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 #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 #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 #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 #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 #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 #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 #1 (Public Review):

The authors investigated how global brain activity varied during reward-based motor learning. During early learning, they found increased covariance between the sensorimotor and dorsal attention networks, coupled with reduced covariance between the sensorimotor and default mode networks; during late learning, they found the opposite pattern. Individual learning performance varied only with changes in the dorsal attention network. The authors certainly used a wide variety of valuable, state-of-the-art techniques to interrogate whole-brain networks and extract the key components of learning behavior. However, the findings are incomplete, tempered by potential confounds in the experimental design. As such, the underlying claim regarding how these networks jointly support reward-based motor learning is unclear.

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 #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 #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 #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 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 #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 #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.

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.

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 #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 #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 #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 #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 #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 #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 #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 #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 #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:

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 #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 #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 #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 #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 #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 #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):

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:

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 #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 #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 #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 #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 #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 #1 (Public Review):

This important study uses a wide variety of convincing, state-of-the-art neuroimaging analyses to characterize whole-brain networks and relate them to reward-based motor learning. During early learning, the authors found increased covariance between the sensorimotor and dorsal attention networks, coupled with reduced covariance between the sensorimotor and default mode networks. During late learning, they observed the opposite pattern. It remains to be seen whether these changes reflect generic changes in task engagement during learning or are specific to reward-based motor learning. This study is highly relevant for researchers interested in reward-based motor learning and decision-making.

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 #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 #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 #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 #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):

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 #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 #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 #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 #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 #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 #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 #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 #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 #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 #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.

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 #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 #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 #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):

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 #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 #1 (Public Review):

Summary:

Organization of cell surface receptors in membrane nanodomains is important for signaling, but how this is regulated is poorly understood. In this study the authors employ TIRFM single-molecule tracking combined with multiple analyses to show that ligand exposure increases diffusion of the immune receptor FLS2 in the plasma membrane and its co-localization with remorin REM1.3 in a manner dependent on the phosphosite S938. They additionally show that ligand increases dwell time of FLS2, and this is linked to FLS2 endocytosis, also in a manner dependent on S938 phosphorylation. The study uncovers a regulatory mechanism of FLS2 localization in the nanodomain crucial for signaling.

Strengths:

TIRFM single-molecule tracking, FRAP, FRET and endocytosis experiments were nicely done. A role of S938 phosphorylation is convincing.

Weaknesses:

In the previous submission, reviewers pointed out multiple issues, which the reviewers believed the authors can address in the revision. The revised version does improve to some extent but still contains many issues in terms of data analysis and writing.

Reviewer #1 (Public Review):

Summary:

In this paper, Misic et al showed that white matter properties can be used to classify subacute back pain patients that will develop persisting pain.

Strengths:

Compared to most previous papers studying associations between white matter properties and chronic pain, the strength of the method is to perform a prediction in unseen data. Another strength of the paper is the use of three different cohorts. This is an interesting paper that provides a valuable contribution to the field.

Weaknesses:

The authors imply that their biomarker could outperform traditional questionnaires to predict pain: "While these models are of great value showing that few of these variables (e.g. work factors) might have significant prognostic power on the long-term outcome of back pain and provide easy-to-use brief questionnaires-based tools, (21, 25) parameters often explain no more than 30% of the variance (28-30) and their prognostic accuracy is limited.(31)". I don't think this is correct; questionnaire-based tools can actually achieve far greater prediction than their model in about half a million individuals from the UK Biobank (Tanguay-Sabourin et al., A prognostic risk score for the development and spread of chronic pain, Nature Medicine 2023).

Moreover, the main weakness of this study is the sample size. It remains small despite having 3 cohorts. This is problematic because results are often overfitted in such a small sample size brain imaging study, especially when all the data are available to the authors at the time of training the model (Poldrack et al., Scanning the horizon: towards transparent and reproducible neuroimaging research, Nature Reviews in Neuroscience 2017). Thus, having access to all the data, the authors have a high degree of flexibility in data analysis, as they can retrain their model any number of times until it generalizes across all three cohorts. In this case, the testing set could easily become part of the training making it difficult to assess the real performance, especially for small sample size studies.

Even if the performance was properly assessed, their models show AUCs between 0.65-0.70, which is usually considered as poor, and most likely without potential clinical use. Despite this, their conclusion was: "This biomarker is easy to obtain (~10 min 18 of scanning time) and opens the door for translation into clinical practice." One may ask who is really willing to use an MRI signature with a relatively poor performance that can be outperformed by self-report questionnaires?

Overall, these criticisms are more about the wording sometimes used and the inference they made. I think the strength of the evidence is incomplete to support the main claims of the paper.

Despite these limitations, I still think this is a very relevant contribution to the field. Showing predictive performance through cross-validation and testing in multiple cohorts is not an easy task and this is a strong effort by the team. I strongly believe this approach is the right one and I believe the authors did a good job.

Minor points:

Methods:

I get the voxel-wise analysis, but I don't understand the methods for the structural connectivity analysis between the 88 ROIs. Have the authors run tractography or have they used a predetermined streamlined form of 'population-based connectome'? They report that models of AUC above 0.75 were considered and tested in the Chicago dataset, but we have no information about what the model actually learned (although this can be tricky for decision tree algorithms).

Minor:<br /> What results are shown in Figure 7? It looks more descriptive than the actual results.

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 #1 (Public Review):

Galanti et al. present an innovative new method to determine the susceptibility of large collections of plant accessions towards infestations by herbivores and pathogens. This work resulted from an unplanned infestation of plants in a greenhouse that was later harvested for sequencing. When these plants were extracted for DNA, associated pest DNA was extracted and sequenced as well. In a standard analysis, all sequencing reads would be mapped to the plant reference genome and unmapped reads, most likely originating from 'exogenous' pest DNA, would be discarded. Here, the authors argue that these unmapped reads contain valuable information and can be used to quantify plant infestation loads.

For the present manuscript, the authors re-analysed a published dataset of 207 sequenced accessions of Thlaspi arvense. In this data, 0.5% of all reads had been classified as exogenous reads, while 99.5% mapped to the T. arvense reference genome. In a first step, however, the authors repeated read mapping against other reference genomes of potential pest species and found that a substantial fraction of 'ambiguous' reads mapped to at least one such species. Removing these reads improved the results of downstream GWAs, and is in itself an interesting tool that should be adopted more widely.

The exogenous reads were primarily mapped to the genomes of the aphid Myzus persicae and the powdery mildew Erysiphe cruciferarum, from which the authors concluded that these were the likely pests present in their greenhouse. The authors then used these mapped pest read counts as an approximate measure of infestation load and performed GWA studies to identify plant gene regions across the T. arvense accessions that were associated with higher or lower pest read counts. In principle, this is an exciting approach that extracts useful information from 'junk' reads that are usually discarded. The results seem to support the authors' arguments, with relatively high heritabilities of pest read counts among T. arvense accessions, and GWA peaks close to known defence genes. Nonetheless, I do feel that more validation would be needed to support these conclusions, and given the radical novelty of this approach, additional experiments should be performed.

A weakness of this study is that no actual aphid or mildew infestations of plants were recorded by the authors. They only mention that they anecdotally observed differences in infestations among accessions. As systematic quantification is no longer possible in retrospect, a smaller experiment could be performed in which a few accessions are infested with different quantities of aphids and/or mildew, followed by sequencing and pest read mapping. Such an approach would have the added benefit of allowing causally linking pest read count and pest load, thereby going beyond correlational associations.

On a technical note, it seems feasible that mildew-infested leaves would have been selected for extraction, but it is harder to explain how aphid DNA would have been extracted alongside plant DNA. Presumably, all leaves would have been cleaned of live aphids before they were placed in extraction tubes. What then is the origin of aphid DNA in these samples? Are these trace amounts from aphid saliva and faeces/honeydew that were left on the leaves? If this is the case, I would expect there to be substantially more mildew DNA than aphid DNA, yet the absolute read counts for aphids are actually higher. Presumably read counts should only be used as a relative metric within a pest organism, but this unexpected result nonetheless raises questions about what these read counts reflect. Again, having experimental data from different aphid densities would make these results more convincing.

Reviewer #1 (Public Review):

Summary:

