盐胁迫下比拉底白刺差异表达基因的转录组分析

[目的]为了挖掘比拉底白刺耐盐相关基因,对其盐胁迫下差异表达基因进行筛选分析。[方法]以比拉底白刺幼苗为材料,用200 mmol·L^-1 NaCl对幼苗处理7 d,并对胁迫处理和对照植株叶片进行转录组测序及生物信息学分析。[结果]有效序列组装共得到应答盐胁迫的168463条unigenes和196个差异表达基因。通过差异基因GO和KEGG功能聚类,分别获得64个GO功能小类和25条KEGG通路。进一步基因相互作用网络分析发现,转录调控、氧化还原以及抗逆相关基因在比拉底白刺应答盐胁迫中发挥重要作用,其中,筛选到3个重要的节点基因,分别是热激同源蛋白基因、L型凝集素类受体激酶基因和Win类蛋白基因。[结论]本研究获得了盐胁迫下比拉底白刺的差异表达基因及功能注释信息,有助于理解其耐盐的分子机制,为后续开发耐盐分子标记及通过基因编辑改良植物耐盐特性提供了科学依据。     

基于高通量转录组测序筛选马尾松抗松材线虫病相关基因

[目的]研究马尾松抗松材线虫病相关基因,为解析马尾松抗性机理和高抗马尾松早期选择提供理论依据。[方法]对高抗和易感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有望被开发为鉴定高抗马尾松的分子标记。     

基于转录组分析椰心叶甲啮小蜂复壮和保种种群差异表达基因

为了明确椰心叶甲啮小蜂Tetrastichus brontispae响应种群退化的关键功能基因,通过对比其复壮和保种种群的转录组数据,对差异表达基因进行GO富集分析,并将差异基因在KEGG数据库中进行对比,获得显著性富集的通路注释信息.结果表明:相较于复壮种群,保种种群共有553个差异表达基因,其中206个基因上调表达...     

白蜡虫雄虫真蛹转录组分析

[目的]蛹期对于白蜡虫雄虫发育和交配具有重要意义,但目前缺乏对于白蜡虫真蛹基因转录情况的分析,本研究旨在获得白蜡虫真蛹转录组数据,了解其转录组特征。[方法]采用Illumina Hi Seq 2000进行高通量测序,并进行生物信息学分析;采用荧光定量PCR(qRT-PCR)检测3个热激蛋白基因(hsp)在不同虫态中的表达动态。[结果]共获得63 957条unigene、63 272个开放阅读框(ORF)、9 561个SSR分子标记,unigene平均长度674 bp,共有14 327条unigene获得了注释,Gene Ontology(GO)注释和Kyoto Encyclopedia of Genes and Genomes(KEGG)注释分析表明,很多基因参与的功能与白蜡虫蛹期的生理特征吻合。[结论]该转录组数据为白蜡虫基因功能鉴定和蛋白组学分析研究奠定了数据基础。     

短日照诱导白花泡桐顶芽死亡过程相关基因的表达

【目的】对短日照诱导的白花泡桐顶芽进行转录组测序分析,探究其顶芽死亡的分子机制,为解决泡桐“冠大干低”提供理论依据。【方法】在25℃恒温条件下,对白花泡桐1年生苗进行短日照(SD)处理,在高生长期(SDa)、高生长停止期(SDb)、顶芽死亡发生期(SDc)3个时期进行转录组测序,并利用荧光定量PCR(qRT-PCR)测定11个差异基因的表达量,同时结合转录组信息进行验证。【结果】白花泡桐顶芽3个时期转录组测序共得到57.47 Gb数据,使用DESeq2软件进行3个时期样品间比较,共筛选出差异基因44 397条,其中注释到7大数据库(NR、NT、GO、EggNOG、KEGG、UniProt、Pfam)的基因为37 076条(83.51%)。3个时期的差异基因在KEGG上主要富集在植物激素信号通路,对顶芽死亡发生的SDc时期植物激素信号通路上差异基因的分析得出：脱落酸(ABA)激素信号的响应因子ABF下调表达,乙烯(ETH)激素信号的响应因子ERF1上调表达,油菜素内酯(BR)激素信号上编码木葡聚糖内糖基转移酶的TCH4和细胞周期特异蛋白基因CYCD3上调表达,这些基因与大多数植物芽休眠过...     

麻竹笋转录组测序及苦涩味物质合成基因差异表达分析

