畜牧兽医学报  2021, Vol. 52 Issue (9): 2452-2463. DOI: 10.11843/j.issn.0366-6964.2021.09.008    PDF    
牦牛StAR基因克隆及其生物信息学与组织表达特性的研究
刘宇, 邬建飞, 李恒, 荆天, 卢建远, 字向东     
西南民族大学 动物科学国家民委重点实验室, 成都 610041
摘要:旨在克隆获得牦牛StAR基因编码序列(CDS)并进行生物信息学分析,探究其mRNA组织表达特性。本研究以屠宰场采集的成年母牦牛心、肝、脾、肺、肾、卵巢、输卵管、子宫组织(n=5),不同年龄(胎牛、1岁、2岁)牦牛的卵巢(n=3),不同发情周期(卵泡期、黄体期)的牦牛卵巢(n=3),黄体期黄牛的卵巢(n=3)及实验室冻存的牦牛颗粒细胞为研究材料。以牦牛黄体期卵巢cDNA为模板,用逆转录PCR克隆StAR基因,并使用MEGA7.0和ExPASy-ProtParam等软件分析其生物信息学特性;采用实时荧光定量PCR技术分析牦牛StAR基因组织表达特性。结果发现,StAR基因CDS区长858 bp,编码285个氨基酸,StAR蛋白总体带正电荷,属于碱性亲水稳定蛋白,无跨膜结构及信号肽,主要存在于细胞质和线粒体;StAR基因具有较高的保守性,符合物种进化规律。牦牛StAR基因在卵巢表达水平最高(P < 0.01),且2岁时卵巢表达水平极显著高于胎牛和1岁龄(P < 0.01),黄体期卵巢表达水平极显著高于卵泡期(P < 0.01);黄体期黄牛卵巢中StAR基因的表达量极显著高于牦牛(P < 0.01);在颗粒细胞的体外培养过程中StAR基因表达量逐渐上升,在培养24 h时达到高峰(P < 0.01),随后显著降低。综上所述,StAR基因序列较为保守,在牦牛卵巢组织中表达最高,且表达水平随年龄与卵巢周期而变化,提示StAR基因可能参与牦牛卵巢及黄体功能相关的繁殖调控。
关键词牦牛    StAR基因    基因克隆    组织表达    卵巢    
Cloning and Analysis of Bioinformatics and Tissue Expression Characteristics of Yak StAR Gene
LIU Yu, WU Jianfei, LI Heng, JING Tian, LU Jianyuan, ZI Xiangdong     
Key Laboratory of Animal Science of State Ethnic Affairs Commission, Southwest Minzu University, Chengdu 610041, China
Abstract: The objective of this study was to obtain the coding sequence (CDS) of yak StAR gene, conduct bioinformatics analysis and investigate its mRNA expression characteristics in different tissues. Yak samples were collected from local slaughter house, including heart, liver, spleen, lung, kidney, ovary, oviduct and uterus of adult female yaks (n=5), ovary of yaks at different ages (fetal, 1-year-old and 2-years-old) (n=3), ovary of yaks during different periods of estrous cycles (follicular phase and luteal phase) (n=3), the luteal phase ovary of cattle (n=3). Granulosa cells of yaks were cryopreserved previously. StAR gene was cloned using yak ovarian cDNA as template during luteal phase by reverse transcriptase-polymerase chain reaction (PCR) and its bioinformatics properties were analyzed by bioinformatics software such as MEGA7.0 and ExPASy-ProtParam etc. The tissue expression characteristics of yak StAR gene were detected by quantitative real-time PCR. The results showed that the CDS of yak StAR gene was 858 bp, encoding 285 amino acids. StAR protein was positively charged and belonged to alkaline hydrophilic stable protein, which had no transmembrane structure and signal peptide. It mainly existed in cytoplasm and mitochondria. StAR gene was highly conserved between different species, which accorded with the law of species evolution. The expression level of yak StAR gene was the highest in ovary (P < 0.01). The expression level of StAR gene in the ovary of 2-years-old was significantly higher than that of fetal and 1-year-old (P < 0.01). Its expression in yak ovary during luteal phase was significantly higher than that during follicular phase (P < 0.01), and was significantly higher in cattle ovary than in yak ovary during luteal phase(P < 0.01). During the in vitro culture (IVC) of granulosa cells, the expression of StAR gene increased gradually and reached the peak level at 24 h IVC (P < 0.01), and then significantly decreased. The results showed that StAR gene was relatively conservative in the course of animal evolution. Its expression level was greatest in the ovary of yak, and changing with age- and estrous cycle-specific characteristics, which indicate that StAR gene may be involved in the regulation of yak reproduction associated with ovary and corpus luteum function.
Key words: yak    StAR gene    gene cloning    tissue expression    ovary    

