基因编辑脱靶效应X基因验证论文

基因编辑脱靶效应X基因验证论文一.摘要

基因编辑技术CRISPR-Cas9作为性的分子工具，在遗传疾病治疗、作物改良及基础生物学研究等领域展现出巨大潜力。然而，脱靶效应作为其固有的技术局限，可能引发非预期基因序列的修饰，从而对生物体健康产生潜在风险。本研究聚焦于某基因编辑案例中脱靶效应的验证，通过结合高通量测序（HTS）和生物信息学分析，系统评估了CRISPR-Cas9系统在目标基因编辑过程中的脱靶位点分布及其影响。研究选取了人类细胞系作为实验模型，利用设计优化的sgRNA对特定基因进行编辑，并通过双重测序技术（ddPCR）精确检测目标序列及潜在脱靶位点的修饰情况。结果表明，尽管主要脱靶事件集中于基因距离较近的区域内，但在远端染色体区域也检测到低频脱靶修饰，这些位点与已知功能性基因存在高度重叠。进一步的功能验证实验揭示了部分脱靶位点可能通过引入移码突变或影响转录调控元件，对细胞表型产生显著干扰。研究构建的脱靶风险评估模型，结合sgRNA序列特征与基因组结构信息，可有效预测潜在高风险脱靶位点。该结果不仅为基因编辑脱靶效应的精准验证提供了实验依据，也为优化CRISPR-Cas9系统设计、降低临床应用风险提供了重要参考。

二.关键词

基因编辑；CRISPR-Cas9；脱靶效应；高通量测序；生物信息学分析；功能验证

三.引言

基因编辑技术，特别是CRISPR-Cas9系统，自问世以来便以其高效、便捷和相对经济的特性，深刻改变了生物学研究和遗传疾病治疗的前景。该技术利用向导RNA（gRNA）引导Cas9核酸酶精准识别并结合目标DNA序列，通过非同源末端连接（NHEJ）或同源定向修复（HDR）等机制实现基因的插入、删除或替换，为攻克镰状细胞贫血、β-地中海贫血等单基因遗传病提供了前所未有的可能性。与此同时，基因编辑在作物培育、畜牧业改良以及基础生命科学问题的探索中同样展现出巨大的应用价值，例如通过编辑关键基因提高作物产量、抗逆性或营养价值，或是在模式生物中解析基因功能网络。据相关统计，全球范围内每年有数百项涉及CRISPR-Cas9的科研论文发表，相关专利申请数量持续攀升，技术迭代速度惊人，其应用范围正从实验室研究逐步向临床试验乃至商业应用拓展。

然而，伴随着基因编辑技术的飞速发展，其固有的局限性，尤其是脱靶效应（off-targeteffects,OTEs），日益成为限制其安全性和有效性的关键瓶颈。脱靶效应是指在基因编辑过程中，Cas9核酸酶除了在预设的目标位点外，还意外地在基因组中其他存在高度相似序列的位点进行切割，引发非预期的基因修饰。这些脱靶位点的存在可能带来多种不良后果：首先，脱靶切割可能导致基因序列的插入突变或移码突变，进而产生有功能的毒性蛋白或破坏关键基因的功能，引发细胞凋亡、基因组不稳定性甚至癌症等严重问题。其次，即使脱靶切割不直接导致功能失活，也可能通过影响染色体重排、基因表达调控元件或非编码RNA等功能区域，间接干扰细胞的正常生理活动。此外，脱靶效应的检测和评估一直是基因编辑领域面临的技术挑战，由于传统测序技术难以全面覆盖整个基因组，且脱靶事件通常以低频发生，因此早期研究往往低估了脱靶效应的普遍性和严重性。随着测序技术的进步，特别是高通量测序（High-ThroughputSequencing,HTS）和生物信息学分析方法的引入，对脱靶位点的检测能力显著增强，使得研究人员能够更系统、更深入地评估特定基因编辑系统的脱靶谱。

尽管近年来在降低脱靶率方面取得了显著进展，例如通过优化sgRNA设计、开发高保真Cas9变体（如HiFiCas9）以及引入脱靶效应检测和修复策略等，但完全消除脱靶效应仍是当前面临的技术难题。因此，建立准确、可靠的脱靶效应验证方法，并对其进行深入的功能分析，对于确保基因编辑技术的安全应用至关重要。目前，脱靶效应的验证主要依赖于生物信息学预测和实验验证相结合的策略。生物信息学预测通常基于算法分析sgRNA与基因组序列的相似性，筛选潜在的脱靶位点。然而，预测结果往往存在假阳性和假阴性的问题，预测的脱靶位点未必全部被实验验证，而未被预测的位点也可能实际发生编辑。因此，后续必须通过实验手段对预测的高风险脱靶位点以及关键功能基因区域进行验证，以确认其是否真正发生修饰以及修饰的类型。实验验证方法主要包括末端修复测序（T7E1assay）、酶切法（如Surveyorassay）、数字PCR（dPCR）和HTS等。T7E1和Surveyorassay操作相对简单，但灵敏度较低，难以检测低频脱靶事件。dPCR技术具有高灵敏度和精确性，特别适用于检测单个碱基位点的突变，能够有效鉴定出由NHEJ产生的移码突变。而HTS则能够一次性检测基因组中成千上万个位点的编辑情况，是目前评估脱靶谱广度和深度的主流方法，但需要复杂的生物信息学流程进行数据处理和位点鉴定。

在本研究中，我们针对一个具体的基因编辑案例，系统性地开展了脱靶效应的验证工作。该案例涉及利用CRISPR-Cas9系统对某个与人类疾病相关的基因进行编辑，旨在修正其致病突变。在初步的生物信息学预测基础上，我们采用了HTS结合dPCR验证的策略，旨在精确绘制该基因编辑系统的脱靶谱，并深入探究潜在脱靶位点的功能影响。研究重点在于明确脱靶位点的分布特征，评估其编辑频率和潜在的生物学效应，特别是关注那些距离目标位点较远、生物信息学预测可能被忽视的位点。同时，通过功能验证实验，我们试揭示脱靶修饰对细胞表型或基因功能的具体影响机制。本研究的意义在于，一方面，通过对特定案例的脱靶效应进行详尽验证，为该基因编辑方案的安全性评估提供了关键实验数据；另一方面，所采用的研究方法和分析策略，包括HTS数据的深度挖掘和结合dPCR的精准验证，可为其他基因编辑研究提供参考，有助于提升脱靶效应检测和验证的整体水平。最终，本研究旨在为优化CRISPR-Cas9系统的设计与应用，降低脱靶风险，推动基因编辑技术在临床和农业等领域的安全、有效转化提供科学依据。研究核心问题是：在特定的基因编辑案例中，CRISPR-Cas9系统产生的脱靶效应的真实范围和程度如何？这些脱靶位点是否具有生物学功能？它们对细胞表型有何具体影响？基于这些问题的回答，我们假设：尽管生物信息学预测可能识别出主要的脱靶区域，但通过综合运用HTS和dPCR验证，可以发现包括远端染色体区域在内的更广泛脱靶位点，部分脱靶位点的修饰可能对细胞功能产生可检测的负面影响。

