猴樟叶片转录组测序及盐碱胁迫响应基因表达模式研究

为了解猴樟响应盐碱胁迫机制,以2个耐盐碱猴樟品系和1个盐碱敏感型猴樟品系3年生种苗叶片为材料,利用Illumina Hiseq^TM平台对其进行转录组测序,对6个精氨酸和脯氨酸通路相关基因进行qRT−PCR分析。结果表明：研究共获得72369条Unigenes,其中39156条Unigenes在NR、NT、KEGG、SwissProt、PFAM、GO、KOG数据库中获得注释;HDHZ vs. HJHZ 组有5415个DEGs,其中4003个得到GO注释;YYHZ vs. HJHZ 组有4665个DEGs,其中331个得到GO注释。此外,KEGG通路分析发现,HDHZ vs. HJHZ组上调基因主要富集在氨基酸代谢和生物碱合成相关通路,下调基因主要富集在倍半萜和三萜的生物合成,以及光合作用相关通路;YYHZ vs. HJHZ组上调基因主要富集在倍半萜和三萜的生物合成及类黄酮生物合成通路;下调基因主要富集在色素代谢和酚类物质代谢相关通路。结合6个差异表达基因的qRT−PCR结果,生物碱类、精氨酸和脯氨酸代谢可能在猴樟响应碱胁迫过程中起主要作用,所得基因的表达模式与RNA−Seq数据结果一致,证实RNA−Seq数据准确。     

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

【目的】鉴定筛选出与卵形鲳鲹卵巢发育相关的候选基因及信号通路,为揭示其卵巢性成熟过程的分子机制打下基础。【方法】挑选卵巢发育处于?期和Ш期的雌性卵形鲳鲹,分别构建卵形鲳鲹卵巢?期和Ш期的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信号通路等与卵形鲳鲹卵巢的发育密切相关,共同调节卵巢的发育与成熟,在卵巢性成熟过程中发挥重要作用。     

盐胁迫对南方泡桐基因表达的影响

为阐明泡桐多倍体的耐盐分子机制,以南方泡桐二倍体及对应的四倍体为试验材料,以白花泡桐基因组为参考序列,利用RNA-Seq测序技术,研究了盐处理前后南方泡桐四倍体和对应二倍体的基因表达变化。通过比对分析,共鉴定出与盐胁迫相关的差异表达基因2 940个。对这些差异基因进行GO功能分类,结果显示差异基因共参与了39个分类,集中在细胞、结合、催化等分类功能区的数量最多。KEGG代谢通路分析发现,这些差异表达基因主要参与了Plant hormone signal transduction,Carbon metabolism,Biosynthesis of amino acid等,其中编码MYB、WRKY、bZIP、ATPase等的基因差异表达显著。表明这些基因可能是潜在的盐胁迫响应基因。采用qRTPCR技术验证了随机挑选的9个差异表达基因,结果与转录组测序数据趋势一致。     

西方蜜蜂工蜂不同虫态发育的转录组学分析

