李昆諭 洪瓊花 陶林 江炎庭 歐陽依娜 儲明星 劉玉芳

摘要 [目的]分析腺苷酸環(huán)化酶1（ADCY1）和胰凝乳蛋白酶（CYM）基因單核苷酸多態(tài)性（SNPs）與山羊產羔性狀的關聯。[方法]通過 MassARRAYSNP分型技術對ADCY1和CYM基因SNPs位點進行分型，檢測其在云上黑山羊（n=544）、濟寧青山羊（n=133）和遼寧絨山羊（n=91）3個群體中的遺傳學特征，并將其多態(tài)性與產羔性能（產羔數、初生窩重、斷奶窩重）進行關聯分析。[結果]群體遺傳學結果表明，ADCY1基因4個突變在云上黑山羊、濟寧青山羊、遼寧絨山羊中均為低度多態(tài)（PIC≤0.25）;CYM基因g.87624254G>A和g.87627082G>A 2個位點在3種山羊中均為中度多態(tài)（0.25<PIC<0.50）;卡方檢驗顯示，除ADCY1 g.43996583C>T在遼寧絨山羊中處于Hardy-Weinberg不平衡狀態(tài)（P<0.05），ADCY1基因其余位點和CYM基因g.87624269 G>A、g.87627082G>A 2個位點在3個群體中均處于Hardy-Weinberg平衡狀態(tài)（P>0.05）。關聯分析結果表明，在云上黑山羊中ADCY1基因g.43979904G>A、g.43974472G>A，CYM基因g.87624254G>A位點多態(tài)性與產羔數存在顯著關聯（P<0.05）;ADCY1基因g.43996550G>A、g.43996583C>T，CYM基因g.87624254G>A、g.87627082G>A位點多態(tài)性與初生窩重存在顯著關聯（P<0.05）。[結論]ADCY1基因g.43974472G>A和g.43979904G>A位點可作為云上黑山羊產羔數選擇的潛在分子標記。ADCY1基因g.43996550G>A、g.43996583C>T位點，CYM基因g.87627082G>A位點可作為云上黑山羊初生窩重選擇的潛在分子標記;CYM基因g.87624254G>A位點可作為云上黑山羊窩重選擇的潛在分子標記。

關鍵詞 山羊;產羔數;窩重;ADCY1基因;CYM基因

中圖分類號 S827? 文獻標識碼 A? 文章編號 0517-6611（2021）24-0141-05

doi：10.3969/j.issn.0517-6611.2021.24.033

Polymorphism of Goat ADCY1 and CYM and Its Association with Litter Size

LI Kun-yu1， 2， HONG Qiong-hua3， TAO Lin2 et al

（1. College of Life Science and Food Engineering， Hebei University of Engineering， Handan， Hebei 056001;2. Key Laboratory of Animal Genetics， Breeding and Reproduction of Ministry of Agriculture and Rural Affairs， Institute of Animal Science， Chinese Academy of Agricultural Sciences， Beijing 100193;3.Yunnan Animal Science and Veterinary Institute， Kunming， Yunnan 650224）

Abstract [Objective]This study aimed to analyze the association between single nucleotide polymorphisms （SNPs） in the adenylate cyclase 1 （ADCY1） and chymotrypsin （CYM） genes and kidding performance in goats.[Method]MassARRAYSNP typing technology was applied to genotype the SNPs of ADCY1 and CYM genes and detect their genetic characteristics in Yunshang black goat （n=544）， Jining grey goat （n=133） and Liaoning cashmere goat （n=91）. Then the association was analyzed between ADCY1 and CYM and kidding performance （litter size， litter weight at birth， litter weight at weaning） in different goats. [Result]The results of population genetics showed that the 4 mutations of ADCY1 gene were low polymorphisms in Yunshang black goat， Jining grey goat， and Liaoning cashmere goat （PIC≤0.25）;CYM gene g.87624254G>A and g.87627082G>A both candidate locus were moderately polymorphic in the three goats （0.25<PIC<0.5）. Chi-square test showed that， except ADCY1 g.43996583C>T in Liaoning cashmere goats， which is in Hardy-Weinberg imbalance （P<0.05）， the other locus are in Hardy-Weinberg equilibrium （P>0.05）;CYM gene g.87624269 two locus G>A and g.87627082G>A were in Hardy-Weinberg equilibrium in the three populations （P>0.05）. The results of association analysis showed that in Yunshang black goats， the ADCY1 gene g.43974472G>A， g.43979904G>A， and CYM gene g.87624254G>A polymorphisms were significantly associated with litter size （P<0.05）;ADCY1 gene g.43996550G>A， g.43996583C>T， CYM gene g.87624254G>A， g.87627082G>A polymorphisms were significantly associated with the primary litter weight （P<0.05）. [Conclusion]In summary， ADCY1 gene g.43974472G>A and g.43979904G>Acan be used as a potential molecular marker for litter size in Yunshang black goats.ADCY1 geneg.43996550G>A and g.43996583C>T locus， CYM gene g.87627082G>A locus can be used as potential molecular markers for litter weight at birth in Yunshang black goats; CYM gene g.87624254G>A locus can be used as potential molecular markers for litter weight selection in Yunshang black goats.