四.文献综述

CRISPR-Cas9基因编辑技术自2012年其潜力被公开以来，便迅速成为生命科学研究领域的核心工具。其基本原理是利用一段与目标DNA序列互补的向导RNA（gRNA）引导Cas9核酸酶识别并结合特定的基因组位点，随后通过非同源末端连接（NHEJ）或同源定向修复（HDR）途径实现基因序列的修改。NHEJ途径是细胞主要的DNA双链断裂修复方式，其特点是效率高但容易引入随机的插入或删除（indels），从而可能导致目标基因失活，这是CRISPR-Cas9在基因功能研究和遗传疾病治疗中最常用的机制。然而，NHEJ的随机性也意味着它极易在基因组中其他存在高度相似序列的位点发生非预期的切割，即脱靶效应，这是限制其安全应用的主要障碍。

早期关于CRISPR-Cas9脱靶效应的研究主要集中在生物信息学预测层面。研究发现，sgRNA的设计是影响脱靶率的关键因素。sgRNA的种子区域（seedregion，通常指前17-20个核苷酸）与基因组序列的匹配度越高，发生非特异性结合和切割的概率也越大。因此，众多研究致力于开发预测脱靶位点的算法。例如，SangerInstitute开发的CHOPCHOP，以及后续的Cas-OFFinder、CRISPR-Ranger、COSMID等工具，通过比较sgRNA种子区域与全基因组序列的相似性，对潜在的脱靶位点进行评分和排序。这些预测工具在不同程度上提高了识别高风险脱靶位点的能力，为实验设计和风险评估提供了重要参考。然而，预测与现实的差距依然存在。一方面，许多预测算法主要关注种子区域的匹配，而忽略了sgRNA的3'末端区域以及PAM位点的序列特征对脱靶的影响。另一方面，算法难以准确评估非目标位点的实际切割效率，因为脱靶事件的发生频率可能极低，且受细胞类型、基因组背景等多种因素影响。此外，预测往往低估了PAM序列附近非完美匹配位点可能发生的切割。

随着测序技术的飞速发展，特别是高通量测序（HTS）技术的应用，研究者开始能够系统性地实验验证CRISPR-Cas9的脱靶谱。早期的研究主要通过末端修复测序（T7E1assay）或限制性酶切结合凝胶电泳（如Surveyorassay）来检测目标附近区域的脱靶indels。这些方法相对简单、快速，但灵敏度有限，且难以检测远距离或低频的脱靶事件。随着NGS技术的普及，基于测序的验证方法成为主流。研究者通常对细胞裂解液进行纯化，提取编辑后的DNA，然后对整个基因组或目标区域进行重测序（re-sequencing），通过与野生型参考序列比较，鉴定出新的indels位点。这类研究揭示了脱靶效应的普遍性，证实即使是设计精良的sgRNA，在多种生物模型和细胞类型中也可能检测到脱靶事件。例如，一部分研究在人类细胞中鉴定出主要脱靶位点集中在距离目标基因5kb以内的区域，这些位点通常与目标位点具有高度序列同源性。另一部分研究则发现了更远距离的脱靶事件，尤其是在基因组重复序列区域或近端着丝粒区域，这些区域由于序列复杂性或特殊结构，难以被常规的生物信息学预测方法有效识别。

对脱靶位点功能性的研究同样重要。早期观点认为，大部分脱靶位点的编辑频率极低，且发生在非关键基因区域，因此其生物学影响可能可以忽略不计。然而，后续研究开始关注低频脱靶事件可能产生的累积效应，尤其是在长期培养或体内环境中。更重要的是，一些研究在距离目标位点较远的位置检测到了可能具有生物学意义的脱靶位点。例如，有研究发现某个用于基因敲除的sgRNA在远端染色体上引发了影响转录启动子或增强子区域的indel，导致下游基因表达水平发生显著改变。另一项研究则报道，一个旨在编辑免疫相关基因的sgRNA在另一个基因的3'非编码区引入了移码突变，干扰了非编码RNA的功能，进而影响了免疫细胞的分化和活化。这些发现表明，即使是低频、远距离的脱靶事件，也可能通过影响基因表达调控、非编码RNA功能或染色体重排等机制，对细胞表型和生物学过程产生不可忽视的干扰，甚至可能引发致癌风险。因此，对检测到的每一个脱靶位点进行功能验证变得至关重要。

近年来，降低脱靶效应的努力主要集中在改进Cas9核酸酶和sgRNA设计两方面。高保真Cas9变体（HiFiCas9）的开发是重要进展之一。通过定向进化或蛋白质工程改造，研究人员获得了如eSpCas9-HF1、SpyCas9-V1、Cpf1等具有更高等位基因特异性（targetsitespecificity）和更低非特异性切割活性的Cas9变体。这些变体在多种实验体系中都显示出显著降低脱靶率的特性。然而，即使是高保真变体，也并非完全消除脱靶，只是将脱靶频率降低至可接受的临床试验标准范围内。另一方面，sgRNA优化策略也在不断发展。除了传统的基于序列相似性筛选，研究者开始考虑sgRNA的二级结构、在细胞内的稳定性、与Cas9的相互作用动力学等因素。一些研究利用机器学习模型来预测sgRNA的脱靶率和编辑效率，从而指导更优sgRNA的设计。此外，双重或多重gRNA的设计策略也被尝试用于提高编辑特异性，通过同时靶向两个邻近位点，只有在两个位点都被成功编辑时才激活标记基因，从而筛选出靶向单一脱靶位点的细胞。

尽管在降低脱靶和验证方法方面取得了长足进步，但脱靶效应的全面治理仍面临诸多挑战和争议。首先，关于脱靶效应的“可接受”阈值是什么，目前尚无统一标准。在基础研究阶段，较低的脱靶率即可满足需求；但在临床应用中，要求则更为严格，通常期望脱靶率低于1×10^-6。然而，如何精确测量和报告脱靶率，特别是远距离、低频的脱靶事件，仍然是一个难题。其次，现有的脱靶验证方法各有优劣，单一方法往往难以全面覆盖。HTS可以提供广泛的脱靶信息，但数据解读复杂，且可能存在假阳性；dPCR可以精确定位和定量特定位点的突变，但通量有限，难以系统性评估整个基因组。如何整合不同方法的优势，建立更全面、更可靠的脱靶验证体系，是当前研究的重要方向。第三，关于脱靶效应的长期影响，尤其是体内环境中的累积效应和潜在致癌风险，仍需更多深入的研究。目前大部分研究集中在体外细胞系，对于脱靶效应在动物模型乃至人体内的真实情况了解有限。第四，不同生物类型、细胞类型和基因组背景对CRISPR-Cas9的脱靶特性有何影响，这些因素如何相互作用，也需要更系统的研究。最后，如何在保证高度特异性的同时，维持CRISPR-Cas9系统的高效编辑能力，也是一个持续存在的平衡难题。针对这些挑战和争议，未来的研究需要更加注重跨学科合作，结合生物信息学、生物化学、细胞生物学、遗传学和动物模型等多种手段，对CRISPR-Cas9的脱靶机制、影响和治理策略进行更深入、更全面的探索。

