Name: 减法基因组学的新序列发现 (视频) | JoVE Journal
Uploaded: 2019-01-25T13:30:00.000+00:00
Duration: 9 min 40 s
Description: 了解减法基因组学方法，利用免费软件从共纯化序列中识别难以分离的新基因组序列。该技术可应用于细菌疫苗靶点和未知病毒序列的鉴定，具有成本效益高、灵活性强的特点。

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

作者感谢 Michelle Biederman、Alyssa Pedersen 和 Colin J. Saldanha 在斑胸草雀基因组项目不同阶段提供的帮助。我们还感谢 Evgeny Bisk 在计算集群系统管理方面的支持，以及美国国立卫生研究院（NIH）提供的资助 1K22CA184297（授予 J.R.B.）和 NIH NS 042767（授予 C.J.S）。

1. 从头 组装起始序列
注意：任何下一代测序（NGS）数据均可使用，只要能够基于这些数据生成组装结果即可。合适的输入数据包括已组装成 fasta 文件的 Illumina、PacBio 或 Oxford Nanopore 测序读段。为具体说明，本节描述的是我们针对斑胸草雀研究所完成的基于 Illumina 的转录组组装7；但需注意，不同项目的具体细节会有所差异。在本示例项目中，原始数据来自 MiSeq 平台，每个样本大约获得 1000 万条成对读段。
<ol>
	<li>使用 Trimmomatic 0.3214 去除Illumina接头序列和低质量碱基。在命令行中输入： 
	java -jar trimmomatic-0.32.jar PE -phred33 正向.fq.gz 反向.fq.gz -baseout 质量与接头修剪 ILLUMINACLIP:TruSeq3-PE.fa:2:30:10 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:20 MINLEN:40</li>
	<li>使用 PEAR15 v. 0.9.6，使用默认参数对 trimmomatic 输出的成对读段进行高质量合并。在命令行中输入： 
	梨 -f &#60;quality_and_adaptor_trimmed_1P.fastq&#62; - r &#60;quality_and_adaptor_trimmed_2P.fastq&#62;</li>
	<li>使用 Reptile v. 1.116 对通过 PEAR 产生的测序读段进行纠错。遵循文中所述的逐步操作方案17.</li>
	<li>使用 Trinity v. 2.4.018 默认模式下组装校正后的序列。对于链特异性文库，使用 -SS_lib_type 参数。输出为一个 fasta 文件（your_assembly.fasta）。在命令行中输入： 
	Trinity --seqType fq --SS_lib_type FR –max_memory 10G –output Trinity_output --left quality_and_adaptor_trimmed_forward_paired_reads.fq –right quality_and_adaptor_trimmed_reverse_paired_reads.fq –CPU 10 
	注意：输出结果将保存在一个名为 Trinity_output 的新目录中，组装结果将被命名为“Trinity.fasta”，可根据需要重命名为 Your_assembly.fasta。更多详细信息请参见 Trinity 官方网站：https://github.com/trinityrnaseq/trinityrnaseq/wiki/Running-Trinity。</li>
</ol>
2. 将组装序列与参考序列进行BLAST比对
注意：当参考序列是组装序列或长读长序列（如 Sanger 测序结果）时，使用此步骤；如果参考序列由原始 Illumina 测序读段组成，请参见下方步骤 3，以将读段比对至查询序列。所有 BLAST 步骤均使用版本 2.2.29+ 完成，但这些命令应在任何较新的 BLAST 版本上均可正常运行。
<ol>
	<li>在命令行中为参考序列（nucleotide_reference.fasta）创建BLAST数据库。在命令行中输入以下命令： 
	makeblastdb -dbtype nucl -in nucleotide_reference.fasta -out nucleotide_reference.db</li>
	<li>将查询组装序列（在步骤1中生成）与参考数据库进行BLAST比对。若要生成输出文件，请使用 [-out BLAST_results.txt]；若要生成表格格式输出（后续使用Python脚本处理所必需），请使用 [-outfmt 6]。这些选项可按任意顺序组合，例如完整的命令为 [blastn -query your_assembly.fasta -db nucleotide_reference.db -out BLAST_results.txt -outfmt 6]。如果需要设置e值，请使用 -evalue 选项并指定合适的数值，例如 [-evalue 1e-6]。但请注意，如讨论部分所述，减法循环实际上会反转e值的筛选效果。</li>
	<li>为提高严谨性，建议使用组装序列中的蛋白质序列作为BLAST查询序列，并采用翻译核苷酸BLAST（tBLASTn），该方法会对（核苷酸）数据库进行六框翻译。此方法推荐用于大多数非模式生物系统，可避免蛋白质注释不完整的问题。
	<ol>
		<li>使用 -db_gencode 选项确保为所研究的生物选择正确的遗传密码。为获得查询用的蛋白质序列，运行TransDecoder包（v. 3.0.1）中的TransDecoder.LongOrfs命令，以识别组装查询序列中最长的开放阅读框。命令为 [TransDecoder.LongOrfs -t your_assembly.fasta]；输出将保存在名为 'transcripts.transdecoder_dir' 的目录中，其中包含一个名为 longest_orfs.pep 的文件，该文件含有来自 your_assembly.fasta 每条序列的最长预测蛋白质序列。</li>
		<li>使用tBLASTn时，运行命令 [tblastn -query longest_orfs.pep -db nucleotide_reference.db -out BLAST_results.txt -outfmt 6]。如果有高质量的蛋白质参考序列，建议使用BLASTp进行蛋白质-蛋白质比对，而非tBLASTn。</li>
		<li>为蛋白质参考序列创建BLAST数据库 [makeblastdb -dbtype prot -in protein_reference.fasta -out protein_reference.db]，然后执行 [blastp -query longest_orfs.pep -db protein_reference.db -out BLAST_results.txt -outfmt 6]。请务必保存结果文件以供下游分析使用，并采用表格格式（outfmt 6），确保Python脚本能够正确解析。</li>
	</ol>
	</li>
