运用高通量技术对杂交鳢进行转录本测序和数据分析,经拼接与组装,最终获得52 065条unigenes,长度范围201～12 407 bp,序列平均长度为710 bp。从长度分布、GC含量和基因表达量等方面对unigenes进行评估,数据显示测序质量较好。进一步的数据分析,共鉴定出20 535个CDS,37 324个SNP位点和7 830个微卫星。用6大数据库Swiss–Prot、Nr、Pfam、KEGG、KOG和GO注释杂交鳢转录组unigenes,分别对应有20 955、28 938、19 339、14 052、19 358和18 635条unigenes获得注释。根据GO数据库的注释,可将这些unigenes分为生物过程、细胞组分和分子功能3大类共50亚类。KEGG代谢通路分析表明,这些unigenes参与了5大类30亚类共268个代谢通路,其中,参与信号转导通路的unigenes数量最多,有1 716条。

【目的】为获得白芷转录组序列信息,了解不同产地间白芷的差异表达基因情况。【方法】以白芷根为研究对象,采用Illumina HiSeq X Ten平台对不同产地白芷进行高通量转录组测序,并对其差异表达基因进行筛选和分析。【结果】川白芷产地间、亳白芷产地间及川白芷与亳白芷间的差异基因总数分别为2686、1704、11 549条,分别被GO富集注释11 038、434、5604条,差异基因的主要功能与细胞组成、细胞、细胞器、结合、催化活性、细胞过程及代谢过程有关。KEGG通路分类显示,差异表达基因被分为4大类,24个小类,所在通路主要与代谢过程相关。与白芷主要化学成分香豆素和挥发油合成相关的4-香豆酸连接酶、丙氨酸氨解酶、3-羟基-3-甲基戊二酰辅酶A还原酶在川白芷中表达量升高,1-脱氧-D-木糖-5-磷酸合成酶、脱氢多萜醇焦磷酸合成酶在川白芷中表达量降低。【结论】通过转录组测序分析了不同产地白芷根基因差异表达状况,为白芷药材道地性品质的形成提供了理论依据。

为获得岩原鲤转录组信息，发掘功能基因，本研究采用ＩｌｌｕｍＩｎａ高通量测序技术对岩原鲤全组织转录组进行测序。结果获得６４ ２５７ ９１８个ＥＳＴ序列，经拼接和组装得到８３ ２５２条单基因序列（ｕｎＩｇＥｎＥ），平均长度７８７ ｂｐ，长度范围２０１～１６ ５７２ ｂｐ。利用ｎＣｂＩ的蛋白质非冗余数据库（ｎｒ）对所有ｕｎＩｇＥｎＥ进行相似性搜索，共有３７ １５７条ｕｎＩｇＥｎＥ（４４．６３％）与数据库中的已知序列同源。利用ｂｌａＳＴ２ｇＯ ｖ２．５软件对ｕｎＩｇＥｎＥ进行注释，共得到２９ ９１９条（３５．９３％）注释基因，根据ｇＯ功能分类将其分为生物过程、细胞组分和分子功能３大类５６亚类。经ＫＯｇ．．．

单细胞转录组测序(scRNA-seq)使研究人员能够在单细胞分辨率下研究基因的表达，了解细胞的动态发育轨迹.然而，单细胞分辨率的基因表达矩阵往往具有高维度和高稀疏性的特点，使后续数据处理面临巨大的挑战.因此，本文综述了scRNA-seq数据分析的一般流程和常用软件，这将有助于研究人员选择更合适的分析软件和流程.

分析梁山慈竹及其体细胞突变体的转录组，挖掘功能基因并对其差异表达基因进行筛选和分析，为梁山慈竹遗传改良提供理论依据。利用ＲＮＡ－Ｓｅｑ技术进行转录组测序，对测序结果进行ｄｅ Ｎｏｖｏ拼接和功能注释；并对差异表达基因进行筛选及ＣｏＧ、Ｇｏ、ＫｅＧＧ数据库中进行比对注释，此外，基于ＳｗｉＳＳ－ＰＲｏｔ功能注释结果，分析纤维素和木质素相关功能基因的表达量差异。测序结果表明，共获得８６ ５７５ ６３１条ＲｅＡｄＳ，ｄｅ Ｎｏｖｏ组装得到８４ ７４１条ｕＮｉＧｅＮｅＳ，共有４９ ８２９条被ＮＲ、ＣｏＧ、Ｇｏ、ＫｅＧＧ、ＳｗｉＳＳ－ＰＲｏｔ注释。从梁山慈竹实生植株（对照）和体细胞突变体Ｎｏ．３０这２个测序．．．

