Reviewer #1 (Public review): <br /> The authors survey the ultrastructural organization of glutamatergic synapses by cryo-ET and image processing tools using two complementary experimental approaches. The first approach employs so-called "ultra-fresh" preparations of brain homogenates from a knock-in mouse expressing a GFP-tagged version of PSD-95, allowing Peukes and colleagues to specifically target excitatory glutamatergic synapses. In the second approach, direct in-tissue (using cortical and hippocampal regions) targeting of the glutamatergic synapses employing the same mouse model is presented. In order to ascertain whether the isolation procedure causes any significant changes in the ultrastructural organization (and possibly synaptic macromolecular organization) the authors compare their findings using both of these approaches. The quantitation of the synaptic cleft height reveals an unexpected variability, while the STA analysis of the ionotropic receptors provides insights into their distribution with respect to the synaptic cleft.

The main novelty of this study lies in the continuous claims by the authors that the sample preservation methods developed here are superior to any others previously used. This leads them as well to systematically downplay or directly ignore a substantial body of previous cryo-ET studies of synaptic structure. Without comparisons with the cryo-ET literature, it is very hard to judge the impact of this work in the field. Furthermore, the data does not show any better preservation in the so-called "ultra-fresh" preparation than in the literature, perhaps to the contrary as synapses with strangely elongated vesicles are often seen. Such synapses have been regularly discarded for further analysis in previous synaptosome studies (e.g. Martinez-Sanchez 2021). Whilst the targeting approach using a fluorescent PSD95 marker is novel and seems sufficiently precise, the authors use a somewhat outdated approach (cryo-sectioning) to generate in-tissue tomograms of poor quality. To what extent such tomograms can be interpreted in molecular terms is highly questionable. The authors also don't discuss the physiological influence of 20% dextran used for high-pressure freezing of these "very native" specimens.

Lastly, a large part of the paper is devoted to image analysis of the PSD which is not convincing (including a somewhat forced comparison with the fixed and heavy-metal staining room temperature approach). Despite being a technically challenging study, the results fall short of expectations.

Reviewer #2 (Public review):

Summary:

The authors set out to visualize the molecular architecture of the adult forebrain glutamatergic synapses in a near-native state. To this end, they use a rapid workflow to extract and plunge-freeze mouse synapses for cryo-electron tomography. In addition, the authors use knockin mice expression PSD95-GFP in order to perform correlated light and electron microscopy to clearly identify pre- and synaptic membranes. By thorough quantification of tomograms from plunge- and high-pressure frozen samples, the authors show that the previously reported 'post-synaptic density' does not occur at high frequency and therefore not a defining feature of a glutamatergic synapse.

Subsequently, the authors are able to reproduce the frequency of post-synaptic density when preparing conventional electron microscopy samples, thus indicating that density prevalence is an artifact of sample preparation. The authors go on to describe the arrangement of cytoskeletal components, membraneous compartments, and ionotropic receptor clusters across synapses.

Demonstrating that the frequency of the post-synaptic density in prior work is likely an artifact and not a defining feature of glutamatergic synapses is significant. The descriptions of distributions and morphologies of proteins and membranes in this work may serve as a basis for the future of investigation for readers interested in these features.

Strengths:

The authors perform a rigorous quantification of the molecular density profiles across synapses to determine the frequency of the post-synaptic density. They prepare samples using two cryogenic electron microscopy sample preparation methods, as well as one set of samples using conventional electron microscopy methods. The authors can reproduce previous reports of the frequency of the post-synaptic density by conventional sample preparation, but not by either of the cryogenic methods, thus strongly supporting their claim.

Reviewer #3 (Public review):

Summary:

The authors use cryo-electron tomography to thoroughly investigate the complexity of purified, excitatory synapses. They make several major interesting discoveries: polyhedral vesicles that have not been observed before in neurons; analysis of the intermembrane distance, and a link to potentiation, essentially updating distances reported from plastic-embedded specimen; and find that the postsynaptic density does not appear as a dense accumulation of proteins in all vitrified samples (less than half), a feature which served as a hallmark feature to identify excitatory plastic-embedded synapses.

Strengths:

(1) The presented work is thorough: the authors compare purified, endogenously labeled synapses to wild-type synapses to exclude artifacts that could arise through the homogenation step, and, in addition, analyse plastic embedded, stained synapses prepared using the same quick workflow, to ensure their findings have not been caused by way of purification of the synapses. Interestingly, the 'thick lines of PSD' are evident in most of their stained synapses.

(2) I commend the authors on the exceptional technical achievement of preparing frozen specimens from a mouse within two minutes.

(3) The approaches highlighted here can be used in other fields studying cell-cell junctions.

(4) The tomograms will be deposited upon publication which will enable neurobiologists and researchers from other fields to carry on data evaluation in their field of expertise since tomography is still a specialized skill and they collected and reconstructed over 100 excellent tomograms of synapses, which generates a wealth of information to be also used in future studies.

(5) The authors have identified ionotropic receptor positions and that they are linked to actin filaments, and appear to be associated with membrane and other cytosolic scaffolds, which is highly exciting.

(6) The authors achieved their aims to study neuronal excitatory synapses in great detail, were thorough in their experiments, and made multiple fascinating discoveries. They challenge dogmas that have been in place for decades and highlight the benefit of implementing and developing new methods to carefully understand the underlying molecular machines of synapses.

Impact on community:

The findings presented by Peukes et al. pertaining to synapse biology change dogmas about the fundamental understanding of synaptic ultrastructure. The work presented by the authors, particularly the associated change of intermembrane distance with potentiation and the distinct appearance of the PSD as an irregular amorphous 'cloud' will provide food for thought and an incentive for more analysis and additional studies, as will the discovery of large membranous and cytosolic protein complexes linked to ionotropic receptors within and outside of the synaptic cleft, which are ripe for investigation. The findings and tomograms available will carry far in the synapse fields and the approach and methods will move other fields outside of neurobiology forward. The method and impactful results of preparing cryogenic, unlabeled, unstained, near-native synapses may enable the study of how synapses function at high resolution in the future.

Reviewer #1 (Public review):

Summary:

This paper is a relevant overview of the currently published literature on low-intensity focused ultrasound stimulation (TUS) in humans, with a meta-analysis of this literature that explores which stimulation parameters might predict the directionality of the physiological stimulation effects.

The pool of papers to draw from is small, which is not surprising given the nascent technology. It seems nevertheless relevant to summarize the current field in the way done here, not least to mitigate and prevent some of the mistakes that other non-invasive brain stimulation techniques have suffered from, most notably the theory- and data-free permutation of the parameter space.

The meta-analysis concludes that there are, at best, weak trends toward specific parameters predicting the direction of the stimulation effects. The data have been incorporated into an open database that will ideally continue to be populated by the community and thereby become a helpful resource as the field moves forward.

Strengths:

The current state of human TUS is concisely and well summarized. The methods of the meta-analysis are appropriate. The database is a valuable resource.

Suggestions:<br /> - The paper remains lengthy and somewhat unfocused, to the detriment of readability. One can understand that the authors wish to include as much information as possible, but this reviewer is sceptical that this will aid the use of the databank, or help broaden the readership. For one, there is a good chunk of repetition throughout. The intro is also somewhat oscillating between TMS, tDCS and TUS. While the former two help contextualizing the issue, it doesn't seem necessary. In the section on clinical applications of TUs and possible outcomes of TUS, there's an imbalance of the content across examples. That's in part because of the difference in knowledge base but some sections could probably be shortened, eg stroke. In any case, the authors may want to consider whether it is worth making some additional effort in pruning the paper

- The terms or concept of enhancement and suppression warrant a clearer definition and usage. In most cases, the authors refer to E/S of neural activity. Perhaps using terms such as "neural enhancement" etc helps distinguish these from eg behavioural or clinical effects. Crucially, how one maps onto the other is not clear. But in any case, a clear statement that the changes outlined on lines 277ff do not

- Re tb-TUS (lines 382ff), it is worth acknowledging here that independent replication is very limited (eg Bao et al 2024; Fong et al bioRxiv 2024) and seems to indicate rather different effects

- The comparison with TPS is troublesome. For one, that original study was incredibly poorly controlled and designed. Cherry-picking individual (badly conducted) proof-of-principle studies doesn't seem a great way to go about as one can find a match for any desired use or outcome.

Moreover, other than the concept of "pulsed" stimulation, it is not clear why that original study would motivate the use of TUS in the way the authors propose; both types of stimulation act in very different ways (if TPS "acts" at all). But surely the cited TPS study does not "demonstrate the capability for TUS for pre-operative cognitive mapping". As an aside, why the authors feel the need to state the "potential for TPS... to enhance cognitive function" is unclear, but it is certainly a non-sequitur. This review feels quite strongly that simplistic analogies such as the one here are unnecessary and misleading, and don't reflect the thoughtful discussion of the rest of the paper. In the other clinical examples, the authors build their suggestions on other TUS studies, which seems more sensible.

Reviewer #1 (Public review):

Summary:

This manuscript presents a compelling and innovative approach that combines Track2p neuronal tracking with advanced analytical methods to investigate early postnatal brain development. The work provides a powerful framework for exploring complex developmental processes such as the emergence of sensory representations, cognitive functions, and activity-dependent circuit formation. By enabling the tracking of the same neurons over extended developmental periods, this methodology sets the stage for mechanistic insights that were previously inaccessible.

Strengths:

(1) Innovative Methodology:<br /> The integration of Track2p with longitudinal calcium imaging offers a unique capability to follow individual neurons across critical developmental windows.

(2) High Conceptual Impact:<br /> The manuscript outlines a clear path for using this approach to study foundational developmental questions, such as how early neuronal activity shapes later functional properties and network assembly.

(3) Future Experimental Potential:<br /> The authors convincingly argue for the feasibility of extending this tracking into adulthood and combining it with targeted manipulations, which could significantly advance our understanding of causality in developmental processes.

(4) Broad Applicability:<br /> The proposed framework can be adapted to a wide range of experimental designs and questions, making it a valuable resource for the field.

Weaknesses:

No major weaknesses were identified by this reviewer. The manuscript is conceptually strong and methodologically sound. Future studies will need to address potential technical limitations of long-term tracking, but this does not detract from the current work's significance and clarity of vision.

Reviewer #2 (Public review):

Summary:

The manuscript by Majnik and colleagues introduces "Track2p", a new tool designed to track neurons across imaging sessions of two-photon calcium imaging in developing mice. The method addresses the challenge of tracking cells in the growing brain of developing mice. The authors showed that "Track2p" successfully tracks hundreds of neurons in the barrel cortex across multiple days during the second postnatal week. This enabled the identification of the emergence of behavioral state modulation and desynchronization of spontaneous network activity around postnatal day 11.

Strengths:

The manuscript is well written, and the analysis pipeline is clearly described. Moreover, the dataset used for validation is of high quality, considering the technical challenges associated with longitudinal two-photon recordings in mouse pups. The authors provide a convincing comparison of both manual annotation and "CellReg" to demonstrate the tracking performance of "Track2p". Applying this tracking algorithm, Majnik and colleagues characterized hallmark developmental changes in spontaneous network activity, highlighting the impact of longitudinal imaging approaches in developmental neuroscience. Additionally, the code is available on GitHub, along with helpful documentation, which will facilitate accessibility and usability by other researchers.

Weaknesses:

(1) The main critique of the "Track2p" package is that, in its current implementation, it is dependent on the outputs of "Suite2p". This limits adoption by researchers who use alternative pipelines or custom code. One potential solution would be to generalize the accepted inputs beyond the fixed format of "Suite2p", for instance, by accepting NumPy arrays (e.g., ROIs, deltaF/F traces, images, etc.) from files generated by other software. Otherwise, the tool may remain more of a useful add-on to "Suite2p" (see https://github.com/MouseLand/suite2p/issues/933) rather than a fully standalone tool.

(2) Further benchmarking would strengthen the validation of "Track2p", particularly against "CaIMaN" (Giovannucci et al., eLife, 2019), which is widely used in the field and implements a distinct registration approach.

(3) The authors might also consider evaluating performance using non-consecutive recordings (e.g., alternate days or only three time points across the week) to demonstrate utility in other experimental designs.

Reviewer #3 (Public review):

Summary:

In this manuscript, Majnik et al. developed a computational algorithm to track individual developing interneurons in the rodent cortex at postnatal stages. Considerable development in cortical networks takes place during the first postnatal weeks; however, tools to study them longitudinally at a single-cell level are scarce. This paper provides a valuable approach to study both single-cell dynamics across days and state-driven network changes. The authors used Gad67Cre mice together with virally introduced TdTom to track interneurons based on their anatomical location in the FOV and AAVSynGCaMP8m to follow their activity across the second postnatal week, a period during which the cortex is known to undergo marked decorrelation in spontaneous activity. Using Track2P, the authors show the feasibility of tracking populations of neurons in the same mice, capturing with their analysis previously described developmental decorrelation and uncovering stable representations of neuronal activity, coincident with the onset of spontaneous active movement. The quality of the imaging data is compelling, and the computational analysis is thorough, providing a widely applicable tool for the analysis of emerging neuronal activity in the cortex. Below are some points for the authors to consider.

Major points:

(1) The authors used 20 neurons to generate a ground truth dataset. The rationale for this sample size is unclear. Figure 1 indicates the capability to track ~728 neurons. A larger ground truth data set will increase the robustness of the conclusions.

(2) It is unclear how movement was scored in the analysis shown in Figure 5A. Was the time that the mouse spent moving scored after visual inspection of the videos? Were whisker and muscle twitches scored as movement, or was movement quantified as the amount of time during which the treadmill was displaced?

(3) The rationale for binning the data analysis in early P11 is unclear. As the authors acknowledged, it is likely that the decoder captured active states from P11 onwards. Because active whisking begins around P14, it is unlikely to drive this change in network dynamics at P11. Does pupil dilation in the pups change during locomotor and resting states? Does the arousal state of the pups abruptly change at P11?

Reviewer #1 (Public review):

Summary:

In this manuscript, Mack and colleagues investigate the role of posttranslational modifications, including lysine acetylation and ubiquitination, in methyltransferase activity of SETD2 and show that this enzyme functions as a tumor suppressor in a KRASG12C-driven lung adenocarcinoma. In contrast to H3K36me2-specific oncogenic methyltransferases, the deletion of SETD2, which is capable of H3K36 trimethylation, increases lethality in a KRASG12C-driven lung adenocarcinoma mouse tumor model. In vitro, the authors demonstrate that polyacetylation of histone H3, particularly of H3K27, H3K14, and H3K23, promotes the catalytic activity of SETD2, whereas ubiquitination of H2A and H2B has no effect.

Strengths:

Overall, this is a well-designed study that addresses an important biological question regarding the functioning of the essential chromatin component. The manuscript contains excellent quality data, and the conclusions are convincing and justified. This work will be of interest to many biochemists working in the field of chromatin biology and epigenetics.

Weaknesses:

A minor comment: labels should be added in the Figures and should be uniform across all Figures (some are distorted).

Reviewer #2 (Public review):

Summary:

Human histone H3K36 methyltransferase Setd2 has been previously shown to be a tumor suppressor in lung and pancreatic cancer. In this manuscript by Mack et al., the authors first use a mouse KRASG12D-driven lung cancer model to confirm in vivo that Setd2 depletion exacerbates tumorigenesis. They then investigate the enzymatic regulation of the Setd2 SET domain in vitro, demonstrating that H2A, H3, or H4 acetylation stimulates Setd2-SET activity, with specific enhancement by mono-acetylation at H3K14ac or H3K27ac. In contrast, histone ubiquitination has no effect. The authors propose that H3K27ac may regulate Setd2-SET activity by facilitating its binding to nucleosomes. This work provides insight into how cross-talk between histone modifications regulates Setd2 function. However, the manuscript lacks a clear discussion on how Setd2's in vivo tumor suppressor role and the in vitro mechanistic regulation findings are connected. Additionally, some experiments require more controls and better data quality for proper interpretation.

Specific comments:

(1) As for Figure 2F, Setd2-SET activity on WT rNuc (H3) appears to be significantly lower compared to what is extensively reported in the literature. This is particularly puzzling given that Figure 2B suggests that using 3H-SAM, H3-nuc are much better substrates than K36me1, whereas in Figure 3F, rH3 is weaker than K36me1. It is recommended for the authors to perform additional experimental repeats and include a quantitative analysis to ensure the consistency and reliability of these findings.

(2) The additional bands observed in Figure 4B, which appear to be H4, should be accompanied by quantification of the intensity of the H3 bands to better assess K36me3 activity. Additionally, the quantification presented in Figure 4C for SAH does not seem accurate as it potentially includes non-specific methylation activity, likely from H4. This needs to be addressed for clarity and accuracy.

(3) In Figure 4E, the differences between bound and unbound substrates are not sufficiently pronounced. Given the modest differences observed, authors might want to consider repeating the assay with sufficient replicates to ensure the results are statistically robust.

(4) Regarding labeling, there are multiple issues that need correction: In the depiction of Epicypher's dNuc, it is crucial to clearly mark H2B as the upper band, rather than ambiguously labeling H2A/H2B together when two distinct bands are evident. In Figure 3B and D, the histones appear to be mislabeled, and the band corresponding to H4 has been cut off. It would be beneficial to refer to Figure 3E for correct labeling to maintain consistency and accuracy across figures.

(5) There are issues with the image quality in some blots; for instance, Figure 2EF and Figure 2D exhibit excessive contrast and pixelation, respectively. These issues could potentially obscure or misrepresent the data, and thus, adjustments in image processing are recommended to provide clearer, more accurate representations.

(6) The authors are recommended to provide detailed descriptions of the materials used, including catalog numbers and specific products, to allow for reproducibility and verification of experimental conditions.

(7) The identification of Setd2 as a tumor suppressor in KrasG12C-driven LUAD is a significant finding. However, the discussion on how this discovery could inspire future therapeutic approaches needs to be more balanced. The current discussion (Page 10) around the potential use of inhibitors is somewhat confusing and could benefit from a clearer explanation of how Setd2's role could be targeted therapeutically. It would be beneficial for the authors to explore both current and potential future strategies in a more structured manner, perhaps by delineating between direct inhibitors, pathway modulators, and other therapeutic modalities.

Reviewer #1 (Public review):

The small conductance calcium-activated potassium channel 2 (SK2) is an important drug target for treating neurological and cardiovascular diseases. However, structural information on this subtype of SK channels has been lacking, and it has been difficult to draw conclusions about activator and inhibitor binding and action in the absence of structural information.

Here the authors set out to (1) determine the structure of the transmembrane regions of a mammalian SK2 channel, (2) determine the binding site of apamin, a historically important SK2 inhibitor whose mode of action is unclear, and (3) use the structural information to generate a novel set of activators/inhibitors that selectively target SK2.

The authors largely achieved all the proposed goals, and they present their data clearly.

Unable to solve the structure of the human SK2 due to excessive heterogeneity in its cytoplasmic regions, the authors create a chimeric construct using SK4, whose structure was previously solved, and use it for structural studies. The data reveal a unique extracellular structure formed by the S2-S3 loop, which appears to directly interact with the selectivity filter and modulate its conductivity. Structures of SK2 in the absence and presence of the activating Ca2+ ions both possess non-K+-selective/conductive selectivity filters, where only sites 3 and 4 are preserved. The S6 gates are captured in closed and open states, respectively. Apamine binds to the S2-S3 loop, and unexpectedly, induces a K+ selective/conductive conformation of the selectivity filter while closing the S6 gate.

Through high-throughput screening of small compound libraries and compound optimization, the group identified a reasonably selective inhibitor and a related compound that acts as an activator. The characterization shows that these compounds bind in a novel binding site. Interestingly, the inhibitor, despite binding in a site different from that of apamine, also induces a K+ selective/conductive conformation of the selectivity filter while the activator induces a non-K+ selective/conductive conformation and an open S6 gate.

The data suggest that the selectivity filter and the S6 gate are rarely open at the same time, and the authors hypothesize that this might be the underlying reason for the small conductance of SK2. The data will be valuable for understanding the mechanism of SK2 channel (and other SK subtypes).

Overall, the data is of good quality and supports the claims made by the authors. However, a deeper analysis of the cryo-EM data sets might yield some important insights, i.e., about the relationship between the conformation of the selectivity filter and the opening of the S6 gate.

Some insight and discussion about the allosteric networks between the SF and the S6 gate would also be a valuable addition.

Reviewer #2 (Public review):

Summary:

The authors have used single-particle cryoEM imaging to determine how small-molecule regulators of the SK channel interact with it and modulate their function.

Strengths:

The reconstructions are of high quality, and the structural details are well described.

Weaknesses:

The electrophysiological data are poorly described. Several details of the structural observations require a mechanistic context, perhaps better relating them to what is known about SK channels or other K channel gating dynamics.

The most pressing point I have to make, which could help improve the manuscript, relates to the selectivity filter (SF) conformation. Whether the two ion-bound state of SK2-4 (Figure 4A) represents a non-selective, conductive SF occluded by F243 or represents a C-type inactivated SF, further occluded by F243, is unclear. It would be important to discuss this. Reconstructions of Kv1.3 channels also feature a similar configuration, which has been correlated to its accelerated C-type inactivation.

Furthermore, binding of a toxin derivative to Kv1.3 restores the SF into a conductive form, though occluded by the toxin. It appears that apamin binding to SK2-4 might be doing something similar. Although I am not sure whether SK channels undergo C-type inactivation like gating, classical MTS accessibility studies have suggested that dynamics of the SF might play a role in the gating of SK channels. It would be really useful (if not essential) to discuss the SF dynamics observed in the study and relate them better to aspects of gating reported in the literature.

The SF of K channels, in conductive states, are usually stabilized by an H-bond network involving water molecules bridged to residues behind the SF (D363 in the down-flipped conformation and Y361). Considering the high quality of the reconstructions, I would suspect that the authors might observe speckles of density (possibly in their sharpened map) at these sites, which overlap with water molecules identified in high-resolution X-ray structures of KcsA, MthK, NaK, NaK2K, etc. It could be useful to inspect this region of the density map.

Reviewer #3 (Public review):

This is a fundamentally important study presenting cryo-EM structures of a human small conductance calcium-activated potassium (SK2) channel in the absence and presence of calcium, or with interesting pharmacological probes bound, including the bee toxin apamin, a small molecule inhibitor, and a small molecule activator. As efforts to solve structures of the wild-type hSK2 channel were unsuccessful, the authors engineered a chimera containing the intracellular domain of the SK4 channel, the subtype of SK channel that was successfully solved in a previous study (reference 13). The authors present many new and exciting findings, including opening of an internal gate (similar to SK4), for the first time resolving the S3-S4 linker sitting atop the outer vestibule of the pore and unanticipated plasticity of the ion selectivity filter, and the binding sites for apamin, one new small molecule inhibitor and another small molecule activator. Appropriate functional data are provided to frame interpretations arising from the structures of the chimeric protein; the data are compelling, the interpretations are sound, and the writing is clear. This high-quality study will be of interest to membrane protein structural biologists, ion channel biophysicists, and chemical biologists, and will be valuable for future drug development targeting SK channels.

The following are suggestions for strengthening an already very strong and solid manuscript:

(1) It would be good to include some information in the text of the results section about the method and configuration used to obtain electrophysiological data and the limitations. It is not until later in the text that the Qube instrument is mentioned in the results section, and it is not until the methods section that the reader learns it was used to obtain all the electrophysiological data. Even there, it is not explicitly mentioned that a series of different internal solutions were used in each cell where the free calcium concentration was varied to obtain the data in Figure1C. Also, please state the concentration of free calcium for the data in Figure 1B.

(2) The authors do a nice job of discussing the conformations of the selectivity filter they observed here in SK as they relate to previous work on NaK and HCN, but from my perspective the authors are missing an opportunity to point out even more striking relationships with slow C-type inactivation of the selectivity filter in Shaker and Kv1 channels. C-type inactivation of the filter in Shaker was seen in 150 mM K using the W434F mutant (PMC8932672) or in 4 mM K for the WT channel (PMC8932672), and similar results have been reported for Kv1.2 (PMC9032944; PMC11825129) and for Kv1.3 (PMC9253088; PMC8812516) channels. For Kv1.3, C-type inactivation occurs even in 150 mM K (PMC9253088; PMC8812516). Not unlike what is seen here with apamin, binding of the sea anemone toxin (ShK) with a Fab attached (or the related dalazatide) inserts a Lys into the selectivity filter and stabilizes the conducting conformation of Kv1.3 even though the Lys depletes occupancy of S1 by potassium (PMC9253088; PMC8812516). What is known about how the functional properties of SK2 channels (where the filter changes conformation) differ from SK4, where the filter remains conducting (reference 13)? Is there any evidence that SK2 channels inactivate? Or might the conformation of the filter be controlled by regulatory processes in SK2 channels? I think connecting the dots here would enhance the impact of this study, even if it remains relatively speculative.

Reviewer #1 (Public Review):

Summary:

In this paper, Li and colleagues overcome solubility problems to determine the structure of FtsEX bound to EnvC from E. coli.

Strengths:

The structural work is well done, and the work is consistent with previous work on the structure of this complex from P. aerugionsa.

Weaknesses:

The model does not take into account all the information that the authors obtained, as well as known in vivo data.

The work lacks a clear comparison to the Pseudomonas structure highlighting new information that was obtained so that it is readily available to the reader.

The authors set out to obtain the structure of FtsEX-EnvC complex from E. coli. Previously, they were unable to do so but were able to determine the structure of the complex from P. aeruginosa. Here they persisted in attacking the E. coli complex since more is known about its involvement in cell division and there is a wealth of mutants in E. coli. The structural work is well done and recapitulates the results this lab obtained with this complex from P. aeruginosa. It would be helpful to compare more directly the results obtained here with the E. coli complex with the previously reported P. aeruginosa complex - are they largely the same or has some insight been obtained from the work that was not present in the previous complex from P. aeruginosa. This is particularly the case in discussing the symmetrical FtsX dimer binding to the asymmetrical EnvC, since this is emphasized in the paper. However, Figures 3C & D of this paper appear similar to Figures 2D & E of the P. aeruginosa structure. Presumably, the additional information obtained and presented in Figure 4 is due to the higher resolution, but this needs to be highlighted and discussed to make it clear to a general audience.

The main issue is the model (Figure 6). In the model ATP is shown to bind to FtsEX before EnvC, however, in Figure 1c, it is shown that ADP is sufficient to promote binding of FtsEX to EnvC.

The work here is all done in vitro, however, information from in vivo needs to be considered. In vivo results reveal that the ATP-binding mutant FtsE(D162N)X promotes the recruitment of EnvC (Proc Natl Acad Sci U S A 2011 108:E1052-60). Thus, even FtsEX in vivo can bind EnvC without ATP (not sure if this mutant can bind ADP).

Perhaps the FtsE protein from E. coli has to have bound nucleotides to maintain its 3D structure.

Comments after revision:

The most interesting aspect of this complex is that it has yet to be determined the order of events in the ATPase cycle as the authors acknowledge. Although the authors have responded quite well to the comments, I am still worried about the significance of the in vitro results compared to the in vivo results reported by others. In vivo ATP binding does not appear required for complex formation (of course it is possible that ADP is responsible in vivo). Have the authors tried to solve the complex with ADP since they suggested that it is sufficient to hold the complex together). If possible, it would confirm the role of ATP binding by comparing the structures. Also, it is not clear if ADP binds to any of the mutants made by the Bernhardt lab (D162N, K41M). If they do not bind ADP then FtsEX without nucleotide is able to bind EnvC as the authors indicate is the case in Pseudomonas. It is also unclear the significance of the ATPase activity of FtsEX in vitro with or without EnvC. Could the activity be some basal activity that is not relevant to the in vivo situation. If EnvC caused FtsEX to hydrolyze ATP it would be a futile cycle as FtsEX and EnvC are localized to the septum long before they are involved septal hydrolysis.

Reviewer #2 (Public Review):

Summary:

Peptidoglycan remodeling, particularly that carried out by enzymes known as amidases, is essential for the later stages of cell division including cell separation. In E. coli, amidases are generally activated by the periplasmic proteins EnvC (AmiA and AmiB) and NlpD (AmiC). The ABC family member, FtsEX, in turn, has been implicated as a modulator of amidase activity through interactions with EnvC. Specifically, how FtsEX regulates EnvC activity in the context of cell division remains unclear.

Strengths:

Li et al. make two primary contributions to the study of FtsEX. The first, the finding that ATP binding stabilizes FtsEX in vitro, enables the second, structural resolution of full-length FtsEX both alone (Figure 2) and in combination with EnvC (Figure 3). Leveraging these findings, the authors demonstrate that EnvC binding stimulates FtsEX-mediated ATP hydrolysis approximately two-fold. The authors present structural data suggesting EnvC binding leads to a conformational change in the complex. Biochemical reconstitution experiments (Figure 5) provide compelling support for this idea.

Weaknesses:

The potential impact of the study is curtailed by the lack of experiments testing the biochemical or physiological relevance of the model which is derived almost entirely from structural data.

Altogether the data support a model in which interaction with EnvC, results in a conformational change stimulating ATP hydrolysis by FtsEX and EnvC-mediated activation of the amidases, AmiA and AmiB. However, the study is limited in both approach and scope. The importance of interactions revealed in the structures to the function of FtsEX and its role in EnvC activation are not tested. Adding biochemical and/or in vivo experiments to fill in this gap would allow the authors to test the veracity of the model and increase the appeal of the study beyond the small number of researchers specifically interested in FtsEX.

Comments after revision:

Although I appreciate the authors' desire to save future biochemical experiments for a separate publication, the lack of in vitro data verifying their model makes it challenging to reconcile with published studies from other groups. The other reviewer's point about EnvC activating FtsEX ATPase activity resulting in a futile cycle since both are recruited to the septum well before constriction, is a good example of the disconnect between the model presented here and in vivo data.

Reviewer #1 (Public review):

Summary:

Melin et al. developed a quantitative assay to measure the fate of paternal mitochondria after fertilization. They combine this assay with C. elegans genetics to show that multiple genes contribute to paternal mitochondrial elimination. However, despite their claims, they unconvincingly place these genes into distinct pathways and fail to determine whether additional unknown genes are involved in the process.

Strengths:

Melin et al. develop a new assay to quantify the fate of paternal mitochondria during embryonic development in C. elegans. They use complex C. elegans genetics to disrupt 5 different genes and nicely measure their contributions to paternal mitochondrial elimination. In an attempt to place these genes into pathways, the authors interrupt genes in various combinations and measure paternal mitochondrial persistence. The authors discovered that disrupting 4 of the genes known to contribute to paternal mitochondrial elimination still resulted in paternal mitochondrial elimination, suggesting that more genes also contribute to this process. Finally, the authors discovered that pink-1, which had previously been discounted, indeed contributes to paternal mitochondrial elimination when the major pathway involving allo-1 is also disrupted.

Weaknesses:

In the introduction, the authors describe the importance of studying the maternal inheritance of mitochondrial DNA. However, the authors mostly study the inheritance of paternally-derived mitochondrial proteins (HSP6::GFP). While the authors do use a PCR approach to measure paternal mitochondrial DNA, their results are not as quantitative and thorough (applied to multiple mutant combinations) as their microscopy assay. Using their microscopy assay, the authors did not combine mutants for all 5 genes. Therefore, they cannot support or discount the possibility that undiscovered paternal mitochondrial elimination mechanisms exist. The author's genetic epistasis experiments are incomplete and occasionally improperly interpreted (as described below). Finally, the authors were unable to achieve paternal mitochondrial transmission to the F2 generation (which admittedly has not been achieved in any experimental system).

Reviewer #2 (Public review):

Summary:

Mitochondrial DNA (mtDNA) is exclusively maternally transmitted in almost all species. Paternal mitochondria, with their mtDNA, must be rapidly degraded after fertilisation to prevent their transmission to progeny, which could lead to subsequent detrimental mito-nuclear incompatibilities. Multiple layers of mechanisms contribute to blocking the transmission of paternal mitochondria and their mtDNA to progeny. Endonuclease activity and mitophagy form a part of these strategies. However, other key regulatory mechanisms remain to be elucidated, as inactivating endonuclease and mitophagy activity only delays the clearance of paternal mitochondria. In this study, the authors mainly focused on genes involved in endonuclease function (csp-6) and autophagy (allo-1) in C. elegans, demonstrating a synergic genetic interaction that potentialize their activity. They also revealed a contribution by pink-1/pink1, in the absence of allo-1.

Strengths:

The majority of data relies on confocal microscopy images and corresponding image analysis and quantification. Images are clear, and quantifications are supported by several biological replicates of >10 n and standard statistical tests. Mutants used were obtained from the Caenorhabditis Genetics Center (CGC) and were previously validated and confirmed by the C. elegans community. The scientific approach is solid and rigorous and in line with state-of-the-art C. elegans methods. Proper controls have been performed to rule out the effect of animal viability on observed results or to confirm the staining validity of TUBES on subcellular structures surrounding paternal mitochondria. Controls validating uaDf5 PCR specificity were conducted.

Weaknesses:

However, the embryonic expression of paternally contributing genes in feminised animals cannot be completely ruled out, as RNAi was used to alter gene expression levels. An issue inherent to RNAi approaches. Also, the impact of pink-1/pink1 is significant, but there is a lack of evidence demonstrating its mitophagic function.

Goal achievements and data supportive of conclusions:

In the first part of the study, the authors strongly and clearly demonstrate the synergistic interaction between the csp-6 and allo-1 in delaying paternal mitochondria degradation and associated mtDNA in the fertilised egg. In wild-type animals, paternal mitochondria are visible (using a mitochondrial HSP-::GFP marker) until the 4-cell stage embryo. In the csp-6; allo-1 double mutant genetic background, paternal mitochondria very significantly perdures until the 2-fold embryonic stage. The uaDf5 mitochondrial deletion, detectable by PCR, that was introduced by crossing with a male, followed the same trend. In addition, loss of fncd-1/fndc1 and phb-2 did not extend the perdurance of paternal mitochondria. In the second part of the study, the authors demonstrate a contribution of the loss of pink-1/pink1, in the absence of allo-1, in delaying paternal mitochondria degradation until the 100-cell stage. Overall, the conclusions are in accordance with the data shown.

Impact on the field:

Endonuclease activity and mitophagy aren't sufficient to prevent the transmission of paternal mitochondria and associated mtDNA to progeny, but they still contribute significantly to regulating the perdurance of paternal mitochondria in early embryos. Understanding how these two functions work in concert to potentialize their activity is important, as they could potentially be manipulated/enhanced to improve paternal mitochondrial degradation in the future. Here, the authors demonstrate a detailed synergistic genetic interaction between these functions. Also, they pointed out a new potential contribution of pink-1/pink1, which may underlie a potentially more complex mitophagic protective function.

Reviewer #3 (Public review):

Summary:

The present study examines the cooperation among four allophagy/mitophagy factors, ALLO-1, CPS-6, FNDC-1, and PHB-2, implicated in the elimination of the sperm-derived mitochondria in C. elegans embryos. The key finding of the cumulative effect of ALLO-1 and CPS-6 inactivation causing delayed sperm mitophagy is significant for the understanding of mitochondrial inheritance in the nematode model and in general. Below are some specific suggestions on how the impact of the article could be elevated:

Abstract:

The authors should shorten the description of previously identified mitophagy factors and provide more detail on the present study results. An impact statement should be added at the end, with significance for understanding mitochondrial inheritance across taxa, all the way to mammals/humans.

Introduction:

The authors should provide more details on ALLO-1 and its interaction with LC-3. Also, it should be specified which of those previously identified allophagy factors are unique to worms and which ones are conserved. See also my comment below about including a diagram and a table of pathways and determinants involved in allophagy/paternal mitophagy.

Results:

If I understand the mtDNA data correctly, paternal mtDNA is maintained throughout the lifespan of the F1 generation but absent from the F2 generation. This is reminiscent of past studies of interspecific Mus musculus/Mus spretus mouse crosses by Kaneda/Shitara in which the paternal mtDNA was maintained F1 generation, resulting in heteroplasmy, but was lost from the F2 generation after back-crossing. Are CPS-6 and ALLO-1 effectors, but not determinants of maternal mtDNA inheritance in the nematode?

The finding that PINK-1 inactivation stabilizes sperm-derived mitochondria in the embryos is interesting. Are the substrates of PINK1 known in C. elegans? This could provide a clue concerning the aforementioned mitophagy determinants acting independently of ALLO-1.

Discussion:

A summary-diagram compiling the intersecting allophagy pathways would be helpful to accompany discussion, in addition to or expanding on the simple diagram presented as Figure 5; also, a table listing all the factors implicated in nematode allophagy next to those implicated in human/mammalian sperm mitophagy, which would highlight the divergences and overlaps between vertebrates and invertebrates.

Is it known how CPC-6 enters/gets imported into the sperm mitochondria inside the embryo? This pathway could potentially be targeted to decipher the allophagy mechanism.

PINK/PARKIN/PACRG and FUNDC1/2 pathways have been implicated in mammalian neurodegeneration as well as in mitophagy, including but not limited to sperm mitophagy after fertilization. These pathways in mammals should be briefly reviewed as they may provide further clues to how the allophagy pathways intersect in C. elegans.

Reviewer #1 (Public review):

Summary:

This study provides new insight into the non-canonicial voltage-gating mechanism of BK channels through prolonged (10 us) MD simulations of the Slo1 transmembrane domain conformation and K+ conduction in response to high imposed voltages (300, 750 mV). The results support previous conclusions based on functional and structural data and MD simulations that the voltage-sensor domain (VSD) of Slo1 undergoes limited conformational changes compared to Kv channels, and predicts gating charge movement comparable in magnitude to experimental results. The gating charge calculations further indicate that R213 and R210 in S4 are the main contributors owing to their large side chain movements and the presence of a locally focused electric field, consistent with recent experimental and MD simulation results by Carrasquel-Ursulaez et al.,2022. Most interestingly, changes in pore conformation and K+ conduction driven by VSD activation are resolved, providing information regarding changes in VSD/pore interaction through S4/S5/S6 segments proposed to underly electromechanical coupling.

Strengths:

Include that the prolonged timescale and high voltage of the simulation allow apparent equilibration in the voltage-sensor domain (VSD) conformational changes and at least partial opening of the pore. The study extends the results of previous MD simulations of VSD activation by providing quantitative estimates of gating charge movement, showing how the electric field distribution across the VSD is altered in resting and activated states, and testing the hypothesis that R213 and R210 are the primary gating charges by steered MD simulations. The ability to estimate gating charge contributions of individual residues in the WT channel is useful as a comparison to experimental studies based on mutagenesis which have yielded conflicting results that could reflect perturbations in structure. Use of dynamic community analysis to identify coupling pathways and information flow for VSD-pore (electromechanical) coupling as well as analysis of state-dependent S4/S5/S6 interactions that could mediate coupling provide useful predictions extending beyond what has been experimentally tested.

Weaknesses:

Weaknesses include that a truncated channel (lacking the C-terminal gating ring) was used for simulations, which is known to have reduced single channel conductance and electromechanical coupling compared to the full-length channel. In addition, as VSD activation in BK channels is much faster than opening, the timescale of simulations was likely insufficient to achieve a fully open state as supported by differences in the degree of pore expansion in replicate simulations, which are also smaller than observed in Ca-bound open structures of the full-length channel. Taken together, these limitations suggest that inferences regarding coupling pathways and interactions in the fully open voltage-activated channel may be only partially supported and therefore incomplete. That said, adequate discussion regarding these limitations are provided together with dynamic community analysis based on the Ca-bound open structure. The latter supports the main conclusions based on simulations, while providing an indication of potential interaction differences between simulated and fully open conformations. Another limitation is that while the simulations convincingly demonstrate voltage-dependent channel opening as evidenced by pore expansion and conduction of K+ and water through the pore, single channel conductance is underestimated by at least an order of magnitude, as in previous studies of other K+ channels. These quantitative discrepancies suggest that MD simulations may not yet be sufficiently advanced to provide insight into mechanisms underlying the extraordinarily large conductance of BK channels.