【目的】基于转录组学对西方蜜蜂工蜂不同虫态间的差异表达基因（DEGs）进行筛选和功能注释分析,揭示与工蜂生长发育相关的信号通路,为深入解析工蜂生长发育的分子调控机理提供基础数据。【方法】以西方蜜蜂工蜂的3日龄幼虫、1日龄蛹和1日龄羽化工蜂3个虫态为研究对象,利用llumina NovaSeq 6000平台进行转录组测序,采用DESeq2筛选不同虫态样品间的表达差异基因,然后分别进行GO功能注释分析及KEGG信号通路富集分析,并通过实时荧光定量PCR进行验证。【结果】经转录组测序,在西方蜜蜂工蜂3日龄幼虫与1日龄蛹间筛选出4823个差异表达基因（51.86%上调,48.14%下调）,在1日龄蛹与1日龄羽化工蜂间筛选出3295个差异表达基因（57.51%上调,42.49%下调）,在3日龄幼虫与1日龄羽化工蜂间筛选出5267个差异表达基因（52.95%上调,47.05%下调）。GO功能注释分析结果显示,3日龄幼虫与1日龄蛹间的差异表达基因注释到43个GO功能条目,1日龄蛹与1日龄羽化工蜂间的差异表达基因注释到45个GO功能条目,3日龄幼虫与1日龄羽化工蜂间的差异表达基因注释到44个GO功能条目,主要涉及细胞过程、细胞部分及结合等。KEGG信号通路富集分析发现,3日龄幼虫与1日龄蛹间有2905个差异表达基因富集到332条KEGG信号通路上,其中17条KEGG信号通路呈显著富集,涉及核糖体、氧化磷酸化和昆虫激素生物合成等;1日龄蛹与1日龄羽化工蜂间有1644个差异表达基因富集到331条KEGG信号通路上,其中45条KEGG信号通路呈显著富集,涉及氧化磷酸化、生热作用和胰岛素分泌等;3日龄幼虫与1日龄羽化工蜂间有2958个差异表达基因富集到337条KEGG信号通路上,其中14条KEGG信号通路呈显著富集,涉及核糖体、蛋白酶体和胰岛素分泌等。6个随机挑选差异表达基因的实时荧光定量PCR检测结果与转录组测序结果相符,即转录组测序结果可靠。【结论】昆虫激素生物合成通路相关差异表达基因调控与西方蜜蜂工蜂各虫态JH滴度变化规律一致,氧化磷酸化信号通路则与各虫态的营养摄入和活动行为相关,而胰岛素分泌通路涉及各虫态的营养调控、脂肪体合成及细胞凋亡。可见,昆虫激素生物合成、胰岛素分泌和氧化磷酸化3种信号通路在西方蜜蜂工蜂幼虫、蛹和成虫的发育调控中发挥着重要作用。     

百合叶烧病发病过程叶片细胞和差异表达基因分析

[目的]明确百合叶烧病发病过程中百合叶片细胞结构的变化,探究百合叶烧病发病的分子调控机理.[方法]通过扫描电镜和透射电镜观察东方百合Tarrango正常叶片、轻度叶烧叶片和重度叶烧叶片超微结构.并通过比较东方百合Tarrango正常叶片、叶烧叶片、正常叶片喷钙和叶烧叶片喷钙4种处理的转录组测序数据,对差异表达基因进行序...     

谷子穗顶端败育突变体sipaa1的表型分析和基因定位

【目的】穗发育对于农作物产量至关重要,而穗顶端败育谷子产量下降的重要原因之一。通过挖掘谷子穗顶端败育的相关基因,探求谷子穗顶端发育的生物学通路,以期为谷子穗发育遗传机理研究提供理论基础。【方法】利用化学诱变剂甲基硫酸乙酯（ethyl methyl sulfonate,EMS）对野生型豫谷一号（Yugu1）进行诱变,在其后代中发现了一个可以稳定遗传的穗顶端败育的突变体,命名为sipaa1,同时对该突变体的农艺性状进行鉴定。以突变体sipaa1母本,SSR41父本构建的F₂定位群体为材料进行遗传分析及图位克隆,确定基因所属染色体以及在该染色体上的位置。对突变体sipaa1和野生型Yugu1的BC₁F₂进行高通量测序,挖掘定位区间内的候选基因,根据候选基因在谷子不同组织部位表达量的差异,找出在穗部高表达的候选基因。对孕穗期的Yugu1和sipaa1进行转录组测序,寻找差异表达基因并分析差异表达基因富集的生物学通路。【结果】与Yugu1相比,突变体sipaa1的平均株高略有增高,增幅不显著,叶长、叶宽分别降低了10.66%和5.08%。突变体的表型变异主要集中在穗部,最突出的表现是穗顶端小花发育异常,谷穗长和谷穗粗分别降低了11.36%和16.12%,单株穗重、谷码数、单穗粒重及千粒重分别降低了30.02%、32.58%、30.55%和18.18%。通过对sipaa1×SSR41的F₂代群体中正常株与突变株的遗传分析表明该突变为隐性单基因控制。经图位克隆将突变基因定位于第1染色体Indel标记1-9.23与1-9.333之间约100 kb的范围内。结合高通量测序数据库,在该定位区间筛选到6个在穗部高表达的候选基因。转录组测序发现,在突变体与野生型之间存在2 768个上调表达基因,507个下调表达基因,且定位区间内有2个差异表达基因主要与激素信号转导、外界胁迫响应、植物-病原互作等生物学通路有关。【结论】谷子穗顶端败育突变体sipaa1由隐性单基因控制,突变基因位于第1染色体Indel标记1-9.23与1-9.333之间,转录组测序与基因功能分析发现了2个在穗部高表达且与植物花器官发育及胁迫响应密切相关的候选基因,候选基因可能通过对激素、胁迫响应,以及细胞程序性死亡等相关通路调控谷子穗顶端败育。     