利用新一代高通量测序技术平台Illumina Solexa对泡桐进行转录组测序和数据denovo组装,并对得到的unigene进行功能注释、分类及代谢通路分析。结果表明,N50值为2278bp、平均长度为1410bp的泡桐unigene共有188019条,将其与NCBI的Nr和Swiss-Prot数据库比对后发现,分别有120808条和85880条unigene与其他物种的基因具有同源性。利用COG数据库可将有关的41914条泡桐unigene分成25类,KEGG数据库分析发现共有43553个unigene参与215种代谢通路。找到与木质素合成相关的泡桐unigene。

为研究饲喂钩吻对仔猪肾脏的影响,将10头健康雄性去势二元杂种(长白猪×大白猪)断奶仔猪,随机分为对照组和给药组,每组5只.对照组仔猪饲喂不含抗菌药物的全价日粮,给药组仔猪饲喂不含抗菌药物的全价日粮+2%钩吻全草粉末,连续饲喂45 d后,分别采集两组猪的肾组织,对其进行高通量测序,并对差异表达基因进行基因本体论(GO)功能分类及京都基因和基因组百科全书(KEGG)数据库分析.结果显示,两组肾脏中共有301个差异表达基因.与对照组相比,给药组有180个基因显著下调,121个基因显著上调,其中与代谢、转运相关的基因,如溶质载体家族2成员13(SLC2A13)、溶质载体家族35成员5(SLC35A5)、含黄素单氧化酶2(FMO2)、醛氧化酶1(AOX1)等出现显著下调(P<0.05).荧光定量PCR对差异表达基因进行验证,其结果与转录组学结果基本一致.表明饲喂钩吻对仔猪肾脏的物质转运和药物代谢有一定影响.

［目的］采用高通量测序技术ＩｌｌｕｍＩｎａ ｍＩ Ｓｅｑ２５０获得蕨类植物芒萁的孢子体转录组数据，以期为芒萁的生长、发育、代谢调控、微进化机制分析等提供重要的分子信息。［方法］应用生物信息学方法对测序获得的大量单基因簇（ｕｎＩｇｅｎｅ）进行基因功能注释、代谢途径及微卫星分析等。［结果］本研究共获得１８ ４６３ ２９６条序列读取片段（ｒｅａｄＳ），总碱基数为４．６２ｇｂｐ序列信息，经序列组装最终得到６３ １６９个ｕｎＩｇｅｎｅ，平均单条ｕｎＩｇｅｎｅ长度为８６３ ｂｐ，ｎ５０为１ ５８７ ｂｐ，其中分布在２００ ５００ ｂｐ长度区间ｕｎＩｇｅｎｅ占总数的５５．４％。数据库中的序列同源性比较表明，．．．

为研究绒山羊绒毛生长周期内皮肤基因表达规律以及皮肤差异表达基因的挖掘,采用Illumina HiSeqTM2000高通量转录组测序平台RNA-Seq技术,对3只成年雌性阿尔巴斯型内蒙古绒山羊全年12个月的皮肤样本进行转录组测序,将无参考基因组de novo拼接组装的unigene运用Nr,GO,COG以及KEGG数据库进行比对注释,并进行差异基因在全年各月份间的变化规律分析。结果显示,测序获得的unigene数为105 854条,比对注释基因55 541个,其中共有51 078条转录本进行了GO注释,21 189条转录本在COG数据库分析得到注释,26 201条预测基因KEGG数据库比对成功。其中皮肤基因GO注释在生物学功能分类中,细胞过程、代谢过程、生物学调控、生物功能调控以及刺激应答方面注释到的转录本比例最大;在分子功能分类中,结合以及催化活性作用方面注释到的转录本比例最大;通过COG数据库的比对注释,看到皮肤表达基因的同源蛋白功能主要集中在基因转录、翻译后修饰、和信号传导机制方面,反而在同工酶转运和代谢、细胞动力、二级代谢生物合成、转运和催化以及核结构方面注释到的基因极少;皮肤转录组数据与KEGG数据库比对,得到皮肤相关转录本较为富集的几条通路分别为MAPK signaling pathway,Wnt signaling pathway,Notch signaling pathway,Hedgehog signaling pathway,TGF-βsignaling pathway,JAK-STAT signaling pathway以及Focal signaling pathway;通过绒山羊不同月份皮肤转录组测序数据进行差异基因分析,得出绒山羊全年皮肤生长周期中基因差异变化规律为绒山羊全年皮肤共有4次较为剧烈的基因差异变化,第1次发生在2月与3月之间,第2次发生在3月与4月之间,第3次发生在6月与7月之间,第4次发生在10月与11月之间,前2次出现的基因差异变化较后2次剧烈得多,说明绒山羊绒毛生长启动时基因变化剧烈,而绒毛生长到休止时的皮肤基因变化较为缓和,是一个基因缓慢变化的过程。

