2. 福建傲芯种业科技有限公司, 漳州 363000;
3. 上杭傲农槐猪产业发展有限公司, 龙岩 364200;
4. 福建傲农生物科技集团股份有限公司, 漳州 363000;
5. 江西农业大学 猪遗传改良与养殖技术国家重点实验室, 南昌 330045
2. Fujian Aoxin Seed Industry Technology Co. Ltd., Zhangzhou 363000, China;
3. Shanghang Aonong Huai Pig Industry Development Co. LTD., Longyan 364200, China;
4. Fujian Aonong BiologicaI Science and Technology Group Co. Ltd., Zhangzhou 363000, China;
5. National Key Laboratory for Swine Genetic Improvement and Production Technology, Jiangxi Agricultural University, Nanchang 330045, China
近交是血缘相关或遗传组成相似个体间的交配繁殖。在畜牧生产中,群体规模有限,近交在所难免。近交水平的增加会导致基因组中有害纯合子的增加[1],也会导致群体遗传多样性的降低[2-3]。近交系数(inbreeding coefficient)表示一个个体中随机选择两个等位基因位点是血缘同源(identical by descent, IBD)的概率[4],近交系数也展示了通过近交降低个体杂合性的比例[5-6]。传统近交系数的计算方法是基于完整的系谱记录信息[7-8]。但在实际生产中,系谱信息的缺乏、不完整、错误率高等都是常见的现象,因此无法准确评估近交系数和亲缘关系[9]。近年来,随着分子生物学技术的高速发展和测序成本的大幅下降,基因组测序、芯片等高通量变异检测技术被广泛应用于家畜的近交程度和亲缘关系的研究中,为家畜育种提供了重要的技术参考[10-11]。
使用基因组数据计算近交系数主要基于SNP纯合性法[12-13],基于SNP间方差法[14-15]和基于连续性纯合片段法[16-17] (runs of homozygosity,ROH)。基于SNP间方差法来自于基因预测中的亲缘关系矩阵,该关系矩阵取代了最佳线性无偏估计(best linear unbiased predicted,BLUP)中的分子亲缘关系矩阵。根据分子亲缘关系矩阵对角线元素等于1加上近交系数,所以一般认为基因型亲缘关系矩阵的对角线为1加上相应个体的基因型近交系数[18-19]。而基于SNP纯合性法与基于SNP间方差法其个体近交系数可能出现负值[20]。ROH是由亲代将来自共同祖先同一段单倍型传递给后代所形成,一般属于IBD单倍型片段。基于ROH的近交系数能同时避免系谱近交系数低估和SNP近交系高估的问题[21]。同一个体不同基因组近交系数的估计结果可能存在较大的差异。因此,比较不同基因组个体近交系数间的差异与相关性,有助于准确评估群体近交程度,维持群体遗传多样性。
槐猪是我国优良地方猪种,具有母性好、繁殖能力强、耐粗饲、营养价值高等优良特性。近年来,随着外来瘦肉型商品猪种大量冲击中国生猪市场以及非洲猪瘟的侵袭,槐猪和绝大多数中国地方猪种一样面临有效群体减少的风险。为此,本试验拟使用槐猪保种群作为研究对象,采用不同的计算方法估计槐猪个体的近交系数,揭示不同参数条件下基因组近交系数的结果及相对优势,以期为评估猪个体近交水平、防止近交衰退和优化选种选配方案提供技术参考。
1 材料与方法 1.1 试验群体本试验所用群体来自福建傲农保种场提供的247头槐猪个体(300日龄),其中包括19头公猪,228头母猪,个体间基本无亲缘相关。使用磁珠法组织基因组提取试剂盒从猪耳组织和公猪原精中提取基因组DNA,检测DNA浓度及琼脂糖凝胶电泳符合分型标准后,在北京康普森生物技术有限公司进行猪“中芯一号”50K芯片分型。
1.2 基因型填充及质控为避免缺失基因型对近交系数计算的影响,本试验使用Beagle v5.4[22]对原始芯片数据进行基因型填充。使用plink v1.90[23]对芯片填充数据进行质控,选择次等位基因频率(minor allele frequency,MAF)≥0.05、哈迪-温伯格平衡(Hardy-Weinberg equilibrium,HWE)检验P < 1×10-5的SNPs进行近交系数计算。
1.3 基因组近交系数计算1.3.1 近交系数计算方式 基于SNP纯合性评估近交系数为Fhom1[23]、Fhom2[24],其计算公式如下:
$ { Fhom } 1=\frac{O(h o m)-E(h o m)}{S-E(h o m)}=1-\frac{\sum _{k=1}^S x_k\left(2-x_k\right)}{\sum _{k=1}^S 2 p_k\left(1-p_k\right)} ; \\ { Fhom } 2=\frac{O(h o m)-E( { hom })}{1-E( { hom })}=1-\frac{1}{S} \sum _{k=1}^S \frac{x_k\left(2-x_k\right)}{2 p_k\left(1-p_k\right)} $ |
在Fhom1中O(hom)、E(hom)、S分别表示纯合子观察数目、纯合子期望数目、SNP总数。在Fhom2中O(hom)、E(hom)分别表示纯合子比例、纯合子期望比例。xk为个体在第k个SNP上的次等位基因剂量数,pk为群体在第k个SNP上的次等位基因频率,下同。
根据Vanraden[18]的基因组亲缘关系矩阵的对角线元素计算得出的近交系数为Fvan1、Fvan2, 其计算公式如下:
$\begin{gathered} F {van} 1=\frac{1}{S} \sum\nolimits_{k=1}^S\left(\frac{\left(x_k-2 p_k\right)^2}{2 p_k\left(1-p_k\right)}-1\right) ; \\ F {van} 2=\frac{\sum\left(x_k-2 p_k\right)^2}{2 \sum p_k\left(1-p_k\right)} \end{gathered} $ |
根据Yang等[24]的基因组亲缘关系矩阵的对角线元素计算得出的近交系数为Fyang, 其计算公式如下:
$ { Fyang }=\frac{1}{S} \sum _{k=1}^s \frac{x_k^2-\left(1+2 p_k\right) x_k+2 p_k^2}{2 p_k\left(1-p_k\right)}。$ |
在基因组某一段区域内,当一定数量一定密度的SNPs表现为纯合时,可以判定该区域存在连续性纯合片段。利用ROH计算近交系数为Froh[25-26], 其公式如下:
$ { Froh }=\sum L_{{ROH}} / L_{ {auto }} {, } $ |
其中,LROH为常染色体上ROH片段的长度,Lauto为常染色体的总长度。ROH检测使用plink v1.90[23]软件,使用滑动窗口的方法对常染色体进行检测,具体的检测参数:至少连续50个SNPs;连续SNP间的距离小于1 Mb;最多允许ROH中有1个杂合;窗口阈值为0.01;密度为0.01 SNP/kb;ROH长度大于1 Mb、5 Mb分别记为Froh1、Froh5。
1.3.2 SNP数目对近交系数评估的影响 为比较不同SNPs数目对近交系数评估的影响,试验从质控后的SNPs中随机选择100、200、500、1 000、2 000、5 000、10 000和20 000个SNPs,计算群体的Fhom1、Fhom2、Fvan1、Fvan2和Fyang,试验重复100次。
1.3.3 不同近交系数统计描述及相关性分析 使用质控后的SNPs(31 545个)计算群体中所有个体Fhom1、Fhom2、Fvan1、Fvan2、Fyang、Froh1和Froh5值,计算不同近交系数评估方式在群体中的最小值、最大值、平均值和方差。为比较不同近交系数计算方式的相关性,本试验使用Person相关系数和Spearman秩相关系数计算不同近交系数评估方式间的相关性。
2 结果 2.1 芯片数据填充与质控本试验使用“中芯一号”对247头槐猪进行基因分型,共包含57 466个SNPs。利用Beagle v5.4对槐猪的芯片数据进行基因型填充,共获得55 156个无缺失SNPs。剔除了829个HWE检验P值小于1×10-5和22 782个MAF < 0.05的SNPs,最终保留31 545个SNPs用于后续分析。
2.2 SNP数目对近交系数评估的影响为探究SNP数目对近交系数评估的影响,试验随机选择了100、200、500、1 000、2 000、5 000、10 000、20 000个SNPs来计算群体的近交系数(Fhom1、Fhom2、Fvan1、Fvan2和Fyang),重复试验100次。结果如图 1所示,虚线为使用质控合格的31 545个SNPs计算的近交系数值。试验发现,随着SNP数目的增加,个体近交系数Fhom1、Fhom2、Fvan1和Fyang的极差逐渐减小,而当SNP数目达到10 000个左右时近交系数趋于稳定,其值与使用所有合格SNPs所计算的近交系数相近。值得注意的是,Fvan2值趋近于0,且随着SNP数目的增加基本无明显变化。
随后,使用合格的SNPs(31 545个)来计算槐猪的近交系数:Fhom1、Fhom2、Fvan1、Fvan2、Fyang、Froh1和Froh2(表 1)。结果显示,Fhom1、Fhom2、Fvan1、Fvan2和Fyang的值分别为-0.216 5~0.344 0、-0.354 8~0.321 7、-0.322 8~ 0.892 1、-0.243 6~0.790 9、-0.189 4~0.543 3,近交系数均存在负值。Froh1和Froh2的值分别为0.008 7~0.406 9和0.002 2~0.402 4。
为衡量不同近交系数评估方式间的相关性,使用Person相关系数进行相关性分析(图 2)。结果显示,Fhom1与Fhom2、Froh1、Froh5间有较强的相关性,分别为0.860 6、0.911 8和0.911 1;Fvan1与Fyang、Fvan2的相关系数分别为0.836 4和0.988 1;Fyang与Fvan2的相关性为0.874 4;Froh1与Froh5之间的相关性高达0.998 2。Fhom2与Fvan1、Fvan2呈现负相关,分别为-0.279 0和-0.192 0。
试验进一步通过Spearman秩相关系数计算不同近交系数评估方式间的相关性(图 3)。结果显示,Fhom1与Fhom2、Froh1、Froh5有较强的正相关,分别为0.828 0、0.825 6、0.821 4;Fvan1与Fvan2间有较强的正相关,为0.985 9;Fyang与Fvan2间有较强的正相关,为0.850 5;Froh1与Froh5间有较强的正相关,为0.995 2。Fvan1与Fhom1、Fhom2呈负相关,分别为-0.010 5和-0.489 0;Fvan2与Fhom2呈负相关,为-0.383 5。
在本研究中,随着SNP数目增加,Fhom1、Fhom2、Fvan1、Fyang均值极差逐渐减少,当SNP数目达10 000时平均近交系数基本稳定,其值与利用所有质控合格的31 545个SNPs计算的近交系数值接近。这与此前Wang[27]发现当基因组较大时,SNP标记达10 000比系谱更能准确的评估个体近交系数结果相符。从Fhom1、Fhom2、Fvan1、Fyang等近交系数评估公式可知,个体近交系数为其不同位点近交系数的均值,当可用位点数较少时,评估个体近交系数值与真实值偏差较大,当可用位点数较多时,评估个体近交系数值与真实值偏差较小。Kardos等[28]发现,当SNP较小时,近交系数评估会出现较大的偏差。而Fvan2值随着SNP数目的增加基本无变化,且Fvan2值一直接近于0。从Fvan2近交系数评估公式可知,SNP取值方式对近交系数评估会产生较大的波动。Silió等[29]使用不同数目的SNP计算系谱近交系数与基因组近交系数间的相关性发现,随着SNP数目增加基因组近交系数与系谱近交系数的相关性增强,这也在一定程度上说明随着SNP数目的增加,分子近交系数的随机性减少。
3.2 不同近交系数评估方式计算个体近交系数的统计在本研究中,槐猪的近交系数Fhom1、Fhom2、Fvan1、Fvan2、Fyang计算结果分别为-0.216 5~0.344 0、-0.354 8~0.321 7、-0.322 8~0.892 1、-0.243 6~0.790 9、-0.189 4~0.543 3,近交系数存在负值。Froh1、Froh2个体近交系数为0.008 7~ 0.406 9、0.002 2~0.402 4,其个体近交系数不存在负值。这可能是因为槐猪群体从不同地方收集组建而成,其来源较广,而在随机配种的大群体中,以当前群体为参照,远交个体间的近交系数为负值[30]。其中,Fhom1、Fhom2的近交系数计算方法会出现-∞~1的值,对于单个位点,其纯合子近交系数为1;对于单个位点,其杂合子近交系数为-∞~ 0[20, 31]。Fvan1、Fvan2、Fyang的近交系数计算方式会出现-1~∞的值,对于单个位点,其纯合子近交系数为-1~∞;对于单个位点,其杂合子近交系数为0~∞[20, 31]。由于赋予稀有位点更高的权重,个体稀有等位基因纯合时,基因组近交系数比普通等位基因更大[18]。而Froh1、Froh5的近交系数计算方式会将值控制在0~1之间[32-34]。Fhom1、Fhom2、Fvan1、Fvan2、Fyang的值可能超出了概率和相关性的允许范围[35-37],但作为近交系数,这可以理解为变异性的损益与基础种群的变异性成正比,负值表示获得了变异性,正值表示失去了变异性[20]。
3.3 不同近交系数评估方式相关性在本研究中,通过Person相关与Spearman秩相关分析可知,Fhom1与Fhom2、Froh1、Froh5间有较强的正相关;Fvan1与Fvan2间有较强的正相关;Fyang与Fvan2间有较强的的正相关;Froh1与Froh5间有较强的正相关;Fhom2与Fvan1、Fvan2呈现负相关。Caballero等[38]比较了100、1 000、10 000等3个不同规模群体的近交系数,均发现Fhom与Froh1、Froh5间有较强的正相关。Zhang[31]比较了50K芯片与序列数据的近交系数,发现Froh与Fhom间有较强的正相关且与Fvan和Fyang的相关性较差。其中,Froh被发现与系谱近交系数的相关性最强,可提供更好的估计效果[39-40]。这与本试验结果互相印证、结果相似。
在比较不同近交系数的取值范围与相关性后,试验推荐Froh作为基因组近交系数的评估方式,一方面是因为Froh的取值范围在0~1之间与系谱近交系数的取值范围接近,符合评估习惯;另一方面,Froh同Fhom1、Fhom2、Fvan1、Fvan2、Fyang等基因组近交系数评估方式有着一定的正相关性,评估结果适用度高。然而Froh受SNP的数目与密度影响较大,对基因分型成本要求较高。
4 结论SNP数目较少时,个体近交系数随机误差较大,SNP数目较多时,个体近交系数与真实近交系数偏差较小。Fhom、Fvan、Fyang等近交系数均会出现负值,Froh近交系数在0 ~ 1之间。Froh同Fhom、Fvan、Fyang等近交系数评估方式有着一定的正相关,Fhom与Fvan间呈现负相关。当SNP数目与密度充足时,可采用Froh作为基因组近交系数度量值,可为基因组近交评估在地方猪的分子保种应用中提供技术参考。
[1] |
DERKS M F L, LOPES M S, BOSSE M, et al. Balancing selection on a recessive lethal deletion with pleiotropic effects on two neighboring genes in the porcine genome[J]. PLoS Genet, 2018, 14(9): e1007661. DOI:10.1371/journal.pgen.1007661 |
[2] |
HULSEGGE I, CALUS M, HOVING-BOLINK R, et al. Impact of merging commercial breeding lines on the genetic diversity of Landrace pigs[J]. Genet Sel Evol, 2019, 51(1): 60. DOI:10.1186/s12711-019-0502-6 |
[3] |
TANG G Q, XUE J, LIAN M J, et al. Inbreeding and genetic diversity in three imported swine breeds in China using pedigree data[J]. Asian-Australas J Anim Sci, 2013, 26(6): 755-765. DOI:10.5713/ajas.2012.12645 |
[4] |
WRIGHT S. Coefficients of inbreeding and relationship[J]. Am Nat, 1922, 56(645): 330-338. DOI:10.1086/279872 |
[5] |
IVERSEN M W, NORDBØ Ø, GJERLAUG-ENGER E, et al. Effects of heterozygosity on performance of purebred and crossbred pigs[J]. Genet Sel Evol, 2019, 51(1): 8. DOI:10.1186/s12711-019-0450-1 |
[6] |
WANG L G, MU Y L, XU L Y, et al. Genomic analysis reveals specific patterns of homozygosity and heterozygosity in inbred pigs[J]. Animals, 2019, 9(6): 314. DOI:10.3390/ani9060314 |
[7] |
VOGES S, DISTL O. Inbreeding trends and pedigree analysis of Bavarian mountain hounds, Hanoverian hounds and Tyrolean hounds[J]. J Anim Breed Genet, 2009, 126(5): 357-365. DOI:10.1111/j.1439-0388.2009.00800.x |
[8] |
BERNARDES P A, GROSSI D A, SAVEGNAGO R P, et al. Population structure of Tabapuã beef cattle using pedigree analysis[J]. Livest Sci, 2016, 187: 96-101. DOI:10.1016/j.livsci.2016.03.002 |
[9] |
TAYLOR H R, KARDOS M D, RAMSTAD K M, et al. Valid estimates of individual inbreeding coefficients from marker-based pedigrees are not feasible in wild populations with low allelic diversity[J]. Conserv Genet, 2015, 16(4): 901-913. DOI:10.1007/s10592-015-0709-1 |
[10] |
WU F, SUN H, LU S X, et al. Genetic diversity and selection signatures within diannan small-ear pigs revealed by next-generation sequencing[J]. Front Genet, 2020, 11: 733. DOI:10.3389/fgene.2020.00733 |
[11] |
WANG Y P, ZHAO X Y, WANG C, et al. High-density single nucleotide polymorphism chip-based conservation genetic analysis of indigenous pig breeds from Shandong Province, China[J]. Anim Biosci, 2021, 34(7): 1123-1133. DOI:10.5713/ajas.20.0339 |
[12] |
MCGIVNEY B A, HAN H G, CORDUFF L R, et al. Genomic inbreeding trends, influential sire lines and selection in the global thoroughbred horse population[J]. Sci Rep, 2020, 10(1): 466. DOI:10.1038/s41598-019-57389-5 |
[13] |
MORALES-GONZÁLEZ E, SAURA M, FERNÁNDEZ A, et al. Evaluating different genomic coancestry matrices for managing genetic variability in turbot[J]. Aquaculture, 2020, 520: 734985. DOI:10.1016/j.aquaculture.2020.734985 |
[14] |
SANTIAGO E, NOVO I, PARDIÑAS A F, et al. Recent demographic history inferred by high-resolution analysis of linkage disequilibrium[J]. Mol Biol Evol, 2020, 37(12): 3642-3653. DOI:10.1093/molbev/msaa169 |
[15] |
ALEMU S W, KADRI N K, HARLAND C, et al. An evaluation of inbreeding measures using a whole-genome sequenced cattle pedigree[J]. Heredity, 2021, 126(3): 410-423. DOI:10.1038/s41437-020-00383-9 |
[16] |
GURGUL A, SZMATOŁA T, TOPOLSKI P, et al. The use of runs of homozygosity for estimation of recent inbreeding in Holstein cattle[J]. J Appl Genet, 2016, 57(4): 527-530. DOI:10.1007/s13353-016-0337-6 |
[17] |
FORUTAN M, MAHYARI S A, BAES C, et al. Inbreeding and runs of homozygosity before and after genomic selection in North American Holstein cattle[J]. BMC Genomics, 2018, 19(1): 98. DOI:10.1186/s12864-018-4453-z |
[18] |
VANRADEN P M. Efficient methods to compute genomic predictions[J]. J Dairy Sci, 2008, 91(11): 4414-4423. DOI:10.3168/jds.2007-0980 |
[19] |
YANG J, BENYAMIN B, MCEVOY B P, et al. Common SNPs explain a large proportion of the heritability for human height[J]. Nat Genet, 2010, 42(7): 565-569. DOI:10.1038/ng.608 |
[20] |
VILLANUEVA B, FERNÁNDEZ A, SAURA M, et al. The value of genomic relationship matrices to estimate levels of inbreeding[J]. Genet Sel Evol, 2021, 53(1): 42. DOI:10.1186/s12711-021-00635-0 |
[21] |
SHI L Y, WANG L G, LIU J X, et al. Estimation of inbreeding and identification of regions under heavy selection based on runs of homozygosity in a Large White pig population[J]. J Anim Sci Biotechnol, 2020, 11(1): 46. DOI:10.1186/s40104-020-00447-0 |
[22] |
BROWNING B L, ZHOU Y, BROWNING S R. A one-penny imputed genome from next-generation reference panels[J]. Am J Human Genet, 2018, 103(3): 338-348. DOI:10.1016/j.ajhg.2018.07.015 |
[23] |
PURCELL S, NEALE B, TODD-BROWN K, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses[J]. Am J Human Genet, 2007, 81(3): 559-575. DOI:10.1086/519795 |
[24] |
YANG J, LEE S H, GODDARD M E, et al. GCTA: a tool for genome-wide complex trait analysis[J]. Am J Human Genet, 2011, 88(1): 76-82. |
[25] |
GOMEZ-RAYA L, RODRÍGUEZ C, BARRAGÁN C, et al. Genomic inbreeding coefficients based on the distribution of the length of runs of homozygosity in a closed line of Iberian pigs[J]. Genet Sel Evol, 2015, 47: 81. DOI:10.1186/s12711-015-0153-1 |
[26] |
WU X D, ZHOU R, ZHANG W, et al. Genome-wide scan for runs of homozygosity identifies candidate genes in Wannan Black pigs[J]. Anim Biosci, 2021, 34(12): 1895-1902. DOI:10.5713/ab.20.0679 |
[27] |
WANG J L. Pedigrees or markers: which are better in estimating relatedness and inbreeding coefficient?[J]. Theor Popul Biol, 2016, 107: 4-13. DOI:10.1016/j.tpb.2015.08.006 |
[28] |
KARDOS M, LUIKART G, ALLENDORF F W. Measuring individual inbreeding in the age of genomics: marker-based measures are better than pedigrees[J]. Heredity, 2015, 115(1): 63-72. |
[29] |
SILIÓ L, RODRÍGUEZ M C, FERNÁNDEZ A, et al. Measuring inbreeding and inbreeding depression on pig growth from pedigree or SNP-derived metrics[J]. J Anim Breed Genet, 2013, 130(5): 349-360. |
[30] |
WANG J. Marker-based estimates of relatedness and inbreeding coefficients: an assessment of current methods[J]. J Evol Biol, 2014, 27(3): 518-530. |
[31] |
ZHANG Q Q, CALUS M P L, GULDBRANDTSEN B, et al. Estimation of inbreeding using pedigree, 50k SNP chip genotypes and full sequence data in three cattle breeds[J]. BMC Genet, 2015, 16(1): 88. |
[32] |
SCHIAVO G, BOVO S, BERTOLINI F, et al. Comparative evaluation of genomic inbreeding parameters in seven commercial and autochthonous pig breeds[J]. Animal, 2020, 14(5): 910-920. |
[33] |
SZMATOŁA T, JASIELCZUK I, SEMIK-GURGUL E, et al. Detection of runs of homozygosity in conserved and commercial pig breeds in Poland[J]. J Anim Breed Genet, 2020, 137(6): 571-580. |
[34] |
WU X D, ZHOU R, WANG Y L, et al. Genome-wide scan for runs of homozygosity in Asian wild boars and Anqing six-end-white pigs[J]. Anim Genet, 2022, 53(6): 867-871. |
[35] |
SCHÄLER J, KRVGER B, THALLER G, et al. Comparison of ancestral, partial, and genomic inbreeding in a local pig breed to achieve genetic diversity[J]. Conserv Genet Resour, 2020, 12(1): 77-86. |
[36] |
DOS SANTOS W B, SCHETTINI G P, FONSECA M G, et al. Fine-scale estimation of inbreeding rates, runs of homozygosity and genome-wide heterozygosity levels in the Mangalarga Marchador horse breed[J]. J Anim Breed Genet, 2021, 138(2): 161-173. |
[37] |
GHOREISHIFAR S M, MORADI-SHAHRBABAK H, FALLAHI M H, et al. Genomic measures of inbreeding coefficients and genome-wide scan for runs of homozygosity islands in iranian river buffalo, Bubalus bubalis[J]. BMC Genet, 2020, 21(1): 16. |
[38] |
CABALLERO A, VILLANUEVA B, DRUET T. On the estimation of inbreeding depression using different measures of inbreeding from molecular markers[J]. Evol Appl, 2021, 14(2): 416-428. |
[39] |
VASILIADIS D, METZGER J, DISTL O. Demographic assessment of the dalmatian dog-effective population size, linkage disequilibrium and inbreeding coefficients[J]. Canine Med Genet, 2020, 7(1): 3. |
[40] |
GAZAL S, SAHBATOU M, PERDRY H, et al. Inbreeding coefficient estimation with dense SNP data: comparison of strategies and application to HapMap Ⅲ[J]. Hum Hered, 2014, 77(1-4): 49-62. |
(编辑 郭云雁)