营养诱导舍饲哈萨克羊非繁殖季节卵巢组织中miRNAs差异表达分析

旨在利用Solexa测序技术分析、挖掘营养诱导非繁殖季节绵羊发情期和乏情期卵巢组织micro RNAs(miRNA),为进一步研究特定miRNA参与非繁殖季节绵羊发情的分子调控机制提供理论基础。以不同营养水平饲喂的哈萨克绵羊作为研究对象,在非繁殖季节前后分别屠宰已鉴定的3只发情(EN)和乏情(AN)母羊,分别采集卵巢(O)组织,构建两个小片段RNA文库(OEN和OAN),进行高通量测序和生物信息学分析,实时荧光定量PCR(qRT-PCR)验证差异表达的miRNAs在哈萨克羊卵巢中表达水平。结果显示,两文库分别获得8 665 978(发情期)和9 071 102条(乏情期)clean reads,且绝大多数小RNA序列长度为21~23nt。两文库相比较,发现有9个已知和104个未知miRNAs表达差异显著(P0.01)。利用qRT-PCR随机选择的6个显著差异表达的miRNAs进行验证,其表达水平和测序分析结果一致。KEGG通路分析表明,这些差异miRNAs的靶基因主要参与了甘氨酸、丝氨酸和苏氨酸代谢,氨基糖和核苷酸糖代谢,精氨酸和鸟氨酸代谢及甘油磷脂代谢等营养调控通路。结合靶基因预测及通路富集分析,推测差异表达的miRNAs是通过营养代谢通路对绵羊的非季节性发情起重要作用。     

Smad2基因在敖汉细毛羊不同发情时期卵巢中的表达分析

Smad2基因是Smads蛋白家族的成员之一,Smad家族中的信号分子对绵羊的繁殖、卵泡发育有重要影响。本实验选取24只3~4岁健康空怀敖汉细毛羊,采用荧光定量PCR、蛋白免疫印迹和免疫组织化学技术,测定Smad2基因在绵羊不同发情时期(乏情期、发情间期、发情前期、发情期)的表达量并分析其在卵泡中表达规律。结果表明:Smad2基因m RNA和蛋白均在绵羊发情前期表达量最高,m RNA水平极显著高于乏情期、发情间期和发情期(P0.01);而发情前期蛋白水平极显著高于乏情期和发情期(P0.01),虽高于发情间期,但差异不显著(P0.05);Smad2蛋白主要在卵巢颗粒细胞、卵泡膜细胞以及卵母细胞表达。据此,Smad2基因可能会促进敖汉细毛羊的发情启动。     

小尾寒羊高繁殖力候选基因ESR在同期发情处理过程中的表达分析

试验选择高繁殖力小尾寒羊为试验动物,采用公认的CIDR结合外源激素处理使其同期发情,与自然发情及乏情的小尾寒羊作对比,利用反转录荧光定量RT-PCR和Western blotting技术研究雌激素受体(estrogen receptor, ESR)基因在卵巢细胞中的表达差异,同时对其启动子区的甲基化岛进行了分析。结果发现,自然发情绵羊和外源激素处理发情的羊在发情表现和时间上无明显差别;在转录水平的检测结果上也无显著性差异(P＞0.05);但是它们显著高于在乏情期绵羊卵巢细胞内的表达(P＜0.05)。尽管激素处理的同期发情羊卵巢细胞内的ESR蛋白表达量显著高于乏情期绵羊,但是仍然显著低于自然发情期绵羊卵巢细胞的表达(P＜0.05)。此外,自然发情和激素处理的绵羊ESR基因启动子区CpG岛甲基化水平分别为0和10.77%,而发情期绵羊ESR基因启动子处于超甲基化状态（92.30%）。以上结果说明,外源激素处理的同期发情小尾寒羊尽管发情表现、时间等与自然发情的绵羊类似,但是ESR基因在蛋白水平上差异显著,且其启动子区CpG岛甲基化水平也存在显著性差异,这也许是目前胚胎移植后着床率低下的原因之一。     

小尾寒羊TGF-β3基因在同期发情处理过程中的表达分析

