牦牛源产细菌素屎肠球菌SWUN5732基因组测序及抑菌相关基因分析

【目的】本研究旨在通过研究1株牦牛源产细菌素屎肠球菌SWUN5732的基因组序列,发掘该菌种的关键功能性特征基因。【方法】通过Illumina PE150测序系统,对屎肠球菌SWUN5732进行基因组测序,并在GO、KEGG、COG和NR数据库对基因组的基本功能特征进行注释,结合NR数据库的注释,挖掘出基因组中潜在的抗菌、抗氧化和抑菌肽相关基因。【结果】屎肠球菌SWUN5732的基因组大小为2 824 168 bp,有效平均GC含量约为38.09%;可检测到的编码基因数量约为2 919个,全部编码基因总长度约为2 387 787 bp,平均长度约为818 bp。在GO数据库中SWUN5732基因组一共有8 851个基因得到注释,分别注释到分子功能、细胞组分和生物过程三大类45个条目;在KEGG数据库中SWUN5732基因组共有1 534个基因被注释,分别为细胞过程、环境信息处理、遗传信息处理、人类疾病、新陈代谢和生物体系统6大功能38个通路;在COG数据库中,SWUN5732基因组共有2 118个基因被注释,在基因组中还包含4个有ABC型氨基酸转运系统(渗透酶组分)的功能序列和具有重复...     

桑品种桂桑优12的转录组生物信息学分析

为了从分子生物学的角度解析桑树的重要功能性状,以广西蚕区主栽桑品种桂桑优12的根部组织为材料提取总RNA进行转录组测序.测序共获得50844314条原始测序数据(Raw data),经过滤后得到50540436条Clean data.对序列进行拼接组装、去冗余后,获得了102254个转录本,总长度为64665080 bp,平均长度为771 bp,N50值为1473 bp,GC含量47.2％.将获得的序列与各大功能数据库比对,进行功能注释,共有86332个转录本获得了注释结果,预测出68218个编码序列(Coding DNA sequence,CDS),总长度34282488 bp,平均每个CDS长度为502 bp,N50值为756 bp,GC含量42.44％.其中有40254个转录本在KEGG获得注释,有26832个转录本在COG获得注释,有27766个转录本在GO获得注释.桂桑优12根部组织的转录本数据分析结果,可为今后研究发掘桑树重要功能性状基因提供一定的数据支持.     

产细菌素乳酸菌LS-1的全基因组测序分析

前期研究从内蒙古手工奶酪中筛选出1株可拮抗金黄色葡萄球菌的乳酸菌Lacticaseibacillus saniviri(LS-1)。为探究LS-1的全基因组信息，挖掘其潜在细菌素编码基因簇，本试验对LS-1进行全基因组测序，对编码基因进行预测，并采用非冗余蛋白质(Nr)和基因本体(GO)数据库进行功能注释；利用BAGEL4数据库预测分析细菌素合成基因簇。结果显示，LS-1的基因组大小为2 356 714 bp, GC含量为47.87%。共预测到2 316个编码基因，其编码基因总长度为2 058 534 bp,平均长度为888 bp。LS-1所分泌细菌素属于Ⅲ类细菌素中的Enterolysin A(EnlA),推测EnlA可能通过裂解细胞壁来实现抑菌作用。本试验通过全基因组测序对LS-1基因组进行全面分析，从基因水平阐明其细菌素编码基因簇，进一步明确了菌株LS-1所产细菌素具有作为新型替抗产品的潜力。     

1株大鲵源多重耐药维氏气单胞菌全基因组解析及其致病性

为明晰重庆地区1株大鲵源多重耐药维氏气单胞菌的生物特征，本试验对分离菌中CQ-AV1株进行了全基因组测序以及致病性分析。利用二代测序技术中Illumina HiSeq测序平台对分离菌进行全基因测序，测序数据经修剪、过滤、去除接头和低质量碱基后，以HGA V1.0软件进行组装并上传至GenBank,并与数据库中已有序列构建系统进化树；利用VFDB(virulence factors of pathogenic bacteria)和CGE(center for genomic epidemiology)数据库对其毒力基因和耐药基因进行注释；将菌株分别腹腔注射健康小鼠和大鲵，观察其对2种实验动物的致病情况，并对发病大鲵的肝脏、脾脏、肺脏和肾脏组织进行病理切片，观察其病理组织学变化。结果显示，CQ-AV1株基因组长度为4 782 590 bp, GC含量为58.48%,共编码4 323个基因，另有假基因62个，GenBank中的注册号为RZII00000000.1;系统进化树显示该菌与猪源维氏气单胞菌菌株具有最近的亲缘关系，与鲈鱼源和罗非鱼源菌株亲缘关系次之，与人源菌株亲缘关系最远；致病性试验...     