Key words Goat;Litter size;Litter weight;ADCY1 gene;CYM gene

基金項目 國家肉羊產業(yè)技術體系專項（CARS-38）;河北省自然科學基金青年項目（C2019402261）;云南省科技創(chuàng)新人才計劃（2018HC017）;云南省萬人計劃——產業(yè)技術領軍人才項目;中國農業(yè)科學院科技創(chuàng)新工程（CAAS-ZDRW202106和ASTIP-IAS13）。

作者簡介 李昆諭（1996—），男，甘肅定西人，碩士研究生，研究方向：動物遺傳育種;洪瓊花（1968—），女，云南鶴慶人，研究員，從事羊遺傳育種研究。李昆諭和洪瓊花為共同第一作者。通信作者：儲明星，研究員，從事羊優(yōu)異繁殖性狀分子機理研究;劉玉芳，博士，從事動物遺傳育種研究。

收稿日期 2021-04-12

山羊（Capra hircus）又稱夏羊，是最古老的馴養(yǎng)動物之一[1]。馴化的山羊可為人類提供肉、奶、皮、毛等主要生活資料，在人類農業(yè)文明和經濟發(fā)展中起著重要作用。產羔數是一個復雜的性狀，受到多種因素的影響，且產羔數是一個遺傳力較低（0.1左右）的性狀，是一個受多個基因控制的復雜的閾性狀[2-3]。利用傳統的育種手段在短時間內難以快速對其進行改良[4-5]。因此，利用分子標記輔助選擇技術，可為快速、高效地改良山羊產羔性狀提供有效位點，顯著提高產羔性狀選擇效率，加快育種進程，為進一步解析山羊高繁殖力形成的分子機理提供理論基礎。

腺苷酸環(huán)化酶1（adenylate cyclase 1，ADCY1）是ADCY蛋白的家族成員，定位于山羊4號染色體上，擁有20個外顯子。在有機體內，ADCY1負責將ATP催化成環(huán)AMP（cAMP），cAMP作為輔助因子參與細胞活動的多種進程[6-7]。Buell等[8-9]研究表明，cAMP可影響卵母細胞減數分裂的過程，其濃度可以影響卵母細胞成熟，進而影響胚胎的質量。Islam等[10]研究表明，在6個中國山羊品種中ADCY1被特異性選擇，并與山羊的繁殖力顯著相關。胰凝乳蛋白酶（chymotrypsin，CYM）具有9個外顯子，定位于山羊3號染色體上，其功能與胰蛋白酶相似，都屬參與蛋白質氨基酸代謝酶，也是一種水解酶。在檢測與胰腺分泌功能不全相關聯的胰腺疾病的篩選試驗方面有一定的應用，然而其對動物生殖的作用尚鮮見報道[11]。Lai等[12]通過對2個品種奶山羊全基因組深度測序篩選到了與產羔性狀顯著相關的重要候選基因CYM 。該研究利用前期對云上黑山羊（云上黑山羊是以云嶺山羊為母本、努比山羊為父本培育而成的我國第一個肉用黑山羊品種[13-14]）的重測序結果，并結合文獻閱讀，推測ADCY1和CYM基因可能對山羊產羔數有一定的影響，在前期數據結果中篩選出與產羔數相關的8個SNPs位點，其中，在ADCY1基因的SNPs分別為g.43974472G>A、g.43979904G>A、g.43996550G>A和g.43996583C>T，CYM基因的SNPs分別為g.87621057G>A、g.87624254G>A、g.87624269G>A和g.87627082G>A，采用 MassARRAYSNP分型技術，在云上黑山羊、濟寧青山羊和遼寧絨山羊群體中對候選基因的基因型進行檢測，分析其多態(tài)性與產羔數性狀之間的關聯性，以期為山羊的高繁殖力分子標記輔助選擇育種提供參考。