[目的]研究避光处理对麻竹笋苦涩味物质合成相关基因表达的影响。[方法]利用Illumina HiSeq~(TM)2500平台对自然生长(CK)和覆土处理(EP)2种类型的麻竹笋进行转录组测序,并对差异表达基因进行分析。[结果]转录组测序共获得36.45 Gb原始数据,经组装去冗余处理得到53 388个Unigene,将所得Unigene比对到Nr、Pfam、COG、Swissprot、GO和KEGG数据库中进行功能注释,发现共有31 462条Unigene与其它物种的基因具有同源性。对CK和EP处理所测得的Unigene进行表达量的比较分析,筛选出1 846个差异基因,其中,在EP处理中上调表达基因998个,下调表达基因848个。由KEGG代谢通路分析发现,差异基因显著地富集在32个代谢途径中,其中,包含与苦涩味物质合成相关的糖酵解途径、苯丙氨酸、酪氨酸和色氨酸合成途径等。进一步研究发现,参与麻竹笋芳香族氨基酸、单宁合成的PFK、ENO、PPY-AT/HPP-AT、LAR等酶基因在EP处理中表达量下降,荧光定量PCR验证结果与测序结果基本一致。[结论]避光处理抑制了苯丙氨酸、酪氨酸、单宁生物合成关键酶基因的表达,可能最终影响麻竹笋苦涩味物质的合成。     

基于马尾松干旱转录组的抗旱功能SSR位点分析

[目的]了解马尾松干旱胁迫转录组序列的功能分布和SSR位点分布特征,并探索与干旱关联的SSR位点。[方法]对马尾松幼苗进行持续干旱胁迫,选取干旱10、15、25 d及正常供水对照的马尾松针叶样品,通过总RNA提取、Illumina测序、raw reads去冗、Trinity拼接获得unigenes。利用Blast比对,对Unigene序列进行GO、KOG、KEGG功能注释及分类;利用Misa软件进行SSR位点批量搜寻,primer 3. 0软件进行规模化SSR引物设计;利用GOSeq(1. 10. 0)、KOBAS(v2. 0. 12)软件对差异表达的含SSR位点Unigene进行GO、KEGG显著性富集分析。[结果]马尾松转录组194 821个Unigene中,101 806个Unigene获得注释,包含64 943个GO功能注释、35 880个KOG功能注释及30 882个KEGG注释。搜寻到6 728个SSR位点,分布于6 367个Unigene中,SSR出现频率为3. 45%;重复类型以单、三、二核苷酸为主,分别占总SSR的35. 82%、33. 03%和25. 22%;重复基序以A/T、AT/AT、AG/CT、AGC/CTG、AAG/CTT为主;基序长度以10 20 bp的短序列SSR为主;基序重复次数以5～10次重复占优势;批量设计13 338对SSR引物。422个含SSR位点Unigene具有差异表达; KEGG富集分析发现,有11个含SSR位点的差异Unigene参与了光合作用、植物激素信号传导及类胡萝卜素合成等3个与干旱响应相关的代谢途径。[结论]从较高质量的马尾松转录组中获得101 806个具有注释的Unigene;从6 367个Unigene中挖掘出6 728个SSR位点; 422个含SSR位点的差异表达Unigene中筛选出11个与干旱关联的SSR功能位点,为抗旱功能基因定位及马尾松抗旱分子机制研究奠定基础。     

光强对无瓣海桑幼苗的生长和光合特性的影响

根据实地测定人工无瓣海桑红树林不同环境下的光照强度,研究了不同光强环境对无瓣海桑幼苗的生长和光合特性的影响.结果表明:随着光照强度的降低,无瓣海桑幼苗的生长和光合特性指标水平均呈下降趋势,且不同光强处理间差异显著.在光照强度为全光照的80％和40％时,幼苗的各项生长指标下降,长势较全光照环境下的幼苗差;在荫蔽初期(2～4个月),弱光环境对幼苗的茎高、叶面积和叶片长宽比的增加有促进作用,随着荫蔽时间的延长,则产生抑制作用,不利于幼苗的生长;当光强达全光照的20％时,幼苗生长极为缓慢.在荫蔽处理180 d时,Sa80、Sa40、Sa20的幼苗死亡率分别为42％、70％、100％;在弱光环境下,无瓣海桑幼苗根、茎、叶的生物量分别下降了90％、80％、91％;随着环境光强水平的下降,无瓣海桑幼苗叶片的叶绿素含量、净光合速率等生理参数均呈下降趋势,各处理间差异显著.试验说明光强不足是造成无瓣海桑幼苗生长和光合水平低下及死亡率较高的主要原因,并对无瓣海桑幼苗的天然更新与自然扩散和人工造林技术进行了讨论.     