</ol>
3. 将测序读段比对到基因组组装结果
注意：如果参考数据集由原始基因组测序读段组成，而非组装后的序列或桑格测序序列，则可使用本方法；在此情况下，请使用 BLAST（步骤 2.1）。
<ol>
	<li>使用 BWA-MEM v. 0.7.1219 或 bowtie220，将下载的原始测序读段（raw_reads.fastq）比对到查询组装序列上。输出结果为 .sam 格式。操作命令如下：首先对组装序列建立索引：[bwa index your_assembly.fasta]，然后比对测序读段 [bwa mem your_assembly.fasta raw_reads.fastq &#62;mapped.sam]。（注意：此处的 '>' 符号并非大于号，而是指示将输出结果写入 mapped.sam 文件中）。</li>
</ol>
4. 使用 Python 脚本去除任何匹配序列
注意：提供的脚本适用于 Python 2.7。
<ol>
	<li>完成步骤2后，使用命令 [./Non-matching_sequences.py your_assembly.fasta BLAST_results.txt] 运行减法Python脚本。在运行脚本前，请确保BLAST输出文件为格式6（表格形式）。该脚本将输出一个包含不匹配序列的fasta格式文件，命名为your_assembly.fasta_non-matching_sequences_BLAST_results.txt.fasta，同时还会生成一个包含匹配序列的文件以供记录，命名为your_assembly.fasta_matching_sequences_BLAST_results.txt.fasta。其中，不匹配序列文件最为重要，可作为潜在T序列的来源，用于后续测试及进一步的减法基因组学循环。</li>
	<li>完成步骤3后，运行Python脚本removeUnmapped.py，输入为步骤3.1生成的.sam文件，该脚本将识别出没有任何匹配读段的查询序列名称，并将其保存至一个新的文本文件中。使用命令 [./removeUnmapped.py mapped.sam]，输出文件名为mapped.sam_contigs_with_no_reads.txt。（该程序还将生成一个已移除所有未比对读段的简化版sam文件；本实验方案中可忽略此文件，但其可能对其他分析有用。）</li>
	<li>由于上一步的输出是一个名为mapped.sam_contigs_with_no_reads.txt的文本文件，其中包含序列名称列表，因此需提取出这些序列以生成一个fasta文件：使用命令 [./getContig.py your_assembly.fasta mapped.sam_contigs_with_no_reads.txt]。输出文件将命名为mapped.sam_contigs_with_no_reads.txt.fasta。</li>
</ol>
5. 为剩余序列设计引物
注意：此时已获得一个包含候选 T 序列的 fasta 文件。本部分描述通过 qPCR 实验验证这些序列是来源于 T，还是来源于此前未知的 R 区域。如果在步骤 4 的扣除过程中移除了所有序列，则可能是因为初始组装未能包含 T，或者扣除条件可能过于严格。
<ol>
	<li>使用 Geneious21 手动确定最佳引物序列。
	<ol>
		<li>选择一段 21-28 bp 的候选序列作为正向引物。避免出现 4 个或更多相同碱基的连续序列。尽量选择碱基组成相对均匀的区域。在 3' 端包含单个 G 或 C 有助于引物的锚定。</li>
		<li>在屏幕右侧点击 Statistics 选项卡，当候选区域被选中时，查看该序列的估计熔解温度（Tm）。目标熔解温度应为 55-60 °C，同时避免重复序列和过长的 G/C 连续序列。</li>
		<li>按照步骤 5.1.1 和 5.1.2 选择反向引物，其位置应位于正向引物 3' 端下游 150-250 个碱基对处。虽然引物长度无需完全一致，但预测的 Tm 应尽可能接近正向引物的 Tm。务必对序列进行反向互补（在 Geneious 中选中序列后右键点击即可找到该选项）。</li>
	</ol>
	</li>
	<li>使用序列窗口顶部工具栏中的 Primer Design 功能。
	<ol>
		<li>点击 Primer Design 按钮，在 Target Region 下输入待扩增的区域。</li>
		<li>在 Characteristics 选项卡下，输入所需的产物大小、熔解温度（Tm）和 %GC 含量（参见步骤 5.1.1）。</li>
		<li>点击 OK 生成引物，并通过定制寡核苷酸服务订购引物。</li>
	</ol>
	</li>
	<li>使用对照 DNA（编码 T 和 R）验证引物，以优化 Tm 和延伸时间。可使用普通 Taq 酶结合凝胶电泳检测条带大小，也可按照步骤 6 的方法通过 qPCR 进行优化。
	<ol>
		<li>将正向和反向引物分别进行 10 倍稀释，使引物终浓度为 10 μM。</li>
		<li>PCR 反应体系为：0.5 μL dNTP、0.5 μL 正向引物、0.5 μL 反向引物、0.1 μL Taq 聚合酶、2 μL 模板、0.75 μL 镁离子、2.5 μL 缓冲液和 18.15 μL 水，每份反应总体积为 25 μL，模板浓度为 5 ng/μL。</li>
		<li>在 PCR 程序中测试不同熔解温度下的引物表现。通常最佳扩增效果出现在略低于引物预测 Tm 的温度，但一般不超过 60 °C。同时根据以下指南测试最佳延伸时间：每 1000 bp 延伸 1 分钟（因此通常根据扩增子长度选择 10-30 秒）。</li>
		<li>进行终点凝胶电泳，确认引物能扩增出预期序列。取 25 μL qPCR 产物与 5 μL 6X 甘油染料混合，在 2% TAE 琼脂糖凝胶上以 200 V 电压电泳 20 分钟。</li>
	</ol>
	</li>
</ol>
6.&#160;qPCR 验证剩余序列
注意：此步骤需要使用在步骤5中已验证的引物和建立的PCR条件。
<ol>
	<li>每个模板设置三个重复，反应体系如下：12.5 μL PowerSYBR Green 主混合液，0.5 μL 浓度为 10 μM 的正向引物，0.5 μL 浓度为 10 μM 的反向引物，10.5 μL 水，以及 1 μL 模板 DNA（浓度为 2 ng/μL），使每孔总体积为 25 μL。</li>
	<li>采用步骤 4 中已验证的退火温度和延伸时间设置 qPCR 程序。我们设计并验证了所有引物均适用于两步循环程序：95 °C 初始变性 10 分钟，随后进行 40 个循环（95 °C 变性 30 秒，60 °C 退火/延伸 1 分钟）。然而，对于某些引物，三步循环程序（变性-退火-延伸）可能更为优化，应根据需要进行调整。我们建议在首次使用引物进行 qPCR 时，至少进行一次最终的熔解曲线分析，以确认仅扩增出单一 DNA 产物。</li>
	<li>通过 Ct 值测定 qPCR/SYBR Green 信号相对于肌动蛋白（或任何其他合适的“R”对照）的表达水平。在所有情况下，计算 2-(目标基因 Ct - β-actin Ct) 的平均值和标准差。</li>
	<li>（可选）进行终点法凝胶电泳，以确认 qPCR 检测到的产物大小正确。具体操作为：将 25 μL qPCR 产物与 5 μL 6x 甘油加样缓冲液混合，在 2% TAE 琼脂糖凝胶上以 200 V 电压电泳 20 分钟。</li>
</ol>
7.&#160;使用新的参考样本进一步筛选数据。
注意：如果第6步已验证了从T中鉴定出的序列，则在此结束循环（图2A）。然而，多种因素可能促使继续该循环，例如当文件中仍存在大量R序列，或在第6步中未通过qPCR验证任何候选T序列时。
<ol>
	<li>获取新的参考序列。此步骤可启动循环的新一轮迭代，可能包括原始基因组数据、原始RNA测序数据或其他已组装的数据集。参考数据的有用资源包括美国国家生物技术信息中心（NCBI）的基因组数据库（https://www.ncbi.nlm.nih.gov/genome），该数据库通过FTP（ftp://ftp.ncbi.nlm.nih.gov/genomes/）提供可访问的已组装基因组；以及基因表达综合数据库（Gene Expression Omnibus，https://www.ncbi.nlm.nih.gov/geo/），其中存储了原始的下一代测序读段。基因组项目也可能通过其他与项目相关的网站和数据库提供其原始序列数据。</li>
</ol>

<ol><li>Barh, D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+Novel+Comparative+Genomics+Analysis+for+Common+Drug+and+Vaccine+Targets+in+Corynebacterium+pseudotuberculosis+and+other+CMN+Group+of+Human+Pathogens%5BTitle%5D)+AND+%22Chemical+Biology+%26+Drug+Design%22%5BJournal%5D)">A Novel Comparative Genomics Analysis for Common Drug and Vaccine Targets in Corynebacterium pseudotuberculosis and other CMN Group of Human Pathogens.</a> Chemical Biology & Drug Design. 78 (1), 73-84 (2011).</li>
<li>Sarangi, A. N., Aggarwal, R., Rahman, Q., Trivedi, N. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Subtractive+Genomics+Approach+for+in+Silico+Identification+and+Characterization+of+Novel+Drug+Targets+in+Neisseria+Meningitides+Serogroup+B%5BTitle%5D)+AND+%22Journal+of+Computer+Science+%26+Systems+Biology%22%5BJournal%5D)">Subtractive Genomics Approach for in Silico Identification and Characterization of Novel Drug Targets in Neisseria Meningitides Serogroup B.</a> Journal of Computer Science & Systems Biology. 2 (5), (2009).</li>
<li>Kaur, N., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Identification+of+Druggable+Targets+for+Acinetobacter+baumannii+Via+Subtractive+Genomics+and+Plausible+Inhibitors+for+MurA+and+MurB%5BTitle%5D)+AND+%22Applied+Biochemistry+and+Biotechnology%22%5BJournal%5D)">Identification of Druggable Targets for Acinetobacter baumannii Via Subtractive Genomics and Plausible Inhibitors for MurA and MurB.</a> Applied Biochemistry and Biotechnology. 171 (2), 417-436 (2013).</li>
<li>Rathi, B., Sarangi, A. N., Trivedi, N. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Genome+subtraction+for+novel+target+definition+in+Salmonella+typhi%5BTitle%5D)+AND+%22Bioinformation%22%5BJournal%5D)">Genome subtraction for novel target definition in Salmonella typhi.</a> Bioinformation. 4 (4), 143-150 (2009).</li>
<li>Epstein, J. H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Identification+of+GBV-D%2C+a+Novel+GB-like+Flavivirus+from+Old+World+Frugivorous+Bats+%28Pteropus+giganteus%29+in+Bangladesh%5BTitle%5D)+AND+%22PLoS+Pathogens%22%5BJournal%5D)">Identification of GBV-D, a Novel GB-like Flavivirus from Old World Frugivorous Bats (Pteropus giganteus) in Bangladesh.</a> PLoS Pathogens. 6 (7), (2010).</li>
<li>Kapoor, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Identification+of+Rodent+Homologs+of+Hepatitis+C+Virus+and+Pegiviruses%5BTitle%5D)+AND+%22MBio%22%5BJournal%5D)">Identification of Rodent Homologs of Hepatitis C Virus and Pegiviruses.</a> MBio. 4 (2), (2013).</li>
<li>Biederman, M. K., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Discovery+of+the+First+Germline-Restricted+Gene+by+Subtractive+Transcriptomic+Analysis+in+the+Zebra+Finch%2C+Taeniopygia+guttata%5BTitle%5D)+AND+%22Current+Biology%22%5BJournal%5D)">Discovery of the First Germline-Restricted Gene by Subtractive Transcriptomic Analysis in the Zebra Finch, Taeniopygia guttata.</a> Current Biology. 28 (10), 1620-1627 (2018).</li>
<li>Readhead, B., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Multiscale+Analysis+of+Independent+Alzheimer%27s+Cohorts+Finds+Disruption+of+Molecular%2C+Genetic%2C+and+Clinical+Networks+by+Human+Herpesvirus%5BTitle%5D)+AND+%22Neuron%22%5BJournal%5D)">Multiscale Analysis of Independent Alzheimer's Cohorts Finds Disruption of Molecular, Genetic, and Clinical Networks by Human Herpesvirus.</a> Neuron. 99, 1-19 (2018).</li>
<li>Carroll, D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+global+virome+project%5BTitle%5D)+AND+%22Science%22%5BJournal%5D)">The global virome project.</a> Science. 359 (6378), 872-874 (2016).</li>
<li>Pigozzi, M. I., Solari, A. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Germ+cell+restriction+and+regular+transmission+of+an+accessory+chromosome+that+mimics+a+sex+body+in+the+zebra+finch.+Taeniopygia+guttata%5BTitle%5D)+AND+%22Chromosome+Research%22%5BJournal%5D)">Germ cell restriction and regular transmission of an accessory chromosome that mimics a sex body in the zebra finch. Taeniopygia guttata.</a> Chromosome Research. 6, 105-113 (1998).</li>
<li>Itoh, Y., Kampf, K., Pigozzi, M. I., Arnold, A. P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Molecular+cloning+and+characterization+of+the+germline-restricted+chromosome+sequence+in+the+zebra+finch%5BTitle%5D)+AND+%22Chromosoma%22%5BJournal%5D)">Molecular cloning and characterization of the germline-restricted chromosome sequence in the zebra finch.</a> Chromosoma. 118, 527-536 (2009).</li>
<li>Warren, W. C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+genome+of+a+songbird%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">The genome of a songbird.</a> Nature. 464, 757-762 (2010).</li>
<li>Balakrishnan, C. N., Lin, Y. C., London, S. E., Clayton, D. F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((RNAseq+transcriptome+analysis+of+male+and+female+zebra+finch+cell+lines%5BTitle%5D)+AND+%22Genomics%22%5BJournal%5D)">RNAseq transcriptome analysis of male and female zebra finch cell lines.</a> Genomics. 100, 363-369 (2012).</li>
<li>Bolger, A. M., Lohse, M., Usadel, B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Trimmomatic%3A+a+flexible+trimmer+for+Illumina+sequence+data%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">Trimmomatic: a flexible trimmer for Illumina sequence data.</a> Bioinformatics. 30 (15), 2114-2120 (2014).</li>
<li>Zhang, J., Kobert, K., Flouri, T., Stamatakis, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((PEAR%3A+a+fast+and+accurate+Illumina+Paired-End+reAd+mergeR%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">PEAR: a fast and accurate Illumina Paired-End reAd mergeR.</a> Bioinformatics. 30, 614-620 (2014).</li>
<li>Yang, X., Dorman, K. S., Aluru, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Reptile%3A+representative+tiling+for+short+read+error+correction%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">Reptile: representative tiling for short read error correction.</a> Bioinformatics. 26, 2526-2533 (2010).</li>
<li>MacManes, M. D., Eisen, M. B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Improving+transcriptome+assembly+through+error+correction+of+high-throughput+sequence+reads%5BTitle%5D)+AND+%22PeerJ%22%5BJournal%5D)">Improving transcriptome assembly through error correction of high-throughput sequence reads.</a> PeerJ. 1 (113), (2013).</li>
<li>Grabherr, M. G., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Full-length+transcriptome+assembly+from+RNA-seq+data+without+a+reference+genome%5BTitle%5D)+AND+%22Nature+Biotechnology%22%5BJournal%5D)">Full-length transcriptome assembly from RNA-seq data without a reference genome.</a> Nature Biotechnology. 29, 644-652 (2011).</li>
<li>Li, H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Aligning+sequence+reads%2C+clone+sequences+and+assembly+contigs+with+BWA-MEM%5BTitle%5D)+AND+%22arXiv%22%5BJournal%5D)">Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM.</a> arXiv. , (2013).</li>
<li>Langmead, B., Salzberg, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Fast+gapped-read+alignment+with+Bowtie+2%5BTitle%5D)+AND+%22Nature+Methods%22%5BJournal%5D)">Fast gapped-read alignment with Bowtie 2.</a> Nature Methods. 9, 357-359 (2012).</li>
<li>Kearse, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Geneious+Basic%3A+an+integrated+and+extendable+desktop+software+platform+for+the+organization+and+analysis+of+sequence+data%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data.</a> Bioinformatics. 28 (12), 1647-1649 (2012).</li>
<li>Peirson, S. N., Butler, J. N. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Quantitative+polymerase+chain+reaction%5BTitle%5D)+AND+%22Methods+in+Molecular+Biology%22%5BJournal%5D)">Quantitative polymerase chain reaction.</a> Methods in Molecular Biology. 362, 349-362 (2007).</li>
<li>Hunt, M., Kikuchi, T., Sanders, M., Newbold, C., Berriman, M., Otto, T. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((REAPR+citation%3A+REAPR%3A+a+universal+tool+for+genome+assembly+evaluation%5BTitle%5D)+AND+%22Genome+Biology%22%5BJournal%5D)">REAPR citation: REAPR: a universal tool for genome assembly evaluation.</a> Genome Biology. 14 (5), (2013).</li>
<li>Meyer, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+mitochondrial+genome+sequence+of+a+hominin+from+Sima+de+los+Huesos%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">A mitochondrial genome sequence of a hominin from Sima de los Huesos.</a> Nature. 505 (7483), 403-406 (2013).</li>
<li>Gunnarsdóttir, E. D., Li, M., Bauchet, M., Finstermeier, K., Stoneking, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((High-throughput+sequencing+of+complete+human+mtDNA+genomes+from+the+Philippines%5BTitle%5D)+AND+%22Genome+Research%22%5BJournal%5D)">High-throughput sequencing of complete human mtDNA genomes from the Philippines.</a> Genome Research. 21 (1), 1-11 (2010).</li>
<li>King, J. L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((High-quality+and+high-throughput+massively+parallel+sequencing+of+the+human+mitochondrial+genome+using+the+Illumina+MiSeq%5BTitle%5D)+AND+%22Forensic+Science+International%3A+Genetics%22%5BJournal%5D)">High-quality and high-throughput massively parallel sequencing of the human mitochondrial genome using the Illumina MiSeq.</a> Forensic Science International: Genetics. 12, 128-135 (2014).</li>
<li>Yao, X., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+First+Complete+Chloroplast+Genome+Sequences+in+Actinidiaceae%3A+Genome+Structure+and+Comparative+Analysis%5BTitle%5D)+AND+%22Plos+One%22%5BJournal%5D)">The First Complete Chloroplast Genome Sequences in Actinidiaceae: Genome Structure and Comparative Analysis.</a> Plos One. 10 (6), (2015).</li>
<li>Zhang, Y., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+Complete+Chloroplast+Genome+Sequences+of+Five+Epimedium+Species%3A+Lights+into+Phylogenetic+and+Taxonomic+Analyses%5BTitle%5D)+AND+%22Frontiers+in+Plant+Science%22%5BJournal%5D)">The Complete Chloroplast Genome Sequences of Five Epimedium Species: Lights into Phylogenetic and Taxonomic Analyses.</a> Frontiers in Plant Science. 7, (2016).</li>
<li>Swart, E. C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+Oxytricha+trifallax+Mitochondrial+Genome%5BTitle%5D)+AND+%22Genome+Biologyogy+and+Evolution%22%5BJournal%5D)">The Oxytricha trifallax Mitochondrial Genome.</a> Genome Biologyogy and Evolution. 4 (2), 136-154 (2011).</li>
<li>Barth, D., Berendonk, T. U. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+mitochondrial+genome+sequence+of+the+ciliate+Paramecium+caudatum+reveals+a+shift+in+nucleotide+composition+and+codon+usage+within+the+genus+Paramecium%5BTitle%5D)+AND+%22BMC+Genomics%22%5BJournal%5D)">The mitochondrial genome sequence of the ciliate Paramecium caudatum reveals a shift in nucleotide composition and codon usage within the genus Paramecium.</a> BMC Genomics. 12 (1), (2011).</li>
<li>Coombe, L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Assembly+of+the+Complete+Sitka+Spruce+Chloroplast+Genome+Using+10X+Genomics%27+GemCode+Sequencing+Data%5BTitle%5D)+AND+%22Plos+One%22%5BJournal%5D)">Assembly of the Complete Sitka Spruce Chloroplast Genome Using 10X Genomics' GemCode Sequencing Data.</a> Plos One. 11 (9), (2016).</li>
<li>Herschleb, J., Ananiev, G., Schwartz, D. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Pulsed-field+gel+electrophoresis%5BTitle%5D)+AND+%22Nature+Protocols%22%5BJournal%5D)">Pulsed-field gel electrophoresis.</a> Nature Protocols. 2 (3), 677-684 (2007).</li>
</ol>