The authors provide solid evidence with a mouse model as well as supporting in vitro and analysis of clinical samples that loss of Fak increases the development of BRAF V600E-induced dysplastic lesions and carcinomas in the cecum via downregulation of Egfr-mediated Erk phosphorylation. This fine-tuning of Erk phosphorylation increases the expression of Lrg4 mRNA expression and promotes Lrg4 stability through downregulation of the E3 ubiquitin ligase Nedd4. The high Lrg4 expression correlates with an increased intestinal stem cell transcriptional signature that the authors suggest drives higher rates of transformation. This provides important insight that factors such as FAK may be able to modulate MAPK-driven tumorigenesis in specific circumstances. The data presented here are largely specific to the cecum. While these specific findings may ultimately have practical implications for human CRC outside the cecum and even therapeutic implications, these remain unexplored and will be a point for future investigations.

Strengths:

The authors use a mouse model (intestinal specific BRAF V600E +/- Fak knockout) as well as supporting in vitro analyses and clinical sample characterization to support their model. For both in vitro and in vivo studies, the authors use a combination of genetic and pharmacologic (including EGFR, FAK, and MEK inhibitors) tools to modulate the MAPK pathway. They also use a combination of transcriptional (RNA-Seq) and protein (IHC and Western blotting) readouts to support their proposed model. Importantly, they use a distinct mouse model (mutant Kras) to demonstrate their findings with Fak loss are specific to instances where EGFR can modulate ERK activation, providing strong evidence for their model. Finally, they also correlate their findings in the murine model with patient samples and with trends in the TCGA database. Collectively, these create a solid and convincing basis for their proposed model.

Weaknesses:

(1) The murine data is largely confined to the cecum. While the analysis of the cecum is appropriate based on the cecum specificity of their phenotype, they often use these findings to make broader generalizations about the nature of tumorigenesis in the intestinal epithelia and in CRC more generally. In my opinion, there was insufficient evidence presented supporting the extension of the proposed model beyond the cecum. While this is a weakness, it could be part of a growing effort to characterize left and right-sided malignancies as related but separate disease processes.

(2) The authors generally do a good job of focusing their analysis on the cecum and supporting their model. For example, Figure 5A examines different colon compartments, including the cecum. However, the authors fail to demonstrate that Fak loss only promotes Lrg4 upregulation in the cecum, where they observe an increase in BRAF V600E dysplasia and carcinoma. This is again seen in Figure 6A, where they only characterize Nedd4 expression in the cecum and not other compartments of the colon.

(3) The authors evaluate a broad range of tissues, including normal colonic mucosa, polyps, pre-cancerous dysplastic lesions, adenocarcinomas, and adenocarcinoma cell lines. While this breadth is a strength of the paper, the authors, at times, equate experimental observations in each of these conditions, despite the difference in the biology of these tissues/cells. For example, in their mouse model, they equate the development of dysplastic lesions and carcinoma lesions. This makes it difficult to accurately interpret their data and conclusions.

(4) In Figure 5i, this experiment was only completed in one cell line (HT29), despite the conclusion that Lrg4 expression is increased by decreased ERK phosphorylation due to protein stabilization. HT29 cells are a transformed human CRC cell line, quite different than a pre-malignant cecum intestinal epithelial cell. While convincing, the authors could have performed this key experiment in non-transformed murine cecal organoids (as they did for other experiments in Figure 5E), which would better recapitulate the mouse and pre-malignant setting to explain their mouse phenotype.

(5) While a large portion of the discussion focusses on the therapeutic implications of these findings, the authors only really investigate tumorigenesis. They likely have additional investigations planned for future manuscripts.

Reviewer #1 (Public Review):

Summary:

The authors establish a recombinant insect cell expression and purification scheme for the antiviral Dicer complex of C. elegans. In addition to Dicer-1, the complex harbors two additional proteins, the RIG-I-like helicase DRH-1 and the dsRNA-binding protein RDE-4. The authors show that the complex prefers blunt-end dsRNA over dsRNAs that contain overhangs. Furthermore, whereas ATP-dependent dsRNA cleavage only exacerbates regular dsRNA cleavage activity, the presence of RDE-4 is essential to ATP-dependent and ATP-independent dsRNA cleavage. Single-particle cryo-EM studies of the ternary C. elegans Dicer complex reveal that the N-terminal domain of DRH-1 interacts with the helicase domain of DCR-1, thereby relieving its autoinhibitory state. Last, the authors show that the ternary complex is able to processively cleave long dsRNA, an activity primarily relying on the helicase activity of DRH-1.

Strengths:

• First thorough biochemical characterization of the antiviral activity of C. elegans Dicer in complex with the RIG-I like helicase DRH-1 and the dsRNA-binding protein RDE-4<br /> • Discovery that RDE-4 is essential to dsRNA processing, whereas ATP hydrolysis is not<br /> • Discovery of an autoinhibitory role of DRH-1's N-terminal domain (in analogy to the CARD domains of RIG-I)<br /> • First structural insights into the ternary complex DCR-1:DRH-1:RDE-4 by cryo-EM to medium resolution<br /> • Trap experiments reveal that the ternary DCR-1 complex cleaves blunt-ended dsRNA processively. Likely, the helicase domain of DRH-1 is responsible for this processive cleavage.

Weaknesses:

• Cryo-EM Structure of the ternary Dicer-1:DRH-1:RED-4 complex to only medium resolution<br /> • High-resolution structure of the C-terminal domain of DRH-1 bound to dsRNA does not reveal the mechanism of how blunt-end dsRNA and overhang-containing one are being discriminated<br /> • The cryo-EM structure of DCR1:DRH-1:RDE-4 in the presence of ATP only reveals the helicase and CTD domains of DRH-1 bound to dsRNA. No information on dsRNA termini recognition is presented. The paragraph seems detached from the general flow of the manuscript.

Reviewer #1 (Public Review):

Summary:

In this work, Odenwald and colleagues show that mutant biotin ligases used to perform proximity-dependent biotin identification (TurboID) can be used to amplify signal in fluorescence microscopy and to label phase-separated compartments that are refractory to many immunofluorescence approaches. Using the parasite Trypanosoma brucei, they show that fluorescent methods such as expansion microscopy and CLEM, which require bright signals for optimal detection, benefit from the elevated signal provided by TurboID fusion proteins when coupled with labeled streptavidin. Moreover, they show that phase-separated compartments, where many antibody epitopes are occluded due to limited diffusion and potential sequestration, are labeled reliably with biotin deposited by a TurboID fusion protein that localizes within the compartment. They show successful labeling of the nucleolus, likely phase-separated portions of the nuclear pore, and stress granules. Lastly, they use a panel of nuclear pore-TurboID fusion proteins to map the regions of the T. brucei nuclear pore that appear to be phase-separated by comparing antibody labeling of the protein, which is susceptible to blocking, to the degree of biotin deposition detected by streptavidin, which is not.

Strengths:

Overall, this study shows that TurboID labelling and fluorescent streptavidin can be used to boost signal compared to conventional immunofluorescence in a manner similar to tyramide amplification, but without having to use antibodies. TurboID could prove to be a viable general strategy for labeling phase-separated structures in cells, and perhaps as a means of identifying these structures, which could also be useful.

Weaknesses:

However, I think that this work would benefit from additional controls to address if the improved detection that is being observed is due to the increased affinity and smaller size of streptavidin/biotin compared to IgGs, or if it has to do with the increased amount of binding epitope (biotin) being deposited compared to the number of available antibody epitopes. I also think that using the biotinylation signal produced by the TurboID fusion to track the location of the fusion protein and/or binding partners in cells comes with significant caveats that are not well addressed here, mostly due to the inability to discern which proteins are contributing to the observed biotin signal.

To dissect the contributions of the TurboID fusion to elevating signal, anti-biotin antibodies could be used to determine if the abundance of the biotin being deposited by the TurboID is what is increasing detection, or if streptavidin is essential for this. Alternatively, HaloTag or CLIP tagging could be used to see if diffusion of a small molecule tag other than biotin can overcome the labeling issue in phase-separated compartments. There are Halo-biotin substrates available that would allow the conjugation of 1 biotin per fusion protein, which would allow the authors to dissect the relative contributions of the high affinity of streptavidin from the increased amount of biotin that the TurboID introduces.

The idea of using the biotin signal from the TurboID fusion as a means to track the changing localization of the fusion protein or the location of interacting partners is an attractive idea, but the lack of certainty about what proteins are carrying the biotin signal makes it very difficult to make clear statements. For example, in the case of TurboID-PABP2, the appearance of a biotin signal at the cell posterior is proposed to be ALPH1, part of the mRNA decapping complex. However, because we are tracking biotin localization and biotin is being deposited on a variety of proteins, it is not formally possible to say that the posterior signal is ALPH1 or any other part of the decapping complex. For example, the posterior labeling could represent a localization of PABP2 that is not seen without the additional signal intensity provided by the TurboID fusion. There are also many cytoskeletal components present at the cell posterior that could be being biotinylated, not just the decapping complex. Similar arguments can be made for the localization data pertaining to MLP2 and NUP65/75. I would argue that the TurboID labeling allows you to enhance signal on structures, such as the NUPs, and effectively label compartments, but you lack the capacity to know precisely which proteins are being labeled.