为了探索胚胎移植过程中妊娠率的根本原因。本研究以高繁殖力小尾寒羊为试验对象,采用公认的CIDR结合外源激素处理使其同期发情,与自然发情及乏情的小尾寒羊作对比,分析了转化生长因子TGF-β3基因在卵巢细胞中的表达差异,同时对其启动子区的甲基化岛进行了分析。结果发现,自然发情绵羊和外源激素处理发情的羊在发情表现和时间上无明显差别,但是TGF-β3基因在转录水平和蛋白水平上自然发情绵羊显著高于外源激素处理发情的羊,而同期处理发情的羊显著高于乏情期绵羊卵巢细胞内的表达。自然发情和激素处理的绵羊TGF-β3基因启动子甲基化水平分别为25.8%和53.3%,而乏情期绵羊TGF-β3基因启动子处于超甲基化水平达到70.0%。以上结果说明外源激素处理的同期发情小尾寒羊尽管在发情表现、时间等与自然发情的绵羊类似,但是胚胎发育和着床重要调控基因TGF-β3的表达和甲基化模式存在显著性差异这也许是目前胚胎移植后着床率低下的原因之一。     

Smad1基因在敖汉细毛羊发情不同时期卵巢中的表达规律研究

为了探索Smad1 mRNA及其蛋白在敖汉细毛羊不同发情时期卵巢中表达量的变化规律,试验采用实时荧光定量PCR、免疫组化方法分析绵羊不同繁殖状态下(乏情期、发情间期、发情前期、发情期)Smad1基因在卵巢中mRNA表达量的变化规律及卵巢中该基因蛋白的表达,并进行图像分析。结果表明:发情期Smad1基因表达量最高,极显著高于其他三个时期(P0.01);发情前期Smad1基因表达量显著高于乏情期和发情间期(P0.05)。在绵羊不同发情时期Smad1蛋白主要表达于乏情期原始卵泡、次级卵泡的颗粒细胞、发情间期次级卵泡的卵泡膜细胞,在发情前期与发情期卵泡中未见表达。卵巢中Smad1基因可能促进黄体溶解,对卵泡发育、卵子的成熟以及释放起作用,与绵羊发情的启动有一定关系。     

发情和乏情初产母猪下丘脑-垂体-卵巢轴差异表达基因筛选及分析

研究旨在分析发情和乏情初产母猪下丘脑-垂体-卵巢轴基因表达的差异,探讨母猪乏情调控的分子机制。选用6头健康的断奶初产母猪,分为发情组（3头）和乏情组（3头）,屠宰后分别采集下丘脑、垂体、卵巢进行转录组测序（RNA-Seq）、GO功能富集、KEGG通路及相关蛋白的Western blotting分析。结果发现,发情组下丘脑、垂体、卵巢分别获得52 432 800、52 573 730和52 209 252条clean reads,与猪参考基因组（Sus scrofa 11.1）的比对率分别为78.69%、81.49%和79.26%;乏情组下丘脑、垂体、卵巢分别获得52 516 724,52 476 820和52 195 962条clean reads,与猪参考基因组的比对率分别为82.38%、83.05%和80.20%。与发情组母猪相比,乏情组母猪下丘脑中共有710个差异表达基因,其中上调基因392个,下调基因318个;垂体中共有707个差异表达基因,其中上调基因283个,下调基因424个;卵巢中共有956个差异表达基因,其中上调基因635个,下调基因321个。3种组织中的共有差异表达基因为36个。与母猪情期调控相关的亲吻素-1（KISS1）、G-蛋白偶联受体54（GPR54）、速激肽3（TAC3）、速激肽3受体（TACR3）、17-α-羟化酶/17,20碳链裂解酶（CYP17A1）、芳香化酶（CYP19A1）、类固醇激素合成急性调节蛋白（STAR）、促性腺激素释放激素受体（GnRHR）和雌激素受体1（ESR1）基因在乏情组母猪体内显著下调（P<0.05;P<0.01）。Western blotting分析结果显示,TAC3、TAC3R、KISS1和GPR54相关蛋白在乏情组母猪下丘脑中也显著下调（P<0.05）。GO和KEGG通路分析发现,差异表达基因显著富集于正调控胆固醇酯化反应、卵巢类固醇生成、GnRH信号通路、FoxO信号等一些与类固醇激素生成和卵泡发育相关的信号通路。本研究结果表明,发情和乏情初产母猪下丘脑-垂体-卵巢轴上与情期调控相关的基因存在差异表达,尤其是Kisspeptin/GPR54和TAC3/TACR3两大系统的mRNA及蛋白表达水平在乏情组母猪下丘脑中均显著下调,推测可能是Kisspeptin/GPR54和TAC3/TACR表达受损导致了下丘脑调节功能障碍,进而导致了初产母猪乏情。该研究为揭示初产母猪乏情分子调控机制提供了重要支持,为今后利用分子技术改善母猪乏情问题提供了重要理论依据。     