1 材料與方法

1.1 供試樣品

共采集768只山羊（2～5歲）的血液，其中云上黑山羊544只，濟寧青山羊133只，遼寧絨山羊91只（表1）。將濟寧青山羊視為高繁殖力群體，遼寧絨山羊視為低繁殖力群體[13-14]。所有試驗用山羊飼養(yǎng)條件環(huán)境一致，云上黑山羊擁有至少1胎產羔數記錄，部分擁有初生窩重和斷奶窩重（3月齡）記錄，所有羊均采用頸靜脈采血（10 mL/只），使用EDTA-K2抗凝，-20 ℃保存。

1.2 血液DNA的提取

試驗樣本使用酚氯仿法提取基因組DNA，采用Nano Drop2000檢測DNA樣本濃度，1.5%瓊脂糖凝膠電泳檢測DNA質量。

1.3 基因分型

采用Sequenom MassARRAYSNP分型技術對ADCY1基因和CYM基因突變位點進行檢測，相關引物信息見表2。分型樣品為DNA，每個樣品需要量為20 μL，DNA濃度為40～80 ng/μL。

1.4 統計分析

應用Microsoft Excel 2020軟件統計山羊ADCY1和CYM基因突變位點的基因型頻率、等位基因頻率、多態(tài)信息含量（PIC）、雜合度（He）和有效等位基因數（Ne），并進行Hardy-Weinberg平衡檢驗。采用 IBM SPSS Statistics 25軟件中的一般線性模型，多重比較采用LSD法，對山羊基因型與產羔表型數據進行關聯分析，所有數據以平均值±標準誤表示（**P<0.01，*P<0.05）。

2 結果與分析

2.1 ADCY1基因多態(tài)性分析

利用Sequenom MassARRAYSNP對768只山羊血液DNA進行基因分型。結果顯示，ADCY1基因4個SNPs位點在云上黑山羊、濟寧青山羊、遼寧絨山羊中存在多態(tài)性，其中g.43974472G>A、g.43979904G>A和g.43996550G>A位點基因型為GG、GA、AA，g.43996583C>T位點基因型為CC、CT、TT（圖1）。

統計濟寧青山羊（高繁）和遼寧絨山羊（低繁）品種ADCY1基因4個SNPs位點的基因型頻率和等位基因頻率（表3），ADCY1基因g.43979904 G>A和43996550 G>A位點的基因型頻率和等位基因頻率在高繁和低繁山羊群體間差異達到顯著水平（P<0.05），且高繁群體和低繁群體的優(yōu)勢基因型均為GG，優(yōu)勢等位基因為G;g.43996583 C>T位點的基因型頻率在高繁和低繁山羊群體間差異均達到極顯著水平（P<0.01）。

群體遺傳學統計表明，ADCY1基因4個突變在云上黑山羊、濟寧青山羊、遼寧絨山羊中均為低度多態(tài)（PIC≤0.25）?？ǚ竭m應性檢驗結果表明，除ADCY1g.43996583C>T在遼寧絨山羊中處于Hardy-Weinberg不平衡狀態(tài)（P<0.05）外，ADCY1其余位點均處于Hardy-Weinberg平衡狀態(tài)（P>0.05）（表4）。

2.2 CYM基因多態(tài)性分析

基因分型結果表明（圖2），CYM基因4個候選位點在3個群體中均存在多態(tài)性。其中，CYM基因g.87624254G>A和g.87627082G>A 2個位點在云上黑山羊、濟寧青山羊、遼寧絨山羊中存在GG、GA 和AA 3種基因型（表4）。

