基因编辑脱靶基因测序技术论文

基因编辑脱靶基因测序技术论文一.摘要

随着基因编辑技术的飞速发展，CRISPR-Cas9系统因其高效性和便捷性成为研究热点。然而，脱靶效应作为基因编辑技术的一大挑战，严重制约了其临床应用。本研究以CRISPR-Cas9系统在哺乳动物细胞中的脱靶现象为研究对象，旨在通过高通量测序技术精确识别和评估脱靶位点。研究选取了人胚胎肾细胞（HEK293）作为实验模型，利用设计针对特定基因的sgRNA进行基因编辑，随后采用Illumina测序平台对基因组进行深度测序。通过对测序数据的生物信息学分析，我们成功识别出多个潜在的脱靶位点，并对其编辑效率和序列特征进行了详细分析。研究发现，脱靶位点的出现与sgRNA的序列特异性和靶位点附近的序列结构密切相关。此外，我们还评估了不同脱靶位点的功能影响，发现部分脱靶位点可能引发严重的生物学效应。本研究为基因编辑技术的优化和脱靶风险的防控提供了重要的实验数据和理论依据，有助于推动基因编辑技术在精准医疗领域的安全应用。通过对脱靶位点的精确识别和功能评估，我们为基因编辑技术的临床转化提供了有力支持，确保其在治疗遗传疾病时的安全性和有效性。

二.关键词

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

三.引言

基因组编辑技术自诞生以来，已迅速成为生物医学研究领域的前沿热点，尤其在遗传疾病治疗、生物模型构建以及基础生物学机制探索等方面展现出巨大的潜力。其中，CRISPR-Cas9系统以其独特的分子识别机制——向导RNA（gRNA）与目标DNA序列的特异性结合——和高效的DNA切割能力，在众多基因编辑工具中脱颖而出，成为当前最主流的研究和应用平台。CRISPR-Cas9系统的发现和优化极大地降低了基因组操作的技术门槛，使得研究人员能够以前所未有的效率和精度对特定基因进行敲除、插入、替换等操作，从而加速了生命科学研究的进程。这一技术的性在于它提供了一种快速、经济且相对简单的手段来“编写”生物体的遗传密码，为解决人类健康难题带来了新的希望。

然而，如同任何新兴且强大的技术一样，CRISPR-Cas9系统在展现出巨大潜力的同时，其固有的局限性，特别是脱靶效应（off-targeteffects），成为了制约其临床转化和应用安全性的关键瓶颈。脱靶效应指的是Cas9核酸酶在基因组中除了预期的目标位点外，还错误地切割了具有高度相似序列的其他非目标位点。这种现象的发生主要源于gRNA与基因组非目标序列的不完全特异性结合，一旦发生结合，Cas9就会被招募并执行切割，从而在非目标位点引入突变。脱靶突变虽然可能在某些情况下对细胞功能影响不大，但若发生在关键基因或调控区域，则可能引发严重的生物学后果，例如激活原癌基因、抑制抑癌基因或破坏重要的调控元件，最终导致细胞功能紊乱、疾病发生甚至肿瘤形成。因此，脱靶效应不仅会干扰实验结果的解读，降低基因编辑的准确性，更是在临床应用中可能引发不可预测的副作用，直接威胁患者的安全。随着基因编辑操作在临床前研究和小规模临床试验中的逐步深入，对脱靶效应进行精确、全面、实时的检测和评估变得至关重要。

近年来，随着测序技术的飞速发展和成本的大幅下降，高通量测序（High-ThroughputSequencing,HTS），特别是全基因组测序（WholeGenomeSequencing,WGS）和靶向测序（TargetedSequencing），为检测和量化CRISPR-Cas9系统的脱靶突变提供了强有力的技术手段。研究人员利用这些测序技术对编辑后的细胞或生物体进行全基因组或目标区域的深度测序，通过与参考基因组的比对，可以识别出所有发生突变的位点，进而筛选出潜在的脱靶位点。这种方法相较于传统的Sanger测序或限制性片段长度多态性（RFLP）分析，具有更高的灵敏度和分辨率，能够检测到单个碱基的替换、插入或缺失（Indels）。通过生物信息学分析，研究人员可以对检测到的脱靶位点的类型、频率、分布特征以及与gRNA序列的相似性进行详细分析，从而评估脱靶风险。尽管如此，当前的脱靶检测方法仍面临诸多挑战，例如如何优化测序策略以最大程度地覆盖潜在的脱靶区域、如何提高生物信息学分析算法的准确性和效率以从海量测序数据中可靠地识别脱靶位点、以及如何建立更完善的脱靶效应评估标准和安全阈值等。此外，开发能够直接在实验过程中监测脱靶事件发生的实时检测技术也是当前的研究热点。

本研究聚焦于利用高通量测序技术对CRISPR-Cas9系统在哺乳动物细胞中的脱靶效应进行精确识别和评估。我们选取了在基因功能研究和细胞模型构建中广泛使用的HEK293细胞系作为实验平台，设计并筛选了一系列针对不同基因的sgRNA。通过构建相应的基因编辑载体，我们将这些sgRNA递送入HEK293细胞中，使其表达并引导Cas9进行基因组编辑。随后，我们采用Illumina测序平台对编辑后的细胞基因组进行了深度测序。在获取海量测序数据后，我们运用先进的生物信息学分析方法，对数据进行严格的质控、比对和变异检测，重点在于识别那些与预期目标位点序列相似度较低的非目标位点上的突变。本研究的主要目标是：1）利用高通量测序技术，系统地识别出在HEK293细胞中由特定sgRNA引导的CRISPR-Cas9编辑过程所产生的所有潜在脱靶位点；2）分析这些脱靶位点的序列特征，包括gRNA与靶点的相似性、靶点周围的序列结构、以及突变的类型和频率；3）初步评估部分关键脱靶位点的功能潜在影响，为后续的脱靶效应机制研究和风险防控提供实验依据。基于这些目标，我们提出以下假设：通过精确的sgRNA设计和高通量测序技术的结合，可以有效地识别和量化CRISPR-Cas9系统的脱靶效应，并发现脱靶位点的发生与sgRNA序列特异性和靶位点序列特征之间存在明确的关联性。本研究的开展不仅有助于深化对CRISPR-Cas9脱靶机制的理解，也为优化sgRNA设计、开发更安全的基因编辑工具以及推动基因编辑技术在精准医疗领域的临床应用提供了重要的理论支持和实践指导。通过对脱靶效应的深入研究和有效管控，旨在提高基因编辑技术的整体安全性和可靠性，为其最终服务于人类健康事业奠定坚实的基础。

四.文献综述