Comments on revisions:

My previous questions and concerns have been adequately addressed.

Reviewer #2 (Public review):

Summary:

The manuscript by Jia and Chen addresses the structural basis of voltage-activation of BK channels using computational approaches. Although a number of experimental studies using gating current and patch-clamp recording have analyzed voltage-activation in terms of observed charge movements and the apparent energetic coupling between voltage-sensor movement and channel opening, the structural changes that underlie this phenomenon have been unclear. The present studies use a reduced molecular system comprising the transmembrane portion of the BK channel (i.e. the cytosolic domain was deleted), embedded in a POPC membrane, with either 0 or 750 mV applied across the membrane. This system enabled acquisition of long simulations of 10 microseconds, to permit tracking of conformational changes of the channel. The authors principal findings were that the side chains of R210 and R213 rapidly moved toward the extracellular side of the membrane (by 8 - 10 Å), with greater displacements than any of the other charged transmembrane residues. These movements appeared tightly coupled to movement of the pore-lining helix, pore hydration, and ion permeation. The authors estimate that R210 and R213 contribute 0.25 and 0.19 elementary charges per residue to the gating current, which is roughly consistent with estimates based on electrophysiological measurements that used the full-length channel.

Strengths:

The methodologies used in this work are sound, and these studies certainly contribute to our understanding of voltage-gating of BK channels. An intriguing observation is the strongly coupled movement of the S4, S5, and S6 helices that appear to underlie voltage-dependent opening. Based on Fig 2a-d, the substantial movements of the R210 and R213 side chains occur nearly simultaneously to the S6 movement (between 4 - 5 usec of simulation time). This seems to provide support for a "helix-packing" mechanism of voltage gating in the so-called "non-domain-swapped" voltage-gated K channels.

Weaknesses:

The main limitation is that these studies used a truncated version of the BK channel, and there are likely to be differences in VSD-pore coupling in the context of the full-length channels that will not be resolved in the present work. Nonetheless, the authors provide a strong rationale for their use of the truncated channel, and the results presented will provide a good starting point for future computational studies of this channel.

Reviewer #1 (Public review):

Summary:

This study provides comprehensive instructions for using the chromatophore tracking software, Chromas, to track and analyse the dynamics of large numbers of cephalopod chromatophores across various spatiotemporal scales. This software addresses a long-standing challenge faced by many researchers who study these soft-bodied creatures, known for their remarkable ability to change colour rapidly. The updated software features a user-friendly interface that can be applied to a wide range of applications, making it an essential tool for biologists focused on animal dynamic signalling. It will also be of interest to professionals in the fields of computer vision and image analysis.

Strengths:

This work provides detailed instructions for this tool kit along with examples for potential users to try. The Gitlab inventory hosts the software package, installation documentation, and tutorials, further helping potential users with a less steep learning curve.

Weaknesses:

The evidence supporting the authors' claims is solid, particularly demonstrated through the use of cuttlefish and squid. However, it may not be applicable to all coleoid cephalopods yet, such as octopuses, which have an incredibly versatile ability to change their body forms.

Comments on revisions:

I am pleased to see the more detailed version of this useful tool along with tutorials designed for diverse users who are interested in animal dynamic colouration. This study provides detailed instructions for using the chromatophore tracking software Chromas to track and analyse the dynamics of large numbers of cephalopod chromatophores across various spatiotemporal scales. The software features a user-friendly interface that is highly compelling and can be applied to a wide range of applications.

Reviewer #2 (Public review):

Summary:

The authors developed a computational pipeline named CHROMAS to track and analyze chromatophore dynamics, which provides a wide range of biological analysis tools without requiring the user to write code.

Strengths:

(1) CHROMAS is an integrated toolbox that provides tools for different biological tasks such as: segment, classify, track and measure individual chromatophores, cluster small groups of chromatophores, analyze full-body patterns, etc.

(2) It could be used to investigate different species. The authors have already applied it to analyze the skin of the bobtail squid Euprymna berryi and the European cuttlefish Sepia officinalis.

(3) The tool is open-source and easy to install. The paper describes in detail the experiment requirements, command to complete each task and provides relevant sample figures, which are easy to follow.

Weaknesses:

(1) There are some known limitations for the current version. The users should read the "Discussion" chapter carefully before preparing their experiments and using CHROMAS.

Reviewer #2 (Public review):

Summary:

The authors show that a combination of arginine methyltransferase inhibitors synergize with PARP inhibitors to kill ovarian and triple negative cancer cell lines in vitro and in vivo using preclinical mouse models.

Strengths and weaknesses

The experiments are well-performed, convincing and have the appropriate controls (using inhibitors and genetic deletions) and use statistics.

They identify the DNA damage protein ERCC1 to be reduced in expression with PRMT inhibitors. As ERCC1 is known to be synthetic lethal with PARPi, this provides a mechanism for the synergy. They use cell lines only for their study in 2D as well as xenograph models.

Joint Public Review:

Summary:

How mechanical forces transmitted by blood flow contribute to cardiac development remains incompletely understood. Using the unique advantages of the zebrafish model, Chen et al make the fundamental discovery that endocardial expression of the transcriptional repressor, Id2b, is maintained in endocardial cells by blood flow. Id1b zebrafish mutants fail to form the valve in the atrioventricular canal (AVC) and show reduced myocardial contractility that they suggest is due to impaired calcium transients. Id2b mutants are largely viable during the first 6 months of life until ~20% display cardiomyopathy characterized by visible edema, structural abnormalities, retrograde blood flow, and reduced systolic function and calcium handling. Mechanistically, the authors suggest that flow-mediated expression of Id2b leads to neuregulin 1 (nrg1) upregulation by physically interacting with and sequestering the Tcf3b transcriptional repressor from conserved tcf3b binding sites upstream of nrg1. Overall, this study advances our understanding of flow-mediated endocardial-myocardial crosstalk during heart development.

Strengths:

The strengths of the study are the significance of the biological question being addressed, use of the zebrafish model, data quality, and use of genetic tools. The text is generally well-written and easy to understand.

Weaknesses:

The main weakness that remains is the lack of rigor surrounding the molecular mechanism where the authors suggest that blood flow induces endocardial expression of Id2b, which binds to Tcf3b and sequesters it from binding the Nrg1 promoter to repress transcription. Although good faith efforts were made to bolster their model, the physical interaction between Id2b and Tcf3b is limited to overexpression of tagged proteins in HEK293 cells. Moreover, no mutagenesis was performed on the tcf3b binding sites identified in the nrg1 promoter to learn their importance in vivo.

Reviewer #1 (Public review):

Summary:

This study uses a cell-based computational model to simulate and study T cell development in the thymus. They initially applied this model to assess the effect of the thymic epithelial cells (TECs) network on thymocyte proliferation and demonstrated that increasing TEC size, density, or protrusions increased the number of thymocytes. They postulated and confirmed that this was due to changes in IL7 signalling and then expanded this work to encompass various environmental and cell-based parameters, including Notch signalling, cell cycle duration, and cell motility. Critical outcomes from the computational model were tested in vivo using medaka fish, such as the role of IL-7 signalling and minimal effect of Notch signalling.

Strengths:

The strength of the paper is the use of computational modelling to obtain unique insights into the niche parameters that control T cell development, such as the role of TEC architecture, while anchoring those findings with in vivo experiments. I can't comment on the model itself, as I am not an expert in modelling, however, the conclusions of the paper seem to be well-supported by the model.

Reviewer #2 (Public review):

Summary:

The authors have worked up a ``virtual thymus' using EPISIM, which has already been published. Attractive features of the computational model are stochasticity, cell-to-cell variability, and spatial heterogeneiety. They seek to explore the role of TECs, that release IL-7 which is important in the process of thymocyte division.