五.正文

1.研究设计与方法

本研究旨在系统评估特定基因编辑案例中CRISPR-Cas9系统的脱靶效应，并对其进行实验验证与功能初步探究。研究主要分为以下几个阶段：sgRNA设计与优化、细胞系构建与基因编辑、脱靶位点高通量测序（HTS）分析、关键脱靶位点数字PCR（dPCR）验证以及潜在脱靶效应功能影响分析。

1.1sgRNA设计与优化

本研究聚焦于一个与特定遗传疾病相关的基因（记为GeneX）的编辑。首先，利用生物信息学工具（如CRISPR-Ranger和CHOPCHOP）筛选出GeneX基因座周围100kb范围内的潜在PAM位点。基于这些位点设计了初始sgRNA库，其中包含针对GeneX关键突变位点的sgRNA，以及针对其上下游调控区域和远端潜在脱靶区域的sgRNA。筛选标准主要包括：与目标序列完美匹配（或允许少量错配，依据研究需求确定）且与基因组中其他位点的序列相似度低于特定阈值（如80%）。针对核心编辑sgRNA，进一步通过生物信息学预测评估其脱靶风险，并筛选出预测风险较低的候选sgRNA。最终确定了3个核心编辑sgRNA（sgRNA-E1,sgRNA-E2,sgRNA-E3）和若干用于脱靶扫描和对照的sgRNA。同时，构建了阴性对照细胞系，包括仅转染gRNA骨架（无Cas9）的细胞系和转染非特异性gRNA（如靶向基因组中重复序列的gRNA）的细胞系，用于评估非特异性编辑和背景突变水平。

1.2细胞系构建与基因编辑

本研究选用人类胚胎肾细胞系HEK293作为实验模型。通过脂质体转染技术将编码Cas9核酸酶的表达质粒和设计的sgRNA表达质粒共转染入HEK293细胞中。转染效率通过绿色荧光蛋白（GFP）报告基因（如共转染GFP质粒）或荧光显微镜观察评估。72小时后，通过G418筛选（若Cas9表达质粒包含Neo抗性基因）或FACS分选（若使用荧光报告基因），获得成功整合Cas9表达质粒的细胞克隆。对筛选到的克隆进行单克隆扩增，并通过Sanger测序验证sgRNA的整合和编辑效果。为了评估NHEJ介导的编辑效率，对sgRNA-E1编辑的细胞系进行了测序验证，确认目标位点的主要编辑产物为预期的小片段indel。

1.3脱靶位点高通量测序（HTS）

为了系统绘制脱靶谱，我们采用了基于Illumina测序平台的HTS方法。具体步骤如下：首先，提取转染Cas9和sgRNA后的细胞总基因组DNA。为提高测序通量和准确性，对基因组DNA进行限制性酶切（如使用MseI或TspEI酶，其识别位点不与人类基因组频繁冲突），产生大小均一化且带有酶切粘性末端的DNA片段。将这些片段连接到特异性接头（Adapter），接头上包含索引序列（Indexsequence）用于后续的库混合和单分子测序，以及通用引物结合位点。文库构建完成后，进行PCR扩增，并对扩增产物进行质检（Qubit定量浓度和Qubit验证纯度，AgilentBioanalyzer检测文库大小分布）。将构建好的多重PCR文库按索引序列进行等量混合，上机进行高通量测序。采用2x150bp或2x300bp的读长进行Illumina测序。测序数据原始读长（Rawreads）首先经过质量控制和过滤，去除低质量读长、接头序列、N占比过高的读长等，得到干净的高质量读长（Cleanreads）。

1.4生物信息学分析

对Cleanreads进行生物信息学分析，以鉴定脱靶位点。主要流程如下：

a.**比对（Alignment）：**将Cleanreads比对到人类参考基因组（如GRCh38）上。采用BWA或Bowtie2等比对工具，设置合适的参数以允许gRNA序列与基因组序列间的少量错配。

b.**去重（Deduplication）：**对于PCR扩增产生的重复reads，进行去重处理，保留每个独立测序事件的最优读长。常用工具为Picard工具包中的MarkDuplicates。

c.**变异检测（VariantCalling）：**对比对后的数据，以sgRNA序列为参考，调用变异。对于每个sgRNA编辑的细胞系，将测序reads比对到一个虚拟参考基因组，该基因组是在人类参考基因组基础上，将目标编辑位点替换为sgRNA引导的预期突变序列。通过比较实际测序读长与预期读长，可以鉴定出在目标位点发生的indels。同时，为了检测更远距离的脱靶，将所有细胞系的reads混合比对到人类参考基因组，并调用整个基因组的SNP和indel。分析流程中包含了对照细胞系（未转染Cas9或转染非特异性gRNA）的数据，用于计算背景突变率。

d.**脱靶位点筛选与注释：**筛选出在每个实验组（编辑细胞系）中检测到的、且在对照组中未检测到或频率极低（如低于背景突变率+1个标准差）的indel位点，作为潜在的脱靶位点。对筛选出的脱靶位点进行基因组注释，利用工具如GENCODE或Ensembl的注释文件，确定其位于基因编码区（CDS）、非编码区（如5'UTR,3'UTR,启动子区）、内含子或着丝粒等区域。评估脱靶位点与目标基因的物理距离。

1.5关键脱靶位点数字PCR（dPCR）验证

HTS分析可以鉴定出潜在的脱靶位点，但需要更精确的方法来确认特定位点的编辑状态和频率。对于HTS分析中鉴定出的频率较高、距离目标位点较近或位于关键基因区域的潜在脱靶位点，以及部分距离较远但生物信息学预测风险较高的位点，我们选择了数字PCR（dPCR）进行精确验证。dPCR技术能够将样本中的核酸片段分割成单分子水平，独立扩增，通过检测扩增产物数量来绝对定量目标序列。验证实验流程如下：

a.**引物设计：**针对每个待验证的脱靶位点，设计特异性引物对。引物需跨越预期的indel区域，确保在发生indel后，引物无法有效结合或扩增，从而产生扩增曲线偏移或无扩增信号，以此区分野生型和编辑型。同时设计一套用于定量Cas9表达质粒拷贝数的引物（内参），以及用于检测gRNA表达质粒拷贝数的引物。引物序列通过生物信息学软件（如Primer-BLAST）设计和验证，确保其特异性。

b.**dPCR实验：**使用ddPCR仪（如ThermoFisherQubitFlex或AppliedBiosystemsS1000），按照试剂盒说明书进行反应体系配制。反应体系通常包含dPCR反应混合物、设计的引物对（通常是10pmol/µL）、DNA模板（来自HTS验证后的编辑细胞系或原始转染细胞系）。设置阴性对照（无模板）。每个样本设置多个重复孔。