干旱胁迫下花生转录组与miRNA测序及相关基因的表达

[目的]探明干旱胁迫下花生转录组与miRNA测序及相关基因的表达,为深入分析花生抗旱分子机制提供理论依据.[方法]采用转录组测序和miRNA测序技术,对桂花37花生在15％ PEG6000模拟干旱胁迫下基因和miRNA的表达模式进行分析.[结果]转录组测序发现有13个显著持续上调的差异基因,其中7个为HSP分子伴侣家族...     

利用RNA-Seq鉴定甘蓝型油菜叶片干旱胁迫应答基因

【目的】利用RNA Sequencing（RNA-Seq）技术比较2种不同生长条件下甘蓝型油菜苗期叶片转录组,鉴定油菜叶片干旱胁迫应答相关基因,从转录组水平揭示油菜适应干旱胁迫环境的分子机制。【方法】提取正常生长（ZY）和自然失水处理（ZY8D）的六叶期甘蓝型油菜中油821的叶片总RNA,以Illumina Hiseq 2000平台进行RNA-Seq分析。利用NGSQCTookit v2.3.3去除低质量和包含模糊碱基的reads。以甘蓝型油菜亲本物种白菜染色体v1.5和甘蓝Scaffold v1.0为参考序列,采用TopHat2-Cufflinks-Cuffmerge-Cuffdiff标准流程进行差异表达基因（differential expressed genes,DEGs）筛选。对上调和下调DEGs分别采用Cytoscape v3.1.0中的BiNGO和KOBAS2.0进行基因本体（gene ontology,GO）和京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes,KEGG）代谢途径富集分析。选择上调和下调DEGs各3个,以实时荧光定量PCR（quantitative real-time PCR,qRT-PCR）验证RNA-Seq结果的可靠性。【结果】过滤低质量reads后,ZY和ZY8D分别保留了26 192 312和28 378 899对高质量reads用于DEGs筛选,其中86.6%和85.8%的reads能准确比对到参考序列上,说明RNA-Seq结果和参考序列可靠。DEGs鉴定结果表明3 657个基因受干旱胁迫诱导差异表达,其中上调表达基因1 431个,下调表达基因2 226个。GO富集分析发现上调表达基因主要与非生物胁迫响应和化学刺激响应相关,其中,参与水分胁迫响应和脱落酸（abscisic acid,ABA）刺激响应的基因分别有127和141个,而下调表达基因与植物病原菌防御、蛋白激酶活性和水杨酸（salicylic acid,SA）刺激相关。KEGG富集分析表明上调表达基因主要富集于苯丙烷和类胡萝卜素的生物合成及淀粉与蔗糖代谢途径,而下调表达基因主要富集于植物-病原菌互作和植物激素ABA、SA和茉莉酸（jasmonic acid,JA）信号转导途径。qRT-PCR检测6个DEGs的表达模式与RNA-Seq分析结果一致,证实了RNA-Seq结果的可靠性。【结论】RNA-Seq分析鉴定出3 657个甘蓝型油菜叶片干旱胁迫应答基因。GO和KEGG代谢途径分析明确了差异表达基因富集的分子功能与代谢途径。     

