Reviewer #3 (Public Review):

Summary:

The receptor tyrosine kinase Anaplastic Lymphoma Kinase (ALK) in humans is nervous system expressed and plays an important role as an oncogene. A number of groups have been studying ALK signalling in flies to gain mechanistic insight into its various roles. In flies, ALK plays a critical role in development, particularly embryonic development and axon targeting. In addition, ALK also was also shown to regulate adult functions including sleep and memory. In this manuscript, Sukumar et al., used a suite of molecular techniques to identify downstream targets of ALK signalling. They first used targeted DamID, a technique that involves a DNA methylase to RNA polymerase II, so that GATC sites in close proximity to PolII binding sites are marked. They performed these experiments in wild type and ALK loss of function mutants (using an Alk dominant negative ALkDN), to identify Alk responsive loci. Comparing these loci with a larval single cell RNAseq dataset identified neuroendocrine cells as an important site of Alk action. They further combined these TaDa hits with data from RNA seq in Alk Loss and Gain of Function manipulations to identify a single novel target of Alk signalling - a neuropeptide precursor they named Sparkly (Spar) for its expression pattern. They generated a mutant allele of Spar, raised an antibody against Spar, and characterised its expression pattern and mutant behavioural phenotypes including defects in sleep and circadian function.

Strengths:

The molecular biology experiments using TaDa and RNAseq were elegant and very convincing. The authors identified a novel gene they named Spar. They also generated a mutant allele of Spar (using CrisprCas technology) and raised an antibody against Spar. These experiments are lovely, and the reagents will be useful to the community. The paper is also well written, and the figures are very nicely laid out making the manuscript a pleasure to read.

Weaknesses:

The manuscript has improved substantially in the revision. Yet, some concerns remain around the genetics and behavioural analysis which is incomplete and confusing. The authors generated a novel allele of Spar - Spar ΔExon1 and examined sleep and circadian phenotypes of this allele and of RNAi knockdown of Spar. The RNAi knockdown is a welcome addition. However, the authors only show one parental control the GAL4 / +, but leave out the other parental control i.e. the UAS RNAi / + e.g. in Fig. 9. It is important to show both parental controls.

Further, the sleep and circadian characterisation could be substantially improved. It is unclear how sleep was calculated - what program was used or what the criteria to define a sleep bout was. In the legend for Fig 8c, it says sleep was shown as "percentage of time flies spend sleeping measured every 5min across a 24h time span". Sleep in flies is (usually) defined as at least 5 min of inactivity. With this definition, I'm not sure how one can calculate the % time asleep in a 5 min bin! Typically people use 30min or 60min bins. The sleep numbers for controls also seem off to me e.g. in Fig. 8H and H' average sleep / day is ~100. Is this minutes of sleep? 100 min / day is far too low, is it a typo? The same applies to Figure 8, figure supplement 2. Other places e.g. Fig 8 figure supplement 1, avg sleep is around 1000 min / day. The numbers for sleep bouts are also too low to me e.g. in Fig 9 number of sleep bouts avg around 4, and in Fig. 8 figure supplement 2 they average 1 sleep bout. There are several free software packages to analyse sleep data (e.g. Sleep Mat, PMID 35998317, or SCAMP). I would recommend that the authors reanalyse their data using one of these standard packages that are used routinely in the field. That should help resolve many issues.

The circadian anticipatory activity analyses could also be improved. The standard in the field is to perform eduction analyses and quantify anticipatory activity e.g. using the method of Harrisingh et al. (PMID: 18003827). This typically computed as the ratio of activity in the 3hrs preceding light transition to activity in the 6hrs preceding light transition. The programs referenced above should help with this.

Finally, in many cases I'm not sure that the appropriate statistical tests have been used e.g. in Fig 8c, 8e, 8h t-tests have been used when are three groups in the figure. The appropriate test here would an ANOVA, followed by post-hoc comparisons.

Reviewer #2 (Public Review):

The authors presented a well-written manuscript describing the comparison of active-learning methods with state-of-art methods for several datasets of pharmaceutical interest. This is a very important topic since active learning is similar to a cyclic drug design campaign such as testing compounds followed by designing new ones which could be used to further tests and a new design cycle and so on. The experimental design is comprehensive and adequate for proposed comparisons.

1) Text in figures still very small and difficult to read. Please redraw the figures increasing the font size: 10-12pt is ideal in comparison with the main text. In my opinion, it seems like the authors drew the Figure properly but there is an automatic resizing by inserting it in the document. Please consider ensuring that the font size will remain legible in the final document.

2) I think this work will benefit from a comparison of obtained models to other models reported in the literature and the interpretability of models (e.g. contribution of molecule groups to the modeled activity) as it is required by OECD guide for QSAR purposes.

Reviewer #3 (Public Review):

Zhang and Lauder characterized both aerobic and anaerobic metabolic energy contributions in schools and solitary fishes in the Giant danio (Devario aequipinnatus) over a wide range of water velocities. By using a highly sophisticated respirometer system, the authors measure the aerobic metabolisms by oxygen uptake rate and the non-aerobic oxygen cost as excess post-exercise oxygen consumption (EPOC). With these data, the authors model the bioenergetic cost of schools and solitary fishes. The authors found that fish schools have a J-shaped metabolism-speed curve, with reduced total energy expenditure per tail beat compared to solitary fish. Fish in schools also recovered from exercise faster than solitary fish. Finally, the authors conclude that these energetic savings may underlie the prevalence of coordinated group locomotion in fish.

The conclusions of this paper are mostly well supported by data.

Reviewer #2 (Public Review):

Summary:<br /> In this study, human induced pluripotent stem cell (hiPSC)-derived liver progenitor cell organoids were transplanted into four different transplantation sites in a mouse model of liver disease, using five organoid delivery methods. Organoids were transplanted into the vascularised chamber device established in the groin, or into the liver, spleen, and subcutaneous fat. Results show that the vascularised chamber had the highest organoid survival, 5.1x higher than the site with the second highest survival (p=0.0002), being the intra-hepatic scaffold approach. Animals with the vascularised chamber also had the highest human albumin levels (0.33 {plus minus} 0.09 ng/mL). No organoid survival was observed when delivered into the liver without a scaffold, or when injected into the spleen. Meager survival occurred in transplantations into subcutaneous fat.

Strengths:<br /> A systematic study with a clear line of experiments and well-presented results. The manuscript is well-written and easy to follow. The results and conclusions drawn are convincing.

Weaknesses:<br /> Although the number of organoids and albumin secretion is visibly higher in the vascularised chamber device, the organoids possess relatively higher Sox9+ cells compared to HNFa4a+ cells suggesting higher biliary differentiation than hepatic differentiation. On the other hand, although the intrahepatic scaffold approach, with a relatively smaller number of organoids and less albumin secretion, showed higher hepatic differentiation (although non-significant) suggesting that better scaffolds could be researched further to assess the clinical application of intrahepatic scaffold-based organoid transplantation.

Reviewer #2 (Public Review):

Summary:<br /> This paper described the role of BRCT repeat 5 in TOPBP1, a DNA damage response protein, in the maintenance of meiotic sex chromosome inactivation (MSCI). By analyzing a Topbp1 mutant mouse with amino acid substitutions in BRCT repeat 5, the authors found reduced phosphorylation of a DNA/RNA helicase, Sentaxin, and decreased localization of the protein to the X-Y sex body in pachynema. Moreover, the authors also found decreased repression of several genes on the sex chromosomes in the male mice.

Strengths:<br /> The works including phospho-proteomics and single-cell RNA sequencing with lots of data have been done with great care and most of the results are convincing.

Weaknesses:<br /> No weakness.

Reviewer #2 (Public Review):

In this manuscript, Libert et al. develop a model to predict an individual's age using physiological traits from multiple organ systems. The difference between the predicted biological age and the chronological age -- ∆Age, has an effect equivalent to that of a chronological year on Gompertz mortality risk. By conducting GWAS on ∆Age, the authors identify genetic factors that affect aging and distinguish those associated with age-related diseases. The study also uncovers environmental factors and employs dropout analysis to identify potential biomarkers and drivers for ∆Age. This research not only reveals new factors potentially affecting aging but also shows promise for evaluating therapeutics aimed at prolonging a healthy lifespan. This work represents a significant advancement in data-driven understanding of aging and provides new insights into human aging. Addressing the points raised would enhance its scientific validity and broaden its implications.

Major points:

1. Enhance the description and clarity of model evaluation.

The manuscript requires additional details regarding the model's evaluation. The authors have stated "To develop a model that predicts age, we experimented with several algorithms, including simple linear regression, Gradient Boosting Machine (GBM) and Partial Least Squares regression (PLS). The outcomes of these approaches were almost identical". It is currently unclear whether the 'almost identical outcomes' mentioned refer to the similarity in top contribution phenotypes, the accuracy of age prediction, or both. To resolve this ambiguity, it would be beneficial to include specific results and comparisons from each of these models.

Furthermore, the authors mention "to test for overfitting, a PLS model had been generated on randomly selected 90% of individuals and tested on the remaining 10% with similar results". To comprehensively assess the model's performance, it is crucial to provide detailed results for both the test and validation datasets. This should at least include metrics such as correlation coefficients and mean squared error for both training and test datasets.

2. External validation and generalization of results

To enhance the robustness and generalizability of the study's findings, it is crucial to perform external validation using an independent population. Specifically, conducting validation with the participants of the 'All of Us' research program offers a unique opportunity. This diverse and extensive cohort, distinct from the initial study group, will serve as an independent validation set, providing insights into the applicability of the study's conclusions across varied demographics.

Reviewer #2 (Public Review):

Kádková, Murach, Pedersen, and colleagues studied how three disease-causing missense mutations in SNAP25 affect synaptic vesicle exocytosis. These mutations have previously been studied by Alten et al., 2021. The authors observed similar impairments in spontaneous and evoked release as Alten et al., 2021, but the measurement of readily releasable pool (RRP) size differed between the two studies. The authors found that the V48F and D166Y mutations affect the interaction with the Ca2+ sensor synaptotagmin-1 (Syt1), but do not entirely phenocopy Syt1 loss-of-function because they also exhibit a gain-of-function. Thus, these mutations affect multiple aspects of the energy landscape for vesicle priming and fusion. The I67N mutation specifically increases the fusion energy barrier without affecting upstream vesicle priming.

The strength of the study includes careful and technically excellent dissection of the synaptic release process and a combination of electrophysiology, biophysics, and modeling approaches. This study gained a deeper mechanistic understanding of these mutations in vesicle exocytosis than the previous study but did not result in a paradigm shift in our understanding of SNAP25-associated encephalopathy because the same spontaneous and evoked release phenotypes were previously identified.

Comments on revised version:

The authors fully addressed the two previous technical concerns and improved the introduction of the paper.

Reviewer #2 (Public Review):

Summary:<br /> In the paper "Inter-regional Delays Fluctuate in the Human Cerebral Cortex," the authors aim to investigate how global changes in the power of brain oscillations affect the latency and strength of cortico-cortical couplings. They measured changes in brain oscillations and inter-regional couplings using human intracranial recordings. Additionally, the authors employed oscillator models to elucidate their empirical findings.

Strengths:<br /> The authors tested their hypotheses using human intracranial data, which provides a direct measurement of brain activity with high spatial and temporal resolution. This offers a unique insight into the interplay between oscillatory power and inter-regional coupling in the human brain.

Weaknesses:<br /> The authors had access to only a subset of brain regions. Although this limitation is common in many intracranial studies, their discussion of global changes in brain oscillations is impacted by the lack of whole-brain coverage, and thus the global nature of these oscillations should be interpreted with caution.

The description of the analysis procedure is not always clear.

Summary of main concerns:<br /> My primary concerns relate to possible circularity in the analysis and the incomplete reporting of statistical results. For instance, correlation values are often provided without associated p-values, making it difficult to assess their significance. Furthermore, in some sections of the text, it is unclear whether specific results are supported by any statistical tests.

Crucial information is buried in the supplemental materials (e.g., the figure showing results for broad-band high-frequency power). Some details about the specific paradigm are missing in the methods section, making it challenging to determine if additional controls are necessary in the analyses. I encourage the authors to clarify certain aspects of the analysis and results to ensure their conclusions are substantiated by the data. Should the results be robust, I believe the study will be significant for researchers interested in brain oscillations and beyond.

Reviewer #2 (Public Review):

Summary:<br /> Both human and non-human animals modulate the frequency of their vocalizations to communicate important information about context and internal state. While regulation of the size of the laryngeal opening is a well-established mechanism to regulate vocal pitch, the contribution of expiratory airflow to vocal pitch is less clear. To consider this question, this study first characterizes the relationship between the dominant frequency contours of adult mouse ultrasonic vocalizations (USVs) and expiratory airflow using whole-body plethysmography. Next, the authors build off of their previous work characterizing intermediate reticular oscillator (iRO) neurons in mouse pups to establish the existence of a genetically similar population of neurons in adults and show that artificial activation of iRO neurons elicits USV production in adults. Third, the authors examine the acoustic features of USV elicited by optogenetic activation of iRO and find that a majority of natural USV types (as defined by pitch contour) are elicited by iRO activation.

Strengths:<br /> Strengths of the study include the novel consideration of expiratory airflow as a mechanism to regulate vocal pitch and the use of intersectional methods to identify and activate the iRO in adult mice. The establishment of iRO neurons as a brainstem population that regulates vocal production across development is an important finding.

Weaknesses:<br /> The study does not include statistical analyses to compare the observed relationships between expiratory airflow and USV pitch to a null model in which expiratory airflow and USV pitch are unrelated. The findings of the study also do not provide clear evidence to support the authors' model in which distinct brainstem populations (iRO and RAm) independently regulate expiratory airflow and laryngeal adduction. Although this study establishes iRO as an important population that regulates USV production in adult mice, the question of whether and how different brainstem populations contribute differentially to vocal production remains an important open question. Lastly, the addition of statistical analyses would help to strengthen the study's conclusion that iRO activation positively biases the relationship between expiratory airflow and USV pitch across multiple USV types.

Reviewer #2 (Public Review):

Summary:<br /> Feng et al. investigated dynamic changes in functional and structural connectivity relationships across a broad age range from childhood to early adulthood (6-22 years) using the large open-source HCP-Development database of multimodal magnetic resonance imaging (MRI). Employing a multilinear model, the study integrates three white-matter structural descriptors derived from diffusion tractography with 'microstructure profile covariance' (MPC) descriptors of relationships between cortical regions in terms of regional T1w/T2w ratio, and evaluates the coupling between these structural connectome (SC) descriptors and functional connectivity (FC) as adjusted coefficients of determination, i.e. how well the structural descriptors correspond to the functional connectivity derived from resting-state functional MRI.

The findings reveal a global increase in SC-FC coupling over development. At a regional level, coupling exhibited distinct profiles of age-related increases and decreases within and between functional networks. Individual variability captured by the presented measures of SC-FC coupling was implicated as a potential marker for the prediction of general intelligence scores. Additionally, the investigation extended to associating changes in SC-FC coupling with age to regional gene expression profiles (derived from Allen Human Brain Atlas that analysed six neurotypical adult brains), suggesting positive associations with oligodendrocyte-related pathways and negative associations with astrocyte-related genes.

Strengths:<br /> Overall, the paper offers an interesting and valuable contribution to our understanding of structure-function relationships in the context of brain development. The commendable efforts to assess robustness across various methodologies, including brain parcellation and tractography, and reproducibility analyses on different data subsets enhance the paper's credibility. Combining cortical MPC with more usual white-matter descriptors of structural connectivity is interesting and provides (potentially) complementary information for the study of structure-function relationships with age. Analysing the changes in SC-FC coupling in relation to profiles of evolutionary expansion and functional principal gradients shows a good effort to position the present observations on SC-FC coupling within the previously described work.

Weaknesses:<br /> Although the paper has many strengths, some aspects of the analysis need to be clarified to further support the proposed conclusions. In particular:

* The authors propose that combining cortical and white-matter connectivity measures yields a more comprehensive descriptor of SC-FC coupling. While this is likely true, the claim is not directly tested by assessing different descriptors separately and then in combination to compare the benefits of incorporating additional information for the description of SC-FC coupling.

* The authors report changes in SC-FC coupling with myelin content (reporting a positive association of coupling with regional myelin) and report positive associations between SC-FC correlation with age and expression of oligodendrocyte-related genes. Given that the computation of SC-FC coupling involves the T1w/T2w ratios within cortical regions (recognised as a myelin marker), it's plausible that these findings may be influenced by potential bias introduced by myelin-related measures in the coupling computation process.

* The authors investigate the predictive power of SC-FC coupling, suggesting non-random (but weak) prediction of individual variability in general intelligence (after age correction). However, again the benefit of using SC-FC coupling measures over using each modality alone is not evaluated. Such comparison might indicate whether the coupling is an informative measure in itself or whether it might be informative only to the extent to which it is a proxy measure of SC and FC (in case the predictive power of each separate modality is much higher).

* Generally, more information on quality assessment of tractography and parcellations (including potential age effects on processing given the wide age range of the participants), additional details on the distribution of cognitive scores used in the predictive section, and further clarifications regarding the design choices and validation strategy would provide the reader with a more detailed understanding of the cohort and proposed analytical pipeline (these minor comments are included in the private recommendations to authors).

Reviewer #2 (Public Review):

Summary:<br /> This study offers a significant advancement in understanding liver innate lymphoid cell (ILC) biology by elucidating the role of the transcription factor Prdm1. It shows that Prdm1 is crucial in maintaining the balance between conventional natural killer (cNK) cells and ILC1s in the liver, with knockout models revealing a vital role in cancer defense mechanisms. Despite not affecting direct cytotoxicity, Prdm1 deficiency leads to increased cancer metastasis and reduced secretion of key molecules like IFN-γ, pointing to its importance in immune regulation. The use of single-cell RNA sequencing further underscores Prdm1's role in cellular communication within the liver's immune milieu. This study is a robust contribution to the field, providing insights that could inform new immunotherapy approaches for liver cancer.

Strengths:<br /> The study's strength lies in its comprehensive approach, combining the specificity of Prdm1 conditional deletion in Ncr1-cre mice with integrative omics analyses and cutting-edge cytometry to delineate Prdm1's role in liver Type 1 ILC biology and its functional implications in tumor immunity. This multifaceted strategy not only clarifies Prdm1's influence on ILC composition and maturation but also conveys potential therapeutic insights for liver cancer immunotherapy.

Weaknesses:<br /> A notable weakness of the study is the limited scope of in vivo disease models, primarily relying on the B16F10 melanoma model, which may not fully capture the complex behavior of Type 1 ILCs across diverse cancer types. Furthermore, the absence of direct human data, such as the effects of PRDM1 deletion in human NK cells or stem cells during their differentiation into NK and ILC1, leaves a gap in translating these findings to clinical settings.

Reviewer #2 (Public Review):

Summary:<br /> Two hypotheses could explain the observation that genes of more complex organisms tend to undergo more alternative splicing. On one hand, alternative splicing could be adaptive since it provides the functional diversity required for complexity. On the other hand, increased rates of alternative splicing could result through nonadaptive processes since more complex organisms tend to have smaller effective population sizes and are thus more prone to deleterious mutations resulting in more spurious splicing events (drift-barrier hypothesis). To evaluate the latter, B́enitiere et al. analyzed transcriptome sequencing data across 53 metazoan species. They show that proxies for effective population size and alternative splicing rates are negatively correlated. Furthermore, the authors find that rare, nonfunctional (and likely erroneous) isoforms occur more frequently in more complex species. Additionally, they show evidence that the strength of selection on splice sites increases with increasing effective population size and that the abundance of rare splice variants decreases with increased gene expression. All of these findings are consistent with the drift-barrier hypothesis.

This study conducts a comprehensive set of separate analyses that all converge on the same overall result and the manuscript is well organized. Furthermore, this study is useful in that it provides a modified null hypothesis that can be used for future tests of adaptive explanations for variation in alternative splicing.

Strengths:<br /> The major strength of this study lies in its complementary approach combining comparative and population genomics. Comparing evolutionary trends across phylogenetic diversity is a powerful way to test hypotheses about the origins of genome complexity. This approach alone reveals several convincing lines of evidence in support of the drift-barrier hypothesis. However, the authors also provide evidence from a population genetics perspective (using resequencing data for humans and fruit flies), making results even more convincing.

The authors are forward about the study's limitations and explain them in detail. They elaborate on possible confounding factors as well as the issues with data quality (e.g. proxies for Ne, inadequacies of short reads, heterogeneity in RNA-sequencing data).

Weaknesses:<br /> The authors primarily consider insects and mammals in their study. This only represents a small fraction of metazoan diversity. Sampling from a greater diversity of metazoan lineages would make these results and their relevance to broader metazoans substantially more convincing. Although the authors are careful about their tone, it is challenging to reconcile these results with trends across greater metazoans when the underlying dataset exhibits ascertainment bias and represents samples from only a few phylogenetic groups. Relatedly, some trends (such as Figure 1B-C) seem to be driven primarily by non-insect species, raising the question of whether some results may be primarily explained by specific phylogenetic groups (although the authors do correct for phylogeny in their statistics). How might results look if insects and mammals (or vertebrates) are considered independently?

Throughout the manuscript, the authors refer to infrequently spliced (mode <5%) introns as "minor introns" and frequently spliced (mode >95%) as "major introns". This is extremely confusing since "minor introns" typically represent introns spliced by the U12 spliceosome, whereas "major introns" are those spliced by the U2 spliceosome. Furthermore, it remains unclear whether the study only considers major introns or both major and minor introns. Minor introns typically have AT-AC splice sites whereas major introns usually have GT/GC-AG splice sites, although in rare cases the U2 can recognize AT-AC (see Wu and Krainer 1997 for example). The authors also note that some introns show noncanonical AT-AC splice sites while these are actually canonical splice sites for minor introns.

Reviewer #2 (Public Review):

The authors used various microscopy techniques, including super-resolution microscopy, to observe the changes that occur in the midpiece of mouse sperm flagella. Previously, it was shown that actin filaments form a double helix in the midpiece. This study reveals that the structure of these actin filaments changes after the acrosome reaction and before sperm-egg fusion, resulting in a thinner midpiece. Furthermore, by combining midpiece structure observation with calcium imaging, the authors show that changes in intracellular calcium concentrations precede structural changes in the midpiece. The cessation of sperm motility by these changes may be important for fusion with the egg. Elucidation of the structural changes in the midpiece could lead to a better understanding of fertilization and the etiology of male infertility. The conclusions of this manuscript are largely supported by the data, but there are several areas for improvement in data analysis and interpretation. Please see the major points below.

1. It is unclear whether an increased FM4-64 signal in the midpiece precedes the arrest of sperm motility. This needs to be clarified in order to argue that structural changes in the midpiece cause sperm motility arrest. The authors should analyze changes in both motility and FM4-64 signal over time for individual sperm.

2. It is possible that sperm stop moving because they die. Figure 1G shows that the FM4-64 signal is increased in the midpiece of immotile sperm, but it is necessary to show that the FM4-64 signal is increased in sperm that are not dead and retain plasma membrane integrity by checking sperm viability with propidium iodide or other means.

3. It is unclear how the structural change in the midpiece causes the entire sperm flagellum, including the principal piece, to stop moving. It will be easier for readers to understand if the authors discuss possible mechanisms.

4. The mitochondrial sheath and cell membrane are very close together when observed by transmission electron microscopy. The image in Figure 9A with the large space between the plasma membrane and mitochondria is misleading and should be corrected. The authors state that the distance between the plasma membrane and mitochondria approaches about 100 nm after the acrosome reaction (Line 330 - Line 333), but this is a very long distance and large structural changes may occur in the midpiece. Was there any change in the mitochondria themselves when they were observed with the DsRed2 signal?

5. In the TG sperm used, the green fluorescence of the acrosome disappears when sperm die. Figure 1C should be analyzed only with live sperm by checking viability with propidium iodide or other means.

Reviewer #2 (Public Review):

Summary:

This interesting paper examines the earliest steps in progesterone-induced frog oocyte maturation, an example of non-genomic steroid hormone signaling that has been studied for decades but is still very incompletely understood. In fish and frog oocytes it seems clear that mPR proteins are involved, but exactly how they relay signals is less clear. In human sperm, the lipid hydrolase ABHD2 has been identified as a receptor for progesterone, and so the authors here examine whether ABHD2 might contribute to progesterone-induced oocyte maturation as well. The main results are:

1. Knocking down ABHD2 makes oocytes less responsive to progesterone, and ectopically expressing ABHD2.S (but not the shorter ABHD2.L gene product) partially rescues responsiveness. The rescue depends upon the presence of critical residues in the protein's conserved lipid hydrolase domain, but not upon the presence of critical residues in its acyltransferase domain.

2. Treatment of oocytes with progesterone causes a decrease in sphingolipid and glycerophospholipid content within 5 min. This is accompanied by an increase in LPA content and arachidonic acid metabolites. These species may contribute to signaling through GPCRs. Perhaps surprisingly, there was no detectable increase in sphingosine-1-phosphate, which might have been expected given the apparent substantial hydrolysis of sphingolipids. The authors speculate that S1P is formed and contributes to signaling but diffuses away.

3. Pharmacological inhibitors of lipid-metabolizing enzymes support, for the most part, the inferences from the lipidomics studies, although there are some puzzling findings. The puzzling findings may be due to uncertainty about whether the inhibitors are working as advertised.

4. Pharmacological inhibitors of G-protein signaling support a role for G-proteins and GPCRs in this signaling, although again there are some puzzling findings.

5. Reticulocyte expression supports the idea that mPR and ABHD2 function together to generate a progesterone-regulated PLA2 activity.

6. Knocking down or inhibiting ABHD2 inhibited progesterone-induced mPRinternalization, and knocking down ABHD2 inhibited SNAP2520-induced maturation.

Strengths:

All in all, this could be a very interesting paper and a nice contribution. The data add a lot to our understanding of the process, and, given how ubiquitous mPR and AdipoQ receptor signaling appear to be, something like this may be happening in many other physiological contexts.

Weaknesses:

I have several suggestions for how to make the main points more convincing.

Main criticisms:

1. The ABHD2 knockdown and rescue, presented in Fig 1, is one of the most important findings. It can and should be presented in more detail to allow the reader to understand the experiments better. E.g.: the antisense oligos hybridize to both ABHD2.S and ABHD2.L, and they knock down both (ectopically expressed) proteins. Do they hybridize to either or both of the rescue constructs? If so, wouldn't you expect that both rescue constructs would rescue the phenotype since they both should sequester the AS oligo? Maybe I'm missing something here.

In addition, it is critical to know whether the partial rescue (Fig 1E, I, and K) is accomplished by expressing reasonable levels of the ABHD2 protein, or only by greatly overexpressing the protein. The author's antibodies do not appear to be sensitive enough to detect the endogenous levels of ABHD2.S or .L, but they do detect the overexpressed proteins (Fig 1D). The authors could thus start by microinjecting enough of the rescue mRNAs to get detectable protein levels, and then titer down, assessing how low one can go and still get rescue. And/or compare the mRNA levels achieved with the rescue construct to the endogenous mRNAs.

Finally, please make it clear what is meant by n = 7 or n = 3 for these experiments. Does n = 7 mean 7 independently lysed oocytes from the same frog? Or 7 groups of, say, 10 oocytes from the same frog? Or different frogs on different days? I could not tell from the figure legends, the methods, or the supplementary methods. Ideally one wants to be sure that the knockdown and rescue can be demonstrated in different batches of oocytes, and that the experimental variability is substantially smaller than the effect size.

2. The lipidomics results should be presented more clearly. First, please drop the heat map presentations (Fig 2A-C) and instead show individual time course results, like those shown in Fig 2E, which make it easy to see the magnitude of the change and the experiment-to-experiment variability. As it stands, the lipidomics data really cannot be critically assessed.

[Even as heat map data go, panels A-C are hard to understand. The labels are too small, especially on the heat map on the right side of panel B. The 25 rows in panel C are not defined (the legend makes me think the panel is data from 10 individual oocytes, so are the 25 rows 25 metabolites? If so, are the individual oocyte data being collapsed into an average? Doesn't that defeat the purpose of assessing individual oocytes?) And those readers with red-green colorblindness (8% of men) will not be able to tell an increase from a decrease. But please don't bother improving the heat maps; they should just be replaced with more informative bar graphs or scatter plots.]

3. The reticulocyte lysate co-expression data are quite important and are both intriguing and puzzling. My impression had been that to express functional membrane proteins, one needed to add some membrane source, like microsomes, to the standard kits. Yet it seems like co-expression of mPR and ABHD2 proteins in a standard kit is sufficient to yield progesterone-regulated PLA2 activity. I could be wrong here - I'm not a protein expression expert - but I was surprised by this result, and I think it is critical that the authors make absolutely certain that it is correct. Do you get much greater activities if microsomes are added? Are the specific activities of the putative mPR-ABHD2 complexes reasonable?

Reviewer #2 (Public Review):

Summary:

Here, Yue et al. set out to determine if the low DNMT3B expression that is observed prior to de novo DNA methylation (before the blastocyst stage) has a function. Re-analyzing existing DNA methylation data from Smith et al. (2012) they find a small DNA methylation gain over a subset of promoters and gene bodies, occurring between the 8-cell and blastocyst stages, and refer to this as "minor de novo DNA methylation". They attempt to assess the relevance/functionality of this minor DNA methylation gain, and report reduced H3K27me3 in Dnmt3b knockdown (KD) trophoblast cells that normally undergo imprinted X-chromosome inactivation (iXCI) before the blastocyst stage. In addition, they assess the proliferation, differentiation, metabolic function, implantation rate, and live birth rate of Dnmt3b KD blastocysts.

Strengths:

Working with early embryos is technically demanding, making the well-designed experiments from this manuscript useful to the epigenetics community. Particularly, the DNMT3B expression and 5-mC staining at different embryonic stages.

Weaknesses:

- Throughout the manuscript, please represent DNA methylation changes as delta DNA methylation instead of fold change.

- Detailed methods on the re-analysis of the DNA methylation data from Smith et al. 2012 are missing from the materials and methods section. Was a minimum coverage threshold used?

- Detailed methods on the establishment and validation of Dnmt3b KO blastocysts and 5-aza-dC treated blastocysts are missing (related to Figure 2).

- Detailed methods on the re-analysis of the ChIPseq data from Liu et al. 2016 are missing from the materials and methods section.

- Some of the data represented in bar graphs does not look convincing/significant. Maybe this data can be better represented differently, such as in box plots or violin plots, which would better represent the data.

- The relevance and rationale for experiments using 5-aza-dC treatment is unclear.

Reviewer #2 (Public Review):

Summary:

The manuscript entitled "N-cadherin mechanosensing in ovarian follicles controls oocyte maturation and ovulation" aimed to investigate the role of N-cadherin in different ovarian physiological processes, including cumulus oocyte expansion, oocyte maturation, and ovulation. The authors performed several in vitro and in vivo mice experiments, using diverse techniques to reinforce their results.

First, they identified two compounds (N-cadherin antagonists) that block the adhesion of periovulatory COCs to fibronectin through screening a small molecule library, using the xCELLigenceTM system, performing proper and complementary controls. Second, the authors showed the presence of N-cadherin adherens junctions between granulosa cells and cumulus cells and at the interface of cumulus cell transzonal projections and the oocyte throughout folliculogenesis. And that these adherens complexes between cumulus cells and oocytes were disrupted when inhibited N-cadherin, as observed by nice representative confocal images. Then, the authors assessed COC expansion and oocyte meiotic maturation to determine whether the loss of oocyte membrane β-catenin and E-cadherin upon N-cadherin inhibitor treatment disrupts the bi-directional communication between cumulus cells and the oocyte. Indeed, N-cadherin antagonists disrupted both processes (cumulus expansion and oocyte meiotic). However, the expression of known mediators of COC expansion (E.g., Areg and Ptgs2) were either increased or unaffected. Nevertheless, RNA-Seq showed consistent effects on cell signaling mRNA genes by the antagonist CRS-066.

In vivo studies using mice were also achieved using stimulated protocols (together with one of the antagonists or vehicle) or granulosa-specific Cdh2 Knockouts to further analyze the role of N-cadherin. N-cadherin antagonist CRS-066 (but not LCRF-0006) significantly reduced mouse ovulation compared to controls. RNA-sequencing data analysis identified distinct gene expression profiles in CRS-066 treated compared to control ovaries. Ovulation in CdhFl/FL; Amhr2Cre mice after stimulation were also significantly reduced; multiple large unruptured follicles were observed in these granulosa-specific Cdh2 mutant ovaries, and the mRNA expression of Areg and Ptgs2 were reduced.

The authors conclude that their study identified N-cadherin as a mechanosensory regulator important in ovarian granulosa cell differentiation able to respond to hormone stimuli both in vivo and in vitro, demonstrating a critical role for N-cadherin in ovarian follicular development and ovulation. They highlighted the potential to inhibit ovulation by targeting this signaling mechanism.

Strengths:<br /> This remarkable manuscript is very well designed, performed, and discussed. The authors analyzed different aspects, and their data supports their conclusions.

Weaknesses:<br /> This study was performed using the mouse as a research model; further studies in larger animals and humans would be interesting and warranted.

Minor comments:

Some results are intriguing. While the AREG y PTGS2 mRNA increased within the COC in vitro by the N-cadherin antagonists, in vivo, the treatment induced a significant increase in both genes when analyzing the whole ovary. What are the authors' ideas that could explain these discrepancies in outcomes?

The authors stated that the ovaries from mice treated in the same manner and collected either before hCG treatment (eCG 44 h) or 11 h after hCG showed equivalent numbers of follicles at each stage of development from primary to antral. However, in Panel l from Figure 5, there is a significant increase in the number of antral follicles in the CRS-066 group (hCG 11 h) compared to the vehicle. Could the author discuss it in the manuscript?

Reviewer #2 (Public Review):

In this study Weinberger et al. investigated cardiac macrophage subsets after ischemia/reperfusion (I/R) injury in mice. The authors studied a ∆FIRE mouse model (deletion of a regulatory element in the Csf1r locus), in which only tissue resident macrophages might be ablated. The authors showed a reduction of resident macrophages in ∆FIRE mice and characterized its macrophages populations via scRNAseq at baseline conditions and after I/R injury. 2 days after I/R protocol ∆FIRE mice showed an enhanced pro inflammatory phenotype in the RNAseq data and differential effects on echocardiographic function 6 and 30 days after I/R injury. Via flow cytometry and histology the authors confirmed existing evidence of increased bone marrow-derived macrophage infiltration to the heart, specifically to the ischemic myocardium. Macrophage population in ∆FIRE mice after I/R injury were only changed in the remote zone. Further RNAseq data on resident or recruited macrophages showed transcriptional differences between both cell types in terms of homeostasis-related genes and inflammation. Depleting all macrophage using a Csf1r inhibitor resulted in a reduced cardiac function and increased fibrosis.

Reviewer #3 (Public Review):

Summary:

In the manuscript by Xiong and colleagues, the roles of TLR2 in hair follicle cycle regulation were investigated. By analyzing published dataset and using immunostaining and transgenic TLR2-GFP reporter mice, the authors showed that TLR2 expression is increased in the late telogen compared to the early telogen, implying that it is important for the transition between telogen to anagen hair cycle. They found that the genetic deletion of Tlr2 in hair follicle stem cells delays hair cycle entry in both homeostatic and wound-induced hair follicle regeneration. In addition, they found that CEP is an endogenous TLR2 activating ligand and triggers the progression of hair cycle in a TLR2-dependent manner. Mechanistically, the activation of TLR2 signaling antagonizes BMP signaling which is critical for the maintenance of hair follicle stem cell quiescence. Clinically, they showed that TLR2 expression is decreased in aging and high-fat diet condition, suggesting that the dysfunctional regulation of TLR2 pathway is responsible for age-related and obesity-related hair thinning and hair loss phenotypes.

Strengths:

Overall, this study presents the role and mechanism of TLR2 in regulating hair follicle regeneration. The functional interrogation parts using HFSC-specific TLR2 genetic deletion is solid, and an endogenous regulator, CEP, is identified.

Weaknesses:<br /> 1)<br /> - In SFig1A, the IF staining of TLR2 and Tlr2-GFP expression seem almost 100% co-localized, which is not usual experimentally.<br /> - In Fig 2J, the relative expression levels of Tlr2 in anagen, telogen, catagen HFSCs were tested. But it is just relative comparison and does not mean whether the expression level is meaningful or not. To make this convincing, adding other cell types such as dermal fibroblasts and immunes to the comparison as negative and positive controls would be a good idea.<br /> - In Fig 2K, the expression of Tlr2 is comparable or a bit lesser in epidermal cells and HFSCs, but the expressions of TLR2 (IF) and Tlr2-GFP in epidermal cells have not been presented at all in the manuscript. As the authors used K15-CrePR1 mice to delete Tlr2 in HFSCs specifically, showing TLR2 IF staining in TLR2-HFSC-KO mice would be nice evidence of significant expression of TLR2 in HFSCs. (still TLR2 expression in epidermis, but no TLR2 expression in HFSCs).<br /> - In Fig 1B, it is still unclear whether TLR2 staining is in epithelial cell or in dermal cells. TLR2 staining patterns in Fig 1B, SFig 1A, and rebuttal seem different. In Fig S1B and rebuttal, TLR2 expression in HFSCs, HG, DP cells, but in Fig 1B, most of HG and DP cells are not TLR2+.<br /> - Together, this reviewer still does not think that there is a clear and solid evidence of Tlr2 expression in HFSCs. Searching the Tlr2 expression in published bulk and single cell RNA-seq dataset would be helpful.