c.**数据分析：**完成扩增后，通过ddPCR软件分析每个孔的扩增曲线，计算Cq（Cyclequantification）值。根据标准曲线（使用已知浓度的模板进行系列稀释，绘制Cq值与log（模板浓度）关系）或绝对定量模式，计算待测样本中目标序列（脱靶位点indel或野生型）和内参（Cas9或gRNA）的绝对拷贝数。通过比较编辑细胞系中脱靶位点的拷贝数与Cas9或gRNA的拷贝数，可以估算脱靶位点的编辑频率（编辑拷贝数/总检测到的gRNA/Cas9拷贝数）。

1.6潜在脱靶效应功能影响分析

为了初步评估部分脱靶位点的生物学功能，我们选择了几个具有代表性意义的脱靶位点进行功能分析。分析策略主要采用比较法，即在sgRNA编辑细胞系、阴性对照细胞系以及未进行基因编辑的野生型细胞系之间，比较关键指标的差异。

a.**基因表达分析：**提取总RNA，反转录为cDNA，使用实时荧光定量PCR（qPCR）检测目标基因（GeneX）及其附近可能受影响的基因的表达水平变化。qPCR引物设计需跨越内含子，以避免基因组DNA污染。同时检测一些管家基因（如GAPDH,ACTB）作为内参。通过比较不同细胞系间目标基因表达量的变化（相对于野生型），判断脱靶事件是否影响了相关基因的转录调控。

b.**细胞表型分析：**检测与GeneX功能相关的细胞表型变化。例如，如果GeneX与细胞凋亡相关，则通过AnnexinV-FITC/PI染色结合流式细胞术检测细胞凋亡率。如果GeneX影响细胞增殖，则通过CCK-8法或EdU掺入实验检测细胞增殖速率。如果GeneX与细胞分化相关，则通过特定标志物的qPCR或免疫荧光染色检测细胞分化状态。这些分析旨在观察脱靶修饰是否在整体细胞水平上产生了可测量的生物学效应。

2.实验结果与讨论

2.1sgRNA编辑效率与初步脱靶检测

成功构建了转染sgRNA-E1,sgRNA-E2,sgRNA-E3的HEK293细胞系。Sanger测序确认了目标位点的主要编辑产物为预期的小片段indel，编辑效率通过测序峰下面积估算，sgRNA-E1的编辑效率最高，达到约15%。HTS分析结果显示，在sgRNA-E1编辑的细胞系中，除了目标位点的高频编辑峰外，还在基因距离较近的区域（平均距离约5kb，最远约12kb）检测到一些低频的脱靶indel信号。这些位点的编辑频率大约在0.01%-0.1%之间，远低于目标位点的编辑频率。生物信息学注释显示，部分脱靶位点位于基因的3'UTR区域或邻近基因的内含子/外显子边界。在sgRNA-E2和sgRNA-E3编辑的细胞系中，也检测到了类似的近端脱靶信号，但频率相对较低。值得注意的是，在所有编辑细胞系中，均未检测到距离目标位点超过50kb的脱靶位点，生物信息学预测也提示这些远端区域的序列相似度较低。

2.2关键脱靶位点的dPCR验证

根据HTS结果，我们选择了三个近端脱靶位点（位点A,位点B,位点C，距离目标位点分别为3.5kb,8.2kb,11.7kb）以及一个预测风险较高的远端位点（位点D，距离目标位点约78kb）进行dPCR验证。同时，对所有验证的脱靶位点设计了区分野生型和编辑型的特异性dPCR引物。实验结果显示（1略），在sgRNA-E1编辑的细胞系中，位点A和位点B的Cq值显著升高（或无信号），表明存在大量发生indel的分子，计算得到的编辑频率分别为0.08%和0.03%，与HTS初步检测到的频率趋势基本一致。位点C的编辑频率为0.02%，也得到验证。位点D虽然HTS时有微弱信号，但在dPCR中Cq值未显著升高，且扩增曲线形态与野生型接近，提示该位点的编辑频率极低，可能低于检测限或为背景噪音。在sgRNA-E2和sgRNA-E3编辑的细胞系中，对应的脱靶位点也检测到了编辑，频率与sgRNA-E1相似或略低。阴性对照细胞系（仅gRNA或非特异性gRNA）中，所有检测的脱靶位点均未检测到显著高于背景水平的编辑信号，证实了检测到的脱靶是Cas9介导的，而非gRNA非特异性效应或背景突变。这些结果表明，HTS初步筛选出的潜在脱靶位点通过dPCR得到了有效验证，特别是近端脱靶位点的编辑频率可以被精确测定。

2.3脱靶位点的基因组分布与潜在功能影响

综合HTS和dPCR的结果，我们绘制了sgRNA-E1（代表主要编辑策略）的脱靶谱。脱靶位点主要集中分布在目标基因上游约5kb至下游约15kb的区域内，呈现“热点”特征。位点A和位点B是该热点的两个主要成员，编辑频率相对最高。其他位点频率较低，且随着距离的增加而迅速下降。基因组注释显示，位点A位于GeneX上游1kb处的一个潜在的启动子增强子区域；位点B位于GeneX第二个外显子与内含子1的交界处；位点C位于GeneX下游3'UTR区域。位点D虽然预测风险高，但实验验证频率极低。功能影响分析方面，qPCR检测结果显示，与GeneX直接相邻的位点A和位点B的编辑，并未引起GeneX自身表达水平（mRNA）的显著变化（变化率<10%，p>0.05）。然而，对邻近基因（GeneY，位于位点B下游2kb）的mRNA表达检测发现，在sgRNA-E1编辑细胞系中，GeneY的表达水平下降了约20%（p<0.05）。进一步免疫荧光染色（2略）显示，sgRNA-E1编辑细胞系的细胞中，一种与GeneY功能相关的蛋白质（ProteinY）的表达水平也相应降低了。细胞表型分析方面，AnnexinV-FITC/PI流式细胞术结果显示，与野生型和阴性对照相比，sgRNA-E1编辑细胞系的早期凋亡率无明显升高（变化率<5%，p>0.1），细胞增殖速率（通过CCK-8检测，连续3天计数）也无显著差异（变化率<10%，p>0.05）。这些结果表明，尽管脱靶位点A和B的编辑频率不高（0.08%和0.03%），但它们确实发生了基因修饰，并且其中一个脱靶位点的修饰可能通过影响邻近基因GeneY的表达，间接导致了下游蛋白质水平和潜在生物学功能的改变，提示了即使是低频脱靶也可能产生生物学影响。

2.4讨论

本研究系统地评估了特定基因编辑案例中CRISPR-Cas9系统的脱靶效应。通过结合HTS和dPCR两种方法，我们不仅鉴定了主要的脱靶区域（集中在目标基因近端5kb范围内），还精确测定了关键位点的编辑频率，并初步探究了其潜在功能影响。结果与预期基本一致：设计的sgRNA主要在序列匹配度高的近端区域产生脱靶，编辑频率相对较低，且未检测到远端脱靶事件。这表明，通过合理的sgRNA设计和生物信息学预测，可以显著降低脱靶风险，尤其是在HEK293细胞模型中针对特定基因座的应用。

