案例1  ChinaMAP 分析 10,588 个人的高深度全基因组序列[1]

研究材料：10588个中国人，随机选择于中国27个省份的8个民族(汉族、回族、满族、苗族、蒙古族、彝族、藏族和壮族) 。平均基线年龄为54岁，女性为64.8%。

分析内容：构建高质量的中国人群遗传变异数据、中国人群体结构分析、基因组特征比较以及变异频谱和致病性变异解析。该研究也通过全基因组关联分析探索了中国人群中2型糖尿病和肥胖遗传相关因素。在血糖相关分析中，结果验证了部分已知的2型糖尿病风险高频基因位点，包括CDKAL1、SLC30A8、SND1-PAX4、IDE-KIF11-HHEX、CDKN2A-CDKN2B、KCNQ1 等，也鉴定和发现了DENND5B、ORM1 、MAFA、PAX6、SOX4 等新位点。

重要研究成果：

图1 ChinaMAP基因变异的分布和模式

图2 体重指数与血糖的关联分析

ChinaMAP研究也通过全基因组关联分析探索了中国人群中2型糖尿病和肥胖遗传相关因素。在体重指数BMI相关分析中，研究团队发现了新的东亚人群特异性CADM2基因位点，CADM2在动物研究中已证实参与调节体重和能量稳态。而FTO等基因在欧美人群中发现的重要肥胖相关基因位点，在ChinaMAP研究结果中并不显著。从这些发现可以提示我们，对大规模中国人群特异性的基因组学的研究，对分子机制和个体化诊治的精准医学体系建立很重要。

案例2 NBT：主流高通量测序仪在人/细菌/宏基因组测序评测结果发布[2]

由生物分子资源设施协会（Association of Biomolecular Resource Facilities, ARBF）主导的ABRF NGS II期研究成果发表于Nature Biotechnology，文章题为“Performance assessment of DNA sequencing platforms in the ABRF Next-Generation Sequencing Study”。研究团队基于来自Illumina、Pacific Biosciences、Thermo Fisher Scientific、BGI、Oxford Nanopore Technologies和Genapsys的多款测序平台，在多个实验室对同一人类基因组家族、三个单独菌株和十种细菌的宏基因组混合物进行测序，并将各平台数据进行全方位、系统性比较，分析各个测序平台的性能差异和测序质量，以提供真实全面的参考证据。

数据显示，在短读长测序平台中，DNBSEQ平台提供了非常低的测序错误率。且SNP/Indel检测的灵敏度和准确度表现也非常优秀。

案例3 DNBSEQ基因组测序揭示肺鳞癌的潜在治疗靶标^[3]

Genomic sequencing and editing revealed the GRM8 signaling pathway as potential therapeutic targets of squamous cell lung cancer

肺腺癌和肺鳞癌（LUSC）是肺癌的主要病理类型，肺鳞癌占原发性肺癌的40%～51%。目前已经有多种靶向药物应用于肺腺癌，但是肺鳞癌的治疗靶点尚没有突破性进展。文章通过外显子重测序（WES），人全基因组重测序（WGS）、靶区域捕获测序（TS）和CRISPR-Cas9基因组编辑技术，利用鳞状细胞肺癌手术肿瘤和对应的源自患者的异种移植瘤（PDX）样本探索和验证肺鳞癌的潜在治疗靶标。

文章亮点：

• Illumina HiSeq X Ten平台WES测序+ DNBSEQ平台WGS测序+ DNBSEQ平台TS测序，多平台数据联合分析
• LUSC PDX模型可广泛应用于潜在治疗目标和策略的验证
• 使用CRISPR系统对PDX肿瘤细胞中的驱动基因进行功能验证

通过基因组测序和CRISPR-Cas9基因组编辑的综合分析，在手术和PDX肿瘤上整合鉴定并验证了GRM8对LUSC肿瘤的促进功能。cAMP活化剂和MEK抑制剂可显著阻断具有GRM8突变的LUSC肿瘤细胞的增殖和存活。因此，GRM8信号传导通路的组成分子可能成为携带GRM8激活突变的鳞状细胞肺癌的治疗靶标。

案例4 BGISEQ-500和HiSeq X Ten全基因组测序鉴定生殖细胞和体细胞变异^[4]

Germline and somatic variant identification using BGISEQ-500 and HiSeq X Ten whole genome sequencing