2)<br /> - In SFig 4B, C, the activation of BMP signaling was hindered by TLR2 signaling activation by PAM3CSK4. But it is in vitro data, and cultured HFSCs are different from in vivo HFSCs, and particularly the changes of HFSCs from quiescence to activation can hardly be recapitulated in vitro.<br /> - In Fig 4H, it is curious that in TLR2-HFSC-KO mice, P21 HFSCs showed no pSMAD1/5/9, but it is increased in P24.<br /> - Also, it is wondered that if ID1 and ID2, key target genes, are increased in TLR2-HFCS-KO.<br /> - The author suggested that BMP7 is a key connection between TLR2 signaling and BMP signaling. It is curious whether BMP7 is a direct target of TLR2 pathway? Are there Nfkb (putative) binding sites in cis-regulatory regions of BMP7?

3)<br /> - In Fig 6C, CEP expression is close to hair follicle in both anagen and telogen. Also, in Telogen, CEP expression is strong and very close to HFSCs. But In rebuttal Fig 2, CEP is localized to sebaceous gland, where MPO, a CEP producing enzyme, is expressed. Which one is correct? Also, if CEP is strongly expressed in Telogen (Fig 6C), how can HFSCs stay in quiescence with decreased BMP signaling?

Reviewer #2 (Public Review):

Summary:

Nonalcoholic fatty liver disease (NASH), recently renamed as metabolic dysfunction-associated steatohepatitis (MASH) is a leading cause of liver-related death. Farnesoid X receptor (FXR) is a promising drug target for treating NASH and several drugs targeting FXR are under clinical investigation for their efficacy in treating NASH. The authors intended to address whether FXR mediates its hepatic protective effects through the regulation of lncRNAs, which would provide novel insights into the pharmacological targeting of FXR for NASH treatment. The authors went from an unbiased transcriptomics profiling to identify a novel enhancer-derived lncRNA FincoR enriched in the liver and showed that the knockdown of FincoR in a murine NASH model attenuated part of the effect of tropifexor, an FXR agonist, namely inflammation and fibrosis, but not steatosis. This study provides a framework for how one can investigate the role of noncoding genes in pharmacological intervention targeting known protein-coding genes. Given that many disease-associated genetic variants are located in the non-coding regions, this study, together with others, may provide useful information for improved and individualized treatment for metabolic disorders.

Strengths:

The study leverages both transcriptional profile and epigenetic signatures to identify the top candidate eRNA for further study. The subsequent biochemical characterization of FincoR using FXR-KO mice combined with Gro-seq and Luciferase reporter assays convincingly demonstrates this eRNA as a FXR transcriptional target sensitive to FXR agonists. The use of in vitro culture cells and the in vivo mouse model of NASH provide multi-level evaluation of the context-dependent importance of the FincoR downstream of FXR in the regulation of functions related to liver dysfunction.

Weaknesses:

As discussed, future work to dissect the mechanisms by which FincoR facilitates the action of FXR and its agonists is warranted. It would be helpful if the authors could base this on the current understanding of eRNA modes of action and the observed biochemical features of FincoR to speculate potential molecular mechanisms explaining the observed functional phenotype. It is unclear if this eRNA is conserved in humans in any way, which will provide relevance to human disease. Additionally, the eRNA knockdown was achieved by deletion of an upstream region of the eRNA transcription. A more direct approach to alter eRNA levels, e.g., overexpression of FincoR in the liver would provide important data to interpret its functional regulation.

Reviewer #2 (Public Review):

Summary:

The objective of authors using metabolomics analysis of primary angle closure glaucoma (PACG) is to demonstrate that serum androstenedione is a novel biomarker that can be used to diagnose PACG and predict visual field progression.

Strengths:

Use of widely targeted and untargeted metabolite detection conditions. Use of liquid chromatography-tandem mass spectrometry and a chemiluminescence method for confirmation of androstenedione.

The authors have incorporated the relevant changes in their manuscript and improved the presentation.

Reviewer #2 (Public Review):

The exact dynamics of responses to volatiles from herbivore-attacked neighbouring plants have been little studied so far. Also, we still lack evidence whether herbivore-induced plant volatiles (HIPVs) induce or prime plant defences of neighbours. The authors investigated the volatile emission patterns of receiver plants that respond to the volatile emission of neighbouring sender plants which are fed upon by herbivorous caterpillars. They applied a very elegant approach (more rigorous than the current state-of-the-art) to monitor temporal response patterns of neighbouring plants to HIPVs by measuring volatile emissions of senders and receivers, senders only and receivers only. Different terpenoids were produced within 2 h of such exposure in receiver plants, but not during the dark phase. Once the light turned on again, large amounts of terpenoids were released from the receiver plants. This may indicate a delayed terpene burst, but terpenoids may also be induced by the sudden change in light. As one contrasting control, the authors also studied the time-delay in volatile emission when plants were just kept under continuous light. Here they also found a delayed terpenoid production, but this seemed to be lower compared to the plants exposed to the day-night-cycle. Another helpful control was now performed for the revision in which the herbivory treatment was started in the evening hours and lights were left on. This experiment revealed that the burst of terpenoid emission indeed shifted somewhat. Circadiane and diurnal processes must thus interact.

Interestingly, internal terpene pools of one of the leaves tested here remained more comparable between night and day, indicating that their pools stay higher in plants exposed to HIPVs. In contrast, terpene synthases were only induced during the light-phase, not in the dark-phase. Moreover, jasmonates were only significantly induced 22 h after onset of the volatile exposure and thus parallel with the burst of terpene release.

An additional experiment exposing plants to the green leaf volatile (glv) (Z)-3-hexenyl acetate revealed that plants can be primed by this glv, leading to a stronger terpene burst. The results are discussed with nice logic and considering potential ecological consequences. All data are now well discussed.

Overall, this study provides intriguing insights in the potential interplay between priming and induction, which may co-occur, enhancing (indirect and direct) plant defence. Follow-up studies are suggested that may provide additional evidence.

Reviewer #2 (Public Review):

The manuscript by Harry and Zakas determined the extent to which gene expression differences contribute to developmental divergence by using a model that has two distinct developmental morphs within a single species. Although the authors did collect a valuable dataset and trends in differential expression between the two morphs of S. benedicti were presented, we found limitations about the methods, system, and resources that the authors should address.

We have two major points:

1. Background information about the biological system needs to be clarified in the introduction of this manuscript. The authors stated that F1 offspring can have intermediate larval traits compared to the parents (Line 81). However, the authors collected F1 offspring at the same time as the mother in the cross. If offspring have intermediate larval traits, their developmental timeline might be different than both parents and necessitate the collection of offspring at different times to obtain the same stages as the parents. Could the authors (1) explain why they collected offspring at the same time as parents given that other literature and Line 81 state these F1 offspring develop at intermediate rates, and (2) add the F1 offspring to Figure 1 to show morphological and timeline differences in development?

Additionally, the authors state (Lines 83-85) that they detail the full-time course of embryogenesis for both the parents and the F1 crosses. However, we do not see where the authors have reported the full-time course for embryogenesis of the F1 offspring. Providing this information would shape the remaining results of the manuscript.

2. We have several concerns about the S. benedicti genome and steps regarding the read mapping for RNA-seq:

The S. benedicti genome used (Zakas et al. 2022) was generated using the PP morph. The largest scaffolds of this assembly correspond to linkage groups, showing the quality of this genome. The authors should point out in the Methods and/or Results sections that the quality of this genome means that PP-specific gene expression can be quantified well. However, the challenges and limitations of mapping LL-specific expression data to the PP genome should be discussed.

It is possible that the authors did not find exclusive gene expression in the LL morph because they require at least one gene to be turned on in one morph as part of the data-cleaning criteria. Because the authors are comparing all genes to the PP morph, they could be missing true exclusive genes responsible for the biological differences between the two morphs. Did they make the decision to only count genes expressed in one stage of the other morph because the gene models and mapping quality led to too much noise?

The authors state that the mapping rates between the two morphs are comparable (Supplementary Figure 1). However, there is a lot of variation in mapping the LL individuals (~20% to 43%) compared to the PP individuals. What is the level of differentiation within the two morphs in the species (pi and theta)? The statistical tests for this comparison should be added and the associated p-value should be reported. The statistical test used to compare mapping rates between the two morphs may be inappropriate. The authors used Salmon for their RNA alignment and differential expression analysis, but it is possible that a different method would be more appropriate. For example, Salmon has some limitations as compared to Kallisto as others have noted. The chosen statistical test should be explained, as well as how RNA-seq data are processed and interpreted.

What about the read mapping rate and details for the F1 LP and PL individuals? How did the offspring map to the P genome? These details should be included in Supplementary Figure 1. Could the authors also provide information about the number of genes expressed at each stage in the F1 LP and PL samples in S Figure 2? How many genes went into the PCA? Many of these details are necessary to evaluate the F1 RNA-seq analyses.

Generally, the authors need to report the statistics used in data processing more thoroughly. The authors need to report the statistics used to (1) process and evaluate the RNA-seq data and (2) determine the significance between the two morphs (Supplementary Figures 1 and 2).

Reviewer #2 (Public Review):

Summary:

Zhou et al report development of a new method, Rec-Seq, that allows rigorous quantitation of the efficiency of 48S ribosomal pre-initiation complex (PIC) formation on messenger RNAs at transcriptome scale in vitro. With a next-generation deep-sequencing approach, Rec-Seq allows precisely targeted dissection of the roles of translation initiation factors in PIC assembly. This level of molecular precision is important to understanding mechanisms of translational control, making Rec-Seq a significant methodological advance. The authors leverage Rec-Seq to investigate the relative roles of two key helicase enzymes, Ded1p and eIF4A. While past work has pointed to differing roles for Ded1p and eIF4A helicase activity in PIC assembly, unambiguous interpretation of prior in-vivo data has been hindered by technical requirements for performing the experiments in cells. Rec-Seq circumvents these challenges, providing robust mechanistic insights. The authors find that Ded1p stimulates PIC formation selectively on mRNAs with long, structured leaders in the Rec-Seq system, while eIF4A provides much more general stimulation across mRNAs. The findings substantiate the past in-vivo results, along with adding new insights. They contrast with evidence that Ded1p promotes translation by suppressing inhibitory upstream initiation through structural remodeling, or through formation of intracellular, phase-separated granules. The conclusions of the study are generally well-supported by the data.

Strengths:

The quantitative nature of Rec-Seq, which uses an internal standard to measure absolute recruitment efficiencies, is an important strength.

The methodology decisively overcomes past experimental limitations, allowing the authors to make clear conclusions with regard to the relative roles of Ded1p and eIF4A in PIC formation. An important and useful addition to the toolbox for studying translation and translational control mechanisms, Rec-Seq substantially expands the throughput and scope of mechanistic analyses for translation initiation.

One significant finding to emerge is that the in-vitro reconstituted system used here recapitulates effects of in-vivo perturbations of translation initiation. Despite the lack of a cellular environment and its components, PIC formation appears to operate much as it does in the cell. Importantly, this highlights an inherent "modularity" to the system that is especially of interest in the context of how regulatory machinery beyond the PIC may control translation.

Weaknesses:

Several findings in this report are quite surprising and may require additional work to fully interpret. Primary among these is the finding that Ded1p stimulates accumulation of PICs at internal site in mRNA coding sequences at an incidence of up to ~50%. The physiological relevance of this is unclear.

A limitation of the methodology is that, as an endpoint assay, Rec-Seq does not readily decouple effects of Ded1p on PIC-mRNA loading from those on the subsequent scanning step where the PIC locates the start codon. Considering that Ded1p activity may influence each of these initiation steps through distinct mechanisms - i.e., binding to the mRNA cap-recognition factor eIF4F, or direct mRNA interaction outside eIF4F - additional studies may be needed to gain deeper mechanistic insights.

As the authors note, the achievable Ded1p concentrations in Rec-Seq may mask potential effects of Ded1p-based granule formation on translation initiation. Additional factors present in the cell could potentially also promote this mechanism. Consequently, the results do not fully rule out granule formation as a potential parallel Ded1p-mediated translation-inhibitory mechanism in cells.

Reviewer #3 (Public Review):

There has been a long-standing link between the biology of sulfur-containing molecules (e.g., hydrogen sulfide gas, the amino acid cysteine, and its close relative cystine, et cetera) and the biology of hypoxia, yet we have a poor understanding of how and why these two biological processes and are co-regulated. Here, the authors use C. elegans to explore the relationship between sulfur metabolism and hypoxia, examining the regulation of cysteine dioxygenase (CDO1 in humans, CDO-1 in C. elegans), which is critical to cysteine catabolism, by the hypoxia inducible factor (HIF1 alpha in humans, HIF-1 in C. elegans), which is the key terminal effector of the hypoxia response pathway that maintains oxygen homeostasis. The authors are trying to demonstrate that (1) the hypoxia response pathway is a key regulator of cysteine homeostasis, specifically through the regulation of cysteine dioxygenase, and (2) that the pathway responds to changes in cysteine homeostasis in a mechanistically distinct way from how it responds to hypoxic stress.

Briefly summarized here, the authors initiated this study by generating transgenic animals expressing a CDO-1::GFP protein chimera from the cdo-1 promoter so that they could identify regulators of CDO-1 expression through a forward genetic screen. This screen identified mutants with elevated CDO-1::GFP expression in two genes, egl-9 and rhy-1, whose wild-type products are negative regulators of HIF-1, raising the possibility that cdo-1 is a HIF-1 transcriptional target. Indeed, the authors provide data showing that cdo-1 regulation by EGL-9 and RHY-1 is dependent on HIF-1 and that regulation by RHY-1 is dependent on CYSL-1, as expected from other published findings of this pathway. The authors show that exogenous cysteine activates cdo-1 expression, reflective of what is known to occur in other systems. Moreover, they find that exogenous cysteine is toxic to worms lacking CYSL-1 or HIF-1 activity, but not CDO-1 activity, suggesting that HIF-1 mediates a survival response to toxic levels of cysteine and that this response requires more than just the regulation of CDO-1. The authors validate their expression studies using a GFP knockin at the cdo-1 locus, and they demonstrate that a key site of action for CDO-1 is the hypodermis. They present genetic epistasis analysis supporting a role for RHY-1, both as a regulator of HIF-1 and as a transcriptional target of HIF-1, in offsetting toxicity from aberrant sulfur metabolism. The authors use CRISPR/Cas9 editing to mutate a key amino acid in the prolyl hydroxylase domain of EGL-9, arguing that EGL-9 inhibits CDO-1 expression through a mechanism that is largely independent of the prolyl hydroxylase activity.

Overall, the data seem rigorous, and the conclusions drawn from the data seem appropriate. The experiments test the hypothesis using logical and clever molecular genetic tools and design. The sample size is a bit lower than is typical for C. elegans papers; however, the experiments are clearly not underpowered, so this is not an issue. The paper is likely to drive many in the field (including the authors themselves) into deeper experiments on (1) how the pathway senses hypoxia and sulfur/cysteine/H2S using these distinct mechanisms/modalities, (2) how oxygen and sulfur/cysteine/H2S homeostasis influence one another, and (3) how this single pathway evolved to sense and respond to both of these stress modalities.

My previous concerns have been addressed. The authors are commended on an excellent body of research.

Reviewer #2 (Public Review):

Summary:<br /> This paper aimed to examine the spatial frequency selectivity of macaque inferotemporal (IT) neurons and its relation to category selectivity. The authors suggest in the present study that some IT neurons show a sensitivity for the spatial frequency of scrambled images. Their report suggests a shift in preferred spatial frequency during the response, from low to high spatial frequencies. This agrees with a coarse-to-fine processing strategy, which is in line with multiple studies in the early visual cortex. In addition, they report that the selectivity for faces and objects, relative to scrambled stimuli, depends on the spatial frequency tuning of the neurons.

Strengths:<br /> Previous studies using human fMRI and psychophysics studied the contribution of different spatial frequency bands to object recognition, but as pointed out by the authors little is known about the spatial frequency selectivity of single IT neurons. This study addresses this gap and they show that at least some IT neurons show a sensitivity for spatial frequency and interestingly show a tendency for coarse-to-fine processing.

Weaknesses and requested clarifications:<br /> 1. It is unclear whether the effects described in this paper reflect a sensitivity to spatial frequency, i.e. in cycles/ deg (depends on the distance from the observer and changes when rescaling the image), or is a sensitivity to cycles /image, largely independent of image scale. How is it related to the well-documented size tolerance of IT neuron selectivity?

2. The authors' band-pass filtered phase scrambled images of faces and objects. The original images likely differed in their spatial frequency amplitude spectrum and thus it is unclear whether the differing bands contained the same power for the different scrambled images. If not, this could have contributed to the frequency sensitivity of the neurons.

3. How strong were the responses to the phase-scrambled images? Phase-scrambled images are expected to be rather ineffective stimuli for IT neurons. How can one extrapolate the effect of the spatial frequency band observed for ineffective stimuli to that for more effective stimuli, like objects or (for some neurons) faces? A distribution should be provided, of the net responses (in spikes/s) to the scrambled stimuli, and this for the early and late windows.

4. The strength of the spatial frequency selectivity is unclear from the presented data. The authors provide the result of a classification analysis, but this is in normalized units so that the reader does not know the classification score in percent correct. Unnormalized data should be provided. Also, it would be informative to provide a summary plot of the spatial frequency selectivity in spikes/s, e.g. by ranking the spatial frequency bands for each neuron based on half of the trials and then plotting the average responses for the obtained ranks for the other half of the trials. Thus, the reader can appreciate the strength of the spatial frequency selectivity, considering trial-to-trial variability. Also, a plot should be provided of the mean response to the stimuli for the two analysis windows of Figure 2c and 2d in spikes/s so one can appreciate the mean response strengths and effect size (see above).

5. It is unclear why such brief stimulus durations were employed. Will the results be similar, in particular the preference for low spatial frequencies, for longer stimulus durations that are more similar to those encountered during natural vision?

6. The authors report that the spatial frequency band classification accuracy for the population of neurons is not much higher than that of the best neuron (line 151). How does this relate to the SNC analysis, which appears to suggest that many neurons contribute to the spatial frequency selectivity of the population in a non-redundant fashion? Also, the outcome of the analyses should be provided (such as SNC and decoding (e.g. Figure 1D)) in the original units instead of undefined arbitrary units.

7. To me, the results of the analyses of Figure 3c,d, and Figure 4 appear to disagree. The latter figure shows no correlation between category and spatial frequency classification accuracies while Figure 3c,d shows the opposite.

8. If I understand correctly, the "main" test included scrambled versions of each of the "responsive" images selected based on the preceding test. Each stimulus was presented 15 times (once in each of the 15 blocks). The LDA classifier was trained to predict the 5 spatial frequency band labels and they used 70% of the trials to train the classifier. Were the trained and tested trials stratified with respect to the different scrambled images? Also, LDA assumes a normal distribution. Was this the case, especially because of the mixture of repetitions of the same scrambled stimulus and different scrambled stimuli?

9. The LDA classifiers for spatial frequency band (5 labels) and category (2 labels) have different chance and performance levels. Was this taken into account when comparing the SNC between these two classifiers? Details and SNC values should be provided in the original (percent difference) instead of arbitrary units in Figure 5a. Without such details, the results are impossible to evaluate.

10. Recording locations should be described in IT, since the latter is a large region. Did their recordings include the STS? A/P and M/L coordinate ranges of recorded neurons?

11. The authors should show in Supplementary Figures the main data for each of the two animals, to ensure the reader that both monkeys showed similar trends.

12. The authors found that the deep nets encoded better the spatial frequency bands than the IT units. However, IT units have trial-to-trial response variability and CNN units do not. Did they consider this when comparing IT and CNN classification performance? Also, the number of features differs between IT and CNN units. To me, comparing IT and CNN classification performances is like comparing apples and oranges.

13. The authors should define the separability index in their paper. Since it is the main index to show a relationship between category and spatial frequency tuning, it should be described in detail. Also, results should be provided in the original units instead of undefined arbitrary units. The tuning profiles in Figure 3A should be in spikes/s. Also, it was unclear to me whether the classification of the neurons into the different tuning profiles was based on an ANOVA assessing per neuron whether the effect of the spatial frequency band was significant (as should be done).

14. As mentioned above, the separability analysis is the main one suggesting an association between category and spatial frequency tuning. However, they compute the separability of each category with respect to the scrambled images. Since faces are a rather homogeneous category I expect that IT neurons have on average a higher separability index for faces than for the more heterogeneous category of objects, at least for neurons responsive to faces and/or objects. The higher separability for faces of the two low- and high-pass spatial frequency neurons could reflect stronger overall responses for these two classes of neurons. Was this the case? This is a critical analysis since it is essential to assess whether it is category versus responsiveness that is associated with the spatial frequency tuning. Also, I do not believe that one can make a strong claim about category selectivity when only 6 faces and 3 objects (and 6 other, variable stimuli; 15 stimuli in total) are employed to assess the responses for these categories (see next main comment). This and the above control analysis can affect the main conclusion and title of the paper.

15. For the category decoding, the authors employed intact, unscrambled stimuli. Were these from the main test? If yes, then I am concerned that this represents a too small number of stimuli to assess category selectivity. Only 9 fixed + 6 variable stimuli = 15 were in the main test. How many faces/ objects on average? Was the number of stimuli per category equated for the classification? When possible use the data of the preceding selectivity test which has many more stimuli to compute the category selectivity.

Reviewer #2 (Public Review):

Summary:<br /> In the manuscript by Oestreicher et al, the authors use patch-clamp electrophysiology, immunofluorescent imaging of the cochlea, auditory function tests, and single-unit recordings of auditory afferent neurons to probe the unique properties of calcium signaling in cochlear hair cells that allow rapid and sustained neurotransmitter release. The calcium-binding proteins (CaBPs) are thought to modify the inactivation of the Cav1.3 calcium channels in IHCs that initiate vesicle fusion, reducing the calcium-dependent inactivation (CDI) of the channels to allow sustained calcium influx to support neurotransmitter release. The authors use knockout mice of Cabp1 and Cabp2 in a double knockout (Cabp1/2 DKO) to show that these molecules are required for enabling sustained calcium currents by reducing CDI and enabling proper IHC neurotransmitter release. They further support their evidence by re-introducing Cabp2 using an injection of AAV containing the Cabp2 sequence into the cochlea, which restores some of the auditory function and reduces CDI in patch-clamp recordings.

Strengths:<br /> Overall the data is convincing that Cabp1/2 is required for reducing CDI in cochlear hair cells, allowing their sustained neurotransmitter release and sound encoding. Figures are well-prepared, recordings are careful and stats are appropriate, and the manuscript is well-written. The discussion appropriately considers aspects of the data that are not yet explained and await further experimentation.

Weaknesses:<br /> There are some sections of the manuscript that pool data from different experiments with slightly different conditions (wt data from a previous paper, different calcium concentrations, different holding voltages, tones vs clicks, etc). This makes the work harder to follow and more complicated to explain. However, the major conclusion, that cabp1 and 2 work together to reduce calcium-dependent inactivation of L-type calcium channels in cochlear inner hair cells, still holds.

Another weakness is that the authors used injections of AAV-containing sequences for Cabp2, but do not present data from sham surgeries. In most cases, the improvement of hearing function with AAV injection is believable and should be attributed to the cabp2 function. However, in at least one instance (Figure 4B), the results of the AAV injection experiments may be overinterpreted - the authors show that upon AAV injection, the hair cells have a much longer calcium current recovery following a large, long depolarization to inactivate the calcium channels. Without comparison to sham surgery, it is not known if this result could be a subtle result of the surgery or indeed due to the Cabp2 expression.<br /> It would be great to see the auditory nerve recordings in AAV-injected animals that have a recovery of ABRs. However, this is a challenging experiment that requires considerable time and resources, so is not required.

Reviewer #2 (Public Review):

Summary:<br /> Chen et al use human embryonic stem cells (ESCs) to determine the impact of wildtype MYC and a point mutant stable form of MYC (MYC-T58A) in the transformation of induced pulmonary neuroendocrine cells (PNEC) in the context of RB1/P53 (RP) loss (tumor suppressors that are nearly universally lost in small cell lung cancer (SCLC)). Upon transplant into immune-deficient mice, they find that RP-MYC and RP-MYC-T58A cells grow more rapidly, and are more likely to be metastatic when transplanted into the kidney capsule, than RP controls. Through single-cell RNA sequencing and immunostaining approaches, they find that these RPM tumors and their metastases express NEUROD1, which is a transcription factor whose expression marks a distinct molecular state of SCLC. While MYC is already known to promote aggressive NEUROD1+ SCLC in other models, these data demonstrate its capacity in a human setting that provides a rationale for further use of the ESC-based model going forward. Overall, these findings provide a minor advance over the previous characterization of this ESC-based model of SCLC published in Chen et al, J Exp Med, 2019.

The major conclusion of the paper is generally well supported, but some minor conclusions are inadequate and require important controls and more careful analysis.

Strengths:<br /> 1. Both MYC and MYC-T58A yield similar results when RP-MYC and RP-MYCT58A PNEC ESCs are injected subcutaneously, or into the renal capsule, of immune-deficient mice, leading to the conclusion that MYC promotes faster growth and more metastases than RP controls.

2. Consistent with numerous prior studies in mice with a neuroendocrine (NE) cell of origin (Mollaoglu et al, Cancer Cell, 2017; Ireland et al, Cancer Cell, 2020; Olsen et al, Genes Dev, 2021), MYC appears sufficient in the context of RB/P53 loss to induce the NEUROD1 state. Prior studies also show that MYC can convert human ASCL1+ neuroendocrine SCLC cell lines to a NEUROD1 state (Patel et al, Sci Advances, 2021); this study for the first time demonstrates that RB/P53/MYC from a human neuroendocrine cell of origin is sufficient to transform a NE state to aggressive NEUROD1+ SCLC. This finding provides a solid rationale for using the human ESC system to better understand the function of human oncogenes and tumor suppressors from a neuroendocrine origin.

Weaknesses:<br /> 1. There is a major concern about the conclusion that MYC "yields a larger neuroendocrine compartment" related to Figures 4C and 4G, which is inadequately supported and likely inaccurate. There is overwhelming published data that while MYC can promote NEUROD1, it also tends to correlate with reduced ASCL1 and reduced NE fate (Mollaoglu et al, Cancer Cell, 2017; Zhang et al, TLCR, 2018; Ireland et al, Cancer Cell, 2020; Patel et al, Sci Advances, 2021). Most importantly, there is a lack of in vivo RP tumor controls to make the proper comparison to judge MYC's impact on neuroendocrine identity. RPM tumors are largely neuroendocrine compared to in vitro conditions, but since RP control tumors (in vivo) are missing, it is impossible to determine whether MYC promotes more or less neuroendocrine fate than RP controls. It is not appropriate to compare RPM tumors to in vitro RP cells when it comes to cell fate. Upon inspection of the sample identity in S1B, the fibroblast and basal-like cells appear to only grow in vitro and are not well represented in vivo; it is, therefore, unclear whether these are transformed or even lack RB/P53 or express MYC. Indeed, a close inspection of Figure S1B shows that RPM tumor cells have little ASCL1 expression, consistent with lower NE fate than expected in control RP tumors.

In addition, since MYC appears to require Notch signaling to induce NE fate (Ireland et al), the presence of DAPT in culture could enrich for NE fate despite MYC's presence. It's important to clarify in the legend of Fig 4A which samples are used in the scRNA-seq data and whether they were derived from in vitro or in vivo conditions (as such, Supplementary Figure S1B should be provided in the main figure). Given their conclusion is confusing and challenges robustly supported data in other models, it is critical to resolve this issue properly. I suspect when properly resolved, MYC actually consistently does reduce NE fate compared to RP controls, even though tumors are still relatively NE compared to completely distinct cellular identities such as fibroblasts.

2. The rigor of the conclusions in Figure 1 would be strengthened by comparing an equivalent number of RP animals in the renal capsule assay, which is n = 6 compared to n = 11-14 in the MYC conditions.

3. Statistical analysis is not provided for Figures 2A-2B, and while the results are compelling, may be strengthened by additional samples due to the variability observed.

4a. Related to Figure 3, primary tumors and liver metastases from RPM or RPM-T58A-expressing cells express NEUROD1 by immunohistochemistry (IHC) but the putative negative controls (RP) are not shown, and there is no assessment of variability from tumor to tumor, ie, this is not quantified across multiple animals.

