发情周期不同阶段绵羊卵巢lncRNA的鉴定和功能分析

doi:10.11843/j.issn.0366-6964.2022.12.009

摘要/Abstract

摘要： 旨在分析不同繁殖周期绵羊卵巢组织中长链非编码RNA (long non-coding RNA，lncRNA)的表达谱，了解lncRNA表达及其调控机理，为绵羊繁育研究提供理论依据。本研究以湖羊为研究对象，选取年龄在1.5～2.5岁的黄体期和卵泡期母羊各3只，通过高通量测序筛选出卵巢组织中的lncRNA，运用生物信息学分析对差异表达的lncRNA进行靶基因预测，通过GO和KEGG富集分析找出与绵羊繁殖相关的通路。结果显示，本研究共获得1 379个差异表达的lncRNAs，其中1 158个表达上调，221个表达下调。GO和KEGG富集分析表明，差异表达的lncRNAs及其靶基因主要参与卵泡发育、排卵周期过程、钙离子信号通路、卵母细胞减数分裂、催产素信号通路、MAPK信号通路、甲状腺激素合成通路、雌激素信号通路。关键的lncRNAs可能通过调控参与这些信号通路和生物学过程的相关基因来调控生殖。其中，LNC_011239、LNC_012847、LNC_003902、LNC_003906、LNC_003907等靶向的MAPK1、ADCY1、ADCY5、PPP3CA和CDC23可能发挥关键的调控作用。qRT-PCR验证证明，随机选取的5个差异lncRNAs定量结果与测序结果基本一致。本研究利用RNA-Seq技术筛选出黄体期和卵泡期卵巢组织中的lncRNAs，并进行差异分析，为揭示绵羊繁殖能力的分子机制提供依据。

关键词: 卵巢, lncRNA, RNA-seq, 绵羊

Abstract: This study aimed to analyze the expression profiles of long non-coding RNA (lncRNA) in sheep ovarian tissues at different reproductive cycles, understand the lncRNA expression and its regulatory mechanism, and provide a theoretical basis for sheep breeding research. In this study, Hu sheep were selected as the research object, and 3 ewes aged from 1.5 to 2.5 years old in luteal phase and follicular phase were selected, respectively. lncRNA in ovarian tissues were screened by high-throughput sequencing. Bioinformatics analysis was used to predict the target genes of differentially expressed lncRNAs. GO and KEGG enrichment analysis were used to identify pathways related to sheep reproduction. The results showed that a total of 1 379 differentially expressed lncRNAs were obtained in this study, of which 1 158 were up-regulated and 221 down-regulated. GO and KEGG enrichment analysis showed that the differentially expressed lncRNAs and their target genes mainly involved in ovarian follicle development, ovulation cycle process, calcium signaling pathway, oocyte meiosis, oxytocin signaling pathway, MAPK signaling pathway, thyroid hormone synthesis signaling pathway, and estrogen signaling pathway. Key lncRNAs might regulate reproduction by regulating genes involved in these signaling pathways and biological processes. Specifically, MAPK1, ADCY1, ADCY5, PPP3CA and CDC23, which are targeted by LNC_011239, LNC_012847, LNC_003902, LNC_003906, LNC_003907 and others, might play key regulatory roles. The results of qRT-PCR showed that the expression levels of 5 randomly selected lncRNAs were consistent with the sequencing results. In this study, RNA-Seq technology was used to screen lncRNAs in luteal and follicular ovarian tissues, and the differences were analyzed to provide the basis for revealing the molecular mechanism of sheep reproductive capacity.

Key words: ovary, lncRNA, RNA-seq, sheep

中图分类号:

S826.3

王梦瑶, 翟振翰, 赵路, 王赛桥, 王玉琴. 发情周期不同阶段绵羊卵巢lncRNA的鉴定和功能分析[J]. 畜牧兽医学报, 2022, 53(12): 4221-4231.

WANG Mengyao, ZHAI Zhenhan, ZHAO Lu, WANG Saiqiao, WANG Yuqin. Identification and Functional Analysis of lncRNA in Ovaries of Sheep at Different Stages of Estrus Cycle[J]. Acta Veterinaria et Zootechnica Sinica, 2022, 53(12): 4221-4231.

参考文献 41