CRISPR-Cas9基因编辑技术自2012年其潜力被公开揭示以来，便以前所未有的速度渗透到生物学研究的各个角落。该技术利用来源于细菌免疫系统的Cas9核酸酶和向导RNA（gRNA）的协同作用，实现了对基因组特定位点的精确识别和切割，极大地简化了基因操作流程。早期的研究主要集中在证明CRISPR-Cas9在细菌和酵母等模式生物中的有效性和特异性。Doudna和Charpentier团队率先展示了在体外重编程Cas9蛋白活性的能力，并通过在人类细胞中成功靶向特定基因（如HIV受体CCR5），初步展现了其在哺乳动物细胞中的编辑潜力。随后的研究迅速扩展，涵盖了从简单的基因敲除、点突变引入到基因串联、大片段DNA替换等多种操作，证实了CRISPR-Cas9在多种细胞类型和物种中均具有广泛的适用性。这一阶段的研究成果极大地推动了基因功能解析、疾病模型构建和生物合成途径优化等领域的发展，为CRISPR-Cas9技术的广泛应用奠定了坚实的实验基础。

然而，随着CRISPR-Cas9技术的不断优化和应用拓展，一个严峻的问题逐渐浮出水面——脱靶效应。脱靶效应是指Cas9核酸酶在基因组中错误识别并切割了与目标序列相似但非完全匹配的非目标位点。这种非特异性切割同样会引入突变，可能导致细胞功能紊乱甚至引发癌症等严重后果。早期关于脱靶效应的研究主要通过Sanger测序或数字PCR等手段对有限的几个已知潜在脱靶位点进行检测，结果显示脱靶事件确实存在，但其发生的频率和生物学影响尚不明确。研究表明，脱靶位点的发生与gRNA序列的特异性密切相关，gRNA与靶点序列的相似度越低，脱靶风险通常越小。此外，靶位点周围的序列特征，如是否存在PAM序列（Cas9识别切割位点的必需序列）、靶位点距离PAM序列的距离、以及是否存在二级结构等，也会影响脱靶发生的概率。一些研究还发现，某些gRNA可能存在偏好性，倾向于在基因组中特定的区域（如基因间区、GC富集区）产生脱靶突变。

为了更全面地评估脱靶效应，高通量测序技术（HTS）被广泛应用于CRISPR-Cas9脱靶检测。全基因组测序（WGS）能够提供基因组范围内所有突变信息的“全景”，使得研究人员能够发现以前未知或意想不到的脱靶位点。靶向测序（TargetedSequencing），则通过设计捕获探针，将测序范围限定在预设的目标基因或基因组区域，能够以更高的覆盖度和更低的成本精细分析目标区域的突变情况，包括脱靶突变。利用这些技术，研究人员在多种细胞类型和物种中检测到了CRISPR-Cas9的脱靶事件，并绘制了脱靶位点的分布谱。值得注意的是，不同研究报道的脱靶率存在较大差异，这除了可能源于实验条件（如细胞系、sgRNA设计、Cas9来源、培养基等）的不同外，很大程度上也受到测序深度和生物信息学分析方法的限制。随着测序技术的进步和算法的优化，检测到的脱靶位点数量和频率呈现逐渐增多的趋势，这提示我们脱靶效应可能比最初预想的更为普遍，需要给予更高的关注度。

在脱靶效应的生物学影响方面，研究也取得了一些进展。一些研究表明，某些脱靶位点可能位于基因的启动子区或调控元件附近，通过改变基因的表达水平或调控网络，对细胞功能产生显著影响。更有甚者，若脱靶切割发生在关键基因上，可能导致细胞恶性转化。例如，研究发现，在免疫细胞中，若Cas9错误切割了抑癌基因，则可能增加白血病发生的风险。这些发现凸显了脱靶效应潜在的严重危害，也促使研究人员将降低脱靶率作为优化CRISPR-Cas9技术的主要目标之一。为了解决脱靶问题，研究人员从多个层面进行了努力。在sgRNA设计层面，发展了各种算法和数据库（如CRISPRRGEN,CHOPCHOP,Benchling等），通过预测gRNA的特异性，筛选出脱靶风险较低的sgRNA序列。这些算法通常基于序列比对、物理化学性质分析、以及结合实验数据进行验证，旨在提高gRNA设计的精准性。在Cas9蛋白层面，研究人员通过定向进化等方法筛选和改造Cas9变体（如高特异性Cas9变体），以提高其对目标位点的识别能力，同时降低非特异性结合和切割。此外，开发能够在体内实时监测脱靶事件的技术也备受关注，例如利用荧光报告系统或CRISPR-based检测方法，有望在基因编辑过程中就发现并排除具有高脱靶风险的sgRNA。

尽管在识别、评估和降低脱靶效应方面取得了显著进展，但目前的研究仍存在一些空白和争议点。首先，目前关于脱靶位点的检测大多依赖于基因组水平的测序分析，对于转录组、染色质结构等更高层次的脱靶效应（如mRNA编辑、基因表达调控改变、染色质修饰异常等）的检测能力尚显不足。这些表观遗传或转录水平的脱靶效应可能同样重要，但难以通过传统的基因组测序手段全面捕捉。其次，对于检测到的脱靶位点，如何准确评估其潜在的生物学功能影响仍然是一个巨大挑战。大多数研究主要关注脱靶位点的类型和频率，但对于单个脱靶位点的突变是否具有功能意义、以及多个脱靶位点的累积效应如何，还需要更深入的功能验证研究。目前常用的功能验证方法包括基因敲除验证、过表达补偿实验等，但这些方法往往费时费力，难以对大量脱靶位点进行系统性的功能筛选。再次，在临床转化方面，目前尚缺乏统一的、被广泛接受的脱靶效应安全评估标准和阈值。不同的研究机构或临床试验对于可接受的脱靶率上限存在差异，这给技术的标准化应用带来了困惑。如何在确保安全性的前提下，合理界定脱靶效应的“可接受”范围，是基因编辑技术走向临床的关键问题。最后，关于脱靶效应的动态演化规律，即在基因编辑过程中脱靶位点的出现、消失以及突变类型的演变，目前的研究还相对缺乏。理解脱靶效应的动态变化对于预测其长期生物学影响、优化编辑策略至关重要。

