论文部分内容阅读
巴西橡胶树（Hevea brasiliensis）是大戟科（Euphorbiaceae）、橡胶树属（Hevea）热带雨林乔木中唯一的栽培种,也是世界上广泛种植的产胶树种,其树皮所产胶乳是天然橡胶的原材料,具有很高经济和战略意义。目前国内外在橡胶树cDNA文库和蛋白组水平上开展了系列分子生物研究,取得了一定的进展,而有关橡胶树全基因组测序的相关研究起步较晚,仅于2013、2016和2017年完成对RRIM600、热研7-33-97和BPM 24三个橡胶树品种的全基因组测序,但也仅组装到Scaffold水平,尚未达到染色体水平。本研究在染色体水平上构建橡胶树单染色体文库,并对其进行测序,通过序列分析对单染色体的一些序列进行功能注释,旨在从单染色体水平上为橡胶树基因组研究提供新思路、新方法,为橡胶树分子辅助育种提供基础,研究取得了如下结果:（1）以巴西橡胶树’热研7-33-97’古铜期的嫩叶为材料,经酶解低渗去壁后制备染色体标本,通过微细玻璃针法从橡胶树中期细胞中共分离出40条单染色体。采用LA-PCR法随机对分离出的18条单染色体进行体外扩增,电泳检测后均有条带,条带大小范围在200～4000bp之间,从中随机选取5条染色体通过Southern杂交检测验证其来源,结果表明该5条染色体均来自橡胶树;将这5条单染色体进行荧光原位杂交验证并结合核型分析,确定其分别是热研7-33-97的第5、6、11、13和14号染色体。（2）对上述5条单染色体的第二轮LA-PCR产物进行纯化,并将纯化后的产物连接、转化、克隆最终构建出热研7-33-97第5、6、11、13和14号染色体文库,其克隆子数目分别是4.5×105个、5.11×105个、3.32×105个、3.01×105个、4.41×105个,插入片段大小范围分别在250～1100bp、250～1200bp、300～1100bp、300～950bp、200～1100bp,平均长度分别为 536 bp、481bp、388 bp、401bp、485 bp,文库的覆盖率分别为2.7倍、3.2倍、1.9倍、3.8倍、2.2倍,空载率分别为2%、3%、3%、1%、2%。（3）从热研7-33-97的第5、6、11、13和14号染色体的文库中,每个随机挑出100～140个左右的阳性克隆子进行测序,并将测序信息与NCBI中的巴西橡胶树RRIM600全基因组数据库（GCA001907995.1）、热研7-33-97全基因组数据库（GCA001654055.1）、EST数据库进行在线比对。比对结果表明,五条染色体文库序列中与橡胶树热研7-33-97全基因组数据库（GCA001654055.1）的高同源性分别有23、30、21、21、18条;与RRIM600全基因组数据库（GCA001907995.1）高同源性分别有23、28、22、20、18条;与橡胶树EST文库比较同源性大于90%的有20、74、40、38、57条,通过功能分析后多数单染色体序列主要与橡胶树耐盐性、耐割性、割面干涸病和抗冻性等有关。（4）对五个染色体文库中与橡胶树的全基因组数据库和EST数据库高同源的19条序列进行了 Blast2GO软件分析,通过比对发现有高同源性的序列有17条。对这17条序列用Blast2GO软件进行GO MAP比对,均得到了 GO术语注释结果,共获得23个GO ID,其中生物学途径10个、细胞组分6个和分子功能7个;对注释结果进行GO分级,共分为5级,有不同的功能,在生物学途径方面主要参与细胞过程、代谢过程、单有机体过程、生物体调控、生物合成过程、初级代谢等;在分子功能方面,主要参与催化活性、运输活性和转移酶活性等;在细胞组件部分主要是参与细胞部分、细胞、细胞器的组成。对经过GO注释后的序列进行Enzyme Code和KEGG代谢途径分析,除了第11号染色体文库序列未能注释到相关酶的功能,在5号、6号、13号及14号染色体文库均得到Enzyme Code 注释,其中 L-5Chr-5、L-5Chr-91 和 L-6Chr-48 序列上都得到一个Enzyme Code注释,所得注释的酶的ID均为2.7.7.49,为RNA-引导的DNA聚合酶;L-13Chr-25序列上都得到一个Enzyme Code注释,所得注释的酶的ID是2.3.1.12,为二氢硫辛酰赖氨酸残留乙酰转移酶;L-14Chr-73所得注释的酶的ID为3.6.1.15,为三磷酸核苷磷酸酶。但进行KEGG代谢途径分析,并未获得代谢途径。