Clark等[1]1994年从小鼠MA-10间质瘤细胞中首次发现类固醇激素合成急性调节蛋白(steroidogenic acute regulatory protein, StAR),并成功地完成了其纯化、前体cDNA的克隆和测序,证明了StAR蛋白对类固醇生物合成有重要调控作用。StAR存在于线粒体中,促进胆固醇由线粒体外膜向内膜转运,胆固醇的这种转运是类固醇合成的限速步骤,胆固醇进入内膜后再被CYP11A1基因编码的细胞色素P450支链裂解酶(P450 scc)催化生成孕烯醇酮[2]。孕烯醇酮作为类固醇激素的合成前体,其离开线粒体后,在不同组织中经过一系列的变化能够合成孕酮、睾酮、雌激素、醛固酮等类固醇激素。StAR基因在肾上腺皮质、睾丸间质细胞、卵巢颗粒细胞和膜细胞以及胎盘合体滋养细胞均有表达[3]StAR基因的突变会导致先天性类脂肾上腺增生,患者肾上腺增大,大量的胆固醇无法转化为孕烯醇酮[4-5]。StAR蛋白的表达水平强烈影响睾丸间质细胞的睾酮生成,而且随着间质细胞老化,StAR蛋白的表达也会减少,衰老间质细胞线粒体内膜的胆固醇供应随之减少[6-7]。卵巢合成的类固醇激素是维持雌性第二性特征和生殖能力所必需的,有研究发现,StAR基因敲除小鼠的卵巢在出生时没有受到影响,但随着小鼠进入青春期,卵巢中的脂质继续积累,导致雌激素水平低,出现乳房发育不良和阴道出血等症状[8]StAR基因在腔前卵泡数多的胎牛卵巢中的表达量要显著高于腔前卵泡数少的卵巢[9],黔北麻羊卵巢中StAR基因在高繁殖力组的表达量显著高于低繁殖力组[10],并且StAR基因在不同繁殖力品种的山羊卵巢中表达差异极显著[11]。目前,在猪[12]、牛[13]、斑马鱼[14]、绵羊[15]和禽[16]等动物上的研究证明,StAR基因可能与性腺的发育和繁殖性状密切相关。

牦牛(Bos grunniens)主要生活在海拔2 500~5 500 m的青藏高原,能适应高原的极端恶劣气候条件,是青藏高原不可替代的主导畜种[17]。但是,与黄牛相比,牦牛繁殖性能低,属于季节性发情动物,主要集中在每年的6~10月发情,黄牛则是四季发情动物,其繁殖力远远高于牦牛[18]StAR基因可能与动物繁殖力密切相关,但目前未见牦牛StAR基因研究的报道。因此,本试验拟通过克隆牦牛StAR基因并进行生物信息学分析,利用实时荧光定量PCR技术对StAR基因在不同组织、不同生长阶段和不同发情周期的卵巢、颗粒细胞中的表达特性进行研究,为进一步研究StAR基因在牦牛繁殖中的调控作用提供理论依据。

1 材料与方法 1.1 试验动物及样品采集

本次试验的样品均采自四川省成都市青白江唐家寺屠宰场,牦牛是来自于阿坝州牧区的麦洼牦牛。电击放血屠宰后,立即采样。样品包括:1)5头5~6岁健康母牦牛心、肝、脾、肺、肾、卵巢、子宫和输卵管组织;2)牦牛胎牛、1岁龄左右和2岁龄左右黄体期卵巢各3头;3)5~6岁卵泡期和黄体期牦牛的卵巢各3头;4)3头5~6岁黄体期黄牛的卵巢;5)牦牛颗粒细胞为本实验室前期冻存的细胞。使用灭菌后的剪刀采集组织,并用0.9%的生理盐水冲洗干净,迅速放入液氮罐中,带回实验室-80 ℃保存备用。

1.2 主要试剂和仪器

TRIzol Reagent和琼脂糖购自美国Invitrogen公司;RNA提取试剂盒购自北京天漠科技开发有限公司;cDNA反转录试剂盒、2×Taq PCR Master mix和PowerUpTM SYBRTM Green均购自Thermo Scientific公司;DNA Marker DL2000购自GenStAR公司;E. coli DH5α感受态细胞和胶回收试剂盒购自成都擎科梓熙生物技术有限公司;克隆载体pMD19-T和氨苄青霉素购自宝生生物工程(大连)有限公司。

PCR仪(ETC811)购自苏州东盛兴业科学仪器有限公司;荧光定量PCR仪(CFX96)购自美国Bio Rad公司;琼脂糖凝胶成像系统(CL1000)购自Invitrogen公司。

1.3 RNA提取和cDNA合成

按RNA提取试剂盒的步骤提取RNA。使用紫外分光光度计检测RNA的质量,1%的琼脂糖凝胶电泳检测RNA完整性。按照反转录试剂盒说明书对RNA进行反转录,合成的cDNA -20 ℃保存备用。

1.4 StAR基因引物设计合成

由于未见牦牛StAR基因序列报道,所以根据NCBI中黄牛(Bos taurus)的StAR基因序列(NM_174189.3)和GAPDH 基因序列(NM_001034034.2),使用Primer Premier 5.0软件设计克隆引物以及内参基因GAPDH荧光定量引物;再根据克隆获得的序列设计StAR基因的荧光定量引物(表 1)。引物均由成都擎科梓熙生物技术有限公司合成。

表 1 PCR和荧光定量PCR引物序列 Table 1 Primer sequences for PCR and quantitative real-time PCR
1.5 牦牛StAR基因克隆与测序

采用卵巢合成的cDNA为模板进行PCR扩增,克隆StAR基因。PCR扩增体系25 μL:2×Taq PCR MasterMix 12.5 μL,上、下游引物各1 μL,cDNA 1 μL,ddH2O 9.5 μL。PCR反应程序:94 ℃预变性2 min;94 ℃变性30 s,退火30 s(退火温度见表 1),72 ℃延伸1 min,共37个循环;72 ℃延伸5 min;4 ℃保存。

