删除或更新信息，请邮件至freekaoyan#163.com(#换成@)

基于有限突变模型和大规模数据的19个常染色体STR的实际突变率研究

本站小编 Free考研考试/2022-01-01
刘志勇2, 任贺3, 陈冲4, 张京晶5, 张晓梦1, 石妍4, 石林玉1, 陈滢4, 程凤1, 贾莉4, 陈曼6, 范庆炜7, 张家榕1, 李万婷1, 王萌春1, 任子林8, 刘雅诚4, 倪铭8, 孙宏钰2, 严江伟,1 1. 山西医科大学法医学院,太原 030001
2. 中山大学中山医学院法医学系,广州 510080
3. 北京警察学院,北京 102202
4. 北京通达首诚司法鉴定所,北京 100192
5. 北京华彦科技有限公司司法鉴定所,北京 100192
6. 南方医科大学法医学院,广州 510515
7. 川北医学院法医学系,南充 637000
8. 北京辐射医学研究所,北京 100850

Actual mutational research of 19 autosomal STRs based on restricted mutation model and big data

Zhiyong Liu2, He Ren3, Chong Chen4, Jingjing Zhang5, Xiaomeng Zhang1, Yan Shi4, Linyu Shi1, Ying Chen4, Feng Cheng1, Li Jia4, Man Chen6, Qingwei Fan7, Jiarong Zhang1, Wanting Li1, Mengchun Wang1, Zilin Ren8, Yacheng Liu4, Ming Ni8, Hongyu Sun2, Jiangwei Yan,1 1. School of Forensic Medicine, Shanxi Medical University, Taiyuan 030001, China
2. Faculty of Forensic Medicine, Zhongshan School of Medicine, Sun Yat-sen University, Guangzhou 510080, China
3. Beijing Police College, Beijing 102202, China
4. Beijing Tongda Shoucheng Institute of Forensic Science, Beijing 100192, China
5. Beijing Huayan Judicial Authentication Institute, Beijing 100192, China
6. School of Forensic Medicine, Southern Medical University, Guangzhou 510515, China
7. Faculty of Forensic Medicine, North Sichuan Medical College, Nanchong 637000, China
8. Beijing Institute of Radiation Medicine, Beijing 100850, China

通讯作者: 严江伟，博士，教授，研究方向：法医基因组学。E-mail:yanjw@sxmu.edu.cn

编委: 朱波峰
收稿日期:2021-06-2修回日期:2021-08-18
基金资助: 国家自然基金重点项目资助编号.82030058


Received:2021-06-2Revised:2021-08-18
Fund supported: Supported by the National Natural Science Foundation of China No.82030058

作者简介 About authors
刘志勇,博士研究生,研究方向：法医遗传学。E-mail:liuzhy255@mail2.sysu.edu.cn;

任贺,硕士,副教授,研究方向：法医遗传学。E-mail:snoopy_fr@sohu.com;刘志勇和任贺并列第一作者。





















摘要
短串联重复序列(short tandem repeat, STR)已广泛用于法医学亲子鉴定和个体识别中,但STR的突变可能会影响其结果的解释。在大多数类似研究中,由于忽略“隐性”突变现象,STR的突变率被低估。鉴于此,为获得更加准确的STR实际突变率,本研究使用Slooten与Ricciardi提出的有限突变模型和大规模数据,对28,313例(78,739个体)中国北京汉族已确认亲生关系的亲子鉴定案的20个常染色体STR基因座(D3S1358、D1S1656、D13S317、Penta E、D16S539、D18S51、D2S1338、CSF1PO、Penta D、TH01、vWA、D21S11、D6S1043、D7S820、D5S818、TPOX、D8S1179、D12S391、D19S433和FGA;由于有限突变模型中未包含D6S1043的矫正参数,因此本文实际计算其余19个STR基因座的突变率)进行了调查。结果发现,所有基因座均存在突变现象,总计发生1665个突变事件,包括1614个一步突变,34个两步突变,8个三步突变和9个非整步突变。基因座特异性的平均实际突变率在三联体中为0.00007700 (TPOX)~0.00459050 (FGA),在二联体中为0.00000000 (TPOX)~0.00344850(FGA)。此外,本研究还分析了表面和实际突变率、三联体和二联体突变率、父源和母源的突变率之间的关系。研究表明,实际突变率多大于表面突变率,而且μ1”/μ2” (表面突变率)的比值通常也大于μ1/μ2 (实际突变率) (μ1”,μ1; μ2”, μ2分别是一步和两步的突变率),即更多的“隐性”突变被释放出来。而且父源和母源的三联体和二联体的突变率也有存在差异。随后,将这些突变率数据与已发表的中国其他汉族人口的相关研究进行比较,展现出了STR突变率的时间与区域差异。由于样本量大,本研究中还报告了一些少见的突变事件,例如同卵双胞胎突变和“假四步突变”等。综上所述,本研究通过大量数据获得了接近真实的STR突变率的估计值,不仅可为中国法医DNA数据库和群体遗传学数据库提供重要的基础数据,也对开展法医学个体识别、亲权鉴定和遗传学研究具有重要的意义。
关键词： 常染色体STR;有限突变模型;亲子鉴定;中国汉族人群;突变分析

Abstract
Short tandem repeat (STR) markers have been widely used in forensic paternity testing and individual identification, but the STR mutation might impact on the forensic result interpretation. Importantly, the STR mutation rate was underestimated due to ignoring the “hidden” mutation phenomenon in most similar studies. Considering this, we use Slooten and Ricciardi’s restricted mutation model based on big data to obtain more accurate mutation rates for each marker. In this paper, the mutations of 20 autosomal STRs loci (D3S1358, D1S1656, D13S317, Penta E, D16S539, D18S51, D2S1338, CSF1PO, Penta D, TH01, vWA, D21S11, D6S1043, D7S820, D5S818, TPOX, D8S1179, D12S391, D19S433, and FGA; The restricted model does not include the correction factor of D6S1043, this paper calculates remaining 19 STR loci mutation rates) were investigated in 28,313 (Total: 78,739 individuals) confirmed parentage-testing cases in Chinese Han population. As a result, total 1665 mutations were found in all loci, including 1614 one-steps, 34 two-steps, 8 three-steps, and 9 nonintegral mutations. The loci-specific average mutation rates ranged from 0.00007700 (TPOX) to 0.00459050 (FGA) in trio’s and 0.00000000 (TPOX) to 0.00344850 (FGA) in duo’s. We analyzed the relationship between mutation rates of the apparent and actual, the trio’s and duo’s, the paternal and maternal, respectively. The results demonstrated that the actual mutation rates are more than the apparent mostly, and the values of μ1”/μ2”(apparent) are also greater than μ1/μ2 (actual) commonly (μ1”, μ1; μ2”, μ2 are the mutation rates of one-step and two-step). Therefore, the “hidden” mutations are identified. In addition, the mutations rates of trio’s and duo’s, the paternal and maternal, exhibit significant difference. Next, those mutation data are used to do a comparison with the studies of other Han populations in China, which present the temporal and regional disparities. Due to the large sample size, some rare mutation events, such as monozygotic (MZ) mutation and “fake four-step mutation”, are also reported in this study. In conclusion, the estimation values of actual mutations are obtained based on big data, they can not only provide basic data for the Chinese forensic DNA and population genetics databases, but also have important significance for the development of forensic individual identification, paternity testing and genetics research.
Keywords：autosomal STR;restricted mutation model;parentage testing;Chinese Han population;mutational analysis


PDF (1741KB)元数据多维度评价相关文章导出EndNote|Ris|Bibtex收藏本文
本文引用格式
刘志勇, 任贺, 陈冲, 张京晶, 张晓梦, 石妍, 石林玉, 陈滢, 程凤, 贾莉, 陈曼, 范庆炜, 张家榕, 李万婷, 王萌春, 任子林, 刘雅诚, 倪铭, 孙宏钰, 严江伟. 基于有限突变模型和大规模数据的19个常染色体STR的实际突变率研究. 遗传[J], 2021, 43(10): 949-961 doi:10.16288/j.yczz.21-197
Zhiyong Liu. Actual mutational research of 19 autosomal STRs based on restricted mutation model and big data. Hereditas(Beijing)[J], 2021, 43(10): 949-961 doi:10.16288/j.yczz.21-197


STR (short tandem repeat)是人类基因组中非常丰富的一类遗传标记[1],由于其高度的多态性[2,3],被广泛应用于个体识别、亲子鉴定和群体遗传学等领域。然而,在大多数STR基因座中都可以观察到突变事件,一般认为复制滑脱是导致STR突变的主要机制[4]。这些突变现象往往会对法医学实践中证据权重评价产生影响,比如在计算累积父权指数(cumulative paternity index, CPI)时需要特别考虑突变现象,过高过低的突变率数据都不利于得出客观结论,因此获得STR基因座的准确突变率是非常重要的。

