Reviewer #1 (Public Review):

This interesting manuscript from the Perozo and Faraldo-Gomez labs investigates the molecular mechanisms underlying the activation of the mechanosensitive ion channel MscS. The authors use a clever combination of cryoEM, coarse-grained (CG) and all-atom (AA) molecular dynamics simulations to determine the first (putatively) open conformation of the WT MscS channel and to show that this channel induces profound deformations of the membrane in the closed but not in the open state. Strikingly, MD simulations reveal that, contrary to what was previously assumed, lipids occupying cavities near the closed pore (hook lipids) come from the outer rather than inner leaflets. On pore opening, the membrane adopts a more relaxed conformation where the lipids contacting the protein are in less strained and tilted conformations. The authors thus propose a mechanism for sensing tension where the equilibrium between the open and closed conformations of the channel is dictated by differences in the membrane morphology in the two states rather than by the association and dissociation of individual lipids with the protein.

Major<br /> The observations on the hook lipids are critical and should be documented better. Based on previous work, it had been proposed that the hook lipids are associated with the inner leaflet and that they leave upon (partial) channel opening. In contrast, the present MD simulations indicate these lipids are associated with the outer leaflet and that their association to the channel persists on opening. These critical observations need to be documented better.<br /> i. Do the authors observe hook lipids in the cryoEM structure of the open channel? If yes, data should be shown. If no, then the discrepancy between MD and EM should be explicitly addressed.<br /> ii. Please show the comparison of the position and coordination of the hook lipids in MD simulations and in the closed (and/or open) structures.<br /> iii. The authors acknowledge that the volume of the cavity where the hook lipids are located decreases on channel opening. How does this not affect the association of the hook lipids with the protein?<br /> iv. Past work revealed several lipids in MscS structures near these cavities besides the hook lipids, and their ordered dissociation from the channel was proposed to be important for gating. Do the simulations show lipids in these cavities?<br /> v. Does the occupancy of the hook lipids in MD simulations change between the open and closed conformations? This should be analyzed.<br /> vi. Is the occupancy of other lipids in the nearby cavity altered upon channel opening?<br /> vii. Is the exchange of lipids near Ile150 affected by the conformational change?

I am a bit confused by the claim that "The comparison clearly highlights the reduction in the width of the transmembrane span of the channel upon opening, and how this changed is well matched by the thickness of the corresponding lipid nanodiscs (approximately from 38 to 23 Å)."<br /> i. How was the nanodisc membrane thickness determined? This should be described.<br /> ii. I do not see a ~15A change in the vertical length of the channel protein or of the nanodisc. While the panels in Fig.2 clearly show a vertical compression of the membrane, it appears that the ~15 A claim might be overstated. Adding a panel with measurements would be helpful to quantify this claim. If this is difficult on the membrane, maybe measurements could be performed on the protein.<br /> iii. What happens to the N-terminal cap structure in the open state? What are the rearrangements that allow the extracellular ends of the TM1 to disassemble the cap.

The data shown in Fig. 6 is cryptic and should be explained better in the main text. As it stands there is a cursory mention in pg. 12 and not much else.<br /> i. It would be helpful if the authors showed the position of Ile150 in the structure.<br /> ii. Does the total number of lipids in proximity of Ile150 change over time? Or the fold change represents ~1:1 exchange of lipids in the pocket?<br /> iii. I am confused by the difference in the maximum possible fold-change in unique lipids, does this reflect the difference in total number of lipids in each leaflet in each system? If so, I am a bit confused as to why there is a ~30% difference in the AA simulations whereas the values are nearly identical for the CG one.<br /> iv. Is it possible to quantify the residence time of the lipids in the pocket of each subunit?

The authors state on Pg. 21 "Nevertheless, we question the prevailing view that density signals of this kind are evidence of regulatory lipid binding sites; that is, we do not concur with the assumption that lipids regulate the gating equilibrium of MscS just like an agonist or antagonist would for a ligand-gated receptor-channel." I am a bit confused by this statement. In principle, binding and unbinding of modulatory ligands can happen on relatively fast time scales, so the observation that in MD simulations lipids exchange on a faster time scale than that of channel gating is not sufficient to make this inference. Indeed, there is ample evidence from other channels (i.e. Trp channels, HCN channels etc) where visualization of similar signals led to the identification of modulatory lipid binding sites. Thus, while I do not necessarily disagree with the authors, I would encourage them to tone down the general portion of the statement.

Reviewer #1 (Public Review):

The authors focused on linking physiological data on theta phase precession and spike-timing-dependent plasticity to the more abstract successor representation used in reinforcement learning models of spatial behavior. The model is presented clearly and effectively shows biological mechanisms for learning the successor representation. Thus, it provides an important step toward developing mathematical models that can be used to understand the function of neural circuits for guiding spatial memory behavior.

However, as often happens in the Reinforcement Learning (RL) literature, there is a lack of attention to non-RL models, even though these might be more effective at modeling both hippocampal physiology and its role in behavior. There should be some discussion of the relationship to these other models, without assuming that the successor representation is the only way to model the role of the hippocampus in guiding spatial memory function.

1. Page 1- "coincides with the time window of STDP" - This model shows effectively how theta phase precession allows spikes to fall within the window of spike-timing-dependent synaptic plasticity to form successor representations. However, this combination of precession and STDP has been used in many previous models to allow the storage of sequences useful for guiding behavior (e.g. Jensen and Lisman, Learning and Memory, 1996; Koene, Gorchetchnikov, Cannon, Hasselmo, Neural Networks, 2003). These previous models should be cited here as earlier models using STDP and phase precession to store sequences. They should discuss in terms of what is the advantage of an RL successor representation versus the types of associative sequence coding in these previous models.

2. On this same point, in the introduction, the successor representation is presented as a model that forms representations of space independent of reward. However, this independence of spatial associations and reward has been a feature of most hippocampal models, that then guide behavior based on interactions between a reward representation and the spatial representation (e.g. Redish and Touretzky, Neural Comp. 1998; Burgess, Donnett, Jeffery, O'Keefe, Phil Trans, 1997; Koene et al. Neural Networks 2003; Hasselmo and Eichenbaum, Neural Networks 2005; Erdem and Hasselmo, Eur. J. Neurosci. 2012). The successor representation should not be presented as if it is the only model that ever separated spatial representations and reward. There should be some discussion of what (if any) advantages the successor representation has over these other modeling frameworks (other than connecting to a large body of RL researchers who never read about non-RL hippocampal models). To my knowledge, the successor representation has not been explicitly tested on all the behaviors addressed in these earlier models.

3. Related to this, successes of the successor representation are presented as showing the backward expansion of place cells. But this was modeled at the start by Mehta and colleagues using STDP-type mechanisms during sequence encoding, so why was the successor representation necessary for that? I don't want to turn this into a review paper comparing hippocampal models, but the body of previous models of the role of the hippocampus in behavior warrants at least a paragraph in each of the introduction and discussion sections. In particular, it should not be somehow assumed that the successor representation is the best model, but instead, there should be some comparison with other models and discussion about whether the successor representation resembles or differs from those earlier models.

4. The text seems to interchangeably use the term "successor representation" and "TD trained network" but I think it would be more accurate to contrast the new STDP trained network with a network trained by Temporal Difference learning because one could argue that both of them are creating a successor representation.

Reviewer #1 (Public Review):

In this manuscript, the authors have assembled a reference transcriptome of the whole head of Loligo vulgaris and used it to perform single cell transcriptomics. With about 20,000 cells, they identify 32 clusters corresponding to a few identifiable cell types - neurons, stem cells, sensory cells, and epidermis. They use select marker genes from these clusters and perform HCR in situs on Loligo heads to describe these cell types. Their in situs describe a region similar to the lateral lip seen in other cephalopods where neural progenitors are found and from where neurons migrate into the brain.

Reviewer #1 (Public Review):

This work aimed at investigating how a BMI decoding performance is impacted by changing the conditions under which a motor task is performed. They recorded motor cortical activity using multielectrode arrays in two monkeys executing a finger flexion and extension task in four conditions: normal (no load, neutral wrist position), loaded (manipulandum attached to springs or rubber bands to resist flexion), wrist (no load, flexed wrist position) or both (loaded and flexed wrist). They found, as expected, that BMI decoders trained and tested on data sets collected during the same conditions performed better at predicting kinematics and muscle activity than others trained and tested across conditions. They also report that the performance of monkeys a BMI task involving the online control of a virtual hand was almost unaffected by changing either the actual manipulandum conditions as above or switching between decoders trained from data collected under different conditions. As for the neuronal activity, they found a mix of changes across task contexts. Interestingly, a principal component analysis revealed that activity in each context falls within well-aligned manifolds, and that the context-dependent variance in neuronal activity strongly correlated to amplitude of muscle activity.

Strengths:

The current study expands on previous findings about BMI decoders generalizability and contributes scientifically in at least three important ways.

First, their results are obtained from monkeys performing a fine finger control task with up to two degrees of freedom. This provides a powerful setting to investigate fine motor control of the hand in primates. The authors use the accuracy of BMI decoders between data sets as a measure of stationarity in the neurons-to-fingers mapping, which provides a reliable assessment. They show that changes in wrist angle or finger load affect the relationship between cortical neurons and otherwise identical movements. Interestingly, this result hold up for both kinematics and muscle activity predictions, albeit being stronger for the latter.

Second, their results confirming that neuronal activity recorded during different task conditions lies effectively within a common manifold is interesting. It supports prior observations, but in the specific context of finger movements.

Third, the dPCA results provide interesting and perhaps unexpected information about the fact that amplitude of muscle activity (or force) is clearly present in the motor cortical activity. This is possibly one of the most interesting findings because extracting a component from neural activity that can related robustly to muscle activity across context would provide great benefits to the development of BMIs for functional electrical stimulation.<br /> Overall, the analyses are well designed and the interpretation of the results is sound.

Weaknesses:

I found the discussion about the possible reasons why offline decoders are more sensitive to context than online decoders very interesting. Nonetheless, as the authors recognize, the possibility that the BMI itself causes a change in context, "in the plant", limits their interpretation. It could mean for the monkeys to switch from one suboptimal decoder to another, causing a ceiling effect occluding generalization errors.

Overall, several new and original results were obtained through these experiments and analyses. Nonetheless, I found it difficult to extract a clear unique and strong take-home message. The study comes short of proposing a new way to improve BMIs generalizability or precisely identifying factors that influence decoders generalizability.

Reviewer #1 (Public Review):

Chakrabarti et al study inner hair cell synapses using electron tomography of tissue rapidly frozen after optogenetic stimulation. Surprisingly, they find a nearly complete absence of docked vesicles at rest and after stimulation, but upon stimulation vesicles rapidly associate with the ribbon. Interestingly, no changes in vesicle size were found along or near the ribbon. This would have indicated a process of compound fusion prior to plasma membrane fusion, as proposed for retinal bipolar cell ribbons. This lack of compound fusion is used to argue against MVR at the IHC synapse. However, that is only one form of MVR. Another form, coordinated and rapid fusion of multiple docked vesicles at the bottom of the ribbon, is not ruled out. Therefore, I agree that the data set provides good evidence for rapid replenishment of the ribbon-associated vesicles, but I do not find the evidence against MVR convincing. The work provides fundamental insight into the mechanisms of sensory synapses.

Reviewer #1 (Public Review):

Hu et al. present findings that extend the understanding of the cellular and synaptic basis of fast network oscillations in the sensory cortex. They developed the ex vivo model system to study synaptic mechanisms of ultrafast (>400Hz) network oscillation ("ripplets") elicited in layer 4 (L4) of the barrel cortex in the mouse brain slice by optogenetically activating thalamocortical axon terminals at L4, which mimic the thalamic transmission of somatosensory information to the cortex. This model allowed them to reproduce extracellular ripplet oscillations in the slice preparation and investigate the temporal relationship of cellular and synaptic response in fast-spiking (FS) inhibitory interneurons and regular spiking (RS) with extracellular ripplet oscillations to common excitatory inputs at these cells. FS cells show precisely timed firing of spike bursts at ripplet frequency, and these spikes are highly synchronized with neighboring FS cells. Moreover, the phase-locked temporal relationship between the ripplets and responses of FS and RS cells, although different phases, to thalamocortical activation are found to closely coincide with EPSCs in RS cells, which suggests that common excitatory inputs to FS and RS cells and their synaptic connectivity are essential to generate reverberating network activity as ripplet oscillations. Additionally, they show that spikes of FS cells in layer 5 (L5) reduced in the slice with a cut between L4 and L5, proposing that recurrent excitation from L4 excitatory cells induced by thalamocortical optogenetic stimulation is necessary to drive FS spike bursts in layer 5 (L5).

Overall, this study helps extend our knowledge of the synaptic mechanisms of ultrafast oscillations in the sensory cortex. However, it would have been nice if the authors had utilized various methodologies and systems.

Although the overall findings are interesting, the conclusion of the study could have been strengthened according to the following points:

1. The authors investigate the temporal relationship between ripplets and FS and RS cells' response elicited by optogenetic activation of TC axon terminals, which is mainly supported by phase-locked responses of FS and RS cells with local ripplets oscillations to optogenetic activation. They also show highly synchronized FS-FS firing by eliminating electrical gap-junction and inhibitory synaptic connections to this synchrony. Based on these findings, the authors suggest that common excitatory inputs to FS and RS cells in L4 would be essential to generate these local ripplets. However, it interferes with the ability to follow the logical flow for biding other findings of phase-locking responses of FS and RS cells in ripplet oscillations in L4.

2. The authors suggest that the optogenetic activation of TC axon terminal elicits local ripplet oscillations via synchronized spike burst of FS inhibitory interneurons and alternating EPSC-IPSC of RS cells in phase-locked with ripplets in L4 barrel cortex, which would be generated by following common excitatory inputs from the local circuits to these cells at the ripple frequency. Thus they intend to investigate the source of these excitatory inputs at this local network of L4 by suppressing the firing of L4 RS cells. However, they show FS spike bursts in L5B, instead of L4, due to the technical limitations of their experimental setup, as described in the manuscript. Although L5 FS spike bursts decrease after cutting the L4/L5 boundary, supposedly inhibiting excitatory input from L4 as depicted in Fig 6D in the author's manuscript, the interpretation of data seems overly extended because it does not necessarily represent cellular and synaptic activities which are phase-locked with the ripplets observed in L4.

3. Authors suggested a circuit model. It would be recommended that the authors try to perform in silico analysis using the suggested model to explore the function of thalamocortical axons on the fast-spiking and regular-spiking neurons to support their circuit model.

Reviewer #1 (Public Review):

The authors had previously developed a method of determining conformational free energy differences between the alternative DFG-in and DFG-out conformational states of kinases using an energy function based on a Potts model. They did this because direct estimates of this free energy change from molecular simulations, while possible in principle, would in practice be hard to do with sufficient accuracy to be useful for such a large conformational transition. Potts model energies have been shown to be correlated with overall protein stability, so it is reasonable that dividing the contacts into DFG-in and DFG-out sets should allow the estimation of a free energy difference between conformational states. In this work they examine the differences between Tyrosine Kinases (TKs) and Serine/Threonine Kinases (STKs) more closely, finding that the model predicts a small free energy change for converting DFG-in to DFG-out for TKs but a significant unfavorable free energy cost to converting to DFG-out for the STKs. The most insightful part of the paper comes in its analysis of how this conformational change may contribute to the overall binding free energies. Calculating binding free energies for Type II inhibitors (which bind DFG-out) by alchemical methods neglects the contribution from any unfavorable conformational change ("reorganization energy") required to adopt the DFG-out conformation. Thus comparing this calculated binding free energy with the total binding free energy estimated from experiment allows an estimate of the conformational reorganization energy. It is found that this estimate is nicely correlated with the free energy change for conformational rearrangement estimated from the Potts model analysis. Thus an important contribution to Type II inhibitor binding is this conformational transition. The different contributions to Type II binding are analyzed in detail by further dissecting the Potts model.

Reviewer #1 (Public Review):

This paper investigates potential mechanisms underlying the generation of hippocampal theta and gamma rhythms using a combination of several modeling approaches. The authors perform new simulation experiments on the existing large-scale biophysical network model previously published by Bezaire et al. Guided by their analysis of this detailed model, they also develop a strongly reduced, rate-based network model, which allows them to run a much larger number of simulations and systematically explore the effects of varying several key parameters. The combined results from these two in silico approaches allow them to predict which cell types and connections in the hippocampus might be involved in the generation and coupling of theta and gamma oscillations.

In my view, several aspects of the general methodology are exemplary. In the current work as well as several earlier papers, the authors are re-using a large-scale network model that was originally developed in a different laboratory (Bezaire et al., 2016) and that still represents the state-of-the-art in detailed hippocampal modeling. Such model reuse is quite rare in computational neuroscience, which is rather unfortunate given the amount of time and effort required to build and share such a complex model. Very often, and also, in this case, the original publication that describes a detailed model provides only limited validation and analysis of model behavior, and the re-use of the same model in later studies represents a great opportunity to further examine and validate the model.

Combining detailed and simplified models can also be a powerful approach, especially when the correspondence between the two is carefully established. Matching results from the two models, in this case, allow strong arguments about key mechanisms of biological phenomena, where the simplified model allows the identification and characterization of necessary and sufficient components, while the detailed model can firmly anchor the models and their predictions to experimental data.

On the other hand, I have several major concerns about the implementation of these approaches and the interpretation of the results in the current study. First of all, the detailed model of Bezaire et al. is considered strictly equivalent, in all of its relevant details, to biological reality, and no attempt is made to verify or even discuss the validity of this assumption, even when particular details of the model are apparently critical for the results presented. I see this as a fundamental limitation of the current work - the fact that the Bezaire et al. model is the best one we have at the moment does not automatically make it correct in all its details, and features of the model that are essential for the new results certainly deserve careful scrutiny (preferably via detailed comparison with experimental data).

An important case in point is the strength of the interactions between specific neuronal populations. This is represented by different quantities in the detailed and simplified model, but the starting point is always the synaptic weight (conductance) values given by Bezaire et al. (2016), also listed in Tables 2 and 3 of the current manuscript. Looking at these parameters, one can identify a handful of connections whose conductance values are much higher than those of the other connections, and also more than an order of magnitude higher (50-100 nS) than commonly estimated values for cortical synapses (normally less than about 5 nS, except for a few very special types of synapse such as the hippocampal mossy fibers). Not surprisingly, several of these connections (such as the pyramidal cell to pyramidal cell connections, and the CCK+BC to PV+BC connections) were found to be critical for the generation and control of theta and gamma oscillations in the model. Given their importance for the conclusions of the paper, it would be essential to double-check the validity of these parameter values. In this context, it is worth noting that, unlike the anatomical parameters (cell numbers and connectivity) that had been carefully calculated and discussed in Bezaire and Soltesz (2013), biophysical parameters (the densities of neuronal membrane conductances and synaptic conductances) in Bezaire et al. (2016) were obtained by relatively simple (partly manual) fitting procedures whose reliability and robustness are mostly unknown. Specifically for synaptic parameters in CA1, a more systematic review and calculation were recently carried out by Ecker et al. (2020); their estimates for the synaptic conductances in question are typically much lower than those of Bezaire et al. (2016) and appear to be more in line with widely accepted values for cortical (hippocampal) synapses.

Furthermore, some key details concerning the construction of the simplified rate model are unclear in the current manuscript. The process of selecting cell types and connections for inclusion in the rate model is described, and the criteria are mostly clear, although the results are likely to be heavily affected by the problems discussed above, and I do not understand why the strength of external input was included among the selection criteria for cell types (especially if the model is meant to capture the internal dynamics of the isolated CA1 region). However, the main issue is that it remains unclear how the parameters of the rate model (the 24 parameters in Table 4) were obtained. The authors simply state that they "found a set of parameters that give rise to theta-gamma rhythms," and no further explanation is provided. Ideally, the parameters of the rate model should be derived systematically from the detailed biophysical model so that the two models are linked as strongly as possible; but even if this was not the case, the methods used to set these parameters should be described in detail.

An important inaccuracy in the presentation of the results concerns the suggested coupling of theta and gamma oscillations in the models. Although the authors show that theta and gamma oscillations can be simultaneously present in the network under certain conditions, actual coupling of the two rhythms (e.g., in the form of phase-amplitude coupling) is not systematically characterized, and it is therefore not clear under what conditions real coupling is present in the two models (although a probable example can be seen in Figure 1C(ii)).

The Discussion of the paper states that gamma oscillations in the model(s) are generated via a pure interneuronal (ING) mechanism. This is an interesting claim; however, I could not find any findings in the Results section that directly support this conclusion.

Finally, although the authors write that they can "envisage designing experiments to directly test predictions" from their modeling work, no such experimental predictions are explicitly identified in the current manuscript.

Reviewer #1 (Public Review):

This study focuses on the role of polo like kinase 1 (PLK-1) during oocyte meiosis. In mammalian oocytes, Plk1 localizes to chromosomes and spindle poles, and there is evidence that it is required for nuclear envelope breakdown, spindle formation, chromosome segregation, and polar body extrusion. However, how Plk1 is targeted to its various locations and how it performs these functions is not well understood. This study uses C. elegans oocytes as a model to explore PLK-1 function during meiosis. They take advantage of an analogue-sensitive allele of plk-1, which enabled them to bypass nuclear envelope breakdown defects that occur following PLK-1 RNAi. This allowed them to dissect later roles of PLK-1 in oocytes, demonstrating that depletion causes defects in spindle organization, chromosome congression, segregation, and polar body extrusion. Moreover, the authors defined mechanisms by which PLK-1 is targeted to chromosomes, showing that CENP-C (HCP-4) is required for localization to chromosome arms and that BUB-1 is required for targeting to the midbivalent region. Finally, they demonstrate that upon removal of PLK-1 from both domains, there are severe meiotic defects. These findings are interesting. However, there is a need for additional analysis to better support some of their conclusions, and to aid in interpretation of particular phenotypes. Specific comments are below.

- For many important claims of the paper, a single representative image is shown but the n is not noted. This is an issue throughout the paper for much of the localization analysis (e.g. Figure 1B, 1C, 1D, 2A, 2B, 3A, 3B, 3C, etc.); in cases like this, numbers should be included to increase the rigor of the presented data. How many images or movies were analyzed that looked like the one shown? For linescans, were they done only on one image? How many independent experiments were done, etc?.

- In the abstract, it is stated that PLK-1 plays a role in spindle assembly/stability (this is also stated elsewhere, e.g. line 101). This phrasing implies that the authors have demonstrated roles in both spindle assembly and stability. However, to distinguish between these roles, they would have to show that removal of PLK-1 before spindle assembly causes defects, and also that removal of PLK-1 from pre-formed spindles causes collapse. I don't think it is necessary to do this, as the spindle roles of PLK-1 are not a focus of the paper. However, the language should be altered so that it does not imply that the paper has demonstrated roles in both. A good place to do this would be in the section from lines 144-147, where they first discuss the spindle defects. It would be straightforward to explain that their approach does not distinguish between spindle assembly and stability, and that PLK-1 could have a role in either or both.

- It is stated that there is kinetochore localization of PLK-1 (and I do see some dim cup-like localization in images after PLK-1 is removed from the chromosome arms via HCP-4 RNAi). However, this cup-like localization is not clear in most wild-type images (e.g. Figure 1B, 1D, 2A, 3A, etc.). Although I recognize that the chromatin staining might be obscuring kinetochore localization, if PLK-1 was truly a kinetochore protein I would also expect it to localize to filaments within the spindle (as many other kinetochore proteins do), especially since the authors state that BUB-1 targets PLK-1 to the kinetochore (and BUB-1 is in the filaments). In fact, the only images where it looks like PLK-1 may be localized to filaments are in Figure 4C and 6A, when HCP-4 has been depleted (though I don't know if this generally true across all HCP-4 RNAi images). For me, this calls into question the conclusion that PLK-1 truly is on the kinetochore in wild type conditions - could it be that PLK-1 only localizes to the kinetochore (and to the filaments) when HCP-4 is depleted? The authors need to resolve this issue and provide better evidence that PLK-1 normally localizes to the kinetochore, if they want to make this claim. Additionally, the observation that PLK-1 is not on the kinetochore filaments (in wild type conditions) should be addressed in the text somewhere - do the authors think that this is a special type of kinetochore protein that does not localize to the filaments?

- The authors should provide a control experiment, treating wild-type worms with 10uM 3-IB-PP1. This would be important to ensure that the spindle defects seen at this concentration in the plk-1as strain are not non-specific effects of the inhibitor. There is a control in Figure 1 - figure supplement 3 using 1uM 3-IB-PP1 but didn't see a control for 10uM (the concentration at which spindle defects are observed).

- In Figure 2F, the gels for BUB-1+PLK-1 look different in the presence and absence of phosphorylation by Cdk1 - for these data, I agree with the authors that it looks as if the complex elutes at a higher volume if BUB-1 is not phosphorylated (lines 200-204). However, Figure 2G has a repeat of the condition with phosphorylated BUB-1, and in this panel, the complex appears to elute at a higher volume than it did on the gel in panel F. The gel in panel G looks much more similar to the unphosphorylated condition in panel F. The authors need to explain this discrepancy (i.e., Is there a reason why the gels cannot be compared between panels? How reproducible are these data?). Ideally, the authors would include a repeat of the unphosphorylated BUB-1 + PLK-1 condition in panel G, done at the same time as the conditions shown in that panel, to avoid the impression that their results may not be reproducible.

- The authors would need to provide convincing evidence that co-depletion of BUB-1 and HCP-4 delocalizes PLK-1 from the chromosomes entirely, and that this co-depletion condition is more severe than either single depletion alone. Additionally, the bub-1T527A and hcp-4T163A alleles are nice tools to, in theory, more specifically delocalize PLK-1 from the midbivalent and chromosome arms, respectively, to explore the functions of chromosome-associated PLK-1. However, I think the authors cannot rule out the possibility that other proteins are also being depleted from the midbivalent and/or chromosome arms in their conditions, and that this delocalization may contribute to the phenotypes observed. For example, hcp-4 depletion was recently shown to delocalize KLP-19 from the chromosome arms (Horton et.al. 2022), so in the experiment shown in Figure 6E (HCP-4 RNAi in the bub-1 mutant), PLK-1 was likely not the only protein missing from the chromosome arms. Therefore, understanding if other proteins are absent from these domains (in the bub-1T527A and hcp-4T16A3 mutants) would help the reader understand and interpret the presented phenotypes (and how specific they are to PLK-1 loss). Consequently, I think that to better understand the co-depletion analysis presented in Figure 6 (and Figure 6 supplement 1), the authors should analyze other midbivalent and chromosome arm proteins, to determine if any are also delocalized (e.g. SUMO, KLP-19, MCAK, etc.). Additionally, instead of performing a combination of mutant and RNAi analysis (i.e. HCP-4 RNAi in the bub-1 mutant (Figure 6) and BUB-1 RNAi in the hcp-4 mutant (Figure 6 figure supplement 1)), it would be more powerful to generate a double mutant - this has a higher chance of being a more specific depletion condition.

Reviewer #1 (Public Review):

This manuscript presents information that will be of great interest to yeast geneticists - standard gene deletions can lead to misleading phenotypes due to effects on adjacent genes. The experiments carefully document this in one case, for the DBP1 gene, and present additional evidence that it can occur at additional genes. An improved version of the standard gene replacement cassette is described, with evidence that it functions in an improved fashion, insulated from affecting adjacent genes.

Reviewer #1 (Public Review):

This study analyzes the detailed chemical mechanics of the formation of a physiologically important protein multimer. The primary strengths of the study are careful analyses of two distinct methods, CG-MALS a direct measure of multimerization, and environment-sensitive tryptophan fluorescence, that each indicates that Ca2+ activation of the C-lobe alone can change the physical interaction with an SK2 C-terminal peptide. An intriguing finding is that while either the N- or C-lobes alone can interact with the C-terminal peptide, only with full-length CaM can the SK C-terminal peptide be bound by two CaM molecules simultaneously. This study also clearly demonstrates that Ca2+ activation of the N-lobe triggers binding to the SK2 C-terminal peptide. Methods descriptions are thorough and excellent. Discussion of relevance to structures and function are nuanced and free of presumptions. The weaknesses of this manuscript are that the physiological implications of these findings are not clear: CaM interacts with regions of SK channels besides the C-terminal peptide studied here, and no evidence is provided here that C-lobe calcium binding alters channel opening. Overall, the evidence for conformational changes of the complex due to Ca2+ binding to the C-lobe alone is very strong, and physiological importance seems likely. The interpretation of data in this manuscript is mostly cautious and logically crystalline, with alternative interpretations discussed at many junctures.

Reviewer #1 (Public Review):

The overarching hypothesis is that cadherin adhesion molecules specify the code that enables the premotor brainstem breathing circuits to innervate the phrenic motor neurons that control the primary breathing muscle, the diaphragm. The authors show that multiple type 1 and 2 cadherins (N-, 6, 9, 10) are expressed by phrenic motor neurons and are necessary for motor neuron development and breathing, and complementarily, that adhesion signaling in medullary breathing circuits are required for normal breathing. The presented data support a model whereby combinations of redundant adhesion molecules create a code to wire the breathing circuit.

Strengths:<br /> 1) The authors first use a complex, rigorous genetic approach to eliminate N, 6, 9, 10 cadherins from motor neurons and discover using whole body plethysmography that neonates do not breath.<br /> 2) Then, the authors provide a thorough description of the anatomy of the mutant motor neurons and discover that the number of motor neurons decreases, the soma anatomical positions and dendritic arborization shift, and there is decreased innervation of the diaphragm breathing muscle.<br /> 3) That Cdh9 medullary expressing neurons are premotor to Cdh9 expressing phrenic motor neurons.<br /> 4) Cadherin signaling is required for normal breathing.

