易 娜 ，黄仁武 ，蒋佳梅 ，孟存溦 ，叶 好 ，李湘萍 ，刘林林 ，胡利平

（1）昆明医科大学法医学院法医学系，云南 昆明 650500;2）昆明医科大学司法鉴定中心，云南 昆明 650500）

近年来，随着法医学鉴定技术的发展，复杂亲缘关系鉴定需求日益增加，全同胞关系作为复杂亲缘关系的一种，指有相同生物学父亲和生物学母亲的多个子代个体之间的关系[1]，目前在已在司法实践的多个领域发挥着重要的作用[2]。因此，根据实际需要，司法部司法鉴定管理局于2014 年发布了《生物学全同胞关系鉴定实施规范》（SF/Z JD0105002—2014，后文简称“旧规范”），规范了基于状态一致性评分（identity by state（IBS）score）的生物学全同胞关系鉴定的技术方法，但由于旧规范中过高的系统效能和对不同群体之间的群体差异考虑的欠缺，导致在实际应用过程中存在一定的误判风险。随后，司法部于2021 年发布并实施了《生物学全同胞关系鉴定技术规范》（SF/T 0117—2021，后文简称“新规范”）以代替旧规范。新规范相较于旧规范更新了累计状态一致性评分（CIBS）的判定阈值并降低了对应检测系统的系统效能，增加了基于等位基因频率计算的全同胞指数（full sibling index，FSI）和累积全同胞指数（CFSI）的判定阈值及相应的系统效能[3]。本文基于新旧规范中相同的19 个必检基因座检测系统探讨了2 份规范中3 种不同全同胞鉴定阈值划分方法在系统效能上的差异以及新规范中两种方法在云南人群生物学全同胞关系鉴定中的适用性。

1 材料与方法

1.1 样本来源

在获得知情同意的前提下，收集1 468 对来自昆明医科大学司法鉴定中心日常案件积累的云南人群样本（口腔拭子、血样），其中全同胞有448 对，无关个体有1 020 对。全同胞对为已通过旧规范中IBS 法认定（包括增加20 个基因座后认定为全同胞的个体对）的兄弟姐妹组合，排除了同卵双生子和STR 突变的家系。无关个体对为实验室随机抽取的无亲缘关系的个体两两随机组合。

1.2 STR 基因座选择及分型检测

案件鉴定用Chelex-100 法提取血样 DNA，采用Promega 公司PowerPlex®21 System 试剂盒进行扩增，扩增产物经 ABI3130XL 型基因分析仪电泳分离，应用ABI GeneMapperID-1.5 软件分析。每个样本均得到D3S1358、D1S1656、D6S1043、D13S317、Penta E、D16S539、D18S51、D2S1338、CSF1PO、PentaD、TH01、vWA、D21S11、D7S820、D5S818、TPOX、D8S1179、D12S391、D19S433、FGA 和 Amelogenin 共21 个STR 基因座的分型结果，包含全同胞鉴定规范中的19 个必检基因座。

1.3 计算累计状态一致性评分

两名有争议个体在同一STR 基因座上出现相同等位基因的个数就是状态一致性评分。当采用包含多个相互独立的常染色体遗传标记分型系统对两名有争议个体进行检测时，各个遗传标记上IBS 之和即称为累计状态一致性评分[3]。利用全同胞鉴定规范中状态一致性评分方法，分别计算每对全同胞和无关个体之间19 个必检基因座的CIBS 值，并对数据进行统计整理。

1.4 计算累积全同胞指数

对于每一个STR 基因座而言，两名有争议个体之间存在全同胞关系时其基因型出现的机率与两名有争议个体之间为无关个体时其基因型出现的机率之比值即为全同胞指数。当采用包含多个相互独立的常染色体遗传标记分型系统对两名有争议个体进行检测时，各个遗传标记上FSI 的乘积即称为常染色体STR 基因座累积全同胞关系指数[3]。根据云南汉族人群中 19 个常染色体STR基因座的等位基因频率资料[4]和新规范中全同胞指数和累积全同胞关系指数的算法[5]，分别计算全同胞组和无关个体组的配对个体之间19 个必检基因座的CFSI 值，并对数据进行统计整理。

