绵羊Y染色体特异性引物及SNPs的筛选

删除或更新信息，请邮件至freekaoyan#163.com(#换成@)

本站小编 Free考研考试/2021-12-26

曹学涛, 裴生伟, 张晋, 李发弟, 李刚, 李万宏, 乐祥鹏. 绵羊Y染色体特异性引物及SNPs的筛选[J]. 中国农业科学, 2018, 51(15): 2990-2999 https://doi.org/10.3864/j.issn.0578-1752.2018.15.014
CAO XueTao, PEI ShengWei, ZHANG Jin, LI FaDi, LI Gang, LI WanHong, YUE XiangPeng. Screening of Y Chromosome Specific Primers and Y-SNPs in Sheep[J]. Scientia Acricultura Sinica, 2018, 51(15): 2990-2999 https://doi.org/10.3864/j.issn.0578-1752.2018.15.014

0 引言

【研究意义】Y染色体基因在研究父系遗传多样性、维持雄性特性,精子发生和雄性繁殖力起着关键作用,Y染色体分子标记已经运用于绵羊父系起源进化研究及父系单倍型构建^[1,2,3]。然而Y染色体测序困难,造成绵羊Y染色体基因序列极度贫乏。因此本研究基于前期牛科动物Y染色体研究成果,进行绵羊Y染色体基因片段鉴定和Y-SNPs的筛选,以期为优秀种公羊的早期选择、胚胎性别的早期鉴定和绵羊品种的父系遗传多样性的保护提供理论依据。【前人研究进展】哺乳动物的性染色体,X和Y染色体,是由一对原始的常染色体进化而来^[4]。Y染色体是雄性动物特有的性染色体,由两部分组成：拟常染色体区域和Y染色体雄性特异区（male specific region of Y chromosome, MSY）组成。MSY区在减数分裂过程中不与X染色体发生重组交换,遵循严格的父系遗传。Y染色体结构非常复杂,包含有大量的重复序列和回文结构,使得其序列的组装极为困难^[5]。到目前为止,绵羊Y染色体只有SRY、ZFY等少数基因的部分序列已公布^[6,7,8]。MEADOWS等利用23对人、牛、绵羊Y染色体特异性引物扩得到绵羊AMELY、DBY、SRY、UTY、ZFY等Y染色体基因的部分片段^[9]。FENG等扩增得到岩羊ZFY基因639 bp的片段^[7]。刘帅兵克隆得到ZFY基因整个mRNA全长^[10]。熊勇等克隆得到藏系绵羊ZFY基因447 bp的片段^[11]。蒋利等利用山羊SRY基因序列设计引物,扩增得到阿勒泰羊、藏绵羊SRY基因整个编码区的雄性特异性片段,研究表明SRY基因在不同绵羊品种间具有高度的同源性^[12]。【本研究切入点】目前只有人类、猩猩、猕猴、小鼠和牛等少数物种的Y染色体完成了测序^[3,5,13-15]。关于绵羊Y染色体,少数研究者只扩增出基因的部分片段,未见克隆得到完整的Y染色体基因序列。【拟解决的关键问题】本研究基于牛科动物现有Y染色体特异引物鉴定绵羊Y染色体单拷贝基因,扫描Y染色体基因在不同绵羊品种中的SNPs,以期找到Y染色体分子标记,以期为筛选不同羊品种Y染色体SNPs提供基础,也为绵羊早期胚胎性别的鉴定提供了分子标记,为进一步研究绵羊Y染色体基因序列奠定基础。

1 材料与方法

试验于2016年10月至2017年5月在兰州大学草地农业生态系统国家重点实验室完成。

1.1 试验材料

496份公绵羊和14份公岩羊的血液样本,包括采自北京澳鑫牧业有限公司的萨福克羊（n=146）、白萨福克羊（n=91）、东弗里升羊（n=6）、特克塞尔羊（n=72）、南非肉用美利奴羊（n=17）、杜泊羊（n=32）,甘肃金昌中天羊业有限公司的湖羊（n=55）,甘肃甘南藏族自治州的藏羊（n=34）,宁夏盐池县的滩羊（n=43）,岩羊采自新疆动物园野生动物保护中心。采集3只湖羊母羊血液（甘肃金昌中天羊业有限公司）及3只岩羊母羊血液,作为阴性对照。同时在收集血样时测定萨福克羊、白萨福克羊、东弗里升羊、特克塞尔羊、南非肉用美利奴羊、杜泊羊周岁（周岁±15—20 d）的阴囊围直径。