作者无任何利益冲突需要披露。

尽管减法基因组学功能强大，但它并非一种千篇一律的方法，需要在多个关键步骤上进行定制化处理，并仔细选择参考序列和测试样本。如果查询组装的质量较差，过滤步骤可能仅能分离出组装过程中的假象。因此，使用适合特定项目的验证方案对从头（de novo）组装结果进行充分验证至关重要。对于RNA-seq数据，Trinity网站提供了相关指南18；而对于DNA数据，则可使用如REAPR23之类的工具。使用BLAST时另一个关键步骤是选择合适的e值，这将决定减法过程是宽松还是严格。然而，该方法中存在一个反向逻辑：与参考序列的匹配越严格，实际上意味着减法越不严格，因为未匹配的序列不会被剔除。因此，为了实现更严格的减法，应在BLAST中使用更大的（即更宽松的）e值。本方案最后的关键步骤是参考序列的选择。为了达到最高效率，参考序列应尽可能完整；但其并不需要完美，因为qPCR检测可以确认剩余序列是否来源于T或R，以及是否需要进一步过滤。在实施该方案过程中，可引入新的参考序列以进一步缩小待验证基因的范围。我们注意到，有时匹配方法可能会发生变化：在最后一次...

该方法的目的是从基因组背景或参考序列（R）中鉴定出新的目标（T）基因组序列，该序列可以是DNA或RNA（图1）。当目标序列无法通过物理手段分离，或分离成本较高时，该方法尤为适用。目前仅有少数生物具有可用于消减分析的完整基因组，因此本方法的关键创新之处在于将计算分析与实验操作相结合，形成一个循环流程，使研究人员能够在参考基因组不完整或仅为非模式生物的草图基因组时，仍能分离出目标序列。每个循环结束后，通过qPCR检测判断是否需要进一步进行消减。经验证的候选T序列应在已知T阳性样本中表现出显著更高的qPCR检出信号。
该方法已被应用于发现不存在宿主同源物的新型细菌药物靶点1,2,3,4，以及从受感染宿主中鉴定新型病毒5,6。除了可用于鉴定目标序列T外，该方法还可提升参考基因组R的完整性：我们近期利用该方法从斑胸草雀参考基因组中鉴定出936个缺失基因，并发现了一个位于生殖系特异性染色体（T）上的新基因7。当目标序列T与已知序列高度分歧，或其身份尚不明确时（例如斑胸草雀的生殖系限制性染色体），减法基因组学尤其具有价值7。
减法基因组学的一个关键优势在于其无偏性，即无需事先对T进行阳性鉴定。在最近的一项研究中，Readhead et al. 探讨了阿尔茨海默病与四种脑区中病毒丰度之间的关系。在病毒鉴定过程中，Readhead et al. 构建了一个包含515种病毒的数据库8，这严重限制了该研究可识别的病毒种类。若采用减法基因组学方法，则可比较健康个体与阿尔茨海默病患者的基因组，从而分离出可能与该疾病相关的新病毒，而无需考虑这些病毒与已知病原体的相似性。目前已知有263种可感染人类的病毒，但据估计，尚有约167万种未被发现的病毒物种存在，其中约63.1万至82.7万种具有感染人类的潜力9。
分离新病毒是减法基因组学特别有效的研究领域，但某些研究可能无需如此严格的方法。例如，已有研究通过非偏向性高通量测序结合逆转录和BLASTx分析来鉴定新病毒序列。5 或富集病毒核酸以提取并反转录病毒序列6. 尽管这些研究采用了 从头合成 测序与组装过程中未使用减法法，因为目标序列已通过BLAST得到明确鉴定。若病毒为完全新型且与其他病毒无关联或仅有远缘关系，则减法基因组学将是一种有用的技术。减法基因组学的优势在于能够获得全新的序列。若宿主生物的基因组已知，可将其序列扣除，从而保留所有病毒序列。例如，在我们已发表的研究中，尽管并非最初目的，但仍通过减法基因组学从斑胸草雀中分离出一条 novel 病毒序列7.
减法基因组学在细菌疫苗靶点的鉴定中也显示出重要价值，这主要是由于抗生素耐药性的急剧上升1,2,3,4。为了最大限度地降低自身免疫反应的风险，研究人员通过剔除在人类宿主中具有同源蛋白的任何蛋白，进一步缩小潜在疫苗靶点的范围。一项针对Corynebacterium pseudotuberculosis的研究中，通过将多种细菌基因组与脊椎动物宿主基因组进行比对并剔除共有序列，以确保可能的药物靶点不会影响宿主中的蛋白，从而避免不良反应1。这类研究的基本工作流程包括：下载细菌的蛋白质组，确定关键蛋白，去除冗余蛋白，利用BLASTp筛选出必需蛋白，并通过与宿主蛋白质组进行BLASTp比对，剔除任何具有宿主同源性的蛋白1,2,3,4。在此背景下，减法基因组学可确保所开发的疫苗在宿主体内不会产生脱靶效应1,2,3,4。
我们利用减法基因组学鉴定了生殖系限制性染色体（GRC）（本例中为T染色体）上的首个蛋白质编码基因，该染色体存在于两性的生殖系中，但不存在于体细胞组织中10。在此项研究之前，关于GRC唯一已知的基因组信息是一个重复区域11。我们对成年斑胸草雀的卵巢和睾丸组织（R+T）进行RNA测序，并进行了从头转录组组装（de novo assembly）。通过比对已发表的体细胞（肌肉）基因组序列（R1）12、其原始（Sanger）测序读段数据（R2）以及体细胞（脑组织）转录组数据（R3）13，在计算层面逐步剔除共有序列。三个参考数据集的依次使用是基于每个循环第5步qPCR检测结果（图2A）的提示，表明需要进一步过滤。所发现的α-SNAP基因通过DNA和RNA的qPCR、克隆及测序得到了验证。本示例表明，该方法具有良好的灵活性：不依赖于核酸类型（DNA与RNA）的匹配，且减法过程可使用由基因组组装结果或原始测序读段构成的参考数据集（R）。