BAD基因在敖汉细毛羊发情不同时期卵巢表达规律及对生殖激素的影响

本实验研究旨在探究BAD基因在敖汉细毛羊中的不同发情时期卵巢中表达量的变化规律以及对生殖激素的影响,为探讨其与绵羊发情的关系奠定基础。首先利用放射免疫方法测定敖汉细毛羊血液中促卵泡激素(FSH)、促黄体生成素(LH)、雌二醇(E_2)和孕激素(P_4)的浓度,然后利用qRT-PCR检测乏情期、发情前期、发情间期和发情期卵巢BAD基因的表达量变化。结果表明:卵巢中BAD基因在发情前期的表达量极显著高于乏情期、发情期和发情间期(P0.01),在发情期的表达量极显著高于乏情期和发情间期(P0.01),乏情期和发情间期的表达量基本相同(P0.05);E_2和LH浓度在BAD基因表达量最高的发情前期维持较高水平,P_4的浓度维持较低水平。综上可知,BAD基因的表达可能促进了E_2和LH浓度升高,抑制了P_4产生,进而影响细毛羊发情,可为绵羊发情的研究提供借鉴。     

绵羊不同繁殖状态下SIX1和TEF基因的表达特征分析

试验旨在探究SIX1和TEF基因在苏尼特羊(Sunite sheep,SNT)和小尾寒羊(Small Tail Han sheep,STH)相关组织中的表达特征,有助于揭示上述2个基因在绵羊季节性发情和繁殖调控中的重要作用。选取季节性发情的SNT母羊在短光照(模拟繁殖季节)和长光照(模拟休情期)条件下及常年发情的STH母羊在不同繁殖时期(卵泡期和黄体期)的下丘脑等10种组织,利用qPCR技术分析上述不同繁殖状态下各组织中SIX1和TEF基因的相对表达量。结果表明,SIX1基因在SNT和STH的垂体组织中均高表达,其它组织中微弱表达;TEF基因在2个绵羊品种的多个组织中广泛表达;SNT垂体中TEF基因在短光照条件下其表达量显著高于长光照条件(P0.05),STH垂体中TEF基因在卵泡期其表达量显著高于黄体期(P0.05);SNT子宫体中TEF在长光照条件下其表达量显著高于短光照条件(P0.05),STH子宫体中TEF在黄体期其表达量极显著高于卵泡期(P0.01)。研究结果显示,SIX1和TEF基因表达在绵羊垂体中发挥重要作用,2个基因在SNT垂体中的表达变化趋势与已知的长光照诱导基因EYA3的表达变化不同,暗示绵羊垂体中SIX1和TEF基因不是通过转录水平变化来参与季节性发情上游基因的调控;TEF基因可能参与绵羊不同繁殖状态下子宫生理变化的调控。本研究为深入探究这2个基因在绵羊繁殖性能调控方面的作用奠定了基础。     

染料木素对舍饲哈萨克羊乏情期诱导发情的作用

为了解染料木素对舍饲哈萨克羊乏情期诱导发情的调理作用,选用体况接近的15月龄半同胞哈萨克母羊18只,随机分成3个组。对照组饲喂基础日粮,高剂量组和低剂量组每天每只分别给予40、20 mg的染料木素,给药24 d。结果显示:高剂量组发情羊所占比例为50.00%,低剂量组16.67%,对照组无发情绵羊;与对照组和低剂量组相比,高剂量组显著提高绵羊体内雌激素水平(P0.05),与对照组相比低剂量组显著提高绵羊体内雌激素水平(P0.05);与对照组相比,试验组绵羊血清中总蛋白、白蛋白和血清钙含量均显著增加(P0.05)。结果表明:一定剂量的染料木素可以诱导舍饲哈萨克羊在乏情期发情。     

甘加藏羊发情周期腺垂体和血浆中GTH含量动态变化

