Reviewer #3 (Public review):

Summary:

The manuscript uses live imaging to study the role of microtubules in the movement of ribeye aggregates in neuromast hair cells in zebrafish. The main findings are that

(1) Ribeye aggregates, assumed to be ribbon precursors, move in a directed motion toward the active zone;<br /> (2) Disruption of microtubules and kif1aa increases the number of ribeye aggregates and decreases the number of mature synapses.

The evidence for point 2 is compelling, while the evidence for point 1 is less convincing. In particular, the directed motion conclusion is dependent upon fitting of mean squared displacement that can be prone to error and variance to do stochasticity, which is not accounted for in the analysis. Only a small subset of the aggregates meet this criteria and one wonders whether the focus on this subset misses the bigger picture of what is happening with the majority of spots.

Strengths:

(1) The effects of Kif1aa removal and nocodozole on ribbon precursor number and size is convincing and novel.<br /> (2) The live imaging of Ribeye aggregate dynamics provides interesting insight into ribbon formation. The movies showing fusion of ribeye spots are convincing and the demonstrated effects of nocodozole and kif1aa removal on the frequency of these events is novel.<br /> (3) The effect of nocodozole and kif1aa removal on precursor fusion is novel and interesting.<br /> (4) The quality of the data is extremely high and the results are interesting.

Weaknesses:

(1) To image ribeye aggregates, the investigators overexpressed Ribeye-a TAGRFP under control of a MyoVI promoter. While it is understandable why they chose to do the experiments this way, expression is not under the same transcriptional regulation as the native protein and some caution is warranted in drawing some conclusions. For example, the reduction in the number of puncta with maturity may partially reflect regulation of the MyoVI promoter with hair cell maturity. Similarly, it is unknown whether overexpression has the potential to saturate binding sites (for example to motors), which could influence mobility. In the revised manuscript, the authors provide evidence to suggest that overexpression is not at unreasonably high levels, which is reasonable. However, I think it remains important to think of these caveats while reading the paper--especially keeping in mind that expression timing is undoubtedly influenced by the transcriptional control of the exogenous promoter .<br /> (2) The examples of punctae colocalizing with microtubules look clear (fig 1 F-G), but the presentation is anecdotal. It would be better and more informative, if quantified.<br /> (3) It appears that any directed transport may be rare. Simply having an alpha >1 is not sufficient to declare movement to be directed (motor driven transport typically has an alpha approaching 2). Due to randomness of a random walk and errors in fits in imperfect data will yield some spread in movement driven by Brownian motion. Many of the tracks in figure 3H look as thought they might be reasonably fit by a straight line (i.e. alpha = 1).<br /> (4) The "directed motion" shown here does not really resemble motor driven transport observed in other systems (axonal transport, for example) even in the subset that have been picked out as examples here. While the role for microtubules and kif1aa in synapse maturation is strong, it seems likely that this role may be something non-canonical (which would be interesting). In the revision, the authors do an excellent job of considering the issues brought up in point 3 and 4. While perhaps no longer a weakness, I am leaving the critiques here for context for the readers to consider. The added taxol results may not completely settle the issue, but are interesting and provide important information.

Reviewer #1 (Public review):

Summary:

In this manuscript, Chua, Daugherty, and Smith analyze a new set of archaeal 20S proteasomes obtained by cryo-EM that illustrate how the occupancy of the HbYX binding pocket induces gate opening. They do so primarily through a V24Y mutation in the α-subunit. These results are supported by a limited set of mutations in K66 in the α subunit, bringing new emphasis to this unit.

Strengths:

The new structure's analysis is comprehensive, occupying the entire manuscript. As such, the scope of this manuscript is very narrow, but the strength of the data is solid, and they offer an interesting and important new piece to the gate-opening literature.

Weaknesses:

Major Concerns

(1) This manuscript rests on one new cryo-EM structure, leading to a single (albeit convincing) experiment demonstrating the importance of occupying the pocket and moving K66. Could a corresponding bulky mutation at K66 not activate the 20S proteasome?

(2) To emphasize the importance of this work, the authors highlight the importance of gate-opening to human 20S proteasomes. However, the key distinctions between these proteasomes are not given sufficient weight.<br /> (a) As the authors note, the six distinct Rpt C-termini can occupy seven different pickets. However, how these differences would impact activation is not thoroughly discussed.<br /> (b) With those other sites, the relative importance of various pockets, such as the one controlling the α3 N-terminus, should be discussed more thoroughly as a potential critical difference.<br /> (c) These differences can lead to eukaryote 20S gates shifting between closed and open and having a partially opened state. This becomes relevant if the goal is to lead to an activated 20S. It would have been interesting to have archaea 20S with a mix of WT and V24Y α-subunits. However, one might imagine the subclassification problem would be challenging and require an extraordinary number of particles.<br /> (d) Furthermore, the conservation of the amino acids around the binding pocket was not addressed. This seems particularly important in the relative contribution of a residue analogous to K66 or V24.

Reviewer #2 (Public review):

Summary:

The manuscript by Chuah et al. reports the experimental results that suggest the occupancy of the HbYX pockets suffices for proteasome gate opening. The authors conducted cryo-EM reconstructions of two mutant archaeal proteasomes. The work is technically sound and may be of special interest in the field of structural biology of the proteasomes.

Strengths:

Overall, the work incrementally deepens our understanding of the proteasome activation and expands the structural foundation for therapeutic intervention of proteasome function. The evidence presented appears to be well aligned with the existing literature, which adds confidence in the presentation.

Weaknesses:

The paper may benefit from some minor revision by making improvements on the figures and necessary quantitative comparative studies.

Reviewer #1 (Public review):

Summary:

In this manuscript, the authors discovered MYL3 of marine medaka (Oryzias melastigma) as a novel NNV entry receptor, elucidating its facilitation of RGNNV entry into host cells through macropinocytosis, mediated by the IGF1R-Rac1/Cdc42 pathway.

Strengths:

In this manuscript, the authors have performed in vitro and in vivo experiments to prove that MnMYL3 may serve as a receptor for NNV via macropinocytosis pathway. These experiments with different methods include Co-IP, RNAi, pulldown, SPR, flow cytometry, immunofluorescence assays and so on. In general, the results are clearly presented in the manuscript.

Comments on revisions:

The authors have addressed all my comments.

Reviewer #2 (Public review):

Summary:

The manuscript offers an important contribution to the field of virology, especially concerning NNV entry mechanisms. The major strength of the study lies in the identification of MmMYL3 as a functional receptor for RGNNV and its role in macropinocytosis, mediated by the IGF1R-Rac1/Cdc42 signaling axis. This represents a significant advance in understanding NNV entry mechanisms beyond previously known receptors such as HSP90ab1 and HSC70. The data, supported by comprehensive in vitro and in vivo experiments, strongly justify the authors' claims about MYL3's role in NNV infection in marine medaka.

Strengths:

(1) The identification of MmMYL3 as a functional receptor for RGNNV is a significant contribution to the field. The study fills a crucial gap in understanding the molecular mechanisms governing NNV entry into host cells.

(2) The work highlights the involvement of IGF1R in macropinocytosis-mediated NNV entry and downstream Rac1/Cdc42 activation, thus providing a thorough mechanistic understanding of NNV internalization process. This could pave the way for further exploration of antiviral targets.

Comments on revisions:

The authors have addressed the concerns from reviewers. This manuscript can be published in the current form.

Reviewer #3 (Public review):

Summary:

The manuscript presents a detailed study on the role of MmMYL3 in the viral entry of NNV, focusing on its function as a receptor that mediates viral internalization through the macropinocytosis pathway. The use of both in vitro assays (e.g., Co-IP, SPR, and GST pull-down) and in vivo experiments (such as infection assays in marine medaka) adds robustness to the evidence for MmMYL3 as a novel receptor for RGNNV. The findings have important implications for understanding NNV infection mechanisms, which could pave the way for new antiviral strategies in aquaculture.

Strengths:

The authors show that MmMYL3 directly binds the viral capsid protein, facilitates NNV entry via the IGF1R-Rac1/Cdc42 pathway, and can render otherwise resistant cells susceptible to infection. This multifaceted approach effectively demonstrates the central role of MmMYL3 in NNV entry.

Reviewer #2 (Public review):

Summary:

Tanaka et al. investigated the role of CCR4 in early atherosclerosis, focusing on the immune modulation elicited by this chemokine receptor under hypercholesterolemia. The study found that Ccr4 deficiency led to qualitative changes in atherosclerotic plaques, characterized by an increased inflammatory phenotype. The authors further analyzed the CD4 T cell immune response in para-aortic lymph nodes and atherosclerotic aorta, showing an increase mainly in Th1 cells and the Th1/Treg ratio in Ccr4-/-Apoe-/- mice compared to Apoe-/- mice. They then focused on Tregs, demonstrating that Ccr4 deficiency impaired their immunosuppressive function in in vitro assays. Authors also states that Ccr4-deficient Tregs had, as expected, impaired migration to the atherosclerotic aorta. Adoptive cell transfer of Ccr4-/- Tregs to Apoe-/- mice mimicked early atherosclerosis development in Ccr4-/-Apoe-/- mice. Therefore, this work shows that CCR4 plays an important role in early atherosclerosis but not in advanced stages.

Strengths:

Several in vivo and in vitro approaches were used to address the role of CCR4 in early atherosclerosis. Particularly, through the adoptive cell transfer of CCR4+ or CCR4- Tregs, the authors aimed to demonstrate the role of CCR4 in Tregs' protection against early atherosclerosis.

Weaknesses:

Flow cytometry experiments are not well controlled. Dead cells and doublets were not excluded from analysis.

Clinical relevance is unclear.

Comments on revisions:

I thank the authors for addressing my suggestions.<br /> I understand that excluding dead cells would require repeating the entire experiment. However, the authors can at least exclude doublets from the existing flow cytometry data.<br /> I also agree with the more cautious claim regarding the role of CCR4 in Treg migration.

Reviewer #3 (Public review):

Summary

Tanaka and colleagues addressed the role of the C-C chemokine receptor 4 (CCR4) in early atherosclerotic plaque development using ApoE-deficient mice on a standard chow diet as a model. Because several CD4+ T cell subsets express CCR4, they examined whether CCR4-deficiency alters the immune response mediated by CD4+ T cells. By histological analysis of aortic lesions, they demonstrated that the absence of CCR4 promoted the development of early atherosclerosis, with heightened inflammation linked to increased macrophages and pro-inflammatory CD4+ T cells, along with reduced collagen content. Flow cytometry and mRNA expression analysis for identifying CD4+ T cell subsets showed that CCR4 deficiency promoted higher proliferation of pro-inflammatory effector CD4+ T cells in peripheral lymphoid tissues and accumulation of Th1 cells in the atherosclerotic lesions. Interestingly, the increased pro-inflammatory CD4+ T cell response occurred despite the expansion of T CD4+ Foxp3+ regulatory cells (Tregs), found in higher numbers in lymphoid tissues of CCR4-deficient mice, suggesting that CCR4 deficiency interfered with Treg's regulatory actions. The findings contrast with earlier studies in a murine model of advanced atherosclerosis, where CCR4 deficiency did not alter the development of the aortic lesions. The authors included a thoughtful discussion about hypothetical mechanisms explaining these contrasting results, including putative differences in the role played by the CCL17/CCL22-CCR4 axis along the stages of atherosclerosis development in this murine model.

Major strengths

• Demonstration of CCR4 deficiency's impact on early atherosclerosis. CCR4 deficiency effects on the early atherosclerosis development in the Apoe-/-mice model were demonstrated by a quantitative analysis of the lesion area, inflammatory cell content and the expression profile of several pro- and anti-inflammatory markers.<br /> • Analysis of the T CD4+ response in various lymphoid tissues (peripheral and para-aortic lymph nodes and spleen) and the atherosclerotic aorta during the early phase of atherosclerosis in the Apoe-/-mice model. This analysis, combining flow cytometry and mRNA expression, showed that CCR4 deficiency enhanced T CD4+ cell activation, favouring the amplification of the typical biased Th1-mediated inflammatory response observed in the lymphoid tissues of hypercholesterolemic mice.<br /> • Treg transference experiments. Transference of Treg from Apoe-/- or Ccr4-/- Apoe-/- mice to Apoe-/- mice under a standard chow diet was useful for addressing the relevance of CCR4 expression on Tregs for the atheroprotective effect of this regulatory T cell subset during early atherosclerosis.

Major weaknesses

• Methodological Limitations: The controls used in the flow cytometry analysis were suboptimal, as neither cell viability nor doublets were assessed. This may have introduced artifacts, particularly when measuring less-represented cell populations within complex samples, such as in assays evaluating Treg migration to the aorta in atherosclerotic mice.<br /> • Incomplete understanding of CCR4-Mediated Mechanisms: The mechanisms by which CCR4 regulates early inflammation and the development of atherosclerosis were not fully clarified.

I have previously addressed the study limitations and their global impact in my earlier reviews.

Reviewer #1 (Public review):

The authors investigate the function and neural circuitry of reentrant signals in visual cortex. Recurrent signaling is thought to be necessary to common types of perceptual experience that are defined by long-range relationships or prior expectation. Contour illusions - where perceptual objects are implied by stimuli characteristics - are a good example of this. The perception of these illusions is thought to emerge as recurrent signals from higher cortical areas feedback onto early visual cortex, to tell early visual cortex that it should be seeing object contours where none are actually present.

The authors test the involvement of reentrant cortical activity in this kind of perception using a drug challenge. Reentrance in visual cortex is thought to rely on NMDAR-mediated glutamate signalling. The authors accordingly employ an NMDA antagonist to stop this mechanism, looking for the effect of this manipulation on visually evoked activity recorded in EEG.

The motivating hypothesis for the paper is that NMDA antagonism should stop recurrent activity, and that this should degrade perceptual activity supporting perception of a contour illusion, but not other types of visual experience. Results in fact show the opposite. Rather than degrading cortical activity evoked by the illusion, memantine makes it more likely that machine learning classification of EEG will correctly infer the presence of the illusion.

On the face of it, this is confusing. But the paper does a good job of providing possible accounts based on specific details of neurochemical signalling and receptor populations.

I broadly find the paper interesting, graceful, and creative. The hypotheses are clear and compelling, the techniques for both manipulation of brain state and observation of that impact are cutting edge and well suited, and the paper draws clear and convincing conclusions that are made necessary by the results. The work sits at the very interesting crux of systems neuroscience, neuroimaging, and pharmacology.

Reviewer #2 (Public review):

This study presents an important finding to the field interested in recurrent processing and the role of NMDA-receptors herein. The evidence for improved decoding under memantine is convincing, while some open questions remain to be followed up in future studies (the lack of a behavioural effect, why is decoding improved rather than decreased?). It is an excellent example of how an unexpected finding can generate novel research ideas to the mechanisms underlying recurrent processing, suggesting that the answer lies in the differences in the effects of ketamine and memantine, rather than their commonalities.

I would like to thank the authors for the great care they have taken in addressing my concerns. I think the revised manuscript is significantly easier to follow now that specific hypothesis have been formulated in the introduction, and the direction of the results is explicitly stated throughout the manuscript. I further appreciate the dampening of some of the claims that are not completely supported by the appropriate interactions.

I think the resulting manuscript is an incredibly exciting contribution to our understanding of NMDA-receptor function, and a great example of how an unexpected finding can raise questions that could potentially drive the field forward. It shows how NMDA's role in recurrent processing is much more complicate than previously assumed, and reveals that it is not the commonalities between memantine and ketamine that are important in understanding recurrent processing, but rather the differences. I look forward to future studies that will target these differences.

Overall great job.

Reviewer #3 (Public review):

Summary:

In this study, Stein and colleagues use a clever masking/attentional blink paradigm using Kanisza stimuli, coupled with EEG decoding and the NMDA antagonist memantine, to isolate putative neural markers of feedforward, lateral, and feedback processing.

In two elegant experiments, they show that memantine selective influences EEG decoding of only illusory Kanisza surfaces (but not contour continuation or raw contrast), only when unmasked, only when attention is available (not when "blinked"), and only when task-relevant.

This neatly implicates NMDA receptors in the feedback mechanisms that are believed to be involved in inferring illusory Kanisza surfaces, and builds a difficult bridge between the large body of human perceptual experiments and pharmacological and neurophysiological work in animals.

Strengths:

Three key strengths of the paper are 1) its elegant and thorough experimental design, which includes internal replication of some key findings, and 2) the clear pattern of results across the full set of experiments, and 3) its clear writing and presentation of results.

The paper effectively reports a 4-way interaction, with memantine only influencing decoding of surfaces (1) that are unmasked (2), with attention available (3) and task-relevant (4). Nevertheless, the results are very clear, with a clear separation between null effects on other conditions and quite a strong (and thus highly selective) effect on this one intersection of conditions. This makes the pattern of findings very convincing.

Weaknesses:

Overall this is an impressive and important paper. However, to my mind there are two minor weaknesses.

First, despite its clear pattern of neural effects, there is no corresponding perceptual effect. Although the manipulation fits neatly within the conceptual framework, and there are many reasons for not finding such an effect (floor and ceiling effects, narrow perceptual tasks etc), this does leave open the possibility that the observation is entirely epiphenomenal, and that the mechanisms being recorded here are not actually causally involved in perception per se.

Second, although it is clear that there is an effect on decoding in this particular condition, what that means is not entirely clear - particularly since performance improves, rather than decreases. It should be noted here that improvements in decoding performance do not necessarily need to map onto functional improvements, and we should all be careful to remain agnostic about what is driving classifier performance. Here too, the effect of memantine on decoding might be epiphenomenal - unrelated to the information carried in the neural population, but somehow changing the balance of how that is electrically aggregated on the surface of the skull. *Something* is changing, but that might be a neurochemical or electrical side-effect unrelated to actual processing (particularly since no corresponding behavioural impact is observed.)

Comments on revisions:

I think the authors responsed fairly to my comments. Even if they weren't really able to add new insight into why behaviour didn't show the same effects as decoding, they discuss this in the revised text.

Reviewer #1 (Public Review):

Many studies reported findings implying that rhizobial infection is associated with cell cycle re-entry and progression, however, our understanding has been fragmented. This study provides exciting new insights as it represents a comprehensive description of the cell cycle progression during early stages of nodulation using fluorescence markers.

To briefly summarize, the authors first monitor H3.1 / H3.3 replacement to distinguish between replicating (S phase) and non-replicating cells to show that M. truncatula cortex cells along the bacterial infection thread are non-replicating (while neighbors enter the S phase). Nuclear size measurements revealed that these non-replicative cells are in the post-replicative stage (G2) rather than in the pre-replicative G1 phase, which the authors confirm with the Plant Cell Cycle Indicator (PlaCCI) fluorescent marker to track cell cycle progression in more detail. Cortex cells in the trajectory of the infection thread did not accumulate the late G2 marker of the PlaCCI nor the G2/M marker KNOLLE, indicating that these cells indeed remain in G2. Because nuclear size measurements indicated that infected cells are polyploid, the authors used the centromere histone marker CENH3 to determine chromosome number. They find that cortex cells giving rise to the nodule primordium are endomitotic and tetraploid, probably because their cell cycle is halted at centromere separation. Although not a focus of this manuscript, the authors also use their fluorescent tools to track cell cycle progression during arbuscular mycorrhiza symbiosis. They confirm that infected cells transition from a replicating to a non-replicating state (H3.1 to H3.3) with progressing development of the arbuscules. In addition, the CENH3 marker confirms previous findings that cortex cells infected by fungi are endocycling (i.e., DNA synthesis without segregation of replicated parts). This represents an important confirmation of previous findings and contrasts with the situation during nodulation symbiosis, where chromosomes separate after replication.

In general, all microscopy images are of very high quality and support the authors' conclusions. While individually each set of fluorescent markers has its limitations, combined they constitute a powerful tool to track various stages of cell cycle progression in individual root cells during symbiosis. Overall, this is a very strong manuscript that comprehensively elucidates root cell cycle changes during microbial infection.

Reviewer #2 (Public Review):

Cell cycle control during nitrogen-fixing symbiosis is an important topic, but our understanding of the process is poor and lacks resolution, as the nodule is a complex organ with many cell types that undergo profound changes. The authors aim to define the cell cycle state of individual plant cells in the emerging nodule primordium, as a transcellular infection thread passes through the meristem to reach cells deep in the incipient nodule and releases bacteria to form symbiosomes. The authors used a number of cell cycle reporters, such as different Histone 3 variants and cyclins, to follow cell cycle progress in exquisite detail. They showed that the host cells in the path of an infection thread exhibit a cell fate distinct from their immediate neighbors: after entering the S phase similar to their neighbors, these cells exit the cell cycle and enter a special differentiated state. This is likely an important shift that allows the proper passage of the infection thread. Although definitive proof needs more investigation, they showed that a pioneering transcription factor, NF-YA1, likely represses these endoreduplicated cells from completing the cell cycle.

Reviewer #1 (Public review):

The structure of a heterohexameric 3:3 LGI1-ADAM22 complex is resolved by Yamaguchi et al. It reveals the intermolecular LGI1 interactions and its role in bringing three ADAM22 molecules together. This may be relevant for the clustering of axonal Kv1 channels and control over their density. While it is currently not clear if the heterohexameric 3:3 LGI1-ADAM22 complex has a physiological role, the detailed structural information presented here allows to pinpoint mutations or other strategies to probe the relevance of the 3:3 complex in future work.

The experimental work is done to a high standard, and all my comments have been addressed. This new version of the manuscript has been improved substantially, and the figures have been enhanced and clarified.

Reviewer #1 (Public review):

This study exploits novel agent (IMT) that inhibits mitochondrial activity in combination with venetoclax. While the concept is not novel, the agent is novel (inhibitor of the mitochondrial RNA polymerase, described in Nature in other tumor models), and quest for safe mitochondrial inhibitors is highly warranted. The strength is in vivo activity data shown in CLDX and in one of the two AML PDX models tested, and apparent safety of the combination. However, the impact on survival is impressive in CLDX but not in PDX, and unclear why Ven-sensitive PDX is resistant to combination (opposite what cell line data show). There is no real evidence that this agent overcome Ven resistance, which could be done for example in primary AML cells. Finally, no on-target pharmacodynamic endpoints are measured in vivo to support the activity of the compound on mitochondrial activity at the doses used (which are safe).

Both Reviewers requested to demonstrate that IMT1 inhibits the target at doses used in vitro or in vivo; while the prior paper showed this for original compound, it is imperative to demonstrate this for this modified agent in a different tumor type such as AML.

These points have not been addressed in the Revision.

Reviewer #2 (Public review):

Summary:

The manuscript by Arabanian and colleagues presents studies showing how inhibition of mitochondrial transcription and replication with a novel inhibitor of the mitochondrial polymerase, IMT, can promote AML cell death in combination with the Bcl2 inhibitor venetoclax. They further show that this combinatorial efficacy is evident in vivo in both the AML cell line MV411 and in a PDX model. Given the multiple studies showing the importance of Oxphos in maintaining AML cell survival, the current studies provide an additional strategy to inhibit Oxphos and thus improve the therapeutic management of AML.

Strengths:

A novel aspect of this work is that IMT is a new class of mitochondrial inhibitor that acts through inhibiting the mitochondrial polymerase. In addition, the demonstration of therapeutic efficacy both in vitro and in vivo (including with PDX), together with some data showing minimal toxicity, adds to the impact of this work. Their overall conclusion that IMT increases the potency of Vex in treating AMLs is supported.

Comments on revisions:

In all, the authors responded to most of the critiques, while two of the major critiques were not experimentally addressed. The work will still have potential impact, but will depend on further studies under more clinically relevant conditions and with a better understanding of drug effects.

Reviewer #1 (Public review):

Summary:

The study by Deng et al reports single cell expression analysis of developing mouse hearts and examines the requirements for cardiac fibroblasts in heart maturation. The work includes extensive gene expression profiling and bioinformatic analysis. The prenatal fibroblast ablation studies show new information on the requirement of these cells on heart maturation before birth.

The strengths of the manuscript are the new single cell datasets and comprehensive approach to ablating cardiac fibroblasts in pre and postnatal development in mice. Extensive data are presented on mouse embryo fibroblast diversity and morphology in response to fibroblast ablation. Histological data support localization of major cardiac cell types and effects of fibroblast ablation on cardiac gene expression at different times of development.

A weakness of the study is that the major conclusions regarding collagen signaling and heart maturation are based on gene expression patterns and are not functionally validated.

Reviewer #2 (Public review):

This study aims to elucidate the role of fibroblasts in regulating myocardium and vascular development through signaling to cardiomyocytes and endothelial cells. This focus is significant, given that fibroblasts, cardiomyocytes, and vascular endothelial cells are the three primary cell types in the heart. The authors employed a Pdgfra-CreER-controlled diphtheria toxin A (DTA) system to ablate fibroblasts at various embryonic and postnatal stages, characterizing the resulting cardiac defects, particularly in myocardium and vasculature development. Single-cell RNA sequencing (scRNA-seq) analysis of the ablated hearts identified collagen as a crucial signaling molecule from fibroblasts that influences the development of cardiomyocytes and vascular endothelial cells.

This is an interesting manuscript; however, there are several major issues, including an over-reliance on the scRNA-seq data, which shows inconsistencies between replicates.

Some of the major issues are described below.

(1) The CD31 immunostaining data (Figure 3B-G) indicate a reduction in endothelial cell numbers following fibroblast deletion using PdgfraCreER+/-; RosaDTA+/- mice. However, the scRNA-seq data show no percentage change in the endothelial cell population (Figure 4D). Furthermore, while the percentage of Vas_ECs decreased in ablated samples at E16.5, the results at E18.5 were inconsistent, showing an increase in one replicate and a decrease in another, raising concerns about the reliability of the RNA-seq findings.

(2) Similarly, while the percentage of Ven_CMs increased at E18.5, it exhibited differing trends at E16.5 (Fig. 4E), further highlighting the inconsistency of the scRNA-seq analysis with the other data.

(3) Furthermore, the authors noted that the ablated samples had slightly higher percentages of cardiomyocytes in the G1 phase compared to controls (Fig. 4H, S11D), which aligns with the enrichment of pathways related to heart development, sarcomere organization, heart tube morphogenesis, and cell proliferation. However, it is unclear how this correlates with heart development, given that the hearts of ablated mice are significantly smaller than those of controls (Figure 3E). Additionally, the heart sections from ablated samples used for CD31/DAPI staining in Figure 3F appear much larger than those of the controls, raising further inconsistencies in the manuscript.

(4) The manuscript relies heavily on the scRNA-seq dataset, which shows inconsistencies between the two replicates. Furthermore, the morphological and histological analyses do not align with the scRNA-seq findings.

(5) There is a lack of mechanistic insight into how collagen, as a key signaling molecule from fibroblasts, affects the development of cardiomyocytes and vascular endothelial cells.

(6) In Figure 1B, Col1a1 expression is observed in the epicardial cells (Figure 1A, E11.5), but this is not represented in the accompanying cartoon.

(7) Do the PdgfraCreER+/-; RosaDTA+/- mice survive after birth when induced at E15.5, and do they exhibit any cardiac defects?

Reviewer #3 (Public review):

Summary:

The authors investigated fibroblasts' communication with key cell types in developing and neonatal hearts, with focus on critical roles of fibroblast-cardiomyocyte and fibroblast-endothelial cells network in cardiac morphogenesis. They tried to map the spatial distribution of these cell types and reported the major pathways and signaling molecules driving the communication. They also used Cre-DTA system to ablate Pdgfra labeled cells and observed myocardial and endothelial cell defects at development. They screened the pathways and genes using sequencing data of ablated heart. Lastly they reported a compensatory collagen expression in long term ablated neonate heart. Overall, this study provides us with important insight on fibroblasts' roles in cardiac development and will be a powerful resource for collagens and ECM focused research.

Strengths:

The authors utilized good analyzing tools to investigate on multiple database of single cell sequencing and Multi-seq. They identified significant pathways, cellular and molecular interactions of fibroblasts. Additionally, they compared some of their analytic findings with human database, and identified several groups of ECM genes with varying roles in mice.

Weaknesses:

This study is majorly based on sequencing data analysis. At the bench, they used very strident technique to study fibroblast functions by ablating one of the major cell population of heart. Also, experimental validation of their analyzed downstream pathways will be required eventually.

Reviewer #1 (Public review):

The authors, Zhang et al., demonstrate the beneficial effects of treating degenerate human primary intervertebral disc (IVD) cells with recombinant human PDGF-AB/BB on the senescence transcriptomic signatures. Utilizing a combination of degenerate cells from elderly humans and experimentally induced senescence in young, healthy IVD cells, the authors show the therapeutic effects on mRNA transcription as well as cellular processes through informatics approaches.

One notable strength of this study is the use of human primary cells and recombinant forms of human PDGF-AB/BB proteins, which increases the translational potential of these in vitro studies. The manuscript is well-written, and the informatics analyses are thorough and clearly presented.

Comments on revisions:

The revised manuscript adds greater clarity, and the impact of the study is greatly enhanced.

Reviewer #2 (Public review):

Summary:

This work highlights a novel role for platelet-derived growth factor (PDGF) in mitigating cellular senescence associated with age-related and painful intervertebral disc degeneration. Prior literature has demonstrated the importance of accumulation of senescent cells in mediating many of the pathological effects associated with the degenerate disc joint, such as inflammation and tissue breakdown. In this study the authors treat clinically relevant human nucleus pulposus and annulus fibrosus cells from patients undergoing discectomy with recombinant PDGF-AB/BB for 5 days and then deep phenotyped the outcomes using bulk RNA sequencing. In addition they irradiated healthy human disc cells which they subsequently treated with PDGF-AB/BB examining the expression of SASP-related markers and also PDGFRA receptor gene expression. Overall PDGF was able to down-regulate many senescent associated pathways and the degenerate phenotype in IVD cells. Altered pathways were associated with neurogenesis, mechanical stimuli, metabolism, cell cycle, reactive oxygen species and mitochondrial dysfunction. Overall the authors achieved their aims and the results by and large support their conclusions although improvements could be made to enhance the rigor of the study and findings

Strengths:

A major strength of this study is the use of human cells from patients undergoing discectomy for disc herniation as well as access to healthy human cells. Investigating the role of PDGF regarding cellular senescence in the degenerate disc joint is novel and an underexplored area of research which is a significant contribution to the field of spine. This study highlights a potential target for addressing cellular senescence where most of the prior focus has been on senolytic drugs. Such studies have broad implications to other age-related diseases where senescence plays a major role. The use of transcriptomics and therefore an unbiased approach to investigating the role of PDGF is also considered a strength as is the follow-up studies involving irradiating healthy human disc cells and treating these cells with PDGF. The combined assessment of both nucleus pulposus and annulus fibrosus cells in the context of these studies adds to the impact.

Weaknesses:

A weakness of these studies relates to qualitative data presented for the B-galactosidase assay. Quantification of such data sets would greatly strengthen the studies and lend further support to the hypotheses. The study in its current form could be strengthened by the inclusion of mechanistic studies probing the downstream PDGF receptor associated pathways for example specifically targeting or modulating the activity of the PDGF receptor PDGFRA.

Reviewer #1 (Public review):

Summary:

Zhao and colleagues employ Drosophila nephrocytes as a model to investigate the effects of a high-fat diet on these podocyte-like cells. Through a highly focused analysis, they initially confirm previous research in their hands demonstrating impaired nephrocyte function and move on to observe the mislocalization of a slit diaphragm-associated protein (pyd) and a knock-in into the locus of the Drosophila nephrin (sns). Employing another reporter construct, they identify activation of the JAK/STAT signaling pathway in nephrocytes. Subsequently, the authors demonstrate the involvement of this pathway in nephrocyte function from multiple angles, using a gain-of-function construct, silencing of an inhibitor, and ectopic overexpression of a ligand. Silencing the effector Stat92E via RNAi or inhibiting JAK/STAT with Methotrexate effectively restored impaired nephrocyte function and slit diaphragm architecture induced by a high-fat diet, while showing no impact under normal dietary conditions.

Strengths:

The findings establish a link between JAK/STAT activity and the impact of a high-fat diet on nephrocytes. This nicely underscores the importance of organ crosstalk for nephrocytes and supports a potential role for JAK/STAT in diabetic nephropathy, as previously suggested by other models.

Weaknesses:

While the analysis provides valuable insights, it appears somewhat over-reliant on tracer uptake in certain instances. Clinical inferences based on a Drosophila model should be interpreted with caution.

Reviewer #2 (Public review):

Summary:

In their manuscript, Zhao et al. describe a link between JAK-STAT pathway activation in nephrocytes upon a high-fat diet. Nephrocytes are the homologs to mammalian podocytes, and it has been previously shown that metabolic syndrome and obesity is associated with worse outcomes for chronic kidney disease. A study from 2021 (Lubojemska et al.) could already confirm a severe nephrocyte phenotype upon feeding Drosophila a high fat diet and also linking lipid overflow by expressing adipose triglyceride lipase in the fat body to nephrocyte dysfunction. In this study, the authors identified a second pathway and mechanism, how lipid dysregulation impact on nephrocyte function. In detail, they show an activation of JAK-STAT signaling in nephrocytes upon feeding a high-fat diet, which was induced by Upd2 expression (a leptin-like hormone) in the fat body, the adipose tissue in Drosophila. Further, they could show genetic and pharmacological interventions can reduce JAK-STAT activation and thereby prevent the nephrocyte phenotype in the high-fat diet model.

Strengths:

The strength of this study is the combination of genetic tools and pharmacological intervention to confirm a mechanistic link between the fat body/adipose tissue and nephrocytes. Inter-organ communication is crucial in the development of several diseases, but the underlying mechanisms are only poorly understood. Using Drosophila, it is possible to investigate several players of one pathway, here JAK-STAT. This was done, by investigating the functional role of Hop, Socs36E and Stat92E in nephrocytes and has also been combined with feeding a high-fat diet, to assess restoration of nephrocyte morphology and function by inhibiting JAK-STAT signaling. Adding a translational approach was done by inhibiting JAK-STAT signaling with methotrexate, which also resulted in attenuated nephrocyte dysfunction. Expression of the leptin-like hormone upd2 in the fat body is a good approach to study inter-organ communication and the impact of other organs/tissue on nephrocyte function and expands their findings from nephrocyte function towards whole animal physiology.

Weaknesses:

Although the general findings of this study are of great interest, the number of flies investigated for the majority of the experiments is very low (6 flies). Also it is not clear whether the 6 flies used are from independent experiments to exclude differences in food/diet.

Reviewer #1 (Public review):

Summary:

The authors use microscopy experiments to track the gliding motion of filaments of the cyanobacteria Fluctiforma draycotensis. They find that filament motion consists of back and forth trajectories along a "track", interspersed with reversals of movement direction, with no clear dependence between filament speed and length. It is also observed that longer filaments can buckle and form plectonemes. A computational model is used to rationalize these findings.

Strengths:

Much work in this field focuses on molecular mechanisms of motility; by tracking filament dynamics this work helps to connect molecular mechanisms to environmentally and industrially relevant ecological behavior such as aggregate formation.

The observation that filaments move on tracks is interesting and potentially ecologically signifiant.

The observation of rotating membrane-bound protein complexes and tubular arrangement of slime around the filament provide important clues to the mechanism of motion.

The observation that long filaments buckle has potential to shed light on the nature of mechanical forces in the filaments, e.g. through study of the length dependence of buckling.

The comparison between motility on agar and on glass is interesting since it shows that filaments have both intrinsic propensity to reverse (that is seen on glass) and mechanically triggered reversal (that is seen on agar when the filament reaches the end of a track).

Weaknesses:

The manuscript makes the interesting statement that the distribution of speed vs filament length is uniform, which would constrain the possibilities for mechanical coupling between the filaments. However Fig 2C does not show a uniform distribution but rather an apparent lack of correlation between speed and filament length, although the statistical degree of correlation is not given. In my view, Fig 2C should not be described as a uniform distribution since mathematically that means something very different than what is shown here. Instead the figure should be described quantitatively with the use of a measured correlation coefficient. This also applies to Fig. S3A.

The statement "since filament speed results from a balance between propulsive forces and drag, these observations of no or positive correlation between filament speed and length show that all (or a fixed proportion of) cells in a filament contribute to propulsive force generation" helps to clarify the important link between Fig 2C and the concept that all cells contribute, but I think this statement is not obvious for many readers, and could be made clearer, e.g. by the use of a simple mathematical model for a chain of bacterial that accounts for drag forces and propulsion forces for each bacterium.

The authors have now clarified that the computational model is 1D and cannot explain the coupling between rotation, slime generation and motion. I find it encouraging and important that model predictions for the dwell time distributions (Fig S12 and S13) are similar to experimental measurements, but I think it would be better to put these results in the main text, together also with Fig S4. If these important results are in the supplement it is harder for the reader to assess the match between model and experiments.

Filament buckling is not analysed in quantitative detail, but the authors have now clarified that this will be the topic of future work with a 2D or 3D computational model.

Reviewer #2 (Public review):

Summary:

The authors combined time-lapse microscopy with biophysical modeling to study the mechanisms and timescales of gliding and reversals in filamentous cyanobacterium Fluctiforma draycotensis. They observed the highly coordinated behavior of protein complexes moving in a helical fashion on cells' surfaces and along individual filaments as well as their de-coordination, which induces buckling in long filaments.

Strengths:

The authors provided concrete experimental evidence of cellular coordination and de-coordination of motility between cells along individual filaments. The evidence is comprised of individual trajectories of filaments that glide and reverse on surfaces as well as the helical trajectories of membrane-bound protein complexes that move on individual filaments and are implicated in generating propulsive forces.

Limitations:

The biophysical model is one-dimensional and thus does not capture the buckling observed in long filaments. I expect that the buckling contains useful information since it reflects the competition between bending rigidity, the speed at which cell synchronization occurs, and the strength of the propulsion forces.

Future directions:

The study highlights the need to identify molecular and mechanical signaling pathways of cellular coordination. In analogy to the many works on the mechanisms and functions of multi-ciliary coordination, elucidating coordination in cyanobacteria may reveal a variety of dynamic strategies in different filamentous cyanobacteria.

Reviewer #3 (Public review):

Summary:

The authors present new observations related to the gliding motility of the multicellular filamentous cyanobacteria Fluctiforma draycotensis. The bacteria move forward by rotating their about their long axis, which causes points on the cell surface to move along helical paths. As filaments glide forward they form visible tracks. Filaments preferentially move within the tracks. The authors device a simple model in which each cell in a filament exerts a force that either pushes forwards or backwards. Mechanical interactions between cells cause neighboring cells to align the forces they exert. The model qualitatively reproduces the tendency of filaments to move in a concerted direction and reverse at the end of tracks.

The authors seek to understand how cells in a filament synchronize their motion to move in a concerted direction. This question connects to the evolution of multicellular life and so is important well beyond the specific field of cyanobacterial locomotion.

Strengths:

The biophysical model used to describe cell-cell coordination of locomotion is clear and reasonable. This model provides a useful phenomenological framework in which to consider the roles of individual cells in the coordinated motion of the group. The qualitative consistency between theory and observation suggests that this model captures some essential qualities of the true system.

The observation that filaments reverse at the ends of tracks is compelling, but difficult to clearly connect to any one microscopic model.

The observations of helical motion of the filament are compelling.

Weaknesses:

The comparison of theory and observation is mainly qualitative. While the authors have done a good job fitting the observations to the theory, it is not possible to systematically vary parameters, independently estimate parameter values, or apply external forces. Consequently, more experiments are needed before the proposed model can the accepted or rejected. This manuscript provides a promising hypothesis but not a compelling justification for it.

Reviewer #1 (Public review):

Summary:

The manuscript by Yamamoto et al. presents a model by which the four main axes of the limb are required for limb regeneration to occur in the axolotl. A longstanding question in regeneration biology is how existing positional information is used to regenerate the correct missing elements. The limb provides an accessible experimental system by which to study the involvement of the anteroposterior, dorsoventral, and proximodistal axes in the regenerating limb. Extensive experimentation has been performed in this area using grafting experiments. Yamamoto et al. use the accessory limb model and some molecular tools to address this question. There are some interesting observations in the study. In particular, one strength is the potent induction of accessory limbs in the dorsal axis with BMP2+Fgf2+Fgf8, which is very interesting.

Strengths:

The manuscript presents some novel phenotypes generated in axolotl limbs due to Wnt signaling. This is generally the first example in which Wnt signaling has provided a gain-of-function in the axolotl limb model. They also present a potent way of inducing limb patterning in the dorsal axis by the addition of just beads loaded with Bmp2+Fgf8+Fgf2.

Weaknesses:

Although interesting, the study makes bold claims about determining the molecular basis of DV positional cues, but the experimental evidence is not definitive and does not take into account the previous work on DV patterning in the amniote limb. Also, testing the hypothesis on blastemas after limb amputation would be needed to support the strong claims in the study. There are several examples of very strong claims, but the evidence lacks support for these claims.

Reviewer #2 (Public review):

Summary:

This study explores how signals from all sides of a developing limb, front/back and top/bottom, work together to guide the regrowth of a fully patterned limb in axolotls, a type of salamander known for its impressive ability to regenerate limbs. Using a model called the Accessory Limb Model (ALM), the researchers created early limb regenerates (called blastemas) with cells from different sides of the limb. They discovered that successful limb regrowth only happens when the blastema contains cells from both the top (dorsal) and bottom (ventral) of the limb. They also found that a key gene involved in front/back limb patterning, called Shh (Sonic hedgehog), is only turned on when cells from both the dorsal and ventral sides come into contact. The study identified two important molecules, Wnt10B and FGF2, that help activate Shh when dorsal and ventral cells interact. Finally, the authors propose a new model that explains how cells from all four sides of a limb, dorsal, ventral, anterior (front), and posterior (back), contribute at both the cellular and molecular level to rebuilding a properly structured limb during regeneration.

Strengths:

The techniques used in this study, like delicate surgeries, tissue grafting, and implanting tiny beads soaked with growth factors, are extremely difficult, and only a few research groups in the world can do them successfully. These methods are essential for answering important questions about how animals like axolotls regenerate limbs with the correct structure and orientation. To understand how cells from different sides of the limb communicate during regeneration, the researchers used a technique called in situ hybridization, which lets them see where specific genes are active in the developing limb. They clearly showed that the gene Shh, which helps pattern the front and back of the limb, only turns on when cells from both the top (dorsal) and bottom (ventral) sides are present and interacting. The team also took a broad, unbiased approach to figure out which signaling molecules are unique to dorsal and ventral limb cells. They tested these molecules individually and discovered which could substitute for actual dorsal and ventral cells, providing the same necessary signals for proper limb development. Overall, this study makes a major contribution to our understanding of how complex signals guide limb regeneration, showing how different regions of the limb work together at both the cellular and molecular levels to rebuild a fully patterned structure.

Weaknesses:

Because the expressional analyses are performed on thin sections of regenerating tissue, they provide only a limited view of the gene expression patterns in their experiments, opening the possibility that they could be missing some expression in other regions of the blastema. Additionally, the quantification method of the expressional phenotypes in most of the experiments does not appear to be based on a rigorous methodology. Therefore, performing alternate expressional analysis, using RNA-seq or qRT-PCR (for example) on the entire blastema would help validate that the authors are not missing something.

Overall, the number of replicates per sample group is quite low (sometimes as low as 3), which is especially risky with challenging techniques like the ones the authors employ. The authors don't appear to have performed a power analysis to calculate the number of animals used in each experiment that are sufficient to identify possible statistical differences between groups. Increasing the sample sizes would substantially increase the rigor of their experiments.

Likewise, the authors' use of an AI-generated algorithm to quantify symmetry on the dorsal/ventral axis, and this approach doesn't appear to account for possible biases due to tissue sectioning angles. They also appear to arbitrarily pick locations in each sample group to compare symmetry measurements. There are other methods, which include using specific muscle groups and nerve bundles as dorsal/ventral landmarks, that would more clearly show differences in symmetry.

Reviewer #3 (Public review):

Summary:

After salamander limb amputation, the cross-section of the stump has two major axes: anterior-posterior and dorsal-ventral. Cells from all axial positions (anterior, posterior, dorsal, ventral) are necessary for regeneration, yet the molecular basis for this requirement has remained unknown. To address this gap, Yamamoto et al. took advantage of the ALM assay, in which defined positional identities can be combined on demand and their effects assessed through the outgrowth of an ectopic limb. They propose a compelling model in which dorsal and ventral cells communicate by secreting Wnt10b and Fgf2 ligands, respectively, with this interaction inducing Shh expression in posterior cells. Shh was previously shown to induce limb outgrowth in collaboration with anterior Fgf8 (PMID: 27120163). Thus, this study completes a concept in which four secreted signals from four axial positions interact for limb patterning. Notably, this work firmly places dorsal-ventral interactions upstream of anterior-posterior, which is striking for a field that has been focussed on anterior-posterior communication. The ligands identified (Wnt10b, Fgf2) are different from those implicated in dorsal-ventral patterning in the non-regenerative mouse and chick models. The results in the context of ALM/ectopic limb engineering are impressive, but the authors do not extend their experiments to assay 'normal' regeneration after amputation.

Strengths:

(1) The ALM and use of GFP grafts for lineage tracing (Figures 1-3) take full advantage of the salamander model's unique ability to outgrow patterned limbs under defined conditions. As far as I am aware, the ALM has not been combined with precise grafts that assay 2 axial positions at once, as performed in Figure 3. The number of ALMs performed in this study deserves special mention, considering the challenging surgery involved.

(2) The authors identify that posterior Shh is not expressed unless both dorsal and ventral cells are present. This echoes previous work in mouse limb development models (AER/ectoderm-mesoderm interaction), but this link between axes was not known in salamanders. The authors elegantly reconstitute dorsal-ventral communication by grafting, finding that this is sufficient to trigger Shh expression (Figure 3 - although see also the Weaknesses section.)

(3) Impressively, the authors discovered two molecules sufficient to substitute dorsal or ventral cells through electroporation into dorsal- or ventral-depleted ALMs (Figure 5). These molecules did not change the positional identity of target cells. The same group previously identified the ventral factor (Fgf2) to be a nerve-derived factor essential for regeneration. In Figure 6, the authors demonstrate that nerve-derived factors, including Fgf2, are alone sufficient to grow out ectopic limbs from a dorsal wound. Limb induction with a 3-factor cocktail without supplementing with other cells is conceptually important for regenerative engineering.

(4) The writing style and presentation of results are very clear.

Weaknesses:

(1) The expression data are the weakest part of this study.

• Despite being a central message, I found the Shh in situs unconvincing (e.g. Figure 2I, 3C, 5C), especially without sense probe controls. An additional assay would be essential to make the Shh data convincing - perhaps like in Figure 5D (qPCR?), RNA-sequencing, or a downstream target gene.

• It is not clear what the n numbers mean for the in situ data (slides analysed / number of biological samples / other?). This is crucial to understanding the reliability of the results.

• The authors do not assay where and when Wnt10b and Fgf2 are expressed beyond the bulk RNA-sequencing (which presumably contains both epidermis and mesenchyme cells). This is a shame, as understanding which cell types express these molecules, and when, would be important for understanding the mechanism.

(2) It is important to consider that the ALM is not 'regeneration', even if the authors have previously argued that ALM bumps and regenerating blastemas are equivalent (PMID: 17959163). The start- and end- points of ALM are different from regeneration, even though there are undoubtedly common principles involved. Thus, I find the word 'regeneration' in the title and last sentence of the abstract unsubstantiated unless evidence is provided that the same mechanisms (Wnt10b/Fgf2/Shh) function during normal limb regeneration.

(3) Drawing the exact boundaries of the Ant/Pos/Dor/Ven BL and grafts in the cartoon in Figure 1 (with respect to anatomical landmarks) would help to better understand the experiments in Figures 3 and 4.

(4) I find the 'positional cue' and 'positional value' terminology confusing, despite the authors' efforts. It is not clear if they refer to cell autonomous or secreted signals, and, as the authors mention, the definitions partially overlap. Lmx1b is defined as a positional value, even though it is necessary and sufficient for dorsal identity (so, isn't it positional information?). Much simpler would be to describe Wnt10b and Fgf2 as what they are: dorsally or ventrally expressed signals that substitute for dorsal or ventral tissue without inducing changes in positional information.

Overall appraisal:

This is a logical and well-executed study that creatively uses the axolotl model to advance an important framework for understanding limb patterning. The reliability of the Shh expression data is a weak point in this otherwise impressive study. The relevance of the mechanisms to normal limb regeneration is not substantiated.

Joint Public Review:

Summary:

The authors investigate how stochastic and deterministic factors are integrated in cell fate decisions, using Dictyostelium discoideum as a model system. They show that cells in different cell cycle phases (a deterministic factor) are predisposed to different fates, albeit with deviations, when exposed to the same environmental stimulus. However, gene expression variability (a stochastic factor) enhances the robustness of cellular responses to environmental cues that disrupt the cell cycle.

Using a simple, tractable mathematical model, the authors demonstrate that cell fate decisions in D. discoideum depend on a combination of deterministic and stochastic factors, i.e., cell cycle phase and gene expression variability, respectively. They then identify Set1 - a key regulator of gene expression variability - indicate the mechanism through which it modulates this variability, and link it to a phenotype in D. discoideum development. Finally, they confirm that gene expression variability contributes to the robustness of the cell's response to environmental disruptions that interfere with the cell cycle.

Strengths:

The authors are careful in the choice of their experiments and in measuring gene expression variability, using methods that account for expected trends with average gene expression.

Weaknesses:

However, in terms of mathematical modelling, it would be important to rule out sources of stochasticity (other than gene expression variability), and also to consider cases where stochastic factors are not necessarily completely independent of the deterministic ones.

Reviewer #1 (Public review):

Summary:

These authors have asked how lytic phage predation impacts antibiotic resistance and virulence phenotypes in methicillin-resistant Staphylococcus aureus (MRSA). They report that staphylococcal phages cause MRSA strains to become sensitized to b-lactams and to display reduced virulence. Moreover, they identify mutations in a set of genes required for phage infection that may impact antibiotic resistance and virulence phenotypes.

Strengths:

Phage-mediated re-sensitization to antibiotics has been reported previously but the underlying mutational analyses have not been described. These studies suggest that phages and antibiotics may target similar pathways in bacteria.

Weaknesses:

One limitation is the lack of mechanistic investigations linking particular mutations to the phenotypes reported here. This limits the impact of the work.

Another limitation of this work is the use of lab strains and a single pair of phages. However, while incorporation of clinical isolates would increase the translational relevance of this work it is unlikely to change the conclusions.

Comments on revisions:

The authors have addressed my concerns.

Reviewer #2 (Public review):

Summary:

The work presented in the manuscript by Tran et al deals with bacterial evolution in the presence of bacteriophage. Here, authors have taken three methicillin-resistant S. aureus strains that are also resistant to beta-lactams. Eventually, upon being exposed to phage, these strains develop beta-lactam sensitivity. Besides this, the strains also show other changes in their phenotype such as reduced binding to fibrinogen and hemolysis.

Strengths:

The experiments carried out are convincing to suggest such in vitro development of sensitivity to the antibiotics. Authors were also able to "evolve" phage in similar fashion thus showing enhanced virulence against the bacterium. In the end, authors carry out DNA sequencing of both evolved bacteria and phage and show mutations occurring in various genes. Overall, the experiments that have been carried out are convincing.

Weaknesses:

None. In the current version of the manuscript, I find the study complete.

Reviewer #1 (Public review):

Summary:

These authors have asked how lytic phage predation impacts antibiotic resistance and virulence phenotypes in methicillin-resistant Staphylococcus aureus (MRSA). They report that staphylococcal phages cause MRSA strains to become sensitized to b-lactams and to display reduced virulence. Moreover, they identify mutations in a set of genes required for phage infection that may impact antibiotic resistance and virulence phenotypes.

Strengths:

Phage-mediated re-sensitization to antibiotics has been reported previously but the underlying mutational analyses have not been described. These studies suggest that phages and antibiotics may target similar pathways in bacteria.

Weaknesses:

One limitation is the lack of mechanistic investigations linking particular mutations to the phenotypes reported here. This limits the impact of the work.

Another limitation of this work is the use of lab strains and a single pair of phages. However, while incorporation of clinical isolates would increase the translational relevance of this work it is unlikely to change the conclusions.

Reviewer #2 (Public review):

Summary:

The work presented in the manuscript by Tran et al deals with bacterial evolution in the presence of bacteriophage. Here, the authors have taken three methicillin-resistant S. aureus strains that are also resistant to beta-lactams. Eventually, upon being exposed to phage, these strains develop beta-lactam sensitivity. Besides this, the strains also show other changes in their phenotype such as reduced binding to fibrinogen and hemolysis.

Strengths:

The experiments carried out are convincing to suggest such in vitro development of sensitivity to the antibiotics. Authors were also able to "evolve" phage in a similar fashion thus showing enhanced virulence against the bacterium. In the end, authors carry out DNA sequencing of both evolved bacteria and phage and show mutations occurring in various genes. Overall, the experiments that have been carried out are convincing.

Weaknesses:

Although more experiments are not needed, additional experiments could add more information. For example, the phage gene showing the HTH motif could be reintroduced in the bacterial genome and such a strain can then be assayed with wildtype phage infection to see enhanced virulence as suggested. At least one such experiment proves the discoveries regarding the identification of mutations and their outcome. Secondly, I also feel that authors looked for beta-lactam sensitivity and they found it. I am sure that if they look for rifampicin resistance in these strains, they will find that too. In this case, I cannot say that the evolution was directed to beta-lactam sensitivity; this is perhaps just one trait that was observed. This is the only weakness I find in the work. Nevertheless, I find the experiments convincing enough; more experiments only add value to the work.

Reviewer #2 (Public review):

Summary:

This study uses in vivo multimodal high-resolution imaging to track how microglia and neutrophils respond to light-induced retinal injury from soon after injury to 2 months post-injury. The in vivo imaging finding was subsequently verified by ex vivo study. The results suggest that despite the highly active microglia at the injury site, neutrophils were not recruited in response to acute light-induced retinal injury.

Strengths:

An extremely thorough examination of the cellular-level immune activity at the injury site. In vivo imaging observations being verified using ex vivo techniques is a strong plus.

Reviewer #3 (Public review):

Summary

This work investigated the immune response in the murine retina after focal laser lesions. These lesions are made with close to 2 orders of magnitude lower laser power than the more prevalent choroidal neovascularization model of laser ablation. Histology and OCT together show that the laser insult is localized to the photoreceptors and spares the inner retina, the vasculature and the pigment epithelium. As early as 1-day after injury, a loss of cell bodies in the outer nuclear layer is observed. This is accompanied by strong microglial proliferation to the site of injury in the outer retina where microglia do not typically reside. The injury did not seem to result in the extravasation of neutrophils from the capillary network, constituting one of the main findings of the paper. The demonstrated paradigm of studying the immune response and potentially retinal remodeling in the future in vivo is valuable and would appeal to a broad audience in visual neuroscience.

Strengths

Adaptive optics imaging of murine retina is cutting edge and enables non-destructive visualization of fluorescently labeled cells in the milieu of retinal injury. As may be obvious, this in vivo approach is a benefit for studying fast and dynamic immune processes on a local time scale - minutes and hours, and also for the longer days-to-months follow-up of retinal remodeling as demonstrated in the article. In certain cases, the in vivo findings are corroborated with histology.

The analysis is sound and accompanied by stunning video and static imagery. A few different sets of mouse models are used: a) two different mouse lines, each with a fluorescent tag for neutrophils and microglia, b) two different models of inflammation - endotoxin-induced uveitis (EAU) and laser ablation are used to study differences in the immune interaction.

One of the major advances in this article is the development of the laser ablation model for 'mild' retinal damage as an alternative to the more severe neovascularization models. This model would potentially allow for controlling the size, depth and severity of the laser injury opening interesting avenues for future study.

The time-course, 2D and 3D spatial activation pattern of microglial activation are striking and provide an unprecedented view of the retinal response to mild injury.

Editor's note: The authors have addressed all the previous concerns raised by the reviewers.

Reviewer #1 (Public review):

Summary:

In this study, authors explored how galanin affects whole-brain activity in larval zebrafish using wide-field Ca2+ imaging, genetic modifications, and drugs that increase brain activity. The authors conclude that galanin has a sedative effect on the brain under normal conditions and during seizures, mainly through the galanin receptor 1a (galr1a). However, acute "stressors(?)" like pentylenetetrazole (PTZ) reduce galanin's effects, leading to increased brain activity and more seizures. Authors claim that galanin can reduce seizure severity while increasing seizure occurrence, speculated to occur through different receptor subtypes. This study confirms galanin's complex role in brain activity, supporting its potential impact on epilepsy.

Strengths:

The overall strength of the study lies primarily in its methodological approach using whole-brain Calcium imaging facilitated by the transparency of zebrafish larvae. Additionally, the use of transgenic zebrafish models is an advantage, as it enables genetic manipulations to investigate specific aspects of galanin signaling. This combination of advanced imaging and genetic tools allows for addressing galanin's role in regulating brain activity.

Weaknesses:

The weaknesses of the study also stem from the methodological approach, particularly the use of whole-brain Calcium imaging as a measure of brain activity. While epilepsy and seizures involve network interactions, they typically do not originate across the entire brain simultaneously. Seizures often begin in specific regions or even within specific populations of neurons within those regions. Therefore, a whole-brain approach, especially with Calcium imaging with inherited limitations, may not fully capture the localized nature of seizure initiation and propagation, potentially limiting the understanding of Galanin's role in epilepsy.

Furthermore, Galanin's effects may vary across different brain areas, likely influenced by the predominant receptor types expressed in those regions. Additionally, the use of PTZ as a "stressor" is questionable since PTZ induces seizures rather than conventional stress. Referring to seizures induced by PTZ as "stress" might be a misinterpretation intended to fit the proposed model of stress regulation by receptors other than Galanin receptor 1 (GalR1).

The description of the EAAT2 mutants is missing crucial details. EAAT2 plays a significant role in the uptake of glutamate from the synaptic cleft, thereby regulating excitatory neurotransmission and preventing excitotoxicity. Authors suggest that in EAAT2 knockout (KO) mice galanin expression is upregulated 15-fold compared to wild-type (WT) mice, which could be interpreted as galanin playing a role in the hypoactivity observed in these animals.

However, the study does not explore the misregulation of other genes that could be contributing to the observed phenotype. For instance, if AMPA receptors are significantly downregulated, or if there are alterations in other genes critical for brain activity, these changes could be more important than the upregulation of galanin. The lack of wider gene expression analysis leaves open the possibility that the observed hypoactivity could be due to factors other than, or in addition to, galanin upregulation.

Moreover, the observation that in double KO mice for both EAAT2 and galanin there was little difference in seizure susceptibility compared to EAAT2 KO mice alone further supports the idea that galanin upregulation might not be the reason to the observed phenotype. This indicates that other regulatory mechanisms or gene expressions might be playing a more pivotal role in the manifestation of hypoactivity in EAAT2 mutants.

These methodological shortcomings and conceptual inconsistencies undermine the perceived strengths of the study, and hinders understanding of Galanin's role in epilepsy and stress regulation.

Comments on revisions:

The revised manuscript and the answers of the authors is appreciated. However, the criticisms were addressed only partially and main weaknesses of the manuscript are still remaining.

Reviewer #2 (Public review):

This revised study is an investigation of galanin and galanin receptor signaling on whole-brain activity in the context of recurrent seizure activity or under homeostatic basal conditions. The authors primarily use calcium imaging to observe whole-brain neuronal activity accompanied by galanin qPCR to determine how manipulations of galanin or the galr1a receptor affect the activity of the whole-brain under non-ictal conditions or when seizure activity occurs. The authors use their eaat2a-/- model (introduced in their Glia 2022 paper, PMID 34716961) that shows recurrent seizure activity as well as suppression of neuronal activity and locomotion interictally. It is compared to the well-known pentylenetetrazole (PTZ) pharmacological model of seizures in zebrafish. Given the literature cited in their Introduction, the authors hypothesize that galanin will exert a net inhibitory effect on brain activity in models of seizures/epilepsy. They were surprised to find that this hypothesis was only moderately supported in their eaat2a-/- model. In contrast, after PTZ, fish with galanin overexpression showed increased seizure number and reduced duration while fish with galanin KO showed reduced seizure number and increased duration.

Previous concerns about sex or developmental biological variables were addressed, as their model's seizure phenotype emerges rapidly and long prior to the establishment of zebrafish sexual maturity. However, in the course of re-review, some additional concerns (below) were detected that, if addressed, could further improve the manuscript. These concerns relate to how seizures were defined from the measurement of fluorescent calcium imaging data. Overall, this study is important and convincing, and carries clear value for understanding the multifaceted functions that neuronal galanin can perform under homeostatic and disease conditions.

Additional Concerns:

- The authors have validated their ability to measure behavioral seizures quantitatively in their 2022 Glia paper but the information provided on defining behavioral seizures was limited. The definition of behavioral seizure activity is not expanded upon in this paper, but could provide detail about how the behavioral seizures relate to a seizure detected via calcium imaging.

- Related to the previous point, for the calcium imaging, the difference between an increase in fluorescence that the authors think reflects increased neuronal activity and the fluorescence that corresponds to seizures is not very clear. This detail is necessary because exactly when the term "seizure" describes a degree of increased activity can be difficult to distinguish objectively.

- The supplementary movies that were added were very useful, but raised some questions. For example, what brain regions were pulsating? What areas seemed to constantly exhibit strong fluorescence and was this an artifact? It seemed that sometimes there was background fluorescence in the body. Perhaps an anatomical diagram could be provided for the readers. In addition, there were some movies with much greater fluorescence changes - are these the seizures? These are some reasons for our request for clarified definitions of the term "seizure".

Reviewer #3 (Public review):

Summary:

The neuropeptide galanin is primarily expressed in the hypothalamus and has been shown to play critical roles in homeostatic functions such as arousal, sleep, stress, and brain disorders such as epilepsy. Previous work in rodents using galanin analogs and receptor-specific knockout have provided convincing evidence for anti-convulsant effects of galanin.

In the present study, the authors sought to determine the relationship between galanin expression and whole-brain activity. The authors took advantage of the transparent nature of larval zebrafish to perform whole-brain neural activity measurements via widefield calcium imaging. Two models of seizures were used (eaat2a-/- and pentylenetetrazol; PTZ). In the eaat2a-/- model, spontaneous seizures occur and the authors found that galanin transcript levels were significantly increased and associated with reduced frequency of calcium events. Similarly, two hours after PTZ galanin transcript levels roughly doubled and the frequency and amplitude of calcium events were reduced.

The authors also used a heat shock protein line (hsp70I:gal) where galanin transcripts levels are induced by activation of heat shock protein, but this line also shows higher basal transcript levels of galanin. Due to problems with whole-brain activity in wild-type larvae, the authors used the line without heat shock. They found higher level of galanin in hsp70I:gal larval zebrafish resulted in a reduction of calcium events and a reduction in amplitude of events. In contrast, galanin knockout (gal-/-) increased calcium activity, indicated by an increased number of calcium events, but a reduction in amplitude and duration. New data in the supplementary figure 2 used antibody staining to confirm the absence of galanin expression in gal-/- knockouts. Knockout of the galanin receptor subtype galr1a via crispants also increased the frequency of calcium events. New data in the revised manuscript reports that galr1aKO did not cause an upregulation of galanin, thereby ruling out genetic compensation effects.

In subsequent experiments in eaat2a-/- mutants were crossed with hsp70I:gal or gal-/- to increase or decrease galanin expression, respectively. These experiments showed modest effects, with eaat2a-/- x gal-/- knockouts showing an increased normalized area under the curve and seizure amplitude.

Lastly, the authors attempted to study the relationship between galanin and brain activity during a PTZ challenge. The hsp70I:gal larva showed increased number of seizures and reduced seizure duration during PTZ. In contrast, gal-/- mutants showed increased normalized area under the curve and a stark reduction in number of detected seizures, a reduction in seizure amplitude, but an increase in seizure duration. The authors then ruled out the role of Galr1a in modulating this effect during PTZ, since the number of seizures was unaffected, whereas the amplitude and duration of seizures was increased.

Strengths:

(1) The gain- and loss-of function galanin manipulations provided convincing evidence that galanin influences brain activity (via calcium imaging) during interictal and/or seizure-free periods. In particular, the relationship between galanin transcript levels and brain activity in figures 1 & 2 was convincing. New antibody staining confirms the absence of galanin in gal-/- mutants. New data also shows galanin transcript levels were unchanged in galr1ako brains.

(2) The authors use two models of epilepsy (eaat2a-/- and PTZ).

(3) Focus on the galanin receptor subtype galr1a provided good evidence for an important role of this receptor in controlling brain activity during interictal and/or seizure-free periods.

(4) The authors have added supplementary video files for calcium imaging to support their observations.

Weaknesses:

(1) Although the relationship between galanin and brain activity during interictal or seizure-free periods was clear, the revised manuscript still lacks mechanistic insight in the role of galanin during seizure-like activity induced by PTZ.

(2) The revised manuscript continues to heavily rely on calcium imaging of different mutant lines. Confirmation of knockouts has been provided with immunostaining in a new supplementary figure. Additional methods could strengthen the data, translational relevance, and interpretation (e.g., acute pharmacology using galanin agonists or antagonists, brain or cell recordings, biochemistry, etc).

Reviewer #1 (Public review):

Summary:

In this manuscript, the authors addressed the previous comments from reviewers.

Strengths:

This study identified that NOLC1 could bind to p53 and decrease its nuclear transcriptional activity, then inhibit p53-mediated ferroptosis in gastric cancer.

Weaknesses:

There are a few Western blot images that were processed with excessive contrast adjustment, such as Figure 2I (Caspase-3 in MKN-45 group), Figure 4H (GPX4 in MKN-45 group), and Figure 5G/5I.

Reviewer #2 (Public review):

Summary:

Shengsheng Zhao et al. investigated the role of nucleolar and coiled-body phosphoprotein 1 (NOLC1) in relegating gastric cancer (GC) development and cisplatin-induced drug resistance in GC. They found a significant correlation between high NOLC1 expression and the poor prognosis of GC. Meanwhile, upregulation of NOLC1 was associated with cis-resistant GC. Experimentally, the authors demonstrate that knocking down NOLC1 increased GC sensitivity to Cis possibly by regulating ferroptosis. Mechanistically, they found NOLC1 suppressed ferroptosis by blocking the translocation of P53 from the cytoplasm to the nucleus and promoting its degradation. In addition, the authors also evaluated the effect of combinational treatment of anti-PD-1 and cisplatin in NOLC1 -knockdown tumor cells, revealing a potential role of NOLC1 in the targeted therapy for GC.

Strengths:

Chemoresistance is considered a major reason causing failure of tumor treatment and death of cancer patients. This paper explored the role of NOLC1 in the regulation of Cis-mediated resistance, which involves a regulated cell death named ferroptosis. These findings provide more evidence highlighting the study of regulated cell death to overcome drug resistance in cancer treatment, which could give us more potential strategies or targets for combating cancer.

Weaknesses:

More evidence supporting the regulation of ferroptosis induced by Cisplatin by NOLC1 should be added. Particularly, the role of ferroptosis in the cisplatin-resistance should be verified and whether NOLC1 regulates ferroptosis induced by additional FINs should be explored. Besides, the experiments to verify the regulation of ferroptosis sensitivity by NOLC1 are sort of superficial. The role of MDM2/p53 in ferroptosis or cisplatin resistance mediated by NOLC1 should be further studied by genetic manipulation of p53, which is the key evidence to confirm its contribution to NOLC1 regulation of GC and relative cell death.

Reviewer #1 (Public review):

Summary:

Gene transfer agent (GTA) from Bartonella is a fascinating chimeric GTA that evolved from the domestication of two phages. Not much is known about how the expression of the BaGTA is regulated. In this manuscript, Korotaev et al noted the structural similarity between BrrG (a protein encoded by the ror locus of BaGTA) to a well-known transcriptional anti-termination factor, 21Q, from phage P21. This sparked the investigation into the possibility that BaGTA cluster is also regulated by anti-termination. Using a suite of cell biology, genetics, and genome-wide techniques (ChIP-seq), Korotaev et al convincingly showed that this is most likely the case. The findings offer the first insight into the regulation of GTA cluster (and GTA-mediated gene transfer) particularly in this pathogen Bartonella. Note that anti-termination is a well-known/studied mechanism of transcriptional control. Anti-termination is a very common mechanism for gene expression control of prophages, phages, bacterial gene clusters, and other GTAs, so in this sense, the impact of the findings in this study here is limited to Bartonella.

Strengths:

convincing results that overall support the main claim of the manuscript.

Weaknesses:

A few important controls are missing.

Comments on revisions:

I am happy with this revised version except for one point, that is a single replicate for ChIP-seq, I don't think that is appropriate.

Reviewer #2 (Public review):

Summary:

In this study, the authors identified and characterized a regulatory mechanism based on transcriptional anti-termination that connects the two gene clusters, capsid and run-off replication (ROR) locus, of the bipartite Bartonella gene transfer agent (GTA). Among genes essential for GTA functionality identified in a previous transposon sequencing project, they found a potential antiterminatior of phage origin within the ROR locus. They employed fluorescence reporter and gene transfer assays of overexpression and knockout strains in combination with ChiPSeq and promoter-fusions to convincingly show that this protein indeed acts as an antiterminator counteracting attenuation of the capsid gene cluster expression.

Impact on the field:

The results provide valuable insights into the evolution of the chimeric BaGTA, a unique example of phage co-domestication by bacteria. A similar system found in the other broadly studied Rhodobacterales/Caulobacterales GTA family suggests that antitermination could be a general mechanism for GTA control.

Strengths:

Results of the selected and carefully designed experiments support the main conclusions.

Weaknesses:

The question why overexpression of the antiterminator does not increase the gene tranfer frequency needs to be answered in further studies.

Comments on revisions:

The authors further improved the already strong manuscript. All my concerns have been addressed. The addition of a summry figure helps to understand the proposed mechanism.

Reviewer #1 (Public review):

Koesters and colleagues investigated the role of the small GTPase Rab3A in homeostatic scaling of miniature synaptic transmission in primary mouse cortical cultures using electrophysiology and immunohistochemistry. The major finding is that TTX incubation for 48 hours does not induce an increase in the amplitude of excitatory synaptic miniature events in neuronal cortical cultures derived from Rab3A KO and Rab3A Earlybird mutant mice. NASPM application had comparable effects on mEPSC amplitude in control and after TTX, implying that Ca2+-permeable glutamate receptors are unlikely modulated during synaptic scaling. Immunohistochemical analysis revealed no significant changes in GluA2 puncta size, intensity, and integral after TTX treatment in control and Rab3A KO cultures. Finally, they provide evidence that loss of Rab3A in neurons, but not astrocytes, blocks homeostatic scaling. Based on these data, the authors propose a model in which neuronal Rab3A is required for homeostatic scaling of synaptic transmission, potentially through GluA2-independent mechanisms.

The major finding - impaired homeostatic up-scaling after TTX treatment in Rab3A KO and Rab3 earlybird mutant neurons - is supported by data of high quality. However, the paper falls short of providing any evidence or direction regarding potential mechanisms. The data on GluA2 modulation after TTX incubation are likely statistically underpowered and do not allow drawing solid conclusions, such as GluA2-independent mechanisms of up-scaling.

The study should be of interest to the field because it implicates a presynaptic molecule in homeostatic scaling, which is generally thought to involve postsynaptic neurotransmitter receptor modulation. However, it remains unclear how Rab3A participates in homeostatic plasticity.

Major (remaining) point:

(1) The current version of the abstract only includes the results on GluA2 immunofluorescence and mEPSC amplitude modulation after TTX treatment in control cultures, and a requirement for Rab3A in neurons instead of astrocytes. The major findings, including the block of the mEPSC amplitude increase upon TTX treatment in Rab3KO/EB mutants, are not mentioned. The abstract should be revised so that it reflects all major findings, potentially at the expense of citing previous work by the authors.

Reviewer #2 (Public review):

First, I would like to thank the authors for the response. I acknowledge that the authors show in previous studies that Rab3A acts from the presynaptic side at the NMJ, and that is, as the authors indicate, their impetus for the current study. However, mechanisms observed at a completely different type of synapses cannot be used as an argument for conclusions here. The authors also acknowledge that they should restrict their conclusions to the data in the current study, and they are merely proposing interpretations. Then perhaps they should restrict these interpretations to the discussion rather than make this claim in the abstract (lines 44-47). Here the authors ask whether Rab3A is involved in the homeostatic increase of postsynaptic AMPARs, am I understanding it correctly that their conclusion for this question is "increase in AMPAR levels in WT cultures is more variable than those in mEPSCs so that it is impossible to determine if Rab3A is involved at all"? If so, then this question has not been answered and should not be regarded as one of the main conclusions with the data presented here. It also remains unclear to me how this piece of inconclusive data serves the main objective of the study.

The authors state at the end that the current study is just an extension of their previous work, and therefore their interpretations here further support the idea that Rab3A is acting presynaptically. I would argue that it is the conclusive data, rather than interpretations that lack concrete evidence, that support ideas and models. I think that we would all agree that immunostaining measurements can be very variable. However, if the authors are determined to use this approach to answer one of their major questions, then perhaps one way to significantly strengthen their conclusions is to find ways to somewhat overcome this technical limitation.

Finally, I thank the authors for addressing other minor concerns of mine.

Reviewer #3 (Public review):

This manuscript presents a number of interesting findings that have the potential to increase our understanding of the mechanism underlying homeostatic synaptic plasticity (HSP). The data broadly support that Rab3A plays a role in HSP, although the site and mechanism of action remain uncertain.

The authors clearly demonstrate the Rab3A plays a role in HSP at excitatory synapses, with substantially less plasticity occurring in the Rab3A KO neurons. There is also no apparent HSP in the Earlybird Rab3A mutation, although baseline synaptic strength is already elevated. In this context, it is unclear if the plasticity is absent, already induced by this mutation, or just occluded by a ceiling effect due the synapses already being strengthened. Occlusion may also occur in the mixed cultures, when Rab3A is missing from neurons but not astrocytes. The authors do appropriately discuss these options. The authors have solid data showing that Rab3A is unlikely to be active in astrocytes, Finally, they attempt to study the linkage between changes in synaptic strength and AMPA receptor trafficking during HSP, and conclude that trafficking may not be solely responsible for the changes in synaptic strength during HSP.

Strengths:

This work adds another player into the mechanisms underlying an important form of synaptic plasticity. The plasticity is likely only reduced, suggesting Rab3A is only partially required and perhaps multiple mechanisms contribute. The authors speculate about some possible novel mechanisms, including whether Rab3A is active pre-synaptically to regulate quantal amplitude.

As Rab3A is primarily known as a pre-synaptic molecule, this possibility is intriguing and novel for this system. However, it is based on the partial dissociation of AMPAR trafficking and synaptic response, and lacks strong support. On average, they saw similar magnitude of change in mEPSC amplitude and GluA2 cluster area and integral, but the GluA2 data was not significant due to higher variability. It is difficult to determine if this is due to biology or methodology - the imaging method involves assessing puncta pairs (GluA2/VGlut1) clearly associated with a MAP2 labeled dendrite. This is a small subset of synapses, with usually less than 20 synapses per neuron analyzed, which would be expected to be more variable than mEPSC recordings averaged across several hundred events. However, when they reduce the mEPSC number of events to similar numbers as the imaging, the mESPC amplitudes are still less variable than the imaging data. The reason for this remains unclear. The pool of sampled synapses is still different between the methods and recent data has shown that synapses have variable responses during HSP. Further, there could be variability in the subunit composition of newly inserted AMPARs, and only assessing GluA2 could mask this (see below). It is intriguing that pre-synaptic changes might contribute to HSP, especially given the likely localization of Rab3A. But it remains difficult to distinguish if the apparent difference in imaging and electrophysiology is a methodological issue rather than a biological one. Stronger data, especially positive data on changes in release, will be necessary to conclude that pre-synaptic factors are required for HSP, beyond the established changes in post-synaptic receptor trafficking. Specific deletion of Rab3A from pre-synaptic neurons would also be highly informative.

Other questions arise from the NASPM experiments, used to justify looking at GluA2 (and not GluA1) in the immunostaining. First, there is a strong frequency effect that is unclear in origin. One would expect NASPM to merely block some fraction of the post-synaptic current, and not affect pre-synaptic release or block whole synapses. But the change in frequency seems to argue (as the authors do) that some synapses only have CP-AMPARs, while the rest of the synapses have few or none. Another possibility is that there are pre-synaptic NASPM-sensitive receptors that influence release probability. Further, the amplitude data show a strong trend towards smaller amplitude following NASPM treatment (Fig 3B). The p value for both control and TTX neurons was 0.08 - it is very difficult to argue that there is no effect. And the decrease on average is larger in the TTX neurons, and some cells show a strong effect. It is possible there is some heterogeneity between neurons on whether GluA1/A2 heteromers or GluA1 homomers are added during HSP. This would impact the conclusions about the GluA2 imaging as compared to the mEPSC amplitude data.

To understand the role of Rab3A in HSP will require addressing two main issues:

(1) Is Rab3A acting pre-synaptically, post-synaptically or both? The authors provide good evidence that Rab3A is acting within neurons and not astrocytes. But where it is acting (pre or post) would aid substantially in understanding its role. The general view in the field has been that HSP is regulated post-synaptically via regulation of AMPAR trafficking, and considerable evidence supports this view. More concrete support for the authors suggestion of a pre-synaptic site of control would be helpful.

(2) Rab3A is also found at inhibitory synapses. It would be very informative to know if HSP at inhibitory synapses is similarly affected. This is particularly relevant as at inhibitory synapses, one expects a removal of GABARs or a decrease in GABA release (ie the opposite of whatever is happening at excitatory synapses). If both processes are regulated by Rab3A, this might suggest a role for this protein more upstream in the signaling; an effect only at excitatory synapses would argue for a more specific role just at those synapses.

Comments on revisions:

The section on TNF is a bit odd. The data on the astrocyte deletion of Rab3A only argues that Rab3A is unlikely to regulate TNF release. But it could easily be downstream of the neuronal TNF receptor. Without any data addressing the TNF response, it seems quite premature to argue that Rab3A is part of a TNF-independent pathway.

The section title (line 506-7) declaring Rab3A as the first presynaptic protein involved in HSP is also premature, as they don't know it is acting pre-synaptically.

Reviewer #1 (Public review):

The authors of this study use electron microscopy and 3D reconstruction techniques to study the morphology of distinct classes of Drosophila sensory neurons *across many neurons of the same class.* This is a comprehensive study attempting to look at nearly all the sensory neurons across multiple sensilla to determine a) how much morphological variability exists between and within neurons of different and similar sensory classes, and 2) identify dendritic features that may have evolved to support particular sensory functions. This study builds upon the authors' previous work, which allowed them to identify and distinguish sensory neuron subtypes in the EM volumes without additional staining so that reconstructed neurons could reliably be placed in the appropriate class. This work is unique in looking at a large number of individual neurons of the same class to determine what is consistent and what is variable about their class-specific morphologies.

This means that in addition to providing specific structural information about these particular cells, the authors explore broader questions of how much morphological diversity exists between sensory neurons of the same class and how different dendritic morphologies might affect sensory and physiological properties of neurons.

The authors found that CO2-sensing neurons have an unusual, sheet-like morphology in contrast to the thin branches of odor-sensing neurons. They show that this morphology greatly increases the surface area to volume ratio above what could be achieved by modest branching of thin dendrites, and posit that this might be important for their sensory function, though this was not directly tested in their study. The study is mainly descriptive in nature, but thorough, and provides a nice jumping-off point for future functional studies. One interesting future analysis could be to examine all four cell types within a single sensilla together to see if there are any general correlations that could reveal insights about how morphology is determined and the relative contributions of intrinsic mechanisms vs interactions with neighboring cells. For example, if higher than average branching in one cell type correlated with higher than average branching in another type, if in the same sensilla. This might suggest higher extracellular growth or branching cues within a sensilla. Conversely, if higher branching in one cell type consistently leads to reduced length or branching in another, this might point to dendrite-dendrite interactions between cells undergoing competitive or repulsive interactions to define territories within each sensilla as a major determinant of the variability.

Reviewer #2 (Public review):

Summary:

The manuscript employs serial block‐face electron microscopy (SBEM) and cryofixation to obtain high‐resolution, three‐dimensional reconstructions of Drosophila antennal sensilla containing olfactory receptor neurons (ORNs) that detect CO2. This method has been used previously by the same lab in Gonzales et. al, 2021. (https://elifesciences.org/articles/69896), which had provided an exemplary model by integrating high-resolution EM with electrophysiology and cell-type-specific labeling. The previous study ended up correlating morphology with activity for multiple olfactory sensillar types. Compared to the 2021 study, this current manuscript appears somewhat incomplete and lacks integration with activity.

In fact older studies have also reported two-dimensional TEM images of the putative CO2 neuron in Drosophila (Shanbhag et al., 1999) and in mosquitoes (McIver and Siemicki, 1975; Lu et al, 2007), and in these instances reported that the dendritic architecture of the CO2 neuron was somewhat different (circular and flattened, lamellated) from other olfactory neurons.

The authors claim that this approach offers an artifact‐minimized ultrastructural dataset compared to earlier. In this study, not only do they confirm this different morphology but also classify it into distinct subtypes (loosely curled, fully curled, split, and mixed). This detailed morphological categorization was not provided in prior studies (e.g., Shanbhag et al., 1999 ). The authors would benefit from providing quantitative thresholds or objective metrics to improve reproducibility and to clarify whether these structural distinctions correlate with distinct functional roles.

Strengths:

The study makes a convincing case that ab1C neurons exhibit a unique, flattened dendritic morphology unlike the cylindrical dendrites found in ab1D neurons. This observation extends previous qualitative TEM findings by not only confirming the presence of flattened lamellae in CO₂ neurons but also quantifying key morphometrics such as dendritic length, surface area, and volume, and calculating surface area-to-volume ratios. The enhanced ratios observed in the flattened segments are speculated to be linked to potential advantages in receptor distribution (e.g., Gr21a/Gr63a) and efficient signal propagation.

Weaknesses:

While the manuscript offers valuable ultrastructural insights and reveals previously unappreciated heterogeneity among CO₂-sensing neurons, several issues warrant further investigation in addition to the points made above.

(1) Although this quantitative approach is robust compared to earlier descriptive reports, its impact is somewhat limited by the absence of direct electrophysiological data to confirm that ultrastructural differences translate into altered neuronal function. A direct comparison or discussion of how the present findings align with the functional data obtained from electrophysiology would strengthen the overall argument.

(2) Clarifying the criteria for dendritic subtype classification with quantitative parameters would enhance reproducibility and interpretability. Moreover, incorporating electrophysiological recordings from ab1C neurons would provide compelling evidence linking structure and function, and mapping key receptor proteins through immunolabeling could directly correlate receptor distribution with the observed morphological diversity.

(3) Even though Cryofixation is claimed to be superior to chemical fixation for generating fewer artifacts, authors need to confirm independently the variation observed in the CO2 neuron morphologies across populations. All types of fixation in TEMs cause some artifacts, as does serial sectioning. Without understanding the error rates or without independent validation with another method, it is hard to have confidence in the conclusions drawn by the authors of the paper.

Addressing these concerns and integrating additional experiments would significantly bolster the manuscript's completeness and advancement.

Reviewer #3 (Public review):

Summary:

In the current manuscript entitled "Population-level morphological analysis of paired CO2- and odor-sensing olfactory neurons in D. melanogaster via volume electron microscopy", Choy, Charara et al. use volume electron microscopy and neuron reconstruction to compare the dendritic morphology of ab1C and ab1D neurons of the Drosophila basiconic ab1 sensillum. They aim to investigate the degree of dendritic heterogeneity within a functional class of neurons using ab1C and ab1D, which they can identify due to the unique feature of ab1 sensilla to house four neurons and the stereotypic location on the third antennal segment. This is a great use of volumetric electron imaging and neuron reconstruction to sample a population of neurons of the same type. Their data convincingly shows that there is dendritic heterogeneity in both investigated populations, and their sample size is sufficient to strongly support this observation. This data proposes that the phenomenon of dendritic heterogenity is common in the Drosophila olfactory system and will stimulate future investigations into the developmental origin, functional implications, and potential adaptive advantage of this feature.

Moreover, the authors discovered that there is a difference between CO2- and odour-sensing neurons of which the first show a characteristic flattened and sheet-like structure not observed in other sensory neurons sampled in this and previous studies. They hypothesize that this unique dendritic organization, which increases the surface area to volume ratio, might allow more efficient Co2 sensing by housing higher numbers of Co2 receptors. This is supported by previous attempts to express Co2 sensors in olfactory sensory neurons, which lack this dendritic morphology, resulting in lower Co2 sensitivity compared to endogenous neurons.

Overall, this detailed morphological description of olfactory sensory neurons' dendrites convincingly shows heterogeneity in two neuron classes with potential functional impacts for odour sensing.

Strength:

The volumetric EM imaging and reconstruction approach offers unprecedented details in single cell morphology and compares dendrite heterogeneity across a great fraction of ab1 sensilla.<br /> The authors identify specific shapes for ab1C sensilla potentially linked to their unique function in CO2 sensing.

Weaknesses:

While the morphological description is highly detailed, no attempts are made to link this to odour sensitivity or other properties of the neurons. It would have been exciting to see how altered morphology impacts physiology in these olfactory sensory cells.

Reviewer #1 (Public review):

In this study, Marocco and colleagues perform a deep characterization of the complex molecular mechanism guiding the recognition of a particular CELLmotif previously identified in hepatocytes in another publication. Having miR-155-3p with or without this CELLmotif as initial focus, the authors identify 21 proteins differentially binding to these two miRNA versions. From these, they decided to focus on PCBP2. They elegantly demonstrate PCBP2 binding to miR-155-3p WT version but not to the CELLmotif-mutated version. miR-155-3p contains a hEXOmotif identified in a different report, whose recognition is largely mediated by another RNA-binding protein called SYNCRIP. Interestingly, mutation of the hEXOmotif contained in miR-155-3p did not only blunt SYNCRIP binding, but also PCBP2 binding despite the maintenance of the CELLmotif. This indicates that somehow SYNCRIP binding is a prerequisite for PCBP2 binding. EMSA assay confirms that SYNCRIP is necessary for PCBP2 binding to miR-155-3p, while PCBP2 is not needed for SYNCRIP binding. The authors aim to extend these findings to other miRNAs containing both motifs. For that, they perform a small-RNA-Seq of EVs released from cells knockdown for PCBP2 versus control cells, identifying a subset of miRNAs whose expression either increases or decreases. The assumption is that those miRNAs containing PCBP2-binding CELLmotif should now be less retained in the cell and go more to extracellular vesicles, thus reflecting a higher EV expression. The specific subset of miRNAs having both the CELLmotif and hEXOmotif (9 miRNAs) whose expressions increase in EVs due to PCBP2 reduction is also affected by knocking down SYNCRIP in the sense that reduction of SYNCRIP leads to lower EV sorting. Further experiments confirm that PCBP2 and SYNCRIP bind to these 9 miRNAs and that knocking down SYNCRIP impairs their EV sorting.

Reviewer #2 (Public review):

Summary:

The author of this manuscript aimed to uncover the mechanisms behind miRNA retention within cells. They identified PCBP2 as a crucial factor in this process, revealing a novel role for RNA-binding proteins. Additionally, the study discovered that SYNCRIP is essential for PCBP2's function, demonstrating the cooperative interaction between these two proteins. This research not only sheds light on the intricate dynamics of miRNA retention but also emphasizes the importance of protein interactions in regulating miRNA behavior within cells.

Strengths:

This paper makes important progress in understanding how miRNAs are kept inside cells. It identifies PCBP2 as a key player in this process, showing a new role for proteins that bind RNA. The study also finds that SYNCRIP is needed for PCBP2 to work, highlighting how these proteins work together. These discoveries not only improve our knowledge of miRNA behavior but also suggest new ways to develop treatments by controlling miRNA locations to influence cell communication in diseases. The use of liver cell models and thorough experiments ensures the results are reliable and show their potential for RNA-based therapies.

Reviewer #1 (Public review):

Summary:

This study uncovers a protective role of the ubiquitin-conjugating enzyme variant Uev1A in mitigating cell death caused by over-expressed oncogenic Ras in polyploid Drosophila nurse cells and by RasK12 in diploid human tumor cell lines. The authors previously showed that overexpression of oncogenic Ras induces death in nurse cells, and now they perform a deficiency screen for modifiers. They identified Uev1A as a suppressor of this Ras-induced cell death. Using genetics and biochemistry, the authors found that Uev1A collaborates with the APC/C E3 ubiquitin ligase complex to promote proteasomal degradation of Cyclin A. This function of Uev1A appears to extend to diploid cells, where its human homologs UBE2V1 and UBE2V2 suppress oncogenic Ras-dependent phenotypes in human colorectal cancer cells in vitro and in xenografts in mice.

Strengths:

(1) Most of the data is supported by a sufficient sample size and appropriate statistics.<br /> (2) Good mix of genetics and biochemistry.<br /> (3) Generation of new transgenes and Drosophila alleles that will be beneficial for the community.

Weaknesses:

(1) Phenotypes are based on artificial overexpression. It is not clear whether these results are relevant to normal physiology.

(2) The phenotype of "degenerating ovaries" is very broad, and the study is not focused on phenotypes at the cellular level. Furthermore, no information is provided in the Materials and Methods on how degenerating ovaries are scored, despite this being the most important assay in the study.

(3) In Figure 5, the authors want to conclude that uev1a is a tumor-suppressor, and so they over-express ubev1/2 in human cancer cell lines that have RasK12 and find reduced proliferation, colony formation, and xenograft size. However, genes that act as tumor suppressors have loss-of-function phenotypes that allow for increased cell division. The Drosophila uev1a mutant is viable and fertile, suggesting that it is not a tumor suppressor in flies. Additionally, they do not deplete human ubev1/2 from human cancer cell lines and assess whether this increases cell division, colony formation, and xenograph growth.

(4) A critical part of the model does not make sense. CycA is a key part of their model, but they do not show CycA protein expression in WT egg chambers or in their over-expression models (nos.RasV12 or bam>RasV12). Based on Lilly and Spradling 1996, Cyclin A is not expressed in germ cells in region 2-3 of the germarium; whether CycA is expressed in nurse cells in later egg chambers is not shown but is critical to document comprehensively.

(5) The authors should provide more information about the knowledge base of uev1a and its homologs in the introduction.

Reviewer #2 (Public review):

Summary:

The authors performed a genetic screen using deficiency lines and identified Uev1a as a factor that protects nurse cells from RasG12V-induced cell death. According to a previous study from the same lab, this cell death is caused by aberrant mitotic stress due to CycA upregulation (Zhang et al.). This paper further reveals that Uev1a forms a complex with APC/C to promote proteasome-mediated degradation of CycA.

In addition to polyploid nurse cells, the authors also examined the effect of RasG12V-overexpression in diploid germline cells, where RasG12V-overexpression triggers active proliferation, not cell death. Uev1a was found to suppress its overgrowth as well.

Finally, the authors show that the overexpression of the human homologs, UBE2V1 and UBE2V2, suppresses tumor growth in human colorectal cancer xenografts and cell lines. Notably, the expression of these genes correlates with the survival of colorectal cancer patients carrying the Ras mutation.

Strength:

This paper presents a significant finding that UBE2V1/2 may serve as a potential therapy for cancers harboring Ras mutations. The authors propose a fascinating mechanism in which Uev1a forms a complex with APC/C to inhibit aberrant cell cycle progression.

Weakness:

The quantification of some crucial experiments lacks sufficient clarity.

Reviewer #1 (Public review):

Summary:

The authors confirmed earlier findings that AVP influences α and β cells differently, depending on glucose concentrations. At substimulatory glucose levels, AVP combined with forskolin - an activator of cAMP -did not significantly stimulate β cells, although it did activate α cells. Once glucose was raised to stimulatory levels, β cells became active, and α cell activity declined, indicating glucose's suppressive effect on α cells and permissive effect on β cells. Under physiological glucose levels (8-9 mM), forskolin enhanced β-cell calcium oscillations, and AVP further modulated this activity. However, AVP's effect on β cells was variable across islets and did not significantly alter AUC measurements (a combined indicator of oscillation frequency and duration). In α cells, forskolin and AVP led to increased activity even at high glucose levels, suggesting that α cells remain responsive despite expected suppression by insulin and glucose.

Experiments with physiological concentrations of epinephrine suggest that AVP does not operate via Gs-coupled V2 receptors in β cells, as AVP could not counteract epinephrine's inhibitory effects. Instead, epinephrine reduced β cell activity while increasing α cell activity through different G-protein-coupled mechanisms. These results emphasize that AVP can potentiate α-cell activation and has a nuanced, context-dependent effect on β cells.

The most robust activation of both α and β cells by AVP occurred within its physiological osmo-regulatory range (~10-100 pM), confirming that AVP exerts bell-shaped concentration-dependent effects on β cells. At low concentrations, AVP increased β cell calcium oscillation frequency and reduced "halfwidths"; high concentrations eventually suppressed β cell activity, mimicking the muscarinic signaling. In α cells, higher AVP concentrations were required for peak activation, which was not blunted by receptor inactivation within physiological ranges.

Attempting to further dissect the role of specific AVP receptors, the authors designed and tested peptide ligands selective for V1b receptors. These included a selective V1b agonist; a V1b agonist with antagonist properties at V1a and oxytocin receptors; and a selective V1a antagonist. In pancreatic slices, these peptides seem to replicate AVP's effects on Ca²⁺ signaling, although responses were highly variable, with some islets showing increased activity and others no change or suppression. The variability was partly attributed to islet-specific baseline activity, and the authors conclude that AVP and V1b receptor agonists can modulate β cell activity in a state-dependent manner, stimulating insulin secretion in quiescent cells and inhibiting it in already active cells.

Strengths:

Overall, the study is technically advanced and provides useful pharmacological tools. However, the conclusions are limited by a lack of direct mechanistic and functional data. Addressing these gaps through a combination of signaling pathway interrogation, functional hormone output, genetic validation, and receptor localization would strengthen the conclusions and reduce the current (interpretive) ambiguity.

Weaknesses:

(1) The study is entirely based on pharmacological tools. Without genetic models, off-target effects or incomplete specificity of the peptides cannot be fully ruled out.

(2) Despite multiple claims about β cell activation or inhibition, the functional output - insulin secretion - is weakly assessed, and only in limited conditions. This aspect makes it very hard to correlate calcium dynamics with physiological outcomes.

(3) Insulin and glucagon secretion assays should be provided; the authors should measure hormone release in parallel with Ca2+ imaging, using perifusion assays, especially during AVP ramp and peptide ligand applications.

Additionally, there is no standardization of the metabolic state of islets. The authors should consider measuring islet NAD(P)H autofluorescence or mitochondrial potential (e.g., using TMRE) to control for metabolic variability that may affect responsiveness.

(4) There is a high degree of variability in response to AVP and V1b agonists across islets (activation, no effect, inhibition). Surprisingly, the authors do not fully explore the cause of this heterogeneity (whether it is due to receptor expression differences, metabolic state, experimental variability, or other conditions).

(5) There is no validation of V1b receptor expression at the protein or mRNA level in α or β cells using in situ hybridization, immunohistochemistry, or spatial transcriptomics.

(6) AVP effects are described in terms of permissive or antagonistic effects on cAMP (especially in relation to epinephrine), but direct measurements of cAMP in α and β cells are not shown, weakening these conclusions. The authors should use Epac-based cAMP FRET sensors in α and β cells to monitor the interaction between AVP, forskolin, and epinephrine more conclusively.

(7) Single-islet transcriptomics or proteomics (also to clarify variability) should be provided to analyze receptor expression variability across islets to correlate with response phenotypes (activation vs inhibition). Alternatively, the authors could perform calcium imaging with simultaneous insulin granule tracking or ATP levels to assess islet functional states.

(8) While the study implies AVP acts through V1b receptors on β cells, the signaling downstream (e.g., PLC activation, IP3R isoforms involved) is simply inferred but not directly shown.

(9) The interpretation that IP3R inactivation (mentioned in the title!) underlies the bell-shaped AVP effect is just hypothetical, without direct measurements. Assays in β (and/or α)-cell-specific V1b KO mice and IP3R KO mice must be provided to support these speculations.

Reviewer #2 (Public review):

Summary:

In this paper, Drs. Kercmar, Murko, and Bombek make a series of observations related to the role of AVP in pancreatic islets. They use the pancreatic slice preparation that their group is well known for. The observations on the slide physiology are technically impressive. However, I am not convinced by the conclusions of this manuscript for a number of reasons. At the core of my concern is perhaps that this manuscript appears to be motivated to resolve 'controversies' surrounding the actions of AVP on insulin and glucagon secretion. This manuscript adds more observations, but these do not move the field forward in improving or solidifying our mechanistic understanding of AVP actions on islets. A major claim in this manuscript is the beta cell expression of the V1b Receptor for AVP, but the evidence presented in this paper falls short of supporting this claim. Observations on the activation of calcium in alpha cells via V1b receptor align with prior observations of this effect.

I have focused my main concerns below. I hope the authors will consider these suggestions carefully - please be assured that they were made with the intent to support the authors and increase the impact of this work.

Strengths:

The main strength of this paper is the technical sophistication of the approach and the analysis and representation of the calcium traces from alpha and beta cells.

Weaknesses:

(1) The introduction is long and summarizes a substantive body of literature on AVP actions on insulin secretion in vivo. There are a number of possible explanations for these observations that do not directly target islet cells. If the goal is to resolve the mechanistic basis of AVP action on alpha and beta cells, the more limited number of papers that describe direct islet effects is more helpful. There are excellent data that indicate that the actions of AVP are mediated via V1bR on alpha cells and that V1bR is a) not expressed by beta cells and b) does not activate beta cell calcium at all at 10 nM - which is the same concentration used in this paper (Figure 4G) for peak alpha cell Ca2+ activation (see https://doi.org/10.1016/j.cmet.2017.03.017; cited as ref 30 in the current manuscript).

(2) We know from bulk RNAseq data on purified alpha, beta, and delta cells from both the Huising and Gribble groups that there is no expression of V2a. I will point you to the data from the Huising lab website published almost a decade ago (http://dx.doi.org/10.1016/j.molmet.2016.04.007) - which is publicly available and can be used to generate figures (https://huisinglab.com/data-ghrelin-ucsc/index.html). They indicate the absence of expression of not only AVP2 receptors anywhere in the islet, but also the lack of expression of V1bra, V1brb, and Oxtr in beta cells. Instead of the detailed list of expression of these 4 receptors elsewhere in the body, it would be more directly relevant to set up their pancreatic slice experiments to summarize the known expression in pancreatic islets that is publicly available. It would also have helped ground the efforts that involved the generation of the V1aR agonist and V2R antagonist, which confirm these known AVP/OXT receptor expression patterns.

(3) Importantly, the lack of V1br from beta cells does not invalidate observations that AVP affects calcium in beta cells, but it does indicate that these effects are mediated a) indirectly, downstream of alpha cell V1br or b) via an unknown off-target mechanism (less likely). The different peak efficacies in Figure 4G would also suggest that they are not mediated by the same receptor.

(4) The rationale for the use of forskolin across almost all traces is unclear. It is motivated by a desire to 'study the AVP dependence of both alpha and beta cells at the same time'. As best as I can determine, the design choice to conduct all studies under sustained forskolin stimulation is related to the permissive actions of AVP on hormone secretion in response to cAMP-generating stimuli. The permissive actions by AVP that are cited are on hormone secretion, which in many cell types requires activation of both calcium and cAMP signaling. Whether the activation of V1br and subsequent calcium response is permitted by cAMP is unclear. I believe the argument the authors are making here is that the activation of beta cell calcium by AVP is permitted by forskolin. i.e., the cAMP stimulated by it in beta cells. However, the design does not account for the elevation of cAMP in alpha cells and subsequent release of glucagon, particularly upon co-stimulation with AVP, which permits glucagon release by activating a calcium response in alpha cells. This glucagon could then activate beta cells. If resolving the mechanism of action is the goal, often less is more. The activation of Gaq-mediated calcium is not cAMP dependent (although the downstream hormone secretion clearly often is). As was shown, AVP does not activate calcium in beta cells in the absence of cAMP. The experiments in Figures 1, 2, and 4 should have been completed in the absence of cAMP first.

(5) It is unexpected that epinephrine in Figure 2 does not activate the alpha cell calcium? A recent paper from the same group (Sluga et al) shows robust calcium activation in alpha cells in a similar prep by 1 nM epinephrine, which is similar to the dose used here.

(6) Figure 8 suggests a pharmacological activation of beta cell V1bR in the low pM range. How do the authors reconcile this comparison with the apparent absence of an effect of AVP stimulation at low pM to low nM doses in beta cells (Figure 4A)? I note that there are changes over time with sustained beta cell stimulation with 8 mM glucose, but these changes are relatively subtle, gradual, and quite likely represent the progression of calcium behaviors that would have occurred under sustained glucose, irrespective of these very low AVP concentrations. I will note that the Kd of the V1bR for AVP is around 1 nM, with tracer displacement starting around 100 pM according to the data in figure 5B, which is hard to reconcile with changes in beta cell calcium by AVP doses that start 10-100-fold lower than this dose at 1 and 10 pM (Figure 8).

Reviewer #3 (Public review):

Summary:

This work aims to better understand the role of arginine vasopressin (AVP) in the control of islet hormone secretion. This builds on previous literature in this area reporting on the actions of AVP to stimulate islet hormones. The gap in literature being addressed by these studies is primarily focused on the glucose-dependency of AVP on both insulin and glucagon secretion. A secondary objective is to explore the role of individual receptors with the use of newly generated peptides and existing tools. The methods include the use of Ca2+ imaging in pancreas slices from mice, with additional outcomes including insulin secretion in some areas. The conclusions presented are that AVP acts through V1b receptors in both alpha- and beta-cells, that this activity occurs in the high cAMP environment, and is glucose dependent.

Strengths:

The area of research is emerging with plenty of room for new contributions. The concept of AVP stimulating islet hormone secretion is important and deserving of further insight. The use of pancreas tissue to image primary cells makes the experiments physiologically relevant. The advancement of novel tools in this area should be helpful to other groups investigating the actions of AVP.

Weaknesses:

The conclusions are only modestly supported by the data and lack experimental depth and rigor. The rationale for only conducting studies at high cAMP conditions is not entirely clear and limits the conclusions that can be made. The use of Ca2+ is helpful, but it is a surrogate for hormone secretion. Additional measurements of hormone secretion are needed to enhance the robustness of these conclusions. Consideration of paracrine effects between alpha- and beta-cells is only superficially made and is likely essential in the context of the experimental design. For instance, there is clear literature that alpha-cells secrete several factors that work in paracrine interactions on beta-cells and autocrine actions back on alpha-cells. Conducting these studies in a high cAMP context only completely overlooks these interactions, skewing the interpretations made by the investigators. Finally, the clarity of the experiments and results could be significantly enhanced.

Reviewer #1 (Public review):

Summary

In this manuscript, the authors introduce Gcoupler, a Python-based computational pipeline designed to identify endogenous intracellular metabolites that function as allosteric modulators at the G protein-coupled receptor (GPCR) - Gα protein interface. Gcoupler is comprised of four modules:

I. Synthesizer - identifies protein cavities and generates synthetic ligands using LigBuilder3

II. Authenticator - classifies ligands into high-affinity binders (HABs) and low-affinity binders (LABs) based on AutoDock Vina binding energies

III. Generator - trains graph neural network (GNN) models (GCM, GCN, AFP, GAT) to predict binding affinity using synthetic ligands

IV. BioRanker - prioritizes ligands based on statistical and bioactivity data

The authors apply Gcoupler to study the Ste2p-Gpa1p interface in yeast, identifying sterols such as zymosterol (ZST) and lanosterol (LST) as modulators of GPCR signaling. Our review will focus on the computational aspects of the work. Overall, we found the Gcoupler approach interesting and potentially valuable, but we have several concerns with the methods and validation that need to be addressed prior to publication/dissemination.

(1) The exact algorithmic advancement of the Synthesizer beyond being some type of application wrapper around LigBuilder is unclear. Is the grow-link approach mentioned in the methods already a component of LigBuilder, or is it custom? If it is custom, what does it do? Is the API for custom optimization routines new with the Synthesizer, or is this a component of LigBuilder? Is the genetic algorithm novel or already an existing software implementation? Is the cavity detection tool a component of LigBuilder or novel in some way? Is the fragment library utilized in the Synthesizer the default fragment library in LigBuilder, or has it been customized? Are there rules that dictate how molecule growth can occur? The scientific contribution of the Synthesizer is unclear. If there has not been any new methodological development, then it may be more appropriate to just refer to this part of the algorithm as an application layer for LigBuilder.

(2) The use of AutoDock Vina binding energy scores to classify ligands into HABs and LABs is problematic. AutoDock Vina's energy function is primarily tuned for pose prediction and displays highly system-dependent affinity ranking capabilities. Moreover, the HAB/LAB thresholds of -7 kcal/mol or -8 kcal/mol lack justification. Were these arbitrarily selected cutoffs, or was benchmarking performed to identify appropriate cutoffs? It seems like these thresholds should be determined by calibrating the docking scores with experimental binding data (e.g., known binders with measured affinities) or through re-scoring molecules with a rigorous alchemical free energy approach.

(3) Neither the Results nor Methods sections provide information on how the GNNs were trained in this study. Details such as node features, edge attributes, standardization, pooling, activation functions, layers, dropout, etc., should all be described in detail. The training protocol should also be described, including loss functions, independent monitoring and early stopping criteria, learning rate adjustments, etc.

(4) GNN model training seems to occur on at most 500 molecules per training run? This is unclear from the manuscript. That is a very small number of training samples if true. Please clarify. How was upsampling performed? What were the HAB/LAB class distributions? In addition, it seems as though only synthetically generated molecules are used for training, and the task is to discriminate synthetic molecules based on their docking scores. Synthetic ligands generated by LigBuilder may occupy distinct chemical space, making classification trivial, particularly in the setting of a random split k-folds validation approach. In the absence of a leave-class-out validation, it is unclear if the model learns generalizable features or exploits clear chemical differences. Historically, it was inappropriate to evaluate ligand-based QSAR models on synthetic decoys such as the DUD-E sets - synthetic ligands can be much more easily distinguished by heavily parameterized ligand-based machine learning models than by physically constrained single-point docking score functions.

(5) Training QSAR models on docking scores to accelerate virtual screening is not in itself novel (see here for a nice recent example: https://www.nature.com/articles/s43588-025-00777-x), but can be highly useful to focus structure-based analysis on the most promising areas of ligand chemical space; however, we are perplexed by the motivation here. If only a few hundred or a few thousand molecules are being sampled, why not just use AutoDock Vina? The models are trained to try to discriminate molecules by AutoDock Vina score rather than experimental affinity, so it seems like we would ideally just run Vina? Perhaps we are misunderstanding the scale of the screening that was done here. Please clarify the manuscript methods to help justify the approach.

(6) The brevity of the MD simulations raises some concerns that the results may be over-interpreted. RMSD plots do not reliably compare the affinity behavior in this context because of the short timescales coupled with the dramatic topological differences between the ligands being compared; CoQ6 is long and highly flexible compared to ZST and LST. Convergence metrics, such as block averaging and time-dependent MM/GBSA energies, should be included over much longer timescales. For CoQ6, the authors may need to run multiple simulations of several microseconds, identify the longest-lived metastable states of CoQ6, and perform MM/GBSA energies for each state weighted by each state's probability.

Reviewer #2 (Public review):

Summary:

Mohanty et al. present a new deep learning method to identify intracellular allosteric modulators of GPCRs. This is an interesting field for e.g. the design of novel small molecule inhibitors of GPCR signalling. A key limitation, as mentioned by the authors, is the limited availability of data. The method presented, Gcoupler, aims to overcome these limitations, as shown by experimental validation of sterols in the inhibition of Ste2p, which has been shown to be relevant molecules in human and rat cardiac hypertrophy models.<br /> They have made their code available for download and installation, which can easily be followed to set up software on a local machine.

Strengths:

- Clear GitHub repository

- Extensive data on yeast systems

Weaknesses:

- No assay to directly determine the affinity of the compounds to the protein of interest.

In conclusion, the authors present an interesting new method to identify allosteric inhibitors of GPCRs, which can easily be employed by research labs. Whilst their efforts to characterize the compounds in yeast cells, in order to confirm their findings, it would be beneficial if the authors show their compounds are active in a simple binding assay.

Reviewer #3 (Public review):

Summary:

In this paper, the authors introduce the Gcoupler software, an open-source deep learning-based platform for structure-guided discovery of ligands targeting GPCR interfaces.<br /> Overall, this manuscript represents a field-advancing contribution at the intersection of AI-based ligand discovery and GPCR signaling regulation.

Strengths:

The paper presents a comprehensive and well-structured workflow combining cavity identification, de novo ligand generation, statistical validation, and graph neural network-based classification. Notably, the authors use Gcoupler to identify endogenous intracellular sterols as allosteric modulators of the GPCR-Gα interface in yeast, with experimental validations extending to mammalian systems. The ability to systematically explore intracellular metabolite modulation of GPCR signaling represents a novel and impactful contribution. This study significantly advances the field of GPCR biology and computational ligand discovery.

Reviewer #1 (Public review):

Summary:

In this study, Tittelmeier et al. explored the role of sphingolipid metabolism in maintaining endolysosomal membrane integrity and its downstream effects on tau aggregation and toxicity, using both worms and human cell models. The authors showed that knockdown of sphingolipid metabolism genes reduced endolysosomal membrane fluidity, as revealed by FRAP and C-Laurdan imaging, leading to increased vesicle rupture. Furthermore, tau aggregates accumulated in endolysosomes and exacerbated membrane rigidity and damage, promoting seeded tau aggregation, likely by enabling tau seed escape into the cytosol. Importantly, unsaturated fatty acid supplementation restored membrane fluidity, suppressed tau propagation, and alleviated neurotoxicity in C. elegans. These findings provide insight into how lipid dysregulation contributes to tau pathology and highlight membrane fluidity restoration as a potential therapeutic avenue for Alzheimer's disease.

Strengths:

The study addresses the connection between sphingolipid metabolism, endolysosomal membrane integrity, and tau pathology, which is a relevant topic in the context of Alzheimer's disease and related tauopathies.

The use of both C. elegans and human cell models provides cross-species perspectives that help frame the findings in a broader biological context.

The combination of FRAP and C-Laurdan dye imaging offers a biophysical approach to investigate changes in membrane properties, which is a technically interesting aspect of the study.

The observation that unsaturated fatty acid supplementation can modulate membrane fluidity and influence tau-related phenotypes adds an element of potential therapeutic interest.

The study presents multiple experimental approaches to address the proposed mechanism, and efforts were made to examine both membrane behavior and tau aggregation dynamics.

Weaknesses:

In Figure 3, the authors used C-Laurdan imaging to assess membrane fluidity and showed that knockdown of SPHK2, the human ortholog of sphk-1, led to increased membrane rigidity. However, the authors did not co-stain with a lysosomal marker, making it unclear whether the observed effect is specific to lysosomal membranes or reflects general membrane changes. Co-staining with LysoTracker or applying segmentation masks to isolate lysosomal signals would significantly improve interpretation.

Line 173 states that Lipofectamine 2000 increases membrane fluidity based on GP index changes, but this is incorrect. A higher GP index indicates increased membrane order (i.e., reduced fluidity), so the statement should be revised. Additionally, Lipofectamine 2000 can itself alter membrane rigidity, posing a risk of false-positive interpretations. To confirm the role of SPHK2 in this phenotype, the authors should use a CRISPR/Cas9 knockout model instead of relying solely on siRNA transfection, which may be confounded by the delivery reagent. Without lysosomal co-staining and SPHK2 KO validation, the authors cannot conclusively claim that SPHK2 loss affects endolysosomal membrane integrity.

The section titled "Fibrillar tau increases membrane rigidity and exacerbates endolysosomal damage" (lines 177-215) requires substantial revision. The narrative jumps abruptly between worms and cell models, making it hard to follow the logic. The use of the F3ΔK281::mCherry strain is introduced without explanation or context. It is unclear whether this strain is relevant to lysosomal membrane rupture, as no reference or justification is provided. The authors should clarify whether this reporter is intended to detect lysosomal membrane permeabilization (LMP). If so, it would be more appropriate to use established LMP reporters, such as lysosome-targeted fluorescent sensors, galectin-based reporters, or dextran leakage assays. Based on the current data in Figure 3G, it is difficult to draw firm conclusions regarding membrane rupture levels.

To support the conclusion that sphingolipid metabolism gene knockdown alters membrane properties, the study would benefit from direct lipidomic analysis. Measuring changes in sphingolipid profiles in both C. elegans and cell models would provide biochemical evidence for the proposed disruption of lipid homeostasis. Given the availability of lipidomics platforms, this type of analysis should be feasible in both worms and human cells and would significantly strengthen the mechanistic claims regarding membrane fluidity and integrity.

The conclusions of the study rely heavily on imaging-based assays, including FRAP, C-Laurdan, and fluorescence microscopy. While these approaches provide valuable spatial and qualitative insights, they are inherently indirect and subject to interpretive limitations. To strengthen the mechanistic claims, the authors should incorporate additional biochemical or quantitative approaches. For example, lipidomics would allow direct measurement of membrane lipid composition changes, and western blotting or quantitative proteomics could assess levels of membrane-associated proteins involved in endolysosomal function or stress responses. Including such data would significantly improve the robustness and reproducibility of the study's conclusions.

The human cell experiments were performed exclusively in HEK293T cells, which are not physiologically relevant for modeling Alzheimer's disease or lysosomal function in neurons. Given that the study aims to draw conclusions related to tau aggregation and lysosomal membrane integrity, the use of a more disease-relevant cellular model is essential. There are several established AD-relevant cell models, including iPSC-derived neurons, neuroblastoma lines expressing tau, or microglial models, which would better reflect the cellular context of tauopathies. Validation of key findings in at least one of these systems would substantially enhance the biological relevance and translational impact of the study.

The authors reported that PUFA supplementation rescues neurotoxic phenotypes by increasing membrane fluidity. However, the data supporting this claim rely entirely on confocal imaging, shown in both the main and supplemental figures. To substantiate the mechanistic link between PUFA treatment and improved lysosomal membrane properties, the authors should include functional assays demonstrating that PUFAs are indeed incorporated into lysosomal membranes. Additionally, lipidomics analysis would be valuable to identify which lipid species are altered upon supplementation and correlate these changes with the observed phenotypic rescue. Furthermore, the conclusion that PUFAs rescue "neurotoxic phenotypes" is not appropriate based on data derived solely from HEK293T cells, which are not neuronal. To make claims about tau-related neurotoxicity, the authors should validate their findings in a more relevant neuronal model, such as SH-SY5Y neuroblastoma cells expressing tau or iPSC-derived neurons. This would better reflect the cellular environment of Alzheimer's disease and provide stronger support for the proposed therapeutic potential of PUFA supplementation.

While the authors demonstrate that ALA supplementation mitigates neurotoxicity in C. elegans expressing aggregated tau (F3ΔK281::mCherry), the current data are not sufficient to conclude that ALA directly rescues tau aggregation toxicity via a lysosome-specific mechanism. It remains unclear how lipid composition is altered upon ALA treatment and whether these changes correlate with functional improvement of lysosomal pathways. The manuscript does not provide mechanistic insight into how ALA enhances lysosomal health or attenuates endolysosomal damage. Moreover, supplementation with PUFAs like ALA can activate a wide range of cellular processes beyond lysosomal function, including alterations in membrane fluidity, signaling cascades, and oxidative stress responses. The authors should clarify how they distinguish the lysosome-related effects from these alternative pathways. For example, did they observe specific lysosomal markers or structural improvements in lysosomes upon ALA treatment? Additional data or controls would be necessary to support a lysosome-specific protective mechanism and to exclude the involvement of other PUFA-responsive pathways in the observed phenotypes.

Reviewer #2 (Public review):

Tittelmeier et al. investigated the role of sphingolipid (SL) metabolism in the maintenance of endolysosomal vesicle integrity. They find that both impaired SL biosynthesis and degradation in C. elegans, decrease the fluidity of endolysosomal membranes and promote their rupture, while it has little effect on plasma membrane fluidity. Endolysosomal membrane fluidity is also negatively affected in human cells upon knockdown (KD) of a gene (SPHK2) involved in the SL degradation pathway. Aggregated forms of tau in both models (C. elegans and human cells) can also cause rigidification of the endolysosomal membrane, with SL homeostasis disruption having an additive effect, exacerbating endolysosomal rupture. Notably, KD of SPHK2 also increased the formation of tau foci, suggesting that compromised endolysosomal integrity may promote tau aggregation. These data provide a clearer understanding of how genetic manipulation of SL metabolism affects endolysosomal membranes and their rigidification in the context of tau aggregation. Supplementation of polyunsaturated fatty acids (PUFAs), which has a beneficial effect on Alzheimer's patients, improved membrane fluidity and reduced tau propagation in human cells and tau-associated neurotoxicity in C. elegans, suggesting a possible mechanism of action.

Overall, the conclusions of this paper are supported by the data, with a few aspects requiring further clarification and elaboration.

(1) A reference to Figure S2E-G, which shows that KD of SL biosynthesis genes do not affect the plasma membrane, is missing from the main text.

(2) In Figure 3C, lipofectamine alone shows that it increases membrane rigidity (increased GP values), not membrane fluidity.

(3) In Figure 3F, the EV cntl condition expressing F3:mCh tau should have increased LGALS3 foci compared to the mCh EV cntl according to Ref (20) and its Figure 2G (at least for Day 5 animals), which would be indicative of the tau spreading in hypodermal tissue. What C. elegans age was examined in Figure 3F? Can the authors provide evidence of the transmission of the F3:mCh tau from the touch receptor neurons to the hypodermis in the EV [similar to Figure 2C & D from Ref (20)] and compare it to the KDs? Otherwise, it seems that KD of SL genes impacts not only endolysosomal rupture but significantly affects tau accumulation/spreading as well (e.g., shown later in HEK cells, where SPHK2 KD increases the formation of tau-Venus foci).

(4) Sphingolipids are essential membrane components and signaling molecules. Does KD of SL genes in C. elegans and the subsequent endolysosomal rupture cause any major, intermediate, or minor defects/phenotypes (in non-aggregation prone models, w/t..)?

Reviewer #3 (Public review):

Summary:

The authors set off with an analysis of the lysosomal integrity upon knockdown of genes of the sphingolipid metabolic pathway that they identified in a previous (yet unpublished) work of an RNA screen using a new C. elegans Tau model. They then used cell culture and C. elegans experiments to study the link between lysosomal rupture and Tau propagation.

Strengths:

The authors use two complementary model systems and use probes to assess membrane rigidity that allow a quick assessment of the membrane dynamics and offer the opportunity to treat the cells with lipids, RNAi. Tau seeds, etc.

Weaknesses:

The main weakness is that this work builds on not-yet-peer-reviewed manuscript that established a new C. elegans Tau model and RNAi screen that aimed to identify genes involved in the propagation of Tau.

This reviewer misses essential information of the C. elegans Tau strain (not included in the method section): e.g., promoter used for the expression, information on the used Tau variant, expression pattern, and aggregation, etc.

Throughout the study, I missed data on:

(1) Effect of the knockdown on Tau expression, localisation (with lysosomal membrane?), aggregation, and proteotoxicity. The effect of the RNAi-mediated knockdown could also simply lead to a reduced expression of Tau that, in turn, leads to suppressed propagation.

(2) A quantification of RNAi knockdown is needed to judge the efficiency of the RNAi, in particular for the combinatorial RNAi experiments involving 2 and even 4 genes in parallel. Ideally, these analyses should be validated with mutants for these genes.

Further:

(3) Figure 4 H, I: Would Tau also aggregate in the absence of externally added Tau?

(4) How specific is the effect for Tau? It would help if the authors could assess other amyloid proteins.

(5) The connection between sphingolipids and AD is not new. See He et al, 2010, Neurobiol. Aging + numerous publications and also not between Tau seeding and lysosomal rupture: Rose et al., PNAS 2024 (that has been cited by the authors).

Reviewer #1 (Public review):

Summary:

It is now increasingly becoming clear that macromolecules and their complexes can form larger structures such as filaments or cages in the cells under certain conditions. These can be beneficial for the cells to promote and coordinate metabolic activity or result in protection against stress. Reactive oxygen species (ROS) can be damaging to macromolecules in cells that grow both aerobically and anaerobically, and they have evolved different mechanisms to cope with ROS. Aerobic organisms have a number of enzymes to combat ROS, while anaerobic organisms have evolved other means, and one such mechanism is described by Song et al in the article.<br /> In Pyrococcus furiosus, a hyperthermophilic anaerobic bacterium, Song et al describe the formation of Oxidative stress-induced tubular structures (OSITs). Using proteomics and electron cryomicroscopy (CryoEM), the authors find that the protein Rubrerythrin is upregulated upon exposure to oxygen, and the tetramer of this protein assembles to form these tubules that are varied in length with a consistent diameter of ~480 Å. They further observe that some of these tubules also have spherical viral-like particles. With enriched fraction of the OSITs from the cells and proteomics, it is shown that the predominant protein is encapsulin, which forms a caged structure and traps ferric iron. The combined structures of OSIT by rubreerythrin and the VLPs of encapsulin protect the cells from oxygen radicals by forming a complex.

Strengths:

The combination of proteomics and electron microscopy with the employment of both tomography of cellular sections and single particle cryoEM of enriched samples.

Weaknesses:

Some description of the methods, in particular the workflow of image processing, is not easy to follow and can be described with more clarity and be easier for non-experts to read/understand.

Reviewer #2 (Public review):

The manuscript entitled "Structure of an oxygen-induced tubular nanocompartment in Pyrococcus furiosus" by Wenfei Song et al. employs whole-cell mass spectrometry and cryo-EM (including tomography, helical reconstruction, and single-particle analysis) to investigate the structure and function of the oxidoreductase Rubrerythrin (Rbr) from Pyrococcus furiosus. The study reports that under oxidative stress, Rbr forms a tubular structure, in contrast to its behaviour under anaerobic conditions. Authors characterized oxidoreductase Rubrerythrin (Rbr) from Pyrococcus furiosus under anaerobic conditions and formed a tubular structure when induced with oxidative stress. This study is well-designed. However, I have several questions related to the experimental design and the results obtained from those experiments, which are listed below.

(1) The authors have mentioned that "Under aerobic conditions, Rbr levels are 3 to 13 times higher compared to anaerobic conditions (Figures 1a-d)." Also, they performed whole-cell mass spec to measure the overexpression of the Rbr enzyme under anaerobic conditions. Thus, from the above statement, I consider the authors' claim that P. furiosus cells were cultured under anaerobic conditions and then exposed to oxidative stress. While cell growth under anaerobic conditions appears perfectly fine, the authors conducted the rest of the experiment under aerobic conditions during mass spectrometry and cryo-EM sample preparation. As a baseline, the author first grew the cells in their preferred anaerobic environment and also imaged the same cells that were exposed to air (aerobic) after anaerobic growth. The cell growth in anaerobic conditions is perfectly fine. But how did authors make sure that during anaerobic conditions, the Rbr enzyme is not expressed or not formed? As a control experiment, authors should demonstrate that during mass spec and cryo-EM sample preparations, cells are not exposed to air or maintained in an anaerobic environment. From anaerobic conditions, whenever cells were selected for spec and cryo-EM, cells were exposed to O2, and definitely controlled cells were not in anaerobic conditions anymore.

The authors collected P. furiosus wild-type or Rbr knockout cells in an anaerobic hood, but after that, they centrifuged the cells and plunged them using a Vitrobot. Are the instrument, centrifuge, and Vitrobot kept in an anaerobic environment? Recently, a few studies (anaerobic plunge-freezing in cryo-electron microscopy, Cook et al. (2024), Hands-Portman and Bakker (2022) DOI: 10.1039/D2FD00060A ) have mentioned the anaerobic plunge freeze setup for protein sample or cell freezing. I guess the authors did not use that setup. In these circumstances, the cell is already exposed to O2 during centrifugation and Vitrobot freezing. How were the control experiments properly performed in anaerobic conditions? A similar argument is true for Lamella grid preparation, where the enzyme was already exposed to O2, and single-particle grid preparation, where the purified enzyme is already exposed to O2. How were the control experiments properly performed in anaerobic conditions?

(2) It is important to provide evidence that the overexpressed protein is actually in an anaerobic condition and is later induced with more O2. Also, authors should confirm biochemically that the overexpressed protein in their desired protein "oxidoreductase Rubrerythrin (Rbr)". No biochemical data were provided in this manuscript. During single-particle analysis, the authors had to purify the protein sample and confirm that these were their desired protein samples. No biochemical or biophysical experiments were performed to confirm that the overexpressed protein is the desired protein.

(3) Figure 3, the atomic model looks different in all four tetramers. However, I have fitted the atomic model into the cryo-EM map, which looks reasonable. However, it will be easier for the reader to evaluate the model if the authors show different orientations of the atomic model, as well as if the authors could show that the atomic model fits the cryo-EM map.

(4) How did the authors select initial particle sets like 24 lakhs when forming helices and not forming isolated particles?

(5) The authors proposed a model for electron transfer upon oxidative stress. However, the data is not convincing that VLP is surrounded by Rbr and forms a tube-like structure. Generally, VLP is a sphere-like structure, and Rbr can form a tube-like structure when it interacts with spherical VLP. Rbr will surround VLP, and it will form a Rbr-decorated sphere-like structure.

(6) It will also be important to comment on the diameter of Oxidative stress-induced tubules (OSITs) and 3D reconstruction and/or helical reconstruction of purified protein samples. The spherical cyan densities within the tube are not very clear. If VLP is surrounded by Rbr (Figure 4), extra Rbr densities will be observed on VLP in the tomogram (in Figure 1). However, in the tomogram, VLP is inside Oxidative stress-induced tubules (OSITs). Figure 1 is a contradicting Figure 4. The authors should explain it properly.

(7) The authors performed helical reconstruction. Where is the Layer line calculation in helical reconstruction, and how do authors identify helical parameters for reconstruction?

(8) The authors used an extremely confusing methodology, which was very difficult to follow. The authors performed tomography, helical reconstruction, and single-particle analysis. Why did the authors need 3 different image processing methods to resolve structures that are not clear to me? The authors should also show the proper fitting between the map and the model. In Supplemental Figure 6c, the overall fitting of the subdomain looks ok. However, many peptide chains and side chains are not fitted properly in the EM density map. It will be helpful to show proper side chain fitting. In Supplementary Fig. 6a, the authors binned the data (Bin 8 or Bin 2) but did not mention when they unbinned the data for data processing. Also, the authors implemented C2 symmetry during local refinement. Why do authors suddenly use C2 symmetry expansion?

Minor Comments:

(1) The authors should properly show a schematic diagram of the enzyme subdomains. It will help to understand interactions or tetrameric assembly.

(2) The introduction is poorly written. It will really be helpful for the reader if the authors provide a proper introduction.

(3) The atomic model did not fit into the cryo-EM, so it was hard to determine the overall fitting.

(4) 17.1A pixel size? It's surprising.

(5) It will be better to calculate local resolution and show the map's angular distribution. It is obvious that resolution at the peripheral region will be poorer than core region. Therefore, it will be better to calculate local resolution. Additionally, authors should show the map to model fitting.

Reviewer #3 (Public review):

Summary:

The manuscript authored by Song et al explores the oxidative stress response of Rubrerythrin in Pyrococcus furiosus and the formation of unique tubules that also encapsulate Encapsulin VLPs. This is an excellent study employing diverse methods to comprehensively study the formation of these assemblies under oxidative stress and lays the foundation of understanding oxidative stress through the formation of tubules among redox-sensing proteins like Rubrerythrin. The authors decipher the molecular structure of the tubules and also present a high-resolution reconstruction of the rubrerythin unit that forms the OSITs.

Strengths:

The study is done thoroughly by employing methods like cryoET, single particle cryoEM, mass spectrometry, and expression analyses of knockout strains to delve into an important mechanism to counter oxidative stress. The authors perform comprehensive analyses, and this study represents a vital contribution to understanding how anaerobic organisms can respond to oxidative stress.

Weaknesses:

Not all encapsulin particles seem to be inside the OSITs. Do the authors have any insights into how the tubules sequester these viral particles? Do the VLPs have a role in nucleating the OSIT assembly, and are there interactions between VLP and OSIT surfaces? These could be points that can be discussed in greater detail by the authors.

Can the authors get a subtomogram averaging done for the encapsulin VLPs? A higher resolution reconstruction may provide potential interaction details with the OSITs, if there are any.

The role of the dense granules observed in the rubrerythrin deletion strain is not very well discussed. Is there a way these granules counter oxidative stress? The EDX scanning seems to show a Phosphate increase similar to Ca and Mg. Are these aggregates therefore likely to be calcium and Mg phosphate aggregates? This section of the paper seems incompletely analysed.

The authors should provide density and coordination distances around the diiron ions and provide a comparison with available crystal structures and highlight differences, if any, in Figure 3. Local resolution for the high-res map may be provided for Supplementary Figure 6.

Overall, this is a well-performed study with clear conclusions. The discussion points need to be improved further.

Reviewer #1 (Public review):

Summary:

This study shows a novel role for SCoR2 in regulating metabolic pathways in the heart to prevent injury following ischemia/reperfusion. It combines a new multi-omics method to determine SCoR2 mediated metabolic pathways in the heart. This paper would be of interest to cardiovascular researchers working on cardioprotective strategies following ischemic injury in the heart.

Strengths:

(1) Use of SCoR2KO mice subjected to I/R injury.

(2) Identification of multiple metabolic pathways in the heart by a novel multi-omics approach.

Weaknesses:

(1) Use of a global SCoR2KO mice is a limitation since the effects in the heart can be a combination of global loss of SCoR2.

(2) Lack of a cell type specific effect.

Reviewer #2 (Public review):

Summary:

This manuscript addresses the gap in knowledge related to the cardiac function of the S-denitrosylase SNO-CoA Reductase 2 (SCoR2; product of the Akr1a1 gene). Genetic variants in SCoR2 have been linked to cardiovascular disease, yet their exact role in the heart remains unclear. This paper demonstrates that mice deficient in SCoR2 show significant protection in a myocardial infarction (MI) model. SCoR2 also affected ketolytic energy production, antioxidant levels, and polyol balance through the S-nitrosylation of crucial metabolic regulators.

Strengths:

(1) Addresses a well-defined gap in knowledge related to the cardiac role of SNO-CoA Reductase 2. Besides the in-depth case for this specific player, the manuscript sheds more light on the links between S-nitrosylation and metabolic reprogramming in the heart.

(2) Rigorous proof of requirement through the combination of gene knockout and in vivo myocardial ischemia/reperfusion.

(3) Identification of precise Cys residue for SNO-modification of BDH1 as SCoR2 target in cardiac ketolysis

Weaknesses:

(1) The experiments with BDH1 stability were performed in mutant 293 cells. Was there a difference in BDH1 stability in myocardial tissue or primary cardiomyocytes from SCoR2-null vs -WT mice? The same question extends to PKM2.

(2) In the absence of tracing experiments, the cross-sectional changes in ketolysis, glycolysis, or polyol intermediates presented in Figures 4 and 5 are suggestive at best. This needs to be stressed while describing and interpreting these results.

(3) The findings from human samples with ischemic and non-ischemic cardiomyopathy do not seem immediately or linearly in line with each other and with the model proposed from the KO mice. While the correlation holds up in the non-ischemic cardiomyopathy (increased SNO-BDH1, SNO-PKM2 with decreased SCoR2 expression), how do the authors explain the decreased SNO-BDH1 with preserved SCoR2 expression in ischemic cardiomyopathy? This seems counterintuitive as activation of ketolysis is a quite established myocardial response to ischemic stress. It may help the overall message clarity to focus the human data part on only NICM patients.

(4) This issue is partially linked to point #(3). Currently, important evidence that is lacking is the demonstration of sufficiency for SCoR2 in S-nytrosylation of targets and cardiac remodeling. Does SCoR2 overexpression in the heart or isolated cardiomyocytes reduce S-nitrosylation of BDH1 and other targets, thus affecting heart function at baseline or under stress?

Reviewer #3 (Public review):

Summary:

This manuscript demonstrates that mice lacking the denitrosylase enzyme SCoR2/AKR1A1 demonstrate robust cardioprotection resulting from reprogramming of multiple metabolic pathways, revealing widespread, coordinated metabolic regulation by SCoR2.

Strengths:

(1) The extensive experimental evidence.

(2) The use of the knockout model.

Weaknesses:

(1) Lack of direct evidence for underlying mechanism(s).

(2) The mouse model used is not tissue-specific.

Reviewer #1 (Public review):

The authors present exciting new experimental data on the antigenic recognition of 78 H3N2 strains (from the beginning of the 2023 Northern Hemisphere season) against a set of 150 serum samples. The authors compare protection profiles of individual sera and find that the antigenic effect of amino acid substitutions at specific sites depends on the immune class of the sera, differentiating between children and adults. Person-to-person heterogeneity in the measured titers is strong, specifically in the group of children's sera. The authors find that the fraction of sera with low titers correlates with the inferred growth rate using maximum likelihood regression (MLR), a correlation that does not hold for pooled sera. The authors then measure the protection profile of the sera against historical vaccine strains and find that it can be explained by birth cohort for children. Finally, the authors present data comparing pre- and post- vaccination protection profiles for 39 (USA) and 8 (Australia) adults. The data shows a cohort-specific vaccination effect as measured by the average titer increase, and also a virus-specific vaccination effect for the historical vaccine strains. The generated data is shared by the authors and they also note that these methods can be applied to inform the bi-annual vaccine composition meetings, which could be highly valuable.

The following points could be addressed in a revision:

(1) The authors conclude that much of the person-to-person and strain-to-strain variation seems idiosyncratic to individual sera rather than age groups. This point is not yet fully convincing. While the mean titer of an individual may be idiosyncratic to the individual sera, the strain-to-strain variation still reveals some patterns that are consistent across individuals (the authors note the effects of substitutions at sites 145 and 275/276). A more detailed analysis, removing the individual-specific mean titer, could still show shared patterns in groups of individuals that are not necessarily defined by the birth cohort.

(2) The authors show that the fraction of sera with a titer below 138 correlates strongly with the inferred growth rate using MLR. However, the authors also note that there exists a strong correlation between the MLR growth rate and the number of HA1 mutations. This analysis does not yet show that the titers provide substantially more information about the evolutionary success. The actual relation between the measured titers and fitness is certainly more subtle than suggested by the correlation plot in Figure 5. For example, the clades A/Massachusetts and A/Sydney both have a positive fitness at the beginning of 2023, but A/Massachusetts has substantially higher relative fitness than A/Sydney. The growth inference in Figure 5b does not appear to map that difference, and the antigenic data would give the opposite ranking. Similarly, the clades A/Massachusetts and A/Ontario have both positive relative fitness, as correctly identified by the antigenic ranking, but at quite different times (i.e., in different contexts of competing clades). Other clades, like A/St. Petersburg are assigned high growth and high escape but remain at low frequency throughout. Some mention of these effects not mapped by the analysis may be appropriate.

(3) For the protection profile against the vaccine strains, the authors find for the adult cohort that the highest titer is always against the oldest vaccine strain tested, which is A/Texas/50/2012. However, the adult sera do not show an increase in titer towards older strains, but only a peak at A/Texas. Therefore, it could be that this is a virus-specific effect, rather than a property of the protection profile. Could the authors test with one older vaccine virus (A/Perth/16/2009?) whether this really can be a general property?

Reviewer #2 (Public review):

This is an excellent paper. The ability to measure the immune response to multiple viruses in parallel is a major advancement for the field, which will be relevant across pathogens (assuming the assay can be appropriately adapted). I only have a few comments, focused on maximising the information provided by the sera.

Firstly, one of the major findings is that there is wide heterogeneity in responses across individuals. However, we could expect that individuals' responses should be at least correlated across the viruses considered, especially when individuals are of a similar age. It would be interesting to quantify the correlation in responses as a function of the difference in ages between pairs of individuals. I am also left wondering what the potential drivers of the differences in responses are, with age being presumably key. It would be interesting to explore individual factors associated with responses to specific viruses (beyond simply comparing adults versus children).

Relatedly, is the phylogenetic distance between pairs of viruses associated with similarity in responses?

Figure 5C is also a really interesting result. To be able to predict growth rates based on titers in the sera is fascinating. As touched upon in the discussion, I suspect it is really dependent on the representativeness of the sera of the population (so, e.g., if only elderly individuals provided sera, it would be a different result than if only children provided samples). It may be interesting to compare different hypotheses - so e.g., see if a population-weighted titer is even better correlated with fitness - so the contribution from each individual's titer is linked to a number of individuals of that age in the population. Alternatively, maybe only the titers in younger individuals are most relevant to fitness, etc.

In Figure 6, the authors lump together individuals within 10-year age categories - however, this is potentially throwing away the nuances of what is happening at individual ages, especially for the children, where the measured viruses cross different groups. I realise the numbers are small and the viruses only come from a small numbers of years, however, it may be preferable to order all the individuals by age (y-axis) and the viral responses in ascending order (x-axis) and plot the response as a heatmap. As currently plotted, it is difficult to compare across panels

Reviewer #3 (Public review):

The authors use high-throughput neutralisation data to explore how different summary statistics for population immune responses relate to strain success, as measured by growth rate during the 2023 season. The question of how serological measurements relate to epidemic growth is an important one, and I thought the authors present a thoughtful analysis tackling this question, with some clear figures. In particular, they found that stratifying the population based on the magnitude of their antibody titres correlates more with strain growth than using measurements derived from pooled serum data. However, there are some areas where I thought the work could be more strongly motivated and linked together. In particular, how the vaccine responses in US and Australia in Figures 6-7 relate to the earlier analysis around growth rates, and what we would expect the relationship between growth rate and population immunity to be based on epidemic theory.

Reviewer #1 (Public review):

Summary:

In this manuscript, Uphoff et al. propose a structural and mechanistic model in which the multidomain ECM protein SVEP1 enables Angiopoietin (ANG) binding to the orphan receptor TIE1, thereby promoting downstream receptor phosphorylation and signaling. Using AlphaFold-based modeling, the authors predict that the CCP20 domain of SVEP1 binds to TIE1, creating a composite surface that facilitates Angiopoietin association and TIE1 activation. The resulting ternary model (SVEP1-TIE1-ANG) offers a structural rationale for how SVEP1 converts TIE1 into a functional, ligand-responsive receptor. Additional models and biological assays suggest roles for other domains of SVEP1, such as CCP5-EGF-L7, although these interactions are predicted with low confidence. The authors interpret these findings as the first structural framework for how SVEP1 enables ANG-TIE1 signaling.

Strengths:

(1) The central hypothesis - that SVEP1 enables ANG binding to the orphan receptor TIE1 - is biologically compelling and addresses an important question in vascular biology.

(2) The AlphaFold-predicted ternary complex (SVEP1-TIE1-ANG) is plausible, high-confidence, and structurally consistent with prior functional data (e.g., poly-Ala scanning from Sato-Nishiuchi et al.).

(3) The authors' model offers a potential explanation for the previously observed role of SVEP1 in enhancing ANG signaling through TIE1, and may represent the first structural insight into TIE1's transition from orphan to ligand-activated receptor.

(4) The potential clinical implication - that a combinatorial ligand (ANG+SVEP1) can activate TIE1- could have translational relevance for vascular leak and inflammatory disease.

Weaknesses:

(1) Lack of structural validation and mechanistic follow-up:
Despite the promising AlphaFold model, there are no figures of the predicted interface, no residue-level interactions shown, no ipTM values reported, and no experimental follow-up to test the interface. PAE plots are incorrectly used as confidence justifications, which is not appropriate for complex predictions.

(2) Biophysical validation is missing:
No surface plasmon resonance (SPR), ITC, or biochemical assays are included to confirm ternary complex formation or quantify binding kinetics. Given the manuscript's structural focus, this is a major gap. For instance, an SPR experiment where ANG is immobilized, and TIE1 binding is measured {plus minus} SVEP1, would directly test the model. And allow direct comparison to ANG-TIE2.

(3) Missed opportunity for mutagenesis-driven validation:
 The manuscript does not include any interface-targeted mutations, despite clear opportunities. For example, mutating T2595 in SVEP1 (to R) or mutating the TIE1-specific residues (residues PL 202-203 to LF) could strongly test the model and potentially reveal dominant-negative behaviors. E.g. A T2595 mutant should block ANG binding but not TIE1 binding.

(4) Overinterpretation of weak models:
The additional AlphaFold model involving the CCP5-EGFL7 domains binding TIE1 has extremely low confidence (ipTM < 0.15) when reexamined by this reader and should not be emphasized. There is no biophysical evidence or binding data (SPR) to support this interaction, and its inclusion detracts from the much stronger CCP20 model.

(5) Language around modeling is overstated and potentially misleading:
Terms like "unequivocal," "high-affinity," or "affirms strong binding" in reference to AlphaFold predictions are inappropriate. These are hypotheses -not confirmations - and must be tested at the biochemical level. This should be clarified throughout the manuscript to ensure non-experts do not misinterpret modeling confidence as binding affinity.

(6) Negative stain EM data is not informative due to low resolution and lack of defined interfaces; unless replaced by higher-resolution Cryo-EM, this should be omitted. Better would be co-gel filtration, AUC, or SEC-MALLs with ANG-SVEP1-TIE1.

(7) Disjointed narrative:
The manuscript presents a compelling mechanism involving CCP20-driven ANG binding to TIE1, but then becomes fragmented by introducing the low-confidence CCP5-EGFL7 model and speculative higher-order polymerization models that are not experimentally supported.

Reviewer #2 (Public review):

Uphoff and colleagues present the results of a study focused on characterizing the binding of SVEP1 to TIE1 along with Angiopoietin-2. Starting with computational prediction of SVEP1 binding to TIE1, the authors identify the region of SVEP1 that serves as a high-affinity ligand for TIE1. Advanced studies identify a weak secondary binding site within SVEP1 that appears to be sufficient but not necessary for its interaction with TIE1 based on in vivo rescue experiments. The most novel contribution of the manuscript seems to be the identification of angiopoietin-1 and -2 as co-factors that seem to enhance the binding of SVEP1 with TIE1 and impact downstream AKT signaling. They propose a complex in which SVEP1 binds to TIE1 and ANG2.

Although the first set of results is essentially confirmatory, the identification of ANG-2 as a "co-factor" enhancing the binding of SVEP1 to TIE1 and associated downstream signaling (i.e., Figures 3 and 4) is novel and is of interest. However, the manuscript and its conclusions would greatly benefit from some clarifying details and additional experiments to ensure rigor and support specific claims.

Reviewer #1 (Public review):

Engineering of AAV replication proteins may provide new insights into Parvoviral replication and potentially enable improved recombinant AAV vector yield when incorporated into the manufacturing process. Silberg and colleagues report an AAV Rep library, that is an interesting and powerful approach, however, the screening design and subsequent experiments lack rigor and ultimately the results are premature. Overall, the manuscript does not accurately describe state-of-the-art in the field and has significant shortcomings with experimental design/data analysis. Key concerns are noted below:

The high enrichment of P19 variants in the library was likely an artifact of the fact they only transfected 20 ng of RepCap into their 5-plate preps. When such little Rep is provided, any boost in Rep expression levels will have a major on yield. When more RepCap is provided, 10 ug in their later evaluation, small changes in Rep expression are unlikely to have major impacts on yield. A more effective strategy would have been to transfect a normal amount of DNA and then utilize serial passaging through infectious cycling to account for cross packaging.

Introduction:<br /> - There are 7 FDA approved AAV gene therapies.<br /> - The description of "shuffling" when citing Mietsczh et al is misapplied. The cited paper discusses rationally designed hybrids.<br /> - The graphic represents a hybrid capsid, but the focus is rep. As such, this should be depicted differently.

Figures 1 and 2 are validation of previously published findings and general optimization of the experimental conditions. These do not provide the reader any new insight or information.

In Figure 3: The experimental approach is limited. It is unclear how the subpooling of different conditions was performed. As mentioned above, their library transfection strategy will significantly bias the results. The enriched variants have not been evaluated - specifically, the enriched non-synonymous mutations have not been shown to yield higher titers when tested individually.<br /> In Figure 4: The claim is made that "several synonymous mutations within the p19 promoter region increase Rep DNA packaging activity." However, Figure 4c does not show statistically significant differences in support of this claim. Additional supporting data is needed. Further, Authors state that the synonymous mutations are near the P19 promoter. However, looking at the sequence shown in figure 4, their annotation of the P19 promoter is not correct and the mutations are actually within the P19 promoter. Relatedly, the authors note that mutations enriched in the p19 region include additional tetranucleotide repeats. No synthetic variants with additional GCTCs have been generated to test this hypothesis. Further, these results would benefit from a Western blot and transcript analysis to confirm Rep52/40, expression levels of constructs.

Reviewer #2 (Public review):

In the present study the authors have investigated the effects of mutations on Rep protein ability to package DNA within the gene therapy vector, AAV. A detailed investigation of Rep mutants selected from a library has been probed for their ability to produce active virions. While the concept is interesting the outcome effects are very limited.

The major issue is the lack of immediate applicability and relevance in the vector production pipeline for AAV. The authors have found that with the synthetic GFP transgene cargo, mutations of the p19 promoter did not lead to enhanced AAV vector packing. Thus the data is quite preliminary and a complete characterization may be necessary to further enhance the translational potential of the approach.

Reviewer #3 (Public review):

While the AAV capsid has long been the target of protein engineering, its Rep proteins have been comparatively less studied. Since Rep plays a variety of roles for genome replication and virion packaging, gaining a deeper mechanistic understanding of their function and/or engineering versions that enable higher packaging productivity would be of interest to the field. This study generates a library of non-synonymous mutations in AAV2 rep (intended to cover all 19 aa changes at all positions, and coming close), packaged an AAV with AAV9 capsid, and sequenced the results to assess which amino acid changes resulted in an enrichment/depletion of genomes containing that variant rep. They found that proline substitutions are disruptive, well known from protein mutagenesis studies. The most significant enrichment sfound, however, were a set of synonymous mutations (unintended members of the library, as the library was designed to contain non-synonymous mutations) that lie within the p19 promoter. However, attempts to package recombinant vector using these individual rep variants in the AAV helper construct did not increase viral titer.

A previous study conducted analogous mutagenesis on Rep: Jain et al., "Comprehensive mutagenesis maps the effect of all single-codon mutations in the AAV2 rep gene on AAV production" eLife 2024 (cited here as reference 19). It is not clear that this current study is a significant advance relative to the prior, quite comprehensive study. Both generated a library of non-synonymous mutations and observed fitness effects on Rep. Because this study sequenced the full rep, rather than barcodes associated with each rep variant, it found the enrichment in the synonymous mutations. However, these should ideally advance a basic understanding of Rep biology and/or result in better AAV production, but they did neither. It is speculated in the Discussion that the mutations generated additional GCTC motifs in p19, elements that mediate protein-DNA interactions. However, the role of GCTC motifs is speculative, and no transcriptional analysis is conducted. Furthermore, as discussed above, the mutations do not result in higher viral titers. Perhaps there's a transcriptional effect at the much lower copy numbers of vector genome present during library selection vs. rAAV packaging. They also found stop codons in Rep domains thought to be required for viral packaging, but functional studies confirming the screening findings are not conducted. As a result, the biological or technical relevance of the findings are extremely unclear, and thus the impact is relatively low.

The description of herring DNA co-transfection and cross-packaging (which is a well-known pitfall) are somewhat technical and arguably don't merit too much main manuscript attention.

Reviewer #1 (Public review):

Summary:

The authors test whether the archerfish can modulate the fast response to a falling target. By manipulating the trajectory of the target, they claim that the fish can modulate the fast response. While it is clear from the result that the fish can modulate the fast response, the experimental support for argument that the fish can do it for a reflex like behavior is inadequate.

Strengths:

Overall, the question that the authors raised in the manuscript is interesting.

Weaknesses:

Major comments:

(1) The argument that the fish can modulate reflex-like behavior relies on the claim that the archerfish makes the decision in 40 ms. There is little support for the 40 ms reaction time. The reaction time for the same behavior in Schlegel 2008, is 60-70 ms and in Tsvilling 2012 about 75 ms, if we take the half height of the maximum as estimated reaction time in both cases. If we take the peak (or average) of the distribution as an estimation of reaction time, the reaction time is even longer. This number is critical for the analysis the authors perform since if the reaction time is longer, maybe this is not a reflex as claimed. In addition, mentioning the 40 ms in the abstract is overselling the result. The title is also not supported by the results.

(2) A critical technical issue of the stimulus delivery is not clear. The frame rate is 120 FPS and the target horizontal speed can be up to 1.775 m/s. This produces target jumping on the screen 15 mm each frame. This is not a continuous motion. Thus, the similarity between the natural system where the target experience ballistic trajectory and the experiment here is not clear. Ideally, another type of stimulus delivery system is needed for a project of this kind that requires fast moving targets (e.g. Reiser, J. Neurosci.Meth. 2008). In addition, the screen is rectangular and not circular, so in some directions the target vanishes earlier than others. It must produce a bias in the fish response but there is no analysis of this type.

(3) The results here rely on the ability to measure the error of response in the case of virtual experiment. It is not clear how this is done since the virtual target does not fall. How do authors validate that the fish indeed perceives the virtual target as falling target? Since the deflection is at a later stage of the virtual trajectory, it is not clear what is the actual physics that governs the world of the experiment. Overall, the experimental setup is not well designed.

Comments on revisions:

The authors handled the comments, and the manuscript has improved accordingly. While some issues could still benefit from further clarification and depth, the current version meets the necessary standards.

Reviewer #2 (Public review):

Summary:

This manuscript studies the prey capture by archer fish, which observe the initial values of motion of aerial prey they made fall by spitting on them, and then rapidly turn to reach the ballistic landing point on the water surface. The question raised by the article is whether this incredibly fast decision-making process is hardwired and thus unmodifiable or can be adjusted by experience to follow a new rule, namely that the landing point is deflected from a certain amount from the expected ballistic landing point. The results show that the fish learn the new rule and use it afterwards in a variety of novel situations that include height, side and speed of the prey, and which preserve the speed of the fish's decision. Moreover, a remarkable finding presented in this work is the fact that fish that have learned to use the new rule can relearn to use the ballistic landing point for an object based on its shape (a triangle) while keeping simultaneously the 'deflected rule' for an object differing in shape (a disc); in other words, fish can master simultaneously two decision-making rules based on the different shape of objects.

Strengths:

The manuscript relies on a sophisticated and clever experimental design that allows changing the apparent landing point of a virtual prey using a virtual reality system. Several robust controls are provided to demonstrate the reliability and usefulness of the experimental setup.

Overall, I like very much the idea conveyed by the authors that even stimuli triggering apparently hardwired responses can be relearned in order to be associated to a different response, thus showing the impressive flexibility of circuits that are sometimes considered as mediating pure reflexive responses. This is the case - as an additional example - of the main component of the Nasanov pheromone of bees (geraniol), which triggers immediate reflexive attraction and appetitive responses, and which can, nevertheless, be learned by bees in association with an electric shock so that bees end up exhibiting avoidance and the aversive response of sting extension to this odorant(1), which is a fully unnatural situation, and which shows that associative aversive learning is strong enough to override preprogrammed responding, thus reflecting an impressive behavioral flexibility.

Weaknesses:

As a general remark, there is some information that I missed and that are mandatory in the analysis of behavioral changes: one is the variability in the performances displayed. The authors mentioned that the results reported come from 6 fish (which is a low sample size). How were the individual performances in terms of consistency? Were all fish equally good in adjusting/learning the new rule? How did errors vary according to individual identity? It seems to me that this kind of information should be available as the authors reported that individual fish could be recognized and tracked (see lines 620-635) and is essential for appreciating the flexibility of the system under study.

The other information that I could not find explained in a proper way is referred to the speed of the learning process. Admittedly, fish learn in an impressive way the new rule and even two rules simultaneously; yet, how long did they need to achieve this? In the article, Fig 2 mention that at least 6 training stages (each defined as a block of 60 evaluated turn decisions, which actually shows that the standard term 'Training Block' would be more appropriate) were required for the fish to learn the 'deflected rule'. While this means 360 trials (turning starts), I was left with the question of how long did this process last? How many hours, days, weeks were needed for the fish to learn? And as mentioned above, were al fish equally fast in learning? I would appreciate explaining this very important point because learning dynamics is relevant to understanding the flexibility of the system.

Comments After Revision:

There was consensus among reviewers that the authors addressed the initial critiques adequately and that the manuscript improved accordingly. The revision clarified several methodological aspects, and the addition of the new Fig. 2 was particularly helpful. It elucidates the experimental approach used in the study and offers essential context for understanding points that may have been unclear in the previous version.

Reviewer #1 (Public review):

In this study the authors aim to understand why decision formation during behavioural tasks is distributed across multiple brain areas. They hypothesize that multiple areas are used in order to implement an information bottleneck (IB). Using neural activity recorded from monkey DLPFC and PMd performing a 2-AFC task, they show that DLPFC represents various task variables (decision, color, target configuration), while downstream PMd primarily represents decision information. Since decision information is the only information needed to make a decision, the authors suggest that PMd has a minimal sufficient representation (as expected from an IB). They then train 3-area RNNs on the same task, and show that activity in the first and third areas resemble the neural representations of DLPFC and PMd, respectively. In order to propose a mechanism, they analyse the RNN and find that area 3 ends up with primarily decision information because feedforward connections between areas primarily propagate decision information.

Overall, the paper reads well and the data analysis and RNN modeling are well done and mostly correct. I agree with the authors that PMd has less information than DLPFC, meaning that some of the target and color information is attenuated. I also agree that this also happens in their multi-area RNN.

However, I find the use of the IB principle here muddles the water rather than clarifying anything. The key problem is that the authors evoke the information bottleneck in a mostly intuitive sense, but they do not actually use it (say, in their modelling). Rather, the IB is simply used to motivate why information will be or should be lost. Since the IB is a generic compressor, however, it does not make any statements about how a particular compression should be distributed or computed across brain areas.

If I ignore the reference to the information bottleneck, I still see a more mechanistic study that proposes a neural mechanism of how decisions are formed, in the tradition of RNN-modelling of neural activity as in Mante et al 2013. Seen through this more limited sense, the present study succeeds at pointing out a good model-data match.

Major points

(1) The IB is a formal, information-theoretic method to identify relevant information. However, in the paper, reference to the information bottleneck method (IB) is only used to motivate why (task-irrelevant) information should be lost in higher areas. The IB principle itself is actually never used. The RNNs are fitted using standard techniques, without reference to the IB. Without a formal link, I think the authors should describe their findings using words (e.g., task-irrelevant information is lost), rather than stating this as evidence for an information-theoretic principle.

(2) The advantage of employing a formal theory is that all assumptions have to be clarified. Since the authors only evoke the IB, but never employ it, they refrain from clarifying some of their assumptions. That is what creates unnecessary confusion.

For instance, the authors cite the following predictions of the IB principle: "(1) There exists a downstream area of cortex that has a minimal and sufficient representation to perform a task ... (2) there exists an upstream area of cortex that has more task information than the minimal sufficient area" - However, since the information bottleneck method is a generic compressor, it does not make any predictions about areas (or neurons). For a given sensory input p(x), a given task output p(y|x), and a given information loss, the IB generates exactly one optimal representation. In other words, the predictions made by the authors relie on other assumptions (e.g. feedforward processing, hierarchy, etc.) and these are not clearly stated.

(3) A corrollary to this problem is that the authors do not formally define task-irrelevant information. It seems the authors simply use the choice or decision as the thing that needs to be computed, and identify all other information as task-irrelevant. That's at least what I glean from the RNN model. However, I find that highly confusing because it suggests the conclusion that color information or target information are task-irrelevant. Surely, that cannot be true, since the decision is based on these quantities!

(4) If we define the output as the only task-relevant information, then any representation that is a pure motor representation would qualify as a minimal sufficient representation to carry out the correct actions. However, it is well-known that sensory information is lost in motor areas. It is not clear to me what exactly we gain by calling motor representations "minimal sufficient representations."

In summary, I think the authors should refrain from evoking the IB - which is a formal, mathematical principle - unless they actually use it formally as well.

Reviewer #2 (Public review):

This study advances our understanding of information encoding in the DLPFC and PMD brain regions. The conclusions are supported with convincing and robust analyses conducted on monkey datasets and trained RNN models. However, there are some concerns regarding the interpretation of findings related to the information bottleneck theory and the mapping of brain areas in the RNN simulations.

The authors' justification regarding mapping between model areas and anatomical areas remains insufficient, in my opinion. However, I recognize that my initial critique may not have been fully clear. The issue I see is this: whichever area is mapped to the first RNN module will trivially exhibit stimulus information, and downstream regions will naturally show a gradual loss of that information if one simply reads out their responses.

Thus, the observed stimulus loss in later modules could be an inevitable consequence of the model's structure, rather than a meaningful analog to the PFC-PMd transition. This point requires more careful justification or a reevaluation of the proposed mapping.

Reviewer #1 (Public review):

Summary:

In this manuscript, Pakula et al. explore the impact of reactive oxygen species (ROS) on neonatal cerebellar regeneration, providing evidence that ROS activates regeneration through Nestin-expressing progenitors (NEPs). Using scRNA-seq analysis of FACS-isolated NEPs, the authors characterize injury-induced changes, including an enrichment in ROS metabolic processes within the cerebellar microenvironment. Biochemical analyses confirm a rapid increase in ROS levels following irradiation and forced catalase expression, which reduces ROS levels, and impairs external granule layer (EGL) replenishment post-injury.

Strengths:

Overall, the study robustly supports its main conclusion and provides valuable insights into ROS as a regenerative signal in the neonatal cerebellum.

Comments on revisions:

The authors have addressed most of the previous comments. However, they should clarify the following response:

*"For reasons we have not explored, the phenotype is most prominent in these lobules, that is why they were originally chosen. We edited the following sentence (lines 578-579):

First, we analyzed the replenishment of the EGL by BgL-NEPs in vermis lobules 3-5, since our previous work showed that these lobules have a prominent defect."*

It has been reported that the anterior part of the cerebellum may have a lower regenerative capacity compared to the posterior lobe. To avoid potential ambiguity, the authors should clarify that "the phenotype" and "prominent defect" refer to more severe EGL depletion at an earlier stage after IR rather than a poorer regenerative outcome. Additionally, they should provide a reference to support their statement or indicate if it is based on unpublished observations.

Reviewer #2 (Public review):

Summary:

The authors have previously shown that the mouse neonatal cerebellum can regenerate damage to granule cell progenitors in the external granular layer, through reprogramming of gliogenic nestin-expressing progenitors (NEPs). The mechanisms of this reprogramming remain largely unknown. Here the authors used scRNAseq and ATACseq of purified neonatal NEPs from P1-P5 and showed that ROS signatures were transiently upregulated in gliogenic NEPs ve neurogenic NEPs 24 hours post injury (P2). To assess the role of ROS, mice transgenic for global catalase activity were assessed to reduce ROS. Inhibition of ROS significantly decreased gliogenic NEP reprogramming and diminished cerebellar growth post-injury. Further, inhibition of microglia across this same time period prevented one of the first steps of repair - the migration of NEPs into the external granule layer. This work is the first demonstration that the tissue microenvironment of the damaged neonatal cerebellum is a major regulator of neonatal cerebellar regeneration. Increased ROS is seen in other CNS damage models, including adults, thus there may be some shared mechanisms across age and regions, although interestingly neonatal cerebellar astrocytes do not upregulate GFAP as seen in adult CNS damage models. Another intriguing finding is that global inhibition of ROS did not alter normal cerebellar development.

Strengths:

This paper presents a beautiful example of using single cell data to generate biologically relevant, testable hypotheses of mechanisms driving important biological processes. The scRNAseq and ATACseq analyses are rigorously conducted and conclusive. Data is very clearly presented and easily interpreted supporting the hypothesis next tested by reduce ROS in irradiated brains.

Analysis of whole tissue and FAC sorted NEPS in transgenic mice where human catalase was globally expressed in mitochondria were rigorously controlled and conclusively show that ROS upregulation was indeed decreased post injury and very clearly the regenerative response was inhibited. The authors are to be commended on the very careful analyses which are very well presented and again, easy to follow with all appropriate data shown to support their conclusions.

Weaknesses:

The authors also present data to show that microglia are required for an early step of mobilizing gliogenic NEPs into the damaged EGL. While the data that PLX5622 administration from P0-P5 or even P0-P8 clearly shows that there is an immediate reduction of NEPs mobilized to the damaged EGL, there is no subsequent reduction of cerebellar growth such that by P30, the treated and untreated irradiated cerebella are equivalent in size. There is speculation in the discussion about why this might be the case. Additional experiments and tools are required to assess mechanisms. Regardless, the data still implicate microglia in the neonatal regenerative response, and this finding remains an important advance.

Reviewer #1 (Public review):

Summary:

Jin et al. investigated how the bacterial DNA damage (SOS) response and its regulator protein RecA affects the development of drug resistance under short-term exposure to beta-lactam antibiotics. Canonically, the SOS response is triggered by DNA damage, which results in the induction of error-prone DNA repair mechanisms. These error-prone repair pathways can increase mutagenesis in the cell, leading to the evolution of drug resistance. Thus, inhibiting the SOS regulator RecA has been proposed as means to delay the rise of resistance.

In this paper, the authors deleted the RecA protein from E. coli and exposed this ∆recA strain to selective levels of the beta-lactam antibiotic, ampicillin. After an 8h treatment, they washed the antibiotic away and allowed the surviving cells to recover in regular media. They then measured the minimum inhibitory concentration (MIC) of ampicillin against these treated strains. They note that after just 8 h treatment with ampicillin, the ∆recA had developed higher MICs towards ampicillin, while by contrast, wild-type cells exhibited unchanged MICs. This MIC increase was also observed in subsequent generations of bacteria, suggesting that the phenotype is driven by a genetic change.

The authors then used whole genome sequencing (WGS) to identify mutations that accounted for the resistance phenotype. Within resistant populations, they discovered key mutations in the promoter region of the beta-lactamase gene, ampC; in the penicillin-binding protein PBP3 which is the target of ampicillin; and in the AcrB subunit of the AcrAB-TolC efflux machinery. Importantly, mutations in the efflux machinery can impact the resistance towards other antibiotics, not just beta-lactams. To test this, they repeated the MIC experiments with other classes of antibiotics, including kanamycin, chloramphenicol, and rifampicin. Interestingly, they observed that the ∆recA strains pre-treated with ampicillin showed higher MICs towards all other antibiotics tested. This suggests that the mutations conferring resistance to ampicillin are also increasing resistance to other antibiotics.

The authors then performed an impressive series of genetic, microscopy, and transcriptomic experiments to show that this increase in resistance is not driven by the SOS response, but by independent DNA repair and stress response pathways. Specifically, they show that deletion of the recA reduces the bacterium's ability to process reactive oxygen species (ROS) and repair its DNA. These factors drive the accumulation of mutations that can confer resistance towards different classes of antibiotics. The conclusions are reasonably well-supported by the data, but some aspects of the data and the model need to be clarified and extended.

Strengths:

A major strength of the paper is the detailed bacterial genetics and transcriptomics that the authors performed to elucidate the molecular pathways responsible for this increased resistance. They systemically deleted or inactivated genes involved in the SOS response in E. coli. They then subjected these mutants to the same MIC assays as described previously. Surprisingly, none of the other SOS gene deletions resulted in an increase in drug resistance, suggesting that the SOS response is not involved in this phenotype. This led the authors to focus on the localization of DNA PolI, which also participates in DNA damage repair. Using microscopy, they discovered that in the RecA deletion background, PolI co-localizes with the bacterial chromosome at much lower rates than wild-type. This led the authors to conclude that deletion of RecA hinders PolI and DNA repair. Although the authors do not provide a mechanism, this observation is nonetheless valuable for the field and can stimulate further investigations in the future.

In order to understand how RecA deletion affects cellular physiology, the authors performed RNA-seq on ampicillin-treated strains. Crucially, they discovered that in the RecA deletion strain, genes associated with antioxidative activity (cysJ, cysI, cysH, soda, sufD) and Base Excision Repair repair (mutH, mutY, mutM), which repairs oxidized forms of guanine, were all downregulated. The authors conclude that down-regulation of these genes might result in elevated levels of reactive oxygen species in the cells, which in turn, might drive the rise of resistance. Experimentally, they further demonstrated that treating the ∆recA strain with an antioxidant GSH prevents the rise of MICs. These observations will be useful for more detailed mechanistic follow-ups in the future.

Weaknesses:

Throughout the paper, the authors use language suggesting that ampicillin treatment of the ∆recA strain induces higher levels of mutagenesis inside the cells, leading to the rapid rise of resistance mutations. However, as the authors note, the mutants enriched by ampicillin selection can play a role in efflux and can thus change a bacterium's sensitivity to a wide range of antibiotics, in what is known as cross-resistance. The current data is not clear on whether the elevated "mutagenesis" is driven by ampicillin selection or by a bona fide increase in mutation rate.

Furthermore, on a technical level, the authors employed WGS to identify resistance mutations in the ampicillin-treated wild-type and ∆recA strains. However, the WGS methodology described in the paper is inconsistent. Notably, wild-type WGS samples were picked from non-selective plates, while ΔrecA WGS isolates were picked from selective plates with 50 μg/mL ampicillin. Such an approach biases the frequency and identity of the mutations seen in the WGS and cannot be used to support the idea that ampicillin treatment induces higher levels of mutagenesis.

Finally, it is important to establish what the basal mutation rates of both the WT and ∆recA strains are. Currently, only the ampicillin-treated populations were reported. It is possible that the ∆recA strain has inherently higher mutagenesis than WT, with a larger subpopulation of resistant clones. Thus, ampicillin treatment might not, in fact, induce higher mutagenesis in ∆recA.

Summary of revised manuscript:

In their revisions, the authors have addressed my major concerns with additional experiments and changes to the text. Thank you!

Reviewer #3 (Public review):

Summary:

In the present work, Zhang et al investigate the involvement of the bacterial DNA damage repair SOS response in the evolution of beta-lactam drug resistance in Escherichia coli. Using a combination of microbiological, bacterial genetics, laboratory evolution, next-generation, and live-cell imaging approaches, the authors propose short-term (transient) drug resistance evolution can take place in RecA-deficient cells in an SOS response-independent manner. They propose the evolvability of drug resistance is alternatively driven by the oxidative stress imposed by accumulation of reactive oxygen species and compromised DNA repair. Overall, this is a nice study that addresses a growing and fundamental global health challenge (antimicrobial resistance).

Strengths:

The authors introduce new concepts to antimicrobial resistance evolution mechanisms. They show short-term exposure to beta-lactams can induce durably fixed antimicrobial resistance mutations. They propose this is due to compromised DNA repair and oxidative stress. Antibiotic resistance evolution under transient stress is poorly studied, so the authors' work is a nice mechanistic contribution to this field.

Weaknesses:

The authors revisions have significantly addressed weaknesses previously identified earlier in the review process.

Reviewer #1 (Public review):

Summary:

This paper shows that maternal high-fat diet during lactation changes microglia morphology in the PVN, potentially to acquire a more active state. Further, the authors reveal that PVN microglia engulf AgRP terminals in the PVN during postnatal development, a previously unrecognized behavior. A notable finding of this paper is that pharmacological reduction of microglial cells can reverse weight gain and terminal loss in the offspring under maternal high fat diet conditions, even though an increase in microglial engulfment of AgRP+ terminals was not observed, suggesting an alternative mechanism. The data support these findings, although questions remain regarding the efficacy and timing of the pharmacological microglial knockdown.

Strengths

(1) The impact of microglia on hypothalamic synaptic pruning is poorly characterized, and thus, the findings herein are especially of interest.

Weaknesses

(1) Most minor concerns were addressed during revisions, including additional details in the methods and results sections that help interpret the data as presented.

(2) The AgRP staining is unclear. For example, in Figure 2, the figure legend says "labeled AgRP terminals (red)" (Fig 2A-D) but then concludes no difference in the number of "AgRP neurons" (Fig 2J). Is this quantification of AgRP+ neurons, terminals, or both?

(3) The PLX experiments are critical to their conclusion that during lactation, microglia in the PVN sculpt AgRP inputs; however, there is no demonstration that PLX treatment effectively eliminated microglia during this postnatal window. Microglia depletion was only assessed at P55, a month past the PLX treatment window making it unclear when and by what percentage the microglia were eliminated.

Reviewer #2 (Public review):

Hypothalamic neural circuits that control body weight develop during the lactation period in rodents. Exposure to maternal high-fat diet during this period (MHFD-L) program has lasting effects on their neuroanatomical organization and function. Microglia sense environmental signals and can sculpt developing circuits by promoting or pruning synaptic connections. Here, the authors examine the contribution of microglia to the effects of MHFD-L to reduce projections from AgRP neurons in the ARH to the PVH, a critical node in circuits regulating energy balance. Using detailed histomorphometric analyses of Iba-1+ cells in the three brain regions (ARH, PVH, and BNST) at two time points (P16 and P30), the authors show that microglial volume and complexity increase, while cell numbers decrease across this period. Exposure to MHFD-L is associated with a transient increase in microglial complexity/volume at P16 in the PVH but not in the other brain regions or time points assessed. Depleting microglia using a pharmacological approach reversed effects of MHD-L on AgRP outgrowth and body weight.

Strengths:

(1) The Introduction is well-written and provides a good overview of what is known about the roles of microglia in sculpting developing circuits in the hippocampus and cortex. This provides a strong rationale for the current investigations in the hypothalamus.

(2) High-quality imaging and detailed 3-D reconstructions of Iba-1 staining in microglia are used to perform unbiased analyses of microglial complexity and to quantify the spatial relationship between microglial processes and AgRP terminals.

Weaknesses:

(1) The central claim of the manuscript is that microglia in the PVH sculpt the density of AgRP inputs to the PVH in a temporally and spatially restricted manner. While the findings of the microglial ablation experiment are consistent with this hypothesis, they do not prove causality, since their manipulations were not limited to the PVH. Further studies are needed to exclude the possibility that increased outgrowth from AgRP neurons results from direct actions in the ARH or indirect consequences of changes in growth rates.

(2) Impacts of microglial depletion were only assessed in adulthood. Given the hypothesized importance of differences in microglia at P16 and not at P30, it would be helpful to demonstrate that PLX5622 does indeed affect microglia at P16, when the circuit is most sensitive to maternal influences.

Reviewer #3 (Public review):

Summary:

The authors interrogated the putative role of microglia in determining AgRP fiber maturation in offspring exposed to a maternal high-fat diet. They found that changes in specific parts of the hypothalamus (but not in others) occur in microglia and that the effect of microglia on AgRP fibers appears to be beyond synaptic pruning, a classical function of these brain-resident macrophages.

Strengths:

The work is very strong in neuroanatomy. The images are clear and nicely convey the anatomical differences. The microglia depletion study adds functional relevance to the paper; however, the pitfalls of the technology regarding functional relevance should be discussed.

Weaknesses:

There was no attempt to functionally interrogate microglia in different parts of the hypothalamus. Morphology alone does not reflect a potential for significant signaling alterations that may occur within and between these and other cell types.

Comments on revised submission: My advice is to change the title by removing "required" and state what is interrogated and found in the paper. A more accurate title would be (for example): Implication of Microglia for Developmental Specification of AgRP Innervation in the Hypothalamus of Offspring Exposed to Maternal High-Fat Diet During Lactation.

I suggest that the authors discuss the limitations of their approach and findings, and propose future directions to address them

Reviewer #1 (Public review):

Summary:

In this manuscript, the authors used a leucine/pantothenate auxotrophic strain of Mtb to screen a library of FDA-approved compounds for their antimycobacterial activity and found significant antibacterial activity of the inhibitor semapimod. In addition to alterations in pathways, including amino acid and lipid metabolism and transcriptional machinery, the authors demonstrate that semapimod treatment targets leucine uptake in Mtb. The work presents an interesting connection between nutrient uptake and cell wall composition in mycobacteria.

Strengths:

(1) The link between the leucine uptake pathway and PDIM is interesting but has not been characterized mechanistically. The authors discuss that PDIM presents a barrier to the uptake of nutrients and shows binding of the drug with PpsB. However it is unclear why only the leucine uptake pathway was affected. We still do not know what PpsB actually does for amino acid uptake - is it a transporter? Does semapimod binding affect its activity? Does the auxotrophic Mtb have lower PDIM levels compared to wild-type Mtb?

(2) The authors show an interesting result where they observed antibacterial activity of semapimod against H37Rv only in vivo and not in vitro. Why do the authors think this is the basis of this observation? It is possible semapimod has an immunomodulatory effect on the host since leucine is an essential amino acid in mice. The authors could check pro-inflammatory cytokine levels in infected mouse lungs with and without drug treatment.

(3) The authors show that the semapimod-resistant auxotroph lacks PDIM. The conclusions would be further strengthened by including validations using PDIM mutants, including del-ppsB Mtb and other genes of the PDIM locus, whether in vivo this mutant would be more susceptible (or resistant) to semapimod treatment.

(4) Prolonged subculturing can introduce mutations in PDIM, which can be overcome by supplementing with propionate (Mullholland et al, Nat Microbiol, 2024). Did the authors also supplement their cultures with propionate? It would be interesting to see what mutations would result in Semr strains with propionate supplementation along with prolonged semapimod treatment.

Weaknesses:

I have summarized the limitations above in my comments. Overall, it would be helpful to provide more mechanistic details to study the connection between leucine uptake and PDIM.

Reviewer #2 (Public review):

Summary

This important study uncovers a novel mechanism for L-leucine uptake by M. tuberculosis and shows that targeting this pathway with 'Semapimod' interferes with bacterial metabolism and virulence. These results identify the leucine uptake pathway as a potential target to design new anti-tubercular therapy.

Strengths

The authors took numerous approaches to prove that L-leucine uptake of M. tuberculosis is an important physiological phenomenon and may be effectively targeted by 'Semapimod'. This study utilizes a series of experiments using a broad set of tools to justify how the leucine uptake pathway of M. tuberculosis may be targeted to design new anti-tubercular therapy.

Weaknesses

The study does not explain how L-leucine is taken up by M. tuberculosis, leaving the mechanism unclear. Even though 'Semapimod' binds to the PpsB protein, the relevant connection between changes in PDIM and amino acid transport remains incomplete. Also, the fact that the drug does not function on WT bacteria makes it a weak candidate to consider its usefulness for a therapeutic option.

Reviewer #3 (Public review):

Agarwal et al identified the small molecule semapimod from a chemical screen of repurposed drugs with specific antimycobacterial activity against a leucine-dependent strain of M. tuberculosis. To better understand the mechanism of action of this repurposed anti-inflammatory drug, the authors used RNA-seq to reveal a leucine-deficient transcriptomic signature from semapimod challenge. The authors then measured a decreased intracellular concentration of leucine after semapimod challenge, suggesting that semapimod disrupts leucine uptake as the primary mechanism of action. Unexpectedly, however, resistant mutants raised against semapimod had a mutation in the polyketide synthase gene ppsB that resulted in loss of PDIM synthesis. The authors believe growth inhibition is a consequence of decreased accumulation of leucine as a result of an impaired cell wall and a disrupted, unknown leucine transporter. This study highlights the importance of branched-chain amino acids for M. tuberculosis survival, and the chemical genetic interactions between semapimod and ppsB indicate that ppsB is a conditionally essential gene in a medium depleted of leucine.

The conclusions regarding the leucine and PDIM phenotypes are moderately supported by experimental data. The authors do not provide experimental evidence to support a specific link between leucine uptake and impaired PDIM production. Additional work is needed to support these claims and strengthen this mechanism of action.

(1) Since leucine uptake and PDIM synthesis are important concepts of the manuscript, experiments would benefit from exploring other BCAAs to know if the phenotypes observed are specific to leucine, and adding additional strains to the 2D TLC experiments to provide confidence in the absence of the PDIM band.

(2) The intriguing observation that wild-type H37Rv is resistant to semapimod but the leucine-auxotroph is sensitive should be further explored. If the authors are correct and semapimod does inhibit leucine uptake through a specific transporter or disrupted cell wall (PDIM synthesis), testing semapimod activity against the leucine-auxotroph in various concentrations of BCAAs could highlight the importance of intracellular leucine. H37Rv is still able to synthesize endogenous leucine and is able to circumvent the effect of semapimod.

Reviewer #1 (Public review):

Summary:

The manuscript finds a negative relationship between tuberculin skin test-induced type I interferon activity with chest X-ray tuberculosis severity in humans. This evidence is between incomplete and solid. It needs a bioinfomatics/transcriptomics reviewer to make a more insightful judgement. The manuscript demonstrates a convincing role for Stat2 in controlling Mycobacterium marinum infection in zebrafish embryos, incomplete data are presented linking reduced leukocyte recruitment to the infection susceptibility phenotype.

Strengths:

(1) An interesting analysis of TST response correlated with chest X-ray pathology.

(2) Novel data on a protective role for Stat2 in a natural host-mycobacterial species infection pairing.

Weaknesses:

(1) The transcriptional modules are very large sets of genes that do not present a clear picture of what is actually being measured relative to other biological pathways.

(2) The link between infection-Stat2-leukocyte recruitment and containment of infection is plausible, but lacks a specific link to the first part of the manuscript.

Major concerns

(1) Line 158: The two transcriptional modules should be placed in the context of other DEG patterns. The macrophage type I interferon module, in particular, is quite large (361 genes). Can this be made more granular in terms of type I IFN ligands and STAT2-dependent genes?

(2) The ifnphi1 injection into mxa:mCherry stat2 crispants is a nice experiment to demonstrate loss of type I IFN responsiveness. Further data is required to demonstrate if important mycobacterial control pathways (IFNy, TNF, il6?, etc) are intact in stat2 crispants before being able to conclude that these phenotypes are specific to type I IFN.

Reviewer #2 (Public review):

Summary:

This study shows that type I interferon (IFN-I) signaling helps protect against mycobacterial infection. Using human gene expression data and a zebrafish model, the authors find that reduced IFN-I activity is linked to more severe disease. They also show that zebrafish lacking the IFN-I signaling gene stat2 are more vulnerable to infection due to poor macrophage migration. These results suggest a protective role for IFN-I in mycobacterial disease, challenging previous findings from other animal models.

Strengths:

Strengths of the manuscript include the use of human clinical samples to support relevance to disease, along with a genetically tractable zebrafish model that enables mechanistic insight.

Weaknesses:

(1) The manuscript presents intriguing human data showing an inverse correlation between IFN-I gene signatures and TB disease, but the findings remain correlative and may be cohort-specific. Given that the skin is not a primary site of TB and is relatively immunotolerant, the biological relevance of downregulated IFN-I-related genes in this tissue to systemic or pulmonary TB is unclear.

(2) The reliance on stat2 CRISPants in zebrafish offers a limited view of IFN-I signaling. Including additional crispant lines targeting other key regulators (e.g., ifnar1, tyk2, irf3, irf7) would strengthen the interpretation and clarify whether the observed effects reflect broader IFN-I pathway disruption.

(3) The conclusion that IFN-I is protective contrasts with established findings from murine and non-human primate models, where IFN-I is often detrimental. While the authors highlight species differences, the lack of functional human data and reliance on M. marinum in zebrafish limit the translational relevance. A more balanced discussion addressing these discrepancies would improve the manuscript.

(4) Quantification of bacterial burden using fluorescence intensity alone may not accurately reflect bacterial viability. Complementary methods, such as qPCR for bacterial DNA, would provide a more robust assessment of antimicrobial activity.

(5) Finally, the authors should clarify whether impaired macrophage recruitment in stat2 crispants results from defects in chemotaxis, differentiation, or survival, and address discrepancies between their human blood findings and prior studies.

Reviewer #3 (Public review):

Summary:

In this manuscript, the authors presented an interesting study providing an insight into the role of Type-I interferon responses in tuberculosis (TB) pathogenesis by combining transcriptome analysis of PBMCs and TST from tuberculosis patients. The zebrafish model was used to identify the changes in the innate immune cell population of macrophages and neutrophils. The findings suggested that Type-I interferon signatures inversely correlated with disease severity in the TST transcriptome data. The authors validated the observations by CRISPR-mediated disruption of stat2 (a critical transcription factor for type I interferon signaling) in zebrafish larvae, showing increased susceptibility to M. marinum infection. Traditionally, type-I interferon responses have been viewed as detrimental in mycobacterial infections, with studies suggesting enhanced susceptibility in certain mouse models. The study tried to identify and further characterize the understanding of the role of type-I interferons in TB.

Strengths:

Traditionally, type-I interferon responses have been viewed as detrimental in mycobacterial infections, with studies suggesting enhanced susceptibility in certain mouse models. The study tried to further understand the role of type-I interferons in TB pathogenesis.

Weaknesses:

Though the study showed an inverse correlation of Type-I interferon with radiological features of TB, the molecular mechanism is largely unexplored in the study, which is making it difficult to understand the basis of the results shown in the manuscript by the authors.

Reviewer #1 (Public review):

Summary:

A fundamental technique for the identification of peptide-specific CD8 T cells is the use of fluorophore-conjugated and peptide loaded MHC tetramers. Classically, refolding of specific peptides with MHC monomers can be labour intensive, and not optimal for screening large numbers of different peptides. Hence, UV-exchanged tetramers have been developed to upscale this, however, still has some associated challenges such as UV-mediated damage to peptide complexes. Here, Pothast, C.R. et al demonstrate the efficacy of using temperature exchanged tetramers for the prevalent alleles HLA-A*03:01, A*11:01, B*07:02, and C*07:02. Building upon their previous work with HLA-A*02:01, H-2Kb, and HLA-E. They first demonstrate the complex stability of tetramers with different affinity peptides at high temperature, showing complex destabilisation can be rescued with higher affinity peptides. This is followed by an optimisation of peptide exchange temperatures, tailored for each allele. The authors then demonstrate successful binding to clonal T cell lines, and then a step further with viral peptides against PBMCs from individuals with confirmed infection history. For the latter they compare to conventional tetramers and demonstrate comparable signal.<br /> Due to the prevalence of these 4 alleles, the ease-of-handling, and short time requirements, these tetramers are likely to show high utility.

Strengths:

The manuscript is well-written and the results are solid, although more detail may add clarity to some of the results, in particular Figures 1 and 2. Other than the points reported below, the study uses accurate controls to demonstrate the specificity of the tetramers, and the data are convincing.

Overall, the interpretation of the results is accurate, and the discussion is thorough. Additional comments may be included to cover potential tetramer batch variability and differences in the stability of different alleles. Specifically, whether certain alleles require higher-affinity peptides to be stable, compared to others.

Weaknesses:

The authors demonstrate the equivalence of temperature-exchanged tetramers to conventional ones, however, as they are an advancement on UV-exchange, it would be useful to show data on how their stability, exchange efficacy, and binding to T cell lines compare to UV-based tetramers. It would be supportive to show that temperature does not impact fluorophore intensity as well.

Reviewer #2 (Public review):

Summary:

The majority of CD8+ T cell responses rely on the proper presentation of antigens through stable MHC-I (but not requiring a stable immunological synapse). This work highlights a new approach to build an array of stable peptide MHC-I using temperature exchange, which can be used to identify antigen-specific CD8+ T cells.

Strengths:

In this work, the authors have proposed an alternative method to reload the peptide MHC-I molecule. Their temperature-exchange approach is distinct from current reloadable peptide MHC technologies involving photolabile peptide, empty MHC-I (Nat Commun 11, 1314 (2020). https://doi.org/10.1038/s41467-020-14862-4), tapasin/TAPBPR chaperone-assisted (eLife 7:e40126.), enzyme exchangeable (WO2020226570) and small alcohol (Curr Res Immunol. 2022 Aug 18;3:167-174. doi: 10.1016/j.crimmu.2022.08.002) approaches.

Weaknesses:

However, the proposed temperature-exchange approach does not substantially improve the quality of antigen-specific T cells that can be identified using the photolabile peptide MHC-I molecules.

The time saved using the temperature-exchange protocol may not be a pull factor as the photolabile peptide MHC-I approach is not unreasonably laborious.

Reviewer #3 (Public review):

Summary:

The study by Pothast and colleagues outlines an extension of their previously described temperature-based MHC-I peptide exchange method on 4 common HLA alleles, to enable the generation of peptide/MCH-I tetramers for characterization of antigen-specific T cells by flow cytometry.

Strengths:

This work outlines a protocol for generating MHC-I tetramers on 4 common HLA allotypes, which can then be applied to monitor T cell responses by flow cytometry studies. The work provides conditional ligands for exchange on each HLA and demonstrates proof of concept studies using clonotypic T cells and CD8+ PBMCs.

The results support that the temperature-exchanged tetramers can perform similarly to conventional tetramers in some settings.

Weaknesses:

Given that there are several proposed methodologies addressing the same task (including UV-mediated, disulfide-bond based stabilization of empty MHC-I conformers, and chaperone-based methods), the relevance of the proposed temperature-mediated technology is questionable.

More specifically, important limitations of the study include:

(1) A lack of quantification of exchanged molecules relative to molecules that retain the original placeholder peptides, or completely empty molecules present in the same sample.

(2) A lack of validation that peptide exchange has occurred in the absence of a reporter T cell line appears to be a significant limitation of the methodology for antigen / T cell discovery.

(3) The sub-optimal exchange efficiency relative to conventional prepared pMHC-I molecules, shown in Figure 4, is a significant limitation of the approach.

(4) There are no data to support that exchange proceeds through the generation of empty molecules during the temperature cycle, or by peptide binding on empty molecules that are already present in the sample. Understanding the mechanism of exchange is important for the necessary improvements to the methodology.

(5) It is possible that the temperature cycle causes protein aggregation or other irreversible changes to the sample - this should be explicitly quantified and addressed in the paper, since misfolded MHC-I molecules can lead to high levels of background staining.

(6) These potential limitations should limit detection of low-affinity/low-avidity interactions between TCRs and their cognate pMHC antigens - this should be addressed explicitly in a model antigen setting.

(7) The approach appears to be limited to the HLAs showing high thermal stability, which have been explored in this study. However, a large fraction of HLAs show sub-optimal thermal stabilities. It seems that explicit validation of peptide exchange would be required for any new HLA allele introduced into this process.

(8) Whether the approach can be used to load suboptimal peptides with lower thermal stabilities that are emerging immunotherapy targets is not addressed in the present study.

Because of these limitations, the present manuscript does not conclusively support the claim that temperature-based exchange can be used as a robust methodology to generate pMHC-I tetramers with desired peptide specificities.

As a result, the scope of applications using these suboptimal exchanged pHLA tetramers is limited, and should be addressed with further improvements of the methodology, including better characterization of exchange efficiency, demonstration of functionality across a broader range of HLA allotypes with varying thermal stability profiles, and validation with clinically relevant low-affinity peptides that would strengthen the potential utility of this approach in immunotherapy development and basic T cell biology research.

Reviewer #1 (Public review):

Summary:

This study focuses on the bacterial metabolite TMA, generated from dietary choline. These authors and others have previously generated foundational knowledge about the TMA metabolite TMAO, and its role in metabolic disease. This study extends those findings to test whether TMAO's precursor, TMA, and its receptor TAAR5 are also involved and necessary for some of these metabolic phenotypes. They find that mice lacking the host TMA receptor (Taar5-/-) have altered circadian rhythms in gene expression, metabolic hormones, gut microbiome composition, and olfactory and innate behavior. In parallel, mice lacking bacterial TMA production or host TMA oxidation have altered circadian rhythms.

Strengths:

These authors use state-of-the-art bacterial and murine genetics to dissect the roles of TMA, TMAO, and their receptor in various metabolic outcomes (primarily measuring plasma and tissue cytokine/gene expression). They also follow a unique and unexpected behavioral/olfactory phenotype. Statistics are impeccable.

Weaknesses:

Enthusiasm for the manuscript is dampened by some ambiguous writing and the presentation of ideas in the introduction, both of which could easily be improved upon revision.

Reviewer #2 (Public review):

Summary:

In the manuscript by Mahen et al., entitled "Gut Microbe-Derived Trimethylamine Shapes Circadian Rhythms Through the Host Receptor TAAR5," the authors investigate the interplay between a host G protein-coupled receptor (TAAR5), the gut microbiota-derived metabolite trimethylamine (TMA), and the host circadian system. Using a combination of genetically engineered mouse and bacterial models, the study demonstrates a link between microbial signaling and circadian regulation, particularly through effects observed in the olfactory system. Overall, this manuscript presents a novel and valuable contribution to our understanding of host-microbe interactions and circadian biology. However, several sections would benefit from improved clarity, organization, and mechanistic depth to fully support the authors' conclusions.

Strengths:

(1) The manuscript addresses an important and timely topic in host-microbe communication and circadian biology.

(2) The studies employ multiple complementary models, e.g., Taar5 knockout mice, microbial mutants, which enhance the depth of the investigation.

(3) The integration of behavioral, hormonal, microbial, and transcript-level data provides a multifaceted view of the observed phenotype.

(4) The identification of olfactory-linked circadian changes in the context of gut microbes adds a novel perspective to the field.

Weaknesses:

While the manuscript presents compelling data, several weaknesses limit the clarity and strength of the conclusions.

(1) The presentation of hormonal, cytokine, behavioral, and microbiome data would benefit from clearer organization, more detailed descriptions, and functional grouping to aid interpretation.

(2) Some transitions-particularly from behavioral to microbiome data-are abrupt and would benefit from better contextual framing.

(3) The microbial rhythmicity analyses lack detail on methods and visualization, and the sequencing metadata (e.g., sample type, sex, method) are not clearly stated.

(4) Several figures are difficult to interpret due to dense layouts or vague legends, and key metabolites and gene expression comparisons are either underexplained or not consistently assessed across models.

(5) Finally, while the authors suggest a causal role for TAAR5 and its ligand in circadian regulation, the current data remain correlative; mechanistic experiments or stronger disclaimers are needed to support these claims.

Reviewer #3 (Public review):

Summary:

Deletion of the TMA-sensor TAAR5 results in circadian alterations in gene expression, particularly in the olfactory bulb, plasma hormones, and neurobehaviors.

Strengths:

Genetic background was rigorously controlled.

Comprehensive characterization.

Weaknesses:

The weaknesses identified by this reviewer are minor.

Overall, the studies are very nicely done. However, despite careful experimentation, I note that even the controls vary considerably in their gene expression, etc, across time (eg, compare control graphs for Cry 1 in IB, 4B). It makes me wonder how inherently noisy these measurements are. While I think that the overall point that the Taar5 KO shows circadian changes is robust, future studies to dissect which changes are reproducible over the noise would be helpful.

Impact:

These data add to the growing literature pointing to a role for the TMA/TMAO pathway in olfaction and neurobehavioral.

Reviewer #1 (Public review):

Summary:

This manuscript analyses primarily the effects of deleting the TgfbR1 and TgfbR2 receptors from endothelial cells at postnatal stages of vascular development and blood-retina barrier maturation in the retina. The authors find that deletion of these receptors affects vascular development in the retina, but importantly, it affects the infiltration of immune cells across the vessels in the retina. The findings demonstrate that Tgfb signaling through TgfbR1/R2 heterodimers regulates primarily the immune phenotypes of endothelial cells in addition to regulating vascular development. The data provided by the authors provide a solid support for their conclusions.

Strengths:

(1) The manuscript uses a variety of elegant genetic studies in mice to analyze the role of TgfbR1 and TgfbR2 receptors in endothelial cells at postnatal stages of vascular development and blood-retina barrier maturation in the retina.

(2) The authors provide a nice comparison of the vascular phenotypes in endothelial-specific knockout of TgfbR1 and TgfbR2 in the retina (and to a lesser degree in the brain) with those from Npd KO mice (loss of Ndp/Fzd signaling) or loss of VEGF-A signaling to dissect the specific roles of Tgf signaling for vascular development in the retina.

(3) The snRNAseq data of vessel segments from the brains of WT versus TgfbR1 -iECKO mice provides a nice analysis of pathways and transcripts that are regulated by Tgfb signaling in endothelial cells.

Weaknesses:

(1) The authors claim that choroidal neovascular tuft phenotypes are similar in TgfbrR1 KO and TgfbrR2 KO mice. However, the phenotypes look more severe in the TgfbrR1 KO rather than TgfbrR2 KO mice. Can the authors show a quantitative comparison of the number of choroidal neovascular tufts per whole eye cross-section in both genotypes?

(2) In the analysis of Sulfo-NHS-Biotin leakage in the retina to assess blood-retina barrier maturation. The authors claim that there is increased vascular leakage in the TgfbR1 KO mice. However, it does not seem like Sulfo-NHS-biotin is leaking outside the vessels. Therefore, it cannot be increased vascular permeability. Can the authors provide a detailed quantification of the leakage phenotype?

(3) The immune cell phenotyping by snRNAseq is premature, as the number of cells is very small. The authors should sort for CD45+ cells and perform single-cell RNA sequencing.

(4) The analysis of BBB leakage phenotype in TgfbR1 KO mice needs to be more detailed and include tracers as well as serum IgG leakage.

(5) A previous study (Zarkada et al., 2021, Developmental Cell) showed that EC-deletion of Alk5 affects the D tip cells. The phenotypes of those mice look very similar to those shown for TgfbrR1 KO mice. Are D-tip cells lost in these mutants by snRNAseq?

Reviewer #2 (Public review):

Summary:

The authors meticulously characterized EC-specific Tgfbr1, Tgfbr2, or double knockout in the retina, demonstrating through convincing immunostaining data that loss of TGF-β signaling disrupts retinal angiogenesis and choroidal neovascularization. Compared to other genetic models (Fzd4 KO, Ndp KO, VEGF KO), the Tgfbr1/2 KO retina exhibits the most severe immune cell infiltration. The authors proposed that TGF-β signaling loss triggers vascular inflammation, attracting immune cells - a phenotype specific to CNS vasculature, as non-CNS organs remain unaffected.

Strengths:

The immunostaining results presented are clear and robust. The authors performed well-controlled analyses against relevant mouse models. snRNA-seq corroborates immune cell leakage in the retina and vascular inflammation in the brain.

Weaknesses:

The causal link between TGF-β loss, vascular inflammation, and immune infiltration remains unresolved. The authors' model posits that EC-specific TGF-β loss directly causes inflammation, which recruits immune cells. However, an alternative explanation is plausible: Tgfbr1/2 KO-induced developmental defects (e.g., leaky vessels) permit immune extravasation, subsequently triggering inflammation. The observations that vein-specific upregulation of ICAM1 staining and the lack of immune infiltration phenotypes in the non-CNS tissues support the alternative model. Late-stage induction of Tgfbr1/2 KO (avoiding developmental confounders) could clarify TGF-β's role in retinal angiogenesis versus anti-inflammation.

Reviewer #1 (Public review):

Summary:

This paper examines how geometric regularities in abstract shapes (e.g., parallelograms, kites) are perceived and processed in the human brain. The manuscript contains multimodal data (behavior, fMRI, MEG) from adults and additional fMRI data from 6-year-old children. The key findings show that (1) processing geometric shapes lead to reduced activity in ventral areas in comparison to complex stimuli and increased activity in intraparietal and inferior temporal regions, (2) the degree of geometric regularity modulates activity in intraparietal and inferior temporal regions, (3) similarity in neural representation of geometric shapes can be captured early by using CNN models and later by models of geometric regularity. In addition to these novel findings, the paper also includes a replication of behavioral data, showing that the perceptual similarity structure amongst the geometric stimuli used can be explained by a combination of visual similarities (as indexed by a feedforward CNN model of the ventral visual pathway) and geometric features.

Strengths:

(1) The study incorporates multi-modal data that uses more than one task and different populations of participants (adults and children).

(2) It replicates behavioral findings of an earlier study in a larger cohort.

(3) The paper comes with openly accessible code in a well-documented GitHub repository, and the data will be published with the paper on OpenNeuro.

Weaknesses:

I wonder how task difficulty and linguistic labels interact with the current findings. Based on the behavioral data, shapes with more geometric regularities are easier to detect when surrounded by other shapes. Do shape labels that are readily available (e.g., "square") help in making accurate and speedy decisions? Can the sensitivity to geometric regularity in intraparietal and inferior temporal regions be attributed to differences in task difficulty? Similarly, are the MEG oddball detection effects that are modulated by geometric regularity also affected by task difficulty?

Reviewer #2 (Public review):

Summary:

The current study seeks to understand the neural mechanisms underlying geometric reasoning. Using fMRI with both children and adults, the authors found that contrasting simple geometric shapes with naturalistic images (faces, tools, houses) led to responses in the dorsal visual stream, rather than ventral regions that are generally thought to represent shape properties. The authors followed up on this result using computational modeling and MEG to show that geometric properties explain distinct variance in the neural response beyond what is captured by a CNN.

Strengths:

These findings contribute much-needed neural and developmental data to the ongoing debate regarding shape processing in the brain and offer additional insights into why CNNs may have difficulty with shape processing. The motivation and discussion for the study are appropriately measured, and I appreciate the authors' use of multiple populations, neuroimaging modalities, and computational models to explore this question.

Weaknesses:

Given that the primary take away from this study is that geometric shape information is found in the dorsal stream, rather than the ventral stream there is very little there is very little discussion of prior work in this area (for reviews, see Freud et al., 2016; Orban, 2011; Xu, 2018). Indeed, there is extensive evidence of shape processing in the dorsal pathway in human adults (Freud, Culham, et al., 2017; Konen & Kastner, 2008; Romei et al., 2011), children (Freud et al., 2019), patients (Freud, Ganel, et al., 2017), and monkeys (Janssen et al., 2008; Sereno & Maunsell, 1998; Van Dromme et al., 2016), as well as the similarity between models and dorsal shape representations (Ayzenberg & Behrmann, 2022; Han & Sereno, 2022).

The presence of activation in aIPS led the authors to interpret their results to mean that geometric reasoning draws on the same processes as mathematical thinking. However, there is not enough evidence in the current study to support this claim.

Reviewer #3 (Public review):

Summary:

The authors report converging evidence from several brain-imaging techniques that geometric figures, notably quadrilaterals, are processed differently in visual (lower activation) and spatial (greater) areas of the human brain than representative figures. Comparison of mathematical models to fit activity for geometric figures shows the best fit for abstract geometric features like parallelism and symmetry. The brain areas active for geometric figures are also active in processing mathematical concepts, even in blind mathematicians, linking geometric shapes to abstract math concepts. The effects are stronger in adults than in 6-year-old Western children. Similar phenomena do not appear in great apes, suggesting that this is uniquely human and developmental.

Strengths:

Multiple converging techniques of brain imaging and testing of mathematical models. Careful reasoning at every step of research and presentation of research, anticipating and addressing possible reservations. Connecting these findings to other findings, brain, behavior, and historical/anthropological, to suggest broad and important fundamental connections between abstract visual-spatial forms and mathematical reasoning, further suggesting visual-spatial origins of mathematical reasoning.

Weaknesses:

Perhaps the manuscript could emphasize that the areas recruited by geometric figures but not objects are spatial, with reduced processing in visual areas. It also seems important to say that the images of real objects are interpreted as representations of 3D objects, as they activate the same visual areas as real objects. By contrast, the images of geometric forms are not interpreted as representations of real objects but rather perhaps as 2D abstractions. The authors use the term "symbolic." That use of that term could usefully be expanded here.

Pigeons have remarkable visual systems. According to my fallible memory, Herrnstein investigated visual categories in pigeons. They can recognize individual people from fragments of photos, among other feats. I believe pigeons failed at geometric figures and also at cartoon drawings of things they could recognize in photos. This suggests they did not interpret line drawings of objects as representations of objects.

Categories are established in part by contrast categories; are quadrilaterals, triangles, and circles different categories?

It would be instructive to investigate stimuli that are on a continuum from representational to geometric, e.g., table tops or cartons under various projections, or balls or buildings that are rectangular or triangular. Building parts, inside and out. like corners. Objects differ from geometric forms in many ways: 3D rather than 2D, more complicated shapes, and internal texture. The geometric figures used are flat, 2-D, but much geometry is 3-D (e. g. cubes) with similar abstract features. The feature space of geometry is more than parallelism and symmetry; angles are important, for example. Listing and testing features would be fascinating. Similarly, looking at younger or preferably non-Western children, as Western children are exposed to shapes in play at early ages.

What in human experience but not the experience of close primates would drive the abstraction of these geometric properties? It's easy to make a case for elaborate brain processes for recognizing and distinguishing things in the world, shared by many species, but the case for brain areas sensitive to processing geometric figures is harder. The fact that these areas are active in blind mathematicians and that they are parietal areas suggests that what is important is spatial far more than visual. Could these geometric figures and their abstract properties be connected in some way to behavior, perhaps with fabrication and construction as well as use? Or with other interactions with complex objects and environments where symmetry and parallelism (and angles and curvature--and weight and size) would be important? Manual dexterity and fabrication also distinguish humans from great apes (quantitatively, not qualitatively), and action drives both visual and spatial representations of objects and spaces in the brain. I certainly wouldn't expect the authors to add research to this already packed paper, but raising some of the conceptual issues would contribute to the significance of the paper.

Reviewer #1 (Public review):

Lipid transfer proteins (LTPs) play a crucial role in the intramembrane lipid exchange within cells. However, the molecular mechanisms that govern this activity remain largely unclear. Specifically, the way in which LTPs surmount the energy barrier to extract a single lipid molecule from a lipid bilayer is not yet fully understood. This manuscript investigates the influence of membrane properties on the binding of Ups1 to the membrane and the transfer of phosphatidic acid (PA) by the LTP. The findings reveal that Ups1 shows a preference for binding to membranes with positive curvature. Moreover, coarse-grained molecular dynamics simulations indicate that positive curvature decreases the energy barrier associated with PA extraction from the membrane. Additionally, lipid transfer assays conducted with purified proteins and liposomes in vitro demonstrate that the size of the donor membrane significantly impacts lipid transfer efficiency by Ups1-Mdm35 complexes, with smaller liposomes (characterized by high positive curvature) promoting rapid lipid transfer.

This study offers significant new insights into the reaction cycle of phosphatidic acid (PA) transfer by Ups1 in mitochondria. Notably, the authors present compelling evidence that, alongside negatively charged phospholipids, positive membrane curvature enhances lipid transfer - an effect that is particularly relevant at the mitochondrial outer membrane. The experiments are technically robust, and my primary feedback pertains to the interpretation of specific results.

(1) The authors conclude from the lipid transfer assays (Figure 5) that lipid extraction is the rate-limiting step in the transfer cycle. While this conclusion seems plausible, it should be noted that the authors employed high concentrations of Ups1-Mdm35 along with less negatively charged phospholipids in these reactions. This combination may lead to binding becoming the rate-limiting factor. The authors should take this point into consideration. In this type of assay, it is challenging to clearly distinguish between binding, lipid extraction, and membrane dissociation as separate processes.

(2) The authors should discuss that variations in the size of liposomes will also affect the distance between them at a constant concentration, which may affect the rate of lipid transfer. Therefore, the authors should determine the average size and size distribution of liposomes after sonication (by DLS or nanoparticle analyzer, etc.).

(3) The authors use NBD-PA in the lipid transfer assays. Does the size of the donor liposomes affect the transfer of NBD-PA and DOPA similarly? Since NBD-labeled lipids are somewhat unstable within lipid bilayers (as shown by spontaneous desorption in Figure 5B), monitoring the transfer of unlabeled PA in at least one setting would strengthen the conclusion of the swap experiments.

(4) The present study suggests that membrane domains with positive curvature at the outer membrane may serve as starting points for lipid transport by Ups1-Mdm35. Is anything known about the mechanisms that form such structures? This should be discussed in the text.

Reviewer #2 (Public review):

Summary:

Lipid transfer between membranes is essential for lipid biosynthesis across different organelle membranes. Ups1-Mdm35 is one of the best-characterized lipid transfer proteins, responsible for transferring phosphatidic acid (PA) between the mitochondrial outer membrane (OM) and inner membrane (IM), a process critical for cardiolipin (CL) synthesis in the IM. Upon dissociation from Mdm35, Ups1 binds to the intermembrane space (IMS) surface of the OM, extracts a PA molecule, re-associates with Mdm35, and moves through the aqueous IMS to deliver PA to the IM. Here, the authors analyzed the early steps of this PA transfer - membrane binding and PA extraction - using a combination of in vitro biochemical assays with lipid liposomes and purified Ups1-Mdm35 to measure liposome binding, lipid transfer between liposomes, and lipid extraction from liposomes. The authors found that membrane curvature, a previously overlooked property of the membrane, significantly affects PA extraction but not PA insertion into liposomes. These findings were further supported by MD simulations.

Strengths:

The experiments are well-designed, and the data are logically interpreted. The present study provides an important basis for understanding the mechanism of lipid transfer between membranes.　

Weaknesses:

The physiological relevance of membrane curvature in lipid extraction and transfer still remains open.

Reviewer #3 (Public review):

The manuscript by Sadeqi et al. studies the interactions between the mitochondrial protein Ups1 and reconstituted membranes. The authors apply synthetic liposomal vesicles to investigate the role of pH, curvature, and charge on the binding of Ups1 to membranes and its ability to extract PA from them. The manuscript is well written and structured. With minor exceptions, the authors provide all relevant information (see minor points below) and reference the appropriate literature in their introduction. The underlying question of how the energy barrier for lipid extraction from membranes is overcome by Ups1 is interesting, and the data presented by the authors could offer a valuable new perspective on this process. It is also certainly a challenging in vitro reconstitution experiment, as the authors aim to disentangle individual membrane properties (e.g., curvature, charge, and packing density) to study protein adsorption and lipid transfer. I have one major suggestion and a few minor ones that the authors might want to consider to improve their manuscript and data interpretation:

Major Comments:

The experiments are performed with reconstituted vesicles, which are incubated with recombinant protein variants and quantitatively assessed in flotation and pelleting assays. According to the Materials and Methods section, the lipid concentration in these assays is kept constant at 5 µM. However, the authors change the size of the vesicles to tune their curvature. Using the same lipid concentration but varying vesicle sizes results in different total vesicle concentrations. Moreover, larger vesicles (produced by freeze-thawing and extrusion) tend to form a higher proportion of multilamellar vesicles, thus also altering the total membrane area available for binding. Could these differences in the experimental system account for the variation in binding? To address this, the authors would need to perform the experiments either under saturation (excess protein) conditions or find an experimental approach to normalize for these differences.

Reviewer #1 (Public review):

Summary:

The nuclear protein SATB1 was originally identified as a protein of the 'nuclear matrix', an aggregate of nuclear components that arose upon extracting nuclei with high salt. While the protein was assumed to have a global function in chromatin organization, it has subsequently been linked to a variety of pathological conditions, notably cancer. The mapping of the factor by conventional ChIP procedures showed strong enrichment in active, accessible chromatin, suggesting a direct role in gene regulation, perhaps in enhancer-promoter communication. These findings did not explain why SATB1-chromatin interaction resisted the 2 M salt extraction during early biochemical fractionation of nuclei.

The authors, who have studied SATB1 for many years, now developed an unusual variation of the ChIP procedure, in which they purify crosslinked chromatin by centrifugation through 8 M urea. Remarkably, while they lose all previously mapped signals for SATB1 in active chromatin, they now gain many binding events in silent regions of the genome, represented by lamin-associated domains (LADs).

SATB1 had previously been shown by the authors and others to bind to DNA with special properties, termed BUR for 'base-unpairing regions'). BURs are AT-rich and apparently enriched in equally AT-rich LADs. The 'urea-ChIP' pattern is essentially complementary to the classical ChIP pattern. The authors now speculate that the previously known SATB1 binding pattern determined by standard ChIP, which does not overlap BURs particularly well, is due to indirect chromatin binding, whereas they consider the urea-ChIP profile, which fits better to the BUR distribution on the chromosome, to be due to direct binding.

Building on the success with urea-ChIP the authors adapted the 4C-procedure of chromosome conformation mapping to work with urea-purified chromatin. The data suggest a model according to which BUR-bound SATB1 mediates long-distance interaction between active loci and some kind of scaffold structure formed by SATB1. Because cell type-specific differences are observed, they suggest that the SATB1 interactions are functionally relevant.

Strengths:

Given the unusual findings of essentially mutually exclusive 'standard ChIP' and 'urea-ChIP' profiles, the authors conducted many appropriate controls. They showed that all SATB1 peaks in urea-ChIP and 96% of peaks in standard-ChIP represent true signals, as they are not observed in a SATB1 knockout cell line. They also show that the urea-ChIP and standard ChIP yield similar profiles for CTCF and polycomb complex subunits. The data appear reproducible judged by at least two replicates and triangulation. The SATB1 KO cells provide a nice control for the specificity of signals, including those that arise from their elaborately modified 4C protocol.

In their revised manuscript the authors provide relevant background information concerning the effect of urea on the denaturation of macromolecules. Importantly, they argue convincingly that urea does not denature DNA under their conditions.

Weaknesses:

Despite the authors' efforts to explain their findings along with a lot of background information, some readers may be left confused due to the complexity of the system. BURs are found enriched in LADs, but are also present in active chromatin. SATB1 binds a subset of BURs, but the reason for discrimination remains unclear. SATB1 appears to bind chromatin in at least two modes with differing diffusion properties and exactly how this relates to the indirect and direct chromatin binding modes is mechanistically unclear.

The authors resort to the term 'SATB1-enriched subnuclear structure' to describe the profile gained through denaturing ChIP, thus avoiding strong statements about involvement of known nuclear structures (such as LADs or heterochromatin) and about functional implications.

The authors acknowledge a potential for RNA to be involved in modulating SATB1 interactions with chromatin, but leave this for future investigation.

Comment on revised version:

The authors revised their manuscript to my satisfaction.

Reviewer #2 (Public review):

Summary:

This study describes the key observation that SATB1 binds directly to so-called BUR elements. This is in contrast to several other reports describing SATB1 binding to promoters and enhancers. This discrepancy is explained by the authors to depend on the features of the ChIP technique being used. Urea-ChIP, innovated by the authors, strips off protein-protein interactions that compound conventional ChIP methods. The authors convincingly make the case that SATB1 and a key genome organiser, CTCF, largely bind different sites, as particularly evident in Figure 2A. In contrast, standard ChIP shows considerable overlap between their sites (Figure 2-figure supplement 1). The report documents convincingly that SATB1 partitions the genome independent of so-called TADs to influence expression patterns. SATB1 controls long-range interactions in thymocytes, and knock down of SATB1 does not affect the TAD patterns.

Strengths:

A new and innovative adaptation of ChIP-seq (urea ChIP-seq) has enabled the authors to successfully question existing data on the patterns of SATB1 binding to the genome. The authors provide a wealth of data to reinforce their claims. This report thus rectifies misconceptions about SATB1 function, which is particularly important given its role in metastasising cancer cells.

Weaknesses:

None

Reviewer #1 (Public review):

Summary:

This work studies representations in a network with one recurrent layer and one output layer that needs to path-integrate so that its position can be accurately decoded from its output. To formalise this problem, the authors define a cost function consisting of the decoding error and a regularisation term. They specify a decoding procedure that, at a given time, averages the output unit center locations, weighted by the activity of the unit at that time. The network is initialised without position information, and only receives a velocity signal (and a context signal to index the environment) at each timestep, so to achieve low decoding error it needs to infer its position and keep it updated with respect to its velocity by path integration.

The authors take the trained network and let it explore a series of environments with different geometries while collecting unit activities to probe learned representations. They find localised responses in the output units (resembling place fields) and border responses in the recurrent units. Across environments, the output units show global remapping and the recurrent units show rate remapping. Stretching the environment generally produces stretched responses in output and recurrent units. Ratemaps remain stable within environments and stabilise after noise injection. Low-dimensional projections of the recurrent population activity forms environment-specific clusters that reflect the environment's geometry, which suggests independent rather than generalised representations. Finally, the authors discover that the centers of the output unit ratemaps cluster together on a triangular lattice (like the receptive fields of a single grid cell), and find significant clustering of place cell centers in empirical data as well.

The model setup and simulations are clearly described, and are an interesting exploration of the consequences of a particular set of training requirements - here: path integration and decodability. But it is not obvious to what extent the modelling choices are a realistic reflection of how the brain solves navigation. Therefore, it is not clear whether the results generalize beyond the specifics of the setup here.

Strengths:

The authors introduce a very minimal set of model requirements, assumptions, and constraints. In that sense, the model can function as a useful 'baseline', that shows how spatial representations and remapping properties can emerge from the requirement of path integration and decodability alone. Moreover, the authors use the same formalism to relate their setup to existing spatial navigation models, which is informative.

The global remapping that the authors show is convincing and well-supported by their analyses. The geometric manipulations and the resulting stretching of place responses, without additional training, are interesting. They seem to suggest that the recurrent network may scale the velocity input by the environment dimensions so that the exact same path integrator-output mappings remain valid (but maybe there are other mechanisms too that achieve the same).

The simulations and analyses in the appendices serve as insightful controls for the main results.

The clustering of place cell peaks on a triangular lattice is intriguing, given there is no grid cell input. It could have something to do with the fact that a triangular lattice provides optimal coverage of 2d space? The included comparison with empirical data is valuable as a first exploration, showing a promising example, but doesn't robustly support the modelling results.

Weaknesses:

The navigation problem that needs to be solved by the model is a bit of an odd one. Without any initial position information, the network needs to figure out where it is, and then path-integrate with respect to a velocity signal. As the authors remark in Methods 4.2, without additional input, the only way to infer location is from border interactions. It is like navigating in absolute darkness. Therefore, it seems likely that the salient wall representations found in the recurrent units are just a consequence of the specific navigation task here; it is unclear if the same would apply in natural navigation. In natural navigation, there are many more sensory cues that help inferring location, most importantly vision, but also smell and whiskers/touch (which provides a more direct wall interaction; here, wall interactions are indirect by constraining velocity vectors). There is a similar but weaker concern about whether the (place cell like) localised firing fields of the output units are a direct consequence of the decoding procedure that only considers activity center locations.

The conclusion that 'representations are attractive' (heading of section 2) is not entirely supported. The authors show 'attractor-like behaviour' within a single context, but there could be alternative explanations for the recovery of stable ratemaps after noise injection. For example, the noise injection could scramble the network's currently inferred position, so that it would need to re-infer its position from boundary interactions along the trajectory. In that case the stabilisation would be driven by the input, not just internal attractor dynamics. Indeed, the useful control analysis in Appendix D suggests such a mechanism: without a velocity signal, only for small noise injections the network returns to a high correlation state. Correlated representations are recovered for larger noise injections due to the same mechanism that allow the network to determine its position upon from an uninformative initial hidden state upon entering a new environment, i.e. boundary interactions.

The authors report empirical data that shows clustering of place cell centers like they find for their output units. They report that 'there appears to be a tendency for the clusters to arrange in hexagonal fashion, similar to our computational findings'. This is an interesting observation on the distribution of place field centres which seems justified based on the example animal shown, but not across the population of animals included.

Reviewer #2 (Public review):

Summary:

The authors proposed a neural network model to explore the spatial representations of the hippocampal CA1 and entorhinal cortex (EC) and the remapping of these representations when multiple environments are learned. The model consists of a recurrent network and output units (a decoder) mimicking the EC and CA1, respectively. The major results of this study are: the EC network generates cells with their receptive fields tuned to a border of the arena; the decoder develops neuron clusters arranged in a hexagonal lattice. Thus, the model accounts for entrohinal border cells and CA1 place cells. It suggests that the remapping of place cells occurs between different environments through state transitions corresponding to unstable dynamical modes in the recurrent network.

Strengths:

The authors found a spatial arrangement of receptive fields similar to their model's prediction in experimental data recorded from CA1. Thus, the model proposes plausible mechanisms to generate hippocampal spatial representations without relying on grid cells. The model also suggests an interesting possibility that path integration is not the speciality of grid cells.

Weaknesses:

The role of grid cells in the proposed view, i.e., the boundary-to-place-to-grid model, remains elusive. The model can generate place cells without generating entorhinal grid cells. Moreover, the model can generate hexagonal grid patterns of place cells in a large arena. Whether and how the proposed model is integrated into the entire picture of the hippocampal-entorhinal meｍory processing remains elusive.

Reviewer #3 (Public review):

Summary:

The authors used recurrent neural network modelling of spatial navigation tasks to investigate border and place cell behaviour during remapping phenomena.

Strengths:

The neural network training seemed for the most part (see comments later) well-performed, and the analyses used to make the points were thorough.

The paper and ideas were well-explained.

Figure 4 contained some interesting and strong evidence for map-like generalisation as environmental geometry was warped.

Figure 7 was striking and potentially very interesting.

It was impressive that the RNN path-integration error stayed low for so long (Fig A1), given that normally networks that only work with dead-reckoning have errors that compound. I would have loved to know how the network was doing this, given that borders did not provide sensory input to the network. I could not think of many other plausible explanations... It would be even more impressive if it was preserved when the network was slightly noisy.

Update:

The analysis of how the RNN remapped, using a context signal to switch between largely independent maps, and the examination of the border like tuning in the recurrent units of the RNN, were both thorough and interesting. Further, in the updated response I appreciated the additional appendix E which helped substantiate the claim that the RNN neurons were border cells.

Weaknesses:

The remapping results were also puzzling. The authors present convincing evidence that the recurrent units effectively form 6 different maps of the 6 different environments (e.g. the sparsity of the code, or fig 6a), with the place cells remapping between environments. Yet, as the authors point out, in neural data the finding is that some cells generalise their co-firing patterns across environments (e.g. grid cells, border cells), while place cells remap, making it unclear what correspondence to make between the authors network and the brain. There are existing normative models that capture both entorhinal's consistent and hippocampus' less consistent neural remapping behaviour (Whittington et al. and probably others), what have we then learnt from this exercise?

Update: see summary below

I felt that the neural data analysis was unconvincing. Most notably, the statistical effect was found in only one of seven animals. Random noise is likely to pass statistical tests 1 in 20 times (at 0.05 p value), this seems like it could have been something similar? Further, the data was compared to a null model in which place cell fields were randomly distributed. The authors claim place cell fields have two properties that the random model doesn't (1) clustering to edges (as experimentally reported) and (2) much more provocatively, a hexagonal lattice arrangement. The test seems to collude the two; I think that nearby ball radii could be overrepresented, as in figure 7f, due to either effect. I would have liked to see a computation of the statistic for a null model in which place cells were random but with a bias towards to boundaries of the environment that matches the observed changing density, to distinguish these two hypotheses.

Update: the authors acknowledge these shortcomings and have appropriately tempered their data related claims.

Some smaller weaknesses:<br /> - Had the models trained to convergence? From the loss plot it seemed like not, and when including regularisors recent work (grokking phenomena, e.g. Nanda et al. 2023) has shown the importance of letting the regularisor minimise completely to see the resulting effect. Else you are interpreting representations that are likely still being learnt, a dangerous business.<br /> Update: I understand that practical limitations make testing this thoroughly impossible, which is fair enough.

- The claim that this work provided a mathematical formalism of the intuitive idea of a cognitive map seems strange, given that upwards of 10 of the works this paper cite also mathematically formalise a cognitive map into a similar integration loss for a neural network.<br /> Update: the introduction of these ideas hasn't changed, and my concerns above remain.

Aim Achieved? Impact/Utility/Context of Work

I think this is a thorough exploration of how this network with these losses is able to path-integrate its position and remap. This is useful, it is good to know how another neural network with slightly different constraints learns to perform these behaviours.

In the updated version of the manuscript I am happy to say that I think there are few claims that are unsubstantiated (see weakness section above that has been significantly updated). The link to neuroscience remains the biggest shortcoming of this work in my view. The authors point to two main results in this direction. First, the ability for interactions only between border-type and place cells to produce many observed place-cell results, providing a new hypothesis. Second, a connection between grid cells, place cells, and border cells, in the production of hexagonal arrangements of place cells.

Regarding the first, as the authors discuss, current evidence suggests border cells are invariant across environments whereas this work finds border cells for specific environments (they use the words rate-remapping boundary-type cells). It seems likely to me that there are many ways a neural network can path-integrate across different environments. In other models where the same base map is re-used (e.g. TEM) grid cells emerge, in this work where the maps for different environments are disjoint these border-like cells that do not match an observed cell type in their tuning to environment are involved. I find this a really interesting alternative (I think what an RNN does is interesting in its own right), but I don't see why I should think it is what the brain does, given that it appears to match observations less well (existence of grid cells, consistent firing patterns of border cells across environments). The smoking gun in favour of the author's hypothesis would be finding these sparse border like cells, or some other evidence of gating like interactions between border and place cells as they discuss. Finding such evidence sounds difficult (so not reasonable to ask for in a rebuttal), and to reiterate, I applaud the authors for clearly outlining an alternative, but I remain unconvinced.

Regarding the second point, while the grid-like placement of field centres was cool, and I applaud the authors for including real neural data comparisons, as the authors say, the data is preliminary, and further evidence would be required to fully substantiate this claim.

As such, in my mind it is an interesting alternative hypothesis. I look forward to seeing experimental predictions or comparisons that can tighten the link, substantiating the claim that what this particular RNN is doing reflects the algorithms at work in the brain.

Reviewer #1 (Public review):

Summary:

In this study, authors explored how galanin affects whole-brain activity in larval zebrafish using wide-field Ca2+ imaging, genetic modifications, and drugs that increase brain activity. The authors conclude that galanin has a sedative effect on the brain under normal conditions and during seizures, mainly through the galanin receptor 1a (galr1a). However, acute "stressors(?)" like pentylenetetrazole (PTZ) reduce galanin's effects, leading to increased brain activity and more seizures. Authors claim that galanin can reduce seizure severity while increasing seizure occurrence, speculated to occur through different receptor subtypes. This study confirms galanin's complex role in brain activity, supporting its potential impact on epilepsy.

Strengths:

The overall strength of the study lies primarily in its methodological approach using whole-brain Calcium imaging facilitated by the transparency of zebrafish larvae. Additionally, the use of transgenic zebrafish models is an advantage, as it enables genetic manipulations to investigate specific aspects of galanin signaling. This combination of advanced imaging and genetic tools allows for addressing galanin's role in regulating brain activity.

Weaknesses:

We have carefully reviewed the revised manuscript and the authors' responses. While the authors have attempted to address the points raised, I find that the revisions and rebuttals are insufficient and not entirely adequate. The authors seem not to have modified the manuscript in any way to take our comments into account.

In particular, many of the methodological and conceptual issues I initially raised remain unresolved. For example, the fundamental concern regarding the use of whole-brain calcium imaging - a method that may not effectively capture the localized and network-specific nature of seizure initiation and propagation - has not been adequately addressed. The authors acknowledge some limitations but do not sufficiently discuss how these affect the interpretation of their findings or propose mitigations. This could be added to the discussion section.

Additionally, the characterization of PTZ as a "stressor" remains problematic. Although the authors have retained this terminology, PTZ is widely understood to act primarily as a proconvulsant agent rather than a general stressor, and framing it otherwise continues to appear like a model-fitting rather than evidence-driven decision. The authors should consider changing the terminology throughout the manuscript and address these concerns when discussing their choice of PTZ as "stressor".

The discussion of the EAAT2 mutant model also remains incomplete. Although the authors mention preliminary transcriptome analyses, no new data were included, and it is stated that the evaluation is ongoing. Without thorough gene expression data, alternative explanations for the hypoactivity phenotype (such as changes in AMPA receptor or other critical neurotransmission-related genes) remain plausible and unaddressed. Moreover, the authors' acknowledgement that galanin upregulation is "at best one of a suite of regulatory mechanisms" further diminishes the centrality of their conclusions without sufficiently reworking the narrative of the study.

Finally, the finding that double knockout animals for EAAT2 and galanin showed little difference in seizure susceptibility compared to EAAT2 knockouts alone suggests that galanin upregulation may not play a dominant functional role, yet this important implication is not adequately reflected in the interpretation of the results.

Conclusion:

In summary, although the authors have made some efforts to respond to the critiques, I do not believe the manuscript has been substantially improved in response in R2, and I do not see reason to change my original assessment made after R1. The major conceptual and methodological concerns remain largely unaddressed, limiting the impact and validity of the study's conclusions. These concerns should be addressed not only in the rebuttal letter but also in the manuscript.

Reviewer #2 (Public review):

This revised paper describes an investigation of galanin and galanin receptor signaling on whole-brain activity in the context of recurrent seizure activity or under homeostatic basal conditions. The authors primarily use calcium imaging to observe whole-brain neuronal activity accompanied by galanin qPCR to determine how manipulations of galanin or the galr1a receptor affect the activity of the whole-brain under non-ictal conditions or when seizure activity occurs. The authors use their eaat2a-/- model (introduced in their Glia 2022 paper, PMID 34716961) that shows recurrent seizure activity as well as suppression of neuronal activity and locomotion interictally. It is compared to the well-known pentylenetetrazole (PTZ) pharmacological model of seizures in zebrafish. Given the literature cited in their Introduction, the authors hypothesize that galanin will exert a net inhibitory effect on brain activity in models of seizures/epilepsy. They were surprised to find that this hypothesis was only moderately supported in their eaat2a-/- model. In contrast, after PTZ, fish with galanin overexpression showed increased seizure number and reduced duration while fish with galanin KO showed reduced seizure number and increased duration.

Previous concerns about sex or developmental biological variables were addressed, as their model's seizure phenotype emerges rapidly and long prior to the establishment of zebrafish sexual maturity. However, it remains unclear whether all seizures detected via calcium imaging alone are also seizures that are detectable at the level of animal behavior. To confirm this, a validation of the threshold used for calcium imaging of "neuronal seizures" would be required to determine if this threshold detects only "neuronal seizures" that co-occur with behavioral seizures. Overall, this study is important and convincing, and carries clear value for understanding the multifaceted functions that neuronal galanin can perform under homeostatic and disease conditions.

Additional Concerns:

- The authors have validated their ability to measure behavioral seizures quantitatively in their 2022 Glia paper but the information provided on defining behavioral seizures as they map onto seizures detected via imaging alone was limited. The definition of behavioral seizure activity as it relates to calcium fluctuations is not expanded upon in this paper, but could provide detail about how the behavioral seizures relate to a seizure detected via calcium imaging alone.

Reviewer #3 (Public review):

Summary:

The neuropeptide galanin is primarily expressed in the hypothalamus and has been shown to play critical roles in homeostatic functions such as arousal, sleep, stress, and brain disorders such as epilepsy. Previous work in rodents using galanin analogs and receptor-specific knockout have provided convincing evidence for anti-convulsant effects of galanin.

In the present study, the authors sought to determine the relationship between galanin expression and whole-brain activity. The authors took advantage of the transparent nature of larval zebrafish to perform whole-brain neural activity measurements via widefield calcium imaging. Two models of seizures were used (eaat2a-/- and pentylenetetrazol; PTZ). In the eaat2a-/- model, spontaneous seizures occur and the authors found that galanin transcript levels were significantly increased and associated with reduced frequency of calcium events. Similarly, two hours after PTZ galanin transcript levels roughly doubled and the frequency and amplitude of calcium events were reduced, while the duration increased.

The authors also used a heat shock protein line (hsp70I:gal) where galanin transcripts levels are induced by activation of heat shock protein, but this line also shows higher basal transcript levels of galanin. Due to problems with whole-brain activity in wild-type larvae, the authors used the line without heat shock. They found higher level of galanin in hsp70I:gal larval zebrafish resulted in a reduction in the number of calcium events and amplitude. In contrast, galanin knockout (gal-/-) significantly increased calcium activity, indicated by an increased number of calcium events, but a reduction in amplitude and duration. Antibody staining confirmed the absence of galanin expression in gal-/- knockouts. Knockout of the galanin receptor subtype galr1a via crispants also increased the frequency of calcium events without influencing amplitude or duration.

In subsequent experiments in eaat2a-/- mutants were crossed with hsp70I:gal or gal-/- to modify galanin expression. These experiments showed modest effects, with eaat2a-/- x gal-/- knockouts showing an increased normalized area under the curve and seizure amplitude.

Lastly, the authors attempted to study the relationship between galanin and brain activity during a PTZ challenge. The hsp70I:gal larva showed increased number of seizures and reduced seizure duration during PTZ. In contrast, gal-/- mutants showed increased normalized area under the curve and a stark reduction in number of detected seizures, a reduction in seizure amplitude, but an increase in seizure duration. The authors then ruled out the role galanin a1 receptor in modulating this effect during PTZ, since the number of seizures was unaffected, whereas the amplitude and duration of seizures was increased in galr1a knockouts.

Strengths:

(1) The gain- and loss-of function galanin manipulations provided convincing evidence that galanin influences brain activity (via calcium imaging) during interictal and/or seizure-free periods. The relationship between galanin transcript levels and brain activity in figures 1 & 2 was convincing. Antibody staining also supports the absence of galanin in gal-/- mutants. Moreover, galanin transcript levels were unchanged in galr1ako brains, suggesting the lack of compensatory effects.

(2) The authors use two models of epilepsy (eaat2a-/- and PTZ).

(3) Supplementary video files for calcium imaging support the observations.

Weaknesses:

(1) I disagree with the idea that PTZ is a 'stressor'. This was raised in previous reviews and has not been acknowledged sufficiently.

(2) Although the relationship between galanin and brain activity during interictal or seizure-free periods was clear, the mechanisms that influence excitability during PTZ remain unclear. The authors show that galr1a does not mediate this effect, since seizure amplitude and duration were more severe in galr1a KO. Therefore, it remains unclear which galanin receptor is modulating this inhibitory effect.

(3) The manuscript is heavily reliant on calcium imaging for interpretation.<br /> Additional methods could strengthen the data, translational relevance, and interpretation (e.g., acute pharmacology using selective galanin agonists or antagonists, brain or cell recordings, biochemistry, etc).

Reviewer #1 (Public review):

Summary:

In this manuscript, Ru and colleagues investigated regulatory gene interactions during osteogenic differentiation. By profiling transcriptomic changes during mesenchymal stem cell differentiation, they identified KLF16 as a key transcription factor that inhibits osteogenic differentiation and mineralization. It was found that overexpression of KLF16 suppressed osteogenesis in vitro, while KLF16⁺/⁻ mice exhibited enhanced bone density, underscoring its regulatory role in bone formation.

Strengths:

(1) Bioinformatics is strong and comprehensive.

(2) Identification of KLF16 in osteoblast differentiation is exciting and innovative.

Weaknesses:

(1) The mechanism of KLF16 function is not studied.

(2) Studies of KLF16 in bone development, from both in vitro and in vivo perspectives, are descriptive.

(3) Findings in bioinformatics analysis are mostly redundant with previous studies in the field, and can be simplified.

Reviewer #2 (Public review):

In their manuscript with the title "Integrated transcriptomic analysis of human induced pluripotent stem cell (iPSC)-derived osteogenic differentiation reveals a regulatory role of KLF16", Ru et al. have analyzed the gene expression changes during the osteogenic differentiation of iPSC-derived mesenchymal stem/stromal cells into preosteoblasts and osteoblasts. As part of the computational analyses, they have investigated the transcription factor regulatory network mediating this differentiation process, which has also led to the identification of the transcription factor KLF16. Overexpression experiments in vitro and the analysis of heterozygous KLF16 knockout mice in vivo indicate that KLF16 is an inhibitor of osteogenic differentiation.

The integrated analysis of iPSC bulk transcriptomic data is a major strength of the study, and it is also great that the authors provide deeper functional characterization of the transcription factor KLF16, one of the newly identified candidate regulators of osteogenic differentiation.

However, characterization of KLF16 expression in the mouse and validation of the knockout model are currently lacking. Alternative explanations for the mutant phenotype should be considered to improve the strength of the conclusions.

If all issues can be addressed, the study would provide an important resource for the field that would facilitate future research on the regulation of osteogenesis in vitro and in vivo, with potential implications for preclinical and clinical research as well as bioengineering.

Reviewer #1 (Public review):

The authors have developed a contextual fear learning (CFC) paradigm in head-fixed mice that produces freezing as the conditioned response. Typically, lick suppression is the conditioned response in such designs, but this 1) introduces a potential confounding influence of reward learning on neural assessments of aversion learning and 2) does not easily allow comparison of head-fixed studies with extensive previous work in freely moving animals, which use freezing as the primary conditioned response. This report describes 3 versions of this virtual reality CFC paradigm, its validation using place-cell remapping, and provides suggestions for further refinement and application.

The first part of this study is a report on the development and outcomes of 3 variations of the CFC paradigm in a virtual reality environment. The fundamental design is strong, with head-fixed mice required to run down a linear virtual track to obtain a water reward. Once trained, the water reward is no longer necessary and mice will navigate virtual reality environments. There are rigorous performance criteria to ensure that mice that make it to the experimental stage show very low levels of inactivity prior to fear conditioning. These criteria do result in only 40% of the mice making it to the experimental stage, but high rates of activity in the VR environment is crucial for detecting learning-related freezing. It is possible that further adjustments to the procedure could improve attrition rates.

Paradigm versions 1 and 2 vary the familiarity of the control context while paradigm versions 2 and 3 vary the inter-shock interval. Version 1 is the most promising, showing the greatest increase in conditioned freezing (~40%) and good discrimination between contexts (delta ~15-20%). Version 2 showed no clear evidence of learning - average freezing at recall day 1 was not different than pre-shock freezing. First lap freezing showed a difference, but this single lap effect is not useful for many of the neural circuit questions for which this paradigm is meant to facilitate. Version 3 produces greater freezing and slower extinction than version 2. While the magnitude of the context discrimination is less than that in version 1, further optimization of the VR CFC is likely to produce robust learning and extinction. The authors discuss several options for further optimization.

The second part of the study is a validation of the head-fixed CFC VR protocol through demonstration that fear conditioning leads to remapping of dorsal CA1 place fields, similar to that observed in freely moving subjects. The results support this aim and largely replicate previous findings in freely moving subjects. One difference from previous work of note is that VR CFC led to remapping of the control environment, not just the conditioning context. The authors present several possible explanations for this lack of specificity to the shock context. While this experiment examined place cell remapping after fear conditioning, it did not attempt to link neural activity to the learned association or freezing behavior.

In summary, this is an important methodological innovation and this study sets the initial parameters and neuronal validation needed to further optimize a head-fixed CFC paradigm that produces freezing. In the discussion, the authors note the limitations of this study, suggest next steps in refinement, and point to several future directions using this protocol to significantly advance our understanding of the neural circuits of threat-related learning and behavior.

Comments on revisions:

The manuscript is much stronger with the additions and revisions the authors provided in their revised submission.

Reviewer #2 (Public review):

Summary:

In this manuscript, Krishnan et al devised three paradigms to perform contextual fear conditioning in head-fixed mice. Each of the paradigms relied on head-fixed mice running on a treadmill through virtual reality arenas. The authors tested the validity of three versions of the paradigms by using various parameters. The authors have addressed some of my initial concerns in their revised manuscript.

Strengths:

The authors have devised three new contextual fear conditioning paradigms in head-fixed mice. The authors tested a number of parameters towards optimization of this approach.

Weaknesses:

While some experimental parameters were tested in the manuscript, it appears that a large amount of additional testing and optimization will be required before reliable behavioral responses can be acquired and ultimately for the paradigm(s) to be useful for answering biological questions. One major factor will be optimizing parameters such that head-fixed mice in this paradigm can (largely) recapitulate what is observed in freely behaving mice. This may be challenging however, as they have previously published one of the three paradigms and the extensive additional testing they did in this current manuscript did not greatly improve the experimental setup. This may indicate limited immediate usefulness for the community as significant work likely remains for optimization.

Achievement of Aims:

The authors have put a significant amount of work in testing the paradigms, and as a result, progress has been made towards their usefulness in the field. However, a significant amount of optimization likely exists.

Impact on the field:

The development of a reliable paradigm for studying contextual fear in head-fixed animals would be a strong contribution to the field as it would enable sophisticated cell and circuit imaging analyses. This study is a good start towards this goal, but significant optimization is required for the paradigm(s) to fully benefit the field - especially to allow those who may have less experience in these approaches to use it in their own research.

Reviewer #3 (Public review):

Summary:

Krishnan et al. present a novel contextual fear conditioning (CFC) paradigm using a virtual reality (VR) apparatus to evaluate whether conditioned context-induced freezing can be elicited in head-fixed mice. By combining this approach with two-photon imaging, the authors aim to provide high-resolution insights into the neural mechanisms underlying learning, memory, and fear. Their experiments demonstrate that head-fixed mice can discriminate between threat and non-threat contexts, exhibit fear-related behavior in VR, and show context-dependent variability during extinction. Supplemental analyses further explore alternative behaviors and the influence of experimental parameters, while hippocampal neuron remapping is tracked throughout the experiments, showcasing the paradigm's potential for studying memory formation and extinction processes.

Strengths:

Methodological Innovation: The integration of a VR-based CFC paradigm with real-time two-photon imaging offers a powerful, high-resolution tool for investigating the neural circuits underlying fear, learning, and memory.

Versatility and Utility: The paradigm provides a controlled and reproducible environment for studying contextual fear learning, addressing challenges associated with freely moving paradigms.

Potential for Broader Applications: By demonstrating hippocampal neuron remapping during fear learning and extinction, the study highlights the paradigm's utility for exploring memory dynamics, providing a strong foundation for future studies in behavioral neuroscience.

Comprehensive Data Presentation: The inclusion of supplemental figures and behavioral analyses (e.g., licking behaviors and variability in extinction) strengthens the manuscript by addressing additional dimensions of the experimental outcomes.

Weaknesses:

Optimization: many parameters remain to be tested in the VR fear conditioning paradigm.

Extended training and attrition rate: the paradigm requires weeks of training and only 40% of mice reach criteria.

Reviewer #1 (Public review):

Summary:

In this manuscript, Guo and colleagues used a cell rounding assay to screen a library of compounds for inhibition of TcdB, an important toxin produced by Clostridioides difficile. Caffeic acid and derivatives were identified as promising leads, and caffeic acid phenethyl ester (CAPE) was further investigated.

Strengths:

Considering the high morbidity rate associated with C. difficile infections (CDI), this manuscript presents valuable research in the investigation of novel therapeutics to combat this pressing issue. Given the rising antibiotic resistance in CDI, the significance of this work is particularly noteworthy. The authors employed a robust set of methods and confirmatory tests, which strengthen the validity of the findings. The explanations provided are clear, and the scientific rationale behind the results is well-articulated. The manuscript is extremely well written and organized. There is a clear flow in the description of the experiments performed. Also, the authors have investigated the effects of CAPE on TcdB in careful detail, and reported compelling evidence that this is a meaningful and potentially useful metabolite for further studies.

Weaknesses:

Although the authors have made changes to the manuscript to address some of my comments, many of the comments were not satisfactorily addressed. Many of the changes are still superficial, and some concerns still need to be addressed. Important details are still missing from the description of some experiments. Authors should carefully revise the manuscript to ascertain that all details that could affect interpretation of their results are presented clearly.

There is still very little discussion (none, really) in the manuscript about the fact that, because the authors observed a significant effect of CAPE on both bacterial growth and spore production, some of the phenotypes observed can no longer be attributed solely to toxin inhibition.

The details about mass spectrometry are still insufficient. It is still unclear whether metabolite identifications were always based on MS1 or MS2. Instead, several details that are really secondary were included. Authors should be unequivocally clear as to how metabolite identities were obtained. They should also indicate which mass spectrometer was used, and there should be a section in the Materials and Methods describing these experiments.

About the removal of carry-over compounds, the authors stated that ultrafiltration centrifugal partition was used. However, although the authors explained this in detail in their response to reviewers file, the details were omitted from the main text. Authors should clearly state in the manuscript text that "Due to the large molecular weight of TcdB, approximately 270 kDa, we selected a 100 kDa molecular weight cutoff ultrafiltration membrane. The centrifugation was performed at 4000 g for 5 min to eliminate the compounds that did not bind to TcdB."

These are important details which need to be included.

Reviewer #2 (Public review):

I appreciate the author's responses to my original review. This is a comprehensive analysis of CAPE on C. difficile activity. It seems like this compound effects all aspects of C. difficile, which could make it effective during infection but also make it difficult to understand the mechanism. Even considering the authors responses, I think it is critical for the authors to work on the conclusions regarding the infection model. There is some protection from disease by CAPE but some parameters are not substantially changed. For instance, weight loss is not significantly different in the C. difficile only group versus the C. difficile + CAPE group. Histology analysis still shows a substantial amount of pathology in the C. difficile + CAPE group. This should be discussed more thoroughly using precise language.

The authors did a good job addressing my concerns regarding the infection model by providing a more accurate descriptions in the Results section for histology. However, the weight loss improvement by CAPE does not look like a significant effect, although it is trending towards improvement. This should be more accurately described.

Another minor concern is that the current Abstract is overstating the amount of disease attenuation. I would replace "remarkably reduces the pathology" with "reduces some of the pathology"

Oh man this website is so freaking useful. dr.loudness-war.info

Gives dynamic range information on albums/tracks.

This guy has Bandcamp as favorite site for buying music.

As second, he has Qobuz.

He also uses 7digital & HDtracks

Reviewer #1 (Public review):

Summary:

Brdar, Osterburg, Munick, et al. present an interesting cellular and biochemical investigation of different p53 isoforms. The authors investigate the impact of different isoforms on the in-vivo transcriptional activity, protein stability, induction of the stress response, and hetero-oligomerization with WT p53. The results are logically presented and clearly explained. Indeed, the large volume of data on different p53 isoforms will provide a rich resource for researchers in the field to begin to understand the biochemical effects of different truncations or sequence alterations.

Strengths:

The authors achieved their aims to better understand the impact/activity of different p53 is-forms, and their data well support their statements. Indeed, the major strengths of the paper lie in its comprehensive characterization of different p53 isoforms and the different assays that are measured. Notably, this includes p53 transcriptional activity, protein degradation, induction of the chaperone machinery, and hetero-oligomerization with wtp53. This will provide a valuable dataset where p53 researchers can evaluate the biological impact of different isoforms in different cell lines. The authors went to great lengths to control and test for the effect of (1) p53 expression level, (2) promotor type, and (3) cell type. I applaud their careful experiments in this regard.

Comments on revised version:

The authors have addressed all of my concerns convincingly, including with a new mass spectrometry experiment to quantify p53 peptides specifically.

Reviewer #1 (Public review):

Batra, Cabrera and Spence et al. present a model which integrates histone posttranslational modification (PTM) data across cell models to predict gene expression with the goal of using this model to better understand epigenetic editing. This gene expression prediction model approach is useful if a) it predicts gene expression in specific cell lines b) it predicts expression values rather than a rank or bin, c) if it helps us to better understand the biology of gene expression or d) it helps us to understand epigenome editing activity. Problematically for points a) and b) it is easier to directly measure gene expression than to measure multiple PTMs and so the real usefulness of this approach mostly relates to c) and d).

Other approaches have been published that use histone PTM to predict expression (e.g. PMID 27587684, 36588793). Is this model better in some way? No comparisons are made, although a claim is made that direct comparisons are difficult. I appreciate that the authors have not used the histone PTM data to predict gene expression levels of an "average cell" but rather that they are predicting expression within specific cell types or for unseen cell types. Approaches that predict expression levels are much more useful, whereas some previous approaches have only predicted expressed or not expressed or a rank order or bin-based ranking. The paper does not seem to have substantial novel insights into understanding the biology of gene expression.

The approach of using this model to predict epigenetic editor activity on transcription is interesting and to my knowledge novel although only examined in the context of a p300 editor. As the author point out the interpretation of the epigenetic editing data is convoluted by things like sgRNA activity scoring and to fully understand the results likely would require histone PTM profiling and maybe dCas9 ChIP-seq for each sgRNA which would be a substantial amount of work.

Furthermore from the model evaluation of H3K9me3 is seems the model is performing modestly for other forms of epigenetic or transcriptional editing- e.g. we know for the best studied transcriptional editor which is CRISPRi (dCas9-KRAB) that recruitment to a locus is associated with robust gene repression across the genome and is associated with H3K9me3 deposition by recruitment of KAP1/HP1/SETDB1 (PMID: 35688146, 31980609, 27980086, 26501517).

One concern overall with this approach is that dCas9-p300 has been observed to induce sgRNA independent off target H3K27Ac (https://www.ncbi.nlm.nih.gov/pmc/articles/PMC8349887/ see Figure S5D) which could convolute interpretation of this type of experiment for the model.

Comments on revisions: This resubmission adds a comparison to existing gene prediction methods, but add no new confirmation experiments with predicting epigenome editing efficiency and had only one minor text edit.

Reviewer #2 (Public review):

Summary:

The authors build a gene expression model based on histone post-translational modifications, and find that H3K27ac is correlated with gene expression. They compare to other gene prediction methods such as DeepChrome. They proceed to perturb H3K27ac at 13 gene promoters in two cell types, and measure gene expression changes to test their model.

Strengths:

The combination of multiple methods to model expression, along with utilizing 6 histone datasets in 13 cell types allowed the authors to build a model that correlates between 0.7-0.79 with gene expression.<br /> They compare three cells types to other prediction models, and this figure should be included in the main figures.<br /> They use dCas9-p300 fusions to perturb H3K27ac and monitor gene expression to test their model. Ranked correlations of the HEK293 data showed some support for the predictions after perturbation of H3K27ac.

Weaknesses:

The authors state in the latest submission that the primary use case of this work is related to predicting epigenome editing outcomes, not predicting gene expression from chromatin. However the first four figures all relate to gene expression prediction. The only main figure that shows epigenome editing prediction is panel 6E. If this authors wish to highlight the use case of this work they should redo figures, including moving panels from current supplemental figures to show this.

The perturbation of 5 genes in K562 with perturb-seq data shows a modest correlation of ~0.5 and is still only shown in supplemental figures, which is odd as this is the true test case of their model in my opinion. The authors are then left to speculate the reasons why the outcome of epigenome editing doesn't fit their predictions, which highlights the limited value in the current version of this method.<br /> As mentioned before, testing genes that were not expressed being most activated by dCas9-p300 weaken the correlations vs. looking at a broad range of different gene expression as the original model was trained on.

If the authors want this method to be used to predict outcomes of epigenome editing, expanding to dCas9-KRAB and other CRISPRa methods (SAM and VPR) would be useful. Those datasets are published and could be analyzed for this manuscript and show how the model holds up across cell types and epigenome editing methods.

The utility of this method as described here, to predict gRNA outcomes seems modest and limited. It is fairly trivial to test 10 or more gRNAs for a single gene to find the best one, and the authors show limited prediction and occasionally no benefit. For example, with CHD8 and CD79 the gRNA with the highest prediction had the lowest actual impact on gene expression of the gRNAs tested. For many other genes the gRNA's prediction and gene expression outcome show no correlation.

Reviewer #1 (Public review):

Summary:

In the present study, the authors examined the possibility of using phosphatidyl-inositol kinase 3-kinase alpha (PI3Ka) inhibitors for heterotopic ossification in fibrodysplasia ossificans progressiva. Administration of BYL719, a chemical inhibitor of PI3Ka, prevented heterotopic ossification in a mouse model of FOP that expressed a mutated ACVR1 receptor. Genetic ablation of PI3Ka also suppressed heterotopic ossification in mice. BYL719 blocked osteo/chondroprogenitor specification and reduced inflammatory responses by reducing the number of fibro-adipogenic progenitors (FAPs) and promoting muscle fibre regeneration in vivo. The authors claimed that inhibition of PI3Ka is a safe and effective therapeutic strategy for heterotopic ossification.

Strengths:

Taking together previous reports on the specificity of BY718 in PI3K, it was suggested that BYL719 inhibits heterotopic ossification by reducing FAPs and promoting muscle regeneration through the PI3K pathway in vivo.

Weaknesses:

In the original manuscript, there was the possibility that BYL719 inhibited heterotopic ossification through non-specific and toxic effects rather than the PI3k pathway.

However, the authors added new data and explanations in the revision to solve the possibility. The findings of the authors would be useful and would provide an additional direction to develop a therapeutic strategy for heterotopic ossification.

Reviewer #2 (Public review):

Summary:

Authors in this study previously reported that BYL719, an inhibitor of PI3Kα, suppressed heterotopic ossification in mice model of a human genetic disease, fibrodysplasia ossificans progressive, which is caused by the activation of mutant ACVR1/R206H by Activin A. The aim of this study is to identify the mechanism of BYL719 for the inhibition of heterotopic ossification. They found that BYL719 suppressed heterotopic ossification in two ways: one is to inhibit the specification of precursor cells for chondrogenic and osteogenic differentiation and the other is to suppress the activation of inflammatory cells.

Strengths:

This study is based on authors' previous reports and the experimental procedures including the animal model are established. In addition, to confirm the role of PI3Kα, authors used the conditional knock-out mice of the subunit of PI3Kα. They clearly demonstrated the evidence indicating that the targets of PI3Kα is not members of TGFBR by a newly established experimental method.

Weaknesses:

Overall, the presented data were closely related to those previously published by authors' group or others and there were very few new findings. The molecular mechanisms through which BYL719 inhibits HO remain unclear, even in the revised manuscript.<br /> Heterotopic ossification in mice model was not stable and inappropriate for the scientific evaluation.<br /> The method for chondrogenic differentiation was not appropriate, and the scientific evidence of successful differentiation was lacking.<br /> The design of gene expression profile comparison was not appropriate and failed to obtain the data for the main aim of this study.<br /> The experiments of inflammatory cells were performed cell lines without ACVR1/R206H mutation, and therefore the obtained data were not precisely related to the inflammation in FOP.

Comments on revisions:

In the R2 version, the authors performed additional experiments using mice with inducible human R206H ACVR1A. BM-MSCs isolated from these mice were used to investigate the effect of Activin-A. The results again suggested that BYL79 inhibited the chondrogenic differentiation of BM-MSCs. However, there are still no data demonstrating the effect of BYL79 on cell growth in these in vitro experiments. In Figures 7A-D, 10 μM BYL79 strongly inhibited the proliferation of inflammatory cells, suggesting that growth inhibition may have contributed to the results shown in Figure 5.

The main point of discussion concerns the significance of the comparisons made. The fundamental disagreement arises from the role of Activin-A in R206H cells and its effect on chondrogenic differentiation. The authors' rebuttal regarding my comments on the RNA-seq analyses should be reconsidered. The core issue lies in the interpretation of Activin-A's role in R206H cells and the distinction between chondrogenic differentiation and ossification.

A key feature of R206H mutant cells is that they respond to Activin-A by activating Smad1/5 signaling-comparable in quality to the signaling induced by BMP6 in WT cells. Another important point, as also acknowledged by the authors, is that Activin-A can transduce Smad2/3 signaling via its canonical receptor, ACVR1B. These dual signaling pathways synergistically contribute to chondrogenic differentiation in precursor cells such as FAPs. Several reports have demonstrated that the combined activation of TGF-β and BMP signaling promotes chondrogenesis more strongly than either pathway alone.

Since the PI3Kα inhibition effect on HO is already known, a critical question in this study is whether BYL79 also inhibits the Smad2/3 pathway. A straightforward experiment would be to compare WT cells treated with Activin-A alone versus Activin-A plus BYL79, and to perform GO term enrichment analyses related specifically to chondrogenic differentiation, not ossification. Additionally, comparing R206H cells treated with Activin-A/BYL79 and WT cells treated with BMP6/BYL79 could help identify gene sets inhibited by BYL79 via Smad2/3 signaling. If these comparisons reveal no specific effect on genes related to chondrogenesis, the effect of BYL79 may be limited to suppression of BMP-mediated osteogenesis. Unfortunately, the authors appear to show little interest in addressing this issue.

Regarding Figure 7, the authors' rebuttal should also be reconsidered. Since the R2 version employed FOP model mice, it would have been possible to evaluate the effects of BYL79 on inflammatory cells harboring the R206H mutation. This could have enabled a more precise assessment of BYL79's influence on inflammatory signaling. While the authors repeatedly claim that BYL79's effect is not specific to any particular ligand or the presence of the FOP mutation, the role of TGF-β signaling in the development of endochondral heterotopic ossification is well recognized. Therefore, the mechanism of BYL79 should be clarified before considering its therapeutic application

Reviewer #1 (Public review):

Summary:

In a previous work Prut and colleagues had shown that during reaching, high frequency stimulation of the cerebellar outputs resulted in reduced reach velocity. Moreover, they showed that the stimulation produced reaches that deviated from a straight line, with the shoulder and elbow movements becoming less coordinated. In this report they extend their previous work by addition of modeling results that investigate the relationship between the kinematic changes and torques produced at the joints. The results show that the slowing is not due to reductions in interaction torques alone, as the reductions in velocity occur even for movements that are single joint. More interestingly, the experiment revealed evidence for decomposition of the reaching movement, as well as an increase in the variance of the trajectory.

Strengths:

This is a rare experiment in a non-human primate that assessed the importance of cerebellar input to the motor cortex during reaching.

Weaknesses:

None

Reviewer #2 (Public review):

This manuscript asks an interesting and important question: what part of 'cerebellar' motor dysfunction is an acute control problem vs a compensatory strategy to the acute control issue? The authors use a cerebellar 'blockade' protocol, consisting of high frequency stimuli applied to the cerebellar peduncle which is thought to interfere with outflow signals. This protocol was applied in monkeys performing center out reaching movements and has been published from this laboratory in several preceding studies. I found the take-home-message broadly convincing and clarifying - that cerebellar block reduces muscle activation acutely particularly in movements that involve multiple joints and therefore invoke interaction torques, and that movements progressively slow down to in effect 'compensate' for these acute tone deficits. The manuscript was generally well written, data were clear, convincing and novel. The key strengths are differentiating acute from sub-acute (within session but not immediate) kinematic consequences of cerebellar block.

Reviewer #3 (Public review):

Summary:

In their revised manuscript, Sinha and colleagues aim to identify distinct causes of motor impairments seen when perturbing cerebellar circuits. This goal is an important one, given the diversity of movement related phenotypes in patients with cerebellar lesion or injury, which are especially difficult to dissect given the chronic nature of the circuit damage. To address this goal, the authors use high-frequency stimulation (HFS) of the superior cerebellar peduncle in monkeys performing reaching movements. HFS provides an attractive approach for transiently disrupting cerebellar function previously published by this group. First, they find a reduction in hand velocities during reaching, which was more pronounced for outward versus inward movements. By modeling inverse dynamics, they find evidence that shoulder muscle torques are especially affected. Next, the authors examine the temporal evolution of movement phenotypes over successive blocks of HFS trials. Using this analysis, they find that in addition to the acute, specific effects on torques in early HFS trials, there was an additional progressive reduction in velocity during later trials, which they interpret as an adaptive response to the inability to effectively compensate for interaction torques during cerebellar block. Finally, the authors examine movement decomposition and trajectory, finding that even when low velocity reaches are matched to controls, HFS produces abnormally decomposed movements and higher than expected variability in trajectory.

Strengths:

Overall, this work provides important insight into how perturbation of cerebellar circuits can elicit diverse effects on movement across multiple timescales.

The HFS approach provides temporal resolution and enables analysis that would be hard to perform in the context of chronic lesions or slow pharmacological interventions. Thus, this study describes an important advance over prior methods of circuit disruption in the monkey, and their approach can be used as a framework for future studies that delve deeper into how additional aspects of sensorimotor control are disrupted (e.g., response to limb perturbations).

In addition, the authors use well-designed behavioral approaches and analysis methods to distinguish immediate from longer-term adaptive effects of HFS on behavior. Moreover, inverse dynamics modeling provides important insight into how movements with different kinematics and muscle dynamics might be differentially disrupted by cerebellar perturbation.

In this revised version of the manuscript, the authors have provided additional analyses and clarification that address several of the comments from the original submission.

Remaining comments:

The argument that there are acute and adaptive effects to perturbing cerebellar circuits is compelling, but there seems to be a lost opportunity to leverage the fast and reversible nature of the perturbations to further test this idea and strengthen the interpretation. Specifically, the authors could have bolstered this argument by looking at the effects of terminating HFS - one might hypothesize that the acute impacts on joint torques would quickly return to baseline in the absence of HFS, whereas the longer-term adaptive component would persist in the form of aftereffects during the 'washout' period. As is, the reversible nature of the perturbation seems underutilized in testing the authors' ideas. While this experimental design was not implemented here, it seems like a good opportunity for future work using these approaches.

The analysis showing that there is a gradual reduction in velocity during what the authors call an adaptive phase is convincing. While it is still not entirely clear why disruption of movement during the adaptive phase is not seen for inward targets, despite the fact that many of the inward movements also exhibit large interaction torques, the authors do raise potential explanations in the Discussion.

The text in the Introduction and in the prior work developing the HFS approach overstates the selectivity of the perturbations. First, there is an emphasis on signals transmitted to the neocortex. As the authors state several times in the Discussion, there are many subcortical targets of the cerebellar nuclei as well, and thus it is difficult to disentangle target-specific behavioral effects using this approach. Second, the superior cerebellar peduncle contains both cerebellar outputs and inputs (e.g., spinocerebellar). Therefore, the selectivity in perturbing cerebellar output feels overstated. Readers would benefit from a more agnostic claim that HFS affects cerebellar communication with the rest of the nervous system, which would not affect the major findings of the study. In the revised manuscript, the authors do provide additional anatomical and evolutionary context and discuss potential limitations in the selectivity of HFS in the Materials and Methods. However, I feel that at least a brief mention of these caveats in the Introduction, where it is stated, "we then reversibly blocked cerebellar output to the motor cortex", would benefit the reader.

Reviewer #2 (Public review):

Summary:

The study characterized the dependence of spike timing-dependent long-term depression (tLTD) on presynaptic NMDA receptors and the intracellular cascade after NMDAR activation possibly involved in the observed decrease in glutamate probability release at L5-L5 synapses of the visual cortex in mouse brain slices.

Strengths:

The genetic and electrophysiological experiments are thorough. The experiments are well reported and mainly support the conclusions. This study confirms and extends current knowledge by elucidating additional plasticity mechanisms at cortical synapses, complementing existing literature.

Weaknesses:

No direct testing for ions passing trough standard NMDAR, mainly sodium and calcium is shown.

Reviewer #3 (Public review):

Summary:

In this manuscript, "Neocortical Layer-5 tLTD Relies on Non-Ionotropic Presynaptic NMDA Receptor Signaling", Thomazeau et al. seek to determine the role of presynaptic NMDA receptors and the mechanism by which they mediate expression of frequency-independent timing-dependent long-term depression (tLTD) between layer-5 (L5) pyramidal cells (PCs) in the developing mouse visual cortex. By utilizing sophisticated methods, including sparse Cre-dependent deletion of GluN1 subunit via neonatal iCre-encoding viral injection, in vitro quadruple patch clamp recordings, and pharmacological interventions, the authors elegantly show that L5 PC->PC tLTD is 1) dependent on presynaptic NMDA receptors, 2) mediated by non-ionotropic NMDA receptor signaling, and 3) is reliant on JNK2/Syntaxin-1a (STX1a) interaction (but not RIM1αβ) in the presynaptic neuron. The study elegantly and pointedly addresses a long-standing conundrum regarding the lack of frequency dependence of tLTD.

Strengths:

The authors did a commendable job presenting a very polished piece of work with high-quality data that this Reviewer feels enthusiastic about. The manuscript has several notable strengths. Firstly, the methodological approach used in the study is highly sophisticated and technically challenging, and successfully produced high-quality data that were easily accessible to a broader audience. Secondly, the pharmacological interventions used in the study targeted specific players and their mechanistic roles, unveiling the mechanism in question step-by-step. Lastly, the manuscript is written in a well-organized manner that is easy to follow. Overall, the study provides a series of compelling evidence that leads to a clear illustration of mechanistic understanding.

Weakness:

No major weaknesses were noted.

Reviewer #1 (Public review):

Summary:

The paper by Graff et al. investigates the function of foxf2 in zebrafish to understand the progression of cerebral small vessel disease. The authors use a partial loss of foxf2 (zebrafish possess two foxf2 genes, foxf2a and foxf2b, and the authors mainly analyze homozygous mutants in foxf2a) to investigate the role of foxf2 signaling in regulating pericyte biology. They find that the number of pericytes is reduced in foxf2a mutants and that the remaining pericytes display alterations in their morphologies. The authors further find that mutant animals can develop to adulthood, but that in adult animals, both endothelial and pericyte morphologies are affected. They also show that mutant pericytes can partially repopulate the brain after genetic ablation.

Strengths:

The paper is well written and easy to follow.

Weaknesses:

The results are mainly descriptive, and it is not clear how they will advance the field at their current state, given that a publication on mice has already examined the loss of foxf2 phenotype on pericyte biology (Reyahi, 2015, Dev. Cell).

(1) Reyahi et al. showed that loss of foxf2 in mice leads to a marked downregulation of pdgfrb expression in perivascular cells. In contrast to expectation, perivascular cell numbers were higher in mutant animals, but these cells did not differentiate properly. The authors use a transgenic driver line expressing gal4 under the control of the pdgfrb promoter and observe a reduction in pericyte (pdgfrb-expressing) cells in foxf2a mutants. In light of the mouse data, this result might be due to a similar downregulation of pdgfrb expression in fish, which would lead to a downregulation of gal4 expression and hence reduced labelling of pericytes. The authors show a reduction of pdgfrb expression also in zebrafish in foxf2b mutants (Chauhan et al., The Lancet Neurology 2016). It would be important to clarify whether, also in zebrafish, foxf2a/foxf2b mutants have reduced or augmented numbers of perivascular cells and how this compares to the data in the mouse. The authors should perform additional characterization of perivascular cells using marker gene expression (for a list of markers, see e.g., Shih et al. Development 2021) and/or genetic lineage tracing.

(2) The authors motivate using foxf2a mutants as a model of reduced foxf2 dosage, "similar to human heterozygous loss of FOXF2". However, it is not clear how the different foxf2 genes in zebrafish interact with each other transcriptionally. Is there upregulation of foxf2b in foxf2a mutants and vice versa? This is important to consider, as Reyahi et al. showed that foxf2 gene dosage in mice appears to be important, with an increase in foxf2 gene dosage (through transgene expression) leading to a reduction in perivascular cell numbers.

(3) Figures 3 and 4 lack data quantification. The authors describe the existence of vascular defects in adult fish, but no quantifiable parameters or quantifications are provided. This needs to be added.

(4) The analysis of pericyte phenotypes and morphologies is not clear. On page 6, the authors state: "In the wildtype brain, adult pericytes have a clear oblong cell body with long, slender primary processes that extend from the cytoplasm with secondary processes that wrap around the circumference of the blood vessel." Further down on the same page, the authors note: "In wildtype adult brains, we identified three subtypes of pericytes, ensheathing, mesh and thin-strand, previously characterized in murine models." In conclusion, not all pericytes have long, slender primary processes, but there are at least three different sub-types? Did the authors analyze how they might be distributed along different branch orders of the vasculature, as they are in the mouse? Which type of pericyte is affected in foxf2a mutant animals? Can the authors identify the branch order of the vasculature for both wildtype and mutant animals and compare which subtype of pericyte might be most affected? Are all subtypes of pericytes similarly affected in mutant animals? There also seems to be a reduction in smooth muscle cell coverage.

(5) Regarding pericyte regeneration data (Figure 7): Are the values in Figure 7D not significantly different from each other (no significance given)?

(6) In the discussion, the authors state that "pericyte processes have not been studied in zebrafish". Ando et al. (Development 2016) studied pericyte processes in early zebrafish embryos, and Leonard et al. (Development 2022) studied zebrafish pericytes and their processes in the developing fin.

Reviewer #2 (Public review):

Summary:

This study investigates the developmental and lifelong consequences of reduced foxf2 dosage in zebrafish, a gene associated with human stroke risk and cerebral small vessel disease (CSVD). The authors show that a ~50% reduction in foxf2 function through homozygous loss of foxf2a leads to a significant decrease in brain pericyte number, along with striking abnormalities in pericyte morphology-including enlarged soma and extended processes-during larval stages. These defects are not corrected over time but instead persist and worsen with age, ultimately affecting the surrounding endothelium. The study also makes an important contribution by characterizing pericyte behavior in wild-type zebrafish using a clever pericyte-specific Brainbow approach, revealing novel interactions such as pericyte process overlap not previously reported in mammals.

Strengths:

This work provides mechanistic insight into how subtle, developmental changes in mural cell biology and coverage of the vasculature can drive long-term vascular pathology. The authors make strong use of zebrafish imaging tools, including longitudinal analysis in transgenic lines to follow pericyte number and morphology over larval development, and then applied tissue clearing and whole brain imaging at 3 and 11 months to further dissect the longitudinal effects of foxf2a loss. The ability to track individual pericytes in vivo reveals cell-intrinsic defects and process degeneration with high spatiotemporal resolution. Their use of a pericyte-specific Zebrabow line also allows, for the first time, detailed visualization of pericyte-pericyte interactions in the developing brain, highlighting structural features and behaviors that challenge existing models based on mouse studies. Together, these findings make the zebrafish a valuable model for studying the cellular dynamics of CSVD.

Weaknesses:

While the findings are compelling, several aspects could be strengthened. First, quantifying pericyte coverage across distinct brain regions (forebrain, midbrain, hindbrain) would clarify whether foxf2a loss differentially impacts specific pericyte lineages, given known regional differences in developmental origin, with forebrain pericytes being neural crest-derived and hindbrain pericytes being mesoderm-derived. Second, measuring foxf2b expression in foxf2a mutants would better support the interpretation that total FOXF2 dosage is reduced in a graded fashion in heterozygote and homozygote foxf2a mutants. Finally, quantifying vascular density in adult mutants would help determine whether observed endothelial changes are a downstream consequence of prolonged pericyte loss. Correlating these vascular changes with local pericyte depletion would also help clarify causality.

Reviewer #3 (Public review):

Summary:

The goal of the work by Graff et al. is to model CSVD in the zebrafish using foxf2a mutants. The mutants show loss of cerebral pericyte coverage that persists through adulthood, but it seems foxf2a does not regulate the regenerative capacity of these cells. The findings are interesting and build on previous work from the group. Limitations of the work include little mechanistic insight into how foxf2a alters pericyte recruitment/differentiation/survival/proliferation in this context, and the overlap of these studies with previous work in fox2a/b double mutants. However, the data analysis is clean and compelling, and the findings will contribute to the field.

Reviewer #1 (Public review):

In this manuscript, Pagano and colleagues test the idea that the protein GMCL1 functions as a substrate receptor for a Cullin RING 3 E3 ubiquitin ligase (CUL3) complex. Using a pulldown approach, they identify GMCL1 binding proteins, including the DNA damage scaffolding protein 53BP1. They then focus on the idea that GMCL1 recruits 53BP1 for CUL3-dependent ubiquitination, triggering subsequent proteasomal degradation of ubiquitinated 53BP1.

In addition to its DNA damage signalling function, in mitosis, 53BP1 is reported to form a stopwatch complex with the deubiquitinating enzyme USP28 and the transcription factor p53 (PMID: 38547292). These 53BP1-stopwatch complexes generated in mitosis are inherited by G1 daughter cells and help promote p53-dependent cell cycle arrest independent from DNA damage (PMID: 38547292). Several studies show that knockout of 53BP1 overcomes G1 cell cycle arrest after mitotic delays caused by anti-mitotic drugs or centrosome ablation (PMID: 27432897, 27432896). In this model, it is crucial that 53BP1 remains stable in mitosis and more stopwatch complex is formed after delayed mitosis.

Pagano and coworkers suggest that 53BP1 levels can sometimes be suppressed in mitosis if the cells overexpress GMCL1. They carry out a bioinformatic analysis of available public data for p53 wild-type cancer cell lines resistant to the anti-mitotic drug paclitaxel and related compounds. Stratifying GMCL1 into low and high expression groups reveals a weak (p = 0.05 or ns) correlation with sensitivity to taxanes. It is unclear on what basis the authors claim paclitaxel-resistant and p53 wild-type cancer cell lines bypass the mitotic surveillance/timer pathway. They have not tested this. Figure 3 is a correlation assembled from public databases but has no experimental tests. Figure 4 looks at proliferation but not cell cycle progression or the length of mitosis. The main conclusions relating to cell cycle progression and specifically the link to mitotic delays are therefore not supported by experimental data. There is no imaging of the cell cycle or cell fate after mitotic delays, or analysis of where the cells arrest in the cell cycle. Most of the cell lines used have been reported to lack a functional mitotic surveillance pathway in the recent work by Meitinger. To support these conclusions, the stability of endogenous 53BP1 under different conditions in cells known to have a functional mitotic surveillance pathway needs to be examined. A key suggestion in the work is that the level of GMCL1 expression correlates with resistance to taxanes. For the mitotic surveillance pathway, the type of drug (nocodazole, taxol, etc) used to induce a delay isn't thought to be relevant, only the length of the delay. Do GMCL1-overexpressing cells show resistance to anti-mitotics in general?

Importantly, if GMCL1 specifically degrades 53BP1 during prolonged mitotic arrests, the authors should show what happens during normal cell divisions without any delays or drug treatments. How much 53BP1 is destroyed in mitosis under those conditions? Does 53BP1 destruction depend on the length of mitosis, drug treatment, or does 53BP1 get degraded every mitosis regardless of length? Testing the contribution of key mitotic E3 ligase activities on mitotic 53BP1 stability, such as the anaphase-promoting complex/cyclosome (APC/C) is important in this regard. One previous study reported an analysis of putative APC/C KEN-box degron motifs in 53BP1 and concluded these play a role in 53BP1 stability in anaphase (PMID: 28228263).

There is no direct test of the proposed mechanism, and it is therefore unclear if 53BP1 is ubiquitinated by a GMCL1-CUL3 ligase in cells, and how efficient this process would be at different cell cycle stages. A key issue is the lack of experimental data explaining why the proposed mechanism would be restricted to mitosis. Indirect effects, such as loss of 53BP1 from the chromatin fraction during M phase upon GMCL1 overexpression, do not necessarily mean that 53BP1 is degraded. PLK1-dependent chromatin-cytoplasmic shuttling of 53BP1 during mitotic delays has been described previously (PMID: 38547292, 37888778). These papers are cited in the text, but the main conclusions of those papers on 53BP1 incorporation into a stopwatch complex during mitotic delays have been ignored. Are the authors sure that 53BP1 is destroyed in mitosis and not simply re-localised between chromatin and non-chromatin fractions? At the very least, these reported findings should be discussed in the text.

The authors use a variety of cancer cell line models throughout their study, most of which have been reported to lack a functional mitotic surveillance pathway. U2OS and HCT116 cells do not respond normally to mitotic delays, despite being annotated as p53 WT. Other studies have used p53 wild-type hTERT RPE-1 cells to study the mitotic surveillance pathway. If the model is correct, then over-expressing GMCL1 in hTERT-RPE1 cells should suppress cell cycle arrest after mitotic delays, and GMCL1 KO should make the cells more sensitive to delays. These experiments are needed to provide an adequate test of the proposed model.

To conclude, while the authors propose a potentially interesting model on how GMCL1 overexpression could regulate 53BP1 stability to limit p53-dependent cell cycle arrest, it is unclear what triggers this pathway or when it is relevant. 53BP1 is known to function in DNA damage signalling, and GMCL1 might be relevant in that context. The manuscript contains the initial description of GMCL1-53BP1 interaction but lacks a proper analysis of the function of this interaction and is therefore a preliminary report.

Reviewer #2 (Public review):

Summary:

This study investigates the role of GMCL1 in regulating the mitotic surveillance pathway (MSP), a protective mechanism that activates p53 following prolonged mitosis. The authors identify a physical interaction between 53BP1 and GMCL1, but not with GMCL2. They propose that the ubiquitin ligase complex CRL3-GMCL1 targets 53BP1 for degradation during mitosis, thereby preventing the formation of the "mitotic stopwatch" complex (53BP1-USP28-p53) and subsequent p53 activation. The authors show that high GMCL1 expression correlates with resistance to paclitaxel in cancer cell lines that express wild-type p53. Importantly, loss of GMCL1 restores paclitaxel sensitivity in these cells, but not in p53-deficient lines. They propose that GMCL1 overexpression enables cancer cells to bypass MSP-mediated p53 activation, promoting survival despite mitotic stress. Targeting GMCL1 may thus represent a therapeutic strategy to re-sensitize resistant tumors to taxane-based chemotherapy.

Strengths:

This manuscript presents potentially interesting observations. The major strength of this article is the identification of GMCL1 as a 53BP1 interaction partner. The authors identified relevant domains and showed that GMCL1 controls 53BP1 stability. The authors further show a potentially interesting link between GMCL1 status and sensitivity to Taxol.

Weaknesses:

However, the manuscript is significantly weakened by unsubstantiated mechanistic claims, overreliance on a non-functional model system (U2OS), and overinterpretation of correlative data. To support the conclusions of the manuscript, the authors must show that the GMCL1-dependent sensitivity to Taxol depends on the mitotic surveillance pathway.

Reviewer #3 (Public review):

Summary:

In this study, Kito et al follow up on previous work that identified Drosophila GCL as a mitotic substrate recognition subunit of a CUL3-RING ubiquitin ligase (CRL3) complex.

Here they characterize mutants of the human ortholog of GCL, GMCL1, that disrupt the interaction with CUL3 (GMCL1E142K) and that lack the substrate interaction domain (GMCL1 BBO). Immunoprecipitation followed by mass spectrometry identified 9 proteins that interacted with wild-type FLAG-GMCL1 and GMCL1 EK but not GMCL1 BBO. These proteins included 53BP1, which plays a well-characterized role in double-strand break repair but also functions in a USP28-p53-53BP1 "mitotic stopwatch" complex that arrests the cell cycle after a substantially prolonged mitosis. Consistent with the IP-MS results, FLAG-GMCL1 immunoprecipitated 53BP1. Depletion of GMCL1 during mitotic arrest increased protein levels of 53BP1, and this could be rescued by wild-type GMCL1 but not the E142K mutant or a R433A mutant that failed to immunoprecipitate 53BP1.

Using a publicly available dataset, the authors identified a relatively small subset of cell lines with high levels of GMCL1 mRNA that were resistant to the taxanes paclitaxel, cabazitaxel, and docetaxel. This type of analysis is confounded by the fact that paclitaxel and other microtubule poisons accumulate to substantially different levels in various cell lines (DOI: 10.1073/pnas.90.20.9552 , DOI: 10.1091/mbc.10.4.947 ), so careful follow-up experiments are required to validate results. The correlation between increased GMCL1 mRNA and taxane resistance was not observed in lung cancer cell lines. The authors propose this was because nearly half of lung cancers harbor p53 mutations, and lung cancer cell lines with wild-type but not mutant p53 showed the correlation between increased GMCL1 mRNA and taxane resistance. However, the other cancer cell types in which they report increased GMCL1 expression correlates with taxane sensitivity also have high rates of p53 mutation. Furthermore, p53 status does not predict taxane response in patients (DOI: 10.1002/1097-0142(20000815)89:4<769::aid-cncr8>3.0.co;2-6 , DOI: 10.1002/(SICI)1097-0142(19960915)78:6<1203::AID-CNCR6>3.0.CO;2-A , PMID: 10955790).

The authors then depleted GMCL1 and reported that it increased apoptosis in two cell lines with wild-type p53 (MCF7 and U2OS) due to activation of the mitotic stopwatch. This is surprising because the mitotic stopwatch paper they cite (DOI: 10.1126/science.add9528 ) reported that U2OS cells have an inactive stopwatch and that activation of the stopwatch results in cell cycle arrest rather than apoptosis in most cell types, including MCF7. Beyond this, it has recently been shown that the level of taxanes and other microtubule poisons achieved in patient tumors is too low to induce mitotic arrest (DOI: 10.1126/scitranslmed.3007965 , DOI: 10.1126/scitranslmed.abd4811 , DOI: 10.1371/journal.pbio.3002339 ), raising concerns about the relevance of prolonged mitosis to paclitaxel response in cancer. The findings here demonstrating that GMCL1 mediates degradation of 53BP1 during mitotic arrest are solid and of interest to cell biologists, but it is unclear that these findings are relevant to paclitaxel response in patients.

Strengths:

This study identified 53BP1 as a target of CRL3GMCL1-mediated degradation during mitotic arrest. AlphaFold3 predictions of the binding interface, followed by mutational analysis, identified mutants of each protein (GMCL1 R433A and 53BP1 IEDI1422-1425AAAA) that disrupted their interaction. Knock-in of a FLAG tag into the C-terminus of GMCL1 in HCT116 cells, followed by FLAG immunoprecipitation, confirmed that endogenous GMCL1 interacts with endogenous CUL3 and 53BP1 during mitotic arrest.

Weaknesses:

The clinical relevance of the study is overinterpreted. The authors have not taken relevant data about the clinical mechanism of taxanes into account. Supraphysiologic doses of microtubule poisons cause mitotic arrest and can activate the mitotic stopwatch. However, in physiologic concentrations of clinically useful microtubule poisons, cells proceed through mitosis and divide their chromosomes on mitotic spindles that are at least transiently multipolar. Though these low concentrations may result in a brief mitotic delay, it is substantially shorter than the arrest caused by high concentrations of microtubule poisons, and the one mimicked here by 16 hours of 0.4 mg/mL nocodazole, which is not used clinically and does not induce multipolar spindles. Resistance to mitotic arrest occurs through different mechanisms than resistance to multipolar spindles. No evidence is presented in the current version of the manuscript that GMCL1 affects cellular response to clinically relevant doses of paclitaxel.

Reviewer #1 (Public review):

Strengths:

Sarpaning et al. provide a thorough characterization of putative Rnt1 cleavage of mRNA in S. cerevisiae. Previous studies have discovered Rnt1 mRNA substrates anecdotally, and this global characterization expands the known collection of putative Rnt1 cleavage sites. The study is comprehensive, with several types of controls to show that Rnt1 is required for several of these cleavages.

Weaknesses:

Formally speaking, the authors do not show a direct role of Rnt1 in mRNA cleavage - no studies were done (e.g., CLIP-seq or similar) to define direct binding sites. Is the mutant Rnt1 expected to trap substrates? Without direct binding studies, the authors rely on genetics and structure predictions for their argument, and it remains possible that a subset of these sites is an indirect consequence of rnt1. This aspect should be addressed in the discussion.

The comprehensive list of putative Rnt1 mRNA cleavage sites is interesting insofar as it expands the repertoire of Rnt1 on mRNAs, but the functional relevance of the majority of these sites remains unknown. Along these lines, the authors should present a more thorough characterization of putative Rnt1 sites recovered from in vitro Rnt1 cleavage.

The authors need to corroborate the rRNA 3'-ETS tetraloop mutations with a northern analysis of 3'-ETS processing to confirm an ETS processing defect (which might need to be done in decay mutants to stabilize the liberated ETS fragment). They state that the tetraloop mutation does not yield a growth defect and use this as the basis for concluding that rRNA cleavage is not the major role of Rnt1 in vivo, which is a surprising finding. But it remains possible that tetraloop mutations did not have the expected disruptive effect in vivo; if the ETS is processed normally in the presence of tetraloop mutations, it would undermine this interpretation. This needs to be more carefully examined.

To support the assertion that YDR514C cleavage is required for normal "homeostasis," and more specifically that it is the major contributor to the rnt1∆ growth defect, the authors should express the YDR514C-G220S mutant in the rDNA∆ strains with mutations in the 3'-ETS (assuming they disrupt ETS processing, see above). This simple experiment should provide a relative sense of "importance" for one or the other cleavage being responsible for the rnt1∆ defect. Given the accepted role of Rnt1 cleavage in rRNA processing and a dogmatic view that this is the reason for the rnt1∆ growth defect, such a result would be surprising and elevate the functional relevance and significance of Rnt1 mRNA cleavage.

Given that some Rnt1 mRNA cleavage is likely nuclear, it is possible that some of these targets are nascent mRNA transcripts, as opposed to mature but unexported mRNA transcripts, as proposed in the manuscript. A role for Rnt1 in co-transcriptional mRNA cleavage would be conceptually similar to Rnt1 cleavage of the rRNA 3'-ETS to enable RNA Pol I "torpedo" termination by Rat1, described by Proudfoot et al (PMID 20972219). To further delineate this point, the authors could e.g., examine the poly-A tails on abundant Rnt1 targets to establish whether they are mature, polyadenylated mRNAs (e.g., northern analysis of oligo-dT purified material). A more direct test would be PARE analysis of oligo-dT enriched or depleted material to determine the poly-A status of the cleavage products. Alternatively, their association with chromatin could be examined.

While laboratory strains of budding yeast have a single RNase III ortholog Rnt1, several other budding yeast have a functional RNAi system with Dcr and Ago (PMID 19745116), and laboratory yeast strains are a derived state due to pressure from the killer virus to lose the RNAi system (PMID 21921191). The current study could provide new insight into the relative substrate preferences of Rnt1 and budding yeast Dicer, which could be experimentally confirmed by expressing Dcr in RNT1 and rnt1∆ strains. In lieu of experiments, discussion of the relevance of Rnt1 cleavage compared to yeast RNAi should be included in the discussion before the "human implications" section.

For SNR84 in Figure S3D, it appears that the TSS may be upstream of the annotated gene model. Does RNA-seq coverage (from external datasets) extend upstream to these additional mapped cleavages? The assertion that the mRNA is uncapped is concerning; an alternative explanation is that the nascent mRNA has a cap initially but is subsequently cleaved by Rnt1. This point should be clarified or reworded for accuracy.

Reviewer #2 (Public review):

The yeast double-stranded RNA endonuclease Rnt1, a homolog of bacterial RNAse III, mediates the processing of pre-rRNA, pre-snRNA, and pre-snoRNA molecules. Cells lacking Rnt1 exhibit pronounced growth defects, particularly at lower temperatures. In this manuscript, Notice-Sarpaning examines whether these growth defects can be attributed at least in part to a function of Rnt1 in mRNA degradation. To test this, the authors apply parallel analysis of RNA ends (PARE), which they developed in previous work, to identify polyA+ fragments with 5' monophosphates in RNT1 yeast that are absent in rnt1Δ cells. Because such RNAs are substrates for 5' to 3' exonucleolytic decay by Rat1 in the nucleus or Xrn1 in the cytoplasm, these analyses were performed in a rat1-ts xrn1Δ background. The data recapitulate known Rtn1 cleavage sites in rRNA, snRNAs, and snoRNAs, and identify 122 putative novel substrates, approximately half of which are mRNAs. Of these, two-thirds are predicted to contain double-stranded stem loop structures with A/UGNN tetraloops, which serve as a major determinant of Rnt1 substrate recognition. Rtn1 resides in the nucleus, and it likely cleaves mRNAs there, but cleavage products seem to be degraded after export to the cytoplasm, as analysis of published PARE data shows that some of them accumulate in xrn1Δ cells. The authors then leverage the slow growth of rnt1Δ cells for experimental evolution. Sequencing analysis of thirteen faster-growing strains identifies mutations predominantly mapping to genes encoding nuclear exosome co-factors. Some of the strains have mutations in genes encoding a larat-debranching enzyme, a ribosomal protein nuclear import factor, poly(A) polymerase 1, and the RNA-binding protein Puf4. In one of the puf4 mutant strains, a second mutation is also present in YDR514C, which the authors identify as an mRNA substrate cleaved by Rnt1. Deletion of either puf4 or ydr514C marginally improves the growth of rnt1Δ cells, which the authors interpret as evidence that mRNA cleavage by Rnt1 plays a role in maintaining cellular homeostasis by controlling mRNA turnover.

While the PARE data and their subsequent in vitro validation convincingly demonstrate Rnt1-mediated cleavage of a small subset of yeast mRNAs, the data supporting the biological significance of these cleavage events is substantially less compelling. This makes it difficult to establish whether Rnt1-mediated mRNA cleavage is biologically meaningful or simply "collateral damage" due to a coincidental presence of its target motif in these transcripts.

(1) A major argument in support of the claim that "several mRNAs rely heavily on Rnt1 for turnover" comes from comparing number of PARE reads at the transcript start site (as a proxy for fraction of decapped transcripts) and at the Rnt1 cleavage site (as a proxy for fraction of Rnt1-cleaved transcripts). The argument for this is that "the major mRNA degradation pathway is through decapping". However, polyA tail shortening usually precedes decapping, and transcripts with short polyA tails would be strongly underrepresented in PARE sequencing libraries, which were constructed after two rounds of polyA+ RNA selection. This will likely underestimate the fraction of decapped transcripts for each mRNA. There is a wide range of well-established methods that can be used to directly measure differences in the half-life of Rnt1 mRNA targets in RNT1 vs rnt1Δ cells. Because the PARE data rely on the presence of a 5' phosphate to generate sequencing reads, they also cannot be used to estimate what fraction of a given mRNA transcript is actually cleaved by Rnt1.

(2) Rnt1 is almost exclusively nuclear, and the authors make a compelling case that its concentration in the cytoplasm would likely be too low to result in mRNA cleavage. The model for Rnt1-mediated mRNA turnover would therefore require mRNAs to be cleaved prior to their nuclear export in a manner that would be difficult to control. Alternatively, the Rnt1 targets would need to re-enter prior to cleavage, followed by export of the cleaved fragments for cytoplasmic decay. These processes would need to be able to compete with canonical 5' to 3' and 3' to 5' exonucleolytic decay to influence mRNA fate in a biologically meaningful way.

(3) The experimental evolution clearly demonstrates that mutations in nuclear exosome factors are the most frequent suppressors of the growth defects caused by Rnt1 loss. This can be rationalized by stabilization of nuclear exosome substrates such as misprocessed snRNAs or snoRNAs, which are the major targets of Rnt1. The rescue mutations in other pathways linked to ribosomal proteins (splicing, ribosomal protein import, ribosomal mRNA binding) support this interpretation. By contrast, the potential suppressor mutation in YDR514C does not occur on its own but only in combination with a puf4 mutation; it is also unclear whether it is located within the Rnt1 cleavage motif or if it impacts Rnt1 cleavage at all. This can easily be tested by engineering the mutation into the endogenous YDR514C locus with CRISPR/Cas9 or expressing wild-type and mutant YDR514C from a plasmid, along with assaying for Rnt1 cleavage by northern blot. Notably, the growth defect complementation of YDR514C deletion in rnt1Δ cells is substantially less pronounced than the growth advantage afforded by nuclear exosome mutations (Figure S9, evolved strains 1 to 5). These data rather argue for a primary role of Rnt1 in promoting cell growth by ensuring efficient ribosome biogenesis through pre-snRNA/pre-snoRNA processing.

Reviewer #1 (Public review):

Summary:

Genome-wide association studies have been an important approach to identifying the genetic basis of human traits and diseases. Despite their successes, for many traits, a substantial amount of variation cannot be explained by genetic factors, indicating that environmental variation and individual 'noise' (stochastic differences as well as unaccounted for environmental variation) also play important roles. The authors' goal was to address whether gene expression variation in genetically identical individuals, driven by historical environmental differences and 'noise', could be used to predict reproductive trait differences.

Strengths:

To address this question, the authors took advantage of genetically identical C. elegans individuals to transcriptionally profile 180 adult hermaphrodite individuals that were also measured for two reproductive traits. A major strength of the paper is its experimental design. While experimenters aim to control the environment that each worm experiences, it is known that there are small differences that each worm experiences even when they are grown together on the same agar plate - e.g. the age of their mother, their temperature, the amount of food they eat, and the oxygen and carbon dioxide levels depending on where they roam on the plate. Instead of neglecting this unknown variation, the authors design the experiment up front to create two differences in the historical environment experienced by each worm: 1) the age of its mother and 2) 8 8-hour temperature difference, either 20 or 25 {degree sign}C. This helped the authors interpret the gene expression differences and trait expression differences that they observed.

Using two statistical models, the authors measured the association of gene expression for 8824 genes with the two reproductive traits, considering both the level of expression and the historical environment experienced by each worm. Their data supports several conclusions. They convincingly show that gene expression differences are useful for predicting reproductive trait differences, predicting ~25-50% of the trait differences depending on the trait. Using RNAi, they also show that the genes they identify play a causal role in trait differences. Finally, they demonstrate an association with trait variation and the H3K27 trimethylation mark, suggesting that chromatin structure can be an important causal determinant of gene expression and trait variation.

Overall, this work supports the use of gene expression data as an important intermediate for understanding complex traits. This approach is also useful as a starting point for other labs in studying their trait of interest.

Weaknesses:

There are no major weaknesses that I have noted. Some important limitations of the work (that I believe the authors would agree with) are worth highlighting, however:

(1) A large remaining question in the field of complex traits remains in splitting the role of non-genetic factors between environmental variation and stochastic noise. It is still an open question which role each of these factors plays in controlling the gene expression differences they measured between the individual worms.

(2) The ability of the authors to use gene expression to predict trait variation was strikingly different between the two traits they measured. For the early brood trait, 448 genes were statistically linked to the trait difference, while for egg-laying onset, only 11 genes were found. Similarly, the total R2 in the test set was ~50% vs. 25%. It is unclear why the differences occur, but this somewhat limits the generalizability of this approach to other traits.

(3) For technical reasons, this approach was limited to whole worm transcription. The role of tissue and cell-type expression differences is important to the field, so this limitation is important.

Reviewer #2 (Public review):

Summary:

This paper measures associations between RNA transcript levels and important reproductive traits in the model organism C. elegans. The authors go beyond determining which gene expression differences underlie reproductive traits, but also (1) build a model that predicts these traits based on gene expression and (2) perform experiments to confirm that some transcript levels indeed affect reproductive traits. The clever study design allows the authors to determine which transcript levels impact reproductive traits, and also which transcriptional differences are driven by stochastic vs environmental differences. In sum, this is a rather comprehensive study that highlights the power of gene expression as a driver of phenotype, and also teases apart the various factors that affect the expression levels of important genes.

Strengths:

Overall, this study has many strengths, is very clearly communicated, and has no substantial weaknesses that I can point to. One question that emerges for me is about the extent to which these findings apply broadly. In other words, I wonder whether gene expression levels are predictive of other phenotypes in other organisms. I think this question has largely been explored in microbes, where some studies (PMID: 17959824) but not others (PMID: 38895328) find that differences in gene expression are predictive of phenotypes like growth rate. Microbes are not the primary focus here, and instead, the discussion is mainly focused on using gene expression to predict health and disease phenotypes in humans. This feels a little complicated since humans have so many different tissues. Perhaps an area where this approach might be useful is in examining infectious single-cell populations (bacteria, tumors, fungi). But I suppose this idea might still work in humans, assuming the authors are thinking about targeting specific tissues for RNAseq.

In sum, this is a great paper that really got me thinking about the predictive power of gene expression and where/when it could inform about (health-related) phenotypes.

Reviewer #3 (Public review):

Summary:

Webster et al. sought to understand if phenotypic variation in the absence of genetic variation can be predicted by variation in gene expression. To this end they quantified two reproductive traits, the onset of egg laying and early brood size in cohorts of genetically identical nematodes exposed to alternative ancestral (two maternal ages) and same generation life histories (either constant 20C temperature or 8-hour temperature shift to 25C upon hatching) in a two-factor design; then they profiled genome-wide gene expression in each individual.

Using multiple statistical and machine learning approaches, they showed that, at least for early brood size, phenotypic variation can be quite well predicted by molecular variation, beyond what can be predicted by life history alone.

Moreover, they provide some evidence that expression variation in some genes might be causally linked to phenotypic variation.

Strengths:

(1) Cleverly designed and carefully performed experiments that provide high-quality datasets useful for the community.

(2) Good evidence that phenotypic variation can be predicted by molecular variation.

Weaknesses:

What drives the molecular variation that impacts phenotypic variation remains unknown. While the authors show that variation in expression of some genes might indeed be causal, it is still not clear how much of the molecular variation is a cause rather than a consequence of phenotypic variation.

Reviewer #1 (Public review):

Summary:

Even though mutations in LRRK2 and GBA1 (which encodes the protein GCase) increase the risk of developing Parkinson's disease (PD), the specific mechanisms driving neurodegeneration remain unclear. Given their known roles in lysosomal function, the authors investigate how LRRK2 and GCase activity influence the exocytosis of the lysosomal lipid BMP via extracellular vesicles (EVs). They use fibroblasts carrying the PD-associated LRRK2-R1441G mutation and pharmacologically modulate LRRK2 and GCase activity.

Strengths:

The authors examine both proteins at endogenous levels, using MEFs instead of cancer cells. The study's scope is potentially interesting and could yield relevant insights into PD disease mechanisms.

Weaknesses:

Many of the authors' conclusions are overstated and not sufficiently supported by the data. Several statistical errors undermine their claims. Pharmacological treatment is very long, leading to potential off-target effects. Additionally, the authors should be more rigorous when using EV markers.

Reviewer #2 (Public review):

Summary:

In this paper, the authors used MEFs expressing the R1441G mutant of leucine-rich repeat kinase 2 (LRRK2), a mutant associated with the early onset of Parkinson's disease. They report that in these cells LAMP2 fluorescence is higher but BMP fluorescence is lower, MVE size is reduced, and that MVEs contain less ILVs. They also report that LAMP2-positive EVs are increased in mutant cells in a process sensitive to LRRK2 kinase inhibition but are further increased by glucocerebrosidase (GCase) inhibition, and that total di-22:6-BMP and total di-18:1-BMP are increased in mutant LRRK2 MEFs compared to WT cells by mass spectrometry. They also report that LRRK2 kinase inhibition partially restores cellular BMP levels, and that GCase inhibition further increases BMP levels, and that in EVs from the LRRK2 mutant, LRRK2 inhibition decreases BMP while GCase inhibition has the opposite effect. Moreover, they report that the BMP increase is not due to increased BMP synthesis, although the authors observe that CLN5 is increased in LRRK2 mutant cells. Finally, they report that GW4869 decreases EV release and exosomal BMP, while bafilomycin A1 increases EV release. They conclude that LRRK2 regulates BMP levels (in cells) and release (via EVs). They also conclude that the process is modulated by GCase in LRRK2 mutant cells, and that these studies may contribute to the use of BMP-positive EVs as a biomarker for Parkinson's disease and associated treatments.

Strengths:

This is an interesting paper, which provides novel insights into the biogenesis of exosomes with exciting biomedical potential. However, I have comments that authors need to address to clarify some aspects of their study.

Weaknesses:

(1) The intensity of LAMP2 staining is increased significantly in cells expressing the R1441G mutant of LRRK2 when compared to WT cells (Figure 1C). Yet mutant cells contain significantly smaller MVEs with fewer ILVs, and the MVE surface area is reduced (Figure 1D-F). This is quite surprising since LAMP2 is a major component of the limiting membrane of late endosomes. Are other proteins of endo-lysosomes (eg, LAMP1, CD63, RAB7) or markers (lysotracker) also decreased (see also below)?

(2) LRRK2 has been reported to interact with endolysosomal membranes. Does the R1441G mutant bind LAMP2- and/or BMP-positive membranes? Does the mutant affect endolysosomes?

(3) Immunofluorescence data indicate that BMP is decreased in mutant LRRK2-expressing cells compared to WT (Figure 1A-B), but mass spec data indicate that di-22:6-BMP and di-18:1-BMP are increased (Figure 3). Authors conclude that the BMP pool detected by mass spec in mutant cells is less antibody-accessible than that present in wt cells, or that the anti-BMP antibody is less specific and that it detects other analytes. This is an awkward conclusion, since the IF signal with the antibody is lower (not higher): why would the antibody be less specific? Could it be that the antibody does not see all BMP isoforms equally well? Moreover, the observations that mutant cells contain smaller MVEs (Figure 1D-F) with fewer ILVs are consistent with the IF data and reduced BMP amounts. This needs to be clarified.

Mass spectrometry data are only shown for two BMP species (di-22:6, di-18:1). What are the major BMP isoforms in WT cells? The authors should show the complete analysis for all BMP species if they wish to draw quantitative conclusions about the amounts of BMP in wt and mutant cells. Finally, BMP and PG are isobaric lipids. Fragmentation of BMPs or PGs results in characteristic fingerprints, but the presence of each daughter ion is not absolutely specific for either lipid. This should be clarified, e.g., were BMP and PG separated before mass spec analysis? Was PG affected? The authors should also compare the BMP data with mass spec data obtained with a control lipid, e.g., PC.

(4) It is quite surprising that the amounts of labeled BMP continue to increase for up to 24h after a short 25min pulse with heavy BMP precursors (Figure 4B).

(5) It is argued that upregulation of CLN5 may be due to an overall upregulation of lysosomal enzymes, as LAMP2 levels were also increased (Figure 2A, C, E). Again, this is not consistent with the observed decrease in MVE size and number (Figure 1D-F). As mentioned above, other independent markers of endo-lysosomes should be analyzed (eg, LAMP1, CD63, RAB7), and/or other lysosomal enzymes (e.g. cathepsin. D).

(6) The authors report that the increase in BMP is not due to an increase in BMP synthesis (Figure 4), although they observe a significant increase in CLN5 (Figure 5A) in LRRK2 mutant cells. Some clarification is needed.

(7) Authors observe that both LAMP2 and BMP are decreased in EVs by GW4869 and increased by bafilomycin (Figure 6). Given my comments above on Figure 1, it would also be nice to illustrate/quantify the effects of these compounds on cells by immunofluorescence.

Reviewer #1 (Public review):

Summary:

This paper seeks to understand the upstream regulation and downstream effectors of glycolysis in retinal progenitor cells, using mouse retinal explants as the main model system. The paper presents evidence that high glycolysis in retinal progenitor cells is required for their proliferation and timely differentiation into photoreceptors. Retinal glycolysis increases after deletion of Pten. The authors suggest that high glycolysis controls cell proliferation and differentiation by promoting intracellular alkalinization, beta-catenin acetylation and stabilization and consequent activation of the canonical Wnt pathway.

Strengths:

- The experiments showing that PFKFB3 overexpression is sufficient to increase proliferation of retinal progenitors (which are already highly dividing cells) and photoreceptor differentiation are striking and the result unanticipated. It suggests that glycolytic flux is normally limiting for proliferation in embryos.<br /> - Likewise the result that an increase in pH from 7.4 to 8.0 is sufficient to increase proliferation implies that pH regulation may have instructive roles in setting the tempo of retinal development and embryonic cell proliferation. Similarly for the results showing that acetate supplementation increases proliferation (I think this result should be moved to the main figures).

Weaknesses:

- Epistatic experiments to test if changes in pH mediate the effects of glycolysis on photoreceptor differentiation, or if Wnt activation is the main downstream effector of glycolysis in controlling differentiation are not presented.<br /> - It is likely that metabolism changes ex vivo vs in vivo, and therefore stable isotope tracing experiments in the explants may not reflect in vivo metabolism.<br /> - The retina at P0 is composed of both progenitors and differentiated cells. It is not clear if the results of the RNA-seq and metabolic analysis reflect changes in the metabolism of progenitors, or of mature cells, or changes in cell type composition rather than direct metabolic changes in a specific cell type.<br /> - The biochemical links between elevated glycolysis and pH and beta-catenin stability are unclear. White et al found that higher pH decreased beta-catenin stability (JCB 217: 3965) in contrast to the results here. Oginuma et al found that inhibition of glycolysis or beta-catenin acetylation does not affect beta-catenin stability (Nature 584:98), again in contrast to these results. Another paper showed that acidification inhibits Wnt signaling by promoting the expression of a transcriptional repressor and not via beta-catenin stability (Cell Discovery 4:37). There are also additional papers showing increased pH can promote cell proliferation via other mechanisms (e.g. Nat Metab 2:1212). It is possible that there is organ-specificity in these signaling pathways however some clarification of these divergent results is warranted.<br /> - The gene expression analysis is not completely convincing. E.g. expression of additional glycolytic genes should be shown in Fig. 1. It is not clear why Hk1 and Pgk1 are specifically shown, and conclusions about changes in glycolysis are difficult to draw from expression of these two genes. The increase in glycolytic gene expression in the Pten-deficient retina is generally small.<br /> - Is it possible that glycolytic inhibition with 2DG slows down development and production of most new differentiated cells rather than specifically affecting photoreceptor differentiation?<br /> - Are the prematurely-born cells caused by PFKFB3 overexpression photoreceptors as assessed by morphology or markers (in addition to position)?

Reviewer #3 (Public review):

Summary:

This study examines the metabolic regulation of progenitor proliferation and differentiation in the developing retina. The authors observe dynamic changes in glycolytic gene expression in retinal progenitors and use various strategies to test the role of glycolysis. They find that elevated glycolysis in Pten-cKO retinas results in alteration of RPC fate, while inhibition of glycolysis has converse effects. They specifically test the role of elevated glycolysis using dominant active cytoPFKB3, which demonstrates the selective effects of elevated glycolysis on progenitor proliferation and rod differentiation. They then show that elevated glycolysis modulates both pHi and Wnt signaling, and provide evidence that these pathways impact proliferation and differentiation of progenitors, particularly affecting rod photoreceptor differentiation.

Strengths:

This is a compelling and rigorous study that provides an important advance in our understanding of metabolic regulation of retina development, addressing a major gap in knowledge. A key strength is that the study utilizes multiple genetic and pharmacological approaches to address how both increased or decreased glycolytic flux affect retinal progenitor proliferation and differentiation. They discover elevated Wnt signaling pathway genes in Pten cKO retina, revealing a potential link between glycolysis and Wnt pathway activation. Altogether the study is comprehensive and adds to the growing body of evidence that regulation of glycolysis plays a key role in tissue development.

Weaknesses:

(1) Following expression of cytoPFKB3, which results in increased glycolytic flux, BrDU labeling was performed at e12.5 and increased labeled cells were detected in the outer nuclear layer. But whether these are cones or rods is not established. The rest of the analysis is focused on the precocious maturation of rhodopsin-labelled outer segments, and the major conclusions emphasize rod photoreceptor differentiation. Therefore it is unclear whether there is an effect on cone differentiation for either Pten cKO or cytoPFKB3 transgenic retina. It is also not established whether rods are born precociously. Presumably this would be best detected by BrDU labeling at later embryonic stages.

(2) The authors find that there is upregulation of multiple Wnt pathway components in Pten cKO retina. They further show that inhibiting Wnt signaling phenocopies the effects of reducing glycolysis. However, they do not test whether pharmacological inhibition of Wnt signaling reverses the effects of high glycolytic activity in Pten cKO retinas. Thus the argument that Wnt is a key downstream effector pathway regulating rod photoreceptor differentiation is weak.

(3) The use of sodium acetate to force protein acetylation is quite non-specific and will have effects beyond beta-catenin acetylation (which the authors acknowledge). Thus it is a stretch to state that "forced activation of beta-catenin acetylation" mimics the impact of Pten<br /> loss/high glycolytic activity in RPCs since the effects could be due to acetylation of other proteins.

Reviewer #1 (Public review):

Summary:

The paper describes the cryoEM structure of RAD51 filament on the recombination intermediate. In the RAD51 filament, the insertion of a DNA-binding loop called the L2 loop stabilizes the separation of the complementary strand for the base-pairing with an incoming ssDNA and the non-complementary strand, which is captured by the second DNA-binding channel called the site II. The molecular structure of the RAD51 filament with a recombination intermediate provides a new insight into the mechanism of homology search and strand exchange between ssDNA and dsDNA.

Strengths:

This is the first human RAD51 filament structure with a recombination intermediate called the D-loop. The work has been done with great care, and the results shown in the paper are compelling based on cryo-EM and biochemical analyses. The paper is really nice and important for researchers in the field of homologous recombination, which gives a new view on the molecular mechanism of RAD51-mediated homology search and strand exchange.

Weaknesses:

The authors need more careful text writing. Without page and line numbers, it is hard to give comments.

Reviewer #2 (Public review):

Summary:

Homologous recombination (HR) is a critical pathway for repairing double-strand DNA breaks and ensuring genomic stability. At the core of HR is the RAD51-mediated strand-exchange process, in which the RAD51-ssDNA filament binds to homologous double-stranded DNA (dsDNA) to form a characteristic D-loop structure. While decades of biochemical, genetic, and single-molecule studies have elucidated many aspects of this mechanism, the atomic-level details of the strand-exchange process remained unresolved due to a lack of atomic-resolution structure of RAD51 D-loop complex.<br /> In this study, the authors achieved this by reconstituting a RAD51 mini-filament, allowing them to solve the RAD51 D-loop complex at 2.64 Å resolution using a single particle approach. The atomic resolution structure reveals how specific residues of RAD51 facilitate the strand exchange reaction. Ultimately, this work provides unprecedented structural insight into the eukaryotic HR process and deepens the understanding of RAD51 function at the atomic level, advancing the broader knowledge of DNA repair mechanisms.

Strengths:

The authors overcame the challenge of RAD51's helical symmetry by designing a minifilament system suitable for single-particle cryo-EM, enabling them to resolve the RAD51 D-loop structure at 2.64 Å without imposed symmetry. This high resolution revealed precise roles of key residues, including F279 in Loop 2, which facilitates strand separation, and basic residues on site II that capture the displaced strand. Their findings were supported by mutagenesis, strand exchange assays, and single-molecule analysis, providing strong validation of the structural insights.

Weaknesses:

Despite the detailed structural data, some structure-based mutagenesis data interpretation lacks clarity. Additionally, the proposed 3′-to-5′ polarity of strand exchange relies on assumptions from static structural features, such as stronger binding of the 5′-arm-which are not directly supported by other experiments. This makes the directional model compelling but contradicts several well-established biochemical studies that support a 5'-to-3' polarity relative to the complementary strand (e.g., Cell 1995, PMID: 7634335; JBC 1996, PMID: 8910403; Nature 2008, PMID: 18256600).

Overall:

The 2.6 Å resolution cryoEM structure of the RAD51 D-loop complex provides remarkably detailed insights into the residues involved in D-loop formation. The high-quality cryoEM density enables precise placement of each nucleotide, which is essential for interpreting the molecular interactions between RAD51 and DNA. Particularly, the structural analysis highlights specific roles for key domains, such as the N-terminal domain (NTD), in engaging the donor DNA duplex.

This structural interpretation is further substantiated by single-molecule fluorescence experiments using the KK39,40AA NTD mutant. The data clearly show a significant reduction in D-loop formation by the mutant compared to wild-type, supporting the proposed functional role of the NTD observed in the cryoEM model.

However, the strand exchange activity interpretation presented in Figure 5B could benefit from a more rigorous experimental design. The current assay measures an increase in fluorescence intensity, which depends heavily on the formation of RAD51-ssDNA filaments. As shown in Figure S6A, several mutants exhibit reduced ability to form such filaments, which could confound the interpretation of strand exchange efficiency. To address this, the assay should either: (1) normalize for equivalent levels of RAD51-ssDNA filaments across samples, or (2) compare the initial rates of fluorescence increase (i.e., the slope of the reaction curve), rather than endpoint fluorescence, to better isolate the strand exchange activity itself.

Based on the structural features of the D-loop, the authors propose that strand pairing and exchange initiate at the 3'-end of the complementary strand in the donor DNA and proceed with a 3'-to-5' polarity. This conclusion, drawn from static structural observations, contrasts with several well-established biochemical studies that support a 5'-to-3' polarity relative to the complementary strand (e.g., Cell 1995, PMID: 7634335; JBC 1996, PMID: 8910403; Nature 2008, PMID: 18256600). While the structural model is compelling and methodologically robust, this discrepancy underscores the need for further experiments.

Reviewer #3 (Public review):

Summary:

Built on their previous pioneer expertise in studying RAD51 biology, in this paper, the authors aim to capture and investigate the structural mechanism of human RAD51 filament bound with a displacement loop (D-loop), which occurs during the dynamic synaptic state of the homologous recombination (HR) strand-exchange step. As the structures of both pre- and post-synaptic RAD51 filaments were previously determined, a complex structure of RAD51 filaments during strand exchange is one of the key missing pieces of information for a complete understanding of how RAD51 functions in the HR pathway. This paper aims to determine the high-resolution cryo-EM structure of RAD51 filament bound with the D-loop. Combined with mutagenesis analysis and biophysical assays, the authors aim to investigate the D-loop DNA structure, RAD51-mediated strand separation and polarity, and a working model of RAD51 during HR strand invasion in comparison with RecA.

Strengths:

(1) The structural work and associated biophysical assays in this paper are solid, elegantly designed, and interpreted.  These results provide novel insights into RAD51's function in HR.

(2) The DNA substrate used was well designed, taking into consideration the nucleotide number requirement of RAD51 for stable capture of donor DNA. This DNA substrate choice lays the foundation for successfully determining the structure of the RAD51 filament on D-loop DNA using single-particle cryo-EM.

(3) The authors utilised their previous expertise in capping DNA ends using monomeric streptavidin and combined their careful data collection and processing to determine the cryo-EM structure of full-length human RAD51 bound at the D-loop in high resolution. This interesting structure forms the core part of this work and allows detailed mapping of DNA-DNA and DNA-protein interaction among RAD51, invading strands, and donor DNA arms (Figures 1, 2, 3, 4). The geometric analysis of D-loop DNA bound with RAD51 and EM density for homologous DNA pairing is also impressive (Figure S5). The previously disordered RAD51's L2-loop is now ordered and traceable in the density map and functions as a physical spacer when bound with D-loop DNA. Interestingly, the authors identified that the side chain position of F279 in the L2_loop of RAD51_H differs from other F279 residues in L2-loops of E, F, and G protomers. This asymmetric binding of L2 loops and RAD51_NTD binding with donor DNA arms forms the basis of the proposed working model about the polarity of csDNA during RAD51-mediated strand exchange.

(4) This work also includes mutagenesis analysis and biophysical experiments, especially EMSA, single-molecule fluorescence imaging using an optical tweezer, and DNA strand exchange assay, which are all suitable methods to study the key residues of RAD51 for strand exchange and D-loop formation (Figure 5).

Weaknesses:

(1) The proposed model for the 3'-5' polarity of RAD51-mediated strand invasion is based on the structural observations in the cryo-EM structure. This study lacks follow-up biochemical/biophysical experiments to validate the proposed model compared to RecA or developing methods to capture structures of any intermediate states with different polarity models.

(2) The functional impact of key mutants designed based on structure has not been tested in cells to evaluate how these mutants impact the HR pathway.

The significance of the work for the DNA repair field and beyond:

Homologous recombination (HR) is a key pathway for repairing DNA double-strand breaks and involves multiple steps. RAD51 forms nucleoprotein filaments first with 3' overhang single-strand DNA (ssDNA), followed by a search and exchange with a homologous strand. This function serves as the basis of an accurate template-based DNA repair during HR. This research addressed a long-standing challenge of capturing RAD51 bound with the dynamic synaptic DNA and provided the first structural insight into how RAD51 performs this function. The significance of this work extends beyond the discovery of biology for the DNA repair field, into its medical relevance. RAD51 is a potential drug target for inhibiting DNA repair in cancer cells to overcome drug resistance. This work offers a structural understanding of RAD51's function with the D-loop and provides new strategies for targeting RAD51 to improve cancer therapies.

Reviewer #1 (Public review):

Summary:

This work contributes several important and interesting observations regarding the heterotolerance of non-growing Escherichia coli and Pseudomonas aeruginosa to the antimicrobial peptide tachyplesin. The primary mechanism of action of tachyplesin is thought to be disruption of the bacterial cell envelope, leading to leakage of cellular contents after a threshold level of accumulation. Although the MIC for tachyplesin in exponentially growing E. coli is just 1 ug/ml, the authors observe that a substantial fraction of a stationary phase population of bacteria survives much higher concentrations, up to 64 ug/ml. By using a fluorescently labelled analogue of tachyplesin, the authors show that the amount of per-cell intracellular accumulation of tachyplesin displays a bimodal distribution, and that the fraction of "low accumulators" correlates with the fraction of survivors. Using a microfluidic device, they show that low accumulators exclude propidium iodide, suggesting that their cell envelopes remain largely intact, while high accumulators of tachyplesin also stain with propidium iodide. They show that this phenomenon holds for several clinical isolates of E. coli with different genetic determinants of antibiotic resistance, and for a strain of Pseudomonas aeruginosa. However, the bimodal distribution does not occur in these organisms for several other antimicrobial peptides, or for tachyplesin in Klebsiella pneumoniae or Staphylococcus aureus, indicating some degree of specificity in the interaction between AMP and bacterial cell envelope. They next explore the dynamics of the fluorescent tachyplesin accumulation and show interestingly that a high degree of accumulation is initially seen in all cells, but that the "low accumulator" subpopulation manages to decrease the amount of intracellular fluorescence over time, while the "high accumulator"subpopulation continues to increase its intracellular fluorescence. Focusing on increased efflux as a hypothesised mechanism for the "low accumulator" phenotype, based on transcriptomic analysis of the two subpopulations, the authors screen putative efflux inhibitors to see if they can block the formation of the low accumulator subpopulation. They find that both the protonophore CCCP and the SSRI sertraline can block the formation of this subpopulation and that a combination of sertraline plus tachyplesin kills a greater fraction of the stationary phase cells than either agent alone, similar to the killing observed when growing cells are treated with tachyplesin.

Strengths:

This study provides new insight into the heterogeneous behaviours of non-growing bacteria when exposed to an antimicrobial peptide, and into the dynamics of their response. The single-cell analysis by FACS and microscopy is compelling. The results provide a much-needed single cell perspective on the phenomenon of tolerance to AMPs and a good starting point for further exploration.

Weaknesses:

The authors have substantially improved the clarity of the manuscript and have added additional experiments to probe further the location of the AMP relative to low and high accumulators, and the physiological states of these sub-populations. These experiments strengthen the assertion that low accumulators keep the AMP at the cell surface while high accumulators permit intracellular access to the AMP.

However, many questions still remain about the physiological characterisation of the "low accumulator" cells. While the evidence presented does support an induced response that removes the AMP from the interior of the cell, no clear mechanism for this is favoured by the experiments presented.

A double deletion of acrA and tolC (two out of the three components of the major constitutive RND efflux pump) reduces the appearance of the low accumulator phenotype, but interestingly, the single deletions have no effect, and a well-characterised inhibitor of RND efflux pumps also has no effect. The authors identify a two-component system, qseCB, that appears necessary for the appearance of low accumulators, but this system has pleiotropic effects on many cellular systems, with only tenuous connections to efflux. The selected pharmacological agents that could prevent the appearance of low accumulators do not offer clear insight into the mechanism by which low accumulators arise, because they have diverse modes of action.

The transcriptomics data collected for low and high accumulator sub-populations are interesting, but in my opinion, the conclusions that can be drawn from these data remain overstated. It is not possible to make any claims about the total amount of "protein synthesis, energy production, and gene expression" on the basis of RNA-Seq data. The reads from each sample are normalised, so there is no information about the total amount of transcript. Many elements of total cellular activity are post-transcriptionally regulated, so it is impossible to assess from transcriptomics alone. Finally, the transcriptomic data are analysed in aggregated clusters of genes that are enriched for biological processes, for example: "Cluster 2 included processes involved in protein synthesis, energy production, and gene expression that were downregulated to a greater extent in low accumulators than high accumulators". However, this obscures the fact that these clusters include genes that are generally inhibitory of the process named, as well as genes that facilitate the process.

The authors have added an experiment to attempt to assess overall metabolic activity in the low accumulator and high accumulator populations, which is a welcome addition. They apply the redox dye resazurin and observe lower resorufin (reduced form) fluorescence in the low accumulator population, which they take to indicate a lower respiration rate. This seems possible, however, an important caveat is that they have shown the low accumulator population to retain substantially lower amounts of multiple different fluorescent molecules (tachyplesin-NBD, propidium iodide, ethidium bromide) intracellularly compared to the high accumulator population. It seems possible that the low accumulator population is also capable of removing resazurin or resorufin from the intracellular space, regardless of metabolic rate. Indeed, it has previously been shown that efflux by RND efflux pumps influences resazurin reduction to resorufin in both P. aeruginosa and E. coli. By measuring only the retained redox dye using flow cytometry, the results may be confounded by the demonstrated ability of the low accumulator population to remove various fluorescent dyes. More work is needed to strongly support broad conclusions about the physiological states of the low and high accumulator populations.

The phenomenon of the emergence of low accumulators, which are phenotypically tolerant to the antimicrobial peptide tachyplesin, is interesting and important even if there is still work to be done to understand the mechanism by which it occurs.

Reviewer #2 (Public review):

Summary:

This study reports on the existence of subpopulations of isogenic E. coli and P. aeruginosa cells that are tolerant to the antimicrobial peptide tachyplesin and are characterized by accumulation of low levels of a fluorescent tachyplesin-NBD conjugate. The authors then set out to address the molecular mechanisms, providing interesting insights even though the mechanism remains incompletely defined: The work demonstrates that increased efflux may cause this phenotype, putatively together with other changes in membrane lipid composition. The authors further demonstrate that pharmacological manipulation can prevent generation of tolerance. The authors are cautious in their interpretation and the claims made are largely justified by the data.

Strengths:

Going beyond the commonly used bulk techniques for studying susceptibility to AMPs , Lee et al. used of fluorescent antibiotic conjugates in combination with flow cytometry analysis to study variability in drug accumulation at the single cell level. This powerful approach enabled the authors to expose bimodal drug accumulation pattern that were condition dependent, but conserved across a variety of E. coli clinical isolates. Using cell sorting in combination with colony-forming unit assays as well as quantitative fluorescence microscopic analysis in a microfludics-setup the authors compellingly demonstrate that low accumulators (where fluorescence signal is mostly restricted to the membrane), can survive antibiotic treatment, whereas high accumulators (with high intracellular fluorescence) were killed.

The relevance of efflux for the ´low accumulator´ phenotype and its survival is convincingly demonstrated by the following lines of evidence: i) A time-course experiment on tachyplesin-NBD pre-loaded cells revealed that all cells initially were high accumulators, before a subpopulation of cells subsequently managed to reduce signal intensity, demonstrating that the ´low accumulator´ phenotype is an induced response and not a pre-existing property. Ii) Double-mutants deficient in the delta acrA delta tolC double-KO, which showed reduced levels of low accumulators´. Interestingly, ´low accumulator´populations were nearly abrogated in bacteria deficient in the qse quorum sensing system, suggesting its centrality for the tachyplesin response. Even though this system may control acrA, the strength of the phenotype may suggest that it may control additional as-of-yet unidenitified factors relevant in the response to tachyplesin. Iii) treatment with efflux pump inhibitor sertraline and verapamil (even though some caution needs to be taken since it is not perfectly selective, see weakness) prevents generation of low accumulators. The observation that sertraline enhances tachyplesin-based killing is an important basis for developing combination therapies.

The study convincingly illustrates how susceptibility to tachyplesin adaptively changes in a heterogeneous way dependent on the growth phases and nutrient availability. This is highly relevant also beyond the presented example of tachyplesin and similar subpopulation-based adaptive changes to the susceptibility towards antimicrobial peptides or other drugs may occur during infections in vivo and they would likely be missed by standardized in vitro susceptibility testing.

Weaknesses:

Some mechanistic questions regarding tachyplesin-accumulation and survival remain. One general shortcoming of the setup of the transcriptomics experiment is that the tachyplesin-NBD probe itself has antibiotic efficacy and induces phenotypes (and eventually cell death) in the ´high accumulator´ cells. As the authors state themselves, this makes it challenging to interpret whether any differences seen between the two groups are causative for the observed accumulation pattern of if they are a consequence of differential accumulation and downstream phenotypic effects.

I have a few minor concerns regarding new data that was added during the revision:

- The statement ´ Moreover, we found that the fluorescence of low accumulators decreased over time when bacteria were treated with 20 μg mL´ is, in my opinion, not supported by the data shown in Figure S4C. That figure shows that the abundance of ´low accumulator´ cells decreases over time. Following the rationale that protease K treatment may cleave surface-associated/extracellular tachyplesin-NDB, this should lead to a shift of ´low accumulator´population to the left, indicating reduced fluorescence intensity per cell. This is not so case, but the population just disappears. However, after 120 min of treatment more cells appear in the ´high accumulator´ state. This result is somewhat puzzling.

- The authors used the metabolic dye resazurin to measure the metabolic activity of low vs. high accumulators. I am not entirely convinced that the lower fluorescence resorufin-fluorescence in tachyplesin-NBD accumulators really indicates lower metabolic activity, since a cell's fluorescence levels would also be affected by the cellular uptake and efflux. It appears plausible that the lower resorufin-fluorescence may result from reduced accumulation/increased efflux in the´low-tachyplesin NBD´ population.

Comment on revisions: All my previous comments have been satisfactorily addressed by the authors.

Reviewer #3 (Public review):

Summary:

This important study shows that stationary phase bacteria survive antimicrobial peptide treatment by switching on efflux pumps, generating low accumulating subpopulations that evade killing-a finding with clear implications for the design of peptide based antibiotics and for researchers studying antimicrobial resistance. The evidence is solid and frequently convincing, as diverse single cell assays, genetics and chemical inhibition coherently link reduced intracellular peptide to survival, even though a few mechanistic details warrant further exploration.

Strengths:

The authors investigate how Escherichia coli (and, to a lesser extent, Pseudomonas aeruginosa) survive exposure to the antimicrobial peptide (AMP) tachyplesin. Because resistance to AMPs is thought to rely heavily on non genetic adaptations rather than on classical mutation based mechanisms, the study focuses on phenotypic heterogeneity and seeks to pinpoint the cellular processes that protect a subset of cells. Using fluorescently labelled tachyplesin, single cell imaging, flow cytometry, transcriptomics, targeted genetics, and chemical perturbations, the authors report that stationary phase cultures harbor two phenotypic states: high accumulating cells that die and low accumulating cells that survive. They further propose and show that inducible efflux activity is the primary driver of survival and show that either efflux inhibition (sertraline, verapamil) or nutrient supplementation prevents the emergence of low accumulators and boosts killing.<br /> The experiments unambiguously reveal that the cells respond to stress heterogeneously, with two distinct subpopulations - one with better survival than the other. This primary phenotype is convincingly shown across various E. coli strains, including clinical isolates. The authors probed the underlying mechanism from several angles, with important additional experiments in the revised version that strengthens the original conclusions in several ways. Newly added efflux assays with ethidium bromide, together with proteinase treatment experiments and ΔacrAΔtolC and ΔqseB/qseC mutant data, illustrate that the low accumulating subpopulation can actively export intracellular compounds. The authors took great care to temper their language to acknowledge other potential alternatives that could explain some of the data such as altered influx, vesicle release or proteolysis, metabolic activity of the cells, indirect effects of sertraline treatment, etc. Additional metabolic dye measurements confirm that low accumulators are less metabolically active, and a new data on nutrient supplementation shows that forcing growth increases peptide uptake and lethality. The authors clarify the crucial point of where antimicrobial peptides actually bind on the cell within the broader survival mechanism and present their conclusions, along with potential caveats, with commendable clarity.

Weaknesses:

Despite these advances, the contribution of efflux may require more direct evidence to further dissect whether efflux is necessary, sufficient, or contributory. The facts that the key low-efflux mutant still retains a small fraction of survivors and that the inhibitors used may cause other physiological changes leading to higher efflux are still unaccounted for. The lipidomic and vesicle findings, while intriguing, remain descriptive, and direct tests of their functional relevance would further solidify the mechanistic models.

Conclusion:

Even with these limitations, the study provides valuable insight into non genetic resistance mechanisms to AMPs and highlights inducible heterogeneity as a critical obstacle to peptide therapeutics. In a much broader context, this study also underscores the importance of efflux physiology even for those antimicrobials that seemingly would not have intracellular targets.

Reviewer #1 (Public review):

Summary:

This is a very well-written paper presenting interesting findings related to the recovery following the end-Permian event in continental settings, from N China. The finding is timely as the topic is actively discussed in the scientific community. The data provides additional insights into the faunal, and partly, floral global recovery following the EPE, adding to the global picture.

Strengths:

The conclusions are supported by an impressive amount of sedimentological and paleontological data (mainly trace fossils) and illustrations.

Reviewer #2 (Public review):

Summary:

The authors made a thorough revision of the manuscript, strengthening the message. They also considered all the comments made by the reviewers and provided appropriate and convincing arguments.

Strengths:

The revised manuscript clarifies all the major points raised by the reviewers, and the way the information is presented (in the text, figures and tables) is clear.

Weaknesses:

The authors provided an appropriate and convincing rebuttal regarding the potential weakness I pointed out in the first review of the manuscript. Therefore, I do not see any major issue in their work.

Reviewer #3 (Public review):

Summary:

The manuscript by Guo and colleagues features the documentation and interpretation of three successions of continental to marginal marine deposits spanning the P/T transition and their respective ichnofaunas. Based on these new data inferences concerning end-Permian mass extinction and Triassic recovery in the tropical realm are discussed.

Strengths:

The manuscript is well written and organized and includes a large amount of new lithological and ichnological data that illuminate ecosystem evolution in a time of large scale transition. The lithological documentations, facies interpretations and ichnotaxonomic assignments look alright (with few exceptions).

Reviewer #1 (Public review):

Summary:

The authors found that IL-1b signaling is pivotal for hypoxemia development and can modulate NETs formation in LPS+HVV ALI model.

Strengths:

They used IL1R1 ko mice and proved that IL1R1 is involved in ALI model proving that IL1b signalling leads towards ARDS. In addition, hypothermia reduces this effect, suggesting a therapeutic option.

Comments on revised version:

The authors have addressed this Reviewer's concerns. The manuscript is much stronger in the current form and can be published.

Reviewer #2 (Public review):

Summary:

The manuscript by Nosaka et al is a comprehensive study exploring the involvement of IL1beta signaling in a 2-hit model of lung injury + ventilation, with a focus on modulation by hypothermia.

Strengths:

The authors demonstrate quite convincingly that interleukin 1 beta plays a role in the development of ventilator-induced lung injury in this model, and that this role includes the regulation of neutrophil extracellular trap formation. The authors use a variety of in vivo animal-based and in vitro cell culture work, and interventions including global gene knockout, cell-targeted knockout and pharmacological inhibition, which greatly strengthen the ability to make clear biological interpretations.

Comments on revised version:

The authors have addressed my concerns/queries.

Reviewer #1 (Public review):

Summary:

The manuscript presents a significant and rigorous investigation into the role of CHMP5 in regulating bone formation and cellular senescence. The study provides compelling evidence that CHMP5 is essential for maintaining endolysosomal function and controlling mitochondrial ROS levels, thereby preventing the senescence of skeletal progenitor cells.

Strengths:

The authors demonstrate that the deletion of Chmp5 results in endolysosomal dysfunction, elevated mitochondrial ROS, and ultimately enhanced bone formation through both autonomous and paracrine mechanisms. The innovative use of senolytic drugs to ameliorate musculoskeletal abnormalities in Chmp5-deficient mice is a novel and critical finding, suggesting potential therapeutic strategies for musculoskeletal disorders linked to endolysosomal dysfunction.

Comments on the latest version:

My concerns were addressed.

Reviewer #3 (Public review):

Summary:

In this study, Zhang et al., reported that CHMP5 restricts bone formation by controlling endolysosome-mitochondrion-mediated cell senescence. Zhang et al., report a novel role of CHMP5 on osteogenesis through affecting cell senescence. Overall, it is an interesting study and provides new insights in the field of cells senescence and bone.

Strengths:

Analyzed the bone phenotype OF CHMP5-periskeletal progenitor-CKO mouse model and found the novel role of senescent cells on osteogenesis and migration.

Weaknesses:

(1) The role and mechanism of CHMP5 gene deletion in enhancing osteogenesis via cellular senescence remain insufficiently elucidated.

(2) The use of the ADTC5 cell line as a skeletal precursor/progenitor model is suboptimal.

Overall, the results support their conclusions.

The impact of this work on the field is its proposal that cellular senescence may exert either inhibitory or promotive effects on osteogenic capacity, depending on cell type and context.

The revised manuscript has addressed most of the concerns raised during the initial review.

Reviewer #1 (Public review):

Summary:

Ferreiro et al. present a method to simulate protein sequence evolution under a birth-death model where sequence evolution is constrained by structural constraints on protein stability. The authors then use this model to explore the predictability of sequence evolution in several viral structural proteins. In principle, this work is of great interest to molecular evolution and phylodynamics, which have struggled to couple non-neutral models of sequence evolution to phylodynamic models like birth-death. Unfortunately, though, the model shows little improvement over neutral models in predicting protein evolution, and this ultimately appears to be due to fundamental conceptual problems with how fitness is modeled and linked to the phylodynamic birth-death model.

Major concerns:

(1) Fitness model: All lineages have the same growth rate r = b-d because the authors assume b+d=1. But under a birth-death model, the growth r is equivalent to fitness, so this is essentially assuming all lineages have the same absolute fitness since increases in reproductive fitness (b) will simply trade off with decreases in survival (d). Thus, even if the SCS model constrains sequence evolution, the birth-death model does not really allow for non-neutral evolution such that mutations can feed back and alter the structure of the phylogeny.

(2) Predictive performance: Similar performance in predicting amino acid frequencies is observed under both the SCS model and the neutral model. I suspect that this rather disappointing result owes to the fact that the absolute fitness of different viral variants could not actually change during the simulations (see comment #1).

(3) Model assessment: It would be interesting to know how much the predictions were informed by the structurally constrained sequence evolution model versus the birth-death model. To explore this, the authors could consider three different models: 1) neutral, 2) SCS, and 3) SCS + BD. Simulations under the SCS model could be performed by simulating molecular evolution along just one hypothetical lineage. Seeing if the SCS + BD model improves over the SCS model alone would be another way of testing whether mutations could actually impact the evolutionary dynamics of lineages in the phylogeny.

(4) Background fitness effects: The model ignores background genetic variation in fitness. I think this is particularly important as the fitness effects of mutations in any one protein may be overshadowed by the fitness effects of mutations elsewhere in the genome. The model also ignores background changes in fitness due to the environment, but I acknowledge that might be beyond the scope of the current work.

(5) In contrast to the model explored here, recent work on multi-type birth-death processes has considered models where lineages have type-specific birth and/or death rates and therefore also type-specific growth rates and fitness (Stadler and Bonhoeffer, 2013; Kunhert et al., 2017; Barido-Sottani, 2023). Rasmussen & Stadler (eLife, 2019) even consider a multi-type birth-death model where the fitness effects of multiple mutations in a protein or viral genome collectively determine the overall fitness of a lineage. The key difference with this work presented here is that these models allow lineages to have different growth rates and fitness, so these models truly allow for non-neutral evolutionary dynamics. It would appear the authors might need to adopt a similar approach to successfully predict protein evolution.

Reviewer #2 (Public review):

Summary:

In this study, "Forecasting protein evolution by integrating birth-death population models with structurally constrained substitution models", David Ferreiro and co-authors present a forward-in-time evolutionary simulation framework that integrates a birth-death population model with a fitness function based on protein folding stability. By incorporating structurally constrained substitution models and estimating fitness from ΔG values using homology-modeled structures, the authors aim to capture biophysically realistic evolutionary dynamics. The approach is implemented in a new version of their open-source software, ProteinEvolver2, and is applied to four viral proteins from HIV-1 and SARS-CoV-2.

Overall, the study presents a compelling rationale for using folding stability as a constraint in evolutionary simulations and offers a novel framework and software to explore such dynamics. While the results are promising, particularly for predicting biophysical properties, the current analysis provides only partial evidence for true evolutionary forecasting, especially at the sequence level. The work offers a meaningful conceptual advance and a useful simulation tool, and sets the stage for more extensive validation in future studies.

Strengths:

The results demonstrate that fitness constraints based on protein stability can prevent the emergence of unrealistic, destabilized variants - a limitation of traditional, neutral substitution models. In particular, the predicted folding stabilities of simulated protein variants closely match those observed in real variants, suggesting that the model captures relevant biophysical constraints.

Weaknesses:

The predictive scope of the method remains limited. While the model effectively preserves folding stability, its ability to forecast specific sequence content is not well supported. Only one dataset (HIV-1 MA) is evaluated for sequence-level divergence using KL divergence; this analysis is absent for the other proteins. The authors use a consensus Omicron sequence as a representative endpoint for SARS-CoV-2, which overlooks the rich longitudinal sequence data available from GISAID. The use of just one consensus from a single time point is not fully justified, given the extensive temporal and geographical sampling available. Extending the analysis to include multiple timepoints, particularly for SARS-CoV-2, would strengthen the predictive claims. Similarly, applying the model to other well-sampled viral proteins, such as those from influenza or RSV, would broaden its relevance and test its generalizability.

It would also be informative to include a retrospective analysis of the evolution of protein stability along known historical trajectories. This would allow the authors to assess whether folding stability is indeed preserved in real-world evolution, as assumed in their model.

Finally, a discussion on the impact of structural templates - and whether the fixed template remains valid across divergent sequences - would be valuable. Addressing the possibility of structural remodeling or template switching during evolution would improve confidence in the model's applicability to more divergent evolutionary scenarios.