<table><tbody><tr><td>Accustart II Taq DNA 聚合酶</td><td>Quanta Bio</td><td>95141</td><td></td></tr><tr><td>BLAST（局部序列比对搜索工具）</td><td></td><td></td><td>https://github.com/trinityrnaseq/trinityrnaseq/wiki/Transcriptome-Assembly-Quality-Assessment</td></tr><tr><td>Bowtie 2</td><td></td><td></td><td>https://www.python.org/download/releases/2.7/</td></tr><tr><td>BWA-MEM v. 0.7.12</td><td></td><td></td><td>https://github.com/BenLangmead/bowtie2</td></tr><tr><td>Geneious</td><td></td><td></td><td>https://blast.ncbi.nlm.nih.gov/Blast.cgi</td></tr><tr><td>PEAR v. 0.9.6</td><td></td><td></td><td>http://www.mybiosoftware.com/reptile-1-1-short-read-error-correction.html</td></tr><tr><td>个人计算机</td><td>Biomatters</td><td></td><td>http://www.geneious.com/</td></tr><tr><td>PowerSYBR qPCR 混合液</td><td>ThermoFisher</td><td>4367659</td><td></td></tr><tr><td>Python v. 2.7</td><td></td><td></td><td>https://sco.h-its.org/exelixis/web/software/pear/</td></tr><tr><td>Reptile v.1.1</td><td></td><td></td><td>https://alurulab.cc.gatech.edu/reptile</td></tr><tr><td>Stratagene Mx3005P</td><td>Agilent Technologies</td><td>401456</td><td></td></tr><tr><td>TransDecoder v. 3.0.1</td><td></td><td></td><td>https://sourceforge.net/projects/bio-bwa/files/</td></tr><tr><td>Trinity v. 2.4.0</td><td></td><td></td><td>https://github.com/TransDecoder/TransDecoder/wiki</td></tr></tbody></table>

