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

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

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

【目的】筛选出凡纳滨对虾生长发育的功能基因及代谢调控网络,揭示其生长发育的分子机制,为后续开展凡纳滨对虾分子生物学研究提供宝贵的基因数据来源。【方法】以快速生长群体和慢速生长群体的凡纳滨对虾肌肉组织为研究材料,通过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个为显著下调表达基因）,主要富集在溶酶体、氨基糖和核苷酸糖新陈代谢及鞘脂类代谢等通路上,在对虾肌肉生长发育过程中发挥重要作用。     

茶树炭疽菌N425全基因组测序与基因功能注释

【目的】对茶树炭疽菌N425进行全基因组测序和基因功能注释并从中初筛出致病相关基因。【方法】通过Illumina HiSeq2500 PE150平台对N425菌株进行全基因组测序和组装,并利用NR、KEGG、KOG等功能基因数据库进行基因预测和功能注释。【结果】茶树炭疽菌N425基因组全长约56.3 Mbp,G+C含量为53.2%,共预测出10 157个蛋白编码基因和若干各类非蛋白编码序列。其中,1 356个基因在病原与宿主互作数据库（PHI）中得到注释,包括140个注释为致病力丧失的基因;720个基因在碳水化合物酶数据库（CAZy）中得到注释;302个基因被预测为次生代谢物相关基因。这些基因中有涉及cAMP-PKA、MAPK等信号通路以及降解宿主细胞壁等致病相关过程,是潜在的致病相关基因。【结论】本研究成功获得茶树炭疽菌N425全基因组序列和基因功能注释,并从中挖掘出潜在致病相关基因,为后续展开对茶树炭疽菌侵染茶树的致病分子机制研究奠定基础。     

一株牦牛源产细菌素植物乳杆菌SWUN5815全基因组测序及序列分析

旨在分析一株牦牛源植物乳杆菌SWUN5815的全基因组序列测序,并对该菌株的功能特性基因进行挖掘.基于PacBio RSII测序平台对乳杆菌SWUN5815进行全基因组测序分析和GO、KEGG、COG和NR数据库进行基因组基本功能注释.结果 表明:1)植物乳杆菌SWUN5815的基因组大小为3.27 Mb,GC含量44...     

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

【目的】对青花菜花蕾进行转录组测序分析，并挖掘与蜡粉合成相关基因，为探明青花菜花球表面蜡粉形成的分子机制提供理论参考。【方法】分别提取野生型和蜡粉缺失型青花菜花球总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家族在调控青花菜蜡粉合成中发挥重要作用。蜡粉合成过程中相关酶基因的差异表达是调控青花菜蜡粉合成的关键，尤其是野生型和蜡粉缺失突变体中特异性表达的差异表达基因，可作为后续研究青花菜花球表面蜡粉形成分子机制的对象。     

生防菌贝莱斯芽孢杆菌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）、杆菌溶素（Bacilysin）、大环内酰亚胺H （Macrolactin H）、杆菌烯（Bacillaene）和艰难菌素（Difficidin）等抑菌物质,其中多数抑菌物质是通过非核糖体途径和聚酮合酶途径合成,且其中6个基因簇在贝莱斯芽孢杆菌中高度保守。【结论】菌株MC2-1的基因组中含有多种次级代谢产物合成基因簇,可编码合成已知的抑菌活性物质,同时还存在降解病原真菌细胞壁的相关基因。推测菌株MC2-1可通过分泌抑菌次生代谢物及相关降解酶达到防治植物病害的效果,在农业生物防治中具有较好的应用前景。     

基于转录组学测序的长白落叶松材性表达基因1）

基于边合成边测序（ Sequencing By Synthesis ,SBS）技术,使用Illumina HiSeq2500高通量测序平台对长白落叶松cDNA文库进行测序。共获得非重复序列基因58683条,总长度为55283938 bp。将得到的转录本使用BLAST软件将非重复序列基因序列与NR、Swiss－Prot、GO、COG、KEGG数据库比对,通过选择BLAST参数E－value不大于10－5,最终获得29350个有注释信息的基因序列,其中19292个转录本注释GO编号,有9112个转录本在COG数据库中被注释并分为25个功能分类,同时注释到120条KEGG代谢途径。在得到注释的转录本中搜索到木质素相关基因共68条,其中与木质素合成相关基因29条,与木质素分解代谢相关基因27条;找到了147条纤维素相关基因,其中与纤维素合成相关基因103条,与纤维素分解代谢相关基因44条。     

