10,000 Matching Annotations
  1. May 2026
    1. On 2024-11-21 16:07:10, user Phillip Gienapp wrote:

      This manuscript has now been published as:

      JJC Ramakers, TE Reed, MP Harris & P Gienapp (2023) Probing variation in reaction norms in wild populations: the importance of reliable environmental proxies. Oikos e09592 doi: 10.1111/oik.09592

    1. On 2024-11-21 06:52:29, user Eva wrote:

      Dear authors,<br /> Thx for sharing the study.

      We have two doubts about the results:<br /> 1. During the RNA/Chimera annealing step, you applied “The samples were incubated at 95℃ for 1 minute and transferred to ice immediately.” However, the original study of Nm-VAQ and other similar experiments usually use gradient and slow cooling to ensure that the RNA and chimera have sufficient time to specifically bind rather than preferentially form secondary structures on their own, thereby affecting the subsequent cleavage of RNase H. This may be the reason for the ‘background’ described in this study and the much higher level of Nm modification seen in Fig 3C than the results in the original article.

      1. In page 7, “PLXNB2, RPL10, and TOP1 were identified as potential candidates for mRNAs with Nm modifications based on various mapping techniques. G498 for PLXNB2, G150 for RPL10, and U877 for TOP1 were the potential sites of 2’-O-methylation”, which only cited reference 29. However, we found ref 29 only reported PLXNB2 as a potential Nm site. Could the authors please point out the provenance of the other two test loci?

      Appreciate your attention.

    1. On 2024-11-20 15:57:30, user Miriam Zagers wrote:

      We recommend that you read our follow-up study as well. In this study, we developed and tested a human embryo culture medium based on these in vivo results, on human preimplantation embryos. The blastocyst formation rate in the culture medium mimicking the uterine environment was similar to that of the control medium. DOI: https://doi.org/10.1101/2024.11.19.623474

    1. On 2024-11-20 15:06:42, user Nicolas Locatelli wrote:

      This manuscript has been published in BMC Genomics w/ open access on 11/20/2024 with the title: "Chromosome-level genome assemblies and genetic maps reveal heterochiasmy and macrosynteny in endangered Atlantic Acropora". The "now published in..." link at the top of this preprint will be coming soon.

    1. On 2024-11-19 16:24:35, user Jesus Valencia wrote:

      An exhaustive analysis of viral diversity in blackflies collected from Onchocerciasis-endemic areas is presented. While it provides valuable data, the study has limitations that affect the coherence between its main objective and results. The title and introduction suggest a focus on neurotropic viruses associated with epilepsy linked to Onchocerca volvulus (OAE), but the results primarily describe viral diversity without providing evidence of these viruses. Possible explanations for the absence of neurotropic viruses are not discussed, nor are methodological alternatives proposed to address this issue. To improve the narrative, it is recommended to reframe the title and present the work as an exploratory study on viral diversity in blackflies.

      Methodologically, the sample size (10 blackflies per pool) and the lack of information about the genus of the collected blackflies compromise the representativeness of the data. Additionally, the criteria for excluding incomplete sequences are inconsistent, as incomplete Riboviria sequences were retained, potentially biasing abundance results.

      Regarding the figures, Figure 2 lacks a clear description of its columns, hindering interpretation. Figure 3 presents phylogenetic trees with barely distinguishable details, complicating analysis. In Figure 4, the point corresponding to Blackfly flavivirus does not clearly cluster with invertebrate hosts, weakening conclusions about its specificity. Analyses based on PCA and dinucleotide composition are interesting but insufficient to infer viral tropism without experimental validation. The supplementary files, although relevant, lack proper descriptions, and the main text does not specify which supplementary figure to consult, limiting their utility.

      Although the study does not successfully address its main objective, it provides an important starting point for future research on the viral diversity of blackflies and highlights the need for further investigation into their relationship with human diseases in endemic areas.

    1. On 2024-11-19 16:05:20, user Jose Morales Perez wrote:

      Dear authors,<br /> I am a biology student with an interest in the evolution of animal communication systems, and out of curiosity, I came across your research titled “Prediction of the Next Elements in Chimpanzee Gesture Sequences.” I would like to congratulate you on your innovative approach and the richness of the information you present. The study is very interesting and makes a significant contribution to the field of animal communication, especially by providing a quantitative approach to chimpanzee behavior. However, I have identified a few aspects that could be improved to make the report clearer and more accessible.

      First, the text refers to a "Table 2," but this table does not appear in the report, which could cause confusion for readers. This might be an issue of incorrect enumeration or an omission, so I would recommend reviewing and correcting this to ensure consistency in data presentation.

      Regarding the figures, while the information is clear and understandable, some of them could benefit from adjustments in formatting to improve their visual presentation. For instance, in Figure 4, the axes are generically labeled as “X” and “Y”; it would be helpful if these were more specifically defined to aid in interpretation. Additionally, the axis titles are in lowercase, which could be standardized for consistency. In Figure 6, the title is not centered, and some axis subtitles and the legend do not follow the proper formatting, which may make them harder to read. In Figure 8, the abbreviation “NB” is used without an explanation, which could cause ambiguity for readers unfamiliar with the term.

      Lastly, I believe it would be useful to include a section for keywords, as they are crucial for indexing and searching for the article, thus improving its accessibility and understanding.

      The work is of high quality, and with these small adjustments, it could further enhance its presentation and clarity. Once again, congratulations on your effort and contribution to the field.

    1. On 2024-10-05 08:38:14, user Martin GIURFA wrote:

      Great work, congratulations! <br /> You may be interested in having a look at the following works, which relate to your findings:

      ° Pheromones modulate reward responsiveness and non-associative learning in honey bees. Baracchi D, Devaud JM, d'Ettorre P, Giurfa M. Sci Rep. 2017 Aug 29;7(1):9875. doi: 10.1038/s41598-017-10113-7

      ° Pheromone components affect motivation and induce persistent modulation of associative learning and memory in honey bees. Baracchi D, Cabirol A, Devaud JM, Haase A, d'Ettorre P, Giurfa M. Commun Biol. 2020 Aug 17;3(1):447. doi: 10.1038/s42003-020-01183-x.

      Good luck with the next steps!

    1. On 2024-11-17 07:22:07, user Velyudhan Mohan Kumar wrote:

      This study, identifying orthologs of human sleep-related genes in the Lokiarchaeota will certainly help in understanding human sleep physiology. The study supports the concept that sleep is a conserved biological function and evolutionarily distinct new features have developed later. The findings that the genes associated with energy metabolism and cellular repair, which play a primordial role of sleep in cellular maintenance in the Lokiarchaeota, suggest the origin and evolutionary purpose of sleep. The results are in line with the previous assumptions regarding the restorative and protective functions of sleep, including cellular maintenance and gene protection. Moreover, such researches can provide valuable information for drug development as conservation and evolutionary background of genes linked to disease could lead to the discovery of novel therapeutic targets.

      This work has used both phylogenomics and bioinformatics to solve a biological problem. This research illustrates the scope of such studies in finding solutions to many biological problems. Why sleep is necessary and is maintained by evolution is not fully answered. From an evolutionary standpoint, phylogenetic approaches would help to unravel some of the unresolved mysteries of nature. The phylogenomics approach to sleep science is rather new, but combining phylogenomics and bioinformatics research on a range of species may help in understanding facts which are still unknown.

    2. On 2024-11-05 08:41:43, user Seithikurippu R Pandi-Perumal wrote:

      We greatly appreciate your constructive feedback and the recognition of Professor Imachi's significant contributions, which we have acknowledged in our manuscript. We fully agree that, ideally, certain wet lab experiments might have strengthened our findings to provide additional insights. However, there was no funding for our ongoing effort. If adequate funding were to be available, we would have pursued a certain set of experiments that we thought would add value to our work. This limitation is transparently discussed in our manuscript, and we consider our current data a valuable "stepping stone" for additional work in this exciting area of research.

      Looking ahead, we are well-positioned to undertake follow-up experiments, contingent upon the availability of appropriate funding and other resources such as access to the strains. However, we do recognize the challenges involved in maintaining and rearing such specialized organisms in the laboratory setting. For example, replicating their native hydrothermal vent conditions and performing wet lab operations under those constraints condition pose additional significant technical and logistical hurdles. These unique challenges underscore the complexity and significance of our research. With support, we would eagerly and carefully advance this aspect of our research in the future.

      Regarding your comment on sleep and circadian cycles in this organism, the identified sleep orthologs underscore that their functional characteristics diverge from those in human sleep processes. From the circadian clock perspectives, organisms whose circadian clocks have periodicities that correspond with those of their cyclic surroundings have survival advantage than those whose periodicities do not. If the development of circadian clocks was driven by circadian resonance, then animals with ~24h periodicity should be more fit in a 24-hour environment than in any other periodic or aperiodic time scale. Alternatively, the persistence of functional clocks and photoentrainment mechanisms implies that circadian clocks provide some "intrinsic adaptive advantage" to organisms in the absence of cyclical environments, where there is no apparent need to synchronize behavioral and metabolic processes with the geophysical cycles ( https://doi.org/10.1186/1740-3391-3-7) . The inherent benefit of clock mechanisms, there, is most likely derived from their ability to help coordinate the organisms' numerous internal metabolic functions.

      In extreme habitats like hydrothermal vents, in the absence of photic cues, organisms likely rely on alternative timing mechanisms, such as temperature or chemical gradients, adapted to their unique environmental niche. In organisms such as Lokiarchaeota, the absence of circadian regulation due to the extreme hydrothermal environment may have driven the evolution of alternative timing mechanisms. Without light cycles, these organisms likely rely on other environmental cues, such as temperature fluctuations or chemical gradients, to regulate physiological processes. This adaptation might provide them with flexible metabolic and behavioural responses suited to their habitat. This suggests that evolutionary forces and selective environmental pressures unique to hydrothermal vents shaped a different evolutionary pathway for metabolic, and chemosynthetic mechanisms for survival than light-regulated, temporal mechanisms seen in surface-dwelling and terrestrial species. It's possible that daily and seasonal variations during evolution further honed the organisms for survival.

    3. On 2024-11-04 12:05:23, user S. Krishnaswamy wrote:

      Sleep in bacteria and archae is usually associated with a state of dormancy or persistency. But that is different from the circadian cycle of sleep. It is therefore surprising to see orthologs of human sleep genes in the deep sea vent dweller like this archae. Unless they had a different function there. It would be interesting to mutate some of the identified genes and see the effect in the cultures grown in the lab of Hiroyuki Imachi

    1. On 2024-11-17 04:42:42, user De Novo Enzyme Enjoyer wrote:

      Fantastic paper, breathtaking results. Just one minor suggestion—when I searched for RFam on GitHub, I found a repo with an almost identical name (Rfam) already exists, for an RNA bioinformatics tool. So to avoid confusion/difficulty locating the Rosetta Flow Atomic Matching codebase (and paper, whenever Ahern et al comes out), it might be a good idea to change the shorthand nickname for the model before publishing the repo (and the papers). RFlam, RoFAM and RoFlam all appear to have more unique namespace than RFam, while still being just two syllables to pronounce. Just a thought!

    1. On 2024-11-15 07:55:21, user Giovanni Bussi wrote:

      Authors of the review

      Olivier Languin-Cattoën, Giovanni Bussi

      Summary

      The authors use Molecular Dynamics with enhanced sampling techniques to gain insight in the directional catch-bond mechanism of Vinculin tail (Vt) interaction with F-actin. They construct two models of the Vt-actin complex that are hypothesized to represent a weak state and a strong state in a force-activated allostery model of the catch bond. They use enhanced sampling techniques to estimate the free-energy landscape of the unbinding process in each state, as well as the unbinding kinetics, in absence and presence of pulling forces of variable intensities and directions along the actin filament axis. Their results demonstrate higher kinetic stability for the strong state with respect to the hypothesized weak one, with unbinding kinetics in range of experimental expectations, confirming the viability of a “3-state” (2 bound states, 4 unbinding pathways) kinetic model of the bidirectional catch bond observed in single-molecule experiments. They additionally observe an increase in bond lifetimes for moderate constant pulling force (10-20 pN) in both directions, indicating an intrinsic catch-bond behavior within each “allosteric” state that may superimpose to the overall allosteric one. They show how an external pulling force affects the positioning of the H1 α-helix believed to act as a regulatory motif of the weak-to-strong transition, providing a compelling structural hypothesis for a force-induced allosteric mechanism. Finally, they provide molecular insight on the difference in stability between both states, highlighting the role of a C-terminal extension (CTE) and the redistribution of Vt-actin contacts under force.

      Comments

      * We wonder if the nomenclature “Holo” can be confusing at first glance for the reader, given the historical usage of the holo- and apo- prefixes to designate protein constructs with and without their constitutive prosthetic groups (usually non-proteic cofactors)

      * We suggest that the authors make clearer the composition of the Holo and Aligned protein sequences given the different numbering in 6UPW and 1QKR (that we guess is due to the presence of Metavinculin instead of Vinculin), that in our understanding are identical except for the absence of the leading H1 helix.

      * In FES calculations, since the constant force (that is, a linear bias) is applied on the same coordinate Q_‖, we expect that the resulting FES could be entirely predicted from the FES at zero force by simple addition of the linear slope, given sufficient exploration of the Q_‖ direction during the OPES-MetaD sampling. This fact could be used by the authors to assess the consistency between the FES computed at different forces. Alternatively, one may want to first aggregate the OPES-MetaD simulations at all forces using appropriate reweighting, and then estimate minimum free-energy paths and free-energy barriers at arbitrary force using the aggregated FES. This approach might lead to a better statistical use of the vast amount of simulations and a smoother estimate of the FES at all forces within the studied range. Finally, the multiple trajectories (20 replicates x 9 force values) could be used in a single bootstrap to assess the statistical uncertainty of the results (see below).

      * Given the coarse nature of the set of chosen CVs (Q_‖ and Q_⊥) it is unclear whether the Vt is able to regain its canonical binding site during OPES-MetaD, notably because of free rotation with respect to the actin filament. <br /> - The authors acknowledge the difficulty in the methods (sec. 5.5) without explicitly stating whether such rebinding events happen at all in their simulations. We believe this is an important piece of information for proper understanding of the FES presented. We would suggest showing time series to clarify how many binding/unbinding events are observed.<br /> - One might expect that the absence of recrossing lead to a poor estimate of the free energy difference between the bound and unbound states as well as the height of the binding free energy barrier. On the bright side, the estimate of the unbinding barrier – which is the one they are the most interested in – can still be expected to be reliable.<br /> - The authors suggest to run multiple (20) separate OPES-MetaD simulations to compensate for this limitation. It should be acknowledged that independent runs starting from the bound state will not correct for the systematic bias caused by the absence of rebinding events. A bootstrap on these replicas would anyway estimate their statistical error.<br /> - We wonder if the use of the more specific CV Q_contact might allow for such recrossing to happen within OPES-MetaD, without the need of aggregating a high number of independent trajectories. In our understanding the authors only used Q_contact to assess the robustness of the free-energy barrier height to the precise choice of the projection space, but did not try to perform OPES-MetaD directly on this CV space, which could be instructive.

      * The authors analyze the FES by determining a minimum free-energy path using the “String method” as a post-processing method. <br /> - The Methods section might benefit from some information about the use of the method, in particular that it is directly applied on the 2D CV-space projected FES (as opposed to a search of a minimum energy path on the full potential energy surface as originally proposed in [50]), and provide details about initialization (choice of end points for the string, number of nodes, initial interpolation) and robustness to these parameters in the converged paths and corresponding barrier estimates.<br /> - Since the FES are aggregated from 20 independent OPES-MetaD runs, it might be relatively straightforward to estimate errors (for example using bootstrapping) and provide error bars on Fig 2c. We believe this would strengthen the significance of the observed barrier difference.

      * One may be concerned about the significance and reliability of the constructed “Aligned” state, since this state was constructed by aligning Vt to another conformation (in what we could refer to as “docking-by-homology”) with little experimental confirmation that such a state is stable in vitro. We understand that the model constitutes the core hypothesis of the whole computational approach, and that the consistency of the computational outcomes with single-molecule experiments themselves validates its plausibility. Nevertheless, it could be argued that lower binding affinity and lifetime are to be expected from a suboptimal binding partner in a suboptimal binding pose. This raises the question of whether the proposed model corresponds to a specific binding mode in reality, or if the results could be reproduced with a different alignment. Is this ruled out by the stability observed in the 500 ns simulations shown in Fig S1?

      * In Sec 3.2 §4, the authors say “these results [...] do not quantitatively explain the observed experimental results, since the experimental changes in lifetime shown in Fig. 1D reflect a net 1.4 kcal/mol change in barrier in the negative direction, and 0.6 kcal/mol in the positive direction (if one assumes a constant prefactor the kinetic rate constant)”. It was somewhat unclear to us where these values come from (Are they computed from the fitted 3-state model in SI S1? At a specific value of the pulling force?) and how exactly they are compared to the computed barriers for Holo and Aligned to conclude to a discrepancy (Overestimated?)

      * In Sec 3.4 §3, the authors convincingly remark that in two out of five simulations of Holo+H1 state pulled towards the barbed end, the conformation of H2–H5 becomes more similar to the Aligned (unbound Vt) structure, suggesting a first step in the strong → weak allosteric transition. We wonder if (i) the specific contacts made with actin and (ii) the specific intra-domain contacts of H1 with the H2–H5 bundle are also indicative of a displacement toward the Aligned state, since this would be an even stronger argument validating the proposed model.

      * Given the relative simplicity of the supposed allosteric motif and the suggestive results of the Holo+H1 simulations, we cannot help but wonder whether the authors also tried to "unfold" the H1 helix in the Aligned model with a N-terminal pulling force since this seems a very natural test to look for equally suggestive indications of a weak → strong allosteric transition under force.

      * Typos or writing remarks:<br /> - Fig 1B: The caption is inconsistent with the figure. In the caption, p₁ denotes the COM of Vt helices H2–H5, p₂ the COM of actin A1/A2 and p₃ the COM of actin A4/A5. On the figure instead (p₁, p₂, p₃) → (p₃, p₁, p₂).<br /> - Some repetitions that might be elegantly avoided<br /> x “capture the difficult-to-capture” (abstract)<br /> x “protein of interest [...] on our molecule of interest” (introduction §1)<br /> x “takes into account a Boltzmann weighted average over all possible configurations [...] a Boltzmann weighted average over all possible configurations” (introduction §4)<br /> - Typos<br /> x “FimH-manose” → “FimH-mannose” (introduction §8)<br /> x “which is the the direction” (3.1 §1)<br /> x “These approximate one-dimensional free energy pathways also give us a way to define when the system has crossed into the unfolded state” → Maybe the authors meant “unbound state” (3.1 §1)<br /> x “not all of the catch bond need come from” → “need to come from” (3.2 §3)<br /> x “we note that these results are suggestive, they do not quantitatively” → We wonder if the authors intended a formulation along the lines of: “we note that despite being suggestive, these results do not quantitatively” (3.2 §4)<br /> x “a constant prefactor the kinetic rate constant” → “a constant prefactor for the kinetic rate constant” (3.2 §4)<br /> x “grafted to our Holo structure random orientation” → “grafted to our Holo structure with a random orientation” (3.4 §2)<br /> x “TIP3 water” → “TIP3P water” (5.1 §2)

      Acknowledgements

      This report was written after a journal club given by OLC in the bussilab group meeting. All the members of the group, including external guests, are acknowledged for participating in the discussion and providing feedback that was useful to prepare this report. The corresponding authors of the original manuscript were consulted before posting this report.

    1. On 2024-11-14 15:48:57, user Mikel Garcia-Marcos wrote:

      Based on information brought up to the attention of the authors, the passage below will be modified in the final published version. More specifically, it has been clarified that the isoformof Rap1GAP used in the EMTA does not affect nucleotide exchange and that additional evidence supports that the this system detects active Galpha(i/o). These conclusions are partially supported by evidence presented in paper that has been published (doi: 10.1126/scisignal.adi4747).

      "For example, EMTA biosensors for G proteins of the Gi/o family are based on using Rap1GAP as a detector module, for which it is unclear whether it affects nucleotide exchange on different Gα subunits of this family or its preference for binding to Gα-GTP or Gα-GDP dissociated from Gβγ (35, 36), thereby raising questions about what is exactly represented by the BRET changes detected by this sensor."

    1. On 2024-11-14 08:45:38, user BioPioneer wrote:

      Another paper had solved this problem and was published about two years ago:

      Najafi, A., Jolly, M. K., & George, J. T. (2023). Population dynamics of EMT elucidates the timing and distribution of phenotypic intra-tumoral heterogeneity. Iscience, 26(7).

    1. On 2024-11-12 14:14:42, user Velyudhan Mohan Kumar wrote:

      In this report, human sleep-related gene orthologs have been studied in the unicellular green alga Chlamydomonas reinhardtii, which had been earlier considered as a model organism for the study of photosynthesis and flagella/cilia. Comparative genomics, which is an essential tool in biological research, is a very potent method that yields biological insights that are not possible with any other approach. It could also help us understand the human genome study, which has the potential to help us comprehend human sleep and sleep disorders. Phylogenomics gives a more profound comprehension of the discoveries made possible by comparative genomics.<br /> The results of this study indicate that the calcium signaling and synaptic transmission involved in sleep regulation are conserved across phyla, and a simple model organisms like Chlamydomonas may be useful in elucidating higher-order, complex behaviors such as sleep.

      Chlamydomonas do not possess the complexity of the sleep seen in vertebrates. However, knowing their genetic makeup can help researchers to understand general molecular processes that may be used to understand sleep in more complex organisms. Comparative phylogenomic studies are also important in understanding human sleep disorders, and they may open new avenues for developing new therapeutic approaches.

    1. On 2024-11-10 18:38:49, user Michele Bonus wrote:

      Hi! DataSAIL seems to be working well for me -- at least on protein structures using the built-in FoldSeek-based clustering.<br /> Regarding the paper: I believe there is something amiss with the coloring in Figure 2. In the lower sub-panel for the identity-based 1D cold split, shouldn't the rows be colored instead of the columns? Right now you are showing the same protein-identity-based split as in the upper sub-panel, and I believe you meant to show sub-optimal partitioning w.r.t. the "ligands" / chemical scaffolds.

    1. On 2024-11-06 08:57:18, user himanshukhandelia wrote:

      Very interesting work. I was curious about the surface tension being applied. For semi-isotropic pressure coupling, the tension is automatically set to zero, and does not need to be set explicitly. Why is this being done?

      We previously demonstrated that even ATP could bind to lipid membranes. 10.1021/acs.langmuir.9b01240. So does UMP (shown by both scattering, fluorescence and simulations (https://doi.org/10.1039/C9S....

      One caveat is that the charge on these nucleotide phosphates are most likely overestimated (see for example https://doi.org/10.1021/acs.jpcb.2c06077 and https://doi.org/10.1016/j.xcrp.2021.100343 . One way to address this is the rescale charges on the P and O atoms, as well as ions. However, a more systematic force field approach may be needed in the future. Pavel Jungwirth is working on this. https://www.developmentaid.org/organizations/awards/view/470320/doing-charges-right-modelling-ioncontrolled-biological-processes-with-the-correct-toolbox-q-scaling

    1. On 2024-11-06 01:23:49, user sun zhen wrote:

      The results of this paper are not credible.

      Based on the gel result, we can see that only introns are removed by oligo dT purification, lots of precursor and nicked RNA still in the products. Obviously, there will be a decrease in dsRNA compared with RNA kit purification which removed nothing from the crude circRNA products.

      When it comes to in vitro expression test, I am surprised to see even the expression of RNA kit circRNA is better than HPLC purified circRNA.We have also tested and compared RP-HPLC purified circRNA (circEGFP and circFLuc from GenScript) before, they got much higher expression than RNA kit purified circRNA or RNase R treatment crude circRNA.

      As the purity of circRNAs is very important for their expression (which can reduce their immunogenicity) based on previous reports, the in vitro or in vivo data in this paper are not convincing.

    1. On 2024-11-05 11:17:02, user Guei-Sheung Liu wrote:

      The paper has been published as below<br /> Proc Natl Acad Sci U S A. 2024 Nov 5;121(45):e2408345121. doi: 10.1073/pnas.2408345121. Epub 2024 Oct 30.<br /> Characterization of RNA editing and gene therapy with a compact CRISPR-Cas13 in the retina

    1. On 2024-11-04 15:12:37, user Jorge Hernández wrote:

      Dear authors, is there a way to access the supplementary data? I am looking for all the Marchantia TFs tried to express as HALO fusions.<br /> Thanks in advance!

    1. On 2024-11-03 15:04:26, user B. Hendrich wrote:

      I think you mean "potency". "Pluripotency" refers to the ability to make all cell types in the embryo. Trophectoderm stem cells cannot do that.

    1. On 2024-11-01 16:14:55, user Haider, Saqlain wrote:

      Dear Authors,<br /> Thank you for sharing this exciting work on the role of ABA and GI in recruiting CO at FT promoter to activate its transcription (and flowering). I appreciate the detailed analysis and valuable insights you've provided. I have a few questions and observations that may help clarify some aspects of the study:<br /> 1. Higher CO in aba1-6 mutants at ZT12<br /> In figure 2A, a slight increase in CO expression could be observed in aba1-6 mutants compared to WT plants though this expression declines at ZT16. Further, a higher CO expression is detected at ZT12 and ZT16 in SUC2::CO:CIT-aba1-6 compared with SUC2::CO:CIT (in WT background). It remains somewhat unclear as to why the lack of ABA is leading to a slight (yet) increase in CO transcription. It is noted that in 9 days old seedlings at ZT16 (figure 2B), there is a decline in CO expression yet aba1-6 mutants transformed with SUC2:CO:CIT in figure 2A have a higher CO at ZT16 than WT transformed with SUC2:CO:CIT. What could be the reason behind that. Finally, it will be valuable if authors could measure the CO-FT expression at ZT12 in 9 days old seedlings (like they did for 13 days old). <br /> 2. Higher Pol-II occupancy at CORE in WT<br /> In figure 4B, a higher enrichment/occupancy of RNA Pol II could be detected at TSS of FT in WT compared with SUC2::HA:CO aba1-6 and SUC2::HA:CO gi100. Yet, these lines flower much earlier than WT (figure 5B). How do authors explain this discrepancy as low Pol II occupancy would lead to a lower FT (and delayed flowering). Moreover, it will be interesting if authors could quantify and compare the FT expression in WT (Col-0), WT (with either SUC2::CO:CIT or SUC2::HA:CO), aba1-6 (with either SUC2::CO:CIT or SUC2::HA:CO), and gi-100 (with either SUC2::CO:CIT or SUC2::HA:CO). For example, do authors expect FT levels to be same in aba1-6, gi-100 (transformed with SUC2::HA:CO) as they flower nearly at the same time (Figure 5B). Moreover, as SUC2::HA:CO aba1-6 gi-100 show a slight delay in flowering compared with either SUC2::HA:CO aba1-6 or SUC2::HA:CO gi-100, it will be interesting to analyse FT transcript accumulation between genotypes. Taken together, addressing these points will increase the readers understanding and make the paper more concise.<br /> 3. ChIP and RNA Pol II occupancy assayed at ZT12 (not ZT16)<br /> In figure 3,4, and 5 the ChIP assays and RNA Pol II occupancy was quantified at ZT12. CO transcript peaks around that time point in long day (LD) conditions but the FT peak is observed at ZT16. Previously, it has been shown that reduced FT peak at ZT16 leads to delay in flowering ( https://doi.org/10.1111/pce.13557) . Authors mention CO binding is observed at ZT12 at CORE (line 329-331) but the delayed flowering phenotype is justified by FT levels so it will be interesting to check CO binding together with RNA Pol II enrichment at ZT16 and then compare both these assays (ZT 12 and ZT16). <br /> 4. RNA Pol II Occupancy at TSS in WT vs. SUC2::HA:CO abi1-6 and SUC2::HA:CO gi-100 with Delayed Flowering in WT <br /> SUC2::HA:CO: gi-100 and SUC2::HA:CO: abi1-6 lines show a drastic decline in Pol II occupancy at TSS of FT, yet these lines flower significantly earlier than WT. Would you consider discussing possible mechanisms by which this delay occurs, despite higher Pol II occupancy, to help clarify this interesting finding? Finally, the SUC2::CO:CIT gi-100 lines flower nearly at same time as SUC2::CO:CIT (figure 5B). Although in supplementary figure 9, an increase in bolting time (flowering in terms of number of days) in SUC2::CO:CIT gi-100 lines could be seen. A few additional details or a discussion of possible mechanisms might help contextualize this observation.

      Thank you again for sharing this preprint. I hope these comments are helpful and supportive of your ongoing work. Please feel free to reach out if you'd like further clarification on any of my points.<br /> Thanks<br /> Saqlain

    1. On 2024-11-01 14:02:04, user Cameron Thrash wrote:

      Hello. May I respectfully recommend a few relevant references for the discussion of SAR11 recombination and speciation?

      Natural variation in SAR11 marine bacterioplankton genomes inferred from metagenomic data<br /> https://biologydirect.biomedcentral.com/articles/10.1186/1745-6150-2-27

      High intraspecific recombination rate in a native population of Candidatus Pelagibacter ubique (SAR11)<br /> https://enviromicro-journals.onlinelibrary.wiley.com/doi/abs/10.1111/j.1462-2920.2007.01361.x

      The Evolutionary Success of the Marine Bacterium SAR11 Analyzed through a Metagenomic Perspective<br /> https://journals.asm.org/do...

      A comparison of homologous recombination rates in bacteria and archaea<br /> https://academic.oup.com/ismej/article/3/2/199/7588207?login=false

    1. On 2024-11-01 09:15:08, user David Žihala wrote:

      Dear Authors,

      Thank you for sharing your impressive work in the form of a preprint. I have a question regarding one aspect of the study that I found a bit confusing. Could you clarify why it was necessary to decontaminate the samples from mouse cells? If two samples were excluded from downstream analysis due to high mouse cell content (over 65%), could you please provide information on the mouse cell content in the remaining samples? Additionally, I am curious about how this type of contamination could have occurred.

      Best regards,<br /> David Zihala

    1. On 2024-10-31 15:54:51, user Josh Mitteldorf wrote:

      The comparison you have made supports the conclusion "that aging is not governed by a conserved universal program" but not the specific alternative you propose, "adaptations to damage and environmental conditions." I suggest that it may be that the ecological need for managing length of life is universal, but that it is implemented differently in different species. Salmon destroy their bodies with glucocorticoids; elephants fail to grow a 7th set of teeth and can't chew their food; mice get cancer; rats get heart disease -- but the underlying adaptive motivations might be universal. The best evidence for this is the role of insulin signaling, which accelerates aging across the animal kingdom.

    1. On 2024-10-31 15:40:16, user James Amos-Landgraf wrote:

      In your methods you state that you purchased C57BL/J mice from the four vendors. These mice are only available from Jackson Lab. The other mice you purchased are C57BL/6N substrains that specific to the vendor. NTac, NCrl, and NHsd. The biggest factor you need to address is the B6J have a mutation in the Nnt (nicotinamide nucleotide transhydrogenase) gene that dramatically influences metabolism. You should address this in your interpretation of bacterial community shifts. There is a good review of these differences that Jax has published on their web site and in various publications. https://www.genengnews.com/sponsored/heres-what-you-need-to-know-about-the-c57bl-6-substrains/

    1. On 2024-10-31 14:51:48, user YC Han wrote:

      This study holds considerable promise, but its findings are significantly undermined by the notable delay in developmental stages observed in embryos injected with CRISPRi/a + sgRNAs, as evident in Figures 1-3. For instance, the embryo injected with CRISPRi + tyr sgRNAs marked with 48 hpf in Figure 1B, appears to be developmentally equivalent to a 30 hpf wild-type embryo, which exhibits substantially fewer pigment cells than embryos at 48 hpf. This discrepancy raises concerns about the direct impact of the CRISPRi/a + sgRNAs injection, whether or not the delay development was caused by an unintended consequence.<br /> I strongly recommend that the authors reexamine their results and conduct a comparative analysis of CRISPRi/a + sgRNAs-injected and control embryos at equivalent conventional developmental stages, as described by Kimmel et al (1995, PMID: 8589427). This crucial consideration will not only enhance the study's rigor but also provide valuable insights to the zebrafish research community.

    1. On 2024-10-31 14:41:16, user CT wrote:

      It’s been known for decades that Blue Jay kents and IBWO kents differ spectrographically… duhhh. And a LOT of other sounds may mimic IBWO kents to the human ear as well (though the paper explores but one single avian species)… all being spectrographically different. This paper adds little to the known literature, except to try to reach further conclusions (actually redundant, re-stated conclusions in previous work) from unsubstantiated assumptions. <br /> I’d be amazed if one can find 2 or 3 (or even one) PhD. level scientists in the entire country who will defend the methods and exaggerated conclusions put forth here.

    1. On 2024-10-30 14:15:51, user Anonymous wrote:

      Dear authors,

      as part of a group activity in our lab we discussed your very interesting paper with the goal to review it. The below review is the result of this exercise and therefor reflects the thoughts and concerns of several people. We hope this helps you with your way forward to publish the paper in a good journal.

      Review:

      The manuscript by De Tito et al. reports a hitherto unknown role for ATG9 in recruitment of PI4K2A to damaged lysosomes to control the levels of PI(4)P, lysosomal membrane contact site formation with the ER, and lipid transfer from ER to damaged lysosomes. They report that this mechanism is controlled by ARFIP2, which shuttles between the Golgi where it anchors ATG9 and lysosomes, and by the AP-3 complex, which mediates retrieval of ATG9 from lysosomes. Finally, the authors show a role for their proposed mechanism in lysosome damage induced by Salmonella infection. The manuscript deals with the timely topic of lysosomal damage repair and PI4K2A was recently emerging as a key player in this process. Mechanisms for PI4K2A recruitment to the sites of damaged lysosomes were so far elusive. Although the manuscript is certainly of interest for the field, the presented evidence on which the claims and the proposed model are based on seem not always substantial enough and further work is required.

      Major concerns<br /> • Novelty and generalization: Although the findings suggest new functions for ATG9A and ARFIP2, they overlap significantly with previous work on autophagic and lysosomal repair pathways. The novelty of these findings, in comparison to earlier research, could be emphasized more clearly.<br /> • The claim that ARFIP2 "modulates lipid transfer for lysosomal repair" could benefit from additional direct evidence linking lipid transfer to lysosomal recovery and ARFIP2's specific role in this process. There is only one in vitro experiment directly showing lipid transfer modulation by ARFIP2 (Fig. 5K). The authors should use the ARFIP2 W99A mutant in this experiment to test whether lipid transfer modulation is specific and according to their model.<br /> • A central claim of the manuscript is a role for ARFIP2 in the repair of damaged lysosomes. The gold standard assay in the field is recovery of lysotracker fluorescence after LLOME-induced damage as shown in Fig. 2G. The way the experiment is presented does not instill confidence. The effects are relatively modest, it is unclear how the statistics were done, the representative images in the supplement only remotely resemble lysotracker stainings, how a reliable lysosome number from these images could be extracted is unclear, it is not shown whether and how much the cells express GFP-ARFIP2, the time points of the images do not match the relevant time points in the quantification, and the lysosome number in the different samples before LLOME treatment is not factored in. <br /> • The major lysosome repair assay in the manuscript is Galectin staining. This is a rather indirect assay as compared to the lysotracker recovery assay as it shows the damage rather than the repair. Could there be repaired lysosomes that still are positive for Galectin? Is the total Galectin expression level the same in the cell lines and conditions used? Furthermore, spot quantification for lysosomal damage markers like LGALS3 (Figure 2f) does not account for cell size or density, potentially leading to misinterpretation of the data. The authors further need to show lysotracker recovery in ATG9 loss of function to substantiate their claims.<br /> • The authors show that ARFIP2 and ATG9 affect PI4K2A localization in cells. The “delivery” of PI4K2A to lysosomes, however, which is a central claim in this manuscript, is insufficiently demonstrated. Same is true for the “retrieval” of ATG9: The authors claim that AP-3 and ARFIP2 are important for the retrieval of ATG9A vesicles. Even though they show effects of these proteins on ATG9A presence at the lysosomes, they never manage to show that retrieval really is impaired. I am not fully sure if it is possible to make videos that convincingly show this. Either they should try that, or they should not hammer so hard on the word retrieval when they never show it specifically. <br /> • The dynamic interactions of the relevant components need to be better demonstrated and characterized. E.g. interactions ATG9-PI4K2A interaction needs to be proven and characterized better (e.g. by co-IP, immunogold or similar techniques). Otherwise the role of ATG9 is much secondary, and should be a ARFIP2-mainly focused paper. The experiments in figure 4H and 4G seem ideally suited for this purpose but the authors fail to show the relevant components in the same blot.<br /> • The increased formation of membrane contact sites is insufficiently demonstrated. The authors need to use electron microscopy, super resolution microscopy, SPLICS or other appropriate techniques to make this claim<br /> • The way of quantifying protein localization to lysosomes (intensity lysosomes/total) seems heavily dependent on the lysosome coverage of the cells. In conditions where to my impression the overall number of lysosomes also changes, I would like to see a negative control that demonstrates that the enhanced lysosome localization is not just by chance, since there are more lysosomes. This applies for instance to Figure 4E+F and 5A+B.<br /> • I am not satisfied with the data analysis they presented. The authors should indicate the statistical test for each piece of data they analysed including an explanation, why did they chose particular test. For instance, in fig. 2g they compare the amount of puncta or lysosomes (it is not clear as well) in the cell in a pairwise manner. It would be more appropriate to implement a statistical test that can compare the curves fitted with the data points or a test which can compare each time point individually, i.e. Kaplan-Meier plot.

      Minor comments<br /> • The zooms generally miss scalebars<br /> • Which part of the picture the zoom is from is often not indicated<br /> • Many of the images are hard to interpret because of low contrast, I would recommend to put their individual channels in grey<br /> • Fig3f. requires a better representative image and although there is a line indicated there is no line plot<br /> • Fig. 1c: show quantification for pS6K<br /> • Fig. 2b: data points are the same as figure 1b<br /> • Fig. 3: Live-imaging of ATG9A with AP-3D1 is needed

    1. On 2024-10-30 08:51:15, user Thomas Munro wrote:

      Given LM189's structural similarity to salmeterol and higher efficacy, I suspect there would be wide interest in whether it shares salmeterol's peculiar time course: activation that persists for hours after washout, but can be reversed by antagonists, and “reasserts” after antagonist washout. The cited patents claim “a particularly long duration of action”, which suggests LM189 may be even more extreme than salmeterol, but don't appear to give any results.<br /> It might be possible to get this data without the dreaded reviewer experiments. Glaxo (now GSK) offers to share clinical trial data. Presumably they would also be willing to share in vitro data for an off-patent compound like this, which has no privacy or IP concerns.

    1. On 2024-10-29 14:01:42, user Amandine Véber wrote:

      In this bioRxiv preprint, the authors question part of the methodology that we have used in Jay et al. (2022) , as well as the conclusions we have drawn from our study. In order to clarify the different points of misunderstanding that have led to these criticisms, in the document available at this link we provide a response which, we hope, will demonstrate the validity of our approach.

      In this document, we first address the main questions raised by our colleagues and then carefully go through our paper to emphasize its main contributions. The format of this response makes it not eligible for publication in standard open archives, hence our non-standard way of making it available online. A link to it can also be found there .

      The current debate which was generated by Jay et al. (2022) is grounded on the very interesting question of which control scenario is appropriate to assess the efficiency of the "sheltering effect" that we investigate there, and we hope to contribute to this discussion with sound logical arguments and relevant biological concepts.

      Amandine Véber, Paul Jay and Tatiana Giraud

    1. On 2024-10-29 04:24:41, user Dan MacNulty wrote:

      The age estimates are impressive. To help readers appreciate the significance of these estimates, consider adding more background about what’s currently known about the age of aspen clones in general. Do estimates exist for other clones? If so, how do they compare to your estimates? How does your aging method compare to previous methods?

    1. On 2024-10-28 15:02:24, user nuala oleary wrote:

      Hi,

      I wanted to point out a small issue with the references to NCBI Datasets. The author names are misformatted. It should read as follows:

      O’Leary NA, Cox E, Holmes JB, Anderson WR, Falk R, Hem V, Tsuchiya MTN, Schuler GD, Zhang X, Torcivia J, Ketter A, Breen L, Cothran J, Bajwa H, Tinne J, Meric PA, Hlavina W, Schneider VA. Exploring and retrieving sequence and metadata for species across the tree of life with NCBI Datasets. Sci Data. 2024 Jul 5;11(1):732. doi: 10.1038/s41597-024-03571-y.

      Thanks!

    1. On 2024-10-28 09:39:36, user Isabella Capellini wrote:

      A revised version of this manuscript is now available in Proceedings B:<br /> Mortlock E, Silovský V, Güldenpfennig J, Faltusová M, Olejarz A, Börger L, Ježek M, Jennings DJ, Capellini I. 2024 Sleep in the wild: the importance of individual effects and environmental conditions<br /> on sleep behaviour in wild boar. Proc. R. Soc. B 291: 20232115.<br /> https://doi.org/10.1098/rspb.2023.2115

    1. On 2024-10-24 16:51:29, user Allisandra Rha wrote:

      Interesting work. The multimodal approach definitely enhanced the integration outcomes and is well-designed. While it is of benefit that AAV efficiently reaches the nucleus, moving forward would you consider evaluating other delivery methods for nuclear translocation to account for patients that are seropositive for AAV antibodies?

      Also, the paper is under-referenced. I would encourage you to look at the text and provide more citations to support your claims.

    1. On 2024-10-23 00:20:34, user CDSL JHSPH wrote:

      I really enjoyed reading your preprint "Universal rules govern plasmid copy number." Your analysis of plasmid copy number across diverse bacterial genera is impressive and the discovery of a scaling law linking plasmid size and PCN is a huge contribution to the understanding of plasmid biology. I was wondering, based on your discussion on the intrinsic variability of PCN in high-copy number plasmids versus low-copy number plasmids, how these variations impact the fitness of bacterial populations, especially under different selective pressures. Additionally, given plasmids' role in gene transfer, how might the scaling law and PCN-size trade-off influence the evolution of plasmid-encoded traits, such as antibiotic resistance genes or virulence factors? <br /> I look forward to seeing how this research progresses! Thank you!

    2. On 2024-10-19 00:01:15, user CDSL JHSPH wrote:

      Thank you for sharing your research! I found your paper on PCN very insightful, especially the discovery of "universal rules" governing PCN. The large-scale analysis of plasmids provides valuable insights into plasmid number variability and regulation, with potential applications in both research and biotechnology. I do have a few questions. Could replicon dominance bring evolutionary advantage, particularly under natural conditions where plasmids may compete for replication resource? Additionally, under external constraints, is there a threshold where dominance shifts between replicons? <br /> I very much enjoyed the article, and look forward to how these findings evolve in future studies.

    1. On 2024-10-23 00:19:08, user CDSL JHSPH wrote:

      I really enjoyed learning about this new Bento Lab because it opens the doors for so many new positive things. I think the paper is written in a way that people will understand the differences between BL and TL and agree that there are so many advantages to the BL. I think your strongest part is how you explain your data from the figures since you talk about the how and identify how to understand the figures and tables. I would say for figure one; I was wondering if maybe also using bar graphs since it shows a clearer way of differentiation like Figure 2. I also think instead of tables, more figures will help the public understand their pros and cons better visually. I also think that the sludge's low read count should be addressed more. While it is mentioned that the Flongle's limitations might have affected the detection of less abundant entities, why do you think happened to that sample, and how could you express that to the audience? I think you could strengthen their argument by referencing specific species, ARGs, or plasmid types they suspect might have been missed due to the lower read count. To end this, I think the authors did a great job explaining the differences, and my only concerns were adding figures and explaining the sludge sample.

    2. On 2024-10-19 18:05:46, user CDSL wrote:

      This article is a good demonstration of the potential of portable DNA sequencing technology for rapid detection of pathogens and antimicrobial resistance in the field, especially for public health emergencies in low-resource or remote areas, and the results section, in particular, is very detailed. However, I feel that there are some deficiencies in the discussion part. It is mentioned in both the introduction and the results that the differences are evaluated from five aspects, but I do not seem to see obvious discussion about DNA yield and purity in the discussion part. Also, have you considered separating the restrictions into a separate section? That may help the reader to read in a more organized way.

    1. On 2024-09-10 18:32:52, user Thomas Sorger wrote:

      Please note that the date of the second reference (2003) is incorrect. The correct date is 1983:<br /> 2. Armstrong E. Relative brain size and metabolism in mammals. Science 220: 1302–1304 (1983).<br /> Tom Sorger

    1. On 2024-10-22 19:01:27, user CDSL JHSPH wrote:

      I believe this is an excellent study and paper, with valuable contributions to the field. It was very well-written, and the methodology very clearly described and reproducible. Given the global clinical significance of malaria, the biological implications of this study can be used as a basis for further studies into the biology of circadian rhythms in malaria transmission, as well as clinical translation for efficient prevention and control strategies and more effective treatments.

      Similar to what previous reviewers have noted, I noticed some inconsistencies between some of the figures and the apparently corresponding results in the text, especially as it relates to the percentages of the transcriptome displaying cyclic expression profiles, for both the mosquito salivary glands and sporozoites. The figures and/or results should therefore be corrected accordingly.

      A possible limitation to the study, which I don't believe was addressed in the paper, was that different mosquito (Anopheles stephensi vs Anopheles gambiae) and parasite (Plasmodium berghei versus Plasmodium falciparum, vivax, ovale, malariae, or knowlesi) species from the ones which are primarily known to cause malaria in humans, were used for analysis in this study. This may limit direct clinical translation of these results for malaria in humans. There should be a statement addressing this in the discussion. Furthermore, the amount of gene homology between the species used in the study and the species causing disease in humans should be stated for reference. If there were any reasons it was not possible or necessary to use the species known to be most associated with disease in humans, this should be stated.

    2. On 2024-10-22 01:49:09, user CDSL JHSPH wrote:

      I thoroughly enjoyed reading your paper, and I believe this study brings valuable new insights into how circadian rhythms influence malaria transmission. This information will undoubtedly impact future strategies for combating malaria, but also other vector-borne diseases. However, I noticed some inconsistencies between the results and the corresponding figures. I might be wrong and misread it but I wanted to bring it to your attention and clarify it on my side as well. In certain sections, figures were cited but did not seem to align with the content discussed. I can mention for example:<br /> - “We found that 27-49% of salivary gland genes displayed a cyclic expression profile, (Fig. 1B-E) representing 5-10x more than what has been previously reported for the mosquito’s head or body”. For C-E, I don’t see how this is actually depicted. Could you please explain; even the next sentence as well I don’t understand how C-E are supposed to back the arguments you are trying to make <br /> - Some figures like 1G and 1H are mentioned in the legend of your figures and in the text as well but are not shown in the figures, why?<br /> Those are not the only inconsistencies I noticed, but before proceeding on citing them, I want to make sure that I am not misreading and not understanding.

      If my assumptions are true, clarifying them would greatly improve the paper’s clarity and presentation of your significant findings. I am looking forward to reading your answer.

      Best regards,<br /> Critical Dissection Class 2024/2025

    3. On 2024-10-19 02:11:16, user Yak Nak wrote:

      The manuscript provides an exciting and valuable look into how circadian rhythms influence malaria transmission by aligning mosquito feeding behavior and parasite activity. The use of RNA-sequencing to uncover rhythmic gene expression in mosquito salivary glands is a significant strength and offers important new insights into the mechanisms behind malaria transmission. The figures are clear and effectively illustrate how these rhythms correlate with mosquito feeding efficiency and parasite infection capabilities, though the figure legends could benefit from more detailed explanations, especially for readers unfamiliar with gene expression data. The introduction is solid but could be improved by providing a more detailed discussion of previous research on circadian rhythms in malaria parasites to better frame the novelty of this study. The discussion section does a good job connecting the findings to broader vector-borne diseases like Zika and dengue, but it would be even stronger with specific examples of how these results could inform practical malaria control strategies, such as optimizing the timing of interventions based on mosquito feeding times. Overall, this is a well-conducted study with important findings, and a few revisions could further enhance its clarity and impact. Two follow-up questions:

      1. Could the authors clarify why specific time intervals (every 4 hours) were chosen for RNA-sequencing, and would more frequent sampling provide additional insights?
      2. Also, how might environmental factors such as temperature or humidity influence these circadian rhythms, and could this affect transmission in different regions?
    4. On 2024-10-18 18:57:00, user CDSL JHSPH wrote:

      I really like this article. Thank you for sharing your article. Malaria remains a major global health challenge, and understanding the biological rhythms of mosquito vectors and malarial parasites is crucial for improving control strategies. However, previous studies have only revealed that mosquitoes have daily rhythmic behaviors. You provided new insights. First, about half of the genes in the mosquito salivary gland transcriptome show 24-hour rhythmic expression, and second, the gene expression of sporozoites in the salivary glands also shows circadian rhythms (parasite movement and infection ability). In addition, you mentioned in the discussion that the parasites and mosquitoes have evolved in coordination with the circadian rhythms, and jointly affect host infection by regulating parasite movement and mosquito blood feeding time. This is very interesting, and your findings provide a new perspective for understanding the temporal regulation mechanism of malaria transmission. At the same time, you mentioned optimizing the insecticide spraying time strategy according to the peak period of mosquito activity. In addition, your findings may also have reference effects on other diseases.

      But I have some questions. First, it is good that you use 12h dark and 12h light to simulate day and night. But in the wild, the temperature and humidity vary between day and night. Perhaps the natural day-night cycle could be better simulated by changing temperature and humidity. Secondly, you mentioned that some observations may apply to uninfected mosquitoes, but you did not specifically discuss the potential differences in rhythms and behaviors between uninfected and infected mosquitoes.

      Finally, thank you again for your contribution and new ideas.

    1. On 2024-10-19 19:13:36, user sam wrote:

      It would be good to see how altered the phosphoproteome is after FACS. Thats a long time for the cells and a lot can change in phosphopatterns in a matter of seconds .

    1. On 2024-10-19 00:19:59, user CDSL JHSPH wrote:

      This article is rigorous. I like the comparison between phages in gut and that in the sea water. And the analysis of phages in hCom2 mice is well designed. Yet I think the equations and the explanations in the result and the method part are not very readable. I’m confused about the meaning of lysogen-host ratio and relative coverage, the role of dilution rate and so on. Plus, figure 2B is not very informative and by just looking at it, I cannot have a clue of what you were doing to get the range of induction rate. But overall your results support your conclusions firmly and clearly.

    1. On 2024-10-18 13:39:11, user Oscar wrote:

      Dear James and co-authors,

      I trust you are already working on submitting this manuscript to a reputable journal. As one of the researchers with the most publications on CAV1 in Ewing Sarcoma (Oscar M. Tirado), I would like to offer some comments and suggestions for your consideration:

      It might be beneficial to mention that CAV1 is a demonstrated direct target of EWS::FLI1, as this could add further context to your findings.

      Regarding the presence of caveolae in Ewing Sarcoma cells, I assume you've reviewed the manuscript (DOI: 10.1016/j.canlet.2016.11.020). In most of the cells studied, including TC71, the levels of Cavin-1, as well as caveolae, are extremely low. As such, the majority of CAV1’s functions in EwS are likely independent of caveolae, which aligns with the results you’ve presented.

      I noticed the focus on the AKT pathway in your work. Have you considered discussing the MAPK pathway as well? It is also affected and likely plays a significant role in the processes you're describing.

      In Figure 4, CAV1 emerges as a potential driver, along with EphA2. I’m sure you are familiar with this study as well (DOI: 10.1002/ijc.31405), which may provide further insights.

      I agree with your assessment that CAV1 plays a key role in the progression of EwS, and your findings align with many of the conclusions we've drawn in our own research. I would be happy to assist further, either through this platform or via personal contact, as we've known each other for many years from various EwS meetings.

      Best regards,<br /> Oscar M. Tirado

    1. On 2024-10-18 01:51:53, user Tushar R. wrote:

      Summary

      Here, the authors sought to apply cryo-EM guided metainterference-based MD simulations to find modeling inaccuracies of flexible helical regions linked to their artificial representation by a single-structure model. They first applied their approach to the group II intron ribozyme from Thermosynechococcus elongatus, one of a few cryo-EM structures available in the PDB that featured mostly RNA and was obtained using a single-structure refinement procedure. They confirmed that remodeling with their approach mostly affected flexible regions at the solvent-exposed stem loops that were not phylogenetically conserved. They found that functional domains that were well-ordered, phylogenetically conserved, and were clearly represented by the cryo-EM density required no remodeling. Extending this analysis through other PDB structures revealed that poor modeling at flexible helical regions was broadly applicable to all RNA-containing cryo-EM-derived structures in the 2.5 - 4 Å resolution range.

      Major Points:

      It is difficult to keep track of the different parameters applied for each simulation. These are scattered in the text. It would be very helpful for the reader to understand these if the authors could make a table of all the MD simulations with specifications on helices restrained, approach, simulation time, force fields present, equilibration time etc.

      Clarification regarding the helical restraints and the simulation time for the initial production run would be helpful—i.e. Why was a simulation time of 2.5 ns chosen, and would this be long enough for your purposes considering the relatively high complexity of the system?

      Were restraints applied to the 3 helices (b,d,i) that unfolded in subsequent simulations or were they only applied to the other 6 helices?

      In the section “Base Pairing Analysis of the Protein Data Bank,” it is mentioned that the analysis likely over-estimates the problematic modeling of helices in cryo-EM derived RNA structures. Given that 15 Å is larger than the expected 8-11 Å (Pietal, M. J., 2012) distance for N1/N9 distance, would this not underestimate the number of problematic helices?

      A detailed analysis of H-bonding within the 6 helices would be useful to get a mechanistic understanding of why the 3 helices (b,d,i) unfold and the other 6 don't—for example, it might be helpful to provide a comparison with the model generated from a single structure approach to know if the same H-bonds exist in it or not as the model forwarded by this paper.

      ERMSD approach: we know relatively less about non-Watson Crick base pairing apart from Hoogsteen and G-U wobble base pairs since so much depends on the context (specific region of the structure, ion concentration, pH etc.) (we’re still thinking about how to phrase this point).<br /> Balance between experimental measurements and MD sampled conformational states: lack of experimental validation when remodeling mismodeled regions for RNA - how do we know when specific interactions are too ideal?

      How is variability between forcefield parameters for divalent metal ions or other ions addressed, especially considering that these are required for proper folding of RNA?

      Dispersive interactions play a huge role in RNA integrity: how confident can we be about these parameters when comparing various nucleic acid specific force fields? Can these differences lead to unfolding of the 3 properly base paired helices?

      With reference to the line on page 4: “The trajectories obtained with metainference simulations were then analyzed by back-calculating the corresponding averaged density map and comparing it to the experimental one”, was the solvent density included in the back-calculated density?

      What proportion of the CC_mask arises just due to fitting of rigid part of the structure and how much of it arises due to the improved fitting of the flexible regions from the 9 loops? Would it be possible to separate these contributions?

      Minor Points:

      Would it be possible to show Fig 2C as a violin or box plot rather than a bar plot? This would allow visualization of the distribution of CC_mask for different conditions with clusters representing conformational states that may agree with the experimental data.

      Is local resolution considered while plotting the RMSF? For this purpose, it would be useful to have a local resolution estimation for the map to help the reader understand whether the unfolding of helices occurs simply because the specific regions were not well defined or if these regions are inherently flexible.

      How effective is DeepFoldRNA in filling gaps in structure as it is normally known for sequence based structure prediction—specifically regarding the modeling of the 38-nucleotide gap?

      For benchmarking purposes, applying this approach to other RNA structures or providing additional validation against independent experimental data (e.g., SAXS or chemical shifts from short RNA motifs) could further strengthen the conclusions.

      Aside from the pre-1r (6ME0) and pre-2r (6MEC) states, Haack et al. (2019) found additional 3D classes that indicated disordered density in conserved regions, which they did not investigate further. They also suggested that the post-2r complex may be captured in one of the 3D classes that yielded a low-resolution 3D reconstruction. Is it possible that one of the structures found by the metainference method could have sampled one of the structures that were not investigated further or the post-2r complex?

      Why is there a sudden fall in the number of helices formed in 32 replicas (Fig 2E)?

      Sentence on page 4: “As a result, 32 replicas appeared to be the best compromise between agreement with experiment and computational cost.” Has it been ensured that this is not a result of overfitting? One of the ways to do this can be to use half-map cross validation i.e. to check if the refined ensemble fits both the half maps equally.

      Have the authors tried to change the weight for 1 μs reference simulation to improve the CC_mask?

      Grammar/spelling mistakes <br /> Results, Test System and Preparation, 1st Paragraph: Change “run” to “ran” to retain past test throughout the paragraph in the sentence beginning “We then run 2.5 ns-long molecular dynamics (MD) simulations in explicit solvent…”<br /> Results, Ensemble Refinement, 1st Paragraph: “Within” is misspelled as “whitin” in the sentence beginning “Conversely, when restraining their helicity, whitin this single-replica refinement approach…”<br /> Results, Ensemble Refinement, 1st Paragraph: “Specifically, all the simulations attempted with a lower number of replicas were crashing reporting missing convergence in enforcing bond constraints, which indicates that the experimental and helical restraints were mutually incompatible.”

      • Tushar Raskar, Sonya Lee, James Fraser
    1. On 2024-10-17 14:05:57, user Christina Warinner wrote:

      It is very nice to see the authors statistically confirm on a large number of samples that oral bacteria contribute to the thanatomicrobiome of archaeological teeth. They may want to note that this pattern has been previously observed and reported twice before: <br /> Mann AE, Sabin S, Ziesemer KA, Vågene Å, Schroeder H, Ozga A, Sankaranarayanan K, Hofman CA, Fellows-Yates J, Salazar Garcia D, Frohlich B, Aldenderfer M, Hoogland M, Read C, Krause J, Hofman C, Bos K, Warinner C. (2018) Differential preservation of endogenous human and microbial DNA in dental calculus and dentin. Scientific Reports 8:9822. DOI: 10.1038/s41598-018-28091-9<br /> Vågene AJ, Campana MG, Robles García N, Warinner C, Spyrou MA, Andrades Valtueña A, Huson D, Tuross N, Herbig A, Bos KI, Krause J. (2018) Salmonella enterica genomes recovered from victims of a major 16th century epidemic in Mexico. Nature Ecology and Evolution 1-9. DOI:10.1038/s41559-017-0446-6.

    1. On 2024-10-15 20:51:50, user Alexa Jennings wrote:

      Overall, this paper is significant to developmental biology. Thank you for your contribution and dedication to the field. As I read through your paper, I have a few comments I would like to offer.

      Fig1: While I was intrigued by the several stainings and experiments carried out by your team, it would be helpful to understand your rationale for choosing the gestational ages of interest. Additionally, since this is a comparative review, it would be useful to visualize the corresponding gestational ages in mice. I also noticed you cited these phases in mice, but it would be much more memorable if there was a clear schematic comparing the two species' gestational ages to go along with the rest of the figure as well.

      Fig2: I noticed a one-way ANOVA was carried out for several of the corresponding data. However, ANOVA's are designed to compare means of groups with samples >3. In the supplemental figures, several groups were either missing data, or had samples <3. Asterisks marking significance would also be useful to visualize significance on the graphs.

      Fig3&4: Quantification of fluorescence would be useful to make statistical comparisons. As is, the evidence is too correlative to make any conclusions.

      Fig5: For clarity, place comparative groups in alphabetical order (ex: 5a and 5c should be 5a and 5b). Additionally, 5k-q show dim fluorescence, and should be quantified to make accurate statistical comparisons. 5f should be compared via a multivariate-ANOVA.

      Final remarks: Further manipulative, molecular experiments are required to make conclusions on similarity. A stronger argument could be made if additional evidence were included. Additionally, sex and cultural differences likely will cause variation among your data. Stating the sex of the samples can help make more accurate comparisons between data. Comparative experiments between sexes and across ethnicities in humans may also be useful.

    1. On 2024-10-14 05:52:34, user PanosMoschou wrote:

      Information missing from the Fig. legend of Fig. S1 panel I, in this version (data not present in the published version of the article in Plos Biology) represents a reconstructed image that fully reproduces the original dataset (in particular, the mating controls DDO and QDO experiments). If you wish to have access to the original dataset, don't hesitate to get in touch with the corresponding author. Furthermore, the group will soon post the original dataset in Zenodo.

    1. On 2024-10-11 22:57:08, user aquape wrote:

      Recent information, google<br /> "GondwanaTalks Verhaegen English"<br /> + see <br /> - M.Vaneechoutte cs 2024 Nat.Anthrop.2,10007 open access https://www.sciepublish.com/article/pii/94 <br /> "Have we been barking up the wrong ancestral tree? Australopithecines are probably not our ancestors” <br /> - id.2024 Nat.Anthrop.2,10008 https://www.sciepublish.com/article/pii/187 <br /> "Reply to Sarmiento E. “Australopithecine Taxonomy and Phylogeny and the Savanna Hypothesis…”"<br /> https://www.gondwanatalks.com/l/the-waterside-hypothesis-wading-led-to-upright-walking-in-early-humans/

    1. On 2024-10-11 09:50:17, user bao zhang wrote:

      I tried to run the script you uploaded on github, and the model.py script reported an error, "NameError: name 'p_SdC' is not defined", how is this p_SdC defined? I am looking forward to your reply, thank you.

    1. On 2024-10-11 08:52:29, user Esther Broset wrote:

      This article is highly informative, and the ionizable lipids show great promise. I was particularly impressed with the targeting strategy and the high expression levels observed in the lungs.

    1. On 2024-09-30 08:04:04, user Ema Nymton wrote:

      Zach's reply does not really address the points in my previous comment. I note that it also does not address the points raised in Bloom 2024.

      It only tangentially addresses 2 points in my post, both of which are apparently misconstrued.<br /> The first point that was addressed was about the p-values- however this response did not acurately portray the argument.<br /> My post never said it was misleading, but rather that the sampling bias towards the stalls would concentrate the highest p-values near the stalls, even more so than a simple relative-risk heatmap would, making the result of elevated p-values near the stalls expected whether or not the stalls were the origin of the outbreak.<br /> Compare: https://i.imgur.com/iDD6D93.jpeg to https://i.imgur.com/2YmYlZp.png

      The second point "an argument that implicitly assumes all expectations for all coronaviruses are identical regardless of their hosts and modes of transmission,"<br /> No, it does not assume they are "identical", but rather that there are general similarities. I am asking for an explanation why the distribution is so different.<br /> Neither mode of transmission nor host explain this - nor does date of sampling.<br /> While not explicitly stated, there is a substantial focus on Raccoon dogs.<br /> You show the distribution of a virus with the same putative host (Raccoon dogs) and same transmission methods (respiratory). The distributions are strikingly different.<br /> Compare https://i.imgur.com/iDD6D93.jpeg to https://i.imgur.com/OFXpCIf.jpeg

      This paper argued (in preprint), and it still argues )in final published form) that the 12 Jan RNA had more time to decay. Considering the market was closed at the same time, why did the SARS-CoV-2 reads decay to such lower values than other CoVs in samples from the same dates, and within the same samples, as shown by Bloom?<br /> See here: https://i.imgur.com/kIOTTYq.jpeg <br /> -even when stratifying by collection date, clear negative correlations still come up.

      Overall, I find that a lot of speculation is offered, but it comes accross as excuses for the data not actually showing that the virus originated from the wildlife stalls.

      I'll also raise additional critiques now:<br /> 1) The heat maps do not colocalize with the potential hosts:<br /> SARS-CoV-2 (RR p-values) and Raccoon dog reads shown here: https://i.imgur.com/AyAiag9.jpeg

      In the paper, instead the focus is on the south-west portion of the market.<br /> As I see it, there are two possibilities:<br /> i) The analysis DOES have sufficient resolution: in which case the resulting heat map not co-localizing with the wildlife stalls would exonerate them<br /> ii) The analysis DOES NOT have sufficient resolution: in which case it must be noted that the eastern and northern parts of the market were very poorly sampled. In this case, finding more positive samples where more samples were taken does not allow any conclusions to be drawn.

      2) The analysis of the SNPs of the Raccoon dogs clearly aligns with the the SNPs of locally caught wildlife (the C14859A + A15304G genome from Hubei, and the C14372T, C15102A, C15252A, T15306C genome from Hubei). Despite the evidence clearly pointing to locally caught raccoon dogs, speculation is offered (again) that maybe racoon dogs with these SNPs are found further south.<br /> Notably, Hubei (the local province) is not a plausible origin of the progenitor. Pekar seemingly agrees: https://www.biorxiv.org/content/10.1101/2023.07.12.548617v1

      In essence, this paper shows:<br /> i) 3 lineage B sequences 1 low quality lineage A sequence have a (very very wide) tMRCA confidence interval that overlaps with a tMRCA calculated from many more A and B lineage genomes: hardly anything surprising or something that conclusions can be drawn from

      ii) Within the heavily sampled area of the market, the areas of elevated relative risk don't actually overlap with the locations of the wildlife stalls<br /> -There is just one stall on the periphery of an area of elevated risk, the one closest to the entrance/mah jong rooms/bathrooms, which is where the higest elevated risk is.

      iii) The raccoon dog genomes from that stall suggest that they are locally caught, and do not come from areas with potential reservoirs of SARS-CoV-2's precursors

      On balance, the sum of this evidence seems to suggest that the HSM wildlife stalls in general, and raccoon dogs specifically, were not the origin of this virus.

      Much speculation and many excuses are offered, but it must be reiterated that these excuses and speculation, not evidence.

      The correlation analysis of Bloom can point right to the hosts of other viruses. The heat maps produced for this paper can point right to the hosts of other viruses. SARS-CoV-2 continuously fails to produce evidence and associations that other animal viruses at the market produced.

      The way these findings are presented in the paper and to the general public through the media seems to be at odds with what the data actually shows.

    2. On 2024-09-21 04:17:05, user Zach Hensel wrote:

      I have a short response to the two comments on this preprint, which, of course, we took into consideration while revising the manuscript, which is now published following peer review here: https://doi.org/10.1016/j.cell.2024.08.010

      One commenter, David Bahry, has taken to social media to call myself and co-authors "frauds" who are "trying to pretend" and made some vulgar comments that can't be repeated in this comments section because, he says, we "ignored" his comments. This is not true, so it's appropriate to respond here.

      Both commenters note that low p-values for relative risk maps (Figure 2 and Figure 4) require sufficient sampling density to obtain a low p-value. Of course, this is correct. Both commenters argue that this is misleading. I disagree. No co-author, peer reviewer, or editor involved was misled. News articles on the paper thus far have almost all portrayed the data and analysis accurately. We do not argue in the paper that SARS-CoV-2 could not have been found in places that were not sampled; I'm unaware of anyone making this argument based on our paper; it's not misleading if no one is misled.

      On top of this, the proportion of positive samples in each location is displayed along with the relative risk p-value. The n=1/1 sample that Bahry complains "shows little heat" is clearly indicated as a sampled location with high positivity in Fig 2A. And the underlying data is available in supplementary figures and tables. It's demonstratively not misleading.

      The other arguments the two commenters make consist of (1) arguing against a different paper published years ago, (2) a demand for a citation of an irrelevant paper, (3) arguing for an alternative analysis method without demonstrating that it would have more statistical power, (4) an argument that implicitly assumes all expectations for all coronaviruses are identical regardless of their hosts and modes of transmission, and lastly (5) a typo in a citation. We addressed the last one -- showing that, in fact, we didn't ignore these comments.

    1. On 2024-10-08 19:41:38, user Reena Sharma wrote:

      The key message of the study is that heavy metals like cadmium (Cd) and mercury (Hg) pose significant threats to plant health, but legumes, including Medicago truncatula, exhibit genetic variation in their ability to tolerate and accumulate these toxic metals. By conducting a transcriptomic analysis of plants with varying levels of tolerance to Cd and Hg, the study identified tissue-specific, genotype-specific, and metal-specific gene expression patterns.

      Notably, plants inoculated with mercury-tolerant rhizobia strains carrying a mercury reductase (Mer) operon experienced less reduction in nodule number, plant biomass, and iron distribution under Hg stress. This suggests that Hg-tolerant rhizobia can mitigate Hg toxicity in plants, enhancing resilience in contaminated environments. These findings highlight the potential to optimize legume-rhizobia interactions for improving plant tolerance to heavy metals and reducing heavy metal transport to edible parts of the plant, which is critical for food safety.

    1. On 2024-10-08 18:44:34, user Paul wrote:

      Please note that the preprint's result is subsumed by that in the published AMB paper. In particular, the AMB paper has modified Algorithm 1 so that it computes S in a right-to-left scan, rather than left-to-right. In this way, S does not need to be stored in memory, eliminating the quadratic space requirement. The space usage of the algorithm is now O(\ell \log \ell).

    1. On 2024-10-08 09:43:09, user Bruno Cenni wrote:

      Very nice and comprehensive dataset and overview across almost all BTKi. A note with regards to the data in Table 2 and Figure 4. For remibrutinib a BTK potency of 1.3 nM as “Kd or IC50” is listed. While the data is correctly referred to Angst et al 2020, this manuscript lists the IC50 for biochemical BTK enzyme inhibition. The same Angst et al 2020 publication also includes the Kd of remibrutinib for BTK (measured in the same assay all the others in the present manuscript) which was 0.63 nM. This is the value that should enter Table 2 and Figure 4.

    1. On 2024-10-08 00:33:10, user Alexis Rohou wrote:

      This comment is based on the version of the article published by Cell.

      I have concerns about the validity of the cryoEM result presented in this paper.

      The authors claim they obtained “a structure of the NINJ1 segment at 4.3 Å”. At this resolution, the following features should be clearly resolved in a cryoEM map of a protein containing alpha helices:<br /> • The pitch of alpha helices (~5.5 Å)<br /> • Large amino acid side residue side chains (e.g. Phe or Tyr, of which there are several in Ninj1)<br /> Neither of these features are demonstrated in the figures prepared by the authors.

      Figure panels S2D and E are consistent with a resolution of ~ 8-10 Å, where alpha helices are resolved as tubular features. No regular indentation or ridging corresponding to the helical pitch is apparent in the figures, or upon visual inspection of the deposited map ( https://www.ebi.ac.uk/emdb/EMD-42301) . No features corresponding to Phe or Tyr side chains are visible. At the resolution claimed by the authors, features corresponding to the side chains of residues Phe100, Phe117, Phe127 and Phe135, which are all located in alpha-helical segments (not in loops), would be expected to be resolved in the map, but they are not.

      In my view, based on this map inspection the authors should not have made this resolution claim.

      To support their resolution claim, the authors present Fourier shell correlation (FSC) curves from the cryoSPARC software in Figure S2C.

      While the FSC curve shown by the authors does cross the 0.143 threshold at ~4.3 Å, the FSC curve exhibits pathologies, which should have alerted the authors to the possibility that the 4.3 Å resolution estimate may be unreliable and that the map should be interpreted with caution. Most notably the curve starts dropping off at around 15 Å, indicating that the signal-to-noise ratio in the map is significantly deteriorated at resolution of ~15 Å and beyond.

      After completion of 3D refinement, cryoSPARC also outputs a second FSC figure, which includes an additional curve (“Corrected”), which accounts for effects of masking on resolution estimation by the FSC. It is unfortunate that the authors didn’t include this in their manuscript.

      The most likely explanation for this pathological FSC resolution estimate and its mismatch with the features resolved in the map is that the 3D refinement failed (due to high noise, or preferred orientation, or other pathologies in the dataset or in the refinement parameters), leading to overfitting to a local minimum in the scoring function.

      Indeed, the validation report for the deposited map and PDB (EMDB: 42301; PDB: 8UIP; https://files.rcsb.org/pub/pdb/validation_reports/ui/8uip/8uip_full_validation.pdf ) contains evidence of overfitting during refinement. The orthogonal projections of the raw map (Section 6.1.2 of the validation report) show overfitting artefacts, as do the orthogonal standard deviation projections in false color of the raw map (Section 6.4.2).

      Such overfitting causes artifically inflated FSC values, which may help explain why the FSC-based estimate of resolution was wrong in this case.

      Given a map of this quality, it’s unclear how the authors built an atomic model of Ninj1. In that respect, the publication’s methods section is not detailed enough. Table S2 indicates that the authors started from a computational structure prediction from AlphaFold2. Given the lack of any features in the map to help place any residue side chains, I assume the location of key residues mentioned in the paper originated from this computational prediction rather than from the cryoEM result itself.

      The lack of support for the modeled atomic coordinates from cryoEM is also made evident in the EMDB/PDB validation report by the unusually poor Q and atom inclusion scores for a map/model of the claimed resolution. A Q-score of 0.04 is unusually low and shows very little map-to-model fit.

      These metrics, together with a visual inspection of the map and model as deposited, suggests that a significant fraction of the sequence was built outside the map and that the cryoEM result does not support the authors’ atomic model beyond the general shape and relative orientation of parts of the alpha-helical segments. The positioning of individual amino acid residues is not directly specified by the cryoEM result and may have come mostly from the AlphaFold model the authors used as a starting point.

      Note: the above observations should be not be taken as having any bearing on other key cryoEM-based observations in the paper, such as the curvatures of Ninj1 assemblies, which are supported by two-dimensional class averages of cryoEM images and not affected by overfitting in the later 3D refinement process.

    1. On 2024-10-07 14:40:02, user BindCraft Enjoyer wrote:

      I like the ‘Design-Until’ architecture of the BindCraft pipeline, but one thing I couldn’t find in the paper is any quantification of BindCraft’s in silico design success rate. In the Introduction, you note that one drawback of RFDiffusion/MPNN-based pipelines is the need to screen thousands to tens of thousands of designs in silico before finding the 10-100 that pass the quality metrics and can be tested experimentally with good success rates. Does BindCraft also require screening of thousands to tens of thousands of designs, or is it more efficient in silico than an RFDiffusion pipeline? You mention in the paper that BindCraft outputs statistics from each design run, and that biasing away from alpha helical binders reduces the in silico design success rate; so it sounds like you have the statistics ready to hand, at least for the targets reported in the paper. I’d love to see these design success rates added to a table, either in the main paper or the SI.

      Another thing I’d like to see is some quantification of the compute time and cost required to run the 4-step pipeline until 100 designs pass the in silico filters. I understand this cost scales with target/binder size and target difficulty, but I would imagine you have the data required to calculate these metrics at least for the design campaigns reported in the paper. I saw on Twitter that you’re working on a direct BindCraft / RFDiffusion pipeline comparison; I hope you’ll include the computational hardware and total CPU/GPU time for each side of that design campaign.

      Great work!

    1. On 2024-10-06 00:56:48, user Rishav Mitra wrote:

      What is the endogenous TDP43 concentration in neuronal cells and how does that change with aging?

      Are there neuron-specific gene regulatory mechanisms like alternative splicing, microRNAs, transcription factors, etc. that control TDP43 expression and consequently the intra-condensate localization inside SGs and liquid-to-solid transition?

      Can microscopy as performed here distinguish between rapid coalescence of distinct TDP43 and G3BP1+ condensates and intracondensate localization?

      Is the phenotype titrable with concentrationof arsenite, i.e., more number of intracondensate TDP43 foci at higher [arsenite]?

      Does demixed TDP43 recruit/co-condense with other SG-resident proteins?

      How can we correlate age-associated oxidative damage with the length of arsenite exposure?

      Does solid-like demixed puncta of TDP-43 having higher local concentrations and slower internal dynamics compared to more liquid-like droplets potentially impact fixing and antibody staining? Is Correlative light and electron microscopy a better approach here than conventional immunofluorescence?

    1. On 2024-10-05 16:00:52, user Walter S Leal wrote:

      We are delighted to share this preprint on BioRxiv. It is the fruit of a couple of years of research collaboration with the citrus growers of the State of Sao Pauli. Fundecitrus has a unique model – they fund research (like most commodity groups do) and have PhD-level research scientists who engage in research activities. Please look at the number of replicates (as well as the raw data in SI) to further appreciate the work's depth. Walter Leal

    1. On 2024-10-04 16:24:29, user Gregory Way wrote:

      We read this paper as part of a journal club, and have decided to compile a collective review and publicly share it with the authors. This was inspired by the Arcadia Science Preprint Review Pizza Party Initiative, and this represents our fourth preprint review.

      Ji et al. present a transformer-based foundation model, called Prophet, which stands for Predictor of Phenotypes. The authors train Prophet on a variety of data modalities including gene expression, cell viability, chemical structures, and cell morphology using publicly-available sources. The authors should be commended for using such vast and disparate resources for such an innovative approach. The task of Prophet is to predict assay endpoints, such as cell viability, and to learn a useful embedding space which can be mined to identify novel, and potentially impactful, relationships. Most often, the phenotype prediction is in the context of some form of perturbation. The authors present a variety of benchmarks comparing Prophet to other methods, and they present both in vitro and in vivo applications to demonstrate potential use-cases. The applications range from looking up untested compounds that are similar to clinically-relevant compounds and predicting zebrafish cell type proportions after gene knockout. Overall, Prophet is methodologically interesting and the applications demonstrate that the method may help generate hypotheses at a low cost. However, we have several major and minor concerns mostly to do with clarity, performance, and software.

      Major concerns:

      1. Unclear and inconsistent terminology and definitions.<br /> a. It is unclear exactly what the authors mean by “phenotype”. It seems that sometimes the term is being used interchangeably with prediction/output but other times it is being used to describe observable physical properties. Additionally, the authors refer to gene expression as a phenotype, and, while technically true, it could be confusing given the authors are also using cell morphology as a phenotype. Furthermore, it is unclear if the authors are describing the collection of genes in the gene expression vector, for example, as the phenotype, or, if it is just a single gene. This confusion is also related to our confusion about model outputs (whether the output is a single element or a vector representation; see below). Given the word is in the article’s title, it seems particularly important.<br /> b. The terminology of “experiments” is also unclear. The authors claim to use 4.7 million experiments, but does this refer to plates, conditions, samples, something else? How did the authors calculate this count?<br /> c. At times, it is unclear what format the input data are and at what level of processing. What are the different possibilities of input data and how might a user decide which to use? Can a user input multiple kinds of data? Did the authors apply any sort of post-processing or quality control?<br /> d. The output of Prophet is ambiguous. Does Prophet predict a single value per input (or different inputs), or, does it predict a cell state vector? The authors describe outputting a one-hot encoding. Does this refer to the output phenotype? What is this structure? The authors write: “we train Prophet to predict cell viability, compound IC50, Cell Painting morphology features, mRNA transcript abundance, and cell type proportion.” Does this mean Prophet will output all of these predictions if the data you have is only morphology features? Does a user have control over these decisions? Furthermore, the authors write "Prophet’s transfer learning capability is not limited to phenotypes seen during the pre-training stage. We did not pre-train Prophet on any morphological measurement, but Prophet fine-tuned on JUMP outperformed both the Prophet-individual model trained only on JUMP and the baseline (Fig. 2b).” What is being predicted from the JUMP data? Morphology feature profiles? Images? Drug class or MOA? This needs to be expanded upon to make these claims. Please clarify the output structure and how a user will interact with the output. Figure 1C does not make this clear.<br /> e. Critical methodological details are discussed without sufficient detail. For example, the data split and validation strategies were ambiguous. How were training, test, and validation splits handled? What partitioning methods, if any, were used? The three-fold cross-validation procedure also lacked clarity. Were all datasets used in cross-validation? How did individual data-set models training differ and influence the full model fine-tuning? What is the specific pseudobulking procedure for RNA?
      2. Unclear justification for Prophet architecture decisions<br /> a. The authors present table S3, which provides hyperparameters. How are these justified? For example, the choice of GeLU over alternatives like ReLU or SiLU. Why use an embedding dimension of 512? Were alternative configurations explored? How would modifying individual architecture decisions impact performance?<br /> b. 20.1 M parameters is a fairly small transformer, will this model need to grow as more perturbations or data types are added? In other words, how long will this model be foundational until the next best model is released? <br /> c. Encoders sequentially relate inputs with a positional embedding. Does this architecture use a positional embedding in the encoder? What does the position represent?
      3. Concerns about model comparisons, baselines, and performance<br /> a. The authors compare Prophet to much simpler machine learning models (random forest, MLP, linear regression), individual Prophet models (trained using only one modality), and a mean baseline representing the average value of that intervention. The authors write: “This approach follows the same strategy as current foundation models (19, 36, 37), which are pre-trained on large amounts of data and then fine-tuned for specific datasets using the pre-trained model as a backbone.” Why not use these current foundation models as benchmarks for Prophet? The authors should also consider comparing different transformer architectures and non-transformer models (e.g., state-space models) as well.<br /> b. Prophet’s performance is low. The highest R2 value is 0.27 with a low of -0.03 and many predictions that perform the same as the mean baseline. Given the low, variable performance, it is difficult to trust Prophet’s output, or, at best, understand which outputs may have incorrect predictions. The authors claim an R2 improvement as low as 0.04 represents a 13x increase in number of hits, but it is unclear how the authors calculate this value. The authors also claim Prophet reduces “the number of experiments needed for viability screens by at least 60x” What statistics are calculating this estimate? <br /> c. The ML model comparisons compared to baseline are incredibly low. Results in Figure 2B suggest that the mean is a better predictor than nearly 100% of ML-based predictions (mean baseline is better in 41/45 comparisons). Our guess is that something might be going wrong in the model training or evaluation procedures.<br /> d. It is confusing if there is only one “final” Prophet model or if there are multiple “final” Prophet models since each must be fine-tuned to single datasets. Why not fine-tune on all datasets? Figure 2A suggests that for each time you perform inference in a new data category (e.g., cell morphology vs. gene expression), fine-tuning is required? If so, this strategy will lead to variable predictions and perhaps unexpected results and this is not a foundational model, but, at best, a foundational architecture.<br /> e. The authors performed a critical analysis, in which they restricted the amount of training data by 50%, 30%, 20%, 10%, and 5%. The authors state: “We found a clear trend: the more treatments and cell states seen by Prophet, the higher the confidence in the predictions (Fig. 2g).” This is an expected result, however, it is unclear how Figure 2G supports this statement. What is an “axis holdout”? Standard deviation of R2 for which predictions? What do the different points represent? How does this show confidence?
      4. Concerns about limited discussion of technical artifacts<br /> a. The authors mention technical variation in the limitations section, but this should be elaborated upon. How might Prophet’s performance be impacted by technical artifacts?<br /> b. Earlier in the manuscript, the authors write: “We universally decompose each experiment into a unique combination of three fundamental elements—the cellular state, the treatments being performed, and the intended phenotypic readout.” Signals from technical artifacts are likely a fourth fundamental element, or, at the least, there should be an experiment to test the impact of technical artifacts.
      5. Concerns about publicly-available source code on Github<br /> a. We provide a full GitHub review following our manuscript review below.

      Minor concerns<br /> - Pg. 3 ln 5, the authors write: “To train it, we collected 9 perturbational datasets to create the largest compendium of publicly available screening datasets to date:...” The language suggests that these datasets were all collected by this group, but the datasets are all publicly available. Also, the Figure S1 reference probably should point to Figure S2. JUMP is listed three times in Figure S2 (one for compound, one for genetic treatment, and then one for both) Why is it listed for both when it is already split into the two perturbation types? Also, the right subplot in S2B (the complexity trade-off) is a bit misleading. JUMP has far higher complexity than PRISM, for instance, but this graph would suggest otherwise. Perhaps the missing piece not described here is readout complexity?<br /> - Figure 1B is a bit difficult to understand. For example, the bullet points for the “intervention” box don’t seem like interventions? It seems like these should be listed elsewhere or the heading should be changed. How is a “gene sequence embedding” an intervention data transformation?<br /> - The authors should tamp down claims. For example, the authors write “We did not pre-train Prophet on any morphological measurement, but Prophet fine-tuned on JUMP outperformed both the Prophet-individual model trained only on JUMP and the baseline”. The performance is only marginally improved (0.01). The authors retrained the zebrafish Prophet model, but it would be helpful to see performance for the original prophet model applied.<br /> - What is a pre-built plate? “To do so, we used a setup with existing pre-built 384-well plates, each with 352 unique drug perturbations applied to 9 cancer cell lines (Table S9). In total, there were 16 pre-built plates.”

      GitHub comments and concerns<br /> 1. Documentation and Usability<br /> a. The README provided by the author is well-structured, offering clear instructions on installation, usage, and licensing, which provides a strong starting point for using Prophet. This level of clarity is especially valuable for researchers and users who are new to the tool. <br /> b. While the README provides a good overview, the documentation around model training is sparse. It would be beneficial to include an explanation in the README on how the model was trained and provide a small explanation of the embeddings captured. More detailed usage/inference examples would enhance comprehension. It’s also unclear how users can apply the method to their own datasets. This will offer quick and easy access for users to understand the functionality and purpose of Prophet.<br /> c. Commented-out code and TODOs are scattered throughout the scripts. Best practices suggest removing unused code to reduce clutter and confusion. Example: prophet/ http://model.py #L12.<br /> 2. Software Environment and Dependencies<br /> a. Users may face reproducibility challenges when trying to set up the software, particularly due to missing environment isolation. Creating an isolated Conda environment or improving instructions around environment setup would help ensure users avoid dependency conflicts.<br /> b. The project uses an older version of Pandas (1.5.x), despite newer versions being available with important fixes. Updating the Pandas version would improve compatibility and performance.<br /> c. A http://setup.py and requirements.txt are both provided but are not used together, creating confusion over proper environment management.<br /> d. No Python version range is specified in the http://setup.py , which led to issues with earlier Python versions (3.8, 3.9). Python 3.10 worked, but this should be clarified for future users.<br /> e. The provided notebook requires a specific version of NumPy that differs from what is stated in http://setup.py . Errors occur with newer versions. NumPy 1.24.4 was found to work, but this should be addressed in the dependencies.<br /> f. We were able to install it into a Linux machine. However, an error occurs when attempting to install the software on macOS. The error reports: “ERROR: No matching distribution found for scipy==1.14.0”<br /> g. A Jupyter notebook is included, but Jupyter is not listed as a dependency in the http://setup.py or requirements.txt, which prevents seamless execution within the provided environment.<br /> 3. Other comments<br /> a. The repository does not include any software tests or automated testing via GitHub Actions or similar tools. Incorporating automated testing would help validate the code’s functionality and improve its robustness.<br /> b. The code does not pass several linting checks (e.g., through dslinter), highlighting the need for improved code quality and adherence to data science best practices.<br /> c. The repository lacks key community health files like http://CONTRIBUTING.md and http://CODE_OF_CONDUCT.md , which are important for guiding open-source contributions and user interactions.<br /> d. Data provided in .xlsx format (e.g., via Figshare) can cause formatting errors and is less open-access friendly than formats like .csv. Switching to a more stable format would improve accessibility and avoid errors.<br /> e. In the tutorial notebook, the code blocks are not executed in a sequential order, which can lead to potential bugs. This lack of sequential execution means that changes made in earlier cells may not be reflected in subsequent cells, resulting in inconsistencies or errors in functionality. <br /> 4. Recommendations for Improvement<br /> a. Address installation issues: Fixing the bugs related to installation and setup should be a priority, as they could deter users from exploring Prophet further. Ensuring an isolated environment setup (e.g., using Conda) would help resolve these issues.<br /> b. Enforce version control for dependencies: Better organization of http://setup.py and requirements.txt, as well as enforcing version control for dependencies, would enhance reliability.<br /> c. Expand the README: Adding a table of contents, additional usage examples, and sections on contributing guidelines and testing procedures would make the README more comprehensive.<br /> d. Adopt best software practices: Implementing clear setup instructions, enforcing dependency version control, and organizing the code more effectively would increase usability and accessibility for a wider scientific audience.

      This is a signed review:<br /> Gregory P. Way, PhD<br /> Erik Serrano<br /> Jenna Tomkinson<br /> Dave Bunten. MEd<br /> Michael J. Lippincott<br /> Cameron Mattson, MSc<br /> University of Colorado Anschutz Medical Campus, Department of Biomedical Informatics

    1. On 2024-10-03 21:49:13, user Francesco Del Carratore wrote:

      At the end of the methods section it is written 'To facilitate reproduction of these findings, all shareable data and code are available in a single structured file, with instructions and links for the non-shareable data, in S1 Data.'. This is great, but where can I find the S1 data as well as the code used for the analysis and figures (S1 code and S2 code)?

    1. On 2024-10-02 20:15:19, user Prof. T. K. Wood wrote:

      1. The authors fail to cite the relevant literature:

      a. Line 59: the first report of ROS inducing persistence for any species is doi:10.1111/j.1751-7915.2011.00327.x (published in 2012, 12,000-fold increase) and should be cited. An erratum was issued to ref 20 by the Conlon group for failure previously to cite this reference.

      b. The first high-throughput screen (10,000 compounds) for waking persister cells was doi:10.1111/1462-2920.14828 (published 2020), which reports persister cells wake by modifying ribosomes with pseudouridine.

      c. The first report showing lower ATP increases persistence is doi:10.1128/AAC.02135-12 (published in 2013) and should be cited.

      1. It is odd that the S. a. strains have 10% persistence (Fig. 1AB) when most strains have persistence far less than 1%.
    1. On 2024-10-02 18:44:54, user Yasas Wijesekara wrote:

      Congratulations on your discovery! This indicates that there is still so much biology waiting to be discovered. I'm interested in examining the genes on these Inocles myself. However, using the provided accession number, I couldn't locate the sequences. Would they be available soon?

    1. On 2024-10-02 16:40:25, user Iva Tolic wrote:

      This manuscript is related to our manuscript "Kinetochore-centrosome feedback linking CENP-E and Aurora kinases controls chromosome congression," doi 10.1101/2024.09.29.614573. This manuscript here expanded significantly during the revision from version 1 to version 2. As a result, we decided to divide it into two separate manuscripts: version 3 (doi: 10.1101/2023.10.19.563150v3) and a new manuscript (doi: 10.1101/2024.09.29.614573).

    1. On 2024-10-02 16:36:02, user Iva Tolic wrote:

      This manuscript is related to our manuscript "CENP-E initiates chromosome congression by opposing Aurora kinases to promote end-on attachments," doi 10.1101/2023.10.19.563150v3, which expanded significantly during the revision from version 1 to version 2. As a result, we decided to divide it into two separate manuscripts: version 3 (doi: 10.1101/2023.10.19.563150v3) and this one (doi: 10.1101/2024.09.29.614573).

    1. On 2024-10-02 09:06:19, user Frédérique Reverchon wrote:

      We discussed this paper in our Journal Club and the comments arising from our Discussion are the following:<br /> 1) Why focus on “rhizosphere soil” rather than “ bulk soil” to study the dynamics of fungal functional groups in fields of different age? The rationale behind that choice is not clearly explained, most of the references used in the Introduction and Discussion sections to explain the effect of time since planting on functional groups come from soil studies rather than rhizosphere studies. The term “field age” is also misleading as it suggests that bulk soil was sampled, and other terms such as “plant age” could maybe be considered.<br /> 2) There were no replicates per field age, as only one field per age was considered, although we understand that different plots were sampled per field. Other parameters than “field age” could cause the observed effects on fungal functional guilds, such as differences in soil nutrients between fields, for example. As seen in the PCoA, fields C and D (both classified as “old fields”) are different in terms of fungal community structure. Including replicate fields would have allowed to avoid possible confounding effects.<br /> 3) Important details are missing from the field description. Considering the Discussion on the possible effects of fertilization and fungicide applications, information should be added on the agrochemical management of these plots. Was it similar across fields? Across sampling years? Information regarding the soil nutrient status, and possible climatic variations between the two years of sampling, should also be included.<br /> 4) Methodological bias was also discussed. For example, most fungal taxa identified with the UNITE database would not be able to be assigned to a functional guild. How did you account for these taxa that could not be included into a functional category? Another possible bias that was discussed is the fact that relative abundance data may not represent what is actually occurring in the rhizosphere: absolute abundance data would also be needed and should be used for the implementation of microbial networks.<br /> 5) The discussion could also benefit from adding some consideration on the “functional group” categorization. How categorical is a “functional group”? Mortierella for example is listed as a saprotroph but it can also act as a plant pathogen. Mycorrhizal fungi included both EMF and AMF, yet strawberry only form symbiosis with AMF. As the focus of this study is on strawberry rhizosphere, why not discuss possible strawberry pathogens and symbionts?<br /> 6) An effect of plant genotype was not detected on relative abundance nor richness of fungal functional guilds. However, significant effects may be found in terms of taxa composition. Why not include such analysis in the study?<br /> Other minor comments are enlisted below.

      Abstract<br /> Results seem to report differences in soil fungi rather than rhizosphere fungi. Why the focus on rhizosphere and not bulk soil? Plant genotype is not mentioned in the Results either, yet it was mentioned in the Background section as an important determinant of microbiomes.

      Introduction<br /> L50: “The essential functions of soil microbiomes, such as pathogenesis, mutualism, and decomposition…”. Is pathogenesis a function? Functions of the soil microbiomes could rather include maintenance of soil structure, fertility, plant productivity… <br /> L66: Research gap: “Despite their significance, our knowledge of the dynamics of different fungal functional groups that co-exist in rhizosphere soil microbiomes is limited (Zanne et al. 2020; Martinović et al. 2021)”. The importance of filling this research gap is not clearly explained. <br /> L69: “Soil fungal functional groups can be dynamic in the plant rhizosphere. For example, fungal pathogens can accumulate over time, especially in monocultures or ecosystems with low plant species diversity (Cook 2006; Li et al. 2014; Peralta et al. 2018; Wang et al. 2023a)”. These references are for soil, not rhizosphere.<br /> Second paragraph of the Introduction: mostly on soil build-up of pathogens, and possible decrease of AMF over time. But this section seems to address the effect of long-term cropping, not plant age (rhizosphere). The rationale behind the selection of rhizosphere as the soil compartment to be studied is needed. You may want to check paper by Sun et al. (2022), published in Applied Microbiology and Biotechnology, on the dynamics of fungal functional guilds in the rhizosphere of wheat.

      Methods<br /> Why establish and sample plots (20 x 20 cm) and not plants if rhizosphere soil was considered? Is there only 1 plant per plot (i.e. 3 plants sampled per genotype and field age)? It is unclear.<br /> Were the same plants sampled in 2021 and 2022? If not, a plant effect could be observed between both years.

      Results<br /> L205: The increase in pathogens in the older fields is observed for one sampling year only (the pattern does not hold in 2022), yet it is described (and discussed) as a general finding.<br /> What happened in 2022 with AMF that decreased dramatically in all fields (both in relative abundance and in richness)? It is briefly discussed as possibly caused by the fertilization regime, which is not described.<br /> Differences in community composition could be assessed statistically to test for significant differences in taxon relative abundance between fields or sampling years. <br /> Regarding your Figures, stating the age of the field in the “x” axis rather than letters “A”, “B”, etc., could be more informative.

      Discussion<br /> The discussion on the accumulation of pathogens over time is also based on references on soil microbiomes, not rhizosphere microbiomes. Were the same plants sampled in 2021 and 2022?<br /> L301: fungicide applications are mentioned to explain the “levelling off” of pathogens in older fields (C and D). Any data on application rates and application dates? This is important if fungal communities are studied. How specific were these agrochemicals? How were they applied?<br /> L311-312: Fusarium is a large genus with many species, not all are pathogens. Making inferences such as “the causal agent of Fusarium wilt in strawberry” is a bit far-stretched. Which Fusarium species are actual strawberry pathogens?<br /> The whole discussion regarding the potential accumulation of pathogens in the rhizosphere calls for symptom incidence data. It is something that was measured?<br /> L387: “the similarity in fungal functional groups across genotypes observed in this study may indicate functional redundancies and stable essential functions of microbiomes.” The microbiome was not studied and differences may be reflected in bacterial communities rather than fungal ones.

    1. On 2024-10-01 18:51:48, user Maria Valle wrote:

      This review was done as part of the SfN Reviewer Mentor Program (Mentor: Joanne Conover, PhD; Mentee: Maria Luisa Valle, PhD)

      Manuscript title: Human iPSC-derived pericyte-like cells carrying APP Swedish mutation overproduce beta-amyloid and induce cerebral amyloid angiopathy-like changes<br /> Journal: bioRxiv

      Overview<br /> Wu et al. characterized human induced pluripotent stem cell (iPSC)-derived pericyte-like cells (iPLCs) to investigate the role of pericytes in Alzheimer’s disease (AD) and cerebral amyloid angiopathy (CAA). First, the authors showed that iPSCs could efficiently differentiate into pericyte-like cells, express pericyte specific markers, and promote angiogenesis, barrier integrity and contractility. They then investigated the differences between iPLCs derived from healthy individuals and those derived from AD patients carrying APPswe mutations. Compared to controls, APPswe iPLCs exhibited a distinct expression of pericyte markers, were able to secrete amyloid beta 1-40 and 1-42 within the media, had an altered transcriptome for key genes involved in cytoskeleton reorganization and metabolic regulation, were more sensitive to mediators of inflammation, and showed compromised angiogenesis, barrier integrity and hypercontractility.<br /> Overall, the manuscript uses a novel approach (iPLCs) to investigate an interesting and sometimes overlooked topic - the specific contribution of pericytes to AD pathology and vascular disfunction. Previous work conducted in in vitro BBB models showed that pericytes play a key role in amyloid clearance contributing to the removal of aggregated Aβ from brain capillaries (Ma et al, Mol Neurodegeneration 13, 57 (2018) and Blanchard et al, Nat Med. 2020 Jun;26(6):952-963). Here, the authors focus on the contribution of pericytes in amyloid secretion, emphasizing the novelty of their research. However, the high variability within datasets and the small number of replicates raises some concern.

      Major comments <br /> • The statement “…overproduce beta-amyloid” in the manuscript title suggests that pericytes have a significant role in Aβ production. Although the authors showed that APPswe iPLCs could secrete 10 times more Aβ1-42 than the control cells, the Aβ1-42 levels are 100 times lower than neurons. Thus, the authors concluded that “contribution of pericytes to total brain amyloid load in AD is limited”. The title should be changed to indicate the main findings of the work and should be supported by the data presented. <br /> • APPswe iPLCs were derive from 3 donors versus iPLCs from 7 healthy controls. Importantly, among the donors, only one had AD, while the others had pre-symptomatic AD or no symptomatology (in this case the mutation was introduced using CRISPR-Cas9 as reported in the methods). The variability in AD cases plus the differences in symptomatology may skew the results and may contribute to the high variability shown in several graphs (Figure 2 A, B, C, F, J).

      Figures<br /> • The authors should be consistent in the number of replicates used: different groups in the graphs show only 1 or 2 replicates, even for control cell lines, which makes the reader question the reproducibility and accuracy of their findings (see Figures 1B, 1I, 1K, 2H, 2J, 4E). <br /> • The authors should clarify the findings reported in Figures 2E and 2H: the figures are similar, but it is not clear if iPLCs in 2H derive from APPswe iPLCs (as reported in the figure legend) or control.<br /> • The authors should correct Figures 1A and 2I as sample labels are missing. The authors should also modify the arrows used in Figure 2I and 2D as it is not clear to what they are pointing. Scale bars should be added on both images since they show different magnification. <br /> • Figures should be arranged in a consistent manner e.g., same format and order should be used consistently.

      Discussion <br /> • An interesting finding is that the HIF1a pathway is downregulated in APPswe iPLCs (Figure 3B). The authors should mention this finding in the discussion. This finding could also support the fact that APPswe cells have decreased VEGF levels and impaired angiogenesis and no change in BACE1 levels (as VEGF and BACE1 are HIF1a target genes). <br /> • For future experiments, the authors should discuss whether APPswe iPLCs exhibit differences in oxidative stress, ROS production and mitochondrial activity compared to controls.<br /> • For future experiments, the authors should use cell lines and human-derived cells as models as they may reveal differences from iPLCs.

      Minor comments<br /> • The authors measured the changes in expression of several pericytes associated genes in Figure 1. However, it is not clear why the authors were not consistent with these specific genes for their further analysis. For example, in Figure 1B they measured PDGFRB, DES, LAMA2, DLC1, and PDE7B while in 1C they measured LAMA2, PDE7B, DES, omitting PDGFRB but adding genes ACTA2 and CD248. Then, all genes were analyzed in Figure 2A-B. Thus, the authors should provide change in expression data for all genes (PDGFRB, LAMA2, DLC1, CD248, PDE7B, ACTA2, DES) in 1B and 1C or provide reasoning for leaving some out. <br /> • Please correct the repeated sentences on page 5: “…which are known to express detectable levels of LRP121 (Figure 2 J). Furthermore, when iPLCs were subjected to pHrodo-conjugated zymosan-coated beads, no uptake of these pathogen-mimicking particles was observed (data not shown). Thus, it appears that the phagocytic activity of these iPLCs is low.”<br /> • Additional edits for word choice and sentence construction are also needed, e.g., pg 10, 2nd to last paragraph, 2nd to last sentence is awkward.

      Decision for the editor: Major revisions<br /> The manuscript presents a novel idea that could advance the AD/CAA field but, at this stage, I have several major concerns regarding reproducibility and possible accuracy of the described findings. I would consider the manuscript for publication only if all major concerns are addressed by the authors.

    1. On 2024-09-30 21:25:22, user Swagatam Barman wrote:

      While the use of a massively parallel combination screen to identify an adjuvant for enhancing the efficacy of rifampicin is innovative, its practical implications are limited. Given that rifampicin is primarily restricted for treating tuberculosis, the likelihood of translating this adjuvant-rifampicin combination into clinical practice is quite low.

      To improve the relevance of this research, it might be more beneficial to explore alternative drug combinations or focus on broader-spectrum agents that could be used in various clinical settings. This approach would enhance the potential impact of the findings and increase the likelihood of clinical applicability.

    1. On 2024-09-30 15:40:53, user Christopher Dunn wrote:

      I left a comment earlier, but it doesn't appear, so I am trying again.

      This is an interesting paper, but I am not sure how wide an interest it will achieve. That aside, I need to re-read this and provide more detailed comments.

      That said, I would note that the official state flower of Connecticut is not the marvel of Peru (Mirabilis jalapa). That species is the State's "Children's flower."

      The official state flower of Connecticut is Kalmia latifolia (mountain-laurel). This is the species that should be used in the analysis.

    1. On 2024-09-27 17:21:46, user Frank wrote:

      After the CcFV-1 paper, which provided rather weak evidence for the virus being capsidated (no infectivity in addition to the rather weak blotting experimental designs), I find it difficult to agree with the wording of "strongly support the notion that polymycoviruses are encapsidated" being used by the authors again this time. Albeit the evidence is stronger here, there remain other plausible interpretations.

    1. On 2024-09-27 11:24:18, user Ruth Berger wrote:

      Important, high quality research. There is a potential alternative explanation for the pattern observed which should be checked: Viola arvensis is a pretty rare wildflower where I live, and in any peri-urban environment must be much, much rarer than garden cultivar violas that usually don't offer any pollinator food at all. Could it be that pollinators avoid them because of their experience with garden variety violas that taught them viola-like flowers are useless? From observation, I have seen various wild and cultivated viola species not getting any pollinator visits at all despite pollinators busily visiting other flowers in the immediate vicinity.

    1. On 2024-09-26 20:35:22, user Trịnh Gia Huy wrote:

      Hi Lacle, thank you for your comment. We fixed the figure already. The completed version and the code will be released later.

    2. On 2024-08-08 13:15:41, user F. Laclé wrote:

      Also, if you can publish your model code in a repository would be great for reproducibility (the model itself is not necessary I reckon). As you know, there are much more system configuration elements to consider, which makes reproducibility efforts complicated. Publishing your model code would allow others to attempt and improve the reproducibility challenges.

    3. On 2024-08-08 13:11:45, user F. Laclé wrote:

      Hi. In Figure 6b, shouldn't the number of filters of the skip connection be equal to 4 times the number of filters in the last Xception module? You have it specified as 16, shouldn't it be 32?

    1. On 2024-09-26 15:09:36, user pLM Enjoyer wrote:

      An important application of pLMs is enhancing the efficiency of protein engineering by adding a classifier top model onto a foundation pLM (with or without fine-tuning), training on a small number (0-96) of experimental sequence/fitness datapoints, and then using this model to score and predict high-fitness sequence variants. This task also provides a good benchmark of pLM quality, since pLMs with ‘better’ embedded representations of sequences produce better variant scores/suggestions. See, for example, Zhou et al 2024 (Enhancing efficiency of protein language models with minimal wet-lab data through few-shot learning), and Jiang et al 2024 (Rapid protein evolution by few-shot learning with a protein language model). I think it would be really useful for you to benchmark AMPLIFY’s performance against ESM/SaProt/etc on few-shot and zero-shot variant fitness prediction using public deep-mutational scanning datasets as described in the papers above. If AMPLIFY really outperformed ESM2-15B on this task, that would be huge!

    1. On 2024-09-26 14:13:00, user Dave Grainger wrote:

      Hi, thanks for the great tool. I was looking through your figures for a journal club. I think there may be a slight error in panel c of the figure. The dotted outline does not highlight the same region in the b2c image as it does for 8um and MERFISH. Worth checking and updating for the final manuscript. BW, Dave

    1. On 2024-09-26 13:55:29, user Paola Casanello wrote:

      Very interesting data!<br /> We showed some years ago that adiponectin had vasodilator effects in the chorionic arteries (wire myography), and that the offspring of women obesity had a limited vasodilatory response to adiponectin.<br /> It would be interesting to discuss these results with your data.<br /> Muñoz-Muñoz et al., J Cell Physiol 2018. doi: 10.1002/jcp.26499.

    1. On 2024-09-26 02:21:06, user Jeff Holst wrote:

      Manuscript has been accepted for publication at EMBO Journal (24th September 2024). We will post a link to the open access revised version once it is available online.

    1. On 2024-09-26 01:06:46, user Kate wrote:

      This is the first author commenting here. Not sure where to discuss after the paper was published, but wanted to add some insight regarding the endogenously purified Pks13 protein used for cryoEM, crosslinking mass spec, and LC-MS (for identification of endogenous lipid). As it's published at NSMB, the paper doesn't address the fact that the endogenously purified Pks13 used for the above mentioned experiments, using the same purification protocol each time, showed variable SEC traces. After a while, we could not reproduce the peak1 and peak2 peaks in SEC (shown in Extended data fig1), but the two peaks either overlapped, or peak1 essentially disappeared. Sometimes, other interacting proteins were also pulled down along with Pks13, which confounded the SEC traces. CryoEM was done with peak1 with TAM inhibitor added, while XL-MS and LC-MS were done with mixed peak1&2 species. Protein used for XL-MS also had "contaminant" species which could actually be functionally-relevant binding partner. Since these protein preps behaved differently in SEC, we are not sure about the functional ramifications of these proteins purified at different times. If others are trying to reproduce our purification results, they may come across these variabilities (which is reasonable, given that this is an endogenous purification!).

    1. On 2024-09-25 14:40:35, user David wrote:

      What a great paper and story! It raises so many questions about the physiological relevance of such a mechanism specifically displayed by the gonadotrope cells! We were particularly interested to see that you identified Neurod1 and Neurod4 as being upregulated during postnatal differentiation of gonadotrope cells.

      As we have recently shown in vitro (PMID: 37658038) that NEUROD1 and NEUROD4 regulate the mobility of immature gonadotrope cells by regulating the expression of NTRK3, we were wondering if you could identify this gene as well? It would be very interesting to know whether the mechanism that we have identified in vitro and in vivo as regulating gonadotrope positioning in the developing pituitary might be relevant to the process that you have just described.<br /> David and Charles

    1. On 2024-09-24 16:44:20, user The Fehr Lab wrote:

      These authors have done a fabulous job at creating new macrodomain inhibitors, which is extremely appreciated. However, having a major conclusion and an implication that Mac1 inhibitors are not antiviral based solely on negative data is misleading. We have published that a macrodomain inhibitor can inhibit virus replication (PMID: 38592023) and will have another story that will soon be available in BioRxiv that describes even more Mac1 targeting compounds that inhibit virus replication. There are some notable problems with the compound described here that could explain its lack of antiviral activity that should be taken into account. Again, I think the novel chemistry identified in this paper is exceptional, but more cautiousness should be taken before making broad claims that this is not a good drug target. Based on genetic data Mac1 and another recent paper on bioRxiv ( https://doi.org/10.1101/2024.08.08.606661 ), it appears that Mac1 is a suitable target for antiviral development, and we are continuing to work to see that dream come to fruition.

    1. On 2024-09-22 21:30:29, user Christian Helker wrote:

      Beautiful work!!! :)<br /> I would like to bring to your attention our publication (“Apelin signaling drives vascular endothelial cells toward a pro-angiogenic state”; https://elifesciences.org/articles/55589) , which explores the function of Apelin on the vasculature. I believe it could provide additional context or complementary insights to your work.

    1. On 2024-09-22 20:05:49, user Fraser Lab wrote:

      Based on: https://www.biorxiv.org/content/10.1101/2024.07.24.604935v1.full

      The manuscript from Lehner and colleagues presents a wealth of mutagenesis information on amyloid aggregation. The central premise of the paper is to use a yeast selection based around the oligmerizaton/aggregation of Sup35 fused to a peptide (in this case abeta) as a proxy for amyloid forming potential. This is cool information on its own and the experimental analysis and computational framework for linking to energies is top notch. The point of using double mutants to enhance the dynamic range is very well explained and will solidify their approach to impactfully link DMS experiments to thermodynamic concepts.

      The major framing of the paper revolves around an analytical protein folding/engineering concept of phi-values that highlights energetic differences in the importance of interactions for forming the transition state vs. the ground state. For the textual interpretation of the results, one must buy into the energetic effects of the Sup35 system as a readout of the transition state (and secondarily for FoldX calculations on various PDBs of abeta polymorphs to be a readout of the ground state). The major issue is that an alternative (and perhaps simpler) explanation is that mutations in APR2 are more disruptive to the Sup35 oligomerization process in the screen and that this reflects amyloid/oligomerization propensities and not strictly TS of initial nucleation. The data from previous studies that is used to draw correlations to justify their interpretation around the TS is buried in extended data figures and is a bit all over the place, especially the deconvolution of primary vs. secondary nucleation. The existence of multiple polymorphs in human cells (and populations), which may or may not have related transition states - and the exact conformational requirements of the Sup35 activation mechanism - further complicate this interpretation.

      In summary - this contains amazing data, but I do not see the language of the interpretations lining up with the strength of the _specificity_ of the claims about the transition state. A fuller discussion of the limitations of the prior low throughput assays that are referenced in extended data 1 and 2 and detailed kinetic characterization of some of the more surprising mutants in a biochemically defined system would greatly improve the match between the data and the claims. These issues should not stop others from building on this beautiful work - but in doing so, other investigators should note that there remains ambiguity as to whether the effects are truly on the TS or on the ground state.

      Avi Samelson and James Fraser

    1. On 2024-09-21 18:49:07, user Flo Débarre wrote:

      Following up on my previous comments about the pangolin datasets featured in Figure 2:

      It had been suggested that NCBI could have tampered with the dates and times shown on their systems. To confirm that a request had been made on June 16, 2021 to make the data public again, I FOIA'd NIH with more search terms than had been done in previous requests. I finally received an email with the user's request, and I can therefore confirm that the pangolin CoV data being public again is unrelated to Jesse Bloom's preprint having been sent to bioRxiv or to NIH on June 18.

    1. On 2024-09-19 14:10:01, user Farhan Feroze wrote:

      Excellent work!<br /> I am curious about the reasons why pulse code EO-151 was preferred over EH-115?<br /> Also, were whole plasmids used as a HDR templates for electroporation? (Since we usually deliver the HDRT as linear dsDNA or ssODN with exposed homology ends)

    1. On 2024-09-18 21:12:56, user Jason Hoskins wrote:

      I am a fan of this approach of pairing multidimensionality reduction of the expression data with biologically meaningful gene sets to generate representative scores of processes or pathways that may be used as quantitative traits in QTL analyses. Approaches like this applied to GWAS variants enable the implication of processes or pathways potentially mediating the genetic effects on complex phenotypes even in the absence of co-localizing cis-eQTL signals, which is unfortunately typical of GWAS signals.

      This pre-print does not include a Discussion section, so it does not yet sufficiently place this work in the context of the relevant published literature. When this is done, I would recommend consideration of our work on expression regulator activity QTLs (aQTLs) that was published in PLOS Computational Biology a few years ago ( Hoskins et al., 2021 ). Our two approaches share a lot conceptually in that we both assume that shared variability among genes often reflects a common underlying regulatory mechanism and that multidimensional reduction of the expression information among such genes can provide a useful metric for summarizing the activity of the gene set. The approach in this pre-print is generalizable for use with any type of gene set deemed potentially relevant to the GWAS trait of interest, while our approach in Hoskins et al. (2021) is more tailored to gene sets representing the target genes of expression regulators based on tissue or cell type-specific co-expression networks where the activity scores are inferred using the VIPER algorithm developed by Andrea Califano’s group (Alvarez et al., 2016). There are also some analyses, results, and speculations in our paper that might suggest potential additional considerations and investigations for this gsQTL study.

    1. On 2024-09-18 14:19:15, user Musaeum Scythia wrote:

      I was curious on which basis the authors wrote the following section:

      "One particularly intriguing finding is the identification of the Y-chromosome haplogroup Q1b1a3a1 in an individual from Santarém_Rua_dos_Barcos_13th-c (PT_23227). This haplogroup is rare in Iberia but more commonly associated with Ashkenazi Jewish populations and Central Asia."

      Were there any sources used for this statement?

      As far as I know, Jewish clades under Q are typically not Q1b1a3a1/Q-L332 (clades under Q-L245 are more typical).

      Q-L332 is a Y-chromosome haplogroup found in several bronze age Siberian populations, and is a fairly prominent lineage across Scythian populations during the iron age and antiquity (hence my interest in the remark above).

      Given the history of the Iberian peninsula, would it not be more plausible to attribute such lineages to the Alans? Sarmatians carried Y-chromosome haplogroup Q-L332, and there seem to be a decent amount of Spaniards and Portuguese which carry Sarmatian-related lineages to this day. Perhaps the R1a-Z94 lineage could also be there through Sarmatians but it would depend on the subclade, as R-Z94 is over 4000 years old and was carried by various peoples - mostly of Indo-Iranian origin however.

      If there were any sources used regarding the Jewish origin of Q-L332 in the Iberian peninsula or perhaps archaeological signs which affirm the Jewish origin of sample PT_23227, I would be interested to have a look.

      Thanks in advance!

    1. On 2024-09-17 10:34:10, user balli wrote:

      The authors state in the intro... "only one clinical study has been published...."<br /> Please look at Jebsen et al 2019 J Med Case Rep / Spicer et al 2021Clin Cancer res.<br /> Also, using short peptides restric their use for intratumoral administration, please provide evidence.<br /> What about immunogenicity and potential ADA responses by longer peptides, please adress.<br /> Given that BOP peptides are strong inducers of ICD, why was immunodefect mouse model used ?

    1. On 2024-09-15 19:00:15, user Dr. M. K.Tiwari wrote:

      This mycobacterium may be similar to M tuberculosis, but it doesn't seem to be very unnatural as M.tuberculosis like organisms are reported from number of aquatic animals, engulfment of mycobacteria by amoeboid cells is well known and these mycobacteria are not digested in food vacuole of amoeba but some time there multiplication in vacuole is also reported.<br /> Sponge may have engulfed the mycobacteria from ablutions of an infected patient and entered in sponge along with incurrent water in spongocoel these mycobacteria were ingested by choanocyte where probably these bacteria were not digested and kept on multiplying. So it may be a possibility that this is not the mycobacteria originally from sponge but entered in sponge along with debris or ablutions of infected patient.

    1. On 2024-09-13 15:32:48, user Ibrahim, Tarhan E wrote:

      Chung et al. (2024) identified a physical interaction between ERC1 and ATG8e, leading them to explore potential ATG8-interaction motifs (AIMs) in ERC1. Using the iLIR database (Jacomin et al., 2016) for AIM prediction, they found the results irrelevant to the ERC1-ATG8e interaction, indicating a false prediction. Through truncated ERC1 variants, they identified a non-canonical AIM undetectable by current prediction tools, which focus on the canonical [W/F/Y]-[X]-[X]-[L/I/V] sequence. They validated this motif with AlphaFold2-multimer (AFM), a method we previously demonstrated (Ibrahim et al., 2023) to accurately predict non-canonical AIMs, as shown with ATG3. Our findings were later confirmed in humans by Farnung et al. (2023) via X-ray crystallography. Despite their similar approach, Chung et al. (2024) did not acknowledge our prior work

    1. On 2024-09-13 14:53:14, user Thibaud Decaens wrote:

      The manuscript has now been published in European Journal of Soil Biology:<br /> Gabriac Q., Ganault P., Barois I., Angeles G., Cortet J., Hedde M., Marchán D.F., Pimentel Reyes J.C., Stokes A., Decaëns T. (2023) Environmental drivers of earthworm communities along an altitudinal gradient in the French Alpes. European Journal of Soil Biology, 116, 103477. https://doi.org/10.1016/j.ejsobi.2023.103477

    1. On 2024-09-13 06:17:44, user Massimo Turina wrote:

      Interesting, but make sure you place it taxonomically in the new family "Konkoviridae".... therefore not really correct to call it a new "Phenuivirus"......but anyway good job.

    1. On 2024-09-12 09:32:02, user Fabienne Jabot-Hanin wrote:

      Could you please share your supplementary table 2 with the 385 index variants which had the same direction of effect on serum creatinine and cystatin C and had significant effects on both biomarkers ?

    1. On 2024-09-12 07:11:20, user Keshava Datta wrote:

      Great study - Extremely important to improve genome annotation as we know it... In one of the first drafts of the human proteome (Nature, 2014), ~16 million spectra that did not match to known proteome were subjected to proteogenomic analysis and ~200 regions with protein coding potential were found. It would have been great if the authors mentioned this and maybe compared these results? As we all know, evidence from multiple groups increases confidence in a finding!!

      (Full disclosure - I was a co-author on the 2014 paper)..

    1. On 2024-09-08 19:36:57, user Cara J. Gottardi wrote:

      Can the authors please confirm use of recombinant human WNT2 from Novus Biologicals (H00007472-P01) for their rescue experiments? The supplier says this protein is not designed to be active, and should not be used for activity-based assays (e.g., the protein is GST-tagged, not ideal for WNT proteins; also wheat germ systems do not allow for glycosylation of secreted proteins). Happy to be wrong if this protein prep really works!

    1. On 2024-09-07 13:32:35, user D_114 wrote:

      The work references Wilkes (2021) as mentioning 80% of glomalin is located within the hyphal network. However, on cross referencing this statement, the article by Wilkes (2021) makes no such claim. Therefore, the statement in this article is misrepresenting research and misleading readers.

    1. On 2024-09-05 16:30:46, user Paolo Ubezio wrote:

      This is a preprint of the following chapter: Ubezio, P., Challenging Age-Structured and First Order Transition Cell Cycle Models of Cell Proliferation, published in Problems in Mathematical Biophysics. SEMA SIMAI Springer Series, vol 38, edited by d'Onofrio, A., Fasano, A., Papa, F., Sinisgalli, C., 2024, Springer, reproduced with permission of Springer Nature Switzerland AG. The final revised and authenticated version is available online at: http://dx.doi.org/10.1007/978-3-031-60773-8_13 .

    1. On 2024-09-05 01:22:44, user GR wrote:

      Hi Authors

      Very interesting paper, I just noticed some small potential mistakes when reading. In Fig 3E, it looks like the beta-tubulin image has been put in place of gamma-tubulin, and in Figure 5A, the GO treatment image appears to be the same as the Figure 4D AA image. Hopefully this comment is helpful!

    1. On 2024-09-04 21:15:50, user Alex Grossman wrote:

      Hello, as a heads up you call Xenorhabdus Gram-positive in your abstract. It appears to be correctly listed as Gram-negative elsewhere in the text.

    1. On 2024-09-03 22:48:26, user Pooja Asthana wrote:

      Summary<br /> The study investigates the human protein DJ-1, which is known for its role in detoxifying the metabolic bioproduct methylglyoxal (MG). There has been an ongoing debate over whether DJ-1 acts directly on MG (direct substrate) or requires a protein intermediate acting as a protein/nucleic acid deglycase (glycated protein substrate). The authors used fixed-target micro-crystallography and mix-and-inject serial crystallography to structurally analyze covalent intermediates in the reaction catalyzed by DJ-1. One of the significant achievements of the study is the successful use of these advanced crystallography methods to determine the structure of key reaction intermediates: hemithioacetal and L-lactoylcysteine, providing new insights into DJ-1's glyoxalase mechanism. However, a major weakness is that the authors' claim refuting the alternative deglycase mechanism are not fully supported by the presented data. Despite this limitation, the study advances our understanding of DJ-1’s enzymatic function by leveraging MISC at synchrotron using the new flow cell injector.

      Major points<br /> Major point 1<br /> The claim made in the discussion that: “These results provide direct structural evidence supporting a growing number of enzymology studies also indicating that DJ-1 is not a deglycase…” is not supported by evidence presented in the manuscript. Although this work elegantly demonstrates that MG covalently modifies the catalytic cysteine of DJ-1 (Cys106), the crystallography experiments presented are unable to test whether the alternative mechanism (with a glycated substrate) occurs. More careful treatment of this logic in the discussion would strengthen the manuscript, and would help the manuscript to be more focused on the compelling X-ray crystallography results. We recognize it is difficult to “prove a negative” however these experiments affirm the primary activity without directly testing the alternative one.

      Major point 2<br /> The authors report compelling evidence that the DJ-1 catalytic cysteine (Cys106) is covalently modified by MG. However, the concentration of MG used was 50 mM, and non-catalytic cysteines might be covalently modified at this concentration. Indeed, it’s possible that one of the DJ-1 surface cysteines is covalently modified (Cys53), based on the large positive difference peak in the FO-FO difference density (Figure 5b, Figure S8) (although it is suggested that this is evidence of allosteric communication). Covalent modification of a surface cysteine leading to lattice disruption is consistent with the observation that MG is known to dissolve DJ-1 crystals. The manuscript could be strengthened by consideration of these points, as well as analysis of difference maps around Cys53 for the fixed target structure (e.g. add panel to Figure S1 showing FO(methylglyoxal)-FO(free) maps around Cys53). Discussion of the differences in modification rates for the catalytic and surface cysteines, and the impact of large versus small crystals, would be helpful.

      Major point 3<br /> Is it plausible that a second, synchronized turnover is captured by the mix-and-inject experiment? This claim might be developed by modeling the concentration profile of the intermediates along the 30 second time course (e.g. similar to Figure 4 in PMID 29848358). To this point, were the occupancies of the covalent adducts refined at each time point? Did the authors consider whether a mixture of species might be present? The evidence supporting the second turnover comes from the featureless difference map calculated between the 3 sec and 20 sec time points (FO(20s)-FO(3s) in Figure S6). Is there an alternative explanation for the decreased occupancy at this time point other than synchronized turnover? E.g. a problem with sample mixing resulting in lower substrate concentration at this time point.

      A related concern is whether the data as presented can discriminate between the two covalent intermediates (HTA or LC). Perhaps Figure S7 would be strengthened by adding the FO-FC difference maps for each of the intermediates modeled with the other species (e.g. the HTA dataset modeled with LC and vice versa). Can the authors comment on the lack of correlated negative (or positive) density in the FO-FO difference map matrix (Figure S5) in panels comparing sp2 and sp3 carbons (e.g. FO(15s)-FO(3s)). In this example, there is a large positive peak in the difference map for the sp2 to sp3 change, but no correlated negative peak.

      Minor points<br /> Minor point 1<br /> Was the covalent adduct observed in the MG-soaked DJ-1 crystals presented in Figure S1c modeled? Is the difference density consistent with the HTA or LC intermediates? Or a mixture of both?<br /> Minor point 2<br /> Is it possible that movement of the active site histidine (His126) away from covalent intermediate (Figure 4a) is consistent with histidine protonation? Or is the geometry such that protonation is unlikely?<br /> Minor point 3<br /> We find it helpful if the figure (or figure legend) includes PDB codes for their quick look up.<br /> Minor point 4<br /> The size of the scale bar in Figure S1a might be increased.

      Review by:<br /> Pooja Asthana, Galen J. Correy & James S. Fraser (UCSF)

    1. On 2024-09-02 14:42:35, user Peer wrote:

      Looks like a very interesting story. Although, since the authors did not include any methodology it's hard to evaluate the quality and reproducibility of this work.

    1. On 2024-09-02 13:32:22, user Adriano R. Lameira wrote:

      What an amazing behaviour, thank you for bringing it to light!<br /> For your interest and to help keep the paper as accurate as possible, here are two former references on movement isochrony in apes:<br /> https://www.nature.com/articles/s41598-019-55360-y <br /> https://www.cell.com/current-biology/fulltext/S0960-9822(22)01601-3 <br /> Hopefully these two papers, together with the current contribution, will help bring some phylogenetic insight into discussions of dance evolution.<br /> Looking forward to the publication of this study, congratulations on this important work.

    1. On 2024-09-01 06:48:28, user Giorgio Cattoretti wrote:

      Dear Colleagues,<br /> I read with interest the manuscript you posted about the Cyclic-multiplex TSA (CmTSA) method.<br /> In the Result subchapter titled “Gentle antibody stripping is necessary for sequential immunostaining with scores of iterations” at line 171 and following you mention several potential Ab stripping reagents, which are not further detailed in the M&M section. Furthermore, you mention beta-ME (line 173) and a low Ph buffer (line 175) as two separate stripping solutions. In the text, there is only one proprietary buffer, IRISKit® HyperView Advanced Ab Stripping Kit , for which a webpage in chinese can be found ( http://m.luminiris.cn/pd.jsp?pid=14):gTX3SSR1ouS7_HrdDaiJq0-yQz4 "http://m.luminiris.cn/pd.jsp?pid=14)") and no information on the composition can be obtained, neither you provide in the manuscript.<br /> To my knowledge, there is no mention in the scientific peer-reviewed literature of beta-ME ALONE as a stripping buffer. Besides the inventor of the 2-ME containing buffer, Laemmli (Nature. 1970;227:680–685), the earliest mention of beta-ME for stripping is Pirici (doi: 10.1369/jhc.2009.953240), after which we did quite some experiments (doi: 10.1369/0022155414536732; doi: 10.1097/pai.0000000000000203; doi: 10.1369/0022155417719419, this latter viewed 23,855 times and quoted in 101 papers). None of these references are among the one you quote.<br /> For reproducibility and to avoid re-inventing the wheel for one’s own advantage, you may wish to quote the relevant previous art and provide evidence of any modification or improvement. If the solution is proprietary, full disclosure is require or mention that the content of the proprietary solution is unknown.<br /> I do appreciate the effort put in the manuscript, which clearly appears an immature draft. Looking forward for improvements.

    1. On 2024-08-30 17:42:18, user David Černý wrote:

      The in-text citations are not properly linked to the corresponding bibliography entries. It would be great to post a v2 that fixes this problem.

    1. On 2024-08-30 12:21:52, user Justin Lemkul wrote:

      An interesting study and nice demonstration of the AMOEBA force field. I would advise the authors to edit the statement "there are no precedents to our knowledge to MD simulations on any HIV-1 forming G-quadruplex" as there have been several studies on LTR GQs in the literature. I suggest citing the following and incorporating them into the discussion:

      https://doi.org/10.1016/j.bpj.2024.03.042 <br /> https://pubs.acs.org/doi/10.1021/acs.jctc.2c00291

    1. On 2024-08-23 06:58:52, user viroid wrote:

      Fantastic work! Integrating your domain meta-clusters (or maybe just the orphaned ones) into the Pfam would be super helpful! A pipeline of FoldMason and HMMer3 shouldn't be too tough?

    1. On 2024-08-22 15:51:55, user Jack M wrote:

      I have two comments:

      1. A recent paper claims that "certain cancers can use ketone bodies as an alternative energy source to sustain tumor fitness, such as pancreatic tumors". The results in this paper, even without a glutamine inhibitor, would disagree with that statement. Maybe you can mention the conflicting approaches to ketogenic diet in the discussion?

      2. Have you looked into the work of Mukherjee et al ? They report a very similar effect with ketogenic diet and DON in late stage experimental glioblastoma. Perhaps this intervention could be a more universal therapy?

      Thank you for this very promising work exploring diet-drug combinations!

    1. On 2024-08-21 16:52:34, user DUPA- Preprint Review wrote:

      The study titled “Temporal dynamics of BMP/Nodal ratio drive tissue-specific gastrulation morphogenesis” investigates how Nodal signaling dynamics and its interactions with BMP signaling affect convergence and extension (C&E) in zebrafish. Emig et.al . employs a combination of methods, including a novel tool for manipulating Nodal activation timing with optogenetics, gene expression profiles, and tissue-specific transplants in embryos.

      This work shows a critical role for the temporal ratio of BMP to Nodal signaling in establishing tissue-specific C&E. High Nodal signaling at only before 4-hours post fertilization promotes mesoderm-driven extension, while a high BMP/Nodal ratio during this time window promotes extension driven by the neural ectoderm. The study efficiently utilizes explants to dissect these signaling pathways, and then supports the findings by demonstrating similar effects in intact embryos. However, the work would be strengthened by a better characterization of the BMP/Nodal ratio and downstream targets in their experimental manipulations.

      The study elucidates how different methods and timings of Nodal signaling activation induce distinct modes of tissue-specific C&E in zebrafish. The interplay between Nodal and BMP signaling, particularly the timing and ratio of these signals, is crucial in determining whether mesoderm or NE-driven morphogenesis occurs. These findings enhance our understanding of early developmental processes. The paper presents compelling evidence on the distinct C&E programs within the mesoderm and ectoderm regulated by BMP/Nodal signaling ratios during a key developmental window.

      Overall, this study provides valuable insights into the interplay between Nodal and BMP signaling in shaping early embryonic development.

      Strengths:

      Clear and Robust Experimental Design: The use of transgenic embryos, mRNA injection, and gene expression analysis provides a strong foundation for the research. The morphological observations coupled with gene expression analysis strengthens the conclusions. The uninjected controls and different concentrations of mRNA were critical in showing the threshold of Nodal signaling. Live imaging provides valuable, real-time data on the cell-autonomous effects of BMP on the ectoderm.

      Time-Sensitive and Tissue-Specific Manipulation: The use of validated optogenetic tools allowed precise control of Nodal signaling onset. The optogenetic BMP tool allowed for precise control over the timing and location of BMP activity. This revealed a critical window for its influence and highlights the concept of BMP/Nodal ratio as a key factor. Furthermore, the tissue-specific activation of BMP in either mesoderm or ectoderm reinforced the link between BMP signaling and distinct tissue behaviors.

      Points of consideration:

      Measuring the BMP/Nodal Ratio: The BMP/Nodal ratio was measured over time using transcriptional profiling in the original explant models. The manipulations that the authors perform in the rest of the paper are meant to alter the BMP/Nodal dynamics throughout this time course. However, it is unclear to the reviewers if the BMP/Nodal ratio was measured in these manipulations. It would be nice to detail the dynamics of this ratio with data rather than the schematics in each figure. The authors could do this with transcriptional profiling or staining for SMAD activity as shown in Figure 2. In addition, staining for SMAD activity would perhaps be more indicative of actual signaling ratios rather than the transcriptional products, which might have multiple inputs.

    1. On 2024-08-21 16:31:36, user DUPA- Preprint Review wrote:

      The manuscript by Siyi Gu and colleagues presents an unbiased approach using well-established APEX2 proximity labeling proteomics and targeted pharmacological experiments that demonstrate different CCR5 chemokine receptor ligands-based induction of distinct receptor signaling responses and trafficking behaviors, including intracellular receptor sequestration, which offers a potential therapeutic strategy for inhibiting CCR5 functions a repertoire of CCR5 functional diversity. The study reveals the molecular basis for receptor sequestration, including information that can be exploited to develop actionable patterns for developing chemokine-based CCR5 targeting molecules that promote retention of the receptor inside the cell. This is particularly noteworthy because CCR5 plays a crucial role in the immune system and is important in numerous physiological and pathological processes such as inflammation, cancer, and HIV transmission. The work has great scope and importance as an alternative therapeutics strategy to address inflammation, cancer, and HIV transmission, and the experiments are well-designed and sound. The manuscript is well-written, clear, and has a reasonable and logical flow.<br /> We have some minor comments and a few methodological suggestions that would improve the quality of the manuscript.

      1. The authors have used the widely accepted APEX proximity labeling technique; however, they have not included a mock or vehicle-treated control to accurately examine the interacting proteins in the different ligand stimulated conditions. This would compensate for any proteins that interact with CCR5 in the absence of a ligand.
      2. The method section sometimes lacks critical information (please read all and add relevant details when needed, which is important for the nonexperts). For example, Fig 1 has missing details on the negative control or vehicle control to compare time-resolved ligand-dependent trafficking along with some agonists, partial agonists, antagonists, and inverse agonists.
      3. The effect of the different ligands on CCR5 functionality and trafficking should be assessed in the physiologically relevant cell line to determine whether the findings from this paper hold true and can be used for therapeutics. Often, using overexpression systems leads to the observation of phenomena that are absent when the receptor is expressed at native levels.
      4. The degradation assays to show CCR5 localization to the lysosome are interesting and relevant but it would be useful to see the whole blot so that the reader can view these as low molecular weight products.
      5. Suppose authors can add time-resolved live confocal microscopy for the trafficking of CCR5 under different treatment conditions and show the colocalization. This would add to the impact of the study.
    1. On 2024-08-19 13:52:54, user Jonathan Rondeau-Leclaire wrote:

      Dear authors, <br /> This study is promising, as you seem to have gathered quality data with a very interesting design that has great potential to generate insights into the impact of microbes on plant productivitiy. I must however venture in a technical comment, as I believe the statistical approach you have chosen is weak and prevents you from truly leveraging all the precious information you have generated with the sequencing experiments.

      Using correlation to find associations between microbial abundance and environmental (or sample) characteristics using data derived from sequencing experiment is generally advised against, as the data is compositional, which means the values are set in a simplex, not a euclidean space. This means that the observed relative abundance values are not independent of each other, as there is an inherent correlation between all taxa: if any microbe increases its true abundance, the relative abundance of everyone else will decrease even if they did not change in true abundance.

      To find bacterial genera associated with plant traits (or any other characteristic), you should use statistical methods developed specifically for handling microbial relative abundance data. You can look up ANCOM-BC, corncob, DESeq, Aldex, and many others that have been developed to work specifically with sequencing data. Some of these even estimate the changes in absolute abundances. There are other reasons to use these methods too, that have to do with special characteristics of sequencing data (which can hardly be ignored), such as sparseness, overdispersion, to name a few. I recommend the following reads:<br /> https://academic.oup.com/gigascience/article/doi/10.1093/gigascience/giz107/5572529 <br /> https://dx.plos.org/10.1371/journal.pcbi.1010467 <br /> https://www.nature.com/articles/s41522-020-00160-w

      Moreover, I do not see any mention of multiple testing correction. As you test multiple genera, it is absolutely essential to correct your p-values for multiple testing, otherwise it is almost certain that some of the genera you identified as significantly correlated were only by pure chance, not for biological reasons. Most differential abundance tests mentioned above do this by default, as it is expected for credible results whenever conducting multiple statistical tests. More on this if you are not familiar with this correction: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6099145/

      I hope this helps, and good luck with your publication!

    1. On 2024-08-19 04:14:33, user Gray Shaw wrote:

      Very interesting observations. Any experiments done as yet to see if a blocking anti-VISTA antibody reduces the luminescence signal generated by the cis interaction of VISTA-SmBiT and PSGL-1-LgBiT at pH6?

    1. On 2024-08-18 14:42:36, user David Ron wrote:

      That mutations abolishing kinase activity (e.g., GCN2 Lys619Ala) also render RAF inhibitors unable to activate the ISR in cells is a powerful argument that GCN2's kinase activity links application of RAF inhibitors to ISR activation.<br /> That the same mutations also abolish GCN2's responsiveness to RAF inhibitors in vitro is less helpful - a dead donkey does not respond to the buggy whip, whether applies its back or to its head. <br /> In this vein, interpreting the inability of RAF inhibitors to activate GCN2 when the gatekeeper residue Met802 is mutated (Figure 9) is key to the conclusion that these agents exert their effect on GCN2 by engaging its kinase active site ATP binding pocket.<br /> I may have missed it, but is there evidence that the Met802 mutations retain reasonable kinase activity? The Met802Phe is not so helpful in this regard as it may lack 'headroom' for further activation.<br /> It might be helpful to measure ISR activation in GCN2∆ cells, reconstituted with GCN2 Met802 mutants in response to histidinol. Preservation of a response to histidinol in face of loss of responsiveness to RAF inhibitors would lend strong support to the authors conclusion.

    1. On 2024-08-12 15:52:52, user Swathi Sukumar wrote:

      Are the IF images in Figure 7, Z-stacks? Max projection of Z-stacks and more technical replicates are required. Additionally, how were the puncta counted?

    1. On 2024-08-12 09:58:10, user John Hawks wrote:

      The second part of the summary of revisions follows below:

      Recognition of pits

      Referee 3 and 4 and several additional commentators have emphasized that the recognition of pit features is necessary to the hypothesis of burial, and questioned whether the data presented in the manuscript were sufficient to demonstrate that pits were present. We have revised the manuscript in several ways to clarify how all the different kinds of evidence from the subsystem test the hypothesis that pits were present. This includes the presentation of a minimal definition of burial to include a pit dug by hominins, criteria for recognizing that a pit was present, and an evaluation of the evidence in each case to make clear how the evidence relates to the presence of a pit and subsequent infill. As referee 3 notes, it can be challenging to recognize a pit when sediment is relatively homogeneous. This point was emphasized in the review by Pomeroy and coworkers (2020b), who reflected on the difficulty seeing evidence for shallow pits constructed by hominins, and we have cited this in the main text. As a result, the evidence for pits has been a recurrent topic of debate for most Pleistocene burial sites. However in addition to the sedimentological and contextual evidence in the cases we describe, the current version also reflects upon other possible mechanisms for the accumulation of bones or bodies. The data show that the sedimentary fill associated with the H. naledi remains in the cases we examine could not have passively accumulated slowly and is not indicative of mass movement by slumping or other high-energy flow. To further put these results into context, we added a section to the Discussion that briefly reviews prior work on distinguishing pits in Pleistocene burial contexts, including the substantial number of sites with accepted burial evidence for which no evidence of a pit is present.

      Extent of articulation and anatomical association

      We have added significantly greater detail to the descriptions of articulated remains and orientation of remains in order to describe more specifically the configuration of the skeletal material. We also provide 14 figures in main text (13 of them new) to illustrate the configuration of skeletal remains in our data. For the Puzzle Box area, this now includes substantial evidence on the individuation of skeletal fragments, which enables us to illustrate the spatial configuration of remains associated with the DH7 partial skeleton, as well as the spatial position of fragments refitted as part of the DH1, DH2, DH3, and DH4 crania. For Dinaledi Feature 1 and the Hill Antechamber Feature we now provide figures that key skeletal parts as identified, including material that is unexcavated where possible, and a skeletal part representation figure for elements excavated from Dinaledi Feature 1.

      Archaeothanatology

      Reviewer 2 suggests that a greater focus on the archaeothanatology literature would be helpful to the analysis, with specific reference to the sequence of joint disarticulation, the collapse of sediment and remains into voids created by decomposition, and associated fragmentation of the remains. In the revised manuscript we have provided additional analysis of the Hill Antechamber Feature with this approach in mind. This includes greater detail and illustration of our current hypothesis for individuation of elements. We now discuss a hypothesis of body disposition, describe the persistent joints and articulation of elements, and examine likely decomposition scenarios associated with these remains. Additionally, we expand our description and illustration of the orientation of remains and degree of anatomical association and articulation within Dinaledi Feature 1. For this feature and for the Hill Antechamber Feature we have revised the text to describe how fracturing and crushing patterns are consistent with downward pressure from overlying sediment and material. In these features, postdepositional fracturing occurred subsequent to the decomposition of soft tissue and partial loss of organic integrity of the bone. We also indicate that the loss by postdepositional processes of most long bone epiphyses, vertebral bodies, and other portions of the skeleton less rich in cortical bone, poses a challenge for testing the anatomical associations of the remaining elements. This is a primary reason why we have taken a conservative approach to identification of elements and possible associations.

      A further aspect of the site revealed by our analysis is the selective reworking of sediments within the Puzzle Box area subsequent to the primary deposition of some bodies. The skeletal evidence from this area includes body parts with elements in anatomical association or articulation, juxtaposed closely with bone fragments at varied pitch and orientation. This complexity of events evidenced within this area is a challenge for approaches that have been developed primarily based on comparative data from single-burial situations. In these discussions we deepen our use of references as suggested by the referee.

      Burial positions

      Reviewer 2 further suggests that illustrations of hypothesized burial positions would be valuable. We recognize that a hypothesized burial position may be an appealing illustration, and that some recent studies have created such illustrations in the context of their scientific articles. However such illustrations generally include a great deal of speculation and artist imagination, and tend to have an emotive character. We have added more discussion to the manuscript of possible primary disposition in the case of the Hill Antechamber Feature as discussed above. We have not created new illustrations of hypothesized burial positions for this revision.

      Carnivore involvement

      Referee 1 suggests that the manuscript should provide further consideration of whether carnivore activity may have introduced bones or bodies into the cave system. The reorganized Introduction now includes a review of previous work, and an expanded discussion within the Supplementary Information (“Hypotheses tested in previous work”). This includes a review of literature on the topic of carnivore accumulation and the evidence from the Dinaledi and Lesedi Chamber that rejects this hypothesis.

      Water transport and mud

      The eLife referees broadly accepted previous work showing that water inundation or mass flow of water-saturated sediment did not occur within the history of Unit 2 and 3 sediments, including those associated with H. naledi remains. However several outside commentators did refer specifically to water flow or mud flow as a mechanism for slumping of deposits and possible sedimentary covering of the remains. To address these comments we have added a section to the Supplementary Information (“Description of the sedimentary deposits of the Dinaledi Subsystem”) that reviews previous work on the sedimentary units and formation processes documented in this area. We also include a subsection specifically discussing the term “mud” as used in the description of the sedimentology within the system, as this term has clearly been confusing for nonspecialists who have read and commented on the work. We appreciate the referees’ attention to the previous work and its terminology.

      Redescription of areas of the cave system

      Reviewer 1 suggests that a detailed reanalysis of all portions of the cave system in and around the Dinaledi Subsystem is warranted to reject the hypothesis that bodies entered the space passively and were scattered from the floor by natural (i.e. noncultural) processes. The referee suggests that National Geographic could help us with these efforts. To address this comment we have made several changes to the manuscript. As noted above, we have added material in Supplementary Information to review the geochronology of the Dinaledi Subsystem and nearby Dragon’s Back Chamber, together with a discussion of the connections between these spaces.

      Most directly in response to this comment we provide additional documentation of the possibility of movement of bodies or body parts by gravity within the subsystem itself. This includes detailed floor maps based on photogrammetry and LIDAR measurement, where these are physically possible, presented in Figures 2 and 3. In some parts of the subsystem the necessary equipment cannot be used due to the extremely confined spaces, and for these areas our maps are based on traditional survey methods. In addition to plan maps we have included a figure showing the elevation of the subsystem floor in a cross-section that includes key excavation areas, showing their relative elevation. All figures that illustrate excavation areas are now keyed to their location with reference to a subsystem plan. These data have been provided in previous publications but the visualization in the revised manuscript should make the relationship of areas clear for readers. The Introduction now includes text that discusses the configuration of the Hill Antechamber, Dinaledi Chamber, and nearby areas, and also discusses the instances in which gravity-driven movement may be possible, at the same time reviewing that gravity-driven movement from the entry point of the subsystem to most of the localities with hominin skeletal remains is not possible.

      Within the Results, we have added a section on the relationship of features to their surroundings in order to assist readers in understanding the context of these bone-bearing areas and the evidence this context brings to the hypothesis in question. We have also included within this new section a discussion of the discrete nature of these features, a question that has been raised by outside commentators.

      Passive sedimentation upon a cave floor or within a natural depression

      Reviewer 3 suggests that the situation in the Dinaledi Subsystem may be similar to a European cave where a cave bear skeleton might remain articulated on a cave floor (or we can add, within a hollow for hibernation), later to be covered in sediment. The reviewer suggests that articulation is therefore no evidence of burial, and suggests that further documentation of disarticulation processes is essential to demonstrating the processes that buried the remains. We concur that articulation by itself is not sufficient evidence of cultural burial. To address this comment we have included a section in the Results that tests the hypothesis that bodies were exposed upon the cave floor or within a natural depression. To a considerable degree, additional data about disarticulation processes subsequent to deposition are provided in our reanalysis of the Puzzle Box area, including evidence for selective reworking of material after burial.

      Postdepositional movement and floor drains

      Reviewer 3 notes that previous work has suggested that subsurface floor drains may have caused some postdepositional movement of skeletal remains. The hypothesis of postdepositional slumping or downslope movement has also been discussed by some external commentators (including Martinón-Torres et al. 2024). We have addressed this question in several places within the revised manuscript. As we now review, previous discussion of floor drains attempted to explain the subvertical orientation of many skeletal elements excavated from the Puzzle Box area. The arrangement of these bones reflects reworking as described in our previous work, and without considering the possibility of reworking by hominins, one mechanism that conceivably might cause reworking was downward movement of sediments into subsurface drains. Further exploration and mapping, combined with additional excavation into the sediments beneath the Puzzle Box area provided more information relevant to this hypothesis. In particular this evidence shows that subsurface drains cannot explain the arrangement of skeletal material observed within the Puzzle Box area. As now discussed in the text, the reworking is selective and initiated from above rather than below. This is best explained by hominin activity subsequent to burial.

      In a new section of the Results we discuss slumping as a hypothesis for the deposition of the remains. This includes discussion of downslope movement within the Hill Antechamber and the idea that floor drains may have been a mechanism for sediment reworking in and around the Puzzle Box area and Dinaledi Feature 1. As described in this section the evidence does not support these hypotheses.

      Hypothesis testing and parsimony

      Referees 1 and 3 and the editorial guidance all suggested that a more appropriate presentation would adopt a null hypothesis and test it. The specific suggestion that the null hypothesis should be a natural sedimentary process of deposition was provided not only by these reviewers but also by some outside commentators. To address this comment, we have edited the manuscript in two ways. The first is the addition of a section to the Discussion that specifically discusses hypothesis testing and parsimony as related to Pleistocene evidence of cultural burial. This includes a brief synopsis of recent disciplinary conversations and citation of work by other groups of authors, none of whom adopted this “null hypothesis” approach in their published work.

      As we now describe in the manuscript, previous work on the Dinaledi evidence never assumed any role for H. naledi in the burial of remains. Reading the reviewer reports caused us to realize that this previous work had followed exactly the “null hypothesis” approach that some suggested we follow. By following this null hypothesis approach, we neglected a valuable avenue of investigation. In retrospect, we see how this approach impeded us from understanding the pattern of evidence within the Puzzle Box area. Thus in the revised manuscript we have mentioned this history within the Discussion and also presented more of the background to our previous work in the Introduction. Hopefully by including this discussion of these issues, the manuscript will broaden conversation about the relation of parsimony to these issues.

      Language and presentation style

      Reviewer 4 criticizes our presentation, suggesting that the text “gives the impression that a hypothesis was formulated before data were collected.” Other outside commentators have mentioned this notion also, including Martinón-Torres et al. (2024) who suggest that the study began from a preferred hypothesis and gathered data to support it. The accurate communication of results and hypotheses in a scientific article is a broader issue than this one study. Preferences about presentation style vary across fields of study as well as across languages. We do not regret using plain language where possible. In any study that combines data and methods from different scientific disciplines, the use of plain language is particularly important to avoid misunderstandings where terms may mean different things in different fields.

      The essential question raised by these comments is whether it is appropriate to present the results of a study in terms of the hypothesis that is best supported. As noted above, we read carefully many recent studies of Pleistocene burial evidence. We note that in each of these studies that concluded that burial is the best hypothesis, the authors framed their results in the same way as our previous manuscript: an introduction that briefly reviews background evidence for treatment of the dead, a presentation of results focused on how each analysis supports the hypothesis of burial for the case, and then in some (but not all) cases discussion of why some alternative hypotheses could be rejected. We do not infer from this that these other studies started from a presupposition and collected data only to confirm it. Rather, this is a simple matter of presentation style.

      The alternative to this approach is to present an exhaustive list of possible hypotheses and to describe how the data relate to each of them, at the end selecting the best. This is the approach that we have followed in the revised manuscript, as described above under the direction of the reviewer and editorial guidance. This approach has the advantage of bringing together evidence in different combinations to show how each data point rejects some hypotheses while supporting others. It has the disadvantage of length and repetition.

      Possible artifact

      We have chosen to keep the description of the possible artifact associated with the Hill Antechamber Feature in the Supplementary Information. We do this while acknowledging that this is against the opinion of reviewer 4, who felt the description should be removed unless the object in question is fully excavated and physically analyzed. The previous version of the manuscript did not rely upon the stone as positive evidence of grave goods or symbolic content, and it noted that the data do not test whether the possible artifact was placed or was intentionally modified. However this did not satisfy reviewer 4, and some outside commentators likewise asserted that the object must be a “geofact” and that it should be removed.

      We have three arguments against this line of thinking. First, we do not omit data from our reporting. Whether Homo naledi shaped the rock or not, used it as a tool or not, whether the rock was placed with the body or not, it is unquestionably there. Omitting this one object from the report would be simply dishonest. Second, the data on this rock are at 16 micron resolution. While physical inspection of its surface may eventually reveal trace evidence and will enable better characterization of the raw material, no mode of surface scanning will produce better evidence about the object’s shape. Third, the position of this possible artifact within the feature provides significant information about the deposition of the skeletal material and associated sediments. The pitch, orientation, and position of the stone is not consistent with slow deposition but are consistent with the hypothesis that the surrounding sediment was rapidly emplaced at the same time as the articulated elements less than 2 cm away.

      In the current version, we have redoubled our efforts to provide information about the position and shape of this stone while not presupposing the intentionality of its shape or placement. We add here that the attitude expressed by referee 4 and other commentators, if followed at other sites, would certainly lead to the loss or underreporting of evidence, which we are trying to avoid.

      Consistency versus variability of behavior

      As described in the revised manuscript, different features within the Dinaledi Subsystem exhibit some shared characteristics. At the same time, they vary in positioning, representation of individuals and extent of commingling. Other localities within the subsystem and broader cave system present different evidence. Some commentators have questioned whether the patterning is consistent with a single common explanation, or whether multiple explanations are necessary. To address this line of questioning, we have added several elements to the manuscript. We created a new section on secondary cultural burial, discussing whether any of the situations may reflect this practice. In the Discussion, we briefly review the ways in which the different features support the involvement of H. naledi without interpreting anything about the intentionality or meaning of the behavior. We further added a section to the Discussion to consider whether variation among the features reflects variation in mortuary practices by H. naledi. One aspect of this section briefly cites variation in the location and treatment of skeletal remains at other sites with evidence of burial.

      Grave goods

      Some commentators have argued that grave goods are a necessary criterion for recognizing evidence of ancient burial. We added a section to the Discussion to review evidence of grave goods at other Pleistocene sites where burial is accepted.

    2. On 2024-08-12 09:57:17, user John Hawks wrote:

      The authors have prepared a summary of revisions to this manuscript, which is available in the Supplementary Information section of the preprint. I have included the text of that summary across the following two comments in two parts (the whole is larger than the maximum comment length).

      Summary of revisions to Evidence for deliberate burial of the dead by Homo naledi

      We extend our sincere thanks to the editor, referees for eLife, and other commentators who have written evaluations of this manuscript, either in whole or in part. Sources of these comments were highly varied, including within the bioRxiv preprint server, social media (including many comments received on X/Twitter and some YouTube presentations and interviews), comments made by colleagues to journalists, and also some reviews of the work published in other academic journals. Some of these are formal and referenced with citations. Others were informal but nonetheless expressed perspectives that helped enable us to revise the manuscript with the inclusion of broader perspectives than the formal review process. It is beyond the scope of this summary to list every one of these, which have often been brought to the attention of different coauthors, but we begin by acknowledging the very wide array of peer and public commentary that have contributed to this work. The reaction speaks to a broad interest in open discussion and review of preprints.

      As we compiled this summary of changes to the manuscript, we recognized that many colleagues made comments about the process of preprint dissemination and evaluation rather than the data or analyses in the manuscript. Addressing such comments is outside the scope of this revised manuscript, but we do feel that a broader discussion of these comments would be valuable in another venue. Many commentators have expressed confusion about the eLife system of evaluation of preprints, which differs from the editorial acceptance or rejection practiced in most academic journals. As authors in many different nations, in varied fields, and in varied career stages, we ourselves are still working to understand how the academic publication landscape is changing, and how best to prepare work for new models of evaluation and dissemination.

      The manuscript and coauthor list reflect an interdisciplinary collaboration. Analyses presented in the manuscript come from a wide range of scientific disciplines. These range from skeletal inventory, morphology, and description, spatial taphonomy, analysis of bone fracture patterns and bone surface modifications, sedimentology, geochemistry, and traditional survey and mapping. The manuscript additionally draws upon a large number of previous studies of the Rising Star cave system and the Dinaledi Subsystem, which have shaped our current work. No analysis within any one area of research stands alone within this body of work: all are interpreted in conjunction with the outcomes of other analyses and data from other areas of research. Any single analysis in isolation might be consistent with many different hypotheses for the formation of sediments and disposition of the skeletal remains. But testing a hypothesis requires considering all data in combination and not leaving out data that do not fit the hypothesis. We highlight this general principle at the outset because a number of the comments from referees and outside specialists have presented alternative hypotheses that may arguably be consistent with one kind of analysis that we have presented, while seeming to overlook other analyses, data, or previous work that exclude these alternatives. In our revision, we have expanded all sections describing results to consider not only the results of each analysis, but how the combination of data from different kinds of analysis relate to hypotheses for the deposition and subsequent history of the Homo naledi remains. We address some specific examples and how we have responded to these in our summary of changes below.

      General organization

      The referee and editor comments are mostly general and not line-by-line questions, and we have compiled them and treated them as a group in this summary of changes, except where specifically noted.

      The editorial comments on the previous version included the suggestion that the manuscript should be reorganized to test “natural” (i.e. noncultural) hypotheses for the situations that we examine. The editorial comment suggested this as a “null hypothesis” testing approach. Some outside comments also viewed noncultural deposition as a null hypothesis to be rejected. We do not concur that noncultural processes should be construed as a null hypothesis, as we discuss further below. However, because of the clear editorial opinion we elected to revise the manuscript to make more explicit how the data and analyses test noncultural depositional hypotheses first, followed by testing of cultural hypotheses. This reorganization means that the revised manuscript now examines each hypothesis separately in turn.

      Taking this approach resulted in a substantial reorganization of the “Results” section of the manuscript. The “Results” section now begins with summaries of analyses and data conducted on material from each excavation area. After the presentation of data and analyses from each area, we then present a separate section for each of several hypotheses for the disposition and sedimentary context of the remains. These hypotheses include deposition of bodies upon a talus (as hypothesized in some previous work), slow sedimentary burial on a cave floor or within a natural depression, rapid burial by gravity-driven slumping, and burial of naturally mummified remains. We then include sections to test the hypothesis of primary cultural burial and secondary cultural burial. This approach adds substantial length to the Results. While some elements may be repeated across sections, we do consider the new version to be easier to take piece by piece for a reader trying to understand how each hypothesis relates to the evidence.

      The Results section includes analyses on several different excavation areas within the Dinaledi Subsystem. Each of these presents somewhat different patterns of data. We conceived of this manuscript combining these distinct areas because each of them provides information about the formation history of the Homo naledi-associated sediments and the deposition of the Homo naledi remains. Together they speak more strongly than separately. In the previous version of the manuscript, two areas of excavation were considered in detail (Dinaledi Feature 1 and the Hill Antechamber Feature), with a third area (the Puzzle Box area) included only in the Discussion and with reference to prior work. We now describe the new work undertaken after the 2013-2014 excavations in more detail. This includes an overview of areas in the Hill Antechamber and Dinaledi Chamber that have not yielded substantial H. naledi remains and that thereby help contextualize the spatial concentration of H. naledi skeletal material. The most substantial change in the data presented is a much expanded reanalysis of the Puzzle Box area. This reanalysis provides greater clarity on how previously published descriptions relate to the new evidence. The reanalysis also provides the data to integrate the detailed information on bone identification fragmentation, and spatial taphonomy from this area with the new excavation results from the other areas.

      In addition to Results, the reorganization also affected the manuscript’s Introduction section. Where the previous version led directly from a brief review of Pleistocene burial into the description of the results, this revised manuscript now includes a review of previous studies of the Rising Star cave system. This review directly addresses referee comments that express some hesitation to accept previous results concerning the structure and formation of sediments, the accessibility of the Dinaledi Subsystem, the geochronological setting of the H. naledi remains, and the relation of the Dinaledi Subsystem to nearby cave areas. Some parts of this overview are further expanded in the Supplementary Information to enable readers to dive more deeply into the previous literature on the site formation and geological configuration of the Rising Star cave system without needing to digest the entirety of the cited sources.

      The Discussion section of the revised manuscript is differentiated from Results and focuses on several areas where the evidence presented in this study may benefit from greater context. One new section addresses hypothesis testing and parsimony for Pleistocene burial evidence, which we address at greater length in this summary below. The majority of the Discussion concerns the criteria for recognizing evidence for burial as applied in other studies. In this research we employ a minimal definition but other researchers have applied varied criteria. We consider whether these other criteria have relevance in light of our observations and whether they are essential to the recognition of burial evidence more broadly.

      Vocabulary

      We introduce the term “cultural burial” in this revised manuscript to refer to the burial of dead bodies as a mortuary practice. “Burial” as an unmodified term may refer to the passive covering of remains by sedimentary processes. Use of the term “intentional burial” would raise the question of interpreting intent, which we do not presume based on the evidence presented in this research. The relevant question in this case is whether the process of burial reflects repeated behavior by a group. As we received input from various colleagues it became clear that burial itself is a highly loaded term. In particular there is a common assumption within the literature and among professionals that burial must by definition be symbolic. We do not take any position on that question in this manuscript, and it is our hope that the term “cultural burial” may focus the conversation around the extent that the behavioral evidence is repeated and patterned.

      Sedimentology and geochemistry of Dinaledi Feature 1

      Reviewer 4 provided detailed comments on the sedimentological and geochemical context that we report in the manuscript. One outside review (Foecke et al. 2024) included some of the points raised by reviewer 4, and additionally addressed the reporting of geochemical and sedimentological data in previous work that we cite.

      To address these comments we have revised the sedimentary context and micromorphology of sediments associated with Dinaledi Feature 1. In the new text we demonstrate the lack of microstratigraphy (supported by grain size analysis) in the unlithified mud clast breccia (UMCB), while such a microstratigraphy is observed in the laminated orange-red mudstones (LORM) that contribute clasts to the UMCB. Thus, we emphasize the presence and importance of a laterally continuous layer of LORM nature occurring at a level that appears to be the maximum depth of fossil occurrence. This layer is severely broken under extensive accumulation of fossils such as Feature 1 and only evidenced by abundant LORM clasts within and around the fossils.

      We have completely reworked the geochemical context associated with Feature 1 following the comments of reviewer 4. We described the variations and trends observed in the major oxides separate from trace and rare-earth elements. We used Harker variations plots to assess relationships between these element groups with CaO and Zn, followed by principal component analysis of all elements analyzed. The new geochemical analysis clearly shows that Feature 1 is associated with localized trace element signatures that exist in the sediments only in association with the fossil bones, which suggests lack of postdepositional mobilization of the fossils and sediments. We additionally have included a fuller description of XRF methods.

      To clarify the relation of all results to the features described in this study, we removed the geochemical and sedimentological samples from other sites within the Dinaledi Subsystem. These localities within the fissure network represent only surface collection of sediment, as no excavation results are available from those sites to allow for comparison in the context of assessing evidence of burial. These were initially included for comparison, but have now been removed to avoid confusion.

      Micromorphology of sediments

      Some referees (1, 3, and 4) and other commentators (including Martinón-Torres et al. 2024) have suggested that the previous manuscript was deficient due to an insufficient inclusion of micromorphological analysis of sediments. Because these commentators have emphasized this kind of evidence as particularly important, we review here what we have included and how our revision has addressed this comment. Previous work in the Dinaledi Chamber (Dirks et al., 2015; 2017) included thin section illustrations and analysis of sediment facies, including sediments in direct association with H. naledi remains within the Puzzle Box area. The previous work by Wiersma and coworkers (2020) used micromorphological analysis as one of several approaches to test the formation history of Unit 3 sediments in the Dinaledi Subsystem, leading to the interpretation of autobrecciation of earlier Unit 1 sediment. In the previous version of this manuscript we provided citations to this earlier work. The previous manuscript also provided new thin section illustrations of Unit 3 sediment near Dinaledi Feature 1 to place the disrupted layer of orange sediment (now designated the laminated orange silty mudstone unit) into context.

      In the new revised manuscript we have added to this information in three ways. First, as noted above in response to reviewer 4, we have revised and added to our discussion of micromorphology within and adjacent to the Dinaledi Feature 1. Second, we have included more discussion in the Supplementary Information of previous descriptions of sediment facies and associated thin section analysis, with illustrations from prior work (CC-BY licensed) brought into this paper as supplementary figures, so that readers can examine these without following the citations. Third, we have included Figure 10 in the manuscript which includes six panels with microtomographic sections from the Hill Antechamber Feature. This figure illustrates the consistency of sub-unit 3b sediment in direct contact with H. naledi skeletal material, including anatomically associated skeletal elements, with previous analyses that demonstrate the angular outlines and chaotic orientations of LORM clasts. It also shows density contrasts of sediment in immediate contact with some skeletal elements, the loose texture of this sediment with air-filled voids, and apparent invertebrate burrowing activity. To our knowledge this is the first application of microtomography to sediment structure in association with a Pleistocene burial feature.

      To forestall possible comments that the revised manuscript does not sufficiently employ micromorphological observations, or that any one particular approach to micromorphology is the standard, we present here some context from related studies of evidence from other research groups working at varied sites in Africa, Europe, and Asia. Hodgkins et al. (2021) noted: “Only a handful of micromorphological studies have been conducted on human burials and even fewer have been conducted on suspected burials from Paleolithic or hunter-gatherer contexts.” In that study, one supplementary figure with four photomicrographs of thin sections of sediments was presented. Interpretation of the evidence for a burial pit by Hodgkins et al. (2021) noted the more open microstructure of sediment but otherwise did not rely upon the thin section data in characterizing the sediments associated with grave fill. Martinón-Torres et al. (2021) included one Extended Data figure illustrating thin sections of sediments and bone, with two panels showing sediments (the remainder showing bone histology). The micromorphological analysis presented in the supplementary information of that paper was restricted to description of two microfacies associated with the proposed “pit” in that study. That study did carry out microCT scanning of the partially-prepared skeletal remains but did not report any sediment analysis from the microtomographic results. Maloney et al. (2022) reported no micromorphological or thin section analysis. Pomeroy et al. (2020a) included one illustration of a thin section; this study may be regarded as a preliminary account rather than a full description of the work undertaken. Goldberg et al. (2017) analyzed the geoarchaeology of the Roc de Marsal deposits in which possible burial-associated sediments had been fully excavated in the 1960s, providing new morphological assessments of sediment facies; the supplementary information to this work included five scans (not microscans) of sediment thin sections and no microphotographs. Fewlass et al. (2023) presented no thin section or micromorphological illustrations or methods. In summary of this research, we note that in one case micromorphological study provided observations that contributed to the evidence for a pit, in other cases micromorphological data did not test this hypothesis, and many researchers do not apply micromorphological techniques in their particular contexts.

      Sediment micromorphology is a growing area of research and may have much to provide to the understanding of ancient burial evidence as its standards continue to develop (Pomeroy et al. 2020b). In particular microtomographic analysis of sediments, as we have initiated in this study, may open new horizons that are not possible with more destructive thin-section preparation. In this manuscript, the thin section data reveals valuable evidence about the disruption of sediment structure by features within the Dinaledi Chamber, and microtomographic analysis further documents that the Hill Antechamber Feature reflects similar processes, in addition to possible post-burial diagenesis and invertebrate activity. Following up in detail on these processes will require further analysis outside the scope of this manuscript.

      Access into the Dinaledi Subsystem

      Reviewer 1 emphasizes the difficulty of access into the Dinaledi Subsystem as a reason why the burial hypothesis is not parsimonious. Similar comments have been made by several outside commentators who question whether past accessibility into the Dinaledi Subsystem may at one time have been substantially different from the situation documented in previous work. Several pieces of evidence are relevant to these questions and we have included some discussion of them in the Introduction, and additionally include a section in the Supplementary Information (“Entrances to the cave system”) to provide additional context for these questions. Homo naledi remains are found not only within the Dinaledi Subsystem but also in other parts of the cave system including the Lesedi Chamber, which is similarly difficult for non-expert cavers to access. The body plan, mass, and specific morphology of H. naledi suggest that this species would be vastly more suited to moving and climbing within narrow underground passages than living people. On this basis it is not unparsimonious to suggest that the evidence resulted from H. naledi activity within these spaces. We note that the accessibility of the subsystem is not strictly relevant to the hypothesis of cultural burial, although the location of the remains does inform the overall context which may reflect a selection of a location perceived as special in some way.

      Stuffing bodies down the entry to the subsystem

      Reviewer 3 suggests that one explanation for the emplacement of articulated remains at the top of the sloping floor of the Hill Antechamber is that bodies were “stuffed” into the chute that comprises the entry point of the subsystem and passively buried by additional accumulation of remains. This was one hypothesis presented in earlier work (Dirks et al. 2015) and considered there as a minimal explanation because it did not entail the entry of H. naledi individuals into the subsystem. The further exploration (Elliott et al. 2021) and ongoing survey work, as well as this manuscript, all have resulted in data that rejects this hypothesis. The revised manuscript includes a section in the results “Deposition upon a talus with passive burial” that examines this hypothesis in light of the data.

    1. On 2024-08-10 09:02:59, user ani1977 wrote:

      Trying to download the data but facing "unzip: inflate error"? Below is the details of file and wget-log

      wget --ftp-user=MSV000095162 --ftp-password=hgu89\!hgb -r " ftp://MSV000095162@massive.ucsd.edu/raw/Raw2/8 --2024-08-09 15:24:16-- ftp://MSV000095162@massive.ucsd.edu/raw/Raw2/80-2.zip => ‘ http://massive.ucsd.edu/raw/Raw2/.listing’ Resolving http://massive.ucsd.edu ... 132.249.211.16 Connecting to http://massive.ucsd.edu |132.249.211.16|:21... connected. Logging in as MSV000095162 ... Logged in! ==> SYST ... done. ==> PWD ... done. ==> TYPE I ... done. ==> CWD (1) /raw/Raw2 ... done. ==> PASV ... done. ==> LIST ... done. http://massive.ucsd.edu/raw/Raw2/.listing `[ <=>

      2024-08-09 15:24:18 (3.41 MB/s) - ‘[http://massive.ucsd.edu/raw/Raw2/.listing’](http://massive.ucsd.edu/raw/Raw2/.listing’:GdLBFgBQICCkUtR2vybJkVCxpvA "http://massive.ucsd.edu/raw/Raw2/.listing’")saved [326]

      Removed ‘[http://massive.ucsd.edu/raw/Raw2/.listing’](http://massive.ucsd.edu/raw/Raw2/.listing’:GdLBFgBQICCkUtR2vybJkVCxpvA "http://massive.ucsd.edu/raw/Raw2/.listing’").<br /> --2024-08-09 15:24:18--[ftp://MSV000095162@massive.ucsd.edu/raw/Raw2/80-2.zip](ftp://MSV000095162@massive.ucsd.edu/raw/Raw2/80-2.zip:Y9SeoGJEmkqa4Bl1VBmCaBwr2E4 "ftp://MSV000095162@massive.ucsd.edu/raw/Raw2/80-2.zip")=> ‘[http://massive.ucsd.edu/raw/Raw2/80-2.zip’](http://massive.ucsd.edu/raw/Raw2/80-2.zip’:O5ImlsutePNNEZ3_IaMs7NMHh6I "http://massive.ucsd.edu/raw/Raw2/80-2.zip’")==> CWD not required.<br /> ==> PASV ... done. ==> RETR[http://80-2.zip](http://80-2.zip:qseKzsDg34VmYze3eF17fnrtQFY "http://80-2.zip")... done.<br /> Length: 173832044544 (162G)[http://massive.ucsd.edu/raw/Raw2/80-2.zip](http://massive.ucsd.edu/raw/Raw2/80-2.zip:eXxGRDBvC8fBOclrVFpPmsM5rs0 "http://massive.ucsd.edu/raw/Raw2/80-2.zip")63%[====================================================================================>

      2024-08-10 03:00:11 (2.54 MB/s) - Data transfer aborted.<br /> Retrying.

      --2024-08-10 03:00:12--[ftp://MSV000095162@massive.ucsd.edu/raw/Raw2/80-2.zip](ftp://MSV000095162@massive.ucsd.edu/raw/Raw2/80-2.zip:Y9SeoGJEmkqa4Bl1VBmCaBwr2E4 "ftp://MSV000095162@massive.ucsd.edu/raw/Raw2/80-2.zip")(try: 2) => ‘[http://massive.ucsd.edu/raw/Raw2/80-2.zip’](http://massive.ucsd.edu/raw/Raw2/80-2.zip’:O5ImlsutePNNEZ3_IaMs7NMHh6I "http://massive.ucsd.edu/raw/Raw2/80-2.zip’")==> CWD not required.<br /> ==> PASV ...<br /> Error in server response, closing control connection.<br /> Retrying.

      --2024-08-10 03:00:14--[ftp://MSV000095162@massive.ucsd.edu/raw/Raw2/80-2.zip](ftp://MSV000095162@massive.ucsd.edu/raw/Raw2/80-2.zip:Y9SeoGJEmkqa4Bl1VBmCaBwr2E4 "ftp://MSV000095162@massive.ucsd.edu/raw/Raw2/80-2.zip")(try: 3) => ‘[http://massive.ucsd.edu/raw/Raw2/80-2.zip’](http://massive.ucsd.edu/raw/Raw2/80-2.zip’:O5ImlsutePNNEZ3_IaMs7NMHh6I "http://massive.ucsd.edu/raw/Raw2/80-2.zip’")Connecting to[http://massive.ucsd.edu](http://massive.ucsd.edu:fnwcjh3df7hMG7pe2bJ2IK5qRBc "http://massive.ucsd.edu")|132.249.211.16|:21... connected.<br /> Logging in as MSV000095162 ... Logged in!<br /> ==> SYST ... done. ==> PWD ... done.<br /> ==> TYPE I ... done. ==> CWD (1) /raw/Raw2 ... done.<br /> ==> PASV ... done. ==> REST 111049517712 ... done.<br /> ==> RETR[http://80-2.zip](http://80-2.zip:qseKzsDg34VmYze3eF17fnrtQFY "http://80-2.zip")... done.<br /> Length: 173832044544 (162G), 62782526832 (58G) remaining[http://massive.ucsd.edu/raw/Raw2/80-2.zip](http://massive.ucsd.edu/raw/Raw2/80-2.zip:eXxGRDBvC8fBOclrVFpPmsM5rs0 "http://massive.ucsd.edu/raw/Raw2/80-2.zip")66%[+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++===>

      2024-08-10 04:49:45 (698 KB/s) - Data connection: Connection timed out; Data transfer aborted.<br /> Retrying.

      --2024-08-10 04:49:48--[ftp://MSV000095162@massive.ucsd.edu/raw/Raw2/80-2.zip](ftp://MSV000095162@massive.ucsd.edu/raw/Raw2/80-2.zip:Y9SeoGJEmkqa4Bl1VBmCaBwr2E4 "ftp://MSV000095162@massive.ucsd.edu/raw/Raw2/80-2.zip")(try: 4) => ‘[http://massive.ucsd.edu/raw/Raw2/80-2.zip’](http://massive.ucsd.edu/raw/Raw2/80-2.zip’:O5ImlsutePNNEZ3_IaMs7NMHh6I "http://massive.ucsd.edu/raw/Raw2/80-2.zip’")==> CWD not required.<br /> ==> PASV ...<br /> Error in server response, closing control connection.<br /> Retrying.

      --2024-08-10 04:49:52--[ftp://MSV000095162@massive.ucsd.edu/raw/Raw2/80-2.zip](ftp://MSV000095162@massive.ucsd.edu/raw/Raw2/80-2.zip:Y9SeoGJEmkqa4Bl1VBmCaBwr2E4 "ftp://MSV000095162@massive.ucsd.edu/raw/Raw2/80-2.zip")(try: 5) => ‘[http://massive.ucsd.edu/raw/Raw2/80-2.zip’](http://massive.ucsd.edu/raw/Raw2/80-2.zip’:O5ImlsutePNNEZ3_IaMs7NMHh6I "http://massive.ucsd.edu/raw/Raw2/80-2.zip’")Connecting to[http://massive.ucsd.edu](http://massive.ucsd.edu:fnwcjh3df7hMG7pe2bJ2IK5qRBc "http://massive.ucsd.edu")|132.249.211.16|:21... connected.<br /> Logging in as MSV000095162 ... Logged in!<br /> ==> SYST ... done. ==> PWD ... done.<br /> ==> TYPE I ... done. ==> CWD (1) /raw/Raw2 ... done.<br /> ==> PASV ... done. ==> REST 115747059452 ... done.<br /> ==> RETR[http://80-2.zip](http://80-2.zip:qseKzsDg34VmYze3eF17fnrtQFY "http://80-2.zip")... done.<br /> Length: 173832044544 (162G), 58084985092 (54G) remaining[http://massive.ucsd.edu/raw/Raw2/80-2.zip](http://massive.ucsd.edu/raw/Raw2/80-2.zip:eXxGRDBvC8fBOclrVFpPmsM5rs0 "http://massive.ucsd.edu/raw/Raw2/80-2.zip")100%[+++++++++++++++++++++++++++++++++++++++++++++++++========================>] 161.89G 6.57MB/s in 4h 43m =

      2024-08-10 09:33:02 (3.26 MB/s) - ‘[http://massive.ucsd.edu/raw/Raw2/80-2.zip’](http://massive.ucsd.edu/raw/Raw2/80-2.zip’:O5ImlsutePNNEZ3_IaMs7NMHh6I "http://massive.ucsd.edu/raw/Raw2/80-2.zip’")saved [173832044544]

      FINISHED --2024-08-10 09:33:02--<br /> Total wall clock time: 18h 8m 45s<br /> Downloaded: 1 files, 54G in 18h 8m 25s (869 KB/s)<br /> ✓

       10/08/2024   09:45.41   /mnt/z/d80  cd[http://massive.ucsd.edu/raw/Raw2/](http://massive.ucsd.edu/raw/Raw2/:PrfWHyVCVLrSn8uZCBTcW9kil6E "http://massive.ucsd.edu/raw/Raw2/")

       10/08/2024   09:45.56   /mnt/z/d80/[http://massive.ucsd.edu/raw/Raw2](http://massive.ucsd.edu/raw/Raw2:wwDJqyDJge16is8KfZcnO2bFthc "http://massive.ucsd.edu/raw/Raw2") unzip[http://80-2.zip](http://80-2.zip:qseKzsDg34VmYze3eF17fnrtQFY "http://80-2.zip")Archive:[http://80-2.zip](http://80-2.zip:qseKzsDg34VmYze3eF17fnrtQFY "http://80-2.zip")creating: BORA_10_16_HEK_80C_C10_RC11_1_22017.d/<br /> creating: BORA_10_16_HEK_80C_C10_RC11_1_22017.d/22017.m/<br /> creating: BORA_10_16_HEK_80C_C10_RC11_1_22017.d/22017.m/backup-2023-12-13.m/<br /> inflating: BORA_10_16_HEK_80C_C10_RC11_1_22017.d/22017.m/backup-2023-12-13.m/diaSettings.diasqlite<br /> inflating: BORA_10_16_HEK_80C_C10_RC11_1_22017.d/22017.m/backup-2023-12-13.m/hystar.method<br /> inflating: BORA_10_16_HEK_80C_C10_RC11_1_22017.d/22017.m/backup-2023-12-13.m/lock.file<br /> inflating: BORA_10_16_HEK_80C_C10_RC11_1_22017.d/22017.m/backup-2023-12-13.m/Maldi.method<br /> inflating: BORA_10_16_HEK_80C_C10_RC11_1_22017.d/22017.m/backup-2023-12-13.m/microTOFQImpacTemAcquisition.method<br /> inflating: BORA_10_16_HEK_80C_C10_RC11_1_22017.d/22017.m/backup-2023-12-13.m/prmSettings.prmsqlite<br /> inflating: BORA_10_16_HEK_80C_C10_RC11_1_22017.d/22017.m/backup-2023-12-13.m/submethods.xml<br /> inflating: BORA_10_16_HEK_80C_C10_RC11_1_22017.d/22017.m/desktop.ini<br /> inflating: BORA_10_16_HEK_80C_C10_RC11_1_22017.d/22017.m/diaSettings.diasqlite<br /> inflating: BORA_10_16_HEK_80C_C10_RC11_1_22017.d/22017.m/hystar.method<br /> inflating: BORA_10_16_HEK_80C_C10_RC11_1_22017.d/22017.m/InstrumentSetup.isset<br /> inflating: BORA_10_16_HEK_80C_C10_RC11_1_22017.d/22017.m/lock.file<br /> inflating: BORA_10_16_HEK_80C_C10_RC11_1_22017.d/22017.m/Maldi.method<br /> inflating: BORA_10_16_HEK_80C_C10_RC11_1_22017.d/22017.m/microTOFQImpacTemAcquisition.method<br /> inflating: BORA_10_16_HEK_80C_C10_RC11_1_22017.d/22017.m/prmSettings.prmsqlite<br /> inflating: BORA_10_16_HEK_80C_C10_RC11_1_22017.d/22017.m/submethods.xml<br /> inflating: BORA_10_16_HEK_80C_C10_RC11_1_22017.d/analysis.tdf<br /> inflating: BORA_10_16_HEK_80C_C10_RC11_1_22017.d/analysis.tdf_bin<br /> inflating: BORA_10_16_HEK_80C_C10_RC11_1_22017.d/chromatography-data-pre.sqlite<br /> inflating: BORA_10_16_HEK_80C_C10_RC11_1_22017.d/chromatography-data.sqlite<br /> inflating: BORA_10_16_HEK_80C_C10_RC11_1_22017.d/chromatography-data.sqlite-journal<br /> inflating: BORA_10_16_HEK_80C_C10_RC11_1_22017.d/SampleInfo.xml<br /> ...<br /> creating: BORA_10_16_HEK_80C_E9_RE10_1_22076.d/<br /> creating: BORA_10_16_HEK_80C_E9_RE10_1_22076.d/22076.m/<br /> creating: BORA_10_16_HEK_80C_E9_RE10_1_22076.d/22076.m/backup-2023-12-14.m/<br /> inflating: BORA_10_16_HEK_80C_E9_RE10_1_22076.d/22076.m/backup-2023-12-14.m/diaSettings.diasqlite<br /> inflating: BORA_10_16_HEK_80C_E9_RE10_1_22076.d/22076.m/backup-2023-12-14.m/hystar.method<br /> inflating: BORA_10_16_HEK_80C_E9_RE10_1_22076.d/22076.m/backup-2023-12-14.m/lock.file<br /> inflating: BORA_10_16_HEK_80C_E9_RE10_1_22076.d/22076.m/backup-2023-12-14.m/Maldi.method<br /> inflating: BORA_10_16_HEK_80C_E9_RE10_1_22076.d/22076.m/backup-2023-12-14.m/microTOFQImpacTemAcquisition.method<br /> inflating: BORA_10_16_HEK_80C_E9_RE10_1_22076.d/22076.m/backup-2023-12-14.m/prmSettings.prmsqlite<br /> inflating: BORA_10_16_HEK_80C_E9_RE10_1_22076.d/22076.m/backup-2023-12-14.m/submethods.xml<br /> inflating: BORA_10_16_HEK_80C_E9_RE10_1_22076.d/22076.m/desktop.ini<br /> inflating: BORA_10_16_HEK_80C_E9_RE10_1_22076.d/22076.m/diaSettings.diasqlite<br /> inflating: BORA_10_16_HEK_80C_E9_RE10_1_22076.d/22076.m/hystar.method<br /> inflating: BORA_10_16_HEK_80C_E9_RE10_1_22076.d/22076.m/InstrumentSetup.isset<br /> inflating: BORA_10_16_HEK_80C_E9_RE10_1_22076.d/22076.m/lock.file<br /> inflating: BORA_10_16_HEK_80C_E9_RE10_1_22076.d/22076.m/Maldi.method<br /> inflating: BORA_10_16_HEK_80C_E9_RE10_1_22076.d/22076.m/microTOFQImpacTemAcquisition.method<br /> inflating: BORA_10_16_HEK_80C_E9_RE10_1_22076.d/22076.m/prmSettings.prmsqlite<br /> inflating: BORA_10_16_HEK_80C_E9_RE10_1_22076.d/22076.m/submethods.xml<br /> inflating: BORA_10_16_HEK_80C_E9_RE10_1_22076.d/analysis.tdf<br /> inflating: BORA_10_16_HEK_80C_E9_RE10_1_22076.d/analysis.tdf_bin<br /> unzip: inflate error`

    1. On 2024-08-08 08:28:33, user Ashleigh Shannon wrote:

      AUTHORS COMMENT / UPDATE <br /> 21 June 2024

      Part of this work, describing the inhibition of both the NiRAN and RdRp domains by the guanosine analogue 5’-triphosphate AT-9010 has now been published in Nature Communications ( http://doi.org/10.1038/s41467-022-28113-1)11:bWV2b7IKjzmu98TO74CkaGwQqk8 "http://doi.org/10.1038/s41467-022-28113-1)1") . This includes the structural characterization of the mechanism of inhibition at both active-sites, which importantly revealed that the SARS-CoV-2 NiRAN domain contains a guanine-specific pocket. This has now been confirmed in other studies 2,3, and represents a promising avenue for future drug development studies.

      Following several contradictory reports on the ability of the NiRAN to NMPylate different cofactor proteins (4–6), we have now carried out additional analysis on the protein-priming activity and NMPylation specificity. This has revealed that our nsp8 product is labeled at the primary amine of a non-native glycine residue, present at the N-terminus of the expression construct (following cleavage with TEV protease), and negating the biological relevance of these findings. Surprisingly, removal of this single residue, exposing the true N-terminal alanine, eliminates all labeling and protein-primed activity, stressing the importance for using proteins with native N- and C- termini. Of note, the current publication does not reflect these new findings, but instead has been left as the initially submitted manuscript.

      Later studies have now revealed nsp9 to be the primary target for both NMPylation5,7, and RNAylation8 – a conclusion that we fully agree with.

      With the exception of our finding that nsp8 can be used to prime-synthesis, albeit by a non-native residue, there was no other evidence for protein-priming until recently. In October, 2023, Schmidt et al., published the finding that nsp9 is covalently linked to the 5' end of positive- and negative-sense RNA produced during SARS-CoV-2 infection9. This linkage was found to be regulated by its interaction with the host protein, staphylococcal nuclease domain-containing protein 1 (SND1), which was found to specifically recognize the 5’ end of negative-sense RNA and be important for viral RNA synthesis. Although it is highly frustrating to have spent several years focused on the wrong protein, this finding supports the notion that protein-priming is occurring in CoVs, and opens up a plethora of options for further mechanistic and structural studies.

      Of note, Schmidt et al. also showed that nsp9-linkage on the (-)sense strand mapped roughly between the genome and poly(A) tail. Intriguingly, we also found a similar specificity for this region. Here we show that the SARS-CoV mRTC (nsp12-nsp7-nsp8) can initiate synthesis through a de novo, NiRAN independent pathway, through the synthesis of a pppGpU primer. This dinucleotide is complementary to the last two nucleotides of the SARS-CoV-2 genome, located precisely at the genome-poly(A) junction. De novo initiation was found to directly compete with the (artificial) nsp8-poly(U) protein-primed synthesis. This shows that the poly(A) tail and 3’ genomic RNA sequence elements guide the positioning of the mRTC to the true 3’-end of the genome.

      The specific details behind initiation of RNA synthesis, including the role of SND1 and potential coordination between protein-priming and/or de novo induction remains to be studied. Synthesis initiation, and the precise role of the NiRAN therefore appears to be a complicated story, which remains to be fully elucidated.

      Ashleigh Shannon and Bruno Canard

    1. On 2024-08-05 20:26:55, user Zach Hensel wrote:

      Summary: Zhan et al argue that nearly identical (>99.9% in the S gene or across a whole genome) sequences sampled within a single species indicate transmission within a single species subsequent to a period of host adaptation, assuming the necessity of adaptive mutations in any new host. They repeatedly claim that similar identical genomes were never sampled in different species in SARS outbreak between 2002 and 2004. This is incorrect. In fact, the data cited in this paper includes at least four examples of 100% sequence identity in samples from different species. Thus, the statement "In the SARS-CoV outbreaks, >99.9% genome or S identity was only observed among isolates collected within a narrow window of time from within the same species" was incorrect.

      Example 1: Figure 5 and the manuscript discuss samples SZ3 and SZ16 collected from palm civets in 2004 in Dongmen market in Shenzhen. However, sample SZ13, reported from the same study (Guan et al, Science 2003, accession AY304487) is not mentioned. It was collected from a raccoon dog and it is 100% identical to SZ16 over 8581 bases in a partial genome sequence inclusive of S.

      Furthermore, Guan et al also reported neutralization of SZ16 infection by sera from masked palm civets (3 of 6 samples), raccoon dog (1 of 1), Chinese ferret badger (1 of 2), and human (12/55; 8/20 for humans in wild animal trade). This indicates widespread infection of diverse species by SARS with >99.9% identical genomes without requiring further host adaptation. It's possible this error is inherited from Song et al (PNAS 2005), which twice states that SZ13 is a sample from a civet.

      This error is perplexing since the primary source for this data is clear on what animals were sampled and sequence identity. What's less well known is that Dongmen market continued to be sampled after May (Yaqing et al, Disease Surveillance 2004). Civets, raccoon dogs, and badgers were still available for sampling and still positive by PCR for SARS in October, November, and December 2003. It wasn't until spillover from civets to humans was definitively shown for the Guangzhou outbreak (Wang et al, Emerg Infect Dis 2005) that Guangdong province acted to close wildlife markets and end the civet trade. Dr. Garry notes similar findings at Xinyuan market in Guangzhou in January 2024.

      Example 2 and 3: Zhan et al write, "The human and civet isolates of the 2003/2004 outbreak, which were collected most closely in time and at the site of cross-species transmission, shared only up to 99.79% S identity.” This is incorrect. Wang et al (Emergency Infect Dis 2005) write, "The 286-bp S gene sequences from isolates from the waitress and the physician were identical to 4 of 5 S gene sequences from palm civets from the restaurant.”

      Example 4: The 99.79% identical sequence described by Wang et al (Emergency Infect Dis 2005) is instead from a different patient without a known link to the civet or other animal trade — a patient identified about one month before the first patient in the restaurant cluster. However, “An isogenic strain with the same SNV pattern, B012G, was detected in a palm civet from the Xinyuan animal market” (Kan et al J Virol 2005).

      Conclusion: Repeated claims of no examples of 99.9% identical viruses in different species during SARS outbreaks were provably incorrect when this manuscript was preprinted. Furthermore, the manuscript appeals to near-identity of sequences obtained within Huanan market. Subsequently, it has been shown that both lineage A and lineage B were found in environmental sampling of Huanan market. I hope that authors will revisit their conclusions given the reality of the facts they based their arguments on.

    1. On 2024-08-05 16:25:30, user Ra Hel wrote:

      Hi, thank you for the thorough article. I would like to comment on : "We then employed GB classifiers in subsequent studies and utilized them to exclude studies that cannot discriminate the disease phenotype based on microbial profile."<br /> It would be really useful to have information on what percentage of studies for each disease passed this criterion - did you need to screen 10 studies per disease to have one that would be able to discriminate the phenotype based on microbial profile or rather vice versa that there were very few odd studies that didn't pass this selection.

    1. On 2024-08-05 08:51:01, user Bayazit Yunusbayev wrote:

      Hi, this sentence is not clear to me: "the probability that a given edge carries at least one mutation is inversely proportional to its<br /> time-length (assuming mutations occur as a Poisson process along the edges)." I might be missing here something, but isn't it so that older edges are expected to encounter more mutations?

    1. On 2024-08-04 17:58:17, user Alex wrote:

      The preprint has now been published in PNAS: Alexander M. Sandercock, Jared W. Westbrook, Qian Zhang, Jason A. Holliday. A genome-guided strategy for climate resilience in American chestnut restoration populations. Proceedings of the National Academy of Sciences, 2024; 121 (30) DOI: 10.1073/pnas.2403505121

    1. On 2024-08-03 15:08:52, user Mario Stanke wrote:

      That should be Supplementary Table S1, rather than Table 1. Thanks for pointing the wrong reference out. The training species are also shown as a tree in Supplementary Figure S1.

    2. On 2024-07-31 13:32:54, user kbseah wrote:

      This looks very promising, congratulations! One point was unclear to me: which genome datasets were used to train the current mammalian model? The text refers to Table 1, but that just shows benchmarking results.

    1. On 2024-08-02 17:16:23, user Lonki wrote:

      The authors have errored by assuming that pregabalin and gabapentin are GABA-receptor ligands. Despite their names, they are not.

      As stated in Wikipedia; “Pregabalin inhibits certain calcium channels, namely, it blocks α2δ subunit-containing voltage-dependent calcium channels (VDCCs).[13][26].

      Also from Wikipedia “Despite the fact that gabapentin is a structural GABA analogue, and in spite of its name, it does not bind to the GABA receptors, does not convert into GABATooltip γ-aminobutyric acid or another GABA receptor agonist in vivo, and does not modulate GABA transport or metabolism within the range of clinical dosing.[82]<br /> “Gabapentin is a ligand of the α2δ calcium channel subunit.[82][83] α2δ is an auxiliary protein connected to the main α1 subunit (the channel-forming protein) of high voltage activated voltage-dependent calcium channels (L-type, N-type, P/Q type, and R-type).[13] Gabapentin is not a direct channel blocker”

      Please read up on these drugs and re-interpret your observations.

    1. On 2024-08-02 16:28:00, user Ana Vasque wrote:

      Dear authors,

      First I'd like to congratulate you on your work. I have read your article on the unique reproductive structure of Euphorbia species with great interest. However, I have specific inquiries regarding the filiform structures that were analyzed to determine the floral identity of the cyathium, and I would appreciate further clarification

      In your study, did you observe results that suggested the upregulation of the B and E genes in these filiform structures, which would indicate a implying reduction in flowering? What is the correlation between these findings and previous research, such as that of Prenner and Rudall (2007), which discusses the presence of these structures associated with individual male flowers in certain species and their absence in others? Additionally, Prenner and Rudall (2007) reference Warming (1870), who observed and interpreted these structures as trichomes, noting their formation subsequent to the initiation of staminate flowers

      Could you provide deeper insights into your interpretation of the formation and functionality of these threadlike structures in your analyses?

      I am grateful for your attention to my inquiries. I look forward to your explanations and further discussion on the outcomes of your study.

      Sincerely,

    1. On 2024-08-01 16:35:16, user John McCusker wrote:

      I just saw ( https://www.science.org/content/article/bad-agar-killing-lab-yeast-around-world-where-it-coming ). Many years ago, I had a similar problem with S. cerevisiae and C. albicans (but not E. coli) growth on agar. (Multiple vendors/suppliers, none of whom found anything wrong.) I eventually found that exposure of agar plates to light while drying (or incubating) caused the problem. Dried/incubated plates in dark and no problem.

    1. On 2024-08-01 12:05:39, user Kenichiro Abe wrote:

      Given the result of transcriptional inhibitor treatment with DRB or THZ1 in Fig.S9 and the report of Abe et al.,2010, NSN-SN transition seems like to be occurred independently of RPB1 degradation <br /> If you have any comments about this point, I would be happy to hear that

    1. On 2024-08-01 07:43:29, user Ben wrote:

      Hi. In the graphical abstract, the mus musculus heart cartoon shows a sagittal section. I would believe you wanted to show a heart with a scar.

    1. On 2024-07-27 13:29:27, user Prof. T. K. Wood wrote:

      Schumacher, M.A., et al., 2009. Molecular Mechanisms of HipA-Mediated Multidrug Tolerance and 531 Its Neutralization by HipB. Science, 323 (5912), 396-401, was found to be false in that HipA is a kinase but has nothing to do with EF-Tu. You should actually read the literature about the TA system you cite.

    2. On 2024-07-26 15:46:59, user Prof. T. K. Wood wrote:

      Line 70: Dy et al. is the 3rd report of phage inhibition by TAs; please cite the seminal ref doi: 10.1128/jb.178.7.2044-2050.1996 that predates it by 18 years.

      Line 69: the TA refs related to persistence in this manuscript have been retracted so the link with persistence is not justified by this ref.

    1. On 2024-07-26 14:37:06, user Julien Racle wrote:

      This is a preprint of the following book chapter: Julien Racle and David Gfeller, How to Predict Binding Specificity and Ligands for New MHC-II Alleles with MixMHC2pred, published in HLA Typing, edited by Sebastian Boegel, 2024, Humana Press reproduced with permission of Springer Science+Business Media, LLC, part of Springer Nature. The final authenticated version is available online at: http://dx.doi.org/10.1007/978-1-0716-3874-3_14 .

    1. On 2024-07-25 15:52:06, user Lynda Delph wrote:

      This is terrific work. The authors might want to cite Lynda F. Delph, Keely E. Brown, Luis Diego Ríos, and John K. Kelly. Sex-specific natural selection on SNPs in Silene latifolia. Evolution Letters 6:308-318.<br /> In that paper they say "Finally, alleles favored by paternity selection tended to reduce female survival (Fig. 3D). The negative correlation here (ρ = −0.10) must be driven by a highly polygenic response."

    1. On 2024-07-23 16:05:44, user Zach Hensel wrote:

      I also have one brief comment regarding the discussion of "no expected interference to biological function" for a recombination breakpoint at position 21,314. The authors suggest that mutation to an isoleucine residue at position nsp16 I219 will "occur without having a significant impact on fitness linked to a disruptive function at a protein level."

      The residue I219 is well resolved in crystal structures (e.g. PDB 6WVN) and makes hydrophobic contacts with other hydrophobic sidechains, suggesting that this is unlikely to be the case. This is testable by mutational fitness analysis -- https://jbloomlab.github.io/SARS2-mut-fitness/nsp16.html -- out of possible non-synonymous mutations that are sufficiently sampled during the pandemic, I219 tolerates substitution to leucine or valine, methionine is slightly disfavored, and mutation to any of arginine, lysine, or threonine is very rare. This suggests a significant fitness impact of substitutions at this position.

    2. On 2024-07-23 16:03:09, user ryhisner wrote:

      This preprint claims that the BA.2.86 lineage, which first appeared in the sequencing record in late July 2023, evolved gradually, while cryptically circulating over the course of about 17 months, beginning in early 2022. I do not find any of the evidence presented here convincing.

      Three clusters of sequences, named C1, C2, and C3, are cited as evidence of cryptic BA.2.86 circulation prior to its emergence on the world stage. Each cluster has a far simpler, more parsimonious explanation than the hypothesis presented by the authors.

      I have documented and analyzed cluster C1 as they have been uploaded, beginning in October 2023, and it appears to be a classic case of a chronic infection in an immunocompromised individual who may have transmitted to one other person. Three of the four C1 sequences have matching metadata—same location, age, and sex—while the fourth has the same location and age but a different sex. Furthermore, there are two additional, closely related sequences collected in 2024 (EPI_ISL_18969735 and EPI_ISL_19259365) whose metadata also match the other three.

      It seems very unlikely the C1 sequences have any relation to BA.2.86. I count approximately 17 spike substitutions and one large deletion (∆138-144) that are in two or more of these four sequences but not in BA.2.86, including P9L, H69D, K77E, T95I, ∆136-144, R158K, Q183E, G213E, D215G, R346T, L452R, F486L, V615A, V642G, H681Y, L841R, D936Y, and D1146N. None of the four C1 sequences have the most distinctive spike mutations of BA.2.86, such as ins16MPLF, the triple-nucleotide F157S-R158G, A264D, I332V, K356T, L452W, or ∆V483. All C1 sequences have ORF1a:L3201F, which is found in BA.2 but not BA.2.86, and none of the C1 sequences possess any of the seven synonymous mutations found in the BA.2.86 branch (C8293T, T13339C, T15756A, A18492G, C21622T, C25207T, C26681T).

      Of the more than 30 spike mutations (relative to baseline BA.2) in BA.2.86, I only see four that are shared between BA.2.86 and the C1 sequences: R403K, A484K, R493Q (reversion), and P621S. All of these are extremely convergent in highly mutated, chronic-infection sequences. I maintain a list of such sequences, and the R493Q reversion is the single most common private mutation, occurring 304 times independently, while R403K appears 103 times, and P621S 60 times. A484K, despite being a two-nucleotide mutation, has independently evolved at least 15 times among the sequences I've recorded. (G446S—very common both in chronic-infection sequences and circulating lineages—is in 2/6 sequences from this cluster).

      Furthermore, the most recently uploaded sequence from this cluster—collected on May 22, 2024, but not listed in the C1 cluster in this paper—contains 16 new spike mutations (at least two of which I suspect can be attributed to sotrovimab treatment). Four of the 16 new spike mutations are also in BA.2.86 (I332V, K356T, L486P, and S939F), a textbook example of how, through convergent intrahost evolution, chronic-infection sequences can come to acquire mutations found in other chronic infections and in unrelated circulating lineages.

      I do not see any resemblance between BA.2.86 and the C1 cluster in the non-spike part of the genome, apart from M:A104V, which is commonly found in chronic-infection sequences (32 independent acquisitions by my count), and is also found in the Pango-designated GS.4/5 lineage (XBB.2.3.11.4/5). It seems to me that this is a case of a chronically infected, immunocompromised individual who developed a few mutations also found in BA.2.86—mutations which are convergent in such chronic infections—and who may have transmitted to one other person (assuming there was not a mistake in the sex assignment of EPI_ISL_18415854, in which case these sequences almost certainly all came from the same patient).

      The phylogenetic relationship between these sequences, as determined by USHER (Ultrafast Sample placement on Existing tRee, created at University of California Santa Cruz and maintained by Angie Hinrichs), can be seen at the following link: <br /> https://nextstrain.org/fetch/raw.githubusercontent.com/ryhisner/posited_BA.2.86_intermediates/main/BA.2.86_C1_posited_intermediates_Mexico_6_seq.json?c=gt-S_841&gmax=25384&gmin=21563&label=id:node_3061533

      The C2 cluster (12 sequences of XBB.1.5.90, 11 from Japan, one from Finland) does not seem to resemble BA.2.86 at all. It is part of a large branch of XBB.1.5.90 (>400 sequences) with S:P621S, also found in BA.2.86, and the only other private mutation I see that it shares with BA.2.86 (but not other hundreds of other XBB.1.5.90 sequences) is C26681T, which is a highly homoplasic synonymous mutation in the coding region of M. Perhaps I am overlooking something, but the C2 cluster looks like a relatively humdrum branch of XBB.1.5.90 to me.

      The Usher tree for these 12 sequences from C2 can be viewed here: https://nextstrain.org/fetch/raw.githubusercontent.com/ryhisner/posited_BA.2.86_intermediates/main/BA.2.86_C2_posited_intermediates_JPN_FIN_12_seq.json?c=userOrOld&label=id:node_4041697

      The C3 cluster consisting of 10 sequences from Sarawak, Malaysia, were all uploaded on the same day (2024-1-25), bear the same collection date (2022-3-11), and have spikes identical to JN.1—including S:L455S—from S:356 to S:681, while the rest of spike is identical to baseline BA.2. The remainder of the genome in these sequences is extremely odd. Two sequences contain the XBB mutation ORF1b:S959P. Seven have the universal BA.2.86 mutation ORF1a:N2526S, while three lack it. One has the BA.2.86.1 mutation ORF1a:K1973R. Seven of the ten have ORF1a:L3201F, which is absent from all BA.2.86. No dropout is indicated in any of the sequences.

      The same Malaysian lab uploaded 321 other sequences (EPI_ISL_18821317-18821647), all from Sarawak, Malaysia, on the same day they uploaded the 10 C3 sequences. The collection dates of these sequences range from 2022-2-27 to 2024-1-9 and include 153 JN.1* sequences and 12 XBB* sequences. As Zach Hensel has noted, six of the ten C3 sequences have G19677T (ORF1b:2070H), which is the defining mutation of BA.2.40, a variant that made up about 60% of all sequences in Sarawak, Malaysia, in mid-March 2022. (Source: https://cov-spectrum.org/explore/Malaysia/AllSamples/from%3D2022-01-15%26to%3D2022-04-28/variants?nextcladeQcSnpClustersScoreTo=55&variantQuery=Nextcladepangolineage%3ABA.2.40*&)

      Sixty-three sequences in this upload are categorized by Nextclade as being BA.2.40 and have 0-3 mutations relative to baseline BA.2.40. Most suspicious of all, 29 of the 153 JN.1* sequences in this same upload also have G19677T. From July 1, 2023 to the present, just 77 sequences categorized by CovSpectrum as BA.2.86* (from 13 different Pango-designated lineages) have had G19677T, with 30 of those coming from Malaysia. (Source: https://cov-spectrum.org/explore/World/AllSamples/from%3D2023-07-01%26to%3D2024-07-14/variants?variantQuery=Nextcladepangolineage%3ABA.2.86*+%26+G19677T&)

      It seems clear that the 10 C2 sequences were BA.2 sequences contaminated by JN.1 sequences from the same upload.

      The authors list 129 sequences they claim shorten the branch leading to BA.2.86, of which I was able to find 128 on GISAID. Ten of these sequences are from the C3 Malaysian cluster described above, along with one additional sequence from the same upload. Apart from these C3 sequences, there are only six sequences with collection dates preceding the first BA.2.86 sequences. All others were collected more than seven weeks after the first BA.2.86 sequences. Six sequences were collected between 7-14 weeks after the first BA.2.86, while the remaining 105 sequences were collected more than 15 weeks afterward.

      It would be surprising if one could not find hundreds of such "hybrid" sequences due purely to contamination. Such sequences have frequently appeared in the sequencing record throughout the pandemic. A few sequences may result from coinfection, but the quality of these sequences, described below, along with the fact that a large proportion of them come from a small number of labs with records of quality-control issues, support the hypothesis that these sequences result from contamination or other lab errors.

      All of the sequences in this list are low quality. They feature a mixture of extensive dropout (particularly in spike), frameshifts, large numbers of mixed nucleotides, clearly artifactual reversions, and mutations from multiple lineages (primarily BA.2.86 and XBB) with no distinct breakpoints. Many of these sequences come from labs known to have frequent quality-control issues. For example, there are 52 sequences from the United States, but none are from the CDC, whose sequences are virtually always first-rate. Instead, they come from smaller local and state labs, whose sequencing quality is often inconsistent.

      The 28 sequences from Texas, for example, come from city hospitals. The average Nextclade qc score of these sequences is 2901 (median 2882). Anything over 100 is designated "bad" by Nextclade. The average number of ambiguous nucleotides per sequences is 17 (median 17), and they average 690 nucleotides of dropout. (EPI_ISL_16599325, EPI_ISL_16599747, EPI_ISL_18546432, EPI_ISL_18690036, EPI_ISL_18690080, EPI_ISL_18690421, EPI_ISL_18690466, EPI_ISL_18690496, EPI_ISL_18743044, EPI_ISL_18743073, EPI_ISL_18743094, EPI_ISL_18743097, EPI_ISL_18743159, EPI_ISL_18743350, EPI_ISL_18743431, EPI_ISL_18743464, EPI_ISL_18743470, EPI_ISL_18743477, EPI_ISL_18743592, EPI_ISL_18816401, EPI_ISL_18816517, EPI_ISL_18816528, EPI_ISL_18816612, EPI_ISL_18816709, EPI_ISL_18816890, EPI_ISL_18816980, EPI_ISL_18874714, EPI_ISL_18908998)

      Similarly, the 17 sequences on this list from Italy are all from the same lab, have an average Nextclade qc score of 2356 (median 2253) and average 1349 nucleotides of dropout. Some sequences on the list are somewhat less bad than these, but none are high-quality. <br /> (EPI_ISL_18496352, EPI_ISL_18674020, EPI_ISL_18677248, EPI_ISL_18721993, EPI_ISL_18722001, EPI_ISL_18722007, EPI_ISL_18722009, EPI_ISL_18755145, EPI_ISL_18792827, EPI_ISL_18792828, EPI_ISL_18792829, EPI_ISL_18792831, EPI_ISL_18820147, EPI_ISL_18820149, EPI_ISL_18820150, EPI_ISL_18820154, EPI_ISL_18820157)

      Finally, I also examined the list of 100 sequences from Supplementary Data 1, Table 5, containing genomes posited by the authors to be recombinants related to ancestors of BA.2.86. These sequences seem to me to fall into five different categories.

      First, there are numerous sequences here that are also listed in one of the C1-C3 clusters—three sequences from C1 (EPI_ISL_18415832, EPI_ISL_18415854, EPI_ISL_18798234), three from C2 (EPI_ISL_18040349, EPI_ISL_18060516, EPI_ISL_18106303, EPI_ISL_18116248), and four from C3 (EPI_ISL_18821484, EPI_ISL_18821485, EPI_ISL_18821487).

      Second, nine of the sequences appear to be fairly unremarkable XDD sequences, which is a designated JN.1/EG.5.1.1 recombinant (EPI_ISL_18617332, EPI_ISL_18706019, EPI_ISL_18706171, EPI_ISL_18553650, EPI_ISL_18653986, EPI_ISL_18531477, EPI_ISL_18569411, EPI_ISL_18695627).

      The third category consists of what seem to me to be relatively normal sequences from a variety of Omicron lineages but with little resemblance to BA.2.86. Some of them have extensive dropout and come from labs known for high rates of artifacts and contamination. (EPI_ISL_18076898, EPI_ISL_17990180, EPI_ISL_18062641, EPI_ISL_18000549, EPI_ISL_18042058, EPI_ISL_18104305, EPI_ISL_18070023, EPI_ISL_18044667, EPI_ISL_18111437, EPI_ISL_15153261, EPI_ISL_17255807, EPI_ISL_15153261, EPI_ISL_16282414, EPI_ISL_16457740)

      The fourth category is BA.2.86 or JN.1 sequences that either don't strike me as very unusual or else have extensive dropout and artifactual reversions. A few of these are from unreliable labs. (EPI_ISL_18097345, EPI_ISL_18556860, EPI_ISL_18567791, EPI_ISL_18682823, EPI_ISL_18705393, EPI_ISL_18635682, EPI_ISL_18503709, EPI_ISL_18400531, EPI_ISL_18717823, EPI_ISL_18446586, EPI_ISL_18584588, EPI_ISL_18631046, EPI_ISL_18700743, EPI_ISL_18675075, EPI_ISL_18659819, EPI_ISL_18713456, EPI_ISL_18636806, EPI_ISL_18704459, EPI_ISL_18686183)

      The fifth and largest category is highly divergent sequences almost certainly deriving from chronic infections, but which appear to me to bear almost no resemblance to BA.2.86 apart from the possession of a few mutations that are widely convergent in such sequences. I've documented most of these, and almost all contain a large number of mutations and deletions not found in BA.2.86 and lack the vast majority of BA.2.86 mutations.

      I hope I haven't misinterpreted any of the authors' hypotheses or data.

      -Ryan Hisner