Advanced
Identification of a SNP in Cattle HGD Gene with its Effect on Economic Trait in Hanwoo
Identification of a SNP in Cattle HGD Gene with its Effect on Economic Trait in Hanwoo
Journal of Life Science. 2014. Nov, 24(11): 1168-1173
Copyright © 2014, Korean Society of Life Science
This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : August 08, 2014
  • Accepted : November 17, 2014
  • Published : November 30, 2014
Download
PDF
e-PUB
PubReader
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
About the Authors
정민 한
홍식 공
kebinkhs@empal.com

Abstract
HGD (homogentisate1,2-dioxygenase)유전자는 소의 1번 염색체에 존재하며, tyrosine과 phenylalanine의 이화 작용을 돕는 효소 중 하나로 알려져 있다. 또한 소에서 도체중과 등심단면적에 영향을 주는 유전자로 보고된 바 있다. 본 연구는 한우 집단을 대상으로 HGD유전자내 변이지역을 탐색하고, 경제형질과의 연관성을 분석하기 위하여 실시하였다. 총 14개의 Exon지역을 바탕으로 변이지역을 탐색한 결과, 10개의 SNPs를 확인하였고, 이 중 genotype의 Frequency (MAF>0.1)를 고려하여 6개의 SNPs (G34256T, G34257C, T34284C, T42333G, T42348C, T42468C)를 발굴하여 경제형질과의 연관성을 분석하였다. 그 결과 G34256T에서 근내지방도와 유의적인 연관성이 확인되었고, G34257C에서 등심단면적과 근내지방도에서 유의적인 연관성이 확인되었다. 따라서 본 연구 결과에서 확인된 HGD 유전자의 SNP유전자형은 한우선발 및 개량을 위한 분자육종의 기초 자료로 이용할 수 있을 것으로 사료된다.
Keywords
서 론
소비자들의 육질에 대한 선호도가 높아짐에 따라 우수한 등급의 쇠고기 생산을 증가시키기 위해 가축의 주요 경제형질 개량을 위한 연구가 활발히 진행 중에 있다. 한우 등급 체계는 육량형질(도체중, 등심단면적, 등지방두께)과 육질형질(근내 지방도, 육색 등)로 구분하고 있으며 [2] , 가축의 주요 경제형질 개량을 위하여 경제형질과 관련된 양적형질좌위(quantitative trait locus, QTL)를 탐색하고, 이를 통해 주요 경제형질들과 연관된 후보 유전자 및 단일염기다형성(single nucleotide polymorphism, SNP)을 발굴하여 주요 경제형질과의 연관성 분석 결과를 토대로 우수 개체를 조기에 선발하는 연구가 진행되고 있다 [5] . 최근 유전체 정보를 활용한 연구결과에 따르면 양적형질에 관여하는 주요 유전자가 존재하고 있으며, 다양한 형질들에 대한 유전자 좌위와 유전자들이 확인되었다. 이 중 육질과 육량에 연관된 대표적인 유전자는 APOE (apolipoprotein E), CAST (calpastatin), DGAT1 (diacylglycerol O-acyltransferase 1), Exostosin-1, GHRH (growth hormone-releasing hormone)등의 유전자가 알려져 있으며 이 외에도 육질과 육량에 연관된 다양한 유전자들이 보고되고 있다 [3 , 13 , 16 , 18] .
Homogentisate1,2-dioxygenase (HGD)유전자는 총 42,632 bp로 소의 1번 염색체에 존재하고 있으며, 14개의 exon으로 구성되어 있으며, tyrosine과 phenylalanine의 이화작용을 돕는 효소 중 하나이다. Tyrosine과 phenylalanine은 호르몬 및 신경 전달 물질 합성에 중요한 역할을 하고 있으며, Homogentisate acid는 phenlyalanine과 이화작용을 통해 생명체에 필요한 당으로 분해된다 [15] . HGD유전자는 인간의 경우 변이(c.87+1G>A)에 따라 tyrosine과 phenylalanine의 분해에 결함이 생겨 나타나 Homogentisate acid가 체내에 축적되어 알캅톤뇨증을 유발하는 것으로 알려져 있으나 [17] , 돼지 경우 Accession No. BE232117의 변이지역인 C>T nt354에 따라 도 체형질에 영향을 주는 것으로 보고되었다 [14] .
본 연구는 도체중과 등심단면적과 연관되었다고 보고된 HGD 유전자를 한우에 적용하여 염기서열을 결정하고, SNP를 발굴하여 각 SNP와 경제형질과의 연관성을 구명하고자 수행하였다.
재료 및 방법
- 공시재료 및 Genomic DNA 추출
본 연구를 위해 사용된 sample은 쇠고기이력제를 위해 축산물 품질평가원에서 2010년도에 도축된 한우를 사용하였다. 한우의 HGD 유전자내 SNP탐색을 위해 아비가 서로 다른 24두에서 근육조직을 채취하여 사용하였으며, 경제형질과의 연관성 분석을 위해 도체 및 성장관련 성적(도체중, 등심 단면적, 등지방 두께, 근내 지방도)을 가진 한우 90두의 근육 조직을QuickGene DNA tissue Kit S (KURABO, Japan)를 이용하여 지침서에 따라 수행하였다(한경대학교 동물실험윤리위원회 2014-6). 추출된 DNA는 ND-1000 UV-Vis Spectrophotometer (NanoDrop Technologies, USA)로 순도 및 농도를 측정한 후 본 연구에 이용하였다.
조사된 형질은 축산물 등급판정 세부기준 제 2004-66호의 측정방법을 이용하여 등지방두께, 등심단면적, 근내지방도 및 도체중 항목을 측정하였다( Table 2 ).
Overall means and standard deviation (SD), minimum (Min) and maximum (Max) of traits analyzed in this study
PPT Slide
Lager Image
Overall means and standard deviation (SD), minimum (Min) and maximum (Max) of traits analyzed in this study
- PCR primer의 설계 및 유전자 증폭
HGD 유전자(Accession No. NC_007299)를 구성하고 있는 모든 Exon 영역을 포함할 수 있도록 7개의 primer를 제작하여 각 primer별로 증폭을 위한 조건을 설정하였다( Table 1 ).
PCR primers information for identification of polymorphisms in bovine HGD gene
PPT Slide
Lager Image
PCR primers information for identification of polymorphisms in bovine HGD gene
PCR 증폭을 위한 반응 조건은 genomic DNA 50 ng, Primer 0.5 μM, dNTP 각 2.5 mM 10X Buffer (10 mM Tris-CL, 50 mM KCL, 1.5 mM MgCl2, pH 8.3), Taq polymerase 0.2 μl (10 unit/μl)를 넣고 최종 반응량이 20 μl가 되게 조정한 후 94℃에서 pre denaturation 5분, 94℃에서 denaturation 30초, primer에 적당한 annealing온도에서 30초, Pre Extension 72℃에서 30초로 총 35 cycle로 실행 후 final Extension을 72℃에서 5분 동안 GeneAmp PCR system 9700 Aplied Biosystems, USA)로 수행하여 증폭된 산물은 EtBr이 포함된 2% agarose gel을 이용하여 증폭여부를 확인하였다.
- SNP 탐색 및 유전자형 결정
HGD 유전자내 SNP를 탐색하기 위하여 증폭된 PCR product를 purification 과정을 거친 후, BigDye ® Terminator v3.1 Cycle Sequencing Kits (Aplied Biosystems, USA)를 활용하여 각각의 염기에 형광 dye를 부착한 후, ABI 3130XL Genetic Analyzer (Applied Biosystems, USA)를 통하여 개체별 염기 서열을 결정하였다. 결정된 염기서열은 seqMAN II (DNA STAR Inc. USA)프로그램을 이용하여 SNP를 탐색하였다.
- Haplotype 및 LD block 분석
HGD유전자의 변이간의 연관불균형(Linkage disequilibrium) 정도를 추정하기 위해 HaploView software package [1] 를 활용하여 분석을 수행하여 D 과 r 2 값을 추정하였다.
- 통계 분석
HGD유전자의 변이지역에 대한 경제형질과의 효과를 추정하기 위해 SAS 9.2 package/PC software (SAS, USA)를 통해 분석을 실시하였다.
  • Yijk=μ+Gj+eijk
  • 상기모형에서,
  • Yijk= 대상형질에 대한 관측치
  • μ= 대상형질의 전체 평균
  • Gj= 유전자형(Genotype) 효과
  • eijk= 임의 오차