荸荠优良品种‘桂蹄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号’球茎的转录组数据库,为进一步研究荸荠球茎淀粉生物合成相关基因的功能及形成的分子机制提供了数据基础。     

卵形鲳鲹卵巢不同发育时期的转录组测序分析

【目的】鉴定筛选出与卵形鲳鲹卵巢发育相关的候选基因及信号通路,为揭示其卵巢性成熟过程的分子机制打下基础。【方法】挑选卵巢发育处于?期和Ш期的雌性卵形鲳鲹,分别构建卵形鲳鲹卵巢?期和Ш期的cDNA文库,采用Illumina HiSeqTM 2500进行转录组测序,经过滤、质量控制及拼接组装后获得的Unigenes在七大数据库（Nr、Nt、Pfam、KOG/COG、Swiss-Prot、KEGG和GO）中进行比对;通过FPKM及DEGseq筛选出差异表达基因,以GOseq和KOBAS对差异表达基因分别进行功能注释及信号通路富集分析,并采用MISA和GATK3进行SSR鉴定及SNP分析。【结果】卵形鲳鲹卵巢组织转录组测序获得的325156432条Raw reads,经过滤筛选得到317206752条Clean reads,拼接组装后得到59554条Unigenes;69.65%的Unigenes在Nr、Nt、Pfam、KOG/COG、Swiss-Prot、KEGG和GO等七大数据库中注释成功,其中有24599条Unigenes被注释到GO数据库,15997条Unigenes被注释到KEGG数据库。在卵形鲳鲹卵巢组织的2个发育时期共鉴定获得56115个基因,经差异表达分析后获得17737个差异基因,其中8169个基因在卵巢Ш期上调表达、9568个基因在卵巢Ш期下调表达。GO功能注释分析发现,卵形鲳鲹卵巢差异表达基因主要注释在细胞过程、氮化合物代谢过程、初级代谢过程、核、核部分、离子结合及水解酶活性等条目上;而KEGG信号通路富集分析结果显示,17737个差异表达基因显著富集在318条代谢途径上,其中前20条KEGG信号通路包括2-氧代羧酸代谢、PI3K-Akt信号通路、甲状腺激素信号通路、磷脂酶D信号通路、Fc εRI信号通路和细胞周期等。卵形鲳鲹卵巢转录组（59554条Unigenes）中共存在30133个SSRs和82490个SNPs。【结论】GnRHR、FSHR、FSHβ、CYP11A、SIRT3和PEG3等差异表达基因及PI3K-Akt信号通路和VEGF信号通路等与卵形鲳鲹卵巢的发育密切相关,共同调节卵巢的发育与成熟,在卵巢性成熟过程中发挥重要作用。     

基于纳米孔全长转录组数据完善东方蜜蜂微孢子虫的基因组注释

【目的】利用已获得的纳米孔全长转录组数据对现有的东方蜜蜂微孢子虫（Nosema ceranae）参考基因组的基因序列和功能注释进行完善。【方法】采用TransDecoder软件预测东方蜜蜂微孢子虫基因的开放阅读框（open reading frame,ORF）及相应的氨基酸。利用gffcompare软件将全长转录本与参考基因组注释的转录本进行比较,对基因组注释基因的非编码区向上游或下游延伸,修正基因的边界。利用MISA软件鉴定长度在500 bp以上的全长转录本的简单重复序列（simple sequence repeat,SSR）位点,包括单核苷酸重复、双核苷酸重复、三核苷酸重复、四核苷酸重复、五核苷酸重复、六核苷酸重复、混合SSR等类型。通过Blast工具将鉴定到的新基因和新转录本比对Nr、KOG、eggNOG、GO和KEGG数据库,从而获得功能注释。【结果】共预测出2 353个完整ORF,其中长度分布在0—100个氨基酸的ORF最多,占总ORF数的72.12%。共对东方蜜蜂微孢子虫的2 340个基因进行了结构优化,其中5′端延长的基因有1 182个,3′端延长的基因有1 158个。共鉴定到1 658个SSR,其中单核苷酸重复、双核苷酸重复、三核苷酸重复、四核苷酸重复的数量分别为1 622、23、7和6个;单核苷酸重复类型的SSR密度最大,达到182.32个/Mb,其次为混合SSR、双核苷酸重复和三核苷酸重复,分别达到6.90、2.78和0.73个/Mb。共鉴定出954个新基因,其中分别有951、333、371、422和321个新基因可注释到Nr、KOG、eggNOG、GO和KEGG数据库。此外,还鉴定出6 164条新转录本,其中分别有6 141、2 808、2 932、3 196和2 585条新转录本可注释到Nr、KOG、eggNOG、GO和KEGG数据库。新基因和新转录本注释数量最多的物种均为东方蜜蜂微孢子虫,其次是蜜蜂微孢子虫（Nosema apis）。【结论】研究结果较好地完善了现有的东方蜜蜂微孢子虫参考基因组已注释基因的序列和功能注释,并补充和注释了大量参考基因组未注释的新基因和新转录本。     