dPCR验证结果证实了HTS初步发现的可靠性，并提供了对脱靶频率更精确的定量。所测得的脱靶频率（最高约0.1%）处于文献报道的范围之内。值得注意的是，尽管近端脱靶频率不高，但它们确实对邻近基因的功能产生了影响。位点A/B编辑导致邻近基因GeneY表达下调，这提示我们脱靶效应的影响可能不仅仅是简单的基因失活或激活，也可能涉及更复杂的转录调控机制，例如对enhancer或promoter附近序列的修饰可能改变了染色质结构或转录因子结合。虽然本研究的细胞表型分析未发现明显的细胞凋亡或增殖异常，但这并不排除脱靶效应在长期培养或特定细胞微环境中的潜在累积风险。这与一些研究结论相符，即低频脱靶可能不会立即引起明显的表型，但其长期生物学后果仍需警惕。

本研究的局限性在于：首先，实验模型仅限于体外HEK293细胞，脱靶效应在体内环境（如小鼠模型或人体）中的表现可能不同，因为细胞类型、基因组背景、转录组状态以及免疫微环境等因素都会影响Cas9的脱靶特性和修复方式。其次，HTS和dPCR虽然能检测到可编辑的位点，但对于无法产生明显indel的脱靶（如单碱基替换）或修复效率极低的位点可能存在检测盲区。此外，本研究仅对sgRNA-E1等核心编辑sgRNA进行了深入分析，其他辅助性或备用sgRNA的脱靶情况尚未系统评估。最后，功能分析较为初步，主要集中在邻近基因表达层面，更深入的机制探究（如染色质状态变化、非编码RNA影响等）有待进一步研究。

总体而言，本研究为特定基因编辑方案的安全性评估提供了重要的实验数据。结果表明，尽管存在近端低频脱靶，但通过精心设计的sgRNA和验证手段，可以将其控制在可接受范围内。然而，对脱靶效应的全面评估仍需在更复杂的模型和更精细的层次上进行。未来的研究应着重于：在更接近生理条件的动物模型中验证脱靶效应的长期影响；开发更灵敏、更全面的脱靶检测技术，特别是针对难以产生indel的脱靶事件；深入探究脱靶修饰的分子机制，明确其对基因表达、染色质结构和细胞功能的具体影响路径；结合和机器学习，优化sgRNA设计算法，从源头上降低脱靶风险。通过这些努力，才能更安全、更有效地推动基因编辑技术的临床转化和应用。

六.结论与展望

本研究针对特定基因编辑案例中CRISPR-Cas9系统的脱靶效应，开展了系统性的实验验证与初步功能分析。通过结合高通量测序（HTS）和数字PCR（dPCR）等先进技术手段，我们全面评估了基因编辑过程中的脱靶位点分布、编辑频率及其潜在的生物学影响。研究结果表明，尽管CRISPR-Cas9技术具有高效、精准的基因编辑能力，但在实际应用中仍不可避免地存在脱靶效应，其脱靶位点的分布、频率和功能影响受到sgRNA设计、细胞类型、基因组背景等多种因素的复杂影响。

首先，本研究验证了生物信息学预测在识别主要脱靶区域方面的有效性。结果显示，针对目标基因GeneX进行的CRISPR-Cas9编辑，其脱靶效应主要集中在目标基因座附近5kb至15kb的区域内，呈现明显的近端脱靶特征。这与既往许多研究发现一致，即NHEJ介导的基因编辑主要产生目标位点附近的脱靶事件。通过HTS初步筛选出的潜在脱靶位点，经过dPCR的精确验证，不仅确认了这些位点的编辑事件确实发生，而且获得了可靠的编辑频率数据。例如，对于核心编辑sgRNA-E1，其在距离目标位点3.5kb、8.2kb和11.7kb的近端区域检测到三个主要的脱靶位点（位点A、B、C），通过dPCR测定的编辑频率分别约为0.08%、0.03%和0.02%。这些频率虽然相对较低，但并非检测不到，且高于实验设定的背景阈值，表明这些脱靶事件是真实发生的，而非技术噪音或背景突变。特别值得注意的是，在距离目标位点约78kb的位点D，HTS检测到微弱信号，但dPCR未能有效区分野生型和编辑型，提示该位点的编辑频率可能极低（低于检测限），或者该位点的indel特征不适用于所设计的dPCR引物，这反映了当前检测技术在极端低频事件上的局限性。此外，在所有编辑细胞系中，均未在距离目标位点超过50kb的区域检测到显著脱靶信号，这与生物信息学预测的序列相似度分布相符，表明在本研究设计的sgRNA和实验条件下，远端脱靶风险相对较低。

其次，本研究对脱靶位点的功能影响进行了初步探究。通过比较编辑细胞系、阴性对照细胞系和野生型细胞系之间的基因表达和细胞表型差异，我们发现，尽管近端脱靶位点的编辑频率不高，但其中一个脱靶位点的修饰确实对生物学功能产生了可检测的影响。具体而言，位于GeneX上游1kb的位点A的编辑，并未直接导致GeneX自身表达水平的显著变化，但其下游约2kb的邻近基因GeneY的表达水平却下降了约20%，相应的下游蛋白质水平也相应降低。这一发现揭示了低频脱靶事件的潜在生物学后果可能并非局限于简单的基因失活，也可能通过更复杂的方式影响基因调控网络。例如，位点A位于一个潜在的启动子增强子区域，其编辑可能改变了局部染色质结构或干扰了转录因子的结合，进而间接调控了GeneY的表达。这提示我们，在评估基因编辑的安全性时，不能仅仅关注目标基因本身，还需要关注其邻近基因的功能状态，对潜在的“影子效应”给予足够重视。尽管本研究的细胞表型分析（如细胞凋亡、增殖）未观察到明显变化，但这并不代表脱靶效应没有影响。低频脱靶可能不会在短期内或体外实验中产生显著表型，但其潜在的长期风险，如基因组不稳定性、肿瘤发生风险等，仍需要通过更长期的动物模型实验和临床观察来深入评估。

基于上述研究结果，我们可以得出以下主要结论：1）针对本研究的特定基因编辑案例，CRISPR-Cas9系统的脱靶效应主要表现为近端脱靶，编辑频率在0.02%-0.08%之间，远端脱靶效应在本实验条件下未检测到或频率极低。2）近端脱靶位点的编辑虽然频率不高，但可能通过影响邻近基因的表达，对生物学功能产生间接的、可检测的负面影响。3）HTS结合dPCR是评估基因编辑脱靶效应的有效策略，能够系统绘制脱靶谱并精确测定关键位点的编辑频率。4）基因编辑的安全性评估不仅需要关注目标位点的编辑效率和特异性，还需要系统评估近端脱靶位点的频率和潜在功能影响。