Reviewer #1 (Public Review):

The manuscript by Wang et al. investigates the role of Rnf220 in hindbrain development and Hox expression. The authors suggest that Rnf220 controls Hox expression in the hindbrain by regulating WDR5 levels. The authors combine in vivo experiments with experiments in P19 cells to demonstrate this mechanism. However, the in vivo data does not provide strong support for the claims the authors make and the role of Rnf in Hox maintenance and pons development is unclear.

Specific concerns with in vivo data:

A major issue throughout the paper is that Hox expression analysis is done exclusively through quantitative PCR, with values ranging from 2-fold to several thousand-fold upregulation, with no antibody validation for any Hox protein (presumably they are all upregulated).

In Figure 1, massive upregulation of most Hox genes in the brainstem is shown after e16.5 but the paper quickly focuses on analysis of PN nuclei. What are the other consequences of this broad upregulation of Hox genes in the brainstem? There is no discussion of the overall phenotype of the mice, the structure of the brainstem, the migration of neurons, etc. The very narrow focus on motor cortex projections to PN nuclei seems bizarre without broad characterization of the mice, and the brainstem in particular. There is only a mention of "severe motor deficits" from previous studies, but given the broad expression of Rnf220, the fact that is a global knockout, and the effects on spinal cord populations shown previously the justification for focusing on PN nuclei does not seem strong.

It is stated that cluster 7 in scRNA-seq corresponds to the PN nuclei. The modest effect shown on Hox3-5 expression in that data in Figure 1 is inconsistent with the larger effect shown in Figure 2.

Presumably, Hox genes are not the only targets of Rnf220 as shown in the microarray/RNA-sequencing data. There is no definitive evidence that any phenotypes observed (which are also not clear) are specifically due to Hox upregulation. The only assay the authors use to look at a Hox-dependent phenotype in the brainstem is the targeting of PN nuclei by motor cortex axons. This is only done in 2 animals and there are no details as to how the data was analyzed and quantified. The only 2 images shown are not convincing of a strong phenotype, they could be taken at slightly different levels or angles. At the very least, serial sections should be shown and the experiment repeated in more animals. There is also no discussion of how these phenotypes, if real, would relate to previous work by the Rijli group which showed very precise mechanisms of synaptic specificity in this system.

The temporal aspect of this regulation in vivo is not clear. The authors show some expression changes begin at e16.5 but are also present at 2 months. Is the presumed effect on neural circuits a result of developmental upregulation at late embryonic stages or does the continuous overexpression in adult mice have additional influence? Are any of the Hox genes upregulated normally expressed in the brainstem, or PN specifically, at 2 months? Why perform single-cell sequencing experiments at 2 months if this is thought to be mostly a developmental effect? Similarly, the significance of the upregulated WRD5 in the pons and pontine nuclei at 2 months in Figure 3 is not clear.

In Figure 3C the levels of RNF220 in wt and het don't seem to be that different.

Based on the single-cell experiments, and the PN nuclei focus, the rescue experiments are confusing. If the Rnf220 deletion has a sustained effect for up to 2 months, why do the injections in utero? If the focus is the PN nuclei why look at Hox9 expression and not Hox3-5 which are the only Hox genes upregulated in PN based on sc-sequencing? No rescue of behavior or any phenotype other than Hox expression by qPCR is shown and it is unclear whether upregulation of Hox9 paralogs leads to any defects in the first place. The switch to the Nes-cre driver is not explained. Also, it seems that wdr5 mRNA levels are not so relevant and protein levels should be shown instead (same for rescue experiments in P19 cells).

Other:<br /> What is the relationship between Retinoic acid and WRD5? In Figure 3E there is no change in WRD5 levels without RA treatment in Rnf KO but an increase in expression with RA treatment and Rnf KO. However, the levels of WRD5 do not seem to change with RA treatment alone. Does Rnf220 only mediate WDR5 degradation in the presence of RA? This does not seem to be the case in experiments in 293 cells in Figure 4.

Why are the levels of Hox upregulation after RA treatment so different in Figure 5 and Figure Supplement 5?

In Figures 4B+C which lanes are input and which are IP? There is no quantitation of Figure 4D, from the blot it does look that there is a reduction in the last 2 columns as well. The band in the WT flag lane seems to have a bubble. Need to quantitate band intensities. Same for E, the effect does not seem to be completely reversed with MG132.

Reviewer #1 (Public Review):

Summary:

The authors were trying to achieve that Tgif1 expression is regulated by EAK1/2 and PTH in a time-dependent manner, and its roles in suppressing Pak3 for facilitating osteoblast adhesion. The authors further tried to achieve that the Tgif1-Pak3 signaling plays a significant role in osteoblast migration to the site of bone repair and bone remodeling.

Strengths:

- In a previous study, they demonstrated that Tgif1 is a target gene of PTH, and the absence of Tgif1 failed to increase bone mass by PTH treatment (Saito et al., Nat Commun., 2019). In this study, they found that Tgif1-Pak3 signaling prompts osteoblast migration through osteoblast adhesion to prompt bone regeneration. This novel finding provides a better understanding of how Tgif1 expression in osteoblasts regulates adherence, spreading, and migration during bone healing and bone remodeling.<br /> - The authors demonstrated that ERK1/2 and PTH regulate Tgif1 expression in a time-dependent manner and its role in suppressing Pak3 through various experimental approaches such as luciferase assay, ChIP assay, and gene silencing. These results contribute to the overall strength of the article.

Weaknesses:

None after substantial revisions especially in vivo parts.

Reviewer #2 (Public Review):

Summary:

in this paper authors show that the degree of pigmentation for RPE cells is not correlated with a level of maturation and function. They suggest that this status could be different in vitro than in vivo but do not provide proper experiments to validate this hypothesis. However, it is the first time that the absence of correlation between pigmentation and function is studied.

Strengths:

The methods are good and experiments very rigorous

Comments on current version:

The authors have modified their title and focus on QC for in vitro process

Reviewer #1 (Public Review):

Summary:

A novel serine protease and a inhibitor pair regulate cell migration in the neural crest.

Strengths:

The reproduction of classical cranial neural crest extirpations and their phenocopy by SerpinE2 morpholino are remarkable. Very scholarly written and data of the highest quality.

Weaknesses:

All were improved upon revision.

Reviewer #1 (Public Review):

Summary:

Building upon their famous tool for the deconvolution of human transcriptomics data (EPIC), Gabriel et al. implemented a new methodology for the quantification of the cellular composition of samples profiled with Assay for Transposase-Accessible Chromatin sequencing (ATAC-Seq). To build a signature for ATAC-seq deconvolution, they first created a compendium of ATAC-seq data and derived chromatin accessibility marker peaks and reference profiles for 21 cell types, encompassing immune cells, endothelial cells, and fibroblasts. They then coupled this novel signature with the EPIC deconvolution framework based on constrained least-square regression to derive a dedicated tool called EPIC-ATAC. The method was then assessed using real and pseudo-bulk RNA-seq data from human peripheral blood mononuclear cells (PBMC) and, finally, applied to ATAC-seq data from breast cancer tumors to show it accurately quantifies their immune contexture.

Strengths:

Overall, the work is of very high quality. The proposed tool is timely; its implementation, characterization, and validation are based on rigorous methodologies and resulted in robust results. The newly-generated, validation data and the code are publicly available and well-documented. Therefore, I believe this work and the associated resources will greatly benefit the scientific community.

Weaknesses:

A few aspects can be improved to clarify the value and applicability of the EPIC-ATAC and the transparency of the benchmarking analysis.

Most of the validation results in the main text assess the methods on all cell types together, by showing the correlation, RMSE, and scatterplots of the estimated vs. true cell fractions. This approach is valuable for showing the overall method performance and for detecting systematic biases and noisy estimates. However, it provides very limited insights regarding the capability of the methods to estimate the individual cell types, which is the ultimate aim of deconvolution analysis. This limitation is exacerbated for rare cell types, which could even have a negative correlation with the ground truth fractions, but not weigh much on the overall RMSE and correlation. I would suggest integrating into the main text and figures an in-depth assessment of the individual cell types. In particular, it should be shown and discussed which cell types can be accurately quantified and which ones are less reliable.

