湖羊PLAG1基因5'调控区多态性及其与早期体重的关联分析

doi:10.11843/j.issn.0366-6964.2021.02.006

摘要/Abstract

摘要： 旨在克隆湖羊PLAG1基因5'调控区序列，明确PLAG1基因5'调控区多态性与湖羊早期体重的关系，寻找用于湖羊生长性状辅助选择的分子标记。本研究利用5'RACE技术鉴定PLAG1基因转录起始位点，以456只断奶湖羊为对象，利用测序法筛选PLAG1 5'调控区SNP位点，使用SPSS 18.0软件分析不同基因型对湖羊初生体重和断奶体重的影响。结果表明，PLAG1基因转录起始位点均位于g.30535位点，测序发现，PLAG1基因5'调控序列存在15个SNPs位点，在湖羊群体中均有3种基因型，均处于Hardy-Weinberg平衡状态。关联分析发现，g.28888A>T位点AA型、g.28901G>T位点GG型、g.30802C>T位点CC型、g.30825A>G位点AA型的初生体重均显著大于相应的杂合基因型（P<0.05），g.30737T>C位点CC型断奶体重显著小于TC型（P<0.05）。连锁分析发现，g.28888A>T、g.28901G>T、g.28913T>A、g.29717G>C、g.29724A>C、g.29725G>A、g.29768T>C、g.29771A>G、g.30141T>A、g.30144A>G、g.30691C>A、g.30694A>T、g.30737T>C、g.30802C>T、g.30825A>G，15个位点连锁且有18种单倍型。将样本量较多的5种单倍型与湖羊早期体重进行关联分析发现，AAGGTAGGAAGGTTAATTAACCAATCCCAA单倍型初生体重显著大于ATGTTAGCACGATCAGTAAGCAATTCCTAG单倍型（P<0.05），断奶体重显著大于AAGGAAGGAAGGTTAATTAACCAACCCCAA单倍型（P<0.05）。上述结果表明，湖羊群体中PLAG1基因5'调控区存在15个连锁的SNPs位点，其中g.28888A>T、g.28901G>T、g.30802C>T、g.30825A>G位点与初生体重显著相关（P<0.05），g.30737T>C位点与断奶体重显著相关（P<0.05），AAGGTAGGAAGGTTAATTAACCAATCCCAA单倍型与初生体重和断奶体重显著相关（P<0.05），说明PLAG1基因可作为湖羊生长性状选择的候选遗传标记。

关键词: 湖羊, PLAG1基因, 5'调控区, SNPs, 初生体重, 断奶体重

Abstract: This experiment aimed to clone the sequence of the 5'regulatory region of Hu sheep's PLAG1 gene, analyze the relationship between its polymorphisms and the early body weight, and look for the molecular markers for assisted selection of growth traits of Hu sheep. The 5'RACE method was used to identify the transcription initiation site of PLAG1 of Hu sheep. Four hundred fifty six weaned Hu sheep were selected as research objects, the nucleotide polymorphisms in the PLAG1 5' regulatory region were screened by cloning and sequencing, SPSS 18.0 was used to analyze the influence of different genotypes on the birth weight and weaning weight of Hu sheep. The results showed that the transcription initiation site was located at g.30535 locus of the PLAG1 gene. Fifteen SNPs were identified with 3 genotypes in PLAG1 5'regulatory region of Hu sheep population, and these SNPs were in Hardy-Weinberg equilibrium. Correlation analysis showed that the birth weight of Hu sheep with AA， GG, CC and AA genotypes at g.28888A>T, g.28901G>T, g.30802C>T and g.30825A>G loci, respectively were significantly higher than those of Hu sheep with corresponding heterozygous genotypes (P<0.05), and weaning weight of Hu sheep with CC genotype at g.30737T>C locus was significantly lower than that of Hu sheep with TC genotype (P<0.05). Linkage analysis showed that g.28888A>T, g.28901G>T, g.28913T>A, g.29717G>C, g.29724A>C, g.29725G>A, g.29768T>C, g.29771A>G, g.30141T>A, G.30144A>G, g.30691C>A, g.30694A>T, g.30737T>C, g.30802C>T, g.30825A>G were closely linked and detected 18 haplotypes. Further association analysis showed that 5 haplotypes with large sample size among them were associated with early body weight, the birth weight of Hu sheep with AAGGTAGGAAGGTTAATTAACCAATCCCAA haplotype was significantly higher than that of Hu sheep with ATGTTAGCACGATCAGTAAGCAATTCCTAG haplotype (P<0.05), the its weaning weight was significantly higher than that of AAGGAAGGAAGGTTAATTAACCAACCCCAA haplotype (P<0.05). The result indicated that there were 15 linkage SNPs loci in the 5' regulatory region of PLAG1 gene in Hu sheep population; g.28888A>T, g.28901G>T, g.30802C>T, g.30825A>G loci were significantly associated with birth weight (P<0.05); g.30737T>C locus was significantly associated with weaning weight (P<0.05); AAGGTAGGAAGGTTAATTAACCAATCCCAA haplotype was significantly associated with birth weight and weaning weight (P<0.05). All these data demonstrated that PLAG1 gene could be used as a candidate genetic marker for the growth traits selection of Hu sheep.