4b. Relatedly, MYC has been shown to be able to push cells beyond NEUROD1 to a double-negative or YAP1+ state (Mollaoglu et al, Cancer Cell, 2017; Ireland et al, Cancer Cell, 2020), but the authors do not assess subtype markers by IHC. They do show subtype markers by mRNA levels in Fig 4B, and since there is expression of ASCL1, and potentially expression of YAP1 and POU2F3, it would be valuable to examine the protein levels by IHC in control RP vs. RPM samples.

5. Given that MYC has been shown to function distinctly from MYCL in SCLC models, it would have raised the impact and value of the study if MYC was compared to MYCL or MYCL fusions in this context since generally, SCLC expresses a MYC family member. However, it is quite possible that the control RP cells do express MYCL, and as such, it would be useful to show.

Reviewer #2 (Public Review):

Summary:<br /> Bolamperti S. et al. 2023 investigate whether the expression of TG-interacting factor (Tgif1) is essential for osteoblastic cellular activity regarding morphology, adherence, migration/recruitment, and repair. Towards this end, germ-line Tgif1 deletion (Tgif1-/-) mice or male mice lacking expression of Tgif1 in mature osteoblastic and osteocytic cells (Dmp1-Cre+; Tgif1fl/fl) and corresponding controls were studied in physiological, bone anabolic, and bone fracture-repair conditions. Both Tgif1-/- and Dmp1-Cre+; Tgif1fl/fl exhibited decreased osteoblasts on cancellous bone surfaces and adherent to collagen I-coated plates. Tgif1-/- mice exhibit impaired healing in the tibial midshaft fracture model, as indicated by decreased bone volume (BV/Cal.V), osteoid (OS/BS), and low osteoblasts (number and surface). Likewise, both Tgif1-/- and Dmp1-Cre+; Tgif1fl/fl show impaired PTH 1-34, (100 µg/kg, 5x/wk for 3 wks) osteoblast activation in vivo, as detected by increases in quiescent bone surfaces. Mechanistic in vitro studies then utilized primary osteoblasts isolated from Tgif1-/- mice and siRNA Tgif1 knockdown OCY454 cells to further investigate and identify the downstream Tgif1 target driving these osteoblastic impairments. In vitro, Tgif1-/- osteoblastic and Tgif1 knockdown OCY454 cells exhibit decreased migration, abnormal morphology, and decreased focal adhesions/cells. Unexpectantly though, localization assays revealed Tgif1 to primarily concentrate in the nucleus and not to co-localize with focal adhesions (paxillin, talin). Also, the expression of major focal adhesion components (paxillin, talin, FAK, Src, etc.) or the Cdc42 family was not altered by loss of Tgif1 expression. In contrast, PAK3 expression is markedly upregulated by loss of Tgif1. In silico analysis followed by mechanistic molecular assays involving ChIP, siRNA (Tgif1, PAK3), and transfection (rat PAK3 promoter) techniques show that Tgif1 physically binds to a specific site in the PAK3 promoter region. Further, the knockdown of PAK3 rescues the Tgif1-deficient abnormal morphology in OCY454 cells. This is the first study to identify the novel transcriptional repression of PAK3 by Tgif1 as well as the specific Tgif1 binding site within the PAK3 promoter.

Strengths:<br /> This work has a plethora of strengths. The co-authors achieved their aim of eliciting the role of Tgif1 expression in osteoblastic cellular functions (morphology, spreading/attachment, migration). Further, this work is the first to depict the novel mechanism of Tgif1 transcriptional repression of PAK3 by a thorough usage of mechanistic molecular assays (in silico analysis, ChIP, siRNA, transfection etc.). The conclusions are well supported and justified by these findings, as the appropriate controls, sample sizes (statistical power), statistics, and assays were fully utilized.

The claims and conclusions are justified by the data.

Weaknesses:<br /> The discussion section could be expanded with a few sentences regarding limitations to the current study and potential future directions.

Reviewer #2 (Public Review):

This manuscript reports several interesting observations that invite follow-up. The notion that hubs, and perhaps condensates that may (or may not embrace them) are functionally and physiologically important is an open issue at this time. The authors note that TFIIIC helps to prune extraneous connections from hubs, but do not comment that the connections that are maintained are also reinforced. At the same time only modest changes in gene expression are associated with expanded or decreased connections and changes in bound proteins. One interesting possibility might be that standard methods for assessing expression miss changes in global or background transcription. It seems that the TFIIIC-MYCN-ER connection has features that would help to suppress such background. The results invite a more global consideration of TFIIIC than as primarily RNAPIII/small RNA transcription factor and of MYCN as an E-box dependent transcription factor. The results use state-of-the-art methods to develop interesting new ideas that have the potential to instruct further studies that may reveal new mechanisms of action for TFIIIC and MYCN

Strengths:<br /> Use of a variety of methods to assess the genomic response to increased MYCN in the presence or absence of TFIIIC. Establishes in vitro and in vivo the TFIIIC-MYCN complex.

Weaknesses:<br /> Dynamic inferences are made without kinetic experiments.

Reviewer #2 (Public Review):

Summary:<br /> To investigate the evolutionary relationship between the RNAi pathway and innate immunity, this study uses biochemistry and structural biology to investigate the trimeric complex of Dicer-1, DRH-1 (a RIGI homologue), and RDE-4, which exists in C. elegans. The three subunits were co-expressed to promote stable purification of the complex. This complex promoted ATP-dependent cleavage of blunt-ended dsRNAs. A detailed kinetic analysis was also carried out to determine the role of each subunit of the trimeric complex in both the specificity and efficiency of cleavage. These studies indicate that RDE-4 is critical for cleavage while DRC-1 is primarily involved in the specificity of the reaction, and DRH-1 promotes ATP hydrolysis. Finally, a moderate density (6-7 angstrom) cryo-EM structure is presented with attempts to position each of the components.

Strengths:<br /> 1. Newly described methods for studying the C. elegans DICER complex.<br /> 2. New structure, albeit only moderate resolution.<br /> 3. Kinetic study of the complex in the presence and absence of individual subunits and mutations, provides detailed insight into the contribution of each subunit.

Weaknesses:<br /> 1. Limited insight due to limited structural resolution.<br /> 2. No attempts to extend findings to other Dicer or RLR systems.

Reviewer #2 (Public Review):

The authors identify a third component in the interaction between myosin Va and melanophilin- an interaction between a 32-residue sequence encoded by exon-g in myosin Va and melanophilin's actin-binding domain. This interaction has implications for how melanosome motility may be regulated.

While this work is largely well done, I believe that additional work would be required to make a more compelling case (e.g. some affinity measurements, necessary controls for the dominant negative experiments).  First, the study provides just one more piece to a well-developed story (the role of exon-F and the GTD in myosin Va: melanophilin (Mlph) interaction), much of which was published 20 years ago by several labs. Second, the study does not demonstrate a physiological significance for their findings other than that exon-G plays an auxiliary role in the binding of myosin Va to Mlph. For example, what dictates the choice between Mlph's actin binding domain (ABD) binding to actin or to exon-G. Is it a PTM or local actin concentration? It is unlikely to be alternative splicing as exon-G is present in all spliced isoforms of myosin Va. And what changes re melanosome dynamics in cells between these two alternatives? Similarly, the paper does not provide any in vitro evidence that binding to exon-G instead of actin effects the processivity of a Rab27a/Myosin Va/Mlph transport complex. For example, if the ABD sticks to exon-G instead of actin, does that block Mlph's ability to promote processivity through its interaction with the actin filament during transport? In summary, given that the authors did not directly test their model either in vitro or in cells, I do not think this story represent a significant conceptual advance.

Reviewer #2 (Public Review):

Summary:<br /> In this work, the authors manage to optimize a simple and rapid protocol using SEC followed by DGCU to isolate sEVs with adequate purity and yield from small volumes of plasma. Isolated fractions containing sEVs using SEC, DGCU, SEC-DGCU, and DGCU-SEC are compared in terms of their yield, purity surface protein profile, and RNA content. Although the combined use of these methodologies has already been evaluated in previous works, the authors manage to adapt them for the use of small volumes of plasma, which allows working in 1.5 mL tubes and reducing the centrifugation time to 2 hours.

The authors finally find that although both the SEC-DGCU and DGCU-SEC combinations achieve isolates with high purity, the SEC-DGCU combination results in higher yields.

This work provides an interesting tool for the rapid obtention of sEVs with sufficient yield and purity for detailed characterization which could be very useful in research and clinical therapy.

Strengths:<br /> -The work is well-written and organized.<br /> -The authors clearly state the problem they want to address, that is, optimizing a method that allows sEV to be isolated from small volumes of plasma.<br /> -Although these methodologies have been tested in previous works, the authors manage to isolate sEVs of high purity and good performance through a simple and fast methodology.<br /> -The characteristics of all isolated fractions are exhaustively analyzed through various state-of-the-art methodologies.<br /> -They present a good interpretation of the results obtained through the methodologies used.

Weaknesses:<br /> -Lack of references that support some of the results obtained.<br /> -Although this work focuses on comparing different techniques and their combinations to find an optimal option, the authors do not use any statistical method that reliably shows the differences between these techniques, except when repeatability is measured.

Reviewer #2 (Public Review):