漳江口红树林国家级自然保护区无瓣海桑的扩散现状研究

通过对漳江口红树林自然保护区无瓣海桑分布现状进行研究,结果表明:(1)漳江口红树林国家级自然保护区内的无瓣海桑尚处于扩散的早期阶段,个体数量呈逐年增长趋势;(2)中游咸淡水交汇区是无瓣海桑分布较集中的区域,且多沿潮沟生长分布;(3)无瓣海桑和乡土红树植物的生态位重叠;(4)保护区内无瓣海桑幼苗主要来源于自然扩散。     

不同氮形态处理条件下杨树根尖差异表达基因的特征分析

[目的]利用高通量转录组测序技术，在硝态氮或铵态氮处理条件下，对杨树根尖差异表达基因进行了筛选和研究，同时分析和描述了差异表达基因对杨树根尖生长发育的影响，为后续开发高氮素吸收利用效率的杨树新种质提供科学依据。[方法]以灰杨幼苗根尖为材料，用0.5mmol·L^-1硝态氮（NO₃^-）和0.5mmol·L^-1铵态氮（NH₄⁺）对幼苗处理10 d，并对植株根尖进行转录组测序以及生物信息学分析。[结果]硝态氮处理下的主根长度几乎是铵态氮处理条件下的一倍。从两种不同氮形态处理杨树根尖转录组文库中，筛选到2 207个差异表达基因。通过差异基因GO和KEGG功能聚类分析，分别获得50个GO功能聚类和20条KEGG通路。进一步利用MapMan分析，筛选出36个氮代谢通路过程、各类氨基酸的生物合成以及代谢过程相关的差异表达基因。对这些差异基因互作调控网络分析发现，硝酸还原酶（Potri.005G172400）基因通过响应不同氮形态，在影响杨树根尖生长发育过程中发挥了重要作用...     

牡丹花型基因的研究

为研究牡丹花型的分子机理,挖掘调控牡丹雌蕊瓣化的候选基因,揭示牡丹不同花型的分子机理,以牡丹‘黑海撒金’(tree peony ‘HEIHAISAJIN’)和‘璎珞宝珠’(‘YINGLUOBAOZHU’)为材料,采用Illumina Hiseq测序技术进行花瓣的转录组测序,通过序列的拼接、组装和过滤,得到牡丹的转录组;比较‘黑海撒金’和‘璎珞宝珠’中转录本特定基因的表达量,获得差异表达基因;将差异表达基因进行GO功能注释和KEGG代谢通路分类,获得与花发育相关的功能基因。通过de novo的拼接和组装,获得了49 191个unigenes;差异显著性分析得到2 735个差异表达基因(DEGs),其中1 082个为上调基因,1 653个差异表达基因下调。GO功能注释发现,DEGs主要分布在25个功能区,最显著富集的为90 S preribosome。KEGG代谢通路中显著富集的通路有13条,最显著富集的通路为Ribosome。对功能基因进行挖掘和比对,发现与开花相关的ABCDEF模型基因有3个,分别为agamous-like MADS-box protein 65(AGL65), floral homeotic protein APETALA 2(AP2)和EMBRYO SAC DEVELOPMENT ARREST 3(EDA3)。认为利用转录组测序技术,可以分析雌雄蕊发育正常和雌蕊瓣化的不同花型的牡丹花瓣的功能基因。3个ABCDEF模型基因中,AGL65和EDA3为首次在牡丹中发现的花器官发育基因。说明这3个基因很可能是调控牡丹雌蕊瓣化的关键基因,参与牡丹花型的发育。该研究系统地挖掘与花型相关的基因,并发现了新的可能调控牡丹花型基因的成员,不仅为调控牡丹花器官发育挖掘了新的重要功能基因,也为揭示牡丹花型发育的分子机理奠定了重要的基础。     

IBA诱导楸树嫩枝扦插不定根发育的转录组分析