Weaknesses: The main conclusion that ablation of the cadherin code decreases synaptic connectivity between the rVRG and phrenic motor neurons is never directly shown. This can only be inferred by the data.<br /> 1) Conclusion that the connectivity between rVRG premotor and phrenic nerve motor neurons is "weaker". This conclusion is inferred from several experiments but is never directly demonstrated. Alternative interpretations of the decreased amplitude of the in vitro phrenic nerve burst is that the rootlet contains fewer axons (as predicted by the fewer motor neurons in S3 and innervation of the diaphragm S2). Additionally, the intrinsic electrophysiological properties of the motor neurons might be different. To show this decisively, the authors could use electrophysiological recordings of phrenic motor neurons to directly measure a change in synaptic input (for example, mEPSPs or EPSPs after optogenetic stimulation of rVRG axon terminals). Without a direct measurement, the synaptic connectivity can only be inferred.<br /> 2) Conclusion that the small phenic nerve burst size in Dbx1 deleted cadherin signaling is due to less synaptic input to the motor neurons. Dbx1 is expressed in multiple compartments of the medullary breathing control circuit, like the breathing rhythm generator (preBötC). The smaller burst size could be due to altered activity between preBötC neurons to create a full burst, the transmission of this burst from the preBötC to the rVRG, etc.<br /> 3) In vitro burst size. The authors use 4 bursts from each animal to calculate the average burst size. How were the bursts chosen? Why did the authors use so few bursts? What is the variability of burst size within each animal? What parameters are used to define a burst? This analysis and the level of detail in the figure legend/methods section is inadequate to rigorously establish the conclusion that burst size is altered in the various genotypes.<br /> 4) The authors state that the in vitro frequency in figure 4 is inaccurate, but then the in vitro frequency is used to claim the preBötC is not impacted in Dbx1 mutants (conclusion section "respiratory motor circuit anatomy and assembly"). To directly assess this conclusion, the bursting frequency of the in vitro preBötC rhythm should be measured.<br /> 5) The burst size in picrotoxin/strychnine is used to conclude that the motor neurons intrinsic physiology is not impacted. The bursts are described, and examples are shown, but this is never quantified across many bursts within in a single recordings nor in multiple animals of each genotype.

Reviewer #1 (Public Review):

This paper describes the results of a MEG study where participants listened to classical MIDI music. The authors then use lagged linear regression (with 5-fold cross-validation) to predict the response of the MEG signal using (1) note onsets (2) several additional acoustic features (3) a measure of note surprise computed from one of several models. The authors find that the surprise regressors predict additional variance above and beyond that already predicted by the other note onset and acoustic features (the "baseline" model), which serves as a replication of a recent study by Di Liberto.

They compute note surprisal using four models (1) a hand-crafted Bayesian model designed to reflect some of the dominant statistical properties of Western music (Temperley) (2) an n-gram model trained on one musical piece (IDyOM stm) (3) an n-gram model trained on a much larger corpus (IDyOM ltm) (4) a transformer DNN trained on a mix of polyphonic and monophonic music (MT). For each model, they train the model using varying amounts of context.

They find that the transformer model (MT) and long-term n-gram model (IDyOM stm) give the best neural prediction accuracy, both of which give ~3% improvement in predicted correlation values relative to their baseline model. In addition, they find that for all models, the prediction scores are maximal for contexts of ~2-7 notes. These neural results do not appear to reflect the overall accuracy of the models tested since the short-term n-gram model outperforms the long-term n-gram model and the music transformer's accuracy improves substantially with additional context beyond 7 notes. The authors replicate all these findings in a separate EEG experiment from the Di Liberto paper.

Overall, this is a clean, nicely-conducted study. However, the conclusions do not follow from the results for two main reasons:

1. Different features of natural stimuli are almost always correlated with each other to some extent, and as a consequence, a feature (e.g., surprise) can predict the neural response even if it doesn't drive that response. The standard approach to dealing with this problem, taken here, is to test if a feature improves the prediction accuracy of a model above and beyond that of a baseline model (using cross-validation to avoid over-fitting). If the feature improves prediction accuracy, then one can conclude that the feature contributes additional, unique variance. However, there are two key problems: (1) the space of possible features to control for is vast, and there will almost always be uncontrolled-for features (2) the relationship between the relevant control features and the neural response could be nonlinear. As a consequence, if some new feature (here surprise) contributes a little bit of additional variance, this could easily reflect additional un-controlled features or some nonlinear relationship that was not captured by the linear model. This problem becomes more acute the smaller the effect size since even a small inaccuracy in the control model could explain the resulting finding. This problem is not specific to this study but is a problem nonetheless.

2. The authors make a distinction between "Gestalt-like principles" and "statistical learning" but they never define was is meant by this distinction. The Temperley model encodes a variety of important statistics of Western music, including statistics such as keys that are unlikely to reflect generic Gestalt principles. The Temperley model builds in some additional structure such as the notion of a key, which the n-gram and transformer models must learn from scratch. In general, the models being compared differ in so many ways that it is hard to conclude much about what is driving the observed differences in prediction accuracy, particularly given the small effect sizes. The context manipulation is more controlled, and the fact that neural prediction accuracy dissociates from the model performance is potentially interesting. However, I am not confident that the authors have a good neural index of surprise for the reasons described above, and this limits the conclusions that can be drawn from this manipulation.

Reviewer #1 (Public Review):

The authors used data from extracellular recordings in mouse piriform cortex (PCx) by Bolding & Franks (2018), they examined the strength, timing, and coherence of gamma oscillations with respiration in awake mice. During "spontaneous" activity (i.e. without odor or light stimulation), they observed a large peak in gamma that was driven by respiration and aligned with the spiking of FBIs. TeLC, which blocks synaptic output from principal cells onto other principal cells and FBIs, abolishes gamma. Beta oscillations are evoked while gamma oscillations are induced. Odors strongly affect beta in PCx but have minimal (duration but not amplitude) effects on gamma. Unlike gamma, strong, odor-evoked beta oscillations are observed in TeLC. Using PCA, the authors found a small subset of neurons that conveyed most of the information about the odor (winner cells). Loser cells were more phase-locked to gamma, which matched the time course of inhibition. Odor decoding accuracy closely follows the time course of gamma power.

I think this is an interesting study that uses a publicly available dataset to good effect and advances the field elegantly, especially by selectively analyzing activity in identified principal neurons versus inhibitory interneurons, and by making use of defined circuit perturbations to causally test some of their hypotheses.

Major:

- The authors show odor-specificity at the time of the gamma peak and imply that the gamma coupling is important for odor coding. Is this because gamma oscillations are important or because gamma is strongest when activity in PCx is strongest (i.e. both excitatory and inhibitory activity, which would cancel each other in the population PSTH, which peaks earlier)? To make this claim, the authors could show that odor decoding accuracy - with a small (~10 ms sliding window) - oscillates at approx. gamma frequencies. As is, Fig. 5 just shows that cells respond at slightly different times in the sniff cycle. What time window was used for computing the Odor Specificity Index? Put another way, is it meaningful that decoding is most accurate when gamma oscillations are strongest, or is this just a reflection of total population activity, i.e., when activity is greatest there is more gamma power, and odor decoding accuracy is best?

- The authors say, "assembly recruitment would depend on excitatory-excitatory interactions among winner cells occurring simultaneously during gamma activity." Can the authors test this prediction by examining the TeLC recordings, in which excitatory-excitatory connections are abolished?

- The authors show that gamma oscillations are abolished in the TeLC condition and use this to claim that gamma arises in the PCx. However, PCx neurons also project back to the OB, where they form excitatory connections onto granule cells. Fukunaga et al (2012) showed that granule cells are essential for generating gamma oscillations in the bulb. Can the authors be sure that gamma is generated in the PCx, per se, rather than generated in the bulb by centrifugal inputs from the PCx, and then inherited from the bulb by the PCx?

Reviewer #1 (Public Review):

In this manuscript, the authors built logistic regression prediction models for linear growth faltering using demographic, socioeconomic, and clinical variables, with the objective of developing a clinical prediction rule that could be applied by healthcare workers to identify and treat high-risk children. A model with 2 variables selected by random forest variable importance performed similarly to a model with 10 variables. Age and HAZ at baseline were selected for the 2-variable model, consistent with existing literature. The authors externally validated the 2-variable model and found similar discriminative ability. Based on typical rule-of-thumb cutoffs, model performance was moderate (AUCs of ~0.65-0.75, depending on model specification); models may still be useful in practice, but this should be further discussed by the authors.

Strengths:

Linear growth faltering is a pressing issue with broad, negative impacts on the health, development, and well-being of children worldwide. In this work, the authors applied clearly explained, thoughtful approaches to variable selection, model specification, and model validation, with large, multi-country cohorts used for training and external validation. Appropriate datasets for external validation can be challenging to find, but the MAL-ED data used here is well-suited to the task, with similar predictor and outcome measurements to the GEMS training data. The well-characterized studies allowed the authors to explore a wide range of potential predictors for stunting, including socioeconomic factors, antibiotic use, and diarrheal etiology.

Weaknesses:

This work would benefit from additional discussion around the clinical relevance of the results. For example, what is the current standard of care for prevention of stunting, and how much would this model improve the status quo? Is specificity of 0.47 in the context of sensitivity of 0.80 an acceptable tradeoff with regards to the interventions that would be used? More discussion around these points is necessary to support the authors' conclusions that these models could potentially be used to support clinical decisions and target resources.

In addition to the external validation, further investigation of model performance in key subpopulations would strengthen the importance and applicability of the work. For example, performance of prediction models may vary widely by setting; it would be valuable to show that the model has similar performance in each country. Another key sensitivity analysis would be to show consistent model performance by HAZ at baseline. The authors note that stunting may be challenging to reverse (p.20), and many of the children are already below the typical cutoff of HAZ<-2 at baseline; it would be valuable to show model performance among the subgroup of children for whom treatment would be most beneficial.

Reviewer #1 (Public Review):

The current study uses microbiology, biochemistry, microscopy, and viral vectors to establish a role for prefrontal cortex expression of the immediate early gene NPAS4 in sucrose preference and dendritic spine morphology in the mouse social defeat stress model. The experimental designs are appropriate and the hypotheses addressed are interesting. The paper is generally very well-written and the figures are clear. Most of the statistical analyses are appropriate, and they are reported in clear and useful tables. Thus, the general potential for the studies is quite high. The authors conclusively show that NPAS4 is induced in mPFC in response to social defeat stress and that NPAS4 is important for stress-induced changes in mPFC dendritic spine number. However, some of the key data regarding reward motivation are difficult to properly interpret and do not convincingly demonstrate a behavioral result of NPAS4 knockdown in mPFC. Moreover, the spine morphology and sequencing analyses lack depth. Most importantly, although the authors explore the effects of reducing NPAS4 expression in mPFC, they do not explore the effects of increasing NPAS4 expression or function, and thus the studies seem incomplete and cannot be fully interpreted.

Reviewer #1 (Public Review):

The study by Osei-Owusu and colleagues addresses the mechanism of desensitization of the proton-activated chloride (PAC) channel. In three recent milestone papers, the authors have cloned the channel, identified its cryo-EM structure under high-pH and low-pH conditions, and addressed the mechanism of its pH-dependent activation. Interestingly, despite dramatic rearrangements in the TM domain, both the high- and the low-pH structures showed a closed pore, suggesting that the latter might represent an inactivated state. In the current study, the authors show that prolonged exposure of PAC to an acidic extracellular solution causes inactivation which is rapidly reversible at high pH. They further show that four mutations (H98R, E107R, D109R, H250R) that are predicted to disrupt interactions that stabilize the low-pH structure reduce PAC inactivation. On the other hand, two mutations that accelerate inactivation (D91R, E94R) are predicted to stabilize the low-pH structure based on MD simulations. The work thus functionally supports the earlier hypothesis that the low-pH cryo-EM structure indeed represents an inactivated state. Moreover, it identifies several key titratable residues that are involved in this process.

The choice of the tested residues is based on strong structural evidence, and the electrophysiological data largely seem to support the conclusions, even though the analysis is not always rigorous. (Time constants seem unreliable as they are extracted from decay time courses that are too short to be reliably fitted, but comparisons of the simple parameter "fractional surviving current after 30 s" seem convincing enough.) Some of the mechanistic conclusions are largely based on MD simulations which I am not qualified to assess.

Reviewer #1 (Public Review):

The authors suggest that there is a long-term periodicity of individual antibody response to influenza A (H3N2). The interesting periodicity may be surely appeared. Though the authors assume that the periodicity is driven by pre-existing antibody responses, the authors could provide more supportive data and discuss some possibilities.

1. The authors can investigate whether the periodicity reflects an epidemic/invasion record of A(H2N3) within Guangzhou or the surrounding city, e.g., the numbers of flu-infected people yearly can be referred to.

2. The authors can consider whether the participants are recently/previously vaccinated and/or infected with flu. The remaining antibodies may reflect a long memory but may show a recent activation.

3. The strains inducing high HI titers may have similar mutations and may be reactive to the same antibodies. What are the mutation frequencies among 21 A(H3N2) strains?

Reviewer #1 (Public Review):

While the circuits underlying the computation of directional motion information in the fly brain are very well described, much less is known about the neurons serving the detection of objects. In a previous publication from the same lab, it has been shown that flies perform body saccades to track a moving object during flight. In the current paper, Frighetto and Frye provide evidence that T3 cells, a population of neurons within the optic lobes, are involved in this task. First, they performed 2-photon Calcium imaging from T3 cells to show that these cells respond to moving bars, which they later use in behavioural experiments. They then silenced T3 cells using genetic tools and tested the behavior of these flies in response to a rotating bar using two different setups. In one, the flies are fixed and bilateral changes in wing stroke amplitude are used as a measure for turning, in the other, flies are magnetically tethered such that they can rotate around the vertical body axis. Silencing T3 cells leads to the abolishment of the steering response induced by object position using a bar that is defined by its motion relative to the surround, but leaves the response to object motion intact. In the magnetically tethered flies, it reduces the number of saccades and thus leads to an impairment of bar-tracking behavior. In another set of experiments they optogenetically activated the whole population of T3 neurons (which supposedly impairs their normal function), which leads to an increase in the number of saccades after the activation (when the light stimulus used to activate the cells is turned off). Silencing the neurons necessary for detection of local motion, T4 and T5 cells, in contrast reduces responses elicited by object motion rather than position, but also has an impact on object tracking saccades. The authors provide a simple model, where speed-dependent signals from multiple T3 cells are integrated and trigger a saccade, when a threshold is reached.

The data generally support the conclusion that T3 cells play a role in detecting bar position and in controlling saccades in response to rotating bars. However, there are some inconsistencies in the data that are not sufficiently explored and discussed.

1. In a previous paper from the lab (Keleş et al., 2020), it was shown that T3 cells respond preferentially to small objects, whereas here they robustly respond to elongated bars and even large-field gratings. This discrepancy is not discussed.

2. In a previous paper, the authors showed that integrated positional error rather than bar position is used to elicit bar-tracking saccades and that saccade amplitude is relatively stereotyped. However, here they show, that T3 cells respond much more strongly to a slowly moving stimulus (18{degree sign}/s) rather than to the fast moving stimuli used for the behavioral experiments (> 90{degree sign}/s). This response property plays an important role for the model they propose. My general concern here is that the findings might not be generalizable to slower moving bars, where more precise, position-dependent responses could play a larger role, and that these fast moving bar stimuli represent an extreme situation, where the flies cannot accurately track bar position any more.

3. The claim that T3 cells are tuned to stimulus velocity is not supported by the data in my view. For the bar stimuli, the authors only tested speeds of 18{degree sign}/s and above 90{degree sign}/s, but nothing in between. For the grating motion there seems to be an influence of temporal frequency for the same stimulus velocity (see e.g. Fig.1_1), but this is not quantified.

4. The results from the optogenetic activation experiments are hard to interpret, as it is unclear how a prolonged activation of all T3 cells would affect the downstream circuitry. It is not clear that this experiment is equivalent to a "loss-of-function perturbation" of T3 cells as the authors claim in the text.

Reviewer #1 (Public Review):

This work provides a new general framework for estimating missing data on cervical cancer epidemiology, including sexual behavior, HPV prevalence, and cervical cancer incidence. These data are useful to determine impact projections of cervical cancer prevention. The authors suggest a three-step approach: 1) a clustering method applied on registries with an intermediate level of data availability to cluster cervical cancer incidence based on a Poisson-regression-based CEM algorithm, 2) a classification method applied on registries with a low level of data availability to classify cervical cancer incidence based on a Random Forest, 3) a projection method applied on missing data based on the mean of available data. The authors use India as a case study to implement this new methodology. Results indicate that two patterns of cervical cancer incidence are identified in India (high and low incidence), classifying all Indian states with missing data to a low incidence. From this classification, missing data is approximated using the mean of the available data within each cluster.

A strength of this approach is that this methodology can be applied to regions with missing data, although a minimum set of information is needed. This makes it possible to have individual data for each unit in the region.

One of the weaknesses of this methodology is the need for a minimum set of epidemiological data to enable impact projections. It is true that when epidemiological cervical cancer data is not available, authors mentioned that general indicators (e.g., human development index, geography) can be used but projections will be probably less realistic. As observed with other techniques, countries with fewer resources have less data available and cannot benefit from these types of techniques to have more adequate guidelines.

Imputation of missing data is always a challenging issue. The technique proposed in this manuscript is an interesting new approach to missing data imputation that could be applied with a minimum set of available data. However, we must focus on obtaining reliable data from each region of the world to help local health authorities implement better preventive measures for the local population.

Reviewer #1 (Public Review):

The manuscript by Masschelin et al. describes how Vitamin B2 deficiency affects body composition, energy expenditure, and glucose metabolism. B2 deficient mice have lower O2 consumption, and locomotor activity, with no difference in food intake. These mice also have lower liver FAD levels, which is expected given that B2 is a necessary cofactor for this coenzyme. Additionally, these mice have lower blood glucose levels following pyruvate injection, implying a lower capacity for gluconeogenesis. Using PPAR KO mice, they show that this effect on pyruvate tolerance is due to PPARα activation, though there is still a minor difference between wild-type and KO mice. Importantly, they show that fenofibrate PPARagonism can improve glucose output following pyruvate injection in the absence of B2. The authors also perform robust metabolomics in each experimental condition and phenotype of the mouse well.

1. The authors have yet to explore other explanations of differences in glucose metabolism under B2D +/-Fenofibrate. The canonical targets of PPARα are involved in fatty acid oxidation, ketogenesis, and VLDL/HDL metabolism, in addition to gluconeogenesis (Bougarne et al. 2018). Gluconeogenesis is more of a fasting response due to CREB, FOXO1/PGC1activation rather than PPAR. In response to B2D, the PPARα KO mice have increased plasma TGs, which may suggest a difference in VLDL TG secretion (Suppl. S3). Perhaps lipid metabolism is more directly affected, and changes in glucose metabolism are secondary to that of triglyceride metabolism. Regarding ketogenesis, the fenofibrate+ B2D fed mice have decreased plasma beta-hydroxybutyrate, suggesting decreased ketogenesis, which is a more canonical PPARα pathway (Suppl. S3). Testing each of these processes would help control that this mechanism is specific to gluconeogenesis and not secondary to something else.

2. Is the effect on ISR dependent on PPARα? Is the mechanism of Fenofibrate on the liver, or on another cell type? In Figure 1, the authors state that Riboflavin deficiency alters body composition and energy expenditure, and then focuses on the liver. However, FAD levels are also increased in the heart and kidneys in addition to the liver. These tissues also respond to PPARα agonism, in addition to the muscle which plays a role in regulating glucose metabolism (B2D mice also have a higher lean mass (Fig 1e)). Additionally, the authors haven't shown specifically if the effects of fenofibrate on electron transport and the ISR are dependent on the presence of PPARα (Figure 5, 6).