该研究使用BGISEQ-500平台对三种恶性胸膜间皮瘤及其对照的正常样本进行全基因组测序，并与Illumina HiSeq X Ten平台测序结果进行评估。两平台数据均使用相同的分析流程，分别比较生殖细胞和体细胞单核苷酸变异（SNP）、小插入或缺失（InDel）。结果表明BGISEQ-500平台通过全基因组测序来鉴定肿瘤样本的体细胞和生殖细胞突变是有潜力的可适用性的平台，这也是该平台首次公开可用的癌症基因组数据。

研究结果：

• 生殖细胞突变：     
  

结果显示，BGISEQ-500平台和HiSeq X Ten平台识别SNP的能力与SNP分型芯片（Infinium Omni2.5–8, Illumina ）是高度一致的（> 99％）。在两个测序平台中鉴定的生殖细胞SNV和indels也是高度一致（分别为86％和81.5％）。

表  SNP芯片数据分别与BGISEQ-500和HiSeq X Ten数据比较，生殖细胞突变基因型一致性的百分比

• 体细胞突变：
  

        三名患者中总共10,890个体细胞SNV，大部分体细胞SNV（72％）在两个平台中被识别，小部分为BGISEQ-500和HiSeq X Ten两平台特有的（分别为14%，14%）。

图  利用BGISEQ-500和HiSeq X Ten的数据鉴定3个胸膜间皮瘤的体细胞突变

案例5 基于DNB的测序平台可有效避免index hopping[5]

Reliable Multiplex Sequencing with Rare Index Mis-Assignment on DNB-Based NGS Platform

本研究使用三种主要的文库制备方法研究了DNB测序平台的Index hopping问题。DNBSEQ测序仪利用独特的DNA纳米球（DNB）技术，基于滚环复制（RCR）进行文库扩增，这种线性扩增可以避免常规PCR带来的错误累积。基于DNB的NGS应用仅使用单个index就实现了前所未有的0.0001％至0.0004％低样本错误分配率。此外，用水代替DNA，加入index，增加空白对照，DNB测序平台发生错误匹配的概率为36 million reads分之一，即0.0000028％。

图  不同测序技术的index hopping比例

研究结果：

• DNA纳米球技术的高index保真度
  

        DNBSEQ平台将DNB加载到规则阵列（patterned arrays）上，并利用组合引物锚定测序技术（cPAS）进行测序。 独特的DNB技术采用具有强链置换活性的Phi29聚合酶和能够进行线性扩增的RCR工艺，每个扩增循环都以原始的单链环状DNA文库为模板，保持每个拷贝子的独立性（图1a）。因此，即使出现寡核苷酸的index hopping等错误，也不会累积错误拷贝，正确的序列总是会在后面的DNA拷贝中复制，保证高的扩增保真度。

图  Index hopping在不同的测序平台产生的机制

• PCR-free文库index hopping污染率极低

        除了常规PCR文库外，文中还对PCR-free文库在DNBSEQ平台的index hopping情况进行调查，未经过任何Q30过滤的99.9998％精度再次证实了DNB可以在很大程度上降低index污染。与上面的常规PCR文库类似，污染率平均约为0.0004%。

表  PCR-free 文库index污染比率

研究意义：

1、  高的检测准确度，保证体细胞低频突变、HPV检测等基因检测的准确性；

2、  Single index避免了繁琐的non-combinatorial dual index带来的额外成本和劳动力浪费；

3、  避免大通量测序中样本数据完整性的丢失。

参考文献

[1]    Cao Y, Li L, Xu M, et al. The ChinaMAP analytics of deep whole genome sequences in 10,588 individuals[J]. Cell research, 2020, 30(9): 717-731.

[2]    Foox, J., Tighe, S.W., Nicolet, C.M. et al. Performance assessment of DNA sequencing platforms in the ABRF Next-Generation Sequencing Study. Nat Biotechnol 39, 1129–1140 （2021）.

[1]    Genomic sequencing and editing revealed the GRM8 signaling pathway as potential therapeutic targets of squamous cell lung cancer.[J]. Cancer letters, 2018.

[2]     Patch A M, Nones K, Kazakoff S H, et al. Germline and somatic variant identification using BGISEQ-500 and HiSeq X Ten whole genome sequencing.[J]. Plos One, 2018, 13(1):e0190264.

[3]     Li Q, Zhao X, Zhang W, et al. Reliable Multiplex Sequencing with Rare Index Mis-Assignment on DNB-Based NGS Platform. bioRxiv, 2018: 343137

数据展示

标准品数据展示

        测试样本选用了“瓶中基因组（Genome in a Bottle）”的人类样本NA12878，这是目前被世界上认为研究最透彻的二倍体人类基因组，并发布了高置信变异集，可作为一个重要工具来了解测序仪和检测结果的表现。