1.2 DNA的提取和DNA池的构建

采用常规的酚-氯仿抽提法提取绵羊血液基因组DNA^[16]。用Nano drop 2000（Thermo,美国）检测提取的DNA浓度和纯度,并用1%的琼脂糖凝胶电泳检测DNA的质量。将质量检测合格绵羊和岩羊DNA统一稀释为50 ng·µL^-1,每个个体DNA取10 µL,以同一个品种尽量在一个DNA池的原则,每10个个体混成一个DNA池,共计51个DNA池,混合均匀后-20℃保存备用。

1.3 绵羊Y染色体特异性引物的筛选

查阅相关文献可知,目前在牛、牦牛、山羊等牛科动物Y染色体研究有29对引物（表1）,分别能扩增出MSY区的SRY、ZFY、USP9Y、UTY、OFD1Y、AMELY等基因。本研究的29对引物均由北京六合华大基因科技有限公司合成。实验以公羊DNA混合池为模板,以母羊DNA为阴性对照,灭菌水为空白对照,通过PCR扩增确定绵羊Y染色体特异性引物。PCR扩增体系为50 µL：模板2 µL,上、下游引物各1 µL,2×EasyTaq酶25 µL和ddH₂O 21µL。扩增反应条件：94℃预变性30 s;94℃变性30 s,退火30 s,72℃延伸1 min/1kb,34个循环;72℃延伸5 min;4℃保存。PCR产物用1％琼脂糖凝胶电泳检测。
Table 1
表1
表1引物序信息
Table 1Primer information

