论文部分内容阅读
北京鸭是世界著名的肉鸭品种,其中” Z型北京鸭瘦肉型配套系”具有饲料转化效率高、瘦肉率高的特点,是研究禽类骨骼肌发育的极好模型,但其骨骼肌发育研究尚不多见。但是,骨骼肌形成与再生的分子机制远非单个调节因子作用的加减,而是所有调节因子相互作用的网络。因此,以大规模的数据分析为基础寻找不同生长阶段北京鸭胸肌发育的差异表达基因及差异表达基因与北京鸭胸肌的关系研究是深度挖掘骨骼肌发育机制的基础。本研究应用抑制消减杂交(SSH)和RNA-seq技术分离鉴定了不同生长阶段北京鸭胸肌发育差异表达基因。分析了经典的骨骼肌发育调控基因(MSTN、IGF1和FOXO3)分子特征、它们在北京鸭胸肌发育过程中的表达变化规律及其与胸肌发育之间的关系。同时,开展了新发现的MUSTN1基因分子特征、时空表达模式及其与北京鸭胸腿肌发育的关系以及北京鸭胸肌内活体转染pcDNA3.1(+)-MUSTN1真核表达载体以探索MUSTN1基因对北京鸭胸肌发育的作用。试验方法及结果如下:1、应用SSH技术分离鉴定不同生长阶段北京鸭胸肌发育差异表达基因应用SSH技术从2周龄和6周龄北京鸭胸肌SSH文库中分离得到257个差异表达基因(165个已知基因和92个未知基因)。多数差异表达基因在6周龄北京鸭胸肌内高表达且与机体的能量代谢和蛋白合成有关。利用qRT-PCR分析了23个差异表达基因在北京鸭2、4、6、8周龄胸肌中的表达模式,通过聚类分析可以将这23个基因分为表达量持续下降、表达量先升高后下降和表达量持续上升等三种类别。2、应用RNA-seq技术分离鉴定不同生长阶段北京鸭胸肌发育的差异表达基因对2、4、6周龄北京鸭胸肌和皮脂组织进行cDNA文库构建及RNA-seq测序与分析,共得到21800万个高质量序列,占原始序列的95.76%。生物信息学分析发现这些高质量序列在鸭子参考基因上随机分布,覆盖了鸭子已注释基因的89.4%,序列的RPKM值与荧光定量值整体趋势一致。相同时期组织间及同一组织不同生长阶段间的差异表达基因数量从几百个到数千个,其中组织间差异表达基因远多于生长阶段间,而皮脂间差异表达基因多于胸肌间。GO和KEGG分析发现,差异表达基因被富集到不同GO条目中,而大量差异表达基因(如IGF1R、IGF1、MSTN、TGFβ3、TGFβ1、TGIF1、FOXO6、FOXO3等)被显著富集到肌肉发育及脂肪沉积的信号通路中,仅MAPK通路就富集了214个差异表达基因。3、肌肉抑制素基因(MSTN)分子特征及其多态性与北京鸭胸肌性状的关联分析应用qRT-PCR技术检测MSTN基因在第1d和2、4、6、8周龄北京鸭胸肌中的表达量发现:从1d到4周龄该基因表达量逐渐升高,但6周龄的表达量显著低于4周龄和8周龄的表达量。分子克隆及生物信息分析发现,北京鸭MSTN基因的CDS区长度为1128bp,其编码蛋白质与多数脊椎动物同源蛋白质具有很高的相似性。MSTN基因具有2个可变剪切体,分别编码375个氨基酸序列和251个氨基酸序列。通过对368只6周龄北京鸭该基因的克隆、SNP扫描、分型及关联分析发现:该基因含有3个SNP位点,其中T129C位点BB基因型鸭子的胸肌厚显著大于AA基因型鸭子的胸肌厚,T952C位点AB基因型鸭子的龙骨长显著长于BB基因型鸭子的龙骨长。4、IGF1和FOXO3表达模式及出雏后北京鸭胸腿肌发育关系研究北京鸭胸腿肌在第2d和2、4、6周龄石蜡切片表明,北京鸭胸肌发育落后于腿肌,其中胸肌在4~6周龄生长最快,而腿肌在2~4周龄生长速度最大。应用qRT-PCR技术检测IGF1和FOXO3在2d和2、4、6周龄胸腿肌中表达情况表明,IGF1分别在第6周和第2周胸肌和腿肌中的表达量最大,而FOXO3分别在第6周和第2周胸肌内表达量最小。5、MUSTN1基因分子特征及其与北京鸭胸肌发育关系研究经过对北京鸭MUSTN1基因cDNA和DNA序列克隆、拼接与分析发现,该基因编码一个含有78个氨基酸的蛋白质,由3个外显子和2个内含子组成,其启动子区包含许多转录因子结合位点。MUSTN1基因在北京鸭胸肌、腿肌中的表达量分别在5周龄和3周龄达到峰值。活体转染pcDNA3.1+MUSTN1表达载体到北京鸭胸肌后,该重组载体能够在北京鸭胸肌中高效表达,同时引起一些骨骼肌发育调控基因表达量的明显变化。综上所述,本研究分离鉴定了北京鸭胸肌发育不同阶段的差异表达基因。对部分经典的骨骼肌发育调控基因及新发现的MUSTN1基因进行了分子特征及其与北京鸭胸肌发育的关系研究。本研究结果将对深入探索北京鸭骨骼肌发育机制提供基础。