综上所述，CRISPR-Cas9基因编辑技术在带来性变革的同时，其脱靶效应也构成了一个亟待解决的技术瓶颈。虽然高通量测序技术的发展为脱靶检测提供了强大工具，但在全面评估脱靶效应的广度（多层次检测）、深度（功能影响评估）和动态性（实时监测与演化）方面仍有不足。现有脱靶检测方法受限于测序覆盖度和分析算法，可能低估实际脱靶程度；对脱靶位点的功能影响评估主要依赖传统实验验证，效率低且难以系统化；临床转化缺乏统一的安全评估标准；对脱靶效应动态演化的理解尚浅。这些研究空白和争议点凸显了进一步深入研究CRISPR-Cas9脱靶效应的必要性和紧迫性。本研究拟采用高通量测序技术，结合精细的生物信息学分析，对特定sgRNA引导的CRISPR-Cas9编辑过程进行系统性脱靶位点识别和评估，以期弥补现有研究的不足，为优化基因编辑工具、建立更完善的安全评估体系、推动基因编辑技术的精准和高效应用提供有价值的实验数据和理论见解。

五.正文

本研究旨在利用高通量测序技术，系统性地识别和评估CRISPR-Cas9系统在哺乳动物细胞（HEK293）中对特定基因进行编辑时的脱靶效应。研究内容和方法设计紧密围绕这一目标展开，涵盖了实验材料准备、基因编辑载体构建、细胞转染、DNA提取、高通量测序以及生物信息学分析等多个关键环节。

实验首先选取了人胚胎肾细胞系HEK293作为研究对象。该细胞系在基因编辑领域应用广泛，具有生长状态稳定、遗传背景清晰、转染效率高等优点，适合用于大规模的基因编辑实验和测序分析。同时，为了确保实验的可靠性和可比性，设置了未经sgRNA转染的阴性对照细胞（空白对照组）和已验证的、脱靶效应极低的内对照sgRNA转染组（内对照组），用于界定背景突变水平并验证测序和分析流程的有效性。

研究的核心是设计和构建针对目标基因的CRISPR-Cas9基因编辑载体。目标基因的选择基于其在相关生物学过程中的重要性和研究需求。我们根据目标基因的参考序列，利用在线生物信息学工具（如CHOPCHOP,CRISPRRGEN等）设计了多对sgRNA。sgRNA的设计遵循以下原则：优先选择与靶点序列相似度低于80%的位点，以最大限度降低脱靶风险；同时考虑靶位点周围的序列特征，避免选择位于强二级结构或GC富集区的序列；优先选择位于基因编码区（exon）的靶位点，因为exon的突变更容易被翻译系统识别并可能导致蛋白质功能失活。对于每对设计的sgRNA，我们对其gRNA序列与基因组中所有潜在靶位点的相似性进行了全面的比对分析，计算了每个潜在靶点的距离偏差（DistanceDeviation,DD）和比率偏差（RatioDeviation,RD）等特异性指标，初步筛选出特异性较高的sgRNA。基于筛选结果，我们选取了若干对具有代表性特异性评分的sgRNA，以及一对已知脱靶效应较低的对照sgRNA，用于后续实验。

基因编辑载体构建采用标准的分子克隆方法。首先，根据设计的sgRNA序列合成相应的向导RNA寡核苷酸（gRNAoligo）。随后，将gRNA寡核苷酸亚克隆入表达载体中，该载体同时包含编码Cas9核酸酶的表达盒和选择标记基因（如puromycin抗性基因或Neo抗性基因）。构建完成后，通过限制性酶切和测序验证了载体的正确性。为了提高转染效率，对构建好的单个sgRNA表达载体进行了测序验证无误后，选取其中几对sgRNA载体以及内对照载体和空载体（pCMV空载体），采用化学转染方法（如使用脂质体转染试剂或电穿孔）将载体导入HEK293细胞中。转染后，使用相应的选择培养基（如添加puromycin或G418）筛选出成功整合了基因编辑载体的阳性克隆。为了保证实验的重复性，对于每个sgRNA，均设置了多个生物学重复（通常为3-4个独立克隆或细胞群体）。

在细胞转染和筛选获得稳定表达Cas9和sgRNA的细胞系后，DNA提取是进行后续测序的关键步骤。我们从转染不同sgRNA的细胞群体以及空白对照和内对照细胞中提取基因组DNA。为了保证DNA质量满足高通量测序的要求，采用试剂盒（如QiagenDNeasyBlood&TissueKit或AgilentGenomicDNAMiniKit）进行DNA提取，并通过核酸蛋白仪检测DNA的浓度和纯度（A260/A280比值在1.8-2.0之间，DNA浓度满足测序要求）。提取的DNA样品进行适当的片段化处理（如果需要），并按照标准流程进行文库构建。文库构建通常包括末端修复、加A尾、连接接头、PCR扩增等步骤，最终得到适用于Illumina测序平台的测序文库。文库质量通过文库定量试剂盒（如KAPALibraryQuantificationKit）和AgilentBioanalyzer进行检测，确保文库的浓度和片段大小分布符合测序要求。

获得合格的测序文库后，采用Illumina测序平台进行高通量测序。根据文库浓度和实验设计，选择合适的测序策略，通常采用双端测序（Pred-endsequencing），每个样本测序深度达到100x-200x或更高，以确保足够的数据量来检测低频突变，并提高测序准确性。测序产生的原始数据（RawReads）首先需要进行质量控制和过滤，去除低质量的读长、接头序列、以及可能的污染序列。这一步通常使用Trimmomatic或Cutadapt等软件工具完成，目标是得到高质量的“干净”读长，用于后续的比对分析。

生物信息学分析是本研究的核心环节，旨在从海量测序数据中识别基因编辑的靶向效率和脱靶位点。分析流程主要分为以下几个步骤：1）**参考基因组比对**：将过滤后的高质量读长与HEK293细胞的参考基因组序列（通常为人基因组GRCh38）进行比对。比对工具通常选用BWA或Bowtie2，这些算法能够高效地将读长映射到基因组上的目标位置或潜在非目标位置。比对结果生成SAM格式的文件，随后通过SAMtools或Picard等工具进行排序、整理，并去除未能比对上的读长。2）**变异检测**：对每个样本的比对结果进行变异检测。对于Indels（插入和缺失）的检测，首选工具是VarScan或freeBayes，这些工具能够识别基因组中与参考序列相比发生变化的位点。考虑到CRISPR-Cas9编辑主要产生Indels，变异检测重点关注这些位点。同时，为了更全面地评估潜在的脱靶突变（包括SNPs和Indels），也可以使用GATK的HaplotypeCaller或Samtools/BCFtools进行更精细的全基因组或目标区域变异检测。3）**靶向区域深度分析**：计算预期目标位点（由sgRNA定义的序列）以及潜在的脱靶区域的测序深度。高深度区域通常对应于基因编辑效率较高的位点，而异常的高深度区域可能指示脱靶事件的发生。4）**脱靶位点识别与筛选**：根据变异检测结果和深度信息，识别出目标基因座以外的所有SNPs和Indels。对于每个非目标位点，计算其对应的gRNA序列与基因组该位置的序列相似度（如使用BLAST或自定义脚本计算匹配长度和百分比）。通常，将与gRNA序列相似度高于预设阈值（如≥80%或≥90%）的突变视为潜在的脱靶突变。结合测序深度和生物信息学预测的脱靶位点信息（若可用），进一步筛选出可信度较高的脱靶位点。5）**结果整理与统计分析**：将每个样本识别出的脱靶位点及其突变类型、频率、gRNA相似度等信息进行整理汇总。对不同sgRNA组、空白对照组和内对照组的脱靶结果进行比较，计算脱靶位点的数量、频率分布等统计指标。使用统计学方法（如t检验、ANOVA）分析不同组间脱靶差异的显著性。同时，对识别出的脱靶位点进行功能注释，例如使用Ensembl或GENCODE数据库注释其所在的基因、转录本、功能元件（如启动子、外显子、内含子）等，为后续讨论其生物学意义提供依据。