Reviewer #1 (Public Review):

In their manuscript, Haenelt et al. investigated the structure-function relationship for cortical columns in the in vivo human brain. The example they used is the thick stripe - pale stripe - thin stripe organisation of secondary visual cortex (V2).

The specific strength of the current study lies in the combination of cutting edge imaging protocols for quantitative measurements of myelin-related signals (qMRI) together with functional activation, both at submillimeter resolution at high field (7T). This allowed the visualisation of the stripy organisation of V2 with regards to colour (thin stripes) and binocular disparity (thick stripes) as well as myelination in individual human subjects. The main results suggest higher myelination for the pale stripe regions. This is in line with some earlier studies, across primate species, but not with others.

One potential issue is that the high myelination signal is associated with the compartment in V2 (pale stripes) which was not functionally defined itself but by the absence of specific functional activations. No difference was reported between those stripes that were defined functionally. Other explanations for the differential pattern of a qMRI signals, e.g. ROI distribution for presumed pale stripes is not evenly distributed (more foveal), ROIs with low activations due to some other factor show higher myelin-related signals, cannot be excluded based on the analysis presented.

Another theoretical and practical issue is the question of "ground truth" for the non-invasive qMRI measures, as the authors - as their starting point - roundly dismiss direct histological tissue studies as conflicting, rather than take a critical look at the merit of the conflicting study results and provide a best hypothesis. If so, they need to explain better how they calibrate their non-invasive MR measurements of myelin.

While this paper makes an important contribution to the question of the association of specific myelination patterns defining the columnar architecture in V2, it is not entirely clear whether the authors can fully resolve it with the data presented.

The highly sophisticated methods and detailed analysis show that high resolution investigation of structure-function relationship of the columnar organisation in human visual cortex are feasible and reliable. V2 stripe patterns can be visualised structurally (with quantitative myelin-related measurements) and functionally (based on functional selectivity, which is of considerable importance for the field. The results indicate that in humans, the pale or inter strip regions might be associated with high patterns of myelination.

Reviewer #1 (Public Review):

This manuscript attempts to disclose new insights into barrel cortex cell class-dependent and cell depth-dependent membrane potential (Vm) dynamics during active whisker sensing. The results highlight similarities but also specific differences between different types of cortical neurons. The approach used is very effective and direct: somatosensory stimulation is performed in awake animals without anesthesia, the neurons are recorded with intracellular whole cell patch clamp recording that can provide fast responses with high resolution, and the identification of various neuron types is achieved by using mice expressing genetically defined selective fluorescent markers. The results support the main conclusions. The work is an extension of previous, similar work performed by this group, However, most previous Vm studies in the mouse barrel cortex during behavior have largely focused on superficial neurons located in the upper ~300 μm of the neocortex since these are more easily targeted through two-photon microscopy. In this study, the authors extend current knowledge by investigating Vm dynamics across a greater range of depths including two-photon targeted whole-cell recordings across the upper ~600 μm of the neocortex. I believe that this manuscript uses a demanding, but excellent approach that will be useful to other researchers in the field. The manuscript is likely to be influential.

Reviewer #1 (Public Review):

Proton Activated Chloride (PAC) channels have been recently identified as important contributors to endosomal acidification, and their activity in the plasma membrane increases under certain pathological conditions and can induce cellular death. There is very limited information on the pharmacology of these ion channels. By recording from endogenous PAC channels stimulated with an acidic extracellular solution in HEK293 cells using the patch-clamp technique, this study finds that PAC channels are inhibited physiological concentrations of the soluble short-chain PIP2 analog dic-8-PIP2. Inhibition is quantified for several PIP2-related lipids with different number of headgroup phosphates and shorter or longer acyl chains, and it is found that an acyl chain with more than 8 carbons and a negative headgroup charge are both required for robust inhibition. Importantly, inhibition appears to result from PIP2 incorporated into the outer membrane leaflet, as treatment of the inner leaflet with PIP2 or poly-lysine to either increase or decrease PIP2, respectively, did not have any effect on channel activity, as opposed to when the lipid is extracellularly applied. A structure of the channel in the presence of PIP2 was obtained using single-particle cryo-electron microscopy - the structure resembles a previously observed conformation for PAC channels that likely represents a non-conducting desensitized state, and it contains densities with a shape that is consistent with a bound PIP2 molecule in the outer leaflet. Mutations to alanine based on the channel-lipid interactions observed in the structure were all found to disrupt inhibition of PAC channels by PIP2, consistent with the location of the lipid binding site proposed in the study. By comparing the amino acid sequence of human PAC channels with those of other species, it is found that the proposed lipid binding site is highly conserved except in zebrafish. Notably, zebrafish PAC channels are less susceptible to inhibition by PIP2, and mutation of residues at the binding site to those present in the human channel increases inhibition, consistent with the proposed location of the binding site for PIP2. Finally, it is found that the kinetics of inhibition by PIP2 are positively correlated with the degree of channel activation and also with the kinetics of desensitization, suggesting that PIP2 binds more favorably to the desensitized state of the channel whereas it does not bind to the closed state, providing a possible mechanism for the inhibition.

Results are clearly reported and findings are generally robust. One concern is that most of the electrophysiological characterization of the inhibition of PAC channels by PIP2 lipids was done using endogenously expressed channels. It is unclear why this was done because mutant channels are studied in a PAC KO cell line that could have been used for all experiments. The effects of acidic pH and acidic pH + PIP2 in cells that do not express PAC channels is therefore not shown, but would be important to establish that the measured effects of the lipid are specific to PAC channels.

Another concern for the study is related to the uncertainty in establishing that the bound lipid is indeed PIP2. Although the mutagenesis results are all consistent with the proposed binding site, it remains a possibility that the mutations affect PIP2 inhibition indirectly by e.g. changing the rate of channel desensitization, which was not measured for any of the mutants on Figure 3E. There is not additional analysis performed to determine whether other types of lipids could occupy the density that is proposed to represent PIP2. Although this might be difficult because no density for the headgroup of the lipid was observed.

A final caveat of the study is that effects of PIP2 on the extracellular leaflet might be non-physiological. If this were the case, however, the identification of a non-physiological binding site that favors desensitization might still be beneficial for drug design in the context of this channel.

Reviewer #1 (Public Review):

In this study, Luo, Han, and Yin et al. conduct a fecal microbiota transplant from MSTN KO pigs exhibiting hypertrophy to recipient antibiotic-depleted B6 mice. The microbiota transplants successfully transferred muscle hypertrophy phenotypes to the mice. Aspects of the pig gut microbiome were recapitulated in the recipient mice, including a higher abundance of Romboutsia and valeric acid. The authors then demonstrate that 5 weeks of daily gavage of valerate, but not isobutyrate or water, was sufficient to increase type IIb myofiber growth and GA muscle mass, and protect mice against dexamethasone-induced muscle atrophy. Taken together, these data neatly demonstrate that genetic disruption of the myostatin gene results in a microbiome-dependent increase in valeric acid, which in turn results in significantly altered skeletal muscle growth.

Reviewer #1 (Public Review):

The authors have used eye image data to create an aging clock of the retina in data from eyePACS with validation in the UK Biobank. They show that the clock predicts mortality independently of chronological age and that it is correlated with phenotypic age. Moreover, a GWAS is conducted in the UK Biobank, which identifies novel genetic loci and a top site located in the ALKAL2 region that is functionally validated in a drosophila model. Overall, the study is interesting with sound methodology and is a nice contribution to the field providing a GWAS summary statistic of the eye clock useful for follow-up analyses.

Reviewer #1 (Public Review):

This work leverages single-cell RNA-sequencing to probe changes in various immune functionings within the ovary in aging. The data provided is the most comprehensive of ovarian immune cells at the resolution of single-cell transcriptomics to-date and will be valuable to other researchers. The authors explore four distinct immune functionings:

- The authors identify macrophages and a unique CD3+CD8-CD4- T-cell subpopulation that change in abundance with aging. While these are interesting findings that align with flow cytometry results, the lack of batch correction and application of single-cell differential abundance tools limit the strength of the claims. The authors also do not further probe gene expression changes specific to these populations.

- The authors also analyze changes in global gene expression across cell types using an enrichment analysis; Figure 3B specifically is an excellent visualization summarizing potential global and cell-type specific changes in gene expression programs during aging.

- The authors infer differences in cell-cell communication mediated by various chemokines and cytokines. In this analysis, they claim a decreased inflammatory response due to aging. Here, the global decrease in gene expression in many cell types is not accounted for. Visualizations and quantitative analyses could benefit from existing, specialized in cell-cell communication tools.

- A discussion of changes to the expression of SASP receptors on immune cell types.

While both the data and biology presented are quite interesting, this study is perhaps too wide in breadth such that no individual result is extensively and rigorously explored.

Reviewer #1 (Public Review):

This study presents an implementation of single-particle tomography within the Bayesian framework of the Relion software package. Similar to previously proposed strategies, the approach leverages single-particle analysis tools and tomographic geometric constraints to improve map resolution. Results on the EMPIAR-10164 benchmark dataset appear to match the performance of previous methods, but no maps were made available or deposited, and no direct comparisons with previous results are shown. Consistent with previously published strategies that use 2D projections instead of sub-volumes, the approach performs favorably in terms of resolution when compared to traditional subvolume averaging.

Strengths

- Use of a Bayesian framework for image refinement and reconstruction requires less parameter tweaking.<br /> - By making the new implementation accessible through a GUI already familiar to many SPA users, this tool will make SPT easier to use.<br /> - The implementation of 3D classification could be potentially beneficial to study sample heterogeneity in situ.<br /> - In cases where high resolution can be achieved (better than 3A), the approach has the potential to correct for higher-order optical aberrations.<br /> - Using two cryo-ET datasets, resolution improvements are shown over traditional subvolume averaging (as implemented in the AV3 Matlab suite of programs [Forster et al., 2007] and Dynamo suite [Castano-Diez, 2012]).

Weaknesses

- The approach recapitulates previously proposed strategies for SPT refinement that use raw tomographic projections instead of sub-volumes to improve resolution. Strategies that leverage the increased SNR of average structures to optimize particle pose and deformation, tilt-series alignment, and CTF refinement, were proposed and validated in earlier studies [1,24].<br /> - Compared to end-to-end pipelines for tomography data analysis such as EMAN2 and Dynamo, this approach only implements the subtomogram averaging step, while still relying on external tools for initial tilt-series alignment, CTF estimation, and particle picking.<br /> - In terms of performance, the HIV-1 Gag maps obtained from the benchmark dataset EMPIAR-10164 do not represent an improvement in resolution over previous methods.<br /> - No validation is provided to support the claim that the tool can correct for higher-order optical aberrations of the microscope from cryo-ET data.<br /> - No results are provided to validate the 3D classification routines to study heterogeneity, and no experiments are shown to support the claim that the new approach is more accurate than previous sub-volume classification strategies that compensate for the missing wedge (such as the approach implemented in the earlier version of Relion [4]).

Overall, this implementation of SPT would be a valuable resource for the cryo-ET community.

Reviewer #1 (Public Review):

The authors present a strong set of experiments to uncover what type of role non-mutant stromal cells might be playing in the development of VM and AST, two vascular lesions that share some similarities.

Questions about experimental design.

1) For quantification of gene expression in VM and AST specimens in Figure 2, the methods say qPCR data were normalized to housekeeping genes, but it would be helpful to normalize to endothelial content. It might be that increased TGFa is due to increased endothelium.

2) The mutant allelic frequency for the HUVEC-PIK3CA WT versus HUVEC-PIK3CA H1047R should be provided. This is critically needed for the interpretation of the results.

3) From Figure 5, it appears that the human primary fibroblasts are not required for the mutant ECs to form perfused vessels (panel H). Is it possible that TGFa from the ECs is sufficient to drive vascular malformation?

Reviewer #1 (Public Review):

The current study melds computational and docking methods with functional measurements in a systematic approach: first, they analyze the mechanism of inhibitor binding to EAAT2; second, they mutate ASCT to resemble EAAT and show that the general binding pocket and inhibition mechanism are conserved; third, they perform an in silico screen to identify compounds that bind to the WT ASCT binding pocket; fourth, they perform electrophysiological assays showing that this novel compound allosterically modulates ASCT function. This is a complete and comprehensive study with extensive experimental support for the major conclusions. The authors identify an allosteric ASCT inhibitor, and although only partial inhibition is achieved, this study serves as proof-of-concept that this site can be targeted in diverse SLC-1 transporters as an allosteric inhibitory site.

Reviewer #1 (Public Review):

Using simultaneous EEG-fMRI, authors asked whether neurovascular coupling is already functional in preterm-born neonates, in whom the underlying physiological mechanisms may still be immature at several levels. The question is very interesting and has implications for the study of brain function development as well as early brain injuries. The manuscript reports a correlation between the "mean duration of EEG microstates" and "fMRI BOLD signal-change", through which authors suggest that such a relationship between the EEG activity and BOLD signal highlights the functionality of neurovascular coupling already at the preterm period. The methodology is interesting, but more (not extensive) analysis is required to support the main conclusion and explain the results.

1. The main finding of the study in support of the conclusion comes from relating the inter-individual variability between EEG microstate duration and fMRI BOLD signal change. Given the few subjects (n=13), small even for neuroimaging in infants, studying effects based on inter-individual variability needs to be done with extra care. It is thus important to check whether interindividual variability can be observed for/accounted for by more basic effects in this population :

- The age range is relatively large (age at scan 31 PMA to 36 PMA - but also the age at birth: 29 to 35 weeks) for the number of included infants. Given the intense age-related changes in brain development at this period, it is important to take this factor into account and study it and to have them perhaps explicitly addressed in the manuscript: a ) Does the duration of EEG microstates depend on the age of the infants? b ) Does the time-to-peak in BOLD decrease with increasing age (Arichi et al., 2012)? c) and eventually does the relation between microstate duration and BOLD signal change holds once controlled for their common dependency on the age (i.e. Once partial correlations are used)?

2. The mean/std for the number of epochs per infant can be detailed more. What was the minimum number of epochs? Did such variability in the number of epochs impact microstate properties such as global explained variance/duration? How variable GEV was across infants and would that relate to the variability in duration?

3. Given that sensory-driven changes in microstates follow a sequential pattern (Hu et al., NeuroImage, 2014), could some "microstate syntax" characterize the underlying brain dynamics during stimulation processing in these neonates? Studying the presence of such syntax could be a way to show structured sensorimotor processing, and to further help quantify the inter-individual variability.

4. Is the sleep state monitored (from the EEG signal itself for example)? Given that the sleep state affects EEG activity and in particular EEG microstate properties in newborns (Khazaei et al., Brain topography, 2021), is there a way to rule out that the variability in microstate duration/BOLD signal change is not due to vigilance states?

5. Some of the conclusions/discussion points could be more cautiously stated and developed. On page 7: "However, our results imply that immature neurovascular coupling may not have a significant role in the pathophysiology of cerebral tissue injuries typically seen in preterm born infants (Volpe, 2009); and even that clinical interventions for perinatal brain injury could account for, accommodate, or capitalize on the presence of neurovascular coupling in the preterm human brain to minimize the severity of the injury and its long-term consequences." With an age range covering very preterm infants to late preterm period, generalizing such a conclusion could be potentially misleading for younger infants for example (fNIRS work in younger preterms does not support neurovascular coupling - Nourhashemi et al, Human brain mapping, 2020). As a group - on average - such a pattern may be reported, but the number of infants at each age does not allow us to draw a conclusion about the developmental stage at which such coupling is truly in place. These points could be more directly discussed with regards to the previous literature.

Reviewer #1 (Public Review):

This manuscript reports a series of studies that evaluate the role of long descending propriospinal neurons arising in the cervical spinal cord that project axons to the lumbar spinal cord in locomotor function recovery after spinal cord injury. The experiment uses several different evaluations of gait including BBB, ladder rung walk tests, and kinematics to compare walking before and after synaptic silencing of long descending propriospinal neurons projecting axons to L2. The data reveal that silencing of these neurons mildly improves walking function. The experiments are carefully described and well-controlled. The use of several different methods to evaluate locomotor function is a strength as is a well-thought-out approach to synaptic silencing. The data support the conclusions proposed by the authors. There are caveats to be considered in interpreting the results which are thoughtfully and thoroughly articulated in the discussion.

Public Review:

The study of Choi and collaborators provide novel information about the microstructural morphology and the crystallographic structure of palaeognathid eggshells.<br /> In terms of format and structure, the work is well organized and the extinction of each section is appropriate. All figures, both those from the main text and Supplementary Information, are of good-quality, informative, and useful, facilitating the understanding of the text. The bibliography is very updated, and all essential references are mentioned.

One of the strongest points of the work, in my opinion, is the designee of the study itself, which included specimens from all living palaeognathid birds and several extinct taxa from a large range of lineages.

The methodology used for analysing the crystallographic nature of the studied specimens (EBSD) is appropriate for the goals of the study. The phylogenetic approaches are also right, which are based on the most recent studies about the phylogenetic relationship of ratites.

Despite their complexity, the results are well presented, being relatively easy to understand for a person not versed in the subject. In fact, the ways in which they are described give them the potential to be used as a guideline to anyone interested in eggshell microstructure.

The discussion of the results seems consistent with the data obtained. Despite the phylogenetic relationship between some palaeognathid taxa remains partially instable, authors present different plausible scenario to explain the variability of the eggshell microstructure within a single monophyletic lineage (homology vs homoplasy). In fact, the homoplastic scenario is, perhaps, the most shocking one to me. In part, it is because it intrinsically suggests that all phylogenetic studies based on eggshell morphological features, and conducted during the last 20 years, are potentially artefacts, and they do not represent real phylogenetic relationships. Far from being a criticism, this interpretation has massive implications, especially for those studies where the taxonomic attribution of a fossil egg is based on phylogenetic results (i.e. Montanoolithus, Cairanoolithus).

Although I do not find negative arguments for any special section of the study, I have a question regarding Triprismatoolithu stephensis:

As mentioned in the text, Triprismatoolithu is analysed by the authors, and several pictures are provided in Fig.S12 alongside a brief description in de Supplementary Tex4. But it seems that it is not included in any of the phylogenetic analyses or figures. Why?

If the specimen has no implication for any of the main analyses, there is no need to be considered as "studied material".

Reviewer #1 (Public Review):

Zeng and colleagues investigated the neural underpinnings of visual-vestibular recalibration. Specifically, they measured changes in three monkeys' perception of unisensory heading cues as well as associated changes in neuronal responses to these cues in three different cortical areas following prolonged exposure to systematic visual-vestibular discrepancies. Behavioral responses in a motion direction discrimination task indicate unisensory perceptual shifts in opposite directions that account for the cross-modal discrepancy the monkeys were exposed to. Neuronal firing patterns, related to motion discrimination judgments by means of neurometric functions indicated analogous shifts in neuronal tuning in areas MSTd and PIVC. In contrast, in area VIP tuning for visual heading stimuli shifted in the same direction as tuning for vestibular stimuli and thus in contradiction to the observed perceptual shifts.

The shifts observed in MSTd and PIVC fit nicely with existing theories and results regarding cross-modal recalibration and substitute claims that activity in these areas might underlie perceptual decisions. The shift of visual tuning in VIP is surprising and will certainly spark further investigation.

Overall the results are really interesting, yet, the manuscript in its current form needs revisions along two dimensions, 1) data analysis and 2) writing.

Reviewer #1 (Public Review):

The manuscript by Coates et al. from the Brown lab adds a fascinating and colorful set of tiles to the growing mosaic of small molecule control of the sterol pathway through strategic employment of different parts of the proteostasis pathway. Dr. Brown is an active and creative leader in this field, and this story brings some new and surprising twists to our understanding of the ways that metabolites, and potentially other small molecules, can alter protein processing and life cycle as part of normal cellular function or pathophysiological states. The data are convincing and thorough, and do a great job of revealing many mechanistic aspects of the intriguing observation that hypoxia changes SM processing and activity by altering its degradative fate. The contributing parts of the whole process include altered MARCH 6 E3 ligase activity, new metabolite-ligand regulators (squalene), and ligand-dependent escape from the proteasome to allow the production of a novel form of SM that is freed from the normal regulation of the full-length protein caused by cholesterol, as the authors have previously described. I particularly appreciate three aspects of this study.

First, they test a lot of hypotheses to gain a very full understanding of the gears that are turning to make this hypoxia response machine run. Importantly, these studies also rule out some oxygen sensing mechanisms that work in other contexts, like proline hydroxylation. Second, the authors go to great lengths to integrate the action of the moving parts in a quantitative way, to ascertain if the effects are explained by the coordinated separate changes that are occurring when hypoxia is imposed. And third, the work includes a very well-thought-out set of ideas about why this sort of response is occurring, both in normal cells experiencing either transient or long-term hypoxia, as well as in cancer cells that seem to prefer this form of truncated and alternatively regulated SM.

There is a growing interest in studying and harnessing small molecules to alter and affect protein stability, and these studies add weight to the idea that there are many evolved mechanisms that can teach us lessons both about foundational biology, and new approaches to drug discovery. These beautiful studies will be an important addition to the literature and will be read and referenced by many.

Reviewer #1 (Public Review):

Using two openly available multi-task fMRI datasets, the authors decompose thalamic activity into a smaller set of components. They show that voxels with higher loadings on the main components (high task hub property) also have a high participation coefficient as derived from resting state data. Cortical activity patterns can be predicted to some degree from thalamic activity patterns, and generally better than from a number of other cortical areas. This prediction relies mainly on the voxels with high task hub scores. The results are valuable and methodological generally solid, with some aspects being incomplete.

1. The finding that thalamic activity exhibits a low dimension structure is in my opinion less of a finding, but rather an assumption that motivates the use of dimensionality reduction techniques. When the authors ask (line 101) "whether thalamic task activity exhibits similar low dimensional structure", what is the alternative hypothesis? I think it is a foregone conclusion that with a restricted number of tasks, and the intrinsic smoothness of fMRI activity data, there are always K<<N components that capture 50,75, 90% of the variance. If you had measured the spiking of the entire population of thalamic neurons or increased the threshold to 99%, the structure of activity would be more high dimensional. So I believe you can either frame this as an assumption going in, or you build carefully an alternative hypothesis of what a "high-dimensional" structure would look like. Generating activity data i.i.d would be the simplest case, but given that both signal and measurement noise in fMRI are reasonably smooth, this would be a VERY trivial null hypothesis.

2. The measure of "task hub" properties that is central to the paper would need to be much better explained and justified. You motivate the measure to be designed to find voxels that are "more flexibly recruited by multiple thalamic activity components", but it is not clear to me at this point that the measure defined on line 634 does this. First, sum_n w_i^2 is constrained to be the variance of the voxel across tasks, correct? Would sum_n abs(w) be higher when the weights are distributed across components? Given that each w is weighted by the variance (eigenvalue) of the component across the thalamus, would the score not be maximal if the voxel only loaded on the most important eigenvector, rather than being involved in a number of components? Also, the measure is clearly not rotational invariant - so would this result change after some rotation PCA solution? Some toy examples and further demonstrations that show why this measure makes sense (and what it really captures) would be essential. The same holds for the participation index for the resting state analysis.<br /> 3. For the activity flow analysis, the null models (which need to be explained better) appear weak (i.e. no differences across tasks?), and it is no small wonder that the thalamus does significantly better. The Pearson correlations are not overwhelmingly impressive either. To give the reader a feel for how good/bad the prediction actually is, it would be essential that the authors would report noise ceilings - i.e. based on the reliability of the cortical activity patterns and thalamic activity patterns, what correlation would the best model achieve (see King et al., 2022, BioRxiv, as an example).<br /> 4. Overall it has not been made clear what the RDM analysis adds to the prediction of the actual activity patterns. If you predicted the activity patterns themselves up to the noise ceiling, you would also hit the RDM correctly. The opposite is not the case, you could predict the correct RDM, but not the spatial location of the activity. However, the two prediction performances are never related to each other and it remains unclear what is learned from the latter (less specific) analysis.

Reviewer #1 (Public Review):

The authors serendipitously discovered that silencing Reln+ stellate neurons from medial entorhinal cortex layer II (mEC2) transiently by hyperpolarizing them causes them to degenerate. They replicate this result with two different tools to hyperpolarize these neurons, as well as with a tool to inhibit synaptic vesicle release at mECII axon terminals. They gain mechanistic insight into the degeneration process by performing a careful time course of axon morphological changes and caspase activation: somatic hyperpolarization causes axon retraction bulbs, while inhibition of glutamate release causes axon fragmentation. Crucially, they find that, unlike mEC2 neurons, neighboring Wfs1+ pyramidal cells or parasubicular cells do not degenerate when silenced in similar ways.

