避光对麻竹笋苦涩味及单宁含量、形态与分布的影响

[目的]探讨不同避光处理对麻竹笋苦涩味强度及单宁含量、形态与分布影响的机理。[方法]运用感官评定、光学显微镜、电子显微镜透射技术和磷钨酸-钨酸钠比色法,采用覆土(EP)、双层不透光套袋(CLPB)和自然生长(CK)3种不同处理,对麻竹鲜笋的口感品质进行感官评定,分析测定麻竹笋笋体单宁物质含量、形态和分布。[结果](1)不同处理的麻竹鲜笋的口感不仅均呈苦涩味,且苦涩味强度均表现为由麻竹笋的基部到笋尖逐渐增强的趋势,其苦涩味强度整体大小顺序为:CKCLPBEP。(2)麻竹竹笋壁中含有单宁物质的细胞(CWT)可被2%氯化亚铁溶液染成黑色,与不含单宁物质的细胞区别明显。(3)单宁大量分布在薄壁细胞内,少量在纤维细胞内,维管束中的筛管和导管细胞中无单宁分布。(4)单宁在薄壁细胞中主要分布于细胞质内,少量分布在液泡中,电镜下CWT中积累的单宁可分为絮状、颗粒状和块状3种类型。(5)不同处理的麻竹笋单宁含量为CK:1.15 2.67 mg·g~(-1),CLPB:1.03 1.43 mg·g~(-1),EP:0.36 1.13 mg·g~(-1),其相应的CWT密度大小顺序为:CKCLPBEP,3种处理不同部位的单宁含量与CWT密度大小均为笋尖中部基部。[结论]避光显著降低麻竹竹笋壁的苦涩味强度及单宁含量,自然光照下生长的麻竹笋苦涩味强度及单宁含量显著高于避光处理;光与竹笋中单宁物质分布和形态相关不显著,单宁形态仅与竹笋部位密切相关。覆土、套袋等避光措施可显著降低笋体单宁物质的含量,降低苦涩味,改善竹笋的口感品质,为培育低苦涩味麻竹笋提供科学的理论依据。     

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

【目的】以自然突变的楸树雄性不育花芽为研究材料,分析与楸树雄性不育相关基因的表达模式,从基因的表达水平上揭示楸树雄性不育的分子机制,为楸树及木本植物的雄性不育研究提供有价值的参考。【方法】采用高通量测序技术对自然状态下楸树的雄性不育的花芽以及可育的花芽进行转录组测序,通过生物信息学对可育花芽与不育花芽的转录本进行比较分析,预测和筛选出与楸树雄性不育有关的基因。【结果】转录组测序共产生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.5mmol·L^-1硝态氮（NO₃^-）和0.5mmol·L^-1铵态氮（NH₄⁺）对幼苗处理10 d，并对植株根尖进行转录组测序以及生物信息学分析。[结果]硝态氮处理下的主根长度几乎是铵态氮处理条件下的一倍。从两种不同氮形态处理杨树根尖转录组文库中，筛选到2 207个差异表达基因。通过差异基因GO和KEGG功能聚类分析，分别获得50个GO功能聚类和20条KEGG通路。进一步利用MapMan分析，筛选出36个氮代谢通路过程、各类氨基酸的生物合成以及代谢过程相关的差异表达基因。对这些差异基因互作调控网络分析发现，硝酸还原酶（Potri.005G172400）基因通过响应不同氮形态，在影响杨树根尖生长发育过程中发挥了重要作用...     

基于马尾松干旱转录组的抗旱功能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功能位点,为抗旱功能基因定位及马尾松抗旱分子机制研究奠定基础。     

赤霉素GA_(4+7)处理下白桦无性系生长及差异基因表达分析

