论文部分内容阅读
课题组前期开展的中国荷斯坦牛产奶性状全基因组关联分析利用传递不平衡(L1-TDT)和回归分析(MMRA)两种统计分析方法同时检测到38个显著SNP位点。本研究旨在基于该GWAS结果进一步筛选和鉴定奶牛产奶性状关键基因并进行功能验证。(1)位置候选基因鉴定及功能注释：基于牛基因组序列草图(Btau4.2),采用生物信息学和比较基因组学方法进行了上述38个显著SNPs的位置候选基因筛查和功能预测,最终鉴定到24个位置候选基因并确定了其物理位置、基因类型及功能注释；基于此,进一步从该24个位置候选基因中挑选了10个基因,利用实时荧光定量RT-PCR技术进行了组织表达谱分析,结果发现EEF1D和PDE9A基因在泌乳期中国荷斯坦牛乳腺组织中的mRNA表达量均明显高于其他7个组织(心肌,子宫,肾,肝,肺,卵巢和小肠)。结合前期GWAS结果,初步说明这2个基因与奶牛泌乳功能相关,将之确定为奶牛产奶性状候选基因,进行细胞水平功能验证。(2)EEF1D和PDE9A基因可变剪接鉴定：通过NCBI和UCSC数据库检索发现EEF1D和PDE9A基因的mRNA序列不完整并且多次改变,同时人和小鼠EEF1D和PDE9A基因mRNA序列的5’端在存在可变剪接,可形成不同转录本,因此推测奶牛的EEF1D和PDE9A基因存在可变剪接可能性较大。通过5’RACE实验,发现EEF1D基因在泌乳期中国荷斯坦牛乳腺组织中确实存在2种转录本,命名为EEF1Da和EEF1D6；EEF1Da(1202bp)和EEF1Db(2195bp)包含不同的第1外显子(exon1a:294bp；exonlb:1287bp),这2个第1外显子没有重叠,其余7个外显子序列和位置一致；利用实时荧光定量RT-PCR技术检测了不同转录本在奶牛乳腺组织的表达水平,结果显示EEF1Da为主要表达的转录本；8个组织的基因表达谱分析表明：EEF1Da在乳腺组织的表达量明显高于EEFIDb并且与EEF1D基因的表达趋势一致。另外,实验结果显不PDE9A基因在荷斯坦牛乳腺组织中只存在一种转录本,通过与已知序列拼接,获得了完整的PDE9A基因mRNA序列。(3)乳腺上皮细胞体外培养：采集健康的泌乳期中国荷斯坦母牛乳腺组织,通过组织块培养法和差酶消化法,获得了纯化的体外培养的乳腺上皮细胞。同时,采用乳汁分离法从收集的泌乳中后期荷斯坦母牛的乳汁中也获得了具有增殖能力的乳腺上皮细胞。经验证,培养的细胞具有生长和增殖能力,可表达上皮特有的角蛋白18基因、部分乳蛋白基因、产奶性状主效基因DAGT1和GHR以及EEF1D和PDE9A基因,适用于目的基因的功能验证。(4)RNA表达量干扰分析：基于EEF1D基因的完整mRNA序列,设计了4条siRNA序列,瞬时转染到奶牛成纤维细胞中,通过实时荧光定量RT-PCR技术检测EEF1D mRNA表达量,结果表明,3条siRNA显著降低EEF1D基因表达量,基于其中2条有效siRNA构建了相应的shRNA载体pGPU6-GFP-NEO-E1357和pGPU6-GFP-NEO-E1893,同时构建了相应shRNA慢病毒表达载体；由于PDE9A基因在牛成纤维细胞中表达量低,未作RNA干扰；另外,人工合成EEF1Da和PDE9A基因的CDS序列并连接到过表达载体,瞬时转染到人293FT细胞,RT-PCR技术检测到了EEF1Da和PDE9A基因表达,说明所构建的载体可用,同时也构建了含EEF1Da和PDE9A基因的CDS序列的慢病毒过表达载体。上述研究结果为深入进行EEF1Da和PDE9A基因的功能验证和奶牛产奶性状关键基因鉴定提供了前期试验基础。