The vulnerability of mEC2 to inactivity is particularly compelling because the authors use three different tools to demonstrate it, two that hyperpolarize neurons (ivermectin-mediated activation of the modified glycine receptor alpha subunit, expressed transgenically; and Kir2.1 overexpression using AAV stereotaxic injection), and one that inhibits synaptic vesicle release at mEC2 terminals (Tetanus toxin overexpression using AAV stereotaxic injection). Each of these tools has its flaws but taken together the findings are very convincing. A few pieces of evidence that the various tools are achieving exactly what the authors say they are achieving are missing. But again, the convergence of the data between the three tools compensates for this to some extent.

I found the significance of the findings really fundamental and the writing of the paper absolutely remarkable - beautifully structured, crystal clear in its articulations and its implications. This paradigm has the potential to reveal crucial biology about plasticity in the adult, and about degeneration and vulnerability mechanisms. Vulnerability is such an important topic common to most neurodegenerative diseases, with absolutely no hints, until now, of what could render some cells more prone to degeneration, and immense potential for the discovery of central disease mechanisms. Even if degeneration relies here on the overexpression of an exogenous protein, it does not rely on the overexpression of a pathological protein directly associated with neurodegeneration, or on the invalidation of an essential protein. There is nothing trivial about the degeneration phenotype observed here, which makes the observations absolutely fascinating. What's more the authors show here evidence for the Grail of vulnerability: the side-by-side comparison of two similar/neighboring cell types treated in the same way, only one of which undergoing degeneration (Reln+ EC2 neurons Wfs1+ EC2 and parasubiculum neurons). The vulnerable cell type here also happens to be the very cell type that is most vulnerable to degeneration in Alzheimer's disease.

These findings are of major importance for a few different reasons:<br /> - Neuronal excitability is clearly an early event occurring in the EC of incipient Alzheimer's disease. This study suggests that the silencing of certain cells by Alzheimer's lesions might contribute to their degeneration.<br /> - A competition-based mechanism for the survival or degeneration of axons and neurons from EC2, is known to operate during development until the end of critical periods. This study suggests that EC2 neurons, which might be particular for their need to be plastic into adulthood, might use these mechanisms as well.<br /> - Again, they establish a paradigm for the mechanistic study of comparative vulnerability between cell types that can be investigated further to understand the molecular underpinnings of degeneration.

Reviewer #1 (Public Review):

By studying the effect of Treg depletion in a CD8+ T cell-dependent diabetes model the group around Ondrej Stepanek described that in the absence of Treg cells antigen-specific CD8+ OT-I T cells show an activated phenotype and accelerate the development of diabetes in mice. These cells - termed KILR cells - express CD8+ effector and NK cell gene signatures and are identified as CD49d- KLRK1+ CD127+ CD8+ T cells. The authors suggest that the generation of these cells is dependent on TCR stimulation and IL-2 signals, either provided due to the absence of Treg cells or by injection of IL-2 complexed to specific anti-IL-2 mAbs. In vivo, these cells show improved target cell killing properties, while the authors report improved anti-tumor responses of combination treatments with doxorubicin combined with IL-2/JES6 complexes. Finally, the authors identified a similar human subset in publicly available scRNAseq datasets, supporting the translational potential of their findings.

The conclusions are mostly well supported, except for the following two considerations:

1) From Fig. 4A and B it is not conclusively shown, that Tregs limit IL-2 necessary for the expansion of OT-I cells and subsequent induction of diabetes. An IL-2 depletion experiment (e.g. with combined injection of the S4B6 and JES6-1 antibodies) would further strengthen this claim. Along these lines, the authors claim "IL-2Rα expression on T cells can be induced by antigen stimulation or by IL-2 itself in a positive feedback loop [20]. Accordingly, downregulation of IL-2Rα in OT-I T cells in the presence of Tregs might be a consequence of the limited availability of IL-2.". The cited reference 20 did observe CD25 upregulation by IL-2 on T cells but the observed effect might only be caused by upregulation of CD25 on Treg cells, which increases the MFI for the whole T cell population. Did the authors observe significant upregulation of CD25 on effector CD4+ and CD8+ T cells in their experiments with IL-2/S4B6 or IL-2/JES6 treatment?

2) The anti-tumor efficacy of KILR cells is intriguing but currently, it is unclear if it is indeed mediated by KILR cells. Have KILR cells been identified by flow cytometry in the BCL1 and B16F10 models treated with doxorubicin and IL-2/JES6? Were specific KILR cell depletion studies conducted, e.g. with an anti-KLRK1 depleting antibody? Additional experiments addressing these questions would be desirable to further support the authors' claims.

Reviewer #1 (Public Review):

The authors are trying to determine how time is valued by humans relative to energy expenditure during non-steady-state walking - this paper proposes a new cost function in an optimal control framework to predict features of walking bouts that start and stop at rest. This paper's innovation is the addition of a term proportional to the duration of the walking bout in addition to the conventional energetic term. Simulations are used to predict how this additional term affects optimal trajectories, and human subjects experiments are conducted to compare with simulation predictions.

I think the paper's key strengths are its simulation and experimental studies, which I regard as cleverly-conceived and well-executed. I think the paper's key weakness is the connection between these two studies, which I regard as tenuous for reasons I will now discuss in detail.

The Title asserts that "humans dynamically optimize walking speed to save energy and time". Directly substantiating this claim would require independently manipulating the (purported) energy and time cost of walking for human subjects, but these manipulations are not undertaken in the present study. What the Results actually report are two findings:<br /> 1. (simulation) minimizing a linear combination of energy and time in an optimal control problem involving an inverted-pendulum model of walking bouts that (i) start and stop at rest and (ii) walk at constant speed yields a gently-rounded speed-vs-time profile (Fig 2A);<br /> 2. (experiment) human subject walking bouts that started and stopped at rest had self-similar speed-vs-time profiles at several bout lengths after normalizing by the average duration and peak speed of each subject's bouts (Fig 4B).<br /> If the paper established a strong connection between (1.) and (2.), e.g. if speed-vs-time trajectories from the simulation predicted experimental results significantly better than other plausible models (such as the 'steady min-COT' and 'steady accel' models whose trajectories are shown in Fig 2A), this finding could be regarded as providing indirect evidence in support of the claim in the paper's Title. Personally, I would regard this reasoning as rather weak evidence - it would be more accurate to assert 'brief human walking bouts look like trajectories of an inverted-pendulum model that minimize a linear combination of energy and time' (of course this phrasing is too wordy to serve as a replacement Title -- I am just trying to convey what assertion I think can be directly substantiated by the evidence in the paper). But unfortunately, the connection between (1.) and (2.) is only discussed qualitatively, and the other plausible models introduced in the Results are not revisited in the Discussion. To my naive eye, the representative 'steady min-COT' trace in Fig 2A seems like a real contender with the 'Energy-Time' trace for explaining the experimental results in Fig 4, but this candidate is rejected at the end of the third-to-last paragraph in the 'Model Predictions' subsection of Results based on the vague rationale that is never revisited.

An additional limitation of the approach not discussed in the manuscript is that a fixed step length was prescribed in the simulations. The 'Optimal control formulation' subsection in the Methods summarizes the results of a sensitivity analysis conducted by varying the fixed step length, but all results reported here impose a constant-step-length constraint on the optimal control problem. Although this is a reasonable modeling simplification for steady-state walking, it is less well-motivated for the walking bouts considered here that start and stop at rest. For instance, the representative trial from a human subject in Figure 8 clearly shows initiation and termination steps that differ in length from the intermediate steps (visually discernable via the slope of the dashed line interpolating the black dots). Presumably different trajectories would be produced by the model if the constant-step-length constraint were removed. It is unclear whether this change would significantly alter predictions from either the 'Energy-Time' or 'steady min-COT' model candidates, and I imagine that this change would entail substantial work that may be out of scope for the present paper, but I think it is important to discuss this limitation.

With my concerns about the paper's framing and through-line noted as above, I want to emphasize that I regard the computational and empirical work reported here to be top-notch and potentially influential. In particular, the experimental study's use of inexpensive wearable sensors (as opposed to more conventional camera-based motion capture) is an excellent demonstration of efficient study design that other researchers may find instructive. To maximize potential impact, I encourage the authors to release their data, simulations, and details about their experimental apparatus (the first two I regard as essential for reproducibility - the third a selfless act of service to the scientific community).

I think the most important point to emphasize is that the bulk of prior work on human walking has focused on steady-state movement - not because of the real-world relevance (since one study reports 50% of walking bouts in daily life are < 16 steps as summarized in Fig 1B), but rather because steady walking is a convenient behavior to study in the laboratory. Significantly, this paper advances both our theoretical and empirical understanding of the characteristics of non-steady-state walking.

It is also significant to note the relationship between this study, where time was incorporated as an additive term in the cost of walking, with previous studies that incorporated time in a multiplicative discount in the cost of eye and arm movements. There is an emerging consensus that time plays a key role in the generation of movement across the body - future studies will discern whether and when additive or multiplicative effects dominate.

Reviewer #1 (Public Review):

Tomanek and Guet describe the results of an evolution experiment where they allowed the bacterium E. coli to adapt to various concentrations of galactose as an additional carbon source. These conditions impose different degrees of demand for the galK enzyme, whose expression level depends on the promoter sequence and on the number of copies of the galK locus. Given that the initial promoter is random and weak, both amplifications of the locus and mutations in the promoter are expected to be adaptive. The experimental strains of E. coli were equipped with a fluorescent reporter system designed to discriminate between these two types of mutations. Furthermore, two strains, IS+ and IS-, were engineered with high and low rates of duplication around the galK locus, respectively. The main result is that at higher concentrations of galactose, where the demand for galK is high, E. coli adapts by acquiring combinations of both types of adaptive mutations, amplifications, and promoter mutations. In contrast, at low concentrations of galactose, where the demand for galK is low but not zero, E. coli appear to adapt by acquiring either an amplification or a promoter mutation but not their combinations. The observation of apparent interference between the acquisition of these two types of mutations is interesting and novel. The authors provide an intuitive explanation for it: when one mutation is sufficient to achieve the optimal expression of the gene, the mutation that is acquired first makes the other mutation obsolete, i.e., there is negative epistasis (possibly even sign epistasis) between these mutations, in the sense that the second mutation is much less adaptive (or possibly even deleterious) in the presence of the first one, in the low-demand environment. The authors discuss the possible implications of this finding for our understanding of molecular evolution and propose a new Amplification Hindrance hypothesis. This hypothesis states that, since amplifications occur at much higher rates than individual point mutations, they can slow down or even prevent sequence divergence. The amplification hindrance hypothesis stands in contrast to the Innovation-Amplification-Divergence hypothesis which is currently the default paradigm and states that amplifications generally accelerate sequence divergence.

STRENGTHS:

The authors designed a powerful reporter system that allows them to monitor the evolutionary dynamics of amplifications and promoter mutations. They ask an important question: how do early evolutionary dynamics of adaptation to environments with different demands for gene expression look like? The phenotypic data they present looks very interesting and shows the existence of interference between amplifications and point mutations in low-demand but not in high-demand conditions. The Amplification Hindrance hypothesis is a novel and useful intellectual contribution to the field.

WEAKNESSES:

In my opinion, the main weakness of the paper is that, while the interference between amplifications and point mutations in the low-demand condition clearly happens (most convincingly shown in Figure 5), its causes remain unclear. In particular, the authors claim that this interference is caused by negative epistasis, but the possibility of clonal interference without epistasis has not been decisively ruled out. The authors mention clonal interference tangentially in the Discussion, but they do not seriously address this alternative explanation. Yet, understanding the cause of this phenomenon is important because clonal interference and negative epistasis have different implications for long-term evolution.

The authors' main hypothesis is that, in the low-demand conditions, expression-increasing point mutations in the promoter provide much lower fitness benefits (or even incur fitness costs) in strains with galK amplifications compared to the ancestral strain without amplifications. The most direct way to test this hypothesis would be to measure the fitness effects of a point mutation in genetic backgrounds with and without amplifications in conditions with low and high demand for galK. This decisive experiment has unfortunately not been done. Instead, the authors construct an indirect argument, whose essence is as follows.

They show that, over the course of the experiment in the low-demand environment, the IS+ populations have acquired fewer point mutations than IS- populations (Figure 5). In addition, the phenotypic data in Figures 2 and 4 demonstrate that IS+ mutations in the low-demand environment contain three phenotypic classes of cells: ancestral, YFP+ and YFP+CFP+. The YFP+ clones are shown to have only one or two promoter mutations. The YFP+CFP+ cells must have duplications, and it is likely (although not quite certain, see below) that they do not have any promoter mutations. These data demonstrate quite convincingly that, whenever adaptation by duplications is possible, the rate at which point mutations segregate and accumulate declines. These data are consistent with the authors' hypothesis based on negative epistasis. However, they also seem to be consistent with the idea that amplifications and point mutations exhibit clonal interference without negative epistasis.

It may be possible to construct an argument against this alternative hypothesis based on the comparison between different environments, but such an argument would have to take into account the fact that clonal interference depends not only on the rates of mutations (which are presumably the same in all environments) but also on their fitness effects which vary across the environment. Another possibility to argue against clonal interference might be by carrying out simulations, although this approach also seems challenging without knowing some key population genetic parameters. The most direct way to resolve this ambiguity would be to demonstrate negative epistasis as discussed above.

Another, less critical but still important, issue mentioned above concerns the authors' claim that the YFP+CFP+ cells have only duplications but no promoter mutations (e.g., LL. 276-277). This is certainly consistent with intuition since these cells have an increased level of both YFP and CFP relative to the ancestor. However, as far as I can tell, there is no evidence to support this claim directly. My understanding is that the authors base this claim on the fact that YFP+CFP+ cells form a cluster of points on the YFP vs CFP plots that is distinct from the cluster of "mixed" cells, which are shown to have both an amplification and a promoter point mutation (Figure 3). But it is still logically possible that the YFP+CFP+ cells have an amplification and a promoter mutation other than the one found in the "mixed" cells (e.g., weaker). The most direct way to show that YFP+CFP+ cells have no promoter mutations would be to sequence a few of them. Another possibility would be to calibrate the YFP/CFP fluorescence measurements against galK copy number.

Reviewer #1 (Public Review):

Dectin-1 is a known C-type lectin receptor that has a role in recognizing pathogen glycans, particularly beta-glucan. Haji et al. present evidence that Dectin-1 also has an endogenous ligand, another C-type lectin that is enriched on platelets called CLEC-2. Using a Dectin-1 reporter line, they identify human platelets as a source for the ligand, raising a panel of mAbs to human platelets and identify one (6D11) that prevents signalling and use this to identify CLEC2 as the Dectin-1 receptor by immunopurification. The authors go on to further characterise this interaction by showing that it occurs between the human but not mouse orthologues, showing that it is the stalk region of human Dectin-1 that contains the ligand which consists of sialylated core 1 O-linked glycans displayed on Thr 105 and adjacent amino acids. Finally, the authors over-express human Dectin-1 in mice within the context of a null background for another CLEC-2 ligand (podoplanin) and show that this rescues embryonic lethality. They show that Dectin-1-dependent CLEC-2 signalling is not sufficient to induce platelet aggregation but can rescue perinatal lethality in podoplanin-deficient mice.

There is a large amount of work in this manuscript and the experiments are well controlled and the results unequivocal and usually verified by several independent techniques. The paper is very well-written making the volume of data accessible. There is novelty here in that this unusually finds that 2 C-type lectin receptors directly interact and the biochemical characterisation of the glycan binding determinants is very well performed. While the authors show that the interaction occurs between the human but not mice orthologues (because mouse Dectin-1 lacks the critical EDxxT motif that is necessary for displaying the o-linked glycan in the correct context) they were able to investigate the functional role of the interaction by generating mice that overexpressed human Dectin-1 in a podoplanin-deficient background. Normally, podoplanin-deficient mice die at birth due to defects in lymphatic vessel development, but this lethality is rescued by overexpression of human Dectin-1. This xeno-overexpression data can be difficult to interpret but suggests that while human Dectin-1 signalling to mouse platelet CLEC-2 is insufficient to drive platelet activation and thrombus formation, it is sufficient to rescue some aspects of platelet function involved in lymph vessel development. They conclude that human Dectin-1 (which is broadly expressed on different tissues) provides a basal level of tonic signalling through platelet-expressed CLEC-2 to establish platelet signalling thresholds.

I thought the manuscript was comprehensive in its approach, well written, and that the experimental data supported the conclusions.

Reviewer #1 (Public Review):

The authors tackle an interesting problem: how do ant colonies regulate foraging in response to their collective hunger? In previous work, the authors related the colony's response to individual ants sensing their own food levels and its temporal dynamics. Looking more carefully at the spatial dynamics of ants, the authors now find that foragers tend to move toward the depth of the nest when their food load is high and toward the nest exit when it is low. This is an elegant and computationally inexpensive set of rules that explains the spatiotemporal dynamics of the system.

Overall, the paper is written clearly, the methods are sound, and I agree with the interpretation of the results.

I do have a few comments and suggestions:

1) How exactly are the inward outward directions defined? Is it simply, away, or towards the entrance? It is not clear from the text, and since this system is not symmetric (cubic with entrance at one of the corners) the authors should clarify.

2) To analyze the biased random walk analysis of the ants, the authors "coarse-grained" the steps as being "inwards" "outwards" and "stay". It's not clear how this level of granulation is justified. Since the authors have access to the actual trajectories and all trophallaxis events, why not just calculate the actual turning angles between consecutive steps the ants take? This would give an actual assessment of both the bias and the noise imposed on the random walks, which the authors could then use directly in their models.

3) It would be important to better connect the author's previous mechanism (relating the colony's response to individual ants sensing their own food levels and its temporal dynamics) to the new mechanism (spatial-temporal dynamics). Are they mutually exclusive? It would be useful to elaborate on this in the Discussion.

4) It would be useful to add a few supplementary movies from the experiments, showing ants moving toward the entrance with low food loads, and moving away from the entrance with high food loads.

Reviewer #1 (Public Review):

This paper addresses an important issue of how humans select foot placement when running on uneven terrain. The authors examine empirical and model-based data to determine whether runners specifically opt for more level locations on such surfaces. The manuscript suggests potential strategies of how humans mitigate fore-aft impulses during foot contact so as to maintain stability in response to changes in terrain.

Overall, the manuscript provides additional insight into lower-limb mechanics of runners on natural surfaces. The model-based analysis of lower limb compliance is especially useful in the context of stability and understanding human motor control. In addition, participant foot placement analysis using empirical and statistical models is very compelling and provides insight into how much planning occurs during running on uneven terrain. However, there are a few concerns of note:

- A central motivation of the study appears to be that past research did not incorporate height and slope variations when studying gait on uneven terrain. Although some of the past work cited by the authors does focus on step-like terrain, the (Voloshina, Ferris 2015) study had both height and slope variations. In that particular study, the terrain consisted of blocks that were smaller than the dimensions of the average human foot. This means that, during each foot-flat phase, the foot had to span at least two blocks of different heights, placing it at a slope in the fore-aft direction. Similarly, the columns of that terrain layout provided slope variations in the medio-lateral direction. Referring to this surface as "step-like" is inaccurate and potentially misleading. Considering that the terrains in the present study and the study from 2015 likely cause very similar types of perturbations to the runner, the motivation behind the current study is not strongly validated. The authors should consider re-evaluating their results in the context of past studies.

- A primary outcome of the study is the analysis of empirical and model-based fore-aft impulses experienced by runners during foot contact. The authors suggest that this measure is related to stability but do not provide extensive explanations. It would be helpful to include additional background information on how impulse analyses have been used previously and why they are particularly fitting in this context.

- The study evaluates participants running back and forth on a 24m track for up to 10 minutes. This means that participants had to perform many turns during each trial. The authors present metabolic energy expenditure data but do not address how these data may be skewed due to the large instances of directional changes. Measuring metabolic data during such tasks is generally noise-prone, potentially leading to an inaccurate representation of energy expenditure. Considering that this is a comparative study and participants had to perform such turns for each terrain trial, this issue could be minor. However, the authors are encouraged to provide more detail on experimental protocol. Addressing whether or not participants stepped off the terrain to switch directions and providing insight into how this experimental approach could potentially affect outcomes could be particularly helpful.

Reviewer #1 (Public Review):

This is a fascinating paper that takes up an important question with a creative and new approach. We have a few suggestions that we hope are constructive for the authors.

1. One nagging concern is that the category structure in the CNN reflects the category structure baked into color space. Several groups (e.g. Regier, Zaslavsky, et al) have argued that color category structure emerges and evolves from the structure of the color space itself. Other groups have argued that the color category structure recovered with, say, the Munsell space may partially be attributed to variation in saturation across the space (Witzel). How can one show that these properties of the space are not the root cause of the structure recovered by the CNN, independent of the role of the CNN in object recognition?

2. In Figure 1, it could be useful to illustrate the central observation by showing a single example, as in Figure 1 B, C, where the trained color is not in the center of the color category. In other words, if the category structure is immune to the training set, then it should be possible to set up a very unlikely set of training stimuli (ones that are as far away from the center of the color category while still being categorized most of the time as the color category). This is related to what is in E, but is distinctive for two reasons: first, it is a post hoc test of the hypothesis recovered in the data-driven way by E; and second, it would provide an illustration of the key observation, that the category boundaries do not correspond to the median distance between training colors. Figure 5 begins to show something of this sort of a test, but it is bound up with the other control related to shape. Similarly, if the claim is that there are six (or seven?) color categories, regardless of the number of colors used to train the data, it would be helpful to show the result of one iteration of the training that uses say 4 colors for training and another iteration of the training that uses say 9 colors for training. The text asserts that Figure 2 reflects training on a range of color categories (from 4 to 9) but doesn't break them out. This is an issue because the average across these iterations could simply be heavily biased by training on one specific number of categories (e.g. the number used in Figure 1). These considerations also prompt the query: how did you pick 4 and 9 as the limits for the tests? Why not 2 and 20? (the largest range of basic color categories that could plausibly be recovered in the set of all languages)?

3. Regarding the transition points in Figure 2A, indicated by red dots: how strong (transition count) and reliable (consistent across iterations) are these points? The one between red and orange seems especially willfully placed.

4. Figure 2E and Figure 5B are useful tests of the extent to which the categorical structure recovered by the CNNs shifts with the colors used to train the classifier, and it certainly looks like there is some invariance in category boundaries with respect to the specific colors uses to train the classifier, an important and interesting result. But these analyses do not actually address the claim implied by the analyses: that the performance of the CNN matches human performance. The color categories recovered with the CNN are not perfectly invariant, as the authors point out. The analyses presented in the paper (e.g. Figure 2E) tests whether there is as much shift in the boundaries as there is stasis, but that's not quite the test if the goal is to link the categorical behavior of the CNN with human behavior. To evaluate the results, it would be helpful to know what would be expected based on human performance.

5. The paper takes up a test of color categorization invariant to luminance. There are arguments in the literature that hue and luminance cannot be decoupled-that luminance is essential to how color is encoded and to color categorization. Some discussion of this might help the reader who has followed this literature. Related, the argument that "neighboring colors in HSV will be neighboring colors in the RGB space" is not persuasive. Surely this is true of any color space?

6. The paper would benefit from an analysis and discussion of the images used to originally train the CNN. Presumably, there are a large number of images that depict man-made artificially coloured objects. To what extent do the present results reflect statistical patterns in the way the images were created, and/or the colors of the things depicted? How do results on color categorization that derive from images (e.g. trained with neural networks, as in Rosenthal et al and presently) differ (or not) from results that derive from natural scenes (as in Yendrikhovskij?).

7. It could be quite instructive to analyze what's going on in the errors in the output of the classifiers, as e.g. in Figure 1E. There are some interesting effects at the crossover points, where the two green categories seem to split and swap, the cyan band (hue % 20) emerges between orange and green, and the pink/purple boundary seems to have a large number of green/blue results. What is happening here?

8. The second experiment using an evolutionary algorithm to test the location of the color boundaries is potentially valuable, but it is weakened because it pre-determines the number of categories. It would be more powerful if the experiment could recover both the number and location of the categories based on the "categorization principle" (colors within a category are harder to tell apart than colors across a color category boundary). This should be possible by a sensible sampling of the parameter space, even in a very large parameter space.

