长期的人工选择使猪的生产性能得到显著提高, 与选择相关的基因组区域也随之发生特定遗传变异表征(选择信号)。不同类型品种所受到选择强度不一, 选择信号亦不相同, 选择性清扫分析已逐渐成为选择信号的主要检测手段。文章基于商用型大白猪(n=45)和地方猪品种通城猪(n=45)的猪60K SNP芯片分型数据, 借助遗传分化系数Fst法进行选择信号检测分析。利用gPLINK软件设定质控标准, 共计34 304个SNPs被筛选出用于统计分析。使用Genepop软件包计算两个猪品种之间的遗传分化参数Fst, 所得Fst平均值为0.3209。选取Fst > 0.7036(即占总Fst值数目的1%), 共计344个SNPs被选择出来。SNP位置注释显示这些位点涉及到79个候选基因(Sus scrofa Build 9)。利用在线软件Ingenuity Pathway Analysis对候选基因的生物学通路进行网络分析, 发现它们多与生长繁殖及免疫应答有关, 如NCOA6、ERBB4、RUNX2和APOB等基因。研究结果为进行猪产肉、抗病等性状候选基因和致因突变深入挖掘提供了有益参考。
The production performance of pigs has been significantly improved due to long-term artificial selection, and the specific variation characterizations (selection signatures) emerged from the selected genome regions. Different types of breeds are subjected to different selection intensities and had different selection signatures. Selective sweep analysis is one of major methods to detect the selection signatures. In this study, based on the 60K BeadChip genotyping data of both commercial Large White (n=45) and local Tongcheng pigs (n=45), genetic differentiation coefficient Fst was applied to detect the selection signatures. Using gPLINK software to set quality control standards, a total of 34 304 SNPs were selected for statistical analysis. Fst values between two breeds were estimated with Genepop package and the average Fst value was 0.3209. Setting Fst>0.7036 (1% of total number of Fst values) as selection threshold, 344 SNPs were obtained and SNP location annotation indicated that there were 79 candidate genes (Sus scrofa Build 9). Furthermore, network analysis was performed using Ingenuity Pathway Analysis and the preliminary results suggested that most genes were in-volved in growth, reproduction, and immune response, such as NCOA6, ERBB4, RUNX2, and APOB genes. The findings from this study will contribute to further identification of candidate genes and causal mutations implying for meat production and disease resistance in pig.
[1] Giuffra E, Kijas JMH, Amarger V, Carlborg Ö, Jeon JT, Andersson L. The origin of the domestic pig: independent domestication and subsequent introgression. Genetics, 2000, 154(4): 1785-1791.
[2] Nielsen R, Hellmann I, Hubisz M, Bustamante C, Clark AG. Recent and ongoing selection in the human genome. Nat Rev Genet, 2007, 8(11): 857-868.
[3] Hudson RR, Kreitman M, Aguadé M. A test of neutral molecular evolution based on nucleotide data. Genetics, 1987, 116(1): 153-159.
[4] Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics, 1989, 123(3): 585-595.
[5] Fu YX, Li WH. Statistical tests of neutrality of mutations. Genetics, 1993, 133(3): 693-709.
[6] Fay JC, Wu CI. Hitchhiking under positive Darwinian selection. Genetics, 2000, 155(3): 1405-1413.
[7] Li WH, Wu CI, Luo CC. A new method for estimating synonymous and nonsynonymous rates of nucleotide sub-stitution considering the relative likelihood of nucleotide and codon changes. Mol Biol Evol, 1985, 2(2): 150-174.
[8] McDonald JH, Kreitman M. Adaptive protein evolution at the Adh locus in Drosophila. Nature, 1991, 351(6328): 652-654.
[9] Sabeti PC, Reich DE, Higgins JM, Levine HZ, Richter DJ, Schaffner SF, Gabriel SB, Platko JV, Patterson NJ, McDonald GJ, Ackerman HC, Campbell SJ, Altshuler D, Cooper R, Kwiatkowski D, Ward R, Lander ES. Detecting recent positive selection in the human genome from hap-lotype structure. Nature, 2002, 419(6909): 832-837.
[10] Lewontin RC, Krakauer J. Distribution of gene frequency as a test of the theory of the selective neutrality of polymorphisms. Genetics, 1973, 74(1): 175-195.
[11] Wright S. Isolation by distance. Genetics, 1943, 28(2): 114-138.
[12] Karlsson EK, Baranowska I, Wade CM, Salmon Hillbertz NH, Zody MC, Anderson N, Biagi TM, Patterson N, Piel-berg GR, Kulbokas EJ 3rd, Comstock KE, Keller ET, Mesirov JP, von Euler H, Kämpe O, Hedhammar Å, Lander ES, Andersson G, Andersson L, Lindblad-Toh K. Efficient mapping of mendelian traits in dogs through genome-wide association. Nat Genet, 2007, 39(11): 1321-1328.
[13] Pariset L, Joost S, Marsan PA, Valentini A, Econogene Consortium (EC). Landscape genomics and biased FST approaches reveal single nucleotide polymorphisms under selection in goat breeds of North-East Mediterranean. BMC Genet, 2009, 10: 7.
[14] Wahlberg P, Carlborg Ö, Foglio M, Tordoir X, Syvänen AC, Lathrop M, Gut IG, Siegel PB, Andersson L. Genetic analysis of an F2 intercross between two chicken lines divergently selected for body-weight. BMC Genomics, 2009, 10(1): 248.
[15] Flori L, Fritz S, Jaffrézic F, Boussaha M, Gut I, Heath S, Foulley JL, Gautier M. The genome response to artificial selection: a case study in dairy cattle. PloS One, 2009, 4(8): e6595.
[16] Kijas JW, Lenstra JA, Hayes B, Boitard S, Porto Neto LR, San Cristobal M, Servin B, McCulloch R, Whan V, Giet-zen K, Paiva S, Barendse W, Ciani E, Raadsma H, McE-wan J, Dalrymple B; International Sheep Genomics Con-sortium Members. Genome-wide analysis of the world’s sheep breeds reveals high levels of historic mixture and strong recent selection. PLoS Biol, 2012, 10(2): e1001258.
[17] Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet, 2007, 81(3): 559-575.
[18] Rousset F. GENEPOP’007: a complete re-implementation of the GENEPOP software for Windows and Linux. Mol Ecol Resour, 2008, 8(1): 103-106.
[19] Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution, 1984, 38(6): 1358-1370.
[20] Wright S. Ev