1.5 3 种方法的阈值划分

基于新旧规范，用IBS 法和FSI 法分别计算1 468对样本的全同胞关系，并对计算结果按照以下三种方法进行划分：方法一以旧规范中 IBS ≥ 22 倾向于认为两名被鉴定人为全同胞，IBS ≤ 13 倾向于认为两名被鉴定人为无关个体和13 ＜ IBS ＜ 22无法给出倾向性意见为标准进行阈值划分；方法二以新规范中IBS＞22 倾向于认为两名有争议个体为全同胞关系，IBS ＜ 12 倾向于认为两名有争议个体为无关个体对和12 ≤ IBS ≤ 22 无法给出倾向性意见为标准进行阈值划分；方法三以新规范中FSI＞1.0×104倾向于认为两名有争议个体为全同胞关系，FSI ＜ 1.0×10-4倾向于认为两名有争议个体为无关个体对和1.0×10-4≤ FSI ≤ 1.0×104无法给出倾向性意见为标准进行阈值划分。

1.6 统计学处理

采用SPSS 25.0 软件对全同胞组和无关个体组的lgCFSI 值和CIBS 值的分布进行正态检验后，统计不同组别lgCFSI 值和CIBS 值的算数均数和标准差，并对其进行独立样本t检验，对计数资料采用χ2检验，均以P＜ 0.05 为差异有统计学意义。

2 结果

2.1 CIBS 值和CFSI 值的分布

根据448 对全同胞和1 020 对无关个体的CIBS 值和对CFSI 值经对数转换后的lgCFSI 值进行正态检验，结果表明全同胞组和无关个体组的CIBS 值及lgFSI 值均符合正态分布（图1，图2）。

图1 全同胞组和无关个体组的CIBS 值拟合所得正态分布曲线Fig.1 The normal distribution curve obtained by fitting the CIBS of the full sibling group and the unrelated individual group

图2 全同胞组和无关个体组的lgCFSI 值拟合所得正态分布曲线Fig.2 The normal distribution curve obtained by fitting the lgCFSI of the full sib group and the unrelated individual group

2.2 无关个体组与全同胞组的差异

统计并描述全同胞组和无关个体组CIBS 值和lgCFSI 值的算数均数、标准差、t值（表1），结果显示无关个体组与全同胞组在CIBS 值与lgCFSI值的差异均有统计学意义，P＜ 0.001。

表1 全同胞组和无关个体组CIBS 值和lgCFSI 值的 和独立样本t 检验结果（）Tab.1 The and independent sample t test results of the CIBS values and lgCFSI values between the full sibling group and the unrelated individual group（）

表1 全同胞组和无关个体组CIBS 值和lgCFSI 值的 和独立样本t 检验结果（）Tab.1 The and independent sample t test results of the CIBS values and lgCFSI values between the full sibling group and the unrelated individual group（）

2.3 3 种方法的统计结果及差异

分别按照3 种方法对全同胞组和无关个体组的计算结果进行统计整理和χ2检验（表2），显示方法一有1 087 对样本可以得到倾向性意见，有381 对样本无法得出倾向性意见，实际检出率为74.05%；方法二有750 对样本可以得到倾向性意见，有718 对样本无法得出倾向性意见，实际检出率为51.09%；方法三有972 对样本可以得到倾向性意见，有496 对样本无法得出倾向性意见，实际检出率为66.21%。3 种方法在19 个常染色体STR 基因座检测水平上的全同胞组检测准确率和无关个体组检测准确率均为100%。3 种方法的χ2检验结果显示在全同胞组和无关个体组中，3 种方法的检出率不完全相同，差异有统计学意义。对2 个组别中3 种方法的检出率进行方法二与方法一和方法三的比较（表3），通过Bonferroni方法进行校正，结果显示在2 个组别中方法二与方法一和方法三的差异均有统计学意义。

表2 3 种方法的统计结果（n）Tab.2 Statistical results of the three methods（n）