通过减法基因组学进行新序列发现

Novel Sequence Discovery by Subtractive Genomics

在旨在鉴定嵌入于较大基因组背景中的基因、蛋白质或一般区域序列的任何研究中，均可使用减法基因组学。减法基因组学使研究人员能够通过全面测序并去除已知遗传元件（参考序列，R），从而分离出目标序列（T）。该方法可用于鉴定线粒体、叶绿体、病毒或生殖系限制性染色体等新序列，尤其适用于T难以从R中轻易分离的情况。该方法从完整的基因组数据（R + T）出发，利用基本局部比对搜索工具（BLAST）比对参考序列或多个参考序列，以去除匹配的已知序列（R），从而留下目标序列（T）。为了使减法效果最佳，R应为相对完整的草图，且不包含T。由于减法后剩余的序列将通过定量聚合酶链式反应（qPCR）进行验证，因此R无需完全完整，该方法仍可有效运行。本文将计算步骤与实验步骤结合为一个可重复迭代的循环过程，依次去除多个参考序列，并逐步优化对T的搜索。减法基因组学的优势在于，即使在物理纯化困难、不可能或成本高昂的情况下，也能鉴定出完全新颖的目标序列。该方法的局限性在于难以找到合适的参考序列，以及获取用于qPCR检测的T阳性与阴性样本。我们在此描述了该方法在鉴定斑胸草雀生殖系限制性染色体上首个基因中的具体应用。在该案例中，计算过滤涉及三个参考序列（R），在三个循环中依次去除：一个不完整的基因组组装、原始基因组数据以及转录组数据。

