10,000 Matching Annotations
  1. Last 7 days
    1. AbstractBackground In recent years, biodiversity data management has emerged as a critical pillar in global conservation efforts. Today, the ability to efficiently collect, structure, and analyze biodiversity data is central to breakthroughs in conservation, drug development, disease monitoring, ecological forecasting, and agri-tech innovation. However, due to the vastness and heterogeneity of biodiversity data, it is often confined to databases for specific research areas in isolated formats and disconnected from other relevant resources. Crucial components of such data in kingdom Plantae comprise of metabolomes - the vast array of compounds produced by plants; traits - measurable characteristics of plants that influence their growth, survival, and reproduction, and that affect ecosystem processes; and biotic interactions - relationships of plants with other living organisms, affecting the ecosystem functions.Results In this work, we present METRIN-KG (MEtabolomes, TRaits, and INteractions-Knowledge Graph) a powerful data resource simplifying the integration of diverse and heterogeneous data resources such as plant metabolomes, traits, and biotic interactions.Conclusions The proposed knowledge graph provides an interface to interactively search for data relating plant metabolomes, traits, and interactions. This, in turn, will facilitate development of research questions in life-sciences. In this context, we provide representative case studies on how to frame queries that can be used to search for relevant data in the knowledge graph.

      This work has been peer reviewed in GigaScience (see https://doi.org/10.1093/gigascience/giag051), which carries out open, named peer-review. These reviews are published under a CC-BY 4.0 license and were as follows:

      Reviewer 2:

      The authors present METRIN-KG, a knowledge graph integrating plant metabolomic, trait, and biotic interaction datasets. The work should have substantial value for multiple plant sciences and ecology domains. The overall effort to harmonize disparate resources into an integrated, semantically coherent resource is impressive. The methodology includes a notable pipeline for ontology alignment across multiple sources. However, despite its technical strengths, several issues regarding data accessibility and manuscript structure and clarity should be addressed. The manuscript is highly technical throughout. While this level of precision is exemplary, it may alienate biologically oriented readers, and effort should be made so that the impact and manuscript is clear to a larger audience.

      Major comments

      1. The online deployment currently presents several problems that must be resolved before publication.

      2. The existence of multiple SPARQL UIs (https://kg.earthmetabolome.org/metrin which redirects to https://qlever.earthmetabolome.org/metrin-kg/ and https://sib-swiss.github.io/sparql-editor/metrin-kg) is not explained, and the redundancy is potentially confusing.

      3. At the time of review, https://qlever.earthmetabolome.org/metrin-kg/ showed expired SSL certificates, making it inaccessible for most users. Automatic certificate renewal (e.g., using certbot) should be implemented.

      4. Attempting to use the SIB UI resulted in an error due to the redirect.

      5. The certificate issue also prevented evaluation of ExpasyGPT querying against METRIN-KG.

      6. Usability represents a significant barrier. Many potential users such as biologists without semantic-web or RDF experience are unlikely to be able to interpret the current figures or formulate SPARQL queries. The manuscript briefly mentions ExpasyGPT, which has strong potential to overcome this barrier by allowing natural language querying. This tool should be emphasized more prominently, potentially with a dedicated subsection and discussion of its role in broadening the resource's accessibility.

      7. To fully demonstrate the relevance of METRIN-KG, the use-case section would benefit from quantitative summaries, visualizations (e.g., distributions, network visualizations), and biological interpretations of query results.

      8. The overall manuscript organisation would benefit from restructuring to improve readability and better expose the impact of the work done.

      9. The methodological content currently in "Mapping of TRY data" and "Mapping of GloBI data" should be moved to the Methods section, while the quantitative outputs (e.g., numbers of records) should be moved into a dedicated Results section.

      10. The current "Data re-use and case studies" section could be reorganised into Results and Discussion sections,

      11. Results include: description of outputs from the methodological steps (e.g. the ontology, the successfulness of the mapping process, the size of the final knowledge graph/number of triples, and other relevant metrics; the user interface, including being able to share and add example questions and write NL questions via ExpasyGPT; examples of SPARQL queries and case studies.

      12. Discussion includes: Reuse potential; case-study interpretations or impact; future directions including planned expansion and enhancing of ontological structure, etc.

      13. An overview and evaluation of the KG should be provided, for example the number of plant species, and the distributions of connections. Any gaps or any (potentially) biases due to the input data needs to be acknowledged. For example, it is very likely certain species may be under/overrepresented in either metabolome or interaction datasets, or possibly geographic skews could exist.

      14. The ontology is variously referred to as the "Earth Metabolome Ontology", "EMI ontology", and "EMI". Consistent naming should be adopted throughout the manuscript and associated repositories. It is also unclear whether the ontology is a result of this work. As written, the "Ontology" section under "Methods" reads more like a result/description than a methodological step. Clarifying what components are original contributions and presenting them in Results would strengthen the manuscript. Additionally, the phrases "our proposed framework" and "our approach" are ambiguous, do these refer to the ontology itself, the metadata-mapping pipeline, or the overall integration process? Finally, referring to METRIN-KG as a "tutorial to build a knowledge graph" appears to be a bit out of place, given the topic of the manuscript.

      15. Given its potential relevance beyond this project, the authors are strongly encouraged to publish the code for the metadata-mapping pipeline. In addition, the following details would strengthen the methodological rigor:

      16. Were any acceptability criteria implemented in the automated step, e.g. a minimum Cosine similarity threshold?

      17. Were the manual corrections systematically documented?

      18. In how many cases were manual corrections needed?

      19. Was any evaluation done on the embedding/model version or the source of errors?

      Minor comments

      1. The manuscript would benefit from proofreading for consistent use of Oxford commas (and an "&" instead of "and").

      2. Reference 190 contains a typo "GiHub"

      3. References should be checked (e.g. citation for [95] references both METRIN-KG Zenodo and GloBI Zenodo)

      4. The METRIN-KG Zenodo link in the article is not to the latest version (Version 5)

      5. Consider improving the figures, e.g. use of colour to better communicate the content and refining layouts.

      6. The authors should deposit a snapshot of the GitHub repository to Zenodo.

      The following are suggestions to the authors, to be followed by their own judgement:

      1. Table 1, Figure 3, Figure 4 could be moved to supplementary material to make space for figures for the case studies.

      2. The dense in-text list of ontologies in the "Metadata mapping" section could be replaced with a summarized table (e.g. by moving Supplementary Table 2, but including references to the main text).

      3. The full SPARQL queries in "Taxonomy mapping" could be moved to supplementary materials, with a high-level description left in the main text.

    2. AbstractBackground In recent years, biodiversity data management has emerged as a critical pillar in global conservation efforts. Today, the ability to efficiently collect, structure, and analyze biodiversity data is central to breakthroughs in conservation, drug development, disease monitoring, ecological forecasting, and agri-tech innovation. However, due to the vastness and heterogeneity of biodiversity data, it is often confined to databases for specific research areas in isolated formats and disconnected from other relevant resources. Crucial components of such data in kingdom Plantae comprise of metabolomes - the vast array of compounds produced by plants; traits - measurable characteristics of plants that influence their growth, survival, and reproduction, and that affect ecosystem processes; and biotic interactions - relationships of plants with other living organisms, affecting the ecosystem functions.Results In this work, we present METRIN-KG (MEtabolomes, TRaits, and INteractions-Knowledge Graph) a powerful data resource simplifying the integration of diverse and heterogeneous data resources such as plant metabolomes, traits, and biotic interactions.Conclusions The proposed knowledge graph provides an interface to interactively search for data relating plant metabolomes, traits, and interactions. This, in turn, will facilitate development of research questions in life-sciences. In this context, we provide representative case studies on how to frame queries that can be used to search for relevant data in the knowledge graph.

      This work has been peer reviewed in GigaScience (see https://doi.org/10.1093/gigascience/giag051), which carries out open, named peer-review. These reviews are published under a CC-BY 4.0 license and were as follows:

      Reviewer 1:

      This study makes an important and timely contribution to plant ecology by combining multiple data sources of plant functional traits, metabolites, and interaction partners. Linking different aspects of a plants' phenotype (morphological and physiological traits, as well as metabolite profiles) with their potential ecological functions (biotic interactions) represents a much-needed step forward in both chemical and functional ecology. The effort the authors have put into compiling these datasets and establishing connections across sources is evident and deserves appreciation. The potential applications of this framework are manifold, and the examples provided on how the database can be used to explore research questions through knowledge paths convincingly demonstrate its value. The Introduction could be strengthened. At present, it feels somewhat long and diffuse in scope, which may make it harder for the reader to quickly grasp the importance of the contribution. Streamlining the text and sharpening the focus on the added value that linking different data sources provides would considerably improve clarity and help the reader more fully appreciate the strength of this novel approach. Minor comments: It would be helpful to specify more clearly what is meant by the "tremendous chemical diversity" that challenges metabolomic data analysis, and to clarify whether the "limited resources available to manage such complexity" refers to analytical limitations, data availability, or both. Similarly, the phrase "multi-level interaction data also lies locked in the guise of pairwise correlation metrics" could benefit from further explanation. As far as I am aware, the GLOBI database reports individual pairwise observations, which could already be seen as an effective way of documenting such interactions.

    1. standard data

      The Qid for descriptive standard / Erschließungsrichtlinie is (Q117460293), now with additional statements and linking items such as Encoded Archival Description (Q1340077), ISAD(G) (Q1654544), Rules for Archival Description (Q22908401) and Records in Contexts–Conceptual Model (Q30216914).

    1. Emotions amplify through a process he called “mutual induction,” where one person’s excitement sparks another’s, escalating rapidly.

      Where do you see this happening?

    1. eLife Assessment

      This is a valuable paper looking at nanoscale organization of the membrane associated periodic cytoskeleton in mouse sciatic nerve axons. Despite previous studies, the precise organisation of the structure remains unclear, especially in vivo, and this manuscript significantly adds to this knowledge with solid data. An unexpected observation is the presence of discrete nanoscale clusters, regularly distributed around sections of axons. However, the paper misses a description of these clusters along the longitudinal axis of the axons.

    2. Reviewer #1 (Public review):

      Summary:

      The article "Nanoscale organization of beta-II spectrin within segments of the membrane-associated periodic skeleton in mouse sciatic nerve axons" by Gazal et al. looks into the organization of the spectrin scaffold in mouse sciatic nerves using super-resolution microscopy. It is now well established that axons, across species, contain a membrane-associated periodic scaffold mainly composed of circumferential actin filaments and longitudinally arranged spectrin tetramers. While super-resolution imaging of neurons in cell culture is relatively easy, exploring the ultrastructure of myelinated axons in intact nerve fibers is a daunting task. Nevertheless, the authors have attempted this by fixing and preparing cross-sections of sciatic nerves. They have then tried to quantify the fluorescence intensity patterns of specific components, especially that of labeled beta-II spectrin and have analysed its distribution.

      One of the main findings is that spectrin is distributed along the axonal periphery and along the outer part of the myelin sheath. By labelling multiple cellular components and using intensity analysis, the authors show the sequence of structural organization of a few key components. They see that, unlike in the case of axons in culture, the axonal cross-sections within the sciatic nerve deviate significantly from a circular shape. They then use 3D-dSTORM to investigate the distribution of beta-II spectrin along the axonal circumference. They see that this distribution is very heterogeneous, both in the sizes of spectrin puncta and their arrangement along the periphery. The amount of spectrin scales linearly with axonal circumference.

      Strengths:

      Super-resolution imaging of axons of intact nerve fibers to investigate the organization of beta-II spectrin.

      Weaknesses:

      While most of the findings, like the spatial distribution of spectrin and related components, are reasonably well supported by data, I have concerns regarding the subsequent claims made in the article. The detection of axial periodicity based on the observation of a peak in the inter-tetramer spacing distribution is not very convincing, and a 3D representation (or a video of 3D reconstruction) would have been better. And so are the claims on characteristic spectrin spacing of 200 nm along the axonal circumference. A peak in the distribution does not imply a periodic arrangement.

    3. Reviewer #2 (Public review):

      Summary:

      This is an interesting paper by the Unsain lab looking at the nanoscale organization of the membrane-associated periodic cytoskeleton in mouse sciatic nerve axons. The precise organization of the structure remains unclear, especially in vivo, and this manuscript significantly adds to our knowledge of this important structure. While some of the findings in the study are somewhat expected (though still valuable to see in an in vivo setting), an interesting observation is the presence of discrete nanoscale clusters that scale up with the size of the axon, which challenges previous assumptions.

      Strengths:

      Strong, convincing data; clever combination of imaging and analytical tools to make novel points; well written; excellent composition of figures.

      Weaknesses:

      (1) Figure 2A/3A: The large and small clusters of spectrin, as seen in cross sections, are unexpected and novel. The authors have done a clever job of combining imaging and analyses, but some things are still unclear. First, the authors should be consistent in their language when they talk about the spectrin clusters. Recommend precise language to define the small and large clusters when they first appear in the text, and then use the definitions consistently throughout the text. Second, based on the data shown, one does not get a clear idea of how the small and large clusters are organized along the longitudinal axis of the axon. In that context, are Figure 2B and C from imaging along the longitudinal axis? If not, it's unclear how the authors can conclude that the spectrin assemblies have a distance of ~170 nm along the linear axis. In general, a perceived limitation of this study is that while the authors have done a good job looking at cross sections, there is no information on the longitudinal distribution of spectrin in these axons. Looking at both cross- and longitudinal sections would also clarify details about the large spectrin clusters. For instance, are they small sausage-like structures, or long rods of spectrin running along the length of the axon? One assumes that all the analyses in Figures 3 and 4 are from the small clusters. Can the authors do a similar analyses of the large clusters? Finally, a schematic model showing both cross- and longitudinal- sections would make things clearer, but the authors would need to show the longitudinal data for that.

      (2) It is interesting to think that the larger spectrin accumulations may be similar to the condensate-like structures seen by Boyer et al., as the authors mention in the discussion. In that context, it is possible that these focal accumulations are local reservoirs of spectrin that are also seen in mature axons (indeed, these accumulations were also seen in mature axons in the Boyer et al. paper, and they also speculated that these accumulations may be local reservoirs). Can the authors check if actin/adducin is also present in these larger spectrin accumulations?

      (3) While talking about the nanoscale clusters, it is important to specify that the authors are talking about circumferential clusters. Though the writing is excellent, one still does not get the precise definition of "clusters" from just reading the abstract, and it would be good if the authors could work on that more (I recognize that this is not easy to do).

    4. Reviewer #3 (Public review):

      Summary:

      In the presented work, the authors investigate spectral staining in axons of the sciatic nerve, where the MPS has been detected before using STED microscopy. They employ 3D-dSTORM in tissue sections and analyze the data, measuring localization of clusters on the axon perimeter and the relative distribution of those. From these data the conclude that large gaps in spectrum localizations exist and that clusters around the axon exist that are spaced at 200nm.

      Major Comments:

      (1) The presented data are at times overinterpreted, and the discussion lacks a critical view of the data. For example, the statement "...Unlike previous suggestions from qualitative evidence in cultured neurons (REfs), βII‑spectrin distribution in MPS segments of peripheral nerves is discontinuous, with extensive stretches of the perimeter lacking βII‑spectrin." is quite strong, given it is based on immunofluorescence staining and dSTORM microscopy in tissue. Absence of evidence of staining is not evidence of absence.

      (2) The authors claim in the abstract that "The number of these clusters scales linearly with the axonal perimeter, maintaining a constant membrane occupancy of ~20% across varying axon diameters." Again, this is from a cut through an axon, while measuring the density of clusters on the perimeter. If they claim area occupancy, an area should be imaged, and the dots (clusters) should be measured in surface coverage in a 2D projection of the axonal surface.

      (3) In general, this reviewer suggests being a bit more moderate in statements such as: "These findings challenge simplified models of the MPS based on cultured systems and demonstrate that the MPS in peripheral nerves is composed of discrete structural units." These statements are bold from the relatively few measurements in a single method and a single viewpoint. Especially when considering that techniques such as dSTORM depend extremely highly on labeling density, and apparent clustering of localization is highly prone to misinterpretation. If the authors desire to make such statements, working with endogenously labeled protein would be warranted. The authors should at least hedge such statements.

      (4) If the authors want to make statements about general organization, why do they not compare adjacent cuts through the axon? If there are continuous spectrin filaments, the clusters should appear at the same site across repeated cuts through the axon.

      Besides this, this reviewer welcomes the effort that has been made to establish dSTORM in tissue sections and to investigate the MPS in native tissue.

    5. Author response:

      We sincerely thank the editors and reviewers for their overall positive assessment and constructive feedback on our manuscript detailing the nanoscale organisation of βII-spectrin of the membrane-associated periodic skeleton (MPS) in mouse sciatic nerve axons. Their perspective and comments will help refining the manuscript.

      A common comment by the reviewers relates to the description of the characteristic longitudinal periodicity of the MPS. We value these comments, which we believe are motivated by the fact that the longitudinal periodicity of the MPS is undoubtedly the most studied and prominent feature of the MPS in cultured neurons. However, the main goal of the present project was to describe how βII-spectrin is organised in the transverse axis of individual segments of the MPS in nerve tissue. This is why we utilised cross-sections of the sciatic nerve, hence achieving the best resolution possible in that plane, at the expense of the resolution in the axial axis. Furthermore, this study clearly shows that the transverse morphology of axons, and thus of the MPS, of neurons in the tissue is highly irregular, in comparison to cultured neurons. This imposes an extra challenge to observe correlated longitudinal structures when the observation length is limited, as in our studies. Nonetheless, to improve this aspect of the manuscript, we will revise our data and previous evidence, clarify the methodological trade-offs made, and make our interpretations more accurate.

      Additionally, we will clarify several imaging- and definition-related inquiries, including tests for insufficient staining, the interpretation of βIII-tubulin staining, the assessment of axon–glia boundaries, and consistency in the use of terms like ‘clusters’ and ‘periodicity’, among others.

      We believe these and other revisions will substantially strengthen the manuscript and comprehensively address the reviewers' feedback.

    1. eLife Assessment

      This fundamental manuscript describes a key role for the integrated stress response-regulated transcription factor CHOP in regulating liver biology in response to endoplasmic reticulum stress through both the downregulation of transcription factors involved in regulating hepatic identity and altering the capacity for integrated stress response and unfolded protein response signaling to induce protective signaling. The data supporting this model is convincing, but including some additional discussion on the mechanism and importance of the work in the context of the published literature would be helpful to better define the complex importance of CHOP signaling. This work will be of interest to a wide range of biologists interested in liver biology, stress-responsive signaling, and ER stress.

    2. Reviewer #1 (Public review):

      Summary:

      The predominant view on CHOP's functions during ER stress is that it promotes cell death. This is in contrast to a handful of reports in the literature that claim that CHOP is a positive regulator of protein synthesis during chronic ER stress, and therefore is part of the adaptation program to ER stress. These previous studies were performed in tissue culture cells. Velarde and co-authors have used a mouse model of induction of mild ER stress to study the function of CHOP in hepatocytes.

      Major strengths and weaknesses of the methods and results:

      The authors use state-of-the-art mice to manipulate (i) CHOP and (ii) ATF6, a protective factor of ER proteostasis, and address the hepatocyte responses to mild ER stress in vivo and in cultures. Validated gene expression programs are well correlated to liver pathology in the mouse models. This is a very well-done study.

      The authors clearly show that CHOP transitions hepatocytes under mild ER stress to a chronic ISR state, which is phenocopied by ATF6-depleted hepatocytes. So the conclusion that CHOP exacerbates ER stress in hepatocytes during mild ER stress is correct. It is also clear that CHOP targets negatively the transcription of hepatocyte identity genes, which opens a new direction of studies on the function of CHOP in secretory cells in general.

      Conclusion:

      This is a significant study that will benefit different research fields, and specifically studies on proteostasis, as was recently highlighted in Nat. Str. Mol. Biol. by experts in the field.

      To this reviewer, the importance of the study is that it links the function of a transcription factor (CHOP) to stress intensity (mild versus severe) in a physiological experimental model (hepatocyte function and pathology).

    3. Reviewer #2 (Public review):

      The Unfolded protein response (UPR) and related integrated stress response (ISR) are critical signaling systems for cell survival in response to acute stresses. While the UPR directs critical adaptive gene expression, certain chronic stresses switch this pathway towards cell death and disease. An important question concerns the mechanisms by which the UPR switches from being adaptive to maladaptive. Prevailing models focus on the transcription factor CHOP (DDIT3 or GADD153), whose levels are enhanced via the UPR, and extended/amplified amounts of CHOP are suggested to boost death-related gene expression. However, the literature and this manuscript point out a number of observations that do not neatly fit with this model, suggesting that there are still unresolved processes by which CHOP adjusts cell outcomes via the UPR.

      This manuscript features a nice hepatocyte-targeted knockout of CHOP to discern the contribution of CHOP in the transition between adaptive and maladaptive outcomes. The key ideas presented in this study are that CHOP-directed gene expression is focused on protein synthesis, metabolism, and hepatocyte identity. In the progression of the UPR, CHOP expression can lead to resumption of protein synthesis, which can assist in the translation of the UPR-directed transcriptome, which includes ATF6/XBP1-directed genes that aid the processing capacity of the endoplasmic reticulum (ER). However, enhanced nascent protein can further stress the ER. CHOP directs gene expression in both the first phase- acute and second phase-chronic in the UPR, and the pivotal decision lies in the transition between the phases.

      Overall, the manuscript includes some new ideas as well as refinements of earlier ones for CHOP-determination of UPR-directed cell fate. The CHOP-hepatocyte knockout mouse model helps to delineate the different tissue functions of CHOP, which has been a problem for some earlier studies. The manuscript progression of experiments is solid, and experimental design and documentation are rigorous. The manuscript text is largely clear, but there are portions that would benefit from fuller explanations of ideas.

      There are three points of concern. First, the manuscript model (Figure 7) lays out a timeline for the progression of the UPR between two phases. The study is not always clear about the times assayed, and there appears to be a single time point for measurements. Second, there is emphasis on protein synthesis changes in the model. It is true that the literature argues that resumption of protein synthesis concurrent with stress damage (i.e., GADD34-directed gene expression) is a key reason for the potentially debilitating effects of CHOP (e.g., Marciniak et al 2004, Han et al 2013). However, the manuscript does not feature protein synthesis measurements. Inclusion of bulk protein synthesis measurements in the context of this model system would strengthen the study and support for the model. Finally, for this reviewer, some of the most interesting ideas center on CHOP-directed transcription of genes that regulate hepatocyte identity. There is solid evidence for direct CHOP regulation of these genes, but the manuscript does not really develop and test the ramifications of these networks on cell fate during ER stress.

      Reviewer Concerns:

      (1) The abstract packs in a lot of information. The ideas would not be clear to a general reader. Furthermore, the UPR and ISR are referred to in the second-to-last sentence, but not defined earlier in the abstract.

      (2) There are some typos/grammar concerns.

      (3) ATF4 diminished with CHOP-depletion (Figure S2A). What is the mechanism here? Does this complicate the analysis of CHOP-directed gene expression? How does this fit with Figure 6J? The timelines for TM treatment are critical. The authors should more fully explain the time courses in the experiments.

      (4) Figures 2 and 3: There is a discussion on enhanced protein synthesis with loss of CHOP (reduced GADD34 expression). What is the time point - 8 hours TM? Emphasize, explain, and justify time points of experiments here and in later panels. It would strengthen the model with direct measurements of protein synthesis. The authors could include GADD34 protein measurements in these panels. Figure 3 - panel D - some abbreviations are not standard.

      (5) Figure 4: One of the most interesting in the manuscript is the transcription factors downstream of CHOP that are linked with hepatocyte differentiation and metabolism. The manuscript would be bolstered by developing some of these target genes into the Figure 7 transition model.

      (6) Figure 6: The comparison of CHOP and ATF6 target genes is a highlight of the manuscript. The literature on this topic is complex, and there are some suggestions that CHOP can be downstream of ATF6. Furthermore, there were some earlier models by Walter and others about extended induction of Perk (death) vs induction of other UPR sensors (survival) (e.g. PMID: 17991856). It would be helpful in the Discussion to delineate between these models and their critical differences.

    4. Reviewer #3 (Public review):

      In this manuscript, the authors aim to understand the function of the transcription factor CHOP, which is known to promote cell death during severe stress in the ER. The authors note that CHOP is induced during less severe stress, but its functional output is not well understood in these cases. Here, they study the effects of conditional knockouts of CHOP in hepatocytes of mice challenged with chemical inducers of ER stress.

      Tunicamycin (an ER stress inducer) injection leads to the upregulation of CHOP and lipid accumulation in the liver, but no significant cell death in the experiments outlined here. Conditional knockout of CHOP results in a number of differences in the way hepatocytes respond to stress, notably resulting in lower steatosis.

      There are two main findings supported by the data presented here. First, the authors show that CHOP suppresses the expression of ONECUT, a master regulator of hepatocyte differentiation and metabolism, during ER stress. They show by ChIP-seq that CHOP binds to the promoter region of this gene, and by RNA-seq that ONECUT expression is suppressed by ER stress in a CHOP-dependent manner. Many predicted targets of ONECUT1 were also suppressed by ER stress in a CHOP-dependent manner, though they were not bound directly by CHOP. The data support a model where CHOP down-regulates hepatocyte metabolism and identity via regulation of ONECUT1. This is a new and interesting finding, perhaps explaining the steatosis phenotype of livers that accompanies ER stress, although this was not tested directly.

      The second main finding of this paper is that CHOP deletion leads to an interesting assortment of effects on genes related to the ER stress response and integrated stress response (ISR). As expected, based on prior work, CHOP deletion led to more phosphorylation of eIF2alpha (CHOP is known to upregulate the phosphatase for this translation factor). However, unexpectedly, this did not cause increased expression of ATF4 (a transcription factor whose upregulation during stress is dependent on eIF2alpha phosphorylation) and its downstream targets; in fact, CHOP deletion had the opposite effect on these. In other words, CHOP seems to both turn off the initiating signal for the ISR (namely, eIF2alpha phosphorylation) and also promote the downstream signaling events that rely on this initiating signal. It makes sense that cells would do this, as restoring translation would be important for realizing the effects of the massive changes in gene expression initiated by ER stress, and yet this would exacerbate stress in the short term, so it would be counterproductive to also turn off the entire stress-regulated program. Having a factor (perhaps CHOP) that coordinates these two events makes sense. It will be interesting in future work to understand the mechanisms behind this regulation.

      Finally, CHOP deletion led to less activity of other aspects of the ER stress response, notably IRE1 (determined through measurement of XBP1 splicing and RIDD of Bloc1s1). This is explained by the continued phosphorylation of eIF2alpha in these knockouts, as the continued attenuation of translation would lessen the burden of misfolded proteins in the ER. Somewhat confusingly, the same pattern is not seen in downstream targets of XBP1. Less splicing, coupled with perhaps less translation of the spliced mRNA, should result in less active transcription factor and lower expression of its target genes in the CHOP KO. This is not observed in Figure 2, although the more global gene expression analysis suggests that all stress-dependent gene expression changes were weaker in the CHOP KO livers.

      The authors characterize the effects of CHOP, promoting restoration of protein synthesis and the accompanying exacerbation of stress while preserving the signaling that should relieve ER stress, as a switch from an acute to chronic phase of ER stress. This is mirrored in their analysis of ATF6 in a similar series of experiments. Although this is an interesting framework for thinking about the stress response, whether CHOP is the key factor or a supporting actor in regulating this transition will require a better understanding of the mechanisms involved.

    1. 🤗 Kernels: Major Updates

      ┌─────────────────────────────────────┐ │ 模型库层:Transformers / Diffusers │ ← 用户写模型代码 ├─────────────────────────────────────┤ │ 🤗 Kernels(内核分发与加载层) │ ← 新增的这一层 │ - 从 Hub 拉取预编译内核 │ │ - 匹配当前硬件/OS/框架版本 │ │ - 安全验证(签名、可信发布者) │ ├─────────────────────────────────────┤ │ 框架层:PyTorch / JAX / CuPy │ ← 张量运算、autograd ├─────────────────────────────────────┤ │ 厂商运行时:CUDA / ROCm / Metal │ ← GPU 编程接口 ├─────────────────────────────────────┤ │ 硬件:NVIDIA / AMD / Apple GPU │ └─────────────────────────────────────┘

    1. Finally, the role of advocacy withinthe queer community and whether queer peoplewith more privilege advocate for those with less.

      The stratification of privilege in the queer community is an important thing to remember. Folx who break away from heteronormative self expression experience more violence than others. The wheel of privilege helps us locate both ourselves and other identities to approximate the priveleges or unearned social benefits of people.

    2. 2000-2010

      It seems so bizarre at first but then if you think about the social climate at the time, it makes perfect sense why folks did not feel safe about publishing at this time period.

      I recently read Malcolm Gladwell's Revenge of the Tipping Point where he discusses the mid-90s resurgence of homophobia around the time that Bill Clinton signed the Defense of Marriage Act (DOMA) in 1996 after the don't ask don't tell policy in 1993 allowing members of military service to serve, so long as they conceal their identity. Gladwell goes on to discuss how the sitcom Will & Grace premiered in 1998 and the intentional characterization of Will as largely heteronormative with a deeply caring relationship with the other main character, Grace. Will's character was a counter-narrative to the pervasive homophobic societal expectation at the time that gay people were not capable of deep relationships in a same sex union. Gladwell's point in the book was that it took only a few years for a critical mass to shift (the tipping point he describes as a critical third) and we see that playing out in 2003 under George Bush, that the United States Supreme Court striking same-sex intimacy. While in 2004, Massachusetts became the first state to permit same-sex marriage. (American Queer History)

      In 1992, Canada ended the prohibition on gay and lesbian service members, and in 1995 The Supreme Court recognized sexual orientation as a protected ground under the Canadian Charter of Rights and Freedoms. By 1996, sexual orientation was added to the Canadian Charter of Rights and Freedoms. By 2005, the Civil Marriage Act provided for the legalization of same-sex marriage across the whole of Canada.

    3. Articles from non-peer-reviewedjournals, magazines, and books

      I understand that this is a journal article and perhaps needed peer-reviewed articles for submission to this journal (???) but I think there can be merit to using non-peer-reviewed articles, especially for a topic like queerness. From a participant perspective, someone might feel more comfortable expressing themselves in a blog or book format, as opposed to a study, due to power dynamics.

    4. white, cisgender,

      Both author's have acknowledged their social positions as white and cis-gender. It seems like culture is largely absent from this literature review and only mentioned a few times. As we have been discussing intersectionality and our positionality this week, I wonder whether acknowledge of our social identities, and perhaps what we cannot personally speak to, warrants excluding articles from authors that can speak to that perspective?

      I personally think that it is not difficult to find authors who have written about the perspectives we cannot personally speak to, but I'm curious to know other folks opinions!

    5. task

      I agree that academic leaders should focus on queer student, but also wonder whether there should also be a focus on queer staff? In my experience, queer students feel more comfortable and empowered when queer staff also feel comfortable to share their identity and experiences. Lee (2022) did a case study highlighting how "sign" and "symbols" of queerness were important in helping queer staff feel comfortable in the workplace, which then allowed them to be support for students. This study specifically mentioned that senior leadership visibility was important to staff empowerment.

    6. White, gay people find lessdifficulty being out (

      There is actually quite a lot of racism, transphobia and sexism in white gay male spaces. White people and white culture are often treated as the default in queer spaces. Because of that, queer and trans people of color frequently feel excluded and have to push back against that logic in both personal interactions and institutions.

    7. In essence, the messagefrom universities was a cold indifferenceto queer people and a genuinereluctance to support or even includethem in programming or policies

      This is happening right now in North America. Withdrawing of support, cutting funding, silencing, looking the other way etc. Trans people in particular are under direct legal, social and cultural attack from the far right. This started in the US, and now is happening in our more conservative provinces.

    1. Crowds are dynamic. They are adaptable. But perhaps their most important and least appreciated feature is that they are a pathway to understanding our social selves.

      Do crowds help us better understand ourselves?

    2. “The collective is so important in forming our everyday identities

      Do you agree with this statement? Who is the main crowd / collective in your life? How are you impacted?

    3. individuals in crowds do not abandon their rationality or surrender their identity to a mob mentality. They do not lose their minds. They do, however, become highly sensitive to what those around them are doing, and become strongly cooperative as a result.

      How is this different from the traditional notion of how a mob behaves? Do you agree with this view?

    4. Most of us seem convinced that crowds inhabit a psychological shadow land of primitive urges and unrestraint, where individuals are stripped of their identity and led unthinking to violent or irrational acts.

      Traditional notion of mob mentality

    1. SQLite 数据库的路径,或 :memory: 使用内存数据库。如果 filename 是空字符串,那么将创建私有的临时磁盘数据库

      三种不同的行为

    1. And so, dear readers, I take theunusual step of inviting you to write inwith your views. After all, you have astake in the issue

      The author says that there doesn't have to be one answer to this question and that us readers should be able to have our say in it.

    2. asPHUHO\HFFHQWULFLQKLVKDELWVbbbLVQRmore sensible than adopting the sameattitude to a child who gets the wronganswer to his sums.”

      Treating a person poorly because they didn't use proper grammar is like treating a kid badly because they didn't get the right answer on a math equation.

    3. Ifyou speak Standard Mandarin, peoplewill always assume that you went to agood university.

      This means that knowing the standard language can actually help someone get a better job, not just speak and write more formally.

    4. “Nonsense,” they collectively said.“Just announce that this class is beingconducted in Standard English, anddemand Standard English.

      At fist I was surprised they said this, but it makes sense, because it is what they've been taught their whole lives.

    5. They can do whatever they wantoutside class,”

      The students think that the the classroom is only for standard english and other dialects are supposed to stay outside.

    6. As children learn Standard English—becoming bidialectal—teachers mustn’tdenigrate whatever dialects are spokenin their homes.

      This essay is already starting to sound similar Young's argument! He also argued for keeping the english that people spoke at home.

    Tags

    Annotators

    1. dministrator diversity gap

      As a vice principal, I have very little power within the board and even less power within the hiring process. I think it's important to address this issue of lack of diversity within our schools because in elementary schools especially, teachers are mostly homogenous group. However, I think when we talk about increasing diversity in school administration, we need to include senior administrations (superintendents) as well. They have more power within the board, and in our board, they make all the hiring decisions. I often find within the school boards, the further up the power chain we go, the more homogenous the group tends to be,

    2. discriminatory practices and microaggressions

      This has been my experience with the local public board. My colleague and I created East Asian Educator Affinity Network at my board to support and network with East Asian educators. The stories shared during our first meeting were heart-breaking. We also felt very invisible as a group within the board. To give an example, the board sends out monthly list of holidays and celebrations from many religious and ethnic groups to all the schools. Yet, they didn't include Chinese New Year, a very obvious but important East Asian holiday.

    3. pervasive nepotism and favouritism

      Given Regulation 274/12 was introduced in 2012, I am curious about how nepotism and favouritism were still allowed to continue when this report was conducted in 2020. When I worked for Peel, teacher hiring was strictly done through seniority list and principals no longer had any input. Even with the Regulation revoked, some boards are still continuing to use the seniority based hiring, not that I am suggesting that solves this issue of teacher diversity. I think lack of teacher diversity comes from a much wider systemic issues.

    4. permanent teachers in Ontario are White,

      At least in TDSB and Peel, there were many BIPOC teachers. However, this was not the case for administrators. I always wondered why this was the case. I took my PQP online asynchronously and everyone in my class came from different parts of Ontario (and some outside of Ontario). Yet, the class was not very diverse. In comparison, I am really enjoying how diverse our cohort is for the PhD program where so many different perspectives are brought together because of the diversity. I just wonder what the difference is between these two programs.

    5. Equity and inclusive education policie

      At times, this EDI policies feel very surface level to me. At Greater Essex County District School Board, they spent an enormous amount of money on surveys and reports to address anti-black racism but beyond the surveys and reports, nothing more came out of it - similar to how so many projects were unaddressed in the Truth and Reconciliation Commission's 94 Calls to Action.

    6. Aspiring Racialized Leaders Mentoring Project

      I am not familiar with this Aspiring Racialized Leaders Mentoring Project but we cannot put the onus of bringing more diverse teachers and administrators to BIPOC...simliar to how Gaudry and Lorenz talked about not requiring Indigenous people to bear the responsibility for change in Indigenous inclusion policy.

    7. economictrends

      I think some of the teacher hiring processes were the results of political and economic climate and its responses from the union. For many years, there were very few teaching jobs which meant that teacher graduates who needed steady income to survive did not stay in teaching. It used to be the case that when teachers were hired at the board, they stayed in the profession for many decades. There was a teacher shortage at the beginning of 2000 when there was a change to retirement requirements (going from 90 factors to 85 factors). So a large number of today's teachers were hired around then where teacher diversity was not considered.

    8. region to region.

      During my master's thesis I learned that I really knew nothing about education system in other provinces. In fact, being an elementary teacher, I really knew very little about what was happening in secondary schools, either. I think this is a disservice to teachers as it encourages teachers to hold a very narrowed view of education. As well, I've learned over the years that the comprehensive and cohesive educational strategies differ not just region to region but swings rather massively based on the political parties and time. For example, for many years, reading was taught through whole language approach and now we are moving back into phonics phase.

    9. that BIPOC educators are underrepresented in per

      In fact, this often leads to Black, Indigenous, and other racialized students to being pushed into “special education” streams and disproportionately ignored, punished, disciplined, or expelled compared with white students. Yet many schools refused to confront the systemic racism driving these patterns, choosing instead to blame the students or their parents for the challenges they face (Gillborn, 2015; James, 2019).

    10. ntario’s lack of teacher diversity h

      The vast majority of teachers in Canada are white women. Gebhard (2020) discusses this as an issue by asserting that “the role of white women in education has historically been to reproduce rather than disrupt colonial epistemologies and discourses of white dominance as these intersect with white womanhood” (p.207). This continued system of ‘keeping the status quo’ may be partly due to ignorance and partly because thinking about or confronting systemic racism and white supremacy is deeply distressing to many people, particularly those in a position of power.

    11. particular Black male studen

      I am going to leave this short piece I wrote in 2023:

      Black children are not scary. They are scared.

      Over and over again, we see news about black children being arrested, beaten or murdered on our streets. There is no respite in the school system either, where black children are neglected, punished, and expelled disproportionally over white children. As Munroe (2021) says so powerfully, “within all of these conversations about racist practices, what often gets left out are stories about the ways school disciplinary practice and policies affect Black students’ emotional well-being and traumatizes them. Black students continue to be judged as inferior and dangerous. Even with ‘progressive discipline,’ Black students who violate rules are seen as offenders rather than teens in need of support” (para. 4).

      It is very clear that black children are not scary. They are scared. Cole (2020) goes on further and very plainly outlines the emotional and systemic consequences of suspending and expelling black children from our schools: The gun violence issues in the city of Toronto have been one of the number one stories in this city for the past two years. Why does no one make the connection to poverty in the city of Toronto and to the expulsion of black students, which further alienates them. We've heard so many times the impact of not giving kids a place at school where they feel welcome or they feel like they have an environment to learn. We push kids out of the education system and we have this huge epidemic of violence and anti-social behaviour (para. 32). This means that keeping kids in school is of utmost importance.

      In addition to keeping them in school, black children need support, care, love and respect. This means teaching black histories, contributions and perspectives and integrating them year-round, including black authors, resources, and contributions in the classroom (Facing History & Ourselves, 2021). This also involves celebrating and appreciating black excellence and cultural contributions in art, music and literature. Love (2014) says ignoring students' culture is an injustice and an oversight. Welcoming black joy into the classroom by allowing students to be themselves and share their culture will enable them to feel seen and appreciated. Emdin (2016) adds to that assertion and further emphasizes that “the classroom that reflects the culture of youth while highlighting how beautiful it could be challenges the ways that the teacher and students engage with each other even before teaching begins” (para. 17).

      References

      Cole, D. (2020, February 6). What’s it like to be black in Canada? Under policing, it’s hell [Interview]. In cbc.ca. https://www.cbc.ca/radio/thecurrent/the-current-for-feb-6-2020-1.5454037/what-s-it-like-to-be-black-in-canada-under-policing-it-s-hell-says-desmond-cole-1.5454039

      Emdin, C. (2016, March 22). For White Folks Who Teach in the Hood (L. Ferlazzo, Interviewer) [Interview]. In Education Week. https://www.edweek.org/teaching-learning/opinion-for-white-folks-who-teach-in-the-hood-an-interview-with-chris-emdin/2016/03

      Facing History & Ourselves. (2021). Black History Month Resources: Approaches, Identities, Histories, Legacies & Inclusion. Facingcanada.facinghistory.org. https://facingcanada.facinghistory.org/black-history-month-resources

      Love, B. (2014). Hip hop, grit, and academic success [Video]. In YouTube. https://www.youtube.com/watch?v=tkZqPMzgvzg

      Munroe, T. (2021, February 17). How to curb anti-Black racism in Canadian schools. The Conversation. https://theconversation.com/how-to-curb-anti-black-racism-in-canadian-schools-150489

    12. hiring

      Hiring processes, though presented as neutral, can reproduce racial bias through credentials, geography, and names. Also, there is a policy contradiction: diversification efforts without anti-racist structural change may reinforce inequality rather than reduce it. Administrator autonomy, then, demands critical self-reflection and accountability.

    13. Whiteness

      I found the emphasis on "Whiteness” and White privilege as a self-reinforcing system a significant idea. According to the article, educational systems validate Whiteness and White privilege and, as a result, promote policies and practices that further oppress BIPOC. This idea also promotes a discussion regarding how the dominant cultural norms that the benefactors are often unaware of, actively manifest obstacles and sustain inequities, rather than being a background that has no impact on the situation.

    14. acially stratified society

      Canada is seen as standing for racially stratified society, where Whiteness is privileged. That's especially the case in resource distribution, for instance in employment and education. Canada is different from other countries as it lacks a national education framework. This means that the construction of educational policies and regulations is the responsibility of the province. This allows each province, including Ontario, to set its own policies in regard to education, teaching, and educational administrators which may lead to different levels of diversity.

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

      Learn more at Review Commons


      Reply to the reviewers

      Reviewer #1:

      Major comments:

      1. Lines 103-116 (first paragraph of the results section) describe mainly published data that is more suitable for the introduction section. It is annoying to refer to different published articles in the Results section to strengthen the results instead of showing them. The same goes for paragraphs two and three. Why mention those data in the Results section if they are already published and known?

      We have reorganized this material by moving some background information to the Introduction. Our intention was not to incorporate published data to strengthen our results, but rather to provide essential context for interpreting our findings. We have therefore left some of this foundational information in the results section to create a clear narrative flow, enabling readers to understand the basis for our experimental design and interpretations without needing to recall details from earlier paragraphs in the Introduction. For example, we considered it crucial to restate the earlier report of the BiP:sfGFP:HDEL phenotype in Atlastin mutants, since our results supporting luminal ER protein displacement contradict the previous fragmentation model.

      The following concept was in line 103 in the Results section, and is now in the introduction in lines 82-92: "Conventional light microscopy, commonly used in studies of neuronal ER structure, lacks the resolution necessary to visualize individual ER tubules in small structures, such as presynaptic terminals. The ER is highly sensitive to fixation, and live imaging experiments in neurons in vivo have been conducted on upright microscopes using water dipping objectives with a typical axial resolution limit of >300 nm, which cannot distinguish the densely packed ER tubules at presynaptic terminals (3,8,21,28,39-42). Electron microscopy offers higher resolution, but cannot be used in live samples and has typically been limited to thin 2D sampling (in which it is difficult to distinguish ER cross-sections from synaptic vesicles) (8,20,22)."

      Figure Legends-(in all Figures): The number of experimental repeats must be mentioned in the figure legends.

      This information is provided in Supplementary Table 1, which contains detailed information about the genotype, statistical analysis, and number of larvae and NMJs analyzed. If the journal requires this information in figure legends, we can move it.

      The way the figures are labeled is worrisome; supplementary figures are not ordered numerically.

      We will be happy to rename supplementary figures according to journal guidelines.

      The tubule extension in Figure 2D is not convincing. Is there a movie showing those changes? Better images are needed. It is essential to show which supplemental movie corresponds to which panel.

      We have now included a corresponding video of the same neuron used as an example of tubule extension. We also added another frame to the figure to provide further information on the tubule event we captured. (Figure 2D, Movie S10)

      This is unnecessary in the results section: "To investigate the relationship between ER structure and function at synapses, we examined mutants of Atlastin, a GTPase that regulates ER tubule fusion. Drosophila has a single homolog while mammals have three Atlastin homologs, with Atlastin-1 enriched in the brain (Rismanchi et al., 2008)."

      This information was moved to the introduction.

      "This reduction in ER membrane marker intensity has also been observed in other HSP mutants, suggesting this is a common feature of ER shaping mutants and could indicate changes in ER membrane composition, integrity, or tubule thickness (Perez-Moreno et al., 2023)." This comparison is important and should be shown in the same settings as for the Atlastin mutant rather than referring to published data.

      We agree with the reviewer that it is important to determine whether other ER-shaping proteins, besides Atlastin, also show a decrease in tdTomato:Sec61b to support our claim that this could be a common feature among ER-shaping mutants. To do this, we examined mutants of another ER-shaping protein, Reticulon 1, which regulates membrane bending and stabilization in ER tubules. These loss-of-function mutants were a gift from Dr. Cahir O'Kane at the University of Cambridge and were used in his lab's Pérez-Moreno et al., 2023 publication. We found that in our hands tdTomato:Sec61b levels were reduced in Reticulon 1 mutants, consistent with the results reported by Pérez-Moreno et al. (2023). These results are in Figure 3E-F. We also examined the synaptic distribution of the luminal ER marker, BiP:sfGFP:HDEL, in Reticulon 1 mutants to see if it is displaced to the cytosol. Notably, it remained ER-associated, unlike in Atlastin mutants. These results are in Figure 6F-G, results lines 267-270, and discussion lines 542-545.

      Does the distribution of the luminal ER marker in Figure 6F diffuse due to mislocalization or reflux after being localized to the ER and then refluxed to the cytosol as was previously shown for the ER to Cytosol signaling (ERCYS) mechanism? Could you assess other ER-luminal protein localization biochemically? It is highly recommended to look at another soluble ER-protein localization in the Atlastin mutant without overexpression, which can be an artifact.

      ER stressors can induce ERCYS, in which some luminal proteins, including PDIA3, DNAJB11, ERp29, and an eroGFP reporter, reflux by 30-70% to the cytoplasm without subsequent degradation (unlike ERAD (ER-associated degradation). This phenomenon has only previously been observed in yeast and glioblastoma tumor cells from mice and human . We believe that our work provide the first suggestion that this may occur in neurons, and particularly in a neurological disease model.

      We do not believe that the reflux phenotype for BiP:sfGFP:HDEL is due to its overexpression for two reasons: (1) we observe reflux in our neuronal Atlastin knockdown experiments, even when the levels of BiP:sfGFP:HDEL are significantly reduced artificially because of titration of the GAL4 between the RNAi and the reporter (Figure 7A), and (2) BiP:sfGFP:HDEL overexpression somewhat suppresses endogenous BiP upregulation ((Figure 10 and see Reviewer 1.10), arguing that the transgene does not induce ER stress). We included a new "limitations of the study" section to be transparent about the caveats of the BiP:sfGFP:HDEL reporter (lines 639-664).

      Identifying potential endogenous neuronal ERCYS substrates in our in vivo preparation poses several challenges. First, biochemical approaches, such as fractionation, are not possible in our complex in vivo sample because neuronal ER proteins would mix with ER from other tissues upon homogenization. Second, detecting endogenous proteins with antibodies requires fixation and permeabilization, which notoriously disrupts ER structure and even causes our reporter BiP:sfGFP:HDEL to collapse from a smooth distribution, as visualized by live imaging and FRAP, to a punctate distribution. Third, using antibodies rather than neuronally restricted transgenes makes it challenging to determine whether the signal originates from the neuron or from dense ER structures in the surrounding muscle. Fourth, some ER luminal proteins can displace as little as 30% in the ERCYS examples cited above, and the sensitivity of our imaging assays may limit our ability to detect these small changes. Finally, the limited availability of tagged transgenes and antibodies specific to Drosophila luminal ER proteins (see next paragraph) poses additional challenges. These limitations highlight the need for future studies to develop novel tools and techniques to more definitively test whether we are indeed observing ERCYS. We have included a paragraph on these future challenges in our discussion in lines 639-664. Identifying endogenous targets of ERCYS in fly neurons is a worthwhile goal, but beyond the scope of the current study. These next steps will particularly benefit from identifying the machinery involved in the reflux of our BiP:sfGFP:HDEL reporter.

      Tools we tested: We investigated several options: (1) a tagged PDI transgene (a gift from Karen Hibbard), which was not detectable at presynaptic terminals, (2) a tagged BiP (FlyORF; F000956) that did not localize to the ER, and (3) full-length endogenous BiP detected by antibody staining. We did not detect obvious reflux of endogenous BiP to the cytoplasm (Figure 9), with the caveat that in fixed samples, the BiP signal was not tightly co-localized with the ER marker even under control conditions. However, we did use this antibody to detect an increase in BiP in Atlastin mutant presynaptic terminals, indicating ER stress (see Reviewer 1.10).

      Though we have not identified endogenous targets, we believe that our studies with the exogenous reporter will be of great interest to the field, as they clarify the previously reported Atlastin phenotype and provide the first report of a new defect in a human disease animal model.

      In comparison to Summerville et al. (2016) in Figure 7, the experiment was not done in the same way. It is important to keep the same settings for comparison

      In Figure 7D-E, we compare the distribution of BiP:sfGFP:HDEL in cell bodies, axons, and muscles between controls and Atlastin mutants. To clarify the experimental approach relative to Summerville et al. (2016): while both our studies examined the same cellular compartments (cell bodies, axons and nerve terminals) using the BiP:sfGFP:HDEL reporter, we employed super-resolution Airyscan microscopy. This enhanced resolution was critical for definitively demonstrating that this is a functional rather than a structural phenotype and that ER displacement is progressive, and repeating this experiment at lower resolution as previously reported does not provide any new information. We identified two distinct distribution phenotypes in Atlastin mutants expressing BiP:sfGFP:HDEL, which were not described in the Summerville et al., 2016 paper. From our manuscript (lines 249-251): "We identified two distinct ER network phenotypes in Atlastin mutants expressing BiP:sfGFP:HDEL: "Partial loss" NMJs retained both diffuse signal and identifiable ER network structures, while "Complete loss" NMJs showed no visible ER network structures. Note that the "Complete loss" phenotype in Atlastin mutants reflects the absence of detectable luminal marker signal in organized ER structures, but not the complete absence of ER membranes, as demonstrated by our ER membrane marker tdTomato:Sec61β results."

      Does the Atlastin mutant induce the unfolded protein response and stress within the ER? It is necessary to look for UPR markers in those settings. It was shown previously that ER stress leads to protein reflux from the ER to the cytosol. Is there a difference in the ER stress markers in the presynaptic terminal?

      The reviewer suggested that Atlastin mutant synapses may exhibit ER stress. To address this, we examined levels of the ER chaperone BiP, a well-established ER stress marker whose expression increases during UPR activation. We first validated that our BiP antibody can detect changes in ER stress by feeding control larvae with 50mM DTT for 24 hours. These results are in the new Figure 10A. Note that we were unable to test sensitivity to ER stress in this way in Atlastin mutant larvae because they did not consume the DTT-treated food, as assessed by blue food coloring in the larvae's guts.

      Using this antibody, we measured baseline BiP levels at NMJs of Atlastin mutants on normal food, and found they were slightly increased compared to controls. We conclude from these experiments that Atlastin mutant synapses have mild ER stress. Notably however, Atlastin mutants co-expressing UAS-BiP:sfGFP:HDEL or UAS-tdTomato:Sec61b did not show significantly increased endogenous BiP levels, suggesting that transgene expression at least partly suppresses the mild ER stress response, even though there is extensive cytosolic displacement. These results argue (1) that the mild ER stress in Atl mutants does not strictly correlate with the reflux phenotype, and (2) that the reflux phenotype is not an artifact of overexpression-induced stress. These results are described on lines 430-436 in the results section and shown in Figure 10B-E, and their implications discussed on lines 585-598.

      We also explored another strategy to detect ER stress by assessing eIF2α phosphorylation, a key event in the Unfolded Protein Response (UPR) pathway. We obtained a phospho-eIF2α antibody (Cell Signaling; #3597) that was reported to work in Drosophila. However, when we tested this antibody by Western blot, we were unable to detect a band at the expected molecular weight for phosphorylated eIF2α, even in positive-control samples treated with DTT to induce ER stress. We therefore concluded that this antibody is not suitable for reliably detecting ER stress in our experimental system. The failure of this antibody highlights the challenges of finding robust tools to measure ER stress in Drosophila.

      It is important to add biochemical experiments to show that no fragmentation of the ER membrane occurred. It can be simply demonstrated by looking at the redox state of the ER, which would change if it were mixed with the reducing cytosol. Moreover, this can be shown by using an ER-targeted redox-sensitive fluorescent protein that is tethered to the ER membrane to follow changes in the redox state of the ER.

      The reviewer asked us to test whether the redox state of the ER is disrupted, which could indicate exchange between the cytosol and ER due to membrane rupture. As noted above, biochemical approaches such as fractionation are not possible in this in vivo sample. We attempted to address this concern by creating a UAS-Sec61β:roGFP construct, using the roGFP sequence from Igbaria et al. (2019) to monitor the ER lumen redox environment in Atlastin mutants. Since Sec61β is membrane-tethered, it should remain in the ER and not undergo reflux, making it an ideal sensor for detecting any mixing between the reducing cytosolic environment and the oxidizing ER lumen that would occur if membrane fragmentation and/or ruptures were present. We tested this approach in wild-type Drosophila S2 cells and used the Gal4-UAS binary expression system to co-express Actin-Gal4 (to drive expression of UAS constructs), UAS-Sec61β:roGFP (redox sensor), and UAS-BiP:Halo:HDEL (as a control reporter insensitive to DTT treatment).

      Our experiments showed no detectable changes in the fluorescent properties of UAS-Sec61β:roGFP following 30 min 10mM DTT treatment compared to DMSO vehicle control, including no increase in 405-nm excitation fluorescence or changes in 488nm/405nm excitation ratios. These results suggest that either the roGFP sensor requires further optimization for sensitivity in this cellular system or that additional controls and calibration steps are needed to establish the dynamic range of the assay. We believe this experiment falls beyond the scope of the current study, given the extensive optimization required. However, it represents an important future direction for testing membrane fragmentation as a mechanism underlying the phenotypes observed in Atlastin mutants. The possibility of ER integrity defects is mentioned in the discussion on lines 547-559.

      Minor comments:

      1. It is important to call figures by order. Figure 2C is called before 2A-B. Figure 2B is called before Figure 2A.

      The revised manuscript has all figures in order of appearance in the text.

      Figure legends (Figure 2): "The same control dataset used in E-G was used in Figure 5 and Figure 5_Supplement." Why is this relevant?

      We wanted to be transparent about reusing the same control dataset across multiple figures to avoid any appearance of data duplication. This notation clarifies that, although the data appear in different contexts (Figures 2 and 5. This version does not contain a Figure 5_Supplement), it represents the same biological samples analyzed for different parameters, ensuring readers understand that these are not independent datasets.

      Figure 4F is called before Figure-4D-E which are not called.

      We revised our manuscript and reorganized Figure 4 to ensure that all figure panels are referenced in sequential order and that panels 4D-E, which were previously not cited in the text, are now properly referenced when discussing their corresponding results.

      Figure 5B is called before the previous ones. Same for Figure 5A supplement.

      We referenced Figure 5A in lines 211-212, which precedes our discussion of Figure 5B. To clarify the figure order, we removed the early references to Figures 2D-G and Movies 7-14, which were mentioned only to indicate that we were analyzing the same dataset in different ways.

      The revised manuscript has all figures in order of appearance in the text.

      Referees cross-commenting

      I agree with the comments raised by reviewer2 and 3. Basically it is highly important to validate those data by genetic rescue. Moreover, it is essential to know the source of the displaced luminal marker to the cytosol. Is it mislocalization or it is a reflux of pre-existing protein to the cytosol after insertion to the ER. It is also recommended by me and the reviewers and me to test the endogenous protein rather than overexpression.

      We have addressed these points in our responses to the following reviewer questions:

      • Genetic rescue: Please see our responses to Reviewer 1/Question #10 and Reviewer 2/Question #1.
      • Source of displaced luminal marker: We provide some evidence addressing this in our response to Reviewer 3/Question #1.
      • Endogenous protein localization: We have examined this and detailed our findings in our responses to Reviewer 1/Question #7 and Reviewer 2/Question #6.

        Reviewer #1 (Significance (Required)):

      General assessment: This interesting paper shows that proteins can escape the ER under special conditions. However, the authors need more evidence to show that and rely less on the overexpression system, especially of BIP-GFP, which can cause proteostasis stress within the ER. Advance: The results have been oversimplified in their explanations, and some points and complexities of the study need to be addressed further to make the most of them. These are often some of the more interesting concepts in the paper. I think many points can be addressed in the text by the authors being clear and concise with their reporting. At the same time, other experiments would turn this paper from an observational one into a very interesting mechanistic one. This paper is based on previously published articles from the group and other groups, and it is a nice progression. However, as mentioned, this paper depends primarily on published data, and the novelty is somehow lost between all the comparisons to other published data instead of emphasizing that. Without a substantial mechanistic improvement, the paper would remain observatory.

      Audience: The microscopy tools can be great addition to researchers in the field to monitor protein trafficking especially Cell biologists (basic research)

      My expertise: ER homeostasis, protein trafficking, cell biology

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

      Summary The endoplasmic reticulum (ER) is a continuous organelle that extends throughout neurons to regulate fundamental processes. The analysis of ER dynamics at synaptic terminals is limited by the challenge of imaging these structures at high resolution. In this manuscript, the authors use super-resolution (~170 nm) live imaging and a combination of membrane and luminal ER markers at the Drosophila larval NMJ, an important model synapse, to investigate dynamic ER architecture in vivo. They report a detailed characterization of the presynaptic ER organization and dynamics at wild-type and GTPase Atlastin mutant NMJs. Their analysis using the ER membrane marker tdTomato:Sec61b reveals the presence of an intact ER network in Atlastin mutants. This contrasts with the apparent ER fragmentation phenotype previously reported and replicated here when using a luminal marker. Their findings instead point to the progressive displacement of luminal proteins to the cytosol in Atlastin mutants specifically at synapses. The authors propose that the disruption of ER protein dynamics at synapses is a compartment-specific ER stress response. The manuscript is well written, results are clearly presented, and experiments are technically rigorous.

      Major comments

      1. The baseline ER phenotypes in Atlastin mutants are mild with complete loss of ER network only observed in terminal boutons. This interesting and unexpected result should be further confirmed by genetic rescue. The authors can use a UAS rescue line previously reported in PMID: 19341724.

      We tested the UAS-Atl-myc rescue line and unfortunately found that even in wild-type neurons, overexpression of Atlastin produced strong ER organization defects that precluded the rescue experiment. Instead, to confirm the cell autonomy of the phenotype and to test it wth an independent tool, we performed a presynaptic knockdown of Atlastin by RNAi and found that BiP:sfGFP:HDEL is displaced, as observed in the Atlastin null mutant. These results are in now shown in Figure 7A-C.

      Lines 204-7: It's not clear how a greater coefficient of variation indicates that the marker is more concentrated in subsynaptic structures or what is meant by 'subsynaptic structures.'

      We added the following text to explain, in lines 181-183: "A higher CoV indicates an uneven distribution of tdTomato:Sec61β within the presynaptic terminal, with some areas showing higher concentrations than others (in contrast to the uniform, diffuse signal expected from fragmentation)." To avoid confusion with postsynaptic structures called the subsynaptic reticulum, we have removed the term "subsynaptic". The intended meaning is distinct structures found within the presynaptic terminal.

      There's a mistake in Figure 6C and the associated text. The summed percentage of the three phenotypic categories adds up to 110% for Atlastin mutants. *

      The reviewer noted that the summed percentage of the three phenotypic categories in Figure 6C adds up to 110% for Atlastin mutants, which appears to be a mathematical error. However, this is not an error, but rather a reflection of our quantification methodology, in which a single bouton can exhibit more than one type of ER dynamics per movie recorded. Our quantification counts each phenotype independently, so boutons displaying multiple phenotypes contribute to more than one category. This approach provides a more comprehensive view of the range of ER dynamics present in Atlastin mutants, as restricting the analysis to mutually exclusive categories would underrepresent the complexity of the phenotypes observed. To make this point clear, we made the following change to the text in lines 257-259: "We note that the sum of these percentages exceeds 100% because one NMJ exhibited multiple phenotypes: one branch had a complete loss, while the other branch had no phenotype. These phenotypes were counted separately."

      Figure 8: the ER looks fragmented in 1st instar controls and mutants. The authors should address this difference from more mature NMJs.

      We would like to clarify that the bulk of experiments in this manuscript (including all ER dynamics, luminal marker redistribution, and membrane marker analyses discussed throughout the Results) were performed in 3rd instar larvae, which are more mature larval NMJ preparations standard in the field. Figure 8 was included specifically to test whether the Atlastin mutant phenotype we describe throughout the paper is also detectable at an earlier developmental stage, not to replace or reinterpret our primary findings.

      Regarding the specific observation that the ER appears more fragmented in Figure 7F-H relative to the more mature NMJs shown elsewhere: this fragmentation, observed similarly in both control and Atlastin mutant 1st instar larvae, likely reflects technical challenges associated with dissecting these smaller, more delicate early-stage specimens rather than a genotype-specific effect. Because fragmentation occurred similarly in both genotypes, we could still reliably assess the redistribution of BiP:sfGFP:HDEL as our primary phenotypic readout in this experiment. We have added the following text (lines 306-309) to clarify this point: "Note that in 1st instar larvae, both normal networks in controls and residual networks in Atlastin mutants appeared more fragmented than in 3rd instar preparations, likely due to the technical challenges of dissecting these smaller, more delicate specimens. Since ER fragmentation occurred similarly in both genotypes, we could still reliably assess the redistribution of BiP:sfGFP:HDEL as our primary phenotypic readout.

      The images in figure 9B do not seem representative of the quantification in Figure 9D. Specifically, the partial loss Atlastin NMJ appears to have recovered as fully as the complete loss Atlastin NMJ.

      The images showed FRAP recovery across the entire bouton, but we photobleached only a small region within each bouton and quantified only this region. We have now added outlines to clearly delineate the specific FRAP regions that were analyzed in each image, which clarify that the partial loss Atlastin showed less recovery than the overall bouton. We have also reordered the figures to more clearly convey our message (Figure 9 is now Figure 8).

      We also made a few changes to the paragraph on lines 347-350 to clarify our experimental reasoning: "We photobleached en passant boutons using a defined region of 6.8 x 7.8 microns (dashed box in Figure 8D) to ensure that BiP:sfGFP:HDEL could recover from the ER networks surrounding the FRAP region (Movies S20-S23)."

      We also added this sentence to the figure legends of Figure 8: "The dashed boxes in (D) indicate areas that were photobleached and analyzed for recovery quantification in (E-F)."

      Optional: An overexpressed luminal marker is displaced to the cytoplasm in Atlastin mutants. It would be interesting to know and increase the significance of the findings if the same is true of endogenous luminal proteins under biological stress conditions.

      As noted in our response to Reviewer #1 suggested that Atlastin mutant synapses may exhibit ER stress. To address this, we examined levels of the ER chaperone BiP, a well-established ER stress marker whose expression increases during UPR activation. We first validated that our BiP antibody can detect changes in ER stress by feeding control larvae with 50mM DTT for 24 hours. We were unable to perform this experiment in Atlastin mutant larvae because they did not consume the DTT-treated food, as assessed by blue food coloring in the larvae's guts. These results are in Figure 10A. In the future, it will be of interest to establish a protocol to examine Atlastin mutants by feeding or treating larval fillets with DTT.

      We measured BiP levels at NMJs of Atlastin mutants and found they were slightly increased compared to controls. Atlastin mutants co-expressing UAS-BiP:sfGFP:HDEL or UAS-tdTomato:Sec61b did not show significantly increased endogenous BiP levels, suggesting that transgene expression suppresses the mild ER stress response. We conclude from these experiments that Atlastin mutant synapses have mild ER stress. These results are in Figure 10B-E).

      Optional: Applying this approach in stimulated conditions (high potassium, increased temperature) might reveal a greater activity-dependent role for Atlastin at synaptic terminals.

      This is a very interesting idea, as we have only examined synapses at rest. However, this is beyond the scope of this paper.

      Minor Comments

      1. Line 16: Atlastin should be italicized.

      Thank you for catching this typo. We have fixed it.

      Figure 5A: Based on the relative intensities, it appears that control and mutant images are not contrast matched but this isn't stated.

      Thank you for catching this omission. We added to the figure legend: "Control and Atlastin mutant images are not contrast matched."

      Line 822: The number of static Atlastin mutant boutons used for analysis is missing.

      Thank you for catching this omission. We have fixed this supplementary table.

      Figure 9: The blue arrows are not annotated in the figure legend.

      Thank you for catching this omission. We have fixed this figure legend.

      Reviewer #2 (Significance (Required)):

      Atlastin is linked to Hereditary Spastic Paraplegia (HSP) and this study changes our understanding of the compartment-specific impacts of its loss. This study reveals the importance of using both membrane and luminal ER markers to accurately interpret phenotypes as well as the importance of considering compartment-specific effects on ER. These findings represent significant mechanistic and conceptual advances. The lack of genetic rescue is a limitation and adding an investigation of an endogenous luminal protein under basal and stress conditions would add significantly to our understanding of Atlastin dysfunction in HSP. Notably, the in vivo imaging approach introduced here can be adapted broadly for live imaging of Drosophila larvae. Thus, this work will be of interest to both neuronal cell biologists and the wider Drosophila community. This review is based on our expertise in neuronal cell biology.

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

      In this manuscript, the authors investigate the structural dynamics of the endoplasmic reticulum (ER) in Drosophila neurons and examine the role of the ER-shaping protein Atlastin in ER morphology. Their discovery on the neuromuscular junction (NMJ)-specific contribution of Atlastin to ER integrity is intriguing and may provide valuable insights into the pathological mechanisms underlying Atlastin mutations associated with hereditary spastic paraplegia (HSP) and hereditary sensory neuropathy. The key observation on ER protein showing an aberrant cytoplasmic localisation in mutant cells appears convincing. Though this phenomenon's characterisation stays at the point of primary observation with its mechanics unclarified, establishing this new and unexpected functional rather than structural Atl effect is important and useful for the field. The observation that ER is structurally preserved in this mutant with absolute lack of Atl are also extremely useful.

      It is unclear if the cytoplasmic localisation affects an exogenous overexpressed ER marker or endogenous protein would also appear in cytoplams, the authors should consider adding an immunostaining data to test that.

      Authors offer speculations on potential reasons for the cyto localisation of the ER marker suggesting that relocation at the cell periphery specifically combined with slow clearance there is the most likely explanation (still unclear what stops the marker from spreading through the entire cell). They suggest that decrease in cotranslational translocation is unlikely as this would result in somatic accumulation of the marker. However, if the clearance in the periphery is less efficient than in soma, the accumulation there might reflect a compromised translocation. Any clarifying experiments, if practical, to directly demonstrate how ER proteins in relocates to the cytoplasm in atl mutant would help understanding better the phenomenon. For example, would proteasomal inhibition make the marker accumulate more across the cell? Authors also suggest links to ER stress. Would stress induction phenocopy the mutant?

      Reviewer #3 asked whether defective proteasomal clearance underlies the cytosolic accumulation of BiP:sfGFP:HDEL in Atlastin mutants. We addressed this directly. First, proteasome function appears intact in the mutants: baseline ubiquitinated protein levels (FK1 antibody) were comparable between control and Atlastin mutants, and MG132 treatment produced a similar increase in ubiquitination in both genotypes, confirming both antibody specificity and normal proteasome activity. We then examined BiP:sfGFP:HDEL directly. In controls, MG132 caused the marker to accumulate at axons and presynaptic terminals, showing that it is normally cleared from these compartments by the proteasome. Critically, this accumulated marker remained associated with intact ER networks: MG132 did not induce diffuse cytosolic BiP:sfGFP:HDEL in any compartment (cell bodies, axons, or presynaptic terminals), even where levels rose substantially. Thus, blocking proteasomal clearance raises ER-localized marker but does not generate the cytosolic pool seen in Atlastin mutants, indicating that impaired clearance is not sufficient to cause the displacement phenotype. We separately noted that BiP:sfGFP:HDEL was already elevated in Atlastin mutant axons without MG132, paralleling the axonal tdTomato:Sec61β accumulation in Figure 4, consistent with reduced baseline clearance specifically in mutant axons, but this does not lead to cytosolic displacement. This experiment is now shown in Figure 11, described in Results (lines 445-475), and discussed in lines 576-581.

      Minor comments:

      Line 146:

      "fast dynamics (Thank you for catching this mistake. We have corrected it.

      Fig. 2D: The data representation of "Tubule displacement" image is unclear. The ER tubule indicated by the red arrow does not seem to show any changes over time (like static). time 0 in stamp appears behind the image.

      Thank you for catching the typo. We have fixed it. Additionally, we added black arrows to highlight a tubule that is not moving, allowing the reader to compare it with the moving tubule. We also included a video of all types of ER tubule dynamics to ensure the reader can also look at the raw data (Movies S9-11).

        • Line 157-158 (and relevant method sections):

      The definition of static and dynamic boutons is ambiguous. The author should describe in more detail this point including how long they observed the structure to define the changes in ER tubule dynamics.

      We provide in the methods (lines 779-791) a detailed explanation of how we categorized boutons as dynamic or static. In addition, we added the following to explain in the results section how we defined static vs dynamic:

      Old sentence: We qualitatively categorized boutons as "static" if we observed no change in ER network structure or "dynamic" if we observed at least one change.

      New sentence in lines 143-147: "We imaged boutons for 40 sec at 0.92 sec intervals to capture ER dynamics over this observation period. Boutons were qualitatively categorized as "static" if we observed no detectable changes in ER network structure throughout the entire 40 sec imaging session, or "dynamic" if we observed at least one of the three defined dynamic events during this time window."

      Fig. 2E: What n=75 and n=29 represent is unclear, are these the number of boutons in en passant and terminal subjected for qualitative analysis?

      We removed these n values from the figure and added this information to the Supplementary Table 1, which contains detailed information about the genotype, statistical analysis, and number of larvae and NMJs analyzed.

      Fig. 2: What the qualitative analysis represents is unclear, are the points pulled from different experiments?

      The data in Fig. 2 E-F comes from movies acquired in the same experiment. The number of independent animals and NMJs imaged is described in Table 1.

      * *Line 231: Regarding "...we found a small but significant reduction in dynamic boutons in Atlastin mutants (76%), ...", how do the authors assess significance. If proportion of static/dynamic ER in boutons was obtained from multiple experiments, it should be presented e.g. as in average {plus minus} standard deviation, or clarify that the proportion is representative of x independent experiments.

      The videos used for this figure were acquired from a single experiment. We use a chi-square test to determine significance relative to the "expected" distribution of dynamics types from controls, as these are categorical rather than continuous data (see PMID 31145670). Information regarding genotype, statistical analysis and number of larvae and NMJs can also be found in Supplementary Table 1.

      Line 267-269 and Fig. 6B: The author's conclusion that "Complete loss of ER network structure in NMJ of BiP:sfGFP:HDEL overexpressing Atl mutant" seem to be based on the lack of signal from luminal marker, which may be undetectable due to changes to tubular volume or marker loss to the cytoplasm, as suggested by the authors, while the membranous ER structure is intact. It would be useful to discuss this point and potentially add ER membrane-stained control.

      We agree with the reviewer that Atlastin mutants categorized as 'complete loss mutants' do not actually lack ER at synapses. We think this is an important point so we added the following to the results in lines 251-254: "Note that the "Complete loss" phenotype in Atlastin mutants reflects the absence of detectable luminal marker signal in organized ER structures, not the complete absence of ER membranes, as demonstrated by our ER membrane marker tdTomato:Sec61β results."

      We attempted to co-label the ER membrane and ER lumen, but these crosses yielded very few live larvae (in either controls or Atlastin mutants, and those that survived had severely deformed NMJs. We added Figure 6-Supplement showing the results of this experiment, and described them on lines 270-273.

      Fig. 6C: In Atl mutant, why does the total of the proportion exceed 100% (10 + 45 + 55)?

      The reviewer noted that the summed percentage of the three phenotypic categories in Figure 6C adds up to 110% for Atlastin mutants. This is not an error, but rather a reflection of our quantification methodology because a single bouton can exhibit more than one type of ER dynamics per movie recorded. Our quantification counts each phenotype independently, so boutons displaying multiple phenotypes contribute to more than one category. This approach provides a more comprehensive view of the range of ER dynamics present in Atlastin mutants, as restricting the analysis to mutually exclusive categories would underrepresent the complexity of the phenotypes observed. To make this point clear, we made the following change to the text in lines 257-259: "We note that the sum of these percentages exceeds 100% because one NMJ exhibited multiple phenotypes: one branch had a complete loss, while the other branch had no phenotype. These phenotypes were counted separately."

      Fig. 9C, line 342-344: In FRAP experiment using CD8, it seems that the Partial loss Atl mutant shows slower recovery that control. There seems to be a mismatch in triangle symbols of Partial loss Atl mutant between legend and plot (one is filled and the other is empty). This should be clarified.

      Thank you for catching this mistake. We have fixed the figure.

      fig. 10 is a clever way to verify the cytoplasmic localization of the ER marker; however, its description and annotation can be improved, and it would be stronger if 4 curves in F for mutant and controls with the trap and normal were shown.

      The reviewer suggested merging our graphs but we believe that keeping them separate is clearer.

      Line 495: Drosophila have ReepA and ReepB, but not Reep1-4. If the authors discuss their speculation based on their observation (using Drosophila), the gene names should be unified in the same species, and explain the corresponding genes to mammalian cells.

      We made the following changes to address the reviewer's concern about gene nomenclature consistency (lines 502-506): "These ER-derived vesicles are likely to involve ReepA and ReepB, the Drosophila orthologs of mammalian REEP1-4, which regulate ER vesicle formation in mammalian cells (67). Notably, while overexpression of Atlastin can regulate REEP vesicle fusion in mammalian systems (67), it is not essential for vesicle formation, suggesting similar regulatory relationships may exist between Atlastin and Reep genes in Drosophila."

      Line 548; should UPR be Unfolded Protein Response?

      Thank you for catching the typo. We have fixed it.

      Reviewer #3 (Significance (Required)):

      This study advances the understanding of how ER morphogens affect neuronal cells specifically, the lack of which limits researchers ability to comprehend the neuronal pathologies associated with ER structure-function. The observation on ER content aberrant localisation caused by the lack of key structural protein should be of a great interest for cell and neuronal biologists and researchers of the associated diseases and shows the field a new direction. Though, mechanistic details remain to be unraveled, it constitutes a fundamental, conceptual advance.

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

      Learn more at Review Commons


      Referee #3

      Evidence, reproducibility and clarity

      In this manuscript, the authors investigate the structural dynamics of the endoplasmic reticulum (ER) in Drosophila neurons and examine the role of the ER-shaping protein Atlastin in ER morphology. Their discovery on the neuromuscular junction (NMJ)-specific contribution of Atlastin to ER integrity is intriguing and may provide valuable insights into the pathological mechanisms underlying Atlastin mutations associated with hereditary spastic paraplegia (HSP) and hereditary sensory neuropathy. The key observation on ER protein showing an aberrant cytoplasmic localisation in mutant cells appears convincing. Though this phenomenon's characterisation stays at the point of primary observation with its mechanics unclarified, establishing this new and unexpected functional rather than structural Atl effect is important and useful for the field. The observation that ER is structurally preserved in this mutant with absolute lack of Atl are also extremely useful.

      It is unclear if the cytoplasmic localisation affects an exogenous overexpressed ER marker or endogenous protein would also appear in cytoplams, the authors should consider adding an immunostaining data to test that.

      Authors offer speculations on potential reasons for the cyto localisation of the ER marker suggesting that relocation at the cell periphery specifically combined with slow clearance there is the most likely explanation (still unclear what stops the marker from spreading through the entire cell). They suggest that decrease in cotranslational translocation is unlikely as this would result in somatic accumulation of the marker. However, if the clearance in the periphery is less efficient than in soma, the accumulation there might reflect a compromised translocation. Any clarifying experiments, if practical, to directly demonstrate how ER proteins in relocates to the cytoplasm in atl mutant would help understanding better the phenomenon. For example, would proteasomal inhibition make the marker accumulate more across the cell? Authors also suggest links to ER stress. Would stress induction phenocopy the mutant?

      Minor comments:

      Line 146: "fast dynamics (<1 sec)" a velocity should be presented as distance/time

      Fig. 2D: The data representation of "Tubule displacement" image is unclear. The ER tubule indicated by the red arrow does not seem to show any changes over time (like static). time 0 in stamp appears behind the image.

      Line 157-158 (and relevant method sections): The definition of static and dynamic boutons is ambiguous. The author should describe in more detail this point including how long they observed the structure to define the changes in ER tubule dynamics.

      Fig. 2E: What n=75 and n=29 represent is unclear, are these the number of boutons in en passant and terminal subjected for qualitative analysis?

      Fig. 2: What the qualitative analysis represents is unclear, are the points pulled from different experiments?

      Line 231: Regarding "...we found a small but significant reduction in dynamic boutons in Atlastin mutants (76%), ...", how do the authors assess significance. If proportion of static/dynamic ER in boutons was obtained from multiple experiments, it should be presented e.g. as in average {plus minus} standard deviation, or clarify that the proportion is representative of x independent experiments.

      Line 267-269 and Fig. 6B: The author's conclusion that "Complete loss of ER network structure in NMJ of BiP:sfGFP:HDEL overexpressing Atl mutant" seem to be based on the lack of signal from luminal marker, which may be undetectable due to changes to tubular volume or marker loss to the cytoplasm, as suggested by the authors, while the membranous ER structure is intact. It would be useful to discuss this point and potentially add ER membrane-stained control.

      Fig. 6C: In Atl mutant, why does the total of the proportion exceed 100% (10 + 45 + 55)?

      Fig. 9C, line 342-344: In FRAP experiment using CD8, it seems that the Partial loss Atl mutant shows slower recovery that control. There seems to be a mismatch in triangle symbols of Partial loss Atl mutant between legend and plot (one is filled and the other is empty). This should be clarified

      fig. 10 is a clever way to verify the cytoplasmic localisatoin of the ER marker, however its description and annotation can be improved, and it would be stronger if 4 curves in F for mutant and controls with the trap and normal were shown.

      Line 495: Drosophila have ReepA and ReepB, but not Reep1-4. If the authors discuss their speculation based on their observation (using Drosophila), the gene names should be unified in the same species, and explain the corresponding genes to mammalian cells.

      Line 548; should UPR be Unfolded Protein Response?

      Significance

      This study advances the understanding of how ER morphogens affect neuronal cells specifically, the lack of which limits researchers ability to comprehend the neuronal pathologies associated with ER structure-function. The observation on ER content aberrant localisation caused by the lack of key structural protein should be of a great interest for cell and neuronal biologists and researchers of the associated diseases and shows the field a new direction. Though, mechanistic details remain to be unraveled, it constitutes a fundamental, conceptual advance.

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

      Learn more at Review Commons


      Referee #2

      Evidence, reproducibility and clarity

      Summary

      The endoplasmic reticulum (ER) is a continuous organelle that extends throughout neurons to regulate fundamental processes. The analysis of ER dynamics at synaptic terminals is limited by the challenge of imaging these structures at high resolution. In this manuscript, the authors use super-resolution (~170 nm) live imaging and a combination of membrane and luminal ER markers at the Drosophila larval NMJ, an important model synapse, to investigate dynamic ER architecture in vivo. They report a detailed characterization of the presynaptic ER organization and dynamics at wild-type and GTPase Atlastin mutant NMJs. Their analysis using the ER membrane marker tdTomato:Sec61b reveals the presence of an intact ER network in Atlastin mutants. This contrasts with the apparent ER fragmentation phenotype previously reported and replicated here when using a luminal marker. Their findings instead point to the progressive displacement of luminal proteins to the cytosol in Atlastin mutants specifically at synapses. The authors propose that the disruption of ER protein dynamics at synapses is a compartment-specific ER stress response. The manuscript is well written, results are clearly presented, and experiments are technically rigorous.

      Major comments

      1. The baseline ER phenotypes in Atlastin mutants are mild with complete loss of ER network only observed in terminal boutons. This interesting and unexpected result should be further confirmed by genetic rescue. The authors can use a UAS rescue line previously reported in PMID: 19341724.
      2. Lines 204-7: It's not clear how a greater coefficient of variation indicates that the marker is more concentrated in subsynaptic structures or what is meant by 'subsynaptic structures.'
      3. There's a mistake in Figure 6C and the associated text. The summed percentage of the three phenotypic categories adds up to 110% for Atlastin mutants.
      4. Figure 8: the ER looks fragmented in 1st instar controls and mutants. The authors should address this difference from more mature NMJs.
      5. The images in figure 9B do not seem representative of the quantification in Figure 9D. Specifically, the partial loss Atlastin NMJ appears to have recovered as fully as the complete loss Atlastin NMJ.
      6. Optional: An overexpressed luminal marker is displaced to the cytoplasm in Atlastin mutants. It would be interesting to know and increase the significance of the findings if the same is true of endogenous luminal proteins under biological stress conditions.
      7. Optional: Applying this approach in stimulated conditions (high potassium, increased temperature) might reveal a greater activity-dependent role for Atlastin at synaptic terminals.

      Minor Comments

      1. Line 16: Atlastin should be italicized.
      2. Figure 5A: Based on the relative intensities, it appears that control and mutant images are not contrast matched but this isn't stated.
      3. Line 822: The number of static Atlastin mutant boutons used for analysis is missing.
      4. Figure 9: The blue arrows are not annotated in the figure legend.

      Significance

      Atlastin is linked to Hereditary Spastic Paraplegia (HSP) and this study changes our understanding of the compartment-specific impacts of its loss. This study reveals the importance of using both membrane and luminal ER markers to accurately interpret phenotypes as well as the importance of considering compartment-specific effects on ER. These findings represent significant mechanistic and conceptual advances. The lack of genetic rescue is a limitation and adding an investigation of an endogenous luminal protein under basal and stress conditions would add significantly to our understanding of Atlastin dysfunction in HSP. Notably, the in vivo imaging approach introduced here can be adapted broadly for live imaging of Drosophila larvae. Thus, this work will be of interest to both neuronal cell biologists and the wider Drosophila community. This review is based on our expertise in neuronal cell biology.

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

      Learn more at Review Commons


      Referee #1

      Evidence, reproducibility and clarity

      In the present manuscript, the authors address an important question related to the ultrastructure and the dynamics of the ER in HSP. In contrast to previous studies, the authors here show (by using a membrane and luminal protein markers) that in the presynaptic terminals, the overexpressed BIP "mislocalizes" to the cytosol without affecting (or with minimal effect) the integrity of the ER membrane. Although they used an artificial system by overexpressing (overexpression) of BIP-sfGFP-HDEL (fused protein), the findings lack validation of the endogenous protein by biochemical and fluorescent tools.

      Concerns:

      I am worried about how the article is presented, mainly in the results section, as most of it refers to published data. The Results section is reserved for presenting new findings without external interpretation or comparison. The paper is written mainly as a comparison paper with other studies or relies on previous studies to strengthen their findings rather than coming up with novel findings. Up to figure 4, I missed the relevance of the new findings. The manuscript needs rewriting to emphasize its novelty and significance without comparing it to previous data. Moreover, the manuscript emphasizes the technology and the findings of the localization of the ER luminal proteins to the cytosol (which is not novel and was previously reported in other settings). Those two aspects were not given enough focus. Here are the main major and minor comments:

      Major comments:

      1. Lines 103-116 (first paragraph of the results section) describe mainly published data that is more suitable for the introduction section. It is annoying to refer to different published articles in the Results section to strengthen the results instead of showing them. The same goes for paragraphs two and three. Why mention those data in the Results section if they are already published and known?
      2. Figure Legends-(in all Figures): The number of experimental repeats must be mentioned in the figure legends.
      3. The way the figures are labeled is worrisome; supplementary figures are not ordered numerically.
      4. The tubule extension in Figure 2D is not convincing. Is there a movie showing those changes? Better images are needed. It is essential to show which supplemental movie corresponds to which panel.
      5. This is unnecessary in the results section: "To investigate the relationship between ER structure and function at synapses, we examined mutants of Atlastin, a GTPase that regulates ER tubule fusion. Drosophila has a single homolog while mammals have three Atlastin homologs, with Atlastin-1 enriched in the brain (Rismanchi et al., 2008)."
      6. "This reduction in ER membrane marker intensity has also been observed in other HSP mutants, suggesting this is a common feature of ER shaping mutants and could indicate changes in ER membrane composition, integrity, or tubule thickness (P.rez-Moreno et al., 2023)." This comparison is important and should be shown in the same settings as for the Atlastin mutant rather than referring to published data.
      7. Does the distribution of the luminal ER marker in Figure 6F diffuse due to mislocalization or reflux after being localized to the ER and then refluxed to the cytosol as was previously shown for the ER to Cytosol signaling (ERCYS) mechanism? Could you assess other ER-luminal protein localization biochemically? It is highly recommended to look at another soluble ER-protein localization in the Atlastin mutant without overexpression, which can be an artifact
      8. "(data not shown)" in line 288. This affects the process of judging those data.
      9. In comparison to Summerville et al. (2016) in Figure 7, the experiment was not done in the same way. It is important to keep the same settings for comparison
      10. Does the Atlastin mutant induce the unfolded protein response and stress within the ER? It is necessary to look for UPR markers in those settings. It was shown previously that ER stress leads to protein reflux from the ER to the cytosol. Is there a difference in the ER stress markers in the presynaptic terminal?
      11. It is important to add biochemical experiments to show that no fragmentation of the ER membrane occurred. It can be simply demonstrated by looking at the redox state of the ER, which would change if it were mixed with the reducing cytosol. Moreover, this can be shown by using an ER-targeted redox-sensitive fluorescent protein that is tethered to the ER membrane to follow changes in the redox state of the ER.

      Minor comments:

      1. It is important to call figures by order. Figure 2C is called before 2A-B. Figure 2B is called before Figure 2A.
      2. Figure legends (Figure 2): "The same control dataset used in E-G was used in Figure 5 and Figure 5_Supplement." Why is this relevant?
      3. Figure 4F is called before Figure-4D-E which are not called.
      4. Figure 5B is called before the previous ones. Same for Figure 5A supplement.

      Referees cross-commenting

      I agree with the comments raised by reviewer2 and 3. Basically it is highly important to validate those data by genetic rescue. Moreover, it is essential to know the source of the displaced luminal marker to the cytosol. Is it mislocalization or it is a reflux of pre-existing protein to the cytosol after insertion to the ER. It is also recommended by me and the reviewers to test the endogenous protein rather than overexpression.

      Significance

      General assessment: This interesting paper shows that proteins can escape the ER under special conditions. However, the authors need more evidence to show that and rely less on the overexpression system, especially of BIP-GFP, which can cause proteostasis stress within the ER.

      Advance: The results have been oversimplified in their explanations, and some points and complexities of the study need to be addressed further to make the most of them. These are often some of the more interesting concepts in the paper. I think many points can be addressed in the text by the authors being clear and concise with their reporting. At the same time, other experiments would turn this paper from an observational one into a very interesting mechanistic one. This paper is based on previously published articles from the group and other groups, and it is a nice progression. However, as mentioned, this paper depends primarily on published data, and the novelty is somehow lost between all the comparisons to other published data instead of emphasizing that. Without a substantial mechanistic improvement, the paper would remain observatory.

      Audience: The microscopy tools can be great addition to researchers in the field to monitor protein trafficking especially Cell biologists (basic research)

      My expertise: ER homeostasis, protein trafficking, cell biology

    1. Research has found that leaders with a high task orientation are likely to emerge in both highly structured contexts like a group that works to maintain a completely automated factory unit and highly unstructured contexts like a group that is responding to a crisis.

      I always found it very interesting to hear from or about someone who stepped up in a crisis or high stress scenario. Especially if that person is not regularly in a leadership role. In a high-pressure group setting finding the nerve to not only take initiative but also direct others is impressive to me. Its also interesting that the others who didn't step up are usually willing and receptive to the directions of someone else even if they are a stranger. It seems the higher the stakes the quicker everyone falls in line and moves once someone has stood up as a leader even if it's a self-appointed role.

    1. Figure 3:

      For all figures:

      I think the font size for the axis labels adn legends are too small.

      I also wonder if it would be more intuitive to show plots of proportions over time rather than cumulative coounts. So on top you'd have proportion uninfects vs infected with susceptible strain vs resistant strain. On the bottom you'd have proportion in the NI-T, NI-NT, I-T, I-NT categories. (I know that wouldn't let you visualize the inter-quartile spread though .. hmm).

      The long and short is that I think you want to decide what the key message is that you want readers to understand from each figure and then plot the results in a way that is as easy as possible to see that result. Also state that key result in the legend. As it stands, I'm not sure that cummulative doses is germane to your key message(s).

    2. First, we increase patient volume to 10 patients per timestep.

      Explain why this is this more realistic - is the idea that there is no update to environmental conditions until 10 patients are simulated?

    3. This creates both between-subpopulation differences (high-risk vs. low-risk mean parameter values) and within-subpopulation variation (σ=0.1 for all attributes)

      It's unclear what aspects of these features are / aren't observed

    4. mixed low-risk and high-risk heterogeneous populations as ground truth

      Not sure what this means. I thought the ground truth is simply which patients are vs aren't infected?

    5. Accurate: true risk levels correctly observed,

      For all three scenarios: I'd be more explicit about what this means. I.e. does it mean the agents know if people are high vs low risk, and that high risk ahs mu of 0.84 while low risk has mu of 0.56?

    6. optimal policies should exploit.

      Can agents tell if patients are in the high vs low risk groups? If so - I'd remind readers of that. If not, then how else can the ground-truth stratification be exploited?

    7. where infection probability followed Gaussian distributions

      If a Gaussian distributio nis used, it seems that the probability of infection could fall outside of [0,1]?

    8. patient information and visible AMR data, but do not update behavior through experience

      Maybe 'patient information and variably-updated clinical practice guidelines, but do not update behavior based on the most current data on clinical performance.'

    1. If some kids in your neighborhood get hold of a little hatchet and chop off thin slices of bark all the way around the base of one of your trees, what will happen to the tree? Why?

      it wil probably dies, because its cut of the xylem and phloem supply connections

    2. What is the most exterior cell layer in an herbaceous stem called? What is the most exterior layer of cells in a five-year-old woody perennial plant stem called.

      The most exterior layer of cells in a five-year-old woody perennial plant stem is called the cork (or phellem)

    3. Describe the origin of annual rings. If a woody dicot is growing in a tropical climate where the weather is the same day in and day out, will you find annual rings in the wood?

      Annual rings originate from the seasonal variation in the activity of the vascular cambium, which alters the size and density of the secondary xylem cells it produces throughout the yea

    4. In a perennial woody dicot, how do the discrete vascular bundles found in the new seedling stem become continuous rings of xylem and phloem in the three-year-old woody stem?

      Trace the boundary between the first, second, and third-year woodIdentify the vascular rays that transport water radiallyLocate the crushed remnants of the primary phloem

    5. If shown a micrograph of an apical meristem, how would you determine whether it is from a root or a shoot?

      look for the external structure and tissue arrangement root apical has a root cap and shoot apiece meristem has the leaf permordia

    1. ARROSTITO: Toda la «organización» éramos doce personas, entre los de Buenos Aires y los de Córdoba. En el operativo jugamos diez.

      De ahi saca Gillespie el nùmero.

    1. the burden ofdiscourse assimilation invariably falls on African American students

      The attempt to assimilate African American students into using only standard English?

    1. When adults validate and celebrate children’s diverse language abilities, they create more productive and engaging learning environments. These abilities include language mixing and innovation.In short, flexible multilingualism contributes to the vitality of diverse mother languages and brings tangible benefits to language learners.

      This is the main idea of the article. Multilingualism should be seen as a strength rather than a weakness. I can use this in my research about language and identity.

    2. My own research with Tibetan families living in urban centres shows that parents encourage children to speak a standard language rather than their regional mother language. Despite the significance of regionally diverse mother languages to adults’ identities, Tibetan communities face pressure to unify heritage language learning around a single standard variety.In this case, language standardization has prevented Tibetan children from accessing forms of linguistic belonging available to adults. It has also unintentionally contributed to a language shift away from Tibetan mother languages and to dominant languages, including Mandarin and English.In such situations, immigrant and minority children face two forms of linguistic marginalization. First, a nation’s official languages exclude their recognized heritage language. Second, a standard language spoken within their community excludes their native, mother languages.

      This connects to my own experience learning English while continuing to speak Spanish. Supporting children's home languages creates better learning environments.

  2. bafybeif5xmfuc2ftek4jt7qzg4mwon2qoeic2pn5ij4j23fvymb24dytaa.ipfs.inbrowser.link bafybeif5xmfuc2ftek4jt7qzg4mwon2qoeic2pn5ij4j23fvymb24dytaa.ipfs.inbrowser.link
    1. Appendix E. Enabling Direct Restore from Object Storage In this article

      It will be better to add some NOTE that this works only when Object Storage is configured as Direct Mode in VBR side.

    1. To increase restore speed

      It doesn't increase restore speed, rather increase time when agent tries to find out which VBR Endpoint is accessible and mounts the backup to local machine.

    1. You can limit bandwidth consumption and restrict network connection usage for Veeam Agent for Mac backup jobs. Limiting bandwidth consumption prevents backup jobs from using all available bandwidth and helps ensure that enough traffic is provided for other network operations.

      Better we align with Windows text

    1. The proband (Patient #20

      Case#: Female, family #5, Patient #20

      DiseaseAssertion: STGD

      FamilyInfo: Proband's sister presented with same clinical prognosis. Sister diagnosed with pattern dystrophy and photoaversion at age 57, with difficulty seeing at night. Sister has nuclear sclerotic and cortical cataracts in both eyes.

      CasePresentingHPOs: HP:0000662, HP:0000603, HP:0000603, HP:0000493

      CaseHPOFreeText: Proband presented with localized blur at age 62, (late onset) in her left eye. BVCA 20/20-3 and 20/20-2 at age 70. Also has macular lesions with stage 2 fundus flecks.

      CaseNotHPOs: n/a

      CaseNotHPOFreeText: n/a

      Genotyping Method: Genotyping performed at Columbia University, sequencing technology used is not disclosed.

      PreviouslyPublished: n/a

      Variant: p.N18681, IVS36:c.5196+1G>A

      ClinVar: M2) 99067, M6) 99351

      CAID: N/A

      SupplementalData: Fig 1: Pedigree illustrating ABCA4 variants and the associated Stargardt phenotype for 5 families. Proband Labeled w/ white arrow for each family. Fig 2: retinal scan measuring melanin in 4 patients of family 2. Panel shows bull's-eye ring of RPE atropy. Fig 3: Macular SD-OCT line profile from b-scans. Reflectivity plotted against function of retinal depth. Table 1: table shows patients with p.N18681 variant, type of mutation, and pathogenicity class. Table 2: Patients, age on-set and first symptom

    1. F17-003

      Case#: Patient 225, Female, age of onset 7 y.o, Poland

      DiseaseAssertion: STGD-1

      FamilyInfo: no given family information.

      CasePresentingHPOs: HP:0007722, HP:0000608, HP:0025158

      CaseHPOFreeText: RPE atrophy, macular degeneration, central hyper-autofluorescence in fundus autofluorescence

      CaseNotHPOs: n/a

      CaseNotHPOFreeText: n/a, non-proband identified HPO's mentioned, but not assignable to individual proband.

      Genotyping Method: DNA isolated from peripheral blood from patients and relatives via MagNA Pure 24, samples screened with MIPs targeting 108 genes involved in pathogensis of IRD's. PCR completed on library, analysed with NGS fragment analysis kit.

      PreviouslyPublished: yes

      Variant: c.[1622T>C;3113C>T]

      ClinVar: 99067, 7894

      CAID: n/a

      gnomeAD 0.0001266 allele frequency

      SupplementalData: Fig1: List of families displaying pseudo-dominant inheritance. Fig2: Number of alleles for most common variants.

    1. 4.4. Disease Course in Patients Harbouring p.(Gly1961Glu) or p.(Asn1868Ile) Allele

      It is known that patients harbouring p.(Gly1961Glu) or p.(Asn1868Ile) allele share some common clinical characteristics and present with a milder disease phenotype than patients carrying other alleles. A typical feature of patients with p.(Gly1961Glu) is BEM, which is otherwise present in around 20% of all STGD1 patients.

    2. 4.4. Disease Course in Patients Harbouring p.(Gly1961Glu) or p.(Asn1868Ile) Allele

      It is known that patients harbouring p.(Gly1961Glu) or p.(Asn1868Ile) allele share some common clinical characteristics and present with a milder disease phenotype than patients carrying other alleles. A typical feature of patients with p.(Gly1961Glu) is BEM, which is otherwise present in around 20% of all STGD1 patients.

    1. 10/22/MClinical, geneticPartial deletion, protein truncating mutation++++−

      PatientID: 10

      KindredID: 10

      Case: M, 22Y0M, Caucasian

      DiseaseAssertion: VHL

      FamilyInfo: Positive family history.

      CohortInfo: Case series of 14 patients with definite or presumed von Hippel-Lindau disease and retinal vascular proliferation.

      CasePresentingHPOs: HP:0002011, HP:0001732, HP:0007850, HP:0012210, HP:0009711 (Morphological central nervous system abnormality, abnormality of the pancreas, retinal vascular proliferation, abnormal renal morphology, retinal capillary hemangioma)

      CaseHPOFreeText: At age 12, his left eye had a large juxtapapillary vascular complex extending from the nerve along the superior arcade that was associated with dense fibrovascular tissue. During the next 29 months, this fibrovascular complex enlarged, extending into the center of the macula and causing a decrease in visual acuity to 20/80. The patient underwent pars plana vitrectomy and membrane peel. When the patient was reexamined 7 years after his surgery, there was no recurrence of the lesion. In the patient's contralateral right eye, 3 typical peripheral RCHs were observed during this 7-year follow-up. At age 14 years, he was also noted to have a small patch of superficial retinal vessels in the inferior retinal quadrant near the equator of the right eye that resembled an arteriovenous anastomosis.

      CaseNotHPOs: N/A

      CaseNotHPOFreeText: N/A

      CasePreviousTesting: N/A

      PreviouslyPublished: N/A

      SupplementalData: N/A

      Variant: Partial Deletion

      LegacyVariant: N/A

      CaseProblemVariantFreeText: N/A

      ClinVarID: N/A

      CAID: N/A

      gnomAD: N/A

      VariantEvidence: N/A

      MutationType: deletion

      CivicName:Null (Partial deletion)

      MultipleGeneVariants: N/A

    1. Complete deletion

      GroupID/ KindredID: 28

      PatientID: II

      Case: 24Y0M, Sex unknown, Polish

      DiseaseAssertion: VHL

      FamilyInfo: Family 28 consisted of 4 members; one relative (age 27Y) was found with multiple cerebellar HABs at 25Y, and multiple spinal HABs, another (40Y) diagnosed at 30Y with multiple cerebellar HABs, and multiple spinal HABs, and a final relative (62Y), diagnosed with a single HAB.

      CasePresentingHPOs: HP:0006880 (cerebellar hemangioblastoma)

      CaseHPOFreeText: Patient was diagnosed with a single cerebellar HAB at 24Y

      GroupNotHPOs: HP:0009713; HP:0009711; HP:0005584 (spinal hemangioblastoma, retinal capillary hemangioma; renal cell carcinoma)

      CaseNotHPOFreeText: N/A

      CasePreviousTesting: N/A

      PreviouslyPublished: N/A

      SupplementalData: N/A

      Variant: complete deletion of VHL gene

      CaseProblemVariantFreeText: N/A

      LegacyVariant: N/A

      ClinVarID: N/A

      CAID: N/A

      gnomAD: N/A

      VariantEvidence: Haplotype analysis information in Table 3

      MutationType: deletion

      CivicName: deletion

      MultipleGeneVariants: N/A

    1. 4 partial deletions

      GroupID/ KindredID: 22

      Case: Age unknown, Sex unknown, Chinese

      DiseaseAssertion: hemangioblastoma

      FamilyInfo: No family history. Genetic testing of their parents confirmed a de novo mutation.

      CasePresentingHPOs: HP:0010797 (hemangioblastoma)

      CaseHPOFreeText: N/A

      CaseNotHPOs: HP:0002666; HP:0005584; HP:0009711; HP:0001732 (pheochromocytoma; renal cell carcinoma; retinal capillary hemangioma; pancreatic lesion)

      CaseNotHPOFreeText: N/A

      CasePreviousTesting: N/A

      PreviouslyPublished: N/A

      SupplementalData: N/A

      Variant: Exon 1 deletion

      LegacyVariant: N/A

      CaseProblemVariantFreeText: N/A

      ClinVar: N/A

      CAID: N/A

      gnomAD: N/A

      VariantEvidence:N/A

      MutationType: exon_loss_variant

      CivicName: Exon 1 Deletion

      MultipleGeneVariants: N/A

    1. In four families, partial deletions of one or more exons were detected by Southern blot analysis

      PatientID: unknown (brother)

      KindredID: A

      Case: M (deceased), Age unknown, Turkish

      DiseaseAssertion: assumed VHL

      FamilyInfo: 5 family members are grouped with this deletion (ages 16-37Y; mean 31Y). VHL was clinically diagnosed in, at minimum, the proband. See Fig. 2 for family pedigree. This patient and his one brother were not tested for the DNA deletion however it is assumed by the authors due to the segregation observed.

      CasePresentingHPOs: HP:0010797 (hemangioblastoma)

      CaseHPOFreeText: N/A

      CaseNotHPOs: HP:0009711; HP:0002666; HP:0005584 (retinal capillary hemangioma; pheochromocytoma; renal cell carcinoma)

      CaseNotHPOFreeText: N/A

      CasePreviousTesting: N/A

      PreviouslyPublished: N/A

      SupplementalData: N/A

      Variant: Deletion of exons 1 and 2

      CaseProblemVariantFreeText: N/A

      LegacyVariant: N/A

      ClinVarID: N/A

      CAID: N/A

      gnomAD: N/A

      VariantEvidence: not tested for the DNA deletion however it is assumed by the authors due to the segregation observed.

      MutationType: exon_loss_variant

      CivicName: exon 1-2 deletion

      MultipleGeneVariants: N/A

    1. GroupID/ KindredID: F28

      Case: Sex unknown, 23Y0M, Spanish

      DiseaseAssertion: assumed VHL

      FamilyInfo: familial antecedents found; unknown of any additional features of index case relatives

      CasePresentingHPOs: HP:0000107; HP:0009711; HP:0006880; HP:0006770; HP:0002666 (renal cysts, retinal capillary hemangioma, cerebellar hemangioblastoma, clear cell renal cell carcinoma, pheochromocytoma)

      CaseHPOFreeText: unilateral pheo, bilateral ccRCC and multiple lesions were found with the patient’s renal cysts. Age of onset is 23Y0M.

      CaseNotHPOs: HP:0001737 (pancreatic cysts)

      CaseNotHPOFreeText: N/A

      CasePreviousTesting: N/A

      PreviouslyPublished: N/A

      SupplementalData: N/A

      Variant: rearrangement (unspecified)

      LegacyVariant: N/A

      CaseProblemVariantFreeText: N/A

      ClinVarID: N/A

      CAID: N/A

      gnomAD: N/A

      VariantEvidence: Southern blot positive

      MutationType: rearrangement_region

      CivicName: Rearrangement

      MultipleGeneVariants: N/A

    1. PatientID: 246

      KindredID: 246

      Case: Sex Unknown, 26Y0M, Ethnicity Unknown

      DiseaseAssertion: VHL

      FamilyInfo: N/A

      CasePresentingHPOs: HP:0010797; HP:0009711; HP:0006770 (Hemangioblastoma, Retinal capillary hemangioma, Clear cell renal cell carcinoma)

      CaseHPOFreeText: CNS and retinal hemangioblastoma and clear cell renal cell carcinoma.

      CaseNotHPOs: HP:0002666; HP:0000107; HP:0001732 (Pheochromocytoma, Renal cyst, Abnormality of the pancreas)

      CaseNotHPOFreeText: N/A

      CasePreviousTesting: N/A

      PreviouslyPublished:N/A

      SupplementalData: Table S1, S2 and S3

      Variant: Partial Deletion (0.7kb)

      LegacyVariant: N/A

      CaseProblemVariantFreeText: N/A

      ClinVarID: N/A

      CAID: N/A

      gnomAD: N/A

      VariantEvidence: N/A

      MutationType: deletion

      CivicName: Partial deletion of 0.7kb

      MultipleGeneVariants: N/A

    1. 16 different families

      PatientID: IX-I

      KindredID: 9

      Case: M, 68Y0M, Ethnicity Unknown

      DiseaseAssertion: VHL

      FamilyInfo: No family history reported.

      CasePresentingHPOs: HP:0006748; HP:0002668; HP:0030405; HP:0001737; HP:0006880 (adrenal pheochromocytoma, paraganglioma, pancreatic endocrine tumor, pancreatic cyst, cerebellar hemangioblastoma)

      CaseHPOFreeText: Patient presented with adrenal pheochromocytoma, paraganglioma, pancreatic endocrine tumor, pancreatic cyst, and cerebellar hemangioblastoma at age 55Y0M and was diagnosed with VHL type 2A (see table 2).

      CaseNotHPOs: HP:0000107; HP:0005584 (renal cyst; renal cell carcinoma)

      CaseNotHPOFreeText: Patient negative for renal cell carcinoma and renal cysts.

      CasePreviousTesting: N/A

      PreviouslyPublished: N/A

      SupplementalData: N/A

      Variant: complete deletion

      LegacyVariant: N/A

      CaseProblemVariantFreeText: N/A

      ClinVarID: N/A

      CAID: N/A

      gnomAD: N/A

      VariantEvidence: N/A

      MutationType: deletion

      CivicName: NULL (deletion)

      MultipleGeneVariants: N/A

    1. Comparably, the result of the CNTNAP2 rs2710102 associations with ASD studies (n = 4+current study) was also not significant [n = 7276; p = 0.26; OR = 1.028 (95 % CI 0.98–1.08), see Suppl. Table S6, Fig. 1b]

      Performed a meta-analysis of all data published on the rs2710102 CNTNAP2 variant up to 2015, and included Anney et al, 2012 (PMID: 20663923), Toma et al., 2013 (PMID: 23277129), Sampath et al, 2013 (PMID: 24147096), Poot et al., 2014 (PMID: 25337070 ), plus there current cohort.

      No significant increase in cases vs controls for CNTNAP2 rs2710102

      While this data does help support Contradictory evidence for CNTNAP2 involvement in autism, the lack of reporting of the total number of controls, as well as the number of cases and controls with the SNP prevents me from curating this information within the ClinGen gene curation interface.

    2. he meta-analysis of the CNTNAP2 rs7794745 associations with ASD studies (n = 5+current study) did not result in significant association with ASD in general [n = 8576; p = 0.112; OR = 1.023 (95 % CI 0.99–1.05); see Suppl. Table S5, Fig. 1a]. Since we detected some heterogeneity for the SNP rs7794745 according to the funnel plot (see Suppl. Figure S10), we performed an additional meta-analysis synthesis using the random effect model. However, no significant association with ASD was observed and OR were very similar as for the fixed-model [OR = 1.081 (95 % CI 0.976–1.196), p = 0.133; for details Suppl. Table S9a].

      Performed a meta-analysis of all data published on the rs7794745 CNTNAP2 variant up to 2015, and included Arking et al., 2008 (PMID:18179894), Li et al 2010 (PMID: 20414140), Anney et al, 2012 (PMID: 20663923), Toma et al., 2013 (PMID: 23277129), Sampath et al, 2013 (PMID: 24147096), plus there current cohort.

      No significant increase in cases vs controls for CNTNAP2 rs7794745

    1. Similar

      CaseControlLabel: Toma 2013 CNTNAP2 rs7794745 analysis

      CaseCohortLabel: 322 Autism patients

      ControlCohortLabel: 524 controls

      CaseCohortDisease: MONDO:0005260 (Autism)

      CaseCohortPhenotype:

      CaseCohortPhenotypesFreeText: Authors indicate that all individuals fulfilled DSM-IV criteria for autism or Asperger disorder. Additionally, some cases included where diagnosed with pervasive developmental disorder they may or may not have been characterized based on ADI-R and ADOS-G.

      CaseControlNOTPhenotype:

      CaseControlNOTPhenotypeFreeText:

      CaseDemographics: 269 men, 53 women, Average age= = 17y.o. of Spanish and/or Caucasian descent.

      ControlDemographics: controls were noted to be sex-matched and unrelated, however no numbers or ages were indicated. Sample obtained from the Blood and Tissues Bank of Hospital Universitari Vali d'Hebron.

      CaseGenotypingMethod: One SNP, rs7794745, was genotyped using PCR-RFLP, from DNA obtained from peripheral blood lymphocytes (salting out method).

      ControlGenotypingMethod: One SNP, rs7794745, was genotyped using PCR-RFLP, from DNA obtained from peripheral blood lymphocytes (salting out method).

      CasePower: 312/322

      ControlPower: 505/524

      CaseAddInfo: Other genetic variation is not noted.

      ControlAddInfo: Other genetic variation is not noted.

      CaseControlStudyType: Single Variant Analysis

      CaseControlDetectionMethod: Cases and controls genotyped for single variant

      CaseControlStats: p-value= 0.87

      CaseControlBias: Cases and controls are ethinically matched, and noted to be sex-matched.

      CaseControlComments:

    1. (GRIN2B, rs7301328)

      Case: Multiple individuals with the same GRIN2B variant (the mean age was 76 ± 11 SD years (Range: 38–101), the mean age at PD onset was 67 ± 12 SD years (Range: 32–94), and 421 patients (63%) were male).

      Disease Assertion: Parkinson's disease

      FamilyInfo: None

      CasePresentingHPOs: HP:0001300

      CaseHPOFreeText: Parkinson's disease

      CaseNotHPOs: HP:0001250

      CasePreviousTesting: Polymorphism genotyping on Sequenom iPLEX ® platform

      PreviouslyPublished: No

      Variant: rs7301328

      ClinVarID: 129205

    1. Table 2.

      Case#:Not assigned

      DiseaseAssertion:MPS type I with evidence of cognitive impairment defined as a score of at least 1 SD below mean on IQ testing or in one domain of neuropsychological function; attenuated MPS I

      FamilyInfo:No information

      CasePresentingHPOs:HP:0100543

      CaseHPOFreeText:

      CaseNotHPOs:HP:0032557

      CaseNotHPOFreeText:Subjects were excluded if they had undergone hematopoietic stem cell transplantation, were unable to comply with study procedures, had significant lumbar pathology precluding access to the intrathecal space via lumbar puncture, or significantly impaired spinal CSF flow as detected on a nuclear medicine flow study as part of screening.

      CaseEnzymeAssay:No information

      CaseUrineGAGs:No information

      CaseERT:Yes

      CaseBMT:No (exclusion criterion)

      Variant1:L238Q

      Variant1ClinVarID:265418

      Variant1CAID:

      Variant2:W402X

      Variant2ClinVarID:

      Variant2CAID:

      AdditionalVariants:N/A

      ParentalGenotype:Not specified

      PreviouslyPublished

    2. Table 2.

      Case#:Not assigned

      DiseaseAssertion:MPS type I with evidence of cognitive impairment defined as a score of at least 1 SD below mean on IQ testing or in one domain of neuropsychological function; attenuated MPS I

      FamilyInfo:No information

      CasePresentingHPOs:HP:0100543

      CaseHPOFreeText:

      CaseNotHPOs:HP:0032557

      CaseNotHPOFreeText:Subjects were excluded if they had undergone hematopoietic stem cell transplantation, were unable to comply with study procedures, had significant lumbar pathology precluding access to the intrathecal space via lumbar puncture, or significantly impaired spinal CSF flow as detected on a nuclear medicine flow study as part of screening.

      CaseEnzymeAssay:No information

      CaseUrineGAGs:No information

      CaseERT:Yes

      CaseBMT:No (exclusion criterion)

      Variant1:L238Q

      Variant1ClinVarID:265418

      Variant1CAID:

      Variant2:63delC

      Variant2ClinVarID:

      Variant2CAID:

      AdditionalVariants:N/A

      ParentalGenotype:Not specified

      PreviouslyPublished

    3. Table 2.

      Case#:Not assigned

      DiseaseAssertion:MPS type I with evidence of cognitive impairment defined as a score of at least 1 SD below mean on IQ testing or in one domain of neuropsychological function; attenuated MPS I

      FamilyInfo:No information

      CasePresentingHPOs:HP:0100543

      CaseHPOFreeText:

      CaseNotHPOs:HP:0032557

      CaseNotHPOFreeText:Subjects were excluded if they had undergone hematopoietic stem cell transplantation, were unable to comply with study procedures, had significant lumbar pathology precluding access to the intrathecal space via lumbar puncture, or significantly impaired spinal CSF flow as detected on a nuclear medicine flow study as part of screening.

      CaseEnzymeAssay:No information

      CaseUrineGAGs:No information

      CaseERT:Yes

      CaseBMT:No (exclusion criterion)

      Variant1:65delC

      Variant1ClinVarID:

      Variant1CAID:

      Variant2:H240R

      Variant2ClinVarID:947028

      Variant2CAID:

      AdditionalVariants:N/A

      ParentalGenotype:Not specified

      PreviouslyPublished

    1. Patient 1

      Case:Patient 1, Male, 4 y/o, German and Thai

      DiseaseAssertion:FOXG1 syndrome

      FamilyInfo: Non-consanguineous parents. De novo. Patient 1 has maternally inherited duplications at 15q21.12(47292250-47309868)x3 and Xp22.31(6451676-8115124)x3 of unknown clinical significance.

      ParentalGenotype:Not provided

      CasePresentingHPOs:HP:0002197, HP:0005484, HP:0011344, HP:0001252, HP:0020045,HP:0000540, HP:0010808, HP:0002019, HP:0003781, HP:0012741, HP:0002421, HP:0000748, HP:0003763, HP:5200017, HP:0025112, HP:0000957

      CaseHPOFreeText: At 19 months, could not sit independently and no speech development. Stereotypic hand and head movements. High pain tolerance. X-ray at age 2 showed broad ribs. Improved head control by age 4.

      CaseNotHPO:HP:0410263, HP:0002376

      CaseNotHPOFreeText:Dysmorphic features. breathing abnormalities

      CasePreviousTesting:No previous testing

      PreviouslyPublished:Not previously published

      GenotypingMethod:Mate-pair sequencing, Sanger sequencing

      Gene:FOXG1

      Variant:46,XY,t(9;14)(q22.3;q11.2)dn

      HGVS:Not provided

      ClinVarID:426096

      gnomAD:Not provided

    2. Patient 3

      Case:Patient 3, Male, 1 y/o, Brazilian

      DiseaseAssertion: FOXG1 syndrome

      FamilyInfo: Non-consanguineous parents

      ParentalGenotype: No family testing mentioned.

      CasePresentingHPOs: HP:0000252, HP:0001252, HP:0000486

      CaseHPOFreeText: Developed microcephaly within 1 year. Developmental delay, unable to sit without assistance at 1 year old. Central Nervous System (CNS) MRI showed a volumetric reduction of the cerebral parenchyma, hypomyelination, and hypoplasia of corpus callosum mainly involving the anterior part and prominent lateral ventricles

      CaseNotHPO: HP:0030917

      CaseNotHPOFreeText: Normal length, weight, HC at birth. No retinal abnormalities.

      CasePreviousTesting: Not asserted.

      PreviouslyPublished: Not previously published.

      GenotypingMethod: Sanger sequencing, mate-pair sequencing

      Gene: FOXG1

      Variant: 46,XY,t(2;14)(q37;q11.2)dn

      HGVS: Not provided.

      ClinVarID:426094

      gnomAD:Not provided

    3. Patient 2

      Case:Patient 2, Female, 5 y/o, Brazilian

      DiseaseAssertion:quadriplegic cerebral palsy with axial hypotonia and distal spasticity

      FamilyInfo: Non-consanguineous parents. Parents had normal karyotypes. Older sister was clinically normal.

      ParentalGenotype: Normal karyotypes for both parents.

      CasePresentingHPOs: HP:0002098, HP:0010959, HP:0008755, HP:0002090, HP:0032661, HP:0001274, HP:0011344, HP:0001249, HP:0001344, HP:0000748, HP:0034435, HP:0000486, HP:0000483, HP:0200136

      CaseHPOFreeText: Dilated superficial blood vessels seen on brain MRI.

      CaseNotHPO:N/A

      CaseNotHPOFreeText: N/A

      CasePreviousTesting: Not asserted

      PreviouslyPublished: Not previously published

      GenotypingMethod:Mate-pair sequencing, Sanger sequencing, chromosome analysis

      Gene: FOXG1

      Variant: 46,XX,t(4;14)(q27;q13)dn

      HGVS: Not provided

      ClinVarID: 426095

      gnomAD: Not provided

    1. (A) With ARL2BP, photoreceptor OS develops normally, with correct DMT formation, axoneme extension, and disk alignment.

      Selection A- Gene with same function implicated in disease

      ExperimentName: Moye et al. 2018

      GOterm: GO:0051179 (localization)

      FunctionalEvidence: It was discovered that without ARL2BP, the photoreceptor OS does not develop normally or have correct DMT formation, axoneme extension, and disk alignment.

      HGNC: HGNC:10295 (RPGR), HGNC:10263 (RP1), HGNC:6816 (MAK)

      GeneEvidence: ARL2BP, RPGR, RP1 and MAK are all localized within the ciliary basal body and rootlet.

      SharedGeneImplication: RPGR was classified as having a "definitive" relationship with RP 15 by the Retina GCEP. Refer to the curation for further details. RP1 was classified as having a "definitive" relationship with RP1 by the Retina GCEP. Refer to the curation for further details. MAK was classified as having a "definitive" relationship with MAK-related retinopathy by the Retina GCEP. Refer to the curation for further details.

    1. RSPH1 and RSPH9, here highlighted, are known to be components of the RS head in human respiratory cilia

      Selection A- Gene with same function implicated in disease

      ExperimentName: Aprea et al. 2023

      GOterm: GO:0001535 (radial spoke head)

      FunctionalEvidence: Taken together, our results indicate that the ruler proteins CCDC39 and CCDC40, the radial spoke head proteins RSPH1 and RSPH9 as well as the central pair associated apparatus proteins SPEF2 and HYDIN have a similar functional role in motile cilia and sperm flagella, as hypothesized in Figure 1.

      HGNC: HGNC:12371

      GeneEvidence: PMID: 24518672 documents how RSPH1 and RSPH9, which encode homologs of components of the ‘head’ structure of ciliary radial spoke complexes identified in Chlamydomonas, cause clinical phenotypes that appear to be indistinguishable except at the gene level. They both also have a relationship with PCD.

      SharedGeneImplication: The gene-disease relationship between RSPH1 and PCD has been established and published in multiple papers. See Motile Ciliopathy GCEP curations of RSPH1 for more information.

    1. (B) Schematic of the axonemal 96 nm repeat unit, including the composition of the IDA subspecies and the location of axonemal components. DNAH7, DNAH6, and DNAH1 are highlighted in purple, orange, and yellow, respectively. The stalks of the IDA dynein heavy chains are highlighted in blue or pink, depending on whether they are associated with DNALI1 or CETN2. The CCDC39/40 ruler complex is depicted in light and dark green.

      Selection A- Gene with same function implicated in disease

      ExperimentName: Wilken et al. 2024

      GOterm: GO:0036156 (inner dynein arm).

      FunctionalEvidence: Both DNAH6 and DNAH1 are both IDA components and localized to the same region.

      HGNC: HGNC:2940 (DNAH1).

      GeneEvidence: Both DNAH6 and DNAH1 are both IDA components and localized to the same region.

      SharedGeneImplication: DNAH1 was classified as having a "limited" relationship with PCD 37 and a “defintive” relationship with spermatogenic failure 18 by the Motile Ciliopathy GCEP. Refer to the curation for further details.

    1. Four genes cause PCD through radial spoke defects - RSPH9 (Chr 6), RSPH4A (Chr 6), RSPH3 (Chr 6), and RSPH1 (Chr 21). These mutations produce various defects at the junction of the radial spoke heads and the central apparatus, resulting in diverse ciliary ultrastructure changes on TEM.

      Selection A - Gene with same function implicated in disease

      ExperimentName: Shapiro 2007 et al.

      GOterm: GO:0051179 (localization)

      FunctionalEvidence: RSPH3 and RSPH9 are both localized to radial spoke heads and have a similar phenotype profile.

      HGNC: HGNC:21057.

      GeneEvidence: PMID: 20301301 documents how RSPH9 and RSPH3, which encode homologs of components of the ‘head’ structure of ciliary radial spoke complexes identified in Chlamydomonas, cause clinical phenotypes that appear to be indistinguishable except at the gene level. They both also have a relationship with PCD.

      SharedGeneImplication: The gene-disease relationship between RSPH9 and PCD has been established and published in multiple papers. See Motile Ciliopathy GCEP curations of RSPH9 for more information.

    1. The source focuses on the clinical, immunological, and molecular features of Severe Combined Immune Deficiency (SCID) based on a multi-institutional experience from India. The specific gene studied in this case is IL2RG. The HGNCID for this gene is HGNC:6016. The sources provide information about a 6-month-old male patient with SCID. The patient's family history suggests a possible X-linked inheritance pattern. The sources note the absence of autoimmunity and malignancies in the patient. Information regarding the patient's variant, previous testing, and publication history is not available in the sources.

    1. 52

      Case#:Patient 52, female, 3 years old

      DiseaseAssertion:Neonatal/Infantile Epileptic Encephalopathy (NIEE)

      FamilyInfo:DeNovo. The family is Chinese

      ParentalGenotype:The authors only conducted singleton and not trio-based exome sequencing so the parents' exomes were not sequenced.

      CasePresentingHPOs:HP:0011344, HP:0002069, HP:0007359, HP:0011097, HP:0100704, HP:0001332, HP:0002072, HP:0012171.

      CaseHPOFreeText:Patient 52 presents with severe global developmental delay and epilepsy.

      Patient 52 has generalized tonic/clonic/tonic-clonic seizures, focal seizures and spasms. Patient 52's seizure onset occurred at 3 months old.

      Patient 52 has cortical visual impairment (CVI), dystonia, chorea, and hand-washing sterotypies.

      Patient History

      @ 3 months - Patient 52 had generalized tonic/clonic/tonic-clonic seizures.

      Patient 52 was on 3 antiepileptic drugs at most recent follow-up visit which reduced seizure frequency by >50%.

      CaseNotHPOs:Not provided

      CaseNotHPOFreeText:Not provided

      CasePreviousTesting:The authors selected a cohort of 31 patients with seizure cryptogenic Neonatal/Infantile Epileptic Encephalopathy (NIEE) and seizure onset before 24 months.

      Exclusion criteria included: (1) Patients with a definite history of brain insult, malformation of cortical development, neurocutaneous and syndromal disorders, and confirmed or highly suspected neurometabolic disorders based on clinical and biochemical markers. (2) Patients with Dravet syndrome and epilepsy at infancy with migrating focal seizure were also excluded because the majority of variants are detected in the SCN1A (>85%) and KCNT1 (approximately 50%) genes.

      Formal neuropsychological testing or best clinical assessment was used to classify patient development or intelligence.

      PreviouslyPublished:Not previously published

      GenotypingMethod:Whole Exome Sequencing (WES) variant results were filtered in a panel of 430 epilepsy-associated genes. After selection of variants from the 430-gene panel, the synonymous variants, variants with variant frequency <10%, and variants with allele frequency >1% were removed.

      Gene:CDKL5

      Variant:NM_003159.2 c. 1849delC (p. Arg617Valfs*4)

      The authors state that the variant is a heterozygous frameshift deletion.

      The authors state that this is a novel variant and is pathogenic.

      HGVS:Not provided

      ClinVarID:Not found

      CAID:Not found.

      gnomAD:Not found

      MultipleGeneVariants:Not provided

    2. 7

      Case#:Patient 7, male, 5 years old

      DiseaseAssertion:Neonatal/Infantile Epileptic Encephalopathy (NIEE)

      FamilyInfo:The family is French/Chinese

      ParentalGenotype:The variant was inherited from Patient 7's asymptomatic mother.

      CasePresentingHPOs:HP:0010864, HP:0012758, HP:0000729, HP:0007359, HP:0002069, HP:0032794, HP:0001250, HP:0000252.

      CaseHPOFreeText:Patient 7 presents with severe intellectual disability, developmental slowdown, and various seizure types. Patient 7 has Autistic Spectrum Disorder (ASD) and microcephaly.

      Patient History

      @ 12 months - Patient 7 presented with seizures.

      Patient 7 developed additional seizure types including: focal seizures with/without generalization, generalized tonic/clonic/tonic-clonic seizures, myoclonic seizures, and hypomotor seizures.

      Patient 7 was on two antiepileptic drugs at most recent followup visit which reduced seizure frequency by >50%.

      CaseNotHPOs:Not provided

      CaseNotHPOFreeText:Not provided

      CasePreviousTesting:The authors selected a cohort of 31 patients with seizure cryptogenic Neonatal/Infantile Epileptic Encephalopathy (NIEE) and seizure onset before 24 months.

      Exclusion criteria included: (1) Patients with a definite history of brain insult, malformation of cortical development, neurocutaneous and syndromal disorders, and confirmed or highly suspected neurometabolic disorders based on clinical and biochemical markers. (2) Patients with Dravet syndrome and epilepsy at infancy with migrating focal seizure were also excluded because the majority of variants are detected in the SCN1A (>85%) and KCNT1 (approximately 50%) genes.

      Formal neuropsychological testing or best clinical assessment was used to classify patient development or intelligence.

      PreviouslyPublished:Not previously published

      GenotypingMethod:Whole Exome Sequencing (WES) variant results were filtered in a panel of 430 epilepsy-associated genes. After selection of variants from the 430-gene panel, the synonymous variants, variants with variant frequency <10%, and variants with allele frequency >1% were removed.

      Gene:SLC9A6

      Variant:Hemizygous splice site NM_001042537.1 c. 794-2A>G was assessed by the authors to be likely pathogenic.

      HGVS:Not provided

      ClinVarID:Not found

      CAID:CA414750320

      gnomAD:Not found

      MultipleGeneVariants:Not provided

    1. Case  1

      Case: Patient 1, female, 2 years old

      DiseaseAssertion: Rett Syndrome

      FamilyInfo: DeNovo. Patient 1 was the second child of a nonconsanguineous Chinese couple. Patient 1 was born at full term with a birth weight of 3.83 kg.

      ParentalGenotype(s):Not provided

      CasePresentingHPOs:HP:0011344, HP:0001250, HP:0001252, HP:0000252, HP:0000748, HP:0003763, HP:0005469, HP:0000486, HP:0012171.

      CaseHPOFreeRext:Patient 1 presents with severe global developmental delay, epilepsy, and hypotonia. The Patient History contains all of Patient 1's phenotypes along with the progression of her condition.

      Patient History

      @ 3 months - Patient 1 had microcephaly (head circumference was less than 3rd percentile with a body weight and body height at 75th percentile). Patient 1 had hypotonia.

      Patient 1 was given a metabolic screening, a muscle enzyme, and a brain tomography with normal results.

      @ 6 months - Patient 1 had microcephaly, flat occiput, right divergent squint, and hypotonia. No syndromal diagnosis could be ascertained at that time.

      @ 2 years - Patient 1 had epilepsy.

      Patient 1 had severe global development delay.

      Patient 1 was given an EEG which showed nonspecific background slowing, but no epileptiform abnormalities. Patient 1 was also given a brain MRI which showed mild thinning of corpus callosum without major structural defect. Patient 1 had no developmental regression but she developed stereotypical hand movements, bruxism, and occasional outburst of laughter.

      Based on these phenotypes, Angelman/Rett Syndrome was suspected.

      CaseNOTHPOs: HP:0002353.

      CaseNOTHPOFreeText: Patient 1 was given an EEG which detected no epileptic abnormalites.

      CasePreviousTesting: Genetic investigations (including methylation-specific multiplex ligation-dependent probe amplification (MS-MLLPA)) was conducted for Angelman Syndrome, UBE3A gene, MECP gene, and array CGH. These studies were negative. A FOXG1 related disease was suspected

      PreviouslyPublished:Not previously published

      GenotypingMethod: A FOXG1 gene test showed a de novo frameshift pathogenic mutation FOXG1 {NM_005249.3} c. 396_397ins26; FOXG1{NP_005240.3} :(p. Gly133TRPfs*68) which confirmed the diagnosis of a FOXG1 related congenital variant of Rett Syndrome.

      Gene:FOXG1

      Variant: (NM_005249.3) c. 396_397ins26 (p. Gly133Trpfs*68)

      HGVS:Not provided

      ClinVarID:Not found

      CAID:Not found

      gnomAD:Not found

    2. Case  2

      Case:Patient 2, female, Chinese

      DiseaseAssertion:Global delay

      FamilyInfo:Patient 2 was the first child of nonconsanguineous Chinese couple born at 38-week gestation. The mother had gestational diabetes mellitus that required insulin therapy.

      ParentalGenotype:Not provided

      CasePresentingHPOs:HP:0000365, HP:0020049, HP:0000252, HP:0001263, HP:0012171, HP:0000154, HP:0000708, HP:0003763, HP:0002376, HP:0012433.

      CaseHPOFreeText:

      Patient History

      @ birth - Patient 2 presented with mild grade bilateral hearing impairment and left divergent squint diagnosed at birth.

      @ Follow-up visit - Patient 2 had microbrachycephaly and global developmental delay. Brain MRI, metabolic screening and array CGH were normal.

      @ 1 year - Patient 2 had stereotypical handwashing movement.

      There was no clinical or electrical seizure.

      Patient 2 has craniofacial features like microbrachycephaly, wide mouth, divergent squint, and behavioral phenotype.

      @ 1.5 years - Patient 2 had bruxism and developmental regression. Patient 2 also had loss of some motor and social skills.

      CaseNotHPOs:HP:0001250

      CaseNotHPOFreeText:Patient 2 has no seizures.

      CasePreviousTesting:Not provided

      PreviouslyPublished:Not previously published

      GenotypingMethod:Not provided

      Gene:MECP2 (MN_004992.3) (NP_004983.1)

      Variant:c. 808C>T (p. Arg270*)

      HGVS:Not provided

      ClinVarID:Not found

      CAID:CA172577

      gnomAD:Not found

      MultipleGeneVariants:NA

    1. first patient

      Case:Patient 1, female, 11 years old

      DiseaseAssertion:Rett Syndrome - Atypical Variant

      FamilyInfo:De Novo. Patient 1's parents were healthy. When Patient 1 was clinically examined, the head circumferences of her mother and father were 53 cm (P25) and 59 cm (P90), respectively. The parents had normal weight.

      ParentalGenotype:Both parents were sequenced for Patient 1's mutation, and the deletion was not detected.

      CasePresentingHPOs:HP:0001249, HP:0001513, HP:0001626, HP:0000256, HP:0010465, HP:0012758, HP:0000750, HP:0012433, HP:0002591, HP:0001250, HP:0020174, HP:0000316, HP:0000336, HP:0000431, HP:0000377, HP:0000470, HP:0001156, HP:0001500.

      CaseHPOFreeText:

      Patient history

      @ birth - Patient 1 was born at 38 weeks of gestation after an uncomplicated pregnancy. Her weight was 3,390 g (P69), height 51 cm (P60), and her head circumference was 36 cm (P90).

      @ 6 months - Patient 1 experienced developmental delay.

      @ 9 months - Patient 1 sat without support.

      @ 2 years - Patient 1 began to walk.

      @ later on - Patient 1 presented with speech delay and behavioral disturbances. The behavioral disturbances were reduced by the drug risperidone.

      @ early childhood - Patient 1 has been obese and appeared to display hyperphagia.

      @ 8 years - Patient 1 developed precocious puberty.

      @ 10 years - Patient 1 had her first epileptic seizure. Treatment with lamotrigine prevented further seizures. The seizures later became refractory to this treatment.

      @ 10 years - Patient 1 presented with a height of 160 cm (P>97) and a head circumference of 59 cm (P>97). She had hypertelorism and prominent eyebrows. Her nasal bridge was broad and auricles were fleshy. In addition, Patient 1's neck was short and the fingers were short and wide.

      Patient 1 has an intellectual disability, metabolic syndrome, and macrocephaly.

      CaseNotHPOs:HP:0002540.

      CaseNotHPOFreeText:Patient 1 sat without support at 9 months old. She walked at 2 years.

      CasePreviousTesting:Patient 1 was given a conventional cytogenetic analysis. The chromosomal analyses (46,XX) and array CGH results (BlueGnome CytoChip ISCA 4×180K v1.0; Agilent Human Genome CGH Microarray 180K) were normal. Patient 1 was also given a methylation-specific MLPA to exclude a Temple syndrome, which is also characterized by weight gain and precocious puberty. In addition, Prader-Willi syndrome was ruled out by methylation testing. This syndrome is another imprinting disease causing obesity and intellectual disability.

      PreviouslyPublished:Not previously published

      GenotypingMethod:Whole-exome sequencing analysis of the entire exome was conducted for Patient 1. Whole-genome sequencing showed heterozygosity in Patient 1, which was confirmed by Sanger sequencing. Macrocephalic syndrome genes including PTEN, NSD1, NFIX, SETBP1, RAI1, and PHF6 were analyzed, and no additional variants of interest (pathogenic, likely pathogenic, or variants of uncertain significance) were observed.

      Gene:MECP2

      Variant:c. 1162_1172del (p. Pro388*), heterozygosity, frameshift.

      HGVS:NM_004992.3

      ClinVarID:Not found

      CAID:CA1139667881

      gnomAD:Not found

      MultipleGeneVariants:Not provided

    1. Patient 2

      Case:Patient 2, female, 2.5 years old

      DiseaseAssertion:CDKL5 Disorder

      FamilyInfo:De Novo with an unremarkable family history.

      ParentalGenotype:Not provided

      CasePresentingHPOs:HP:0032792, HP:0007359, HP:0011154, HP:0002194, HP:0010862, HP:0000750, HP:0012434, HP:0000710, HP:0100023, HP:0000252, HP:0009062, HP:0010845, HP:0020174, HP:0010841.

      CaseHPOFreeText:

      Patient history

      @ 2 months to 2.5 years - Patient 2 experienced tonic and focal seizures with autonomic symptoms.

      @ 2 years - Patient 2 was diagnosed with CDKL5 disorder.

      Patient 2 had severe delayed gross motor development, severe delayed fine motor development, severe delayed language development and delayed social development. Patient 2 had hyperoral and hand flapping stereotypies. She also had microcephaly (< 2 SD) and axial hypotonia.

      Patient 2 was given a brain EEG which detected Delta slowing of the background followed by generalized attenuation during seizures and multifocal interictal epileptiform abnormalities.

      @ 2.5 years - Patient 2 experienced seizures in clusters.

      Treatments

      Antiepleptic treatments included:Phenobarbital, Topirimate, Clobazam, Valproate, Keppra, Vitamin B6, Phenytoin, Lamotrigine, Nitrazepam, Oxcarbazepine, Mirtazipine, KCI, Levetiracetam, ant ketogenic diet.

      Patient 2 was also treated with Carnitine. Patient 2's seizures seem to have changed in type over time but continued.

      Other testing

      Tests were conducted on Patient 2 to obtain data in the following areas: NBS, lactate, lipoprotein profile, plasma and urine amino acids, urine organic acids, acylcarnitine profile, total and free serum carnitine levels, plasma ammonia, total plasma homocysteine, serum CK, liver enzymes, urine alpha-AASA, creatine, biotinidase, VLCFA, Batten disease screen, CSF analysis (amino acids, lactate, glucose, protein, cell count, neurotransmitters), MRI-brain with spectroscopy, and karyotype.

      CaseNotHPOs:Not provided

      CaseNotHPOFreeText:Not provided

      CasePreviousTesting:Patient 2 was given a gene sequence test for the genes SCNIA and MECP2. No gene mutation was found for these two genes.

      PreviouslyPublished:Not previously published

      GenotypingMethod:Patient 2 was given a CDKL5 gene sequence test. The test detected a mutation in the CDKL5 gene.

      Gene:CDKL5

      Variant:c. 2480_2486dupCAGATCT. frameshift

      The authors state that CDKL5 gene sequencing detected a de novo duplication in exon 17, c. 2480_2486dupCAGATCT, resulting in a frameshift (Boston University School of Medicine, Center for Human Genetics, Boston, MA). This is a novel change that has not been reported before in ExAC. Only pathogenic point mutations in exon 17 have previously been reported. In a patient with a previously reported frameshift mutation in exon 18, a truncated CDKL5 transcript was detected. A truncated protein would lack the C-terminus and would not localize correctly in the cell, as demonstrated in vitro. Accordingly, any reading frame altering mutations proximal to exon 18 are null-variants. Therefore, this mutation is classified as pathogenic according to ACMG criteria.

      HGVS:Not provided

      ClinVarID:547188

      NM_001323289.2 (CDKL5): c. 2480_2486dup (p. Gln830fs)

      Allele ID: 538303

      CAID:Not found

      gnomAD:Not found

      MultipleGeneVariants:Not provided

    1. In this study, we report a 4-month-old boy with T−B+NK− SCID due to an unreported nonsense mutation in exon 2 of the IL2RG gene. The patient was derived from a twin pregnancy, and his twin brother was asymptomatic with no immune defects. In order to confirm the pathogenic effect of the detected novel variant on the protein structure, a modeling process was performed.

      Case: Patient, Male, 4 months old <br /> DiseaseAssertion: SCID <br /> FamilyInfo: third child of non-consanguineous parents; has twin brother that is asymptomatic with no immune defects; no family history of primary immunodeficiencies <br /> CasePresentingHPOs: HP:0002014, HP:0020099, HP:0030148 <br /> CaseHPOFreeText: diarrhea, norovirus infection, heart murmur <br /> CaseNotHPOs: increased CRP <br /> CasePreviousTesting: N/A <br /> Gene: IL2RG <br /> Variant: NM_000206(IL2RG): <br /> ClinVar: <br /> CAID: <br /> gnomAD: <br /> SupplementalData:

    1. Sixty-six individuals representing 54 families were studied (Supplementary Material, Table S1). All individuals were found to harbor two ABCA4 variants likely to cause the retinal disease (18,20–26). In 40 families (74%), independent segregation of the two alleles was demonstrated. The ages at the time of their first visit ranged from 9 to 74 years (mean = 35.9, median = 35.2 years); in the majority of individuals (36/66=55%), data were available from a second visit that occurred on average 8.7 years (range=2–20 years, median = 6.9 years) after the first visit.

      Case#: Patient #35, male, 35yo at report, 14yo at onset,

      DiseaseAssertion: STGD

      FamilyInfo: family 30, segregation was noted as "yes" but no other details provided

      CasePresentingHPOs:

      CaseHPOFreeText:

      CaseNotHPOs:

      CaseNotHPOFreeText:

      GenotypingMethod:

      PreviouslyPublished: n/a

      Variant: allele 1: A1038V;L541P allele 2: G818E

      ClinVar: 99135

      CAID: CA227000

      SupplementalData: supplemental table 1

    1. STGD87 2588G→C Q1750X Yes

      Case#: STGD87, 10-14yo at onset, German

      DiseaseAssertion: STGD

      FamilyInfo: segregation in family

      CasePresentingHPOs:

      CaseHPOFreeText: "The diagnosis of STGD was based on the demonstration of bilateral impairment of central vision and the appearance of perimacular and/or peripheral yellow-white flecks, with or without atrophy of the central retinal-pigment epithelium and a normal or only mildly abnormal flash electroretinogram when recorded in early stages of the disease."

      CaseNotHPOs:

      CaseNotHPOFreeText:

      GenotypingMethod: denaturing gradient gel electrophoresis, dHPLC, and SSCP analysis, PCR amplification of individual coding exons and flanking intron sequences, direct DNA sequencing

      PreviouslyPublished: n/a

      Variant: Q1750X; 2588G→C in trans "Correct segregation of disease alleles was demonstrated in all 39 cases in which family samples were available for study"

      ClinVar: 7879

      CAID: CA119128

      SupplementalData: n/a

    2. STGD47/164 IVS13+1G→A 2588G→C Yes

      Case#: STGD47/164, 10-14yo at onset, German

      DiseaseAssertion: STGD

      FamilyInfo: segregation in family

      CasePresentingHPOs:

      CaseHPOFreeText: "The diagnosis of STGD was based on the demonstration of bilateral impairment of central vision and the appearance of perimacular and/or peripheral yellow-white flecks, with or without atrophy of the central retinal-pigment epithelium and a normal or only mildly abnormal flash electroretinogram when recorded in early stages of the disease."

      CaseNotHPOs:

      CaseNotHPOFreeText:

      GenotypingMethod: denaturing gradient gel electrophoresis, dHPLC, and SSCP analysis, PCR amplification of individual coding exons and flanking intron sequences, direct DNA sequencing

      PreviouslyPublished: n/a

      Variant: IVS13+1G→A; 2588G→C in trans "Correct segregation of disease alleles was demonstrated in all 39 cases in which family samples were available for study"

      ClinVar: 7879

      CAID: CA119128

      SupplementalData: n/a

    1. Case 4A 52-year-old male was examined for declining vision OS over the past few months. He was previously clinically diagnosed with STGD 7 years before presentation. Family history was not significant for ocular disease. Best-corrected visual acuity measured 20/100 OD and 20/70 OS. Spherical refractive error measured −3.00 OD and −3.25 OS. Anterior segment examination was unremarkable and applanation tonometry measured 17 mmHg OD and 14 mmHg OS. Posterior segment examination was significant for central atrophy and classic peripheral pisciform flecks sparing the peripapillary regions OU (Figure 4, A and B). Autofluorescence imaging demonstrated inner atrophic flecks and outer hyperautofluorescent flecks. Moderate peripapillary hypoautofluorescence, but not atrophy, was present, likely secondary to the patient’s myopia (Figure 4, C and D). Genotyping revealed two heterozygous ABCA4 mutations, P1380L and S1696N.Open in a separate windowFig. 4Case 4. STGD mutation IVS40 + 5G>A. A, Color Photo OU. B, Red-Free Photo OU reveal central atrophy and classic peripheral pisciform flecks sparing the peripapillary regions OU. C, Autofluorescence OD. D, Autofluorescence OS show that the innermost flecks are hypoautofluorescent, consistent with atrophy, whereas the outermost flecks are hyperautofluorescent, demonstrating excess lipofuscin. There is moderate peripapillary hypoautofluorescence that is not as dark as this patient’s central atrophy or the peripapillary atrophy of Case 1. This finding may thus be due to the patient’s myopia.

      Case#: Hwang Case 4, male, 52yo at report, 45yo at onset

      DiseaseAssertion: Stargardt

      FamilyInfo: Family history was not significant for ocular disease.

      CasePresentingHPOs: HP:0000545

      CaseHPOFreeText: declining vision OS, BCVA was 20/100 OD and 20/70 OS. Spherical refractive error measured −3.00 OD and −3.25 OS. Posterior segment examination was significant for central atrophy and classic peripheral pisciform flecks sparing the peripapillary regions OU (Figure 4, A and B). Autofluorescence imaging demonstrated inner atrophic flecks and outer hyperautofluorescent flecks. Moderate peripapillary hypoautofluorescence, but not atrophy, was present (Figure 4, C and D).

      CaseNotHPOs: HP:0500087

      CaseNotHPOFreeText:

      GenotypingMethod: Genotyping was performed by the ABCR400 microarray followed by direct sequencing to confirm identified variants.

      PreviouslyPublished: n/a

      Variant: P1380L and S1696N

      ClinVar: 7904

      CAID: CA129033

      SupplementalData: n/a

    2. Case 1A 55-year-old male was examined for long-standing central visual impairment since age 18. Family history was not significant for ocular disease. His best-corrected visual acuity of 20/350 OD and 20/200 OS was consistent with measurements over the last 20 years. Spherical refractive error measured −3.5 OD and −2.0 OS. Anterior segment examination was unremarkable and applanation tonometry measured 17 mmHg OD and 14 mmHg OS. Posterior segment examination and autofluorescence imaging were significant for sharply demarcated central and peripapillary zones of atrophy and the absence of fleck lesions (Figure 1, A–D). Humphrey visual fields demonstrated bilateral central scotomas with eccentric fixation at the inferior border. ERG examination was subnormal and similar to results obtained 22 years ago.Open in a separate windowFig. 1Case 1. STGD with peripapillary atrophy and mutations P1380L and IVS40 + 5G>A. A, Autofluorescence OD. B, Color Photo OD. C, Autofluorescence OS. D, Color Photo OS. All show marked peripapillary and macular atrophy with a sharply demarcated zone of sparing between them. These characteristics caused initial diagnostic confusion with choroidal sclerosis.Genetic testing was employed for further diagnostic information and two heterozygous ABCA4 mutations, P1380L and IVS40 + 5G>A, were identified and classified as disease-causing alleles, thereby confirming the diagnosis of STGD.

      Case#: Hwang Case 1, US, male, 55yo at report, 18yo at onset

      DiseaseAssertion: Stargardt disease

      FamilyInfo: Family history was not significant for ocular disease

      CasePresentingHPOs: HP:0007663, HP:0500087, HP:0000603, HP:0000512

      CaseHPOFreeText: BCVA of 20/350 OD and 20/200 OS. Spherical refractive error measured −3.5 OD and −2.0 OS. Posterior segment examination and autofluorescence imaging were significant for sharply demarcated central and peripapillary zones of atrophy and the absence of fleck lesions (Figure 1, A–D). Humphrey visual fields demonstrated bilateral central scotomas with eccentric fixation at the inferior border.

      CaseNotHPOs:

      CaseNotHPOFreeText:

      GenotypingMethod: Genotyping was performed by the ABCR400 microarray followed by direct sequencing to confirm identified variants.

      PreviouslyPublished: n/a

      Variant: P1380L and IVS40 + 5G>A

      ClinVar: 7904

      CAID: CA129033

      SupplementalData: n/a

    1. 13/8 35 0.017 163 0 – 3 – 3 4 L541P R1098C

      Case#: Patient 13, 35yo

      DiseaseAssertion: STGD

      FamilyInfo: Family 8

      CasePresentingHPOs:

      CaseHPOFreeText: Visual acuity=0.017. OCT ft (μm)=163. MP (dB)=0. Fundus=3(extensive atrophic-appearing RPE changes). ERG=3(abnormal responses involving both rods and cones). mfERG=4(subnormal mfERG in the entire test field (0°–30°) plus pathologic Ganzfeld ERG).

      CaseNotHPOs:

      CaseNotHPOFreeText:

      GenotypingMethod: PCR of coding regions, intron/exon boundaries, and 5′ and 3′ regions of ABCA4; Standard cycle-sequencing reactions with BigDye Terminator

      PreviouslyPublished:

      Variant: L541P; R1098C

      CAID: CA226911;

      SupplementalData: n/a

    1. 3 39 6/6 1 RCD Val552Ile 6/6

      Case#: Case 3, 39yo

      DiseaseAssertion: BEM, RCD

      FamilyInfo: n/a

      CasePresentingHPOs:

      CaseHPOFreeText: a ring of increased AF surrounding decreased foveal AF, visual acuity= 6/6, 6/6

      CaseNotHPOs:

      CaseNotHPOFreeText: acquired toxic aetiology

      GenotypingMethod: The entire coding sequence (50 exons), including exon–intron boundaries, of the ABCA4 gene of each patient was screened using single‐stranded conformational polymorphism (SSCP) analysis and direct sequencing.

      PreviouslyPublished: n/a

      Variant: Val552Ile heterozygous

      CAID: CA239745

      SupplementalData: n/a

    1. An uncommon case of retinitis pigmentosa patients basedon clinical and genetic studyAyudha Bahana Bahana Ilham Perdamaian, MSc2, Dewi Kartikawati Paramita, PhD3, Riris Istighfari Jenie,PhD4, Supanji Supanji, PhD11Universitas Gadjah Mada Fakultas Kedokteran Kesehatan Masyarakat dan Keperawatan, 2Doctorate Program of Health andMedicine Science, Faculty of Medicine, Public Health, and Nurse, Universitas Gadjah Mada, Yogyakarta, Indonesia. Departmentof Ophthalmology, Faculty of Medicine, Public Health and Nursing, Universitas Gadjah Mada, 3Department of Histology andMolecular Biology, Faculty of Medicine, Public Health and Nursing, Universitas Gadjah Mada, Yogyakarta, Indonesia,Integrated Research Laboratory, Faculty of Medicine, Public Health and Nursing, Universitas Gadjah Mada, Yogyakar,4Department of Pharmaceutical Chemistry, Faculty of Pharmacy, Gadjah Mada University, Yogyakarta, IndonesiaCASE REPORTThis article was accepted: 25 August 2024Corresponding Author: Supanji SupanjiEmail: supanji@ugm.ac.id19-An uncommon00304.qxp_3-PRIMARY.qxd 29/08/2024 3:47 PM Page 98

      PMID:39215425

      Gene: ABCA4

      HGNC ID: 34

      case27-year-old male, the brother of case 1

      DiseaseAssertion: Table 1 The summary of the clinical assessment of IRD patients’ family in this research fro there down they did a whole pannel on the family

      Pedigree one can be fore form the beggginnings of case presention section?

      CasePresentingHPOs: Case 2, a 27-year-old male, the brother of case 1 had blurry vision which was not corrected with an eyeglass and inconveniences under bright light starting from 14 years ago. Case 2 also underwent a fundus examination after finding that case 1 was RP. In further examination of those patients and their family members found that case 1 was confirmed as RP and case 2

      CaseHPOFreeText:NA

      CaseNotHPOs:NA

      CaseNotHPOFreeText:NA

      Genotyping Method:NA

      PreviouslyPublished:NA

      Variant:NA

      ClinVar:

      CAID:NA

      SupplementalData:NA

      Inheritance pattern Autosomal Recessive

    1. (http://genetics.bwh.harvard.edu/pph2/). In addition, mutation taster predicted both L168F and L168S variant as disease-causing with PROVEAN predictions of L168F (-2.767) and L168S (-4.083) as deleterious (https://www.mutationtaster.org/). As such, it was not surprising that the L168S variant patient had much more severe disease onset and rapid progression compared to other SCA34-causing ELOVL4 variants. For example, a patient carrying the T233M ELOVL4 variant was reported to develop ataxia starting at 15 years of age [10]. However, at the time of examination of this patient at 60 years of age, an MRI of the brain showed only subtle flattening of the ventral pons and mild cerebellar atrophy [10]. Another patient carrying the Q180P ELOVL4 variant developed ataxia in his mid-20 s and showed cerebellar and pontine atrophy [11]. Japanese patients also carrying the W256G variant developed gait ataxia between 13–56 years of age [12]. However, disease progression was reported to be very slow, and patients did not require assistance with walking with a walker or cane until the age of 60 years or older [12]. Taken together, it looks like the nature of the mutation and its effect on normal ELOVL4 function most likely through defects in VLC-FA biosynthesis or conformational changes in protein structure are critical to disease onset and severity of the pathologies.

      SupplementalData:

    1. WDR19-associated retinopathy presenting with adult-onset Stargardt-likephenotype

      PMID:39967245

      Gene: ABCA4

      HGNC ID: 34

      Case#:39 man

      DiseaseAssertion:NA

      FamilyInfo:NA

      CasePresentingHPOs:Snellen in both eye, visual impairment with night blindnessisual acuity was 20/20Snellen in both eyes, with a minor correction for astig-matism. The anterior segment and intraocular pressurewere within normal limits. On fundus examination, dif-fuse fleck-like lesions were scattered both inside and out-side the arcades, while sharply demarcated areas ofmacular atrophy with foveal sparing, more pronouncedin the left eye, were visible.

      CaseHPOFreeText:NA

      CaseNotHPOs:NA

      CaseNotHPOFreeText:

      Genotyping Method:Next-Generation Sequencing (NGS), using theTruSight One Clinical Exome sequencing panel on anIllumina NexSeq500 platform, enriching for 4800 genesincluding ABCA4, CNGB3, ELOVL4, PROM1, and PRPH2

      PreviouslyPublished:Under refernces?

      Variant:WDR19 variants:the novel deletion at c.1777 + 1 within the donor splicingsite (class 4) and the rare c.1430 G>T variant causing theamino-acid substitution p.(Arg477Leu) (class 3) in the putative protein. Additionally, a heterozygous c.1793A>G(class 3) variant in the CDH23 gene was found, though it was deemed as not contributive to the patient’s clinical phenotype. All reported variants were confirmed throughSanger sequencing

      ClinVar:NA

      CAID:NA

      SupplementalData:NA

    1. The STGD patient from Family 12 is a compound heterozygous with p.Val931Met and a novel nonsense mutation at exon 33 (p.Glu1574X; Figure 1B). Disease onset for this patient was at age 43. Ophthalmic examination revealed moderate central retinal changes, decreased mfERG responses exclusively in the central 15 degrees, and decreased visual acuity.

      Case#: Family 12 Proband, male, 43yo at onset, Portuguese

      DiseaseAssertion: Stargardt

      FamilyInfo: no affected family members in pedigree (Fig. 1)

      CasePresentingHPOs: HP:0007663

      CaseHPOFreeText: "The criteria for STGD phenotype included bilateral central vision loss and pigmentary macular lesions, normal caliber of retinal vessels, absence of pigmented bone spicules, and compatibility with recessive mode of inheritance." Moderate central retinal changes, decreased mfERG responses exclusively in the central 15 degrees

      CaseNotHPOs:

      CaseNotHPOFreeText:

      PreviouslyPublished: n/a

      Variant: p.Glu1574X; p.Val931Met. Several other polymorphisms also reported. ABCR400 gene chip microarray, DHPLC

      ClinVar: 1460063

      CAID: CA341283936

      SupplementalData: n/a

    1. Mutation scanning and direct DNA sequencing of all 50 exons of ABCR were completed for 150 families segregating recessive Stargardt disease (STGD1). ABCR variations were identified in 173 (57%) disease chromosomes, the majority of which represent missense amino acid substitutions. These ABCR variants were not found in 220 unaffected control individuals (440 chromosomes) but do cosegregate with the disease in these families with STGD1, and many occur in conserved functional domains. Missense amino acid substitutions located in the amino terminal one-third of the protein appear to be associated with earlier onset of the disease and may represent misfolding alleles. The two most common mutant alleles, G1961E and A1038V, each identified in 16 of 173 disease chromosomes, composed 18.5% of mutations identified. G1961E has been associated previously, at a statistically significant level in the heterozygous state, with age-related macular degeneration (AMD). Clinical evaluation of these 150 families with STGD1 revealed a high frequency of AMD in first- and second-degree relatives. These findings support the hypothesis that compound heterozygous ABCR mutations are responsible for STGD1 and that some heterozygous ABCR mutations may enhance susceptibility to AMD.

      Annotating here since the full text is a PDF.

      Case#: Family AR321 proband, US, 6yo at onset

      DiseaseAssertion: Stargardt

      FamilyInfo: proband and two other siblings are affected

      CasePresentingHPOs: The essential and defining features of STGD were (1) pedigrees with at least one living affected individual compatible with autosomal recessive inheritance; (2) an ophthalmoscopically characteristic retinal disorder in families with both parents living; (3) bilateral central visual loss with both “beaten metal” elliptical foveal dystrophy and temporal pallor of the optic discs, documented by retinal color photography, with or without yellow-pigment epithelial flecks in the macular and/or retinal “near periphery”; and (4) the characteristic fluorescein angiographic feature of a dark choroid (Blacharski 1988).

      CaseHPOFreeText:

      CaseNotHPOs:

      CaseNotHPOFreeText: (1) evidence of autosomal dominant inheritance; (2) any history of night blindness, loss of peripheral vision, or "retinitis pigmentosa"; (3) cataracta complicata or cells in the vitreous; (4) substantially abnormal electroretinographic or electrooculographic responses; (5) no fluorescein angiography performed or no dark choroid documented; (6) neurological disease (including loss of cognition or seizures); 7) drug exposures (especially to antimalarial and agents known to cause crystalline retinopathies); or (8) any "atypical" maculopathies in which a unique diagnosis of STGD could not be established.

      PreviouslyPublished: PMID: 8533764

      Variant: c.3113C>T p.A1038V; c.1715G>C p.R572P . Heteroduplex and SSCP analyses were used to screen the 50 exons of ABCA4. Linkage analysis and haplotype analysis were previously performed

      ClinVar: 99073

      CAID: CA226919

      SupplementalData: n/a

    1. JB260 Stargardt ABCA4 c.6119G>A p.Arg2040Gln rs148460146 Zernant et al (2014)50 c.2879del p.Ala960Aspfs*17 N/A

      Case#: Bryant Subject JB260, US

      DiseaseAssertion: Stargardt

      FamilyInfo:

      CasePresentingHPOs: "Stargardt disease is a childhood-onset macular degeneration and is most commonly caused by mutations in ABCA4. Characteristic yellow flecks are typically seen under the macula during a fundus exam."

      CaseHPOFreeText:

      CaseNotHPOs:

      CaseNotHPOFreeText:

      GenotypingMethod: WES; previously screened using arrayed primer extension (APEX) multigene panels for the relevant disease and no disease-causing variants had been identified; PCR and Sanger for verification

      PreviouslyPublished: n/a

      Variant: c.6119G>A p.Arg2040Gln; c.2879del p.Ala960Aspfs*17

      CAID: CA232815

      SupplementalData:

    1. See Supplementary Table S2 for a complete genotypic glossary of the cohort.

      Case#: Patients were identified from the inherited retinal disease (IRD) database at UC San Diego (UCSD).

      DiseaseAssertion: RP with macular edema

      FamilyInfo:

      CasePresentingHPOs:

      CaseHPOFreeText: Dx of RP based on "a history of progressive peripheral vision loss or nyctalopia, and ocular examination findings of RP including bone spicule pigmentation, disc pallor and attenuated vessels and genetic confirmation."

      CaseNotHPOs:

      CaseNotHPOFreeText:

      GenotypingMethod: Next-generation sequencing (NGS), exome sequencing, and/or targeted Sanger sequencing were the primary genetic testing approaches.

      PreviouslyPublished: PMID:10206579 is referenced but it seems a reference to the variant and not the proband

      Variant: c.6383A>G (p.His2128Arg); c.3G>T (p.Met1?). phase unknown

      ClinVar: 99455

      CAID: CA227399

      SupplementalData: Variant is found in table S2

    1. Protoderm will differentiate or mature into what type of tissue? How about procambium? And ground meristem?

      the 3 foundation system: protoderm, procambium, ground meristem

    2. Explain how roots elongate and increase in diameter via primary and secondary meristems.

      Roots grow in length through primary growth (elongation) and increase in thickness through secondary growth (widening), driven by specialized tissues called meristems.

    1. eLife Assessment

      This valuable study shows the impact of the metabolic state of bacteria on phage infection. The experimental results, based on various phages infecting E. coli, are convincing and consistent with a two-step adsorption mathematical model. This study should be of interest to the communities working on cell metabolism and on host-pathogen interactions.

    2. Reviewer #1 (Public review):

      In the wild, bacteria can be found in a wide range of metabolic states, including states in which they are resource limited. Because phages heavily rely on the infected cell's molecular machinery to replicate, it is natural to wonder how phage-bacteria interactions depend on the metabolic state of the cell. In this work, Marantos et al. investigate specifically how the rate of infection of 5 different phages changes between cells grown in energy-rich conditions and cells grown in energy-depleted conditions. Their results clearly show that 4 out of the 5 phages studied display a significant reduction in infection rate in cells that are energetically depleted and provide a potential explanation for this observation by looking into the mechanisms that these phages use to irreversibly infect their host cells.

      The work also tries to explain the observation using a mathematical/mechanistic model that describes infection as the sequence of two steps, where a phage first needs to bind to a cell receptor, from which it can potentially unbind, and then irreversibly infects by injecting its genome. The mechanistic interpretation offered by the model highlights an interesting trade-off between adsorbing to a metabolically active host and discriminating between active and inactive hosts that, somehow, a phage has to optimize. It would be interesting, in the future, to investigate how different phages optimize this task.

      Comments on revised version.

      I am happy with how the authors have addressed all the comments. The manuscript is much clearer and more readable and the previous overstated claims have been removed/clarified.

    3. Reviewer #2 (Public review):

      Summary:

      The authors investigate the dependence of phage adsorption rates on host metabolic state, using 5 coliphages that differ in their infection cycles and host receptors. They find that four of the 5 phages showed significantly reduced infection under low metabolic states, with phage that generally have weaker adsorption being more strongly affected by low metabolism. The authors complement their findings with a 2-step infection model where phages can disengage from their hosts after initial adsorption. The paper illustrates the power of standardized experimental protocols for quantitative trait comparisons and highlights the dependence of phage infection success on host physiology.

      Strengths:

      The paper is well written and clearly structured.

      The experiments are well designed and particularly commendable is the diligent use of control scenarios to allow for quantitative comparison between phages. This standardized protocol will be valuable for the entire phage community.

      The authors convincingly show the impact of host physiology on phage adsorption success. This dependence has so far mainly been considered for intracellular phage replication and the paper shows that host physiology has to be taken into account at all steps of phage infection.

    4. Reviewer #3 (Public review):

      Marantos et al. showed that for some coliphages, the energetic state of the bacterial host cell has a strong impact on whether phage infection is initiated. The authors drew this conclusion from the observation that there are more free phages remaining in the medium after infection of arsenate-azide-treated cells as compared to after infection of untreated cells. These data were analyzed and reported both as ratios of the treated vs. untreated conditions and using a mass-action kinetic model of phage-cell collision in the infection mixture. The data supported the findings that for four phages infecting Escherichia coli bacteria, namely, phages λ, 𝜙80, m13, and T6, the phages are less likely to initiate infection if the host bacteria are energy depleted. However, for phage T5, the authors found that their infection propensity is not impacted.

      As I have stated in the first submission of this manuscript, the data presented by the authors clearly supported the principal conclusion of the study. The five phages chosen by the authors represent different viral lifestyles and infection mechanisms, highlighting the potential applicability to other Escherichia coli phages. Finally, the authors successfully use a classic mass-action model of phage-cell collision to interpret their data. The simplicity of their experimental assay, combined with the use of this mathematical model, offers other investigators who study phage-bacterial interactions in other contexts a potentially useful toolkit to examine infection in general, and specifically, the dependence of phage infection on the host's metabolic state.

      Comments on revised version.

      In this revised version, the authors have successfully resolved all of my comments. I appreciate that the main text has been majorly revamped, which greatly helps the readers follow the motivation behind the experiment and analyses, and interpret the data. I agree that the revised terminology choice "commitment to infection", instead of the previous interchangeably used "adsorption"/"entry", is much more logical, considering the experimental data. I also commend the authors for writing the modeling part in a very clear, pedagogical, and instructive manner. Overall, I believe that this manuscript will be valuable to those who are interested in phage-bacterial interactions.

    5. Author response:

      The following is the authors’ response to the original reviews.

      Public Reviews:

      Reviewer #1 (Public review):

      In the wild, bacteria can be found in a wide range of metabolic states, including states in which they are resource-limited. Because phages heavily rely on the infected cell's molecular machinery to replicate, it is natural to wonder how phage-bacteria interactions depend on the metabolic state of the cell. In this work, Marantos et al. investigate specifically how the rate of infection of 5 different phages changes between cells grown in energy-rich conditions and cells grown in energy-depleted conditions. Their results clearly show that 4 out of the 5 phages studied display a significant reduction in infection rate in cells that are energetically depleted and provide a potential explanation for this observation by looking into the mechanisms that these phages use to irreversibly infect their host cells.

      The work also tries to explain the observation using a mathematical/mechanistic model that describes infection as the sequence of two steps, where a phage first needs to bind to a cell receptor, from which it can potentially unbind, and then irreversibly infects by injecting its genome. While the model is sensible from a mechanistic perspective, the experimental evidence that supports how each model's rate is affected by the cell metabolic state is weak, as only ratios of these rates can be inferred from the data.

      Reviewer #2 (Public review):

      Summary:

      The authors investigate the dependence of phage adsorption rates on host metabolic state, using 5 coliphages that differ in their infection cycles and host receptors. They find that four of the 5 phages showed significantly reduced infection under low metabolic states, with phages that generally have weaker adsorption being more strongly affected by low metabolism. The authors complement their findings with a 2-step infection model where phages can disengage from their hosts after initial adsorption. The paper illustrates the power of standardized experimental protocols for quantitative trait comparisons and highlights the dependence of phage infection success on host physiology.

      Strengths:

      The paper is well written and clearly structured.

      The experiments are well-designed, and particularly commendable is the diligent use of control scenarios to allow for quantitative comparison between phages. This standardized protocol will be valuable for the entire phage community.

      The authors convincingly show the impact of host physiology on phage adsorption success. This dependence has so far mainly been considered for intracellular phage replication, and the paper shows that host physiology has to be taken into account at all steps of phage infection.

      Weaknesses:

      There are some concerns about the experimental setup and which conclusions can be drawn from it:

      Before phage infection, bacterial cultures are grown to exponential growth, washed, and then resuspended with glucose or arsenate-azide for 10min. It is however, questionable that 10 minutes is enough to simulate high and low metabolic states realistically. 10 minutes seems to be quite short to go from exponential growth to a low metabolic state, given the transcriptional memory of previous environments. It seems more likely that the population will be quite heterogeneous, with cells in various states of transition towards low metabolic states.

      While we agree with the reviewer that during metabolic transitions there may be a period in which the population is heterogeneous, with cells in different stages of transition toward a low metabolic state, the 10-minute treatment used here was chosen based on prior work showing that arsenate–azide rapidly inhibits cellular energy metabolism and is sufficient to eliminate the hyper diffusion of the λ receptor (Winther et al., Biophysical Journal 2009, http://dx.doi.org/10.1016/j.bpj.2009.06.027). We have also corrected the DOI for this reference in the manuscript. Furthermore, the ATP pool of log-phase E. coli turns over several times per second (Holms et al., Arch. Mikrobiol. 1972, http://dx.doi.org/10.1007/BF00425016). We therefore assumed the bacteria were energy depleted after 10 minutes. We have clarified this point in the revised manuscript.

      Given that arsenate and azide inhibit cellular metabolism, i.e., have antimicrobial effects, cells might not just downregulate metabolism but also activate the stress response, and this causes some of the observed effects on phage adsorption. Therefore, the 'low metabolic state' of the cells in this paper could mean that cells are starved or that they are stressed or both.

      The reviewer is correct. We don’t exclude indirect effects. However, as nutrients were removed from the bacteria by washing and energy metabolism was inhibited by the addition of arsenate and azide, we assumed a stress response requiring biosynthesis would be unlikely to occur.

      The abundance of receptors could change between the high and low metabolic media conditions and contribute to the observed differences in adsorption, while the authors seem to assume in their model that the initial adsorption rate always remains the same.

      We do not think that the observed differences in adsorption are explained by a change in receptor abundance. In a previous study using the same experimental protocol as in the present work, phage λ was compared to the metabolically insensitive mutant λh (Brown et al., PNAS 2022, http://dx.doi.org/10.1073/pnas.2106005119). If the lower adsorption in the low-metabolic condition were caused by a reduced number of receptors, then λh should also have shown a lower adsorption rate under the same condition. Instead, no measurable effect on λh adsorption rate was observed. We therefore conclude that the effect is not explained by changes in receptor number on the timescale of the experiment. We have clarified this point in the revised manuscript.

      Reviewer #3 (Public review):

      Summary:

      Marantos et al. showed that for some coliphages, the energetic state of the bacterial host cell has a strong impact on whether phage infection is initiated. The authors drew this conclusion from the observation that there are more free phages remaining in the medium after infection of arsenate-azide-treated cells as compared to after infection of untreated cells. These data were analyzed and reported both as ratios of the treated vs. untreated conditions and using a mass-action kinetic model of phage-cell collision in the infection mixture. The data supported the findings that for four phages infecting Escherichia coli bacteria, namely, phages λ, ɸ80, m13, and T6, the phages are less likely to initiate infection if the host bacteria are energy-depleted. However, for phage T5, the authors found that their infection propensity is not impacted.

      Strengths:

      The data presented by the authors clearly supported the principal conclusion of the study ("Viral commitment to infection depends on host metabolism"). The five phages chosen by the authors represent different viral lifestyles and infection mechanisms, highlighting the potential applicability to other Escherichia coli phages. Finally, the authors successfully used a classic mass-action model of phage-cell collision to interpret their data. The simplicity of their experimental assay, combined with the use of this mathematical model, offers other investigators who study phage-bacterial interactions in other contexts a potentially useful toolkit to examine infection in general, and specifically, the dependence of phage infection on the host's metabolic state.

      Weaknesses:

      (1) The authors isolated and measured the numbers of free phages in the medium after infection of bacteria under different treatments. These measurements were analyzed in two different ways: (1) simply as ratios (corrected/normalized using different controls), and (2) fitted using a simple mathematical model. I have concerns regarding both analyses.

      (1.1) For the first method, having different time points at which the sample of each phage is collected critically complicates data interpretation. As one incubates the phage-bacteria mixture for a longer time, more infection occurs, and the number of phages collected from the mixture decreases. Therefore, the different incubation time forfeits the goal of "a systematic and quantitative comparison across different phages [...]", just as the authors self-criticized. Conceivably, the authors could have used the shortest measurement time for all phages (i.e., 10 minutes, as for phage λ). Alternatively, the authors could have applied a systematic criterion such as half (or any other fraction) of the latent period of each phage, which would still "maximize the incubation period while ensuring that manipulations were completed before the first infection cycle concluded". In my view, the seemingly arbitrary measurement time for each phage renders the entire first analysis very challenging to interpret. It also goes against the author's proposition that the protocol was "standardized" or "consistent". It is not clear what the readers are supposed to take away from this first analysis, or rather, which evidence, finding, or conclusion the manuscript would lose if the authors only presented the modeling-based analysis.

      (1.2) The second method of analysis sought to remove the dependence of the measurements on time. I completely agree with this goal, and the findings extracted from this analysis significantly contributed to the merits of this manuscript. However, the authors achieved this goal using a single time point for each phage to calculate the infection rate (η). As shown in Figure S3, each of the phage depletion curves is anchored by only one data point (note that the P(t)/P(0) = 1 at t = 0 is assumed, not measured). This goes against the typical way this collision model is used in the literature, where a time series is measured and used to fit the model (e.g., DOI 10.1007/978-1-60327-164-6 18, or more recently, PMID 39700139). This practice in the current manuscript reduced the robustness of the inferred η values. This problem is exacerbated by assumptions used by the authors in formulating this model. For instance, the authors used a constant value for the bacterial concentration, B, because "bacterial growth and lysis were negligible" (lines 135-136). However, considering that the bacteria were cultured at 37oC in a very rich medium (first in YT broth, then in 2% glucose), the measurement times of 20, 30, and 55 minutes are most likely one or a few generations of bacterial growth and division.

      Related note: I suggest that one of the panels in Figure S3 should be moved to the main text, since it is critical to the second method of analysis.

      We would like to clarify that the manuscript does not present two separate methods, but rather one method presented in two steps: a first step with results that are directly tied to the experimental measurements and show whether the effect is present for each phage, followed by a second, analytical step that makes the results comparable across phages.

      The first step presents the ratios because they directly reflect the measurements performed in the experiment and allow the reader to see the effect of the metabolic state for each phage in contrast to its control. We agree that these ratios are time-dependent and therefore not suitable for quantitative comparison between phages. Their purpose is to illustrate the experimental outcome and to show that the effect is present (or absent) on a per-phage basis not to compare magnitudes across phages.

      We then follow this with the second step, allowing the reader to follow the logic of the analysis. The analytical step that follows does not represent a second method, but a continuation of the same analysis. Here, we remove the time-dependence specifically in order to make comparison of the effect across phages possible, by connecting our results to standard measures such as the adsorption rate η. Importantly, P(0) is measured for every phage in every experiment. The only modeling assumption used (a standard one in the field) is the exponential form for the decay in free phage number, which naturally yields P(t)/P(0) = 1 at t = 0.

      Regarding the reviewer’s concern that bacterial growth may not have been negligible over the relevant time window, we note that recent work on rich-to-minimal growth lags in E. coli reports substantial delays before growth resumes after nutrient downshift. One 2023 study (Wu et al., Nature Microbiology 2023, https://doi.org/10.1038/s41564-022-01310-w) considering wild-type E. coli shows in Fig. 2c a lag of up to about 2 hours after a shift from MOPS minimal medium with 0.2% glucose plus 18 amino acids to the same medium without amino acids. Another 2023 study (Zhu and Dai, Nature Communications 2023, https://doi.org/10.1038/s41467-023-36254-0) examining both rel+ and rel− strains reports a growth lag of about 49 minutes for rel+ and more than 5 hours for the relA deletion strain. While these conditions are not identical to ours, they support the general point that growth does not immediately resume after such shifts. We therefore think it is unlikely that, following transfer from YT, the cells underwent one or a few full generations during the time window of our adsorption measurements.

      On the related note: Following the comments of all reviewers on Figure S3, we have decided to remove it to avoid confusion.

      (2) The data were able to distinguish phages that successfully infected bacteria and those that remained free in the medium, and the authors appropriately interpreted the data as such throughout the Results section. However, in the Discussion (starting from the very first sentence, line 172), the authors used terms that include "adsorption" and "entry" more interchangeably (for example, see the three sentences in lines 310-313, for "viral entry efficiency is shaped by [...]", then "adsorption kinetics modeling"). I do not see how the authors' data could distinguish between adsorption (the phage particles attaching to the outside of the cell) and entry (the phage DNA being injected into the cell). Conceivably, any phage particles that irreversibly attach to a cell but do not yet inject their genome into the cell would still be removed from the medium and therefore not quantified. Another example: in lines 189-191, the authors interpreted that "[...] when the bacterium is in a low metabolic state, the phage does not bind irreversibly to the host", but how do the authors eliminate the case of no phage binding (i.e., the reversible step) to begin with?

      We agree with the reviewer that our use of the terms adsorption, entry, and infection should have been more careful. Our experiment can only identify the irreversible commitment of phage to a host cell. We have therefore revised the text to refer consistently to phage commitment.

      Similarly, in lines 283-293, how do the authors delineate whether energy depletion would increase the k_off term or decrease the k_inj term, because either would result in more free phages in the medium as observed in the data? I believe that the writing of the Discussion, as it stands now, is doing a disservice to the conclusions presented in the Results section.

      We thank the reviewer for this important point. We agree that the model would work either by k_off or k_inj being dependent on the host metabolic state, and that our original wording was therefore too restrictive. The data do not distinguish between these possibilities; they only constrain the ratio k_off/k_inj. In the revised text, we therefore formulate the argument in terms of this ratio: if energy depletion leads to reduced commitment, this can arise either because k_off increases, because k_inj decreases, or because both change, as long as k_off/k_inj becomes larger in the inactive case. Put differently, what matters is not which individual rate changes, but that the balance between leaving and committing shifts in a way that disfavors commitment to inactive cells. This also leads to the trade-off now discussed in the revised manuscript: efficient commitment to active hosts requires a small k_off/k_inj, whereas strong discrimination against inactive hosts requires this ratio to become significantly larger in the inactive case. Depending on whether this is achieved through changes in k_off or k_inj, the cost of discrimination appears either as slower commitment or as additional energy dissipation. We agree that the previous wording overstated the mechanistic interpretation, and we have revised the Discussion accordingly to bring it in line with what the Results actually support. Based on the comments from all reviewers, we have also revised the terminology throughout the manuscript: instead of error correction, we now refer to this as a discrimination process, and we replaced k_inj by k_com to reflect that our assay resolves irreversible phage commitment rather than DNA injection specifically.

      (3) The authors presented an argument that performing infection of all five phages in the same condition is an advantage, allowing for comparison across different phages. While this goal is a completely valid one, it is difficult to reconcile that with the fact that different phages require different optimal conditions for successful infection. For instance, phage T5 famously requires Ca2+ for successful infection into the host bacterium (and later successful replication); see PMID 13174489. However, all infections were performed in TMG, which lacks Ca2+. Perhaps the absence of T5 dependence on the host metabolism is because the infection condition used by the authors was not optimal for T5 to begin with? Similar arguments could be made for other phages.

      Our study alone cannot eliminate that possibility. However, we have cited multiple previous studies, for example references citing Braun et al., showing that T5 remains insensitive to the host metabolic state under different buffer conditions. We therefore believe it is unlikely that the lack of metabolic dependence we observe for T5 is simply due to suboptimal infection conditions.

      (4) Whereas the manuscript examined five coliphages, only phage T5 and phage λ were discussed extensively. I believe some discussion points for these two phages need clarification.

      We focused our discussion on the phages T5, λ and φ80 because these are the phages for which similar effects have been reported previously in the literature. This allowed us to connect our findings directly to existing work and to discuss mechanistic hypotheses in a meaningful comparative framework. For the remaining phages, to our knowledge no prior studies have examined their behavior under comparable metabolic conditions, and therefore a similarly detailed discussion would have been speculative. Nevertheless, all five phages are treated equally in the presentation of the experimental results and in the quantitative comparison of adsorption rates.

      (4.1) Phage T5: The data obtained by the authors show that the infection rate of phage T5 is not impacted by the metabolic state of the host cell. Considering that the authors used the terms "infection", "adsorption", and "entry" interchangeably to refer to the irreversible commitment of a phage to a host cell (see point 2), this discussion regarding phage T5 lacks one critical literature context: DNA entry of phage T5 is known to occur in two phases (first-step transfer and second-step transfer). Critically, the second step can only occur if phage proteins encoded by the phage DNA transferred in the first step are expressed (see PMID 10577483 and the cited papers therein). In that context, metabolic poisoning of the host bacteria should have impeded T5 infection. The authors should comment on this point.

      As the reviewer pointed out, our usage of the terms infection, adsorption, and entry should have been more careful. Our experiment can only identify irreversible commitment of phage to a host cell. For T5, we expect that this irreversible commitment already occurs upon first-step transfer of phage DNA. As a result, even if second-step transfer is impeded under metabolic poisoning, our method would not resolve that effect. We have added this clarification to the revised manuscript.

      (4.2) Phage λ: The experiment using phage λ in this current study shares many resemblances to that in Brown et al. 2022. That feature alone is not a problem, but at many places in the text, the writing is ambiguous as to whether it is discussing the results in Brown et al. 2022 or in the current manuscript. I am giving three examples below, but this is not exhaustive: (i) Lines 67-69, there is no Brown et al. 2022 reference immediately after "a mutant phage variant (λh) could bypass this dependency [...]" (not just in the previous sentence); (ii) Line 228 should clearly say "Our previous findings suggested that phage λ is capable of [...]", since it concerns Brown et al., 2022, not the current study; and (iii) Lines 245-246, there is no Brown et al., 2022 reference immediately after "we observed that a mutant variant [...] even energy-depleted host" (without a reference, it reads like the authors "observed" that finding in this current manuscript).

      The reviewer is right. In those places, the text was ambiguous as to whether it referred to the present study or to Brown et al. (2022). We have now inserted the reference at the relevant points and revised the wording where needed to make this distinction explicit.

      Also, regarding phage λ: The discussion between line 230 and line 249 is very interesting, but since it concerns the differences between λ PaPa and Ur-λ, the authors should consider mentioning and discussing a very relevant recent study, PMCID: PMC6312755.

      We agree that the study by Guan et al. is very relevant and interesting. However, our point in this part of the Discussion is only to clarify that we used λ PaPa and not the originally isolated λ strain. We have therefore limited the discussion here to that distinction.

      (5) Control experiments, or references to prior studies, are needed to support that the As/Az treatment at this concentration and duration (at least 10 minutes) is sufficient to deplete the metabolic state of the cell. For instance, this can be shown by impeded or null cell growth, arrested motility (using a standard swimming assay), or a fluorescent reporter for the energetic state of the cell.

      The 10-minute treatment used here was chosen based on prior work showing that arsenate–azide rapidly inhibits cellular energy metabolism and is sufficient to eliminate the hyperdiffusion of the λ receptor (Winther et al., Biophysical Journal 2009, http://dx.doi.org/10.1016/j.bpj.2009.06.027) where the effect was assessed by monitoring the rate of movement of the λ receptor on the bacterial surface. We have clarified this point in the revised manuscript.

      Recommendations for the authors:

      Reviewer #1 (Recommendations for the authors):

      As mentioned earlier, I found the paper interesting and addressed an important and significant knowledge gap.

      My biggest concern is about the interpretation of the experimental data in light of the two-step model. In particular, around line 286, it is stated "k_inj is more sensitive to metabolic state than k_off". Assuming k does not depend on metabolic state, which is a fair assumption, the equation for eta only depends on the ratio between k_inj and k_off and not on the individual parameters separately. Consequently, there is no way of saying which one of the two is more affected by metabolic state, unless the model already assumes that k_off is not influenced by metabolic state. The results could equally be explained by k_inj decreasing in metabolically depleted cells, or k_off increasing in such cells. If this is an assumption of the model, this should be clearly stated and not reported as a consequence of the data, as it is at the moment. Also, how does this mathematical model connect to the fitting function used in Figure 2b?

      We thank the reviewer for this important point. We agree that the model would work either by k_off or k_inj being dependent on the host metabolic state, and that our original wording was therefore too restrictive. The data do not distinguish between these possibilities; they only constrain the ratio k_off/k_inj. In the revised text, we therefore formulate the argument in terms of this ratio: discrimination requires that k_off/k_inj be larger for inactive hosts than for active hosts, such that commitment is specifically reduced in the inactive case. Put differently, what matters is not which individual rate changes, but that the balance between leaving and committing shifts in a way that disfavors commitment to inactive cells. This introduces a trade-off: efficient commitment to active hosts requires a small k_off/k_inj, whereas strong discrimination requires this ratio to become significantly larger for inactive hosts. If this is achieved through changes in k_off, discrimination comes at the cost of slower commitment by allowing more time to leave; if it is achieved through changes in k_inj, it can preserve fast commitment to active hosts but requires additional energy dissipation in order to actively modulate commitment. We have therefore revised the text accordingly to frame the argument in terms of this trade-off, rather than attributing the effect specifically to k_inj. Based on the comments from all reviewers, we have also revised the terminology throughout the manuscript: instead of error correction, we now refer to this as a discrimination process, and we replaced k_inj by k_com to reflect that our assay resolves irreversible phage commitment rather than DNA injection specifically.

      I have a related experimental criticism. The kinetic model presented assumes an exponential decay of free phage, which is a commonly used assumption in the phage literature. Given that the phage types used in this study lyse relatively slowly, it would be good to actually see adsorption curves, in which free phage is measured at different time points between inoculation and lysis. This data would not only provide useful evidence for the kinetic model, but it should also replace what is now in Figure S3, which consists of fitting one experimental point with one line. As it currently stands, Figure S3 is not useful actually misleading.

      We appreciate the reviewer’s point. We agree that adsorption curves, in which free phage is measured at different time points between inoculation and lysis, would provide a stronger basis for evaluating the kinetic model. However, we do not have the resources to perform these additional experiments within the scope of the present study. Following the comments of all reviewers on this point, we have therefore decided to remove Figure S3 to avoid confusion.

      Finally, it is not clear to me why the quantity "Ratio" has been chosen to be presented in Figure 1, rather than the ratio of estimated adsorption rates eta'/eta, which is much more intuitive for a phage study and contains the same information. I would recommend switching to this choice, unless there is a clear rationale for why the quantity "Ratio" is more useful/effective. Showing eta'/eta would also increase the readability of Figure 1, as it would move the y-axis to a logarithmic scale and better visualize values around 1.

      We used “Ratio” in Figure 1 to illustrate the experimental design, controls, and measured quantities directly, as it more transparently reflects the data collected. In the second part of the analysis, where we compare time-independent adsorption rate estimates, we have presented the corresponding values of η′/η as suggested.

      Minor comments:

      (1) Introduction

      Line 31: "... such as nutrient limitation, fluctuating temperatures, and variable energy availability" - if drawing a distinction between energy availability and nutrient limitation, please make explicit what this distinction is. Energy availability seems like a natural consequence of nutrient availability.

      While energy and nutrient availability are often linked in E. coli, they represent distinct physiological constraints. Nutrient limitation refers to the lack of essential biosynthetic precursors such as nitrogen, phosphorus, or amino acids. Energy availability, in contrast, reflects the cell’s ability to generate ATP and reducing equivalents through metabolic processes. For example, under anaerobic conditions, E. coli may have ample nutrients but limited energy production due to the lower efficiency of fermentation compared to aerobic respiration. Thus, energy limitation can occur independently of nutrient limitation.

      (2) Results

      (a) Whole Section: Please label equations.

      All equations have now been labelled in the revised manuscript.

      (b) Lines 105 to 114: As stated in Major Comments, I think the clarity of the paper would be improved by introducing the relative adsorption rate here and dropping the concept of Ratio entirely. However, if the authors wish to use Ratio, I would recommend the following:

      Lines 105 to 109 are confusing to read because of the number of connectives: "... ratio of free viruses from permissive AND resistant hosts respectively TO the free viruses in buffer under energy-depleted AND energy competent conditions". This would be clearer if each quantity were given an algebraic symbol, and RP, RR, and Ratio were defined through formal algebra, rather than mixed mathematical and sentence notation.

      This section has been rewritten for clarity. We now introduce explicit algebraic symbols and define the quantities formally, which removes the ambiguity present in the sentence-only description while retaining the intended meaning.

      The chemical names "arsenate" and "azide" should appear in the body of the text before they appear abbreviated in an equation. Please state at this point that these are both metabolic inhibitors, as it is not immediately clear what role they play or why you are using them.

      The text has been updated to introduce arsenate and azide by name before the abbreviations are used, and we now explicitly note that they act as metabolic inhibitors.

      On line 114, the authors helpfully provide an interpretation of Ratio = 1. It would be useful to provide at the same time interpretations of Ratio >1 and <1, perhaps 2 and 0.5 specifically?

      We have added brief explanations illustrating the interpretation of Ratio values greater than and less than 1, including examples of 2 and 0.5.

      I would consider giving this quantity a more interpretable name than Ratio. This quantity represents how much a bacteriophage preferentially adsorbs to metabolically active cells, so perhaps "Selectivity" or "Adsorption Bias"?

      We intentionally retained the generic term “Ratio”, as this quantity reflects an intermediate experimental measure used to describe the process rather than a newly defined metric. Its purpose is to bridge the experimental observations and the subsequent quantification of effects on the adsorption rate (η).

      (c) Lines 117 to 122: the authors sometimes refer to ratios explicitly, "average ratio of around 1.6" and other times say e.g., "a greater than 3 times increase in viral particles". Using more consistent language (saying "Ratio" every time) would be clearer.

      We have standardized the terminology in this section and now refer to all fold-changes consistently using “Ratio” to avoid ambiguity.

      (d) Figure 1

      Phages λ and T6 look like they have ratios less than 1 for resistant cells? If this is true / if the ratio is statistically significantly below 1, please comment.

      The ratios for λ and T6 are not statistically different from 1. The apparent deviation is within the standard error of the mean. To make this clearer, we have added the corresponding p-values to Table S2 in the Supplementary Information.

      Ratios near 1 are difficult to distinguish from 1, especially in panels A and D. Using a logarithmic scale on the y-axis would make the plots more readable.

      Because the values in these panels are not statistically different from 1, changing to a logarithmic scale would not alter the interpretation. We therefore retained the current axis scaling to reflect that there is no meaningful deviation from 1 in these cases.

      The data corresponding to individual experiments have no error bars. Given that the number of free virions was determined by plaque assay, which carries an intrinsic sampling error, this uncertainty should be reflected in the plots.

      We thank the reviewer for this important comment. Because plaque assays have compound sources of stochastic variation, assigning a per-measurement error bar would risk implying false precision. For this reason, we present the values from each biological replicate directly, and the uncertainty is represented in the statistical summary across replicates. Specifically, for each phage and condition we show the three independent experimental measurements and report the mean along with the standard error of the mean. This approach allows us to represent biological variability without implying a precision that cannot be accurately quantified at the level of single plaque counts.

      Similarly, the average value does show error bars, but it is not stated what these error bars correspond to: standard error in the mean, standard deviation of the sample, or combined uncertainty?

      The caption has been updated to state that the error bars represent the standard error of the mean.

      The resistant bacteria seemed to have ratios close to 1 in all cases. Is this because very few virions adsorbed under both energy conditions?

      Resistance is commonly associated with a lack of a surface receptor for the phage (or generally an entry pathway). We use the resistant bacteria as a control group for the effect of the conditions on adsorption. For resistant bacteria, the Ratio should be 1 since virions do not adsorb under both energy conditions. Any slight variations from 1 should come from sampling errors or small heterogeneity in the population.

      (e) Figure 2

      Please comment on what the error bars here represent. Error bars in Figure 2 A seem to permit negative (or at least zero) values of relative adsorption rate for phages m13 and T6, possibly implying an overestimate of the error? If it is the case that multiple values used to calculate the mean are far apart, possibly showing the values individually through a superimposed swarm plot would be clearer.

      This point is now addressed in the Supplementary Information, where we clarify how the error bars were calculated.

      (3) Discussion

      (a) Line 189: "high metabolic state" is imprecise. Say "energy-competent" to be consistent with earlier language.

      To maintain continuity with earlier terminology, we now include “energy-competent” in parentheses alongside “high metabolic state,” while retaining the original phrasing for readability.

      (b) Figure 3, population level

      Show adsorbed virions physically attached to bacteria, rather than removing them completely from the image, as currently, the implication is that at a high metabolic state, there are fewer virions total, not fewer virions remaining in solution because more are adsorbed. You could go as far as to add a third "after centrifuging" row, showing the adsorbed phages stuck in the pellet and the unadsorbed phages remaining in solution.

      Thank you for this suggestion. Figure 3 has been updated to depict adsorbed virions attached to bacterial cells, clarifying that the decrease represents adsorption rather than loss of total particles. This change improves the accuracy and interpretability of the schematic.

      (4) Methods and Materials

      (a) Figure 5

      The step "estimate cell numbers from OD" appears to follow incubating plates overnight. If the cells you are counting come from the pellet produced by centrifuging 3 steps prior, you could add a fork into the black line connecting the steps, with one branch corresponding to the supernatant and phages, and the other to the pellet and cells?

      Thank you for pointing this out. The order in the figure has been corrected: cell numbers are estimated from OD before overnight incubation. This resolves the confusion without the need for branching in the workflow diagram.

      (a) Line 332

      You allow as much time as possible for adsorption without the possibility of lysis. Did you determine the lysis times / latent periods of these phages through one-step-growth-curves, or use published results, in which case please cite? Having obtained the lysis time by either method, what fraction of the lysis time did you allow for adsorption? Also, please add supplementary tables with lysis times used for the different phages.

      We thank the reviewer for this comment. We used published latent-period values as guides and verified compatibility with our own system when selecting incubation times. We have clarified this in the text and added the relevant citations. We did not use a common fixed fraction of the lysis time for all phages; instead, incubation times were chosen to allow sufficient time for adsorption but not for completion of the first lytic cycle. For λ, productive lytic development was blocked in the host background used, as in Brown et al., PNAS 2022, http://dx.doi.org/10.1073/pnas.2106005119. For ϕ80 and T5, we used published latent-period values as guides and verified their compatibility with our own system (De Paepe and Taddei, PLoS Biology 2006, http://dx.doi.org/10.1371/journal.pbio.0040193). M13 is a chronic filamentous phage and therefore does not have a standard lytic latent period; in our host–phage combination, it required more than 1 h before phage release. For T6, we relied primarily on the kinetics observed in our own system, since adsorption was unusually slow for this phage–host pair under our assay conditions. Although literature reports describe shorter T6 latent periods under specific assay conditions (Foster and Johnson, Journal of General Physiology 1951, http://dx.doi.org/10.1085/jgp.34.5.529), this is consistent with published work showing that adsorption and infection kinetics can vary substantially with host background, surface structure, and experimental conditions (Heller and Braun, Journal of Bacteriology 1979, http://dx.doi.org/10.1128/jb.139.1.32-38.1979; Storms et al., Biochemical Engineering Journal 2012, http://dx.doi.org/10.1016/j.bej.2012.02.010).

      (5) Supplementary

      Figure S1

      This data is useful in understanding the main body of the paper, and I think this should form part of a main figure (possibly with the individual experimental data points superimposed over the bars). This could come before or as part of Figure 1?

      We thank the reviewer for this suggestion. We have explored including these data directly in the main figure but found that doing so substantially reduced the readability of the figure, as the underlying table is visually dense. For this reason, we chose to summarize the results in Figure 1 and present the detailed data separately in Figure S1 of the Supplementary Material, along with the Ratio analysis, which more effectively conveys the trends without overloading the main figure.

      Reviewer #2 (Recommendations for the authors):

      Minor comments:

      (1) L16-18: This sentence could be made more accessible as 'error correction' is not an intuitive term in the phage field.

      We have updated the overall theory section including the terminology. Instead of error correction, we now refer to it as a discrimination process.

      (2) L96-98: Does this potentially indicate a trade-off where evolution for stronger binding cannot evolve at the same time as responsiveness to metabolic activity?

      We agree that this sentence made a stronger evolutionary claim than our data support. Since we only tested four laboratory phages, we cannot conclude that there is an evolutionary trade-off between stronger binding and responsiveness to host metabolic activity. We have therefore removed this sentence to avoid making an unsupported evolutionary interpretation.

      (3) L102: What does 'post-cellular' mean?

      Postcellular supernatant is simply the liquid that remains after cells have been removed. During centrifugation, the cells pellet at the bottom, and the liquid above (which can contain viruses) is the postcellular supernatant.

      (4) L105-107: Worth splitting into two sentences as it is a bit unclear if ratios are built between permissible and resistant hosts or between buffers or both.

      Thank you for the suggestion. We have rewritten this section into two sentences to clarify how the ratios are constructed, and we hope the revised wording improves readability.

      (5) L110-122: Figures S1 and S2 could be referenced here.

      References to Figures S1 and S2 have now been added in this section.

      (6) L137: As P(0) is the viral concentration in buffer, I am assuming that the phage lysate has been diluted in buffer and phages have been added to cultures from the same dilution tube to guarantee equal starting numbers, but I couldn't find this in the methods.

      This clarification has been added to the Methods and Media section of the Supplementary Information.

      (7) L243: It would be worth defining what 'hyperdiffusion' means.

      We have added a brief definition of “hyperdiffusion”.

      (8) L253-256: I do not entirely follow this explanation.

      We thank the referee for pointing out this lack of clarity. This was also raised by Reviewer #3. The point we intended to convey is that λ behaves differently toward E. coli LamB depending on whether it is on a living cell or isolated in buffer, but makes no such distinction for Shigella LamB, binding it in both contexts. More specifically, previous work showed that wild-type E. coli extracts could only inactivate λ in the presence of added solvents, whereas control extracts prepared similarly from Shigella did not require added solvent for λ inactivation. This observation is consistent with E. coli LamB requiring a specific state to irreversibly bind λ. We therefore meant to suggest that the capacity for metabolic-state sensing is not simply a function of phage identity, but also depends on receptor-specific properties that differ between the two bacterial species.

      We have rephrased it as follows: Notably, wild-type λ is inactivated by E. coli K-12 extracts only when solvents are added, whereas Shigella extracts inactivate λ without this requirement (Randall-Hazelbauer and Schwartz, J. Bacteriol. 1973; Schwartz, J. Mol. Biol. 1975; Schwartz and Le Minor, J. Virol. 1975). This suggests that E. coli LamB requires a specific state for irreversible binding, a conditionality absent in Shigella LamB, indicating that the capacity for metabolic-state sensing may depend on receptor-specific properties.

      (9) L284: Why is k_inj necessarily more sensitive to the metabolic state than k_off? Could membrane changes under stress increase k_off?

      We thank the reviewer for this important point. We agree that the model would work either by k_off or k_inj being dependent on the host metabolic state, and that our original wording was therefore too restrictive. The data do not distinguish between these possibilities; they only constrain the ratio k_off/k_inj. In the revised text, we therefore formulate the argument in terms of this ratio: reduced commitment in inactive cells can arise through an increase in k_off, a decrease in k_inj, or both, as long as k_off/k_inj becomes larger in the inactive case. What matters is therefore not which individual rate changes, but that the balance between leaving and committing shifts in a way that disfavors commitment to inactive cells. This also underlies the trade-off now discussed in the manuscript: efficient commitment to active hosts requires a small k_off/k_inj, whereas strong discrimination against inactive hosts requires this ratio to become much larger in the inactive case. We have revised the Discussion accordingly to bring it in line with what the Results actually support. Based on the comments from all reviewers, we have also revised the terminology throughout the manuscript: instead of error correction, we now refer to this as a discrimination process, and we replaced k_inj by k_com to reflect that our assay resolves irreversible phage commitment rather than DNA injection specifically.

      (10) Figure 1: There seems to be more variation between replicates in phage Lambda than in other phages. Is this caused by receptor number heterogeneity in the population?

      Unfortunately we do not have a way to compare receptor number heterogeneity across the different phage receptors in our experiments. We therefore cannot conclude that the larger variation observed for phage λ is caused by receptor number heterogeneity in the population.

      (11) Figure S1: There seems to be a significant difference between phage Lambda viability in the two buffers - do the authors have an idea where this comes from?

      There is no difference in λ viability between the two buffers. The apparent difference in the figure is due to sampling variability.

      (12) Figure S3: Last sentence of the legend probably shouldn't say 'upper'.

      Following the suggestions from all of the reviewers we have removed Figure S3 as it created more confusion than clarity.

      Reviewer #3 (Recommendations for the authors):

      (1) The text reads as incomplete in some places. Can the authors please provide clarifications on the following points?

      (1.1) Lines 235-256: How do the authors draw a conclusion that "a phage can detect host metabolic status" from a study that used purified LamB receptors (i.e., no live cells with any metabolism) extracted from two different bacterial species (i.e., not a difference in metabolic states)?

      We thank the referee for pointing out this lack of clarity. This was also raised by Reviewer #2. The point we intended to convey is that λ behaves differently toward E. coli LamB depending on whether it is on a living cell or isolated in buffer, but makes no such distinction for Shigella LamB, binding it in both contexts. More specifically, previous work showed that wild-type E. coli extracts could only inactivate λ in the presence of added solvents, whereas control extracts prepared similarly from Shigella did not require added solvent for λ inactivation. This observation is consistent with E. coli LamB requiring a specific state to irreversibly bind λ. We therefore meant to suggest that the capacity for metabolic-state sensing is not simply a function of phage identity, but also depends on receptor-specific properties that differ between the two bacterial species.

      We have rephrased it as follows: Notably, wild-type λ is inactivated by E. coli K-12 extracts only when solvents are added, whereas Shigella extracts inactivate λ without this requirement (Randall-Hazelbauer and Schwartz, J. Bacteriol. 1973; Schwartz, J. Mol. Biol. 1975; Schwartz and Le Minor, J. Virol. 1975). This suggests that E. coli LamB requires a specific state for irreversible binding, a conditionality absent in Shigella LamB, indicating that the capacity for metabolic-state sensing may depend on receptor-specific properties.

      (1.2) Line 270, in the abstract, and in the caption of Figure 4: The authors described the model using terms such as "an error-correction mechanism" or "standard error correction", but there is little explanation. Can the authors clarify what kind of "error" is discussed here, and how it is "corrected"? In the "standard error correction" model, what determines which method of correction is "standard"? If "error correction" is a standard term in phage-bacterial interaction modeling, please provide references.

      We agree with the reviewer that our use of the term error correction was not appropriate in this context. The proper term is discrimination process rather than error correction. We have now corrected this terminology throughout the manuscript and clarified the underlying logic in the relevant sections.

      (1.3) Line 301: The authors speculated that phage T5 is "better suited to ecological niches", but I am not sure how that is consistent with their data showing T5 is more rampant, that they infect both energy-competent and energy-depleted cells, not just depleted cells. Why "niches", and why are T5 better suited to environments "where energy-limited cells dominate", not just any environment?

      We agree that this point was not stated clearly enough. What we intended to convey is that T5 would be at a net disadvantage in a niche containing a mixture of energy-competent and energy-deficient hosts. We have updated the main text accordingly.

      (1.4) Line 303, and related to point 6.3. above: Phage λ can also infect and replicate in "starved bacterial cells" (shown in Kourilsky 1974 and Geng et al. 2024, both of which were cited in this manuscript). How do the authors reconcile these reports with the discussion point in line 303, and their data that only phage T5, but not λ, shows insensitivity to the host metabolic state?

      Our data do not imply that phage λ is unable to infect starved bacteria. As shown in Kourilsky (1974) and Geng et al. (2024), λ can indeed infect and replicate in nutrient-limited cells. Our results specifically indicate that λ infection under starvation proceeds with a reduced adsorption rate, while T5 maintains the same adsorption rate even when the host is starved. Thus, our conclusion is that T5 is insensitive to the host metabolic state at the level of adsorption, whereas λ is not. We acknowledge that the wording in line 303 may have unintentionally led to confusion, and we have revised this part of the text to avoid that.

      (2) The following comments relate to the text and figures in the manuscript. There are many places in the manuscript that could use fine proofreading and copy-editing for clarity and consistency. For example:

      (2.1) If I understand it correctly, the equation in between lines 109 and 110 should be clarified using terms such as "Free viral particles after mixing with bacteria in Arsenate and Azide" and "Free viral particles in bacteria-free buffer with Arsenate and Azide". As it stands, it is not clear which terms correspond to conditions where bacteria are present.

      The equation has been updated to explicitly indicate which terms refer to mixtures containing bacteria and which refer to bacteria-free controls, so that the correspondence between conditions is now clear.

      (2.2) Equations in between line 276 and 283, and elsewhere: Some concentration terms are enclosed in brackets ("[BP]"), while most are not.

      This notation has been clarified. We now use “[PB]” specifically to denote the transient phage–bacterium complex, distinguishing it from the product P⋅B. All other concentration terms are written without brackets for consistency.

      (2.3) Figure 4 and in equations: "BP" or "PB"?

      The notation has been made consistent throughout; we now use “PB” exclusively to denote the phage–bacterium complex.

      (2.4) Line 284 and line 286: The "inj" in "k_inj" is sometimes italicized, sometimes not.

      The notation has been standardized so that k_inj is now formatted consistently throughout the manuscript, without italicizing “inj.” Also we have replaced k_inj by k_com to reflect that our assay resolves irreversible phage commitment rather than DNA injection specifically.

      (2.5) Figure 5: Was the step "Estimate cell numbers from OD" really performed on the next day after the experiment (i.e., >12 hours after infection and phage plating), not immediately after cell washing?

      Thank you for pointing this out. The figure has been updated to reflect the correct order of steps: cell numbers are estimated from OD immediately after washing, followed by overnight incubation of the plates.

      (2.6) Figure S1: As it stands now, the x-axis of each panel can be read either as "Permissive, Resistant bacteria, Buffer" (missing "bacteria" for the first pair of bars), or "Permissive (bacteria), Resistant (bacteria), Buffer (bacteria)" (extra "bacteria" for the last pair of bars).

      The intended interpretation is the second one (permissive bacteria, resistant bacteria, buffer).

      (2.7) Figure S3: The panel letters "A" and "B" are missing in the figure. Also, it is not clear why the legend for the five phages and the legend for the measurement times are not combined.

      Following the suggestions from all of the reviewers we have removed Figure S3 as it created more confusion than clarity.

      (2.8) Strain table in the Methods and Materials: Please write genotypes with italicization, and consistently indicate mutations and deletions with the minus sign superscript or the Δ prefix. Also, for the S3222 strain: Is it really the entire Mal regulon mutated ("Mal-"), or just lamB-? In Brown et al. 2022, it was only the latter.

      Genotypes have been reformatted with consistent notation. For S3222, the correct designation is Mal-, as in the SI of Brown et al. 2022. In this case, Mal- is intended as a phenotypic designation rather than a specific genotype, and we have therefore formatted it accordingly, i.e. neither italicized nor written in lower case.

    1. the sequence of the entire human genome was completed in 2003, just 50 years after the description of the double-helical nature of DNA by Watson and Crick
      • by 2003 the sequence of the entire human genome was completed
      • this was only 50 years after the discovery that a human gene was a double helix
    2. late 1990s in sequencing the human genome led in the mid-2000s to the announcement that over 90% of the genome had been sequenced
      • from the early 1990s to the mid 2000s, much of the genome had been figured out
    3. Apart from infectious organisms and environmental pollutants, many diseases are manifestations of abnormalities in genes, proteins, chemical reactions, or biochemical processes, each of which can adversely affect one or more critical biochemical functions
      • along side of having a biochemical bases affected by the environment, much of diseases often times manifest from gene abnormalities
    4. Recent increasing emphasis on systematic attempts to maintain health and forestall disease, or preventive medicine, includes nutritional approaches to the prevention of diseases such as atherosclerosis and cancer.
      • recent approaches to maintain health, and prevent the use of medicine includes nutritional approaches in orders to prevent chronic diseases
    5. The maintenance of health requires optimal dietary intake of vitamins, certain amino acids and fatty acids, various minerals, and water.
      • the components of what is needed to maintain health
    6. The World Health Organization (WHO) defines health as a state of “complete physical, mental, and social well-being and not merely the absence of disease and infirmity.”
      • definition of health by WHO
    1. Dicots typically have a pith while monocots do not. Why?

      because their vascular bundles are scattered throughout the stem, consuming the central region with connective tissue

    2. Herbaceous perennials die back to the ground in the spring. Where do they get the energy to grow the next spring, and from what tissue do new shoots emerge?

      They store energy for the following spring as carbohydrates (starches) within specialized underground structures, and new shoots emerge from dormant renewal buds

    1. eLife Assessment

      This study makes a valuable contribution by broadening the range of eukaryotic model systems and establishing Blastocystis, the most prevalent microeukaryote in the human gut, tractable for reverse-genetics investigations. The presented imaging data are convincing and informative, although confirmation by molecular methods would further strengthen the study. The work should interest readers studying host-microbe interactions in the human gut, as well as those developing new systems for eukaryotic research.

    2. Reviewer #1 (Public review):

      Summary:

      This paper presents a toolkit for the transformation of Blastocystis. The authors have screened a number of selectable agents, promoters and reporter genes and present their findings. This resource will be of immense use to those in Blastocystsis field, as well as those seeking to establish transformation tools in other species where such tools do not yet exist. Establishing new transformation tools is extremely challenging, and the authors have done an excellent job.

      Strengths:

      The authors have carried out a systematic screen of promoters, reporter genes and selectable agents. They have screened numerous for each, and all the data is presented. It is good to see when things did not work as well as when things did - so this data set is extremely useful indeed.

      Weaknesses:

      The findings are reported by reporter gene assay (microscopy). No evidence is given using genetics. The authors claim that the DNA is maintained episomally. However, could it be possible that there is integration? No PCRS/RT-PCRs are shown (although it can safely be assumed that the DNA/RNA is present where the transformation was successful), nor are any Western blots. These would have been useful to show that the P2A ribosomal skipping had occurred, and that proteins were expressed individually rather than as a polyprotein.

      Comments on revised version.

      The authors have revised their manuscript to clarify that molecular analyses have not yet occurred and have resolved the technical/publication issues with the figures. I look forward to seeing these tools used in future publications to answer important questions in Blastocystsis research.

    3. Reviewer #3 (Public review):

      Summary:

      The primary objective of this study was to establish a practical and functional framework for propagation of stable transgenic cell lines of Blastocystis, a common animal gut microeukaryote. Although the work focused on Blastocystis ST7-B, a subtype with relatively low prevalence in humans, this choice is justified by its association with more frequent negative health effects. Beyond their relevance to the medical field, the methodological advances described here have the potential to also expand cell biology studies of this anaerobic organism, including its unusual mitochondria and redox metabolism.

      Strengths:

      Prior to this work, genetic tools for Blastocystis were very limited, relying on a single strong promoter-terminator combination. The authors successfully expanded the available promoter set across a range of expression strengths by testing two dozen variants in luciferase-based assays. Critically, they developed an integrated workflow from a modular transgenic construct design to an expanded inventory of molecular components (promoters, reporters), optimized DNA delivery, stepwise antibiotic resistance-mediated clonal selection and propagation, and to reporter validation. The evaluation of several anaerobiosis-compatible labeling strategies for live (and fixed) cell optical imaging will be particularly useful, with the SNAP-tag system appearing especially promising for Blastocystis.

      Weaknesses:

      The presented data generally provide a solid support for the conclusions that the work reached, but clarification of reasoning and several inconsistencies, as well as amendments to visual presentation of the data would be highly beneficial, as detailed below.

      (1) Episomal persistence of the construct:

      The manuscript repeatedly assumes, including in its title, that constructs persist in Blastocystis in their episomal form, but no direct evidence is provided. Although this interpretation is plausible, it should be identified more clearly as provisional. Nuclear genomic integration (e.g., via NHEJ) remains a possible explanation unless supporting evidence or rationale is provided to exclude it. Testing whether the phenotype persists without drug-mediated selection in the generated transgenic cell lines would help strengthen the case for episomal maintenance.

      (2) Promoters and terminators:

      (2.1) There is a discrepancy between the claimed number of loci (14), from which promoters used to drive luciferase expression were derived, and those detailed as having been actually generated in Table 1 (11). This inconsistency should be corrected or explained, as it creates uncertainty around the accuracy of the dataset.

      (2.2) Based on the presented evidence, constructs benchmarked in bioluminescence assays differed only in their promoter composition. Although terminator selection is mentioned in the Methods section, no additional details are provided; for instance, Table 1 and Figure 2 only list 23 promoters in total. Figure 2A likewise shows only promoter-dependent variation. If the terminator was held constant (LeguP1?), this should be stated explicitly. The authors may then consider revising the wording of having tested "23 promoter-terminator pairs" to better reflect that only promoters varied.

      (2.3) Promoter benchmarking was done with a plasmid lacking a selection marker, so it is unclear how the maintenance of the luciferase construct was ensured. Without selection, the observed reporter intensity could reflect differential or stochastic plasmid retention rather than promoter strength alone. The luminescence assay was performed 16-18 hours after transfection, but the rationale for this particular timeframe should be explained. In this context, the authors should explicitly state whether the experiments shown in Fig.2A represent biological triplicates or technical triplicates from a single transfection.

      (3) Figure 2:

      (3.1) Several aspects of the current design may lead to ambiguity for the reader. The boxplots are colour-coded, but it is unclear whether the colours carry meaning or are purely decorative. Because the data are already spatially separated into bins, additional random colouring is redundant and may suggest distinctions that are not intended. In addition, the part A of Figure 2 is split into two panels with the scale for the left panel shown in the right panel and some of the boxplot colours falling in the range of the scale, but not in line with their counterparts in the left panel. Because the colour use is not consistent, it is difficult to tell whether the same scale should be applied to both panels or how it should be interpreted.

      (3.2) The left panel of the part A uses a diverging blue-white-red colour scheme, which is most appropriate when the midpoint represents a meaningful central value such as zero. Because the values shown in this graph are only positive, a non-diverging 2-colour scale or a colour palette such as 'viridis' would make the plot easier to interpret.

      (3.3) A black background should be avoided: 'B' and 'C' labels are invisible and it draws attention to a distracting design feature rather to the data themselves.

      (4) Figure 3:

      (4.1) Individual snapshots should be separated more clearly, either by using a white background or by adding visible borders to make the overall composition clearer. As currently displayed, some boundaries between fluorescent channels resemble image artifacts rather than intentional panel divisions.

      (4.2) In the parts B-D, the legend should explain more clearly what each image shows and the figure itself would benefit from annotations. There seem to be three sub-panels in each 'condition' of part B (as well as C and D): while the middle and rightmost panel can be easily inferred to represent the fluorescent protein and bright-field image, what the leftmost panels represent is not specified. If DAPI was used to dye DNA, an explanation why mostly multiple labelled regions are visible should be provided.

      (4.3) Cell morphology and appearance differ markedly between UnaG/smURFP and SNAP-tag images, which should be explained. A microscope issue is mentioned in the main text, but if that was the cause, the authors should consider replacing the images as the current distortions complicate interpretation.

      Comments on revised version.

      The revised version provides sufficient clarity and appropriate visual presentation. Some confusion evidently arose due to my misunderstanding, so I thank the authors for their comprehensive clarifications and patience.

    4. Author response:

      The following is the authors’ response to the original reviews.

      In revising the manuscript, we have focused on three main priorities raised during review: (1) improving precision around evidential claims, particularly concerning vector maintenance and P2A-mediated protein separation; (2) substantially improving figure quality, accessibility, and legend clarity; and (3) correcting inconsistencies and expanding methodological detail where requested.

      This study was intended as a foundational genetic toolkit and methodological framework for Blastocystis ST7-B, establishing practical workflows for DNA delivery, endogenous regulatory-element benchmarking, antibiotic-selected recovery, clonal propagation, and reporter-based analysis in a genetically challenging anaerobic microbial eukaryote. The central evidence presented is therefore functional in nature: reproducible transgene delivery, selectable recovery and propagation of colony-derived transgenic lines, and detectable reporter expression using multiple anaerobic-compatible reporter systems.

      We agree with the reviewers that several additional experiments, including Western blot analysis of P2A-containing constructs, outward-facing PCR, plasmid rescue assays, and selection-withdrawal experiments, would further strengthen the mechanistic interpretation of the system and help distinguish episomal persistence from genomic integration. We have therefore revised the manuscript throughout to clearly separate what is directly demonstrated from what remains a plausible working interpretation or important future direction.

      Importantly, the revised manuscript no longer presents episomal maintenance or complete P2A-mediated protein separation as demonstrated conclusions. Instead, these are now discussed explicitly as unresolved mechanistic questions requiring future molecular analysis. Nevertheless, the central methodological conclusion remains unchanged: stable selectable transgene expression, recovery of colony-derived transgenic lines, and reporter-positive Blastocystis ST7-B transformants can now be reproducibly obtained.

      Reviewer #1 (Public review):

      Summary:

      This paper presents a toolkit for the transformation of Blastocystis. The authors have screened a number of selectable agents, promoters and reporter genes and present their findings. This resource will be of immense use to those in the Blastocystis field, as well as those seeking to establish transformation tools in other species where such tools do not yet exist. Establishing new transformation tools is extremely challenging, and the authors have done an excellent job.

      Strengths:

      The authors have carried out a systematic screen of promoters, reporter genes and selectable agents. They have screened numerous for each, and all the data is presented. It is good to see when things did not work as well as when things did, so this data set is extremely useful indeed.

      Weaknesses:

      The findings are reported by reporter gene assay (microscopy). No evidence is given using genetics. The authors claim that the DNA is maintained episomally. However, could it be possible that there is integration? No PCRS/RT-PCRs are shown (although it can safely be assumed that the DNA/RNA is present where the transformation was successful), nor are any Western blots. These would have been useful to show that the P2A ribosomal skipping had occurred, and that proteins were expressed individually rather than as a polyprotein.

      We thank the reviewer for the positive assessment of the manuscript and for recognising both the technical difficulty and broader utility of establishing genetic tools in Blastocystis and other experimentally challenging microbial eukaryotes. We also appreciate the reviewer’s identification of the main evidential limitations in the original manuscript, particularly regarding vector maintenance and P2A-mediated protein separation.

      First, regarding the question of vector topology and the interpretation of episomal maintenance.

      We agree that the original manuscript presented episomal persistence too strongly relative to the evidence currently available. We have therefore revised the manuscript throughout to clarify that episomal maintenance should presently be regarded as a plausible working model rather than a directly demonstrated conclusion.

      The transfection system used here was adapted from Li et al. (2019), including use of the pXS2-P<sub>Legumain</sub>-derived plasmid framework. Importantly, the construct used in the present study does not contain the original Trypanosoma brucei tubulin-targeting region associated with homologous integration in the original pXS2 system. Complete plasmid sequencing confirmed that the constructs function here as heterologous expression plasmids carrying Blastocystis ST7-B regulatory elements and transgenes. While this does not demonstrate episomal persistence, it also means that genomic integration cannot be inferred from the historical pXS2 vector architecture alone.

      We further note that comparative genomic analyses by Gentekaki et al. (2017) suggest that Blastocystis lacks components of the canonical non-homologous end-joining (NHEJ) machinery, implying that homologous recombination is likely to represent the principal route for double-stranded DNA repair. Because the constructs used here did not contain Blastocystis homology arms, there is currently no obvious mechanism favouring targeted homologous integration. Nevertheless, we fully agree that genomic integration cannot presently be excluded.

      To reflect this appropriately, the revised manuscript now explicitly separates the demonstrated functional outcomes from unresolved mechanistic questions concerning vector maintenance. We also identify several future approaches that would help distinguish episomal persistence from genomic integration, including outward-facing PCR, plasmid rescue followed by full plasmid sequencing, Southern blotting, FISH, selection-withdrawal experiments, and long-read sequencing approaches.

      We have revised the manuscript throughout to remove statements implying demonstrated episomal maintenance and now present episomal persistence only as a plausible working interpretation.

      In the Methods section under Cloning, the following text has been added:

      Lines 202–206: “The constructs used in this study were derived from the pXS2-P<sub>Legumain</sub> vector described by Li et al. (2019), which adapted a heterologous expression-vector backbone for transient plasmid-based expression in Blastocystis ST7-B. Here, the same molecular backbone was used as a plasmid scaffold carrying Blastocystis-derived regulatory elements and transgenes.”

      In the Discussion, the following text has been added/edited:

      Lines 665–673: “The molecular maintenance state of the introduced constructs remains unresolved: episomal maintenance is a plausible working model, but genomic integration cannot be formally excluded. The constructs used here lack Blastocystis homology arms, and comparative genomic analyses suggest that Blastocystis lacks canonical non-homologous end-joining components (Gentekaki et al., 2017), making targeted integration by standard repair routes unlikely but not impossible. Direct assays such as outward-facing PCR, plasmid rescue followed by full plasmid sequencing, FISH, or selection-withdrawal experiments will be required to distinguish episomal persistence from integration.”

      Second, regarding P2A-mediated protein separation.

      We agree that Western blotting would provide the most direct biochemical assessment of P2A-mediated ribosomal skipping efficiency in Blastocystis ST7-B and would help determine the extent of any residual uncleaved fusion product. We have therefore revised the manuscript to avoid implying that complete protein-level separation was directly demonstrated.

      The revised manuscript now states only what is directly supported by the current data: that P2A-containing bicistronic constructs supported antibiotic-selected recovery of transgenic lines together with detectable downstream reporter expression. The microscopy data therefore support functional downstream reporter expression, but do not by themselves exclude residual uncleaved fusion products.

      We selected P2A because it is a compact and well-characterised peptide with high reported separation efficiency across multiple eukaryotic systems, including microbial eukaryotes. However, we agree that P2A performance can be context-dependent, and we now explicitly identify biochemical validation of P2A cleavage efficiency as an important future direction.

      Importantly, these revisions do not alter the central methodological conclusion of the study, namely that selectable transgene expression, propagation of reporter-positive lines, and recovery of colony-derived Blastocystis ST7-B transformants can now be reproducibly achieved.

      Text inserted in the Results:

      Lines 394–396: “The P2A peptide is expected to promote ribosomal skipping during translation, allowing two separate polypeptides to be produced from a single open reading frame.”

      Lines 403–404: “However, protein-level separation was not directly tested, and the extent of any residual uncleaved fusion product remains unresolved.”

      Text inserted in the Discussion:

      Lines 619–629: “P2A was selected because it is a well-characterised peptide with high reported separation efficiency in human cell lines, zebrafish embryos, and mice (Kim et al., 2011). It also has precedent across microbial eukaryotes, including the protest Dictyostelium discoideum (Zhu et al., 2023), the fungi Aspergillus niger (Schuetze and Meyer, 2017) and Ustilago maydis (Müntjes et al., 2020), and the apicomplexan parasites Toxoplasma gondii (Markus et al., 2019) and Plasmodium falciparum (Dans et al., 2024). However, P2A performance is context-dependent, and the evidence presented here is functional rather than biochemical. P2A-containing constructs support antibiotic-selected recovery and downstream reporter expression in Blastocystis ST7-B, but ribosomal skipping efficiency and any residual uncleaved product will require direct protein-level validation.”

      Reviewer #1 (Recommendations for the authors):

      (1) Please could you show a Western blot to confirm if P2A has worked? It could be that the proteins are being expressed as a polyprotein.

      We agree that Western blotting would provide the most direct biochemical assessment of P2A-mediated ribosomal skipping efficiency in Blastocystis ST7-B and would help determine the extent of any residual uncleaved fusion product. This is an important point, and we have revised the manuscript accordingly to avoid implying that complete protein-level separation was directly demonstrated.

      The current study was designed as a first-generation functional genetic toolkit for Blastocystis ST7-B, focused primarily on establishing reproducible workflows for selectable transgene expression, reporter recovery, and propagation of transgenic lines in this experimentally challenging anaerobic microbial eukaryote. The toolkit is therefore validated here through functional outcomes, including antibiotic-selected survival, stable propagation through extended passaging (>15 passages) and cryopreservation, and detectable reporter fluorescence above wild-type autofluorescence.

      P2A was selected because it is a compact and well-characterised peptide with high reported ribosomal skipping efficiency across multiple eukaryotic systems, including microbial eukaryotes, as discussed above. Nevertheless, we fully agree that direct biochemical validation would strengthen the mechanistic interpretation of the bicistronic system in Blastocystis ST7-B. We therefore now explicitly identify Western blot analysis, ideally using epitope-tagged upstream and downstream products, as an important future direction for quantitative assessment of P2A cleavage efficiency and any residual uncleaved fusion products.

      Relevant manuscript revisions are described above under the general response to Reviewer 1.

      (2) Something has gone wrong with figure formatting. Figure 2 is nearly illegible and I cannot read the text in section A. Sections B, C, and D have lost their labels and are fuzzy and surrounded by black. A similar issue affects Figure 3. Everything is just black with a few cells. It is illegible when printed.

      We thank the reviewer for highlighting these presentation issues and agree that the submitted figure quality significantly impaired readability and interpretation. The problems appear to have arisen primarily during manuscript compilation and export, particularly affecting image resolution, contrast, and panel labelling in the review PDF.

      To address this, Figures 2 and 3 have been completely reformatted and replaced with revised high-resolution versions. We have also improved typography, panel separation, colour scaling, and legend clarity throughout. In response to additional reviewer suggestions, individual data points have now been added to Figures 2B and 2C to improve transparency and interpretability of the underlying data distributions.

      Figures 2 and 3 have been replaced with fully revised high-resolution versions with improved panel labelling, accessibility, typography, and figure legends.

      (3) The data from Figure 2B would be better placed in Table 1 with a column for robust/moderate/intermediate/weak/very weak. This would be much easier for the reader.

      We thank the reviewer for this helpful suggestion. We believe the comment refers to the promoter activity data shown in Figure 2A rather than the voltage optimisation data in Figure 2B. To improve readability and accessibility of these data, we have revised Figure 2A extensively to make the promoter activity tiers more legible and easier to interpret directly from the heat map and accompanying box plots.

      We considered incorporating simplified activity classifications into Table 1. However, activity patterns were construct-specific rather than simply locus-specific. In several cases, multiple promoter fragments derived from the same locus produced substantially different reporter outputs, and activity did not scale monotonically with promoter fragment length. We therefore felt that assigning a single categorical activity label at the locus level would oversimplify the dataset and reduce the construct-level resolution that is central to the toolkit value of the study.

      Instead, we addressed the reviewer’s concern by substantially improving the presentation and readability of Figure 2A, allowing readers to identify robust, moderate, intermediate, weak, and very weak expression constructs more directly while preserving the underlying construct-specific information.

      Figure 2A has been revised to improve clarity, accessibility, and legibility of the promoter activity tiers, allowing construct-level expression classes to be interpreted more directly from the heat map and accompanying boxplots.

      (4) How do you know if the constructs are maintained as episomes? Have you done an outward-facing PCR?

      We agree that direct molecular evidence distinguishing episomal persistence from genomic integration is currently lacking, and we appreciate the reviewer highlighting this important limitation. We have therefore revised the manuscript throughout to avoid presenting episomal maintenance as a demonstrated conclusion and now describe it only as a plausible working interpretation based on the current evidence and vector design.

      We have not performed outward-facing PCR in the present study. As discussed in the general response above, we now explicitly identify outward-facing PCR, plasmid rescue followed by full plasmid sequencing, selection-withdrawal assays, FISH, and long-read sequencing approaches as important future directions for resolving the molecular maintenance state of the constructs.

      The revised manuscript now clearly separates the demonstrated functional outcomes, including selectable transgene expression, recovery of colony-derived transgenic lines, and stable reporter-positive propagation under selection, from the unresolved mechanistic question of vector topology.

      This issue has been addressed throughout the revised manuscript, including in the Methods and Discussion sections, where episomal maintenance is now presented as a plausible but unconfirmed interpretation rather than a demonstrated conclusion.

      Minor Comments

      Line 66: is this one to two billion individuals with Blastocystis, or one to two billion Blastocystis cells per gut?

      The intended meaning was colonised individuals globally. We agree that the original phrasing was ambiguous and have corrected it for clarity.

      Lines 66–67 revised to: “…microorganisms in the human gut, and is estimated to colonise approximately one to two billion people globally (Scanlan and Stensvold, 2013).”

      Line 148: Supplier of IMDM?

      The supplier information was already present in the original manuscript as IMDM L0191 (Biowest).

      No additional manuscript change required.

      Line 157: Who annotated the dataset, the 2017 paper or the present study?

      The dataset annotation derives from Armengaud et al. (2017). We agree that the original wording was unclear and have revised this section substantially to improve clarity regarding the rationale and workflow used for promoter and terminator candidate selection.

      “The relevant Methods section has been extensively revised for clarity and expanded detail” (Lines 156–189).

      Line 166: Who predicted the 3′ UTR, the 2017 paper?

      This information derives from the NCBI annotation associated with the Blastocystis ST7-B genome based on Denoeud et al. (2011). This has now been clarified in the Methods section.

      Clarified in revised Methods section.

      Line 237: How long did it take in days?

      Approximately 2 days.

      Line 270 revised to: “…turned yellow without drug treatment, usually within 2 days post-transfection.”

      Line 325: Typo, missing gap between Figure and 1A.

      Corrected in revised manuscript.

      Reviewer #2 (Public review):

      This manuscript presents a substantial technical advance for the genetic manipulation of Blastocystis by establishing an integrated workflow for stable episomal transgenesis, antibiotic selection, clonal recovery, and reporter-based imaging in the ST7-B subtype. The study is particularly valuable because it combines multiple previously fragmented approaches into a coherent and practically applicable toolkit, including endogenous regulatory elements, optimized electroporation conditions, selectable markers, and anaerobic compatible fluorescent reporters. This methodological work greatly expands the molecular toolbox and future studies focused on both basic and infection biology can now build on the ability to express and localize proteins in fixed as well as live cells.

      The microscopy data are convincing and clearly demonstrate functional reporter expression and successful recovery of stable transgenic lines. Nevertheless, because this is primarily a methodological paper, the study would be further strengthened by the inclusion of Western blot validation of reporter expression and bicistronic constructs. In particular, biochemical analysis of the P2A-containing constructs would help assess the efficiency of ribosomal skipping and exclude the possible presence of uncleaved fusion proteins, thereby providing stronger support for the interpretation of the imaging data and the functionality of the expression system.

      We thank the reviewer for this thoughtful and positive assessment of the manuscript and for recognising the value of integrating previously fragmented approaches into a coherent and practically usable genetic toolkit for Blastocystis ST7-B. We particularly appreciate the reviewer’s recognition that the system expands the currently available molecular toolbox for both cell biological and infection-related studies in this experimentally challenging anaerobic microbial eukaryote.

      We also appreciate the reviewer’s comments regarding biochemical validation of the P2A-containing bicistronic constructs. We agree that Western blot analysis would strengthen the mechanistic interpretation of the reporter system by directly assessing ribosomal skipping efficiency and the possible presence of residual uncleaved fusion products. In response, we have revised the manuscript throughout to ensure that the conclusions remain appropriately evidence-based and do not imply that complete protein-level separation was directly demonstrated.

      The revised manuscript now explicitly distinguishes the demonstrated functional outcomes, including selectable transgene expression, stable propagation of reporter-positive lines, and detectable downstream reporter expression, from unresolved mechanistic questions concerning P2A cleavage efficiency and vector maintenance state. We now also identify biochemical validation of P2A-mediated protein separation as an important future direction for further refinement of the system.

      Relevant manuscript revisions addressing these points are described above under the response to Reviewer 1.

      Reviewer #2 (Recommendations for the authors):

      The quality of images could be better. The figures lacked resolution — possibly a conversion artefact.

      We agree that the figure quality in the submitted review PDF significantly reduced readability and visual interpretation. The issues appear to have arisen primarily during manuscript compilation and export, particularly affecting image resolution, typography, panel labelling, and contrast rendering.

      To address this, Figures 2 and 3 have been completely reformatted and replaced with revised high-resolution versions. We have also improved panel separation, typography, colour scaling, contrast settings, and figure legends to improve accessibility and interpretability both on screen and in print. In addition, the export workflow and file formatting have been updated to improve compatibility with journal production requirements and reduce the likelihood of compression-related rendering artefacts during manuscript compilation.

      Figures 2 and 3 have been replaced with revised high-resolution versions with improved typography, panel labelling, contrast settings, and accessibility.

      Reviewer #3 (Public review):

      Summary:

      The primary objective of this study was to establish a practical and functional framework for the propagation of stable transgenic cell lines of Blastocystis, a common animal gut microeukaryote. Although the work focused on Blastocystis ST7-B, a subtype with relatively low prevalence in humans, this choice is justified by its association with more frequent negative health effects. Beyond their relevance to the medical field, the methodological advances described here have the potential to also expand cell biology studies of this anaerobic organism, including its unusual mitochondria and redox metabolism.

      Strengths:

      Prior to this work, genetic tools for Blastocystis were very limited, relying on a single strong promoter-terminator combination. The authors successfully expanded the available promoter set across a range of expression strengths by testing two dozen variants in luciferase-based assays. Critically, they developed an integrated workflow from a modular transgenic construct design, to an expanded inventory of molecular components (promoters, reporters), optimized DNA delivery, stepwise antibiotic resistance-mediated clonal selection and propagation, and to reporter validation. The evaluation of several anaerobiosis-compatible labeling strategies for live (and fixed) cell optical imaging will be particularly useful, with the SNAP-tag system appearing especially promising for Blastocystis.

      Weaknesses:

      The presented data generally provide solid support for the conclusions that the work reached, but clarification of reasoning and several inconsistencies, as well as amendments to the visual presentation of the data, would be highly beneficial, as detailed below.

      (1) Episomal persistence of the construct:

      The manuscript repeatedly assumes, including in its title, that constructs persist in Blastocystis in their episomal form, but no direct evidence is provided. Although this interpretation is plausible, it should be identified more clearly as provisional. Nuclear genomic integration (e.g., via NHEJ) remains a possible explanation unless supporting evidence or rationale is provided to exclude it. Testing whether the phenotype persists without drug-mediated selection in the generated transgenic cell lines would help strengthen the case for episomal maintenance.

      We thank the reviewer for this important point and agree that the original manuscript presented episomal persistence too strongly relative to the currently available evidence. In particular, we agree that the title and several sections of the manuscript implied a level of mechanistic certainty that was not directly demonstrated.

      We have therefore revised the manuscript throughout to clarify that episomal maintenance should presently be regarded as a plausible working interpretation rather than a demonstrated conclusion. The revised text now explicitly distinguishes the demonstrated functional outcomes, including selectable transgene expression, recovery and propagation of colony-derived transgenic lines, and stable reporter-positive maintenance under selection, from the unresolved mechanistic question of vector topology.

      As discussed in our response to Reviewer 1, the constructs used here do not contain Blastocystis homology arms, and comparative genomic analyses suggest that Blastocystis lacks canonical non-homologous end-joining components, making targeted integration by standard repair routes less strongly supported mechanistically, although genomic integration cannot presently be excluded.

      We agree that selection-withdrawal experiments would provide useful additional evidence regarding construct persistence and have now explicitly identified such assays, together with outward-facing PCR, plasmid rescue, FISH, and long-read sequencing approaches, as important future directions for resolving the molecular maintenance state of the transgenes.

      The manuscript has been revised throughout to remove wording implying demonstrated episomal maintenance. Episomal persistence is now discussed only as a plausible working interpretation pending direct molecular validation.

      (2) Promoters and terminators:

      (2.1) There is a discrepancy between the claimed number of loci (14), from which promoters used to drive luciferase expression were derived, and those detailed as having been actually generated in Table 1 (11). This inconsistency should be corrected or explained, as it creates uncertainty around the accuracy of the dataset.

      We thank the reviewer for this careful reading and for identifying this inconsistency. We agree that the distinction between candidate loci and successfully generated promoter constructs was not sufficiently clear in the original manuscript and could create uncertainty regarding the dataset.

      The original candidate set comprised 14 loci selected for promoter and terminator discovery. However, only 11 loci yielded successfully cloned and experimentally tested promoter constructs. The remaining three loci were retained in Table 1 for completeness and transparency, as repeated cloning attempts were unsuccessful despite two independent efforts.

      We have revised the manuscript to make this distinction explicit and to clarify that the reported NanoLuc benchmarking experiments were ultimately performed using constructs derived from 11 successfully cloned loci.

      Lines 361–364: “To expand the available regulatory parts, we screened 23 NanoLuc reporter constructs containing putative endogenous promoter–terminator pairs from 11 of 14 candidate loci; three loci could not be cloned after two independent attempts and are indicated in Table 1.”

      (2.2) Based on the presented evidence, constructs benchmarked in bioluminescence assays differed only in their promoter composition. Although terminator selection is mentioned in the Methods section, no additional details are provided; for instance, Table 1 and Figure 2 only list 23 promoters in total. Figure 2A likewise shows only promoter-dependent variation. If the terminator was held constant (LeguP1?), this should be stated explicitly. The authors may then consider revising the wording of having tested “23 promoter-terminator pairs” to better reflect that only promoters varied.

      We thank the reviewer for the opportunity to clarify this point. We agree that the original presentation may have created the impression that promoter and terminator regions were independently varied and benchmarked, whereas the experimental design was primarily focused on construct-level comparison of endogenous regulatory modules.

      As described in the Methods, each construct contained a candidate endogenous upstream promoter region together with the corresponding endogenous downstream terminator region derived from the same locus. For consistency and to keep the cloning and screening strategy experimentally tractable, a fixed 500 bp downstream terminator fragment was used for each locus rather than systematically varying terminator length or independently testing terminator activity.

      We therefore retain the description “endogenous promoter–terminator pairs,” since each construct contains both endogenous upstream and downstream regulatory regions from the same genomic locus. However, we agree that the assay was not designed to independently dissect promoter versus terminator contributions to reporter output. We have revised the manuscript accordingly to make this distinction explicit and avoid ambiguity regarding the scope of the benchmarking analysis.

      Lines 365–368: “Each construct paired a candidate upstream promoter region with the corresponding downstream terminator region from the same locus, defined here as the native 500 bp sequence immediately downstream of the stop codon. Where multiple promoter lengths were tested for the same locus, the terminator fragment was kept constant (Table 1; Figure 1A).”

      This design allowed construct-level benchmarking of paired promoter–terminator modules but did not test promoter strength or terminator activity independently.

      (2.3) Promoter benchmarking was done with a plasmid lacking a selection marker, so it is unclear how the maintenance of the luciferase construct was ensured. Without selection, the observed reporter intensity could reflect differential or stochastic plasmid retention rather than promoter strength alone. The luminescence assay was performed 16-18 hours after transfection, but the rationale for this particular timeframe should be explained. In this context, the authors should explicitly state whether the experiments shown in Fig.2A represent biological triplicates or technical triplicates from a single transfection.

      We thank the reviewer for these important methodological points. We agree that the original manuscript did not sufficiently clarify the transient nature of the NanoLuc benchmarking assay or the rationale underlying the assay design and timing.

      The promoter benchmarking assay was designed as an early transient-expression screen adapted from the NanoLuc-based workflow of Li et al. (2019), with modifications, rather than as a stable-maintenance assay. No selectable marker was included because the objective was to compare relative early reporter output across constructs shortly after DNA delivery, before prolonged culture effects became dominant.

      The 16–18 h post-electroporation time point was selected based on the NanoLuc expression kinetics reported by Li et al. (2019) and empirical optimisation during assay development. This window allowed robust transient reporter detection while limiting confounding effects arising from prolonged plasmid loss, differential outgrowth, variable recovery, or later culture-level changes.

      We agree that, in the absence of selection, the observed NanoLuc signal cannot be interpreted as an absolute measure of promoter strength independent of DNA uptake efficiency, early plasmid retention, or post-transfection recovery dynamics. We have therefore revised the manuscript to clarify that Figure 2A reports relative transient reporter output under standardized early post-transfection conditions rather than isolated promoter activity alone.

      We now also explicitly state that the data shown in Figure 2A derive from three independent electroporation experiments per construct, each assayed in technical duplicate.

      Lines 241–248: “Promoter–terminator activity was assessed 16–18 h after electroporation using a transient NanoLuc assay adapted from Li et al. (2019), with modifications. This early time point was selected to capture reporter output within the transient-expression window after DNA delivery, before prolonged plasmid loss, differential outgrowth, or culture-level changes could dominate the readout. Because the constructs did not contain a selectable marker, the measured NanoLuc signal reflects early transient reporter output rather than promoter strength independent of DNA uptake, early plasmid retention, or post-transfection recovery.”

      Additional clarification added to Figure 2 legend stating that measurements derive from three independent electroporation experiments, each assayed in technical duplicate.

      (3) Figure 2:

      (3.1) Several aspects of the current design may lead to ambiguity for the reader. The boxplots are colour-coded, but it is unclear whether the colours carry meaning or are purely decorative. Because the data are already spatially separated into bins, additional random colouring is redundant and may suggest distinctions that are not intended. In addition, part A of Figure 2 is split into two panels, with the scale for the left panel shown in the right panel and some of the boxplot colours falling in the range of the scale, but not in line with their counterparts in the left panel. Because the colour use is not consistent, it is difficult to tell whether the same scale should be applied to both panels or how it should be interpreted.

      (3.2) The left panel of part A uses a diverging blue-white-red colour scheme, which is most appropriate when the midpoint represents a meaningful central value such as zero. Because the values shown in this graph are only positive, a non-diverging 2-colour scale or a colour palette such as 'viridis' would make the plot easier to interpret.

      (3.3) A black background should be avoided: 'B' and 'C' labels are invisible, and it draws attention to a distracting design feature rather than the data themselves.

      We thank the reviewer for these detailed comments regarding figure design and visual interpretation. We agree that the original presentation of Figure 2 introduced unnecessary visual ambiguity through inconsistent colour usage, the use of a diverging colour scale for strictly positive values, and poor readability associated with the dark background and low-resolution export.

      In response, Figure 2 has been extensively redesigned to improve clarity, accessibility, and interpretability. The previous blue–white–red diverging heatmap has been replaced with a sequential colour palette appropriate for positive-only expression data. Boxplot colouring has also been simplified and harmonised with the heatmap scheme to avoid implying unsupported categorical distinctions. In addition, panel organisation, typography, scaling, and legend structure have all been revised to improve readability and reduce ambiguity regarding interpretation of the plotted values.

      We also agree that the black background distracted from the data presentation and impaired visibility of panel labels and image boundaries. The revised figures therefore use white backgrounds together with clearer panel separation and improved label visibility throughout.

      Figure 2 has been completely reformatted using a sequential colour scale in panel A, simplified and harmonised boxplot colouring, larger typography, improved panel separation, revised legends, and white backgrounds throughout. Corrected high-resolution source figures have been provided.

      (4) Figure 3:

      (4.1) Individual snapshots should be separated more clearly, either by using a white background or by adding visible borders to make the overall composition clearer. As currently displayed, some boundaries between fluorescent channels resemble image artifacts rather than intentional panel divisions.

      We thank the reviewer for this helpful comment regarding figure composition and panel separation. We agree that the original presentation made it difficult to distinguish intentional panel boundaries from imaging artefacts, particularly in the low-resolution review PDF generated during manuscript compilation.

      To improve clarity, Figure 3 has been reformatted using white backgrounds, clearer panel spacing, and more explicit separation between individual snapshots and imaging channels. High-resolution source images have also been provided to ensure that fluorescence patterns, image boundaries, and panel organisation remain clearly interpretable both on screen and in print.

      Figure 3 has been reformatted with improved panel separation, white backgrounds, clearer image boundaries, and revised high-resolution source figures.

      (4.2) In parts B-D, the legend should explain more clearly what each image shows, and the figure itself would benefit from annotations. There seem to be three sub-panels in each 'condition' of part B (as well as C and D): while the middle and rightmost panel can be easily inferred to represent the fluorescent protein and bright-field image, what the leftmost panels represent is not specified. If DAPI was used to dye DNA, an explanation why mostly multiple labelled regions are visible should be provided.

      We thank the reviewer for these helpful suggestions regarding figure annotation and legend clarity. We agree that the original presentation did not sufficiently explain the composition of the imaging panels, particularly under the low-resolution conditions of the review PDF.

      To improve interpretability, the revised Figure 3 now includes clearer panel organisation, improved annotations, and expanded figure legends explicitly identifying the individual imaging channels and staining conditions shown in each subpanel. The leftmost panels in parts B–D are now more clearly identified in both the figure and legend, together with the corresponding fluorescence or staining conditions used in each experiment.

      As mentioned in the Methods sections we used Hoechst 33342 to visualise DNA; but we agree that the Hoechst 33342-labelled structures required additional clarification. The revised legend section now explains that multiple Hoechst 33342-positive regions are commonly observed because Blastocystis cells can contain multiple nuclei depending on cell stage and subtype-specific morphology.

      In addition, high-resolution source images have been provided to ensure that fluorescent signals, panel boundaries, and imaging features remain clearly interpretable both on screen and in print.

      Figure 3 legends and annotations have been revised to clarify imaging channels, staining conditions, and panel organisation. The figure caption was also edited to include: “DNA was visualised using Hoechst 33342. Most cells contained two nuclei, and smaller Hoechst 33342-positive signals consistent with mitochondrial DNA were also observed in some instances.”

      (4.3) Cell morphology and appearance differ markedly between UnaG/smURFP and SNAP-tag images, which should be explained. A microscope issue is mentioned in the main text, but if that was the cause, the authors should consider replacing the images, as the current distortions complicate interpretation.

      We thank the reviewer for this important observation and agree that the apparent morphological differences between the UnaG/smURFP and SNAP-tag panels required additional clarification.

      The images shown for the different reporter systems were acquired under different imaging conditions and microscope configurations following an instrument-related issue during part of the imaging workflow, as noted in the Methods section. As a result, direct visual comparison of cell morphology between reporter systems is not appropriate. The primary purpose of these panels is instead to demonstrate reporter detectability, live-cell labelling capability, and the characteristic fluorescence patterns obtained with the different anaerobiosis-compatible reporter systems.

      In particular, the SNAP-tag panels were included to demonstrate successful live-cell labelling without permeabilisation together with the expected increase in fluorescence signal at higher substrate concentrations, rather than to support quantitative comparison of cell morphology across imaging conditions.

      We considered replacing the affected images. However, equivalent replacement datasets acquired under directly comparable conditions are not currently available. We have therefore retained the original images but revised the figure legend to clarify the intended interpretation and limitations of these panels explicitly.

      Figure 3 legend revised to include:

      “Because images for the different reporter systems were acquired under different imaging conditions, they are presented to demonstrate reporter detectability and labelling pattern and should not be used for quantitative comparison of cell morphology across reporter systems.”

      Reviewer #3 (Recommendations for the authors):

      The reader may find the current order confusing starting with construct design before testing which drug to use for selection. The narrative would work better if it started with antibiotic selection as the first logical step for generating stable cell lines.

      We thank the reviewer for this thoughtful suggestion regarding narrative structure and agree that multiple organisational strategies are possible for presenting a methodological workflow of this type.

      We considered reorganising the Results section to begin with antibiotic selection and drug sensitivity profiling. However, we ultimately retained the overall structure because the manuscript is organised as a toolkit-development framework rather than as a strictly chronological experimental protocol. The Results therefore begin with regulatory-element discovery and construct design, which form the conceptual and experimental foundation of the toolkit, before progressing to DNA delivery optimisation, drug sensitivity profiling, clonal recovery, and reporter validation.

      We felt that this structure most clearly reflects the dependency relationships within the system: regulatory elements are required before constructs can be assembled, constructs are required before electroporation conditions can be evaluated, and selectable constructs are required before stable selection and clonal recovery can be meaningfully assessed.

      (2) The text states that the screen 'focused on the 1,000 most abundant proteins to establish a preliminary library capable of supporting varying levels of transcription.' Since the genome has ~6,000 protein-coding genes, the top 1,000 cover the most abundant proteins — not a wide expression range.

      We thank the reviewer for this important clarification. We agree that the original wording could incorrectly imply that the screen was intended to sample broadly across the full transcriptional range of the Blastocystis genome. This was not the case, and we have revised the manuscript accordingly.

      Our strategy was instead designed to enrich for candidate loci with a higher prior likelihood of supporting detectable transgene expression. Because no genome-wide promoter map, transcription start site dataset, or experimentally validated regulatory annotation was available for Blastocystis ST7-B at the inception of this work, we used the abundance-ranked Blastocystis ST4-WR1 proteomic dataset of Armengaud et al. (2017) as a practical starting point for candidate discovery.

      Importantly, the Blastocystis ST4-WR1 proteome is highly skewed, with 193 proteins contributing approximately 50% of the detected proteome and the 13 most abundant proteins contributing approximately 10% (Armengaud et al., 2017). We therefore selected the top 1,000 proteins not as a representation of the genome-wide expression range, but as a proteomics-guided enrichment strategy to identify loci more likely to contain active endogenous regulatory regions suitable for initial toolkit development.

      We have revised the relevant Methods section substantially to clarify both the rationale and the workflow used for candidate selection, homolog identification, and promoter/terminator definition.

      The Methods section (Lines 156–189) has been extensively revised to clarify the rationale underlying candidate regulatory-element selection. The revised text now explicitly states that the strategy was designed to enrich for likely active loci for toolkit development rather than to systematically survey the full range of promoter strengths across the Blastocystis genome.

      Additional methodological detail has also been added regarding:

      Use of the Armengaud et al. (2017) proteomic and proteogenomic datasets,

      Homolog identification in Blastocystis ST7-B,

      Locus selection criteria,

      Promoter boundary definition,

      And operational definition of candidate terminator regions.

      (3) The Methods contain an inconsistency: cells were left in 0.5 mL, then 1 mL was added, but then only 0.5 mL is apparently used for transfection. What happened to the 1 mL?

      We thank the reviewer for identifying this ambiguity in the transfection workflow description. The apparent inconsistency arose because the protocol description moved from bulk cell resuspension to preparation of individual electroporation reactions without explicitly stating how the intermediate suspension was used.

      After washing, approximately 0.5 mL of cytomix buffer remained above the pellet, and 1 mL of complete cytomix buffer was then added to generate an approximately 1.5 mL cell suspension. Cells were counted from this pooled suspension, after which the volume corresponding to 5 × 10<sup>7</sup> cells was transferred into each individual electroporation reaction. Following addition of DNA, each electroporation reaction was adjusted to a final volume of 500 µL with complete cytomix buffer. The remaining cell suspension was retained for additional transfections or control reactions.

      We agree that the original wording could be misinterpreted and have revised the Methods section to clarify the sequential handling steps more explicitly.

      Lines 225-229 revised to read: “The resulting approximately 1.5 mL pooled cell suspension was used for total viable cell counting using a hemacytometer.”

      “After counting, the volume corresponding to 5 x 10<sup>7</sup> cells was transferred to each electroporation reaction and combined with 25 µg of plasmid DNA. The total electroporation volume was adjusted to 500 µL with complete cytomix buffer.”

      (4) Figures 2 and 3 are too low-resolution for the font size used and for clearly viewing the microscopy images.

      We thank the reviewer for highlighting these readability issues. As noted in our responses above regarding Figures 2 and 3, the low-resolution appearance primarily resulted from manuscript compilation and PDF export artefacts affecting typography, image rendering, and panel clarity in the review version.

      To address this, Figures 2 and 3 have been completely reformatted and replaced with revised high-resolution versions featuring improved typography, panel labelling, contrast, accessibility, and image clarity for both on-screen viewing and print reproduction.

      Revised high-resolution versions of Figures 2 and 3 have been provided as described above. No additional manuscript changes were required beyond the figure revisions already outlined.

      (5) Figure 4 is confusing because the left and right panels appear inconsistent, with much higher concentrations required for growth inhibition in the culture-based assay than the resazurin assay indicated. The rationale for the resazurin assay should be explained, and the complete growth inhibition (CGI) concentration should be highlighted in the right panel.

      We thank the reviewer for highlighting this potential source of confusion. We agree that the distinction between the two assay endpoints was not sufficiently emphasised in the original figure presentation and legend.

      The apparent discrepancy arises because the two assays measure different biological endpoints under different assay conditions. The resazurin assay was used to estimate IC<sub>50</sub> values, corresponding to the concentration at which metabolic activity was reduced by approximately 50% under the assay conditions. In contrast, the small-culture assay was designed to determine complete growth inhibition (CGI), defined operationally as the concentration at which no detectable culture outgrowth occurred after incubation, using phenol red acidification as a culture-level readout.

      Because these assays measure partial metabolic inhibition versus complete suppression of detectable culture outgrowth, the corresponding concentration ranges are not expected to coincide directly. The higher concentrations observed in the right-hand panels therefore reflect the more stringent endpoint associated with complete growth inhibition rather than inconsistency between the assays.

      We agree that this distinction should have been explained more clearly in the original manuscript. We have therefore substantially revised the Figure 4 legend to clarify the rationale underlying both assays, explicitly distinguish IC<sub>50</sub> and CGI endpoints, and explain how the CGI values were used to guide subsequent antibiotic selection conditions for Blastocystis ST7-B transformants. The CGI transition range has also been made more visually explicit in the revised figure presentation.

      Figure 4 caption revised to: “Antibiotic potency and selection-window determination in Blastocystis ST7-B. Dose–response curves for puromycin, trimethoprim, and WR99210 were estimated from a resazurin-based viability assay (n = 3 independent replicates per drug per concentration). Points show mean ± SD, and the insets list the estimated IC50 values with R<sup>2</sup>-values > 0.75 for all fitted curves. The IC<sub>50</sub> estimates represent the drug concentrations that reduced resazurin-based metabolic activity by 50% under the assay conditions.”

      Right panels: “small-culture complete growth inhibition assay using 1 × 10<sup>7</sup> WT Blastocystis ST7-B cells per culture, assayed in triplicate across a wide range of concentrations. Cultures were incubated for 2 days, and outgrowth was assessed using phenol red acidification of the medium as a culture-level readout, with yellow indicating growth and red indicating no detectable growth. The yellow-to-red transition was used to estimate the concentration required for complete growth inhibition and to guide the subsequent antibiotic selection strategy for Blastocystis ST7-B transformants.”

      “IC<sub>50</sub> and CGI represent distinct assay endpoints: the former measures partial reduction in metabolic activity, whereas the latter identifies the concentration at which no detectable culture outgrowth occurs under the small-culture assay conditions.”

      (6) In Figure 3B, the unexpected UnaG fluorescence pattern could be due to protein sequestration because the protein is mildly toxic to the cell. This should be discussed in addition to the reasons already provided.

      We thank the reviewer for this thoughtful suggestion and agree that protein sequestration or reporter-associated cellular stress represent plausible alternative interpretations of the observed UnaG fluorescence pattern.

      We considered the possibility of UnaG-associated toxicity during interpretation of these data. However, under the conditions tested, we did not observe clear evidence of a substantial toxic effect: UnaG-expressing Blastocystis ST7-B cells could be recovered as stable lines, maintained under antibiotic selection, and propagated through continued culture. We therefore felt that direct attribution of the observed fluorescence pattern to reporter toxicity would currently remain speculative.

      At present, we consider the biochemical properties of the UnaG system itself to provide a more parsimonious explanation for the observed localisation pattern. In particular, unconjugated bilirubin is highly hydrophobic and would be expected to partition preferentially into lipid-rich cellular environments. This interpretation is consistent with the lipid-rich peripheral and intracellular structures previously reported in Blastocystis ST7-B (Liao et al., 2023).

      We have therefore revised the Discussion to acknowledge that the observed UnaG fluorescence pattern may reflect a combination of reporter-specific biochemical behaviour, bilirubin partitioning, local intracellular environment, or possible sequestration phenomena. At the same time, we avoid assigning toxicity as a demonstrated mechanism in the absence of direct measurements of cell fitness, reporter abundance, or bilirubin distribution. Such experiments would be required to evaluate this possibility rigorously.

      Lines 642-648: “Consistent with this, lipid-rich peripheral and intracellular structures have been reported in Blastocystis ST7-B, potentially providing favourable microenvironments for BR partitioning and contributing to the punctate UnaG fluorescence pattern (Liao et al., 2023). An alternative possibility is that the observed signal pattern reflects reporter sequestration or reporter-associated cellular stress. However, because UnaG-expressing lines were recovered, maintained under selection, and propagated through continued culture, toxicity remains a possible but untested explanation rather than a demonstrated mechanism.”

      Minor Comments

      Figure 2: Parts B and C should also show individual datapoints for better reader assessment.

      We agree that inclusion of individual data points improves transparency and interpretability of the underlying data distributions.

      Individual data points have now been overlaid on the boxplots in Figures 2B and 2C.

      Figure 3A: Separate channels (fluorescence, bright-field, merge) should be shown rather than only the merge. The current overlay is difficult to interpret, especially for colour-blind readers.

      We appreciate the reviewer’s concern regarding accessibility and interpretability. We considered separating the fluorescence, bright-field, and merged channels for Figure 3A. However, this panel was intended primarily as an overview demonstrating reporter detectability within the bicistronic construct context, while the detailed fluorescence distribution is explored more extensively in the subsequent UnaG panels. We therefore retained the merged presentation for Figure 3A. Importantly, the image is not dependent on red–green discrimination, as it combines a greyscale bright-field background with a high-contrast green/cyan fluorescence signal that remains distinguishable through brightness and contrast differences. In addition, colour-blind-friendly lookup tables (LUTs) were used throughout the revised figure set.

      To further improve accessibility, the original red annotation arrow has been replaced with a colour-blind-friendly annotation colour.

      Briefly define system components (P2A, UnaG, smURFP, SNAP-tag) and add an abbreviation list.

      We agree that brief contextual definitions improve accessibility for readers less familiar with these reporter systems. Rather than adding a separate abbreviation list, we have added short explanatory descriptions at the points where these components are first introduced in the manuscript.

      Lines 394–396: “The P2A peptide is expected to promote ribosomal skipping during translation, allowing two separate polypeptides to be produced from a single open reading frame.” Line 515–516: “UnaG, a bilirubin-binding fluorescent protein originally isolated from the muscle of the Japanese eel (Kumagai et al., 2013)…” Line 527: “smURFP (small ultra-red fluorescent protein)…”

      Abstract: “among the most prevalent microbial eukaryote” should be “eukaryotes”.

      Corrected in revised manuscript.

      Conclusion (2nd sentence): unclear what “endogenous regulatory part discovery” means.

      We agree that this phrase required clarification. The intended meaning was the identification and benchmarking of native Blastocystis ST7-B promoter and terminator elements for construct design and toolkit development. We have clarified this directly in the revised Conclusion section.

      Lines 682–683 revised to: “By bringing endogenous regulatory part discovery, namely the identification of native promoter and terminator elements, …”

      Author contributions: “critical advise” should be “advice”.

      Corrected in revised manuscript.

      Again, we thank the reviewers for their careful evaluation, constructive criticism, and thoughtful feedback on the manuscript. The review process has substantially strengthened the manuscript by helping us clarify the distinction between what is directly demonstrated experimentally and what remains mechanistically unresolved.

      The central methodological conclusions of the study remain unchanged: the toolkit enables selectable transgene expression, recovery of colony-derived lines, and propagation of reporter-positive transgenic Blastocystis ST7-B lines, extending genetic accessibility in this organism substantially beyond the previous transient transfection framework.

      At the same time, the revised manuscript now more explicitly acknowledges important unresolved mechanistic questions, including vector topology, P2A-mediated protein separation efficiency, and persistence in the absence of selection. These are now discussed transparently together with the future experimental approaches that will be required to address them directly.

      We believe the revised manuscript now presents a clearer, more rigorous, and more accessible description of a practical genetic toolkit for Blastocystis ST7-B and hope that the revisions and clarifications satisfactorily address the reviewers’ concerns.

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

      Learn more at Review Commons


      Reply to the reviewers

      1. __ General Statements__ We thank the reviewers for their thoughtful and constructive evaluations of our work. We are particularly encouraged that both recognize the value of this study as a scalable and systematic framework for the functional exploration of the human KZFP family and agree that the resource generated here will be of broad interest to the KZFP, transposable element, and genome regulation communities. Reviewer 1 explicitly notes that "the screening framework itself represents a potentially useful resource for prioritizing candidate KZFPs for downstream study" and that "the study may nonetheless serve as a useful starting point for future investigations into KZFP biology and transcriptional regulation." Reviewer 2 similarly emphasizes that "the authors provide an efficient and valuable screening platform that can identify promising candidates for further investigation" and that "the methodological advance represents the primary contribution of the work."

      We can only concur with these assessments. The principal goal of this study was not to elucidate the physiological roles of all or even a subset of individual KZFPs, but rather to provide a scalable framework that enables their systematic prioritization and generates experimentally testable hypotheses regarding their functions. To support our argument, we ventured into some mechanistic analyses, but these could not pretend to be complete and definitive. In that respect, we hear the reviewers when they note that the original manuscript does not always sufficiently distinguish candidate discovery from mechanistic validation. In its revised version, we will therefore more clearly frame the inducible K562 overexpression assay as a standardized and sensitive readout of regulatory potency rather than as a direct surrogate of physiological function. Within this framework, K562 fitness defects are interpreted as a quantitative measure of the extent to which ectopic KZFP expression perturbs transcriptional homeostasis in a controlled cellular context, while the direct targets and transcriptional networks identified through our integrative analyses are presented as hypotheses to be tested in more physiologically relevant systems. Accordingly, the revised manuscript preserves the broad scope and resource aspect of the study while incorporating additional experimental validation, expanded methodological descriptions, and a more cautious interpretation of the proposed biological functions of the selected KZFPs.

      __Although this document is submitted as a Revision Plan, we have already incorporated a substantial number of revisions into the transferred manuscript. In particular, we have implemented most of the presentation, methodological, and conceptual modifications requested by the reviewers, including clarification of the scope of the study, extensive revisions of the Results and Discussion, expanded Materials and Methods, and numerous figure and text corrections. These revisions are detailed in Section 3 ("Description of the revisions that have already been incorporated into the transferred manuscript"). __

      The remaining points requiring additional experimentation or more extensive analyses are described in Section 2 ("Description of the planned revisions").

      __ Description of the planned revisions__

      Reviewer 1 Major comment 1

      “Finally, several aspects of the data presentation are currently difficult to reconcile. In Fig. 1D, the meaning of the purple category is unclear, and the percentage scaling on the x-axis is difficult to reconcile with the cumulative values displayed. For instance, the sum of all the bars would not reach 100%, as the values of the bars span percentages up to 4% at most (for 105 MYO KZFPs) according to this plot. Similarly, the reported numbers of TE-binding KZFPs in Fig. 1E-F and Fig. S1D appear internally inconsistent and should be clarified. Specifically, 53+14=67 KZFPs are reported to bind TEs in total, yet a larger number of KZFPs appears associated with individual TE families (e.g., 86 for LTR.ERV1). If the values shown correspond to percentages rather than absolute counts, this should be explicitly clarified in both the figure and legend. In addition, Fig. S1D appears inconsistent with the counts reported in Fig. 1E-F, as only 5 out of the 53 toxic KZFPs displayed in the plot show no enrichment for any of the highlighted TE families.”

      We thank the reviewer for this insightful comment, which has helped us identify areas where the presentation of our data can be substantially improved. We agree that the current presentation of the TE-binding analyses could be clearer and that revising these figures will improve both their readability and the overall consistency of the manuscript. In the revised manuscript, we will clarify the apparent inconsistencies in the presentation of the TE-binding KZFP analyses and revise the corresponding figures and legends accordingly. Importantly, these inconsistencies do not arise from errors in the underlying data but rather from an insufficient explanation of the statistical enrichment analyses and the way the results are represented. We will therefore redesign the relevant figures and expand their legends to more clearly describe the analytical approach, the enrichment criteria, and the interpretation of the results. We believe that these revisions will improve the clarity, transparency, and internal consistency of the manuscript, allowing readers to more readily interpret the TE-binding analyses. Minor comments of the reviewer 1 were extremely useful to detect mistakes and we are grateful for that. All the modifications that were asked see below were included in the manuscript.

      Reviewer 1 Major comment 2

      “Finally, while the proteomics results aimed at identifying SCAN-dependent interactors are of interest, several aspects of the experimental design and data analysis remain unclear. In particular, it is not specified whether the experiment was performed in biological replicates or as a single measurement. This is important, as it directly affects how the data can be interpreted and how stringent downstream filtering can be. In the Results section, the authors state that "we identified a set of SCAN-dependent interactors, i.e., proteins that co-immunoprecipitated with the full-length construct but were absent in controls and lost upon deletion of the SCAN domain," which suggests a relatively binary, "presence/absence" filtering strategy. However, this description does not specify whether any quantitative threshold (e.g., enrichment ratio) was applied when comparing full-length constructs to deletion mutants. In contrast, the Methods section states that "proteins lacking signal above background were excluded and proteins were additionally required to show stronger signal in at least one bait condition than in GFP controls based on heatmap clustering (see script)," which instead suggests that a threshold-based criterion was used to define enrichment relative to controls and deletion mutants. If this is the case, the exact criteria and thresholds used for filtering should be clearly stated and consistently reported between the Results and Methods sections. If replicate measurements were not performed, this should be explicitly acknowledged, as peptide-level variability may substantially influence the identification of high-confidence interactors, particularly if the applied cutoffs are not highly stringent.”

      We agree that a more detailed description of the experimental design and analysis strategy, together with additional validation, will strengthen the interpretation of the proteomic data. In the revised manuscript, we expanded the Results and Materials and Methods sections to provide a clearer and more quantitative description of the filtering strategy, including the enrichment criteria and thresholds used to define SCAN-dependent interactors. To further strengthen these findings, we propose to perform an independent biological replicate of the co-immunoprecipitation mass spectrometry experiment. This additional experiment will increase confidence in the identified SCAN-dependent interactors and further support the conclusions drawn from the proteomic analysis.

      Reviewer 1 Minor comments

      • In Fig. 2A, readability could be improved by adjusting the layering of points, as the darker dots (in particular the red ones) are currently obscured by lighter ones. Alternatively, removing the outline of the points (which is not transparent) may also improve visibility, but in that case the legend for point size would need to be updated accordingly.

      Thank you for this helpful suggestion. We will revise Figure 2A to improve its readability by reworking the layering of the points in accordance with the reviewer's recommendation. We will also evaluate the point outlines and, if appropriate, remove them and update the point-size legend accordingly to ensure the figure is clear and easy to interpret.

      Reviewer 2 – Major comment

      “- The authors looked at available chromatin data in either K562 cells or HEK293 cells, which I think is a very good way of utilizing publicly available data. Since the authors showed that different KZFPs might be functionally relevant in different cell types/tissues, I was wondering if they checked if there is available ChIP Seq or CUT&RUN data in those specific cell types/tissues. If yes, that data should be included in the manuscript.”

      We agree that integrating KZFP binding data generated in biologically relevant cell types or tissues would further strengthen the proposed regulatory models. As described in the revised manuscript, we have already adopted this approach for ZNF43 by integrating chromatin landscape data from thymus and liver, where suitable datasets were available.

      To further address this point, we propose to systematically explore publicly available ChIP-seq, CUT&RUN, CUT&Tag, and related chromatin profiling datasets for the other KZFPs investigated in this study. Where suitable datasets are available, these analyses will be incorporated into the revised manuscript to further support the proposed tissue-specific regulatory models and provide additional biological context for the identified target genes.

      __ Description of the revisions that have already been incorporated in the transferred manuscript__

      Reviewer 1 Major comment 1

      “The large-scale overexpression screen represents the foundation of the manuscript and provides a potentially valuable resource for prioritizing candidate KZFPs for downstream study. However, several aspects of the experimental setup and data presentation currently limit the interpretation of the reported proliferation defects. First, key details regarding the screening workflow remain unclear. While the Methods section describes the overall procedure, it is difficult to determine when cells were seeded relative to doxycycline induction, in which plate format the cells were maintained throughout the experiment, and whether medium exchange was performed during the 9-day assay. These points are particularly relevant given the use of suspension K562 cells (which can complicate medium exchange in a 96-well plate format and make long-term culture more difficult to control) and a metabolic viability readout (PrestoBlue), as differences in nutrient depletion or overgrowth could also influence the signal independently of reduced proliferation or toxicity. Additional clarification regarding seeding density, timing of induction, plate format, culture handling throughout the assay, and whether cell morphology/density was visually monitored would substantially improve interpretability and reproducibility. Second, it is unclear whether the observed proliferation phenotypes may be influenced by differences in transgene expression levels or integration effects. Were all constructs validated for comparable expression following induction? In the absence of such controls, it remains difficult to determine whether the reported phenotypes reflect specific KZFP activities or differences in overexpression efficiency. While it may not be possible to conclusively distinguish KZFP-specific effects from toxicity associated with high transgene expression levels, this limitation should at least be acknowledged. In addition, the possibility that some phenotypes may be influenced by transgene integration effects should also be considered. Unless independent transductions were validated for the KZFPs classified as toxic, it remains difficult to exclude integration-site-specific contributions to the observed proliferation defects. Third, the normalization strategy would benefit from additional clarification. In Fig. S1A, the LacZ control appears variably affected by doxycycline treatment across plates, whereas the GFP control appears more stable. Since normalization relies on the mean behavior of both controls within each batch and condition, the authors should clarify whether this variability could influence hit calling.”

      We agree that additional methodological details improve the clarity and reproducibility of the screening assay. Accordingly, we substantially expanded the Materials and Methods section to describe the experimental workflow, quality controls, data normalization, and hit-calling criteria. The revised paragraph is reproduced below.

      Arrayed overexpression screen

      To systematically assess the effect of human KZFP overexpression on cellular fitness, K562 cells were individually transduced with doxycycline-inducible lentiviral vectors encoding 366 human KZFPs. Lentiviral particles were produced as described above and used to transduce cells without MOI calculation. __Instead, a fixed volume of viral supernatant (200µL per × 104 cells in 48 well plate filled with 200ul of RPMI) was used for all transductions to ensure comparable experimental conditions. Transduced cells were selected with puromycin before doxycycline induction. Following puromycin selection (1µg/mL for 3 days), cells were seeded at 20 000 cells per well in 24-well plates filled with 1ml of medium in technical triplicate for each KZFP. Following puromycin selection and prior to doxycycline induction, cell survival was visually assessed as a quality control metric for each KZFP construct (Supp __Table 2____). Doxycycline (1µg/mL) was added immediately after cell seeding to induce expression of the HA-tagged KZFPs. At each time point, metabolic activity was measured using PrestoBlue™ reagent according to the manufacturer's instructions (10µL reagent added to 100µL culture medium, incubated for 3h in a 96 plates). Absorbance was recorded at 570 nm and 600 nm using a plate reader (Hidex Sense Microplate Reader), GFP- and LacZ-expressing control wells were included on every plate to account for plate-to-plate and batch-to-batch variability. Peripheral wells were filled with culture medium to minimize evaporation-induced edge effects. Cells were maintained in RPMI supplemented with 10% fetal bovine serum (FBS) and 1× penicillin–streptomycin, and splited (1/10) with aspiration of the surface medium every three days throughout the assay while maintaining doxycycline at 1µg/mL. Cell proliferation was assessed after 4, 7, and 9 days of induction. and the A570/A600 ratio was used as a surrogate measure of viable cell number and proliferative capacity. For computational normalization, raw A570/A600 values were first background-corrected by subtracting the signal from medium-only controls and then normalized in two steps. First, each value was divided by the mean signal obtained from the GFP and LacZ control wells from the corresponding batch and induction condition to correct for inter-batch variability. Second, the resulting value was normalized to the corresponding −Dox condition for the same KZFP and time point to correct for seeding variability, yielding a relative proliferation score that reflects the effect of KZFP induction. KZFPs with a normalized proliferation score ≤ 0.85 at day 9 were arbitrarily classified as proliferation-impairing hits in this screening framework.

      After doxycycline induction, dot blot analysis using anti-HA and anti-actin antibodies was systematically performed to assess KZFP expression and sample loading, respectively (Supplementary DotBlot.pdf). The HA signal following doxycycline induction (HA_Dox) and actin signal following doxycycline induction (Actin_Dox) were visually scored from the dot blot signals (__Supp __Table 2).

      In addition, to strengthen the methodological description and address these concerns more directly, we will:

      1/ Include a supplementary table summarizing our experimental observations for each individual KZFP throughout the screening process (See preliminary Supp Table 2). -> See header here:

      2/ Perform and include Dot Blot analyses, to assess and compare transgene expression levels across KZFP constructs. (Supplementary File DotBlot.pdf____). Generation of these files is in progress, with a few missing dot blots still being completed (we have done 303 over 366 already). However, preliminary versions have already been submitted. -> See header of the .pdf here:

      In addition, we agree that a more explicit discussion of the limitations of our screening approach improves the interpretation of our findings. Accordingly, we expanded the Discussion to address the limitations associated with variable transgene integration, heterogeneous transgene expression, potential toxicity due to ectopic KZFP overexpression, and the use of K562 cells as a standardized rather than physiological cellular model.

      “Several methodological considerations should be taken into account when interpreting these results. As with any lentiviral overexpression screen, three potential sources of technical variability may influence the observed phenotypes: differences in transgene integration sites, heterogeneity in transgene expression levels, and non-specific toxicity resulting from ectopic overexpression. Variable integration sites are unlikely to represent a major source of bias in the present study because all analyses were performed on polyclonal populations of transduced cells rather than individual clones, thereby averaging integration-site effects across many independent events. In contrast, heterogeneity in transgene expression levels is expected, as the abundance of each KZFP depends not only on transduction efficiency but also on intrinsic differences in mRNA stability, translational efficiency, and protein stability. To minimize these sources of variability, all constructs underwent systematic quality control, including assessment of cell survival following puromycin selection and evaluation of transgene expression by HA dot blot after doxycycline induction. Although transgene expression levels varied across KZFPs (Supplementary File DotBlot.pdf), this variability showed no systematic relationship with the proliferation phenotypes, suggesting that differences in overexpression efficiency are unlikely to be the primary determinant of toxicity. Nevertheless, ectopic expression exposes cells to supraphysiological concentrations of KZFPs capable of generating non-physiological interactions or regulatory effects. Therefore, while the screening strategy is well suited for identifying candidate functional regulators, independent validation under endogenous expression conditions remains essential to confirm KZFP-specific functions.”

      Reviewer 1 Major comment 2:

      “A central conceptual issue throughout the manuscript is that the downstream functional analyses of the selected KZFPs remain largely disconnected from the original screening phenotype. The four candidates were prioritized based on proliferation defects observed upon overexpression in K562 cells; however, the subsequent analyses (with the only exception being a more in-depth experimental analysis of ZNF498 in ciliogenesis, which stands out as comparatively more directly supported by experimental evidence) primarily rely on correlative expression patterns and KZFP ChIP-seq datasets to infer potential biological functions in unrelated cellular contexts. As a result, it remains unclear whether the proposed transcriptional programs are mechanistically linked to the proliferation phenotypes that motivated candidate selection in the first place. This issue is evident across multiple sections of the manuscript. For example, the proposed role of ZNF43 in regulating fatty acid metabolism and detoxification pathways is primarily inferred from tissue-level expression correlations. While these analyses focus on genes identified as potential ZNF43 targets, the underlying ChIP-seq datasets were themselves generated under ZNF43 overexpression conditions. Therefore, the current analyses do not establish whether ZNF43 regulates these pathways under physiological expression levels or within a relevant cellular context, nor how such regulation relates to the proliferation defect observed in K562 cells. Moreover, several proposed target genes remain substantially expressed in tissues where ZNF43 expression is not particularly low (e.g., kidney and heart muscle), suggesting that additional regulators are likely involved. Similarly, the proposed model of ZNF257-mediated regulation of MAGEA genes during spermatogenesis is intriguing but does not fully account for the expression behavior of all MAGEA family members, particularly MAGEA2B, which displays strong expression in spermatocytes despite high ZNF257 expression. This expression pattern should be acknowledged in the main text and reflected in Fig. 3K. In addition, the labels for MAGEA6 and MAGEA2B in Fig. 3C appear to be inverted. More broadly, the proposed regulatory model is difficult to reconcile with the generally restricted expression pattern of MAGEA genes across adult tissues, as their expression does not appear to consistently correlate with ZNF257 levels outside the germline context. Related concerns also apply to the analyses of ZNF498 and ZNF18, where the proposed functions in cilium formation and sperm maturation remain disconnected from the proliferation defects identified in the initial screen.”

      We agree that this comment raises an important conceptual point and has helped us clarify the scope of the study and the interpretation of our findings. In the revised manuscript, we explicitly distinguish hypothesis generation from mechanistic validation by clarifying that the proliferation phenotype observed in K562 cells reflects the regulatory potential of ectopically expressed KZFPs rather than their physiological functions. We also adopted a more cautious interpretation of the functional analyses, emphasizing that the proposed regulatory networks are hypothesis-generating and that individual KZFPs are unlikely to act as sole regulators. More broadly, we emphasize that the primary objective of this study is to establish a scalable screening platform for prioritizing KZFPs and identifying biologically relevant contexts for future investigation, rather than to provide a comprehensive functional characterization of individual KZFPs. We agree that this comment highlights an important limitation of our proposed regulatory model. In the revised manuscript, we adopted a more nuanced interpretation by presenting ZNF257 as a contributor to, rather than the sole regulator of, the MAGEA transcriptional program, and by explicitly discussing the exceptions identified by the reviewer.

      Modification in the revised manuscript:

      1/

      “Integrative transcriptomic, chromatin and proteomic analyses reveal diverse mechanisms, including transposable element–linked repression (ZNF43), promoter-proximal regulation (ZNF257), and SCAN domain–dependent transcriptional activation (ZNF498/ZSCAN25 and ZNF18).”

      Is now:

      “Integrative transcriptomic, chromatin and proteomic analyses identify distinct regulatory properties and generate testable hypotheses regarding diverse mechanisms, including transposable element-associated repression (ZNF43), promoter-proximal regulation (ZNF257), and SCAN domain-dependent transcriptional activation (ZNF498/ZSCAN25 and ZNF18).”

      2/

      “Detailed follow-up of four such candidates, ZNF43, ZNF257, ZNF498 and ZNF18, revealed as hypothesized distinct modes of action, ranging from TE-linked transcriptional repression to promoter-proximal gene silencing and SCAN domain-mediated transcriptional activation. These findings reinforce the view that KZFPs, while often viewed as a homogeneous family of TE-repressive TFs, are rather functionally diverse regulators with wide-ranging impacts on human biology.”

      Is now:

      “Detailed follow-up of four such candidates, ZNF43, ZNF257, ZNF498 and ZNF18, identified distinct regulatory properties and generated hypotheses regarding their physiological functions. By integrating overexpression-induced transcriptional responses, chromatin occupancy, proteomic analyses and tissue-specific expression data, we propose candidate biological contexts in which these KZFPs may operate. These hypotheses now provide a framework for future mechanistic studies performed under physiological conditions. Together, these findings reinforce the view that KZFPs, while often viewed as a homogeneous family of TE-repressive transcription factors, comprise functionally diverse regulators with broad potential roles in human biology.”

      3/

      “We conclude from these data that ZNF43 regulates a transcriptional program related to fatty acid metabolism and detoxification, allowing for the preferential expression of its effectors in the liver (Fig. 2G). Interestingly, neither expression nor chromatin state followed the same pattern at the functionally unrelated DNAI4 locus, indicating that this gene is subjected to other dominant regulators.”

      Is now:

      “Together, these observations identify a small set of candidates ZNF43 target genes involved in fatty acid metabolism and detoxification and suggest that ZNF43 may contribute to the regulation of these transcriptional programmes in appropriate physiological contexts (Fig. 2G). However, these conclusions are derived from overexpression-based datasets and tissue-level expression analyses and should therefore be considered hypothesis-generating. Interestingly, neither expression nor chromatin state followed the same pattern at the functionally unrelated DNAI4 locus, indicating that additional regulatory mechanisms contribute to the control of these genes.”

      4/

      “It strongly suggests that ZNF257 contributes to initiating the transcriptional repression of these two MAGEA genes during early spermiogenesis, after which their silencing may be stabilized through stable epigenetic mechanisms such as DNA methylation.”

      Is now:

      “These observations suggest that ZNF257 may contribute to the initiation of transcriptional repression of a subset of MAGEA genes during the spermatogonia-to-spermatocyte transition, after which their silencing may be stabilized through epigenetic mechanisms such as DNA methylation.”

      5/

      “Together, these results identify ZNF498 as a transcriptional activator of gene modules controlling cytoskeleton-dependent processes and suggest that this TF may act as a regulator of neuronal cytoskeletal architecture, warranting investigation in relevant neural models.”

      Is now:

      “Together, these results indicate that ZNF498 functions as a transcriptional activator in our overexpression system and support the hypothesis that it contributes to transcriptional programmes controlling cytoskeleton-dependent processes in physiologically relevant neural contexts, warranting further investigation in dedicated neural models.”

      6/

      “The co-expression of ZNF18 and its target genes at the spermatid stage suggests that ZNF18 activates a transcriptional program supporting these processes.”

      Is now:

      “The co-expression of ZNF18 and its candidate target genes at the spermatid stage is consistent with the hypothesis that ZNF18 contributes to transcriptional programmes supporting these processes.”

      7/

      “The four KZFPs characterised here illustrate this diversity. ZNF43 represses a coherent set of genes involved in fatty acid metabolism and detoxification through binding to nearby LTR/ERV1 integrants, with its expression anticorrelating that of its targets: i.e., highly expressed in thymus and bone marrow, where these metabolic genes are silent, and lowly expressed in liver, where they are most active. This represents a clear example of host genomes coopting TE-derived sequences and shaping their regulatory activities in a cell-type specific manner by the differential expression of KZFPs. ZNF257, by contrast, acts as a promoter-proximal repressor whose targets show accelerated sequence evolution at their promoters, consistent with integration into a KZFP-orchestrated GRN through rapid promoter diversification, a feature previously described for KZFPs (Farmiloe et al., 2023). Its regulation of the MAGEA gene cluster exemplifies a distinct evolutionary mechanism: an ancestral intronic binding site, present in MAGEA6 gene body, before ZNF257 emerged, was propagated across the cluster through tandem duplication, enabling coordinated regulation of multiple paralogs. Temporal expression analysis during spermatogenesis further suggests that ZNF257 initiates MAGEA repression at the spermatogonia-to-spermatocyte transition, after which silencing may be maintained through epigenetic mechanisms such as DNA methylation. ZNF498 and ZNF18, both SCAN-containing KZFPs with variant KRAB domains, on the other hand acted as transcriptional activators. ZNF498 activates a programme centred on microtubule cytoskeleton organisation, as demonstrated by the disruption of ciliogenesis upon its overexpression, and both ZNF498 and its targets are broadly expressed in the central nervous system, particularly in excitatory neurons where microtubule dynamics are essential for axonal architecture. ZNF18 similarly activates genes involved in chromatin remodelling and cytoskeletal reorganisation at the spermatid stage, processes that are hallmarks of spermiogenesis. Together, these case studies demonstrate that even within a single screen, KZFPs with fundamentally different regulatory logics can be identified through a single unifying phenotype and then mechanistically dissected to uncover their unique properties.”

      Is now:

      “The four KZFPs characterized here illustrate the functional diversity that can be uncovered using this screening strategy. For ZNF43, integration of overexpression transcriptomics with ChIP-exo binding data identified a small set of candidate direct target genes located near LTR/ERV1 elements. Their tissue-specific expression patterns are consistent with the hypothesis that ZNF43 contributes to transcriptional programmes associated with fatty acid metabolism and detoxification, although these analyses, which rely on overexpression-derived datasets and tissue-wide correlations, do not establish physiological regulation or causality. Rather, they identify a candidate regulatory network whose functional relevance will require investigation in appropriate biological models. More generally, these observations support the concept that host genomes may exploit TE-derived regulatory sequences in a tissue-specific manner through differential KZFP expression, while recognizing that additional transcription factors almost certainly participate in controlling these gene expression programmes. Similarly, ZNF257 emerged as a promoter-associated transcriptional repressor in our overexpression system. Evolutionary analyses suggest that tandem duplication propagated an ancestral ZNF257-binding sequence across the MAGEA locus, generating the hypothesis that ZNF257 may contribute to coordinated regulation of this gene cluster during spermatogenesis. The temporal expression profiles of ZNF257 and the MAGEA genes are compatible with such a model but remain correlative and therefore require direct functional validation. ZNF498 and ZNF18, two SCAN-containing KZFPs with variant KRAB domains, displayed transcriptional activation rather than repression following overexpression. For ZNF498, the integration of transcriptomic analyses with expression profiling pointed to microtubule cytoskeleton organization as a candidate biological process, a prediction that was further supported experimentally by the marked impairment of ciliogenesis following ZNF498 overexpression in hTERT-RPE1 cells. This represents the strongest functional validation presented in this study and supports the biological relevance of the analytical framework developed here. For ZNF18, the co-expression of the KZFP and its candidate target genes during spermatogenesis is consistent with the hypothesis that it contributes to transcriptional programmes involved in chromatin remodelling and cytoskeletal reorganization during spermatid differentiation. Together, these case studies illustrate how a standardized overexpression screen can identify KZFPs with distinct regulatory properties and generate biologically coherent hypotheses regarding their physiological functions. Rather than establishing definitive functions for individual KZFPs, this framework prioritizes candidates, proposes relevant cellular contexts, and provides a foundation for future mechanistic studies performed under physiological conditions.”

      “In addition, interpretation of the SCAN-deletion experiments is complicated by the reduced expression levels of the deletion constructs relative to the corresponding full-length proteins, making it difficult to determine whether the observed proliferation phenotypes are pathway-specific or partially driven by differential expression.”

      We thank the reviewer for this important observation and agree that differences in expression levels between the full-length and ΔSCAN constructs could complicate the interpretation of the observed phenotypes. To address this concern, we performed a quantitative comparison of the expression levels of full-length and ΔSCAN proteins using both western blotting and transgene expression using RNAseq, while accounting for differences in transgene length. This result are now added in (Fig S6C, D).

      With modification of the legend:

      • HA signal after OE of HA-tagged ZNF18, ZNF18∆SCAN, ZNF498, ZNF498∆SCAN or GFP in K562 cells. Actin as control.
      • Quantification of ZNF18, ZNF18∆SCAN, ZNF498, ZNF498∆SCAN It appears that the difference is small (Minor comments of the reviewer 1

      “- In the Abstract and in the "Limitations of the study" section, the term "annotation" is used. It would be preferable to specify "functional characterization" instead of "annotation".

      Done as suggested by the reviewer.

      • In the Introduction, there may be a minor citation confusion. Following the sentence: "Characterized by an N-terminal KRAB domain and a C-terminal tandem array of C2H2 zinc fingers, KZFPs primarily target transposable element (TE)-embedded sequences," the cited references are predominantly experimental studies supporting this statement. However, the inclusion of the review "Bruno, Mahgoub and Macfarlan, 2019" appears less appropriate in this context, as it does not directly present ChIP-seq data supporting this claim. More relevant primary studies from the same research area include "Wolf et al. 2020" and "Bruno et al. 2025.".

      Done as suggested by the reviewer.

      • In Fig. 1A, "D10" appears inconsistent with the text and other figures (Fig. 1B, 1G, 1H), which refer to 9 days post-induction.

      Done as suggested by the reviewer.

      • In Fig. S1, there may be a mismatch in the highlighted plate: the zoomed image appears to correspond to the first plate from the top. The correct plate should be highlighted for consistency.

      Done as suggested by the reviewer.

      • In Fig. 1B, there is a typographical error ("K ZFPs" instead of "KZFPs").

      Done as suggested by the reviewer.

      • In Fig. S1E, it is unclear what "other" refers to. Please clarify whether this represents the mean of all remaining KZFPs or a defined subset, ideally in the figure description.

      Done as suggested by the reviewer.

      • In Fig. S2E, "SetDB1" should be corrected to "SETDB1".

      Done as suggested by the reviewer.

      • In Fig. 3B, it is unclear what distinguishes the upper and lower "Diverse REs". A brief clarification in the figure legend would improve interpretability, particularly regarding the transposable element families included.

      Done as suggested by the reviewer.

      • In Fig. S3C, the x-axis labels appear slightly misaligned and shifted to the right.

      Done as suggested by the reviewer.

      • In Fig. 3C, the labels for MAGEA6 and MAGEA2B appear to be inverted.

      Done as suggested by the reviewer.

      • In Fig. 3K, "MAGE3" should be corrected to "MAGEA3".

      Done as suggested by the reviewer.

      • In the ZNF498 section, line 4, the punctuation should be corrected so that the period appears after the figure reference ("promoters (Fig. S1E).").

      Done as suggested by the reviewer.

      • In the final sentence of the ZNF498 section, a noun appears to be missing after "cytoskeleton-dependent," possibly "processes".

      Done as suggested by the reviewer.

      • In the last section of the Results and corresponding figures and their descriptions, "SCAN dependant" should be corrected to "SCAN-dependent".”

      Done as suggested by the reviewer.

      Major comments of the reviewer 2

      “- The authors chose four KZFPs to study in detail, but why they chose these 4 candidates is unlcear to me. It would be nice to add a more detailed description of the process by which they chose the four candidates.”

      We agree that the rationale for selecting the four KZFPs should be presented more explicitly. Accordingly, we revised the manuscript to clarify the selection criteria.

      “However, a modest correlation was noted between the number of transcription start sites (TSS) bound by KZFPs and the drop in PrestoBlue signal induced by their overexpression (Fig. 1G), and SCAN-containing KZFPs (SKZFPs) tended to induce proliferation defects more frequently than family members lacking this domain (Fig. 1H).”

      Is now:

      “However, a modest correlation was noted between the number of transcription start sites (TSS) bound by KZFPs and the drop in PrestoBlue signal induced by their overexpression (Fig. 1G), and SCAN-containing KZFPs (SKZFPs) tended to induce proliferation defects more frequently than family members lacking this domain (Fig. 1H). These observations indicated that KZFPs affecting proliferation do not constitute a homogeneous functional group, prompting us to select representative candidates spanning the evolutionary, structural, and genomic diversity of the KZFP family for mechanistic characterization.____”

      “- The materials and methods part of the manuscript is not detailed enough for other researchers to reproduce the study. They should add more details to both experiments and data analysis part of this section. Below I highlight some examples for sake of clarity, but the authors should revise the whole materials and methods section and add more details keeping these examples in mind:

      • The authors do not state the titer of lentiviral vectors they generate nor the MOI or amount of virus they use to transduce the cells

      • In many cases, the specific softwares and the software version is not stated e.g., the analysis of the Gene Ontology Biological Processes

      • It would be beneficial for the readers to get more details about the construct they used, for example a map of the plasmid.

      • It is unclear how many cells were used for RNA extraction

      • It is unclear which microscopes were used for imaging.

      • The concentration of antibodies used for staining and the product number, and provider of the antibody is not always depicted.”

      We agree that the additional methodological details requested by the reviewer will improve the reproducibility and transparency of the study. Accordingly, we have expanded the Methods section to provide a more detailed description of the experimental procedures and data analysis workflow.

      “Lentiviral particles were produced in HEK293T cells by transient co-transfection of transfer, packaging and envelope plasmids. Cells were transfected at approximately 70–80% confluence using a standard lipid-based transfection reagent. Viral supernatants were collected 48 h after transfection, cleared by centrifugation, filtered through 0.22-µm membranes, and used fresh or stored appropriately until use. Recipient K562 or hTERT-RPE1 cells were transduced under conditions optimized for efficient gene delivery.”

      Is now:

      “Lentiviral particles were produced in HEK293T cells. 105 cells were seeded in 24 well plates filled with 1ml DMEM the day before transfection. Cells were co-transfected individually with 0.15ug of each plasmids encoding KZFPs tagged with HA (pTRE-KZFPX-HA-PGK-puro), 0.1ug of the packaging plasmid (pR8.74) and 0.07ug of the envelope plasmid (pMD2G) using TransIT®-LT1 Transfection Reagent (MIR 2306), according to the manufacturer's instructions. Viral supernatants were harvested 24h after transfection, clarified by centrifugation, filtered through 0.45-µm filters and used immediately.”

      “Coding sequences were cloned into doxycycline-inducible lentiviral transfer vectors designed to express N-terminally HA-tagged proteins.”

      Is now:

      “Coding sequences corresponding to 366 human KZFP open reading frames were codon-optimized for human expression and cloned into doxycycline-inducible lentiviral transfer vectors expressing C-terminal HA-tagged proteins under the control of a tetracycline-responsive promoter pTRE-KZFPX-HA-PGK-puro. All expression constructs used in the primary overexpression screen have been deposited and are publicly available (De Tribolet et al., 2023). A schematic representation of the lentiviral expression cassette, including the promoter, HA tag, cloning site, antibiotic resistance cassette, and regulatory elements, is provided in Supplementary file. Selected constructs encoding ZNF43, ZNF257, ZNF498 and ZNF18 were used for follow-up mechanistic studies. For SCAN-domain functional analyses, deletion constructs lacking the SCAN domain (ΔSCAN) were generated for ZNF18 and ZNF498 in the same lentiviral backbone. Deletion were done using In-Fusion cloning with specific primers. PCR was performed with high-fidelity polymerase, followed by gel purification and recombination with the linearized plasmid using the In-Fusion HD Cloning Kit (Takara Bio©) according to the manufacturer’s protocol. The product was transformed into HB101 Escherichia coli cells, and colonies were screened by PCR. Positive clones were verified by Sanger sequencing, and confirmed plasmids were propagated and purified for further use.”

      “Total RNA was extracted...”

      Is now:

      “For each biological replicate, approximately 1 × 10⁶ K562 cells were harvested 72 h after doxycycline induction. Total RNA was extracted…”

      “Images were acquired by fluorescence microscopy under identical conditions across samples.”

      Is now:

      “Images were acquired using a confocal microscope Leica-SP8 (Leica Biosystems) with an objective HC PL APO 63x/1.40 and a pinhole size of 1 AU, using identical acquisition settings for all conditions. Images were processed using Fiji/ImageJ (version 2.9.0) without nonlinear intensity adjustments.”

      “Cells were fixed and stained with antibodies against ciliary markers (ARL13B)”

      Is now:

      “Cells were fixed in 4% paraformaldehyde, permeabilized with 0.1% Triton X-100, blocked with 2% BSA, and incubated with rabbit anti-ARL13B (Proteintech, Cat. No. 17711-1-AP, 1:200) followed by Alexa Fluor 568-conjugated donkey anti-rabbit IgG (Thermo Fisher Scientific, Cat. No. A-10042, 1:1000). Nuclei were stained with Hoechst (1 µg/mL).”

      “- The authors mention that KZFPs are usually expressed at a low level in the K562 cell line they use, but there is no figure showing the expression level of KZFPs in this cell type. It would be important to see the baseline KZFP expression in these cells, the level of overexpression and compare it to the endogenous expression levels they show in different cell types/tissues, at least for the four candidates studied more in depth. This would help to understand whether this level of activity is something that could occur naturally in a physiologically relevant context.”

      We thank the reviewer for this insightful suggestion and fully agree that providing additional context regarding endogenous and ectopic KZFP expression levels will help readers better assess the physiological relevance of our findings. As suggested, we included data showing the baseline expression levels of the four selected KZFPs in K562 cells together with the expression levels achieved following doxycycline-induced overexpression. We also compared these values with publicly available transcriptomic data from cell lines. Importantly, only cell lines are assessed as we need ground through (K562) to estimate transgene expression. We modified Fig. S2, Fig. S3, Fig. S4 and Fig. S5 to add the results of these analysis. Here is ZNF43 as an example:

      With the following legend:

      “(C) Distribution of endogenous expression levels, (using GFP control cells), of all expressed genes (light grey) and all KZFPs (dark grey) in K562 cells. The solid red line indicates endogenous ZNF43 expression in GFP control cells, whereas the dashed red line indicates the corrected transgene expression following doxycycline induction.

      (D) Endogenous ZNF43 expression across Human Protein Atlas cell lines, (https://www.proteinatlas.org/about/download#cell_line), following normalization to the local RNA-seq dataset. K562 cells are highlighted in red. The dashed red line indicates the corrected transgene level measured following doxycycline-induced overexpression in K562 cells overexpressing ZNF43.”

      Modified the result section:

      “ZNF43 is a ~43-million-year-old KZFP with a canonical TRIM28-recruiting KRAB domain and 19 zinc fingers that preferentially recognize an LTR/ERV1-embedded sequence (Fig. S1F). We first verified that ZNF43 overexpression impaired the growth of K562 cells (Fig. S2A, B). Endogenous ZNF43 expression was readily detectable in K562 cells and across human cell lines (Fig. S2C, D). Following doxycycline induction, transcript abundance markedly increased and exceeded the highest endogenous expression level observed among the analyzed cell lines (Fig. S2C, D).”

      We also updated the Methods section:

      Quantification of endogenous and transgene expression levels

      Endogenous KZFP expression in K562 cells was estimated from GFP control RNA-seq samples using normalized mean expression values obtained from the differential expression analyses. For ZNF18, whose transgene sequence is identical to the endogenous coding sequence (i.e., not codon-optimized), transgene-derived expression was estimated directly by subtracting the endogenous transcript abundance measured in GFP controls from the total transcript abundance measured following doxycycline induction (OE − GFP). For ZNF43, ZNF257 and ZNF498, the overexpression constructs were synthesized using codon-optimized coding sequences. RNA-seq reads were therefore additionally aligned against the codon-optimized transgene reference sequences to specifically quantify exogenous transcripts without interference from endogenous reads. Because these codon-specific counts are generated through an independent alignment strategy, they are not directly comparable to the endogenous RNA-seq expression values. To calibrate these measurements, a scaling factor was derived from the ZNF18 dataset by comparing the codon-specific read counts with the transgene abundance estimated from the differential expression analysis (OE − GFP). This empirically determined correction factor was subsequently applied to all codon-optimized constructs, thereby expressing transgene abundance on the same scale as the endogenous RNA-seq measurements. Corrected transgene expression values were then used for all downstream comparisons. To compare endogenous expression across physiological contexts, publicly available RNA-seq datasets from the Human Protein Atlas (cell lines) were downloaded and normalized to the local RNA-seq scale. A normalization factor was calculated from the median expression ratio of KZFPs detected in both the Human Protein Atlas K562 dataset and the local K562 GFP control RNA-seq dataset, and subsequently applied uniformly to all Human Protein Atlas datasets. This normalization enabled direct comparison of endogenous expression across biological contexts with the corrected transgene expression values. Global KZFP expression was calculated as the median normalized expression of all annotated KZFPs within each biological context. For the four KZFPs selected for detailed characterization, endogenous expression across Human Protein Atlas cell lines was compared with corrected transgene expression following doxycycline induction. Expression distributions of all genes and KZFPs were visualized using ranked expression plots and density histograms. All analyses were performed in R using the tidyverse package.”

      We fully acknowledge that the overexpression system used in this study was primarily designed as a discovery platform to identify candidate functions, targets, and interaction partners of KZFPs that are otherwise expressed at lower levels in K562 cells. As the reviewer correctly points out, determining whether these regulatory effects occur at endogenous expression levels in physiologically relevant cellular contexts represents an important next step. We Thereby also clarified this in the “Limitations to this study” paragraph:

      “To better place our experimental system into a physiological context, we compared endogenous KZFP expression in K562 cells with publicly available transcriptomic datasets from the Human Protein Atlas. These analyses showed that K562 cells do not exhibit unusually low global KZFP expression compared with other human cell lines. However, consistent with the restricted expression patterns of this protein family, KZFPs as a whole are expressed at substantially lower levels than the average human gene. For the four KZFPs characterized in detail, doxycycline induction produced transcript levels that exceeded the highest endogenous expression observed across the analyzed human cell lines. Accordingly, the overexpression system used in this study was not designed to recapitulate physiological expression levels but rather to maximize the identification of candidate target genes, interacting partners, and regulatory pathways for KZFPs that are otherwise expressed at low endogenous levels. Consequently, the molecular interactions identified here should be considered as hypotheses requiring validation under endogenous expression conditions in physiologically relevant cellular models.”

      “- RNA seq analysis: It is unclear how many cells were used in the RNA seq analysis, I would like to ask the authors to clarify that. Moreover, from my understanding the RNA seq analysis was done on day 3, while the Presto Blue analysis was done on days 4, 7 and 9. I would like to kindly ask the authors to motivate their choice for the day of the RNA sequencing analysis.”

      We agree that this information required clarification. The Methods section has been revised to specify the number of cells used for RNA-seq library preparation and to explain the rationale for performing RNA-seq after 3 days of doxycycline induction, before measurable proliferation defects emerge, in order to capture primary transcriptional responses to KZFP overexpression. The corresponding modification has also been added to the Results section when introducing the RNA-seq analyses.

      “For transcriptome profiling, K562 cells expressing the indicated inducible constructs were treated with doxycycline for 72 h before harvest. Total RNA was extracted using the NucleoSpin RNA plus kit (Macherey-Nagel) according to the manufacturer’s recommendations. RNA quantity and purity were assessed by spectrophotometry, and RNA integrity was evaluated before library preparation.”

      Is now:

      “For transcriptome profiling, 1 × 10⁶ K562 cells expressing the indicated inducible constructs were treated with doxycycline for 72 h before harvest. RNA was collected after 3 days of induction to capture the primary transcriptional responses to KZFP overexpression before substantial differences in proliferation became apparent. This early time point was chosen to minimize secondary transcriptional changes resulting from altered cell growth, cell-cycle distribution, or cellular stress, which become detectable in the proliferation assays performed after 4, 7, and 9 days of induction. Total RNA was extracted using the NucleoSpin RNA plus kit (Macherey-Nagel) according to the manufacturer’s recommendations. RNA quantity and purity were assessed by spectrophotometry, and RNA integrity was evaluated before library preparation.”

      “We then profiled the transcriptome of K562 cells overexpressing ZNF43 by deep RNA sequencing (RNA-seq)”

      Is now:

      “We then profiled the transcriptome of K562 cells overexpressing ZNF43 by deep RNA sequencing (RNA-seq) after 3 days of doxycycline induction, a time point selected to capture primary transcriptional responses before the onset of measurable proliferation defects.”

      Minor comments of the reviewer 2

      “- Figure S1D is not mentioned in the text before figure S1E. The order of the panels should be changed in the figure.

      Done as suggested by the reviewer.

      • "We selected genes that were downregulated upon ZNF43 overexpression and harboured a ZNF43 binding site within 10kb of their TSS (Fig. 1A) - don't the authors mean Fig. 2A?

      Done as suggested by the reviewer.

      • In Figure 4D, the GO terms cannot be read, as the sentences seem to be cut.

      Done as suggested by the reviewer.

      • All figures and figure legends need to be revised. In some cases, the letter size is too small, or the legend and explanation of colours is missing. Please see some examples below: Fig. S6C, Fig 6C, Fig S4C, Fig S5C (letter size too small) Fig S6G, Fig 4E (label/scale is missing)”

      Homogenized to Arial 6 by default as requested by most of journal guidelines

      __ Description of analyses that authors prefer not to carry out__

      We think that by proceeding as described above we will have addressed all major conceptual issues raised by the reviewers.

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

      Learn more at Review Commons


      Referee #2

      Evidence, reproducibility and clarity

      Summary

      Foley et al establishes a scalable framework to probe KZFP function. They performed an array of inducible overexpression screen of 366 human KZFPs in K562 cells. This screen, together with the analysis of transcriptomic and available chromatin and proteomic datasets revealed that KZFPs regulate many different mechanisms, highlighting the functional diversity of KZFPs. Understanding this functional diversity is a very interesting, timely and relevant question, but it is also challenging to study. Therefore, the approach the authors develop is promising. While the quality of the experiments and data analysis is high, the weakness I see in the manuscript is the lack of major biological insights in relevant model systems. Please see my detailed comment in the significance part.

      Major comments

      • The authors chose four KZFPs to study in detail, but why they chose these 4 candidates is unlcear to me. It would be nice to add a more detailed description of the process by which they chose the four candidates.
      • The materials and methods part of the manuscript is not detailed enough for other researchers to reproduce the study. They should add more details to both experiments and data analysis part of this section. Below I highlight some examples for sake of clarity, but the authors should revise the whole materials and methods section and add more details keeping these examples in mind:
        • The authors do not state the titer of lentiviral vectors they generate nor the MOI or amount of virus they use to transduce the cells
        • In many cases, the specific softwares and the software version is not stated e.g. the analysis of the Gene Ontology Biological Processes
        • It would be beneficial for the readers to get more details about the construct they used, for example a map of the plasmid.
        • It is unclear how many cells were used for RNA extraction
        • It is unclear which microscopes were used for imaging.
        • The concentration of antibodies used for staining and the product number, and provider of the antibody is not always depicted.
      • The authors looked at available chromatin data in either K562 cells or HEK293 cells, which I think is a very good way of utilizing publicly available data. Since the authors showed that different KZFPs might be functionally relevant in different cell types/tissues, I was wondering if they checked if there is available ChIP Seq or CUT&RUN data in those specific cell types/tissues. If yes, that data should be included in the manuscript.
      • The authors mention that KZFPs are usually expressed at a low level in the K562 cell line they use, but there is no figure showing the expression level of KZFPs in this cell type. It would be important to see the baseline KZFP expression in these cells, the level of overexpression and compare it to the endogenous expression levels they show in different cell types/tissues, at least for the four candidates studied more in depth. This would help to understand whether this level of activity is something that could occur naturally in a physiologically relevant context.
      • RNA seq analysis: It is unclear how many cells were used in the RNA seq analysis, I would like to ask the authors to clarify that. Moreover, from my understanding the RNA seq analysis was done on day 3, while the Presto Blue analysis was done on days 4, 7 and 10. I would like to kindly ask the authors to motivate their choice for the day of the RNA sequencing analysis.

      Minor comments

      • Figure S1D is not mentioned in the text before figure S1E. The order of the panels should be changed in the figure.
      • "We selected genes that were downregulated upon ZNF43 overexpression and harboured a ZNF43 binding site within 10kb of their TSS (Fig. 1A) - don't the authors mean Fig. 2A?
      • In Figure 4D, the GO terms cannot be read, as the sentences seem to be cut.
      • All figures and figure legends need to be revised. In some cases, the letter size is too small, or the legend and explanation of colours is missing. Please see some examples below: Fig. S6C, Fig 6C, Fig S4C, Fig S5C (letter size too small) Fig S6G, Fig 4E (label/scale is missing)

      Significance

      Understanding the diverse roles of KZFPs is an important and interesting research question. However, studying KZFPs is challenging, as many KZFP-mediated effects appear to be highly cell type- and tissue-specific. This complexity is also highlighted by the findings of the current manuscript.

      A major strength of this study is the development of a scalable system that enables the simultaneous investigation of the entire KZFP family. Performing such analyses on an individual basis would be extremely time-consuming. Therefore, the authors provide an efficient and valuable screening platform that can identify promising candidates for further investigation. In this regard, the methodological advance represents the primary contribution of the work.

      At the same time, the study lacks a clear biological conclusion. While the screen identifies KZFPs with potential functional effects, it would substantially increase the impact of the manuscript if the authors selected at least one candidate for in-depth characterization in a biologically relevant cellular context. The current study is still of high quality and importance without these experiments, but such follow-up analyses would greatly strengthen the biological significance of the findings.

      Another limitation is that the experiments were performed in a cell type in which many of the investigated KZFPs are not normally expressed. As a result, the forced overexpression strategy may not accurately reflect physiological conditions and could potentially generate false-positive results. This concern is particularly relevant in light of the authors' statement that "KZFPs with sufficient regulatory potency to perturb cellular fitness outside of their normal setting are strong candidates for playing important roles within it." While this may indeed be true for some KZFPs, it is also possible that certain observed phenotypes simply arise from ectopic expression in an inappropriate cellular environment.

      More generally, the observation that KZFPs can have functions beyond TE repression is already established in the literature. Therefore, the manuscript provides limited new biological insight into this concept. The authors could potentially strengthen the novelty of the study by placing greater emphasis on specific KZFP subfamilies, such as SCAN-containing zinc finger proteins, which are a novel direction and have been implicated in non-canonical regulatory roles.

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

      Learn more at Review Commons


      Referee #1

      Evidence, reproducibility and clarity

      Summary

      In the manuscript "An overexpression platform reveals the functional diversity of human KRAB-Zinc Finger Proteins in maintaining cellular homeostasis", the authors describe a scalable framework aimed at prioritizing individual KZFPs for functional characterization. The study is centered on a large-scale screening approach in which human erythroleukemia K562 cells were transduced with inducible constructs enabling the overexpression of 366 individual KZFPs. The effects of KZFP overexpression on cellular proliferation were then assessed using a cell viability assay comparing induced and non-induced conditions.

      Based on the results of this primary screen, the authors selected four KZFPs for further investigation among those whose overexpression was associated with proliferation defects. Follow-up analyses integrated tissue- and cell type-specific expression patterns of these KZFPs with expression analyses of putative target genes identified from previously published ChIP-seq datasets, with the aim of inferring potential biological functions. On the basis of these analyses, the authors propose that the selected KZFPs may regulate distinct gene networks involved in processes including fatty acid metabolism, spermatogenesis, and other aspects of cellular homeostasis.

      The study spans multiple biological contexts, but the evidence supporting several of the individual conclusions remains relatively preliminary, and the breadth of the manuscript often comes at the expense of mechanistic depth.

      Major comments

      1. Interpretation and robustness of the initial overexpression screen

      The large-scale overexpression screen represents the foundation of the manuscript and provides a potentially valuable resource for prioritizing candidate KZFPs for downstream study. However, several aspects of the experimental setup and data presentation currently limit the interpretation of the reported proliferation defects. First, key details regarding the screening workflow remain unclear. While the Methods section describes the overall procedure, it is difficult to determine when cells were seeded relative to doxycycline induction, in which plate format the cells were maintained throughout the experiment, and whether medium exchange was performed during the 9-day assay. These points are particularly relevant given the use of suspension K562 cells (which can complicate medium exchange in a 96-well plate format and make long-term culture more difficult to control) and a metabolic viability readout (PrestoBlue), as differences in nutrient depletion or overgrowth could also influence the signal independently of reduced proliferation or toxicity. Additional clarification regarding seeding density, timing of induction, plate format, culture handling throughout the assay, and whether cell morphology/density was visually monitored would substantially improve interpretability and reproducibility. Second, it is unclear whether the observed proliferation phenotypes may be influenced by differences in transgene expression levels or integration effects. Were all constructs validated for comparable expression following induction? In the absence of such controls, it remains difficult to determine whether the reported phenotypes reflect specific KZFP activities or differences in overexpression efficiency. While it may not be possible to conclusively distinguish KZFP-specific effects from toxicity associated with high transgene expression levels, this limitation should at least be acknowledged. In addition, the possibility that some phenotypes may be influenced by transgene integration effects should also be considered. Unless independent transductions were validated for the KZFPs classified as toxic, it remains difficult to exclude integration-site-specific contributions to the observed proliferation defects. Third, the normalization strategy would benefit from additional clarification. In Fig. S1A, the LacZ control appears variably affected by doxycycline treatment across plates, whereas the GFP control appears more stable. Since normalization relies on the mean behavior of both controls within each batch and condition, the authors should clarify whether this variability could influence hit calling. Finally, several aspects of the data presentation are currently difficult to reconcile. In Fig. 1D, the meaning of the purple category is unclear, and the percentage scaling on the x-axis is difficult to reconcile with the cumulative values displayed. For instance, the sum of all the bars would not reach 100%, as the values of the bars span percentages up to 4% at most (for 105 MYO KZFPs) according to this plot. Similarly, the reported numbers of TE-binding KZFPs in Fig. 1E-F and Fig. S1D appear internally inconsistent and should be clarified. Specifically, 53+14=67 KZFPs are reported to bind TEs in total, yet a larger number of KZFPs appears associated with individual TE families (e.g. 86 for LTR.ERV1). If the values shown correspond to percentages rather than absolute counts, this should be explicitly clarified in both the figure and legend. In addition, Fig. S1D appears inconsistent with the counts reported in Fig. 1E-F, as only 5 out of the 53 toxic KZFPs displayed in the plot show no enrichment for any of the highlighted TE families. 2. Relationship between the screening phenotype and the proposed biological functions of selected KZFPs

      A central conceptual issue throughout the manuscript is that the downstream functional analyses of the selected KZFPs remain largely disconnected from the original screening phenotype. The four candidates were prioritized based on proliferation defects observed upon overexpression in K562 cells; however, the subsequent analyses (with the only exception being a more in-depth experimental analysis of ZNF498 in ciliogenesis, which stands out as comparatively more directly supported by experimental evidence) primarily rely on correlative expression patterns and KZFP ChIP-seq datasets to infer potential biological functions in unrelated cellular contexts. As a result, it remains unclear whether the proposed transcriptional programs are mechanistically linked to the proliferation phenotypes that motivated candidate selection in the first place.

      This issue is evident across multiple sections of the manuscript. For example, the proposed role of ZNF43 in regulating fatty acid metabolism and detoxification pathways is primarily inferred from tissue-level expression correlations. While these analyses focus on genes identified as potential ZNF43 targets, the underlying ChIP-seq datasets were themselves generated under ZNF43 overexpression conditions. Therefore, the current analyses do not establish whether ZNF43 regulates these pathways under physiological expression levels or within a relevant cellular context, nor how such regulation relates to the proliferation defect observed in K562 cells. Moreover, several proposed target genes remain substantially expressed in tissues where ZNF43 expression is not particularly low (e.g. kidney and heart muscle), suggesting that additional regulators are likely involved.

      Similarly, the proposed model of ZNF257-mediated regulation of MAGEA genes during spermatogenesis is intriguing but does not fully account for the expression behavior of all MAGEA family members, particularly MAGEA2B, which displays strong expression in spermatocytes despite high ZNF257 expression. This expression pattern should be acknowledged in the main text and reflected in Fig. 3K. In addition, the labels for MAGEA6 and MAGEA2B in Fig. 3C appear to be inverted. More broadly, the proposed regulatory model is difficult to reconcile with the generally restricted expression pattern of MAGEA genes across adult tissues, as their expression does not appear to consistently correlate with ZNF257 levels outside the germline context.

      Related concerns also apply to the analyses of ZNF498 and ZNF18, where the proposed functions in cilium formation and sperm maturation remain disconnected from the proliferation defects identified in the initial screen. In addition, interpretation of the SCAN-deletion experiments is complicated by the reduced expression levels of the deletion constructs relative to the corresponding full-length proteins, making it difficult to determine whether the observed proliferation phenotypes are pathway-specific or partially driven by differential expression. Finally, while the proteomics results aimed at identifying SCAN-dependent interactors are of interest, several aspects of the experimental design and data analysis remain unclear. In particular, it is not specified whether the experiment was performed in biological replicates or as a single measurement. This is important, as it directly affects how the data can be interpreted and how stringent downstream filtering can be. In the Results section, the authors state that "we identified a set of SCAN-dependent interactors, i.e. proteins that co-immunoprecipitated with the full-length construct but were absent in controls and lost upon deletion of the SCAN domain," which suggests a relatively binary, "presence/absence" filtering strategy. However, this description does not specify whether any quantitative threshold (e.g. enrichment ratio) was applied when comparing full-length constructs to deletion mutants. In contrast, the Methods section states that "proteins lacking signal above background were excluded and proteins were additionally required to show stronger signal in at least one bait condition than in GFP controls based on heatmap clustering (see script)," which instead suggests that a threshold-based criterion was used to define enrichment relative to controls and deletion mutants. If this is the case, the exact criteria and thresholds used for filtering should be clearly stated and consistently reported between the Results and Methods sections. If replicate measurements were not performed, this should be explicitly acknowledged, as peptide-level variability may substantially influence the identification of high-confidence interactors, particularly if the applied cutoffs are not highly stringent.

      Overall, many of the proposed biological functions are currently supported primarily by correlative analyses and would benefit either from additional functional validation or from a more cautious framing of the conclusions.

      Minor comments

      • In the Abstract and in the "Limitations of the study" section, the term "annotation" is used. It would be preferable to specify "functional characterization" instead of "annotation".
      • In the Introduction, there may be a minor citation confusion. Following the sentence: "Characterized by an N-terminal KRAB domain and a C-terminal tandem array of C2H2 zinc fingers, KZFPs primarily target transposable element (TE)-embedded sequences," the cited references are predominantly experimental studies supporting this statement. However, the inclusion of the review "Bruno, Mahgoub and Macfarlan, 2019" appears less appropriate in this context, as it does not directly present ChIP-seq data supporting this claim. More relevant primary studies from the same research area include "Wolf et al. 2020" and "Bruno et al. 2025.".
      • In Fig. 1A, "D10" appears inconsistent with the text and other figures (Fig. 1B, 1G, 1H), which refer to 9 days post-induction.
      • In Fig. S1, there may be a mismatch in the highlighted plate: the zoomed image appears to correspond to the first plate from the top. The correct plate should be highlighted for consistency.
      • In Fig. 1B, there is a typographical error ("K ZFPs" instead of "KZFPs").
      • In Fig. S1E, it is unclear what "other" refers to. Please clarify whether this represents the mean of all remaining KZFPs or a defined subset, ideally in the figure description.
      • In Fig. 2A, readability could be improved by adjusting the layering of points, as the darker dots (in particular the red ones) are currently obscured by lighter ones. Alternatively, removing the outline of the points (which is not transparent) may also improve visibility, but in that case the legend for point size would need to be updated accordingly.
      • In Fig. S2E, "SetDB1" should be corrected to "SETDB1".
      • In Fig. 3B, it is unclear what distinguishes the upper and lower "Diverse REs". A brief clarification in the figure legend would improve interpretability, particularly regarding the transposable element families included.
      • In Fig. S3C, the x-axis labels appear slightly misaligned and shifted to the right.
      • In Fig. 3C, the labels for MAGEA6 and MAGEA2B appear to be inverted.
      • In Fig. 3K, "MAGE3" should be corrected to "MAGEA3".
      • In the ZNF498 section, line 4, the punctuation should be corrected so that the period appears after the figure reference ("promoters (Fig. S1E).").
      • In the final sentence of the ZNF498 section, a noun appears to be missing after "cytoskeleton-dependent," possibly "processes".
      • In the last section of the Results and corresponding figures and their descriptions, "SCAN dependant" should be corrected to "SCAN-dependent".

      Significance

      General assessment

      This study presents a large-scale inducible overexpression platform aimed at systematically exploring the functional diversity of human KZFPs. The screening framework itself represents a potentially useful resource for prioritizing candidate KZFPs for downstream investigation and may be of interest to researchers studying KZFPs, transcriptional regulation, and transposable element biology. A notable strength of the study is the breadth of the screening effort and the attempt to integrate multiple orthogonal datasets to generate functional hypotheses for relatively understudied KZFPs. The more in-depth experimental analysis of ZNF498-mediated ciliogenesis also provides an example of how the platform may be used to identify biologically relevant candidates for further characterization. At the same time, many of the functional conclusions currently remain preliminary and are supported primarily by correlative analyses integrating overexpression phenotypes, expression datasets, and KZFP binding preferences also determined by overexpressing KZFP constructs. In several cases, the proposed biological functions remain insufficiently connected to the original proliferation phenotype used for candidate prioritization. As a result, the current study is best viewed as an exploratory and hypothesis-generating framework rather than a definitive functional characterization of the selected KZFPs. Clarifying these limitations and moderating some of the broader conclusions would substantially strengthen the manuscript.

      Advance

      The principal advance of the work is therefore primarily technical and resource-oriented, providing a scalable experimental framework for systematic KZFP prioritization and downstream functional exploration. While the study does not yet provide extensive mechanistic validation for most proposed functions, it may nonetheless serve as a useful starting point for future investigations into KZFP biology and transcriptional regulation.

      Audience

      The manuscript will likely be of greatest interest to a specialized basic research audience working on KZFPs, transposable element regulation, epigenetic regulation, and transcriptional control. Researchers interested in large-scale functional screening approaches may also find the methodological framework useful.

      Field of expertise: transposable elements, KRAB-zinc finger proteins, epigenetic regulation, functional genomics, genome regulation.

    1. eLife Assessment

      This important study explores whether complex structures that are lost during evolution can re-evolve, which is a long-standing debate in evolutionary and developmental biology. The authors demonstrate that re-evolution can occur if the gene regulatory network that underlies the development of complex traits is maintained. The evidence supporting its conclusions is convincing and the work will be of interest to those studying the evolution and development of complex traits.