西藏亚东牦牛线粒体DNA分析及系统进化研究

应用基因测序技术对西藏亚东牦牛线粒体基因组进行测序,并对测序数据质控与修剪及基因组组装、精细注释、功能注释与分析。结果表明:亚东牦牛线粒体基因组大小约为15 389 bp,包含13个PCGs、21个tRNAs、2个rRNAs及1个D-loop,其A+T碱基含量为60.86%,G+C碱基含量为39.14%;编码基因(PCGs、tRNAsrRNAs)片段总长度为14 495 bp,占基因组总长度的94.19%;同时,分析了亚东牦牛与其他11个物种间的进化关系,发现亚东牦牛与甘南牦牛归为一类,亲缘关系最近。西藏亚东牦牛线粒体DNA分析为研究牦牛群体遗传、起源与分子进化提供了技术参考。     

1株利用乳酸的猪源丁酸梭菌的分离鉴定及基因组序列分析

【目的】 研究1株利用乳酸的猪源丁酸梭菌的基因组信息。【方法】 采用厌氧培养法, 通过乳酸分解培养基(LADM)初筛、梭菌增殖培养基(RCM)复筛, 从饲喂乳酸菌发酵饲料的健康仔猪肠道内容物中分离利用乳酸快、丁酸转化率高的菌株, 对所筛选的菌株进行形态鉴定、16S rDNA序列分析、乳酸转化特性试验、基因组重测序和框架图测序, 并对其编码蛋白进行功能注释。【结果】 分离得到1株分离菌LY33, 基于形态和16S rDNA序列分析鉴定为丁酸梭菌。菌株LY33对乳酸具有较好的利用能力, 在生长第6天可将66 mmol的乳酸基本转化完全, 生成17 mmol左右的丁酸。LY33菌株的重测序表明, 其编码基因整体变异较小, 与参考基因组相比, LY33菌株全基因组杂合比率为0.022‰, 单核苷酸多态性(SNP)位点变异总数21 590个(占整个基因组0.4666%), 插入缺失(InDel)突变总和594个(占0.0128%), 不存在拷贝数变异(CNV), 结构变异(SV)注释的变异总数103个(占0.0003%)。突变基因主要为与菌株生长、能量代谢调节、维生素合成(特别是生物素)有关的基因, 分布于2条染色体的多条序列上。LY33菌株框架图测序及其编码基因蛋白的功能注释表明, 其基因组序列长度4 627 127 bp, GC含量28.60%, 含有4 171个基因, 编码区总长度占比84.12%, 参与了糖代谢(carbohydrate metabolism)、膜转运(membrane transport)和氨基酸代谢(amino acid metabolism)等多种代谢途径转化过程, 基因组中抗大环内酯类药物、抗杆菌肽、VanI糖肽抗性基因个数较多。【结论】 本研究从饲喂4%乳酸菌发酵饲料的健康仔猪肠道内容物中分离到1株利用乳酸快、丁酸转化率高的LY33菌株, 经鉴定为丁酸梭菌, 其具有良好的应用前景, 可作为一种新的微生物资源用于微生态制剂产品。     

副猪嗜血杆菌血清4型全基因组序列草图测定与分析

本研究应用Illumina Hiseq 2000测序系统对副猪嗜血杆菌(HPS)血清4型gx033株进行全基因组序列测定,经过DNA文库构建、测序、数据过滤和组装,绘制了全基因组草图.全基因组草图显示,HPS血清4型的基因组大小为2 155 493 bp,G+C含量为39.79％,含有编码基因2 189个,基因总长度1 881 801 bp.基因功能注释表明,849个基因功能尚不明确,1 340个基因被分为能量产生和转化、转录和细胞周期调控等19个功能组.该HPS血清4型全基因组草图的绘制为HPS病的致病机制等研究提供依据.     

1株大鲵源葡萄牙柠檬酸杆菌的全基因组测定及毒力基因和耐药基因的分析