[1]	甘春艳,李军,金海.2021年我国肉羊产业发展概况、未来发展趋势及建议[J].中国畜牧杂志,2022,58(3):258-263.GAN C Y,LI J,JIN H.Development situation,future development trend and suggestions of meat sheep industry in China in 2021[J].Chinese Journal of Animal Science,2022,58(3):258-263.(in Chinese)
[2]	周书杰,张应利,张传辉.全国肉羊养殖成本与收益分析[J].合作经济与科技,2022(3):100-101.ZHOU S J,ZHANG Y L,ZHANG C H.Cost and benefit analysis of mutton sheep breeding in China[J].Co-Operative Economy & Science,2022(3):100-101.(in Chinese)
[3]	LAURENT G S,WAHLESTEDT C,KAPRANOV P.The Landscape of long noncoding RNA classification[J].Trends Genet,2015,31(5):239-251.
[4]	TUCK A C,TOLLERVEY D.A transcriptome-wide atlas of RNP composition reveals diverse classes of mRNAs and lncRNAs[J]. Cell,2013,154(5):996-1009.
[5]	DERRIEN T,JOHNSON R,BUSSOTTI G,et al.The GENCODE v7 catalog of human long noncoding RNAs:analysis of their gene structure,evolution,and expression[J].Genome Res,2012,22(9):1775-1789.
[6]	KNOWLING S,MORRIS K V.Non-coding RNA and antisense RNA.Nature’s trash or treasure?[J].Biochimie, 2011,93(11): 1922-1927.
[7]	HOMBACH S,KRETZ M.Non-coding RNAs:classification,biology and functioning[J].Adv Exp Med Biol,2016,937:3-17.
[8]	STATELLO L,GUO C J,CHEN L L,et al.Gene regulation by long non-coding RNAs and its biological functions[J].Nat Rev Mol Cell Biol, 2021,22(2):96-118.
[9]	MATTICK J S.Long noncoding RNAs in cell and developmental biology[J].Semin Cell Dev Biol,2011,22(4):327.
[10]	CHEN Y G,SATPATHY A T,CHANG H Y.Gene regulation in the immune system by long noncoding RNAs[J].Nat Immunol, 2017, 18(9):962-972.
[11]	李耀坤,许香萍,陆婷婷,等.影响山羊产羔数lncRNA的筛选及功能分析[J].畜牧兽医学报,2020,51(8):1853-1865.LI Y K,XU X P,LU T T,et al.Screening and functional analysis of lncRNA affecting litter size in goats[J].Acta Veterinaria et Zootechnica Sinica,2020,51(8):1853-1865.(in Chinese)
[12]	FENG X,LI F Z,WANG F,et al.Genome-wide differential expression profiling of mRNAs and lncRNAs associated with prolificacy in Hu sheep[J].Biosci Rep,2018,38(2):BSR20171350.
[13]	SU T,YU H L,LUO G,et al.The interaction of lncRNA XLOC-2222497,AKR1C1,and progesterone in porcine endometrium and pregnancy[J].Int J Mol Sci,2020,21(9):3232.
[14]	CHEN S,GUO X F,HE X Y,et al.Transcriptome analysis reveals differentially expressed genes and long non-coding RNAs associated with fecundity in sheep hypothalamus with different FecB genotypes[J].Front Cell Dev Biol,2021,9:633747.
[15]	PERTEA M,KIM D,PERTEA G M,et al.Transcript-level expression analysis of RNA-seq experiments with HISAT,stringtie and ballgown[J].Nat Protoc,2016,11(9):1650-1667.
[16]	PERTEA M,PERTEA G M,ANTONESCU C M,et al.StringTie enables improved reconstruction of a transcriptome from RNA-seq reads[J].Nat Biotechnol, 2015,33(3):290-295.
[17]	YOUNG M D,WAKEFIELD M J,SMYTH G K,et al.Gene ontology analysis for RNA-seq:accounting for selection bias[J]. Genome Biol, 2010,11(2):R14.
[18]	KANEHISA M,ARAKI M,GOTO S,et al.KEGG for linking genomes to life and the environment[J].Nucleic Acids Res,2008, 36:D480-D484.
[19]	MAO X Z,CAI T,OLYARCHUK J G,et al.Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary[J].Bioinformatics,2005,21(19):3787-3793.
[20]	马骏杰,汪亮,刘娣,等.冷应激后民猪背脂中差异lncRNA的筛选与分析[J].畜牧兽医学报,2021,52(4):932-942.MA J J,WANG L,LIU D,et al.Screening and analysis of differentially expressed lncRNA in Min pig back fat after cold stress[J].Acta Veterinaria et Zootechnica Sinica,2021,52(4):932-942.(in Chinese)
[21]	王瑞雪,刘佳森,李蕴华,等.苏尼特羊不同生长时期肌肉组织lncRNA的差异表达分析[J].中国农业大学学报,2021, 26(1): 51-61.WANG R X,LIU J S,LI Y H,et al.Differential expression analysis of lncRNA in muscle tissue of Sunite sheep at different growth stages[J].Journal of China Agricultural University,2021,26(1):51-61.(in Chinese)
[22]	NAKAGAWA S,SHIMADA M,YANAKA K,et al.The lncRNA Neat1 is required for corpus luteum formation and the establishment of pregnancy in a subpopulation of mice[J].Development,2014,141(23):4618-4627.
[23]	LI W B,NOTANI D,MA Q,et al.Functional roles of enhancer RNAs for oestrogen-dependent transcriptional activation[J]. Nature,2013,498(7455):516-520.
[24]	LI J,CAO Y X,XU X F,et al.Increased new lncRNA-mRNA gene pair levels in human cumulus cells correlate with oocyte maturation and embryo development[J].Reprod Sci,2015,22(8):1008-1014.
[25]	GAO W L,LIU M,YANG Y,et al.The imprinted H19 gene regulates human placental trophoblast cell proliferation via encoding miR-675 that targets Nodal Modulator 1 (NOMO1)[J].RNA Biol,2012,9(7):1002-1010.
[26]	JAFFE L A,EGBERT J R.Regulation of mammalian oocyte meiosis by intercellular communication within the ovarian follicle[J].Annu Rev Physiol,2017,79:237-260.
[27]	RUSNAK F,MERTZ P.Calcineurin:form and function[J].Physiol Rev,2000,80(4):1483-1521.
[28]	ISLAM R,LIU X X,GEBRESELASSIE G,et al.Genome-wide association analysis reveals the genetic locus for high reproduction trait in Chinese Arbas Cashmere goat[J].Genes Genomics,2020,42(8):893-899.
[29]	BAI Y Y,LI J,ZHU H J,et al.Deletion mutation within the goat PPP3CA gene identified by GWAS significantly affects litter size[J].Reprod Fertil Dev,2021,33(7):476-483.
[30]	ZHOU Q,LI J,YUE W,et al.Cell division cycle 23 is required for mouse oocyte meiotic maturation[J].FASEB J,2020,34(7): 8990-9002.
[31]	LIU A J,LIU M H,LI Y X,et al.Differential expression and prediction of function of lncRNAs in the ovaries of low and high fecundity Hanper sheep[J].Reprod Domest Anim,2021,56(4):604-620.
[32]	KANG H M,JEONG C B,LEE Y H,et al.Cross-reactivities of mammalian MAPKs antibodies in rotifer and copepod:Application in mechanistic studies in aquatic ecotoxicology[J].Mar Pollut Bull,2017,124(2):614-623.
[33]	刘永刚,李文哲,苏建民.哺乳动物卵母细胞成熟相关信号的调控机制研究进展[J].动物医学进展,2019,40(9):93-97.LIU Y G,LI W Z,SU J M.Progress on regulation mechanisms of mammalian oocyte maturation signals[J].Progress in Veterinary Medicine,2019,40(9):93-97.(in Chinese)
[34]	DEVI Y S,SEIBOLD A M,SHEHU A,et al.Inhibition of MAPK by prolactin signaling through the short form of its receptor in the ovary and decidua:involvement of a novel phosphatase[J].J Biol Chem,2011,286(9):7609-7618.
[35]	SALVADOR L M,MAIZELS E,HALES D B,et al.Acute signaling by the LH receptor is independent of protein kinase C activation[J].Endocrinology,2002,143(8):2986-2994.
[36]	SELA-ABRAMOVICH S,CHOREV E,GALIANI D,et al.Mitogen-activated protein kinase mediates luteinizing hormone-induced breakdown of communication and oocyte maturation in rat ovarian follicles[J].Endocrinology, 2005,146(3): 1236-1244.
[37]	YOSHIMURA T,YASUO S,WATANABE M,et al.Light-induced hormone conversion of T4 to T3 regulates photoperiodic response of gonads in birds[J].Nature,2003,426(6963):178-181.
[38]	IKEGAMI K,YOSHIMURA T.The hypothalamic-pituitary-thyroid axis and biological rhythms:The discovery of TSH’s unexpected role using animal models[J].Best Pract Res Clin Endocrinol Metab,2017,31(5):475-485.
[39]	任巧玲,马强,白红杰,等.雌激素及其对动物生殖功能的影响[J].河南农业科学,2006(4):114-115.REN Q L,MA Q,BAI H J,et al.Estrogen and its effect on animal reproductive function[J].Journal of Henan Agricultural Sciences, 2006(4):114-115.(in Chinese)
[40]	CAI Z X,GULDBRANDTSEN B,LUND M S,et al.Prioritizing candidate genes for fertility in dairy cows using gene-based analysis,functional annotation and differential gene expression[J].BMC Genomics,2019,20(1):255.
[41]	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):938.