亲权鉴定:从STR分型到亲权指数(PI)的计算逻辑

栏目:法医百科 0 人浏览

先看一个真实工作场景:一家司法鉴定机构接到委托,离婚诉讼里男方要确认自己和孩子有没有血缘关系。实验室采了三份血样,跑了二十多个STR基因座,出具的意见书里写着一行字:累计亲权指数3.2×107,支持被检男子是孩子的生物学父亲。这个数字怎么来的,凭什么一个数就能定"是不是亲生的",这篇文章把计算逻辑讲透。

一、两个先要弄明白的定义:亲权鉴定和STR

国家标准GB/T 37223-2018《亲权鉴定技术规范》对亲权鉴定的定义很朴素:通过对人类遗传标记的检测,根据遗传规律分析,对个体之间血缘关系的鉴定。拆开看有三个要件:检测的对象是遗传标记,分析的依据是遗传规律,回答的是血缘关系有没有。

先把标准的沿革交代一句。这个国家标准2018年12月28日发布,2019年4月1日实施,取代了司法部技术规范SF/Z JD0105001-2016。司法部以司发通〔2019〕33号通知,自2019年4月1日起废止那批技术规范,对应表里写得很清楚:旧规范SF/Z JD0105001-2016,新标准GB/T 37223-2018。眼下鉴定机构执行的是国家标准,下文说的"国家标准"都指GB/T 37223-2018。

短串联重复序列(short tandem repeat, STR)是今天用得最广的那类遗传标记。它的核心序列由2到6个碱基组成,串联重复的次数不同,就形成长度不等的等位基因。拿D13S317基因座说,重复7次是等位基因7,重复11次是等位基因11。每个人在每个常染色体STR基因座上有两个等位基因,一个来自父亲,一个来自母亲。这就是孟德尔分离定律在分子水平上的样子。

国家标准对选用的STR基因座还有一条推荐性要求,条款用的是"宜":串联重复单位宜为四核苷酸或五核苷酸。四核苷酸重复电泳时"影子带"少,分型干净。实验室按两种模式受理。三联体是:被检男子、生母、孩子。二联体是:被检男子和孩子,生母不参加。

用STR做亲权鉴定,第一层逻辑是排除。孩子身上找不到来源的等位基因,只能来自生父,这个等位基因叫生父基因。被检男子提供不了生父基因,这个位点就排除他。一个遗传标记系统对非父的排除能力叫排除概率,多个基因座合起来是累计非父排除率(CPE)。国家标准要求,亲权鉴定用的STR组合,CPE不应小于0.9999。意思是:如果他不是亲爹,这套系统有99.99%以上的机会把他揪出来。没别的。

二、亲权指数的定义:一道比值,两种假设

排除只回答"是不是",回答不了"有多像"。司法上需要的不止是排除,还要在支持时给出证据强度。亲权指数(paternity index, PI)就是干这个的。

国家标准给亲权指数的定义是:亲权关系鉴定中判断遗传证据强度的指标。它是两个条件概率的似然比率,写成公式:PI=X/Y。

X:假设被检男子是孩子的生物学父亲,观察到当前分型结果的概率。Y:假设一个与案件无关的随机男子是孩子的生物学父亲,观察到当前分型结果的概率。

PI大于1,说明这组分型在"被检男子是生父"的假设下更可能出现,数值越大,对父权的支持越强。PI小于1,结果偏向随机男子假设。PI等于0,这个位点直接排除。

这个思路有年头了。1938年,瑞典学者Essen-Möller提出用似然比的框架处理父权问题。他还给出把比值折成概率的W值,W=X/(X+Y)。Gürtler主张直接报告比值X/Y,这就是今天PI的来历。Hummel在1984年的综述里把这段历史理得很清楚。W值后来演变成亲权概率(也叫父权相对机会,RCP),作为PI的另一种表达保留下来。

二联体和三联体的PI公式不一样。二联体缺了生母,孩子另一个等位基因的来源得按群体频率去估计,所以二联体的PI普遍比三联体小。这也是鉴定实践里条件允许就尽量做三联体的原因,证据强度更扎实。