葡萄牙柠檬酸杆菌是一种新发现的条件性致病菌,为明确重庆地区大鲵源多重耐药葡萄牙柠檬酸杆菌的基因组及其耐药特征,本试验对患病大鲵分离菌进行鉴定、药敏试验、全基因组测序及毒力、耐药基因的分析。分离菌首先通过16S rDNA进行初步鉴定,腹腔注射健康小鼠观察其致病性;利用Kirby-Bauer(K-B)纸片法进行9类共35种抗生素的敏感性试验;再利用二代测序技术进行全基因组测序,通过生物软件和数据库对分离菌的毒力基因和耐药基因进行分析。结果显示,分离菌鉴定为葡萄牙柠檬酸杆菌,将其命名为CQ-CP1,该菌对小鼠具有较强的致病性。药敏试验显示CQ-CP1对多种抗生素不敏感,为多重耐药菌株。全基因组测序该菌株共有5 006 960 bp, G+C含量为51.8%,共编码有4 722个结构蛋白,基因组信息在GenBank中的登录号为RZIH00000000.1。病原毒力因子数据库(VFDB)比对发现CQ-CP1中共有49个与毒力相关的基因,与细菌的运动、黏附、转录等相关;耐药基因数据库比对结果显示,该菌携带6个耐药基因(ARG),分别是氨基糖苷类aadA1,喹诺酮类qnrS1、qnrB32,β内酰胺...     

1株西藏牦牛巴氏杆菌全基因组测序与比较基因组学分析

为进一步对西藏牦牛巴氏杆菌病基因组学研究提供理论基础,从西藏林芝某地病死牦牛肺脏中分离鉴定出1株荚膜A型牦牛巴氏杆菌,并通过SMRT法对该菌株进行全基因组序列测定,对测序数据组装后基因组组分及比较基因组学进行分析。结果显示,获得大小为2.3 Mb基因组,GC含量40.30%。同时预测2 096个编码基因数量,编码基因序列总长度2 047 410 bp。预测出总长度为4 739 bp的重复序列,重复序列含量0.21%。预测非编码rRNA数量为19、tRNA数量为59。假基因数量为3,假基因总长度为736 bp。此外,对测序的1株菌株和2株参考菌株(Mannheimia haeemolytica、Pasteurella multocica)进行基因组共线性分析,显示测序菌株与参考菌株之间共线性良好。对测序的1株菌株和2株参考菌株进行家族分类,共得到2 105个基因族,其中,3株菌共有的基因族类(核心基因族)有1 483个,占总基因组的70.45%。进化分析研究显示,测序菌株与参考菌株Mannheimia haeemolytica较远。表明本研究对西藏牦牛巴氏杆菌分离菌株进行全基因测序,同时进行基因组组分及比较基因组学分析,为牦牛巴氏杆菌病的进一步研究提供数据支持。     

利用全基因组重测序数据检测8个鸭品种基因组拷贝数变异

旨在鉴别与鸭重要经济性状潜在相关的拷贝数变异(copy number variations, CNVs),为解析CNVs对鸭经济性状的影响提供前期研究基础。本研究利用从美国国家生物技术信息中心(National Center for Biotechnology Information, NCBI)公共数据库中下载的8个鸭品种共78个个体的全基因组重测序数据，采用CNVnator和CNVcaller软件进行全基因组CNVs检测，同时只保留两个软件检测结果中存在至少1 bp重叠的同类型CNVRs,以消除假阳性对试验结果的影响。结果显示，8个鸭品种的CNVs合并后共检测出7 550个CNV regions (CNVRs),其中包括7 098个duplications和452个deletions。这些CNVRs在鸭29条常染色体上呈不均匀分布，总长度为16 111.2 kb,平均长度为2 134 bp,约占鸭基因组的1.51%。此外，本研究在8个鸭品种共筛选到4 304个潜在的品种特异性CNVRs,覆盖1 230个注释基因。通过基因功能Gene Ontology(GO)富集分析，鉴别到38个可...     

1株鸡滑液囊支原体的分离鉴定与全基因组分析