9. Finally, the paper sets itself up as taking "a different approach by evaluating whether color categorization could be a side effect of learning object recognition", as distinct from the approach of studying "communicative concepts". But these approaches are intimately related. The central observation in Gibson et al. is not the discovery of warm-vs-cool categories (these as the most basic color categories have been known for centuries), but rather the relationship of these categories to the color statistics of objects-those parts of the scene that we care about enough to label. This idea, that color categories reflect the uses to which we put our color-vision system, is extended in Rosenthal et al., where the structure of color space itself is understood in terms of categorizing objects versus backgrounds (u') and the most basic object categorization distinction, animate versus inanimate (v'). The introduction argues, rightly in our view, that "A link between color categories and objects would be able to bridge the discrepancy between models that rely on communicative concepts to incorporate the varying usefulness of color, on the one hand, and the experimental findings laid out in this paragraph on the other". This is precisely the link forged by the observation that the warm-cool category distinction in color naming correlates with object-color statistics (Gibson, 2017; see also Rosenthal et al., 2018). The argument in Gibson and Rosenthal is that color categorization structure emerges because of the color statistics of the world, specifically the color statistics of the parts of the world that we label as objects, which is the same approach adopted by the present work. The use of CNNs is a clever and powerful test of the success of this approach.

Reviewer #1 (Public Review):

This paper uses 2 cohorts from the UK and links a) risk factors associated with low antibody levels after vaccination and b) risk factors for infection. The paper makes the important point that following the third vaccination, risk factors associated with low antibody response after the first vaccination, are less likely to lead to sub-protective levels. This highlights the importance of obtaining a booster shot. Though it is not a primary finding of the paper, the observed discordance between self-reported infection and anti-nucleocapsid positivity is an important finding. While these findings are potentially useful, the presentation of the data is somewhat unfocused, and the message is presented in a diffuse fashion. Moreover, certain key components of the analysis such as the assay threshold and timing of samples after the 2nd vaccine are a bit confusing and require clarification. The use of univariate analyses can be misleading. Finally, the relevance of the findings relative to our current stage of the pandemic with multiple new VOCs requires a clearer explication.

Reviewer #1 (Public Review):

The core conclusions in this manuscript are well supported. First, the genomic data clear support a recent origin of the Baltic and Arctic ecotypes from an Atlantic-like ancestor, without major bottlenecks and with gene exchange at least on the one sampled sympatric location but probably also more widely. The genome scans strongly suggest the involvement of multiple loci in divergence. This is supported by the quantitative trait locus analysis but this support is relatively weak because the number of markers used was small and markers were rather unevenly distributed on the genetic map, leaving little power to detect QTL in some areas. An analysis of the power of the QTL analysis is lacking and there might also be an issue about whether the measured trait truly captures rhythmicity.

The presence of segregating variants across all sampled populations for the loci that appear to underlie local adaptation, plus the absence of strong sweep signatures, is good evidence that adaptation to the Baltic environment was based on standing genetic variation. This is supported by the known presence of local adaptation to tidal regimes within the Atlantic ecotype, providing a mechanism for the maintenance of standing variation. Clustering of putatively adaptive loci in one region of chromosome 1 seems clear, although it is not formally compared to a random distribution.

The authors make an interesting point that it is only when many genes are involved that Gene Ontology enrichment is expected to be informative. Here, they also had clear a priori expectations which are neatly fulfilled, further suggesting that differentiated loci are good candidates for a role in adaptation.

The broader context provided for the analysis of the fascinating marine midge system is brief and could usefully be expanded. It should make clear that, despite some focus on major gene effects that can be assigned to individual loci, there is widespread evidence for polygenic adaptation. Even where structural variants are known to have major effects, many other regions of the genome typically contribute to local adaptation. It would also be helpful to refer to theoretical expectations regarding distributions of effect sizes and clustering of locally-adaptive loci, with or without gene flow.

Reviewer #1 (Public Review):

The authors define regulatory networks across 77 tissue contexts using software they have previously published (PECA2, Duren et al. 2020). Each regulatory network is a set of nodes (transcription factors (TF), target genes (TG), and regulatory elements (RE)) and edges (regulatory scores connecting the nodes). For each context, the authors define context-specific REs, as those that do not overlap REs from any of the other 76 contexts, and context-specific regulatory networks as the collection of TFs, TGs, and REs connected to at least one context-specific RE. This approach essentially creates annotations that are aggregated across genes, elements, and specific contexts. For each tissue, the authors use linkage disequilibrium score regression (LDSC) to calculate enrichment for complex trait heritability within the set of all REs from the corresponding context-specific regulatory network. Heritability enrichments in context-specific regulatory network REs are compared with heritability enrichments in regions defined using other approaches.

Reviewer #1 (Public Review):

Yuan et al. propose that the magnesium transporter unextended (uex) controls Drosophila sleep via Mg2+ efflux, Ca2+-dependent CREB signaling, and a CNK-ERK pathway. UEX protein levels display daily oscillations in fly heads, whereas UEX-depleted flies show long sleep with low levels of Ca2+ and synaptic plasticity. Transgenic expression of wild-type UEX or the mammalian uex homolog CNNM1 possibly rescues the uex mutant sleep, supporting their evolutionary conservation. UEX forms a protein complex with the ERK signaling suppressor CNK, and UEX depletion appears to de-repress ERK activation. As expected, pan-neuronal CNK depletion phenocopies uex mutant sleep. Taken together, the authors suggest a novel mechanism whereby magnesium shapes animal sleep through the specific MAPK pathway. Overall, the authors reported quite a few exciting phenotypes in UEX-depleted flies using a range of analyses (e.g., behaviors, gene expression, Mg2+/neural imaging, and biochemistry). However, evidence for their causal link to sleep regulation is missing, and some key conclusions remain justified by more rigorous analyses. It is noteworthy that Wu et al. have previously demonstrated the Mg2+ efflux transporter activity of UEX and mapped its memory-enhancing function to a specific group of adult neurons in the fly brain (https://pubmed.ncbi.nlm.nih.gov/33242000/). Given the intimate interactions between sleep and memory, these findings may have broad implications in sleep-relevant physiology and disorders. However, a number of important control studies and statistical assessments are not included, but are necessary to conclusively interpret the data.

Reviewer #1 (Public Review):

Modified rabies viral vectors allow high throughput mapping of neuronal circuits with cell type specificity. However, lack of standardization in this field limits extrapolation of useful information, beyond the identity of anatomically connected regions, such as differential input densities and connectivity motifs. in this manuscript, Tran-Van-Minh et al., attempt to develop a statistical approach which will allow consolidation of new, as well as previously-acquired datasets, to yield biologically significant insights into the logic underlying rabies vectors' expansion from single starter cells.

This question the authors address is of high importance and the presentation of this manuscript is timely, as rabies tracing experiments increase at an exponential rate. The authors provide a largely complete description of the main pitfalls and caveats in current analysis approaches and common misconceptions in the interpretation of results. In addition, they correctly diagnose the potential of this methodology to extract pertinent information from such experiments and provide the reader with useful tips on how to better design, analyze and interpret them.

While such a paper is undoubtedly called for, and has the potential to substantially improve circuit mapping experiments, with little cost to the experimenter, there are a few critical flaws in their premises, which will limit a widespread adoption of their analysis approach. While the authors discuss caveats in design and interpretation of results, they do not implement these suggestions into their own experiments, such as the need for accurate estimation of starter cells and automated cell counting which is not biased by strong signal coming from dendritic arbors/axons in densely-labeled regions. While the authors claim that the reduction in residuals and relative conformity across experiments following log-transformation of n(i) and n(s) shows that their approach is robust, the large degree of variability across experiments, the datasets of some are biologically implausible, showing a decrease in n(i) as a function of n(s), suggests that this transformation disposes of useful information, which might help detect anomalies in data acquisition. The fact that the most rigorously-acquired dataset which was presented by the Allen Institute (BRAIN Initiative Cell Census Network, Nature, 2021) is the only one which, does not comply with the transformation, attests further to this caveat. This is probably the reason why the estimated number of individual neurons retrogradely labeled by a single starter cell (~1400) is more than an order of magnitude higher than any previously-published estimation, including those which include tracing from single-cell, and is highly unrealistic.

In addition, the model selected by the authors to fit the various datasets does not take into proper account saturation of n(i), as all proposed functions are growing functions. Here, the most suitable model which should be used to describe the nature of RV expansion is a cumulative distribution function, while the rest are either private cases (e.g. linear fit given low n(s) or zero divergence) or are biologically unrealistic (e.g. exponential fit). Again here, the only dataset which appears to fit this function the best comes from the BRAIN Initiative Cell Census Network (Nature, 2021).

In conclusion, while such work is called for and has the potential of becoming a staple in the connectivity toolbox, many of the premises presented here will need to be significantly adjusted, before the approach could be put into widespread use.

Reviewer #1 (Public Review):

The manuscript by Chen and colleagues contains a number of important findings. First, they showed with careful imaging, cell cycle characterization, and models that replisomes of fast-growing E. coli form factories and super factories (factories involving cousin replisomes), they estimated that factories split after the replication of 1/3 of the chromosome. Second, they used a replication fork block to monitor the interdependence of replisomes of factories. Third, they show that orphaned replisomes replicate slowly and frequently require a repair process to complete the replication cycle. This is the first characterization of an advantage of replicating bacterial genomes in a factory rather than as independent replisomes. This is an important discovery that simultaneously confirms the existence of DNA factories, a much-debated topic, and reveals one of their functions. The existence and characterization of factories were mostly addressed through imaging methods; here the combination of imaging with molecular genomics assays is a real asset for the manuscript.

Reviewer #1 (Public Review):

In this study, the authors investigate multi-attribute economic decisions in humans. More specifically, they investigate the impact of adding a decoy option to a binary choice, the effect of which has been debated in the previous literature with recent studies finding effects in different directions. By re-analyzing large datasets from some of these studies and by applying computational modeling to the data, the authors propose that a subjective encoding rule comparing attributes separately, and not computed as one multiplicative value, explains participants' behavior on different levels of granularity. Context effects were well captured by a selective integration model. The work has many positive features, including praising open science, the analysis of the potential confounds of previously used designs, and both quantitative and qualitative comparisons of several well-justified models.

Reviewer #1 (Public Review):

This manuscript provides the first cellular analysis of how neuronal activity in axons (in this case the optic nerve) regulates the diameter of nearby blood vessels and hence the energy supply to neuronal axons and their associated cells. This is an important subject because, in a variety of neurological disorders, there is damage to the white matter that may result from a lack of sufficient energy supply, and this paper will stimulate work on this important subject.

Axonal spiking is suggested to release glutamate which activates NMDA receptors on myelin-making oligodendrocytes wrapped around the axons: the oligodendrocytes - either directly or indirectly via astrocytes - then generate prostaglandin E2 which relaxes pericytes on capillaries, thus decreasing the resistance of the vascular bed and (presumably) increasing blood flow in the nerve.

Strengths of the paper

The paper identifies some important characteristics of axon-vascular coupling, notably its slow temporal development and long-lasting nature, the involvement of PgE2 in an oxygen-dependent manner, and a role for NMDARs. Rigorous criteria (constriction and dilation of capillaries by pharmacological agents) are used to select functioning pericytes for analysis.

Weaknesses of the paper

The study focuses exclusively on pericytes. It would have been interesting to assess whether arteriolar SMCs also contribute to regulating blood flow.

The slow (~10 minutes) time scale of the responses seen is remarkable when compared to grey matter neurovascular coupling which occurs in seconds. The authors suggest that this reflects movement of messengers along the nerve, but the action potential moves so rapidly down the nerve that it will activate the release of vasodilators essentially at the same time everywhere along the nerve, and there will be no concentration (or voltage) gradient to drive such a movement of messengers (unless there is a spatially localized unique site in the vascular bed where propagating responses are generated). Thus it remains unclear why the responses are so slow.

The activity-evoked dilation is thought to reflect PgE2 release at what is probably a physiological O2 concentration, but not at the hyperoxic 95% O2 used in most of the experiments. The involvement of NMDA receptors is apparently only shown (using OL-specific receptor subunit deletion) at the unphysiological high O2 level. This raises two questions: are NMDA receptors also involved in the response at low O2 (as schematized in Fig 6), and what is the messenger downstream of NMDARs at high O2 (NO would be an obvious candidate, previously shown to contribute to NVC in the grey matter).

Reviewer #1 (Public Review):

The goal of this work is to study whether sleep and wake regulate cerebellar structural plasticity. Scanning electron microscopy and 3D reconstruction were used to characterize structural changes in spines and synapses in the mouse cerebellar cortex. The net number and size of synapses did not change between sleep and wake. However, the number of small portion of spines ("naked" spines without synapses) increased during sleep, and the number of branched spines (contains more than one synapse) increased during wake.<br /> The methodological approach (scanning electron microscopy) is laborious but adequate. However, it cannot follow the same synapses longitudinally, and live imaging, even only in superficial regions, is an ideal approach to achieve the goals.

The results did not find changes in the total number of spines or synapses during sleep and wake. Structural changes were found in small number of spines lacking a synapse. Whether these changes are functionally important is unclear. The results support circuit-specific, rather than global, effect of sleep on synaptic plasticity.

The role of sleep range from cellular maintenance to memory consolidation. Multiple evidence suggests that sleep regulates synaptic plasticity. Although this work did not attempt to understand the functional role of sleep in regulating learning and memory, it discovered that sleep promotes synaptic pruning in specific circuits of the cerebellum. Future similar anatomical studies in additional brain regions, including excitatory and inhibitory circuits, combine with physiological and behavioral assays are expected to provide insights to the role of sleep in regulating synaptic plasticity, learning and memory.

Reviewer #1 (Public Review):

Zebrafish have become a well-established model organism for studies of damage and repair in hair cell sensory organs due to organ accessibility and robust mechanisms for repair and regeneration. While lateral line organs in zebrafish are widely used to study hair cell physiology, zebrafish also contain hair cell organs in their inner ears that have the potential to be more comparable to mammalian counterparts. Using single-cell RNA seq in embryonic through adult stages, this study found that hair cells and supporting cells of the zebrafish inner ear are distinct from those of the lateral line. Additionally, they identified distinct cellular subtypes within the maculae and cristae of the ear.

This work provides some novel findings, including identifying distinct markers in the macula and crista hair cells and supporting cells as well as detecting a domain in the zebrafish utricle that coincides with features of the striolar region in mouse utricle. However, while much of the focus of in situ expression analysis was the spatial separation of calcium-binding proteins capb1b and capb2b in the maculae, it was not clear how Capb1 & 2 expression corresponds to striolar and extrastriolar regions in the mammalian utricle, where Capb2 is expressed throughout. In addition, the authors assumed the zebrafish utricle and saccule perform similar functions (i.e. hearing) and contain hair cells with similar frequency tuning. It would improve this study to consider the unique functions and frequency sensitivities of the utricle and the saccule, given that different expression of voltage-gated channels may give insight into the specific physiology of these two sensory organs.

Reviewer #1 (Public Review):

The authors Rem et al., examine the mechanism of action of APP, a protein implicated in Alzheimer's disease pathology, on GABAB receptor function. It has been reported earlier that soluble APP (sAPP) binds to the Sushi domain 1 of the GABAB1a subunit. In the current manuscript, authors examine this issue in detail and report that sAPP or APP17 interacts with GABABR with nano Molar affinity. However, binding of APP to GABAB receptor does not influence any of the canonical effects such as receptor function, K+ channel currents, spontaneous release of glutamate, or EPSC in vivo. The experimental evidence provided to support the conclusions is thorough and statistically sound. The range of techniques used to address each of the aims has been carefully curated to draw meaningful conclusions.

The authors use HEK293T heterologous cell line to confirm the affinity of APP17 for the receptor, ligand displacement, and receptor activation. They also use this method to study PKA activation downstream of the GPCR. They use slice electrophysiology to measure changes in glutamatergic transmission EPSC and then in vivo 2-photon microscopy to measure functional changes in vivo.

The work is significant for the field of Alzheimer's and also GABAB receptor biology, as it has been assumed for sAPP acts via GABAB receptors to influence neurotransmission in the brain. The results presented here open up the question yet again, what is the physiological function of sAPP in the brain?

The manuscript is clearly written and easy to follow. The main criticism would be that the manuscript fails to identify the mechanism downstream of APP17 interaction with GB1a SD1.

Reviewer #1 (Public Review):

This paper by Hubmann et al investigates the role of a large ECM protein in lymphangiogenesis using the zebrafish model to test genetic interactions and with state-of-the-art reporter readouts to evaluate expression patterns and vascular behaviors. The experiments examine the effects of svep1 deletion on facial lymphatics and produce the surprising result that svep1 and vegfc appear to be required in non-overlapping areas of this lymphatic vascular bed. They then use these tools and mutants for tie1 or tie2 to examine genetic interactions in this area and in parachordal lymphangioblast migration in the trunk, showing that tie1 but not tie2 shows epistasis with svep1. The data are overall rigorous and quantified to show the statistical significance of the stated results. The work provides important insights into a complex signaling pathway that is widely utilized in vascular development. The evidence is convincing in supporting the findings and is predicted to be of interest to vascular biologists and others interested in Tie/Tek signaling such as cancer biologists.

Reviewer #1 (Public Review):

Gormley et al. conduct a comprehensive Mendelian randomisation (MR) analysis to study the causal effect of obesity and metabolic disorder on oral and oropharyngeal cancer. This work follows on from observational studies that found evidence of associations between these variables and continued on from a previous MR study that focused on the effect of circulating lipid traits on these cancer outcomes. In this study, the authors found limited evidence for an effect of obesity-related traits on either cancer type, contradicting the previous observational studies and thus implying that the latter may have been affected by confounding. Sensitivity analyses that adjusted for potential pleiotropy and outlying instrumental variants agreed with the main study findings. Risk factors including smoking, risk-taking (a proxy for sexual behavior), and alcohol consumption were also stratified for, and only evidence was found for smoking as a mediator of the observed effect.

Strengths:<br /> Strengths of this study include the thorough MR analysis conducted, with all required sensitivity analyses and checks conducted and summarised appropriately, including the use of recent MR methods such as SIMEX where required. For the disease outcomes, oral and oropharyngeal cancer, the authors used summary statistics from the single largest trans-ethnic published meta-analysis GWAS available. Known risk factors were stratified for in sensitivity analyses to follow up any results of interest. Where measures of risk factors were not directly available (e.g. sexual behaviour), genetically-correlated proxy measures were used.

Weaknesses:<br /> While this study has several strengths, there are also a number of limitations present, many of which the authors outline in the paper. In particular, the sample size of the outcome GWAS used was rather small which may hinder the power. Also, some of the cancer subtypes of interest that had previously been investigated in observational studies were not available as GWAS, and so could not be included in this study. The authors only used European or trans-ancestry GWAS, as these are the only ancestries presently available, but a future investigation into other ethnicities is warranted. They also used risk-taking (self-defined via questionnaire) as a proxy for sexual behaviour, which despite having a strong genetic correlation, may not be the best substitute.

Reviewer #1 (Public Review):

This is an elegant article where the authors define valuable criteria to identify and classify high-confidence miRNA in fungi. The data supporting the conclusions are solid, and the results are essential to unify the annotation criteria in these organisms. Interestingly, the paper shows that miRNAs in fungi look and position within the genome more similarly to their animal counterpart but may act like plant miRNAs.

The conclusions of this paper are well supported by data, but some aspects need to be clarified and extended.

Reviewer #1 (Public Review):

In this manuscript, the author characterizes the lattice of kinesin-decorated microtubule reconstituted from porcine tubulins in vitro and Xenopus egg extract using cryo-electron tomography and subtomogram averaging. Using the SSTA, they looked at the transition in the lattice of individual microtubules. The authors found that the lattice is not always uniform but contains transitions of different types of lattices. The finding is quite interesting and probably will lead to more investigation of the microtubule lattice inside the cells later on for this kind of lattice transition.

The manuscript is easy to read and well-organized. The supporting data is very well prepared.

Overall, it seems the conclusion of the author is justified. However, the manuscript appears to show a lack of data. Only 4 tomograms are done for the porcine microtubules. Increasing the data number would make the manuscript statistically convincing. In addition, having the same transition with the missing wedge orientation randomly from different subtomograms will allow a better average of transition without the missing wedge artifact.

Another thing that I found lacking is the mapping of the transition region/alignment in the raw data. However, it is not easy for me or the reader to understand how each segment is oriented relative to each other apart from the simplified seam diagrams in the figures, and also the orientation of the seam corresponding to the missing wedge in the average. With these improvements, I think the conclusion of the manuscript will be better justified.

Reviewer #1 (Public Review):

In this manuscript, Garratt and collaborators describe exposure to male vs. female olfactory cues as an intervention modulating mouse lifespan. They use urine exposure (early life) and soiled bedding exposure (after weaning) to determine the impact of sex-specific olfactory cues on lifespan. They use wild-type and Gnao1 neural-specific mutants to determine potential dependency on vomeronasal function as well. They also recorded information about body temperature, body weight, glucose levels, grip strength, and balance. Overall, they identify female-olfactory cues as increasing the lifespan of female but not male mice, regardless of genotype.

This study reports an intriguing new intervention with an impact on mouse lifespan (with a female-specific effect). However, some of the data with negative results are not plotted/shown, and although all experiments were performed in wild-type and mutant background, all the shown data is pooled and not split by genotype. Overall, this study will be valuable to the field provided that a few analytical points are addressed for clarity and reproducibility and that all methodological details are included.

Joint Public Review:

In this work Malis et al introduce a novel spin-labeling MRI sequence to measure cerebrospinal fluid (CSF) outflow. The glymphatic system is of growing interest in a range of diseases, but few studies have been conducted in humans due to the requirement for and invasiveness of contrast injections. By labeling one hemisphere of the brain the authors attempt to assess outflow through the superior sagittal sinus (SSS), one of the major drainage pathways for CSF, signal changes across time were assessed to extract commonly used metrics. Additionally, correlations with age are explored in their cohort of healthy volunteers. The authors report the movement of labeled CSF from the subarachnoid space to the dura mater, parasagittal dura, and ultimately SSS, evidence of leakage from the subarachnoid space to the SSS, and decreases in CSF outflow metrics with older age.

1. I don't think that the description of Parasagittal dura in figure 1 is correct. There is no anatomical structure at the top of SSS that is known as PSD. The location of the lymphatic structures is also incorrect. Please review "Anatomic details of intradural channels in the parasagittal dura: a possible pathway for flow of cerebrospinal fluid" Neurosurgery 1996 Fox at al. There is usually no obvious tissue between the upper wall of the SSS and the calvarium, which can also be seen in the authors' fig 2A and 2B. All of the tissues located lateral to the SSS are known as PSD. Also, the SSS wall is not as thick as the authors stated and is known as PSD in this region. For this reason, the authors need to revise Fig 1 and it should be changed to PSD in the areas referred to as the SSS wall in the article.

2. The authors described tagged CSF in two pathways: from dura mater to PSD and SAS into the SSS and directly from SAS to SSS. Flow from dura mater to PSD and SAS in the main and supplement cannot be seen. Only a flow from PSD to SSS can be seen. Also, regular dura cannot carry flow-collagen-rich fibrous tissue, except parasagittal dura. There is no flow from dura to the CSF in the figures.

3. The authors have conducted many tests to prevent venous contamination. However, measurements were made based on SSS flow rates in all tests. Small parenchymal venous structures, and small cortical-SAS veins might be tagged due to different flow patterns and T2- Relaxation times.

4. The rate of CSF formation in humans is 0.3 - 0.4 ml min-1. ( Brinker et al 2014. Fluids Barriers CNS). We can assume that the absorption rate is also similar to the CSF formation for the entire system brain and Spine. Therefore, the absorption rate of this very small amount of CSF by SSS is very low in seconds. It is hard to detect by MR and especially CSF flow from the PSD to SSS. The authors concluded that using this technique the rate averaged less than a couple of seconds, rather than on the order of hours or days as previously reported with the use of intrathecal administration of GBCA (Ringstad et al., 2020).

5. Overall, I think that the CSF flow from the PSD to the CNS described by the authors - the CSF flow, might be the venous flow that drains into the SSS slowly, predominantly in the rich venous channels, venous lacunae, and previously described channels in the PSD. Additional explanations are needed.

6. The study is generally well described and to the best of my knowledge an innovative approach. The findings are broadly consistent with what might be expected from the literature and the authors make a good argument in support of their findings. However, the lack of validation is a major limitation of the presented work. In introducing a novel technique a comparison with an existing approach, such as Gd enhanced contrast techniques, or phase contrast would have been expected. Several considerations could have been mentioned/addressed in more detail e.g. what effect labeling efficiency, tortuosity of vessels, lack of gating, the effectiveness of the intensity thresholding to remove the signal from blood, etc may have on the quantification, etc. Without a more thorough validation, it is difficult to evaluate the findings. While scans were conducted on two volunteers to assess reproducibility this is a very small sample and it is notable that scans were conducted consecutively, which might be expected to reduce variance relative to scans further apart e.g. on different dates, scanned by a different operator and no information is provided on how the two scans were positioned (i.e. separately vs copied from the first to the second scan), some metrics showed large percentage differences, which were more pronounced in one subject than the other. Without further data, it is difficult to interpret the reproducibility results. No assessment of the effect of physiological parameters e.g. breathing, cardiac pulsations, or factors affecting glymphatic clearance e.g. amount of sleep the evening before was given.

7. Given these limitations it is hard to adequately assess the likely impact or utility. In recent years several groups have published work e.g. doi.org/10.1038/s41467-020-16002-4 , doi.org/10.1016/j.neuroimage.2021.118755 assessing the blood-CSF barrier. However, previous work has generally focused on larger structures, and by labeling in the oblique-sagittal plane it is unclear how drainage and blood flow rates may affect the presented values here.

8. Some validation data would greatly increase the value of the reported work. I would therefore encourage the authors to consider acquiring some additional datasets to compare measures of CSF draining against another method e.g. 2-D or 4-D phase contrast, or Gd-based contrast-enhanced techniques. Some additional points to consider are noted below.

8. Abstract

CSF outflow may also be imaged with phase contrast MRI (albeit in a limited way).<br /> Demographics would fit better in Results, breakdown could be given for the young and old groups i.e. n, ages, sex.<br /> Conclusion - unless further validation can be provided I think some of the claims should be toned down.

9. Introduction

The authors emphasise the role of Nedergaard, however, there was some relevant earlier work (e.g. Rennels et al, PMID: 2396537).

10. Methods

It would be more conventional to summarise the volunteer characteristics in the Results.<br /> Given the age difference between the two groups, and the fact that for conventional ASL we know of differences in labelling efficiency and the need for a different post-labelling duration in more elderly patients how did the authors account for this?<br /> More broadly what would the effect of differences in labeling efficiency be, given the labeling plane is unlikely to be perpendicular to the draining vessels?<br /> While the authors mention circadian effects there is no mention of controlling for other factors before the scan e.g. caffeine consumption, smoking, etc.<br /> Various mechanisms have been hypothesised to drive glymphatic pulsations. Assessing how physiological signals correlated with the flow may have been a useful proof of concept. Why was it not considered necessary to use a gated acquisition? Did the authors consider the potential impact of respiratory and cardiac pulsations on their measurements?<br /> ROI segmentation - manually selected by two raters, was this done independently and blinded? How were consensus ROIs agreed?<br /> Intensity values outwith MEAN +/- 2 SD were excluded from further analyses. This is justified as removing pulsatile blood. However, was this done independently for tag-on and tag-off? Does this mean slight differences were present in the number of voxels between the two?<br /> The starting points and parameter ranges are given in Eq'n 3, how were the ranges defined? Was there a reason for constraining the fit to positive values only, is there a risk of bias from this?<br /> While the main results appear to have a reasonable sample size n=2 for the reproducibility analysis is very limited. Additional datasets would be useful in properly interpreting the results.

11. Results<br /> While the authors have taken some measures to reduce potential contamination from blood I would be concerned about the risk of surface vessels affecting the signal, and there does not seem to be an evaluation of how effective their measures are.<br /> The labeling pulse is applied in the oblique sagittal orientation, but in tandem with differing rates of blood flow and CSF drainage from the labeling plane does that not risk circulating flow from other slices potentially affecting the values?<br /> Figure 4. The authors focus on the parasagittal dura, but in both the subtraction image and panel C showing different slices at TI=1250 ms some movement appears visible in the opposing hemisphere. Similarly in S2 If the signal does represent CSF movement then this seems counterintuitive and should be explained.<br /> In Figures 4 and 5 the angulation of the TIME-SLIP tag pulse seems quite different. What procedure was used to standardise this, and what effect may this have on the results?

12. Discussion<br /> Phrasing error 'which will be assessed in future studies'.<br /> I would suggest that some of the claims of novelty be moderated e.g. 'may facilitate establishment of normative values for CSF outflow' seems a stretch given multiple pathways exist and this is only considered one.<br /> More consideration should be given to some of the points mentioned in the results. The lack of validation should be properly discussed.

Reviewer #1 (Public Review):

The authors use what is potentially a novel method for bootstrapping sequence data to evaluate the extent to which SARS-CoV-2 transmissions occurred between regions of the world, between France and other European countries, and between some distinct regions within France. Data from the first two waves of SARS-CoV-2 in Europe were considered, from 2020 into January 2021. The paper provides more detail about the specific spread of the virus around Europe, specifically within France, than other work in this area of which I am aware.

An interesting facet of the methodology used is the downsampling of sequence data, generating multiple bootstraps each of around 500-1000 sequences and conducting analysis on each one. This has the strength of sampling, in total, a large number of sequences, while reducing the overall computational cost of analysis on a database that contains in total several hundred thousand sequences. A question I had about the results concerns the extent of downsampling versus the rate of viral migration: If between-country movements are rapid, a reduced sample could be misleading, for example characterising a transmission path from A to B to C as being from A to C by virtue of missing data. I acknowledge that this would be a problem with any phylogeographic analysis relying on limited data. However, in this case, how does the rate of migration between locations compare to the length of time between samples in the reduced trees? Along these lines, I was unclear to what extent the reported proportions of intra- versus inter-regional transmissions (e.g. line 223) would be vulnerable to sampling effects.

A further question around the methodology was the use of an artificially high fixed clock rate in the phylogenetic analysis so as to date the tree in an unbiased way. Although I understood that the stated action led to the required results, given the time available for review I was unable to figure out why this should be so. Is this an artefact of under-sampling, or of approximations made in the phylogenetic inference? Is this a well-known phenomenon in phylogenetic inference?

The value of this kind of research is highlighted in the paper, in that genomic data can be used to assess and guide public health measures (line 64). This work elucidates several facts about the geographical spread of SARS-CoV-2 within France and between European countries. The more clearly these facts can be translated into improved or more considered public health action, through the evaluation of previous policy actions, or through the explication of how future actions could lead to improved outcomes, the more this work will have a profound and ongoing impact.

Reviewer #1 (Public Review):

Ge et. al., examined sodium-glucose cotransporter-2 inhibitors (SGLT2i) in Alport syndrome (AS), and demonstrate that it was beneficial in AS through reduced lipotoxicity in podocytes as a key mechanism of action. The SGLT2i empagliflozin has been previously shown to have positive effects on hyperglycemia control, as well as on cardiovascular and renal outcomes of type II diabetes mellitus through tubuloglomerular feedback, but its effect on glomerular diseases such as AS are unknown to date. The authors have previously identified that cholesterol efflux in podocytes plays a critical pathogenic role in a diabetic kidney disease setting. The evidence that authors provide in favor of their hypothesis in a disease of non-metabolic origin such as AS, was supported as the SGLT2i was effective in reducing the deleterious effects of lipotoxicity in podocytes, ameliorated glomerular injury and proteinuria, and extending the life span of Col4a3 knockout mice. They further show that empagliflozin treatment mitigated AS podocytes from cell death through apoptosis, but did not impact the cell's cytotoxicity. These results support the notion that empagliflozin affects the regulation of important metabolic switch in mouse kidneys, perhaps through decreasing lipid accumulation in podocytes.

However, the authors solely rely on one IHC staining image of a human biopsy to demonstrate SGLT2 expression in podocytes in vivo. Although the authors have done several experiments which greatly increase the confidence in their findings that empagliflozin is beneficial in AS and would have clinical significance, their data does not rule out the possibility that empagliflozin has beneficial effects through the other glomerular cells in AS, or limited to impacting lipids in podocytes in AS.

Reviewer #1 (Public Review):

This manuscript describes the generation and characterisation of a mouse knockout model of Cep78, which codes for a centrosomal protein previously implicated in cone-rod dystrophy (CRD) and hearing loss in humans. Previous work in cultured mammalian cells (including patient fibroblasts) also indicated roles for CEP78 in primary cilium assembly and length control, but so far no animal models for CEP78 were described. Here, the authors first use CRISPR/Cas9 to knock out Cep78 in the mouse and convincingly demonstrate loss of CEP78 protein in lysates of retina and testis of Cep78-/- animals. Next, by careful phenotypic analysis, the authors demonstrate significant defects in photoreceptor structure and function in these mutant animals, which become more severe over a 9 (or 18) month period. Specifically, TEM analysis demonstrates ultrastructural defects of the connecting cilium and photoreceptor outer segments in the Cep78 mutants, which is in line with previously reported roles for CEP78 in CRD and in regulating primary cilia assembly in humans. In addition to a CRD-like phenotype, the authors also convincingly show that male Cep78-/- animals are infertile and exhibit severe defects in spermatogenesis, sperm flagella structure and manchette formation (MMAF phenotype). Furthermore, the authors provide evidence for an MMAF phenotype from a male individual carrying a previously reported CEP78 c.1629-2A>G mutation, substantiating that CEP78 is required for sperm development and function in mammals and supporting previously published work (Ascari et al. 2020).

Finally, to identify the underlying molecular mechanism by which CEP78 loss causes MMAF, the authors perform some biochemical analyses, which suggest that CEP78 physically interacts with IFT20 and TTC21A (an ortholog of Chlamydomonas IFT139) and might regulate their stability. The authors conclude that CEP78 directly binds IFT20 and TTC21A in a trimeric complex and that disruption of this complex underlies the MMAF phenotype observed in Cep78-/- male mice. However, this conclusion is not fully justified by the data provided, and the mechanism by which CEP78 affects spermatogenesis therefore remains to be clarified.

Specific strengths are weaknesses of the manuscript are listed below.

Strengths:

Overall, the phenotypic characterisation of the Cep78-/- animals appears convincing and provides new evidence supporting that CEP78 plays an important role in the development and function of photoreceptors and sperm cells in vertebrates.

Weaknesses:

1) The immunoprecipitation experiments of mouse testis extracts that were used for the mass spectrometry analysis in Table S4 were performed with an antibody against endogenous CEP78 (although antibody details are missing). One caveat with this approach is that the antibody might block binding of CEP78 to some of its interactors, e.g. if the epitope recognized by the antibody is located within one or more interactor binding sites in CEP78. This could explain why the authors did not identify some of the previously identified CEP78 interactors in their IP analysis, such as CEP76 and the EDD-DYRK2-DDB1-VprBP complex (Hossain et al. 2017) as well as CEP350 (Goncalves et al. 2021).