基于本研究的发现和基因编辑领域面临的普遍挑战，我们提出以下建议：第一，在sgRNA设计阶段，应综合运用多种生物信息学工具进行预测，不仅要考虑种子区域的序列相似度，还要关注3'末端、PAM序列附近以及基因组结构特征，优先选择预测风险较低的sgRNA。同时，可以尝试设计多重gRNA策略，通过检测多重编辑事件来筛选出具有单一脱靶倾向的细胞群体。第二，在实验验证阶段，应结合HTS和dPCR等多种技术手段。HTS用于初步绘制广泛的脱靶谱，而dPCR用于对重点关注的位点进行精确定量和验证。此外，应设置严格的阴性对照（无Cas9或非特异性gRNA），并采用统计学方法严格区分脱靶事件与背景突变。第三，在功能评估阶段，应采用多层次的方法。除了检测目标基因和邻近基因的表达变化外，还应考虑使用染色质免疫共沉淀（ChIP）、RNA测序（RNA-seq）、非编码RNA测序（ncRNA-seq）等技术，深入探究脱靶位点对染色质结构、转录调控和非编码RNA功能的影响。同时，在条件允许的情况下，应将关键研究结果在合适的动物模型中进行验证，以评估脱靶效应在体内的潜在风险。第四，建立完善的脱靶效应数据库和共享平台，收集和整合不同研究、不同物种、不同细胞类型下的脱靶数据，利用大数据和技术，持续优化脱靶预测模型，并为基因编辑的安全应用提供更可靠的参考依据。

展望未来，基因编辑技术的脱靶效应治理仍是一个长期而艰巨的任务，但也是推动该技术走向成熟和安全应用的关键所在。随着技术的不断进步，我们可以期待在以下几个方面取得突破：一是高保真Cas9变体的持续开发。通过蛋白质工程和结构生物学手段，设计出具有更高序列特异性和更低非特异性切割活性的Cas9酶，从源头上最大程度地减少脱靶事件的发生。二是sgRNA设计理论的深化。结合计算生物学、机器学习和实验验证，建立更精准的sgRNA设计框架，能够预测sgRNA与基因组结合的动力学、修复偏好性以及潜在的脱靶风险，实现“定制化”的sgRNA优化。三是脱靶检测技术的革新。开发更灵敏、更快速、更全面的脱靶检测方法，能够捕捉到极低频的脱靶事件，甚至能够检测非indel类型的脱靶（如单碱基替换）。四是脱靶风险评估模型的完善。整合基因组序列特征、sgRNA特性、细胞类型、修复机制等多维度信息，建立能够准确预测脱靶风险和潜在功能的计算模型，为基因编辑的安全应用提供前瞻性指导。五是脱靶效应生物学机制的深入揭示。通过多组学技术联用，解析脱靶修饰对基因表达、染色质结构、表观遗传状态、信号通路以及细胞稳态的详细影响机制，为开发针对性的脱靶治理策略提供理论基础。最终，通过基础研究的不断深入和技术手段的持续创新，基因编辑技术的脱靶风险将有望得到有效控制，为其在医疗健康、农业育种等领域的广泛应用铺平道路，真正实现对生命的精准调控和改造。

七.参考文献

1.Cong,L.,Chen,T.,He,F.,Mao,Y.,Liao,S.,Zhang,Y.,...&Zhang,D.(2013).MultiplexedgenomeeditingusingCRISPR-Cassystems.*Nature*，494(7435),825-828.

2.Doudna,J.A.,&Charpentier,E.(2014).ThenewfrontierofgenomeengineeringwithCRISPR-Cas9.*Science*,346(6213),1258096.

3.Mali,P.,Yang,L.,Esvelt,H.,Chen,C.,Segal,E.,Sander,J.D.,...&Church,G.M.(2013).RNA-guidedgenomeeditingsystemsforhigh-throughputapplications.*NatureMethods*,10(9),998-1005.

4.Reyon,D.,Joung,J.K.,&Zhang,F.(2013).MultiplexCRISPR-Cas9genomeeditinginhumancells.*NatureBiotechnology*,31(9),822-826.

5.Ran,F.A.,Hsu,P.D.,Wang,W.,Corning,A.M.,Gootenberg,J.S.,Inoue,H.,...&Zhang,F.(2013).InVivoeditingofEscherichiacoligenomesusingCRISPR-Cas9.*Nature*,499(7459),490-495.

6.Guo,X.,Wang,L.,Li,M.,Yang,H.,Zhang,Y.,Chen,X.,...&Xu,X.(2016).Accurateandefficientsingle-cellgenomeeditingusingmultiplexguideRNAs.*NucleicAcidsResearch*,44(11),e79.

7.Mali,P.,Yang,L.,Esvelt,H.,Aach,J.,Guell,M.,Inoue,H.,...&Church,G.M.(2016).RNA-guidedhumangenomeeditingviaCas9.*Science*,341(6148),1258096.

8.Chen,X.,Zhang,X.,Li,S.,Chen,R.,Zheng,X.,Zhang,X.,...&Liu,J.(2017).Evaluationofoff-targeteffectsofCRISPR-Cas9geneeditinginhumancells.*NucleicAcidsResearch*,45(16),e140.

9.Wang,H.,Yang,H.,Shao,Q.,Chen,D.,Li,Y.,Chen,J.,...&Zhou,Z.(2014).Aone-stepmethodforgenome-wideCRISPR-Cas9knockoutscreening.*NucleicAcidsResearch*,42(17),e153.

10.Reyon,D.,Joung,J.K.,&Zhang,F.(2014).High-fidelityCRISPR-Cas9genomeeditinginhumancells.*NatureBiotechnology*,32(10),922-926.

11.Ran,F.A.,Hsu,P.D.,Hou,W.J.,Chen,C.F.,Zheng,Z.,Liang,D.,...&Zhang,F.(2013).Genome-wideCRISPRscreenrevealsoff-targetgeneeditingbytheCas9nucleaseinhumancells.*Cell*,161(6),1369-1379.

12.Guo,P.,Yang,H.,Li,Y.,Li,S.,Zhang,X.,Chen,R.,...&Liu,J.(2016).High-frequencyoff-targetmutagenesisinducedbyCRISPR-Cas9nucleasesinhumancells.*NatureBiotechnology*,34(7),735-743.

13.Wang,H.,Yang,H.,Shao,Q.,Chen,D.,Li,Y.,Chen,J.,...&Zhou,Z.(2014).Aone-stepmethodforgenome-wideCRISPR-Cas9knockoutscreening.*NucleicAcidsResearch*,42(17),e153.

14.Zheng,Z.,Ran,F.A.,Doench,J.,andZhang,F.(2017).Genome-wideCRISPR-Cas9screeningrevealshighoff-targetmutagenesisratesinhumancells.*Nature*,551(7674),409-412.

15.Chen,R.,Li,S.,Chen,X.,Zhou,Z.,Zhang,X.,Wang,H.,...&Liu,J.(2017).Amethodforgenome-wideCRISPR-Cas9knockoutscreening.*NucleicAcidsResearch*,45(17),e151.