福清山羊与努比亚黑山羊发情期卵巢组织RNA-Seq分析

【目的】 研究福清山羊与努比亚黑山羊发情期卵巢组织转录组差异表达水平,一方面为山羊繁殖性状形成的分子机制提供理论依据,另一方面为利用努比亚黑山羊杂交改良福清山羊提供可加快遗传进展的分子标记。【方法】 利用转录组测序方法对福清山羊和努比亚黑山羊发情期卵巢组织进行研究,筛选品种间的差异表达基因（differentially expressed genes,DEGs）,并对DEGs功能进行注释和若干基因荧光定量PCR（quantitative real-time PCR,qRT-PCR）验证;同时通过与参考基因组比对,分析和筛选测序样品中存在的SNP/InDel。【结果】 6个样品共得到46.68Gb Clean Data,DESeq分析发现了福清山羊和努比亚黑山羊发情期卵巢组织的DEGs149个（包含25个新转录本）,其中表达上调53个,表达下调96个;初步认为输卵管素基因（oviductin,OVN）、类固醇合成快速调节蛋白基因(steroidogenic acute regulatory protein,STAR)、早期生长应答1基因(early growth response 1,EGR1)可作为福清山羊繁殖性能的候选基因。149个DEGs中的108个基因能被GO（gene ontology）数据库注释,30个DEGs能够被COG（Cluster of Orthologous Groups of proteins）数据库注释,91个DEGs能够被KEGG（kyoto encyclopediaof genes and genomes）数据库注释。KEGG通路分析表明DEGs共富集到21条信号通路中,6条通路显著富集。经过进一步的与参考基因组序列比对分析,共发掘1 506个新转录本。经qRT-PCR验证,所选基因（转录本）表达变化模式与转录组测序结果一致,表明测序结果可靠。【结论】 在转录组水平上筛选出了福清山羊和努比亚黑山羊发情期卵巢组织的149个DEGs,发掘了1 506个新转录本,初步揭示了OVN、STAR和EGR1在山羊繁殖过程中可能发挥重要作用,为进一步探索山羊繁殖性状相关机理提供参考依据。     

基于转录组学测序的长白落叶松材性表达基因

基于边合成边测序(Sequencing By Synthesis,SBS)技术,使用Illumina Hi Seq2500高通量测序平台对长白落叶松c DNA文库进行测序。共获得非重复序列基因58 683条,总长度为55 283 938 bp。将得到的转录本使用BLAST软件将非重复序列基因序列与NR、Swiss-Prot、GO、COG、KEGG数据库比对,通过选择BLAST参数Evalue不大于10-5,最终获得29 350个有注释信息的基因序列,其中19 292个转录本注释GO编号,有9 112个转录本在COG数据库中被注释并分为25个功能分类,同时注释到120条KEGG代谢途径。在得到注释的转录本中搜索到木质素相关基因共68条,其中与木质素合成相关基因29条,与木质素分解代谢相关基因27条;找到了147条纤维素相关基因,其中与纤维素合成相关基因103条,与纤维素分解代谢相关基因44条。     

菊叶香藜转录组单核苷酸多态性(SNP)信息挖掘及功能注释分析

