1. Sep 2026
    1. Note: This response was posted by the corresponding author to Review Commons. The content has not been altered except for formatting.

      Learn more at Review Commons


      Reply to the reviewers

      Manuscript number: RC-2026-03654

      Corresponding author(s): Yusuke, Kishi

      1. General Statements [optional]

      Reviewer #1 (Evidence, reproducibility and clarity (Required))

      The manuscript by You et al. investigates changes in gene expression and histone modifications after juvenile social isolation (jSI) in the nucleus accumbens (NAc). They find many differentially expressed genes and find overlap with activating or repressive histone marks. They then go on to compare their data to other published datasets to support their findings. This is an interesting study, and the authors use creative approaches using unique and published data to identify epigenetic mechanisms underlying the changes in gene expression within the NAc. However, there are many points of clarification that are needed to fully evaluate the manuscript and several experimental details.

      Reviewer #1 (Significance (Required))

      This is an interesting study and the authors use creative approaches using unique and published data to identify epigenetic mechanisms underlying the changes in gene expression within the NAc. However there are many points of clarification that are needed to fully evaluate the manuscript and several experimental details are missing or unclear.

      We thank the reviewer for the positive assessment of our study and for the constructive comments, which have helped us to improve the manuscript considerably. The reviewer's points centered in particular on the framing of the Introduction, especially the conflation of adult and adolescent social isolation, the use of nominal rather than FDR-corrected p-values to define DEGs, and the lack of clarity in several parts of the Methods.

      We have addressed each of these points in our responses below, in part through revisions already made to the manuscript. Briefly, the Introduction and the Methods have been substantially revised, and we will carry out additional threshold-free analyses to support our conclusions. We are grateful to the reviewer for identifying the points that most needed clarification.

      Reviewer #2 (Evidence, reproducibility and clarity (Required))

      This paper profiles NeuN positive NAc nuclei after juvenile social isolation. The authors report RNA seq changes and CUT and Tag for H3K4me1, H3K4me3, H3K27ac, and H3K27me3. They then link DEGs to public datasets for Kdm6b, Brd4, and Setd1a. The neuron enriched design is useful. The main claim remains correlative. Causal support is thin.

      Reviewer #2 (Significance (Required))

      This study provides a useful neuron-enriched transcriptomic and histone modification resource from the nucleus accumbens following juvenile social isolation. The integration of RNA-seq and CUT&Tag data adds value for researchers studying epigenetic regulation and stress-related neurobiology. However, the advance is primarily descriptive rather than mechanistic, as the conclusions rely largely on correlative analyses without functional validation. The manuscript will be of interest to the neuroepigenetics and psychiatric neuroscience communities, but the conceptual advance is incremental, and the mechanistic claims should be moderated.

      We thank the reviewer for recognizing the value of our neuron-enriched dataset and for the constructive comments, which have helped us to improve the manuscript considerably. The reviewer's points centered in particular on the correlative nature of our findings and the need to moderate causal language, the permissive thresholds used to define DEGs and the absence of correction for multiple testing, and the interpretation of the comparisons with published datasets.

      We address each of these points in detail in our responses below, and have already implemented a substantial part of the revisions. Briefly, we have removed wording implying that the identified epigenetic factors mediate jSI-induced transcriptional changes and now describe these relationships as associations that remain to be tested, and the comparisons with published datasets are explicitly framed as exploratory. We will also carry out threshold-free analyses and apply correction for multiple testing. We are grateful to the reviewer for these suggestions, which we believe have made the manuscript more accurate about what our data can and cannot support.

      Reviewer #3 (Evidence, reproducibility and clarity (Required))

      Summary:

      In this study, the authors were investigating the effect of juvenile social isolation (jSI) on the nucleus accumbens (NAc) transcriptome in female mice. They used P21 wild type (C57BL6) female, isolated at P21 or group housed, and reunited at P35 to generate their samples for RNAseq on NAc punches. They also performed FACS sorting on the nuclei from NAc lysates to select only neuronal nuclei (NeuN staining). They studied the differentially expressed genes (DEGs) between group housed and jSI by RNAseq. Then, they studied the protein/protein interaction via stringDB on the DEGs identified (up or down) and perform GO analysis on them. They identified Ntrk2, Grin3a, Grik1 and Bcl2; associated with the neuronal function or transcription regulation terms. They also studied the histones modifications (H3K4me1, H3K4me3, H3K27ac, and H3K27me3) after jSI and identified neuronal function and transcription regulation terms again on the DDR (Cut and Tag method). They found that these histones modifications could play a role in jSI-induced adaptations and neuronal function. Finally, they reanalyzed public datasets of RNAseq data to identify histone modifications associated with their DEGs of interest, and compare their DEGs to differences between genotypes in the public datasets (in conditional KO models of some genes of interest, as Kdm6b cKO, BET inhibitor, or Setd1a +/- mice with 3 different mutations. They conclude that histone modification could be involved in jSI-induced gene expression alteration.

      Reviewer #3 (Significance (Required))

      Globally, this study is showing transcriptomic and histone modifications occurring after juvenile social isolation in female mice. The authors identified differentially expressed genes linked to neuronal development, regulation of transcription and chromatin remodeling. The reanalyzed public datasets to identify histone modification on the DEG identified and observed the impact of some already published mutations on the gene expression to compare it to their data. The limits of this issue are caused by the reanalysis parts, since the datasets used are cortical and cerebellum samples, in diverse development stage (embryonic, early juvenile, adults - both sexes) that is quite different compared to their paradigm (juvenile females). They also never display the name of the principal DEGs identified (text or plots) which leads to difficult understanding of the findings of this paper. A more focused analysis on NAc or striatal, female only, juvenile stage datasets would be more helpful in this situation. The text should be more precise sometimes and a specific explanation on the exclusion of male mice should be introduce early in the methods.

      This study finds its place in the current research on the role of NAc function/dysfunction in behavioral abnormalities induced by social isolation and try to understand the mechanisms behind the abnormal behaviors induced by separation. The audience could be composed of researchers from several domains where social isolation is the cause or the consequence of pathological behaviors, including studies on loss, depression, ASD, Alzheimer, Schizophrenia, etc). The context is quite broad. The results from this paper could help find new molecular targets to alleviate the effects of social isolation and perhaps ameliorate the behavior for mouse models of several diseases or later in patients. A better understanding of the effects of social isolation in female is interesting, but being able to compare both sexes would be even better: identifying sex-differences and common defect is of great interest nowadays in several domains.

      The present reviewer has expertise in behavior in mice (both male and female) from juvenile to adult stages, has studied neuronal circuits including prefrontal cortex and striatum (mainly NAc) in behavioral abnormalities in mice in a model of ASD and more recently in an addiction model. The reviewer is interested particularly in sex-differences in neuronal circuits defects and behavior expression in diseases. Finally, the reviewer has recently focused on spatial transcriptomic approaches in addiction models.

      We thank the reviewer for the detailed and constructive assessment of our study, and for the many specific suggestions, which have helped us to improve the manuscript considerably. The reviewer's points centered in particular on the discrepancy between the published datasets we reanalyzed and our own experimental paradigm, the absence of gene names in the text and figures, which made our findings difficult to follow, and the need for greater precision in the text, including an explicit explanation of why only female mice were used.

      We address each of these points in detail in our responses below, and have already implemented a substantial part of the revisions. Briefly, the rationale for using female mice is now stated in both the Introduction and the Methods, the text has been revised throughout for precision, and the comparisons with published datasets are explicitly framed as exploratory. We will also reanalyze the striatal dataset from Chen et al. (Sci. Adv., 2020), which is considerably closer to our samples than the prefrontal cortex datasets, and we will label the relevant genes in the text and figures. We are grateful to the reviewer for the care taken in reading the manuscript.

      2. Description of the planned revisions

      1-2

      The authors cite several studies from the Nestler lab on early life stress that link early life stress to histone modifications but failed to cite the manuscripts that investigated the transcriptional changes in response to jSI. These studies also highlight sex differences in jSI. This is important given that this study only uses females. Many of the effects described might not be comparable simply because of the sex of the animals.

      We thank the reviewer for this helpful comment. We agree that studies investigating transcriptional changes in response to jSI, including those reporting sex differences, should be cited, and we will add these references in the revised Introduction.

      The paragraph highlighted by the reviewer (the third paragraph of the Introduction) was intended to summarize the relationship between stress and epigenetic regulation in the NAc, which is why studies from the Nestler laboratory are prominently represented; transcriptional responses to social isolation in the NAc were summarized in the preceding paragraph. As the reviewer notes, however, other important studies have since been reported, and we will incorporate them accordingly.

      We also thank the reviewer for raising the issue of sex differences, which has been highlighted repeatedly by the other reviewers as well and is indeed an important point. We have revised the Introduction to specify the sex of the animals used in the rodent studies, and to make the reported sex differences explicit (lines 75–78).

      1-3

      Can the authors clarify if they used an adjusted p-value or nominal p-value. If they are using a nominal p-value the authors should explain their reasoning and provide information regarding if any of the transcripts survived a p-value correction. The addition of threshold-free approaches are more appropriate (GSEA) rather than focusing on transcripts with a nominal p-value. If they are going to present data using a nominal p-value, this should be justified and the cut off should be explained and every interpretation should include a caveat.

      We thank the reviewer for this important comment. The reviewer is correct that nominal p-values, rather than adjusted p-values (FDR), were used in this study, and we agree that this point requires both clarification and revision.

      In the revised manuscript, we will take the following steps. First, as the reviewer suggests, we will perform GSEA as a threshold-free approach and replace the current GO analysis of DEGs with the GSEA results. Second, for analyses that necessarily require a defined gene set, namely the PPI analysis, the ChIP-Atlas analysis, and the overlap analyses shown as Venn diagrams, we will continue to use gene sets defined by nominal p-values, but we will state explicitly in the text that these analyses are exploratory and hypothesis-generating rather than confirmatory.

      Regarding our original choice of threshold, we note that the number of biological replicates in this study (n = 4–6) is small relative to the number of genes tested, and that transcriptional responses to stress in the nervous system are typically of small effect size with substantial inter-individual variability. Under these conditions, FDR correction is highly conservative, and previous studies in this field have reported findings based on nominal p-values (e.g. Torres-Berrío et al., Nat. Neurosci., 2024). Since the aim of the present study was to generate hypotheses rather than to establish definitive causal factors, we adopted a permissive threshold. We will make this rationale explicit in the revised Methods and Results.

      Finally, we will add a statement to the Limitations section noting that the findings reported here require validation in future work. We hope that these revisions adequately address the reviewer's concern.

      1-4

      It is unclear how the authors confirmed that input RNA or neurons were similar across samples for the library prep. This is especially important given the top genes that are differentially expression. The finding that beta actin (Actb) is up in jSI vs GH animals. The results could be due to differences in input rather than actual differences in expression.

      We thank the reviewer for raising this important point. We understand the concern to be that the observed change in Actb, a gene commonly regarded as a housekeeping gene, may reflect differences in input material between samples rather than a genuine difference in expression.

      We would like to clarify that the same number of sorted nuclei was used for every sample, and we have stated this explicitly in the Methods section of the revised manuscript. Consistent with this, library quality assessed by fluorometry (Qubit) and capillary electrophoresis (TapeStation) was comparable across all samples, and the FACS profiles showed the similar pattern in every case; we will present these data in the point-by-point response accompanying the revised manuscript. In addition, we will report the TMM normalization factors calculated in edgeR, together with PCA and/or MDS analyses, to confirm that no sample deviated substantially from the others and that normalization was applied appropriately.

      Regarding housekeeping gene expression, Actb is known to undergo activity-dependent changes in expression in neurons and to contribute to synaptic function and plasticity. We consider this an interesting observation in its own right, particularly as another reviewer noted that the “actin reorganization” GO term among our differentially expressed genes warrants further description. In the revised manuscript, we will cite the relevant literature and discuss this point. We will also confirm that the expression of other housekeeping genes was largely unchanged, supporting the consistency of RNA-seq quality across samples.

      1-5

      There are aspects of the methods are difficult to understand. For example, under RNA-seq, the authors mention "frozen nuclei were thawed and centrifuged.....the supernatant was discarded and nuclei were centrifuged again under the same conditions" Can the authors clarify what was done here? Were the nuclei resuspended in STEM CellBANKER or something else? While this is a concrete example, there are many other places in the methods the are like this, meaning that steps seem to be skipped and it then becomes difficult to assess the approach. It is recommended that the authors work to clarify the methods. Another example, how the DNA was treated in the CUT&TAG and how much DNA was added to the library prep.

      We apologize for the lack of clarity in the Methods section. We will review and revise the entire Methods section to ensure that each step is described unambiguously.

      Regarding the specific example raised by the reviewer, sorted nuclei were resuspended in STEM CELLBANKER and frozen in this solution; after thawing, they were centrifuged directly without resuspension in any other buffer. We will make this explicit in the revised text.

      For the CUT&Tag experiments, the same number of nuclei was used for every sample, and we have stated this explicitly in the revised Methods. Regarding the amount of DNA used for library preparation, we did not quantify the DNA prior to PCR amplification, as the yield at this step is below the range of reliable quantification and measurement is not included in the original CUT&Tag protocol. Instead, all libraries were amplified with the same number of PCR cycles (12 cycles, as stated in the Methods), ensuring that they were prepared under identical conditions throughout.

      2-major-2

      DEG thresholds are loose. DEGs are defined as "p-value 1.2." There is no clear FDR cutoff. With ~1250 DEGs from n = 5 to 6, many hits may be noise. Please report FDR filtered lists or justify the uncorrected p value choice. Re run key GO and overlap tests on a stricter set.

      We thank the reviewer for this important comment, and we agree that the thresholds used to define DEGs require clarification and revision.

      In the revised manuscript, we will take the following steps. First, we will perform GSEA as a threshold-free approach and replace the current GO analysis of DEGs with the GSEA results, so that our functional conclusions do not depend on an arbitrary cutoff. Second, for analyses that necessarily require a defined gene set, namely the PPI analysis, the ChIP-Atlas analysis, and the overlap analyses shown as Venn diagrams, we will continue to use gene sets defined by nominal p-values, but we will state explicitly in the text that these analyses are exploratory and hypothesis-generating rather than confirmatory.

      Regarding our original choice of threshold, the number of biological replicates in this study (n = 4–6) is small relative to the number of genes tested, and transcriptional responses to stress in the nervous system are typically of small effect size with substantial inter-individual variability. Under these conditions, FDR correction is highly conservative, and previous studies in this field have reported findings based on nominal p-values (e.g. Torres-Berrío et al., Nat. Neurosci., 2024). Since the aim of the present study was to generate hypotheses rather than to establish definitive causal factors, we adopted a permissive threshold. We will make this rationale explicit in the revised Methods and Results, and we will add a statement to the Limitations section noting that the findings reported here require validation in future work.

      2-major-3

      Public data overlaps are hard to interpret. Kdm6b data are from cerebellum. Brd4 data are from cultured cortical neurons treated with JQ1. Setd1a data are mostly PFC. The authors note that "different brain regions, cell types, and experimental conditions... may contribute to... false negative or false positive results." That caveat is important. Overlaps should be framed as hypothesis generating only. Do not treat them as evidence that these enzymes act in NAc under jSI.

      We fully agree with the reviewer on this point, and we have revised the manuscript accordingly.

      First, we have explicitly framed all comparisons with published datasets as exploratory and hypothesis-generating, and we have removed any wording implying that these enzymes act as mediators of jSI-induced transcriptional changes in the NAc.

      Second, as pointed out by Reviewer #3, the study by Chen et al. (Sci. Adv., 2022) also performed scRNA-seq using the striatum of Setd1a heterozygous mice. Since striatal tissue is far closer to our NAc samples than the prefrontal cortex, we will reanalyze the striatal dataset and compare it with our jSI DEGs. Depending on the outcome of this analysis, we will reorganize the figures so that the most relevant comparison is presented in the main text and the remaining reanalyses are moved to the supplementary material.

      Third, we have explained our rationale for dataset selection in the revised text. The three factors examined here were nominated by our own data: Brd4 and Setd1a emerged from the ChIP-Atlas promoter analysis, and Kdm6b was itself upregulated in our RNA-seq data. We then searched for publicly available RNA-seq datasets in which these factors had been perturbed in the nervous system. No such dataset exists for the NAc or striatum for Kdm6b or Brd4, and we therefore selected the datasets that were closest to our system among those available. We have stated this limitation in the Limitation section, together with the differences in brain region, cell type, developmental stage, and sex between these datasets and our own.

      2-major-5

      CUT and Tag analysis is coarse for promoter claims. Signals are quantified in "all 5 kbp bins." That bin size can blur promoters, enhancers, and neighboring genes. Please add peak calling or TSS centered analyses for key loci such as Grik1, Bcl2, and Dkk3. Also show more browser tracks beyond one example.

      We thank the reviewer for this comment, and we agree that quantification in 5 kbp bins is too coarse to support claims about promoter-level regulation.

      In the revised manuscript, we will perform higher-resolution analyses centered on transcription start sites for the key loci highlighted by the reviewer, including Grik1, Bcl2, and Dkk3, so that promoter signals can be evaluated separately from those of surrounding regions. We will also try to perform peak calling to define regions of enrichment more precisely.

      We will also increase the number of browser tracks shown, so that examples are provided for each histone modification rather than for a single locus.

      2-major-6

      Multiple testing for overlaps needs attention. Many Fisher tests compare DEGs with DDRs and with several public DEG lists. Report whether p values were corrected across tests. Some reported overlaps are small in absolute numbers even when p values look significant.

      We thank the reviewer for this comment. The reviewer is correct that we did not correct for multiple testing across the repeated Fisher's exact tests.

      In the revised manuscript, we will apply the Benjamini-Hochberg procedure and report adjusted p-values in the figures. Correction will be performed within each analysis category rather than across all tests, namely the overlaps between DEGs and DDRs, and the overlaps between our jSI DEGs and each of the published datasets for Kdm6b, Brd4, and Setd1a. We note that these tests are not fully independent, since they share the jSI DEG list and involve mutually exclusive up- and down-regulated gene sets; the Benjamini-Hochberg procedure remains valid under such positive dependence.

      We also agree that p-values alone can be misleading when the absolute number of overlapping genes is small. For every comparison, we will additionally report the observed and expected numbers of overlapping genes together with the odds ratio or fold enrichment, so that the magnitude of each overlap can be assessed independently of its p-value.

      2-minor-2

      Down DEG GO terms include "Chondrocyte differentiation" and "Positive regulation of cartilage development." These look odd for NAc neurons. Check annotation quality and whether these terms survive stricter DEG filters.

      We thank the reviewer for this observation. Genes involved in developmental processes are frequently shared across tissues, and GO annotation assigns such pleiotropic genes to multiple terms, which can produce enrichment for categories that appear unrelated to the tissue under study.

      In our case, the genes driving the enrichment of "chondrocyte differentiation" and "positive regulation of cartilage development" include Sox5, Bmp4, and Nfib, all of which have established roles in nervous system development. Whether these genes are functionally important in the NAc remains unknown, but their appearance in these terms reflects annotation overlap between chondrocyte and neural developmental programs rather than an implausible result.

      We nevertheless acknowledge the concern regarding our DEG thresholds, as discussed in our response to the Reviewer #2's comment (#2-major-2). We will perform GSEA as a threshold-free approach and use it to confirm the functional categories identified by the current GO analysis.

      3-major-5(OPTIONAL)

      Datasets selection: The public datasets used through this study are far from the original experimental design proposed in this study (cerebellum at P14, cortex at E16.5, nonspecific inhibitor, whole adult PFC = 12-14weeks old). It is important to note that the dataset used, from Chen et al 2022 (whole PFC) also studied the striatum of heterozygous Setd1a mice in the same paper. Why did the author reanalyzed PFC data instead of striatum, which would probably look more like their NAc-restricted samples? Restrict the analysis of the public dataset on NAc data, if possible on females only, and/or try to obtain data from similar experimental design (social isolation, stressed mice). Note here, that the development stage during which the mice have been isolated/regrouped and sample taken will probably of importance. The reanalyzes mix embryonic, early juvenile and adult samples, none of which is consistent nor look like their set up (P31-35).

      We thank the reviewer for this thoughtful comment, and we agree that the discrepancy between the published datasets and our own experimental design is a substantial limitation.

      Regarding the Chen et al. (Sci. Adv., 2022) dataset, we are grateful to the reviewer for pointing out that the same study also profiled the striatum. Since the NAc is part of the striatum, this dataset is considerably closer to our samples than the prefrontal cortex, and we will reanalyze it and compare it with our jSI DEGs in the revised manuscript.

      We will also explain our rationale for dataset selection explicitly. The three factors examined here were nominated by our own data: Brd4 and Setd1a emerged from the ChIP-Atlas promoter analysis, and Kdm6b was itself upregulated in our RNA-seq data. We then searched for publicly available RNA-seq datasets in which these factors had been perturbed in the nervous system. To our knowledge, no dataset combines perturbation of these enzymes with the NAc or striatum (except for the Setd1a striatal dataset noted above), with female animals only, or with a social isolation or stress paradigm. Among the datasets available, we therefore selected those closest to our system, prioritizing perturbation of the factor of interest, since this was the specific question the analysis was designed to address. We will state this constraint clearly in the revised text, together with the differences in brain region, cell type, developmental stage, and sex between these datasets and our own.

      Finally, in line with the comments from this reviewer and from Reviewer #2 (#2-major-3), we will frame all of these comparisons as exploratory and hypothesis-generating, and will remove wording implying that these enzymes mediate jSI-induced transcriptional changes.

      3-major-7

      Assumptions are made based on 2 sets of reanalyses. These parts should be displayed in the supplementary to help the authors target some genes of interest rather than the principal figures. These analyses didn't seem convincing due to too much shift from the original issue of the paper (which is jSI in the NAc in female mice).

      We thank the reviewer for this comment, and we agree that the reanalyses of published datasets are exploratory in nature and are considerably removed from the central question of this study.

      As described in our response to the Reviewer #3's comment (#3-major-5), we will reanalyze the striatal dataset from Chen et al. (Sci. Adv., 2022), which is far closer to our NAc samples than the prefrontal cortex datasets used previously. Depending on the outcome of this analysis, we will reorganize the figures so that the most relevant comparison is retained in the main text and the remaining reanalyses are moved to the supplementary material.

      Throughout the revised manuscript, we will present these comparisons explicitly as a means of narrowing down candidate genes for future investigation, rather than as evidence that these enzymes act in the NAc under jSI. We have also added a statement at the beginning of this section noting that the datasets were obtained under conditions different from ours, and the specific differences will be described in the Limitations section.

      3-major-8

      Volcano plots throughout the study: Should display the genes names (at least top 10 up and top 10 down DEGs) on the graph, otherwise the volcano plots are unreadable.

      Thank you for your comment. We will present our top annotated genes in the figure.

      3-minor-intro-10

      General comment: Since the paper is focused on female, it could be of interest to state in the introduction if/how sex differences exist in relation to social isolation and human diseases showing isolation as a phenotype.

      We thank the reviewer for this suggestion, which we have adopted. We have added a statement to the Introduction describing what is known about sex differences in the effects of jSI in rodents, noting that while some outcomes are shared between sexes, others such as sociability and aggression differ. We will also describe what is known about sex differences in the human conditions in which isolation or loneliness features as a phenotype. We have also specified the sex of the animals used in the rodent studies we cite, as requested by Reviewer #1 (#1-2). We have kept this description focused on social isolation rather than surveying sex differences in psychiatric disease more broadly, so that it remains relevant to the present study.

      3-minor-results-1

      Section 1: Only 1 or 2 GO terms (down / up DEGs) are described, but top5 is represented in the figure. Description of the others would be of interest (for example actin reorganization could be particularly interesting). Also, citing some DEGs from the top UP and DOWN, representative of the GO terms, could be of huge interest here.

      We thank the reviewer for this comment. We agree that describing the GO terms in more detail would improve the readability of this section. In the revised manuscript, we will describe all of the top-ranked GO terms shown in the figure, rather than only one or two, including terms such as actin reorganization. We will also name representative differentially expressed genes belonging to these terms in the text, so that the reader can appreciate which genes underlie each enrichment without consulting the supplementary tables.

      3-minor-results-2

      Figure 1, E/F/G: Some genes in Fig1G are not from top10 nodes in DEGs (Slc17a8, Hcrtr2, Dkk3, Dact1, Nr4a1, Fosl2, Htr5a). They are stated as "potentially important genes" in the legends and are found later in the study as important using other methods than RNAseq. Since Fig1G is displaying RNAseq results, it would be better to stick to the Top10 genes displayed here and put in another figure the other "potentially important genes".

      What means "potentially important"? Why these? What are the criterion? Also, the fig1G is unclear visually: separate it in two for top10 down and top10 up, it would be easier to understand and navigate.

      Legends: Statistics used for 1G are not stated.

      Volcano plot: Top 10 genes UP/down could be displayed on the graph.

      We thank the reviewer for these comments.

      We agree that including genes in Fig. 1G that were not among the top nodes made the figure difficult to follow. These genes will be moved to the section corresponding to the previous Fig. 5, where they are first identified as candidates, and Fig. 1G will show only the top node genes.

      We also agree that the term "potentially important genes" was unclear. Since this panel shows the top nodes identified by the PPI analysis, we have replaced this wording with "top nodes" in the figure legend.

      Figure 1G will be divided into separate panels for up-regulated and down-regulated DEGs, as suggested.

      We will state the statistical method used in the figure legend, and we will label the top 10 genes on the volcano plot.

      3-minor-results-4

      Figure 2: A/B/C: only one mention of the NAc in the data; Fig D: 5/8 NAc datasets. The analysis here seems unbalanced. Why not take into account only NAc datasets? The composition of cortical area or retina is highly different than the NAc (mostly Glutamatergic vs GABAergic populations). If doable, the analysis focused on NAc datasets would be better.

      We thank the reviewer for this suggestion, which we agree would improve the specificity of this analysis.

      For the histone modification analysis, a sufficient number of NAc-derived datasets is available in ChIP-Atlas to support the analysis on its own, and we will therefore repeat this part using only NAc datasets in the revised manuscript.

      For the transcription factor analysis, however, the number of NAc-derived ChIP-seq datasets is too small for a comparable enrichment analysis, and restricting the analysis in this way would leave most candidate factors untested. We will therefore retain the "Neural" category for this part, and we will state explicitly that the underlying datasets derive from a range of neural tissues whose cellular composition differs from that of the NAc, as the reviewer notes. As described in our response to the Reviewer #3's comment (#3-minor-results-3), we will also make clear that this analysis was intended to generate candidates rather than to identify regulatory relationships operating in the NAc.

      3-minor-results-5

      Section 3: "First, neuronal development-related genes, such as 'nervous system development', were found in all DDRs of the four histone modifications" line 193-194: sentence is unclear, the author probably meant "term". No gene have been cited here in any histone modification experiment (nor visible in the figure, only dots without names, top 10 up/down could be displayed on the volcano plot). It would be of interest to state at least some of the genes identified here (DDRs, closest loci) and to see if/how many common genes from RNAseq data were found again here. GO terms: again, only one or two examples are described, but figure shows the top5. They all could be at least stated.

      The conclusion of the first paragraph states: "These results were consistent with the transcriptome analysis that neuronal function and transcription-related genes were affected, and with the transcription factor analysis that epigenetic regulators were predicted to bind to the promoter regions of these genes." line 197-199: Since the author did not state any genes, we can only believe that the result are consistent based on two vague GO terms "neuronal system development" and "regulation of transcription by RNApol II/chromatin remodeling". If the lector has to read itself every gene table to know which genes are dysregulated in jSI, this study will be really time consuming.

      We thank the reviewer for these comments, and we apologize that the description of "nervous system development" was inaccurate. We have rewritten this sentence so that it refers to genes functionally related to this GO term being enriched among the DDRs, rather than to the term itself being found among the DDRs (lines 218–219).

      We will also describe all of the top-ranked GO terms shown in the figure, rather than only one or two, and we will name the genes associated with the DDRs both in the text and on the volcano plots. In addition, we will state how many of these genes overlap with the DEGs identified in our RNA-seq analysis, so that the correspondence between the two datasets is apparent without consulting the supplementary tables.

      3-minor-results-6

      Figure 3: Volcano could display the top10 names of DDRs.

      "The results indicated that down-DEGs were associated with H3K4me1, H3K4me3, and H3K27ac." line 203: In which direction are altered H3K4me1, H3K4me3, and H3K27ac? This is important to know.

      "Consistent with their active roles in transcription, downregulation of H3K4me1 and H3K27ac was more relevant to down-DEGs than up-DEGs." line 205: Why? Unclear statement.

      "Considering the composite roles of H3K4me3 (an active histone modification) and H3K27me3 (a repressive histone modification), we hypothesize a major role of H3K27me3 in these up-DEGs, and the contribution of H3K4me3 to gene expression alteration by jSI might be small, though we cannot exclude the possibility that it regulates certain genes locally or plays a repressive role." line 208-211: It is very unclear here, why H3K27me3 should play a major role while H3K4me3 alteration "might be small". This has to be further discussed.

      "For top 10 nodes among up-DEGs, we didn't find any significant alterations in any of the four histone modifications around their gene loci, except for downregulated H3K4me3 around Aldh18a1, Lamp1, and Gnb4" line 217-219: Formulation is clumsy here, reformulate.

      "Some of these genes were marked by multiple altered histone modifications." Line 223: Which ones? Only Bcl2 displayed, but the authors state "some of these genes" right after writing "H3K27ac was found to be downregulated around Grin3a, Grik1, and Adgre1". Are these the other genes showing several histone modifications? It is unclear.

      We thank the reviewer for these comments, which have helped us to clarify this section.

      Regarding the direction of the histone modification changes, we have revised the text to state explicitly in which direction each modification was altered, rather than referring only to an association.

      Regarding the roles of H3K4me3 and H3K27me3 in up-DEGs, we agree that our reasoning was not adequately explained. We have rewritten this passage to make the logic explicit: the reduction of the repressive mark H3K27me3 is consistent with the upregulation of these genes, whereas the concurrent reduction of the active mark H3K4me3 is not, which is why we consider H3K27me3 the more likely contributor at these loci (lines 237–241).

      We have also rewritten the sentence describing the top 10 nodes among up-DEGs, which was awkwardly constructed, so that it now states positively which genes showed an alteration and in which direction (lines 243–253). Similarly, we have revised the sentence referring to genes marked by multiple altered histone modifications, so that it specifies which gene is being described rather than referring vaguely to "some of these genes".

      For the volcano plots, we will label the top-ranked DDRs, as we will also do for the volcano plots elsewhere in the manuscript.

      3-minor-results-7

      Figure 4B: only Grik1 as an example. Why only this one and not Bcl2 that moreover show several modifications? Could be helpful to show an example of each modification.

      Thank you for your comment. We will present the modification enrichment of specific genes that we mentioned.

      3-minor-results-8

      Section 4, Kdm6b: "To examine the possible contribution of Kdm6b to jSI, we re-analyzed the RNA-seq data from Kdm6b-knockout in the published study (Ramesh et al., 2023)." line 240-241: this study is about conditional Kdm6b KO in the cerebellum, on naive P14 male and female mice's neurons in culture. The authors extrapolate the results from a completely different neuronal population/region and sex to justify the potential effect jSI could have on their adolescent female mice. This sentence is misleading for the reader, since the model used (not stressed) and experimental conditions are far from what they are studying. This sentence needs some reformulation to better explain their goal. They show the DEGs (up/down) from reanalyzed data and overlap between these DEGs and the one from Figure1, but the conditions are far from each other here. One could ask what the specificity of their overlap demonstrated here.

      "Fosl2 and Nr4a1 are immediate early genes (IEGs) in response to neuronal activation in many brain regions (Dave et al., 2025; Shi et al., 2024), and these two genes have been reported to be involved in memory maintenance (McNulty et al., 2012; Mizuno et al., 2020) and Parkinson's disease (PD) (Fan et al., 2020; Rouillard et al., 2018). In addition, Htr5a, which encodes serotonin receptor 5A, was also upregulated in jSI and downregulated by Kdm6b KO (Fig. 1G, 5G, Table S1). And Htr5a has been reported to be a risk factor of human schizophrenia (Guan et al., 2016)." line 254-260: This part of the result paragraph is about introduction/discussion again. This should be moved appropriately.

      We thank the reviewer for these comments.

      Regarding the Ramesh et al. (Elife, 2023) dataset, we agree that the experimental conditions differ substantially from ours, and that our original wording did not make this clear. We have added a statement at the beginning of this section explaining how the datasets were selected and noting that they were obtained under conditions different from ours, so that the reader understands from the outset that these comparisons were intended to narrow down candidate genes rather than to test whether these enzymes act in the NAc after jSI. As described in our response to the Reviewer #3's comment (#3-major-5), we will also state the specific differences in brain region, cell type, developmental stage, and sex in the Limitations section, and we will reanalyze the striatal dataset from Chen et al. (Sci. Adv., 2022), which is considerably closer to our samples.

      Regarding the passage describing Fosl2, Nr4a1, and Htr5a, we agree that the discussion of their roles in memory and disease belongs in the Discussion rather than the Results. We have removed this material from the Results, retaining only the minimal information needed to follow why these genes were of interest, and we will incorporate the remainder into the Discussion. We have applied the same principle to the corresponding passage in the transcription factor section, as described in our response to the Reviewer #3's comment (#3-minor-results-3).

      3-minor-results-10

      Section 4, Setd1a: Here, they used 3 separate datasets: whole PFC of Setd1a heterozygous mice (exon 4 LacZ/Neo cassette insertion), whole PFC from loss of function Setd1a heterozygous mice and FoxP2+ nuclei from PFC of Setd1a +/- mice (frameshift in the 15th exon). These datasets are quite different between themselves and compared to NAc samples from jSI mice. The authors stated that the first two datasets had a low DEG overlap with their samples but continued with the third which showed a significant overlap for Hcrtr2, Dkk3 and Dact1.

      They finally conclude that: "Taken together, these results suggest that epigenetic factors, such as Kdm6b, Brd4, and Setd1a, may mediate jSI-induced gene expression alterations." Nothing in these datasets is comparable to what they want to prove here, it is a huge stretch to propose these genes as mediators of jSI. Reformulate. These results could be exploited as exploratory, to reduce the number of potential targets, but needs to be investigated on their own.

      We thank the reviewer for these comments, with which we largely agree.

      Regarding the differences between the three Setd1a datasets and our own samples, we have stated these in the Limitations section (lines 514–519), as described in our response to the Reviewer #3's comment (#3-major-5). We will also reanalyze the striatal dataset from the same study by Chen et al. (Sci. Adv., 2022), which is considerably closer to our NAc samples than the prefrontal cortex datasets, and we will reorganize this section accordingly.

      Regarding the difference in overlap between the three datasets, we would note that all three comparisons were performed and reported, and that the low overlap with the Mukai and Nagahama datasets was described in the original manuscript rather than omitted. One possible explanation for this difference is that the Chen dataset was generated from sorted Foxp2-positive nuclei, whereas the other two were derived from whole prefrontal cortex, in which signals from neurons may be diluted by non-neuronal cell types. We have added this to the text as a possible interpretation rather than a demonstrated explanation, and we have stated explicitly that all three datasets were compared in the same way (lines 318–322).

      Finally, we agree that proposing these enzymes as mediators of jSI-induced gene expression changes overstates what our data support. We have removed such wording and reframed these results as exploratory analyses that narrow down candidate genes for future investigation (lines 328–332), as described in our response to the Reviewer #2's comment (#2-major-1).

      3-minor-results-11

      Figure 5: volcano plots: top10 genes visible could be useful. This section of the results would fit better displayed in the supplementary, since they reanalyzed datasets far from their experimental conditions. These genes of interest should be further investigated in their jSI model.

      As described in our response to the Reviewer #3's comment (#3-major-7), we will reorganize this section in light of the reanalysis of the striatal dataset, retaining the most relevant comparison in the main text and moving the remaining reanalyses to the supplementary material. We have also revised the text so that each of these sections concludes by identifying the genes concerned as candidates requiring further investigation in our jSI model, rather than as established targets. We will label the top differentially expressed genes on the volcano plots, as described in our response to the Reviewer #3's comment (#3-minor-results-2).

      3-minor-discussion-7

      "Besides these two main shared functions affected by jSI, our results suggest that other biological processes are potentially mediated by one or more histone modifications. For example, some DDRs of H3K27ac and H3K27me3 are functionally enriched around cell adhesion-associated genes, and this is consistent with previous papers suggesting that cell adhesion is affected by isolation (Santiago et al., 2023; Wu et al., 2022)..." line 377-382: This GO term appeared in the figure, but has never been mentioned clearly in the results. The explanation goes on for a full paragraph. It could be better to introduce it before if it is of interest. Also, which cell-adhesion genes have been found in the RNAseq / cut&tag experiments for this family of genes (never stated)?

      We thank the reviewer for this comment. We agree that discussing cell adhesion at length in the Discussion is inappropriate when the corresponding GO term was never described in the Results, and that the genes underlying this enrichment were not identified.

      We will introduce this GO term in the Results section, where the functional enrichment of the DDRs is described, and we will name the cell adhesion-associated genes identified in our RNA-seq and CUT&Tag analyses both there and in the Discussion. This will be done together with the more comprehensive description of the top-ranked GO terms that we will add in response to the Reviewer #3's comments (#3-minor-results-1 and #3-minor-results-5).

      3. Description of the revisions that have already been incorporated in the transferred manuscript

      1-1

      The text in the introduction conflates adult and adolescent social isolation which have very different effects on behavior. Additionally, there is evidence that isolation during adolescence can have permanent effects on behavior but the behavioral effects of adult isolation in rodents are transient. It is recommended that the authors restructure the intro to be more specific to describing the adolescent period and why epigenetic mechanisms would be expected to regulate changes induced by jSI.

      We thank the reviewer for pointing out that adult and adolescent social isolation are conflated in the Introduction. We agree that the effects of social isolation differ between adulthood and adolescence, and we have revised the Introduction to clearly distinguish between the two, with a specific focus on the adolescent period. Specifically, we now note that isolation during adolescence can produce lasting behavioral alterations, whereas the effects of adult isolation are largely transient, and we have added a statement explaining why the adolescent period may therefore be particularly susceptible to epigenetic reprogramming (lines 69–74). In addition, we now explain why epigenetic regulation is a plausible candidate mechanism for the changes induced by jSI: since the effects of jSI persist long after the isolation period has ended, environmental stress during this window is likely to leave a lasting molecular trace within affected cells, and epigenetic regulation can stably maintain altered transcriptional states (lines 93–96).

      1-6

      Can the authors please explain why only females were used for these experiments? In addition, can the authors please comment on potential caveats in the interpretation by only including females in the study?

      We thank the reviewer for raising this point, which was also noted by the other reviewers. We apologize that this rationale was not stated explicitly in the original manuscript, and that our previous work was not cited in this context.

      Our focus on female mice was not arbitrary but followed from our previous work using the same isolation paradigm as in the present study. In Sazhina et al. (Neuroimage, 2025), in which mice were isolated from P21 to P35 and regrouped from P35 to P49, we found that jSI produced a heightened fear response in female but not male mice. Since the NAc has been implicated in scaling fear responses to threat intensity, and since this region undergoes a critical period around P28, we hypothesized that isolation during this window disrupts NAc development in a manner that leads to inappropriate fear responses in adulthood, and that this underlies the female-specific phenotype we had observed. We therefore designed the present study to examine molecular changes in the NAc of female mice. We have stated this rationale explicitly in the Introduction and Methods of the revised manuscript (lines 73–74, 542-543), and added the relevant references.

      We also agree that restricting the study to females limits the interpretation of our findings. Because jSI is known to produce sex-dependent effects on both behavior and gene expression, the alterations reported here cannot be assumed to occur in males, and comparisons with published datasets derived from male or mixed-sex animals must be made with this in mind. We have discussed these caveats explicitly in the Limitations section (lines 524–527), noting that a parallel analysis in males would be required to distinguish shared mechanisms of jSI from sex-specific ones.

      1-7

      Were females shipped to the facility on P21? It is unclear.

      We apologize that this was not clearly described in the original manuscript. Female mice were shipped from the breeder and arrived at our animal facility at P21, at which point isolation was started directly. We have stated this explicitly in the revised Methods section (lines 540–541). Group-housed control animals were shipped and received on the same day and under the same conditions, so that both groups experienced identical transport.

      1-8

      Were any animals used for multiple endpoints or was each endpoint a separate cohort? Were any samples pooled?

      We apologize for not describing this clearly. Nuclei were isolated from the NAc of a single animal and divided into five aliquots, each of which was used as one sample for RNA-seq or for one of the four CUT&Tag experiments. Each biological replicate therefore corresponds to a single animal, and no samples were pooled; all endpoints were derived from the same set of animals rather than from separate cohorts. Both our RNA-seq and CUT&Tag protocols have been optimized for use with small numbers of nuclei, which allowed all five libraries to be prepared from a single animal. We have stated this explicitly in the revised Methods section (lines 554–563).

      2-major-1

      Causality is not shown. The Discussion states: "Although we didn't show the molecular mechanism of histone modification alteration regulating gene expression, we revealed the association between transcriptome and histone modifications." That limit should shape the Abstract and title more clearly. Phrases like "epigenetic alterations may also play a role" are fine. Stronger wording about mediation should be toned down until NAc specific perturbation is done.

      We agree with the reviewer that this study is hypothesis-generating and does not demonstrate causality. We were mindful of this in preparing the original manuscript, but we acknowledge that language implying mediation remained in several places. We have gone through the entire manuscript, including the Abstract, and revised the wording so that it accurately reflects the correlative nature of our findings. In particular, we have removed expressions implying that the identified epigenetic factors mediate jSI-induced transcriptional changes, and now describe these relationships as associations that remain to be tested by NAc-specific perturbation (lines 34-35, 40-42, 197-200, 330-332, 433-434, 465-467, 481-483, 493-496).

      Regarding the title, we would prefer to retain the current wording. The title states that jSI is accompanied by alterations in gene expression and in histone modifications, and does not assert that the latter mediates the former; we therefore believe it does not overstate our findings. We note that the title has been modified to specify the sex of the animals used, as requested by Reviewer #3.

      2-minor-1

      Only female mice were used. State this early and discuss sex limits. Juvenile isolation effects often differ by sex.

      We thank the reviewer for this comment, and we agree on both points.

      As described in our response to the Reviewer #1's comment (#1-6), our focus on female mice followed from our previous work using the same isolation paradigm, in which jSI produced a heightened fear response in female but not male mice. We have stated this rationale explicitly in the Introduction and at the beginning of the Methods (lines 524-527, 542-543), and we have made clear in the Abstract and the title that this study was performed in female mice (line 42).

      We also agree that the sex-specific limitations of our findings require explicit discussion. Since jSI is known to produce sex-dependent effects on both behavior and gene expression, our results cannot be assumed to generalize to males, and comparisons with published datasets derived from male or mixed-sex animals must be interpreted with this in mind. We have addressed these points in the Limitations section of the revised manuscript.

      2-minor-3

      Figure 1 lists "II2ra" in the top nodes table. That is likely Il2ra. Please correct.

      We thank the reviewer for catching this error. The reviewer is correct that the gene name in the top nodes table should read Il2ra rather than "II2ra", and we have corrected this in Figure 1.

      2-minor-4

      Sample sizes differ a lot across marks. H3K4me1 and H3K27me3 have n = 4 in jSI. Discuss power and why replicates differ.

      We thank the reviewer for this comment. In this study, 10 control and 9 jSI animals were prepared, and all analyses were performed on nuclei derived from each of these animals. However, a subset of CUT&Tag libraries failed to pass our post-sequencing quality criteria and was excluded from the analysis, which resulted in the differing numbers of replicates across histone modifications. We have described this in the revised Methods (line 554), together with the quality criteria used for exclusion, so that the basis for these differences is transparent.

      We also agree that the reduced number of replicates for H3K4me1 and H3K27me3 lowers the statistical power for these marks relative to the others, and that the number of differentially distributed regions detected for them may therefore be underestimated. We have stated this explicitly in the Limitations section.

      2-minor-5

      The isolation protocol includes regrouping from P35 to P49. Make clear that effects are lasting post isolation effects, not acute isolation effects.

      We thank the reviewer for this comment, which correctly identifies our intent. The regrouping period was included precisely because our interest is in the effects of juvenile isolation that persist into adulthood, rather than in the acute consequences of isolation itself. Our previous work using the same protocol demonstrated a lasting fear phenotype in female mice after the regrouping period (Sazhina et al., Neuroimage, 2025), and a central aim of the present study is to ask whether epigenetic regulation contributes to the persistence of such environmentally induced changes.

      We have stated this rationale explicitly in the Methods (lines 548–550), and the Introduction now notes that the behavioral effects of jSI persist long after the isolation period has ended (lines 69–72, 93–94).

      2-minor-6

      Methods say "GPT-5.4... and Claude Sonnet 4.6... was used." Fix subject verb agreement.

      We thank the reviewer for pointing this out. We have corrected the subject-verb agreement in this sentence of the Methods (lines 682–683), which now reads "GPT-5.4 ... and Claude Sonnet 4.6 and Opus 5 ... were used".

      2-minor-7

      Data Availability lists "GSE3508789." Confirm this accession. It looks malformed.

      We thank the reviewer for catching this. The accession number was indeed malformed; the correct accession is GSE123652, and we have corrected it in the Data Availability section (lines 672, 717).

      2-minor-8

      Abstract keywords include "Loneliness." The mouse work is social isolation. Keep that distinction clear, as the Introduction already does.

      We thank the reviewer for this comment. We agree that including "loneliness" as a keyword was inappropriate given that this study examines social isolation in mice, and we have removed it from the keyword list (lines 44–45).

      3-major-1

      Title: Add the sex: "in female mice" since it is specific.

      We agree with the reviewer and have revised the title to specify that this study was performed in female mice (lines 1–3).

      3-major-2

      Referencing: Biorender.com has been used to generate some schematics in this study, but it is never acknowledged or referenced.

      We thank the reviewer for pointing out this omission. The schematic in Figure 1 was created using BioRender, and we have added the citation to the figure legend in the format specified by BioRender, together with the corresponding publication license (line 737).

      3-major-6

      Wording: Formulation throughout the current study is vague, sometimes misleading, with some unclear sentences (see minor comments for the sentences showing problems). This issue needs to be checked again.

      Thank you for your comment. The responses could be checked in minor comments part.

      3-minor-intro-1

      "These negative effects are further supported by evidence from Covid-19 during the last few years" line 56/57: reformulate.

      We thank the reviewer for this comment. We agree that the original sentence was awkwardly phrased, since it referred to evidence "from Covid-19" rather than to the studies conducted during that period, and since "the last few years" was vague. We have rewritten it to state that studies conducted during the COVID-19 pandemic, when social contact was widely restricted, provided further evidence for these negative effects (lines 56–58).

      3-minor-intro-2

      "In addition, isolation contributes to severe social issues, such as increased human suicide risk" line 57/58: issues is plural, but only one example is given; moreover, the term "social issue" associated to suicide is poorly-worded.

      We thank the reviewer for this comment. We agree that the plural "issues" was not supported by the single example given, and that describing suicide as a "social issue" was poorly worded. We have rewritten the sentence so that social isolation is described as being associated with adverse outcomes, including the risk of suicide (lines 59–60).

      3-minor-intro-3

      "In the case of rodents, socially isolated animal models are proposed to be associated with various human diseases" line 59/60: clumsy sentence, the animal models are associated to human disease? This could be reformulated.

      We thank the reviewer for this comment. We agree that the original sentence was awkwardly constructed, since it stated that the animal models themselves were associated with human diseases. We have rewritten it so that the socially isolated animals are described as exhibiting a range of behavioral abnormalities, and these abnormalities, rather than the models themselves, are described as modelling aspects of human diseases such as depression and schizophrenia (lines 60–63).

      3-minor-intro-4

      "and the effects of juvenile social isolation (jSI) on motor, emotional, learning, and sociability-related behaviors in rodents have been widely reported (Li et al., 2021; Powell & Swerdlow, 2023; Walker et al., 2019), which further provides evidence of the pathogenesis and molecular mechanisms of human mental disorders." line 63-66: Examples of behavioral dysfunctions would be appreciated here.

      We thank the reviewer for this suggestion. We have added examples of the behavioral dysfunctions reported after jSI, namely hyperactivity, elevated anxiety, impaired spatial learning, and altered social play (lines 66–68). In the same paragraph we have also added a description of the sex differences reported for these effects, and a reference to our own previous work using the same isolation paradigm, as described in our responses to the Reviewer #1's comments (#1-2 and #1-6).

      3-minor-intro-5

      "The nucleus accumbens (NAc) is a critical component of the brain reward circuitry, and dysfunction of the NAc is associated with drug addiction (Zinsmaier et al., 2022), impaired social interaction (Pomrenze et al., 2022; Shan et al., 2022), and abnormal emotion expression (Gebara et al., 2021)" line 67-70: here, the statement reads as dysfunction of the NAc is responsible of abnormal behaviors (addiction, social behavior or emotional expression), but the articles show that NAc is dysregulated in models of these pathologies. Is the dysfunction in the NAc responsible of or a consequence of the pathologies? This could be better formulated.

      We thank the reviewer for this comment. We agree that the original wording could be read as asserting that NAc dysfunction causes these conditions, whereas the cited studies show that the NAc is dysregulated in models of them. We have rewritten the sentence so that NAc dysfunction is described as having been reported in animal models of these conditions, without implying a direction of causation (lines 81–85).

      3-minor-intro-6

      "An fMRI study showed that activity of the human NAc is associated with the sense of loss (Cooper et al., 2009; O'Connor et al., 2008)." line 70-72: Sentence says one study, but two references are used. The sentence refers to O'Connor only. Cooper is about reward/effort and NAc activity, not grief/loss, reformulate. Also, "activity" is unclear, the authors could be more precise with "hyperactivity of the NAc has been found in people suffering from loss".

      We thank the reviewer for pointing out these problems. We have removed the citation to Cooper et al. (2009), which concerns reward and effort rather than grief, so that the sentence now refers only to O’Connor et al. (Neuroimage, 2008) (lines 85–86). We have also replaced the vague reference to "activity" with a statement that hyperactivity of the NAc was reported in individuals experiencing loss, as the reviewer suggested.

      3-minor-intro-7

      "The NAc from lonely individuals showed key differentially expressed genes (DEGs) that are associated with both neurodegenerative and neuropsychiatric diseases" line 74-75: Unclear, give examples of the pathologies here to be consistent with the next sentence about female rats (Alzheimer, Parkinson, Huntington).

      We thank the reviewer for this suggestion. We have added the specific diseases identified in that study, namely Alzheimer's disease, Parkinson's disease, and major depression disorder, so that this sentence is consistent with the following sentence describing the findings in female rats (lines 88–92).

      3-minor-intro-8

      The next paragraph (line 79-95) about DNA methylation and histone modifications is missing a general conclusion: what is interesting or needs to be more studied? Also, H3K9 and H3K79 have been introduced but unused in the paper, while H3K27ac/me3 have not been introduced. What is known about them?

      We thank the reviewer for these comments. We agree that this paragraph lacked a conclusion and that the histone modifications discussed did not match those examined in this study.

      We have shortened the description of H3K79 methylation and removed the statement concerning H3K9, neither of which is examined here. In their place we have added a description of the four modifications we analyse, namely H3K4me1, H3K4me3, H3K27ac, and H3K27me3, together with what is known about their roles in the NAc (lines 108–118). The paragraph now concludes by noting that little is known about H3K4me3 and H3K27ac in the NAc under stress, and that how any of these modifications are altered after jSI remains unknown, which motivates the present study.

      3-minor-intro-9

      "How do epigenetic elements mediate gene expression dysfunction under jSI stress? In this study, we aimed to reveal the alterations in gene expression and histone modifications induced by jSI, and to elucidate their roles in the context of psychiatric disorders promoted by jSI" line 96-99: Statement is too general, it is missing the term "NAc" here.

      We thank the reviewer for this comment. We agree that our statement of aims was too general and omitted the brain region under study. We have rewritten both the question that opens this paragraph and the statement of aims so that they specify the NAc, and we have also indicated that the study was performed in female mice (lines 119–122). In addition, we have removed the wording implying that epigenetic elements mediate gene expression dysfunction, in line with the request from Reviewer #2 that causal language be moderated (#2-major-1).

      3-minor-methods-1

      Female mice only have been used in this study. It is never explained why so. Knowing that many diseases show sex differences in prevalence or symptom expression, it would have been helpful to include also male mice in the study, to identify common mechanism linked to jSI versus sex-specific alterations.

      We thank the reviewer for this comment, and we apologize that our rationale was not stated in the original manuscript. As described in our response to the Reviewer #1's comment (#1-6), our focus on female mice followed from our previous finding that the same isolation paradigm produced a behavioral phenotype in females but not males (Sazhina et al., Neuroimage, 2025). We have added a statement to this effect at the beginning of the Animals and sample collection section (lines 542–543), and we have discussed the resulting limitations in the Limitations section (lines 524–527).

      3-minor-methods-2

      Cut & Tag: Dilution of antibodies is not displayed ("1µL") and the reference for antibodies is unclear: "primary antibodies (H3K4me1, MABI, 536 MABI0302; H3K4me3, abcam, ab8580; H3K27me3, CST, 9733S; H3K27ac, CST, 537 8173S; 1 µL per reaction)". This should be adapted to look like the FACS antibody description: "anti-NeuN-488 conjugated antibody (Millipore, #MAB377X, 1:400 dilution)".

      We thank the reviewer for pointing this out. We have revised the description of the CUT&Tag antibodies so that it follows the same format as the FACS antibody description, specifying the host species and clonality, the supplier, the catalogue number, and the dilution for each antibody, in place of the volume per reaction given previously (lines 611–615).

      3-minor-methods-3

      "Frozen nuclei were thawed and bound to concanavalin A (ConA)-coated magnetic beads (BioMag®Plus Concanavalin A, 10 µL per reaction) for 10-60 min on a rotator at room temperature." Why so much difference in incubation time here?

      We thank the reviewer for pointing this out. We have checked our experimental records and confirmed that the incubation was performed for 10–20 min in all experiments reported here, and we have corrected the text accordingly (line 608). The wider range given in the original manuscript reflected our general protocol, in which we have confirmed that binding is satisfactory anywhere between 10 and 60 min, but this was not the range actually used in the present study.

      3-minor-methods-4

      Number of animals: While reading the manuscript, it was unclear that 5 different groups of mice have been used. The number of animals should be stated in the methods in the "nucleus extraction" part or "animals" section to facilitate understanding the methods.

      We thank the reviewer for this comment. We have renamed this section "Animals and sample collection" and added a paragraph stating the number of animals analysed in each group, that nuclei from each animal were divided between the RNA-seq and the four CUT&Tag experiments, and that no samples were pooled (lines 539, 554–563). We have also explained that the number of replicates is smaller than the number of animals for some datasets because a subset of libraries did not pass our quality criteria, and we note that the replicate number for each dataset is given in the corresponding figure legend.

      3-minor-methods-5

      Data analysis: "For RNA-seq data, p-value 1.2 were used as the threshold for identifying DEGs. For CUT&Tag data, p-value 2 were selected as the standard for identifying DDRs." Cut&Tag p-value and threshold is written in RNAseq section, move it to its proper part hereafter "CUT&Tag data analysis".

      Thank you for your comment. We revised it (lines 664–665).

      3-minor-methods-6

      Suggestion DDR: acronym is present in the methods but not explained. It is however explained in the results. This depends on the order in the publication, but if methods appear first, it would be helpful to understand what stands for DDR.

      Thank you for your comment. We revised it (line 665).

      3-minor-methods-7

      Cut&Tag data analysis: "The procedures of quality check and trimming were the same as RNA-seq data analysis." line 591; "and the removal of blacklisted regions was the same as RNA-seq" line 594; "GO analysis was the same as RNA-seq analysis." line 599: These parts could be ameliorated to avoid repetition. Since the preparation of nuclei and most of the analysis are the same, the methods could be more straightforwardly explained separating common preparation from specific analysis.

      Thank you for your comment. We revised it (lines 659–660).

      3-minor-methods-8

      Methods explaining how the authors performed the reanalysis of the public datasets is missing.

      We thank the reviewer for pointing out this omission. We have added a "Public data analysis" section to the Methods, describing how the raw data were retrieved from the DDBJ and GEO databases, and how they were processed (lines 667–671). Steps shared with our own datasets are indicated as such rather than repeated in full.

      3-minor-methods-9

      Animals have been separated at P21 and regrouped at P35: It is not stated if they were regrouped together or with a group of unstressed WT never separated, which could influence their behavior and stress levels.

      We thank the reviewer for this comment. We have clarified that the isolated mice were regrouped with other previously isolated animals, rather than with mice that had never been separated (lines 547–550). As described in our response to the Reviewer #2's comment (#2-minor-5), we have also stated why the regrouping period was included.

      3-minor-methods-10

      No behavioral test has been performed on these mice. It would have been appreciated to see that 2-weeks social isolation was efficient to generate stress in these animals (anxiety test, sociability at least). And if this protocol has been previously used in their lab, at least to explain briefly what behavior abnormalities the jSI was inducing.

      We thank the reviewer for this comment, and we apologize that this information was not included in the original manuscript.

      This isolation protocol has been used previously in our laboratory, and the resulting behavioral phenotype was reported in Sazhina et al. (Neuroimage, 2025), in which the same paradigm produced a heightened fear response in female mice. We have added this finding to the Introduction (lines 73–74), and we have also cited it in the Methods as the basis for our use of female animals (lines 542–543), so that the behavioral consequences of the paradigm are documented.

      As described in our response to the Reviewer #2's comment (#2-major-4), behavioral testing was not performed on the cohorts used for molecular analysis, in order to avoid introducing transcriptional and epigenetic changes unrelated to isolation.

      3-minor-results-3

      Section 2: The first paragraph here is about which Transcription Factor (TF) is predicted to participate in their DEG's expression. Half of this paragraph is introduction about the function of several TF. This is not part of results and should be moved appropriately in the introduction or discussion section, or shortened significantly, since it is now longer than the result part. Importantly here, the analysis is done on ChIP Atlas (public datasets).

      They state: "Taken together, promoter analysis of DEGs suggests that potential epigenetic mechanisms may act upstream of jSI-induced transcriptional dysregulation in the NAc." line 176, but the database has never been stated to be NAc-only data nor data from jSI animals. If not, this sentence has to be modified. It was unclear globally if this part was based on their work or data mining on a first read, it should be more clearly stated at the beginning that this is exploratory.

      We thank the reviewer for these comments.

      We agree that the introductory description of the transcription factors was disproportionately long for a Results section. We have shortened it substantially, retaining only the information required to follow why these factors were of interest, namely that Setd1a and Brd4 act through the histone modifications examined in this study (lines 184–190). The remaining background, including the association of these factors with neurological and psychiatric disease, has been moved to the Discussion.

      We also agree that the exploratory nature of this analysis was not stated clearly. The ChIP-Atlas database is compiled from published ChIP-seq experiments and does not contain data from the NAc of socially isolated animals. We have stated this at the outset of the section (lines 177–180), so that the reader understands from the beginning that the analysis was intended to generate candidates rather than to identify regulatory relationships operating in our system, and we have revised the concluding sentence so that it no longer implies that these mechanisms were demonstrated in the NAc under jSI (lines 199–200).

      3-minor-results-9

      Section 4, Brd4: Here, the authors reanalyzed data from E16.5 cortical neuronal culture treated with or without BET family inhibitor, which they state is not selective of Brd4 (even if it is part of the BET family). The crossover between this embryonic cortical neuronal population treated with nonspecific inhibitor and their model (juvenile Social Isolation, NAc) is a bit of a stretch. What do the overlap in DEGs really mean here?

      "We also examined the possible downstream Brd4 target genes within the gene sets of down-DEGs by jSI and down-DEGs by JQ1 treatment, and we identified Hcrtr2 and Dkk3 in these gene sets (Fig. 1G, 5H, Table S1). Hcrtr2 encodes an orexin receptor, and it has been reported to be involved in altered arousal levels through dopamine neurons (Bandarabadi et al., 2024). Dkk3 inhibits Wnt signaling and is reported to be related to anxiety and memory formation (X. Chen et al., 2025; Flores et al., 2024)." line 275-280: What is the conclusion on these results?

      We thank the reviewer for these comments.

      Regarding the Korb et al. (Nat. Neurosci., 2015) dataset, we have added a statement at the beginning of this section explaining how the datasets were selected and noting that they were obtained under conditions different from ours, so that the reader understands from the outset that these comparisons were intended to narrow down candidate genes rather than to test whether these enzymes act in the NAc after jSI (lines 261–265). As described in our response to the Reviewer #3's comment (#3-major-5), the specific differences are also stated in the Limitations section (lines 514–519).

      Regarding the meaning of the overlaps, we agree that our original wording did not convey this clearly, and in particular that opening with "as expected" was misleading given that we also observed a significant overlap in the opposite direction. We have rewritten this passage so that the overlap in the unexpected direction is stated as an independent observation rather than as a qualifying clause, and we now state explicitly that the two gene sets are related but that the direction of change does not correspond in a simple manner (lines 295–301).

      Regarding the conclusion of this section, we agree that none was previously given. We have removed the description of the roles of Hcrtr2 and Dkk3 in arousal, anxiety, and memory, which belongs in the Discussion, and have instead concluded the section by identifying these two genes as candidates whose expression may be regulated by Brd4 in the context of jSI (lines 304–306). We have applied the same principle to the corresponding passages in the Kdm6b and Setd1a sections, as described in our response to the Reviewer #3's comment (#3-minor-results-8).

      3-minor-discussion-1

      " For example, the expression of glutamate receptors is reduced in the NAc, prefrontal cortex, and hippocampus under isolation stress (Hermes et al., 2011; Mao et al., 2022; Sestito et al., 2011)." line 315-317: Which GluR are reduced here? It needs to be more precise for the reader here, and to state if some genes have been found in common between this literature and their DEGs.

      We thank the reviewer for this comment. We agree that the original sentence was too vague, and we have revised it to specify which glutamate receptors were reduced in the cited studies (lines 343–351), namely GluA1 and GluA3 in the NAc and caudate putamen under chronic social isolation stress.

      We now also state explicitly how these findings relate to our own data. Gria1 and Gria3, encoding GluA1 and GluA3, respectively, were not among our DEGs, but we identified other glutamatergic synapse-associated genes, including Grin3a and Grik1, and we note that isolation may therefore affect glutamatergic signaling in the NAc through a partly distinct set of genes under our conditions.

      3-minor-discussion-2

      "The NAc is a key component of the brain reward circuit, and it is involved in drug addiction and social behavior (Pomrenze et al., 2022; Zinsmaier et al., 2022). NAc neurons receive glutamatergic inputs from the PFC, basolateral amygdala (BLA), hippocampus, and ventral tegmental area (VTA) (Arrondeau et al., 2024; Dieterich et al., 2021; Elam et al., 2025; Le Borgne et al., 2025; Zinsmaier et al., 2022), and neurons in the NAc output the information to the ventral pallidum (VP) (Liu et al., 2022), VTA (Qi et al., 2022), and other areas of the basal ganglia (Lanciego et al., 2012)." line 318-324: These lines are describing the circuitry of the NAc, some of its inputs (no mention of dopamine afferences from the VTA) and outputs. No use of this information is used after, since they conclude the paragraph with: "Thus, deficits in glutamatergic synapses possibly mediate jSI-induced behavioral abnormalities, including impaired social interaction, anxiety, and an increased risk of substance abuse." line 325-326: What is the point of describing the circuit, if it is not interpreted regarding their results? What is their hypothesis on the circuit dysfunction in jSI female mice? They were discussing the DEGs from RNAseq result before this paragraph. What is the link/hypothesis between their DEGs and the glutamatergic circuits of the NAc? Is the NAc directly responsible of jSI-induced behavioral abnormalities for them or cortical/amygdal/hippocampal/VTA glutamatergic projection neurons are dysregulated, creating DEGs at the synapse in the NAc? This part of the discussion should be more specific on what they mean.

      We thank the reviewer for this comment. We agree that the description of the NAc circuitry was not connected to our own findings, and we have substantially shortened it, retaining only a single sentence summarizing the glutamatergic inputs to the NAc and its outputs to downstream regions (lines 351–356).

      We have also revised the concluding sentence of this paragraph so that it follows from the preceding discussion of our DEGs (lines 351–356). Since our data indicate that glutamatergic synapse-associated genes are downregulated in NAc neurons after jSI, and since the NAc integrates glutamatergic inputs from several regions implicated in social and emotional behavior, we now state that altered glutamatergic signaling at these synapses may contribute to the behavioral abnormalities induced by jSI, rather than asserting that such deficits mediate them.

      3-minor-discussion-3

      "Deficiencies in these proteins are associated with various behavioral abnormalities (Araujo et al., 2017; Chasse et al., 2024; Guo et al., 2020; Huang et al., 2021; Mukai et al., 2019)." line 334-335: what proteins and what behavioral abnormalities? This is not precise enough and needs reformulation/conclusions.

      We thank the reviewer for this comment. We agree that the original sentence was not sufficiently specific, since it referred to deficiencies in several proteins and to behavioral abnormalities without indicating which protein was associated with which phenotype.

      We have revised this passage to describe the reported phenotypes individually for each factor (lines 362–370). We have also incorporated here the background material on Setd1a and Brd4 that we removed from the Results section, as described in our response to the Reviewer #3's comment (#3-minor-results-3), so that the association of these factors with neurological and psychiatric conditions is presented in the Discussion rather than interrupting the presentation of our findings.

      3-minor-discussion-4

      "A previous report suggests that histone modifications such as H3K4me3 in the hippocampus respond to an enriched environment (Schaffner et al., 2023), and our results indicate that these histone modifications may influence gene expression in the NAc under jSI stress as well." line 342-344: In which direction is the modification in the hippocampus in enriched environment? Is it opposite to what the authors have found in jSI (which would be interesting, since one could see a more social environment as an enriched condition too)? The idea behind this sentence needs to be precised.

      We thank the reviewer for this comment. We have revised this sentence to describe the reported finding more precisely (lines 376–379), namely that the loss of H3K4me1 observed in SNCA transgenic mice was partially dampened by environmental enrichment.

      Regarding the reviewer's suggestion that enrichment might be viewed as the converse of isolation, we agree this is an interesting possibility, but we have chosen not to develop the comparison in the text. Social isolation and environmental enrichment are not straightforwardly opposite conditions, since the absence of social contact is not simply the inverse of enrichment relative to standard housing, and the two may engage distinct circuits and cell populations. We therefore refer to this study only as evidence that these histone modifications are responsive to the housing environment, rather than drawing a directional comparison with our own data.

      3-minor-discussion-6

      "Since neural development relies on the regulation of gene expression (Jain et al., 2001; Xiang et al., 2020), we hypothesize that these terms reflect altered gene expression regulation mechanisms under jSI stress." "However, these hypotheses need to be further validated by additional experiments." line 369-371 & 375-377: The authors are not integrating their results, they are being cautious, but the message stays unclear to the reader. What is the message here?

      We thank the reviewer for this comment. We agree that our original wording was cautious to the point of leaving the message unclear, and we have rewritten this passage.

      We now state explicitly what we wish to propose: that genes involved in transcriptional and chromatin regulation carried altered histone modifications, that changes at such loci may have consequences extending beyond the genes themselves through their downstream targets, and that this may be particularly relevant during the developmental window examined here (lines 404–416). The paragraph now concludes by presenting this as a hypothesis, namely that histone modification changes at regulatory genes act as an upstream event after jSI whose consequences are amplified through the targets of those regulators, together with a statement that this remains to be tested experimentally, rather than ending with a general remark that further validation is required.

      3-minor-discussion-8

      "To determine whether histone modifications regulate specific genes, we focused on potentially important genes. Grik1, for example, exhibits reduced H3K27ac levels. It encodes a subunit of ionotropic glutamate receptors, and its deficiency has been found in mental diseases, such as schizophrenia and ADHD (Chatterjee et al., 2022; Hirata et al., 2012). The inactivation of Grik1 in rodents promotes anxiety-like behaviors via glutamatergic transmission (Englund et al., 2021)." line 387-392: What is the conclusion/hypothesis on Grik1's role?

      We thank the reviewer for this comment. We agree that the original passage described what is known about Grik1without stating what we ourselves wished to conclude.

      We have added a statement of our hypothesis at the end of this passage: that the reduction in H3K27ac around the Grik1 locus contributes to the downregulation of Grik1 after jSI, which in turn may contribute to the anxiety-like phenotypes associated with isolation (lines 432–437). We also indicate what would be required to test this, namely manipulating H3K27ac at the Grik1 locus and assessing the resulting transcriptional and behavioral changes.

      3-minor-discussion-9

      "Bcl2, for example, has downregulated H3K4me1, H3K4me3, and upregulated H3K27me3 levels. Bcl2 is an apoptosis-related gene that determines neuronal survival under stress. A previous study suggests that chronic social defeat stress decreases the Bcl-2/Bax ratio in NeuN+ neurons in the hippocampus (Zhu et al., 2024). Our data suggest that jSI is another type of stress that suppresses Bcl2 expression, and that the epigenetic factors are possible upstream regulatory mechanisms." line 393-399: No links or clear hypothesis have been made here. The authors proposed to go deeper in this direction later. It would be interesting to conclude on the hypothesis on these two genes (Grik1/Bcl2) in their model.

      We thank the reviewer for this comment. We agree that this passage describes our observations concerning Bcl2 without stating what we conclude from them.

      We have revised it to present our hypothesis explicitly, namely that the coordinated reduction of H3K4me1 and H3K4me3 together with the increase in H3K27me3 around the Bcl2 locus contributes to its downregulation after jSI, and that this may in turn affect neuronal survival under stress (lines 443–446). We have also indicated what would be required to test this, in the same way as for Grik1, as described in our response to the Reviewer #3's comment (#3-minor-discussion-8).

      3-minor-discussion-10

      "We found that Kdm6b may be at least partially involved in the gene expression changes in jSI mice, especially for potentially important genes like Nr4a1. Nr4a1 encodes a transcription factor that regulates dopamine metabolism, and it is reported to be involved in drug addiction, which is probably mediated by histone modification alterations." line 408-411: The authors are very cautious about their conclusion, maybe too much. Since they base their hypothesis on P14 cerebellum neurons in culture, more work is needed here to establish clearly the role of Kdm6b. It feels like the authors are dropping clues for the reader, without concluding themselves on their hypothesis.

      We thank the reviewer for this comment, which we found particularly helpful. We recognize that our repeated use of hedging language left the reader uncertain as to what we were actually proposing.

      We have restructured these passages so that our hypothesis is stated explicitly, followed by a statement of what would be required to test it (lines 455–469). For Kdm6b, we now propose that Kdm6b-mediated H3K27 demethylation contributes to the upregulation of Nr4a1 in the NAc after jSI, and that testing this hypothesis will require NAc-specific manipulation of Kdm6b together with assessment of both transcriptional and behavioral outcomes. The limitations of the published dataset from which this hypothesis derives, including its origin in cerebellar tissue at P14, are stated in the Limitations section, as described in our response to the Reviewer #3's comment (#3-major-5).

      We believe this approach addresses the reviewer's concern without overstating our findings, and it is consistent with the request from Reviewer #2 that causal language be moderated (#2-major-1).

      3-minor-discussion-11

      "Our experiments suggest that Dkk3 is regulated by jSI stress, and this is probably mediated by Brd4." / "Our analysis also implies that Dkk3 and Hcrtr2 are probably regulated by other epigenetic regulators such as Setd1a." line 429-430 & 431-432: "Probably". No proof is given on that statement. Everything is "potential" or "possible" or "probable" here.

      We agree with the reviewer. As described in our response to the Reviewer #3's comment (#3-minor-discussion-10), we have gone through the Discussion and removed or replaced the repeated use of "probably", "possible", and "potential", restructuring the relevant passages so that each hypothesis is stated explicitly, together with a statement of what would be required to test it (lines 476–483, 484-496).

      3-minor-discussion-12

      "To summarize, we revealed the changes in gene expression and histone modification levels under jSI stress." line 444: Here they are missing the term "in the NAc" and "in female mice" and maybe need to add "some changes" in the sentence.

      We thank the reviewer for this comment. We have revised the summary sentence to specify both the brain region and the sex of the animals studied, and to indicate that we identified a subset of changes rather than a comprehensive account (lines 498–499). It now reads: "To summarize, we revealed some changes in gene expression and histone modification levels in the NAc of female mice under jSI stress."

      3-minor-discussion-13

      Major issue of this study is stated as "Another limitation related to this is that we inferred candidate epigenetic factors based on previously published public data. However, different brain regions, cell types, and experimental conditions between our analyses and public datasets may contribute to the gene expression differences, leading to false negative or false positive results in this study." line 456-460: Indeed, as they highlight here, the differences between the datasets chosen and their initial context (jSI, NAc) is huge. Maybe the authors could explain better why did they choose these datasets instead of more similar ones.

      Also, in the RNAseq / cut&tag analysis, they "were unable to separate these subtypes (D1+ or D2+) for this current analysis". This analysis could be interesting to see in the supplementary or a brief description of what they have tried in the discussion, since the composition of the NAc is made of about 90-95% of MSNs and a plethora of interneurons (Parvalbumin, Cholinergic, etc), leading to different circuit connectivity and function. The striatum also contains a third population of D1/D2 hybrid neurons, maybe this population (about 5-10% of the whole striatum) made the identification of the neuronal subtypes harder.

      We thank the reviewer for these comments.

      Regarding the choice of public datasets, we have revised the Limitations section to state that these datasets were selected as the closest available to our system, since no dataset in which these factors had been perturbed exists for the NAc or striatum, and we have added developmental stage and sex to the list of differences between those datasets and our own (lines 506–527).

      Regarding the separation of neuronal subtypes, our nuclei were sorted using an anti-NeuN antibody, which labels neuronal nuclei broadly and does not distinguish D1- from D2-positive medium spiny neurons. We did not attempt to separate these subtypes, and we have now stated this reason explicitly in the Limitations section. As the reviewer notes, this is a meaningful limitation, since changes restricted to one subtype, or occurring in opposite directions between subtypes, would be underestimated or missed entirely in our analysis. We agree that subtype-resolved approaches will be required to address this in future work.

      4. Description of analyses that authors prefer not to carry out

      2-major-4

      Behavior is missing from this study. The Introduction says jSI affects "motor, emotional, learning, and sociability-related behaviors." This manuscript does not show that the molecular changes track those phenotypes in the same cohort. Without that link, the psychiatric disease framing stays speculative.

      We thank the reviewer for this comment, and we apologize that our reasoning was not made clear in the original manuscript.

      Behavioral phenotypes produced by this isolation paradigm have been characterized in our previous study (Sazhina et al., Neuroimage, 2025), which used exactly the same protocol as the present work, with isolation from P21 to P35 and regrouping from P35 to P49. We have described these findings in the Introduction so that the behavioral consequences of our paradigm are explicit (lines 554–557), rather than referring only to the literature in general terms.

      We did not perform behavioral testing on the animals used for molecular analysis, because behavioral testing itself constitutes a substantial stimulus. Exposure to a novel environment, handling, learning experience, and aversive stimuli such as the foot shock used in fear conditioning all induce transcriptional and epigenetic changes in the brain, including the induction of immediate early genes such as Nr4a1 and Fosl2, which are among the candidate genes discussed in this study. Had the same animals been subjected to behavioral testing, we would not have been able to attribute the observed changes to jSI rather than to the testing procedure. We therefore used dedicated cohorts for molecular profiling, and we have stated this rationale explicitly in the revised Methods.

      We acknowledge, nevertheless, that this design means we cannot demonstrate a correspondence between molecular changes and behavioral phenotypes within the same individuals. We have stated this limitation in the Limitation section and have moderated the framing of our findings in relation to psychiatric disease accordingly (lines 524–527), as described in our response to the Reviewer #2's comment (#2-major-1).

      3-major-3(OPTIONAL)

      Males: Experiments are only focused on females here. The authors didn't state why. These experiments (RNAseq, Cut&Tag) could be performed independently in male mice. The introduction refers to several pathologies that show sex-bias in human or animal models and even articles with sex differences. If the authors had a reason to select female only, it should be clearly stated. This proposition would take the same amount of resources and time that for this issue but would increase the knowledge about the (sex differences in the) effect of juvenile social isolation greatly.

      We thank the reviewer for this suggestion, and we agree that a parallel analysis in male mice would substantially extend the value of this work.

      As described in our response to the Reviewer #1's comment (#1-6), our focus on female mice was not arbitrary. In our previous study using the same isolation paradigm (Sazhina et al., Neuroimage, 2025), jSI produced a heightened fear response in female but not male mice, and the present study was designed to examine the molecular basis of that female-specific phenotype in the NAc. We apologize that this rationale was not stated in the original manuscript, and we have described it explicitly in the Introduction and Methods of the revised version (lines 73–74, 524-527, 542-543).

      We would nevertheless like to explain why we are unable to perform the proposed experiments within the scope of this revision. Generating a comparable dataset in males would require a new cohort of animals, isolation and regrouping over four weeks, nuclear isolation and sorting, and the preparation and sequencing of five libraries per animal, followed by the full analysis pipeline. As the reviewer notes, this would require resources and time comparable to those invested in the present study, and it is not feasible within the revision period. We have stated in the Limitations section that a parallel analysis in males is required to distinguish shared mechanisms of jSI from sex-specific ones (lines 524–527), and we intend to pursue this in future work.

      3-major-4(OPTIONAL)

      Behavior abnormalities: Showing behavior abnormalities after jSI of female mice or, if it has been published somewhere else previously, a general description of the phenotypes observed in juvenile and adult female mice. This experiment could be performed in one batch of female separated at P21 and regrouped at P35, with anxiety (openfield or elevated plus maze, 1 day each), sociability (direct or 3-Chambers, 1 day each) and eventually depressive-like behaviors/anhedonia (sucrose preference test, 1 week including habituation to the two bottles and/or tail suspension/forced swimming, 1 day) tested. Same tests would be performed in juveniles or in adults.

      We thank the reviewer for this suggestion, and for noting that a description of previously published phenotypes would be an acceptable alternative to new behavioral experiments.

      Behavioral phenotypes produced by this isolation paradigm have been characterized in our previous study (Sazhina et al., Neuroimage, 2025), which used exactly the same protocol as the present work, with isolation from P21 to P35 and regrouping from P35 to P49. In that study, jSI produced a heightened fear response in female but not male mice. We have described these findings explicitly in the Introduction (lines 73–74), so that the behavioral consequences of our paradigm are stated rather than left to the reader to infer from the general literature.

      As described in our response to the Reviewer #2's comment (#2-major-4), we deliberately did not perform behavioral testing on the animals used for molecular profiling, since behavioral testing itself induces transcriptional and epigenetic changes in the brain and would have confounded the changes attributable to isolation. We have stated this rationale explicitly in the revised Methods (lines 554–557).

      3-major-9(OPTIONAL)

      Include more mechanistic experiments (as suggested in the discussion) on some of the factors identified (Kdm6b, Brd4, Setd1a) to better confirm their involvement in the jSI-induced changes in transcriptome (and behavior). Conditional KO in the NAc with viral infection (mouse line Kdm6b-flox or Brd4-flox or Setd1a-flox existing + AAV-cre) or AAV Crispr for specific knockdown could hardly be performed in the time window used in this study (at least 3week expression of the Cre/Cas9). But systemic or intracerebral pharmacological approach (ip injection of an inhibitor or cannula implantation, local intra-accumbens injection) could be performed in juvenile (1 week recovery post-surgery only, can be done at P21 before isolation).

      We thank the reviewer for this thoughtful suggestion, and we agree that functional analysis of Kdm6b, Brd4, and Setd1a would substantially strengthen the conclusions of this study.

      As the reviewer notes, conditional knockout approaches requiring viral expression are not feasible within the developmental window used here. Regarding the pharmacological alternatives proposed, systemic administration would not allow us to attribute any resulting changes specifically to the NAc, and local intra-accumbens administration, while addressing this point, would require establishing a new surgical and behavioral pipeline in juvenile animals alongside the molecular analyses. Neither is achievable within the revision period.

      We have therefore stated explicitly in the Discussion that functional validation of these candidate factors is required to test the hypotheses raised here (lines 509–513), and we intend to pursue this in future work. We are grateful to the reviewer for these constructive suggestions, which we will take up in designing those experiments.

      3-minor-discussion-5

      "Since our mice were isolated from P21 to P35, a period critical for the maturation of neurons (Makinodan et al., 2012; Walker et al., 2019; Yamaguchi et al., 2024), it is possible that the neuronal development process is affected by the isolated housing environment." line 351-354: This sentence states that juvenile social isolation during development could impair development. It is probable, since early life stress (even maternal stress) can have prolonged effect on the offspring behaviors. It is probable that the time-window of jSI is important for development.

      OPTIONAL: The authors would benefit of more experiments here: 1) They could perform SI in the same conditions but in adult females and compare the DEGs observed in that case (less development related genes probably). 2) "Neurons during adolescence mainly experience synaptic pruning and elimination (Afroz et al., 2016; Germann et al., 2021; Watanabe & Kano, 2024), and isolation stress probably impairs such processes." Here, the author could check in the NAc of their jSI female mice the state of dendritic arborization of MSNs (number, length, etc) on NAc brain slices.

      We thank the reviewer for these suggestions, both of which we agree would strengthen the study.

      Regarding social isolation in adult females, a comparison with adult-isolated animals would indeed help to establish whether the changes we observed are specific to the adolescent period. This would require a new cohort, a full isolation and regrouping schedule, nuclear isolation and sorting, and library preparation and sequencing, and is not feasible within the revision period. We will note this as a direction for future work. We have, however, revised the Introduction to cite studies defining P21 to P35 as a critical period for maturation in this system (lines 64–65), so that the rationale for focusing on this window is better supported.

      Regarding dendritic morphology, we agree that examining the arborization of NAc MSNs would provide a useful structural correlate of the developmental processes discussed in this section. We did not collect tissue suitable for morphological analysis from these animals, since the entire NAc punch was used for nuclear isolation, and this analysis would therefore also require a new cohort. We will likewise note this as a direction for future work, and we have revisedthe corresponding statements in the Discussion so that they are presented as hypotheses rather than as established consequences of isolation (lines 389–391).

    2. Note: This preprint has been reviewed by subject experts for Review Commons. Content has not been altered except for formatting.

      Learn more at Review Commons


      Referee #3

      Evidence, reproducibility and clarity

      Summary:

      In this study, the authors were investigating the effect of juvenile social isolation (jSI) on the nucleus accumbens (NAc) transcriptome in female mice. They used P21 wild type (C57BL6) female, isolated at P21 or group housed, and reunited at P35 to generate their samples for RNAseq on NAc punches. They also performed FACS sorting on the nuclei from NAc lysates to select only neuronal nuclei (NeuN staining). They studied the differentially expressed genes (DEGs) between group housed and jSI by RNAseq. Then, they studied the protein/protein interaction via stringDB on the DEGs identified (up or down) and perform GO analysis on them. They identified Ntrk2, Grin3a, Grik1 and Bcl2; associated with the neuronal function or transcription regulation terms. They also studied the histones modifications (H3K4me1, H3K4me3, H3K27ac, and H3K27me3) after jSI and identified neuronal function and transcription regulation terms again on the DDR (Cut and Tag method). They found that these histones modifications could play a role in jSI-induced adaptations and neuronal function. Finally, they reanalyzed public datasets of RNAseq data to identify histone modifications associated with their DEGs of interest, and compare their DEGs to differences between genotypes in the public datasets (in conditional KO models of some genes of interest, as Kdm6b cKO, BET inhibitor, or Setd1a +/- mice with 3 different mutations. They conclude that histone modification could be involved in jSI-induced gene expression alteration.

      Major comments:

      • Title: Add the sex : "in female mice" since it is specific.
      • Referencing: Biorender.com has been used to generate some schematics in this study, but it is never acknowledged or referenced.
      • OPTIONNAL: Males: Experiments are only focused on females here. The authors didn't state why. These experiments (RNAseq, Cut&Tag) could be performed independently in male mice. The introduction refers to several pathologies that show sex-bias in human or animal models and even articles with sex differences. If the authors had a reason to select female only, it should be clearly stated. This proposition would take the same amount of resources and time that for this issue but would increase the knowledge about the (sex differences in the) effect of juvenile social isolation greatly.
      • OPTIONNAL: Behavior abnormalities: Showing behavior abnormalities after jSI of female mice or, if it has been published somewhere else previously, a general description of the phenotypes observed in juvenile and adult female mice. This experiment could be performed in one batch of female separated at P21 and regrouped at P35, with anxiety (openfield or elevated plus maze, 1 day each), sociability (direct or 3-Chambers, 1 day each) and eventually depressive-like behaviors/anhedonia (sucrose preference test, 1 week including habituation to the two bottles and/or tail suspension/forced swimming, 1 day) tested. Same tests would be performed in juveniles or in adults.
      • OPTIONNAL: Datasets selection: The public datasets used through this study are far from the original experimental design proposed in this study (cerebellum at P14, cortex at E16.5, nonspecific inhibitor, whole adult PFC = 12-14weeks old). It is important to note that the dataset used, from Chen et al 2022 (whole PFC) also studied the striatum of heterozygous Setd1a mice in the same paper. Why did the author reanalyzed PFC data instead of striatum, which would probably look more like their NAc-restricted samples? Restrict the analysis of the public dataset on NAc data, if possible on females only, and /or try to obtain data from similar experimental design (social isolation, stressed mice). Note here, that the development stage during which the mice have been isolated/regrouped and sample taken will probably of importance. The reanalyzes mix embryonic, early juvenile and adult samples, none of which is consistent nor look like their set up (P31-35).
      • Wording: Formulation throughout the current study is vague, sometimes misleading, with some unclear sentences (see minor comments for the sentences showing problems). This issue needs to be checked again.
      • Assumptions are made based on 2 sets of reanalises. These parts should be displayed in the supplementary to help the authors target some genes of interest rather than the principal figures. These analyses didn't seem convincing due to too much shift from the original issue of the paper (which is jSI in the NAc in female mice).
      • Volcano plots throughout the study: Should display the genes names (at least top 10 up and top 10 down DEGs) on the graph, otherwise the volcano plots are unreadable.
      • OPTIONNAL: Include more mechanistic experiments (as suggested in the discussion) on some of the factors identified (Kdm6b, Brd4, Setd1a) to better confirm their involvement in the jSI-induced changes in transcriptome (and behavior). Conditional KO in the NAc with viral infection (mouse line Kdm6b-flox or Brd4-flox or Setd1a-flox existing + AAV-cre) or AAV Crispr for specific knockdown could hardly be performed in the time window used in this study (at least 3week expression of the Cre/Cas9). But systemic or intracerebral pharmalogical approach (ip injection of an inhibitor or cannula implantation, local intra-accumbens injection) could be performed in juvenile (1 week recovery post-surgery only, can be done at P21 before isolation).
      • Some of the top genes differentially expressed should be clearly cited in the results.

      Minor comments:

      Introduction:

      • Throughout the introduction, some parts could be ameliorated, since some sentences are unclear or confusing. Here are some examples and explanations on the issue detected (find in red comments/suggestions from the reviewer): "These negative effects are further supported by evidence from Covid-19 during the last few years » line 56/57 : reformulate.

      "In addition, isolation contributes to severe social issues, such as increased human suicide risk » line 57/58: issues is plural, but only one example is given; moreover, the term "social issue" associated to suicide is poorly-worded.

      "In the case of rodents, socially isolated animal models are proposed to be associated with various human diseases » line 59/60: clumsy sentence, the animal models are associated to human disease? This could be reformulated.

      "and the effects of juvenile social isolation (jSI) on motor, emotional, learning, and sociability-related behaviors in rodents have been widely reported (Li et al., 2021; Powell & Swerdlow, 65 2023; Walker et al., 2019), which further provides evidence of the pathogenesis and 66 molecular mechanisms of human mental disorders. » line 63-66: Examples of behavioral dysfunctions would be appreciated here.

      "The nucleus accumbens (NAc) is a critical component of the brain reward 68 circuitry, and dysfunction of the NAc is associated with drug addiction (Zinsmaier et al., 69 2022), impaired social interaction (Pomrenze et al., 2022; Shan et al., 2022), and 70 abnormal emotion expression (Gebara et al., 2021) » line 67-70: here, the statement reads as dysfunction of the NAc is responsible of abnormal behaviors (addiction, social behavior or emotional expression), but the articles show that NAc is dysregulated in models of these pathologies. Is the dysfunction in the NAc responsible of or a consequence of the pathologies? This could be better formulated.

      "An fMRI study showed that activity of the human NAc is associated with the sense of loss (Cooper et al., 2009; O'Connor et al., 2008). » line 70-72: Sentence says one study, but two rfereneces are used. The sentence refers to O'Connor only. Cooper is about reward/effort and NAc activity, not grief/loss, reformulate. Also, "activity" is unclear, the authors could be more precise with "hyperactivity of the NAc has been found in people suffering from loss".

      "The NAc from lonely individuals showed key differentially expressed genes (DEGs) that are associated with both neurodegenerative and neuropsychiatric diseases » line 74-75: Unclear, give examples of the pathologies here to be consistent with the next sentence about female rats (Alzheimer, Parkinson, Huntington). The next paragraph (line79-95) about DNA methylation and histone modifications is missing a general conclusion: what is interesting or needs to be more studied? Also, H3K9 and H3K79 have been introduced but unused in the paper, while H3K27ac/me3 have not been introduced. What is known about them?

      « How do epigenetic elements mediate gene expression dysfunction under jSI stress? In this study, we aimed to reveal the alterations in gene expression and histone modifications induced by jSI, and to elucidate their roles in the context of psychiatric disorders promoted by jSI » line 96-99: Statement is too general, it is missing the term "NAc" here. - General comment: Since the paper is focused on female, it could be of interest to state in the introduction if/how sex differences exist in relation to social isolation and human diseases showing isolation as a phenotype.

      Methods:

      • Female mice only have been used in this study. It is never explained why so. Knowing that many diseases show sex differences in prevalence or symptom expression, it would have been helpful to include also male mice in the study, to identify common mechanism linked to jSI versus sex-specific alterations.
      • Cut & Tag: Dilution of antibodies is not displayed ("1µL") and the reference for antibodies is unclear: "primary antibodies (H3K4me1, MABI, 536 MABI0302; H3K4me3, abcam, ab8580; H3K27me3, CST, 9733S; H3K27ac, CST, 537 8173S; 1 µL per reaction) ». This should be adapted to look like the FACS antibody description: "anti-NeuN-488 conjugated antibody (Millipore, #MAB377X, 1:400 dilution) ».
      • "Frozen nuclei were thawed and bound to concanavalin A (ConA)-coated magnetic beads (BioMag®Plus Concanavalin A, 10 µL per reaction) for 10-60 min on a rotator at room temperature. » Why so much difference in incubation time here?
      • Number of animals: While reading the manuscript, it was unclear that 5 different groups of mice have been used. The number of animals should be stated in the methods in the "nucleus extraction" part or "animals" section to facilitate understanding the methods.
      • Data analysis: "For RNA-seq data, p-value < 0.05 and Fold Change > 1.2 were used as the threshold for identifying DEGs. For CUT&Tag data, p-value < 0.05 and Fold Change > 2 were selected as the standard for identifying DDRs. » Cut&Tag p-value and threshold is written in RNAseq section, move it to its proper part hereafter "CUT&Tag data analysis ».
      • Suggestion DDR: acronym is present in the methods but not explained. It is however explained in the results. This depends on the order in the publication, but if methods appear first, it would be helpful to understand what stands for DDR.
      • Cut&Tag data analysis: "The procedures of quality check and trimming were the same as RNA-seq data analysis. » line 591; "and the removal of blacklisted regions was the same as RNA-seq » line 594; "GO analysis was the same as RNA-seq analysis. » line 599 : These parts could be ameliorated to avoid repetition. Since the preparation of nuclei and most of the analysis are the same, the methods could be more straightforwardly explained separating common preparation from specific analysis.
      • Methods explaining how the authors performed the reanalysis of the public datasets is missing.
      • Animals have been separated at P21 and regrouped at P35: It is not stated if they were regrouped together or with a group of unstressed WT never separated, which could influence their behavior and stress levels.
      • No behavioral test has been performed on these mice. It would have been appreciated to see that 2-weeks social isolation was efficient to generate stress in these animals (anxiety test, sociability at least). And if this protocol has been previously used in their lab, at least to explain briefly what behavior abnormalities the jSI was inducing.

      Results:

      • Section 1: Only 1 or 2 GO terms (down / up DEGs) are described, but top5 is represented in the figure. Description of the others would be of interest (for example actine reorganization could be particularly interesting). Also, citing some DEGs from the top UP and DOWN, representative of the GO terms, could be of huge interest here.
      • Figure 1, E/F/G: Some genes in Fig1G are not from top10 nodes in DEGs (Slc17a8, Hcrtr2, Dkk3, Dact1, Nr4a1, Fosl2, Htr5a). They are stated as "potentially important genes" in the legends and are found later in the study as important using other methods than RNAseq. Since Fig1G is displaying RNAseq results, it would be better to stick to the Top10 genes displayed here and put in another figure the other "potentially important genes". What means "potentially important" ? Why these ? What are the criterion? Also, the fig1G is unclear visually: separate it in two for top10 down and top10 up, it would be easier to understand and navigate. Legends: Statistics used for 1G are not stated. Volcano plot: Top 10 genes UP/down could be displayed on the graph.
      • Section 2: The first paragraph here is about which Transcription Factor (TF) is predicted to participate in their DEG's expression. Half of this paragraph is introduction about the function of several TF. This is not part of results and should be moved appropriately in the introduction or discussion section, or shortened significantly, since it is now longer than the result part.

      Importantly here, the analysis is done on ChIP Atlas (public datasets). They state : "Taken together, promoter analysis of DEGs suggests that potential epigenetic mechanisms may act upstream of jSI-induced transcriptional dysregulation in the NAc. » line 176, but the database has never been stated to be NAc-only data nor data from jSI animals. If not, this sentence has to be modified. It was unclear globally if this part was based on their work or data mining on a first read, it should be more clearly stated at the beginning that this is exploratory. - Figure2 : A/B/C: only one mention of the NAc in the data; Fig D: 5/8 NAc datasets. The analysis here seams unbalanced. Why not take into account only NAc datasets ? The composition of cortical area or retina is highly different than the NAc (mostly Glutamatergic vs GABAergic populations). If doable, the analysis focused on NAc datasets would be better. - Section 3: "First, neuronal development-related genes, such as "nervous system development", were found in all DDRs of the four histone modifications » line 193-194: sentence is unclear, the author probably meant "term".

      No gene have been cited here in any histone modification experiment (nor visible in the figure, only dots without names, top 10 up/down could be displayed on the volcano plot). It would be of interest to state at least some of the genes identified here (DDRs, closest loci) and to see if / how many common genes from RNAseq data were found again here.

      GO terms: again, only one or two examples are described, but figure shows the top5. They all could be at least stated.

      The conclusion of the first paragraph states : "These results were consistent with the transcriptome analysis that neuronal function and transcription-related genes were affected, and with the transcription factor analysis that epigenetic regulators were predicted to bind to the promoter regions of these genes. » line 197-199: Since the author did not state any genes, we can only believe that the result are consistent based on two vague GO terms "neuronal system development" and "regulation of transcription by RNApol II/chromatin remodeling". If the lector has to read itself every gene table to know which genes are dysregulated in jSI, this study will be really time consuming. - Figure3: Volcano could display the top10 names of DDRs.

      "The results indicated that down-DEGs were associated with H3K4me1, H3K4me3, and H3K27ac. » line 203: In which direction are altered H3K4me1, H3K4me3, and H3K27ac ? This is important to know. "Consistent with their active roles in transcription, downregulation of H3K4me1 and 205 H3K27ac was more relevant to down-DEGs than up-DEGs. » line 205: Why ? Unclear statement.

      "Considering the composite roles of H3K4me3 (an active histone modification) and H3K27me3 (a repressive histone modification), we hypothesize a major role of H3K27me3 in these up-DEGs, and the contribution of H3K4me3 to gene expression alteration by jSI might be small, though we cannot exclude the possibility that it regulates certain genes locally or plays a repressive role. » line 208-211: It is very unclear here, why H3K27me3 should play a major role while H3K4me3 alteration "might be small". This has to be further discussed.

      "For top 10 nodes among up-DEGs, we didn't find any significant alterations in any of the four histone modifications around their gene loci, except for downregulated H3K4me3 around Aldh18a1, Lamp1, and Gnb4 » line 217-219: Formulation is clumsy here, reformulate.

      "Some of these genes were marked by multiple altered histone modifications." Line 223: Which ones? Ony Bcl2 displayed, but the authors state "some of these genes" right after writing "H3K27ac was found to be 222 downregulated around Grin3a, Grik1, and Adgre1 ». Are these the other genes showing several histone modifications? It is unclear.<br /> - Figure 4B : only Grik1 as an example. Why only this one and not Bcl2 that moreover show several modifications? Could be helpful to show an example of each modification. -Section 4:

      Kdm6b:

      "To examine the possible contribution of Kdm6b to jSI, we re-analyzed the RNA-seq data from Kdm6b-knockout in the published study (Ramesh et al., 2023). » line 240-241: this study is about conditional Kdm6b KO in the cerebellum, on naive P14 male and female mice's neurons in culture. The authors extrapolate the results from a completely different neuronal population/region and sex to justify the potential effect jSI could have on their adolescent female mice. This sentence is misleading for the reader, since the model used (not stressed) and experimental conditions are far from what they are studying. This sentence needs some reformulation to better explain their goal. They show the DEGs (up/down) from reanalyzed data and overlap between these DEGs and the one from Figure1, but the conditions are far from each other here. One could ask what the specificity of their overlap demonstrated here.

      "Fosl2 and Nr4a1 are immediate early genes (IEGs) in response to neuronal activation in many brain regions (Dave et al., 2025; Shi et al., 2024), and these two genes have been reported to be involved in memory maintenance (McNulty et al., 2012; Mizuno 257 et al., 2020) and Parkinson's disease (PD) (Fan et al., 2020; Rouillard et al., 2018). In addition, Htr5a, which encodes serotonin receptor 5A, was also upregulated in jSI and downregulated by Kdm6b KO (Fig. 1G, 5G, Table S1). And Htr5a has been reported to be a risk factor of human schizophrenia (Guan et al., 2016). » line 254-260: This part of the result paragraph is about introduction/discussion again. This should be move appropriately.

      Brd4:

      Here, the authors reanalyzed data from E16.5 cortical neuronal culture treated with or without BET family inhibitor, which they state is not selective of Brd4 (even if it is part of the BET family). The crossover between this embryonic cortical neuronal population treated with nonspecific inhibitor and their model (juvenile Social Isolation, NAc) is a bit of a stretch. What do the overlap in DEGs really means here? "We also examined the possible downstream Brd4 target genes within the gene sets of down-DEGs by jSI and down-DEGs by JQ1 treatment, and we identified Hcrtr2 and Dkk3 in these gene sets (Fig. 1G, 5H, Table S1). Hcrtr2 encodes an orexin receptor, and it has been reported to be involved in altered arousal levels through dopamine neurons (Bandarabadi et al., 2024). Dkk3 inhibits Wnt signaling and is reported to be related to anxiety and memory formation (X. Chen et al., 2025; Flores et al., 2024). » line 275-280: What is the conclusion on these results?

      Setd1a:

      Here, they used 3 separate datasets: whole PFC of Setd1a heterozygous mice (exon 4 LacZ/Neo cassette insertion), whole PFC from loss of function Setd1a heterozygous mice and FoxP2+ nuclei from PFC of Setd1a +/- mice (frameshift in the 15th exon). These datasets are quite different between themselves and compared to NAc samples from jSI mice. The authors stated that the first two datasets had a low DEG overlap with their samples but continued with the third which showed a significant overlap for Hcrtr2, Dkk3 and Dact1. They finally conclude that: "Taken together, these results suggest that epigenetic factors, such as Kdm6b, Brd4, and Setd1a, may mediate jSI-induced gene expression alterations. » Nothing in these datasets is comparable to what they want to prove here, it is a huge stretch to propose these genes as mediators of jSI. Reformulate. These results could be exploited as exploratory, to reduce the number of potential targets, but needs to be investigated on their own.

      Figure 5: volcano plots: top10 genes visible could be useful. This section of the results would fit better displayed in the supplementary, since they reanalyzed datasets far from their experimental conditions. These genes of interest should be further investigated in their jSI model.

      Discussion:

      "For example, the expression of glutamate receptors is reduced in the NAc, prefrontal cortex, and hippocampus under isolation stress (Hermes et al., 2011; Mao et al., 2022; Sestito et al., 2011). » line 315-317: Which GluR are reduced here ? It needs to be more precise for the reader here, and to state if some genes as been found in common between this literature and their DEGs.

      « The NAc is a key component of the brain reward circuit, and it is involved in drug addiction and social behavior (Pomrenze et al., 2022; Zinsmaier et al., 2022). NAc neurons receive glutamatergic inputs from the PFC, basolateral amygdala (BLA), hippocampus, and ventral tegmental area (VTA) (Arrondeau et al., 2024; Dieterich et al., 2021; Elam et al., 2025; Le Borgne et al., 2025; Zinsmaier et al., 2022), and neurons in the NAc output the information to the ventral pallidum (VP) (Liu et al., 2022), VTA (Qi et al., 2022), and other areas of the basal ganglia (Lanciego et al., 2012). » line 318-324: These lines are describing the circuitry of the NAc, some of its inputs (no mention of dopamine afferences from the VTA) and outputs. No use of this information is used after, since they conclude the paragraph with: "Thus, deficits in glutamatergic synapses possibly mediate jSI-induced behavioral abnormalities, including impaired social interaction, anxiety, and an increased risk of substance abuse ». line 325-326: What is the point of describing the circuit, if it is not interpreted regarding their results? What is their hypothesis on the circuit dysfunction in jSI female mice? They were discussing the DEGs from RNAseq result before this paragraph. What is the link/hypothesis between their DEGs and the glutamatergic circuits of the NAc? Is the NAc directly responsible of jSI-induced behavioral abnormalities for them or cortical/amygdal/hippocampal/VTA glutamatergic projection neurons are dysregulated, creating DEGs at the synapse in the NAc? This part of the discussion should be more specific on what they mean.

      "Deficiencies in these proteins are associated with various behavioral abnormalities (Araujo et al., 2017; Chasse 335 et al., 2024; Guo et al., 2020; Huang et al., 2021; Mukai et al., 2019). » line 334-335: what proteins and what behavioral abnormalities? This is not precise enough and needs reformulation/conclusions. "A previous report suggests that histone modifications such as H3K4me3 in the hippocampus respond to an enriched environment (Schaffner et al., 2023), and our results indicate that these histone modifications may influence gene expression in the NAc under jSI stress as well. » line 342-344: In which direction is the modification in the hippocampus in enriched environment? Is it opposite to what the authors have found in jSI (which would be interesting, since one could see a more social environment as an enriched condition too)? The idea behind this sentence needs to be precised.

      "Since our mice were isolated from P21 to P35, a period critical for the maturation of neurons (Makinodan et al., 2012; Walker et al., 2019; Yamaguchi et al., 2024), it is possible that the neuronal development process is affected by the isolated housing environment. » line 351-354: This sentence states that juvenile social isolation during development could impair development. It is probable, since early life stress (even maternal stress) can have prolonged effect on the offspring behaviors. It is probable that the time-window of jSI is important for development. OPTIONNAL: The authors would beneficiate of more experiments here: 1) They could perform SI in the same conditions but in adult females and compare the DEGs observed in that case (less development related genes probably). 2) "Neurons during adolescence mainly experience synaptic pruning and elimination (Afroz et al., 2016; Germann et al., 2021; Watanabe & Kano, 2024), and isolation stress probably impairs such processes. ". Here, the author could check in the NAc of their jSI female mice the state of dendritic arborization of MSNs (number, length, etc) on NAc brain slices.

      ". Since neural development relies on the regulation of gene expression (Jain 370 et al., 2001; Xiang et al., 2020), we hypothesize that these terms reflect altered gene expression regulation mechanisms under jSI stress » "However, these hypotheses need to be further validated by additional experiments. » line 369-371 & 375-377: The author are not integrating their results, they are being cautious, but the message stays unclear to the reader. What is the message here ?

      « Besides these two main shared functions affected by jSI, our results suggest that other biological processes are potentially mediated by one or more histone modifications. For example, some DDRs of H3K27ac and H3K27me3 are functionally enriched around cell adhesion-associated genes, and this is consistent with previous papers suggesting that cell adhesion is affected by isolation (Santiago et al., 2023; Wu et 382 al., 2022)... » line 377-382: This GO term appeared in the figure, but has never been mentioned clearly in the results. The explanation goes on for a full paragraph. It could be better to introduce it before if it is of interest. Also, which cell-adhesion genes have been found in the RNAseq / cut&tag experiments for this family of genes (never stated)?

      "To determine whether histone modifications regulate specific genes, we focused on potentially important genes. Grik1, for example, exhibits reduced H3K27ac levels. It encodes a subunit of ionotropic glutamate receptors, and its deficiency has been found in mental diseases, such as schizophrenia and ADHD (Chatterjee et al., 2022; Hirata et al., 2012). The inactivation of Grik1 in rodents promotes anxiety-like behaviors via glutamatergic transmission (Englund et al., 2021). » line 387-392: What is the conclusion/hypothesis on Grik1's role ?

      "Bcl2, for example, has downregulated H3K4me1, H3K4me3, and upregulated H3K27me3 levels. Bcl2 is an apoptosis-related gene that determines neuronal survival under stress. A previous study suggests that chronic social defeat stress decreases the Bcl-2/Bax ratio in NeuN+ neurons in the hippocampus (Zhu et al., 2024). Our data suggest that jSI is another type of stress that suppresses Bcl2 expression, and that the epigenetic factors are possible upstream regulatory mechanisms » line 393-399: No links or clear hypothesis have been made here. The authors proposed to go deeper in this direction later. It would be interesting to conclude on the hypothesis on these two genes (Grik1/Bcl2) in their model.

      "We found that Kdm6b may be at least partially involved in the gene expression changes in jSI mice, especially for potentially important genes like Nr4a1. Nr4a1 encodes a transcription factor that regulates dopamine metabolism, and it is reported to be involved in drug addiction, which is probably mediated by histone modification alterations » line 408-411: The authors are very cautious about their conclusion, maybe too much. Since they base their hypothesis on P14 cerebellum neurons in culture, more work is needed here to establish clearly the role of Kdm6b. It feels like the authors are dropping clues for the reader, without concluding themselves on their hypothesis.

      "Our experiments suggest that Dkk3 is regulated by jSI stress, and this is probably mediated by Brd4.» / "Our analysis also implies that Dkk3 and Hcrtr2 are probably regulated by other epigenetic regulators such as Setd1a. "line 429-430 & 431-432: "Probably". No proof is given on that statement. Everything is "potential" or "possible" or "probable" here.

      Summary paragraph: "To summarize, we revealed the changes in gene expression and histone modification levels under jSI stress. » line 444: Here they are missing the term "in the NAc" and "in female mice" and maybe need to add "some changes" in the sentence.

      Limitations paragraph: Major issue of this study is stated as "Another limitation related to this is that we inferred candidate epigenetic factors based on previously published public data. However, different brain regions, cell types, and experimental conditions between our analyses and public datasets may contribute to the gene expression differences, leading to false negative or false positive results in this study. » line 456-460: Indeed, as they highlight here, the differences between the datasets chosen and their initial context (jSI, NAc) is huge. Maybe the authors could explain better why did they choose these datasets instead of more similar ones. Also, in the RNAseq / cut&tag analysis, they "were unable to separate these subtypes (D1+ or D2+) for this current analysis ». This analysis could be interesting to see in the supplementary or a brief description of what they have tried in the discussion, since the composition of the NAc is made of about 90-95% of MSNs and a plethora of interneurons (Parvalbumin, Cholinergic, etc), leading to different circuit connectivity and function. The striatum also contains a third population of D1/D2 hybrid neurons, maybe this population (about 5-10% of the whole striatum) made the identification of the neuronal subtypes harder.

      Significance

      Globally, this study is showing transcriptomic and histone modifications occurring after juvenile social isolation in female mice. The authors identified differentially expressed genes linked to neuronal development, regulation of transcription and chromatin remodeling. The reanalyzed public datasets to identify histone modification on the DEG identified and observed the impact of some already published mutations on the gene expression to compare it to their data. The limits of this issue are caused by the reanalysis parts, since the datasets used are cortical and cerebellum samples, in diverse development stage (embryonic, early juvenile, adults - both sexes) that is quite different compared to their paradigm (juvenile females). They also never display the name of the principal DEGs identified (text or plots) which leads to difficult understanding of the findings of this paper. A more focused analysis on NAc or striatal, female only, juvenile stage datasets would be more helpful in this situation. The text should be more precise sometimes and a specific explanation on the exclusion of male mice should be introduce early in the methods.

      This study finds its place in the current research on the role of NAc function/dysfunction in behavioral abnormalities induced by social isolation and try to understand the mechanisms behind the abnormal behaviors induced by separation. The audience could be composed of researchers from several domains where social isolation is the cause or the consequence of pathological behaviors, including studies on loss, depression, ASD, Alzheimer, Schizophrenia, etc). The context is quite broad. The results from this paper could help find new molecular targets to alleviate the effects of social isolation and perhaps ameliorate the behavior for mouse models of several diseases or later in patients. A better understanding of the effects of social isolation in female is interesting, but being able to compare both sexes would be even better: identifying sex-differences and common defect is of great interest nowadays in several domains.

      The present reviewer has expertise in behavior in mice (both male and female) from juvenile to adult stages, has studied neuronal circuits including prefrontal cortex and striatum (mainly NAc) in behavioral abnormalities in mice in a model of ASD and more recently in an addiction model. The reviewer is interested particularly in sex-differences in neuronal circuits defects and behavior expression in diseases. Finally, the reviewer has recently focused on spatial transcriptomic approaches in addiction models.

    3. Note: This preprint has been reviewed by subject experts for Review Commons. Content has not been altered except for formatting.

      Learn more at Review Commons


      Referee #2

      Evidence, reproducibility and clarity

      This paper profiles NeuN positive NAc nuclei after juvenile social isolation. The authors report RNA seq changes and CUT and Tag for H3K4me1, H3K4me3, H3K27ac, and H3K27me3. They then link DEGs to public datasets for Kdm6b, Brd4, and Setd1a. The neuron enriched design is useful. The main claim remains correlative. Causal support is thin.

      Major points

      Causality is not shown. The Discussion states: "Although we didn't show the molecular mechanism of histone modification alteration regulating gene expression, we revealed the association between transcriptome and histone modifications." That limit should shape the Abstract and title more clearly. Phrases like "epigenetic alterations may also play a role" are fine. Stronger wording about mediation should be toned down until NAc specific perturbation is done.

      DEG thresholds are loose. DEGs are defined as "p-value < 0.05 and Fold Change (FC) > 1.2." There is no clear FDR cutoff. With ~1250 DEGs from n = 5 to 6, many hits may be noise. Please report FDR filtered lists or justify the uncorrected p value choice. Re run key GO and overlap tests on a stricter set.

      Public data overlaps are hard to interpret. Kdm6b data are from cerebellum. Brd4 data are from cultured cortical neurons treated with JQ1. Setd1a data are mostly PFC. The authors note that "different brain regions, cell types, and experimental conditions... may contribute to... false negative or false positive results." That caveat is important. Overlaps should be framed as hypothesis generating only. Do not treat them as evidence that these enzymes act in NAc under jSI.

      Behavior is missing from this study. The Introduction says jSI affects "motor, emotional, learning, and sociability-related behaviors." This manuscript does not show that the molecular changes track those phenotypes in the same cohort. Without that link, the psychiatric disease framing stays speculative.

      CUT and Tag analysis is coarse for promoter claims. Signals are quantified in "all 5 kbp bins." That bin size can blur promoters, enhancers, and neighboring genes. Please add peak calling or TSS centered analyses for key loci such as Grik1, Bcl2, and Dkk3. Also show more browser tracks beyond one example.

      Multiple testing for overlaps needs attention. Many Fisher tests compare DEGs with DDRs and with several public DEG lists. Report whether p values were corrected across tests. Some reported overlaps are small in absolute numbers even when p values look significant.

      Minor points

      Only female mice were used. State this early and discuss sex limits. Juvenile isolation effects often differ by sex. Down DEG GO terms include "Chondrocyte differentiation" and "Positive regulation of cartilage development." These look odd for NAc neurons. Check annotation quality and whether these terms survive stricter DEG filters.

      Figure 1 lists "II2ra" in the top nodes table. That is likely Il2ra. Please correct. Sample sizes differ a lot across marks. H3K4me1 and H3K27me3 have n = 4 in jSI. Discuss power and why replicates differ.

      The isolation protocol includes regrouping from P35 to P49. Make clear that effects are lasting post isolation effects, not acute isolation effects.

      Methods say "GPT-5.4... and Claude Sonnet 4.6... was used." Fix subject verb agreement. Data Availability lists "GSE3508789." Confirm this accession. It looks malformed.

      Abstract keywords include "Loneliness." The mouse work is social isolation. Keep that distinction clear, as the Introduction already does.

      Overall

      Solid descriptive resource on NAc neuron transcriptome and histone marks after jSI. Not yet strong enough for firm mechanistic claims about Kdm6b, Brd4, or Setd1a. Tighten statistics, soften causal language, and add locus level epigenomic detail. Functional tests in NAc would raise impact a lot. At present I see this as useful but preliminary.

      Significance

      This study provides a useful neuron-enriched transcriptomic and histone modification resource from the nucleus accumbens following juvenile social isolation. The integration of RNA-seq and CUT&Tag data adds value for researchers studying epigenetic regulation and stress-related neurobiology. However, the advance is primarily descriptive rather than mechanistic, as the conclusions rely largely on correlative analyses without functional validation. The manuscript will be of interest to the neuroepigenetics and psychiatric neuroscience communities, but the conceptual advance is incremental, and the mechanistic claims should be moderated.

      My expertise: Single-cell and bulk transcriptomics, epigenomics, neuropsychiatric disorders, and computational genomics.

    4. Note: This preprint has been reviewed by subject experts for Review Commons. Content has not been altered except for formatting.

      Learn more at Review Commons


      Referee #1

      Evidence, reproducibility and clarity

      The manuscript by You et al. investigates changes in gene expression and histone modifications after juvenile social isolation (jSI) in the nucleus accumbens (NAc). They find many differentially expressed genes and find overlap with activating or repressive histone marks. They then go on to compare their data to other published datasets to support their findings. This is an interesting study, and the authors use creative approaches using unique and published data to identify epigenetic mechanisms underlying the changes in gene expression within the NAc. However, there are many points of clarification that are needed to fully evaluate the manuscript and several experimental details.

      1) The text in the introduction conflates adult and adolescent social isolation which have very different effects on behavior. Additionally, there is evidence that isolation during adolescence can have permanent effects on behavior but the behavioral effects of adult isolation in rodents are transient. It is recommended that the authors restructure the intro to be more specific to describing the adolescent period and why epigenetic mechanisms would be expected to reguate changes induced by jSI.

      2) The authors cite several studies from the Nestler lab on early life stress that link early life stress to histone modifications but failed to cite the manuscripts that investigated the transcriptional changes in response to jSI. These studies also highlight sex differences in jSI. This is important given that this study only uses females. Many of the effects described might not be comparable simply because of the sex of the animals.

      3) Can the authors clarify if they used an adjusted p-value or nominal p-value. If they are using a nominal p-value the authors should explain their reasoning and provide information regarding if any of the transcripts survived a p-value correction. The addition of threshold-free approaches are more appropriate (GSEA) rather than focusing on transcripts with a nominal p-value. If they are going to present data using a nominal p-value, this should be justified and the cut off should be explained and every interpretation should include a caveat.

      4) It is unclear how the authors confirmed that input RNA or neurons were similar across samples for the library prep. This is especially important given the top genes that are differentially expression. The finding that beta actin (Actb) is up in jSI vs GH animals. The results could be due to differences in input rather than actual differences in expression.

      5) There are aspects of the methods are difficult to understand. For example, under RNA-seq, the authors mention "frozen nuclei were thawed and centrifuged.....the supernatant was discarded and nuclei were centrifuged again under the same conditions" Can the authors clarify what was done here? Were the nuclei resuspended in STEM CellBANKER or something else? While this is a concrete example, there are many other places in the methods the are like this, meaning that steps seem to be skipped and it then becomes difficult to assess the approach. It is recommended that the authors work to clarify the methods. Another example, how the DNA was treated in the CUT&TAG and how much DNA was added to the library prep.

      6) Can the authors please explain why only females were used for these experiments? In addition, can the authors please comment on potential caveats in the interpretation by only including females in the study?

      7) Were females shipped to the facility on P21? It is unclear.

      8) Were any animals used for multiple endpoints or was each endpoint a separate cohort? Were any samples pooled?

      Significance

      This is an interesting study and the authors use creative approaches using unique and published data to identify epigenetic mechanisms underlying the changes in gene expression within the NAc. However there are many points of clarification that are needed to fully evaluate the manuscript and several experimental details are missing or unclear.

    1. eLife Assessment

      This manuscript presents important work on macrostructure formation in a freshwater filamentous cyanobacterium, focusing on its ability to aggregate and buckle. The authors employ a wide range of experiments and approaches, including time-lapse imaging and 3D theoretical modelling, to establish the physical dynamics of filaments, such as gliding motility and buckling. They demonstrate that, in addition to gliding motility, filament length and flexibility are essential for the ability of cyanobacteria to capture particles and form aggregates. Overall, the study provides convincing evidence and advances our understanding of the role of microbial motility in the environment. It will be of broad interest to biophysicists and environmental microbiologists.

    2. Reviewer #1 (Public review):

      Summary:

      In this manuscript, the authors investigate the mechanisms underlying macrostructure formation in a freshwater filamentous cyanobacterium strain, F. draycotensis, focusing on how its ability to aggregate and form these structures depends on the physical properties of the filaments. Using experimental observations, they demonstrate that the cyanobacterium actively captures and surrounds particles, a process driven primarily by gliding motility.

      To explain these physical dynamics, the authors present a 3D model indicating that particle collection relies on filament length, as well as a specific mechanical response, namely, filament buckling and the subsequent formation of loops of bundles of filaments. While the authors have previously documented the buckling and looping characteristics of this strain, this study provides new insight by demonstrating that these physical phenomena are essential for particle capture and collection.

      Strengths:

      This manuscript benefits from a rigorous and detailed quantitative analysis of video recordings, which clearly documents the motility, buckling behaviour, and particle collection dynamics of the filaments.

      The authors effectively validate their hypothesis by using a naturally shorter filamentous strain, which fails to collect particles, suggesting that filament length is indeed a critical parameter.

      To further confirm the length dependency within the same species, the authors experimentally generated shorter filaments of F. draycotensis. The fact that these shortened filaments also lose the capacity to collect particles provides strong evidence supporting their proposed mechanism.

      Weaknesses:

      There is a conceptual concern. The authors linked the specific physical properties of this strain to evolutionary data, highlighting that the studied lineages diverged approximately two billion years ago. This creates a misleading impression that particle collection via flexible looping filaments is a recent evolutionary adaptation. However, particle collection has been observed in other cyanobacteria, such as Trichodesmium, which features short, rigid filaments. Therefore, the term "emerging" does not seem appropriate for the title and text. The capacity to collect particles in the studied strain F. draycotensis appears to be primarily a function of physical characteristics (filament length and flexibility) rather than evolutionary age. Any cyanobacterial strain possessing similar physical properties is likely to exhibit comparable behaviour, rendering the evolutionary timeframe largely irrelevant to the core mechanism. In addition, the phylogenetic tree presented in Figure S5 does not reflect the current consensus on cyanobacterial evolution and systematics and does not align with modern phylogenomic frameworks (see, for example, Strunecky et al., 2023 https://doi.org/10.1111/jpy.13304). There is also no such order Cyanobacteriales, which has been mentioned in a few older publications but is clearly outdated.

      Another concern is that the authors nearly completely ignore the role of type IV pili in the gliding motility of cyanobacteria, including filamentous strains. For a long time, there was a misconception that the gliding motility of cyanobacteria was due to slime protrusion. Slime plays a role in this process. However, several studies have shown that filamentous strains also use type IV pili to glide on surfaces. The authors should discuss this and include it in their model. In addition, the authors concluded that gliding motility is responsible for particle collection by Fluctiforma draycotensis. Although I believe that their conclusion is correct, there might be several limitations to the experiments which allow for other reasons to be considered. Their conclusions were based on the use of a non-motile strain and an unspecified community without the motile Fluctiforma draycotensis strain. The problem I see here is that it is not clear why this strain is not motile; it could be because of the lack of type IV pili, mutations which alter their functionality, defects in slime secretion, any other mutation (e.g. in chemoreceptors), cellular structure, metabolism, or combinations of these. Furthermore, it is possible that the community changes its composition and behaviour when it lives without the cyanobacterium with a rich carbon source (glucose) or with a non-motile cyanobacterium which may not secrete slime or, for example, a signalling component which controls behaviour of the bacteria in the community. For that reason, the authors should be more cautious with their conclusion that solely motility behaviour of Fluctiforma draycotensis is responsible for particle collection. Additional factors might be responsible for these effects.

    3. Reviewer #2 (Public review):

      Summary:

      The authors studied aggregation, buckling, and particle collection by the filamentous cyanobacterium Fluctiforma draycotensis, as well as by the filamentous Pseudanabaena sp. (order Pseudoanabenales). They performed a range of experiments, from imaging individual gliding filaments to multiple-day experiments showing the formation of large aggregates around a particle formed from a precipitate. They also developed a model of buckling filaments to argue that the ability of elastic filaments to collect particles and form macrostructures is confined to a part of the filament phase space in terms of length and flexibility, meaning that gliding combined with certain filament length and flexibility naturally reproduces the observations.

      Strengths:

      This is an impressive study that uses multiple tools to connect macrostructure formation with filaments' gliding motility and buckling. It adds an important perspective on the biological and physical factors at play in the emergence of aggregates.

      Weaknesses:

      The authors ignore the possibility that filament behavior plays an important role in the emergence of the observed patterns. Cyanobacteria have been shown to control their gliding motility (Pfreundt et al Science 2023; Kurjahn et al Nature Comm 2024), and their molecular motors are known to be regulated by chemotaxis-like signaling pathways (Risser ARM 2025). As far as I know, how the coordination between the pulling agents along an individual filament works is actively debated, but there seems to be little doubt that it exists. 

      To illustrate this point better, note that the aggregation observed by the authors is consistent with the length-dependent ability of filaments to coordinate gliding (I'm not saying this is how it works in Fluctiforma draycotensis; I'm saying it's consistent). Suppose the coordination requires sufficiently long filaments, which could be the case when signaling molecules travel along the filament, propagating information about when individual pulling agents should reverse. In such a model, short filaments act randomly because they fail to coordinate gliding by the time they glide off nascent aggregates, whereas longer filaments can perform informed reversals because they have more time for coordination. Such behavior then explains the lack of aggregation in Pseudanabaena sp. (via behavior, not lack of stiffness). Note that Trichodesmium is stiff; its filaments do not buckle, yet Trichodesmium forms organized aggregates via tightly controlled motility. Note also that, as the authors report, since Pseudanabaena sp. is both shorter and faster, its filaments have relatively (to the time needed to glide the filaments' length) little time to coordinate reversals. In my opinion, whether the observed patterns passively emerge from gliding and buckling or result from active behavior remains an open question.

      I also have a small suggestion regarding this statement on model novelty:

      The essential novelty of this model is that the filament itself is active and out of equilibrium, and additionally, the forces and torques are applied locally along its centreline, and not at its extremities as in previous steady-state mechanical studies of elastic, twistable filaments such as DNA [31-33] (see Methods and SI).

      This statement needs to be revised as it ignores a substantial body of work on self-organization of active filaments: (R. E. Isele-Holder, J. Elgeti, G. Gompper, Soft Matter 2015; Pfreudnt et al, Science 2023; Faluweki et al PRL 2023; Kurjahn et al Nature Comm 2024).

      Last point: the authors often say that their observations are reproducible ('...reproducibly forms macroscopic granules...'). What is meant? Different experiments on different days, different aliquots?

    4. Reviewer #3 (Public review):

      Summary:

      The authors report and characterize the formation of aggregate microstructures by the motile filamentous cyanobacterium Fluctiforma draycotensis, which exhibits gliding motility accompanied by rotation along the long axis while excreting EPS. In experiments with motile F. draycotensis cultures, they observed the formation of granular structures composed of cyanobacteria and other material (iron, polystyrene beads, etc.), with macrostructures on the scale of 1mm within 24 hours. The structures were motile at speeds comparable to that of the cyanobacteria filaments, resulting in their growth through coalescence over time. Notably, such macrostructures were absent in nonmotile F. draycotensis, pointing to the role of filament motility in their formation. Through experiments examining the micro-scale dynamics, inert material such as small polystyrene beads was found to be transported by the gliding, buckling, and plectoneme dynamics of the filaments, pointing to the underlying mechanism by which particles are collected into larger-scale microgranule structures.

      To interrogate the properties that drive the cyanobacteria filament buckling, plectoneme formation, and entanglement, the authors develop a mechanical model for filaments as nearly inextensible, slender bodies with resistance to twisting and bending under active gliding forces and torques and responding to fluid flows and surface adhesion. They derive expressions for the thresholds for buckling and twisting instabilities, which are additionally demonstrated and interrogated through simulation via the Immersed Boundary Method. Most importantly, bending and plectoneme formation only occur with sufficiently long filaments, and the threshold is shorter for bending than for plectoneme formation. Experimental observations with wild-type filaments agree with the model-predicted thresholds. The authors perform additional experiments with shorter filaments below both thresholds, including the filamentous bacterium Pseudanabaena, which fail to collect particles (though can in principle form macrostructures).

      Strengths:

      This work appears to be novel (notably, the discovery and characterization of the particle collection behavior of a filamentous cyanobacterium) and has interesting implications for both naturally observed cyanobacterial macrostructures as well as the controllable parameters in engineering them. The experimental and modeling work is well motivated, contributing to the broader understanding of macrostructure formation and material aggregation through active filament dynamics (not exclusive to cyanobacteria), as well as the underlying physical properties governing important filamentous cyanobacterium dynamics. As such, I would expect the results of this paper to be of broad interest to both biophysicists and microbiologists. Generally, the manuscript is well written with clear, compelling figures that illustrate the important conclusions of this study.

      Weaknesses:

      In the section on "Shorter gliding filaments cannot collect particles nor form granule macrostructures", the filamentous cyanobacteria considered *all* fall below the predicted thresholds for bending and twisting. The "long" F. draycotensis are 60 microns in length, notably less than the 120 and 320 micron thresholds derived in the previous section as well as the lengths of filaments considered in Figure 3D, yet these "long" 60 micron filaments form macrostructures. How can this be understood in the context of the model predictions? Is the nature of the macrostructures in Figure 4B, the microscale parameters, or the collection of particles somehow different than those with filaments an order of magnitude longer in earlier parts of the paper? The paper would be stronger if these sorts of questions were addressed in the text and/or with supplementary figures.

    5. Author response:

      Public Reviews:

      Reviewer #1 (Public review):

      Summary.

      In this manuscript, the authors investigate the mechanisms underlying macrostructure formation in a freshwater filamentous cyanobacterium strain, F. draycotensis, focusing on how its ability to aggregate and form these structures depends on the physical properties of the filaments. Using experimental observations, they demonstrate that the cyanobacterium actively captures and surrounds particles, a process driven primarily by gliding motility.

      To explain these physical dynamics, the authors present a 3D model indicating that particle collection relies on filament length, as well as a specific mechanical response, namely, filament buckling and the subsequent formation of loops of bundles of filaments. While the authors have previously documented the buckling and looping characteristics of this strain, this study provides new insight by demonstrating that these physical phenomena are essential for particle capture and collection.

      Strengths:

      This manuscript benefits from a rigorous and detailed quantitative analysis of video recordings, which clearly documents the motility, buckling behaviour, and particle collection dynamics of the filaments.

      The authors effectively validate their hypothesis by using a naturally shorter filamentous strain, which fails to collect particles, suggesting that filament length is indeed a critical parameter.

      To further confirm the length dependency within the same species, the authors experimentally generated shorter filaments of F. draycotensis. The fact that these shortened filaments also lose the capacity to collect particles provides strong evidence supporting their proposed mechanism.

      We thank the reviewer for the accurate summary of our work and for their identified strengths of the study.

      Weaknesses:

      There is a conceptual concern. The authors linked the specific physical properties of this strain to evolutionary data, highlighting that the studied lineages diverged approximately two billion years ago. This creates a misleading impression that particle collection via flexible looping filaments is a recent evolutionary adaptation. However, particle collection has been observed in other cyanobacteria, such as Trichodesmium, which features short, rigid filaments. Therefore, the term ”emerging” does not seem appropriate for the title and text. The capacity to collect particles in the studied strain F. draycotensis appears to be primarily a function of physical characteristics (filament length and flexibility) rather than evolutionary age. Any cyanobacterial strain possessing similar physical properties is likely to exhibit comparable behaviour, rendering the evolutionary timeframe largely irrelevant to the core mechanism.

      We would like to first clarify our use of the term “emergent”. It seems that the referee took this in an evolutionary context, whereas we are using this term in the context of its use in systems dynamics, and referring to: “a complex entity displaying behaviors that its components do not have on their own, and emerge only when they interact in a wider whole”. Here, particle collection and dynamic aggregate formation “emerges” from the buckling and interaction of many filaments.

      With regards to the evolution of particle collection behavior, our comment on the evolutionary distance between F. draycotensis and Pseudoanabena sp. was meant to highlight the point that particle collection seems to be a function of physical characteristics and motility: Despite a large evolutionary distance, and possibly many biological differences, a physics-based argument is capturing the difference between the particle collection ability of these two organisms. Thus, we are in agreement here with the reviewer. We did not intend to make any arguments about “evolutionary age” of the particle collection behavior.

      We see that the short, evolutionary comment in the Introduction has confused the reviewer and potentially is confusing to other readers too. We will therefore remove this evolutionary comment from the Introduction section of the revised manuscript and make the point in more detail in the Discussion section.

      In addition, the phylogenetic tree presented in Figure S5 does not reflect the current consensus on cyanobacterial evolution and systematics and does not align with modern phylogenomic frameworks (see, for example, Strunecky et al., 2023 https://doi.org/10.1111/jpy.13304). There is also no such order Cyanobacteriales, which has been mentioned in a few older publications but is clearly outdated.

      We thank the reviewer for this comment, as it has made us realise that we never explained our choice of taxonomic framework in the manuscript, and perhaps this is the source of the confusion.

      The order Cyanobacteriales does exist: it is the order-level name applied in the Genome Taxonomy Database (GTDB) [5, 6], currently the most comprehensive and actively curated genome-based taxonomy of prokaryotes. GTDB classifies taxonomic groups algorithmically, as monophyletic groups in a concatenated marker-protein phylogeny with ranks normalised by relative evolutionary divergence. This has produced a number of re-groupings and new names relative to the older, morphology-derived classifications; many of these have since been formally proposed under the International Code of Nomenclature of Prokaryotes and the SeqCode [2], and are progressively being adopted by the NCBI. The placement of Cyanobacteriales, and of the other orders shown in Figure S5, can be inspected directly on the GTDB “Taxonomy Tree” (see here for the orders within the class Cyanobacteriia).

      We would also like to note that we do not see our tree and the framework of Strunecky et al. as being in conflict. Strunecky et al. constructed their phylogenomic backbone using GTDB-Tk and the same 120-marker concatenated alignment that the GTDB itself uses. What differs between the two schemes is therefore not the underlying phylogeny but the nomenclature applied to the resulting clades: Strunecky et al. work within the botanical tradition and combine the phylogenomic tree with phenotypic characteristics, thereby proposing ten new orders and fifteen new families, whereas GTDB assigns rank boundaries purely by evolutionary divergence and so draws broader order limits. In practice, the GTDB order Cyanobacteriales spans several of the families (e.g. Oscillatoriales and Coleofasciculales) and orders (e.g. Chroococcales and Nostocales), that are proposed within the Strunecky et al. work. Our reason for adopting the GTDB nomenclature is for practical reasons specific to this study. F. draycotensis is a recently described organism [3] that is not included in Strunecky et al. and has no placement in their tree. In GTDB it falls within a family-level lineage (placeholder name JAAUUE01) inside the Cyanobacteriales, with the sequenced members of the Coleofasciculaceae as its closest relatives. We could not have assigned it to one of the Strunecky orders without inventing a placement. The same applies to some of the other, recent metagenomically described cyanobacteria [10], which similarly have no assigned names in the literature. GTDB, by contrast, provides a reproducible, algorithmic assignment for all of these genomes, and is now widely used for this reason in genome- and metagenome-based studies of cyanobacteria (e.g. [1]). We therefore used it consistently throughout.

      Finally, with regards to the reviewer’s point about the tree itself, we would like to note that Figure S5 was intended only to convey the evolutionary distance between F. draycotensis and Pseudanabaena sp., and it was built from a modest set of six concatenated ribosomal protein markers using an approximate maximum-likelihood method with SH-like local support values. This is considerably less rigorous than the 120-marker RAxML and Bayesian analysis of Strunecky et al., and we agree that a stronger tree may be preferable. For the revised manuscript we are recomputing the tree from a substantially larger set of concatenated single-copy marker genes, using IQTREE with model selection and non-parametric bootstrap support. We would note, however, that the specific conclusion drawn from this figure — that the two strains we use for our experimental work, namely F. draycotensis and Pseudoanabena sp. belong to deeply divergent cyanobacterial lineages — is supported by the deep backbone of the cyanobacterial tree, which is stable across marker sets and inference methods, and is equally supported by the tree of Strunecky et al.

      We will make these points clearer in the Methods and Discussion sections of the revised manuscript, as well as the Figure S5 legend.

      Another concern is that the authors nearly completely ignore the role of type IV pili in the gliding motility of cyanobacteria, including filamentous strains. For a long time, there was a misconception that the gliding motility of cyanobacteria was due to slime protrusion. Slime plays a role in this process. However, several studies have shown that filamentous strains also use type IV pili to glide on surfaces. The authors should discuss this and include it in their model.

      The reviewer is correct that we did not include molecular details of gliding motility in our biophysical model. They are also correct to point out that pili and slime biosynthesis genes are shown to be involved in gliding motility [8]. It is, however, still unclear how these factors interact to produce mechanical gliding forces that can result in filament rotation (observed only in some filamentous cyanobacteria), filament reversal, as well as decoordination during such reversals, which we have previously shown in F. draycotensis [9]. Therefore, we have chosen to keep the biophysical model at a coarse-grained, phenomenological level. Instead of explicitly modelling the detailed molecular mechanisms behind force generation, we model only the minimum necessary physical forces and torques needed to reproduce the observed rotation and translation of the filament during gliding under de-coordinated conditions. This model is able to reproduce the experimentally observed buckling and twisting of filaments, and is therefore sufficient and useful to achieve a coarse-grained understanding of mechanical forces and their relationship to buckling, twisting and entanglement, which are the main processes we focus on here. As molecular details behind force generation in rotating, filamentous cyanobacteria become available, more detailed physical models can be constructed. We also note, in this context, that the two filamentous cyanobacteria we compare both encode the type IV pilus machinery, so the presence of a pilus motor does not by itself distinguish a particle-collecting from a non-collecting strain (see our response to the reviewer’s next point).

      We will make these points clearer in the Methods and Discussion sections of the revised manuscript.

      In addition, the authors concluded that gliding motility is responsible for particle collection by Fluctiforma draycotensis. Although I believe that their conclusion is correct, there might be several limitations to the experiments which allow for other reasons to be considered. Their conclusions were based on the use of a non-motile strain and an unspecified community without the motile Fluctiforma draycotensis strain. The problem I see here is that it is not clear why this strain is not motile; it could be because of the lack of type IV pili, mutations which alter their functionality, defects in slime secretion, any other mutation (e.g. in chemoreceptors), cellular structure, metabolism, or combinations of these. Furthermore, it is possible that the community changes its composition and behaviour when it lives without the cyanobacterium with a rich carbon source (glucose) or with a non-motile cyanobacterium which may not secrete slime or, for example, a signalling component which controls behaviour of the bacteria in the community. For that reason, the authors should be more cautious with their conclusion that solely motility behaviour of Fluctiforma draycotensis is responsible for particle collection. Additional factors might be responsible for these effects.

      Our conclusion that gliding motility is the main factor underpinning particle collection is based on several observations.

      Firstly, on the macroscopic scale we present several control experiments where we did not observe particle collection: (i) in the community featuring a non-motile F. draycotensis, and with mostly the same other bacterial species as the community featuring the motile F. draycotensis, (ii) in a bacterial community derived from the original F. draycotensis community but lacking any cyanobacteria, (iii) in the original community with physically shortened F. draycotensis, and (iv) in another cyanobacterial community featuring different bacteria and a naturally shorter, filamentous gliding cyanobacteria Pseudanabaena sp. A straightforward, parsimonious explanation that satisfies all these observations is that particle collection is underpinned by physical characteristics of gliding filamentous cyanobacteria.

      Secondly and more directly, in time-lapse microscopy imaging we repeatedly observe clusters of beads being moved by gliding filaments, and thereby being collected into larger clusters. Thus, whilst factors such as slime secretion also contribute, the primary mechanism driving the observed particle motion seems to be that particles stick to filaments and are carried around with them as they glide. We cannot rule out a contribution of pili to bead attachment and transport. We note, however, that both cyanobacteria compared here encode the type IV pilus machinery. In a homology survey of the two genomes, Pseudanabaena sp. and F. draycotensis both carry orthologues of the core T4P components — the assembly ATPase PilB, the retraction ATPase PilT, the inner-membrane platform protein PilC, the prepilin peptidase PilD, and the alignment-complex proteins PilM and PilF — together with the hormogonium-associated hmpD, hmpF and hmpG. Pseudanabaena sp. is therefore not pilus-deficient, and it does glide, yet it does not collect particles. The difference between the two organisms consequently cannot be attributed to the presence or absence of the pilus motor, which we would argue supports the physical argument we make here. Consistent with this, we have not identified mutations in pilus-related genes in the mutant, non-motile F. draycotensis.

      We are currently in the process of preparing another manuscript describing the mutations that led to motility loss in the non-motile F. draycotensis, as well as the proteins that are differentially expressed in the motile and non-motile F. draycotensis. These analyses will shed more light on the molecular mechanisms abolishing motility and how they might be influencing particle collection.

      In the revised manuscript, we will make these points clearer in the Discussion section.

      Reviewer #2 (Public review):

      Summary:

      The authors studied aggregation, buckling, and particle collection by the filamentous cyanobacterium Fluctiforma draycotensis, as well as by the filamentous Pseudanabaena sp. (order Pseudoanabenales). They performed a range of experiments, from imaging individual gliding filaments to multiple-day experiments showing the formation of large aggregates around a particle formed from a precipitate. They also developed a model of buckling filaments to argue that the ability of elastic filaments to collect particles and form macrostructures is confined to a part of the filament phase space in terms of length and flexibility, meaning that gliding combined with certain filament length and flexibility naturally reproduces the observations.

      Strengths:

      This is an impressive study that uses multiple tools to connect macrostructure formation with filaments’ gliding motility and buckling. It adds an important perspective on the biological and physical factors at play in the emergence of aggregates.

      We thank the reviewer for the accurate summary of our work and highlighting the strengths of the study.

      Weaknesses:

      The authors ignore the possibility that filament behavior plays an important role in the emergence of the observed patterns. Cyanobacteria have been shown to control their gliding motility (Pfreundt et al Science 2023; Kurjahn et al Nature Comm 2024), and their molecular motors are known to be regulated by chemotaxislike signaling pathways (Risser ARM 2025). As far as I know, how the coordination between the pulling agents along an individual filament works is actively debated, but there seems to be little doubt that it exists.

      To illustrate this point better, note that the aggregation observed by the authors is consistent with the length-dependent ability of filaments to coordinate gliding (I’m not saying this is how it works in Fluctiforma draycotensis; I’m saying it’s consistent). Suppose the coordination requires sufficiently long filaments, which could be the case when signaling molecules travel along the filament, propagating information about when individual pulling agents should reverse. In such a model, short filaments act randomly because they fail to coordinate gliding by the time they glide off nascent aggregates, whereas longer filaments can perform informed reversals because they have more time for coordination. Such behavior then explains the lack of aggregation in Pseudanabaena sp. (via behavior, not lack of stiffness). Note that Trichodesmium is stiff; its filaments do not buckle, yet Trichodesmium forms organized aggregates via tightly controlled motility. Note also that, as the authors report, since Pseudanabaena sp. is both shorter and faster, its filaments have relatively (to the time needed to glide the filaments’ length) little time to coordinate reversals. In my opinion, whether the observed patterns passively emerge from gliding and buckling or result from active behavior remains an open question.

      We appreciate the comment by the reviewer. We certainly agree that behavioral responses exist in filamentous cyanobacteria and will interplay with the physical aspects to produce exciting, complex dynamics. Besides the exemplar ideas that the reviewer provides, there can be many other scenarios involving behavioral responses, such as responses to light and to quorum sensing molecules or photosynthesis-generated radicals. For example, in F. draycotensis we have observed photo-responses at the aggregate level, which we are are currently studying. Photoresponses are also observed in Trichodesmium aggregates [7]. In general, a full understanding of the interaction of the biological (i.e. behavioral) and the physical aspects will require several future studies.

      In the current study, however, we focus on characterising the physical aspects of gliding motility alone, combined with experimental observations. We believe that this approach is important to establish a form of “null expectation” from the physics of gliding, elastic filaments alone. Currently, the molecular mechanisms responsible for coordinating the reversal behaviour of multiple filaments are still unclear, so it is difficult to experimentally demonstrate behavioural contributions to aggregate formation, e.g. via experiments where such behaviour is switched off. In the meantime, simulations such as those presented here allow us to test more precisely the potential role of activity, coordinated reversals and the elastic properties of the filament. In future it will be interesting to scale up the presented model to include multiple interacting filaments, and to systematically test the respective roles of active coordination behaviour for one individual filament (reversals) and for multiple interacting filaments (where contacts modulate activity), as well as the physical properties (length and flexibility). Such modelling studies can then identify if a ‘purely physical’ model can or cannot generate realistic aggregates, and pinpoint whether additional coordination mechanisms are needed to regulate aggregation. By testing the combination of different physical and biological coordination mechanisms, it would then help to indicate how much of a role is played by various potential active coordination behaviours.

      We will bring out this point more clearly in the Discussion section of the revised manuscript.

      I also have a small suggestion regarding this statement on model novelty:

      The essential novelty of this model is that the filament itself is active and out of equilibrium, and additionally, the forces and torques are applied locally along its centreline, and not at its extremities as in previous steady-state mechanical studies of elastic, twistable filaments such as DNA [31-33] (see Methods and SI).

      This statement needs to be revised as it ignores a substantial body of work on self-organization of active filaments: (R. E. Isele-Holder, J. Elgeti, G. Gompper, Soft Matter 2015; Pfreudnt et al, Science 2023; Faluweki et al PRL 2023; Kurjahn et al Nature Comm 2024).

      We agree with the reviewer that there is a significant literature on active filaments, some of which we have already cited and will now discuss in more details, as well as adding and discussing the suggested additional references. Our statement on “model novelty” refers to the analysis of buckling instabilities of biological filaments, and in particular DNA, due to a combination of forces and torques. To our knowledge, this has only be studied explicitely by [4], and only in the local (resistive force theory) limit. The elastohydrodynamic simulations coupled to local active forces and torques, as we implemented here, are therefore novel and will expand the analysis of both microbial filaments and other biological polymers. We will clarify these points in the Methods and Discussion sections of the revised manuscript.

      Last point: the authors often say that their observations are reproducible (’...reproducibly forms macroscopic granules...’). What is meant? Different experiments on different days, different aliquots?

      The “replicability” statement was in reference to different experiments started on different days using cultures obtained from serial transfer experiments, as well as cultures re-initiated from cyrostocks. This point will be made clear in the revised manuscript.

      Reviewer #3 (Public review):

      Summary:

      The authors report and characterize the formation of aggregate microstructures by the motile filamentous cyanobacterium Fluctiforma draycotensis, which exhibits gliding motility accompanied by rotation along the long axis while excreting EPS. In experiments with motile F. draycotensis cultures, they observed the formation of granular structures composed of cyanobacteria and other material (iron, polystyrene beads, etc.), with macrostructures on the scale of 1mm within 24 hours. The structures were motile at speeds comparable to that of the cyanobacteria filaments, resulting in their growth through coalescence over time. Notably, such macrostructures were absent in nonmotile F. draycotensis, pointing to the role of filament motility in their formation. Through experiments examining the micro-scale dynamics, inert material such as small polystyrene beads was found to be transported by the gliding, buckling, and plectoneme dynamics of the filaments, pointing to the underlying mechanism by which particles are collected into larger-scale microgranule structures.

      To interrogate the properties that drive the cyanobacteria filament buckling, plectoneme formation, and entanglement, the authors develop a mechanical model for filaments as nearly inextensible, slender bodies with resistance to twisting and bending under active gliding forces and torques and responding to fluid flows and surface adhesion. They derive expressions for the thresholds for buckling and twisting instabilities, which are additionally demonstrated and interrogated through simulation via the Immersed Boundary Method. Most importantly, bending and plectoneme formation only occur with sufficiently long filaments, and the threshold is shorter for bending than for plectoneme formation. Experimental observations with wild-type filaments agree with the model-predicted thresholds. The authors perform additional experiments with shorter filaments below both thresholds, including the filamentous bacterium Pseudanabaena, which fail to collect particles (though can in principle form macrostructures).

      Strengths:

      This work appears to be novel (notably, the discovery and characterization of the particle collection behavior of a filamentous cyanobacterium) and has interesting implications for both naturally observed cyanobacterial macrostructures as well as the controllable parameters in engineering them. The experimental and modeling work is well motivated, contributing to the broader understanding of macrostructure formation and material aggregation through active filament dynamics (not exclusive to cyanobacteria), as well as the underlying physical properties governing important filamentous cyanobacterium dynamics. As such, I would expect the results of this paper to be of broad interest to both biophysicists and microbiologists. Generally, the manuscript is well written with clear, compelling figures that illustrate the important conclusions of this study.

      We thank the reviewer for the accurate summary of our work and recognising the broad relevance of the study.

      Weaknesses:

      In the section on “Shorter gliding filaments cannot collect particles nor form granule macrostructures”, the filamentous cyanobacteria considered “all” fall below the predicted thresholds for bending and twisting. The “long” F. draycotensis are 60 microns in length, notably less than the 120 and 320 micron thresholds derived in the previous section as well as the lengths of filaments considered in Figure 3D, yet these “long” 60 micron filaments form macrostructures. How can this be understood in the context of the model predictions? Is the nature of the macrostructures in Figure 4B, the microscale parameters, or the collection of particles somehow different than those with filaments an order of magnitude longer in earlier parts of the paper? The paper would be stronger if these sorts of questions were addressed in the text and/or with supplementary figures.

      We thank the reviewer for this point. Indeed as we mention in the text, the ‘long’ population has a mean length of 60 micron. However, as shown in the length distribution plot in Fig 4A, the maximum filament lengths observed in these populations (within the samples used for microscopy) are 560 microns for the long filaments, versus 240 microns for the short filaments. Thus, we expect the long population to contain multiple filaments that can buckle and a few that can form plectonemes, whilst the short population might have some buckling filaments and none that form plectonemes. We stress that Fig 4A only shows the length distribution for what we believe to be a representative sample taken from the long and short populations, not the full data from the entire population.

      We will revise the main text to include the maximum filament lengths of the two populations as well as the mean values. We will also add lines to Fig 4A to indicate the buckling and plectoneme threshold lengths from the analytical estimate for the F. draycotensis filaments (same values as in Fig 3), to make it clear that the long population contains more buckling/plectoneming filaments than the short population.

      References:

      (1) E. S. Cameron et al. “Diversity and specificity of molecular functions in cyanobacterial symbionts”. In: Sci Rep 14.1 (2024), p. 18658. issn: 2045-2322 (Electronic) 2045-2322 (Linking). doi: 10.1038/s41598-024-69215-8. url: https: //www.ncbi.nlm.nih.gov/pubmed/39134591.

      (2) M. Chuvochina et al. “Proposal of names for 329 higher rank taxa defined in the Genome Taxonomy Database under two prokaryotic codes”. In: FEMS Microbiol Lett 370 (2023). issn: 1574-6968 (Electronic) 0378-1097 (Print) 0378-1097 (Linking). doi: 10.1093/femsle/fnad071. url: https://www.ncbi.nlm.nih.gov/ pubmed/37480240.

      (3) S. J. N. Duxbury et al. “Niche formation and metabolic interactions contribute to stable diversity in a spatially structured cyanobacterial community”. In: ISME J (2025). issn: 1751-7370 (Electronic) 1751-7362 (Linking). doi: 10.1093/ismejo/ wraf126. url: https://www.ncbi.nlm.nih.gov/pubmed/40577531.

      (4) Raymond E. Goldstein, Thomas R. Powers, and Chris H. Wiggins. “Viscous Nonlinear Dynamics of Twist and Writhe”. In: Physical Review Letters 80.23 (June 1998), pp. 5232–5235. issn: 1079-7114. doi: 10.1103/physrevlett.80.5232.

      (5) D. H. Parks et al. “A standardized bacterial taxonomy based on genome phylogeny substantially revises the tree of life”. In: Nat Biotechnol 36.10 (2018), pp. 996–1004. issn: 1546-1696 (Electronic) 1087-0156 (Linking). doi: 10.1038/ nbt.4229. url: https://www.ncbi.nlm.nih.gov/pubmed/30148503.

      (6) D. H. Parks et al. “GTDB release 10: a complete and systematic taxonomy for 715 230 bacterial and 17 245 archaeal genomes”. In: Nucleic Acids Res 54.D1 (2026), pp. D743–D754. issn: 1362-4962 (Electronic) 0305-1048 (Print) 03051048 (Linking). doi: 10.1093/nar/gkaf1040. url: https://www.ncbi.nlm.nih. gov/pubmed/41123020.

      (7) U. Pfreundt et al. “Controlled motility in the cyanobacterium Trichodesmium regulates aggregate architecture”. In: Science 380.6647 (2023), pp. 830–835. issn: 1095-9203 (Electronic) 0036-8075 (Linking). doi: 10.1126/science.adf2753.

      (8) Douglas D Risser. “Motility in Filamentous Cyanobacteria”. In: Annual Review of Microbiology 79 (2025).

      (9) Jerko Rosko et al. “Cellular coordination underpins rapid reversals in gliding filamentous cyanobacteria and its loss results in plectonemes”. In: eLife 13 (2025), RP100768.

      (10) A. Scarampi et al. “Enrichment of convergent metabolic functions in microbial communities through imposed and emergent environmental niches”. In: bioRxiv (2026). doi: 10.64898/2026.02.11.705344.

    1. the lack of social validationsuggests a failure to achieve status as a social peer

      why does lack of active engagement entail grok’s failure to achieve status as a social peer, esp given that passive engagement is at par?

    2. objective judge of credibility and competence.

      isn’t them percieving that as an objective “judge” contingent on how they uptake the response it produces, as opposed to just invoking it for such questions? an alternative explaination could be that they are looking for confirmation, which does not entail that it is a “judge” with immense credibility, but just another actor with baseline cognitive sophistication that is one of many supporters of their positions. in that case they are considering grok to be as much of an “objective judge with credibilty and competence” as much as me asking a random person on the street their opinion when im disagreeing with a friend to bolster my point somewhat.

    3. Large Language Models (LLMs) are increas-ingly deployed as active participants on pub-lic social media platforms,

      what are some prominent examples?

    4. ts responses have the potentialto shape user opinions, validate existing beliefs, orintroduce new frames into sensitive debates

      why do we not test for this using our dataset.

    5. Yetthis very neutrality becomes a high-stakes featurewhen applied to questions with no neutral answe

      hmm interesting characterisation — “no neutral answer”

    1. Whole exome sequencing identifies a novel splice-site mutation in IMPG2gene causing Stargardt-like juvenile macular dystrophy in a northIndian family

      PMID:35973334

      Gene: ABCA4

      HGNC ID: 34

      Case#: eldest sister II.2 aged 22 year

      Variant splice-site variant NC_000003.11(NM_016247.3):c.1239 + 1G > T [Chr3:100972539C > A

      FammilyInfo two-generation north Indian family with three members affected with Stargardt-like macular dys trophy

      CasePresentingHPOs:ow vision and difficulty in night vision, with symptoms starting in the early second decade of life, which progressed slowly over time

      PedrigreeIn the results section is mentioned

      CaseHPOFreeText:NA

      CaseNotHPOs:Na

      CaseNotHPOFreeText:NA

      Genotyping Method:2.3. Validation of identified variant by Sanger sequencing

      PreviouslyPublished:NA

    1. Case 3

      Case#:3, 34–year-old woman

      DiseaseAssertion:NR

      FamilyInfo: no family history of an ocular disease

      CasePresentingHPOs: HP:0007663, HP:0030506, HP:0030528, HP:0000603

      CaseHPOFreeText:bilateral markedly decreased vision (logMar BCVA OD:0.93, OS:0.95), bilateral atrophic lesions of the macula accompanied by yellow-white stellate flecks at the level of the retinal pigment epithelium, Atrophic lesions and flecks were also extending to the mid-periphery of both retinae, bilateral absolute central scotoma and relative paracentral scotomas as well in both eyes, PERG was significantly reduced, while scotopic and photopic amplitudes were also lower than normal.

      CaseNotHPOs:

      CaseNotHPOFreeText:NR

      Genotyping Method: “Analysis of the ABCA4 gene”

      PreviouslyPublished: No

      Variant: NM_000350.3:c.5882G>A, NM_000350.3:c.6709dup

      ClinVar: ​​7888, 99485

      CAID: CA119132, CA227437

      SupplementalData:

    1. Interestingly, we identified one 30-year-old patient (ARDM-247), double heterozygous for the p.Arg1129Leu and p.Cys2137Tyr alleles, who presented a CRD phenotype. This p.Cys2137Tyr change was located more towards the amino terminus. Moreover, in other study, the results showed that the changes located in this zone appear to result in altered processing of the protein and to be associated with an earlier onset of disease.16 The p.Cys2137Tyr change in combination with the p.Arg1129Leu allele produced a CRD phenotype. Therefore, we speculate that the novel p.Cys2137Tyr variant could be a severe allele which is modifying the patient’s phenotype.

      Case#: Family MD-0247/ARDM-247 Proband, 12yo at onset

      DiseaseAssertion: AR cone rod dystrophy

      FamilyInfo:

      CasePresentingHPOs:

      CaseHPOFreeText: STGD diagnosed based on "bilateral central vision loss; fundus presenting with a beaten-bronze appearance and/or the presence of orange-yellow flecks in the retina from the posterior pole to the mid-periphery; fluorescein angiography showing typical dark choroid; and normal to subnormal electroretinogram (ERGs)." VA loss, VF loss, BCVA=0.05/0.1, cone-pattern on ERG

      CaseNotHPOs:

      CaseNotHPOFreeText:

      PreviouslyPublished: n/a

      Variant: c.6410G>A (p.Cys2137Tyr); c.3386G>T (p.Arg1129Leu) found by ABCR400 + dHPLC + HRM

      ClinVar: 2202779

      CAID: CA341277358

      SupplementalData: n/a

    1. 18 5709 Mi 9 2/10 / 2/10 c.32T>C(1) / c.1804C<T(13) p.Leu11Pro/p.Arg602Thr

      Case#: Maia-Lopes Family 18 Proband 5709, Portuguese, 9yo at onset

      DiseaseAssertion: STGD

      FamilyInfo: Family 18

      CasePresentingHPOs:

      CaseHPOFreeText: mild central fundus changes, vision: 2/10 / 2/10

      CaseNotHPOs:

      CaseNotHPOFreeText:

      GenotypingMethod: ABCR400 microarray, dHPLC

      PreviouslyPublished: n/a

      Variant: c.1804C<T/ c.32T>C(1); p.Arg602Thr/p.Leu11Pro

      ClinVar: 99217

      CAID: CA227106

      SupplementalData: n/a

    1. ABCA4

      In supplemental table S2, Case LL291 is a female listed as "likely solved" with a clinical diagnosis of CRD. Proband is compound heterozygous for c.32T>C p.(Leu11Pro) and c.5196+1137G>A.

      In Supplemental table S1, this proband is a female, 54yo at report, 50yo at onset. Nyctalopia, visual acuity=0.05/0.6, OD: extremely high myopia OS: high myopia, no family information

    1. Patient 2 (P2)

      Case#: Somali origin, 11 years of age, onset age 8

      DiseaseAssertion: STGD

      FamilyInfo: Parents were heterozygous for Arg212Cys, asymptoamtic 32 year old father was found to be homozygous for Gly1961Glu

      CasePresentingHPOs: HP:0011504, ORPHA:827, HP:0000608

      CaseHPOFreeText: BCVA 20/70 and 20/80, atrophic zone of outer retinal and RPE atrophy, Bull's Eye Maculopathy, Macular atrophy

      CaseNotHPOs: N/a

      CaseNotHPOFreeText: N/a

      Genotyping Method: Whole genome sequencing

      PreviouslyPublished: N/a

      Variant: NM_000350.3:c.5882G>A p.(Gly1961Glu) ; NM_000350.3:c.634C>T p.(Arg212Cys)

      ClinVar: 7888; 7898

      CAID: CA119132, CA203216

      SupplementalData: N/a

    1. In humans, clinical studies have implicated mutations in 19 of the 48 known ABC transporters in diseases such as cystic fibrosis and adrenoleukodystrophy.

      Annotating here since the article is a PDF.

      This variant is mentioned in table 2 as a "disease associated mutation", but only as being present in the NBD/NBD interface. No further details are provided.

    1. Patients older than 60 years or with ocular comorbidities such as diabetic retinopathy, uveitis, or glaucoma were excluded. From the remaining list, subjects for whom high-resolution SD-OCT imaging was available were selected. A review of the patient imaging and medical records was performed to identify those who received a clinical diagnosis of Stargardt macular dystrophy based on their clinical phenotype, including color fundus, infrared, FAF, and fluorescein angiography images and electroretinographic findings. 12

      Case#: P3, male, 16yo at report, 9yo at dx, US with Indian ethnicity

      DiseaseAssertion: Stargardt

      FamilyInfo: n/a

      CasePresentingHPOs:

      CaseHPOFreeText: BCVA (logMAR)= OD=20/160 (0.90), OS=20/125 (0.80)

      CaseNotHPOs:

      CaseNotHPOFreeText: ocular comorbidities such as diabetic retinopathy, uveitis, or glaucoma

      GenotypingMethod: "genetic testing"

      PreviouslyPublished: n/a

      Variant: c.2453G>A; c.4532C>A (p.Pro1511His)

      ClinVar: 99135

      CAID: CA227000

      SupplementalData: table 1

    1. Finally, we examined whether the phenotype‐associated known/candidate pathogenic variants could explain the patient's disease, andif the MAF in population‐matched control data (8.3kJPN) was relatedto disease prevalence. Patients were classified as “Solved” if theirgenotype was consistent with their clinical phenotype. Patients wereclassified as “Partially solved” when a heterozygous known/candidatepathogenic variant was detected in a recessive allele, but without anadditional variant in trans. Patients were categorized as “Unsolved” iftheir genotypes exhibited either no candidate pathogenic variants ormultiple heterozygous pathogenic variants that did not explain thephenotype clearly. Variants annotated as causal for solved patients arelisted in Supporting Information: Table S2. Novel variants identified inthis study are listed in the second sheet of Table S2. SupportingInformation: Table S3 shows the phenotypes and genotypes of solvedpatients.

      This variant is listed in supplementary tables S2 and S3. Proband KN-187 is a "solved" patient, meaning the phenotype matches the genotype. Compound heterozygous (c.6290C>T p.P2097L; c.6445C>T p.R2149X) male with Stargardt disease- all that is provided.

    1. STGD1 was determined according to initial symptoms of VA loss; fundus images showing orange-yellow flecks in the retina, a beaten-bronze appearance; and normal or cone-altered ffERG results

      Case#: MD-0790, Spanish

      DiseaseAssertion: STGD1

      FamilyInfo: n/a

      CasePresentingHPOs:

      CaseHPOFreeText: "STGD1 was determined according to initial symptoms of VA loss; fundus images showing orange-yellow flecks in the retina, a beaten-bronze appearance; and normal or cone-altered ffERG results"

      CaseNotHPOs:

      CaseNotHPOFreeText:

      GenotypingMethod: Index cases were studied by different next-generation sequencing (NGS) strategies, including targeted gene panels, clinical exome, and/or whole-exome sequencing

      PreviouslyPublished: n/a

      Variant: c.1715G>C p.(Arg572Pro); c.4918C>T p.(Arg1640Trp)

      ClinVar: 99073

      CAID: CA226919

      SupplementalData: Table S1

    1. that theater is not “about” dialogue after all, but what might be called fields of dramatic tension; force-fields of human relationships beneath the level of language

      Q: Oates says theater isn't really "about" the dialogue, the real meaning is something underneath it. Does she mean subtext or something else entirely?

    2. What shimmers with life on the page may die within minutes in the theater precisely because prose is a language to be spoken to an individual, recreated in an individual reader’s consciousness, usually in solitude, while dramatic dialogue is a special language spoken by living actors to one another, a collective audience overhearing

      O: While theater may be one of a few places that can simply "display" the character without direction or added context. descriptions, or remarks, written prose that can seem rather complex and alive also exist, if the character is demonstrated through actions, background or incentives, internal thinking or self-muses, which brings out a more vivid character, whose actions are realistically dependent on previous influences, instead of being a "character sheet" character or Mary Sue, or of that who follows a classical trope. The character is not created as a plot driven device, and acts with intention, but also unconscious decisions, much like what happens in real life, which illustrates and adds "noise" that creates more complex and compelling characters, continuous instead of something discrete, free from defined tropes, vast and liberating. This type of fiction is often called literary fiction.

    1. Adieu, adieu, adieu.⌜

      Meant to submit this annotation on thursday's class, completely forgot to submit

      This entire scene and before is handled really well in the hoffman version, and compared to the BBC version it's night and day. I enjoyed the shift between humor and seriousness in Flute's acting in the Hoffman version, and didn't like theseus' death distracting from that in the BBC version.

    1. фигурка представляет собой прямоугольник шириной в W i W_i Wi​ клеток и высотой в одну клетку; самая левая клетка фигурки имеет абсциссу a i a_i ai

      Мы никак не можем шевелить эти плитки.

      только порядок можем выбрать. никакие манипуляций не разрешены с ними. ни перемещение не перевороты.

      нужно сделать получившуюся фигуру минимальной высоты

    1. Figure 3. The same five configurations measured on one workload, so the levers read as numbers rather than directions. Scaling replicas dp=1 to dp=8 - 8 GPUs to 64 - multiplies decode throughput by 5.7, not 8. On a fixed 8 GPUs, changing only the parallelism to tp=4/pp=2 gives 1.85×. TTFT here is dominated by queueing at 5,000 arrivals, so read the fall from ~956 s to ~117 s as backlog draining, not as model latency.

      same Q as before.. 16 > 32 > 64 - wy would throughput not increase in the same proportiton? you are effetively adding replicas

    2. The knobs that decide it Each knob trades one outcome against another. The direction, the mechanism, and where the direction stops holding: Tensor parallelism, tensor_parallel_size. Each GPU holds a shard of every layer, so weight bytes, KV-cache bytes and arithmetic per GPU per token all fall roughly as 1/TP - and because decode is bandwidth-bound, less to stream means a shorter step. That benefit has a ceiling. Two all-reduces run per layer, one after attention and one after the MLP; the bytes each GPU sends level off as TP grows, but the number of synchronisation steps does not, and the cost jumps whenever the group stops fitting inside one fast interconnect domain. KV sharding also stops at the model’s KV-head count. Meanwhile total GPUs = tensor parallelism × pipeline parallelism × replicas, so at a fixed GPU count raising TP costs you replicas. Figure 3 shows TP=8 sitting past the useful point for this model: against TP=4 with PP=2 on the same 8 GPUs it delivers 46% less throughput, at 48 ms p99 TPOT against 47 ms - no latency gain at all. That comparison moves PP as well, so it bounds the benefit of higher TP rather than isolating it - but it is enough to say there is no single direction here. TP has a useful range, and 8 is above it. Replicas, data_parallel_size. Each replica is a full copy of the model with its own GPUs, so aggregate throughput rises with replica count - sublinearly, and only while arrivals keep every replica batched. Nothing pools across replicas: weights are duplicated per replica, the prefix cache is per-replica so a 70% hit rate inside one replica does not become 70% across eight, and a request routed to a busy replica cannot borrow a quiet one’s free KV blocks. In Figure 3, dp=1 to dp=8 gives 5.7×, not 8×. Every configuration in that sweep is backlogged - average TTFT is 117 seconds even at dp=8 - so it is a saturated-throughput comparison rather than a latency one, and the gap against linear scaling is not attributable to any single one of those causes from this data alone. Batch size, capped by the token budget, by max_num_seqs, and by free KV blocks - whichever binds first. Larger batches raise decode throughput because the weight bytes streamed per iteration are amortised across more tokens. While weight bytes dominate, step time is nearly flat as the batch grows, so that throughput is close to free. It stops being free in two ways. Once the live KV cache outweighs the weights - the reasoning case above, at 1.5 GB per request - every added request adds bytes to stream and step time grows with batch size directly. And once kernels turn compute-bound, step time grows with arithmetic instead. Either way TPOT is what pays. Net: throughput up, TPOT up, cheap only in the weight-dominated regime. Precision and KV-cache format. FP8 weights instead of BF16 halve weight bytes, which frees HBM for KV cache and cuts what decode must stream, relieving capacity and bandwidth together. An FP8 KV cache does the same to the cache itself. What it costs is output quality, and not equally: quantising the KV cache degrades more than quantising weights at the same width, because the cache is read back as attention state for the rest of the request. Kernel support also varies by GPU, so the gain does not carry across the hardware axis of the grid.

      reduce the size of this section too. by 70%.

    3. very request runs prefill once, then decode in a loop. Each panel shows what those phases do to one resource, drawn to scale. Coding fills the token budget: prefill is 96% of the work, and a prompt too large for one pass is split into chunks that then compete with decode. Chat fills the latency budget: 48 ms of a 50 ms target, so the batch stops growing while HBM sits largely unused. Reasoning fills HBM, and the two mirrored bars are the mechanism - cache per request grows 6× as the request runs, so requests that fit fall 6×, from 295 to 49, while they are still running.

      keep the terminology consistent.. latency bound vs kv cache bound is not in the same group. keep it flops vs vram vs b/w

    4. A coding assistant sends huge prompts - the whole file, the repo context - and gets short edits back. At ISL ~8,000 and OSL ~300, prefill is 96% of the token work. What binds is the token budget: a prompt larger than one forward pass can hold is admitted in chunks, one chunk per iteration - 2 chunks at a 6,144-token budget, 12 at 682. Each of those iterations spends most of its budget on prefill and has little left for decode, so decode tokens per iteration drop and aggregate decode throughput falls for the whole batch while that prompt is being read. Every request already decoding waits longer for its next token, so TPOT rises for all of them. The chunked request needs 12 iterations before its first output token instead of 2, so its TTFT rises with the chunk count. A chatbot sends short prompts and needs a steady stream of words back. Prefill is over quickly, so the request is almost entirely decode, and one iteration advances every active request by exactly one token - so TPOT is iteration time. Admitting one more request lengthens that iteration, because the GPU must compute one more token and stream one more request’s KV cache through HBM within the same forward pass. With a 50 ms per-token SLO against the ~48 ms p99 TPOT measured for this model on an H100 node, 2 ms of margin remain: the next admission consumes it and the one after breaches the SLO. So the SLO caps batch size, and because throughput scales with batch size it caps tokens per second per GPU. Serving the same traffic then takes more replicas - you buy GPUs to hold a latency line rather than to do work. Memory is not the constraint here: each request holds ~125 MB of KV cache, so 72 GB would hold roughly 590 concurrent requests, far more than the latency target allows. A reasoning or agent workload writes far more than it reads - long chains of thought, thousands of output tokens. Every one appends to the KV cache, so a request that begins with 2,000 tokens of cache ends holding 12,000: 250 MB becomes 1.5 GB. Capacity: the same 72 GB holds 295 of these requests at the start but 49 once grown, so as the live requests age fewer new ones fit, the batch drains faster than it refills, and aggregate throughput decays over the run. Bandwidth: decode re-reads the entire live cache every iteration, so bytes streamed per iteration grow as the caches grow and TPOT rises as requests age.

      this section has too much detail. try to keep it simple with fewer numbers and state the constraint clearly towards the end - flops or vram or b/w.

    1. Give me that boy and I will go with thee.

      In the BBC version, this was changed into "Give me your hand". No kid is mentioned, but still Titania gets mad, refuses. This adaptation changes a lot of things: we go from a dispute about an exterior factor that they disagree on, to a couple fighting. Oberon is a humiliated husband in this scene, which makes it way more personal.

    2. Do you amend it, then. It lies in you. 0486  Why should Titania cross her Oberon?

      In Davies version, Oberon, speechless, cannot speak. Pucks intervenes and take these lines, which shows his importance to his master

    3. At a fair vestal thronèd by ⌜the⌝ west, 0529 165 And loosed his love-shaft smartly from his bow 0530  As it should pierce a hundred thousand hearts. 0531  But I might see young Cupid’s fiery shaft 0532  Quenched in the chaste beams of the wat’ry moon, 0533  And the imperial vot’ress passèd on

      Elizabeth I. The compliment is that Cupid's arrow missed her and she stayed unmarried by choice. The arrow lands on a flower instead and causes everything that goes wrong for the next three acts.

    4. Well, go thy way. Thou shalt not from this grove 0516  Till I torment thee for this injury.—

      The BBC version stands out for many changs, and this one is one too: Oberon threatens Titania even though she is still here, compared to after she leaves in the original version

    5. If you will patiently dance in our round 0510  And see our moonlight revels, go with us. 0511  If not, shun me, and I will spare your haunts.

      In Bennett's version, Titania keeps the invitational, but in the BBC one, same thing: it is cut. Its now Oberon's turn to try and make the move

    6. Set your heart at rest: 0490  The Fairyland buys not the child of me. 0491  His mother was a vot’ress of my order, 0492  And in the spicèd Indian air by night 0493  Full often hath she gossiped by my side 0494 130 And sat with me on Neptune’s yellow sands, 0495  Marking th’ embarkèd traders on the flood, 0496  When we have laughed to see the sails conceive 0497  And grow big-bellied with the wanton wind; 0498  Which she, with pretty and with swimming gait, 0499 135 Following (her womb then rich with my young 0500  squire), 0501  Would imitate and sail upon the land 0502  To fetch me trifles and return again, 0503  As from a voyage, rich with merchandise. 0504 140 But she, being mortal, of that boy did die, 0505  And for her sake do I rear up her boy, 0506  And for her sake I will not part with him

      This passage is completely cut in the BBC version, because the kid is not mentioned. This is the main reason why the fight in this version is more personal: no exterior element is included

    1. Listing tasks may not be all that you need to do. There may be so many tasks that you must group them so that readers can find individual ones more easily.

      I enjoyed this point as it made me aware of how to properly convey my voice when in situation where I may need to explain multiple things.

    1. Overall, highlighting specific results in each résumé category will increase your chances of getting your résumé noticed

      I think this is an important and true point. Having specific examples of experience tailored to individual sections can show your strength across multiple areas.

    1. Quoting the work of others in your writing is fine, provided that you credit the source fully enough that your readers can find it on their own. If you fail to take careful notes, or the sentence is present in your writing but later fails to get accurate attribution, it can have a negative impact on you and your organization.

      This is a good summarizing point and also important to consider who you're doing work for.

    1. When you begin picking team members for a writing project in a technical writing course, you should consider choosing people with different backgrounds and interests. Just as a diverse, well-rounded background for an individual writer is an advantage, a group of diverse individuals makes for a well-rounded writing team.

      I agree with this point and think it is very important to consider when forming a team for a writing project.

    1. “Since Berniehas taken a salary from taxpayers for 40 years, he would know nothing about this.”

      I think this is an example of mudslinging or ad hominem because they are attacking Bernie Sanders personally instead of just talking about whether the 32-hour workweek would work. They bring up how long he has been getting a salary from taxpayers and use that to question his knowledge about the issue. This takes attention away from the actual argument.

    2. Get ready for a long weekend at Bernie’s!

      I think this is an example of sensationalism because the article starts with “Get ready for a long weekend at Bernie’s!” instead of just introducing the bill. The phrase is meant to catch the reader’s attention and make the story sound more entertaining. It also connects Bernie Sanders to the idea of a long weekend before the reader even learns the details of the proposal.

    3. Ben Shapiro believes Americans should work until they drop dead,and Bernie Sanders believes Americans deserve more time off.

      I think this is an example of word choice because the author says Ben Shapiro believes Americans should “work until they drop dead.” That wording makes Ben Shapiro’s position sound much more extreme then simply saying he supports raising the retirement age. The author then describes Bernie Sanders as believing Americans “deserve more time off,” which gives Sanders a more positive description. The way its wording makes one seem bad and one good which can influence bias.

    4. It’s a dismal situation.

      I think that “dismal” shows subjective bias because it makes the situation sound very negative. The author could have just explained that workers have less free time and more unpredictable schedules, but instead they call it a “dismal situation.” This makes the reader see the situation in a negative way before they can form their own opinion.

    1. Frontier AI companies within democratic countries coordinate to establish common safety standards as well as limits on the rate of unchecked AI progress.

      Utilitarian Lens: The purpose of safety regulations in general is to efficiently manage and prevent hazards and/or risks, and in turn, any resulting incidents. The purpose of setting restrictions in this scenario is to more easily keep track of and manage progress, and the frequency at which it is occurring. Preventing any significant losses (for instance, inefficient time and/or resource usage), that could potentially occur due to a lack of awareness. In the process of decision-making, regulations can also be beneficial in determining which decisions would result in the best possible outcome for a business.

    2. More generally, society must have a say in how this technology is used, and more time for the necessary public deliberations — which pacing the frontier would bring us — is surely a good thing.

      Common Good Lens: Assuming that Amodei is referring to the general public, the idea of reliance on "society" (Amodei) in (at least, partially) determining the extent of AI may be somewhat complicated due to the inevitability of divisive responses. Though alternatively, this could be in reference to the idea of constructive feedback and/or criticism, in that responses from the general public will contribute to the decisions regarding and improvement of AI in the near future. Regardless, it is most likely the businesses responsible for creating and/or holding ownership of AI technologies will ultimately make the final decision, which could be based on feedback, and/or a variety of other factors (legality, ethics and values, etc.).

    3. Thus, a key part of pacing within democracies is to keep democracies’ AI lead over autocracies as large as possible, to give us the breathing room we need in order to pace effectively.

      Justice Lens: The implications in this part of the essay focus on keeping the US, the "good guys", ahead of China, as the "bad guys". The argument here is right versus wrong based on the idea that we need to punish or restrict the "bad guy", based on the US view of the CCP.

    4. Any cooperation we are able to achieve with China will extend the amount of time we have to spend on pacing the frontier within the democratic nations. We should aim for the higher levels while seeing the lower levels as much more likely and realistic.

      Utilitarian Lens: I find it interesting that so much of the focus is on China as the most likely to match the US in the race to build the biggest and best AI. It seems short-sighted to ignore that there are other authoritarian powers that would be willing to do as much as China but not so much in the public eye. I appreciate that his focus is on the biggest, most obvious "enemy", but he is also simplifying his argument a lot. This argument does put most of the focus on the highest, fastest "common good", and would be the least expensive way to do it as well, in a highly utilitarian way.

    5. I believe that AI could cure most major diseases in the next 5–10 years, greatly accelerate economic growth rates, create a world of abundance and empowerment, and usher in a renaissance of democracy and freedom

      While all of this sounds good especially when looking at it through a Common Good lens. It doesn't really consider any of the consequences. Curing diseases is great and something that society has been aiming towards for a long time. We have made a lot of progress when it comes to medicine but there are still diseases that are uncurable. If we can speed up the rate in which medicine is developed more people can be saved. If most diseases are cured, then more people will also live. When looking at the same phrase with a Utilitarian lens things look a bit different. AI has consequences it greatly affects the environment by needing a lot of water to cool systems. With less water, plants and animals suffer. A lot of disease treatments have been found within plants and nature. While AI might be able to quickly figure out what is needed to cure a disease, we may not have the recourses to make the cure because they were put towards AI. If the environment is significantly destroyed to further AI development then we can't even use the information that AI gives us (curing diseases). There is also the idea of living rather than surviving. AI could help cure diseases and allow us to survive longer but if the world is seemingly destroyed is there really any point? Do the pros outweigh the cons? If AI is handled carefully and more research is done to be able to develop it in a more sustainable way, and the pros outweigh the cons, then AI could be beneficial.

    6. In other words, to create a race to the top. We have always devoted a substantial fraction of our efforts to studying, addressing, and informing the public about these AI risks, as well as advocating for well-considered regulation of AI, even when this gets us accused of hype, “doomerism”, or regulatory capture.

      Utilitarian Lens: Dario Amodei shows us how he thinks our world with AI is a race to the top and how we must pace the frontier. I think that AI is being used by a lot of companies and is replacing jobs that humans had before. We need to think about how this is affecting the population and the jobs being lost. Slowing down AI and regulating it is important before there is more harm or damage. AI is growing and evolving fast, so continuing to prioritize caution over speed and prudence over product is vital. While AI can have the benefits of advanced technology, regulating the use of AI before it grows too quickly is important.

    7. Carefully wielded, AI can be the latest in a long line of technological miracles that have uplifted and ennobled humanity.

      Applying the Utilitarian Lens, pacing AI serves as both a detriment and a positive for us. AI is at a weird stage where everyone in the world is competing to develop a more advanced system; however, the rate at which things are moving raises an eyebrow regarding public safety and privacy concerns. Slowing down the process of developing AI seems to be the obvious response when it comes to taking everyone into consideration—because the attention is literally going off of AI and onto the people. However, the reason it would be nice to speed up AI, for the population's sake, would be to help solve world issues such as disease. I don't know if we're quite at that because if we're weighing it out, the developing/understanding ratio is skewed; in other words, can we progress enough with AI at this speed in order to truly stay intellectually on par with it? If it becomes too advanced, would we know what to do? Applying a different lens like the Virtue, we arrive at a very similar conclusion. Honesty is something that I think most people can value in other people. If it's clear that AI is developing to the nth degree, why does the general public not know more about it? If these organizations were more candid about how little they truly understand as far as safety is concerned, then maybe we would take a step back. That would be the responsible thing to do, which is another trait people value.

    8. Today, however, the picture is totally different. The current models are an almost endless gold mine of insight into both how to build AI well and what can sometimes go wrong with it if it isn’t built well. I believe that if slowing down bought us even an extra year or two before models reach critical levels of capability, and we used that time to advance alignment, we could greatly reduce the risk that something goes seriously wrong.

      When applying the Utilitarian lens, Amodei points out how it would be beneficial to slow the production of AI so that we can better align with it so that there are safety protocols in place. This way would cause the least harmful affect because if AI companies were to slow down to better understand all aspects of what they are creating, they can better ensure that they will not get any unwanted or unfriendly outcome from their AI creations. This helps to keep everyone safe while still advancing the AI field, albeit a little slower than it is currently being developed. With that, while it may be beneficial to society in order to ensure its safety, it does come with the pitfall that it will take longer to find solution to different problems that they are using AI for which could cause harm along the way.

    9. I myself survived an early-stage cancer that would not have been treatable even fifty years ago

      Utilitarian Lens: One of the potential benefits that we may lose in pacing or delaying AI development is the creation of new drugs to prevent disease. Research shows it takes on average 10-15 years to develop a new drug that can be ready for patient use, and by slowing the development of AI, we could be slowing the development of these drugs when it already takes so long in the first place. On the contrary, slowing the growth of AI could be beneficial to saving the enviroment. The constant creation of new AI technology needs large data centers to help power it, and these data centers take tons of water and invade the forests and parks that were once there. If we slow the progress of AI, we could potentially cut back of the number of data centers needed until we can develop more laws and regulations for them. Before I can make a conclusive statement that slowing development causes more good than harm, I'd have to do research into the major fields that benefit from AI technology like medicine, and ask professionals in those fields on their opinion of slowing the progress of the AI they use. I highlighted what I believe is an assumption Amodei is making that also affects my judgement aswell. He actually should have a bias towards speeding up the development of AI because of its help in curing him of an early stage of cancer, yet he still takes the opposite stance. It seems he's not letting his personal bias get in the way of what he believes to be right. Hearing how AI has helped save lives does affect how I feel about this specific branch of its use as seen in my earlier comment, yet like Amodei, I'm trying to look at things objectively. If you apply the Care Ethics Lens to this same argument, you have to be wary of the patients this decision may affect if you slow AI's development, as it does have a positive use in the medical field.

    1. Columbus embarked on his historic voyage that failed to discover the route to Asia but instead discovered the New World

      Why did he fail to discover resources from asia, did he not find any resources or any other people that lived their?

    1. Bodily outlines are also dynamic within individuals, with whatindicates sex within-bodily outline changing over time, need, andage. Sex changes for everyone; for example, no one is born withbreasts. Instead, adults and younger people who have them do sobecause of endogenous hormonal processes, exogenous hormoneuse, surgical processes, and/or material in clothing items like bras.

      I found this interesting because it makes a good point. In America, people tend to say the female sex qualifications involve the identification of breasts. However, young adolescents do not have breasts until puberty, and even then, breasts don't always form for all females. Does that make them less female for lacking a certain quality that ties them to their sex?

    2. . Yet binaries are notoriously inept atdescribing the world, most poignantly because they exclude real-ities that exist beyond them (Fausto-Sterling, 2020; Waters, 2001).

      This part of the passage is almost, if not exactly criticizing the viewpoint of this concept. It states the fact that the sex-gender binary and nature-nuture/culture binary isn't considering the outside viewpoints of what sex and gender are across different cultures and enviorments. That they aren't considering that there might not be a definite universal way of viewing gender.

    3. Human biobehavioral research is often rooted in and supports two interrelated binaries: one that divides andessentializes humans into two sexes and values maleness over femaleness, and a second that divides behaviorinto two sides, assigning more truth to one—nature and sex—over the other—gender, nurture, and culture.

      Though I'm not sure this counts as a valid bias annotation; I can identify that two biases are being addressed. It is talking about two concepts of sex and gender. One group believes in a more conservative way of viewing the topic, in which there are only two possible sexes and that male qualities are valued over females. While the other view explores a more open-viewed way of seeing the concept of sex and gender across different variations.

    1. The No. 1 reason why crude oil prices and product prices have not been higher than they are is that China has cut its crude oil imports

      This is such a key point and speaks to China adaptability. In a time when crud oil prices spike. China is able to cut consumption of jet fuel. How was China able to cut it's use and is this something the United States can do.

    1. history is living and relevant to current concerns, not the “dead hand of the past,”

      The author is arguing that studying history helps us understand problems we face today. New events and new research can change what questions we ask about the past.

    2. humans are changing, overwhelming, or displacing other global processes of nature

      This shows us how environmental history has become more relevant over time. The author is stating that humans do not exist separately from nature, but that our actions can affect entire ecosystems.

    3. “The Rise of the West.” The new scholarship on Asia—which Jack Goldstone dubbed “The California School” because so many of us lived, worked, or published in California—raised questions about how and why the modern world came to have its essential characteristics: politically organized into nation-states and economically centered around industrial capitalism. Our findings that Asian societies had many of the characteristics others had seen as exclusively European and thus “causes” of the “European miracle” led us to argue that similarities cannot cause differences and to look for alternative explanations for how and why the world came to be the way it is. Andre Gunder Frank and Kenneth Pomeranz pulled this scholarship into two important books that changed the way we now understand how the world works, decentering Eurocentric explanations of history.

      This is important because it challenges the typical way history is taught that Europe is the main force behind global development.

    4. The past may not change, but our history of it does.

      This stands out to me because it explains that history is not just a big collection of facts together, history is constantly being reexamined rather than one simple answer.

    5. A growing human population requires additional food and energy supplies to support it, and given the agricultural technology available in 1400, those increases could come from but three sources: bringing more land under cultivation, increasing the labor inputs on a given plot of land (including selecting better seed), or increasing the amount of water or fertilizer.

      This causes some problems as it becomes naturally harder to support more and more people.

    6. Moreover, those densely populated regions of Earth corresponded to just fifteen highly developed civilizations, the most notable being (from east to west) Japan, Korea, China, Indonesia, Indochina, India, Islamic West Asia, Europe (both Mediterranean and West), Aztec, and Inca.

      People mainly lived in numbers back then, as it was not wise to live without community and not close to people.

    7. As an environmental historian myself, I found that observation compelling and incorporated an ecological theme in my narrative.

      This shows the authors credentials by stating that he is a primarily environmental historian.

    8. Later, during the Enlightenment of the seventeenth and eighteenth centuries, they attributed their superiority to a Greek heritage of secular, rationalistic, and scientific thought.

      I find it interesting that they changed their religion as time passed.

    9. discovery

      The authors mention of Columbus made me think about how much of a bad guy was. Many historical events show us the true colors of people and it makes me think: were there any good people? I suppose history focuses on the big, rich figures, and rich people can often be greedy. There were probably great men who were lost to time.

    10. world peace would prevail

      Today this is a crazy thought to have because there is controversy in almost everything. I guess back then world peace didn’t seem unattainable. It was a simpler time.

    1. eLife Assessment

      This study presents valuable findings regarding the role of the P2X7 receptor in retinoic acid signaling-mediated retinal remodeling following photoreceptor degeneration. However, the evidence supporting the main conclusions is currently incomplete, as key comparisons are confounded by differences in genetic background, and the cellular source of receptor expression is not sufficiently resolved. While the work offers mechanistic insights of interest to researchers studying retinal degeneration, additional experimental controls are needed to fully support the authors' claims.

    2. Reviewer #1 (Public review):

      In this study, Telias et al. identify the P2X7 receptor as a key component of retinoic acid signaling-mediated remodeling in the degenerating rd1 retina. The authors report increased P2X7R expression in the inner retina following photoreceptor loss and link P2X7R signaling to retinal ganglion cell hyperactivity, membrane hyperpermeability, and altered calcium homeostasis. Genetic deletion of p2rx7 abolishes ganglion cell hyperpermeability and reduces several features of pathological remodeling in the rd1 retina. The study provides valuable mechanistic insight into retinal remodeling following photoreceptor degeneration. However, several issues currently limit the strength of the conclusions. In particular, key comparisons are confounded by differences in genetic background, the cellular source of P2X7R expression is not sufficiently resolved, and several experiments require additional controls and more cautious interpretation.

      Major comments:

      (1) The Methods state that C57BL/6J mice were used as wild-type controls, whereas rd1 mice were maintained on a C3H/HeJ background, and the rd1-p2rx7 knockout line is on a mixed background. Direct comparisons among these groups are therefore potentially confounded by strain-specific differences. Authors should use littermate rd1 het mice as healthy controls in all their experiments.

      (2) The authors do not demonstrate P2X7R expression in RGCs. The P2X7R signal in the rd1 retina shown in Figure 1B appears saturated and is therefore difficult to compare directly with the WT image in Figure 1A. Furthermore, Figure 1E indicates that overall P2X7R fluorescence in the GCL is not significantly different between WT and rd1 retinas, whereas the representative images appear to suggest a marked increase. The authors should provide images acquired and displayed under identical settings and consider including retinal whole-mount staining with an RGC-specific marker. P2X7R abundance should then be quantified specifically within identified RGCs in both healthy and rd1 retinas. As mentioned above, het rd1 mice should be included. As an additional control, the authors also should include the staining of rd1 p2x7r KO retinas.

      (3) Does the increase in P2X7R abundance correlate with photoreceptor loss? Please include staining of younger rd1 mice along with het rd1 littermates.

      (4) In Figure 1G-H, the description of the reporter is internally inconsistent: the Results refer to an artificial mini-Pax6 promoter, whereas the Figure 1 legend describes the Ple344 neuronal mini-promoter derived from Tubb3. Please clarify which promoter is used and what cell population the ECFP signal labels. Figure 1G should explain the function of each reporter element and how RAR activity is inferred. Figure 1H should include an RGC marker such as RBPMS and provide quantitative analysis of RBPMS-positive, reporter-positive, and Yo-Pro-positive cells. Additional controls are needed to exclude effects of viral transduction or retinal inflammation on Yo-Pro uptake. They should include rd1 retinas without AAV, rd1 het retinas with and without the reporter AAV.

      (5) In addition, lines 143-144 state that two experiments were performed, but only one is described in that paragraph; the text should be reorganized or clarified.

      (6) Yo-Pro-1-positive cell density in rd1 mice between Figure 1H and Figure 2D is different. Why?

      (7) Constitutive deletion of p2rx7 may cause developmental or compensatory changes that could contribute to the observed phenotype in Figures 2 and 3. Additional controls are therefore needed to distinguish acute effects of p2rx7 loss from developmental consequences. The authors should assess whether p2rx7 deletion alters retinal cell-type composition, including RGC density, or affects the timing or extent of photoreceptor degeneration. Perform a rescue experiment to determine whether overexpression of p2rx7 in the knockout background restores the phenotype. For all experiments, het rd1 control mice should be included.

      (8) Figure 3D and E experiments should include control AAV expression such as GFP.

      (9) Figure 5 experiments should include control rd1 het mice.

    3. Reviewer #2 (Public review):

      Summary:

      In this study, the authors used genetic, transcriptomic, imaging, and electrophysiological approaches to investigate the role of P2X7R in retinal remodeling in rd1 mice. The authors propose that RA signaling upregulates P2X7R, leading to altered Ca²⁺ signaling, HCN1 expression, and spontaneous RGC hyperactivity.

      Strengths:

      Multiple lines of experiments were conducted, focusing on an important question in the field.

      Weaknesses:

      However, some aspects of the experimental design, statistical analysis, and interpretation require clarification. In particular, the genetic controls and causal evidence should be strengthened before the proposed RA-P2X7R-Ca²⁺-HCN1 pathway can be fully supported.

      Some major issues:

      (1) The Abstract states that P2X7R deletion prevents the upregulation of RA-responsive genes, whereas the Results state that RAR-dependent genes were not consistently changed by P2rx7 deletion. This central statement should be corrected and clarified.

      (2) The P2rx7 knockout model requires further discussion. JAX strain 005576 targets exon 13 and has previously been reported to retain truncated P2X7 transcripts with residual activity (Masin et al., 2012). The authors should avoid describing this allele as complete loss of the entire P2rx7 gene unless additional isoform-specific validation is provided.

      (3) The evidence that RAR directly regulates P2rx7 transcription remains incomplete. Predicted promoter motifs and reduced P2X7R protein after BMS-493 treatment are supportive, but they do not establish direct transcriptional regulation. Measurement of P2rx7 mRNA, RAR promoter occupancy, or promoter mutagenesis would strengthen this conclusion.

      (4) Several statistical results require verification. In Figure 2H, four paired eyes are analyzed using an unpaired Mann-Whitney test, and the reported p<0.001 is difficult to reconcile with n=4. Similarly, the significance levels in Figure 5A-D are not compatible with a two-sided Wilcoxon rank-sum test using n=3 mice per group. The statistical unit, exact P values, and number of biological replicates should be rechecked.

      (5) The overall causal pathway remains partially inferential. The study does not directly show that increased Ca²⁺ causes HCN1 upregulation or that HCN1 is required for RGC hyperactivity. These steps should either be experimentally tested or presented as a proposed model rather than an established mechanism.

    1. Help make your website accessible

      Best practice. A button with a distinctive call to action and high visual contrast, making it easier for people with low vision to locate.

    2. Summarize full blog with:

      Good practice. It is connected to other tools—such as ChatGPT or Perplexity—that facilitate reading or text comprehension for readers who might otherwise struggle with the material.

    3. 1. Web accessibility increases your client base and boosts your reputation.

      Good practice. Organizing information using headings and numbered lists or markers makes it possible to process the complexity of the webpage in a clear and predictable way.

    4. Written by: accessiBe Team

      Good practice. Images must include descriptive labels so that people with visual impairments can perceive the visual content using screen readers usch as the picture of the writer "AccessiBe Team".

    1. Typeface Identification

      It's a lot easier to identify typefaces of typewriters if you include a few things in your photos: <br /> - Full frontal picture of the machine with clear shot of one or more of the rulers (so one can attempt to determine pitch as well as more easily identify the make/model for cross referencing).<br /> - A full typed sample of all the characters, especially some of the more unique ones like Q, 3, 5, 4, 7, &, G, A, E, and W in both upper and lower case.<br /> - A close up shot of the slugs which clearly shows the foundry marks (generally found in the middle of the slug) of both letter slugs and special character/number slugs as sometimes the latter are different.

    2. Princess 300 with RaRo Elite Imperial typeface

      The "AR" foundry mark is a shorthand for Albert Rodrian, one half of RaRo. If you ever see an R inside an octagon that's the mark for the other half, Alfred Ransmayer.

      The "84" on the foundry mark is typical of RaRo's 12 pitch Elite Imperial which can be found with a sample at: https://typecast.munk.org/2023/02/04/1961-adler-typewriter-type-font-styles/

      You can also compare other known Princess typefaces at the database: https://typewriterdatabase.com/twdb.0.typefaces?mfr_search=70&model_search=100#

      If you put a ruler up to your type, you should find that there are 12 characters per inch.

      You can also cross reference the following: - the Hass catalog, page 3, row 1:<br /> https://archive.org/details/typewriter-foundry-marks-haas/page/n1/mode/2up<br /> - RaRo Catalogue Page B3-4 "Elite-Schatten Ro 84":<br /> https://www.sommeregger.name/typewriters/files/RaRo/rarotype_katalog_1980_OCR_20160523.pdf


      Reply to u/CaptainCulinary at https://old.reddit.com/r/typewriters/comments/1wjb0gk/princess_300_with_typeface/

    1. These issues must be addressed urgently, in the mathematical community, by the companies developing these technologies and, more broadly, by a society that will confront similar problems in many other forms of intellectual work.

      【局限】作者指出这些问题需要在数学界、AI开发公司以及更广泛的社会层面紧急解决,反映了这一问题的复杂性和紧迫性。这一局限表明AI对数学领域的影响只是更广泛社会变革的前奏。

    2. Whether these changes ultimately benefit the field or have a destructive effect will in large part be determined by the decisions of the humans in control of this new technology.

      【非共识】作者强调AI对数学领域的影响最终取决于控制这项技术的人类决策,而非技术本身。这一观点挑战了技术决定论,强调了人类选择和价值观在塑造AI发展方向中的核心作用。

    3. We are witnessing a general threat to intellectual work, with misalignment between the outcome of the use of AI and its initial purpose.

      【局限】作者指出AI使用结果与其初始目的之间存在错位,这对智力工作构成了普遍威胁。这一局限不仅限于数学领域,而是扩展到所有创造性工作,反映了当前AI技术应用的系统性问题。

    4. In many fields and activities, years of training have traditionally served not only to produce a final answer or product, but also to develop understanding and the ability to formulate new questions and ideas.

      【方法】作者阐述了传统训练的双重目的:不仅产生最终答案或产品,还发展理解能力和提出新问题的能力。这一方法论视角强调了学习过程本身的价值,而非仅关注结果,这与当前AI只关注结果输出的模式形成对比。

    5. Without the willing mathematicians who must take care of their development and integration into the mathematical canon, AI-conceived ideas would never become fully alive and the crucial human transmission chain between mathematicians would be lost.

      【非共识】作者认为没有数学家的积极参与,AI产生的想法永远不会真正"活起来",人类数学家之间的关键传承链将会断裂。这一观点挑战了AI可以独立推动数学进步的假设,强调了人类在数学知识传承中的不可替代性。

    6. Often these solutions are announced in a rush, leaving no time for a proper writeup, the isolation of new methods and ideas, and citing relevant previous work of others.

      【局限】作者指出AI解决方案常常仓促公布,缺乏适当的文档编写、新方法与思想的提炼以及对他人相关工作的引用。这一局限揭示了当前AI生成数学内容的质量控制和学术规范缺失问题。

    7. Solving one of these problems has been a certain sign of new insights and interesting methods, which would then be studied by a community of mathematicians, through a long and arduous process of talks, discussions, simplifications.

      【方法】作者详细描述了数学问题解决的理想过程:解决新问题带来新见解和方法,然后通过学术社区的长久讨论和简化过程最终形成教科书级内容。这一方法论强调的是数学作为人类集体知识的演进过程,而非单纯的答案获取。

    8. The goals of the AI companies and the goals of the mathematical community are severely misaligned.

      【非共识】作者明确指出AI公司与数学社区的目标存在严重错位,这是一个非共识观点。在技术乐观主义者看来,AI应该辅助而非阻碍数学发展,而本文认为商业利益与学术追求的根本冲突可能导致数学研究的质变。

    9. Over the last few months, the mathematical capabilities of LLMs have improved dramatically, to the point that they can solve major outstanding problems in many fields of mathematics.

      【数据】作者指出LLMs在数学能力上取得了戏剧性进步,能够解决许多数学领域的主要未解决问题。这一陈述暗示AI已经达到可以独立解决重大数学难题的水平,但未提供具体数据支持,这可能是一个需要验证的关键主张。

    10. I am proud to be among the list of 25 initial signatories — all Fields Medallists — to the declaration below

      【非共识】这篇声明由25位菲尔兹奖获得者共同签署,代表了数学界最高权威对AI的集体担忧。这种顶级数学家的集体发声本身就是对AI在数学领域应用的强烈质疑,与主流科技界对AI能力的乐观态度形成鲜明对比。

    1. eLife Assessment

      This important work considerably advances our understanding of the nature of ongoing thought across individuals and across task contexts. Using a combination of fMRI brain data and experience sampling, they find that activation of the MDN leads to greater stability of deliberate task-focused thoughts during task performance. The evidence supporting the claims is convincing and both analytically strong and creative, offering a novel approach to capture thought stability and its relationship to independent brain maps. The findings will be of broad interest to cognitive scientists and neuroimagers.

    2. Reviewer #1 (Public review):

      Summary:

      This paper suggests an alternative model for the function of the multiple demand network. Specifically, its role is not necessarily to sustain cognitive control and maintain task sets, but rather to "stabilize task-appropriate modes of thought". Evidence for this would be that the MD network is responsible for maintaining a particular thought state during a task. To investigate this, they use a combination of fMRI brain data during a set of 14 tasks and experience sampling in a different set of participants performing the same tasks. Using dimensionality reduction, they reduced the space of task features (and brain systems) to a smaller, more tractable set of dimensions and examined whether stability in specific thought components was related to recruitment of specific brain systems during the task.

      Strengths:

      Overall, this is an interesting and creative study with strong analytic methods that do a good job accounting for confounds or alternative explanations (save one I mention below).

      Weaknesses:

      I have mostly minor comments and one major one.

      Major:

      The principal finding is that tasks that evoke brain activity patterns that resemble the MDN also had more stable "deliberate task focus" features. While all the analysis and controls are impressive, I'm still left with the sense that this is reifying something we already know or that alternative explanations are more parsimonious than the MDN induces stability in "thought".

      I thought an example might be easiest to understand my point: If I gave participants a series of working-memory-related tasks. Some of these are the crème de la crème, and others are sloppy and poorly designed. Then suppose I assess them on measures related to deliberate thought; I'd likely find the "good" tasks elicit more consistent/reliable deliberate task focus. I also would bet money that these same tasks would evoke canonical WM and MDN activity patterns more than the sloppy tasks. This isn't evidence of MDN stabilizing patterns, but rather that both stable thought patterns and activity in the MDN share a common cause. Thus, I would predict that with my thought experiment, your analysis would find the same result. So it strikes me as a strong alternative possibility for these results is that tasks that reliably evoke deliberate task focus are also those that more strongly and consistently evoke working-memory demand (i.e., Figure 2 shows that they are primarily driven by the executive/WM tasks).

      Minor:

      (1) The PCA was reviewed previously, and I don't want to relitigate a prior method, but I had one minor concern. It would be useful to know how the principal results are based merely on the "deliberate" or "focus" items specifically. Is the thought space necessary, or do the individual items that likely drive the "deliberate task focus" PC essentially replicate the main result?

      (2) I struggled with the motivation for projecting the task data onto a resting state FC analysis that focuses on "gradients". I understand that with 14 tasks activation maps, data reduction is a good thing. But as someone who isn't as enmeshed in this work, I didn't follow why this specific "atlas" was chosen over any other (parcellations, meta-analytic maps of canonical networks, etc.). Maybe a brief sentence saying why this and not that would help readers who find themselves in my shoes.

    3. Reviewer #2 (Public review):

      Summary:

      The study's aim was to establish whether stability in thought patterns relates to the particular thought pattern, the task context, or their interaction. And further, whether stable thought patterns could be linked with distinct brain patterns.

      Strengths:

      The core reliability framing is novel, and the trait/state/interaction decomposition is a fruitful way to pose the question, leading to the finding that stability is neither a pure trait nor a pure task property, but emerges from their interaction, which is a solid contribution.

      The aim to characterise aspects of stability across individuals and across tasks was achieved and is well supported.

      Weaknesses:

      While the paper makes excellent use of existing data sets, the independent samples, i.e., one study sample for the cognitive/thought-sampling data, and several different study samples to generate the brain maps, do limit the brain-behaviour conclusions that can be drawn.

    1. eLife Assessment

      This important study presents compelling evidence for dopamine-based representation learning in the mouse olfactory tubercle. While this work suggests that gradient descent in deep networks serves as a helpful framework for understanding neural representations, the empirical evidence provided remains incomplete. Nonetheless, this work will be of significant interest to systems and computational neuroscientists, cognitive scientists, and machine learning theorists.

    2. Reviewer #1 (Public review):

      Summary:

      This manuscript reports on simultaneous neural recordings in the olfactory tubercle (OTu) and ventral tegmental area in head-fixed mice performing a simple go-nogo odor-guided reversal task. In the task, there were 3 odors that predicted lick-spout water at 0%, 50%, and 100% probability, with 0% and 100% odors reversing at some point each session. The authors fitted this neural data to a value function approximator in which reward prediction errors fed back onto state representations, allowing the optimal set of representations to be learned. They found that model predictions of adjustments to state representations were correlated with trial-by-trial changes in OTu neural activity, from which the authors conclude that such a system, with dopaminergic errors feeding back onto OTu representations of states, which then generate reward predictions, is biologically plausible.

      Strengths:

      This is a novel and creative modeling approach that has important implications. It seems to be showing the biological plausibility of a model that can learn state representations rather than relying on a fixed set of states that are programmed into the model. This makes tremendous sense, because the real world is much less well-defined than the kinds of tasks conventionally used by neuroscientists to probe reinforcement learning. As such, it is an important demonstration.

      Weaknesses:

      The task used in this study is quite simple in its state space, and in particular in how it maps sensory stimuli (odors) onto states, such that the task would seem not to require a system that can learn state representations, or at least that it would not be ideal for testing such a model. This mismatch raises some questions about why this model would perform as well as it appears to be doing here.

      The model has two updating functions, both using dopaminergic RPE's. One of these maps raw stimuli to state representations using a parameter termed theta; the second maps state representations to a value prediction, using a parameter termed w. The interaction of these two updating functions seems to be giving the model its interesting characteristics. But it is critical to test what the first of these updates is doing in this model, given that odor stimuli appear to map straightforwardly onto states. For example, one might test the effect of ablating this part of the model, leaving only updates of what the authors term w. Relatedly, one might test the extent to which OTu neurons show simple odor selectivity before and after reversals, asking whether these neurons reflect state representations in this model merely by being selective for a particular odor, or if they develop a more complex kind of responsiveness.<br /> The authors compare the full model, which uses a gradient descent update, to a series of alternatives. The fact that the full model performs significantly better than any of the alternatives leads to the conclusion that the brain is using something like this model in this task. But these alternative models are all reduced or simplified versions of the primary model. This suggests that the full model is the best version within the basic framework posed by the authors. But to draw the conclusion that this model is capturing what is occurring in the brain, one would want to test how this model would perform compared to a different class of model, in particular one that assumes a fixed set of states.<br /> A second weakness of the paper is that the authors do not show the behavioral or raw neural data, which would be important to summarize for the sake of transparency and to help readers get an intuitive sense of what is going on in the task and why the model performs as well as it does. One essential issue is: how much training do mice receive before neural data used in the analysis are collected? Do mice get pre-exposure to contingency reversals before analyzed neural data is collected? Relatedly, how quickly (i.e., in how many trials) do mice show behavioral evidence of having learned initial contingencies and then reversed contingencies? What is the behavioral criterion? With regard to the number of trials mice take to learn the reversals, this can change enormously over training, and such changes could have a big effect on how the model performs. Regarding the neural data, one would want to show some measure of odor selectivity of SPN's and DAN's and how each population responds to delivery and omission of reward in different conditions.

    3. Reviewer #2 (Public review):

      Summary:

      In this paper, the authors use electrophysiological recordings from the olfactory tubercle (OTu) and ventral tegmental area (VTA) of mice learning an olfactory Pavlovian conditioning task to demonstrate the consistency of changes in OTu odour responses with gradient descent updates minimising reward prediction error (RPE). The paper is clearly written and the work well motivated. The authors address a gap in the literature on animal reinforcement learning by providing neural evidence for gradient-based representation learning, something that had been proposed but not yet tested. The results are convincing and the limitations comprehensively addressed. Of particular interest is the proposal that OTu SPNs could solve the weight transport problem through knowledge of the sign of their downstream connections from their expression of either D1/D2 receptors. This makes a concrete experimental prediction that future research could test.

      Strengths:

      The paper provides one of the first demonstrations backed by neural recordings that representation learning in the brain is consistent with gradient descent. It shows how, although weight transport may be biologically implausible, the brain appears to find other ways to compute a gradient in multilayered networks for efficient learning. The study builds nicely on recent work in systems neuroscience and provides evidence for a concrete implementation in the OTu-VTA circuitry of mice.

      The paper demonstrates that changes in OTu striatal projection neuron (SPN) activity over trials are proportional not only to the RPE relayed by VTA dopamine neurons, but also account for the influence of each particular SPN on the RPE. If an SPN decreases the RPE when active, its activity will increase on the next trial after a positive RPE. The activity will instead decrease for an SPN that increases the RPE. This relation is encapsulated in the update rule of Equation 5.

      Weaknesses:

      The main weakness of the paper is that it provides only indirect evidence for the update mechanism by inferring synaptic weights based on the (justified) assumption of VTA dopamine neurons encoding RPE. This is still a substantial contribution to understanding representation learning in the brain, though I do think that the authors could provide some additional evidence to further convince readers. One idea could be to analyse the distribution of inferred weights and validate whether it agrees with known statistics of connectivity between OTu SPNs and their downstream projections (e.g. fraction of D1/D2 SPNs).

      Further, as dopamine neurons are known to have asymmetrical responses for positive and negative RPEs (with the dips in activity related to negative RPEs being generally smaller) I'd expect an improvement in the correlation of the learning updates particularly after the reversal if the authors account for this in the model.

    4. Reviewer #3 (Public review):

      Most models of reinforcement learning in the brain treat the question of how the external world is represented as an afterthought. In "Error driven representation learning in the mesolimbic system", the authors begin by calling attention to this limitation of previous work, then proceed to show that neural activity in part of the ventral striatum evolves in a manner consistent with dopaminergic RPE-driven representation learning. Overall, the question is interesting, the modeling is well done, and the claims are bold. While I am not convinced that olfactory tubercle outputs _mainly_ reflect state features acquired through error-driven learning (see main points below), I am now more willing to believe that they might. I am confident that this work will spark discussion.

      Strengths:

      The latent weight trajectory inference approach is a nice application of Kalman filtering, the model validation and comparison steps are well executed and thorough, the discussion is nicely written and includes an appropriate caveat about the weight transport problem, and the general idea of striatal output being value-like but also having state-like aspects that evolve over time is thought-provoking.

      Weaknesses:

      In my view, the main limitation of this work is that the authors do not clearly rule out (1) non-representation learning and (2) non-error-driven learning. This fits with the narrow research question stated at the end of the introduction (L71-72), which is confirmatory in nature and does not claim to exclude alternatives, but other aspects of the framing are less consistent.

      Main points of criticism

      (1) Representation vs. value learning

      The first two pages of the manuscript gave me the impression that the authors wish to draw a clear distinction between representation and value learning, and that they would squarely position this paper as a study of representation learning. The substance of the work does not seem consistent with this positioning, a mismatch that could be addressed either by incorporating new analysis and discussion or by changing the framing.

      Examples of emphasis on representation:

      - Abstract L14-17 defines value and representation learning and states why they are different.

      - First three paragraphs of introduction explain the power of learned representations.

      - Results L103-108 attribute state representations to OTu and value to downstream regions.

      For this level of emphasis, it would be good to see a convincing argument that OTu MSN activity is better understood as a state representation than as the output of a value function.

      Options for establishing OTu as state and not value:

      - Explicitly claim in the introduction that previous work has established this, and explain why evidence of stimulus valence being encoded in this area (citations on L68-70) does not favour the value interpretation. Discussing how OTu differs from other parts of ventral striatum that are canonically seen as value coding would also help.

      - Directly compare the extent of state vs. value encoding in the present data. The $w=1$ control is a good step in this direction, but its connection to value coding is mentioned only in passing.

      (2) Error driven learning vs. other types of learning

      Similar to my previous point, the authors seem to claim that the representational changes they study are specifically error driven. While the authors include a good number of controls and ablations, it was not obvious to me that any of them correspond to a form of non-error driven learning that could plausibly generate useful representations. Adding Hebbian learning or a sparsifying learning rule would strengthen this aspect of the work.

    1. Welcome to my portfolio. Here you’ll find a selection of my work. Explore my projects to learn more about what I do.Camping Planner AppArchitecture App UI DesignZodiac T-shirtWood engraving reliefTwenty-Four Heavenly Guardians (二十四诸天, Èrshísì Zhūtiān)"Imperial Plumage" IntaglioAnatomy collection Watermelon print practise Green_pepper#1Cell_Phone

      I think there is a lot of good work on this page, but everything is shown together, so it feels a little messy. I would organize the projects into different categories, such as App Design, UI Design, Sketches, T-shirt Design, and Chinese-inspired Design. This would make the portfolio easier to understand and help users find specific projects more quickly.

    1. Climate Feedbacks Climate feedbacks are natural processes that respond to global warming by offsetting or further increasing change in the climate system. Feedbacks that offset the change in climate are called negative feedbacks. Feedbacks that amplify changes are called positive feedbacks. Water vapor appears to cause the most important positive feedback. As the earth warms, the rate of evaporation and the amount of water vapor in the air both increase. Because water vapor is a greenhouse gas, this leads to further warming. The melting of Arctic sea ice is another example of a positive climate feedback. As temperatures rise, sea ice retreats. The loss of ice exposes the underlying sea surface, which is darker and absorbs more sunlight than ice, increasing the total amount of warming. Less snow cover during warm winters has a similar effect. Clouds can have both warming and cooling effects on climate. They cool the planet by reflecting sunlight during the day, and they warm the planet by slowing the escape of heat to space (this is most apparent at night, as cloudy nights are usually warmer than clear nights).Climate change can lead to changes in the coverage, altitude, and reflectivity of clouds. These changes can then either amplify (positive feedback) or dampen (negative feedback) the original change. The net effect of these changes is likely an amplifying, or positive, feedback due mainly to increasing altitude of high clouds in the tropics, which makes them better able to trap heat, and reductions in coverage of lower-level clouds in the mid-latitudes, which reduces the amount of sunlight they reflect. The magnitude of this feedback is uncertain due to the complex nature of cloud/climate interactions.9

      Lastly, the webpage also contains distinguishable and readable content. Almost the entirety of the article is done with black text over a white background making it clearly readable, whereas the images and videos are placed in a proper manner far apart from the text without making it all clumped up.

    2. Key Greenhouse Gases

      The webpage also includes a video titled "The Greenhouse Effect" (Right above the highlighted title). This video shows commendable web accessbility practice as the video has closed captioning on videos, and the important sounds and contexts are communicated nicely.

    3. The greenhouse effect helps trap heat from the sun, which keeps the temperature on earth comfortable. But people’s activities are increasing the amount of heat-trapping greenhouse gases in the atmosphere, causing the earth to warm up.

      The article also has images and videos. These are very important to consider when looking at web accessiblity because these types of media need to be accessible to different users. (Images and videos can make it accessble to a larger visual audience).

    4. Basics of Climate Change Learn about some of the key concepts related to climate change: The Greenhouse Effect Key Greenhouse Gases Other Greenhouse Gases Aerosols Climate Feedbacks

      The webpage has the information organized using clear and proper heading and subheadings. The linked subheadings take the user to their preferred topic about climate change, and make it very easy for them to navigate and understand the entire article. Nice web accessbility practice!

    1. Welcome to my professional Art Portfolio. One of my defining principles is establishing an open dialogue between form and function in all aspects of my artwork. Look through my work to see this in action, and get in touch for more information.

      The white text is placed over a busy background image, which can make some parts of the text difficult to read. Adding a darker overlay could improve the contrast. I also think the background image should be one of the artist’s most representative works because it is the first image visitors see. This could help people understand the artist’s style more quickly.

    2. HomeCamping app concept appSelf-ordering web appFirst-home field guidePortfolioExperiencesQualificationContact

      The navigation titles are clear and specific, so users can understand what each page is about before clicking on it. This also helps screen reader users navigate the website more easily because each link clearly describes where it will take them.

    3. Home

      The color of the “Home” link is light purple, whereas all other navigation links have a black color. In my opinion, this is a good solution since a different color enables users to easily see that the current page is the Home one. Moreover, this color makes the navigation unique. On the other hand, it may be difficult to see the light purple color against a white background for users with vision problems or color blindness. Therefore, I would suggest making the “Home” text slightly bold so that it is easier to read while still keeping the purple colour.

    1. I have never heard of a finer ship  fitted with the weapons and armor of war,  swords and harnesses. In its embrace lay  a multitude of treasures, which were to go with him  far off, into the dominion of the sea.

      this line emphasizing that no ship had ever been more magnificently outfitteds for war and the after life.

    2. here were many treasures  from faraway lands, such precious things loaded there.

      This is similar to the dicussion we had in class classifying what treasure is.

    3. They laid down their dear king,  giver of rings, in the bosom of the ship,  mighty by the mainmast.

      the traditional Anglo-Saxon bond of loyalty between a lord and his loyal thanes through generous gift-giving

    4. His beloved companions carried him then to the  water’s edge, as he himself had instructed  when he still governed, that much-loved Scylding friend,  their beloved land-prince held power a long time.3

      This Sheild sheafson grand ship burial one of the most vividly detailed funerals in all of world literature.

    5. Through glorious deeds  a man shall prosper among peoples everywhere.

      In a warrior society immortality was achieved not through a long life but through praise and fame generous leadership.

    6. He consumed honors until each of the other surrounding  tribes over the whale’s road were forced to obey him and pay tribute.

      Forcing neighboring nations across the whales road to pay tribute wad the ultimate proof of a kings strength

    7. The Lord of life,16b ruler of glory, gifted worldly honour:  Beowulf was famed with widespread renown,  son of Scyld, in the northern lands.

      this line introduces Beowulf the son of Sheafson being an early Danish king. The foundation anglo-saxon is ideal of a young prince should earn loyalty through gift

    1. Canada’s prairie provinces export the most honey of all of Canada’s provinces. About 56 percent of the honey Canada exported last year went to the United States; it is the top foreign market for Canadian honey,

      At what point is the imposed tariff enough to stop Canadian producers from exporting their products? How detrimental is that to their overall income, and is only selling domestically sustainable for the Canadian economy?

    2. “If you’re picking fights with Canadians, then you don’t have a friend in the world,”

      To what extent can Canada position itself on the world stage as a victim to United States aggression, and will that have an impact on other nations trading with the United States?

    3. “We’re standing our ground not to be bullied on the world stage and pushed around,”

      I find it interesting how the tariff war is just as much about national identity and pride as it is about money and trade. Does Canada want to uphold their national pride and image on a world stage to incentivize trade from other nations? I wonder how Canada can stand their ground against America tarriffs without taking too large of a hit to their economy

    1. Hello my name is Alexander Reyes, I chose communications because previously I was a Psych major, who really only enjoyed the aspects of Psych that involved communication and interaction between people. I then decided it would be best if I took a step in a direction that was finally for me, and landed on the Communications Major. One thing I think our major could improve on is class sizes. Yes please add 500 people into one intro to argumentation class without giving the Prof. a raise that is how you do it! Oh and my high school mascot was a Cougar lol.