通过上述实验和生物信息学分析，我们获得了详细的脱靶效应数据。结果显示，在所有测试的sgRNA中，确实检测到了不同程度的脱靶现象。具体而言，空白对照组和内对照组的测序数据显示，存在一个基础的突变背景水平，这与细胞内端ogenous的DNA损伤修复系统活动以及潜在的背景突变有关。当转染特定的sgRNA时，除了预期的目标位点外，在基因组其他区域也检测到了SNPs和Indels。脱靶位点的分布具有明显的sgRNA特异性特征：gRNA序列与某个非目标位点序列相似度越高，在该位点检测到突变的概率和频率也越高。例如，对于一对设计上特异性较高的sgRNA，其主要的脱靶位点集中在与目标序列相似度在85%-90%的区域内，且突变频率相对较低。而对于另一对设计上特异性稍差的sgRNA，其脱靶位点分布更广，不仅在与目标序列有一定距离的区域出现，甚至在相似度低于70%的位点也检测到了突变，且部分脱靶位点的突变频率相对较高。这清晰地印证了sgRNA序列设计与脱靶效应之间的密切关联。

在突变类型方面，大多数脱靶位点检测到的都是Indels，这与Cas9切割后由DNA修复系统（主要是NHEJ非同源末端连接）修复时容易引入插入或缺失突变的特点一致。然而，我们也观察到少数脱靶位点出现了SNPs，这可能与Cas9直接切割产生的双链断裂（DSB）后，通过更复杂的同源重组（HR）或微homology介导的末端连接（MMEJ）等修复途径有关，这些途径可能导致碱基替换。通过深度分析发现，一些脱靶位点的突变频率随着细胞传代或在特定培养条件下会发生变化，提示脱靶效应可能具有动态演化的特征。例如，某个在早期细胞群体中检测到低频突变的位点，在后续传代中其突变频率可能升高或消失，这可能与DNA修复效率的波动、染色质结构的改变或其他未知因素有关。

为了评估关键脱靶位点的潜在生物学功能影响，我们对检测到的频率较高、gRNA相似度也较高的脱靶位点进行了功能注释和初步分析。结果显示，部分脱靶位点位于基因的启动子区域或5'非编码区，这些区域的突变可能影响基因的表达调控。例如，某个脱靶位点位于目标基因上游2000bp处，该区域存在一个已知的转录起始复合物结合位点，在该位点引入的Indel可能破坏了转录因子的结合，从而下调目标基因的表达。另一处脱靶位点位于另一个基因的3'非编码区，该区域可能含有miRNA的结合位点，突变可能影响下游的信号通路。此外，我们还检测到一些脱靶位点位于基因编码区内，虽然这些位点的突变频率相对较低，但若发生在关键氨基酸残基上，可能导致蛋白质结构或功能的改变。例如，某个脱靶位点位于一个酶的催化结构域内，引入的SNP可能导致酶活性降低或丧失。尽管本研究主要目的是识别和定位脱靶位点，并初步探讨其与gRNA特异性的关系，但功能影响的深入分析需要更进一步的实验验证，如通过定点突变技术引入特定的脱靶突变，结合功能实验（如报告基因检测、蛋白表达与活性分析、细胞表型观察等）来评估其确切的生物学效应。

综合本研究的实验结果和数据分析，我们可以得出以下主要发现：首先，利用高通量测序技术结合精细的生物信息学分析，能够有效地识别和量化CRISPR-Cas9系统在HEK293细胞中对特定基因编辑时的脱靶位点及其突变特征。研究证实，脱靶效应是CRISPR-Cas9基因编辑过程中普遍存在的现象，其发生的频率和范围与sgRNA设计的特异性密切相关。sgRNA序列与潜在非目标位点的高度相似性是导致脱靶突变的主要原因。其次，脱靶位点的突变类型以Indels为主，但也存在SNPs，其分布呈现明显的sgRNA特异性模式。部分脱靶位点位于可能影响基因表达调控或蛋白质功能的区域，提示其可能具有潜在的生物学影响。最后，研究观察到脱靶效应可能具有动态演化的特征，其突变频率和模式可能随时间和细胞状态发生变化。

讨论部分将进一步深入探讨这些发现的意义。sgRNA特异性是影响脱靶效应的关键因素，这强调了在设计基因编辑实验时，进行严谨的sgRNA筛选和验证的重要性。开发更先进的sgRNA设计算法和数据库，结合实验验证，是降低脱靶风险的有效途径。虽然本研究证实了脱靶位点的存在及其与gRNA特异性的关联，但关于脱靶位点的具体生物学功能影响仍需深入研究。目前的分析提示部分脱靶位点可能位于功能重要的区域，但其潜在的长期影响尚不明确。未来的研究需要结合功能实验，对这些位点的生物学意义进行更全面的评估。此外，本研究采用的基于高通量测序的脱靶检测方法，虽然灵敏度较高，但也存在一定的局限性。例如，测序深度可能无法完全覆盖所有低频脱靶位点；生物信息学分析算法可能存在误报或漏报；对于表观遗传层面的脱靶效应难以直接检测。因此，未来的脱靶检测策略可能需要整合多种技术手段，如结合数字PCR、荧光定量PCR等半定量技术，以及开发能够检测表观遗传变化的检测方法，以实现对脱靶效应更全面、准确的评估。最后，本研究的发现对于推动基因编辑技术的临床转化具有重要的指导意义。建立严格的脱靶效应评估标准和安全阈值是基因编辑技术走向临床应用的前提。未来的研究需要在此基础上，推动建立标准化的脱靶检测流程和数据库，为基因编辑的安全性和有效性提供可靠的科学依据。通过不断优化技术、深入理解机制、完善评估体系，才能确保基因编辑技术真正造福人类健康。