引物 Primers	序列 Sequences	退火温度 Anneal temperatures (℃)	参考文献 References
AMELY1	GCCTGCTTGACAGACTTA	—	[17]
AMELY1	AATGCTGTTAGGGATTTG	—	[17]
AMELY2	CCTGCCATCTACCTCTAT	—	[17]
AMELY2	GAACCGCTTTAACTATCC	—	[17]
AMELY4	TGGTTGGCTCTGGCGTAG	—	[17]
AMELY4	TCCTGCCCTCCAAACTCC	—	[17]
OFD1Y1	ACTCCTCACCCACAGACACC	—	[17]
OFD1Y1	CTGGATCGGGAAAGGGTAAT	—	[17]
OFD1Y9	CCTTCTCGCTTAGCATTTGG	—	[17]
OFD1Y9	CACATCTCATGGGAGGGACT	—	[17]
SRY1	ACAATATACAGTTGTTACTTCCA	—	[18]
SRY1	CAGAATTTGTGAATGTGC	—	[18]
SRY2	AAATAACTTCATAATGACACTT	—	[18]
SRY2	TGTGACCGGCTTAATTGGCTT	—	[18]
SRY3	ACTCTATGTCAACTTTCAAG	—	[18]
SRY3	GGCAATTTCAGATACAACTA	—	[18]
SRY4	GTCTGCTGCACCTTCATC	—	[18]
SRY4	CTTATTGTGGCCCAGGCTTGTC	—	[18]
引物 Primers	序列 Sequences	退火温度 Anneal temperatures (℃)	参考文献 References
SRY5	CCGGGCTATAAATATCGACCTC	—	[18]
SRY5	GATGAAACCTTGGGTCTCACAG	—	[18]
SRY6	GATGAAACCTTGGGTCTCACAG	55.4	[18]
SRY6	CTCACAAATTCATGGTACAGAG	55.4	[18]
SRY7	CAACTTTCAAGTTTGCCTTATGG	—	[19]
SRY7	ACAGCCCAATCCTGTTATATA	—	[19]
SRY8	GTATCTGAGATTGCCTGT	—	[19]
SRY8	CATCGTCGTTCAATACTC	—	[19]
SRY9	TTCATAACCGTAGACAAT	—	[19]
SRY9	CATGCATCGGGTTGCATA	—	[19]
SRY10	ACAGTCCAGCTGTGGTAC	—	[19]
SRY10	TTGTTACAGGGAAAGTCC	—	[19]
SRY11	GTCTGCTGCACCTTCATCCT	62.6	[20]
SRY11	GTTCATGGGTCGCTTGACGT	62.6	[20]
SRY12	CCACAGAAATCGCTTCCTGC	—	[18]
SRY12	TATGACTAAAGGAACTGAG	—	[18]
SRY13	GCTGACTGCCAGGACGTA	—	[19]
SRY13	TTGTTACAGGGAAAGTCC	—	[19]
SRY14	GCATGTAGAGACATTGCA	—	[19]
SRY14	TACAGTTTCACCATGGAC	—	[19]
USP9Y	AAACCCTTCAAGGTAATAAAACAAAA	55.4	[21]
USP9Y	CACAGCTCCTCAAAACCAGA	55.4	[21]
UTY	TGATGCCTATATTAGCCATTGACA	57.4	[22]
UTY	AAAATTCTTTATGATGTTCCATCC	57.4	[22]
ZFY1	CAGGTGAGGGCACATGAG	—	[18]
ZFY1	ATCACATTCGATGGCCTT	—	[18]
ZFY2	ATATGCTTGAAGAGACGACAAC	—	[18]
ZFY2	AGTCAGAAGACAAATGTCACA	—	[18]
ZFY3	TTCTAATTTGAAGACGCATGTG	50.8	[18]
ZFY3	CAACTTCTTTATGGTGTC GTG	50.8	[18]
ZFY4	GATCTCAGAGGAAATTCCAGTA	—	[18]
ZFY4	GATTGTGTAAACTTATCTTCTT	—	[18]
ZFY5	TCATTGGATCACTTATGCTC	—	[18]
ZFY5	TATGGATTCGCATGTGCTT	—	[18]
ZFY6	ATATGCTTGAAGAGACGACAAC	57.6	[17]
ZFY6	AGTCAGAAGACAAATGTCACA	57.6	[17]
ZFY7	TGAAGCAGGCACTTTGACAG	—	[21]
ZFY7	GGAGCATAAGTGATCCAATGAA	—	[21]
ZFY8	TCACATTGCAGCTTTAGGATTG	—	[23]
ZFY8	CCTTCACTTGGCAGATGGAT	—	[23]

The annealing temperatures were listed for those ovine Y-specific primers表中仅列出绵羊Y染色体特异性引物的退火温度

新窗口打开

1.4 绵羊Y染色体特异性基因片段SNPs筛选

将上述鉴定为绵羊Y染色体特异性引物,以DNA混合池为模板的PCR产物送北京六合华大基因科技有限公司进行测序。DNAstar和Chromas软件用于测序结果分析,DNAMAN用于多序列比对及同源性分析。经过混合池测序,在SRY11片段内鉴定得到一个G>A突变,利用Watcut软件分析（http://watcut. uwaterloo.ca/template.php?act=snp_new）发现该突变位点可被限制性内切酶Hpy188I识别。酶切反应体系：10×NEB Buffer 3 µL,PCR产物4 µL和5U Hpy188I内切酶,加灭菌水至25 µL。酶切反应条件为37℃温育3 h。取酶切消化液10 µL加样于2.5%琼脂糖凝胶上进行电泳。

1.5 数据处理与统计分析

利用SPSS 19.0软件分析基因位点与绵羊睾丸大小的相关性。先对数据进行描述性统计分析,确定是否存在离群值,根据数据特征,利用单因素方差分析（ANOVA）探究基因型效应。

2 结果

2.1 绵羊Y染色体特异性引物

通过PCR扩增,本研究发现29对引物中6对引物,包括ZFY3、SRY6、USP9Y、UTY、SRY11、ZFY6确定为Y染色体特异性引物（图1）,只在公羊DNA上出现扩增条带,在母羊DNA及水中没有扩增产物。此外17对引物未能出现扩增条带或非特异性扩增,6对引物母羊DNA能扩增出条带（电子版附图S1）。