由表3可知，CYM基因g.87621057G>A和g.87624269G>A位點的基因型頻率和等位基因頻率在高繁和低繁山羊群體間差異均達到顯著水平（P<0.05），但g.87627082G>A位點的基因型頻率和等位基因頻率在高繁和低繁山羊群體間差異均不顯著。其中g.87624254G>A 和g.87627082G>A位點在高繁和低繁山羊群體中的優(yōu)勢基因型為GA，g.87621057G>A和g.87624269G>A位點在高繁和低繁山羊群體中的優(yōu)勢基因型為GG，優(yōu)勢等位基因均為G。

由表4可知，CYM基因g.87624254G>A和g.87627082G>A 2個候選位點在3種山羊中均為中度多態(tài)（0.25<PIC<0.50），CYM基因g.87624269G>A在3種山羊中均為低度多態(tài)（PIC≤0.25）?？ǚ綑z驗表明，CYM基因g.87624269 G>A和g.87627082G>A 2個位點在3個群體中均處于Hardy-Weinberg平衡狀態(tài)（P>0.05）。

2.3 候選SNPs位點與云上黑山羊繁殖性能的關聯分析

ADCY1基因g.43974472G>A位點AA型的產羔數顯著高于GG型（P<0.05），ADCY1基因g.43979904G>A位點GG和GA型的產羔數顯著高于AA型（P<0.05）（表5）。CYM基因g.87624254G>A位點GG型的產羔數顯著高于GA型（P<0.05）。ADCY1基因g.43996550G>A、g.43996583C>T，CYM基因g.87627082G>A 位點多態(tài)性與云上黑山羊初生窩重存在顯著關聯（P<0.05）。CYM基因g.87624254G>A位點多態(tài)性與云上黑山羊初生窩重和斷奶窩重存在顯著關聯（P<0.05）。

3 討論

3.1 ADCY1基因與生殖性狀的關系

cAMP作為第二信使的一種，是細胞內信號轉導的啟動組成部件之一，與細胞活動及生殖密切相關[8，15-17]。Chang等[18-19]研究發(fā)現，抑制cAMP表達可使芳香化酶的表達降低，從而影響顆粒細胞的生長發(fā)育，ADCY1作為cAMP形成所必需的物質，對于細胞活動尤為重要[20]。因此，ADCY1基因與動物的生殖過程息息相關，或可作為研究動物生殖過程的一個標記。

3.2 ADCY1和CYM基因多態(tài)性及其與云上黑山羊產羔數、窩重之間的關系

目前，關于ADCY1和CYM基因與山羊產羔數之間關系的報道較少。該試驗發(fā)現，ADCY1基因g.43974472G>A位點AA型個體產羔數高于GG型個體0.32只，GA型個體產羔數高于GG型0.07只，筆者推測該位點G>A的突變或是一個優(yōu)良的突變。ADCY1基因g.43979904G>A位點G>A的突變導致其產羔數下降，推測該位點G>A的突變或是一個不利的突變。ADCY1基因g.43996550G>A和g.43996583C>T位點的突變性都與云上黑山羊的初生窩重顯著相關，或可作為篩選云上黑山羊初生窩重的潛在分子標記。

CYM基因g.87624254G>A位點GG型個體相比GA型產羔數高出0.17只，而GA型個體初生窩重顯著高于GG型，因此該突變位點還可能需要進行進一步的驗證。CYM基因g.87627082G>A位點初生窩重雜合型個體顯著高于純合型個體，或可作為篩選云上黑山羊初生窩重的潛在分子標記。綜上所述，ADCY1基因g.43974472G>A和g.43979904G>A位點或可作為云上黑山羊產羔數選擇的潛在分子標記。ADCY1基因g.43996550G>A、g.43996583C>T位點，CYM基因g.87627082G>A位點或可作為云上黑山羊初生窩重的潛在分子標記，g.87624254G>A位點可作為云上黑山羊窩重選擇的潛在分子標記。

4 結論

ADCY1基因g.43974472G>A和g.43979904G>A位點與云上黑山羊產羔數存在顯著關聯，可作為云上黑山羊產羔數選擇的潛在分子標記。ADCY1基因g.43996550G>A、g.43996583C>T位點，CYM基因g.87627082G>A位點與云上黑山羊初生窩重存在顯著關聯，可作為云上黑山羊初生窩重選擇的潛在分子標記。CYM基因g.87624254G>A位點與云上黑山羊初生窩重和斷奶窩重存在顯著關聯，可作為云上黑山羊窩重選擇的潛在分子標記。

參考文獻