【目的】利用RNA-seq技术对楸树易生根品种‘豫楸1号’的嫩枝扦插生根的4个发育阶段进行转录组测序和分析,为阐明楸树不定根发育的分子机制奠定理论基础。【方法】首先利用浓度为2 g·L-1IBA处理‘豫楸1号’嫩枝,然后分别提取处理后0天(对照)、1天(激活期)、15天(愈伤形成期)的插穗基部的RNA及25天(不定根形成期)和35天(不定根伸长期)的根部的RNA进行转录组测序;接着测序得到raw reads,通过去除接头、重复序列、低质量的序列,得到高质量的clean reads,使用Trinity软件进行拼接获得contigs;随后,利用BLASTx、RSEM、Cytoscape及Me V4.9.0分子生物学软件对测序结果进行注释和差异基因鉴定;最后随机选取15个差异表达的RNA-Seq数据利用q PCR验证基因的表达。【结果】在62 955个unigenes中,31 646(50.26%)个unigenenes获得了注释;差异基因分析发现了11 100个差异基因,其中包括10 200个特异的和900个共同的差异基因;GO富集分析发现‘细胞骨架’仅在激活期显著富集,而‘DNA代谢过程’仅在愈伤组织形成期显著富集;KEGG富集分析显示参与丙烷合成、糖酵解和植物激素代谢等途径的基因对楸树不定根的形成具有重要的作用,并且随着楸树不定根发育进程的发展,参与糖酵解途径的基因数目逐渐减少而参与苯丙素合成的基因数量却在持续增加,表明在IBA刺激后,代谢通路的动态变化响应了不定根的发育进程。【结论】通过对楸树‘豫楸1号’不定根发育的4个阶段的转录组分析,发现细胞分裂素和乙烯间的互作促进楸树愈伤形成而生长素和油菜素内酯间的互作则促进了楸树不定根的伸长。虽然本研究不能完全解释‘豫楸1号’不定根形成的分子机制,但它可以作为进一步探索与这一复杂过程相关的候选途径和基因的有力工具,可为解析楸树不定根发育的分子机制和其他近缘物种的研究提供帮助。     

无瓣海桑人工混交林结构及幼苗天然更新能力

为探究无瓣海桑人工混交林在河口海湾的生态效应问题,采用样方法对福建省厦门市海沧湾16年生无瓣海桑与其它种类红树的人工混交林结构和幼苗天然更新扩散能力进行实地调查。结果表明:①该群落外貌简单、林相整齐、空间层次分层明显。位于乔木一层的是无瓣海桑,位于乔木二层的有无瓣海桑、秋茄、红海榄;位于灌木-小树层的有秋茄、红海榄和白骨壤;位于幼苗层的只有秋茄。②无瓣海桑在乔木层的重要值为35.1%,秋茄在乔木层、灌木-小树层和幼苗层3个不同层次的重要值均为最大。③该群落各多样性指数值均较低,灌木-小树层的各多样性指数值均比乔木层和幼苗层的多样性指数值高。④在4种不同类型样地调查幼苗扩散情况,均未发现无瓣海桑幼苗或历年生的幼树,而乡土红树植物秋茄在4种类型样地中均有幼苗自然生长,白骨壤在光滩中有发现幼苗生长。     

刺槐幼苗对NaCl胁迫的生理生化响应

[目的 ]研究刺槐幼苗对不同浓度NaCl的生理生化响应,为建立刺槐抗逆性生态生理指标体系提供理论依据。[方法 ]以1年生刺槐幼苗为材料,施加不同浓度NaCl(0、50、100 mmol·L~(-1)),分析其生长、生理及水通道蛋白和Na+/H+逆向转运蛋白基因转录表达变化特征。[结果 ]在NaCl胁迫下,刺槐生长被抑制,叶片相对含水量、光合色素含量、光合作用参数降低,水分利用效率和气孔限制值升高;根和叶中MDA含量、抗氧化酶(CAT、APX和GR)活性及渗透调节物质(游离脯氨酸、氨基酸和可溶性蛋白质)含量升高;NaCl胁迫改变了刺槐体内Na~+、K~+、Mg~(2+)的离子平衡状态,但是对Ca~(2+)的含量影响不显著。NaCl胁迫诱导了刺槐根中水通道蛋白基因(TIP1;1、PIP1;1和PIP2;1)、叶中水通道蛋白基因(TIP1;1和PIP1;1)及根和叶中Na+/H+逆向转运蛋白基因NHX1的转录表达。[结论 ]刺槐通过增加叶片水分利用效率、合成并积累抗氧化酶和渗透调节物质及诱导相关基因转录表达,增强了盐适应能力。     