显示原图|下载原图ZIP|生成PPT
图S1绵羊Y染色体非特异性PCR扩增结果 A：ZFY2 PCR扩增结果;B：SRY4 PCR扩增结果;C：SRY5 PCR扩增结果;D：ZFY9 PCR扩增结果;E：OFD1Y1 PCR扩增结果;F：OFD1Y9 PCR扩增结果;1-3公羊DNA,4-6母羊DNA,7-9ddH₂O,M.DNA相对分子质量标准DL2000
-->Fig. S1PCR products of Y nonspecific fragment in sheepA.PCR products of ZFY3 fragment; B.PCR products of SRY6 fragment; C.PCR products of USP9Y fragment; D.PCR products of UTY fragment; E.PCR products of SRY11 fragment; F.PCR products of ZFY3 fragment;1-3 rams DNA;4-6 ewes DNA;7-9 ddH₂O;M.DL2000 marker
-->

显示原图|下载原图ZIP|生成PPT
图1绵羊Y染色体特异性引物的PCR扩增 A：ZFY3基因PCR扩增结果;B：SRY6基因PCR扩增结果;C：USP9Y基因PCR扩增结果;D：UTY基因PCR扩增结果;E：SRY11基因PCR扩增结果;F：ZFY6基因PCR扩增结果;1-3公羊DNA,4-6母羊DNA,7-9 ddH2O,M.DNA相对分子质量标准DL2000
-->Fig. 1PCR products of ovine Y specific primersA.PCR products of ZFY3 genes fragment; B.PCR products of SRY6 genes fragment; C.PCR products of USP9Y genes fragment; D.PCR products of UTY genes fragment; E.PCR products of SRY11 genes fragment; F.PCR products of ZFY3 genes fragment;1-3 rams DNA;4-6 ewes DNA;7-9 ddH₂O;M.DL2000 marker
-->

2.2 绵羊Y染色体特异性DNA片段与其他物种的同源性分析

通过对上述6对特异性引物的PCR产物进行测序,分析发现ZFY3片段大小为614 bp,该片段与山羊ZFY第11外显子和部分CDS区相似度较高,相似度为98.84%。SRY6片段大小为461 bp,该片段与山羊SRY部分CDS区相似度较高,相似度为98.84%。USP9Y片段大小为444 bp,该片段与牛USP9Y第26外显子相似度较高,相似度为84.96%。UTY片段大小为256 bp,该片段与牛第19外显子相似度较高,相似度为81.51%。SRY 11片段大小为440 bp,该片段与山羊SRY基因序列相似度较高,相似度为90.41%,该基因片段的部分片段（237 bp）与NCBI公布的绵羊的SRY基因CDS区（AY604735.1）相似度达100%。ZFY3片段大小为621 bp,该片段与牛ZFY基因第11外显子和部分CDS区相似度较高,相似度为97.67%。
用DNAMAN生物软件,将绵羊ZFY3、SRY6、USP9Y、UTY、SRY11、ZFY6基因序列与山羊、牛、牦牛和岩羊相应基因序列进行序列比对（电子版附图S2—S7）,计算不同物种基因序列间的同源性大小（表2）。结果表明,绵羊ZFY3基因片段与岩羊、山羊、牛、牦牛各物种间同源性大小分别为98.84%、97.89%、97.08%和95.79%;绵羊SRY6基因片段与岩羊、山羊、牛、牦牛各物种间同源性大小分别为95.25%、93.93%、90.54%和90.32%;绵羊USP9Y基因片段与牛、牦牛各物种间同源性大小分别为84.96%、84.67%;绵羊UTY基因片段与岩羊、黄牛、牦牛各物种间同源性同为96.77%、81.51%和81.51%;绵羊SRY11基因片段与岩羊、山羊、牛、牦牛各物种间同源性大小分别为92.63%、90.41%、84.97%和84.75%;绵羊ZFY6基因片段与岩羊、山羊、牛、牦牛各物种间同源性大小分别为99.06%、89.11%、97.67%和86.76%。由于USP9Y基因是从岩羊母羊DNA中扩增出,表明该片段并不是岩羊Y染色体片段。