In the benchmarking analysis, EPIC-ATAC is compared to several deconvolution methods, most of which were originally developed for transcriptomics data. This comparison is not completely fair unless their peculiarities and the limitations of tweaking them to work with ATAC-seq data are discussed. For instance, some methods (including the original EPIC) correct for cell-type-specific mRNA bias, which is not present in ATAC-seq data and might, thus, result in systematic errors.

On a similar note, it could be made more explicit which adaptations were introduced in EPIC, besides the ad-hoc ATAC-seq signature, to make it applicable to this type of data.

Given that the final applicability of EPIC-ATAC is on real bulk RNA-seq data, whose characteristics might not be completely recapitulated by pseudo-bulk samples, it would be interesting to see EPIC and EPIC-ATAC compared on a dataset with matched, real bulk RNA-seq and ATAC-seq, respectively. It would nicely complement the analysis of Figure 7 and could be used to dissect the commonalities and peculiarities of these two approaches.

Reviewer #1 (Public Review):

Cellulose is the major component of the plant cell wall and as such is a major component of all plant biomass on the planet. It is made at the cell surface by a large membrane-bound complex known as the cellular synthase complex. It is the structure of the cellulose synthase complex that determines the structure of the cellulose microfibril, the unit of cellulose found in nature. Consequently, while understanding the molecular structure of individual catalytic subunits that synthesise individual beta 1-4 glucose chains is important, to really understand cellulose synthesis it is necessary to understand the structure of the entire complex.

In higher plants, cellulose is synthesised by a large membrane-bound complex composed of three different CESA proteins. During cellulose synthesis in the primary cell wall, this is composed of members of groups CESA1, CESA3, and CESA6. While the authors have previously presented structural data on CESA8, required for cellulose synthesis in the secondary cell wall, here they provide structural and enzymatic analysis of CESA1, CESA3, and CESA6 from soya beans.

The authors have utilised their established protocol to purify trimers for all three classes of CESA proteins and obtain structural information using electron microscopy. The structures reveal some subtle, but interesting differences between the structures obtained in this study and that previously obtained for CESA8. In particular, they identify a change in the position of transmembrane helices 7 that in previous structures formed part of the transmembrane channel. In the structure of CESA1 TM7 is shifted laterally to a position more towards the periphery of the protomer, where it is stabilised by inter-protomer interactions. This creates a large lipid-exposed channel opening that is likely encountered by the growing cellulose chain. In the discussion, the authors speculate this channel might facilitate lateral movement of cellulose chains in the membrane which would allow them to associate to form the microfibril. There is, however, no explanation for why this might be different for CESA proteins involved in primary and secondary cell wall CESA proteins.

Interactions within the trimer as stabilised by the plant conserved regions (PCR), while in common with previous studies that class-specific regions (CSR) are not resolved, are likely highly disordered as has been suggested in previous studies. As the name suggests these regions are likely to be important for determining how different CESA proteins interact, but it remains to be seen how they achieve this. Similarly, the N-terminal domain (NTD) remains rather intriguing. In the CESA3 structure, the NTD forms a stalk that protrudes into the cytoplasm that was previously observed for CESA8, while it remains unresolved in CESA1 and CESA6. The authors suggest the inability to resolve this region is likely the result of the NTD being able to form multiple conformations. Loss of the NTD does not prevent the formation of trimers and CESA1 and CESA3 are still able to interact. Previous bioinformatic studies suggest that the CSR part of the NTD is also highly class-specific (Carrol et al. 2011 Frontiers in Plant Science 2, 5-5) suggesting it is also likely to participate in interactions between different CESA proteins. This analysis provides little new information on the structure of the NTD or how it functions as part of the cellulose synthase complex.

The other important point regarding cellulose synthesis is how the different CESA trimers function during cellulose synthesis and complex assembly. The authors provide biochemical evidence that mixed complexes of two different CESA proteins are able to synergistically increase the rate of cellulose synthesis. This increase is not dramatic, around 2-fold as it is unclear what brings about this increase and whether it results from the ability to form larger complexes favouring greater rates of cellulose synthesis.

It is clear however from electron microscopy that mixing of CESA proteins can lead to the formation of large aggregates not seen with single CESA proteins. The aggregates observed do not form rosette-type shapes but appear to be much more random aggregates of different CESA trimers. The authors suggest that this is likely a result of the fact that the complexes are not constrained in two dimensions by the membrane, however, if these are biologically relevant interactions that form aggregates it is somewhat surprising that they do not form hexameric structures, particularly since they are essentially forming as a single layer.

Overall the study provides some important data and raises a number of important questions.

Reviewer #1 (Public Review):

Summary:

This manuscript successfully established a new short-term model of diabetic retinopathy by treating zebrafish embryos with high concentrations of monosaccharides, resembling the hyperangiogenic characteristics observed in proliferative diabetic retinopathy in patients. They found that excessive angiogenesis induced by glucose and noncaloric monosaccharides was achieved by activating the quiescent endothelial cells into proliferating tip cells. Importantly, the authors further confirmed the effects of monosaccharides on inducing excessive angiogenesis were mediated by the foxo1a-marcksl1a pathway.

Strengths:

These results showed the potentially detrimental effects of the noncaloric monosaccharides on blood vessel function and provided novel insights into the underlying mechanisms.

Weaknesses:

The mechanism of noncaloric monosaccharides inducing excessive sprouting angiogenesis is not solid enough.

Reviewer #1 (Public Review):

Summary:

This impressive study by Bandet and Winship uses 2-photon imaging in awake behaving mice to examine long-term changes in neural activity and functional connectivity after focal ischemic stroke. The authors discover that there are long-lasting perturbations in neural activity and functional connectivity, specifically within peri-infarct cortex but not more distant cortical regions. Overall I thought the study provided important new findings that were supported by compelling data.

Strengths:

This is a technically challenging study and the experiments appear to be well done. The manuscript was written in a concise manner, and the figures were clearly presented. The analytic tools were rigorous and appropriate, leading to novel insights regarding neural activity patterns during movement or rest. The discovery of long-lasting impairments in neural activity/functional connectivity is important (and often overlooked) given that future stroke studies need to recognize what problems exist in order to properly rectify them. The authors also question the spatial extent to which functional changes occur after stroke, at least at the single cell level. Overall, I think this was a well-executed study whose primary conclusions were justified by the data presented.

Weaknesses:

I found very little in the way of weaknesses. The authors addressed my comments about the methodology, statistical analysis, normalization of data and discussion points about cortical plasticity during stroke recovery.

Reviewer #1 (Public Review):

Summary:

This is an interesting study by Xu et al showing the effects of infection with the Treponema pallidum virus (which causes syphilis disease) on neuronal development using iPSC-derived human brain organoids as a model and single-cell RNA sequencing. This work provides an important insight into the impact of the virus on human development, bridging the gap between the phenomena observed in studies using animal models as well as non-invasive human studies showing developmental abnormalities in fetuses infected with the virus in utero through maternal vertical transmission.

Using single-cell RNAseq in combination with qPCR and immunofluorescence techniques, the authors show that T. pallidum infected organoids are smaller in size, in particular during later growth stages, contain a larger number of undifferentiated neuronal lineage cells, and exhibit decreased numbers of specific neuronal subcluster, which the authors have identified as undifferentiated hindbrain neurons.

The study is an important first step in understanding how T. pallidum affects human neuronal development and provides important insight into the potential mechanisms that underlie the neurodevelopmental abnormalities observed in infected human fetuses. Several important weaknesses have also been noted, which need to be addressed to strengthen the study's conclusions.

Strengths:

(1) The study is well written, and the data quality is good for the most part.

(2) The study provides an important first step in utilizing human brain organoids to study the impact of T. pallidum infection on neuronal development.

(3) The study's conclusions may provide important insight to other researchers focused on studying how viral infections impact neuronal development.

Weaknesses:

(1) It is unclear how T. pallidum infection was validated in the organoids. If not all cells are infected, this could have important implications for the study's conclusions, in particular the single-cell RNAseq experiments. Were only cells showing the presence of the virus selected for sequencing? A detailed description of how infection was validated and the process of selection of cells for RNAseq would strongly support the study's conclusions.