1）有效深度30X数据展示

• 高测序数据质量

        测试数据有至少96%的碱基识别准确率高达99%，至少89%的碱基识别准确率高达99.9%。

表1 测序质量 

        测序质量值可衡量碱基未正确检出的概率。测序技术，一种类phred算法[1,2]会为片段中的每个碱基分配一个质量分值，与最初因桑格测序实验而开发的算法类似。

        一个给定碱基的测序质量分值Q定义为下面的等式：

         Q = -10log10(e)

        其中，e为预计碱基检出不正确的概率。

        • Q分值较高表示出错的概率较小。

        • Q分值较低可能会导致相当大一部分的片段不可用，还可能导致假阳性的变异检出增加，以致得出不准确的结论。

        如下所示，质量分值20表示错误率为1/100，相应的检出精确度为99%。

表2 标测序质量分值与碱基检出精确度的关系

• 高比对率和覆盖度

        应用人的标准品 NA12878，分别在 DNBSEQ 和 N平台分别测了 100Gb Raw data 数据。从下表可以看出，由于 N平台 Duplicate rate 较高，需多测一定比例的数据，才有和 DNBSEQ 平台相当的有效深度。

表3 DNBSEQ 与 N平台比对数据比较

Clean bases：过滤掉接头，低质量和含N的reads后剩下的碱基数量；

Mapping rate：碱基比对率，比对到参考基因组的碱基数目除以clean data的碱基数目，如果测序样本存在污染或者与参考基因组差异较大，比对率偏低会影响后续的信息分析；

Unique rate：比对到基因组上唯一位置的base比率，一条reads在相同数量的容错时会有两个或者两个以上的位点都吻合，那么，它的比对结果不唯一。对于某些下游分析，需要去除比对多个位点的reads，只保留唯一比对的reads；

Duplicate reads：重复的 reads 所占比例，为了保证后续变异分析的准确性，会去掉duplicate reads后进行下游信息分析，相同数据量重复率越低，后续可用的数据量越多；

Mismatch rate：碱基的错配率；

Average sequencing depth：有效平均深度（不计算duplication），比对到参考基因组的碱基数目除以基因组的大小；目前行业对外承诺的30X（90G）、40X（120G）等深度只是测序量的简单换算，并不是指有效深度。

Coverage at least 1X（4X、10X、20X）：覆盖率，指测序深度达到1X、4X、10X、20X以上的全基因组占比。

• 高灵敏度和精准度

        高灵敏度（Sensitivity）和高精准度（Precision）意味着 DNBSEQ 平台检测发现变异的能力更强，并且结果中为真的突变的概率也高。

表4 DNBSEQ 与 N平台 SNP 精准度和敏感度对比

表5  DNBSEQ 与 N平台 InDel 精准度和敏感度对比

Sensitivity：灵敏度，又叫真阳性率（TPR），计算公式：灵敏度=真阳性/（真阳性+假阴性）。是指实际为阳性的样本中，判断为阳性的比例。例如，真正突变中，被判断为有突变的比例，它反映筛检发现变异的能力，灵敏度越高，假阴性越低；

Precision：精准度，也叫阳性预测值（PPV），计算公式：精准度=真阳性/（真阳性+假阳性），指筛检试验检出的全部阳性变异中，真正“变异”的例数（真阳性）所占的比例，反映筛检变异结果阳性中为真的突变的可能性，精准度越高，假阳性越低。

        *上述分析结果由华大信息分析流程所得，本结果不代表交付指标，最终解释权归深圳华大基因股份有限公司所有。

参考文献

[1] Li Q, Zhao X, Zhang W, et al. Reliable Multiplex Sequencing with Rare Index Mis-Assignment on DNB-Based NGS Platform. bioRxiv, 2018: 343137

[2] Illumina. Effects of Index Misassignment on Multiplexing and Downstream Analysis (white paper). 4 (2017). doi:10.1101/125724 


表1 DNA样本送样建议

注意事项：

1) 务必附上凝胶电泳、NanoDropTM、Qubit®、Agilent Bioanalyzer等其中至少一种的检测结果，电泳图需标明所用marker的条带大小。

2) 样品质量以BGI的质检结论为准，望合作伙伴理解，检测结果可能会由于检测地点，仪器设备和操作者等不同造成固有差异。因质检有一定的消耗量，合作伙伴寄送的样本量必须高于各产品样品标准至少50ng以上。强烈建议根据2倍以上标准制备样品，否则很可能会导致大量样本质检未能达标，延误项目进展。