显示原图|下载原图ZIP|生成PPT
图S2不同物种ZFY3基因序列比对结果
-->Fig. S2Comparison of ovine ZFY3 gene sequence with other species
-->

显示原图|下载原图ZIP|生成PPT
图S3不同物种SRY6基因序列比对结果
-->Fig. S3Comparison of ovine SRY6 gene sequence with other species
-->

显示原图|下载原图ZIP|生成PPT
图S4不同物种USP9Y基因序列比对结果
-->Fig. S4Comparison of ovine USP9Y gene sequence with other specie
-->

显示原图|下载原图ZIP|生成PPT
图S5不同个物种UTY基因序列比对结果
-->Fig. S5Comparison of ovine UTY gene sequence with other specie
-->

显示原图|下载原图ZIP|生成PPT
图S6不同物种SRY11基因序列比对结果
-->Fig. S6Comparison of ovine SRY11 gene sequence with other specie
-->

显示原图|下载原图ZIP|生成PPT
图S7不同物种ZFY6基因序列比对结果
-->Fig. S7Comparison of ovine ZFY6 gene sequence with other specie
-->

Table 2
表2
表2绵羊Y染色体基因与不同物种同源基因的比较
Table 2Comparison of homologous genes between sheep Y chromosome and others species (%)

基因片段 Gene fragments	物种 Species	绵羊 Sheep	岩羊 Bharal	山羊 Goat	牛 Cattle
ZFY3	岩羊Bharal	98.84	—	—	—
	山羊Goat	97.89	96.1	—	—
	牛Cattle	97.08	94.63	98.21	—
	牦牛Yak	95.79	95.88	96.59	97.72
SRY6	岩羊Bharal	95.25	—	—	—
	山羊Goat	93.93	95.74	—	—
	牛Cattle	90.54	89.45	91.97	—
	牦牛Yak	90.32	89.24	92.19	99.35
USP9Y	黄牛Cattle	84.96	—	—	—
USP9Y	牦牛Yak	84.67	—	—	83.41
UTY	岩羊Bharal	96.77		—	—
	牛Cattle	81.51	80.46	—	—
	牦牛Yak	81.51	80.77	—	86.77
SRY11	岩羊Bharal	92.63	—	—	—
	山羊Goat	90.41	97.19	—	—
	牛Cattle	84.97	89.39	90.33	—
	牦牛Yak	84.75	89.18	90.11	99.32
ZFY6	岩羊Bharal	99.06	—	—	—
	山羊Goat	89.11	88.24	—	—
	牛Cattle	97.67	96	88.34	—
	牦牛Yak	86.76	86.2	84.09	87.87

新窗口打开

2.3 绵羊Y染色体SNPs的筛选及基因效应分析

通过测序,ZFY3、SRY6、USP9Y、UTY、ZFY6基因序列均未发现突变位点。在萨福克和白萨福克羊群体中,SRY11片段混合池测序出现双峰,进一步对其中一个混池的每个个体分别进行测序,验证为G>A的突变（图2）。对萨福克及白萨福克的个体进行PCR-RFLP检测,PCR产物经Hpy188I酶切后产生AA和GG 2种基因型（图3）。AA基因型个体的扩增产物能被Hpy188I酶切后将产生4个片段,分别为54、63、151和172 bp;GG基因型个体的扩增产物能被Hpy188I酶切后将产生3个片段,分别为63、172和205bp。SRY11片段G>A突变位点在萨福克羊、白萨福克羊、东弗里升羊、特克塞尔羊、南非肉用美利奴羊、杜泊羊、湖羊、藏羊、滩羊公羊群体中的基因型频率、等位基因频率见表3。结合混合池测序及酶切结果发现,在萨福克羊、白萨福克羊中存在AA和GG基因型,其中在白萨福克羊群体中AA基因型频率为0.253,GG基因型频率为0.747;在萨福克羊群体中AA基因型频率为0.014,GG基因型频率为0.986。结果显示GG基因型在萨福克羊和白萨福克羊群体中为优势基因型;而东弗里升羊、特克塞尔羊、南非肉用美利奴羊、杜泊羊、湖羊、藏羊、滩羊公羊群体中只有GG基因型。通过对岩羊个体的混合池测序,在SRY11片段也未发现突变位点。性状关联分析表明,在白萨福克羊群体中,GG基因型个体的睾丸大小显著的高于AA基因型（表4）。

