生防菌贝莱斯芽孢杆菌MC2-1全基因组测序分析

【目的】分析贝莱斯芽孢杆菌MC2-1的全基因组序列信息,挖掘该菌株拮抗物质合成的相关基因,为其生防机理研究和开发应用提供数据基础。【方法】采用二代Illumina Hiseq与三代Pacbio平台相结合的测序技术,对从美洲大蠊肠道内分离获得的烟草疫霉拮抗菌MC2-1进行全基因组测序,并对测序数据进行基因组装、预测、功能注释、共线性分析及次级代谢产物合成基因簇预测。【结果】菌株MC2-1的全基因组大小为3929792 bp,平均GC含量为46.5%,共编码4015个基因;包含tRNA基因86个、rRNA基因27个、sRNA基因81个;含有基因组岛2个、前噬菌体1个、CRISPRCas 9个。在NR、Swiss-Prot、COG、GO和KEGG数据库中分别注释到4015、3347、2996、2841和2223个基因。在CAZyme数据库中注释到几丁质酶、纤维素酶和淀粉酶等可降解植物病原真菌细胞壁的相关酶基因。预测到菌株MC2-1中有12个次级代谢产物合成基因簇,编码表面活性素（Surfactin）、丰原素（Fengycin）、杆菌素（Bacillibactin）、杆菌溶素（Bacilysi...     

嗜麦芽窄食单胞菌CA1定殖及促生长特性相关基因挖掘

【目的】挖掘嗜麦芽窄食单胞菌CA1定殖和促生长特性相关基因,为探究其促生长作用机制提供理论参考。【方法】利用Nanopore测序技术平台对CA1菌株进行全基因组测序,将其基因组序列与GO、NR、KEGG、COG等数据库进行比对,获得基因功能注释信息,挖掘分析其功能特异性基因;使用基质培法收集甘蔗根系分泌物,并分析根系分泌物对菌株CA1的趋化性、生物膜形成和定殖相关基因表达的影响。【结果】嗜麦芽窄食单胞菌CA1基因组大小为4371327 bp,携带6个质粒,GC含量64.58%,其中4696个基因被注释预测,包括4525个蛋白编码基因、87个t RNA、13个rRNA和1个tmRNA。GO、COG、KEGG和NR数据库分别成功注释2795、2078、2494和2584个基因。该菌基因组含有细菌趋化性（cheA、cheV、cheY）、生物膜形成（flhA、flk、flgG）、细菌黏附（pgaA、mrkD）、嗜铁（fes、fepA、fecA）、合成氨（gltP、gltB、gltD）和磷酸盐代谢（phoR、phoB、pstA）等促进植物生长和定殖相关的基因。趋化性试验结果显示,与对照组相比,根...     

牛分枝杆菌CVCC68002全基因组测序与功能分析

【目的】为进一步深入研究牛分枝杆菌CVCC68002基因组特征,致病性和遗传进化关系,为牛分支杆菌感染疾病的预防和深入研究提供参考。【方法】利用Illumina和Pacbio测序技术平台对CVCC68002菌株进行全基因组测序建库,采用生物信息学和系统进化树等方法,对该菌株基因组功能注释、致病性和遗传进化关系进行了分析。【结果】通过分析发现该菌株基因组全长4 357 100 bp,GC含量为65.64%,编码基因达4 145个;在GO、KEGG、COG、NR数据库基因组特征显示其具有高氨基酸和碳水化合物代谢能力。通过PHI,VFDB,ARDC和CARD数据库对其致病性进行相关分析,得到456个毒力基因注释,主要涉及细菌黏附与定殖、细菌逃避宿主免疫防御系统、分泌系统蛋白和转录因子等。并找到相关包括万古霉素,青霉素,多胺类抗生素、异烟肼、吡嗪酰胺、氟喹诺酮等的抗药性基因,另外,外排泵的抗性蛋白注释量显著,推测牛分枝杆菌通过外排抗性蛋白从而达到抗药自我保护的目的。【结论】获得牛分枝杆菌CVCC68002株完整基因组信息,完善系统进化信息,明确了其毒力因子和有效地靶标抗药基因,为后续对该菌株分...     