三、计算逻辑:从一个基因座到累计亲权指数

三联体的公式表

三联体情形下,生母的基因型能帮我们锁定孩子哪个等位基因是生父基因。用p、q、r表示等位基因P、Q、R的群体频率,常见情形的公式节选如下:

生母基因型孩子基因型被检父基因型PI计算公式
PPPPPP1/p
PPPPPQ1/(2p)
PPPQQQ1/q
PPPQQR1/(2q)
PQPQPP1/(p+q)
PQPQPQ1/(p+q)
PQPQPR1/[2(p+q)]

看第一行:生母PP、孩子PP、被检男子PP。被检男子若为生父,必然传递P,概率1;随机男子传递P的概率是p。PI=1/p。p越小,这份证据越值钱。

再看倒数第三行:生母PQ、孩子PQ、被检男子PP。孩子从母亲拿P还是拿Q定不下来,生父基因可能是P也可能是Q。被检男子是生父时传递P,概率1;随机男子传递P或Q的概率是p+q。PI=1/(p+q)。

二联体的公式表

二联体没有生母参与,公式里要多乘一项"母亲随机贡献",数值自然变小:

孩子基因型被检父基因型PI计算公式
PPPP1/p
PPPQ1/(2p)
PQPP1/(2p)
PQPQ(p+q)/(4pq)
PQPR1/(4p)

逐位点相乘

一个基因座给出一个PI。不同STR基因座在染色体上相距足够远,遗传上相互独立,各个基因座的PI可以直接相乘,得到累计亲权指数(combined paternity index, CPI):CPI=PI1×PI2×…×PIn。

美国国家司法研究所(NIJ)的教材里有个4基因座的小例子,可以完整走一遍:

基因座生母孩子被检父等位基因频率PI
D3S1358161615,16p16=0.23151/(2p16)=2.16
vWA17,1817,1818,20p17=0.2627,p18=0.22191/[2(p17+p18)]=1.03
FGA25,2723,2521,23p23=0.15811/(2p23)=3.16
D18S5114,181818,20p18=0.09181/(2p18)=5.45

四个基因座的PI分别是2.16、1.03、3.16、5.45,CPI=2.16×1.03×3.16×5.45=38.4。意思是:观察到这组分型,在"被检男子是生父"假设下比在"随机男子是生父"假设下可能性高38倍。

4个位点就有38倍。中国实验室常规做20个以上位点,文献报告的案例里CPI从104到1011都有。CPI还能换算成亲权概率:W=CPI/(CPI+1)。在前概率取0.5的前提下,CPI超过10000,对应的W就高于0.9999。这就是"亲权概率大于99.99%"这句话的出处。

四、不符合孟德尔遗传:突变的机制和算法

真实检案里常有这种分型表现:孩子某个等位基因既不像父亲也不像母亲,图谱上冒出一个"陌生等位基因"。这种不符合孟德尔遗传的现象,绝大多数是STR突变。

突变的发生机制,学界接受度最高的是DNA复制滑动(也叫滑链错配):DNA合成时新链和模板链错位,重复单位多一个或少一个。所以STR突变以一步突变为主,占90%以上。Brinkmann团队1998年发表的数据显示,突变率和核心序列的结构、重复长度有关,重复越长突变率越高。法医界给常染色体STR的平均突变率常取0.2%这个量级。

国内数据印证了这条规律。内蒙古1360例鉴定里,23个基因座平均突变率0.15%,D12S391最高到0.38%。广东10271例家系里,母源突变率在0.003%到0.05%之间,95.4%是一步突变。父源突变明显多于母源,比例大约3:1到4:1。原因不复杂:精子经历的分裂次数比卵子多得多,复制出错的机会跟着多。

突变位点不能丢掉不算。国家标准明确:任何情况下都不能为了获得较高的累计亲权指数,把不符合遗传规律的遗传标记删除。突变位点用突变模型算PI。以μ表示这个基因座的平均突变率,p为突变等位基因频率:三联体一步突变PI=μ/(4p),两步突变PI=μ/(40p);二联体一步突变PI=μ/(8p),两步突变PI=μ/(80p)。男性突变率μ取0.002,女性取0.0005。