鸡滑液囊支原体（Mycoplasma synoviae，MS）能引起家禽呼吸道疾病、传染性滑膜炎、蛋壳顶端异常及生长迟缓等，给全球养禽业带来巨大的经济损失。然而不同国家和地区MS的耐药性和致病性存在差异，因此本研究旨在挖掘福建MS菌株的特征。作者采集福建省疑似患MS的病鸡关节腔内容物，对菌株进行分离鉴定，并利用三代联合二代测序技术对其进行全基因组测序。结果显示，该分离株病原特性稳定，已保藏至中国典型培养物保藏中心，保藏编号为CCTCC No.：M 2021210，并将其命名为MS-FJ01。MS-FJ01基因组全长为795 381 bp，GC含量为28.39%；预测编码基因数为704个，占整个基因组的90.48%；并在COG、GO和KEGG数据库分别注释到476、382、425个基因。通过各数据库的分析，注释了40个耐药相关基因，4个碳水化合物相关酶基因，136个病原与宿主互作相关基因，47个毒力因子相关基因，106个膜转运蛋白相关基因，25个限制修饰相关基因。通过与MS HN01和MS 5-9基因组进一步比较发现，MS-FJ01与这两个中国菌株均存在良好的共线性关系，且与MS 5-9的亲缘关系更近。本研究阐明了福建MS菌株的基因组结构、通过相关数据库的分析，预测和注释了其基因的功能，为深入研究MS提供参考依据。     

布鲁菌疫苗株S2全基因组测序及功能分析

为找出布鲁菌疫苗株S2基因组分泌性蛋白、参与生物学途径基因、细胞组件基因、分子功能基因、毒力相关因子及耐药基因,对布鲁菌疫苗株S2全基因组测序并进行Ⅲ型分泌系统效应蛋白预测、GO注释、VFDB注释、ARDB注释。全基因组测序结果表明,S2基因组大小约为3.3 Mb,杂合度和重复度较低,GC含量分布未出现异常。通过Ⅲ型分泌系统效应蛋白预测未找出S2分泌性蛋白。通过GO数据库注释,发现参与生物学途径基因有3 823个、细胞组件基因有1 949个、分子功能基因有2 585个。通过VFDB数据库注释,找出S2基因组有79个毒力相关基因,其中最大基因为8 604bp,最小基因为180bp。通过ARDB数据库注释,找出S2基因组有13个耐药基因,最大基因为3 225bp,最小基因为333bp;其中功能注释的有GL002835基因,是枯草杆菌肽类基因,此基因耐受抗微生物类药物。本研究为进一步了解布鲁菌S2基因组基本功能奠定了基础。     

纳米孔测序技术对一份禽流感病毒阳性样品的分子溯源

H5亚型禽流感为人兽共患病,可严重威胁人类和动物健康。直接对样品进行禽流感病毒靶向全基因组测序、系统进化分析,对于揭示禽流感病毒分子特征具有重要作用。本研究采用第三代靶向纳米孔测序技术,对1份禽流感病毒阳性拭子样品进行了全基因组测序。序列鉴定结果显示,病毒为H5N5亚型禽流感病毒,序列覆盖度达96.44%,仅PB1基因约480 bp的序列未被测出;HA蛋白裂解位点含有5个连续的碱性氨基酸,具备高致病性禽流感病毒分子特征,未发现向SAα2-6 Gal受体基因转变的氨基酸突变;HA基因进化分析结果显示,病毒为2.3.4谱系,NA基因为欧亚谱系。本研究采用靶向纳米孔测序技术,成功对灭活禽泄殖腔拭子中的禽流感病毒进行了基因组测序,证实该技术快速、便捷,可直接用于拭子样品的病毒全基因组测序,适用于动物流感病毒的快速分子鉴定与溯源。     

1株跨种裂解李氏杆菌噬菌体的分离鉴定