基于全基因组预测和分析花生果腐病原菌分泌蛋白

【目的】为分析花生果腐病病原菌-新孢镰刀菌的侵染机制并筛选可用基因。【方法】以全基因组测序结果为基础,根据分泌蛋白与效应因子的特点,采用SignalP、TMHMM、WoLF PSORT及NetGPI等软件一步步筛选,对新孢镰刀菌全基因组12 756条蛋白序列进行了分泌蛋白预测和功能分析。【结果】新孢镰刀菌全基因组编码蛋白中有701条序列具有典型的分泌蛋白特征,占蛋白总数的5.49%,其长度主要分布在101～700个氨基酸范围内。信号肽的分析结果表明,信号肽-3和-1位置的氨基酸相对保守,其切割位点属于典型的A-X-A型。对分泌蛋白功能预测结果表明,656个分泌蛋白获得了功能注释,其功能主要涉及碳水化合物的运输、代谢过程,蛋白翻译后修饰和氨基酸代谢、运输过程。分泌蛋白中效应蛋白预测结果表明,新孢镰刀菌分泌蛋白中共有161个潜在的效应蛋白。【结论】对新孢镰刀菌分泌蛋白的有效预测和功能分析,为筛选新效应蛋白,明确其致病机理,筛选寄主抗性基因奠定了基础。同时为花生果腐病抗性品种选育和综合防治提供理论依据。     

2 个烟草青枯病菌菌株的致病力差异与基因组比较分析

【目的】探索比较 2 个不同青枯病菌菌株之间致病力差异和基因组差异,分析影响菌株致病力差异的关键基因,为烟草青枯病的防治提供理论参考。【方法】从烟草青枯病发病烟株上分离获得不同的青枯病菌菌株并进行鉴定,通过侵染试验比较 2 个菌株之间的致病力差异,并对 2 个菌株的基因组进行测序,比较二者之间基因组序列差异,分析影响菌株致病力差异的基因。【结果】从江西省抚州市广昌县和南丰县分离到2 株青枯病菌菌株 RsTP2 和 RsTY2,egl 测序表明菌株 RsTP2 和 RsTY2 均属于亚洲分支演化型Ⅰ中的序列变种13,其中菌株 RsTP2 对云烟 87 致病力较强,而菌株 RsTY2 致病力较弱。通过对 2 个菌株进行基因组测序,获得二者的基因组框架图,其基因组大小分别为 5.43 MB 和 5.23 MB。将获得的基因和蛋白序列与相关数据库进行比对,获得基因的功能及物种注释信息。通过分析发现,菌株 RsTY2 有 52 个蛋白和菌株 RsTP2 的 44 个蛋白比对后序列存在差异,其相似度在 50.00%~99.88%,其中存在较大差异的编码蛋白主要是凝集素类蛋白,其次是细胞溶素分泌激活蛋白和穿孔素家族蛋白。【结论】通过比较不同致病力青枯病菌株 RsTP2 和 RsTY2 的基因组序列差异,发现二者基因组中编码凝集素类蛋白、细胞溶素分泌激活蛋白和穿孔素家族蛋白等蛋白序列差异较大,这类编码蛋白主要参与细菌的黏附、膜溶解、穿透等功能,与细菌侵入宿主细胞的能力关系密切。这几类编码蛋白的差异推测是导致菌株 RsTP2 和 RsTY2 致病力差异的重要因素之一。     

新暗色柱节孢H1全基因组测序和注释

