上海发布
判断“人or猪or狗”最可靠的方式,不是观察文件名、样本名称或单条基因,而是对经过质控的测序数据分别比对到人、猪、狗的同版本参考基因组,再结合唯一比对比例、覆盖范围、线粒体与核基因组信号、污染情况共同判断。DNA数据下载则应先确认样本或参考基🤔因组的🎆 accession、文件类型、基因组版本和校验值,再下载并验证文件完整性。
质控工具的输出应保存为独立报告。fastp、FastQC等工具可以⚡帮助发现质量问题,但它们只能说明序列是否适合分析,不能直接回答样本来自人、猪还是狗。
k-mer分类结果适合做快速筛查。Kraken类工具、Centrifuge类工具或基于Mash的相似性分▶️析可以先给出候选物种🌅,但分类数据库的构建版本会改变结果,因此必须记录数据库来源、版本和过滤参数,并用独立比对复核。
DNA样本数据和参考基因组承担的任务不同,下载前必须先区分两者。样本数据通常来自测序项目,常见文件是FASTQ,里面保存实际测得的读段及质量值;参考基因组通常是FASTA,里面保存经过组装的染色体或 contig 序列,常配套GFF、GTF、索引和染色体长度文件。
人类测序数据还需要额外检查访问权限。公开数据可以直接下载👍,受控访问数据则必须按照数据仓库的🔥授权要求获取,不能通过改文件名、复制他人凭证或绕过权限取得个人基因组信息。
实际处理时,建议按照“明确数据类型—核对元数据—下载原始文件—校验压缩包—质量控制—物种分类—候选参考比对”的☀️顺序操作。只有一份短序列时可以先使用线粒体标记或保守基因进🌺行初筛;拥有全基因组或全外显子测序数据时,应优先采用核基因组证据,避免把污染、样本混合或错误注释误判成物种差异。
人or猪or狗的物种判定必须建立在可用读段之上,低质量碱基、接头残留和过短序列会明显增加错误匹配。短读段分析可以先检查每个样本的总读段数、平均质量、N碱基比例、接头污染和重复率💪,再决定是否进行剪切。
短序列的物种判断需要考虑片段长度和保守程度。保守基因的一小段序列可能同时匹配多个哺乳动物,只有包含足够特异位点、且正反向比对结果一致时,结论才更可信。
DNA数据下载流程的第一步是核对 accession 和样本元数据。搜索结果中可能同时出现项目编号、样本编号、实验编号、运行编号和文件编号,下载前应确认这些编号是否属于同一个样本,以及测序平台、读长、单双端类型和物种描述是否一致。
压缩格式不能替代完整性验证。文件能够解压,只说明压缩结构基本可读,不能证明文件从头到尾没有缺失;校验值一致、读段数量合✨理、配对关系正确,才构成可用的下载结果。