《锈菌类真菌基因组构造分析.docx》由会员分享，可在线阅读，更多相关《锈菌类真菌基因组构造分析.docx（12页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、锈菌类真菌基因组构造分析锈菌种群庞大，能够引起很多重要经济作物和林木病害，严重威胁全球粮食和林业生产安全。全基因组分析为锈菌基因功能研究、毒性变异研究及锈菌演化规律研究提供了重要基础，为制定锈病有效防控策略和创制抗锈新材料提供理论根据。本文综述了目前锈菌全基因组分析领域的进展，对锈菌的基因组构造、基因组成、基因组变异等特征进行了归纳分析，对基因组变异与其专性寄生特性的关系、基因组变异对其毒性变异的可能影响等进行了阐述。基因组学将为最终揭示锈菌生活史复杂性和毒性高度变异性的根本成因提供有力工具。关键词：锈菌，基因组，毒性，变异，基因组测序，进化锈菌属于担子菌门真菌，已发现130多属，5000余种
2、，是高等专性活体寄生真菌。寄主植物被锈菌侵染后，在叶片、叶鞘、茎或果实上长出大量的橙色、褐色或者红色锈状孢子堆，有的还可在枝干上引起肿瘤、粗皮、丛枝、曲枝等症状，或造成落叶、焦梢、生长不良等。严重时孢子堆密集成片，植株因体内水分大量蒸发而迅速枯死。锈病是世界范围广泛分布的农作物和森林作物重大灾祸性病害Petersen1974；Aniksteretal.2004。近一个世纪以来，锈病中最具代表性的小麦锈病频繁爆发流行，造成包括我国在内的全球小麦主产国毁灭性损失。研究表明，锈菌的毒性变异是导致小麦品种抗病性的“丧失，进而造成病害大规模流行的根本原因。因而，揭示锈菌毒性变异的基础成因是解决品种抗性“
3、丧失，实现病害持久控制的必由之路。生物遗传变异的重要表现是其各种途径导致的基因组变异，包括基因组序列变异和构造变异。序列变异主要有单核苷多态性、短序列插入与缺失、微卫星或简单重复序列和转座子等类型；构造变异主要包括大片段插入与缺失、基因有/无变异和基因拷贝数变异，表现为基因组区域大范围的插入、缺失、复制、倒置inversion和移位translocation。研究表明，植物病原菌基因组水平的各类变异能够最终导致病原菌毒性变异。因而，全基因组测序分析和比拟基因组分析是揭示锈菌类真菌演化规律和毒性变异分子机制的有效途径。近年来，基因组学的迅猛发展为最终揭示锈菌生活史复杂性和毒性高度变异性的根本成因
4、提供了有力工具。本文综述了目前锈菌全基因组测序分析领域的进展，对锈菌的基因组构造、基因组成、基因组变异等特征进行了归纳分析，对基因组变异与其专性寄生特性的关系、基因组变异对其毒性变异的可能影响等进行了阐述。1锈菌的全基因组测序分析锈菌为专性寄生、难以在人工培养基上培养、生活史具复杂多型现象和遗传转化困难等使得锈菌的基因组学及基因功能研究发现相对滞后。在第一个植物病原真菌水稻大角间座壳稻瘟菌Magnaportheoryzae基因组草图Deanetal.2005公布6年后，基于Sanger法测序的两个锈菌禾柄锈菌小麦专化型和落叶松杨栅锈菌全基因组测序分析发表。随后，基于二代高通量测序技术的5个锈菌
5、基因组草图先后发表表1。在锈菌基因组测序中，用于DNA提取的菌体都是相对易于收获的夏孢子一种双核的厚壁孢子。基因组DNA序列存在高度的杂合性，杂合度是人类基因组的1.6倍Zhengetal.2013，给锈菌基因组拼接带来严重障碍，作为基因组拼接完好性重要参数的CotigN50和ScalffoldN50均大大低于水稻大角间座壳等其他真菌基因组草图。这在基于二代测序技术的锈菌基因组拼接中更为明显，条形柄锈菌小麦专化型Pucciniastriiformisf.sp.tritici130基因组采用全基因组鸟枪法测序拼接，ScalffoldN50甚至只到达5.2kb表1。条形柄锈菌小麦专化型Puccin
6、iastriiformisf.sp.triticiCY32基因组通过采用FosmidtoFosmid二次拼接策略将这一参数提高到125kbZhengetal.2013。在已公布的锈菌基因组中普遍包含相对较大数量的预测蛋白编码基因，基因数量在14445个25288个之间表2。对两个小麦锈菌Pst和Pgt和落叶松杨栅锈菌Mlp编码蛋白基因的比照分析表明，两种小麦锈菌的同源基因有5560个，而Pst和Mlp的同源基因只要1937个，相比于Mlp、Pst和Pgt的同源关系更近。Pst、Pgt和Mlp种特异基因数分别为4469、5827和4922个。在Pst特异性的4469个基因中，只要14.43%10
7、66和3.38%250在GO和PAMGO数据库中有已知功能的同源基因Zhengetal.2013。这表明在基因组水平上，锈菌物种之间的差异特别显著，这或许与其生活史的复杂性及转主侵染的适应性需求相关。2锈菌基因组大小与锈菌演化基因组承载着生物进化的印迹，基因组大小是物种的重要特征。基因组的大小和动态变化对物种的进化合适度、繁衍和产生多样性的非有性机制都有直接的影响Dhondtetal.2011。真核生物基因组的大小变化范围很大，真菌是真核生物中基因组相对较小的类群。利用流式细胞技术对1850种真菌基因组大小进行评估，发现所检测真菌平均基因组大小仅为44.2Mbp，担子菌门真菌平均基因组大小为7
8、0.4Mbp，而锈菌平均基因组大小为305.5Mbp。近期发表的鬼针草单孢锈菌基因组到达了2489Mbp，是目前已知最大的真菌基因组Ramosetal.2015a。锈菌和其他真菌相比明显表现出基因组增大的进化趋向，这种扩增趋势可能与其寄生专化性和繁衍复杂性无性、有性及稀少有性的进化需求有关。具有大的、可塑性强的基因组可能有利于锈菌在分支选择cladeselection中避免灭绝，在与寄主共进化中以更高的频率产生新的毒性类型偶尔地扩展到新的寄主群体或新寄主。不同锈菌基因组大小差异能够到达30倍，从栎柱锈菌梭形专化型Cronartiumquercuumf.sp.fusiforme基因组大小76.6
9、Mbp到鬼针草单孢锈菌U.bidentis基因组大小为2489Mbp。锈菌属内种间的基因组大小变动范围也明显大于其他真菌，讲明基因组大小的变异可能在锈菌进化、遗传分化和寄主专化性演化中发挥活跃的作用Ramosetal.2015b。从已发表的植物病原真菌基因组比照结果能够看出专性活体营养型真菌如锈菌、白粉菌等基因组明显大于非专性活体营养、半活体营养和死体营养型真菌图1。尽管目前尚不能完好解释这一现象的成因，但通过增大基因组转座子作用为主要动力能够增加遗传多样性，进而减少锈菌作为活体营养的专性寄生菌在群体遗传漂移drift中灭绝的可能性。3锈菌基因组的转座子和重复序列组分锈菌类真菌基因组大小的显著
10、膨胀主要是其转座子和重复序列成分的增加造成的，其基因数量并没有同比例显著增加。转座子成分的活性可能是锈菌群体遗传和寄生专化性毒性保持多样性的重要动力，这对未发现转主寄主的不完全锈菌无性繁衍为主可能更为重要，它们也往往保持着较大的基因组Ramosetal.2015b。落叶松杨栅锈菌和禾柄锈菌小麦专化型全基因组序列分析结果表明，其基因组中不存在大片段的重复，但是含有大量的转座子和重复序列，约占整个基因组的45%。比拟不同营养型真菌基因组转座子和重复序列含量发现，与非专性活体营养真菌相比，专性活体营养锈菌类基因组转座子和重复序列含量明显更高图2。在不同锈菌基因组中，不同类型的转座子含量存在明显差异图
11、3。在落叶松杨栅锈菌基因组中，类末端反向重复序列TIR，14%转座子的比例较类长末端重复序列反转录LTR，8%转座子的比例更高。与之相反，在禾柄锈菌小麦专化型基因组中类TIR转座子和类LTR转座子所占基因组的比例则分别为9.8%和12.4%。在条形柄锈菌小麦专化型基因组中类长末端重复序列反转录转座子LTR和类末端反向重复序列TIR转座子含量最高，分别占基因组的27%和17%Zhengetal.2013。咖啡驼孢锈菌Hemileiavastatrix的类长末端重复序列反转录转座子LTR最多，占总基因组大小的38.7%。同时在咖啡驼孢锈菌基因组中发现了四种新的类长末端重复序列反转录转座子LTR。由
12、于基因组中大量转座子的存在造成基因组的重组，使得近源种间直系同源家族基因的共线性发生了很大的变化。目前转座子含量的差异对不同锈菌基因组以及致病性的影响还有待研究。转座子是基因组的基本组成部分，是促使基因组防御系统发展的主要推动力，并且能够多种方式对基因组的可塑性和演变作奉献。4锈菌基因家族扩增与丢失基因注释及比拟基因组分析发现，在进化经过中锈菌基因家族出现大量的扩增或缺失，以及家族内基因数量的增加或减少。相比于拥有4583个基因家族的原始祖先，禾柄锈菌小麦专化型Pgt目前有5413个基因家族，有1241个基因家族的增加和411个基因家族的减少。而落叶松杨栅锈菌Mlp分别增加了909个基因家族，
13、减少了188个基因家族，目前拥有5304个基因家族表2。一些膨胀的基因家族有物种特异性，可能在进化经过中这些基因对寄主专化性起重要作用。在禾柄锈菌小麦专化型Pgt和落叶松杨栅锈菌Mlp中，共有5045个同源基因。但是，由于扩张的转座子对基因组大量的改组，使得同源基因的共线性很低。Duplessisetal2011利用silico对Mlp和Pgt基因组分泌蛋白SP进行基因预测和手动注释，分别鉴定出了1184和1106个SP，其中谱系特异性基因分别占74%和84%。Mlp和Pgt基因组SP基因家族的膨胀非常显著，占整个扩张基因组家族的10%表2。对接种叶片后表达量最高的前100个基因分析发现：落叶
14、松杨栅锈菌Mlp和禾柄锈菌小麦专化型Pgt中分泌蛋白的数量分别占到50%和28%。在侵染的植物组织中大多数上调表达的SP基因是谱系特异性的，只要16%在两种锈菌中有同源基因。Zhengetal.2013通过比拟3种锈菌Pst、Pgt、Mlp基因组成分发现，条形柄锈菌小麦专化型特异的基因中多数已知功能的基因类型属于转运蛋白，激酶、碳水化合物激活酶类和分泌蛋白SPs，其中分泌蛋白所占比例最高。表明这些分泌蛋白基因进化速率很高，保守性很低，可能在毒性变异中发挥着重要作用图4。Nemrietal2014通过对4种锈菌Pst、Pgt、Mlp、Mli基因组候选分泌蛋白家族预测发现，编码几丁质酶的家族总共有
15、61个成员，其中32个是分泌蛋白。这表明一些效应蛋白可能从一些非分泌原始祖先进化而来，进而在真菌细胞壁以外的寄主细胞壁或者在寄主细胞内部执行类似的生物学功能。对已知的亚麻栅锈菌无毒基因、蚕豆锈菌效应蛋白、以及玉米黑粉菌Ustilagomaydis效应蛋白分支酸变位酶1进行生物信息学分析，在Pst、Pgt、Mlp中也未发现同源基因，或者同源基因不具有无毒基因特性Persoonsetal.2014。这表明分泌蛋白基因家族具有较高的进化速率以及较低的保守性。子囊菌纲的黄枝孢菌效应蛋白能够分泌到胞外被寄主免疫受体所识别DeWitetal.2009。然而，锈菌候选的胞外分泌蛋白能否也介入寄主免疫响应有待
16、进一步研究。Petreetal2016通过将Mlp候选的效应蛋白注射本氏烟的方法研究其亚细胞定位，结果显示六个效应蛋白定位于细胞核、核仁、叶绿体、线粒体和细胞外。其中一个效应蛋白叶绿体靶蛋白1能够模拟寄主转运肽进入叶绿体和线粒体，表明这些效应蛋白可能操纵线粒体和叶绿体来完成其专性寄生。5锈菌的基因丢失与生活方式演化锈菌是活体营养的专性寄生菌，不能人工培养。揣测锈菌在进化的经过中可能丢失了部分腐生性营养担子菌类所特有的某些基因。比拟分析21种已测序不同营养型真菌基因组中碳水化合物激活酶类基因成分发现，Mlp和Pgt具有相对较少的糖苷水解酶编码基因，数量和担子菌纲的双色蜡蘑较为类似，但是相比半活体
17、营养型、死体营养型和腐生型生物则要少很多。在进化为活体专性营养型的经过中，锈菌类真菌已经失去了一些编码分泌糖苷水解酶和多聚糖裂解酶的基因。但是，编码一些去除植物纤维素和半纤维素的糖苷水解酶的基因家族，例如GH7、GH10、GH12、GH26和GH27均出现了稍微的扩增。这些酶可能在锈菌侵染寄主的初始阶段起主要作用，例如蛋白酶、酯酶和一些碳水化合物激活酶类基因在被侵染的植物体内表达量高度上调，表明这些侵入菌丝能够通过水解酶类穿透进入寄主细胞ElGueddarietal.2002。锈菌通过侵入的菌丝在寄主体内构成吸器来获取糖类和氨基酸等养分对于成功建立起活体专性营养型互作具有重要意义。Mlp、Pg
18、t和Mli缺失了一些介入硝酸盐同化吸收的基因，例如硝酸盐和亚硝酸盐转运蛋白以及亚硝酸盐复原酶。Mlp具有基本的硫素同化吸收相关基因，而Pgt、Pst则没有。Pgt缺少了亚硫酸盐复原酶SiR的和亚基，而Mlp的SiR的亚基缺失，导致失去转酮醇酶活性，但是在亚麻栅锈菌Mli中存在完好的硫素代谢途径。锈菌缺失部分氮素和硫素同化吸收相关途径与他们专性活体寄生方式一致，它们能够通过扩张的转运蛋白家族直接从寄主植物细胞获取氮素NH4+或氨基酸和硫素Spanu，2012图4。锈菌可能通过丢失一些非必需基因来平衡其大进化macroevolution压力。6锈菌基因组变异与毒性演化在锈菌中，还未发现像水稻大角间
19、座壳那样由于大片段缺失造成分泌蛋白缺失而引起的毒性变异。甚至在比拟毒性差异最大的不同生理小种分泌蛋白组时，也没有在任何一个小种中发现超过15个预测的分泌蛋白基因缺失。这表明基因的缺失可能不是造成锈菌毒性改变的主要原因，而不同生理小种的分泌蛋白等位基因的变异可能是造成毒性变异的重要原因。通过重测序4个不同毒性的条形柄锈菌小麦专化型生理小种，生物信息学分析鉴定不同小种候选效应蛋白基因的单核苷酸多态性SNP变异，在2999个分泌蛋白中发现5个候选多态性分泌蛋白基因。这些多态性分泌蛋白可能反响不同生理小种条形柄锈菌小麦专化型对特定基因型寄主小麦的快速适应性。Zhengetal2013以PstCY32参
20、照基因组，对5个世界范围的重要流行菌系CY23、PKCDRD、Hu092、104E137A、Pst78进行重测序分析，检测到超过100000个SNPs。同时还发现不同菌系还存在丰富SV、InDel等遗传变异。其中，分布于编码区的cSNPs占总共SNP数量的35%，这些SNPs可能影响锈菌小种毒性和寄生特性的演化。以此全基因组SNP数据构建的6个菌系的进化关系显示，Pst分离系之间有明显的遗传分化，有性生殖可能在其毒性演化经过中发挥了重要作用。通过对15个落叶松杨栅锈菌不同毒性分离系进行了重测序，结果显示与Mlp98AG31基因组相比，不同菌系存在大量的单核苷酸变异SNVs、多核苷酸变异MNVs
21、；连续的SNVs和短序列的插入/缺失变异InDels。同时，落叶松杨栅锈菌16%的基因上游1000bp区段有超过10个SNPs，表明这些区段可能与基因表达调控有关。重复序列诱导的点突变是真菌所特有的一种真核生物基因沉默机制，1987年在粗糙脉孢菌Neurosporacrassa中初次发现。Hornsetal2012通过对来自担子菌门3个不同亚门伞菌亚门、黑粉菌亚门、柄锈菌亚门的8种担子菌转座子的变异分析发如今柄锈菌亚门的落叶松杨栅锈菌Mlp和禾柄锈菌小麦专化型Pgt中存在大量的RIP情况。在Mlp和Pgt中RIP主要发生在Gypsy类转座子的TpCpG三核苷酸区域，使TpCpG转换成为TpTp
22、G，进而构成富含T+A的片段。而在伞菌亚门和黑粉菌亚门中均未发现此类现象。重复序列诱导点突变RIP在减数分裂经过中通过突变多拷贝DNA，能最大限度地减少转座子的影响，这一特性被以为是为转座子传播提供的一种防御机制。此外，发生RIP的序列多集中在着丝粒区域，主要是转座子甲基化后的痕迹。锈菌类真菌基因组宏大、杂合度高，菌系地理分布广泛、生态位多样，蕴含着诸多基因组进化待解之谜。目前，高通量测序的低成本化和第三代测序技术的成熟化为锈菌构造基因组学提供了史无前例的机遇，通过基因组精细测序、群体重测序以及比拟基因组学分析能够更深化地揭示锈菌寄生专化性进化和毒性变异的遗传基础。对基因组变异的深化研究将为最终在分子水平揭示其毒性变异机制提供理论根据，为预测新的生理小种的产生及锈病防治提供有力工具。