【目的】测序并组装火龙果溃疡病病原菌新暗色柱节孢(Neoscytalidium dimidiatum)H1的基因组序列，为了解其基因组构成、遗传变异、致病机理和制定防控策略奠定基础。【方法】提取贵州火龙果溃疡病病原菌N.dimidiatum H1的基因组DNA,利用Illumina NovaSeq和Nanopore PromethION测序平台进行基因组测序、组装和注释分析。【结果】通过测序和组装获得N.dimidiatum H1的基因组规模为43.78 Mb, Gap数量为0;组装获得14个contig, N50长度为3.92 Mb,平均长度为3.13 Mb,最长长度为5.63 Mb。使用真核和真菌直系同源基因集进行BUSCO评估，完整的单拷贝基因比例分别为98.60%和98.50%。N.dimidiatum H1的基因组含有12 962个编码基因和205个非编码RNA。N.dimidiatum H1的基因组中重复序列含量为1.62%,简单重复为主要类型。转座子含量为1.59%,以LTR类Gypsy亚族和LINE类为主要类型。共计12 703个(98.00%)编码基因被注释，与同属于...     

山麻鸭开产期和产蛋高峰期卵巢组织转录组分析

【目的】探明山麻鸭开产期和产蛋高峰期卵巢组织的转录组差异,丰富蛋鸭转录组数据信息。【方法】选取6份山麻鸭卵巢组织样品（开产期和产蛋高峰期各3份）,利用Illumina HiSeq^TM 2000进行高通量测序,构建山麻鸭开产期和产蛋高峰期卵巢组织转录组文库,并用测序评估、基因功能注释等生物信息学方法进行分析。【结果】经过测序获得质量不低于20的碱基比例（Q20）均高于93%;开产期获得了43 489 798条reads,4 380 847 061 bp数据量;高峰期获得了42 782 676条reads,4 307 499 083 bp数据量;分别有70.92%和72.70%的reads能比对到鸭参考基因组序列上。对开产期与产蛋高峰期转录组进行比较,发现开产期有26 808个表达基因,产蛋高峰期有27 013个表达基因,与开产期为参考,在高峰期卵巢组织中共获得1 929个差异表达基因,其中上调基因989个,下调基因940个;进一步将这些差异基因在Nr数据库进行注释,发现获得注释的基因有1 423个,其中上调基因695个,下调基因728个;COG功能注释分析发现这些差异表达基因共获得663个功能注释,涉及25个功能分类;GO 功能注释分析表明,有1 122个基因获得GO 功能注释,可以分为61个功能分类,这些分类主要涉及到分子绑定、催化活性、细胞过程、生物调节等诸多生理生化过程,其中涉及到发育繁殖生物学过程的相关注释基因有406个;KEGG分析发现共有425个基因注释到160个代谢通路中,其中GnRH信号通路、GPI-anchor生物合成、TGF-β信号通路、核糖体生物合成等通路显著富集。GnRH信号通路和TGF-β信号通路在卵巢的生理生殖活动发挥了重要作用,而GPI-anchor生物合成和核糖体生物合成,这两类代谢通路在卵巢生理生殖活动的作用尚不明确。【结论】利用高通量测序技术对山麻鸭开产期和产蛋高峰期卵巢组织的转录组进行测序和分析,揭示了山麻鸭不同生理状态下卵巢组织差异表达基因的数量,获得了差异表达基因的功能、分类和代谢通路。为丰富蛋鸭卵巢组织转录组信息,为开展蛋鸭产蛋性状相关基因的研究及分子调控机制研究奠定基础。     

苹果果实大小相关的ARF-Aux/IAA互作组合筛选