Key words: Hu sheep, PLAG1 gene, 5'regulatory region, SNPs, birth weight, weaning weight

中图分类号:

郭潇潇, 李隐侠, 王悦, 张含, 张俊, 钱勇, 孟春花, 王慧利, 仲跻峰, 曹少先. 湖羊PLAG1基因5'调控区多态性及其与早期体重的关联分析[J]. 畜牧兽医学报, 2021, 52(2): 331-343.

GUO Xiaoxiao, LI Yinxia, WANG Yue, ZHANG Han, ZHANG Jun, QIAN Yong, MENG Chunhua, WANG Huili, ZHONG Jifeng, CAO Shaoxian. Polymorphisms in the 5'Regulatory Region of PLAG1 Gene and Their Association with Early Body Weight of Hu Sheep[J]. Acta Veterinaria et Zootechnica Sinica, 2021, 52(2): 331-343.

参考文献 37

[1]	KAS K,VOZ M L,HENSEN K,et al.Transcriptional activation capacity of the novel PLAG family of zinc finger proteins[J].J Biol Chem,1998,273(36):23026-23032.
[2]	KAS K,VOZ M L,RÖIJER E,et al.Promoter swapping between the genes for a novel zinc finger protein and β-catenin in pleiomorphic adenomas with t(3;8)(p21;q12) translocations[J].Nat Genet,1997,15(2):170-174.
[3]	LANDRETTE S F,KUO Y H,HENSEN K,et al.Plag1 and Plagl2 are oncogenes that induce acute myeloid leukemia in cooperation with Cbfb-MYH11[J].Blood,2005,105(7):2900-2907.
[4]	CAO Y J,TAO Q S,KAO X M,et al.Hsa-circRNA-103809 promotes hepatocellular carcinoma develop-ment via microRNA-1270/PLAG1 like zinc finger 2 axis[J].Dig Dis Sci,2020,doi:10.1007/s10620-020-06416-x.
[5]	BANĚKOVČ M,URO-COSTE E,PTÁKOVÁ N,et al.What is hiding behind S100 protein and SOX10 positive oncocytomas? Oncocytic pleomorphic adenoma and myoepithelioma with novel gene fusions in a subset of cases[J].Hum Pathol,2020,103:52-62.
[6]	HENSEN K,VAN VALCKENBORGH I C C,KAS K,et al.The tumorigenic diversity of the three PLAG family members is associated with different DNA binding capacities[J].Cancer Res,2002,62(5):1510-1517.
[7]	TANG Q Q,WU W,XU X,et al.MiR-141 contributes to fetal growth restriction by regulating PLAG1 expression[J].PLoS One,2013,8(3):e58737.
[8]	QUEIMADO L,LOPES C,DU F,et al.Pleomorphic adenoma gene 1 is expressed in cultured benign and malignant salivary gland tumor cells[J].Lab Invest,1999,79(5):583-589.
[9]	HENSEN K,BRAEM C,DECLERCQ J,et al.Targeted disruption of the murine Plag1 proto-oncogene causes growth retardation and reduced fertility[J].Dev Growth Differ,2004,46(5):459-470.
[10]	MADISSOON E,DAMDIMOPOULOS A,KATAYA-MA S,et al.Pleomorphic adenoma gene 1 is needed for timely zygotic genome activation and early embryo development[J].Sci Rep,2019,9(1):8411.
[11]	WONG J,JUMA A R,TRAN S C,et al.Deficiency of the transcription factor PLAG1 results in aberrant coiling and morphology of the epididymis[J].Asian J Androl,2020,22(4):342-347.
[12]	ABI HABIB W,BRIOUDE F,EDOUARD T,et al.Genetic disruption of the oncogenic HMGA2-PLAG1-IGF2 pathway causes fetal growth restriction[J].Genet Med,2018,20(2):250-258.
[13]	KARIM L,TAKEDA H,LIN L,et al.Variants modulating the expression of a chromosome domain encompassing PLAG1 influence bovine stature[J].Nat Genet,2011,43(5):405-413.
[14]	FINK T,TIPLADY K,LOPDELL T,et al.Functional confirmation of PLAG1 as the candidate causative gene underlying major pleiotropic effects on body weight and milk characteristics[J].Sci Rep,2017,7(1):44793.
[15]	SONG Y X,XU L Y,CHEN Y,et al.Genome-wide association study reveals the PLAG1 gene for knuckle,biceps and shank weight in Simmental Beef Cattle[J].PLoS One,2016,11(12):e0168316.
[16]	NISHIMURA S,WATANABE T,MIZOSHITA K,et al.Genome-wide association study identified three major QTL for carcass weight including the PLAG1-CHCHD7 QTN for stature in Japanese Black cattle[J].BMC Genet,2012,13(1):40.
[17]	LEE S H,CHOI B H,LIM D,et al.Genome-wide association study identifies major loci for carcass weight on BTA14 in Hanwoo (Korean cattle)[J].PLoS One,2013,8(10):e74677.
[18]	OKADA Y,KAMATANI Y,TAKAHASHI A,et al.A genome-wide association study in 19633 Japanese subjects identified LHX3-QSOX2 and IGF1 as adult height loci[J].Hum Mol Genet,2010,19(11):2303-2312.
[19]	CHO Y S,GO M J,KIM Y J,et al.A large-scale genome-wide association study of Asian populations uncovers genetic factors influencing eight quantitative traits[J].Nat Genet,2009,41(5):527-534.
[20]	QIAO R M,GAO J,ZHANG Z Y,et al.Genome-wide association analyses reveal significant loci and strong candidate genes for growth and fatness traits in two pig populations[J].Genet Sel Evol,2015,47(1):17.
[21]	ZHONG J L,XU J W,WANG J,et al.A novel SNP of PLAG1 gene and its association with growth traits in Chinese cattle[J].Gene,2019,689:166-171.
[22]	LI Z,WU M L,ZHAO H D,et al.The PLAG1 mRNA expression analysis among genetic variants and relevance to growth traits in Chinese cattle[J].Anim Biotechnol,2020,31(6):504-511.
[23]	GUO Y M,HOU L J,ZHANG X F,et al.A meta analysis of genome-wide association studies for limb bone lengths in four pig populations[J].BMC Genet,2015,16(1):95.
[24]	ZHAO X D,REN W H,YANG W J,et al.Wnt pathway is involved in pleomorphic adenomas induced by overexpression of PLAG1 in transgenic mice[J].Int J Cancer,2006,118(3):643-648.
[25]	DECLERCQ J,VAN DYCK F,VAN DAMME B,et al.Upregulation of IGF and Wnt signalling associated genes in pleomorphic adenomas of the salivary glands in PLAG1 transgenic mice[J].Int J Oncol,2008,32(5):1041-1047.
[26]	HUANG W,LI B R,FENG H.PLAG1 silencing promotes cell chemosensitivity in ovarian cancer via the IGF2 signaling pathway[J].Int J Mol Med,2020,45(3):703-714.
[27]	ZHENG Y B,XU L,HASSAN M,et al.Bayesian modeling identifies PLAG1 as a key regulator of proliferation and survival in rhabdomyosarcoma cells[J].Mol Cancer Res,2020,18(3):364-374.
[28]	HENSEN K,BRAEM C,DECLERCQ J,et al.Targeted disruption of the murine Plag1 proto-oncogene causes growth retardation and reduced fertility[J].Dev Growth Differ,2004,46(5):459-470.
[29]	张天,李祥龙,周荣艳,等.不同毛色山羊皮肤组织Agouti基因剪接体类型研究[J].畜牧兽医学报,2015,46(11):1934-1943.ZHANG T,LI X L,ZHOU R Y,et al.Study on Agouti spliceosome types in goat skin with different coat color[J].Acta Veterinaria et Zootechnica Sinica,2015,46(11):1934-1943.(in Chinese)
[30]	JIA L F,MAO Y H,JI Q Q,et al.Decoding mRNA translatability and stability from the 5'UTR[J].Nat Struct Mol Biol,2020, 27(9):814-821.
[31]	CAPELL A,FELLERER K,HAASS C.Progranulin transcripts with short and long 5'untranslated regions (UTRs) are differentially expressed via post-transcriptional and translational repression[J].J Biol Chem,2014,289(37):25879-25889.
[32]	EISENHUT P,MEBRAHTU A,MORADI BARZADD M,et al.Systematic use of synthetic 5'-UTR RNA structures to tune protein translation improves yield and quality of complex proteins in mammalian cell factories[J].Nucl Acids Res,2020,doi:10.1093/nar/gkaa847.
[33]	BABU D,VANNELLI S,FANELLI A,et al.Variants in the 5'UTR reduce SHOX expression and contribute to SHOX haploinsufficiency[J].Eur J Hum Genet,2020,doi:10.1038/s41431-020-0676-y.
[34]	KARIM L,TAKEDA H,LIN L,et al.Variants modulating the expression of a chromosome domain encompassing PLAG1 influence bovine stature[J].Nat Genet,2011,43(5):405-413.
[35]	LITTLEJOHN M,GRALA T,SANDERS K,et al.Genetic variation in PLAG1 associates with early life body weight and peripubertal weight and growth in Bos taurus[J].Anim Genet,2012,43(5):591-594.
[36]	范家萌.NR6A1、PLAG1、LCORL基因在民猪群体内的单倍型及连锁不平衡分析[D].哈尔滨:东北农业大学,2015.FAN J M.The studies on the haplotypes of NR6A1、PLAG1 and LCORL and the association on Min pig population[D]. Harbin:Northeast Agricultural University,2015.(in Chinese)
[37]	ZHOU Z H,HUANG B Z,LAI Z Y,et al.The distribution characteristics of a 19-bp indel of the PLAG1 gene in Chinese cattle[J].Animals,2019,9(12):1082.