[1] MACHUGH D E，BRADLEY D G.Livestock genetic origins：Goats buck the trend[J].PNAS，2001，98（10）：5382-5384.

[2]? 狄冉，潘章源，儲明星，等.魯中肉羊GTF2A1基因多態(tài)性及其與產羔數的關聯分析[J].中國草食動物科學，2019，39（4）：15-17，80.

[3]? 湯繼順.利用轉錄組測序和蛋白質組學分析篩選綿羊多羔候選基因的研究[D].北京：中國農業(yè)科學院，2019.

[4]? 何煒.分子標記在山羊遺傳育種中的應用[J].青海畜牧獸醫(yī)雜志，2005，35（5）：35-37.

[5]? 蘭蓉，朱蘭，姚新榮，等.山羊產羔數全基因組關聯分析[J].畜牧獸醫(yī)學報，2015，46（4）：549-554.

[6]? ZOU T，LIU J Y，SHE L，et al.A perspective profile of ADCY1 in cAMP signaling with drug-resistance in lung cancer[J].J Cancer，2019，10（27）：6848-6857.

[7]? SIMONDS W F.G protein regulation of adenylate cyclase[J].Trends Pharmacol Sci，1999，20（2）：66-73.

[8] BUELL M，CHITWOOD J L，ROSS P J.cAMP modulation during sheep in vitro oocyte maturation delays progression of meiosis without affecting oocyte parthenogenetic developmental competence[J].Anim Reprod Sci，2015，154：16-24.

[9]? ROSE R D，GILCHRIST R B，KELLY J M，et al.Regulation of sheep oocyte maturation using cAMP modulators[J].Theriogenology，2013，79（1）：142-148.

[10]? ISLAM R，LI Y F，LIU X X，et al.Genome-wide runs of homozygosity，effective population size，and detection of positive selection signatures in six Chinese goat breeds[J].Genes（Basel），2019，10（11）：1-25.

[11]? 史艷霞，張永軍，竇寶峰，等.Bt抗性和敏感棉鈴蟲幼蟲中腸類胰蛋白酶和類胰凝乳蛋白酶基因克隆及半定量測定[J].應用與環(huán)境生物學報，2009，15（5）：630-633.

[12]? LAI F N，ZHAI H L，CHENG M，et al.Whole-genome scanning for the litter size trait associated genes and SNPs under selection in dairy goat（Capra hircus）[J].Sci Rep，2016，6：1-12.

[13]? 陶林，楊紅遠，江炎庭，等.云上黑山羊體尺性狀對體重影響的通徑分析[J].家畜生態(tài)學報，2020，41（7）：18-22.

[14]? 陶林，楊紅遠，江炎庭，等.云上黑山羊體重與成年性狀的典型相關分析[J].中國農學通報，2020，36（11）：108-112.

[15]? MEHLMANN L M.Stops and starts in mammalian oocytes：Recent advances in understanding the regulation of meiotic arrest and oocyte maturation[J].Reproduction，2005，130（6）：791-799.

[16] ALBUZ F K，SASSEVILLE M，LANE M，et al.Simulated physiological oocyte maturation（SPOM）：A novel in vitro maturation system that substantially improves embryo yield and pregnancy outcomes[J].Hum Reprod，2010，25（12）：2999-3011.

[17] TIAN X Z，WANG F，ZHANG L，et al.Beneficial effects of melatonin on the in vitro maturation of sheep oocytes and its relation to melatonin receptors[J].Int J Mol Sci，2017，18（4）：1-15.

[18] CHANG H M，KLAUSEN C，LEUNG P C K.Antimüllerian hormone inhibits follicle-stimulating hormone-induced adenylyl cyclase activation，aromatase expression，and estradiol production in human granulosa-lutein cells[J].Fertil Steril，2013，100（2）：585-592.

[19]? GROSSMAN M P，NAKAJIMA S T，FALLAT M E，et al.Müllerian-inhibiting substance inhibits cytochrome P450 aromatase activity in human granulosa lutein cells culture[J].Fertil Steril，2008，89（S5）：1364-1370.

[20]? CROSBY I M，MOOR R M，HESLOP J P，et al.cAMP in ovine oocytes：Localization of synthesis and its action on protein synthesis，phosphorylation，and meiosis[J].J Exp Zool，1985，234（2）：307-318.