【目的】通过转录组学和生物信息学,在苹果全基因组中对可能互作的MdARF和MdIAA进行鉴定,为明确相关基因功能和解析生长素调控苹果果实大小的分子机理奠定基础。【方法】对野生大果型‘皇家嘎啦’和miRNA172p过表达的转基因小果型‘皇家嘎啦’进行不同发育时期和不同组织材料的转录组测序,对测序结果进行基因的功能注释及差异表达分析。利用转基因小果和野生型大果的转录组数据,筛选在果实发育中表达的苹果MdARFs和MdIAAs基因家族成员,通过逐一在两个家族间计算基因时空表达的相关系数,筛选可能互作的MdARF-MdIAA组合,将从拟南芥基因组中下载的23个ARFs和34个Aux/IAAs、番茄基因组中下载的21个ARFs和25个Aux/IAAs,分别与互作候选MdARFs和MdIAAs进行比对,并进一步构建系统发育树。使用MEME和TBtools对苹果候选互作对中的MdARFs和MdIAAs蛋白进行Motif分析。利用STRING蛋白互作预测数据库进行同源映射,构建苹果中的蛋白-蛋白互作网络,进一步的确认候选互作对,最终得到苹果中通过互作参与果实发育可能性最高的MdARF-MdIAA组合。...     

林麝肺源ST1249型铜绿假单胞菌的分离鉴定及其全基因组序列分析

【目的】对分离自林麝化脓肺脏的1株疑似铜绿假单胞菌(PA)进行鉴定,并分析其致病和耐药机制,为林麝PA感染化脓性疾病的防控奠定基础。【方法】将病原菌分离纯化后,依次进行生化试验、16S rRNA序列分析、药敏试验和小鼠致病性试验,并通过全基因组测序,对分离菌株进行群体进化和物种分型分析以及基因功能注释。【结果】分离自林麝化脓肺脏的1株疑似铜绿假单胞菌经鉴定与PA相符,命名为FMDP001。药敏试验显示其对阿莫西林、头孢曲松、氨曲南、多粘菌素B和林可霉素耐药;对小鼠半数致死量为9.4×10~5 CFU/mL。全基因组序列分析显示该菌基因组大小为6 955 100 bp,序列类型为ST1249,与B136-33株同源性最高,且两菌株基因组平均核苷酸一致性(ANI)值达98.93%;全基因组中共有357个序列编码FMDP001与致病性相关的基因,根据功能分为黏附、铁摄取、胞外毒性蛋白和调控系统;84个序列编码耐药基因,其中多药耐药外排泵为主要成分。【结论】从林麝化脓肺脏中分离到1株致病性较强的PA,并获得ST1249型林麝源PA的全基因组序列,序列显示该菌携带大量药物外排泵及生物膜形成相关基因,决定其具有多重耐药特性,哌拉西林等可作为该类型PA感染的临床用药。     

青花菜花蕾转录组测序分析及蜡粉合成相关基因挖掘

【目的】对青花菜花蕾进行转录组测序分析，并挖掘与蜡粉合成相关基因，为探明青花菜花球表面蜡粉形成的分子机制提供理论参考。【方法】分别提取野生型和蜡粉缺失型青花菜花球总RNA，采用Illumina HiSeqTM 2500平台进行转录组测序，获得高质量Clean reads，采用Trinity进行序列组装后获得青花菜Unigene库，将获得的Unigene序列与Nr、Nt、KEGG、Pfam、KOG/COG、Swiss-Prot和GO数据库比对，获得基因功能注释信息；使用DESeq2进行差异表达分析。【结果】共获得44.68 Gb Clean data，De novo组装得到41244条Unigenes，N50长度为1847 bp。从所获得的Unigenes中筛选出8685个差异表达基因（DEGs）（上调基因5747个，下调基因2938个），共有8038个基因被注释到不同数据库，其中，5220个基因注释到Pfam数据库； 2066个基因注释到COG数据库，3866个基因注释到KOG数据库； 2580个差异表达基因被注释到75个转录因子家族中，注释最多的是MYB家族（235个）；GO数据库中6095个差异表达基因注释到细胞组分、分子功能和生物学过程三大类的52个功能分类； KEGG数据库中，1671个差异表达基因富集到138条代谢通路，其中13个差异表达基因与脂肪酸合成有关，7个差异表达基因与蜡粉生物合成途径有关。【结论】转录因子MYB家族在调控青花菜蜡粉合成中发挥重要作用。蜡粉合成过程中相关酶基因的差异表达是调控青花菜蜡粉合成的关键，尤其是野生型和蜡粉缺失突变体中特异性表达的差异表达基因，可作为后续研究青花菜花球表面蜡粉形成分子机制的对象。     