2) Figure 7A-D and page 18-25: based on IPs performed on cell or tissue lysates the authors conclude that CEP78 directly binds IFT20 and TTC21A in a "trimeric complex". However, this conclusion is not justified by the data provided, nor by the previous studies that the authors are referring to (Liu et al. 2019 and Zhang et al. 2016). The reported interactions might just as well be indirect. Indeed, IFT20 is a known component of the IFT-B2 complex (Taschner et al., 2016) whereas TTC21A (IFT139) is part of the IFT-A complex, which suggests that they may interact indirectly. In addition, the IPs shown in Figure 7A-D are lacking negative controls that do not coIP with CEP78/IFT20/TTC21A. It is important to include such controls, especially since IFT20 and CEP78 are rich in coiled coils that tend to interact non-specifically with other proteins.

3) In Figure 7D, the input blots show similar levels of TTC21A and IFT20 in control and Cep78-/- mouse testicular tissue. This is in contrast to panels E-G in the same figure where TTC21A and IFT20 levels look reduced in the mutant. Please explain this discrepancy.

4) The efficiency of the siRNA knockdown shown in 7J-M was only assessed by qPCR (Figure S4), but this does not necessarily mean the corresponding proteins were depleted. Western blot analysis needs to be performed to show depletion at the protein level. Furthermore, it would be desirable with rescue experiments to validate the specificity of the siRNAs used.

5) Figure 7I: the resolution of the IFM is not very high and certainly not sufficient to demonstrate that CEP78, IFT20 and TTC21A co-localize to the same region on the centrosome, which one would have expected if they directly interact.

6) It is not really clear what information the authors seek to obtain from the global proteomic analysis of elongating spermatids shown in Figure 3N, O and Tables S2 and S3. Also, in Table S2, why are the numbers for CEP78 in columns P, Q and R so high when Cep78 is knocked out in these spermatid lysates? Please clarify.

7) Figure 1F and Figure 4K: the data needs to be quantified.

8) Figure 2A: It is difficult to see a difference in connecting cilium length in control and Cep78-/- mutant retinas based on the images shown here.

Reviewer #1 (Public Review):

The idea that a passive living being can improve the wind dispersal of its seeds by passively changing their drag is enticing. The manuscript shows that high wind events in Scotland are inversely correlated with the ambient humidity. The dandelion pappus morphs with the ambient humidity, being more open in dry conditions, which is associated with stronger wind events. This passive morphing of the shape of the pappi thus leads to a dispersal of the seeds further away from their origin.

The analysis and discussion in the paper is focused on "distance", i.e., how far the pappus will fly. Could the notion of time be relevant too? In wet conditions, perhaps it's better for a seed to hit the ground quickly and start germinating, whereas if its dry, staying up in the air for longer to travel farther might be a better strategy.

Reviewer #1 (Public Review):

This paper shows that sibling neurons in the zebrafish spinal cord have different inputs and outputs, and do not show interconnectivity - somewhat surprising considering their very similar development. The differences in sibling neuron connectivity are strongly correlated with the level of Notch signaling, suggesting that Notch signaling regulates circuit assembly.

Reviewer #1 (Public Review):

3' UTR-derived sRNAs are increasingly recognized as post-transcriptional regulators of diverse bacterial processes. While initially assumed to be highly specific regulators with single targets, global RNA interactome approaches challenge this view by suggesting 3'-derived sRNAs can bind and regulate multiple target mRNAs, reminiscent of the regulons of intergenic sRNAs (PMID: 35388892). In Enterobacteriaceae, GlnZ is an sRNA derived from the 3' UTR of the glnA mRNA (PMID: 15718303) that encodes glutamine synthetase. However, the role and function of GlnZ have not previously been determined. In the present work, the authors set out to investigate GlnZ in E. coli and Salmonella enterica. In doing so, they uncover the mechanism of GlnZ biogenesis, namely RNase E-mediated release from the parental glnA mRNA. Additionally, they identify several GlnZ targets (involved in carbon/nitrogen metabolism), some of which are conserved between E.coli and Salmonella while others are species-specific. Downstream mutational characterization of sRNA variants and experiments with target reporter constructs allow them to map sRNA-target interaction sites at nucleotide resolution. Together, their findings further support the idea that 3'-derived sRNAs, too, can act as more global post-transcriptional regulators with multiple direct targets.