PCR产物用1%琼脂糖凝胶电泳检测,按照胶回收试剂盒说明书将胶回收和纯化。然后将回收产物与pMD19-T Vector连接并在16 ℃过夜,后转化DH5α感受态细胞,均匀涂于LB固体培养基(AMP+)上,挑取经37 ℃培养24 h过后呈白色的单菌落,接种于LB液体培养基(AMP+)中摇床6 h。菌液进行PCR鉴定,将鉴定结果呈阳性的样品菌液装袋送至成都擎科梓熙生物技术有限公司测序。

1.6 StAR基因生物信息学分析

牦牛StAR基因生物学分析工具及功能参见表 2

表 2 生物信息学分析工具及相关功能 Table 2 Bioinformatics analysis tools and corresponding functions
1.7 牦牛StAR基因表达特性分析

根据克隆获得的StAR基因设计的引物(表 1),使用qPCR方法检测StAR基因在各个组织(心、肝、脾、肺、肾、卵巢、输卵管、子宫)及颗粒细胞中的表达水平。颗粒细胞由实验室前期冷冻保存,将其复苏并参考洪磊等[20]的方法进行体外培养(in vitro culture,IVC),分别收集不同时期(0、12、24、36、48、60 h)的颗粒细胞。每个样品重复3次,用GAPDH作为内参基因校正基因的相对表达水平。qPCR反应体系15 μL:PowerUpTM SYBRTM Green MasterMix 7.5 μL,cDNA 1 μL,上、下游引物各0.5 μL,ddH2O 5.5 μL。反应条件:95 ℃预变性3 min;95 ℃变性10 s,60 ℃退火30 s,39个循环;读取熔解曲线。

1.8 统计分析

本试验用2-△△Ct法计算基因相对表达量,使用SPSS 21.0软件的独立样本t检验和单因素方差分析法进行显著性检验,当P < 0.05时,表示差异有统计学意义。

2 结果 2.1 牦牛StAR基因克隆

以卵巢cDNA为模板,PCR扩增后,产物经1%的琼脂糖凝胶电泳检测,获得了与预期目的产物大小相符的片段(图 1),表明本试验成功的扩增了牦牛StAR基因。测序得到StAR基因长度为1 156 bp,其中CDS长度为858 bp,共编码285个氨基酸。序列提交至NCBI,获得登录号:MW495047。

M. DNA相对分子质量标准; 1~3. StAR引物PCR扩增产物 M. DL2000 DNA Marker; 1-3. PCR amplification products of StAR primers 图 1 牦牛StAR基因的PCR扩增结果 Fig. 1 PCR amplification result of yak StAR gene
2.2 牦牛StAR基因生物信息学分析

2.2.1 牦牛StAR基因核苷酸序列分析   使用DNAMAN 6.0软件将测定序列的CDS与黄牛StAR基因的CDS进行比较,发现了8个碱基突变,其中有3个导致了氨基酸的改变(表 3)。牦牛与黄牛的StAR基因CDS同源性为99.07%。

表 3 牦牛与黄牛StAR基因编码区序列差异 Table 3 Coding sequence difference of StAR gene between yak and cattle

2.2.2 StAR基因编码氨基酸序列同源性分析   使用DNAMAN 6.0软件对不同物种StAR蛋白氨基酸序列进行比对,结果显示,牦牛与黄牛(NP_776614.2)、绵羊(NP_001009243.1)、人(NP_ 000340.2)、猪(NP_998920.2)、原鸡(NP_990017.1)、斑马鱼(NP_571738.1)的StAR氨基酸序列相似性分别为98.95%、96.84%、87.37%、90.18%、65.96%、58.54%,说明StAR基因在进化过程中有较高的保守性(图 2)。

黑色.同源性=100%;粉色.同源性>75%;蓝色.同源性>50%;白色.同源性 < 50% Black. Homology=100%;Pink. Homology>75%;Blue. Homology>50%;White. Homology < 50% 图 2 牦牛和其他物种间StAR氨基酸序列比对 Fig. 2 Alignment of amino acid sequence of StAR between yak and other species

2.2.3 牦牛StAR蛋白理化性质分析   使用ExPASy-ProtParam在线工具对牦牛的StAR蛋白进行分析,发现其蛋白分子式为C1398H2283N409O407S16,分子质量为31 845.96 u;StAR蛋白的氨基酸组成中亮氨酸(Leu)、丙氨酸(Ala)、谷氨酸(Glu)所占比例较高, 分别为11.9%、8.4%、8.1%;带正电荷的(Arg+Lys) 氨基酸总数为38,带负电荷的(Asp+Glu)的氨基酸总数为32,说明牦牛StAR蛋白带正电荷;理论等电点为9,表明StAR蛋白是碱性蛋白;在哺乳动物的红细胞中半衰期为30 h;脂肪族系数为93.75;平均亲水系数为-0.288,不稳定系数为39.99,表明StAR蛋白属于亲水稳定蛋白。

根据TMHMM2.0在线工具分析发现,StAR蛋白没有跨膜结构;使用Signal P 5.0预测,结果表明, StAR蛋白不存在信号肽;使用PSORTII在线工具分析发现,StAR蛋白主要存在于细胞质(47.8%),也存在于线粒体(30.4%)、细胞核(13%)、细胞骨架(4.3%)、内质网(4.3%)。使用Net Phos 3.1在线工具分析,发现牦牛StAR蛋白有22个磷酸化位点,其中丝氨酸(Ser)磷酸化位点12个、苏氨酸(Thr)磷酸化位点6个、酪氨酸(Tyr)磷酸化位点4个。