运行 BLAST 后，输出文件将包含查询序列中与数据库匹配的序列列表。经过 Python 减法分析后，将获得一批不匹配的序列，并通过 qPCR 进行验证。相关结果及后续步骤如下所述。
阴性结果。在将序列比对到参考序列后，可能出现两种阴性结果。一种情况是没有比对结果，意味着该总序列与参考序列之间没有任何相似序列。这可能是由于为测序样本选择了错误的参考序列所致。另一种可能是起始组装中不存在独特序列（所有序列均被扣除），因此无法找到目标序列对应的基因。请核查参考序列的来源，并确保其并非来自与查询组装相同的组织。
经过计算过滤后，qPCR 可能产生阴性结果，示例见图 3A、3B、C，其中在不同鸟类组织间的检测结果无...

本实验方案的目的是结合计算研究与实验研究，寻找那些难以从共纯化序列中分离出来的新型序列，这些共纯化序列可能仅有部分被知晓。

减法基因组学可用于任何旨在鉴定嵌入于……中的基因、蛋白质或一般区域序列的研究。

本方案的意义在于可用于鉴定那些无法从共纯化序列中分离出来的基因组序列，而这些共纯化序列本身可能也仅有部分信息被知晓。该技术的主要优势在于成本低廉，主要使用可免费下载的软件，且具有高度灵活性。该方法可应用于多种生物学问题的研究。潜在应用包括鉴定细菌疫苗靶点，以及鉴定序列与已知微生物差异极大的病毒。只要能够获得不含基因组靶标的参考序列，该方法即可应用于任何无法通过实验手段将未知成分分离的体系。该技术需要一定的尝试与调整过程，因此耐心至关重要。您可能需要对程序进行故障排查。您可以使用与本文所述不同的程序。请尽可能参考用户手册。该方法的可视化演示非常有帮助，因为计算工作依赖于对命令行编程结构的基本理解。首先，使用 Trimmomatic 0.32 去除 illumina 接头序列和低质量碱基。然后，使用 Pear 0.9.11 版本，以默认参数将 Trimmomatic 输出的成对读段合并，生成高质量的拼接读段。然后使用 Reptile 1.1 版本对通过 Pear 产生的测序读段进行纠错。最后在默认模式下使用 Trinity 2.4.0 版本组装纠错后的序列。对于链特异性文库，需使用 SS_lib_type 参数。输出为一个 FASTA 文件，将被保存至名为 trinity_output 的新目录中。在命令行中为参考序列 nucleotide_reference.fasta 创建一个 BLAST 数据库。BLAST 将查询序列与参考数据库进行比对。要获得输出文件，请使用 BLAST_results.txt。若需生成后续 Python 脚本处理步骤所需的表格格式输出，请使用 outfmt 6。为提高比对的严格性，可使用组装得到的蛋白质序列作为 BLAST 查询序列，并采用翻译核酸 BLAST（translated nucleotide BLAST），后者会对核酸数据库进行六框翻译。为了获得查询序列的蛋白质序列，运行 TransDecoder.LongOrfs 命令以识别从组装好的查询序列中得到的最长开放阅读框。然后运行 tblastn。如有必要，使用 db_gencode 并选择适当的编码选项，确保为所研究的生物体选择正确的遗传密码。如果有高质量的蛋白质参考序列可用，应使用 blastp 进行蛋白质-蛋白质比对，而不是使用 tblastn。建立该蛋白质参考序列的 BLAST 数据库。务必将结果保存为文件以供后续分析，并使用表格形式输出，以确保 Python 脚本能够正确解析结果。现在，使用减法 Python 脚本去除任何匹配的序列。为了将测序读段比对到组装结果上，使用 BWA-MEM Version 0.7.12 或 bowtie 2 将下载的原始测序读段比对到查询组装序列上。首先，对组装序列建立索引，然后进行读段比对。输出结果为 SAM 格式。使用 SAM 文件作为输入，运行 Python 脚本 removeUnmapped.py。该步骤可识别出没有匹配读段的查询序列名称，并将其保存至一个新的文本文件中。上一步骤的输出结果是一个包含序列名称的 txt 文件。请从中提取这些序列，生成一个 FASTA 文件。输出结果将是一个fasta文件。使用Genius手动确定最佳引物序列。为正向引物选择一段21至28个碱基对的候选序列，避免出现四个或更多相同碱基的连续序列。尽量选择一个碱基对组成相对均匀的区域。在3'端有一个单独的G或C是有利的，有助于引物的结合。在屏幕右侧点击“Statistics”标签，当候选区域被高亮显示时，可查看该序列的估计熔解温度。引物的熔解温度应控制在55至60摄氏度之间，避免出现重复序列以及连续的G或C碱基。反向引物的设计方法相同，应位于正向引物3'端下游150至250个碱基对范围内。引物长度无需完全一致，但预测的熔解温度应尽可能与正向引物接近。在Genius软件中选中序列后，通过右键菜单中的选项对序列进行反向互补操作。另一种方法是使用序列窗口顶部工具栏中的引物设计（Primer Design）功能。在“目标区域”（Target Region）下插入需要扩增的区域。在“特性”（characteristics）选项卡下，输入所需的引物长度、熔解温度（melting temperature）以及 G-C 含量百分比（per cent G C），然后点击“确定”（OK）以生成引物。为了对剩余序列进行定量PCR验证，首先为每个模板准备三份反应混合液，包含SYBR Green Master Mix、正向引物和反向引物以及水，总体积为25 µL。根据先前已验证的退火温度和延伸时间设置qPCR扩增程序。首次使用引物进行qPCR时，至少应生成最终的熔解曲线，以确认扩增产物为单一DNA片段。通过Ct值将qPCR SYBR Green信号相对于肌动蛋白（Actin）进行测量。在所有情况下，计算相对于肌动蛋白的2的Ct次方的平均值和标准偏差。进行终点凝胶电泳，以确认qPCR检测到的产物大小正确。在此，将25微升qPCR产物与5微升6X Glitterol染料混合，于2% TAE琼脂糖凝胶中在200伏电压下电泳20分钟。如果您的 qPCR 结果显示未鉴定到目标序列，请使用新的参考序列重复整个流程，新参考序列可从在线数据库中获取。本案例中的消减文库项目始于对成年雄性和雌性斑胸草雀生殖系组织 RNA 的测序。最终，共鉴定出 935 个此前未被纳入全基因组注释的体细胞基因。经过计算过滤后，定量 PCR 可能产生阴性结果，即在不同鸟类组织间的检测无差异。相反，当基因组 DNA 的 qPCR 显示目标组织中的检测值相对于参考样本具有统计学上显著更高的水平时，则确认为阳性结果，代表真实目标序列被成功鉴定。在此案例中，alpha snap 基因被验证为生殖系限制性基因，因其在体细胞组织中的含量相对于睾丸 DNA 显著降低，而在睾丸 DNA 中其含量与肌动蛋白（Actin）相当。在进行每一步计算时，务必注意使用正确的输入数据。可能需要多次重复循环减法步骤，才能获得目标序列或多个目标序列。可对发现的基因进行多种系统发育、结构和功能分析。这些附加方法有助于深入了解基因的进化和功能作用。我们鉴定了鸣禽生殖系限制性染色体上的首个基因，这一发现使人们对这种出人意料的基因组元件产生了更广泛的关注。后续研究显示，在多种鸣禽物种中存在类似的染色体，且这些染色体含有大量其他基因。

