Reviewer #1 (Public review):

Summary:

Cook et al. have presented an important study on the transcriptomic and epigenomic signature underlying craniofacial development in marsupials. Given the lack of a dunnart genome, the authors also prepared long and short-read sequence datasets to assemble and annotate a novel genome to allow for the mapping of RNAseq and ChIPseq data against H3K4me3 and H3K27ac, which allowed for the identification of putative promoter and enhancer sites in dunnart. They found that genes proximal to these regulatory loci were enriched for functions related to bone, skin, muscle and embryonic development, highlighting the precocious state of newborn dunnart facial tissue. When compared with mouse, the authors found a much higher proportion of promoter regions aligned between species than for enhancer regions, and subsequent profiling identified regulatory elements conserved across species and are important for mammalian craniofacial development. In contrast, the identification of dunnart-specific enhancers and patterns of RNA expression further confirm the precocious state of muscle development, as well as for sensory system development, in dunnart suggesting that early formation of these features are critical for neonate marsupials likely to assist with detecting and responding to cues that direct the joeys to the mother's teat after birth. This is one of the few epigenomic studies performed in marsupials (of any organ) and the first performed in fat-tailed dunnart (also of any organ). Marsupials are emerging as an important model for studying mammalian development and evolution and the authors have performed a novel and thorough analysis, impressively including the assembly of a new marsupial reference genome that will benefit many future studies.

Strengths:

The study provides multiple pieces of evidence supporting the important role enhancer elements play in mammalian phenotypic evolution, namely the finding of a lower proportion of peaks present in both dunnart and mouse for enhancers than for promoters, and dunnart showing more genes uniquely associated with it's active enhancers than any other combination of mouse and dunnart samples, whereas this pattern was less pronounced than for promoter-associated genes. In addition, rigorous parameters were used for the cross-species analyses to identify the conserved regulatory elements and the dunnart-specific enhancers. For example, for the results presented in Figure 1, I agree that it is a little surprising that the average promoter-TSS distance is greater than that for enhancers, but that this could be related to the possible presence of unannotated transcripts between genes. The authors addressed this well by examining the distribution of promoter-TSS distances and using proximal promoters (cluster #1) as high confidence promoters for downstream analyses.

The genome assembly method was thorough, using two different long read methods (Pacbio and ONT) to generate the long reads for contig and scaffold construction, increasing the quality of the final assembled genome.

Weaknesses:

Biological replicates of facial tissue were collected at a single developmental time point of the fat-tailed dunnart within the first postnatal day (P0), and analysed this in the context of similar mouse facial samples from the ENCODE consortium at six developmental time points, where previous work from the authors have shown that the younger mouse samples (E11.5-12.5) approximately corresponds to the dunnart developmental stage (Cook et al. 2021). However, it would be useful to have samples from at least one older dunnart time point, for example, at a developmental stage equivalent to mouse E15.5. This would provide additional insight into the extent of accelerated face development in dunnart relative to mouse, i.e. how long do the regulatory elements that activated early in dunnart remain active for and does their function later influence other aspects of craniofacial development?

The authors refer to the development of the CNS being delayed in marsupials relative to placental mammals, however, evidence shows how development of the dunnart brain (whole brain or cortex) is protracted compared to mouse, by a factor of at least 2 times, rather than delayed per se (Workman et al. 2013; Paolino et al. 2023). In addition, there is evidence that cortical formation and cell birth may begin at approximately the same stage across species equivalent to the neonate period in dunnart (E10.5 in mouse), and that shortly after this at the stage equivalent to mouse E12.5, the dunnart cortex shows signs of advanced neurogenesis followed by a protracted phase of neuronal maturation (Paolino et al. 2023). Therefore, it is possible that marsupial CNS development appears delayed relative to mouse but instead begins at the same stage and then proceeds to develop on a different timing scale.

Reviewer #1 (Public review):

Summary:

This study by Howe and colleagues investigates the role of the posterolateral cortical amygdala (plCoA) in mediating innate responses to odors, specifically attraction and aversion. By combining optogenetic stimulation, single-cell RNA sequencing, and spatial analysis, the authors identify a topographically organized circuit within plCoA that governs these behaviors. They show that specific glutamatergic neurons in the anterior and posterior regions of plCoA are responsible for driving attraction and avoidance, respectively, and that these neurons project to distinct downstream regions, including the medial amygdala and nucleus accumbens, to control these responses.

Strengths:

The major strength of the study is the thoroughness of the experimental approach, which combines advanced techniques in neural manipulation and mapping with high-resolution molecular profiling. The identification of a topographically organized circuit in plCoA and the connection between molecularly defined populations and distinct behaviors is a notable contribution to understanding the neural basis of innate motivational responses. Additionally, the use of functional manipulations adds depth to the findings, offering valuable insights into the functionality of specific neuronal populations.

Weaknesses:

There are some weaknesses in the study's methods and interpretation. The lack of clarity regarding the behavior of the mice during head-fixed imaging experiments raises the possibility that restricted behavior could explain the absence of valence encoding at the population level. Furthermore, while the authors employ chemogenetic inhibition of specific pathways, the rationale for this choice over optogenetic inhibition is not fully addressed, and this could potentially affect the interpretation of the results. Additionally, the choice of the mplCoA for manipulation, rather than the more directly implicated anterior and posterior subregions, is not well-explained, which could undermine the conclusions drawn about the topographic organization of plCoA.

Despite these concerns, the work provides significant insights into the neural circuits underlying innate behaviors and opens new avenues for further research. The findings are particularly relevant for understanding the neural basis of motivational behaviors in response to sensory stimuli, and the methods used could be valuable for researchers studying similar circuits in other brain regions. If the authors address the methodological issues raised, this work could have a substantial impact on the field, contributing to both basic neuroscience and translational research on the neural control of behavior.

Reviewer #1 (Public review):

Summary:

MHC (Major Histocompatibility Complex) genes have long been mentioned as cases of trans-species polymorphism (TSP), where alleles might have their most recent common ancestor with alleles in a different species, rather than other alleles in the same species (e.g., a human MHC allele might coalesce with a chimp MHC allele, more recently than the two coalesce with other alleles in either species). This paper provides a more complete estimate of the extent and ages of TSP in primate MHC loci. The data clearly support deep TSP linking alleles in humans to (in some cases) old world monkeys, but the amount of TSP varies between loci.

Strengths:

The authors use publicly available datasets to build phylogenetic trees of MHC alleles and loci. From these trees they are able to estimate whether there is compelling support for Trans-species polymorphisms (TSPs) using Bayes Factor tests comparing different alternative hypotheses for tree shape. The phylogenetic methods are state-of-the-art and appropriate to the task.

The authors supplement their analyses of TSP with estimates of selection (e.g., dN/dS ratios) on motifs within the MHC protein. They confirm what one would suspect: classical MHC genes exhibit stronger selection at amino acid residues that are part of the peptide binding region, and non-classical MHC exhibit less evidence of selection. The selected sites are associated with various diseases in GWAS studies.

Weaknesses:

An implication drawn from this paper (and previous literature) is that MHC has atypically high rates of TSP. However, rates of TSP are not estimated for other genes or gene families, so readers have no basis of comparison. No framework to know whether the depth and frequency of TSP is unusual for MHC family genes, relative to other random genes in the genome, or immune genes in particular. I expect (from previous work on the topic), that MHC is indeed exceptional in this regard, but some direct comparison would provide greater confidence in this conclusion.

Given the companion paper's evidence of genic gain/loss, it seems like there is a real risk that the present study under-estimates TSP, if cases of TSP have been obscured by the loss of the TSP-carrying gene paralog from some lineages needed to detect the TSP. Are the present analyses simply calculating rates of TSP of observed alleles, or are you able to infer TSP rates conditional on rates of gene gain/loss?

Figure 5 (and 6) provide regression model fits (red lines in panel C) relating evolutionary rates (y axis not labeled) to site distance from the peptide binding groove, on the protein product. This is a nice result. I wonder, however, whether a linear model (as opposed to non-linear) is the most biologically reasonable choice, and whether non-linear functions have been evaluated. The authors might consider generalized additive models (GAMs) as an alternative that relaxes linearity assumptions.

The connection between rapidly evolving sites, and disease associations (lines 382-3) is very interesting. However, this is not being presented as a statistical test of association. The authors note that fast-evolving amino acids all have at least one association: but is this really more disease-association than a random amino acid in the MHC? Or, a randomly chosen polymorphic amino acid in MHC? A statistical test confirming an excess of disease associations would strengthen this claim.

Reviewer #1 (Public review):

Summary:

Meissner et al describe an update on the collection of split-GAL4 lines generated by a consortium led by Janelia Research Campus. This follows the same experimental pipeline described before and presents as a significant increment to the present collection. This will strengthen the usefulness and relevance of "splits" as a standard tool for labs that already use this tool and attract more labs and researchers to use it.

Strengths:

This manuscript presents a solid step to establish Split-GAL4 lines as a relevant tool in the powerful Drosophila toolkit. Not only the raw number of available lines contribute to the relevance of this tool in the "technical landscape" of genetic tools, but additional features of this effort contribute to the successful adoption. These include:

(1) A description of expression patterns in the adult and larvae, expanding the "audience" for these tools<br /> (2) A classification of line combination according to quality levels, which provides a relevant criterion while deciding to use a particular set of "splits".<br /> (3) Discrimination between male and female expression patterns, providing hints regarding the potential role of these gender-specific circuits.<br /> (4) The search engine seems to be user-friendly, facilitating the retrieval of useful information.<br /> (5) An acknowledgement of the caveats and challenges that splits (like any other genetic tool) can carry.<br /> Overall, the authors employed a pipeline that maximizes the potential of the Split-GAL4 collection to the scientific community.

Weaknesses:

My concerns were resolved regarding the existence of caveats while using these tools that researchers should be aware of, particularly those using them for the first time.

Reviewer #1 (Public review):

Summary:

This study provides convincing data showing that expression of the PIK3R1(deltaExon11) dominant negative mutation in Activated PI3K Delta Syndrome 1/2 (APDS1/2) patient-derived cells reduces AKT activation and p110δ protein levels. Using a 3T3-L1 model cell system, the authors show that overexpressed p85α(deltaExon 11) displays reduced association with the p110α catalytic subunit but strongly interacts with Irs1/2. Overexpression of PIK3R1 dominant negative mutants inhibit AKT phosphorylation and reduce cellular differentiation of preadipocytes. The experimental design, interpretation, and quantification broadly support the authors' conclusions, which establishes a new paradigm that warrants future studies.

Strengths:

The strength of this study is the clear results derived from Western blots analysis of cell signaling markers (e.g. pAKT1), and co-immunoprecipitation of PI3K holoenzyme complexes and associated regulatory factors (e.g. Irs1/2). The authors analyze a variety of PIK3R1 mutants (i.e. deltaExon11, E489K, R649W, and Y657X), which reveals a range of phenotypes that support the proposed model for dominant negative activity. The use of clonal cell lines with doxycycline induced expression of the PIK3R1 mutants (deltaExon 11, R649W, and Y657X) provides convincing experimental data concerning the relationship between p85α mutant expression and AKT phosphorylation in vivo. This approach for overexpression is excellent and should be utilized more broadly by cell biologists. The authors convincingly show that p85α(deltaExon11, R649W, or Y657X) is unable to associate with p110α but instead more strongly associates with Irs1/2 compared to wild type p85α. Overall, this article does a great job of motivating future studies of SHORT and APDS2 PIK3R1 mutants expressed from their endogenous loci (e.g. knock-in mice).

Weaknesses:

The limitations for this study lie in the complexity of the cell signaling pathway under investigation, rather than a lack of rigor by the authors. Future experimentation will help reconcile the cell type specific differences (e.g. APDS2 patient derived cells vs. the 3T3-L1 cell model system) in PIK3R1 mutant behavior reported by the authors. This is also intimately linked to variable expression of PIK3R1 mutants and cell-type specific regulatory factors. Although beyond the scope of this work, an unbiased proteomic study that broadly evaluates the cell signaling landscape could provide a more holistic understanding of the APDS2 and SHORT mutants compared to a candidate-based approach. Additional structural biochemistry of the p110α/p85α(deltaExon 11) complex is needed to explain why PIK3R1 mutant regulatory subunits do not strongly associate with the p110 catalytic subunit. A more comprehensive biochemical analysis of p110α/p85α, p110β/p85α, and p110δ/p85α mutant protein complexes will also be necessary to explain various cell signaling phenotypes. A minor limitation of this study is the use of single end point assays to measure PI3K lipid kinase activity in the presence of one regulatory input (i.e. RTK-derived pY peptide). An expanded biochemical analysis of purified mutant PI3K complexes across the canonical membrane signaling landscape will be important for deciphering how competition between wild-type and mutant regulatory subunits are regulated in different cell signaling contexts.

Reviewer #1 (Public review):

This study presents an investigation into the physiological functions of RIPK1 within the context of liver physiology, particularly during short-term fasting. Through the use of hepatocyte-specific Ripk1-deficient mice (Ripk1Δhep), the authors embarked on an examination of the consequences of Ripk1 deficiency in hepatocytes under fasting conditions. They discovered that the absence of RIPK1 sensitized the liver to acute injury and hepatocyte apoptosis during fasting, a finding of significant interest given the crucial role of the liver in metabolic adaptation. Employing a combination of transcriptomic profiling and single-cell RNA sequencing techniques, the authors uncovered intricate molecular mechanisms underlying the exacerbated proinflammatory response observed in Ripk1Δhep mice during fasting. While the investigation offers valuable insights into the consequences of Ripk1 deficiency in hepatocytes during fasting conditions, there appears to be a primarily descriptive nature to the study with a lack of clear connection between the experiments. Thus, a stronger focus is warranted, particularly on understanding the dialogue between hepatocytes and macrophages. Moreover, the data would benefit from reinforcement through additional experiments such as Western blotting, flow cytometry, and rescue experiments, which would offer a more quantitative aspect to the findings. By incorporating these enhancements, the study could achieve a more comprehensive understanding of the underlying mechanisms and ultimately strengthen the overall impact of the research.

Comments on revision:

The authors have addressed my comments accordingly.

Reviewer #1 (Public Review):

Summary:

The emergence of Drosophila EM connectomes has revealed numerous neurons within the associative learning circuit. However, these neurons are inaccessible for functional assessment or genetic manipulation in the absence of cell-type-specific drivers. Addressing this knowledge gap, Shuai et al. have screened over 4000 split-GAL4 drivers and correlated them with identified neuron types from the "Hemibrain" EM connectome by matching light microscopy images to neuronal shapes defined by EM. They successfully generated over 800 split-GAL4 drivers and 22 split-LexA drivers covering a substantial number of neuron types across layers of the mushroom body associative learning circuit. They provide new labeling tools for olfactory and non-olfactory sensory inputs to the mushroom body; interneurons connected with dopaminergic neurons and/or mushroom body output neurons; potential reinforcement sensory neurons; and expanded coverage of intrinsic mushroom body neurons. Furthermore, the authors have optimized the GR64f-GAL4 driver into a sugar sensory neuron-specific split-GAL4 driver and functionally validated it as providing a robust optogenetic substitute for sugar reward. Additionally, a driver for putative nociceptive ascending neurons, potentially serving as optogenetic negative reinforcement, is characterized by optogenetic avoidance behavior. The authors also use their very large dataset of neuronal anatomies, covering many example neurons from many brains, to identify neuron instances with atypical morphology. They find many examples of mushroom body neurons with altered neuronal numbers or mistargeting of dendrites or axons and estimate that 1-3% of neurons in each brain may have anatomic peculiarities or malformations. Significantly, the study systematically assesses the individualized existence of MBON08 for the first time. This neuron is a variant shape that sometimes occurs instead of one of two copies of MBON09, and this variation is more common than that in other neuronal classes: 75% of hemispheres have two MBON09's, and 25% have one MBON09 and one MBON08. These newly developed drivers not only expand the repertoire for genetic manipulation of mushroom body-related neurons but also empower researchers to investigate the functions of circuit motifs identified from the connectomes. The authors generously make these flies available to the public. In the foreseeable future, the tools generated in this study will allow important advances in the understanding of learning and memory in Drosophila.

Strengths:

(1) After decades of dedicated research on the mushroom body, a consensus has been established that the release of dopamine from DANs modulates the weights of connections between KCs and MBONs. This process updates the association between sensory information and behavioral responses. However, understanding how the unconditioned stimulus is conveyed from sensory neurons to DANs, and the interactions of MBON outputs with innate responses to sensory context remains less clear due to the developmental and anatomic diversity of MBONs and DANs. Additionally, the recurrent connections between MBONs and DANs are reported to be critical for learning. The characterization of split-GAL4 drivers for 30 major interneurons connected with DANs and/or MBONs in this study will significantly contribute to our understanding of recurrent connections in mushroom body function.

(2) Optogenetic substitutes for real unconditioned stimuli (such as sugar taste or electric shock) are sometimes easier to implement in behavioral assays due to the spatial and temporal specificity with which optogenetic activation can be induced. GR64f-GAL4 has been widely used in the field to activate sugar sensory neurons and mimic sugar reward. However, the authors demonstrate that GR64f-GAL4 drives expression in other neurons not necessary for sugar reward, and the potential activation of these neurons could introduce confounds into training, impairing training efficiency. To address this issue, the authors have elaborated on a series of intersectional drivers with GR64f-GAL4 to dissect subsets of labeled neurons. This approach successfully identified a more specific sugar sensory neuron driver, SS87269, which consistently exhibited optimal training performance and triggered ethologically relevant local searching behaviors. This newly characterized line could serve as an optimized optogenetic tool for sugar reward in future studies.

(3) MBON08 was first reported by Aso et al. 2014, exhibiting dendritic arborization into both ipsilateral and contralateral γ3 compartments. However, this neuron could not be identified in the previously published Drosophila brain connectomes. In the present study, the existence of MBON08 is confirmed, occurring in one hemisphere of 35% of imaged flies. In brains where MBON08 is present, its dendrite arborization disjointly shares contralateral γ3 compartments with MBON09. This remarkable phenotype potentially serves as a valuable resource for understanding the stochasticity of neurodevelopment and the molecular mechanisms underlying mushroom body lobe compartment formation.

Reviewer #1 (Public review):

Summary:

Knudstrup et al. use two-photon calcium imaging to measure neural responses in the mouse primary visual cortex (V1) in response to image sequences. The authors presented mice with many repetitions of the same four-image sequence (ABCD) for four days. Then on the fifth day, they presented unexpected stimulus orderings where one stimulus was either omitted (ABBD) or substituted (ACBD). After analyzing trial-averaged responses of neurons pooled across multiple mice, they observed that stimulus omission (ABBD) caused a small, but significant, strengthening of neural responses but observed no significant change in the response to stimulus substitution (ACBD). Next, they performed population analyses of this dataset. They showed that there were changes in the correlation structure of activity and that many features about sequence ordering could be reliably decoded. This second set of analyses is interesting and exhibited larger effect sizes than the first results about predictive coding. However, concerns about the design of the experiment temper my enthusiasm.

The most recent version of this manuscript makes a few helpful changes (entirely in supplemental figures--the main text figures are unchanged). It does not resolve any of the larger weaknesses of the experimental design, or even perform single-neuron tracking in the one case where it was possible (between similar FOVs shown in Supplemental Figure 1).

Strengths:

(1) The topic of predictive coding in the visual cortex is exciting, and this task builds on previous important work by the senior author (Gavornik and Bear 2014) where unexpectedly shuffling sequence order caused changes in LFPs recorded from visual cortex.

(2) Deconvolved calcium responses were used appropriately here to look at the timing of the neural responses.

(3) Neural decoding results showing that the context of the stimuli could be reliably decoded from trial-averaged responses were interesting. But I have concerns about how the data was formatted for performing these analyses.

Weaknesses:

(1) All analyses were performed on trial-averaged neural responses that were pooled across mice (except for Supplementary Figure 6, see below). Owing to differences between subjects in behavior, experimental preparation quality, and biological variability, it seems important to perform most analyses on individual datasets to assess how behavioral training might differently affect each animal.

In the most recent draft, a single-mouse analysis was added for Figure 4C (Supplementary Figure 6). This effect of "representational drift" was not statistically quantified in either the single-mouse results or in the main text figure panel. Moreover, the apparent correlational drift could be accounted for by a reduction in SNR as a consequence of photobleaching.

(2) The correlation analyses presented in Figure 3 (labeled the second Figure 2 in the text) should be conducted on a single-animal basis. Studying population codes constructed by pooling across mice, particularly when there is no behavioral readout to assess whether learning has had similar effects on all animals, appears inappropriate to me. If the results in Figure 3 hold up on single animals, I think that is definitely an interesting result.

In the most recent draft, this analysis was still not performed on single mice. I was referring to the "decorrelation of responses" analysis in Figure 3, not the "representational drift" analysis in Figure 4. See my comments on Supplementary Figure 6 above.

(3) On Day 0 and Day 5, the reordered stimuli are presented in trial blocks where each image sequence is shown 100 times. Why wasn't the trial ordering randomized as was done in previous studies (e.g. Gavornik and Bear 2014)? Given this lack of reordering, did neurons show reduced predictive responses because the unexpected sequence was shown so many times in quick succession? This might change the results seen in Figure 2, as well as the decoder results where there is a neural encoding of sequence order (Figure 4). It would be interesting if the Figure 4 decoder stopped working when the higher order block structure of the task were disrupted.

In the rebuttal letter for the most recent draft, the authors refer to recent work in press (Hosmane et al. 2024) suggesting that because sleep may be important for plastic changes between sessions, they do not expect much change to be apparent within a session. However, they admit that this current study is too underpowered to know for sure--and do not cite or mention this yet unpublished work in the manuscript itself.

As a control, I would be interested to at least know how much variance in neural responses is observed between intermediate "training" sessions with identical stimuli, e.g. between Day 1 and Day 4, but this is not possible as imaging was not performed on these days.

Despite being referred to as "similar" I do not think early and late responses are clearly shown--aside from the histograms comparing "early traces" to "all traces" which include early traces in Figure 5B and Figure 6A. Showing variance in single-cell responses would be helpful to add in Supplementary Figure 3 and Supplementary Figure 4.

(4) A primary advantage of using two-photon calcium imaging over other techniques like extracellular electrophysiology is that the same neurons can be tracked over many days. This is a standard approach that can be accomplished by using many software packages-including Suite2P (Pachitariu et al. 2017), which is what the authors already used for the rest of their data preprocessing. The authors of this paper did not appear to do this. Instead, it appears that different neurons were imaged on Day 0 (baseline) and Day 5 (test). This is a significant weakness of the current dataset.

In the most recent draft, this concern has not been mitigated. Despite Supplementary Figure 1 showing similar FOVs, mostly different neurons were still extracted. In all other sessions, it is not reported how far apart the other recorded FOVs were from each other.

The rebuttal comment that the PE statistic is computed on an individual cell within-session basis is reasonable. Moreover, the bootstrapped version of the PE analysis in Supplementary Figure 8 is an improvement of the main analysis in the paper. As a control, it would have been helpful to compute the stability of the PE ratio statistics between training days (e.g. between day 1 and day 4). How much change would have been observed when none is expected? Unfortunately, imaging was not performed on these training days so this analysis will not be readily possible to perform. Moreover, the PE statistic requires averaging across cells and trials and is therefore very likely to wash out many interesting effects. Even if it is the population response that is changing, why would it be the arithmetic mean that changes in particular vs. some other projection of the population activity? The experimental and analysis design of the paper here remains weak in my mind.

Reviewer #2 (Public review):

Summary:

In this manuscript, the authors provide a method aiming to accurately reflect the individual deviation of longitudinal/temporal change compared to the normal temporal change characterized based on pre-trained population normative model (i.e., a Bayesian linear regression normative model), which was built based on cross-sectional data. This manuscript aims at solving a recently identified problem of using normative models based on cross-sectional data to make inferences about longitudinal change.

Strengths:

The efforts of this work make a good contribution to addressing an important question of normative modeling. With the greater availability of cross-sectional studies for normative modeling than longitudinal studies, and the inappropriateness of making inferences about longitudinal subject-specific changes using these cross-sectional data-based normative models, it's meaningful to try to address this gap from the aspect of methodological development.

Reviewer #1 (Public review):

Summary:

This work aims at understanding the role of thalamus POm in dorsal lateral striatum (DLS) projection in learning a sensorimotor associative task. The authors first confirm that POm forms "en passant" synapses with some of the DLS neuronal subtypes. They then perform a go/no-go associative task that consists of the mouse learning to discriminate between two different textures and to associate one of them with an action. During this task they either record the activity of the POm to DLS axons using endoscopy or silence their activity. They report that POm axons in the DLS are activated around the sensory stimulus but that the activity is not modulated by the reward. Last, they showed that silencing the POm axons at the level of DLS slows down learning the task.

The authors show convincing evidence of projections from POm to DLS and that POm inputs to DLS code for whisking whatever the outcome of the task is. However, their results do not allow to conclude if more neurones are recruited during the learning process or if the already activated fibres get activated more strongly. Last, because POm fibres in the DLS are also projecting to S1, silencing the POm fibres in the DLS could have affected inputs in S1 as well and therefore, the slowdown in acquiring the task is not necessarily specific to the POm to DLS pathway.

Strengths:

One of the main strengths of the paper is to go from slice electrophysiology to behaviour to get an in-depth characterization of one pathway. The authors did a comprehensive description of the POm projections to the DLS using transgenic mice to unambiguously identify the DLS neuronal population. They also used a carefully designed sensorimotor association task, and they exploited the results in depth.

It is a very nice effort to have measured the activity of the axons in the DLS not only after the mice have learned the task but throughout the learning process. It shows the progressive increase of activity of POm axons in the DLS, which could imply that there is a progressive strengthening of the pathway. The results show convincingly that POm axons in the DLS are not activated by the outcome of the task but by the whisker activity, and that this activity in average increases with learning.

Weaknesses:

One of the main targets of the striatum from thalamic input are the cholinergic neurons that weren't investigated here, is there information that could be provided?

It is interesting to know that the POm projects to all neuronal types in the DLS, but this information is not used further down the manuscript so the only take-home message of Figure 1 is that the axons that they image or silence in the DLS are indeed connected to DLS neurons and not just passing fibres. In this line, are these axons the same as the ones projecting to S1? If this is the case, why would we expect a different behaviour of the axon activity at the DLS level compared to S1?

The authors used endoscopy to measure the POm axons in the DLS activity, which makes it impossible to know if the progressive increase of POm response is due to an increase of activity from each individual neurons or if new neurons are progressively recruited in the process.

The picture presented in Figure 4 of the stimulation site is slightly concerning as there are hardly any fibres in neocortical layer 1 while there seems to be quite a lot of them in layer 4, suggesting that the animal here was injected in the VB. This is especially striking as the implantation and projection sites presented in Figure 1 and 2 are very clean and consistent with POm injection.

Comment after review: The weaknesses remain as concerns have not been addressed. The dataset is interesting but the interpretation, due partly to the lack of control (especially relative to VPM contamination), is difficult.

Reviewer #1 (Public review):

Summary:

In this series of studies, Locantore et al. investigated the role of SST-expressing neurons in the entopeduncular nucleus (EPNSst+) in probabilistic switching tasks, a paradigm that requires continued learning to guide future actions. In prior work, this group had demonstrated EPNSst+ neurons co-release both glutamate and GABA and project to the lateral habenula (LHb), and LHb activity is also necessary for outcome evaluation necessary for performance in probabilistic decision-making tasks. Previous slice physiology works have shown that the balance of glutamate/GABA co-release is plastic, altering the net effect of EPN on downstream brain areas and neural circuit function. The authors used a combination of in vivo calcium monitoring with fiber photometry and computational modelling to demonstrate that EPNSst+ neural activity represents movement, choice direction and reward outcomes in their behavioral task. However, viral-genetic manipulations to synaptically silence these neurons or selectively eliminate glutamate release had no effect on behavioral performance in well-trained animals. The authors conclude that despite their representation of task variables, EPN Sst+ neuron synaptic output is dispensable for task performance.

Strengths and Weaknesses:

Overall, the manuscript is exceptionally scholarly, with a clear articulation of the scientific question and a discussion of the findings and their limitations. The analyses and interpretations are careful and rigorous. This review appreciates the thorough explanation of the behavioral modelling and GLM for deconvolving the photometry signal around behavioral events, and the transparency and thoroughness of the analyses in the supplemental figures. This extra care has the result of increasing the accessibility for non-experts, and bolsters confidence in the results. To bolster a reader's understanding of results, we suggest it would be interesting to see the same mouse represented across panels (i.e. Fig 1 F-J, Supp 1 F,K etc i.e via inclusion of faint hash lines connecting individual data points across variables. Additionally, Fig 3E demonstrates that eliminating the 'reward' and 'choice and reward' terms from the GLM significantly worsens model performance; to demonstrate the magnitude of this effect, it would be interesting to include a reconstruction of the photometry signal after holding out of both or one of these terms, alongside the 'original' and 'reconstructed' photometry traces in panel D. This would help give context for how the model performance degrades by exclusion of those key terms. Finally, the authors claimed calcium activity increased following ipsilateral movements. However, figure 3C clearly shows that both SXcontra and SXisi increase beta coefficients. Instead, the choice direction may be represented in these neurons, given that beta coefficients increase following CXipsi and before SEipsi, presumably when animals make executive decisions. Could the authors clarify their interpretation on this point? Also, it is not clear if there is a photometry response related to motor parameters (i.e. head direction or locomotion, licking), which could change the interpretation of the reward outcome if it is related to a motor response; could the authors show photometry signal from representative 'high licking' or 'low licking' reward trials, or from spontaneous periods of high. Vs low locomotor speeds (if the sessions are recorded) to otherwise clarify this point?

There are a few limitations with the design and timing of the synaptic manipulations that would improve the manuscript if discussed or clarified. The authors take care to validate the intersectional genetic strategies: Tetanus Toxin virus (which eliminates synaptic vesicle fusion) or CRISPR editing of Slc17a6, which prevents glutamate loading into synaptic vesicles. The magnitude of effect in the slice physiology results are striking. However, this relies on co-infection of a second AAV to express channelrhodopsin for the purposes of validation, and it is surely the case that there will not be 100% overlap between the proportion of cells infected. Alternative means of glutamate packaging (other VGluT isoforms, other transporters etc) could also compensate for the partial absence of VGluT2, which should be discussed. The authors do not perform a complimentary experiment to delete GABA release (i.e. via VGAT editing), which is understandable, given the absence of an effect with the pan-synaptic manipulation. A more significant concern is the timing of these manipulations as the authors acknowledge. The manipulations are all done in well-trained animals, who continue to perform during the length of viral expression. Moreover, after carefully showing that mice use different strategies on the 70/30 version vs the 90/10 version of the task, only performance on the 90/10 version is assessed after the manipulation. Together, the observation that EPNsst activity does not alter performance on a well learned, 90/10 switching task decreases the impact of the findings, as this population may play a larger role during task acquisition or under more dynamic task conditions. Additional experiments could be done to strengthen the current evidence, although the limitations is transparently discussed by the authors.