天台铁线莲花色变异的转录组分析

天台铁线莲是中国特有的铁线莲属植物,为大花平展型藤本花卉,具有较高观赏价值,花色通常白色.项目组在天台铁线莲引种驯化过程中发现了花萼呈现白色和紫色相间以及花萼全部为蓝紫色的多个花色变异类型.为进一步探究影响花色变异机制,以白色原种天台铁线莲为对照进行了转录组测序.测序结果共获得136.03 Gb的clean data....     

谷子光敏色素家族基因的生物信息学及表达模式分析

光敏色素(PHY)是植物体重要的光感受器,在光形态建成中发挥重要作用.为挖掘新型模式植物谷子的光敏色素基因,通过BLAST从xiaomi基因组中共鉴定出4个SiPHYs,命名为SiPHY-1、SiPHY-2、SiPHY-3、SiPHY-4.利用生物信息学方法对SiPHYs的基因结构、理化性质、系统进化、启动子顺式作用元件和不同时期不同组织的表达谱进行分析,并采用RT-qPCR分析6个品种(系)在不同光周期下SiPHYs差异表达.结果表明,4个SiPHYs分布在8号和9号染色体上,除SiPHY-2外其余3个成员均为酸性亲水性蛋白;系统进化显示,SiPHY-1和SiPHY-2为PHYA亚族,SiPHY-3为PHYB亚族,SiPHY-4为PHYC亚族,且SiPHYs与狗尾草SvPHYs、玉米Zm-PHYs、高粱SbPHYs的亲缘关系近,与拟南芥AtPHYs亲缘关系较远;SiPHYs的顺式作用元件中光响应元件占比最多.表达谱分析结果显示,SiPHY-1可能与谷子的种子萌发有关,SiPHY-3可能与谷子形态建成有关,SiPHY-4可以调控谷子的抽穗开花期;进一步对6个品种(系)不同光周期下的差异表达显示,SiPHY-1在晋谷21和大同红谷中差异显著,SiPHY-3在大同红谷中差异显著,SiPHY-4在晋谷21、xiaomi、黄粟和xiaomi4中差异显著,而Si-PHY-2在不同光周期和不同材料中表达量均很低或不表达,说明SiPHY表达量受光周期和基因型的影响.     

谷子SiARGOS1的克隆、表达分析和功能标记开发