显示原图|下载原图ZIP|生成PPT
图2SRY11片段测序 A. AA基因型; B. GG基因型
-->Fig. 2Sequencing results of SRY11 fragment A. AA genotype; B. GG genotype
-->

显示原图|下载原图ZIP|生成PPT
图3绵羊SRY11片段酶切1,3：AA基因型;2,4：GG基因型;M为DNA相对分子质量标准DL2000
-->Fig. 3PCR-RFLP pattern for SRY11 fragment with Hpy188I digestion1, 3: AA genotype; 2, 4: GG genotype; M: DL2000 marker
-->

Table 3
表3
表3SRY11片段G>A 位点在不同绵羊品种中的基因型和等位基因频率
Table 3Genotypic and allelic frequencies of G>A in SRY11 fragment among different sheep breeds

品种 Species	总数（只） Total	AA基因型（只） Genotype AA	GG基因型（只） Genotype BB	AA基因型频率 AA genotype frequency	GG基因型频率 GG genotype frequency
白萨福克White Suffolk sheep	91	23	68	0.253	0.747
萨福克Suffolk sheep	146	2	144	0.014	0.986
特克赛尔Texel sheep	72	0	72	0	1
杜泊Dorper sheep	32	0	32	0	1
东弗里升East Friesian sheep	6	0	6	0	1
南非肉用美利奴South African Mutton Merino	17	0	17	0	1
湖羊Hu sheep	55	0	55	0	1
藏羊Tibetan sheep	34	0	34	0	1
滩羊Tan sheep	43	0	43	0	1
岩羊Bharal	14	0	14	1	1
总计Total	510	25	485	0.049	0.951

新窗口打开
Table 4
表4
表4 SRY11片段G>A位点与白萨福克羊睾丸大小之间的关联分析
Table 4Correlation analysis of the G>A in SRY11 and testis size in White Suffolk sheep

繁殖性状 Reproductive traits	基因型（平均值±标准差） Genotypes (Mean ± SD)		P值 P value
繁殖性状 Reproductive traits	AA	GG	P值 P value
睾丸大小 Testis size	25.79±4.49	27.06±3.58	0.029

新窗口打开

3 讨论

