论文部分内容阅读
明晰适应性表型变异的分子基础是理解生物多样性形成与维持的核心。根据Gommans等的理论，遗传变异分子机制如果能提高生物体产生新表型以适应环境变化的能力，该机制便可能会增加物种的可进化性，如此，高遗传变异水平或表型可塑性的体系往往会更易被选择。而相较于基因组变异“全或无”的特性，转录和转录后水平修饰的多样度和灵活性更高，它们可能正是真核生物的可进化能力显著强于原核生物的重要体现。家鸡自其野生祖先——红原鸡(red junglefowl，RJF)被驯化后，演化出相当丰富的品种，是表型多样度最高的鸟类。但鸡基因组很小(不到1G)，要在仅仅8，000年的分歧时间里产生如此高的表型多样度，若只依靠基因组变异，显然是不够的。高通量测序技术的迅速发展，显著降低了分析的成本和周期。故此，本论文基于大范围基因组和转录组数据，构建鸡基因间长链非编码RNA(long intergenic non-coding RNA，lincRNA)和RNA编辑谱，分析它们的组织特性、功能以及其在家鸡驯化过程中的作用，并探讨了云南特有品种鸡——瓢鸡的无尾性状的分子遗传机制。 基于821个RNA-Seq(RNA Sequencing)样本，我们鉴定得到了4，754个lincRNA基因，并分析其诸如表达、组织特异性以及转录本长度和外显子数等特征，发现lincRNA的表达水平要比蛋白编码基因低得多，但组织特异性更高，转录本短而外显子数少。通过比较转录组分析，我们找到了68个在不同条件下差异表达的lincRNA。而利用已有的全基因组关联研究（genome-wide association studies，GWAS）数据，我们发现47个与特定表型显著相关的lincRNA。我们的家鸡与红原鸡间的大范围群体基因组比较分析结果揭示，419个lincRNA可能在家鸡驯化过程中受到强烈的人工选择压力，它们或许通过影响鸡的行为和繁殖力等来参与新品种的形成。 利用96个鸡RNA-Seq样本，我们构建了鸡的RNA编辑谱。分析不同组织的整体编辑水平模式发现，鸡里，肝脏中RNA编辑水平最低，骨骼肌中最高，这与哺乳动物肌肉中编辑水平相对较低的结论迥然相异。家鸡和红原鸡间RNA编辑水平比较分析揭示，两个群体间具有非常相似的编辑模式，群体间差异编辑的位点不到500个，仅约占所有A-to-I编辑位点的1％。两个群体的基因组比较分析结果亦表明，编辑位点在家鸡驯化过程中受到的选择压力很弱。而PhastCons值显示RNA编辑位点在鸟类中是保守的，但在其它脊椎动物中发生了快速进化。ADAR酶的RNA编辑贡献性分析也得到了与哺乳动物截然不同的现象:鸡里，ADARB1和ADARB2基因的表达与编辑水平呈显著正相关，而ADAR基因与编辑水平呈显著负相关;哺乳动物中，ADAR1和ADAR2酶是RNA编辑的主要贡献者，而ADAR3是RNA编辑的主要抑制子。不过，与哺乳动物中相一致的是，我们在鸡中也找到许多除ADAR酶以外的RNA编辑调节基因。例如，编码小亚基加工体的非核糖体蛋白基因UTP18可能通过影响核糖体的组装来抑制RNA编辑。 基于21个不同胚胎发育时期的尾部转录组的比较分析，我们发现了许多与鸡尾部发育显著相关的差异表达基因和共表达模块，例如，MYL3和MYH7B参与的肌发生，或PAX9和PAX1涉及的脊柱形成。通过对20个成年瓢鸡与98个成年正常有尾鸡的共118个基因组重测序数据的群体基因组比较分析，我们找到两个在瓢鸡和正常有尾鸡间呈现显著遗传分化的区域，chr2:84.5Mb-89Mb和chr24:6.12Mb-6.18Mb。它们包含了诸如TERT,PDCD6等许多与骨发生、凋亡及钙信号通路等相关的基因。而尾部胚胎发育显微观察结果显示，瓢鸡胚胎尾部缺陷在胚胎发育第四天前就已形成。因此，我们推断，调节元件所受选择压力可能导致这些区域基因在胚胎发育早期时的活性改变，并通过促凋亡、BMP或者NF-κB等信号通路破坏尾牙的间质维持，最终导致不正常的尾部发育。 本论文通过整合大范围转录组和基因组数据分析，揭示了转录及转录后水平修饰在家鸡驯化过程中的重要作用，完善了现有lincRNA和RNA编辑谱及其功能注释，并探讨了瓢鸡无尾性状发生的可能遗传机制。我们的工作将为今后复杂表型适应性进化机制的相关研究提供思路和奠定基础。