【目的】从谷子中分离受激素诱导表达、参与器官大小控制的拟南芥ARGOS(Auxin-regulated gene involved in organ size)基因家族的同源基因,进行生物信息学分析,明确其在不同组织器官及其受植物激素诱导的表达模式,分析基因编码区及其启动子序列差异,开发功能标记,为谷子产量性状相关基因的改良提供依据。【方法】通过对已有ARGOS蛋白保守结构域进行BLAST,明确谷子ARGOS家族成员数目并进行蛋白序列分析,采用同源克隆方法获得谷子ARGOS家族成员之一——SiARGOS1编码区及其启动子序列,用生物信息学方法分析SiARGOS1启动子的顺式作用元件,通过实时荧光定量PCR分析该基因在谷子各器官中以及不同植物激素条件下的诱导表达模式,利用基因编码区及启动子序列的SNP和插入缺失序列开发分子标记,同时利用85份谷子品种的穗重(panicle weight,PW)、穗粒重(grain weight,GW)和千粒重(thousand-grain weight,TGW)等产量性状数据进行基因型间的差异显著性分析,挖掘用于检测该基因与谷子产量性状相关优异等位变异的功能标记。【结果】获得6个谷子ARGOS家族成员,均具有典型的保守OSR(organ size related)结构域,包含2个跨膜螺旋结构和1个高度保守富含亮氨酸区域,克隆了与拟南芥AtARGOS同源的家族成员之一——SiARGOS1编码区及其启动子序列,该基因位于谷子第8染色体上,开放阅读框为342 bp,无内含子,编码113个氨基酸,启动子区域为2 109 bp,含有与生长素、乙烯、茉莉酸和赤霉素等多种植物激素调控有关的元件。表达分析发现,SiARGOS1在谷子根、茎、叶和穗等器官中均有表达,在根中表达量最高,其次为茎和叶,穗中表达量最低。SiARGOS1对生长素吲哚乙酸(indole-3-acetic acid,IAA)不敏感,但受乙烯利(ethephon,ETH)上调表达。不同基因型谷子SiARGOS1序列分析发现,SiARGOS1编码区151 bp(起始密码子83 bp)处存在1个SNP(C/G),导致该基因第28个氨基酸发生突变(Ala/Gly),据此设计一个CAPS-AccⅡ标记;另外,启动子区存在19个SNP和2个In Del,根据-1 652—-1 651处(TA)_(2/3)和-1 165—-1 163处(TCA)_(1/2)的序列差异分别设计SSR引物AP-1和AP-2。同时,用这些标记对85份谷子品种进行检测,CPAS-AccⅡ和AP-1检测到不同基因型的穗重、穗粒重和千粒重的差异均不显著,而AP-2检测的2种基因型间除千粒重差异不显著外,穗重和穗粒重在2015和2016两年间差异均达到显著水平。【结论】在谷子中发现6个ARGOS家族成员,均具有保守OSR结构域,其中,谷子SiARGOS1开放阅读框为342 bp,无内含子,与拟南芥AtARGOS同源,该基因对生长素不敏感,但受乙烯利上调表达。在该基因启动子区开发的SSR标记AP-2可作为功能标记,用于谷子穗重和穗粒重等产量性状相关优异等位变异的鉴定和筛选。     

不同基因型谷子叶片衰老的研究

以谷子品种‘05-61'和‘金谷3号'为材料,研究不同基因型谷子开花至成熟期叶片的叶绿素含量、酶促防御系统的保护酶SOD、CAT活性和MDA积累量的动态变化规律.结果表明,两个参试品种各功能叶片的叶绿素含量,SOD、CAT等保护性酶活性呈由高到低的变化趋势,MDA含量随着叶片的衰老而升高;与‘金谷3号'谷子相比,‘05-61'叶片防御活性氧毒害的保护酶活性高且下降缓慢,丙二醛含量低且增长缓慢,在籽粒灌浆中后期表现更为明显.研究认为,在旱作农业生产中,延缓叶片衰老,对促进糜子籽粒充实,提高谷子产量和品质具有重要作用.     

谷子GRAS转录因子家族的全基因组鉴定、表达分析及标记开发

为了鉴定谷子GRAS家族基因并揭示SiGRASs响应外源植物激素和逆境胁迫过程的表达规律,本研究采用生物信息学方法对谷子GRAS转录因子家族进行全基因组鉴定,利用实时荧光定量PCR(qRT-PCR)进行4种激素和2种胁迫处理后的表达分析,并根据SiGRAS23序列差异开发分子标记。结果表明,谷子全基因组共包含52个GRAS转录因子基因,大部分编码亲水性蛋白;82.69%的基因编码酸性蛋白,长度为362~734 aa,分子量为39.81~100.09 kDa,等电点为4.85~9.53。系统发育分析将谷子GRAS家族划分为10个亚家族。组织表达量分析表明,各亚族基因具有明显的组织表达特异性,LISCL、DELLA和SHR亚族基因分别在叶、茎和根中有较高的表达量,PAT1和HAM亚族大部分基因为组成型表达,但在叶中的表达量最高。SiGRASs启动子区含有多种植物激素、逆境响应相关的顺式作用元件;qRT-PCR结果显示,SiGRASs在不同激素和非生物胁迫下的表达水平存在较大差异,其中PAT1亚族中Seita.2G369400对6种不同的处理响应最为敏感;少数SiGRASs基因在各组织和各种激素和非生物胁迫处理下,表达量均在极低水平。DELLA亚族中的SiGRAS23在遗传群体AJF₅的双亲矮宁黄和晋谷21号间存在序列差异,据此开发的分子标记D8-1与该群体株高性状紧密连锁。本研究为解析谷子GRAS家族基因参与激素信号转导及逆境胁迫响应的功能研究奠定了基础,SiGRAS23的分子标记可用于今后谷子种质资源株高等位变异的筛选。     

