绵羊不同妊娠时期卵巢转录组差异表达比较分析

doi:10.12160/j.issn.1672-5190.2021.04.007

摘要/Abstract

摘要： [目的]探究随着绵羊体内胎儿的发育,母体卵巢转录组表达水平的变化。[方法]试验绵羊为苏格兰黑面羊和特克赛尔羊F₁代杂交品种,3个妊娠时间分别为妊娠23 d、妊娠35 d和妊娠100 d。从NCBI的高通量测序数据SRA存储库下载3个妊娠时期绵羊卵巢转录组的18个样本数据。使用Z-Score方法对转录组数据进行标准化处理。将妊娠23 d、妊娠35 d的绵羊卵巢组织进行比对,妊娠35 d、妊娠100 d的绵羊卵巢组织进行比对,使用R语言的limma包进行基因差异表达分析,筛选DEGs,以P-Value<0.05和|logFC|≥2作为筛选DEGs条件。利用DAVID在线软件对DEGs进行注释及功能富集分析,并以P-Value<0.05作为信号通路显著富集标志。[结果]在妊娠23 d与妊娠35 d组绵羊卵巢中筛选出54个DEGs,妊娠35 d与妊娠100 d组绵羊卵巢中筛选出68个DEGs,两组均包含HSP90AA1、HSP90AB1、MMP2和HSP90B1基因。妊娠23 d与35 d组绵羊卵巢的54个DEGs显著富集到包括内质网蛋白质加工信号通路、核糖体信号通路、雌激素信号通路、癌症蛋白聚糖信号通路、扩张型心肌病信号通路、肥厚型心肌病信号通路6条通路;妊娠35 d、妊娠100 d组绵羊卵巢的68个DEGs显著富集到包括癌症相关信号通路、黏着斑信号通路、抗原加工与提呈信号通路、细胞外基质受体相互作用信号通路、细菌侵袭上皮细胞信号通路、内质网蛋白质加工信号通路、核糖体信号通路、雌激素信号通路、PI3K-Akt信号通路、癌症蛋白聚糖信号通路和前列腺癌信号通路11条通路。[结论]在妊娠期间雌激素信号通路上的HSP90AA1、HSP90AB1、MMP2和HSP90B1基因表达有很大变化,这些基因可能对绵羊妊娠维持有着重要的作用。

关键词: 绵羊, 转录组, 妊娠, 卵巢

Abstract: [Objective] To transcriptomically identify the differentially expressed genes （DEGs） in maternal ovary with the development of fetus in sheep. [Method] F₁ hybrids of Scottish black-faced sheep and Texel sheep were used as experimental animals, and the transcriptomic analysis profiles were collected on 23^th day, 35^th day and 100^th day of gestation, respectively. A total of 18 transcriptomic data samples of sheep ovaries in the above three periods of gestation were downloaded from SRA in NCBI, a high-throughput sequencing data storing database. The transcriptomic data were standardized using Z-Score method. Comparative transcriptomic analysis was carried out between 23^th day and 35^th day of gestation as well as between 35^th day and 100^th day of gestation. The limma package of R language was used for analysis of differential gene expression to screen DEGs, and P-Value < 0.05 and |logFC| ≥ 2 were used as screening conditions for DEGs. DAVID online software was used to annotate DEGs and analyze their functional enrichment, and P-Value < 0.05 was used as a significant enrichment marker of signaling pathway. [Result]A total of 54 DEGs in sheep ovaries were screened between 23^th day and 35^th day of gestation, and a total of 68 DEGs were screened between 35^th day and 100^th day of gestation. The genes of HSP90AA1,HSP90AB1,MMP2 and HSP90B1 were observed in both groups of the DEGs. The 54 DEGs identified between 23^th day and 35^th day of gestation were found to be significantly enriched in six signaling pathways, including protein processing in endoplasmic reticulum, ribosome, estrogen, proteoglycans in cancer, dilated cardiomyopathy, and hypertrophic cardiomyopathy. The 68 DEGs identified between 35^th day and 100^th day of gestation were found to be significantly enriched in eleven signaling pathways, including cancer, focal adhesion, antigen processing and presentation, extracellular matrix （ECM）-receptor interaction, bacterial invasion of epithelial cells, protein processing in endoplasmic reticulum, ribosome, estrogen, PI3K-Akt, proteoglycans in cancer, and prostate cancer. [Conclusion] The gene expressions of HSP90AA1、HSP90AB1、MMP2 and HSP90B1 in estrogen signaling pathway during gestation have significant changes, indicating that these genes may play important roles in sheep gestation maintenance.

Key words: sheep, transcriptome, gestation, ovary

中图分类号:

S826

史静茹, 郭丽丽, 刘在霞, 刘永斌, 张家新, 张文广. 绵羊不同妊娠时期卵巢转录组差异表达比较分析[J]. 畜牧与饲料科学, 2021, 42(4): 37-41.

SHI Jing-ru, GUO Li-li, LIU Zai-xia, LIU Yong-bin, ZHANG Jia-xin, ZHANG Wen-guang. Comparative Transcriptomic Analysis to Identify Differentially Expressed Genes in Sheep Ovaries during Different Periods of Gestation[J]. Animal Husbandry and Feed Science, 2021, 42(4): 37-41.

参考文献