表3 3 种方法两两比较Tab.3 Pairwise comparison of three methods

3 讨论

对于同样的19 个常染色体STR 基因座检测系统，以上不同全同胞鉴定阈值划分方法的系统效能在两份规范中分别为0.750 0，0.556 5 和0.662 5；在云南人群中的系统效能分别为0.740 5，0.510 9 和0.662 1。比较新旧规范中IBS 法的应用，因CIBS 阈值划分标准不同，两种方法的系统效能相差约0.193 5，在云南人群中2 种方法的系统效能相差约0.226 9。在新规范中减少的这部分系统效能是由于更新的规范中将个体对的得分为“CIBS=12”、“CIBS=13”、“CIBS=22”时判定为“无法给出倾向性意见”，同时靠近阈值附近的个体对数也明显较多，从而导致两种方法的系统效能相差较大。新规范降低了IBS 法判定全同胞关系的系统效能，降低了误判率，使得全同胞鉴定时给出肯定意见的准确性有所提高[6]。

FSI 法对于全同胞关系指数的计算由ITO 法提供[7]，用以判定两名有争议个体是否为全同胞关系的方法。该方法的计算过程较IBS 法稍显复杂，但FSI 法考虑到了人群群体遗传学差异，基于不同人群等位基因频率在某一检测系统下进行计算，因此能有效避免全同胞关系鉴定中的误判，增加结果的可信度[8]。本研究基于新规范以及文献报道的云南人群的等位基因频率，在19 个常染色体STR 基因座检测系统下，通过FSI 法计算的系统效能均高于IBS 法。同时，比较全同胞组和无关个体组的正态分布曲线，lgCFSI 值的重叠范围和比例明显小于CIBS 值的重叠范围和比例，也说明FSI 法能对全同胞和无关个体进行更有效的区分。

本研究考虑到实际鉴定中，不同人群同一基因座上等位基因频率分布的差异会使得同一基因座的个体识别能力有所不同[9]，从而导致同一检测系统在不同人群中的系统效能有一定差异的情况。应用新规范中IBS 法和FSI 法在19 个常染色体STR 检测系统上对1 468 对云南人群样本进行了全同胞关系鉴定，计算所得在云南人群中两种方法的系统效能分别为0.510 9 和0.662 1 与新规范中同一检测系统下对应方法的系统效能0.556 5和0.662 5 相近，说明新规范在云南人群中有很高的应用价值。对于云南地区少数民族聚集，不同的少数民族间或有部分等位基因频率存在较大差异的情况[10]，新规范规定，当某个基因座的个体识别能力与规范中存在较大差异时，可以在准确性不低于99.99%的前提下重新计算并制定全同胞关系鉴定的CIBS 阈值，这也能有效解决新规范在各地区人群中实际应用的局限。

对于无法给出倾向性意见可补充检测基因座，以增加能给出倾向性意见的可能性[11]。新规范中也给出了基于19 个常染色体STR 基因座依次补充36 个与其不存在连锁不平衡的常染色体基因座后相应的检测系统效能，对于IBS 法和FSI 法其最高系统效能分别可达0.998 6 和0.996 5[3]，能有效区分全同胞和无关个体。对于特殊的鉴定案件，还可补充检验X 染色体、Y 染色体或线粒体DNA遗传标记的情况[12]，但应该根据其遗传规律采用文字描述的方式进行分析说明。

本研究发现，基于19 个常染色体STR 基因座检测系统，新规范在云南人群全同胞关系鉴定中具有较高的法医学应用价值；相较于旧规范，新规范中的IBS 法和FSI 法能有效降低全同胞关系鉴定的错判风险。由于IBS 法和FSI 法在计算原理上的差异，本文建议实际应用中将IBS 法和FSI 法相结合，对结果进行更为科学和准确的判断。需要注意的是，无论新规范还是旧规范均只适用于同父同母的全同胞关系鉴定，对于半同胞关系的判别并不适用。实际案件中如遇到不确定为全同胞的情况，可通过ITO 法对全同胞和半同胞关系进行判别[13-14]，从而得出相应结论。