山羊5个组织的全长转录组构建与分析

【目的】构建山羊（Capra hircus）全长转录组,旨在改进参考基因组的功能注释。【方法】利用PacBio长读段测序技术对南江黄羊肝、脾、肾周脂肪、背最长肌和睾丸5个组织的混合样本进行高通量测序;对原始数据进行校正、聚类去冗余,最后用KEGG数据库对预测转录本序列进行功能注释。【结果】质控后获得28 432条高质量一致性全长序列,其中99.49%的序列被成功比对到山羊基因组;得到源自9 879个基因的19 115个转录本。进一步分析得到表达基因总数的41.39%发生了可变剪接;其中起始外显子可变剪接占比最高（36.47%）,外显子互斥占比最低（1.78%）。另外,144条序列未被比对到参考基因组上,这些序列的平均GC含量高达56.50%;利用KEGG数据库对其进行功能注释,揭示28个转录本与免疫系统相关。【结论】获得的测序数据全面涵盖混合样本中表达的转录本,证明了山羊基因组转录的复杂性。研究结果总体上提高了每个基因平均表达的转录本数量,证明可变剪接的广泛存在;未成功比对的转录本序列GC含量较高且许多转录本来自免疫基因的表达。     

基于银耳转录组测序的多糖代谢途径分析

【目的】对银耳进行转录组测序分析,分析多糖生物合成途径和挖掘关键功能基因。【方法】采用Illumina Hiseq 2500测序平台对银耳菌丝体进行测序,组装得到unigene序列,并与公共数据库进行对比分析和注释,挖掘果糖和甘露糖代谢通路相关基因。【结果】本研究共得到4.72 Gb raw data,拼接后获得17 008条unigene序列,N50为2073 bp;注释结果表明,76.47%序列能够注释到Uniprot公共数据库,注释到GO分类和KEGG通路的unigene分别为10 152和5671条。进一步分析银耳果糖和甘露糖代谢途径相关基因,挖掘得到74条unigene;其中编码L-iditol 2-dehydrogenase的unigene最多,其次是butanol dehydrogenase。【结论】基于转录组测序数据解析及预测,为研究银耳多糖和其它产物的生物合成途径及分子机制提供了数据支持,也为银耳品质的形成提供理论依据。     

辣椒疫霉拮抗菌Y4-39全基因组测序与分析

【目的】解析辣椒疫霉拮抗菌Y4-39的全基因组序列信息,阐明其防病促生机制,为其开发和应用提供理论依据。【方法】采用三代牛津纳米孔测序（ONT）和二代测序平台（BGISEQ）相结合的方法,对从黑水虻肠道分离获得的对辣椒疫霉具有较强抑制活性的贝莱斯芽孢杆菌（Bacillus velezensis） Y4-39进行全基因组测序、组装、基因预测和功能注释,并进行比较基因组学分析,从全基因组水平挖掘菌株Y4-39的次生代谢产物合成基因簇。【结果】菌株Y4-39全基因组大小为3891759 bp,GC含量46.67%,编码3713个蛋白基因,27个rRNA和86个tRNA基因,不含质粒。基于全基因组序列构建的B.velezensis系统发育进化树可分为5个亚群,亚群间的平均核苷酸一致率（ANI）大于97%,亚群内各菌株之间ANI大于98%。预测得到12个潜在的次生代谢产物合成基因簇,包括表面活性素、大环内酰亚胺H、bacillaene、芬枯草菌素、地非西丁、杆菌素和杆菌溶素等抑菌活性物质,同时还存在5个产物未知的次生代谢产物合成基因簇。【结论】贝莱斯芽孢杆菌种内存在一定程度的分化。菌株Y4-39基因组含有多个次生代谢产物合成基因簇,具有合成多种抑菌活性物质的能力,具有较好的应用前景。     