novel-sequence-discovery-by-subtractive-genomics

Research

科研

JoVE Journal

Genetics

遗传学

减法基因组学的新序列发现

我们已开发出全面且无偏倚的全基因组筛选方法，以研究基因-药物和基因-环境之间的相互作用。本文介绍了对这些突变体文库进行筛选的方法。...

条形码酵母文库的竞争性基因组筛选

混合DNA测序是一种快速且经济高效的策略，可用于在大样本队列中检测与复杂表型相关的罕见变异。本文介绍了利用SPLINTER软件包对32个癌症相关基因进行混合的下一代测序数据进行计算分析的方法。该方法具有可扩展性，适用于任何感兴趣的表型研究。...

使用 SPLINTER 检测池化测序中的稀有基因组变异

表征微生物群落一直是环境微生物学中的长期目标。新一代测序技术如今能够以极低的成本和工作量，以前所未有的深度对微生物群落进行表征。本文详细介绍了我们使用台式测序仪对细菌16S核糖体RNA基因进行测序的方法。...

16S核糖体RNA基因扩增子的高通量测序

亚硫酸氢盐扩增子测序（BSAS）是一种用于定量目标基因组区域中胞嘧啶甲基化水平的方法。该方法在下一代测序之前，结合亚硫酸氢盐转化与目标区域的PCR扩增，从而在单碱基水平上实现DNA甲基化的绝对定量。...