谷子种质资源遗传多样性研究

利用谷子基因组序列开发出SSR引物205对,其中171对可以稳定扩增出目的条带,149对有多态性。利用30个多态性SSR标记对41份谷子种质资源进行遗传多样性分析,共检测出291个等位变异,平均每个位点9.7个。谷子地方品种基因多样性指数平均为0.7907,育成品种基因多样性指数平均为0.7571。对41个谷子品种进行聚类,在遗传相似系数为0.164时聚为6组,与生态型基本一致。     

不同短日照对不同叶龄谷子生长发育的影响

谷子是短日照作物,短日照处理可以减少单株叶片数,促进出叶速度,从而提前抽穗,缩短育种进程。谷子在五叶期时,给予8h短日照光周期处理效果最显著。     

谷子类甜蛋白基因家族的鉴定与密码子偏性分析

类甜蛋白在植物的生长发育和抵御胁迫过程中发挥作用。利用生物信息学方法对谷子类甜蛋白家族基因组成、结构、启动子顺式作用元件、密码子使用偏性等进行分析。结果表明,谷子类甜蛋白家族包含43个成员,分布在9条染色体上,分为3种基因结构类型,其中16个成员仅包含1个外显子,14个包含3个外显子成员的内含子相位均为1-2型。根据系统发育分析归为12个聚类组,聚类组5中的基因主要来自Ⅲ号染色体,聚类组6中的基因主要来自Ⅰ号染色体,其内含子相位多为1-2型,这些基因可能来自同一祖先基因,与进化过程中的染色体重组事件有关。88.4%的基因参与应激反应,多个基因启动子区含有激素和胁迫响应元件,说明该家族基因在植物抵御胁迫过程中发挥作用。多数基因有效密码子数(ENC)较小,密码子的第3位的G+C含量(GC3s)值分布集中,包含10个最优密码子,均以G或C结尾,表明谷子类甜蛋白家族基因密码子使用偏性强,基因表达潜力高,进化过程中主要受自然选择压力影响。     

谷子GATA基因家族的鉴定及表达分析

本研究从谷子基因组水平鉴定出33个SiGATA成员,命名为SiGATA1～SiGA TA33.利用生物信息学方法对谷子SiGATA家族的系统发育、基因结构、染色体分布、保守基序、蛋白三级结构和顺式调控元件进行预测和分析.结果 表明,33个SiGATA转录因子被划分为3组,不均匀的分布在8条染色体上,多数含有CX2CX1...     

谷子萌发期响应干旱胁迫的基因表达谱分析

