水通道蛋白(AQPs)是一类具有高选择性的跨膜通道蛋白,其主要功能是通过多种途径实现细胞运输水和促进一些低分子量的溶质通过[1-2]。AQPs由13个成员组成(AQP0~AQP12)[3],广泛分布于各种组织和器官中,其中水通道蛋白7(Aquaporins 7,AQP7)是一种水/甘油转运蛋白[4],允许水和小的中性电荷溶质,如甘油、尿素、嘌呤以及亚砷酸盐通过[5-6],在机体糖代谢和脂代谢过程中起着重要的作用,影响脂质和葡萄糖的稳态。AQP7已成为脂肪细胞代谢的潜在调节因子,在脂肪组织中受到胰岛素的反向调节[7-8],胰岛素通过PI3K途径影响AQP7基因启动子区的胰岛素反应元件,与成人肥胖有着密切的联系[9-10]。除此之外,其还参与了脑缺血、充血性心力衰竭、高血压和血管生成等相关疾病的心血管功能和病理过程[11]。牦牛是青藏高原上的关键物种,生长在寒冷低氧等极端环境,对青藏高原生态系统的协调和稳定发展起着重要的作用[12-13],所以牦牛的健康养殖和生产性能的提高极其关键。目前,AQP7基因主要在小鼠[14]、藏山羊[15]和人类[16]中进行研究,在高原牦牛上的研究还未见报道,其作为参与代谢的心血管功能的关键基因,对探讨牦牛在高原低氧环境下的适应性和生长有着重大的意义。本研究对九龙牦牛AQP7基因进行克隆及表达分析,探讨其在牦牛组织器官中的表达模式,旨在为后续探讨该基因的功能提供参考资料。
本试验所用牦牛来自四川九龙县,随机选择3头健康成年雄性九龙牦牛作为研究对象,采集心脏、肝脏、脾脏、肺脏、肾脏、小肠、肌肉及瘤胃8个组织,一部分组织装入冻存管中储存于液氮,另取一部分用装有4%多聚甲醛固定液的离心管保存。
AQP7兔多克隆抗体购自北京博奥森生物技术有限公司;RNAiso Plus试剂盒、PrimeScriptTM RT Reagent Kit试剂盒、pMDTM19-T Vector Cloning Kit试剂盒和DH5α感受态细胞均购于宝日医生物技术(北京)有限公司。AQP7基因扩增及荧光定量引物由上海生物工程有限公司合成(表1)。
表1 引物信息
Tab.1 Primer information
基因Gene引物序列(5′-3′)Primer sequence产物长度/bp Product length退火温度/℃Tm应用ApplicationAQP7F:ATGACTCAGTCAGACAGGAAGAGGC99348.0PCRR:TTAAGGGTTGTACAAGGCTGAGGGTAQP7-QF:GCCATCATCTACTTGTTCTTCATT17949.0q-PCRR:AAGGCTGAGGGTCTATCTTCGGGGTβ-actinF:CT TCGAGCAGGAGATGGC16754.9q-PCRR:CCGTGT TGGCGTAGAGGT
称取0.1 g牦牛心脏、肝脏、脾脏、肺脏、肾脏、小肠、肌肉和瘤胃,将组织剪碎后加入液氮至研钵中进行研磨,按照RNAiso Plus试剂盒说明书对8种组织RNA进行提取,1%琼脂糖电泳和Nano Drop 2000分光光度计对其进行检测。获得质量完好的RNA按照PrimeScriptTM RT Reagent Kit说明书进行cDNA的制备,获得的cDNA于-20 ℃保存备用。
牦牛心脏cDNA为模板扩增AQP7基因,对应的PCR扩增体系为50 μL:上游和下游引物各2 μL(0.5 μmol/L),1.1×T3 Super PCR Mix酶44 μL(1×),cDNA模板2 μL(100 ng/μL)。反应程序:98 ℃预变性2 min;98 ℃变性10 s,54 ℃退火10 s,72 ℃延伸15 s,共35个循环;72 ℃再延伸2 min,4 ℃保存。获得的PCR产物用1%琼脂糖电泳检测,按照DNA回收试剂盒说明书对目的片段进行胶回收,胶回收产物按照pMDTM19-T Vector Cloning Kit说明书进行连接,转化到DH5α中,将转化菌液涂布于含有氨苄(Amp+)的LB固体培养基上,37 ℃培养过夜,挑取单菌落,PCR进行鉴定,阳性菌液送生工生物工程(上海)股份有限公司进行测序。
测序得到的序列与目的基因序列进行比对,确定其为目的基因后对其基本性质进行分析,首先使用ORF Finder寻找序列的开放阅读框,同时得到相应的氨基酸序列。用在线软件ExPASy、ProtScal、TMHMM、SignalP和SMART对蛋白质理化性质、疏水性、跨膜结构、信号肽、结构功能域和蛋白互作进行分析。利用SOPMA软件对氨基酸二级结构进行预测,用MEGA 6.0软件进行系统进化树的构建。
利用RT-qPCR技术检测牦牛AQP7基因在心脏、肝脏、脾脏、肺脏、肾脏、小肠、肌肉和瘤胃组织中的表达差异,以β-actin作为内参矫正基因的相对表达水平。10 μL反应体系:上下游引物均为0.8 μL(0.5 μmol/L),cDNA模板1 μL(50 ng/μL),TB GreenTMPremix Ex TaqTMⅡ 5.2 μL(2×),ddH2O 2.2 μL。反应程序:95 ℃预变性3 min;95 ℃变性10 s,59.4 ℃退火10 s,共40个循环。3个生物学重复,每个样本重复检测3次,用2-ΔΔCt法进行分析,用SPSS 26.0软件分析其显著性。
用多聚甲醛固定的8种牦牛组织经脱水、包埋和切片后,进行染色,3%甲醇双氧水中室温放置10 min,用PBS将切片进行清洗,浸入pH值6.0、0.01 mol/L的柠檬酸盐缓冲液中,将其进行加热至沸腾,冷却后用PBS洗涤2次,随后向其滴加山羊血清封闭液后室温放置20 min,滴加一抗(1∶100),滴加生物素化二抗,将其放置于37 ℃中孵育30 min,用PBS清洗后室温下进行DAB显色,最后利用苏木素进行轻度复染,脱水透明后用中性树胶封片。通过显微镜观察染色结果,采集3张400倍显微图像,用Image-Pro Plus 6.0图像分析系统对结果图片进行分析,用SPSS 26.0软件分析其显著性。
心脏cDNA为模板扩增后获得与预期目的片段大小一致的基因片段(图1),经测序获得AQP7基因片段长度为993 bp,ORF为993 bp,序列共编码330个氨基酸,将获得的序列提交到NCBI,获得登录号为MW045313。
图1 AQP7基因PCR扩增
Fig.1 PCR amplification of AQP7
获得的牦牛AQP7基因序列比对后,下载相应的同属及其他动物的序列进行同源性分析,发现与黄牛的同源性最高为99.7%,其次为水牛(98.7%)、山羊(95.0%)和瘤牛(90.6%),与小鼠的同源性最低,为67.8%。系统进化树的构建发现该基因与黄牛和瘤牛聚为一类,亲缘关系最近,与水牛的物种亲缘关系次之(图2)。九龙牦牛AQP7蛋白理化性质分析结果显示,其分子式为C1651H2548N424O445S16,分子质量为35.97 ku,等电点9.32。编码的氨基酸有17个带负电残基(Asp+Glu),24个带正电残基(Arg+Lys),脂肪系数为93.94,不稳定系数为33.43,亲水性平均值为0.318,推测该蛋白为稳定疏水蛋白。功能预测发现,牦牛AQP7蛋白有典型的疏水区,此外有6个典型的跨膜螺旋区,含有1个MIP结构功能域。蛋白质互作分析发现其与同家族的AQP8、AQP11和AQP12B等存在相互作用。AQP7有α螺旋(27.27%)、β折叠(2.42%)、延伸片段(24.55%)和无规卷曲(43.33%)4种二级结构,其中无规卷曲占比最高。
图2 九龙牦牛AQP7基因进化树
Fig.2 Phylogenetic tree of AQP7 gene of Jiulong yak
牦牛AQP7基因在不同组织中的表达结果显示(图3),AQP7基因在心脏的表达量最高,肌肉和肾脏次之,表达量相对较低的组织为肝脏、脾脏、肺脏、小肠和瘤胃。AQP7基因在心脏和肌肉中的表达量极显著高于肾脏、肝脏和脾脏等其他组织(P<0.01),且在心脏中的表达量极显著高于肌肉组织(P<0.01)。
不同大写字母表示差异极显著(P<0.01);
相同字母表示差异不显著(P>0.05)。
The different capital letters show extremely significant different(P<0.01);
The same letter show no significant different between stages(P>0. 05).
图3 AQP7基因在九龙牦牛不同组织中的表达
Fig.3 Expression of AQP7 in different tissues of Jiulong yak
通过免疫组织化学方法检测,发现AQP7免疫阳性细胞在九龙牦牛心脏、肝脏、脾脏、肺脏、肾脏、小肠、肌肉和瘤胃组织中均检测到,阳性颗粒主要分布在细胞质中。AQP7蛋白表达于多种组织器官,在心脏和肌肉的肌细胞、肾脏的近曲小管内皮细胞和集合管内皮细胞中表达较为丰富。除此之外,AQP7蛋白在毛细血管内皮细胞、肝细胞、淋巴细胞、小肠固有层细胞核中也都有表达(图4)。平均光密度值分析可知(图5),牦牛AQP7蛋白在心脏中的表达量最高,肾脏和肌肉次之,在肝脏、脾脏、肺脏、小肠和瘤胃中的表达量相对较低,且AQP7蛋白在心脏、肾脏和肌肉中的表达显著高于肝脏、脾脏、肺脏、小肠和瘤胃(P<0.05)。
A1~H1.牦牛AQP7蛋白在心脏、肝脏、脾脏、肺脏、肾脏、小肠、肌肉、瘤胃中的免疫组织染色;A2~H2.阴性对照。MC.肌细胞;KC.枯否细胞;HCs.肝细胞;L.淋巴细胞;PC.上皮细胞;G.肾小球;CT.集合管;PT.近曲小管;IG.肠腺; IP.固有层; SG.颗粒层。
A1-H1.The immunohistochemical staining of yak AQP7 protein in heart,liver,spleen,lung,kidney,intestine,muscle,rumen;A2-H2.The negative control.MC.Muscle cells;KC.Kupffer cell;HCs.Hepatocyte;L.Lymphocyte;PC.Epithelial cells;G.Glomerular;CT.Collecting tubes;PT.Proximal convoluted tubule;IG.Intestinal gland;IP.Lamina propria;SG.Stratum granulosum.
图4 AQP7在九龙牦牛组织的分布(400×)
Fig.4 Distribution of AQP7 in tissues of Jiulong yak(400 ×)
不同小写字母表示差异显著(P<0.05);
相同字母表示差异不显著(P>0.05)。
The different small letters show significant difference(P<0.05);
The same letter show no significant difference between stages(P>0. 05).
图5 九龙牦牛AQP7阳性反应平均光密度值分析
Fig.5 Average integral optical density analysis of AQP7 expression of Jiulong yak
目前,关于水通道蛋白转运水和其他小分子物质的功能研究已相当完善,但对于其通过转运功能参与其他维持机体相对稳态的机制还不清楚。近年来,水通道蛋白在皮肤病[17]、肾脏疾病[18]和癌症[19]等疾病中的研究越来越多,关于其低氧适应性和生长机制的研究相对较少,牦牛作为高原上耐寒和耐低氧的物种,了解其天然的适应性机制有利于解决目前适应高原环境的问题以及提升牦牛生产性能。前人的研究表明,AQP4和AQP9基因与牦牛的低氧适应性相关[20-21],对AQP7基因的研究有利于加强水通道蛋白家族在低氧适应性过程中的作用机制研究,并且参与糖代谢和脂代谢[22-23],在一定程度上有利于进一步加强探讨其低氧适应性机制。因此,本研究克隆获得九龙牦牛编码区为993 bp 的AQP7基因,编码330个氨基酸,与张亚楠等[15]克隆获得的藏山羊的AQP7基因片段长度一致,同源性也相对较高,其具有的跨膜结构和MIP结构域与藏山羊一致,但蛋白分子式和氨基酸的序列有一定的差异。牦牛和藏山羊均生活在高海拔地区具有较强的抗逆性,所以在一定程度上有相似性。除此之外,蛋白基本理化性质分析和结构也为进一步探讨牦牛AQP7蛋白在机体当中可能的功能作用提供参考资料。
前人研究表明,人AQP7基因在肾脏[24]、脂肪组织[25]、睾丸[26]和心脏[27]中均有表达,本研究通过RT-qPCR方法发现,AQP7基因在九龙牦牛心脏、肝脏和脾脏等8种组织中均有表达,且在心脏、肾脏和肌肉中的表达量最高,与藏山羊的组织表达相对一致[15]。AQPs在心血管系统中举足轻重,AQPs参与了一些相关疾病的心血管功能和病理过程[28],九龙牦牛AQP7在心脏中高表达说明了其在心脏的功能作用中起着重要的作用,这可能是由于牦牛心脏对高原环境有一定的适应性所导致的,可以通过对心脏血液循环系统的调节来促进机体在高原环境的适应性;在肾脏中的高表达可能是由于肾脏作为生物体的泌尿系统,参与了水的过滤和重吸收等过程,水的运输和其他小分子物质的运输均依赖于水通道蛋白的功能作用;肌肉当中AQP7的高表达,可能是因为肌肉参与了机体的能量代谢过程,且肌肉组织中含有脂肪成分。已有报道表明,人AQP7参与机体的脂代谢过程,调节机体的能量稳态[29]。
组织定位研究发现,AQP7蛋白在九龙牦牛心脏、肝脏和脾脏等8个组织当中均有阳性颗粒,和RT-qPCR结果基本一致。AQP7定位于小鼠脂肪组织、心脏和横纹肌的毛细血管[23],这与本研究结果相似。AQP7基因敲除小鼠心脏的甘油和ATP含量降低,说明AQP7对阐明心脏生理和病理的过程具有重要意义[28]。此外,AQPs在肾脏特异性细胞类型和液体稳态中的生理和病理功能已被深入研究[30],并可能成为一种新型的肾脏疾病的理想生物标志物。综上所述,可以推断AQP7蛋白在牦牛心脏、肾脏和肌肉中起着至关重要的作用。牦牛在生理学和组织学上都具有独特的性状,以适应青藏高原低氧高海拔的恶劣生活环境,具有与黄牛相比更敏感的低氧感知能力和更高的代谢能力[31]。心脏和肌肉均作为能量代谢的重要器官,在对环境的适应性方面起着重要的作用,AQP7通过参与牦牛心脏和肌肉的能量代谢过程来影响牦牛对极端环境的适应性,且AQP7在脂代谢的过程中有着不可替代的作用,本研究结果对于了解AQP7在牦牛生长过程中对脂代谢的调控机制提供了重要的参考数据,并且为探讨牦牛通过AQP7参与能量代谢来影响高原低氧环境的适应性机制奠定了基础。
[1] Prieto-Martínez N,Vilagran I,Morat R,Rivera del
lamo M M,Rodríguez-Gil J E,Bonet S,Yeste M. Relationship of aquaporins 3(AQP3),7(AQP7),and 11(AQP11)with boar sperm resilience to withstand freeze-thawing procedures[J].Andrology,2017,5(6):1153-1164. doi:10.1111/andr.12410.
[2] Skowronski M T,Skowronska A,Rojek A,Oklinski M K,Nielsen S. Prolonged starvation causes up-regulation of AQP1 in adipose tissue capillaries of AQP7 knock-out mice[J].International Journal of Molecular Sciences,2016,17(7):1101. doi:10.3390/ijms17071101.
[3] Iena F M,Lebeck J. Implications of aquaglyceroporin 7 in energy metabolism[J].International Journal of Molecular Sciences,2018,19(1):154. doi:10.3390/ijms19010154.
[4] Fu X H,Zhu J,Zhang L,Shu J. Long non-coding RNA NEAT1 promotes steatosis via enhancement of estrogen receptor alpha-mediated AQP7 expression in HepG2 cells[J].Artificial Cells,Nanomedicine,and Biotechnology,2019,47(1):1782-1787.doi:10.1080/21691401.2019.1604536.
[5] Liu Z J,Shen J,Carbrey J M,Mukhopadhyay R,Agre P,Rosen B P. Arsenite transport by mammalian aquaglyceroporins AQP7 and AQP9[J].Proceedings of the National Academy of Sciences of the United States of America,2002,99(9):6053-6058. doi:10.1073/pnas.092131899.
[6] Rojek A,Praetorius J,Frøkiaer J,Nielsen S,Fenton R A. A current view of the mammalian aquaglyceroporins[J].Annual Review of Physiology,2008,70:301-327. doi:10.1146/annurev.physiol.70.113006.100452.
[7] Lebeck J,Østergård T,Rojek A,Füchtbauer E M,Lund S,Nielsen S,Praetorius J. Gender-specific effect of physical training on AQP7 protein expression in human adipose tissue[J].Acta Diabetologica,2012,49(S1):S215-S226. doi:10.1007/s00592-012-0430-1.
[8] Kishida K,Kuriyama H,Funahashi T,Shimomura I,Kihara S,Ouchi N,Nishida M,Nishizawa H,Matsuda M,Takahashi M,Hotta K,Nakamura T,Yamashita S,Tochino Y,Matsuzawa Y. Aquaporin adipose,a putative glycerol channel in adipocytes[J].Journal of Biological Chemistry,2000,275(27):20896-20902. doi:10.1074/jbc.M001119200.
[9] Zheng D,Peng M Y,Zhu X N,Wang J. Effects of qutanhuoxue decoction on AQP7 and AQP9 expression in nonalcoholic fatty liver model rats[J].Evidence Based Complementary and Alternative Medicine,2019,2019:5709626. doi:10.1155/2019/5709626.
[10] Fu X H,Xing L L,Xu W H,Shu J. Treatment with estrogen protects against ovariectomy-induced hepatic steatosis by increasing AQP7 expression[J].Molecular Medicine Reports,2016,14(1):425-431.doi:10.3892/mmr.2016.5236.
[11] Gomes D,Agasse A,Thiébaud P,Delrot S,Gers H,Chaumont F. Aquaporins are multifunctional water and solute transporters highly divergent in living organisms[J].Biochimicaet Biophysica Acta,2009,1788(6):1213-1228. doi:10.1016/j.bbamem.2009.03.009.
[12] 毛进彬,毛旭东,代舜尧,涂永强,邓启红,沈常欢.亚丁牦牛生长发育性能研究[J].中国草食动物科学,2020,40(5):74-76. doi:10.3969/j.issn.2095-3887.2020 .05.017.
Mao J B,Mao X D,Dai S Y,Tu Y Q,Deng Q H,Shen C H. Research on the growth and development performance of Yading yak[J]. China Herbivore Science,2020,40(5):74-76.
[13] 贾功雪,丁路明,徐尚荣,方有贵,付弘赟,杨其恩. 青藏高原牦牛遗传资源保护和利用:问题与展望[J].生态学报,2020,40(18):6314-6323. doi:10.5846/stxb201912232763.
Jia G X,Ding L M,Xu S R,Fang Y G,Fu H Y,Yang Q E. Conservation and utilization of yak genetic resources in Qinghai-Tibet plateau:Problems and perspectives[J].Acta Ecologica Sinica,2020,40(18):6314-6323.
[14] Lebeck J,Søndergaard E,Nielsen S. Increased AQP7 abundance in skeletal muscle from obese men with type 2 diabetes[J].American Journal of Physiology Endocrinology and Metabolism,2018,315(3): 367-373. doi:10.1152/ajpendo.00468.2017.
[15] 张亚楠,王永,朱江江,许晴,林亚秋.藏山羊AQP7基因克隆和不同组织器官差异表达分析[J].西南农业学报,2019,32(8):1956-1960. doi:10.16213/j. cnki. scjas. 2019. 8.042.
Zhang Y N,Wang Y,Zhu J J,Xu Q,Lin Y Q. Analysis of cloning and expression of AQP7 gene in Tibetan goat[J].Southwest China Journal of Agricultural Sciences,2019,32(8):1956-1960.
[16] 郭榕. AQP7介导的人参皂苷Rb1改善肥胖的作用及其机制研究[D].福州:福建医科大学,2018.
Guo R. AQP7 involved in ginsenoside-Rb1-mediated anti-obestiy through the activation of PPAR-γ pathway [D].Fuzhou:Fujian Medical University,2018.
[17] Hara-Chikuma M,Satooka H,Watanabe S,Honda T,Miyachi Y,Watanabe T,VerkmanA S. Aquaporin-3-mediated hydrogen peroxide transport is required for NF-κB signalling in keratinocytes and development of psoriasis[J].Nature Communications,2015,6:7454. doi:10.1038/ncomms8454.
[18] Chen Y H,Jin S,Teng X,Hu Z J,Zhang Z H,Qiu X,Tian D Y,Wu Y M. Hydrogen sulfide attenuates LPS-induced acute kidney injury by inhibiting inflammation and oxidative stress[J].Oxidative Medicine and Cellular Longevity,2018,2018:6717212. doi:10.1155/2018/6717212.
[19] Vieceli Dalla Sega F,Zambonin L,Fiorentini D,Rizzo B,Caliceti C,Landi L,Hrelia S,Prata C. Specific aquaporins facilitate Nox-produced hydrogen peroxide transport through plasma membrane in leukaemia cells[J].Biochimicaet Biophysica Acta,2014,1843(4):806-814. doi:10.1016/j.bbamcr.2014.01.011.
[20] 宋国强. 牦牛AQP4基因的克隆及其对高原脑水肿抗性的表达机理[D].兰州:兰州大学,2011.
Song G Q. AQP4 gene cloning and its expression mechanism of resistance to high altitude cerebral edema in the yak(Bos grunniens)[D].Lanzhou:Lanzhou University,2011.
[21] 汪迁昊. 牦牛脑AQP9的功能特征[D].兰州:兰州大学,2017.
Wang Q H. Functional characteristics of AQP9 in the yak brain [D].Lanzhou:Lanzhou University,2017.
[22] Mourelatou R,Kostopoulou E,Rojas-Gil A P,Kehagias I,Linos D,Kalfarentzos F E,Spiliotis B E. Decreased adipocyte glucose transporter 4(GLUT4)and aquaglyceroporin-7(AQP7)in adults with morbid obesity:Possible early markers of metabolic dysfunction[J].Hormones,2019,18(3):297-306. doi:10.1007/s42000-019-00130-8.
[23] Skowronski M T,Lebeck J,Rojek A,Praetorius J,Füchtbauer E M,Frøkiaer J,Nielsen S. AQP7 is localized in capillaries of adipose tissue,cardiac and striated muscle:Implications in glycerol metabolism[J].American Journal of Physiology Renal-Physiology,2007,292(3):956-965. doi:10.1152/ajprenal.00314.2006.
[24] Li Y J,Wang W L,Jiang T,Yang B X. Aquaporins in urinary system[J].Advances in Experimental Medicine and Biology,2017,969:131-148. doi:10.1007/978-94-024-1057-0_9.
[25] Guo M,Chen F,Lin T F,Peng Y Q,Li W P,Zhu X X,Lin L,Chen Y S. Apelin-13 decreases lipid storage in hypertrophic adipocytes in vitro through the upregulation of AQP7 expression by the PI3K signaling pathway[J].Medical Science Monitor,2014,20:1345-1352. doi:10.12659/MSM.890124.
[26] Lebeck J. Metabolic impact of the glycerol channels AQP7 and AQP9 in adipose tissue and liver[J].Journal of Molecular Endocrinology,2014,52(2):165-178. doi:10.1530/JME-13-0268.
[27] Karaca A,Palaby
k O,Ta ʂtekin E,Turan F N,Vardar S A. High fructose diet suppresses exercise-induced increase in AQP7 expression in the in vivo rat heart[J].Anatolian Journal of Cardiology,2016,16(12):916-922. doi:10.14744/anatoljcardiol.2016.6958.
[28] Tie L,Wang D,Shi Y D,Li X J. Aquaporins in cardiovascular system[J].Advances in Experimental Medicine and Biology,2017,969:105-113. doi:10.1007/978-94-024-1057-0_6.
[29] Oikonomou E,Kostopoulou E,Rojas-Gil A P,Georgiou G,Spiliotis B E. The metabolic implications of aquaporin 7(AQP7)promoter variants in lean children and children with obesity[J].Hormones,2020,19(2):187-195. doi:10.1007/s42000-020-00184-z.
[30] He J Z,Yang B X. Aquaporins in renal diseases[J].International Journal of Molecular Sciences,2019,20(2):366. doi:10.3390/ijms20020366.
[31] 贺大芳.牦牛和黄牛心脏、肺脏组织microRNA转录组的鉴定与差异表达分析[D].雅安:四川农业大学,2017.
He D F.Identification and differential expression analysis of microRNAs in heart and lung tissues between yak and cattle[D].Ya'an:Sichuan Agricultural University,2017.