为了探讨甘加藏羊发情周期腺垂体和血浆中FSH和LH含量动态变化规律及排卵的关系。选取2.5~5岁,健康未孕雌性甘加藏羊30只,采集发情周期各阶段和乏情期腺垂体和血浆,采用酶联免疫吸附试验(ELISA)测定发情前期、发情期、发情后期、间情期及乏情期腺垂体和血浆中FSH和LH的含量。结果表明,甘加藏羊发情周期内腺垂体FSH浓度在发情期(6.873IU/L)和发情后期(6.863IU/L)较低,与间情期最大(8.228IU/L)差异显著(P<0.05),发情前期又降低(8.096IU/L);腺垂体LH浓度发情期(9.39IU/L)和发情后期(9.585IU/L)较低,间情期最大(10.629IU/L),发情前期(9.7325IU/L)降低,各时期差异不显著(P>0.05);发情周期内血浆中FSH浓度在发情前期(2.599IU/L)达到峰值,其他时期均较低,发情前期与其他时期差异显著(P<0.05);血浆LH浓度在发情期(3.866IU/L)达到峰值,发情后期(3.179IU/L)和间情期(3.292IU/L)降低,各时期差异不显著(P>0.05);并且FSH和LH在发情周期内均出现4个波峰:乏情期腺垂体和血浆中FSH和LH的浓度均低于发情周期内FSH和LH各自的基础值。FSH和LH含量在垂体和血浆的这些差异性应该是两者协同调控藏羊发情和排卵的主要原因。     

多浪羊与小尾寒羊皮下脂肪组织转录组分析

旨在探索多浪羊（D组）与小尾寒羊（X组）皮下脂肪组织中的特异性表达基因,并研究其潜在的作用,为理解绵羊脂肪组织发育规律以及对脂代谢相关疾病的预防和治疗研究提供依据。本研究选取脂肪沉积能力存在差异、健康无病、体况良好、种内个体体重相近（约50 kg）的雌性成年多浪羊和小尾寒羊为试验材料,分为D组（试验组）和X组（对照组）,每组3个重复,采集位于背最长肌的皮下脂肪组织,应用RNA-Seq技术和生物信息学方法进行转录组测序并对结果进行分析。以|Fold change|≥2,P adjust≤0.05为标准筛选差异表达基因,通过对差异表达基因进行GO功能注释和KEGG通路富集分析,得到与脂肪沉积和脂代谢有关的差异基因。为了验证测序数据的可靠性,本研究随机选取了6个差异表达基因进行qRT-PCR验证。结果显示,在6个样本中共检测到38 672个已知的mRNAs,新的mRNAs为1 606个,在两组中共有839个差异表达基因,其中有320个差异基因在多浪羊组中上调表达,有519个差异基因在多浪羊组中下调表达。通过GO功能注释分析发现,差异表达基因主要参与脂质分解代谢过程、脂质生物合成过程、脂质分解代谢负调控过程、MAPK级联反应调控、对甘油三酯的反应等生物学过程。KEGG通路富集结果显示,差异表达基因显著富集到了PI3K-Akt、MAPK、胰岛素以及PPAR等信号通路中。qRT-PCR结果与测序结果一致,表明测序结果可靠。通过对多浪羊和小尾寒羊皮下脂肪组织进行转录组测序以及生物信息学分析,筛选到与脂肪沉积和脂代谢相关的差异表达基因,这些基因主要参与脂质生物合成、脂质代谢等过程,其中COL1A1、AKT2、SCD、LPL、PCK1与PPP2R5A可能在多浪羊与小尾寒羊的皮下脂肪组织的沉积与代谢中发挥重要作用。     

不同繁殖季节绵羊卵巢组织中piRNAs序列特征及表达差异分析

为了分析绵羊卵巢组织中piRNAs的序列特征及不同繁殖季节绵羊卵巢组织中piRNAs的表达差异,本试验采用Solexa测序技术对休情季节（休情期）和发情季节（卵泡期和黄体期）滩羊卵巢组织进行高通量测序,采用生物信息学方法筛选绵羊卵巢中的piRNAs,对piRNAs的长度分布、首位碱基偏好、基因组来源、染色体分布、链特异性、ping-pong循环特征、piRNAs比对重复序列类型进行分析,并对发情季节与休情季节绵羊卵巢之间差异表达的piRNAs进行靶基因预测及其靶基因的GO功能和KEGG Pathway富集分析。结果显示,绵羊卵巢中不同长度的piRNAs首位碱基偏好性不一致;3个时期（休情期、黄体期和卵泡期）的piRNAs均主要比对到基因内含子、编码区及lncRNA区域,而比对到重复序列的piRNAs相对较少;绵羊卵巢piRNAs在染色体上呈不均匀分布;绵羊piRNAs簇存在很强的链特异性,但首位碱基偏好性不明显;3个时期的piRNAs均无明显的ping-pong循环特征。KEGG富集分析发现,差异piRNA靶基因主要富集在RNA运输通路（RNA transport pathway）、嘌呤代谢通路（purine metabolism pathway）、黏着斑通路（focal adhesion）、PI3K-Akt信号通路（PI3K-Akt signaling pathway）、Hippo信号通路（Hippo signaling pathway）等与绵羊繁殖相关的通路。本研究揭示了不同繁殖季节绵羊卵巢组织中的piRNAs序列特征,暗示绵羊卵巢组织中的piRNAs可能通过BAD等基因和黏着斑通路等调控绵羊季节性繁殖,为深入解析绵羊季节性发情分子机制提供一定参考。     