[1] 吕桂萍,文亚洲.提高绵羊繁殖力的关键措施[J].养殖与饲料,2017(6):37-38.
[2] MCNATTY K P, GALLOWAY S M, WILSON T, et al.Physiological effects of major genes affecting ovulation rate in sheep[J].Genetics Selection Evolution,2005,37(S1):25-38.
[3] NOSATI M,NANAEI H A,GHANATSAMAN Z A,et al.Whole genome sequence analysis to detect signatures of positive selection for high fecundity in sheep[J].Reproduction in Domestic Animals,2019,54(2):358-364.
[4] BODIN L,PASQUALE D E,FABRE S,et al.A novel mutation in the bone morphogenetic protein 15 gene causing defective protein secretion is associated with both increased ovulation rate and sterility in Lacaune sheep[J].Endocrinology,2007,148(1):393-400.
[5] 田志龙,刘秋月,王翔宇,等.绵羊BMP15基因多态性及其与产羔数的关联分析[J].中国畜牧杂志,2018,54(7):23-27.
[6] 田志龙. 绵羊繁殖候选基因的组织表达、SNP及其与产羔数的关联分析[D].洛阳:河南科技大学,2019.
[7] YANG H, LIN S, LEI X,et al.Identification and profiling of microRNAs from ovary of estrous Kazakh sheep induced by nutritional status in the anestrous season[J].Animal Reproduction Science,2016(175):18-26.
[8] 安春利. 孕激素与早期妊娠[J].幸福家庭,2020(11):8.
[9] 田晓,孙振高.早孕期雌激素水平预测妊娠结局的研究进展[J].医学综述,2020,26(5):872-877.
[10] 赵梦圆. 绵羊早期妊娠诊断关键生殖激素及差异蛋白筛选[D].晋中:山西农业大学,2019.
[11] 庄晨. 孕酮相关基因在绵羊早期妊娠肝脏的表达及在生产上的应用[D].邯郸:河北工程大学,2019.
[12] CLARK E L,BUSH S J, MCCULLOCH M E B,et al. A high resolution atlas of gene expression in the domestic sheep (Ovis aries)[J].PLoS Genet,2017,13(9):e1006997.
[13] RITCHIE M E,BELINDA P,DI W,et al.Limma powers differential expression analyses for RNA-sequencing and microarray studies[J].Nucleic Acids Research,2015,43(7):e47.
[14] HUANG D W,SHERMAN B T,TAN Q,et al.DAVID bioinformatics resources: expanded annotation database and novel algorithms to better extract biology from large gene lists[J].Nucleic Acids Research,2007,35(S2):169-175.
[15] 任巧玲,马强,白红杰,等. 雌激素及其对动物生殖功能的影响[J].河南农业科学, 2006(4):114-115.
[16] FILIS P,NAGRATH N,FRASER M,et al.Maternal smoking dysregulates protein expression in second trimester human fetal livers in a sex-specific manner[J].The Journal of Clinical Endocrinology and Metabolism,2015,100(6):861-870.
[17] TESFAYE D,GHANEM N,CARTER F,et al.Gene expression profile of cumulus cells derived from cumulus-oocyte complexes matured either in vivo or in vitro[J].Reproduction,Fertility and Development,2009,21(3):451-461.
[18] TORRES-SALZAR Q,MARTINEZ-LOPEZ Y,REYES-ROMERO M,et al.Differential methylation in promoter regions of the genes NR3C1 and HSP90AA1, involved in the regulation, and bioavailability of cortisol in leukocytes of women with preeclampsia[J].Frontiers in Medicine,2020,7:206.
[19] NIKISHIN D A,FILATOV M A,KISELEVA M V,et al.Selection of stable expressed reference genes in native and vitrified/thawed human ovarian tissue for analysis by qRT-PCR and Western blot[J].Journal of Assisted Reproduction and Genetics,2018,35(10):1851-1860.
[20] CHEN J J,KHALIL R A .Matrix Metalloproteinases in normal pregnancy and preeclampsia[J].Progress in Molecular Biology and Translational science, 2017,148:87-165.
[21] KIZAKI K,USHIZAWA K,TAKAHASHI T,et al.Gelatinase (MMP-2 and -9) expression profiles during gestation in the bovine endometrium[J].Reproductive Biology and Endocrinology,2008,6(1):66.
[22] RAMU D,VENKATESAN V,PAUL S F D,et al. Genetic variation in matrix metalloproteinase MMP2 and MMP9 as a risk factor for idiopathic recurrent spontaneous abortions in an Indian population[J].Journal of Assisted Reproduction and Genetics,2017,34(7):945-949.
[23] YAN Y,FANG L L,LI Y X,et al.Association of MMP2 and MMP9 gene polymorphisms with the recurrent spontaneous abortion:A meta-analysis[J].Gene,2021,767:145173.
[24] PEREZA N,OSTOJI S,VOLK M,et al.Matrix metalloproteinases 1, 2, 3 and 9 functional single-nucleotide polymorphisms in idiopathic recurrent spontaneous abortion[J].Reproductive Biomedicine Online,2012,24(5):567-575.
[25] GREMLICH S,NGUYEN D,REYMONDIN D,et al.Fetal MMP2/MMP9 polymorphisms and intrauterine growth restriction risk[J].Journal of Reproductive Immunology,2007,74(1/2):143-151.
[26] ALAM S M K,KONNO T,SOARES M J. Identification of target genes for a prolactin family paralog in mouse decidua[J].Reproduction,2015,149(6):625-632.