无瓣海桑引入对龙海九龙江口乡土红树林影响的调查

由于龙海九龙江口引种的无瓣海桑在树体高度、生长速度等方面都显著超过乡土红树植物秋茄、桐花树和白骨壤,引起对外来生物入侵的担忧,因此加强对无瓣海桑的生态监测十分必要。通过样方调查,研究了无瓣海桑引入对秋茄、桐花树的影响。结果表明：无瓣海桑林冠下光照良好,为秋茄和桐花树的生长和更新提供了充足的空间。无瓣海桑的引入,提高了红树林群落结构的复杂性和生物多样性。无瓣海桑林冠下秋茄种群密度提高,桐花树种群密度下降。无瓣海桑对林下秋茄的幼苗数量没有影响,对桐花树的幼苗数量略有促进。所有样地内均未发现无瓣海桑的幼树或幼苗,未发现生物入侵迹象。     

汕头市海桑、无瓣海桑冻害调查初报

海桑(Sonneratiacaseolaris(L.)Engl.)、无瓣海桑(S.apetalaBuch-Ham)是汕头市近几年成功引种的2个优良红树林(mangroveforests)树种。报告了汕头市海桑、无瓣海桑的苗木和林木的冻害情况,并就海桑和无瓣海桑的抗寒性及其在汕头市育苗和防寒的技术措施进行讨论。     

偏肿革裥菌在木质环境下转录组构建与相关基因表达分析

【目的】对白腐菌偏肿革裥菌在木质和非木质环境下的转录组进行测序,为白腐菌降解木材的机制研究提供支持。【方法】采用高通量测序技术对木屑和非木屑处理条件下的菌丝样本进行转录组测序。利用eggNOG、GO、KEGG等数据库对转录本进行注释和比较分析,预测和筛选出偏肿革裥菌与木材降解有关的基因。应用荧光定量PCR(qPCR)检测mnp2等11个与木质素降解相关的差异基因在木屑处理5天条件下的表达量,结合转录组数据对这些基因的表达量变化趋势进行分析。【结果】偏肿革裥菌转录组测序共得到38.9 Gb Clean Data,各样品Clean Data平均达到6.49 Gb,各样品的Reads与参考基因组的比对效率在71.23%~74.25%之间。使用edgeR软件进行差异表达分析,得到差异表达基因898个,其中上调351个、下调547个。差异基因被注释到GO数据库的有251个,注释到KEGG数据库的有223个,注释到eggNOG数据库的有704个。eggNOG分析表明,差异基因表达多聚集在能量产生和转换、转录后修饰、蛋白质代谢、伴侣关系、次生代谢物的生物合成、碳水化合物的运输和分解代谢等功能分类下。GO分析表明,显著性富集与频率较高的生物过程是丙酮酸代谢过程、类异戊二烯生物合成过程、天冬氨酸家族氨基酸代谢过程和木质素分解过程。在KEGG通路分析中,差异基因分布到86个不同的生物途径,差异基因显著富集的代谢通路主要为:碳代谢、柠檬酸循环、芳香化合物降解、乙醛酸代谢。qPCR结果表明在木屑处理条件下基因mnp2、mnp3、lip9、lip2、laccase1和nadp的表达量明显上调;mnp10 s、mrp、cyp450-1、GroES1和GroES2的表达量明显下调,qPCR与转录组测序结果在表达量差异倍数上存在较小偏差,但整体趋势一致,可以证明转录组测序结果是正确可靠的。【结论】偏肿革裥菌降解木材与碳代谢、芳香化合物降解等通路密切相关;与木质素分解等生物过程密切相关。根据基因功能注释的结果得到11个与白腐菌降解木质素相关的重要差异表达基因。     

粤东沙质滩涂6种红树林树种造林试验研究

选择无瓣海桑、秋茄、白骨壤、木榄、桐花和拉关木6个红树林树种,在粤东沿海沙质滩涂地带进行造林试验。结果表明:(1)涨退潮风浪大的流动沙质滩涂不适宜种植红树林,避风、涨退潮海浪较小、沙土流动较缓慢的沙质滩涂可以成功造林;(2)种植后14个月,6个红树林树种造林存活率排序为白骨壤(31%)拉关木(28%)无瓣海桑(25%)秋茄(22%)桐花(18%)木榄(14%),平均保存率为23%;(3)在风浪较小、沙土流动较缓慢的沙质滩涂营造红树林,可选择无瓣海桑、白骨壤和拉关木3个树种。     

