家畜全基因组关联分析研究进展

doi:10.16431/j.cnki.1671-7236.2018.02.022

摘要/Abstract

摘要：

全基因组关联分析（genome-wide association study，GWAS）是一种研究经济性状候选基因的分析方法。近年来，随着家畜全基因组测序的完成，大量的单核苷酸多态性（single nucleotide polymorphisms，SNPs）被标识，GWAS也越来越多地应用于家畜重要性状的研究领域中，在动物遗传育种中，通过对家畜基因组进行GWAS分析研究，找到控制家畜主要经济性状的重要SNPs，从而挖掘重要经济性状的候选基因。作者详细综述了GWAS的分析方法及其在重要家畜育种中的研究进展。GWAS分析方法包括基因组控制法（genommic control）、分层分析法（stratification analysis）、主成分分析法（principal components analysis，PAC）和混合线性模型分析法（mixed-linear-model association，MLMA），通路分析方法包括非核算法（基因功能富集分析（gene set enrichment analysis，GSEA）和分层贝叶斯优取（hierarchical Bayes prioritization，HBP））和核算法。依据不同的目标性状选择合理的分析方法，提高GWAS分析结果的准确性，为进一步利用GWAS分析各种性状的遗传基础提供合理的借鉴。

关键词: 家畜; 全基因组关联分析(GWAS); 统计方法; 通路分析

Abstract:

Genome-wide association study (GWAS) is an analytical method, has been used to study candidate of economical traits in animals. In recent years, with the completion of the whole genome sequencing for livestock, a large number of single nucleotide polymorphisms (SNPs) have been identified. GWAS is also more and more applied in the field of important traits for livestock. Finding important SNPs that control the important economic traits of livestock by GWAS analysis to unearth candidate genes for important economic traits in animal genetics and breeding. This paper reviews the GWAS analysis methods and it's research progress in important livestock breeding. Statistical analysis methods include genommic control, stratification analysis, principal components analysis (PAC) and mixed-linear-model association (MLMA). The methods of pathway analysis include non-nuclear methods (gene enrichment analysis (GSEA) and hierarchical Bayes prioritization (HBP)) and algorithms. This review will provide reference for further research on genetic background of the important traits of livestock by GWAS.

Key words: livestock; genome-wide association study (GWAS); statistical methods; pathway analysis

中图分类号:

S813.3

赵德胜, 张崇志. 家畜全基因组关联分析研究进展[J]. 《中国畜牧兽医》, 2018, 45(2): 463-470.

ZHAO Desheng, ZHANG Chongzhi. Research Progress on Genome-wide Association in Livestock[J]. , 2018, 45(2): 463-470.

参考文献