2.2.4 StAR蛋白二、三级结构预测   使用SOPMA在线工具预测StAR蛋白二级结构,结果显示,α-螺旋占比为39.3%、延伸链占比为17.54%、β-转角占比为4.21%、无规则卷曲占比为38.95%(图 3)。

h. α-螺旋;e. 延伸链;t. β-转角;c. 无规卷曲 h. α-helix; e. Extended chain; t. β-turn; c. Random coil 图 3 牦牛StAR蛋白二级结构预测 Fig. 3 Secondary structure prediction of yak StAR protein

使用SWISS-MODEL在线工具预测StAR蛋白三级结构,预测结果与二级结构预测结果基本一致(图 4)。

图 4 牦牛StAR蛋白三级结构预测 Fig. 4 The tertiary structure prediction of yak StAR protein

2.2.5 蛋白相互作用分析   根据STRING11.0在线工具分析发现,StAR蛋白可能与电压依赖性阴离子通道蛋白1(voltage dependent anion channel protein 1, VDAC1)、VDAC2、VDAC3、FIGLA、NR5A1、HSD17B6、POMC等蛋白有相互作用(图 5)。

图 5 StAR蛋白网络互作分析图 Fig. 5 StAR protein network interaction analysis

2.2.6 牦牛StAR基因系统进化树分析   从NCBI下载黄牛(NM_174189.3)、绵羊(NM_ 001009243.1)、山羊(XM_013975437.2)、小鼠(NM_ 011485.5)、大鼠(NM_031558.3)、猪(NM_213755.2)、马(NM_001081800.3)、狗(NM_001097542.1)、人(NM_000349.3)、斑马鱼(NM_131663.1)、原鸡(NM_204686.2)和黑猩猩(XM_001170375.5)12个物种的StAR基因序列。利用Mega 6.0构建基因进化树,发现牦牛与黄牛首先聚为一类,其次与绵羊、山羊聚为一类,与斑马鱼亲缘关系最远,进化比较保守,且符合物种进化规律(图 6)。

图 6 StAR基因系统进化树 Fig. 6 Phylogenetic tree of StAR gene
2.3 牦牛StAR基因表达差异分析

2.3.1 牦牛StAR基因在各个组织中的表达差异分析   使用qPCR检测StAR基因在心、肝、脾、肺、肾、卵巢、输卵管、子宫组织中的表达水平,结果显示,StAR基因在卵巢中高表达,极显著高于其他组织(P < 0.01),在其他组织中表达均较低(图 7)。

不同大写字母表示差异极显著(P < 0.01)。下同 Different capital letters represent the extremely significant differences (P < 0.01). The same as below 图 7 StAR基因在不同组织中的相对表达水平 Fig. 7 Relative expression level of StAR gene in different tissues

2.3.2 StAR基因在卵巢的表达分析   使用qPCR检测StAR基因在胎牛、1岁左右、2岁左右(黄体期)牦牛卵巢的表达情况,结果显示,随着年龄的增大,StAR基因在卵巢中的相对表达量逐渐增加,各个时期之间差异极显著(P < 0.01,图 8A)。StAR基因在黄体期卵巢的表达要极显著高于卵泡期(P < 0.01,图 8B),在黄牛黄体期卵巢的表达要极显著高于在牦牛黄体期卵巢的表达(P < 0.01,图 8C)。

图 8 卵巢StAR基因的相对表达水平 Fig. 8 Relative expression level of StAR in ovary

2.3.3 StAR基因在颗粒细胞中的时序表达分析   利用qPCR检测StAR基因在体外培养颗粒细胞不同阶段的相对表达水平,结果发现,StAR基因的相对表达水平呈先上升后下降的趋势,在24 h达到最高,极显著高于其他时间段(P < 0.01,图 9)。

图 9 StAR基因在颗粒细胞不同时序的相对表达水平 Fig. 9 Relative expression level of StAR gene in granulosa cells at different time
3 讨论

StAR是类固醇生物合成的限速步骤,StAR蛋白在肾上腺、卵巢、睾丸、胎盘、脑的类固醇生成细胞中合成,影响一系列发育和生理过程,特别是在繁殖相关性状中具有重要的调控作用[4, 21-23]StAR基因敲除鼠在卵巢上会出现脂质沉积聚集、基质细胞黄素化和卵泡的不完全成熟,最终导致卵巢早衰[8]。这些研究说明StAR基因在繁殖性状上有重要的作用。

为了探究牦牛StAR基因与繁殖性状的相关性及与黄牛StAR基因的差异性,本试验克隆了牦牛StAR基因,并获得了其完整的CDS,发现CDS全长858 bp,编码285个氨基酸。本研究发现,牦牛StAR蛋白共有22个磷酸化位点,早期的研究也表明,仓鼠StAR是一种磷酸化蛋白[24]。Artemenko等[25]还发现,StAR蛋白存在多种形式,StAR首先被合成为一个37 ku的前体,它被蛋白激酶A磷酸化,并在线粒体内被切割成30 ku大小的蛋白,并且磷酸化后的蛋白提供了主要的胆固醇转运机制。本研究发现,StAR蛋白主要存在于细胞质(47.8%)和线粒体(30.4%),这与Bose等[26]发现的StAR蛋白具有从线粒体外膜向内膜转运胆固醇的功能相符合。牦牛StAR蛋白氨基酸同源性比对和基因进化树分析结果显示,StAR基因与黄牛、绵羊、山羊具有较高的同源性,在进化上是比较保守的。Bauer等[27]对鱼类、两栖动物、鸟类和哺乳动物StAR的氨基酸序列进行了比较,也发现这种蛋白质在不同的脊椎动物群体中具有很强的保守性,与本研究结果一致。StAR基因的突变会导致类固醇缺乏症,影响睾丸和卵巢功能[5, 8, 23]。本研究发现。牦牛StAR基因的编码区与黄牛相比有8个碱基突变,其中有3个导致了氨基酸的改变,这种差异是否影响到牦牛StAR蛋白的功能,从而影响牦牛的繁殖性能需要进一步研究。另外有研究发现,在三级结构当中StAR蛋白存在一个StART结构域,在本研究结果的三级结构图中也能明显看到了该结构域,该结构是α/β类型的,含有1个九链反平行的β片层和4个α-螺旋,形成的长疏水通道可以结合1个胆固醇分子[28-29]

