伊犁鹅产蛋前后各期卵巢转录谱的构建及卵泡发育相关基因的分析

【目的】构建产蛋前后各期伊犁鹅卵巢转录组文库，结合生物信息学分析揭示不同产蛋期伊犁鹅卵巢组织差异表达基因，并鉴定出影响鹅卵巢发育的关键基因，为伊犁鹅的繁殖调控提供理论参考。【方法】选择处于开产期（KL组）、产蛋期（CL组）及休产期（XL组）的伊犁鹅各4只，屠宰后取其卵巢组织，用于构建卵巢转录组文库。通过新基因挖掘、基因差异表达、基因注释以及蛋白互作网络分析筛选出与鹅卵泡发育相关的候选基因；随机挑选8个差异表达基因，应用实时荧光定量PCR验证其表达情况。【结果】伊犁鹅卵巢组织学结果表明，在开产期时，伊犁鹅卵巢表面有大量初级卵泡，而产蛋期卵巢则显示出卵泡的层级性，在休产期卵巢中可观察到卵泡出现向内凹陷、闭锁的现象。通过转录组测序（RNA-Seq），从构建的12个伊犁鹅卵巢cDNA文库中获得有效读数57 811 186~85 328 377条，Q30值均>93.38%，每个样品所产的测序读数于鹅参考基因组上的比对率在82.79%~89.24%。在卵巢中注释的新基因共有1 112个，单核苷酸多态性（SNP）位点总数为1 642 273~2 425 069个；SNP和插入缺失（InDel）均主要注释于内含子区域；各时期的可变剪切类型主要集中为最后1个外显子可变剪切（TSS）及第1个外显子可变剪切（TTS）。在KL vs CL、XL vs CL及XL vs KL组中分别有337、1 136和525个差异表达基因，共有差异表达基因为α2A肾上腺素能受体（ADRA2A）、表皮蛋白（CP）、非转移性黑色素瘤糖蛋白B （GPNMB）及α-1-抗胰蛋白酶样（LOC106033756）。GO功能富集分析发现，差异表达基因主要富集在肽基酪氨酸磷酸化的正调控、细胞黏附及质膜外侧等过程。KEGG通路富集分析发现，差异表达基因主要显著富集于神经活性配体-受体相互作用、ECM-受体相互作用及类固醇生物合成等。结合蛋白互作网络分析，筛选到与鹅卵巢发育相关的潜在调控因子Bruton’s酪氨酸激酶（BTK）、血小板源性生长因子受体α（PDGFRA）、整合素β3（ITGB3）等。实时荧光定量PCR结果显示，RNA-Seq结果准确可靠。【结论】本研究揭示了产蛋前后各期伊犁鹅卵巢组织中的基因表达差异，筛选到影响伊犁鹅卵巢发育的神经活性配体-受体相互作用、ECM-受体相互作用、类固醇生物合成等重要通路与BTK、PDGFRA和ITGB3等关键候选基因，为了解鹅卵巢组织调控产蛋性能的分子机制提供了理论支撑。

Construction of Ovarian Transcriptional Profiles and Analysis of Follicle Development-Related Genes in Yili Geese at Various Stages Before and After Egg Laying