Summary:<br /> Recent advances in single cell profiling of gene expression (RNA) permit the analysis of specialized cell types, an approach that has great value in the nervous system which is characterized by prodigious neuronal diversity. The novel data in this study focus primarily on genetic profiling to compare autonomic neurons from ganglia associated with the cranial parasympathetic outflow (sphenopalatine (also known as pteropalatine), the thoraco-lumbar sympathetic outflow (stellate, coeliac) and the sacral parasympathetic outflow (pelvic). Using statistical methods to reduce the dimensionality of the data and map gene expression, the authors provide interesting evidence that cranial parasympathetic and sacral sympathetic ganglia differ from each other and from sympathetic ganglia (Figures 1, S1 - S4). The authors interpret the mapping analysis as evidence that the cranial and sacral outflows differ so that calling them both parasympathetic is unjustified. Based on anatomical localization of markers (Figure 2 ) (mainly transcription factors) the authors show a similarity between the sympathetic and pelvic ganglion. In Figure 3 they present evidence that some pelvic ganglionic neurons are dually innervated by sympathetic preganglionic neurons and sacral preganglionic neurons. These observations are interpreted to mean that the pelvic ganglion is not parasympathetic, but rather a modified sympathetic ganglion - hence the title of the manuscript.

Strengths:<br /> The extensive use of single cell profiling in this work is both interesting and exciting. Although still in its early stages, it holds promise for a deepened understanding of autonomic development and function. As noted in the introduction, this study extends previous work by Professor Brunet and his associates.

Weaknesses:<br /> This work further documents differences between the cranial and sacral parasympathetic outflows that have been known since the time of Langley - 100 years ago. The approach taken by Brunet et al. has focused on late neonatal and early postnatal development, a time when autonomic function is still maturing. In addition, the sphenopalatine and other cranial ganglia develop from placodes and the neural crest, while sympathetic and sacral ganglia develop from the neural crest alone. How then do genetic programs specifying brainstem and spinal development differ and how can this account for kinship that Brunet documents between spinal and sacral ganglia? One feature that seems to set the pelvic ganglion apart is the mixture of 'sympathetic' and 'parasympthetic' ganglion cells and the convergence of preganglionic sympathetic and parasympathetic synapses on individual ganglion cells (Figure 3). This unusual organization has been reported before using microelectrode recordings (see Crowcroft and Szurszewski, J Physiol (1971) and Janig and McLachlan, Physiol Rev (1987)). Anatomical evidence of convergence in the pelvic ganglion has been reported by Keast, Neuroscience (1995). It should also be noted that the anatomy of the pelvic ganglion in male rodents is unique. Unlike other species where the ganglion forms a distributed plexus of mini-ganglia, in male rodents the ganglion coalesces into one structure that is easier to find and study. Interestingly the image in Figure 3A appears to show a clustering of Chat-positive and Th-positive neurons. Does this result from the developmental fusion of mini ganglia having distinct sympathetic and parasympathetic origins. In addition, Brunet et al dismiss the cholinergic and noradrenergic phenotypes as a basis for defining parasympathetic and parasympathetic neurons. However, see the bottom of Figure S4 and further counterarguments in Horn (Clin Auton Res (2018)). What then about neuropeptides, whose expression pattern is incompatible with the revised nomenclature proposed by Brunet et al.? Figure 1B indicates that VIP is expressed by sacral and cranial ganglion cells, but not thoracolumbar ganglion cells. The authors do not mention neuropeptide Y (NPY). The immunocytochemistry literature indicates that NPY is expressed by a large subpopulation of sympathetic neurons but never by sacral or cranial parasympathetic neurons.

The title of this paper is misleading because it implies a conclusion that is not adequately supported by the data and that is difficult for a general reader to parse. Independent assessments by two referees both agreed on title's problematic message. If one can get beyond the title, then the paper does contain data that is of interest. The authors compared single cell gene expression in neurons from the cranial sphenopalatine ganglion, paravertebral chain ganglia (stellate and lumbar), the prevertebral coeliac ganglion and the bladder ganglion. The cranial and pelvic ganglia are parasympathetic, while the paravertebral and prevertebral ganglia are sympathetic. The gene expression data identified differences between the cranial, sympathetic, and pelvic ganglia. Based primarily on this finding the authors concluded that the sacral bladder ganglion is not parasympathetic. Since some genes suggest a kinship between the pelvic and sympathetic neurons, the authors conclude that the pelvic neurons are pelvo-sympathetic - hence the title. This nomenclature does little to improve understanding of the autonomic motor system and it ignores important anatomical and functional properties that underlie existing definitions of the sympathetic and parasympathetic systems. The idea that the cranial and sacral autonomic outflows have some differences is not new (see for example Nilsson, 1983 and Janig, 2022). Since many of the genes identified in the present study are HOX genes and other transcription factors that specify the rostro-caudal axis during development, it is also not surprising that these genes suggest a kinship between sacral parasympathetic neurons and sympathetic neurons, all of which derive from the neural crest and are supplied by the spinal cord. The different profile of cranial parasympathetic neurons is also not surprising given that they derive from a mixture of placodal and neural crest progenitors and are supplied by the brainstem. (see my previous comments for anatomical and functional criteria that further support the existing nomenclature for the sympathetic and parasympathetic motor systems.

Reviewer #2 (Public Review):

The authors make the interesting observation that the developmental refinement of apical M/T cell dendrites into individual glomeruli proceeds normally even when the majority of neighboring M/T cells are ablated. At later stages, the remaining neurons develop additional dendrites that invade multiple glomeruli ectopically, and similarly, OSN inputs to glomeruli lose projection specificity as well. The authors conclude that the normal density of M/T neurons is not required for developmental refinement, but rather for maintaining specific connectivity in adults.

The observations are indeed quite striking; however, the authors' conclusions are not entirely supported by the data.

1. It is unclear whether the expression of diphtheria toxin that eventually leads to the ablation of the large majority of M/T neurons compromises the cell biology of the remaining ones.

2. The authors interpret the growth of ectopic dendrites later in life as a lack of maintenance of dendrite structure; however, maybe the observed changes reflect actually adaptations that optimize wiring for extremely low numbers of M/T neurons. The finding that olfactory behavior was less affected than predicted supports this interpretation.

3. The number of remaining M/T neurons is much higher at P10 than later. Can the relatively large number of remaining neurons (or their better health status) be the reason that dendrites refine normally at the early developmental stages rather than a (currently unknown) developmental capacity that preserves refinement?

4. While the effect of reduced M/T neuron density on both M/T dendrites and OSN axons is described well, the relationship between both needs to be characterized better: Is one effect preceding the other or do they occur simultaneously? Can one be the consequence of the other?

5. Page 7: the observation that not all neurons develop additional dendrites is not a sign of differences between cell types, it may be purely stochastic.

6. Page 8: the fact that activity blockade did not affect the formation of ectopic dendrites does not suggest that the process is not activity-dependent: both manipulations have the same effect and may just mask each other.

7. It remains unclear how the observed structural changes can explain the behavioral effects.

Reviewer #2 (Public Review):

Carla de la Fuente et al., utilize a diversity of approaches to understand which plant traits contribute to the stress resilience of pearl millet in the Sahelian desert environment. By comparing data resulting from crop modeling of pearl millet growth and meteorological data from a span of 20 years, the authors clearly determined that early season drought resilience is contributed by accelerated growth of the seedling primary root, which confirms a hypothesis generated in a previous study, Passot et al., 2016. To determine the genetic basis for this trait, they performed a combination of GWAS, QTL analysis, and RNA sequencing and identified a previously unannotated coding sequence of a glutaredoxin C9-like protein, PgGRXC9, as the strongest candidate. Phenotypic analysis using a mutant of the closest Arabidopsis homolog AtROXY19 suggests the broad conservation of this pathway. Comparisons between the transcript of PgGRXC9 by in situ hybridization (this work) and AtROXY19 pattern expression (Belin et al., 2014) support the hypothesis that this pathway acts in the elongation zone of the root. Additional analysis of cell production and elongation rates in root apex in both pearl millet and A. thaliana suggests that PgGRXC9 specifically regulates primary root through the promotion of cell elongation. While several studies have established the connection between redox status of cells and root growth, the current study represents an important contribution to the field because of the agricultural importance of the plant studied, and the connection made between this developmental trait and stress resilience in a specific and stressful environmental context of the Sahelian desert.

Reviewer #2 (Public Review):

Summary:<br /> The evolution of resistance to antimalarial drugs follows a seemingly counterintuitive pattern, in which resistant strains typically originate in regions where malaria prevalence is relatively low. Previous investigations have suggested that frequent exposures in high-prevalence regions produce high levels of partial immunity in the host population, leading to subclinical infections that go untreated. These subclinical infections serve as refuges for sensitive strains, maintaining them in the population. Prior investigations have supported this hypothesis; however, many of them excluded important dynamics, and the results cannot be generalized. The authors have taken a novel approach using a deterministic model that includes both general and adaptive immunity. They find that high levels of population immunity produce refuges, maintaining the sensitive strains and allowing them to outcompete resistant strains. While general population immunity contributed, adaptive immunity is key to reproducing empirical patterns. These results are robust across a range of fitness costs, treatment rates, and resistance efficacies. Given sufficient antigenic diversity and high transmission, sensitive parasites remain in circulation even when there is no cost to resistance. This work demonstrates that future investigations cannot overlook adaptive immunity and antigenic diversity.

Strengths:<br /> Overall, this is a very nice paper that makes a significant contribution to the field. It is well-framed within the body of literature and achieves its goal of providing a generalizable, unifying explanation for otherwise disparate investigations. The model is innovative. The approach is elegant and rigorous, with results that are supported across a broad range of parameters when considered within an equilibrium setting. Their exploration of geographical patterns of resistance makes the results of their simulations even more compelling. As such, this work will likely serve as a foundation for many future investigations.

Weaknesses:

Although the authors model resistance invasion, it does not align with empirical observations of the spread of resistance. For example, Plasmodium's mutation rate and population size mean that mutations providing chloroquine resistance should arise repeatedly even within a single infection. Nevertheless, Africa remained free of chloroquine resistant strains until a lineage was introduced from Asia. Upon introduction, it spread across the continent within ten years. The difference between the fate of chloroquine resistance originating in Africa versus chloroquine resistance originating in Asia cannot be attributed to changes in population immunity and treatment.

The source of this disparity may be in part attributable to the use of a deterministic, compartmental model, as the authors mention in the discussion. Strains are not explicitly modeled. This means that in terms of the distribution of strain diversity, the resistant and the sensitive compartments are identical, and the locus determining resistance is equally distributed across all strain backgrounds. However, substantial rates of linkage disequilibrium and clonal reproduction are found even in high transmission settings. The model assumptions may be met at equilibrium, but are not appropriate for most scenarios involving the invasion of a rare mutation.

Reviewer #2 (Public Review):

Summary:<br /> The landmark publication of the "Fly Atlas" in 2022 provided a single cell/nuclear transcriptomic dataset from 15 individually dissected tissues, the entire head, and the body of male and female flies. These data led to the annotation of more than 250 cell types. While certainly a powerful and data-rich approach, a significant step forward relies on mapping these data back to the organism in time and space. The goal of this manuscript is to map 150 transcripts defined by the Fly Atlas by FISH and in doing so, provide, for the first time, a spatial transcriptomic dataset of the adult fly. Using this approach (Molecular Cartography with Resolve Biosciences), the authors, furthermore, distinguish different RNA localizations within a cell type. In addition, they seek to use this approach to define previously unannotated clusters found in the Fly Atlas. As a resource for the community at large interested in the computational aspects of their pipeline, the authors compare the strengths and weaknesses of their approach to others currently being performed in the field.

Strengths:<br /> 1. The authors use Resolve Biosciences and a novel bioinformatics approach to generate a FISH-based spatial transcriptomics map. To achieve this map, they selected 150 genes (50 body; 100 head) that were highly expressed in the single nuclear RNA sequencing dataset and were used in the 2022 paper to annotate specific cell types; moreover, the authors chose several highly expressed genes characteristic of unannotated cell types. Together, the approach and generated data are important next steps in translating the transcriptomic data to spatial data in the organism.<br /> 2. Working with Resolve, the authors developed a relatively high throughput approach to analyze the location of transcripts in Drosophila adults. This approach confirmed the identification of particular cell types suggested by the FlyAtlas as well as revealed interesting subcellular locations of the transcripts within the cell/tissue type. In addition, the authors used co-expression of different RNAs to unbiasedly identify "new cell types". This pipeline and data provide a roadmap for additional analyses of other time points, female flies, specific mutants, etc.<br /> 3. The authors show that their approach reveals interesting patterns of mRNA distribution (e.g alpha- and beta-Trypsin in apical and basal regions of gut enterocytes or striped patterns of different sarcomeric proteins in body muscle). These observations are novel and reveal unexpected patterns. Likewise, the authors use their more extensive head database to identify the location of cells in the brain. They report the resolution of 23 clusters suggested by the single-cell sequencing data, given their unsupervised clustering approach. This identification supports the use of spatial cell transcriptomics to characterize cell types (or cell states).<br /> 4. Lastly, the authors compare three different approaches --- their own described in this manuscript, Tangram, and SpaGE - which allow integration of single cell/nuclear RNA-seq data with spatial localization FISH. This was a very helpful section as the authors compared the advantages and disadvantages (including practical issues, like computational time).

Weaknesses:<br /> 1. Experimental setup. It is not clear how many and, for some of the data, the sex of the flies that were analyzed. It appears that for the body data, only one male was analyzed. For the heads, methods say male and female heads, but nothing is annotated in the figures. As such, it remains unclear how robust these data are, given such a limited sample from one sex. As such, the claims of a spatial atlas of the entire fly body and its head ("a rosetta stone") are overstated. Also, the authors should clearly state in the main text and figure legends the sex, the age, how many flies, and how many replicates contributed to the data presented (not just the methods). What also adds to the confusion is the use of "n" in para 2 of the results. " ... we performed coronal sections at different depths in the head (n=13)..." 13 sections in total from 1 head or sections from 13 heads? Based on the body and what is shown in the figure, one assumes 13 sections from one head. Please clarify.<br /> 2. Probes selected: Information from the methods section should be put into the main text so that it is clear what and why the gene lists were selected. The current main text is confusing. If the authors want others to use their approach, then some testing or, at the very least, some discussion of lower expressed genes should be added. How useful will this approach be if only highly expressed genes can be resolved? In addition, while it is understood that the company has a propriety design algorithm for the probes, the authors should comment on whether the probes for individual genes detect all isoforms or subsets (exons and introns?), given the high level of splicing in tissues such as muscle.<br /> 3. Imaging: it isn't clear from the text whether the repeated rounds of imaging impacted data collection. In many of what appear to be "stitched" images, there are gradients of signal (eg, figure 2F); please comment. Also, since this a new technique, could a before and after comparison of the original images and the segmented images be shown in the supplemental data so that the reader can better appreciate how the authors assessed/chose/thresholded their data? More discussion of the accuracy of spot detection would be helpful.<br /> 4. The authors comment on how many RNAs they detected (first paragraph of results). How do these numbers compare to the total mRNA present as detected by single-cell or single-nuclear sequencing?<br /> 5. Using this higher throughput method of spatial transcriptomics, the authors discern different cell types and different localization patterns within a tissue/cell type.<br /> a. The authors should comment on the resolution provided by this approach, in terms of the detection of populations of mRNAs detected by low throughput methods, for example, in glia, motor neuron axons, and trachea that populate muscle tissue. Are these found in the images? Please show.<br /> b. The authors show interesting localization patterns in muscle tissue for different sarcomere protein-coding mRNAs, including enrichment of sls in muscle nuclei located near the muscle-tendon attachment sites. As this high throughput approach is newly being applied to the adult fly, it would increase confidence in these data, if the authors would confirm these data using a low throughput FISH technique. For example, do the authors detect such alternating "stripes" ( Act 88F, TpnC4, and Mhc) or enriched localization (sls) using FISH that doesn't rely on the repeated colorization, imaging, decolorization of the probes?<br /> 6. The authors developed an unbiased method to identify "new cell types" which relies on co-expression of different transcripts. Are these new cell types or a cell state? While expression is a helpful first step, without any functional data, the significance of what the authors found is diminished. The authors need to soften their statements.

Appraisal:<br /> The authors' goal is to map single cell/nuclear RNAseq data described in the 2022 Fly Atlas paper spatially within an organism to achieve a spatial transcriptomic map of the adult fly; no doubt, this is a critical next step in our use of 'omics approaches. While this manuscript does the hard work of trying to take this next step, including developing and testing a new pipeline for high throughput FISH and its analysis, it falls short, in its present form, in achieving this goal. The authors discuss creating a robust spatial map, based on one male fly. Moreover, they do not reveal principles of mRNA localization, as stated in the abstract; they show us patterns, but nothing about the logic or function of these patterns. This same criticism can be said of the identification of "new cell types, just based on RNA colocalization. In both cases (mRNA subcellular localization or cell type identification), further data in the form of validation with traditional low throughput FISH and genetic manipulations to assess the relation to cell function are required for the authors to make such claims.

Discussion of likely impact:<br /> If revised, these data, and importantly the approach, would impact those working on Drosophila adults as well as those working in other model systems where single cell/nuclear sequencing is being translated to the spatial localization within the organism. The subcellular localization data - for example, the size of transcripts and how that relates to localization or the patterns of sarcomeric protein localization in muscle - are intriguing, and would likely impact our thinking on RNA localization, transport, etc if confirmed. Lastly, the authors compare their computational approaches to those available in the field; this is valuable as this is a rapidly evolving field and such considerations are critical for those wishing to use this type of approach.

Reviewer #2 (Public Review):

Summary: This work presented by Kudo and colleagues is of great importance to strengthen our understanding of electrophysiological changes in the course of AD. Although the main conclusions regarding functional connectivity and spectral power change through the course of the disease are not new and have been largely studied and theorised on, this article offers an innovative approach that certainly consolidates previous knowledge on the topic. Not only that, this article also broadens our knowledge presenting useful and important details on the specificity of frequency and cortical distribution of these early alterations. The main take-home message of this work is the early disruption of electrophysiological signatures that precedes detectable alterations in other more commonly used pathology markers (i.e. gray matter atrophy and cognitive impairment). More specifically, these signatures include long-range connectivity in the alpha and beta bands, and local synchrony (spectral power) in the same frequency bands.

Strengths: The present work has some major strengths that make it paramount for the advance of our understanding of AD electrophysiology. It is a very well written manuscript that, despite the complexity of the analyses employed, runs the reader through the different steps of the analysis in a pedagogic and clever way, making the points raised by the results easy to grasp. The methodology itself is carefully chosen and appropriate to the nature of the question posed by the researchers, as event-based models are well-suited for cross-sectional data.

The quality of the figures is outstanding; not only are they aesthetic but, more importantly, the figures convey information exceptionally well and facilitate comprehension of the main results.<br /> The conclusions of the paper are, in general, well described and discussed, and consider the state-of-the-art works of AD electrophysiology. Furthermore, even though the conclusions themselves are not groundbreaking at all (synaptic damage preceding structural and cognitive impairment is one of the epitomes of the pathological cascading model proposed by Jack in 2010), this article is innovative and groundbreaking in the way they address with clever analyses in a relatively large sample for neuroimaging standards.

Weaknesses: The authors increased the clarity of sample description after revisions (particularly control group characterization). However, even though it is true that a certain percentage of AB positivity is to be expected amongst cognitively healthy individuals, that doesn´t discard they are not expressing preclinical AD to some extent. I still feel that including only biomarker negative participants in the control group would increase the quality of the work. However, the sample is relatively well characterized as a whole and the results are interesting and in line with previous literature, thus limiting the apparent impact of these possible confounds.

Reviewer #2 (Public Review):

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

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

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

Major Weaknesses<br /> - There is no evidence of actual transparency of the entire mouse brain across different treatments. The suggested protocol is very good at removing lipids (as assessed by DiD staining) and by results of fluorescence registration deep within the brain. BUT, since in many places of the manuscript authors speak of "transparency" the reader will expect the typical picture in which control and processed brains are on top of a white graphical pattern that would evidence transparency (see as an example Figure 1 and 2 of Wan et al. 2018 (Neurophotonics. 2018 Jul;5(3):035007. doi: 10.1117/1.NPh.5.3.035007.)

- The manuscript lacks clarity on the applicability of MOCAT to regular formalin-fixed tissue and tissues other than the brain.

- Insufficient information is provided on the "epoxy treatment" or "hydrogel," and a more detailed explanation is warranted.

- The differences between passive and active immunolabeling, as well as photobleaching data, should be addressed for a comprehensive understanding.

- The assertion that MOCAT can be rapidly applied in hospital pathology departments seems overstated due to the limited availability of light-sheet microscopes outside research labs.

- The compatibility of MOCAT with genetically encoded fluorescent proteins remains unclear and warrants further investigation.

- The control of equivalent depths in cryosections for evaluating the intensity of DiD staining should be elaborated upon.

- The composition of NFC1 and NFC2 solutions for refractive index matching should be provided.

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

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

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

Reviewer #2 (Public Review):

Summary:<br /> The goals of this study were to develop a genetic approach that would specifically and comprehensively target axo-axonic cells (AACs) throughout the brain and then to describe the patterns and characteristics of the targeted AACs in multiple, selected brain regions. The investigators have been successful in providing the most complete description of the regional distribution of putative (pAACs) throughout the brain to date. The supporting evidence is convincing, even though incomplete in some brain regions. The findings should serve as a guide for more detailed studies of AACs within each brain region and lead to new insights into the connectivity and functional organization of this important group of GABAergic interneurons.

Strengths:<br /> The study has numerous strengths. A major strength is the development of a unique intersectional genetic strategy that uses cell lineage (Nkx2.1) and molecular (Unc5b or Pthlh) markers to identify axo-axonic AACs specifically and, apparently, nearly completely throughout the mouse brain. While AACs have been described previously in the cerebral cortex, hippocampus, and amygdala, there has been no specific genetic marker that selectively identifies all AACs in these regions.

The current genetic strategy has labeled pAACs in a large number of additional brain regions, including the claustrum-insular complex, extended amygdala, and several olfactory centers. In general, the findings provide support for the specificity of the methods for targeting AACs, and include some examples of labeling near markers of axon initial segments. However, the Investigators are careful to refer to labeled neurons as "putative AACs" as they have not been fully characterized and their identity verified.

The descriptions and numerous low-magnification images of the brain provide a roadmap for subsequent, detailed studies of AACs in numerous brain regions. The overview and summaries of the findings in the Abstract, Introduction, and Discussion are particularly clear and helpful in placing the extensive regional descriptions of AACs in context.

Weaknesses:<br /> One weakness of the study is the lack of an illustration of the high-resolution cell labeling that can be achieved with the methods, including labeling of numerous rows of axon terminals in contact with axon initial segments. The initial images of the brain-wide distribution of putative AACs are necessarily presented at low magnification. Although the authors indicate that the cells have "highly characteristic AAC labeling patterns throughout the neocortex, hippocampus and BLA", these morphological details cannot be visualized by the reader at the current magnification, even when the images are enlarged on the computer screen. Some of the details become evident in later Figures, but an initial illustration of single cell labeling with confocal microscopy, or tracing of their characteristic axonal arbors, would support the specificity of the labeling in the low magnification images.

Table 1 indicates that the AAC identity of the cells has been validated in many brain regions but not in all. The methods used for validation have not been described and should be included for completeness. The authors are careful to acknowledge that labeled cells in some regions have not been validated and refer to such cells as pAACs.

The intersectional genetic methods included the use of the lineage marker Nkx2.1 with either Unc5b or Pthlh as the molecular marker. As described, the mice with intersectional targeting of Nkx2.1 and Unc5b appear to show the most specific brain-wide labeling for AACs, and the majority of the descriptions are from these mice. The targeting with Nkx2.1 and Pthlh is less convincing. The title for Figure 1 Supplemental Figure 3 suggests a similar AAC distribution in the Pthlh;Nkx2.1 mouse compared to the Unc5b;Nkx2.1 mouse. However, the descriptions of the individual panels suggest a number of inconsistencies and non-AAC labeling. The heavy labeling in the caudate and cells in layer 4 is particularly problematic. Based on the data presented, it appears that heavy labeling achieved in these mice could not be relied on for specific labeling of all AACs, although specific labeling could be achieved under some conditions, such as following tamoxifen administration at select ages.

The methods described for dense labeling and single-cell labeling are described briefly in the methods. Some discussion of the development of the methods would be useful, including how it was determined that methods for heavy labeling identified AACs specifically and completely.

Reviewer #2 (Public Review):

Summary:

Jeong & Choi (2023) use a semi-naturalistic paradigm to tackle the question of how the activity of neurons in the mPFC might continuously encode different functions. They offer two possibilities: either there are separate dedicated populations encoding each function, or cells alter their activity depending on the current goal of the animal. In a threat-avoidance task rats procured sucrose in an area of a chamber where, after remaining there for some amount of time, a 'Lobsterbot' robot attacked. To initiate the next trial rats had to move through the arena to another area before returning to the robot encounter zone. Therefore the task has two key components: threat avoidance and navigating through space. Recordings in the IL and PL of the mPFC revealed encoding that depended on what stage of the task the animal was currently engaged in. When animals were navigating, neuronal ensembles in these regions encoded distance from the threat. However, whilst animals were directly engaged with the threat and simultaneously consuming reward, it was possible to decode from a subset of the population whether animals would evade the threat. Therefore the authors claim that neurons in the mPFC switched between two functional modes: representing allocentric spatial information, and representing egocentric information pertaining to the reward and threat.

Strengths:

As the authors point out, whilst these multiple functions of activity in the mPFC have generally been observed in tasks dedicated to the study of a singular function, less work has been done in contexts where animals continuously switch between different modes of behaviour in a more natural way. Being able to assess whether previous findings of mPFC function apply in natural contexts is very valuable to the field, even outside of those interested in the mPFC directly. This also speaks to the novelty of the work; although mixed selectivity encoding of threat assessment and action selection has been demonstrated in some contexts (e.g. Grunfeld & Likhtik, 2018) understanding the way in which encoding changes on-the-fly in a self-paced task is valuable for verifying whether current understanding holds true.

The authors are also generally thoughtful in their analyses and use a variety of approaches to probe the information encoded in the recorded activity. In particular, they also use relatively close analysis of behaviour as well as manipulating the task itself by removing the threat to verify their own results. The use of such a rich task also allows them to draw comparisons, e.g. in different zones of the arena or different types of responses to threats, that a more reduced task would not otherwise allow.

Weaknesses:

The central question the paper seeks to answer is whether 'individual cells are dedicated to spatial representation and emotional stimuli processing or if they adapt their function to the current goal'. However, there does not seem to be a direct analysis that answers this question. It is not clear what proportion of each of the ensembles recorded is necessary for decoding distance from the threat, and whether it is these same neurons that directly 'switch' to responding to head entry or withdrawal in the encounter phase within the total population. The PCA gets closest to answering this question by demonstrating that activity during the encounter is different from activity in the nesting or foraging zones, but in principle this could be achieved by neurons or ensembles that did not encode spatial parameters. The population analyses are focused on neurons sensitive to behaviours relating to the threat encounter, but even before dividing into subtypes etc., this is at most half of the recorded population. And again it is difficult to ascertain how the final ensemble analysis of the avoidance response relates to the prior spatial encoding. As a result, the model of the results proposed in Fig. 7 cannot be validated by the data as is.

A second concern is also illustrated by Fig. 7: in the data presented, separate reward and threat encoding neurons were not shown - in the current study design, it is not possible to dissociate reward and threat responses as the data without the threat present were only used to study spatial encoding integrity. To be able to claim this working model, a key additional analysis is to compare PETHs around head entry and withdrawal for sucrose without attack. Alternatively, a small proportion of probe trials could have been added where rats did not receive any reward for being in the encounter zone. This would allow the authors to ascertain whether the elevated response of the Type 2 neurons in particular is partially driven by reward receipt.

Thirdly, the findings of this work are not mechanistic or functional but are purely correlational. For example, it is claimed that analysing activity around the withdrawal period allows for ascertaining their functional contributions to decisions. But without a direct manipulation of this activity, it is difficult to make such a claim. The authors later discuss whether the elevated response of Type 2 neurons might simply represent fear or anxiety motivation or threat level, or whether they directly contribute to the decision-making process. As is implicit in the discussion, the current study cannot differentiate between these possibilities. However, the language used throughout does not reflect this.

Fourthly, the authors mention the representation of different functions in 'distinct spatiotemporal regions' but the bulk of the analyses, particularly in terms of response to the threat, do not compare recordings from PL and IL although - as the authors mention in the introduction - there is prior evidence of functional separation between these regions.

Reviewer #2 (Public Review):

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

To isolate and study the neural crest, the authors use embryos collected at E8.5 from two transgenic mouse lines. Wnt1-Cre;RosaeYFP labels Wnt1-positive neuroepithelial cells in the dorsolateral neural plate, which includes pre-migratory neural crest that resides in the dorsal neuroectoderm and neural plate border before induction (as well as some other lineages). Mef2c-F10N-LacZ leverages a neural crest cell-specific enhancer of Mef2c to control LacZ expression in the predominantly migratory neural crest. This dual genetic approach that allows the authors to distinguish and compare pre-migratory and migratory neural crest cells is a strength of the work. However, one potential weakness needing to be addressed is that some workers (e.g., Lewis et al., 2013) have reported phenotypic effects of Wnt1-Cre transgene expression including ectopic Wnt pathway activation, abnormal neuroepithelial development, and increases in CyclinD1 expression and cell proliferation. The authors should discuss the extent to which the results of their study were or were not influenced by these potentially confounding effects, especially since Wnt canonical signaling is known to regulate the G1/S transition and promote delamination of the neural crest.

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

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

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

As a functional test of their hypothesis that one of the genes indicative of an EMT intermediate stage (i.e., Dlc1) is essential for neural crest migration, the authors use a lentivirus-mediated knockdown strategy. A strength is that the authors include appropriate scramble and cell death controls as part of their experimental design. However, a weakness is that the authors do not justify why they chose a knockdown strategy, which has its limitations including its systemic injection into the amniotic cavity, its likely global and more variable effects, and its need to be conducted in culture. Why the authors did not instead use a Wnt1-Cre-mediated deletion of Dlc1, which would have been "cleaner" and more specific to the neural crest, is not clear (maybe so they could specifically target different Dcl1 isoforms?). Also, the authors use Sox10 as a marker to count neural crest cells, but Sox10 may only label a subset of neural crest cells and thus some unaffected lineages may not have been counted. The authors should mention what is known about the regulation of Dcl1 by Sox10 in the neural crest. Although the data are persuasive, a second marker for counting neural crest cells following knockdown would make the analysis more robust. Can the authors explain why they did not simply use the Mef2c-F10N-LacZ line and count LacZ-positive cells (if fluorescence signal was required for the quantification workflow, then could they have used an anti-beta Galactosidase antibody to label cells)?

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

Reviewer #2 (Public Review):

Summary:<br /> Chromosome organization in E. coli and related species ('transversal') deviates starkly from the pattern more commonly found in bacteria ('longitudinal'). The underlying mechanisms and the physiological roles, however, are not well understood. This manuscript by Seba et al. investigates the activity and regulation of MukBEF in chromosome folding in E. coli. Using a construct for inducible expression of MukBEF, the authors first demonstrate that the initiation of long-range chromosome contacts (likely by loop extrusion) is not restricted to few positions on the chromosome and rather widely distributed but excluding the replication terminus region. Using ChIP-Seq, the authors show that the distribution of MukBEF over the chromosome is consistent with widely distributed loading and moreover indicate a connection of chromosome folding and DNA replication with newly replicated DNA shower an increased tendency for MukBEF binding. To dissect this further, they then redistribute matS sites on the chromosome by a clever strategy based on large-scale transpositions. The results reveal that matS-free DNA segments undergo MukBEF dependent folding regardless of their position relative to the origin of replication, being consistent with a broad distributed loading of MukBEF. By fine-mapping with smaller transposition events, they show that few matS sites are sufficient to impede MukBEF activity. Surprisingly, however, E. coli and most related genomes harbor many matS sites, which are particularly highly concentrated near the chromosome dimer resolution dif site (Fig. 5).

This is a well-executed and well-presented study. The findings show that the MatP/matS system acts locally and independent of DNA replication to restrict MukBEF in the replication terminus region. Few of the many matS sites are sufficient for MukBEF restriction. The main conclusions of the work are clear and well supported by the data.

Reviewer #2 (Public Review):

Summary<br /> The authors seek to characterize the role of splicing factor SRSF1 during spermatogenesis. Using a conditional deletion of Srsf1 in germ cells, they find that SRSF1 is required for male fertility. Via immunostaining and RNA-seq analysis of the Srsf1 conditional knockout (cKO) testes, combined with SRSF1 CLIP-seq and IP-MS data from the testis, they ultimately conclude that Srsf1 is required for homing of precursor spermatogonial stem cells (SCCs) due to alternative splicing.

Strengths<br /> The overall methods and results are robust. The histological analysis of the Srsf1 cKO traces the origins of the fertility defect to the postnatal testis, and the authors have generated interesting datasets characterizing SRSF1's RNA targets and interacting proteins specifically in the testis.

Ultimately, the authors have shown that SRSF1's effects on alternative splicing are required to establish spermatogenesis. In the absence of Srsf1, the postnatal gonocytes do not properly mature into spermatogonia and consequently never initiate spermatogenesis.

Reviewer #2 (Public Review):

The manuscript of Duewell et al has made critical observations that help to understand the mechanisms of activation of the class IA PI3Ks. By using single-molecule kinetic measurements, the authors have made outstanding progress toward understanding how PI3Kbeta is uniquely activated by phosphorylated tyrosine kinase receptors, Gbeta/gamma heterodimers and the small G protein Rac1. While previous studies have defined these as activators of PI3Kbeta, the current manuscript makes clear the quantitative limitations of these previous observations. Most previous quantitative in vitro studies of PI3Kbeta activation have used soluble peptides derived from bis-phosphorylated receptors to stimulate the enzyme. These soluble peptides stimulate the enzyme, and even stimulate membrane interaction. Although these previous studies showed that the release of p85-mediated autoinhibition unmasks an intrinsic affinity of the enzyme for lipid membranes, they ignored what would be the consequence of these peptide sequences being present in the context of intrinsic membrane proteins. The current manuscript shows that the effect of membrane-conjugated peptides on the enzyme activity is profound, in terms of recruiting the enzyme to membranes. In this context, the authors show that G proteins associated with the membranes have an important contribution to membrane recruitment, but they also have a profound allosteric effect on the activity on the membrane, These are observations that would not have been possible with bulk measurements, and they do not simply recapitulate observations that were made for other class IA PI3Ks.

An important observation that the authors have made is that Gbeta/gamma heterodimers and RAc1 alone have almost no ability to recruit PI3Kbeta to the membranes that they are using, and this is central to one of the most profoundly novel activation mechanisms offered by the manuscript. The authors propose that the nSH2- and Gbeta/gamma binding sites partially overlap, so that Gbeta/gamma can only bind once the nSH2 domain releases the p110beta subunit. This mechanism would mean that once the nSH2 is engaged by membrane-congugated pY, the Gbg heterodimer can bind and increase the association of the enzyme with membranes. Indeed, this increased membrane association is observed by the authors. However, the authors also show that this increased recruitment to membranes accounts for relatively little increase in activity, and that the far greater component of activation is due to an allosteric effect of the membrane association on the activity of the enzyme. The proposal for competition between Gbg binding and the nSH2 is consistent with the behavior of an nSH2 mutant that cannot bind to pY and which, consequently, does not vacate the Gbg-binding site. In addition to the outstanding contribution to understanding the kinetics of activation of PI3Kbeta, the authors have offered the first structural interpretation for the kinetics of Gbg activation in synergy with pY activation. The proposal for an overlapping nSH2/Gbg binding site is supported by predictions made by John Burke, using alphafold multimer. Although there is no experimental structure to support this structural model, it is consistent with HDX-MS analyses that were published previously.

Reviewer #2 (Public Review):

In the study by Hreich et al, the potency of P2RX7-specific positive modulator HEI3090, developed by the authors, for the treatment of Idiopathic pulmonary fibrosis (IPF) was investigated. Recently, the authors have shown that HEI3090 can protect against lung cancer by stimulating dendritic cell P2RX7, resulting in IL-18 production that stimulates IFN-γ production by T and NK cells (DOI: 10.1038/s41467-021-20912-2). Interestingly, HEI3090 increases IL-18 levels only in the presence of high eATP. Since the treatment options for IPF are limited, new therapeutic strategies and targets are needed. The authors first show that P2RX7/IL-18/IFNG axis is downregulated in patients with IPF. Next, they used a bleomycin-induced lung fibrosis mouse model to show that the use of a positive modulator of P2RX7 leads to the activation of the P2RX7/IL-18 axis in immune cells that limits lung fibrosis onset or progression. Mechanistically, treatment with HEI3090 enhanced IL-18-dependent IFN-γ production by lung T cells leading to a decreased production of IL-17 and TGFβ, major drivers of IPF. The major novelty is the use of the small molecule HEI3090 to stimulate the immune system to limit lung fibrosis progression by targeting the P2RX7, which could be potentially combined with current therapies available. Overall, the study was well performed and the manuscript is clear.

Reviewer #2 (Public Review):

Summary:<br /> Gaikwad et al. investigated the role of eIF2A in translational response to stress in yeast. For this purpose, the authors conducted ribosome profiling under SM treatment in eIF2A-depleted strain. Data analysis revealed that eIF2A did not influence translation from mRNAs bearing uORFs or cellular IRESes, in the stress condition, broadly. The authors found that only a small number of mRNAs were supported by eIF2A. The data should be helpful for researchers in the fields.

Major points:<br /> 1. The weakness of this work is the lack of clarification on the function of eIF2A in general. The novelty of this study was limited.

2. Related to this, it would be worth investigating common features in mRNAs selectively regulated (surveyed in Figure 3A). Also, it would be worth analyzing the effect of eIF2A deletion on elongation (ribosome occupancy on each codon and/or global ribosome footprint distribution along CDS) and termination/recycling (footprint reads on stop codon and on 3′ UTR).

3. Regarding Figure 3D, the reporters were designed to include promoter and 5′ UTR of the target genes. Thus, it should be worth noting that reporter design was based on the assumption that eIF2A-dependency in translation regulation was not dependent on 3′ UTR or CDS region. The reason why the effects on ribosome profiling-supported mRNAs could not be recapitulated in reporter assay may originate from this design. This should be also discussed.

4. Related to the point above, the authors claimed that eIF2A affects "possibly only one" (HKR1) mRNA. However, this was due to the reporter assay which is technically variable and could not allow some of the constructs to pass the authors' threshold. Authors may be worth considering better wording for this point.

5. For Figure 3D, it would be worth considering to test all the #-marked genes (in Figure 3C) in this set up.

6. In box plots, the authors should provide the statistical tests, at least where the authors explained in the main text.

Reviewer #2 (Public Review):

Summary:

In this manuscript, the goal of the authors is to understand the process of mature sprout formation from mini-sprouts to develop new blood vessels during angiogenesis. For this, they use their earlier experimental setup of engineered blood vessels in combination with a modified spatio-temporal model for Notch signalling. The authors first study the role of VEGF on Tip (Delta-rich) and Stalk (Notch-rich) patterning. The Tip cells are further examined for their space-time dynamics as Mini-sprouts and mature Sprouts. The Notch signalling model is later supplemented with a phenomenological _random uniform model_ for Sprout selection as a plausible mechanism for Sprout formation from Mini-Sprouts. Finally, the authors look into the role of fibronectin in the Sprout formation process. Overall, the authors propose that VEGF interacts with Notch signalling in blood vessels to generate spatially disordered and co-localized Tip cells. VEGF and fibronectin then provide external cues to dynamically modulate mature Sprout formation from Mini-Sprouts that could control the location and density of developing blood vessels with a process that is consistent with a Turing-like mechanism.

Strengths and Weaknesses

In this manuscript, work motivation, problem definition, experimental procedures, analysis techniques, mathematical methods (including the parameters), and findings are all presented quite clearly. Moreover, the authors carefully indicate whenever they make any assumptions, and do not mix unproven hypothesis with deduced or known facts. The experimental techniques and most of the mathematical methods used in this paper are borrowed from the earlier works of the corresponding authors, and thus are not completely novel. However, the use of these ideas to provide a simple elucidation of the role of VEGF and fibronectin in Sprout formation, in an otherwise complex system, is very interesting and useful. Some of the data analysis methods presented in the paper - (i) quantification of Tip spatial patterns (Fig. 3) and (ii) Sprout temporal dynamics using Sankey diagram (Fig. 4) - seem quite novel to me in the context of Notch signalling literature. Similarly, the authors also provide a new mechanism (VEGF) to obtain disordered Delta-Notch patterning without explicitly including _noise_ in the system (Fig. 2 and Fig. S1). The authors also systematically quantify the statistics of spacing between the Sprouts and show that the Sprouts have a tendency to be away from each other, something that they could also partially recapitulate by additionally including a novel _random uniform model_ for Sprout selection (Fig. 5). Although the association between fibronectin and angiogenesis is known in the literature, in this manuscript, the authors could clearly demonstrate that fibronectin is present in high and low levels, respectively, around Sprouts and Mini-sprouts (Fig. 6). A combination of these findings could then motivate the authors to hypothesize, as mentioned above, a Turing-like mechanism for Sprout formation, something that I find interesting.

Although I find the relative simplicity of the experimental system and theoretical model and the clear findings they generate appealing, some aspects raise a few questions. The authors experimentally find 20 +- 0.08 percent of Tip cells in the model blood-vessels that is consistent with the salt-and-pepper pattern seen in Notch signalling model (~25 %). However, it is not clear to me if the reverse is true, i.e., 25% of Tip cells automatically imply a salt-and-pepper pattern - the authors do not seem to provide a direct experimental evidence. Furthermore, the authors use their Notch signalling model on a regular hexagonal lattice, but there is a large variability in the cell sizes (Fig. 3) in the experimental system. Since it is observed in the literature that signalling depends on the contact area between the neighbouring cells, it is not clear how that would affect the findings presented in this paper. Similarly, since some of the cells are quite small compared to the others, I worry how appropriate it is to express the distance between the Tip cells in terms of _cell numbers_ (Fig. 3). Regarding Sprout classification, as per Table 1, a bridge of two cells is formed as per early-stage-I mechanism for Sprout. On the other hand, the entire data interpretation of experiments seems to be based on early Stage II and matured stage in that same table (also Figs. 3 and 4) in which only one Tip cell seems to be counted per mature Sprout. However, if some Sprouts are formed via early stage-I mechanism, a projection in 2D for analysis would give a count of __two__ adjacent Tip cells, but corresponding to a __single__ Sprout. It could be possible that the presence of such two-cell Sprouts affects the statistics of inter-Sprout distances (Fig. 5). Finally, I find the proposed mechanism of Sprout formation dynamics to be somewhat unsatisfactory. Other than the experimental evidence regarding the spacing of Sprouts and the fibronectin levels around Sprouts and Mini-sprouts (Figs. 4 and 5), there is very little evidence to support the hypothesis about a Turing-like mechanism for Sprouting. Moreover, it seems to me that Turing patterns can appear in a wide variety of settings and could be applied to the current problem in an abstract manner without making any meaningful connections with the system variables. Also, from a modeling point of view, cell migration and mechanics, are expected to take a major part in Sprout formation, while cell division and inclusion would most likely influence Tip-Stalk cell formation. However, it seems that in the present work, these effects are coarse-grained into Notch signalling parameters and the Sprout selection model, thus making any experimental connection quite vague.

Overall Assessment

I feel that the authors, on the whole, do achieve their main goals. Although I have a few concerns that I have raised above, overall, I find the work presented in this manuscript to be a solid addition to the broad field of collective cell dynamics. The authors use well established experimental and mathematical methods while adding a few novel analysis techniques and modeling ideas to provide a compelling, albeit incomplete, picture of Sprout formation during angiogenesis. While the direct application of this work in the context of angiogenesis is obvious, the broad set of ideas and techniques (discussed above) in this work would also be useful to researchers who work on Notch signalling in morphogenesis, collective cell migration, and epithelial-mesenchymal-transition.

Reviewer #2 (Public Review):

This is an excellent study. It starts with the identification of two bactofilins in H. neptunium, a demonstration of their important role for the determination of cell shape and discovery of an associated endopeptidase to provide a convincing model for how these two classes of proteins interact to control cell shape. This model is backed up by a quantitative characterisation of their properties using high-resolution imaging and image analysis methods.

Overall, all evidence is very convincing and I do not have many recommendations on how to improve the manuscript.

In my opinion, there are only two issues that I have with the paper:

1. The single particle dynamics of BacA is presented and analysed and I would like to give some suggestions on how to maybe extract even more information from the already acquired data:

1.1. Presentation: Figure 5A is only showing projections of single particle time-lapse movies. To convince the reader that it was indeed possible to detect single molecules it would be helpful if the authors present individual snapshots and intensity traces. In case of single molecules these will show step wise bleaching<br /> 1.2. Analysis: Figure 5B and Supplement Figure 1 are showing the single particle tracking results, revealing that there are two populations of BacA-YFP in the cell. However, this data does not show if individual BacA particles transition between these two populations or not. A more detailed analysis of the existing data, where one can try to identify confinement events in single particle trajectories could be very revealing and help to understand the behaviour of BacA in more detail.

2. The title of Fig. 3 says that BacA and BacD copolymerise, however, the data presented to confirm this conclusions is actually rather weak. First, the Alphafold prediction does not show the co-polymer, and second, the in vitro polymerisation experiments were only done with BacA in the absence of BacD. Accordingly, the only evidence that supports this is their colocalization in fluorescence microscopy. I suggest to either weaken the statement or change the title and add more evidence.

Finally, did the authors think about biochemical experiments to study the interaction between the cytoplasmic part of LmdC and the bactofilins? These could further support their model.

Reviewer #2 (Public Review):

This is the first comprehensive study aimed at assessing the impact of landscape modification on the prevalence of P. knowlesi malaria in non-human primates in Southeast Asia. This is a very important and timely topic both in terms of developing a better understanding of zoonotic disease spillover and the impact of human modification of landscape on disease prevalence.

This study uses the meta-analysis approach to incorporate the existing data sources into a new and completely independent study that answers novel research questions linked to geospatial data analysis. The challenge, however, is that neither the sampling design of previous studies nor their geospatial accuracy are intended for spatially-explicit assessments of landscape impact. On the one hand, the data collection scheme in existing studies was intentionally opportunistic and does not represent a full range of landscape conditions that would allow for inferring the linkages between landscape parameters and P. knowlesi prevalence in NHP across the region as a whole. On the other hand, the absolute majority of existing studies did not have locational precision in reporting results and thus sweeping assumptions about the landscape representation had to be made for the modeling experiment. Finally, the landscape characterization was oversimplified in this study, making it difficult to extract meaningful relationships between the NHP/human intersection on the landscape and the consequences for P. knowlesi malaria transmission and prevalence.

Despite study limitations, the authors point to the critical importance of understanding vector dynamics in fragmented forested landscapes as the likely primary driver in enhanced malaria transmission. This is an important conclusion particularly when taken together with the emerging evidence of substantially different mosquito biting behaviors than previously reported across various geographic regions.

Another important component of this study is its recognition and focus on the value of geospatial analysis and the availability of geospatial data for understanding complex human/environment interactions to enable monitoring and forecasting potential for zoonotic disease spillover into human populations. More multi-disciplinary focus on disease modeling is of crucial importance for current and future goals of eliminating existing and preventing novel disease outbreaks.

Reviewer #2 (Public Review):

Summary:<br /> In eukaryotes, sterols are crucial for signaling and regulating membrane fluidity, however, the mechanism governing cholesterol production and transport across the cell membrane in bacteria remains enigmatic. The manuscript by Zhai et al. sheds light on this topic by uncovering three potential cholesterol transport proteins. Through comprehensive bioinformatics analysis, the authors identified three genes bstA, bstB, and bstC encoding proteins which share homology with transporters, periplasmic binding proteins, and periplasmic components superfamily, respectively. Furthermore, the authors confirmed the specific interaction between these three proteins and C-4 methylated sterols and determined the structures of BstB and BstC. Combining these structural insights with molecular dynamics simulation, they postulated several plausible substrate binding sites within each protein.

Strengths:<br /> The authors have identified 3 proteins that seem likely to be involved in sterol transport between the inner and outer membrane. The structures are of high quality, and the sterol binding experiments support a role for these proteins in sterol transport.

Weaknesses:<br /> While the author's model is very plausible, direct evidence for a role of BstABC in transport, or that the 3 proteins function together in a single pathway, is limited.

Reviewer #2 (Public Review):

Summary:<br /> The authors describe the structure-functional relationship of domains in S. pombe CAF-1, which promotes DNA replication-coupled deposition of histone H3-H4 dimer. The authors nicely showed that the ED domain with an intrinsically disordered structure binds to histone H3-H4, that the KER domain binds to DNA and that, in addition to a PIP box, the KER domain also contributes to the PCNA binding. The ED and KER domains as well as the WHD domain are essential for nucleosome assembly in vitro. The ED, KER domains and the PIP box are important for the maintenance of heterochromatin.

Strengths:<br /> The combination of structural analysis using NMR and Alphafold2 modeling with biophysical and biochemical analysis provided strong evidence on the role of the different domain structures of the large subunit of SpCAF-1, spPCF-1 in the binding to histone H3-H4, DNA as well as PCNA. The conclusion was further supported by genetic analysis of the various pcf1 mutants. The large amounts of data provided in the paper support the authors' conclusion very well.

Weaknesses:

Reviewer #2 (Public Review):

Summary:<br /> This study from Bamgbose et al. identifies a new and important interaction between H4K20me and Parp1 that regulates inducible genes during development and heat stress. The authors present convincing experiments that form a mostly complete manuscript that significantly contributes to our understanding of how Parp1 associates with target genes to regulate their expression.

Strengths:<br /> The authors present 3 compelling experiments to support the interaction between Parp1 and H4K20me, including:

1) PR-Set7 mutants remove all K4K20me and phenocopy Parp mutant developmental arrest and defective heat shock protein induction.

2) PR-Set7 mutants have dramatically reduced Parp1 association with chromatin and reduced poly-ADP ribosylation.

3) Parp1 directly binds H4K20me in vitro.

Weaknesses:<br /> 1) The histone array experiment in Fig1 strongly suggests that PARP binds to all mono-methylated histone residues (including H3K27, which is not discussed). Phosphorylation of nearby residues sometimes blocks this binding (S10 and T11 modifications block binding to K9me1, and S28P blocks binding to K27me1). However, H3S3P did not block H3K4me1, which may be worth highlighting. The H3K9me2/3 "blocking effect" is not nearly as strong as some of these other modifications, yet the authors chose to focus on it. Rather than focusing on subtle effects and the possibility that PARP "reads" a "histone code," the authors should consider focusing on the simple but dramatic observation that PARP binds pretty much all mono-methylated histone residues. This result is interesting because nucleosome mono-methylation is normally found on nucleosomes with high turnover rates (Chory et al. Mol Cell 2019)- which mostly occurs at promoters and highly transcribed genes. The author's binding experiments could help to partially explain this correlation because PARP could both bind mono-methylated nucleosomes and then further promote their turnover and lower methylation state.

2) The RNAseq analysis of Parp1/PR-Set7 mutants is reasonable, but there is a caveat to the author's conclusion (Line 251): "our results indicate H4K20me1 may be required for PARP-1 binding to preferentially repress metabolic genes and activate genes involved in neuron development at co-enriched genes." An alternative possibility is that many of the gene expression changes are indirect consequences of altered development induced by Parp1 or PR-Set7 mutants. For example, Parp1 could activate a transcription factor that represses the metabolic genes that they mention. The authors should consider discussing this possibility.

3) The section on the inducibility of heat shock genes is interesting but missing an important control that might significantly alter the author's conclusions. Hsp23 and Hsp83 (group B genes) are transcribed without heat shock, which likely explains why they have H4K20me without heat shock. The authors made the reasonable hypothesis that this H4K20me would recruit Parp-1 upon heat shock (line 270). However, they observed a decrease of H4K20me upon heat shock, which led them to conclude that "H4K20me may not be necessary for Parp1 binding/activation" (line 275). However, their RNA expression data (Fig4A) argues that both Parp1 and H40K20me are important for activation. An alternative possibility is that group B genes indeed recruit Parp1 (through H4K20me) upon heat shock, but then Parp1 promotes H3/H4 dissociation from group B genes. If Parp1 depletes H4, it will also deplete H4K20me1. To address this possibility, the authors should also do a ChIP for total H4 and plot both the raw signal of H4K20me1 and total H4 as well as the ratio of these signals. The authors could also note that Group A genes may similarly recruit Parp1 and deplete H3/H4 but with different kinetics than Group B genes because their basal state lacks H4K20me/Parp1. To test this possibility, the authors could measure Parp association, H4K20methylation, and H4 depletion at more time points after heat shock at both classes of genes.

Reviewer #2 (Public Review):

Summary:<br /> In contrast to the recent findings reported by Schuster S et al., this brief paper presents evidence suggesting that the stumpy form of T. brucei is likely the most pre-adapted form to progress through the life cycle of this parasite in the tsetse vector.

Strengths:<br /> One significant experimental point is that all fly infection experiments are conducted in the absence of "boosting" metabolites like GlcNAc or S-glutathione. As a result, flies infected with slender trypanosomes present very low or nonexistent infection rates. This provides important experimental evidence that the findings of Schuster S and colleagues may need to be revisited.

Weaknesses:<br /> However, I believe the authors should have included their own set of experiments demonstrating that the presence of these metabolites in the infectious bloodmeal enhances infection rates in flies receiving blood meals containing slender trypanosomes. Considering the well-known physiological variabilities among flies from different facilities, including infection rates, this would have strengthened the experimental evidence presented by the authors.

Reviewer #2 (Public Review):

Summary:<br /> This work provides a new framework, "GPsite" to predict DNA, RNA, peptide, protein, ATP, HEM, and metal ions binding sites on proteins. This framework comes with a webserver and a database of annotations. The core of the model is a Geometric featurizer neural network that predicts the binding sites of a protein. One major contribution of the authors is the fact that they feed this neural network with predicted structure from ESMFold for training and prediction (instead of native structure in similar works) and a high-quality protein Language Model representation. The other major contribution is that it provides the public with a new light framework to predict protein-ligand interactions for a broad range of ligands.

The authors have demonstrated the interest of their framework with mostly two techniques: ablation and benchmark.

Strengths:<br /> The performance of this framework as well as the provided dataset and web server make it useful to conduct studies.

The ablations of some core elements of the method, such as the protein Language Model part, or the input structure are very insightful and can help convince the reader that every part of the framework is necessary. This could also guide further developments in the field. As such, the presentation of this part of the work can hold a more critical place in this work.

Weaknesses:<br /> Overall, we can acknowledge the important effort of the authors to compare their work to other similar frameworks. Yet, the lack of homogeneity of training methods and data from one work to the other makes the comparison slightly unconvincing, as the authors pointed out. Overall, the paper puts significant effort into convincing the reader that the method is beating the state of the art. Maybe, there are other aspects that could be more interesting to insist on (usability, interest in protein engineering, and theoretical works).

Reviewer #2 (Public Review):

Summary:<br /> Glioblastoma is a common primary brain cancer, that is difficult to treat and has a low survival rate. The lack of genetically tractable and immunocompetent vertebrate animal models has prevented the discovery of new therapeutic targets and limited efforts for screening pharmaceutical agents for the treatment of the disease. Here Weiss et al., express oncogenic variants frequently observed in human glioblastoma within zebrafish lacking the tumor suppressor TP53 to generate a patient-relevant in vivo model. The authors demonstrate that loss of TP53 and overexpression of EGFR, PI3KCA, and mScarlet (p53EPS) in neural progenitors and radial glia leads to visible fluorescent brain lesions in live zebrafish. The authors performed RNA expression analysis that uncovered a molecular signature consistent with human mesenchymal glioblastoma and identified gene expression patterns associated with inflammation. Live imaging revealed high levels of immune cell infiltration and associations between microglia/macrophages and tumor cells. To define functional roles for regulators of inflammation on specific immune-related responses during tumorigenesis, transient CRISPR/Cas9 gene targeting was used to disrupt interferon regulator factor proteins and showed Inflammation-associated irf7 and irf8 are required to inhibit p53EPS tumor formation. Further, experiments to deplete the macrophages using clodronate liposomes suggest that macrophages contribute to the suppression of tumor engraftment following transplantation. The authors' conclusions are largely supported by the data and the experiments are thoroughly controlled throughout. Taken together, these results provide new insights into the regulation of glioblastoma initiation and growth by the surrounding microenvironment and provide a novel in vivo platform for the discovery of new molecular mechanisms and testing of therapeutics.