根据蛋白互作网络图显示,StAR蛋白可能与VDAC1、VDAC2、VDAC3、FIGLA、NR5A1、HSD17B6、POMC等蛋白相互作用。VDAC是线粒体孔道蛋白,参与代谢物/离子的跨膜转运和许多细胞过程[30]。Bose等[31]发现,在肾上腺细胞中磷酸化的StAR与线粒体外膜上的电压依赖阴离子通道1(voltage-dependent anion channel 1, VDAC1)相互作用,从而促进了37 ku磷酸化StAR到32 ku中间体的加工。生殖细胞因子(factor in the germline alpha, FIGLA)是一种在生殖细胞中特异表达的转录因子,可以调控卵泡发育和编码透明带的相关基因,FIGLA基因的缺失明显阻碍了减数分裂的进展,并最终导致卵母细胞凋亡,另外干扰睾丸细胞中的FIGLA基因后,会上调StAR基因的表达[32-34]。核受体亚家族5组A成员1(nuclear receptor subfamily 5 group A member 1, NR5A1)基因编码的类固醇生成因子1(steroidogenic factor 1,SF-1)是一种转录激活剂,在调节正常生殖生理和内分泌功能上具有重要作用, 控制肾上腺和性腺发育[35]。Sandhoff等[36]发现,大鼠StAR基因的启动子区域存在一个SF-1结合位点,对StAR基因的转录是至关重要。17β羟基类固醇脱氢酶6型(hydroxysteroid 17-beta dehydrogenase 6, HSD17B6)基因编码的蛋白质既具有氧化还原酶活性,又具有表异构酶活性,参与雄激素分解代谢,有研究发现,大鼠间质细胞在化疗药物的作用下会促进StAR基因的转录,抑制HSD17B6的转录,并伴随睾酮水平降低[37-38]。阿片黑素促皮质激素原(pro-opiomelanocortin, POMC)基因编码一种前原蛋白,可以产生多达10种的生物活性肽,促肾上腺皮质激素(ACTH)是其中一种,有研究发现,ACTH可以通过cAMP/PKA途径调控类固醇合成,促进StAR基因的表达[39]。从蛋白互作网络的分析结果可以发现,与StAR蛋白存在相互作用的蛋白主要参与类固醇激素合成和性腺发育的调控,猜测StAR基因与繁殖相关性状有很强的关联性。

对多个物种的研究表明,动物StAR基因的正常表达与繁殖性能密切相关[12-16]。本试验采用qPCR检测StAR基因的组织表达情况,发现StAR基因在牦牛的卵巢中表达量最高,有研究也发现,StAR基因主要在雌性动物的卵巢中表达,比如猪[12]、鸡[40]、斑马鱼[14]、鲈鱼[41]等。这些研究结果与StAR蛋白的功能及卵巢是类固醇激素主要来源的内分泌特性相吻合。牦牛StAR基因在脾、肺、输卵管、子宫有极其微量的表达,而在其他组织中不表达,但是刘倩等[42]在中华鳖中发现,StAR基因在脾脏和心脏中高表达,欧阳宏佳等[16]在水禽中检测到StAR基因在胸肌和腿肌中表达,说明该基因的表达具有一定的种属特异性。本研究发现,StAR基因在2岁左右牦牛卵巢表达最高,极显著高于胎牛和1岁左右的牦牛。欧阳宏佳等[16]发现,鹅睾丸中StAR基因的表达量在3岁时极显著高于52和1日龄,并且Michalovic等[43]发现,奶牛卵巢颗粒细胞中StAR基因的表达量在成年时极显著高于初情期前,说明随着年龄的增长,性腺逐渐发育完全其功能逐渐完善,睾丸与卵巢合成与分泌类固醇激素的能力也增强。StAR基因在牦牛黄体期的表达量极显著高于卵泡期,Jahromi等[44]在大鼠中也发现,StAR基因在发情间期(黄体期)的表达水平极显著高于其他时期,说明黄体期卵巢孕酮等类固醇激素的合成分泌增加与StAR基因的高表达有关。本研究发现,在黄牛卵巢黄体期StAR基因的表达量极显著高于牦牛。Prakash等[45]发现,在发情周期的8~12天(黄体期)黄牛的血浆孕酮水平显著高于牦牛,而毛宁等[46]发现,高产香猪组卵巢StAR基因表达量显著高于低产组,杨沛方等[10]报道,山羊在高繁殖力组卵巢中StAR基因的表达量显著高于低繁殖力组,推测StAR基因在牦牛和黄牛卵巢表达的差异性可能是导致两物种繁殖力差异的重要原因之一。StAR基因在牦牛颗粒细胞IVC过程中的相对表达水平呈先上升后下降的趋势,在培养24 h表达量最高(图 9),这与牦牛卵母细胞体外培养18~24 h时达到成熟的结果一致[47-49],说明颗粒细胞StAR基因的表达与卵母细胞成熟有一定的相关性,对卵母细胞成熟具有重要的作用。