[目的]为研究赤霉素GA4+7对白桦生长的调控机制。[方法]试验利用1个白桦无性系的分株材料,连续2年分别叶面喷施浓度为0.2 g·L~(-1)(T1)、0.4 g·L~(-1)(T2)的GA4+7溶液,对实验苗木的生长情况进行跟踪调查,并对相应顶端生长组织进行转录组测序和差异基因分析。[结果]对2个年份的苗高、地径等性状分析发现,不同处理间苗高与地径的差异均达到显著或极显著水平,GA4+7对高生长有明显的促进作用,2年后的T2白桦苗高较CK提高了19.65%。于6月28日取GA4+7连续处理2年的白桦顶芽开展RNA-seq分析,结果显示:0.4 g·L~(-1)GA4+7(T2)处理与对照(CK)间的差异基因数量最多,即上调表达的基因有181个,下调表达的基因有55个,这些差异基因在叶绿体类囊体膜、红光、远红光及蓝光的细胞响应等方面富集明显。Pathway富集分析显示:氧化磷酸化途径与光合途径中的8条基因呈上调表达,认为施加0.4 g·L~(-1)GA4+7处理后,上述基因的上调表达加快电子传递进程、催化ATP的合成,从而促进白桦光合作用增强,提高苗期生长量。[结论]对白桦外源施加赤霉素GA4+7后,可上调氧化磷酸化途径和光合途径相关基因的表达,进而促进苗木的高和地径生长。研究结果可为白桦生长相关途径基因的克隆提供参考。     

蕨类植物芒萁幼孢子体转录组高通量测序及特征分析

[目的]采用高通量测序技术Illumina Mi Seq250获得蕨类植物芒萁的孢子体转录组数据,以期为芒萁的生长、发育、代谢调控、微进化机制分析等提供重要的分子信息。[方法]应用生物信息学方法对测序获得的大量单基因簇(Unigene)进行基因功能注释、代谢途径及微卫星分析等。[结果]本研究共获得18 463 296条序列读取片段(reads),总碱基数为4.62Gbp序列信息,经序列组装最终得到63 169个Unigene,平均单条Unigene长度为863 bp,N50为1 587 bp,其中分布在200 500 bp长度区间Unigene占总数的55.4%。数据库中的序列同源性比较表明,26 826个Unigene与其他物种的已知基因具有不同程度的同源性。芒萁转录组中的Unigene根据GO功能大致可分为细胞组成、分子功能和生物学过程3大类47个分支,其中有大量的Unigene与细胞进程、绑定活性、代谢过程和催化活性相关。将Unigene与COG数据库进行比对,根据其功能大致可分为26类。以KEGG数据库作为参考,依据代谢途径可将Unigene定位到257个代谢途径分支。此外,利用MISA软件检索2 6碱基微卫星,共找到13 286个SSR。在不同长度微卫星中,三核苷重复数量最多,占总数的40.41%。在各重复基序类型中出现频率最高的为AG/CT(14.45%)与AAG/CTT(12.39%)。利用重复基序开发的多态性SSR标记,可应用于芒萁不同个体的基因型分型鉴定。[结论]本研究获得了较高质量的芒萁转录组数据库,揭示了芒萁孢子体生长发育过程中表达基因的功能总体特征,可为芒萁进一步的功能基因挖掘和分子标记规模化开发奠定基础。     

无瓣海桑根响应盐胁迫的转录组分析

[目的]初步探究无瓣海桑的耐盐分子机制,筛选出无瓣海桑抗盐候选基因,为后续功能验证实验及林木抗盐性遗传育种奠定分子基础。[方法]以1年生无瓣海桑幼苗为材料,用500 mmol·L^-1 NaCl分别处理0 d(对照组)和10 d(处理组),取不同条件下的根部组织进行转录组测序,并结合三代全长转录组数据进行后续生物信息学分析。[结果](1)与对照组相比,NaCl处理10 d后,无瓣海桑幼苗根系中共有14 401个差异表达基因,其中,7 153个上调,7 248个下调。(2)GO分析发现,共有11 068个差异基因在47个GO条目得到注释。(3)在KEGG富集分析中,共有6 189个差异基因富集到134条通路,其中,共有14条通路显著富集(P值<0.01,Q值<0.05)。(4)通过进一步对差异基因进行功能注释分析,共筛选出抗盐候选基因89个,其中,抗氧化基因24个,渗透调节物质基因22个,植物激素基因19个,蛋白激酶基因10个,转录因子基因14个。[结论]活性氧清除、渗透调节、植物激素、蛋白激酶及转录因子相关基因参与调控无瓣海桑盐逆境适应过程。     