【目的】挖掘菊叶香藜转录组数据中的SNP信息,对具有SNP位点的unigene(SNP-unigene)进行功能注释分析,为菊叶香藜SNP分子标记的开发提供理论依据。【方法】利用Samtools和VarScan v.2.2.7软件寻找候选的SNP位点,并对SNP-unigene进行GO注释、COG注释和KEGG注释。【结果】在菊叶香藜花和叶转录组中,分别鉴定出了889个和673个SNP位点,转换分别占63.0%和62.7%,颠换分别占37.0%和37.3%,转换和颠换的比值分别为1.70和1.68,非编码区分别占21.15%和21.10%,编码区分别占67.27%和67.46%。菊叶香藜所有SNP位点分布于643条SNP-unigene上,功能注释结果显示,总共有343条SNP-unigene具有GO注释,370条SNP-unigene具有COG注释,232条SNP-unigene具有KEGG注释,这些SNP-unigene涉及较多的功能主要与代谢、核糖体、次生代谢产物的生物合成相关。【结论】本研究通过大规模地筛选菊叶香藜中的候选SNP位点,可以为研究菊叶香藜基因分型、构建遗传图谱等方面奠定基础,对于开发利用菊叶香藜植物资源具有重要的价值。     

枯草芽孢杆菌NCD-2对棉花根系分泌物L-脯氨酸响应的转录-蛋白质组学联合分析

【目的】棉花根系分泌物L-脯氨酸是影响生防菌株定殖的关键因素之一,前期研究发现L-脯氨酸能够提高枯草芽孢杆菌（Bacillus subtilis）NCD-2生物膜形成能力。通过高通量测序技术分析L-脯氨酸调控枯草芽孢杆菌NCD-2生物膜形成和生防潜力相关的基因,为深入认识棉花根系分泌物与生防菌株的分子互作关系打下基础。【方法】以枯草芽孢杆菌NCD-2菌株为供试材料,外源添加浓度为10 mg·mL^-1的L-脯氨酸共培养24 h后,分别进行转录组（RNA-seq）和同位素标记相对定量蛋白质组技术（iTRAQ）分析,对所获得的转录组和蛋白质组数据进行生物信息学分析,利用实时荧光定量PCR（RT-qPCR）验证不同代谢通路中部分差异基因的表达情况。 【结果】转录组分析发现,L-脯氨酸和NCD-2菌株共培养后,获得1 071个差异表达基因（DEG）,其中602个基因上调,469个基因下调。GO分析发现在生物过程、细胞组分和分子功能方面分别有49、14和30个功能条目显著富集。KEGG代谢通路主要富集在化合物代谢、鞭毛组装、细菌运动和趋化作用。蛋白质组学分析发现,筛选到211个差异表达蛋白（DEP）,其中118个蛋白上调,93个蛋白下调。GO分析发现在生物过程和分子功能方面分别有13和8个功能条目显著富集。KEGG代谢通路主要富集在氨基酸代谢、碳水化合物类代谢、鞭毛组装和ABC转运蛋白。进一步转录-蛋白质组学联合分析发现,在L-脯氨酸作用下,检测到共有的显著差异表达基因（或蛋白）112个,其中38个基因（或蛋白）下调,74个基因（或蛋白）上调。GO功能显著富集主要集中在营养库活性、催化活性、细胞膜、定位、细胞脂质代谢过程、氧化还原过程、sigma因子活性、转运活性、芽孢形成9个方面。KEGG代谢通路主要富集在能量代谢、ABC转运蛋白、抗生素生物合成、鞭毛组装、运动或趋化作用和双组分系统方面。RT­qPCR验证了26个差异表达显著基因,结果发现在表达量上存在一定的差异,但表达趋势与RNA­seq和iTRAQ组学分析结果基本一致。【结论】棉花根系分泌物L-脯氨酸与枯草芽孢杆菌NCD-2之间的互作存在一个复杂的生物过程,依赖于不同代谢通路网络中的多个基因。双组分系统、抗生素生物合成、能量代谢、运动或趋化作用、鞭毛组装和ABC转运蛋白通路中的差异基因（或蛋白）可能在棉花根系分泌物与枯草芽孢杆菌互作过程方面发挥着重要作用。     

金鱼源维氏气单胞菌JY01全基因组测序及生物信息分析