STR基因座突变率的实际评估计算中,往往会受抽样群体大小的影响[5]。尽管已有研究者发表大量关于STR基因座突变率的论文,但基于较少的样本量,其STR突变率的计算结果可能存在偏差。除了以上因素外,突变率的不同计算方法也会对STR实际突变率计算产生影响。在已报道的多数研究中STR突变率一般是通过直接计数法来计算的[5,6,7,8],并且将三联体和二联体的数据进行了合并。然而,该方法没有考虑“隐性突变(hidden mutation) ”现象,因为人们普遍认为当发生突变时,突变步数越少越真实,但实际并非总是如此,即某些突变会被隐藏。鉴于通过三联体和二联体的计算得到的STR基因座突变率有一定差异,因此合并三联体和二联体的突变数据可能是不合理的。即,直接计数法对于STR的实际突变率的估计有一定偏差。

针对以上抽样群体尺度与计算方法的问题,本研究采用了大量亲子鉴定数据和Slooten与Ricciardi[9]的更为严谨的有限突变模型计算方法。基于28,313例已确认亲生关系的中国汉族人群亲子鉴定案例,分别计算了三联体和二联体两种情形的19个常染色体STR基因座实际突变率。本文详细比较了表面突变率(直接计数法)和实际突变率(有限突变模型),并与其他多个群体的STR突变率进行了对比分析。最后,本文还讨论了影响STR突变率的几个因素,比如突变年龄与来源、等位基因长度、基因座杂合度等。

1 材料与方法

1.1 数据来源

本研究涉及的基因座突变数据来自2009~2019年间中国北京汉族人群的常规亲子鉴定案例的统计结果,该数据已获得北京通达首诚司法鉴定所许可使用。所有数据均采用匿名化处理,不涉及个体相关除突变信息以外的其他STR基因座具体分型,也不涉及家系情况。本次统计分析的亲子鉴定案例的CPI值均大于10,000,达到亲子鉴定的确认标准(未能确认的案例被除外),且突变的STR基因座已经过多个试剂盒验证。总计28,313例(78,739个体)亲子鉴定数据被纳入,包括22,113例三联体(trio's)和6200例二联体(duo's) (父亲-孩子：2342例;母亲-孩子：3858例)。该数据涉及的STR基因座包括CODIS 系统(D3S1358、D13S317、D16S539、D18S51、CSF1PO、TH01、vWA、D21S11、D7S820、D5S818、TPOX、D8S1179、FGA)和D1S1656、D6S1043、Penta E、Penta D、D2S1338、D12S391、D19S433。由于有限突变模型[9]中未包含D6S1043基因座的计算矫正参数,因此本研究中共包括除D6S1043之外其余19个STR基因座的突变率计算。

1.2 突变率计算

在已确认亲生关系的亲子鉴定案例中,当发现STR等位基因违反孟德尔遗传定律时,认为该STR基因座发生突变事件[10]。在计算过程中,突变步数与来源的判定根据Brinkmann等[10]和Weber等[11]的描述,本研究中所有可疑的STR突变事件均再次通过同一试剂盒和其他试剂盒进行检测,确保分型正确。所有STR基因座的三联体和二联体父源母源突变事件的统计结果见附表1。在突变率计算过程中,一方面使用直接计数法公式μ”=n/N计算表面突变率,其中n是观察到的突变事件的数目,N为减数分裂的发生次数;另一方面在得到表面突变率后,通过概率矫正因子矩阵A矫正得到实际突变率(见方程(I))。其中A_{k, l}代表观测到k步突变,实际是l步突变的概率。对于任一A_{k, l},Slooten和Ricciardi[9]基于NFI的参考数据集(2085个体)的等位基因型频率,随机组合一对夫妇,模拟生成子代的基因型,其中从母亲遗传一个没有突变的等位基因,但从父亲遗传得到一个实际有l步突变的等位基因,即已知实际突变步数,然后通过观察得到表面突变步数。上述过程通过计算机模拟重复100,000次,计算不同突变步数的频率,进而估计最终突变转移概率A_{k, l}。该模型假设：(1)当发生不确定来源的突变时,认为父源与母源发生的概率相同,均为0.5;(2)当观察到的突变事件可以同时解释为相同的突变增加或者突变减少时,该模型方法不予区分。本文采用模型中下列矩阵方程(I)计算实际突变率μ。

$\left[ \begin{matrix} \mu 0” \\ \mu 1” \\ \mu 2” \\ \vdots \\ \mu \text{k}” \\ \end{matrix} \right]=\left[ \begin{matrix} 1 & {{A}_{1,0}} & {{A}_{2,0}} & \cdots & {{A}_{\text{k},0}} \\ 0 & {{A}_{1,1}} & {{A}_{2,1}} & \cdots & {{A}_{\text{k},1}} \\ 0 & 0 & {{A}_{2,2}} & \cdots & {{A}_{\text{k},2}} \\ \vdots & \vdots & \vdots & \vdots & \vdots \\ 0 & 0 & 0 & 0 & {{A}_{\text{k,k}}} \\ \end{matrix} \right]\cdot \left[ \begin{matrix} \mu 0 \\ \mu 1 \\ \mu 2 \\ \vdots \\ \mu \text{k} \\ \end{matrix} \right]$ 方程(I)

在“方程(I) ”中,μk”为表面为k步的突变率,μ0”为表面不发生突变的概率;μk代表实际为k步的突变率,μ0为实际不发生突变的概率;Ak,l是表面突变步数为k,而实际步数为l的概率;μ”为基因座特异性表面突变率; μ为基因座特异性实际突变率。为了便于理解该模型,方程(Ⅱ)列出了方程(I)的展开式。以TH01为例,统计计算得到了表面突变率μk”后,再根据Slooten和Ricciardi[9]提供的矫正矩阵A_{k, l},可以结合方程(Ⅱ)反解出实际突变率μk,其具体的计算过程见附表2。本文采用在线计算器(https://www.numberempire.com/matrixbinarycalculator.php)进行矩阵求逆。

$\left\{ \begin{align} & \mu 0+\mu 1{{A}_{1,0}}+\mu 2{{A}_{2,0}}+\mu 3{{A}_{3,0}}=\mu 0” \\ & \mu 1{{A}_{1,1}}+\mu 2{{A}_{2,1}}+\mu 3{{A}_{3,1}}=\mu 1” \\ & \mu 2{{A}_{2,2}}+\mu 3{{A}_{3,2}}\text{=}\mu 2” \\ & \mu 3{{A}_{3,3}}=\mu 3” \\ \end{align} \right. 方程(Ⅱ)$

2 结果与分析

2.1 突变分析

本项研究中,在28,313例确认的亲子鉴定案例(22,113例三联体和6200例二联体)的20个常染色体STR基因座中共发现了1665个突变事件,所有基因座均观察到突变现象。基于有限突变模型计算得到的三联体和二联体的基因座特异性实际突变率和总突变率数据显示(表1,不包含D6S1043基因座)：三联体基因座特异性平均实际突变率范围为0.00007700 (TPOX)~0.00459050 (FGA),二联体为0.00000000 (TPOX)~0.00344850 (FGA);FGA基因座具有较高的父源突变率(0.00767100)和总突变率(0.00459050),Penta E具有较高的母源突变率(0.00237900),与Jin等[12]基于直接计数法的研究结果一致。

Table 1
表1
表1三联体与二联体STR基因座特异性实际突变率
Table 1The total STR loci-specific actual mutation rates in trio's and duo's
基因座三联体二联体
FμMμ平均μFμMμ平均μ
CSF1PO0.004491000.000898000.002694500.003795000.001152000.00247350
D12S3910.006745000.001885000.004315000.005328000.000000000.00266400
D13S3170.001972000.000899000.001435500.003105000.000943000.00202400
D16S5390.002461000.000739000.001600000.001642000.000498000.00107000
D18S510.005249000.001404000.003326500.003094000.000751000.00192250
D19S4330.002223000.001233000.001728000.001581000.000960000.00127050
D1S16560.002643000.001202000.001922500.001959000.000000000.00097950
D21S110.002557000.002012390.002284700.002033000.001234000.00163350
D2S13380.002738000.000477000.001607500.001203000.000000000.00060150
D3S13580.001605390.000734000.001169700.001496000.000978000.00123700
D5S8180.002747000.000926000.001836500.006100000.000617000.00335850
D7S8200.003250000.000509000.001879500.001498000.000455000.00097650
D8S11790.004438000.000810000.002624000.003681000.000000000.00184050
FGA0.007671000.001510000.004590500.005736000.001161000.00344850
Penta D0.001825000.000494000.001159500.002431000.000489000.00146000
Penta E0.005131000.002379000.003755000.005256000.001587000.00342150
TH010.000226110.000100000.000163060.000776000.000471000.00062350
TPOX0.000077000.000077000.000077000.000000000.000000000.00000000
vWA0.006006000.001088000.003547000.000749000.000909000.00082900
Total0.064055500.019376390.041715950.051463000.012205000.03183400
“μ”表示实际突变率;“F”表示父源;“M”表示母源。

新窗口打开|下载CSV

同卵双胞胎(monozygotic, MZ)是由单个精子和单个卵子结合形成受精卵后,分裂成两个胚胎形成的,其发生率约为1/250[13]。从理论上讲,同卵双胞胎拥有相同的基因组,然而研究者在STR[14]、拷贝数变异(CNV)[15]和单核苷酸多态性(SNP)[16]中观察到了一些差异。本次研究中,在10对MZ中发现了相同突变,包括7例男性和3例女性(附表3),突变基因座为D6S1043 (2例)、Penta E (2例)、D18S51 (2例)、FGA、D21S11、D12S391和CSF1PO,且均为一步突变。在已公开发表的论文中,尚未报道同时在MZ中发现不一致的STR突变,但Wang等[17]报道了一对MZ中其中一个个体的vWA基因座出现三等位基因模式的案例。在亲子鉴定实践中,一般较少会在同一个体观察到多个突变。如果观察到,其可能会干扰鉴定人得出亲子关系的正确结论。在本次研究的28,313例亲子鉴定案件中,共发现27例孩子的基因型中可观察到2个突变的事件,发生率为0.00095363 (27/28,313,20个常染色体STR基因座) (附表4A);同时也发现3例孩子的基因型中可观察到3个突变的事件(附表4B)。在这30例的多基因座突变案例中,本研究使用了总计超过39个基因座的多个STR试剂盒对亲子关系进行了确认。

2.2 表面与实际突变率的比较

基于直接计数法和有限突变模型,计算得到每个STR基因座的表面突变率(μ”)与实际突变率(μ),三联体的相关突变率数据见附表5,二联体相关数据见附表6。在使用有限突变模型的矩阵计算(方程(I))中,本研究发现有三个负值的实际突变率,即三联体组的D13S317和D19S433基因座的Mμ2;二联体组的D3S1358基因座的Fμ1(M,F分别代表母源和父源突变)。为了查找原因,本研究分析了D13S317和D19S433基因座在三联体中的突变特征,发现母源突变中可观察到一步和三步突变,但未观察到两步突变;同理,二联体组的D3S1358基因座仅在父源中观察到两步突变,但未观察到一步突变。即表明,实际突变率负值可能是在某STR基因座可观察到高步数突变,但未观察到某个低步数突变的情况下发生的(即非连续的突变步数现象)。换言之,即这几个负值可能在提示潜在低步数突变应该会被观察到,但实际中并没有观察到,这已经超过了有限突变模型的可解释范畴。因此在数据处理过程中,把三个负值实际突变率设为零。如图1所示,在父源和母源突变之间有一个明显的界限,特别是在三联体的一侧(图1A),这可能提示三联体与二联体的计算得到的STR突变率是有差异的,尽管这种差异可能来自于三联体与二联体本身,或者来自于样本量的差异,但可间接说明以往直接合并三联体与二联体数据进行计算存在不合理之处。此外,从该图也可发现在三联体和二联体中,TH01和TPOX基因座的表面和实际突变率都较低。与之相反的是,D12S391和FGA基因座在三联体中的表面和实际突变率均较高,而D21S11和Penta E基因座在二联体中的表面和实际突变率较高。为了进一步分析有限突变模型的校正效果,绘制了实际突变率与表面突变率比值图(μ/μ”) (图2)。该图表明,除父源的D3S1358和TH01(μ/μ”≈1)和母源的D21S11基因座(μ/μ”≈1)外,其余基因座的突变率均得到了较大程度的矫正(μ/μ”>1)。就矫正程度而言,二联体组的基因座的效应表现得更为明显(与Slooten和Ricciardi[9]的研究结果相似),尤其是CSF1PO和D5S818基因座。也就是说有限突变模型通过参数矫正的方式,将更多的隐性的突变挖掘出来。然而,在亲子鉴定过程中,需要计算有突变率参与的父权指数时,目前一般采用直接计数法计算得到平均STR突变率(比如男性取0.002;女性取0.001~0.0005)来进行。若为了更严谨的计算,可采用三联体和二联体的父源、母源或者STR基因座特异性突变率,当然也可以使用本次研究基于有限突变模型计算获得的STR基因座特异性实际突变率。从理论上讲,计算单个基因座父权指数时,采用上述两种算法得到的突变率,结果是会存在差异,但对总体的鉴定结论的影响有限。

图1

新窗口打开|下载原图ZIP|生成PPT
图1三联体和二联体表面和实际突变率热图

A：三联体情形;B：二联体情形;“μ””表示表面突变率;“μ”表示实际突变率。
Fig. 1The heatmap of apparent and actual mutation rates in trio's and duo's



图2

新窗口打开|下载原图ZIP|生成PPT
图2表面突变率与实际突变率比较

“A”和“B”代表三联体组父源(F)和母源(M)的μ/μ”比值;“C”和“D”分别代表二联体组父源(F)和母源(M)的μ/μ”比值;空白处表示比式的分母为零的情况。
Fig. 2The comparisons of apparent and actual mutation rates



2.3 三联体和二联体组突变率的比较

在以往基于直接计数法研究中,一般将三联体和二联体组的数据进行合并计算突变率。然而,三联体和二联体得到的突变率水平是有一定差异,其在2.1和2.2中已进行部分阐述。在大多数STR基因座中,三联体父源的实际突变率值大于二联体的父源相关值(图3,绿色>黑色);三联体母源的实际突变率值大于二联体的母源相关值(图3,粉色>蓝色);它们的比值表现出相似的规律(图3),特别是在vWA基因座。与此同时,本研究也发现存在相反趋势的基因座,如三联体组的D13S317、D5S818和二联体组的CSF1PO、D3S1358、TH01等,因此将三联体和二联体的数据合并是不合理的,它们之间的差异是存在且不宜忽略的。

图3

新窗口打开|下载原图ZIP|生成PPT
图3三联体和二联体中Fμ、Mμ及其比值的比较

“trio's”表示三联体组;“duo's”表示二联体组;空白处表示比式的分母为零的情况。
Fig. 3The comparisons of Fμ, Mμ and their ratios in trio's and duo's



2.4 父源和母源突变率的比较

在1665个突变事件中,共观察到1178例父源突变和293例母源突变事件,同时有194例突变来源不明(父源或母源),其中父源突变占多数,为70.75% (图4A)。在三联体组中,总Fμ为0.06405550,Mμ为0.01937639 (表1);二联体中,总Fμ为0.05146300,Mμ为0.01220500 (表1)。从以上数据可以得知,父源突变率一般较母源突变率高,在三联体和二联体组中分别约为3.3倍和4.2倍,其与以往的研究结果一致[11,12]。Jónsson等[18]对冰岛1548个三联体进行的新生突变(de novo mutations, DNMs)研究也表明父源的突变数量约是母源的4倍,这也侧面反映了本次研究的结果。本研究还比较了不同突变来源的基因座特异性突变率(图5)。如图5A显示,三联体中绝大多数基因座不同突变来源计算得到的突变率具有一定差异性,除D21S11和TPOX基因座外,父源突变率是母源突变率2倍以上。而二联体组也显示了相似的趋势(图5B)。这可能是男性与女性形成配子时经历的细胞分裂次数不同,而男性需要产生更多的配子的缘故[10]。可以预知,使用突变来源特异性基因座突变率将有助于亲子关系等法医学鉴定得出的更加客观而准确结论。另外,为了说明不同突变来源的年龄范围,本研究制作了父源(图4B)和母源(图4C)突变年龄的扇形图,从中可以看到父源突变的年龄范围比母源突变的年龄范围更广,这可能是由于男性在较长时间内可以保持一定的生殖能力,而女性不可以的原因。

图4

新窗口打开|下载原图ZIP|生成PPT
图4本次研究中突变事件来源、年龄和步数的比较

A：突变来源分布;B：父源突变年龄分布;C：母源突变年龄分布;D：突变步数分布。
Fig. 4The comparisons of mutation origin, age, and step in this research



图5

新窗口打开|下载原图ZIP|生成PPT
图5不同突变来源的基因座特异性突变率比较

A：三联体情形;B：二联体情形;空白处表示比式的分母为零的情况。
Fig. 5The comparisons of loci-specific mutation rates in different mutation origins



2.5 突变步长与突变方向分析

在所有的突变事件中,共观察到1614个一步突变、34个两步突变、8个三步突变和9个非整步突变。一步突变占总突变事件的96.94%,是主要的突变类型,逐步突变模型(stepwise mutation model, SMM)可以解释一步突变为主的现象[19,20]。在其他非一步突变中,两步突变占2.04%,三步突变占0.48%,而非整步突变占0.54% (图4D)。发生非整步突变的基因座分别为：D19S433 (3例)、D21S11 (2例)、D1S1656、D6S1043、TH01和D12S391;突变来源为父亲的6例,母亲的1例,未能确定突变来源的2例;非整步的突变方向,其中7例突变减少,1例突变增加,另有1例不能确定突变方向;关于非整步突变的步数,0.1、0.2、0.3步均有观察到,其中有5例含有0.2步突变,2例有0.3步突变,1例有0.1步突变,另有1例为0.3或者0.1步突变。就本研究观察到的9例非整步突变而言,突变减少和包含0.2步突变是占多数的,详细信息展现在附表7中。为了进一步分析有限突变模型对于STR基因座等位基因不同突变步数的突变率影响,本研究计算了三联体中发生一步突变与两步突变的比值(图6)。在父源(图6A)中,μ1”/μ2”的值普遍大于μ1/μ2,特别是在vWA基因座(本次研究：μ1”/μ2”≈104;μ1/μ2≈73),这表明由于有限突变计算模型的使用,μ1/μ2之间的差异正在减小。这些规律与Slooten与Ricciardi[9]的研究相似(μ1”/μ2”≈75;μ1/μ2≈58),也就是说更多的高步数隐性突变被挖掘出来。在母源(图6B)基本上也显示了这一规律,但发现了具有相反趋势的基因座,如D21S11、FGA和Penta E,根据该模型作者的解释[9],这可能是由于模型本身不完美的原因造成的。

图6

新窗口打开|下载原图ZIP|生成PPT
图6三联体父源与母源的突变步长分析

A：父源情形;B：母源情形;空白处表示比式的分母为零的情况。
Fig. 6The mutation step analysis of the paternal and maternal origins in trio's



尽管有限突变计算模型中不考虑突变方向(突变减少/增加),但本研究仍然统计了等位基因减少和增加的分布情况。在1665次突变事件中(附表8),总计有706次突变减少(mutation loss),703次突变增加(mutation gain),256次为不确定突变方向或者分组;突变减少与增加(loss/gain)之比约为1∶1,与Lan等[21]的研究相似。?Brinkmann等[10]报道认为突变的发生与等位基因的长度存在一定关联,长的等位基因比短的等位基因更易于发生突变;Shao等[5]认为较长的等位基因倾向于发生突变减少事件,而短的等位基因则更倾向于发生突变增加事件。为了探究等位基因长度与突变方向之间的关系,参照Ge等[22]和Xu等[23]的方法,将每个STR基因座的等位基因根据长度分为三组,即较短(short group)、中等(medium group)、较长(long group)等位基因组,详细的统计数据见附表8。从表中可计算得知在以上三个分组中,总的突变次数分别为105、1042、397,其中中等等位基因长度组发生的突变事件最多(约为67.49%),较长等位基因长度组次之(约为25.71%),而较短等位基因长度组最少(约为6.80%),也即中等及较长等位基因更易于发生突变,较短的等位基因发生突变较少,这与Xu等[23]的研究结论一致;较短等位基因组突变减少与增加(loss/gain)的比值约为0.49 (<1),中等组为0.86,较长组为1.72 (>1),本研究结果也表现出等位基因“短变长,长变短”的规律,这也得到Shao等[5]、Ge等[22]、Xu等[23]、Dupuy等[24]的支持,其可能与STR的序列结构有关。

2.6 杂合度与实际突变率的关系

一般认为,杂合度高的基因座其突变率也较高,也即突变率高可能是基因座杂合度高的一个原因。据此根据已有文献报道的杂合度数据[25],分析了每个STR基因座的杂合度和实际突变率的关系,使用Origin软件绘制线性拟合图(图7),Pearsonʼs r= 0.68454。结果也相似地表明,杂合度较低的STR基因座突变率普遍较低,如TH01、TPOX基因座,这与前人的研究相一致[12]。

图7

新窗口打开|下载原图ZIP|生成PPT
图719个常染色体STR基因座杂合度与实际突变率的关系

Fig. 7The relationship between heterozygosity and actual mutation rates for the 19 autosomal STR loci



3 讨论

不同的突变率计算方法结果可能存在潜在的偏差,比如传统的直接计数法在处理突变数据时未考虑到“隐性突变”现象[26,27]。此处以两个确认的三联体亲子鉴定为例进行说明(下文“a”为正整数)。案例一：某STR基因座母亲的分型是M=(a, a+1),孩子是C=(a, a+1),父亲是F=(a-1, a+1),假设真实的遗传情况为孩子C的“a+1”来自母亲M的“a+1”,“a”来自父亲F的“a-1”或者“a+1”一步突变后的“a”,即实际发生一步突变。从表面上看,可能会认为孩子C的“a+1”遗传自父亲F的“a+1”,而“a”遗传自母亲M的“a”,此时表面上会认为没有突变发生,即“突变发生有无认识偏差”型隐性突变(Type Ⅰ型); 案例二：某STR基因座母亲的分型是M=(a+1, a+1),孩子是C=(a, a-1),父亲是F=(a, a-1),假设真实遗传情况为孩子C的“a-1”来自母亲M的“a+1”两步突变后的“a-1”,“a”来自父亲F的“a”,即实际发生两步突变。从表面上看,可能会认为孩子C的“a-1”遗传自父亲F的“a-1”,“a”遗传自母亲M的“a+1”一步突变后的“a”,此时表面上会认为发生一步突变,因为当不知道突变真相时,一般判断为较短的突变步数,即“真实突变步数认识偏差”型隐性突变(Type Ⅱ型)。类似地,隐性突变也可以出现在二联体中。在上述案例中,表面观察到的突变被定义为表面突变(apparent mutation, ApM),而不可见的真实突变被定义为实际突变(actual mutation, AcM)。因此,如果不考虑隐性突变(比如直接计数法),突变率可能就会被低估。尽管研究者无法直接从表面的STR分型突变现象中洞察真实的突变,但Slooten和Ricciardi[9]通过突变模拟的方式研究了该隐性突变现象,并提供了可接近真实突变率的计算模型和矫正参数。该模型根据突变的等位基因是否超过STR等位基因分型标准品(ladder)的范围L,可分为有限突变模型(restricted mutation model)和无限突变模型(unrestricted mutation model),而真实的生物学过程可能介于这两个模型之间。当不知道突变的等位基因是否超过Ladder时,无限突变模型可能是合适的。然而,基于大量数据研究调查发现[9],突变的等位基因不包括在Ladder中是一种非常少见的现象,本次研究的数据也证明这一点。综合上述,鉴于使用简单的有限突变模型是具有一定合理性的,本研究拟采用该模型计算实际STR突变率。

基于以上有限突变模型计算获得的STR基因座实际突变率,本研究一方面为了观察其是否有时间差异,收集了Yan等[28]于2006年报道的与本研究具有相同来源的北京汉族突变率数据,这组数据调查了6532例亲子鉴定案例中的19,754个样本的17个常用STR基因座特异性突变率, 这些基因座均包含在本次研究中;另一方面为了观察本研究的北京汉族与中国其他汉族人群突变率的地域差异,收集了中国其他地区或省(华北、华东、华南地区和河南、上海、贵州、河北、云南、福建、广东等省)汉族群体的突变率数据(附表9)。如各基因座突变率比较的散点图所示(图8),在时间变化方面：共享的17个STR基因座中,除D21S11、TH01、TPOX和Penta D外,其余13个基因座的突变率在14年间中增幅较大(红色方框与黑色圆点),尤其是D18S51和Penta E的突变率增加将近一倍;在地域变化方面,不同汉族群体间的STR基因座特异性突变率具有一定差异,特别是在vWA、FGA、D21S11、D8S1179、Penta E和D12S391基因座,这些差异可能来自于多个地域来源样本的遗传结构不同[29,30],或是使用群体样本量不同,或使用的计算模型方法不同。

图8

新窗口打开|下载原图ZIP|生成PPT
图8中国10个不同地区或省汉族人群突变率比较

Fig. 8The Han population mutation rate comparisons of 10 different regions or provinces in China



亲子鉴定实践中,有时可能遇到假突变现象(fake mutation),其应引起足够的重视。例如,本次研究中在一例确认的三联体亲子鉴定中D13S317基因座发现其基因型不符合遗传规律的现象,称之为“假四步突变(fake four-steps mutation) ”：假设父AF=8,母亲M=12和孩子C=8 (图9)。当使用四个STR扩增试剂盒(PowerPlex®21 (Promega,美国)、MicroreaderTM21 (阅微基因,北京)、MicroreaderTM23sp-B (阅微基因,北京)、AGCU 17+1 system (中德美联,无锡))确认基因型后,得到相同的结果。起初,怀疑母亲和孩子之间可能发生了四步突变,即母亲的等位基因“12”突变为“8”。为了进一步证实这个猜想,首先使用ForenSeq DNA Signature Prep Kit (Verogen, 美国)在MiSeq FGx仪器(Illumina,美国)中进行测序,并使用ForenSeq通用分析软件(Verogen,美国)分析原始数据。测序结果显示(图9),母亲的等位基因“8”发生了丢失,其序列与孩子的等位基因“8”是一致的。为了继续寻找母亲等位基因丢失的原因,使用primer5软件设计了D13S317基因座的引物(F：TGGAAGCTGCATAGAC;R：ACTCCATACAAATGCAAGAGA),并进行了克隆测序(图9)。如红色方框所示(图9),母亲的等位基因“8”的核心区上游第17 bp处发生单碱基突变,由正常的“C”突变为“T”。以上研究表明,D13S317基因座引物结合区的突变可能是多个试剂盒未能成功扩增出母亲等位基因“8”的原因。实际中,真正的四步突变是非常少见的,当遇到可疑的四步突变时,测序有时是很必要的。

图9

新窗口打开|下载原图ZIP|生成PPT
图9基于三种方法对D13S317基因座的测序结果

Fig. 9The sequencing results of D13S317 locus using three methods



本研究在28,313例确认的亲子鉴定案例中,共发现1665个突变事件。基于以上突变数据和有限突变模型,分别计算了三联体和二联体案例中的父源与母源的19个常染色体STR基因座的特异性实际突变率。与直接计数法计算突变率相比,本次的计算结果揭示了更多的隐性突变现象,因此获得的突变率数值理论上更加接近真实情况。本研究不仅可为中国法医DNA数据库和群体遗传学数据库提供重要的基础数据,也对开展法医学个体识别、亲权鉴定和遗传学研究具有重要的意义。

(责任编委: 朱波峰)

附录：

附加材料详见文章电子版www.chinagene.cn。

Table 1 Supplementary
表1附
表1附中国汉族20个常染色体STR基因座的突变数和案例数
Table 1 SupplementaryThe number of mutations and cases in the 20 STR loci studied in Chinese Han Population
基因座三联体突变数三联体
案例数		二联体突变数二联体案例数
父源母源父源或母源总计父源母源总计父亲孩子二联体母亲孩子二联体
CSF1PO657158722,11342623423858
D12S391113242115820,56280820573406
D13S3173011105122,11342623423858
D16S539388105622,11321323423858
D18S5196221012822,11352723423858
D19S433361956022,11322423423858
D1S1656291184813,71220213612149
D21S11444179222,11333623423858
D2S133850596422,11320223423858
D3S1358296134822,11322423423858
D6S104354856720,56202220573406
D5S818401185922,11361723423858
D7S820513136722,11321323423858
D8S1179751099422,11350523423858
FGA135181917222,113931223423858
Penta D30734020,56231420573406
Penta E92401014220,562841220573406
TH01412722,11311223423858
TPOX111322,11300023423858
vWA97111612422,11312323423858
Total1109264194156722,11369299823423858
由于有限突变模型中未包含D6S1043基因座的计算矫正参数,因此我们的实际突变率计算共包括除D6S1043之外其余19个STR基因座。

新窗口打开|下载CSV

Table 2 Supplementary
表2附
表2附基于有限突变模型三联体组TH01基因座实际突变率计算实例
Table 2 Supplementary The detail mutation rate calculation example of the TH01 locus in trio's
基因座孩子父亲母亲突变
来源		突变
步数		F1”F2”F3”M1”M2”M3”
基因型
TH018, 999父亲10.5000.500
母亲
TH018,999父亲10.5000.500
母亲
TH01986, 9父亲1100000
TH016, 119, 106, 9父亲1100000
TH01979父亲2010000
TH0186,98, 9父亲1100000
TH017, 1079母亲1000100
总计(n)-410200
突变”0.0001810.00004500.0000900
F1”, F2” F3”:代表突变来源为父亲,突变步数为1,2,3 的情况. M1”, M2” M3”: 代表突变来源为母亲,突变步数为1,2,3的情况。

新窗口打开|下载CSV



Table 3 Supplementary
表3附
表3附10例同卵双生子STR基因座突变描述
Table 3 SupplementaryThe STR mutated description of monozygotic twins from 10 cases
案例孩子试剂盒父亲年龄母亲年龄孩子性别突变基因座孩子分型父亲分型母亲分型突变来源突变步数
案例1MZ1AB-Identifiler, AGCU21+13030男FGA22,2421,2322,24未知±1
MZ23030男FGA22,2421,2322,24未知±1
案例2MZ1Goldeneye 20A, AGCU21+13328男D6S104313, 2013, 1818, 21母亲-1
MZ23328男D6S104313, 2013, 1818, 21母亲-1
案例3MZ1Goldeneye 20A2123女Penta E15,2517,2412,15父亲+1
MZ22123女Penta E15,2517,2412,15父亲+1
案例4MZ1Goldeneye 20A2624男D18S5112,1614,1512,15父亲+1
MZ22624男D18S5112,1614,1512,15父亲+1
案例5MZ1Goldeneye 20A, AGCU21+14039女D21S1130,33.231,32.230,31.2父亲+1
MZ24039女D21S1130,33.231,32.230,31.2父亲+1
案例6MZ1Goldeneye 20A4440男D6S104318,21.312,22.311,18父亲-1
MZ24440男D6S104318,21.312,22.311,18父亲-1
案例7MZ1PowerPlex®213231女Penta E16,171616,18未知±1
MZ23231女Penta E16,171616,18未知±1
案例8MZ1PowerPlex®214224男D12S39117,2218,2218,19母亲-1
MZ24224男D12S39117,2218,2218,19母亲-1
案例9MZ1PowerPlex®213832男CSF1PO12,1410,1312父亲+1
MZ23832男CSF1PO12,1410,1312父亲+1
案例10MZ1PowerPlex®213131男D18S5115,1714,1614,15父亲+1
MZ23131男D18S5115,1714,1614,15父亲+1
“MZ”表示同卵双胞胎(monozygotic twins);“-”表示突变减少;“+”表示突变增加;“±”表示突变增加或者减少。

新窗口打开|下载CSV

Table 4A Supplementary
表4A附
表4A附同一个孩子观察到2个突变事件的27例案件
Table 4A SupplementaryThe 27 cases of 2 mutation events observed in one child
案例试剂盒父亲年龄母亲年龄孩子性别突变基因座孩子分型父亲分型母亲分型突变来源突变步数
案例1AB-Identifiler2422女D3S13581515,1616母亲-1
D19S43312,14.213,1414,14.2父亲-1
案例2Goldeneye 20A2522女D18S511613,1513,16父亲+1
D12S39118,2118,2018,20未知+1
案例3Goldeneye 20A3129男D7S8209, 12119, 10父亲+1
CSF1PO12, 1310, 1313父亲-1
案例4Goldeneye 20A2627男CSF1PO13,1412,149,14未知±1
FGA22,2423,2624,24.2父亲-1
案例5Goldeneye 20A5525女D3S135816,1816,1715,16父亲+1
D12S39119,2018,2119父亲-1
案例6Goldeneye 20A2622女D6S10431211,1312,13父亲±1
D8S11791414,1513,15母亲±1
案例7Goldeneye 20A3738女D12S39120,232118,23父亲-1
FGA22,2322,2422,25父亲±1
案例8Goldeneye 20A2419男D5S81811,1410,1410,12母亲±1
D6S104314,1713,1917,18父亲+1
案例9PowerPlex®213027女D21S1130,3131,31.231,33未知-1
D2S133823,2418,2419,24未知-1
案例10PowerPlex®213523女D21S1130,32.22931.2,32.2父亲+1
D12S39118,2220,2118,19父亲+1
案例11PowerPlex®214034男vWA17,181914,17父亲-1
FGA22,2421,2321,24父亲±1
案例12PowerPlex®2133NA女penta E10,2216,23NA父亲-1
D5S8181311,12NA父亲+1
案例13PowerPlex®212925女D8S11791211,1312,13父亲±1
FGA23,2422.2,2520,23父亲-1
案例14PowerPlex®213728女D18S5113,2113,1519,22母亲-1
vWA1816,1918父亲-1
案例15PowerPlex®2139NA女Penta E16,1710,18NA父亲-1
D18S5114,1918,20NA父亲±1
案例16PowerPlex®212425女D13S3179,118,109,10父亲+1
Penta E12,1513,1616母亲+1
案例17PowerPlex®213032男D18S51141413,17母亲+1
FGA23,2521,2318,22父亲-1
案例18PowerPlex®212421女TPOX8,118,108父亲+1
D8S117913,1515,1715,18父亲-2
案例19PowerPlex®212828女D6S104314,171713,18母亲+1
D16S53911,1410,119,12父亲+1
案例20PowerPlex®213733男Penta E11,2713,2611,20父亲+1
Penta D10910,13父亲+1
案例21PowerPlex®212825女D13S31710,1210,118,10父亲+1
D21S1130,3130,3129,32母亲+1
案例22PowerPlex®212623女D1S165610,151511母亲-1
D13S31718,2619,2718,19父亲-1
案例23PowerPlex®213928女D16S53910,129,1312父亲+1
D21S1130.2,31.228,30.229,32.2母亲-1
案例24PowerPlex®214541女D16S539119,1211,12父亲-1
D8S117913,1813,1712,13父亲+1
案例25PowerPlex®212829女Penta D11,148,912,14父亲+2
D21S1130,31.229,3030未知+1.2
案例26PowerPlex®213028女D13S317998母亲+1
D5S81812,1410,139,12父亲+1
案例27PowerPlex®214033男D1S165614,1615,17.312,14父亲+1
D13S3178,139,128,12父亲+1

新窗口打开|下载CSV

Table 4B Supplementary
表4B附
表4B附同一个孩子观察到3个突变事件的3例案件
Table 4B SupplementaryThe 3 cases of 3 mutation events observed in one child
案例试剂盒父亲年龄母亲年龄孩子性别突变基因座孩子分型父亲分型母亲分型突变来源突变步数
案例1Goldeneye 20A, MicroreaderTM23sp
ID Sysrem		4739男D8S117912,141314父亲-1
D17S129015,1915,2015,18未知±1
FGA22,232422,24父亲-1
案例2PowerPlex®21, MicroreaderTM23sp
ID Sysrem		3926男D6S104312,1717,2013,14母亲-1
D8S113218,2118,1922-1
FGA19,2420,2519,25父亲-1
案例3PowerPlex®21, MicroreaderTM23sp
ID Sysrem		3433男D2S133819,211919,20母亲+1
D10S1435131312,14±1
FGA21,2522,2421,23父亲+1
“-”表示突变减少;“+”表示突变增加;“±”表示突变增加或者减少;“NA”表示无相关数据。

新窗口打开|下载CSV

Table 5A Supplementary
表5A附
表5A附三联体组的表面突变率
Table 5A SupplementaryThe apparent mutation rates in trio's
基因座Fμ”Fμ1”Fμ2”Fμ3”Mμ”Mμ1”Mμ2”Mμ3”
CSF1PO0.003278610.003278610.000000000.000000000.000655720.000655720.000000000.00000000
D12S3910.006006230.005908960.000097270.000000000.001677850.001677850.000000000.00000000
D13S3170.001582780.001537560.000045220.000000000.000723560.000678330.000000000.00004522
D16S5390.001944560.001944560.000000000.000000000.000587890.000542670.000045220.00000000
D18S510.004567450.004431780.000090440.000045220.001221000.001175780.000045220.00000000
D19S4330.001741060.001695840.000045220.000000000.000972280.000927060.000000000.00004522
D1S16560.002406650.002260790.000145860.000000000.001093930.001093930.000000000.00000000
D21S110.002148060.002148060.000000000.000000000.002012390.001921950.000090440.00000000
D2S13380.002464610.002464610.000000000.000000000.000429610.000429610.000000000.00000000
D3S13580.001605390.001605390.000000000.000000000.000565280.000565280.000000000.00000000
D5S8180.001989780.001944560.000045220.000000000.000678330.000633110.000045220.00000000
D7S8200.002600280.002600280.000000000.000000000.000407000.000407000.000000000.00000000
D8S11790.003595170.003595170.000000000.000000000.000655720.000655720.000000000.00000000
FGA0.006534620.006263280.000180890.000090440.001288830.001198390.000045220.00004522
Penta D0.001531950.001337420.000145900.000048630.000413380.000364750.000048630.00000000
Penta E0.004717440.004498590.000218850.000000000.002188500.002066920.000072950.00004863
TH010.000226110.000180890.000045220.000000000.000090440.000090440.000000000.00000000
TPOX0.000067830.000067830.000000000.000000000.000067830.000067830.000000000.00000000
vWA0.004748340.004703120.000045220.000000000.000859220.000859220.000000000.00000000

新窗口打开|下载CSV

Table 5B Supplementary
表5B附
表5B附三联体组的实际突变率
Table 5B SupplementaryThe actual mutation rates in trio's
基因座FμFμ1Fμ2Fμ3MμMμ1Mμ2Mμ3
CSF1PO0.004491000.004491000.000000000.000000000.000898000.000898000.000000000.00000000
D12S3910.006745000.006604000.000141000.000000000.001885000.001885000.000000000.00000000
D13S3170.001972000.001891000.000081000.000000000.000899000.00083200-0.000043000.00011000
D16S5390.002461000.002461000.000000000.000000000.000739000.000660000.000079000.00000000
D18S510.005249000.005046000.000119000.000084000.001404000.001334000.000070000.00000000
D19S4330.002223000.002152000.000071000.000000000.001233000.00118200-0.000024000.00007500
D1S16560.002643000.002451000.000192000.000000000.001202000.001202000.000000000.00000000
D21S110.002557000.002557000.000000000.000000000.002012390.001921950.000090440.00000000
D2S13380.002738000.002738000.000000000.000000000.000477000.000477000.000000000.00000000
D3S13580.001605390.001605390.000000000.000000000.000734000.000734000.000000000.00000000
D5S8180.002747000.002669000.000078000.000000000.000926000.000848000.000078000.00000000
D7S8200.003250000.003250000.000000000.000000000.000509000.000509000.000000000.00000000
D8S11790.004438000.004438000.000000000.000000000.000810000.000810000.000000000.00000000
FGA0.007671000.007252000.000238000.000181000.001510000.001373000.000047000.00009000
Penta D0.001825000.001521000.000214000.000090000.000494000.000414000.000080000.00000000
Penta E0.005131000.004847000.000284000.000000000.002379000.002222000.000077000.00008000
TH010.000226110.000180890.000045220.000000000.000100000.000100000.000000000.00000000
TPOX0.000077000.000077000.000000000.000000000.000077000.000077000.000000000.00000000
vWA0.006006000.005925000.000081000.000000000.001088000.001088000.000000000.00000000
“Fμ”表示父源突变率;“Mμ”表示母源突变率;“Fμ”与“Mμ”后的数字表示突变步数;标红的数字表示计算出的实际突变率为负值,详见文中。

新窗口打开|下载CSV

Table 6A Supplementary
表6A附
表6A附二联体组表面突变率
Table 6A SupplementaryThe apparent mutation rates in duo's
基因座Fμ”Fμ1”Fμ2”Fμ3”Mμ”Mμ1”Mμ2”Mμ3”
CSF1PO0.001707940.001707940.000000000.000000000.000518400.000518400.000000000.00000000
D12S3910.003889160.003889160.000000000.000000000.000000000.000000000.000000000.00000000
D13S3170.001707940.001707940.000000000.000000000.000518400.000518400.000000000.00000000
D16S5390.000853970.000853970.000000000.000000000.000259200.000259200.000000000.00000000
D18S510.002134930.002134930.000000000.000000000.000518400.000518400.000000000.00000000
D19S4330.000853970.000853970.000000000.000000000.000518400.000518400.000000000.00000000
D1S16560.001469510.001469510.000000000.000000000.000000000.000000000.000000000.00000000
D21S110.001280960.001280960.000000000.000000000.000777600.000777600.000000000.00000000
D2S13380.000853970.000853970.000000000.000000000.000000000.000000000.000000000.00000000
D3S13580.000853970.000426990.000426990.000000000.000518400.000518400.000000000.00000000
D5S8180.002561910.002561910.000000000.000000000.000259200.000259200.000000000.00000000
D7S8200.000853970.000853970.000000000.000000000.000259200.000259200.000000000.00000000
D8S11790.002134930.002134930.000000000.000000000.000000000.000000000.000000000.00000000
FGA0.003842870.003842870.000000000.000000000.000777600.000777600.000000000.00000000
Penta D0.001458430.001458430.000000000.000000000.000293600.000293600.000000000.00000000
Penta E0.003889160.003889160.000000000.000000000.001174400.001174400.000000000.00000000
TH010.000426990.000426990.000000000.000000000.000259200.000259200.000000000.00000000
TPOX0.000000000.000000000.000000000.000000000.000000000.000000000.000000000.00000000
vWA0.000426990.000426990.000000000.000000000.000518400.000518400.000000000.00000000

新窗口打开|下载CSV

Table 6B Supplementary
表6B附
表6B附二联体组的实际突变率
Table 6B SupplementaryThe actual mutation rates in duo's
基因座FμFμ1Fμ2Fμ3MμMμ1Mμ2Mμ3
CSF1PO0.003795000.003795000.000000000.000000000.001152000.001152000.000000000.00000000
D12S3910.005328000.005328000.000000000.000000000.000000000.000000000.000000000.00000000
D13S3170.003105000.003105000.000000000.000000000.000943000.000943000.000000000.00000000
D16S5390.001642000.001642000.000000000.000000000.000498000.000498000.000000000.00000000
D18S510.003094000.003094000.000000000.000000000.000751000.000751000.000000000.00000000
D19S4330.001581000.001581000.000000000.000000000.000960000.000960000.000000000.00000000
D1S16560.001959000.001959000.000000000.000000000.000000000.000000000.000000000.00000000
D21S110.002033000.002033000.000000000.000000000.001234000.001234000.000000000.00000000
D2S13380.001203000.001203000.000000000.000000000.000000000.000000000.000000000.00000000
D3S13580.00149600-0.001554000.003050000.000000000.000978000.000978000.000000000.00000000
D5S8180.006100000.006100000.000000000.000000000.000617000.000617000.000000000.00000000
D7S8200.001498000.001498000.000000000.000000000.000455000.000455000.000000000.00000000
D8S11790.003681000.003681000.000000000.000000000.000000000.000000000.000000000.00000000
FGA0.005736000.005736000.000000000.000000000.001161000.001161000.000000000.00000000
Penta D0.002431000.002431000.000000000.000000000.000489000.000489000.000000000.00000000
Penta E0.005256000.005256000.000000000.000000000.001587000.001587000.000000000.00000000
TH010.000776000.000776000.000000000.000000000.000471000.000471000.000000000.00000000
TPOX0.000000000.000000000.000000000.000000000.000000000.000000000.000000000.00000000
vWA0.000749000.000749000.000000000.000000000.000909000.000909000.000000000.00000000
“Fμ”表示父源突变率;“Mμ”表示母源突变率;“Fμ”与“Mμ”后的数字表示突变步数;标红的数字表示计算出的实际突变率为负值,详见文中。

新窗口打开|下载CSV

Table 7 Supplementary
表7附
表7附非整步突变事件汇总
Table 7 SupplementaryThe total nonintegral mutations events
基因座基因型突变来源可能的突变步数
孩子父亲母亲
D19S4331314.213, 15.2父亲-1.2
D19S43313, 1615, 16.213, 14-0.2
D19S4331315.2--2.2
D1S16561616.3, 1715, 16-0.3
D6S104317, 23.318, 2411, 17-0.1
D21S1130, 3229, 30.232, 32.2-0.2
TH0196, 99.3母亲-0.3
D21S1130, 31.229, 3030未知+1.2
D12S39119.3, 2019, 2020, 21+0.3/-0.1
加粗的字体表示突变的等位基因。

新窗口打开|下载CSV

Table 8 Supplementary
表8附
表8附等位基因突变增加与减少事件汇总
Table 8 SupplementaryThe total mutation events of allele mutation gain and loss
基因座较短组 Short group中等组 Medium group较长组 Long group未确定分组
Undeterminable
group		Total
ARNGNLNUUMARNGNLNUUMARNGNLNUUM
CSF1PO?7-9.11001?9.2-11.164010?12-1531283622093
D12S391?14-18.3811120?19-2336648108?24-288701523166
D13S317?5-70000?8-122211235?13-15512017557
D16S539?6-80000?9-122110435?13-16611118659
D18S51?10-139009?13.2-205432591?21-2515180332135
D19S433?9-11.22002?12-162425251?16.2-18.238011064
D1S1656?10-120101?13-17.32020242?18-19.31405250
D21S11?27-2987015?29.2-33.12631360?33.2-35.2212014998
D2S1338?16-19.25106?19.3-23.3217028?24-26919129366
D3S1358?12-130000?14-182122952?19-200000052
D6S1043?9-122002?12.2-181315432?18.2-22.31319234169
D5S818?7-90000?9.2-1327161356?13.2-1528010066
D7S820?7-92114?9.1-1158215?11.1-1412193341770
D8S1179?8-104004?11-1534191164?16-18722130199
FGA13-18.2110219-24.251531812225-28152934713184
Penta E?5-15167528?16-21.35437697?22-33.17140218154
PentaD?5-95106?10-181812535?19-221001244
TH01?5-71001?8-122608?13-15000009
TPOX?7-80101?8.1-11.21001?12-13010103
vWA?14-151113?16-19414019100?19.3-214110159127
Total -65328105-4974321131042-141242143971211665
AR: allele range (等位基因长度范围); NG: number of mutation gain (突变增加的数量); NL: number of mutation loss (突变减少的数量); NU: number of uncertain mutation (不能确定突变方向的数量); NM: number of total mutation (突变总数)。

新窗口打开|下载CSV

Table 9 Supplementary
表9附
表9附中国其他省或区汉族人群的突变率汇总
Table 9 SupplementaryThe Han populations mutation rates from other provinces or regions in China
基因座本次
研究		Yan等
研究[1]		华北
汉族[2]		华东
汉族[2]		华南
汉族[3]		河南
汉族[4]		上海
汉族[5]		贵州
汉族[6]		河北
汉族[7]		云南
汉族[8]		福建
汉族[9]		广东
汉族[9]
n=78739n=19754n=16612n=4494n=6303n=2475n=12413n=3879n=5006n=10941n=3279n=4094
CSF1PO0.00270.00110.00140.00080.00050.00240.00260.00040.00100.00090.00080.0014
D12S3910.0043-0.00290.00000.00280.00120.00170.00130.00240.00270.00210.0037
D13S3170.00140.00060.00070.00080.00100.00000.00170.00090.00070.00030.00000.0017
D16S5390.00160.00060.00070.00000.00030.00180.00060.00090.00070.00050.00040.0017
D18S510.00330.00060.00220.00170.00080.00180.00290.00090.00100.00150.00170.0017
D19S4330.00170.00040.00040.00170.00080.00000.00090.00040.00030.00140.00000.0010
D1S16560.0019---0.00000.0012-0.0013-0.00080.00000.0003
D21S110.00230.00210.00070.00040.00180.00300.00200.00090.00070.00210.00120.0010
D2S13380.00160.00040.00170.00040.00050.00000.00110.00040.00070.00090.00080.0014
D3S13580.00120.00050.00120.00040.00050.00060.00150.00180.00030.00090.00000.0027
D5S8180.00180.00060.00100.00080.00080.00000.00110.00090.00070.00120.00040.0003
D7S8200.00190.00080.00140.00000.00050.00060.00030.00040.00030.00060.00040.0020
D8S11790.00260.00100.00170.00120.00130.00060.00290.00040.00140.00140.00080.0027
FGA0.00460.00220.00240.00460.00280.00490.00460.00130.00310.00240.00250.0027
Penta D0.00120.00090.00020.00080.00100.00060.00060.00040.00070.00060.00040.0007
Penta E0.00380.00090.00320.00120.00210.00000.00320.00180.00140.00200.00250.0020
TH010.00020.00020.00000.00000.00000.00060.0002-0.00000.00050.00000.0003
TPOX0.00010.00010.00020.00000.00000.00000.0003-0.00000.00030.00040.0003
vWA0.00350.00160.00180.00040.00180.00300.00150.00270.00100.00270.00210.0027
Total0.0013-0.00130.0008--0.0016-0.01630.00120.00080.0014
本次研究中显示为三联体平均突变率;“n”表示样本量;“-”表示无相关数据。

新窗口打开|下载CSV


参考文献 原文顺序
文献年度倒序
文中引用次数倒序
被引期刊影响因子

Subramanian S, Mishra RK, Singh L . Genome-wide analysis of microsatellite repeats in humans: their abundance and density in specific genomic regions
Genome Biol, 2003,4(2):R13.
[本文引用: 2]

B?r W, Brinkmann B, Budowle B, Carracedo A, Gill P, Lincoln P, Mayr W, Olaisen B . DNA recommendations. Further report of the DNA Commission of the ISFH regarding the use of short tandem repeat systems. International Society for Forensic Haemogenetics,
Int J Legal Med, 1997,110(4):175-176.
[本文引用: 3]

Butler JM . Short tandem repeat typing technologies used in human identity testing
Biotechniques, 2007,43(4):2-5.
[本文引用: 2]

Kornberg A, Bertsch LL, Jackson JF, Khorana HG . Enzymatic synthesis of deoxyribonucleic acid, XVI. oligonucleotides as templates and the mechanism of their replication
Proc Natl Acad Sci USA, 1964,51(2):315-323.
[本文引用: 2]

Shao CC, Lin MX, Zhou ZH, Zhou YQ, Shen YW, Xue AM, Zhou HG, Tang QQ, Xie JH . Mutation analysis of 19 autosomal short tandem repeats in Chinese Han population from Shanghai
Int J Legal Med, 2016,130(6):1439-1444.
[本文引用: 5]

Li HX, Peng D, Wang Y, Wu RG, Zhang YM, Li R, Sun HY . Evaluation of genetic parameters of 23 autosomal STR loci in a southern Chinese Han population
Ann Hum Biol, 2018,45(4):359-364.
[本文引用: 2]

Wang HD, Kand B, Su N, He M, Zhang B, Guo YX, Zhu BF, Liao SX, Zeng ZS . Evaluation of the genetic parameters and mutation analysis of 22 STR loci in the central Chinese Han population,
Int J Legal Med, 2017,131(1):103-105.
[本文引用: 2]

Qu N, Zhang XC, Liang H, Ou XL . Analysis of genetic polymorphisms and mutations at 23 autosomal STR loci in Guangdong Han population
Forensic Sci Int Genet, 2019,38:e16-e17.
[本文引用: 2]

Slooten K, Ricciardi F . Estimation of mutation probabilities for autosomal STR markers
Forensic Sci Int Genet, 2013,7(3):337-344.
[本文引用: 11]

Brinkmann B, Klintschar M, Neuhuber F, Hühne J, Rolf B . Mutation rate in human microsatellites: influence of the structure and length of the tandem repeat
Am J Hum Genet, 1998,62(6):1408-1415.
[本文引用: 4]

Weber JL, Wong C . Mutation of human short tandem repeats
Hum Mol Genet, 1993,2(8):1123-1128.
[本文引用: 2]

Jin B, Su Q, Luo HB, Li YB, Wu J, Yan J, Hou YP, Liang WB, Zhang L . Mutational analysis of 33 autosomal short tandem repeat (STR) loci in southwest Chinese Han population based on trio parentage testing
Forensic Sci Int Genet, 2016,23:86-90.
[本文引用: 3]

Taylor MJ, Fisk NM . Prenatal diagnosis in multiple pregnancy
Baillieres Best Pract Res Clin Obstet Gynaecol, 2000,14(4):663-675.
[本文引用: 1]

Wurmb-Schwark NV, Schwark T, Christiansen L, Lorenz D, Oehmichen M . The use of different multiplex PCRs for twin zygosity determination and its application in forensic trace analysis
Leg Med(Tokyo), 2004,6(2):125-130.
[本文引用: 1]

Xu YJ, Li TT, Pu T, Cao RX, Long F, Chen S, Sun K, Xu R . Copy number variants and exome sequencing analysis in six pairs of Chinese monozygotic twins discordant for congenital heart disease
Twin Res Hum Genet, 2017,20(6):521-532.
[本文引用: 1]

Weber Lehmann J, Schilling E, Gradl G, Richter DC, Wiehler J, Rolf B . Finding the needle in the haystack: differentiating “identical” twins in paternity testing and forensics by ultra-deep next generation sequencing
Forensic Sci Int Genet, 2014,9:42-46.
[本文引用: 1]

Wang LF, Yang Y, Zhang XN, Quan XL, Wu YM . Tri-allelic pattern of short tandem repeats identifies the murderer among identical twins and suggests an embryonic mutational origin
Forensic Sci Int Genet, 2015,16:239-245.
[本文引用: 1]

Jónsson H, Sulem P, Kehr B, Kristmundsdottir S, Zink F, Hjartarson E, Hardarson MT, Hjorleifsson KE, Eggertsson HP, Gudjonsson SA, Ward LD, Arnadottir GA, Helgason EA, Helgason H, Gylfason A, Jonasdottir A, Jonasdottir A, Rafnar T, Frigge M, Stacey SN, Magnusson OT, Thorsteinsdottir U, Masson G, Kong A, Halldorsson BV, Helgason A, Gudbjartsson DF, Stefansson K . Parental influence on human germline de novo mutations in 1,548 trios from Iceland
Nature, 2017,549(7673):519-522.
[本文引用: 1]

Valdes AM, Slatkin M, Freimer NB . Allele frequencies at microsatellite loci: the stepwise mutation model revisited
Genetics, 1993,133(3):737-749.
[本文引用: 1]

Klintschar M, Dauber EM, Ricci U, Cerri N, Immel UD, Kleiber M, Mayr WR . Haplotype studies support slippage as the mechanism of germline mutations in short tandem repeats
Electrophoresis, 2010,25(20):3344-3348.
[本文引用: 1]

Lan Q, Wang HD, Shen CM, Guo YX, Yin CY, Xie T, Fang YT, Zhou YS, Zhu BF . Mutability analysis towards 21 STR loci included in the AGCU 21+1 kit in Chinese Han population
Int J Legal Med, 2018,132(5):1287-1291.
[本文引用: 1]

Ge JY, Budowle B, Aranda XG, Planz JV, Eisenberg AJ, Chakraborty R . Mutation rates at Y chromosome short tandem repeats in Texas populations,
Forensic Sci Int Genet, 2009,3(3):179-184.
[本文引用: 2]

Xu W, Wang YQ, Zhang DD, Wang DX, Zhou L, Ye XL, Zhu CG, Shi YZ . Mutation analysis of 21 autosomal short tandem repeats in Han population from Hunan, China
Ann Hum Biol. 2019,46(3):254-260.
[本文引用: 3]

Dupuy BM, Stenersen M, Egeland T, Olaisen B . Y-chromosomal microsatellite mutation rates: differences in mutation rate between and within loci,
Hum Mutat, 2004,23(2):117-124.
[本文引用: 1]

Xie BB, Chen L, Yang YR, Lv YX, Chen J, Shi Y, Chen C, Zhao HY, Yu ZL, Liu YC, Fang XD, Yan JW . Genetic distribution of 39 STR loci in 1027 unrelated Han individuals from northern China
Forensic Sci Int Genet, 2015,19:205-206.
[本文引用: 1]

Chakraborty R, Stivers DN, Zhong YX . Estimation of mutation rates from parentage exclusion data: applications to STR and VNTR loci
Mutat Res, 1996,354(1):41-48.
[本文引用: 1]

Vicard P, Dawid AP . A statistical treatment of biases affecting the estimation of mutation rates
Mutat Res, 2004,547(1-2):19-33.
[本文引用: 1]

Yan JW, Liu YC, Tang H, Zhang QX, Huo ZY, Hu SN, Yu J . Mutations at 17 STR loci in Chinese population
Forensic Sci Int, 2006,162(1-3):53-54.
[本文引用: 1]

Lu DJ, Liu QL, Wu WW, Zhao H . Mutation analysis of 24 short tandem repeats in Chinese Han population
Int J Legal Med, 2012,126(2):331-335.
[本文引用: 1]

Qian XQ, Yin CY, Ji Q, Li K, Fan HT, Yu YF, Bu FL, Hu LL, Wang JW, Mu HF, Haigh S, Chen F . Mutation rate analysis at 19 autosomal microsatellites
Electrophoresis, 2015,36(14):1633-1639.
[本文引用: 1]

附录涉及的参考文献(References):



Yan JW, Liu YC, Tang H, Zhang QX, Huo ZY, Hu SN, Yu J . Mutations at 17 STR loci in Chinese population
Forensic Sci Int, 2006,162(1-3):53-54.
[本文引用: 2]

Qian XQ, Yin CY, Ji Q, Li K, Fan HT, Yu YF, Bu FL, Hu LL, Wang JW, Mu HF . Mutation rate analysis at 19 autosomal microsatellites
Electrophoresis, 2015,36(14):1633-1639.
[本文引用: 3]

Li HX, Peng D, Wang Y, Wu RG, Zhang YM, Li R, Sun HY . Evaluation of genetic parameters of 23 autosomal STR loci in a Southern Chinese Han population,
Ann Hum Biol, 2018,45(4):359-364.
[本文引用: 2]

Wang. HD, Kang B, Su N, He M, Zhang B, Guo YX, Zhu BF, Liao. SX, Zeng ZS. Evaluation of the genetic parameters and mutation analysis of 22 STR loci in the central Chinese Han population
Int J Legal Med, 2016,131(1):103-105.
[本文引用: 2]

Shao CC, Lin MX, Zhou ZH, Zhou YQ, Shen YW, Xue AM, Zhou HG, Tang QQ, Xie JH . Mutation analysis of 19 autosomal short tandem repeats in Chinese Han population from Shanghai
Int J Legal Med, 2016,130(6):1439-1444.
[本文引用: 5]

Yang MQ, Ren Z, Ji JY, Zhou H, Zhang HL, Dai JL, Wang J, Huang J . Population genetic data and mutations of 22 autosomal STR loci in Guizhou Han population
Forensic Sci Int Genet, 2017,29:e29-e30.
[本文引用: 2]

Sun LJ, Shi K, Tan L, Zhang Q, Fu LH, Zhang XJ, Fu GP, Li SJ, Cong B . Analysis of genetic polymorphisms and mutations at 19 STR loci in Hebei Han population
Forensic Sci Int Genet, 2017,31:e50-e51.
[本文引用: 2]

Zhang XF, Liu LL, Xie RF, Wang GY, Shi Y, Gu T, Hu LP, Nie SJ . Population data and mutation rates of 20 autosomal STR loci in a Chinese Han population from Yunnan Province, Southwest China
Int J Legal Med, 2018,132(4):1083-1085.
[本文引用: 2]

Qu N, Zhang XC, Liang H, XL Ou . Analysis of genetic polymorphisms and mutations at 23 autosomal STR loci in Guangdong Han population
Forensic Sci Int Genet, 2019,38:e16-e17.
[本文引用: 11]

相关话题/基因 计算 数据 突变 北京

Free考研考试FreeKaoYan.Com
欢迎来到Free考研考试，"为实现人生的Free而奋斗"

© 2020 FreeKaoYan! . All rights reserved.