This is a thoroughly conducted study and both, important and timely. As the corresponding author points out in the cover letter, this is an instance of similar findings being simultaneously reported by more than one group (see preprint from Gisela Storz' lab: https://www.biorxiv.org/content/10.1101/2022.04.01.486790v1). The results of these two, independently conducted studies largely agree and complement one another.

Reviewer #1 (Public Review):

In this article, Susswein and colleagues use SafeGraph mobile device location data to characterize seasonal trends in indoor activity in the United States at the county level with relevance for respiratory disease transmission. They find substantial variation in indoor activity over the course of the year, ranging from roughly 25% (summer trough) to 200% (winter peak) of the average/baseline indoor activity in each county. Additionally, they identify two main regions with distinct seasonal trends in indoor activity: one in the north, where indoor activity follows a roughly standard sinusoidal trend, and one in the south where indoor activity may feature an additional summer peak. They also identify a third minor region with spring and fall peaks in indoor activity, corresponding to mountainous areas that are hubs for winter tourism. Using a simple mathematical disease transmission model, they demonstrate that using different seasonal forcing terms as inputs can yield substantially different epidemic curves.

This study's main strength is the volume and resolution of the data. Because of this, the authors are able to provide convincing evidence that seasonal variation in indoor activity exists, that it is substantial, and that it varies geographically across the US. Another important strength is the approach that the authors used to identify regions with similar seasonal trends in indoor activity. By using a network community detection algorithm, they were able to avoid making a priori assumptions about the number, size, and geographic connectivity of the regions, allowing them to make better use of the data itself to inform the delineation between regions.

Despite the volume of the underlying dataset, it is geographically limited to the United States and only captures the locations of mobile devices for which their users have opted in to sharing location data. This calls into question the generalizability of the findings to other countries and to other populations within the US that may have reduced access to mobile devices or may be less likely to share location data. The assessment of between-county differences in seasonal indoor activity trends and the assessment of the impact of the COVID-19 pandemic on indoor activity could benefit from greater detail, as they currently rest mainly on visual inspection of the trends.

Overall, the authors have largely achieved their aims of characterizing indoor seasonal activity in the United States at a fine geographic resolution. This work will be immediately useful for the construction of more evidence-based infectious disease transmission models. The authors have made available their estimates of seasonal deviation in indoor activity at the county level, which can be incorporated directly into disease transmission models. Their descriptions are also sufficient for building models that do not incorporate the full county-level detail but nevertheless account for important regional differences in indoor activity across the US.

Reviewer #1 (Public Review):

This paper elucidates the developmental-genetic mechanisms that generate the winged and wingless form (morph) of female pea aphids (Acyrthosiphon pisum). Pea aphids reproduce parthenogenetically generating genetically identical offspring, and so the difference between the winged and wingless morphs is environmentally induced (referred to as a polyphenism). Previous studies have shown that the crowding of mothers is sufficient to induce the winged phenotype in the offspring. The authors develop a technique so that they can reliably generate wing-destined (WD) or wingless-destined (WLD) offspring. This allowed them to examine the early development of WD and WLD offspring during the 1st nymphal stage, before wing development can be observed externally, in the 3rd nymphal stage. They find that the wing primordia are apparent in both WD and WLD 1st instars immediately after birth, but that the primordia degrades 30-36h after birth in WLD nymphs. They then demonstrate that this degeneration is due to autophagy rather than apoptosis, evident through the increased expression of autophagy-related genes in WLD nymphs, but not pro-apoptotic genes. The authors next ask what is responsible for inducing autophagy in the WLD nymphs and so transcriptomics to look for genes that are differentially up- or down-regulated in 1st insar WLD versus WD nymphs. One gene, REPTOR2, is markedly down-regulated in WLD versus WD nymphs. REPTOR is an established target of TOR-signaling, which is in turn an established regulator of autophagy. The authors, therefore, focus on REPTOR2 and show that it has arisen through gene duplication of REPTOR in the A. pisum lineages, that it is differentially upregulated in the thorax of WLD versus WD nymphs, and that knock-down of REPTOR2 both reduces levels of the autophagic protein ATG8 in the wing primordia of 1st instar nymph and increases the proportion of winged offspring. Finally, the authors demonstrate that TOR, which canonically represses RAPTOR, negatively regulates autophagy of the wing primordia and positively regulates the generation of the winged morph.

The strength of this paper is that it cleanly implicates a novel gene, REPTOR2, that has arisen through gene duplication, in the generation of alternative morphs in a polyphenism. The paper also provides compelling evidence that the degeneration of the wing primordia in wingless aphids is through autophagy rather than apoptosis. Further, the paper provides another example of how signaling pathways known to be involved in the generation of reaction norms (continuous phenotypic responses to environmental variation) are also implicated in the generation of polyphenisms (discrete phenotypic responses to environmental variation). The paper uses a reliable and reproducible technique to generate wing and wingless forms of aphids upon which developmental studies can be conducted. The results of the paper are straightforward and convincing. A weakness of the paper is that, while it implicates both RAPTOR2 and TOR-signaling in the generation of winged and wingless morphs, it does not provide a causal link between the two. REPTOR (Repressed by TOR) is known to be a regulator of TOR-signaling in Drosophila, activating transcriptional stress response upon TOR inhibition, and the authors argue that REPTOR2 serves to exert a negative effect of TOR signaling on autophagy initiation in wingless aphids. Nevertheless, their data do not unambiguously show this. Specifically, they do not demonstrate that REPTOR2 is downstream of TOR in the signaling pathway that regulates winglessness.

Reviewer #1 (Public Review):

In a very interesting and technically advanced study, the authors measured the force production of curved protofilaments at depolymerizing mammalian microtubule ends using an optical trap assay that they developed previously for yeast microtubules. They found that the magnesium concentration affects this force production, which they argue based on a theoretical model is due to affecting the length of the protofilament curls, as observed previously by electron microscopy. Comparing with their previous force measurements, they conclude that mammalian microtubules produce smaller force pulses than yeast microtubules due to shorter protofilament curls. This work provides new mechanistic insight into how shrinking microtubules exert forces on cargoes such as for example kinetochores during cell division. The experiments are sophisticated and appear to be of high quality, conclusions are well supported by the data, and language is appropriate when conclusions are drawn from more indirect evidence. Given that the experimental setup differs from the previous optical trap assay (antibody plus tubulin attached to bead versus only antibody attached to bead), a control experiment could be useful with yeast microtubules using the same protocol used in the new variant of the assay, or at least a discussion regarding this issue. One open question may be whether the authors can be sure that measured forces are only due to single depolymerizing protofilaments instead of two or more protofilaments staying laterally attached for a while. How would this affect the interpretation of the data?<br /> This work will be of interest to cell biologists and biophysicists interested in spindle mechanics or generally in filament mechanics.

Reviewer #1 (Public Review):

In this work, Li & Meister provide an extensive description of functional cell types within the posterior-medial part of the mouse superior colliculus, corresponding to the upper lateral visual field. Presenting a battery of visual stimuli to head-fixed wild-type mouse lines, they use calcium imaging and subsequent clustering of functional responses to identify 24 functional cell types. Besides the comprehensive sampling of SC cell types, a major strength of the manuscript in my view is the direct comparison with the previously published in vitro RGC data. Overall, the manuscript and the associated data promise to be a valuable resource to experimentalists and computational neuroscientists comparing visual processing across the major processing stages of the mouse visual system. However, in the current form the manuscript still comes with some limitations. In my view, these are related to some parts, where statistical justifications for the conclusions are still missing, where the findings should be more strongly embedded into the current and past literature, and where more efforts should be made to relate the findings in the different mouse lines to those for the overall population.

Reviewer #1 (Public Review):

The authors have generated a set of seven nanobody tools against two of the largest Drosophila proteins, which are related to vertebrate titin and essential for muscle function. The study of such gigantic proteins is a challenge. They show that each of these nanobodies recognizes their epitope with high affinity (as expected from antibodies), fails to generate a signal after immune-fixation of a mutant for the cognate protein, do not cross-react with each other, and generates a signal in the muscle that makes sense with what one would anticipate for fly titin homologs. In addition, they show that these nanobodies have better penetration and labeling efficiency than conventional antibodies in thick tissues after classical paraformaldehyde fixation. Using these nanobodies, they could deduce the organization of the epitopes in different muscle types and propose a model for Sallimus and Projectin arrangement in muscles, including in larvae which are difficult to label with traditional antibodies due to their impermeable chitin skeleton. Finally, they could fuse the gene encoding one of the nanobodies to the open reading frame of NeonGreen and express the corresponding fusion protein in animals to use the probe in FRAP assays.

The work is very well performed and convincing. However, given its significant redundancy in terms of biological conclusions with the companion study "Nanobodies combined with DNA-PAINT super-resolution reveal a staggered titin nano-architecture in flight muscles" by the same authors, and other published papers, I recommend the authors further prove the use of their nanobodies in live assays. In particular, the authors should test whether they can use the nanobodies to induce protein degradation either permanently or conditionally.

Reviewer #1 (Public Review):

The authors found that the IDR in Cdc15 gets phosphorylated by multiple kinases, Pom1/Shk1/Pck1/Kin1, and the phosphorylation on IDR inhibits the phase separation of the Cdc15 protein. The phosphorylation was demonstrated in the cell as well as in vitro. Moreover, the phosphorylation sites were identified by mass spectrometry. The phospho-regulation of Cdc15 LLPS was demonstrated by in vitro assay using recombinant proteins. The significance of the phosphorylation on contractile actomyosin ring (CAR) was demonstrated by using a cdc15 mutant carrying 31 Ala-substitutions at the phosphorylation sites (cdc15 31A). The CAR assembled comparable to cdc15+, but maturation and contraction of the ring were faster in the cdc15 31A mutant, suggesting the contribution of the phosphorylation for delaying cytokinesis. This could be one of the mechanisms to ensure the completion of chromosome segregation before the cytokinesis. In this paper, the authors showed over-accumulation of type-II myosin regulatory light chain Rlc1 on CAR in the cdc15 31A mutant during the CAR assembly and its contraction. In addition, the kinases for the Cdc15 IDR phosphorylation are identified as polarity kinases, which restrict the assembly of the CAR formation in the middle. Indeed, inhibition of the kinases increases the ratio of septa formation at the cell tip in the mid1 knockout mutant, which lacks a major positive polarity cue during the mitotic phase. However, in this manuscript, this phenotype is not solely explained by the phosphorylation of the cdc15 31A, because the authors did not show the tip septa formation using cdc15 31A.

Overall, the data supports their conclusion, Cdc15 forms LLPS, and the process is inhibited by the phosphorylation of amino acid residues in the IDR in Cdc15 by polarity kinases. It is still unclear whether LLPS formation is a reversible process regulated by the protein kinases. In vitro experiments showed condensate formation by dephosphorylation of Cdc15 IDR but not diffusion of the LLPS by phosphorylation. I wonder if incubation of the kinases and the Cdc15 IDR condensates induces demolition of the LLPS.

The transition of the Cdc15 IDR phosphorylation and LLPS formation through the cell cycle progression is unclear. In asynchronous cells (most of the cells may be in the G2 phase) and nda3 or cps1 mutants, Cdc15 was still highly phosphorylated. This indicates that the Cdc15 is phosphorylated and the LLPS formation is inhibited throughout the cell cycle. The transition of the phosphorylation status for individual residues could be the next challenge for this research. In addition, currently, there is no approach to monitor the LLPS in wild-type cells. Therefore, it is still unclear if LLPS formation is the physiological mechanism regulating cell division in wild-type cells.

Reviewer #1 (Public Review):

This paper tests whether people vary their reliance on episodic memory vs. incremental learning as a function of the uncertainty of the environment. The authors posit that higher uncertainty environments should lead to more reliance on episodic memory, and they find evidence for this effect across several kinds of analyses and across two independent samples.

The paper is beautifully written and motivated, and the results and figures are clear and compelling. The replication in an independent sample is especially useful. I think this will be an important paper of interest to a broad group of learning, memory, and decision-making researchers. I have only two points of concern about the interpretation of the results:

1. My main concern regards the indirect indicator of participants' use of episodic memory on a given trial. The authors assume that episodic memory is used if the value of the chosen object (as determined by its value the last time it was presented) does not match the current value of the deck it is presented in. They find that these mismatch choices happen more often in the high-volatility environment. But if participants simply choose in a more noisy/exploratory way in the high volatility environment, I believe that would also result in more mismatched judgments. What proportion of the trials labeled as episodic should we expect to be a result of noise or exploration? It seems conceivable that a judgment to explore could take longer, and result in the observed RT effects. Perhaps it could be useful to match up putative episodic trials with later recognition memory for those particular items. The across-subjects correlations are an indirect version of this, but could potentially be subject to a related concern if participants who explore more (and are then judged as more episodic) also simply have a better memory.

2. The paper is framed as tapping into a trade-off between the use of episodic memory vs. incremental learning, but it is not clear why participants would not use episodic memory in this particular task setup whenever it is available to them. The authors mention that there is "computational expense" to episodic memory, but retrieval of an already-established strong episodic memory could be quite effortless and even automatic. Why not always use it, since it is guaranteed in this task to be a better source of information for the decision? If it is true that RT is higher when using episodic memory, that is helpful toward establishing the trade-off, so this links to the concern above about how confident we can be about the use of episodic memory in particular trials.

Reviewer #1 (Public Review):

Trudel and colleagues aimed to uncover the neural mechanisms of estimating the reliability of the information from social agents and non-social objects. By combining functional MRI with a behavioural experiment and computational modelling, they demonstrated that learning from social sources is more accurate and robust compared with that from non-social sources. Furthermore, dmPFC and pTPJ were found to track the estimated reliability of the social agents (as opposed to the non-social objects).

The strength of this study is to devise a task consisting of the two experimental conditions that were matched in their statistical properties and only differed in their framing (social vs. non-social). The novel experimental task allows researchers to directly compare the learning from social and non-social sources, which is a prominent contribution of the present study to social decision neuroscience.

One of the major weaknesses is the lack of a clear description about the conceptual novelty. Learning about the reliability/expertise of social and non-social agents has been of considerable concern in social neuroscience (e.g., Boorman et al., Neuron 2013; and Wittmann et al., Neuron 2016). The authors could do a better job in clarifying the novelty of the study beyond the previous literature.

Another weakness is the lack of justifications of the behavioural data analyses. It is difficult for me to understand why 'performance matching' is suitable for an index of learning accuracy. I understand the optimal participant would adjust the interval size with respect to the estimated reliability of the advisor (i.e., angular error); however, I am wondering if the optimal strategy for participants is to exactly match the interval size with the angular error. Furthermore, the definitions of 'confidence adjustment across trials' and 'learning index' look arbitrary.

As the authors assumed simple Bayesian learning for the estimation of reliability in this study, the degree/speed of the learning should be examined with reference to the distance between the posterior and prior belief in the optimal Bayesian inference.

Reviewer #1 (Public Review):

The Voeltz lab in previous work has established a physical connection between the endoplasmic reticulum and mitochondria as an organizing principle in mitochondrial fission and fusion. A key cytoplasmic protein, Drp1 and a mitochondrial outer membrane protein, Mfn1, are known to localize to nodes of interaction between the ER and mitochondria but until now, no ER membrane proteins required in this process have been described. Using a Turbo ID fused to Mfn1, the authors have identified a known ER membrane protein that interacts with functional Mfn. This ER membrane protein, which they call Aphyd1, contains putative acyl hydrolase and acyltransferase domains. Further, in compelling work combining high resolution fluorescence microscopy and gene function analysis they have described the role of this protein in the recruitment of Drp1 and Mnn1 to nodes of interaction between the ER and mitochondria and in the sequential processes of mitochondrial constriction, fission and fusion. The work is clear and nicely quantitative and adds a new molecular element to the important question of how the ER serves to organize the division and fusion of mitochondria.

Reviewer #1 (Public Review):

This manuscript clearly demonstrates that murine malaria infection with Plasmodium chabaudi impairs B cells' interaction with T cells, rather than DCs interaction with T cells. The authors elegantly showed that DCs were activated, capable of acquiring antigens and priming T cells during P. chabaudi infection. B cells are the main APC to capture particulate antigens such as infected RBC (iRBC), while DCs preferentially take up soluble antigens. This study is important to understand how ongoing infections such as malaria may negatively affect heterologous immunizations.

Overall, the experimental designs are straightforward, and the manuscript is well-written. However, there were several limitations in this study.

Specific comments:

1. The mechanism of how the prior capture of iRBC by B cells lead to the impairment of B-T interaction was not understood. It is unclear whether the impairment of B-T cell interaction is due to direct BCR interaction with iRBC, or an indirect response to extrinsic factors induced by malaria infection.

2. Would malaria infection in MD4 mouse that carries transgenic BCR that does not recognize malaria parasite impair subsequent B cell response to HEL immunization? This may clarify whether the impairment of subsequent B cell response is BCR-specific. If malaria impairs subsequent B cell response to HEL in MD4 mouse, it might suggest that other cell types and B cell-extrinsic factors might be involved in causing the impaired B cell responses, instead of malaria affecting B cells directly.

3. MD4 mice were mentioned in the Methods in vitro RBC binding, although none of the figures described the usage of MD4 mice. This experiment data might be important to show whether RBC binding to B cells is mediated through BCR.

4. Does P. chabaudi infection have any effects on B cell uptake of subsequent antigens, such as soluble antigen PE or particulate antigen CFSE-labeled P. yoelii iRBC?

5. Is this phenomenon specific to malaria infection? Does malaria-irrelevant particulate immunization affect T-B interaction of subsequent heterologous immunization?

6. Despite the impaired Tfh and GC 8 days after immunization following malaria infection, Fig. 5F showed GP-specific IgG eventually increased to the same level as the uninfected immunized mice on day 23. Did the authors check whether these mice had a delayed Tfh and GC response that eventually increase on day 23? Are these antibody responses derived from GC, or GC-independent response?

7. Does recovery from malaria infection by antimalarial treatment rescue the B cell response to subsequent heterologous immunization?

8. Fig. 1C shows more nRBC was taken up than iRBC in B cells, but Line 142 states that "B cells bound significantly more iRBC than nRBC. Is there a mistake in the figure arrangement? Why do B cells take up for naïve RBC than iRBC?

9. Fig. S1 C and D are confusing. CD45.1+ CD45.2+ mouse did not receive labeled iRBC, but why iRBC was detected as much as 40% in the spleen of this naïve mouse?

Reviewer #1 (Public Review):

This paper details the creation and data behind the website http://pandemics.okstate.edu/covid19/. The authors attempt to explore if there is a cause and effect between the detection of unusually increased mutation activity in the genomic surveillance databases and subsequent near-term surges in SARS-CoV-2 case numbers.

Overall the premise is interesting as other than following case numbers reported to health authorities and observing what is happening in another country, there is no reliable way to predict when a surge is going to occur. Unfortunately, the data demonstrate that there was no reliable metric that could be used to predict surge events. Interestingly, the website has issued a "surge alert" currently for the month of September. It will be interesting to observe whether their model indeed has predictive power or whether the current analysis is merely coincidental with the surges but not necessarily predictive of them.

Reviewer #1 (Public Review):

This interesting manuscript by Goto and Miyamichi analyses calcium dynamics in the kisspeptin neurons of the arcuate nucleus of the hypothalamus during the estrous cycle and during reproductive aging in female mice. In particular, the authors succeed in tracking arcuate kisspeptin calcium activity in the same mice over 10 months, which is quite impressive and brings highly valuable information. The authors demonstrate that the frequency and the amplitude and waveforms of individual synchronous episodes of arcuate kisspeptin neuronal activation vary across the estrous cycle. Unexpectedly, however, aging does not appear to alter markedly calcium dynamics in these neurons.

Reviewer #1 (Public Review):

Ferroportin (Fpn) is a Fe2+/H+ antiporter that extrudes Fe2+ from cells and is important for iron homeostasis. Using a combination of proteoliposome assays, mutagenesis and structural studies by cryo EM the authors are aiming to demonstrate that the Fpn-transporter is also capable of Ca2+ influx, indicating a novel route for Ca2+ entry into cells.

Strengthens: The paper combines a number of different methods to robustly demonstrate the interaction of Ca2+ with the iron transporter and to show translocation of Ca2+ is not pH dependent.

Weaknesses: Fpn uses proton-gradient to drive Fe2+ efflux. The proposal is that the antiporter can also passively uptake Ca2+. This means that after Ca2+ release on the inside, Fpn would need to spontaneously rest to the outside again, which it has not evolved to do in the absence of Fe2+. To provide further support for Ca2+ uptake it is important to show that there are mutiple turnovers of the transporter, i.e., more kinetic information is needed,

The impact of this paper is the demonstration that transporters (exchangers) can also operate as facilitative transporters for other substrates. The study also implies that Ca2+ can enter cells by this pathway, but if so the physiological context of this entry route needs further investigation and/or justification.

Reviewer #1 (Public Review):

Here the authors set out to disentangle neural responses to acoustic and linguistic aspects of speech. Participants heard a short story, which could be in a language they understood or did not (French vs. Dutch stories, presented to Dutch listeners). Additional predictors included a combination of acoustic and linguistic factors: Acoustic, Phoneme Onsets, Phoneme Surprisal, Phoneme Entropy and, Word Frequency. Accuracy of reconstruction of the acoustic amplitude envelope was used as an outcome measure.

The use of continuous speech and the use of comprehended vs. uncomprehended speech are both significant strengths of the approach. Overall the analyses are largely appropriate to answer the questions posed.

The reconstruction accuracies (e.g., R^2 values Figure 1) seem lower perhaps than might be expected - some direct comparisons with prior literature would be welcome here. Specifically, the accuracies in Figure 1A are around .002-.003 whereas the range seen in some other papers is about an order of magnitude or more larger (e.g. Broderick et al. 2019 J Neurosci; Ding and Simon 2013 J Neurosci).

One theoretical point relevant to this and similar studies concerns the use of acoustic envelope reconstruction accuracy as the dependent measure. On the one hand, reconstruction accuracy provides an objective measure of "success", and a satisfying link between stimulus and brain activity. On the other hand, as the authors point out, envelope reconstruction is probably not the primary goal of listeners in a conversation: comprehension is. Some discussion of the implications of envelope reconstruction accuracy might be useful in guiding interpretation of the current work, and importantly, helping the field as a whole grapple with this issue.

Overall, the results support the authors' conclusions that acoustic edges and phoneme features are treated differently depending on whether a listener comprehends the language being spoken. In particular, phoneme features contribute to a greater degree when language is comprehended, whereas acoustic edges contribute similarly regardless of comprehension. These findings are important in part because of prior work suggesting that acoustic edges are critically important for "chunking" continuous speech into linguistic units; the current results re-center language units (phonemes) as critical to comprehension.

Reviewer #1 (Public Review):

In this genetic and imaging based analysis of stem-cell maintenance and organ initiation, two phases important for continued production of shoot organs in plants, the authors tested whether SHR and targets/partners (SCR, SCL23, JWD) provide the circuitry to maintain stem cell pool and contribute to the production of lateral organs. Finding that these factors are indeed expressed in and required for SAM activities, and furthermore, behaviors of SHR and SCR in the root are recapitulated in the meristem, including mobility of SHR (here to epidermis from internal layers), activation of SCR by SHR, and "trapping" of SHR movement by complexing with SCR. Strengths include high quality imaging of reporters and FRET-FLIM measurement to assess in vivo complex formation. The analysis is then extended to link SHR and SCR to shoot-specific factors and auxin, again by testing expression, genetic dependencies and physical interaction. This is repeated for a number of factors and individually, each is well done experiment. Conclusions about causal relationships are somewhat overstated (for example, the idea that SHR-SCR act through CYCD6 to alter cell division is based on expression patterns, not a functional analysis of cycd6).

In general, there are many high-quality studies included in this paper, and the presentation of imaging data (both the images themselves and quantification of data) is excellent. There is also a lot of data, and while each section was presented in a logical way, connections between sections, and the overarching developmental questions were sparse. Because the authors found that many of the relationships defined in the root were recapitulated in the shoot, the present organization leaves one with somewhat of a sense that little new was learned, and yet, the shoot meristem IS different and there are shoot specific inputs into the core regulatory factors. Rewriting to highlight the different activities (and thus expectation about regulation) could make the finding of the same network more interesting and creating a summary figure that highlights the input of shoot specific signals would bring the unique analysis to the forefront.

Reviewer #1 (Public Review):

In this manuscript, the authors present a suite of statistical models that can identify patterns of somatic single base mutational signatures (SBS), insertions, deletions, and structural rearrangements predictive of DNA damage response (DDR) gene deficiency. A similar approach (HRDetect) has already proved successful in identifying BRCA1/BRCA2 deficiency in breast cancers and other tumor types. To generate their models, Sørensen and colleagues consider over 700 DDR gene deficiencies across more than 6,000 patients enrolled through the Hartwig Medical Foundation and PCAWG studies. The authors also consider the full set of COSMIC SBS reference signatures. The models recapitulate known associations between BRCA1/2,TP53, and CDK12 and mutational patterns but also characterize previously undescribed associations involving ATRX, PTEN, CDKN2A, and SMARCA4. Many of these novel models generalize across different tumour types and also primary and metastatic cancers. The authors also consider negative coefficient features in their models which is worthwhile and present hypotheses supporting how negative features might arise. One of the further strengths of the study is its innovative reuse of large-scale cancer genome data sets leading to predictive models with potential use for clinical intervention and demonstrating the potential of using WGS mutational signatures to guide cancer treatment. Many of the findings and observations presented in the paper have the collateral potential to enhance our understanding of the aetiology of SBS signatures. While I don't think the paper presents any major conceptual and technical advances in terms of methodology, the manuscript is important, interesting, and timely.

Reviewer #1 (Public Review):

This work by Olesen et al provides a clear and thorough examination of participation in organised colorectal cancer screening in Denmark over 2018-2021, with a particular focus on potential impacts of the COVID-19 pandemic. There is also an analysis by population subgroups that offers additional insight into how the pandemic may have affected participation.

The key strength of this manuscript is access to the presumably complete data from Danish Colorectal Cancer Screening Database, with only a small proportion of individuals being excluded for sensible reasons. Combined with a standard statistical analysis, the results are clear and inarguable.

A weakness is that the manuscript is wholly quantitative. Therefore some of the observations made, particularly for population subgroups, cannot be explained within the scope of this manuscript. The authors have provided hypotheses to explain these, but without further research no firm conclusion can be made.

Another weakness is that, due to the nature of screening, no conclusions can be made on the health impact of COVID-related changes to screening. It is unclear whether small reductions in screening and colonoscopy follow-ups is likely to lead to additional cancers, or later-stage diagnoses.

This manuscript, together with similar analyses in other settings worldwide, provides both an overview of how the COVID-19 pandemic has affected screening, but also potentially guidance for how future potential disruptions should be managed in the cancer screening space. By considering the potential downstream impact of changes to screening, and balancing these against potential harms and resource restrictions associated with (eg) attending health services during a pandemic, policymakers can make informed decisions to manage future shutdowns.

Reviewer #1 (Public Review):

Liu, et al. describe an essential role for prostaglandin signaling during neonephrogenesis in the zebrafish kidney following acute kidney injury. They identified that renal interstitial cells are the source of prostaglandin, and demonstrated which components of the prostaglandin biosynthesis pathway as well as which receptor is involved in the signaling. Further, they determined the mechanism by which PGE2 stimulates the proliferation of renal progenitor cells which make the new nephrons, namely to regulate B-catenin levels. This work is systematic, thorough, and well-controlled. The applications of these findings may have a profound impact on the formulation of regenerative medicine treatments for kidney disease.

Reviewer #1 (Public Review):

Tan and colleagues examine the role of additional genetic removal of Munc13 in murine cultured synapses deficient for RIM and ELKS. They utilize a comprehensive set of morphological, ultrastructural and functional experiments to conclude that Munc13 is a nonredundent factor in synaptic vesicle priming. In addition, the results contribute to the ongoing discussion whether synaptic vesicle docking represents synaptic vesicle priming.

The main strength of the work is the use of a sophisticated genetic model, and the quality of the performed experiments. The results strongly support the conclusion of the nonredundant role of Munc13 in synaptic vesicle priming.

The weakness of the paper is that the findings from these study has limited impact, as the genetic and functional interaction of Munc13 and RIM has been extensively analyzed on a qualitative and quantitative level (Kaeser 2021, Zarebidaki 2020). While this paper benefits from the additional deletion of ELKS, the specific contribution of ELKS in comparison to the older studies is not in the focus of the study.

While the genetic removal of all 6 genes involved clearly require the use of conditional KO mice, the resulting outcome of the experimental design is a hypomorphic phenotype, as neurotransmitter release is still detected. This complicates the interpretation of the findings and weakens the strength of the conclusions.

Reviewer #1 (Public Review):

This manuscript addressed an important question regarding an obstacle to hair reprogramming in older mice that is not present in newborn mice. The conclusions were justified by experimental results. However, the scientific novelty is limited, and there is a lack of functional characterization of the newly formed hair cells.

There are several strengths of this study: 1. It addressed a significant question as hearing loss is an important public health issue. 2. Well-designed genetic approaches. 3. Experiments were well designed and justified. 4. Experimental results are convincing. 5. Conclusions were well justified by experimental results. There are also several weaknesses:

1. The scientific novelty is limited. It is known that overexpression of several transcription factors simultaneously can reprogram hair cells and non-hair cells with hair cell characteristics.

2. Transcription factor Atoh1 and downstream GFI1and POU4F3 have been used to reprogram embryonic stem cells and chick otic epithelial cells in vitro to cells expressing several hair cell genes and displaying key hair cell features.

3. There is no functional characterization of newly reprogramed hair cells in adult mice although FM 1-43 dye was used for characterizing reprogramed hair cells in neonatal mice.

4. It is not understood why the changes in transduction channel protein expression were not highlighted in gene analysis.

5. It will be nice if hair cell-like electrophysiological properties can be found in newly reprogramed hair cells.

Reviewer #1 (Public Review):

The authors investigate whether and how PFA fixation affects the structures formed by some proteins that undergo LLPS. They do that by over-expressing a number of constructs in cells and imaging them by live cell fluorescence microscopy, after which they fix the cells and image the same cell after fixation. Their results clearly show that, for different proteins and with different tags, there is a non-systematic alteration of the LLPS structures.

In parallel to this experimental work, the authors also present and analyze a dynamic computational model in which they investigate how different fixation rates for those proteins inside and outside the condensates can lead to alterations in the overall condensate organization after fixation. Their model shows that if the fixation rate inside the condensate is different than outside the condensate, and if the dynamics of protein exchange in/out of the condensates are fast enough, fixation artifacts are to be expected.

It remains to be seen whether the alterations in condensate structures after fixation (as seen experimentally) are caused by different fixation rates (as shown computationally).

Overall, this manuscript puts forward an important observation, on how chemical fixation can alter cellular structures, such as those of membraneless organelles.

Reviewer #1 (Public Review):

There are five major claims. First is a replication of lower spatial frequency representation in V4. This is based on two examples shown in Fig. 3. The differences look clear but should be analyzed statistically.

The second claim is that, on the large scale of the visual field representation in V2 and V4, spatial frequency is mapped from high to low going from fovea to periphery, here estimated as lateral to medial, as in V1. The analysis for this is to plot the geometric centroids of 6 different spatial frequency band responses, for orientation contrasts (Fig. 4A) and color contrasts (Fig. 4B), and show that they progress in position from lateral to medial for high to low frequencies. This seems like an unusual analysis that obscures most of the original data concerning the relationship of spatial frequency response profiles to the two-dimensional imaging area. And, the centroids do not show a continuous map, but rather a bunching of points except for the extremes. The additional data are 4 supplemental examples with three, different frequency ranges. These data are not analyzed statistically.

The third and most important claim is that spatial frequency and orientation are mapped orthogonally (and recursively) in V2 and V4, as seen in Fig. 5 and the Fig. 5 supplement. Together these figures present two regions in V4 and two regions in V2. If these are the only analyzed regions, the authors need to specify more clearly how they were selected. Presumably, though, other regions were analyzed, and the authors should present results from all analyzable regions, and use statistical analyses to establish significance.

The fourth claim is that color-sensitive regions in V4 are more associated with low spatial frequencies. The one significant example (the analysis and statistical tests need to be explained), shown in Fig. 6, shows a weak relationship to color for both spatial frequency bands, and the other examples presented in the supplementary are not significant and have even lower absolute relationships. These results, if presented, should be considered inconclusive.

The fifth claim is for stripe-like periodicity of spatial frequency representation in V2, related to color tuning. This is supported by ostention to binary maps of spatial frequency tuning in Fig. 7 and supplement. Establishing this periodicity would require statistical analysis, and in any case, seems impossible since only a sliver of V2 is visible in these brain surface images, so stripes orthogonal to the V1/V2 boundary (i.e. CO stripes) cannot be distinguished from other patterns of spatial frequency tuning. In fact, Fig. 5E and S5I do not appear to have iso-frequency contours biased toward that orientation.

Reviewer #1 (Public Review):

Auwerx et al. have taken a new approach to mine large existing datasets of intermediary molecular data between GWAS and phenotype, with the aim of uncovering novel insight into the molecular mechanisms which lead a GWAS hit to have a phenotypic effect. The authors show that you can get additional insight by integrating multiple omics layers rather than analyzing only a single molecular type, including a handful of specific examples, e.g. that the effect of SNPs in ANKH on calcium are mediated by citrate. Such additional data is necessary because, as the authors' point out, while we have thousands of SNPs with significant impact on phenotypes of interest, we often don't know at all the mechanism, given that the majority of significant SNPs found through GWAS are in non-coding (and often intergenic) regions.

This paper shows how one can mine large existing datasets to better estimate the cellular mechanism of significant, causal SNPs, and the authors have proven that by providing insight into the links between a couple of genes (e.g. FADS2, TMEM258) and metabolite QTLs and consequent phenotypes. There is definitely a need and utility for this, given how few significant SNPs (and even fewer recently-discovered ones) hit parts of the DNA where the causal mechanism is immediately obvious and easily testable through traditional molecular approaches.

I find the paper interesting and it provides useful insight into a still relatively new approach. However, I would be interested in knowing how well this approach scales to the general genetics community: would this method work with a much smaller N (e.g. n = 500)? Being able to make new insights using cohorts of nearly 10,000 patients is great, but the vast majority of molecular studies are at least an order of magnitude smaller. While sequencing and mass spectrometry are becoming exponentially cheaper, the issue of sample size is likely to remain for the foreseeable future due to the challenges and expenses of the initial sample collection.

Reviewer #1 (Public Review):

The paper studied spatial-temporal characteristics of dominant T cell clones in juvenile idiopathic arthritis. The authors found that the composition and functional characteristics of immune infiltrates are strikingly similar between joints within one patient, and observed a strong overlap between dominant T cell clones, especially Treg. Moreover, in localized autoimmune disease there is auto-antigen driven expansion of both T effector and Treg clones, that are highly persistent and are (re)circulating.

Reviewer #1 (Public Review):

Tafenoquine is an important 8-aminoquinoline antimalarial, mostly aimed at the management of Plasmodium vivax malaria. Through the retrospective analysis of several previously performed efficacy trials, the authors aimed to better understand the drugs mechanism of action, while exploring the possibility of improved efficacy through dose increment.

Strengths: robust analysis approaches unlocked three main messages with the potential of improving the clinical practice:<br /> i. P. vivax recurrency is positively associated with tafenoquine terminal half-life and D7 methemoglobin levels.<br /> ii. The methemoglobin levels support the current view that tafenoquine, acts through its metabolites, similar to what is believed for primaquine.<br /> ii. Most importantly, the therapeutic window of tafenoquine is wider than previously considered, allowing the suggestion of a significant increase in dosing, from 300 mg to 450 mg, leading to significantly increased efficacy.

Weaknesses: being a retrospective analysis, the work is limited to the available data. In particular, and as referred by the authors, no drug levels are reported. Additionally, there are some aspects that in my view need more detailed analysis and discussion, in particular, what seems to be a lack of exploration as to the importance (or lack of it) of the patient CYP2D6 status in Tafenoquine T1/2, methemoglobin levels, and overall efficacy. These mild weaknesses do not change the overall conclusions of the study.

Reviewer #1 (Public Review):

This paper is a continuation of other research by this group and represents another step back in time for peptide preservation in eggshells. It is exciting to see Miocene age peptides and that they overlap so completely with both extant ostrich struthiocalcin as well as the previously described Pliocene peptides. The biggest weakness is the lack of tables showing both the de novo peptides as well as those detected by database searching.

Reviewer #1 (Public Review):

The manuscript by Heckman and Doe describes a nice set of experiments that extend previous studies of the plasticity in wiring and growth of a sensory axon terminal (Dbd) and its connection to a partner interneuron (A08a) in Drosophila embryonic/larval ventral nerve cord. The authors confirm and extend prior studies in the lab that showed misrouting of Dbd axons cause changes in its site of connection on medial versus lateral dendrites of A08a. The authors show using misrouting and ablation of Dbd that the site of axonal innervation plays a role in promoting dendrite outgrowth in that specific domain of the A08a dendritic field, suggesting a contact-dependent dendrite growth mechanism regulates early connections. The authors then describe a second mechanism where activity of the Dbd sensory neuron regulates a separate aspect of early connectivity, whereby reduced activity leads to increase A08a dendrite growth globally and increased activity suppresses overall A08a dendrite growth. This study fits with other work in the field on the role of activity in synaptic wiring while highlighting the opposing roles of contact versus activity in establishing early connection patterns. They also identify a brief developmental window where Dbd ablation causes A08a dendritic undergrowth, suggesting an early critical period for contact-dependent dendritic growth modulation, similar to that observed for activity-dependent plasticity. Overall, this is a nice study that provides important advances in our understanding of the plasticity of early neuronal wiring.

Reviewer #1 (Public Review):

This study examines how the COVID-19 pandemic impacted cervical cancer screening participation to invitations sent through the organized cervical cancer screening program of Denmark. I think the results are particularly enlightening in the context of pandemic recovery, as they show that while the short-term participation (90 days) dropped due to public health messaging emphasizing staying at home, the long-term participation (365 days) did not drop; this suggests that women did not completely miss the opportunity to screen during the pandemic, but simply postponed their screening to a later point in time. I think this has implications, especially for modeling the impact of the pandemic on cancer incidence, as many screening models have made the assumption that screenings missed during the pandemic would not be "caught up" later leading to higher cancer incidence in the long term; however, this study suggests that this is not the case and that there is a natural 'catch up' of screening that occurs over time. This is reassuring, as a short delay in cervical cancer screening would not be expected to lead to overly important long-term negative health outcomes.

Particular strengths of this study include the population-based registry covering the whole target population, and the ability to link the data to socioeconomic variables of interest to examine whether there were particular groups of women which were more impacted than others. The models also accounted for seasonal and long-term trends in cancer screening participation, which bolsters the confidence that their results are not the result of trends in cervical screening participation over time and are most likely attributable to the COVID-19 pandemic. However, as the statistical methods do not include an interaction test for the overall effect of each socioeconomic variable, it is not clear whether the differences that are observed between women by age and socioeconomic status are significant.

Reviewer #1 (Public Review):

This study explores the mechanisms responsible for reduced steroidogenesis of adrenocortical cells in a mouse model of systemic inflammation induced by LPS administration. Working from RNA and protein profiling data sets in adrenocortical tissue from LPS-treated mice they report that LPS perturbs the TCA cycle at the level of succinate dehydrogenase B (SDHB) impairing oxidative phosphorylation. Additional studies indicate these events are coupled to increased IL-1β levels which inhibit SDHB expression through DNA methyltransferase-dependent DNA methylation of the SDHB promoter.

In general, these are interesting studies with some novel implications. I do, however, have concerns with some of the author's rather broad conclusions given the limitations of their experimental approach. The paper could be improved by addressing the following points:

1. The limitations of using LPS as the model for systemic inflammation need to be explicitly described.<br /> 2. The initial in vivo findings, which support the proposed metabolic perturbation, are based on descriptive profiling data obtained at one time point following a single dose of LPS. The author's conclusion that the ultimate transcriptional pathway identified hinges critically on knowledge of the time course of this effect following LPS, which is not adequately addressed in the paper. How was this time and dose of LPS established and are there data from different dose and time points?<br /> 3. Related to the point above, the authors data supporting a break in the TCA cycle would be strengthened direct biochemical assessment (metabolic flux analysis) of step kin the TCA cycle process impacted.<br /> 4. The proposed connection of DNMT and IL1 signaling to systemic inflammation and reduced steriodogenesis could be more firmly established by additional studies in adrenal cortical cells lacking these genes.

Reviewer #1 (Public Review):

Neverov and colleagues present a large-scale computational investigation of epistatic interactions between substitutions in the spike protein. The analysis is based on an improved version of their previous approach that has been applied to other organisms to the Influenza A virus. They find several sets of interacting sites that tend to change in concert.

The approach is sensible and the work seems well executed. A systematic investigation of epistatic interactions is important to better understand the constraints and drivers of future SC2 evolution. This work is hence an important contribution to the field and a nice complement to experimental work by Jesse Bloom's group and others.

The authors uncover several groups of residues that seem to change in concert. The identified groups make sense, but further validation and comparison with experimental or other computational approaches would strengthen the conclusions.

Reviewer #1 (Public Review):

Summary:

It is widely known that lesioning the vHP produces anxiogenic effects, and cells in the vHP increase their firing rates in the anxiogenic location. This paper aims to investigate the neural dynamics of the vHP when a non-anxiogenic location changes to an anxiogenic one during spatial navigation. For testing, the authors removed half of the side walls of the elevated linear maze or track in the middle of the session. They reported that cells in the vHP remapped and overrepresented anxiogenic places as the walls were removed. Also, the authors claim that single-cell activities recorded before entering the open portion of the track could be used to predict how far the rat would explore along the open segment of the track.

Strength and weakness:

The experimental paradigm was well-designed to examine the main research question. It is novel in that the authors recorded single cells electrophysiologically in the vHP, as this has been done only in very few studies. However, in the current version of the manuscript, the main argument (tied to Figure 5) is not supported by detailed neural and behavioral data. Specifically, the authors did not provide the basic firing properties of the electrophysiological data to verify the quality of single-cell recording data. Also, they did not compare the velocity and position data before entering the open arm between the proximal and distal exploration. Thus, it is hard to reject the alternative hypothesis that the difference in neural activities between the exploration types stems from behavioral differences, not necessarily based on prediction signals.

Significance of the work:

This study should contribute to the single-cell-level understanding of the vHP with only very few experimental data available in the literature on the topic. Since some of the previous studies that claimed to record the ventral hippocampus actually targeted the intermediate portion of the hippocampus, this study would set a new standard for investigating the true ventral hippocampus.

Reviewer #1 (Public Review):

The authors set out to develop an in vitro model of multiple species representing diversity in the CF airway as a platform for a range of studies on why polymicrobial communities resist therapy. The rationale for their design is sound and the methods appear justifiable and reproducible. The major strength of this work is in producing a method for a range of future work, ideally for multiple groups in the field. The primary findings are interesting but not groundbreaking. One weakness in the method of reporting interspecies interactions and another in evaluating alternative causes of lasR advantages present opportunities for a stronger research contribution beyond this terrific method.

Reviewer #1 (Public Review):

The authors sought to explore the brain age paradigm in the early stages of Alzheimer's disease, focusing on the combination of different MRI modalities (brain structure derived from T1-weighted MRI and functional connectivity derived from resting state fMRI). Their goal was to understand how different unimodal brain ages and a combined multi-modal brain age related to risk factors related to Alzheimer's disease, namely hippocampal volume, cognitive performance, amyloid or tau positivity from PET scans or CSF data and neurofilaments. As part of this, they aimed to ascertain which brain age models performed more accurately.

The major strength of the methods is the novel combination of different MRI modalities using Gaussian Processes and stacking to predict brain age. Another strength is the use of multiple data sources for both model training and testing, reducing the reliance on a single site and decreasing the likelihood of overfitting, which should improve generalisability. A weakness is the poor fit of the functional connectivity model to the data, whereby the vast majority of test participants were shown to have younger appearing brains, even those with cognitive impairment. This indicates that an alternative fMRI processing pipeline could have been beneficial, however, no experimentation on this important facet of the analysis was included. Another weakness is the relatively limited sample size compared too much of the brain age literature and the failure to report the R^2 metric, which is important for the comparison of this study with previously published reports. Potentially, more accurate models would have led to clearer results, as there are a number of borderline findings which hinder clear interpretation.

In general, the study did meet its stated goals and was able to generate a multi-modality brain-age model and this model did show older appearing brains in people with cognitive impairment. This model also showed that people with older appearing brains had poorer cognition, lower hippocampal volume, and greater amyloid deposition. In people who met the criteria for being cognitively impaired, greater tau deposition on PET scans was associated with an older appearing brain. One claim of the study is that the multi-modality brain-age model was more accurate than the brain-volume model, however, it is unclear from the report whether appropriate statistics were used for this. The authors need to clarify exactly what procedure they undertook to compare the models, as they potentially employed an erroneous method (determining statistical significance based on the number of bootstraps instead of the number of observations) which may have led them to mistakenly claim better performance.

Given the relatively poor or equivocal performance of the brain age models and the relatively small sample sizes available, it is not clear that the modelling or dataset will have a big impact on the field. More accurate modelling methods are openly available, as are larger datasets. Nevertheless, the study is well-motivated and scientifically rigorous, so the results themselves are informative regarding the interrelationships of key Alzheimer's biomarkers and risk factors.

Reviewer #1 (Public Review):

This is an important, well-written and easily comprehended quantitative imaging study that analyzes the motion of endo-lysosomal compartments within axons in vivo using simultaneous multiphoton imaging in the mammalian brain. The simultaneous dual two-photon imaging is well-executed and represents a substantive advance in a field that relies heavily on in vitro neuronal culture preparations. This work opens the door to neurons that have aged appropriately and done so in the context of normal synaptic and neuromodulatory input, without an excess of added factors that occurs with in vitro cell culture. The authors solve an issue of cell polarity, providing strong support for their ability to determine directional movement (anterograde versus retrograde). In principle, this could become a generalized approach, opening this type of experiment up to other investigators. Finally, interesting differences in motion are observed, including activity-dependent and calcium-dependent changes that differ from measurements made in vitro. This is a significant technical advance with interesting observations that substantively move the field forward.

Reviewer #1 (Public Review):

The authors develop and freely disseminate the THINGS-data collection, a large-scale dataset incorporating MRI, MEG, eye-tracking, and 4.7 million similarity ratings for 1,854 object concepts. Demonstrating the reliability of their data, the authors replicate nearly a dozen previous neuroimaging papers. This "big data" approach significantly advances our ability to link behavioral measures with neuroimaging at scale, with the potential to spark future insights into how the mind represents objects.

I thought that the article was well-written, with a sound methodological approach, high-quality results, and well-supported conclusions. I am overall enthusiastic about this work, and I think THINGS will provide an important benchmark for future big data approaches in cognitive and computational neuroscience.

However, I thought it was also important to articulate more directly the potential insights this dataset can offer to the field. Although the authors mentioned that they "provided five examples for potential research directions", it was not clear to me what these new research directions were, given that the authors entirely describe replications in the results.

Reviewer #1 (Public Review):

The authors of this paper have used emerging environmental DNA (eDNA) methods to profile depth-biodiversity in two deep ocean trenches.

Strengths:<br /> Working in deep ocean habitats is challenging and this paper provides data from not one but two deep ocean trenches and provides new perspectives on biodiversity distributions in these habitats. The most interesting finding is that deep ocean habitats appear to contribute much more biodiversity, as measured using relatively novel eDNA techniques than previously thought. The comparison of these two trenches also illustrates the variable nature of these biodiversity patterns and suggests trench-specific features (e.g. unique habitats and/or productivity) influence deep ocean biodiversity.

Weaknesses:<br /> While eDNA methods are becoming more established, there remains skepticism by many in the scientific community about the origins of the detected DNA (e.g. does it drift in from other areas or water layers?). If these concerns aren't addressed (i.e. by citing supporting literature on the fate of eDNA), the different biodiversity profiles between trenches could possibly be explained by differing oceanography. There is also some important methodological information that is missing from this manuscript. For example, sampling volumes will affect the amount of biodiversity detected, but it is not clear if sample volumes are consistent across depths and study areas. It was also not indicated whether field controls (blanks) were taken to assess the potential contamination of samples. Lastly, the literature in the eDNA field is progressing rapidly and there are some missing papers (e.g Thomsen et al. 2016, Canals et al. 2021, McClenaghan et al. 2020, Govindarajan et al. 2021, etc.) that are relevant to the technique used in this manuscript and the habitat studied.

Reviewer #1 (Public Review):

This group previously demonstrated that trisomy 21 causes an increase in PCNT levels, and this increase leads to pericentrosomal crowding and inhibition of ciliogenesis in fibroblasts. The authors here use trisomy and tetrasomy 21 retinal pigment epithelium cells generated by microcell-mediated chromosome transfer (MMCT) and previously generated mouse models of human trisomy 21. The well-quantified data and well-reasoned paper compellingly demonstrate that modestly increased PCNT levels can attenuate ciliogenesis and may result in trisomy 21-associated phenotypes such as cerebellar growth defects.

Reviewer #1 (Public Review):

The authors describe a feeding system for killifish that allows high precision control of feeding amount and schedule on a per-tank basis. The system permits automation of this task using open-source and affordable components and software. Due to this emphasis, the system appears amenable to manufacture by individual research groups and the approach appears very scalable (although more detailed build, programming and assembly instructions and videos might be useful for groups with little experience with microcontrollers and manufacturing). An exciting aspect of the system is the possibility to modify the system for different purposes. For example, it might be possible to reduce the minimum feeding amount, thereby allowing more fine grained exploration of effects related to feeding shedule. I am very enthusiastic about the open-source "maker" aspects of this work.

The authors next explore two interesting applications of the system. First, they show that precise control of food allows automated investigation of lifespan extension under calorie restriction (CR) conditions. This is an important use case for a system of this type and showing that it is fit for this application is important.

Secondly, the authors show an exciting modification of the system that involves only addition of a simple red light LED. This modification allows use of the system in a associative learning / conditioning paradigm.

Finally, they show that there is an age-dependent decline in learning as evaluated by this conditioning paradigm. I am very enthusiastic about this additional function and, again, this example demonstrates the flexibly and open nature of the technology, suggesting that others can likely modify and expand the system to suite their own questions and applications. In summary, I am enthusiastic about the technology described and about the approach by which the system was developed.

However, at the current stage, the biological applications are essentially validation experiments - e.g. showing that CR can be implemented and that the system can be used for learning and memory experiments. Neither of these aspects is pursued beyond the basic validation experiments (showing that lifespan extension can be achieved and that there is age-dependent decline in associative learning).

Reviewer #1 (Public Review):

Leukemic cells are known to remodel bone marrow niche to promote their expansion and to suppress normal hematopoiesis. However, molecular mechanisms remain largely unknown. In this manuscript, authors developed new experimental models in mice to address this issue, using mouse BCR-ABL-driven ALL cells marked with YFP, or DOX-inducible MLL-AF9 AML cells. After transplantation of either of these cells, authors discovered suppression of host hematopoiesis. Using these systems, authors tested their hypothesis on lymphotoxin receptor-mediated interaction of the leukemic cells and stroma cells.

The main conclusions here are: 1) lymphotoxin signaling through its receptor mediates IL7 down regulation and alters gene expression related to inflammation etc. in stroma cells, 2) IL7 down regulation leads to reduction in B lymphoid cells but not myeloid cells, 3) lymphotoxin expression in leukemic cells is induced by DNA damage response, and 4) CXCR4, which is known to be induced in B cells in response to stroma cells, collaborates with DNA damage in induction of lymphotoxin in leukemic cells. Taken together, authors suggest that a positive feedback loop of leukemic cells and stroma cells for leukemic cell proliferation and normal hematopoietic suppression, involving lymphotoxin and CXCR4 in leukemic cells and lymphotoxin receptor in stroma cells. Generally, these conclusions and the model of the positive feedback regulation are supported, to a reasonable level, by the experimental results provided in the manuscript. However, some of the results show small effects of manipulations, leaving the pathological significance of the feedback model as a future issue.