결과 및 고찰
- HGD 유전자내 SNP탐색
NCBI에 등록된 염기서열(Accession No. NC_007299)의 Exon지역을 바탕으로 7개의 Primer를 제작하여 혈연관계가 없는 한우 24두를 대상으로 염기서열 분석을 통해 SNP지역을 확인하였다. 분석한 결과 총 10개의 SNPs (T11187C, T11301A, T11398G, G29833A, G34256T, G34257C, T34284C, T42333G, T42348C, T42468C)가 확인되었다( Fig. 1 ). 발견된 10개의 SNPs 중 T11187C (rs43239703), G29833A (rs109637397), T42348C (rs137531856), T42468C (rs133513769)은 NCBI SNP DB (http://www.ncbi.nlm.nih.gov/SNP/snp_ref.cgi?locusId=507496 )에 등록되어 있었으며 6개는 신규 SNPs 지역으로 확인되었다 ( Table 3 ). Zhou 등(2010) [20] 에 의해 보고된 9개의 변이지역 중 A219G, A282C, A11206G, A29614G, G29859A, C34233, A42460G은 한우에서는 발견되지 않았다. 그리고 Exon 10번의 SNP (G29833A)는 Arginine이 Histidine으로 바뀌는 nonsynonymous mutation으로 확인되었다.
Genotype and allele frequencies of the identified HGD gene in Hanwoo (24 samples)
PPT Slide
Lager Image
MAF : Minor Allele Frequency, AA: Amino Acid Change
PPT Slide
Lager Image
SNP map of Hanwoo HGD gene. The coding exon is marked by black, and 5 and 3 UTRs are white blocks.
- HGD 유전자형 분석
발굴된 10개의 SNP 중 genotype의 Frequency (MAF > 0.1)를 고려하여 6개 지역(G34256T, G34257C, T34284C, T42333G, T42348C, T42468C)을 선별하여 90두의 개체별 유전자형을 확인하였다. 유전자형의 빈도는 G34256T는 GG유전자형(0.5843), GT유전자형(0.3708) 그리고 TT유전자형(0.0449)이 확인되었다. G34257C는 GG유전자형(0.3596), GC유전자형(0.3708) 그리고 CC유전자형(0.2697)으로 확인되었다. T34284C는 TT유전자형(0.1011), TC유전자형(0.5056) 그리고 CC유전자형(0.3933)으로 확인되었다. T42333G는 TT유전자형(0.3708), TG유전자형(0.5281) 그리고 GG유전자형(0.1011)으로 확인되었다. T42348C는 TT유전자형(0.1000), TC유전자형(0.5444) 그리고 CC유전자형(0.3556)으로 확인되었다. T42468C는 TT유전자형(0.1222), TC유전자형(0.5556), CC유전자형(0.3222)으로 확인되었다( Table 4 ).
Genotype and allele frequencies of the identified HGD gene in Hanwoo (90 samples)
PPT Slide
Lager Image
Genotype and allele frequencies of the identified HGD gene in Hanwoo (90 samples)
- LD 분석 및 경제형질과 연관성 분석
6개의 SNPs를 대상으로 연관불균형(Linkage disequilibrium, LD) 분석 및 경제형질과의 연관성 분석을 수행하여 Table 5 에 제시하였다. T42333G, T42348C 그리고 T42468C 변이들 간의 D 값의 경우 0.784에서 1로 나타났으며, r-square 값은 0.819에서 0.952로 확인되어 강한 LD가 형성되어 있음이 확인되었다( Table 5-a ; Fig. 2 ). 그리고 Haplotype은 총8개로 확인되었으며, 빈도수가 0.01이상인4개(-ACC-, -GTT-, -ACT-, -ATT-)의 haplotype빈도를 Table5-b 에 나타내었다. HGD유전자의 Ht 1번(-ACC-)의 빈도는 0.595, Ht 2번(-GTT-)의 빈도는 0.354로 확인되었다. Ht 3번(-ACT-)의 빈도는 0.034, Ht 4번(-ATT-)의 빈도는 0.01로 매우 낮게 확인되었다.
Linkage disequilibrium and haplotype structure in the HGD gene
PPT Slide
Lager Image
Linkage disequilibrium and haplotype structure in the HGD gene
PPT Slide
Lager Image
Linkage disequilibrium plots with the single nucleotide polymorphism within the HGD of this study.
또한 각각의 haplotype과 경제형질간의 연관성 분석한 결과 통계적으로 유의적이지 못한 것( p >0.05)으로 확인되었다.
- 경제형질과의 연관성 분석
90두를 대상으로 6개의 변이지역에 대한 개체별 유전자형 분석결과를 바탕으로 유전자형과 경제형질간의 연관성 분석을 실시하였다. G34256T지역과 경제형질간의 연관성 분석 결과 근내지방도에서 유의적인 연관성이 검출되었다. GG(3.82±0.27), TT(3.25±0.98) 유전자형에 비해 GT(5.09±0.344)유전자형이 우수한 근내지방도를 가지고 있는 것으로 확인되었다 ( p <0.05). 또한 G34257C에서는 등심단면적과 근내지방도에서 유의적인 연관성이 검출되었다( Table 6 ). GG(79.96±1.59)유전자형이 CC(75.83±1.83), GC(74.54±1.56) 유전자형에 비해 우 수한 등심단면적을 가지는 것으로 확인되었으며, GC(5.09±0.34), GG(4.03±0.48)유전자형이 CC(3.45±0.40)유전자형 보다 우수한 근내지방도를 가지는 것으로 확인되었다( p <0.05).
Least squares means and standard errors for economic traits of genotypes in Hanwoo HGD gene
PPT Slide
Lager Image
BF : Backfat Thickness, LMA : Longissimus Muscle Area, CW : Carcass Weight, MS : Marbling Score
소의 1번 염색체는 약 158,337,000 bp로 육질등급과 육량등급에 영향을 주는 도체중(61.852~79.15 cM), 근내지방도(51.01~76.50 cM), 등심단면적(65.30 cM) 등 다양한 QTL을 가지고 있다 [6 , 10 , 11] . HGD 유전자는 58.5cM에 위치하고 있어 한우의 등급판정에 중요한 요소인 도체중, 근내지방도, 등심단면적의 QTL영역 내에 위치하고 있다 [20] .
본 연구 결과 한우에서HGD유전자는 총 6개의 변이지역이 확인되었으며. 9개의 SNPs가 확인된 Chinese red cattle 종과는 1개의 SNP (T42333G)가 공통적으로 나타나는 것을 확인하였다. 또한 Ponsuksili 등(2005) [14] 과 Zhou 등(2010) [20] 에 따르면 Chinese red cattle과 돼지에서 도체중과 관련된 후보 유전자로 보고되었으나 본 연구 결과 한우에서는 등심단면적과 근내지방도 형질과의 연관성이 나타나는 것을 확인할 수 있었다.
따라서 본 연구 결과에서 확인된 근내지방도와 등심단면적과 연관된 한우의 HGD 유전자 내 변이지역인 G34256T, G34257C는 한우의 육질과 관련된 주요 형질에 영향을 미칠것으로 사료된다. 또한 추후 HGD유전자를 포함하고 있는 QTL 영역 내의 변이지역을 확인한다면 한우에서 성장과 지방대사 작용에 관련된 후보유전자 연구에 활용 가능할 것으로 사료된다.
References
Barrett J. C. , Fry B. , Daly M. J. 2005 Haploview : analysis and visualization of LD and haplotype maps Bioinformatics 21 263 - 265    DOI : 10.1093/bioinformatics/bth457
Cahyadi M. , Maharani D. , Ryoo S. H. , Lee S. H. , Lee J. H. 2012 Identification of SNPs in TG and EDG1 genes and their relationships with carcass traits in Korean cattle (Hanwoo) CNU J Agric Sci 39 (3) 349 - 355    DOI : 10.7744/cnujas.2012.39.3.349
Cheong H. S. , Yoon D. H. , Kim L. H. , Park B. L. , Choi Y. H. , Chung E. R. , Cho Y. M. , Park E. W. , Cheong I. C. , Oh S. J. , Yi S. G. , Park T. , Shin H. D. 2006 Growth hormone-releasing hormone (GHRH) polymorphisms associated with carcass traits of meat in Korean cattle BMC Genetics 7 35 -
Comings D. E. , MacMurray J. P. 2000 Molecular heterosis : a review Mol Genet Metab 71 19 - 31    DOI : 10.1006/mgme.2000.3015
Kim B. S. , Kim N. K. , Lee S. , Cho Y. M. , Heo H. N. , Park E. W. , Yang B. K. , Yoon D. 2011 Association study between the polymorphisms of exostosin-1 gene and economic traits in Hanwoo J Anim Sci Technol 53 (1) 7 - 13    DOI : 10.5187/JAST.2011.53.1.7
Kim J. J. , Farnir F. , Savell J. , Taylor J. F. 2003 Detection of quantitative trait loci for growth and beef carcass fatness traits in a cross between bos taurus (Angus) and bos indicus (Brahman) cattle J Anim Sci 81 1933 - 1942
Laschi M. , Tinti L. , Braconi D. , Millucci L. , Ghezzi L. , Amato L. , Selvi E. , Spreafico A. , Bernardini G. , Santucci A. 2013 Homogentisate 1,2 dioxygenase is expressed in human osteoarticular cells: implications in alkaptonuria J Cell Physiol 227 (9) 3254 - 3257
Lee S. H. , Cho Y. M. , Lee S. H. , Kim B. S. , Kim N. K. , Choy Y. H. , Kim K. H. , Yoon D. , Im S. K. , Oh S. J. , Park E. W. 2008 Identification of marbling-related candidate genes in M. longissimus dorsi of high- and low marbled Hanwoo (Korean Native Cattle) steers BMB Reports 41 (12) 846 - 851    DOI : 10.5483/BMBRep.2008.41.12.846
Lee S. H. , Lim D. , Dang C. G. , Chang S. S. , Kim H. C. , Jeon G. J. , Yeon S. H. , Jang G. W. , Park E. W. , Oh J. D. , Lee H. K. , Lee J. H. , Kang H. S. , Yoon D. 2013 Genome wide association test to identity QTL for dressing percentage in Hanwoo CNU Journal of Agricultural Science 40 (2) 155 - 162    DOI : 10.7744/cnujas.2013.40.2.155
McClure M. C. , Morsci N. S. , Schnabel R. D. , Kim J. W. , Yao P. , Rolf M. M. , McKay S. D. , Gregg S. J. , Chapple R. H. , Northcutt S. L. , Taylor J. F. 2010 A genome scan for quantitative trait loci influencing carcass, post-natal growth and reproductive traits in commercial Angus cattle Anim Genet 41 (6) 597 - 607    DOI : 10.1111/j.1365-2052.2010.02063.x
Nalaila S. M. , Stothard P. , Moore S. S. , Li C. , Wang Z. J. 2012 Whole-genome QTL scan for ultrasound and carcass merit traits in beef cattle using Bayesian shrinkage method J Anim Breed Genet 129 (2) 107 - 119    DOI : 10.1111/j.1439-0388.2011.00954.x
Oh J. D. , Lee K. W. , Jung I. J. , Jeon G. J. , Lee H. K. , Kong H. S. 2011 Association of a single nucleotide polymorphism with economic traits in porcine uncoupling protein 3 gene J Life Sci 21 155 - 158    DOI : 10.5352/JLS.2011.21.1.155
Oh J. D. , Lee J. A. , Lee K. W. , Park D. K. , Cho B. W. , Jeon G. J. , Lee H. K. , Kong H. S. 2010 Identification of polymorphisms in CAST gene associated with economic traits in Hanwoo (Bos taurus coreanae) J Life Sci 20 1498 - 1504    DOI : 10.5352/JLS.2010.20.10.1498
Ponsuksili S. , Chomdej S. , Murani E. , Bläser U. , Schreinemachers H. J. , Schellander K. , Wimmers K. 2005 SNP detection and genetic mapping of porcine genes encoding enzymes in hepatic metabolic pathways and evaluation of linkage with carcass traits Anim Genet 36 (6) 477 - 483
Rodríguez J. M. , Timm D. E. , Titus G. P. , Beltrán-Valero De Bernabé D. , Criado O. , Mueller H. A. , Rodríguez De Córdoba S. , Peñalva M. A. 2000 Structural and functional analysis of mutations in alkaptonuria Hum Mol Genet 9 (15) 2341 - 2350    DOI : 10.1093/oxfordjournals.hmg.a018927
Shin K. H. , Shin S. C. , Chung K. Y. , Chung E. R. 2009 Effects of SNP markers of the apolipoprotein E (APOE) gene on meat quantity and quality traits in Korea cattle Korean J Food Sci Ani Resour 29 (1) 108 - 113    DOI : 10.5851/kosfa.2009.29.1.108
Sakthivel S. , Zatkova A. , Nemethova M. , Surovy M. , Kodasi L. , Saravanan S. P. 2014 Mutation screening of the HGD gene identifies a novel alkaptonuria mutation with significant founder effect and high prevalence Ann Hum Genet 78 155 - 164    DOI : 10.1111/ahg.12055
Thaller G. , Kramer W. , Winter A. , Kaupe B. , Erhardt G. , Fries R. 2003 Effects of DGAT1 variants on milk production traits in German cattle breeds J Anim Sci 81 1911 - 1918
Zhou G. L. , Cao Y. , Li M. , Zhang L. C. , Yu Y. S. , Jin H. G. 2010 Meat quality and carcass traits in relation to HGD-BstXI and HGD-HaeIII PCR-RFLP polymorphism in Chinese red cattle Meat Sci 85 (2) 270 - 273    DOI : 10.1016/j.meatsci.2010.01.011
Zhou G. , Dudgeon C. , Li M. , Cao Y. , Zhang L. , Jin H. 2010 Molecular cloning of the HGD gene and association of SNPs with meat quality traits in Chinese red cattle Mol Biol Rep 37 (1) 603 - 611    DOI : 10.1007/s11033-009-9860-4