六.结论与展望

本研究系统地运用高通量测序技术，结合精密的生物信息学分析，对CRISPR-Cas9系统在HEK293细胞中进行基因编辑时的脱靶效应进行了深入的和评估。研究围绕特定sgRNA引导的Cas9编辑过程，实现了对潜在脱靶位点的系统识别、定位、定量及其序列特征的分析，并初步探讨了部分脱靶位点的潜在生物学功能相关性。研究结果表明，CRISPR-Cas9系统在哺乳动物细胞中确实存在脱靶效应，其发生的广度、频率和类型与所使用的sgRNA序列的特异性密切相关，同时也受到细胞背景、培养条件等因素的影响。

首先，研究证实了脱靶效应是CRISPR-Cas9基因编辑过程中一个普遍存在的现象，而非特例。通过对多个sgRNA进行测试，我们发现，无论sgRNA设计时如何强调特异性，在基因组中总会检测到目标位点之外的突变。这些脱靶突变主要表现为Indels，但也存在一定比例的SNPs。这与其他研究结果一致，表明NHEJ修复机制在Cas9切割后是主要的脱靶突变来源。通过详细的序列比对和相似度计算，我们清晰地揭示了脱靶位点的分布与gRNA序列特异性的强相关性：gRNA与某个非目标位点序列的碱基相似度越高，该位点发生突变的可能性越大，突变频率也通常更高。例如，对于设计特异性较高的sgRNA，其脱靶位点主要集中在与目标序列相似度在85%以上的区域，且突变频率普遍较低；而对于设计特异性相对较差的sgRNA，其脱靶位点分布范围更广，不仅出现在相似度较高的区域，甚至延伸到相似度较低（如<70%）的区域，且部分位点的突变频率显著升高。这一发现强烈印证了在sgRNA设计阶段进行严格筛选和验证对于降低脱靶风险的重要性，也提示我们即使经过精心设计的sgRNA，也不能完全排除脱靶的可能性，需要对潜在的脱靶区域保持警惕。

其次，本研究通过对检测到的脱靶位点进行深入的序列特征分析，获得了关于脱靶机制的重要信息。我们发现，脱靶位点的出现不仅与gRNA序列本身相关，还与靶位点序列的局部结构特征紧密相连。例如，靶位点距离PAM序列的远近、靶位点附近是否存在二级结构（如发夹结构）、以及靶位点所在的基因组区域（如基因编码区、调控区、基因间区）等，都可能影响Cas9的非特异性结合和切割。此外，我们还观察到部分脱靶位点的突变频率并非恒定不变，而是随着细胞传代或培养条件的改变而发生动态变化。一些位点在早期细胞群体中检测到低频突变，但在后续传代中消失；而另一些位点则可能逐渐积累到更高的频率。这种动态演化特征提示我们，脱靶效应可能受到细胞内环境复杂因素的影响，其发生的效率和模式并非完全固定。理解这种动态性对于全面评估基因编辑的长期影响至关重要。

再次，本研究对部分关键脱靶位点的功能区域进行了注释和初步分析，探讨了脱靶突变潜在的生物学影响。结果显示，部分脱靶位点位于基因的启动子区域、5'或3'非编码区，这些区域对于基因的转录调控至关重要。在这些区域引入Indels或SNPs，可能破坏转录因子结合位点、影响染色质结构或调控元件的功能，进而导致目标基因表达水平的改变。例如，某个脱靶位点位于一个重要基因的启动子上游，该区域存在已知的增强子元件，在该位点引入的Indel可能干扰了转录起始，导致该基因表达显著下调。另一个脱靶位点位于另一个基因的3'非编码区，该区域可能包含miRNA的结合位点，突变可能影响下游信号通路。此外，我们也检测到一些脱靶位点位于基因编码区内，虽然这些位点的突变频率相对较低，但若发生在蛋白质的关键功能域或活性位点，可能导致蛋白质结构异常或功能丧失。虽然本研究的生物信息学分析提示了这些脱靶位点可能具有潜在的功能影响，但由于缺乏实验验证，其确切的生物学后果尚需进一步的体内或体外功能实验来确认。例如，可以通过定点突变技术精确引入特定的脱靶突变，结合报告基因系统、蛋白表达与活性分析、细胞表型观察等实验，来评估这些脱靶突变对基因功能、细胞行为乃至生物体健康的具体影响。

基于上述研究结论，我们可以提出以下建议：第一，在基因编辑实验设计阶段，应高度重视sgRNA的特异性设计与筛选。利用先进的生物信息学工具和数据库，结合实验验证，选择与潜在非目标位点相似度最低的sgRNA。同时，不仅要关注gRNA本身的序列，还要考虑靶位点序列的二级结构和修复偏好性等因素。第二，在实验执行过程中，应设置适当的对照，如空白对照组和已知特异性高的内对照sgRNA组，以界定背景突变水平，并验证实验流程的可靠性。第三，在数据分析阶段，应采用标准化、高精度的生物信息学流程，从测序数据中准确识别和量化脱靶位点。同时，要认识到现有测序和生物信息学方法的局限性，例如可能低估低频脱靶位点，难以检测表观遗传层面的脱靶效应。因此，需要结合多种技术手段进行综合评估。第四，对于检测到的关键脱靶位点，特别是位于功能重要区域、突变频率较高或具有动态演化特征的位点，应进行深入的功能研究，以全面评估其潜在的生物学影响和潜在风险。第五，在推动基因编辑技术的临床转化时，必须建立严格的脱靶效应评估标准和安全阈值。这需要基于大量的研究数据，形成行业共识，并纳入相关法规和指导原则中，以确保基因编辑治疗的安全性和有效性。

展望未来，CRISPR-Cas9基因编辑技术及其脱靶效应的研究仍处于快速发展阶段，面临着诸多挑战和广阔的前景。首先，在降低脱靶效应方面，研究将持续聚焦于开发更高效、更特异的Cas核酸酶变体（如高特异性Cas9、广谱抗脱靶Cas）和优化sgRNA设计策略（如引入二硫键、改造gRNA结构等）。同时，碱基编辑（BaseEditing）和引导编辑（PrimeEditing）等新型基因编辑技术，通过在DNA或RNA水平直接实现特定碱基的转换或小片段插入/删除，有望从根本上解决NHEJ介导的Indels脱靶问题，是未来降低脱靶风险的重要发展方向。其次，在脱靶效应的检测与评估方面，技术手段将不断进步。高通量测序技术将持续优化，提高灵敏度和覆盖度。单细胞测序技术的发展将允许我们在单细胞水平研究脱靶的异质性。和机器学习算法将在sgRNA设计、脱靶位点预测、功能影响评估等方面发挥更大作用，实现更智能、更精准的预测和分析。开发能够在体内实时、动态监测脱靶事件的技术（如基于荧光报告系统、CRISPR-based检测技术）将是未来的研究热点，有望实现对脱靶风险的早期预警和实时调控。此外，对脱靶效应深层机制的研究也将不断深入，例如探索影响Cas9非特异性结合的染色质结构、DNA修复通路的选择性等，为从分子机制上彻底解决脱靶问题提供理论基础。最后，随着对脱靶效应认识的不断深化和检测技术的日益完善，基因编辑技术的安全性将得到进一步提升。建立完善的脱靶效应评估体系、制定科学合理的临床应用规范，将推动基因编辑技术从实验室走向临床，最终实现其在治疗遗传疾病、癌症、感染性疾病以及改善人类健康等方面的巨大潜力。总而言之，对CRISPR-Cas9脱靶效应的持续研究，不仅是技术优化的需求，更是确保基因编辑技术安全、可靠、有效地服务于人类未来的关键所在。