黄体期和卵泡期小尾寒羊DIO2与DIO3组织表达

为分析DIO2和DIO3基因在小尾寒羊不同繁殖时期(黄体期和卵泡期)下丘脑-垂体-性腺轴各组织中的表达差异,阐明这2个基因在绵羊发情转换中的表达模式,实验采用实时荧光定量PCR(Real-time PCR)技术对比分析DIO2和DIO3基因在黄体期和卵泡期小尾寒羊下丘脑-垂体-性腺轴各组织中的表达差异。结果表明:DIO2和DIO3基因在下丘脑、垂体、松果体、大脑、小脑、卵巢、子宫、输卵管、肾脏、肾上腺10种组织中均表达;DIO2基因在黄体期和卵泡期的垂体组织中表达量显著高于其他组织(P<0.05),其中黄体期垂体、子宫、下丘脑、松果体、输卵管和卵巢DIO2的表达量显著高于卵泡期(P<0.05);DIO3在卵泡期松果体、下丘脑、子宫、垂体和输卵管的表达量显著高于黄体期(P<0.05)。综上,DIO2能抑制绵羊发情,而DIO3促进发情。     

非繁殖季节中国美利奴羊血清E2、P4、FSH、LH激素变化规律研究

为了在非繁殖季节判断绵羊的发情状况,试验以新疆高海拔地区中国美利奴羊为研究对象,采用酶联免疫测定(ELISA)法测定绵羊血清中雌二醇(E2)、孕酮(P4)、促卵泡素(FSH)、促黄体生成素(LH)的含量,研究激素诱导发情羊、发情症状不明显羊、不发情羊体内激素分泌规律的差异。结果表明:激素诱导发情羊的E2含量比不发情羊高,差异显著(P<0.05);但是在P4、FSH、LH水平上,激素诱导发情羊要低于不发情羊,差异不显著(P>0.05)。发情羊与发情不明显羊在4种激素含量上差异不显著(P>0.05),说明公羊试情存在一定缺陷,而通过测定E2含量判定绵羊发情状态将成为未来研究的重点。     

牦牛泌乳期垂体组织比较转录组学分析

试验旨在了解不同产奶量牦牛泌乳期垂体组织间的转录组差异,为进一步阐述牦牛垂体组织调控泌乳的机制提供参考。提取牦牛垂体组织总RNA,并应用RNA-Seq技术对高产奶量牦牛和低产奶量牦牛泌乳期间垂体组织转录组进行高通量测序及分析。结果显示,通过比较分析高产奶量牦牛和低产奶量牦牛垂体组织转录组数据,筛选出114个差异表达基因,其中55个表现为上调,59个表现为下调。进一步功能分析表明,这些差异基因涉及多种GO分类及KEGG通路,其中GO分析显示,差异基因与许多和氨基酸代谢相关的生物学过程存在紧密关联;KEGG通路分析显示,最为富集的是细胞黏附分子通路,且金黄色葡萄球菌感染和利什曼病等大量与免疫或疾病相关的通路也出现显著富集。本研究对比了高产奶量牦牛和低产奶量牦牛泌乳期垂体组织转录组数据,筛选并分析了相关差异表达基因,为提高牦牛产奶量的研究提供新的思路。     

基于转录组数据挖掘藏羊立毛肌发生的关键基因

