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
Reply to the reviewers
We thank the reviewers for your thoughtful assessment of our study and for recognizing its technical advance. We appreciate the reviewers’ comments and suggestions, all of which we have carefully addressed in the revised manuscript.
Reviewer #1
1. In the first section of results, a self-cleaving 2A peptide is inserted upstream of the selection marker to improve selection of correctly tagged genes, as reported in the cited literature. To demonstate this in their PCR-tagging method, the authors tag the HNRNPA1 gene with the mNG-2A-NeoR strategy and find a tagging efficiency of 70% and more, in contrast with the 32% efficiency found with the original method developed by Fueller et al. (2020) (page 5). In support of this conclusion, only a representative image of cells showing mNG and HOECHST signal is presented, with no quantification (Fig1b). As this is a central result of the manuscript that allowed the authors to maintain the 2A-NeoR selection strategy for further investigation of sequence-level data, it is advisable to perform a more thorough quantification of the tagging efficiency by counting (with the help of image analyis software such as Fiji) the mNG+ cells and the total number of cells (HOECHST+), and calculate the % of mNG+ cells across different biological replicates (ideally, cells transfected in different days). A direct comparison (statistical analysis) with the % of mNG+ cells obtained by the original methodology (Fueller et al) would be desirable. In this way, the robustness of the improved tagging efficiency with the optimized method could be also evaluated within one single gene.
Re____sponses
- We have now quantified the fraction of mNG+ cells across three biological replicates for three different genes, using cells transfected in independent experiments. The representative images and quantitative analysis have been added to the revised Figure 1b.
- We also revised the corresponding text in the revised manuscript on page 5 (lines 143-151) to avoid overinterpreting the comparison with the tagging efficiency previously reported by Fueller et al. (2020). The previous experiments were performed in parental HEK293T cells using Zeocin-based selection, whereas the present experiments were performed in HEK293T NeoR-KO cells and G418 selection. Since several experimental parameters were not fully matched, we agree that a direct statistical comparison would not be appropriate. The revised text is:” Using the 2A-linked selection cassette, we tagged HNRNPA1, DDX21, and SLC16A1 in HEK293T NeoR-KO cells (HEK293T cells in which the neomycin-resistance (NeoR) gene linked to the T-antigen has been inactivated), with three independent biological replicates for each gene. … (Figure 1b). Quantification across biological replicates confirmed the reproducibility of the 2A-linked tagging workflow and support its use in subsequent characterization and optimization.”
- In addition, we changed the title of this Results section to “2A-linked selection and characterization of integration outcomes of CRISPR-Cas12a-assisted PCR tagging” to avoid the claim of direct improvement relative to the previous method. Nevertheless, the high and reproducible tagging efficiencies observed here support the use of the 2A-linked selection strategy for subsequent experiments.
- In the following line (page 5), authors talk about "correctly localized" signals, although no clear criterium of "correct" signal localization is explained, nor it is indicated anywhere in the example image. While in the following section(s) it is mentioned that diffuse signal may come from unspecific mNG expression of concatenated PCR cassettes, it would be still important to quantify the "correct" vs "uncorrect" mNG signal (i.e: nuclear vs cytoplasmatic) via image analysis.
__Re____sponse: __Thank you for pointing out the lack of localization definitions. In the revised Figure 1b, examples of diffuse cytoplasmic fluorescence are indicated by white arrows, and the fractions of the two localization categories (“correct” vs “uncorrect”) are quantified. Accordingly, the Results text was revised as “Following G418 selection, fluorescence microscopy showed that most mNeonGreen-positive cells displayed signals consistent with the expected gene-specific subcellular localization, whereas only a small fraction exhibited diffuse mNeonGreen fluorescence (Figure 1b)” on page 5 (lines 146-149). And the criteria used for classification are now described in the Methods and figure legend.
- Of particular significance would be to clarify whether there is a relationship between the % of concatemer product and the presence of diffuse cytoplasmic signal. In several instances, the authors mention their previous finding that diffuse cytoplasmic signal usually results from unspecific expression of the fluorescent reporter within a concatenated PCR cassette. However, from Fig1b there is no evident correlation between % of concatemer reads and % of diffuse cytoplasmic signal. Moreover, Fig2c is cited but no arrows are indicated in the figure in contrast with what mentioned in the text (page 10). Moreover, it is unclear what the authors would like to highlight in the following paragraph when talking about the percentage of diffuse signal (the discrepancy of the concatemers vs diffuse signal? Please make the paragraph more clear). Overall, I think that showing actual microscopy data and image quantification of diffuse vs localized signal would be ideal instead of only mentioning the result in the text. Moreover, a clearer and more detailed discussion about the authors' interpretation of the discrepancy observed between the % of concatemers and % of diffuse cytoplasmic signal is needed, and how it relates to previous findings by Fueller et al. (2020).
Re____sponses
- We apologize for the confusion. In the revised Figure 1b we included the example of diffuse signals and also quantified the fraction of cells with diffuse signal across three biological replicates. We also show example of diffuse signals in revised Fig. S2d. In addition, we revised the paragraphs on page 22 (lines 494-514) to make the description of concatemers and diffuse signals clearer: “Since these reads represent donor-donor junctions rather than intact alleles, their abundance cannot be interpreted directly as concatemer copy number or as the fraction of cells harboring concatemeric integrations.
- In our previous study (Fueller et al., 2020), we proposed that concatemer formation may contribute to diffuse cytoplasmic fluorescence (examples indicated by white arrows in Figure 1b and Figure S2d) through unintended expression of mNG from downstream cassette copies. Although an SV40 PolyA terminator sequence was included downstream of NeoR, the U6 Pol III promoter driving crRNA expression has been reported to mediate Pol II-driven transcription that can contribute to the expression of downstream mNG sequences (Gao et al., 2018; Rumi et al., 2006; Figure 4a). We therefore quantified the fraction of cells displaying diffuse signal at eight loci under control conditions with unmodified donors and compared with the corresponding fractions of concatemer reads, to examine the relationship between concatemer formation and diffuse fluorescence.
- We found that the frequency of diffuse signal varied substantially across loci and did not closely track the fraction of concatemer reads. No diffuse signal was observed at the CLTC, DDX21, and HNRNPA1 loci. Only a few cells (
- It is expectable that the efficiency of PCR-tagging might vary if different preparations of the PCR products are used. The authors should add technical details in the method section that are currently missing regarding kits and protocols used for PCR, PCR purification/extraction, verification of PCR by sequencing (within the "Generation of plasmids and PCR cassettes"paragraph). This includes kits, polymerase type (Q5? Phusion?), extraction protocols (was the PCR gel-extracted?), sequence verification (was the extracted band sent for sequencing prior to perform PCR-tagging?) and so on. In particular, sequence verification by Sanger sequencing would help understand whether some small nucleotide changes were generated by PCR prior to cell transfection and should be added to the discussion part as source of possible mutations.
Responses
- The PCR cassettes were amplified using Velocity DNA polymerase (Bioline), following our previously published protocol (Fueller et al., 2020). As described there, other commercia high-fidelity polymerases like Phusion (Thermo Fisher Scientific), as well as self-purified Pfu-Sso7d DNA polymerase, can also be used.
- Sanger sequencing of PCR cassettes before transfection was not performed, as the cassette sequence is defined by the known template plasmid and the primer sequences. Template plasmid was verified by sanger sequencing and the primers were synthesized and purified by cartridge to ensure high fidelity/purity. Although small nucleotide changes introduced by PCR cannot be excluded as a potential source of mutations, but their overall contribution is expected to be minor given the high-fidelity nature of the polymerases.
- We have now added the experimental details to the Methods on page 33 (lines 817-827):” The PCR cassettes were amplified from either pMaCTag-2A-N07 (mNeonGreen) or pMaCTag-2A-N11 (mScarlet-i) template plasmid based on our previously published protocol (Fueller et al., 2020). In brief, PCR cassettes were amplified using a Velocity polymerase (Bioline) in a reaction mixture containing self-made HiFi buffer (200 mM Tris-HCl pH 8.8, 100 mM (NH₄)₂SO₄, 500 mM KCl, 1% [vol/vol] Triton X-100, 1 mg/ml BSA, 20 mM MgCl₂), 10 mM dNTPs (Bioline), 5 M betaine (Sigma-Aldrich), 200 ng template plasmid, and 10 µM each of M1 and M2 tagging oligos. The following cycling conditions were used: 3 min at 95°C; 30 cycles of 20 s at 95°C, 30 s at 64°C, and 45 s per kb at 72°C; 5 min at 72°C; and 4°C hold. After PCR, 0.4 µl DpnI was added and the reaction was incubated for 1 h at 37°C to digest the plasmid template. PCR products were verified by agarose gel electrophoresis to confirm product size. As our PCR reactions consistently yielded a single specific band of the expected size, column-based purification (Macherey-Nagel, NucleoSpin Gel and PCR Clean-up) was performed without gel extraction.”
- [OPTIONAL] A significant improvement of the study in terms of reproducibility would be to analyze, via microscopy image analysis (NOT by sequencing), the % of tagging efficiency as suggested above (see point 1 under the microscopy section) for the same gene (i.e. HNRNPA1 or other) with biological replicates consisting of cells transfected with PCR reactions prepared in different days, or [OPTIONAL] even analyzing how the tagging efficiency varies with PCR reactions purified with different kits or protocols (example: gel extraction vs PCR elution). Image analysis should include signal quantification with estimation of the % of tagging efficiency and statistical analysis. These experiment, not particularly demanding in terms of reagents and time (~3 weeks in total due to puromycin selection according to the workflow in Fig2b) would greatly help the reader understanding the robustness of tagging efficiency within single genes.
Re____sponses
- As described in our response to the comment 1, we have now quantified tagging efficiency by microscopy for three different genes across three independent biological replicates. For each replicate, cells were transfected with independently prepared PCR products (generated on different days, in some cases by different lab members). Fraction of mNG+ cells was quantified using Fiji, and the statistical analysis has been added to the revised Fig. 1b.
- We agree that systematically comparing different PCR purification methods could provide additional technical information. However, our PCR reactions consistently yielded one specific band of the expected size, we therefore performed a simple column-based purification to minimize the processing time. We have nevertheless clarified the experimental procedures in the Methods on page 34 (lines 824-827): “PCR products were verified by agarose gel electrophoresis to confirm product size. As our PCR reactions consistently yielded a single specific band of the expected size, column-based purification (Macherey-Nagel, NucleoSpin Gel and PCR Clean-up) was performed without gel extraction. All the primers and oligos used for cloning, validation and PCR cassettes are listed in Table S1.”
- The concentrations of pharmacological inhibitors used (stated in the result section at page 8) should be included in the method section (in the "Transfection" section) and in FigS2, along with details of their use in cell culture. "If NHEJ inhibitors were applied, cell culture media containing small molecules (KU0060648, NU7441, VE-822, M3814) were added after 8 hours of transfection and cells were incubated for 16 hours. Two days after transfection, cells were selected with 500 μg/ml neomycin (Sigma-Aldrich)". Was the drug kept in the medium during neomycin selection?
Response: We have revised Fig. S2a to highlight the selected concentrations of inhibitors in red on x-axis. We have also added the concentrations to the Methods section and revised the text accordingly on page 34 (lines 834-837): “Briefly, cells at ~70% confluency in one well of a 24-well plate were transfected with 100 ng of Cas12a helper plasmid (either impLbCas12a (Addgene #137715) or impLbCas12a-hGEM) and 900 ng PCR cassette. After 8 h, the transfection medium was replaced with either fresh culture medium. Where indicated, NHEJ inhibitors (200nM KU0060648, 1μM NU7441, 1μM VE-822, and 500nM M3814) were added at this time, and cells were incubated for an additional 16 hours before the inhibitor-containing media was replaced with fresh culture medium. One day after transfection, cells were expanded to a 6-well plate. Two days after transfection, cells were selected with 500 µg/ml G418 (Sigma-Aldrich).”
- While the authors mention the absence of technical/biological replication within single genes, some graphs are presented with error bars (Fig1F, Fig2d, FigS2b) although no details are given in the text. What replication has been done? What do the error bars represent? What is the sample size?
Re____sponse: Thank you for pointing out this ambiguity. For the samples shown in the original Fig. 1f and Fig. 2d, sequencing libraries were prepared independently in different batches from the same population of tagged cells (from the same transfection), and were sequenced on separate batches. We initially displayed the variation between these technical measurements as error bars. But since they do not represent independent biological replicates, we agree that the error bars could be misleading, so we now show the results from a single sequencing batch for each condition (revised Fig. 1f and revised Fig. 2e). For Fig. S2b, the results of 8h siRNA knockdown were generated from two biological replicate, we have also clarified this point in the revised figure legends.
- (pages 9-10, Result section): "We quantified the ratio of the HDR band intensity to the total intensity of both bands, defining this as "HDR purity". At the HNRNPA1 locus, we observed that both inhibition of c-NHEJ and synchronization of Cas12a increased HDR purity (Figure 3a)." How were the bands quantified to obtain the reported HDR purity %? Which software was used for band intensity quantification? Similarly, the information about how HDR purity was calculated with the anchor-seq data is lacking. These pieces of information are not present in the Method section and should be added.
Response: This information has been included in the Method section on page 36 (lines 903-905): ” The PCR products were analyzed by gel electrophoresis on 1% agarose/TAE, where the upper band corresponds to junctions generated by integration via c-NHEJ and the lower via HR. The intensity of both bands was quantified using Fiji (ImageJ 2.14.0/1.54f). HDR purity was quantified as the ratio of the HDR band intensity to the total intensity of both bands.” and on page 38 (line 966): ” In addition, HDR purity was calculated as the ratio of the HDR / (HDR + Indel) reads. All custom scripts were integrated into a Snakemake workflow, allowing automated execution of the analysis pipeline.”
Minor comments:
- A few typos are present throughout the text. Examples: "because the probability of perfect insertion at both ends of the reporter are very low" (page 2); "With this 'mammalian PCR tagging' method constitutes a" (page 3); "In this work we address these gaps and employ selection markers fused in-frame to the tag of choice (e.g. mNeonGreen) using a 2A peptide spacing the two protein domains ensures that" (page 4); "but increased the off-target frequency for TOMM20 by 0.6-fold? (page 9); "from the wide-type and tagged cells" (page 17); "Together, our improvements outline an robust strategy" (abstract). Some other might be present, please check carefully.
Re____sponse: We apologize for the typos. These have been fixed in the revised manuscript. We also checked to make sure other typos are corrected.
- The Acknowledgement part seems incomplete, and the author contributions are not listed.
Re____sponse: The Acknowledgement part has been revised and author contributions have been added.
- Can the authors explain the concept/meaning of "index hopping" in the text (page 6)?
Response: We have revised the text on page 6 (lines 207-213): “In Tn5-Anchor-Seq, the customized barcode strategy enables the pooling of a large number of samples immediately after Tn5 tagmentation, but this pre-PCR pooling can introduce index hopping — the misassignment of reads to an incorrect sample barcode through cross-amplification. To assess the extent of index hopping from pre-PCR pooling, we designed a dedicated experiment in which reads assigned to barcode-gene combinations not expected in the pool were counted as misassigned, and their fraction relative to total reads was calculated as the index hopping rate. The observed rates ranged from 0.05% to 6%.”
- The FigS2g does not exist (cited at page 8). Did the authors mean f?
Re____sponse: We apologize for the incorrect figure citation. It has been corrected in the revised manuscript.
- There is a discrepancy between Figures and labels regarding FigS3 at page 10: "The PS&Btn modifications on PCR cassettes reduced the frequency of concatemeric reads by 0.5- to 1.7-fold for most genes, while no effect was seen for CANX and SLC16A1 (Figure 4c and Figure S3d). As expected, inhibition of NHEJ with KU, NU, or M3814 decreased the concatemeric integration, while treatment with VE or the use of Cas12a-hGEM increased concatemer frequency at specific loci, such as CANX, DDX21, and SLC16A1 (Figure 4d, 4e, and Figure S3e). Inhibition of MMEJ by knocking down Polθ reduced concatemeric integration at the EMD and SLC16A1 loci but not at the LMNA locus (Figure S3e)." FigS3e does not exist, instead S3d seems to contain what S3e was supposed to show.
Re____sponse: We apologize for these mistakes. We have carefully reviewed the entire manuscript to ensure that all Figures and labels are consistent.
- Figure 5b is never cited in the main text, nor elsewhere in the manuscript.
Re____sponse: We apologize for the typo. Fig 5b was mistakenly referred to as Fig 5a. This has now been fixed.
- In the FigS2c caption, remove the sentence "All PCR cassettes used in this experiment were modified with PS and Btn." since both modified and unmodified PCR cassettes are shown in the panel.
Re____sponse: We apologize for this mistake. The Fig. S2c caption has been revised.
- The scale bar in the microscopy images (Fig1b, Fig5c) is not clearly visible, please fix the thickness.
Re____sponse: Thanks for pointing this out. We have increased the thickness of the scale bars in the revised Figure 1b and Figure 5c.
- In Fig5d, some lanes are highlighted in red although no explanation is given in the figure caption or in the text of what red color lanes represent.
Re____sponse: The red lanes indicate the clones that were subjected for sequencing. The Fig. 5d caption has been revised accordingly.
- In FigS1, different graphs of different panels have the same legend, that is placed only within panel f. The same applies for FigS3. It would help readability to state where to find the legend in the figure caption or to rearrange the panels so that the legend appears next to each graph.
Re____sponse: Thank you for pointing this out. The plots in Fig. S1b, d, and f use the same symbols and colors for the integration outcome categories, with a shared legend next to panel f. We have clarified this in the revised Fig. S1b legend as well as the caption to improve the readability. We also make sure the other figure legends are clear in the revised manuscript.
11.Throughout the schematics used to exemplify the method workflow or cassette architecture there are several abbreviations that are not expanded in the relative figure caption. Examples: Figure 1a (HA, mNG, NeoR, U6, crRNA), Figure 1b (ONT, BC, Tn5) [...]. While some of them can be evinced in the main text, expanding each abbreviation in the respective figure caption can improve the figure readability.
Re____sponse: We have included the abbreviations in the respective figure captions throughout the revised manuscript.
- In Fig1e, several symbols are present in the IGV snapshots whose meaning can only be guessed (what are purple vertical lines? What do numbers and letters indicate and what do the color indicate?). This information should be added as legend in the panel or explained in the figure caption. Similarly, it is not clear what the IGV snapshot indicates, it would be useful to add vertical lines in the IGV snapshot indicating, for example, the 5'-HA/3'-HA portion or something else that could help the reader understand what portion of the cassette is indicated.
Re____sponse: Thank you for the suggestions. We have added vertical lines to the IGV snapshots in the revised Fig. 1e to indicate the approximate positions of 5’-HA junctions. We have also expanded the figure caption to explain the IGV annotations: “(e) Schematics illustrating genomic integration scenarios following PCR tagging and corresponding IGV snapshots showing sequencing-read alignments: i. HDR-mediated integration; ii. NHEJ-mediated deletions and insertions; iii. Off-target integration; iv. Concatemeric integration. Black dashed lines connect each schematic to the genomic interval displayed in the corresponding IGV snapshot. In the IGV snapshots, grey bars represent aligned sequencing reads, whereas mismatches to the reference sequence are shown as colored vertical marks according to the standard IGV nucleotide color scheme (A, green; C, blue; G, brown; T, red). Purple vertical marks denote insertions, with numbers indicating the base-pair length of the insertion. Gaps or thin connecting lines within individual reads indicate deletions, with numbers indicating the size of deletion. Blue vertical dashed lines indicate the approximate positions of the 5’-HA junctions.”
- Please specify that FACS was used only for the seamless tagging experiment, as it's unclear in the text. Were sorted cells seeded as one cell per well, or were they collected in one single tube and then separated? Why is Fig5b showing the sorting of mNG fluorescent cells while the result and method only mScarlet-i is mentioned (possibly wrong figure)?
Re____sponses
- FACS was only used for the seamless tagging experiment, and single cells were seeded as one cell per well. We have added this information to the Methods text on page 35 (lines 877-884): “During seamless PCR tagging, following G418 selection, cells were harvested, washed with PBS…… Single mScarlet-i-positive cells were sorted directly into separate wells of 96-well plates containing……”
- We have also revised Fig. 5b, in which the mScarlet-i was used as the fluorescent protein.
- In FigS2, panels e-f show the same variable (fraction of HDR reads). For EMD locus, however, there is a possible discrepancy on the fraction shown in e (yellow bar, PS&Btn) and f (grey bar, no treatment, modified oligos). Can the authors explain this apparent discrepancy?
Re____sponse: Thank you for pointing this out. Panels e and f (correspond to the revised Fig. 2d and S2e) were generated from different batches of transfection. The difference of EMD values is attributed to batch-to-batch transfection variability. Therefore, comparisons within each panel were performed within the same transfection batch. We have clarified this in Methods text on page 34 (lines 844-847): “Independent experiments shown in separate figure panels were performed in distinct transfection batches unless otherwise indicated. Because absolute editing efficiencies can vary between transfections, treatment effects were evaluated relative to the corresponding controls within each experimental batch. ”
- Looking at some figures, such as Fig2, it is sometimes hard to relate the results to a "bigger picture". For example, Fig2d depicts "% of HDR integrations", while Fig2f shows "Fraction of HDR reads (%)". Besides possible confusion regarding the name of the variable analyzed, the caption of Fig2f explains that the relative panel contains the data in Fig2d. As Fig2d contains % of HDR-mediated integration of unmodified cassettes ranging 60-85%, one would expect to find the same percentages in Fig2f, under unmodified, control cassette; however, on average the fraction of HDR reads in these samples is only 10-20%. How is this possible? Is Fig2f including unmapped reads or reads mapped to off-target loci?
Re____sponse: We apologize for the confusion. We agree that the original description was not reader-friendly. We have restructured Results section 2 (pages 7-9). We now first describe the full spectrum of tagging outcomes across all gene loci and treatment conditions using the metrics “fraction of tagging outcomes”. Then we present in detail the “HDR efficiency” across conditions. In addition, we have clarified these metrics in Methods on page 37 (lines 961-966):” The total read count for tagging outcomes was defined as HDR + Indel + off-target + concatemer reads, whereas genomic integrations (Figure 1c) were defined as HDR + Indel + off-target reads, excluding concatemer reads. Thus, outcome fractions in the holistic analysis (Figure 2d) were normalized to all classified tagging outcomes, while event efficiencies/frequency (HDR, Indel, off-target) were calculated relative to genomic integrations only.”
- The graph in Fig3d displays Log2FoldChange of indel frequencies. To have a better understanding of the data, it is advisable to show an additional graph with "raw" indel frequencies for each gene, as side-by-side bars (similar to the graph in Fig2c). The same applies for the fold change shown in Fig3g.
Re____sponse: Thank you for this suggestion. We have revised the figure to include the additional graphs showing the raw frequencies under the control conditions (revised Fig. 3d, 3g), and moved the Log2FoldChange graph to revised Fig. S3c, S3e. We have also included the side-by-side bars showing the raw frequencies under all conditions in the revised Fig. S3b, S3e.
- [SUGGESTION TO IMPROVE CLARITY] Relative to the comments above (15-16), it would be advisable to add one or more graphs that summarize the read composition, similar to what has been done for Fig1f, for all genes analyzed in the different conditions. For example, it would be advisable to add a supplementary figure with one graph per gene analyzed (8 graphs in total) in the form of the one displayed in Fig2f, with stacked bars instead of boxplots as shown in Fig1f. This will help the reader understanding better and compare side-by-side the effect of the treatments applied in read composition for each gene.
Re____sponse: Thank you for this suggestion. We agree that a holistic overview of the outcome composition is important and helpful to illustrate the overall effects of different conditions. We have therefore showed the full spectrum of tagging outcomes across all gene loci and treatment conditions in the revised Fig. 2d and S2e. We have also revised the results text accordingly on pages 12-13 (lines 333-355).
Reviewer #____2
- The most important shortcoming of this manuscript is the lack of biological replicates. The authors present results from multiple loci, but in this setting, variation will be largely driven by locus-specific effects rather than true experimental variation, which can be substantial for these types of methods. This makes it difficult for readers and potential users to evaluate how robust this pipeline is. The authors must present biological replicates for at least a small number of loci with their finalized method.
o Significant time and resources will be needed.
Re____sponse: Thank you for this valuable suggestion. We agree that biological replicates are critical for demonstrating the reproducibility and robustness of our workflow. We have therefore revised Fig. 1b, to include quantification of the fraction of mNG+ cells across three biological replicates for three different genes. These replicates were obtained from independent transfection experiments and were subjected to statistical analysis. The text on page 5 (lines 143-151) has also been revised accordingly.
- The authors characterize clonal genotypes following seamless PCR tagging at two loci. They report that 2/6 isolated clones of HNRNPA1 have the desired edit, while 3/3 clones of SLC16A1 have the desired edit. However, the small number of amplicons sequenced makes it difficult to draw robust conclusions about clonal genotype outcomes. The authors should perform additional amplicon sequencing to strengthen these conclusions. In addition, they should sequence the incorrectly sized products generated from the junction-PCRs to gain further insight into alternative genotypes and incorrect editing that may occur with this method.
o Additional sequencing can easily be accomplished if cell lines are still available. An estimated time of about 1-2 weeks for cell thawing, genome extraction, amplicon sequencing, PCR purification, and outsourced Plasmidsaurus sequencing.
Re____sponse: Thank you for this important comment. We agree that sequencing more clones would robustly determine the frequencies of different genomic integration outcomes following seamless PCR tagging. However, our intention was not to quantitatively mapping full spectrum of editing outcomes, rather demonstrate that the two-cassette strategy can generate seamless integrations at different gene loci. We have therefore revised the Results text on pages 26 (lines 603-614) to avoid overclaiming regarding the genotype frequencies.In addition, we have revised the corresponding Discussion text on page 30 (lines 704-713) accordingly.
- It is not obvious how the authors were able to accurately quantify concatemeric integrations with the Tn5-anchor-seq tool. Specifically, the schematic in Figure 4a makes it clear that concatemers have multiple binding sites for the PCR primer. If these binding sites give rise to multiple products from the same template, the quantification would be inaccurate. The authors should comment on their considerations in suppressing the formation of secondary products, such as the use of a non-strand displacing polymerase, and in ensuring that their quantifications are accurate, including any bioinformatic steps.
o The authors can clarify methods and adjust results presented in the manuscript in 1-2 days.
Re____sponse: Thank you for pointing this out. The scenario where a single template generates multiple PCR products would require amplification of an intact concatemer. However, in our Tn5-anchor-seq, Tn5 tagmentation precedes PCR so that concatemers are fragmented before PCR. Primer sites from a N-copy concatemer therefore reside on separate tagmented fragments, which are amplified independently. So, the concatemer-derived read counts scale with the amount of concatemeric sequence present in the library. Potential secondary products are further minimized by the nested two-round PCR design, biotin-streptavidin enrichment of cassette-primed amplicons, limited first-round PCR (15 cycles), hot-start polymerase, and size selection of products > 500bp. In bioinformatics analysis, each read is counted once after demultiplexing, mapping, filtering and removal of secondary and supplementary alignments. Importantly, this metric reports the fraction of concatemeric sequence, rather than the number of concatemeric integration events per allele. We have clarified the interpretation of concatemer reads in the revised manuscript on page 22 (lines 494-497).
- The authors begin the results section by describing how their revised tagging cassette improves selection efficiency. It appears that they compare their newer data to data that were obtained for their previous manuscript, which were likely obtained at different points in time, potentially by different individuals, and using different sets of reagents, which makes direct comparisons impossible. It would be beneficial to include at least one side-by-side experiment to compare the authors' new method with their previous method to conclusively demonstrate any improvement.
Re____sponses
- We agree that our original text overinterpreted the comparison with the tagging efficiency previously reported by Fueller et al. (2020). The previous experiments were performed in parental HEK293T cells using Zeocin-based selection, whereas the present experiments were performed in HEK293T NeoR-KO cells using G418 selection. Since several experimental parameters were not fully matched, we recognize that a direct quantitative comparison is not be appropriate.
- We have therefore removed the sentence “compared to 32% when employing a cassette in which the selection marker was inserted as a separate gene (Fueller et al., 2020). ” in the revised Results text on page 5. Instead, we now report the quantification across biological replicates to show the reproducibility of the workflow using the new 2A-linked selection cassette, thereby supporting the use of the 2A-linked selection strategy for subsequent experiments.
- In addition, we changed the title of this Results section to “2A-linked selection and characterization of integration outcomes of CRISPR-Cas12a-assisted PCR tagging” to avoid the claim of direct improvement relative to the previous method.
- Data and method presentation: For much of the manuscript, the authors provide detailed descriptions of experiments and methods. However, the authors should provide more detail in the figure legends about the concentration and timing of chemical inhibitor treatments.
Re____sponse: Thank you for pointing this out. We have added the concentrations and treatment details to the Methods section and revised the text on page 34 (lines 834-837).
6.Replication and statistical analysis: While the authors mention their lack of replicates, and reason multiple loci and conditions compensate for this, replication of their experiments is still necessary in order to draw meaningful conclusions of their data (see comment above).
Re____sponse: As described in the response to the first comment, we have revised Fig. 1b, to include quantification of the fraction of mNG+ cells across three biological replicates for three different genes. The text on page 5 has also been revised accordingly.
Minor Comments:
- The authors clearly identify and cite the prior work on which their method is based.
Re____sponse: We thank you for this kind comment.
- The discussion focuses on the findings and limitations of the authors' specific study. Readers would benefit from a conceptual comparison of the strengths and weakness of their method to other similar technologies in this space, which would allow researchers to make better informed decisions about the appropriateness of this method for their own work.
Re____sponse: Thank you for this suggestion. We have added a paragraph in Discussion on page 31 (lines 718-727) to compare with similar technologies.
- The text and line of the scale bar on all microscopy images needs to be enlarged for readability.
Re____sponse: All the scale bars of the microscopy images have been enlarged in the revised manuscript.
- Text has inappropriate question mark in "... off-target frequency for TOMM20 by 0.6 fold?".
Re____sponse: The question mark has been removed.
- Fold-changes in the text are described in an unconventional and unintuitive manner. For example, the authors note a 0.2 to 0.7-fold reduction in off target frequency contrasted against a 0.6-fold increase, which appear to be two different interpretations of fold-changes. These results might be more accurately represented by using a 1-fold change as the baseline for no change. Alternatively, the authors could provide this measure as percent decreases and increases.
Re____sponse: Thank you for pointing this out. We have revised the text to express decreases and increases as percent changes relative to the control condition, while retaining log2 fold-change values in the heatmaps. The revised text show on lines 436, 447, 448, and 517.
- The authors mention that an example of diffuse cytoplasmic fluorescent signals observed in microscopy can be seen in Figure 2c, which are labeled with white arrows. A figure with this description is not present in the manuscript.
Re____sponse: We apologize for the confusion. In the revised Fig. 1b and S2d, examples of diffuse cytoplasmic fluorescence are indicated by white arrows. Accordingly, the Results text was revised on page 22 (lines 498-500).
- Schematics in Figure 5 show mNeonGreen, however, Figure 5c is labeled as mScarlet. The Sauthors should clarify which fluorescence tag was used.
Re____sponse: We have revised Fig. 5a by replacing “mNG” with the generic word “tag”, as this panel is a general schematic to compare the two protein-tagging strategies. We also corrected Fig. 5b and 5d to indicate mScarlet-i, which is the fluorescent tag used in this part.
- Figure S3e is mentioned in the text, however, the supplemental material does not contain a 3e.
Re____sponse: The typo has been corrected in the revised manuscript.
- The authors have developed a fantastic pipeline to quantify tagging outcomes holistically, and in Figure 1f they visualize the results of the pipeline in an intuitive and informative manner. However, for the remainder of the manuscript, they proceed not to use this visualization and instead focus on the effects of their experimental parameters on individual outcomes or a ratio of two outcomes ("HDR purity"). This choice makes it difficult to clearly evaluate the effects of the parameters on all outcomes, and the "HDR purity" metric, although useful to compare against the gel-based readout, is difficult to interpret. We would recommend including the holistic visualization from Figure 1f for all experiments, in the supplement if necessary, for example as a large grid that depicts the frequency of all outcomes for each experimental condition by locus. This will make it much easier for readers to evaluate how a given experimental parameter affects the frequency of different editing outcomes.
Re____sponse: Thank you for this important suggestion. We agree that a holistic overview of the outcome composition is important to illustrate the overall effects of different conditions. We have therefore restructured the Results section 2 by first showing the full spectrum of tagging outcomes across all gene loci and treatment conditions in the revised Fig. 2d and S2e. We then described the HDR efficiency in the revised Fig. 2e, 2f, and S2g. The corresponding text has also been revised accordingly (pages 12-13, lines 333-355). We have also revised the description about HDR purity results on page 17 (lines 414-423) to make it clearer.
- While the authors detail their rationale in selecting loci to tag in a previous publication, it is important to include the information in the current manuscript as well.
__Re____sponse: __We have included this information in the revised manuscript on page 11 (lines 305-311).
- The authors present many figures that look similar but display different quantities. To improve readability, they should use more descriptive figure titles and axis labels so that the main findings and quantities are immediately clear. This is particularly needed in Fig.2d, where the variable plotted is not clear.
Re____sponse: We have restructured Figure 2 to provide a more holistic overview of the full spectrum of tagging outcomes. The original Fig. 2d has been revised as Fig. 2e with the y-axis labeled “HDR efficiency” and the condition indicated as “Unmodified, no treatment” on top. We have also revised the axis labels and results description on page 12 (lines 349-355) for greater clarity.
- Two aspects of the final method may (or may not) represent limitations that would be helpful for the authors to discuss. First, the authors use a Cas12a with expanded PAM range, which has the advantage of increasing the number of available cut sites but also has the potential to lead to more off-target cutting, which can cause secondary mutations in the selected cell populations. Second, tagging the 3′ UTR may lead to differences in transcript and protein stability (especially due to the remainder of the 2A site). It would be helpful for the authors to include their perspectives on these issues and if they took efforts to address them.
Re____sponses
- Thank you for raising these two important points. Regarding the expanded-PAM Cas12a variant, we have added this point explicitly to the limitations paragraph in Discussion on page 32 (lines 757-759): “The use of an expanded-PAM Cas12a variant may also alter off-target activity relative to wild-type Cas12a (Kleinstiver et al., 2019). Thus, comprehensive off-target profiling is needed in applications where genomic integrity is critical. ” Systematic off-target characterization was beyond the scope of the current study, but we consider this an important direction for future work, particularly for therapeutic applications.
- Regarding the 2A site remainder and 3’-UTR effects, we partially address the concerns about protein stability/integrity by tagging mNG-FKBP degron at four loci and confirmed efficient inducible depletion of the tagged protein following degron induction (revised Fig. S4). This demonstrates that the integrated fusion proteins are expressed and that the degron domain remains structurally accessible. However, we acknowledge that direct measurement of endogenous protein activity will be important in future studies. We have therefore added discussion of this point in limitations paragraph on page 32 (lines 760-764).
Reviewer #____3
- All optimization experiments are performed in HEK293T-N cells, which are aneuploid/polyploid, permissive for genome editing and have atypical DNA repair dynamics. While acceptable for method development, this limits generalizability, especially given the strong locus and pathway dependent effects observed. DNA repair pathway balance (HDR vs NHEJ vs MMEJ, etc.) can be cell-type specific. The improvements reported here may not generalize to non-transformed cells, stem cells, cancer cell lines, or other therapeutically relevant cell types. Selecting a small, representative gene set (2-3 genes) and performing controlled tagging experiments across a panel including cancer cell lines, non-transformed cells (like hTERT immortalized RPE-1 cells), and embryonic stem cells (H1 or H9) could provide readers with a clearer view of integration and tagging efficiencies across diverse cellular contexts. Additionally, assessing the effects of pathway modulation in these cell types would help determine how broadly generalizable these strategies are increasing the impact of the method detailed herein.
Re____sponses
- Thank you for raising this important point. We fully agree that cell-type generalizability is very important. HEK293T NeoR-KO cells were selected for systematic optimization because their high editing permissiveness enables robust detection of differences across conditions and loci, rather than because they are representative of the broader range of cell types.
- To assess the portability of our workflow, we tested it in two additional cell lines. First, we applied the optimized standard PCR tagging protocol — using PS- and Btn- modified donors in combination with 200 nM KU0060648 — to HeLa cells. Tagging of HNRNPA1, DDX21, and TOMM20 was evaluated by fluorescence microscopy and flow cytometry. The fraction of mNG+ cells at expected localization were comparable to the observations in HEK293T-N cells. Second, we applied the seamless tagging strategy to HepG2 cells at HNRNPA1 locus and confirmed the comparable tagging efficiency via fluorescence microscopy. All these data are now summarized in Supplementary Note 1.
- We note that sequencing-based integration outcome characterization was not performed in either HeLa or HepG2 cells. We have clarified this in the revised Discussion on page 32 (754-757).
- While the manuscript presents a clear and well-supported optimization of CRISPR-Cas12a assisted PCR tagging, it would benefit from a more explicit contextual comparison to other contemporary endogenous tagging approaches, such as Cas9-mediated HDR with plasmid donors. Ideally, a single or small panel of genes could be compared to assess the benefits and drawbacks of each technique. Even a conceptual discussion outlining relative strengths, limitations, and use cases would help readers better position the proposed method within the broader genome-editing landscape. Such a comparison would clarify when PCR tagging is the most appropriate choice and where alternative strategies may offer advantages, thereby enhancing the interpretability and practical impact of the study.
Re____sponse: Thank you for this suggestion. We have revised the Discussion to include a dedicate paragraph for comparison with other similar technology on page 31 (lines 718-727).
- Lack of functional validation beyond localization. Localization is used as the primary functional readout. Including limited functional validation (e.g., protein abundance, stability, or activity) for a representative target would strengthen the claim that optimized conditions preserve protein function, since donor modifications and DNA repair pathway inhibition could affect protein dosage or integrity. Alternatively, explicitly acknowledging this limitation would better frame the scope of the functional conclusions. At a minimum this needs to be discussed carefully.
Re____sponses:
- Thank you for this important suggestion. We agree that subcellular localization alone provides limited evidence for functional protein integrity. Thus, we applied mNG-degron tagging at HNRNPA1, DDX21, CANX, and KIF11 loci to confirm protein integrity beyond localization. Following degron induction, efficient depletion of the mNG-tagged proteins was confirmed by microscopy, flow cytometry, and western blot. The combination of correct subcellular localization and successful inducible degradation indicates that the tagged fusion proteins are expressed and structurally intact at the endogenous locus.
- The results are presented in revised Fig. S4, and a brief description has been added to the Results section on page 23 (lines 530-539).
- We noted that these results are not full functional validation, the degradation efficiency does not directly show whether native protein activity is preserved. We have revised the Discussion on page 32 (lines 764-767) to acknowledge this limitation.
- The effects of donor cassette size and sequence composition on integration efficiency and concatemer formation are not systematically explored. Because larger or more complex payloads may engage repair pathways differently, a brief discussion of how donor size might influence the reported improvements would help contextualize the results. If extra curious, a single reliably editable target gene could be used to assess and explore the effect of varying payload composition.
Response: Thank you for this suggestion. The mNG-degron cassette described in the revised Fig. S4 and also our response to Major point 3 is larger and compositionally more complex, as it contains an additional degron sequence and an HA epitope tag. With this degree of payload size expansion, we didn’t observe any obvious reduction in performance with respect to PCR amplification or tagging efficiency assessed via microscopy. However, we acknowledge that the effects on HDR efficiency and concatemer formation were not evaluated specifically. We have now discussed this limitation in the revised manuscript on page 31 (lines 728-735).
Minor ____Comments:
- The frequency of concatemeric integrations is inferred from junction-centric Tn5-Anchor-Seq data rather than direct measurement of intact alleles, which limits the ability to determine concatemer copy number, allelic structure, or the presence of complex rearrangements within a single genomic locus. While this limitation is acknowledged, it should be more clearly emphasized when interpreting the magnitude of concatemer suppression and improvements in genomic homogeneity, particularly given the polyploid nature of HEK cells. Clarifying that junction counts do not directly equate to allele counts would help prevent overinterpretation of these results and better frame the conclusions drawn from the sequencing data.
__Response____: __Thank you for this suggestion. We have revised the corresponding text on page 22 (lines 495-497) to avoid overinterpretation of the concatemer results.
- Inhibition of DNA-PKcs emerges as one of the most effective strategies for increasing HDR-mediated integration and suppressing undesired byproducts in this study. However, the broader practicality of this approach is not explored in depth. DNA-PKcs inhibitors can introduce cytotoxicity, alter cell-cycle progression, or induce stress responses that may be incompatible with certain experimental systems or long-term phenotypic analyses. A more explicit discussion of these limitations, including guidance on dosing, duration of treatment, and scenarios in which donor modification alone may be preferable, would help readers/researcher assess the applicability of this strategy across different cell types and use cases.
__Response____: __We have expanded the discussion of DNA-PKcs inhibitors in the revised manuscript on page 29 (lines 667-670).
- The study focuses exclusively on C terminal tagging, and the applicability of the optimized workflow to N terminal or internal tagging configurations is not examined. Because tagging position can substantially influence repair outcomes and expression, a brief discussion of whether the observed improvements are expected to generalize beyond C terminal insertions would help clarify the scope of the method.
__Response____: __We have expanded the limitation paragraph in Discussion on page 31 (lines 736-743).
- Clarify the rationale for inhibitor concentrations in the main text, not only in supplement.
__Response____: __We have revised the text on page 11 (lines 295-297).
- Define all abbreviations at first use (e.g., c-NHEJ, MMEJ) for clarity.
__Response____: __The abbreviations have been defined in the revised manuscript.
- Explicitly state sample sizes for microscopy-based quantifications.
__Response____: __This has been addressed.
7.Clarify how index hopping was quantified and controlled in pooled sequencing analysis.
__Response____: __We have clarified this in the revised manuscript on page 6 (lines 207-213).
- Briefly describe computational resource requirements for the analysis pipeline.
__Response____: __Thank you for this suggestion. We have included this information in the Methods text on page 38 (lines 971-975).
- State whether the "seamless" tagging method was tested with donor modifications and inhibitors (Assumption is no from your drawings in Figure 5).
__Response____: __We apologize for the missing information. Seamless tagging was performed using PS&Btn-modified donor in combination with KU treatment, as this was the best-performing condition identified in the standard PCR-tagging experiments. We have now clarified this in revised Fig. 5b and Results text on page 26 (lines 593-595).
- Clarify whether donor-end modifications affect PCR yield or fidelity. Also, discuss cost implications of biotinylated/modified primers for scalability
__Response____: __We did not observe any obvious reduction in PCR performance when using modified primers, although PCR yield and fidelity were not quantified as separate endpoints. We have therefore clarified this as well as the primer cost in Discussion on page 29 (lines 660-664).
- State whether sequencing depth was normalized across conditions.
__Response____: __Thank you for pointing this out. Sequencing depth was not explicitly normalized across conditions as all reported metrics are per-sample fractions. However, we ensured that sequencing depth was sufficient for all samples, such that fraction estimates were not biased by low coverage. This is now stated in the Methods on page 38 (lines 966-969).
- Abstract: replace 'an robust strategy' with 'a robust strategy'.
__Response____: __This has been fixed.
- Results section: ensure consistent use of 'Cas12a-hGEM' vs 'hGEM-Cas12a'.
__Response____: __It’s now consistently “Cas12a-hGEM” in the revised manuscript.
- Figure 3g: remove stray question mark after '0.6-fold?'.
__Response____: __It has been removed.
- Methods section: standardize units (µM vs uM) throughout.
__Response____: __We now use µM throughout the revised manuscript.
- Figure references: inconsistent capitalization of figure labels (Figure vs Fig.).
__Response____: __We now use consistently “Figure” in the revised manuscript.