(2) The authors show that T. pallidum infection results in impaired development of hindbrain neurons. How does this finding compare to what has already been shown in animal studies? Is a similar deficit in this brain region observed with this specific virus? It would be useful to strengthen the study's conclusions if the authors added a discussion about the observed deficits in hindbrain neuronal development, and prior literature on similar studies conducted in animal models or human patients. Does T. pallidum preferentially target these neurons, or is this a limitation of the current organoid model system?

(3) The authors show that T. pallidum-infected organoids are smaller in size by measuring organoid diameter during later stages of organoid growth, with no change during early stages. Does that represent insufficient infection at the early stages? Is this due to increased cell death or lack of cell division in the infected organoids? Experiments using IHC to quantify levels of cleaved caspase and/or protein markers for cell proliferation would be able to address these questions.

4) In Figure 1D authors show differences in rosette-like structure in the infected organoids. The representative images do not appear to be different in any of the discussed components (e.g., the sox2 signal looks fairly similar between the two conditions). No quantification of these structures was presented. Authors should provide quantification or a more representative image to support their statement.

5) The IHC images shown in Figures 3E, G, and Figure 4E look very similar between the two conditions despite the discussed decrease in the text. A more suitable representative image should be presented, or the analysis should be amended to reflect the observed results.

Reviewer #1 (Public Review):

Summary:

This manuscript aimed to investigate the emergence of emotional sensitivity and its relationship with gestational age. Using an oddball paradigm and event-related potentials, the authors conducted an experiment in 120 healthy neonates with a gestational age range of 35 to 40 weeks. A significant developmental milestone was identified at 37 weeks gestational age, marking a crucial juncture in neonatal emotional responsiveness.

Strengths:

This study has several strengths, by providing profound insights into the early development of social-emotional functioning and unveiling the role of gestational age in shaping neonatal perceptual abilities. The methodology of this study demonstrates rigor and well-controlled experimental design, particularly involving matched control sounds, which enhances the reliability of the research. Their findings not only contribute to the field of neurodevelopment, but also showcase potential clinical applications, especially in the context of autism screening and early intervention for neurodevelopmental disorders.

Weaknesses:

More details should be provided in terms of inclusion and exclusion criteria for the participants, as well as missing data due to the non-cooperation of newborns during the experimental process. Potential differences between preterm and full-term infants are worth exploring. Several aspects of EEG data analyses and data interpretation should be better clarified.

Reviewer #1 (Public Review):

Summary:

In this manuscript, the authors describe a new pipeline to measure changes in vasculature diameter upon opt-genetic stimulation of neurons.

The work is interesting and the topic is quite relevant to better understand the hemodynamic response on the graph/network level.

Strengths:

The manuscript provides a pipeline that allows for the detection of changes in the vessel diameter as well as simultaneously allowing for the location of the neurons driven by stimulation.

The resulting data could provide interesting insights into the graph-level mechanisms of regulating activity-dependent blood flow.

The interesting findings include that vessel radius changes depend on depth from the cortical surface and that dilations on average happen closer to the activated neurons.

Weaknesses:

The utility of a pipeline depends on the generalization properties.

While the proposed pipeline seems to work for the data the authors acquired, it is unclear if this pipeline will actually generalize to novel data sets possibly recorded by a different microscope (e.g. different brand), or different imagining conditions (e.g. illumination or different imagining artifacts) or even to different brain regions or animal species, etc.

The authors provide a 'black-box' approach that might work well for their particular data sets and image acquisition settings but it is left unclear how this pipeline is actually widely applicable to other conditions as such data is not provided.

In my experience, without well-defined image pre-processing steps and without training on a wide range of image conditions pipelines typically require significant retraining, which in turn requires generating sufficient amounts of training data, partly defying the purpose of the pipeline.

It is unclear from the manuscript, how well this pipeline will perform on novel data possibly recorded by a different lab or with a different microscope.

Analysis

Some of the chosen analysis results seem to not fully match the shown data, or the visualization of the data is hard to interpret in the current form. Additionally, some measures seem not fully adapted to the current situation (e.g. the efficiency measure does not consider possible sources or sinks). Thus, some additional analysis work might be required to account for this.

Reviewer #2 (Public Review):

The authors set out to identify CAPs (Candidate Adaptive Polymorphyisms), i.e., simply put mutations that carry a potential functional advantage, and utilize computational methods based on the perturbation of C-alpha positions with an Elastic Network Model to determine if dynamics of CAP residues are different in any way.

The authors have addressed the main methodological concerns.

However, one point remains. The specific comparison of which CAPs have been previously identified by other means, particularly with other computational methods that look into dynamics is still lacking. It is unfortunate that the authors do not present such analysis, particularly with respect to single point mutational analysis of Teruel et al. in Plos Comp. Bio. If CAP positions were previously identified by other means it adds strength to the methodology used by the authors. The authors also do not discuss their results in light of the work of Lam et al. (Sci. Comm, 2020) where an evolutionary analysis of Spike/ACE2 binding across homologues is performed. I believe that such deeper discussion of the current results in light of previous work, adds strength to the analysis presented in this manuscript as the methodology is different. Even if all results were not new, with the method being different from the other means by which such results were obtained, it would be still a worthy contribution to the field. Furthermore, for the community at large trying to understand the importance of particular positions in Spike, knowing that a particular position identified here was also identified by works X, Y, Z adds a lot of to the field. I can only think that the authors may imagine that if one of their CAPs was identified by other means previously, it takes away from the merit of their work, but it is actually the opposite. I urge the authors to not brush away this. In fact, more important than any methodological aspect of the present work, this strengthening of evidence for particular positions by several independent methods is the most important evidence that the authors can contribute to the field.

Reviewer #1 (Public Review):

The work by Ginatt et al. uses genome-scale metabolic modeling to identify and characterize trophic interactions between rhizosphere-associated bacteria. Beyond identifying microbial species associated with specific host and soil traits (e.g., disease tolerance), a detailed understanding of the interactions underlying these associations is necessary for developing targeted microbiome-centered interventions for plant health. It has nonetheless remained challenging to define the roles of specific organisms and metabolic species in natural rhizobiomes. Here, the authors combine microbial compositional data obtained through metagenomic sequencing with a new collection of genome-scale models to predict interactions in the native rhizosphere communities of apple rootstocks. To do this, they have established processes to integrate these sources of data and model specific trophic exchanges, which they use to obtain testable hypotheses for targeted modulation of microbiota members in situ.

The authors carry out a careful model curation process based on metagenomic sequencing data and existing model generation tools, which, together with basing the in silico medium composition on known root exudates, strengthens their predictions of interaction network features. Moreover, its reliance on genome-scale models provides a broader basis for linking sequence-based information to predictions of function on a multispecies level beyond rhizosphere microbiomes.

Having generated a set of predicted trophic interactions, the authors carried out a detailed analysis linking features of these interactions to organism taxonomy and broader ecosystem properties. Intriguingly, the organisms predicted to grow in the first iteration of their framework (i.e., on only root exudates) broadly correspond to taxonomic groups experimentally shown to benefit from these compounds. Additionally, the simulations predicted some patterns of vitamin and amino acid secretion that are known to form the basis for interactions in the rhizosphere. Together, these outcomes underscore the applicability of this method to help disentangle trophic interaction networks in complex microbiomes.

The methodology described in this paper represents a useful and promising framework to better understand the complexity of microbial interaction networks in situ. However, the degree to which the predictions can vary according to environmental composition remains difficult to quantify, and the work does not address the sensitivity of the modeling predictions beyond a simulated medium containing 33 root exudates. I find this especially important given that relatively few (84 of 243) species were predicted to grow even after cross-feeding, suggesting that a richer medium could lead to different interaction network structures. While the authors do state the importance of environmental composition and have carefully designed an in silico medium, I believe that simulating a broader set of resource pools would add necessary insight into both the predictive power of the models themselves and trophic interactions in the rhizosphere more generally.

Reviewer #1 (Public Review):

Summary:

Gain-of-function mutations and amplifications of PPM1D are fond across several human cancers and are associated with advanced tumor stage, worse prognosis, and increased lymph node metastasis. This manuscript presents important findings that SOD1 inhibition is a potential strategy to achieve therapeutic synergism for PPM1D-mutant leukemia; and demonstrates the redox landscape of PPM1D-mutant cells.

Strengths:

In this manuscript, Zhang and colleagues investigate the synthetic-lethal dependencies of PPM1D (protein phosphatase, Mg2+/Mn2+ dependent 1D) in leukemia cells using CRISPR/Cas9 screening. They identified that SOD1 (superoxide dismutase-1) as the top hit, whose loss reduces cellular growth in PPM1D-mutant cells, but not wildtype (WT) cells. Consistently, the authors demonstrate that PPM1D-mutant cells are more sensitive to SOD1 inhibitor treatment. By performing different in vitro studies, they show that PPM1D-mutant leukemia cells have elevated level of reactive oxygen species (ROS), decreased basal respiration, increased genomic instability, and impaired non-homologous end-joining repair. These data highlight the potential of SOD1 inhibition as a strategy to achieve therapeutic synergism for PPM1D-mutant leukemia; and demonstrates the redox landscape of PPM1D-mutant cells.

Weaknesses:

While the current study has identified synthetic lethality of PPM1D-mutant leukemia cells upon SOD1 inhibition, the underlying mechanism remains elusive. Although ROS levels have been assessed between wild-type (WT) and PPM1D-mutant leukemia cells, the specific redox alterations induced by SOD1 inhibition in PPM1D mutant versus WT cells have not been elucidated. To address this gap, direct comparisons of ROS levels using various probes should be conducted between PPM1D mutant and WT cells under conditions of SOD1 inhibition.

Reviewer #1 (Public Review):

Warming and precipitation regime change significantly influences both above-ground and below-ground processes across Earth's ecosystems. Soil microbial communities, which underpin the biogeochemical processes that often shape ecosystem function, are no exception to this, and although research shows they can adapt to this warming, population dynamics and ecophysiological responses to these disturbances are not currently known. The Qinghai-Tibet Plateau, the Third Pole of the Earth, is considered among the most sensitive ecosystems to climate change. The manuscript described an integrated, trait-based understanding of these dynamics with the qSIP data. The experimental design and methods appear to be of sufficient quality. The data and analyses are of great value to the larger microbial ecological community and may help advance our understanding of how microbial systems will respond to global change. There are very few studies in which the growth rates of bacterial populations from multifactorial manipulation experiments on the Qinghai-Tibet Plateau have been investigated via qSIP, and the large quantity of data that comprises the study described in this manuscript, will substantially advance our knowledge of bacterial responses to warming and precipitation manipulations.

Showing the parallel between him and other babies of his age, he describes that he did not shout "Miez, Miez, Miez, Miez", which translates to "no, no, no, no" in old germanic dialects like the other children. But instead, in a polar opposite manor, yelled "high, sturdy and big voice" shouted drink, drink, drink. Showing from an early age leadership, confidence, and "greatness".

Reviewer #1 (Public Review):

Osteoarthritis (OA) is associated with painful, chronic inflammation that often leads to severe joint pain and joint stiffness for people over the age of 55. There is no effective therapeutic drug in the treatment of osteoarthritis. The authors found that mice without Cbfβ in their chondrocytes develop spontaneous OA. Authors uncovered that the deficiency of Cbfβ caused increased canonical Wnt signaling and inflammatory response, and decreased Hippo/YAP signaling and TGF-β signaling in articular cartilage. Authors showed that ACLT surgery-induced OA decreased Cbfβ and Yap expression and increased active β-catenin expression in articular cartilage, while local AAV-mediated Cbfβ overexpression promoted Yap expression, diminished active β-catenin expression in OA lesions. The authors demonstrated that AAV-mediated Cbfβ overexpression in knee joints of mice with OA showed the significant protective effect of Cbfβ on articular cartilage in the ACLT OA mouse model. The results from the study demonstrated Cbfβ maintains articular cartilage homeostasis through inhibiting Wnt/β-catenin signaling and increasing Hippo/Yap, and TGFβ signaling. Importantly, the authors proved that local Cbfβ overexpression could be an effective strategy for treatment of OA. The data shown in the study demonstrated that the findings are novel and very significant, and the authors' claims and conclusions are justified by their data. The paper is generally excellent with an interesting scientific premise and strong scientific rigor. The findings in this manuscript are novel, the manuscript is clearly written, and the findings will make a significant impact in the field.

Reviewer #1 (Public Review):

In this paper, Nikolaou et al. demonstrated that CYRI-B expression is upregulated in a mouse model of pancreatic ductal adenocarcinoma (PDAC). Interestingly, they found that, while CYRI-B KO promotes the early stages of tumour progression, it prevents the formation of metastasis at later stages. Focusing on the latter, the authors highlight a role for CYRI-B in controlling the membrane availability of the LPA receptor LPAR1, which is required to support PDAC cell chemotaxis towards serum or LPA.

Strengths: the in vivo and imaging data are very solid, and convincingly support the authors' conclusions. The KPC model is well-established in PDAC research and is a very powerful tool to investigate disease onset and progression. The imaging approaches used are of a very high standard. Good data presentation with the use of super-plots.

Weaknesses: the authors focused on chemotaxis, but did not present any evidence with regard to the role of CYRI-B in 3D cell invasion, which is a key process associated with cell invasion. The data presented clearly show a specific effect towards liver metastasis, while diaphragm and bowel metastasis were not affected by CYRI-B deletion. It would be beneficial to include a discussion about this, providing some potential explanation behind this observation.

This work is of interest to cell biologists not only working in pancreatic cancer but also more broadly to researchers interested in vesicular trafficking, plasma membrane receptor dynamics and cell migration.

Reviewer #1 (Public Review):

Hoving and colleagues investigated the mechanisms of contact inhibition of locomotion (CIL) in Schwann cells using cell migration assays, in combination with siRNA as well as an ex-vivo model for collective cell migration of the peripheral nervous system. They found that N-cadherin is needed for proper cell repulsion during CIL. Schwann Cells depleted of N-cadherin failed CIL when encountering other Schwann cells depleted of N-cadherin, however they maintained CIL when encountering Schwann cells expressing N-cadherin. Depletion of alpha-catenin and to some degree p120 did not have the same effect as N-cadherin depletion. Further, they determined that the extracellular domain is needed for CIL as well as an interaction with Glypican-4. Glypicans often act as co-receptors for other signaling molecules, and so the authors further narrowed CIL's dependence to Slit signaling. N-cadherin was needed for proper Slit surface expression, again, dependent on the extracellular domain, and depletion of both Slit2 and 3 lead to a cell clumping and rounding phenotype. Finally, using an ex-vivo model of Schawnn cell migration they showed that rSlit lead to a similar cell rounding and clumping phenotype, ultimately leading to an inhibition of cell migration.

Strengths

This was a very methodical examination of what is needed for CIL in cultured Schwann cells. The data presented largely supports the findings and the linking of N-cad to glypican-4 to Slit signaling further illuminates this process helping to define the molecular players. The mechanistic insight goes further in that they demonstrate the Slit does not get to the cell surface without the expression of the extracellular domain of N-cad.

Weaknesses

The conclusions that can be drawn from this study remain a little narrow since only Schwann cells were used. This is not so much a weakness in that authors were indeed investigating the periphery nervous system regeneration but it does limit their findings. The experiments carried out in the ex-vivo system only touch on one aspect of their cell culture work, the mechanism of Slit. No other aspects of their cell culture system was tested ex-vivo which

Reviewer #1(Public Review):

In this manuscript the authors report an experiment to assess how training on a perceptual task may not only increase performance on that task but impact on the appearance of the trained stimuli. They compare discrimination performance, coherence thresholds, and estimation biases for random dot motion direction relative to horizontal rightward in three groups of observers before and after 3 days in which they either trained on a discrimination task, an estimation task, or did not train. The authors report significant increases in discrimination performance post training compared to not training. They also report increases in estimation biases when assessed as the average estimate (over a bimodal distribution that crosses 0) but not when assessed as the mode of the bimodal distribution. They conclude that training resulted in "increases in already-large estimation biases away from horizontal".

The methods and results are strengthened by the combination of classical psychophysical techniques and sophisticated computational modelling. One weakness is the possibility is misleading summary statistics when dealing with bimodal distributions. Convincing evidence that observers perceived stimulus directions as further from horizontal (in the absolute sense) following training is not presented in the current manuscript. Irrespective, this work is likely to impact the field.

Reviewer #1 (Public Review):

Identifying individual BCR/Ab chain sequences that are members of the same clone is a long-standing problem in the analysis of BCR/Ab repertoire sequencing data. The authors propose a new method designed to be scalable for application to huge repertoire data sets without sacrificing accuracy. Their approach utilizes Hamming Distance between CDR3 sequences followed by clustering for a fast, high-precision approach to classifying pairs of sequences as related or not, and then refines the classification using mutation information from germline-encoded regions. They compare their method with other state-of-the-art methods using synthetic data.