【目的】揭示天麻种植土壤真菌群落组成结构，天麻种植前后土壤真菌群落结构的变化及种完天麻种方竹后土壤真菌群落结构的变化，为研究天麻根际土壤微生态提供理论依据并为解决天麻连作障碍提供技术支持。【方法】基于高通量测序方法对采集的3组土样（原壤charlle 001、一麻壤charlle 002、方竹壤charlle 006）进行基因组DNA提取，PCR扩增，产物检测、纯化和定量，建库测序，数据分析，物种组成多样性分析及差异分析。【结果】3组土样中真菌主要类群为Ascomycota、Basidiomycota和Mortierellomycota三门真菌，其中Ascomycota真菌丰度普遍最大；3组样本特有真菌OTU：一麻壤>方竹壤>原壤；种水平下分析原壤charlle 001中主要优势真菌为Clitocybe sp.、Solicoccozyma sp.和Spirosphaera sp.，一麻壤charlle 002中主要优势真菌为Clavulinopsis sp.和Rosellinia sp.，方竹壤charlle 006中主要优势真菌为Trichoderma sp.；Clavulinopsis和Rosellinia在一麻壤charllle002中丰度极大在charllle001和charllle006中丰度极小；charlle001和charlle002组间真菌Mortierella相对丰度差异显著，charlle001和charlle006组间真菌Penicillium相对丰度差异显著；charlle002和charlle006组间物种差异Penicillium相对丰度差异显著。种植收获一茬天麻后土壤中真菌种类数量递增，群落结构发生变化，一麻壤中真菌ASV、特有真菌数目明显增大；一麻壤种植方竹后木霉菌显著增多；Clavulinopsis和Rosellinia部分真菌易引起天麻病害，而方壤中Trichoderma多用于病原菌生物防治。【结论】天麻不同种植土壤的真菌群落结构不同，种过一年天麻后土壤中病原菌增多，再种植方竹呈现修复作用，但具体修复机制还有待后续深入研究。

为全面了解虾夷扇贝(Patinopecten yessoensis)微卫星分布频率和数量,深化对虾夷扇贝基因组的认识,本研究运用第二代高通量测序技术,进行虾夷扇贝简化基因组测序(RAD-seq),从基因组水平阐明虾夷扇贝基因组微卫星特征。结果显示,简化基因组测序共获得序列总长为92,551,435 bp,经过滤筛选,获得259,535个contig,其中,包含微卫星序列3618条,经引物设计共获得3460对微卫星引物。统计微卫星序列的重复类型,其中,三核苷酸重复单元数量最多(1587个,45.87%),其

