Reviewer #3 (Public Review):

Summary:

This research article reports that a greater number of senescent osteoclasts (SnOCs), which produce Netrin-1 and NGF, are responsible for innervation in the LSI and aging animal models.

Strengths:

The research is based on previous findings in the authors' lab and the fact that the IVD structure was restored by treatment with ABT263. The logic is clear and clarifies the pathological role of SnOCs, suggesting the potential utilization of senolytic drugs for the treatment of LBP. Generally, the study is of good quality and the data is convincing.

Weaknesses:

All my concerns have been well addressed, no further comments.

Reviewer #1 (Public Review):

Here, using an organoid system, Wong et al generated a new model of hereditary diffuse leukoencephalopathy with axonal spheroids, with which they investigated how CSF1R-mutaions affect the phenotypes of microglia/macrophages, and revealed metabolic changes in microglia/macrophages associated with a proinflammatory phenotype.

In general, this paper is interesting and well-written, and tackles important issues to be addressed.

This study suffers from several major concerns and limitations that dampen the value of the study. As the authors also mentioned, models that perfectly recapitulate the complexity of the HDLS brain the models would be required to better understand the molecular mechanisms of the disease. In this regard, it is unclear how nicely the organoid system in this study can recapitulate the condition in patients with HDLS (e.g. reduced microglia density, downregulated expression of P2YR12, pathological alterations). In addition, the authors used two different models with distinct mutations that could produce different readouts in CSF1R-mediated cellular responses.

Although the reviewer does understand the importance of providing several options/tools to study rare diseases like HDLS and the difficulty of generating stable organoids with less variation, it is unclear if the different outcomes between HD1 and HD2 are generated through different mutations or simply due to different differentiation efficiency from iMacs (e.g. Figure 2B), which needs to be confirmed. Lastly, there is an over-interpretation regarding the results in Figure 6A. There is no difference between isoHD1 iMac control and HD1 Mut iMac.

Reviewer #2 (Public Review):

Summary:

This paper investigates a rare and severe brain disease called Hereditary Diffuse Leukoencephalopathy with Axonal Spheroids (HDLS). The authors aimed to understand how mutations in the gene CSF-1R affect microglia, the resident immune cells in the brain, and which alterations and factors lead to the specific pathophysiology. To model the human brain with the pathophysiology of HDLS, they used the human-specific model system of induced pluripotent stem cell (iPSC)-derived forebrain organoids with integrated iPSC-derived microglia (iMicro) from patients with the HDLS-causing mutation and an isogenic cell line with the corrected genome. They found that iPSC-derived macrophages (iMac) with HDLS mutations showed changes in their response, including increased inflammation and altered metabolism. Additionally, they studied these iMacs in forebrain organoids, where they differentiate into iMicro, and showed transcriptional differences in isolated iMicro when carrying the HDLS mutation. In addition, the authors described the influence of the mutation within iMicro on the transcriptional level of neurons and neural progenitor cells (NPCs) in the organoid. They observed that the one mutation showed implications for impaired development of neurons, possibly contributing to the progression of the disease.

Overall, this study provides valuable insights into the mechanisms underlying HDLS and emphasizes the importance of studying diseases like these with a suitable model system. These findings, while promising, represent only an initial step towards understanding HDLS and similar neurodegenerative diseases, and thus, their direct translation into new treatment options remains uncertain.

Strengths:

The strength of the work lies in the successful reprogramming of two HDLS patient-derived induced pluripotent stem cells (iPSCs) with different mutations, which is crucial for the study of HDLS using human forebrain organoid models. The use of corrected isogenic iPSC lines as controls increases the validity of the mutation-specific observations. In addition, the model effectively mimics HDLS, particularly concerning deficits in the frontal lobe, mirroring observations in the human brain. Obtaining iPSCs from patients with different CSF1R mutations is particularly valuable given the limitations of rodent and zebrafish models when studying adult-onset neurodegenerative diseases. The study also highlights significant metabolic changes associated with the CSF1R mutation, particularly in the HD2 mutant line, which is confirmed by the HD1 line. In addition, the work shows transcriptional upregulation of the proinflammatory cytokine IL-1beta in cells carrying the mutation, particularly when they phagocytose apoptotic cells, providing further insight into disease mechanisms.

Weaknesses:

The authors have not elucidated the significance of the increased CSF1 dosage in Figure 2F, aside from its effect on cell viability, lacking a thorough discussion of this result. Additionally, while transcriptomic and metabolic alterations related to the mutation were demonstrated in iMac models, similar investigations in iMicros are absent, necessitating further experiments to validate the findings across cell models. The conclusion drawn regarding cytokine levels lacks robust support from the data, particularly considering the varied responses observed in different mutant lines. Further analysis of the secretome (e.g. via ELISA) could provide additional insights. Moreover, the characterization of iMicros is incomplete, with limited protein-level analysis (e.g. validate RNA-seq via flow cytometry). Additionally, the claim of microglial-like morphology lacks adequate evidence, as the provided image is insufficient for such an assessment. RNA-seq experiments should be represented better, it is not possible to read the legends or gene names in the figures. Maybe the data sets can be combined into PCAone and one overall analysis, e.g. via WGCNA-like analyses? This would make it easier for the reader to compare the two cell lines side by side. Furthermore, inaccuracies and omissions in the figure legends compromise the clarity of data representation. Statistical test information is missing. Finally, inconsistent terminology usage throughout the paper may confuse readers (iMac versus iMicros).

Reviewer #2 (Public Review):

Summary:

Fertilization is a crucial event in sexual reproduction, but the molecular mechanisms underlying egg-sperm fusion remain elusive. Elofsson A et al. used AlphaFold to explore possible synapse-like assemblies between sperm and egg membrane proteins during fertilization. Using a systematic search of protein-protein interactions, the authors proposed a pentameric complex of three sperm (IZUMO1, SPACA6, and TMEM81) and two egg (JUNO and CD9) proteins, providing a new structural model to be used in future structure-function studies.

Strengths:

(1) The study uses the AlphaFold algorithm to predict higher-order assemblies. This approach could offer insights into a highly transient protein complex, which are challenging to detect experimentally.<br /> (2) The article predicts a pentameric complex between proteins involved in fertilization, shedding light on the architectural aspects of the egg-sperm fusion synapse.

Weaknesses:

The proposed model, which is a prediction from a modeling algorithm, lacks experimental validation of the identity of the components and the predicted contacts.

It is noteworthy that in an independent study, Deneke et al. provides experimental evidence of the interaction between IZUMO1/SPACA6/TMEM81 in zebrafish. This is an important element that supports the findings presented in this manuscript

Regarding the authors response on the question of a global search:<br /> I understand that a global search might be difficult to interpret because a large number of putative false positives. But it is this type of information that is needed to assess the validity of the model and the scoring power in the absence of any experimental validation. At minimum, the search should include a negative control set of proteins known to be unrelated to sperm fertilization or homologous egg-sperm fusion complexes from incompatible species to account for species-specific interactions.

I acknowledge that experimentally validating highly transient complexes presents technical hurdles. However, a high-confidence structural model could enable the design of point mutations specifically disrupting the predicted interactions. Subsequent rescue experiments could then validate the directionality of these interactions. Ultimately, such experiments are crucial for robust model validation.

Reviewer #3 (Public Review):

Summary:

Sperm-egg fusion is a critical step in successful fertilization. Although several proteins have been identified in mammals that are required for sperm-egg adhesion and fusion, it is still unclear whether there are other proteins involved in this process and how the reported proteins complex co-operate to complete the fusion process. In this study, the authors first identified TMEM81 as a structural homologue of IZUMO1 and SPACA6, and predicted the interactions with a pool of human proteins associated with gamete fusion, using AlphaFold-Multimer, a recent advance in protein complex structure prediction. The prediction is compelling and well discussed, and the experimental evidence to verify this interaction is lacking in this study but supported by a complementary and independent study by another group.

Strengths:

The authors present a pentameric complex formation of four previously reported proteins involved in egg/sperm interaction together with TMEM181 using a deep learning tool, AlphaFold-Multimer.

Weaknesses:

It is intriguing to see that some of the proteins involved in sperm-egg interaction are successfully predicted to be assembled into a single multimeric structure by AlphaFold-Multimer. The experimental validation of the interactions is not directly supported in this study. As there are more candidate proteins in the process, testing other possible protein interactions more comprehensively will provide more rationale for the current 3D multi-protein modeling.

Reviewer #1 (Public Review):

Summary:

In this study, the authors generated a novel transgenic mouse line OpalinP2A-Flpo-T2A-tTA2 to specifically label mature oligodendrocytes, and at the same time their embryonic origins by crossing with a progenitor cre mouse line. With this clever approach, they found that LGE/CGE-derived OLs make minimum contributions to the neocortex, whereas MGE/POA-derived OLs make a small but lasting contribution to the cortex. These findings are contradictory to the current belief that LGE/CGE-derived OPCs make a sustained contribution to cortical OLs, whereas MGE/POA-derived OPCs are completely eliminated. Thus, this study provides a revised and more comprehensive view on the embryonic origins of cortical oligodendrocytes. To specifically label mature oligodendrocytes, and at the same time their embryonic origins by crossing with a progenitor cre mouse line. With this clever approach, they found that LGE/CGE-derived OLs make minimum contributions to the neocortex, whereas MGE/POA-derived OLs make a small-but-lasting contribution to to cortex. These findings are contradictory to the current belief that LGE/CGE-derived OPCs make a sustained contribution to cortical OLs, whereas MGE/POA-derived OPCs are completely eliminated. Thus, this study has provided a revised and updated view on the embryonic origins of cortical oligodendrocytes.

Strengths:

The authors have generated a novel transgenic mouse line to specifically label mature differentiated oligodendrocytes, which is very useful for tracing the final destiny of mature myelinating oligodendrocytes. Also, the authors carefully compared the distribution of three progenitor cre mouse lines and suggested that Gsh-cre also labeled dorsal OLs, contrary to the previous suggestion that it only marks LGE-derived OPCs. In addition, the author also analyzed the relative contributions of OLs derived from three distinct progenitor domains in other forebrain regions (e.g. Pir, ac). Finally, the new transgenic mouse lines and established multiple combinatorial genetic models will facilitate future investigations of the developmental origins of distinct OL populations and their functional and molecular heterogeneity.

Weaknesses:

Since OpalinP2A-Flpo-T2A-tTA2 only labels mature oligodendrocytes but not OPCs, the authors can not suggest that the lack of LGE/CGE-derived-OLs in the neocortex is less likely caused by competitive postnatal elimination, but more likely due to limited production and/or allocation (line 118-9). It remains possible that LGE/CGE-derived OPCs migrate into the cortex but are later eliminated.

Reviewer #2 (Public Review):

Summary:

In this manuscript, Cai et al use a combination of mouse transgenic lines to re-examine the question of the embryonic origin of telencephalic oligodendrocytes (OLs). Their tools include a novel Flp mouse for labelling mature oligodendrocytes and a number of pre-existing lines (some previously generated by the last author in Josh Huang's lab) that allowed combinatorial or subtractive labelling of oligodendrocytes with different origins. The conclusion is that cortically-derived OLs are the predominant OL population in the motor and somatosensory cortex and underlying corpus callosum, while the LGE/CGE generates OLs for the piriform cortex and anterior commissure rather than the cerebral cortex. Small numbers of MGE-derived OLs persist long-term in the motor, somatosensory and piriform cortex.

Strengths:

The strength and novelty of the manuscript lies in the elegant tools generated and used and which have the potential to elegantly and accurately resolve the issue of the contribution of different progenitor zones to telencephalic regions.

Weaknesses:

(1) Throughout the manuscript (with one exception, lines 76-78), the authors quantified OL densities instead of contributions to the total OL population (as a % of ASPA for example). This means that the reader is left with only a rough estimation of the different contributions.

(2) All images and quantifications have been confined to one level of the cortex and the potential of the MGE and the LGE/CGE to produce oligodendrocytes for more anterior and more posterior cortical regions remains unexplored.

(3) Hence, the statement that "In summary, our findings significantly revised the canonical model of forebrain OL origins (Figure 4A) and provided a new and more comprehensive view (Figure 4B )." (lines 111, 112) is not really accurate as the findings are neither new nor comprehensive. Published manuscripts have already shown that (a) cortical OLs are mostly generated from the cortex [Tripathi et al 2011 (https://doi.org/10.1523/JNEUROSCI.6474-10.2011), Winker et al 2018 (https://doi.org/10.1523/JNEUROSCI.3392-17.2018) and Li et al (https://doi.org/10.1101/2023.12.01.569674)] and (b) MGE-derived OLs persist in the cortex [Orduz et al 2019 (https://doi.org/10.1038/s41467-019-11904-4) and Li et al 2024 (https://doi.org/10.1101/2023.12.01.569674)]. Extending the current study to different rostro-caudal regions of the cortex would greatly improve the manuscript.

Reviewer #3 (Public Review):

In the manuscript entitled "Embryonic Origins of Forebrain Oligodendrocytes Revisited by Combinatorial Genetic Fate Mapping," Cai et al. used an intersectional/subtractional strategy to genetically fate-map the oligodendrocyte populations (OLs) generated from medial ganglionic eminence (NKX2.1+), lateral ganglionic eminences, and dorsal progenitor cells (EMX1+). Specifically, they generated an OL-expressing reporter mouse line OpalinP2A-Flpo-T2A-tTA2 and bred with region-specific neural progenitor-expressing Cre lines EMX1-Cre for dOL and NKX2.1-Cre for MPOL. They used a subtractional strategy in the OpalinFlp::Emx1Cre::Nkx2.1Cre::RC::FLTG mouse line to predict the origins of OLs from lateral/caudal ganglionic eminences (LC). With their genetic tools, the authors concluded that neocortical OLs primarily consist of dOLs. Although the populations of OLs (dOLs or MP-OLs) from Emx1+ or Nkx2.1+ progenitors are largely consistent with previous findings, they observed that MP-OLs contribute minimally but persist into adulthood without elimination as in the previous report (PMID: 16388308).

Intriguingly, by using an indirect subtraction approach, they hypothesize that both Emx1-negative and Nkx2.1-negative cells represent the progenitors from lateral/caudal ganglionic eminences (LC), and conclude that neocortical OLs are not derived from the LC region. This is in contrast to the previous observation for the contribution of LC-expressing progenitors (marked by Gsx2-Cre) to neocortical OLs (PMID: 16388308). The authors claim that Gsh2 is not exclusive to progenitor cells in the LC region (PMID: 32234482). However, Gsh2 exhibits high enrichment in the LC during early embryonic development. The presence of a small population of Gsh2-positive cells in the late embryonic cortex could originate/migrate from Gsh2-positive cells in the LC at earlier stages (PMID: 32234482). Consequently, the possibility that cortical OLs derived from Gsh2+ progenitors in LC could not be conclusively ruled out. Notably, a population of OLs migrating from the ventral to the dorsal cortical region was detected after eliminating dorsal progenitor-derived OLs (PMID: 16436615).

The indirect subtraction data for LC progenitors drawn from the OpalinFlp-tdTOM reporter in Emx1-negative and Nkx2.1-negative cells in the OpalinFlp::Emx1Cre::Nkx2.1Cre::RC::FLTG mouse line present some caveats that could influence their conclusion. The extent of activity from the two Cre lines in the OpalinFlp::Emx1Cre::Nkx2.1Cre::RC::FLTG mice remains uncertain. The OpalinFlp-tdTOM expression could occur in the presence of either Emx1Cre or Nkx2.1Cre, raising questions about the contribution of the individual Cre lines. To clarify, the authors should compare the tdTOM expression from each individual Cre line, OpalinFlp::Emx1Cre::RC::FLTG or OpalinFlp::Nkx2.1Cre::RC::FLTG, with the combined OpalinFlp::Emx1Cre::Nkx2.1Cre::RC::FLTG mouse line. This comparison is crucial as the results from the combined Cre lines could appear similar to only one Cre line active.

Overall, the authors provided intriguing findings regarding the origin and fate of oligodendrocytes from different progenitor cells in embryonic brain regions. However, further analysis is necessary to substantiate their conclusion about the fate of LC-derived OLs convincingly.

Reviewer #1 (Public Review):

This is an interesting, informative, and well-designed study that combines theoretical and experimental methodologies to tackle the phenomenon of higher-resolution structures/substructures in model biomolecular condensates. However, there is significant room for improvement in the presentation and interpretation of the results. As it stands, the precise definition of "frustration," which is a main theme of this manuscript (as emphasized in the title), is not sufficiently well articulated. This situation should be rectified to avoid "frustration" becoming a "catch-all" term without a clear perimeter of applicability rather than a precise, informative description of the physical state of affairs. There are also a few other concerns, e.g., regarding interpretation of correlation of phase-separation critical temperature and transfer free energy of amino acid residues as well as the difference between critical temperature and onset temperature, and the way the simulated configurations are similar to that of gyroids. Accordingly, the manuscript should be revised to address the following:

(1) It is accurately pointed out on p.4 that elastin-like polypeptides (ELPs) undergo heat-induced phase separation and therefore exhibit lower critical solution temperatures (LCSTs). But it is not entirely clear how this feature is reproduced by the authors' simulation. A relationship between simulated surface tension and "transition temperature" is provided in Fig.1C; but is the "transition temperature" (authors cited ref.41 by Urry) the same as critical temperature? Apparently, Urry's Tt is "critical onset temperature", the temperature when phase separation happens at a given polymer concentration. This is different from the (global) critical temperature LCST - though the two may be correlated-or not-depending on the shape of the phase boundary. Moreover, is the MOFF coarse-grained forcefield (first step in the multi-scale simulation), by itself, capable of reproducing heat-induced phase separation in a way similar to the forcefield of Dignon et al., ACS Cent Sci 5, 821-230 (2019)? Or, is this temperature-dependent effect appearing only subsequently, after the implementation of the MARTINI and/or all-atom steps? Clarification is needed. To afford a more informative context for the authors' introductory discussion, the aforementioned Dignon et al. work and the review by Cinar et al. [Chem Eur J 25, 13049-13069 (2019)], both touching upon the physical underpinning of the LCST feature of elastin, should also be cited along with refs.41-43.

(2) "Frustration" and "frustrated" are used prominently in the manuscript to characterize certain observed molecular configurations (11 times total, in both the title and in the abstract). Apparently, it is the most significant conceptual pronouncement of this work, hence its precise meaning is of central importance to the authors' thesis. Whereas one should recognize that the theoretical and experimental observations are striking without invocation of the "frustration" terminology, usage of the term can be useful if it offers a unifying conceptual framework. However, as it stands, a clear definition of the term "frustration" is lacking, leaving readers to wonder what molecular configurations are considered "frustrated" and what are not (i.e.,is the claim of observation of frustration falsifiable?). For instance, "frustrated microphase separation" appears in both the title and abstract. A logical question one may ask is: "Are all microphase separations frustrated"? If the answer is in the affirmative, does invocation of the term "frustration" add anything to our physical insight? If the answer is not in the affirmative, then how does one distinguish between microphase separations that are frustrated from those that are not frustrated? Presumably all simulated and experimental molecular configurations in the present study are those of lowest free energy for the given temperature. In other words, they are what they are. In the discussion about frustrated phase separation on p.13, for example, the authors appear to refer to the fact that chain connectivity is preventing hydrophobic residues to come together in a way to achieve the most favorable interactions as if there were no chain connectivity (one may imagine in that case all the hydrophobic residues will form a large cluster without microphase separation). Is this what the authors mean by "frustration"? If that's true, isn't that merely stating the obvious, at least for the observed microphase separation? In general, does "frustration" always mean deviation of actual, physical molecular configurations from certain imagined/hypothetical/reference molecular configurations, and therefore dependent upon the choice of the imagined reference configuration? If this is how the authors apply the term "frustration" in the present work, what is the zero-frustration reference state/configuration for microphase separation? And, similarly, what is the zero-frustration reference state/configuration when frustrated EPS-water interactions are discussed (~p.14-p.15, Fig.5)? How do non-frustrated water-protein interactions look like? Is the classic clathrate-like organization of water hydrogen bonds around small nonpolar solute "frustrated"?

(3) In the discussion about the correlation of various transfer free energy scales for amino acids and Urry's critical onset temperature (ref.41) on p.11 and Fig.4, is there any theoretical relationship to be expected between the interactions among amino acids of ELPs and their critical onset temperatures? While a certain correlation may be intuitively expected if the free energy scale "is working", is there any theoretical insight into the mathematical form of this relationship? A clarifying discussion is needed because it bears logically on whether the observed correlation or lack thereof for different transfer energy scales is a good indication of the adequacy of the energy scales in describing the actual physical interactions at play. This question requires some prior knowledge of the expected mathematical relationship between interaction parameters and onset temperature.

(4) To provide a more comprehensive context for the present study, it is useful to compare the microphase separation seen in the authors' simulation with the micelle-like structures observed in recent simulated condensed/aggregated states of hydrophobic-polar (HP) model sequences in Statt et al., J Chem Phys 152, 075101 (2020) [see esp. Fig.6] and Wessén et al., J Phys Chem B 126, 9222-9245 (2022) [see, e.g., Fig.10].

(5) "Gyroid-like morphology" is mentioned several times in the manuscript (p.4, p.8, p.17, Fig.S3). This is apparently an interesting observation but a clear explanation is lacking. A more detailed and specific discussion, perhaps with additional graphical presentations, should be provided to demonstrate why the simulated condensed-phase ELP configurations are similar to the classical description of gyroid as in, e.g., Terrones & Mackay, Chem Phys Lett 207, 45-50 (1993) and Lambert et al., Phil Trans R Soc A 354, 2009-2023 (1996).

Comments on the revised manuscript:

The authors have adequately addressed my previous concerns.

Reviewer #2 (Public Review):

Summary:

Latham A.P. et al. apply simulations and FLIM to analyse several di-block elastin-like polypetides and connect their sequence to the micro-structure of coacervates resulting from their phase-separation.

Strengths:

Understanding the molecular grammar of phase separating proteins and the connection with mesoscale properties of the coacervates is highly relevant. This work provides insights into micro-structures of coacervates resulting from di-block polypetides.

Weaknesses:

The results apply to a very specific architecture (di-block polypetides) with specific sequences.

Reviewer #1 (Public Review):

Most amino acids are stereoisomers in the L-enantiomer, but natural D-serine has also been detected in mammals and its levels shown to be connected to a number of different pathologies. Here, the authors convincingly show that D-serine is transported in the kidney by the neutral amino acid transporter ASCT2 and as a non-canonical substrate for the sodium-coupled monocarboxylate transporter SMCTs. Although both transport D-serine, this important study further shows in a mouse model for acute kidney injury that ASCT2 has the dominant role.

Strengths:

The paper combines proteomics, animal models, ex vivo transport analyses and in vitro transport assays using purified components. The exhaustive methods employed provide compelling evidence that both transporters can translocate D-serine in the kidney.

Weakness:

In the model for acute kidney injury the SMCTs proteins were not showing a significant change in expression levels and were rather analysed based on other, circumstantial evidence. Although its clear SMCTs can transport D-serine its physiological role is less obvious compared to ASCT2.

Reviewer #2 (Public Review):

Summary:

The manuscript "A multi-hierarchical approach reveals D-1 serine as a hidden substrate of sodium-coupled monocarboxylate transporters" by Wiriyasermkul et al. is a resubmission of a manuscript, which focused first on the proteomic analysis of apical membrane isolated from mouse kidney with early ischemia- reperfusion injury (IRI), a well-known acute kidney injury (AKI) model. In a second part, the transport of D-serine by Asct2, Smct1, and Smct2 has been characterized in detail in different model systems, such as transfected cells and proteoliposomes.

Strengths:

A major problem with the first submission was the explanation of the link between the two parts of the manuscript: it was not very clear why the focus on Asct2, Smct1 and Smct2 was a consequence of the proteomic analysis. In the present version of the manuscript, the authors have focused on the expression of membrane transporters in the proteome analysis, thus making the reason for studying Asct2, Smct1 and Smct2 transporters more clear. In addition, the authors used 2D-HPLC to measure plasma and urinary enantiomers of 20 amino acids in plasma and urine samples from sham and ischaemia-reperfusion injury (IRI) mice. The results of this analysis demonstrated the value of D-serine as a potential marker of renal injury. These changes have greatly improved the manuscript and made it more convincing.

Weaknesses:

More than weakness I would speak of discussion points: I have a few suggestions that may help to make the paper more accessible to a general audience.<br /> (1) In the Introduction, when the authors introduce the term "micromolecules", it would be beneficial to provide a precise definition or clarification of what they mean by this term. Adding a brief explanation may help the reader to better understand the context.<br /> (2) In line 91, I suggest specifying that this is a renal IRI model.<br /> (3) Lines 167-168 state that Asct2 is localised to the apical side of the renal proximal tubules. Is there any expression of Asct2 in other nephron segments?<br /> (4) Lines 225-226: Have the authors expressed the candidate genes in HEK293 cells with ASCT2 knockdown?<br /> (5) lines 254-255: why was D-serine transport enhanced by ASCT2 knockdown in FlpInTR-SMCT1 or 2 cells?<br /> (6) line 265: The low affinity of SMCT1 for D-serine alone makes it an unlikely transporter for urinary D-serine.<br /> (7) line 316: The authors state that there is a high tubular D-serine reabsorption in IRI and in line 424 that there is an inactivation of DAAO during the pathology. This suggests that there is a reabsorption of D-serine mediated by a transport system in the basolateral membrane domain of proximal tubular cells. Do the authors have any information about this transporter?<br /> (8) in lines 462-463, the authors state: "It is suggested that PAT1 is less active at the apical membrane where the luminal pH is neutral". However, the pH of urine in the proximal tubules is normally acidic due to the high activity of NH3. I suggest rewording this sentence.

Reviewer #3 (Public Review):

Summary:

The main objective of this work has been to delve into the mechanisms underlying the increment of D-serine in serum, as a marker of renal injury.

Strengths:

With a multi-hierarchical approach, the work shows that Ischemia reperfusion injury in kidney causes a specific increment in renal reabsorption of D-serine that, at least in part, is due to the increased expression of the apical transporter ASCT2. In the way, the authors revealed that SMCT1 also transports D-serine.

The manuscript also supports that increased expression of ASCT2, even together with the parallel decreased expression of SMCT1, in renal proximal tubules underlies the increased reabsorption of D-serine responsible of the increment of this enantiomer in serum in a murine model of ischemia reperfusion injury.

Weaknesses:

Remains to be clarified whether ASCT2 has substantial stereospecificity in favor of D- versus L-serine to sustain a ~10-fold decreased in the ratio D-serine/L-serine in the urine of mouse under ischemia reperfusion injury (IRI).<br /> It is not clear how the increment in the expression of ASCT2, in parallel with the decreased expression of SMCT1, results in increased renal reabsorption of D-serine in IRI.

I am satisfied with the changes the authors have introduced in the text of the revised version of their manuscript.

Reviewer #1 (Public Review):

Summary:

LRRK2 protein is familially linked to Parkinson's disease by the presence of several gene variants that all confer a gain-of-function effect on LRRK2 kinase activity.

The authors examine the effects of BDNF stimulation in immortalized neuron-like cells, cultured mouse primary neurons, hIPSC-derived neurons, and synaptosome preparations from the brain. They examine an LRRK2 regulatory phosphorylation residue, LRRK2 binding relationships, and measures of synaptic structure and function.

Strengths:

The study addresses an important research question: how does a PD-linked protein interact with other proteins, and contribute to responses to a well-characterized neuronal signalling pathway involved in the regulation of synaptic function and cell health?

They employ a range of good models and techniques to fairly convincingly demonstrate that BDNF stimulation alters LRRK2 phosphorylation and binding to many proteins. Some effects of BDNF stimulation appear impaired in (some of the) LRRK2 knock-out scenarios (but not all). A phosphoproteomic analysis of PD mutant Knock-in mouse brain synaptosomes is included.

Weaknesses:

The data sets are disjointed, conclusions are sweeping, and not always in line with what the data is showing. Validation of 'omics' data is very light. Some inconsistencies with the major conclusions are ignored. Several of the assays employed (western blotting especially) are likely underpowered, findings key to their interpretation are addressed in only one or other of the several models employed, and supporting observations are lacking.

As examples to aid reader interpretation:

(a) pS935 LRRK2 seems to go up at 5 minutes but goes down below pre-stimulation levels after (at times when BDNF-induced phosphorylation of other known targets remains very high). This is ignored in favour of discussion/investigation of initial increases, and the fact that BDNF does many things (which might indirectly contribute to initial but unsustained changes to pLRRK2) is not addressed.

(b) Drebrin coIP itself looks like a very strong result, as does the increase after BDNF, but this was only demonstrated with a GFP over-expression construct despite several mouse and neuron models being employed elsewhere and available for copIP of endogenous LRRK2. Also, the coIP is only demonstrated in one direction. Similarly, the decrease in drebrin levels in mice is not assessed in the other model systems, coIP wasn't done, and mRNA transcripts are not quantified (even though others were). Drebrin phosphorylation state is not examined.

(c) The large differences in the CRISPR KO cells in terms of BDNF responses are not seen in the primary neurons of KO mice, suggesting that other differences between the two might be responsible, rather than the lack of LRRK2 protein.

(d) No validation of hits in the G2019S mutant phosphoproteomics, and no other assays related to the rest of the paper/conclusions. Drebrin phosphorylation is different but unvalidated, or related to previous data sets beyond some discussion. The fact that LRRK2 binding occurs, and increases with BDNF stimulation, should be compared to its phosphorylation status and the effects of the G2019S mutation.

Reviewer #2 (Public Review):

Taken as a whole, the data in the manuscript show that BDNF can regulate PD-associated kinase LRRK2 and that LRRK2 modifies the BDNF response. The chief strength is that the data provide a potential focal point for multiple observations across many labs. Since LRRK2 has emerged as a protein that is likely to be part of the pathology in both sporadic and LRRK2 PD, the findings will be of broad interest. At the same time, the data used to imply a causal throughline from BDNF to LRRK2 to synaptic function and actin cytoskeleton (as in the title) are mostly correlative and the presentation often extends beyond the data. This introduces unnecessary confusion. There are also many methodological details that are lacking or difficult to find. These issues can be addressed.

(1) The writing/interpretation gets ahead of the data in places and this was confusing. For example, the abstract highlights prior work showing that Ser935 LRRK2 phosphorylation changes LRRK2 localization, and Figure 1 shows that BDNF rapidly increases LRRK2 phosphorylation at this site. Subsequent figures highlight effects at synapses or with synaptic proteins. So is the assumption that LRRK2 is recruited to (or away from) synapses in response to BDNF? Figure 2H shows that LRRK2-drebrin interactions are enhanced in response to BDNF in retinoic acid-treated SH-SY5Y cells, but are synapses generated in these preps? How similar are these preps to the mouse and human cortical or mouse striatal neurons discussed in other parts of the paper (would it be anticipated that BDNF act similarly?) and how valid are SH-SY5Y cells as a model for identifying synaptic proteins? Is drebrin localization to synapses (or its presence in synaptosomes) modified by BDNF treatment +/- LRRK2? Or do LRRK2 levels in synaptosomes change in response to BDNF? The presentation requires re-writing to stay within the constraints of the data or additional data should be added to more completely back up the logic.

(2) The experiments make use of multiple different kinds of preps. This makes it difficult at times to follow and interpret some of the experiments, and it would be of great benefit to more assertively insert "mouse" or "human" and cell type (cortical, glutamatergic, striatal, gabaergic) etc.

(3) Although BDNF induces quantitatively lower levels of ERK or Akt phosphorylation in LRRK2KO preps based on the graphs (Figure 4B, D), the western blot data in Figure 4C make clear that BDNF does not need LRRK2 to mediate either ERK or Akt activation in mouse cortical neurons and in 4A, ERK in SH-SY5Y cells. The presentation of the data in the results (and echoed in the discussion) writes of a "remarkably weaker response". The data in the blots demand more nuance. It seems that LRRK2 may potentiate a response to BDNF that in neurons is independent of LRRK2 kinase activity (as noted). This is more of a point of interpretation, but the words do not match the images.

(4) Figure 4F/G shows an increase in PSD95 puncta per unit length in response to BDNF in mouse cortical neurons. The data do not show spine induction/dendritic spine density/or spine morphogenesis as suggested in the accompanying text (page 8). Since the neurons are filled/express gfp, spine density could be added or spines having PSD95 puncta. However, the data as reported would be expected to reflect spine and shaft PSDs and could also include some nonsynaptic sites.

(5) Experimental details are missing that are needed to fully interpret the data. There are no electron microscopy methods outside of the figure legend. And for this and most other microscopy-based data, there are few to no descriptions of what cells/sites were sampled, how many sites were sampled, and how regions/cells were chosen. For some experiments (like Figure 5D), some detail is provided in the legend (20 segments from each mouse), but it is not clear how many neurons this represents, where in the striatum these neurons reside, etc. For confocal z-stacks, how thick are the optical sections and how thick is the stack? The methods suggest that data were analyzed as collapsed projections, but they cite Imaris, which usually uses volumes, so this is confusing. The guide (sgRNA) sequences that were used should be included. There is no mention of sex as a biological variable.

(6) For Figures 1F, G, and E, how many experimental replicates are represented by blots that are shown? Graphs/statistics could be added to the supplement. For 1C and 1I, the ANOVA p-value should be added in the legend (in addition to the post hoc value provided).

(7) Why choose 15 minutes of BDNF exposure for the mass spec experiments when the kinetics in Figure 1 show a peak at 5 mins?

(8) The schematic in Figure 6A suggests that iPSCs were plated, differentiated, and cultured until about day 70 when they were used for recordings. But the methods suggest they were differentiated and then cryopreserved at day 30, and then replated and cultured for 40 more days. Please clarify if day 70 reflects time after re-plating (30+70) or total time in culture (70). If the latter, please add some notes about re-differentiation, etc.

(9) When Figures 6B and 6C are compared it appears that mEPSC frequency may increase earlier in the LRRK2KO preps than in the WT preps since the values appear to be similar to WT + BDNF. In this light, BDNF treatment may have reached a ceiling in the LRRK2KO neurons.

(10) Schematic data in Figures 5A and C and Figures 5B and E are too small to read/see the data.

Reviewer #1 (Public Review):

Amason et al. investigated the formation of granulomas in response to Chromobacterium violaceum infection, aiming to uncover the cellular mechanisms governing the granuloma response. They identify spatiotemporal gene expression of chemokines and receptors associated with the formation and clearance of granulomas, with a specific focus on those involved in immune trafficking. By analyzing the presence or absence of chemokine/receptor RNA expression, they infer the importance of immune cells in resolving infection. Despite observing increased expression of neutrophil-recruiting chemokines, treatment with reparixin (an inhibitor of CXCR1 and CXCR2) did not inhibit neutrophil recruitment during infection. Focusing on monocyte trafficking, they found that CCR2 knockout mice infected with C. violaceum were unable to form granulomas, ultimately succumbing to infection.

The spatial transcriptomics data presented in the figures could be considered a valuable resource if shared, with the potential for improved and clarified analyses. The primary conclusion of the paper, that C. violaceum infection in the liver cannot be contained without macrophages, would benefit from clarification.

While the spatial transcriptomic data generated in the figures are interesting and valuable, they could benefit from additional information. The manual selection of regions of granulomas for analysis could use additional context - was the rest of the liver not sequenced, or excluded for other reasons? Including a healthy liver in the analysis could serve as a control for any lasting effects at the final time point of 21 days. Providing more context for the scalebars throughout the spatial analyses, such as whether the data are raw counts or normalized based on the number of reads per spatial spot, would be helpful for interpretation, as changes in expression could signal changes in the numbers of cells or changes in the gene expression of cells.

In Figure 4, qualitative measurements are valuable, but having an idea of the raw data for a few of the pursued chemokines/receptors would aid interpretation. It would also be beneficial to clarify whether the reported values are across all clusters and consider focusing on clusters with the greatest change in expression. Figures 5E and F would benefit from clarification regarding the x-axis units and whether the expression levels are summed across all clusters for each time point. Additionally, information on the sequencing depth of the samples would be helpful, particularly as shallow sequencing of RNA can result in poor capture of low-expression transcripts.

Regarding the conclusion of the essentiality of macrophages in granuloma formation, it may be prudent to further investigate the role of macrophages versus CCR2. Analyzing total cell counts in the liver after infection could provide insight into whether the decrease in the fraction of macrophages is due to decreased numbers or infiltration of other cell types. Consideration of experiments deleting macrophages directly, instead of CCR2, could provide more definitive evidence of the necessity of macrophage migration in containing infections.

Reviewer #2 (Public Review):

Summary:

In this study, Amason et al employ spatial transcriptomics and intervention studies to probe the spatial and temporal dynamics of chemokines and their receptors and their influence on cellular dynamics in C. violaceum granulomas. As a result of their spatial transcriptomic analysis, the authors narrow in on the contribution of neutrophil- and monocyte-recruiting pathways to host response. This results in the observation that monocyte recruitment is critical for granuloma formation and infection control, while neutrophil recruitment via CXCR2 may be dispensable.

Strengths:

Since C. violaceum is a self-limiting granulomatous infection, it makes an excellent case study for 'successful' granulomatous inflammation. This stands in contrast to chronic, unproductive granulomas that can occur during M. tuberculosis infection, sarcoidosis, and other granulomatous conditions, infectious or otherwise. Given the short duration of C. violaceum infection, this study specifically highlights the importance of innate immune responses in granulomas.

Another strength of this study is the temporal analysis. This proves to be important when considering the spatial distribution and timing of cellular recruitment. For example, the authors observe that the intensity and distribution of neutrophil- and monocyte-recruiting chemokines vary substantially across infection time and correlate well with their previous study of cellular dynamics in C. violaceum granulomas.

The intervention studies done in the last part of the paper bolster the relevance of the authors' focus on chemokines. The authors provide important negative data demonstrating the null effect of CXCR1/2 inhibition on neutrophil recruitment during C. violaceum infection. That said, the authors' difficulty with solubilizing reparixin in PBS is an important technical consideration given the negative result. On the other hand, monocyte recruitment via CCR2 proves to be indispensable for granuloma formation and infection control. I would hesitate to agree with the authors' interpretation that their data proves macrophages are serving as a physical barrier from the uninvolved liver. It is possible and likely that they are contributing to bacterial control through direct immunological activity and not simply as a structural barrier.

Weaknesses:

There are several shortcomings that limit the impact of this study. The first is that the cohort size is very limited. While the transcriptomic data is rich, the authors analyze just one tissue from one animal per time point. This assumes that the selected individual will have a representative lesion and prevents any analysis of inter-individual variability. Granulomas in other infectious diseases, such as schistosomiasis and tuberculosis, are very heterogeneous, both between and within individuals. It will be difficult to assert how broadly generalizable the transcriptomic features are to other C. violaceum granulomas. Furthermore, this undermines any opportunity for statistical testing of features between time points, limiting the potential value of the temporal data.

Another caveat to these data is the limited or incompletely informative data analysis. The authors use Visium in a more targeted manner to interrogate certain chemokines and cytokines. While this is a great biological avenue, it would be beneficial to see more general analyses considering Visum captures the entire transcriptome. Some important questions that are left unanswered from this study are:

What major genes defined each spatial cluster?

What were the top differentially expressed genes across time points of infection?

Did the authors choose to focus on chemokines/receptors purely from a hypothesis perspective or did chemokines represent a major signature in the transcriptomic differences across time points?

In addition to the absence of deep characterization of the spatial transcriptomic data, the study lacks sufficient quantitative analysis to back up the authors' qualitative assessments. Furthermore, the authors are underutilizing the spatial information provided by Visium with no spatial analysis conducted to quantify the patterning of expression patterns or spatial correlation between factors.

Impact:

The author's analysis helps highlight the chemokine profiles of protective, yet host protective granulomas. As the authors comment on in their discussion, these findings have important similarities and differences with other notable granulomatous conditions, such as tuberculosis. Beyond the relevance to C. violaceum infection, these data can help inform studies of other types of granulomas and hone candidate strategies for host-directed therapy strategies.

Reviewer #1 (Public Review):

Pineda et al investigate the association of the hypothesis that Dux4, an embryonic transcription factor, expression in tumor cells is associated with immune evasion and resistance to immunotherapy. They analyze existing cohorts of bulk RNAseq sequenced tumors across cancer types to identify Dux4 expression and association with survival. They find that Dux4 expression is detected in a higher proportion of metastatic tumors compared to primary tumors, is associated with decreased immune infiltrate and a variety of immune metrics and previously nominated immune signatures, and do an in depth evaluation of a cohort of metastatic urothelial cell carcinoma, finding that Dux4 expression is associated with a more immunodeficient tumor microenvironment (desert or excluded microenvironment) and worse survival in this aPDL1 treated cohort. They then find that Dux4 expression is a major independent predictor of survival in this cohort using different types of survival analyses (KM, Cox PH, and random survival forests). With prior existing biological data supporting the hypothesis (in prior work, the senior author has demonstrated Dux4 expression causally suppresses MHC-I expression in interferon-gamma treated cell lines), the current work links Dux4 expression with less immune activity in clinical tumor samples and with survival in ICI treated urothelial carcinomas, and demonstrates that Dux4 expression provides independent information towards survival including other molecular and clinical characteristics (TMB, ECOG PS as the other strongest markers), and provides interesting resolution on landmark analyses with TMB and Dux4 expression providing greater informativeness at later survival landmarks (e.g. 1 year and later), while ECOG PS has strong informativeness already at earlier time points. This work provides impetus towards more mechanistic and functional dissection of the mechanism of Dux4-associated changes with the tumor microenvironment (e.g. in vivo mouse studies) as well as potential interventional studies (e.g. Dux4 as a target in combination therapies). What the work does not provide is additional resolution on the mechanism of how Dux4 may be associated with a more immunodeficient microenvironment.

The conclusions are generally well supported, but there are issues that would benefit from clarification and extension:

- The finding that Dux4 expression is detected in a higher proportion of metastatic tumors and at higher levels compared to TCGA samples (Fig 1BC) is striking. However, at least for one tumor type (melanoma), the TCGA cohort is comprised of mostly locoregional metastatic (n=81 primary and 367 metastatic tumors in the PanCan Atlas). Since there are annotations for primary and (locoregional) metastatic samples in TCGA, an analysis of the primary vs. locoregional metastasis vs distant metastatic samples seems reasonable and likely informative. The analysis of tumors with matched FFPE and flash frozen samples with hybrid probe capture and polyA sequencing, respectively is a nice validation to show that the difference in Dux4 expression is not due to differences in preservation of starting material/sequencing in the metastatic samples vs TCGA samples (S1BC).<br /> - The findings that Dux4 expression in the metastatic urothelial carcinoma setting is associated with a more immunodeficient microenvironment (Figure 2) is clear and unambiguous using multiple lines of data and analyses (bulk RNAseq, DUX4-positive vs DUX4-negative tumors, different immune cell and cytokine signatures; IHC showing an association with immune deserts and immune excluded phenotypes). However, this is an association and does not demonstrate causality.<br /> - The survival analyses (Fig 3,4,5) show fairly convincingly that Dux4 provide independent predictive information beyond clinical variables and TMB towards survival in the aPDL1 treated metastatic urothelial carcinoma cohort. However, the choice to split the cohort into Dux4 negative (defined as < 0.25 TPM) and Dux4 positive (> 1 TPM) while excluding a large number of patients (n=126 pts) that fall in between has significant impact on the rigor of conclusions. This would benefit from showing all the data (e.g. including the 3rd group of in-betweens in the survival analyses as a separate group).<br /> - The authors demonstrate that adding Dux4 to clinical markers and TMB results in an improved predictive model for survival, but there are a few questions regarding this model as a clinical biomarker<br /> o Is Dux4 expression better than other correlated immune signatures/markers (e.g. interferon gamma, T effector signature, overall immune infiltrate) in providing additional information?<br /> - The use of random survival forests to quantify the (predictive) marginal effect of Dux4+ vs Dux4- expression on survival in a non-parametric model as well as shed light on association with survival at different landmark times using Shapley values is quite interesting and well conducted.

Reviewer #2 (Public Review):

Summary:

This article takes an expansive look at the potential role of DUX4 in cancer treatment and prognosis, including its correlation with other key biomarkers, the potential for cancer to be resistant to treatment, and risk prediction.

Strengths:

The primary strength of this work is the breadth of the analyses. The authors have linked DUX4 to not just one but multiple points in the trajectory of cancer, which increases the face validity of their conclusion that DUX4 is meaningfully related to the course of a cancer as well as the prognosis for a patient.

Statistically, the authors have taken care to properly validate their findings using appropriate bootstrapping and testing strategies.

Weaknesses:

Several weaknesses are noted. First, there is little-to-no description of the underlying sample population. It is only stated that "several large cohorts of patients with different metastatic cancers" were analyzed, and that a cohort of patients with advanced urothelial cancer was used for estimating associations with clinical outcomes. Lacking is information on the sampling mechanism, inclusion/exclusion criteria, treatment modalities, the definition of 'time = 0', the number of events observed, or even the sample size. Knowledge about the underlying study design would help explain some counterintuitive results, e.g. that the hazard of death among patients with Stage IV cancer is half that of those with Stage I cancer (Table 1); presumably this is not because Stage IV is actually protective but rather an artifact of the sampling scheme for these data. Second, the definition of negative versus positive DUX4 expression varies throughout the paper. In Figure 2B, Figure 3A, and Figure 3C, it is defined as >1 TPM vs. <= 1 TPM; in Figure 4A and Figure 5A, it is defined as >1 TPM vs. < 0.25 TPM; in Figure S1C it is partitioned into four groups, with boundaries defined at 0.25 TPM, 1 TPM, and 5 TPM. If categorization is needed, a rationale should be provided (ideally prospectively and not based upon the observed data, so as to avoid the perception of forking paths analyses), and it should be consistently applied. Third and finally, data seem to be occasionally excluded without rationale. For example, as mentioned above, the Cox model presented in Figure 4A seems to exclude all patients with DUX4 TPM between 0.25 and 1. Figure 3C excludes patients with TMB in the lowest quartile (although the decision was ostensibly to control for TMB confounding, there are more appropriate ways to do so that don't result in loss of data, e.g. a stratified KM plot). Excluding patients based upon a particular region of the covariate space makes interpreting the resulting model awkward.

Reviewer #1 (Public Review):

The authors have made a novel and important effort to distinguish and include different sources of active deformations for fitting C elegans embryo development: cyclic muscle contractions and actomyosion circumferential stresses. The combination and synchronisation of both contributions are, according to the model, responsible for different elongation rates, and can induce bending and torsion deformations, which are a priori not expected from purely contractile forces. The model can be applied to other growth processes in initially cylindrical shapes.

The tilt of the fibers is an important assumption of the model. However, fiber direction in Figure 3B is not fully clear for explaining the tilting. The fiber in 3B has not very much in common with the fibers in the color part of the figure. Also, is vector m supposed to be tangent to the fiber? In the figure does not seem to be so. It should be expected that alpha is a consequence of the deformation, not as an input parameter, as it seems in the tests of Figure 6A. How is the value of alpha chosen? According to Figure 6, torsion is expected for alpha>0, but for beta=pi/2 and alpha>0 no torsion may be obtained. In fact, it seems that torsion should appear when cos(beta)*sin(alpha)>0. As a consequence, value of beta should be given in Figure 6. Can the amount of torsion be tested as a function of alpha and beta?

The transfer of energy and deformation is a very interesting aspect of the paper, and also crucial for the model and predicting elongation. However, the modelling of this transfer remains very obscure and only explained in the Appendix. Some more details on how the transfer is selected should be given in the main text. Can the transfer of energy interpreted as a change of the relaxed reference configuration? Once a ratio of the energy transferred is fixed, the assumption on elongation distribution should be stated. (Uniformly? ) The authors should also define in the main text the factor g_a1, and explain how this value is computed from condition W_c=W_r .

Given the convoluted shape of the embryo in the egg, contact may be a crucial mechanism for determining growth and torsion. The model does not include this contact, and this limitation should be reflected in the article.

Minor comment:<br /> -Line 300: "we determine the optimal values for the activation parameters". the optimal with respect to which objective? Norm of difference between experimental and computational displacements? How this is quantified needs to be specified.

Reviewer #2 (Public Review):

Summary:

During C. elegans development, embryos undergo elongation of their body axis in absence of cell proliferation or growth. This process relies in an essential way on periodic contractions of two pairs muscles that extend along the embryo's main axis. How contraction can lead to extension along the same direction is unknown.

To address this question, the authors use a continuum description of a multicomponent elastic solid. The various components are the interior of the animal, the muscles, and the epidermis. The different components form separate compartments and are described as hyperelastic solids with different shear moduli. For simplicity, a cylindrical geometry is adopted. The authors consider first the early elongation phase, which is driven by contraction of the epidermis, and then late elongation, where contraction of the muscles injects elastic energy into the system, which is then transferred into elongation. The authors get elongation that can be successfully fitted to the elongation dynamics of wild type worms and two mutant strains.

Strengths:

The work proposes a physical mechanism underlying a puzzling biological phenomenon. The framework developed by the authors could be used to explain phenomena in other organisms and could be exploited in the design of soft robots.

Weaknesses:

(1) The manuscript is hard to read without being very familiar with continuum descriptions of elastic media. This might make the work difficult to access for biologists. This is a real pity because the findings are potentially of great interest to developmental biologists and engineers alike.

(2) The discussion of the worm's mechanical properties could go deeper. The authors hardly justify their assumptions.

Reviewer #1 (Public Review):

Summary:

This study presents careful biochemical experiments to understand the relationship between LRRK2 GTP hydrolysis parameters and LRRK2 kinase activity. The authors report that incubation of LRRK2 with ATP increases the KM for GTP and decreases the kcat. From this they suppose an autophosphorylation process is responsible for enzyme inhibition. LRRK2 T1343A showed no change, consistent with it needing to be phosphorylated to explain the changes in G-domain properties. The authors propose that phosphorylation of T1343 inhibits kinase activity and influences monomer-dimer transitions.

Strengths:

Strengths of the work are the very careful biochemical analyses and interesting result for wild type LRRK2.

Weaknesses:

The conclusions related to involvement of a monomer-dimer transition are to this reviewer, premature and an independent method needs to be utilized to bolster this aspect of the story.

Reviewer #2 (Public Review):

As discussed in the original review, this manuscript is an important contribution to a mechanistic understanding of LRRK2 kinase. Kinetic parameters for the GTPase activity of the ROC domain have been determined in the absence/presence of kinase activity. A feedback mechanism from the kinase domain to GTP/GDP hydrolysis by the ROC domain is convincingly demonstrated through these kinetic analyses. However, a regulatory mechanism directly linking the T1343 phospho-site and a monomer/dimer equilibrium is not fully supported. The T1343A mutant has reduced catalytic activity and can form similar levels of dimer as WT. The revised manuscript does point out that other regulatory mechanisms can also play a role in kinase activity and GTP/GDP hydrolysis (Discussion section). The environmental context in cells cannot be captured from the kinetic assays performed in this manuscript, and the introduction contains some citations regarding these regulatory factors. This is not a criticism, the detailed kinetics here are rigorous, but it is simply a limitation of the approach. Caveats concerning effects of membrane localization, Rab/14-3-3 proteins, WD40 domain oligomers, etc... should be given more prominence than a brief (and vague) allusion to 'allosteric targeting' near the end of the Discussion.

Specific comments

(1) The revised version is better organized with respect to the significance of monomer/dimer equilibrium and the relevance of the GTP-binding region of ROC domain that encompasses the T1343 phospho-site. The relevance of monomers/dimers of LRRK2 from previous studies is better articulated and readers are able to follow the reasoning for the various mutations.

(2) As a suggestion I would change the following on page 6 to clarify for readers:<br /> "...would show no change in kcat and KM values upon in vitro ATP treatment" to:<br /> "...would show no change in kcat and KM values for GTP hydrolysis upon in vitro ATP treatment"

(3) The levels of dimer in WT (+ATP) and T1343A (+/- ATP) are the same, about 40-45%. These data are cited when the authors state that ATP-induced monomerization is 'abolished' (page 6). My suggestion is to re-phrase this conclusion for consistency with data (Fig 5). For example, one can state that 'ATP incubation does not affect the percentage of dimer for the T1343A variant of LRRK2'. This would be similar to the authors' description of these data on page 8 - 'no difference in dimer formation upon ATP treatment'.

Reviewer #1 (Public Review):

The author studies a family of models for heritable epigenetic information, with a focus on enumerating and classifying different possible architectures. The key aspects of the paper are:

- Enumerate all 'heritable' architectures for up-to 4 constituents.<br /> - A study of whether permanent ("genetic") or transient ("epigenetic") perturbations lead to heritable changes<br /> - Enumerated the connectivity of the "sequence space" formed by these heritable architectures<br /> - Incorporating stochasticity, the authors explore stability to noise (transient perturbations)<br /> - A connection is made with experimental results on C elegans.

The study is timely, as there is a renewed interest in the last decade in non-genetic, heritable heterogeneity (e.g., from single-cell transcriptomics). Consequently, there is a need for a theoretical understanding of the constraints on such systems. There are some excellent aspects of this study: for instance, the attention paid to how one architecture "mutates" into another. Unfortunately, the manuscript as a whole does not succeed in formalising nor addressing any particular open questions in the field. Aside from issues in presentation and modelling choices (detailed below), it would benefit greatly from a more systematic approach rather than the vignettes presented.

## Terminology

The author introduces a terminology for networks of interacting species in terms of "entities" and "sensors" -- the former being nodes of a graph, and the latter being those nodes that receive inputs from other nodes. In the language of directed graphs, "entities" would seem to correspond to vertices, and "sensors" those vertices with positive indegree and outdegree. Unfortunately, the added benefit of redefining accepted terminology from the study of graphs and networks is not clear.

## Model

The model seems to suddenly change from Figure 4 onwards. While the results presented here have at least some attempt at classification or statistical rigour (i.e. Fig 4 D), there are suddenly three values associated with each entity ("property step, active fraction, and number"). Furthermore, the system suddenly appears to be stochastic. The reader is left unsure what has happened, especially after having made the effort to deduce the model as it was in Figs 1 through 3. No respite is to be found in the SI, either, where this new stochastic model should have been described in sufficient detail to allow one to reproduce the simulation.

## Perturbations

Inspired especially by experimental manipulations such as RNAi or mutagenesis, the author studies whether such perturbations can lead to a heritable change in network output. While this is naturally the case for permanent changes (such as mutagenesis), the author gives convincing examples of cases in which transient perturbations lead to heritable changes. Presumably, this is due the the underlying multistability of many networks, in which a perturbation can pop the system from one attractor to another.

Unfortunately, there appears to be no attempt at a systematic study of outcomes, nor a classification of when a particular behaviour is to be expected. Instead, there is a long and difficult-to-read description of numerical results that appear to have been sampled at random (in terms of both the architecture and parameter regime chosen). The main result here appears to be that "genetic" (permanent) and "epigenetic" (transient) perturbations can differ from each other -- and that architectures that share a response to genetic perturbation need not behave the same under an epigenetic one. This is neither surprising (in which case even illustrative evidence would have sufficed) nor is it explored with statistical or combinatorial rigour (e.g. how easy is it to mistake one architecture for another? What fraction share a response to a particular perturbation?)

As an additional comment, many of the results here are presented as depending on the topology of the network. However, each network is specified by many kinetic constants, and there is no attempt to consider the robustness of results to changes in parameters.

## DNA analogy

At two points, the author makes a comparison between genetic information (i.e. DNA) and epigenetic information as determined by these heritable regulatory architectures. The two claims the author makes are that (i) heritable architectures are capable of transmitting "more heritable information" than genetic sequences, and (ii) that, unlike DNA, the connectivity (in the sense of mutations) between heritable architectures is sparse and uneven (i.e. some architectures are better connected than others).

In both cases, the claim is somewhat tenuous -- in essence, it seems an unfair comparison to consider the basic epigenetic unit to be an "entity" (e.g., an entire transcription factor gene product, or an organelle), while the basic genetic unit is taken to be a single base-pair. The situation is somewhat different if the relevant comparison was the typical size of a gene (e.g., 1 kb).

Reviewer #1 (Public Review):

Summary:

This describes the molecular identity of the intermediate status of cranial neural crest cells (NCCs) during the initial delamination process. Taking advantage of single-cell RNA seq, the authors identify new populations of cells during EMT characterized by a specific set of gene expressions, including Dlc1. Promigratory cranial NCCs differentiate through different trajectories depending on their cell cycle phases but converge into a common progenitor, then differentiate into mesenchymal cells expressing region-specific genes.

Strengths:

Single-cell RNA seq data convincingly support what the authors claim. This is the first time to identify intermediate states between premigratory and migratory cranial NCCs. Silencing one of the marker genes, Dlc1, reduces the migratory activity of cranial NCCs. These findings deepen our understanding of the mechanism of EMT in general.

Comments on revised version:

Weaknesses:

None after substantial revision.

Reviewer #2 (Public Review):

Zhao et al., focus on mechanisms through which cells convert from epithelium to mesenchyme and become migratory. This phenomenon of epithelial-to-mesenchymal transition (EMT) occurs during both embryonic development and cancer progression. During cancer progression, EMT seemingly includes cells at intermediate states as defined by the combinatorial expression of epithelial and mesenchymal markers. But the importance of these markers and the role of these intermediate states remains unclear. Moreover, whether EMT during development also involves equivalent intermediate cell states is not known. To address this gap in knowledge, the authors devise a strategy to identify and characterize changes that an embryonic population of cells called the cranial neural crest undergo as they delaminate from the neuroepithelium and become a highly migratory population of mesenchymal cells that ultimately give rise to a broad range of derivatives.

To isolate and study the neural crest, the authors use embryos collected at E8.5 from two transgenic mouse lines. Wnt1-Cre;RosaeYFP labels Wnt1-positive neuroepithelial cells in the dorsolateral neural plate, which includes pre-migratory neural crest that reside in the dorsal neuroectoderm and neural plate border before induction (as well as some other lineages). Mef2c-F10N-LacZ leverages a neural crest cell-specific enhancer of Mef2c to control LacZ expression in predominantly migratory neural crest. This dual genetic approach that allows the authors to distinguish and compare pre-migratory and migratory neural crest cells is a strength of the work.

To assay for the differential expression of genes involved in the EMT and migration of cranial neural crest, the authors perform single cell RNA sequencing (scRNA-seq) using current methods. A strength is a large sample size per mouse line, and relatively high numbers of single cells analyzed. The authors identify six major cell/tissue types present in mouse E8.5 cranial tissues using known markers, which they then segregate into a cranial neural crest cluster using a well-reasoned bioinformatic strategy. The cranial neural crest cluster contains pre-migratory and migratory cells that they partition further into five subclusters and then characterize using the differential expression and combinatorial patterns of neural crest specifier genes, markers of pre-migratory neural crest, markers of early versus late migratory neural crest, markers of undifferentiated versus differentiated neural crest, tissue-specific markers, and region-specific markers. One weakness is that there is little attempt to map potential novel genes and/or pathways that also distinguish these clusters.

The authors then go on to subdivide the five cranial neural crest subclusters into almost two dozen smaller subclusters, again using the combinatorial expression of known markers (e.g., neural crest genes, cell junction genes, and cell cycle genes). A weakness is that the marker analysis and accompanying interpretation of the results relies heavily on the purported roles of different genes as described in the published work of others, which potentially introduces some untested assumptions and a bit of hand-waving into the study. Moreover, the limited correlation between mRNA and protein abundance for cell cycle markers is well documented in the literature but the authors rely heavily on gene expression to determine cell cycle status. Even though the authors add a compelling Edu/pHH3 double-labeling experiment and cell cycle inhibition studies, the work would be strengthened by including some analysis of protein expression to see if the cell cycle correlations hold up. Nonetheless, the subcluster and cell cycle analyses lead the authors to conclude that there are a series of intermediate cell states between neural crest EMT and delamination, and that cell cycle regulation is a defining feature and necessary component of those states. These novel findings are generally well supported by the data.

To test if there are spatiotemporal differences in the localization of neural crest cells during EMT in vivo, the authors apply a cutting-edge technique called signal amplification by exchange reaction for multiplexed fluorescent in situ hybridization (SABER-FISH), which they validate using standard in situ hybridization. The authors select specific marker genes that seem justified based on their scRNA-seq dataset, and they generate a series of convincing images and quantitative data that add valuable depth to the story.

As a functional test of their hypothesis that one of the genes indicative of an EMT intermediate stage (i.e., Dlc1) is essential for neural crest migration, the authors use a lentivirus-mediated knockdown strategy. A strength is that the authors include appropriate scramble and cell death controls as part of their experimental design.

The authors use Sox10 as a marker to count neural crest cells, but Sox10 may only label a subset of neural crest cells and thus some unaffected lineages may not have been counted. Although the data are persuasive, a second marker for counting neural crest cells following knockdown would make the analysis more robust.

Overall, this is a first-rate study with many more strengths than weaknesses. The authors generate high quality data, and their interpretations are reasonable and balanced. Another strength is the writing, which is clear and well organized, and the figures (including supplemental), which are excellent and provide unambiguous visualization of some very complex data sets. The methods are state-of the art and are effectively executed, and they will be useful to the broader cell and developmental biology community. The work contains well-substantiated findings and supports the conclusion that EMT is a highly dynamic, multi-step process, which was previously thought to be more-or-less binary. Such findings will alter the way the field thinks about EMT in neural crest and the work will likely serve as an important example alongside cancer metastasis.

Reviewer #3 (Public Review):

Summary:

Zhao et al. address the question of whether intermediate states of the epithelial-to-mesenchymal transition (EMT) exist in a natural developmental context as well as in cancer cells. This is important not only for our understanding of these developmental systems but also for their development as resources for new anti-cancer approaches. Guided by single-cell RNA sequencing analysis of delaminating mouse cranial neural crest cells, they identify two distinct populations with transcriptional signatures intermediate between neuroepithelial progenitors and migrating crest. Both clusters are also spatially intermediate and are actively cycling, with one in S-phase and one in G2/M. They show that blocking progression through S phase prior to the onset of delamination and knockdown of intermediate state marker Dlc1 both reduce the number of migratory cells that have completed EMT. Overall, the work provides a modern take and new insights into the classical developmental process of neural crest delamination.

Strengths:

• Deep analysis of the scRNAseq dataset revealed previously unappreciated cell populations intermediate between premigratory and migratory crest.<br /> • The observation that delaminating/intermediate neural crest cells appear to be in S or G2/M phase is interesting and worth reporting, though the ultimate significance remains unclear, given that they do not make distinct derivatives depending on their cycle state.<br /> • The authors employ new methods for multiplex spatial imaging to more accurately define their populations of interest and their relative positions.<br /> • The authors present evidence that intermediate state gene Dlc1 (a Rho GAP) is not just a marker but functionally required for neural crest delamination in mouse, as previously shown in chicken.

Weaknesses:

• Similar experiments involving blockade of cell cycle progression and Dlc1 dose manipulation were previously performed in chick models, as noted in the discussion. The newly-defined intermediate states give added context to the results, but they are not entirely novel.

Reviewer #1 (Public Review):

In this manuscript, the authors report a molecular mechanism for recruiting syntaxin 17 (Syn17) to the closed autophagosomes through the charge interaction between enriched PI4P and the C-terminal region of Syn17. How to precisely control the location and conformation of proteins is critical for maintaining autophagic flux. Particularly, the recruitment of Syn17 to autophagosomes remains unclear. In this paper, the author describes a simple lipid-protein interaction model beyond previous studies focusing on protein-protein interactions. This represents conceptual advances.

Reviewer #2 (Public Review):

Summary:

Syntaxin17 (STX17) is a SNARE protein that is recruited to mature (i.e., closed) autophagosomes, but not to immature (i.e., unclosed) ones, and mediates the autophagosome-lysosome fusion. How STX17 recognizes the mature autophagosome is an unresolved interesting question in the autophagy field. Shinoda and colleagues set out to answer this question by focusing on the C-terminal domain of STX17 and found that PI4P is a strong candidate that causes the STX17 recruitment to the autophagosome.

Strengths:

The main findings are: 1) Rich positive charges in the C-terminal domain of STX17 are sufficient for the recruitment to the mature autophagosome; 2) Fluorescence charge sensors of different strengths suggest that autophagic membranes have negative charges and the charge increases as they mature; 3) Among a battery of fluorescence biosensors, only PI4P-binding biosensors distribute to the mature autophagosome; 4) STX17 bound to isolated autophagosomes is released by treatment with Sac1 phosphatase; 5) By dynamic molecular simulation, STX17 TM is shown to be inserted to a membrane containing PI4P but not to a membrane without it. These results indicate that PI4P is a strong candidate that STX17 binds to in the autophagosome.

Weaknesses:

• It was not answered whether PI4P is crucial for the STX17 recruitment in cells because manipulation of the PI4P content in autophagic membranes was not successful for unknown reasons.<br /> • The question that the authors posed in the beginning, i.e., why is STX17 recruited to the mature (closed) autophagosome but not to immature autophagic membranes, was not answered. The authors speculate that the seemingly gradual increase of negative charges in autophagic membranes is caused by an increase in PI4P. However, this was not supported by the PI4P fluorescence biosensor experiment that showed their distribution to the mature autophagosome only. Here, there are at least two possibilities: 1) The increase of negative charges in immature autophagic membranes is derived from PI4P. However the fluorescence biosensors do not bind there for some reason; for example, they are not sensitive enough to recognize PI4P until it reaches a certain level, or simply, their binding does not occur in a quantitative manner. 2) The negative charge in immature membranes is not derived from PI4P, and PI4P is generated abundantly only after autophagosomes are closed. In either case, it is not easy to explain why STX17 is recruited to the mature autophagosome only. For the first scenario, it is not clear how the PI4P synthesis is regulated so that it reaches a sufficient level only after the membrane closure. In the second case, the mechanism that produces PI4P only after the autophagosome closure needs to be elucidated (so, in this case, the question of the temporal regulation issue remains the same).

Reviewer #3 (Public Review):

Summary:

In this study, the authors set out to address the question of how the SNARE protein Syntaxin 17 senses autophagosome maturation by being recruited to autophagosomal membranes only once autophagosome formation and sealing is complete. The authors discover that the C-terminal region of Syntaxin 17 is essential for its sensing mechanism that involves two transmembrane domains and a positively charged region. The authors discover that the lipid PI4P is highly enriched in mature autophagosomes and that electrostatic interaction with Syntaxin 17's positively charged region with PI4P drives recruitment specifically to mature autophagosomes. The temporal basis for PI4P enrichment and Syntaxin 17 recruitment to ensure that unsealed autophagosomes do not fuse with lysosomes is a very interesting and important discovery. Overall, the data are clear and convincing, with the study providing important mechanistic insights that will be of broad interest to the autophagy field, and also to cell biologists interested in phosphoinositide lipid biology. The author's discovery also provides an opportunity for future research in which Syntaxin 17's c-terminal region could be used to target factors of interest to mature autophagosomes.

Strengths:

The study combines clear and convincing cell biology data with in vitro approaches to show how Syntaxin 17 is recruited to mature autophagosomes. The authors take a methodical approach to narrow down the critical regions within Syntaxin 17 required for recruitment and use a variety of biosensors to show that PI4P is enriched on mature autophagosomes.

Weaknesses:

There are no major weaknesses, overall the work is highly convincing. It would have been beneficial if the authors could have shown whether altering PI4P levels would affect Syntaxin 17 recruitment. However, this is understandably a challenging experiment to undertake and the authors outlined their various attempts to tackle this question.

Reviewer #1 (Public Review):

Li et al report that upon traumatic brain injury (TBI), Pvr signalling in astrocytes activates the JNK pathway and up-regulates the expression of the well-known JNK target MMP1. The FACS sort astrocytes, and carry out RNAseq analysis, which identifies pvr as well as genes of the JNK pathway as particularly up-regulated after TBI. They use conventional genetics loss of function, gain of function and epistasis analysis with and without TBI to verify the involvement of the JNK-MMP1 signalling pathway downstream of PVR. They also show that blocking endocytosis prolongs the involvement of this pathway in the TBI response.

The strengths are that multiple experiments are used to demonstrate that TBI in their hands damaged the BBB, induced apoptosis and increased MMP1 levels. The RNAseq analysis on FACS sorted astrocytes is nice and will be valuable to scientists beyond the confines of this paper. The functional genetic analysis is conventional, yet sound, and supports claims of JNK and MMP1 functioning downstream of Pvr in the TBI context.

For this revised version the authors have removed all the unsupported claims. This renders their remaining claims more solid. However, it has resulted in the loss of important cellular aspects of the response to TBI, limiting the scope and value of the work.

The main weakness is that novelty and insight are both rather limited. Others had previously published that both JNK signalling and MMP1 were activated upon injury, in multiple contexts (as well as the articles cited by the authors, they should also see Losada-Perez et al 2021). That Pvr can regulate JNK signalling was also known (Ishimaru et al 2004). The authors claim that the novelty was investigating injury responses in astrocytes in Drosophila. However, others had investigated injury responses by astrocytes in Drosophila before. It had been previously shown that astrocytes - defined as the Prospero+ neuropile glia, and also sharing evolutionary features with mammalian NG2 glia - respond to injury both in larval ventral nerve cords and in adult brains, where they proliferate regenerating glia and induce a neurogenic response (Kato et al 2011; Losada-Perez et al 2016; Harrison et al 2021; Simoes et al 2022). The authors argue that the novelty of the work is the investigation of the response of astrocytes to TBI. However, this is of somewhat limited scope. The authors mention that MMP1 regulates tissue remodelling, the inflammatory process and cancer. Exploring these functions further would have been an interesting addition, but the authors did not investigate what consequences the up-regulation of MMP1 after injury has in repair or regeneration processes.

The statistical analysis is incorrect in places, and this could affect the validity of some claims.

Altogether, this is an interesting and valuable addition to the repertoire of articles investigating neuron-glia communication and glial responses to injury in the Drosophila central nervous system (CNS). It is good and important to see this research area in Drosophila grow. This community together is building a compelling case for using Drosophila and its unparalleled powerful genetics to investigate nervous system injury, regeneration and repair, with important implications. Thus, this paper will be of interest to scientists investigating injury responses in the CNS using Drosophila, other model organisms (eg mice, fish) and humans.

Reviewer #3 (Public Review):

In this study, authors used the Drosophila model to characterize molecular details underlying traumatic brain injury (TBI). Authors used the transcriptomic analysis of astrocytes collected by FACS sorting of cells derived from Drosophila heads following brain injury and identified upregulation of multiple genes, such as Pvr receptor, Jun, Fos, and MMP1. Additional studies identified that Pvr positively activates AP-1 transciption factor (TF) complex consisting of Jun and Fos, of which activation leads to the induction of MMP1. Finally, authors found that disruption of endocytosis and endocytotic trafficking facilitates Pvr signaling and subsequently leads to induction of AP-1 and MMP1.

Overall, this study provides important clues to understanding molecular mechanisms underlying TBI. The identified molecules linked to TBI in astrocytes could be potential targets for developing effective therapeutics. The obtained data from transcriptional profiling of astrocytes will be useful for future follow-up studies. The manuscript is well-organized and easy to read.

However, the connection suggested by the authors between Pvr and AP-1, potentially mediated through the JNK pathway, lacks strong experimental support in my view. It's important to recognize that AP-1 activity is influenced by multiple upstream signaling pathways, not just the JNK pathway, which is the most well-characterized among them. Therefore, assuming that AP-1 transcriptional activity solely reflects the activity of the JNK pathway without additional direct evidence is unwarranted. To strengthen their argument, the study could benefit from direct evidence implicating the JNK pathway in linking Pvr to AP-1. This could be achieved through genetic studies involving mutants or transgenes targeting key components of the JNK pathway, such as Bsk and Hep, the Drosophila homologues of JNK and JNKK, respectively. Alternatively, employing p-JNK antibody-based techniques like Western blotting, while considering the potential challenges associated with p-JNK immunohistochemistry, could provide further validation. This important criticism regarding the molecular link between Pvr and AP-1 has been overlooked.

Reviewer #1 (Public Review):

Dasguta et al. have dissected the role of Sema7a in fine tuning of a sensory microcircuit in the posterior lateral line organ of zebrafish. They attempt to also outline the different roles of a secreted verses membrane-bound form of Sema7a in this process. Using genetic perturbations and axonal network analysis, the authors show that loss of both Sema7a isoforms causes abnormal axon terminal structure with more bare terminals and fewer loops in contact with presynaptic sensory hair cells. Further, they show that loss of Sema7a causes decreased number and size of both the pre- and post-synapse. Finally, they show that overexpression of the secreted form of Sema7a specifically can elicit axon terminal outgrowth to an ectopic Sema7a expressing cell. Together, the analysis of Sema7a loss of function and overexpression on axon arbor structure is fairly thorough and revealed a novel role for Sema7a in axon terminal structure. However, the connection between different isoforms of Sema7a and the axon arborization needs to be substantiated. Furthermore, the effect of loss of Sema7a on the presynaptic cell is not ruled out as a contributing factor to the synaptic and axon structure phenotypes. These issues weaken the claims made by the authors including the statement that they have identified dual roles for the GPI-anchored verses secreted forms of Sema7a on synapse formation and as a chemoattractant for axon arborization respectively.

Reviewer #2 (Public Review):

In this work, Dasgupta et al. investigates the role of Sema7a in the formation of peripheral sensory circuit in the lateral line system of zebrafish. They show that Sema7a protein is present during neuromast maturation and localized, in part, to the base of hair cells (HCs). This would be consistent with pre-synaptic Sema7a mediating formation and/or stabilization of the synapse. They use sema7a loss-of-function strain to show that lateral line sensory terminals display abnormal arborization. They provide highly quantitative analysis of the lateral line terminal arborization to show that a number of specific topological parameters are affected in mutants. Next, they ectopically express a secreted form of Sema7a to show that lateral line terminals can be ectopically attracted to the source. Finally, they also demonstrate that the synaptic assembly is impaired in the sema7a mutant. Overall, the data are of high quality and properly controlled. The availability of Sema7a antibody is a big plus, as it allows to address the endogenous protein localization as well to show the signal absence in the sema7a mutant. The quantification of the arbor topology should be useful to people in the field who are looking at the lateral line as well as other axonal terminals. I think some results are overinterpreted though. The authors state: "Our findings demonstrate that Sema7A functions both as a juxtracrine and as a secreted cue to pattern neural circuitry during sensory organ development." However, they have not actually demonstrated which isoform functions in HCs (also see comments below). In addition, they have to be careful in interpreting their topology analysis, as they cannot separate individual axons. Thus, such analysis can generate artifacts. They can perform additional experiments to address these issues or adjust their interpretations.

Reviewer #3 (Public Review):

The data reported here demonstrate that Sema7a defines the local behavior of growing axons in the developing zebrafish lateral line. The analysis is sophisticated and convincingly demonstrates effects on axon growth and synapse architecture. Collectively, the findings point to the idea that the diffusible form of sema7a may influence how axons grow within the neuromast and that the GPI-linked form of sema7a may subsequently impact how synapses form, though additional work is needed to strongly link each form to its' proposed effect on circuit assembly.

Comments on revised submission:

The revised manuscript is significantly improved. The authors comprehensively and appropriately addressed most of the reviewers' concerns. In particular, they added evidence that hair cells express both Sema7A isoforms, showed that membrane bound Sema7A does not have long range effects on guidance, demonstrated how axons behave close to ectopic Sema7A, and analyzed other features of the hair cells that revealed no strong phenotypes. The authors also softened the language in many, but not all places. Overall, I am satisfied with the study as a whole.

Reviewer #4 (Public Review):

This study provides direct evidence showing that Sema7a plays a role in the axon growth during the formation of peripheral sensory circuits in the lateral-line system of zebrafish. This is a valuable finding because the molecules for axon growth in hair-cell sensory systems are not well understood. The majority of the experimental evidence is convincing, and the analysis is rigorous. The evidence supporting Sema7a's juxtracrine vs. secreted role and involvement in synapse formation in hair cells is less conclusive. The study will be of interest to cell, molecular and developmental biologists, and sensory neuroscientists.

Reviewer #1 (Public Review):

Summary:

Through an unbiased genomewide KO screen, the authors identified loss of DBT to suppress MG132-mediated death of cultured RPE cells. Further analyses suggested that DBT reduces ubiquitinated proteins by promoting autophagy. Mechanistic studies indicated that DBT loss promotes autophagy via AMPK and its downstream ULK and mTOR signaling. Furthermore, loss of DBT suppresses polyglutamine- or TDP-43-mediated cytotoxicity and/or neurodegeneration in fly models. Finally, the authors showed that DBT proteins are increased in ALS patient tissues, compared to non-neurological controls.

Strengths:

The idea is novel, the evidence is convincing, and the data are clean. The findings have implications for human diseases.

Weaknesses:

None.

Reviewer #2 (Public Review):

Summary:

Hwang, Ran-Der et al utilized a CRISPR-Cas9 knockout in human retinal pigment epithelium (RPE1) cells to evaluate for suppressors of toxicity by the proteasome inhibitor MG132 and identified that knockout of dihydrolipoamide branched chain transacylase E2 (DBT) suppressed cell death. They show that DBT knockout in RPE1 cells does not alter proteasome or autophagy function at baseline. However, with MG132 treatment, they show a reduction in ubiquitinated proteins but with no change in proteasome function. Instead, they show that DBT knockout cells treated with MG132 have improved autophagy flux compared to wildtype cells treated with MG132. They show that MG132 treatment decreases ATP/ADP ratios to a greater extent in DBT knockout cells, and in accordance causes activation of AMPK. They then show downstream altered autophagy signaling in DBT knockout cells treated with MG132 compared to wild-type cells treated with MG132. Then they express the ALS mutant TDP43 M337 or expanded polyglutamine repeats to model Huntington's disease and show that knockdown of DBT improves cell survival in RPE1 cells with improved autophagic flux. They also utilize a Drosophila models and show that utilizing either a RNAi or CRISPR-Cas9 knockout of DBT improves eye pigment in TDP43M337V and polyglutamine repeat-expressing transgenic flies. Finally, they show evidence for increased DBT in postmortem spinal cord tissue from patients with ALS via both immunoblotting and immunofluorescence.

Strengths:

This is a mechanistic and well-designed paper that identifies DBT as a novel regulator of proteotoxicity via activating autophagy in the setting of proteasome inhibition. Major strengths include careful delineation of a mechanistic pathway to define how DBT is protective. These conclusions are well-justified.

Weaknesses:

None

Reviewer #1 (Public Review):

Summary:

In this study, the authors distinguished afferent inputs to different cell populations in the VTA using dimensionality reduction approaches and found significantly distinct patterns between normal and drug treatment conditions. They also demonstrated negative correlations of the inputs induced by drugs with gene expression of ion channels or proteins involved in synaptic transmission and demonstrated the knockdown of one of the voltage-gated calcium ion channels caused decreased inputs.

Weaknesses:

(1) For quantifications of brain regions in this study, boundaries were based on the Franklin-Paxinos (FP) atlas according to previous studies (Beier KT et al 2015, Beier KT et al 2019). It has been reported significant discrepancies exist between the anatomical labels on the FP atlas and the Allen Brain Atlas (ref: Chon U et al., Nat Commun 2019). Although a summary of conversion is provided as a sheet, the authors need to describe how consistent or different the brain boundaries they defined in the manuscript with Allen Brain Atlas by adding histology images. Also, I wonder how reliable the annotations were for over a hundred of animals with manual quantification. The authors should briefly explain it rather than citing previous studies in the Material and Methods Section.

(2) Regarding the ellipsoids in the PC, although it's written in the manuscript that "Ellipsoids were centered at the average coordinate of a condition and stretched one standard deviation along the primary and secondary axes", it's intuitively hard to understand in some figures such as Figure 2O, P and Figure S1. The authors need to make their data analysis methods more accessible by providing source code to the public.

(3) In histology images (Figure 1B and 3K), the authors need to add dashed lines or arrows to guide the reader's attention.

(4) In Figure 2A and G, apparently there are significant differences in other brain regions such as NAcMed or PBN. If they are also statistically significant, the authors should note them as well and draw asterisks(*).

(5) In Figure 2N about the spatial distribution of starter cells, the authors need to add histology images for each experimental condition (i.e. saline, fluoxetine, cocaine, methamphetamine, amphetamine, nicotine, and morphine) as supplement figures.

(6) In the manuscript, it is necessary to explain why Cacna1e was selected among other calcium ion channels.

Reviewer #2 (Public Review):

The application of rabies virus (RabV)-mediated transsynaptic tracing has been widely utilized for mapping cell-type-specific neural connectivities and examining potential modifications in response to biological phenomena or pharmacological interventions. Despite the predominant focus of studies on quantifying and analyzing labeling patterns within individual brain regions based on labeling abundance, such an approach may inadvertently overlook systemic alterations. There exists a considerable opportunity to integrate RabV tracing data with the global connectivity patterns and the transcriptomic signatures of labeled brain regions. In the present study, the authors take an important step towards achieving these objectives.

Specifically, the authors conducted an intensive reanalysis of a previously generated large dataset of RabV tracing to the ventral tegmental area (VTA) using dimension reduction methods such as PCA and UMPA. This reaffirmed the authors's earlier conclusion that different cell types in the VTA, namely dopamine neurons (DA) and GABAergic neurons, exhibit quantitatively distinct input patterns, and a single dose of addictive drugs, such as cocaine and morphine, induced altered labeling patterns. Additionally, the authors illustrate that distinct axes of PCA can discriminate experimental variations, such as minor differences in the injection site of viral tracers, from bona fide alternations in labeling patterns caused by drugs of abuse. While the specific mechanisms underlying altered labeling in most brain regions remain unclear, whether involving synaptic strength, synaptic numbers, pre-synaptic activities, or other factors, the present study underscores the efficacy of an informatics approach in extracting more comprehensive information from the RabV-based circuit mapping data.

Moreover, the authors showcased the utility of their previously devised bulk gene expression patterns inferred by the Allen Gene Expression Atlas (AGEA) and "projection portrait" derived from bulk axon mapping data sourced from the Allen Mouse Brain Connectivity Atlas. The utilization of such bulk data rests upon several limitations. For instance, the collection of axon mapping data involves an arbitrary selection of both cell type-specific and non-specific data, which might overlook crucial presynaptic partners, and often includes contamination from neighboring undesired brain regions. Concerns arise regarding the quantitativeness of AGEA, which may also include the potential oversight of key presynaptic partners. Nevertheless, the authors conscientiously acknowledged these potential limitations associated with the dataset.

Notably, building on the observation of a positive correlation between the basal expression levels of Ca2+ channels and the extent of drug-induced changes in RabV labeling patterns, the authors conducted a CRISPRi-based knockdown of a single Ca2+ channel gene. This intervention resulted in a reduction of RabV labeling, supporting that the observed gene expression patterns have causality in RabV labeling efficiency. While a more nuanced discussion is necessary for interpreting this result (see below), overall I commend the authors for their efforts to leverage the existing dataset in a more meaningful way. This endeavor has the potential to contribute significantly to our understanding of the mechanisms underlying alterations in RabV labeling induced by drugs of abuse.

Finally, drawing upon the aforementioned reanalysis of previous data, the authors underscored that a single administration of ketamine/xylazine anesthesia could induce enduring modifications in RabV labeling patterns for VTA DA neurons, specifically those projecting to the nucleus accumbens and amygdala. Given the potential impact of such alterations on motivational behaviors at a broader level, I fully agree that prudent consideration is warranted when employing ketamine/xylazine for the investigation of motivational behaviors in mice.

Specific Points:

(1) Beyond advancements in bioinformatics, readers may find it insightful to explore whether the PCA/UMPA-based approach yields novel biological insights. For example, the authors are encouraged to discuss more functional implications of PBN and LH in the context of drugs of abuse, as their labeling abundance could elucidate the PC2 axis in Fig. 2M.

2) While I appreciate the experimental data on Cacna1e knockdown, I am unclear about the rationale behind specifically focusing on Cacna1e. The logic behind the statement, "This means that expression of this gene is not inhibitory towards RABV transmission," is also unclear. Loss-of-function experiments only signify the necessity or permissive functions of a gene. In this context, Cacna1e expression levels are required for efficient RabV labeling, but this neither supports nor excludes the possibility that this gene expression instructively suppresses RabV labeling/transmission, which could be assessed through gain-of-function experiments.

Reviewer #3 (Public Review):

Summary:

Authors mapped monosynaptic inputs to dopamine, GABA, and glutamate neurons in VTA under different anesthesia methods, and under drugs (cocaine, morphine, methamphetamine, amphetamine, nicotine, fluoxetine). They found that input patterns under different conditions are separated, and identified some key brain areas to contribute to such separation. They also searched a database for gene expression patterns that are common across input brain areas with some changes by anesthesia or drug administration.

Strengths:

The whole-brain approach to address drug effects is appealing and their conclusion is clear. The methodology and motivation are clearly explained.

Weaknesses:

While gene expression analyses may not be related to their findings on the anatomical effects of drugs, this will be a nice starting point for follow-up studies.

Reviewer #2 (Public Review):

The authors analysed functional MRI recordings of brain activity at rest, using state-of-the-art methods that reveal the diverse ways in which the information can be integrated in the brain. In this way, they found brain areas that act as (synergistic) gateways for the 'global workspace', where conscious access to information or cognition would occur, and brain areas that serve as (redundant) broadcasters from the global workspace to the rest of the brain. The results are compelling and consisting with the already assumed role of several networks and areas within the Global Neuronal Workspace framework. Thus, in a way, this work comes to stress the role of synergy and redundancy as complementary information processing modes, which fulfill different roles in the big context of information integration.<br /> In addition, to prove that the identified high-order interactions are relevant to the phenomenon of consciousness, the same analysis was performed in subjects under anesthesia or with disorders of consciousness (DOC), showing that indeed the loss of consciousness is associated with a deficient integration of information within the gateway regions.

However, there is still a standing issue that could be the basis for an improved analysis: the concepts of gateways and broadcasters allude to a directionality in the information flow. In fact, Figure 1 depicts Stage (i) and Stage (iii) as one-way processes. However, the identification of gateway and broadcaster regions relies on matrices that are symmetrical, i.e. they are not directed. Would it be possible to assess the gateway or broadcaster nature of a region taking into account the directionality of the information flow? In other words, if region X is a gateway, I would expect a synergistic relationship between the past of X,Y and present of Y (Y not being a gateway) towards the present of X; but not necessarily the other way around (i.e. the present of Y being less dependent on the past/present of X). A similar reasoning can be made for broadcasters.

Although regional differences in haemodynamics complicate attempts to map directed information flow from fMRI recordings, perhaps the IID framework could be leveraged to extract directed data (i.e., there are many atoms that are explicitly directed). As an avenue for future research, it would be interesting to discuss the feasibility or limitations of such analysis.

Also, there is something confusing in Figure 4B-C and its description. Awake should be similar to recovery (they are both awake, aren't they? Not much info is given, anyway); thus it seems counterintuitive that anesthesia minus awake looks so different than anesthesia minus recovery. The first is mostly blue-ish and the second is mostly red. Is it possible that Figure 4C is actually recovery minus anesthesia? That would make much more sense, also for Figure 4D. Please correct me if I am wrong.

Reviewer #3 (Public Review):

The work proposes a model of neural information processing based on a 'synergistic global workspace,' which processes information in three principal steps: a gatekeeping step (information gathering), an information integration step, and finally, a broadcasting step. They provided an interpretation of the reduced human consciousness states in terms of the proposed model of brain information processing, which could be helpful to be implemented in other states of consciousness. The manuscript is well-organized, and the results are important and could be interesting for a broad range of literature, suggesting interesting new ideas for the field to explore.

Comments on revised version:

The authors have addressed all my comments made in the previous revision.

Reviewer #3 (Public Review):

Strengths:<br /> There are many reports on the effect of chemical properties of foods on feeding in fruit flies, but only limited studies reported how physical properties of food affect feeding especially pharyngeal mechanosensory neurons. First, they found that mechanosensory mutants, including nompC, Tmc, and Piezo, showed impaired swallowing, mainly the emptying process. Next, they identified cibarium multidendritic mechanosensory neurons (md-C) are responsible for controlling swallowing by regulating motor neuron (MN) 12 and 11, which control filling and emptying, respectively.

Weaknesses:<br /> While the involvement of md-C and mechanosensory channels in controlling swallowing is convincing, it is not yet clear which stimuli activate md-C. Can it be an expansion of cibarium or food viscosity, or both? In addition, if rhythmic and coordinated contraction of muscles 11 and 12 is essential for swallowing, how can simultaneous activation of MN 11 and 12 by md-C achieve this? Finally, previous reports showed that food viscosity mainly affects the filling rather than the emptying process, which seems different from their finding.

Reviewer #1 (Public Review):

Qin et al. set out to investigate the role of mechanosensory feedback during swallowing and identify neural circuits that generate ingestion rhythms. They use Drosophila melanogaster swallowing as a model system, focusing their study on the neural mechanisms that control cibarium filling and emptying in vivo. They find that pump frequency is decreased in mutants of three mechanotransduction genes (nompC, piezo, and Tmc), and conclude that mechanosensation mainly contributes to the emptying phase of swallowing. Furthermore, they find that double mutants of nompC and Tmc have more pronounced cibarium pumping defects than either single mutants or Tmc/piezo double mutants. They discovered that the expression patterns of nompC and Tmc overlap in two classes of neurons, md-C and md-L neurons. The dendrites of md-C neurons warp the cibarium and project their axons to the subesophageal zone of the brain. Silencing neurons that express both nompC and Tmc leads to severe ingestion defects, with decreased cibarium emptying. Optogenetic activation of the same population of neurons inhibited filling of the cibarium and accelerated cibarium emptying. In the brain, the axons of nompC∩Tmc cell types respond during ingestion of sugar but do not respond when the entire fly head is passively exposed to sucrose. Finally, the authors show that nompC∩Tmc cell types arborize close to the dendrites of motor neurons that are required for swallowing and that swallowing motor neurons respond to the activation of the entire Tmc-GAL4 pattern.

Strengths:<br /> -The authors rigorously quantify ingestion behavior to convincingly demonstrate the importance of mechanosensory genes in the control of swallowing rhythms and cibarium filling and emptying<br /> -The authors demonstrate that a small population of neurons that express both nompC and Tmc oppositely regulate cibarium emptying and filling when inhibited or activated, respectively<br /> -They provide evidence that the action of multiple mechanotransduction genes may converge in common cell types

Weaknesses:<br /> -A major weakness of the paper is that the authors use reagents that are expressed in both md-C and md-L but describe the results as though only md-C is manipulated<br /> -Evidence that the defects they see in pumping can be specifically attributed to md-C is based on severing the labellum and allowing md-L neurons to degrade.<br /> -GRASP is known to be non-specific and prone to false positives when neurons are in close proximity but not synaptically connected. A positive GRASP signal supports but does not confirm direct synaptic connectivity between md-C/md-L axons and MN11/MN12.<br /> -MN11/MN12 LexA lines are not included in the manuscript and their expression patterns (shared with the reviewers in the author response) do not appear to contain any motor neurons. Double labeling with previously described MN11 and MN12 motor neuron Gal4 lines is needed to support the claim that these LexA lines in fact label MN11 and MN12.<br /> -As seen in Figure Supplement 2, the expression pattern of Tmc-GAL4 is broader than md-C alone. Therefore, the functional connectivity the authors observe between Tmc expressing neurons and MN11 and 12 cannot be traced to md-C alone<br /> -Example traces of md-C calcium imaging during ingestion in vivo are not included, and evidence that md-C neurons respond to mechanical force is lacking<br /> -A positive control (perhaps demonstrating that sugar sensory neurons respond to sucrose in this preparation) is needed to assess whether the lack of response to sucrose ex vivo in Figure 4K is informative<br /> -Proximity between md-C neurons and muscles is not evidence that they sense stretch<br /> -Reporting of posthoc tests needs to be improved throughout the manuscript, as it is not clear which comparisons are noted with asterisks in the figures.

Overall, this work convincingly shows that swallowing and swallowing rhythms are dependent on several mechanosensory genes. Qin et al. also characterize a candidate neuron, md-C, that is likely to provide mechanosensory feedback to pumping motor neurons, but the results they present here are not sufficient to assign this function to md-C alone. This work will have a positive impact on the field by demonstrating the importance of mechanosensory feedback to swallowing rhythms and providing a potential entry point for future investigation of the identity and mechanisms of swallowing central pattern generators.

Reviewer #3 (Public Review):

Summary:

The authors aimed to develop an automated tool to easily collect, process, and annotate the biomedical literature for higher efficiency and better reproducibility.

Strengths:

Two charms coming with the efforts made by the team are Pubget (for efficient and reliable grabbing articles from PubMed) and labelbuddy (for annotating text). They make text-mining of the biomedical literature more accessible, effective, and reproducible for streamlined text-mining and meta-science projects. The data were collected and analyzed using solid and validated methodology and demonstrated a very promising direction for meta-science studies.

Weaknesses:

More developments are needed for different resources of literature and strengths of AI-powered functions.

Reviewer #1 (Public Review):

Summary:

In this paper, the authors present new tools to collect and process information from the biomedical literature that could be typically used in a meta-analytic framework. The tools have been specifically developed for the neuroimaging literature. However, many of their functions could be used in other fields. The tools mainly enable to downloading of batches of paper from the literature, extracting relevant information along with meta-data, and annotating the data. The tools are implemented in an open ecosystem that can be used from the command line or Python.

Strengths:

The tools developed here are really valuable for the future of large-scale analyses of the biomedical literature. This is a very well-written paper. The presentation of the use of the tools through several examples corresponding to different scientific questions really helps the readers to foresee the potential application of these tools.

Weaknesses:

The tools are command-based and store outcomes locally. So users who prefer to work only with GUI and web-based apps may have some difficulties. Furthermore, the outcomes of the tools are constrained by inherent limitations in the scientific literature, in particular, here the fact that only a small portion of the publications have full text openly available.

Reviewer #2 (Public Review):

Summary:

In this manuscript, the authors described the litmining ecosystem that can flexibly combine automatic and manual annotation for meta-research.

Strengths:

Software development is crucial for cumulative science and of great value to the community. However, such works are often greatly under-valued in the current publish-or-perish research culture. Thus, I applaud the authors' efforts devoted to this project. All the tools and repositories are public and can be accessed or installed without difficulty. The results reported in the manuscript are also compelling that the ecosystem is relatively mature.

Weaknesses:

First and foremost, the logic flow of the current manuscript is difficult to follow.

The second issue is the results from the litmining ecosystem were not validated and the efficiency of using litmining was not quantified. To validate the results, it would be better to directly compare the results of litmining with recognized ground truth in each of the examples. To prove the efficiency of the current ecosystem, it would be better to use quantitative indices for comparing the litmining and the other two approaches (in terms of time and/or other costs in a typical meta-research).

The third family of issues is about the functionality of the litmining ecosystem. As the authors mentioned, the ecosystem can be used for multiple purposes, however, the description here is not sufficient for researchers to incorporate the litmining ecosystem into their meta-research project. Imagine that a group of researchers are interested in using the litmining ecosystem to facilitate their meta-analyses, how should they incorporate litmining into their workflow? I have this question because, in a complete meta-analysis, researchers are required to (1) search in more than one database to ensure the completeness of their literature search; (2) screen the articles from the searched articles, which requires inspection of the abstract and the pdf; (3) search all possible pdf file of included articles instead of only relying on the open-access pdf files on PMC database. That said, if researchers are interested in using litmining in a meta-analysis that follows reporting standards such as PRISMA, the following functionalities are crucial:<br /> (a) How to incorporate the literature search results from different databases;<br /> (b) After downloading the meta-data of articles from databases, how to identify whose pdf files can be downloaded from PMC and whose pdf files need to be searched from other resources;<br /> (c) Is it possible to also annotate pdf files that were not downloaded by pubget?<br /> (d) How to maintain and update the meta-data and intermediate data for a meta-analysis by using litmining? For example, after searching in a database using a specific command and conducting their meta-analysis, researchers may need to update the search results and include items after a certain period.

Reviewer #1 (Public Review):

Original Review:

Bischoff et al present a carefully prepared study on a very interesting and relevant topic: the role of ion channels (here a Ca2+-activated K+ channel BK) in regulating mitochondrial metabolism in breast cancer cells. The potential impact of these and similar observations made in other tumor entities has only begun to be appreciated. That being said, the authors pursue in my view an innovative approach to understanding breast cancer cell metabolism.

Considering the following points would further strengthen the manuscript:

Methods:

(1) The authors use an extracellular Ca2+ concentration (2 mM) in their Ringer's solutions that is almost twice as high as the physiologically free Ca2+ concentration (ln 473). Moreover, the free Ca2+ concentration of their pipette solution is not indicated (ln 487).

(2) Ca2+I measurements: The authors use ATP to elicit intracellular Ca2+ signals. Is this then physiological stimulus for Ca2+ signaling in breast cancer? What is the rationale for using ATP? Moreover, it would be nice to see calibrated baseline values of Ca2+i

(3) Membrane potential measurements: It would be nice to see a calibration of the potential measurements; this would allow to correlate IV relationship with membrane potential. Without calibration it is hard to compare unless the identical uptake of the dye is shown.<br /> Do paxilline or IbTx also induce a depolarization?

(4) mito-potential measurements: Why did the authors use such a long time course and preincubated cells mit channel blockers overnight? Why did they not perform paired experiments and record the immediate effect of the BK channel blockers in the mito potential?

(5) MTT assay are also based on mitochondrial function - since modulation of mito function is at the core of this manuscript, an alternative method should be used.

Results:

(1) Fig. 5G: The number of BK "positive" mitoplasts is surprisingly low - how does this affect the interpretation? Did the authors attempt to record mitoBK current in the "whole-mitoplast" mode? How does the mitoBK current density compare with that of the plasma membrane? Is it possible to theoretically predict the number of mitoBK channels per mitochondrium to elicit the observed effects? Can these results be correlated with immuno-localization of mitoBK channels?

(2) There are also reports about other mitoK channels (e.g. Kv1.3, KCa3.1, KATP) playing an important role in mitochondrial function. Did the authors observe them, too? Can the authors speculate on the relative importance of the different channels? Is it known whether they are expressed organ-/tumor-specifically?

Comments on revised version:

The authors responded to all of my comments - except for one - in a satisfactory way so that I have no further concerns. The authors have prepared a very interesting piece of work that advances the field.

However, I disagree with respect to their interpretation of statistics. Individually analyzed cells are not the best biological replicate per se. In my view a true replicate requires the use of an independent batch of cells derived from a new passage. The statistical analysis can only based on the total number of n cells, if each replicate contributes the same number of cells. If this is not the case, the authors will have to calculate the average of each replicate first so that they are equally weighted.

Reviewer #2 (Public Review):

Summary:

The large-conductance Ca2+ activated K+ channel (BK) has been reported to promote breast cancer progression, but it is not clear how. The present study, carried out in breast cancer cell lines, concludes that BK located in mitochondria reprograms cells towards the Warburg phenotype, one of the metabolic hallmarks of cancer.

Strengths:

The use of a wide array of modern complementary techniques, including metabolic imaging, respirometry, metabolomics and electrophysiology. On the whole experiments are astute and well designed, and appear carefully done. The use of a BK knock out cells to control for the specificity of the pharmacological tools is a major strength. The manuscript is clearly written. There are many interesting original observations that may give birth to new studies.

Weaknesses: The main conclusion regarding the role of a BK channel located in mitochondria appears is not sufficiently supported. Other perfectible aspects are the interpretation of co-localization experiments and the calibration of Ca2+ dyes. These points are discussed in more detail in the following paragraphs:

(1) May the metabolic effects be ascribed to a BK located in mitochondria? Unfortunately not, at least with the available evidence. While it is clear these cells have a BK in mitochondria (characteristic K+ currents detected in mitoplasts) and it is also well substantiated that the metabolic effects in intact cells are explained by an intracellular BK (paxilline effects absent in the BK KO), it does not follow that both observations are linked. Given that ectopic BK-DEC appeared at the surface, a confounding factor is the likely expression of BK in other intracellular locations such as ER, Golgi, endosomes, etc. To their credit authors acknowledge this limitation several times throughout the text ("...presumably mitoBK...") but not in other important places, particularly in title and abstract.

(2) mitoBK subcellular location. Pearson correlations of 0.6 and about zero were obtained between the locations of mitoGREEN on one side, and mRFP or RFP-GPI on the other (Figs. 1G and S1E). These are nice positive and negative controls. For BK-DECRFP however the Pearson correlation was about 0.2. What is the Z resolution of apotome imaging? Assuming an optimum optical section of 600 nm, as obtained a 1.4 NA objective with a confocal, that mitochondria are typically 100 nm in diameter and that BK-DECRFP appears to stain more structures that mitoGREEN, the positive correlation of 0.2 may not reflect colocalization. For instance, it could be that BK-DECRFP in not just in mitochondria but in a close underlying organelle e.g. the ER. Along the same line, why did BK-RFP also give a positive Pearson? Isn´t that unexpected? Considering that BK-DEC was found by patch clamping at the plasma membrane, the subcellular targeting of the channel is suspect. Could it be that the endogenous BK-DEC does actually reside exclusively in mitochondria (a true mitoBK), but overflows to other membranes upon overexpression? Regarding immunodetection of BK in the mitochondrial Percoll preparation (Fig. S5), absence of NKA demonstrates absence of plasma membrane contamination, but does not inform about contamination by other intracellular membranes.

(3) Calibration of fluorescent probes. The conclusion that BK blockers or BK expression affects resting Ca2+ levels should be better supported. Fluorescent sensors and dyes provide signals or ratios that need be calibrated if comparisons between different cell types or experimental conditions are to be made. This is implicitly acknowledged here when monitoring ER Ca2+, with an elaborate protocol to deplete the organelle in order to achieve a reading at zero Ca2+.

(4) Line 203. "...solely by the expression of BKCa-DECRFP in MCF-7 cells". Granted, the effect of BKCa-DECRFP on the basal FRET ratio appears stronger than that of BK-RFP, but it appears that the latter had some effect. Please provide the statistics of the latter against the control group (after calibration, see above).

The revised version of the manuscript has incorporated my suggestions to a very reasonable degree, in several cases with new experiments. The details of these improvements can be found in the correspondence.

Reviewer #3 (Public Review):

The original research article, titled "mitoBKCa is functionally expressed in murine and human breast cancer cells and promotes metabolic reprogramming" by Bischof et al, has demonstrated the underlying molecular mechanisms of alterations in the function of Ca2+ activated K+ channel of large conductance (BKCa) in the development and progression of breast cancer. The authors also proposed that targeting mitoBKCa in combination with established anti-cancer approaches, could be considered as a novel treatment strategy in breast cancer treatment.

The paper is modified according to the reviewer's comments. Most of the queries raised by this reviewer were answered. However, the preclinical implication of this study can also be manifested in combinatorial treatment with known chemotherapeutic drugs which is lacking in this manuscript. Hopefully, the authors will consider this in their future study.

Joint Public Review:

Neuropeptide signaling is an important component of nervous systems, where neuropeptides typically act via G-protein coupled receptors (GPCRs) to regulate many physiological and behavioral processes. Neuropeptides and their cognate GPCRs have been extensively characterized in bilaterian animals, revealing that a core set of neuropeptide signaling systems originated in common ancestors of extant Bilateria. Neuropeptides have also been identified in cnidarians, which are a sister group to the Bilateria. However, the GPCRs that mediate the effects of neuropeptides in cnidarians have not been identified.

In this paper the authors perform a phylogenetic analysis of GPCRs in metazoans and report that the orthologs of bilaterian neuropeptide receptors are not found in cnidarians. This indicates that neuropeptide signaling systems have largely evolved independently in cnidarians and bilaterians. To accomplish this, they generated a library of putative and known neuropeptides computationally identified in the genome of the cnidarian sea anemone Nematostella vectensis. These peptides were systematically screened for their ability to activate any of the 161 putative Nematostella GPCRs.

This work identified 31 validated GPCRs. These, together with GPCRs that cluster with them, were then used to demonstrate the independent expansion of GPCRs in cnidarian and bilaterian lineages. The authors then mapped validated receptors and ligands to the Nematostella single cell data to provide an overview of the cell types expressing these signaling genes. In addition, the authors have begun to analyze neuropeptide signaling networks in N. vectensis by showing potential signaling connections between cell types expressing neuropeptides and cell types expressing cognate receptors.

This work is the most extensive pharmacological characterization of neuropeptide GPCRs in a cnidarian to date and thus represents an important accomplishment, and is one that will improve our understanding of how peptidergic signaling evolved in animals and its impact on evolution of nervous systems. In addition, this impressive work transforms our knowledge of neuropeptide signaling systems in cnidarians and provides the foundations for extensive functional characterization neuropeptide systems in the context of nervous systems that exhibit radial symmetry, contrasting with the bilaterally symmetrical architecture of the majority of bilaterian nervous systems.

The reviewers did not detect any weaknesses in the work but asked that the authors comment on the following points, which they have done in the revised version.

(1) Clearly, other neuropeptide signaling systems in cnidarians remain to be discovered but this paper represents a huge step forward.

(2) There are limitations in what can be interpreted from single cell transcriptomic data but the data nevertheless provide the foundations for future studies involving i). detailed anatomical analysis of neuropeptide and neuropeptide receptor expression in N. vectensis using mRNA in situ hybridization and/or immunohistochemical methods and ii). functional analysis of the physiological/behavioral roles of neuropeptide signaling systems in N. vectensis.

Reviewer #1 (Public Review):

Summary:

In the resubmitted manuscript by Chen et al. entitled, "Retinal metabolism: Evidence for uncoupling of glycolysis and oxidative phosphorylation via Cori-, Cahill-, and mini-Krebs-cycle", the authors look to provide insight on retinal metabolism and substrate utilization but using a murine explant model with various pharmacological treatments in conjunction with metabolomics. The authors conclude that photoreceptors, a specific cell within the explant, which also includes retinal pigment epithelium (RPE) and many other types of cells, are able to uncouple glycolytic and Krebs-cycle metabolism via three different pathways: 1) the mini-Krebs-cycle, fueled by glutamine and branched-chain amino acids; 2) the alanine-generating Cahill-cycle; and 3) the lactate-releasing Cori-cycle. While the authors have toned down some of their bold conclusions made in the original manuscript, they did very little in the way of providing additional well-controlled experiments, including cell-specific treatments, genetic knockouts, or stable isotope tracing to support their conclusions. Rather, the authors proceed to speculate more without additional data. The major issues raised by this reviewer were not adequately addressed. As such, the conclusions continue to be highly speculative and not well supported with evidence.

Strengths of resubmission:

The resubmission toned down some of its bold statements.

Weaknesses of resubmission:

Major weaknesses of this study persist including lack of in vivo supporting data. Also, retinal explant culture metabolomics are done in neuroretina with RPE attached, which are metabolically active and can be altered by the treatments investigated herein, further confounding the claims made regarding the neuroretina. While including the RPE in the explant model is commended, it needs to be separated from the retina prior to metabolomics to get a better sense of each tissues' metabolism. Also, melanin within RPE will hinder immunofluorescence signal, so one cannot state that RPE do not express certain enzymes based solely on immunofluorescence. Pharmacologic treatments are not cell-specific as the enzymes are expressed in numerous cells within the retina and RPE, and/or the treatments have significant off-target effects (such as shikonin). So, it is difficult to ascertain that the metabolic changes are secondary to the effects on photoreceptors alone, which the authors claim. Additionally, the explants are taken at a very early age when photoreceptors are known to still be maturing. No mention or data is presented on how these metabolic changes are altered in retinal explants after photoreceptors have fully matured. Likewise, significant assumptions are made based on a single metabolomics experiment with no stable isotope tracing to support the pathways suggested. In vivo, stable-isotope retinal metabolomics are being done and have been done, so stating this technology is beyond our field is false. Therefore, the conclusions reached in this manuscript are still not supported.

Reviewer #2 (Public Review):

Summary:

The authors aim to learn about retinal cell specific metabolic pathways, which could substantially improve the way retinal diseases are understood and treated. They culture ex vivo mouse retinas for 6 days with 2 - 4 days of various drug treatments targeting different metabolic pathways or by removing the RPE/choroid tissue from the neural retina. They then look at photoreceptor survival, stain for various metabolic enzymes/transporters and quantify a broad panel of metabolites. While this is an important question to address, the results are not sufficient to support the conclusions.

Strengths:

The questions the authors are exploring at extremely valuable and I commend the authors and working to learn more about retina metabolism. The different sensitivity of the cones to various drugs is interesting and may suggest key differences between rod and cones. The authors also provide a thoughtful discussion of various metabolic pathways in the context of previous publications.

Weaknesses:

As the authors point out, ex vivo culture models allow for control over multiple aspects of the environment (such as drug delivery) not available in vivo. Ex vivo cultures can provide good hints as to what pathways are available between interacting tissues. However, there are many limitations to ex vivo cultures, including shifting to a very artificial culture media condition that is extremely different than the native environment of the retina. It is well appreciated that cells have flexible metabolism and will adapt to conditions provided. Therefore, observations of metabolic responses obtained under culture conditions need to be interpreted with caution, they indicate what the tissue is doing under those specific conditions (which include cells adapting and dying).

Chen et al use pharmacological interventions are to the impact of various metabolic pathways on photoreceptor survival and "long term" metabolic changes. The dose and timing of these drug treatments are not examined though. It is also hard to know how these drugs penetrate the tissue and it needs to be validated that they intended targets are being accurately hit. These relatively long term treatments should be causing numerous downstream changes to metabolism, cell function and survival, which makes looking at a snap shot of metabolite levels hard to interpret. It would be more valuable to look at multiple time points after drug treatment, especially easy time points (closer to 1 hr). the authors use metabolite ratios to make conclusions about pathway activity. It would be more valuable to directly measure pathway activity by looking a metabolite production rates in the media and/or with metabolic tracers again in time scales closer to minutes and hours instead of days.

While the data is interesting and may give insights into some rod and cone specific metabolic susceptibility, more work is needed to validate these conclusions. Given the limitations of the model the authors have over interpreted their findings and the conclusions are not supported by the results. They need to either dramatically limit the scope of their conclusions or validate these hypotheses with additional models and tools.

Reviewer #3 (Public Review):

Summary:

The neural retina is one of the most energetically active tissues in the body and research into retinal metabolism has a rich history. Prevailing dogma in the field is that the photoreceptors of the neural retina (rods and cones) are heavily reliant on glycolysis, and as oxygen tension at the level of photoreceptors is very low, these specialized sensory neurons carry out aerobic glycolysis, akin to the Warburg effect in cancer cells. It has been found that this unique metabolism changes in many retinal diseases, and targeting disease-altered retinal metabolism may be a viable treatment strategy. The neural retina is composed of 11 different cell types, and many research groups over the past century have contributed to our current understanding of cell-specific metabolism of retinal cells. More recently, it has been shown in mouse models and co-culture of the mouse neural retina with human RPE cultures that photoreceptors are reliant on the underlying retinal pigment epithelium for supplying nutrients. Chen and colleagues add to this body of work by studying an ex vivo culture of the developing mouse retina that maintained contact with the retinal pigment epithelium. They exposed such ex vivo cultures to small molecule inhibitors of specific metabolic pathways, performing targeted metabolomics on the tissue and staining tissue with key metabolic enzymes to lay the groundwork for what metabolic pathways may be active in particular cell types of the retina. The authors conclude that rod and cone photoreceptors are reliant on different metabolic pathways to maintain their cell viability - in particular, that rods rely on oxidative phosphorylation and cones rely on glycolysis. Further, their data suggest multiple mechanisms whereby glycolysis may occur simultaneously with anapleurosis to provide abundant energy to photoreceptors. The data from metabolomics revealed several novel findings in retinal metabolism, including the use of glutamine to fuel the mini-Krebs cycle, the utilization of the Cahill cycle in photoreceptors, and a taurine/hypotaurine shuttle between the underlying retinal pigment epithelium and photoreceptors to transfer reducing equivalents from the RPE to photoreceptors. In addition, this study provides quantitative metabolomics datasets that can be compared across experiments and groups. The use of this platform will allow for rapid testing of novel hypotheses regarding the metabolic ecosystem in the neural retina.

Strengths:

The data on differences in susceptibility of rods and cones to mitochondrial dysfunction versus glycolysis provides novel hypothesis-generating conjectures that can be tested in animal models. The multiple mechanisms that allow anapleurosis and glycolysis to run side-by-side add significant novelty to the field of retinal metabolism, setting the stage for further testing of these hypotheses as well.

Weaknesses:

Almost all of the conclusions from the paper are preliminary, based on data showing enzymes necessary for a metabolic process are present and the metabolites for that process are also present. However, to truly prove whether these processes are happening (rather than speculation of the possibility they are happening), further experiments are necessary. As it currently stands, results from this study contradict results from other studies - in particular that cones, not rods, are most reliant of glycolysis. The authors attempt to address these contradictions, but without further experimentation, logical arguments carry only so much weight. At a minimum, the authors have argued that the small molecules they use are exquisitely specific for their intended targets, but validating results with a second small molecule that hits the same target but is structurally different would bolster their claims. Genetically knocking down the intended targets with interfering RNA technology would also be possible, as would explant cultures from knock-out animals. Without these studies to confirm target specificity, combined with the fact that conclusions from this study contradict existing studies in the literature, the results have to be categorized as speculative and hypothesis-generating rather than conclusive.

Joint Public Review:

Bolivar et al. set out to explore whether four distinct neuronal subtypes within the peripheral nervous system exhibit varying potentials for axon regeneration following nerve injury. To investigate this question, they harnessed the power of four distinct reporter mouse models featuring fluorescent labeling of these neuronal subtypes. Their findings reveal that axons of nociceptor neurons exhibit faster regeneration than those of motor neurons, with mechanoreceptors, and proprioceptors displaying the slowest regeneration rate.

To delve into the molecular mechanisms underlying this divergence in regeneration potential, the authors employed the Ribotag technique in mice. This approach enabled them to dissect the differential translatomes of these four neuronal populations after nerve injury, comparing them to uninjured neurons. Their comprehensive expression profiling data uncovers a remarkably heterogeneous response among these neuron subtypes to axon injury.

To focus on one identified target with a mechanistic experiment as a proof of concept, their analysis highlights a striking upregulation of MED12 in proprioceptors, leading to the hypothesis that this molecule may play an inhibitory role, contributing to the comparatively slower regeneration of proprioceptor axons when compared to other neuronal subtypes. This hypothesis gains support from their in vitro model, where siRNA-mediated downregulation of MED12 results in a significant increase in neurite outgrowth in proprioceptive neurons after plating in cell culture dishes.

Overall, this is an interesting study, and in conjunction with similar work from others will be highly valuable for neurobiologists studying how to modulate the regeneration of axons from distinct neuronal subtypes. The quality of data presentation appears to be very good in general, and the manuscript is appropriately written.

Comments on revised version:

Because there are multiple explanations for the differential regeneration responses, the authors have provided further discussion about how regeneration may be regulated in vitro and in vivo. The detection of a gene, Med12, which is unregulated in proprioceptive neurons, but not nociceptive and mechanoceptors, gives support to the existence of specific programs of responses in the peripheral nervous system after injury. Further investigation is needed to define this responsiveness in detail.

Another response is the role of neurotrophins and their receptors. The authors have considered outcomes as a result of different Trk receptor signaling and also the effect of TGFbeta and IL6 as cytokine modulators. Add to this list is the possibility that axon guidance molecules and downstream substates may also play a role.

The original title was considered to be too broad and did not explain all the mechanistic aspects of this study. Therefore a revised title "Neuron-specific RNA-sequencing reveals different response in peripheral neurons after nerve injury" was used. It is appropriately suitable for the results reported in this manuscript.

Reviewer #1 (Public Review):

Summary:

The authors are interested in the developmental origin of the neurons of the cerebellar nuclei. They identify a population of neurons with a specific complement of markers that originate in a distinct location from where cerebellar nuclear precursor cells have been thought to originate that show distinct developmental properties. The cerebellar nuclei have been well studied in recent years both to understand their development and through an evolutionary lens, which supports the importance of this study. The discovery of a new germinal zone giving rise to a new population of CN neurons is an exciting finding, and it enriches our understanding of cerebellar development, which has previously been quite straightforward, where cerebellar inhibitory cells arise from the ventricular zone and the excitatory cells arise from the rhombic lip.

Strengths:

One of the strengths of the manuscript is that the authors use a wide range of technical approaches, including transgenic mice that allow them to disentangle the influence of distinct developmental organizers such as ATOH.<br /> Their finding of a novel germinal zone and a novel population of CN neurons is important for developmental neuroscientists, and cerebellar neuroscientists.

Weaknesses:

One important question raised by this work is what these newly identified cells eventually become in the adult cerebellum. Are they excitatory or inhibitory? Do they correspond to a novel cell type or perhaps one of the cell classes that have been recently identified in the cerebellum (e.g. Fujita et al., eLife, 2020)? Understanding this would significantly bolster the impact of this manuscript.

The major weakness of the manuscript is that it is written for a very specialized reader who has a strong background in cerebellar development, making it hard to read for a general audience. It's challenging to follow the logic of some of the experiments as well as to contextualize these findings in the field of cerebellar development.

Reviewer #2 (Public Review):

Summary:

Canonically cerebellar neurons are derived from 2 primary germinal zones within the anterior hindbrain (dorsal rhombomere 1). This manuscript identifies an important, previously underappreciated origin for a subset of early cerebellar nuclei neurons - the dorsal mesencephalon. This is an exciting finding. While the conclusions are generally supported, several of the figure panels are of inferior quality and do not readily convey the results the authors assert.

Strengths:

The authors have identified a novel early population of cerebellar neurons with likely novel origin in the midbrain. They have used multiple assays to support their conclusions, including immunohistochemistry and in situ analyses of a number of markers of this population which appear to stream from the midbrain into the dorsal anterior cerebellar anlage.

The inclusion of Otx2-GFP short-term lineage analyses and analysis of Atoh1 -/- animals also provide considerable support for the midbrain origin of these neurons as streams of cells seem to emanate from the midbrain. However, without live imaging, there remains the possibility that these streams of cells are not actually migrating, and rather, gene expression is changing in static cells. Hence the authors have conducted midbrain diI labelling experiments of short-term and long-term cultured embryos showing di-labelled cells in the developing cerebellum. These studies confirm the migration of cells from the midbrain into the early cerebellum.

Reviewer #1 (Public Review):

Summary:

This interesting study explores the mechanism behind an increased susceptibility of daf-18/PTEN mutant nematodes to paralyzing drugs that exacerbate cholinergic transmission. The authors use state-of-the-art genetics and neurogenetics coupled with locomotor behavior monitoring and neuroanatomical observations using gene expression reporters to show that the susceptibility occurs due to low levels of DAF-18/PTEN in developing inhibitory GABAergic neurons early during larval development (specifically, during the larval L1 stage). DAF-18/PTEN is convincingly shown to act cell-autonomously in these cells upstream of the PI3K-PDK-1-AKT-DAF-16/FOXO pathway, consistent with its well-known role as an antagonist of this conserved signaling pathway. The authors exclude a role for the TOR pathway in this process and present evidence implicating selectivity towards developing GABAergic neurons. Finally, the authors show that a diet supplemented with a ketogenic body, β-hydroxybutyrate, which also counteracts the PI3K-PDK-1-AKT pathway, promoting DAF-16/FOXO activity, partially rescues the proper development (morphology and function) of GABAergic neurons in daf-18/PTEN mutants, but only if the diet is provided early during larval development. This strongly suggests that the critical function of DAF-18/PTEN in developing inhibitory GABAergic neurons is to prevent excessive PI3K-PDK-1-AKT activity during this critical and particularly sensitive period of their development in juvenile L1 stage worms. Whether or not the sensitivity of GABAergic neurons to DAF-18/PTEN function is a defining and widespread characteristic of this class of neurons in C. elegans and other animals, or rather a particularity of the unique early-stage GABAergic neurons investigated remains to be determined.

Strengths:

The study reports interesting and important findings, advancing the knowledge of how daf-18/PTEN and the PI3K-PDK-1-AKT pathway can influence neurodevelopment, and providing a valuable paradigm to study the selectivity of gene activities towards certain neurons. It also defines a solid paradigm to study the potential of dietary interventions (such as ketogenic diets) or other drug treatments to counteract (prevent or revert?) neurodevelopment defects and stimulate DAF-16/FOXO activity.

Weaknesses:

(1 )Insufficiently detailed methods and some inconsistencies between Figure 4 and the text undermine the full understanding of the work and its implications.

The incomplete methods presented, the imprecise display of Figure 4, and the inconsistency between this figure and the text, make it presently unclear what are the precise timings of observations and treatments around the L1 stage. What exactly do E-L1 and L1-L2 mean in the figure? The timing information is critical for the understanding of the implications of the findings because important changes take place with the whole inhibitory GABAergic neuronal system during the L1 stage into the L2 stage. The precise timing of the events such as neuronal births and remodelling events are well-described (e.g., Figure 2 in Hallam and Jin, Nature 1998; Fig 7 in Mulcahy et al., Curr Biol, 2022). Likewise, for proper interpretation of the implication of the findings, it is important to describe the nature of the defects observed in L1 larvae reported in Figure 1E - at present, a representative figure is shown of a branched commissure. What other types of defects, if any, are observed in early L1 larvae? The nature of the defects will be informative. Are they similar or not to the defects observed in older larvae?

(2) The claim of proof of concept for a reversal of neurodevelopment defects is not fully substantiated by data.

The authors state that the work "constitutes a proof of concept of the ability to revert a neurodevelopmental defect with a dietary intervention" (Abstract, Line 56), however, the authors do not present sufficient evidence to distinguish between a "reversal" or prevention of the neurodevelopment defect by the dietary intervention. This clarification is critical for therapeutic purposes and claims of proof-of-concept. From the best of my understanding, reversal formally means the defect was present at the time of therapy, which is then reverted to a "normal" state with the therapy. On the other hand, prevention would imply an intervention that does not allow the defect to develop to begin with, i.e., the altered or defective state never arises. In the context of this study, the authors do not convincingly show reversal. This would require showing "embryonic" GABAergic neuron defects or showing convincing data in newly hatched L1 (0-1h), which is unclear if they do so or not, as I have failed to find this information in the manuscript. Again, the method description needs to be improved and the implications can be very different if the data presented in Figure 2D-E regard newly born L1 animals (0-1h) or L1 animals at say 5-7h after hatching. This is critical because the development of the embryonically-born GABAergic DD neurons, for instance, is not finalized embryonically. Their neurites still undergo outgrowth (albeit limited) upon L1 birth (see DataS2 in Mulcahy et al., Curr Biol 2022), hence they are susceptible to both committing developmental errors and to responding to nutritional interventions to prevent them. In contrast to embryonic GABAergic neurons, embryonic cholinergic neurons (DA/DB) do not undergo neurite outgrowth post-embryonically (Mulcahy et al., Curr Biol 2022), a fact which could provide some mechanistic insight considering the data presented. However, neurites from other post-embryonically-born neurons also undergo outgrowth post-embryonically, but mostly during the second half of the L1 stage following their birth up to mid-L2, with significant growth occurring during the L1-L2 transition. These are the cholinergic (VA/VB and AS neurons) and GABAergic (VD) neurons. The fact that AS neurons undergo a similar amount of outgrowth as VD neurons is informative if VD neurons are or are not susceptible to daf-18/PTEN activity. Independently, DD neurons are still quite unique on other aspects (see below), which could also bring insight into their selective response.

Finally, even adjusting the claim to "constitutes a proof-of-concept of the ability of preventing a neurodevelpmental defect with a dietary intervention" would not be completely precise, because it is unclear how much this work "constitutes a proof of concept". This is because, unless I misunderstood something, dietary interventions are already applied to prevent neurodevelopment defects, such as when folic acid supplementation is recommended to pregnant women to prevent neural tube defects in newborns.

(3) The data presented do not warrant the dismissal of DD remodeling as a contributing factor to the daf-18/PTEN defects.

Inhibitory GABAergic DD neurons are quite unique cells. They are well-known for their very particular property of remodeling their synaptic polarity (DD neurons switch the nature of their pre- and post-synaptic targets without changing their wiring). This process is called DD remodeling. It starts in the second half of the L1 stage and finishes during the L2 stage. Unfortunately, the fact that the authors find a specific defect in early GABAergic neurons (which are very likely these unique DD neurons) is not explored in sufficient detail and depth. The facts that these neurons are not fully developed at L1, that they still undergo limited neurite growth, and that they are poised for striking synaptic plasticity in a few hours set them apart from the other explored neurons, such as early cholinergic neurons, which show a more stable dynamics and connectivity at L1 (see Mulcahy et al., Curr Biol 2022).

The authors use their observation that daf-18/PTEN mutants present morphological defects in GABAergic neurons prior to DD remodeling to dismiss the possibility that the DAF-18/PTEN-dependent effects are "not a consequence of deficient rearrangement during the early larval stages". However, DD remodeling is just another cell-fate-determined process and as such, its timing, for instance, can be affected by mutations in genes that affect cell fates and developmental decisions, such as daf-18 and daf-16, which affect developmental fates such as those related with the dauer fate. Specifically, the authors do not exclude the possibility that the defects observed in the absence of either gene could be explained by precocious DD remodeling. Precocious DD remodeling can occur when certain pathways, such as the lin-14 heterochronic pathway, are affected. Interestingly, lin-14 has been linked with daf-16/FOXO in at least two ways: during lifespan determination (Boehm and Slack, Science 2005) and in the L1/L2 stages via the direct negative regulation of an insulin-like peptide gene ins-33 (Hristova et al., Mol Cell Bio 2005). It is likely that the prevention of DD dysfunction requires keeping insulin signaling in check (downregulated) in DD neurons in early larval stages, which seems to coincide with the critical timing and function of daf-18/PTEN. Hence, it will be interesting to test the involvement of these genes in the daf-18/daf-16 effects observed by the authors.

Discussion on the impact of the work on the field and beyond:

The authors significantly advance the field by bringing insight into how DAF-18/PTEN affects neurodevelopment, but fall short of understanding the mechanism of selectivity towards GABAergic neurons, and most importantly, of properly contextualizing their findings within the state-of-the-art C. elegans biology.

For instance, the authors do not pinpoint which type of GABAergic neuron is affected, despite the fact that there are two very well-described populations of ventral nerve cord inhibitory GABAergic neurons with clear temporal and cell fate differences: the embryonically-born DD neurons and the post-embryonically-born VD neurons. The time point of the critical period apparently defined by the authors (pending clarifications of methods, presentation of all data, and confirmation of inconsistencies between the text and figures in the submitted manuscript) could suggest that DAF-18/PTEN is required in either or both populations, which would have important and different implications. An effect on DD neurons seems more likely because an image is presented (Figure 2D) of a defect in an L1 daf-18/PTEN mutant larva with 6 neurons (which means the larva was processed at a time when VD neurons were not yet born or expressing pUnc-47, so supposedly it is an image of a larva in the first half of the L1 stage (0-~7h?)). DD neurons are also likely the critical cells here because the neurodevelopment errors are partially suppressed when the ketogenic diet is provided at an "early" L1 stage, but not later (e.g., from L2-L3, according to the text, L2-L4 according to the figure? ).

This study brings important contributions to the understanding of GABAergic neuron development in C. elegans, but unfortunately, it is justified and contextualized mostly in distantly-related fields - where the study has a dubious impact at this stage rather than in the central field of the work (post-embryonic development of C. elegans inhibitory circuits) where the study has stronger impact. This study is fundamentally about a cell fate determination event that occurs in a nutritionally-sensitive developmental stage (post-embryonic L1 larval stage) yet the introduction and discussion are focused on more distantly related problems such as excitatory/inhibitory (E/I) balance, pathophysiology of human diseases, and treatments for them. Whereas speculation is warranted in the discussion, the reduced in-depth consideration of the known biology of these neurons and organisms weakens the impact of the study as redacted. For instance, the critical role of DAF-18/PTEN seems to occur at the early L1 larval stage, a stage that is particularly sensitive to nutritional conditions. The developmental progression of L1 larvae is well-known to be sensitive to nutrition - eg, L1 larvae arrest development in the absence of food, something that is explored in nematode labs to synchronize animals at the L1 stage by allowing embryos to hatch into starvation conditions (water). Development resumes when they are exposed to food. Hence, the extensive postembryonic developmental trajectory that GABAergic neurons need to complete is expected to be highly susceptible to nutrition. Is it? The sensitivity towards the ketogenic diet intervention seems to favor this. In this sense, the attribution of the findings to issues with the nutrition-sensitive insulin-like signaling pathway seems quite plausible, yet this possibility seems insufficiently considered and discussed.

Finally, the fact that imbalances in excitatory/inhibitory (E/I) inputs are linked to Autism Spectrum Disorders (ASD) is used to justify the relevance of the study and its findings. Maybe at this stage, the speculation would be more appropriate if restricted to the discussion. In order to be relevant to ASD, for instance, the selectivity of PTEN towards inhibitory neurons should occur in humans too. However, at present, the E/I balance alteration caused by the absence of daf-18/PTEN in C. elegans could simply be a coincidence due to the uniqueness of the post-embryonic developmental program of GABAergic neurons in C. elegans. To be relevant, human GABAergic neurons should also pass through a unique developmental stage that is critically susceptible to the PI3K-PDK1-AKT pathway in order for DAF-18/PTEN to have any role in determining their function. Is this the case? Hence, even in the discussion, where the authors state that "this study provides universally relevant information on.... the mechanisms underlying the positive effects of ketogenic diets on neuronal disorders characterized by GABA dysfunction and altered E/I ratios", this claim seems unsubstantiated as written particularly without acknowledging/mentioning the criteria that would have to be fulfilled and demonstrated for this claim to be true.

Reviewer #2 (Public Review):

Summary:

Disruption of the excitatory/inhibitory (E/I) balance has been reported in Autism Spectrum Disorders (ASD), with which PTEN mutations have been associated. Giunti et al choose to explore the impact of PTEN mutations on the balance between E/I signaling using as a platform the C. elegans neuromuscular system where both cholinergic (E) and GABAergic (I) motor neurons regulate muscle contraction and relaxation. Mutations in daf-18/PTEN specifically affect morphologically and functionally the GABAergic (I) system, while leaving the cholinergic (E) system unaffected. The study further reveals that the observed defects in the GABAergic system in daf-18/PTEN mutants are attributed to reduced activity of DAF-16/FOXO during development.

Moreover, ketogenic diets (KGDs), known for their effectiveness in disorders associated with E/I imbalances such as epilepsy and ASD, are found to induce DAF-16/FOXO during early development. Supplementation with β-hydroxybutyrate in the nematode at early developmental stages proves to be both necessary and sufficient to correct the effects on GABAergic signaling in daf-18/PTEN mutants.

Strengths:

The authors combined pharmacological, behavioral, and optogenetic experiments to show the GABAergic signaling impairment at the C. elegans neuromuscular junction in DAF-18/PTEN and DAF-16/FOXO mutants. Moreover, by studying the neuron morphology, they point towards neurodevelopmental defects in the GABAergic motoneurons involved in locomotion. Using the same set of experiments, they demonstrate that a ketogenic diet can rescue the inhibitory defect in the daf-18/PTEN mutant at an early stage.

Weaknesses:

The morphological experiments hint towards a pre-synaptic defect to explain the GABAergic signaling impairment, but it would have also been interesting to check the post-synaptic part of the inhibitory neuromuscular junctions such as the GABA receptor clusters to assess if the impairment is only presynaptic or both post and presynaptic.

Moreover, all observations done at the L4 stage and /or adult stage don't discriminate between the different GABAergic neurons of the ventral nerve cord, ie the DDs which are born embryonically and undergo remodeling at the late L1 stage, and VDs which are born post-embryonically at the end of the L1 stage. Those additional elements would provide information on the mechanism of action of the FOXO pathway and the ketone bodies.

Conclusion:

Giunti et al provide fundamental insights into the connection between PTEN mutations and neurodevelopmental defects through DAF-16/FOXO and shed light on the mechanisms through which ketogenic diets positively impact neuronal disorders characterized by E/I imbalances.

Reviewer #3 (Public Review):

Summary:

This is a conceptually appealing study by Giunti et al in which the authors identify a role for PTEN/daf-18 and daf-16/FOXO in the development of inhibitory GABA neurons, and then demonstrate that a diet rich in ketone body β-hydroxybutyrate partially suppresses the PTEN mutant phenotypes. The authors use three assays to assess their phenotypes: (1) pharmacological assays (with levamisole and aldicarb); (2) locomotory assays and (3) cell morphological assays. These assays are carefully performed and the article is clearly written. While neurodevelopmental phenotypes had been previously demonstrated for PTEN/daf-18 and daf-16/FOXO (in other neurons), and while KB β-hydroxybutyrate had been previously shown to increase daf-16/FOXO activity (in the context of aging), this study is significant because it demonstrates the importance of KB β-hydroxybutyrate and DAF-16 in the context of neurodevelopment. Conceptually, and to my knowledge, this is the first evidence I have seen of a rescue of a developmental defect with dietary metabolic intervention, linking, in an elegant way, the underpinning genetic mechanisms with novel metabolic pathways that could be used to circumvent the defects.

Strengths:

What their data clearly demonstrate, is conceptually appealing, and in my opinion, the biggest contribution of the study is the ability of reverting a neurodevelopmental defect with a dietary intervention that acts upstream or in parallel to DAF-16/FOXO.

Weaknesses:

The model shows AKT-1 as an inhibitor of DAF-16, yet their studies show no differences from wildtype in akt-1 and akt-2 mutants. AKT is not a major protein studied in this paper, and it can be removed from the model to avoid confusion, or the result can be discussed in the context of the model to clarify interpretation.

When testing additional genes in the DAF-18/FOXO pathway, there were no significant differences from wild type in most cases. This should be discussed. Could there be an alternate pathway via DAF-18/DAF-16, excluding the PI3K pathway or are there variations in activity of PI3K genes during a ketogenic diet that are hard to detect with current assays?

The consequence of SOD-3 expression in the broader context of GABA neurons was not discussed. SOD-3 was also measured in the pharynx but measuring it in neurons would bolster the claims.

If they want to include AKT-1, seeing its effect on SOD-3 expression could be meaningful to the model.

Reviewer #1 (Public Review):

This study presents valuable observations of white matter organisation from diffusion MRI and two types of synchrotron imaging in both monkeys and mice. Cross-modality comparisons are interesting as the different methods are able to probe anatomical structures at different length scales, from single axons in high-resolution synchrotron (ESRF) imaging, to clusters of axons in lower-resolution synchrotron (DEXY) data, to axon populations at the mm-scale in diffusion MRI. By acquiring all modalities in monkey and mouse ex vivo samples, the authors can observe principles of fibre organisation, and characterise how fibre characteristics, such as tortuosity and micro-dispersion, vary across select brain regions and in healthy tissue versus a demyelination model. The results are solid, though some statements (in the abstract/discussion) do not appear to be fully supported, and statistical tests would help confirm whether tissue characteristics are similar/different between different conditions.

One very interesting result is the observation of apparent laminar organisation of fibres in ex vivo monkey white matter samples. DESY data from the corpus callosum shows fibres with two dominant orientations (one L-R, one slightly inclined), clustered in laminar structures within this major fibre bundle. Thanks to the authors providing open data, I was able to look through the raw DESY volume and observe regions with different "textures" (different orientations) in the described laminar arrangement. That this organisation can be observed by eye, as well as by structure tensor, is fairly convincing. As not all readers will download the data themselves, the manuscript could benefit from additional figures/videos to demonstrate (1) the quality of the DESY data and (2) a more 3D visualisation of the laminar structures (where the coronal plane shows convincing columnar structure or stripes). Similarly in Figure 5A, though this nicely depicts two populations with different orientations, it is somewhat difficult to see the laminar structure in the current image.

ESRF data of the centrum semiovale (CS) contributes evidence for similar laminar structures in a crossing fibre region, where primarily AP fibres are shown to cluster in 3 laminar structures. As above, further visualisations of the ESRF volume in the CS (as shown in Figure 4E) would be of value (e.g. showing consistency across the 4 volumes, 2D images showing stripey/columnar patterns along different axes, etc).

A key limitation of this result is that, though the DESY data from the CC seems convincing, the same structures were not observed in high-resolution synchrotron (ESRF) data of the same tissue sample in the corpus callosum. This seems surprising and the manuscript does not provide a convincing explanation for this inconsistency. The authors argue that this is due to the limited FOV of the ESRF data (~200x200x800 microns). However, the observed laminar structures in DESY are ~40 microns thick, and ERSF data from the CST suggests laminar thicknesses in the range of 5-40 microns with a similar FOV. This suggests the ERSF FOV would be sufficient to capture at least a partial description of the laminar organisation. Further, the DESY data from the CC shows columnar variations along the LR axis, which we might expect to be observed along the long axis of the ESFR volume of the same sample. Additional analyses or explanations to reconcile these apparently conflicting observations would be of value. For example, the authors could consider down-sampling the ESRF data in an appropriate manner to make it more similar to the DESY data, and running the same analysis, to see if the observed differences are related to resolution (i.e. the thinner laminar structures cluster in ways that they look like a thicker laminar structure at lower resolution), or crop the DESY data to the size of the ESRF volume, to test whether the observed differences can be explained by differences in FOV.

Laminar structures were not observed in mouse data, though it is unclear if this is due to anatomical differences or somewhat related to differences in data quality across species.

The authors further quantify various other characteristics of the white matter, such as micro-dispersion, tortuosity, and maximum displacement. Notably, the microscopic FA calculated via structure tensor is fairly consistent across regions, though not modalities. When fibre orientations are combined across the sample, they are shown to produce similar FODs to dMRI acquired in the same tissue, which is reassuring. As noted in the text, the estimates of tortuosity and max displacement are dependent on the FOV over which they are calculated. Calculating these metrics over the same FOV, or making them otherwise invariant to FOV, could facilitate more meaningful comparisons across samples and/or modalities.

Though the results seem solid, some statements, particularly in the abstract and discussion, do not seem to be fully supported by the data. For example, the abstract states "Our findings revealed common principles of fibre organisation in the two species; small axonal fasciculi and major bundles formed laminar structures with varying angles, according to the characteristics of major pathways.", though the results show "no strong indication within the mouse CC of the axonal laminar organisation observed in the monkey". Similarly, the introduction states: "By these means, we demonstrated a new organisational principle of white matter that persists across anatomical length scales and species, which governs the arrangement of axons and axonal fasciculi into sheet-like laminar structures." Further comments on the text are provided below.

One observation not notably discussed in the paper is that the spherical histograms of Figure 3E/H appear to have an anisotropic spread of the white points about 0,0. It would be interesting if the authors could comment on whether this could be interpreted as the FOD having asymmetric dispersion and if so, whether the axis of dispersion relates to the fibre orientations of the laminar structures.

A limitation of the study is that it considers only small ex vivo tissue samples from two locations in a single postmortem monkey brain and slightly larger regions of mouse brain tissue. Consequently, further evidence from additional brain regions and subjects would be required to support more generalised statements about white matter organisation across the brain.

Given the monkey results, the mouse study (section 2.5 onwards) lacks some motivation. In particular, it is unclear why a demyelination model was studied and if/how this would link to the laminar structure observed in the monkey data. Further, it is unclear how comparable tortuosity/max deviation values are across species, considering the differences in data quality and relative resolution, given that the presented results show these values are very modality-dependent.

The paper introduces a new method of "scale-space" parameters for structure tensors. Since, to my understanding, this is the first description of the method, some simple validation of the method would be welcomed. Further, the same scale parameters are not used across monkeys and mice, with a larger kernel used in mice (Table 2) which is surprising given their smaller brain size. Some explanation would be helpful.

Reviewer #2 (Public Review):

Summary:<br /> In this work, the authors combine diffusion MRI and high-resolution x-ray synchrotron phase-contrast imaging in monkey and mouse brains to investigate the 3D organization of brain white matter across different scales and species. The work is at the forefront of the anatomical investigation of the human connectome and aligns with several current efforts to bridge the resolution gap between what we can see in vivo at the millimeter scale and the complexity of the human brain at the sub-micron scale. The authors compare the 3D white matter organization across modalities within 2 small regions in one monkey brain (body of the corpus callosum, centrum semiovale) and within one region (splenium of the corpus callosum) in healthy mice and in one murine model of focal demyelination. The study compares measures of tissue anisotropy and fiber orientations across modalities, performs a qualitative comparison of fasciculi trajectories across brain regions and tissue conditions using streamlined tractography based on the structure tensor, and attempts to quantify the shape of fasciculi trajectories by measuring the tortuosity index and the maximum deviation for each reconstructed streamline. Results show measures of anisotropy and fiber orientations largely agree across modalities, especially for larger FOV data. The high-resolution data allows us to explore the fiber trajectories in relation to tissue complexity and pathology. The authors claim the study reveals new common organization principles of white matter fibers across species and scales, for which axonal fasciculi arrange into sheet-like laminar structures.

Strengths:

The aim of the study is of central importance within present efforts to bridge the gap between macroscopic structures observable in vivo in humans using conventional diffusion MRI and the microscopic organization of white matter tissue. Results obtained from this type of study are important to interpret data obtained in vivo, inform the development of novel methodologies, and expand our knowledge of the structural and thus functional organization of brain circuits.

Multi-scale data acquired across modalities within the same sample constitute extremely valuable data that is often hard to acquire and represent a precious resource for validation of both diffusion MRI tractography and microstructure methods.

The inclusion of multi-species data adds value to the study, allowing the exploration of common organization principles across species.

The addition of data from a murine cuprizone model of focal demyelination adds interesting opportunities to study the underlying biological changes that follow demyelination and how these impact tissue anisotropy and fiber trajectories. These data can inform the interpretation and development of diffusion MRI microstructure models.

Weaknesses:

The main claim of a newly discovered laminar organization principle that is consistent across scales and species is not supported strongly enough by the data. The main evidence in support of the claim comes from the larger FOV data obtained from the body of the corpus callosum in the monkey brain. A laminar organization principle is partially shown in the centrum semiovale in the monkey brain and it is not shown in mice data. Additionally, the methods lack details to help the correct interpretation of these findings (e.g., how were these fasciculi defined?; how well do they represent different axonal populations?; what is the effect of blood vessels on the structure tensor reconstruction?; how was laminar separation quantified?) and the discussion does not provide a biological background for this organization. The corpus callosum sample suggests axons within a bundle of fibers are organized in a sheet-like fashion, while data from the centrum semiovale suggest fibers belonging to different fiber bundles are organized in a sheet-like arrangement. While I acknowledge the challenges in acquiring such high-resolution data, additional samples from different regions in the same animals and from different animals would help strengthen this claim.

The main goal of the study is to bridge the organization of white matter across anatomical length scales and species. However, given the substantial difference in FOVs between the two imaging modalities used, and the absence of intermediate-resolution data, it remains difficult to effectively understand how these results can be used to inform conventional diffusion MRI. In this sense, the introduction does not do a good enough job of building a strong motivation for the scientific questions the authors are trying to answer with these experiments and for the specific methodology used.

The cuprizone data represent a unique opportunity to explore the effect of demyelination on white matter tissue. However, this specific part of the study is not well motivated in the introduction and seems to represent a missed opportunity for further exploration of the qualitative and quantitative relationship between diffusion MRI and sub-micron tissue information (although unfortunately not within the same brain sample). This is especially true considering the diffusion MRI protocol for mice would allow extrapolation of advanced measures from different tissue compartments.

Reviewer #1 (Public Review):

Summary:

The work in the manuscript titled " Altered firing output of VIP interneurons and early dysfunctions in CA1 hippocampal circuits in the 3xTg mouse model of Alzheimer's disease" utilized patch-clamp techniques to explore the electrophysiological characteristics of VIP interneurons in the early stages of AD using the 3xTg mouse model. The study revealed that VIP interneurons exhibited prolonged action potentials and reduced firing rates. These changes could not be attributed to modifications in input signals or morphological transformations. The authors attributed aberrant VIP activity to the accumulation of beta-amyloid in those interneurons.

The decreased frequency of VIP inhibitory events was associated with no observed changes in excitatory drive to these interneurons. Consequently, heightened activity in the general population of CA1 interneurons was observed during a decision-making task and an object recognition test. In light of these findings, the authors concluded that the altered firing patterns of VIP interneurons may initiate early-stage dysfunction in hippocampal CA1 circuits, potentially influencing the progression of AD pathology.

Strengths:

Overall the work is novel and moves the field of Alzheimer's disease forward in a significant way. The manuscript reports a novel concept of aberrant activity in VIP interneurons during the early stages of AD thus contributing to dysfunctions of the CA1 microcircuit. This results in the enhancement of the inhibitory tone on the primary cells of CA1. Thus, the disinhibition by VIP interneurons of Principal Cells is dampened. The manuscript was skillfully composed, and the study was of strong scientific rigor featuring well-designed experiments. Necessary controls were present. Both sexes were included.

Limitations:

(1) The authors attributed aberrant circuit activity to the accumulation of "Abeta intracellularly" inside IS-3 cells. That is problematic. 6E10 antibody recognizes amyloid plaques in addition to Amyloid Precursor Protein (APP) as well as the C99 fragment. There are no plaques at the ages 3xTg mice were examined. Thus, the staining shown in Figure 1a is of APP/C99 inside neurons, not abeta accumulations in neurons. At the ages of 3-6 months, 3xTg starts producing abeta oligomers and potentially tau oligomers as well (Takeda et al., 2013 PMID: 23640054; Takeda et al., 2015 PMID: 26458742 and others). Emerging literature suggests that abeta and tau oligomers disrupt circuit function. Thus, a more likely explanation of abeta and tau oligomers disrupting the activity of VIP neurons is plausible.

(2) Authors suggest that their animals do not exhibit loss of synaptic connections and show Figure 3d in support of that suggestion. However, imaging with confocal microscopy of 70-micron thick sections would not allow the resolution of pre- and post-synaptic terminals. More sensitive measures such as electron microscopy or array tomography are the appropriate techniques to pursue. It is important for the authors to either remove that data from the manuscript or address the limitations of their technique in the discussion section. There is a possibility of loss of synaptic connections in their mouse model at the ages examined.

Reviewer #2 (Public Review):

Summary:

The submitted manuscript by Michaud and Francavilla et al., is a very interesting study describing early disruptions in the disinhibitory modulation exerted by VIP+ interneurons in CA1, in a triple transgenic model of Alzheimer's disease. They provide a comprehensive analysis at the cellular, synaptic, network, and behavioral level on how these changes correlate and might be related to behavioral impairments during these early stages of the disease.

Main findings:

- 3xTg mice show early Aß accumulation in VIP-positive interneurons.

- 3xTg mice show deficits in a spatially modified version of the novel object recognition test.

- 3xTg mice VIP cells present slower action potentials and diminished firing frequency upon current injection.

- 3xTg mice show diminished spontaneous IPSC frequency with slower kinetics in Oriens / Alveus interneurons.

- 3xTg mice show increased O/A interneuron activity during specific behavioral conditions.

- 3xTg mice show decreased pyramidal cell activity during specific behavioral conditions.

Strengths:

This study is very important for understanding the pathophysiology of Alzheimer´s disease and the crucial role of interneurons in the hippocampus in healthy and pathological conditions.

Weaknesses:

Although results nicely suggest that deficits in VIP physiological properties are related to the differences in network activity, there is no demonstration of causality.

Reviewer #1 (Public Review):

Summary:

The paper of Mao et al. expands the genetic toolset that was previously developed by the Rao lab (Denfg et al 2019) to introduce the conditional KO or downregulation of neurotransmission components in Drosophila. The authors then use these tools to investigate neurotransmission in the the clock neurons of the Drosophila brain. They first test some known components and then analyze the contribution of the CNMa neuropeptide and its receptor to the circadian behavior. The results indicate that CNMA acts from a subset of DN1ps (dorsal clock neurons) to set the phase of the morning peak of locomotor activity in light:dark cycles, with an advanced morning activity in the absence of the neuropeptide. Interestingly, the receptor for the PDF neuropeptide appears to be acting in some of the CNMa neurons to control morning activity.

Strengths/weaknesses:

This is clearly a very useful new set of tools to restrict the manipulation of these components to specific neuronal populations, and overall (see specific points below), the paper is convincing to show that the tools indeed allow to efficiently and specifically eliminate neuropeptides/receptors from subsets of neurons. The analysis of the CNMa function in the clock network reveals a new and interesting function for CNMa in the control of morning anticipation in LD conditions. This function appears to depend on CNMA_expressing DN1ps.

Comment on revised version:

I believe that the authors properly addressed the main points that were raised in my comment on version 1.

Reviewer #2 (Public Review):

Original Review:

In this study Mao and co-workers deliver a substantial suite of genetic tools in support of the senior author's recent proposal to create a "chemoconnectomic" tool kit for the expression mapping and conditional disruption of specific neurotransmitter systems with fly neurons of interest. Specifically, they describe the creation of two toolsets for recombination-based and CRISPR/Cas9-based conditional knockouts of genes supporting neurotransmitter and neuromodulator function and Flp-Out and Split-LexA toolkit for the examination of gene expression within defined subsets of neurons. The authors report the creation of conditional genetic tools for the disruption/mapping of approximately 200 chemoconnectomic gene products, an examination of the general effectiveness of these tools in the fly brain and apply them to the circadian clock network in an attempt to reveal new information regarding the transmitter/modulator systems involved in daily behavioral timing. The authors provide clear evidence of the effectiveness of the new methods along with a transparent assessment of the variability of the tools. In addition, they present evidence that the neuro peptide CNMa influences the morning peak of daily activity in the fly by regulating the timing of activity increases in anticipation of dawn.

A major strength of the study is the transparent assessment of the effectiveness and variability of the conditional genetic approaches developed by the authors. The authors have largely achieved their aims and the study therefore represents a major delivery on the promise of chemoconnectomics made by the senior author in 2019 (Neuron, Vol. 101, p. 876). Though there are some concerns about the variability of knockout effectiveness, off target effects of the knockout strategies, and (especially) the accuracy of the gene expression approach, the tools created for this study will almost certainly be useful for the field and support a great deal of future work.

Comments on revised version:

The authors have responded to each of my concerns. Most importantly, they have made the discrepancies within the study and between the study and previously published work clearer to the reader. they have also corrected statements that are not consistent with the current state of the field. The issue regarding opposing effects of PDF signaling and CNMa, which was also raised by Reviewer One still stands, notwithstanding the edits made to the text.

Reviewer #3 (Public Review):

Summary:

Mao and colleagues generated powerful reagents to genetically analyse chemical communication (CCT) in the brain, and in the process uncovered a function for the CNMa neuropeptide expressed in a subset of DN1p neurons that contributes to the temporal organization of locomotor activity, i.e., the timing of morning anticipation.

Strengths:

The strength of the manuscript relies in the generation/characterization of new tools for conditional targeting a well-defined set of CCT genes along with the design and testing of improved versions of Cas9 for efficient knock out. Such invaluable resources will be of interest to the whole community. The authors employed these tools and intersectional genetics to provide an alternative profiling of clock neurons, which is complementary to the ones already published. Furthermore, they uncovered a role for CNMamide, expressed in two DN1ps, in the timing of morning anticipation.

Weaknesses:

All prior concerns have been addressed.

Reviewer #1 (Public Review):

Summary:

In this work, Xie, Prescott and colleagues have reevaluated the role of Nav1.7 in nociceptive sensory neurons excitability. They find that nociceptors can make use of different sodium channel subtypes to reach equivalent excitability. The existence of this degeneracy is critical to understanding the neuronal physiology under normal and pathological conditions and could explain why Nav subtype-selective drugs have failed in clinical trials. More concretely, nociceptor repetitive spiking relies on Nav1.8 at DIV0 (and probably under normal conditions in vivo), but on Nav1.7 and Nav1.3 at DIV4-7 (and after inflammation in vivo).

The conclusions of this paper are mostly well supported by data, and these findings should be of broad interest to scientists working on pain, drug development, neuronal excitability and ion channels.

Strengths:

The authors have employed elegant electrophysiology experiments (including specific pharmacology and dynamic clamp) and computational simulations to study the excitability of a subpopulation of DRGs that would very likely match with nociceptors (they take advantage of using transgenic mice to detect Nav1.8-expressing neurons). They make a strong point showing the degeneracy that occurs at the ion channel expression level in nociceptors, adding this new data to previous observations in other neuronal types. They also demonstrate that the different Nav subtypes functionally overlap and are able to interchange their "typical" roles in action potential generation. As Xie, Prescott and colleagues argue, the functional implications of the degenerate character of nociceptive sensory neurons excitability need to be seriously taken into account regarding drug development and clinical trials with Nav subtype-selective inhibitors.

In this revised version, the quality of the manuscript has been visibly improved. In my opinion, the questions and concerns raised by reviewers have been addressed clearly. After a detailed reading of this version and the comments to the reviewers, I have no additional comments or criticisms.

Reviewer #2 (Public Review):

Summary:

The authors have noted in preliminary work that tetrodotoxin (TTX), which inhibits NaV1.7 and several other TTX-sensitive sodium channels, has differential effects on nociceptors, dramatically reducing their excitability under certain conditions but not under others. Partly because of this coincidental observation, the aim of the present work was to re-examine or characterize the role of NaV1.7 in nociceptor excitability and the effects on drug efficacy. The manuscript demonstrates that a NaV1.7-selective inhibitor produces analgesia only when nociceptor excitability is based on NaV1.7. More generally and comprehensively, the results show that nociceptors can achieve equivalent excitability through changes in differential NaV inactivation and NaV expression of different NaV subtypes (NaV 1.3/1.7 and 1.8). This can cause widespread changes in the role of a particular subtype over time. The degenerate nature of nociceptor excitability shows functional implications that make the assignment of pathological changes to a particular NaV subtype difficult or even impossible.

Thus, the analgesic efficacy of NaV1.7- or NaV1.8-selective agents depends essentially on which NaV subtype controls excitability at a given time point. These results explain, at least in part, the poor clinical outcomes with the use of subtype-selective NaV inhibitors and therefore have major implications for the future development of Nav-selective analgesics.

Strengths:

The results are clearly and impressively supported by the experiments and data shown. During the revision, the manuscript was consistently improved and the concerns of the first reviews were resolved. All methods are described in detail, and presumably, allow good reproducibility and were suitable to address the scientific question.

The results showing that nociceptors can achieve equivalent excitability through changes in differential NaV inactivation and expression of different NaV subtypes are of great importance in the fields of basic and clinical pain research and sodium channel physiology and pharmacology, but also for a broad readership and community. The degenerate nature of nociceptor excitability, which is clearly shown and well supported by data has large functional implications. The results are of great importance because they may explain, at least in part, the poor clinical outcomes with the use of subtype-selective NaV inhibitors and therefore have major implications for the future development of Nav-selective analgesics.

In summary, the authors achieved their overall aim to enlighten the role of the NaV1.7 in nociceptor excitability and the effects on drug efficacy. The data support the conclusions and clinical implications are highlighted as far as is currently justifiable due to the still limited experience in translation. This appears well-considered, not too speculative, and ultimately appropriate.

The main weaknesses of the first version were fixed during the revision:

(i) After revising the manuscript, the initial weakness that the computer model was described superficially has been fixed. Important information was added to the main text and additional information, including the full code and equations and values are deposited on ModelDB or are given in the Supplementary information (Suppl. Table 5 & 6).

(ii) The authors now comment that corresponding studies on protein levels or e.g. neuroinflammatory changes could support the characterization of the time course of membrane expression and cellular changes, but this should be addressed in future studies, as these analyses would also raise new questions, such as about membrane trafficking, post-translational modifications, etc. This is plausible and has now been appropriately addressed in the text.

(iii) During the initial review the authors were asked to discuss the promising role of NaV1.7 in the light of clinical results. In their response, the authors confidently state that they „wish to avoid speculating on which particular clinical results are better explained because our study was not designed for that." They, however, emphasize their take-home message, which is well supported "Instead, our take-home message (which is well supported; see Discussion on lines 309-321) is that NaV1.7-selective drugs may have a variable clinical effect because nociceptors' reliance on NaV1.7 is itself variable - much more than past studies would have readers believe. ... The challenge (as highlighted in the Abstract, lines 21-22) is that identifying the dominant Nav subtype to predict drug efficacy is difficult."

Against the background of this argumentation, it must be admitted that the decision not to present as yet unproven speculations is probably appropriate from a scientific point of view and that this ultimately proves the critical assessment of one's own data and the limitations of the study. This is undoubtedly acceptable and - in retrospect - probably the right way to go.

Reviewer #3 (Public Review):

Summary:

In this study the authors used patch-clamp to characterize the implication of various voltage-gated Na+ channels in the firing properties of mouse nociceptive sensory neurons. They claim that depending on the culture conditions NaV1.3, NaV1.7, and NaV1.8 have distinct contributions to action potential firing and that similar firing patterns can result from distinct relative roles of these channels.

Strengths:

The paper addresses the important issue of understanding the lack of success of therapeutic strategies targeting NaV channels in the context of pain. Specifically, the authors test the hypothesis that different NaV channels contribute in a plastic manner to action potential firing, which may be the reason why it is difficult to target pain by inhibiting these channels.

Weaknesses:

(1) - The main claim of this paper is that "nociceptors can achieve equivalent excitability using different combinations of NaV1.3, NaV1.7, and NaV1.8". From this, they allude to the manifestation of "degeneracy", a concept implying that a biological process can occur via distinct sets of underlying components.<br /> In my opinion, the analyses of the data is biased towards the author's interpretation.<br /> - First, when comparing the excitability across neurons one should relate the response (in this case mean firing frequency) to the absolute size of the stimulus, not to the size of the stimulus normalized to the rheobase (see e.g., Figs. 1A). From this particular figure the authors conclude that the excitability is similar in the culture stages DIV0 and DIV4-7, but these data were not directly compared.<br /> - Second, the authors reach their conclusion from the comparison of the (average) firing rate determined over 1 s current stimulation in distinct conditions. However, this is not the only parameter that determines how sensory neurons might convey information. For instance, the time dependence of the instantaneous frequency, the actual firing pattern, maybe also important.<br /> - Third, the use of 1 s of current stimulation might not be sufficient to characterize the firing pattern if one wants to obtain conclusions that could translate to clinical settings (i.e., sustained pain).<br /> - Fourth, out of principle, the gating properties of NaV1.7 and NaV1.8 channels are not identical, and therefore their contributions to excitability should not be the same. A neuron in which NaV1.7 is the main contributor is expected to have a damping firing pattern due to cumulative channel inactivation, whereas another depending mainly on NaV1.8 is expected to display more sustained firing. This is actually seen in the results of the modelling.

(2) - The quality of some recordings is dubious. The currents shown as TTX-sensitive in Fig. 1D look very strange (not like the ones at Baseline DIV4-7). These traces show abnormally fast inactivation and even transient deflections above zero current line. These are obvious artifacts of the subtraction procedure, probably due to unstable current amplitudes along the recording time. Similar odd-looking traces are shown in Fig. 3A.

(3) - I would like to point out that the main Significance Statement of the manuscript reads "The analgesic efficacy of subtype-selective drugs hinges on which subtype controls excitability". I would like to point out that, in addition of being extremely obvious for anyone knowing a bit about pain signaling, the authors did not test the analgesic efficacy of any drug in this study.

(4) - A critical issue in the manuscript is the unnecessary use of phrases that imply that biological entities have some sort of willpower, flirting with anthropomorphism and teleological language.<br /> Sentences such as "Nociceptive sensory neurons convey pain signals to the CNS using action potentials" (see the Abstract) should be avoided. Neurons do not really "use" action potentials, they have no will to do so. Action potentials are not tools or means to be "used" by neurons. There are many other examples of misuse of the verb "use" in many other sentences. These were pointed out during the revision phase, but unfortunately the authors refused to correct them.

Reviewer #1 (Public Review):

Summary:

In Ryu et al., the authors use a cortical mouse astrocyte culture system to address the functional contribution of astrocytes to circadian rhythms in the brain. The authors' starting point is transcriptional output from serum-shocked culture, comparative informatics with existing tools and existing datasets. After fairly routine pathway analyses, they focus on the calcium homeostasis machinery and one gene, Herp, in particular. They argue that Herp is rhythmic at both mRNA and protein levels in astrocytes. They then use a calcium reporter targeted to the ER, mitochondria, or cytosol and show that Herp modulates calcium signaling as a function of circadian time. They argue that this occurs through the regulation of inositol receptors. They claim that the signaling pathway is clock-controlled by a limited examination of Bmal1 knockout astrocytes. Finally, they switch to calcium-mediated phosphorylation of the gap junction protein Connexin 43 but do not directly connect HERP-mediated circadian signaling to these observations. While these experiments address very important questions related to the critical role of astrocytes in regulating circadian signaling, the mechanistic arguments for HERP function, its role in circadian signaling through inositol receptors, the connection to gap junctions, and ultimately, the functional relevance of these findings is only partially substantiated by experimental evidence.

Strengths:

- The paper provides useful datasets of astrocyte gene expression in circadian time.

- Identifies HERP as a rhythmic output of the circadian clock.

- Demonstrates the circadian-specific sensitivity of ATP -> calcium signaling.

- Identifies possible rhythms in both Connexin 43 phosphorylation and rhythmic movement of calcium between cells.

Weaknesses:

- It is not immediately clear why the authors chose to focus on Ca2+ homeostasis or Herp from their initial screens as neither were the "most rhythmic" pathways in their primary analyses.

- It would have been interesting (and potentially important) to know whether various methods of cellular synchronization would also render HERP rhythmic (e.g., temperature, forskolin, etc). If Herp is indeed relatively astrocyte-specific and rhythmic, it should be easy to assess its rhythmicity in vivo.

- The authors show that Herp suppression reduces ATP-mediated suppression of calcium whereas it initially increases Ca2+ in the cytosol and mitochondria and then suppresses it. The dynamics of the mitochondrial and cytosolic responses are not discussed in any detail and it is unclear what their direct relationship is to Herp-mediated ER signaling. What is the explanation for Herp (which is thought to be ER-specific) to calcium signaling in other organelles?

- What is the functional significance of promoting ATP-mediated suppression of calcium in ER?

- The authors then nicely show that the effect of ATP is dependent on intrinsic circadian timing but do not explain why these effects are antiphase in cytosol or mitochondria. Moreover, the ∆F/F for calcium in mitochondria and cytosol both rise, cross the abscissa, and then diminish - strongly suggesting a biphasic signaling event. Therefore, one wonders whether measuring the area under the curve is the most functionally relevant measurement of the change.

- Why are mitochondrial and cytosolic calcium not also demonstrated for Bmal1 KO astrocytes?

- The authors claim that Herp acts by regulating the degradation of ITPRs but this hypothesis - rather central to the mechanisms proposed in this study - is not experimentally substantiated.

- There is no clear demonstration of the functional relevance of the circadian rhythms of ATP-mediated calcium signaling.

Reviewer #2 (Public Review):

Summary:

The article entitled "Circadian regulation of endoplasmic reticulum calcium response in mouse cultured astrocytes" submitted by Ryu and colleagues describes the circadian control of astrocytic intracellular calcium levels in vitro.

Strengths:

The authors used a variety of technical approaches that are appropriate

Weaknesses:

Statistical analysis is poor and could lead to a misinterpretation of the data

Several conceptual issues have been identified.

Overinterpretation of the data should be avoided. This is a mechanistic paper done completely in vitro, all references to the in vivo situation are speculative and should be avoided.

Reviewer #3 (Public Review):

Astrocyte biology is an active area of research and this study is timely and adds to a growing body of literature in the field. The RNA-seq, Herp expression, and Ca2+ release data across wild-type, Bmal1 knockout, and Herp knockdown cellular models are robust and lend considerable support to the study's conclusions, highlighting their importance. Despite these strengths, the manuscript presents a gap in elucidating the dynamics of HERP and the involvement of ITPR1/2 in modulating Ca2+ release patterns and their circadian variations, which remains insufficiently supported and characterized. While the Connexin data underscore the importance of rhythmic Ca2+ release triggered by ATP, the relationship here appears correlational and the role of HERP and ITPR in Cx function remains to be characterized. Moreover, enhancing the manuscript's clarity and readability could significantly benefit the presentation and comprehension of the findings.

Reviewer #1 (Public Review):

This is an interesting and well-written paper reporting on a novel approach to studying cerebellar function based on the idea of selective recruitment using fMRI. The study is well-designed and executed. Analyses are sound and results are properly discussed. The paper makes a significant contribution to broadening our understanding of the role of the cerebellum in human behavior.

- While the authors provide a compelling case for the link between BOLD and the cerebellar cortical input layer, there remains considerable unexplained variance. Perhaps the authors could elaborate a bit more on the assumption that BOLD signals mainly reflect the input side of the cerebellum (see for example King et al., elife. 2023 Apr 21;12:e81511).

- The current approach does not appear to take the non-linear relationships between BOLD and neural activity into account.

- The authors may want to address a bit more the issue of closed loops as well as the underlying neuroanatomy including the deep cerebellar nuclei and pontine nuclei in the context of their current cerebello-cortical correlational approach. But also the contribution of other brain areas such as the basal ganglia and hippocampus.

- What about the direct projections of mossy fibers to the DCN that actually bypasses the cerebellar cortex?

Reviewer #2 (Public Review):

Summary:

Shahshahani and colleagues used a combination of statistical modelling and whole-brain fMRI data in an attempt to separate the contributions of cortical and cerebellar regions in different cognitive contexts.

Strengths:

* The manuscript uses a sophisticated integration of statistical methods, cognitive neuroscience, and systems neurobiology.

* The authors use multiple statistical approaches to ensure robustness in their conclusions.

* The consideration of the cerebellum as not a purely 'motor' structure is excellent and important.

Weaknesses:

* Two of the foundation assumptions of the model - that cerebellar BOLD signals reflect granule cells > purkinje neurons and that corticocerebellar connections are relatively invariant - are still open topics of investigation. It might be helpful for the reader if these ideas could be presented in a more nuanced light.

* The assumption that cortical BOLD responses in cognitive tasks should be matched irrespective of cerebellar involvement does not cohere with the idea of 'forcing functions' introduced by Houk and Wise.

Reviewer #1 (Public Review):

In this work, the authors set out to ask whether the MYRF family of transcription factors, represented by myrf-1 and myrf-2 in C. elegans, have a role in the temporally controlled expression of the miRNA lin-4. The precisely timed onset of lin-4 expression in the late L1 stage is known to be a critical step in the developmental timing ("heterochronic") pathway, allowing worms to move from the L1 to the L2 stage of development. Despite the importance of this step of the pathway, the mechanisms that control the onset of lin-4 expression are not well understood.

Overall, the paper provides convincing evidence that MYRF factors have a key role in promoting lin-4 expression in young larvae. Using state-of-the-art techniques (knock-in reporters and conditional alleles), the authors show that MYRF factors are essential for lin-4 activation and act cell-autonomously. Results using some unusual gain-of-function alleles are supported by consistent results using other approaches. The authors also provide evidence supporting the idea that MYRF factors activate lin-4 by directly activating its promoter. Because these results are indirect test of this, further experiments will be necessary to conclusively determine whether lin-4 is indeed a direct target of MYRF factors. myrf-1 and myrf-2 likely function redundantly to activate lin-4; potential complex interactions between these two genes will be an interesting area for future work.

Overall, the paper's results are convincing. The important findings on miRNA regulation and the control of developmental timing will make this work of interest to a broad range of developmental biologists.

Reviewer #2 (Public Review):

Summary:

In this manuscript, the authors examine how temporal expression of the lin-4 microRNA is transcriptionally regulated.

Comments on revised version:

In the revised manuscript, the authors have suitably addressed my original concerns.

Aims achieved: The aims of the work are now achieved.

Impact: This study shows that a single transcription factor (MYRF-1) is important for the regulation of multiple microRNAs that are expressed early in development to control developmental timing.

Reviewer #2 (Public Review):

Summary:

Larouche et al show that TEs are broadly expressed in thymic cells, especially in mTECs and pDCs. Their data suggest a possible involvement of TEs in thymic gene regulation and IFN-alpha secretion. They also show that at least some TE-derived peptides are presented by MHC-I in the thymus.

Strengths:

The idea of high/broad TE expression in the thymus as a mechanism for preventing TE-mediated autoimmunity is certainly an attractive one, as is their involvement in IFN-alpha secretion therein. The analyses and experiments presented here are therefore a very useful primer for more in-depth experiments, as the authors point out towards the end of the discussion.

Weaknesses:

There are many dangers about analysing RNA-seq data at the subfamily level. Outputs may be greatly confounded by pervasive transcription, DNA contamination, and overlap of TEs with highly expressed genes. Whether TE transcripts are independent units or part of a gene also has important implications for the conclusions drawn. The authors have tried to mitigate against some of these issues, but they have not been completely ruled out.

Reviewer #1 (Public Review):

Summary:

Transposable Elements (TEs) are exogenously acquired DNA regions that have played important roles in the evolutional acquisition of various biological functions. TEs may have been important in the evolution of the immune system, but their role in thymocytes has not been fully clarified.

Using the human thymus scRNA dataset, the authors suggest the existence of cell type-specific TE functions in the thymus. In particular, it is interesting to show that there is a unique pattern in the type and expression level of TEs in thymic antigen-presenting cells, such as mTECs and pDCs, and that they are associated with transcription factor activities. Furthermore, the authors suggested that TEs may be non-redundantly regulated in expression by Aire, Fezf2, and Chd4, and that some TE-derived products are translated and present as proteins in thymic antigen-presenting cells. These findings provide important insights into the evolution of the acquired immune system and the process by which the thymus acquires its function as a primary lymphoid tissue.

Strengths:

(1) By performing single-cell level analysis using scRNA-seq datasets, the authors extracted essential information on heterogeneity within the cell population. It is noteworthy that this revealed the diversity of expression not only of known autoantigens but also of TEs in thymic antigen-presenting cells.

(2) The attempt to use mass spectrometry to confirm the existence of TE-derived peptides is worthwhile, even if the authors did not obtain data on as many transcripts as expected.

(3) The use of public data sets and the clearly stated methods of analysis improved the transparency of the results.

Weaknesses:

(1) The authors sometimes made overstatements largely due to the lack or shortage of experimental evidence.

For example in Figure 4, the authors concluded that thymic pDCs produced higher copies of TE-derived RNAs to support the constitutive expression of type-I interferons in thymic pDCs, unlike peripheral pDCs. However, the data was showing only the correlation between the distinct TE expression pattern in pDCs and the abundance of dsRNAs. We are compelled to say that the evidence is totally too weak to mention the function of TEs in the production of interferon. Even if pDCs express a distinct type and amount of TE-derived transcripts, it may be a negligible amount compared to the total cellular RNAs. How many TE-derived RNAs potentially form the dsRNAs? Are they over-expressed in pDCs?<br /> The data interpretation requires more caution to connect the distinct results of transcriptome data to the biological significance.

We contend that our manuscript combines the attributes of a research article (novel concepts) and a resource article (datasets of TEs implicated in various aspects of thymus function). The critical strength of our work is that it opens entirely novel research perspectives. We are unaware of previous studies on the role of TEs in the human thymus. The drawback is that, as with all novel multi-omic systems biology studies, our work provides a roadmap for a multitude of future mechanistic studies that could not be realized at this stage. Indeed, we performed wet lab experiments to validate some but not all conclusions: i) presentation of TE-derived MAPs by TECs and ii) formation of dsRNAs in thymic pDCs. In response to Reviewer #1, we performed supplementary analyses to increase the robustness of our conclusions. Also, we indicated when conclusions relied strictly on correlative evidence and clarified the hypotheses drawn from our observations. Regarding the Reviewer's questions about TE-derived dsRNAs, LINE, LTR, and SINE elements all have the potential to generate dsRNAs, given their highly repetitive nature and bi-directional transcription (1). As ~32% of TE subfamilies are overexpressed in pDCs, we hypothesized that these TE sequences might form dsRNA structures in these cells. To address the Reviewer's concerns regarding the amount of TE-derived RNAs among total cellular RNAs, we also computed the percentage of reads assigned to TEs in the different subsets of thymic APCs (see Reviewer 1 comment #4).<br /> ------

I appreciate the authors' efforts to improve the quality of this valuable paper. The additional data proposed by the authors enhanced the possibility that the non-negligible amount of RNAs in pDCs is derived from TE elements. Their biological roles and significance will be demonstrated in future research.

(2) Lack of generality of specific examples. This manuscript discusses the whole genomic picture of TE expression. In addition, one good way is to focus on the specific example to clearly discuss the biological significance of the acquisition of TEs for the thymic APC functions and the thymic selection.

In Figure 2, the authors focused on ETS-1 and its potential target genes ZNF26 and MTMR3, however, the significance of these genes in NK cell function or development is unclear. The authors should examine and discuss whether the distinct features of TEs can be found among the genomic loci that link to the fundamental function of the thymus, e.g., antigen processing/presentation.

We thank the Reviewer for this highly relevant comment. We investigated the genomic loci associated with NK cell biology to determine if ETS1 peaks would overlap with TE sequences in protein-coding genes' promoter region. Figure 2h illustrates two examples of ETS1 significant peaks overlapping TE sequences upstream of PRF1 and KLRD1. PRF1 is a protein implicated in NK cell cytotoxicity, whereas KLRD1 (CD94) dimerizes with NKG2 and regulates NK cell activation via interaction with the nonclassical MHC-I molecule HLA-E (2, 3). Thus, we modified the section of the manuscript addressing these results to include these new analyses: "Finally, we analyzed publicly available ChIP-seq data of ETS1, an important TF for NK cell development (4), to confirm its ability to bind TE sequences. Indeed, 19% of ETS1 peaks overlap with TE sequences (Figure 2g). Notably, ETS1 peaks overlapped with TE sequences (Figure 2h, in red) in the promoter regions of PRF1 and KLRD1, two genes important for NK cells' effector functions (2, 3)."<br /> ------

I am convinced by the authors' explanation that TE elements may contribute to the functions of NK cells.<br /> However, since I have understood that the main topic of this paper is about the thymus and thymic antigen-presenting cells, the mention of NK cells seems abrupt and unconnected to me. NK cells are a type of innate lymphocyte that arise in the bone marrow, and thymus is dispensable for their development and function. The readers might expect to find something more fundamental regarding the function of the thymus and immunological tolerance.

(3) Since the deep analysis of the dataset yielded many intriguing suggestions, why not add a discussion of the biological reasons and significance? For example, in Figure 1, why is TE expression negatively correlated with proliferation? cTEC-TE is mostly postnatal, while mTEC-TE is more embryonic. What does this mean?

We thank the Reviewer for this comment. To our knowledge, the relationship between cell division and transcriptional activity of TEs has not been extensively studied in the literature. However, a recent study has shown that L1 expression is induced in senescent cells. We therefore added the following sentences to our Discussion: "The negative correlation between TE expression and cell cycle scores in the thymus is coherent with recent data showing that transcriptional activity of L1s is increased in senescent cells (5). A potential rationale for this could be to prevent deleterious transposition events during DNA replication and cell division." We also added several discussion points regarding the regulation of TEs by KZFPs to answer concerns raised by Reviewer 2 (see Reviewer 2 comment #1).<br /> ------

I agree on the possibility suggested by the authors.

(4) To consolidate the experimental evidence about pDCs and TE-derived dsRNAs, one option is to show the amount of TE-derived RNA copies among total RNAs. The immunohistochemistry analysis in Figure 4 requires additional data to demonstrate that overlapped staining was not caused by technical biases (e.g. uneven fixation may cause the non-specifically stained regions/cells). To show this, authors should have confirmed not only the positive stainings but also the negative staining (e.g. CD3, etc.). Another possible staining control was showing that non-pDC (CD303- cell fractions in this case) cells were less stained by the ds-RNA probe.

We thank the Reviewer for this suggestion. We computed the proportion of reads in each cell assigned to two groups of sequences known to generate dsRNAs: TEs and mitochondrial genes (1). These analyses showed that the proportion of reads assigned to TEs is higher in pDCs than other thymic APCs by several orders of magnitude (~20% of all reads). In contrast, reads derived from mitochondrial genes had a lower abundance in pDCs. We included these results in Figure 4 - figure supplement 2 and included the following text in the Results section "To evaluate if these dsRNAs arise from TE sequences, we analyzed in thymic APC subsets the proportion of the transcriptome assigned to two groups of genomic sequences known as important sources of dsRNAs, TEs and mitochondrial genes (1). Strikingly, whereas the percentage of reads from mitochondrial genes was typically lower in pDCs than in other thymic APCs, the proportion of the transcriptome originating from TEs was higher in pDCs (~22%) by several orders of magnitude (Figure 4 - figure supplement 2)." As a negative control for the immunofluorescence experiments, we used CD123- cells. Indeed, flow cytometry analysis of the magnetically enriched CD303+ fraction was around 90% pure, as revealed by double staining with CD123 and CD304 (two additional markers of pDCs): CD123- cells were also CD304-/lo, showing that these cells are non- pDCs. Thus, we decided to compare the dsRNA signal between CD123+ cells (pDCs) and CD123- cells (non-pDCs). The difference between CD123+ and CD123- cells was striking (Figure 4d).<br /> ------

Although the technical concerns about immunostaining were not resolved, it is understandable that it would be difficult to rerun the experiment since the authors used the precious human thymi as the experimental material. Immunostaining co-staining requires careful interpretation so that careful experimental setup is needed.

Reviewer #2 (Public Review):

Summary:

Shotgun data have been analysed to obtain fungal and bacterial organisms abundance. Through their metabolic functions and through co-occurrence networks, a functional relationship between the two types of organisms can be inferred. By means of metabolomics, function-related metabolites are studied in order to deepen the fungus-bacteria synergy.

Strengths:

Data obtained in bacteria correlate with data from other authors.<br /> The study of metabolic "interactions" between fungi and bacteria is quite new.<br /> The inclusion of metabolomics data to support the results is a great contribution.

Weaknesses:

All my concerns have been clarified

Reviewer #1 (Public Review):

The manuscript by Mullen et al. investigated the gene expression changes in cancer cells treated with the DHODH inhibitor brequinar (BQ), to explore the therapeutic vulnerabilities induced by DHODH inhibition. The study found that BQ treatment causes upregulation of antigen presentation pathway (APP) genes and cell surface MHC class I expression, mechanistically which is mediated by the CDK9/PTEFb pathway triggered by pyrimidine nucleotide depletion. The combination of BQ and immune checkpoint therapy demonstrated a synergistic (or additive) anti-cancer effect against xenografted melanoma, suggesting the potential use of BQ and immune checkpoint blockade as a combination therapy in clinical therapeutics.

The interesting findings in the present study include demonstrating a novel cellular response in cancer cells induced by DHODH inhibition. However, whether the increased antigen presentation by DHODH inhibition actually contributed to the potentiation of the efficacy of immune-check blockade (ICB) is not directly examined is the limitation of the study. Moreover, the mechanism of the increased antigen presentation pathway by pyrimidine depletion mediated by CDK9/PTEFb was not validated by genetic KD or KO targeting by CDK9/PTEFb pathways. Finally, high concentrations of BQ have been reported to show off-target effects, sensitizing cancer cells to ferroptosis, and the authors should discuss whether the dose used in the in vivo study reached the ferroptotic sensitizing dose or not.

Comment on the revised version:

In their response letter, the authors appropriately addressed the reviewer's comments.

However, it is unfortunate that these comments are not reflected in the main text. Consequently, readers may encounter the same questions. Therefore, the reviewer recommends mentioning them in the discussion or limitations of the study, even if briefly, to address readers' concerns. Especially, addressing the comments such as the dosage of BQ being lower than the reported pro-ferroptotic dose (PMID 37407687), and the lack of examining potential impact of immune cell depletion on the efficacy of BQ treatment would be necessary for considering the proposed mechanism. The latter limitation is also raised by the other reviewer.

Reviewer #2 (Public Review):

In their manuscript entitled "DHODH inhibition enhances the efficacy of immune checkpoint blockade by increasing cancer cell antigen presentation", Mullen et al. describe an interesting mechanism of inducing antigen presentation. The manuscript includes a series of experiments that demonstrate that blockade of pyrimidine synthesis with DHODH inhibitors (i.e. brequinar (BQ)) stimulates the expression of genes involved in antigen presentation. The authors provide evidence that BQ mediated induction of MHC is independent of interferon signaling. A subsequent targeted chemical screen yielded evidence that CDK9 is the critical downstream mediator that induces RNA Pol II pause release on antigen presentation genes to increase expression. Finally, the authors demonstrate that BQ elicits strong anti-tumor activity in vivo in syngeneic models, and that combination of BQ with immune checkpoint blockade (ICB) results in significant lifespan extension in the B16-F10 melanoma model. Overall, the manuscript uncovers an interesting and unexpected mechanism that influences antigen presentation and provides an avenue for pharmacological manipulation of MHC genes, which is therapeutically relevant in many cancers. However, a few key experiments are needed to ensure that the proposed mechanism is indeed functional in vivo.

Major Points:

(1) According to the proposed model, BQ mediated induction of antigen presentation is a contributing factor to the efficacy of this therapeutic strategy. If this is true, then depletion of immune cells should reduce the therapeutic efficacy of BQ in vivo. The authors should perform the B16-F10 transplant experiments in either Rag null mice (if available) or with CD8/CD4 depletion. The expectation would be that T cell depletion (or MHC loss with genetic manipulation) should reduce the efficacy of BQ treatment. Absent this critical experiment, it is difficult to confidently conclude that induction of antigen presentation is a fundamental component of the in vivo response to DHODH inhibition.

(2) Does BQ treatment induce antigen presentation in non-malignant cells? APCs? If the induction of antigen presentation is not cancer specific and related to a pyrimidine depletion stress response, then there is a possibility that healthy tissues will also exhibit a similar phenotype, raising concerns about the specificity of a de novo immune response. The authors should examine antigen presentation genes in healthy tissues treated with BQ.

(3) In the title, the authors claim that DHODH enhances the efficacy of ICB. However, the experiment shown in Figure 5D does not demonstrate this. The Kaplan Meier curves reflect more of an additive response versus a synergistic combination. Furthermore, the concurrent treatment of BQ and ICB seems to inhibit the efficacy of ICB due to BQ toxicity in immune cells. When concurrently administered, the survival of the mice is the same as with brequinar alone, suggesting that the efficacy of ICB was diminished. However, if ICB is administered following an initial dose of BQ, there is an added survival benefit of a magnitude that is similar to ICB alone. This result seems to contradict the title. Furthermore, the authors should show the longitudinal growth curves of these tumors.

(4) Related to Point 3, the temporal separation of BQ and ICB raises the question of whether the induction of antigen presentation with BQ is persistent during the course of delayed ICB treatment. One explanation for the results is that BQ treatment reduces tumor burden, and then a subsequent course of ICB also reduces tumor burden but not that the two therapies are functioning in synergy. To address this, the authors should measure the duration of BQ mediated induction of antigen presentation after stopping treatment.

(5) In Figure 1, the authors show that DHODH inhibition induces expression of both MHC-I and MHC-II genes at the RNA level. However, they only validate MHC-I by flow cytometry. A simple experiment to evaluate the effect of BQ treatment on MHC-II surface expression would provide important additional mechanistic insight into the immunomodulatory effects of DHODH inhibition, especially given recent literature reinforcing the importance of MHC-II expression on epithelial cancers, including melanoma (Oliveira et al. Nature 2022).

Minor Points:

(1) The authors show ChIP-seq tracks from Tan et al. for HLA-B. However, given the pervasive effect of Ter treatment across many HLA genes, the authors should either show tracks at additional loci, or provide a heatmap of read density across more loci. This would substantiate the mechanistic claim that RNA Pol II occupancy and activity across antigen presentation genes is the major driver of response to DHODH inhibition as opposed to mRNA stabilization/increased translation.

(2) A compelling way to demonstrate a change in antigen presentation is through mass spectrometry based immunopeptidomics. Performing immunopeptidomic analysis of BQ treated cell lines would provide substantial mechanistic insight into the outcome of BQ treatment. While this approach may be outside the scope of the current work, the authors should speculate on how this treatment may specifically alter the antigenic landscape where future directions would include empirical immunopeptidomics measurements.

(3) While the signaling through CDK9 seems convincing, it still does not provide a mechanistic link between depleted pyrimidines and CDK9 activity. The authors should speculate on the mechanism that signals to CDK9.

(4) Related to minor point 2, the authors should consider a genetic approach to confirm the importance of CDK9. While the pharmacological approach, including multiple mechanistically distinct CDK9 inhibitors provides strong evidence, an additional experiment with genetic depletion of CDK9 (CRISPR KO, shRNA, etc) would provide compelling mechanistic confirmation.

(5) The authors should comment in the discussion on how this strategy may be particularly useful in patients harboring genetic or epigenetic loss of interferon signaling, a known mechanism of ICB resistance. Perhaps DHODH inhibition could rescue MHC expression in cells that are deficient in interferon sensing.

Overall, the paper is clearly written and presented. With the additional experiments described above, especially in vivo, this manuscript would provide a strong contribution to the field of antigen presentation in cancer. The distinct mechanisms by which DHODH inhibition induces antigen presentation will also set the stage for future exploration into alternative methods of antigen induction.

Comments on latest version:

The authors address the majority of the points raised in my previous review. However, no additional in vivo experiments were performed, which seems necessary for the major conclusions of the paper.

I disagree with the authors' assessment of Major Point 3 in my review. I have updated the text of Major Point 3 in my public review to further clarify my position.

My final assessment is that if the authors want to claim that DHODH inhibition potentiates immune checkpoint blockade, as is stated in the title, then further in vivo experimentation is needed.

Reviewer #3 (Public Review):

Mullen et al present an important study describing how DHODH inhibition enhances efficacy of immune checkpoint blockade by increasing cell surface expression of MHC I in cancer cells. DHODH inhibitors have been used in the clinic for many years to treat patients with rheumatoid arthritis and there has been a growing interest in repurposing these inhibitors as anti-cancer drugs. In this manuscript, the Singh group builds on their previous work defining combinatorial strategies with DHODH inhibitors to improve efficacy. The authors identify an increased expression of genes in the antigen presentation pathway and MHC I after BQ treatment which is mediated strictly by pyrimidine depletion and CDK9/P-TEFb. The authors rationalize that increased MHC I expression induced by DHODH inhibition might favor efficacy of dual immune checkpoint blockade. In fact, this combinatorial treatment prolonged survival in an immunocompetent B16F10 melanoma model.

Previous studies have shown that DHODH inhibitors can increase expression of innate immunity-related genes but the role of DHODH and pyrimidine nucleotides in antigen presentation has not been previously reported. A strength of the manuscript is the solid in vitro mechanistic data supported by analysis in multiple cell lines. The in vivo data show compelling additive effects of DHODH inhibitors and ICB. However, more controls and experiments would be required to define the nature of these effects and to confirm that the mechanistic in vitro data is conserved in vivo.

This is a relevant manuscript proposing a mechanistic link between pyrimidine depletion and MHC I expression and a novel therapeutic approach combining DHODH inhibitors with dual checkpoint blockade. These results might be relevant for the clinical development of DHODH inhibitors in the treatment of solid tumors, a setting where these have not shown optimal efficacy yet.

Comments on revised version:

The authors have addressed my questions regarding validation of gene expression in other cell lines. They have also provided an explanation about why in vivo evaluations could not be performed for the experiment in Figure 5E.

Reviewer #1 (Public Review):

Summary:

This manuscript by Warfvinge et al. reports the results of CITE-seq to generate single-cell multi-omics maps from BM CD34+ and CD34+CD38- cells from nine CML patients at diagnosis. Patients were retrospectively stratified by molecular response after 12 months of TKI therapy using European Leukemia Net (ELN) recommendations. They demonstrate heterogeneity of stem and progenitor cell composition at diagnosis, and show that compared to optimal responders, patients with treatment failure after 12 months of therapy demonstrate increased frequency of molecularly defined primitive cells at diagnosis. These results were validated by deconvolution of an independent previously published dataset of bulk transcriptomes from 59 CML patients. They further applied a BCR-ABL-associated gene signature to classify primitive Lin-CD34+CD38- stem cells as BCR:ABL+ and BCR:ABL-. They identified variability in the ratio of leukemic to non-leukemic primitive cells between patients, showed differences in expression of cell surface markers and determined that a combination of CD26 and CD35 cell surface markers could be used to prospectively isolate the two populations. The relative proportion of CD26-CD35+ (BCR:ABL-) primitive stem cells was higher in optimal responders compared to treatment failures, both at diagnosis and following 3 months of TKI therapy.

Strengths:

The studies are carefully conducted and the results are very clearly presented. The data generated will be a valuable resource for further studies. The strengths of this study are the application of single-cell multi-omics using CITE-Seq to study individual variations in stem and progenitor clusters at diagnosis that are associated with good versus poor outcomes in response to TKI treatment. These results were confirmed by deconvolution of a historical bulk RNAseq data set. Moreover, they are also consistent with a recent report from Krishnan et al. and are a useful confirmation of those results. The major new contribution of this study is the use of gene expression profiles to distinguish BCR-ABL+ and BCR-ABL- populations within CML primitive stem cell clusters and then applying antibody-derived tag (ADT) data to define molecularly identified BCR:ABL+ and BCR-ABL- primitive cells by expression of surface markers. This approach allowed them to show an association between the ratio of BCR-ABL+ vs BCR-ABL- primitive cells and TKI response and study dynamic changes in these populations following short-term TKI treatment.

Weaknesses:

The number of samples studied by CITE-Seq is limited. However, the authors have confirmed their key observations in additional samples. The BCR-ABL+ versus BCR-ABL- status of cells was not confirmed by direct sequencing for BCR-ABL. However, we recognize that the methodologies to perform these analyses on single cells is still evolving and the authors have shown that CD26 and CD35 expression can consistently identify BCR-ABL+ versus BCR-ABL- cells. It will be of interest to learn whether the GEP and surface markers identified here can distinguish BCR-ABL+ primitive stem cells later in the course of TKI treatment.

Reviewer #3 (Public Review):

Summary:

In this study, Warfvinge and colleagues use CITE-seq to interrogate how CML stem cells change between diagnosis and after one year of TKI therapy. This provides important insight into why some CML patients are "optimal responders" to TKI therapy while others experience treatment failure. CITE-seq in CML patients revealed several important findings. First, substantial cellular heterogeneity was observed at diagnosis, suggesting that this is a hallmark of CML. Further, patients who experienced treatment failure demonstrated increased numbers of primitive cells at diagnosis compared to optimal responders. This finding was validated in a bulk gene expression dataset from 59 CML patients, in which it was shown that the proportion of primitive cells versus lineage-primed cells correlates to treatment outcome. Even more importantly, because CITE-seq quantifies cell surface protein in addition to gene expression data, the authors were able to identify the BCR/ABL+ and BCR/ABL- CML stem cells express distinct cell surface markers (CD26+/CD35- and CD26-/CD35+, respectively). In optimal responders, BCR/ABL- CD26-/CD35+ CML stem cells were predominant, while the opposite was true in patients with treatment failure. Together, these findings represent a critical step forward for the CML field and may allow more informed development of CML therapies, as well as the ability to predict patient outcomes prior to treatment.

Strengths:

This is an important, beautifully written, well-referenced study that represents a fundamental advance in the CML field. The data are clean and compelling, demonstrating convincingly that optimal responders and patients with treatment failure display significant differences in the proportion of primitive cells at diagnosis, and the ratio of BCR-ABL+ versus negative LSCs. The finding that BCR/ABL+ versus negative LSCs display distinct surface markers is also key and will allow for more detailed interrogation of these cell populations at a molecular level.

Weaknesses:

CITE-seq was performed in only 9 CML patient samples and 2 healthy donors. Additional samples would greatly strengthen the very interesting and notable findings.

Reviewer #1 (Public Review):

This is an important manuscript that links gene expression to genetic variants and regions of open chromatin. The mechanisms of genetic gene regulation are essential to understanding how standing genetic variation translates to function and phenotype. This data set has the ability to add substantial insight into the field. In particular, the authors show how the relationships between variants, chromatin, and genes are spatially constrained by topologically associated domains.

Reviewer #2 (Public Review):

The experiments described in the manuscript are well designed and executed. Most of the data presented are of high quality, convincing, and in general support the conclusions made in the manuscript. This manuscript should be of great interest to the field of mammalian gene regulation and the approaches used here can have broader applications in studying genetic and epigenetic regulations of gene expression. The key finding reported here, the importance of 3D chromatin structure in controlling gene expression, although not unexpected, offers a better understanding of the physiological roles of TADs.

Comments on revised version:

I think the authors have substantially addressed reviewers' concerns. I have no further comments to add.

Reviewer #1 (Public Review):

Cell type deconvolution is one of the early and critical steps in the analysis and integration of spatial omic and single cell gene expression datasets, and there are already many approaches proposed for the analysis. Sang-aram et al. provide an up-to-date benchmark of computational methods for cell type deconvolution.

In doing so, they provide some (perhaps subtle) additional elements that I would say are above the average for a benchmarking study: i) a full Nextflow pipeline to reproduce their analyses; ii) methods implemented in Docker containers (which can be used by others to run their datasets); iii) a fairly decent assessment of their simulator compared to other spatial omics simulators. A key aspect of their results is that they are generally very concordant between real and synthetic datasets. And, it is important that they authors include an appropriate "simpler" baseline method to compare against and surprisingly, several methods performed below this baseline. Overall, this study also has the potential to also set the standard of benchmarks higher, because of these mentioned elements.

The only weakness of this study that I can readily see is that this is a very active area of research and we may see other types of data start to dominate (CosMx, Xenium) and new computational approaches will surely arrive. The Nextflow pipeline will make the the prospect of including new reference datasets and new computational methods easier.

Reviewer #2 (Public Review):

In this manuscript Sangaram et al provide a systematic methodology and pipeline for benchmarking cell type deconvolution algorithms for spatial transcriptomic data analysis in a reproducible manner. They developed a tissue pattern simulator that starts from single-cell RNA-seq data to create silver standards and used spatial aggregation strategies from real in situ-based spatial technologies to obtain gold standards. By using several established metrics combined with different deconvolution challenges they systematically scored and ranked 12 different methods and assessed both functional and usability criteria. Altogether, they present a reusable and extendable platform and reach very similar conclusions to other deconvolution benchmarking paper, including that RCTD, SpatialDWLS and Cell2location typically provide the best results. Major strengths of the simulation engine include the ability to downsample and recapitulate several cell and tissue organization patterns.

More specifically, the authors of this study sought to construct a methodology for benchmarking cell type deconvolution algorithms for spatial transcriptomic data analysis in a reproducible manner. The authors leveraged publicly available scRNA-seq, seqFISH, and STARMap datasets to create synthetic spatial datasets modeled after that of the Visium platform. It should be noted that the underlying experimental techniques of seqFISH and STARMap (in situ hybridization) do not parallel that of Visium (sequencing), which could potentially bias simulated data. Furthermore, to generate the ground truth datasets cells and their corresponding count matrix are represented by simple centroids. Although this simplifies the analysis it might not necessarily accurately reflect Visium spots where cells could lie on a boundary and affect deconvolution results.

The authors thoroughly and rigorously compare methods while addressing situational discrepancies in model performance, indicative of a strong analysis. The authors make a point to address both inter- and intra- dataset reference handling, which has a significant impact on performance, as the authors note in the text and conclusions. Indeed, supplying optimal reference data is - potentially most - important to achieve the best performance and hence it's important to understand that experimental design or sample matching is at least equally important to selecting the ideal deconvolution tool.

Similarly, the authors conclude that many methods are still outperformed by bulk deconvolution methods (e.g. Music or NNLS), however, it needs to be noted that these 'bulk' methods are also among the most sensitive when using an external (inter) dataset (S10), which likely resembles the more realistic scenario for most labs.

As the authors also discuss it's important to realize that deconvolution approaches are typically part of larger exploratory data analysis (EDA) efforts and require users to change parameters and input data multiple times. Thus, running time, computing needs, and scalability are probably key factors that researchers would like to consider when looking to deconvolve their datasets.

The authors achieve their aim to benchmark different deconvolution methods and the results from their thorough analysis support the conclusions that creating cell type deconvolution algorithms that can handle both cell abundance and rarity throughout a given tissue sample are challenging.

The reproducibility of the methods described will have significant utility for researchers looking to develop cell type deconvolution algorithms, as this platform will allow simultaneous replication of the described analysis and comparison to new methods.

Reviewer #1 (Public Review):

The authors primary objective in this study was to identify differences between patients with preeclampsia and normal patients with respect to the placental syncytiotrophoblast extracellular vesicle proteome.

A strength of this study is that the authors identified novel STB-EV protein markers that are more abundant in the placenta of patients with preeclampsia compared with normal controls. This contributes a little more to what is already known about STB-EV markers and preeclampsia. If these markers can be shown to be more abundant in maternal plasma of preeclampsia patients, it would be very useful for identifying patients who are at high risk for developing early-onset preeclampsia.

Weaknesses include:<br /> (1) The small sample size. There were only 6 patients in the study group and 6 normal controls. However, this can be considered as a pilot study.<br /> (2) The normal controls were not matched with the study patients and the authors did not state how the controls were selected.<br /> (3) The authors state that the placenta samples were obtained at the time of elective cesarean section. However, it is likely that all the preeclampsia patients were delivered for clinical indications rather than electively. This should be clarified.

Reviewer #2 (Public Review):

Summary:

Preeclampsia is a disorder of pregnancy that affects 4-5% of pregnancies worldwide. Identifying this condition early is clinically relevant as it will help clinicians to make management decisions to prevent adverse outcomes. The placenta holds a key to many pregnancy-related pathologies including preeclampsia and studies have shown many differences in the placenta of women with preeclampsia as compared to controls. However as the placenta cannot be collected directly during pregnancy, the exosomes secreted by it are considered a good alternative to tissue biopsy. In this study, the authors have compared the proteins in different sizes of exosomes from the placenta of women with and without preeclampsia. The idea is to eventually use these as biomarkers for early detection of preeclampsia.

Strengths:

The novelty factor of this study is the use of two different-sized exosomes which has not been achieved earlier.

Weaknesses:

The study measured the proteins at only a single time point after the disease has already occurred. However, the placenta is an ever-changing tissue throughout pregnancy and different proteins can come up at different times in pregnancy. Thus serial measurements are necessary and a single time point measurement. The has not validated the identified biomarkers in plasma or circulating placental exosomes from women with and without preeclampsia. Thus the utility of these findings in real-life situations can not be judged from this work.

Reviewer #1 (Public Review):

The authors set out to use structural biology (cryo-em), SPR and complement convertase assays to understand the mechanism(s) by which ISG65 dampens the cytotoxicity/cellular clearance to/of trypanosmes opsonised with C3b by the innate immune system.

The cryo-EM structure adds significantly the the author's previous crystallographic data because the latter was limited to the C3d sub-domain of C3b. Further, the in vitro convertase assay adds an additional functional dimension to this study.

The authors have achieved their aims and the results support their conclusions.

The role of complement in immunity to T. brucei (or lack thereof) has been a significant question in molecular parasitology for over 30 years. The identification of ISG65 as the C3 receptor and now this study providing mechanistic insights represents a major advance in the field.

The authors have appropriately put their results into perspective with other recent reports on the role of ISG65.

Reviewer #3 (Public Review):

The authors investigate the mechanisms by which ISG65 and C3 recognize and interact with each other. The major strength is the identification of exo-site by determining the cryoEM structure of the complex, which suggests new intervention strategies. This is a solid body of work that has an important impact in parasitology, immunology, and structural biology.

Comments on revised version:

The authors have addressed all the previous concerns.

Reviewer #1 (Public Review):

The goal of this study is to understand the allosteric mechanism of overall activity regulation in an anaerobic ribonucleotide reductase (RNR) that contains an ATP-cone domain. Through cryo-EM structural analysis of various nucleotide-bound states of the RNR, the mechanism of dATP inhibition is found to involve order-disorder transitions in the active site. These effects appear to prevent binding of substrate and a radical transfer needed to initiate the reaction.

Strengths of the manuscript include the comprehensive nature of the work - including both numerous structures of different forms of the RNR and detailed characterization of enzyme activity to establish the parameters of dATP inhibition. The manuscript has been improved in a revision by performing additional experiments to help corroborate certain aspects of the study. But these new experiments do not address all of the open questions about the structural basis for mechanism. Additionally, some questions about the strength of biochemical data and fit of binding or kinetic curves to data that were raised by other referees still remain. Some experimental observations are not consistent with the proposed model. For example, why does dATP enhance Gly radical formation when the proposed mechanism of dATP inhibition involves disorder in the Gly radical domain?

The work is impactful because it reports initial observations about a potentially new mode of allosteric inhibition in this enzyme class. It also sets the stage for future work to understand the molecular basis for this phenomenon in more detail.

Reviewer #3 (Public Review):

The manuscript by Bimai et al describes a structural and functional characterization of an anaerobic ribonucleotide reductase (RNR) enzyme from the human microbe, P. copri. More specifically, the authors aimed to characterize the mechanism by how (d)ATP modulates nucleotide reduction in this anaerobic RNR, using a combination of enzyme kinetics, binding thermodynamics, and cryo-EM structural determination, complemented by hydrogen-deuterium exchange (HDX). One of the principal findings of this paper is the ordering of a NxN 'flap' in the presence of ATP that promotes RNR catalysis and the disordering (or increased protein dynamics) of both this flap and the glycyl radical domain (GRD) when the inhibitory effector, dATP, binds. The latter is correlated with a loss of substrate binding, which is the likely mechanism for dATP inhibition. It is important to note that the GRD is remote (>30 Ang) from the binding site of the dATP molecule, suggesting long-range communication of the structural (dis)ordering. The authors also present evidence for a shift in oligomerization in the presence of dATP. The work does provide evidence for new insights/views into the subtle differences of nucleotide modulation (allostery) of RNR, in a class III system, through long-range interactions.

The strengths of the work are the impressive, in-depth structural analysis of the various regulated forms of PcRNR by (d)ATP using cryo-EM. The authors present seven different models in total, with striking differences in oligomerization and (dis)ordering of select structural features, including the GRD that is integral to catalysis. The authors present several, complementary biochemical experiments (ITC, MST, EPR, kinetics) aimed at resolving the binding and regulatory mechanism of the enzyme by various nucleotides. The authors present a good breadth of the literature in which the focus of allosteric regulation of RNRs has been on the aerobic orthologues.

The addition of hydrogen-deuterium exchange mass spectrometry (HDX-MS) complements the results originating from cryo-EM data. Most notably, is the observation of the enhanced exchange (albeit quite subtle) of the GRD domain in the presence of dATP that matches the loss of structural information in this region in the cryo-EM data. The most pronounced and compelling HDX results are seen in the form of dATP-induced protection of peptides immediately adjacent to the b-hairpin at the s-site, where dATP is expected to bind based on cryo-EM. It is clear that the presence of dATP increases the rigidity of this region.

Weaknesses: The discussion of the change in peptide mobility in the N-terminal region is complicated by the presence of bimodal mass spectral features and this may prevent detailed interpretation of the data, especially for select peptide region that shows opposite trends upon nucleotide association. Further, the HDX data in the NxN flap is unchanged upon nucleotide binding (ATP, dATP, or CTP), despite changes observed in the cryo-EM data.

Reviewer #3 (Public Review):

In the revised manuscript by Maio et al, the authors examined the bioenergetic mechanisms involved in the delayed migration of DC's during Mtb infection. The authors performed a series of in vitro infection experiments including bioenergetic experiments using the Agilent Seahorse XF, and glucose uptake and lactate production experiments. Also, data from SCENITH is included in the revised manuscript as well as some clinical data. This is a well written manuscript and addresses an important question in the TB field. A remaining weakness is the use of dead (irradiated) Mtb in several of the new experiments and claims where iMtb data were used to support live Mtb data. Another notable weakness lies in the author's insistence on asserting that lactate is the ultimate product of glycolysis, rather than acknowledging a large body of historical data in support of pyruvate's role in the process. This raises a perplexing issue highlighted by the authors: if Mtb indeed upregulates glycolysis, one would expect that inhibiting glycolysis would effectively control TB. However, the reality contradicts this expectation. Lastly, the examination of the bioenergetics of cells isolated from TB patients undergoing drug therapy, rather than studying them at their baseline state is a weakness.

Reviewer #1 (Public Review):

Summary:

The authors use insights into the dynamics of the PKA kinase domain, obtained by NMR experiments, to inform MD simulations that generate an energy landscape of PKA kinase domain conformational dynamics.

Strengths:

The authors integrate strong experimental data through the use of state-of-the-art MD studies and derive detailed insights into allosteric communication in PKA kinase. Comparison of wt kinase with a mutant (F100A) shows clear differences in the allosteric regulation of the two proteins. These differences can be rationalized by NMR and MD results. During the revision process, the authors have addressed the reviewers' comments adequately and have improved the accessibility of the manuscript to a wider audience.

Reviewer #3 (Public Review):

Summary:

Combining several MD simulation techniques (NMR-constrained replica-exchange metadynamics, Markov State Model, and unbiased MD) the authors identified the aC-beta4 loop of PKA kinase as a switch crucially involved in PKA nucleotide/substrate binding cooperatively. They identified a previously unreported excited conformational state of PKA (ES2), this switch controls and characterized ES2 energetics with respect to the ground state. Based on translating the simulations into chemical shits and NMR characterizing of PKA WT and an aC-beta4 mutant, the author made a convincing case in arguing that the simulation-suggested excited state is indeed an excited state observed by NMR, thus giving the excited state conformational details.

Strengths:

This work incorporates extensive simulation works, new NMR data, and in vitro biochemical analysis. It stands out in its comprehensiveness, and I think it made a great case.

Weaknesses:

The manuscript is somewhat difficult to read even for kinase experts, and even harder for the layman. The difficulty partially arises from mixing the technical description of the simulations with the structural interpretation of the results, which is more intuitive, and partially arises from the assumption that readers are familiar with kinase architecture and its key elements (the aC helix, the APE motif, etc).

Reviewer #1 (Public Review):

Summary:

The authors have implemented Optimal Transport algorithm in GromovMatcher for comparing LC/MS features from different datasets. This paper gains significance in the proteomics field for performing meta-analysis of LC/MS data.

Strengths:

The main strength is that GromovMatcher acheives significant performance metrics compared to other existing methods. The authors have done extensive comparisons to claim that GromovMatcher performs well.

Weaknesses:

The authors might need to add the limitation of datasets and thus have tested/validated their tool using simulated data in the abstract as well.

Reviewer #2 (Public Review):

Summary

The goal of untargeted metabolomics is to identify differences between metabolomes of different biological samples.

Untargeted metabolomics identifies features with specific mass-to-charge-ratio (m/z) and retention time (RT). Matching those to specific metabolites based on the model compounds from databases is laborious and not always possible, which is why methods for comparing samples on the level of unmatched features are crucial.<br /> The main purpose of the GromovMatcher method presented here is to merge and compare untargeted metabolomes from different experiments. These larger datasets could then be used to advance biological analyses, for example, for identification of metabolic disease markers.

The main problem that complicates merging different experiments is that m/z and RT vary slightly for the same feature (metabolite).

The main idea behind the GromovMatcher is built on the assumption that if two features match between two datasets (that feature i from dataset 1 matches feature j from dataset 2, and feature k from dataset 1 matches feature l from dataset 2), then the correlations or distances between the two features within each of the datasets (i and k, and j and l) will be similar. The authors then use the Gromov-Wasserstein method to find the best matches matrix from these data.

The variation in m/z between the same features in different experiments is a user-defined value and it is initially set to 0.01 ppm. There is no clear limit for RT deviations, so the method estimates a non-linear deviation (drift) of RT between two studies. GromovMatcher estimates the drift between two studies, and then discards the matching pairs where the drift would deviate significantly from the estimate. It learns the drift from a weighted spline regression.

The authors validate the performance of their GromovMatcher method using a dataset of cord blood. They use 20 different splits and compare the GromovMatcher (both its GM and GMT iterations, whereby GMT version uses the deviation from estimated RT drift to filter the matching matrix) with two other matching methods: M2S and metabCombiner.

The second validation was done using a (scaled and centered) dataset of metabolics from cancer datasets from the EPIC cohort that were manually matched by an expert. This dataset was also used to show that using automated methods can identify more features that are associated with a particular group of samples than what was found by manual matching. Specifically, the authors identify additional features connected to alcohol consumption.

Strengths:

I see the main strength of this work in its combination of all levels of information (m/z, RT, and higher-order information on correlations between features) and using each of the types of information in a way that is appropriate for the measure. The most innovative aspect is using the Gromov-Wasserstein method to match the features based on distance matrices.

The authors of the paper identify two main shortcomings with previously established methods that attempt to match features from different experiments: a) all other methods require fine-tuning of user-defined parameters, and, more importantly, b) do not consider correlations between features. The main strength of the GromovMatcher is that it incorporates the information on distances between the features (in addition to also using m/z and RT).

Weaknesses:

The main weakness is that there seem not to be enough manually curated datasets that could be used for validation. It will, therefore, be important, for the authors, and the field in general to keep validating and improving their methods if more datasets become available.

The second weakness, as emphasized by the authors in the discussion is that the method as it is set up now can be directly used only to compare two datasets. I am confident that the authors will successfully implement novel algorithms to address this issue in the future.

Reviewer #1 (Public Review):

Theoretical principles of viscous fluid mechanics are used here to assess likely mechanisms of transport in the ER. A set of candidate mechanisms is evaluated, making good use of imaging to represent ER network geometries. Evidence is provided that the contraction of peripheral sheets provides a much more credible mechanism than the contraction of individual tubules, junctions, or perinuclear sheets.

The work has been conducted carefully and comprehensively, making good use of underlying physical principles. There is a good discussion of the role of slip; sensible approximations (low volume fraction, small particle size, slender geometries, pragmatic treatment of boundary conditions) allow tractable and transparent calculations; clear physical arguments provide useful bounds; stochastic and deterministic features of the problem are well integrated.

There are just a couple of areas where more discussion might be warranted, in my view.

(1) The energetic cost of tubule contraction is estimated, but I did not see an equivalent estimate for the contraction of peripheral sheets. It might be helpful to estimate the energetic cost of viscous dissipation in generated flows at higher frequencies. The mechanism of peripheral sheet contraction is unclear: do ATP-driven mechanisms somehow interact with thermal fluctuations of membranes?

(2) Mutations are mentioned in the abstract but not (as far as I could see) later in the manuscript. It would be helpful if any consequences for pathologies could be developed in the text.

Reviewer #2 (Public Review):

Summary:

This study explores theoretically the consequences of structural fluctuations of the endoplasmic reticulum (ER) morphology called contractions on molecular transport. Most of the manuscript consists of the construction of an interesting theoretical flow field (physical model) under various hypothetical assumptions. The computational modeling is followed by some simulations

Strengths:

The authors are focusing their attention on testing the hypothesis that a local flow in the tubule could be driven by tubular pinching. We recall that trafficking in the ER is considered to be mostly driven by diffusion at least at a spatial scale that is large enough to account for averaging of any random flow occurring from multiple directions [note that this is not the case for plants].

Weaknesses:

The manuscript extensively details the construction of the theoretical model, occupying a significant portion of the manuscript. While this section contains interesting computations, its relevance and utility could be better emphasized, perhaps warranting a reorganization of the manuscript to foreground this critical aspect.

Overall, the manuscript appears highly technical with limited conclusive insights, particularly lacking predictions confirmed by experimental validation. There is an absence of substantial conclusions regarding molecular trafficking within the ER.

Reviewer #1 (Public Review):

Summary:

In this paper, proteomics analysis of the plasma of human subjects that underwent an exercise training regime consisting of a combination of endurance and resistance exercise led to the identification of several proteins that were responsive to exercise training. Confirming previous studies, many exercise-responsive secreted proteins were found to be involved in the extra-cellular matrix. The protein CD300LG was singled out as a potential novel exercise biomarker and the subject of numerous follow-up analyses. The levels of CD300LG were correlated with insulin sensitivity. The analysis of various open-source datasets led to the tentative suggestion that CD300LG might be connected with angiogenesis, liver fat, and insulin sensitivity. CD300LG was found to be most highly expressed in subcutaneous adipose tissue and specifically in venular endothelial cells. In a subset of subjects from the UK Biobank, serum CD300LG levels were positively associated with several measures of physical activity - particularly vigorous activity. In addition, serum CD300LG levels were negatively associated with glucose levels and type 2 diabetes. Genetic studies hinted at these associations possibly being causal. Mice carrying alterations in the CD300LG gene displayed impaired glucose tolerance, but no change in fasting glucose and insulin. Whether the production of CD300LG is changed in the mutant mice is unclear.

Strengths:

The specific proteomics approach conducted to identify novel proteins impacted by exercise training is new. The authors are resourceful in the exploitation of existing datasets to gain additional information on CD300LG.

Weaknesses:

While the analyses of multiple open-source datasets are necessary and useful, they lead to relatively unspecific correlative data that collectively insufficiently advance our knowledge of CD300LG and merely represent the starting point for more detailed investigations. Additional more targeted experiments of CD300LG are necessary to gain a better understanding of the role of CD300LG and the mechanism by which exercise training may influence CD300LG levels. One should also be careful to rely on external data for such delicate experiments as mouse phenotyping. Can the authors vouch for the quality of the data collected?

Reviewer #2 (Public Review):

Summary:

This manuscript from Lee-Odegard et al reports proteomic profiling of exercise plasma in humans, leading to the discovery of CD300LG as a secreted exercise-inducible plasma protein. Correlational studies show associations of CD300LG with glycemic traits. Lastly, the authors query available public data from CD300LG-KO mice to establish a causal role for CD300LG as a potential link between exercise and glucose metabolism. However, the strengths of this manuscript were balanced by the moderate to major weaknesses. Therefore in my opinion, while this is an interesting study, the conclusions remain preliminary and are not fully supported by the experiments shown so far.

Strengths:

(1) Data from a well-phenotyped human cohort showing exercise-inducible increases in CD300LG.

(2) Associations between CD300LG and glucose and other cardiometabolic traits in humans, that have not previously been reported.

(3) Correlation to CD300LG mRNA levels in adipose provides additional evidence for exercise-inducible increases in CD300LG.

Weaknesses:

(1) CD300LG is by sequence a single-pass transmembrane protein that is exclusively localized to the plasma membrane. How CD300LG can be secreted remains a mystery. More evidence should be provided to understand the molecular nature of circulating CD300LG. Is it full-length? Is there a cleaved fragment? Where is the epitope where the o-link is binding to CD300LG? Does transfection of CD300LG to cells in vitro result in secreted CD300LG?

(2) There is a growing recognition of specificity issues with both the O-link and somalogic platforms. Therefore it is critical that the authors use antibodies, targeted mass spectrometry, or some other methods to validate that CD300LG really is increased instead of just relying on the O-link data.

(3) It is insufficient simply to query the IMPC phenotyping data for CD300LG; the authors should obtain the animals and reproduce or determine the glucose phenotypes in their own hands. In addition, this would allow the investigators to answer key questions like the phenotype of these animals after a GTT, whether glucose production or glucose uptake is affected, whether insulin secretion in response to glucose is normal, effects of high-fat diet, and other standard mouse metabolic phenotyping assays.

(4) I was unable to find the time point at which plasma was collected at the 12-week time point. Was it immediately after the last bout of exercise (an acute response) or after some time after the training protocol (trained state)?

Reviewer #3 (Public Review):

Summary:

This manuscript by Liu et al. presents a case that CAPSL mutations are a cause of familial exudative vitreoretinopathy (FEVR). Attention was initially focused on the CAPSL gene from whole exome sequence analysis of two small families. The follow-up analyses included studies in which CAPSL was manipulated in endothelial cells of mice and multiple iterations of molecular and cellular analyses. Together, the data show that CAPSL influences endothelial cell proliferation and migration. Molecularly, transcriptomic and proteomic analyses suggest that CAPSL influences many genes/proteins that are also downstream targets of MYC and may be important to the mechanisms.

Strengths:

This multi-pronged approach found a previously unknown function for CAPSLs in endothelial cells and pointed at MYC pathways as high-quality candidates in the mechanism.

Weaknesses:

Two issues shape the overall impact for me. First, the unreported population frequency of the variants in the manuscript makes it unclear if CAPSL should be considered an interesting candidate possibly contributing to FEVR, or possibly a cause. Second, it is unclear if the identified variants act dominantly, as indicated in the pedigrees. The studies in mice utilized homozygotes for an endothelial cell-specific knockout, leaving uncertainty about what phenotypes might be observed if mice heterozygous for a ubiquitous knockout had instead been studied.

In my opinion, the following scientific issues are specific weaknesses that should be addressed:

(1) Please state in the manuscript the number of FEVR families that were studied by WES. Please also describe if the families had been selected for the absence of known mutations, and/or what percentage lack known pathogenic variants.

(2) A better clinical description of family 3104 would enhance the manuscript, especially the father. It is unclear what "manifested with FEVR symptoms, according to the medical records" means. Was the father diagnosed with FEVR? If the father has some iteration of a mild case, please describe it in more detail. If the lack of clinical images in the figure is indicative of a lack of medical documentation, please note this in the manuscript.

(3) The TGA stop codon can in some instances also influence splicing (PMID: 38012313). Please add a bioinformatic assessment of splicing prediction to the assays and report its output in the manuscript.

(4) More details regarding utilizing a "loxp-flanked allele of CAPSL" are needed. Is this an existing allele, if so, what is the allele and citation? If new (as suggested by S1), the newly generated CAPSL mutant mouse strain needs to be entered into the MGI database and assigned an official allele name - which should then be utilized in the manuscript and who generated the strain (presumably a core or company?) must be described.

(5) The statement in the methods "All mice used in the study were on a C57BL/6J genetic background," should be better defined. Was the new allele generated on a pure C57BL/6J genetic background, or bred to be some level of congenic? If congenic, to what generation? If unknown, please either test and report the homogeneity of the background, or consult with nomenclature experts (such as available through MGI) to adopt the appropriate F?+NX type designation. This also pertains to the Pdgfb-iCreER mice, which reference 43 describes as having been generated in an F2 population of C57BL/6 X CBA and did not designate the sub-strain of C57BL/6 mice. It is important because one of the explanations for missing heritability in FEVR may be a high level of dependence on genetic background. From the information in the current description, it is also not inherently obvious that the mice studied did not harbor confounding mutations such as rd1 or rd8.

(6) In my opinion, more experimental detail is needed regarding Figures 2 and 3. How many fields, of how many retinas and mice were analyzed in Figure 2? How many mice were assessed in Figure 3?

(7) I suggest adding into the methods whether P-values were corrected for multiple tests.

Reviewer #2 (Public Review):

Summary:

The authors have developed marker selection and k-means (k=2) based binary clustering algorithm for the first-level supervised clustering of the CyTOF dataset. They built a seamless pipeline that offers the multiple functionalities required for CyTOF data analysis.

Strengths:

The strength of the study is the potential use of the pipeline for the CyTOF community as a wrapper for multiple functions required for the analysis. The concept of the first line of binary clustering with known markers can be practically powerful.

Weaknesses:

The weakness of the study is that there's little conceptual novelty in the algorithms suggested from the study and the benchmarking is done in limited conditions.

Reviewer #1 (Public Review):

Summary:

This manuscript presented a useful toolkit designed for CyTOF data analysis, which integrates 5 key steps as an analytical framework. A semi-supervised clustering tool was developed, and its performance was tested in multiple independent datasets. The tool was compared to human experts as well as supervised and unsupervised methods.

Strengths:

The study employed multiple independent datasets to test the pipeline. A new semi-supervised clustering method was developed.

Weaknesses:

The examination of the whole pipeline is incomplete. Lack of descriptions or justifications for some analyses.

Reviewer #3 (Public Review):

Summary:

ImmCellTyper is a new toolkit for Cytometry by time-of-flight data analysis. It includes BinaryClust, a semi-supervised clustering tool (which takes into account prior biological knowledge), designed for automated classification and annotation of specific cell types and subpopulations. ImmCellTyper also integrates a variety of tools to perform data quality analysis, batch effect correction, dimension reduction, unsupervised clustering, and differential analysis.

Strengths:

The proposed algorithm takes into account the prior knowledge.<br /> The results on different benchmarks indicate competitive or better performance (in terms of accuracy and speed) depending on the method.

Weaknesses:

The proposed algorithm considers only CyTOF markers with binary distribution.

Reviewer #1 (Public Review):

Summary:

This study identifies new types of interactions between Drosophila gustatory receptor neurons (GRNs) and shows that these interactions influence sensory responses and behavior. The authors find that HCN, a hyperpolarization-activated cation channel, suppresses the activity of GRNs in which it is expressed, preventing those GRNs from depleting the sensillum potential, and thereby promoting the activity of neighboring GRNs in the same sensilla. HCN is expressed in sugar GRNs, so HCN dampens the excitation of sugar GRNs and promotes the excitation of bitter GRNs. Impairing HCN expression in sugar GRNs depletes the sensillum potential and decreases bitter responses, especially when flies are fed on a sugar-rich diet, and this leads to decreased bitter aversion in a feeding assay. The authors' conclusions are supported by genetic manipulations, electrophysiological recordings, and behavioral assays.

Strengths:

(1) Non-synaptic interactions between neurons that share an extracellular environment (sometimes called "ephaptic" interactions) have not been well-studied, and certainly not in the insect taste system. A major strength of this study is the new insight it provides into how these interactions can impact sensory coding and behavior.

(2) The authors use many different types of genetic manipulations to dissect the role of HCN in GRN function, including mutants, RNAi, overexpression, ectopic expression, and neuronal silencing. Their results convincingly show that HCN impacts the sensillum potential and has both cell-autonomous and nonautonomous effects that go in opposite directions. There are a couple of conflicting or counterintuitive results, but the authors discuss potential explanations.

(3) Experiments comparing flies raised on different food sources suggest an explanation for why the system may have evolved the way that it did: when flies live in a sugar-rich environment, their bitter sensitivity decreases, and HCN expression in sugar GRNs helps to counteract this decrease.

Weaknesses/Limitations:

(1) The genetic manipulations were constitutive (e.g. Ih mutations, RNAi, or misexpression), and depleting Ih from birth could lead to compensatory effects that change the function of the neurons or sensillum. Using tools to temporally control Ih expression could help to confirm the results of this study.

(2) The behavioral experiment shows a striking loss of bitter sensitivity, but it was only conducted for one bitter compound at one concentration. It is not clear how general this effect is. The same is true for some of the bitter GRN electrophysiological experiments that only tested one compound and concentration.

(3) Several experiments using the Gal4/UAS system only show the Gal4/+ control and not the UAS/+ control (or occasionally neither control). Since some of the measurements in control flies seem to vary (e.g., spiking rate), it is important to compare the experimental flies to both controls to ensure that any observed effects are in fact due to the transgene expression.

(4) I was surprised that manipulations of sugar GRNs (e.g. Ih knockdown, Gr64a-f deletion, or Kir silencing) can impact the sensillum potential and bitter GRN responses even in experiments where no sugar was presented. I believe the authors are suggesting that the effects of sugar GRN activity (e.g., from consuming sugar in the fly food prior to the experiment) can have long-lasting effects, but it wasn't entirely clear if this is their primary explanation or on what timescale those long-lasting effects would occur. How much / how long of a sugar exposure do the flies need for these effects to be triggered, and how long do those effects last once sugar is removed?

(5) The authors mention that HCN may impact the resting potential in addition to changing the excitability of the cell through various mechanisms. It would be informative to record the resting potential and other neuronal properties, but this is very difficult for GRNs, so the current study is not able to determine exactly how HCN affects GRN activity.

Reviewer #2 (Public Review):

Summary:

In this manuscript, the authors start by showing that HCN loss-of-function mutation causes a decrease in spiking in bitter GRNs (bGRN) while leaving sweet GRN (sGRN) response in the same sensillum intact. They show that a perturbation of HCN channels in sweet-sensing neurons causes a similar decrease while increasing the response of sugar neurons. They were also able to rescue the response by exogenous expression. Ectopic expression of HCN in bitter neurons had no effect. Next, they measure the sensillum potential and find that sensillum potential is also affected by HCN channel perturbation. These findings lead them to speculate that HCN in sGRN increases sGRN spiking which in turn affects bGRNs. To test this idea that carried out multiple perturbations aimed at decreasing sGRN activity. They found that decreasing sGRN activity by either using receptor mutant or by expressing Kir (a K+ channel) in sGRN increased bGRN responses. These responses also increase the sensillum potential. Finally, they show that these changes are behaviorally relevant as conditions that increase sGRN activity decrease avoidance of bitter substances.

Strengths:

There is solid evidence that perturbation of sweet GRNs affects bitter GRN in the same sensillum. The measurement of transsynaptic potential and how it changes is also interesting and supports the authors' conclusion.

Weaknesses:<br /> The ionic basis of how perturbation in GRN affects the transepithelial potential which in turn affects the second neuron is not clear.

Reviewer #3 (Public Review):

Ephaptic inhibition between neurons housed in the same sensilla has been long discovered in flies, but the molecular basis underlying this inhibition is underexplored. Specifically, it remains poorly understood which receptors or channels are important for maintaining the transepithelial potential between the sensillum lymph and the hemolymph (known as the sensillum potential), and how this affects the excitability of neurons housed in the same sensilla.

Lee et al. used single-sensillum recordings (SSR) of the labellar taste sensilla to demonstrate that the HCN channel, Ih, is critical for maintaining sensillum potential in flies. Ih is expressed in sugar-sensing GRNs (sGRNs) but affects the excitability of both the sGRNs and the bitter-sensing GRNs (bGRNs) in the same sensilla. Ih mutant flies have decreased sensillum potential, and bGRNs of Ih mutant flies have a decreased response to the bitter compound caffeine. Interestingly, ectopic expression of Ih in bGRNs also increases sGRN response to sucrose, suggesting that Ih-dependent increase in sensillum potential is not specific to Ih expressed in sGRNs. The authors further demonstrated, using both SSR and behavior assays, that exposure to sugars in the food substrate is important for the Ih-dependent sensitization of bGRNs. The experiments conducted in this paper are of interest to the chemosensory field. The observation that Ih is important for the activity in bGRNs albeit expressed in sGRNs is especially fascinating and highlights the importance of non-synaptic interactions in the taste system.

Despite the interesting results, this paper is not written in a clear and easily understandable manner. It uses poorly defined terms without much elaboration, contains sentences that are borderline unreadable even for those in the narrower chemosensory field, and many figures can clearly benefit from more labeling and explanation. It certainly needs a bit of work.

Below are the major points:

(1) Throughout the paper, it is assumed that Ih channels are expressed in sugar-sensing GRNs but not bitter-sensing GRNs. However, both this paper and citation #17, another paper from the same lab, contain only circumstantial evidence for the expression of Ih channels in sGRNs. A simple co-expression analysis, using the Ih-T2A-GAL4 line and Gr5a-LexA/Gr66a-LexA line, all of which are available, could easily demonstrate the co-expression. Including such a figure would significantly strengthen the conclusion of this paper.

(2) Throughout this paper, it is often unclear which class of labellar taste sensilla is being recorded. S-a, S-b, I-a, and I-b sensilla all have different sensitivities to bitters and sugars. Each figure should clearly indicate which sensilla is being recorded. Justification should be provided if recordings from different classes of sensilla are being pooled together for statistics.

(3) In many figures, there is a lack of critical control experiments. Examples include Figures 1C-F (lacking UAS control), Figure 2I-J (lacking UAS control), Figure 4E (lacking the UAS and GAL4 control, and it is also strange to compare Gr64f > RNAi with Gr66a > RNAi, instead of with parental GAL4 and UAS controls.), and Figure 5D (lacking UAS control). Without these critical control experiments, it is difficult to evaluate the quality of the work.

(4) Figure 2A could benefit from more clarification about what exactly is being recorded here. The text is confusing: a considerable amount of text is spent on explaining the technical details of how SP is recorded, but very little text about what SP represents, which is critical for the readers. The authors should clarify in the text that SP is measuring the potential between the sensillar lymph, where the dendrites of GRNs are immersed, and the hemolymph. Adding a schematic figure to show that SP represents the potential between the sensillar lymph and hemolymph would be beneficial.

(5) The sGRN spiking rate in Figure 4B deviates significantly from previous literature (Wang, Carlson, eLife 2022; Jiao, Montell PNAS 2007, as examples), and the response to sucrose in the control flies is not dosage-dependent, which raises questions about the quality of the data. Why are the responses to sucrose not dosage-dependent? The responses are clearly not saturated at these (10 mM to 100 mM) concentrations.

(6) In Figure 4C, instead of showing the average spike rate of the first five seconds and the next 5 seconds, why not show a peristimulus time histogram? It would help the readers tremendously, and it would also show how quickly the spike rate adapts to overexpression and control flies. Also, since taste responses adapt rather quickly, a 500 ms or 1 s bin would be more appropriate than a 5-second bin.

(7) Lines 215 - 220. The authors state that the presence of sugars in the culture media would expose the GRNs to sugar constantly, without providing much evidence. What is the evidence that the GRNs are being activated constantly in flies raised with culture media containing sugars? The sensilla are not always in contact with the food.

(8) Line 223. To show that bGRN spike rates in Ih mutant flies "decreased even more than WT", you need to compare the difference in spike rates between the sorbitol group and the sorbitol + sucrose group, which is not what is currently shown.

(9) To help readers better understand the proposed mechanisms here, including a schematic figure would be helpful. This should show where Ih is expressed, how Ih in sGRNs impacts the sensillum potential, how elevated sensillum potential increases the electrical driving force for the receptor current, and affects the excitability of the bGRNs in the same sensilla, and how exposure to sugar is proposed to affect ion homeostasis in the sensillum lymph.

Reviewer #2 (Public Review):

The authors used a whole genome CRISPR screen to identify targetable synthetic lethalities associated with PPM1D mutations, known poor prognosis and currently undruggable factors in leukemia. The authors identified the cytosolic superoxide dismutase (SOD1, Cu/Zn SOD) as a major protective factor in PPMD1 mutant vs. wt cells, and their study investigates associated mechanisms of this protection. Using both genetic depletion and small molecule inhibitors of SOD1, the authors conclude that SOD1 loss exacerbates mitochondrial dysfunction, ROS levels and DNA damage phenotypes in PPM1D mutant cells, decreasing cell growth in AML cells. The data strongly support that PPMD1 mutant cells have high levels of total peroxides and elevated DNA breaks, and that genetic depletion of SOD1 decreases cell growth in two AML cell lines. However, the authors don't explain how superoxide radical (which is not damaging by itself) induces such damage, the on-target effects of the SOD1 inhibitors at the concentrations is not clear, the increase in total hydroperoxides is not supported by loss of SOD1, the changes in mitochondrial function are small, and there is no assessment of how the mitochondrial SOD2 expression or function, which dismutates mitochondrial superoxide, is altered. Overall these studies do not distinguish between signal vs. damaging aspects of ROS in their models and do not rule out an alternate hypothesis that loss of SOD1 increases superoxide production by cytosolic NADPH activity which would significantly alter ROS-driven regulation of kinase/phosphatase signal modulation, affecting cell growth and proliferation as well as DNA repair. Additionally, with the exception of growth defects demonstrated with sgSOD1, the majority of data are acquired using two chemical inhibitors, LCS1 and ATN-224, without supporting evidence that these inhibitors are acting in an on-target manner.

Overall, the authors address an important problem by seeking targetable vulnerabilities in PPM1D mutant AML cells, it is clear SOD1 deletion induces strong growth defects in the AML cell lines tested, most of the approaches are appropriate for the outcomes being evaluated, and the data are technically solid and well-presented. The major weakness lies in which redox pathways and ROS species are evaluated, how the resulting data are interpreted, and gaps in the follow-up experiments. Due to these omissions, as currently presented, the broader impact of these findings are unclear.

These specific concerns are outlined in detail below and I offer some suggestions regarding how to clarify the mechanisms underlying their initial observation of SOD1 synthetic lethality:

(1) Fig. 1 - SOD1 appears to be clustered with several other genes in the volcano plot (including FANC proteins). Did any other ROS-detoxifying enzymes show similar fitness scores? The effects of the SOD1 sgRNA are striking, however it would be useful to see qPCR or immunoblot data confirming robust depletion.

Does SOD1 co-expression in PPM1-mutant patient AML correspond to poorer disease outcomes? This can be evaluated in publicly available patient datasets and would support the idea of SOD1 synthetic lethality.

It would also be useful to know (given the subsequent results) whether expression of the SOD2, the mitochondrial superoxide dismutase, is altered in response to SOD1 loss.

(2) Fig. 2 - What are the relative SOD1 levels in the mutant PPM1D vs. wt. cell lines? The effects of the chemical inhibitors are stronger in MOLM-13 than the other two lines. These data could also point to whether LCS-1 and ATN-224 cytotoxicity is on-target or off-target at these concentrations, which is a key issue not currently addressed in these studies. This is a particular concern as the OCI-AML2 line shows a stronger growth defect with CRISPR SOD1 KO (in Fig 1) but the smallest effects with these chemical inhibitors.

While endogenous mitochondrial superoxide levels are elevated in PPM1D mutant lines, it is entirely unclear why SOD1 inhibition should affect mitochondrial superoxide as it detoxifies cytosolic superoxide. Also unclear why DCFDA signal (which measures total hydroperoxides) is *increased* under SOD1 inhibition - SOD1 dismutates superoxide radicals into hydrogen peroxide, therefore unless SOD2 is compensating for SOD1 loss, one might expect hydroperoxides to be lower (unless some entirely different oxidase is increasing their levels). None of these outcomes appear to be considered. Finally, it is not explained how lipid peroxidation, which requires production of hydroxyl or similarly high potency radicals, is being caused by increased superoxide or peroxides. One possibility is there is an increase in labile iron, in which case this phenotype would be rescued by the iron chelator desferal, and by the lipophilic antioxidant, ferrostatin.

Do the sgSOD1 cells also show similar increases in MitoSox green, DCFDA and BODIPY signal? These experiments would clarify whether the effects with the inhibitors are directly related directly to SOD1 loss or if they represent off-target effects from the inhibitors and/or compensatory changes in SOD2.

(3) Fig. 3 - the effects on mitochondrial respiratory parameters, while statistically significant, do not seem biologically striking. Also, these data are shown for OCI-AML2 cells which show the smallest cytotoxic effects with the SOD1 inhibitors among the 3 lines tested. They do however show the most robust growth defect with sgSOD1. This discrepancy could suggest that mitochondrial dysfunction does not underlie the observed growth defect and/or the inhibitor cytotoxicity is not on-target. Ideally mitochondrial profiling should also be carried out on this cell line with inducible SOD1 depletion. Have the authors assessed whether the mitochondrial Bcl family proteins are affected by the inhibitors?

(4) Fig. 4 - Currently the data in this figure do not support the authors claim that PPM1D-mutant cells have impaired antioxidant defense mechanisms, leading to an elevation in ROS levels and reliance on SOD1 for protection. It should be noted that oxidative stress specifically refers to adverse cellular effects of increasing ROS, not baseline levels of various redox parameters. Ideally levels of GSSG/GSH would be a better measure of potential redox stress tolerance than the total antioxidant capacity assay. Finally, oxidative stress can be assessed by challenging the wt and mutant PPM1D cell lines with oxidant stressors such as paraquat which elevates superoxide or drugs like erastin which elevate mitochondrial ROS. The immunoblot shows negligible changes in the antioxidant proteins assayed. Again, this blot should include SOD2 which is the most relevant antioxidant in the context of mitochondrial superoxide.

(5) Fig. 5 - These data support that DNA breaks are elevated in PPM1D mutant vs. wt cells. However, the data with the chemical SOD1 inhibitor again do not convince that the enhanced levels are due to on-target effects on SOD1. Use of the alkaline comet assay is appropriate for these studies and the 8-oxoguanine data do indicate contributions from oxidative DNA base damage. But these are unlikely to result directly from altered superoxide levels, as this species cannot directly oxidize DNA bases or cause DNA strand breaks.

The following points summarize my specific experimental and textual recommendations:

(1) These studies require an assessment of on-target efficacy of the inhibitors at the relevant concentration ranges. Ideally, they should have minimal effects against SOD1 knockout cell lines (acute challenge at a time point before the growth defects become apparent) and show better efficacy in SOD1-overexpressing lines. Key experiments (changes in superoxide, OCR profiling, DNA alkaline comet assay) would be more convincing if they are carried out with SOD1 knockout lines to compare against the inhibitor effects (3-4 days after introducing sgSOD1 when growth defects are not apparent).

(2) Instead of using NAC, which elevates glutathione synthesis but also has several known side-effects, the authors may want to determine whether Tempol, a SOD mimetic can rescue the effects of SOD1 knockout or inhibition. This would directly prove that SOD1 functional loss underlies the observed growth defect and cytotoxicity from genetic SOD1 knockdown or chemical inhibition.

(3) The complete lack of consideration of SOD2 in these studies is a missed opportunity as it reduces mitochondrial superoxide levels but elevates hydrogen peroxide levels. It would be very interesting to see whether SOD1 inhibition leads to compensatory increases in SOD2. SOD2 can be easily measured by immunoblot. Furthermore, measuring total superoxide via hydroethidium in a flow cytometric assay vs. mitochondrial ROS in PPM1D mut vs. wt cells and under SOD1 knockout would enable a determination of which species dominates (cytosolic or mitochondrial). These experiments are required to fill some logical gaps in interpretation of their redox data.

(4) Given the DNA breaks observed in PPM1D mutant cells, it is highly recommended the authors assess whether iron levels are elevated in mut vs. wt cells and whether desferal can rescue observed SOD1 inhibition defects.

(5) The authors may want to assess whether Rac1 or NADPH oxidase activity is altered in the SOD1 KO in wt vs. PPM1D cells. Their results may be the consequence of compromised ROS-driven survival signaling or DNA repair rather than direct ROS-induced damage, which is not caused directly by superoxide (or hydrogen peroxide).

(6) It is recommended the discussion focus more strongly on how the signaling function of superoxide vs. its reactions with other molecular entities to induce genotoxic outcomes could be contributing to the observed phenotypes. The discussion of FANC proteins, which were targets with similar fitness scores but not experimentally investigated at all, is an unwarranted digression.

Reviewer #1 (Public Review):

Summary:

The authors demonstrate that it is possible to carry out eQTL experiments for the model eukaryote S. cerevisiae, in "one pot" preparations, by using single-cell sequencing technologies to simultaneously genotype and measure expression. This is a very appealing approach for investigators studying genetic variation in single-celled and other microbial systems, and will likely inspire similar approaches in non-microbial systems where comparable cell mixtures of genetically heterogeneous individuals could be achieved.

Strengths:

While eQTL experiments have been done for nearly two decades (the corresponding author's lab are pioneers in this field), this single-cell approach creates the possibility for new insights about cell biology that would be extremely challenging to infer using bulk sequencing approaches. The major motivating application shown here is to discover cell occupancy QTL, i.e. loci where genetic variation contributes to differences in the relative occupancy of different cell cycle stages. The authors dissect and validate one such cell cycle occupancy QTL, involving the gene GPA1, a G-protein subunit that plays a role in regulating the mating response MAPK pathway. They show that variation at GPA1 is associated with proportional differences in the fraction of cells in the G1 stage of the cell cycle. Furthermore, they show that this bias is associated with differences in mating efficiency.

Weaknesses:

While the experimental validation of the role of GPA1 variation is well done, the novel cell cycle occupancy QTL aspect of the study is somewhat underexploited. The cell occupancy QTLs that are mentioned all involve loci that the authors have identified in prior studies that involved the same yeast crosses used here. It would be interesting to know what new insights, besides the "usual suspects", the analysis reveals. For example, in Cross B there is another large effect cell occupancy QTL on Chr XI that affects the G1/S stage. What candidate genes and alleles are at this locus? And since cell cycle stages are not biologically independent (a delay in G1, could have a knock-on effect on the frequency of cells with that genotype in G1/S), it would seem important to consider the set of QTLs in concert.

Reviewer #2 (Public Review):

Boocock and colleagues present an approach whereby eQTL analysis can be carried out by scRNA-Seq alone, in a one-pot-shot experiment, due to genotypes being able to be inferred from SNPs identified in RNA-Seq reads. This approach obviates the need to isolate individual spores, genotype them separately by low-coverage sequencing, and then perform RNA-Seq on each spore separately. This is a substantial advance and opens up the possibility to straightforwardly identify eQTLs over many conditions in a cost-efficient manner. Overall, I found the paper to be well-written and well-motivated, and have no issues with either the methodological/analytical approach (though eQTL analysis is not my expertise), or with the manuscript's conclusions.

I do have several questions/comments.

393 segregant experiment:<br /> For the experiment with the 393 previously genotyped segregants, did the authors examine whether averaging the expression by genotype for single cells gave expression profiles similar to the bulk RNA-Seq data generated from those genotypes? Also, is it possible (and maybe not, due to the asynchronous nature of the cell culture) to use the expression data to aid in genotyping for those cells whose genotypes are ambiguous? I presume it might be if one has a sufficient number of cells for each genotype, though, for the subsequent one-pot experiments, this is a moot point.

Figure 1B:<br /> Is UMAP necessary to observe an ellipse/circle - I wouldn't be surprised if a simple PCA would have sufficed, and given the current discussion about whether UMAP is ever appropriate for interpreting scRNA-Seq (or ancestry) data, it seems the PCA would be a preferable approach. I would expect that the periodic elements are contained in 2 of the first 3 principal components. Also, it would be nice if there were a supplementary figure similar to Figure 4 of Macosko et al (PMID 26000488) to indeed show the cell cycle dependent expression.

Aging, growth rate, and bet-hedging:<br /> The mention of bet-hedging reminded me of Levy et al (PMID 22589700), where they saw that Tsl1 expression changed as cells aged and that this impacted a cell's ability to survive heat stress. This bet-hedging strategy meant that the older, slower-growing cells were more likely to survive, so I wondered a couple of things. It is possible from single-cell data to identify either an aging, or a growth rate signature? A number of papers from David Botstein's group culminated in a paper that showed that they could use a gene expression signature to predict instantaneous growth rate (PMID 19119411) and I wondered if a) this is possible from single-cell data, and b) whether in the slower growing cells, they see markers of aging, whether these two signatures might impact the ability to detect eQTLs, and if they are detected, whether they could in some way be accounted for to improve detection.

AIL vs. F2 segregants:<br /> I'm curious if the authors have given thought to the trade-offs of developing advanced intercross lines for scRNA-Seq eQTL analysis. My impression is that AIL provides better mapping resolution, but at the expense of having to generate the lines. It might be useful to see some discussion on that.

10x vs SPLit-Seq<br /> 10x is a well established, but fairly expensive approach for scRNA-Seq - I wondered how the cost of the 10x approach compares to the previously used approach of genotyping segregants and performing bulk RNA-Seq, and how those costs would change if one used SPLiT-Seq (see PMID 38282330).

Reviewer #1 (Public Review):

The manuscript investigates the role of the membrane-deforming cytoskeletal regulator protein Abba in cortical development and its potential implications for microcephaly. It is a valuable contribution to the understanding of Abba's role in cortical development. The strengths and weaknesses identified in the manuscript are outlined below:

Clinical Relevance:

The authors identified a patient with microcephaly and a patient with an intellectual disability harboring a mutation in the Abba variant (R671W) adding a clinically relevant dimension to the study.

Mechanistic Insights:

The study offers valuable mechanistic insights into the development of microcephaly by elucidating the role of Abba in radial glial cell proliferation, radial fiber organization, and the migration of neuronal progenitors. The identification of Abba's involvement in the cleavage furrow during cell division, along with its interaction with Nedd9 and positive influence on RhoA activity, adds depth to our understanding of the molecular processes governing cortical development. Though the reported results establish the novel interaction between Abba and Nedd9, the authors have not addressed whether the mutant protein loses this interaction and whether that results in the observed effects.

In Vivo Validation:

The overexpression of mutant Abba protein (R671W) resulting in phenotypic similarities to Abba knockdown effects supports the significance of Abba in cortical development.

Reviewer #2 (Public Review):

Summary:

Carabalona and colleagues investigated the role of the membrane-deforming cytoskeletal regulator protein Abba (MTSS1L/MTSS2) in cortical development to better understand the mechanisms of abnormal neural stem cell mitosis. The authors used short hairpin RNA targeting Abba20 with a fluorescent reporter coupled with in-utero electroporation of E14 mice to show changes to neural progenitors. They performed flow cytometry for in-depth cell cycle analysis of Abba-shRNA impact on neural progenitors and determined an accumulation in the S phase. Using culture rat glioma cells and live imaging from cortical organotypic slides from mice in utero electroporated with Abba-shRNA, the authors found Abba played a prominent role in cytokinesis. They then used a yeast-two-hybrid screen to identify three high-confidence interactors: Beta-Trcp2, Nedd9, and Otx2. They used immunoprecipitation experiments from E18 cortical tissue coupled with C6 cells to show Abba's requirement for Nedd9 localization to the cleavage furrow/cytokinetic bridge. The authors performed a shRNA knockdown of Nedd9 by in-utero electroporation of E14 mice and observed similar results as with the Abba-shRNA. They tested a human variant of Abba using in-utero electroporation of cDNA and found disorganized radial glial fibers and misplaced, multipolar neurons, but lacked the impact of cell division seen in the shRNA-Abba model.

Strengths:

A fundamental question in biology about the mechanics of neural stem cell division.

Directly connecting effects in Abba protein to downstream regulation of RhoA via Nedd9.

Incorporation of human mutation in ABBA gene.

Use of novel technologies in neurodevelopment and imaging.

Weaknesses:

Unexplored components of the pathway (such as what neurogenic populations are impacted by Abba mutation) and unleveraged aspects of their data (such as the live imaging) limit the scope of their findings and leave significant questions about the effect of ABBA on radial glia development.

(1) The claim of disorganized radial glial fibers lacks quantifications.<br /> On page 11, the authors claim that knockdown of Abba leads to changes in radial glial morphology observed with vimentin staining. Here they claim misoriented apical processes, detached end feet, and decreased number of RGP cells in the VZ. However, they do not provide quantification of process orientation to better support their first claim. Measurements of radial glia fiber morphology (directionality, length) and angle of division would be metrics that can be applied to data. Some of these analyses could be done in their time-lapse microscopy images, such as to quantify the number of cell divisions during their period of analysis (though that is short-15 hours).

(2) It is unclear where the effect is:

-In RG or neuroblasts? Is it in cell cleavage that results in the accumulation of cells at VZ (as sometimes indicated by their data like in Figure 2A or 4D)? Interrogation of cell death (such as by cleaved caspase 3) would also help. Given their time-lapse, can they identify what is happening to the RG fiber? The authors describe a change in "migration" but do not show evidence for this for either progenitor or neuroblast populations. Given they have nice time-lapse imaging data, could they visualize progenitor versus young neuron migration? Analysis of neuroblasts (such as with doublecortin expression in the tissue) would also help understand any issues in migration (of neurons v stem cells).

-At cleavage furrow? In abscission? There is high-resolution data that highlights the cleavage furrow as the location of interest (Figure 3A), however, there is also data (Figure 3B) to suggest Abba is expressed elsewhere as well and there is an overall soma decrease. More detail of the localization of Abba during the division process would be helpful for example, could cleavage furrow proteins, such as Aurora B, co-localization (and potentially co-IP) help delineate subpopulations of Abba protein? Furthermore, the FRET imaging is a unique way to connect their mutation with function - could they measure/quantify differences at furrow compared to the rest of soma to further corroborate that the Abba-associated RhoA effect was furrow-enriched?

-The data highlights nicely that a furrow doesn't clearly form when ABBA expression and subsequent RhoA activity are decreased (in Figure 3 or 5A). Does this lead to cells that can't divide because of poor abscission, especially since "rounding" still occurs? Or abnormal progenitors (with loss of fiber or inability to support neuroblast migration)? Or abnormal progression of progenitors to neuroblasts?

(3) Limited to a singular time point of mouse cortical development

On page 13, the authors outline the results of their Y2H screen with the identification of three high-confidence interactors. Notably, they used an E10.5-E12.5 mouse brain embryo library rather than one that includes E14, the age of their in-utero electroporation mice. Many of the authors' claims focus on in-utero electroporation of shRNA-Abba of E14 mice that are then evaluated at E16-18. Justification for the focus on this age range should be included to support that their findings can then be applied to all mouse corticogenesis.

(4) Detail of the effect of the human variant of the ABBA mutation in mice is lacking.

Their identification of the R671W mutation is interesting and the IUE model warrants more characterization, as they did with their original KD experiments.

-Could they show that Abba protein levels are decreased (in either cell lines or electroporated tissue)?

-While time-lapse morphology might not have been performed, more analysis on cell division phenotype (such as plane of division and radial glia morphology) would be helpful.

Reviewer #1 (Public Review):

Summary:

This work proposes a new method, DyNetCP, for inferring dynamic functional connectivity between neurons from spike data. DyNetCP is based on a neural network model with a two-stage model architecture of static and dynamic functional connectivity.

This work evaluates the accuracy of the synaptic connectivity inference and shows that DyNetCP can infer the excitatory synaptic connectivity more accurately than a state-of-the-art model (GLMCC) by analyzing the simulated spike trains. Furthermore, it is shown that the inference results obtained by DyNetCP from large-scale in-vivo recordings are similar to the results obtained by the existing methods (jitter-corrected CCG and JPSTH). Finally, this work investigates the dynamic connectivity in the primary visual area VISp and in the visual areas using DyNetCP.

Strengths:

The strength of the paper is that it proposes a method to extract the dynamics of functional connectivity from spike trains of multiple neurons. The method is potentially useful for analyzing parallel spike trains in general, as there are only a few methods (e.g. Aertsen et al., J. Neurophysiol., 1989, Shimazaki et al., PLoS Comput Biol 2012) that infer the dynamic connectivity from spikes. Furthermore, the approach of DyNetCP is different from the existing methods: while the proposed method is based on the neural network, the previous methods are based on either the descriptive statistics (JSPH) or the Ising model.

Weaknesses:

Although the paper proposes a new method, DyNetCP, for inferring the dynamic functional connectivity, its strengths are neither clear nor directly demonstrated in this paper. That is, insufficient analyses are performed to support the usefulness of DyNetCP.

First, this paper attempts to show the superiority of DyNetCP by comparing the performance of synaptic connectivity inference with GLMCC (Figure 2). However, the improvement in the synaptic connectivity inference does not seem to be convincing. While this paper compares the performance of DyNetCP with a state-of-the-art method (GLMCC), there are several problems with the comparison. For example:

(1) This paper focused only on excitatory connections (i.e., ignoring inhibitory neurons).

(2) This paper does not compare with existing neural network-based methods (e.g., CoNNECT: Endo et al. Sci. Rep. 2021; Deep learning: Donner et al. bioRxiv, 2024).

(3) Only a population of neurons generated from the Hodgkin-Huxley model was evaluated.

Thus, the results in this paper are not sufficient to conclude the superiority of DyNetCP in the estimation of synaptic connections. In addition, this paper compares the proposed method with the standard statistical methods Jitter-corrected CCG (Figure 3) and JPSTH (Figure 4). Unfortunately, these results do not show the superiority of the proposed method. It only shows that the results obtained by the proposed method are consistent with those obtained by the existing methods (CCG or JPSTH). This paper also compares the proposed method with standard statistical methods, such as jitter-corrected CCG (Figure 3) and JPSTH (Figure 4). It only shows that the results obtained by the proposed method are consistent with those obtained by the existing methods (CCG or JPSTH), which does not show the superiority of the proposed method.

In summary, although DyNetCP has the potential to infer synaptic connections more accurately than existing methods, the paper does not provide sufficient analysis to make this claim. It is also unclear whether the proposed method is superior to the existing methods for estimating functional connectivity, such as jitter-corrected CCG and JPSTH. Thus, the strength of DyNetCP is unclear.

Reviewer #2 (Public Review):

Summary:

Here the authors describe a model for tracking time-varying coupling between neurons from multi-electrode spike recordings. Their approach extends a GLM with static coupling between neurons to include dynamic weights, learned by a long-short-term-memory (LSTM) model. Each connection has a corresponding LSTM embedding and is read out by a multi-layer perceptron to predict the time-varying weight.

Strengths:

This is an interesting approach to an open problem in neural data analysis. I think, in general, the method would be interesting to computational neuroscientists.

Weaknesses:

It is somewhat difficult to interpret what the model is doing. I think it would be worthwhile to add some additional results that make it more clear what types of patterns are being described and how.

Major Issues:

Simulation for dynamic connectivity. It certainly seems doable to simulate a recurrent spiking network whose weights change over time, and I think this would be a worthwhile validation for this DyNetCP model. In particular, I think it would be valuable to understand how much the model overfits, and how accurately it can track known changes in coupling strength. If the only goal is "smoothing" time-varying CCGs, there are much easier statistical methods to do this (c.f. McKenzie et al. Neuron, 2021. Ren, Wei, Ghanbari, Stevenson. J Neurosci, 2022), and simulations could be useful to illustrate what the model adds beyond smoothing.

Stimulus vs noise correlations. For studying correlations between neurons in sensory systems that are strongly driven by stimuli, it's common to use shuffling over trials to distinguish between stimulus correlations and "noise" correlations or putative synaptic connections. This would be a valuable comparison for Figure 5 to show if these are dynamic stimulus correlations or noise correlations. I would also suggest just plotting the CCGs calculated with a moving window to better illustrate how (and if) the dynamic weights differ from the data.

Reviewer #2 (Public Review):

Summary:

In this manuscript, using Staphylococcus aureus as a model organism, Panda et al. aim to understand how organic acids inhibit bacterial growth. Through careful characterization and interdisciplinary collaboration, the authors present valuable evidence that acetic acid specifically inhibits the activity of Ddl enzyme that converts 2 D-alanine amino acids into D-ala-D-ala dipeptide, which is then used to generate the stem pentapeptide of peptidoglycan (PG) precursors in the cytoplasm. Thus, a high concentration of acetic acid weakens the cell wall by limiting PG-crosslinking (which requires a D-ala portion). However, S. aureus maintains a high intracellular D-ala concentration to circumvent acetate-mediated growth inhibition.

Strengths:

The authors utilized a well-established transposon mutant library to screen for mutants that struggle to grow in the presence of acetic acid. This screen allowed authors to identify that a strain lacking intact alr1, which encodes for alanine racemase (converts L-ala to D-ala), is unable to grow well in the presence of acetic acid. This phenotype is rescued by the addition of external D-ala. Next, the authors rule out the contribution of other pathways that could lead to the production of D-ala in the cell. Finally, by analyzing D-ala and D-ala-D-ala concentrations, as well as muropeptide intermediates accumulation in different mutants, the authors pinpoint Ddl as the specific target of acetic acid. In fact, the synthetic overexpression of ddl alone overcomes the toxic effects of acetic acid. Using genetics, biochemistry, and structural biology, the authors show that Ddl activity is specifically inhibited by acetic acid and likely by other biologically relevant organic acids. Interestingly, this mechanism is different from what has been reported for other organisms such as Escherichia coli (where methionine synthesis is affected). It remains to be seen if this mechanism is conserved in other organisms that are more closely related to S. aureus, such as Clostridioides difficile and Enterococcus faecalis.

Weaknesses:

Although the authors have conclusively shown that Ddl is the target of acetic acid, it appears that the acetic acid concentration used in the experiments may not truly reflect the concentration range S. aureus would experience in its environment. Moreover, Ddl is only significantly inhibited at a very high acetate concentration (>400 mM). Thus, additional experiments showing growth phenotypes at lower organic acid concentrations may be beneficial. Another aspect not adequately discussed is the presence of D-ala in the gut environment, which may be protective against acetate toxicity based on the model provided.

Reviewer #1 (Public Review):

Summary:

The manuscript entitled "Staphylococcus aureus counters organic acid anion-mediated inhibition of peptidoglycan cross-linking through robust alanine racemase activity" by Panda, S et al. reports an extensive biochemical analysis of the result from a Tn screen that identified alr1 as being required for acetic acid tolerance. In the end, they demonstrate that reduced D-Ala pools in the ∆alr1 mutant lead to a drastic reduction in D-Ala-D-Ala dipeptide. They show that this is due to the ability of organic acid anions to limit the D-Ala-D-Ala ligase enzyme Ddl. They demonstrate that:

(1) Acetate exposure in the ∆alr1 results in reduced D-Ala-D-Ala dipeptide, but not the monomers.

(2) Acetate can bind to purified Ddl in vitro.

(3) This binding results in reduced enzyme activity.

(4) Other organic acid anions such as lactate, proprionate, and itaconitate can also inhibit Ddl.

The experiments are clearly described and logically laid out. I have only a few minor comments to add.

Strengths:

The most significant strength is the exceptional experimental data that supports the authors' hypotheses.

Weaknesses:

Only minor weaknesses were identified by this reviewer.

(1) Which allele is alr1, the one upstream of MazEF or the one in the Lysine biosynthetic operon?

(2) Figure 3B. Where does the C3N2 species come from in the WT and why is it absent in the mutants? It is about 25% of the total dipeptide pool.

(3) Figure 3D could perhaps be omitted. I understand that the authors attained statistical significance in the fitness defect, but biologically this difference is very minor. One would have to look at the isotopomer distribution in the Dat overexpressing strain to make sure that increased flux actually occurred since there are other means of affecting activity (e.g. allosteric modulators).

(4) In Figure 4A, why is the complete subunit UDP-NAM-AEKAA increasing in each strain upon acetate challenge if there was such a stark reduction in D-Ala-D-Ala, particularly in the ∆alr1 mutant? For that matter, why are the levels of UDP-NAM-AEKAA in the ∆alr1 mutant identical to that of WT with/out acetate?

(5) Figure 4B. Is there no significant difference between ddl and murF transcripts between WT and ∆alr1 under acetate stress? This comparison was not labeled if the tests were done.

(6) Although tricky, it is possible to measure intracellular acetate. It might be of interest to know where in the Ddl inhibition curve the cells actually are.

Reviewer #1 (Public Review):

A subclass of inhibitory heterotrimeric guanine nucleotide-binding protein subunits, GNAI, has been implicated in sensory hair cell formation, namely the establishment of hair bundle (stereocilia) orientation and staircase formation. However, the former role of hair bundle orientation has only been demonstrated in mutants expressing pertussis toxin, which blocks all GNAI subunits, but not in mutants with a single knockout of any of the Gnai genes, suggesting that there is a redundancy among various GNAI proteins in this role. Using various conditional mutants, the authors concluded that GNAI3 is the primary GNAI proteins required for hair bundle morphogenesis, whereas hair bundle orientation requires both GNAI2 and GNAI3.

Strength

Various compound mutants were generated to decipher the contribution of individual GNAI1, GNAI2, GNAI3 and GNAIO in the establishment of hair bundle orientation and morphogenesis. The study is thorough with detailed quantification of hair bundle orientation and morphogenesis, as well as auditory functions.

The revised manuscript has clarified the phenotypic differences raised between the Gnai2/3 double mutants and Ptx mutant phenotypes and resolved the weakness pointed out in the previous submission. These results further illustrate the dynamic requirement of Gnai/O in hair bundle establishment and is an important contribution to the field.

Reviewer #2 (Public Review):

Jarysta and colleagues set out to define how similar GNAI/O family members contribute to the shape and orientation of stereocilia bundles on auditory hair cells. Previous work demonstrated that loss of particular GNAI proteins, or inhibition of GNAIs by pertussis toxin, caused several defects in hair bundle morphogenesis, but open questions remained which the authors sought to address. Some of these questions include whether all phenotypes resulting from expression of pertussis toxin stemmed from GNAI inhibition; which GNAI family members are most critical for directing bundle development; whether GNAI proteins are needed for basal body movements that contribute to bundle patterning. These questions are important for understanding how tissue is patterned in response to planar cell polarity cues.

To address questions related to the GNAI family in auditory hair cell development, the authors assembled an impressive and nearly comprehensive collection of mouse models. This approach allowed for each Gnai and Gnao gene to be knocked out individually or in combination with each other. Notably, a new floxed allele was generated for Gnai3 because loss of this gene in combination with Gnai2 deletion was known to be embryonic lethal. Besides these lines, a new knockin mouse was made to conditionally express untagged pertussis toxin following cre induction from a strong promoter. The breadth and complexity involved in generating and collecting these strains makes this study unique, and likely the authoritative last word on which GNAI proteins are needed for which aspect of auditory hair bundle development.

Appropriate methods were employed by the authors to characterize auditory hair bundle morphology in each mouse line. Conclusions were carefully drawn from the data and largely based on excellent quantitative analysis. The main conclusions are that GNAI3 has the largest effect on hair bundle development. GNAI2 can compensate for GNAI3 loss in early development but incompletely in late development. The Gnai2 Gnai3 double mutant recapitulates nearly all the phenotypic effects associated with pertussis toxin expression and also reveals a role for GNAIs in early movement of the basal body. This comprehensive study builds on earlier reports, both uncovering new functions and putting previously putative functions on solid ground.

Reviewer #2 (Public Review):

In this manuscript, the authors analyze the shapes of cerebral cortices from several primate species, including subgroups of young and old humans, to characterize commonalities in patterns of gyrification, cortical thickness, and cortical surface area. The authors state that the observed scaling law shares properties with fractals, where shape properties are similar across several spatial scales. One way the authors assess this is to perform a "cortical melting" operation that they have devised on surface models obtained from several primate species. The authors also explore differences in shape properties between brains of young (~20 year old) and old (~80) humans. A challenge the authors acknowledge struggling with in reviewing the manuscript is merging "complex mathematical concepts and a perplexing biological phenomenon." This reviewer remains a bit skeptical about whether the complexity of the mathematical concepts being drawn from are justified by the advances made in our ability to infer new things about the shape of the cerebral cortex.

(1) The series of operations to coarse-grain the cortex illustrated in Figure 1 produces image segmentations that do not resemble real brains. The process to assign voxels in downsampled images to cortex and white matter is biased towards the former, as only 4 corners of a given voxel are needed to intersect the original pial surface, but all 8 corners are needed to be assigned a white matter voxel. The reason for introducing this bias (and to the extent that it is present in the authors' implementation) is not provided. The authors provide an intuitive explanation of why thickness relates to folding characteristics, but ultimately an issue for this reviewer is, e.g., for the right-most panel in Figure 2b, the cortex consists of several 4.9-sided voxels and thus a >2 cm thick cortex. A structure with these morphological properties is not consistent with the anatomical organization of typical mammalian neocortex.

(2) For the comparison between 20-year-old and 80-year-old brains, a well-documented difference is that the older age group possesses more cerebral spinal fluid due to tissue atrophy, and the distances between the walls of gyri becomes greater. This difference is born out in the left column of Figure 4b. It seems this additional spacing between gyri in 80 year olds requires more extensive down-sampling (larger scale values in Figure 4a) to achieve a similar shape parameter K as for the 20 year olds. The authors assert that K provides a more sensitive measure (associated with a large effect size) than currently used ones for distinguishing brains of young vs. old people. A more explicit, or elaborate, interpretation of the numbers produced in this manuscript, in terms of brain shape, might make this analysis more appealing to researchers in the aging field.

(3) In the Discussion, it is stated that self-similarity, operating on all length scales, should be used as a test for existing and future models of gyrification mechanisms. Given the lack of association between the abstract mathematical parameters described in this study and explicit properties of brain tissue and its constituents, it is difficult to envision how the coarse-graining operation can be used to guide development of "models of cortical gyrification."

(4) There are several who advocate for analyzing cortical mid-thickness surfaces, as the pial surface over-represents gyral tips compared to the bottoms of sulci in the surface area. The authors indicate that analyses of mid-thickness representations will be taken on in future work, but this seems to be a relevant control for accepting the conclusions of this manuscript.

Reviewer #3 (Public Review):

Summary: Through a rigorous methodology, the authors demonstrated that within 11 different primates, the shape of the brain followed a universal scaling law with fractal properties. They enhanced the universality of this result by showing the concordance of their results with a previous study investigating 70 mammalian brains, and the discordance of their results with other folded objects that are not brains. They incidentally illustrated potential applications of this fractal property of the brain by observing a scale-dependant effect of aging on the human brain.

Strengths:<br /> - New hierarchical way of expressing cortical shapes at different scales derived from previous report through implementation of a coarse-graining procedure<br /> - Investigation of 11 primate brains and contextualisation with other mammals based on prior literature<br /> - Proposition of tool to analyse cortical morphology requiring no fine tuning and computationally achievable<br /> - Positioning of results in comparison to previous works reinforcing the validity of the observation.<br /> - Illustration of scale-dependance of effects of brain aging in the human.

Weaknesses:<br /> - The notion of cortical shape, while being central to the article, is not really defined, leaving some interpretation to the reader<br /> - The organization of the manuscript is unconventional, leading to mixed contents in different sections (sections mixing introduction and method, methods and results, results and discussion...). As a result, the reader discovers the content of the article along the way, it is not obvious at what stages the methods are introduced, and the results are sometimes presented and argued in the same section, hindering objectivity.<br /> To improve the document, I would suggest a modification and restructuring of the article such that: 1) by the end of the introduction the reader understands clearly what question is addressed and the value it holds for the community, 2) by the end of the methods the reader understands clearly all the tools that will be used to answer that question (not just the new method), 3) by the end of the results the reader holds the objective results obtained by applying these tools on the available data (without subjective interpretations and justifications), and 4) by the end of the discussion the reader understands the interpretation and contextualisation of the study, and clearly grasps the potential of the method depicted for the better understanding of brain folding mechanisms and properties.

Reviewer #1 (Public Review):

In this study, Lin et al developed a protocol termed HIF-Clear, to perform tissue clearing and labelling on large-scale FFPE mouse brain specimens. They have optimized protocols for dewaxing and adequate delipidation of FFPE tissues to enable deep immunolabelling, even for whole mouse brains. This was useful for the study of disease models such as in an astrocytoma model to evaluate spatial architecture of the tumour and its surrounding microenvironment. It was also used in a traumatic brain injury model to quantify changes in vasculature density and differences in monoaminergic innervation. They have also demonstrated the potential of multi-round immunolabelling using photobleaching, as well as expansion microscopy with FFPE samples using Hif Clear.

Comments on revised version:

The revised manuscript by Lin et al is much improved with a more detailed methods description. There are only a few minor comments for the authors that are still valid:

- Some procedures, including the basic HIF-Clear protocol, seem to produce marked tissue expansion that is not mentioned in the manuscript. Users should take this fact into consideration when making measurements.<br /> - The authors have provided a comparison between mouse and human brain samples in Figure S12. However, it is misleading to mention that the "fluorescent signals are comparable at varying depth" as the figure clearly showed a lack of continuous staining especially for SMI312 at 900um depth, and human brain tissue showed considerably increased background signal (likely due to endogenous lipofuscin which has autofluorescent properties). Also, This is difficult to assess in the present design of the experiment because, at different depths, the tissue and the antigen may change themselves... making it difficult to make a direct staining comparison with other depths.

Reviewer #2 (Public Review):

The manuscript details an investigation aimed at developing a protocol to render centimeter-scale formalin-fixed paraffin-embedded specimens optically transparent and suitable for deep immunolabeling. The authors evaluate various detergents and conditions for epitope retrieval such as acidic or basic buffers combined with high temperatures in entire mouse brains that had been paraffin-embedded for months. They use various protein targets to test active immunolabeling and light-sheet microscopy registration of such preparations to validate their protocol. The final procedure, called MOCAT pipeline, briefly involves 1% Tween 20 in citrate buffer, heated in a pressure cooker at 121 {degree sign}C for 10 minutes. The authors also note that part of the delipidation is achieved by the regular procedure.

Major Strengths<br /> - The simplicity and ease of implementation of the proposed procedure using common laboratory reagents distinguish it favorably from more complex methods.

- Direct comparisons with existing protocols and exploration of alternative conditions enhance the robustness and practicality of the methodology.

Final considerations<br /> The evidence presented supports the effectiveness of the proposed method in rendering thick FFPE samples transparent and facilitating repeated rounds of immunolabeling.

The developed procedure holds promise for advancing tissue and 3D-specific determination of proteins of interest in various settings, including hospitals, basic research, and clinical labs, particularly benefiting neuroscience research.

The methodological findings suggest that MOCAT could have broader applications beyond FFPE samples, differentiating it from other tissue-clearing approaches in that the equipment and chemicals needed are broadly accessible.

Reviewer #1 (Public Review):

Summary:

Working memory is imperfect - memories accrue errors over time and are biased towards certain identities. For example, previous work has shown memory for orientation is more accurate near the cardinal directions (i.e., variance in responses is smaller for horizontal and vertical stimuli) while being biased towards diagonal orientations (i.e., there is a repulsive bias away from horizontal and vertical stimuli). The magnitude of errors and biases increase the longer an item is held in working memory and when more items are held in working memory (i.e., working memory load is higher). Previous work has argued that biases and errors could be explained by increased perceptual acuity at cardinal directions. However, these models are constrained to sensory perception and do not explain how biases and errors increase over time in memory. The current manuscript builds on this work to show how a two-layer neural network could integrate errors and biases over a memory delay. In brief, the model includes a 'sensory' layer with heterogenous connections that lead to the repulsive bias and decreased error in the cardinal directions. This layer is then reciprocally connected with a classic ring attractor layer. Through their reciprocal interactions, the biases in the sensory layer are constantly integrated into the representation in memory. In this way, the model captures the distribution of biases and errors for different orientations that have been seen in behavior and their increasing magnitude with time. The authors compare the two-layer network to a simpler one-network model, showing that the one-model network is harder to tune and shows an attractive bias for memories that have lower error (which is incompatible with empirical results).

Strengths:

The manuscript provides a nice review of the dynamics of items in working memory, showing how errors and biases differ across stimulus space. The two-layer neural network model is able to capture the behavioral effects as well as relate to neurophysiological observations that memory representations are distributed across the sensory cortex and prefrontal cortex.

The authors use multiple approaches to understand how the network produces the observed results. For example, analyzing the dynamics of memories in the low-dimensional representational space of the networks provides the reader with an intuition for the observed effects.

As a point of comparison with the two-layer network, the authors construct a heterogenous one-layer network (analogous to a single memory network with embedded biases). They argue that such a network is incapable of capturing the observed behavioral effects but could potentially explain biases and noise levels in other sensory domains where attractive biases have lower errors (e.g., color).

The authors show how changes in the strength of Hebbian learning of excitatory and inhibitory synapses can change network behavior. This argues for relatively stronger learning in inhibitory synapses, an interesting prediction.

The manuscript is well-written. In particular, the figures are well done and nicely schematize the model and the results.

Weaknesses:

Despite its strengths, the manuscript does have some weaknesses.

First, as far as we can tell, behavioral data is only presented in schematic form. This means some of the nuances of the effects are lost. It also means that the model is not directly capturing behavioral effects. Therefore, while providing insight into the general phenomenon, the current manuscript may be missing some important aspects of the data.

Relatedly, the models are not directly fit to behavioral data. This makes it hard for the authors to exclude the possibility that there is a single network model that could capture the behavioral effects. In other words, it is hard to support the authors' conclusion that "....these evolving errors...require network interaction between two distinct modules." (from the abstract, but similar comments are made throughout the manuscript). Such a strong claim needs stronger evidence than what is presented. Fitting to behavioral data could allow the authors to explore the full parameter space for both the one-layer and two-layer network architectures.

In addition, directly comparing the ability of different model architectures to fit behavioral data would allow for quantitative comparison between models. Such quantitative comparisons are currently missing from the manuscript.

To help broaden the impact of the paper, it would be helpful if the authors provided insight into how the observed behavioral biases and/or network structures influence cognition. For example, previous work has argued that biases may counteract noise, leading to decreased variance at certain locations. Is there a similar normative explanation for why the brain would have repulsive biases away from commonly occurring stimuli? Are they simply a consequence of improved memory accuracy? Why isn't this seen for all stimulus domains?

Previous work has found both diffusive noise and biases increase with the number of items in working memory. It isn't clear how the current model would capture these effects. The authors do note this limitation in the Discussion, but it remains unclear how the current model can be generalized to a multi-item case.

The role of the ring attractor memory network isn't completely clear. There is noise added in this stage, but how is this different from the noise added at the sensory stage? Shouldn't these be additive? Is the noise necessary? Similarly, it isn't clear whether the memory network is necessary - can it be replaced by autapses (self-connections) in the sensory network to stabilize its representation? In short, it would be helpful for the authors to provide an intuition for why the addition of the memory network facilitates the repulsive bias.

Overall:

Overall, the manuscript was successful in building a model that captured the biases and noise observed in working memory. This work complements previous studies that have viewed these effects through the lens of optimal coding, extending these models to explain the effects of time in memory. In addition, the two-layer network architecture extends previous work with similar architectures, adding further support to the distributed nature of working memory representations.

Reviewer #2 (Public Review):

In this manuscript, Yang et al. present a modeling framework to understand the pattern of response biases and variance observed in delayed-response orientation estimation tasks. They combine a series of modeling approaches to show that coupled sensory-memory networks are in a better position than single-area models to support experimentally observed delay-dependent response bias and variance in cardinal compared to oblique orientations. These errors can emerge from a population-code approach that implements efficient coding and Bayesian inference principles and is coupled to a memory module that introduces random maintenance errors. A biological implementation of such operation is found when coupling two neural network modules, a sensory module with connectivity inhomogeneities that reflect environment priors, and a memory module with strong homogeneous connectivity that sustains continuous ring attractor function. Comparison with single-network solutions that combine both connectivity inhomogeneities and memory attractors shows that two-area models can more easily reproduce the patterns of errors observed experimentally. This, the authors take as evidence that a sensory-memory network is necessary, but I am not convinced about the evidence in support of this "necessity" condition. A more in-depth understanding of the mechanisms operating in these models would be necessary to make this point clear.

Strengths:

The model provides an integration of two modeling approaches to the computational bases of behavioral biases: one based on Bayesian and efficient coding principles, and one based on attractor dynamics. These two perspectives are not usually integrated consistently in existing studies, which this manuscript beautifully achieves. This is a conceptual advancement, especially because it brings together the perceptual and memory components of common laboratory tasks.

The proposed two-area model provides a biologically plausible implementation of efficient coding and Bayesian inference principles, which interact seamlessly with a memory buffer to produce a complex pattern of delay-dependent response errors. No previous model had achieved this.

Weaknesses:

The correspondence between the various computational models is not fully disclosed. It is not easy to see this correspondence because the network function is illustrated with different representations for different models and the correspondence between components of the various models is not specified. For instance, Figure 1 shows that a specific pattern of noise is required in the low-dimensional attractor model, but in the next model in Figure 2, the memory noise is uniform for all stimuli. How do these two models integrate? What element in the population-code model of Figure 2 plays the role of the inhomogeneous noise of Figure 1? Also, the Bayesian model of Figure 2 is illustrated with population responses for different stimuli and delays, while the attractor models of Figures 3 and 4 are illustrated with neuronal tuning curves but not population activity. In addition, error variance in the Bayesian model appears to be already higher for oblique orientations in the first iteration whereas it is only first shown one second into the delay for the attractor model in Figure 4. It is thus unclear whether variance inhomogeneities appear already at the perceptual stage in the attractor model, as it does in the population-code model. Of course, correspondences do not need to be perfect, but the reader does not know right now how far the correspondence between these models goes.

The manuscript does not identify the mechanistic origin in the model of Figure 4 of the specific noise pattern that is required for appropriate network function (with higher noise variance at oblique orientations). This mechanism appears critical, so it would be important to know what it is and how it can be regulated. In particular, it would be interesting to know if the specific choice of Poisson noise in Equation (3) is important. Tuning curves in Figure 4 indicate that population activity for oblique stimuli will have higher rates than for cardinal stimuli and thus induce a larger variance of injected noise in oblique orientations, based on this Poisson-noise assumption. If this explanation holds, one wonders if network inhomogeneities could be included (for instance in neural excitability) to induce higher firing rates in the cardinal/oblique orientations so as to change noise inhomogeneities independently of the bias and thus control more closely the specific pattern of errors observed, possibly within a single memory network.

The main conclusion of the manuscript, that the observed patterns of errors "require network interaction between two distinct modules" is not convincingly shown. The analyses show that there is a quantitative but not a qualitative difference between the dynamics of the single memory area compared to the sensory-memory two-area network, for specific implementations of these models (Figure 7 - Figure Supplement 1). There is no principled reasoning that demonstrates that the required patterns of response errors cannot be obtained from a different memory model on its own. Also, since the necessity of the two-area configuration is highlighted as the main conclusion of the manuscript, it is inconvenient that the figure that carefully compares these conditions is in the Supplementary Material.

The proposed model has stronger feedback than feedforward connections between the sensory and memory modules. This is not a common assumption when thinking about hierarchical processing in the brain, and it is not discussed in the manuscript.

Reviewer #3 (Public Review):

Summary:

The present study proposes a neural circuit model consisting of coupled sensory and memory networks to explain the circuit mechanism of the cardinal effect in orientation perception which is characterized by the bias towards the oblique orientation and the largest variance at the oblique orientation.

Strengths:

The authors have done numerical simulations and preliminary analysis of the neural circuit model to show the model successfully reproduces the cardinal effect. And the paper is well-written overall. As far as I know, most of the studies on the cardinal effect are at the level of statistical models, and the current study provides one possibility of how neural circuit models reproduce such an effect.

Weaknesses:

There are no major weaknesses and flaws in the present study, although I suggest the author conduct further analysis to deepen our understanding of the circuit mechanism of the cardinal effects. Please find my recommendations for concrete comments.

Reviewer #1 (Public Review):

Summary:

A cortico-centric view is dominant in the study of the neural mechanisms of consciousness. This investigation represents the growing interest in understanding how subcortical regions are involved in conscious perception. To achieve this, the authors engaged in an ambitious and rare procedure in humans of directly recording from neurons in the subthalamic nucleus and thalamus. While participants were in surgery for the placement of deep brain stimulation devices for the treatment of essential tremor and Parkinson's disease, they were awakened and completed a perceptual-threshold tactile detection task. The authors identified individual neurons and analyzed single-unit activity corresponding with the task phases and tactile detection/perception. Among the neurons that were perception-responsive, the authors report changes in firing rate beginning ~150 milliseconds from the onset of the tactile stimulation. Curiously, the majority of the perception-responsive neurons had a higher firing rate for missed/not perceived trials. In summary, this investigation is a valuable addition to the growing literature on the role of subcortical regions in conscious perception.

Strengths:

The authors achieved the challenging task of recording human single-unit activity while participants performed a tactile perception task. The methods and statistics are clearly explained and rigorous, particularly for managing false positives and non-normal distributions. The results offer new detail at the level of individual neurons in the emerging recognition of the role of subcortical regions in conscious perception.

Weaknesses:

"Nonetheless, it remains unknown how the firing rate of subcortical neurons changes when a stimulus is consciously perceived." (lines 76-77) The authors could be more specific about what exactly single-unit recordings offer for interrogating the role of subcortical regions in conscious perception that is unique from alternative neural activity recordings (e.g., local field potential) or recordings that are used as proxies of neural activity (e.g., fMRI).

Related comment for the following excerpts:

"After a random delay ranging from 0.5 to 1 s, a "respond" cue was played, prompting participants to verbally report whether they felt a vibration or not. Therefore, none of the reported analyses are confounded by motor responses." (lines 97-99).

"These results show that subthalamic and thalamic neurons are modulated by stimulus onset, irrespective of whether it was reported or not, even though no immediate motor response was required." (lines 188-190).

"By imposing a delay between the end of the tactile stimulation window and the subjective report, we ensured that neuronal responses reflected stimulus detection and not mere motor responses." (lines 245-247).

It is a valuable feature of the paradigm that the reporting period was initiated hundreds of milliseconds after the stimulus presentation so that the neural responses should not represent "mere motor responses". However, verbal report of having perceived or not perceived a stimulus is a motor response and because the participants anticipate having to make these reports before the onset of the response period, there may be motor preparatory activity from the time of the perceived stimulus that is absent for the not perceived stimulus. The authors show sensitivity to this issue by identifying task-selective neurons and their discussion of the results that refer to the confound of post-perceptual processing. Still, direct treatment of this possible confound would help the rigor of the interpretation of the results.

"When analyzing tactile perception, we ensured that our results were not contaminated with spurious behavior (e.g. fluctuation of attention and arousal due to the surgical procedure)." (lines 118-117).

Confidence in the results would be improved if the authors clarified exactly what behaviors were considered as contaminating the results (e.g., eye closure, saccades, and bodily movements) and how they were determined.

The authors' discussion of the thalamic neurons could be more precise. The authors show that only certain areas of the thalamus were recorded (in or near the ventral lateral nucleus, according to Figure S3C). The ventral lateral nucleus has a unique relationship to tactile and motor systems, so do the authors hypothesize these same perception-selective neurons would be active in the same way for visual, auditory, olfactory, and taste perception? Moreover, the authors minimally interpret the location of the task, sensory, and perception-responsive neurons. Figure S3 suggests these neurons are overlapping. Did the authors expect this overlap and what does it mean for the functional organization of the ventral lateral nucleus and subthalamic nucleus in conscious perception?

"We note that, 6 out of 8 neurons had higher firing rates for missed trials than hit trials, although this proportion was not significant (binomial test: p = 0.145)." (lines 215-216).

It appears that in the three example neurons shown in Figure 4, 2 out of 3 (#001 and #068) show a change in firing rate predominantly for the missed stimulations. Meanwhile, #034 shows a clear hit response (although there is an early missed response - decreased firing rate - around 150 ms that is not statistically significant). This is a counterintuitive finding when compared to previous results from the thalamus (e.g., local field potentials and fMRI) that show the opposite response profile (i.e., missed/not perceived trials display no change or reduced response relative to hit/perceived trials). The discussion of the results should address this, including if these seemingly competing findings can be rectified.

The authors report 8 perception-responsive neurons, but there are only 5 recording sites highlighted (i.e., filled-in squares and circles) in Figures S3C and 4D. Was this an omission or were three neurons removed from the perception-responsive analysis?

Could the authors speak to the timing of the responses reported in Figure 4? The statistically significant intervals suggested both early (~160-200ms) to late responses (~300ms). Some have hypothesized that subcortical regions are early - ahead of cortical activation that may be linked with conscious perception. Do these results say anything about this temporal model for when subcortical regions are active in conscious perception?

Reviewer #2 (Public Review):

The authors have studied subpopulations of individual neurons recorded in the thalamus and subthalamic nucleus (STN) of awake humans performing a simple cognitive task. They have carefully designed their task structure to eliminate motor components that could confound their analyses in these subcortical structures, given that the data was recorded in patients with Parkinson's Disease (PD) and diagnosed with an Essential Tremor (ET). The recorded data represents a promising addition to the field. The analyses that the authors have applied can serve as a strong starting point for exploring the kinds of complex signals that can emerge within a single neuron's activity. Pereira et. al conclude that their results from single neurons indicate that task-related activity occurs, purportedly separate from previously identified sensory signals. These conclusions are a promising and novel perspective for how the field thinks about the emergence of decisions and sensory perception across the entire brain as a unit.

Despite the strength of the data that was obtained and the relevant nature of the conclusions that were drawn, there are certain limitations that must be taken into consideration:

(1) The authors make several claims that their findings are direct representations of consciousness identifiable in subcortical structures. The current context for consciousness does not sufficiently define how the consciousness is related to the perceptual task.

(2) The current work would benefit greatly from a description and clarification of what all the neurons that have been recorded are doing. The authors' criteria for selecting subpopulations with task-relevant activity are appropriate, but understanding the heterogeneity in a population of single neurons is important for broader considerations that are being studied within the field.

(3) The authors have omitted a proper set of controls for comparison against the active trials, for example, where a response was not necessary. Please explain why this choice was made and what implications are necessary to consider.

Reviewer #3 (Public Review):

Summary:

This important study relies on a rare dataset: intracranial recordings within the thalamus and the subthalamic nucleus in awake humans, while they were performing a tactile detection task. This procedure allowed the authors to identify a small but significant proportion of individual neurons, in both structures, whose activity correlated with the task (e.g. their firing rate changed following the audio cue signalling the start of a trial) and/or with the stimulus presentation (change in firing rate around 200 ms following tactile stimulation) and/or with participant's reported subjective perception of the stimulus (difference between hits and misses around 200 ms following tactile stimulation). Whereas most studies interested in the neural underpinnings of conscious perception focus on cortical areas, these results suggest that subcortical structures might also play a role in conscious perception, notably tactile detection.

Strengths:

There are two strongly valuable aspects in this study that make the evidence convincing and even compelling. First, these types of data are exceptional, the authors could have access to subcortical recordings in awake and behaving humans during surgery. Additionally, the methods are solid. The behavioral study meets the best standards of the domain, with a careful calibration of the stimulation levels (staircase) to maintain them around the detection threshold, and an additional selection of time intervals where the behavior was stable. The authors also checked that stimulus intensity was the same on average for hits and misses within these selected periods, which warrants that the effects of detection that are observed here are not confounded by stimulus intensity. The neural data analysis is also very sound and well-conducted. The statistical approach complies with current best practices, although I found that, in some instances, it was not entirely clear which type of permutations had been performed, and I would advocate for more clarity in these instances. Globally the figures are nice, clear, and well presented. I appreciated the fact that the precise anatomical location of the neurons was directly shown in each figure.

Weaknesses:

Some clarification is needed for interpreting Figure 3, top rows: in my understanding the black curve is already the result of a subtraction between stimulus present trials and catch trials, to remove potential drifts; if so, it does not make sense to compare it with the firing rate recorded for catch trials.

I also think that the article could benefit from a more thorough presentation of the data and that this could help refine the interpretation which seems to be a bit incomplete in the current version. There are 8 stimulus-responsive neurons and 8 perception-selective neurons, with only one showing both effects, resulting in a total of 15 individual neurons being in either category or 13 neurons if we exclude those in which the behavior is not good enough for the hit versus miss analysis (Figure S4A). In my opinion, it should be feasible to show the data for all of them (either in a main figure, or at least in supplementary), but in the present version, we get to see the data for only 3 neurons for each analysis. This very small selection includes the only neuron that shows both effects (neuron #001; which is also cue selective), but this is not highlighted in the text. It would be interesting to see both the stimulus-response data and the hit versus miss data for all 13 neurons as it could help develop the interpretation of exactly how these neurons might be involved in stimulus processing and conscious perception. This should give rise to distinct interpretations for the three possible categories. Neurons that are stimulus-responsive but not perception-selective should show the same response for both hits and misses and hence carry out indifferently conscious and unconscious responses. The fact that some neurons show the opposite pattern is particularly intriguing and might give rise to a very specific interpretation: if the neuron really doesn't tend to respond to the stimulus when hits and misses are put together, it might be a neuron that does not directly respond to the stimulus, but whose spontaneous fluctuations across trials affect how the stimulus is perceived when they occur in a specific time window after the stimulus. Finally, neuron #001 responds with what looks like a real burst of evoked activity to stimulation and also shows a difference between hits and misses, but intriguingly, the response is strongest for misses. In the discussion, the interesting interpretation in terms of a specific gating of information by subcortical structures seems to apply well to this last example, but not necessarily to the other categories.

Reviewer #1 (Public Review):

Summary:

People with Parkinson's disease often experience a variety of nonmotor symptoms, the biological bases of which remain poorly understood. Johansson et al began to study potential roles of the dorsal raphe nucleus (DRN) degeneration in the pathophysiology of neuropsychiatric symptoms in PD.

Strengths:

Boi et al validated a transgenic reporter mouse line that can reliably label dopaminergic neurons in the DRN. This brain region shows severe neurodegeneration and has been proposed to contribute to the manifestation of neuropsychiatric symptoms in PD. Using this mouse line (and others), Boi and colleagues characterized electrophysiological and morphological phenotypes of dopaminergic and serotoninergic neurons in the raphe nucleus. This study involved very careful topographical registration of recorded neurons to brain slices for post hoc immunohistochemical validation of cell identity, making it an elegant and thorough piece of work.

In relevance to PD pathophysiology, the authors evaluated the physiological and morphological changes of DRN serotoninergic and dopaminergic neurons after a partial loss of nigrostriatal dopamine neurons, which serves as a mouse model of early parkinsonian pathology. Moreover, the authors identified a series of physiological and morphological changes of subtypes of DRN neurons that depend on nigral dopaminergic neurodegeneration, LC noradrenergic neurodegeneration, or both. Indeed this works highlights the importance of LC noradrenergic degeneration in PD pathophysiology.

Overall, this is a well-designed study with high significance to the Parkinson's research field.

Reviewer #2 (Public Review):

In this paper, Boi et al. thoroughly classified the electrophysiological and morphological characteristics of serotonergic and dopaminergic neurons in the DRN and examined the alterations of these neurons in the 6-OHDA-induced mouse PD model. Using whole-cell patch clamp recording, they found that 5-HT and dopamine (DA) neurons in the DRN are electrophysiologically well-distinguished from each other. In addition, they characterized distinct morphological features of 5-HT and DA neurons in the DRN. Notably, these specific features of 5-HT and DA neurons in the DRN exhibited different changes in the 6-OHDA-induced PD model. Then the authors utilized desipramine (DMI) to separate the effects of nigrostriatal DA depletion and noradrenalin (NA) depletion which are induced by 6-OHDA. Interestingly, protection from NA depletion by DMI pretreatment reversed the changes in 5-HT neurons, while having a minor impact on the changes in DA neurons in the DRN. These data indicate that the role of NA lesion in the altered properties of DRN 5-HT neurons by 6-OHDA is more critical than the one of DA lesion.

Overall, this study provides foundational data on the 5-HT and DA neurons in the DRN and their potential involvement in PD symptoms. Given the defects of the DRN in PD, this paper may offer insights into the cellular mechanisms that may underlie non-motor symptoms associated with PD. Despite the importance of the primary claim proposed by the authors, however, the interpretation of the authors on some DMI experiments is not explained well.

Reviewer #3 (Public Review):

Summary:

Using ex vivo electrophysiology and morphological analysis, Boi et al. investigate the electrophysiological and morphological properties of serotonergic and dopaminergic subpopulations in the dorsal raphe nucleus (DRN). They performed labor-intensive and rigorous electrophysiology with posthoc immunohistochemistry and neuronal reconstruction to delineate the two major cell classes in the DRN: DRN-DA and DRN-5HT, named according to their primary neurotransmitter machinery. They find that the dopaminergic (DRN-DA) and serotonergic (DRN-5HT) neurons are electrophysiologically and morphologically distinct, and are altered following striatal injection of the toxin 6-OHDA. However, these alterations were largely prevented in DRN-5HT neurons by pre-treatment with desipramine. These findings suggest an important interplay between catecholaminergic systems in healthy and parkinsonian conditions, as well as a relationship between neuronal structure and function.

Strengths:

Large, well-validated dataset that will be a resource for others.<br /> Complementary electrophysiological and anatomical characterizations.<br /> Conclusions are justified by the data.<br /> Relevant for basic scientists interested in DRN cell types and physiology<br /> Relevant for those interested in serotonin and/or DRN neurons in Parkinson's Disease

Weaknesses:

Given the scope of the author's questions and hypotheses, I did not identify any major weaknesses.

Reviewer #1 (Public Review):

Summary:

This is a study on the role of the retrosplenial cortex (RSC) and the hippocampus in working memory. Working memory is a critical cognitive function that allows temporary retention of information for task execution. The RSC, which is functionally and anatomically connected to both primary sensory (especially visual) and higher cognitive areas, plays a key role in integrating spatial-temporal context and in goal-directed behaviors. However, the specific contributions of the RSC and the hippocampus in working memory-guided behaviors are not fully understood due to a lack of studies that experimentally disrupt the connection between these two regions during such behaviors.

In this study, researchers employed eArch3.0 to silence hippocampal axon terminals in the RSC, aiming to explore the roles of these brain regions in working memory. Experiments were conducted where animals with silenced hippocampal axon terminals in the RSC performed a delayed non-match to place (DNMP) task. The results indicated that this manipulation impaired memory retrieval, leading to decreased performance and quicker decision-making in the animals. Notably, the authors observed that the effects of this impairment persisted beyond the light-activation period of the opsin, affecting up to three subsequent trials. They suggest that disrupting the hippocampal-RSC connection has a significant and lasting impact on working memory performance.

Strengths:

They conducted a study exploring the impact of direct hippocampal inputs into the RSC, a region involved in encoding spatial-temporal context and transferring contextual information, on spatial working memory tasks. Utilizing eArch3.0 expressed in hippocampal neurons via the viral vector AAV5-hSyn1-eArch3.0, they aimed to bilaterally silence hippocampal terminals located at the RSC in rats pre-trained in a DNMP task. They discovered that silencing hippocampal terminals in the RSC significantly decreased working memory performance in eArch+ animals, especially during task interleaving sessions (TI) that alternated between trials with and without light delivery. This effect persisted even in non-illuminated trials, indicating a lasting impact beyond the periods of direct manipulation. Additionally, they observed a decreased likelihood of correct responses following TI trials and an increased error rate in eArch+ animals, even after incorrect responses, suggesting an impairment in error-corrective behavior. This contrasted with baseline sessions where no light was delivered, and both eArch+ and control animals showed low error rates.

Weaknesses:

While I agree with the authors that the role of hippocampal inputs to the RSC in spatial working memory is understudied and merits further investigation, I find that the optogenetic experiment, a core part of this manuscript that includes viral injections, could be improved. The effects were rather subtle, rendering some of the results barely significant and possibly too weak to support major conclusions. Additionally, no mechanistic investigation was conducted beyond referencing previous reports to interpret the core behavioral phenotypes.

Reviewer #2 (Public Review):

The authors examine the impact of optogenetic inhibition of hippocampal axon terminals in the retrosplenial cortex (RSP) during the performance of a working memory T-maze task. Performance on a delayed non-match-to-place task was impaired by such inhibition. The authors also report that inhibition is associated with faster decision-making and that the effects of inhibition can be observed over several subsequent trials. The work seems reasonably well done and the role of hippocampal projections to retrosplenial cortex in memory and decision-making is very relevant to multiple fields. However, the work should be expanded in several ways before one can make firm conclusions on the role of this projection in memory and behavior.

(1) The work is very singular in its message and the experimentation. Further, the impact of the inhibition on behavior is very moderate. In this sense, the results do not support the conclusion that the hippocampal projection to retrosplenial cortex is key to working memory in a navigational setting.

(2) There are no experiments examining other types of behavior or working memory. Given that the animals used in the studies could be put through a large number of different tasks, this is surprising. There is no control navigational task. There is no working memory test that is non-spatial. Such results should be presented in order to put the main finding in context.

(3) The actual impact of the inhibition on activity in RSP is not provided. While this may not be strictly necessary, it is relevant that the hippocampal projection to RSP includes, and is perhaps dominated by inhibitory inputs. I wonder why the authors chose to manipulate hippocampal inputs to RSP when the subiculum stands as a much stronger source of afferents to RSP and has been shown to exhibit spatial and directional tuning of activity. The points here are that we cannot be sure what the manipulation is really accomplishing in terms of inhibiting RSP activity (perhaps this explains the moderate impact on behavior) and that the effect of inhibiting hippocampal inputs is not an effective means by which to study how RSP is responsive to inputs that reflect environmental locations.

(4) The impact of inhibition on trials subsequent to the trial during which optical stimulation was actually supplied seems trivial. The authors themselves point to evidence that activation of the hyperpolarizing proton pump is rather long-lasting in its action. Further, each sample-test trial pairing is independent of the prior or subsequent trials. This finding is presented as a major finding of the work, but would normally be relegated to supplemental data as an expected outcome given the dynamics of the pump when activated.

(5) In the middle of the first paragraph of the discussion, the authors make reference to work showing RSP responses to "contextual information in egocentric and allocentric reference frames". The citations here are clearly deficient. How is the Nitzan 2020 paper at all relevant here?

(6) The manuscript is deficient in referencing and discussing data from the Smith laboratory that is similar. The discussion reads mainly like a repeat of the results section.

Reviewer #2 (Public Review):

The authors aimed at elucidating the development of high altitude polycythemia which affects mice and men staying in a hypoxic atmosphere at high altitude (hypobaric hypoxia; HH). HH causes increased erythropoietin production which stimulates the production of red blood cells. The authors hypothesize that increased production is only partially responsible for exaggerated red blood cell production, i.e. polycythemia, but that decreased erythrophagocytosis in the spleen contributes to high red blood cells counts.

The main strength of the study is the use of a mouse model exposed to HH in a hypobaric chamber. However, not all of the reported results are convincing due to some smaller effects which one may doubt to result in the overall increase in red blood cells as claimed by the authors. Moreover, direct proof for reduced erythrophagocytosis is compromised due to a strong spontaneous loss of labelled red blood cells, although effects of labelled E. coli phagocytosis are shown.

Comments on latest version:

The authors have partly addressed my comments.

(1) The response to my question regarding unchanged MCH is a kind of "hand waiving" - maybe it would require substantially more extensive work to clarify this issue

(2) The moderate if not marginal difference in normal vs splenectomy argues against a significant role of the spleen - even if the difference was slightly larger in HH

(3) There is still overinterpretation of data. My Q was: Is the reduced phagocytic capacity in Fig 4B significant? Response: "This is indicative of a diminished phagocytic capacity, particularly when contrasted<br /> with NN conditions." I guess that is a "no"

(4) I assume my question with respect to bi- or trivalent iron chelators was misunderstood.

In general, as indicated above, it is an interesting hypothesis which is corroborated by data in several instances. Maybe the scientific community should decide whether it is all in all conclusive.

Reviewer #3 (Public Review):

The manuscript by Yang et al. investigated in mice how hypobaric hypoxia can modify the RBC clearance function of the spleen, a concept that is of interest. Via interpretation of their data, the authors proposed a model that hypoxia causes an increase in cellular iron levels, possibly in RPMs, leading to ferroptosis, and downregulates their erythrophagocytic capacity.

Comments on revised version:

The manuscript has now improved with all the new data, supporting the model proposed by the authors. However, it remains not very easy to follow for the conclusions and experimental details. Some of the most important remaining comments are listed below:

(1) Lines 401-406 - The conclusions in this new fragment sound a bit overstated - the authors do not directly measure erytrophagocytosis capacity, only the total RBC parameters in the circulation. The increase is also very mild biologically between sham and splenectomized mice in HH conditions.

(2) scRNA seq data are still presented in a way that is very difficult to understand. The readers could not see from the graphics that macrophages are depleted. The clusters are not labelled - some clusters in the bin 'macrophahes+DC' seem actually to be more represented in Fig. 3E; Fig. 3F does not correspond to Fig. 3D. It would be maybe more informative to present like in Figure D side by side NN versus HH? The authors could consider moving the data from supplements that relate to RPMs to the main figure and making it consistent for the Clusters - eg, the authors show data for Cluster 0 in the supplement, and the same Cluster is not marked as macrophages in the main figure. This is quite difficult to follow.

(3) Figure 3G has likely mislabeled axis for F4/80 and CD11b - such mistakes should be avoided in a second revised version of the manuscript, and this data is now redundant with the data shown as new Figure 5A.

(4) The data from new Figure 4 should be better mentioned in the main body of the manuscript - all panels are mentioned twice in the text, first speaking about the decline of labelled RBCs and second referring to phagocytic capacity, whereas this figure only illustrates the decline of labelled RBCs, not directly phagocytic capacity of RPMs. What is lacking, as opposed to typical RBC life span assay, is the time '0' ('starting point') - this is particularly important as we can observe a big drop in labelled RBCs for eg 7 days between NN and HH group, actually implying increased removal of labelled RBCs within the first days of hypoxia exposure. What should be better labelled in this figure is that the proportion of RBCs are labelled RBCs not all RBCs (Y axis in individual panels). Overall, the new Figure 4 brings new data to the study, but how it is presented and discussed is not at the 'state-of-the-art' level (eg, missing the time '0') and is not very straightforward to the reader.

(5) In Figure 7, the experiments with Tuftsin are not very easy to follow, especially for the major conclusions. In panels A and B, the focus is the drug itself under NN conditions, with RBC removal as a readout. Then, in the next panels, the authors introduce HH, and then look at the F4/80 and iron staining. What was exactly the major point the authors wanted to make here?

(6) The data from Figure 8 are informative but do not address the individual cell types - eg, a drop in HO1 or FT may be due to the depletion of RPMs. An increase of TFR1 could be due to the retention of RBCs, the same as maybe labile iron. The data from PBMC are only very loosely linked to these phenotypes observed in the total spleen, and the reason for the regulation of the same proteins in PBMC might be different. It goes back to the data in Figure 3A-C, where also total splenocytes are investigated for their viability.

(7) Can the authors provide the data for the purity (eg cell surface markers) of their primary splenic macrophage cultures? Only ensuring that these are macrophages or addressing the readouts from Figure 8 in RPMs could link ferroptosis to RPMs under HH conditions.

(8) All the data are not presented as individual data points which is not widely applied in papers.

(9) No gating strategies are nicely illustrated or described.

Reviewer #1 (Public Review):

A descriptive manuscript investigating the ability of a peptide, implicated in development, to induce signalling responses indicative of immunity. The work clearly documents the ability of the synthetic peptide to induce these responses, and open future work to link this back to physiology.

Comments on revised version:

Congratulations to the authors for the improvements to the manuscript.

I still have reservations, as raised by other reviewers, about whether the outputs observed can definitively be classified as immune/defence outputs without assaying an impact upon microbial growth. Indeed, this is challenging to address as many of the outputs are shared by multiple pathways. This is especially the case here as the peptide could have different effects in different tissues or cells with different expression levels of the receptors (e.g. hypothetically - no expression = no effect; weak expression - cell wall loosening and susceptibility; high expression - strong response and 'defence' response). I do however appreciate that the authors have toned down some of the conclusions regarding the defence response and also they included further reference to outputs also being from developmental pathways.

Reviewer #3 (Public Review):

Previously, it has been shown the essential role of IDA peptide and HAESA receptor families in driving various cell separation processes such as abscission of flowers as a natural developmental process, of leaves as a defense mechanism when plants are under pathogenic attack or at the lateral root emergence and root tip cell sloughing. In this work, Olsson et al. show for the first time the possible role of IDA peptide in triggering plant innate immunity after the cell separation process occurred. Such an event has been previously proposed to take place in order to seal open remaining tissue after cell separation to avoid creating an entry point for opportunistic pathogens. The elegant experiments in this work demonstrate that IDA peptide is triggering the defense-associated marker genes together with immune specific responses including release of ROS and intracellular CA2+. Thus, the work highlights an intriguing direct link between endogenous cell wall remodeling and plant immunity. Moreover, the upregulation of IDA in response to abiotic and especially biotic stimuli are providing a valuable indication for potential involvement of HAE/IDA signalling in other processes than plant development.

Comments on revised version:

We thank the authors for addressing our previous comments. Overall, we are satisfied with the improvements and appreciate the hard work that has gone into this manuscript. We wish you all the best on the further publication pathway.

Reviewer #1 (Public Review):

This study presents valuable data on effector proteins (=virulence factors) used by the bacterial pathogen Legionella pneumophila that target host vesicle trafficking GTPases during infection. The evidence supporting the claims of the authors is robust, and the data suggest a sophisticated interplay between multiple effectors with the goal of temporarily exploiting host cell Rab10 during infection.

The authors have done a nice job addressing my earlier concerns. I have no further criticism about the revised paper.

Reviewer #2 (Public Review):

This manuscript explores the interplay between Legionella Dot/Icm effectors that modulate ubiquitination of the host GTPase Rab10. Rab10 undergoes phosphoribosyl-ubiquitination (PR-Ub) by the SidE family of effectors which is required for its recruitment to the Legionella containing vacuole (LCV). Through a series of elegant experiments using effector gene knockouts, co-transfection studies and careful biochemistry, Kubori et al further demonstrate that:

(1) The SidC family member SdcB contributes to the polyubiquitination (poly-Ub) of Rab10 and its retention at the LCV membrane.

(2) The transglutaminase effector, MavC acts as an inhibitor of SdcB by crosslinking ubiquitin at Gln41 to lysine residues in SdcB.

Reviewer #1 (Public Review):

Summary:

This study brings new information about the function of serotonin-gated ion channels 5-HT3AR, by describing the conformational changes undergoing during ligands binding. These results can be potentially extrapolated to other members of the Cys-loop ligand-gated ion channels. By combining fluorescence microscopy with electrophysiological recordings, the authors investigate structural changes inside and outside the orthosteric site elicited by agonists, partial agonists, and antagonists. The results are convincing and correlate well with the observations from cryo-EM structures. The work will be of important significance and broad interest to scientists working on channel biophysics but also drug development targeting ligand-gated ion channels.

Strengths:

The authors present an elegant and well-designed study to investigate the conformational changes on 5-HT3AR where they combine electrophysiological and fluorometry recordings. They determined four positions suitable to act as sensors for the conformational changes of the receptor: two inside and two outside the agonist binding site. They make a strong point showing how antagonists produce conformational changes inside the orthosteric site similarly as agonists do but they failed to spread to the lower part of the ECD, in agreement with previous studies and Cryo-EM structures. They also show how some loss-of-function mutant receptors elicit conformational changes (changes in fluorescence) after partial agonist binding but failed to produce measurable ionic currents, pointing to intermediate states that are stabilized in these conditions. The four fluorescence sensors developed in this study may be good tools for further studies on characterizing drugs targeting the 5-HT3R. The major conclusions of the manuscript seem well justified.

Weaknesses:

Weaknesses have been very well addressed during the review process.

Reviewer #2 (Public Review):

Summary:

This study focuses on the 5-HT3 serotonin receptor, a pentameric ligand-gated ion channel important in chemical neurotransmission. There are many cryo-EM structures of this receptor with diverse ligands bound, however assignment of functional states to the structures remains incomplete. The team applies voltage-clamp fluorometry to measure, at once, both changes in ion channel activity, and changes in fluorescence. Four cysteine mutants were selected for fluorophore labeling, two near the neurotransmitter site, one in the ECD vestibule, and one at the ECD-TMD junction. Agonists, partial agonists, and antagonists were all found to yield similar changes in fluorescence, a proxy for conformational change, near the neurotransmitter site. The strength of the agonist correlated to a degree with propagation of this fluorescence change beyond the local site of neurotransmitter binding. Antagonists failed to elicit a change in fluorescence in the vestibular of the ECD-TMD junction sites. The VCF results further turned up evidence supporting intermediate (likely pre-active) states.

Strengths:

The experiments appear rigorous, the problem the team tackles is timely and important, the writing and the figures are for the most part very clear. We sorely need approaches orthogonal to structural biology to annotate conformational states and observe conformational transitions in real membranes- this approach, and this study, get right to the heart of what is missing.

Weaknesses:

The weaknesses in the study itself are overall minor, I only suggest improvements geared toward clarity. What we are still missing is application of an approach like this to annotate the conformation of the part of the receptor buried in the membrane; there is an important debate about which structure represents which state, and that is not addressed in the current study.

Reviewer #3 (Public Review):

Summary:

The authors have examined the 5-HT3 receptor using voltage clamp fluorometry, which enables them to detect structural changes at the same time as the state of receptor activation. These are ensemble measurements, but they enable an impressive scheme of the action of different agonists and antagonists to be built up. The growing array of structural snapshots of 5-HT3 receptors is used to good effect to understand the results.

Strengths:

The combination of rigorously tested fluorescence reporters with oocyte electrophysiology across a large panel of ligands is a solid development for this receptor type.

Weaknesses:

In their revision, the authors corrected all the weaknesses of the original submission.

Joint Public Review:

Summary:

This important manuscript investigates a subpopulation of glutamatergic neurons in the suprammamillary nucleus that projects to the pre-optic hypothalamus area (SuM-VGLUT2+::POA). First, they define the neural circuitry of these neurons, which contact many stress/threat-associated brain regions. Then they employ fibre photometry to measure the activity of these neurons during various threatening tasks and find the responses correlate well with threat stimuli. Finally, they stimulate these neurons and find multiple lines of evidence that mice find this aversive and will act to avoid receiving this stimulation. In sum, they provide solid evidence that this neuronal population represents a new node in stress response circuitry that allows the animal to produce flexible behaviours in response to stress, which will be of interest to neuroscientists across several sub-fields.

Strengths:

Overall this is a solid manuscript tackling an important question. Coping with stress by an animal in danger is essential for survival. This manuscript identifies a novel population of neurons in the murine supramamillary nucleus (SuM) projecting to the pre-optic hypothalamus area among other regions that is involved in this important process. The evidence to support the conclusions is solid.

Specific strengths:

• The topic is novel.

• The manuscript follows a logical structure and neatly moves through the central story. Several potential alternate interpretations are well-controlled for.

• The manuscript employs an array of different tasks to provide converging evidence for their conclusions.

• The authors provide excellent evidence of the specificity of the function of this neuronal population, both from anatomical studies and from behavioural studies (e.g. demonstrating that activity of gabaergic neurons in the same region does not correlate with behaviours in the same way).

• The study is well-powered (sample sizes are good) and the effects are convincing.

Weaknesses:

* Not all of the reviewer comments were addressed in the manuscript itself, although this was acknowledged in the author's responses to reviewers. One key example is as follows:

* The authors did not entirely address comments related to rigor but they at least acknowledged it. For example, in multiple places they argue that WT, purchased mice are probably not different in baseline behavior compared to Vgltu2-IRES-Cre because it is unlikely that adding the IRES-Cre will change behavior. However, they do not acknowledge that transgenic lines are not from the exact same genetic background and generation number, and there is ample evidence in the literature that transgenic mice on a B6J background can differ in basal phenotypes from one another and B6J. In one place they show some basal behavior, at least in heat map form though not quantified. Had the authors decided to apply this more pervasively, it would have made the story even more compelling in terms of a stress/threat-induced phenotype.

Comments on revised version from the Reviewing Editor:

The authors have done a thorough job of answering the reviewer queries, and a good job of explaining why they have not answered a particular point. Indeed, there is so much additional information in response to the reviewers that I hope readers of the manuscript will read the reviews and responses as well! I think they add a lot.

Reviewer #3 (Public Review):

Summary:

Chang et al. investigated the mechanisms governing collagen fibrillogenesis, firstly demonstrating that cells within tail tendons are able to uptake exogenous collagen and use this to synthesize new collagen-1 fibrils. Using an endocytic inhibitor, the authors next showed that endocytosis was required for collagen fibrillogenesis and that this process occurs in a circadian rhythmic manner. Using knockdown and overexpression assays, it was then demonstrated that collagen fibril formation is controlled by vacuolar protein sorting 33b (VPS33b), and this VPS33b-dependent fibrillogenesis is mediated via Integrin alpha-11 (ITGA11). Finally, the authors demonstrated increased expression of VPS33b and ITGA11 at the gene level in fibroblasts from patients with idiopathic pulmonary fibrosis (IPF), and greater expression of these proteins in both lung samples from IPF patients and in chronic skin wounds, indicating that endocytic recycling is disrupted in fibrotic diseases.

Strengths:

The authors have performed a comprehensive functional analysis of the regulators of endocytic recycling of collagen, providing compelling evidence that VPS33b and ITGA11 are crucial regulators of this process.

Weaknesses:

Throughout the study, several different cell types have been used (immortalised tail tendon fibroblasts, NIHT3T cells, and HEK293T cells). In general, it is not clear which cells have been used for a particular experiment, and the rationale for using these different cell types is not explained. In addition, some experimental details are missing from the methods.

There is also a lack of functional studies in patient-derived IPF fibroblasts which means the link between endocytic recycling of collagen and the role of VPS33b and ITGA11 cannot be fully established.

Reviewer #2 (Public Review):

Summary:

In this manuscript, the authors describe a mechanism, by which fluorescently-labelled Collagen type I is taken up by cells via endocytosis and then incorporated into newly synthesized fibers via an ITGA11 and VPS33B-dependent mechanism. The authors claim the existence of this collagen recycling mechanism and link it to fibrotic diseases such as IPF and chronic wounds.

Strengths:

The manuscript is well-written, and experimentally contains a broad variation of assays to support their conclusions. Also, the authors added data of IPF patient-derived fibroblasts, patient-derived lung samples, and patient-derived samples of chronic wounds that highlight a potential in vivo disease correlation of their findings.

The authors were also analyzing the membrane topology of VPS33B and could unravel a likely 'hairpin' like conformation in the ER membrane.

Weaknesses:

Experimental evidence is missing that supports the non-degradative endocytosis of the labeled collagen.

The authors show and mention in the text that the endocytosis inhibitor Dyngo®4a shows an effect on collagen secretion. It is not clear to me how specific this readout is if the inhibitor affects more than endocytosis. This issue was unfortunately not further discussed. The authors use commercial rat tail collagen, it is unclear to me which state the collagen is in when it's endocytosed. Is it fully assembled as collagen fiber or are those single heterotrimers or homotrimers?

The Cy-labeled collagen is clearly incorporated into new fibers, but I'm not sure whether the collagen is needed to be endocytosed to be incorporated into the fibers or if that is happening in the extracellular space mediated by the cells.

In general for the collagen blots, due to the lack of molecular weight markers, what chain/form of collagen type I are you showing here?

Besides the VPS33B siRNA transfected cells the authors also use CRISPR/Cas9-generated KO. The KO cells do not seem to be a clean system, as there is still a lot of mRNA produced. Were the clones sequenced to verify the KO on a genomic level? For the siRNA transfection, a control blot for efficiency would be great to estimate the effect size. To me it is not clear where the endocytosed collagen and VPS33B eventually meet in the cells and whether they interact. Or is ITGA11 required to mediate this process, in case VPS33B is not reaching the lumen?

The authors show an upregulation of ITGA11 and VPS33B in IPF patients-derived fibroblasts, which can be correlated to an increased level of ColI uptake, however, it is not clear whether this increased uptake in those cells is due to the elevated levels of VPS33B and/or ITGA11.

Reviewer #1 (Public Review):

Summary:

The authors describe that the endocytic pathway is crucial for ColI fibrillogenesis. ColI is endocytosed by fibroblasts, prior to exocytosis and formation of fibrils, which can include a mixture of endogenous/nascent ColI chains and exogenous ColI. ColI uptake and fibrillogenesis are regulated by circadian rhythm as described by the authors in 2020, thanks to the dependence of this pathway on circadian-clock-regulated protein VPS33B. Cells are capable of forming fibrils with recently endocytosed ColI when nascent chains are not available. Previously identified VPS33B is demonstrated not to have a role in endocytosis of ColI, but to play a role in fibril formation, which the authors demonstrate by showing the loss of fibril formation in VPS33B KO, and an excess of insoluble fibrils - along-side a decrease in soluble ColI secretion - in VPS33B overexpression conditions. A VPS33B binding protein VIPAS39 is also shown to be required for fibrillogenesis and to colocalise with ColI. The authors thus conclude that ColI is internalised into endosomal structures within the cell, and that ColI, VPS33B, and VIPA39 are co-trafficked to the site of fibrillogenesis, where along with ITGA11, which by mass spectrometric analysis is shown to be regulated by VPS33B levels, ColI fibrils are formed. Interestingly, in involved human skin sections from idiopathic pulmonary fibrosis (IPF) patients, ITGA11 and VPS33B expression is increased compared to healthy tissue, while in patient-derived fibroblasts, uptake of fluorescently-labelled ColI is also increased. This suggests that there may be a significant contribution of endocytosis-dependent fibrillogenesis in the formation of fibrotic and chronic wound-healing diseases in humans.

Strengths:

This is an interesting paper that contributes an exciting novel understanding of the formation of fibrotic disease, which despite its high occurrence, still has no robust therapeutic options. The precise mechanisms of fibrillogenesis are also not well understood, so a study devoted to this complex and key mechanism is well appreciated. The dependence of fibrillogenesis on VPS33B and VIPA39 is convincing and robust, while the distinction between soluble ColI secretion and insoluble fibrillar ColI is interesting and informative.

Weaknesses:

There are a number of limitations to this study in its current state. Inhibition of ColI uptake is performed using Dyngo4a, which although proposed as an inhibitor of Clathrin-dependent endocytosis is known to be quite un-specific. This may not be a problem however, as the endocytic mechanism for ColI also does not seem to be well defined in the literature, in fact, the principle mechanism described in the papers referred to by the authors is that of phagocytosis. It would be interesting to explore this important part of the mechanism further, especially in relation to the intracellular destination of ColI. The circadian regulation does not appear as robust as the authors' last paper, however, there could be a larger lag between endocytosis of ColI and realisation of fibrils. The authors state that the endocytic pathway is the mechanism of trafficking and that they show ColI, VPS33B, and VIPA39 are co-trafficked. However, the only link that is put forward to the endosomes is rather tenuously through VPS33B/VIPA39. There is no direct demonstration of ColI localisation to endosomes (ie. immunofluorescence), and this is overstated throughout the text. Demonstrating the intracellular trafficking and localisation of ColI, and its actual relationship to VPS33B and VIPA39, followed by ITGA11, would broaden the relevance of this paper significantly to incorporate the field of protein trafficking. Finally, the "self-formation" of ColI fibrils is discussed in relation to the literature and the concentration of fluorescently-tagged ColI, however as the key message of the paper is the fibrillogenesis from exocytosed colI, I do not feel like it is demonstrated to leave no doubt. Specific inhibition of intracellular trafficking steps, or following the progressive formation of ColI fibrils over time by immunofluorescence would demonstrate without any further doubt that ColI must be endocytosed first, to form fibrils as a secondary step, rather than externally-added ColI being incorporated directly to fibrils, independent of cellular uptake.

Reviewer #1 (Public Review):

Summary:

The authors quantitatively describe the complex binding equilibria of BRAF and its inhibitors resulting in some cases in the paradoxical activation of BRAF dimer when bound to ATP competitive inhibitors. The authors use a biophysical tour de force involving FRET binding assays, NMR, kinase activity assays, and DEER spectroscopy.

Strengths:

The strengths of the study are the beautifully conducted assays that allow for a thorough characterization of the allostery in this complex system. Additionally, the use of F-NMR and DEER spectroscopy provides important insights into the details of the process.

The resulting model for binding of inhibitors and dimerization (Figure 4) is very helpful.

Weaknesses:

This is a complex system and its communication is inherently challenging. It might be of interest to the broader readership to understand the implications of the model for drug development and therapy.

Reviewer #2 (Public Review):

Summary:

This manuscript uses FRET, 19F-NMR, and DEER/EPR solution measurements to examine the allosteric effects of a panel of BRAF inhibitors (BRAFi). These include first-generation aC-out BRAFi, and more recent Type I and Type II aC-in inhibitors. Intermolecular FRET measurements quantify Kd for BRAF dimerization and inhibitor binding to the first and second subunits. Distinct patterns are found between aC-in BRAFi, where Type I BRAFi binds equally well to the first and second subunits within dimeric BRAF. In contrast, Type II BRAFi shows stronger affinity for the first subunit and weaker affinity for the second subunit, an effect named "allosteric asymmetry". Allosteric asymmetry has the potential for Type II inhibitors to promote dimerization while favoring occupancy of only one subunit (BBD form), leading to the enrichment of an active dimer.

Measurements of in vitro BRAF kinase activity correlate amazingly well with the calculated amounts of the half-site-inhibited BBD forms with Type II inhibitors. This suggests that the allosteric asymmetry mechanism explains paradoxical activation by this class of inhibitors. DEER/EPR measurements further examine the positioning of helix aC. They show systematic outward movement of aC with Type II inhibitors, relative to the aC-in state with Type I inhibitors, and further show that helix aC adopts multiple states and is therefore dynamic in apo BRAF. This makes a strong case that negative cooperativity between sites in the BRAF dimer can account for paradoxical kinase activation by Type II inhibitors by creating a half-site-occupied homodimer, BBD. In contrast, Type I inhibitors and aC-out inhibitors do not fit this model, and are therefore proposed to be explained by previously proposed models involving negative allostery between subunits in BRAF-CRAF heterodimers, RAS priming, and transactivation.

Strengths:

This study integrates orthogonal spectroscopic and kinetic strategies to characterize BRAF dynamics and determine how it impacts inhibitor allostery. The unique combination of approaches presented in this study represents a road map for future work in the important area of protein kinase dynamics. The work represents a worthy contribution not only to the field of BRAF regulation but to protein kinases in general.

Weaknesses:

Some questions remain regarding the proposed model for Type II inhibitors and its comparison to Type I and aC-out inhibitors that would be useful to clarify. Specifically, it would be helpful to address whether the activation of BRAF by Type II inhibitors, while strongly correlated with BBD model predictions in vitro, also depends on CRAF via BRAF-CRAF in cells and therefore overlaps with the mechanisms of paradoxical activation by Type I and aC-out inhibitors.

Reviewer #1 (Public Review):

Summary:

Szathmary and colleagues explore the parabolic growth regime of replicator evolution. Parabolic growth occurs when nucleic acid strain separation is the rate limiting step of the replication process which would have been the case for non-enzymatic replication of short oligonucleotide that could precede the emergence of ribozyme polymerases and helicases. The key result is that parabolic replication is conducive to the maintenance of genetic diversity, that is, coexistence of numerous master sequences (the Gause principle does not apply). Another important finding is that there is no error threshold for parabolic replication except for the extreme case of zero fidelity.

Strengths:

I find both the analytic and the numerical results to be quite convincing and well described. The results of this work are potentially important because they reveal aspects of a realistic evolutionary scenario for the origin of replicators.

Weaknesses:

There are no obvious technical weaknesses. It can be argued that the results represent an incremental advance because many aspects of parabolic replication have been explored previously (the relevant publications are properly cited). Obviously, the work is purely theoretical, experimental study of parabolic replication is due. In the opinion of this reviewer, though, these are understandable limitations that do not actually detract from the value of this work.

Reviewer #2 (Public Review):

Summary:

A dominant hypothesis concerning the origin of life is that, before the appearance of the first enzymes, RNA replicated non-enzymatically by templating. However, this replication was probably not very efficient, due to the propensity of single strands to bind to each other, thus inhibiting template replication. This phenomenon, known as product inhibition, has been shown to lead to parabolic growth instead of exponential growth. Previous works have shown that this situation limits competition between alternative replicators and therefore promotes RNA population diversity. The present work examines this scenario in an agent-based model of RNA replication, taking into account finite population size, mutations and differences in GC content. The main results are (1) confirmation that parabolic growth promotes diversity, but that when the population size is small enough, sequences least efficient at replicating may nevertheless go extinct; (2) the observation that fitness is not only controlled by the replicability of sequences, but also by their GC content ; (3) the observation that parabolic growth attenuates the impact of mutations and, in particular, that the error threshold to which exponentially growing sequences are subject can be exceeded, enabling sequence identity to be maintained at higher mutation rates.

Strengths:

The analyses are sound and the observations intriguing. Indeed, while it has been noted previously that parabolic growth promotes coexistence, this is the first work to show that it can also mitigate the error threshold catastrophe, which is often presented as a major obstacle to our understanding of the origin of life.

Weaknesses:

A general weakness, which can however be seen as inherent in an agent-based model that aims to be more realistic than earlier, more phenomenological models, is the proliferation of parameters. The choice and values of these parameters are generally justified and, in many cases, several values are tested to assess the robustness of the results, but it can be difficult for the reader to identify the modeling choices that are truly critical from those that are less so.

Reviewer #4 (Public Review):

Summary:

Masala N et al showed interesting aberrant calcium microwaves in the hippocampus when synapsin promoter driven GCaMPs were expressed for a long period of time. These aberrant hippocampal Ca2+ micro-waves depend on the viral titre of the GECI. The microwave of Ca2+ was not observed when GECI was expressed only a sparse set of neurons.

Strengths:

These findings are important to wide neuroscience community especially when considering a great number of investigators are using similar approaches. Results look convincing and are consistent across several laboratories.

Weaknesses:

Synapsin promoter labels both excitatory pyramidal neurons and inhibitory neurons. To avoid aberrant Ca2+ microwave, a combination of Flex virus and CaMKII-Cre or Thy-1-GCaMP6s and 6f mice were tested. However, all these approaches limit the number of infected pyramidal neurons. While the comprehensive display of these results is appreciated, one additional important test would be more informative. To distinguish whether the microwave of Ca2+ is sufficiently caused via the expression of GCaMP in interneurons, or just a matter of pyramidal neuron density, testing Flex-GCaMP6 in interneuron specific mouse lines such as PV-Cre and SOM-Cre will provide further clarifications.

Reviewer #2 (Public Review):

Summary:

The authors describe and quantify a phenomenon in the CA1 and CA3 of the hippocampus that they call aberrant Ca2+ micro-waves. Micro-waves are sometimes seen during 2-photon calcium imaging of populations of neurons under certain conditions. They are spatially confined slow calcium events that start in a few cells and slowly spread to neighboring groups of cells. This phenomenon has been uttered between researchers in the field at conferences, but no one has taken the time to carefully capture and quantify micro-waves and pin down the causes. The authors show that micro-waves are dependent on the viral titre of the genetically encoded calcium indicators (GECIs), the genetic promoter (synapsin), the neuronal subtype (granule cells in the dentate gyrus do not produce micro-waves and they are not seen in the neocortex), and the density of GECI expression. The authors should be commended for their work and for raising awareness to all labs doing any form of calcium imaging in populations of neurons. The authors also come up with alternative approaches to avoid artifactual micro-waves such as reducing the transduction titre (1:2 dilution of virus) and a transduction method employing sparser and cre-dependent GECI expression in principal cells using a CaMKII promoter.

Strengths:

The micro-waves reported in the paper were robustly observed across 4 laboratories and 3 different countries with various experimenters and calcium imaging set-ups. This adds significant strength to the work.

The age of mice used covered a broad range (from 6 to 43 weeks). This is a strength because it covers most ages that are used in labs that regularly do calcium imaging.

Another strength is they used different GCaMP variants (GCaMP6m, GCaMP6s, GCaMP7f), as well as a red indicator: RCaMP. This shows the micro-waves are not an issue with any particular GECI, as the authors suggest.

The authors include many movies of micro-waves. This is extremely useful for researchers in the field to view them in real-time so they can identify them in their own data.

They provide a useful table with specific details of the virus injected, titre, dilution, and other information along with the incidence of micro-waves. A nice look-up table for researchers to see if their viral strategy is associated with a high or low incidence of micro-waves.

Weaknesses:

The effect of mico-waves on single cell function was not analyzed. It would be useful, for example, if we knew the influence of micro-waves on place fields. Can a place cell still express a place field in a hippocampus that produces micro-waves? What effect might a microwave passing over a cell have on its place field? Mice were not trained in these experiments, so the authors do not have the data. However, they do briefly discuss these ideas.

Reviewer #3 (Public Review):

Summary:

The work by Masala and colleagues highlights a striking artifact that can result from a particular viral method for expressing genetically encoded calcium indicators (GECIs) in neurons. In a cross-institutional collaboration, the authors find that viral transduction of GECIs in the hippocampus can result in aberrant slow-traveling calcium (Ca2+) micro-waves. These Ca2+ micro-waves are distinct from previously described ictal activity but nevertheless are likely a pathological consequence of overexpression of virally transduced proteins. Ca2+ micro-waves will most-likely obscure the physiology that most researchers are interested in studying with GECIs, and their presence indicates that the neural circuit is in an unintended pathological state. Interestingly this pathology was not observed using the same viral transduction methods in other brain regions. The authors recommend several approaches that may help other experimenters avoid this confound in their own data such as reducing the titer of viral injections or using recombinase-dependent expression. The intent of this manuscript is to raise awareness of the potential unintended consequences of viral overexpression, particularly for GECIs. A rigorous investigation into the exact causes of Ca2+ micro-waves or the mechanisms supporting them are beyond the authors' intended scope.

Strengths:

The authors clearly demonstrate that Ca2+ micro-waves occur in the CA1 and CA3 regions of the hippocampus following large volume, high titer injections of adeno-associated viruses (AAV1 and AAV9) encoding GECIs. The supplementary videos provide undeniable proof of their existence.

By forming an inter-institutional collaboration, the authors demonstrate that this phenomenon is robust to changes in surgical techniques or imaging conditions.

Weaknesses:

I believe that the weaknesses of the manuscript are appropriately highlighted by the authors themselves in the discussion. The manuscript does not attempt to exhaustively characterize the conditions under which calcium micro-waves occur. Rather, the authors raise awareness of this problem.

Reviewer #1 (Public Review):

This paper identifies GABA cells in the preoptic hypothalamus and others in the posterior hypothalamus which are involved in REM sleep rebound (the increase in REM sleep) after selective REM sleep deprivation. By calcium photometry, these preoptic cells are most active during REM, and show more calcium signals during REM deprivation, suggesting they respond to "REM pressure". Inhibiting these cells ontogenetically diminishes REM sleep. The optogenetic and photometry work is carried out to a high standard, the paper is well written, and the findings are interesting and enhance our understanding of REM sleep regulation. The new findings make it clear that as for the circuitry that regulates NREM sleep, REM sleep circuitry is also quite distributed in the brain. It is unclear if there is a true "REM center". The study of mechanisms of catching up on lost sleep (sleep homeostasis), has previously focused on NREM sleep, where various circuits have been identified. That there is a special mechanism that also tracks time awake and compensates with REM sleep is intriguing.

In a broader context, the existence of REM rebound suggests that REM sleep must have a function, otherwise why catch up on it. There is a lot of literature that suggests REM contributes to emotional processing, for example. The new findings deepen our appreciation of REM regulation. As REM sleep is often disturbed in stress (e.g. post-traumatic stress disorder) and in depression, understanding more about REM regulation could ultimately aid treatments for people living with these conditions.

Reviewer #2 (Public Review):

Maurer et al investigated the contribution of GAD2+ neurons in the preoptic area (POA), projecting to the tuberomammillary nucleus (TMN), to REM sleep regulation. They applied an elegant design to monitor and manipulate activity of this specific group of neurons: a GAD2-Cre mouse, injected with retrograde AAV constructs in the TMN, thereby presumably only targeting GAD2+ cells projecting to the TMN. Using this set-up in combination with technically challenging techniques including EEG with photometry and REM sleep deprivation, the authors found that this cell-type studied becomes active shortly (≈40sec) prior to entering REM sleep and remains active during REM sleep. Moreover, optogenetic inhibition of GAD2+ cells inhibits REM sleep by a third, and also impairs the rebound in REM sleep in the following hour. Thus, the data makes a convincing case for a role of GAD2+ neurons in the POA projecting to the TMN in REM sleep regulation.

Reviewer #2 (Public Review):

Summary:

Roy et al. investigated the role of non-canonical DNA structures called G-quadruplexes (G4s) in long-range chromatin interactions and gene regulation. Introducing a G4 array into chromatin significantly increased the number of long-range interactions, both within the same chromosome (cis) and between different chromosomes (trans). G4s functioned as enhancer elements, recruiting p300 and boosting gene expression even 5 megabases away. The study proposes a mechanism where G4s directly influence 3D chromatin organization, facilitating communication between regulatory elements and genes.

Strength:

The findings are valuable for understanding the role of G4-DNA in 3D genome organization and gene transcription.

Weaknesses:

The study would benefit from more robust and comprehensive data, which would add depth and clarity.

(1) Lack of G4 Structure Confirmation: The absence of direct evidence for G4 formation within cells undermines the study's foundation. Relying solely on in vitro data and successful gene insertion is insufficient.

(2) Alternative Explanations: The study does not sufficiently address alternative explanations for the observed results. The inserted sequences may not form G4s or other factors like G4-RNA hybrids may be involved.

(3) Limited Data Depth and Clarity: ChIP-qPCR offers limited scope and considerable variation in some data makes conclusions difficult.

(4) Statistical Significance and Interpretation: The study could be more careful in evaluating the statistical significance and magnitude of the effects to avoid overinterpreting the results.

Reviewer #3 (Public Review):

Summary:

This paper aims to demonstrate the role of G-quadruplex DNA structures in the establishment of chromosome loops. The authors introduced an array of G4s spanning 275 bp, naturally found within a very well-characterized promoter region of the hTERT promoter, in an ectopic region devoid of G-quadruplex and annotated gene. As a negative control, they used a mutant version of the same sequence in which G4 folding is impaired. Due to the complexity of the region, 3 G4s on the same strand and one on the opposite strand, 12 point mutations were made simultaneously (G to T and C to A). Analysis of the 3D genome organization shows that the WT array establishes more contact within the TAD and throughout the genome than the control array. Additionally, a slight enrichment of H3K4me1 and p300, both enhancer markers, was observed locally near the insertion site. The authors tested whether the expression of genes located either nearby or up to 5 Mb away was up-regulated based on this observation. They found that four genes were up-regulated from 1.5 to 3-fold. An increased interaction between the G4 array compared to the mutant was confirmed by the 3C assay. For in-depth analysis of the long-range changes, they also performed Hi-C experiments and showed a genome-wide increase in interactions of the WT array versus the mutated form.

Strengths:

The experiments were well-executed and the results indicate a statistical difference between the G4 array inserted cell line and the mutated modified cell line.

Weaknesses:

The control non-G4 sequence contains 12 point mutations, making it difficult to draw clear conclusions. These mutations not only alter the formation of G4, but also affect at least three Sp1 binding sites that have been shown to be essential for the function of the hTERT promoter, from which the sequence is derived. The strong intermingling of G4 and Sp1 binding sites makes it impossible to determine whether all the observations made are dependent on G4 or Sp1 binding. As a control, the authors used Locked Nucleic Acid probes to prevent the formation of G4. As for mutations, these probes also interfere with two Sp1 binding sites. Therefore, using this alternative method has the same drawback as point mutations. This major issue should be discussed in the paper. It is also possible that other unidentified transcription factor binding sites are affected in the presented point mutants.

Reviewer #1 (Public Review):

Summary:

In this manuscript, Chowdhury and co-workers provide interesting data to support the role of G4-structures in promoting chromatin looping and long-range DNA interactions. The authors achieve this by artificially inserting a G4-containing sequence in an isolated region of the genome using CRISPR-Cas9 and comparing it to a control sequence that does not contain G4 structures. Based on the data provided, the authors can conclude that G4-insertion promotes long-range interactions (measured by Hi-C) and affects gene expression (measured by qPCR) as well as chromatin remodelling (measured by ChIP of specific histone markers).

Whilst the data presented is promising and partially supports the authors' conclusion, this reviewer feels that some key controls are missing to fully support the narrative. Specifically, validation of actual G4-formation in chromatin by ChIP-qPCR (at least) is essential to support the association between G4-formation and looping. Moreover, this study is limited to a genomic location and an individual G4-sequence used, so the findings reported cannot yet be considered to reflect a general mechanism/effect of G4-formation in chromatin looping.

Strengths:

This is the first attempt to connect genomics datasets of G4s and HiC with gene expression. The use of Cas9 to artificially insert a G4 is also very elegant.

Weaknesses:

Lack of controls, especially to validate G4-formation after insertion with Cas9. The work is limited to a single G4-sequence and a single G4-site, which limits the generalisation of the findings.

Reviewer #3 (Public Review):

Summary:

The manuscript uses a combination of evolutionary approaches and structural/dynamics observations to provide mechanistic insights in the adaptation of the Spike protein during the evolution of SARS-COV-2 variants. The conclusion that CAP sites should be taken in particular account when considering the impact of the emergence of new strains and mutations is warranted.

Strengths:

The results presented in this work are very well outlined with well-written text, pleasant and well-described pictures, didactical and clear description of the methods, e.g. the discussion of how the MD equilibration procedure is applied and evaluated is clear and well argument.<br /> The citation of relevant similar results with different approaches strengthens the reasoning; in particular, comparing the calculated scores with previous experimentally obtained data is one of the strongest points of the manuscript.

Weaknesses:

There are no replicas of the molecular dynamics (MD) simulations, understandable since it's not a MD-focused paper. However, the comparison of multiple replicas could enhance the reliability of the findings.

Reviewer #1 (Public Review):

The manuscript by Oleh et al. uses in vitro electrophysiology and compartmental modeling (via NEURON) to investigate the expression and function of HCN channels in mouse L2/3 pyramidal neurons. The authors conclude that L2/3 neurons have developmentally regulated HCN channels, the activation of which can be observed when subjected to large hyperpolarizations. They further conclude via blockade experiments that HCN channels in L2/3 neurons influence cellular excitability and pathway-specific EPSP kinetics, which can be neuromodulated. While the authors perform a wide range of slice physiology experiments, concrete evidence that L2/3 cells express functionally relevant HCN channels is limited. There are serious experimental design caveats and confounds that make drawing strong conclusions from the data difficult. Furthermore, the significance of the findings is generally unclear, given modest effect sizes and a lack of any functional relevance, either directly via in vivo experiments or indirectly via strong HCN-mediated changes in known operations/computations/functions of L2/3 neurons.

Specific points:

(1) The interpretability and impact of this manuscript are limited due to numerous methodological issues in experimental design, data collection, and analysis. The authors have not followed best practices in the field, and as such, much of the data is ambiguous and/or weak and does not support their interpretations (detailed below). Additionally, the authors fail to appropriately explain their rationale for many of their choices, making it difficult to understand why they did what they did. Furthermore, many important references appear to be missing, both in terms of contextualizing the work and in terms of approach/method. For example, the authors do not cite Kalmbach et al 2018, which performed a directly comparable set of experiments on HCN channels in L2/3 neurons of both humans and mice. This is an unacceptable omission. Additionally, the authors fail to cite prior literature regarding the specificity or lack thereof of Cs+ in blocking HCN. In describing a result, the authors state "In line with previous reports, we found that L2/3 PCs exhibited an unremarkable amount of sag at 'typical' current commands" but they then fail to cite the previous reports.

(2) A critical experimental concern in the manuscript is the reliance on cesium, a nonspecific blocker, to evaluate HCN channel function. Cesium blocks HCN channels but also acts at potassium channels (and possibly other channels as well). The authors do not acknowledge this or attempt to justify their use of Cs+ and do not cite prior work on this subject. They do not show control experiments demonstrating that the application of Cs+ in their preparation only affects Ih. Additionally, the authors write 1 mM cesium in the text but appear to use 2 mM in the figures. In later experiments, the authors switch to ZD7288, a more commonly used and generally accepted more specific blocker of HCN channels. However, they use a very high concentration, which is also known to produce off-target effects (see Chevaleyre and Castillo, 2002). To make robust conclusions, the authors should have used both blockers (at accepted/conservative concentrations) for all (or at least most) experiments. Using one blocker for some experiments and then another for different experiments is fraught with potential confounds.

(3) A stronger case could be made that HCN is expressed in the somatic compartment of L2/3 cells if the authors had directly measured HCN-isolated currents with outside-out or nucleated patch recording (with appropriate leak subtraction and pharmacology). Whole-cell voltage-clamp in neurons with axons and/or dendrites does not work. It has been shown to produce erroneous results over and over again in the field due to well-known space clamp problems (see Rall, Spruston, Williams, etc.). The authors could have also included negative controls, such as recordings in neurons that do not express HCN or in HCN-knockout animals. Without these experiments, the authors draw a false equivalency between the effects of cesium and HCN channels, when the outcomes they describe could be driven simply by multiple other cesium-sensitive currents. Distortions are common in these preparations when attempting to study channels (see Williams and Womzy, J Neuro, 2011). In Fig 2h, cesium-sensitive currents look too large and fast to be from HCN currents alone given what the authors have shown in their earlier current clamp data. Furthermore, serious errors in leak subtraction appear to be visible in Supplementary Figure 1c. To claim that these conductances are solely from HCN may be misleading.

(4) The authors present current-clamp traces with some sag, a primary indicator of HCN conductance, in Figure 2. However, they do not show example traces with cesium or ZD7288 blockade. Additionally, the normalization of current injected by cellular capacitance and the lack of reporting of input resistance or estimated cellular size makes it difficult to determine how much current is actually needed to observe the sag, which is important for assessing the functional relevance of these channels. The sag ratio in controls also varies significantly without explanation (Figure 6 vs Figure 7). Could this variability be a result of genetically defined subgroups within L2/3? For example, in humans, HCN expression in L2/3 varies from superficial and deep neurons. The authors do not make an effort to investigate this. Regardless of inconsistencies in either current injection or cell type, the sag ratio appears to be rather modest and similar to what has already been reported previously in other papers.

(5) In the later experiments with ZD7288, the authors measured EPSP half-width at greater distances from the soma. However, they use minimal stimulation to evoke EPSPs at increasingly far distances from the soma. Without controlling for amplitude, the authors cannot easily distinguish between attenuation and spread from dendritic filtering and additional activation and spread from HCN blockade. At a minimum, the authors should share the variability of EPSP amplitude versus the change in EPSP half-width and/or stimulation amplitudes by distance. In general, this kind of experiment yields much clearer results if a more precise local activation of synapses is used, such as dendritic current injection, glutamate uncaging, sucrose puff, or glutamate iontophoresis. There are recording quality concerns here as well: the cell pictured in Figure 3a does not have visible dendritic spines, and a substantial amount of membrane is visible in the recording pipette. These concerns also apply to the similar developmental experiment in 6f-h, where EPSP amplitude is not controlled, and therefore, attenuation and spread by distance cannot be effectively measured. The outcome, that L2/3 cells have dendritic properties that violate cable theory, seems implausible and is more likely a result of variable amplitude by proximity.

(6) Minimal stimulation used for experiments in Figures 3d-i and Figures 4g-h does not resolve the half-width measurement's sensitivity to dendritic filtering, nor does cesium blockade preclude only HCN channel involvement. Example traces should be shown for all conditions in 3h; the example traces shown here do not appear to even be from the same cell. These experiments should be paired (with and without cesium/ZD). The same problem appears in Figure 4, where it is not clear that the authors performed controls and drug conditions on the same cells. 4g also lacks a scale bar, so readers cannot determine how much these measurements are affected by filtering and evoked amplitude variability. Finally, if we are to believe that minimal stimulation is used to evoke responses of single axons with 50% fail rates, NMDA receptor activation should be minimal to begin with. If the authors wish to make this claim, they need to do more precise activation of NMDA-mediated EPSPs and examine the effects of ZD7288 on these responses in the same cell. As the data is presented, it is not possible to draw the conclusion that HCN boosts NMDA-mediated responses in L2/3 neurons.

(7) The quality of recordings included in the dataset has concerning variability: for example, resting membrane potentials vary by >15-20 mV and the AP threshold varies by 20 mV in controls. This is indicative of either a very wide range of genetically distinct cell types that the authors are ignoring or the inclusion of cells that are either unhealthy or have bad seals.

(8) The authors make no mention of blocking GABAergic signaling, so it must be assumed that it is intact for all experiments. Electrical stimulation can therefore evoke a mixture of excitatory and inhibitory responses, which may well synapse at very different locations, adding to interpretability and variability concerns.

(9) The investigation of serotonergic interaction with HCN channels produces modest effect sizes and suffers the same problems as described above.

(10) The computational modeling is not well described and is not biologically plausible. Persistent and transient K channels are missing. Values for other parameters are not listed. The model does not seem to follow cable theory, which, as described above, is not only implausible but is also not supported by the experimental findings.

Taken together, there are serious methodological and analytical concerns that need to be addressed before the authors' claims can be supported. Combined with the small effect sizes and high data variability throughout the paper, this makes it hard to see how the manuscript could make a strong contribution to advancing our understanding of L2/3 cortical pyramidal neuron function.

Reviewer #2 (Public Review):

Summary:

This paper by Olah et al. uncovers a previously unknown role of HCN channels in shaping synaptic inputs to L2/3 cortical neurons. The authors demonstrate using slice electrophysiology and computational modeling that, unlike layer 5 pyramidal neurons, L2/3 neurons have an enrichment of HCN channels in the proximal dendrites. This location provides a locus of neuromodulation for inputs onto the proximal dendrites from L4 without an influence on distal inputs from L1. The authors use pharmacology to demonstrate the effect of HCN channels on NMDA-mediated synaptic inputs from L4. The authors further demonstrate the developmental time course of HCN function in L2/3 pyramidal neurons. Taken together, this a well-constructed investigation of HCN channel function and the consequences of these channels on synaptic integration in L2/3 pyramidal neurons.

Strengths:

The authors use careful, well-constrained experiments using multiple pharmacological agents to asses HCN channel contributions to synaptic integrations. The authors also use a voltage clamp to directly measure the current through HCN channels across developmental ages. The authors also provide supplemental data showing that their observation is consistent across multiple areas of the cerebral cortex.

Weaknesses: