全基因组关联分析鉴别白色杜洛克×二花脸F2资源家系中猪阴囊疝易感基因位点的鉴定

doi:10.3864/j.issn.0578-1752.2014.14.017

摘要/Abstract

摘要： 【目的】通过基于单倍型的病例-对照全基因组关联分析(GWAS)鉴别白色杜洛克×二花脸F2资源家系中影响猪阴囊疝的易感位点及位置功能候选基因，并在远缘纯种猪阴囊疝三联体核心家系（Trio）群体中进行重复验证。【方法】利用Illumina porcine 60K SNP芯片对1 020头白色杜洛克×二花脸F2资源家系阴囊疝群体(包含19个阴囊疝患病个体)进行扫描获取基因型。通过Plink v1.07软件对基因型数据进行质量控制；剔除个体基因型检出率< 90%的个体，剔除SNP位点检出率< 90%、哈迪温伯格平衡卡方检验P≤10-3和最小等位基因频率< 0.05及性染色体上和无法定位的SNP位点。质控合格的SNP位点用PHASEBOOK构建出F2资源家系中每个个体的单倍型，并采用隐马尔可夫模型将其归类到数目预先定义的祖先单倍型中。使用基于广义线性混合模型的GLASCOW软件进行利用祖先单倍型的病例-对照全基因组关联分析。采用保守的Bonferroni校正方法得到基因组显著水平和染色体显著水平的阈值。对F2资源家系中达到染色体显著水平的易感SNP位点在包含237个患病个体的远缘纯种猪阴囊疝三联体核心家系(Trio)群体中进行同样的质控，并利用Haploview软件对质控合格的SNP位点分别展开基于单点和基于单倍型的传递不平衡分析(TDT)，进行重复验证。【结果】白色杜洛克×二花脸F2资源家系中全部个体的38 033个SNP标记通过质控。在GWAS分析中，共发现108个达染色体显著水平(P<2.63×10-5)的猪阴囊疝关联SNP位点，分别位于染色体2、8和17上，最强相关的SNP位点位于17号染色体上。在远缘三联体核心家系群中，724个个体的96个SNP位点通过质控。对质控合格的SNP位点进行基于单点的TDT分析，结果发现5个显著相关的SNP位点得到重复验证(P<0.05)，其中2号染色体上有1个SNP位点和17号染色体上有4个SNP位点，分别位于IQGAP2，CHMP4B，SERINC3，ZNF334 等4个基因内或其上下游。基于单倍型的TDT验证分析发现两个易感单倍型框，其中与基于单点TDT验证分析有重合的仅有SSC17上CHMP4B内及下游的两个易感位点。结合疝气发生机制及TDT分析结果，推测IQGAP2和CHMP4B可能是影响猪阴囊疝发生的两个重要的位置功能候选基因。【结论】本研究利用基于小样本的GWAS分析和较大样本群验证分析，在SSC2和SSC17上分别鉴别1个和4个阴囊疝易感位点，在SSC17上鉴别到两个易感单倍型。易感位点附近的2个基因IQGAP2和CHMP4B可能是影响猪阴囊疝发生的位置功能候选基因，值得进一步研究。

关键词: 猪 , 全基因组关联分析 , 传递不平衡分析 , 阴囊疝 , 易感基因

Abstract: 【Objective】Isolation of susceptible loci and positional candidate genes associated with pig scrotal hernia through case-control genome-wide association study (GWAS) based on haplotypes in a F2 White Duroc × Erhualian intercross and it’s replication in outbred trio families.【Method】Phenotypic data and Illumina porcine 60K SNP genotypes were collected from 1020 individuals in the White Duroc × Erhualian intercross (containing 19 affected pigs) and outbred trio families (containing 237 scrotal hernia pigs). Quality control was carried out using PLINK v1.07 for each population， separately. Samples of low (<90%) call rate were removed. SNP markers were removed if they had genotype-missing rates ＞0.1, or Hardy-Weinberg P≤10-3(based on Chi-squared test) or minor allele frequencies (MAF) ＜0.05, and SNPs on sex chromosomes or with unknown positions were also excluded. PHASEBOOK package was used to infer haplotypes from the genotyped individuals and to assign them to K ancestral haplotype clusters by hidden monto carlo markov models. The software GLASCOW based on generalized linear mixed models was used to carry out GWAS for pig scrotal hernia in the F2 intercross by the ancestral haplotype and significant SNPs were chosen to perform transmission disequilibrium test (TDT) based on both single SNP markers and haplotypes in the outbred population (containing 237 affected pigs) by Haploview. 【Result】After quality control, all samples in F2 population and 724 samples in outbred populations passed the filter and a final set of 38 033 SNPs and 96 SNPs were selected for GWAS and TDT analysis. In total, 108 SNPs with Bonferroni chromosome-wise significance (P＜2.63×10-5) were identified to be associated with scrotal hernia in F2 intercross, and located on SSC2, SSC8 and SSC17, respectively. The most promising SNP was located on SSC17. Five significant SNPs, located in on nearby IQGAP2，CHMP4B，SERINC3，ZNF334 genes, were replicated (located on SSC2 and SSC17) in outbred trio families by TDT analysis. TDT analysis based on haplotype revealed that two haplotypes contribute to risk of scrotal hernia, one of which was overlapped with two significant SNP on SSC17 in the TDT analysis based on single marker. IQGAP2 and CHMP4B genes close to the SNPs most associated with pig scrotal hernia were selected as positional candidate genes according to its biological and physiological functions.【Conclusion】A GWAS with small samples was conducted and variation and replication test was performed in a larger populations, which revealed 5 susceptible loci for pig scrotal hernia both in the F2 intercross and the outbred populations. IQGAP2 and CHMP4B may play a role in the risk of scrotal hernia, and deeper analysis should be carried out in the future study.

Key words: pig , genome-wide association study , transmission disequilibrium test , scrotal hernia , susceptibility gene

宿英1, 阮国荣2, 龙毅1, 杨斌1, 张志燕1, 邓玮芸1, 吴丽花1, 吕显山1, 艾华水1, 肖石军1, 任军1, 黄路生1, 丁能水1. 全基因组关联分析鉴别白色杜洛克×二花脸F2资源家系中猪阴囊疝易感基因位点的鉴定[J]. 中国农业科学, 2014, 47(14): 2872-2880.

SU Ying-1, RUAN Guo-Rong-2, LONG Yi-1, YANG Bin-1, ZHANG Zhi-Yan-1, DENG Wei-Yun-1, WU Li-Hua-1, 吕Xian-Shan-1 , AI Hua-Shui-1, XIAO Shi-Jun-1, REN Jun-1, HUANG Lu-Sheng-1, DING Neng-Shui-1. Genome-Wide Association Study Reveals Candidate Susceptibility Loci for Pig Scrotal Hernia Using Both F2 Intercross and Outbred Populations[J]. Scientia Agricultura Sinica, 2014, 47(14): 2872-2880.

0
/ / 推荐

导出引用管理器 EndNote|Reference Manager|ProCite|BibTeX|RefWorks

链接本文: https://www.chinaagrisci.com/CN/10.3864/j.issn.0578-1752.2014.14.017

https://www.chinaagrisci.com/CN/Y2014/V47/I14/2872

参考文献

[1]Walters J R. Have we forgotten about inherited disease//AGBU Pig Genetics Workshop, 2010: 27-28.

[2]孙淑东, 殷凤振, 刘文明, 来长海. 仔猪腹股沟阴囊疝的预防及治疗措施. 山东畜牧兽医, 2011, 32(8): 30-31.

Sun S D, Yin F Z, Liu W M, Lai C H. Piglets scrotal/inguinal hernia prevention and treatment measures. Shandong Journal of Animal Husbandry and Veterinary Science, 2011, 32 (8): 30-31. (in Chinese)

[3]Chang K S, Chang J H. Animal models of pediatric surgical diseases. Pediatric Surgery International 1994, 9(5/6): 307-322.

[4]Gatphayak K, Chongkasikit N, Charoensook R, Laenoi W, Vearasilp T, Sardsud V, Knorr C, ter Meulen U, Brenig B. Present situation of porcine hernia inguinalis/scrotalis in Thailand. The Global Food & Product Chain-Dynamics, Innovalions, Conficts, Strategies, Tropentag Dctober 11-13, 2005.

[5]Grindflek E, Moe M, Taubert H, Simianer H, Lien S, Moen T. Genome-wide linkage analysis of inguinal hernia in pigs using affected sib pairs. BMC Genetics, 2006, 7(1): 25.

[6]Ding N S, Mao H R, Guo Y M, Ren J, Xiao S J, Wu G Z, Shen H Q, Wu L H, Ruan G F, Brenig B. A genome-wide scan reveals candidate susceptibility loci for pig hernias in an intercross between White Duroc and Erhualian. Journal of Animal Science, 2009, 87(8): 2469-2474.

[7]Du F X, Mathialagan N, Dyer C, Grosz M, Messer L, Clutter A, Wang T, Lohuis M, Byatt J. Mapping genes affecting scrotal hernia condition in domestic pigs. Journal of Dairy Science, 2004, 87: 453-453.

[8]Du Z Q, Zhao X, Vukasinovic N, Rodriguez F, Clutter AC, Rothschild M F. Association and haplotype analyses of positional candidate genes in five genomic regions linked to scrotal hernia in commercial pig lines. PloS One, 2009, 4(3): e4837.

[9]Tanyel F, Müftüoglu S, Dagdeviren A, Kaymaz F, Büyükpamukçu N. Myofibroblasts defined by electron microscopy suggest the dedifferentiation of smooth muscle within the sac walls associated with congenital inguinal hernia. BJU International, 2001, 87(3): 251-255.

[10]Tanyel F C, Erdem S, Büyükpamukçu N, Tan E. Smooth muscle within incomplete obliterations of processus vaginalis lacks apoptotic nuclei. Urologia Internationalis, 2002, 69(1): 42-45.

[11]Rodrigues Junior A J, Rodrigues C J, Cunha A C, Jin Y. Quantitative analysis of collagen and elastic fibers in the transversalis fascia in direct and indirect inguinal hernia. Revista do Hospital das Clínicas, 2002, 57(6): 265-270.

[12]Rodrigues C J, Yoo J H, Rodrigues Junior A J. Elastin (ELN) gene point mutation in patients with inguinal hernia. Genetics and Molecular Biology, 2006, 29(1): 45-46.

[13]Andersson L, Georges M. Domestic-animal genomics: deciphering the genetics of complex traits. Nature Reviews Genetics, 2004, 5(3): 202-212.

[14]兰旅涛, 郭源梅, 陈从英, 杨斌, 毛辉荣, 任军, 周利华. 在白色杜洛克× 二花脸 F2 资源家系中定位影响猪 210 日龄 8 个体尺性状的 QTL. 中国农业科学, 2010, 43(15): 3214-3220.

Lan L T, Guo Y M, Chen C Y, Yang B, Mao H R, Ren J, Zhou L H. QTL Mapping for eight traits related to porcine body dimensions at 210 days in a White Duroc × Erhualian F2 resource population. Scientia Agricultura Sinica, 2010, 43(15): 3214-3220. (in Chinese)

[15]Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira M A, Bender D, Maller J, Sklar P, De Bakker P I, Daly M J. PLINK: a tool set for whole-genome association and population-based linkage analyses. The American Journal of Human Genetics, 2007, 81(3): 559-575.

[16]Druet T, Georges M. A hidden Markov model combining linkage and linkage disequilibrium information for haplotype reconstruction and quantitative trait locus fine mapping. Genetics, 2010, 184(3): 789-798.

[17]Zhang Z, Guillaume F, Sartelet A, Charlier C, Georges M, Farnir F, Druet T. Ancestral haplotype-based association mapping with generalized linear mixed models accounting for stratification. Bioinformatics, 2012, 28(19): 2467-2473.

[18]Yang Q, Cui J, Chazaro I, Cupples L A, Demissie S. Power and type I error rate of false discovery rate approaches in genome-wide association studies. BMC Genetics, 2005, 6(Suppl 1): 134.

[19]Barrett J C, Fry B, Maller J, Daly M. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 2005, 21(2): 263-265.

[20]Horvath S, Xu X, Laird N M. The family based association test method: strategies for studying general genotype–phenotype associations. European Journal of Human Genetics, 2001, 9(4): 301-306.

[21]Klein R J, Zeiss C, Chew E Y, Tsai J Y, Sackler R S, Haynes C, Henning A K, SanGiovanni J P, Mane S M, Mayne S T. Complement factor H polymorphism in age-related macular degeneration. Science, 2005, 308(5720): 385-389.

[22]McCarthy M I, Hirschhorn J N. Genome-wide association studies: potential next steps on a genetic journey. Human Molecular Genetics, 2008, 17(R2): R156-R165.

[23]Brown M D, Sacks D B. IQGAP1 in cellular signaling: bridging the GAP. Trends in Cell Biology, 2006, 16(5): 242-249.

[24]Katoh K, Shibata H, Suzuki H, Nara A, Ishidoh K, Kominami E, Yoshimori T, Maki M. The ALG-2-interacting protein Alix associates with CHMP4b, a human homologue of yeast Snf7 that is involved in multivesicular body sorting. Journal of Biological Chemistry, 2003, 278(40): 39104-39113.

[25]Vito P, Pellegrini L, Guiet C, D’Adamio L. Cloning of AIP1, a novel protein that associates with the apoptosis-linked gene ALG-2 in a Ca2+-dependent reaction. Journal of Biological Chemistry, 1999, 274(3): 1533-1540.

[26]Mahul-Mellier A-L, Hemming F J, Blot B, Fraboulet S, Sadoul R. Alix, making a link between apoptosis-linked gene-2, the endosomal sorting complexes required for transport, and neuronal death in vivo. The Journal of Neuroscience, 2006, 26(2): 542-549.

[27]Avila M, Kirchhoff M, Marle N, Hove H D, Chouchane M, Thauvin‐Robinet C, Masurel A, Mosca‐Boidron A L, Callier P, Mugneret F. Delineation of a new chromosome 20q11. 2 duplication syndrome including the ASXL1 gene. American Journal of Medical Genetics Part A, 2013, 161(7): 1594-1598.

[28]Li X, Leo BM, Beck G, Balian G, Anderson G D. Collagen and proteoglycan abnormalities in the GDF-5-deficient mice and molecular changes when treating disk cells with recombinant growth factor. Spine, 2004, 29(20): 2229-2234.