[1] HIRSCHHORN J N, DALY M J. Genome-wide association studies for common diseases and complex traits[J].Nature Reviews Genetics, 2005, 6(2):95-108.
[2] RISCH N, MERIKANGAS K. The future of genetic studies of complex human diseases[J].Science, 1996, 273(5281):1516-1517.
[3] KLEIN R J, ZEISS C, CHEW E Y, et al. Complement factor H polymorphism in age-related macular degeneration[J].Science, 2005, 308(5720):385-389.
[4] 严卫丽. 复杂疾病全基因组关联研究进展——遗传统计分析[J]. 遗传, 2008, 30(5):543-549.YAN W L. Genome-wide association study on complex diseases:Genetic statistical issues[J]. Heredias, 2008, 30(5):543-549.(in Chinese)
[5] SPIELMAN R S, MCGINNIS R E, EWENS W J. Transmission test for linkage disequilibrium:The insulin gene region and insulin-dependent diabetes mellitus (IDDM)[J]. American Journal of Human Genetics, 1993, 52(3):506.
[6] 郝兴杰, 胡林, 张淑君. 全基因组关联分析方法的研究进展[J]. 畜牧兽医学报, 2016, 47(2):213-217.HAO S J, HU L, ZHANG S J. Progresse in research of genome-wide association study method[J]. Acta Veterinaria et Zootechnica Sinica, 2016, 47(2):213-217.(in Chinese)
[7] DEVLIN B, ROEDER K, WASSERMAN L. Genomic control, a new approach to genetic-based association studies[J].Theoretical Population Biology, 2001, 60(3):155.
[8] PRITCHARD J K, ROSENBERG N A. Use of unlinked genetic markers to detect population stratification in association studies[J].The American Journal of Human Genetics, 1999, 65(1):220-228.
[9] PRITCHARD J K, STEPHENS M, ROSENBERG N A, et al. Association mapping in structured populations[J].The American Journal of Human Genetics, 2000, 67(1):170-181.
[10] CAPELLI C, REDHEAD N, ROMANO V, et al. Population structure in the Mediterranean basin:AY chromosome perspective[J]. Annals of Human Genetics, 2006, 70(2):207-225.
[11] PROGRAM A L. Population structure and eigenanalysis[J].PLoS Genetics, 2006, 2(12):e190.
[12] PRICE A L, PATTERSON N J, PLENGE R M, et al. Principal components analysis corrects for stratification in genome-wide association studies[J]. Nature Genetics, 2006, 38(8):904-909.
[13] ZHANG Z, BUCKLER E S, CASSTEVENS T M, et al. Software engineering the mixed model for genome-wide association studies on large samples[J].Briefings in Bioinformatics,2009,10(6):664-675.
[14] LAM A C, MICHAEL S, AULCHENKO Y S, et al. Rapid and robust association mapping of expression quantitative trait loci[J]. BMC Proceedings, 2007, 1(1):S144.
[15] SACCONE S F, SACCONE N L, SWAN G E, et al. Systematic biological prioritization after a genome-wide association study:An application to nicotine dependence[J].Bioinformatics, 2008, 24(16):1805.
[16] ELEFTHEROHORINOU H, WRIGHT V, HOGGART C, et al. Pathway analysis of GWAS provides new insights into genetic susceptibility to 3 inflammatory diseases[J].PLoS One, 2009, 4(11):e8068.
[17] JIA P, WANG L, MELTZER H Y, et al. Common variants conferring risk of schizophrenia:A pathway analysis of GWAS data[J]. Schizophrenia Research, 2010, 122(1-3):38.
[18] LEWINGER J P, CONTI D V, BAURLEY J W, et al. Hierarchical Bayes prioritization of marker associations from a genome-wide association scan for further investigation[J].Genetic Epidemiology, 2007, 31(8):871-882.
[19] SUBRAMANIAN A, TAMAYO P, MOOTHA V K, et al. Gene set enrichment analysis:A knowledge-based approach for interpreting genome-wide expression profiles[J]. Proceedings of the National Academy of Sciences, 2005, 102(43):15545-15550.
[20] KULESHOV M V, JONES M R, ROUILLARD A D, et al. Enrichr:A comprehensive gene set enrichment analysis web server 2016 update[J]. Nucleic Acids Research, 2016, 44(W1):W90.
[21] TAMAYO P, STEINHARDT G, LIBERZON A, et al. The limitations of simple gene set enrichment analysis assuming gene independence[J].Statistical Methods in Medical Research, 2016, 25(1):472.
[22] XU Z, PAN W. Binomial mixture model based association testing to account for genetic heterogeneity for GWAS[J].Genetic Epidemiology, 2016, 40(3):202-209.
[23] SCHAID D J. Genomic similarity and kernel methods Ⅱ:Methods for genomic information[J]. Human Heredity, 2010, 70(2):132-140.
[24] LARSON N B, SCHAID D J. A kernel regression approach to gene-gene interaction detection for case-control studies[J].Genetic Epidemiology, 2013, 37(7):695.
[25] DAETWYLER H D, SCHENKEL F S, SARGOLZAEI M, et al. A genome scan to detect quantitative trait loci for economically important traits in Holstein cattle using two methods and a dense single nucleotide polymorphism map[J].Journal of Dairy Science, 2008, 91(8):3225-3236.
[26] SCHOPEN G C B, VISKER M, KOKS P D, et al. Whole-genome association study for milk protein composition in dairy cattle[J].Journal of Dairy Science, 2011, 94(6):3148-3158.
[27] HUANG W, KIRKPATRICK B W, ROSA G J M, et al. A genome-wide association study using selective DNA pooling identifies candidate markers for fertility in Holstein cattle[J]. Animal Genetics, 2010, 41(6):570-578.
[28] SODELAND M, KENT M P, OLSEN H G, et al. Quantitative trait loci for clinical mastitis on chromosomes 2, 6, 14 and 20 in Norwegian Red cattle[J].Animal Genetics, 2011, 42(5):457.
[29] CRISPIM A C, KELLY M J, GUIMARÃES S E, et al. Multi-trait GWAS and new candidate genes annotation for growth curve parameters in Brahman cattle[J]. PLoS One, 2015, 10(10):e0139906.
[30] BOICHARD D, SANCHEZ M P, BARBAT A, et al. Identification of candidate causal variants underlying QTL in dairy cattle through GWAS and Bayesian approach at the sequence level[A].PAG XXIV-Plant and Animal Genome Conference[C]. 2016.
[31] VAN DEN BERG I, BOICHARD D, LUND M S. Multi-breed GWAS and meta-analysis using sequences of five dairy cattle breeds improve QTL mapping[A].Annual Meeting of the European Association for Animal Production[C]. 2015.
[32] KIJAS J W, TOWNLEY D, DALRYMPLE B P, et al. A genome wide survey of SNP variation reveals the genetic structure of sheep breeds[J].PLoS One, 2009, 4(3):e4668.
[33] TOSSER-KLOPP G, BARDOU P, BOUCHEZ O, et al. Design and characterization of a 52K SNP chip for goats[J].PLoS One, 2014, 9(1):e86227.
[34] BECKER D, TETENS J, BRUNNER A, et al. Microphthalmia in Texel sheep is associated with a missense mutation in the paired-like homeodomain 3(PITX3) gene[J].PLoS One, 2010, 5(1):e8689.
[35] DONG Y, XIE M, JIANG Y, et al. Sequencing and automated whole-genome optical mapping of the genome of a domestic goat (Capra hircus)[J]. Nature Biotechnology, 2013, 31(2):135-141.
[36] JIANG Y, XIE M, CHEN W, et al. The sheep genome illuminates biology of the rumen and lipid metabolism[J].Science, 2014, 344(6188):1168-1173.
[37] GREYVENSTEIN O F, REICH C M, VAN MARLE-KOSTER E, et al. Polyceraty (multi-horns) in Damara sheep maps to ovine chromosome 2[J]. Animal Genetics, 2016, 47(2):263-266.
[38] REN X, YANG G L, PENG W F, et al.Erratum:A genome-wide association study identifies a genomic region for the polycerate phenotype in sheep (Ovis aries)[J]. Scientific Reports, 2016, 7:25322.
[39] LI M H, TⅡRIKKA T, KANTANEN J. A genome-wide scan study identifies a single nucleotide substitution in ASIP associated with white versus non-white coat-colour variation in sheep (Ovis aries)[J]. Heredity, 2014, 112(2):122-131.
[40] MARTIN P M, PALHIÈRE I, RICARD A, et al. Genome wide association study identifies new loci associated with undesired coat color phenotypes in Saanen goats[J]. PLoS One, 2016, 11(3):e0152426.
[41] GHOLIZADEH M, RAHIMIMIANJI G, NEJATIJAVAREMI A. Genomewide association study of body weight traits in Baluchi sheep[J].Journal of Genetics, 2015, 94(1):143.
[42] ALMAMUN H A, KWAN P, CLARK S A, et al. Genome-wide association study of body weight in Australian Merino sheep reveals an orthologous region on OAR6 to human and bovine genomic regions affecting height and weight[J].Genetics Selection Evolution, 2015, 47(1):66.
[43] MATIKA O, RIGGIO V, ANSELMEMOIZAN M, et al. Genome-wide association reveals QTL for growth, bone andin vivo carcass traits as assessed by computed tomography in Scottish Blackface lambs[J]. Genetics Selection Evolution, 2016, 48(1):11.
[44] ELBELTAGY A R. Heat stress in desert sheep and goats:Signatures of natural selection and initiation of GWAS[A].Plant and Animal Genome XXⅡ Conference[C]. 2014.
[45] CHO I C, ZHONG T, SEO B Y, et al. Whole-genome association study for the roan coat color in an intercrossed pig population between Landrace and Korean native pig[J]. Genes & Genomics, 2011, 33(1):17-23.
[46] DUIJVESTEIJN N, KNOL E F, MERKS J W, et al. A genome-wide association study on androstenone levels in pigs reveals a cluster of candidate genes on chromosome 6[J].BMC Genetics, 2010, 11(1):1-11.
[47] ZHANG C, BRUCE H, YANG T, et al. Genome wide association studies (GWAS) identify QTL on SSC2 and SSC17 affecting loin peak shear force in crossbred commercial pigs[J].PLoS One, 2016, 11(2):e0145082.
[48] VERARDO L L, SILVA F F, VARONA L, et al. Bayesian GWAS and network analysis revealed new candidate genes for number of teats in pigs[J]. Journal of Applied Genetics, 2015, 56(1):123-132.
[49] TANG J, ZHANG Z, YANG B, et al. Identification of loci affecting teat number by genome-wide association studies on three pig populations[J].Asian-Australasian Journal of Animal Sciences, 2017, 30(1):1-7.
[50] WANG X, LIU X, DENG D, et al. Genetic determinants of pig birth weight variability[J].BMC Genetics, 2016, 17(Suppl 1):15.
[51] GUO X, SU G, CHRISTENSEN O F, et al. Genome-wide association analyses using a Bayesian approach for litter size and piglet mortality in Danish Landrace and Yorkshire pigs[J].BMC Genomics, 2016, 17(1):1-12.
[52] 李学武, 王瑞军, 王志英,等. 绒山羊绒毛品质遗传研究进展[J]. 中国畜牧兽医, 2016, 43(12):3275-3284.LI X W, WANG R J, WANG Z Y, et al. The genetic research progress of fiber and wool quality in cashmere goat[J]. China Animal Husbandry & Veterinary Medicine, 2016, 2016, 43(12):3275-3284.(in Chinese)