通过下一代测序进行靶向DNA甲基化分析

我们描述了一种将流式细胞术与高通量测序相结合的技术，用于鉴定基因组中晚期复制的区域。...

G2-seq：一种基于高通量测序技术鉴定基因组晚期复制区域的方法

下一代测序（NGS）是一种强大的基因组表征工具，但其应用受到该平台较高错误率（约0.5–2.0%）的限制。我们介绍了错误校正测序的方法，可消除NGS的错误率影响，检测到低至0.0001变异等位基因分数的突变。...

利用纠错型DNA和RNA测序进行稀有事件检测

鉴定基因与调控元件之间的物理相互作用具有挑战性，但染色体构象捕获技术为此提供了便利。该4C-seq方案的改进通过减少PCR模板的过度扩增，减轻了PCR偏好性，同时通过增加限制性内切酶消化步骤，提高了测序读段的可比对性。...

利用下一代测序技术对染色体构象捕获环化文库（4C-seq）进行基因调控序列的高通量鉴定

本文建立并描述了一种由质谱引导的基因组挖掘方案。该方案基于基因组序列信息和液相色谱-串联质谱（LC-MS/MS）分析，旨在促进从复杂的微生物和植物提取物中鉴定分子。...

质谱引导的基因组挖掘作为发现新型天然产物的工具

基于转录组学的低成本药物筛选

在本文中，我们描述了一种用于功能性检测转录过程以及针对细菌转录的抗菌剂抑制活性的实验方案。...

体外 转录实验及其在药物发现中的应用

了解减法基因组学方法，利用免费软件从共纯化序列中识别难以分离的新基因组序列。该技术可应用于细菌疫苗靶点和未知病毒序列的鉴定，具有成本效益高、灵活性强的特点。

通过减法基因组学进行新序列发现

探索更多视频