基于中麻黄萌发种子转录组的黄酮类化合物合成途径基因的挖掘

应用新一代高通量测序技术(Illumina Solexa),对中麻黄(Ephedra intermedin)同萌期的种子进行转录组测序,数据重头(denovo)组装后,共获得52 007条Unigene序列,总长为25 043 596 bp。COG功能注释、GO分类及KEGG代谢通路分析后,获得了64 751个GO功能注释,17 701个COG功能注释以及16 942个KEGG注释,并从KEGG通路中找到参与黄酮类化合物合成途径的关键基因片段283个,其中包含了查尔酮合酶基因、细胞色素相关基因和花青素还原酶等基因。中麻黄转录组测序工作的完成,极大地扩充了麻黄的基因资源,为麻黄属植物分子系统进化分析、抗性和药用功能基因的发掘与利用、遗传改良等研究奠定基础。     

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

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

普通油茶花发育过程数字表达谱测序及分析

为了解普通油茶花发育过程的分子机制及成花关键基因,选择花芽分化前、萼片形成期、花瓣形成期、雌雄蕊形成期、子房及花药形成期和雌雄蕊成熟期的花芽作为试验材料,对此6个时期样品进行表达谱测序并进行分析。结果表明,测序共获得原始reads数为103 249 386条,平均Cycle Q20均为100%;各样品测序reads与转录组部分测序构建的Unigene库的序列对比,对比效率超过60%,共获得差异表达基因26 861个,7月3日与6月15日差异基因达到19 596个,上调基因数19 408个;对差异基因进行生物功能分析、翻译,核糖体结构和生物合成为245条,转录因子为372条,复制、重组和修复为320条,细胞信号转导机制为285条,翻译后修饰、蛋白质周转和伴侣蛋白为283条,能源生产和转换为114条,碳水化合物运输和代谢为136条,氨基酸的转运和代谢为122条;与花器官形成有关A类基因Unigene数量最多,其中Unigene67783在不同的发育阶段差异较明显,表达水平也较高;B类基因有9条Unigene,其中Unigene56059表达水平较高,各发育阶段差异也较明显;属于C类基因的Unigene数量有9条,总体表达水平较低且差异不明显。     

重金属胁迫下白骨壤数字基因表达谱分析

[目的]从转录组水平研究重金属污水胁迫处理下白骨壤的基因表达变化,以揭示白骨壤响应重金属污染胁迫的机制。[方法]采用基于高通量测序的数字基因表达谱(DGE)技术对重金属污水胁迫处理组和对照组白骨壤材料进行转录组测序分析。[结果]与对照组相比,重金属处理下白骨壤共有10 459个差异表达基因,其中,8 685个表达量上升,1 774个表达量下降。Me V聚类表明:大部分基因在重金属处理样本中的表达量受到极大的诱导。GO功能注释发现,差异表达基因主要定位于叶绿体以及细胞内各种膜结构,参与"细胞代谢"、"细胞组分的组合和生物合成"、"对刺激的响应"等过程。KEGG代谢通路分析显示,差异表达基因分布于126条Pathways,涉及"核糖体"、"光合作用"、"乙醛酸盐代谢"、"丙酮酸代谢"等途径。重金属胁迫还促进萜类、黄酮类化合物生物合成的关键酶基因(牻牛儿基焦磷酸合酶(0009238)、黄酮醇合成酶(0001833))的表达,进而促进白骨壤有效成分的积累;显著诱导细胞分裂素水解酶编码基因CKX7(0057124)、油菜素内脂信号转导组分基因BSK(0043741)等的表达,进而提高白骨壤对重金属胁迫的适应能力。此外,转录因子分析发现,b HLH、NAC、MYB-related和WRKY在重金属胁迫中发挥重要作用。实验选取8个与环境刺激响应密切相关的基因,通过qRT-PCR验证了它们在重金属污水胁迫处理下的表达差异,结果与数字基因表达谱分析的结果较一致,证实了差异表达基因数据的有效性。[结论]重金属处理影响了白骨壤大量的新陈代谢、光合作用、小分子酸合成、激素合成与信号转导及次生代谢产物合成相关基因的表达。     

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