Strengths/Weaknesses:<br /> The authors convincingly show that co-injection of activated human EGFRviii, PI3KCAH1047R, and mScarlet into TP53 null zebrafish promotes the formation of fluorescent brain lesions and glioblastoma-like tumor formation. The authors state that oncogenic MAPK/AKT pathway activation drives this glial-derived tumor formation. It would be important to include a wild-type or uninjected control for the pERK and pAKT staining shown in Fig1 I-K to aid in the interpretation of these results. Likewise, quantification of the pERK and pAKT staining would be useful to demonstrate the increase over WT, and would also serve to facilitate comparison with the similar staining in the KPG model (Supp Fig 2D).

The authors use a transplantation assay to further test the tumorigenic potential of dissociated cells from glial-derived tumors. Listing the percentage of transplants that generate fluorescent tumor would be helpful to fully interpret these data. Additionally, it was not clear based on the description in the results section that the transplantation assay was an "experimental surrogate" to model the relapse potential of the tumor cells. This is first mentioned in the discussion. The authors may consider adding a sentence for clarity earlier in the manuscript as it helps the reader better understand the logic of the assay.

The authors nicely show high levels of immune cell infiltration and associations between microglia/macrophages and tumor cells. However, a quantification of the emergence of macrophages over time in relation to tumor initiation and growth would provide significant support to the observations of tumor suppressive activity of the phagocytes. Along these lines, the inclusion of a statement about when leukocytes emerge during normal development would be informative for those not familiar with the zebrafish model.

From the data provided in Figure 4G and Supp Fig 7b, the authors suggest that "increased p53EPS tumor initiation following Ifr gene knock-down is a consequence of irf7 and irf8 loss-of-function in the TME". Given the importance of the local microenvironment highlighted in this study, spatial information in the form of in situ hybridization to identify the relevant location of the expression change would be important to support this conclusion.

The authors used neutral red staining that labels lysosomal-rich phagocytes to assess enrichment at the early stages of tumor initiation. The images in Figure 3 panel A should be labeled to denote the uninjected controls to aid in the interpretation of the data. In Supplemental Figure 6, the neutral red staining in the irf8 CRISPR-injected larvae looks to be increased, counter to the quantification. Can the authors comment if the image is perhaps not representative?

Reviewer #2 (Public Review):

Summary:<br /> This study tests a plausible and intriguing hypothesis that one cause of the differences in the neural underpinnings of concrete and abstract words is differences in their grounding in the current sensory context. The authors reasoned that, in this case, an abstract word presented with a relevant visual scene would be processed in a more similar way to a concrete word. Typically, abstract and concrete words are tested in isolation. In contrast, this study takes advantage of naturalistic movie stimuli to assess the neural effects of concreteness in both abstract and concrete words (the speech within the film), when the visual context is more or less tied to the word meaning (measured as the similarity between the word co-occurrence-based vector for the spoken word and the average of this vector across all present objects). This novel approach allows a test of the dynamic nature of abstract and concrete word processing, and as such could extend the literature and add a useful perspective accounting for differences in processing these word types.

The critical contrasts needed to test the key hypothesis are not presented or not presented in full within the core text. To test whether abstract processing changes when in a situated context, the situated abstract condition would first need to be compared with the displaced abstract condition as in Supplementary Figure 6. Then to test whether this change makes the result closer to the processing of concrete words, this result should be compared to the concrete result. The correlations shown in Figure 6 in the main text are not focused on the differences in activity between the situated and displaced words or comparing the correlation of these two conditions with the other (concrete/abstract) condition. As such they cannot provide conclusive evidence as to whether the context is changing the processing of concrete/abstract words to be closer to the other condition. Additionally, it should be considered whether any effects reflect the current visual processing only or more general sensory processing.

Overall, the study would benefit from being situated in the literature more, including a) a more general understanding of the areas involved in semantic processing (including areas proposed to be involved across different sensory modalities and for verbal and nonverbal stimuli), and b) other differences between abstract and concrete words and whether they can explain the current findings, including other psycholinguistic variables which could be included in the model and the concept of semantic diversity (Hoffman et al.,). It would also be useful to consider whether difficulty effects (or processing effort) could explain some of the regional differences between abstract and concrete words (e.g., the language areas may simply require more of the same processing not more linguistic processing due to their greater reliance on word co-occurrence). Similarly, the findings are not considered in relation to prior comparisons of abstract and concrete words at the level of specific brain regions.

The authors use multiple methods to provide a post hoc interpretation of the areas identified as more involved in concrete, abstract, or both (at different times) words. These are designed to reduce the interpretation bias and improve interpretation, yet they may not successfully do so. These methods do give some evidence that sensory areas are more involved in concrete word processing. However, they are still open to interpretation bias as it is not clear whether all the evidence is consistent with the hypotheses or if this is the best interpretation of individual regions' involvement. This is because the hypotheses are provided at the level of 'sensory' and 'language' areas without further clarification and areas and terms found are simply interpreted as fitting these definitions. For instance, the right IFG is interpreted as a motor area, and therefore sensory as predicted, and the term 'autobiographical memory' is argued to be interoceptive. Language is associated with the 'both' cluster, not the abstract cluster, when abstract >concrete is expected to engage language more. The areas identified for both vs. abstract>concrete are distinguished in the Discussion through the description as semantic vs. language areas, but it is not clear how these are different or defined. Auditory areas appear to be included in the sensory prediction at times and not at others. When they are excluded, the rationale for this is not given. Overall, it is not clear whether all these areas and terms are expected and support the hypotheses. It should be possible to specify specific sensory areas where concrete and abstract words are predicted to be different based on a) prior comparisons and/or b) the known locations of sensory areas. Similarly, language or semantic areas could be identified using masks from NeuroSynth or traditional meta-analyses. A language network is presented in Supplementary Figure 7 but not interpreted, and its source is not given. Alternatively, there could be a greater interpretation of different possible explanations of the regions found with a more comprehensive assessment of the literature. The function of individual regions and the explanation of why many of these areas are interpreted as sensory or language areas are only considered in the Discussion when it could inform whether the hypotheses have been evidenced in the results section.

Additionally, these methods attempt to interpret all the clusters found for each contrast in the same way when they may have different roles (e.g., relate to different senses). This is a particular issue for the peaks and valleys method which assesses whether a significantly larger number of clusters is associated with each sensory term for the abstract, concrete, or both conditions than the other conditions. The number of clusters does not seem to be the right measure to compare. Clusters differ in size so the number of clusters does not represent the area within the brain well. Nor is it clear that many brain regions should respond to each sensory term, and not just one per term (whether that is V1 or the entire occipital lobe, for instance). The number of clusters is therefore somewhat arbitrary. This is further complicated by the assessment across 20 time points and the inclusion of the 'both' category. It would seem more appropriate to see whether each abstract and concrete cluster could be associated with each different sensory term and then summarise these findings rather than assess the number of abstract or concrete clusters found for each independent sensory term. In general, the rationale for the methods used should be provided (including the peak and valley method instead of other possible options e.g., linear regression).

The measure of contextual situatedness (how related a spoken word is to the average of the visually presented objects in a scene) is an interesting approach that allows parametric variation within naturalistic stimuli, which is a potential strength of the study. This measure appears to vary little between objects that are present (e.g., animal or room), and those that are strongly (e.g., monitor) or weakly related (e.g., science). Additional information validating this measure may be useful, as would consideration of the range of values and whether the split between situated (c > 0.6) and displaced words (c < 0.4) is sufficient.

Finally, the study assessed the relation of spoken concrete or abstract words to brain activity at different time points. The visual scene was always assessed using the 2 seconds before the word, while the neural effects of the word were assessed every second after the presentation for 20 seconds. This could be a strength of the study, however almost no temporal information was provided. The clusters shown have different timings, but this information is not presented in any way. Giving more temporal information in the results could help to both validate this approach and show when these areas are involved in abstract or concrete word processing. Additionally, no rationale was given for this long timeframe which is far greater than the time needed to process the word, and long after the presence of the visual context assessed (and therefore ignores the present visual context).

Reviewer #2 (Public Review):

This is a valuable contribution that will facilitate brain transcriptomic analyses and the joint analyses of gene expression and structural and functional imaging. The methods used are solid, and the authors conducted a wide range of analyses to demonstrate the value of the dense gene expression data.

Reviewer #2 (Public Review):

This is a very interesting paper about the coupling of Slack and Nav1.6 and the insight this brings to the effects of quinidine to treat some epilepsy syndromes.

Slack is a sodium-activated potassium channel that is important to hyperpolarization of neurons after an action potential. Slack is encoded by KNCT1 which has mutations in some epilepsy syndromes. These types of epilepsy are treated with quinidine but this is an atypical antiseizure drug, not used for other types of epilepsy. For sufficient sodium to activate Slack, Slack needs to be close to a channel that allows robust sodium entry, like Nav channels or AMPA receptors. but more mechanistic information is not available. Of particular interest to the authors is what allows quinidine to be effective in reducing Slack.

In the manuscript, the authors show that Nav, not AMPA receptors, are responsible for Slack's sensitization to quinidine blockade, at least in cultured neurons (HeK293, primary cortical neurons). Most of the paper focuses on the evidence that Nav1.6 promotes Slack sensitivity to quinidine.

Reviewer #2 (Public Review):

Although the trans-synaptic tracing method mediated by the rabies virus (RV) has been widely utilized to infer input connectivity across the brain to a genetically defined population in mice, the analysis of labeled pre-synaptic neurons in terms of cell-type has been primarily reliant on classical low-throughput histochemical techniques. In this study, the authors made a significant advance toward high-throughput transcriptomic (TC) cell typing by both dissociated single-cell RNAseq and the spatial TC method known as BARseq to decode a vast array of molecularly labeled ("barcoded") RV vector library. First, they demonstrated that a barcoded RV vector can be employed as a simple retrograde tracer akin to AAVretro. Second, they provided a theoretical classification of neural networks at the single-cell resolution that can be attained through barcoded-RV and concluded that the identification of the vast majority (ideally 100%) of starter cells (the origin of RV-based trans-synaptic tracing) is essential for the inference of single-cell resolution neural connectivity. Taking this into consideration, the authors opted for the BARseq-based spatial TC that could, in principle, capture all the starter cells. Finally, they demonstrated the proof-of-concept in the somatosensory cortex, including infrared connectivity from 381 putative pre-synaptic partners to 31 uniquely barcoded-starter cells, as well as many insightful estimations of input convergence at the cell-type resolution in vivo. Collectively, this work will establish a cornerstone for future advancements in rabies-barcode technology.

This revised version incorporates imaging data to assess the stringency of identifying the starter cells in comparison with conventional protein-based detection methods. Additionally, it encompasses insightful discussions concerning potential limitations and offers perspectives on future improvements. The method section is systematically subdivided with subsection numbers, facilitating the cross-referencing of the corresponding sections in the main text and figure legends. I posit that adopting this stylistic approach as the standard for manuscripts delineating innovative methodological strides would be prudent. The clarity of the figure legends has been significantly enhanced, contributing to a more accessible understanding of the figure panels. In sum, this manuscript is articulate and thorough, epitomizing scientific rigor.

Reviewer #3 (Public Review):

The present study presents a comprehensive exploration of the distinct impacts of Isoflurane and Ketamine on c-Fos expression throughout the brain. To understand the varying responses across individual brain regions to each anesthetic, the researchers employ principal component analysis (PCA) and c-Fos-based functional network analysis. The methodology employed in this research is both methodical and expansive. Notably, the utilization of a custom software package to align and analyze brain images for c-Fos positive cells stands out as an impressive addition to their approach. This innovative technique enables effective quantification of neural activity and enhances our understanding of how anesthetic drugs influence brain networks as a whole.

The primary novelty of this paper lies in the comparative analysis of two anesthetics, Ketamine and Isoflurane, and their respective impacts on brain-wide c-Fos expression. The study reveals the distinct pathways through which these anesthetics induce loss of consciousness. Ketamine primarily influences the cerebral cortex, while Isoflurane targets subcortical brain regions. This finding highlights the differing mechanisms of action employed by these two anesthetics-a top-down approach for Ketamine and a bottom-up mechanism for Isoflurane. Furthermore, this study uncovers commonly activated brain regions under both anesthetics, advancing our knowledge about the mechanisms underlying general anesthesia.

Reviewer #2 (Public Review):

Summary<br /> The authors report an extensive series of neuroimaging experiments (at both 3T and 7T) to provide evidence for a scene-selective visual area in the human posterior parietal cortex (PIGS) that is distinct from the main three (parahippocampal place area, PPA; occipital place area, OPA; medial place area, MPA) typically reported in the literature. Further, they argue that in comparison with the other three, this region may specifically be involved in representing ego-motion in natural contexts. The characterization of this scene-selective region provides a useful reference point for studies of scene processing in humans.

Strengths<br /> One of the major strengths of the work is the extensive series of experiments reported, showing clear reproducibility of the main finding and providing functional insight into the region studied. The results are clearly presented and for the most part, convincing.

Weaknesses<br /> One of the major weaknesses of the work is the failure to relate the current results to other findings in the literature, making it hard to assess whether it is is a "previously undescribed scene-selective site".

First, the scene-selective region identified appears to overlap with regions that have previously been identified in terms of their retinotopic properties. In particular, it is unclear whether this region overlaps with V7/IPS0 and/or IPS1. This is particularly important since prior work has shown that OPA often overlaps with v7/IPS0 (Silson et al, 2016, Journal of Vision). The findings would be much stronger if the authors could show how the location of PIGS relates to retinotopic areas (other than V6, which they do currently consider). I wonder if the authors have retinotopic mapping data for any of the participants included in this study. If not, the authors could always show atlas-based definitions of these areas (e.g. Wang et al, 2015, Cerebral Cortex).

Second, recent studies have reported a region anterior to OPA that seems to be involved in scene memory (Steel et al, 2021, Nature Communications; Steel et al, 2023, The Journal of Neuroscience; Steel et al, 2023, biorXiv). Is this region distinct from PIGS? Based on the figures in those papers, the scene memory-related region is inferior to V7/IPS0, so characterizing the location of PIGS to V7/IPS0 as suggested above would be very helpful here as well.

If PIGS overlaps with either of V7/IPS0 or the scene memory-related area described by Steel and colleagues, then arguably it is not a newly defined region (although the characterization provided here still provides new information).

Another reason that it would be helpful to relate PIGS to this scene memory area is that this scene memory area has been shown to have activity related to the amount of visuospatial context (Steel et al, 2023, The Journal of Neuroscience). The conditions used to show the sensitivity of PIGS to ego-motion also differ in the visuospatial context that can be accessed from the stimuli. Even if PIGS appears distinct from the scene memory area, the degree of visuospatial context is an alternative account of what might be represented in PIGS.

Reviewer #2 (Public Review):

Summary:<br /> Tian et al. explore the developmental organs of cortical reorganization in blindness. Previous work has found that a set of regions in the occipital cortex show different functional responses and patterns of functional correlations in blind vs. sighted adults. In this paper, Tian et al. ask: how does this organization arise over development? Is the "starting state" more like the blind pattern, or more like the adult pattern? Their analyses reveal that the answer depends on the particular networks investigated; some functional connections in infants look more like blind than sighted adults; other functional connections look more like sighted than blind adults; and others fall somewhere in the middle, or show an altogether different pattern in infants compared with both sighted and blind adults.

Strengths:<br /> The question raised in this paper is extremely important: what is the starting state in development for visual cortical regions, and how is this organization shaped by experience? This paper is among the first to examine this question, particularly by comparing infants not only with sighted adults but also blind adults, which sheds new light on the role of visual (and cross-modal) experience. Another clear strength lies in the unequivocal nature of many results. Many results have very large effect sizes, critical interactions between regions and groups are tested and found, and infant analyses are replicated in split halves of the data.

Weaknesses:<br /> A central claim is that "infant secondary visual cortices functionally resemble those of blind more than sighted adults" (abstract, last paragraph of intro). I see two potential issues with this claim. First, a minor change: given the approaches used here, no claims should be made about the "function" of these regions, but rather their "functional correlations". Second (and more importantly), the claim that the secondary visual cortex in general resembles blind more than sighted adults is still not fully supported by the data. In fact, this claim is only true for one aspect of secondary visual area functional correlations (i.e., their connectivity to A1/M1/S1 vs. PFC). In other analyses, the infant secondary visual cortex looks more like sighted adults than blind adults (i.e., in within vs. across hemisphere correlations), or shows a different pattern from both sighted and blind adults (i.e., in occipito-frontal subregion functional connectivity). It is not clear from the manuscript why the comparison to PFC vs. non-visual sensory cortex is more theoretically important than hemispheric changes or within-PFC correlations (in fact, if anything, the within-PFC correlations strike me as the most important for understanding the development and reorganization of these secondary visual regions). It seems then that a more accurate conclusion is that the secondary visual cortex shows a mix of instructive effects of vision and reorganizing effects of blindness, albeit to a different extent than the primary visual cortex.

Relatedly, group differences in overall secondary visual cortex connectivity are particularly striking as visualized in the connectivity matrices shown in Figure S1. In the results (lines 105-112), it is noted that while the infant FC matrix is strongly correlated with both adult groups, the infant group is nonetheless more strongly correlated with the blind than sighted adults. I am concerned that these results might be at least partially explained by distance (i.e., local spread of the bold signal), since a huge portion of the variance in these FC matrices is driven by stronger correlations between regions within the same system (e.g., secondary-secondary visual cortex, frontal-frontal cortex), which are inherently closer together, relative to those between different systems (e.g., visual to frontal cortex). How do results change if only comparisons between secondary visual regions and non-visual regions are included (i.e., just the pairs of regions within the bold black rectangle on the figure), which limits the analysis to long-rang connections only? Indeed, looking at the off-diagonal comparisons, it seems that in fact there are three altogether different patterns here in the three groups. Even if the correlation between the infant pattern and blind adult pattern survives, it might be more accurate to claim that infants are different from both adult groups, suggesting both instructive effects of vision and reorganizing effects of blindness. It might help to show the correlation between each group and itself (across independent sets of subjects) to better contextualize the relative strength of correlations between the groups.

It is not clear that differences between groups should be attributed to visual experience only. For example, despite the title of the paper, the authors note elsewhere that cross-modal experience might also drive changes between groups. Another factor, which I do not see discussed, is possible ongoing experience-independent maturation. The infants scanned are extremely young, only 2 weeks old. Although no effects of age are detected, it is possible that cortex is still undergoing experience-independent maturation at this very early stage of development. For example, consider Figure 2; perhaps V1 connectivity is not established at 2 weeks, but eventually achieves the adult pattern later in infancy or childhood. Further, consider the possibility that this same developmental progression would be found in infants and children born blind. In that case, the blind adult pattern may depend on blindness-related experience only (which may or may not reflect "visual" experience per se). To deal with these issues, the authors should add a discussion of the role of maturation vs. experience and temper claims about the role of visual experience specifically (particularly in the title).

The authors measure functional correlations in three very different groups of participants and find three different patterns of functional correlations. Although these three groups differ in critical, theoretically interesting ways (i.e., in age and visual/cross-modal experience), they also differ in many uninteresting ways, including at least the following: sampling rate (TR), scan duration, multi-band acceleration, denoising procedures (CompCor vs. ICA), head motion, ROI registration accuracy, and wakefulness (I assume the infants are asleep).

Addressing all of these issues is beyond the scope of this paper, but I do feel the authors should acknowledge these confounds and discuss the extent to which they are likely (or not) to explain their results. The authors would strengthen their conclusions with analyses directly comparing data quality between groups (e.g., measures of head motion and split-half reliability would be particularly effective).

Reviewer #2 (Public Review):

Summary:

This study investigates the neural substrates of syntax variation in Bengalese finch songs. Here, the authors tested the effects of bilateral lesions of mMAN, a brain area with inputs to HVC, a premotor area required for song production. Lesions in mMAN induce variability in syntactic elements of song specifically through increased transition entropy, variability within stereotyped song elements known as chunks, and increases in the repeat number of individual syllables. These results suggest that mMAN projections to HVC contribute to multiple aspects of song syntax in the Bengalese finch. Overall the experiments are well-designed, the analysis excellent, and the results are of high interest.

Strengths:

The study identifies a novel role for mMAN, the medial magnocellular nucleus of the anterior nidopallium, in the control of syntactic variation within adult Bengalese finch song. This is of particular interest as multiple studies previously demonstrated that mMAN lesions do not affect song structure in zebra finches. The study undertakes a thorough analysis to characterise specific aspects of variability within the song of lesioned animals. The conclusions are well supported by the data.

Weaknesses:

The study would benefit from additional mechanistic information. A more fine-grained or reversible manipulation, such as brain cooling, might allow additional insights into how mMAN influences specific aspects of syntax structure. Are repeat number increases and transition entropy resulting from shared mechanisms within mMAN, or perhaps arising from differential output to downstream pathways (i.e. projections to HVC)? Similarly, unilateral manipulations would allow the authors to further test the hypothesis that mMAN is involved in inter-hemispheric synchronization.

Reviewer #2 (Public Review):

This paper seeks to determine whether the human visual system's sensitivity to causal interactions is tuned to specific parameters of a causal launching event, using visual adaptation methods. The three parameters the authors investigate in this paper are the direction of motion in the event, the speed of the objects in the event, and the surface features or identity of the objects in the event (in particular, having two objects of different colors).

The key method, visual adaptation to causal launching, has now been demonstrated by at least three separate groups and seems to be a robust phenomenon. Adaptation is a strong indicator of a visual process that is tuned to a specific feature of the environment, in this case launching interactions. Whereas other studies have focused on retinotopically-specific adaptation (i.e., whether the adaptation effect is restricted to the same test location on the retina as the adaptation stream was presented to), this one focuses on feature-specificity.

The first experiment replicates the adaptation effect for launching events as well as the lack of adaptation event for a minimally different non-causal 'slip' event. However, it also finds that the adaptation effect does not work for launching events that do not have a direction of motion more than 30 degrees from the direction of the test event. The interpretation is that the system that is being adapted is sensitive to the direction of this event, which is an interesting and somewhat puzzling result given the methods used in previous studies, which have used random directions of motion for both adaptation and test events.

The obvious interpretation would be that past studies have simply adapted to launching in every direction, but that in itself says something about the nature of this direction-specificity: it is not working through opposed detectors. For example, in something like the waterfall illusion adaptation effect, where extended exposure to downward motion leads to illusory upward motion on neutral-motion stimuli, the effect simply doesn't work if motion in two opposed directions is shown (i.e., you don't see illusory motion in both directions, you just see nothing). The fact that adaptation to launching in multiple directions doesn't seem to cancel out the adaptation effect in past work raises interesting questions about how directionality is being coded in the underlying process. In addition, one limitation of the current method is that it's not clear whether the motion-direction-specificity is also itself retinotopically-specific, that is, if one retinotopic location were adapted to launching in one direction and a different retinotopic location adapted to launching in the opposite direction, would each test location show the adaptation effect only for events in the direction presented at that location?

The second experiment tests whether the adaptation effect is similarly sensitive to differences in speed. The short answer is no; adaptation events at one speed affect test events at another. Furthermore, this is not surprising given that Kominsky & Scholl (2020) showed adaptation transfer between events with differences in speeds of the individual objects in the event (whereas all events in this experiment used symmetrical speeds). This experiment is still novel and it establishes that the speed-insensitivity of these adaptation effects is fairly general, but I would certainly have been surprised if it had turned out any other way.

The third experiment tests color (as a marker of object identity), and pits it against motion direction. The results demonstrate that adaptation to red-launching-green generates an adaptation effect for green-launching-red, provided they are moving in roughly the same direction, which provides a nice internal replication of Experiment 1 in addition to showing that the adaptation effect is not sensitive to object identity. This result forms an interesting contrast with the infant causal perception literature. Multiple papers (starting with Leslie & Keeble, 1987) have found that 6-8-month-old infants are sensitive to reversals in causal roles exactly like the ones used in this experiment. The success of adaptation transfer suggests, very clearly, that this sensitivity is not based only on perceptual processing, or at least not on the same processing that we access with this adaptation procedure. It implies that infants may be going beyond the underlying perceptual processes and inferring genuine causal content. This is also not the first time the adaptation paradigm has diverged from infant findings: Kominsky & Scholl (2020) found a divergence with the object speed differences as well, as infants categorize these events based on whether the speed ratio (agent:patient) is physically plausible (Kominsky et al., 2017), while the adaptation effect transfers from physically implausible events to physically plausible ones. This only goes to show that these adaptation effects don't exhaustively capture the mechanisms of early-emerging causal event representation.

One overarching point about the analyses to take into consideration: The authors use a Bayesian psychometric curve-fitting approach to estimate a point of subjective equality (PSE) in different blocks for each individual participant based on a model with strong priors about the shape of the function and its asymptotic endpoints, and this PSE is the primary DV across all of the studies. As discussed in Kominsky & Scholl (2020), this approach has certain limitations, notably that it can generate nonsensical PSEs when confronted with relatively extreme response patterns. The authors mentioned that this happened once in Experiment 3 and that a participant had to be replaced. An alternate approach is simply to measure the proportion of 'pass' reports overall to determine if there is an adaptation effect. I don't think this alternate analysis strategy would greatly change the results of this particular experiment, but it is robust against this kind of self-selection for effects that fit in the bounds specified by the model, and may therefore be worth including in a supplemental section or as part of the repository to better capture the individual variability in this effect.

In general, this paper adds further evidence for something like a 'launching' detector in the visual system, but beyond that, it specifies some interesting questions for future work about how exactly such a detector might function.

Kominsky, J. F., & Scholl, B. J. (2020). Retinotopic adaptation reveals distinct categories of causal perception. Cognition, 203, 104339. https://doi.org/10.1016/j.cognition.2020.104339

Kominsky, J. F., Strickland, B., Wertz, A. E., Elsner, C., Wynn, K., & Keil, F. C. (2017). Categories and Constraints in Causal Perception. Psychological Science, 28(11), 1649-1662. https://doi.org/10.1177/0956797617719930

Leslie, A. M., & Keeble, S. (1987). Do six-month-old infants perceive causality? Cognition, 25(3), 265-288. https://doi.org/10.1016/S0010-0277(87)80006-9

Reviewer #2 (Public Review):