淅川乌骨鸡全基因组转座子的鉴定与分析

【目的】通过分析淅川乌骨鸡全基因组重测序数据,获取淅川乌骨鸡转座子的基本信息和特征分布,探究转座子相关基因参与的通路。不仅对研究淅川乌骨鸡转座子的生物学功能具有重要的意义,而且为探索基因组扩增、基因组功能及进化研究提供重要基础数据。【方法】对淅川乌骨鸡血液DNA进行全基因组重测序,采用双末端短序列比对基因组,运用RepeatModeler、TEclass、RepeatMasker等使用流程对转座子进行鉴定、构建、校正、分类和注释,获得淅川乌骨鸡整个基因组所有的转座子,对转座子的特征、分布及和基因的关系等方面进行分析。并对转座子插入的所有基因进行GO和KEGG数据库富集,分别结合GO和KEGG注释结果对功能进行描述。【结果】经鉴定、分类和注释,淅川乌骨鸡共鉴定到370 252个转座子序列,分为19个超家族,主要为CR1、TcMar-Mariner、ERVL、ERV1等超家族,进一步说明淅川乌骨鸡转座子类型主要是逆转录转座子;转座子的数目和染色体长度有关,染色体越长,转座子数目越多。在基因组中转座子和基因的密度成反比,基因密集的区域中转座子密度较低;基因组内转座子的插入并非随机的,LTR/...     

大丽轮枝菌（Verticillium dahliae VdLs.17）分泌组预测及分析

 【目的】预测并分析大丽轮枝菌基因组范围内的分泌蛋白,为大丽轮枝菌分泌蛋白致病机理的研究奠定基础。【方法】利用已公布的大丽轮枝菌全基因组序列,组合使用生物信息学软件SignalP、TargetP、TMHMM、Big-pi和PROSITE,预测大丽轮枝菌基因组范围内所有分泌蛋白,定义为分泌组。统计分析分泌组中蛋白N-端信号肽特点;应用碳水化合物活性酶类数据库和病原菌-寄主互作蛋白数据库对分泌组蛋白进行注释,预测分泌组中潜在果胶酶、纤维素酶和病原菌寄主互作蛋白;利用真菌激发子的保守结构域,预测分泌组中潜在的激发子蛋白集;应用BLASTP程序比较分析大丽轮枝菌和黑白轮枝菌分泌组,获得大丽轮枝菌相对于黑白轮枝菌特异的分泌蛋白。【结果】大丽轮枝菌分泌组共有922个蛋白。信号肽分析表明,以19个氨基酸为信号肽的蛋白数目最多,非极性氨基酸丙氨酸的出现频率最高,而有带电侧链的氨基酸天冬氨酸和谷氨酸的出现频率最低,信号肽的-3和-1位置上的氨基酸相对保守。大丽轮枝菌分泌组含有158个潜在的碳水化合物活性酶类,其中,包括10个果胶水解酶和14个果胶裂解酶;190个潜在的病原菌-寄主互作蛋白、97个含有RxLx[EDQ]模体的蛋白和52个富含半胱氨酸的小分子量分泌蛋白;58个相对于黑白轮枝菌分泌组特异的蛋白。【结论】本文建立了预测大丽轮枝菌分泌组蛋白的方法。分泌组蛋白信号肽长度具有高度的变异性,氨基酸组成多为脂肪族氨基酸,序列在C-端结构域较为保守。分泌组中包含大量潜在的果胶降解酶、病原菌-寄主互作蛋白、RxLx[EDQ]模体蛋白和富含半胱氨酸的小分子量蛋白等致病相关蛋白。     

荸荠优良品种‘桂蹄3号’球茎发育过程转录组测序分析