维氏气单胞菌(Aeromonasveronii JY01)是从福州某养殖场内金鱼烂尾、腐皮和溃疡等病灶处分离出的1株革兰氏阴性优势菌。为了解其基因组结构及功能特征,基于IlluminaPE150和Pacbio测序平台对维氏气单肥菌菌株JY01进行全基因组测序,并对测序数据进行组装和组分分析,利用生物信息学手段进行基因预测与功能注释、同时预测致病因子和耐药基因。结果表明:维氏气单肥菌菌株JY01基因组序列长度为4.97Mb,GC含量为58.17%;共预测到4 601个编码基因、31个rRNA、16个基因岛及124个tRNA;其中,3 461个基因在COG中获得注释,3 082个基因聚集到GO聚类分析中,4 345个基因在KEGG代谢通路中富集,预测维氏气单胞菌菌株JY01携带13种毒力因子包含240个相关毒力基因及9大类抗生素耐药性包含175个抗生素耐药基因。     

不同糖源对葡萄试管苗蛋白激酶相关基因表达的影响

【目的】探究不同外源糖对葡萄试管苗生长发育及蛋白激酶基因转录调控的影响,应用转录组测序挖掘蛋白质磷酸化过程中的基因,为葡萄蛋白激酶相关基因功能的验证奠定一定基础。【方法】在基本培养基中分别添加2%的蔗糖、葡萄糖和果糖,以无糖为对照,分别命名为S20、G20、F20和CK,经过37 d培养后,测定不同处理的地上和地下鲜重,并采用Illumina HiSeq ^TM 2000对各处理叶片进行转录组测序,通过综合生物信息学分析（参考基因组比对、差异基因（DEGs）筛选、COG（Cluster of Orthologous Groups of proteins）注释、GO（Gene Ontology）注释等）筛选出蛋白激酶相关基因,通过qRT-PCR分析该蛋白激酶相关基因的表达特性。【结果】F20、G20和S20处理下的葡萄（‘红地球’）试管苗与CK相比,地上鲜重具有明显差异,且F20最高,而G20地下鲜重最高。SNP统计发现,转换是主要的变异类型,颠换次之,且发生在基因间区的SNP数量最多,其次是下游;剪接位点供体和同义终止发生的基因数量最少且相等。4个样品中共获得了2 633个差异基因,3个处理与CK相比,共有差异基因180个且被聚类为3组,第一组中127个基因仅在CK中高表达,第二组19个基因仅在G20下高表达,而第三组34个基因在3个处理下表达模式不尽相同。这些共有的差异基因在COG中注释到了26个基因并分在11个功能类别中,且主要注释在一般的功能类别中。在GO分类中,共有的基因分别被注释在分子功能、生物学过程和细胞组分的14、22和13个功能类别中。共筛选出7种蛋白激酶,分别为葡萄糖激酶（Glucokinase,GK）、丝裂原活化蛋白激酶（Mitogen-activated protein kinases,MAPKs）、钙调蛋白激酶（Calcineurin protein kinase,CBL）、蛋白磷酸酶2（Protein phosphatase 2,PP2A）、己糖激酶（Hexokinase,HXK）、组氨酸蛋白激酶（Histidine protein kinase,HPK）和酪氨酸激酶（Tyrosine kinase,TK）,其不同激酶的基因在不同处理中具有各自的表达模式,经qRT-PCR验证,选择的20个差异基因中有17个基因表达与转录组测序结果相一致。【结论】在葡萄试管苗培养中,果糖较葡萄糖和蔗糖相比对生长较好。测序得出180个差异基因对3种不同糖均作出响应,这些基因在COG数据库中主要富集在膜酯转运和代谢、次级代谢物和碳水化合物的合成、转运和分解;GO中大多注释在蛋白激酶和氧化还原酶的活性中;筛选出了7种蛋白激酶,这些差异基因在数量、功能分类和代谢通路上对糖的响应各不相同。     

基于转录组测序对紫花苜蓿细胞质雄性不育系相关代谢通路的鉴定

【目的】研究紫花苜蓿(Medicago sative L.)细胞质雄性不育(CMS)的机理,为紫花苜蓿利用不育系的育种工作提供理论基础。【方法】本研究以紫花苜蓿细胞质雄性不育系MSGN1A及其保持系MSGN1B为材料,进行转录组测序(RNA-Seq)以及基因功能注释,筛选与不育系机理相关的差异表达基因及代谢通路。【结果】Illumina测序共得到紫花苜蓿95 679条功能基因,经过差异表达分析筛选出187个显著差异表达的基因(DEGs),其中包括18个上调基因,169个下调基因。采用GO、KEGG、COG注释库分别对187个DEGs进行功能注释,其主要涉及催化活性、代谢活动、信号传导机制、合成、膜功能、碳水化合物运输和代谢、能量代谢等相关通路。【结论】筛选出与紫花苜蓿细胞质雄性不育性相关的代谢通路,初步探讨了紫花苜蓿细胞质雄性不育系的分子机理。     