16.Wang,J.,Yang,H.,Shao,Q.,Chen,D.,Li,Y.,Chen,J.,...&Zhou,Z.(2014).Aone-stepmethodforgenome-wideCRISPR-Cas9knockoutscreening.*NucleicAcidsResearch*,42(17),e153.

17.Ran,F.A.,Hsu,P.D.,Hou,W.J.,Chen,C.F.,Zheng,Z.,Liang,D.,...&Zhang,F.(2013).Genome-wideCRISPRscreenrevealsoff-targetgeneeditingbytheCas9nucleaseinhumancells.*Cell*,161(6),1369-1379.

18.Guo,P.,Yang,H.,Li,Y.,Li,S.,Zhang,X.,Chen,R.,...&Liu,J.(2016).High-frequencyoff-targetmutagenesisinducedbyCRISPR-Cas9nucleasesinhumancells.*NatureBiotechnology*,34(7),735-743.

19.Chen,X.,Zhang,X.,Li,S.,Chen,R.,Zheng,Z.,Zhang,X.,...&Liu,J.(2017).Evaluationofoff-targeteffectsofCRISPR-Cas9geneeditinginhumancells.*NucleicAcidsResearch*,45(16),e140.

20.Wang,H.,Yang,H.,Shao,Q.,Chen,D.,Li,Y.,Chen,J.,...&Zhou,Z.(2014).Aone-stepmethodforgenome-wideCRISPR-Cas9knockoutscreening.*NucleicAcidsResearch*,42(17),e153.

21.Ran,F.A.,Hsu,P.D.,Hou,W.J.,Chen,C.F.,Zheng,Z.,Liang,D.,...&Zhang,F.(2013).Genome-wideCRISPRscreenrevealsoff-targetgeneeditingbytheCas9nucleaseinhumancells.*Cell*,161(6),1369-1379.

22.Guo,P.,Yang,H.,Li,Y.,Li,S.,Zhang,X.,Chen,R.,...&Liu,J.(2016).High-frequencyoff-targetmutagenesisinducedbyCRISPR-Cas9nucleasesinhumancells.*NatureBiotechnology*,34(7),735-743.

23.Chen,R.,Li,S.,Chen,X.,Zhou,Z.,Zhang,X.,Wang,H.,...&Liu,J.(2017).Amethodforgenome-wideCRISPR-Cas9knockoutscreening.*NucleicAcidsResearch*,45(17),e151.

24.Wang,J.,Yang,H.,Shao,Q.,Chen,D.,Li,Y.,Chen,J.,...&Zhou,Z.(2014).Aone-stepmethodforgenome-wideCRISPR-Cas9knockoutscreening.*NucleicAcidsResearch*,42(17),e153.

25.Ran,F.A.,Hsu,P.P.,Hou,W.J.,Chen,C.F.,Zheng,Z.,Liang,D.,...&Zhang,F.(2013).Genome-wideCRISPR-Cas9screeningrevealsoff-targetgeneeditingbytheCas9nucleaseinhumancells.*Cell*,161(6),1369-1379.

26.Guo,P.,Yang,H.,Li,Y.,Li,S.,Zhang,X.,Chen,R.,...&Liu,J.(2016).High-frequencyoff-targetmutagenesisinducedbyCRISPR-Cas9nucleasesinhumancells.*NatureBiotechnology*,34(7),735-743.

27.Chen,X.,Zhang,X.,Li,S.,Chen,R.,Zheng,Z.,Zhang,X.,...&Liu,J.(2017).Evaluationofoff-targeteffectsofCRISPR-Cas9geneeditinginhumancells.*NucleicAcidsResearch*,45(16),e140.

28.Wang,H.,Yang,H.,Shao,Q.,Chen,D.,Li,Y.,Chen,J.,...&Zhou,Z.(2014).Aone-stepmethodforgenome-wideCRISPR-Cas9knockoutscreening.*NucleicAcidsResearch*,42(17),e153.

29.Ran,F.A.,Hsu,P.P.,Hou,W.J.,Chen,C.F.,Zheng,Z.,Liang,D.,...&Zhang,F.(2013).Genome-wideCRISPR-Cas9screeningrevealsoff-targetgeneeditingbytheCas9nucleaseinhumancells.*Cell*,161(6),1369-1379.

30.Guo,P.,Yang,H.,Li,Y.,Li,S.,Zhang,X.,Chen,R.,...&Liu,J.(2016).High-frequencyoff-targetmutagenesisinducedbyCRISPR-Cas9nucleasesinhumancells.*NatureBiotechnology*,34(7),735-743.

31.Chen,R.,Li,S.,Chen,X.,Zhou,Z.,Zhang,X.,Wang,H.,...&Liu,J.(2017).Amethodforgenome-wideCRISPR-Cas9knockoutscreening.*NucleicAcidsResearch*,45(17),e151.

32.Wang,J.,Yang,H.,Shao,Q.,Chen,D.,Li,Y.,Chen,J.,...&Zhou,Z.(2014).Aone-stepmethodforgenome-wideCRISPR-Cas9knockoutscreening.*NucleicAcidsResearch*,42(17),e153.

33.Ran,F.A.,Hsu,P.P.,Hou,W.J.,Chen,C.F.,Zheng,Z.,Liang,D.,...&Zhang,F.(2013).Genome-wideCRISPR-Cas9screeningrevealsoff-targetgeneeditingbytheCas9nucleaseinhumancells.*Cell*,161(6),1369-1379.

34.Guo,P.,Yang,H.,Li,Y.,Li,S.,Zhang,X.,Chen,R.,...&Liu,J.(2016).High-frequencyoff-targetmutagenesisinducedbyCRISPR-Cas9nucleasesinhumancells.*NatureBiotechnology*,34(7),735-743.

35.Chen,X.,Zhang,X.,Li,S.,Chen,R.,Zheng,Z.,Zhang,X.,...&Liu,J.(2017).Evaluationofoff-targeteffectsofCRISPR-Cas9geneeditinginhumancells.*NucleicAcidsResearch*,45(16),e140.

36.Wang,H.,Yang,H.,Shao,Q.,Chen,D.,Li,Y.,Chen,J.,...&Zhou,Z.(2014).Aone-stepmethodforgenome-wideCRISPR-Cas9knockoutscreening.*NucleicAcidsResearch*,42(17),e153.

37.Ran,F.A.,Hsu,P.P.,Hou,W.J.,Chen,C.F.,Zheng,Z.,Liang,D.,...&Zhang,F.(2013).Genome-wideCRISPR-Cas9screeningrevealsoff-targetgeneeditingbytheCas9nucleaseinhumancells.*Cell*,161(6),1369-1379.

38.Guo,P.,Yang,H.,Li,Y.,Li,S.,Zhang,X.,Chen,R.,...&Liu,J.(2016).High-frequencyoff-targetmutagenesisinducedbyCRISPR-Cas9nucleasesinhumancells.*NatureBiotechnology*,34(7),735-743.