在目前研究绵羊起源进化及单倍型分析中,更多的是以线粒体基因遗传多样性为基础的母系遗传研究,而利用以Y-SNPs为基础的父系遗传的研究较少。绵羊的父系及母系遗传相结合才能更加全面、系统的分析绵羊起源进化、单倍型及分类。本研究筛选得到了6对绵羊Y染色体特异性引物,均为首次扩增得到,与MEADOWS等筛选得到9对绵羊Y染色体特异性引物无重复^[9]。SRY（sex-determining region on Y chromosome, SRY）位于哺乳动物Y染色体,是性别决定直接相关的基因。该基因从人的Y染色体上首次分离得到,并被确定为睾丸决定因子^[24]。MEADOWS等通过测序在SRY基因5′启动区发现8个SNPs,结合Y染色微卫星SRYM18鉴定绵羊及野生绵羊单倍型^[25]。MEADOWS等发现得到绵羊SRY第一外显子与牛的同源性为94.7%,并在5′启动区鉴定得到了一个A>G的突变^[9]。在本次研究扩增得到了两个SRY基因片段,其与牛的同源性分别为90.54%和84.97%,这与前人的研究相似,说明绵羊和牛SRY基因具有较高的同源性。本研究在SRY11片段鉴定到1个G>A的同义突变,通过Hpy188I酶切分析,发现只有在萨福克羊、白萨福克羊中有AA和GG 2种基因型,而在其他7个绵羊品种和岩羊中只有GG基因型,暗示该位点可能只存在于萨福克羊和白萨福克羊的血统中,但还需要在更多的绵羊品种群体中进行验证。此外该突变位点与白萨福克羊周岁的睾丸大小紧密相关,即GG基因型个体的睾丸大小显著大于AA型个体。本研究为首次发现Y染色体突变与绵羊繁殖性状相关。该位点在更大的白萨福克羊群体进行验证后可用于白萨福克种公羊的选育。
ZFY基因即Y连锁的锌指蛋白因子（zinc finger protein, Y-linked）,位于Y染色体短臂（Yp11.3）上,由11个外显子和1个随机重复区域的13个“锌-指”结构组成,在雄性动物睾丸的生长发育过程中发挥作用^[26,27]。毛德才研究发现牦牛、猩猩、人、牛等9个物种ZFX/ZFY间的核苷酸序列有较高同源性^[28]。熊勇等研究表明藏系绵羊ZFY基因与猪、摩弗伦羊、狼、人类等其他物种之间具有很高的同源性^[11]。LAWSON等通过进化树分析发现：绵羊ZFY与岩羊、山羊、牛ZFY的亲缘关系较近^[8]。ARKADI克隆得到了整个牛ZFY基因整个编码区及马和猪ZFY基因编码区部分序列,通过比对发现ZFY基因编码区在3个物种的同源性较高^[29]。本研究发现绵羊ZFY3和ZFY6片段与岩羊、山羊、牛、牦牛等物种间同源性大小在89.11%—99.6%,表明该基因在不同物种间也具有较高的同源性,与前人的研究基本一致。AASEN等^[30]利用ZFY片段的扩增及酶切鉴别人、牛、绵羊和山羊的胚胎性别,该研究的结果也可用于绵羊胚胎性别的早期鉴定。
SATOU等^[31]研究发现USP9Y基因的突变影响小鼠睾丸重量。有****发现在牛USP9Y基因第26内含子发现了一个81 bp的插入缺失,并用于构建牛Y染色体单倍型及区分瘤牛和普通牛^[21,32-33]。本研究在USP9Y基因片段中没有检测到遗传变异,今后需要在更多的绵羊品种进行Y-SNPs扫描。先前的研究利用BAC的方法在绒山羊上鉴定到USP9Y基因348 bp的片段,其与家猫、人、恒河猴、黑猩猩等的同源性在89％以上^[34]。本研究扩增得到的绵羊USP9Y基因与牛、牦牛各物种间同源性大小分别为84.96%、84.67%,表明该基因在进化过程中具有一定的保守性。USP9Y基因在岩羊母羊DNA中扩增出,表明该片段并不是岩羊Y染色体片段,因此本研究的结果可作为区分岩羊和绵羊的分子标记。
UTY基因即Y连锁的泛转录三角形四肽重复因子,在人上的研究发现UTY基因有6个不同的mRNA转录本,UTY蛋白被认为是未成熟的组织适合性抗原,可启动男性特异性细胞毒性T淋巴细胞反应^[35,36]。MEADOWS等扩增得到了UTY基因外显子1,通过比对发现与牛的同源性为98.1%^[9]。本研究扩增得到的绵羊UTY基因与岩羊、牛、牦牛各物种间同源性达81.51%—96.77%,表明该基因在进化过程中相对保守。GINJA等扩增牛UTY基因19外显子,发现有一个C>A的突变^[23]。然而在本研究中在UTY基因上并未检测到突变位点。
本研究首次在SRY发现一个G>A的突变,在ZFY3、SRY6、USP9Y、UTY、ZFY6片段均未发现突变,其原因可能与样本量大小及品种差异有关。但是本研究的结果为今后我国绵羊品种父系遗传多样性的研究奠定了基础。当然,为了系统和全面研究绵羊Y染色体,必须依赖绵羊Y染色体完整的序列,完成绵羊Y染色体测序是研究绵羊父系遗传多样性和绵羊雄性繁殖力育种亟待解决的难题。

4 结论

本研究首次扩增得到绵羊ZFY、SRY、USP9Y和UTY 4个基因的6个片段。通过比对发现绵羊ZFY3、SRY6、USP9Y、UTY、ZFY6片段与岩羊、牛、山羊、牦牛的具有较高的同源性。首次在SRY11片段发现G>A的突变,该突变位点只存在于萨福克和白萨福克公羊群体中,并与白萨福克羊睾丸大小紧密相关。
（责任编辑林鉴非）
The authors have declared that no competing interests exist.