[目的]研究马尾松抗松材线虫病相关基因,为解析马尾松抗性机理和高抗马尾松早期选择提供理论依据。[方法]对高抗和易感2种马尾松基因型接种松材线虫,在接种后第1、15和30天取样进行转录组高通量测序,通过对高抗和易感马尾松差异基因识别以及富集分析,以筛选与抗松材线虫病相关的组成型差异基因。[结果]对高抗和易感马尾松接种松材线虫后的转录组进行比较,在接种后第1、15和30天分别获得2 866、679和1 657个差异基因,且在接种松材线虫后不同时间点间,共有差异基因相对较少。对差异基因进行GO富集分析,结果显示:在接种后第1、15和30天,最显著的生物学进程是氧化-还原过程,而GO项刺激反应、转录调节以及香叶基二磷酸代谢过程也显著富集。对这几类GO项相关基因进一步分析,发现接种后3个时间点GO项刺激反应中共包括26个R基因,除了2个R基因外,其它R基因表达皆为组成型,在接种松材线虫和对照间差异不显著,而GO项转录调节中仅2个差异基因被nr注释为ERF转录因子,其余为未知基因。接种松材线虫后,ERF转录因子在高抗马尾松上的表达量始终高于易感马尾松上的表达量,且在接种松材线虫和对照间表达量变化不显著,也为组成型基因表达。GO项香叶基二磷酸代谢过程中,GGPPS 3个基因也在高抗马尾松上具有更高表达量,为组成型基因表达。这些基因中TIR-NBS-LRR基因(c65785.graph_c0)、ERF转录因子(c78073.graph_c0)和3个GGPPS基因在高抗马尾松中表达量较高,而在易感马尾松中表达量较低,甚至为0,可作为开发分子标记的候选基因用于鉴定抗性马尾松。[结论]高抗和易感马尾松在接种松材线虫后第1天其表达量差异达到显著水平的基因最多,其中,LRR基因、ERF转录因子和GGPPS与马尾松的抗性相关,部分TIR-NBS-LRR基因、ERF转录因子和GGPPS有望被开发为鉴定高抗马尾松的分子标记。

动物基因组学研究是进行动物遗传资源保护和利用以及分子育种的一项重要基础工作，高通量测序技术的出现为动物基因组学研究带来了革命性飞跃。文章对基于高通量测序技术的动物基因组从头测序、重测序、简化基因组测序等基因组学研究进行了综述，总结了相关的生物信息分析内容，并对不同分析工具及方法进行了比较分析。最后，讨论了当前高通量测序技术存在的问题，对该技术未来发展方向进行了展望。

为深入了解柔鱼食性，探明其在食物网中的位置和生态功能价值，本研究根据2022年6—8月于西北太平洋公海采集的15尾柔鱼样本，采用高通量测序技术分析其食物组成。结果显示，15尾样品的高质量序列为646980条，通过聚类分析共得到624个运算分类单元（OTU），经筛选比对检测到的饵料类别共计65种，分属6门11纲22目37科51属。研究表明，柔鱼主要摄食头足类、鱼类以及多种浮游动物（包括磷虾类、介形类、水母类等）。优势饵料包括拟沙丁鱼、日本爪乌贼、北方拟黵乌贼、萤乌贼、尾棘背灯鱼等。柔鱼摄食范围广，饵料物种与海区的天然饵料有关；不同胴长组柔鱼表现出不同的摄食偏好。本研究通过柔鱼胃含物的全组成分析，对深度挖掘摄食与柔鱼生长、洄游、繁殖等生活史的关系有重要作用，同时为基于生态系统角度的柔鱼资源可持续开发和养护提供科学依据。

【目的】建立一种批量分析T-DNA插入位点的简单、有效的方法。【方法】提供一种基于高通量测序技术的T-DNA插入位点的分析方法，将其命名为PSORA:Parallel sequencing of one round amplicons。首先对一轮交错式热不对称PCR(TAIL-PCR)的扩增产物进行高通量测序，随后通过生物信息学分析T-DNA插入位点，该方法降低了TAIL-PCR过程中对特异性扩增的要求。PSORA中使用的引物一侧为简并引物，另一侧为T-DNA特异性引物。在特异性引物的5’端设计6 nt的样品标签（Barcode），用于标记不同的转化事件。所使用的5个转基因烟草株系（L1、L6、L9、L15和L19）由农杆菌介导质粒pBI121转化获得。此外，通过标准PCR对PSORA的结果进行验证。【结果】利用PSORA对5个转基因株系的T-DNA插入位点进行分析，结果显示，L6含2个插入位点（NW＿015801367的36 316 bp处和NW＿015950898的42 202 bp处），L9、L15和L19各含1个插入位点（L9的插入位点为NW＿015943682的235 969 bp处；L15的插入位点为NW＿015802951的60 529 bp处；L19的插入位点为NW＿015863435的12 188 bp处），L1的插入位点未能成功获取。对生物信息学分析结果进行PCR验证，含不同插入位点的转基因株系之间可互为阴性对照，野生型（WT）作为空白对照，结果表明，在各转基因株系中均得到与预期一致的特异性扩增，该结果验证了PSORA的有效性。【结论】PSORA是一种分析T-DNA插入位点的有效方法，可以同时分析多个插入事件，相对于传统的染色体步移方法更简便、快速。