【目的】开展荸荠球茎发育过程转录组测序研究,为研究荸荠球茎发育过程相关基因表达信息提供参考。【方法】运用高通量测序技术,对荸荠球茎不同发育时期进行转录组测序研究。【结果】组装共得到223 182条转录本和90 542条Unigene,平均长度为809 bp,N50为1119。组装完整性较高,效果较好。对所得Unigene进行不同数据库注释,共有50 583条Unigene成功注释到7个数据库(NR、GO、NT、Pfam、KEGG、KOG及Swiss-prot)。对差异表达基因分析,发现球茎膨大初期的差异表达基因数目最多,为最活跃阶段。GO功能富集分析结果表明,33 205个基因获得功能注释,分为分子功能、细胞组分和生物学过程等3大类和54个亚类。COG功能分类结果表明,17 743个基因分布于25个功能区域,其中碳水化合物代谢占重要地位。KEGG代谢通路注释结果表明有20 667个基因获得功能注释,共有116条代谢途径,其中淀粉-蔗糖代谢占主要作用。【结论】利用高通量转录组测序技术首次建立了荸荠优良品种‘桂蹄3号’球茎的转录组数据库,为进一步研究荸荠球茎淀粉生物合成相关基因的功能及形成的分子机制提供了数据基础。     

红橘响应褐斑病菌侵染的转录组学分析

【目的】明确红橘(Citrus reticulata Blanco, cv. Hongjv)接种褐斑病致病菌——链格孢菌橘致病型(Alternaria alternata tangerine pathotype)后基因种类和表达量在转录水平的变化规律,确定红橘响应该致病型侵染的关键基因。【方法】采用链格孢菌橘致病型接种红橘离体叶片,28 h后选取感病叶片和未接种叶片提取RNA,进行转录组高通量测序,然后利用生物信息学分析,以甜橙基因组为参考,以|log2 fold change|≥1,q-value≤0.01为阈值选取感病和健康红橘叶片转录组的差异表达基因,应用GO数据库对差异表达基因(differentially expressed gene,DEG)进行功能分类,KEGG分析代谢途径,MapMan软件分析生物胁迫信号通路相关基因的表达变化。采用qRT-PCR方法对测序结果进行验证。【结果】红橘接种链格孢菌橘致病型28 h后产生大量与胁迫相关的差异表达基因,获得上调差异基因5 173个,下调差异基因6 555个。GO功能分类显示差异基因主要与蛋白结合、膜、氧化还原过程等相关。通过KE...     

基于简化基因组测序技术的番茄雄性不育基因定位

【目的】获得番茄雄性不育基因的候选基因,提高番茄雄性不育基因定位的准确性,为该雄性不育基因的克隆及功能机制研究提供理论基础。【方法】以一个与苗期绿茎基因紧密连锁的花粉败育型材料为母本,以一个苗期紫茎可育系为父本,通过集群分离(BSA)法分别建立30个F2单株的不育和可育DNA池,基于简化基因组测序(SLAF-Seq)技术检测足够量的简化基因组扩增片段,通过关联分析,获得目标雄性不育基因的候选区域,并利用GO数据库对关联区域内的基因进行功能注释。【结果】通过SLAF测序共获得20.27 Mb的序列量;共获得103 250个SLAF标签,标签整体平均深度达160.39倍;完成了2个混合池间SLAF标签标记的多态性分析,共获得多态性标记6 892个;利用2个混合池进行关联分析,共得到8个与目标雄性不育基因相关的候选区域,区域平均大小为0.29Mb,区域内共有27个差异标记,146个候选基因。【结论】利用SLAF测序技术对番茄雄性不育基因进行了初步定位。     