3) BGI原则上只接收1.5mL/2.0mL EP管，要求每管样品体积在15-100μL之间（推荐30μL），根据实验要求，如果样品体积小于15μL，BGI可能会在检测之前稀释原始样品。

表2 组织样本送样建议

注意事项：

1) 2 mL 螺纹旋盖保存管

2) 组织样品保存方法选择:首选液氮速冻;没有液氮条件的，可直接放入-80°C冰箱冻存; 环境条件限制的，可使用商业核酸保护液保存，并严格按相应试剂说明操作。

3) 长年保存的组织:保存时间超过一年的组织不建议送样。 

Q1：DNBSEQ 人全基因组重测序的数据格式是否与 Illumina 平台的一致？

           是一致的。所以信息分析流程都一样。

Q2：想对我们 DNBSEQ 产出的数据先分析确定一下格式和数据质量，现在是否有测试过的数据可以提供？

           华大智造官网，测试数据下载：https://www.mgi-tech.com/Resource/demo/ping_id/3 

Q3：如何实现基因组变异可视化？

           基因组可视化软件 IGV (Integrative Genomics Viewer) 是高性能的基因组数据可视化工具，能够帮助使用者同时合并分析不同类型的基因组数据，并能灵活放大基因组上的某个特定区域。IGV 软件免费下载地址: http://www.broadinstitute.org/igv. IGV 可查看 SAM / BAM 比对文件和 VCF 变异检测文件，下图显示的是 IGV 可视化窗口。

Q4：如何寻找候选变异？

           可以使用“人基因组变异交互分析系统”进行后续变异位点筛选，寻找候选变异位点时，可利用变异注释结果，关注非同义突变、剪接突变、移码突变。

           1)去除千人基因组数据库中 MAF >=1% 的变异

           2)去除 NHLBI-ESP6500 European American 群体数据库中 MAF >=1% 的 变异 

           3)去除 NHLBI-ESP6500 African American 群数据库中 MAF >=1%的变异

           4) 推 测 变 异 的 致 病 性 。 利 用 SIFT/PolyPhen2/Mutation assessor/Condel/FATHMM 进行打分，预测某个变异和氨基酸置换是否影响蛋白 功 能 。

            如 果 score<=0.05 或 PolyPhen2>=0.909 或 MA score>=1.9 或 Condel = deleterious 或 FATHMM=deleterious，就推测该变异可能是有害变异。

Q5：一般用什么方法来验证 call SNP 准确率?

           华大炎黄计划是用 Sanger 测序的方法和芯片分型两种方法来验证 SNP 的准确性的， 因为 Sanger 测序被认为是测序中的“金标准”。

Q6：GWAS后期验证一般需要用什么方案？客户已经做过GWAS，找出的位点有一些在内含子部分，想要进行后期验证。

1、在更大的群体上针对候选位点质谱验证；

2、动物模型上验证。

Q7：FFPE样本可以承诺什么指标？

           FFPE样本由于保存年份和降解程度的差异，不同样本的情况，测序质量、duplication比率、覆盖度等非常不同，属于风险建库，只承诺Raw data的数据量，其他都不能承诺。

Q8：用人唾液测全基因组，效果如何？

           可以做，因为会含有口腔里的微生物，存在污染，一些商业样本的结果显示，比对率、覆盖度等均比血液样本略低。

Q9：突变位点为有效位点时使用的 depth 阈值是多少?

           GATK在call变异时SNP和InDel均要求depth大于等于4 。

Q10：数据中的 Duplication 指什么?如何定义?有何影响？

           生物学意义:由同一个序列经过PCR扩增而产生的不同的reads，处理的时候去掉这些不同的reads只留一条。

           生物信息学上意义:跟比对软件有关; 严格定义:起始和终止坐标都一样，mismatch 的位置和类型都一样的不同的reads，算duplication; 宽松定义:起始坐标一样的不同的reads，定义为 duplication。位点一样，碱基也一样的就是确定的 duplication; 而那些去完 adaptor 之后的信息(起始终止坐标，mismatch 位置，类型等)一样的不同 reads 不一定就是duplication，因为没有算上去掉的那一部分。与测序深度可能有关系，理论上，测序深度越高，得到的 duplicate 的 reads 也会越多，因为这些 reads 也会比对到基因组上，在 call SNP 的时候，会对局部的覆盖深度有影响， 甚至于影响到 SNP 的进一步过滤，所以对于变异的检测也是有干扰的。对于InDel 和 SV 也是类似的。