pH对解淀粉芽胞杆菌TR2生长和在土壤中定殖的影响

【目的】为了明确解淀粉芽胞杆菌TR2菌株在不同pH值土壤中的定殖能力及运动性、生物膜形成等相关的影响。【方法】采用平板计数法、结晶紫染色法、运动性平板法等研究不同pH条件对菌株的生长速度、趋化性、生物膜形成、涌动性及在土壤中的定殖能力。【结果】解淀粉芽胞杆菌TR2菌株适合生长和定殖的pH值范围为pH5～10,其中最适pH为7～8,且在pH为7～8的环境中，菌株有较好的涌动性运动能力。解淀粉芽孢杆菌TR2在pH7的中性环境中向不同pH环境的趋向中，向pH 8的趋向最为明显。解淀粉芽孢杆菌TR2在pH8～9形成良好复杂的生物膜，在pH8的土壤中的定殖能力好于其他处理。【结论】中性或弱碱性土壤更利于解淀粉芽胞杆菌TR2菌株的生长与定殖，为解淀粉芽胞杆菌TR2菌株的开发和利用提供了试验依据。     

2种不同生长规格墨瑞鳕肌肉组织转录组测序分析

【目的】掌握不同生长规格墨瑞鳕个体间差异表达基因的表达特点,为其功能相关基因深度挖掘及分子遗传育种提供科学依据。【方法】挑选同一养殖条件下极大个体和极小个体的墨瑞鳕,构建肌肉组织cDNA文库后,采用Illumina HiSeqTM 4000测序平台对存在生长差异的墨瑞鳕肌肉组织进行转录组测序分析,获得的Unigenes在Nr、Nt、Pfam、KOG/COG、Swiss-Prot、KEGG和GO等数据库中进行比对;通过FPKM及DEGseq筛选出差异表达基因,以GOseq和KOBAS对差异表达基因分别进行GO功能注释及KEGG信号通路富集分析,并采用MISA进行SSR鉴定分析。【结果】从墨瑞鳕肌肉组织中共测序获得39749条Unigenes,其长度范围在301~55230 bp,平均长度为1705 bp。注释到Nt、Nr、Swiss-Prot、Pfam数据库的Unigenes分别有27046、20824、18268和17772条,在7个数据库中均得到注释的Unigenes共计6742条,占Unigenes总数的16.96%。根据差异表达基因筛选条件P<0.05且|log₂ Fold Change|>1,共筛选出722个差异表达基因,其中上调基因308个、下调基因414个。差异表达基因GO功能注释分析结果表明,注释基因数目较多的GO功能条目包括细胞过程、代谢过程、膜、细胞器及结合等;KEGG信号通路富集分析发现,差异表达基因被成功富集到234条信号通路上,主要涉及磷脂酰肌醇3激酶/蛋白激酶信号通路、MAPK信号通路、胰岛素信号通路及FoxO信号通路等。在39749条Unigenes中鉴定筛选出22120个SSRs,占Unigene总数的55.65%,SSR的平均间距为3063 bp。【结论】基于转录组测序分析获得的墨瑞鳕肌肉组织差异表达基因以发挥结合、细胞过程及代谢过程等功能为主,且主要富集在PI3K-Akt信号通路、核糖体信号通路、FoxO信号通路及细胞凋亡等能量代谢相关通路上,通过共同协调而对墨瑞鳕的生长发育起调控作用。     

旱柳转录组测序及生物学分析

对盐敏感旱柳沿江柳和耐盐旱柳9901的杂交F_1代根部进行了RNA-Seq测序和分析,共获取了107 950条Unigene,平均长度为1 076.96 bp。通过与COG、GO等8个数据库比对,60 848个基因获得注释信息,其中38 182条Unigene在GO数据库中获得注释,24 101条Unigene在KEGG数据库中获得注释。GO和KEGG富集分析结果表明,差异表达基因主要调节核糖体代谢、植物激素信号转导等生物学功能。