[目的]蛹期对于白蜡虫雄虫发育和交配具有重要意义,但目前缺乏对于白蜡虫真蛹基因转录情况的分析,本研究旨在获得白蜡虫真蛹转录组数据,了解其转录组特征。[方法]采用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)注释分析表明,很多基因参与的功能与白蜡虫蛹期的生理特征吻合。[结论]该转录组数据为白蜡虫基因功能鉴定和蛋白组学分析研究奠定了数据基础。     

云南干热河谷地区余甘子转录组分析

[目的]对云南干热河谷地区余甘子转录组特征进行描述,旨在为余甘子微卫星标记的开发和功能基因的挖掘提供较全面的背景信息。[方法]采用Illumina Hiseq 4000测序平台对余甘子叶片进行转录组测序,对原始数据进行过滤、de novo组装及聚类去冗余等处理后,再与公共数据库进行比对,对Unigenes进行基本功能注释、CDS预测、TF编码能力预测及R-Gene预测等分析。[结果]本研究共获得10. 52 Gb的Clean reads,Q20、Q30分别为98. 47%、95. 28%。组装并去冗余后获得76 881条Unigenes,平均长度、N50分别为713、1 257 nt。通过与NR、COG、KEGG和Swiss Prot数据库进行比对,44 768条Unigenes获得功能注释。余甘子转录组Unigenes根据COG功能注释信息大致分为25类;按GO功能注释信息划分为生物学过程、细胞组分和分子功能3大类47亚类;参考KEGG注释信息,可归为6大代谢通路、21类代谢途径,其中约3/5为代谢相关通路。根据以上注释结果共检测出42 953个CDS,其余未比对上的Unigenes用ESTScan预测后得到2 058个CDS。同时,预测到56个TF家族以及18种RGene。[结论]本研究获得的余甘子转录组Unigenes序列的组装质量较高、完整性较好、基因丰富、功能多样,极大地扩充了余甘子基因信息库,为今后余甘子乃至叶下珠属植物功能基因挖掘、抗性机理分析、分子标记开发、分子辅助育种等研究提供了重要的基础数据。     

毛竹笋呈味物质种类、含量与辛辣味强度的关系

[目的]对3种类型毛竹笋进行品尝和分析,探索毛竹笋辛辣味的呈味物质,为提高毛竹笋食用品质和经济价值提供参考。[方法]试验运用感官评定、气质联用(GC-MS)、高效液相色谱(HPLC)等方法,对出土毛竹笋(AS)、未出土毛竹笋(US)和毛竹鞭笋(RS)3种新鲜竹笋及AS竹笋水煮后的笋渣、笋汤进行感官评定、成分分析,开展呈味物质含量与相应滋味强度的相关性分析。[结果](1)经感官评定,AS和US 2种类型毛竹笋辛辣味均表现为从笋基部到笋尖依次增强的趋势,而AS竹笋从笋基部到笋尖辛辣味强度值为2.0 9.5,明显高于US的1 6,RS竹笋滋味则表现为微甜,笋基部到笋尖无明显变化。AS竹笋水煮后笋渣滋味为辛辣味,且辛辣味强度随水煮时间的延长而降低;笋汤滋味则表现为苦涩味,滋味强度随水煮时间的延长而增强。(2)通过GC-MS方法,共检测出毛竹鲜笋43种成分,其中,酯类9种、醇类6种、酸类5种、酚类3种、酮类2种、醛类3种、烃类4种和其它11种,其中,相对含量较高且与辛辣味有关的物质为草酸、单宁和麻黄碱。(3)相关性分析表明,与AS和US竹笋辛辣味显著相关的呈味物质为氰化物、氰苷、单宁和麻黄碱,与AS竹笋辛辣味强度的相关系数分别为0.812、0.850、0.939、0.818(P0.01);与US竹笋辛辣味强度的相关系数分别为0.905、0.857、0.939、0.937(P0.01)。AS水煮后笋渣中氰化物、氰苷和单宁与辛辣味强度显著相关,相关系数分别为0.859、0.861、0.933(P0.01),笋汤中与苦涩味强度显著相关的有氰化物、单宁和草酸,相关系数分别为0.982、0.954、0.976(P0.01)。综合表明:毛竹笋辛辣味的呈味物质是氰化物、氰苷和单宁。[结论]影响毛竹笋辛辣味的主要呈味物质为氰化物、氰苷和单宁,其辛辣味是三者相互作用的结果。     