In this work, the authors set out to identify time-varying subspaces in the premotor cortical activity of monkeys as they executed/observed a reach-grasp-hold movement of 4 different objects. Then, they projected the neural activity to these subspaces and found evidence of shifting subspaces in the time course of a trial in both conditions, executing and observing. These shifting subspaces appear to be distinct in execution and observation trials. However, correlation analysis of neural dynamics reveals the similarity of dynamics in these distinct subspaces. Taken together, Zhao and Schieber speculate that the condition-dependent activity studied here provides a representation of movement that relies on the actor.<br /> This work addresses an interesting question. The authors developed a novel approach to identify instantaneous subspaces and decoded the object type from the projected neural dynamics within these subspaces. As interesting as these results might be, I have a few suggestions and questions to improve the manuscript:<br /> 1- Repeating the analyses in the paper, e.g., in Fig5, using non-MN units only or the entire population, and demonstrating that the results are specific to MNs would make the whole study much more compelling.<br /> 2- The method presented here is similar and perhaps related to principal angles (https://doi.org/10.2307/2005662). It would be interesting to confirm these results with principal angles. For instance, instead of using the decoding performance as a proxy for shifting subspaces, principal angles could directly quantify the 'shift' (similar to Gallego et al, Nat Comm, 2018). Relatedly, why the decoding of the 'object type' is used to establish the progressive shifting of the subspaces? I would be interested to see the authors' argument. The object type should be much more decodable during movement or hold, than instruction, which is probably why the chance-level decoding performance (horizontal lines) is twice the instruction segment for the movement segment.<br /> 3- Why aren't execution and observation subspaces compared together directly? Especially given that there are both types of trials in the same session with the same recorded population of neurons. Using instantaneous subspaces, or the principal angles between manifolds during exec trials vs obs trials.<br /> 4- The definition of the instantaneous subspaces is a critical point in the manuscript. I think it is slightly unclear: based on the Methods section #715-722 and the main text #173-#181, I gather that the subspaces are based on trial averaged neural activity for each of the 4 objects, separately. So for each object and per timepoint, a vector of size (1, n) -n neurons- is reduced to a vector of (1, 2 or 3 -the main text says 2, methods say 3-) which would be a single point in the low-d space. Is this description accurate? This should be clarified in the manuscript.<br /> 5- Isn't the process of projecting segments of neural dynamics and comparing the results equivalent to comparing the projection matrices in the first place? If so, that might have been a more intuitive avenue to follow.<br /> 6- Lines #385-#389: This process seems unnecessarily complicated. Also, given the number of trials available, this sometimes doesn't make sense. E.g. Monkey R exec has only 8 trials of one of the objects, so bootstrapping 20 trials 500 times would be spurious. Why not, as per Gallego et al, Nat Neurosci 2020 and Safaie et al, Nat 2023 which are cited, concatenate the trials?<br /> 7- Related to the CCA analysis, what behavioural epoch has been used here, the same as the previous analyses, i.e. 100ms? how many datapoint is that in time? Given that CCA is essentially a correlation value, too few datapoints make it rather meaningless. If that's the case, I encourage using, let's say, one window combined of I and G until movement, and one window of movement and hold, such that they are both easier to interpret. Indeed low values of exec-exec in CC2 compared to Gallego et al, Nat Neurosci, 2020 might be a sign of a methodological error.

Reviewer #2 (Public Review):

In this paper, the authors presented a compelling rationale for investigating the role of UBCs in prolonging and diversifying signals. Based on the two types of UBCs known as ON and OFF UBC subtypes, they have highlighted the existing gaps in understanding UBCs connectivity and the need to investigate whether UBCs target UBCs of the same subtype, different subtypes, or both. The importance of this knowledge is for understanding how sensory signals are extended and diversified in the granule cell layer.

The authors designed very interesting approaches to study UBCs connectivity by utilizing transgenic mice expressing GFP and RFP in UBCs, Brainbow approach, immunohistochemical and electrophysiological analysis, and computational models to understand how the feed-forward circuits of interconnected UBCs transform their inputs.

This study provided evidence for the existence of distinct ON and OFF UBC subtypes based on their electrophysiological properties, anatomical characteristics, and expression patterns of mGluR1 and calretinin in the cerebellum. The findings support the classification of GRP UBCs as ON UBCs and P079 UBCs as OFF UBCs and suggest the presence of synaptic connections between the ON and OFF UBC subtypes. In addition, they found that GRP and P079 UBCs form parallel and convergent pathways and have different membrane capacitance and excitability. Furthermore, they showed that UBCs of the same subtype provide input to one another and modify the input to granule cells, which could provide a circuit mechanism to diversify and extend the pattern of spiking produced by mossy fiber input. Accordingly, they suggested that these transformations could provide a circuit mechanism for maintaining a sensory representation of movement for seconds.

Overall, the article is well written in a sound detailed format, very interesting with excellent discovery and suggested model.

Reviewer #2 (Public Review):

The focus of this manuscript was to investigate whether Kv1.8 channels, which have previously been suggested to be expressed in type I hair cells of the mammalian vestibular system, are responsible for the potassium conductance gK,L. This is an important study because gK,L is known to be crucial for the function of type I hair cells, but the channel identity has been a matter of debate for the past 20 years. The authors have addressed this research topic by primarily investigating the electrophysiological properties of the vestibular hair cells from Kv1.8 knockout mice. Interestingly, gK,L was completely abolished in Kv1.8-deficient mice, in agreement with the hypothesis put forward by the authors based on the literature. The surprising observation was that in the absence of Kv1.8 potassium channels, the outward potassium current in type II hair cells was also largely reduced. Type II hair cells express the largely inactivating potassium conductance g,K,A, but not gK,L. The authors concluded that heteromultimerization of non-inactivating Kv1.8 and the inactivating Kv1.4 subunits could be responsible for the inactivating gK,A. Overall, the manuscript is very well written and most of the conclusions are supported by the experimental work. The figures are well described, and the statistical analysis is robust.

My only comment relates to the statement regarding the results providing "evidence" that Kv1.4 form heteromultimers with Kv1.8 channels (see Discussion). The only data I can see from the results is that Kv1.4 channels are expressed in the membrane of type II hair cells, which is not sufficient evidence for the above claim. Is the distribution of Kv1.8 and Kv1.4 overlapping in type II hair cells? Have the authors attempted to perform some pharmacological studies on Kv1.4? For example, would gK,A be completely blocked by a Kv1.4 antagonist? Addressing at least some of these questions would strengthen your argument.

Reviewer #2 (Public Review):

Summary:<br /> This manuscript provides a detailed analysis of B-cell lymphomagenesis in mice lacking an alternative exon in the region encoding the C-terminal (regulatory) domain of the p53 protein and thus enable to assemble the so-called p53AS isoform. This isoform differs from canonical p53 by the replacement of roughly 30 c-terminal residues by about 10 residues encoded by the alternative exon. There is biochemical and biological evidence that p53AS retains strong transcriptional and somewhat enhanced suppressive activities, with mouse models expressing protein constructs similar to p53AS showing signs of increased p53 activity leading to rapid and lethal anemia. However, the precise role of the alternative p53AS variant has not been addressed so far in a mouse model aimed at demonstrating whether the lack of this particular p53 isoform (trp53ΔAS/ΔAS mice) may cause a specific pathological phenotype.

Results show that lack of AS expression does not noticeably affect p53 transcriptional activity but reveals a subtle pathogenic phenotype, with trp53ΔAS/ΔAS males, but not females, tending to develop more frequently and earlier B-cell lymphoma than WT. Next, the authors then introduced ΔAS in transgenic Eμ-Myc mice that show accelerated lymphomagenesis. They show that lack of AS caused increased lethality and larger tumor lymph nodes in p53ΔAS Eμ-Myc males compared to their p53WT Eμ-Myc male counterparts, but not in females. Comparative transcriptomics identified a small set of candidate, differentially expressed genes, including Ackr4 (atypical chemokine receptor 4), which was significantly less expressed in the spleens of ΔAS compared to WT controls. Ackr4 encodes a dummy receptor acting as an interceptor for multiple chemokines and thus may negatively regulate a chemokine/cytokine signalling axis involved in lymphomagenesis, which is down-regulated by estrogen signalling. Using in vitro cell models, the authors provide evidence that Ackr4 is a transcriptional target for p53 and that its p53-dependent activation is repressed by 17b-oestradiol. Finally, seeking evidence for a relevance for this gene in human lymphomagenesis, the authors analyse Burkitt lymphoma transcriptomic datasets and show that high ACKR4 expression correlated with better survival in males, but not in females

Strengths:<br /> A convincing demonstration of a subtle, gender-specific pathogenic phenotype associated with the lack of p53AS. The characterization of trp53ΔAS/ΔAS is well described and the data presented are convincing. This represents a significant achievement since, as mentioned, in vivo data establishing the relevance of p53AS isoform remains scarce. Based on this initial observation, the authors provide strong correlative evidence that this particular phenotype is associated by differential expression of Ackr4.

Weaknesses:<br /> The study does not demonstrate how p53AS may specifically and differentially contribute to the regulation of Ackr4, nor whether restoring Ackr4 expression may nullify the observed phenotype.

Reviewer #3 (Public Review):

Summary:<br /> This paper describes transcriptomes from three tardigrade species with or without treatment with ionizing radiation (IR). The authors show that IR produces numerous single-strand and double-strand breaks as expected and that these are substantially repaired within 4-8 hours. Treatment with IR induces strong upregulation of transcripts from numerous DNA repair proteins including Dsup specific to the Hypsobioidea superfamily. Transcripts from the newly described protein TDR1 with homologs in both Hypsibioidea and Macrobiotoidea supefamilies are also strongly upregulated. They show that TDR1 transcription produces newly translated TDR1 protein, which can bind DNA and co-localizes with DNA in the nucleus. At higher concentrations, TDR appears to form aggregates with DNA, which might be relevant to a possible function in DNA damage repair. When introduced into human U2OS cells treated with bleomycin, TDR1 reduces the number of double-strand breaks as detected by gamma H2A spots. This paper will be of interest to the DNA repair field and to radiobiologists.

Strengths:<br /> The paper is well-written and provides solid evidence of the upregulation of DNA repair enzymes after irradiation of tardigrades, as well as upregulation of the TRD1 protein. The reduction of gamma-H2A.X spots in U2OS cells after expression of TRD1 supports a role in DNA damage.

Weaknesses:<br /> Genetic tools are still being developed in tardigrades, so there is no mutant phenotype to support a DNA repair function for TRD1, but this may be available soon.

Reviewer #2 (Public Review):

Summary:

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

Strengths:

The analyses are sound and the observations are intriguing. Indeed, it has been noted previously that parabolic growth promotes coexistence, its role in mitigating the error threshold catastrophe - which is often presented as a major obstacle to our understanding of the origin of life - had not been examined before.

Weaknesses:

Although all the conclusions are interesting, most are not very surprising for people familiar with the literature. As the authors point out, parabolic growth is well known to promote diversity (Szathmary-Gladkih 89) and it has also been noted previously that a form of Darwinian selection can be found at small population sizes (Davis 2000). Given that under parabolic growth, no sequence is ever excluded for infinite populations, it is also not surprising to find that mutations have a less dramatic exclusionary impact.

A general weakness is the presentation of models and parameters, whose choices often appear arbitrary. Modeling choices that would deserve to be further discussed include the association of the monomers with the strands and the ensuing polymerization, which are combined into a single association/polymerization reaction (see also below), or the choice to restrict to oligomers of length L = 10. Other models, similar to the one employed here, have been proposed that do not make these assumptions, e.g. Rosenberger et al. Self-Assembly of Informational Polymers by Templated Ligation, PRX 2021. To understand how such assumptions affect the results, it would be helpful to present the model from the perspective of existing models.

The values of the (many) parameters, often very specific, also very often lack justifications. For example, why is the "predefined error factor" ε = 0.2 and not lower or higher? How would that affect the results? Similarly, in equation (11), where does the factor 0.8 come from? Why is the kinetic constant for duplex decay reaction 1.15e10−8? Are those values related to experiments, or are they chosen because specific behaviors can happen only then?

The choice of the model and parameters potentially impact the two main results, the attenuation of the error threshold and the role of GC content:

Regarding the error threshold, it is also noted (lines 379-385) that it disappears when back mutations are taken into account. This suggests that overcoming the error threshold might not be as difficult as suggested, and can be achieved in several ways, which calls into question the importance of the particular role of parabolic growth. Besides, when the concentration of replicators is low, product inhibition may be negligible, such that a "parabolic replicator" is effectively growing exponentially and an error catastrophe may occur. Do the authors think that this consideration could affect their conclusion? Can simulations be performed?

Regarding the role of the GC content, GC-rich oligomers are found to perform the worst but no rationale is provided. One may assume that it happens because GC-rich sequences are comparatively longer to release the product. However, it is also conceivable that higher GC content may help in the polymerization of the monomers as the monomers attach longer on the template (as described in Eq.(9)). This is an instance where the choice to pull into a single step the association and polymerization reactions are pulled into a single step independent of GC content may be critical. It would be important to show that the result arises from the actual physics and not from this modeling choice.

Some more specific points that would deserve to be addressed:

- Line 53: it is said that p "reflects how easily the template-reaction product complex dissociates". This statement is not correct. A reaction order p<1 reflects product inhibition, the propensity of templates to bind to each other, not slow product release. Product release can be limiting, yet a reaction order of 1 can be achieved if substrate concentrations are sufficiently high relative to oligomer concentrations (von Kiedrowski et al., 1991).

- Population size is a key parameter, and a comparison is made between small (10^3) and large (10^5) populations, but without explaining what determines the scale (small/large relative to what?).

- In the same vein, we might expect size not to be the only important parameter, but also concentration.

- Lines 543-546: if understanding correctly, the quantitative result is that the error threshold rises from 0.1 in the exponential case to 0.196 in the parabolic. Are the authors suggesting that a factor of 2 is a significant difference?

- Figure 3C: this figure shows no statistically significant effect?

- line 542: "phase transition-like species extension (Figure 4B)": such a clear threshold is not apparent.

Reviewer #2 (Public Review):

Summary:<br /> This work describes the structure of Heparan-alpha-glucosaminide N-acetyltransferase (HGSNAT), a lysosomal membrane protein that catalyzes the acetylation reaction of the terminal alpha-D-glucosamine group required for the degradation of heparan sulfate (HS). HS degradation takes place during the degradation of the extracellular matrix, a process required for restructuring tissue architecture, regulation of cellular function, and differentiation. During this process, HS is degraded into monosaccharides and free sulfate in lysosomes.

HGSNAT catalyzes the transfer of the acetyl group from acetyl-CoA to the terminal non-reducing amino group of alpha-D-glucosamine. The molecular mechanism by which this process occurs has not been described so far. One of the main reasons to study the mechanism of HGSNAT is that multiple mutations spanning the entire sequence of the protein, such as nonsense mutations, splice-site variants, and missense mutations lead to dysfunction that causes abnormal accumulation of HS within the lysosomes. This accumulation is a cause of mucopolysaccharidosis IIIC (MPS IIIC), an autosomal recessive neurodegenerative lysosomal storage disorder, for which there are no approved drugs or treatment strategies.

This paper provides a 3.26A structure of HGSNAT, determined by single-particle cryo-EM. The structure reveals that HGSNAT is a dimer in detergent micelles and a density assigned to acetyl-CoA. The authors speculate about the molecular mechanism of the acetylation reaction, map the mutations known to cause MPS IIIC on the structure and speculate about the nature of the HGSNAT disfunction caused by such mutations.

Strengths:<br /> The description of the architecture of HGSNAT is the highlight of the paper since this corresponds to the first description of the structure of a member of the transmembrane acyl transferase (TmAT) superfamily. The high resolution of an HGSNAT bound to acetyl-CoA is an important leap in our understanding of the HGSNAT mechanism. The density map is of high quality, except for the luminal domain. The location of the acetyl-CoA allows speculation about the mechanistic role of multiple residues surrounding this molecule. The authors thoroughly describe the architecture of HGSNAT and map the mutations leading to MPS IIIC. The description of the dimeric interphase is a novel result, and future studies are left to confirm the importance of oligomerization for function.

Weaknesses:<br /> Apart from the cryo-EM structure, the article does not provide any other experimental evidence to support or explain a molecular mechanism. Due to the complete absence of functional assays, mutagenesis analysis, or other structures such as a ternary complex or an acetylated enzyme intermediate, the mechanistic model depicted in Figure 5 should be taken with caution.

The authors discuss that H269 is an essential residue that participates in the acetylation reaction, possibly becoming acetylated during the process. However, there is no solid experimental evidence, e.g. mutagenesis analysis or structural analysis, in this or previous articles, that demonstrates this to be the case.

In the discussion part, the authors mention previous studies in which it was postulated that the catalytic reaction can be described by a random order mechanistic model or a Ping Pong Bi Bi model. However, the authors leave open the question of which of these mechanisms best describes the acetylation reaction. The structure presented here does not provide evidence that could support one mechanism or the other.

Although the authors map the mutations leading to MPS IIIC on the structure and use FoldX software to predict the impact of these mutations on folding and fold stability, there is no experimental evidence to support FoldX's predictions.

Reviewer #2 (Public Review):

Summary:

The manuscript from Bekeova et al. entitled "Acyl-CoA thioesterase-2 facilitates P-oxidation in glycolytic skeletal muscle in a lipid supply dependent manner" examines whether loss of acyl-CoA thioesterase-2 (ACOT2) in the mitochondrial matrix of skeletal muscle alters mitochondrial fatty acid metabolism. The authors generate data demonstrating that under normal chow conditions, loss of ACOT2 increases mitochondrial respiration of long-chain fatty acid, but also increases susceptibility to the build-up of metabolic intermediates. However, during short-term high-fat feeding (7 days), mice with knockout of skeletal muscle ACOT2 had better glucose and insulin tolerance. Interestingly, skeletal muscle ACOT2 knockout mice on chow and high-fat diet utilized more glucose during the active (dark cycle) portion of the day. These data suggest that ACOT2 may be a potential therapeutic target to improve glucose homeostasis.

Strengths:

The use of creatine kinase cre recombinase to specifically target striated muscle localizes the genetic manipulation, thus increasing the rigor of these experiments by limiting potential off-site changes in ACOT2 expression. Also, the assessment of mitochondrial respiration and response to changes in energy change via the creatine kinase clamp technique is a strength. These data provide a measurement of isolated mitochondrial respiration at physiologically relevant concentrations of ATP and ADP, while also allowing for assessment of how these mitochondria respond to changes in free energy (Fisher-Wellman et al. 2018). The indirect calorimetry data provides systemic physiological context to the striated muscle-specific genetic manipulation, while also allowing for the examination of how this change in skeletal muscle ACOT2 impacts systemic responses to different energy challenges. Finally, the extensive metabolomics, transcriptomics, and lipidomics analysis, not only provides a wealth of data but is used to further the authors' investigation of skeletal muscle ACOT2 activity in mitochondrial fatty acid oxidation and glucose homeostasis.

Weaknesses:

Several general confounding factors exist in the experimental design that could potentially impact the interpretation of the observed outcomes. First, all mice were housed at housing temperatures (22C) below the thermoneutral zone, which has been well described by many investigators to result in dramatically increased energy expenditure. Changes in total and resting energy expenditure could alter the skeletal muscle and systemic utilization of lipids, response to high-fat diet, and glucose homeostasis. Second, no dietary control was observed in these experiments. While this did not impact outcomes when the diets were not compared, once the authors began to compare normal chow to high-fat diet, numerous differences in the composition of these diets could impact the outcomes. Third, the extended food withdrawal before the glucose- and insulin tolerance tests puts the mouse in a state of extreme energy stress more akin to starvation than fasting, which can negatively impact outcomes (Ayala et al. 2010, Virtue & Vidal-Puig 2021). Fourth, the use of the Seahorse platform for the assessment of respiration of isolated mitochondria is highly debatable (Schmidt et al. 2021), particularly when the investigators also used high-resolution respirometry specifically designed for the purpose of measuring isolated mitochondrial oxygen consumption. Importantly, the use of the Seahorse platform to assess cellular respiration in this investigation is quite appropriate. Finally, while the authors present data demonstrating that ACOT2 expression is highest in Type I fibers compared to the various Type II fiber types, a large number of the experiments are performed in a muscle that is primarily composed of Type II fibers. The authors briefly acknowledge this limitation. But, is important for the reader to keep this in mind when trying to consider how these findings would translate to humans.

Impact:

The authors have generated data that implicates skeletal muscle mitochondrial coenzyme A handling as a therapeutic target in the improvement of glucose homeostasis. While the exact role of increased tissue lipid burden on insulin action, glucose uptake, and substrate metabolism is still debated, the association between increased tissue lipid and impaired tissue- and systemic glucose handling is very strong. The data herein suggest that ACOT2 represents a pharmaceutical target to improve systemic glucose homeostasis in the population with obesity.

Reviewer #2 (Public Review):

Summary:<br /> In the present study, van Gerwen et al. perform deep phosphoproteomics on muscle from saline or insulin-injected mice from 5 distinct strains fed a chow or HF/HS diet. The authors follow these data by defining a variety of intriguing genetic, dietary or gene-by-diet phosphor-sites which respond to insulin accomplished through application of correlation analyses, linear mixed models and a module-based approach (WGCNA). These findings are supported by validation experiments by intersecting results with a previous profile of insulin-responsive sites (Humphrey et al, 2013) and importantly, mechanistic validation of Pfkfb3 where overexpression in L6 myotubes was sufficient to alter fatty acid-induced impairments in insulin-stimulated glucose uptake. To my knowledge, this resource provides the most comprehensive quantification of muscle phospho-proteins which occur as a result of diet in strains of mice where genetic and dietary effects can be quantifiably attributed in an accurate manner. Utilization of this resource is strongly supported by the analyses provided highlighting the complexity of insulin signaling in muscle, exemplified by contrasts to the "classically-used" C57BL6/J strain. As it stands, I view this exceptional resource as comprehensive with compelling strength of evidence behind the mechanism explored. I raised several comments in the last round of assessment but all of them have now been thoughtfully addressed.

Strengths: Generation of a novel resource to explore genetic and dietary interactions influencing the phospho-proteome in muscle. This is accompanied by elegant application of in silico tools to highlight the utility

Weaknesses: none noted

Reviewer #3 (Public Review):

The study by Jain et al. on recombinant adeno-associated viruses (rAAVs) represents a valuable contribution to the fields of virus genetics and gene therapy. As non-pathogenic vectors, rAAVs have become a popular choice for delivering gene therapies. The authors have previously investigated the effects of all possible single codon substitutions, deletions, and insertions in the AAV2 cap gene on AAV production. In this study, they extend their analysis to the AAV2 rep gene and rep genes in two additional capsid serotypes, establishing a genotype-phenotype landscape that enhances our understanding of Rep protein function and offers potential strategies for improving Rep function in gene therapy applications. The experimental design is rigorous, the analyses well-executed, and the interpretations of the data are convincing. While I have a few suggestions to further refine the study, I believe it is overall an excellent piece of research.

One aspect that may warrant further consideration is the assumption, as mentioned in Figure 2's legend, that synonymous mutations are neutral and can serve as controls for normalizing the production rate. However, Figures S5-6 and Figures S11-12 suggest that synonymous mutations are not necessarily neutral, as their distribution is similar to that of nonsynonymous mutations. Thus, it may be beneficial to more thoroughly examine the potential effects of synonymous mutations on the genotype-phenotype landscape.

Additionally, previous research by Jeff Collar and others has reported that synonymous mutations can affect mRNA levels through mRNA degradation rate. It would be interesting to determine if the 20-bp barcodes located at the 3' end are positioned within the untranslated regions and could thus be employed to quantify the mRNA levels of individual variants. This information could offer insight into another potential mechanism by which single codon mutations impact the production rate of rAAV.

The authors discovered several novel mutations that enhance AAV production yet are absent in natural occurrences. This intriguing finding could benefit from further elaboration, particularly with regard to the distribution of these mutations within the protein structure and the nature of the amino acid transitions involved. It would also be informative if the authors could provide a brief discussion as to why these mutations have not been observed in nature. For instance, could it be that optimal viral fitness necessitates an intermediate production rate rather than an excessively rapid one? Expanding on these points may further enrich the paper and offer valuable insights for readers.

The authors have taken commendable steps to address the concerns I raised in my previous evaluation. They have provided comprehensive clarifications, performed necessary revisions, and expanded upon certain key points in the manuscript.

Reviewer #3 (Public Review):

In the revised version the authors have addressed some of the reviewers' concerns, but, despite the new explanatory paragraph on page 16, the paper remains confusing because as shown in Figure 7 at the end of the Results the PARG KO 293A cells that were analyzed at the beginning of the Results are not true PARG knockouts. The authors stated that they did not rewrite the Results because they wanted to describe the experiments in the order in which they were carried out, but there is no imperative for the experiments to be described in the order in which they were done, and it would be much easier for the uninitiated reader to appreciate the significance of these studies if the true PARG KO cell data were presented at the beginning, as all three of the original reviewers proposed.

While the authors have to some extent clarified the nature of the PARG KO alleles, they have not been able to identify the source of the residual PARG activity in the PARG KO cells, in part because different commercial PARG antibodies give different and conflicting immunoblotting results. Additional sequence characterization of PARG mRNAs expressed in the PARG cKO cells, and also in-depth proteomic analysis of the different PARG bands could provide further insight into the origins and molecular identities of the various PARG proteins expressed from the different KO PARG alleles, and determine which of them might retain catalytic activity.

The authors have made no progress in identifying which are the key PARG substrates required for S phase progression, although they suggest that PARP1 itself may be an important target.

Reviewer #2 (Public Review):

In this manuscript, the authors have meticulously constructed a comprehensive atlas delineating hematopoietic stem/progenitor cell (HSPC) and immune-cell types within the zebrafish kidney, employing single-cell transcriptome profiling analysis. Notably, these cell populations exhibited distinctive responses to viral infection. Intriguingly, the investigation revealed that HSPCs manifest positive reactivities to viral infection, indicating the effective induction of trained immunity in select HSPCs. Furthermore, the study unveiled the capacity for the generation of antigen-stimulated adaptive immunity within the kidney, suggesting a role for the zebrafish kidney as a secondary lymphoid organ. This research elucidates the distinctive features of the fish immune system and underscores the multifaceted biology of the kidney in ancient vertebrates.

Reviewer #2 (Public Review):

Summary:<br /> The aim of this manuscript is to use molecular dynamics (MD) simulations to describe the conformational changes of the neurotransmitter binding site of a nicotinic receptor. The study uses a simplified model including the alpha-delta subunit interface of the extracellular domain of the channel and describes the binding of four agonists to observe conformational changes during the weak-to-strong affinity transition.

Strength:<br /> The 200 ns-long simulations of this model suggest that the agonist rotates about its centre in a 'flip' motion, while loop C 'flops' to restructure the site. The changes appear to be reproduced across simulations and different ligands and are thus a strong point of the study.

Weaknesses:<br /> After carrying out all-atom molecular dynamics, the authors revert to a model of binding using continuum Poisson-Boltzmann, surface area, and vibrational entropy. The motivations for and limitations associated with this approximate model for the thermodynamics of binding, rather than using modern atomistic MD free energy methods (that would fully incorporate configurational sampling of the protein, ligand, and solvent) could be provided. Despite this, the authors report a correlation between their free energy estimates and those inferred from the experiment. This did, however, reveal shortcomings for two of the agonists. The authors mention their trouble getting correlation to experiment for Ebt and Ebx and refer to up to 130% errors in free energy. But this is far worse than a simple proportional error, because -24 Vs -10 kcal/mol is a massive overestimation of free energy, as would be evident if the authors were to instead express results in terms of KD values (which would have an error exceeding a billion fold). The MD analysis could be improved with better measures of convergence, as well as a more careful discussion of free energy maps as a function of identified principal components, as described below. Overall, however, the study has provided useful observations and interpretations of agonist binding that will help understand pentameric ligand-gated ion channel activation.

Main points:<br /> Regarding the choice of model, some further justification of the reduced 2 subunit ECD-only model could be given. On page 5 the authors argue that, because binding free energies are independent of energy changes outside the binding pocket, they could remove the TMD and study only an ECD subunit dimer. While the assumption of distant interactions being small seems somewhat reasonable, provided conformational changes are limited and localised, how do we know the packing of TMD onto the ECD does not alter the ability of the alpha-delta interface to rearrange during weak or strong binding? They further write that "fluctuations observed at the base of the ECD were anticipated because the TMD that offers stability here was absent.". As the TMD-ECD interface is the "gating interface" that is reshaped by agonist binding, surely the TMD-ECD interface structure must affect binding. It seems a little dangerous to completely separate the agonist binding and gating infrastructure, based on some assumption of independence. Given the model was only the alpha and delta subunits and not the pentamer with TMD, I am surprised such a model was stable without some heavy restraints. The authors state that "as a further control we carried out MD simulation of a pentamer docked with ACh and found similar structural changes at the binding pocket compared to the dimer." Is this sufficient proof of the accuracy of the simplified model? How similar was the model itself with and without agonist in terms of overall RMSD and RMSD for the subunit interface and the agonist binding site, as well as the free energy of binding to each model to compare?

Although the authors repeatedly state that they have good convergence with their MD, I believe the analysis could be improved to convince us. On page 8 the authors write that the RMSD of the system converged in under 200 ns of MD. However, I note that the graph is of the entire ECD dimer, not a measure for the local binding site region. An additional RMSD of local binding site would be much more telling. You could have a structural isomerisation in the site and not even notice it in the existing graph. On page 9 the authors write that the RMSF in Figure S2 showed instability mainly in loops C and F around the pocket. Given this flexibility at the alpha-delta interface, this is why collecting those regions into one group for the calculation of RMSD convergence analysis would have been useful. They then state "the final MD configuration (with CCh) was well-aligned with the CCh-bound cryo-EM desensitized structure (7QL6)... further demonstrating that the simulation had converged." That may suggest a change occurred that is in common with the global minimum seen in cryo EM, which is good, but does not prove the MD has "converged". I would also rename Figure S3 accordingly.

The authors draw conclusions about the dominant states and pathways from their PCA component free energy projections that need clarification. It is important first to show data to demonstrate that the two PCA components chosen were dominant and accounted for most of the variance. Then when mapping free energy as a function of those two PCA components, to prove that those maps have sufficient convergence to be able to interpret them. Moreover, if the free energies themselves cannot be used to measure state stability (as seems to be the case), that the limitations are carefully explained. First, was PCA done on all MD trajectories combined to find a common PC1 & PC2, or were they done separately on each simulation? If so, how similar are they? The authors write "the first two principal components (PC-1 and PC-2) that capture the most pronounced C. displacements". How much of the total variance did these two components capture? The authors write the changes mostly concern loop C and loop F, but which data proves this? e.g. A plot of PC1 and PC2 over residue number might help.

The authors map the -kTln rho as a free energy for each simulation as a function of PC1 & PC2. It is important to reveal how well that PC1-2 space was sampled, and how those maps converged over time. The shapes of the maps and the relative depths of the wells look very different for each agonist. If the maps were sampled well and converged, the free energies themselves would tell us the stabilities of each state. Instead, the authors do not even mention this and instead talk about "variance" being the indicator of stability, stating that m3 is most stable in all cases. While I can believe 200ns could not converge a PC1-2 map and that meaningful delta G values might not be obtained from them, the issue of lack of sampling must be dealt with. On page 12 they write "Although the bottom of the well for 3 energy minima from PCA represent the most stable overall conformation of the protein, they do not convey direct information regarding agonist stability or orientation". The reasons why not must be explained; as they should do just that if the two order parameters PC1 and PC2 captured the slowest degrees of freedom for binding and sampling was sufficient. The authors write that "For all agonists and trajectories, m3 had the least variance (was most stable), again supporting convergence by 200 ns." Again the issue of actual free energy values in the maps needs to be dealt with. The probabilities expressed as -kTln rho in kcal/mol might suggest that m2 is the most stable. Instead, the authors base stability only on variance (I guess breadth of the well?), where m3 may be more localised in the chosen PC space, despite apparently having less preference during the MD (not the lowest free energy in the maps).

The motivations and justifications for the use of approximate PBSA energetics instead of atomistic MD free energies should be dealt with in the manuscript, with limitations more clearly discussed. Rather than using modern all-atom MD free energy methods for relative or absolute binding free energies, the author selects clusters from their identified states and does Poisson-Boltzmann estimates (electrostatic, vdW, surface area, vibrational entropy). I do believe the following sentence does not begin to deal with the limitations of that method: "there are limitations with regard to MM-PBSA accurately predicting absolute binding free energies (Genheden & Ryde, 2015; Hou et al., 2011) that depends on the parameterization of the ligand (Oostenbrink et al., 2004)." What are the assumptions and limitations in taking continuum electrostatics (presumably with parameters for dielectric constants and their assignments to regions after discarding solvent), surface area (with its assumptions and limitations), and of course assuming vibration of a normal mode can capture entropy. On page 30, regarding their vibrational entropy estimate, they write that the "entropy term provides insights into the disorder within the system, as well as how this disorder changes during the binding process". It is important that the extent of disorder captured by the vibrational estimate be discussed, as it is not obvious that it has captured entropy involving multiple minima on the system's true 3N-dimensional energy surface, and especially the contribution from solvent disorder in bound Vs dissociated states.

As discussed above, errors in the free energy estimates need to be more faithfully represented, as fractional errors are not meaningful. On page 21 the authors write "The match improved when free energy ratios rather than absolute values were compared." But a ratio of free energies is not a typical or expected measure of error in delta G. They also write "For ACh and CCh, there is good agreement between.Gm1 and GLA and between.Gm3 and GHA. For these agonists, in silico values overestimated experimental ones only by ~8% and ~25%. The agreement was not as good for the other 2 agonists, as calculated values overestimated experimental ones by ~45%(Ebt) and ~130% (Ebt). However, the fractional overestimation was approximately the same for GLA and GHA." See the above comment on how this may misrepresent the error. On page 21 they write, in relation to their large fractional errors, that they "do not know the origin of this factor but speculate that it could be caused by errors in ligand parameterization". However the estimates from the PBSA approach are, by design, only approximate. Both errors in parameterisation (and their likely origin) and the approximate model used, need discussion.

Reviewer #4 (Public Review):

The manuscript entitled "Ebola Virus Sequesters IRF3 in Viral Inclusion Bodies to Evade Host Antiviral Immunity" mainly describes that the function of IBs formed by the viral proteins VP35 and NP in evading host antiviral immunity. They proved that Ebola virus VP35 protein can interact with STING, but not IRF3, to sequester IRF3 into inclusion bodies and thereby inhibit type-I interferon production. This work will be of some interest to readers in the Ebola Virus field, however, the current data do not clearly explain the relationship of VP35 protein and IRF3.

Reviewer #2 (Public Review):

Summary:

Bian et al. calculated Phenotypic Age Acceleration (PhenoAgeAccel) via a linear model regressing Phenotypic Age on chronological age. They examined the associations between PhenoAgeAccel and cancer incidence using 374,463 individuals from the UK Biobank and found that older PhenoAge was consistently related to an increased risk of incident cancer, even among each risk group defined by genetics.

Strengths:

The study is well-designed, and uses a large sample size from the UK biobank.

Weaknesses:

Since the UK biobank has a large sample size, it should have enough power to split the dataset into discovery and validation sets. Why did the authors use 10-fold cross-validation instead of splitting the dataset?

Reviewer #2 (Public Review):

Summary:<br /> Drosophila hematopoiesis has been shown to be governed by a number of signaling pathways such as JAK/STAT and Dpp. This important study shows the role of nutrient sensing and autophagy in determining blood cell differentiation. The authors show that General control non-derepressible 5 (Gcn5), a histone acetyltransferase affects blood cell differentiation. Gcn5 also negatively regulates autophagy through its effector TFEB which directly regulates autophagy genes. The authors also show that mTORC1 modulates Gcn5 levels and through it, TFEB activity thus acting as a fine-tuning mechanism that maintains optimal levels of autophagy.

Strengths:<br /> The main strength of the work lies in the interesting finding that cellular metabolic processes such as autophagy have a direct role in blood cell differentiation and has the potential to be of interest to those working on vertebrate haematopoiesis as well. The report has generated intriguing data, using promoters specific for sub-sections of the lymph gland, that different cellular subsets of the lymph gland contribute differently towards haematopoiesis, but this is not followed up in detail and the final conclusions are derived from a combination of whole lymph gland perturbations as well as those from specific promoters.

Weaknesses:<br /> 1. Gc5 seems to be expressed throughout the lymph gland but modulating it in the subsections does not have the same result. It is very striking that the knockdown of Gcn5 in the prohemocyte population does not have an effect on differentiation whereas overexpression does. The modulations of Gcn5 in PSC also have variable effects across hemocyte subpopulations which is not explored in the manuscript. Interestingly, also the domain deletion constructs show a differential effect on blood cell differentiation when altered solely in the prohemocytes which is not explained. While Gcn5 can be seen in all sections of the lymph gland in the first figure, under the HHLT-Gal4 and Hml-Gal4, Gcn5 looks cytoplasmic and almost completely excluded from the nucleus strikingly unlike Gcn5 expression under the Collier-Gal4 and Dome-Gal4. The rest of the experiments in the manuscript are done with multiple promoters, with autophagy flux measured by modulating Gcn5 with a pan hemocyte promoter, but the mTORC1-Gcn5 axis is explored using chemical modulators which affect the whole of the lymph gland (Fig7) or using two pro-hemocyte promoters (Fig8).

2. The knockdown of Gcn5 seems to affect the gland size (A compared to B and C). Since mTORC1 is a central regulator of cell size, it is possible that some of the effects seen in these knockdowns are potentially through mTORC1 affecting size suggesting that the signalling axis between mTORC1 and Gcn5 might not be a one-way axis as suggested in Figure 9. Also, this would mean that in experiments where absolute cell counts of crystal cells or niche cells are used to assess blood cell differentiation, further analysis to consider total cell numbers in the lymph gland would strengthen the manuscript.

3. A genetic manipulation of mTORC1 specifically in the pro hemocytes would strengthen the role of mTORC1 in the pathway rather than the chemical modulation which affects the whole of the lymph gland.

Reviewer #2 (Public Review):

Based on the corresponding author's response, the questions I raised were not addressed for various reasons. This is not necessarily a negative. The authors indicated that most of the points raised will be addressed in a separate manuscript. Specifically, the Cdu1 targeting of IkBa. They mentioned intriguing findings regarding IkBa in cells infected with a cdu1-null strain C. trachomatis in their response to reviewers. Similar to this, there appears to be a planned manuscript that will address the question of the timing of CTL0480's function in inclusion extrusion.

The lack of more direct infection-related evidence of Cdu1 interaction with various type III effectors was raised; and the authors attributed this to technical difficulties and low abundance of starting materials. It was not clear if they tried other approaches to demonstrate interaction.

Another suggestion was the quantitation of the three target effectors of Cdu1 in wild type and cdu1-null background. The authors provided western blot data and immunofluorescence images that revealed potential differences in stability/turnover kinetics. The authors might want to discuss the implications of the different kinetics of stability/turnover. For example, if all three proteins are necessary for optimal extrusion of inclusions, and concertedly act to mediate this process, all three would need to be present at the required levels. Could this be a temporal regulation strategy? Does acetylation also regulate function, interactions, etc.?

In short, the response to some of the questions is forthcoming in the form of follow-up manuscripts. New observations on the different stability profiles could be elaborated in the Discussion section, with a brief discussion on functional and/or regulatory implications.

Reviewer #2 (Public Review):

The authors succeed in generalizing the pre-alignment procedure for their cell identification method to allow it to work effectively on data with only small subsets of cells labeled. They convincingly show that their extension accurately identifies head angle, based on finding auto florescent tissue and looking for a symmetric l/r axis. They demonstrate method works to allow the identification of a particular subset of neurons. Their approach should be a useful one for researchers wishing to identify subsets of head neurons in C. elegans, and the ideas might be useful elsewhere.

The authors also assess the relative usefulness of several atlases for making identity predictions. They attempt to give some additional general insights on what makes a good atlas, but here insights seem less clear as available data does not allow for experiments that cleanly decouple: 1. the number of examples in the atlas 2. the completeness of the atlas. and 3. the match in strain and imaging modality discussed. In the presented experiments the custom atlas, besides the strain and imaging modality mismatches discussed is also the only complete atlas with more than one example. The neuroPAL atlas, is an imperfect stand in, since a significant fraction of cells could not be identified in these data sets, making it a 60/40 mix of Openworm and a hypothetical perfect neuroPAL comparison. This waters down general insights since it is unclear if the performance is driven by strain/imaging modality or these difficulties creating a complete neuroPal atlas. The experiments do usefully explore the volume of data needed. Though generalization remains to be shown the insight is useful for future atlas building that for the specific (small) set of cells labeled in the experiments 5-10 examples is sufficient to build a accurate atlas.

Reviewer #2 (Public Review):

The authors aimed to examine the role of a group of neurons expressing Foxb1 in behaviors through projections to the dlPAG. Standard chemogenetic activation or inhibition and optogentic terminal activation or inhibition at local PAG were used and results suggested that, while activation led to reduced locomotion and breathing, inhibition led to a small degree of increased locomotion.

The observed effects on breathing are evident and dramatic. However, due to the circumstance that does not permit to perform additional experiments, the conclusion is not as strong as it could be.

Reviewer #2 (Public Review):

Summary<br /> The study used eye tracking with a focus on pupillometry to examine how infants can learn to distinguish between informative and uninformative visual cues. Infants (n = 30, mean age = 8.2-months-old) viewed displays consisting of a sequence of stimuli: a fixation point, a central cue that predicted a subsequent informative or uninformative signal, the signal itself, and the target event (a cartoon animal, referred to as the reward). The key results are that: (1) pupil size differs depending on whether the infants anticipated an informative or uninformative signal, (2) this difference develops across trials, consistent with a slow learning process, and (3) there is rapid generalization when new shapes were introduced that shared features with the informative vs uninformative cues. The study complements a rich literature, including from this same group, showing that children are sensitive to information gains, and is interesting and important in revealing that pupil size is a physiological marker of information anticipation. We have several comments and concerns and believe that addressing them would substantially strengthen the manuscript.

Major points are related to interpretation, statistical robustness, and clarity

1. There is a tendency to overinterpret the findings.<br /> a. Throughout, the authors interpret the findings as meaning that pupil size tracks the "value" of information; however, the results do not demonstrate conclusively whether, or what kind of value information has in this task. A natural hypothesis is that infants are intrinsically motivated to predict - i.e., value the ability to predict the target event as early as possible. In a supplementary figure, the authors present evidence that infants indeed fixate on the target event sooner after seeing informative vs uninformative cues, consistent with the idea that they use the information for improving predictions. However, those results are not fully convincing, as we detail in point 2. Most importantly, the analysis is not integrated or even mentioned in the main analyses analysis. Making the link between the pupil reaction and the use of the information would greatly strengthen the paper (whether or not the supplementary findings hold up to more thorough scrutiny). Either this link should be made and discussed, or the authors should soften their conclusions about the utility of the informative cues.

b. On line 236, the text states that the evidence "...supports the growing body of evidence indicating that infants are proactive in shaping their learning environment by searching for and focusing on information-rich stimuli". The results do not show that the infants search for information, only that they have a pupil reaction that differentiates between informative and uninformative stimuli.

c. On lines 248-249, it seems a stretch to relate the changes in pupil dilation to a shift in information value onto the cue. Without some other measure (e.g., EEG), this remains speculative. While I believe the suggestion is plausible, the language should be softened to highlight this as a follow-up research question that the present research cannot directly speak to.

2. Several findings are statistically weak and several analyses are insufficiently controlled.

a. The analysis in Supplementary Figure 2, which shows that the latencies of target fixations are shorter after informative vs uninformative cues, raises several questions.<br /> i. We were unable to fully test these analyses as the OSF project seems to only contain latency data for 33 participants (including 22 of the 30 that remain in the final sample).<br /> ii. The results are described as revealing a significant difference, but the 89% confidence interval of the difference contains 0. How did the authors establish significance here?<br /> iii. How do the authors distinguish incidental fixations (which just happened to land near the target) from true predictive gaze shifts? Fixations were pooled if they occurred from 1.25 seconds before to 1 second after target onset. This is sufficient time for the eye to move in and out of the window several times. The authors should analyse the distributions of fixation durations to rule out various artifacts unrelated to target prediction.<br /> iv. Latencies to fixation were standardized, bringing the mean across each participant to 0, and yet the statistical model includes a random intercept; is there a justification for this?<br /> v. Standardizing removes information about whether fixations were proactive or reactive. It would be very interesting to see if/how information affects these two differently.<br /> vi. Since informativeness was learned across trials, it seems desirable that the model should include as random effects a trial number and an interaction between trial number and informativeness. This would allow a comparison between learning to predict and the pupil reaction. Are infants who have a stronger (or earlier) pupil reaction also more likely to show stronger learning to anticipate?

b. The main finding that pupil size differs between informative and uninformative cues is based on a 3-second analysis window. This long window most likely spans many saccades, which can affect pupil size on its own or by bringing the eye on or off visual stimuli. There is no analysis to show that the statistics of saccades or fixation locations are equivalent between the two trial types - but this is necessary to convincingly rule out a spurious artifact.

c. The second main finding that the effect of informativeness grows across trials seems statistically weak. The text (line 138) states that the interaction had a beta of 0.002, which was equal to the lower border of the 89%HDI ([0.002, 0.003]). For the second claim that pupil size decreased across informative trials, the beta is -0.002, and 89% HID is non-existent - i.e., [-0.002, -0.002]. (In general, the authors should check their numbers more carefully and make sure they are presented with a degree of precision that allows the reader to interpret them meaningfully.

d. The analyses do not indicate how well the TD model fits; we are shown only that it fits better than a linear model. On line 177 a correlation analysis is mentioned between the data and model, but the statistic cited for this test on line 179 is a mean beta coefficient, so it is impossible to know what this means. An analysis of goodness of fit or, at the very least, a figure superimposing the model and data, would be much more convincing.

3. The descriptions are very unclear in some key parts of the paper

a. The description of the TD model applied to pupil learning (starting on line 391) is very unclear. The model has to include some measure of informativeness - i.e., the match between the cued and true target location - but it is unclear how this was formalized. It is also very unclear how time within the trial is incorporated (the meaning of the TDE equation).

b. The description of the generalization analysis (Fig. 5) is also very unclear. Every single sentence in it evoked some confusion, so I will go through them one by one. "A Bayesian additive model showed that infants' pupil dilation was reduced for novel cues." Reduced relative to what? "This was specific to those novel cues that shared the features of the familiar informative cues (estimated mean difference = -0.05, 89%HDI = [-0.062, -0.038])." All the novel cues shared features with the informative cues; do the authors mean the novel cues that had the critical feature indicative of the informative cue? "The size of this effect approximated the difference between conditions that were observed for familiar stimuli (estimated mean difference = -0.067, 89% HDI = [-201 0.077, -0.057])." What is "this effect"? "Crucially, this difference was not observable at the start of the task, when the familiar stimuli were first introduced (estimated mean difference = -0.007, 89%HDI = [-0.015, 0.001])." At the start of the task, the stimuli were novel, and not familiar.

Reviewer #2 (Public Review):

Summary:<br /> The authors investigate sub-skin surface deformations to a number of different, relevant tactile stimuli, including pressure and moving stimuli. The results demonstrate and quantify the tension and compression applied from these types of touch to fingerprint ridges, where pressure flattens the ridges. Their study further revealed that on lateral movement, prominent vertical shearing occurred in ridge deformation, with somewhat inconsistent horizontal shear. This also shows how much the deeper skin layers are deformed in touch, meaning the activation of all cutaneous mechanoreceptors, as well as the possibility of other deeper non-cutaneous mechanoreceptors.

Strengths:<br /> The paper has many strengths. As well as being impactful scientifically, the methods are sound and innovative, producing interesting and detailed results. The results reveal the intricate workings of the skin layers to pressure touch, as well as sliding touch over different conditions. This makes it applicable to many touch situations and provides insights into the differential movements of the skin, and thus the encoding of touch in regards to the function of fingerprints. The work is very clearly written and presented, including how their work relates to the literature and previous hypotheses about the function of fingerprint ridges. The figures are very well-presented and show individual and group data well. The additional supplementary information is informative and the video of the skin tracking demonstrates the experiments well.

Weaknesses:<br /> There are very few weaknesses in the work, rather the authors detail well the limitations in the discussion. Therefore, this opens up lots of possibilities for future work.

Impact/significance:<br /> Overall, the work will likely have a large impact on our understanding of the mechanics of the skin. The detail shown in the study goes beyond current understanding, to add profound insights into how the skin actually deforms and moves on contact and sliding over a surface, respectively. The method could be potentially applied in many other different settings (e.g. to investigate more complex textures, and how skin deformation changes with factors like dryness and aging). This fundamental piece of work could therefore be applied to understand skin changes and how these impact touch perception. It can further be applied to understand skin mechanoreceptor function better and model these. Finally, the importance of fingertip ridges is well-detailed, demonstrating how these play a role in directly shaping our touch perception and how they can shape the interactions we have with surfaces.

Reviewer #3 (Public Review):

The authors utilize chimpanzee-human hybrid cell lines to assess cis-regulatory evolution. These hybrid cell lines offer a well-controlled environment, enabling clear differentiation between cis-regulatory effects and environmental or other trans effects.<br /> In their research, Wang et al. expand the range of chimpanzee-human hybrid cell lines to encompass six new developmental cell types derived from all three germ layers. This expansion allows them to discern cell type-specific cis-regulatory changes between species from more pleiotropic ones. Although the study investigates only two iPSC clones, the RNA- and ATAC-seq data produced for this paper is a valuable resource.

The authors begin their analysis by examining the relationship between allele-specific expression (ASE) as a measure of species divergence and cell type specificity. They find that cell-type-specific genes exhibit more divergent expression. By integrating this data with measures of constraint within human populations, the authors conclude that the increased divergence of tissue-specific genes is, at least in part, attributable to positive selection. A similar pattern emerges when assessing allele-specific chromatin accessibility (ASCA) as a measure of divergence of cis-regulatory elements (CREs) in the same cell lines.

By correlating these two measures, the authors identify 95 CRE-gene pairs where tissue-specific ASE aligns with tissue-specific ASCA. Among these pairs, the authors select two genes of interest for further investigation. Notably, the authors employ an intriguing machine learning approach in which they compare the inferred chromatin state of the human sequence with that of the chimpanzee sequence to pinpoint putatively causal variants.

Overall, this study delves into the examination of gene expression and chromatin accessibility within hybrid cell lines, showcasing how this data can be leveraged to identify potential causal sequence differences underlying between-species expression changes.

All in all most conclusions appear solid, with the exception of the interpretation of a cell type/state identification machine learning model to pinpoint putatively causal variants. The described variants lack any functional validation and there is no data that measure the certainty of the results.

Reviewer #2 (Public Review):

In this work, the authors elaborate on an analytically tractable, continuous-attractor model to study an idealized neural network with realistic spiking phase precession/procession. The key ingredient of this analysis is the inclusion of a mechanism for slow firing-rate adaptation in addition to the otherwise fast continuous-attractor dynamics. The latter continuous-attractor dynamics classically arises from a combination of translation invariance and nonlinear rate normalization.

For strong adaptation/weak external input, the network naturally exhibits an internally generated, travelling-wave dynamics along the attractor with some characteristic speed. For small adaptation/strong external stimulus, the network recovers the classical externally driven continuous-attractor dynamics. Crucially, when both adaptation and external input are moderate, there is a competition with the internally generated and externally generated mechanisms leading to an oscillatory tracking regime. In this tracking regime, the population firing profile oscillates around the neural field tracking the position of the stimulus. The authors demonstrate by a combination of analytical and computational arguments that oscillatory tracking corresponds to realistic phase precession/procession. In particular the authors can account for the emergence of unimodal and bimodal cells, as well as some other experimental observations with respect the dependence of phase precession/procession on the animal's locomotion.

The strengths of this work are at least three-fold: 1) Given its simplicity, the proposed model has a surprisingly large explanatory power of the various experimental observations. 2) The mechanism responsible for the emergence of precession/procession can be understood as a simple yet rather illuminating competition between internally driven and externally driven dynamical trends. 3) Amazingly, and under some adequate simplifying assumptions, a great deal of analysis can be treated exactly, which allows for a detailed understanding of all parametric dependencies. This exact treatment culminates with a full characterization of the phase space of the network dynamics, as well as the computation of various quantities of interest, including characteristic speeds and oscillating frequencies.