4 结论

本研究克隆得到牦牛StAR基因CDS全长858 bp,编码285个氨基酸,主要在线粒体膜上发挥作用,在进化中和物种间比较保守。牦牛StAR基因主要在卵巢中表达,而且黄体期表达高于卵泡期,并且在体外颗粒细胞培养至24 h表达水平最高。提示,StAR基因可能在牦牛卵巢周期中参与繁殖调控。

参考文献
[1]
CLARK B J, WELLS J, KING S R, et al. The purification, cloning, and expression of a novel luteinizing hormone-induced mitochondrial protein in MA-10 mouse Leydig tumor cells: characterization of the steroidogenic acute regulatory protein (StAR)[J]. J Biol Chem, 1994, 269(45): 28314-28322. DOI:10.1016/S0021-9258(18)46930-X
[2]
TUCKEY R C, HEADLAM M J, BOSE H S, et al. Transfer of cholesterol between phospholipid vesicles mediated by the steroidogenic acute regulatory protein (StAR)[J]. J Biol Chem, 2002, 277(49): 47123-47128. DOI:10.1074/jbc.M206965200
[3]
MILLER W L, BOSE H S. Early steps in steroidogenesis: intracellular cholesterol trafficking[J]. J Lipid Res, 2011, 52(12): 2111-2135. DOI:10.1194/jlr.R016675
[4]
LIN D, SUGAWARA T, STRAUSS Ⅲ J F, et al. Role of steroidogenic acute regulatory protein in adrenal and gonadal steroidogenesis[J]. Science, 1995, 267(5205): 1828-1831. DOI:10.1126/science.7892608
[5]
KANG E, KIM Y M, KIM G H, et al. Mutation spectrum of STAR and the founder effect of p.Q258* in Korean patients with congenital lipoid adrenal hyperplasia[J]. Mol Med, 2017, 23: 149-154.
[6]
LEERS-SUCHETA S, STOCCO D M, AZHAR S. Down-regulation of steroidogenic acute regulatory (StAR) protein in rat Leydig cells: implications for regulation of testosterone production during aging[J]. Mech Ageing Dev, 1999, 107(2): 197-203. DOI:10.1016/S0047-6374(98)00149-3
[7]
STOCCO D M. StAR protein and the regulation of steroid hormone biosynthesis[J]. Annu Rev Physiol, 2001, 63: 193-213. DOI:10.1146/annurev.physiol.63.1.193
[8]
HASEGAWA T, ZHAO L P, CARON K M, et al. Developmental roles of the steroidogenic acute regulatory protein (StAR) as revealed by StAR knockout mice[J]. Mol Endocrinol, 2000, 14(9): 1462-1471. DOI:10.1210/mend.14.9.0515
[9]
WELLER M M D C A, FORTES M R S, MARCONDES M I, et al. Effect of maternal nutrition and days of gestation on pituitary gland and gonadal gene expression in cattle[J]. J Dairy Sci, 2016, 99(4): 3056-3071. DOI:10.3168/jds.2015-9673
[10]
杨沛方, 周志楠, 敖叶, 等. 黔北麻羊StAR基因的克隆、生物信息学分析及在单、多羔黔北麻羊不同性腺组织中的表达[J]. 中国畜牧杂志, 2020, 56(11): 78-84.
YANG P F, ZHOU Z N, AO Y, et al. Cloning, bioinformatics analysis and expression of StAR gene in Qianbei brown goat[J]. Chinese Journal of Animal Science, 2020, 56(11): 78-84. (in Chinese)
[11]
ZI X D, LU J Y, ZHOU H, et al. Comparative analysis of ovarian transcriptomes between prolific and non-prolific goat breeds via high-throughput sequencing[J]. Reprod Domest Anim, 2018, 53(2): 344-351. DOI:10.1111/rda.13111
[12]
PILON N, DANEAU I, BRISSON C, et al. Porcine and bovine steroidogenic acute regulatory protein (StAR) gene expression during gestation[J]. Endocrinology, 1997, 138(3): 1085-1091. DOI:10.1210/endo.138.3.5003
[13]
HARTUNG S, RUST W, BALVERS M, et al. Molecular cloning and in vivo expression of the bovine steroidogenic acute regulatory protein[J]. Biochem Biophys Res Commun, 1995, 215(2): 646-653. DOI:10.1006/bbrc.1995.2513
[14]
SHANG G H, PENG X Y, JI C, et al. Steroidogenic acute regulatory protein and luteinizing hormone are required for normal ovarian steroidogenesis and oocyte maturation in zebrafish[J]. Biol Reprod, 2019, 101(4): 760-770. DOI:10.1093/biolre/ioz132
[15]
JUENGEL J L, MEBERG B M, TURZILLO A M, et al. Hormonal regulation of messenger ribonucleic acid encoding steroidogenic acute regulatory protein in ovine corpora lutea[J]. Endocrinology, 1995, 136(12): 5423-5429. DOI:10.1210/endo.136.12.7588291
[16]
欧阳宏佳, 孙敬帅, 江丹莉, 等. 水禽StAR基因克隆、表达及其对睾丸发育的影响[J]. 畜牧兽医学报, 2020, 51(12): 3013-3022.
OUYANG H J, SUN J S, JIANG D L, et al. Cloning and expression of waterfowl StAR gene and its effect on testicular development[J]. Acta Veterinaria et Zootechnica Sinica, 2020, 51(12): 3013-3022. DOI:10.11843/j.issn.0366-6964.2020.12.011 (in Chinese)
[17]
QIU Q, ZHANG G J, MA T, et al. The yak genome and adaptation to life at high altitude[J]. Nat Genet, 2012, 44(8): 946-949. DOI:10.1038/ng.2343
[18]
ZI X D. Reproduction in female yaks (Bos grunniens) and opportunities for improvement[J]. Theriogenology, 2003, 59(5-6): 1303-1312. DOI:10.1016/S0093-691X(02)01172-X
[19]
KUMAR S, STECHER G, TAMURA K. MEGA7:molecular evolutionary genetics analysis version 7.0 for bigger datasets[J]. Mol Biol Evol, 2016, 33(7): 1870-1874. DOI:10.1093/molbev/msw054
[20]
洪磊, 陈祥, 唐文, 等. 黔北麻羊SRD5A2基因干扰载体构建及其对产羔相关基因表达的影响[J]. 畜牧兽医学报, 2020, 51(12): 2980-2990.
HONG L, CHEN X, TANG W, et al. Construction of SRD5A2 gene interference vector and its effect on the expression of lambing related genes in Qianbei Ma goats[J]. Acta Veterinaria et Zootechnica Sinica, 2020, 51(12): 2980-2990. DOI:10.11843/j.issn.0366-6964.2020.12.008 (in Chinese)
[21]
MANNA P R, STETSON C L, SLOMINSKI A T, et al. Role of the steroidogenic acute regulatory protein in health and disease[J]. Endocrine, 2016, 51(1): 7-21. DOI:10.1007/s12020-015-0715-6
[22]
MILLER W L. Steroidogenic acute regulatory protein (StAR), a novel mitochondrial cholesterol transporter[J]. Biochim Biophys Acta, 2007, 1771(6): 663-676. DOI:10.1016/j.bbalip.2007.02.012
[23]
CHEN H, ZHANG Q R, CHEN R M, et al. Lipoid congenital adrenal hyperplasia due to steroid acute regulatory protein (STAR) variants in Three Chinese patients[J]. J Steroid Biochem Mol Biol, 2020, 200: 105635. DOI:10.1016/j.jsbmb.2020.105635
[24]
FLEURY A, MATHIEU A P, DUCHARME L, et al. Phosphorylation and function of the hamster adrenal steroidogenic acute regulatory protein (StAR)[J]. J Steroid Biochem Mol Biol, 2004, 91(4-5): 259-271. DOI:10.1016/j.jsbmb.2004.04.010
[25]
ARTEMENKO I P, ZHAO D, HALES D B, et al. Mitochondrial processing of newly synthesized steroidogenic acute regulatory protein (StAR), but not total StAR, mediates cholesterol transfer to cytochrome P450 side chain cleavage enzyme in adrenal cells[J]. J Biol Chem, 2001, 276(49): 46583-46596. DOI:10.1074/jbc.M107815200
[26]
BOSE H S, LINGAPPA V R, MILLER W L. Rapid regulation of steroidogenesis by mitochondrial protein import[J]. Nature, 2002, 417(6884): 87-91. DOI:10.1038/417087a
[27]
BAUER M P, BRIDGHAM J T, LANGENAU D M, et al. Conservation of steroidogenic acute regulatory (StAR) protein structure and expression in vertebrates[J]. Mol Cell Endocrinol, 2000, 168(1-2): 119-125. DOI:10.1016/S0303-7207(00)00316-6
[28]
TSUJISHITA Y, HURLEY J H. Structure and lipid transport mechanism of a StAR-related domain[J]. Nat Struct Biol, 2000, 7(5): 408-414. DOI:10.1038/75192
[29]
SCHRICK K, BRUNO M, KHOSLA A, et al. Shared functions of plant and mammalian StAR-related lipid transfer (START) domains in modulating transcription factor activity[J]. BMC Biol, 2014, 12: 70. DOI:10.1186/s12915-014-0070-8
[30]
CATERINO M, RUOPPOLO M, MANDOLA A, et al. Protein-protein interaction networks as a new perspective to evaluate distinct functional roles of voltage-dependent anion channel isoforms[J]. Mol Biosyst, 2017, 13(12): 2466-2476. DOI:10.1039/C7MB00434F
[31]
BOSE M, WHITTAL R M, MILLER W L, et al. Steroidogenic activity of StAR requires contact with mitochondrial VDAC1 and phosphate carrier protein[J]. J Biol Chem, 2008, 283(14): 8837-8845. DOI:10.1074/jbc.M709221200
[32]
BAYNE R A L, DA SILVA S J M, ANDERSON R A. Increased expression of the FIGLA transcription factor is associated with primordial follicle formation in the human fetal ovary[J]. Mol Hum Reprod, 2004, 10(6): 373-381. DOI:10.1093/molehr/gah056
[33]
WANG Z P, LIU C Y, ZHAO Y G, et al. FIGLA, LHX8 and SOHLH1 transcription factor networks regulate mouse oocyte growth and differentiation[J]. Nucleic Acids Res, 2020, 48(7): 3525-3541. DOI:10.1093/nar/gkaa101
[34]
LI H L, XU W T, ZHANG N, et al. Two Figla homologues have disparate functions during sex differentiation in half-smooth tongue sole (Cynoglossus semilaevis)[J]. Sci Rep, 2016, 6: 28219. DOI:10.1038/srep28219
[35]
NA X X, MAO Y, TANG Y M, et al. Identification and functional analysis of fourteen NR5A1 variants in patients with the 46 XY disorders of sex development[J]. Gene, 2020, 760: 145004. DOI:10.1016/j.gene.2020.145004
[36]
SANDHOFF T W, HALES D B, HALES K H, et al. Transcriptional regulation of the rat steroidogenic acute regulatory protein gene by steroidogenic factor 1[J]. Endocrinology, 1998, 139(12): 4820-4831. DOI:10.1210/endo.139.12.6345
[37]
MUTHUSAMY S, ANDERSSON S, KIM H J, et al. Estrogen receptorβ and 17β-hydroxysteroid dehydrogenase type 6, a growth regulatory pathway that is lost in prostate cancer[J]. Proc Natl Acad Sci U S A, 2011, 108(50): 20090-20094. DOI:10.1073/pnas.1117772108
[38]
AL-BADER M, KILARKAJE N. Effects of bleomycin, etoposide and cisplatin treatment on Leydig cell structure and transcription of steroidogenic enzymes in rat testis[J]. Eur J Pharmacol, 2015, 747: 150-159. DOI:10.1016/j.ejphar.2014.12.006
[39]
CLARK B J. ACTH action on StAR biology[J]. Front Neurosci, 2016, 10: 547.
[40]
SECHMAN A, GRZEGORZEWSKA A K, GRZESIAK M, et al. Nitrophenols suppress steroidogenesis in prehierarchical chicken ovarian follicles by targeting STAR, HSD3B1, and CYP19A1 and downregulating LH and estrogen receptor expression[J]. Domest Anim Endocrinol, 2020, 70: 106378. DOI:10.1016/j.domaniend.2019.07.006
[41]
PRUCHA M S, MARTYNIUK C J, DOPERALSKI N J, et al. Steroidogenic acute regulatory protein transcription is regulated by estrogen receptor signaling in largemouth bass ovary[J]. Gen Comp Endocrinol, 2020, 286: 113300. DOI:10.1016/j.ygcen.2019.113300
[42]
刘倩, 马晓, 岑双双, 等. 中华鳖StAR基因的克隆、表达及多克隆抗体制备[J]. 水产学报, 2020, 44(11): 1813-1823.
LIU Q, MA X, CEN S S, et al. Molecular cloning, expression analysis and polyclonal antibody preparation of StAR gene in Chinese soft-shelled turtle (Pelodiscus sinensis)[J]. Journal of Fisheries of China, 2020, 44(11): 1813-1823. (in Chinese)
[43]
MICHALOVIC L, CURRIN L, GUTIERREZ K, et al. Granulosa cells of prepubertal cattle respond to gonadotropin signaling and upregulate genes that promote follicular growth and prevent cell apoptosis[J]. Mol Reprod Dev, 2018, 85(12): 909-920. DOI:10.1002/mrd.23066
[44]
JAHROMI M S, TEHRANI F R, NOROOZZADEH M, et al. Elevated expression of steroidogenesis pathway genes; CYP17, GATA6 and StAR in prenatally androgenized rats[J]. Gene, 2016, 593(1): 167-171. DOI:10.1016/j.gene.2016.07.067
[45]
PRAKASH B S, SARKAR M, PAUL V, et al. Postpartum endocrinology and prospects for fertility improvement in the lactating riverine buffalo (Bubalus bubalis) and yak (Poephagus grunniens L.)[J]. Livest Prod Sci, 2005, 98(1-2): 13-23. DOI:10.1016/j.livprodsci.2005.10.014
[46]
毛宁, 王嘉福, 张福平, 等. 香猪卵巢StARCYP11A1基因的差异表达研究[J]. 中国畜牧兽医, 2018, 45(5): 1137-1144.
MAO N, WANG J F, ZHANG F P, et al. Study on differential expression of StAR and CYP11A1 genes in Xiang Pig ovary[J]. China Animal Husbandry & Veterinary Medicine, 2018, 45(5): 1137-1144. (in Chinese)
[47]
兰道亮, 熊显荣, 陈亚冰, 等. 牦牛体外成熟卵母细胞差异转录组学研究[J]. 中国科学: 生命科学, 2017, 47(10): 1099-1112.
LAN D L, XIONG X R, CHEN Y B, et al. Differential transcriptome analysis of yak oocytes in vitro maturation[J]. Scientia Sinica Vitae, 2017, 47(10): 1099-1112. (in Chinese)
[48]
何翃闳, 潘阳阳, 张慧珠, 等. FSH对牦牛卵母细胞EGFEGFR表达及其细胞凋亡的影响[J]. 畜牧兽医学报, 2018, 49(9): 1899-1907.
HE H H, PAN Y Y, ZHANG H Z, et al. The effects of follicle-stimulating hormone (FSH) on the expression of EGF and EGFR in yak oocytes and apoptosis[J]. Acta Veterinaria et Zootechnica Sinica, 2018, 49(9): 1899-1907. (in Chinese)
[49]
字向东, 罗斌, 夏威, 等. 基于RNA-Seq技术的牦牛体外受精胚胎发育转录组分析[J]. 中国农业科学, 2018, 51(8): 1577-1589.
ZI X D, LUO B, XIA W, et al. Transcriptomic analysis of IVF embryonic development in the yak (Bos grunniens) via RNA-Seq[J]. Scientia Agricultura Sinica, 2018, 51(8): 1577-1589. (in Chinese)

(编辑   郭云雁)