The authors address an important problem in an interesting, innovative, and rigorous way, using probabilistic representations of CDR3 differences, frequencies of shared and not-shared mutations, and the relationships between the two under hypotheses of related pairs and unrelated pairs, and from these develop an approach for determining thresholds for classification and lineage assignment. Benchmarking shows that the proposed method, the complete method including both steps, outperforms other methods.

Strengths of the method include its theoretical underpinnings which are consistent with an immunologist's intuition about how related and unrelated sequences would compare with each other in terms of the metrics to use and how those metrics are related to each other.

I have two high-level concerns:<br /> (1) It isn't clear how the real and synthetic data are being used to estimate parameters for the classifier and evaluate the classifier to avoid circularity. It seems like the approach is used to assign lineages in the data from [1], and then properties of this set of lineages are used to estimate parameters that are then used to refine the approach and generate synthetic data that is used to evaluate the approach. This may not be a problem with the approach but rather with its presentation, but it isn't entirely clear what data is being used and where for what purpose. An understanding of this is necessary in order to truly evaluate the method and results.<br /> (2) Regarding the data used for benchmarking - given the intertwined fashion by which the classification approach and synthetic data generation approach appear to have been developed, it is not surprising that the proposed approach outperforms the other methods when evaluated on the synthetic data presented here. It would be better to include in the benchmark the data used by the other methods to benchmark themselves or also generate synthetic data using their data generation procedures.

An improved method for BCR/Ab sequence lineage assignment would be a methodologic advancement that would enable more rigorous analyses of BCR/Ab repertoires across many fields, including infectious disease, cancer, autoimmune disease, etc., and in turn, enable advancement in our understanding of humoral immune responses. The methods would have utility to a broad community of researchers.

Reviewer #1 (Public Review):

Summary:

Doxorubincin has long been known to cause bone loss by increasing osteoclast and suppressing osteoblast activities. The study by Wang et al. reports a comprehensive investigation into the off-target effects of doxorubicin on bone tissues and potential mechanisms.. They used a tumor-free model with wild type mice and found that even a single dose of doxorubicin has a major influence by increasing leukopenia and DAMPs and inflammasomes in macrophages and neutrophils, and inflammation-related cell death (pyroptosis and NETosis). The gene knockout study shows that AIM2 and NLRP3 are the major contributors to bone loss. Overall, the study confirmed previous findings regarding the impact of doxorubicin on tissue inflammation and expands the research further into bone tissue. The presented data presented are consistent; however, a major question remains regarding whether doxorubicin drives inflammation and its related events. Most in vitro study showed that the effect of doxorubincin cannot be demonstrated without LPS priming. This observation raises the question of whether doxorubincin itself could activate the inflammasome and the related events. In vivo study, on the other hand, suggested that it doesn't require LPS. The inconsistency here was not explained further. Moreover, a tumor-free mouse model was used for the study; however, immune responses in tumor bearing models would likely be distinct from tumor-free ones. The justification for using tumor-free models is not well-established.

Strengths:

The paper includes a comprehensive study that shows the effects of doxorubincin on cytokine levels in serum, release of DAMPs and NETosis, and leukopenia using both in vivo and in vitro models. Bone marrow cells, macrophages and neutrophils were isolated from the bone marrow, and the levels of cytokines in serum were also determined.

They employed multiple knockout models with deficiency in Aim 2, Nlirp3, and double deficiencies to dissect the functional involvement of these two inflammasomes.

The experiments in general are well designed. The paper is also logically written, and figures were clearly labeled.

Weaknesses:

Most of the data presented are correlative, and there is not much effort to dissect the underlying molecular mechanism.

It is not entirely clear why a tumor free model is chosen to study immune responses, as immune responses can differ significantly with or without tumor-bearing.

Immune responses in isolated macrophages, neutrophils and bone marrow cells require priming with LPS, while such responses are not observed in vivo. There is no explanation for these differences.

The band intensities on Western blots in Fig. 4 and Fig. 5 are not quantified, and the numbers of repeats are also not provided.

Many abbreviations are used throughout the text, and some of the full names are not provided.

Fig. 5B needs a label on X axis.

Reviewer #1 (Public Review):

In this manuscript by Buchanan and colleagues, the authors set out to determine if mutations associated with resistance to the Plasmodium apicoplast inhibitor azithromycin (AZ) had a measurable impact on the fitness of Plasmodium berghei and P. falciparum parasites as they traverse both the mosquito host and vertebrate liver.

The Plasmodium endosymbiotic organelles - the mitochondrion and apicoplast - are attractive drug targets as they (1) possess essential functions across the multi-host multi-compartment life cycle of these parasites, and (2) are of bacterial origin and thus are vulnerable to inhibition both to extant antibiotics, and novel drugs with high parasite specificity.

Historically however the high resistance propensity of drug targets encoded in the organellar genomes (most notably atovaquone and doxycycline) has precluded the use of these drugs in an endemic setting, limiting these potent compounds to use in prophylaxis for travelers from non-endemic countries. Several studies in the last decade now fairly definitively show that mutations conferring resistance to atovaquone in the mitochondrial gene cytochrome b are, in a mutation-dependent manner, totally or near-totally compromised in their ability to infect, grow, and escape the mosquito host, leading to a reexamination of the potential utility of this extraordinarily potent drug in endemic settings. Symmetries exist between the Plasmodium mitochondrion and apicoplast, which both appear to have highly fexpanded roles in the mosquito and liver relative to the blood stages. Thus, the authors set out to explore whether mutations in essential apicoplast genes were, in a similar manner to mutations in cytochrome B, associated with fitness effects in the mosquito and/or liver.

Towards this, the authors selected for several AZ-resistant parasite populations, all of which acquired mutations in the apicoplast genome-encoded ribosomal protein Rpl4. Interestingly, the authors observed contrasting fitness effects caused by these mutations, both between mutants within Plasmodium species, and between species. In P. berghei, AZ mutants were compromised in their ability to form oocysts and sporozoites, and a large proportion of sporozoites lacked an intact apicoplast and displayed aberrant gliding behaviour. Similarly, in the liver, Rpl4 mutant P. berghei liver schizonts were smaller, had fewer nuclei, and appeared extremely limited in their ability to cause a patent infection - crucially in particular via mosquito bites. Surprisingly, a P. falciparum Rpl4 mutant (notably in a different position of the protein) had no impact on sporogony but appeared to have a strong impact on liver schizont development in a liver-humanized mouse model, suggesting that establishment of blood stage infection in a subsequent human host would be less likely for mutant parasites.

This is a well-executed study, that presents novel and noteworthy findings. The impact of drug-resistance-conferring mutations in Plasmodium outside of the blood stage is woefully understudied, primarily due to significant challenges associated with studying Plasmodium, especially P. falciparum, in both the mosquito and liver which the authors navigate commendably. The results presented in this manuscript leverage state-of-the-art techniques and clearly support the authors' conclusion that AZ-conferring resistance mutations have a strong negative effect on the ability of Plasmodium parasites to both reinfect and cause symptomatic infection in a subsequent vertebrate host. This could indicate that apicoplast-targeted inhibitors are more attractive as co-drugs for malaria treatment than previously thought, due to the reduced probability of the spread of resistance, which has been a perennial issue in malaria therapeutic care.

Reviewer #1 (Public Review):

Mohseni and Elhaik's article offers a critical evaluation of Geometric Morphometrics (GM), a common tool in physical anthropology for studying morphological differences and making phylogenetic inferences. I read their article with great interest, although I am not a geneticist or an expert on PCA theory since the problem of morphology-based classification is at the core of paleoanthropology.

The authors developed a Python package for processing superimposed landmark data with classifier and outlier detection methods, to evaluate the adequacy of the standard approach to shape analysis via modern GM. They call into question the accuracy, robustness, and reproducibility of GM, and demonstrate how PCA introduces statistical artefacts specific to the data, thus challenging its scientific rigor. The authors demonstrate the superiority of machine learning methods in classification and outlier detection tasks. The paper is well-written and provides strong evidence in support of the authors' argument. Thus, in my opinion, it constitutes a major contribution to the field of physical anthropology, as it provides a critical and necessary evaluation of what has become a basic tool for studying morphology, and of the assumptions allowing its application for phylogenetic inferences. Again, I am not an expert in these statistical methods, nor a geneticist, but the authors' contribution is of substantial relevance to our field (physical anthropology). The examples of NR fossils and HLD 6 are cases in point, in line with other notable examples of critical assessment of phylogenetic inferences made on the basis of PCA results of GM analysis. For example, see Lordkipanidze et al.'s (2014) GM analyses of the Dmanisi fossils, suggesting that the five crania represent a single regional variant of Homo erectus; and see Schwartz et al.'s (2014) comment on their findings, claiming that the dental, mandibular, and cranial morphology of these fossils suggest taxic diversity. Schwartz et al. (2014) ask, "Why did the GMA of 78 landmarks not capture the visually obvious differences between the Dmanisi crania and specimens commonly subsumed H. erectus? ... one wonders how phylogenetically reliable a method can be that does not reflect even easily visible gross morphological differences" (p. 360).

