中国农业科学 ›› 2015, Vol. 48 ›› Issue (1): 140-153.doi: 10.3864/j.issn.0578-1752.2015.01.14

• 畜牧·兽医 • 上一篇    下一篇

两个贵州山羊品种GHSR和GHRL基因遗传变异及其与生长性状的关联性

宋桃伟,蔡惠芬,罗卫星,刘若余,张依裕,孙岩岩,刘彬   

  1. 贵州大学动物科学学院/高原山地动物遗传育种与繁殖教育部重点实验室,贵阳 550025
  • 收稿日期:2014-01-23 出版日期:2015-01-01 发布日期:2015-01-01
  • 通讯作者: 罗卫星,E-mail:ef.wxluo@gzu.edu.cn
  • 作者简介:宋桃伟,E-mail:songtw400407@126.com
  • 基金资助:
    贵州省农业科技攻关项目(黔科合NY字(2012)3006号)

Association of GHSR and GHRL Gene Genetic Variation with Growth Traits in Two Guizhou Goat Breeds

SONG Tao-wei, CAI Hui-fen, LUO Wei-xing, LIU Ruo-yu, ZHANG Yi-yu, SUN Yan-yan, LIU Bin   

  1. Key Laboratory of Animal Genetics, Breeding and Reproduction in the Plateau Mountainous Region, Ministry of Education/College of Animal Science, Guizhou University , Guiyang 550025
  • Received:2014-01-23 Online:2015-01-01 Published:2015-01-01

摘要: 【目的】探讨GHSR和GHRL基因作为山羊体重体尺性状候选基因的可能性,寻找与山羊生长性状相关的分子遗传标记。【方法】以贵州白山羊和贵州黑山羊为研究素材,采用DNA混合池结合PCR产物直接测序及PCR-SSCP技术检测GHSR和GHRL基因的单核苷酸多态性,利用SPSS (18.0)软件GLM统计模型分析基因SNPs不同基因型及合并基因型与贵州山羊生长性状的关联性,同时采用在线软件对GHSR基因进行生物信息学分析。【结果】在GHSR基因外显子2和3′UTR分别检测到G200A同义突变和T628C位点,而在5′UTR区和外显子1则未发现多态性。设计一对特异性引物对G200A位点进行PCR-SSCP分析,表现为3种基因型,依次命名为GG、GA和AA。GG基因型为优势基因型,在贵州白山羊和贵州黑山羊母羊群体中等位基因G为优势等位基因,其等位基因频率分别为0.6731和0.5243。而在贵州黑山羊公羊群体中A则为优势等位基因,其频率为0.5122。多态信息含量均表现为中度多态(0.25<PIC<0.50)。在GHRL基因内含子2处发现1个G141A突变,外显子2和外显子4处分别检测到2个同义突变(T78C和C14T),设计一对特异性引物对C14T位点进行SSCP分析,显示为2种基因型CT和CC。在所有试验群体中,CC基因型为优势基因型,其频率分别为0.7692、0.9417和0.9390,等位基因C为优势等位基因,其频率依次为0.8846、0.9709和0.9695,多态信息含量均显示为低度多态(PIC<0.25)。χ2检验显示所有试验群体G200A和C14T位点基因型分布均符合Hardy-Weinberg平衡(P>0.05)。关联分析表明,GHSR基因G200A位点与山羊体重、体高、体斜长和管围显著相关 (P<0.05),纯合子GG基因型优势较为明显,贵州白山羊GG基因型个体的体高显著高于GA基因型个体和管围显著高于AA基因型个体(P<0.05)。C14T位点贵州黑山羊 (公羊) CC基因型个体的体重、体高和胸围显著高于CT基因型 (P<0.05)。组合基因型对体高和胸围指标的影响显著(P<0.05),表现为CCGG合并基因型个体的体高显著高于CCAA合并基因型个体 (P<0.05)。生物信息学分析揭露,GHSR基因5’UTR处未检测到核心启动子区和CpG岛,仅发现1个CAAT-box和部分潜在的转录因子结合位点,G200A导致mRNA二级结构发生明显变化。GHSR蛋白二级结构以α-螺旋为主 (48.90%),三级结构及跨膜螺旋结构预测表明该蛋白是膜蛋白。【结论】研究初步推测GHSR基因和GHRL基因多态性及合并基因型对山羊生长性状具有较显著的影响,G200A和C14T位点可作为山羊生长性状的有效遗传标记用于指导山羊的分子育种。

关键词: GHSR基因, GHRL基因, 遗传变异, 生长性状, 关联分析

Abstract: 【Objective】 The aim of this study was to explore the relationships between GHSR and GHRL gene and body weight & size traits, and search the molecular genetic markers related to goat growth traits.【Method】Guizhou white goat and Guizhou black goat were chosen as subjects, the DNA pooling was constructed and the technology of direct sequencing of PCR products and PCR-SSCP were used to detect the single nucleotide polymorphism of GHSR and GHRL genes, and its association of polymorphisms and different genotypes and combined genotypes with growth traits was analyzed by SPSS(18.0). Meanwhile, the bioinformatics analysis of GHSR gene was conducted by online software.【Result】The result showed that two SNPs (G200A and T628C) were detected in exon 2 and 3′UTR but no SNPs in exon1 and 5′UTR of GHSR gene, respectively. G200A was a sense mutation site and devided into three genotypes GG, GA and AA by PCR-SSCP. Genotype GG was the dominant genotype and allele G was dominant allele with the frequencies of 0.6731 and 0.5243 in Guizhou white goat and Guizhou black goat (♀) groups instead of allele A in Guizhou black goat (♂). All populations were moderate polymorphic in G200A site (0.25<PIC<0.50). G141A was found in intron 2 and two sense mutation sites (T78C and C14T) were identified in exon2 and exon4 of GHRL gene, respectively. There were two genotypes (CT and CC) in C14T site, genotype CC was the dominant genotype, which the frequency were 0.7692, 0.9417 and 0.9390, allele C was dominant allele with the frequency of 0.8846, 0.9709 and 0.9695, respectively. These populations were low polymorphic (PIC<0.25). χ2 test indicated that SNPs(G200A and C14T) fit with Hardy-Weinberg equilibrium in the population (P>0.05). Correlation analysis revealed that G200A site was significantly associated with body weight, body height, body length and hucklebone width (P<0.05), GG genotype was superior genotype. Additionally, body weight, body height and chest girth of individuals of Guizhou black goat (♂)with genotype CC was significantly higher than those with genotype CT(P<0.05) in C14T site. Meanwhile, the combined genotype was significantly associated with body height and chest girth, the individuals with CCGG genotype had significantly better body height than that of genotype CCAA (P<0.05). Bioinformatics analysis indicated that no core promoter regions and CpG island were found but a CAAT-box and several important transcription factors in 5′UTR region of GHSR gene. G200A had led to the change of secondary structure of mRNA. The secondary structure of protein was mainly α-helix (48.90%). Tertiary structure and transmembrane helix prediction showed that GHSR protein was a membrane protein.【Conclusion】The results revealed that SNPs and combined genotype of GHSR and GHRL gene had effect on growth traits in goat, G200A and C14T sites could be used as effective genetic markers for goat molecular breeding.

Key words: GHSR gene, GHRL gene, genetic variation, growth traits, correlation analysis