姜科基因组复杂性的比较基因组学研究.pdf

以菖蒲为外类群，对姜科物种进行比较基因组学分析，结果表明，姜、草果及砂仁在进化过程中经历并共享了至少4次全基因组加倍事件,且各物种基因组在分歧之后没有再发生加倍事件，各物种基因组受多倍化影响发生了大量的基因丢失。对姜科基因组复杂性的比较研究对其物种遗传基因功能创新具有重要意义，可为单子叶植物进化史研究提供参考。关键词：姜科；基因组复杂性；多倍化中图分类号：Q943.2Comparative Genomics Study on the Genome Complexity of Zingiberaceae(1.North China University of Science and Technology,Tangshan 063210,China;2.Genomics and ComputationalBiology Research Center,North China University of Science and Technology,Tangshan 063210,China)Abstract:The study does comparative genomic analysis of Zingiberaceae with calamus as an outgroup.The results showginger,grass and amomum have experienced and shared at least 4 genome-wide doubling events in the course ofevolution,and the doubling events havent occur again after the divergence of the genome of each species,and a largenumber of genes are lost due to the influence of polyploidy.The comparative study of genomic complexity ofZingiberaceae is of great significance for the innovation of genetic function of its species,and can provide reference forthe study of the evolutionary history of monocotyledonous plants.Key words:Zingiberaceae;Genome complexity;Polyploidization0引言姜科（Zingiberaceae）是单子叶植物姜目（Zingibe-rales)的一个重要分支,在生产生活中占据重要地位。其下属植物具有重要的经济价值及药用价值，如姜（Zi n g i b e r o f f i c i n a l e）、草果（Amomum tsao-ko）及砂仁（W u r f b a i n i a v i l l o s a），作为中药应用了几千年 。多倍化是物种进化与分歧的重要推动力 2 。研究表明,单子叶植物在进化过程中均经历了全基因组加倍（W G D，w h o l e g e n o me d u p l i c a t i o n），包括姜科在内的大部分单子叶植物共同拥有一次全基因加倍事件，即TWGD3-5。全基因组加倍导致的基因丢失、易位等现象对基因组结构及复杂性造成了重大影响 5-7 。研究表明，姜在进化过程中从单子叶植物祖先到现在共经历了3 次全基因组加倍事件,草果基因组在文献标志码：B文章编号：16 7 4-8 6 46(2 0 2 3)14-0 0 6 8-0 4Hu Fubo,Wang Xiyin?近期没有经历过全基因组加倍事件，砂仁基因组的研究表明,砂仁和姜的共同祖先可能经历了1次全基因组加倍事件。但现有研究对姜科植物古多倍化的认识还很模糊。随着单子叶植物基部物种的测序完成，为姜科植物在进化过程中经历的更古老加倍事件的研究提供了条件。WGD收稿日期:2 0 2 3 -0 5-10基金项目：国家自然科学基金项目“基于基因共线性的植物基因组学数据库建设研究”（3 2 0 7 0 6 6 9）作者简介：胡福博（19 9 6），男，硕士。研究方向：比较与进化基因学。通讯作者：王希胤（19 7 2）,男,博士，教授。研究方向：比较与进化基因组学。68Atz ZofAtswvi CnuAta代表菖蒲，Zof代表姜，Ats代表草果，Wvi代表砂仁,Cnu代表椰子图1待研究植物系统发育树Fig.1 Plant phylogenetic tree to be studied对菖蒲（Acorus tatarinowii）基因组的研究表明，它只经历过1次全基因组加倍事件 8 。椰子（Cocos nu-cifera)在进化过程中除与姜科祖先物种共享了古老的TWGD以外，只单独经历了1次全基因组加倍事件,两者基因组相对保守。4 这有助于进一步了解姜科植物基因组结构的复杂性。1材料与方法1.1物种基因组数据姜基因组数据下载自公共数据库Genebank（I n-dex of/genomes/genbank/plant/Zingiber_officinale/lat-est_assembly_versions（n i h.g o v)/G CA _0 18 446 3 8 5.1_Zo_vl.1/）。草果基因组数据来自国家基因库生命大数据平台 CNGBdb（h t t p s:/d b.c n g b.o r g/s e a r c h/p r o-ject/CNP0003772/）。砂仁基因组数据来自公共数据库 Refseq(https:/ftp.ncbi.nlm.nih.gov/genomes/ref-seq/plant/Elaeis_guineensis/all_assembly_versions/GCF_000442705.1_EG5/。菖蒲基因组数据来自中国国家基因库(https:/gb.org/pub/CNSA/data4/CNP0001708/CNS0456199/CNA0036157/）。编写python脚本,将下载得到的原始数据进行数据预处理得到所需的注释文件（gff）、染色体长度文件（lens）、蛋白序列文件(pep）及蛋白编码序列文件（cds）。1.2系统发育及共线性基因推断使用蛋白序列比对工具BLAST+9,对研究物种蛋白序列（pep)文件进行种内及种间同源基因搜索（E=value100)。使用 orthofinder 提取待研究物种的单拷贝基因,构建系发育物种树（见图1）。根据blast得到的结果，运行生信分析流程软件WGDI的“d模块,绘制基因组内及基因组间的同源基因点阵图10 ,运用共线性分析软件ColinearScan，提取所物种AtaCnuZofAtsWvi2.2Ks 分析同义核苷酸替换（ks）是蛋白质编码序列的核苷酸变异不引起氨基酸改变 13 。有研究利用椰子基因组内共线性基因计算了ks 分布并进行数字拟合，发现有两个明显的峰,认为ks 峰值大的峰代表椰子经历的WGD。提取了那部分共线性基因绘制了其ks柱状图,计算了姜、砂仁及菖蒲基因组内与组间共线性基因的ks分布并绘制ks柱状图（见图2）。发现菖蒲和姜研究物种基因组内及基因组之间的共线性基因对(共线性片段的基因对=5 个)。1.3Ks 计算为了区分不同全基因组加倍事件产生的共线基因,计算了待研究物种间同源性基因对的核苷酸同义替代(Ks）,以估计共线基因之间的分歧水平。利用clustalW将基因对的编码序列进行比对，调用PAML包的Nei-Gojobori 的方法进行Ks分布计算 12 。利用WGDI对ks分布进行数学拟合，得到ks峰值。1.4?染色体同源区域深度推断编写Python 脚本,将利用软件colinearscan提取的姜基因组与菖蒲同源的染色体片段投影到蒲染色体上，构建一个展示姜基因组染色体同源区域深度的列表,结合上述处理得到的染色体长度（lens）文件并运行软件WGDI的“-ci模块将列表可视化。2结果分析2.1基因共线性分析对共线性分析软件colinearscan获得的共线性片段进行分析，分别统计了菖蒲、姜、草果、砂仁的基因组内及基因组间的共线性基因对数量5个以上的共线性片段的共线性基因数（见表1）。姜基因组内共线性基因对有1118 5对,砂仁有6 8 15对,草果有 6 6 8 对。菖蒲与姜基因组间的共线性基因对有48 0 9 对,与砂仁有3323对，与草果有6 6 8 对。同理,椰子与姜基因组间的共线性基因对有1442 5对,与砂仁有10 18 9 对,与草果有9 9 3 5对。结果显示,姜基因组进化过程中的同源共线性基因对保留最多,砂仁次之,草果最少，故后续分析以姜作为主要研究对象，探究姜科基因组的复杂性。表1同源基因统计Tab.1 Statistics of homologous geneAtaCnu58511971323 92129 865480916 07366810535332310596Zof398514 42511 18520 78124 255基因组间的共线性基因对的ks峰值为1.43,对应菖蒲和姜发生分歧的时间。姜和砂仁基因组间共线性基因对的ks峰值为0.3 6,对应姜和砂仁发生分歧的时间。椰子共线性基因对的ks峰值为0.9 9,对应椰子经历-WGD的时间。姜基因组内的同源共线性基因对的ks分布柱状图显示在0.3 6 1.43 有很多小峰,推测在与蒲分歧之后，姜科各物种在发生分歧之前可能经历了包括-WGD在内至少4次全基因组加倍事件。69Ats448993520 41466814 206Wvi286910 18915 55511 5366815越多,说明姜经历的全基因组加倍事件次数较多,受古老的加倍事件影响产生的重复基因随物种进化而大量丢失。Ata_Zof Ata3Ata2Ata4Ata5Cnu_CnuAta6Ata7Ata8 Ata9Zof_WviAtalAta12AtallAtal0Ch#123456789101112Zof_Zof-0900.99图2 菖蒲、椰子、姜、砂仁共线性基因同义核苷酸替换Fig.2Collinear gene synonym nucleotide substitution ofcalamus,coconut,Zingiberaceae and amomum kernel2.3染色体同源区域深度推断在一定程度上,染色体同源区域深度可以大致反映物种进化过程中经历的全基因组加倍次数。将姜与菖蒲基因组同源性区域映射到菖蒲染色体并绘制圈图（见图3）。结果显示，姜映射到菖蒲染色体上的同源区域深度最高达到13，由此推测，姜在进化过程中可能至少受到4次全基因组二倍乘事件的影响,越靠近外圈,姜的同源结构空白占比越大，代表同源基因丢失a7891011121314151617181920212223240.51.0A311.5Ats62.0ks外圈代表菖蒲映射的姜同源性基因图3 菖蒲映射姜的染色体同源片段深度Fig.3 Chromosomal homologous fragment depth of acorusmapping Zingiberaceae2.4基因组间同源基因点阵图为探究姜科物种的基因组结构，绘制了菖蒲、椰子、姜等物种基因组间的同源基因点阵图。物种间的同源点阵图中的线性片段是由物种分化产生的同源基因片段。如姜与草果、砂仁的点图中（见图4）,草果的2号染色体同源最好匹配姜的1号染色体，砂仁的2号染色体同源最好匹配姜的1号染色体，说明姜科祖先基因组在分化形成现存物种之后并没有再经历全基因组加倍事件。b7Wvi9101112131415161718192021222324234JOZ6JoZ7R91011点图中红点表示同源性最好，蓝点表示同源性次好，灰点表示同源性较差。图4姜与草果、砂仁基因组间的同源点阵图Fig.4 Homologous dot map of the genome of Zingiberaceae,grass fruit and amomum kernel7091011菖蒲与姜的同源基因点图呈现的点非常散乱，不过部分菖蒲染色体对应的姜的同源片段数量显示大于8条（见图5）。椰子与姜的点图显示，椰子6 号染色体对应姜的同源片段数量达到15 2 2 条，再次说明姜基因组在进化过程中经历的多倍化事件至少有4次，甚至可能更多。至于椰子对应的姜的同源片段数量不同,可能是由于椰子受2 次多倍化事件的影响,造成染色体基因的缺失。a16b61416图5#姜与菖蒲、椰子基因组间局部同源基因点阵图Fig.5ILocal homologous gene dot map of Zingiberaceae,calamus and coconut genome3结论与讨论以菖蒲为外类群,对姜科下属植物姜基因组进行比较研究，了解姜科经历的全基因组加倍事件可能多于3 次,有可能是4次也可能更多,且姜科下属植物草果在进化过程中受多次加倍事件产生的重复基因大量丢失,导致其更为复杂的基因组结构,令人们对姜科基因组的深入挖掘难度增加。随着测序技术的不断更新换代，会有更多姜科乃至单子叶下属植物基因组完成更为精确的测序工作,有了更多精确的物种基因组数据做参考，会令姜科下属植物基因组复杂性结构的分析与挖掘得到长足进步,结果更为清晰、准确。生信分析软件及算法的开发，可为未来姜科基因组比较分析研究提供更强大的技术支撑。参考文献：1 Yang P,Zhao HY,Wei JS,et al.Chromosome-level genome assemblyand functional characterization of terpene synthases provide insights intothe volatile terpenoid biosynthesis of Wurfbainia villosa J.Plant J,2022,112(03):630-645.2 Wang J,Yu J,Sun P,et al.Paleo-polyploidization in lycophytes J.Genomics Proteomics Bioinformatics,2020,18(03):333-340.3 Cheng SP,Jia KH,Liu H,et al.Haplotype-resolved genome assemblyand allele-specific gene expression in cultivated ginger J.HorticZofRes,2021,8(01):188.23456789.1011口Zof12345564 Wang S,Xiao Y,Zhou ZW,et al.High-quality reference genomesequences of two coconut cultivars provide insights into evolution ofmonocot chromosomes and differentiation of fiber content and plantheightJ.Genome Biol,2021,22(01):304.89.10115 Paterson AH,Bowers JE,Chapman BA.Ancient polyploidizationpredating divergence of the cereals,and its consequences forcomparative genomics J.Proceedings of the National Academy ofSciences of the United States of America,2004,101(26):9903-9908.6 Wang X,Guo H,Wang J,et al.Comparative genomic de-convolution ofthe cotton genome revealed a decaploid ancestor and widespreadchromosomal fractionationJ.New Phytol,2016,209(03):1252-1263.7 Marques A,Klemme S,Houben A.Evolution of plant B chromosomeenriched sequencesJ.Genes(Basel),2018,9(10):115-118.8 Shi T,Huneau C,Zhang Y,et al.The slow-evolving Acorus tatarinowiigenome sheds light on ancestral monocot evolution J.Nat Plants,2022,8(07):764 777.9 Camacho C,Coulouris G,Avagyan V,et al.BLAST+:architecture andapplicationsJ.BMC Bioinformatics,2009(10):421.1o Sun P,Jiao B,Yang Y,et al.WGDI:a user-friendly toolkit forevolutionary analyses of whole-genome duplications and ancestralkaryotypesJ.BioRxiv,2021(05):337-340.11 Wang X,Shi X,Li Z,et al.Statistical inference of chromosomalhomology based on gene colinearity and applications to Arabidopsisand riceJ.BMC Bioinformatics,2006(07):447.12 Yang Z.PAML 4:phylogenetic analysis by maximum likelihoodJ.Mol Biol Evol,2007,24(08):1586-1591.13 Pan ZX,Xu D,Zhang JB,et al.Reviews in comparative genomicresearch based on orthologsJ.Hereditas,2009,31(05):457-463.14 Wang J,Sun P,Li Y,et al.An overlooked paleotetraploidization incucurbitaceaeJ.Mol Biol Evol,2018,35(01):16-26.(上接第6 7 页)参考文献：1】韩艳,董嘉雯，王龙龙，等.承插式装配桥墩抗震性能拟静力试验与数值模拟 J.工程抗震与加固改造，2 0 2 0,42（0 5）：8-12.2徐艳，曾增,葛继平，等.承插式预制拼装桥墩的最小合理承插深度J.同济大学学报（自然科学版）,2 0 19,47（12）：17 0 6-17 11.3 Guo A,Li H,Ba X.et al.Experimental investigation on the cyclicperformance of reinforced concrete piers with chloride-induced corrosionin marine environmentJ.Eng Struct,2015,105:1-11.4 Ou Y,Fan H,Nguyen N.Long-term seismic performance of reinforcedconcrete bridges under steel reinforcement corrosion due to chlorideattackJ.Earthq Eng Struct Dyn,2013(42):2113-2127.5】中华人民共和国住房和城乡建设部.水泥基灌浆材料应用技术规范M.北京：中国建筑工业出版社，2 0 15.6刘建鑫.混凝土结构设计规范（GB50010-2010）修订动态及主要内容分析 J.呼伦贝尔学院学报,2 0 11,19（0 1：4.7吴元，杨晓婧，王凯，等.灌浆料单轴受压应力应变曲线试验研究 J.工业建筑,2 0 14(S1）：9 0 9-9 13.71