小麦叶锈菌与TcLr19非亲和互作cDNA文库的EST分析

 【目的】建立小麦叶锈菌与小麦非亲和互作的基因表达数据库,为从分子水平阐明小麦抗叶锈病机制奠定基础。【方法】对叶锈菌诱导的TcLr19小麦叶片cDNA文库随机测序,利用BLASTx对所得序列进行功能注释,按照Bevan的植物基因功能分类标准进行分类。采用RT-PCR技术,以Actin为参照,对2条信号转导相关基因和3条抗病与防御相关基因进行表达分析。【结果】随机对cDNA文库中756个阳性克隆测序,获得649条高质量EST。对EST序列聚类拼接获得472条非重复序列(Unisequence)。功能分类结果表明,25.2%为未知功能蛋白,21.0%与能量代谢基因相关,17.8%与抗病防御基因及信号转导基因相关,其余EST分别与转录、转运、蛋白代谢等功能基因相关。RT-PCR分析结果表明同一基因在叶锈菌侵染后不同时间点的小麦叶片中表达量不一致,且各基因有不同的表达模式。【结论】推测钙转运ATP酶,谷胱甘肽-S-转移酶,蛋白激酶Pti1,NBS-LRR抗病蛋白和分子伴侣DnaJ等参与了小麦与叶锈菌的非亲和互作过程。该cDNA文库可用于了解小麦与叶锈菌非亲和互作过程中的功能基因,为建立研究病原菌基因及小麦抗病基因数据库奠定基础。     

凡纳滨对虾转录组测序分析及肌肉生长发育相关基因的筛选

【目的】筛选出凡纳滨对虾生长发育的功能基因及代谢调控网络,揭示其生长发育的分子机制,为后续开展凡纳滨对虾分子生物学研究提供宝贵的基因数据来源。【方法】以快速生长群体和慢速生长群体的凡纳滨对虾肌肉组织为研究材料,通过Illumina HiSeqTM2500平台对构建的cDNA文库进行高通量测序分析,以StringTie进行拼接组装后利用DESeq2筛选差异表达基因,并基于KOG、GO、nr、COG、Swiss-Pro、KEGG和Pfam等数据库进行差异表达基因功能注释分析。【结果】经拼接组装共获得53458844条Clean reads,各样品Clean reads的Q30均在93.00%以上;Illumina测序获得的Clean reads与参考基因组的比对效率在87.75%~87.80%,说明转录组测序数据真实可靠。采用SnpEff进行SNP/InDel变异注释分析,结果显示,SNP位点中以A>G、G>A、C>T和T>C等4种类型的数量较多（14950~21562个）,InDel位点则以SYNONYMOUS_CODING的数量最多（33630个）。使用StringTie对Mapped reads（比对到参考基因组的Reads）进行拼接,共发掘到4607个新基因,分别输入COG、GO、KEGG、KOG、Pfam、Swiss-Prot、eggNOG和nr数据库中进行序列比对,最终发现共有1098个新基因被注释,以被nr数据库注释的新基因数量最多（1077个）,而被COG数据库注释的新基因数量最少（416个）。基于Q-value<0.05且Fold Change>2的筛选条件,共获得1408个差异表达基因（661个为显著上调表达基因,747个为显著下调表达基因）;1408个差异表达基因被注释到53个GO功能条目中,其中,22条被注释到生物学过程（Biological process）,16条被注释到细胞组分（Cellular component）,15条被注释到分子功能（Molecular function）;KEGG信号通路富集分析发现3条重要的信号通路,分别是溶酶体通路（Lysosome）、氨基糖和核苷酸糖新陈代谢（Amino sugar and nucleotide sugar metabolism）及鞘脂类代谢（Sphingolipid metabolism）。在溶酶体通路中,CTSL基因、Nramp基因、MyoG基因和Myf5基因在凡纳滨对虾快速生长群体肌肉组织中呈上调表达,而Trypsin基因呈下调表达。【结论】通过转录组测序分析从快速生长群体和慢速生长群体的凡纳滨对虾肌肉组织中筛选出1408个差异表达基因（661个为显著上调表达基因,747个为显著下调表达基因）,主要富集在溶酶体、氨基糖和核苷酸糖新陈代谢及鞘脂类代谢等通路上,在对虾肌肉生长发育过程中发挥重要作用。