As an alternative to the PCA step in GM, the authors tested eight leading supervised learning classifiers and outlier detection methods on three-dimensional datasets. The authors demonstrated inconsistency of PCA clustering with the taxonomy of the species investigated for the reconstruction of their phylogeny, by analyzing a database comprising landmarks of 6 known species that belong to the Old World monkeys tribe Papionini, using PCA for classification. The authors also demonstrated that high explained variance should not be used as an estimate of high accuracy (reliability). Then, the authors altered the dataset in several ways to simulate the characteristic nature of paleontological data.

The authors excluded taxa from the database to study how PCA and alternative classifiers are affected by partial sampling, and the results presented in Figures 4 and 5, among others, are quite remarkable in showing the deviations from the benchmark data. These results expose the perils of applying PCA and GM for interpreting morphological data. Furthermore, they provide evidence showing that the alternative classifiers are superior to PCA, and that they are less susceptible to experimenter intervention. Similar results, i.e., inconsistencies in the PC plots, were obtained in examinations of the effect of removing specimens from the dataset and in the interesting test of removing landmarks to simulate partial morphological data, as is often the case with fossils. To test the combined effect of these data alterations, the authors combined removal of taxa, specific samples, and landmarks from the dataset. In this case, as well, the PCA results indicate deviation from the benchmark data. However, the ML classifiers could not remedy the situation. The authors discuss how these inconsistencies may lead to different interpretations of the data, and in turn, different phylogenetic conclusions. Lastly, the authors simulated the situation of a specimen of unknown taxonomy using outlier detection methods, demonstrating LOF's ability to identify a novelty in the morphospace.

References<br /> Bookstein FL. 1991. Morphometric tools for landmark data: geometry and biology [Orange book]. Cambridge New York: Cambridge University Press.<br /> Cooke SB, and Terhune CE. 2015. Form, function, and geometric morphometrics. The Anatomical Records 298:5-28.<br /> Lordkipanidze D, et al. 2013. A complete skull from Dmanisi, Georgia, and the evolutionary biology of early Homo. Science 342: 326-331.<br /> Schwartz JH, Tattersall I, and Chi Z. 2014. Comment on "A complete skull from Dmanisi, Georgia, and the evolutionary biology of Early Homo". Science 344(6182): 360-a.

Joint Public Review:

Roget et al. build on their previous work developing a simple theoretical model to examine whether ageing can be under natural selection, challenging the mainstream view that ageing is merely a byproduct of other biological and evolutionary processes. The authors propose an agent-based model to evaluate the adaptive dynamics of a haploid asexual population with two independent traits: fertility timespan and mortality onset. Through computational simulations, their model demonstrates that ageing can give populations an evolutionary advantage. Notably, this observation arises from the model without invoking any explicit energy tradeoffs, commonly used to explain this relationship.

Additionally, the theoretical model developed here indicates that mortality onset is generally selected to start before the loss of fertility, irrespective of the initial values in the population. The selected relationship between the fertility timespan and mortality onset depends on the strength of fertility and mortality effects, with larger effects resulting in the loss of fertility and mortality onset being closer together. By allowing for a trans-generational effect on ageing in the model, the authors show that this can be advantageous as well, lowering the risk of collapse in the population despite an apparent fitness disadvantage in individuals. Upon closer examination, the authors reveal that this unexpected outcome is a consequence of the trans-generational effect on ageing increasing the evolvability of the population (i.e., allowing a more effective exploration of the parameter landscape), reaching the optimum state faster.

The simplicity of the proposed theoretical model represents both the major strength and weakness of this work. On one hand, with an original and rigorous methodology, the logic of their conclusions can be easily grasped and generalised, yielding surprising results. Using just a handful of parameters and relying on direct competition simulations, the model qualitatively recapitulates the negative correlation between lifespan and fertility without requiring energy tradeoffs. This alone makes this work an important milestone for the rapidly growing field of adaptive dynamics, opening many new avenues of research, both theoretically and empirically.

On the other hand, the simplicity of the model also makes its relationship with living organisms difficult to gauge, leaving open questions about how much the model represents the reality of actual evolution in a natural context. In particular, a more explicit discussion on how the specifics of the model can impact the results and their interpretation is needed. For example, the lack of mechanistic details on the trans-generational effect on ageing makes the results difficult to interpret. Even if analytical results are obtained, most of the observations appear derived from simulations as they are currently presented. Also, the choice of parameters for the simulations shown in the paper and how it relates to our biological knowledge is not fully addressed by the authors. Finally, the conclusions of evolvability are insufficiently supported, as the authors do not show if the wider genotypic variability in populations with the ageing trans-generational effect is, in fact, selected.

Reviewer #1 (Public Review):

Summary:

In this study, Faniyan and colleagues build on their recent finding that renal Glut2 knockout mice display normal fasting blood glucose levels despite massive glucosuria. Renal Glut2 knockout mice were found to exhibit increased endogenous glucose production along with decreased hepatic metabolites associated with glucose metabolism. Crh mRNA levels were higher in the hypothalamus while circulating ACTH and corticosterone was elevated in this model. While these mice were able to maintain normal fasting glucose levels, ablating afferent renal signals to the brain resulted in substantially lower blood glucose levels compared to wildtype mice. In addition, the higher CRH and higher corticosterone levels of the knockout mice were lost following this denervation. Finally, acute phase proteins were altered, plasma Gpx3 was lower, and major urinary protein MUP18 and its gene expression were higher in renal Glut2 knockout mice. Overall, the main conclusion that afferent signaling from the kidney is required for renal glut2 dependent increases in endogenous glucose production is well supported by these findings.

Strengths:

An important strength of the paper is the novelty of the identification of kidney to brain communication as being important for glucose homeostasis. Previous studies had focused on other functions of the kidney modulated by or modulating brain function. This work is likely to promote interest in CNS pathways that respond to afferent renal signals and the response of the HPA axis to glucosuria. Additional strengths of this paper stem from the use of incisive techniques. Specifically, the authors use isotope enabled measurement of endogenous glucose production by GC-MS/MS, capsaicin ablation of afferent renal nerves, and multifiber recording from the renal nerve. The authors also paid excellent attention to rigor in the design and performance of these studies. For example, they used appropriate surgical controls, confirmed denervation through renal pelvic CGRP measurement, and avoided the confounding effects of nerve regrowth over time. These factors strengthen confidence in their results. Finally, humans with glucose transporter mutations and those being treated with SGLT2 inhibitors show a compensatory increase in endogenous glucose production. Therefore, this study strengthens the case for using renal Glut2 knockout mice as a model for understanding the physiology of these patients.

Weaknesses:

A few weaknesses exist. Most concerns relate to the interpretation of this study's findings. The authors state that loss of glucose in urine is sensed as a biological threat based on the HPA axis activation seen in this mouse model. This interpretation is understandable but speculative. Importantly, whether stress hormones mediate the increase in endogenous glucose production in this model and in humans with altered glucose transporter function remains to be demonstrated conclusively. For example, the paper found several other circulating and local factors that could be causal. This model is also unable to shed light on how elevated stress hormones might interact with insulin resistance, which is known to increase endogenous glucose production. That issue is of substantial clinical relevance for patients with T2D and metabolic disease. Finally, while findings from the Glut2 knockout mice are of scientific interest, it should be noted that the Glut2 receptor is critical to the function of pancreatic islets and as such is not a good candidate for pharmacological targeting

Reviewer #1 (Public Review):

Summary:

The study investigated how root cap cell corpse removal affects the ability of microbes to colonize Arabidopsis thaliana plants. The findings demonstrate how programmed cell death and its control in root cap cells affect the establishment of symbiotic relationships between plants and fungi. Key details on molecular mechanisms and transcription factors involved are also given. The study suggests reevaluating microbiome assembly from the root tip, thus challenging traditional ideas about this process. While the work presents a key foundation, more research along the root axis is recommended to gain a better understanding of the spatial and temporal aspects of microbiome recruitment.