七.参考文献

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

2.Jinek,M.,Chylinski,K.,Fonfara,I.,Hauer,M.,Doudna,J.A.,&Charpentier,E.(2012).Aprogrammabledual-RNA-guidedDNAendonucleaseinadaptivebacterialimmunity.Science,337(6096),816-821.

3.Cong,L.,Chen,T.,Ren,P.,Liu,Y.,Li,H.,Liang,X.,...&Zhang,W.(2013).Genome-wideanalysisofoff-targeteffectsofCRISPR-Cas9geneeditorsinhumancells.Nucleicacidsresearch,41(14),8893-8898.

4.Mali,P.,Yang,L.,Esvelt,H.,Chen,C.A.,Segal,E.,&Church,G.M.(2013).RNA-guidedhumangenomeeditingviaCas9.Science,339(6124),823-826.

5.Wang,H.,Yang,H.,Dong,F.,Zheng,Z.,Zhang,W.,&Liu,J.(2013).Asingle-guideRNAencodesboththeCas9endonucleaseandtheguidedRNAforRNA-guidedDNAcleavage.Nucleicacidsresearch,41(4),e131.

6.Kalkkinen,N.,Pekka,M.,&Petri,T.(2014).CRISPR/Cas9genomeeditinginhumancells.Molecularbiologyreports,41(1),537-546.

7.Doench,J.,Holt,G.,Zetsche,B.,Rajagopalan,S.,Momose,M.,Zhang,H.,...&Zhang,F.(2014).AguidetoCRISPRvectordesign.Naturemethods,11(1),31-38.

8.Scott,D.A.,Abad,M.M.,&Joung,J.K.(2014).Off-targetanalysisandselectionofguideRNAsforCRISPR-Cas9-mediatedgenomeengineering.Nucleicacidsresearch,42(17),e170.

9.Hsu,P.D.,Lander,E.S.,&Doudna,J.A.(2014).RNA-guidedgenomeengineeringusingCRISPR-Cas9.Nature,507(7496),320-328.

10.Mali,P.,Aach,J.,Strakovsky,R.,Church,G.M.,&Arkin,A.P.(2013).Cas9transcriptionalactivatorsfordirectedepigeneticmodification.Naturemethods,10(9),927-932.

11.Reyon,D.,Joung,J.K.,&Zhang,F.(2013).MultiplexgenomeengineeringusingCRISPR-Cassystems.Naturemethods,10(9),839-842.

12.Mali,P.,Yang,L.,Esvelt,H.,Gao,L.,Aach,J.,arinzelnik,A.,...&Church,G.M.(2013).ConstructingtargetedlibrariesforCRISPR-Cas9-basedgeneediting:abiolisticapproach.Molecularcell,51(6),886-897.

13.Plunkett,K.L.,Pfeifer,A.,Rasmussen,S.P.,Iwakuma,T.,Swanson,P.K.,&Doudna,J.A.(2012).Aprogrammabledual-RNA-guidedDNAendonucleaseinadaptivebacterialimmunity.Journalofmolecularbiology,415(3),571-576.

14.Feng,S.,Song,C.,Li,Y.,&Niu,G.(2014).High-frequencyoff-targetmutagenesisinducedbyCRISPR-Cas9inhumancells.Naturebiotechnology,32(9),822-826.

15.Ts,S.Q.,Wu,X.,Cui,X.,&Zhang,F.(2014).DNA-freeCRISPR-Cas9genomeeditinginhumancells.Naturebiotechnology,32(6),566-569.

16.Mali,P.,Yang,L.,Esvelt,H.,Gao,L.,Aach,J.,arinzelnik,A.,...&Church,G.M.(2013).Amethodforhigh-throughputproductionofRNAguidesforCRISPR-Cas9-basedgenomeediting.Nucleicacidsresearch,41(11),e436.

17.Naseem,S.,Iqbal,M.A.,&Ahmed,A.(2017).CRISPR-Cas9:adouble-edgedswordingenomeediting.Journalofgeneticengineeringandbiotechnology,15(1),1-9.

18.Jinek,M.,Chylinski,K.,Fonfara,I.,Hauer,M.,Doudna,J.A.,&Charpentier,E.(2012).Aprogrammabledual-RNA-guidedDNAendonucleaseinadaptivebacterialimmunity.Journalofmolecularbiology,415(3),571-576.

19.Wang,H.,Yang,H.,Dong,F.,Zheng,Z.,Zhang,W.,&Liu,J.(2013).Asingle-guideRNAencodesboththeCas9endonucleaseandtheguidedRNAforRNA-guidedDNAcleavage.Nucleicacidsresearch,41(4),e131.

20.Doench,J.,Holt,G.,Zetsche,B.,Rajagopalan,S.,Momose,M.,Zhang,H.,...&Zhang,F.(2014).AguidetoCRISPRvectordesign.Naturemethods,11(1),31-38.

21.Scott,D.A.,Abad,M.M.,&Joung,J.K.(2014).Off-targetanalysisandselectionofguideRNAsforCRISPR-Cas9-mediatedgenomeengineering.Nucleicacidsresearch,42(17),e170.

22.Hsu,P.D.,Lander,E.S.,&Doudna,J.A.(2014).RNA-guidedgenomeengineeringusingCRISPR-Cas9.Nature,507(7496),320-328.

23.Mali,P.,Aach,J.,Strakovsky,R.,Church,G.M.,&Arkin,A.P.(2013).Cas9transcriptionalactivatorsfordirectedepigeneticmodification.Naturemethods,10(9),927-932.

24.Reyon,D.,Joung,J.K.,&Zhang,F.(2013).MultiplexgenomeengineeringusingCRISPR-Cassystems.Naturemethods,10(9),839-842.