【目的】谷子是一种耐旱作物,通过二代测序技术获得大量的谷子萌发期响应干旱胁迫的差异基因,进而挖掘谷子萌发期抵御干旱的关键基因及其相关的分子机制。【方法】以晋谷45为材料,谷子萌发期分别用18%PEG-6000干旱胁迫（处理组）和蒸馏水（对照组）处理种子并测定1、10和18 h种子的SOD、POD和CAT活性。SOD活性用氮蓝四唑（NBT）法测定,POD活性用愈创木酚法测定,CAT活性用比色法测定;对萌发10和18 h种子的对照组和处理组构建cDNA文库并进行差异基因表达谱分析;利用Bowtie将reads比对到参考基因组,采用RSEM对bowtie的比对结果进行表达量估计;使用DESeq进行差异表达基因分析;利用NR、Swiss-Prot、KEGG、COG和GO在线数据库对差异基因进行功能注释,挖掘调控谷子萌发的关键基因;利用qRT-PCR验证测序结果的可靠性。【结果】处理组的SOD活性整体比对照组高,而POD活性和CAT活性与之相反;随着萌发时间的变化,SOD活性在不断地增加,但CAT和POD活性逐渐减小。基因表达谱序列与所选参考基因组序列高度一致,基因表达呈现出高度不均一性。通过高通量测序最后获得35 470个基因,以RPKM≥0.01为筛选标准,对照样本中分别筛选出24 030和24 486个表达基因,PEG干旱胁迫处理10和18 h的样本分别筛选出24 019和23 877个表达基因;差异表达基因分析表明,谷子萌发10和18 h分别筛选出456和545个差异基因,其中87和267个上调表达基因,369和278个下调表达基因;GO功能显著性富集分析表明,差异基因主要涉及代谢过程,细胞进程和响应刺激;KEGG富集分析表明,差异基因参与到苯丙烷代谢和植物激素信号转导过程;通过qRT-PCR对5个差异基因在干旱胁迫下种子萌发时的表达分析表明,其表达趋势与表达谱分析结果基本一致。【结论】差异表达基因广泛涉及到糖、蛋白质、核酸等生物大分子代谢、次生代谢和能量代谢等过程;SnRK2和PAL可能在干旱胁迫下调节种子的萌发。     

谷子抗锈病反应相关MYB转录因子的鉴定与表达

【目的】谷子锈病是影响其产量和品质的重要因素之一。鉴定谷子抗锈病相关的MYB转录因子,为谷子抗锈病机理研究奠定基础。【方法】 利用real-time PCR技术,检测9个SiMYBs转录因子基因在（1）谷子孕穗期根、茎、叶和穗4个部位的表达情况,（2）在谷子抗（resistance,R）、感（susceptibility,S）锈病反应120 h内的表达丰度差异,（3）在植株外接水杨酸（salicylic acid,SA）和茉莉酸甲酯（methyl jasmonate,MeJA）后24 h内的表达变化;比较SiMYBs基因在谷子R、S、SA与MeJA 4种反应中的表达模式;选择抗病相关SiMYBs基因进行转录激活活性检测和亚细胞定位。【结果】 SiMYB100在叶部表达量最高,其余8个基因均在根部表达量最高,SiMYB074最显著。5个基因的表达与抗病相关：SiMYB074和SiMYB202受锈菌侵染诱导表达,抗病反应早期的表达量明显高于感病反应;SiMYB041和SiMYB177在抗病反应早期下调表达,后期上升至接种前,而感病反应中保持低水平表达;SiMYB100在接种后的48 h内,抗病、感病反应表达模式相反。在响应外源激素SA和MeJA反应中,SiMYBs基因的表达量均发生不同程度的变化。4个基因（SiMYB074、SiMYB100、SiMYB174和SiMYB202）在R、SA与MeJA反应中的表达模式一致,不同于S反应。5个抗病相关SiMYBs基因具有转录激活活性,其编码蛋白质定位于细胞核中。【结论】 鉴定到SiMYB041、SiMYB074、SiMYB100、SiMYB177和SiMYB202 5个基因的表达与谷子抗锈病相关;SiMYB074与SiMYB100在谷子生长发育和抗病过程中都发挥一定的功能;SiMYB074、SiMYB100、SiMYB174和SiMYB202 4个基因可能通过SA和JA信号途径参与谷子的早期抗病反应。