泡泡刺高通量转录组鉴定及其黄酮类代谢途径初步分析

[目的]为了更好的了解泡泡刺黄酮类生物合成途径及其潜在的抗逆作用,[方法]应用高通量RNA-seq测序技术对泡泡刺当年生新叶进行了转录组测序及相关生物信息学分析。[结果]表明:泡泡刺当年生叶中获得了13 013 444 total reads,总核苷酸数为1 171 209 960 bp(1.09 Gb),组装拼接后得到48 921条unigene序列;Nr、Swiss-Prot、COG、GO、KEGG等数据库分析显示,有30 407个Nr注释,19 671个Swissprot注释,9 273个COG功能注释,46 153个GO功能注释,13 654个KEGG注释,并从KEGG通路中找到参与黄酮类化合物合成途径的关键基因片段186个。[结论]本研究获得了较好的泡泡刺转录组序列信息,并且其中富含黄酮类化合物代谢途径的相关基因,这为其抗逆机制、药用价值的研究及开发利用提供宝贵资源。     

杨树PtROP家族基因的表达分析与功能预测

[目的]为探究ROP基因在林木中的功能。[方法]本研究以杨树为模式,在全基因组水平上对ROP基因的家族成员、基因结构、保守结构域、氨基酸序列相似性及表达模式进行了分析。[结果]显示,杨树PtROP基因家族包含13个成员,不同成员间在进化上相对保守,均存在与GTP结合与水解相关的结构域。这些基因在不同组织、不同胁迫条件下的表达具有明显差异,说明它们参与不同的生物学过程。通过构建功能基因网络(functional gene network for poplar)发现PtROP基因主要参与了信号转导过程。[结论]杨树PtROP基因家族包含13个成员参与不同的生物学过程,可能主要参与了信号转导过程。     

11条白桦BpSPL家族基因的生物信息学和表达分析

[目的]研究白桦中SPL转录因子的基因序列特征及其在不同时期不同组织中的表达规律。[方法]依据白桦45个转录组测序结果,共获得12条全长SPL基因,依次命名为白桦BpSPL1-BpSPL12,对其中11条BpSPL进行了生物信息学及基因表达特征分析。[结果]生物信息学分析发现,11条BpSPL均含有1个高度保守的SBP结构域,且基因长度差异较大,含有2 10个不等数目的外显子,系统进化分析发现11条白桦SPL分属于六大类SPL蛋白;qRT-PCR分析结果显示,白桦BpSPL基因在不同时期的叶、顶芽、茎、雄花序中表达变化显著,多数SPL基因在顶芽中7月5日和8月20日至9月20日时期表达较高,除BpSPL1基因外在雄花序从6月份到9月份的生长阶段中,SPL的表达水平呈现逐渐升高的总体趋势。[结论]BpSPL基因可能参与到了白桦顶芽和雄花序的生长发育过程。