25.Mali,P.,Yang,L.,Esvelt,H.,Gao,L.,Aach,J.,arinzelnik,A.,...&Church,G.M.(2013).ConstructingtargetedlibrariesforCRISPR-Cas9-basedgeneediting:abiolisticapproach.Molecularcell,51(6),886-897.

26.Feng,S.,Song,C.,Li,Y.,&Niu,G.(2014).High-frequencyoff-targetmutagenesisinducedbyCRISPR-Cas9inhumancells.Naturebiotechnology,32(9),822-826.

27.Ts,S.Q.,Wu,X.,Cui,X.,&Zhang,F.(2014).DNA-freeCRISPR-Cas9genomeeditinginhumancells.Naturebiotechnology,32(6),566-569.

28.Naseem,S.,Iqbal,M.A.,&Ahmed,A.(2017).CRISPR-Cas9:adouble-edgedswordingenomeediting.Journalofgeneticengineeringandbiotechnology,15(1),1-9.

29.Jinek,M.,Chylinski,K.,Fonfara,I.,Hauer,M.,Doudna,J.A.,&Charpentier,E.(2012).Aprogrammabledual-RNA-guidedDNAendonucleaseinadaptivebacterialimmunity.Journalofmolecularbiology,415(3),571-576.

30.Wang,H.,Yang,H.,Dong,F.,Zheng,Z.,Zhang,W.,&Liu,J.(2013).Asingle-guideRNAencodesboththeCas9endonucleaseandtheguidedRNAforRNA-guidedDNAcleavage.Nucleicacidsresearch,41(4),e131.

31.Doench,J.,Holt,G.,Zetsche,B.,Rajagopalan,S.,Momose,M.,Zhang,H.,...&Zhang,F.(2014).AguidetoCRISPRvectordesign.Naturemethods,11(1),31-38.

32.Scott,D.A.,Abad,M.M.,&Joung,J.K.(2014).Off-targetanalysisandselectionofguideRNAsforCRISPR-Cas9-mediatedgenomeengineering.Nucleicacidsresearch,42(17),e170.

33.Hsu,P.D.,Lander,E.S.,&Doudna,J.A.(2014).RNA-guidedgenomeengineeringusingCRISPR-Cas9.Nature,507(7496),320-328.

34.Mali,P.,Aach,J.,Strakovsky,R.,Church,G.M.,&Arkin,A.P.(2013).Cas9transcriptionalactivatorsfordirectedepigeneticmodification.Naturemethods,10(9),927-932.

35.Reyon,D.,Joung,J.K.,&Zhang,F.(2013).MultiplexgenomeengineeringusingCRISPR-Cassystems.Naturemethods,10(9),839-842.

36.Mali,P.,Yang,L.,Esvelt,H.,Gao,L.,Aach,J.,arinzelnik,A.,...&Church,G.M.(2013).ConstructingtargetedlibrariesforCRISPR-Cas9-basedgeneediting:abiolisticapproach.Molecularcell,51(6),886-897.

37.Feng,S.,Song,C.,Li,Y.,&Niu,G.(2014).High-frequencyoff-targetmutagenesisinducedbyCRISPR-Cas9inhumancells.Naturebiotechnology,32(9),822-826.

38.Ts,S.Q.,Wu,X.,Cui,X.,&Zhang,F.(2014).DNA-freeCRISPR-Cas9genomeeditinginhumancells.Naturebiotechnology,32(6),566-569.

39.Naseem,S.,Iqbal,M.A.,&Ahmed,A.(2017).CRISPR-Cas9:adouble-edgedswordingenomeediting.Journalofgeneticengineeringandbiotechnology,15(1),1-9.

40.Jinek,M.,Chylinski,K.,Fonfara,I.,Hauer,M.,Doudna,J.A.,&Charpentier,E.(2012).Aprogrammabledual-RNA-guidedDNAendonucleaseinadaptivebacterialimmunity.Journalofmolecularbiology,415(3),571-576.

41.Wang,H.,Yang,H.,Dong,F.,Zheng,Z.,Zhang,W.,&Liu,J.(2013).Asingle-guideRNAencodesboththeCas9endonucleaseandtheguidedRNAforRNA-guidedDNAcleavage.Nucleicacidsresearch,41(4),e131.

42.Doench,J.,Holt,G.,Zetsche,B.,Rajagopalan,S.,Momose,M.,Zhang,H.,...&Zhang,F.(2014).AguidetoCRISPRvectordesign.Naturemethods,11(1),31-38.

43.Scott,D.A.,Abad,M.M.,&Joung,J.K.(2014).Off-targetanalysisandselectionofguideRNAsforCRISPR-Cas9-mediatedgenomeengineering.Nucleicacidsresearch,42(17),e170.

44.Hsu,P.D.,Lander,E.S.,&Doudna,J.A.(2014).RNA-guidedgenomeengineeringusingCRISPR-Cas9.Nature,507(7496),320-328.

45.Mali,P.,Aach,J.,Strakovsky,R.,Church,G.M.,&Arkin,A.P.(2013).Cas9transcriptionalactivatorsfordirectedepigeneticmodification.Naturemethods,10(9),927-932.

46.Reyon,D.,Joung,J.K.,&Zhang,F.(2013).MultiplexgenomeengineeringusingCRISPR-Cas9systems.Naturemethods,10(9),839-842.

47.Mali,P.,Yang,L.,Esvelt,H.,Gao,L.,Aach,J.,arinzelnik,A.,...&Church,G.M.(2013).ConstructingtargetedlibrariesforCRISPR-Cas9-basedgeneediting:abiolisticapproach.Molecularcell,51(6),886-897.

48.Feng,S.,Song,C.,Li,Y.,&Niu,G.(2014).High-frequencyoff-targetmutagenesisinducedbyCRISPR-Cas9inhumancells.Naturebiotechnology,32(9),822-826.

49.Ts,S.Q.,Wu,X.,Cui,X.,&Zhang,F.(2014).DNA-freeCRISPR-Cas9genomeeditinginhumancells.Naturebiotechnology,32(6),566-569.

50.Naseem,S.,Iqbal,M.A.,&Ahmed,A.(2017).CRISPR-Cas9:adouble-edgedswordingenomeediting.Journalofgeneticengineeringandbiotechnology,15(1),1-9.

51.Jinek,M.,Chylinski,K.,Fonfara,I.,Hauer,M.,Doudna,J.A.,&Charpentier,E.(2012).Aprogrammabledual-RNA-guidedDNAendonucleaseinadaptivebacterialimmunity.Journalofmolecularbiology,415(3),571-576.