楸树雄性不育花芽转录组测序及分析

【目的】以自然突变的楸树雄性不育花芽为研究材料,分析与楸树雄性不育相关基因的表达模式,从基因的表达水平上揭示楸树雄性不育的分子机制,为楸树及木本植物的雄性不育研究提供有价值的参考。【方法】采用高通量测序技术对自然状态下楸树的雄性不育的花芽以及可育的花芽进行转录组测序,通过生物信息学对可育花芽与不育花芽的转录本进行比较分析,预测和筛选出与楸树雄性不育有关的基因。【结果】转录组测序共产生27.18 Gb数据,组装并去冗余后得到86 076个Unigene,然后将Unigene比对到7大功能数据库(NR,NT,GO,COG,KEGG,Swissprot,Interpro)进行功能注释,最终被7大数据库中任意一个数据库注释上的Unigene总数为64 600(75.05%)。将试验组(雄性不育花芽,SL)与对照组(可育花芽,FL)所测得的Unigene进行表达量的分析后,筛选出表达量有差异且可信度高的差异表达基因。在噪音分布中(差异表达倍数在2以上,可信度在0.8以上)共筛选出试验组中有6 915个基因上调表达,3 504个基因下调表达。在泊松分布中(差异表达倍数在2以上,错误发生率在0.001以下),SL-1 vs FL-1、SL-2 vs FL-2、SL-3 vs FL-3三个生物学重复中得出差异上调表达基因分别有13 979,13 513,13 055个,差异下调表达基因分别有12 170,13 807,10 411个。差异基因的GO功能分析表明,在生物过程中显著性富集的条目集群频率较高的是生殖过程、生殖发育过程、生殖系统发育、生殖结构发育,在分子功能中只有生长素流出跨膜转运蛋白活性显著性富集。差异基因KEGG通路分析中,差异基因映射到127个不同的生物途径,显著富集代谢通路中差异基因的生物途径主要为:代谢途径、次生代谢产物的生物合成、剪接体、RNA转运以及甘油磷脂及淀粉和蔗糖的代谢。在用差异基因与同源基因比对分析中,共比对出246个同源性高的Unigene,其COG功能分类多聚集在RNA加工和修饰,细胞周期控制、细胞分裂、染色体分区,转录等,同时将这些同源性高的差异表达基因映射到了丙酮酸代谢、植物激素信号转导途径中。【结论】楸树雄性不育的形成与生殖发育的多个过程、丙酮酸代谢途径、生长素流出跨膜转运蛋白活性以及油菜素内酯介导的信号转导通路密切相关。根据分析的结果和已经完成的相关细胞学观察,推测楸树雄性不育的形成可能是楸树体内丙酮酸代谢过程异常,抑制油菜素内酯的合成,使绒毡层发育异常并进一步影响到小孢子减数分裂,最终导致不育花粉的形成。     

翠亨湿地无瓣海桑人工林土壤有机碳分布特征及与土壤理化指标相关性

[目的]通过对比不同潮位带无瓣海桑人工林土壤有机碳的含量和分布规律,探讨与无瓣海桑人工林土壤有机碳分布有关的主要土壤理化性质指标的相关关系。[方法]以中山翠亨国家湿地公园无瓣海桑人工林为例,选择高潮带和中潮带无瓣海桑人工林以及低潮带光滩区域,测定土壤有机碳含量及主要土壤理化指标,并计算土壤有机碳密度,再进一步分析其与土壤理化指标的相关性。[结果]研究表明:在垂直方向上,土壤有机碳平均含量和平均密度均表现出显著差异,从大到小依次为:高潮带〉中潮带〉低潮带;其中,在高潮带,土壤有机碳含量最大值出现在0~20 cm土层,有机碳密度最大值出现在20~40 cm土层。土壤有机碳含量与含盐量显著负相关,与土壤全氮、全磷极显著正相关,与土壤pH值极显著负相关。[结论]高潮带无瓣海桑人工林土壤有机碳的密度和含量均高于中潮带和低潮带;土壤全氮、全磷含量、pH值与土壤有机碳含量关系最密切,这些指标可以用来预测无瓣海桑人工林的土壤有机碳分布。