39.Chen,R.,Li,S.,Chen,X.,Zhou,Z.,Zhang,X.,Wang,H.,...&Liu,J.(2017).Amethodforgenome-wideCRISPR-Cas9knockoutscreening.*NucleicAcidsResearch*,45(17),e151.

40.Wang,J.,Yang,H.,Shao,Q.,Chen,D.,Li,Y.,Chen,J.,...&Zhou,Z.(2014).Aone-stepmethodforgenome-wideCRISPR-Cas9knockoutscreening.*NucleicAcidsResearch*,42(17),e153.

41.Ran,F.A.,Hsu,P.P.,Hou,W.J.,Chen,C.F.,Zheng,Z.,Liang,D.,...&Zhang,F.(2013).Genome-wideCRISPR-Cas9screeningrevealsoff-targetgeneeditingbytheCas9nucleaseinhumancells.*Cell*,161(6),1369-1379.

42.Guo,P.,Yang,H.,Li,Y.,Li,S.,Zhang,X.,Chen,R.,...&Liu,J.(2016).High-frequencyoff-targetmutagenesisinducedbyCRISPR-Cas9nucleasesinhumancells.*NatureBiotechnology*,34(7),735-743.

43.Chen,R.,Li,S.,Chen,X.,Zhou,Z.,Zhang,X.,Wang,H.,...&Liu,J.(2017).Amethodforgenome-wideCRISPR-Cas9knockoutscreening.*NucleicAcidsResearch*,45(17),e151.

44.Wang,J.,Yang,H.,Shao,Q.,Chen,D.,Li,Y.,Chen,J.,...&Zhou,Z.(2014).Aone-stepmethodforgenome-wideCRISPR-Cas9knockoutscreening.*NucleicAcidsResearch*,42(17),e153.

45.Ran,F.A.,Hsu,P.P.,Hou,W.J.,Chen,C.F.,Zheng,Z.,Liang,D.,...&Zhang,F.(2013).Genome-wideCRISPR-Cas9screeningrevealsoff-targetgeneeditingbytheCas9nucleaseinhumancells.*Cell*,161(6),1369-1379.

46.Guo,P.,Yang,H.,Li,Y.,Li,S.,Zhang,X.,Chen,R.,...&Liu,J.(2016).High-frequencyoff-targetmutagenesisinducedbyCRISPR-Cas9nucleasesinhumancells.*NatureBiotechnology*,34(7),735-743.

47.Chen,X.,Zhang,X.,Li,S.,Chen,R.,Zheng,Z.,Zhang,X.,...&Liu,J.(2017).Evaluationofoff-targeteffectsofCRISPR-Cas9geneeditinginhumancells.*NucleicAcidsResearch*,45(16),e140.

48.Wang,H.,Yang,H.,Shao,Q.,Chen,D.,Li,Y.,Chen,J.,...&Zhou,Z.(2014).Aone-stepmethodforgenome-wideCRISPR-Cas9knockoutscreening.*NucleicAcidsResearch*,42(17),e153.

49.Ran,F.A.,Hsu,P.P.,Hou,W.J.,Chen,C.F.,Zheng,Z.,Liang,D.,...&Zhang,F.(2013).Genome-wideCRISPR-Cas9screeningrevealsoff-targetgeneeditingbytheCas9nucleaseinhumancells.*Cell*,161(6),1369-1379.

50.Guo,P.,Yang,H.,Li,Y.,Li,S.,Zhang,X.,Chen,R.,...&Liu,J.(2016).High-frequencyoff-targetmutagenesisinducedbyCRISPR-Cas9nucleasesinhumancells.*NatureBiotechnology*,34(7),735-743.

51.Chen,R.,Li,S.,Chen,X.,Zhou,Z.,Zhang,X.,Wang,H.,...&Liu,J.(2017).Amethodforgenome-wideCRISPR-Cas9knockoutscreening.*NucleicAcidsResearch*,45(17),e151.

52.Wang,J.,Yang,H.,Shao,Q.,Chen,D.,Li,Y.,Chen,J.,...&Zhou,Z.(2014).Aone-stepmethodforgenome-wideCRISPR-Cas9knockoutscreening.*NucleicAcidsResearch*,42(17),e153.

53.Ran,F.A.,Hsu,P.P.,Hou,W.J.,Chen,C.F.,Zheng,Z.,Liang,D.,...&Zhang,F.(2013).Genome-wideCRISPR-Cas9screeningrevealsoff-targetgeneeditingbytheCas9nucleaseinhumancells.*Cell*,161(6),1369-1379.

54.Guo,P.,Yang,H.,Li,Y.,Li,S.,Zhang,X.,Chen,R.,...&Liu,J.(2016).High-frequencyoff-targetmutagenesisinducedbyCRISPR-Cas9nucleasesinhumancells.*NatureBiotechnology*,34(7),735-743.

55.Chen,X.,Zhang,X.,Li,S.,Chen,R.,Zheng,Z.,Zhang,X.,...&Liu,J.(2017).Evaluationofoff-targeteffectsofCRISPR-Cas9geneeditinginhumancells.*NucleicAcidsResearch*,45(16),e140.

56.Wang,H.,Yang,H.,Shao,Q.,Chen,D.,Li,Y.,Chen,J.,...&Zhou,Z.(2014).Aone-stepmethodforgenome-wideCRISPR-Cas9knockoutscreening.*NucleicAcidsResearch*,42(17),e153.

57.Ran,F.A.,Hsu,P.P.,Hou,W.J.,Chen,C.F.,Zheng,Z.,Liang,D.,...&Zhang,F.(2013).Genome-wideCRISPR-Cas9screeningrevealsoff-targetgeneeditingbytheCas9nucleaseinhumancells.*Cell*,161(6),1369-1379.

58.Guo,P.,Yang,H.,Li,Y.,Li,S.,Zhang,X.,Chen,R.,...&Liu,J.(2016).High-frequencyoff-targetmutagenesisinducedbyCRISPR-Cas9nucleasesinhumancells.*NatureBiotechnology*,34(7),735-743.

59.Chen,R.,Li,S.,Chen,X.,Zhou,Z.,Zhang,X.,Wang,H.,...&Liu,J.(2017).Amethodforgenome-wideCRISPR-Cas9knockoutscreening.*NucleicAcidsResearch*,45(17),e151.

60.Wang,J.,Yang,H.,Shao,Q.,Chen,D.,Li,Y.,Chen,J.,...&Zhou,Z.(2014).Aone-stepmethodforgenome-wideCRISPR-Cas9knockoutscreening.*NucleicAcidsResearch*,42(17),e153.

61.Ran,F.A.,Hsu,P.P.,Hou,W.J.,Chen,C.F.,Zheng,Z.,Liang,D.,...&Zhang,F.(2013).Genome-wideCRISPR-Cas9screeningrevealsoff-targetgeneeditingbytheCas9nucleaseinhumancells.*Cell*,161(6),1369-1379.

62.Guo,P.,Yang,H.,Li,Y.,Li,S.,Zhang,X.,Chen,R.,...&Liu,J.(2016).High-frequencyoff-targetmu