分子是突变发生的概率,分母是随机男子碰巧提供这个等位基因的概率。逻辑没变:证据还是两个假设的比值。说实话,突变位点怎么算,是亲权鉴定里最容易翻车的地方。算错了,一个本来支持亲权的案子可能掉进"不排除也不支持"的灰色区间。

五、判定标准与法医学意义

国家标准定的判定标准,一张表说清楚:

指标数值含义
累计非父排除率CPE不应小于0.9999检测体系的前提条件
累计亲权指数CPI大于10000支持亲权关系,对应亲权概率高于99.99%
累计亲权指数CPI小于0.0001排除亲权关系
累计亲权指数CPI0.0001至10000之间(含两端)增加检测的遗传标记,直到结论明确

两个阈值放一起看很有意思。支持线在104之上,排除线在10-4之下,中间隔着8个数量级的灰色地带。落在灰色地带的案子,国家标准要求增加检测的遗传标记,把结论推向两端,不许模棱两可地出具意见。

还有几条红线值得记住:任何情况下都不能仅凭一个遗传标记不符合遗传规律就作出排除意见;鉴定意见以不考虑双胞胎和近亲为前提,近亲案件要用另外的公式;对接受过外周血干细胞移植的当事人不能采血,要采口腔拭子或毛发;Y-STR、X-STR和线粒体DNA只能做补充检验和描述性分析,不推荐跟常染色体STR的指数联合计算。

从法医学意义上讲,PI和CPI把"亲子关系"这个高度情感化、容易各执一词的问题,压缩成了一个可复核、可质证的数字。法庭上双方争议的不再是"像不像",而是分型对不对、频率用哪个人群、公式用得对不对。这对鉴定意见的质量控制是好事。

它不是万能钥匙。突变、无效等位基因、近亲干扰、嵌合体,都可能让数字偏离事实。近年学界还在发展基于SNP高通量测序的无创产前亲权鉴定,以及全同胞、半同胞、祖孙等复杂亲缘关系指数的通用算法。方向没变:把更复杂的亲缘问题,装进同样的似然比框架里。

这件事的底线是:数字必须来自标准,不能来自感觉。

参考文献

  1. 国家市场监督管理总局、国家标准化管理委员会. GB/T 37223-2018《亲权鉴定技术规范》(现行有效,2018-12-28发布,2019-04-01实施,主管部门:司法部). 国家标准全文公开系统
  2. 司法部. 关于废止《印章印文鉴定规范》等15项司法鉴定技术规范的通知(司发通〔2019〕33号,2019年3月8日):自2019年4月1日起废止SF/Z JD0105001-2016《亲权鉴定技术规范》,对应替代标准为GB/T 37223-2018. 通知全文及废止对应表
  3. 侯一平. 法医物证学. 第4版. 北京:人民卫生出版社,2016. ISBN 978-7-117-22236-5. 人卫智慧服务商城
  4. Hummel K. On the theory and practice of Essen-Möller's W value and Gürtler's paternity index (PI). Forensic Science International, 1984, 25: 1-17. DOI链接
  5. Gjertson DW, Brenner CH, Baur MP, et al. ISFG: Recommendations on biostatistics in paternity testing. Forensic Science International: Genetics, 2007, 1(3-4): 223-231. DOI链接
  6. Brinkmann B, Klintschar M, Neuhuber F, et al. Mutation rate in human microsatellites: influence of the structure and length of the tandem repeat. American Journal of Human Genetics, 1998, 62(6): 1408-1415. DOI链接
  7. 兰菲菲, 梁杰, 胡听听, 等. 常染色体STR基因座母源突变的观察分析与亲权指数计算. 中国生育健康杂志, 2022, 33(4). 期刊官网全文
  8. 郭振民, 武海鹏, 霍志鑫, 等. 亲子鉴定中23个常染色体STR基因座突变分析. 中国医科大学学报, 2024, 53(4): 363-366. 期刊官网全文
  9. National Institute of Justice. Population Genetics and Statistics for Forensic Analysts: Paternity Calculations. NIJ官网