As mentioned by the authors themselves, the main limitation of this work is that it deals with a very idealized model and it remains to see how the proposed dynamical behaviors would persists in more realistic models. For example, the model is based on a continuous attractor model that assumes perfect translation-invariance of the network connectivity pattern. Would the oscillating tracking behavior persist in the presence of connection heterogeneities? Another limitation is that the system needs to be tuned to exhibit oscillation within the theta range and that this tuning involves a priori variable parameters such as the external input strength. Is the oscillating-tracking behavior overtly sensitive to input strength variations? The author mentioned that an external pacemaker can serve to drive oscillation within the desired theta band but there is no evidence presented supporting this. A final and perhaps secondary limitation has to do with the choice of parameter, namely the time constant of neural firing which is chosen around 3ms. This seems rather short given that the fast time scale of rate models (excluding synaptic processes) is usually given by the membrane time constant, which is typically about 15ms. I suspect this latter point can easily be addressed.

Reviewer #2 (Public Review):

Hong and collaborators investigated variations in the amount of synaptic proteins in plasma extracellular vesicles (EV) in Parkinson's Disease (PD) patients on one-year follow-up. Their findings suggest that plasma EV synaptic proteins may be used as clinical biomarkers of PD progression.

It is a preliminary study using semi-quantitative analysis of synaptic proteins.

The authors have a cohort of PD patients with clinical examination and a know-how on EV purification. Regarding this latter part, they may improve their description of EV purification. EV may be broken into smaller size EV after freezing. Does it explain the relatively small size in their EV preparation? Do the authors refer to the MISEV guidelines for EV purity? Regarding synaptic protein quantification, the choice of western blotting may not be the best one. ELISA and other multiplex arrays are available. How the authors do justify their choice? Do the authors try to sort plasma EV by membrane-associated neuronal EV markers using either vesicle sorting or immunoprecipitation?

Many technical aspects may be improved. Such technical questions weakened the authors' conclusions.

Reviewer #2 (Public Review):

The authors made an atlas of single-cell transcriptome of on a pure population of leukocytes isolated from the brain of adult Tg(cd45:DsRed) transgenic animals by flow cytometry. Seven major leukocyte populations were identified, comprising microglia, macrophages, dendritic-like cells, T cells, natural killer cells, innate lymphoid-like cells, and neutrophils. Each cluster was analyzed to characterize subclusters. Among lymphocytes, in addition to 2 subclusters expressing typical T cell markers, a group of il4+ il13+ gata3+ cells was identified as possible ILC2. This hypothesis is supported by the presence of this population in rag2KO fish, in which the frequency of lck and zap70+ cells is strongly reduced. The use of KO lines for such validations is a strength of this work (and the zebrafish model).

The subcluster analysis of mpeg1.1 + myeloid cells identified 4 groups of microglial cells, one novel group of macrophage-like cells (expressing s100a10b, sftpbb, icn, fthl27, anxa5b, f13a1b and spi1b), and several groups of DC like cells expressing the markers siglec15l, ccl19a.1, ccr7, id2a, xcr1a.1, batf3, flt3, chl1a and hepacam2. Combining these new markers and transgenic reporter fish lines, the authors then clarified the location of leukocyte subsets within the brain, showing for example that DC-like cells stand as a parenchymal population along with microglia. Reporter lines were also used to perform a detailed analysis of cell subsets, and cross with a batf3 mutant demonstrated that DC-like cells are batf3 dependent, which was similar to mouse and human cDC1. Finally, analysis of classical mononuclear phagocyte deficient zebrafish lines showed they have reduced numbers of microglia but exhibit distinct DC-like cell phenotypes. A weakness of this study is that it is mainly based on FACS sorting, which might modify the proportion of different subtypes.

This atlas of zebrafish brain leukocytes is an important new resource for scientists using the zebrafish models for neurology, immunology, and infectiology, and for those interested in the evolution of the brain and immune system.

Reviewer #2 (Public Review):

Summary<br /> This study deepens the former authors' investigations of the mechanisms involved in gating the long-term consolidation of an associative memory (LTM) in Drosophila melanogaster. After having previously found that LTM consolidation 1. costs energy (Plaçais and Préat, Science 2013) provided through pyruvate metabolism (Plaçais et al., Nature Comm 2017) and 2. is gated by the increased tonic activity in a type of dopaminergic neurons ('MP1 neurons') following only training protocol relevant for LTM, i.e. interspaced in time (Plaçais et al., Nature Neuro 2012), they here dig into the intra-cell signalling triggered by dopamine input and eventually responsible for the increased mitochondria activity in Kenyon Cells. They identify a particular PKC, PKCδ, as a major molecular interface in this process and describe its translocation to mitochondria to promote pyruvate metabolism, specifically after spaced training.

Methodological approach<br /> To that end, they use RNA interference against the isozyme PKCδ, in a time-controlled way and in the whole Kenyon cell populations or in the subpopulation forming the α/β lobe. This knock-down decreased the total PKCδ mRNA level in the brain by ca. 30%, and is enough to observe decreased in flies performances for LTM consolidation. Using Pyronic, a sensor for pyruvate for in vivo imaging, and pharmacological disruption of mitochondrial function, the authors then show that PKCδ knock-down prevents a high level of pyruvate from accumulating in the Kenyon cells at the time of LTM consolidation, pointing towards a role of PKCδ in promoting pyruvate metabolism. They further identify the PDH kinase PDK as a likely target for PKCδ since knocking down both PKCδ and PDK led to normal LTM performances, likely counterbalancing PKCδ knock-down alone.

To understand the timeline of PKCδ activation and to visualise its mitochondrial translocation in a subpart of Mushroom body lobes they imported in fruitfly the genetically-encoded FRET reporters of PKCδ, δCKAR, and mitochondria-δCKAR (Kajimoto et al 2010). They show that PKCδ is activated to the sensor's saturation only after spaced training, and not other types of training that are 'irrelevant' for LTM. Further, adding thermogenetic activation of dopaminergic neurons and RNA interference against Gq-coupled dopamine receptor to FRET imaging, they identify that a dopamine-triggered cascade is sufficient for the elevated PKCδ-activation.

Strengths and weaknesses<br /> The authors use a combination of new fluorescent sensors and behavioral, imaging, and pharmacological protocols they already established to successfully identify the molecular players that bridge the requirement for spaced training/dopaminergic neurons MP1 oscillatory activity and the increased metabolic activity observed during long-term memory consolidation.

The study is dense in new exciting findings and each methodological step is carefully designed. Almost all possible experiments one could think of to make this link have been done in this study, with a few exceptions that do not prevent the essential conclusions from being drawn.

The discussion is well conducted, with interesting parallels with mammals, where the possibility that this process takes place as well is yet unknown.

Impact<br /> Their findings should interest a large audience:<br /> They discover and investigate a new function for PKCδ in regulating memory processes in neurons in conjunction with other physiological functions, making this molecule a potentially valid target for neuropathological conditions. They also provide new tools in drosophila to measure PKCδ activation in cells. They identify the major players for lifting the energetic limitations preventing the formation of a long-term memory.

Reviewer #2 (Public Review):

Summary:<br /> In this manuscript, Shore et al. investigate the consequent changes in excitability and synaptic efficacy of diverse neuronal populations in an animal model of juvenile epilepsy. Using electrophysiological patch-clamp recordings from dissociated neuronal cultures, the authors find diverging changes in two major populations of inhibitory cell types, namely somatostatin (SST)- and parvalbumin (PV)-positive interneurons, in mice expressing a variant of the KCNT1 potassium channel. They further suggest that the differential effects are due to a compensatory increase in the persistent sodium current in PV interneurons in pharmacological and in silico experiments.

Strengths:<br /> 1) Heterozygous KCNT1 gain of function variant was used which more accurately models the human disorder.<br /> 2) The manuscript is clearly written, and the flow is easy to follow. The authors explicitly state the similarities and differences between the current findings and the previously published results in the homozygous KCNT1 gain of function variant.<br /> 3) This study uses a variety of approaches including patch clamp recording, in silico modeling, and pharmacology that together make the claims stronger.<br /> 4) Pharmacological experiments are fraught with off-target effects and thus it bolsters the authors' claims when multiple channel blockers (TTX and VU170) are used to reconstruct the sodium-activated potassium current. Having said that, it would be helpful to see the two drug manipulations be used in the same experiment. Notably, does the more selective blocker VU170 mimic the results of TTX for NFS GABAergic cells in Figure 2? And does it unmask a genotype difference for FS GABAergic cells like the one seen in PV interneurons in Figure 5C3.

Weaknesses:<br /> 1) This study relies on recordings in dissociated cortical neurons. Although specific WT interneurons showed intrinsic membrane properties like those reported for acute brain slices, it is unclear whether the same will be true for those cells expressing KCNT1 variants. This reviewer highly recommends confirming some of the key findings using an ex vivo slice preparation. This is especially important given the discrepant result of reduced excitability of PV cells reported by Gertler et al., 2022 (cited here in the manuscript but not discussed in this context) in acute hippocampal slices for a different KCTN1 gain of function variant.<br /> 2) It is unclear how different pieces of results fit together to form a story about the disease pathophysiology. For example, hyperexcitability of PV cells would suggest more inhibition which would counter seizure propensity. However, spontaneous inhibitory postsynaptic currents show no change in pyramidal neurons. Moreover, how do the authors reconcile that the reductions in synaptic inputs onto interneurons in Figure 3B with the increases in Figure 8? This should be discussed.<br /> 3) Similarly, the results in this work are not entirely internally consistent. For example, given the good correspondence between FS and NFS GABAergic cells with PV and SST expression, why are FS GABAergic cells hyperexcitable in Figure 1? If anything, there is a tendency to show reduced excitability like the NFS GABAergic cells. Also, why do the WT I-V curves look so different between Figures 2 and 5? This reviewer suggests at least a brief explanation in the discussion.<br /> 4) Given the authors' claim that the KCNT1 activation curve is a major contributor to the observed excitability differences in specific GABA cell subtypes, it would be helpful to directly measure the activation curve in the variants experimentally as was done for WT KCNT1 in Figure 6A and use the derived kinetics in the compartmental model.

Reviewer #2 (Public Review):

Summary:

In this study, Zhao, Nern, et al. investigate a population of neurons in the optic lobe of Drosophila melanogaster that process optic flow, relative motion between the insect's eyes and its environment that is generated during flight and provides useful information to the fly about its own self-motion. Although a sample of these Lobula Plate Tangential (LPT) neurons has been studied across Diptera in prior work, the full population has not been exhaustively and thoroughly cataloged in a single species, limiting our understanding of how LPT tuning properties across the population convey features of optic flow fields relevant to downstream motor regions.

Through extensive manual reconstructions in a fly electron microscopy volume, the authors of this study identify 58 LPT neurons in the fruit fly encompassing previously studied Horizontal and Vertical cells and novel cells that have not been previously characterized. Using the detailed anatomy of each cell and knowledge of upstream T4/T5 selectivity, the authors derive the predicted motion pattern map (PMPM) of each neuron. To understand how optic flow field tunings of individual LPTs align with global optic flow patterns flies are expected to encounter during flight such as translation and rotation, the authors compute the average angular difference between each PMPM and idealized rotation and translation optic flow fields. The authors also map individual LPTs to their counterparts in a second fly brain to explore LPT-LPT connectivity and downstream connectivity to central brain neuropils. They find that distinct groupings of LPTs have diverse downstream connectivity patterns and that downstream neurons align more closely to global optic flow fields that are expected during flight. This study is a valuable resource to researchers studying motion vision in the insect brain and is of interest to researchers studying sensorimotor processing by providing hypotheses for how optic flow information is integrated downstream to guide fly behavior.

Strengths:

A key strength of this study is the thoroughness with which the authors comprehensively identify individual LPT neurons in the FAFB volume. They not only conduct an impressive number of careful manual reconstructions to recover individual LPTs, but they also attempt, and often succeed, to map each individual neuron to its counterpart in light microscopy, studies across Diptera, and available auto-segmented connectome datasets such as FlyWire, FAFB-FFN1, and Hemibrain. The authors are similarly thorough when surveying individual LPT properties such as neurotransmitter identities, in some cases using multiple datasets to reconcile ambiguous neurotransmitter predictions. The care with which the complete LPT population has been identified establishes this study as a useful resource for future studies of insect motion.

In addition to providing a comprehensive catalog of individual LPTs, the authors also contextualize their findings within broader sensorimotor circuitry by considering connectivity between LPTs and from LPTs to downstream regions. Exploration of structure in downstream connectivity suggests that optic flow information is directed to various central brain neuropils through specific groups of LPTs. With some additional analyses, these results broaden the scope of this study by providing useful hypotheses for sensorimotor circuit organization.

Weaknesses:

A novel method introduced in this study is the derivation of individual LPT-predicted motion pattern maps (PMPMs) using T4 preferred directions and LPT morphology. Although this method underlies core findings in this study, such as alignment to global optic flow fields and properties of downstream integration, aspects of the methods used to derive PMPMs are not explained sufficiently well, particularly in the main text. For example, in the Methods, the authors briefly describe the process of computing a weighted sum of T4 preferred directions to obtain the PMPM for each LPT, but a detailed understanding of these preferred directions combined is missing in Figure 2 or the associated descriptions in the main text. It is also not clear how PMPMs are derived in cases where LOP layer coverages are overlapping (for example VS 13-1 in Figure 3) to yield smooth PMPMs. In addition, it is not clear how the PMPMs of bilateral LPTs such as LPT-45 and LPT-50 in Figure 4 were integrated to compute downstream target composite PMPMs. Finally, all the PMPMs were derived from the T4 preferred direction that relies on the ommatidial viewing directions ("Eyemap") introduced in Zhao et al. 2022. It is also important for the current study to give an indication of how sensitive their results are to possible inaccuracies in this map and derived T4/T5 direction selectivities.

Although the authors explore some features of connectivity from LPT to downstream partners (Figure 6), there is a lack of reconciliation of these findings with individual LPT properties explored earlier in the study, such as those presented in Figures 2-4. In that sense, there is a disconnect between the two parts of the manuscript (and a missed opportunity). For instance, an important follow-up analysis would be to use knowledge about LPT-LPT connectivity to better predict effective PMPMs of LPTs taking into account network effects. This extension would lead to a better understanding of how LPT-LPT interactions shape optic flow responses in the LOP. In addition, in Figure 6 Supplement 2 (which I recommend to move to the main figures), the authors show that LPTs can be grouped together based on similarity of output connectivity (Panel B-D) and that this structure corresponds to output synapses located in different groups of central brain neuropils. However, they do not attempt to explicitly link these groupings with individual LPT PMPMs, alignment to global optic flow patterns, LPT layer enervation, cell morphologies, and input connectivity patterns. Such an analysis would be an important step to bring the manuscript together and to get a better understanding of the organization of the whole system.

Reviewer #2 (Public Review):

Summary:

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

Strengths:

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

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

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

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

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

Weaknesses:

Whether micro-waves are associated with the age of mice was not quantified. This would be good to know and the authors do have this data.

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

The CaMKII-Cre approach for flexed-syn-GCaMP expression shows no micro-waves and is convincing, but it is only from 2 animals, even though both had no micro-waves.

The authors state in their Discussion that even without observable microwaves, a syn-Ca2+-indicator transduction strategy could still be problematic. This may be true, but they do not check this in their analysis, so it remains unknown.

Reviewer #2 (Public Review):

Summary:<br /> Non-canonical Wnt signaling plays an important role in morphogenesis, but how different components of the pathway are required to regulate different developmental events remains an open question. This paper focuses on elucidating the overlapping and distinct functions of dact1 and dact2, two Dishevelled-binding scaffold proteins, during zebrafish axis elongation and craniofacial development. By combining genetic studies, detailed phenotypic analysis, lineage tracing, and single-cell RNA-sequencing, the authors aimed to understand (1) the relative function of dact1/2 in promoting axis elongation, (2) their ability to modulate phenotypes caused by mutations in other non-canonical wnt components, and (3) pathways downstream of dact1/2.

Strong qualitative evidence was provided to support dact1/2's role in genetically modulating non-canonical wnt signaling to regulate body axis elongation and the morphology of the anterior neurocranium (ANC). However, there is currently insufficient evidence supporting the author's claim that suppression of calpain 8 by dact1/2 is important for craniofacial development and that "embryonic fields determined during gastrulation affect the CNCC ability to contribute to the craniofacial skeleton".

Strengths:<br /> (1) The generation of dact1/2 germline mutants and the use of genetic approaches to dissect their genetic interactions with wnt11f2 and gpc4 provide unambiguous and consistent results that inform the relative functions of dact1 and dact2, as well as their combined effects.

(2) Because the anterior neurocranium exhibits a spectrum of phenotypes in different genetic mutants, it is a useful system for studying how tissue morphology can be modulated by different components of the same pathway, as demonstrated in this study.

(3) The authors leveraged lineage tracing by photoconversion to dissect how dact1/2 differentially impacts the ability of different cranial neural crest populations to contribute to the anterior neurocranium. This revealed that distinct mechanisms can lead to similar phenotypes in different mutants.

Weaknesses:<br /> (1) While the qualitative data show altered morphologies in each mutant, quantifications of these phenotypes are lacking in several instances, making it difficult to gauge reproducibility and penetrance, as well as to assess the novel ANC forms described in certain mutants.

(2) Germline mutations limit the authors' ability to study a gene's spatiotemporal functional requirement. They therefore cannot concretely attribute nor separate early-stage phenotypes (during gastrulation) to/from late-stage phenotypes (ANC morphological changes).

(3) Given that dact1/2 can regulate both canonical and non-canonical wnt signaling, this study did not specifically test which of these pathways is altered in the dact1/2 mutants, and it is currently unclear whether disrupted canonical wnt signaling contributes to the craniofacial phenotypes, even though these phenotypes are typical non-canonical wnt phenotypes.

(4) The use of single-cell RNA sequencing unveiled genes and processes that are uniquely altered in the dact1/2 mutants, but not in the gpc4 mutants during gastrulation. However, how these changes lead to the manifested ANC phenotype later during craniofacial development remains unclear. The authors showed that calpain 8 is significantly upregulated in the mutant, but the fact that only 1 out of 142 calpain-overexpressing animals phenocopied dact1/2 mutants indicates the complexity of the system.

(5) Craniofacial phenotypes observed in this study are attributed to convergent extension defects but convergent extension cell movement itself was not directly examined, leaving open if changes in other cellular processes, such as cell differentiation, proliferation, or oriented division, could cause distinct phenotypes between different mutants.

Reviewer #2 (Public Review):

George and colleagues present a novel open-source toolbox to model rodent locomotor patterns and electrophysiological responses of spatially modulated neurons, such as hippocampal "place cells". The present manuscript describes a comprehensive Python package ("RatInABox") with powerful capabilities to simulate a variety of environments, exploratory behaviors and concurrent responses of a variety of cell types. In addition, they provide the tools to expand these basics functions and potentially multiple different model designs, new cell types or more complex neural network architectures. The manuscript also illustrated several simple application cases. The authors have also created a comprehensive GitHub repository with more detailed explanations, tutorials and example scripts. Overall, I found both the manuscript and associated repository very clear, well written and easy the scrips easy to follow and implement, to a superior level of many commercial software packages. RatInABox fills several existing gaps in the literature and features important improvements over previous approaches; for example, the implementation of continuous 2D environments instead of tabularized state spaces. I believe this toolbox will be of great interest for many researchers in the field of spatial navigation and beyond and provide them with a remarkably powerful and flexible tool. I don't have any major issues with the manuscript. However, the manuscript can be further improved by clarifying some aspects of the toolbox, discussing its limitations and biological plausibility.

Reviewer #2 (Public Review):