旨在探究调控藏羊皮肤内立毛肌发生的分子机制，挖掘影响立毛肌发生的关键基因。本研究采集了30个75~110日龄健康情况良好的藏羊胚胎，收集背部皮肤组织，制作石蜡组织切片，利用HE染色观察背部皮肤内不同日龄立毛肌的形态学变化特征，推测出立毛肌发生的关键时期为胚龄75~85天（E75~E85）；分别挑选胚龄大约在75和85天各3个背部皮肤组织样品，分样品提取RNA，利用Illumina Hiseq平台进行转录组测序，对测序数据进行质控、过滤、比对，筛选得到差异表达基因，并对差异基因进行功能注释、富集分析、蛋白网络互作和RT-qPCR验证，挖掘调控立毛肌早期发育的关键候选基因。测序结果共得到1 159个差异表达的基因（P<0.05），其中900个基因上调，259个基因下调。随机选择6个差异表达的基因进行RT-qPCR验证，表达趋势与转录组测序结果一致，表明测序结果可靠。通过对差异基因进行GO功能注释和KEGG富集分析，筛选到与立毛肌发育相关的5个生物学过程条目和5个信号通路共47个非冗余基因（如BAMBI、TNNI3、HOXA9等）。本研究表明，藏羊立毛肌开始发育的关键时期约为E75~E85天，该时期内SPP1、BAMBI、TNNI3、HOXA9、SOX15可能为影响藏羊立毛肌发生的关键候选基因。     

绵羊EYA3基因表达及其多态性与季节性发情之间的关系

为探究绵羊EYA3基因表达及其多态性与季节性发情之间的关系,本实验通过采用实时荧光定量PCR(qPCR)检测EYA3基因在常年发情绵羊(小尾寒羊)和季节性发情绵羊(苏尼特羊)的大脑、下丘脑和垂体组织中的表达情况,同时利用Sequenom MassARRAY~?SNP技术检测2组绵羊(常年发情组:小尾寒羊、策勒黑羊和湖羊;季节性发情组:滩羊、苏尼特羊和草原型藏羊)EYA3基因SNP位点(g.238191128G>C)的多态性,并与绵羊季节性发情性状进行关联分析。结果表明:无论在小尾寒羊还是苏尼特羊,EYA3基因在各组织中均广泛表达,在苏尼特羊垂体中表达较高(P<0.01);且长光照条件下季节性发情的苏尼特羊垂体中EYA3基因表达量显著高于小尾寒羊(P<0.01)。绵羊EYA3基因g.238191128G>C位点均存在GG、GC和CC3种基因型。g.238191128G>C位点在小尾寒羊、湖羊、策勒黑羊、滩羊、苏尼特羊和草原型藏羊中表现为低度多态(PIC<0.25);在小尾寒羊、策勒黑羊、滩羊、苏尼特羊和草原型藏羊5个群体中均处于哈代温伯格平衡状态(P>0.05);并且该位点基因型频率和等位基因频率在常年发情和季节性发情绵羊组间差异均达到极显著水平(P<0.01)。综上,EYA3基因表达与g.238191128G>C位点与绵羊季节性发情性状存在一定关联。     

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

旨在分析不同繁殖周期绵羊卵巢组织中长链非编码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,并进行差异分析,为揭示绵羊繁殖能力的分子机制提供依据。     

云南半细毛羊精子冷冻前后转录组表达差异分析

【目的】 本研究旨在对云南半细毛羊精子冷冻前后差异表达的基因与精子冷冻损伤的关系进行研究。【方法】 采用电刺激法采集3只云南半细毛羊精液,各分成2份,一份直接用于提取精子RNA,另一份经过冷冻保存7 d后再进行精子RNA提取。利用获得的精子RNA,采用单细胞转录组测序技术分别检测新鲜精子和冻融精子的mRNA转录组,构建mRNA转录组文库,用Hiseq2500测序仪进行测序,对测序结果做进一步过滤处理,然后对差异表达基因进行GO功能注释和KEGG通路分析。【结果】 6份精液样本共产生623 399 754条原始测序序列,过滤处理后得到514 313 802条（82.50%）高质量的Clean reads。生物信息学分析结果表明,冷冻前后差异表达显著的基因共1 213个。其中,解冻后表达量下调的有1 208个,上调的有5个。通过GO功能注释和KEGG通路富集分析,差异表达基因分别归类于56个GO分类,参与40个信号通路。其中与精子功能密切相关的主要包括精子细胞过程、代谢过程、应激反应、生殖、质膜和氧化还原反应;精子冷冻前后差异表达基因参与的信号通路主要有炎症、细胞因子-细胞因子受体相互作用和细胞凋亡等通路。【结论】 云南半细毛羊精子的冷冻损伤可能与获得的差异表达基因参与的生物过程有关。此外,获得的差异表达基因也可作为分子标记物应用于绵羊冻精品质的评价。