Finally, intersectional strategies target LHb-projecting neurons, although in the original characterization it is not entirely clear that the LHb is the only projection target of EPNsst neurons. A projection map would help clarify this point.

Overall, the authors used a pertinent experimental paradigm and common cell-specific approaches to address a major gap in the field, which is the functional role of glutamate/GABA co-release from the major basal ganglia output nucleus in action selection and evaluation. The study is carefully conducted, their analyses are thorough, and the data are often convincing and thought-provoking. However, the limitations of their synaptic manipulations with respect to the behavioral assays reduces generalizability and to some extent the impact of their findings.

Comments on the latest version:

Specifically, they have included more thorough analyses to address several concerns related to interpreting activity patterns of EPSst+ neurons. The authors clearly point out that calcium activity increased during ipsilateral movements, and the increase was statistically larger during the choice phase (Figure 2 supplement 1F-G), indicating that these neurons may represent movement and additional factors (e.g. executive decision-making). Correspondingly, we appreciate the thorough explanation of using a GLM model to determine which behavioural variables contribute to observed physiological signals and adding the example reconstructed signal with direction and reward variables omitted in Figure 3 supplements 1 and 2.

Although no new manipulation experiment is added to the manuscript, the authors respond to common critiques related to testing the behavioural effect after the manipulations in well-trained mice. The discussion related to technical limitations, possible compensatory mechanisms and alternative interpretations is thorough and overall satisfying. Based on the behaviour modeling results, the authors speculate that animals need to integrate more evidence from the past to guide choice in a more uncertain environment (70/30 version), instead of adopting a 'win-stay, lose-shift' strategy in the more deterministic 90/10 version. The authors expand the discussion, but the possibility that EPNSst+ neurons contribute to task performance in well-trained animals under uncertainty is not directly tested. Along with other alternative explanations discussed in the manuscript, we think the paper is valuable literature for future studies to understand the basal ganglia circuits in learning and decision-making.

Reviewer #2 (Public review):

This study by Bell et al. focuses on understanding the roles of two alternatively spliced exons in the single Drosophila Cav2 gene cac. The authors generate a series of cac alleles in which one or the other mutually exclusive exons are deleted to determine the functional consequences at the neuromuscular junction. They find alternative splicing at one exon encoding part of the voltage sensor impacts the activation voltage as well as localization to the active zone. In contrast, splicing at the second exon pair does not impact Cav2 channel localization, but it appears to determine the abundance of the channel at active zones. Together, the authors propose that alternative splicing at the Cac locus enables diversity in Cav2 function generated through isoform diversity generated at the single Cav2 alpha subunit gene encoded in Drosophila.

Overall this is an excellent, rigorously validated study that defines unanticipated functions for alternative splicing in Cav2 channels. The authors have generated an important toolkit of mutually exclusive Cac splice isoforms that will be of broad utility for the field, and show convincing evidence for distinct consequences of alternative splicing of this single Cav2 channel at synapses. Importantly, the authors use electrophysiology and quantitative live sptPALM imaging to determine the impacts of Cac alternative splicing on synaptic function. There remain some questions regarding the mechanisms underlying the changes in Cac localization to somatodendritic compartments. Nonetheless, this is a compelling investigation of alternative splicing in Cav2 channels that should be of interest to many researchers.

Reviewer #1 (Public review):

Summary:

Authors benchmarked 5 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 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 detected by each of the methods. They also tracked the computational efficiency of these methods. The authors work is valuable for the tested species and the analyses presented appear to support their claim that users should be cautious calling IBD when SNP density is low and recombination rate is high.

Strengths:

The study design was solid. The authors set up their reasoning for using P. falciparum very well. The high recombination rate and similar mutation rate to human is indeed an interesting case. Further, they chose methods that were developed explicitly for each species. This was a strength of the work, as well as incorporating both simulated and empirical data to support their goal that IBD detection should be benchmarked in P. falciparum.

Weaknesses:

The scope of the optimization and application of results from the work are narrow, in that everything is fine-tuned for Plasmodium. Some of the results were not entirely unexpected for users of any of the tested software that was developed for humans. For example, it is known that Refined IBD is not going to do well with the combination of short IBD segments and low SNP density. Lastly, it appears the authors only did one large-scale simulation (there are no reported SDs).

Reviewer #1 (Public review):

Summary:

This Tanzanian study focused on the relationship between human genetic ancestry, Mycobacterium tuberculosis complex (MTBC) diversity, and tuberculosis (TB) disease severity. The authors analyzed the genetic ancestry of 1,444 TB patients and genotyped the corresponding MTBC strains isolated from the same individuals. They found that the study participants predominantly possess Bantu-speaking genetic ancestry, with minimal European and Asian ancestry. The MTBC strains identified were diverse and largely resulted from introductions from South or Central Asia. Unfortunately, no associations were identified between human genetic ancestry, the MTBC strains, or TB severity. The authors suggest that social and environmental factors are more likely to contribute to TB severity in this setting.

Strengths:

In comparison to other studies investigating the role of human genetics in TB phenotypes, this study is relatively large, with more than 1,400 participants.

The matched human-MTBC strain collection is valuable and offers the opportunity to address questions about human-bacterium co-evolution.

Weaknesses:

Although the authors had genome-wide genotyping and whole genome sequencing data, they only compared the associations between human ancestry and MTBC strains. Given the large sample size, they had the opportunity to conduct a genome-wide association study similar to that of Muller et al. (https://doi.org/10.1016/j.ygeno.2021.04.024).

The authors tested whether human genetic ancestry is associated with TB severity. However, the basis for this hypothesis is unclear. The studies cited as examples all focused on progression to active TB (from a latent infection state), which should not be conflated with disease severity. It is difficult to ascertain whether the role of genetic ancestry in disease severity would be detectable through this study design, as some participants might simply have been sicker for longer before being diagnosed (despite the inquiry about cough duration). This delay in diagnosis would not be influenced solely by human genetics, which is the conclusion of the study.

Additionally, the study only included participants who attended the TB clinic.

Including healthy controls from the general population would have provided an interesting comparison to see if ancestry proportions differ.

Although the authors suggest that social and environmental factors contribute to TB severity, only age, smoking, and HIV status were characterised in the study.

Reviewer #1 (Public review):

Summary:

This work introduces the differentiable Gillespie algorithm, DGA, which is a differentiable variant of the celebrated (and exact) Gillespie algorithm commonly used to perform stochastic simulations across numerous fields, notably in the life sciences. The proposed DGA approximates the exact Gillespie algorithm using smooth functions yielding a suitable approximate differentiable stochastic system as a proxy for the underlying discrete stochastic system, where DGA stochastic reactions have continuous reaction index and the species abundances. To illustrate their methodology, the authors specifically consider in detail the case of a well-studied two-state promoter gene regulation system that they analyze using a machine learning approach, and by combining simulation data with analytical results. For the two-state promoter gene system, the DGA is benchmarked by accurately reproducing the results of the exact Gillespie algorithm. For this same simple system, the authors also show that how the DGA can be used for estimating kinetic parameters of both simulated and real noisy experimental data. This let them argue convincingly that the DGA can become a powerful computation tool for applications in quantitative and synthetic biology. In order to argue that the DGA can be employed to design circuits with ad-hoc input-output relations, these considerations are then extended to a more complex four-state promoter model of gene regulation.

Strengths:

The main strength of the paper is its clarity and its pedagogical presentation of the simulation methods.

Weaknesses:

It would have been useful to have a brief discussion, based on a concrete example, of what can be achieved with the DGA and is totally beyond the reach of the Gillespie algorithm and the numerous existing stochastic simulation methods. A more comprehensive and quantitative analysis of the limitations of the DGA, e.g. for rare events, would have also been helpful.

Reviewer #1 (Public review):

Summary:

The authors provide a resource to the systems neuroscience community, by offering their Python-based CLoPy platform for closed-loop feedback training. In addition to using neural feedback, as is common in these experiments, they include a capability to use real-time movement extracted from DeepLabCut as the control signal. The methods and repository are detailed for those who wish to use this resource. Furthermore, they demonstrate the efficacy of their system through a series of mesoscale calcium imaging experiments. These experiments use a large number of cortical regions for the control signal in the neural feedback setup, while the movement feedback experiments are analyzed more extensively.

Strengths:

The primary strength of the paper is the availability of their CLoPy platform. Currently, most closed-loop operant conditioning experiments are custom built by each lab and carry a relatively large startup cost to get running. This platform lowers the barrier to entry for closed-loop operant conditioning experiments, in addition to making the experiments more accessible to those with less technical expertise.

Another strength of the paper is the use of many different cortical regions as control signals for the neurofeedback experiments. Rodent operant conditioning experiments typically record from the motor cortex and maybe one other region. Here, the authors demonstrate that mice can volitionally control many different cortical regions not limited to those previously studied, recording across many regions in the same experiment. This demonstrates the relative flexibility of modulating neural dynamics, including in non-motor regions.

Finally, adapting the closed-loop platform to use real-time movement as a control signal is a nice addition. Incorporating movement kinematics into operant conditioning experiments has been a challenge due to the increased technical difficulties of extracting real-time kinematic data from video data at a latency where it can be used as a control signal for operant conditioning. In this paper they demonstrate that the mice can learn the task using their forelimb position, at a rate that is quicker than the neurofeedback experiments.

Weaknesses:

There are several weaknesses in the paper that diminish the impact of its strengths. First, the value of the CLoPy platform is not clearly articulated to the systems neuroscience community. Similarly, the resource could be better positioned within the context of the broader open-source neuroscience community. For an example of how to better frame this resource in these contexts, I recommend consulting the pyControl paper. Improving this framing will likely increase the accessibility and interest of this paper to a less technical neuroscience audience, for instance by highlighting the types of experimental questions CLoPy can enable.

While the dataset contains an impressive amount of animals and cortical regions for the neurofeedback experiment, and an analysis of the movement-feedback experiments, my excitement for these experiments is tempered by the relative incompleteness of the dataset, as well as its description and analysis in the text. For instance, in the neurofeedback experiment, many of these regions only have data from a single mouse, limiting the conclusions that can be drawn. Additionally, there is a lack of reporting of the quantitative results in the text of the document, which is needed to better understand the degree of the results. Finally, the writing of the results section could use some work, as it currently reads more like a methods section.

Suggestions for improved or additional experiments, data or analyses:

Not necessary for this paper, but it would be interesting to see if the CLNF group could learn without auditory feedback.

There are no quantitative results in the results section. I would add important results to help the reader better interpret the data. For example, in: "Our results indicated that both training paradigms were able to lead mice to obtain a significantly larger number of rewards over time," You could show a number, with an appropriate comparison or statistical test, to demonstrate that learning was observed.

For: "Performing this analysis indicated that the Raspberry Pi system could provide reliable graded feedback within ~63 {plus minus} 15 ms for CLNF experiments." The LED test shows the sending of the signal, but the actual delay for the audio generation might be longer. This is also longer than the 50 ms mentioned in the abstract.

It could be helpful to visualize an individual trial for each experiment type, for instance how the audio frequency changes as movement speed / calcium activity changes.

The sample sizes are small (n=1) for a few groups. I am excited by the variety of regions recorded, so it could be beneficial for the authors to collect a few more animals to beef up the sample sizes.

I am curious as to why 60 trials sessions were used. Was it mostly for the convenience of a 30 min session, or were the animals getting satiated? If the former, would learning have occurred more rapidly with longer sessions?

Figure 4 E is interesting, it seems like the changes in the distribution of deltaF was in both positive and negative directions, instead of just positive. I'd be curious as to the author's thoughts as to why this is the case. Relatedly, I don't see Figure 4E, and a few other subplots, mentioned in the text. As a general comment, I would address each subplot in the text.

For: "In general, all ROIs assessed that encompassed sensory, pre-motor, and motor areas were capable of supporting increased reward rates over time," I would provide a visual summary showing the learning curves for the different types of regions.

Relatedly, I would further explain the fast vs slow learners, and if they mapped onto certain regions.

Also I would make the labels for these plots (e.g. Supp Fig3) more intuitive, versus the acronyms currently used.

The CLMF animals showed a decrease in latency across learning, what about the CLNF animals? There is currently no mention in the text or figures.

Reviewer #1 (Public review):

Summary:

Fluorescence imaging has become an increasingly popular technique for monitoring neuronal activity and neurotransmitter concentrations in the living brain. However, factors such as brain motion and changes in blood flow and oxygenation can introduce significant artifacts, particularly when activity-dependent signals are small. Yogesh et al. quantified these effects using GFP, an activity-independent marker, under two-photon and wide-field imaging conditions in awake behaving mice. They report significant GFP responses across various brain regions, layers, and behavioral contexts, with magnitudes comparable to those of commonly used activity sensors. These data highlight the need for robust control strategies and careful interpretation of fluorescence functional imaging data.

Strengths:

The effect of hemodynamic occlusion in two-photon imaging has been previously demonstrated in sparsely labeled neurons in V1 of anesthetized animals (see Shen and Kara et al., Nature Methods, 2012). The present study builds on these findings by imaging a substantially larger population of neurons in awake, behaving mice across multiple cortical regions, layers, and stimulus conditions. The experiments are extensive, the statistical analyses are rigorous, and the results convincingly demonstrate significant GFP responses that must be accounted for in functional imaging experiments. However, whether these GFP responses are driven by hemodynamic occlusion remains less clear, given the complexities associated with awake imaging and GFP's properties (see below).

Weaknesses:

(1) The authors primarily attribute the observed GFP responses to hemodynamic occlusion. While this explanation is plausible, other factors may also contribute to the observed signals. These include uncompensated brain movement (e.g., axial-direction movements), leakage of visual stimulation light into the microscope, and GFP's sensitivity to changes in intracellular pH (see e.g., Kneen and Verkman, 1998, Biophysical Journal). Although the correlation between GFP signals and blood vessel diameters supports a hemodynamic contribution, it does not rule out significant contributions from these (or other) factors. Consequently, whether GFP fluorescence can reliably quantify hemodynamic occlusion in two-photon microscopy remains uncertain.

(2) Regardless of the underlying mechanisms driving the GFP responses, these activity-independent signals must be accounted for in functional imaging experiments. However, the present manuscript does not explore potential strategies to mitigate these effects. Exploring and demonstrating even partial mitigation strategies could have significant implications for the field.

(3) Several methodology details are missing from the Methods section. These include: (a) signal extraction methods for two-photon imaging data (b) neuropil subtraction methods (whether they are performed and, if so, how) (c) methods used to prevent visual stimulation light from being detected by the two-photon imaging system (d) methods to measure blood vessel diameter/area in each frame. The authors should provide more details in their revision.

Reviewer #1 (Public review):

Summary:

This study aims to investigate the links between social behaviors observed in free-moving situations and behavioral performances measured in well-controlled, laboratory settings. The authors assessed general social tendencies and dyadic relationships among four monkeys in a group by scoring agonistic (aggression) and affiliative (grooming and proximity) behaviors in each pair. By measuring the saccadic reaction time in a classic social interference task, the authors reported that the monkeys with higher SEIs (i.e., more social individuals) were less distracted by the faces of other monkeys. These effects were enhanced when the distractors were out-group monkey faces rather than in-group ones. Lastly, oxytocin administration increased the impact of the out-group monkey faces in the social interference task, while reducing the magnitude of general social tendencies measured with SEI.

Strengths:

(1) The combination of behavioral data obtained in a colony room and in a laboratory environment is rare and important.<br /> (2) The evaluation of social interactions were successfully performed based on an automated target detection algorithm. The resulting multi-dimensional, complicated social interactions were summarized into simple indices (SEI and IEI). These indices provide a good measure for the social tendencies of each monkey.<br /> (3) Well-designed and robust experiments in the laboratory environment that are linked nicely with the general social tendencies observed in spontaneous behaviors.

Weaknesses:

(1) While the overall results are interesting, I am somewhat left confused about how to interpret the difference in the scores derived from different conditions. For example, the authors stated "Comparing the weights for in-group and out-group distractors, the effect of proximity was larger than that of aggression and grooming" in p.8. Does this mean that the proximity is indeed the type of behavior most affected in the out-group condition compared to the in-group condition? The out-group effects are difficult to examine with actual behavioral data, but some in-group effects such as those involving OT can be tested, which possibly provides good insights into interpreting the differences of the weights observed across the experimental conditions.

(2) I think it is important to provide how variable spontaneous social interactions were across sessions and how impactful the variability of the interactions is on the SEI and IEI, as it helps to understand how meaningful the differences of weights are across the conditions, but such data are missing. In line with this point, although the conclusions still hold as those data were obtained during the same experimental periods, shouldn't the weights in Fig. 3f and Figs. 4g and 4h (saline) be expected to be similar, if not the same?

Reviewer #1 (Public review):

Summary:

In this manuscript, Kondo et al. developed a method to suppress somatic action potentials while recording spine calcium signals using two-photon imaging in the L2/3 visual cortex in response to visual stimuli. The authors identified different patterns of dendritic spine activation by visual stimuli and analyzed how the different patterns of spine responses may contribute to somatic visual responses. Their analysis results suggest that spines on dendrites with a clustered arrangement can potentially generate sharply tuned output.

Strengths:

This is an interesting study addressing a standing question of how previously reported pepper-and-salt-like distributed sensory inputs on individual spines may give rise to somatic sensory selectivity. The method of somatic inhibition to prevent bAPs appears new and effective. The measurements of spine activity are carefully done. The finding that a small number of spines located in the same branch with similar tuning properties would predict the somatic tuning is consistent with local dendritic nonlinear integration mechanisms.

Weaknesses:

(1) The demonstration of the effectiveness of soma-specific inhibition is inadequate. Figure 1 only provides a single example trace showing the inhibition of somatic visual responses. The authors should provide statistical analysis over grouped data. For the effect of soma-specific inhibition on spine activity, the authors provided mostly negative results, lacking effects on spine responses for both soma inhibition and bAP subtraction. This is confusing. One possible explanation is that bAPs normally have little influence on spine activity. However, this would conflict with the known fact that somatic APs can easily invade spines in L2/3 neurons (e.g., Chen et al., Nature 2011). Another possibility is that under the current experimental conditions, somatic APs were rarely evoked by the visual stimulus. The authors should also rule out the possibility that the spines they imaged are from different neurons than the ones with somatic inhibition. The authors may consider identifying those cases where somatic APs have a significant impact on spine activity or spine tuning and show how bAP inhibition influences the dendritic and spine responses.

(2) Figure 4 shows that the proportion of spines with a preferred orientation similar to the soma (ΔOri {less than or equal to} 30{degree sign}) was 60%, which is surprisingly high. It is intriguing that without somatic AP invasion, there could be such a high degree of similarity between spine activity and somatic tuning. What is the ratio without soma inhibition? One could reason that with bAP invasion, there should be even more spines showing visual responses similar to those of the soma. Moreover, with such a high proportion of spines showing similar sensory tuning to the soma, it is inevitable that many branches contain more spines with similar tuning as the soma, exhibiting an apparent branch-specific clustering. While such apparent clustering may well predict somatic tuning, it primarily reflects a correlational relationship rather than a causal synaptic integration mechanism.

(3) There has been extensive work studying how the integration of spine activity or sub-branch activity gives rise to somatic output. The proposed main contribution of this study is to use an improved method to inhibit somatic activity in order to more confidently measure spine-specific activity and examine the integration mechanisms. However, the results showed that the measured spine-specific activity under soma inhibition was not significantly different from that measured under normal conditions (see point 1). It becomes unclear how this new method contributes to obtaining new insights into the synaptic integration mechanism.

(4) Figure 6 shows how the tuning similarity between spines depends on the distance between them. It is unclear what new information was acquired regarding the functional clustering of spines. This result can be largely explained by the overall higher proportion of similarly tuned spines (60%) compared to the soma's preferred orientations. Moreover, the authors did not demonstrate how such clustering may contribute to nonlinear synaptic integration.

(5) The results shown in Figure 7 can again be largely explained by the static property of a higher proportion of spines tuned similarly to the soma. These results do not reveal any active dendritic integration mechanisms.

Reviewer #1 (Public review):

In this manuscript, Chang et al. investigated the cell type-specific role of the integrin activator Shv in activity-dependent synaptic remodeling. Using the Drosophila larval neuromuscular junction as a model, they show that glial-secreted Shv modulates synaptic plasticity by maintaining the extracellular balance of neuronal Shv proteins and regulating ambient extracellular glutamate concentrations, which in turn affects postsynaptic glutamate receptor abundance. Furthermore, they report that genetic perturbation of glial morphogenesis phenocopies the defects observed with the loss of glial Shv. Altogether, their findings propose a role for glia in activity-induced synaptic remodeling through Shv secretion. While the conclusions are intriguing, several issues related to experimental design and data interpretation merit further discussion.

Reviewer #1 (Public review):

Summary:

As TDP-43 mislocalization is a hallmark of multiple neurodegenerative diseases, the authors seek to identify pathways that modulate TDP-43 levels. To do this, they use a FACS based genome wide CRISPR KD screen in a Halo tagged TDP-43 KI iPSC line. Their screen identifies a number of genetic modulators of TDP-43 expression including BORC which plays a role in lysosome transport.

Strengths:

Genome wide CRISPR based screen identifies a number of modulators of TDP-43 expression to generate hypotheses regarding RNA BP regulation and perhaps insights into disease.

Weaknesses:

It is unclear how altering TDP-43 levels may relate to disease where TDP-43 is not altered in expression but mislocalized. This is a solid cell biology study, but the relation to disease is not clear without providing evidence of BORC alterations in disease or manipulation of BORC reversing TDP-43 pathology in disease.

The mechanisms by which BORC and lysosome transport modulate TDP-43 expression are unclear. Presumably, this may be through altered degradation of TDP protein but this is not addressed.

Previous studies have demonstrated that TDP-43 levels can be modulated by altering lysosomal degradation so the identification of lysosomal pathways is not particularly novel.

It is unclear whether this finding is specific to TDP-43 levels or whether lysosome localization may more broadly impact proteostasis in particular of other RNA BPs linked to disease.

Unclear whether BORC depletion alters lysosome function or simply localization.

Reviewer #1 (Public review):

Summary of what is achieved: This manuscript validates and extends upon the sigh generating circuit between the NMB/GRP+ RTN/parafacial neurons and the NMBR/GRPR+ preBötC neurons established in Li et al., 2016. The authors generate multiple transgenic lines that enable selective targeting of these various sub-populations of cells and demonstrate the sufficiency of each type in generating a sigh breath. Additionally, they show that NMBR and GPRP preBötC neurons are glutamatergic, have overlapping and distinct expression, and do not express SST. Beyond this validation, the authors show that ectopic stimulation of SST neurons is sufficient to evoke sighs and that they are necessary for NMB/GRP induced sighing. This data is the first time that preBötC neurons downstream of NMBR/GRPR neurons have been identified that transform a eupneic breath into a sign breath. The five conclusions stated at the end of the introduction are supported by the data.

Summary of a primary weakness: A strong emphasis throughout the manuscript is the identification of an unsubstantiated slow sigh rhythm that is produced by NMBR/GRPR neurons. It is even suggested that this is an intrinsic property of these neurons. However, to make such a novel (and quite surprising) claim requires many more studies and the conclusion is dependent on how the authors have defined a sigh. Moreover, some data within the paper conflicts with this idea. The resubmitted manuscript does not contain any revisions and the rebuttal does not sufficiently address the critiques.

In summary, the optogenetic and chemogenetic characterization of the neuropeptide pathway transgenic lines nicely aligns with and provides important validation of the previous study by Li et. al., 2016 and the SST neuron studies provide a new mechanism for the transformation of NMBR/GRPR neuropeptide activation into a sigh. These are important findings, and they should be the points emphasized. The proposal of a slow sigh rhythm should be more rigorously established with new experiments and analysis or should be more carefully described and discussed.

Reviewer #1 (Public review):

Summary:

Dong et al here have studied the impact of the small Ras-like GTPase Rab10 on the exocytosis of dense core vesicles (DVC), which are important mediators of neuropeptide signaling in brain. They use optical imaging to show that lentiviral depletion of Rab10 in mouse hippocampal neurons in culture independent of the established defects in neurite outgrowth hamper DCV exocytosis. They further demonstrate that such defects are paralleled by changes in ER morphology and defective ER-based calcium buffering as well as reduced ribosomal protein expression in Rab10-depleted neurons. Re-expression of Rab10 or supplementation of exogenous L-leucine to restore defective neuronal protein synthesis rescues impaired DCV secretion. Based on these results they propose that Rab10 regulates DCV release by maintaining ER calcium homeostasis and neuronal protein synthesis.

Strengths:

This work provides interesting and potentially important new insights into the connection between ER function and the regulated secretion of neuropeptides via DCVs. The authors combine advanced optical imaging with light and electron microscopy, biochemistry and proteomics approaches to thoroughly assess the effects of Rab10 knockdown at the cellular level in primary neurons. The proteomic dataset provided may be valuable in facilitating future studies regarding Rab10 function. This work will thus be of interest to neuroscientists and cell biologists.

Weaknesses:

Whether and how the phenotypes of Rab10 reported in this study are linked remains an open question. Likewise, a possible role of Rab10 in exocytosis cannot be excluded at this stage.

Comments on revisions:

My previous questions and concerns have been satisfactorily addressed by the authors.

Reviewer #2 (Public review):

This work aggregates data across 5 openly available stopping studies (3 at 7 tesla and 2 at 3 tesla) to evaluate activity patterns across the common contrasts of Failed Stop (FS) > Go, FS > stop success (SS), and SS > Go. Previous work has implicated a set of regions that tend to be positively active in one or more of these contrasts, including the bilateral inferior frontal gyrus, preSMA, and multiple basal ganglia structures. However, the authors argue that upon closer examination, many previous papers have not found subcortical structures to be more active on SS than FS trials, bringing into question whether they play an essential role in (successful) inhibition. In order to evaluate this with more data and power, the authors aggregate across five datasets and find many areas that are *more* active for FS than SS, including bilateral preSMA, GPE, thalamus, and VTA. They argue that this brings into question the role of these areas in inhibition, based upon the assumption that areas involved in inhibition should be more active on successful stop than failed stop trials, not the opposite as they observed.

Comments on revisions:

The authors have been responsive to the feedback of both reviewers and they have significantly improved the manuscript. I now judge the work as valuable and solid. The authors have achieved their aims to characterize subcortical BOLD activation in the stop-signal paradigm.

Reviewer #1 (Public review):

The authors set out to define the molecular basis for LP as the origin of BRCA1-deficient breast cancers. They showed that LPs have the highest level of replicative stress, and hypothesise that this may account for their tendency to transform. They went on to identify ELF3 as a candidate driver of LP transformation and showed that ELF3 expression is up-regulated in response to replicative stress as well as BRCA1 deficiency. They went on to show that ELF3 inactivation led to a higher level of DNA damage, which may result from compromised replicative stress responses.

While the manuscript supports the interesting idea wherein ELF3 may fuel LP cell transformation, it remains obscure how ELF3 promotes cell tolerance to DNA damage. Interestingly the authors proposed that ELF3 suppresses excessive genomic instability, but in my opinion, I do not see any evidence that supports this claim. In fact, one might think that genomic instability is key to cell transformation.

Comments on revisions:

The authors have addressed most of my concerns.

This being said, the one major criticism raised by both Reviewers is the lack of evidence to support ELF3 as a driver of transformation of and in LP cells. The authors appear to have invested much resource and time but were not successful in isolating LP cells for experimentations. I would therefore suggest that the authors tone down their claims throughout the manuscript.

Reviewer #1 (Public review):

Summary:

This study analyzed biomarker data from 28 subjects with geographic atrophy (GA) in a Phase I/II clinical trial of PPY988, a subretinal AAV2 complement factor I (CFI) gene therapy, to evaluate pharmacokinetics and pharmacodynamics. Post-treatment, a 2-fold increase in vitreous humor (VH) FI was observed, correlating with a reduction in FB breakdown product Ba but minimal changes in other complement factors. The aqueous humor (AH) was found to be an unreliable proxy for VH in assessing complement activation. In vitro assays showed that the increase in FI had a minor effect on the complement amplification loop compared to the more potent C3 inhibitor pegcetacoplan. These findings suggest that PPY988 may not provide enough FI protein to effectively modulate complement activation and slow GA progression, highlighting the need for thorough biomarker review to determine optimal dosing in future studies.

Strengths:

This manuscript provides critical data on the efficacy of gene therapy for the eye, specifically introducing complement FI expression. It presents the results from a halted clinical trial, making the publication of this data essential for understanding the outcomes of this gene therapy approach. The findings offer valuable insights and lessons for future gene therapy attempts in similar contexts.

Weaknesses:

No particular weaknesses. The study was carefully performed and limitations are discussed.

I have just some concerns about the methodology used. The authors use the MILLIPLEX assays, which allow for multiplexed detection of complement proteins and they mention extensive validation. How are the measurements with this assay correlating with gold standard methods? Is the specificity and the expected normal ranges preserved with this assay? This also stands for the Olink assay. Some of the proteins are measured by both assay and/or by standard ELISA. How do these measurements correlate?

Comments on revisions:

The authors answered part of my comments. Only one remained - please provide a comparison between ELISA/Multiplex and Olink data to judge the robustness of the Olinkl assay for complement.

Reviewer #1 (Public review):

Summary:

Chlamydia spp. has a biphasic developmental cycle consisting of an extracellular, infectious form called an elementary body (EB) and an intracellular, replicative form known as a reticular body (RB). The structural stability of EBs is maintained by extensive cross linking of outer membrane proteins while the outer membrane proteins of RBs are in a reduced state. The overall redox state of EBs is more oxidized than RBs. The authors propose that redox state may be a controlling factor in the developmental cycle. To test this, alkyl hydroperoxide reductase subunit C (ahpC) was overexpressed or knocked down to examine effects on developmental gene expression. KD of ahpC induced increased expression of EB-specific genes and accelerated EB production. Conversely, overexpression of phpC delayed differentiation to EBs. The results suggest that chlamydial redox state may play a role in differentiation.

Strengths:

Uses modern genetic tools to explore the difficult area of temporal gene expression throughout the chlamydial developmental cycle.

Weaknesses:

The environmental signals triggering ahpC expression/activity are not determined.

Comments on revisions:

I am satisfied with the modifications made to the manuscript.

Reviewer #2 (Public review):

Summary:

In this manuscript the authors evaluate the attenuation, immunogenicity, and protection efficacy of a live-attenuated SARS-CoV-2 vaccine candidate (BK2102) against SARS-CoV-2.

Strengths:

The authors demonstrate that intranasal inoculation of BK2102 is safe and able to induce humoral and cellular immune responses in hamsters, without apparent signs of damage in the lungs, that protects against homologous SARS-CoV-2 and Omicron BA.5 challenge. Safety of BK2102 was further confirmed in a new hACE2 transgenic mouse model generated by the authors.

Weaknesses:

The authors have addressed my previous comments on the first submission of the document.

Reviewer #1 (Public Review):

Summary:

The present study's main aim is to investigate the mechanism of how VirR controls the magnitude of MEV release in Mtb. The authors used various techniques, including genetics, transcriptomics, proteomics, and ultrastructural and biochemical methods. Several observations were made to link VirR-mediated vesiculogenesis with PG metabolism, lipid metabolism, and cell wall permeability. Finally, the authors presented evidence of a direct physical interaction of VirR with the LCP proteins involved in linking PG with AG, providing clues that VirR might act as a scaffold for LCP proteins and remodel the cell wall of Mtb. Since the Mtb cell wall provides a formidable anatomical barrier for the entry of antibiotics, targeting VirR might weaken the permeability of the pathogen along with the stimulation of the immune system due to enhanced vesiculogenesis. Therefore, VirR could be an excellent drug target. Overall, the study is an essential area of TB biology.

Strengths:

The authors have done a commendable job of comprehensively examining the phenotypes associated with the VirR mutant using various techniques. Application of Cryo-EM technology confirmed increased thickness and altered arrangement of CM-L1 layer. The authors also confirmed that increased vesicle release in the mutant was not due to cell lysis, which contrasts with studies in other bacterial species.

Another strength of the manuscript is that biochemical experiments show altered permeability and PG turnover in the mutant, which fits with later experiments where authors provide evidence of a direct physical interaction of VirR with LCP proteins.

Transcriptomics and proteomics data were helpful in making connections with lipid metabolism, which the authors confirmed by analyzing the lipids and metabolites of the mutant.

Lastly, using three approaches, the authors confirm that VirR interacts with LCP proteins in Mtb via the LytR_C terminal domain.

Altogether, the work is comprehensive, experiments are designed well, and conclusions were made based on the data generated after verification using multiple complementary approaches.

Weaknesses:

The major weakness is that the mechanism of VirR-mediated EV release remains enigmatic. Most of the findings are observational and only associate enhanced vesiculogenesis observed in the VirR mutant with cell wall permeability and PG metabolism. Authors suggest that EV release occurs during cell division when PG is most fragile. However, this has yet to be tested in the manuscript-the AFM of the VirR mutant, which produces thicker PG with more pore density, displays enhanced vesiculogenesis. No evidence was presented to show that the PG of the mutant is fragile, and there are differences in cell division to explain increased vesiculogenesis. These observations, counterintuitive to the authors' hypothesis, need detailed experimental verification.

Transcriptomic data only adds a little substantial. Transcriptomic data do not correlate with the proteomics data. It remains unclear how VirR deregulates transcription. TLCs of lipids are not quantitative. For example, the TLC image of PDIM is poor; quantitative estimation needs metabolic labeling of lipids with radioactive precursors. Further, change in PDIMs is likely to affect other lipids (SL-1, PAT/DAT) that share a common precursor (propionyl- CoA).

The connection of cholesterol with cell wall permeability is tenuous. Cholesterol will serve as a carbon source and contribute to the biosynthesis of methyl-branched lipids such as PDIM, SL-1, and PAD/DAT. Carbon sources also affect other aspects of physiology (redox, respiration, ATP), which can directly affect permeability and import/export of drugs. Authors should investigate whether restoration of the normal level of permeability and EV release is not due to the maintenance of cell wall lipid balance upon cholesterol exposure of the VirR mutant.

Finally, protein interaction data is based on experiments done once without statistical analysis. If the interaction between VirR and LCP protein is expected on the mycobacterial membrane, how SPLIT_GFP system expressed in the cytoplasm is physiologically relevant. No explanation was provided as to why VirR interacts with the truncated version of LCP proteins and not with the full-length proteins.

Comments on revisions:

The authors have addressed my comments. I have no further issues.

Reviewer #1 (Public review):

Summary:

In this preprint, the authors systematically and rigorously investigate how specific classes of residue mutations alter the critical temperature as a proxy for the driving forces for phase separation. The work is well executed, the manuscript well-written, and the results reasonable and insightful.

Strengths:

The introductory material does an excellent job of being precise in language and ideas while summarizing the state of the art. The simulation design, execution, and analysis are exceptional and set the standard for large-scale simulation studies. The results, interpretations, and Discussion are largely nuanced, clear, and well-motivated, and the pedagogical nature with which sampling convergence is discussed is greatly appreciated. Finally, the underlying data are shared in a clear and accessible manner. Overall, the manuscript is a model

Weaknesses:

The simplicity of a one-bead-per-residue model parameterized to capture UCST-type phase behavior does perhaps impact some aspects of the generality of this work. That said, the authors carefully acknowledge these limitations, and overall, this is not seen as a major weakness of the conclusions drawn or the manuscript, given those conclusions are appropriately couched.

Joint Public Review:

This manuscript by Tao et al. reports on an effort to better specify the underlying interactions driving the effects of biodiversity on productivity in biodiversity experiments. The authors are especially concerned with the potential for competitive interactions to drive positive biodiversity-ecosystem functioning relationships by driving down the biomass of subdominant species. The authors suggest a new partitioning schema that utilizes a suite of partial density treatments to capture so-called competitive ability.

Readers are encouraged to consider the original reviews in full, which outline the strengths and weaknesses of the work:

First version: https://elifesciences.org/reviewed-preprints/98073v1/reviews

Second version: https://elifesciences.org/reviewed-preprints/98073v2/reviews

There are no further reviews for this version because the authors declined to make further improvements to their manuscript.

Reviewer #1 (Public review):

Summary

In this human neuroimaging and electrophysiology study, the authors aimed to characterise effects of a period of visual deprivation in the sensitive period on excitatory and inhibitory balance in the visual cortex. They attempted to do so by comparing neurochemistry conditions ('eyes open', 'eyes closed') and resting state, and visually evoked EEG activity between ten congenital cataract patients with recovered sight (CC), and ten age-matched control participants (SC) with normal sight.

First, they used magnetic resonance spectroscopy to measure in vivo neurochemistry from two locations, the primary location of interest in the visual cortex, and a control location in the frontal cortex. Such voxels are used to provide a control for the spatial specificity of any effects, because the single-voxel MRS method provides a single sampling location. Using MR-visible proxies of excitatory and inhibitory neurotransmission, Glx and GABA+ respectively, the authors report no group effects in GABA+ or Glx, no difference in the functional conditions 'eyes closed' and 'eyes open'. They found an effect of group in the ratio of Glx/GABA+ and no similar effect in the control voxel location. They then perform multiple exploratory correlations between MRS measures and visual acuity, and report a weak positive correlation between the 'eyes open' condition and visual acuity in CC participants.

The same participants then took part in an EEG experiment. The authors selected two electrodes placed in the visual cortex for analysis and report a group difference in an EEG index of neural activity, the aperiodic intercept, as well as the aperiodic slope, considered a proxy for cortical inhibition. Control electrodes in the frontal region did not present with the same pattern. They report an exploratory correlation between the aperiodic intercept and Glx in one out of three EEG conditions.

The authors report the difference in E/I ratio, and interpret the lower E/I ratio as representing an adaptation to visual deprivation, which would have initially caused a higher E/I ratio. Although intriguing, the strength of evidence in support of this view is not strong. Amongst the limitations are the low sample size, a critical control cohort that could provide evidence for higher E/I ratio in CC patients without recovered sight for example, and lower data quality in the control voxel. Nevertheless, the study provides a rare and valuable insight into experience-dependent plasticity in the human brain.

Strengths of study

How sensitive period experience shapes the developing brain is an enduring and important question in neuroscience. This question has been particularly difficult to investigate in humans. The authors recruited a small number of sight-recovered participants with bilateral congenital cataracts to investigate the effect of sensitive period deprivation on the balance of excitation and inhibition in the visual brain using measures of brain chemistry and brain electrophysiology. The research is novel, and the paper was interesting and well written.

Limitations

Low sample size. Ten for CC and ten for SC, and further two SC participants were rejected due to lack of frontal control voxel data. The sample size limits the statistical power of the dataset and increases the likelihood of effect inflation.

In the updated manuscript, the authors have provided justification for their sample size by pointing to prior studies and the inherent difficulties in recruiting individuals with bilateral congenital cataracts. Importantly, this highlights the value the study brings to the field while also acknowledging the need to replicate the effects in a larger cohort.

Lack of specific control cohort. The control cohort has normal vision. The control cohort is not specific enough to distinguish between people with sight loss due to different causes and patients with congenital cataracts with co-morbidities. Further data from a more specific populations, such as patients whose cataracts have not been removed, with developmental cataracts, or congenitally blind participants, would greatly improve the interpretability of the main finding. The lack of a more specific control cohort is a major caveat that limits a conclusive interpretation of the results.

In the updated version, the authors have indicated that future studies can pursue comparisons between congenital cataract participants and cohorts with later sight loss.

MRS data quality differences. Data quality in the control voxel appears worse than in the visual cortex voxel. The frontal cortex MRS spectrum shows far broader linewidth than the visual cortex (Supplementary Figures). Compared to the visual voxel, the frontal cortex voxel has less defined Glx and GABA+ peaks; lower GABA+ and Glx concentrations, lower NAA SNR values; lower NAA concentrations. If the data quality is a lot worse in the FC, then small effects may not be detectable.

In the updated version, the authors have added more information that informs the reader of the MRS quality differences between voxel locations. This increases the transparency of their reporting and enhances the assessment of the results.

Because of the direction of the difference in E/I, the authors interpret their findings as representing signatures of sight improvement after surgery without further evidence, either within the study or from the literature. However, the literature suggests that plasticity and visual deprivation drives the E/I index up rather than down. Decreasing GABA+ is thought to facilitate experience dependent remodelling. What evidence is there that cortical inhibition increases in response to a visual cortex that is over-sensitised to due congenital cataracts? Without further experimental or literature support this interpretation remains very speculative.

The updated manuscript contains key reference from non-human work to justify their interpretation.

Heterogeneity in patient group. Congenital cataract (CC) patients experienced a variety of duration of visual impairment and were of different ages. They presented with co-morbidities (absorbed lens, strabismus, nystagmus). Strabismus has been associated with abnormalities in GABAergic inhibition in the visual cortex. The possible interactions with residual vision and confounds of co-morbidities are not experimentally controlled for in the correlations, and not discussed.

The updated document has addressed this caveat.

Multiple exploratory correlations were performed to relate MRS measures to visual acuity (shown in Supplementary Materials), and only specific ones shown in the main document. The authors describe the analysis as exploratory in the 'Methods' section. Furthermore, the correlation between visual acuity and E/I metric is weak, not corrected for multiple comparisons. The results should be presented as preliminary, as no strong conclusions can be made from them. They can provide a hypothesis to test in a future study.

This has now been done throughout the document and increases the transparency of the reporting.

P.16 Given the correlation of the aperiodic intercept with age ("Age negatively correlated with the aperiodic intercept across CC and SC individuals, that is, a flattening of the intercept was observed with age"), age needs to be controlled for in the correlation between neurochemistry and the aperiodic intercept. Glx has also been shown to negatively correlates with age.

This caveat has been addressed in the revised manuscript.

Multiple exploratory correlations were performed to relate MRS to EEG measures (shown in Supplementary Materials), and only specific ones shown in the main document. Given the multiple measures from the MRS, the correlations with the EEG measures were exploratory, as stated in the text, p.16, and in Fig.4. yet the introduction said that there was a prior hypothesis "We further hypothesized that neurotransmitter changes would relate to changes in the slope and intercept of the EEG aperiodic activity in the same subjects." It would be great if the text could be revised for consistency and the analysis described as exploratory.

This has been done throughout the document and increases the transparency of the reporting.

The analysis for the EEG needs to take more advantage of the available data. As far as I understand, only two electrodes were used, yet far more were available as seen in their previous study (Ossandon et al., 2023). The spatial specificity is not established. The authors could use the frontal cortex electrode (FP1, FP2) signals as a control for spatial specificity in the group effects, or even better, all available electrodes and correct for multiple comparisons. Furthermore, they could use the aperiodic intercept vs Glx in SC to evaluate the specificity of the correlation to CC.

This caveat has been addressed. The authors have added frontal electrodes to their analysis, providing an essential regional control for the visual cortex location.

Comments on the latest version:

The authors have made reasonable adjustments to their manuscript that addressed most of my comments by adding further justification for their methodology, essential literature support, pointing out exploratory analyses, limitations and adding key control analyses. Their revised manuscript has overall improved, providing valuable information, though the evidence that supports their claims is still incomplete.

Reviewer #1 (Public review):

Contractile Injection Systems (CIS) are versatile machines that can form pores in membranes or deliver effectors. They can act extra or intracellularly. When intracellular they are positioned to face the exterior of the cell and hence should be anchored to the cell envelope. The authors previously reported the characterization of a CIS in Streptomyces coelicolor, including significant information on the architecture of the apparatus. However, how the tubular structure is attached to the envelope was not investigated. Here they provide a wealth of evidence to demonstrate that a specific gene within the CIS gene cluster, cisA, encodes a membrane protein that anchors the CIS to the envelope. More specifically, they show that:

- CisA is not required for assembly of the structure but is important for proper contraction and CIS-mediated cell death<br /> - CisA is associated to the membrane (fluorescence microscopy, cell fractionation) through a transmembrane segment (lacZ-phoA topology fusions in E. coli)<br /> - Structural prediction of interaction between CisA and a CIS baseplate component<br /> - In addition they provide a high-resolution model structure of the >750-polypeptide Streptomyces CIS in its extended conformation, revealing new details of this fascinating machine, notably in the baseplate and cap complexes.

All the experiments are well controlled including trans-complemented of all tested phenotypes.

One important information we miss is the oligomeric state of CisA.

While it would have been great to test the interaction between CisA and Cis11, to perform cryo-electron microscopy assays of detergent-extracted CIS structures to maintain the interaction with CisA, I believe that the toxicity of CisA upon overexpression or upon expression in E. coli render these studies difficult and will require a significant amount of time and optimization to be performed. It is worth mentioning that this study is of significant novelty in the CIS field because, except for Type VI secretion systems, very few membrane proteins or complexes responsible for CIS attachment have been identified and studied.

Reviewer #1 (Public review):

Summary:

Carter et al. present the eduWOSM imaging platform, a promising development in open-source microscopy for educational purposes. The paper outlines the construction and setup of this versatile microscope, demonstrating its capabilities through three key examples: single fluorophore tracking of tubulin heterodimers in gliding microtubules, 4D deconvolution imaging and tracking of chromosome movements in dividing human cells, and automated single-particle tracking in vitro and in live cells, with motion classified into sub-diffusive, diffusive, and super-diffusive behaviors.

The paper is well-written and could be strengthened by providing more empirical data on its performance, addressing potential limitations, and offering detailed insights into its educational impact. The project holds great potential and more discussion on long-term support and broader applications would provide a more comprehensive view of its relevance in different contexts.

Strengths:

(1) The eduWOSM addresses a crucial need in education, providing research-quality imaging at a lower cost (<$10k). The fact that it is open-source adds significant value, enabling broad accessibility even in under resourced areas.<br /> (2) There is availability of extensive resources, including a dedicated website, YouTube channel, and comprehensive tutorial guides to help users replicate the microscope.<br /> (3) The compact, portable, and stable design makes it easy to build multiple systems for use in diverse environments, including crowded labs and classrooms. This is further enhanced by the fact multiple kind of imaging experiments can be run on the system, from live imaging to super-resolution imaging.<br /> (4) The paper highlights the user-friendly nature of the platform, with the imaging examples in the paper being acquired by undergrad students.

Weaknesses:

(1) The paper mentions the microscope is suitable not just for education but even for research purposes. This claim needs validation through quantitative comparison to existing research-grade microscopes in terms of resolution, signal-to-noise ratio, and other key metrics. Adding more rigorous comparisons would solidify its credibility for research use, which would immensely increase the potential of the microscope.<br /> (2) The open-source microscope field is crowded with various options catering to hobby, educational, and research purposes (e.g., openFLexure, Flamingo, Octopi, etc.). The paper would benefit from discussing whether any aspects set the eduWOSM platform apart or fulfill specific roles that other microscopes do not.<br /> (3) While the eduWOSM platform is designed to be user-friendly, the paper would benefit from discussing whether the microscope can be successfully built and operated by users without direct help from the authors. It's important to know if someone with basic technical knowledge, relying solely on the provided resources (website, YouTube tutorials, and documentation), can independently assemble, calibrate, and operate the eduWOSM.<br /> (4) Ensuring long-term support and maintenance of the platform is crucial. The paper would benefit from addressing how the eduWOSM developers plan to support updates, improvements, or troubleshooting.

Reviewer #1 (Public review):

Summary:

The authors show that early life experience of juvenile bats shape their outdoor foraging behaviors. They achieve this by raising juvenile bats either in an impoverished or enriched environment. They subsequently test the behavior of bats indoors and outdoors. The authors show that behavioral measures outdoors were more reliable in delineating the effect of early life experiences as the bats raised in enriched environments were more bold, active and exhibit higher exploratory tendencies.

Strengths:

The major strength of the study is providing a quantitative study of animal "personality" and how it is likely shaped by innate and environmental conditions. The other major strength is the ability to do reliable long term recording of bats in the outdoors giving researchers the opportunity to study bats in their natural habitat. To this point, the study also shows that the behavioral variables measured indoors do not correlate to that measured outdoors, thus providing a key insight into the importance of testing animal behaviors in their natural habitat.

Weaknesses:

It is not clear from the analysis presented in the paper how persistent those environmentally induced changes, do they remain with the bats till the end of their lives.

Reviewer #1 (Public review):

Summary:

The manuscript aims to elucidate the impact of a prophage within the genome of Shewanella fidelis on its interaction with the marine tunicate Ciona robusta. The authors made a deletion mutant of S. fidelis that lacks one of its two prophages. This mutant exhibited an enhanced biofilm phenotype, as assessed through crystal violet staining, and showed reduced motility. The authors examined the effect of prophage deletion on several genes that could modulate cyclic-diGMP levels. While no significant changes were observed under in vitro conditions, the gene for one protein potentially involved in cyclic-diGMP hydrolysis was overexpressed during microbe-host interactions. The mutant was retained more effectively within a one-hour timeframe, whereas the wild-type (WT) strain became more abundant after 24 hours. Fluorescence microscopy was used to visualize the localization patterns of the two strains, which appeared to differ. Additionally, a significant difference in the expression of one immune protein was noted after one hour, but this difference was not evident after 23 hours. An effect of VCBC-C addition on the expression of one prophage gene was also observed.

Strengths:

I appreciate how the authors integrate diverse expertise and methods to address questions regarding the impact of prophages on gut microbiome-host interactions. The chosen model system is appropriate, as it allows for high-throughput experimentation and the application of simple imaging techniques.

Weaknesses:

My primary concern is that the manuscript primarily describes observations without providing insight into the molecular mechanisms underlying the observed differences. It is particularly unclear how the presence of the prophage leads to the phenotypic changes related to bacterial physiology and host-microbe interactions. Which specific prophage genes are critical, or is the insertion at a specific site in the bacterial genome the key factor? While significant effects on bacterial physiology are reported under in vitro conditions, there is no clear attribution to particular enzymes or proteins. In contrast, when the system is expanded to include the tunicate, differences in the expression of a cyclic-diGMP hydrolase become apparent. Why do we not observe such differences under in vitro conditions, despite noting variations in biofilm formation and motility? Furthermore, given that the bacterial strain possesses two prophages, I am curious as to why the authors chose to target only one and not both.

Regarding the microbe-host interaction, it is not clear why the increased retention ability of the prophage deletion strain did not lead to greater cell retention after 24 hours, especially since no differences in the immune response were observed at that time point.

Concerning the methodological approach, I am puzzled as to why the authors opted for qPCR instead of transcriptomics or proteomics. The latter approaches could have provided a broader understanding of the prophage's impact on both the microbe and the host.

Reviewer #1 (Public review):

Summary:

The authors showed the presence of Mtb in human liver biopsy samples of TB patients and reported that chronic infection of Mtb causes immune-metabolic dysregulation. Authors showed that Mtb replicates in hepatocytes in a lipid rich environment created by up regulating transcription factor PPARγ. Authors also reported that Mtb protects itself from anti-TB drugs by inducing drug metabolising enzymes.

Strengths:

It has been shown that Mtb induces storage of triacylglycerol in macrophages by induction of WNT6/ACC2 which helps in its replication and intracellular survival, however, creation of favorable replicative niche in hepatocytes by Mtb is not reported. It is known that Mtb infects macrophages and induces formation of lipid-laden foamy macrophages which eventually causes tissue destruction in TB patients. In a recent article it has been reported that "A terpene nucleoside from M. tuberculosis induces lysosomal lipid storage in foamy macrophages" that shows how Mtb manipulates host defense mechanisms for its survival. In this manuscript, authors reported the enhancement of lipid droplets in Mtb infected hepatocytes and convincingly showed that fatty acid synthesis and triacylglycerol formation is important for growth of Mtb in hepatocytes. The authors also showed the molecular mechanism for accumulation of lipid and showed that the transcription factor associated with lipid biogenesis, PPARγ and adipogenic genes were upregulated in Mtb infected cells.

The comparison of gene expression data between macrophages and hepatocytes by authors is important which indicates that Mtb modulates different pathways in different cell type as in macrophages it is related to immune response whereas, in hepatocytes it is related to metabolic pathways.

Authors also reported that Mtb residing in hepatocytes showed drug tolerance phenotype due to up regulation of enzymes involved in drug metabolism and showed that cytochrome P450 monooxygenase that metabolize rifampicin and NAT2 gene responsible for N-acetylation of isoniazid were up regulated in Mtb infected cells.

Weaknesses:

There are reports of hepatic tuberculosis in pulmonary TB patients especially in immune-compromised patients, therefore finding granuloma in human liver biopsy samples is not surprising.<br /> Mtb infected hepatic cells showed induced DME and NAT and this could lead to enhanced metabolism of drug by hepatic cells as a result Mtb in side HepG2 cells get exposed to reduced drug concentration and show higher tolerance to drug. The authors mentioned that " hepatocyte resident Mtb may display higher tolerance to rifampicin". In my opinion higher tolerance to drugs is possible only when DME of Mtb inside is up regulated or the target is modified. Although, in the end authors mentioned that drug tolerance phenotype can be better attributed to host intrinsic factors rather than Mtb efflux pumps. It may be better if the Drug tolerant phenotype section can be rewritten to clarify the facts.

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<br /> (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.<br /> (3) It is not entirely clear why specific length distribution for VCG oligomers has to be assumed rather than emerged from simulations.<br /> (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.

Weaknesses:

The conclusions are somewhat convoluted and could be presented better.

Reviewer #1 (Public review):

Summary:

Persistence is a phenomenon by which genetically susceptible cells are able to survive exposure to high concentrations of antibiotics. This is especially a major problem when treating infections caused by slow growing mycobacteria such as M. tuberculosis and M. abscessus. Studies on the mechanisms adopted by the persisting bacteria to survive and evade antibiotic killing can potentially lead to faster and more effective treatment strategies.

To address this, in this study, the authors have used a transposon mutagenesis based sequencing approach to identify the genetic determinants of antibiotic persistence in M. abscessus. To enrich for persisters they employed conditions, that have been reported previously to increase persister frequency - nutrient starvation, to facilitate genetic screening for this phenotype. M.abs transposon library was grown in nutrient rich or nutrient depleted conditions and exposed to TIG/LZD for 6 days, following which Tn-seq was carried out to identify genes involved in spontaneous (nutrient rich) or starvation-induced conditions. About 60% of the persistence hits were required in both the conditions. Pathway analysis revealed enrichment for genes involved in detoxification of nitrosative, oxidative, DNA damage and proteostasis stress. The authors then decided to validate the findings by constructing deletions of 5 different targets (pafA, katG, recR, blaR, Mab_1456c) and tested the persistence phenotype of these strains. Rather surprisingly only 2 of the 5 hits (katG and pafA) exhibited a persistence defect when compared to wild type upon exposure to TIG/LZD and this was complemented using an integrative construct. The authors then investigated the specificity of delta-katG susceptibility against different antibiotic classes and demonstrated increased killing by rifabutin. The katG phenotype was shown to be mediated through the production of oxidative stress which was reverted when the bacterial cells were cultured under hypoxic conditions. Interestingly, when testing the role of katG in other clinical strains of Mab, the phenotype was observed only in one of the clinical strains demonstrating that there might be alternative anti-oxidative stress defense mechanisms operating in some clinical strains.

Strengths:

While the role of ROS in antibiotic mediated killing of mycobacterial cells have been studied to some extent, this paper presents some new findings with regards to genetic analysis of M. abscessus susceptibility, especially against clinically used antibiotics, which makes it useful. Also, the attempts to validate their observations in clinical isolates is appreciated.

Weaknesses:

- Fig. 3 - 5 of the hits from the transposon screen were reconstructed as clean deletion strains and tested for persistence. However, only 1 (katG) gave a strong and 1 (Mab_1456c) exhibited a minor defect. Two of the clones did not show any persistence phenotype (blaR and recR) and one (pafA) showed a minor phenotype, however it was not clear if this difference was really relevant as the mutant exhibited differences at Day 0, prior to the addition of antibiotics. Considering these results from the validation, the conclusion would be that the Tn-seq approach to screen persistence defects is not reliable and is more likely to result in misses than hits.

- Fig 3 - Why is there such a huge difference in the extent of killing of the control strain in media, when exposed to TIG/LZD, when compared to Fig. 1C and Fig. 4. In Fig. 1C, M. abs grown in media decreases by >1 log by Day 3 and >4 log by Day 6, whereas in Fig. 3, the bacterial load decreases by <1 log by Day 3 and <2 log by Day 6. This needs to be clarified, if the experimental conditions were different, because if comparing to Fig. 1C data then the katG mutant strain phenotype is not very different.

Reviewer #1 (Public review):

Summary:

The imaging pipeline presented in this paper is a useful tool for visualizing and dynamically tracking bacterial colony formation at the individual worm level, enabling the study of microbiome colonization's association with host physiology, including lifespan, infection severity, and genetic mutations in real-time. This technique allows for certain biological information to be obtained that was previously missed such as pmk-1 mutants exhibiting a higher rate of colonization by E. coli OP50 than wild-type animals. Overall, this platform could be of interest to many labs studying C. elegans interactions with their microbiome and with bacterial pathogens.

Strengths:

This platform allows for unbiased quantifications of microbe colonization of bacteria at scale. This is particularly important in a field studying dynamic responses or potentially more subtle or variable phenotypes.

Platform could be adapted for multiple uses or potentially other species of nematodes for evolutionary comparisons.

The platform allows researchers to correlate bacterial colonization with predicted lifespan.

Weaknesses:

Platform will require optimization for any given bacteria species which restricts its ease of use for researchers that won't regularly be studying the same bacteria.

Requires the bacteria to be genetically tractable so cannot be easily adapted to microbes that do not have established ways of expressing GFP or other reporters.

This platform requires the use of relatively older adult animals that are more prone to larger gut colonies of bacteria. Thus, studies using this platform are restricted to studying older populations.

The relationship between bacterial colonization and host lifespan requires further investigation. The current SICKO platform and experimentation cannot fully address whether animals in poorer health are more susceptible to colonization, or whether colonization casually contributes to a decline in health. Furthermore, while such effects are statistically significant their effect size in some cases is modest.

Reviewer #3 (Public Review):

In multiple cancers, the key roles of B cells are emerging in the tumor microenvironment (TME). The authors of this study appropriately introduce that B cells are relatively under-characterised in the TME and argue correctly that it is not known how the B cell receptor (BCR) repertoires across tumor, lymph node and peripheral blood relate. The authors therefore supply a potentially useful study evaluating the tumor, lymph node and peripheral blood BCR repertoires and site-to-site as well as intra-site relationships. The authors employ sophisticated analysis techniques, although the description of the methods is incomplete.

Major strengths:

(1) The authors provide a unique analysis of BCR repertoires across tumor, dLN, and peripheral blood. The work provides useful insights into inter- and intra-site BCR repertoire heterogeneity. While patient-to-patient variation is expected, the findings with regard to intra-tumor and intra-dLN heterogeneity with the use of fragments from the same tissue are of importance, contribute to the understanding of the TME, and will inform future study design.

(2) A particular strength of the study is the detailed CDR3 physicochemical properties analysis which leads the authors to observations that suggest a less-specific BCR repertoire of TIL-B compared to circulating B cells.

Comments on revisions:

Your efforts in addressing concerns related to methodological details, narrative clarity, and data representation are commendable. The expanded descriptions of Fig. 1A and the experimental design, as well as the restructuring of the discussion, have greatly enhanced the manuscript's clarity and coherence.

Reviewer #1 (Public review):

Summary:

In this work, a screening platform is presented for rapid and cost-effective screening of candidate genes involved in Fragile Bone Disorders. The authors validate the approach of using crispants, generating FO mosaic mutants, to evaluate the function of specific target genes in this particular condition. The design of the guide RNAs is convincingly described, while the effectiveness of the method is evaluated to 60% to 92% of the respective target genes being presumably inactivated. Thus, injected F0 larvae can be directly used to investigate the consequences of this inactivation.

Skeletal formation is then evaluated at 7dpf and 14dpf, first using a transgenic reporter line revealing fluorescent osteoblasts, second using alizarin-red staining of mineralized structures. In general, it appears that the osteoblast-positive areas are more often affected in the crispants compared to the mineralized areas, an observation that appears to correlate with the observed reduced expression of bglap, a marker for mature osteoblasts, and the increased expression of col1a1a in more immature osteoblasts.

Finally, the injected fish (except two lines that revealed a high mortality) are also analyzed at 90dpf, using alizarin red staining and micro-CT analysis, revealing an increased incidence of skeletal deformities in the vertebral arches, fractures, as well as vertebral fusions and compressions for all crispants except those for daam2. Finally, the Tissue Mineral Density (TMD) as determined by micro-CT is proposed as an important marker for investigating genes involved in osteoporosis.<br /> Taken together, this manuscript is well presented, the data are clear and well analyzed, and the methods well described. It makes a compelling case for using the crispant technology to screen the function of candidate genes in a specific condition, as shown here for bone disorders.

Strengths:

Strengths are the clever combination of existing technologies from different fields to build a screening platform. All the required methods are comprehensively described.

Weaknesses:

One may have wished to bring one or two of the crispants to the stage of bona fide mutants, to confirm the results of the screening, however, this is done for some of the tested genes as laid out in the discussion.

Comments on latest version:

All my issues were resolved.

Reviewer #1 (Public review):

Bacterial effectors that interfere with the inner molecular workings of eukaryotic host cells are of great biological significance across disciplines. On the one hand they help us to understand the molecular strategies that bacteria use to manipulate host cells. On the other hand they can be used as research tools to reveal molecular details of the intricate workings of the host machinery that is relevant for the interaction/defence/symbiosis with bacteria. The authors investigate the function and biological impact of a rhizobial effector that interacts with and modifies, and curiously is modified by, legume receptors essential for symbiosis. The molecular analysis revealed a bacterial effector that cleaves a plant symbiosis signaling receptor to inhibit signaling and the host counterplay by phosphorylation via a receptor kinase. These findings have potential implications beyond bacterial interactions with plants.

Bao and colleagues investigated how rhizobial effector proteins can regulate the legume root nodule symbiosis. A rhizobial effector is described to directly modify symbiosis-related signaling proteins, altering the outcome of the symbiosis. Overall, the paper presents findings that will have a wide appeal beyond its primary field.

Out of 15 identified effectors from Sinorhizobium fredii, they focus on the effector NopT, which exhibits proteolytic activity and may therefore cleave specific target proteins of the host plant. They focus on two Nod factor receptors of the legume Lotus japonicus, NFR1 and NFR5, both of which were previously found to be essential for the perception of rhizobial nod factor, and the induction of symbiotic responses such as bacterial infection thread formation in root hairs and root nodule development (Madsen et al., 2003, Nature; Tirichine et al., 2003; Nature). The authors present evidence for an interaction of NopT with NFR1 and NFR5. The paper aims to characterize the biochemical and functional consequences of these interactions and the phenotype that arises when the effector is mutated.

Evidence is presented that in vitro NopT can cleave NFR5 at its juxtamembrane region. NFR5 appears also to be cleaved in vivo. and NFR1 appears to inhibit the proteolytic activity of NopT by phosphorylating NopT. When NFR5 and NFR1 are ectopically over-expressed in leaves of the non-legume Nicotiana benthamiana, they induce cell death (Madsen et al., 2011, Plant Journal). Bao et al., found that this cell death response is inhibited by the coexpression of nopT. Mutation of nopT alters the outcome of rhizobial infection in L. japonicus. These conclusions are well supported by the data.

The authors present evidence supporting the interaction of NopT with NFR1 and NFR5. In particular, there is solid support for cleavage of NFR5 by NopT (Figure 3) and the identification of NopT phosphorylation sites that inhibit its proteolytic activity (Figure 4C). Cleavage of NFR5 upon expression in N. benthamiana (Figure 3A) requires appropriate controls (inactive mutant versions) that have been provided, since Agrobacterium as a closely rhizobia-related bacterium might increase defense related proteolytic activity in the plant host cells.

Key results from N. benthamiana appear consistent with data from recombinant protein expression in bacteria. For the analysis in the host legume L. japonicus transgenic hairy roots were included. To demonstrate that the cleavage of NFR5 occurs during the interaction in plant cells the authors build largely on western blots. Regardless of whether Nicotiana leaf cells or Lotus root cells are used as the test platform, the Western blots indicate that only a small proportion of NFR5 is cleaved when co-expressed with nopT, and most of the NFR5 persists in its full-length form (Figures 3A-D). It is not quite clear how the authors explain the loss of NFR5 function (loss of cell death, impact on symbiosis), as a vast excess of the tested target remains intact. It is also not clear why a large proportion of NFR5 is unaffected by the proteolytic activity of NopT. This is particularly interesting in Nicotiana in the absence of Nod factor that could trigger NFR1 kinase activity.

Comments on latest version:

The presentation of the figures and the language has greatly improved and the specific mistakes pointed out in the last review have been corrected. I especially appreciate the new images used to illustrate the observed mutant phenotypes, which are much clearer and easier to understand. The pictures used to illustrate the mutant phenotypes seem to be of more comparable root regions than before. Overall, the requested changes have been implemented, with some exceptions described below.

• Figure 1: New representative images are shown for BAX1 and CERK1. These pictures are more consistent with the phenotype seen in other treatments, but since the data has not changed, I presume the data from leaf discs (where the leaf discs for these treatments looked very different) previously shown is still included. The criteria for what was considered cell death is in my opinion still not described in the legend. The cell death/total ratio has been added for all leaf discs, as requested.<br /> • Figure 2: the discussion of the figure now emphasizes direct protein interaction. There is still no size marker in 2D or a description of size in the figure legend, making it difficult to compare the result to Figure 3. If I understand the rebuttal comments correctly, there are other bands on the blot, including non-specific bands. This does not negate the need to include the full blot as a supplemental figure to show cleaved NFR5 as well as other bands. I do not see any other clarifications on this subject in the manuscript.<br /> • Figure 5: From the pictures, it is now easier to understand what is meant by "infection foci". Although there is no description in the methods of how these were distinguished from infection threads, I believe the images are clear enough.<br /> • Figure 6: The changes in the discussion are appreciated, but panel E still misrepresents the evidence in the paper, as from the drawing it still seems that the cleaved NFR5 is somehow directly responsible for suppressing infection when this was not shown

Reviewer #1 (Public review):

Summary:

The authors introduced their previous paper with the concise statement that "the relationships between lineage-specific attributes and genotypic differences of tumors are not understood" (Chen et al., JEM 2019, PMID: 30737256). For example, it is not clear why combined loss of RB1 and TP53 is required for tumorigenesis in SCLC or other aggressive neuroendocrine (NE) cancers, or why the oncogenic mutations in KRAS or EGFR that drive NSCLC tumorigenesis are found so infrequently in SCLC. This is the main question addressed by the previous and current papers.

One approach to this question is to identify a discrete set of genetic/biochemical manipulations that are sufficient to transform non-malignant human cells into SCLC-like tumors. One group reported transformation of primary human bronchial epithelial cells into NE tumors through a complex lentiviral cocktail involving inactivation of pRB and p53 and activation of AKT, cMYC and BCL2 (PARCB) (Park et al., Science 2018, PMID: 30287662). The cocktail previously reported by Chen and colleagues to transform human pluripotent stem-cell (hPSC)-derived lung progenitors (LPs) into NE xenografts was more concise: DAPT to inactivate NOTCH signaling combined with shRNAs against RB1 and TP53. However, the resulting RP xenografts lacked important characteristics of SCLC. Unlike SCLC, these tumors proliferated slowly and did not metastasize, and although small subpopulations expressed MYC or MYCL, none expressed NEUROD1.

MYC is frequently amplified or expressed at high levels in SCLC, and here, the authors have tested whether inducible expression of MYC could increase the resemblance of their hPSC-derived NE tumors to SCLC. These RPM cells (or RPM T58A with stabilized cMYC) engrafted more consistently and grew more rapidly than RP cells, and unlike RP cells, formed liver metastases when injected into the renal capsule. Gene expression analyses reveled that RPM tumor subpopulations expressed NEUROD1, ASCL1 and/or YAP1.

The hPSC-derived RPM model is a major advance over the previous RP model. This may become a powerful tool for understanding SCLC tumorigenesis and progression and for discovering gene dependencies and molecular targets for novel therapies. However, the specific role of cMYC in this model needs to be clarified.

Recommended Revision:

cMYC can drive proliferation, tumorigenesis or apoptosis in a variety of lineages depending on concurrent mutations. For example, in the Park et al., study, normal human prostate cells could be reprogrammed to form adenocarcinoma-like tumors by activation of cMYC and AKT alone, without manipulation of TP53 or RB1. In their previous manuscript, the authors carefully showed the role of each molecular manipulation in NE tumorigenesis. DAPT was required for NE differentiation of LPs to PNECs, shRB1 was required for expansion of the PNECs, and shTP53 was required for xenograft formation. cMYC expression could influence each of these steps, and importantly, could render some steps dispensable. For example, shRB1 was previously necessary to expand the DAPT-induced PNECs, as neither shTP53 nor activation of KRAS or EGFR had no effect on this population, but perhaps cMYC overexpression could expand PNECs even in the presence of pRB, or even induce LPs to become PNECs without DAPT. Similarly, both shRB1 and shTP53 were necessary for xenograft formation, but maybe not if cMYC is overexpressed. If a molecular hallmark of SCLC, such as loss of RB1 or TP53, has become dispensable with the addition of cMYC, this information is critically important in interpreting this as a model of SCLC tumorigenesis.

To interpret the role of cMYC expression in hPSC-derived RPM tumors, we need to know what this manipulation does without manipulation of pRB, p53 or NOTCH, alone or in combination. There are 7 relevant combinations that should be presented in this manuscript: (1) cMYC alone in LPs, (2) cMYC + DAPT, (3) cMYC + shRB1, (4) cMYC + DAPT + shRB1, (5) cMYC + shTP53, (6) cMYC + DAPT + shTP53, and (7) cMYC + shRB1 + shTP53. Wild-type cMYC is sufficient; further exploration with the T58A mutant would not be necessary.

Please present the effects of these combinations on LP differentiation to PNECs, expansion of PNECs as well as other lung cells, xenograft formation and histology, and xenograft growth rate and capacity for metastasis. If this could be clarified experimentally, and the results discussed in the context of other similar approaches such as the Park et al., paper, this study would be a major addition to the field.

Reviewer #1 (Public review):

Summary:

Here the authors convincingly identify and characterize the SERBP1 interactome and further define its role in the nucleus, where it is associated with complexes involved in splicing, cell division, chromosome structure, and ribosome biogenesis. Many of the SERBP1-associated proteins are RNA-binding proteins and SERBP1 exerts its impact, at least in part, through these players. SERBP1 is mostly disordered but along with its associated proteins displays a preference for G4 binding and can can bind to PAR and be PARylated. They present data that strongly suggest that complexes in which SERBP1 participates are assembled through G4 or PAR binding. The authors suggest that because SERBP1 lacks traditional functional domains yet is clearly involved in distinct regulatory complexes, SERBP1 likely acts in the early steps of assembly through the recognition of interacting sites present in RNA, DNA, and proteins.

Strengths:

The data is very convincing and demonstrated through multiple approaches.

Weaknesses:

None. The authors have adequately addressed earlier reviewer concerns.

Reviewer #1 (Public review):

Summary:

In cells undergoing Flavivirus infection, cellular translation is impaired but the viruses themselves escape this inhibition and are efficiently translated. In this study, the authors use very elegant and direct approaches to identify the regions in the 5' and 3' UTRs that are important for this phenomenon and then use them to retrieve two cellular proteins that associate with them and mediate translational shutoff evasion (DDX3 and PABP1). A number of experimental approaches are used with a series of well-controlled experiments that fully support the authors' conclusions.

Strengths:

The work identifies the regions in the 5' and 3' UTRs of the viral genome that mediate the escape of JEV from cellular transcriptional shutoff, they evaluate the infectivity of the mutant viruses bearing or not these structures and even explore their pathogenicity in mice. They then identify the cellular proteins that bind to these regions (DDX3 and PABP1) and determine their role in translation blockade escape, in addition to examining and assessing the conservation of the stem-loop identified in JEV in other Flaviviridae.

In almost all of their systematic analyses, translational effects are put in parallel with the replication kinetics of the different mutant viruses. The experimental thread followed in this study is rigorous and direct, and all experiments are truly well-controlled, fully supporting the authors' conclusions

Reviewer #1 (Public review):

Summary:

The authors revealed the cellular heterogeneity of companion cells (CCs) and demonstrated that the florigen gene FT is highly expressed in a specific subpopulation of these CCs in Arabidopsis. Through a thorough characterization of this subpopulation, they further identified NITRATE-INDUCIBLE GARP-TYPE TRANSCRIPTIONAL REPRESSOR 1 (NIGT1)-like transcription factors as potential new regulators of FT. Overall, these findings are intriguing and valuable, contributing significantly to our understanding of florigen and the photoperiodic flowering pathway. However, there is still room for improvement in the quality of the data and the depth of the analysis. I have several comments that may be beneficial for the authors.

Strengths:

The usage of snRNA-seq to characterize the FT-expressing companion cells (CCs) is very interesting and important. Two findings are novel: 1) Expression of FT in CCs is not uniform. Only a subcluster of CCs exhibits high expression level of FT. 2) Based on consensus binding motifs enriched in this subcluster, they further identify NITRATE-INDUCIBLE GARP-TYPE TRANSCRIPTIONAL REPRESSOR 1 (NIGT1)-like transcription factors as potential new regulators of FT.

Weaknesses:

(1) Title: "A florigen-expressing subpopulation of companion cells". It is a bit misleading. The conclusion here is that only a subset of companion cells exhibit high expression of FT, but this does not imply that other companion cells do not express it at all.<br /> (2) Data quality: Authors opted for fluorescence-activated nuclei sorting (FANS) instead of traditional cell sorting method. What is the rationale behind this decision? Readers may wonder, especially given that RNA abundance in single nuclei is generally lower than that in single cells. This concern also applies to snRNA-seq data. Specifically, the number of genes captured was quite low, with a median of only 149 genes per nucleus. Additionally, the total number of nuclei analyzed was limited (1,173 for the pFT:NTF and 3,650 for the pSUC2:NTF). These factors suggest that the quality of the snRNA-seq data presented in this study is quite low. In this context, it becomes challenging for the reviewer to accurately assess whether this will impact the subsequent conclusions of the paper. Would it be possible to repeat this experiment and get more nuclei?<br /> (3) Another disappointment is that the authors did not utilize reporter genes to identify the specific locations of the FT-high expressing cells (cluster 7 cells) within the CC population in vivo. Are there any discernible patterns that can be observed?<br /> (4) The final disappointment is that the authors only compared FT expression between the nigtQ mutants and the wild type. Does this imply that the mutant does not have a flowering time defect particularly under high nitrogen conditions?

Reviewer #1 (Public review):

Summary

Chabukswar et al analysed endogenous retrovirus (ERV) Env variation in a set of primate genomes using consensus Env sequences from ERVs known to be present in hominoids using a Blast homology search with the aim of characterising env gene changes over time. The retrieved sequences were analysed phylogenetically, and showed that some of the integrations are LTR-env recombinants.

Strengths

The strength of the manuscript is that such an analysis has not been performed yet for the subset of ERV Env genes selected and most of the publicly available primate genomes.

Weaknesses

Unfortunately, the weaknesses of the manuscript outnumber its strengths. Especially the methods section does not contain sufficient information to appreciate or interpret the results. The results section contains methodological information that should be moved, while the presentation of the data is often substandard. For instance, the long lists of genomes in which a certain Env was found could better be shown in tables. Furthermore, there is no overview of the primate genomes, or accession numbers, used. It is unclear whether the analyses, such as the phylogenetic trees, are based on nucleotide or amino acid sequences since this is not stated. tBLASTn was used in the homology searches, so one would suppose aa are retrieved. In the Discussion, both env (nt?) and Env (aa?) are used.

For the non-hominoids, genome assembly of publicly available sequences is not always optimal, and this may require Blasting a second genome from a species. Which should for instance be done for the HML2 sequences found in the Saimiri boliviensis genome, but not in the related Callithrix jacchus genome. Finally, the authors propose to analyse recombination in Env sequences but only retrieve env-LTR recombinant Envs, which should likely not have passed the quality check.

Since the Methods section does not contain sufficient information to understand or reproduce the results, while the Results are described in a messy way, it is unclear whether or not the aims have been achieved. I believe not, as characterisation of env gene changes over time is only shown for a few abberrant integrations containing part of the LTR in the env ORF.

Reviewer #1 (Public review):

Summary:

In this study, the authors sequenced emm89 serotype genomes of clinical isolates from patients in Japan, where the number of invasive Group A Streptococcus (GAS), especially those of the emm89 serotype, has drastically increased over the past 10-15 years. The sequences from this cohort were compared against a large collection of publicly available global isolates, yielding a total of almost 1000 genomes in the analysis. Because the researchers focused on the emm89 serotype, they could construct a common core genome, with subsequent ability to analyze genomic differences in accessory genes and intergenic regions that contributed to the invasive phenotype using multiple types of GWAS analysis (SNP, k-mer). Their analysis demonstrates some mutations responsible for invasiveness are specific to the Japanese strains, and that multiple independent virulence factors can contribute to invasiveness. None of the invasive phenotypes were correlated with new gene acquisition. Together, the data support that synergy between bacterial survival and upregulation of virulence factors contributes to the development of severe infection.

Strengths:

• The authors verify their analysis by confirming that covS is one of the more frequently mutated genes in invasive strains of GAS, as has been shown in other publications.

• A mutation in one of the SNPs attributed to invasiveness (SNP fhuB) was introduced into an invasive strain. The authors demonstrate that this mutant strain survives less well in human blood. Therefore, the authors have experimental data to support their claims that their analysis uncovered a new mutation/SNP that contributed to invasiveness.

Weaknesses:

• It would be helpful for the authors to highlight why their technique (large scale analysis of one emm type) can yield more information than a typical GWAS analysis of invasive vs. non-invasive strains. Are SNPs easier to identify using a large-scale core genome? Is it more likely evolutionarily to find mutations in non-coding regions as opposed to the core genome and accessory genes, and this is what this technique allows? Did the analysis yield unexpected genes or new genes that had not been previously identified in other GWAS analyses? These points may need to be made more apparent in the results and deserve some thought in the discussion section.

• The Alpha-fold data does not demonstrate why the mutations the authors identified could contribute to the invasive phenotype. It would be helpful to show an overlay of the predicted structures containing the different SNPs to demonstrate the potential structural differences that can occur due to the SNP. This would make the data more convincing that the SNP has a potential impact on the function of the protein. Similarly, the authors discuss modification of the hydrophobicity of the side chain in the ferrichrome transporter (lines 317-318) due to a SNP, but this is not immediately obvious in the figure (Fig. 5).

Reviewer #1 (Public review):

Overall, the manuscript reveals the role of actin polymerization to drive the fusion of myoblasts during adult muscle regeneration. This pathway regulates fusion in many contexts, but whether it was conserved in adult muscle regeneration remained unknown. Robust genetic tools and histological analyses were used to support the claims convincingly.

There are a few interpretations that could be adjusted.

The beginning of the results about macrophages traversing ghost fibers after regeneration was a surprise given the context in the abstract and introduction. These results also lead to new questions about this biology that would need to be answered to substantiate the claims in this section. Also, it is unclear the precise new information learned here because it seems obvious that macrophages would need to extravasate the basement membrane to enter ghost fibers and macrophages are known to have this ability. Moreover, the model in Figure 4D has macrophages and BM but there is not even mention of this in the legend. The authors may wish to consider removing this topic from the manuscript.

Which Pax7CreER line was used? In the methods, the Jax number provided is the Gaka line but in the results, Lepper et al 2009 are cited, which is not the citation for the Gaka line.

Did the authors assess regeneration in the floxed mice that do not contain Cre as a control? Or is it known these alleles do not perturb the function of the targeted gene?

The authors comment: 'Interestingly, expression of the fusogenic proteins, MymK and MymX, was up-regulated in the TA muscle of these mice (Fig. S4F), suggesting that fusogen overexpression is not able to rescue the SCM fusion defect resulted from defective branched actin polymerization.' It is unclear if fusogens are truly overexpressed because the analysis is performed at dpi 4 when the expression of fusogens may be decreased in control mice because they have already fused. Also, only two animals were analyzed and it is unclear if MymX is definitively increased. The authors should consider adjusting the interpretation to SCM fusion defect resulting from defective branched actin polymerization is unlikely to be caused by a lack of fusogen expression.

Reviewer #1 (Public review):

Summary:

The authors use FIB SEM methods to generate 3D volumes of almost all cells comprising a miniature wasp eye and describe the anatomy of each cell type in detail. The function of each cell type is determined through comparisons with descriptions using other methods from larger insect species.

Strengths:

The data show that, despite the small size, many elements of the eye are consistent with those found in larger insects. In addition, the powerful FIB-SEM technique revealed a hitherto unknown case of ectopic photoreceptors.

Weaknesses:

As this paper only uses anatomical analyses, no functional interpretations of cell function are tested.

The aim of this paper was to describe the ultrastructural organization of compound eyes in the extremely small wasp Megaphragma viggianii. The authors successfully achieved this aim and provided an incredibly detailed description of all cell types with respect to their location, volume, and dimensions. As this is the first of its kind, the results cannot easily be compared with previous work. The findings are likely to be an important reference for future work that uses similar techniques to reconstruct the eyes of other insect species. The FIB-SEM method used is being used increasingly often in structural studies of insect sensory organs and brains and this work demonstrates the utility of this method.

Reviewer #1 (Public review):

Summary:

Mallimadugula et al. combined Molecular Dynamics (MD) simulations, thiol-labeling experiments, and RNA-binding assays to study and compare the RNA-binding behavior of the Interferon Inhibitory Domain (IID) from Viral Protein 35 (VP35) of Zaire ebolavirus, Reston ebolavirus, and Marburg marburgvirus. Although the structures and sequences of these viruses are similar, the authors suggest that differences in RNA binding stem from variations in their intrinsic dynamics, particularly the opening of a cryptic pocket. More precisely, the dynamics of this pocket may influence whether the IID binds to RNA blunt ends or the RNA backbone.

Overall, the authors present important findings to reveal how the intrinsic dynamics of proteins can influence their binding to molecules and, hence, their functions. They have used extensive biased simulations to characterize the opening of a pocket which was not clearly seen in experimental results - at least when the proteins were in their unbound forms. Biochemical assays further validated theoretical results and linked them to RNA binding modes. Thus, with the combination of biochemical assays and state-of-the-art Molecular Dynamics simulations, these results are clearly compelling.

Strengths:

The use of extensive Adaptive Sampling combined with biochemical assays clearly points to the opening of the Interferon Inhibitory Domain (IID) as a factor for RNA binding. This type of approach is especially useful to assess how protein dynamics can affect its function.

Weaknesses :

Although a connection between the cryptic pocket dynamics and RNA binding mode is proposed, the precise molecular mechanism linking pocket opening to RNA binding still remains unclear.

Reviewer #1 (Public review):

Summary:

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

Strengths:

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

Weaknesses:

Below are specific comments and concerns:

(1) The diversity of cell types recovered from scRNA-seq libraries of sorted Nes-CFP cells is unexpected, especially the inclusion of minor types such as microglia, meninges, and ependymal cells. The authors should validate whether Nes and CFP mRNAs are enriched in the sorted cells; if not, they should discuss the potential pitfalls in sampling bias or artifacts that may have affected the dataset, impacting interpretation.<br /> (2) The authors should de-emphasize that ROS signaling and related gene upregulation exclusively in gliogenic NEPs. Genes such as Cdkn1a, Phlda3, Ass1, and Bax are identified as differentially expressed in neurogenic NEPs and granule cell progenitors (GCPs), with Ass1 absent in GCPs. According to Table S4, gene ontology (GO) terms related to ROS metabolic processes are also enriched in gliogenic NEPs, neurogenic NEPs, and GCPs.<br /> (3) The authors need to justify the selection of only the anterior lobe for EGL replenishment and microglia quantification.<br /> (4) Figure 1K: The figure presents linkages between genes and GO terms as a network but does not depict a gene network. The terminology should be corrected accordingly.<br /> (5) Figure 1H and S2: The x-axis appears to display raw p-values rather than log10(p.value) as indicated. The x-axis should ideally show -log10(p.adjust), beginning at zero. The current format may misleadingly suggest that the ROS GO term has the lowest p-values.<br /> (6) Genes such as Ppara, Egln3, Foxo3, Jun, and Nos1ap were identified by bulk ATAC-seq based on proximity to peaks, not by scRNA-seq. Without additional expression data, caution is needed when presenting these genes as direct evidence of ROS involvement in NEPs.<br /> (7) The authors should annotate cell identities for the different clusters in Table S2.<br /> (8) Reiterative clustering analysis reveals distinct subpopulations among gliogenic and neurogenic NEPs. Could the authors clarify the identities of these subclusters? Can we distinguish the gliogenic NEPs in the Bergmann glia layer from those in the white matter?<br /> (9) In the Methods section, the authors mention filtering out genes with fewer than 10 counts. They should specify if these genes were used as background for enrichment analysis. Background gene selection is critical, as it influences the functional enrichment of gene sets in the list.<br /> (10) Figure S1C: The authors could consider using bar plots to better illustrate cell composition differences across conditions and replicates.<br /> (11) Figures 4-6: It remains unclear how the white matter microglia contribute to the recruitment of BgL-NEPs to the EGL, as the mCAT-mediated microglia loss data are all confined to the white matter.

Reviewer #1 (Public review):

Summary:

This manuscript investigates genes that escape X-Chromosome Inactivation (XCI) across human tissues, using females that exhibit skewed or non-random XCI. The authors identified 2 female individuals with skewed XCI in the GTex database, in addition to the 1 female skewed sample in this database that has been described in a previous publication (Ref.16). The authors also determined the genes that escape XCI for 380 X-linked genes across 30 different tissues.

Strengths:

The novelty of this manuscript is that the authors have identified the XCI expression status for a total of 380 genes across 30 different human tissues, and also discovered the XCI status (escape, variable escape, or silenced) for 198 X-linked genes, whose status was previously not determined. This report is a good resource for the field of XCI, and would benefit from additional analyses and clarification of their comparisons of XCI status.

Weaknesses:

Specific comments:

(1) The authors state that they have reclassified the allelic expression status of 32 genes (shown in Table S5, Supplementary Figure 3). The concern is the source of the tissue or cell line which was originally used to make the classification of XCI status, and whether the comparisons are equivalent. For example, if cell lines (and not tissues) were used to define the XCI status for EGFL6, TSPAN6, and CXorf38, then how can the authors be sure that the escape status in whole tissues would be the same? Also along these lines, the authors should consider whether escape status in previous studies using immortalized/cancer cell lines (such as the meta analyses done in Balaton publication) would be different compared to healthy tissues (seems like it should be). Therefore making comparisons between healthy whole tissues and cancer cell lines doesn't make sense.

(2) The authors note that skewed XCI is prevalent in the human population, and cite some publications (references 8, 10-12). If RNAseq data is available from these female individuals with skewed XCI (such as ref 12), the authors should consider using their allelic expression pipeline to identify XCI status of more X-linked genes.

(3) It has been well established that the human inactive X has more XCI escape genes compared to the mouse inactive X. In light of the author's observations across human tissues, how does the XCI status compare with the same tissues in mice?

Reviewer #1 (Public review):

Summary<br /> In this beautiful paper the authors examined the role and function of NR2F2 in testis development and more specifically on fetal Leydig cells development. It is well known by now that FLC are developed from an interstitial steroidogenic progenitors at around E12.5 and are crucial for testosterone and INSL3 production during embryonic development, which in turn shapes the internal and external genitalia of the male. Indeed, lack of testosterone or INSL3 are known to cause DSD as well as undescended testis, also termed as cryptorchidism.<br /> The authors first characterized the expression pattern of the NR2R2 protein during testis development and then used two cKO systems of NR2F2, namely the Wt1-creERT2 and the Nr5a1-cre to explore the phenotype of loss of NR2F2. They found in both cases that mice are presenting with undescended testis and major reduction in FLC numbers. They show that NR2F2 has no effect on the amount and expression of the progenitor cells but in its absence, there are less FLC and they are immature.<br /> The effect of NR2F2 is cell autonomous and does not seem to affect other signalling pathways implemented in Leydig cell development as the DHH, PDGFRA and the NOTCH pathway.

Overall, this paper is excellent, very well written, fluent and clear. The data is well presented, and all the controls and statistics are in place. I think this paper will be of great interest to the field and paves the way for several interesting follow up studies as stated in the discussion

Reviewer #1 (Public review):

Summary:

This study examines the role of host blood meal source, temperature, and photoperiod on the reproductive traits of Cx. quinquefasciatus, an important vector of numerous pathogens of medical importance. The host use pattern of Cx. quinquefasciatus is interesting in that it feeds on birds during spring and shifts to feeding on mammals towards fall. Various hypotheses have been proposed to explain the seasonal shift in host use in this species but have provided limited evidence. This study examines whether the shifting of host classes from birds to mammals towards autumn offers any reproductive advantages to Cx. quinquefasciatus in terms of enhanced fecundity, fertility, and hatchability of the offspring. The authors found no evidence of this, suggesting that alternate mechanisms may drive the seasonal shift in host use in Cx. quinquefasciatus.

Strengths:

Host blood meal source, temperature, and photoperiod were all examined together.

Weaknesses:

The study was conducted in laboratory conditions with a local population of Cx. quinquefasciatus from Argentina. I'm not sure if there is any evidence for a seasonal shift in the host use pattern in Cx. quinquefasciatus populations from the southern latitudes.

Comments on the revision:

Overall, the manuscript is much improved. However, the introduction and parts of the discussion that talk about addressing the question of seasonal shift in host use pattern of Cx. quin are still way too strong and must be toned down. There is no strong evidence to show this host shift in Argentinian mosquito populations. Therefore, it is just misleading. I suggest removing all this and sticking to discussing only the effects of blood meal source and seasonality on the reproductive outcomes of Cx. quin.

Reviewer #1 (Public review):

The aim of this paper is to describe a novel method for genetic labelling of animals or cell populations, using a system of DNA/RNA barcodes.

Strengths:

• The author's attempt at providing a straightforward method for multiplexing Drosophila samples prior to scRNA-seq is commendable. The perspective of being able to load multiple samples on a 10X Chromium without antibody labelling is appealing.<br /> • The authors are generally honest about potential issues in their method, and areas that would benefit from future improvement.<br /> • The article reads well. Graphs and figures are clear and easy to understand.

Weaknesses:

• The usefulness of TaG-EM for phototaxis, egg laying or fecundity experiments is questionable. The behaviours presented here are all easily quantifiable, either manually or using automated image-based quantification, even when they include a relatively large number of groups and replicates. Despite their claims (e.g., L311-313), the authors do not present any real evidence about the cost- or time-effectiveness of their method in comparison to existing quantification methods.<br /> • Behavioural assays presented in this article have clear outcomes, with large effect sizes, and therefore do not really challenge the efficiency of TaG-EM. By showing a T-maze in Fig 1B, the authors suggest that their method could be used to quantify more complex behaviours. Not exploring this possibility in this manuscript seems like a missed opportunity.<br /> • Experiments in Figs S3 and S6 suggest that some tags have a detrimental effect on certain behaviours or on GFP expression. Whereas the authors rightly acknowledge these issues, they do not investigate their causes. Unfortunately, this question the overall suitability of TaG-EM, as other barcodes may also affect certain aspects of the animal's physiology or behaviour. Revising barcode design will be crucial to make sure that sequences with potential regulatory function are excluded.<br /> • For their single-cell experiments, the authors have used the 10X Genomics method, which relies on sequencing just a short segment of each transcript (usually 50-250bp - unknown for this study as read length information was not provided) to enable its identification, with the matching paired-end read providing cell barcode and UMI information (Macosko et al., 2015). With average fragment length after tagmentation usually ranging from 300-700bp, a large number of GFP reads will likely not include the 14bp TaG-EM barcode. When a given cell barcode is not associated with any TaG-EM barcode, then demultiplexing is impossible. This is a major problem, which is particularly visible in Figs 5 and S13. In 5F, BC4 is only detected in a couple of dozen cells, even though the Jon99Ciii marker of enterocytes is present in a much larger population (Fig 5C). Therefore, in this particular case, TaG-EM fails to detect most of the GFP-expressing cells. Similarly, in S13, most cells should express one of the four barcodes, however many of them (maybe up to half - this should be quantified) do not. Therefore, the claim (L277-278) that "the pan-midgut driver were broadly distributed across the cell clusters" is misleading. Moreover, the hypothesis that "low expressing driver lines may result in particularly sparse labelling" (L331-333) is at least partially wrong, as Fig S13 shows that the same Gal4 driver can lead to very different levels of barcode coverage.<br /> • Comparisons between TaG-EM and other, simpler methods for labelling individual cell populations are missing. For example, how would TaG-EM compare with expression of different fluorescent reporters, or a strategy based on the brainbow/flybow principle?<br /> • FACS data is missing throughout the paper. The authors should include data from their comparative flow cytometry experiment of TaG-EM cells with or without additional hexameric GFP, as well as FSC/SSC and fluorescence scatter plots for the FACS steps that they performed prior to scRNA-seq, at least in supplementary figures.<br /> • The authors should show the whole data described in L229, including the cluster that they chose to delete. At least, they should provide more information about how many cells were removed. In any case, the fact that their data still contains a large number of debris and dead cells despite sorting out PI negative cells with FACS and filtering low abundance barcodes with Cellranger is concerning.

Overall, although a method for genetic tagging cell populations prior to multiplexing in single-cell experiments would be extremely useful, the method presented here is inadequate. However, despite all the weaknesses listed above, the idea of barcodes expressed specifically in cells of interest deserves more consideration. If the authors manage to improve their design to resolve the major issues and demonstrate the benefits of their method more clearly, then TaG-EM could become an interesting option for certain applications.

Comments on revisions:

The authors have addressed many important points, providing reassurances about the initial weaknesses of their work. Although the TaG-EM is unlikely to have a significant influence on the field due to its limited benefits, the results are now sound and provide the reader with an unbiased view of the possibilities and limitations of the method.

Reviewer #1 (Public review):

Summary:

BMP signaling is, arguably, best known for its role in the dorsoventral patterning, but not in nematodes, where it regulates body size. In their paper, Vora et al. analyze ChIP-Seq and RNA-Seq data to identify direct transcriptional targets of SMA-3 (Smad) and SMA-9 (Schnurri) and understand the respective roles of SMA-3 and SMA-9 in the nematode model Caenorhabditis elegans. The Authors use SMA-3 and SMA-9 ChIP-Seq data and RNA-Seq data from SMA-3 and SMA-9 mutants, and bioinformatic analyses to identify the genes directly controlled by these two transcription factors (TFs) and find approximately 350 such targets for each. They show that all SMA-3-controlled targets are positively controlled by SMA-3 binding, while SMA-9-controlled targets can be either up- or downregulated by SMA-9. 129 direct targets were shared by SMA-3 and SMA-9, and, curiously, the expression of 15 of them was activated by SMA-3 but repressed by SMA-9. In case of such opposing effects, the SMA-9 appears to act epistatically to SMA-3. Since genes responsible for cuticle collagen production were eminent among the SMA-3 targets, the Authors focused on trying to understand the body size defect known to be elicited by the modulation of BMP signaling. Vora et al. provide compelling evidence that this defect is likely to be due to problems with the BMP signaling-dependent collagen secretion necessary for cuticle formation.

Strengths:

Vora et al. provide a valuable analysis of ChIP-Seq and RNA-Seq datasets, which will be very useful for the community. They also shed light on the mechanism of the BMP-dependent body size control by identifying SMA-3 target genes regulating cuticle collagen synthesis and by showing that downregulation of these genes affects body size in C. elegans.

Weaknesses:

(1) Although the analysis of the SMA-3 and SMA-9 ChIP-Seq and RNA-Seq data is extremely useful, the goal "to untangle the roles of Smad and Schnurri transcription factors in the developing C. elegans larva", has not been reached. While the role of SMA-3 as a transcriptional activator appears to be quite straightforward, the function of SMA-9 in the BMP signaling remains obscure.

(2) The Authors clearly show that both TFs can bind independently of each other, however, by using distances between SMA-3 and SMA-9 ChIP peaks, they claim that when the peaks are close these two TFs likely act as complexes. In the absence of proof that SMA-3 and SMA-9 physically interact (e.g. that they co-immunoprecipitate - as they do in Drosophila), this is an unfounded claim, which still has to be experimentally substantiated. In the revised version of the manuscript, the authors acknowledge this.

(3) The second part of the results (the collagen story) is loosely connected the first part. dpy-11 encodes an enzyme important for cuticle development, and it is a differentially expressed direct target of SMA-3. dpy-11 can be bound by SMA-9, but it is not affected by this binding according to RNA-Seq. Thus, technically, this part of the paper does not require any information about SMA-9. However, this can likely be improved by addressing the function of the 15 genes, with the opposing mode of regulation by SMA-3 and SMA-9.

Comments on revisions:

In comparison to the first version of the manuscript, the authors have significantly improved the "readability" of the paper, made the Discussion much better, and toned down some of the less supported arguments.

Reviewer #1 (Public review):

Summary:

The study aims to create a comprehensive repository about the changes in protein abundance and their modification during oocyte maturation in Xenopus laevis.

Strengths:

The results contribute meaningfully to the field.

Weaknesses:

The manuscript could have benefitted from more comprehensive analyses and clearer writing. Nonetheless, the key findings are robust and offer a valuable resource for the scientific community.

Reviewer #1 (Public review):

Summary:

In this article, Gupta and colleagues explore the parameters that could promote the elimination of active Ras cells when surrounded by WT cells. The elimination of active Ras cells by surrounding WT cells was previously described extensively and associated with a process named cell competition, a context dependant elimination of cells. Several mechanisms have been associated with competition, including more recently elimination processes based on mechanical stress. This was explored theoretically and experimentally and was either associated with differential growth and sensitivity to pressure and/or differences in homeostatic density/pressure. This was extensively validated for the case of Scribble mutant cells which are eliminated by WT MDCK cells due to their higher homeostatic density. However, there has been so far very little systematic characterisation of the mechanical parameters and properties of these different cell types and how this could contribute to mechanical competition.

Here, the authors used the context of active Ras cells in MDCK cells (with some observations in vivo in mice gut which are a bit more anecdotal) to explore the parameters causal to Ras cell elimination. Using for the first time traction force microscopy, stress microscopy combined with Bayesian inference, they first show that clusters of active Ras cells experience higher pressure compared to WT. Interestingly, this occurs in absence of differences in growth rate, and while Ras cells seems to have lower homeostatic density, in contractions with the previous models associated with mechanical cell competition. Using a self-propelled Voronoi model, they explored more systematically the conditions that will promote the compression of transformed cells, showing globally that higher Area compressibility and/or lower junctional tension are associated with higher compressibility. Using then an original and novel experimental method to measure bulk compressibility of cell populations, they confirmed that active Ras cells are globally twice more compressible than WT cells. This compressibility correlates with a disruption of adherens junctions. Accordingly, the higher pressure near transformed Ras cells can be completely rescued by increasing cell-cell adhesion through E-cad overexpression, which also reduces the compressibility of the transformed cells. Altogether, these results go along the lines of a previous theoretical work (Gradeci et al. eLife 2021) which was suggesting that reduced stiffness/higher compressibility was essential to promote loser cell elimination. Here, the authors provide for the first time a very convincing experimental measurement and validation of this prediction. Moreover, their modelling approach goes far beyond what was performed before in terms of exploration of conditions promoting compressibility, and their experimental data point at alternative mechanisms that may contribute to mechanical competition.

Strengths:

- Original methodologies to perform systematic characterisation of mechanical properties of Ras cells during cell competition, which include a novel method to measure bulk compressibility.<br /> - A very extensive theoretical exploration of the parameters promoting cell compaction in the context of competition.

Weaknesses:

- Most of the theoretical focus is centred on the bulk compressibility, but so far does not really explain the final fate of the transformed cells. Classic cell competition scenario (including the one involving active Ras cells) lead to the elimination of one cell population either by cell extrusion/cell death or global delamination. This aspect is absolutely not explored in this article, experimentally or theoretically, and as such it is difficult to connect all the observables with the final outcome of cell competition. For instance, higher compressibility may not lead to loser status if the cells can withstand high density without extruding compared to the WT cells (and could even completely invert the final outcome of the competition). Down the line, and as suggested in most of the previous models/experiments, the relationship between pressure/density and extrusion/death will be the key factor that determine the final outcome of competition. However, there is absolutely no characterisation of cell death/cell extrusion in the article so far.

- While the compressibility measurement are very original and interesting, this bulk measurement could be explained by very different cellular processes, from modulation of cell shape, to cell extrusion and tissue multilayering (which by the way was already observed for active Ras cells, see for instance https://pubmed.ncbi.nlm.nih.gov/34644109/). This could change a lot the interpretation of this measurement and to which extend it can explain the compression observed in mixed culture. This compressibility measurement could be much more informative if coupled with an estimation of the change of cell aspect ratio and the rough evaluation of the contribution of cell shape changes versus alternative mechanisms.

- So far, there is no clear explanation of why transformed Ras cells get more compacted in the context of mixed culture compared to pure Ras culture. Previously, the compaction of mutant Scribble cells could be explained by the higher homeostatic density of WT cells which impose their prefered higher density to Scribble mutant (see Wagstaff et al. 2016 or Gradeci et al 2021), however that is not the case of the Ras cells (which have even slightly higher density at confluency). If I understood properly, the Voronoid model assumes some directional movement of WT cell toward transformed which will actively compact the Ras cells through self-propelled forces (see supplementary methods), but this is never clearly discussed/described in the results section, while potentially being one essential ingredient for observing compaction of transformed cells. In fact, this was already described experimentally in the case of Scribble competition and associated with chemoattractant secretion from the mutant cells promoting directed migration of the WT (https://pubmed.ncbi.nlm.nih.gov/33357449/). It would be essential to show what happens in absence of directional propelled movement in the model and validate experimentally whether there is indeed directional movement of the WT toward the transformed cells. Without this, the current data does not really explain the competition process.

- Some of the data lack a bit of information on statistic, especially for all the stress microscopy and traction forces where we do no really know how representative at the stress patterns (how many experiment, are they average of several movies ? integrated on which temporal window ?)

Reviewer #1 (Public review):

Summary:

In this article, Gupta and colleagues explore the parameters that could promote the elimination of active Ras cells when surrounded by WT cells. The elimination of active Ras cells by surrounding WT cells was previously described extensively and associated with a process named cell competition, a context dependant elimination of cells. Several mechanisms have been associated with competition, including more recently elimination processes based on mechanical stress. This was explored theoretically and experimentally and was either associated with differential growth and sensitivity to pressure and/or differences in homeostatic density/pressure. This was extensively validated for the case of Scribble mutant cells which are eliminated by WT MDCK cells due to their higher homeostatic density. However, there has been so far very little systematic characterisation of the mechanical parameters and properties of these different cell types and how this could contribute to mechanical competition.

Here, the authors used the context of active Ras cells in MDCK cells (with some observations in vivo in mice gut which are a bit more anecdotal) to explore the parameters causal to Ras cell elimination. Using for the first time traction force microscopy, stress microscopy combined with Bayesian inference, they first show that clusters of active Ras cells experience higher pressure compared to WT. Interestingly, this occurs in absence of differences in growth rate, and while Ras cells seems to have lower homeostatic density, in contractions with the previous models associated with mechanical cell competition. Using a self-propelled Voronoi model, they explored more systematically the conditions that will promote the compression of transformed cells, showing globally that higher Area compressibility and/or lower junctional tension are associated with higher compressibility. Using then an original and novel experimental method to measure bulk compressibility of cell populations, they confirmed that active Ras cells are globally twice more compressible than WT cells. This compressibility correlates with a disruption of adherens junctions. Accordingly, the higher pressure near transformed Ras cells can be completely rescued by increasing cell-cell adhesion through E-cad overexpression, which also reduces the compressibility of the transformed cells. Altogether, these results go along the lines of a previous theoretical work (Gradeci et al. eLife 2021) which was suggesting that reduced stiffness/higher compressibility was essential to promote loser cell elimination. Here, the authors provide for the first time a very convincing experimental measurement and validation of this prediction. Moreover, their modelling approach goes far beyond what was performed before in terms of exploration of conditions promoting compressibility, and their experimental data point at alternative mechanisms that may contribute to mechanical competition.

Strengths:

- Original methodologies to perform systematic characterisation of mechanical properties of Ras cells during cell competition, which include a novel method to measure bulk compressibility.<br /> - A very extensive theoretical exploration of the parameters promoting cell compaction in the context of competition.

Weaknesses:

- Most of the theoretical focus is centred on the bulk compressibility, but so far does not really explain the final fate of the transformed cells. Classic cell competition scenario (including the one involving active Ras cells) lead to the elimination of one cell population either by cell extrusion/cell death or global delamination. This aspect is absolutely not explored in this article, experimentally or theoretically, and as such it is difficult to connect all the observables with the final outcome of cell competition. For instance, higher compressibility may not lead to loser status if the cells can withstand high density without extruding compared to the WT cells (and could even completely invert the final outcome of the competition). Down the line, and as suggested in most of the previous models/experiments, the relationship between pressure/density and extrusion/death will be the key factor that determine the final outcome of competition. However, there is absolutely no characterisation of cell death/cell extrusion in the article so far.

- While the compressibility measurement are very original and interesting, this bulk measurement could be explained by very different cellular processes, from modulation of cell shape, to cell extrusion and tissue multilayering (which by the way was already observed for active Ras cells, see for instance https://pubmed.ncbi.nlm.nih.gov/34644109/). This could change a lot the interpretation of this measurement and to which extend it can explain the compression observed in mixed culture. This compressibility measurement could be much more informative if coupled with an estimation of the change of cell aspect ratio and the rough evaluation of the contribution of cell shape changes versus alternative mechanisms.

- So far, there is no clear explanation of why transformed Ras cells get more compacted in the context of mixed culture compared to pure Ras culture. Previously, the compaction of mutant Scribble cells could be explained by the higher homeostatic density of WT cells which impose their prefered higher density to Scribble mutant (see Wagstaff et al. 2016 or Gradeci et al 2021), however that is not the case of the Ras cells (which have even slightly higher density at confluency). If I understood properly, the Voronoid model assumes some directional movement of WT cell toward transformed which will actively compact the Ras cells through self-propelled forces (see supplementary methods), but this is never clearly discussed/described in the results section, while potentially being one essential ingredient for observing compaction of transformed cells. In fact, this was already described experimentally in the case of Scribble competition and associated with chemoattractant secretion from the mutant cells promoting directed migration of the WT (https://pubmed.ncbi.nlm.nih.gov/33357449/). It would be essential to show what happens in absence of directional propelled movement in the model and validate experimentally whether there is indeed directional movement of the WT toward the transformed cells. Without this, the current data does not really explain the competition process.

- Some of the data lack a bit of information on statistic, especially for all the stress microscopy and traction forces where we do no really know how representative at the stress patterns (how many experiment, are they average of several movies ? integrated on which temporal window ?)

Reviewer #1 (Public review):

Summary:

Olfaction is fundamental to the survival and reproduction of animals, as they rely on olfactory sensory neurons (OSNs) in the olfactory epithelium (OE) to detect volatile chemical cues in their environment. Most mature OSNs adhere to the 'one neuron one receptor' rule, wherein each neuron selects a single receptor for expression from a large repertoire of olfactory receptor genes. The precise regulation of olfactory receptor expression is critical for accurate odorant recognition. Since the seminal discovery of olfactory receptors by Linda Buck and Richard Axel in 1991, substantial efforts have been made to elucidate the mechanisms underlying OSN differentiation and receptor expression. However, these processes remain incompletely understood. The development of in vitro olfactory epithelium organoids offers a promising platform to address these fundamental questions. The in vivo OE is composed of a complex array of cell types, which has posed a significant challenge for recapitulating its structure and function in vitro. Previous attempts to generate olfactory organoids from adult human or mouse OE cells yielded tissue containing OSNs, but these constructs were structurally distinct from the in vivo OE and lacked the characteristic pseudostratified epithelium.

In this study, Kazuya et al. successfully established olfactory epithelium organoids from E13.5 mouse embryonic OE stem cells, which developed into a pseudostratified structure closely resembling the native OE. They further examined the influence of different cultural conditions on OE differentiation, confirming the pivotal role of niche factors in promoting OSN development. Through immunofluorescence staining and single-cell RNA sequencing, they demonstrated that the organoids encompass a diverse range of cell types analogous to those present in the in vivo OE. Notably, calcium imaging revealed that the organoids were functionally responsive to odorants, and single-cell transcriptomic analysis showed that the majority of mature OSNs conformed to the 'one neuron one receptor' rule. Using these organoids, the authors performed a preliminary investigation into the developmental trajectories of OSNs, developed a tool to predict subpopulations of mature OSNs, and identified novel markers associated with OSN maturation. Collectively, the data provide compelling evidence for the reliability and utility of this olfactory organoid model. Further in-depth analyses may enable readers to better assess and utilize this tool to advance the study of olfactory biology.

Strengths:

The authors developed and established olfactory epithelium organoids, with immunofluorescence imaging confirming the presence of a pseudostratified structure similar to that of the in vivo olfactory epithelium, representing a significant advancement. Single-cell sequencing and calcium imaging further demonstrated the utility of these organoids, as they contain multiple cell types analogous to the in vivo olfactory epithelium. Importantly, they are physiologically functional, capable of responding to odor stimuli.

Weakness:

Although the authors have made significant progress in the technique, there are some gaps in understanding its underlying principles. First, it remains unclear what specific characteristics of E13.5 embryonic olfactory stem cells enable them to generate organoids in vitro that more closely resemble the in vivo olfactory epithelium, compared to adult mouse olfactory stem cells. Second, it is not clearly defined which specific cell type(s) from the embryonic olfactory epithelium give rise to these organoids, and the efficiency of organoid formation from the isolated cells also warrants further clarification.

Reviewer #1 (Public review):

Summary:

van der Linden et al. report on the development of a new green-fluorescent sensor for calcium, following a novel rational design strategy based on the modification of the cyan-emissive sensor mTq2-CaFLITS. Through a mutational strategy similar to the one used to convert EGFP into EYFP, coupled with optimization of strategic amino acids located in proximity of the chromophore, they identify a novel sensor, G-CaFLITS. Through a careful characterization of the photophysical properties in vitro and the expression level in cell cultures, the authors demonstrate that G-CaFLITS combines a large lifetime response with a good brightness in both the bound and unbound states. This relative independence of the brightness on calcium binding, compared with existing sensors that often feature at least one very dim form, is an interesting feature of this new type of sensors, which allows for a more robust usage in fluorescence lifetime imaging. Furthermore, the authors evaluate the performance of G-CaFLITS in different subcellular compartments and under two-photon excitation in Drosophila. While the data appears robust and the characterization thorough, the interpretation of the results in some cases appears less solid, and alternative explanations cannot be excluded.

Strengths:

- The approach is innovative and extends the excellent photophysical properties of the mTq2-based to more red-shifted variants. While the spectral shift might appear relatively minor, as the authors correctly point out, it has interesting practical implications, such as the possibility to perform FLIM imaging of calcium using widely available laser wavelengths, or to reduce background autofluorescence, which can be a significant problem in FLIM.<br /> - The screening was simple and rationally guided, demonstrating that, at least for this class of sensors, a careful choice of screening positions is an excellent strategy to obtain variants with large FLIM responses without the need of high-throughput screening.<br /> - The description of the methodologies is very complete and accurate, greatly facilitating the reproduction of the results by others, or the adoption of similar methods. This is particularly true for the description of the experimental conditions for optimal screening of sensor variants in lysed bacterial cultures.<br /> - The photophysical characterization is very thorough and complete, and the vast amount of data reported in the supporting information is a valuable reference for other researchers willing to attempt a similar sensor development strategy. Particularly well done is the characterization of the brightness in cells, and the comparison on multiple parameters with existing sensors.<br /> - Overall, G-CaFLITS displays excellent properties for a FLIM sensor: very large lifetime change, bright emission in both forms and independence from pH in the physiological range.

Weaknesses:

- The paper demonstrates the application of G-CaFLITS in various cellular sub-compartments without providing direct evidence that the sensor's response is not affected by the targeting. Showing at least that the lifetime values in the saturated state are similar in all compartments would improve the robustness of the claims.<br /> - In some cases, the interpretation of the results is not fully convincing, leaving alternative hypotheses as a possibility. This is particularly the case for the claim of the origin of the strongly reduced brightness of G-CaFLITS in Drosophila. The explanation of the intensity changes of G-CaFLITS also shows some inconsistency with the basic photophysical characterization.<br /> - While the claims generally appear robust, in some cases they are conveyed with a lack of precision. Several sentences in the introduction and discussion could be improved in this regard. Furthermore, the use of the signal-to-noise ratio as a means of comparison between sensors appears to be imprecise, since it is dependent on experimental conditions.

Joint Public Review:

Summary of the work:

This manuscript defines the differential stress response signaling induced by nuclear and cytoplasmic protein misfolding. To accomplish this, the authors used superfolder GFP fused to a destabilized FKBP protein-bearing targeting signal for cytosolic or nuclear localization. When cells were grown in the presence of the ligand Shield-1, this protein was stable, allowing fluorescence of the GFP protein. Upon removal of Shield-1, the FKBP protein is unfolded targeting the entire fusion protein to proteasomal degradation. Using this approach, they performed RNAseq to probe similarities and differences in transcriptional responses to the accumulation of unfolded proteins in the cytosol or nucleus. As expected, many of the pathways upregulated in both datasets involved protein homeostasis pathways such as the proteasome and cytosolic chaperones. The increase in proteasome subunits correlated with the stabilization of Nrf1 under these conditions, suggesting that protein misfolding might induce proteasome subunits through an Nrf1-dependent mechanism, but this was not explicitly tested. In contrast, the authors report that the p53-dependent transcriptional response was selectively induced by protein misfolding stress in the nucleus, but not the cytosol. Deletion of p53 blocked this increase, indicating that this response is attributable to p53 stabilization. The increased p53 transcriptional activity corresponded with the stabilization of p53 and its target p21 in cells subjected to nuclear but not cytosolic protein misfolding stress. Using a reporter of nuclear proteasome activity, they show that nuclear proteasome activity is reduced in cells following protein misfolding stress in the nucleus, indicating that the stabilization of p53 (and other transcription factors such as NRF1) might be attributed to reduced proteasomal degradation. Additionally, the authors showed that nuclear misfolding stress also induces cell cycle arrest. However, this effect was not dependent on p53 deletion, indicating that this is mediated by other unknown mechanisms.

Major strengths and weaknesses of the methods and results:

The findings reported here define specific transcriptional outputs induced by targeted protein misfolding stress in the nucleus and cytosol, revealing new insights into the organelle-specific stress signaling. The approach is interesting and effective at revealing cellular responses induced by compartment-specific protein misfolding stress.

One major weakness of the study is the lack of mechanistic follow-up for the transcriptional study. For example, what is the mechanistic basis for p53 stabilization by nuclear-destabilized domain (Nuc DD)? Is this entirely caused by diminished nuclear degradation activity as shown in Figure 6 or are there additional factors to be considered? If limited proteasome degradation capacity is the main reason for p53 upregulation, wouldn't the authors also see stabilization of other short-lived transcription factors? The fact that Nrf1 and Nrf2 are also stabilized by Nuc DD is consistent with the authors' hypothesis. On the other hand, if Nuc DD also affects other short-lived transcription factors such as c-fos or c-myc via proteasome inhibition, why did the gene expression analysis only pick up the p53 pathway as the one differentially regulated by Nuc DD? Would this imply that only p53 is specifically targeted by the nuclear proteasome, whereas other short-lived transcription factors are degraded either by the cytosolic proteasome or by both nuclear and cytosolic proteasome like Nrf1? Is there any evidence in the literature that supports this speculation? Additionally, how does Nuc DD affect the UPS system in the nucleus? Does it clog the proteasome directly or affect other assisting factors like chaperones or ubiquitinating enzymes? Lastly, it isn't clear what the functional implications of p53 stabilization would be for cells subjected to nuclear protein misfolding stress, particularly as the small effect on cell cycle arrest is not dependent on p53. In the end, the lack of mechanistic and/or functional follow-up reduces the overall importance of this manuscript. While the reviewers do not expect the authors to answer all these questions by experiments, additional work/clarifications/discussions along these lines would significantly improve the paper (see the recommendations).

Another major weakness is the lack of statistical analysis (SA) to better support their conclusions. In fact, no SA was provided for many figures even though the authors tried to make many comparisons.

The failure of the DD reporter to mount a significant heat shock response was puzzling. The presence of non-native proteins is the primary trigger for the heat shock response, but the authors acknowledge that inducible chaperones such as Hspa1a/b and Hsp90aa1 were not significantly changed in their system (page 8). Could this suggest a problem with the approach? What exactly is the nature of the stress mounted by Nuc DD?

The cell cycle data presented in Figure 5 is less robust, particularly as the p53 data in panels C and D was collected only once.

The Western blot data shown in Figure 6 does not have quantification to show how representative the blot is and how robust the changes in protein levels are over time. Western blots are known to be variable with different replicates and therefore the authors need to mention the number of biological repeats represented by the blot.

Reviewer #1 (Public review):

This manuscript discusses from a theory point of view he mechanisms underlying the formation of specialized or mixed factories. To investigate this, a chromatin polymer model was developed to mimic the chromatin binding-unbinding dynamics of various complexes of transcription factors (TFs).

The model revealed that both specialized (i.e., demixed) and mixed clusters can emerge spontaneously, with the type of cluster formed primarily determined by cluster size. Non-specific interactions between chromatin and proteins were identified as the main factor promoting mixing, with these interactions becoming increasingly significant as clusters grow larger.

These findings, observed in both simple polymer models and more realistic representations of human chromosomes, reconcile previously conflicting experimental results. Additionally, the introduction of different types of TFs was shown to strongly influence the emergence of transcriptional networks, offering a framework to study transcriptional changes resulting from gene editing or naturally occurring mutations.

Overall I think this is an interesting paper discussing a valuable model of how chromosome 3D organisation is linked to transcription. I would only advise the authors to polish and shorten their text to better highlight their key findings and make it more accessible to the reader.

Reviewer #1 (Public review):

Summary:

This study has preliminarily revealed the role of ACVR2A in trophoblast cell function, including its effects on migration, invasion, proliferation, and clonal formation, as well as its downstream signaling pathways.

Strengths:

The use of multiple experimental techniques, such as CRISPR/Cas9-mediated gene knockout, RNA-seq, and functional assays (e.g., Transwell, colony formation, and scratch assays), is commendable and demonstrates the authors' effort to elucidate the molecular mechanisms underlying ACVR2A's regulation of trophoblast function. The RNA-seq analysis and subsequent GSEA findings offer valuable insights into the pathways affected by ACVR2A knockout, particularly the Wnt and TCF7/c-JUN signaling pathways.

Weaknesses:

The molecular mechanisms underlying this study require further exploration through additional experiments. While the current findings provide valuable insights into the role of ACVR2A in trophoblast cell function and its involvement in the regulation of migration, invasion, and proliferation, further validation in both in vitro and in vivo models is needed. Additionally, more experiments are required to establish the functional relevance of the TCF7/c-JUN pathway and its clinical significance, particularly in relation to pre-eclampsia. Additional techniques, such as animal models and more advanced clinical sample analyses, would help strengthen the conclusions and provide a more comprehensive understanding of the molecular pathways involved.

Reviewer #1 (Public review):

Summary:

Meissner-Bernard et al present a biologically constrained model of telencephalic area of adult zebrafish, a homologous area to the piriform cortex, and argue for the role of precisely balanced memory networks in olfactory processing.

This is interesting as it can add to recent evidence on the presence of functional subnetworks in multiple sensory cortices. It is also important in deviating from traditional accounts of memory systems as attractor networks. Evidence for attractor networks has been found in some systems, like in the head direction circuits in the flies. However, the presence of attractor dynamics in other modalities, like sensory systems, and their role in computation has been more contentious. This work contributes to this active line of research in experimental and computational neuroscience by suggesting that, rather than being represented in attractor networks and persistent activity, olfactory memories might be coded by balanced excitation-inhibitory subnetworks.

Strengths:

The main strength of the work is in: (1) direct link to biological parameters and measurements, (2) good controls and quantification of the results, and (3) comparison across multiple models.

(1) The authors have done a good job of gathering the current experimental information to inform a biological-constrained spiking model of the telencephalic area of adult zebrafish. The results are compared to previous experimental measurements to choose the right regimes of operation.<br /> (2) Multiple quantification metrics and controls are used to support the main conclusions, and to ensure that the key parameters are controlled for - e.g. when comparing across multiple models.<br /> (3) Four specific models (random, scaled I / attractor, and two variant of specific E-I networks - tuned I and tuned E+I) are compared with different metrics, helping to pinpoint which features emerge in which model.

In the revised manuscript, the authors have also:<br /> (a) made a good effort to provide a mechanistic explanation of their results (especially on the mechanism underlying medium amplification in specific E/I network models);<br /> (b) performed a systematic analysis of the parameter space by changing different parameters of E and I neurons (specifically showing that different time constants of E and I neurons do not change the results and therefore the main effects result from connectivity);<br /> (c) added further analysis and discussion on the potential functional and computational significance of balanced specific E-I subnetworks.

These additions substantially strengthen the study, presenting compelling evidence for how networks with specific E-I structure can underpin olfactory processing and memory representations. The findings have potential implications that extend beyond the olfactory system and may be applicable to other neural systems and species.

Reviewer #1 (Public review):

Summary:

It is evident that studying leukocyte extravasation in vitro is a challenge. One needs to include physiological flow, culture cells and isolate primary immune cells. Timing is of utmost importance and a reproducible setup is essential. Extra challenges are met when extravasation kinetics in different vascular beds is required, e.g., across the blood-brain barrier. In this study, the authors describe a reliable and reproducible method to analyze leukocyte TEM under physiological flow conditions, including this analysis. That the software can also detect reverse TEM is a plus.

Strengths:

It is quite a challenge to get this assay reproducible and stable, in particular as there is flow included. Also for the analysis, there is currently no clear software analysis program, and many labs have their own methods. This paper gives the opportunity to unify the data and results obtained with this assay under label-free conditions. This should eventually lead to more solid and reproducible results.

Also, the comparison between manual and software analysis is appreciated.

Reviewer #1 (Public review):

Summary:

In this manuscript, Shibata describes a method to assess rapidly fluctuating CpG sites (fCpGs) from single-cell methylation sequencing (sc-MeSeq) data. Assuming that fCpGs are largely consistent over time with changes induced by inheritable events during replication, the author infers lineage relationships in available brain-derived sc-MeSeq. Supplementing current lineage tracing through genomic and mitochondrial mosaic variants is an interesting concept that could supplement current work or allow additional lineage analysis in existing data.

However, the author failed to convincingly show the power of fCpG analysis to determine lineages in the human brain. While the correlation with cellular division and distinction of cell types appears plausible and strong, the application to detect specific lineages is less convincing. Aspects of this might be due to a lack of clarity in presentation and erroneous use of developmental concepts. However, without addressing these problems it is challenging for a reader to come to the same conclusions as the author.

On the flip side, this novel application of fCpGs will allow the re-use of existing sc-MeSeq to infer additional features that were previously unavailable, once the biological relevance has been further elucidated.

Strengths:

• Novel re-analysis application of methylation data to infer the status of fCpGs and the use as a lineage marker<br /> • Application of this method to an innovative existing data set to benchmark this framework against existing developmental knowledge

Weaknesses:

• Inconsistent or erroneous use of neurodevelopmental concepts which hinders appropriate interpretation of the results.<br /> • Somewhat confusing presentation at times which makes it hard to judge the value of this novel approach.

Reviewer #1 (Public review):

Summary:

The authors intended to investigate the earliest mechanisms enabling self-prioritization, especially in the attention. Combining a temporal order judgement task with computational modelling based on the Theory of Visual Attention (TVA), the authors suggested that the shapes associated with the self can fundamentally alter the attentional selection of sensory information into awareness. This self-prioritization in attentional selection occurs automatically at early perceptual stages. Furthermore, the processing benefits obtained from attentional selection via self-relatedness and physical salience were separated from each other.

Strengths:

The manuscript is written in a way that is easy to follow. The methods of the paper are very clear and appropriate.

Comments on revisions:

The authors clearly showed the relationship between attention and self-prioritization.

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.

Weaknesses:

The robust data showing the quality of this model at the transcriptomic level can be strengthened with confirmation at protein and functional levels. The authors were unable to investigate the effects of Adar1-KO using ER-Hoxb8 cells and instead had to rely on a mouse model with a patient-specific Adar1 mutation (Adar1 D1113H). Additionally, ER-Hoxb8-derived microglia do not express Sall1, a key marker of microglia, which limits their fidelity as a full microglial replacement, as has been rightfully pointed out in the discussion.

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. While Adar1-KO macrophages do not engraft well, the success of TLR4-KO line highlights the model's potential for investigating other genes. 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.

Joint Public Review:

Summary:

Microfossils from the Paleoarchean Eon represent the oldest evidence of life, but their nature has been strongly debated among scientists. To resolve this, the authors reconstructed the lifecycles of Archaean organisms by transforming a Gram-positive bacterium into a primitive lipid vesicle-like state and simulating early Earth conditions. They successfully replicated all morphologies and life cycles of Archaean microfossils and studied cell degradation processes over several years, finding that encrustation with minerals like salt preserved these cells as fossilized organic carbon. Their findings suggest that microfossils from 3.8 to 2.5 billion years ago were likely liposome-like protocells with energy conservation pathways but without regulated morphology.

Strengths:

The authors have crafted a compelling narrative about the morphological similarities between microfossils from various sites and proliferating wall-deficient bacterial cells, providing detailed comparisons that have never been demonstrated in this detail before. The extensive number of supporting figures is impressive, highlighting numerous similarities. While conclusively proving that these microfossils are proliferating protocells morphologically akin to those studied here is challenging, we applaud this effort as the first detailed comparison between microfossils and morphologically primitive cells.

Summary of reviewer comments on this revision:

Each of the original reviewers evaluated the revised manuscript and were complimentary about how the authors addressed their original concerns. One reviewer added: "It is a thought-provoking manuscript that will be well received." We encourage readers of this version of the paper to consider the original reviewer comments and the authors' responses: https://elifesciences.org/reviewed-preprints/98637/reviews

Reviewer #1 (Public review):

Summary:

In their manuscript, the authors provide compelling evidence that stimulus-frequency otoacoustic emission (SFOAE) phase-gradient delays predict the sharpness (quality factors) of auditory-nerve-fiber (ANF) frequency tuning curves in budgerigars. In contrast with mammals, neither SFOAE- nor ANF-based measures of cochlear tuning match the frequency dependence of behavioral tuning in this species of parakeet. Although the reason for the discrepant behavioral results (taken from previous studies) remains unexplained, the present data provide significant and important support for the utility of otoacoustic estimates of cochlear tuning, a methodology previously explored only in mammals.

Strengths:

* The OAE and ANF data appear solid and believable. (The behavioral data are taken from previous studies.)

* No other study in birds (and only a single previous study in mammals) has combined behavioral, auditory-nerve, and otoacoustic estimates of cochlear tuning in a single species.

* SFOAE-based estimates of cochlear tuning now avoid possible circularity and were are obtained by assuming that the tuning ratio estimated in chicken applies also to the budgerigar.

Weaknesses:

* In mammals, accurate prediction of neural Q_ERB from otoacoustic N_SFOAE involves the application of species-invariance of the tuning ratio combined with an attempt to compensate for possible species differences in the location of the so-called apical-basal transition (for a review, see Shera & Charaziak, Cochlear frequency tuning and otoacoustic emissions. Cold Spring Harb Perspect Med 2019; 9:pii a033498. doi: 10.1101/cshperspect.a033498; in particular, the text near Eq. 2 and the value of CFa|b).

Despite this history, the manuscript makes no mention of the apical-basal transition, its possible role in birds, or why it was ignored in the present analysis. As but one result, the comparative discussion of the tuning ratio (paragraph beginning on lines 383) is incomplete and potentially misleading. Although the paragraph highlights differences in the tuning ratio across groups, perhaps these differences simply reflect differences in the value of CFa|b. For example, if the cochlea of the budgerigar is assumed to be entirely "apical" in character (so that CFa|b is around 7-8 kHz), then the budgerigar tuning ratios appear to align remarkably well with those previously obtained in mammals (see Shera et al 2010, Fig 9).

* For the most part, the authors take previous behavioral results in budgerigar at face value, attributing the discrepant behavioral results to hypothesized "central specializations for the processing of masked signals". But before going down this easy road, the manuscript would be stronger if the authors discussed potential issues that might affect the reliability of the previous behavioral literature. For example, the ANF data show that thresholds rise rapidly above about 5 kHz. Might the apparent broadening of the behavioral filters arise as<br /> a consequence of off-frequency listening due to the need to increase signal levels at these frequencies? Or perhaps there are other issues. Inquiring readers would appreciate an informed discussion.

Reviewer #1 (Public review):

Tu, Wen, et al. investigated the activity of mPFC putative glutamatergic neurons during a probabilistic threat discrimination and avoidance learning task using miniaturized GRIN lens implantation and single-photon calcium imaging in freely moving mice. In conjunction with this cellular recording, they employed channelrhodopsin-mediated optogenetic excitation of terminals from basal forebrain cholinergic projection neurons coupled to the delivery of an air puff on either of two maze paths with differential threat probability. The authors found that the optogenetic manipulation altered mPFC encoding of outcomes and disrupted animals' behavioral adaptation. Over the course of multiple learning sessions, optogenetically stimulated mice lagged behind control animals in resolving the differential threat probabilities on the two paths and making adaptive choices. In particular, the animals with optogenetic stimulation of cholinergic terminals were significantly more likely to switch to the path with higher threat probability after having just gotten a rare air puff on the generally "safer" path. Combined with data from a deterministic version of the task showing that optogenetically stimulated mice could behaviorally discriminate between the paths appropriately under such circumstances, these results suggest an impairment in the experimental animals' ability to make use of threat history over multiple trials. This comparison of probabilistic and deterministic versions of the same task is a highlight of this paper, representing a thoughtfulness about what information can be gleaned from such variations in the design of behavioral experiments that is all too often lacking. These data are timely in contributing to an ongoing discussion in the field about the role of phasic cholinergic signaling to the cortex, about which relatively little is known.

While the ensemble recording of mPFC neurons during the task appears to be reliable and well-designed and the behavioral effects of the optogenetic stimulation are convincing, some major weaknesses of the paper limit its usefulness to others in the field:

(1) Optogenetic excitation of presynaptic terminals can lead to antidromic action potentials that alter the firing properties of the target cell (see the excellent review on challenges of and strategies for presynaptic optogenetic experiments Rost et al., Nat Neurosci 2022). To their credit, the authors explicitly acknowledge this fact, but they believe that the only alternative possibility is that their intervention could lead to increased acetylcholine release at collateral projections in other prefrontal subregions. In fact, we do not know that the mechanism mediating the behavioral changes observed involves acetylcholine at all, as many ChAT+ basal forebrain neurons co-transmit using GABA (Saunders et al., Nature, 2015; Saunders et al., eLife, 2015; Granger et al., Neuropharmacology, 2016). A very useful internal control, which is recommended by Rost et al. for such presynaptic excitation experiments, would be to locally infuse nicotinic or muscarinic cholinergic antagonists into the mPFC in an attempt to reverse the optogenetically induced deficit; this would resolve whether the effect is indeed mediated by cholinergic neurotransmission and if it is specific to the mPFC.

(2) In a similar vein, the fact that LED illumination in the no-opsin control group appears to increase activity in prefrontal neurons (Figure 2C) and, moreover, has a functional effect in disrupting location-selective cellular activity to a similar extent as in the ChrimsonR group (Figure S3) is inadequately explained and cause for concern. Although the authors argue that the degree or "robustness" of puff-evoked activity was significantly greater in the ChrimsonR group as compared to fluorophore-only controls, their statistical test for demonstrating this is the Kolmogorov-Smirnov test (Figure 2D), thus showing that the two samples likely are drawn from different distributions but little else.

(3) Throughout the paper, the authors rely heavily on the Kolmogorov-Smirnov and binomial tests (Figures 2D, 3, 4D, S3, S4) to compare distributions in this manner, but it is unclear to me why these would be the most appropriate statistical tests for what they seek to demonstrate. Given the holistic nature of these tests in comparing the shape and spread of distributions, I am concerned that they might be inflating the significance of the differences between groups. Even if the authors were seeking a nonparametric statistical test, which most likely would be quite appropriate, there are nonparametric versions of ANOVA that they could use (e.g. Kruskal-Wallis, Friedman). Indeed, in much of this data set a repeated measures statistical analysis would seem to be called for, whereas the Kolmogorov-Smirnov test assumes that the two samples must be independent of each other. The most notable example of this premise being violated is in Figure 3, where data from the same cell populations in the same animals are being compared between experimental days and across various trial types.

Reviewer #1 (Public review):

This manuscript presents an interesting exploration of the potential activation mechanisms of DLK following axonal injury. While the experiments are beautifully conducted and the data are solid, I feel that there is insufficient evidence to fully support the conclusions made by the authors.

In this manuscript, the authors exclusively use the puc-lacZ reporter to determine the activation of DLK. This reporter has been shown to be induced when DLK is activated. However, there is insufficient evidence to confirm that the absence of reporter activation necessarily indicates that DLK is inactive. As with many MAP kinase pathways, the DLK pathway can be locally or globally activated in neurons, and the level of DLK activation may depend on the strength of the stimulation. This reporter might only reflect strong DLK activation and may not be turned on if DLK is weakly activated. The results presented in this manuscript support this interpretation. Strong stimulation, such as axotomy of all synaptic branches, caused robust DLK activation, as indicated by puc-lacZ expression. In contrast, weak stimulation, such as axotomy of some synaptic branches, resulted in weaker DLK activation, which did not induce the puc-lacZ reporter. This suggests that the strength of DLK activation depends on the severity of the injury rather than the presence of intact synapses. Given that this is a central conclusion of the study, it may be worthwhile to confirm this further. Alternatively, the authors may consider refining their conclusion to better align with the evidence presented.

As noted by the authors, DLK has been implicated in both axon regeneration and degeneration. Following axotomy, DLK activation can lead to the degeneration of distal axons, where synapses are located. This raises an important question: how is DLK activated in distal axons? The authors might consider discussing the significance of this "synapse connection-dependent" DLK activation in the broader context of DLK function and activation mechanisms.

Reviewer #1 (Public review):

This study identifies two behavioral processes that underlie learned pathogen avoidance behavior in C. elegans: exiting and re-entry of pathogenic bacterial lawns. Long-term behavioral tracking indicates that animals increase the prevalence of both behaviors over long-term exposure to the pathogen Pseudomonas aeruginosa. Using an optogenetic silencing screen, the authors identify groups of neurons, whose activity regulates lawn occupancy. Surprisingly, they find that optogenetic inhibition of neurons during only the first two hours of pathogen exposure can establish subsequent long-term changes in pathogen aversion. By leveraging a compressed sensing approach, the authors define a set of neurons involved in either lawn exit or lawn re-entry behavior using a constrained set of transgenic lines that drive Arch-3 expression in overlapping groups of neurons. They then measure the calcium activity of the candidate neurons involved in lawn re-entry in freely moving animals using GCaMP, and observe a reduction in their neural activity after exposure to pathogen. Optogenetic inhibition of AIY and SIA neurons during acute pathogen exposure in naïve animals delays lawn entry whereas activating these neurons in animals previously exposed to pathogen enhances lawn entry, albeit transiently.

This work is missing experiments and analyses that are necessary to substantiate their claims. Although the authors convincingly show that neuronal inhibition experiments during pathogen exposure reveal separable groups of neurons controlling pathogenic lawn exiting and re-entry, their methods to validate these results at single neuron cell-type resolution lack rigor.

In Figure 4, the authors claim that the reduction in calcium activity in cells of interest following pathogen exposure encodes pathogen experience. However, they make no effort to correlate the observed decreased activity with concomitant shifts in increased immobility (decreased forward locomotion) or the increased age of the worms since pathogen exposure began (24 hours have elapsed), either of which could easily explain these results. A better comparison would be between age-matched naive animals and animals exposed to pathogen. More to the point, we are interested in the involvement of these neurons' activity patterns with the behavioral motifs associated with lawn exits and re-entries, so examining these activity patterns in the absence of any pathogen before or after long-term pathogen exposure yields little insight into their relevant signaling roles. To substantiate the authors' claims, a better experiment would measure these neurons' calcium activity during lawn exits and re-entries in naive and post-exposed age-matched worms.

In Figure 5, the authors attempt to show that manipulating AIY and SIA/SIB neuronal activity controls pathogenic lawn re-entry behavior. Although they show that inhibiting these neurons in naive animals increases latency to enter pathogenic lawns, they never test the effect of neuronal inhibition in post-exposed animals. Instead they activate these neurons using channelrhodopsin, whereby they observe an increase in lawn entry and exit behavior, indicative of high forward locomotion speed. Although suggestive, neither of these experiments prove these neurons' involvement in pathogenic lawn re-entry behavior following pathogen exposure. To rigorously test the hypothesis that AIY and SIA/SIB neurons are required to sustain higher latency to lawn re-entry following pathogen exposure, the authors should perform neuronal inhibition experiments in post-pathogen-exposed animals as well and compare the results. The interpretation of this figure is further complicated by the fact that Npr-4::ChR2 animals express ChR2 in AIY in addition to SIA/SIB neurons: experiments that calculated lawn re-entry rates in Npr-4::ChR2 activation in post-exposed animals may include the known effect of stimulating AIY alone (Fig. 5J) since no discernible attempt at structured illumination to limit excitation to SIA/SIB neurons was made in these animals (Fig. 5 K, L).

This work raises the interesting possibility that different sets of neurons control lawn exit and lawn re-entry behaviors following pathogen exposure. However, the authors never directly test this claim. To rigorously show this, the authors would need to show that lawn-exit promoting neurons (CEPs, HSNs, RIAs, RIDs, SIAs) are dispensable for lawn re-entry behavior and that lawn re-entry promoting neurons (AVK, SIA, AIY, MI) are dispensable for lawn exit behavior in pathogen-exposed animals. The authors identify AVK neurons as important for modulating lawn re-entry behavior by brief inhibition at the start of pathogen exposure but fail to find that these neurons are required for increased latency to re-entry in naïve animals (Fig. 5D). Recent work from Marquina-Solis et al (2024) shows that chronic silencing of these neurons delays pathogen lawn leaving, due to impaired release of flp-1 neuropeptide. Authors may wish to connect their work more closely with the existing literature by investigating the behavioral process by which AVK contributes to lawn evacuation.

Reviewer #1 (Public review):

Summary:

In this manuscript by Lopez-Blanch and colleagues, 21 microexons are selected for a deep analysis of their impacts on behavior, development, and gene expression. The authors begin with a systematic analysis of microexon inclusion and conservation in zebrafish and use these data to select 21 microexons for further study. The behavioral, transcriptomic, and morphological data presented are for the most part convincing. Furthermore, the discussion of the potential explanations for the subtle impacts of individual microexon deletions versus loss-of-function in srrm3 and/or srrm4 is quite comprehensive and thoughtful. One major weakness: data presentation, methods, and jargon at times affect readability / might lead to overstated conclusions. However, overall this manuscript is well-written, easy to follow, and the results are of broad interest.

Strengths:

(1) The study uses a wide variety of techniques to assess the impacts of microexon deletion, ranging from assays of protein function to regulation of behavior and development.

(2) The authors provide comprehensive analyses of the molecular impact of their microexon deletions, including examining how host-gene and paralog expression is affected.

Weaknesses / Major Points:

(1) According to the methods, it seems that srrm3 social behavior is tested by pairing a 3mpf srrm3 mutant with a 30dpf srrm3 het. Is this correct? The methods seem to indicate that this decision was made to account for a slower growth rate of homozygous srrm3 mutant fish. However, the difference in age is potentially a major confound that could impact the way that srrm3 mutants interact with hets and the way that srrm3 mutants interact with one another (lower spread for the ratio of neighbour in front value, higher distance to neighbour value). This reviewer suggests testing het-het behavior at 3 months to provide age-matched comparisons for del-del, testing age-matched rather than size-matched het-del behavior, and also suggests mentioning this in the main text / within the figure itself so that readers are aware of the potential confound.

(2) Referring to srrm3+/+; srrm4-/- controls for double mutant behavior as "WT for simplicity" is somewhat misleading. Why do the authors not refer to these as srrm4 single mutants?

(3) It's not completely clear how "neurally regulated" microexons are defined / how they are different from "neural microexons"? Are these terms interchangeable?

(4) Overexpression experiments driving srrm3 / srrm4 in HEK293 cells are not described in the methods.

(4) Suggest including more information on how neurite length was calculated. In representative images, it appears difficult to determine which neurites arise from which soma, as they cross extensively. How was this addressed in the quantification?

Reviewer #1 (Public review):

Summary:

The authors aimed to investigate the cellular mechanisms underlying place field formation (PFF) in hippocampal CA1 pyramidal cells by performing in vivo two-photon calcium imaging in head-restrained mice navigating a virtual environment. Specifically, they sought to determine whether BTSP-like (behavioral time scale synaptic plasticity) events, characterized by large calcium transients, are the primary mechanism driving PFFs or if other mechanisms also play a significant role. Through their extensive imaging dataset, the authors found that while BTSP-like events are prevalent, a substantial fraction of new place fields are formed via non-BTSP-like mechanisms. They further observed that large calcium transients, often associated with BTSP-like events, are not sufficient to induce new place fields, indicating the presence of additional regulatory factors (possibly local dendritic spikes).

Strengths

The study makes use of a robust and extensive dataset collected from 163 imaging sessions across 45 mice, providing a comprehensive examination of CA1 place-cell activity during navigation in both familiar and novel virtual environments. The use of two-photon calcium imaging allows the authors to observe the detailed dynamics of neuronal activity and calcium transients, offering insights into the differences between BTSP-like and non-BTSP-like PFF events. The study's ability to distinguish between these two mechanisms and analyze their prevalence under different conditions is a key strength, as it provides a nuanced understanding of how place fields are formed and maintained. The paper supports the idea that BTSP is not the only driving force behind PFF, and other mechanisms are likely sufficient to drive PFF, and BTSP events may also be insufficient to drive PFF in some cases. The longer-than-usual virtual track used in the experiment allowed place cells to express multiple place fields, adding a valuable dimension to the dataset that is typically lacking in similar studies. Additionally, the authors took a conservative approach in classifying PFF events, ensuring that their findings were not confounded by noise or ambiguous activity.

Weaknesses

Despite the impressive dataset, there are several methodological and interpretational concerns that limit the impact of the findings. Firstly, the virtual environment appears to be poorly enriched, relying mainly on wall patterns for visual cues, which raises questions about the generalizability of the results to more enriched environments. Prior studies have shown that environmental enrichment can significantly influence spatial coding, and it would be important to determine how a more immersive VR environment might alter the observed PFF dynamics. Secondly, the study relies on deconvolution methods in some cases to infer spiking activity from calcium signals without in vivo ground truth validation. This introduces potential inaccuracies, as deconvolution is an estimate rather than a direct measure of spiking, and any conclusions drawn from these inferred signals should be interpreted with caution. Thirdly, the figures would benefit from clearer statistical annotations and visual enhancements. For example, several plots lack indicators of statistical significance, making it difficult for readers to assess the robustness of the findings. Furthermore, the use of bar plots without displaying underlying data distributions obscures variability, which could be better visualized with violin plots or individual data points. The manuscript would also benefit from a more explicit breakdown of the proportion of place fields categorized as BTSP-like versus non-BTSP-like, along with clearer references to figures throughout the results section. Lastly, the authors' interpretation of their data, particularly regarding the sufficiency of large calcium transients for PFF induction, needs to be more cautious. Without direct confirmation that these transients correspond to actual BTSP events (including associated complex spikes and calcium plateau potentials), concluding that BTSP is not necessary or sufficient for PFF formation is speculative.

Reviewer #1 (Public review):

Summary:

Zhang et al. addressed the question of whether advantageous and disadvantageous inequality aversion can be vicariously learned and generalized. Using an adapted version of the ultimatum game (UG), in three phases, participants first gave their own preference (baseline phase), then interacted with a "teacher" to learn their preference (learning phase), and finally were tested again on their own (transfer phase). The key measure is whether participants exhibited similar choice preferences (i.e., rejection rate and fairness rating) influenced by the learning phase, by contrasting their transfer phase and baseline phase. Through a series of statistical modeling and computational modeling, the authors reported that both advantageous and disadvantageous inequality aversion can indeed be learned (Study 1), and even be generalised (Study 2).

Strengths:

This study is very interesting, it directly adapted the lab's previous work on the observational learning effect on disadvantageous inequality aversion, to test both advantageous and disadvantageous inequality aversion in the current study. Social transmission of action, emotion, and attitude have started to be looked at recently, hence this research is timely. The use of computational modeling is mostly appropriate and motivated. Study 2, which examined the vicarious inequality aversion in conditions where feedback was never provided, is interesting and important to strengthen the reported effects. Both studies have proper justifications to determine the sample size.

Weaknesses:

Despite the strengths, a few conceptual aspects and analytical decisions have to be explained, justified, or clarified.

INTRODUCTION/CONCEPTUALIZATION<br /> (1) Two terms seem to be interchangeable, which should not, in this work: vicarious/observational learning vs preference learning. For vicarious learning, individuals observe others' actions (and optionally also the corresponding consequence resulting directly from their own actions), whereas, for preference learning, individuals predict, or act on behalf of, the others' actions, and then receive feedback if that prediction is correct or not. For the current work, it seems that the experiment is more about preference learning and prediction, and less so about vicarious learning. The intro and set are heavily around vicarious learning, and later the use of vicarious learning and preference learning is rather mixed in the text. I think either tone down the focus on vicarious learning, or discuss how they are different. Some of the references here may be helpful: Charpentier et al., Neuron, 2020; Olsson et al., Nature Reviews Neuroscience, 2020; Zhang & Glascher, Science Advances, 2020

EXPERIMENTAL DESIGN<br /> (2) For each offer type, the experiment "added a uniformly distributed noise in the range of (-10 ,10)". I wonder what this looks like? With only integers such as 25:75, or even with decimal points? More importantly, is it possible to have either 70:30 or 90:10 option, after adding the noise, to have generated an 80:20 split shown to the participants? If so, for the analyses later, when participants saw the 80:20 split, which condition did this trial belong to? 70:30 or 90:10? And is such noise added only to the learning phase, or also to the baseline/transfer phases? This requires some clarification.

(3) For the offer conditions (90:10, 70:30, 50:50, 30:70, 10:90) - are they randomized? If so, how is it done? Is it randomized within each participant, and/or also across participants (such that each participant experienced different trial sequences)? This is important, as the order especially for the learning phase can largely impact the preference learning of the participants.

STATISTICAL ANALYSIS & COMPUTATIONAL MODELING<br /> (4) In Study 1 DI offer types (90:10, 70:30), the rejection rate for DI-AI averse looks consistently higher than that for DI averse (ie, the blue line is above the yellow line). Is this significant? If so, how come? Since this is a between-subject design, I would not anticipate such a result (especially for the baseline). Also, for the LME results (eg, Table S3), only interactions were reported but not the main results.

(5) I do not particularly find this analysis appealing: "we examined whether participants' changes in rejection rates between Transfer and Baseline, could be explained by the degree to which they vicariously learned, defined as the change in punishment rates between the first and last 5 trials of the Learning phase." Naturally, the participants' behavior in the first 5 trials in the learning phase will be similar to those in the baseline; and their behavior in the last 5 trials in the learning phase would echo those at the transfer phase. I think it would be stronger to link the preference learning results to the change between the baseline and transfer phase, eg, by looking at the difference between alpha (beta) at the end of the learning phase and the initial alpha (beta).

(6) I wonder if data from the baseline and transfer phases can also be modeled, using a simple Fehr-Schimdt model. This way, the change in alpha/beta can also be examined between the baseline and transfer phase.

(7) I quite liked Study 2 which tests the generalization effect, and I expected to see an adapted computational modeling to directly reflect this idea. Indeed, the authors wrote, "[...] given that this model [...] assumes the sort of generalization of preferences between offer types [...]". But where exactly did the preference learning model assume the generalization? In the methods, the modeling seems to be only about Study 1; did the authors advise their model to accommodate Study 2? The authors also ran simulation for the learning phase in Study 2 (Figure 6), and how did the preference update (if at all) for offers (90:10 and 10:90) where feedback was not given? Extending/Unpacking the computational modeling results for Study 2 will be very helpful for the paper.

Reviewer #1 (Public review):

Summary:

Optical blur is characterized by contrast losses and phase shifts that alter the local relationship between the component spatial frequencies in the image. The eye experiences optical blur on several occasions - for instance, physiologically, when the focus state of the eye does not match the optical vergence demand and, in cases of pathologies like keratoconus where the cornea gets progressively distorted leading to degraded retinal image quality. Recalibration of the visual system to suprathreshold contrast losses arising from the optical blur and the mechanisms that may underlie such a recalibration have been well-researched. This study by Barbot et al presents convincing evidence that the visual system could also recalibrate itself to the phase distortions experienced with optical blur. This was demonstrated, in principle, on a small number of participants with normal vision but with induced blur (?? experienced psychophysical observers) and in a few keratoconic patients using their state-of-the-art adaptive optics apparatus. In the former cohort, known magnitudes of radially asymmetric blur from a vertical coma were induced while participants judged the position of a compound grating target that shifted in predictable ways with the induction of blur. Immediate exposure to images blurred with such higher-order aberrations resulted in position shifts that were consistent with optical theory, but prolonged exposure to such blur resulted in the position shift returning to veridical perception (albeit, not completely). When the blur was removed after the adaptation phase, after effects of the position offset were noticed. In the keratoconic cohort, such position offsets were observed even when the eye was completely corrected for optical degradation. These results are discussed in the context of the perception of real-world targets, the underlying neurophysiology, and what it means to space perception in disease conditions like keratoconus.

Strengths:

A clear hypothesis, a parameterized experimental space, rigor of optical correction and psychophysical judgements, and clarity in the explanation of results are the major strengths of the paper. Additional strengths include the control experiments to address confounders and the additional analyses shown in the supplementary section to rule out analytical inconsistencies in explaining the results.

Weaknesses:

The small sample size (especially in the keratoconic cohort) may be a limitation of the study. While the experiments conducted in this study are meant to demonstrate a basic visual phenomenon, that only 6 keratoconic patients were included in the study precludes the results from being extrapolated to the heterogeneity of disease presentation. It must, however, be noted that these are difficult experiments to conduct, and getting multiple participants to agree to such an experiment is not an easy task.

Second, the analysis shown in Figure 6C relating the magnitude of habitual higher-order RMS to the absolute PSE shift is not convincing. The PSE's were both positive and negative in the KC patients. The direction of the phase shift experienced by the patient (i.e., positive or negative shift in the PSE) should also be determined by the pattern of HOA's in their eyes. Simply comparing the absolute magnitudes does not make sense. Would it be possible to convolve the compound grating with the PSF obtained from each patient and predict which direction should the PSE shift? This prediction can then be compared with the observed shift in the PSE's.

A third weakness of the study may be the assumption that the phase recalibration in keratoconic cohort may be eye-specific. That is, if the participant has dissimilar severities of keratoconus, the probed eye's aberration profile may determine the phase profile that the eye is calibrated to. I am not sure to what extent this assumption is valid. Further, under natural viewing, the pupil size will change with light intensity and accommodative state and this will, in turn, determine the optical quality of the eye. Given this, it is not clear what will the visual system recalibrate itself to, when the phase shifts in the retinal image may keep changing from the underlying blur profile in the retina. Also, if the disease is progressive in nature (in their cohort, the authors indicate that the disease did not progress), the calibration state should also constantly change. What is the time scale of such a calibration and could there be multiple states of such adaptation remains to be explored. This, of course, is not a weakness of the present study, but an open question for the future.

Finally, one additional experiment could have been performed (this is good to have information and certainly not a necessity). What is the wavefront profile of a few keratoconic patients that participated in the study, used as the adaptation profile in the 2nd experiment (as opposed to a fixed level of coma)? Would a 60-min paradigm result in adapted states that will result in phase shifts matching what is experienced by keratoconic eyes (see Marella et al., Vis Res, 2024 for a similar induced experiment for studying the impact of phase shifts on visual and stereoacuities)?

Reviewer #1 (Public review):

Summary:

This study uses information from the UK Biobank and aims to investigate the role of BMI on various health outcomes, with a focus on differences by sex. They confirm the relevance of many of the well-known associations between BMI and health outcomes for males and females and suggest that associations for some endpoints may differ by sex. Overall their conclusions appear supported by the data. The significance of the observed sex variations will require confirmation and further assessment.

Strengths:

This is one of the first systematic evaluations of sex differences between BMI and health outcomes.

The hypothesis that BMI may be associated with health differentially based on sex is relevant and even expected. As muscle is heavier than adipose tissue, and as men typically have more muscle than women, as a body composition measure BMI is sometimes prone to classifying even normal weight/muscular men as obese, while this measure is more lenient when used in women.

Confirmation of the many well-known associations is as expected and attests to the validity of their approach.

Demonstration of the possible sex differences is interesting, with this work raising the need for further study.

Weaknesses:

Many of the statistical decisions appeared to target power at the expense of quality/accuracy. For example, they chose to use self-reported information rather than doctor diagnoses for disease outcomes for which both types of data were available.

Despite known problems and bias arising from the use of one sample approach, they chose to use instruments from the UK Biobank instead of those available from the independent GIANT GWAS, despite the difference in sample size being only marginally greater for UKB for the context. With the way the data is presented, it is difficult to assess the extent to which results are compatible across approaches.

The approach to multiple testing correction appears very lenient, although the lack of accuracy in the reporting makes it difficult to know what was done exactly. The way it reads, FDR correction was done separately for men, and then for women (assuming that the duplication in tests following stratification does not affect the number of tests). In the second stage, they compared differences by sex using Z-test, apparently without accounting for multiple testing.

Presentation lacks accuracy in a few places, hence assessment of the accuracy of the statements made by the authors is difficult.

Conclusion "These findings highlight the importance of retaining a healthy BMI" is rather uninformative, especially as they claim that for some attributes the effects of BMI may be opposite depending on sex/gender.

Reviewer #1 (Public review):

Summary:

This study highlights the strengths of using predictive computational models to inform C. elegans screening studies of compounds' effects on aging and lifespan. The authors primarily focus on all-trans retinoic acid (atRA), one of the 5 compounds (out of 16 tested) that extended C. elegans lifespan in their experiments. They show that atRA has positive effects on C. elegans lifespan and age-related health, while it has more modest and inconsistent effects (i.e., some detrimental impacts) for C. briggsae and C. tropicalis. In genetic experiments designed to evaluate contributing mediators of lifespan extension with atRA exposure, it was found that 150 µM of atRA did not significantly extend lifespan in akt-1 or akt-2 loss-of-function mutants, nor in animals with loss of function of aak-2, or skn-1 (in which atRA had toxic effects); these genes appear to be required for atRA-mediated lifespan extension. hsf-1 and daf-16 loss-of-function mutants both had a modest but statistically significant lifespan extension with 150 µM of atRA, suggesting that these transcription factors may contribute towards mediating atRA lifespan extension, but that they are not individually required for some lifespan extension. RNAseq assessment of transcriptional changes in day 4 atRA-treated adult wild-type worms revealed some interesting observations. Consistent with the study's genetic mutant lifespan observations, many of the atRA-regulated genes with the greatest fold-change differences are known regulated targets of daf-2 and/or skn-1 signaling pathways in C. elegans. hsf-1 loss-of-function mutants show a shifted atRA transcriptional response, revealing a dependence on hsf-1 for ~60% of the atRA-downregulated genes. On the other hand, RNAseq analysis in aak-2 loss-of-function mutants revealed that aak-2 is only required for less than a quarter of the atRA transcriptional response. All together, this study is proof of the concept that computational models can help optimize C. elegans screening approaches that test compounds' effects on lifespan, and provide comprehensive transcriptomic and genetic insights into the lifespan-extending effects of all-trans retinoic acid (atRA).

Strengths:

(1) A clearly described and well-justified account describes the approach used to prioritize and select compounds for screening, based on using the top candidates from a published list of computationally ranked compounds (Fuentealba et al., 2019) that were cross-referenced with other bioinformatics publications to predict anti-aging compounds, after de-selecting compounds previously evaluated in C. elegans as per the DrugAge database. 16 compounds were tested at 4-5 different concentrations to evaluate effects on C. elegans lifespan.

(2) Robust experimental design was undertaken evaluating the lifespan effects of atRA, as it was tested on three strains each of C. elegans, C. briggsae, and C. tropicalis, with trial replication performed at three distinct laboratories. These observations extended beyond lifespan to include evaluations of health metrics related to swimming performance.

(3) In-depth analyses of the RNAseq data of whole-worm transcriptional responses to atRA revealed interesting insights into regulator pathways and novel groups of genes that may be involved in mediating lifespan-extension effects (e.g., atRA-induced upregulation of sphingolipid metabolism genes, atRA-upregulation of genes in a poorly-characterized family of C. elegans paralogs predicted to have kinase-like activity, and disproportionate downregulation of collagen genes with atRA).

Weaknesses:

(1) The authors' computational-based compound screening approach led to a ~30% prediction success rate for compounds that could extend the median lifespan of C. elegans. However, follow-up experiments on the top compounds highlighted the fact that some of these observed "successes" could be driven by indirect, confounding effects of these compounds on the bacterial food source, rather than direct beneficial effects on C. elegans physiology and lifespan. For instance, this appeared to be the case for the "top" hit of propranolol; other compounds were not tested with metabolically inert or killed bacteria. In addition, there are no comparative metrics provided to compare this study's ~30% success rate to screening approaches that do not use computational predictions.

(2) Transcriptomic analyses of atRA effects were extensive in this study, but evaluations and discussions of non-transcriptional effects of key proposed regulators (such as AMPK) were limited. For instance, non-transcriptional effects of aak-2/AMPK might account for its requirement for mediating lifespan extension effects, since aak-2 was not required for a major proportion of atRA transcriptional responses.

Reviewer #1 (Public review):

Summary:

The authors demonstrate that female Spodoptera littoralis moths prefer to oviposit on well-watered tomato plants and avoid drought-stressed plants. The study then recorded the sounds produced by drought-stressed plants and found that they produce 30 ultrasonic clicks per minute. Thereafter, the authors tested the response of female S. littoralis moths to clicks with a frequency of 60 clicks per minute in an arena with and without plants and in an arena setting with two healthy plants of which one was associated with 60 clicks per minute. These experiments revealed that in the absence of a plant, the moths preferred to lay eggs on the side of the area in which the clicks could be heard, while in the presence of a plant the S. littoralis females preferred to oviposit on the plant where the clicks were not audible. In addition, the authors also tested the response of S. littoralis females in which the tympanic membrane had been pierced making the moths unable to detect the click sounds. As hypothesised, these females placed their eggs equally on both sites of the area. Finally, the authors explored whether the female oviposition choice might be influenced by the courtship calls of S. littoralis males which emit clicks in a range similar to a drought-stressed tomato plant. However, no effect was found of the clicks from ten males on the oviposition behaviour of the female moths, indicating that the females can distinguish between the two types of clicks. Besides these different experiments, the authors also investigated the distribution of egg clusters within a longer arena without a plant, but with a sugar-water feeder. Here it was found that the egg clusters were mostly aggregated around the feeder and the speaker producing 60 clicks per minute. Lastly, video tracking was used to observe the behaviour of the area without a plant, which demonstrated that the moths gradually spent more time at the arena side with the click sounds.

Strengths:

This manuscript is very interesting to read and the possibility that female moths might use sound as an additional sensory modality during host-searching is exciting and very relevant to the field of insect-plant interactions.

Weaknesses:

The study addresses a very interesting question by asking whether female moths incorporate plant acoustic signals into their oviposition choice, unfortunately, I find it very difficult to judge how big the influence of the sound on the female choice really is as the manuscript does not provide any graphs showing the real numbers of eggs laid on the different plants, but instead only provides graphs with the Bayesian model fittings for each of the experiments. In addition, the numbers given in the text seem to be relatively similar with large variations e.g. Figure 1B3: 1.8 {plus minus} 1.6 vs. 1.1 {plus minus} 1.0. Furthermore, the authors do not provide access to any of the raw data or scripts of this study, which also makes it difficult to assess the potential impact of this study. Hence, I would very much like to encourage the authors to provide figures showing the measured values as boxplots including the individual data points, especially in Figure 1, and to provide access to all the raw data underlying the figures.

Regarding the analysis of the results, I am also not entirely convinced that each night can be taken as an independent egg-laying event, as the amount of eggs and the place were the eggs are laid by a female moth surely depends on the previous oviposition events. While I must admit that I am not a statistician, I would suggest, from a biological point of view, that each group of moths should be treated as a replicate and not each night. I would therefore also suggest to rather analyse the sum of eggs laid over the different consecutive nights than taking the eggs laid in each night as an independent data point.

Furthermore, it did not become entirely clear to me why a click frequency of 60 clicks per minute was used for most experiments, while the plants only produce clicks at a range of 30 clicks per minute. Independent of the ecological relevance of these sound signals, it would be nice if the authors could provide a reason for using this frequency range. Besides this, I was also wondering about the argument that groups of plants might still produce clicks in the range of 60 clicks per minute and that the authors' tests might therefore still be reasonable. I would agree with this, but only in the case that a group of plants with these sounds would be tested. Offering the choice between two single plants while providing the sound from a group of plants is in my view not the most ecologically reasonable choice. It would be great if the authors could modify the argument in the discussion section accordingly and further explore the relevance of different frequencies and dB-levels.

Finally, I was wondering how transferable the findings are towards insects and Lepidopterans in general. Not all insects possess a tympanic organ and might therefore not be able to detect the plant clicks that were recorded. Moreover, I would imagine that generalist herbivorous like Spodoptera might be more inclined to use these clicks than specialists, which very much rely on certain chemical cues to find their host plants. It would be great if the authors would point more to the fact that your study only investigated a single moth species and that the results might therefore only hold true for S. littoralis and closely related species, but not necessary for other moth species such as Sphingidae or even butterflies.

Reviewer #1 (Public review):

The origin recognition complex (ORC) is an essential loading factor for the replicative Mcm2-7 helicase complex. Despite ORC's critical role in DNA replication, there have been instances where the loss of specific ORC subunits has still seemingly supported DNA replication in cancer cells, endocycling hepatocytes, and Drosophila polyploid cells. Critically, all tested ORC subunits are essential for development and proliferation in normal cells. This presents a challenge, as conditional knockouts need to be generated, and a skeptic can always claim that there were limiting but sufficient ORC levels for helicase loading and replication in polyploid or transformed cells. That being said, the authors have consistently pushed the system to demonstrate replication in the absence or extreme depletion of ORC subunits.

Here, the authors generate conditional ORC2 mutants to counter a potential argument with prior conditional ORC1 mutants that Cdc6 may substitute for ORC1 function based on homology. They also generate a double ORC1 and ORC2 mutant, which is still capable of DNA replication in polyploid hepatocytes. While this manuscript provides significantly more support for the ability of select cells to replicate in the absence or near absence of select ORC subunits, it does not shed light on a potential mechanism.

The strengths of this manuscript are the mouse genetics and the generation of conditional alleles of ORC2 and the rigorous assessment of phenotypes resulting from limiting amounts of specific ORC subunits. It also builds on prior work with ORC1 to rule out Cdc6 complementing the loss of ORC1.

The weakness is that it is a very hard task to resolve the fundamental question of how much ORC is enough for replication in cancer cells or hepatocytes. Clearly, there is a marked reduction in specific ORC subunits that is sufficient to impact replication during development and in fibroblasts, but the devil's advocate can always claim minimal levels of ORC remaining in these specialized cells.

The significance of the work is that the authors keep improving their conditional alleles (and combining them), thus making it harder and harder (but not impossible) to invoke limiting but sufficient levels of ORC. This work lays the foundation for future functional screens to identify other factors that may modulate the response to the loss of ORC subunits.

This work will be of interest to the DNA replication, polyploidy, and genome stability communities.

Reviewer #1 (Public review):

Summary:

The authors in this manuscript performed scRNA-seq on a cohort of 15 early-stage cervical cancer patients with a mixture of adeno- and squamous cell carcinoma, HPV status, and several samples that were upstaged at the time of surgery. From their analyses they identified differential cell populations in both immune and tumour subsets related to stage, HPV status, and whether a sample was adenocarcinoma or squamous cell. Putative microenvironmental signaling was explored as a potential explanation for their differential cell populations. Through these analyses the authors also identified SLC26A3 as a potential biomarker for later stage/lymph node metastasis which was verified by IHC and IF. The dataset is likely useful for the community. The accuracy and clarity have been improved from the previous version, and additional immunofluorescence supporting the existence of their proposed cluster is now present. That said, there remain some issues with the strength of some claims (particularly in the abstract and results sections) and some of the cell type definitions.

Strengths

The dataset could be useful for the community<br /> SLC26A3 could potentially be a useful marker to predict lymph node metastasis with further study

Weaknesses

Casual language is used in the abstract around immunosuppressive microenvironment and signal cross-talk between Epi_10_CYSTM1 cluster and Tregs. The data show localization that supports a possible interaction and probable cytokines, but functional experiments would be needed to establish causality.

In the description of the single cell data processing there is no mention of batch effect correction. Given that many patients were analyzed, and no mention was made of pooling or deconvolution, it must be assumed these were run separately which invariably leads to batch effects. Given the good overlays across patients some batch correction must have been performed. How was batch effect correction performed?

While statistics were added to the clinical correlates, it would appear that single variables are being assessed one at a time by chi-squared analysis. This ignores the higher order structure of the data and the correlations between some variables resulting in potentially spurious findings. This is compounded as some categories had below 5 observations violating the assumptions of a chi-squared test.

The description of all analytical steps remains quite truncated. While the inclusion of annotated code is useful, a full description of which tools were used, with which settings, and why each were chosen, is a minimum needed to properly interpret the results. This is as important in a mainly analytical paper as the experimental parameters.

Validation of the clustering results remains a problem. The only details provided are that FindClusters was used. This depends on a manual choice of multiple parameters including the k-nearest neighbours included, whether Louvain or Leiden clustering is used, the resolution parameter, and others (how many variable genes/PCs etc...). Why were these parameters selected, how do you know that you're not over or under-clustering.

The cluster Epi_10_CYSTM1 remains somewhat problematic. None of the additional data supports its existence outside of the single patient who has cells from that population. Additionally, it falls well outside of any of the other Epithelial cells to the point that drawing it as part of a differentiation order doesn't even make sense. Indeed, most of the upregulated pathways in this cluster appear to be related to class II antigen presentation which would fit better with a dendritic cell/macrophage than an epithelial cell. While the IF at the end does support the existence of the cluster, numbers are still very limited, and this doesn't have data on the antigen presenting function. At the least a strong disclaimer should be included in the text that this population is essentially exclusive to one sample in the scRNA data.

The linkage between the cluster types and IHC for prediction of lymph node metastasis is tenuous. Most of the strongly cluster associated markers were not predictive despite their clusters being theoretically enriched. This inability to recognize the clusters in additional samples using alternative methods does not give confidence that these clusters are robust. SLC26A3 being associated with upstaging may very well be a useful marker, however, given the lack of association of the other markers, it may be premature to say this is due to the same Epi_10_CYSTM1 cluster.

There are multiple issues in the classification of T cells and neutrophils. In the analysis of T cell subset, all CD4+ T cells are currently scored as Tregs, what happened to the T-helper cells? Additionally, Activated T and Cytotoxic T both seem to contain CD8+ cells, but all their populations have equivalent expression of the activation marker CD69. Moreover, the "Cytotoxic" ones also express TIGIT, HAVCR2 and LAG3 which are generally exhaustion markers. For neutrophils, several obviously different clusters have been grouped together (Neu_1 containing two diametrically opposite cell clouds being an obvious example).

Again in the CellChat section of the results causal language is being repeatedly used. These are just possible interactions, not validated ones. While the co-localization in the provided IF images certainly supports the co-localization, this still is only correlative and doesn't prove causality.

Minor Issues<br /> The sentence "However, due to the low morbidity of ADC, in-depth investigations are insufficient" could be misinterpreted. Morbidity generally refers to the severity or health burden rather than the frequency of cases, though it's true in some studies prevalence is used for the overall impact of the disease on a population and referred to as morbidity. In this instance though, "incidence" or "prevalence" would be clearer word choices.

The previous rebuttal states that clusters/cell type calls were refined to eliminate issues such as epithelial cells creeping into the T cell cluster, however, the cell %s have not been altered according to the change tracking. Shouldn't all the %s have been altered even if only slightly?

Joint Public Review:

Reviewers thought that the authors addressed some, but not all the concerns raised in the previous round of a review.

Strengths: The authors employed a battery of next-generation sequencing and crosslinking techniques (e.g., Quick-irCLIP, APA-Seq, and Ribo-Seq) to describe a previously unappreciated binding of eIF3 to the 3'UTRs of the mRNAs. It is also shown that eIF3:3'UTR binding occurs in the vicinity of poly(A) tail of mRNAs that are actively translated in neuronal progenitor cells derived from human pluripotent stem cells. Collectively, these findings provide evidence for the role of eIF3 in regulating translation from the 3'UTR end of the mRNA.

Weaknesses: In addition to these clear strengths of the article, some weaknesses were observed pertinent to the lack of mechanistic data. It was therefore thought that the experiments aiming to dissect the mechanisms of eIF3 binding to 3'UTRs and their impact on translation warrant future studies. Finally, establishing the impact of the proposed eIF3:3'UTR binding mechanism of translational regulation on cellular fate is required to further support the biological importance of the observed phenomena. It was found that this should also be addressed in the follow up studies.

Reviewer #1 (Public Review):

Summary:

The authors performed two-sample MR combined with sensitivity analyses and colocalization to test the effect of PUFA on cerebral aneurysms. They found that genetically predicted omega-3 and DHA decreased the risk for intracranial aneurysm (IA) and subarachnoid haemorrhage (SAH) but not for unruptured IA (uIA).

Strengths:

PUFA on the risk of cerebral aneurysms is of clinical importance; the authors performed multiple sensitivity analyses to ensure MR fulfils its assumptions.

Reviewer #1 (Public review):

The findings of Ziolkowska and colleagues show that a specific projection from the nucleus reuniens of the thalamus (RE) to dorsal CA1 of the hippocampus plays an important role in fear extinction learning in male and female mice. In and of itself, this is not a new finding. Yet, the potential novelty and excitement comes from the authors' identification of structural alterations from RE projecting neurons to the specific stratum lacunosum moleculare subregion of CA1 after learning. The authors use a range of anatomical and functional approaches to demonstrate structural synaptic changes in dorsal CA1 that parallel the necessary role of RE inputs in modulating extinction learning. The significance of these findings was previously hampered by several technical shortcomings in the experimental design and interpretation. The authors adequately addressed some of the design concerns raised in the previous round, along with the interpretive critique that they couldn't localize the timing of effects to consolidation as originally claimed. Nevertheless, the authors provided an inadequate response to the concern regarding their misapplication of Ns and missing controls in one experiment.

In the previous review, a major methodological weakness in the experimental design involved the widespread misapplication of Ns used for the statistical analyses. Much of the anatomical analyses of structural synaptic changes in the RE-CA1 pathway used N = number of axons (Figs. 1, 2), N = number of dendrites (Figs. 3, 4), and N = number of sections (Fig. 7). In each instance it was recommended that N = animal number should be used. Reasons for this are as follows: this is standard practice in neuroanatomical research; using N = branch/ dendrite/ bouton/ spine number artificially inflates the statistical power and this incorrectly assumes independence of observations; using N = number of sections, etc., doesn't account for imbalances in the number of observations that vary from animal to animal that may skew group results.

In the authors' response, they generally concurred, but then they followed up with the defense that the number of items was too few in some cases, or absent in others, to permit using N = animal number. While they changed some of their data to N = animal numbers, other aspects of their data remained as-is. The description of the statistics in the figure legend is also dense and difficult to follow in places. Ns should be checked in the legend and figure to make sure they're correct, as at least one error was noted (e.g., see Fig. 2C). Overall, the authors' response falls short of the standard of rigor that helps to reinforce scientific findings from reliability and reproducibility concerns when generating more data to increase Ns (i.e., the number of animals) would have been the better choice.

Another persistent concern from the previous review is that, in the electron microscopic analyses of dendritic spines (Fig. 5), the authors only compared fear acquisition versus extinction training. One critique was that the lack of inclusion of a naïve control group made it difficult to understand how these structural synaptic changes are occurring relative to baseline. It was also noted that the authors appropriately included naïve controls in other experiments in the paper. In the revised submission the authors simply added in naïve control data to their previous histogram. It is not considered good practice to collect, process, or analyze data one group at a time, as this would be prone to cohort effects or experimental bias. These data should be discarded and the experiment should be run correctly with randomized cases in each group, or instead these data should be eliminated from the report since there is a key control group missing. Again, the nature of the authors' response perpetuates the aforementioned concern that data collection and analysis in this report may fall short of an acceptable standard of rigor.

Reviewer #1 (Public review):

Summary:

The authors are trying to develop a microscopy system that generates data output exceeding the previous systems based on huge objectives.

Strengths:

They have accomplished building such a system, with a field of view of 1.5x1.0 cm2 and a resolution of up to 1.2 um. They have also demonstrated their system performance on samples such as organoids, brain sections, and embryos.

Weaknesses:

To be used as a volumetric imaging technique, the authors only showcase the implementation of multi-focal confocal sectioning. On the other hand, most of the real biological samples were acquired under the wide-field illumination, and processed with so-called computational sectioning. Despite the claim that it improves the contrast, sometimes I felt that the images were oversharpened and the quantitative nature of these fluorescence images may be perturbed.

Reviewer #1 (Public review):

Summary:

Shen et al. conducted three experiments to study the cortical tracking of the natural rhythms involved in biological motion (BM), and whether these involve audiovisual integration (AVI). They presented participants with visual (dot) motion and/or the sound of a walking person. They found that EEG activity tracks the step rhythm, as well as the gait (2-step cycle) rhythm. The gait rhythm specifically is tracked superadditively (power for A+V condition is higher than the sum of the A-only and V-only condition, Experiments 1a/b), which is independent of the specific step frequency (Experiment 1b). Furthermore, audiovisual integration during tracking of gait was specific to BM, as it was absent (that is, the audiovisual congruency effect) when the walking dot motion was vertically inverted (Experiment 2). Finally, the study shows that an individual's autistic traits are negatively correlated with the BM-AVI congruency effect.

Strengths:

The three experiments are well designed and the various conditions are well controlled. The rationale of the study is clear, and the manuscript is pleasant to read. The analysis choices are easy to follow, and mostly appropriate.

Weaknesses:

There is a concern of double-dipping in one of the tests (Experiment 2, Figure 3: interaction of Upright/Inverted X Congruent/Incongruent). I raised this concern on the original submission, and it has not been resolved properly. The follow-up statistical test (after channel selection using the interaction contrast permutation test) still is geared towards that same contrast, even though the latter is now being tested differently. (Perhaps not explicitly testing the interaction, but in essence still testing the same.) A very simple solution would be to remove the post-hoc statistical tests and simply acknowledge that you're comparing simple means, while the statistical assessment was already taken care of using the permutation test. (In other words: the data appear compelling because of the cluster test, but NOT because of the subsequent t-tests.)

Reviewer #2 (Public review):

Summary:

Griesius et al. investigate the dendritic integration properties of two types of inhibitory interneurons in the hippocampus: those that express NDNF+ and those that express somatostatin. They found that both neurons showed supralinear synaptic integration in the dendrites, blocked by NMDA receptor blockers but not by blockers of Na+ channels. These experiments are critically overdue and very important because knowing how inhibitory neurons are engaged by excitatory synaptic input has important implications for all theories involving these inhibitory neurons.

Comments on revisions:

The authors have addressed the reviewers' comments, but haven't resolved most of the key issues.

Specifically, performing only a single uncaging experiment at a single dendritic location per cell prevents a detailed biophysical analysis of NDNF and OLM cell integration properties. A more extended exploration would have potentially addressed several of the reviewers' questions. It is particularly worrying that the authors cite cell health, dendritic blebbing, and changes in input resistance as the reason for terminating experiments after a single uncaging event. This suggests that the uncaging laser may be damaging the dendrite, potentially affecting the membrane potential directly, and overall cell health, beyond simply uncaging glutamate.

While the authors' qualitative conclusions about supra-linear integration and NMDA receptor dependency seem plausible, the limited data and potential methodological issues weaken any quantitative interpretations and comparisons between the two cell types.

Similarly, the absence of dendritic Na-spikes remains unexplained, despite reports of strong dendritic Na-currents in these cells.

Reviewer #1 (Public review):

In recent years, our understanding of the nuclear steps of the HIV-1 life cycle has made significant advances. It has emerged that HIV-1 completes reverse transcription in the nucleus and that the host factor CPSF6 forms condensates around the viral capsid. The precise function of these CPSF6 condensates is under investigation, but it is clear that the HIV-1 capsid protein is required for their formation. This study by Tomasini et al. investigates the genesis of the CPSF6 condensates induced by HIV-1 capsid, what other co-factors may be required, and their relationship with nuclear speckels (NS). The authors show that disruption of the condensates by the drug PF74, added post-nuclear entry, blocks HIV-1 infection, which supports their functional role. They generated CPSF6 KO THP-1 cell lines, in which they expressed exogenous CPSF6 constructs to map by microscopy and pull down assays of the regions critical for the formation of condensates. This approach revealed that the LCR region of CPSF6 is required for capsid binding but not for condensates whereas the FG region is essential for both. Using SON and SRRM2 as markers of NS, the authors show that CPSF6 condensates precede their merging with NS but that depletion of SRRM2, or SRRM2 lacking the IDR domain, delays the genesis of condensates, which are also smaller.

The study is interesting and well conducted and defines some characteristics of the CPSF6-HIV-1 condensates. Their results on the NS are valuable. The data presented are convincing.

I have two main concerns. Firstly, the functional outcome of the various protein mutants and KOs is not evaluated. Although Figure 1 shows that disruption of the CPSF6 puncta by PF74 impairs HIV-1 infection, it is not clear if HIV-1 infection is at all affected by expression of the mutant CPSF6 forms (and SRRM2 mutants) or KO/KD of the various host factors. The cell lines are available, so it should be possible to measure HIV-1 infection and reverse transcription. Secondly, the authors have not assessed if the effects observed on the NS impact HIV-1 gene expression, which would be interesting to know given that NS are sites of highly active gene transcription. With the reagents at hand, it should be possible to investigate this too.

Reviewer #1 (Public review):

Summary:

The study examines how pyruvate, a key product of glycolysis that influences TCA metabolism and gluconeogenesis, impacts cellular metabolism and cell size. It primarily utilizes the Drosophila liver-like fat body, which is composed of large post-mitotic cells that are metabolically very active. The study focuses on the key observations that over-expression of the pyruvate importer MPC complex (which imports pyruvate from the cytoplasm into mitochondria) can reduce cell size in a cell-autonomous manner. They find this is by metabolic rewiring that shunts pyruvate away from TCA metabolism and into gluconeogenesis. Surprisingly, mTORC and Myc pathways are also hyper-active in this background, despite the decreased cell size, suggesting a non-canonical cell size regulation signaling pathway. They also show a similar cell size reduction in HepG2 organoids. Metabolic analysis reveals that enhanced gluconeogenesis suppresses protein synthesis. Their working model is that elevated pyruvate mitochondrial import drives oxaloacetate production and fuels gluconeogenesis during late larval development, thus reducing amino acid production and thus reducing protein synthesis.

Strengths:

The study is significant because stem cells and many cancers exhibit metabolic rewiring of pyruvate metabolism. It provides new insights into how the fate of pyruvate can be tuned to influence Drosophila biomass accrual, and how pyruvate pools can influence the balance between carbohydrate and protein biosynthesis. Strengths include its rigorous dissection of metabolic rewiring and use of Drosophila and mammalian cell systems to dissect carbohydrate:protein crosstalk.

Weaknesses:

However, questions on how these two pathways crosstalk, and how this interfaces with canonical Myc and mTORC machinery remain. There are also questions related to how this protein:carbohydrate crosstalk interfaces with lipid biosynthesis. Addressing these will increase the overall impact of the study.

Reviewer #1 (Public review):

Summary:

This manuscript by Toledo and colleagues describes the generation and characterization of Y220C mice (Y217C in the mouse allele). The authors make notable findings: Y217C mice that have been backcrossed to C57Bl/6 for five generations show decreased female pup births due to exencephaly, a known defect in p53 -/- mice, and they show a correlation with decreased Xist expression, as well increased female neonatal death. They also noted similar tumor formation in Y217C/+ and p53 +/- mice, suggesting that Y217C may not function as a dominant negative. Notably, the authors find that homozygous Y217C mice die faster than p53 -/- mice and that the lymphomas in the Y217C mice were more aggressive and invasive. The authors then perform RNA seq on thymi of Y217C homozygotes compared to p53 -/-, and they suggest that these differentially expressed genes may explain the increased tumorigenesis in Y217C mice.

Strengths:

Overall, the study is well controlled and quite well done and will be of interest to a broad audience, particularly given the high frequency of the Y220C mutation in cancer (1% of all cancers, 4% of ovarian cancer).

Weaknesses:

No weaknesses were noted by this reviewer.

Reviewer #1 (Public review):

Summary:

In this report, the authors made use of a murine cell life derived from a MYC-driven liver cancer to investigate the gene expression changes that accompany the switch from normoxic to hypoxia conditions during 2D growth and the switch from 2D monolayer to 3D organoid growth under normoxic conditions. They find a significant (ca. 40-50%) overlap among the genes that are dysregulated in response to hypoxia in 2D cultures and in response to spheroid formation. Unsurprisingly, hypoxia-related genes were among the most prominently deregulated under both sets of conditions. Many other pathways pertaining to metabolism, splicing, mitochondrial electron transport chain structure and function, DNA damage recognition/repair, and lipid biosynthesis were also identified.

Major comments:

(1) Lines 239-240: The authors state that genes involved in DNA repair were identified as being necessary to maintain survival of both 2D and 3D cultures (Figure S6A). Hypoxia is a strong inducer of ROS. Thus, the ROS-specific DNA damage/recognition/repair pathways might be particularly important. The authors should look more carefully at the various subgroups of the many genes that are involved in DNA repair. They should also obtain at least a qualitative assessment of ROS and ROS-mediated DNA damage by staining for total and mitochondrial-specific ROS using dyes such as CM-H2-DCFDA and MitoSox. Actual direct oxidative damage could be assessed by immunostaining for 8-oxo-dG and related to the sub-types of DNA damage-repair genes that are induced. The centrality of DNA damage genes also raises the question as to whether the previously noted prominence of the TP53 pathway (see point 5 below) might represent a response to ROS-induced DNA damage.

(2) Because most of the pathway differences that distinguish the various cell states from one another are described only in terms of their transcriptome variations, it is not always possible to understand what the functional consequences of these changes actually are. For example, the authors report that hypoxia alters the expression of genes involved in PDH regulation but this is quite vague and not backed up with any functional or empirical analyses. PDH activity is complex and regulated primarily via phosphorylation/dephosphorylation (usually mediated by PDK1 and PDP2, respectively), which in turn are regulated by prevailing levels of ATP and ADP. Functionally, one might expect that hypoxia would lead to the down-regulation of PDH activity (i.e. increased PDH-pSer392) as respiration changes from oxidative to non-oxidative. This would not be appreciated simply by looking at PDH transcript levels. This notion could be tested by looking at total and phospho-PDH by western blotting and/or by measuring actual PDH activity as it converts pyruvate to AcCoA.

(3) Line 439: Related to the above point: the authors state: "It is likely that blockade of acetyl-CoA production by PDH knockout may force cells to use alternative energy sources under hypoxic and 3D conditions, averting the Warburg effect and promoting cell survival under limited oxygen and nutrient availability in 3D spheroids." This could easily be tested by determining whether exogenous fatty acids are more readily oxidized by hypoxic 2D cultures or spheroids than occurs in normoxic 2D cultures.

(4) Line 472: "Hypoxia induces high expression of Acaca and Fasn in NEJF10 cells indicating that hypoxia promotes saturated fatty acid synthesis...The beneficial effect of Fasn and Acaca KO to NEJF10 under hypoxia is probably due to reduction of saturated fatty acid synthesis, and this hypothesis needs to be tested in the future.". As with the preceding comment, this supposition could readily be supported directly by, for example, performing westerns blots for these enzymes and by showing that incubation of hypoxic 2D cells or spheroids converted more AcCoA into lipid.

(5) In Supplementary Figure 2B&C, the central hub of the 2D normoxic cultures is Myc (as it should well be) whereas, in the normoxic 3D, the central hub is TP53 and Myc is not even present. The authors should comment on this. One would assume that Myc levels should still be quite high given that Myc is driven by an exogenous promoter. Does the centrality of TP53 indicate that the cells within the spheroids are growth-arrested, being subjected to DNA damage and/or undergoing apoptosis?

(6) In the Materials and Methods section (lines 711-720), the description of how spheroid formation was achieved is unclear. Why were the cells first plated into non-adherent 96 well plates and then into non-adherent T75 flasks? Did the authors actually utilize and expand the cells from 144 T75 flasks and did the cells continue to proliferate after forming spheroids? Many cancer cell types will initially form monolayers when plated onto non-adherent surfaces such as plastic Petri dishes and will form spheroid-like structures only after several days. Other cells will only aggregate on the "non-adherent" surface and form spheroid-like structures but will not actually detach from the plate's surface. Have the authors actually documented the formation of true, non-adherent spheroids at 2 days and did they retain uniform size and shape throughout the collection period? The single photo in Supplementary Figure 1 does not explain when this was taken. The authors include a schematic in Figure 2A of the various conditions that were studied. A similar cartoon should be included to better explain precisely how the spheroids were generated and clarify the rationale for 96 well plating. Overall, a clearer and more concise description of how spheroids were actually generated and their appearance at different stages of formation needs to be provided.

(7) The authors maintained 2D cultures in either normoxic or hypoxic (1% O2) states during the course of their experiments. On the other hand, 3D cultures were maintained under normoxic conditions, with the assumption that the interiors of the spheroids resemble the hypoxic interiors of tumors. However, the actual documentation of intra-spheroid hypoxia is never presented. It would be a good idea for the authors to compare the degree of hypoxia achieved by 2D (1% O2) and 3D cultures by staining with a hypoxia-detecting dye such as Image-iT Green. Comparing the fluorescence intensities in 2D cultures at various O2 concentrations might even allow for the construction of a "standard curve" that could serve to approximate the actual internal O2 concentration of spheroids. This would allow the authors to correlate the relative levels of hypoxia between 2D and 3D cultures.

(8) Related to the previous 2 points, the authors performed RNAseq on spheroids only 48 hours after initiating 3D growth. I am concerned that this might not have been a sufficiently long enough time for the cells to respond fully to their hypoxic state, especially given my concerns in Point 6. Might the results have been even more robust had the authors waited longer to perform RNA seq? Why was this short time used?

(9) What happens to the gene expression pattern if spheroids are re-plated into standard tissue culture plates after having been maintained as spheroids? Do they resume 2D growth and does the gene expression pattern change back?

(10) Overall, the paper is quite descriptive in that it lists many gene sets that are altered in response to hypoxia and the formation of spheroids without really delving into the actual functional implications and/or prioritizing the sets. Some of these genes are shown by CRISPR screening to be essential for maintaining viability although in very few cases are these findings ever translated into functional studies (for example, see points 1-4 above). The list of genes and gene pathways could benefit from a better explanation and prioritization of which gene sets the authors believe to be most important for survival in response to hypoxia and for spheroid formation.

(11) The authors used a single MYC-driven tumor cell line for their studies. However, in their original paper (Fang, et al. Nat Commun 2023, 14: 4003.) numerous independent cell lines were described. It would help to know whether RNAseq studies performed on several other similar cell lines gave similar results in terms of up & down-regulated transcripts (i.e. representative of the other cell lines are NEJF10 cells).

Reviewer #1 (Public review):

Summary:

In this manuscript, the authors identified that NOLC1 was upregulated in gastric cancer samples, which promoted cancer progression and cisplatin resistance. They further found that NOLC1 could bind to p53 and decrease its nuclear transcriptional activity, then inhibit p53-mediated ferroptosis. There are several major concerns regarding the conclusions.

Strengths:

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

Weaknesses:

The major conclusions were not sufficiently supported by the results. The experiments were not conducted in a comprehensive manner.

Reviewer #1 (Public review):

Summary:

In this study, Tiang et al. explore the role of ubiquitination of non-structural protein 16 (nsp16) in the SARS-CoV-2 life cycle. nsp16, in conjunction with nsp10, performs the final step of viral mRNA capping through its 2'-O-methylase activity. This modification allows the virus to evade host immune responses and protects its mRNA from degradation. The authors demonstrate that nsp16 undergoes ubiquitination and subsequent degradation by the host E3 ubiquitin ligases UBR5 and MARCHF7 via the ubiquitin-proteasome system (UPS). Specifically, UBR5 and MARCHF7 mediate nsp16 degradation through K48- and K27-linked ubiquitination, respectively. Notably, degradation of nsp16 by either UBR5 or MARCHF7 operates independently, with both mechanisms effectively inhibiting SARS-CoV-2 replication in vitro and in vivo. Furthermore, UBR5 and MARCHF7 exhibit broad-spectrum antiviral activity by targeting nsp16 variants from various SARS-CoV-2 strains. This research advances our understanding of how nsp16 ubiquitination impacts viral replication and highlights potential targets for developing broadly effective antiviral therapies.

Strengths:

The proposed study is of significant interest to the virology community because it aims to elucidate the biological role of ubiquitination in coronavirus proteins and its impact on the viral life cycle. Understanding these mechanisms will address broadly applicable questions about coronavirus biology and enhance our overall knowledge of ubiquitination's diverse functions in cell biology. Employing in vivo studies is a strength.

Weaknesses:

While the conclusions are generally well-supported by the data, additional work is needed to confirm that NSP16 is ubiquitinated in a biologically relevant context and to better define the roles of the reported E3 ligases. Clarifications regarding aspects of data acquisition, data analysis, and text editing could notably strengthen the manuscript and its conclusions.

Reviewer #1 (Public review):

Summary:

The study "Monitoring of Cell-free Human Papillomavirus DNA in Metastatic or Recurrent Cervical Cancer: Clinical Significance and Treatment Implications" by Zhuomin Yin and colleagues focuses on the relationship between cell-free HPV (cfHPV) DNA and metastatic or recurrent cervical cancer patients. It expands the application of cfHPV DNA in tracking disease progression and evaluating treatment response in cervical cancer patients. The study is overall well-designed, including appropriate analyses.

Strengths:

The findings provide valuable reference points for monitoring drug efficacy and guiding treatment strategies in patients with recurrent and metastatic cervical cancer. The concordance between HPV cfDNA fluctuations and changes in disease status suggests that cfDNA could play a crucial role in precision oncology, allowing for more timely interventions. As with similar studies, the authors used Droplet Digital PCR to measure cfDNA copy numbers, a technique that offers ultrasensitive nucleic acid detection and absolute quantification, lending credibility to the conclusions.

Weaknesses:

Despite including 28 clinical cases, only 7 involved recurrent cervical cancer, which may not be sufficient to support some of the authors' conclusions fully. Future studies on larger cohorts could solidify HPV cfDNA's role as a standard in the personalized treatment of recurrent cervical cancer patients.

Reviewer #1 (Public review):

This work presents a self-supervised method for the segmentation of 3D cells in microscopy images, an annotated dataset, as well as a napari plugin. While the napari plugin is potentially useful, there is insufficient evidence in the manuscript to support the claim that the proposed method is able to segment cells in other light-sheet microscopy image datasets than the four specific ones used here.

I acknowledge that the revision is now more upfront about the scope of this work. However, my main point still stands: even with the slight modifications to the title, this paper suggests to present a general method for self-supervised 3D cell segmentation in light-sheet microscopy data. This claim is simply not backed up.

I still think the authors should spell out the assumptions that underlie their method early on (cells need to be well separated and clearly distinguishable from background). A subordinate clause like "often in cleared neural tissue" does not serve this purpose. First, it implies that the method is also suitable for non-cleared tissue (which would have to be shown). Second, this statement does not convey the crucial assumptions of well separated cells and clear foreground/background differences that the method is presumably relying on.

It does appear that the proposed method works very well on the four investigated datasets, compared to other pre-trained or fine-tuned models. However, it still remains unclear whether this is because of the proposed method or the properties of those specific datasets (namely: well isolated cells that are easily distinguished from the background). I disagree with the authors that a comparison to non-learning methods "is unnecessary and beyond the scope of this work". In my opinion, this is exactly what is needed to proof that CellSeg3D's performance can not be matched with simple image processing.

As I mentioned in the original review, it appears that thresholding followed by connected component analysis already produces competitive segmentations. I am confused about the authors' reply stating that "[this] is not the case, as all the other leading methods we fairly benchmark cannot solve the task without deep learning". The methods against which CellSeg3D is compared are CellPose and StarDist, both are deep-learning based methods. That those methods do not perform well on this dataset does not imply that a simpler method (like thresholding) would not lead to competitive results. Again, I strongly suggest the authors include a simple, non-learning based baseline method in their analysis, e.g.:<br /> * comparison to thresholding (with the same post-processing as the proposed method)<br /> * comparison to a normalized cut segmentation (with the same post-processing as the proposed method)

Regarding my feedback about the napari plugin, I apologize if I was not clear. The plugin "works" as far as I tested it (i.e., it can be installed and used without errors). However, I was not able to recreate a segmentation on the provided dataset using the plugin alone (see my comments in the original review). I used the current master as available at the time of the original review and default settings in the plugin.

Reviewer #1 (Public review):

Summary:

Dalal and Haddad investigated how neurons in the olfactory bulb are synchronized in oscillatory rhythms at gamma frequency. Temporal coordination of action potentials fired by projection neurons can facilitate information transmission to downstream areas. In a previous paper (Dalal and Haddad 2022, https://doi.org/10.1016/j.celrep.2022.110693), the authors showed that gamma frequency synchronization of mitral/tufted cells (MTCs) in the olfactory bulb enhances the response in the piriform cortex. The present study builds on these findings and takes a closer look at how gamma synchronization is restricted to a specific subset of MTCs in the olfactory bulb. They combined odor and optogenetic stimulations in anesthetized mice with extracellular recordings.

The main findings are that lateral synchronization of MTCs at gamma frequency is mediated by granule cells (GCs), independent of the spatial distance, and strongest for MTCs with firing rates close to 40 Hz. The authors conclude that this reveals a simple mechanism by which spatially distributed neurons can form a synchronized ensemble. In contrast to lateral synchronization, they found no evidence for the involvement of GCs in lateral inhibition of nearby MTCs.

Strengths:

Investigating the mechanisms of rhythmic synchronization in vivo is difficult because of experimental limitations for the readout and manipulation of neuronal populations at fast timescales. Using spatially patterned light stimulation of opsin-expressing neurons in combination with extracellular recordings is an elegant approach. The paper provides evidence for an activity-dependent synchronization of MTCs in gamma frequency that is mediated by GCs.

Weaknesses:

The study provides several results showing the firing of MTCs in gamma frequency range, however, direct evidence for the synchronization of MTCs in gamma frequency is missing.

Reviewer #1 (Public review):

The authors aimed to investigate how the probability of a reversal in a decision-making task is represented in cortical neurons. They analyzed neural activity in the prefrontal cortex of monkeys and units in recurrent neural networks (RNNs) trained on a similar task. Their goal was to understand how the dynamical systems that implement computation perform a probabilistic reversal learning task in RNNs and nonhuman primates.

Major strengths and weaknesses:

Strengths:

(1) Integrative Approach: The study exemplifies a modern approach by combining empirical data from monkey experiments with computational modeling using RNNs. This integration allows for a more comprehensive understanding of the dynamical systems that implement computation in both biological and artificial neural networks.

(2) The focus on using perturbations to identify causal relationships in dynamical systems is a good goal. This approach aims to go beyond correlational observations.

Weaknesses:

(1) The description of the RNN training procedure and task structure lacks detail, making it difficult to fully evaluate the methodology.

(2) The conclusion that the representation is better described by separable dynamic trajectories rather than fixed points on a line attractor may be premature.

(3) The use of targeted dimensionality reduction (TDR) to identify the axis determining reversal probability may not necessarily isolate the dimension along which the RNN computes reversal probability.

Appraisal of aims and conclusions:

The authors claim that substantial dynamics associated with intervening behaviors provide evidence against a line attractor. The conclusion that this representation is better described by separable dynamic trajectories rather than fixed points on a line attractor may be premature. The authors found that the state was translated systematically in response to whether outcomes were rewarded, and this translation accumulated across trials. This is consistent with a line attractor, where reward input bumps the state along a line. The observed dynamics could still be consistent with a curved line attractor, with faster timescale dynamics superimposed on this structure.

Likely impact and utility:

This work contributes to our understanding of how probabilistic information is represented in neural circuits and how it influences decision-making. The methods used, particularly the combination of empirical data and RNN modeling, provide a valuable approach for investigating neural computations. However, the impact may be limited by some of the methodological concerns raised.

The data and methods could be useful to the community, especially if the authors provide more detailed descriptions of their RNN training procedures and task structure. However, reverse engineering of the network dynamics was minimal. Most analyses didn't take advantage of the full access to the RNN's update equations.

Reviewer #1 (Public review):

Summary:

The authors use high-throughput gene editing technology in larval zebrafish to address whether microexons play important roles in the development and functional output of larval circuits. They find that individual microexon deletions rarely impact behavior, brain morphology, or activity, and raise the possibility that behavioral dysregulation occurs only with more global loss of microexon splicing regulation. Other possibilities exist: perhaps microexon splicing is more critical for later stages of brain development, perhaps microexon splicing is more critical in mammals, or perhaps the behavioral phenotypes observed when microexon splicing is lost are associated with loss of splicing in only a few genes.

A few questions remain:

(1) What is the behavioral consequence for loss of srrm4 and/or loss-of-function mutations in other genes encoding microexon splicing machinery in zebrafish?

(2) What is the consequence of loss-of-function in microexon splicing genes on splicing of the genes studied (especially those for which phenotypes were observed).

(3) For the microexons whose loss is associated with substantial behavioral, morphological, or activity changes, are the same changes observed in loss-of-function mutants for these genes?

(4) Do "microexon mutations" presented here result in the precise loss of those microexons from the mRNA sequence? E.g. are there other impacts on mRNA sequence or abundance?

(5) Microexons with a "canonical layout" (containing TGC / UC repeats) were selected based on the likelihood that they are regulated by srrm4. Are there other parallel pathways important for regulating the inclusion of microexons? Is it possible to speculate on whether they might be more important in zebrafish or in the case of early brain development?

Strengths:

(1) The authors provide a qualitative analysis of splicing plasticity for microexons during early zebrafish development.

(2) The authors provide comprehensive phenotyping of microexon mutants, addressing the role of individual microexons in the regulation of brain morphology, activity, and behavior.

Weaknesses:

(1) It is difficult to interpret the largely negative findings reported in this paper without knowing how the loss of srrm4 affects brain activity, morphology, and behavior in zebrafish.

(2) The authors do not present experiments directly testing the effects of their mutations on RNA splicing/abundance.

(3) A comparison between loss-of-function phenotypes and loss-of-microexon splicing phenotypes could help interpret the findings from positive hits.

Reviewer #1 (Public review):

Summary:

Previous studies have shown that treatment with 17α-estradiol (a stereoisomer of the 17β-estradiol) extends lifespan in male mice but not in females. The current study by Li et al, aimed to identify cell-specific clusters and populations in the hypothalamus of aged male rats treated with 17α-estradiol (treated for 6 months). This study identifies genes and pathways affected by 17α-estradiol in the aged hypothalamus.

Strengths:

Using single-nucleus transcriptomic sequencing (snRNA-seq) on hypothalamus from aged male rats treated with 17α-estradiol they show that 17α-estradiol significantly attenuated age-related increases in cellular metabolism, stress, and decreased synaptic activity in neurons.<br /> Moreover, sc-analysis identified GnRH as one of the key mediators of 17α-estradiol's effects on energy homeostasis. Furthermore, they show that CRH neurons exhibited a senescent phenotype, suggesting a potential side effect of the 17α-estradiol. These conclusions are supported by supervised clustering by neuropeptides, hormones, and their receptors.

Weaknesses:

However, the study has several limitations that reduce the strength of the key claims in the manuscript. In particular:

(1) The study focused only on males and did not include comparisons with females. However, previous studies have shown that 17α-estradiol extends lifespan in a sex-specific manner in mice, affecting males but not females. Without the comparison with the female data, it's difficult to assess its relevance to the lifespan.

(2) It's not known whether 17α-estradiol leads to lifespan extension in male rats similar to male mice. Therefore, it is not possible to conclude that the observed effects in the hypothalamus, are linked to the lifespan extension.

(3) The effect of 17α-estradiol on non-neuronal cells such as microglia and astrocytes is not well described (Fig.1). Previous studies demonstrated that 17α-estradiol reduces microgliosis and astrogliosis in the hypothalamus of aged male mice. Current data suggest that the proportion of oligo, and microglia were increased by the drug treatment, while the proportions of astrocytes were decreased. These data might suggest possible species differences, differences in the treatment regimen, or differences in drug efficiency. This has to be discussed.

A more detailed analysis of glial cell types within the hypothalamus in response to drug should be provided.

(4) The conclusion that CRH neurons are going into senescence is not clearly supported by the data. A more detailed analysis of the hypothalamus such as histological examination to assess cellular senescence markers in CRH neurons, is needed to support this claim.

Comments on revisions:

Some of the concerns were addressed in this revised version, and the authors responded and addressed study design limitations in both sexes/ages.

However, there are still some concerns that were not sufficiently addressed:

While the term "senescent" was changed to "stressed," some histological/ cellular validation of this phenotype is still needed.

Some discussion on the sex-specific effects of 17α-estradiol in the hypothalamus is still required. Previous studies in mice demonstrated that 17α-estradiol reduced hypothalamic microgliosis and astrogliosis in male but not female UM-HET3 mice.

Additionally, the provided analysis on astrocytes and microglia is superficial.

Reviewer #1 (Public review):

The manuscript examines the role of Naa10 in cKO animals, in immortalized neurons, and in primary neurons. Given that Naa10 mutations in humans produce defects in nervous system function, the authors used various strategies to try to find a relevant neuronal phenotype and its potential molecular mechanism.

This work contains valuable findings that suggest that the depletion of Naa10 from CA1 neurons in mice exacerbates anxiety-like behaviors. Using neuronal-derived cell lines authors establish a link between N-acetylase activity, Btbd3 binding to CapZb, and F-actin, ultimately impinging on neurite extension. The evidence demonstrating this is in most cases incomplete, since some key controls are missing and clearly described or simply because claims are not supported by the data. The manuscript also contains biochemical, co-immunoprecipitation, and proteomic data that will certainly be of value to our knowledge of the effects of protein N--acetylation in neuronal development and function.

Reviewer #1 (Public review):

Summary:

In this manuscript, Bisht et al address the hypothesis that protein folding chaperones may be implicated in aggregopathies and in particular Tau aggregation, as a means to identify novel therapeutic routes for these largely neurodegenerative conditions.

The authors conducted a genetic screen in the Drosophila eye, which facilitates the identification of mutations that either enhance or suppress a visible disturbance in the nearly crystalline organization of the compound eye. They screened by RNA interference all 64 known Drosophila chaperones and revealed that mutations in 20 of them exaggerate the Tau-dependent phenotype, while 15 ameliorated it. The enhancer of the degeneration group included 2 subunits of the typically heterohexameric prefoldin complex and other co-translational chaperones.

The authors characterized in depth one of the prefoldin subunits, Pfdn5, and convincingly demonstrated that this protein functions in the regulation of microtubule organization, likely due to its regulation of proper folding of tubulin monomers. They demonstrate convincingly using both immunohistochemistry in larval motor neurons and microtubule binding assays that Pfdn5 is a bona fide microtubule-associated protein contributing to the stability of the axonal microtubule cytoskeleton, which is significantly disrupted in the mutants.

Similar phenotypes were observed in larvae expressing Frontotemporal dementia with Parkinsonism on chromosome 17-associated mutations of the human Tau gene V377M and R406W. On the strength of the phenotypic evidence and the enhancement of the TauV377M-induced eye degeneration, they demonstrate that loss of Pfdn5 exaggerates the synaptic deficits upon expression of the Tau mutants. Conversely, the overexpression of Pfdn5 or Pfdn6 ameliorates the synaptic phenotypes in the larvae, the vacuolization phenotypes in the adult, and even memory defects upon TauV377M expression.

Strengths:

The phenotypic analyses of the mutant and its interactions with TauV377M at the cell biological, histological, and behavioral levels are precise, extensive, and convincing and achieve the aims of characterization of a novel function of Pfdn5.

Regarding this memory defect upon V377M tau expression. Kosmidis et al (2010) pmid: 20071510, demonstrated that pan-neuronal expression of TauV377M disrupts the organization of the mushroom bodies, the seat of long-term memory in odor/shock and odor/reward conditioning. If the novel memory assay the authors use depends on the adult brain structures, then the memory deficit can be explained in this manner.

If the mushroom bodies are defective upon TauV377M expression does overexpression of Pfdn5 or 6 reverse this deficit? This would argue strongly in favor of the microtubule stabilization explanation.

The discovery that Pfdn5 (and 6 most likely) affect tauV377M toxicity is indeed a novel and important discovery for the Tauopathies field. It is important to determine whether this interaction affects only the FTDP-17-linked mutations, or also WT Tau isoforms, which are linked to the rest of the Tauopathies. Also, insights on the mode(s) that Pfdn5/6 affect Tau toxicity, such as some of the suggestions above are aiming at, will likely be helpful towards therapeutic interventions.

Weaknesses:

What is unclear however is how Pfdn5 loss or even overexpression affects the pathological Tau phenotypes.

Does Pfdn5 (or 6) interact directly with TauV377M? Colocalization within tissues is a start, but immunoprecipitations would provide additional independent evidence that this is so.

Does Pfdn5 loss exacerbate TauV377M phenotypes because it destabilizes microtubules, which are already at least partially destabilized by Tau expression?<br /> Rescue of the phenotypes by overexpression of Pfdn5 agrees with this notion.

However, Cowan et al (2010) pmid: 20617325 demonstrated that wild-type Tau accumulation in larval motor neurons indeed destabilizes microtubules in a Tau phosphorylation-dependent manner.

So, is TauV377M hyperphosphorylated in the larvae?? What happens to TauV377M phosphorylation when Pfdn5 is missing and presumably more Tau is soluble and subject to hyperphosphorylation as predicted by the above?

Expression of WT human Tau (which is associated with most common Tauopathies other than FTDP-17) as Cowan et al suggest has significant effects on microtubule stability, but such Tau-expressing larvae are largely viable. Will one mutant copy of the Pfdn5 knockout enhance the phenotype of these larvae?? Will it result in lethality? Such data will serve to generalize the effects of Pfdn5 beyond the two FDTP-17 mutations utilized.

Does the loss of Pfdn5 affect TauV377M (and WTTau) levels?? Could the loss of Pfdn5 simply result in increased Tau levels? And conversely, does overexpression of Pfdn5 or 6 reduce Tau levels?? This would explain the enhancement and suppression of TauV377M (and possibly WT Tau) phenotypes. It is an easily addressed, trivial explanation at the observational level, which if true begs for a distinct mechanistic approach.

Finally, the authors argue that TauV377M forms aggregates in the larval brain based on large puncta observed especially upon loss of Pfdn5. This may be so, but protocols are available to validate this molecularly the presence of insoluble Tau aggregates (for example, pmid: 36868851) or soluble Tau oligomers as these apparently differentially affect Tau toxicity. Does Pfdn5 loss exaggerate the toxic oligomers and overexpression promotes the more benign large aggregates??

Reviewer #1 (Public review):

Summary:

This study utilized publicly available Hi-C data to ensemble a comprehensive set of breast cancer cell lines (luminal, Her2+, TNBC) with varying metastatic features to answer whether breast cancer cells would acquire organ-specific features at the 3D genome level to metastasize to that specific organ. The authors focused on lung metastasis and included several controls as the comparison including normal mammary lines, normal lung epithelial lines, and lung cancer cell lines. Due to the lower resolution at 250KB binning size, the authors only addressed the compartments (A for active compartment and B for inactive compartment) not the other 3D organization of the genome. They started by performing clustering and PCA analysis for the compartment identity and discovered that this panel of cell lines could be well separated based on Her2 and epithelial-mesenchymal features according to the compartment identity. While correlating with the transcriptomic changes, the authors noticed the existence of concordance and divergence between the compartment changes and transcriptomic changes. The authors then switched gears to tackle the core question of metastatic organotropism to the lung. They discovered a set of "lung permissive compartment changes" and concluded that "lung metastatic breast cancer cell lines acquire lung-like genome architecture" and "organotropic 3D genome changes match target organ more than an unrelated organ". To prove the latter point, the authors enlisted an additional non-breast cancer cell line (prostate cancer) in the setting of brain metastasis. This is a piece of pure dry computational work without wet bench experiments.

Strengths:

The authors embarked on an ambitious journey to seek the answer regarding 3D genome changes predisposing to metastatic organotropism. The authors succeeded in the assembly of a comprehensive panel of breast cancer cell lines and the aggregation of the 3D genome structure data to conduct a hypothesis-driven computation analysis. The authors also achieved in including proper controls representing normal non-cancerous epithelium and the end organ of interest. The authors did well in the citation of relevant references in 3D genome organization and EMT.

Weaknesses:

(1) The authors should clearly indicate how they determine the patterns of spread of the breast cancer cell lines being utilized in this manuscript. How did the authors arrive at the conclusion that certain cell lines would be determined as "localized spread" and "metastatic tropism to the lung"? This definition is crucial, and I will explain why.

Todd Golub's team from the Broad Institute of MIT and Harvard published "A metastasis map of human cancer cell lines" to exhaustively create a first-generation metastasis map (MetMap) that reveals organ-specific patterns of metastasis. (By the way, this work was not cited in the reference in this manuscript.) The MetMap Explorer (https://depmap.org/metmap/vis-app/index.html) is a public resource that could be openly accessed to visualize the metastatic potential of each cell line as determined by the in vivo barcoding approach as described in the MetMap paper in the format of petal plots. 5 organs were tested in the MetMap paper, including brain, lung, liver, kidney, and bone. The authors would discover that some of the organ-specific metastasis patterns defined in the MetMap Explorer would be different from the authors' classification. For example, the authors defined MCF7 as a line as lung metastatic, and rightly so the MetMap charted a signal towards lung with low penetrance and low metastatic potential. The authors defined ZR751 as a line with localized spread, however, the MetMap charted a signal towards the kidney with low penetrance and low metastatic potential, the signal strength similar to the lung metastasis in MCF7. A similar argument could be made for T47D. The TNBC line MDA-MB-231 is indeed highly metastatic, however, in MetMap data, its metastasis is not only specific to the lung but towards all 5 organs with high penetrance and metastatic potential. The 2 lung cancer cell lines mentioned in this study, A549 and H460, the authors defined them as localized spread to the lung. However, the MetMap data clearly indicated that A549 and H460 are highly metastatic to all 5 organs with high penetrance and high metastatic potential.

Since results will vary among different experimental models testing metastatic organotropism, (intra-cardiac injection was the metastasis model being adopted in the MetMap), the authors should state more clearly which experimental model system served as the basis for their definition of organ-specific metastasis. In my opinion, this is the most crucial first step for this entire study to be sound and solid.

(2) Figure 1b: The authors found that "MDA-MB-231 cells were grouped with the lung carcinoma cells. This implies that the genome organization of this cell line is closer to that of lung cells than to other breast epithelial cell lines.". In fact, another TNBC line BT549 was also clustered under the same clade. So this clade consisted of normal-like and highly metastatic lines. Therefore, the authors should be mindful of the fact that the compartment features might not directly link to metastasis (or even metastatic organotropism).

(3) Figure 3: In the text, the authors stated, "To further investigate this result, we examined the transcription status of genes that changed compartment across the EMT spectrum and, conversely, the compartment status of genes that changed transcription (Fig. 3b, c, and d)". However, it was not apparent in the figure that the cell lines were arranged according to an EMT spectrum. Also, the clustering heatmaps did not provide sufficient information regarding the genes with concordant/divergent compartments vs transcription changes. It would be more informative if the authors could spend more effort in annotating these genes/pathways.

(4) Figure 4: The title of the subheading of this section was 'Lung metastatic breast cancer cell lines acquire lung-like genome architecture". Echoing my comments in point 1, I am a bit hesitant to term it as "lung metastatic" but rather "metastatic' in general since cell lines such as MDA-MD-231 do metastasize to other organs as well. However, I do get the point that the definition of "lung metastasis" is derived from the common metastasis features among the cell lines here (MCF7, T47D, SKBR3, MDA-MB-231).

There might be another argument about whether the "lung" carcinoma cell lines can be considered "localized" since they are also capable of metastasizing to other organs. In a way, what the authors probably were trying to leverage here is the "tissue" identity of that organ. Having said this, in addition to showing the "lung permissive changes", the authors should show the "breast identity conservation" as well. Because this section started to deal with the concept of "tissue/lineage identify", the authors should also clarify whether these breast cancer cell lines capable of making lung metastasis are also preserving their original tissue identity from the compartment features (which would most likely be the case).

(5) Rest of the sections: The authors started to claim that the organ-specific metastasis permissive compartmental features mimic the destinated end organ. The authors utilized additional non-breast cancer cell lines (prostate cancer cell lines LNCaP as localized and DU145 as brain metastatic) in brain metastasis to strengthen this claim. (DU145 in MetMap again is highly metastatic to lung, brain, and kidney). However, this makes one wonder that for cell lines that are capable of metastasizing to multiple organ sites (eg. MDA-MB-231, DU145, A459, H460), does it mean that they all acquire the permissive features for all these organs? This scenario is clinically relevant in Stage 4 patients who often present with not only one metastatic lesion in one single organ but multiple metastatic lesions in more than one organ (eg. concomitant liver and lung metastasis). Do the authors think that there might be different clones having different tropism-permissive 3D genome features or there might be evolutionary trajectory in this?

In my opinion, to further prove this point, the authors might need to consider doing in vivo experiments to collect paired primary and organ-specific metastatic samples to look at the 3D genome changes.

(6) Technically, the study utilized public Hi-C data without generating new Hi-C data. The resolution of the Hi-C data for compartments was set at 250KB as the binning size indicating that the Hi-C data was at lower resolution so it might not be ideal to address other 3D genome architecture changes such as TADs or long-range loops. It is therefore unknown whether there might be permissive TAD/loop changes associated with organotropism and this is the limitation of this study.

(7) In the final sentence of the discussion the authors stated "Overall, our results suggest that genome spatial compartment changes can help encode a cell state that favors metastasis (EMT)". The "metastasis (EMT)" was in fact not clearly linked inside the manuscript. The authors did not provide a strong link between metastasis and EMT in their result description. It is also unclear whether the EMT-associated compartment identity would also correlate with the organotropic compartment identity.