In this study, Lopes and colleagues provide convincing evidence to support the potential for gene therapy to restore expression of heparanase-2 (Hpse2) in mice mutant for this gene, as occurs in urofacial syndrome. Beyond symptomatic relief for the consequences of outlet obstruction that results from Hpse2 mutation, no treatments exist. Building on prior studies describing the nature of urinary tract dysfunction in Hpse2 mutant mice, the authors applied a gene therapy approach to determine whether gene replacement could be achieved, and if so, whether restoration of HPSE2 expression could mitigate the urinary tract dysfunction and present a potential cure. Using an AAV9 viral vector encoding HPSE2, the authors performed gene replacement in neonatal wild-type or Hpse2 mutant mice and determined gene and protein expression as well as the impact on bladder outflow tract and bladder body physiology in juvenile mice. In addition to dose-dependent transduction of liver and pelvic ganglia (that innervate the bladder) with HPSE2, and demonstration of increased HPSE2 protein in Hpse2 mutant mice, the authors showed restoration of nerve-evoked outflow tract relaxation and bladder body contraction, both of which were deficient in mutant mice. They also showed that the viral vector-based approach was not deleterious to weight gain or to liver morphology. Based on these findings the authors concluded that AAV9-based HPSE2 replacement is feasible and safe, mitigates the physiological deficits in outflow tract and bladder tissue from Hpse2 mutant mice, and provides a foundation for gene replacement approaches for other genes implicated in lower urinary tract disorders.

Strengths include a rigorous experimental design, solid data in support of the conclusions, and a discussion of the limitations of the approach.

Weaknesses include a lack of discussion of the basis for differences in carbachol sensitivity in Hpse2 mutant mice, limited discussion of bladder tissue morphology in Hpse2 mutant mice, some questions over the variability of the functional data, and a need for clarification on the presentation of statistical significance of functional data

Reviewer #3 (Public Review):

Summary: The authors used cross-linking of a known P-gp substrate in combination with single particle cyro-EM to investigate the translocation pathway of this important ABC transporter. Based on the results of this study, a new translocation mechanism is proposed that is supported by the data. While only one substrate was used, the data obtained are convincing. In addition, the proposed model will stimulate new experiments from other laboratories to proof or disproof the model.

Strengths: the combination of cross-linking and structural biology allowed novel insights in the translocation pathway of ABCB1

Weaknesses: While only one substrate was used, the data obtained are convincing. In addition, the proposed model will stimulate new experiments from other laboratories to proof or disproof the model.

Reviewer #2 (Public Review):

Summary:<br /> This important study by LaBerge and co-authors seeks to understand the causal drivers of faculty gender demographics by quantifying the relative importance of faculty hiring and attrition across fields. They leverage historical data to describe past trends and develop models that project future scenarios that test the efficacy of targeted interventions. Overall, I found this study to be a compelling and important analysis of gendered hiring and attrition in US institutions, and one that has wide-reaching policy implications for the academy. The authors have also suggested a number of fruitful future avenues for research that will allow for additional clarity in understanding the gendered, racial, and socioeconomic disparities present in US hiring and attrition, and potential strategies for mitigating or eliminating these disparities.

Strengths:<br /> In this study, LaBerge et al use data from over 268,000 tenured and tenure-track faculty from over 100 fields at more than 12,000 PhD-granting institutions in the US. The period they examine covers 2011-2020. Their analysis provides a large-scale overview of demographics across fields, a unique strength that allows the authors to find statistically significant effects for gendered attrition and hiring across broad areas (STEM, non-STEM, and topical domains).

LaBerge et al. find gendered disparities in attrition-using both empirical data and their counterfactual model-that account for the loss of 1378 women faculty across all fields between 2011 and 2020. It is true that "this number is both a small portion of academia... and a staggering number of individual careers," as ." - as this loss of women faculty is comparable to losing more than 70 entire departments. I appreciate the authors' discussion about these losses-they note that each of these is likely unnecessary, as women often report feeling that they were pushed out of academic jobs.

LaBerge et al. also find-by developing a number of model scenarios testing the impacts of hiring, attrition, or both-that hiring has a greater impact on women's representation in the majority of academic fields in spite of higher attrition rates for women faculty relative to men at every career stage. Unlike many other studies of historical trends in gender diversity, which have often been limited to institution-specific analyses, they provide an analysis that spans over 100 fields and includes nearly all US PhD-granting institutions. They are able to project the impacts of strategies focusing on hiring or retention using models that project the impact of altering attrition risk or hiring success for women. With this approach, they show that even relatively modest annual changes in hiring accumulate over time to help improve the diversity of a given field. They also demonstrate that, across the model scenarios they employ, changes to hiring drive the largest improvement in the long-term gender diversity of a field.

Future work will hopefully - as the authors point out - include intersectional analyses to determine whether a disproportionate share of lost gender diversity is due to the loss of women of color from the professoriate. I appreciate the author's discussion of the racial demographics of women in the professoriate, and their note that "the majority of women faculty in the US are white" and thus that the patterns observed in this study are predominately driven by this demographic. I also highly appreciate their final note that "equal representation is not equivalent to equal or fair treatment," and that diversifying hiring without mitigating the underlying cause of inequity will continue to contribute to higher losses of women faculty.

Weaknesses<br /> First, and perhaps most importantly, it would be beneficial to include a distinct methods section. While the authors have woven the methods into the results section, I found that I needed to dig to find the answers to my questions about methods. I would also have appreciated additional information within the main text on the source of the data, specifics about its collection, inclusion and exclusion criteria for the present study, and other information on how the final dataset was produced. This - and additional information as the authors and editor see fit - would be helpful to readers hoping to understand some of the nuance behind the collection, curation, and analysis of this important dataset.

I would also encourage the authors to include a note about binary gender classifications in the discussion section. In particular, I encourage them to include an explicit acknowledgement that the trends assessed in the present study are focused solely on two binary genders - and do not include an analysis of nonbinary, genderqueer, or other "third gender" individuals. While this is likely because of the limitations of the dataset utilized, the focus of this study on binary genders means that it does not reflect the true diversity of gender identities represented within the professoriate.

In a similar vein, additional context on how gender was assigned on the basis of names should be added to the methods section.

I do think that some care might be warranted regarding the statement that "eliminating gendered attrition leads to only modest changes in field-level diversity" (Page 6). while I do not think that this is untrue, I do think that the model scenarios where hiring is "radical" and attrition is unchanged from present (equal representation of women and men among hires (ER) + observed attrition (OA)) shows that a sole focus on hiring dampens the gains that can otherwise be addressed via even modest interventions (see, e.g., gender-neutral attrition (GNA) + increasing representation of women among hires (IR)). I am curious as to why the authors did not include an additional scenario where hiring rates are equal and attrition is equalized (i.e., GNA + ER). The importance of including this additional model is highlighted in the discussion, where, on Page 7, the authors write: "In our forecasting analysis, we find that eliminating the gendered attrition gap, in isolation, would not substantially increase representation of women faculty in academia. Rather, progress towards gender parity depends far more heavily on increasing women's representation among new faculty hires, with the greatest change occurring if hiring is close to gender parity." I believe that this statement would be greatly strengthened if the authors can also include a comparison to a scenario where both hiring and attrition are addressed with "radical" interventions.

Reviewer #2 (Public Review):

Summary:<br /> The manuscript by Herneisen et al. examines the Toxoplasma SPARK kinase orthologous to mammalian PDK1 kinase. The extracellular signals trigger cascades of the second messengers and play a central role in the apicomplexan parasites' survival. In Toxoplasma, these cascades regulate active replication of the tachyzoites, which manifests as acute toxoplasmosis, or the development into drug-resilient bradyzoites characteristic of the chronic stage of the disease. This study focuses on the poorly understood signaling mechanisms acting upstream of such second messenger kinases as PKA and PKG. The authors showed that similar to PDK1, Toxoplasma SPARK appears to regulate several AGC kinases.

Strengths:<br /> The study demonstrated a strong association of the SPARK kinase with an elongin-like SPARKEL factor and an uncharacterized AGC kinase. Using a set of standard assays, the authors determined the SPARK/SPARKEL role in parasite egress and invasion. Finally, the study presented evidence of the SPARK/SPARKEL involvement in the bradyzoite differentiation.

Weaknesses:<br /> Although the study can potentially uncover essential sensing mechanisms operating in Toxoplasma, the evidence of the SPARK/SPARKEL mechanisms is weak. Specifically, due to incomplete data analysis, the SPARK/SPARKEL-dependent phosphoregulation of AGC kinases cannot be evaluated. The manuscript requires better organization and lacks guidance on the described experiments. Although the study is built on advanced genetics, at times, it is unnecessarily complicated, raising doubts rather than benefiting the study.

Reviewer #2 (Public Review):

Lentiviral infection of primate species has been linked to the rapid mutational evolution of numerous primate genes that interact with these viruses, including genes that inhibit lentiviruses as well as genes required for viral infection. In this manuscript, Warren et al. provide further support for the diversification of CD4, the lentiviral entry receptor, to resist lentiviral infection in great ape populations. This work builds on their prior publication (Warren et al. 2019, PMCID: PMC6561292 ) and that of other groups (e.g., Russell et al. 2021, PMCID: PMC8020793; Bibollet-Ruche et al. 2019, PMCID: PMC6386711) documenting both sequence and functional diversity in CD4, specifically within (1) the CD4 domain that binds to the lentiviral envelope and (2) great ape populations with endemic lentiviruses. Thus, the paper's finding that gorilla populations exhibit diverse CD4 alleles that differ in their susceptibility to lentiviral infection is well demonstrated both here and in a prior publication.

To bolster the argument that lentiviruses are indeed the causative driver of this diversification, which seems likely from a logical perspective but is difficult to prove, Warren et al. pursue two novel lines of evidence. First, the authors reconstruct ancestral CD4 genes that predate lentiviral infection of hominid populations. They then demonstrate that resistance to lentiviral infection is a derived trait in chimpanzees and gorillas, which have been co-evolving with endemic lentiviruses, but not in humans, which only recently acquired HIV. Nevertheless, the derived resistance could be stochastic or due to drift. This argument would be strengthened by demonstrating that bonobo and orangutan CD4, which also do not have endemic lentiviruses, resemble the ancestral and human susceptibility to great-ape-infecting lentiviruses.

Second, Warren et al. provide a population genetic argument that only endemically infected primates exhibit diversifying selection, again arguing for endemic lentiviruses being the evolutionary driver. The authors compare SNP occurrence in CD4 to neighboring genes, demonstrating that non-synonymous SNP frequency is only elevated in endemically infected species. Moreover, these amino-acid-coding changes are significantly concentrated in the CD4 domain that binds the lentiviral envelope. This is a creative analysis to overcome the problem of very small sample sizes, with very few great ape individuals sequenced. The additional small number of species compared (2-3 in each group) also limits the power of the analysis; the authors could consider expanding their analysis to Old World Monkey species that do or do not have endemic lentiviruses, as well as great apes.

Overall, this manuscript lends additional support to a well-documented example of a host-virus arms race: that of lentiviruses and the viral entry receptor.

Reviewer #2 (Public Review):

In this paper, Wang and collaborators characterize the rapid evolution of the X-linked miR-506 cluster in mammals and characterize the functional reference of depleting a few or most of the miRNAs in the cluster. The authors show that the cluster originated from the MER91C DNA transposon and provide some evidence that it might have expanded through the retrotransposition of adjacent LINE1s. Although the animals depleted of most miRNAs in the cluster show normal sperm parameters, the authors observed a small but significant reduction in litter size. The authors then speculate that the depletion of most miRNAs in the cluster could impair sperm competitiveness in polyandrous mating. Using a successive mating protocol, they show that, indeed, sperm lacking most X-linked miR-506 family members is outcompeted by wild-type sperm. The authors then analyze the evolution of the miR-506 cluster and its predicted targets. They conclude that the main difference between mice and humans is the expansion of the number of target sites per transcript in humans.

The conclusions of the paper are, in most cases, supported by the data; however, a more precise and in-depth analysis would have helped build a more convincing argument in most cases.

1) In the abstracts and throughout the manuscript, the authors claim that "... these X-linked miRNA-506 family miRNA [...] have gained more targets [...] " while comparing the human miRNA-506 family to the mouse. An alternative possibility is that the mouse has lost some targets. A proper analysis would entail determining the number of targets in the mouse and human common ancestor.

2) The authors claim that the miRNA cluster expanded through L1 retrotransposition. However, the possibility of an early expansion of the cluster before the divergence of the species while the MER91C DNA transposon was active was not evaluated. Although L1 likely contributed to the diversity within mammals, the generalization may not apply to all species. For example, SINEs are closer on average than L1s to the miRNAs in the SmiR subcluster in humans and dogs, and the horse SmiR subcluster seems to have expanded by a TE-independent mechanism.

3) Some results are difficult to reconcile and would have benefited from further discussion. The miR-465 sKO has over two thousand differentially expressed transcripts and no apparent phenotype. Also, the authors show a sharp downregulation of CRISP1 at the RNA and protein level in the mouse. However, most miRNAs of the cluster increase the expression of Crisp1 on a reporter assay. The only one with a negative impact has a very mild effect. miRNAs are typically associated with target repression; however, most of the miRNAs analyzed in this study activate transcript expression.

4) More information is required to interpret the results of the differential RNA targeting by the murine and human miRNA-506 family. The materials and methods section needs to explain how the authors select their putative targets. In the text, they mention the use of four different prediction programs. Are they considering all sites predicted by any method, all sites predicted simultaneously by all methods, or something in between? Also, what are they considering as a "shared target" between mice and humans? Is it a mRNA that any miR-506 family member is targeting? Is it a mRNA targeted by the same miRNA in both species? Does the targeting need to occur in the same position determined by aligning the different 3'UTRs?

5) The authors highlight the particular evolution of the cluster derived from a transposable element. Given the tendency of transposable elements to be expressed in germ cells, the family might have originated to repress the expression of the elements while still active but then remained to control the expression of the genes where the element had been inserted. The authors did not evaluate the expression of transcripts containing the transposable element or discuss this possibility. The authors proposed an expansion of the target sites in humans. However, whether this expansion was associated with the expansion of the TE in humans was not discussed either. Clarifying whether the transposable element was still active after the divergence of the mouse and human lineages would have been informative to address this outstanding issue.

Post-transcriptional regulation is exceptionally complex in male haploid cells, and the functional relevance of many regulatory pathways remains unclear. This manuscript, together with recent findings on the role of piRNA clusters, starts to clarify the nature of the selective pressure that shapes the evolution of small RNA pathways in the male germ line.

Reviewer #2 (Public Review):

Summary:<br /> The manuscript "Cohesin still drives homologous recombination repair of DNA double-strand breaks in late mitosis" by Ayra-Plasencia et al. investigates regulations of HR repair in conditional cdc15 mutants, which arrests the cell cycle in late anaphase/telophase. Using a non-competitive MAT switching system of S. cerevisiae, they show that a DSB in telophase-arrested cells elicits a delayed DNA damage checkpoint response and resection. Using a degron allele of SMC3 they show that MATa-to-alpha switching requires cohesin in this context. The presence of a DSB in telophase-arrested cells leads to an increase in the kleisin subunit Scc1 and a partial rejoining of sister chromatids after they have separated in a subset of cells.

Strengths:<br /> The experiments presented are well-controlled. The induction systems are clean and well thought-out.

Weaknesses:<br /> The manuscript is very preliminary, and I have reservations about its physiological relevance. I also have reservations regarding the usage of MAT to make the point that inter-sister repair can occur in late mitosis.

Reviewer #2 (Public Review):

Summary:

The ability of cardiac cells to regenerate has been the object of intense (and sometimes controversial) research in biology. While lower organisms can robustly undergo cardiac regeneration by reactivation of the embryonic cardiogenic pathway, this ability is strongly reduced in mice, both temporally and qualitatively. Finding a way to derive precursor cells with regenerative ability from differentiated cells in mammals has been challenging.

Zhou, He, and colleagues hypothesized that ISL-1-positive cells would show regenerative capacity and developed a small molecules screen to dedifferentiate cardiomyocytes (CM) to ISL1-positive precursor cells. Using hESC-derived CM, the authors found that the combination of both, WNT activation (CHIR99021) and p300 acetyltransferase inhibition (A-485) (named 2C protocol) induces CM dedifferentiation to regenerative cardiac cells (RCCs). RCCs are proliferative and re-express embryonic cardiogenic genes while decreasing the expression of more mature cardiac genes, bringing them towards a more precursor-like state. RCCs were able to differentiate into CM, smooth muscle cells, and endothelial cells, highlighting their multipotent property. In vivo, administration of 2C in rats and mice had protective effects on myocardial infarction.

Mechanistically, the authors report that the 2C protocol drives CM-specific transcriptional and epigenetic changes.

Strengths:

The authors made a great effort to validate their data using orthogonal ways, and several hESC lines. The use of lineage tracing convincingly showed a dedifferentiation from CM. They translate their findings into an in vivo model of myocardial injury, and show functional cardiac regeneration post-injury. They also showed that 2C could surprisingly be used as a preventive treatment. Together their data may suggest a regenerative effect of 2C both in vitro and in vivo settings. If confirmed, this study might unlock therapeutic strategy for cardiac regeneration.

Weaknesses:

Several points remain puzzling to me and some aspects of this study need to be clarified and extended:

General comments:

* Experimental design & Interpretation*

1) The main hypothesis (line 50) that Isl1 cells have regenerative properties is not extremely novel (10.1172/jci.insight.80920, doi.org/10.1038/nature03215,10.1016/s1534-5807(03)00363-0).

2) Based on Table S1, concentration of A-485 used in the screen is 10uM but used throughout this study at 0.5um. Could the authors provide a rationale for this 20x reduction of concentration? It would be useful to get a titration of this compound for the effects tested.

3) It is confusing to clearly understand what proportion of CMs dedifferentiate to become RCCs. The lineage tracing data suggests only 0.6%-1.5% of cells undergo this transition. It is difficult to understand how such a small fraction can have wide effects in their different experimental settings. This is specifically true when the author quantified nuclear and cytosolic area on brightfield pictures - would the same effect on nuclear/cytosolic area be observed in Isl1 KO cells?

4) The authors totally disregard the effect of i-BET-762 that gives a very similar percentage of Isl1-positive cells when combined to 2C (Supp. 1E). What is the effect of CHIR + I-BET-762 alone?

5) It is really hard to understand the contradictory effects of A-485 on acetylation status. The authors mentioned that A-485 only has an effect on H3K27Ac and not on H3K9Ac (line 221) to later (line 226) contradict themselves by saying it also has an effect on H3K9Ac. To explain this discrepancy, the authors vaguely mentioned "further analyses" without giving any other details. It would be transparent to explain what led to this radical change in interpretation.

6) The difference in the ChIP peak height is rather minimal for the H3K9Ac data. Were the peaks normalized to the sequencing depth? What does the y-axis represent on these ChIP traces (number of counts?)

7) Would it be possible to test this 2C protocol on mESC and see if similar changes occur? How transcriptionally different would these mouse RCCs be to Isl1+ progenitors isolated from neonatal mice (P1-P5)?

Statistics & Data acquisition

1) The authors mentioned experiments were done at least 3 times and each dot plotted on a graph is an average of technical repeat for one biological repeat. My understanding would be that if I see 9 dots, it means the experiment has been done 9 times - What would be the rationale for such a high number of repeats? It is an "artificial" way to increase the power of a test and might lead to misinterpretation of the data. This becomes relevant for some figures where biological difference is minimal and they still show statistical differences (e.g. Supp 2E, Supp 3A, Supp 9C,...). This is also true for in vivo section (Fig. 4G).

It would help to have a precise clarification between technical and biological repeats in the figure legends (e.g. n=3 biological repeat (aka 3 dots on a graph) obtained from averaging XX technical repeats), as well as the specific test stated the legend in addition to the general paragraph in the methods. Providing raw numerical data so readers can re-test them independently would also be a transparent way to do it.

2) Does the author test for normality before applying a specific test? Please clarify and justify either way.

3) If each dot represents a biological repeat as stated in the method section, why do some datasets have different numbers of repeats between NC and 2C if obtained in parallel? Have repeats been excluded?

Reviewer #2 (Public Review):

Summary:<br /> Starting from an AlphaFold2 model of the outward-facing conformation of the GAT1 transporter, the authors primarily use state-of-the-art MD simulations to dissect the role of the two Na+ ions that are known to be co-transported with the substrate, GABA (and a co-transported Cl- ion). The simulations indicated that Na+ binding to OF GAT depends on the electrostatic environment. The authors identify an extracellular recruiting site including residues D281 and E283 which they hypothesized to increase transport by locally increasing the available Na+ concentration and thus increasing binding of Na+ to the canonical binding sites NA1 and NA2. The charge-neutralizing double mutant D281A-E283A showed decreased binding in simulations. The authors performed GABA uptake experiments and whole-cell patch clamp experiments that taken together validated the hypothesis that the Na+ staging site is important for transport due to its role in pulling in Na+.

Detailed analysis of the MD simulations indicated that Na+ binding to NA2 has multiple structural effects: The binding site becomes more compact (reminiscent of induced fit binding) and there is some evidence that it stabilizes the outward-facing conformation.

Binding to NA1 appears to require the presence of the substrate, GABA, whose carboxylate moiety participates in Na+ binding; thus the simulations predict cooperativity between binding of GABA and Na+ binding to NA1.

Strengths:<br /> - MD simulations were used to propose a hypothesis (the existence of the staging Na+ site) and then tested with a mutant in simulations AND in experiments. This is an excellent use of simulations in combination with experiments.

- A large number of repeat MD simulations are generally able to provide a consistent picture of Na+ binding. Simulations are performed according to current best practices and different analyses illuminate the details of the molecular process from different angles.

- The role of GABA in cooperatively stabilizing Na+ binding to the NA1 site looks convincing and intriguing.

Weaknesses:<br /> - Assessing the effects of Na+ binding on the large-scale motions of the transporter is more speculative because the PCA does not clearly cover all of the conformational space and the use of an AlphaFold2 model may have introduced structural inconsistencies. For example, it is not clear if movements of the inner gate are due to an AF2 model that's not well packed or really a feature of the open outward conformation.

- Quantitative analyses are difficult with the existing data; for example, the tICA "free energy" landscape is probably not converged because unbinding events haven't been observed.

Reviewer #2 (Public Review):

Summary:<br /> In this study by Bendzunas et al, the authors show that the formation of intra-molecular disulfide bonds involving a pair of Cys residues near the catalytic HRD motif and a highly conserved T-Loop Cys with a BRSK-specific Cys at an unusual CPE motif at the end of the activation segment function as repressive regulatory mechanisms in BSK1 and 2. They observed that mutation of the CPE-Cys only, contrary to the double mutation of the pair, increases catalytic activity in vitro and drives phosphorylation of the BRSK substrate Tau in cells. Molecular modeling and molecular dynamics simulations indicate that oxidation of the CPE-Cys destabilizes a conserved salt bridge network critical for allosteric activation. The occurrence of spatially proximal Cys amino acids in diverse Ser/Thr protein kinase families suggests that disulfide-mediated control of catalytic activity may be a prevalent mechanism for regulation within the broader AMPK family. Understanding the molecular mechanisms underlying kinase regulation by redox-active Cys residues is fundamental as it appears to be widespread in signaling proteins and provides new opportunities to develop specific covalent compounds for the targeted modulation of protein kinases.

The authors demonstrate that intramolecular cysteine disulfide bonding between conserved cysteines can function as a repressing mechanism as indicated by the effect of DTT and the consequent increase in activity by BSK-1 and -2 (WT). The cause-effect relationship of why mutation of the CPE-Cys only increases catalytic activity in vitro and drives phosphorylation of the BRSK substrate Tau in cells is not clear to me. The explanation given by the authors based on molecular modeling and molecular dynamics simulations is that oxidation of the CPE-Cys (that will favor disulfide bonding) destabilizes a conserved salt bridge network critical for allosteric activation. However, no functional evidence of the impact of the salt-bridge network is provided. If you mutated the two main Cys-pairs (aE-CHRD and A-loop T+2-CPE) you lose the effect of DTT, as the disulfide pairs cannot be formed, hence no repression mechanisms take place, however when looking at individual residues I do not understand why mutating the CPE only results in the opposite effect unless it is independent of its connection with the T+2residue on the A-loop.

Strengths:<br /> This is an important and interesting study providing new knowledge in the protein kinase field with important therapeutic implications for the rationale design and development of next-generation inhibitors.

Weaknesses:<br /> There are several issues with the figures that this reviewer considers should be addressed.

Reviewer #2 (Public Review):

Summary:<br /> In this paper, the authors investigated the admixture history of domestic cattle since they were introduced into Iberia, by studying genomic data from 24 ancient samples dated to ~2000-8000 years ago and comparing them to modern breeds. They aimed to (1) test for introgression from (local) wild aurochs into domestic cattle; (2) characterize the pattern of admixture (frequency, extent, sex bias, directionality) over time; (3) test for correlation between genetic ancestry and stable isotope levels (which are indicative of ecological niche); and (4) test for the hypothesized higher aurochs ancestry in a modern breed of fighting bulls.

Strengths:<br /> Overall, this study collects valuable new data that are useful for testing interesting hypotheses, such as admixture between domestic and wild populations, and correlation between genome-wide aurochs ancestry and aggressiveness.

Weaknesses:<br /> Most conclusions are partially supported by the data presented. The presence of admixed individuals in prehistorical periods supports the hypothesized introgression, although this conclusion needs to be strengthened with an analysis of potential contamination. The frequency, sex-bias, and directionality of admixture remain highly uncertain due to limitations of the data or issues with the analysis. There is considerable overlap in stable isotope values between domestic and wild groups, indicating a shared ecological niche, but variation in classification criteria for domestic vs wild groups and in skeletal elements sampled for measurements significantly weakens this claim. Lastly, the authors presented convincing evidence for relatively constant aurochs ancestry across all modern breeds, including the Lidia breed which has been bred for aggressiveness for centuries. My specific concerns are outlined below.

Contamination is a common concern for all ancient DNA studies. Contamination by modern samples is perhaps unlikely for this specific study of ancient cattle, but there is still the possibility of cross-sample contamination. The authors should estimate and report contamination estimates for each sample (based on coverage of autosomes and sex chromosomes, or heterozygosity of Y or MT DNA). Such contamination estimates are particularly important to support the presence of individuals with admixed ancestry, as a domestic sample contaminated with a wild sample (or vice versa) could appear as an admixed individual.

A major limitation of this study is uncertainty in the "population identity" for most sampled individuals (i.e., whether an individual belonged to the domesticated or wild herd when they were alive). Based on chronology, morphology, and genetic data, it is clear the Mesolithic samples from the Artusia and Mendandia sites are bona fide aurochs, but the identities of individuals from the other two sites are much less certain. Indeed, archeological and morphological evidence from El Portalon supports the presence of both domestic animals and wild aurochs, which is echoed by the inter-individual heterogeneity in genetic ancestry. Based on results shown in Fig 1C and Fig 2 it seems that individuals moo017, moo020, and possibly moo012a are likely wild aurochs that had been hunted and brought back to the site by humans. Although the presence of individuals (e.g., moo050, moo019) that can only be explained by two-source models strongly supports that interbreeding happened (if cross-contamination is ruled out), it is unclear whether these admixed individuals were raised in the domestic population or lived in the wild population and hunted.

Such uncertainty in "population identity" limits the authors' ability to make conclusions regarding the frequency, sex bias, and directionality of gene flow between domestic and wild populations. For instance, the wide range of ancestry estimates in Neolithic and Chalcolithic samples could be interpreted as evidence of (1) frequent recent gene flow or (2) mixed practices of herding and hunting and less frequent gene flow. Similarly, the statement about "bidirection introgression" (on pages 8 and 11) is not directly supported by data. As the genomic, morphological, and isotope data cannot confidently classify an individual as belonging to the domesticated or wild population, it seems impossible to conclude the direction of gene flow (if by "bidirection introgression" the authors mean something other than "bidirectional gene flow", they need to clearly explain this before reaching the conclusion.)

The f4 statistics shown in Fig 3B are insufficient to support the claim regarding sex-biased hybridization, as the f4 statistic values are not directly comparable between the X chromosome and autosomes. Because the effective population size is different for the X chromosome and autosomes (roughly 3:4 for populations with equal numbers of males and females), the expected amount of drift is different, hence the fraction of allele sharing (f4) is expected to be different. In fact, the observation that moo004 whose autosomal genome can be modeled as 100% domestic ancestry still shows a higher f4 value for the X chromosome than autosomes hints at this issue. A more robust metric to test for sex-biased admixture is the admixture proportion itself, which can be estimated by qpAdm or f4-ratio (see Patterson et al 2012). However, even with this method, criticism has been raised (e.g., Lazaridis and Reich 2017; Pfennig and Lachance, 2023). In general, detecting sex-bias admixture is a tough problem.

In general, the stable isotope analysis seems to be very underpowered, due to the issues of variation in classification criteria and skeletal sampling location discussed by the authors in supplementary material. The authors claimed a significant difference in stable nitrogen isotope between (inconsistently defined) domestic cattle and wild aurochs, but no figures or statistics are presented to support this claim. Please describe the statistical method used and the corresponding p-values. The authors can consider including a figure to better show the stable isotope results.

Reviewer #2 (Public Review):

Summary:<br /> Wang and collaborators have evaluated the impact of inflammation on bone loss induced by Doxorubicin, which is commonly used in chemotherapy to treat various cancers. In mice, they show that a single injection of Doxorubicin induces systemic inflammation, leukopenia, and significant bone loss associated with increased bone-resorbing osteoclast numbers. In vitro, the authors show that Doxorubicin activates the AIM2 and NLRP3 inflammasomes in macrophages and neutrophils. Importantly, they show that the full knockouts (germline deletions) of AIM2 (Aim2-/-) and NLRP3 (Nlrp3-/-) and Caspase 1 (Casp1-/-) limit (but do not completely abolish) bone loss induced 4 weeks after a single injection of Doxorubicin in mice. From these results, they conclude that Doxorubicin activates inflammasomes to cause inflammation-associated bone loss.

Strengths:<br /> While numerous studies have reported that Doxorubicin activates the inflammasome in myeloid cells and various other cell types, that Doxorubicin induces systemic inflammation, and that both the systemic inflammation and Doxorubicin treatment leads to bone loss, functional experiments demonstrating that NRLP3 and/or AIM2 loss-of-functions, and thus the systemic impairment of the inflammatory response, may prevent bone-loss induced by Doxorubicin were lacking. The strength of this manuscript is that it provides these missing data.

Weaknesses:<br /> However, one could argue that most of the conclusions drawn from the data presented here have been previously reported and that it was very much expected that reducing systemic inflammation in treated animals (in Aim2-/- and/or Nlrp3-/- mice) would preserve bone homeostasis to some extent, similarly to what has been reported in the context of cardiotoxicity induced by Doxorubicin.

Since the manuscript focuses on therapeutic considerations aiming to preserve bone homeostasis in animals treated with Doxorubicin, additional experiments evaluating and comparing various therapeutic options could improve the impact of the study. Drugs targeting the inflammasomes could be tested in addition to the genetic mouse models. Since increased osteoclast numbers (and likely bone resorption) are associated with Doxorubicin-induced bone loss, antiresorptive drugs such as Bisphosphonates or anti-RANKL antibodies could be tested and compared to anti-inflammatory drugs. Since autophagy and senescence have been shown to contribute to bone loss induced by Doxorubicin, it would be interesting to use the pharmacologic inhibitors (targeting autophagy or senescence) used in these previous studies to evaluate the relative impact of these different cellular mechanisms, on bone loss induced by Doxorubicin.

Moreover, the cellular and molecular mechanisms by which Doxorubicin induces bone loss in vivo could be further evaluated. Doxorubicin has been reported to directly affect bone-making osteoblasts and bone-resorbing osteoclasts. It would be important to determine the relative importance of the activation of the AIM2 and NLRP3 inflammasomes in these cells compared to macrophages and neutrophils. Floxed mouse lines exist for both Aim2 and Nlrp3, as well as relevant cell-specific Cre lines. Thus, cell-specific conditional knockouts could have been used in the current study, instead of using global knockout animals. Genetic tools also exist to induce the specific ablation of macrophages or neutrophils and could be used. Furthermore, it is unclear whether local inflammation is induced in the bone marrow of Doxorubicin-treated mice, and what is the relative impact of local versus systemic inflammation in bone loss in these mice. Markers of the inflammasomes, pyroptosis, and NETosis could be evaluated on bone sections, and on bone and bone marrow samples. The effect of Doxorubicin on osteoblast numbers in vivo and on bone resorption (not just osteoclast numbers) should be evaluated as well. These mechanistic aspects are important and needed to better understand the cytotoxic mechanisms triggered by Doxorubicin, and define the best therapeutic approaches to preserve bone integrity in chemotherapy.

Finally, it would be important to assess the bone mass of Doxorubicin-treated control, Aim2-/-, Nlrp3-/- and Cas1-/- mice at a later time point than 4 weeks post-injection. Nlrp3 knockout has been reported to increase the density of the cortical and trabecular bones. The bone mass of Aim2-/-, Nlrp3-/- and Cas1-/- mice at baseline may be higher than that of control mice, and it may take slightly longer for Doxorubicin to reduce bone mass to the same extent than in controls. It would be also interesting to do similar experiments using animals treated multiple times with Doxorubicin instead of using a single injection, since patients receive their chemotherapy multiple times.