Reviewer #1 (Public Review):

Polarization in cells and organs is often dictated by opposing polarity domains. In grass subsidiary cells, several proteins (including PAN1) were previously found to polarize in a discrete patch prior to asymmetric division. Zhang et al. identify POLAR via transcriptional profiling of Bdmute, a mutant that lacks subsidiary cells The authors effectively show that Bdpolar mutants have defective subsidiary cells. A distinctive and exciting localization pattern of POLAR is demonstrated, which is opposite to PAN1. This localization pattern is further contextualized by showing that PAN1 and MUTE are both required for POLAR's distinctive localization; however, PAN1 polarization is unaffected in both polar and mute. The integration of MUTE, POLAR, and PAN1 is particularly important as it integrates how polarity proteins and fate factors interact with each other.

Bdpolar mutants have defects in subsidiary cells that lead to defects in stomatal function. The authors carefully and quantitatively compare the phenotypes of pan1 and polar and conclude distinct roles for the two proteins based on differences in phenotypes including nuclear polarization, division site specification, and repeated rounds of cell division. The discovery and localization of POLAR are very exciting, but the comparison between single alleles of pan1 and polar and the extrapolation requires scrutiny. In particular, the data on division site specification in pan1 seem inconsistent with the % defective subsidiary cells and nuclear migration defects. However, these are addressable and given the exciting nature of the localization and pathway determination, the paper's impact stands.

Reviewer #1 (Public Review):

The research investigates the genetic basis for resistance to high CO2 levels in the human pathogenic fungus Cryptococcus neoformans. Screening collections of over 5,000 gene deletion strains revealed 96 with impaired growth, including a set of genes all related to the same RAM signaling pathway. Further genetic dissection was able compellingly to place where this pathway lies relative to upstream inputs and through the isolation of suppressor mutants as potential downstream targets of the pathway. Given the high levels of CO2 encountered by fungi in the human host, this work may provide new directions for the control of disseminated fungal disease.

The research presents both strengths and weaknesses.

Strengths include:

(1) One of the largest scale analyses of genes involved in growth under high CO2 concentrations in a fungus, revealing a set of just under 100 mutants with impaired growth.<br /> (2) Elegant genetic epistasis analysis to show where different components fit within a pathway of transmission of CO2 exposure. For example, over expression of one of the kinases, Cbk1, can overcome the CO2-sensitivity of mutations in the CDC24 or CNA1 genes (but not in the reciprocal overexpression direction).<br /> (3) Isolation of suppressor mutations in the cbk1 background, now able to grow at high CO2 levels, was able to lead to the identification of two genes. Follow up characterization, which included examining in vitro phenotypes, gene expression analysis, and impact during mouse infection was able to reveal that the two suppressors restore a subset of the phenotypes impacted by mutation of CBK1. Indeed, one conclusion from this careful work is that the reduced virulence of the cbk1 mutant is not due to its sensitivity to high levels of CO2, perhaps an unexpected finding given the original goals of the study towards linking CO2 sensitivity with decreased virulence.

Weaknesses include:

(1) What is the rationale for examining gene expression using the NanoString technology of 118 genes rather than a more genome-wide approach such as RNA-sequencing?<br /> (2) Without additional species examined, some of the conclusions about differences in impact between ascomycetes and basidiomycetes might instead reflect differences between species. For example, RAM mutants in other strains of C. neoformans do not exhibit so strong a temperature sensitive phenotype. Or to extend the comparison further, one might assume given the use of CO2 for Drosophila manipulations that the RAM pathway components in an insect would not be required for surviving high CO2.<br /> (3) Given the relative ease of generate progeny of this species, it would have been informative to explore if the suppressors of cbk1 also suppressed the loss of genes like CDC24, CNA1, etc, equivalent to the experiment performed of overexpression of CBK1 in those backgrounds.

Reviewer #1 (Public Review):

This paper makes an important contribution to the current debate on whether the diversity of a microbial community has a positive or negative effect on its own diversity at a later time point. In my view, the main contribution is linking the diversity-begets-diversity patterns, already observed by the same authors and others, to genomic signatures of gene loss that would be expected from the Black Queen Hypothesis, establishing an eco-evolutionary link. In addition, they test this hypothesis at a more fine-grained scale (strain-level variation and SNP) and do so in human microbiome data, which adds relevance from the biomedical standpoint. The paper is a well-written and rigorous analysis using state-of-the-art methods, and the results suggest multiple new experiments and testable hypotheses (see below), which is a very valuable contribution.

That being said, I do have some concerns that I believe should be addressed. First of all, I am wondering whether gene loss could also occur because of environmental selection that is independent of other organisms or the diversity of the community. An alternative hypothesis to the Black Queen is that there might have been a migration of new species from outside and then loss of genes could have occurred because of the nature of the abiotic environment in the new host, without relationship to the community diversity. Telling the difference between these two hypotheses is hard and would require extensive additional experiments, which I don't think is necessary. But I do think the authors should acknowledge and discuss this alternative possibility and adjust the wording of their claims accordingly.

Another issue is that gene loss is happening in some of the most abundant species in the gut. Under Black Queen though, we would expect these species to be most likely "donors" in cross-feeding interactions. Authors should also discuss the implications, limitations, and possible alternative hypotheses of this result, which I think also stimulates future work and experiments.

Regarding Figure 5B, there is a couple of questions I believe the authors should clarify. First, How is it possible that many species have close to 0 pathways? Second, besides the overall negative correlation, the data shows some very conspicuous regularities, e.g. many different "lines" of points with identical linear negative slope but different intercept. My guess is that this is due to some constraints in the pathway detection methods, but I struggle to understand it. I think the authors should discuss these patterns more in detail.

Finally, I also have some conceptual concerns regarding the genomic analysis. Namely, genes can be used for biosynthesis of e.g. building blocks, but also for consumption of nutrients. Under the Black Queen Hypothesis, we would expect the adaptive loss of biosynthetic genes, as those nutrients become provided by the community. However, for catabolic genes or pathways, I would expect the opposite pattern, i.e. the gain of catabolic genes that would allow taking advantage of a more rich environment resulting from a more diverse community (or at least, the absence of pathway loss). These two opposing forces for catabolic and biosynthetic genes/pathways might obscure the trends if all genes are pooled together for the analysis. I believe this can be easily checked with the data the authors already have, and could allow the authors to discuss more in detail the functional implications of the trends they see and possibly even make a stronger case for their claims.

Reviewer #1 (Public Review):

This is an interesting paper that presents a novel idea for the identification of risk factors amongst highly correlated traits in a Mendelian randomization paradigm - a previous investigation (https://www.ncbi.nlm.nih.gov/pmc/articles/PMC8438050/) has considered PCA, but not sparse PCA. There are clear conceptual reasons why sparse PCA may be an improvement, as detailed in this paper. Overall, the paper does a good job in terms of motivating this work and comparing the methods. A large chunk of the motivation for the method is conceptual (rather than empirical), and it's unlikely that any method would outperform others in all circumstances, but the authors do a good job of illustrating differences and giving a clear and qualified recommendation.

Reviewer #1 (Public Review):

There were two parts to this paper. The first was to build a network model with parameters carefully adjusted to match those seen in the turtle cortex. The second was to simulate the circuit, and show that it could produce reasonably repeatable patterns of activity in response to a single, externally added, spike.

As a model of the turtle cortex, the paper was pretty convincing. And the explanation for the repeatable patterns of activity - a small number of very strong connections and a very low background firing rate - seemed eminently reasonable. This paper should serve as a very good starting point for understanding computing in the turtle cortex.

However, average firing rates in the turtle are extremely low - 0.1 Hz, at least in these simulations. Their model is unlikely, therefore, to account for activity in the mammalian cortex, which exhibits a much higher background firing rate, and for which there's not a lot of evidence for the extremely strong connections seen in the turtle.