本研究旨在获得天然李氏杆菌噬菌体,提供防治李氏杆菌病新型生物制剂,减少抗微生物药物的使用,抑制病原菌耐药性的产生。利用产单核细胞李氏杆菌作为宿主菌对屠宰场污水进行双层琼脂平板法筛选,获得噬菌体,并对其进行透射电镜观察、生长特性检测(温度、pH、一步生长曲线、最佳感染复数、有机溶剂影响)、基因组酶切鉴定和全基因组测序分析。分离出的产单核细胞李氏杆菌噬菌体中,选取1株裂解性强,遗传稳定的噬菌体进行后续试验,并命名为LP8;经电镜观察为肌尾科噬菌体,可以跨种裂解18株产单核细胞李氏杆菌和5株威尔斯李氏杆菌,确定LP8的最佳感染复数为1、最适生长温度为45 ℃、最适pH为7;在6种有机溶剂中,仅异戊醇可导致LP8活性丧失;对提取的基因组进行酶切鉴定,确定为双链DNA;LP8全基因组测序结果表明,基因组大小为87 038 bp、含有120个编码基因、编码基因的累计长度为76 326 bp、编码基因的平均长度为636 bp、编码区域长度占基因组的比例为87.69%。本试验分离出产单核细胞李氏杆菌噬菌体,并对其噬菌能力和应用价值进行鉴定。分离出的LP8噬菌体相较于李氏杆菌的噬菌体,细菌裂解能力更强,适应环境范围更广。本研究为实验室后续建立产单核细胞李氏杆菌噬菌体库和产单核细胞李氏杆菌噬菌体的其他应用提供了良好的基础。     

上海白猪(上系)基因组遗传变异检测与功能注释分析

拟通过全基因组范围内遗传变异的检测和注释分析,深入了解上海白猪(上系)群体的遗传现状,以便更好地进行提纯复壮与开发利用。本研究应用基因组简化测序平台—GGRS,对上海白猪(上系)群体(99头)进行基因组测序,并综合实验室前期太湖地方猪种和西方引进品种共计447头测序数据进行SNP和InDel等遗传变异检测和功能注释分析。结果显示,上海白猪基因组测序平均覆盖度为2.87%,平均测序深度为3.9倍,共检测到328 586个SNPs、693 220个InDels。进一步分析表明,SNP和Indel在染色体上的分布相对均一,但在不同染色体上数量和密度的分布存在一定差异。结构注释显示,SNP和InDel对应基因的数量分别为11 496个和13 216个,按照基因的结构区域分类,SNP和InDel在各类功能基因区间的分布特点一致,绝大多数的SNP和InDel变异都分布在基因间区,分别为74.61%和83.38%。内含子中次之,分别为22.76%和14.98%。基因富集分析结果表明,SNP主要与蛋白代谢过程和高分子代谢过程等相关,InDel多在组织、器官发育和疾病相关的通路中发生富集。此外还发现,含有高密度SNP和InDel等遗传变异的QTL主要集中在肉质与胴体性状和健康性状。本研究利用先进的简化基因组测序技术对上海白猪(上系)全基因组范围内的遗传变异进行了普查,并通过对太湖流域地方猪种和西方引进品种的合并比较分析,阐释了这些遗传变异的染色体分布特征,基因区间分布规律和功能注释信息,为后续的开发利用奠定了坚实的基础。     

鸭源马尿气球菌的鉴定及全基因组测序分析

为鉴定一株鸭源分离菌,并探究其生物信息学特性,本研究分别通过Pacbio Sequel三代测序技术和Illumina NovaSeq二代测序技术对该株鸭源分离菌进行全基因组测序;基于该分离菌16S rRNA基因序列构建进化树;通过GTDB基于基因组之间的平均核苷酸相似性及采用infer软件构建120个单拷贝基因序列进化树并对其分类;利用NR数据库比对分离菌全基因组蛋白编码基因序列;采用VFDB数据库分析该分离菌相关毒力基因;采用K-B纸片法测定分离菌对21种药物的药敏性试验;采用CARD数据库分析分离菌的耐药基因;利用GO、KEGG、COG数据库对分离菌进行基因功能、信号通路及基因的注释。测序结果显示分离菌由一个环形染色体和一个环形质粒组成,大小分别为2 133 811 bp和37 505 bp,GC含量分别为39.04%和35.03%;通过GTDB分类及NR数据库比对分析显示,该菌为马尿气球菌亚种;VFDB数据库比对发现6种致病毒力基因;药敏试验结果显示分离菌对红霉素、多西环素、庆大霉素、卡那霉素等多种药物耐药,CARD数据库比对发现6种耐药基因;基因功能分析显示,1 277个基因参...     

不同生长时期梅花鹿鹿茸转录组分析

