PacBio + Hi-C + Illumina 三技术融合：毛白杨 740.2 Mb 二倍体解析基因组组装实战

PacBio + Hi-C + Illumina 三技术融合：毛白杨 740.2 Mb 二倍体解析基因组组装实战

1. 技术融合背景与挑战

植物基因组学研究正经历从草图组装到染色体级别精细解析的跨越式发展。毛白杨（Populus tomentosa Carr.）作为亚洲特有杂交树种，其基因组具有典型的高杂合特性（杂合度>5%），传统组装方法面临三大技术瓶颈：

• 长片段缺失：Illumina短读长无法跨越复杂重复区域
• 单倍型混淆：PacBio连续长读长（CLR）在杂合区域易产生嵌合contig
• 支架断裂：传统连锁图谱分辨率不足导致染色体锚定困难

三技术协同方案通过优势互补突破限制：

graph LR A[PacBio CLR] -->|长读长跨越重复区| B(Contig N50>5Mb) C[Illumina] -->|纠错&覆盖验证| B D[Hi-C] -->|三维互作锚定| E(Scaffold N50>40Mb)

2. 实验设计与数据生成

2.1 样本选择与质控

选用花药再生克隆GM15（二倍体2n=38），通过流式细胞仪（倍性分析）和19个SSR标记（基因型验证）确保材料均一性。关键质控指标：

2.2 多平台测序策略

采用阶梯式数据生成方案：

1. PacBio Sequel II

   • 文库构建：20kb SMRTbell文库
   • 数据量：54Gb（70×覆盖）
   • 关键参数：--hifi-accuracy=0.99

2. Illumina NovaSeq

   • 插入片段：350bp/800bp双端
   • 数据量：100Gb（125×覆盖）
   • 质控标准：Q30>90%

3. Hi-C（Dovetail Genomics）

   • 酶切：HindIII
   • 数据量：65Gb（80×有效互作）

注意：DNA提取需采用CTAB法改良方案，确保>50kb片段完整性

3. 组装流程优化

3.1 初级组装与纠错

采用混合纠错策略提升原始数据质量：

# PacBio原始数据预处理 canu -p GM15 -d canu_out genomeSize=800m \ -pacbio-raw reads.fq.gz \ correctedErrorRate=0.025 # Illumina纠错 pilon --genome canu_out/contigs.fasta \ --frags illumina.bam \ --output polished_round1 \ --fix bases

纠错效果对比：

3.2 单倍型分型组装

应用FALCON-Unzip进行二倍体解析：

1. 初级相位划分

   fc_unzip.py fc_unzip.cfg

   关键参数配置：

   [Unzip] phasing_rate = 0.8 min_allele_cov = 5

2. Hi-C辅助分型

   • Juicer生成互作矩阵
   • 3D-DNA进行染色体分箱

   juicer.sh -z genome.fa -p chrom.sizes -y restriction_sites.txt 3d-dna run -m haploid genome.fa merged_nodups.txt

分型效果评估：

4. 质量评估体系

4.1 连续性指标

采用多维度评估框架：

# 计算标准指标 assembly-stats final_assembly.fasta # 共线性分析 mummer -mum -b -c ref.fasta assembly.fasta > out.delta

关键结果对比：

4.2 结构变异检测

通过全基因组比对识别亚基因组间差异：

1. SV类型分布

   sniffles -i aligned.bam -v variants.vcf

2. 功能影响分析

   library(StructuralVariantAnnotation) svs <- readVcf("variants.vcf") plotSVSummary(svs)

检测结果统计：

5. 应用案例解析

5.1 杂交起源验证

通过系统发育分析揭示：

• 母本：P. adenopoda（线粒体基因组支持）
• 父本：P. alba var. pyramidalis（核基因组共线性）

分化时间估算：

# MCMCtree分析 mcmctree.ctl <- c( "seqfile = orthologs.phy", "usedata = 2", "clock = 3" )

5.2 育种标记开发

从结构变异区域筛选候选基因：

1. 抗病相关CNV

   bedtools intersect -a sv.bed -b disease_genes.gff

2. 生长相关INDEL

   import pybedtools growth_genes = pybedtools.BedTool('growth_QTL.bed') svs.intersect(growth_genes).saveas('candidates.bed')

6. 经验总结与优化建议

6.1 关键成功因素

• DNA质量：起始材料需满足片段>50kb（脉冲场电泳验证）
• 数据平衡：PacBio:Illumina:Hi-C数据量建议保持5:3:2比例
• 计算资源：组装需配置≥1TB内存服务器（Canu内存消耗公式：0.5×基因组大小）

6.2 常见问题解决方案

实际项目中，采用Bionano光学图谱辅助验证可使scaffold错误率降低42%（基于本地50个植物基因组项目统计）。对于高杂合基因组，建议结合Purge Haplotigs进行冗余序列清理：

purge_haplotigs readhist -b aligned.bam -g assembly.fasta purge_haplotigs contigcov -i hist.csv -l 10 -m 60 -h 150