In the model, ordinary clones have IL7R levels chosen from a distribution, while `lesioned' clones have an IL7R value set to the maximum. The observation is that the lesioned clones are larger families, but the difference is not dramatic. This might be called a cell-intrinsic mechanism. One promising cell-extrinsic mechanism is mentioned: if a lesioned clone happens to be near a source of IL-7 and begins to proliferate, the progeny can crowd out cells of other clones and monopolise the IL-7 source. The effect will be more noticeable if sources are rare, so is seen when the TEC network is sparse.

Strengths:

Thymic disfunctions are of interest, not least because of T-ALL. New cells are added, one at a time, to simulate the conveyor belt of thymocytes on a background of stationary cells. They are thus able to follow cell lineages, which is interesting because one progenitor can give rise to many progeny.

There are some experimental results in Figures 4,5 and 6. For example, il7 crispant embryos have fewer thymocytes and smaller thymii; but increasing IL-7 availability produces large thymii.

Reviewer #3 (Public review):

Summary:

Tsingos et al. seek to advance beyond the current paradigm that proliferation of malignant cells in T-cell acute lymphoblastic leukemia occurs in a cell-autonomous fashion. Using a computational agent-based model and experimental validation, they show instead that cell proliferation also depends on interaction with thymic epithelial cells (TEC) in the thymic niche. One key finding is that a dense TEC network inhibits the proliferation of malignant cells and favors the proliferation of normal cells, whereas a sparse TEC network leads to rapid expansion of malignant thymocytes.

Strengths:

A key strength of this study is that it combines computational modeling using an agent-based model with experimental work. The original modeling and novel experimental work strengthen each other well. In the agent-based model, the authors also tested the effects of varying a few key parameters of cell proliferation.

Reviewer #1 (Public review):

Petrovic et al. investigate CCR5 endocytosis via arrestin2, with a particular focus on clathrin and AP2 contributions. The study is thorough and methodologically diverse. The NMR titration data are particularly compelling, clearly demonstrating chemical shift changes at the canonical clathrin-binding site (LIELD), present in both the 2S and 2L arrestin splice variants.

To assess the effect of arrestin activation on clathrin binding, the authors compare: truncated arrestin (1-393), full-length arrestin, and 1-393 incubated with CCR5 phosphopeptides. All three bind clathrin comparably, whereas controls show no binding. These findings are consistent with prior crystal structures showing peptide-like binding of the LIELD motif, with disordered flanking regions. The manuscript also evaluates a non-canonical clathrin binding site specific to the 2L splice variant. Though this region has been shown to enhance beta2-adrenergic receptor binding, it appears not to affect CCR5 internalization.

Similar analyses applied to AP2 show a different result. AP2 binding is activation-dependent and influenced by the presence and level of phosphorylation of CCR5-derived phosphopeptides. These findings are reinforced by cellular internalization assays.

In sum, the results highlight splice-variant-dependent effects and phosphorylation-sensitive arrestin-partner interactions. The data argue against a (rapidly disappearing) one-size-fits-all model for GPCR-arrestin signaling and instead support a nuanced, receptor-specific view, with one example summarized effectively in the mechanistic figure.

Reviewer #2 (Public review):

Summary:

Based on extensive live cell assays, SEC, and NMR studies of reconstituted complexes, these authors explore the roles of clathrin and the AP2 protein in facilitating clathrin-mediated endocytosis via activated arrestin-2. NMR, SEC, proteolysis, and live cell tracking confirm a strong interaction between AP2 and activated arrestin using a phosphorylated C-terminus of CCR5. At the same time, a weak interaction between clathrin and arrestin-2 is observed, irrespective of activation.

These results contrast with previous observations of class A GPCRs and the more direct participation by clathrin. The results are discussed in terms of the importance of short and long phosphorylated bar codes in class A and class B endocytosis.

Strengths:

The 15N,1H, and 13C, methyl TROSY NMR and assignments represent a monumental amount of work on arrestin-2, clathrin, and AP2. Weak NMR interactions between arrestin-2 and clathrin are observed irrespective of the activation of arrestin. A second interface, proposed by crystallography, was suggested to be a possible crystal artifact. NMR establishes realistic information on the clathrin and AP2 affinities to activated arrestin, with both kD and description of the interfaces.

Weaknesses:

This reviewer has identified only minor weaknesses with the study.

(1) Arrestin-2 1-418 resonances all but disappear with CCR5pp6 addition. Are they recovered with Ap2Beta2 addition, and is this what is shown in Supplementary Figure 2D?

(2) I don't understand how methyl TROSY spectra of arrestin2 with phosphopeptide could look so broadened unless there are sample stability problems.

(3) At one point, the authors added an excess fully phosphorylated CCR5 phosphopeptide (CCR5pp6). Does the phosphopeptide rescue resolution of arrestin2 (NH or methyl) to the point where interaction dynamics with clathrin (CLTC NTD) are now more evident on the arrestin2 surface?

(4) Once phosphopeptide activates arrestin-2 and AP2 binds, can phosphopeptide be exchanged off? In this case, would it be possible for the activated arrestin-2 AP2 complex to re-engage a new (phosphorylated) receptor?

(5) Did the authors ever try SEC measurements of arrestin-2 + AP2beta2+CCR5pp6 with and without PIP2, and with and without clathrin (CLTC NTD? The question becomes what the active complex is and how PIP2 modulates this cascade of complexation events in class B receptors.

Reviewer #3 (Public review):

Summary:

Overall, this is a well-done study, and the conclusions are largely supported by the data, which will be of interest to the field.

Strengths:

(1) The strengths of this study include experiments with solution NMR that can resolve high-resolution interactions of the highly flexible C-terminal tail of arr2 with clathrin and AP2. Although mainly confirmatory in defining the arr2 CBL 376LIELD380 as the clathrin binding site, the use of the NMR is of high interest (Figure 1). The 15N-labeled CLTC-NTD experiment with arr2 titrations reveals a span from 39-108 that mediates an arr2 interaction, which corroborates previous crystal data, but does not reveal a second area in CLTC-NTD that in previous crystal structures was observed to interact with arr2.

(2) SEC and NMR data suggest that full-length arr2 (1-418) binding with the 2-adaptin subunit of AP2 is enhanced in the presence of CCR5 phospho-peptides (Figure 3). The pp6 peptide shows the highest degree of arr2 activation and 2-adaptin binding, compared to less phosphorylated peptides or not phosphorylated at all. It is interesting that the arr2 interaction with CLTC NTD and pp6 cannot be detected using the SEC approach, further suggesting that clathrin binding is not dependent on arrestin activation. Overall, the data suggest that receptor activation promotes arrestin binding to AP2, not clathrin, suggesting the AP2 interaction is necessary for CCR5 endocytosis.

(3) To validate the solid biophysical data, the authors pursue validation experiments in a HeLa cell model by confocal microscopy. This requires transient transfection of tagged receptor (CCR5-Flag) and arr2 (arr2-YFP). CCR5 displays a "class B"-like behavior in that arr2 is rapidly recruited to the receptor at the plasma membrane upon agonist activation, which forms a stable complex that internalizes into endosomes (Figure 4). The data suggest that complex internalization is dependent on AP2 binding, not clathrin (Figure 5).

Weaknesses:

The interaction of truncated arr2 (1-393) was not impacted by CCR5 phospho-peptide pp6, suggesting the interaction with clathrin is not dependent on arrestin activation (Figure 2). This raises some questions.

Overall, the data are solid, but for added rigor, can these experiments be repeated without tagged receptor and/or arr2? My concern stems from the fact that the stability of the interaction between arr2 and receptor may be related to the position of the tags.

Reviewer #1 (Public review):

Summary:

The study shows, perhaps surprisingly, that human fecal homogenates enhance the invasiveness of Salmonella typhimurium into cells of a swine colonic explant. This effect is only seen with chemotactic cells that express the chemoreceptor Tsr. However, two molecules sensed by Tsr that are present at significant concentrations in the fecal homogenates, the repellent indole and the attractant serine, do not, either by themselves or together at the concentrations in which they are present in the fecal homogenates, show this same effect. The authors then go on to study the conflicting repellent response to indole and attractant response to serine in a number of different in vitro assays.

Strengths:

The demonstration that homogenates of human feces enhance the invasiveness of chemotactic Salmonella Typhimurium in a colonic explant is unexpected and interesting. The authors then go on to document the conflicting responses to the repellent indole and the attractant serine, both sensed by the Tsr chemoreceptor, as a function of their relative concentration and the spatial distribution of gradients.

Weaknesses:

The authors do not identify what is the critical compound or combination of compounds in the fecal homogenate that gives the reported response of increased invasiveness. They show it is not indole alone, serine alone, or both in combination that have this effect, although both are sensed by Tsr and both are present in the fecal homogenates. Some of the responses to conflicting stimuli by indole and serine in the in vitro experiments yield interesting results, but they do little to explain the initial interesting observation that fecal homogenates enhance invasiveness.

Reviewer #3 (Public review):

Summary:

In this manuscript, Franco and colleagues describe careful analyses of Salmonella chemotactic behavior in the presence of conflicting environmental stimuli. By doing so, the authors describe that this human pathogen integrates signals from a chemoattractant and a chemorepellent into an intermediate "chemohalation" phenotype.

Strengths:

The study was clearly well-designed and well-executed. The methods used are appropriate and powerful. The manuscript is very well written, and the analyses are sound. This is an interesting area of research, and this work is a positive contribution to the field.

Weaknesses:

No significant weaknesses noted.

Reviewer #1 (Public review):

Summary:

The authors developed a tool for simulating multiplexed single-cell RNA-seq and ATAC-seq data with various adjustable settings like ambient RNA/DNA rate and sequencing depth. They used the simulated data with different settings to evaluate the performance of many demultiplexing methods. They also proposed a new metric at single-cell level that correlates with the RNA/DNA contamination level.

Strengths:

The simulation tool has a straightforward design and provides adjustability in multiple parameters that have practical relevance, such as sequencing depth and ambient contamination rate. With the growing use of multiplexing in single-cell RNAseq and ATACseq experiments, the tools and results in this paper can guide the experimental design and tool selection for many researchers. The simulation tool also provides a platform for benchmarking newly developed demultiplexing tools.

Weaknesses:

The usefulness of the proposed new metric of "variant consistency" and how it can guide users in selecting demultiplexing methods seems a little unclear. It correlates with the level of ambient RNA/DNA contamination, which makes it look like a metric on data quality. However, it does depend on the exact demultiplexing method, yet it's not clear how it directly connects to the "accuracy" of each demultiplexing method, which is the most important property that users of these methods care about. Since the simulated data has ground truth of donor identities available, I would suggest using the simulated data to show whether "variant consistency" directly indicates the accuracy of each method, especially the accuracy within those "C2" reads.

I also think the tool and analyses presented in this paper need some further clarification and documentation on the details, such as how the cell-type gene and peak probabilities are determined in the simulation, and how doublets from different cell types are handled in the simulation and analysis. A few analyses and figures also need a more detailed description of the exact methods used.

Reviewer #2 (Public review):

Li et al. describe ambisim, a tool with the goal of creating realistic synthetic single-nucleus RNA/ATAC sequencing datasets. It has become standard to pool multiple genetically distinct donors when using single-cell sequencing followed by genotype-based demultiplexing (i.e., using donor single-nucleotide variants to identify specific donor origin). A plethora of tools exist to accomplish this demultiplexing, but advanced tools to create synthetic datasets, and therefore provide definitive benchmarking, are lacking. Ambisim is a well-thought-out simulator that improves upon previous tools available by allowing for modeling of variable ambient contamination proportions and doing so in a genotype-aware fashion. This provides more realistic synthetic datasets that provide challenging scenarios for future demultiplexing tools. The authors use ambisim to benchmark a large number of available and commonly used genotype-free and -dependent demultiplexing tools. They identify the strengths and weaknesses of these tools. They also go on to define a new metric, variant consistency, to further assess demultiplexing performance across tools. Overall, this manuscript provides a useful framework to more thoroughly evaluate future demultiplexing tools, as well as provides rationale for tool selection depending on a user's experimental conditions.

The authors provide measured conclusions that are supported by their findings. There are some aspects that are unclear.

(1) Throughout the manuscript, the figure legends are difficult to understand, and this makes it difficult to interpret the graphs.

(2) Since this is both a new tool and a benchmark, it would be worthwhile in the Discussion to comment on which demultiplexing tools one may want to choose for their dataset, especially given the warning against ensemble methods. From this extensive benchmarking, one may want to choose a tool based on the number of donors one has pooled, the modalities present, and perhaps even the ambient RNA (if it has been estimated previously).

(3) What are the minimal computational requirements for running ambisim? What is the time cost?

Reviewer #1 (Public review):

Summary:

The authors had previously found that a brief social isolation could increase the activity of these neurons, and that manipulation of these neurons could alter social behavior in a social rank dependent fashion. This manuscript explored which of the outputs were responsible for this, identifying the central nucleus of the amygdala as the key output region. The authors identified some discrete behavior changes associated with these outputs, and found that during photostimulation of these outputs, neuronal activity appeared altered in 'social response' neurons. In the revised manuscript, the authors address the comments in a rigorous fashion.

Strengths:

Rigorous analysis of the anatomy. Careful examination of the hetergenous effects on cell activity due to stimulation, linking the physiology with the behavior via photostimulation during recording in vivo.

Weaknesses:

The authors have responded to all of my comments.

Reviewer #2 (Public review):

Summary:

The authors perform a series of studies to follow up on their previous work, which established a role for dorsal raphe dopamine neurons (DRN) in the regulation of social-isolation-induced rebound in mice. In the present study, Lee et. al, use a combination of modern circuit tools to investigate putatively distinct roles of DRN dopamine transporting containing (DAT) projections to the bed nucleus of the stria terminalis (BNST), central amygdala (CeA), and posterior basolateral amygdala (BLP). Notably, they reveal that optogenetic stimulation of distinct pathways confers specific behavioral states, with DRNDAT-BLP driving aversion, DRNDAT-BNST regulating non-social exploratory behavior, and DRNDAT-CeA promoting socialability. A combination of electrophysiological studies and in situ hybridization studies reveal heterogenous dopamine and neuropeptide expression and different firing properties, providing further evidence of pathway-specific neural properties. Lastly, the authors combine optogenetics and calcium imaging to resolve social encoding properties in the DRNDAT-CeA pathway, which correlates observed social behavior to socially engaged neural ensembles.

Collectively, these studies provide an interesting way of dissecting out separable features of a complex multifaceted social-emotional state that accompanies social isolation and the perception of 'loneliness.' The main conclusions of the paper provide an important and interesting set of findings that increase our understanding of these distinct DRN projections and their role in a range of social (e.g., prosocial, dominance), non-social, and emotional behaviors. However, as noted below, the examination of these circuits within a homeostatic framework is limited given that a number of the datasets did not include an isolated condition. The DRNDAT-CeA pathway was investigated with respect to social homeostatic states in the present study for some of the datasets.

Strengths:

(1) The authors perform a comprehensive and elegant dissection of the anatomical, behavioral, molecular, and physiological properties of distinct DRN projections relevant to social, non-social, and emotional behavior, to address multifaceted and complex features of social state.

(2) This work builds on prior findings of isolation-induced changes in DRN neurons and provides a working framework for broader circuit elements that can be addressed across social homeostatic state.

(3) This work characterizes a broader circuit implicated in social isolation and provides a number of downstream targets to explore, setting a nice foundation for future investigation.

(4) The studies account for social rank and anxiety-like behavior in several of the datasets, which are important consideration to the interpretation of social motivation states, especially in male mice with respect to dominance behavior.

Weaknesses:

(1) The conceptual framework of the study is based on the premise of social isolation and perceived 'loneliness' under the framework of social homeostasis, analogous to hunger. In this framework, social isolation should provoke an aversive state and compensatory social contact behavior. In the authors' prior work, they demonstrate synaptic changes in DRN neurons and social rebound following acute social isolation. Thus, the prediction would be that downstream projections also would show state dependent changes as a function of social isolation state (e.g., grouped/socially engaged vs. isolated). In the current paper, a social isolation condition was included for some but not all experiments, which should be considered in the interpretation of the data, specifically within the context of dynamic isolation states.

(2) Figure 1 confirms co-laterals in the BNST and CeA via anatomical tracing studies. The goal of the optogenetic studies is to dissociate functional/behavioral roles of distinct projections. One limitation of optogenetic projection targeting is the possibility of back-propagating action potentials (stimulation of terminals in one region may back-propagate to activate cell bodies, and then afferent projections to other regions), and/or stimulation of fibers of passage. However, this is addressed in the discussion and the present data are convincing, which minimizes the concern.

(3) Sex as a biological variable should be considered in the present data, as included in the discussion.

Reviewer #3 (Public review):

Summary:

The authors investigated the role of dopaminergic neurons (dopamine transporter expressing, DAT) in the dorsal raphe nucleus (DRN) in regulating social and affective behavior through projections to the central nucleus of the amygdala (CeA), bed nucleus of the stria terminalis (BNST), and the posterior subdivision of the basolateral amygdala. The largest effect observed was in the DRN-DAT projections to the CeA. Augmenting previously published results from this group (Matthews et al., 2016), the comprehensive behavioral analysis relative to social dominance, gene expression analysis, electrophysiological profiling, and in vivo imaging provides novel insights into how DRN-DAT projections to the CeA influence the engagement of social behavior in the contexts of group housed and socially isolated mice.

Strengths:

Correlational analysis with social dominance is a nice addition to the study. The overall computational analyses performed are well-designed and rigorous.

Weaknesses:

(1) Analysis of dopamine receptor expression did not include Drd3, Drd4, or Drd5 which may provide more insights into how dopamine modulates downstream targets. This is particularly relevant to the BNST projection in which the densest innervation did not robustly co-localize with the expression of either Drd1 or Drd2. It is also possible that dopamine release from DRN-DAT neurons in any or all of these structures in modulating neurotransmitter release from inputs to these regions that contain D2 receptors on their terminals.

(2) Although not the focus of this study, without pharmacological blockade of dopamine receptors, it is not possible to assess what the contribution of dopamine is to the behavioral outcomes. Given the co-release of glutamate and GABA from these neurons it is possible that dopamine plays only a marginal role in the functional connectivity of DRN-DAT neurons.

(3) Photostimulation parameters used during the behavioral studies (8 pulses of light delivered at 30 Hz for several minutes) could lead to confounding results limiting data interpretation. As shown in Figure 6J, 8 pulses of light delivered at 30 Hz results in a significant attenuation of the EPSC amplitude in the BLP and CeA projection. Thus, prolonged stimulation could lead to significant synaptic rundown resulting in an overall suppression of connectivity in the later stages of the behavioral analyses.

Comments on revisions:

No further issues have been identified.

Reviewer #1 (Public review):

Summary:

This paper reports an intracranial SEEG study of speech coordination, where participants synchronize their speech output with a virtual partner that is designed to vary its synchronization behavior. This allows the authors to identify electrodes throughout the left hemisphere of the brain that have activity (both power and phase) that correlates with the degree of synchronization behavior. They find that high-frequency activity in secondary auditory cortex (superior temporal gyrus) is correlated to synchronization, in contrast to primary auditory regions. Furthermore, activity in inferior frontal gyrus shows a significant phase-amplitude coupling relationship that is interpreted as compensation for deviation from synchronized behavior with the virtual partner.

Strengths:

(1) The development of a virtual partner model trained for each individual participant, which can dynamically vary its synchronization to the participant's behavior in real time, is novel and exciting.

(2) Understanding real-time temporal coordination for behaviors like speech is a critical and understudied area.

(3) The use of SEEG provides the spatial and temporal resolution necessary to address the complex dynamics associated with the behavior.

(4) The paper provides some results that suggest a role for regions like IFG and STG in the dynamic temporal coordination of behavior both within an individual speaker and across speakers performing a coordination task.

Reviewer #2 (Public review):

Summary:

This paper investigates the neural underpinnings of an interactive speech task requiring verbal coordination with another speaker. To achieve this, the authors recorded intracranial brain activity from the left (and to a lesser extent, the right) hemisphere in a group of drug-resistant epilepsy patients while they synchronised their speech with a 'virtual partner'. Crucially, the authors were able to manipulate the degree of success of this synchronisation by programming the virtual partner to either actively synchronise or desynchronise their speech with the participant, or else to not vary its speech in response to the participant (making the synchronisation task purely one-way). Using such a paradigm, the authors identified different brain regions that were either more sensitive to the speech of the virtual partner (primary auditory cortex), or more sensitive to the degree of verbal coordination (i.e. synchronisation success) with the virtual partner (left secondary auditory cortex and bilateral IFG). Such sensitivity was measured by (1) calculating the correlation between the index of verbal coordination and mean power within a range of frequency bands across trials, and (2) calculating the phase-amplitude coupling between the behavioural and brain signals within single trials (using the power of high-frequency neural activity only). Overall, the findings help to elucidate some of the brain areas involved in interactive speaking behaviours, particularly highlighting high-frequency activity of the bilateral IFG as a potential candidate supporting verbal coordination.

Strengths:

This study provides the field with a convincing demonstration of how to investigate speaking behaviours in more complex situations that share many features with real-world speaking contexts e.g. simultaneous engagement of speech perception and production processes, the presence of an interlocutor and the need for inter-speaker coordination. The findings thus go beyond previous work that has typically studied solo speech production in isolation, and represent a significant advance in our understanding of speech as a social and communicative behaviour. It is further an impressive feat to develop a paradigm in which the degree of cooperativity of the synchronisation partner can be so tightly controlled; in this way, this study combines the benefits of using pre-recorded stimuli (namely, the high degree of experimental control) with the benefits of using a live synchronisation partner (allowing the task to be truly two-way interactive, an important criticism of other work using pre-recorded stimuli). A further key strength of the study lies in its employment of stereotactic EEG to measure brain responses with both high temporal and spatial resolution, an ideal method for studying the unfolding relationship between neural processing and this dynamic coordination behaviour.

Weaknesses:

One limitation of the current study is the relatively sparse coverage of the right hemisphere by the implanted electrodes (91 electrodes in the right compared to 145 in the left). Of course, electrode location is solely clinically motivated, and so the authors did not have control over this. In a previous version of this article, the authors therefore chose not to include data from the right hemisphere in reported analyses. However, after highlighting previous literature suggesting that the right hemisphere likely has high relevance to verbal coordination behaviours such as those under investigation here, the authors have now added analyses of the right hemisphere data to the results. These confirm an involvement of the right hemisphere in this task, largely replicating left hemisphere results. Some hemispheric differences were found in responses within the STG; however, interpretation should be tempered by an awareness of the relatively sparse coverage of the right hemisphere meaning that some regions have very few electrodes, resulting in reduced statistical power.

Reviewer #3 (Public review):

Summary:

In microbiology, accurately characterizing microbial populations and communities is essential. One widely used approach is to measure the absolute or relative abundance of microbial species. Recent research in microbial ecology, for instance, has shown that even genetically identical hosts exposed to the same microbial pool can develop very different gut microbiota, largely due to random colonization events. This study builds on that idea but adds a valuable layer: it suggests that some of the observed variability might actually result from experimental noise, specifically the randomness introduced by dilution and plate counting techniques. To address this, the authors introduce REPOP, a new tool designed to improve the quantification of microbial populations by explicitly accounting for the inherent stochasticity in these methods. They test REPOP using both simulated and experimental datasets, showing how it can help recover meaningful trends.

Strengths:

Overall, this paper is a good contribution to the field. The motivation is clear: improving our ability to quantify microbial populations is crucial for many research areas. The authors make a strong case that ignoring experimental noise is no longer acceptable, and they offer a well-argued solution. The manuscript is well-written and easy to follow, and the logic behind REPOP is convincingly laid out. The use of simulated data is especially valuable, as it allows the authors to test whether the method can recover known inputs, an important validation step. Even with experimental data, where true values are unknown, the method seems to behave in a reasonable and expected way, which is reassuring. All in all, this is an important step forward in how we quantify microbial populations.

Weaknesses:

While the study is promising, there are a few areas where the paper could be strengthened to increase its impact and usability. First, the extent to which dilution and plating introduce noise is not fully explored. Could this noise significantly affect experimental conclusions? And under what conditions does it matter most? Does it depend on experimental design or specific parameter values? Clarifying this would help readers appreciate when and why REPOP should be used. Second, more practical details about the tool itself would be very helpful. Simply stating that it is available on GitHub may not be enough. Readers will want to know what programming language it uses, what the input data should look like, and ideally, see a step-by-step diagram of the workflow. Packaging the tool as an easy-to-use resource, perhaps even submitting it to CRAN or including example scripts, would go a long way, especially since microbiologists tend to favor user-friendly, recipe-like solutions. Third, it would be great to see the method tested on existing datasets, such as those from Nic Vega and Jeff Gore (2017), which explore how colonization frequency impacts abundance fluctuation distributions. Even if the general conclusions remain unchanged, showing that REPOP can better match observed patterns would strengthen the paper's real-world relevance. Lastly, it would be helpful for the authors to briefly discuss the limitations of their method, as no approach is without its constraints. Acknowledging these would provide a more balanced and transparent perspective.

Reviewer #2 (Public review):

Summary:

Microbial population abundances are regularly estimated by multiplying plate counts by dilution factors, with inferences made about sample heterogeneity without taking into account heterogeneity generated through dilution and plating methods. The authors have developed REPOP, a method for disentangling methodological stochasticity from ecological heterogeneity using a Bayesian framework. They present three models: a unimodal distribution, a multimodal distribution, and a multimodal distribution that incorporates a colony count cutoff. They use a combination of simulated and experimental data to show the effectiveness of the REPOP method in resolving true microbial population distributions.

Strengths:

Overall, this paper addresses a significant issue in microbial ecology and reliably demonstrates that the REPOP method improves upon current methods of estimating microbial population heterogeneity, particularly with simulation data. The three models presented build upon each other and are discussed in a way that is fairly accessible to a broad audience. The authors also show that leveraging the information provided by non-countable plates is important. Additionally, the authors address the potential for extending this method to other sources of methodological stochasticity that may occur in microbial plating. However, it does seem that they could extend this further by discussing ways that this method could be applied to non-microbial systems, allowing this work to appeal to a broader audience.

Weaknesses:

A more thorough discussion of when and by how much estimated microbial population abundance distributions differ from the ground truth would be helpful in determining the best practices for applying this method. Not only would this allow researchers to understand the sampling effort necessary to achieve the results presented here, but it would also contextualize the experimental results presented in the paper. Particularly, there is a disconnect between the discussion of the large sample sizes necessary to achieve accurate multimodal distribution estimates and the small sample sizes used in both experiments.

Reviewer #1 (Public review):

Summary:

The authors developed a novel theoretical/computational procedure to count bacterial populations without introducing artificial randomness effects due to dilution. Surprisingly, this very important aspect of studies of bacterial systems has been overlooked. The proposed method provides a simple and transparent approach to eliminate the randomness of bacterial accounting procedures, allowing now to fully concentrate on the intrinsic effects of the studied systems.

Strengths:

A very simple and clear procedure is introduced and explained in full detail. This elegant approach finds an excellent compromise between mathematical rigor and computational efficiency, which is important for practical applications. The provided examples are convincing beyond a doubt, clearly indicating the potential strong impact of the proposed framework. Various complications and possible issues are also discussed and analyzed. This seems to be a very powerful novel method that should significantly advance the analysis of complex biological systems.

Weaknesses:

The only minor weakness that I found is the assumption of independence of bacterial species, which is expressed as the well-stirred approximation. One could imagine that bacterial species might cooperate, leading to non-uniform distributions that are real. How to distinguish such situations?

I believe that this method can be extended to determine if this is the case or not before the application. For example, if the bacteria species are independent of each other and one can use the binomial distributions, then the Fano factor would be proportional to the overall relative fraction of bacterial species. Maybe a simple test can be added to test it before the application of REPOP. However, I believe that this is a minor issue.

Reviewer #1 (Public review):

The authors conducted an fMRI study to investigate the neural effects of sustaining attention to areas of different sizes. Participants were instructed to attend to alphanumeric characters arranged in a circular array. The size of attention field was manipulated in four levels, ranging from small (18 deg) to large (162 deg). They used a model-based method to visualize attentional modulation in early visual cortex V1 to V3, and found spatially congruent modulations of the BOLD response, i.e., as the attended area increased in size, the neural modulation also increased in size in the visual cortex. They suggest that this result is a neural manifestation of the zoom-lens model of attention and that the model-based method can effectively reconstruct the neural modulation in the cortical space.

The study is well-designed with sophisticated and comprehensive data analysis. The results are robust and show strong support for a well-known model of spatial attention, the zoom-lens model. Overall, I find the results interesting and useful for the field of visual attention research.

Comments on revisions:

The authors have addressed my previous comments satisfactorily. I would encourage the authors to make data and code publicly available, which appears to be the custom in this era.

Reviewer #2 (Public review):

Summary:

The study in question utilizes functional magnetic resonance imaging (fMRI) to dynamically estimate the locus and extent of covert spatial attention from visuocortical activity. The authors aim to address an important gap in our understanding of how the size of the attentional field is represented within the visual cortex. They present a novel paradigm that allows for the estimation of the spatial tuning of the attentional field and demonstrate the ability to reliably recover both the location and width of the attentional field based on BOLD responses.

Strengths:

(1) Innovative Paradigm: The development of a new approach to estimate the spatial tuning of the attentional field is a significant strength of this study. It provides a fresh perspective on how spatial attention modulates visual perception.

(2) Refined fMRI Analysis: The use of fMRI to track the spatial tuning of the attentional field across different visual regions is methodologically rigorous and provides valuable insights into the neural mechanisms underlying attentional modulation.

(3) Clear Presentation: The manuscript is well-organized, and the results are presented clearly, which aids in the reader's comprehension of the complex data and analyses involved.

Weaknesses:

(1) Lack of Neutral Cue Condition: The study does not include a neutral cue condition where the cue width spans 360{degree sign}, which could serve as a valuable baseline for assessing the BOLD response enhancements and diminishments in both attended and non-attended areas.

(2) Clarity on Task Difficulty Ratios: The explicit reasoning for the chosen letter-to-number ratios for various cue widths is not detailed. Ensuring clarity on these ratios is crucial, as it affects the task difficulty and the comparability of behavioral performance across different cue widths. It is essential that observed differences in behavior and BOLD signals are attributable solely to changes in cue width and not confounded by variations in task difficulty.

Comments on revisions:

(1) Please standardize the naming of error metrics across Figures 4-6 to improve clarity (e.g., "angular error" (Figure 4), "|angular error|" (Figure 5), and "absolute error" (Figure 6) appear to refer to the same measure). This inconsistency is also present in the main text.

(2) Consider briefly mentioning the baseline offset in Lines 179-186. It is included in Figures 4-7 and serves as a reference for interpreting attentional modulation alongside gain. Introducing it with other model parameters would improve clarity.

(3) It may be valuable to examine BOLD responses in unattended visual regions. As shown in Figure 2a, suppression patterns (e.g., the most negative responses) appear to vary in extent and distribution with attentional cue width. Analyzing these unattended regions may offer a more complete view of how attention shapes the spatial profile of cortical activity.

Reviewer #1 (Public review):

Henshall et al. delete the highly abundant merozoite surface protein PfMSP2 from two Plasmodium falciparum laboratory lines (3D7 and Dd2) using CRISPR-Cas9. Parasites lacking MSP2 replicate and invade red cells normally, opposing the experimental history that suggests MSP2 is essential. Unexpectedly, the knock-outs become more susceptible to several inhibitory antibodies - most strikingly those that target the apical antigen AMA1-while antibodies to other surface or secreted proteins are largely unaffected. Recombinant MSP2 added in vitro can dampen AMA1-antibody binding, supporting a "conformational masking" model. The reported data suggest that MSP2 helps shield key invasion ligands from host antibodies and may itself be a double-edged vaccine target.

Reviewer #2 (Public review):

Summary:

The authors were trying to establish the role of Plasmodium falciparum surface protein 2 in merozoite biology, specifically the process of erythrocyte invasion.

Strengths:

The major strengths of the manuscript are in the Plasmodium falciparum genetic and phenotyping approaches. PfMSP2 knockouts are made in two different strains, which is important as it is known that invasion pathways can vary between strains, but is a level of comprehensiveness that is not always delivered in P. falciparum genetic studies. The knockout strains are characterised very thoroughly using multiple different assays, and the authors should be commended for publishing a good deal of negative data, where no phenotype was detected. This is not always done, but is very helpful for the field and reduces the potential for experimental redundancy, i.e., others repeating work that has already been performed but never published. The quality of the writing, referencing, and figures is also generally strong, although a few minor typos and technical comments on presentation have been communicated to the authors.

Weaknesses:

There are, however, some areas that are weaker.

(1) The section describing Laverania and avian Plasmodium MSP2 comparison is a lengthy section and could be told much more concisely for clarity in delivering the key message, i.e., that conservation in distantly related Plasmodium species could indicate an important function. The identification of MSP2-like genes in avian Plasmodium species was highlighted previously in the referenced Escalante paper, so it is not entirely novel, although this paper goes into more detailed characterisation of the extent of conservation. Overall, this section takes up much more space in the manuscript than is merited by the novelty and significance of the findings.

(2) Characterisation of the knockout strains is generally thorough, though relatively few interactions were followed by live microscopy (Figures 3E-H). A minimum of 30 merozoites were followed in each assay (although the precise number is not specified in the figure or legend), but there are intriguing trends in the data that could potentially have become significant if n was increased.

(3) The comparative RNAseq data is interesting, but is not followed up to any significant degree. Multiple transcripts are up-regulated in the absence of PfMSP2, but they are largely dismissed because they are genes of unknown function, not previously linked to invasion, or lack an obvious membrane anchor. Having gone to the lengths of exploring potentially compensatory changes in gene expression, it is disappointing not to validate or explore the hits that result.

(4) Given the abundance of PfMSP2 on the merozoite surface, it would have been interesting to see whether the knockout lines have any noticeable difference in surface composition, as viewed by electron microscopy, although, of course, this experiment relies on access to the appropriate facilities.

(5) One of the key findings is that deletion of PfMSP2 increases inhibition by some antibodies/nanobodies (some anti-CSS2, some anti-AMA1) but not others (anti-EBA/RH, anti-EBA175, anti-Rh5, anti-TRAMP, some anti-CSS2, some anti-AMA1). The data supporting these changes in inhibition are solid, but the selectivity of the effect (only a few antibodies, and generally those targeting later stages in invasion) is not really discussed in any detail. Do the authors have a hypothesis for this selectivity? The authors make attempts to explore the mechanisms for this antibody-masking (Figure 7), but the data is less solid. Surface Plasmon Resonance was non-conclusive, while an ELISA approach co-incubating MSP2 and anti-AMA1 antibodies to wells coated with AMA1 lacks appropriate controls (eg, including other merozoite proteins in similar experiments).

Overall, the claim that PfMSP2 is non-essential for in vitro growth is well justified and is an important contribution to the field. The impact of PfMSP2 deletion on antibody inhibition (which is highlighted in the title of the manuscript) and the mechanism behind it is much less definitive, but does open up an interesting area for further investigation, with more work to be done.

Reviewer #3 (Public review):

The central issue for evaluating the overfilling hypothesis is the identity of the mechanism that causes the very potent (>80% when inter pulse is 20 ms), but very quickly reverting (< 50 ms) paired pulse depression (Fig 1G, I). To summarize: the logic for overfilling at local cortical L2/3 synapses depends critically on the premise that probability of release (pv) for docked and fully primed vesicles is already close to 100%. If so, the reasoning goes, the only way to account for the potent short-term enhancement seen when stimulation is extended beyond 2 pulses would be by concluding that the readily releasable pool overfills. However, the conclusion that pv is close to 100% depends on the premise that the quickly reverting depression is caused by exocytosis dependent depletion of release sites, and the evidence for this is not strong in my opinion. Caution is especially reasonable given that similarly quickly reverting depression at Schaffer collateral synapses, which are morphologically similar, was previously shown to NOT depend on exocytosis (Dobrunz and Stevens 1997). Note that the authors of the 1997 study speculated that Ca2+-channel inactivation might be the cause, but did not rule out a wide variety of other types of mechanisms that have been discovered since, including the transient vesicle undocking/re-docking (and subsequent re-priming) reported by Kusick et al (2020), which seems to have the correct timing.

In an earlier round of review, I suggested raising extracellular Ca2+, to see if this would increase synaptic strength. This is a strong test of the authors' model because there is essentially no room for an increase in synaptic strength. The authors have now done experiments along these lines, but the result is not clear cut. On one hand, the new results suggest an increase in synaptic strength that is not compatible with the authors' model; technically the increase does not reach statistical significance, but, likely, this is only because the data set is small and the variation between experiments is large. Moreover, a more granular analysis of the individual experiments seems to raise more serious problems, even supporting the depletion-independent counter hypothesis to some extent. On the other hand, the increase in synaptic strength that is seen in the newly added experiments does seem to be less at local L2/3 cortical synapses compared to other types of synapses, measured by other groups, which goes in the general direction of supporting the critical premise that pv is unusually high at L2/3 cortical synapses. Overall, I am left wishing that the new data set were larger, and that reversal experiments had been included as explained in the specific points below.

Specific Points:

(1) One of the standard methods for distinguishing between depletion-dependent and depletion-independent depression mechanisms is by analyzing failures during paired pulses of minimal stimulation. The current study includes experiments along these lines showing that pv would have to be extremely close to 1 when Ca2+ is 1.25 mM to preserve the authors' model (Section "High double failure rate ..."). Lower values for pv are not compatible with their model because the k1 parameter already had to be pushed a bit beyond boundaries established by other types of experiments. The authors now report a mean increase in synaptic strength of 23% after raising Ca to 2.5 mM. The mean increase is not quite statistically significant, but this is likely because of the small sample size. I extracted a 95% confidence interval of [-4%, +60%] from their numbers, with a 92% probability that the mean value of the increase in the full population is > 5%. I used the 5% value as the greatest increase that the model could bear because 5% implies pv < 0.9 using the equation from Dodge and Rahamimoff referenced in the rebuttal. My conclusion from this is that the mean result, rather than supporting the model, actually undermines it to some extent. It would have likely taken 1 or 2 more experiments to get above the 95% confidence threshold for statistical significance, but this is ultimately an arbitrary cut off.

(2) The variation between experiments seems to be even more problematic, at least as currently reported. The plot in Figure 3-figure supplement 3 (left) suggests that the variation reflects true variation between synapses, not measurement error. And yet, synaptic strength increased almost 2-fold in 2 of the 8 experiments, which back extrapolates to pv < 0.2. If all of the depression is caused by depletion as assumed, these individuals would exhibit paired pulse facilitation, not depression. And yet, from what I can tell, the individuals depressed, possibly as much as the synapses with low sensitivity to Ca2+, arguing against the critical premise that depression equals depletion, and even arguing - to some extent - for the counter hypothesis that a component of the depression is caused by a mechanism that is independent of depletion. I would strongly recommend adding an additional plot that documents the relationship between the amount of increase in synaptic strength after increasing extracellular Ca2+ and the paired pulse ratio as this seems central.

(3) Decrease in PPR. The authors recognize that the decrease in the paired-pulse ratio after increasing Ca2+ seems problematic for the overfilling hypothesis by stating: "Although a reduction in PPR is often interpreted as an increase in pv, under conditions where pv is already high, it more likely reflects a slight increase in pocc or in the number of TS vesicles, consistent with the previous estimates (Lin et al., 2025)." I looked quickly, but did not immediately find an explanation in Lin et al 2025 involving an increase in pocc or number of TS vesicles, much less a reason to prefer this over the standard explanation that reduced PPR indicates an increase in pv. The authors should explain why the most straightforward interpretation is not the correct one in this particular case to avoid the appearance of cherry picking explanations to fit the hypothesis.

(4) The authors concede in the rebuttal that mean pv must be < 0.7, but I couldn't find any mention of this within the manuscript itself, nor any explanation for how the new estimate could be compatible with the value of > 0.99 in the section about failures.

(5) Although not the main point, comparisons to synapses in other brain regions reported in other studies might not be accurate without directly matching experiments. As it is, 2 of 8 synapses got weaker instead of stronger, hinting at possible rundown, but this cannot be assessed because reversibility was not evaluated. In addition, comparing axons with and without channel rhodopsins might be problematic because the channel rhodopsins might widen action potentials.

(6) Perhaps authors could double check with Schotten et al about whether PDBu does/does not decrease the latency between osmotic shock and transmitter release. This might be an interesting discrepancy, but my understanding is that Schotten et al didn't acquire information about latency because of how the experiments were designed.

(7) The authors state: "These data are difficult to reconcile with a model in which facilitation is mediated by Ca2+-dependent increases in pv." However, I believe that discarding the premise that depression is always caused by depletion would open up wide range of viable possibilities.

Reviewer #1 (Public review):

Summary:

Authors benchmarked five IBD detection methods (hmmIBD, isoRelate, hap-IBD, phasedIBD, and Refined IBD) in Plasmodium falciparum using simulated and empirical data. Plasmodium falciparum has a mutation rate similar to that of humans but a much higher recombination rate and lower SNP density. Thus, the authors evaluated how recombination rate and marker density affect IBD segment detection. Next, they performed parameter optimization for Plasmodium falciparum and benchmarked the robustness of downstream analyses (selection detection and Ne inference) using IBD segments detected by each method. They also tracked the computational efficiency of these methods. The authors' work is valuable for the tested species, and the analyses presented support their claim that users should be cautious when calling IBD in contexts of low SNP density and high recombination rate.

Strengths:

The study design is convincing and well-structured. The authors chose to use P. falciparum, which presents an interesting case due to its high recombination rate and a mutation rate similar to that of humans. The authors note that despite the widespread use of IBD for genomic surveillance, comprehensive evaluation of these methods in high-recombination, low-marker-density contexts has been lacking. Furthermore, they also examined the performance of IBD detection methods developed specifically for P. falciparum, and evaluated it with phased data which broadened the applicability of the work.

Weaknesses:

The authors thoughtfully addressed our prior concerns by 1) expanding the simulations; 2) updating figures and methods for clarity; and 3) more clearly framing the broader utility of their benchmarking effort. These updates strengthen the manuscript and make the relevance of their findings beyond Plasmodium falciparum more apparent.

More specifically:

The authors added three full replicates per simulation scenario and updated figures to reflect uncertainty at relevant levels, which addresses earlier concerns about reproducibility. The limited number of replicates is due to computational intensity. In the future, broader generalizability and deeper exploration of variance in benchmarking accuracy across parameter space would further strengthen the conclusions/generalizability. The author's also emphasized that, while the study is centered on Plasmodium falciparum, the benchmarking framework, not the parameters, are broadly applicable to other sexually recombining species. Lastly, they extensively updated multiple figures to include simulation models, results from simulation replicates, and improved the figures from the previous version of the manuscript.

Reviewer #2 (Public review):

Summary:

Guo et al. benchmarked and optimized methods for detecting Identity-By-Descent (IBD) segments in Plasmodium falciparum (Pf) genomes, which are characterized by high recombination rates and low marker density. Their goal was to address the limitations of existing IBD detection tools, which were primarily developed for human genomes and do not perform well in the genomic context of highly recombinant genomes. They first analysed various existing IBD callers, such as hmmIBD, isoRelate, hap-IBD, phased-IBD, and refinedIBD. They focused on the impact of recombination on the accuracy, which was calculated based on two metrics, the false negative rate and the false positive rate. The results suggest that high recombination rates significantly reduce marker density, leading to higher false negative rates for short IBD segments. This effect compromises the reliability of IBD-based downstream analyses, such as effective population size (Ne) estimation.<br /> They showed that the best tool for IBD detection in Pf is hmmIBD, because it has relatively low FN/FP error rates and is less biased for relatedness estimates. However, this method is less computationally efficient.<br /> Their suggestion is to optimize human-oriented IBD methods and use hmmIBD only for the estimation of Ne.

Strengths:

Although I am not an expert on Plasmodium falciparum genetics, I believe the authors have developed a valuable benchmarking framework tailored to the unique genomic characteristics of this species. Their framework enables a thorough evaluation of various IBD detection tools for non-human data, such as high recombination rates and low marker density, addressing a key gap in the field.

This study provides a comparison of multiple IBD detection methods, including probabilistic approaches (hmmIBD, isoRelate) and IBS-based methods (hap-IBD, Refined IBD, phased IBD). This comprehensive analysis offers researchers valuable guidance on the strengths and limitations of each tool, allowing them to make informed choices based on specific use cases. I think this is important beyond the study of Pf.

The authors highlight how optimized IBD detection can help identify signals of positive selection, infer effective population size (Ne), and uncover population structure.

They demonstrate the critical importance of tailoring analytical tools to suit the unique characteristics of a species. Moreover, the authors provide practical recommendations, such as employing hmmIBD for quality-sensitive analyses and fine-tuning parameters for tools originally designed for non-P. falciparum datasets before applying them to malaria research.

Overall, this study represents a meaningful contribution to both computational biology and malaria genomics, with its findings and recommendations likely to have an impact on the field.

Weaknesses:

One weakness of the study is the lack of emphasis on the broader importance of studying Plasmodium falciparum as a critical malaria-causing organism. Malaria remains a significant global health challenge, causing hundreds of thousands of deaths annually.

While the study provides a thorough technical evaluation of IBD detection methods and their application to Pf, it does not adequately connect these findings to the broader implications for malaria research and control efforts. Additionally, the discussion on malaria and its global impact could have framed the study in a more accessible and compelling way, making the importance of these technical advances clearer to a broader audience, including researchers and policymakers in the fight against malaria. In the revised version, the authors have placed greater emphasis on this aspect, while still maintaining the methodological focus of the paper.

Reviewer #1 (Public review):

Summary:

Gekko, Nomura et al., show that Drp1 elimination in zygotes using the Trim-Away technique leads to mitochondrial clustering and uneven mitochondrial partitioning during the first embryonic cleavage, resulting in embryonic arrest. They monitor organellar localization and partitioning using specific targeted fluorophores. They also describe the effects of mitochondrial clustering in spindle formation and the detrimental effect of uneven mitochondrial partitioning to daughter cells.

Strengths:

The authors have gathered solid evidence for the uneven segregation of mitochondria upon Drp1 depletion through different means: mitochondrial labelling, ATP labelling and mtDNA copy number assessment in each daughter cell. Authors have also characterised the defects in cleavage mitotic spindles upon Drp1 loss

Weaknesses:

This study convincingly describes the phenotype seen upon Drp1 loss. Further studies should be conducted to elucidate the mechanism by which Drp1 ensures even mitochondrial partitioning during the first embryonic cleavage.

Reviewer #2 (Public review):

Gekko et al investigate the impact of perturbing mitochondrial during early embryo development, through modulation of the mitochondrial fission protein Drp1 using Trim-Away technology. They aimed to validate a role for mitochondrial dynamics in modulating chromosomal segregation, mitochondrial inheritance and embryo development and achieve this through the examination of mitochondrial and endoplasmic reticulum distribution, as well as actin filament involvement, using targeted plasmids, molecular probes and TEM in pronuclear stage embryos through the first cleavages divisions. Drp1 deletion perturbed mitochondrial distribution, leading to asymmetric partitioning of mitochondria to the 2-cell stage embryo, prevented appropriate chromosomal segregation and culminated in embryo arrest. Resultant 2-cell embryos displayed altered ATP, mtDNA and calcium levels. Microinjection of Drp1 mRNA partially rescued embryo development. A role for actin filaments in mitochondrial inheritance is described, however the actin-based motor Myo19 does not appear to contribute.

Overall, this study builds upon their previous work and provides further support for a role of mitochondrial dynamics in mediating chromosomal segregation and mitochondrial inheritance. In particular, Drp1 is required for the redistribution of mitochondria to support symmetric partitioning and ongoing development.

Strengths:

The study is well designed, the methods are appropriate and the results are clearly presented. The findings are nicely summarised in a schematic.

The addition of further quantification, including mitochondrial cluster size, elongation/aspect ratio and ROS, as requested by the reviewers, has provided further evidence for the impact of Drp1 depletion on mitochondrial morphology and function.

Understanding the role of mitochondria in binucleation and mitochondrial inheritance is of clinical relevance for patients undergoing infertility treatment, particularly those undergoing mitochondrial replacement therapy.

Weaknesses:

The only remaining weakness is that the authors have not undertaken additional experiments to clarify any role for mitochondrial transport following Drp1 depletion.

Reviewer #3 (Public review):

Why mitochondria are finely maintained in the female germ cell (oocyte), zygotes, and preimplantation embryos? Mitochondrial fusion seems beneficial in somatic cells to compensate for unhealthy mitochondria, for example, mitochondria with mutated mtDNA that potentially defuel the respiratory activity if accumulated above a certain threshold. However, in the germ cells, it may rather increase the risk of transmitting mutated mtDNA to the next generation. Also, finely maintained mitochondria would also be beneficial for efficient removal when damaged, as the authors briefly discussed. Due in part to the limited suitable model, physiological role of mitochondrial fission in embryos were obscure. In this study, authors demonstrated that mitochondrial fission prevents multiple adverse outcomes, especially including the aberrant demixing of parental genome (a clinical phenotype of human embryos) in zygotic stage. Thus, this study would be also of clinical importance that could contribute by proposing a novel mechanism.

The authors have adequately indicated the limitations at each of the specific points. The revisions the authors made have consolidated their conclusion, thus still, making this study an excellent one.

Reviewer #1 (Public review):

Summary:

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

Strengths:

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

Weaknesses:

Extrapolating from the core HbYX activating motif shared by Archae and Eukaryotes to the specific operations of gate opening, which is more elaborate in eukaryotes, may prove challenging.

Reviewer #2 (Public review):

Summary:

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

Strengths:

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

Comments on revisions:

The authors have addressed all my questions.

Reviewer #1 (Public review):

N'Guessan et al have improved the report of their study of expression QTL (eQTL) mapping in yeast using single cells. The authors make use of advances in single cell RNAseq (scRNAseq) in yeast to increase the efficiency with which this type of analysis can be undertaken. Building on prior research led by the senior author that entailed genotyping and fitness profiling of almost 100,000 cells derived from a cross between two yeast strains (BY and RM) they performed scRNAseq on a subset of ~5% (n = 4,489) individual cells. To address the sparsity of genotype data in the expression profiling they used a Hidden Markov Model (HMM) to infer genotypes and then identify the most likely known lineage genotype from the original dataset. To address the relationship between variance in fitness and gene expression the authors partition the variance to investigate the sources of variation. They then perform eQTL mapping and study the relationship between eQTL and fitness QTL identified in the earlier study.

This paper seeks to address the question of how quantitative trait variation and expression variation are related. scRNAseq represents an appealing approach to eQTL mapping as it is possible to simultaneously genotype individual cells and measure expression in the same cell. As eQTL mapping requires large sample sizes to identify statistical relationships, the use of scRNAseq is likely to dramatically increase the statistical power of such studies. However, there are several technical challenges associated with scRNAseq and the authors' study is focused on addressing those challenges. The authors have successfully demonstrated their stated goal of developing, and illustrating the benefit of, a one-pot scRNA-seq experiment and analysis for eQTL mapping.

Reviewer #2 (Public review):

This work describes the single-cell expression profiling of thousands of cells of recombinant genotypes from a model natural-variation system, a cross between two divergent yeast strains.

I appreciate the addition of lines 282-291, which now makes the authors' point about one advantage of the single-cell technique for eQTL mapping clearly: the authors don't need to normalize for culture-to-culture variation the way standard bulk methods do (e.g. in Albert et al., 2018 for the current yeast cross), and without this normalization, they can integrate analyses of expression with those of estimates of growth behaviors from the abundance of a genotype in the pool. The main question the manuscript addresses with the latter, in Figure 3, is how much variation in growth appears to have nothing to do with expression, for which the answer the authors given is 30%. I agree that this represents a novel finding. The caveats are (1) the particular point will perhaps only be interesting to a small slice of the eQTL research community; (2) the authors provide no statistical controls/error estimate or independent validation of the variance partitioning analysis in Figure 3, and (3) the authors don't seem to use the single-cell growth/fitness estimates for anything else, as Figure 4 uses loci mapped to growth from a previously published, standard culture-by-culture approach.

Reviewer #1 (Public review):

Summary:

This is an interesting follow-up to a paper published in Human Molecular Genetics reporting novel roles in corticogenesis of the Kif7 motor protein that can regulate the activator as well as the repressor functions of the Gli transcription factors in Shh signalling. This new work investigates how a null mutation in the Kif7 gene affects the formation of corticofugal and thalamocortical axon tracts and the migration of cortical interneurons. It demonstrates that Kif7 null mutant embryos present with ventriculomegaly and heterotopias as observed in patients carrying KIF7 mutations. The Kif7 mutation also disrupts the connectivity between cortex and thalamus and leads to an abnormal projection of thalamocortical axons. Moreover, cortical interneurons show migratory defects that are mirrored in cortical slices treated with the Shh inhibitor cyclopamine suggesting that the Kif7 mutation results in a down-regulation of Shh signalling. Interestingly, these defects are much less severe at later stages of corticogenesis.

Strengths/weaknesses:

The findings of this manuscript are clearly presented and are based on detailed analyses. Using a compelling set of experiments, especially the live imaging to monitor interneuron migration, the authors convincingly investigate Kif7's roles and their results support their major claims. The migratory defects in interneurons and the potential role of Shh signalling present novel findings and provide some mechanistic insights but rescue experiments would further support Kif7's role in interneuron migration. Similarly, the mechanism underlying the misprojection which has previously been reported in other cilia mutants remains unexplored. Taken together, this manuscript makes novel contributions to our understanding of the role of primary cilia in forebrain development and to the aetiology of the neural symptons in ciliopathy patients.

Comments on revisions:

The authors addressed most of the points I raised in my original review. However, I am not convinced by the figures the authors present on Shh protein expression. The "bright tiny dots" of Shh protein in the cortex are not visible on the images in Figure 7. I wonder whether the authors could present higher magnification and/or black and white images with increased contrast.

Reviewer #2 (Public review):

Summary:

This study investigates the role of KIF7, a ciliary kinesin involved in the Sonic Hedgehog (SHH) signaling pathway, in cortical development using Kif7 knockout mice. The researchers examined embryonic cortex development (mainly at E14.5), focusing on structural changes and neuronal migration abnormalities.

Strengths:

(1) The phenotype observed is interesting, and the findings provide neurodevelopmental insight into some of the symptoms and malformations seen in patients with KIF7 mutations.

(2) The authors assess several features of cortical development, including structural changes in layers of the developing cortex, connectivity of the cortex with thalamus, as well as migration of cINs from CGE and MGE to cortex.

Weaknesses:

(1) The Kif7 null does have phenotype differences from individual mutations seen in patients. It would be interesting to add more thoughts about how the null differs from these mutants in ciliary structure and SHH signaling via the cilium.

(2) The description of altered cortex development at E14.5 is perhaps rather descriptive. It would be useful to assess more closely the changes occurring in different cell types and stages. For this it seems very important to have a time course of cortical development and how the structural organization changes over time. This would be easy to assess with the addition of serial sections from the same mice. It might also be interesting to see how SHH signaling is altered in different cortical cell types over time with a SHH signaling reporter mouse.

(3) Abnormal neurodevelopmental phenotypes have been widely reported in the absence of other key genes affecting primary cilia function (Willaredt et al., J Neurosci 2008; Guo et al., Nat Commun 2015). It would be interesting to have more discussion of how the Kif7 null phenotype compares to some of these other mutants.

(4) The authors see alterations in cIN migration to the cortex and observe distinct differences in the pattern of expression of Cxcl12 as well as suggest cell intrinsic differences within cIN in their ability to migrate. The slice culture experiments though make it a little difficult to interpret the cell intrinsic effects on cIN of loss of Kif7, as the differences in Cxcl12 patterns still exist presumably in the slice cultures. It would be useful to assess their motility in an assay where they were isolated, as well as assess transcriptional changes in cINs in vivo lacking KIF7 for expression patterns that may affect motility or other aspects of migration.

Comments on revisions:

The authors have made significant and thoughtful responses as well as experimental additions to the authors comments. Their efforts are appreciated and the manuscript is much improved.

Reviewer #1 (Public review):

Summary:

In this study by Li et al., the authors re-investigated the role of cDC1 for atherosclerosis progression using the ApoE model. First, the authors confirmed the accumulation of cDC1 in atherosclerotic lesions in mice and humans. Then, in order to examine the functional relevance of this cell type, the authors developed a new mouse model to selectively target cDC1. Specifically, they inserted the Cre recombinase directly after the start codon of the endogenous XCR1 gene, thereby avoiding off-target activity. Following validation of this model, the authors crossed it with ApoE-deficient mice and found a striking reduction of aortic lesions (numbers and size) following a high-fat diet. The authors further characterized the impact of cDC1 depletion on lesional T cells and their activation state. Also, they provide in-depth transcriptomic analyses of lesional in comparison to splenic and nodal cDC1. These results imply cellular interactions between lesion T cells and cDC1. Finally, the authors show that the chemokine XCL1, which is produced by activated CD8 T cells (and NK cells), plays a key role in the interaction with XCR1-expressing cDC1 and particularly in the atherosclerotic disease progression.

Strengths:

The surprising results on XCL1 represent a very important gain in knowledge. The role of cDC1 is clarified with a new genetic mouse model.

Weaknesses:

My criticism is limited to the analysis of the scRNAseq data of the cDC1. I think it would be important to match these data with published data sets on cDC1. In particular, the data set by Sophie Janssen's group on splenic cDC1 might be helpful here (PMID: 37172103; https://www.single-cell.be/spleen_cDC_homeostatic_maturation/datasets/cdc1). It would be good to assign a cluster based on the categories used there (early/late, immature/mature, at least for splenic DC).

Reviewer #2 (Public review):

This study investigates the role of cDC1 in atherosclerosis progression using Xcr1Cre-Gfp Rosa26LSL-DTA ApoE-/- mice. The authors demonstrate that selective depletion of cDC1 reduces atherosclerotic lesions in hyperlipidemic mice. While cDC1 depletion did not alter macrophage populations, it suppressed T cell activation (both CD4+ and CD8+ subsets) within aortic plaques. Further, targeting the chemokine Xcl1 (ligand of Xcr1) effectively inhibits atherosclerosis. The manuscript is well-written, and the data are clearly presented. However, several points require clarification:

(1) In Figure 1C (upper plot), it is not clear what the Xcr1 single-positive region in the aortic root represents, or whether this is caused by unspecific staining. So I wonder whether Xcr1 single-positive staining can reliably represent cDC1. For accurate cDC1 gating in Figure 1E, Xcr1+CD11c+ co-staining should be used instead.

(2) Figure 4D suggests that cDC1 depletion does not affect CD4+/CD8+ T cells. However, only the proportion of these subsets within total T cells is shown. To fully interpret effects, the authors should provide:<br /> a) Absolute numbers of total T cells in aortas.<br /> b) Absolute counts of CD4+ and CD8+ T cells.

(3) How does T cell activation mechanistically influence atherosclerosis progression? Why was CD69 selected as the sole activation marker? Were other markers (e.g., KLRG1, ICOS, CD44) examined to confirm activation status?

(4) Figure 7B: Beyond cDC1/2 proportions within cDCs, please report absolute counts of: Total cDCs,cDC1, and cDC2 subsets. Figure 7D: In addition to CD4+/CD8+ T cell proportions, the following should be included:<br /> a) Total T cell numbers in aortas<br /> b) Absolute counts of CD4+ and CD8+ T cells.

(5) cDC1 depletion reduced CD69+CD4+ and CD69+CD8+ T cells, whereas Xcl1 depletion decreased Xcr1+ cDC1 cells without altering activated T cells. How do the authors explain these different results? This discrepancy needs explanation.

Reviewer #1 (Public review):

Summary:

This study presents a compelling strategy for ultrasound-mediated immunomodulation in HCC, supported by robust scRNA-Seq data. While the mechanistic depth and translational validation require further refinement, the work significantly advances the field of noninvasive cancer immunotherapy. Addressing the major concerns, particularly regarding calcium signaling specificity and STNDs@Ca²⁺ safety, will strengthen the manuscript's impact.

Strengths:

(1) Innovative Approach:<br /> The integration of noninvasive ultrasound with calcium-targeted nanotechnology (STNDs@Ca²⁺) represents a significant advancement in cancer immunotherapy. The dual mechanism - direct immunomodulation via FUS and calcium delivery via nanoparticles - is both novel and promising.

(2) Comprehensive Mechanistic Insights:<br /> The use of scRNA-seq and flow cytometry provides a detailed map of immune cell dynamics, highlighting key pathways (TNF, NFκB, MAPK) and cellular transitions (e.g., MDSC suppression, CD8⁺/NK cell activation).

(3) Robust Preclinical Validation:<br /> The study validates findings in two distinct HCC models (H22 and Hepal-6), demonstrating consistent tumor suppression (>70-90%) and prolonged survival, which strengthens translational relevance.

Weaknesses:

Major Issues:

(1) Mechanistic Specificity of Calcium Influx:<br /> While the study attributes immunomodulation to ultrasound-induced calcium influx, the exact mechanism (e.g., involvement of mechanosensitive channels like Piezo1 or TRP families) remains underdiscussed. The qRT-PCR data shows no changes in TRP channels, but the upregulation of Piezo1 warrants deeper exploration.

Suggestion: The authors should include experiments to inhibit Piezo1 or other calcium channels to confirm their role in FUS-mediated effects.

(2) STNDs@Ca²⁺ Biodistribution and Safety:<br /> Although biodistribution data show splenic accumulation, potential off-target effects (e.g., liver/lung uptake) and long-term toxicity are not fully addressed. The serum biochemical analysis (Table 2) lacks critical markers like inflammatory cytokines or immune cell counts.

Suggestion: The authors should provide longitudinal toxicity data (e.g., histopathology beyond 3 hours) and assess systemic immune activation/inflammation.

(3) Statistical and Technical Clarifications:<br /> The statistical methods for multi-group comparisons (e.g., ANOVA vs. t-test) are inconsistently described. For instance, Figure 1 labels significance without specifying correction for multiple comparisons.

Suggestion: the authors should clarify statistical methods in figure legends and the Methods section; apply Bonferroni or FDR correction where appropriate.

(4) Interpretation of scRNA-seq Data:<br /> The clustering of MDSCs using surface markers (Itgam/Ly6c2/Ly6g) overlaps with conventional myeloid populations (Supplementary Figure 16), raising questions about subset specificity.

Suggestion: The authors should validate MDSC identity using functional assays (e.g., T cell suppression) or additional markers (e.g., Arg1, iNOS).

Reviewer #2 (Public review):

Summary:

Wang et al. studied the therapeutic potential of focused ultrasound noninvasively stimulating the spleen (FUS sti. spleen) to modulate the splenic immunity, with an aim to exert anti-tumor effect. They found that the treatment enhanced antitumor capability in CD8 T/ NK cells and reduced the immunosuppression, facilitating the inhibition of HCC growth in vivo.

Strengths:

They have utilized bulk RNA sequencing, single cell RNA sequencing, and flow cytometry to investigate the immune and tumor cell profiling in the mouse models upon FUS sti. spleen. Moreover, they highlighted the importance of combining FUS with spleen-targeted nanodroplets encapsulating bioavailable calcium ions (STND@Ca2+), which facilitated the calcium influx into the murine spleen and further enhanced the therapeutic efficacy of FUS.

Weaknesses:

While the study is interesting and potentially clinically impactful, the mechanism of action of the therapy is not fully elucidated. It would benefit from more rigorous approaches. With the theoretical part strengthened, this paper would be of interest to cancer cell biologists and clinician scientists working on the oncology field.

Reviewer #1 (Public review):

Summary:

In this fMRI study, the authors wished to assess neural mechanisms supporting flexible "temporal construals". For this, human participants learned a story consisting of fifteen events. During fMRI, events were shown to them and they were instructed to consider the event from "an internal" or from "an external" perspective. The authors found opposite patterns of brain activity in the posterior parietal cortex and the anterior hippocampus for the internal and the external viewpoint. They conclude that allocentric sequences are stored in the hippocampus, whereas egocentric sequences are used in the parietal cortex. The claims align with previous fMRI work addressing this question.

Strengths:

The research topic is fascinating, and very few labs in the world are asking the question of how time is represented in the human brain. Working hypotheses have been recently formulated, and this work seems to want to tackle some of them.

Weaknesses:

The current writing is fuzzy both conceptually and experimentally. I cannot provide a sufficiently well-informed assessment of the quality of the experimental work because there is a paucity of details provided in the report. Any future revisions will likely improve transparency.

(1) Improving writing and presentation:

The abstract and the introduction make use of loaded terms such as "construals", "mental timeline", "panoramic views" in very metaphoric and unexplained ways. The authors do not provide a comprehensive and scholarly overview of these terms, which results in verbiage and keywords/name-dropping without a clear general framework being presented. Some of these terms are not metaphors. They do refer to computational concepts that the authors should didactically explain to their readership. This is all the more important that some statements in the Introduction are misattributed or factually incorrect; some statements lack attributions (uncited published work).

Once the theory, the question, and the working hypothesis are clarified, the authors should carefully explain the task.

(2) The experimental approach lacks sufficient details to be comprehensible to a general audience. In my opinion, the results are thus currently uninterpretable. I highlight only a couple of specific points (out of many). I recommend revision and clarification.

a) No explanation of the narrative is being provided. The authors report a distribution of durations with no clear description of the actual sequence of events. The authors should provide the text that was used, how they controlled for low-level and high-level linguistic confounds.

b) The authors state, "we randomly assigned 15 phrases to the events twice". It is impossible to comprehend what this means. Were these considered stimuli? Controls? IT is also not clear which event or stimulus is part of the "learning set" and whether these were indicated to be such to participants.

c) The left/right counterbalancing is not being clearly explained. The authors state that there is counterbalancing, but do not sufficiently explain what it means concretely in the experiment. If a weak correlation exists between sequential position and distance, it also means that the position and the distance have not been equated within. How do the authors control for these?

d) The authors used two tasks. In the "external perspective" one, the authors asked participants to report whether events were part of the same or a different part of the day. In the "internal perspective one", the authors asked participants to project themselves to the reference event and to determine whether the target event occurred before or after the projected viewpoint. The first task is a same/different recognition task. The second task is a temporal order task (e.g., Arzy et al. 2009). These two asks are radically different and do not require the same operationalization. The authors should minimally provide a comprehensive comparison of task requirements, their operationalization, and, more importantly, assess the behavioral biases inherent to each of these tasks that may confound brain activity observed with fMRI.

e) The authors systematically report interpreted results, not factual data. For instance, while not showing the results on behavioral outcomes, the authors directly interpret them as symbolic distance effects.

Crucially, the authors do not comment on the obvious differences in task difficulty in these two tasks, which demonstrates a substantial lack of control in the experimental design. The same/different task (task 1 called "external perspective") comes with known biases in psychophysics that are not present in the temporal order task (task 2 called " internal perspective"). The authors also did not discuss or try to match the performance level in these two tasks. Accordingly, the authors claim that participants had greater accuracy in the external (same/different) task than in the internal task, although no data are shown and provided to support this report. Further, the behavioral effect is trivialized by the report of a performance accuracy trade-off that further illustrates that there is a difference in the task requirements, preventing accurate comparison of the two tasks.

All fMRI contrasts are also confounded by this experimental shortcoming, seeing as they are all reported at the interaction level across a task. For instance, in Figure 4, the authors report a significant beta difference between internal and external tasks. It is impossible to disentangle whether this effect is simply due to task difference or to an actual processing of the duration that differs across tasks, or to the nature of the representation (the most difficult to tackle, and the one chosen by the authors).

Conclusion:

In conclusion, the current experimental work is confounded and lacks controls. Any behavioral or fMRI contrasts between the two proposed tasks can be parsimoniously accounted for by difficulty or attentional differences, not the claim of representational differences being argued for here.

Reviewer #2 (Public review):

Summary:

Xu et al. used fMRI to examine the neural correlates associated with retrieving temporal information from an external compared to internal perspective ('mental time watching' vs. 'mental time travel'). Participants first learned a fictional religious ritual composed of 15 sequential events of varying durations. They were then scanned while they either (1) judged whether a target event happened in the same part of the day as a reference event (external condition); or (2) imagined themselves carrying out the reference event and judged whether the target event occurred in the past or will occur in the future (internal condition). Behavioural data suggested that the perspective manipulation was successful: RT was positively correlated with sequential distance in the external perspective task, while a negative correlation was observed between RT and sequential distance for the internal perspective task. Neurally, the two tasks activated different regions, with the external task associated with greater activity in the supplementary motor area and supramarginal gyrus, and the internal condition with greater activity in default mode network regions. Of particular interest, only a cluster in the posterior parietal cortex demonstrated a significant interaction between perspective and sequential distance, with increased activity in this region for longer sequential distances in the external task, but increased activity for shorter sequential distances in the internal task. Only a main effect of sequential distance was observed in the hippocampus head, with activity being positively correlated with sequential distance in both tasks. No regions exhibited a significant interaction between perspective and duration, although there was a main effect of duration in the hippocampus body with greater activity for longer durations, which appeared to be driven by the internal perspective condition. On the basis of these findings, the authors suggest that the hippocampus may represent event sequences allocentrically, whereas the posterior parietal cortex may process event sequences egocentrically.

Strengths:

The topic of egocentric vs. allocentric processing has been relatively under-investigated with respect to time, having traditionally been studied in the domain of space. As such, the current study is timely and has the potential to be important for our understanding of how time is represented in the brain in the service of memory. The study is well thought out, and the behavioural paradigm is, in my opinion, a creative approach to tackling the authors' research question. A particular strength is the implementation of an imagination phase for the participants while learning the fictional religious ritual. This moves the paradigm beyond semantic/schema learning and is probably the best approach besides asking the participants to arduously enact and learn the different events with their exact timings in person. Importantly, the behavioural data point towards successful manipulation of internal vs. external perspective in participants, which is critical for the interpretation of the fMRI data. The use of syllable length as a sanity check for RT analyses, as well as neuroimaging analyses, is also much appreciated.

Weaknesses/Suggestions:

Although the design and analysis choices are generally solid, there are a few finer details/nuances that merit further clarification or consideration in order to strengthen the readers' confidence in the authors' interpretation of their data.

(1) Given the known behavioural and neural effects of boundaries in sequence memory, I was wondering whether the number of traversed context boundaries (i.e., between morning-afternoon, and afternoon-evening) was controlled for across sequential length in the internal perspective condition? Or, was it the case that reference-target event pairs with higher sequential numbers were more likely to span across two parts of the day compared to lower sequential numbers? Similarly, did the authors examine any potential differences, whether behaviourally or neurally, for day part same vs. day part different external task trials?

(2) I would appreciate further insight into the authors' decision to model their task trials as stick functions with duration 0 in their GLMs, as opposed to boxcar functions with varying durations, given the potential benefits of the latter (e.g., Grinband et al., 2008). I concur that in certain paradigms, RT is considered a potential confound and is taken into account as a nuisance covariate (as the authors have done here). However, given that RTs appear to be critical to the authors' interpretation of participant behavioural performance, it would imply that variations in RT actually reflect variations in cognitive processes of interest, and hence, it may be worth modelling trials as boxcar functions with varying durations.

(3) The activity pattern across tasks and sequential distance in the posterior parietal cortex appears to parallel the RT data. Have the authors examined potential relationships between the two (e.g., individual participant slopes for RT across sequential distance vs. activity betas in the posterior parietal cortex)?

(4) There were a few places in the manuscript where the writing/discussion of the wider literature could perhaps be tightened or expanded. For instance:

i) On page 16, the authors state 'The negative correlation between the activation level in the right PPC and sequential distance has already been observed in a previous fMRI study (Gauthier & van Wassenhove, 2016b). The authors found a similar region (the reported MNI coordinate of the peak voxel was 42, -70, 40, and the MNI coordinate of the peak voxel in the present study was 39, -70, 35), of which the activation level went up when the target event got closer to the self-positioned event. This finding aligns with the evidence suggesting that the posterior parietal cortex implements egocentric representations.' Without providing a little more detail here about the Gauthier & van Wassenhove study and what participants were required to do (i.e., mentally position themselves at a temporal location and make 'occurred before' vs. 'occurred after' judgements of a target event), it could be a little tricky for readers to follow why this convergence in finding supports a role for the posterior parietal cortex in egocentric representations.

ii) Although the authors discuss the Lee et al. (2020) review and related studies with respect to retrospective memory, it is critical to note that this work has also often used prospective paradigms, pointing towards sequential processing being the critical determinant of hippocampal involvement, rather than the distinction between retrospective vs. prospective processing.

iii) The authors make an interesting suggestion with respect to hippocampal longitudinal differences in the representation of event sequences, and may wish to relate this to Montagrin et al. (2024), who make an argument for the representation of distant goals in the anterior hippocampus and immediate goals in the posterior hippocampus.

Reviewer #1 (Public review):

Summary:

This is a rigorous data-driven modeling study, extending the authors' previous model of spinal locomotor central pattern generator (CPG) circuits developed for the mouse spinal cord and adapted here to the rat to explore potential circuit-level changes underlying altered speed-dependent gaits, due to asymmetric (lateral) thoracic spinal hemisection and symmetric midline contusion. The model reproduces key features of the rat speed-dependent gait-related experimental data before injury and after recovery from these two different thoracic spinal cord injuries and suggests injury-specific mechanisms of circuit reorganization underlying functional recovery. There is much interest in the mechanisms of locomotor behavior recovery after spinal cord injury, and data-driven behaviorally relevant circuit modeling is an important approach. This study represents an important advance in the authors' previous experimental and modeling work on locomotor circuitry and in the motor control field.

Strengths:

(1) The authors use an advanced computational model of spinal locomotor circuitry to investigate potential reorganization of neural connectivity underlying locomotor control following recovery from symmetrical midline thoracic contusion and asymmetrical (lateral) hemisection injuries, based on an extensive dataset for the rat model of spinal cord injury.

(2) The rat dataset used is from an in vivo experimental paradigm involving challenging animals to perform overground locomotion across the full range of speeds before and after the two distinct spinal cord injury models, enabling the authors to more completely reveal injury-specific deficits in speed-dependent interlimb coordination and locomotor gaits.

(3) The model reproduces the rat gait-related experimental data before injury and after recovery from these two different thoracic spinal cord injuries, which exhibit roughly comparable functional recovery, and suggests injury-specific, compensatory mechanisms of circuit reorganization underlying recovery.

(4) The model simulations suggest that recovery after lateral hemisection mechanistically involves partial functional restoration of descending drive and long propriospinal pathways. In contrast, recovery following midline contusion relies on reorganization of sublesional lumbar circuitry combined with altered descending control of cervical networks.

(5) These observations suggest that symmetrical (contusion) and asymmetrical (lateral hemisection) injuries induce distinct types of plasticity in different spinal cord regions, suggesting that injury symmetry partly dictates the location and type of neural plasticity supporting recovery.

(6) The authors suggest that therapeutic strategies may be more effective by targeting specific circuits according to injury symmetry.

Weaknesses:

The recovery mechanisms implemented in the model involve circuit connectivity/connection weights adjustment based on assumptions about the structures involved and compensatory responses to the injury. As the authors acknowledge, other factors affecting locomotor patterns and compensation, such as somatosensory afferent feedback, neurochemical modulator influences, and limb/body biomechanics, are not considered in the model.

Reviewer #2 (Public review):

Summary:

In this paper, the authors present a detailed computational model and experimental data concerning overground locomotion in rats before and after recovery from spinal cord injury. They are able to manually tune the parameters of this physiologically based, detailed model to reproduce many aspects of the observed animals' locomotion in the naive case and in two distinct injury cases.

Strengths:

The strengths are that the model is driven to closely match clean experimental data, and the model itself has detailed correspondence to proposed anatomical reality. As such, this makes the model more readily applicable to future experimental work. It can make useful suggestions. The model reproduces a large number of conditions across frequencies, and with the model structure changed by injury and recovery. The model is extensive and is driven by known structures, with links to genetic identities, and has been extensively validated across multiple experiments and manipulations over the years. It models a system of critical importance to the field, and the tight coupling to experimental data is a real strength.

Weaknesses:

A downside is that, scientifically, here, the only question tackled is one of sufficiency. By manually tuning parameters in a manner that aligns with the field's understanding from experimental work, the detailed model can accurately reproduce the experimental findings. One of the benefits of computational models is that the counterfactual can be tested to provide evidence against alternative hypotheses. That isn't really done here. I'm fairly certain that there are competing theories regarding what happens during recovery from a hemi-section injury and a contusion injury. The model could be used to make predictions for some alternative hypotheses, supporting or rejecting theories of recovery. This may be part of future plans. Here, the focus is on showing that the model is capable of reproducing the experimental results at all, for any set of parameters, however tuned.

Reviewer #3 (Public review):

Summary:

This study describes a computational model of the rat spinal locomotor circuit and how it could be reconfigured after lateral hemisection or contusion injuries to replicate gaits observed experimentally.

The model suggests the emergence of detour circuits after lateral hemisection, whereas after a midline contusion, the model suggests plasticity of left-right and sensory inputs below the injury.

Strengths:

The model accurately models many known connections within and between forelimb and hindlimb spinal locomotor circuits.

The simulation results mirror closely gait parameters observed experimentally. Many gait parameters were studied, as well as variability in these parameters in intact versus injured conditions.

Weaknesses:

The study could provide some sense of the relative importance of the various modified connectivities after injury in setting the changes in gait seen after the two types of injuries.

Overall, the authors achieved their aims, and the model provides solid support for the changes in connectivity after the two types of injuries were modelled. This work emphasizes specific changes in connectivity after lateral hemisection or after contusion that could be investigated experimentally. The model is available for public use and could serve as a tool to analyze the relative importance of various highlighted or previously undiscovered changes in connectivity that may underlie the recovery of locomotor function in spinalized rats.

Reviewer #1 (Public review):

In this manuscript, Saeb et al reported the mechanistic roles of the flexible stalk domain in sTREM2 function using molecular dynamics simulations. They have reported some interesting molecular bases explaining why sTREM2 shows protective effects during AD, such as partial extracellular stalk domain promoting binding preference and stabilities of sTREM2 with its ligand even in the presence of known AD-risk mutation, R47H. Furthermore, they found that the stalk domain itself acts as the site for ligand binding by providing an "expanded surface", known as 'Expanded Surface 2' together with the Ig-like domain. Also, they observed no difference in the binding free energy of phosphatidyl-serine with wild TREM2-Ig and mutant TREM2-Ig, which is a bit inconsistent with the previous report with experiment studies by Journal of Biological Chemistry 293, (2018), Alzheimer's and Dementia 17, 475-488 (2021), Cell 160, 1061-1071 (2015).

Reviewer #2 (Public review):

Significance:

TREM2 is an immunomodulatory receptor expressed on myeloid cells and microglia in the brain. TREM2 consists of a single immunoglobular (Ig) domain that leads into a flexible stalk, transmembrane helix, and short cytoplasmic tail. Extracellular proteases can cleave TREM2 in its stalk and produce a soluble TREM2 (sTREM2). TREM2 is genetically linked to Alzheimer's disease (AD), with the strongest association coming from an R47H variant in the Ig domain. Despite intense interest, the full TREM2 ligand repertoire remains elusive, and it is unclear what function sTREM2 may play in the brain. The central goal of this paper is to assess the ligand-binding role of the flexible stalk that is generated during the shedding of TREM2. To do this, the authors simulate the behavior of constructs with and without stalk. However, it is not clear why the authors chose to use the isolated Ig domain as a surrogate for full-length TREM2. Additionally, experimental binding evidence that is misrepresented by the authors contradicts the proposed role of the stalk.

Summary and strengths:

The authors carry out MD simulations of WT and R47H TREM2 with and without the flexible stalk. Simulations are carried out for apo TREM2 and for TREM2 in complex with various lipids. They compare results using just the Ig domain to results including the flexible stalk that is retained following cleavage to generate sTREM2. The computational methods are well-described and should be reproducible. The long simulations are a strength, as exemplified in Figure 2A where a CDR2 transition happens at ~400-600 ns. The stalk has not been resolved in structural studies, but the simulations suggest the intriguing and readily testable hypothesis that the stalk interacts with the Ig domain and thereby contributes to the stability of the Ig domain and to ligand binding. I suspect biochemists interested in TREM2 will make testing this hypothesis a high priority.

Comments on latest version:

The authors have addressed my critiques and carried out additional simulations, as requested. I would upgrade my assessment of the evidence to "solid."

Reviewer #1 (Public review):

Summary:

Odor- and taste-sensing are mediated by two different systems, the olfactory and gustatory systems, and have different behavioral roles. In this study, Wei et al. challenge this dichotomy by showing that odors can activate gustatory receptor neurons (GRNs) in Drosophila to promote feeding responses, including the proboscis extension response (PER) that was previously thought to be driven only by taste. While previous studies suggested that odors can promote PER to appetitive tastants, Wei et al. go further to show that odors alone cause PER, this effect is mediated through sweet-sensing GRNs, and sugar receptors are required. The study also shows that odor detection by bitter-sensing GRNs suppresses PER. The authors' conclusions are supported by behavioral assays, calcium imaging, electrophysiological recordings, and genetic manipulations. The observation that both attractive and aversive odors promote PER leaves an open question as to why this effect is adaptive. Overall, the study sheds new light on chemosensation and multimodal integration by showing that odor and taste detection converge at the level of sensory neurons, a finding that is interesting and surprising while also being supported by another recent study (Dweck & Carlson, Sci Advances 2023).

Strengths:

(1) The main finding that odors alone can promote PER by activating sweet-sensing GRNs is interesting and novel.

(2) The study uses video tracking of the proboscis to quantify PER rather than manual scoring, which is typically used in the field. The tracking method is less subjective and provides a higher-resolution readout of the behavior.

(3) The study uses calcium imaging and electrophysiology to show that odors activate GRNs. These represent complementary techniques that measure activity at different parts of the GRN (axons versus dendrites, respectively) and strengthen the evidence for this conclusion.

(4) Genetic manipulations show that odor-evoked PER is primarily driven by sugar GRNs and sugar receptors rather than olfactory neurons. This is a major finding that distinguishes this work from previous studies of odor effects on PER and feeding (e.g., Reisenman & Scott, 2019; Shiraiwa, 2008) that assumed or demonstrated that odors were acting through olfactory neurons.

Weaknesses/Limitations:

(1) Many of the odor effects on behavior or neuronal responses were only observed at very high concentrations. Most effects seemed to require concentrations of at least 10^-2 (0.01 v/v), which is at the high end of the concentration range used in olfactory studies (e.g., Hallem et al., 2004), and most experiments in the paper used a far higher concentration of 0.5 v/v. It is unclear whether these are concentrations that would be naturally encountered by flies. In addition, it is difficult to compare the concentrations used for electrophysiology and behavior given that they are presented in solution versus volatile form.

(2) The timecourse of GRN activation by odors seems quite prolonged (and possibly delayed, depending on the exact timing of odor onset to the fly), and this timecourse is not directly compared with activation by tastes to determine whether it is a property of the calcium sensor or a real difference.

(3) While the overall effect of different conditions is tested using appropriate statistical methods, post-hoc tests are not always used to determine which specific groups are different from each other (e.g., which odors and concentrations elicit significant PER compared to air or mineral oil controls in Fig. 1; which odors show impaired responses without olfactory organs in Fig. 2A).

Discrepancies with previous studies:

These discrepancies are important to note but should not necessarily be considered "weaknesses" of the present study.

(1) It is not entirely clear why PER to odors alone has not been previously reported, especially as this study shows that it is a broad effect evoked by many different odors. Previous studies (Oh et al., 2021; Reisenman & Scott, 2019; Shiraiwa, 2008) tested the effect of odors on PER and only observed enhancement of PER to sugar rather than odor-evoked PER; some of these studies explicitly show no effect of odor alone or odor with low sugar concentration. In the Response to Reviewers, the authors propose that genetic background may explain discrepancies, but this is not discussed much in the paper itself. Differences in behavioral quantification (automated vs. manual scoring, quantification of PER duration versus probability) may also contribute.

(2) The calcium imaging data showing that sugar GRNs respond to a broad set of odors contrasts with results from Dweck & Carlson (Sci Adv, 2023) who recorded sugar neurons with electrophysiology and observed responses to organic acids, but not other odors. This discrepancy is mentioned in the Discussion but the underlying reason is not clear.

Reviewer #3 (Public review):

Summary:

Using flies, Kazama et al. combined behavioral analysis, electrophysiological recordings, and calcium imaging experiments to elucidate how odors activate gustatory receptor neurons (GRNs) and elicit a proboscis extension response, which is interpreted as a feeding response.

The authors used DeepLabCut v2.0 to estimate the extension of the proboscis, which represents an unbiased and more precise method for describing this behavior compared to manual scoring.

They demonstrated that the probability of eliciting a proboscis extension increases with higher odor concentrations. The most robust response occurs at a 0.5 v/v concentration, which, despite being diluted in the air stream, remains a relatively high concentration. Although the probability of response is not particularly high it is higher than control stimuli. Notably, flies respond with a proboscis extension to both odors that are considered positive and those regarded as negative.

The authors used various transgenic lines to show that the response is mediated by GRNs. Specifically, inhibiting Gr5a reduces the response, while inhibiting Gr66a increases it in fed flies. Additionally, they find that odors induce a strong positive response in both types of GRNs, which is abolished when the labella of the proboscis are covered. This response was also confirmed through electrophysiological tip recordings.

Finally, the authors demonstrated that the response increases when two stimuli of different modalities, such as sucrose and odors, are presented together, suggesting clear multimodal integration

Strengths:

The integration of various techniques, which collectively supports the robustness of the results.<br /> The assessment of electrophysiological recordings in intact animals, preserving natural physiological conditions.

Weaknesses:

Only highly concentrated odours are capable of evoking positive responses and, even then, the proportion remains relatively low.

The authors have incorporated my suggestions.

Reviewer #1 (Public review):

Summary:

The authors attempted to dissect the function of a long non-coding RNA, lnc-FANCI-2, in cervical cancer. They profiled lnc-FANCI-2 in different cell lines and tissues, generated knockout cell lines, and characterized the gene using multiple assays.

Strengths:

A large body of experimental data has been presented and can serve as a useful resource for the scientific community, including transcriptomics and proteomics datasets. The reported results also span different parts of the regulatory network and open up multiple avenues for future research.

Weaknesses:

The write-up is somewhat unfocused and lacks deep mechanistic insights in some places.

Comments on revisions:

The manuscript is much improved. I am satisfied with the authors' responses.

Reviewer #3 (Public review):

Summary:

A long noncoding RNA, lnc-FANCI-2, was reported to be regulated by HPV E7 oncoprotein and a cell transcription factor, YY1 by this group. The current study focuses on the function of lnc-FANCI-2 in HPV-16 positive cervical cancer is to intrinsically regulate RAS signaling, thereby facilitating our further understanding additional cellular alterations during HPV oncogenesis. Authors used the advanced technical approaches such as KO, transcriptome and (IRPCRP) and LC- MS/MS analyses in the current study and concluded that KO Inc-FANCI-2 significantly increase RAS signaling, especially phosphorylation of Akt and Erk1/2.

Strengths:

(1) HPV E6E7 are required for full immortalization and maintenance of malignant phenotype of cervical cancer, but they are NOT sufficient for full transformation and tumorigenesis. This study helps further the understanding of other cellular alterations in HPV oncogenesis.<br /> (2) lnc-FANCI-2 is upregulated in cervical lesion progression from CIN1, CIN2-3 to cervical cancer, cancer cell lines and HPV transduced cell lines.<br /> (3) Viral E7 of high-risk HPVs and host transcription factor YY1 are two major factors promoting lnc-FANCI-2 expression.<br /> (4) Proteomic profiling of cytosolic and secreted proteins showed inhibition of MCAM, PODXL2 and ECM1 and increased levels of ADAM8 and TIMP2 in KO cells.<br /> (5) RNA-seq analyses revealed that KO cells exhibited significantly increased RAS signaling but decreased IFN pathways.<br /> (6) Increased phosphorylated Akt and Erk1/2, IGFBP3, MCAM, VIM, and CCND2 (cyclin D2) and decreased RAC3 were observed in KO cells.

Comments on revisions:

The revised manuscript has been significantly improved. The authors addressed all my concerns.

Reviewer #1 (Public review):

Summary

Detecting unexpected epistatic interactions between multiple mutations requires a robust null expectation-or neutral function-that predicts the combined effects of multiple mutations on phenotype based on the individual effects of single mutations. This study evaluated the relevance of the product neutrality function, where double-mutant fitness is represented as a multiplicative combination of single-mutant fitness in the absence of epistatic interactions. The authors used a recent large dataset on fitness, specifically yeast colony size, to analyze epistatic interactions.

The study confirmed that the product function outperformed other neutral functions in predicting double-mutant fitness, showing no bias between negative and positive epistatic interactions. Additionally, in the theoretical portion of the study, the authors employed a previously established theoretical model of bacterial cell growth to simulate growth rates of both single- and double-mutants under multiple parameters. The simulations similarly demonstrated that the product function was superior to other functions in predicting the fitness of hypothetical double-mutants. Based on these findings, the authors concluded that the product function is a robust tool for analyzing epistatic interactions in growth fitness and effectively reflects how growth rates depend on the combination of multiple biochemical pathways.

Strength

By leveraging a previously published large dataset of yeast colony sizes for single- and double-knockout mutants, this study validated the relevance of the product function, which has frequently been used in genetics to analyze epistatic interactions. The confirmation that the product function provides a more reliable prediction of double-mutant fitness compared to other neutral functions is valuable for researchers analyzing epistatic interactions, particularly those working with the same dataset.<br /> Notably, this dataset has been previously used in studies exploring epistatic interactions with the product neutrality function. This study's findings affirm the validity of using the product function, which could enhance confidence in the conclusions drawn by those earlier studies. Consequently, both researchers utilizing this dataset and readers of prior research will benefit from the confirmation provided by this study.

Weakness

This study contains several serious problems, primarily stemming from the following issues: ignoring the substantial differences in the mechanisms regulating cell growth between prokaryotes and eukaryotes and adopting an overly specific and unrealistic set of assumptions in the mutation model. Below, the details are discussed.

(1) Misapplication of prokaryotic growth models

The mechanistic origin of the multiplicative model observed in yeast colony fitness is explained using a bacterial cell growth model. However, there is no valid justification for linking these two systems. The bacterial growth model, the Scott-Hwa model, heavily rely on specific molecular mechanisms, such as ppGpp-mediated regulation, which adjusts ribosome expression and activity during translation. In particular, this mechanism is critical to ensure growth-dependency of the fraction of ribosome in proteome in the Scott-Hwa model [https://doi.org/10.1111/j.1462-2920.2010.02357.x; https://doi.org/10.1073/pnas.2201585119]. Yeast cells lack this regulatory mechanism, making it inappropriate to directly apply bacterial growth models to yeast.<br /> The Weiße model is based on a larger set of underlying equations and involves more parameters than the Scott-Hwa model. In the original paper by Weiße et al. (PNAS, 2015), however, the model parameters were fitted solely to experimental data from E. coli, and the model's applicability to yeast was never assessed. In summary, for neither the Scott-Hwa model nor the Weiße model has it been demonstrated that the entire model quantitatively fits experimental data from yeast. A positive correlation between growth rate and RNA/protein ratio, often observed in yeast, supports only a limited portion of either model, and does not constitute validation of the models as a whole.

(2) Overly specific assumptions in the theoretical model

The theoretical model assumes that two mutations affect only independent parameters of specific biochemical processes. However, this overly restrictive assumption weakens the model's validity in explaining the general occurrence of the multiplicative model in mutations. Furthermore, experimental evidence suggests limitations of this approach. For example, in most viable yeast deletion mutants with reduced growth rates, the expression of ribosomal proteins remained largely unchanged, contrary to the predictions of the Scott-Hwa model [https://doi.org/10.7554/eLife.28034]. This discrepancy highlights that the Scott-Hwa model and its derivatives cannot reliably explain mutants' growth rates based on current experimental evidence.

(3) Limited reliability of the mechanistic origin of the multiplicative model

The authors seem to regard growth-optimizing feedback as the mechanistic origin of the multiplicative model. However, the importance of growth-optimizing feedback in explaining product neutrality heavily depends on the very specific framework of the Scott-Hwa model. As I pointed out above, the Scott-Hwa model is a bacterial growth model that considers only a narrowly defined set of biochemical reactions. Using such a narrow model to explore the mechanistic origin of product neutrality observed on a genome-wide scale appears to be inappropriate. Arguments based on either the Scott-Hwa model or the Weiße model fail to account for the generality of product neutrality across diverse genetic perturbations. These models, in their current form, do not explain the broader patterns of product neutrality observed experimentally.

Reviewer #2 (Public review):

The paper deals with the important question of gene epistasis, focusing on asking what is the correct null model for which we should declare no epistasis.

In the first part, they use the Synthetic Genetic Array dataset to claim that the effects of a double mutation on growth rate is well predicted by the product of the individual effects (much more than e.g. the additive model). The second (main) part shows this is also the prediction of two simple, coarse-grained models for cell growth.

I find the topic interesting, the paper well written, and the approach innovative.

Comments on revisions:

The authors have adequately addressed the comments raised in the review below, and I find that the paper has improved.

Reviewer #1 (Public review):

Summary:

This work provides a new potential tool to manipulate Tregs function for therapeutic use. It focuses on the role of PGAM in Tregs differentiation and function. The authors, interrogating publicly available transcriptomic and proteomic data of human regulatory T cells and CD4 T cells, state that Tregs express higher levels of PGAM (at both message and protein levels) compared to CD4 T cells. They then inhibit PGAM by using a known inhibitor ECGC and show that this inhibition affects Tregs differentiation. This result was also observed when they used antisense oligonucleotides (ASOs) to knockdown PGAM1.

PGAM1 catalyzes the conversion of 3PG to 2PG in the glycolysis cascade. However, the authors focused their attention on the additional role of 3PG: acting as starting material for the de novo synthesis of serine.

They hypothesized that PGAM1 regulates Tregs differentiation by regulating the levels of 3PG that are available for de novo synthesis of serine, which has a negative impact on Tregs differentiation. Indeed, they tested whether the effect on Tregs differentiation observed by reducing PGAM1 levels was reverted by inhibiting the enzyme that catalyzes the synthesis of serine from 3PG.

The authors continued by testing whether both synthesized and exogenous serine affect Tregs differentiation and continued with in vivo experiments to examine the effects of dietary serine restriction on Tregs function.

In order to understand the mechanism by which serine impacts Tregs function, the authors assessed whether this depends on the contribution of serine to one-carbon metabolism and to DNA methylation.

The authors therefore propose that extracellular serine and serine whose synthesis is regulated by PGAM1 induce methylation of genes Tregs associated, downregulating their expression and overall impacting Tregs differentiation and suppressive functions.

Strengths:

The strength of this paper is the number of approaches taken by the authors to verify their hypothesis. Indeed, by using both pharmacological and genetic tools in in vitro and in vivo systems they identified a potential new metabolic regulation of Tregs differentiation and function.

Reviewer #2 (Public review):

Summary:

The authors have tried to determine the regulatory role of Phosphoglycerate mutate (PGAM), an enzyme involved in converting 3-phosphoglycerate to 2-phosphoglycerate in glycolysis, in differentiation and suppressive function of regulatory CD4 T cells through de novo serine synthesis. This is done by contributing one carbon metabolism and eventually epigenetic regulation of Treg differentiation.

Strengths:

The authors have rigorously used inhibitors and antisense RNA to verify the contribution of these pathways in Treg differentiation in-vitro. This has also been verified in an in-vivo murine model of autoimmune colitis. This has further clinical implications in autoimmune disorders and cancer.

[Editors' note: The authors addressed important comments by the reviewers.]

Reviewer #1 (Public review):

Summary:

The authors strived for an inventory of GPCRs and GPCR pathway component genes within the genomes of 23 choanoflagellates and other close relatives of metazoans.

Strengths:

The authors generated a solid phylogenetic overview of the GPCR superfamily in these species. Intriguingly, they discover novel GPCR families, novel assortments of domain combinations, and novel insights into the evolution of those groups within the Opisthokonta clade. A particular focus is laid on adhesion GPCRs, for which the authors discover many hitherto unknown subfamilies based on Hidden Markov Models of the 7TM domain sequences, which were also reflected by combinations of extracellular domains of the homologs. In addition, the authors provide bioinformatic evidence that aGPCRs of choanoflagellates also contain a GAIN domain, which is self-cleavable, thereby reflecting the most remarkable biochemical feat of aGPCRs.

Weaknesses:

The chosen classification scheme for aGPCRs may require reassessment and amendment by the authors in order to prevent confusion with previously issued classification attempts of this family.

Reviewer #2 (Public review):

Summary:

The authors set out to characterise the GPCR family in choanoflagellates (and other unicellular holozoans). GPCRs are the most abundant gene family in many animal genomes, playing crucial roles in a wide range of physiological processes. Although they are known to evolve rapidly, GPCRs are an ancient feature of eukaryotic biology. Identifying conserved elements across the animal-protist boundary is therefore a valuable goal, and the increasing availability of genomes from non-animal holozoans provides new opportunities to explore evolutionary patterns that were previously obscured by limited taxon sampling. This study presents a comprehensive re-examination of GPCRs in choanoflagellates, uncovering examples of differential gene retention and revealing the dynamic nature of the GPCR repertoire in this group. As GPCRs are typically involved in environmental sensing, understanding how these systems evolved may shed light on how our unicellular ancestors adapted their signalling networks in the transition to complex multicellularity.

Strengths:

The paper combines a broad taxonomic scope with the use of both established and recently developed tools (e.g., Foldseek, AlphaFold), enabling a deep and systematic exploration of GPCR diversity. Each family is carefully described, and the manuscript also functions as an up-to-date review of GPCR classification and evolution. Although similar attempts to understand GPCR evolution were made over the last decade, the authors build on this foundation by identifying new families and applying improved computational methods to better predict structure and function. Notably, the presence of Rhodopsin-like GPCRs in some choanoflagellates and ichthyosporeans is intriguing, even though they do not fall within known animal subfamilies. The computational framework presented here is broadly applicable, offering a blueprint for surveying GPCR diversity in other non-model eukaryotes (and even in animal lineages), potentially revealing novel families relevant to drug discovery or helping revise our understanding of GPCR evolution beyond model systems.

Weaknesses:

While the study contributes several interesting observations, it does not radically revise the evolutionary history of the GPCR family. However, in an era increasingly concerned with the reproducibility of scientific findings, this is arguably a strength rather than a weakness. It is encouraging to see that previously established patterns largely hold, and that with expanded sampling and improved methods, new insights can be gained, especially at the level of specific GPCR subfamilies. Then, no functional follow-ups are provided in the model system Salpingoeca rosetta, but I am sure functional work on GPCRs in choanoflagellates is set to reveal very interesting molecular adaptations in the future.

Reviewer #1 (Public review):

Summary:

This is an interesting theoretical study examining the viability of Virtual Circular Genome (VCG) model, a recently proposed scenario of prebiotic replication in which a relatively long sequence is stored as a collection of its shorter subsequences (and their compliments). It was previously pointed out that VCG model is prone to so-called sequence scrambling which limits the overall length of such a genome. In the present paper, additional limitations are identified. Specifically, it is shown that VCG is well replicated when the oligomers are elongated by sufficiently short chains from "feedstock" pool. However, ligation of oligomers from VCG itself results in a high error rate. I believe the research is of high quality and well written. However, the presentation could be improved and the key messages could be clarified.

(1) It is not clear from the paper whether the observed error has the same nature as sequence scrambling

(2) The authors introduce two important lengths LS1 and LS2 only in the conclusions and do not explain enough which each of them is important. It would make sense to discuss this early in the manuscript.

(3) It is not entirely clear why specific length distribution for VCG oligomers has to be assumed rather than emerged from simulations.

(4) Furthermore, the problem has another important length, L0 that is never introduced or discussed: a minimal hybridization length with a lifetime longer than the ligation time. From the parameters given, it appears that L0 is sufficiently long (~10 bases). In other words, it appears that the study is done is a somewhat suboptimal regime: most hybridization events do not lead to a ligation. Am I right in this assessment? If that is the case, the authors might want to explore another regime, L0<br /> Strengths:

High-quality theoretical modeling of an important problem is implemented.

Reviewer #2 (Public review):

Summary:

This important theoretical and computational study by Burger and Gerland attempts to set environmental, compositional, kinetic, and thermodynamic constraints on the proposed virtual circular genome (VCG) model for the early non-enzymatic replication of RNA. The authors create a solid kinetic model using published kinetic and thermodynamic parameters for non-enzymatic RNA ligation and (de)hybridization, which allows them to test a variety of hypotheses about the VCG. Prominently, the authors find that the length (longer is better) and concentration (intermediate is better) of the VCG oligos have an outsized impact on the fidelity and yield of VCG production with important implications for future VCG design. They also identify that activation of only RNA monomers, which can be achieved using environmental separation of the activation and replication, can relax the constraints on the concentration of long VCG component oligos by avoiding the error-prone oligo-oligo ligation. Finally, in a complex scenario with multiple VCG oligo lengths, the authors demonstrate a clear bias for the extension of shorter oligos compared to the longer ones. This effect has been observed experimentally (Ding et al., JACS 2023) but was unexplained rigorously until now. Overall, this manuscript will be of interest to scientists studying the origin of life and the behavior of complex nucleic acid systems.

Strengths:

- The kinetic model is carefully and realistically created, enabling the authors to probe the VCG thoroughly.<br /> - Fig. 6 outlines important constraints for scientists studying the origin of life. It supports the claim that the separation of activation and replication chemistry is required for efficient non-enzymatic replication. One could easily imagine a scenario where activation of molecules occurs, followed by their diffusion into another environment containing protocells that encapsulate a VCG. The selective diffusion of activated monomers across protocell membranes would then result in only activated monomers being available to the VCG, which is the constraint outlined in this work. The proposed exclusive replication by monomers also mirrors the modern biological systems, which nearly exclusively replicate by monomer extension.<br /> - Another strength of the work is that it explains why shorter oligos extend better compared to the long ones in complex VCG mixtures. This point is independent of the activation chemistry used (it simply depends on the kinetics and thermodynamics of RNA base-pairing) so it should be very generalizable.

Reviewer #1 (Public review):

Summary:

The authors have assembled a cohort of 10 SiNET, 1 SiAdeno, and 1 lung MiNEN samples to explore the biology of neuroendocrine neoplasms. They employ single-cell RNA sequencing to profile 5 samples (siAdeno, SiNETs 1-3, MiNEN) and single-nuclei RNA sequencing to profile seven frozen samples (SiNET 4-10).

They identify two subtypes of siNETs, characterized by either epithelial or neuronal NE cells, through a series of DE analyses. They also report findings of higher proliferation in non-malignant cell types across both subtypes. Additionally, they identify a potential progenitor cell population in a single lung MiNEN samples.

In the revised study, they have addressed my points and I have no further comments.

Reviewer #2 (Public review):

Summary:

The research identifies two main SiNET subtypes (epithelial-like and neuronal-like) and reveals heterogeneity in non-neuroendocrine cells within the tumor microenvironment. The study validates findings using external datasets and explores unexpected proliferation patterns. While it contributes to understanding SiNET oncogenic processes, the limited sample size and depth of analysis present challenges to the robustness of the conclusions.

Strengths:

The studies effectively identified two subtypes of SiNET based on epithelial and neuronal markers. Key findings include the low proliferation rates of neuroendocrine (NE) cells and the role of the tumor microenvironment (TME), such as the impact of Macrophage Migration Inhibitory Factor (MIF).

Weaknesses:

However, the analysis faces challenges such as a small sample size and lack of clear biological interpretation in some analyses.

Reviewer #3 (Public review):

This study profiles small intestine NETs and one mixed lung NET at single cell resolution and identifies two subtypes of neuroendocrine cells, as well as explores the proliferation patterns in malignant and nonmalignant cell types, identifying MIF as a potential factor that promotes proliferation of B and plasma cells in siNETs. Furthermore, they explore the single-cell landscape of a mixed LCNEC and squamous cell carcinoma, from which they identify a putative stem cell population with expression of features from both lineages.

Strengths:

This work showcases single-cell profiling of a rare tumor type, which is very informative for the field of NETs. The authors highlight very interesting observations, including the identification of the epithelial and neuronal subtype of siNETs, which they validated with an independent bulk RNA sequencing cohort. Furthermore, the observation of low cycling in malignant cells and high cycling in nonmalignant cells is an interesting one which may be applicable to other NETs.

Weaknesses:

• The authors do not connect their findings to clinical outcome. For example, is the epithelial or neuronal subtype enriched in tumors with worse or better prognosis or high grade vs. low grade siNETs or in patients who metastasize vs. who don't? As the authors show they can identify epithelial vs. neuronal subtypes in bulk RNA seq, perhaps they can take advantage of these other studies with larger sample sizes to investigate this. Additionally, the authors identify that the phenomenon of higher B/plasma cell proliferation is particular to epithelial siNETs and write that "The implications of high B/plasma cell turnover, and of other downstream effects of high MIF expression, are unclear, but raise the possibility that MIF-CD74 interaction may constitute a relevant target for the epithelial-like SiNET subtype." However, if this interaction contributes to survival in these patients, targeting this interaction may not be beneficial. Thus, it is important for the authors to try to connect their finding to clinical outcomes to enhance the translational relevance of this paper.

• The generalizability of this study would be enhanced if the authors analyzed other available single cell studies of NETs and found a similar phenomenon of high proliferating nonmalignant cell types. Although these studies are also very limited in sample size, seeing concordance in findings across independent cohorts and different experimental techniques would help to strengthen the findings. While the authors rationalize that these other studies are too distinct from their own due to enrichment for immune cells, this limitation should be noted but does not prevent such an analysis from being attempted.

• On page 3, the authors claim that "Technical effects (e.g. single cell analysis of fresh samples vs. single nuclei analysis of frozen samples) could also impact the capture of distinct cell types, although we did not observe a clear pattern of such bias." Can the authors show that cell type frequencies are not significantly different between the samples profiled with these two methods?

• Why did siNET3 and siNET9 have much lower recovery of neuroendocrine cells compared to other samples? It would be interesting to see how similar or different the transcriptional profiles are of the samples that were obtained from the same patient, considering that multifocal siNETs are found to derive from distinct clones, although this analysis is understandably not possible in this case due to the lack of neuroendocrine cells in one of two samples from the same patient.

• It should be more clearly stated in the text that these samples were previously treated with somatostatin analogues, as this impacts the interpretation of the findings.

• The identification of a potential progenitor subtype in the miNEN is very intriguing, albeit a case study and represents a distinct cancer from the lowly proliferating siNETs. While the authors mention this in the text, the case study feels rather tangential to the other parts of the paper.

• How the authors compared the DE genes to known signatures for the fibroblast and endothelial cells should be clarified and discussed in the Methods section.

Reviewer #2 (Public review):

Summary:

In this manuscript by Peto et al., the authors describe the impact of different antimicrobials on gut microbiota in a prospective observational study of 225 participants (healthy volunteers, inpatients and outpatients). Both cross-sectional data (all participants) and longitudinal data (subset of 79 haematopoietic cell transplant patients) were used. Using metagenomic sequencing, they estimated the impact of antibiotic exposure on gut microbiota composition and resistance genes. In their models, the authors aim to correct for potential confounders (e.g. demographics, non-antimicrobial exposures and physiological abnormalities), and for differences in the recency and total duration of antibiotic exposure. I consider these comprehensive models an important strength of this observational study. Yet, the underlying assumptions of such models may have impacted the study findings and residual confounding is likely. Other strengths include the presence of both cross-sectional and longitudinal exposure data and presence of both healthy volunteers and patients. Together, these observational findings expand on previous studies (both observational and RCTs) describing the impact of antimicrobials on gut microbiota.

Weaknesses:

(1) The main weaknesses result from the observational design. This hampers causal interpretation and makes correction for potential confounding necessary. While the authors have used comprehensive models to correct for potential confounders and for differences between participants in duration of antibiotic exposure and time between exposure and sample collection, I believe residual confounding is likely (which is mentioned as a limitation in the discussion).<br /> For their models, the authors found a disruption half-life of 6 days to be the best fit based on Shannon diversity. Yet, the disruption caused by antimicrobials may be longer than represented in this model - as highlighted in the discussion.

(2) Another consequence of the observational design of this study is the relatively small number of participants available for some comparisons (e.g. oral clindamycin was only used by 6 participants). Care should be taken when drawing any conclusions from such small numbers.

Comments on revisions:

The authors have adequately addressed all of my comments.

Reviewer #1 (Public review):

Summary:

The study by Cao et al. provides a compelling investigation into the role of mutational input in the rapid evolution of pesticide resistance, focusing on the two-spotted spider mite's response to the recent introduction of the acaricide cyetpyrafen. This well-documented introduction of the pesticide-and thus a clearly defined history of selection-offers a powerful framework for studying the temporal dynamics of rapid adaptation. The authors combine resistance phenotyping across multiple populations, extensive resequencing to track the frequency of resistance alleles, and genomic analyses of selection in both contemporary and historical samples. These approaches are further complemented by laboratory-based experimental evolution, which serves as a baseline for understanding the genetic architecture of resistance across mite populations in China. Their analyses identify two key resistance-associated genes, sdhB and sdhD, within which they detect 15 mutations in wild-collected samples. Protein modeling reveals that these mutations cluster around the pesticide's binding site, suggesting a direct functional role in resistance. The authors further examine signatures of selective sweeps and their distribution across populations to infer the mechanisms-such as de novo mutation or gene flow-driving the spread of resistance, a crucial consideration for predicting evolutionary responses to extreme selection pressure. Overall, this is a well-rounded, thoughtfully designed and well-written manuscript. It shows significant novelty, as it is relatively rare to integrate broad-scale evolutionary inference from natural populations with experimentally informed bioassays, however, follow up work will be needed to fully resolve haplotype structure and the functional effects of resistance mutations in the system.

Strengths:

One of the most compelling aspects of this study is its integration of genomic time-series data in natural populations with controlled experimental evolution. By coupling genome sequencing of resistant field populations with laboratory selection experiments, the authors tease apart the individual effects of resistance alleles along with regions of the genome where selection is expected to occur, and compare that to the observed frequency in the wild populations over space and time. Their temporal data clearly demonstrates the pace at which evolution can occur in response to extreme selection. This type of approach is a powerful roadmap for the rest of the field of rapid adaptation.

The study effectively links specific genetic changes to resistance phenotypes. The identification of sdhB and sdhD mutations as major drivers of cyetpyrafen resistance is well supported by allele frequency shifts in both field and experimental populations. The scope of their sampling clearly facilitated the remarkable number of observed mutations within these target genes, and the authors provide a careful discussion of the likelihood of these mutations from de novo or standing variation. Furthermore, the discovered cross-resistance that these mutations confer to other mitochondrial complex II inhibitors highlights the potential for broader resistance management and evolution.

Weaknesses:

(1) Pleiotropy without pesticide modes of action (cyflumetofen and cyetpyrafen) may also play a role in the rapid response to the focal pesticide in this study<br /> (2) Other aspects of the environment that might influence selection were not considered in the structure of resistance alleles (i.e. climate, elevation)<br /> (3) Very little data were used for haplotype reconstruction, only 8 SNPs, and this excluded all heterozygous alleles, which could dramatically influence the complexity of these inferred haplotype networks.<br /> (4) Single Mutations and Their Effects:<br /> - Allelic effects were not estimated in isogenic lines, so the effects presented also include heterogeneity from allelic interactions with the genomic background<br /> - The authors see populations that segregate for resistance mutations but that have no survival to pesticides. This suggests either not all of the resistance mutations studied here actually have functional effects or that dominance is playing a role in masking their effects in the heterozygous state.

Reviewer #2 (Public review):

Summary:

This paper investigates the evolution of pesticide resistance in the two-spotted spider mite following the introduction of an SDHI acaricide, cyatpyrafen, in China. The authors make use of cyatpyrafen-naive populations collected before that pesticide was first used, as well as more recent populations (both sensitive and resistant) to conduct comparative population genomics. They report 15 different mutations in the insecticide target site from resistant populations, many reported here for the first time, and look at the mutation and selection processes underlying the evolution of resistance, through GWAS, haplotype mapping, and testing for loss of diversity indicating selective sweeps. None of the target site mutations found in resistant populations was found in pre-exposure populations, suggesting that the mutations may have arisen de novo rather than being present as standing variation, unless initially present at very low frequencies; a de novo origin is also supported by evidence of selective sweeps in some resistant populations. Furthermore, there is no significant evidence of migration of resistant genotypes between the sampled field populations indicating multiple origins of common mutations. Overall, this indicates a very high mutation rate and a wide range of mutational pathways to resistance for this target site in this pest species. The series of population genomic analyses carried out here, in addition to the evolutionary processes that appear to underly resistance development in this case, could have implications for the study of resistance evolution more widely.

Strengths:

This paper combines phenotypic characterisation with extensive comparative population genomics, made possible by the availability of multiple population samples (each with hundreds of individuals) collected before as well as after then introduction of the pesticide cyatpyrafen, as well as lab-evolved lines. This resuts in findings of mutation and selection processes that can be related back to the pesticide resistance trait of concern. Large numbers of mites were tested phenotypically to show the levels of resistance present, and the authors also made near-isogenic lines to confirm the phenotypic effects of key mutations. The population genomic analyses consider a range of alternative hypotheses, including mutations arising by de novo mutation or selection from standing genetic variation; and mutations in different populations arising independently or arriving by migration. The claim that mutations most likley arose by multiple repeated de novo mutations is therefore supported by multiple lines of evidence: the direct evidence of none of the mutations being found in over 2000 individuals from naive populations, and the indirect evidence from population genomics showing evidence of selective sweeps but not of significant migration between the sampled populations.

Weaknesses:

As acknowledged within the discussion, whilst evidence supports a de novo origin of the resistance associated mutations, this cannot be proven definitively as mutations may have been present at a very low frequency and therefore not found within the tested pesticide-naive population samples.

Near-isofemale lines were made to confirm the resistance levels associated with five of the 15 mutations, but otherwise the genotype-phenotype associations are correlative as confirmation by functional genetics was beyond the scope of this study.

Comments on revisions:

My recommendations have all been addressed in the revised version.

Reviewer #1 (Public review):

In this manuscript, Azlan et al. identified a novel maternal factor called Sakura that is required for proper oogenesis in Drosophila. They showed that Sakura is specifically expressed in the female germline cells. Consistent with its expression pattern, Sakura functioned autonomously in germline cells to ensure proper oogenesis. In sakura KO flies, germline cells were lost during early oogenesis and often became tumorous before degenerating by apoptosis. In these tumorous germ cells, piRNA production was defective and many transposons were derepressed. Interestingly, Smad signaling, a critical signaling pathway for the GSC maintenance, was abolished in sakura KO germline stem cells, resulting in ectopic expression of Bam in whole germline cells in the tumorous germline. A recent study reported that Bam acts together with the deubiquitinase Otu to stabilize Cyc A. In the absence of sakura, Cyc A was upregulated in tumorous germline cells in the germarium. Furthermore, the authors showed that Sakura co-immunoprecipitated Otu in ovarian extracts. A series of in vitro assays suggested that the Otu (1-339 aa) and Sakura (1-49 aa) are sufficient for their direct interaction. Finally, the authors demonstrated that the loss of otu phenocopies the loss of sakura, supporting their idea that Sakura plays a role in germ cell maintenance and differentiation through interaction with Otu during oogenesis.

Latest comments:

The reviewer acknowledges the importance of sharing the observed defects in Sakura mutant ovaries and the possible physiological significance of the Sakura-Out interaction with the research community, as this information could lay the groundwork for future functional analysis research.

Reviewer #2 (Public review):

In this study, the authors identified CG14545 (named it sakura), as a key gene essential for Drosophila oogenesis. Genetic analyses revealed that Sakura is vital for both oogenesis progression and ultimate female fertility, playing a central role in the renewal and differentiation of germ stem cells (GSC).

The absence of Sakura disrupts the Dpp/BMP signaling pathway, resulting in abnormal bam gene expression, which impairs GSC differentiation and leads to GSC loss. Additionally, Sakura is critical for maintaining normal levels of piRNAs. Also, the authors convincingly demonstrate that Sakura physically interacts with Otu, identifying the specific domains necessary for this interaction, suggesting a cooperative role in germline regulation. Importantly, the loss of otu produces similar defects to those observed in sakura mutants, highlighting their functional collaboration.

The authors provide compelling evidence that Sakura is a critical regulator of germ cell fate, maintenance, and differentiation in Drosophila. This regulatory role is mediated through modulation of pMad and Bam expression. However, the phenotypes observed in the germarium appear to stem from reduced pMad levels, which subsequently trigger premature and ectopic expression of Bam. This aberrant Bam expression could lead to increased CycA levels and altered transcriptional regulation, impacting piRNA expression. In this revised manuscript, the authors further investigated whether Sakura affects the function of Orb, a binding partner they identified, in deubiquitinase activity when Orb interacts with Bam.

This elaborate study will be embraced by both germline-focused scientists and the developmental biology community.

Latest comments:

The authors answered all my persistent concerns and made changes according to the recommendations I incorporated for the revised version of the manuscript.

Reviewer #3 (Public review):

In this very thorough study, the authors characterize the function of a novel Drosophila gene, which they name Sakura. They start with the observation that sakura expression is predicted to be highly enriched in the ovary and they generate an anti-sakura antibody, a line with a GFP-tagged sakura transgene, and a sakura null allele to investigate sakura localization and function directly. They confirm the prediction that it is primarily expressed in the ovary and, specifically, that it is expressed in germ cells, and find that about 2/3 of the mutants lack germ cells completely and the remaining have tumorous ovaries. Further investigation reveals that Sakura is required for piRNA-mediated repression of transposons in germ cells. They also find evidence that sakura is important for germ cell specification during development and germline stem cell maintenance during adulthood. However, despite the role of sakura in maintaining germline stem cells, they find that sakura mutant germ cells also fail to differentiate properly such that mutant germline stem cell clones have an increased number of "GSC-like" cells. They attribute this phenotype to a failure in the repression of Bam by dpp signaling. Lastly, they demonstrate that sakura physically interacts with otu and that sakura and otu mutants have similar germ cell phenotypes. Overall, this study helps to advance the field by providing a characterization of a novel gene that is required for oogenesis. The data are generally high-quality and the new lines and reagents they generated will be useful for the field.

Latest comments:

As with my previous assessment, I remain supportive of publication of this manuscript. Though I agree with the other reviewers that additional experimentation would increase the value of this study even further, I feel it will also be a useful contribution to the field as is.

Reviewer #2 (Public review):

Summary:

The authors report the discovery of a population of gingival fibroblasts displaying the expression of cellular senescence markers P21 and P16 in human periodontitis samples and a murine ligature-induced periodontitis (LIP) model. They support this finding in the murine model through bulk RNA-sequencing and show that differentially expressed genes are significantly enriched in the SenMayo cellular senescence in aging dataset. They then show that Ligature-Induced Periodontitis (LIP) mice treated with the senomorphic drug metformin display overall diminished bone damage, reduced histomorphic alterations, and a reduction in P21 and P16 immunostaining signal. To explore the cell types expressing cellular senescence markers in periodontitis, the authors make use of a combination of bioinformatic analyses on publicly available scRNA-seq data, immunostainings on patient samples and their LIP model; as well as in vitro culture of healthy human gingival fibroblasts treated with LPS. They found that fibroblasts are a cell population expressing P16 in periodontitis which are also enriched for SenMayo genes, suggesting they have a senescent phenotype. They then point to a subgroup of fibroblasts expressing CD81+ with the highest enrichment for a SASP geneset in periodontitis. They also show that treatment of LIP mice and human LPS-treated gingival fibroblasts with metformin leads to a reduction of P21 and P16-positive cells, as well as the senescence-associated beta-galactosidase (SA-beta-gal) marker. Finally, they show evidence suggesting that CD81+ senescent fibroblasts are the source of C3 complement protein, which they stipulate signals through the C3AR1 receptor present in neutrophils in periodontitis. The authors observed that both CD81+ fibroblast and C3AR1+ neutrophil populations are expanded in periodontitis, that both populations appear to be in close contact, and that treatment with metformin reduced both C3 and the neutrophil marker MPO in their mouse LIP model.

After a round of revision, the authors have made significant improvements to their manuscript, such as improving the quality of the data/evidence and also included new data from experiments using a well-known senolytic and the senomorphic metformin, which all together provide a solid support to their main claims.

Strengths:

The study implements several different techniques and tools on human samples, mouse models, fibroblast cultures, and publicly available data to support their conclusions. In summary, they provide solid evidence showing that in the context of periodontitis, there is an expansion of cells expressing senescence markers P21, and P16, as well as members of the SASP, and that this includes CD81+ fibroblasts.

Weaknesses:

The fact that in this study the periodontitis samples belonged to patients with a significantly higher median age (all older than 50 years of age) and the healthy samples belonged to young adults (all younger than 35 years of age), raises the need for caution in interpretation due to a possible effect of aging in the accumulation of CD81+ senescent fibroblasts. However, the recruitment of similar age groups in this case is of course difficult due to the higher prevalence of periodontitis in older adults. In this regard it is important to note that the authors still support their findings using a mouse ligature model. Similar studies comparing healthy and periodontitic patients from similar age groups will be of great importance in the future.

Reviewer #1 (Public review):

Summary:

In this paper, Wang et al show that differentiated peridermal cells of the zebrafish epidermis extend cytoneme-like protrusions toward the less differentiated, intermediate layer below. They present evidence that expression of a dominant-negative cdc42, inhibits cytoneme formation and leads to elevated expression of a marker of undifferentiated keratinocytes, krtt1c19e, in the periderm layer. It is demonstrated that Delta-Notch signaling is involved in keratinocyte differentiation and that loss of cytonemes correlates with a loss of Notch signaling. Finally, changes in expression of the inflammatory cytokine IL-17 and its receptors is shown to affect cytoneme number and periderm structure in a manner similar to Notch and cdc42 perturbations.

Strengths:

Overall, the idea that differentiated cells signal to underlying undifferentiated cells via membrane protrusions in skin keratinocytes is interesting and novel, and it is clear that periderm cells send out thin membrane protrusions that contain a Notch ligand. Further, and perturbations that affect cytoneme number, Notch signaling and IL-17 expression clearly lead to changes in periderm structure and gene expression.

Weaknesses:

The mechanisms by which IL-17 affects cytoneme formation requires further investigation.

Reviewer #2 (Public review):

Summary:

The aim of the study was to understand how cells of the skin communicate across dermal layers. The research group has previously demonstrated that cellular connections called airinemes contribute to this communication. The current work builds upon this knowledge by showing that differentiated keratinocytes also use cytonemes, specialized signaling filopodia, to communicate with undifferentiated keratinocytes. They show that cytonemes are the more abundant type of cellular extension used for communication between the differentiated keratinocyte layer and the undifferentiated keratinocytes. Disruption of cytoneme formation led to expansion of the undifferentiated keratinocytes into the periderm, mimicking skin diseases like psoriasis. The authors go on to show that disruption of cytonemes results in perturbations in Notch signaling between the differentiated keratinocytes of the periderm and the underlying proliferating undifferentiated keratinocytes. Further the authors show that Interleukin-17, also known to drive psoriasis, can restrict formation of periderm cytonemes, possibly through the inhibition of Cdc42 expression. This work suggests that cytoneme mediated Notch signaling plays a central role in normal epidermal regulation. The authors propose that disruption of cytoneme function may be an underlying cause of various human skin diseases.

Strengths:

The authors provide strong evidence that periderm keratinocytes cytonemes contain the notch ligand DeltaC to promote Notch activation in the underlying intermediate layer to regulate accurate epidermal maintenance.

Weaknesses:

The impact of the study would be increased if the mechanism by which Interlukin-17 and Cdc42 collaborate to regulate cytonemes was defined. Experiments measuring Cdc42 activity, rather than just measuring expression, would strengthen the conclusions.

Comments on revisions:

The authors have sufficiently addressed my critiques from the initial round of evaluation. They have included useful representative images, clarified how they scored cytonemes and provided additional controls/experimental conditions that improve the rigor of the study. The results provided now support the key conclusions of the study.

Reviewer #1 (Public review):

Summary:

This fundamental study presents a practical modification of the orthogonal hybridization chain reaction (HCR) technique, a promising yet underutilized method with broad potential for future applications across various fields. The authors advance this technique by integrating peptide ligation technology and nanobody-based antibody mimetics-cost-effective and scalable alternatives to conventional antibodies-into a DNA-immunoassay framework that merges oligonucleotide-based detection with immunoassay methodologies. They demonstrate this with compelling evidence that this approach facilitates a modified ELISA platform capable of simultaneously quantifying multiple target protein expression levels within a single protein mixture sample.

Strengths:

The hybridization chain reaction (HCR) technique was initially developed to enable the simultaneous detection of multiple mRNA expression levels within the same tissue. This method has since evolved into immuno-HCR, which extends its application to protein detection by utilizing antibodies. A key requirement of immuno-HCR is the coupling of oligonucleotides to antibodies, a process that can be challenging due to the inherent difficulties in expressing and purifying conventional antibodies.<br /> In this study, the authors present an innovative approach that circumvents these limitations by employing nanobody-based antibody mimetics, which recognize antibodies, instead of directly coupling oligonucleotides to conventional antibodies. This strategy facilitates oligonucleotide conjugation-designed to target the initiator hairpin oligonucleotide of HCR-through peptide ligation and click chemistry.

Weaknesses:

The sandwich-format technique presented in this study, which employs a nanobody that recognizes primary IgG antibodies, may have limited scalability compared to existing methods that directly couple oligonucleotides to primary antibodies. This limitation arises because the C-region types of primary antibodies are relatively restricted, meaning that the use of nanobody-based detection may constrain the number of target proteins that can be analyzed simultaneously. In contrast, the conventional approach of directly conjugating oligonucleotides to primary antibodies allows for a broader range of protein targets to be analyzed in parallel.

Additionally, in the context of HCR-based protein detection, the number of proteins that can be analyzed simultaneously is inherently constrained by fluorescence wavelength overlap in microscopy, which limits its multiplexing capability. By comparison, direct coupling oligonucleotides to primary antibodies can facilitate the simultaneous measurement of a significantly greater number of protein targets than the sandwich-based nanobody approach in the barcode-ELISA/NGS-based technique.

Comments on revisions:

The previous suggestions were well incorporated in the revised manuscript.

for - paper - climate crisis - rebound effect - paper - title - Energy efficiency and economy-wide rebound effects: A review of the evidence and its implications - from - post - LinkedIn - rebound effect - https://hyp.is/yz4m_ldBEfC18Bfg0RPf2w/www.linkedin.com/feed/update/urn:li:activity:7346027213776953344/

Reviewer #2 (Public review):

Summary:

This study explores how a simple sensorimotor circuit in the nematode C. elegans enables it to navigate salt gradients based on past experiences. Using computational simulations and previously described neural connections, the study demonstrates how a single neuron, ASER, can change its signaling behavior in response to different salt conditions, with which the worm is able to "remember" prior environments and adjust its navigation toward "preferred" salinity accordingly.

Strengths:

The key novelty and strength of this paper is the explicit demonstration of computational neurobehavioral modeling and evolutionary algorithms to elucidate the synaptic plasticity in a minimal neural circuit that is sufficient to replicate memory-based chemotaxis. In particular, with changes in ASER's glutamate release and sensitivity of downstream neurons, the ASER neuron adjusts its output to be either excitatory or inhibitory depending on ambient salt concentration, enabling the worm to navigate toward or away from salt gradients based on prior exposure to salt concentration.

Weaknesses:

While the model successfully replicates some behaviors observed in previous experiments, some key assumptions of the work still need to be verified by biological validation of further experiments.

Comments on revisions:

Thank you for the authors' response. The revision and their response have substantially addressed my concerns.

Reviewer #1 (Public review):

Summary:

In the present study, Chen et al. investigate the role of Endophilin A1 in regulating GABAergic synapse formation and function. To this end, the authors use constitutive or conditional knockout of Endophilin A1 (EEN1) to assess the consequences on GABAergic synapse composition and function, as well as the outcome for PTZ-induced seizure susceptibility. The authors show that EEN1 KO mice show a higher susceptibility to PTZ-induced seizures, accompanied by a reduction in the GABAergic synaptic scaffolding protein gephyrin as well as specific GABAAR subunits and eIPSCs. The authors then investigate the underlying mechanisms, demonstrating that Endophilin A1 binds directly to gephyrin and GABAAR subunits, and identifying the subdomains of Endophilin A1 that contribute to this effect. Overall, the authors state that their study places Endophilin A1 as a new regulator of GABAergic synapse function.

Strengths:

Overall, the topic of this manuscript is very timely, since there has been substantial recent interest in describing the mechanisms governing inhibitory synaptic transmission at GABAergic synapses. The study will therefore be of interest to a wide audience of neuroscientists studying synaptic transmission and its role in disease. The manuscript is well written and contains a substantial quantity of data. In the revised version of the manuscript, the authors have increased the number of samples analyzed and have significantly improved the statistical analysis, thereby substantially strengthening the conclusions of their study.

Comments on revised version:

The authors have addressed all of my concerns, and the manuscript has been substantially improved.

Reviewer #2 (Public review):

Summary:

The function of neural circuits relies heavily on the balance of excitatory and inhibitory inputs. Particularly, inhibitory inputs are understudied when compared to their excitatory counterparts due to the diversity of inhibitory neurons, their synaptic molecular heterogeneity, and their elusive signature. Thus, insights into these aspects of inhibitory inputs can inform us largely on the functions of neural circuits and the brain.

Endophilin A1, an endocytic protein heavily expressed in neurons, has been implicated in numerous pre- and postsynaptic functions, however largely at excitatory synapses. Thus, whether this crucial protein plays any role in inhibitory synapse, and whether this regulates functions at the synaptic, circuit, or brain level remains to be determined.

New Findings:

(1) Endophilin A1 interacts with the postsynaptic scaffolding protein gephyrin at inhibitory postsynaptic densities within excitatory neurons.

(2) Endophilin A1 promotes the organization of the inhibitory postsynaptic density and the subsequent recruitment/stabilization of GABA A receptors via Endophilin A1's membrane binding and actin polymerization activities.

(3) Loss of Endophilin A1 in CA1 mouse hippocampal pyramidal neurons weakens inhibitory input and leads to susceptibility to epilepsy.

(4) Thus the authors propose that via its role as a component of the inhibitory postsynaptic density within excitatory neurons, Endophilin A1 supports the organization, stability, and efficacy of inhibitory input to maintain the excitatory/inhibitory balance critical for brain function.

(5) The conclusion of the manuscript is well supported by the data but will be strengthened by addressing our list of concerns and experiment suggestions.

Comments on revised version:

The authors addressed the concerns adequately. The three remaining concerns are:

(1) The use of one-way ANOVA is not well justified.

(2) The use of superplots to show culture to culture variability would make it more transparent.

(3) Change EEN1 in Figure 8B to EndoA1.

Reviewer #3 (Public review):

Summary:

The authors investigated a possible role of Endophilin A1 in the inhibitory postsynaptic density.

Strengths:

The authors used a broad array of experimental approaches to investigate this, including tests of seizure susceptibility, electrophysiology, biochemistry, neuronal culture and image analysis.

Weaknesses:

Many results are difficult to interpret, and data quality is not always convincing, unfortunately. The basic premise of the study, that gephyrin and endophilin A1 interact, requires more robust analysis to be convincing.

Specific comments:

The authors have made a substantial effort to improve their manuscript. A number of issues, related to numbers of observations mentioned by the reviewers, are clarified in the revised manuscript. The authors have also clarified some of the other questions from the reviewers. The long list of issues brought up by the reviewers and the many corrections needed still raise questions about data quality in this manuscript.<br /> In response to my comments (Point 2), the added experiment with PSD95.FingR and GPN.FingR in cultured neurons (Fig. S5A-D) is a good addition; the in vivo data using FingRs in Figure S3 look less convincing however. In response to my Point 5, the authors have added a cell-free binding assay (Figure 5I). This is a useful addition, but to convincingly make the point of interaction between Gephyrin and EndoA1, more rigorous biophysical quantitation of binding is needed. The legend in Figure 5I states that 4 independent experiments were performed, but the graph only shows 3 dots. This needs to be corrected.

Reviewer #1 (Public review):

Summary:

The manuscript by Ross, Miscik, and others describes an intriguing series of observations made when investigating the requirement for podxl during hepatic development in zebrafish. Podxl morphants and CRISPants display a reduced number of hepatic stellate cells (HSCs), while mutants are either phenotypically wild type or display an increased number of HSCs.

The absence of observable phenotypes in genetic mutants could indeed be attributed to genetic compensation, as the authors postulate. However, in my opinion, the evidence provided in the manuscript at this point is insufficient to draw a firm conclusion. Furthermore, the opposite phenotype observed in the two deletion mutants is not readily explainable by genetic compensation and invokes additional mechanisms.

Major concerns:

(1) Considering discrepancies in phenotypes, the phenotypes observed in podxl morphants and CRISPants need to be more thoroughly validated. To generate morphants, authors use "well characterized and validated ATG Morpholino" (lines 373-374). However, published morphants, in addition to kidney malformations, display gross developmental defects including pericardial edema, yolk sack extension abnormalities, and body curvature at 2-3 dpf (reference 7 / PMID: 24224085). Were these gross developmental defects observed in the knockdown experiments performed in this paper? If yes, is it possible that the liver phenotype observed at 5 dpf is, to some extent, secondary to these preceding abnormalities? If not, why were they not observed? Did kidney malformations reproduce? On the CRISPant side, were these gross developmental defects also observed in sgRNA#1 and sgRNA#2 CRISPants? Considering that morphants and CRISPants show very similar effects on HSC development and assuming other phenotypes are specific as well, they would be expected to occur at similar frequencies. It would be helpful if full-size images of all relevant morphant and CRISPant embryos were displayed, as is done for tyr CRISPant in Figure S2. Finally, it is very important to thoroughly quantify the efficacy of podxl sgRNA#1 and sgRNA#2 in CRISPants. The HRMA data provided in Figure S1 is not quantitative in terms of the fraction of alleles with indels. Figure S3 indicates a very broad range of efficacies, averaging out at ~62% (line 100). Assuming random distribution of indels among cells and that even in-frame indels result in complete loss of function (possible for sgRNA#1 due to targeting the signal sequence), only ~38% (.62*.62) of all cells will be mutated bi-allelically. That does not seem sufficient to reliably induce loss-of-function phenotypes. My guess is that the capillary electrophoresis method used in Figure S3 underestimates the efficiency of mutagenesis, and that much higher mutagenesis rates would be observed if mutagenesis were assessed by amplicon sequencing (ideally NGS but Sanger followed by deconvolution analysis would suffice). This would strengthen the claim that CRISPant phenotypes are specific.

(2) In addition to confidence in morphant and CRISPant phenotypes, the authors' claim of genetic compensation rests on the observation that podxl (Ex1(p)_Ex7Δ) mutants are resistant to CRISPant effect when injected with sgRNA#1 (Figure 3L). Considering the issues raised in the paragraph above, this is insufficient. There is a very straightforward way to address both concerns, though. The described podxl(-194_Ex7Δ) and podxl(-319_ex1(p)Δ) deletions remove the binding site for the ATG morpholino. Therefore, deletion mutants should be refractive to the Morpholino (specificity assessment recommended in PMID: 29049395, see also PMID: 32958829). Furthermore, both deletion mutants should be refractive to sgRNA#1 CRISPant phenotypes, with the first being refractive to sgRNA#2 as well.

Reviewer #2 (Public review):

In this manuscript, Ross and Miscik et. al described the phenotypic discrepancies between F0 zebrafish mosaic mutant ("CRISPants") and morpholino knockdown (Morphant) embryos versus a set of 5 different loss-of-function (LOF) stable mutants in one particular gene involved in hepatic stellate cells development: podxl. While transient LOF and mosaic mutants induced a decrease of hepatic stellate cells number stable LOF zebrafish did not. The authors analyzed the molecular causes of these phenotypic differences and concluded that LOF mutants are genetically compensated through the upregulation of the expression of many genes. Additionally, they ruled out other better-known and described mechanisms such as the expression of redundant genes, protein feedback loops, or transcriptional adaptation.

While the manuscript is clearly written and conclusions are, in general, properly supported, there are some aspects that need to be further clarified and studied.

(1) It would be convenient to apply a method to better quantify potential loss-of-function mutations in the CRISPants. Doing this it can be known not only percentage of mutations in those embryos but also what fraction of them are actually generating an out-of-frame mutation likely driving gene loss of function (since deletions of 3-6 nucleotides removing 1-2 aminoacid/s will likely not have an impact in protein activity, unless that this/these 1-2 aminoacid/s is/are essential for the protein activity). With this, the authors can also correlate phenotype penetrance with the level of loss-of-function when quantifying embryo phenotypes that can help to support their conclusions.

(2) It is unclear that 4.93 ng of morpholino per embryo is totally safe. The amount of morpholino causing undesired effects can differ depending on the morpholino used. I would suggest performing some sanity check experiments to demonstrate that morpholino KD is not triggering other molecular outcomes, such as upregulation of p53 or innate immune response.

(3) Although the authors made a set of controls to demonstrate the specificity of the CRISPant phenotypes, I believe that a rescue experiment could be beneficial to support their conclusions. Injecting an mRNA with podxl ORF (ideally with a tag to follow protein levels up) together with the induction of CRISPants could be a robust manner to demonstrate the specificity of the approach. A rescue experiment with morphants would also be good to have, although these are a bit more complicated, to ultimately demonstrate the specificity of the approach.

(4) In lines 314-316, the authors speculate on a correlation between decreased HSC and Podxl levels. It would be interesting to actually test this hypothesis and perform RT-qPCR upon CRISPant induction or, even better and if antibodies are available, western blot analysis.

(5) Similarly, in lines 337-338 and 342-344, the authors discuss that it could be possible that genes near to podxl locus could be upregulated in the mutants. Since they already have a transcriptomic done, this seems an easy analysis to do that can address their own hypothesis.

(6) Figures 4 and 5 would be easier to follow if panels B-F included what mutants are (beyond having them in the figure legend). Moreover, would it be more accurate and appropriate if the authors group all three WT and mutant data per panel instead of showing individual fish? Representing technical replicates does not demonstrate in vivo variability, which is actually meaningful in this context. Then, statistical analysis can be done between WT and mutant per panel and per set of primers using these three independent 3-month-old zebrafish.

Reviewer #3 (Public review):

Summary:

Ross et al. show that knockdown of zebrafish podocalyxin-like (podxl) by CRISPR/Cas or morpholino injection decreased the number of hepatic stellate cells (HSC). The authors then generated 5 different mutant alleles representing a range of lesions, including premature stop codons, in-frame deletion of the transmembrane domain, and deletions of the promoter region encompassing the transcription start site. However, unlike their knockdown experiment, HSC numbers did not decrease in podxl mutants; in fact, for two of the mutant alleles, the number of HSCs increased compared to the control. Injection of podxl CRISPR/Cas constructs into these mutants had no effect on HSC number, suggesting that the knockdown phenotype is not due to off-target effects but instead that the mutants are somehow compensating for the loss of podxl. The authors then present multiple lines of evidence suggesting that compensation is not exclusively due to transcriptional adaptation - evidence of mRNA instability and nonsense-mediated decay was observed in some but all mutants; expression of the related gene endoglycan (endo) was unchanged in the mutants and endo knockdown had no effect on HSC numbers; and, expression profiling by RNA sequencing did not reveal changes in other genes that share sequence similarity with podxl. Instead, their RNA-seq data showed hundreds of differentially expressed genes, especially ECM-related genes, suggesting that compensation in podxl mutants is complex and multi-genic.

Strengths:

The data presented is impressively thorough, especially in its characterization of the 5 different podxl alleles and exploration of whether these mutants exhibit transcriptional adaptation.

Weaknesses:

RNA sequencing expression profiling was done on adult livers. However, compensation of HSC numbers is apparent by 6 dpf, suggesting compensatory mechanisms would be active at larval or even embryonic stages. Although possible, it's not clear that any compensatory changes in gene expression would persist to adulthood.

Reviewer #1 (Public review):

This study presents cryoEM-derived structures of the Trypanosome aquaporin AQP2, in complex with its natural ligand, glycerol, as well as two trypanocidal drugs, pentamidine and melarsoprol, which use AQP2 as an uptake route. The structures are high quality, and the density for the drug molecules is convincing, showing a binding site in the centre of the AQP2 pore.

The authors then continue to study this system using molecular dynamics simulations. Their simulations indicate that the drugs can pass through the pore and identify a weak binding site in the centre of the pore, which corresponds with that identified through cryoEM analysis. They also simulate the effect of drug resistance mutations, which suggests that the mutations reduce the affinity for drugs and therefore might reduce the likelihood that the drugs enter into the centre of the pore, reducing the likelihood that they progress through into the cell.

While the cryoEM and MD studies are well conducted, it is a shame that the drug transport hypothesis was not tested experimentally. For example, did they do cryoEM with AQP2 with drug resistance mutations and see if they could see the drugs in these maps? They might not bind, but another possibility is that the binding site shifts, as seen in Chen et al. Do they have an assay for measuring drug binding? I think that some experimental validation of the drug binding hypothesis would strengthen this paper. Without this, I would recommend the authors to soften the statement of their hypothesis (i.e, lines 65-68) as this has not been experimentally validated.

Reviewer #2 (Public review):

Summary:

The authors present 3.2-3.7 Å cryo-EM structures of Trypanosoma brucei aquaglyceroporin-2 (TbAQP2) bound to glycerol, pentamidine, or melarsoprol and combine them with extensive all-atom MD simulations to explain drug recognition and resistance mutations. The work provides a persuasive structural rationale for (i) why positively selected pore substitutions enable diamidine uptake, and (ii) how clinical resistance mutations weaken the high-affinity energy minimum that drives permeation. These insights are valuable for chemotherapeutic re-engineering of diamidines and aquaglyceroporin-mediated drug delivery.

My comments are on the MD part.

Strengths:

The study

(1) Integrates complementary cryo-EM, equilibrium, applied voltage MD simulations, and umbrella-sampling PMFs, yielding a coherent molecular-level picture of drug permeation.

(2) Offers direct structural rationalisation of long-standing resistance mutations in trypanosomes, addressing an important medical problem.

Weaknesses:

Unphysiological membrane potential. A field of 0.1 V nm⁻¹ (~1 V across the bilayer) was applied to accelerate translocation. From the traces (Figure 1c), it can be seen that the translocation occurred really quickly through the channel, suggesting that the field might have introduced some large changes in the protein. The authors state that they checked visually for this, but some additional analysis, especially of the residues next to the drug, would be welcome.

Based on applied voltage simulations, the authors argue that the membrane potential would help get the drug into the cell, and that a high value of the potential was applied merely to speed up the simulation. At the same time, the barrier for translocation from PMF calculations is ~40 kJ/mol for WT. Is the physiological membrane voltage enough to overcome this barrier in a realistic time? In this context, I do not see how much value the applied voltage simulations have, as one can estimate the work needed to translocate the substrate on PMF profiles alone. The authors might want to tone down their conclusions about the role of membrane voltage in the drug translocation.

Pentamidine charge state and protonation. The ligand was modeled as +2, yet pKa values might change with the micro-environment. Some justification of this choice would be welcome.

I don't follow the RMSD calculations. The authors state that this RMSD is small for the substrate and show plots in Figure S7a, with the bottom plot being presumably done for the substrate (the legends are misleading, though), levelling off at ~0.15 nm RMSD. However, in Figure S7a, we see one trace (light blue) deviating from the initial position by more than 0.2 nm - that would surely result in an RMSD larger than 0.15, but this is somewhat not reflected in the RMSD plots.

Reviewer #3 (Public review):

Summary:

Recent studies have established that trypanocidal drugs, including pentamidine and melarsoprol, enter the trypanosomes via the glyceroaquaporin AQP2 (TbAQP2). Interestingly, drug resistance in trypanosomes is, at least in part, caused by recombination with the neighbouring gene, AQP3, which is unable to permeate pentamidine or melarsoprol. The effect of the drugs on cells expressing chimeric proteins is significantly reduced. In addition, controversy exists regarding whether TbAQP2 permeates drugs like an ion channel, or whether it serves as a receptor that triggers downstream processes upon drug binding. In this study the authors set out to achieve three objectives:<br /> (1) to determine if TbAQP2 acts as a channel or a receptor,<br /> (2) to understand the molecular interactions between TbAQP2 and glycerol, pentamidine, and melarsoprol, and<br /> (3) to determine the mechanism by which mutations that arise from recombination with TbAQP3 result in reduced drug permeation.

Indeed, all three objectives are achieved in this paper. Using MD simulations and cryo-EM, the authors determine that TbAQP2 likely permeates drugs like an ion channel. The cryo-EM structures provide details of glycerol and drug binding, and show that glycerol and the drugs occupy the same space within the pore. Finally, MD simulations and lysis assays are employed to determine how mutations in TbAQP2 result in reduced permeation of drugs by making entry and exit of the drug relatively more energy-expensive. Overall, the strength of evidence used to support the author's claims is solid.

Strengths:

The cryo-EM portion of the study is strong, and while the overall resolution of the structures is in the 3.5Å range, the local resolution within the core of the protein and the drug binding sites is considerably higher (~2.5Å).

I also appreciated the MD simulations on the TbAQP2 mutants and the mechanistic insights that resulted from this data.

Weaknesses:

(1) The authors do not provide any empirical validation of the drug binding sites in TbAQP2. While the discussion mentions that the binding site should not be thought of as a classical fixed site, the MD simulations show that there's an energetically preferred slot (i.e., high occupancy interactions) within the pore for the drugs. For example, mutagenesis and a lysis assay could provide us with some idea of the contribution/importance of the various residues identified in the structures to drug permeation. This data would also likely be very valuable in learning about selectivity for drugs in different AQP proteins.

(2) Given the importance of AQP3 in the shaping of AQP2-mediated drug resistance, I think a figure showing a comparison between the two protein structures/AlphaFold structures would be beneficial and appropriate.

(3) A few additional figures showing cryo-EM density, from both full maps and half maps, would help validate the data.

(4) Finally, this paper might benefit from including more comparisons with and analysis of data published in Chen et al (doi.org/10.1038/s41467-024-48445-4), which focus on similar objectives. Looking at all the data in aggregate might reveal insights that are not obvious from either paper on their own. For example, melarsoprol binds differently in structures reported in the two respective papers, and this may tell us something about the energy of drug-protein interactions within the pore.

Reviewer #1 (Public review):

The manuscript by Chiu et al describes the modification of the Zwitch strategy to efficiently generate conditional knockouts of zebrafish betapix. They leverage this system to identify a surprising glia-exclusive function of betapix in mediating vascular integrity and angiogenesis. Betapix has been previously associated with vascular integrity and angiogenesis in zebrafish, and betapix function in glia has also been proposed. However, this study identifies glial betapix in vascular stability and angiogenesis for the first time.

The study derives its strength from the modified CRISPR-based Zwitch approach to identify the specific role of glial betapix (and not neuronal, mural, or endothelial). Using RNA-in situ hybridization and analysis of scRNA-Seq data, they also identify delayed maturation of neurons and glia and implicate a reduction in stathmin levels in the glial knockouts in mediating vascular homeostasis and angiogenesis. The study also implicates a betapix-zfhx3/4-vegfa axis in mediating cerebral angiogenesis.

There is both technical (the generation of conditional KOs) and knowledge-related (the exclusive role of glial betapix in vascular stability/angiogenesis) novelty in this work that is going to benefit the community significantly.<br /> While the text is well written, it often elides details of experiments and relies on implicit understanding on the part of the reader. Similarly, the figure legends are laconic and often fail to provide all the relevant details.

Specific comments:

(1) While the evidence from cKO's implicating glial betapix in vascular stability/angiogenesis is exciting, glia-specific rescue of betapix in the global KOs/mutants (like those performed for stathmin) would be necessary to make a water-tight case for glial betapix.

(2) Splice variants of betapix have been shown to have differential roles in haemorrhaging (Liu, 2007). What are the major glial isoforms, and are there specific splice variants in the glial that contribute to the phenotypes described?

(3) Liu et al, 2012 demonstrated reduced proliferation of endothelial cells in bbh fish and linked it to deficits in angiogenesis. Are there proliferation/survival defects in endothelial cells in the glial KOs?

Reviewer #2 (Public review):

Summary:

Using a genetic model of beta-pix conditional trap, the authors are able to regulate the spatio-temporal depletion of beta-pix, a gene with an established role in maintaining vascular integrity (shown elsewhere). This study provides strong in vivo evidence that glial beta-pix is essential to the development of the blood-brain barrier and maintaining vascular integrity. Using genetic and biochemical approaches, the authors show that PAK1 and Stathmins are in the same signaling axis as beta-pix, and act downstream to it, potentially regulating cytoskeletal remodeling and controlling glial migration. How exactly the glial-specific (beta-pix driven-) signaling influences angiogenesis or vascular integrity is not clear.

Strengths:

(1) Developing a conditional gene-trap genetic model which allows for tracking knockin reporter driven by endogenous promoter, plus allowing for knocking down genes. This genetic model enabled the authors to address the relevant scientific questions they were interested in, i.e., a) track expression of beta-pix gene, b) deletion of beta-pix gene in a cell-specific manner.

(2) The study reveals the glial-specific role of beta-pix, which was unknown earlier. This opens up avenues for further research. (For instance, how do such (multiple) cell-specific signaling converge onto endothelial cells which build the central artery and maintain the blood-brain barriers?)

Weaknesses:

Major:

(1) The study clearly establishes a role of beta-pix in glial cells, which regulates the length of the central artery and keeps the hemorrhages under control. Nevertheless, it is not clear how this is accomplished.<br /> a. Is this phenotype (hemorrhage) a result of the direct interaction of glial cells and the adjacent endothelial cells? If direct, is the communication established through junctions or through secreted molecules?<br /> b. The authors do not exclude the possibility that the effects observed on endothelial cells (quantified as length of central artery) could be secondary to the phenotype observed with deletion of glial beta-pix. For instance, can glial beta-pix regulate angiogenic factors secreted by peri-vascular cells, which consequently regulate the length of the central artery or vascular integrity?<br /> c. The pictorial summary of the findings (Figure 7) does not include Zfhx or Vegfa. The data do not provide clarity on how these molecules contribute (directly or indirectly) to endothelial cell integrity. Vegfaa is expressed in the central artery, but the expression of the receptor in these endothelial cells is not shown. Similarly, all other experimental analyses for Zfhx and Vegfa expression were performed in glial cells. More experimental evidence is necessary to show the regulation of angiogenesis (of endothelial cells) by glial beta-pix. Is the Vegfaa receptor present on central arteries, and how does glial depletion of beta-pix affect its expression or response of central artery endothelial cells (both pertaining to angiogenesis and vascular integrity).

(2) Microtubule stabilization via glial beta-pix, claimed in Figure 5M, is unclear. Magnified images for h-betapix OE and h-stmn-1 glial cells are absent. Is this migration regulated by beta-pix through its GEF activity for Cdc42/Rac?

(3) Hemorrhages are caused by compromised vascular integrity, which was not measured (either qualitatively or quantitatively) throughout the manuscript. The authors do measure the length of the central artery in several gene deletion models (2I, 3C. 5F/J, 6G/K), which is indicative of artery growth/ angiogenesis. How (if at all) defects in angiogenesis are an indication of hemorrhage should be explained or established. Do these angiogenic growth defects translate into junctional defects at later developmental timepoints? Formation and maintenance of endothelial cell junctions within the hemorrhaging arteries should be assessed in fish with deleted beta-pix from astrocytes.

(4) More information is required about the quality control steps for 10X sequencing (Figure 4, number of cells, reads, etc.). What steps were taken to validate the data quality? The EC groups, 1 and 2-days post-KO are not visible in 4C. One appreciates that the progenitor group is affected the most 2 days post-KO. But since the effects are expected to be on the endothelial cell group as well (which is shown in in vivo data), an extensive analysis should be done on the EC group (like markers for junctional integrity, angiogenesis, mesenchymal interaction, etc.). Are Stathmins limited to glial cells? Are there indicators for angiogenic responses in endothelial cells?

Reviewer #1 (Public review):

Summary:

This study employs spatial transcriptomics to explore the molecular architecture of the adult mouse adrenal gland and the adjacent adipose tissue. The research aimed to identify zonation-specific genetic markers, elucidate cellular differentiation patterns, and investigate inter- and intra-zone communication within the adrenal gland. The findings support the centripetal differentiation model, highlighting the transition of cell populations across different cortical zones. The study also integrates ligand-receptor interaction analysis to uncover the adrenal gland's role in endocrine and neuroendocrine signaling, particularly in stress response. This high-resolution spatial transcriptomic map provides novel insights into adrenal gland biology and is a resource for further investigations.

Strengths:

The study, using the latest technologies and methods such as Visium CytAssist technology, UMAP & Seurat analysis, Gene Ontology (GO) & KEGG pathway enrichment analysis, Monocle3, and CellChat analysis, performed three-dimensional analysis, which has been challenging to achieve using the two-dimensional transcriptomics that have been commonly used up until now.

The unique gene expression patterns were demonstrated for each adrenal zone. Spatial transcriptomics confirmed unique gene expression patterns for each adrenal zone (ZG, ZF, ZX, medulla). The centripetal differentiation model shows the migration of the progenitor cells from the adrenal capsule towards the inner cortex. Key genetic markers were identified in each adrenal zone and adjacent adipose tissues. In addition, CellChat analysis identified major signaling pathways, including Wnt signaling, Hedgehog signaling, IGF2-IGF2R interactions, and Neuropeptide Y (NPY) signaling in the medulla. All these results offer a valuable dataset for future adrenal biology research, with potential applications in disease modeling and therapeutic target identification.

The results, high-resolution mapping of adrenal gland zonation, validation of the centripetal differentiation model, perspective on cell-cell communication, and potential translational impact on human adrenal gland function and disorders, are quite noble.

Weaknesses:

The reviewer requests that the following issues be addressed in the text:

(1) The study focuses only on adult male mice, which limits insights into developmental and sex-specific differences. What do the authors predict about the gender and age difference?

(2) Despite advanced methodologies, single-cell heterogeneity may not be fully captured, as Visium technology has limited spatial resolution.

(3) While the study suggests that ZX might have a role in androgen synthesis, further functional validation is required.

(4) The study is primarily descriptive, lacking in-depth mechanistic experiments to validate cell-cell communication interactions. It is quite interesting to suggest cell-cell communication, but the authors are still required to provide some evidence to support it.

(5) The data supports the conclusions, particularly in validating the centripetal differentiation model using Monocle3 trajectory analysis. However, functional validation experiments (e.g., gene knockout studies) would strengthen the findings, especially regarding ZX function and ligand-receptor interactions.

Reviewer #2 (Public review):

This study by M. Blatkiewicz et al. seeks to define the spatial gene expression pattern of the adult male mouse adrenal gland using current spatial transcriptomic techniques. They propose new zone-specific gene markers and specific intra- and inter-zonal signaling pathways based on receptor-ligand expression patterns. Their web tool is user-friendly and will be helpful for adrenal scientists. The manuscript is easy to follow, but validation of crucial results of the large dataset is missing. There are also several contradictory results/interpretations, and the opportunity to dissect the sexually dimorphic gene expression pattern and mouse-human interspecies differences is a missed opportunity.

(1) The authors used 10-week-old CD1 male mouse adrenal glands to assess the spatial transcriptomics of the adrenal gland. As they also mentioned, male mice typically lose their zone-X after puberty (around 6-8 weeks of age). However, their analysis in 10-week-old mice suggests that zone-X covers most of the adrenal cortex. As shown in Figure 3A, the dots between the zona glomerulosa and the medulla are mostly positive for zone-X, which would suggest that the zona fasciculata represents a relative minority of the overall adult adrenal cortex. Is this correct? Is the presence of zone-X in sexually mature adult male mice unique to the CD1 strain? Providing histology data in support of this conclusion, using zone-specific markers combined with RNA in situ hybridization or immunofluorescence techniques in the CD1 male adrenal gland, would help to interpret these data further. Given the relatively low resolution of their gene expression profiles, it is possible there is overlap between the zona fasciculata and the zone-X.

(2) The pseudotime trajectory analysis confirms prior reports in the literature showing zonal transdifferentiation but does not provide novel insight. It would be nice to know what gene expression patterns correlate (positively or negatively) based on an unbiased analysis.

(3) The authors suggest that they identified new zonal markers, but it would be nice to see confirmation of some of these markers (e.g., Frmpd4, Oca2, Sphkap for the ZG or Cited1, Nat8f5 for the ZF, etc. ) with in situ or immunofluorescence combined with known markers such as Dab2, Cyp11b2, or Cyp11b1.

(4) The authors mention a gradual transition between the zones. It would be interesting to know whether transition zones exist between the zona glomerulosa and the zona fasciculata or the zona fasciculata and the zone-X.

(5) The authors note using Visium cyst assist, but they do not discuss the advantages of this system compared to other systems. Explanation of the approximate resolution of their analysis (e.g., how many cells were pooled in the wells) would help readers to interpret their data. It would also be nice to compare it to other spatial transcriptomic analyses of human adrenals, given the differences between the zonation of human and mouse adrenals.

(6) Interestingly, CellChat analysis suggests possible communication between the medulla and the zona fasciculata and zona glomerulosa. How do the authors explain the transfer of these molecules from the medulla to the outer zones given centripetal blood flow in the adrenal? Also, how does the fact that Igf2 expression has been shown to be expressed in the capsule (PMID: 22266195) affect the interpretation of their data?

(7) The study misses the opportunity to dissect sexually dimorphic gene expression patterns in the mouse adrenal. For example, the authors could have focused on the role of stem cells between male and female mouse adrenals, which have been reported to differ (PMID: 31104943). In addition, the authors could have focused on the sexually dimorphic zone-X and its regulation by sex hormone signaling.

(8) The capsule is classified as a connective tissue, which may be misleading given its important role as a signaling center in the adrenal. Genes enriched in typical connective tissues do not include many of the genes that seem to define the adrenal capsule. Also, some of the capsule markers appear to be found in the zona glomerulosa. Is this a result of not being able to fully resolve the small layer of zG cells and the even smaller layer of capsular cells? Guided reclustering of the cells based on known markers and separation of capsule and connective tissue might help to present their data on adrenal zonation more clearly.

Reviewer #3 (Public review):

Summary:

In summary, the scientists used Visium spatial transcriptomics technology to create a thorough spatial transcriptomic atlas of the adult male mouse adrenal gland and the adipose tissues that surround it. Their primary goals were to map the cell communication network, determine the differentiation direction of various cell types, and find marker genes for various adrenal zones.

Strengths:

(1) Undoubtedly, one of the biggest strengths of the manuscript is a spatial transcriptomic o mouse adrenal gland tissue, which, to my knowledge, has not been done before.

(2) Comprehensive Zonal Characterization: Seven distinct clusters were identified, corresponding to known anatomical and functional regions (ZG, ZF, ZX, medulla, connective tissue, brown and white adipose tissue), each with robust marker gene sets.

(3) The authors manage to integrate advanced bioinformatical tools such as CellChatDB, Monocle3, and CARD to study the relationship between cell types and differentiation of the tissue.

(4) The authors manage to identify novel marker genes for some adrenal zones.

Weaknesses:

(1) The study focused only on one adult male CD1 IGS mouse, which is a limiting factor for other strains, ages, or females, especially given the sexual dimorphism of the ZX. Although the authors claim that four slices of the adrenal gland have been processed on Visium and sequenced, for "clarity," they show only one, which might bias the results.

(2) Lack of detailed QC analysis of the Visium slide.

(3) The study misses the functional validation of the novel marker genes - this needs to be addressed.

(4) What worries me a lot is the fact that, actually, there might be more than one cell present within a Visium spot, so the only way to define zones is by anatomical observation rather than cellular composition.

(5) In cell chat analysis, the authors show the strength of the interactions, but miss out on the number of interactions.

Conclusions:

The authors' stated goals were mostly accomplished:

By mapping the mouse adrenal gland's molecular landscape, they were able to clearly establish unique molecular signatures for every anatomical zone.

Pseudotime study of the cell progression from the capsule through ZG, ZF, and ZX demonstrates that the data strongly support the centripetal differentiation concept. Conclusions on the functional importance of newly discovered marker genes are conjectural and need additional experimental support.

Nevertheless, several findings are still tentative and will need more experimental support, especially when it comes to the significance of ZX persistence and the functional involvement of recently discovered marker genes.

Reviewer #1 (Public review):

Summary:

This study addresses the issue of rapid skill learning and whether individual sequence elements (here: finger presses) are differentially represented in human MEG data. The authors use a decoding approach to classify individual finger elements, and accomplish an accuracy of around 94%. A relevant finding is that the neural representations of individual finger elements dynamically change over the course of learning. This would be highly relevant for any attempts to develop better brain machine interfaces - one now can decode individual elements within a sequence with high precision, but these representations are not static but develop over the course of learning.

Strengths:

The work follows from a large body of work from the same group on the behavioural and neural foundations of sequence learning. The behavioural task is well established a neatly designed to allow for tracking learning and how individual sequence elements contribute. The inclusion of short offline rest periods between learning epochs has been influential because it has revealed that a lot, if not most of the gains in behaviour (ie speed of finger movements) occur in these so-called micro-offline rest periods.

The authors use a range of new decoding techniques, and exhaustively interrogate their data in different ways, using different decoding approaches. Regardless of the approach, impressively high decoding accuracies are observed, but when using a hybrid approach that combines the MEG data in different ways, the authors observe decoding accuracies of individual sequence elements from the MEG data of up to 94%.

Reviewer #2 (Public review):

Summary:

The current paper consists of two parts. The first part is the rigorous feature optimization of the MEG signal to decode individual finger identity performed in a sequence (4-1-3-2-4; 1~4 corresponds to little~index fingers of the left hand). By optimizing various parameters for the MEG signal, in terms of (i) reconstructed source activity in voxel- and parcel-level resolution and their combination, (ii) frequency bands, and (iii) time window relative to press onset for each finger movement, as well as the choice of decoders, the resultant "hybrid decoder" achieved extremely high decoding accuracy (~95%).

In the second part of the paper, armed with the successful 'hybrid decoder,' the authors asked how neural representation of individual finger movement that is embedded in a sequence, changes during a very early period of skill learning and whether and how such representational change can predict skill learning. They assessed the difference in MEG feature patterns between the first and the last press 4 in sequence 41324 at each training trial and found that the pattern differentiation progressively increased over the course of early learning trials. Additionally, they found that this pattern differentiation specifically occurred during the rest period rather than during the practice trial. With a significant correlation between the trial-by-trial profile of this pattern differentiation and that for accumulation of offline learning, the authors argue that such "contextualization" of finger movement in a sequence (e.g., what-where association) underlies the early improvement of sequential skill. This is an important and timely topic for the field of motor learning and beyond.

Strengths:

The use of temporally rich neural information (MEG signal) has a significant advantage over previous studies testing sequential representations using fMRI. This allowed the authors to examine the earliest period (= the first few minutes of training) of skill learning with finer temporal resolution. Through the optimization of MEG feature extraction, the current study achieved extremely high decoding accuracy (approx. 94%) compared to previous works. The finding of the early "contextualization" of the finger movement in a sequence and its correlation to early (offline) skill improvement is interesting and important. The comparison between "online" and "offline" pattern distance is a neat idea.

Weaknesses:

One potential weakness, in terms of the generality, is that the study assessed the single sequence, the "41324" across all participants. Future confirmation test of using different sequences would be important.

Reviewer #3 (Public review):

Summary:

One goal of this paper is to introduce a new approach for highly accurate decoding of finger movements from human magnetoencephalography data via dimension reduction of a "multi-scale, hybrid" feature space. Following this decoding approach, the authors aim to show that early skill learning involves "contextualization" of the neural coding of individual movements, relative to their position in a sequence of consecutive movements. Furthermore, they aim to show that this "contextualization" develops primarily during short rest periods interspersed with skill training, and correlates with a performance metric which the authors interpret as an indicator of offline learning.

Strengths:

A strength of the paper is the innovative decoding approach, which achieves impressive decoding accuracies via dimension reduction of a "multi-scale, hybrid space". This hybrid-space approach follows the neurobiologically plausible idea of concurrent distribution of neural coding across local circuits as well as large-scale networks.

Weaknesses:

A clear weakness of the paper lies in the authors' conclusions regarding "contextualization". Several potential confounds, which partly arise from the experimental design, and which are described below, question the neurobiological implications proposed by the authors, and offer a simpler explanation of the results. Furthermore, the paper follows the assumption that short breaks result in offline skill learning, while recent evidence casts doubt on this assumption.

Specifically:

The authors interpret the ordinal position information captured by their decoding approach as a reflection of neural coding dedicated to the local context of a movement (Figure 4). One way to dissociate ordinal position information from information about the moving effectors is to train a classifier on one sequence, and test the classifier on other sequences that require the same movements, but in different positions (Kornysheva et al., Neuron 2019). In the present study, however, participants trained to repeat a single sequence (4-1-3-2-4). As a result, ordinal position information is potentially confounded by the fixed finger transitions around each of the two critical positions (first and fifth press). Across consecutive correct sequences, the first keypress in a given sequence was always preceded by a movement of the index finger (=last movement of the preceding sequence), and followed by a little finger movement. The last keypress, on the other hand, was always preceded by a ring finger movement, and followed by an index finger movement (=first movement of the next sequence). Figure 3 - supplement 5 shows that finger identity can be decoded with high accuracy (>70%) across a large time window around the time of the keypress, up to at least {plus minus}100 ms (and likely beyond, given that decoding accuracy is still high at the boundaries of the window depicted in that figure). This time window approaches the keypress transition times in this study. Given that distinct finger transitions characterized the first and fifth keypress, the classifier could thus rely on persistent (or "lingering") information from the preceding finger movement, and/or "preparatory" information about the subsequent finger movement, in order to dissociate the first and fifth keypress. Currently, the manuscript provides little evidence that the context information captured by the decoding approach is more than a by-product of temporally extended, and therefore overlapping, but independent neural representations of consecutive keypresses that are executed in close temporal proximity - rather than a neural representation dedicated to context.

During the review process, the authors pointed out that a "mixing" of temporally overlapping information from consecutive keypresses, as described above, should result in systematic misclassifications and therefore be detectable in the confusion matrices in Figures 3C and 4B, which indeed do not provide any evidence that consecutive keypresses are systematically confused. However, such absence of evidence (of systematic misclassification) should be interpreted with caution. The authors also reported that there was only a weak relation between inter-press intervals and "online contextualization" (Figure 5 - figure supplement 6), however, their analysis suprisingly includes a keypress transition that is shared between OP1 and OP5 ("4-4"), rather than focusing solely on the two distinctive transitions ("2-4" and "4-1").

Such temporal overlap of consecutive, independent finger representations may also account for the dynamics of "ordinal coding"/"contextualization", i.e., the increase in 2-class decoding accuracy, across Day 1 (Figure 4C). As learning progresses, both tapping speed and the consistency of keypress transition times increase (Figure 1), i.e., consecutive keypresses are closer in time, and more consistently so. As a result, information related to a given keypress is increasingly overlapping in time with information related to the preceding and subsequent keypresses. Furthermore, learning should increase the number of (consecutively) correct sequences, and, thus, the consistency of finger transitions. Therefore, the increase in 2-class decoding accuracy may simply reflect an increasing overlap in time of increasingly consistent information from consecutive keypresses, which allows the classifier to dissociate the first and fifth keypress more reliably as learning progresses, simply based on the characteristic finger transitions associated with each. In other words, given that the physical context of a given keypress changes as learning progresses - keypresses move closer together in time, and are more consistently correct - it seems problematic to conclude that the mental representation of that context changes. During the review process, authors pointed at absence of evidence of a relation between tapping speed and "ordinal coding" (Figure 5 - figure supplement 7). However, a rigorous test of the idea that the mental representation of context changes would require a task design in which the physical context remains constant.

A similar difference in physical context may explain why neural representation distances ("differentiation") differ between rest and practice (Figure 5). The authors define "offline differentiation" by comparing the hybrid space features of the last index finger movement of a trial (ordinal position 5) and the first index finger movement of the next trial (ordinal position 1). However, the latter is not only the first movement in the sequence, but also the very first movement in that trial (at least in trials that started with a correct sequence), i.e., not preceded by any recent movement. In contrast, the last index finger of the last correct sequence in the preceding trial includes the characteristic finger transition from the fourth to the fifth movement. Thus, there is more overlapping information arising from the consistent, neighbouring keypresses for the last index finger movement, compared to the first index finger movement of the next trial. A strong difference (larger neural representation distance) between these two movements is, therefore, not surprising, given the task design, and this difference is also expected to increase with learning, given the increase in tapping speed, and the consequent stronger overlap in representations for consecutive keypresses.

A further complication in interpreting the results stems from the visual feedback that participants received during the task. Each keypress generated an asterisk shown above the string on the screen. It is not clear why the authors introduced this complicating visual feedback in their task, besides consistency with their previous studies. The resulting systematic link between the pattern of visual stimulation (the number of asterisks on the screen) and the ordinal position of a keypress makes the interpretation of "contextual information" that differentiates between ordinal positions difficult. While the authors report the surprising finding that their eye-tracking data could not predict asterisk position on the task display above chance level, the mean gaze position seemed to vary systematically as a function of ordinal position of a movement - see Figure 4 - figure supplement 3.

The authors report a significant correlation between "offline differentiation" and cumulative micro-offline gains. However, to reach the conclusion that "the degree of representational differentiation -particularly prominent over rest intervals - correlated with skill gains.", the critical question is rather whether "offline differentiation" correlates with micro-offline gains (not with cumulative micro-offline gains). That is, does the degree to which representations differentiate "during" a given rest period correlate with the degree to which performance improves from before to after the same rest period (not: does "offline differentiation" in a given rest period correlate with the degree to which performance has improved "during" all rest periods up to the current rest period - but this is what Figure 5 - figure supplements 1 and 4 show).

The authors follow the assumption that micro-offline gains reflect offline learning. However, there is no compelling evidence in the literature, and no evidence in the present manuscript, that micro-offline gains (during any training phase) reflect offline learning. Instead, emerging evidence in the literature indicates that they do not (Das et al., bioRxiv 2024), and instead reflect transient performance benefits when participants train with breaks, compared to participants who train without breaks, however, these benefits vanish within seconds after training if both groups of participants perform under comparable conditions (Das et al., bioRxiv 2024). During the review process, the authors argued that differences in the design between Das et al. (2024) on the one hand (Experiments 1 and 2), and the study by Bönstrup et al. (2019) on the other hand, may have prevented Das et al. (2024) from finding the assumed (lasting) learning benefit by micro-offline consolidation. However, the Supplementary Material of Das et al. (2024) includes an experiment (Experiment S1) whose design closely follows the early learning phase of Bönstrup et al. (2019), and which, nevertheless, demonstrates that there is no lasting benefit of taking breaks for the acquired skill level, despite the presence of micro-offline gains.

Along these lines, the authors argue that their practice schedule "minimizes reactive inhibition effects", in particular their short practice periods of 10 seconds each. However, 10 seconds are sufficient to result in motor slowing, as report in Bächinger et al., elife 2019, or Rodrigues et al., Exp Brain Res 2009.

An important conceptual problem with the current study is that the authors conclude that performance improves, and representation manifolds differentiate, "during" rest periods. However, micro-offline gains (as well as offline contextualization) are computed from data obtained during practice, not rest, and may, thus, just as well reflect a change that occurs "online", e.g., at the very onset of practice (like pre-planning) or throughout practice (like fatigue, or reactive inhibition).

The authors' conclusion that "low-frequency oscillations (LFOs) result in higher decoding accuracy compared to other narrow-band activity" should be taken with caution, given that the critical decoding analysis for this conclusion was based on data averaged across a time window of 200 ms (Figure 2), essentially smoothing out higher frequency components.

Reviewer #1 (Public review):

Summary:

Diarrheal diseases represent an important public health issue. Among the many pathogens that contribute to this problem, Salmonella enterica serovar Typhimurium is an important one. Due to the rise in antimicrobial resistance and the problems associated with widespread antibiotic use, the discovery and development of new strategies to combat bacterial infections is urgently needed. The microbiome field is constantly providing us with various health-related properties elicited by the commensals that inhabit their mammalian hosts. Harnessing the potential of these commensals for knowledge about host-microbe interactions as well as useful properties with therapeutic implications will likely to remain a fruitful field for decades to come. In this manuscript, Wang et al use various methods, encompassing classic microbiology, genomics, chemical biology, and immunology, to identify a potent probiotic strain that protects nematode and murine hosts from S. enterica infection. Additionally, authors identify gut metabolites that are correlated with protection, and show that a single metabolite can recapitulate the effects of probiotic administration.

Strengths:

The utilization of varied methods by the authors, together with the impressive amount of data generated, to support the claims and conclusions made in the manuscript is a major strength of the work. Also, the ability the move beyond simple identification of the active probiotic, also identifying compounds that are at least partially responsible for the protective effects, is commendable.

Weaknesses:

No major weaknesses noted.

Reviewer #2 (Public review):

Summary:

In this work, the investigators isolated one Lacticaseibacillus rhamnosus strain (P118), and determined this strain worked well against Salmonella Typhimurium infection. Then, further studies were performed to identify the mechanism of bacterial resistance, and a list of confirmatory assays were carried out to test the hypothesis.

Strengths:

The authors provided details regarding all assays performed in this work, and this reviewer trusted that the conclusion in this manuscript is solid. I appreciate the efforts of the authors to perform different types of in vivo and in vitro studies to confirm the hypothesis.

Reviewer #1 (Public review):

Summary:

Walton et al. set out to isolate new phages targeting the opportunistic pathogen Pseudomonas aeruginosa. Using a double ∆fliF ∆pilA mutant strain, they were able to isolate 4 new phages, CLEW-1. -3, -6 and -10, that were unable to infected the parental PAO1F Wt strain. Further experiments showed that the 4 phages were only able to infect a ∆fliF strain, indicating a role of the MS-protein in the flagellum complex. Through further mutational analysis of the flagellum apparatus, the authors were able to identify the involvement of c-di-GMP in phage infection. Depletion of c-di-GMP levels by an inducible phosphodiesterase render the bacteria resistant to phage infection, while elevation of c-di-GMP through the Wsp system made the cells sensitive to infection by CLEW-1. Using TnSeq, the authors were able to not only reaffirm the involvement of c-di-GMP in phage infection but also able to identify the exopolysaccharide PSL as a downstream target for CLEW-1. C-di-GMP is a known regulator of PSL biosynthesis. The authors show that CLEW-1 binds directly to PSL on the cell surface and that deletion of the pslC gene resulted in complete phage resistance. The authors also provide evidence that the phage - PSL interaction happens during the biofilm mode of growth and that the addition of the CLEW-1 phage specifically resulted in a significant loss of biofilm biomass. Lastly, the authors set out to test if CLEW-1 could be used to resolve a biofilm infection using a mouse keratitis model. Unfortunately, while the authors noted a reduction in bacterial load assessed by GFP fluorescence, the keratitis did not resolve under the tested parameters.

Strengths:

The experiments carried out in this manuscript are thoughtful and rational, and sufficient explanation is provided for why the authors chose each specific set of experiments. The data presented strongly supports their conclusions and they give present compelling explanations for any deviation. The authors have not only developed a new technique for screening for phages targeting P. aeruginosa, but also highlights the importance of looking for phages during the biofilm mode of growth, as opposed to the more standard techniques involving planktonic cultures.

Weaknesses:

The authors did not include host-range testing or resistance development in this study, which would have strengthened the paper. Additionally, further characterisation of the CLEW-1 interaction with PSL at the molecular level would also have been welcomed. However, this will likely be the subject of future studies.

Reviewer #2 (Public review):

This manuscript by Walton et al. suggests that they have identified a new bacteriophage that uses the exopolysaccharide Psl from Pseudomonas aeruginosa (PA) as a receptor. As Psl is an important component in biofilms, the authors suggest that this phage (and others similarly isolated) may be able to specifically target biofilm-growing bacteria.

Comments on revised version:

The authors have generally responded well to the reviewers' comments. This has served to improve this manuscript that has identified a new bacteriophage that uses the exopolysaccharide Psl from Pseudomonas aeruginosa as a receptor.

Reviewer #1 (Public review):

Summary:

In this paper, Behruznia and colleagues use long-read sequencing data for 339 strains of the Mycobacterium tuberculosis complex to study genome evolution in this clonal bacterial pathogen. They use both a "classical" pangenome approach that looks at the presence and absence of genes, and a pangenome graph based on whole genomes in order to investigate structural variants in non-coding regions. The comparison of the two approaches is informative and shows that much is missed when focussing only on genes. The two main biological results of the study are that 1) the MTBC has a small pangenome with few accessory genes, and that 2) pangenome evolution is driven by genome reduction. In the revised article, the description of the data set and the methods is much improved, and the comparison of the two pangenome approaches is more consistent. I still think, however, that the discussion of genome reduction suffers from a basic flaw, namely the failure to distinguish clearly between orthologs and homologs/paralogs.

Strengths:

The authors put together the so-far largest data set of long-read assemblies representing most lineages of the Mycobacterium tuberculosis context, and covering a large geographic area. They sequenced and assembled genomes for strains of M. pinnipedi, L9, and La2, for which no high-quality assemblies were available previously. State-of-the-art methods are used to analyze gene presence-absence polymorphisms (Panaroo) and to construct a pangenome graph (PanGraph). Additional analysis steps are performed to address known problems with misannotated or misassembled genes.

Weaknesses:

The revised manuscript has gained much clarity and consistency. One previous criticism, however, has in my opinion not been properly addressed. I think the problem boils down to not clearly distinguishing between orthologs and paralogs/homologs. As this problem affects a main conclusion - the prevalence of deletions over insertions in the MTBC - it should be addressed, if not through additional analyses, then at least in the discussion.

Insertions and deletions are now distinguished in the following way: "Accessory regions were further classified as a deletion if present in over 50% of the 192 sub-lineages or an insertion/duplication if present in less than 50% of sub-lineages." The outcome of this classification is suspicious: not a single accessory region was classified as an insertion/duplication. As a check of sanity, I'd expect at least some insertions of IS6110 to show up, which has produced lineage- or sublineage-specific insertions (Roychowdhury et al. 2015, Shitikov et al. 2019). Why, for example, wouldn't IS6110 insertions in the single L8 strain show up here?

In a fully clonal organism, any insertion/duplication will be an insertion/duplication of an existing sequence, and thus produce a paralog. If I'm correctly understanding your methods section, paralogs are systematically excluded in the pangraph analysis. Genomic blocks are summarized at the sublineage levels as follows (l.184 ): "The DNA sequences from genomic blocks present in at least one sub-lineage but completely absent in others were extracted to look for long-term evolution patterns in the pangenome." I presume this is done using blastn, as in other steps of the analysis.

So a sublineage-specific copy of IS6110 would be excluded here, because IS6110 is present somewhere in the genome in all sublineages. However, the appropriate category of comparison, at least for the discussion of genome reduction, is orthology rather than homology: is the same, orthologous copy of IS6110, at the same position in the genome, present or absent in other sublineages? The same considerations apply to potential sublineage-specific duplicates of PE, PPE, and Esx genes. These gene families play important roles in host-pathogen interactions, so I'd argue that the neglect of paralogs is not a finicky detail, but could be of broader biological relevance.

Reviewer #2 (Public review):

Summary:

The authors attempted to investigate the pangenome of MTBC by using a selection of state-of-the-art bioinformatic tools to analyse 324 complete and 11 new genomes representing all known lineages and sublineages. The aim of their work was to describe the total diversity of the MTBC and to investigate the driving evolutionary force. By using long read and hybrid approaches for genome assembly, an important attempt was made to understand why the MTBC pangenome size was reported to vary in size by previous reports. This study provides strong evidence that the MTBC pangenome is closed and that genome reduction is the main driver of this species evolution.

Strengths:

A stand-out feature of this work is the inclusion of non-coding regions as opposed to only coding regions which was a focus of previous papers and analyses which investigated the MTBC pangenome. A unique feature of this work is that it highlights sublineage-specific regions of difference (RDs) that was previously unknown. Another major strength is the utilisation of long-read whole genomes sequences, in combination with short-read sequences when available. It is known that using only short reads for genome assembly has several pitfalls. The parallel approach of utilizing both Panaroo and Pangraph for pangenomic reconstruction illuminated limitations of both tools while highlighting genomic features identified by both. This is important for any future work and perhaps alludes to the need for more MTBC-specific tools to be developed. Lastly, ample statistical support in the form of Heaps law and genome fluidity calculations for each pangenome to demonstrate that they are indeed closed.

Weaknesses:

There are no major weaknesses in the revised version of this manuscript.

Reviewer #2 (Public review):

Summary:

The authors aimed to explore the role of climbing fibers (CFs) in cerebellar learning, with a focus on optokinetic reflex (OKR) adaptation. Their goal was to understand how CF activity influences memory acquisition, memory consolidation, and memory retrieval by optogenetically suppressing CF inputs at various stages of the learning process.

Strengths:

The study addresses a significant question in the cerebellar field by focusing on the specific role of CFs in adaptive learning. The authors use optogenetic tools to manipulate CF activity. This provides a direct method to test the causal relationship between CF activity and learning outcomes.

Weaknesses:

Despite shedding light on the potential role of CFs in cerebellar learning, the study is hampered by significant methodological issues that question the validity of its conclusions. The absence of detailed evidence on the effectiveness of CF suppression and concerns over tissue damage from optogenetic stimulation weaken the argument that CFs are not essential for memory consolidation. These challenges make it difficult to confirm whether the study's objectives were fully met or if the findings conclusively support the authors' claims. The research commendably attempts to unravel the temporal involvement of CFs in learning but also underscores the difficulties in pinpointing specific neural mechanisms that underlie the phases of learning. Addressing these methodological issues, investigating other signals that might instruct consolidation, and understanding CFs' broader impact on various learning behaviors are crucial steps for future studies.

Comments on revisions:

In this revision, the authors provide new data regarding the effect of eNpHR on CF-evoked complex spiking in vivo but fails to address overall concerns showing the functional effect that explains their causal results. Additionally, the paper has a narrow "CF-or-nothing" framing that leaves unanswered the central question of which signal instructs consolidation if CFs do not. Substantial new experiments and tighter logic are required before the work can serve as a definitive test of CF involvement in different memory processes.

Reviewer #3 (Public review):

Summary:

The authors attempted to study connections with the inferior olive to the cerebellar cortex and analyze impacts on optokinetic reflex using optogenetics to perturb the pathway. This is a commendable effort as these methods are very challenging due to the location of the inferior olive and recording methods.

Strengths:

The authors have shown that climbing fiber activity was altered due to the optogenetic perturbation. They have added an additional figure to show that complex spikes disappear with inhibitory optogenetics and the impacts on behavior are interesting.

Weaknesses:

The images provided to show injection region are difficult to see and specific cell types are not co-labeled. The data and strength of the results would benefit from high-resolution images demonstrating selectivity and expression, in particular for Figure 2A and 3A. In addition, while the processed recording data looks very striking, including the raw data, as done in Figure 2, would again support the conclusions.

One major concern is that the viruses chosen are non-specific to the cell targets and a cre-based approach is lacking to draw conclusions on only the targeted pathway of interest. It is unclear based on the figures provided if the AAVs labeled only the pathway of interest. It would be interesting to know if typical memory acquisition returns in the same animals if inhibition stops and if animal movement was impacted by the perturbation.

Reviewer #1 (Public review):

This paper describes technically-impressive measurements of calcium signals near synaptic ribbons in goldfish bipolar cells. The data presented provides high spatial and temporal resolution information about calcium concentrations along the ribbon at various distances from the site of entry at the plasma membrane. This is important information. Important gaps in the data presented mean that the evidence for the main conclusions is currently inadequate.

Strengths

• The technical aspects of the measurements are impressive. The authors use calcium indicators bound to the ribbon and high speed line scans to resolve changes with a spatial resolution of ~250 nm and temporal resolution of less than 10 ms. These spatial and temporal scales are much closer to those relevant for vesicle release than previous measurements.

• The use of calcium indicators with very different affinities and of different intracellular calcium buffers helps provide confirmation of key results.

Weaknesses

• Multiple key points of the paper lack a statistical test or summary data from populations of cells. For example, the text states that the proximal and distal calcium kinetics in Figure 2A differ. This is not clear from the inset to Figure 2A - where the traces look like scaled versions of each other. Values for time to half-maximal peak fluorescence are given for one example cell but no statistics or summary are provided. Figure 8 shows examples from one cell with no summary data. This issue comes up in other places as well.

• The rise time measurements in Figure 2 are very different for low and high affinity indicators, but no explanation is given for this difference. Similarly, the measurements of peak calcium concentration in Figure 4 are very different with the two indicators. That might suggest that the high affinity indicator is strongly saturated, which raises concerns about whether that is impacting the kinetic measurements.

Reviewer #2 (Public review):

Summary:

The study introduces new tools for measuring intracellular Ca2+ concentration gradients around retinal rod bipolar cell (rbc) synaptic ribbons. This is done by comparing the Ca2+ profiles measured with mobile Ca2+ indicator dyes versus ribbon-tethered (immobile) Ca2+ indicator dyes. The Ca2+ imaging results provide a straightforward demonstration of Ca2+ gradients around the ribbon and validate their experimental strategy. This experimental work is complemented by a coherent, open-source, computational model that successfully describes changes in Ca2+ domains as a function of Ca2+ buffering. In addition, the authors try to demonstrate that there is heterogeneity among synaptic ribbons within an individual rbc terminal.

Strengths:

The study introduces a new set of tools for estimating Ca2+ concentration gradients at ribbon AZs, and the experimental results are accompanied by an open-source, computational model that nicely describes Ca2+ buffering at the rbc synaptic ribbon. In addition, the dissociated retinal preparation remains a valuable approach for studying ribbon synapses. Lastly, excellent EM.

Comments on revisions:

Specific minor comments:

(1) Rewrite the final sentence of the Abstract. It is difficult to understand.

(2) Add a definition in the Introduction (and revisit in the Discussion) that delineates between micro- and nano-domain. A practical approach would be to round up and round down. If you round up from 0.6 um, then it is microdomain which means ~ 1 um or higher. Likewise, round down from 0.3 um to nanodomain? If you are using confocal, or even STED, the resolution for Ca imaging will be in the 100 to 300 nm range. The point of your study is that your new immobile Ca2-ribbon indicator may actually be operating on a tens of nm scale: nanophysiology. The Results are clearly written in a way that acknowledges this point but maybe make such a "definition" comment in the intro/discussion in order to: 1) demonstrate the power of the new Ca2+ indicator to resolve signals at the base of the ribbon (effectively nano), and 2) (Discussion) to acknowledge that some are achieving nanoscopic resolution (50 to 100nm?) with light microscopy (as you ref'd Neef et al., 2018 Nat Comm).

(3) Suggested reference: Grabner et al. 2022 (Sci Adv, Supp video 13, and Fig S5). Here rod Cav channels are shown to be expressed on both sides the ribbon, at its base, and they are within nanometers from other AZ proteins. This agrees with the conclusions from your imaging work.

(4) In the Discussion, add a little more context to what is known about synaptic transmission in the outer and inner retina.. First, state that the postsynaptic receptors (for example: mGluR6-OnBCs vs KARs-Off-BCs, vs. AMPAR-HCs), and possibly the synaptic cleft (ground squirrel), are known to have a significant impact on signaling in the outer retina. In the inner retina, there are many more unknowns. For example, when I think of the pioneering Palmer JPhysio study, which you sight, I think of NMDAR vs AMPAR, and uncertainty in what type postsynaptic cell was patched (GC or AC....). Once you have informed the reader that the postsynapse is known to have a significant impact on signaling, then promote your experimental work that addresses presynaptic processes: "...the new tool and results allow us to explore release heterogeneity, ribbon by ribbon in dissociated preps, which we eventually plan to use at ribbon synapses within slices......to better understand how the presynapse shapes signaling......".

Reviewer #3 (Public review):

Summary:

In this study, the authors have developed a new Ca indicator conjugated to the peptide, which likely recognizes synaptic ribbons and have measured microdomain Ca near synaptic ribbons at retinal bipolar cells. This interesting approach allows one to measure Ca close to transmitter release sites, which may be relevant for synaptic vesicle fusion and replenishment. Though microdomain Ca at the active zone of ribbon synapses has been measured by Hudspeth and Moser, the new study uses the peptide recognizing synaptic ribbons, potentially measuring the Ca concentration relatively proximal to the release sites.

Strengths:

The study is, in principle, technically well done, and the peptide approach is technically interesting, which allows one to image Ca near the particular protein complexes. The approach is potentially applicable to other types of imaging.

Weaknesses:

Peptides may not be entirely specific, and genetic approach tagging particular active zone proteins with fluorescent Ca indicator proteins may well be more specific. Although the authors are aware of this and the peptide approach is generally used for ribbon synapses, the authors should be aware of this, when interpreting the results.

Reviewer #1 (Public review):

Summary:

In this manuscript, the authors attempt to reconstitute some active zone properties by introducing synaptic ribbon proteins into HEK cells. This "ground-up" approach can be valuable for assessing the necessity of specfic proteins in synaptic function. Here, the authors co-transfect a membrane-targeted bassoon, RBP2, calcium channel subunits and Ribeye to generate what they call "synthetic ribbons". The resultant structures show an ability to cluster calcium channels (Figure 4B) and a modest ability to concentrate calcium entry locations (figure 7J). At the light level, the ribeye aggregates look spherical and localize to the membrane through its interaction with the membrane-targeted bassoon and at the EM level the structures resemble those observed when Ribeye is overexpressed alone. It is a nice proof-of-principle in establishing a useful experimental system for studying calcium channel localization and with expression of other proteins perhaps a means to understanding structure and function of the ribbon. The paper does establish that previously described protein-interactions can be reconstituted in a heterologous system to and that the addition of Ribeye can increase the size of calcium channel patches via indirect interactions.

Strengths:

(1) The authors establish a new experimental system for the study of calcium channel localization to active zones.<br /> (2) The clustering of calcium channels to bassoon via RBP2 is a nice confirmation of a previously-described interaction between bassoon and calcium channels in a cell-based system<br /> (3) The "ground-up" approach is an attractive one and theoretically allows one learn a lot about the essential interactions for building a ribbon structure.<br /> (4) The finding that introducing Ribeye can enhance the size of calcium channel patches is a novel finding that is interesting.

Weaknesses:

(1) The addition of EM is welcome, but the structures seem to resemble those created by overexpression of Ribeye alone, albeit at the membrane. It is unclear to me whether the interaction with Bsn or indirect interactions with other proteins has any effect on these structures. Also, while the abstract mentions that the size and shape are similar to ribbons, the EM seems to show that the size and shape are quite variable.<br /> (2) The clustering of channels is accomplished by taking advantage of previously described interactions between RBP2, Ca channels and bassoon. While it is nice to see that it can be reconstituted in a naive cell, the interactions were previously described. The localization of Ribeye to bassoon takes advantage of a previously described interaction between the two and the membrane localization of the complexes required introduction of a membrane-anchoring motif. These factors limit the novelty of the findings.<br /> (3) The difference in Ca imaging between SyRibbons and other locations is subtle. While there are reasonable explanations for why this could be the case, it may limit the utility of this system for studying Ca-channel-ribbon dynamics moving forward.

Reviewer #2 (Public review):

Summary:

The authors show that co-expression of bassoon, RIBEYE, Cav1.3-alpha1, Cav-beta3, Cav-alpha2delta1, and RBP2 in a heterologus system (HEK293 cells) is sufficient to generate a protein complex resembling a presyanptic ribbon-type active zone both in morphology and in function (in clustering voltage-gated Ca channels and creating sites for localized Ca2+ entry). If the 3 separate Cav gene products are taken as a single protein (i.e. a Ca channel), the conclusion is that the core of a ribbon synapse comprises 4 proteins: bassoon holds the RIBEYE-containing ribbon to the plasma membrane, and RPB2 binds to bassoon and Ca channels, tethering the Ca channels to the presynaptic active zone.

Strengths:

(1) Good use of a heterologous system with generally appropriate controls provides convincing evidence that a presynaptic ribbon-type active zone (without the ability to support exocytosis), with the ability to support localized Ca2+ entry (a key feature of ribbon-type pre-synapses) can be assembled from a few proteins.<br /> (2) In the revised manuscript, the authors do a good job of addressing the limitations of their cultured cell-system.

Weaknesses:

(1) Relies on over-expression, which almost certainly diminishes the experimentally-measured parameters (e.g. pre-synapse clustering, localization of Ca2+ entry).<br /> (2) Are HEK cells the best model? HEK cells secrete substances and have a studied-endocytitic pathway, but they do not create neurosecretory vesicles. Initially, I asked why didn't the authors did not try to reconstitute a ribbon synapse in a cell that makes neurosecretory vesicles like a PC12 cell, and the authors addressed this question in their revision.<br /> (3) Related to 1 and 2: the Ca channel localization observed is significant but not so striking given the presence of Cav protein and measurements of Ca2+ influx distributed across the membrane. Presumably, this is the result of overexpression and an absence of pathways for pre-synaptic targeting of Ca channels. But, still, it was surprising that Ca channel localization was so diffuse. I suppose that the authors tried to reduce the effect of over-expression by using an inducible Cav1.3? Even so, the accessory subunits were constitutively over-expressed.

Reviewer #3 (Public review):

Summary:

Ribbon synapses are complex molecular assemblies responsible for synaptic vesicle trafficking in sensory cells of the eye and the inner ear. The Ca2+-dependent exocytosis occurs at the active zone (AZ), however, the molecular mechanisms orchestrating the structure and function of the AZs of ribbon synapses are not well understood. To advance in the understanding of those mechanisms, the authors present a novel and interesting experimental strategy pursuing the reconstitution of a minimal active zone of a ribbon synapse within a synapse-naïve cell line: HEK293 cells. The authors have used stably transfected HEK293 cells that express voltage-gated Ca2+ channels subunits (constitutive -CaV beta3 and CaV alpha2 beta1- and inducible CaV1.3 alpha1). They have expressed in those cells several proteins of the ribbon synapse active zone: (1) RIBEYE, (2) a modified version of Bassoon that binds to the plasma membrane through artificial palmitoylation (Palm-Bassoon) and (3) RIM-binding protein 2 (RBP2) to induce the formation of a minimal active zone that they called SyRibbons. The formation of such structures is convincing, however, the evidence of such structures having a functional impact (for example enhancing Ca2+-currents), as the authors claim, is weak. In conclusion, the novel approach shows that expression of a multiprotein complex partially reproduces properties, especially structural properties, of ribbon-type active zones in a heterologous system. Although the approach opens interesting possibilities for further experiments, the evidence supporting the functional properties of the so called "synthetic ribbon synapses" is incomplete.

Strengths of the study:

(1) The study is carefully carried out using a remarkable combination of (1) superresolution, correlative light microscopy and cryo-electron tomography, to analyze the formation and subcellular distribution of molecular assemblies and (2) functional assessment of voltage-gated Ca2+ channels using patch-clamp recording of Ca2+-currents and fluorometry to correlate Ca2+ influx with the molecular assemblies formed by AZ proteins. The results are of high quality and are in general accompanied of required control experiments.<br /> (2) The method opens new opportunities to further investigate the minimal and basic properties of AZ proteins that are difficult to study using in vivo systems. The cells that operate through ribbon synapses (e.g. photoreceptors and hair cells) are particularly difficult to manipulate, so setting up and validating the use of a heterologous system more suitable for molecular manipulations is highly valuable.<br /> (3) The structures formed by RIBEYE and Palm-Bassoon in HEK293 cells identified by STED nanoscopy and cryo-electron microscopy share relevant similarities similar to the AZs of ribbon synapses found in rat inner hair cells.

Weaknesses of the study:

(1) The evidence of the functional properties of the "synthetic ribbon-type active zones" has been only assessed by its effect on the modulation of Ca2+-channel function, and that effect is rather weak. The authors provide reasonable explanations regarding such a weak effect but, however, it is difficult to conclude that indeed the "synthetic ribbon-type active zones" are bona fide functional multiprotein complexes.

Reviewer #1 (Public review):

Summary:

The authors, Dalal, et. al., determined cryo-EM structures of open, closed, and desensitized states of the pentameric ligand-gated ion channel ELIC reconstituted in liposomes, and compared them to structures determined in varying nanodisc diameters. They argue that the liposomal reconstitution method is more representative of functional ELIC channels, as they were able to test and recapitulate channel kinetics through stopped-flow thallium flux liposomal assay. The authors and others have described channel interactions with membrane scaffold proteins (MSP), initially thought to be in a size-dependent manner. However, the authors reported their cryo-EM ELIC structure interacts with the large nanodisc spNW25, contrary to their original hypotheses. This suggests that the channels interactions with MSPs might alter its structure, possibly influencing the functional states of the channel. Thus, the authors describe reconstitution in liposomes are more representative of the native structure and can recapitulate all channel states.

Strengths:

Cryo-EM structural determination from proteoliposomes is promising methodology within the ion channel field due to their large surface area and lack of MSP or other membrane memetics that could alter channel structure. The authors succeeded in comparing structures determined in liposomes to those in a wide range of nanodisc diameters. This comparison gives rise to important discussions for other membrane protein structural studies when deciding the best method for individual circumstances.

Weaknesses:

As the overarching goal of the study was to determine structural differences of ELIC in detergent nanodiscs and liposomes. The authors stated they determined open, closed, and desensitized states of ELIC reconstituted in liposomes and suggest the desensitization gate is at the 9' region of the pore. However, limited functional data was provided when determining the functional states of the channel with most of the evidence deriving from structures, which only provides snapshots of channels.

Reviewer #2 (Public review):

Summary

The report by Dalas and colleagues introduces a significant novelty in the field of pentameric ligand-gated ion channels (pLGICs). Within this family of receptors, numerous structures are available, but a widely recognised problem remains in assigning structures to functional states observed in biological membranes. Here, the authors obtain both structural and functional information of a pLGIC in a liposome environment. The model receptor ELIC is captured in the resting, desensitised and open states. Structures in large nanodiscs, possibly biased by receptor-scaffold protein interactions, are also reported. Altogether these results set the stage for the adoption of liposomes as a proxy for the biological membranes, for cryoEM studies of pLGICs and membrane proteins in general.

Strengths

The structural data is comprehensive, with structures in liposomes in the 3 main states (and for each, both inward-facing and outward-facing), and an agonist-bound structure in the large spNW25 nanodisc (and a retreatment of previous data obtained in a smaller disc). It adds up to a series of work from the same team that constitutes a much-needed exploration of various types of environment for the transmembrane domain of pLGICs. The structural analysis is thorough.<br /> The tone of the report is particularly pleasant, in the sense that the authors' claims are not inflated. For instance, a sentence such as "By performing structural and functional characterization under the same reconstitution conditions, we increase our confidence in the functional annotation of these structures." is exemplary.

Weakness

All the details necessary to reproduce the work are present in the Methods. Nevertheless, the biochemistry might have been shown and discussed in greater details. While I do believe that liposomes will be in most cases better than, say, nanodiscs, the process that leads from the protein in its membrane down to the liposome will play a big role in preserving the native structure.

Reviewer #1 (Public review):

Summary:

The authors investigated the elasticity of controllability by developing a task that manipulates the probability of achieving a goal with a baseline investment (which they refer to as inelastic controllability) and the probability that additional investment would increase the probability of achieving a goal (which they refer to as elastic controllability). They found that a computational model representing the controllability and elasticity of the environment accounted better for the data than a model representing only the controllability. They also found that prior biases about the controllability and elasticity of the environment was associated with a composite psychopathology score. The authors conclude that elasticity inference and bias guide resource allocation.

Strengths:

This research takes a novel theoretical and methodological approach to understanding how people estimate the level of control they have over their environment and how they adjust their actions accordingly. The task is innovative and both it and the findings are well-described (with excellent visuals). They also offer thorough validation for the particular model they develop. The research has the potential to theoretically inform understanding of control across domains, which is a topic of great importance.

Weaknesses:

In its revised form, the manuscript addresses most of my previous concerns. The main remaining weakness pertains to the analyses aimed at addressing my suggesting of Bayesian updating as an alternative to the model proposed by the authors. My suggestion was to assume that people perform a form of function approximation to relate resource expenditure to success probability. The authors performed a version of this where people were weighing evidence for a few canonical functions (flat, step, linear), and found that this model underperformed theirs. However, this Bayesian model is quite constrained in its ability to estimate the function relating resources. A more robust test would be to assume a more flexible form of updating that is able to capture a wide range of distributions (e.g., using basis functions, gaussian processes, or nonparametric estimators); see, e.g., work by Griffiths on human function learning). The benefit of testing this type of model is that it would make contact with a known form of inference that individuals engage in across various settings and therefore could offer a more parsimonious and generalizable account of function learning, whereby learning of resource elasticity is a special case. I defer to the authors as to whether they'd like to pursue this direction, but if not I think it's still important that they acknowledge that they are unable to rule out a more general process like this as an alternative to their model. This pertains also to inferences about individual differences, which currently hinge on their preferred model being the most parsimonious.

Reviewer #2 (Public review):

Summary:

In this paper, the authors test whether controllability beliefs and associated actions/resource allocation are modulated by things like time, effort, and monetary costs (what they call "elastic" as opposed to "inelastic" controllability). Using a novel behavioral task and computational modeling, they find that participants do indeed modulate their resources depending on whether they are in an "elastic," "inelastic," or "low controllability" environment. The authors also find evidence that psychopathology is related to specific biases in controllability.

Strengths:

This research investigates how people might value different factors that contribute to controllability in a creative and thorough way. The authors use computational modeling to try to dissociate "elasticity" from "overall controllability," and find some differential associations with psychopathology. This was a convincing justification for using modeling above and beyond behavioral output and yielded interesting results. Notably, the authors conclude that these findings suggest that biased elasticity could distort agency beliefs via maladaptive resource allocation. Overall, this paper reveals important findings about how people consider components of controllability.

Weaknesses:

The authors have gone to great lengths to revise the manuscript to clarify their definitions of "elastic" and "inelastic" and bolster evidence for their computational model, resulting in an overall strong manuscript that is valuable for elucidating controllability dynamics and preferences. One minor weakness is that the justification for the analysis technique for the relationships between the model parameters and the psychopathology measures remains lacking given the fact that simple correlational analyses did not reveal any significant associations.

Reviewer #3 (Public review):

A bias in how people infer the amount of control they have over their environment is widely believed to be a key component of several mental illnesses including depression, anxiety, and addiction. Accordingly, this bias has been a major focus in computational models of those disorders. However, all of these models treat control as a unidimensional property, roughly, how strongly outcomes depend on action. This paper proposes---correctly, I think---that the intuitive notion of "control" captures multiple dimensions in the relationship between action and outcome. In particular, the authors identify one key dimension: the degree to which outcome depends on how much *effort* we exert, calling this dimension the "elasticity of control". They additionally argue that this dimension (rather than the more holistic notion of controllability) may be specifically impaired in certain types of psychopathology. This idea has the potential to change how we think about several major mental disorders in a substantial way and can additionally help us better understand how healthy people navigate challenging decision-making problems. More concisely, it is a very good idea.

Unfortunately, my view is that neither the theoretical nor empirical aspects of the paper really deliver on that promise. In particular, most (perhaps all) of the interesting claims in the paper have weak empirical support.

Starting with theory, the authors do not provide a strong formal characterization of the proposed notion of elasticity. There are existing, highly general models of controllability (e.g., Huys & Dayan, 2009; Ligneul, 2021) and the elasticity idea could naturally be embedded within one of these frameworks. The authors gesture at this in the introduction; however, this formalization is not reflected in the implemented model, which is highly task-specific. Moreover, the authors present elasticity as if it is somehow "outside of" the more general notion of controllability. However, effort and investment are just specific dimensions of action; and resources like money, strength, and skill (the "highly trained birke") are just specific dimensions of state. Accordingly, the notion of elasticity is necessarily implicitly captured by the standard model. Personally, I am compelled by the idea that effort and resource (and therefore elasticity) are particularly important dimensions, ones that people are uniquely tuned to. However, by framing elasticity as a property that is different in kind from controllability (rather than just a dimension of controllability), the authors only make it more difficult to integrate this exciting idea into generalizable models.

Turning to experiment, the authors make two key claims: (1) people infer the elasticity of control, and (2) individual differences in how people make this inference are importantly related to psychopathology.

Starting with claim 1, there are three subclaims here; implicitly, the authors make all three. (1A) People's behavior is sensitive to differences in elasticity, (1B) people actually represent/track something like elasticity, and (1C) people do so naturally as they go about their daily lives. The results clearly support 1A. However, 1B and 1C are not strongly supported.

(1B) The experiment cannot support the claim that people represent or track elasticity because effort is the only dimension over which participants can engage in any meaningful decision-making. The other dimension, selecting which destination to visit, simply amounts to selecting the location where you were just told the treasure lies. Thus, any adaptive behavior will necessarily come out in a sensitivity to how outcomes depend on effort.

Notes on rebuttal: The argument that vehicle/destination choice is not trivial because people occasionally didn't choose the instructed location is not compelling to me-if anything, the exclusion rate is unusually low for online studies. The finding that people learn more from non-random outcomes is helpful, but this could easily be cast as standard model-based learning very much like what one measures with the Daw two-step task (nothing specific to control here). Their final argument is the strongest, that to explain behavior the model must assume "a priori that increased effort could enhance control." However, more literally, the necessary assumption is that each attempt increases the probability of success-e.g. you're more likely to get a heads in two flips than one. I suppose you can call that "elasticity inference", but I would call it basic probabilistic reasoning.

For 1C, the claim that people infer elasticity outside of the experimental task cannot be supported because the authors explicitly tell people about the two notions of control as part of the training phase: "To reinforce participants' understanding of how elasticity and controllability were manifested in each planet, [participants] were informed of the planet type they had visited after every 15 trips." (line 384).

Notes on rebuttal: The authors try to retreat, saying "our research question was whether people can distinguish between elastic and inelastic controllability." I struggle to reconcile this with the claim in the abstract "These findings establish the elasticity of control as a distinct cognitive construct guiding adaptive behavior". That claim is the interesting one, and the one I am evaluating the evidence in light of.

Finally, I turn to claim 2, that individual differences in how people infer elasticity are importantly related to psychopathology. There is much to say about the decision to treat psychopathology as a unidimensional construct (the authors claim otherwise, but see Fig 6C). However, I will keep it concrete and simply note that CCA (by design) obscures the relationship between any two variables. Thus, as suggestive as Figure 6B is, we cannot conclude that there is a strong relationship between Sense of Agency (SOA) and the elasticity bias---this result is consistent with any possible relationship (even a negative one). As it turns out, Figure S3 shows that there is effectively no relationship (r=0.03).

Notes on rebuttal: The authors argue for CCA by appeal to the need to "account for the substantial variance that is typically shared among different forms of psychopathology". I agree. A simple correlation would indeed be fairly weak evidence. Strong evidence would show a significant correlation after *controlling for* other factors (e.g. a regression predicting elasticity bias from all subscales simultaneously). CCA effectively does the opposite, asking whether-with the help of all the parameters and all the surveys-one can find any correlation between the two sets of variables. The results are certainly suggestive, but they provide very little statistical evidence that the elasticity parameter is meaningfully related to any particular dimension of psychopathology.

There is also a feature of the task that limits our ability to draw strong conclusions about individual differences about elasticity inference. In the original submission, the authors stated that the study was designed to be "especially sensitive to overestimation of elasticity". A straightforward consequence of this is that the resulting *empirical* estimate of estimation bias (i.e., the gamma_elasticity parameter) is itself biased. This immediately undermines any claim that references the directionality of the elasticity bias (e.g. in the abstract). Concretely, an undirected deficit such as slower learning of elasticity would appear as a directed overestimation bias.

When we further consider that elasticity inference is the only meaningful learning/decision-making problem in the task (argued above), the situation becomes much worse. Many general deficits in learning or decision-making would be captured by the elasticity bias parameter. Thus, a conservative interpretation of the results is simply that psychopathology is associated with impaired learning and decision-making.

Notes on rebuttal: I am very concerned to see that the authors removed the discussion of this limitation in response to my first review. I quote the original explanation here:

- In interpreting the present findings, it needs to be noted that we designed our task to be especially sensitive to overestimation of elasticity. We did so by giving participants free 3 tickets at their initial visits to each planet, which meant that upon success with 3 tickets, people who overestimate elasticity were more likely to continue purchasing extra tickets unnecessarily. Following the same logic, had we first had participants experience 1 ticket trips, this could have increased the sensitivity of our task to underestimation of elasticity in elastic environments. Such underestimation could potentially relate to a distinct psychopathological profile that more heavily loads on depressive symptoms. Thus, by altering the initial exposure, future studies could disambiguate the dissociable contributions of overestimating versus underestimating elasticity to different forms of psychopathology.

The logic of this paragraph makes perfect sense to me. If you assume low elasticity, you will infer that you could catch the train with just one ticket. However, when elasticity is in fact high, you would find that you don't catch the train, leading you to quickly infer high elasticity-eliminating the bias. In contrast, if you assume high elasticity, you will continue purchasing three tickets and will never have the opportunity to learn that you could be purchasing only one-the bias remains.

The authors attempt to argue that this isn't happening using parameter recovery. However, they only report the *correlation* in the parameter, whereas the critical measure is the *bias*. Furthermore, in parameter recovery, the data-generating and data-fitting models are identical-this will yield the best possible recovery results. Although finding no bias in this setting would support the claims, it cannot outweigh the logical argument for the bias that they originally laid out. Finally, parameter recovery should be performed across the full range of plausible parameter values; using fitted parameters (a detail I could only determine by reading the code) yields biased results because the fitted parameters are themselves subject to the bias (if present). That is, if true low elasticity is inferred as high elasticity, then you will not have any examples of low elasticity in the fitted parameters and will not detect the inability to recover them.

Minor comments:

Below are things to keep in mind.

The statistical structure of the task is inconsistent with the framing. In the framing, participants can make either one or two second boarding attempts (jumps) by purchasing extra tickets. The additional attempt(s) will thus succeed with probability p for one ticket and 2p - p^2 for two tickets; the p^2 captures the fact that you only take the second attempt if you fail on the first. A consequence of this is buying more tickets has diminishing returns. In contrast, in the task, participants always jumped twice after purchasing two tickets, and the probability of success with two tickets was exactly double that with one ticket. Thus, if participants are applying an intuitive causal model to the task, they will appear to "underestimate" the elasticity of control. I don't think this seriously jeopardizes the key results, but any follow-up work should ensure that the task's structure is consistent with the intuitive causal model.

The model is heuristically defined and does not reflect Bayesian updating. For example, it over-estimates maximum control by not using losses with less than 3 tickets (intuitively, the inference here depends on what your beliefs about elasticity). Including forced three-ticket trials at the beginning of each round makes this less of an issue; but if you want to remove those trials, you might need to adjust the model. The need to introduce the modified model with kappa is likely another symptom of the heuristic nature of the model updating equations.

Reviewer #1 (Public review):

Summary:

This is an interesting study characterizing and engineering so-called bathy phytochromes, i.e., those that respond to near infrared (NIR) light in the ground state, for optogenetic control of bacterial gene expression. Previously, the authors have developed a structure-guided approach to functionally link several light-responsive protein domains to the signaling domain of the histidine kinase FixL, which ultimately controls gene expression. Here, the authors use the same strategy to link bathy phytochrome light-responsive domains to FixL, resulting in sensors of NIR light. Interestingly, they also link these bathy phytochrome light-sensing domains to signaling domains from the tetrathionate-sensing SHK TtrS and the toluene-sensing SHK TodS, demonstrating the generality of their protein engineering approach more broadly across bacterial two-component systems.

This is an exciting result that should inspire future bacterial sensor design. They go on to leverage this result to develop what is, to my knowledge, the first system for orthogonally controlling the expression of two separate genes in the same cell with NIR and Red light, a valuable contribution to the field.

Finally, the authors reveal new details of the pH-dependent photocycle of bathy phytochromes and demonstrate that their sensors work in the gut - and plant-relevant strains E. coli Nissle 1917 and A. tumefaciens.

Strengths:

(1) The experiments are well-founded, well-executed, and rigorous.

(2) The manuscript is clearly written.

(3) The sensors developed exhibit large responses to light, making them valuable tools for ontogenetic applications.

(4) This study is a valuable contribution to photobiology and optogenetics.

Weaknesses:

(1) As the authors note, the sensors are relatively insensitive to NIR light due to the rapid dark reversion process in bathy phytochromes. Though NIR light is generally non-phototoxic, one would expect this characteristic to be a limitation in some downstream applications where light intensities are not high (e.g., in vivo).

(2) Though they can be multiplexed with Red light sensors, these bathy phytochrome NIR sensors are more difficult to multiplex with other commonly used light sensors (e.g., blue) due to the broad light responsivity of the Pfr state. This challenge may be overcome by careful dosing of blue light, as the authors discuss, but other bacterial NIR sensing systems with less cross-talk may be preferred in some applications.

Reviewer #2 (Public review):

Summary:

In this manuscript, Meier et al. engineer a new class of light-regulated two-component systems. These systems are built using bathy-bacteriophytochromes that respond to near-infrared (NIR) light. Through a combination of genetic engineering and systematic linker optimization, the authors generate bacterial strains capable of selective and tunable gene expression in response to NIR stimulation. Overall, these results are an interesting expansion of the optogenetic toolkit into the NIR range. The cross-species functionality of the system, modularity, and orthogonality have the potential to make these tools useful for a range of applications.

Strengths:

(1) The authors introduce a novel class of near-infrared light-responsive two-component systems in bacteria, expanding the optogenetic toolbox into this spectral range.

(2) Through engineering and linker optimization, the authors achieve specific and tunable gene expression, with minimal cross-activation from red light in some cases.

(3) The authors show that the engineered systems function robustly in multiple bacterial strains, including laboratory E. coli, the probiotic E. coli Nissle 1917, and Agrobacterium tumefaciens.

(4) The combination of orthogonal two-component systems can allow for simultaneous and independent control of multiple gene expression pathways using different wavelengths of light.

(5) The authors explore the photophysical properties of the photosensors, investigating how environmental factors such as pH influence light sensitivity.

Weaknesses:

(1) The expression of multi-gene operons and fluorescent reporters could impose a metabolic burden. The authors should present data comparing optical density for growth curves of engineered strains versus the corresponding empty-vector control to provide insight into the burden and overall impact of the system on host viability and growth.

(2) The manuscript consistently presents normalized fluorescence values, but the method of normalization is not clear (Figure 2 caption describes normalizing to the maximal fluorescence, but the maximum fluorescence of what?). The authors should provide a more detailed explanation of how the raw fluorescence data were processed. In addition, or potentially in exchange for the current presentation, the authors should include the raw fluorescence values in supplementary materials to help readers assess the actual magnitude of the reported responses.

(3) Related to the prior point, it would be useful to have a positive control for fluorescence that could be used to compare results across different figure panels.

(4) Real-time gene expression data are not presented in the current manuscript, but it would be helpful to include a time-course for some of the key designs to help readers assess the speed of response to NIR light.

Reviewer #3 (Public review):

Summary:

This paper by Meier et al introduces a new optogenetic module for the regulation of bacterial gene expression based on "bathy-BphP" proteins. Their paper begins with a careful characterization of kinetics and pH dependence of a few family members, followed by extensive engineering to produce infrared-regulated transcriptional systems based on the authors' previous design of the pDusk and pDERusk systems, and closing with characterization of the systems in bacterial species relevant for biotechnology.

Strengths:

The paper is important from the perspective of fundamental protein characterization, since bathy-BphPs are relatively poorly characterized compared to their phytochrome and cyanobacteriochrome cousins. It is also important from a technology development perspective: the optogenetic toolbox currently lacks infrared-stimulated transcriptional systems. Infrared light offers two major advantages: it can be multiplexed with additional tools, and it can penetrate into deep tissues with ease relative to the more widely used blue light-activated systems. The experiments are performed carefully, and the manuscript is well written.

Weaknesses:

My major criticism is that some information is difficult to obtain, and some data is presented with limited interpretation, making it difficult to obtain intuition for why certain responses are observed. For example, the changes in red/infrared responses across different figures and cellular contexts are reported but not rationalized. Extensive experiments with variable linker sequences were performed, but the rationale for linker choices was not clearly explained. These are minor weaknesses in an overall very strong paper.

Reviewer #1 (Public review):

Summary:

Can a plastic RNN serve as a basis function for learning to estimate value. In previous work this was shown to be the case, with a similar architecture to that proposed here. The learning rule in previous work was back-prop with an objective function that was the TD error function (delta) squared. Such a learning rule is non-local as the changes in weights within the RNN, and from inputs to the RNN depends on the weights from the RNN to the output, which estimates value. This is non-local, and in addition, these weights themselves change over learning. The main idea in this paper is to examine if replacing the values of these non-local changing weights, used for credit assignment, with random fixed weights can still produce similar results to those obtained with complete bp. This random feedback approach is motivated by a similar approach used for deep feed-forward neural networks.

This work shows that this random feedback in credit assignment performs well but is not as well as the precise gradient-based approach. When more constraints due to biological plausibility are imposed performance degrades. These results are consistent with previous results on random feedback.

Strengths:

• The authors show that random feedback can approximate well a model trained with detailed credit assignment.<br /> • The authors simulate several experiments including some with probabilistic reward schedules and show results similar to those obtained with detailed credit assignments as well as in experiments.<br /> • The paper examines the impact of more biologically realistic learning rules and the results are still quite similar to the detailed back-prop model.

Weaknesses:

• The impact of the article is limited by using a network with discrete time-steps, and only a small number of time steps from stimulus to reward. They assume that each time step is on the order of hundreds of ms. They justify this by pointing to some slow intrinsic mechanisms, but they do not implement these slow mechanisms is a network with short time steps, instead they assume without demonstration that these could work as suggested. This is a reasonable first approximation, but its validity should be explicitly tested.

• As the delay between cue and reward increases the performance decreases. This is not surprising given the proposed mechanism, but is still a limitation, especially given that we do not really know what a is the reasonable value of a single time step.

Reviewer #2 (Public review):

Summary:

Tsurumi et al. show that recurrent neural networks can learn state and value representations in simple reinforcement learning tasks when trained with random feedback weights. The traditional method of learning for recurrent network in such tasks (backpropogation through time) requires feedback weights which are a transposed copy of the feed-forward weights, a biologically implausible assumption. This manuscript builds on previous work regarding "random feedback alignment" and "value-RNNs", and extends them to a reinforcement learning context. The authors also demonstrate that certain non-negative constraints can enforce a "loose alignment" of feedback weights. The author's results suggest that random feedback may be a powerful tool of learning in biological networks, even in reinforcement learning tasks.

Strengths:

The authors describe well the issues regarding biologically plausible learning in recurrent networks and in reinforcement learning tasks. They take care to propose networks which might be implemented in biological systems and compare their proposed learning rules to those already existing in literature. Further, they use small networks on relatively simple tasks, which allows for easier intuition into the learning dynamics.

Weaknesses:

The principles discovered by the authors in these smaller networks are not applied to larger networks or more complicated tasks with long temporal delays (>100 timesteps), so it remains unclear to what degree these methods can scale or can be used more generally.

Comments on revisions: I would still want to see how well the network learns tasks with longer time delays (on the order of 100 or even 1000 timesteps). Previous work has shown that random feedback struggles to encode longer timescales (see Murray 2019, Figure 2), so I would be interested to see how that translates to the RL context in your model.

Reviewer #3 (Public review):

Summary:

The paper studies learning rules in a simple sigmoidal recurrent neural network setting. The recurrent network has a single layer of 10 to 40 units. It is first confirmed that feedback alignment (FA) can learn a value function in this setting. Then so-called bio-plausible constraints are added: (1) when value weights (readout) is non-negative, (2) when the activity is non-negative (normal sigmoid rather than downscaled between -0.5 and 0.5), (3) when the feedback weights are non-negative, (4) when the learning rule is revised to be monotic: the weights are not downregulated. In the simple task considered all four biological features do not appear to impair totally the learning.

Strengths:

(1) The learning rules are implemented in a low-level fashion of the form: (pre-synaptic-activity) x (post-synaptic-activity) x feedback x RPE. Which is therefore interpretable in terms of measurable quantities in the wet-lab.

(2) I find that non-negative FA (FA with non negative c and w) is the most valuable theoretical insight of this paper: I understand why the alignment between w and c is automatically better at initialization.

(3) The task choice is relevant, since it connects with experimental settings of reward conditioning with possible plasticity measurements.

Weaknesses:

(4) The task is rather easy, so it's not clear that it really captures the computational gap that exists with FA (gradient-like learning) and simpler learning rule like a delta rule: RPE x (pre-synpatic) x (post-synaptic). To control if the task is not too trivial, I suggest adding a control where the vector c is constant c_i=1.

(5) Related to point 3), the main strength of this paper is to draw potential connection with experimental data. It would be good to highlight more concretely the prediction of the theory for experimental findings. (Ideally, what should be observed with non-negative FA that is not expected with FA or a delta rule (constant global feedback) ?).

(6a) Random feedback with RNN in RL have been studied in the past, so it is maybe worth giving some insights how the results and the analyzes compare to this previous line of work (for instance in this paper [1]). For instance, I am not very surprised that FA also works for value prediction with TD error. It is also expected from the literature that the RL + RNN + FA setting would scale to tasks that are more complex than the conditioning problem proposed here, so is there a more specific take-home message about non-negative FA? or benefits from this simpler toy task?

(6b) Related to task complexity, it is not clear to me if non-negative value and feedback weights would generally scale to harder tasks. If the task in so simple that a global RPE signal is sufficient to learn (see 4 and 5), then it could be good to extend the task to find a substantial gap between: global RPE, non-negative FA, FA, BP. For a well chosen task, I expect to see a performance gap between any pair of these four learning rules. In the context of the present paper, this would be particularly interesting to study the failure mode of non-negative FA and the cases where it does perform as well as FA.

(7) I find that the writing could be improved, it mostly feels more technical and difficult than it should. Here are some recommendations:<br /> 7a) For instance, the technical description of the task (CSC) is not fully described and requires background knowledge from other paper which is not desirable.<br /> 7b) Also the rationale for the added difficulty with the stochastic reward and new state is not well explained.<br /> 7c) In the technical description of the results I find that the text dives into descriptive comments of the figures but high-level take home messages would be helpful to guide the reader. I got a bit lost, although I feel that there is probably a lot of depth in these paragraphs.

(8) Related to the writing issue and 5), I wished that "bio-plausibility" was not the only reason to study positive feedback and value weights. Is it possible to develop a bit more specifically what and why this positivity is interesting? Is there an expected finding with non-negative FA both in the model capability? or maybe there is a simpler and crisp take-home message to communicate the experimental predictions to the community would be useful?

[1] https://www.nature.com/articles/s41467-020-17236-y

Comments on revisions:

Thank you for addressing all my comments in your reply.

Reviewer #2 (Public review):

Summary:

In this study, the authors characterized population genetic variation in the MHC locus across primates and looked for signals of long-term balancing selection (specifically trans-species polymorphism, TSP) in this highly polymorphic region. To carry out these tasks, they used Bayesian methods for phylogenetic inference (i.e. BEAST2) and applied a new Bayesian test to quantify evidence supporting monophyly vs. transspecies polymorphism for each exon across different species pairs. Their results, although mostly confirmatory, represent the most comprehensive analyses of primate MHC evolution to date and novel findings or possible discrepancies are clearly pointed out. However, as the authors discuss, the available data are insufficient to fully capture primates' MHC evolution.

Strengths of the paper include: using appropriate methods and statistically rigorous analyses; very clear figures and detailed description of the results methods that make it easy to follow despite the complexity of the region and approach; a clever test for TSP that is then complemented by positive selection tests and the protein structures for a quite comprehensive study.

That said, weaknesses include: lack of information about how many sequences are included and whether uneven sampling across taxa might results in some comparisons without evidence for TSP; frequent reference to the companion paper instead of summarizing (at least some of) the critical relevant information (e.g., how was orthology inferred?); no mention of the quality of sequences in the database and whether there is still potential effects of mismapping or copy number variation affecting the sequence comparison.

Comments on revisions:

The authors have sufficiently addressed the reviewers' comments or provided additional details justifying their work. In particular, expansion of the discussion section on limitations of the analysis and clearer reference to how this relates to their companion paper represent improvements. Remaining suggestions are to still make clearer how much sparsity of sequences in the database may impact the conclusions (e.g., is this more of a problem for some genes or taxa than others? Is it a small problem or a large problem?). The data summary tables are a bit hard to read and seem to contain some information not used in the article - maybe the presentation of these could be improved or the full details, or a shorter table summer in the main paper and full details only in the supplement.

Reviewer #3 (Public review):

Summary:

The study uses publicly available sequences of classical and non-classical genes from a number of primate species to assess the extent and depth of TSP across the primate phylogeny. The analyses were carried out in a coherent and, in my opinion, robust inferential framework and provide evidence for ancient (even > 30 million years) TSP at several classical class I and class II genes. The authors also characterise evolutionary rates at individual codons, map these rates onto MHC protein structures, and find that the fastest evolving codons are extremely enriched for autoimmune and infectious disease associations.

Strengths:

The study is comprehensive, relying on a large data set, state-of-the-art phylogenetic analyses and elegant tests of TSP. The results are not entirely novel, but a synthesis and re-analysis of previous findings is extremely valuable and timely.

Weaknesses:

Following the revision by the Authors I see mostly one weakness - Older literature on the subject is duly cited, but the discussion of the findings the context of this literature is limited.

Comments on revisions:

Lines 441-452 - In this section, you discuss an apparent paradox between long-lived balancing selection and strong directional selection, referencing elevated substitution rates. However, this issue is more nuanced and may not be best framed in terms of substitution rates. That terminology is common in phylogenetic analyses, where differences between sequences-or changes along phylogenetic branches-are often interpreted as true substitutions in the population genetic sense. In the case of MHC trees and the rates you're discussing here, the focus is more accurately on the rate at which new mutations become established within particular allelic lineages. So while this still concerns evolutionary rates at specific codons, equating them directly with substitution rates may be misleading. A more precise term or framing might be warranted in this context.

Reviewer #1 (Public review):

Summary and Strengths:

This work shows that the gene encoding Layilin is expressed preferentially in human skin Tregs, and that the fraction of Tregs expressing Layilin may overexpress genes related to T cell activation and adhesion. Expression of Layilin on Tregs would have no impact on activation markers or in vitro suppressive function. However, activation of Layilin either with a cross-linking antibody or collagen IV, its natural ligand, would promote cell adhesion via LFA1 activation. The in vivo functional role of Layilin in Tregs is studied in a conditional KO mouse model in a model of skin inflammation. Deletion of Layilin in Tregs led to an attenuation of the disease score and a reduction in the cutaneous lymphocyte infiltrate. This work is clearly innovative, but a number of major points limit its interest.

Weakness and major points:

(1) The number of panels and figures suggests that this story is quite complete but several data presented in the main figures do not provide essential information for a proper understanding of Layilin's role in Tregs.

Figures 1I, 1J, and the whole of Figure 2 could be placed as supplementary figures. Also, for Figure 3E, it would be preferable to show the percentage of cells expressing cytokines rather than their absolute numbers. In fact, the drop in the numbers of cytokine-producing cells is probably due solely to the drop in total cell numbers and not to a decrease in the proportion of cells expressing cytokines. If this is the case, these data should be shown in supplementary figures. Finally, Figures 4 and 5 could be merged.

(2) Some important data are not shown or not mentioned.

(a) It would be important to show the proportion of Treg, Tconv, and CD8 expressing Layilin in healthy skin and in patients developing psoriasis, as well as in the blood of healthy subjects.<br /> (b) We lack information to be convinced that there is enrichment for migration and adhesion genes in Layilin+ Tregs in the GSEA data. The authors should indicate what geneset libraries they used. Indeed, it is tempting to show only the genesets that give results in line with the message you want to get across. If these genesets come from public banks, the bank used should be indicated, and the results of all gene sets shown in an unbiased way. In addition, it should be indicated whether the analyses were performed on untransformed or pseudobulk scRNAseq data analyses. Finally, it would be preferable to confirm the GSEA data with z-score analyses, as Ingenuity does, for example. Indeed, in GSEA-type analyses, there are genes that have activating but also inhibiting effects on a pathway in a given gene set.<br /> (c) For all FACS data, the raw data should be shown as histograms or dot plots for representative samples.<br /> (d) For Figure 5B, the number of samples analyzed is insufficient to draw clear conclusions.

(3) For Figs. 4 and 5, the design of the experiment poses a problem. Indeed, the comparison between Layn+ and Layn- cells may, in part, not be directly linked to the expression or absence of expression of this protein. Indeed, Layn+ and Layn- Tregs may constitute populations with different biological properties, beyond the expression of Layn. However, in the experiment design used here, a significant fraction of the sorted Layn- Tregs will be cells belonging to the population that has never expressed this protein. It would have been preferable to sort first the Layn+ Tregs, then knock down this protein and re-sort the Layn- Tregs and Layn+ Tregs. If this experiment is too cumbersome to perform, I agree that the authors should not do it. However, it would be important to mention the point I have just made in the text.

Reviewer #2 (Public review):

Summary:

In their manuscript, Gouirand et al. report on the role of Layilin expression for the motility and suppressive capacity of regulatory T cells (Tregs). In previous studies, the authors had already demonstrated that Layilin is expressed on Tregs, that it acts as a negative regulator of their suppressive capacity, that it functions to anchor Tregs in non-lymphoid tissues, and that it enhances the adhesive properties of Layilin-expressing cells by co-localization with the integrin αLβ2 (LFA-1). Building on these published data, the authors now show that Layilin is highly expressed on a subset of clonally expanded effector Tregs in both healthy and psoriatic skin and that deletion of Layilin in Tregs in vivo resulted in significantly attenuated skin inflammation. Furthermore, the authors addressed the molecular mechanism by which Layilin affects the suppressive capacity of Tregs and showed that Layilin increased Treg adhesion via modulation of LFA-1, resulting in distinct cytoskeletal changes.

Strengths:

Certainly, the strength of this study lies in the combination of data from mouse and human models.

Weaknesses:

Some of the conclusions drawn by the authors must be treated with caution, as the experimental conditions were not always appropriate, leading to a risk of misinterpretation.

Reviewer #3 (Public review):

Summary:

Gouirand et al explore the function of Layilin on Treg in the context of psoriasis using both patient samples and a conditional mutant mouse model. They perform functional analysis in the patient samples using Cas9-mediated deletion. The authors suggest that Layilin works in concert with integrins to bind collagen IV to attenuate cell movement.

The work is well done and built on solid human data. The report is a modest advance from the authors' previous report in 2021 that focused on tumor responses, with this report focusing on psoriasis. There are some experimental concerns that should be considered.

Strengths:

(1) Good complementation of patient and animal model data.

(2) Solid experimentation using state-of-the-art approaches.

(3) There is clearly a biological effect of LAYN deficiency in the mouse model.

(4) The report adds some new information to what was already known from the previous reports.

Weaknesses:

(1) It is not clear that the assays used for functional analysis of the patient samples were optimal.

(2) Several conclusions are not fully substantiated.

(3) The report is lacking some experimental details.

Reviewer #1 (Public review):

Summary:

Aicardi-Goutières Syndrome (AGS) is a genetic disorder that primarily affects the brain and immune system through excessive interferon production. The authors sought to investigate the role of microglia in AGS by first developing bone-marrow-derived progenitors in vitro that carry the estrogen-regulated (ER) Hoxb8 cassette, allowing them to expand indefinitely in the presence of estrogen and differentiate into macrophages when estrogen is removed. When injected into the brains of Csf1r-/- mice, which lack microglia, these cells engraft and resemble wild-type (WT) microglia in transcriptional and morphological characteristics, although they lack Sall1 expression. The authors then generated CRISPR-Cas9 Adar1 knockout (KO) ER-Hoxb8 macrophages, which exhibited increased production of inflammatory cytokines and upregulation of interferon-related genes. This phenotype could be rescued using a Jak-Stat inhibitor or by concurrently mutating Ifih1 (Mda5). However, these Adar1-KO macrophages fail to successfully engraft in the brain of both Csf1r-/- and Cx3cr1-creERT2:Csf1rfl/fl mice. To overcome this, the authors used a mouse model with a patient-specific Adar1 mutation (Adar1 D1113H) to derive ER-Hoxb8 bone marrow progenitors and macrophages. They discovered that Adar1 D1113H ER-Hoxb8 macrophages successfully engraft the brain, although at lower levels than WT-derived ER-Hoxb8 macrophages, leading to increased production of Isg15 by neighboring cells. These findings shed new light on the role of microglia in AGS pathology.

Strengths:

The authors convincingly demonstrate that ER-Hoxb8 differentiated macrophages are transcriptionally and morphologically similar to bone marrow-derived macrophages. They also show evidence that when engrafted in vivo, ER-Hoxb8 microglia are transcriptomically similar to WT microglia. Furthermore, ER-Hoxb8 macrophages engraft the Csf1r-/- brain with high efficiency and rapidly (2 weeks), showing a homogenous distribution. The authors also effectively use CRISPR-Cas9 to knock out TLR4 in these cells with little to no effect on their engraftment in vivo, confirming their potential as a model for genetic manipulation and in vivo microglia replacement.

Overall, this paper demonstrates an innovative approach to manipulating microglia using ER-Hoxb8 cells as surrogates. The authors present convincing evidence of the model's efficacy and potential for broader application in microglial research, given its ease of production and rapid brain engraftment potential in microglia-deficient mice. Using mouse-derived cells for transplantation reduces complications that can come with the use of human cell lines, highlighting the utility of this system for research in mouse models.

Reviewer #2 (Public review):

Summary:

Microglia have been implicated in brain development, homeostasis, and diseases. "Microglia replacement" has gain tractions in recent years, using primary microglia, bone marrow or blood-derived myeloid cells, or human iPSC-induced microglia. Here, the authors extended their previous work in the area and provide evidence to support: (1) Estrogen-regulated (ER) homeobox B8 (Hoxb8) conditionally immortalized macrophages from bone marrow can serve as stable, genetically manipulated cell lines. These cells are highly comparable to primary bone marrow-derived (BMD) macrophages in vitro, and, when transplanted into a microglia-free brain, engraft the parenchyma and differentiate into microglia-like cells (MLCs). Taking advantage of this model system, the authors created stable, Adar1-mutated ER-Hoxb8 lines using CRISPR-Cas9 to study the intrinsic contribution of macrophages to Aicardi-Goutières Syndrome (AGS) disease mechanism.

Strengths:

The studies are carefully designed and well-conducted. The imaging data and gene expression analysis are carried out at a high level of technical competences and the studies provide strong evidence that ER-Hoxb8 immortalized macrophages from bone marrow are a reasonable source for "microglia replacement" exercise. The findings are clearly presented, and the main message will be of general interest to the neuroscience and microglia communities.

Reviewer #2 (Public review):

In the manuscript, Ruhling et al propose a rapid uptake pathway that is dependent on lysosomal exocytosis, lysosomal Ca2+ and acid sphingomyelinase, and further suggest that the intracellular trafficking and fate of the pathogen is dictated by the mode of entry. Overall, this is manuscript argues for an important mechanism of a 'rapid' cellular entry pathway of S.aureus that is dependent on lysosomal exocytosis and acid sphingomyelinase and links the intracellular fate of bacterium including phagosomal dynamics, cytosolic replication and host cell death to different modes of uptake.

Key strength is the nature of the idea proposed, while continued reliance on inhibitor treatment combined with lack of phenotype for genetic knock out is a major weakness. While the authors argue a role for undetectable nano-scale Cer platforms on the cell surface caused by ASM activity, results do not rule out a SM independent role in the cellular uptake phenotype of ASM inhibitors.

The authors have attempted to address many of the points raised in the previous revision. While the new data presented provide partial evidence, the reliance on chemical inhibitors and lack of clear results directly documenting release of lysosomal Ca2+, or single bacterial tracking, or clear distinction between ASM dependent and independent processes dampen the enthusiasm.

I acknowledge the author's argument of different ASM inhibitors showing similar phenotypes across different assays as pointing to a role for ASM, but the lack of phenotype in ASM KO cells is concerning. The author's argument that altered lipid composition in ASM KO cells could be overcoming the ASM-mediated infection effects by other ASM-independent mechanisms is speculative, as they acknowledge, and moderates the importance of ASM-dependent pathway. The SM accumulation in ASM KO cells does not distinguish between localized alterations within the cells. If this pathway can be compensated, how central is it likely to be ?

The authors allude to lower phagosomal escape rate in ASM KO cells compared to inhibitor treatment, which appears to contradict the notion of uptake and intracellular trafficking phenotype being tightly linked. As they point out, these results might be hard to interpret. Could an inducible KD system recapitulate (some of) the phenotype of inhibitor treatment ? If S. aureus does not escape phagosome in macrophages, could it provide a system to potentially decouple the uptake and intracellular trafficking effects by ASM (or its inhibitor treatment) ?

The role of ASM on cell surface remains unclear. The hypothesis proposed by the authors that the localized generation of Cer on the surface by released ASM leads to generation of Cer-enriched platforms could be plausible, but is not backed by data, technical challenges to visualize these platforms notwithstanding. These results do not rule out possible SM independent effects of ASM on the cell surface, if indeed the role of ASM is confirmed by controlled genetic depletion studies.

The reviewer acknowledges technical challenges in directly visualizing lysosomal Ca2+ using the methods outlined. Genetically encoded lysosomal Ca2+ sensor such as Gcamp3-ML1 might provide better ways to directly visualize this during inhibitor treatment, or S. aureus infection.

Reviewer #1 (Public review):

Summary:

This important study by Chen et. al. describes a novel approach for optogentically evoking seizures in an etiologically relevant mouse model of epilepsy. The authors developed a model that can trigger seizures "on demand" using optogenetic stimulation of CA1 principal cells in mice rendered epileptic by an intra-hippocampal kainate (IHK) injection into CA3. The authors discuss their model in the context of the limitations of current animal models used in epilepsy drug development. In particular, their model addresses concerns regarding existing models where testing typically involves inducing acute seizures in healthy animals or waiting on infrequent, spontaneous seizures in epileptic animals.

Strengths:

A strength of this manuscript is that this approach may facilitate the evaluation of novel therapeutics since these evoked seizures, despite having some features that were significantly different from spontaneous seizures, are suggested to be sufficiently similar to spontaneous seizures which are more laborious to analyze. The data demonstrating the commonality of pharmacology and EEG features between evoked seizures and spontaneous seizures in epileptic mice, while also being different from evoked seizures in naïve mice, are convincing. The structural, functional, and behavioral differences between a seizure-naïve and epileptic mouse, which emerge due to the enduring changes occurring during epileptogenesis, are complex and important. Accordingly, this study highlights the importance of using mice that have underwent epileptogenesis as model organisms for testing novel therapeutics. Furthermore, this study positively impacts the wider epilepsy research community by investigating seizure semiology in these populations.

Weaknesses:

This study convincingly demonstrates that the feature space measurements for stimulus-evoked seizures in epileptic mice were significantly different from those in naïve mice; this result allows the authors to conclude that "seizures induced in chronically epileptic animals differed from those in naïve animals". However, the authors also conclude that "induced seizures resembled naturally occurring spontaneous seizures in epileptic animals" despite their own data demonstrating similar, albeit fewer, significant differences in feature space measurements. It is unclear if and what the threshold is whereby significant differences in these feature space measurements lead to the conclusion that the differences are meaningful, as in the comparison of epileptic and naïve mice, or not meaningful, as in the comparison of evoked and spontaneous seizures.

Reviewer #2 (Public review):

The authors aimed to develop an animal model of temporal lobe epilepsy (TLE) that will generate "on-demand" seizures and an improved platform to advance our ability to find new anti-seizure drugs (ASDs) for drug-resistant epilepsy (DRE). Unlike some of the work in this field, the authors are studying actual seizures, and hopefully events that are similar to actual epileptic seizures. To develop an optimized screening tool, however, one also needs high-throughput systems with actual seizures as a quantitative, rigorous, and reproducible outcome measures. The authors aim to provide such a model; however, this approach may be over-stated here and seems unlikely to address the critical issue of drug resistance, which is their most important claim.

Strengths:

- The authors have generated an animal model of "on demand" seizures, which could be used to screen new ASDs and potentially other therapies. The authors and their model make a good-faith effort to emulate the epileptic condition and to use seizure susceptibility or probability as a quantitative output measure.

- The events considered to be seizures appear to be actual seizures, with some evidence that the seizures are different from seizures in the naïve brain. Their effort to determine how different ASDs raise seizure probability or threshold to an optogenetic stimulus to the CA1 area of the rodent hippocampus is focused on an important problem, as many if not most ASD screening uses surrogate measures that may not be as well linked to actual epileptic seizures.

- Another concern is their stimulation of dorsal hippocampus, while ventral hippocampus would seem more appropriate.

- Use of optogenetic techniques allows specific stimulation of the targeted CA1 pyramidal cells, and it appears that this approach is reproducible and reliable with quantitative rigor.

- The authors have taken on a critically important problem, and have made a good-faith effort to address many of the technical concerns raised in the reviews, but the underlying problem of DRE remains.

Weaknesses:

- Although the model has potential advantages, it also has disadvantages. As stated by the authors, the pre-test work-load to prepare the model may not be worth the apparent advantages. And most important, the paper frequently mentions DRE but does not directly address it, and yet drug resistance is the critical issue in this field.

- Although the paper shows examples of actual seizures, there remains some concern that some of the events might not be seizures - or a homogeneous population of seizures. More quantitative assessment of the electrical properties (e.g., duration) of the seizures and their probability is likely to be more useful than the proposed quantification in the future of the behavioral seizure stages, because the former could be both more objective and automated, while the behavioral analysis of the seizures will likely be more subjective and less reliable (and also fraught with subjectivity and analytical problems). Nonetheless, the authors point that the presence of "Racine 3 or above" behavioral seizures (in addition to their electrical data) is a good argument that many (if not all) of the "seizures" are actual epileptic seizures.

- Optogenetic stimulation of CA1 provides cell-specificity for the stimulation, but it is not clear that this method would actually be better than electrical stimulation of a kindled rodent with superimposed hippocampal injury. The reader is unfortunately left with the concern of whether this model would be easier and more efficacious than kindling.

- Although the authors have taken on a critically important problem, and have combined a variety of technologies, this approach may facilitate more rapid screening of ASDs against actual seizures (beneficial), but it does not really address the fundamentally critical yet difficult problem of DRE. A critical issue for DRE that is not well-addressed relates to adverse effects, which is often why many ASDs are not well tolerated by many patients (e.g., LEV). Thus, we are left with: how does this address anti-seizure DRE?

- The focus of this paper seems to be more on seizures more than on epilepsy. In the absence of seizure spontaneity, the work seems to primarily address the issues of seizure spread and duration. Although this is useful, it does not seem to be addressing the question of what trips the system to generate a seizure.

An appraisal of whether the authors achieved their aims, and whether the results support their conclusions:

- The authors seem to have developed a new and useful model; however, it is not clear how this will address that core problem of DRE, which was their stated aim.

- A discussion of the likely impact of the work on the field, and the utility of the methods and data to the community.

- As stated before in the original review, the potential impact would primarily be aimed at the ETSP or a drug-testing CRO; however, much more work will be required to convince the epilepsy community that this approach will actually identify new ASDs for DRE. The approach is potentially time-consuming with a steep and potentially difficult optimization curve, and thus may not be readily adaptable to the typical epilepsy-models neuroscience laboratory.

Any additional context you think would help readers interpret or understand the significance of the work:

- The problem of DRE is much more complicated than described by the authors here; however, the paper could end up being more useful than is currently apparent. Although this work could be seen as technically - and maybe conceptually - elegant and a technical tour de force, will it "deliver on the promise"? Is it better than kindling for DRE? In attempting to improve the discovery process, how will the model move us to another level? Will this model really be any better than others, such as kindling?

Reviewer #3 (Public review):

This revised paper develops and characterizes a new approach for screening drugs for epilepsy. The idea is to increase the ability to study seizures in animals with epilepsy because most animal models have rare seizures. Thus, the authors use the existing intrahippocampal kainic acid (IHKA) mouse model, which can have very unpredictable seizures with long periods of time between seizures. This approach is of clear utility to researchers who may need to observe many seizure events per mouse during screening of antiseizure medications. A key strength is also that more utility can be derived from each individual mouse. The authors modified the IHKA model to inject KA into CA3 instead of CA1 in order to preserve the CA1 pyramidal cells that they will later stimulate. To express the excitatory opsin channelrhodopsin (ChR2) in area CA1, they use a virus that expresses ChR2 in cells that express the Thy-1 promoter. The authors demonstrate that CA3 delivery of KA can induce a very similar chronic epilepsy phenotype to the injection of KA in CA1 and show that optical excitation of CA1 can reliably induce seizures. The authors evaluate the impact of repeated stimulation on the reliability of seizure induction and show that seizures can be reliably induced by CA1 stimulation, at least for the short term (up to 16 days). These are strengths of the study.

However, there are several limitations: the seizures are evoked, not spontaneous. It is not clear how induced seizures can be used to investigate if antiseizure medication can reduce spontaneous seizures. Although seizure inducibility and severity can be assessed, the lack of spontaneous seizures is a limitation. To their credit, the authors show that electrophysiological signatures of induced vs spontaneous seizures are similar in many ways, but the authors also show several differences. Notably, the induced seizures are robustly inhibited by the antiseizure medication levetiracetam and variably but significantly inhibited by diazepam, similar to many mouse models with chronic recurrent seizure activity. One also wonders if using a mouse model with numerous seizures (such as the pilocarpine model) might be more efficient than using a modified IHKA protocol.

In this revised manuscript, the authors address some previous concerns related to definitions of seizures and events that are trains of spikes, sex as a biological variable, and present new images of ChR2 expression (but these images could be improved to see the cells more clearly). A few key concerns remain unaddressed, however. For example, it is still not clear that evoked seizures triggered by stimulating CA1 are similar to spontaneous seizures, regardless of the idea that CA1 plays a role in seizure disorders. It also remains unclear whether repeated activation of the hippocampal circuit will result in additional alterations to this circuit that affect the seizure phenotype over prolonged intervals (after 16 days). Furthermore, the use of SVM with the number of seizures being used as replicates (instead of number of mice) is inappropriate. Another theoretical concern is whether the authors are correct in suggesting that one will be able to re-use the mice for screening multiple drugs in a row.

Strengths:<br /> - The authors show that the IHKA model of chronic epilepsy can be modified to preserve CA1 pyramidal cells, allowing optogenetic stimulation of CA1 to trigger seizures.<br /> - The authors show that repeated optogenetic stimulation of CA1 in untreated mice can promote kindling and induce seizures, indeed generating two mouse models in total.<br /> - Many electrophysiological signatures are similar between the induced and spontaneous seizures, and induced seizures reliably respond to treatment with antiseizure medications.<br /> - Given that more seizures can be observed per mouse using on-demand optogenetics, this model enhances the utility of each individual mouse.<br /> - Mice of each sex were used.

Weaknesses:<br /> - Evaluation of seizure similarity using the SVM modeling and clustering is not sufficiently justified when using number of seizures as the statistical replicate (vs mice).<br /> - Related to the first concern, the utility of increasing number of seizures for enhancing statistical power is limited because standard practice is for sample size to be numbers of mice.<br /> - The term "seizure burden" usually refers to the number of spontaneous seizures per day, not the severity of the seizures themselves. Because the authors are evoking the seizures being studied, this study design precludes assessment of seizure burden.<br /> - It seems likely that repeatedly inducing seizures will have a long-term effect, especially in light of the downward slope at day 13-16 for induced seizures seen in Figure 4C. A duration of evaluation that is longer than 16 days is warranted.<br /> - Human epilepsy is extensively heterogeneous in both etiology and individual phenotype, and it may be hard to generalize the approach.

Reviewer #1 (Public review):

Strengths:

This work adds another mouse model for LAMA2-MD that re-iterates the phenotype of previously published models. Such as dy3K/dy3K; dy/dy and dyW/dyW mice. The phenotype is fully consistent with the data from others.

One of the major weaknesses of the manuscript initially submitted was the overinterpretation and the overstatements. The revised version is clearly improved as the authors toned-down their interpretation and now also cite the relevant literature of previous work.

Comments on revisions:

This is the second revision of a paper focusing on the generation of a CRISPR/Cas9-engineered mouse model for LAMA2-MD. I have reviewed the initial submission, the first revision, and now this second revision. While there have been improvements, several issues still need to be addressed by the authors. I will outline these points without dividing them into major and minor categories:

Introduction:

The statement regarding existing mouse models requires correction: The claim, "They were established in the pre-gene therapy era, leaving trace of engineering, such as bacterial elements in the Lama2 gene locus, thus unsuitable for testing various gene therapy strategies," is inaccurate. Current mouse models can indeed be used for testing gene therapy strategies, regardless of whether they contain elements in the Lama2 locus. The primary consideration is whether or not they express laminin-alpha2. Please revise this statement.<br /> Results Section:

scRNA-seq:

The authors note that they analyzed "a total of 8,111 cells from the dyH/dyH mouse brain and 8,127 cells from the WT mouse brain were captured using the 10X Genomics platform (Figure supplement 4A, B)." This is too few cells to support firm conclusions. Furthermore, there is a discrepancy in the referred figure S4, which indicates that 10,094 cells were analyzed for dyH/dyH mice and 10,496 for wild-type mice. Please correct this inconsistency.

Figure 5C displays differences in cell populations between wild-type and dyH/dyH mice. Given the low number of cells analyzed and the lack of replicates, these differences cannot be considered reliable. More samples should be analyzed to support these findings.

The data suggest a defect in the BBB for dyH/dyH mice, but this conclusion is based on minimal cell counts and remains purely correlative. If BBB issues exist, experimental validation is necessary, such as injecting dyes into the bloodstream to detect any leakage. I have previously highlighted this in my comments on earlier manuscript versions.

Bulk RNA-seq:

The number of samples analyzed here is substantial, making the data potentially more robust. These data could serve as a valuable resource for other researchers. However, it is important to note that all data are correlative and do not provide functional insights.

Overall:

The manuscript still lacks significant insights, partly because existing mouse models for LAMA2-MD have been extensively analyzed. While the bulk RNA-seq data offer some value as a resource, I recommend that the authors re-assess their writing and further temper their interpretations of the findings.

Reviewer #1 (Public review):

Translating discoveries from model organisms to humans is often challenging, especially in neuropsychiatric diseases, due to the vast gaps in the circuit complexities and cognitive capabilities. Kajtor et al. propose to bridge this gap in the fly models of Parkinson's disease (PD) by developing a new behavioural assay where flies respond to a moving shadow by modifying their locomotor activities. The authors believe the flies' response to the shadow approximates their escape response to an approaching predator. To validate this argument, they tested several PD-relevant transgenic fly lines and showed that some of them indeed have altered responses in their assay.

Strengths:

This single-fly-based assay is easy and inexpensive to set up, scalable and provides sensitive, quantitative estimates to probe flies' optomotor acuity. The behavioural data is detailed, and the analysis parameters are well-explained.

Weaknesses:

The authors have yet to link cellular physiology to behaviour. It will be interesting to see how future use of this assay helps uncover connections between cellular pathology and behavioural changes.

Reviewer #2 (Public review):

The manifestation and progression of neurodegenerative disorders is poorly understood. Many of the neuronal disorders start by presenting subtle changes in neuronal circuit and quantification and measurement of these subtle behavior responses could help one delineate the mechanisms involved. The present study very nicely uses the flies' behavioral response to predator-mimicking passing shadows to measure subtle changes in their behavior. The data from various fly genetic models of Parkinson's disease supports their claim. This single trial method is useful to capture the individual animal's response to the threatening stimuli but stops short of capturing the fine ambulatory responses which could provide further information on an individual's behavioral response. By capturing the fine features, the authors could get detailed observations, such as posture, gait or wing positioning for a better understanding the behavioral response to the passing shadow.

Reviewer #1 (Public review):

Summary:

This manuscript explores the multiple cell types present in the wall of murine collecting lymphatic vessels with the goal of identifying cells that initiate the autonomous action potentials and contractions needed to drive lymphatic pumping. Through the use of genetic models to delete individual genes or detect cytosolic calcium in specific cell types, the authors convincingly determine that lymphatic muscle cells are the origin of the action potential that triggers lymphatic contraction.

Strengths:

The experiments are rigorously performed, the data justify the conclusions and the limitations of the study are appropriately discussed.

There is a need to identify therapeutic targets to improve lymphatic contraction and this work helps identify lymphatic muscle cells as potential cellular targets for intervention.

Comments on revisions: The authors have addressed all of the reviewer comments. They should be congratulated on their precise and comprehensive study.

Reviewer #2 (Public review):

Summary:

This is a well written manuscript describing studies directed at identifying the cell type responsible for pacemaking in murine collecting lymphatics. Using state of the art approaches, the authors identified a number of different cell types in the wall of these lymphatics and then using targeted expression of Channel Rhodopsin and GCaMP, the authors convincingly demonstrate that only activation of lymphatic muscle cells produces coordinated lymphatic contraction and that only lymphatic muscle cells display pressure-dependent Ca2+ transients as would be expected of a pacemaker in these lymphatics.

Strengths:

The use of targeted expression of channel rhodopsin and GCaMP to test the hypothesis that lymphatic muscle cells serve as the pacemakers in musing lymphatic collecting vessels.

Weaknesses:

The only significant weakness was the lack of quantitative analysis of most of the imaging data shown in Figures 1-11. In particular the colonization analysis should be extended to show cells not expected to demonstrate colocalization as a negative control for the colocalization analysis that the authors present. These weaknesses have been resolved by revision and addition of new and novel RNAseq data, additional colocalization data and membrane potential measurements.

Comments on revisions: No additional concerns.

Reviewer #3 (Public review):

Summary:

Zawieja et al. aimed to identify the pacemaker cells in the lymphatic collecting vessels. Authors have used various Cre-based expression systems and optogentic tools to identify these cells. Their findings suggest these cells are lymphatic muscle cells that drive the pacemaker activity in the lymphatic collecting vessels.

Strengths:

The authors have used multiple approaches to test their hypothesis. Some findings are presented as qualitative images, while some quantitative measurements are provided.

Weaknesses:<br /> - More quantitative measurements.<br /> - Possible mechanisms associated with the pacemaker activity.<br /> - Membrane potential measurements.

Comments on revisions: I do not have any additional comments.

Reviewer #2 (Public review):

Summary:

Sukhina et al. uses a chronic murine dietary restriction model to investigate the cellular mechanisms underlying nutritionally acquired immunodeficiency as well as the consequences of a refeeding intervention. The authors report a substantial impact of undernutrition to the myeloid compartment, which is not rescued by refeeding despite rescue of other phenotypes including lymphocyte levels, and which is associated with maintained partial susceptibility to bacterial infection.

Strengths:

Overall, this is a nicely executed study with an appropriate number of mice, robust phenotypes, and interesting conclusions, and the text is very well written. The authors' conclusions are generally well-supported by their data.

Weaknesses:

There is little evaluation of known critical drivers of myelopoiesis (e.g. PMID 20535209, 26072330, 29218601) over the course of the 40% diet, which would be of interest with regard to comparing this chronic model to other more short-term models of undernutrition.

Further, the microbiota, well-established to be regulated by undernutrition (e.g. PMID 22674549, 27339978, etc.), and also well-established to be a critical regulator of hematopoiesis/myelopoiesis (e.g. PMID 27879260, 27799160, etc.), should be studied in any future explorations using this model.

The authors have recognized these limitations to the study in their discussion.

Reviewer #3 (Public review):

This communication from Sukhina et al argues that a period of malnutrition (modeled by caloric restriction) causes lasting immune deficiencies (myelopoesis) not rescued by re-feeding. This is a potentially important paper exploring the effects of malnutrition on immunity, which is a clinically important topic. The revised study adds some details with respect to kinetics of immune compartment and body weight changes, but most aspects raised by the referees were deferred experimentally. Several textual changes have been made to avoid over-interpreting their data. My overall assessment of this revised study is similar to my impression before, which is that while the observations are interesting, there is both a lack of mechanistic understanding of the phenomena and a lack of resolution/detail about the phenomena itself.

Reviewer #1 (Public review):

Summary:

Desingu et al. show that JEV infection reduces SIRT2 expression. Upon JEV infection, 10-day-old SIRT2 KO mice showed increased viral titer, more severe clinical outcomes, and reduced survival. Conversely, SIRT2 overexpression reduced viral titer, clinical outcomes, and improved survival. Transcriptional profiling shows dysregulation of NF-KB and expression of inflammatory cytokines. Pharmacological NF-KB inhibition reduced viral titer. The authors conclude that SIRT2 is a regulator of JEV infection.

Strengths:

This paper is novel because sirtuins have been primarily studied for aging, metabolism, stem cells/regeneration. Their role in infection has not been explored until recently. Indeed, Barthez et al. showed that SIRT2 protects aged mice from SARS-CoV-2 infection (Barthez, Cell Reports 2025). Therefore, this is a timely and novel research topic. Mechanistically, the authors showed that SIRT2 suppresses the NF-KB pathway. Interestingly, SIRT2 has also been shown recently to suppress other major inflammatory pathways, such as cGAS-STING (Barthez, Cell Reports 2025) and the NLRP3 inflammasome (He, Cell Metabolism 2020; Luo, Cell Reports 2019). Together, these findings support the emerging concept that SIRT2 is a master regulator of inflammation.

Weaknesses:

(1) Figures 2 and 3. Although SIRT2 KO mice showed increased viral titer, more severe clinical outcomes, and reduced survival upon JEV infection, the difference is modest because even WT mice exhibited very severe disease at this viral dose. The authors should perform the experiment using a sub-lethal viral dose for WT mice, to allow the assessment of increased clinical outcomes and reduced survival in KO mice.

(2) Figure 5K-N, the authors examined the expression of inflammatory cytokines in WT and SIRT2 KO cells upon JEV infection, in line with the dysregulation of NF-kB. It has been shown recently that SIRT2 also regulates the cGAS-STING pathway (Barthez, Cell Reports 2025) and the NLRP3 inflammasome (He, Cell Metabolism 2020; Luo, Cell Reports 2019). Do you also observe increased IFNb, IL1b, and IL18 in SIRT2 KO cells upon JEV infection? This may indicate that SIRT2 regulates systemic inflammatory responses and represents a potent protection upon viral infection. This is particularly important because in Figure 7F, the authors showed that SIRT2 overexpression reduced viral load even when NF-KB is inhibited, suggesting that NF-KB is not the only mediator of SIRT2 to suppress viral infection.

Reviewer #2 (Public review):

The manuscript by Desingu et al., explores the role of SIRT2 in regulating Japanese Encephalitis Virus (JEV) replication and disease progression in rodent models. Using both an in vitro and an in vivo approach, the authors demonstrate that JEV infection leads to decreased SIRT2 expression, which they hypothesize is exploited by JEV for viral replication. To test this hypothesis, the authors utilize SIRT2 inhibition (via AGK2 or genetic knockout) and demonstrate that it leads to increased viral load and worsens clinical outcomes in JEV-infected mice. Conversely, SIRT2 overexpression via an AAV delivery system reduces viral replication and improves survival among infected mice. The study proposes a mechanism in which SIRT2 suppresses JEV-induced autophagy and inflammation by deacetylating NF-κB, thereby reducing Beclin-1 expression (an NF-κB-dependent gene) and autophagy, which the authors consider a pathway that JEV exploits for replication. Transcriptomic analysis further supports that SIRT2 deficiency leads to NF-κB-driven cytokine hyperactivation. Additionally, pharmacological inhibition of NF-κB using Bay 11 (an IKK inhibitor) results in reduced viral load and improved clinical pathology in WT and SIRT2 KO mice. Overall, the findings from Desingu et al. are generally supported by the data and suggest that targeting SIRT2 may serve as a promising therapeutic approach for JEV infection and potentially other RNA viruses that SIRT2 helps control. However, the paper does fall short in some areas. Please see below for our comments to help improve the paper.

Reviewer #1 (Public review):

Fombellida-Lopez and colleagues describe the results of an ART intensification trial in people with HIV infection (PWH) on suppressive ART to determine the effect of increasing the dose of one ART drug, dolutegravir, on viral reservoirs, immune activation, exhaustion, and circulating inflammatory markers. The authors hypothesize that ART intensification will provide clues about the degree to which low-level viral replication is occurring in circulation and in tissues despite ongoing ART, which could be identified if reservoirs decrease and/or if immune biomarkers change. The trial design is straightforward and well-described, and the intervention appears to have been well tolerated. The investigators observed an increase in dolutegravir concentrations in circulation, and to a lesser degree in tissues, in the intervention group, indicating that the intervention has functioned as expected (ART has been intensified in vivo). Several outcome measures changed during the trial period in the intervention group, leading the investigators to conclude that their results provide strong evidence of ongoing replication on standard ART. The results of this small trial are intriguing, and a few observations in particular are hypothesis-generating and potentially justify further clinical trials to explore them in depth. However, I am concerned about over-interpretation of results that do not fully justify the authors' conclusions.

(1) Trial objectives: What was the primary objective of the trial? This is not clearly stated. The authors describe changes in some reservoir parameters and no changes in others. Which of these was the primary outcome? No a priori hypothesis / primary objective is stated, nor is there explicit justification (power calculations, prior in vivo evidence) for the small n, unblinded design, and lack of placebo control. In the abstract (line 36, "significant decreases in total HIV DNA") and conclusion (lines 244-246), the authors state that total proviral DNA decreased as a result of ART intensification. However, in Figures 2A and 2E (and in line 251), the authors indicate that total proviral DNA did not change. These statements are confusing and appear to be contradictory. Regarding the decrease in total proviral DNA, I believe the authors may mean that they observed transient decrease in total proviral DNA during the intensification period (day 28 in particular, Figure 2A), however this level increases at Day 56 and then returns to baseline at Day 84, which is the source of the negative observation. Stating that total proviral DNA decreased as a result of the intervention when it ultimately did not is misleading, unless the investigators intended the day 28 timepoint as a primary endpoint for reservoir reduction - if so, this is never stated, and it is unclear why the intervention would then be continued until day 84? If, instead, reservoir reduction at the end of the intervention was the primary endpoint (again, unstated by the authors), then it is not appropriate to state that the total proviral reservoir decreased significantly when it did not.

(2) Intervention safety and tolerability: The results section lacks a specific heading for participant safety and tolerability of the intervention. I was wondering about clinically detectable viremia in the study. Were there any viral blips? Was the increased DTG well tolerated? This drug is known to cause myositis, headache, CPK elevation, hepatotoxicity, and headache. Were any of these observed? What is the authors' interpretation of the CD4:8 ratio change (line 198)? Is this a significant safety concern for a longer duration of intensification? Was there also a change in CD4% or only in absolute counts? Was there relative CD4 depletion observed in the rectal biopsy samples between days 0 and 84? Interestingly, T cells dropped at the same timepoints that reservoirs declined... how do the authors rule out that reservoir decline reflects transient T cell decline that is non-specific (not due to additional blockade of replication)?

(3) The investigators describe a decrease in intact proviral DNA after 84 days of ART intensification in circulating cells (Figure 2D), but no changes to total proviral DNA in blood or tissue (Figures 2A and 2E; IPDA does not appear to have been done on tissue samples). It is not clear why ART intensification would result in a selective decrease in intact proviruses and not in total proviruses if the source of these reservoir cells is due to ongoing replication. These reservoir results have multiple interpretations, including (but not limited to) the investigators' contention that this provides strong evidence of ongoing replication. However, ongoing replication results in the production of both intact and mutated/defective proviruses that both contribute to reservoir size (with defective proviruses vastly outnumbering intact proviruses). The small sample size and well-described heterogeneity of the HIV reservoir (with regard to overall size and composition) raise the possibility that the study was underpowered to detect differences over the 84-day intervention period. No power calculations or prior studies were described to justify the trial size or the duration of the intervention. Readers would benefit from a more nuanced discussion of reservoir changes observed here.

(4) While a few statistically significant changes occurred in immune activation markers, it is not clear that these are biologically significant. Lines 175-186 and Figure 3: The change in CD4 cells + for TIGIT looks as though it declined by only 1-2%, and at day 84, the confidence interval appears to widen significantly at this timepoint, spanning an interquartile range of 4%. The only other immune activation/exhaustion marker change that reached statistical significance appears to be CD8 cells + for CD38 and HLA-DR, however, the decline appears to be a fraction of a percent, with the control group trending in the same direction. Despite marginal statistical significance, it is not clear there is any biological significance to these findings; Figure S6 supports the contention that there is no significant change in these parameters over time or between groups. With most markers showing no change and these two showing very small changes (and the latter moving in the same direction as the control group), these results do not justify the statement that intensifying DTG decreases immune activation and exhaustion (lines 38-40 in the abstract and elsewhere).

(5) There are several limitations of the study design that deserve consideration beyond those discussed at line 327. The study was open-label and not placebo-controlled, which may have led to some medication adherence changes that confound results (authors describe one observation that may be evidence of this; lines 146-148). Randomized/blinded / cross-over design would be more robust and help determine signal from noise, given relatively small changes observed in the intervention arm. There does not seem to be a measurement of key outcome variables after treatment intensification ceased - evidence of an effect on replication through ART intensification would be enhanced by observing changes once intensification was stopped. Why was intensification maintained for 84 days? More information about the study duration would be helpful. Table 1 indicates that participants were 95% male. Sex is known to be a biological variable, particularly with regard to HIV reservoir size and chronic immune activation in PWH. Worldwide, 50% of PWH are women. Research into improving management/understanding of disease should reflect this, and equal participation should be sought in trials. Table 1 shows differing baseline reservoir sizes betweenthe control and intervention groups. This may have important implications, particularly for outcomes where reservoir size is used as the denominator.

(6) Figure 1: the increase in DTG levels is interesting - it is not uniform across participants. Several participants had lower levels of DTG at the end of the intervention. Though unlikely to be statistically significant, it would be interesting to evaluate if there is a correlation between change in DTG concentrations and virologic / reservoir / inflammatory parameters. A positive relationship between increasing DTG concentration and decreased cell-associated RNA, for example, would help support the hypothesis that ongoing replication is occurring.

(7) Figure 2: IPDA in tissue- was this done? scRNA in blood (single copy assay) - would this be expected to correlate with usCaRNA? The most unambiguous result is the decrease in cell-associated RNA - accompanying results using single-copy assay in plasma would be helpful to bolster this result. The use of the US RNA / Total DNA ratio is not helpful/difficult to interpret since the control and intervention arms were unmatched for total DNA reservoir size at study entry.

Reviewer #2 (Public review):

Summary:

An intensification study with a double dose of 2nd generation integrase inhibitor with a background of nucleoside analog inhibitors of the HIV retrotranscriptase in 2, and inflammation is associated with the development of co-morbidities in 20 individuals randomized with controls, with an impact on the levels of viral reservoirs and inflammation markers. Viral reservoirs in HIV are the main impediment to an HIV cure, and inflammation is associated with co-morbidities.

Strengths:

The intervention that leads to a decrease of viral reservoirs and inflammation is quite straightforward forward as a doubling of the INSTI is used in some individuals with INSTI resistance, with good tolerability.

This is a very well documented study, both in blood and tissues, which is a great achievement due to the difficulty of body sampling in well-controlled individuals on antiretroviral therapy. The laboratory assays are performed by specialists in the field with state-of-the art quantification assays. Both the introduction and the discussion are remarkably well presented and documented.

The findings also have a potential impact on the management of chronic HIV infection.

Weaknesses:

I do not think that the size of the study can be considered a weakness, nor the fact that it is open-label either.