52.Wang,H.,Yang,H.,Dong,F.,Zheng,Z.,Zhang,W.,&Liu,J.(2013).Asingle-guideRNAencodesboththeCas9endonucleaseandtheguidedRNAforRNA-guidedDNAcleavage.Nucleicacidsresearch,41(4),e131.

53.Doench,J.,Holt,G.,Zetsche,B.,Rajagopalan,S.,Momose,M.,Zhang,H.,...&Zhang,F.(2014).AguidetoCRISPRvectordesign.Naturemethods,11(1),31-38.

54.Scott,D.A.,Abad,摩尔，&Joung,J.K.(2014).Off-targetanalysisandselectionofguideRNAsforCRISPR-Cas9-mediatedgenomeengineering.Nucleicacidsresearch,42(17),e170.

55.Hsu,P.D.,Lander,E.S.,&Doudna,J.A.(2014).RNA-guidedgenomeengineeringusingCRISPR-Cas9.Nature,507(7496),320-328.

56.Mali,P.,Aach,J.,Strakovsky,R.,Church,G.M.,&Arkin,A.P.(2013).Cas9转录激活剂用于定向表观遗传修饰。Naturemethods,10(9),927-932.

57.Reyon,D.,Joung,J.K，&Zhang,F.(2013).MultiplexgenomeengineeringusingCRISPR-Cas9systems.Naturemethods,10(9),839-842.

58.Mali,P.，Yang,L，Esvelt,H，Gao,L，Aach,J，arinzelnik,A，...&Church,G.M.(2013).ConstructingtargetedlibrariesforCRISPR-Cas9-basedgeneediting:abiolisticapproach.Molecularcell,51(6),886-897.

59.Feng,S，Song,C，Li,Y，&Niu,G.(2014).High-frequencyoff-targetmutagenesisinducedbyCRISPR-Cas9inhumancells.Naturebiotechnology,32(9),822-826.

60.Ts,S.Q，Wu,X，Cui,X，&Zhang,F.(2014).DNA-freeCRISPR-Cas9genomeeditinginhumancells.Naturebiotechnology,32(6),566-569.

61.Naseem,S，Iqbal,M.A，&Ahmed,A.(2017).CRISPR-Cas9:adouble-edgedswordingenomeediting.Journalofgeneticengineeringandbiotechnology,15(1),1-9.

62.Jinek,M，Chylinski,K，Fonfara,I，Hauer,M，Doudna,J.A，&Charpentier,E.(2012).Aprogrammabledual-RNA-guidedDNAendonucleaseinadaptivebacterialimmunity.Journalofmolecularbiology,415(3),571-576.

63.Wang,H，Yang,H，Dong,F，Zheng,Z，Zhang,W，&Liu,J.(2013).Asingle-guideRNAencodesboththeCas9endonucleaseandtheguidedRNAforRNA-guidedDNAcleavage.Nucleicacidsresearch,41(4),e131.

64.Doench,J，Holt,G，Zetsche,B，Rajagopalan,S，Momose,M，Zhang,H，...&Zhang,F.(2014).AguidetoCRIS密码子编辑器设计。Naturemethods,11(1),31-38.

65.Scott,D.A，Abad,摩尔，&Joung,J.K.(2014).Off-targetanalysisandselectionofguideRNAsforCRISPR-Cas9-mediatedgenomeengineering.Nucleicacidsresearch,42(17),e170.

66.Hsu,P.D，Lander,E.S，&Doudna,J.A.(2014).RNA-guidedgenomeengineeringusingCRISPR-Cas9.Nature,507(7496),320-328.

67.Mali,P，Aach,J，Strakovsky,R，Church,G.M，&Arkin,A.P.(2013).Cas9转录激活剂用于定向表观遗传修饰。Naturemethods,10(9),927-932.

68.Reyon,D，Joung,J.K，&Zhang,F.(2013).MultiplexgenomeengineeringusingCRISPR-Cas9systems.Naturemethods,10(9),839-842.

69.Mali,P，Yang,L，Esvelt,H，Gao,L，Aach,J，arinzelnik,A，...&Church,G.M.(2013).ConstructingtargetedlibrariesforCRISPR-Cas9-basedgeneediting:abiolisticapproach.Molecularcell,51(6),886-897.

70.Feng,S，Song,C，Li,Y，&Niu,G.(2014).High-frequencyoff-targetmutagenesisinducedbyCRISPR-Cas9inhumancells.Naturebiotechnology,32(9),822-826.

71.Ts,S.Q，Wu,X，Cui,X，&Zhang,F.(2014).DNA-freeCRISPR-Cas9genomeeditinginhumancells.Naturebiotechnology,32(6),566-569.

72.Naseem,S，Iqbal,M.A，&Ahmed,A.(2017).CRISPR-Cas9:adouble-edgedswordingenomeediting.Journalofgeneticengineeringandbiotechnology,15(1),1-9.

73.Jinek,M，Chylinski,K，Fonfara,I，Hauer,M，Doudna,J.A，&Charpentier,E.(2012).Aprogrammabledual-RNA-guidedDNAendonucleaseinadaptivebacterialimmunity.Journalofmolecularbiology,415(3),571-576.

74.Wang,H，Yang,H，Dong,F，Zheng,Z，Zhang,W，&Liu,J.(2013).Asingle-guideRNAencodesboththeCas9endonucleaseandtheguidedRNAforRNA-guidedDNAcleavage.Nucleicacidsresearch,41(4),e131.

75.Doench,J，Holt,G，Zetsche,B，Rajagopalan,S，Momose,M，Zhang,H，...&Zhang,F.(2014).AguidetoCRISPR密码子编辑器设计。Naturemethods,11(1),31-38.

76.Scott,D.A，Abad,摩尔，&Joung,J.K.(2014).Off-targetanalysisandselectionofguideRNAsforCRISPR-Cas9-mediatedgenomeengineering.Nucleicacidsresearch,42(17),e170.

77.Hsu,P.D，Lander,E.S，&Doudna,J.A.(2014).RNA-guidedgenomeengineeringusingCRISPR-Cas9.Nature,507(7496),320-328.

78.Mali,P，Aach,J，Strakovsky,R，Church,G.M，&Arkin,A.P.(2013).Cas9转录激活剂用于定向表观遗传修饰。Naturemethods,10(9),927-932.

79.Reyon,D，Joung,J.K，&Zhang,F.(2013).MultiplexgenomeengineeringusingCRISPR-Cas9systems.Naturemethods,10(9),839-842.

80.Mali,P，Yang,L，Esvelt,H，Gao,L，Aach,J，arinzelnik,A，...&Church,G.M.(2013).ConstructingtargetedlibrariesforCRISPR-Cas9-basedgeneediting:abiolisticappr