本研究旨在探究不同生长时期梅花鹿鹿茸的转录组差异,丰富梅花鹿鹿茸转录组信息。选取5个生长时期(10、20、28、44、66d)的梅花鹿鹿茸组织作为试验样品,利用Illumina HiSeqTM2500高通量测序,并用测序评估、基因注释等生物信息学方法进行分析。结果,经过测序、转录本拼接获得了375 657条contigs,平均长度为688bp,329 946个unigenes,平均长度为534bp。通过比对Nr、Nt、Pfam、KOG/COG、Swiss-prot、KEGG、GO等7大数据库,90.07%unigenes得到了成功注释。其中39 674个(12.02%)基因成功注释到GO数据库,在level 2水平上共分为41个类别;17 732个(5.37%)基因成功注释到将KOG的26个类别中。对5个生长时期的鹿茸转录本文库进行比较分析,共发现了509个差异基因,其中407个差异表达基因成功注释到GO数据库中,主要为信号转导、氧化还原、转录调节、蛋白酶降解等生物学过程。其中转录调节基因PER1、EGR1的表达随着鹿茸的生长而增加,而GAS1则相反,随着鹿茸的生长而降低,推测PER1、EGR1、GAS1等转录因子在鹿茸生长过程中可能起着重要的调节作用。本研究利用高通量测序技术对不同生长时期的梅花鹿鹿茸组织进行了转录组测序和分析,筛选到了梅花鹿鹿茸在不同生长时期下的差异表达基因,并获得了差异表达基因的功能。     

牙釉质基因鉴定山羊性别的研究

根据牙釉质(amelogenin,AMEL)同源基因在山羊X、Y染色体上存在不同程度碱基缺失的特点,设计1对引物,对山羊血液基因组DNA混合样进行PCR扩增,将PCR产物纯化、克隆并测序;然后对46个山羊血液DNA样本(♀:22个,♂:24个)进行性别鉴定。结果显示,雌性山羊样品经扩增只产生一条非特异性条带,雄性山羊样品产生一条非特异性条带和一条特异性条带;非特异性片段长度为349 bp,特异性片段为289 bp;46个山羊血液DNA样本鉴定结果与实际性别对比,雌性准确率为100%(22/22),雄性准确率为100%(24/24)。     

长白山野杂猪和东北民猪背最长肌的mRNA转录组比较分析

本研究通过RNA-seq技术对3个长白山野杂猪与3个东北民猪背最长肌样本进行了mRNA测序和差异分析,从而为筛选长白山野杂猪肉质性状功能基因奠定基础。结果显示:在长白山野杂猪与东北民猪中均有100%以上的序列能完全比对到猪的基因组上,而所转录的基因占总基因的比例在3头长白山野杂猪中平均为59.67%,3头东北民猪为58.73%,而在不同个体间SNP位点数为14 299～19 671;在2个品种猪背最长肌中共表达基因有15 842个,仅在长白山野杂猪中表达的基因有619个,仅在东北民猪中表达的基因有710个;在2个品种猪背最长肌中共获得显著差异表达的基因153个,其中在长白山野杂猪中上调表达的基因有55个,下调表达的基因有98个;GO分析显示,153个差异表达基因被注释到33条GO条目中,其中11条被注释到生物学过程中,11条被注释到细胞组分,11条被注释到分子功能;KEGG分析显示,这些差异表达基因被富集到130条通路中,其中主要与脂肪代谢相关的信号通路有关。结果表明,长白山野杂猪与东北民猪背最长肌中mRNA表达存在差异,这些差异表达基因将是将来研究2个品种猪间肉质差异的主要靶标之一。     

牦牛皱胃全长转录组测序分析

【目的】获得牦牛(Bos grunniens)皱胃全长转录组数据库,深入挖掘牦牛皱胃功能基因。【方法】采用PacBio Sequel高通量测序系统,使用单分子实时(single molecular real time, SMRT)测序技术对成年牦牛皱胃全长转录组进行测序,对原始数据进行质控和去冗余分析,再比对参考基因组获取过滤后的非重复序列(Unigenes),使用多种生物信息软件对牦牛皱胃全长转录本数据进行功能注释、转录因子注释、编码区预测、简单重复序列(simple sequence repeats, SSR)分析及可变剪接分析。【结果】通过测序共获得14 467 420条子序列,平均子序列长度为3 344.23 bp,质控得到循环一致性序列(CCS)有296 840条,全长非嵌合(FLNC)序列有277 402条,过滤和去冗余后对比参考基因组最终获得8 556条Unigenes。通过与NR、Swiss-Prot、KEGG、KOG、eggNOG、GO和Pfam数据库比对,对Unigenes进行注释,其中NR数据库注释了8 544条Unigenes; Swiss-Prot数据库注释了...