miR-181a-5p靶基因预测及其在鹅卵泡颗粒层中的表达分析

为深入探究miR-181a-5p在鹅颗粒细胞中可能的功能及其作用机制奠定基础。以多个网络数据库资源为基础,通过生物信息学的方法预测可能的靶基因及其信号通路,然后利用q PCR检测miR-181a-5p及其靶基因在不同阶段颗粒层中的表达水平。结果显示,miR-181a-5p在脊椎动物中高度保守;靶基因的GO和KEGG分析发现,miR-181a-5p主要参与调控细胞的增殖、分化、凋亡等过程,多数靶基因主要富集到FOXO、MAPK、PI3K-Akt等信号通路;根据靶基因互作网络图,筛选出10个核心靶基因VEGFA、MAPK1、FOS、KRAS、HRAS、SIRT1、BCL2、ESR1、PTEN和CDKN1A,其中多数靶基因均与细胞增殖凋亡相关。结果表明,miR-181a-5p在2～4 mm表达量最高,但在4～6 mm急剧下降,之后在4～6 mm、8～10 mm、F5这3个阶段均呈现显著上升趋势。miR-181a-5p可能通过靶向调控SIRT1、BCL2、CDKN1B、PTEN等基因,参与FOXO、MAPK、PI3K-Akt等信号通路来调控颗粒细胞的增殖凋亡过程。     

亚砷酸盐胁迫下水稻microRNA的应答功能及调控通路分析

为明确砷胁迫下microRNA的应答功能及调控通路,本研究利用miRNA芯片检测水稻幼苗在亚砷酸盐胁迫后miRNA的表达差异。结果发现11个miRNA表达显著下调,7个显著上调,并采用RT-PCR验证了9个miRNA的差异表达。差异miRNA靶基因多编码转录因子以及重金属响应网络中的蛋白。基因本体(GO)和KEGG通路富集表明,砷胁迫下miRNA主要介导调控激素信号转导、硫代谢、氨基酸代谢和抗氧化防御。利用PLACE数据库发现miRNA上游启动子区密集分布多种金属应答元件。本研究有助于阐明水稻中砷转运和积累的分子机制,为水稻抗逆育种提供理论参考。     

杨树响应细菌性溃疡病菌侵染的miRNA靶基因预测及分析

欧美杨细菌性溃疡病是由Lonsdalea quercina subsp. Populi引起的枝干病害,给杨树的生长和生存带来极大危害。miRNA是植物响应逆境胁迫的重要调控分子,为鉴定杨树中响应欧美杨细菌性溃疡病菌侵染的miRNA靶基因,以期为林木抗病分子育种提供新的候选基因,本研究以miRBase数据库中杨树全部miRNA成熟序列为探针,对‘中林46’杨接种细菌性溃疡病后的差异表达基因进行了miRNA的靶基因预测。共筛选出杨树127个miRNA家族的276个miRNA对应的566个靶基因,涉及植物病原互作、植物激素信号传导和苯丙素生物合成等多个途径,表明miRNA参与了杨树对细菌性溃疡病菌侵染的响应。结合靶基因的GO分析、KEGG通路富集和基因功能注释,发现miR482、miR6459、miR7812、miR7835和miR396等通过靶定RPS2、NBS-LRR、AUX1、CCR、Ca2+跨膜运输蛋白编码基因,参与调控杨树对病菌侵染的响应。miR7835靶定的9个差异表达基因中有6个基因均编码苯丙素生物合成途径通路中的CCR酶,推测其可能影响杨树木质素合成,进而参与杨树对病菌侵染的防御反应。本研究发现的miRNA:mRNA调控基因对,丰富了林木抗病分子育种的基因资源,为林木抗病分子调控机制的研究提供了新的线索。     

荣昌猪、长白猪先天免疫候选基因筛选

旨在探究影响荣昌猪和长白猪机体先天免疫功能的候选基因.以新生荣昌猪、长白猪胸腺为研究对象,利用RNA-seq分析荣昌猪、长白猪胸腺组织中的差异表达基因,同时对差异表达基因进行GO和KEGG富集分析,并通过qRT-PCR验证RNA-seq结果的准确性.结果表明:在新生荣昌猪和长白猪胸腺中,差异倍数在2倍以上的基因共421个,其中,荣昌猪相对长白猪上调的有231个,下调的有190个.GO富集分析发现,这些差异表达基因主要富集在清道夫受体活性、G蛋白偶联受体结合、T细胞活化、细胞分泌调节、胞吐作用调节、SNARE复合物组装调节等生物学过程.KEGG通路分析发现,这些差异表达基因主要参与细胞因子-细胞因子受体互作、cAMP信号、代谢(磷脂酶D、半胱氨酸和蛋氨酸)及一些疾病(非洲锥形病、癌症)相关信号通路.6个差异表达基因的RT-qPCR结果表明,此次RNA-seq的结果准确可靠.综合差异表达基因、GO和KEGG信号通路结果发现,TMSB15A和STAB2、SAA1和HP分别是荣昌猪和长白猪胸腺先天免疫的候选基因.     

辽宁绒山羊皮肤毛囊mir-1298-5p靶基因预测及表达载体构建

探讨辽宁绒山羊皮肤毛囊不同发育时期mir-1298-5p的调控作用及其与靶基因的潜在关系,为mir-1298-5p与其靶基因对皮肤毛囊发育作用提供理论依据。以辽宁绒山羊皮肤毛囊组织为材料,通过miRNA的分离、表达、靶基因预测筛选及靶位点3'UTR表达载体构建,采用RT-PCR进行表达检测,在线靶基因预测软件预测mir-1298-5p的靶基因,并对其靶基因TGF-βR1的靶位点进行克隆测序、缺失突变。结果表明：靶基因预测发现TGF-β R1的3'UTR存在mir-1298-5p种子区结合位点,且mir-1298-5p及其靶基因TGF-β R1均在皮肤毛囊不同发育时期呈现差异性表达。mir-1298-5p对皮肤毛囊周期性发育可能起到一定的调控作用,并成功构建了TGF-βR1 3'UTR表达载体,为后期过表达转染体系的建立和功能基因的验证奠定了基础。     

基于转录组测序筛选黄牛低氧适应性相关差异基因

为探明藏黄牛对低氧适应的相关差异表达基因并分析差异表达基因功能和代谢通路,以藏黄牛和三江黄牛心脏组织为研究对象,提取总RNA,建立测序文库后利用Illumina HiSeqTM 4000进行测序,经生物信息学分析筛选出差异表达转录本,并对差异表达基因进行GO富集和KEGG通路分析,最后采用荧光定量PCR(RT-qPCR)验证测序数据的准确性。结果表明:藏黄牛相比于三江黄牛得到显著差异表达基因608个,其中上调基因467个,下调基因141个,推测表达量较高的B2M、COX7C、NFATC1、ECE1、CAST基因是与低氧适应性相关的基因;差异转录本GO功能富集分析可知,在分子功能、细胞组分、生物过程分别富集到327,253,2 191个条目,其中富集程度最高的条目分别是细胞发育过程、结合与转录因子活性、细胞前缘。KEGG分析结果表明,差异转录本显著富集在8条通路,筛选低氧适应性候选基因NFATC1是MP-PKG信号通路与T细胞受体信号通路成员,推测其通过cGMP-PKG通路与T细胞受体信号通路参与低氧适应性调控。采用高通量测序技术对藏黄牛和三江黄牛心脏组织进行转录组测序分析,获得大量差异表达基因,为研究动物低氧适应性机制奠定基础。     

基于高通量测序分析中华辣椒果实发育期相关的miRNA

辣椒的成熟过程是一个受基因调控表达以及多种通路协作完成的复杂的过程,其中microRNA在调控基因表达方面具有重要的作用。本研究为了进一步探索辣椒果实中miRNA及其靶基因在果实成熟过程中的功能,构建了中华辣椒(Capsicum chinense)绿熟期、转色期、红熟期不同成熟时期果实的小RNA文库,并利用生物信息学分析探讨其调控过程中的信号通路。结果表明,三个时期分别获得22 138 507、17 692 000和14 480 729条miRNA序列,共鉴定出149条保守的miRNA。其中,绿熟期到转色期下调表达和上调表达的miRNA分别有3和7条,绿熟期到红熟期下调表达和上调表达的miRNA分别有26和22条,说明microRNA参与了果实成熟的过程。靶基因预测及功能富集分析显示它们可能参与了乙烯、生长素、细胞分裂素、脱落酸、茉莉酸介导的信号途径。本研究结果为今后辣椒成熟中的激素调节变化的研究提供了可能的理论依据。     

ssc-miR-486靶基因预测及生物信息学分析

通过对ssc-miR-486进行靶基因预测及生物信息学分析,以探索其影响猪脂肪沉积的作用机制。首先利用实时荧光定量PCR技术对脂肪沉积能力相差较大的莱芜猪和大白猪皮下脂肪组织中ssc-miR-486表达情况进行比较研究,然后利用miRBase、Ensembl、NCBI等数据库及miRanda、PITA、RNAhybrid、Cytoscape、JASPAR、Promoter Scan等软件对ssc-miR-486进行生物信息学分析,包括保守性分析、转录因子结合位点预测、靶基因预测、GO富集分析和KEGG信号通路富集分析。结果表明,ssc-miR-486在脂肪沉积能力较强的莱芜猪皮下脂肪组织中发生显著上调。生物信息学分析结果发现,miR-486在各物种间非常保守,ssc-miR-486启动子区域有SREBP1、SREBP2等多个转录因子结合位点,获得的229个靶基因显著富集于脂质代谢、脂质生物合成、细胞代谢等生物学过程及类固醇生物合成、磷脂肌醇代谢、PPAR等信号转导通路中。结果表明,ssc-miR-486参与了猪脂肪沉积或脂肪代谢的调控,为今后研究脂肪沉积提供基础理论支撑。     

基于转录组测序筛选鸡皮肤毛囊性状相关基因和关键通路

旨在通过对崇仁麻鸡背部皮肤毛囊组织的转录组学分析，筛选关键的差异表达基因和信号通路，探讨鸡皮肤毛囊性状的分子调控机制。采用TRIzol法提取皮肤毛囊组织总RNA,通过Illumina平台进行转录组高通量测序(RNA-Seq)。最后采用实时荧光定量的方法对转录组测序进行验证。研究结果如下：共筛选出了3 856个差异表达基因，在公鸡(BGB vs YGB)比较组合中筛选出了971个差异表达基因，其中274个表达上调，697个表达下调。在母鸡(BMB vs YMB)比较组合中筛选出了3 529个差异表达基因，其中1 477个上调，2 052个下调。在公鸡、母鸡共同的差异表达基因中筛选出了3个与皮肤毛囊性状有关的高表达基因KRT75、KRT6A和KRT14。在GO功能显著富集中，BGB vs YGB比较组合富集到了516个GO term, BMB vs YMB比较组合富集到了1 020个GO term,基因主要集中参与细胞黏附等活动。显著富集了多条KEGG通路，在BGB vs YGB比较组合中富集了8条显著的KEGG通路，在BMB vs YMB比较组合中富集了20条显著的KEGG通路，筛选出了...     

2种黑籽南瓜响应枯萎病菌侵染的转录组学研究

为了筛选白皮黑籽南瓜和绿皮花纹黑籽南瓜抗枯萎病的差异表达基因,挖掘参与抗性相关的代谢通路。以尖孢镰刀菌侵染2 d后和未侵染的2种幼苗叶片作为材料,利用Illumina Hiseq 2500测序技术进行转录组测序分析。以未侵染的幼苗作为对照,白皮和绿皮花纹黑籽南瓜差异表达基因分别有2082和1116个。其中141个基因在2种黑籽南瓜中均呈现出差异表达,62个基因具有相同表达趋势。对差异表达基因进行GO功能注释,结果显示2种黑籽南瓜差异表达基因主要富集在生物学过程中。KEGG代谢通路分析发现白皮有22个差异表达基因显著富集在植物激素信号通路,绿皮花纹有30个差异表达基因富集在苯丙烷生物合成通路。本研究挖掘出黑籽南瓜抗枯萎病差异表达基因,预测了抗性通路,为瓜类抗病机制的深入研究奠定了理论基础。     

初步探究LncRNA在甘蓝型油菜生态型分化中的作用

甘蓝型油菜是我国重要的油料作物之一,根据其成花转变过程中对低温春化时间需求的不同分为3种生态型:春性、半冬性和冬性。前人研究发现,长链非编码RNA (LncRNA)可在多层面上调控基因的表达,参与对植物生长发育的调控。在拟南芥中, LncRNA可以通过调控春化途径相关基因的表达来影响开花。本研究以3种生态型甘蓝型油菜为材料,利用高通量测序技术进行苗期叶片的mRNA和LncRNA测序,初步探究LncRNA在油菜生态型分化及适应性形成中的作用。3种生态型甘蓝型油菜共差异表达基因的GO及KEGG富集分析表明,不同生态型甘蓝型油菜之间存在大量基础化合物合成代谢差异,特别是脂质类化合物。3种生态型甘蓝型油菜中共鉴定获得3775个LncRNA,其中285个在2个及以上生态型组合中存在差异表达,涉及到1517个候选靶基因。这些差异表达LncRNA涉及到的靶基因也富集到大量基础化合物合成代谢途径。通过mRNA-LncRNA联合分析,我们预测到了一个开花基因的调控网络,包含8个开花基因和23个LncRNA,涉及到温度和光信号调控通路。通过比较鉴定得到的LncRNA和972个甘蓝型油菜重要农艺性状QTL位...     

金花葵开花前后差异基因及代谢物分析

为深入研究金花葵花朵开花前后基因和代谢物的变化,利用转录组学和代谢组学技术相结合的方法对金花葵花蕾和花朵进行检测。结果显示,通过转录组分析鉴定了206 636个Unigenes,筛选出42 618个差异表达Unigenes,其中包括63个差异表达转录因子家族,24个转录调节因子家族。GO分析结果显示,差异表达基因主要富集于囊泡介导的逆行运输等生物学过程。KOG功能注释显示,差异表达基因功能以通用功能预测居多,其次是信号转导机制、翻译后修饰蛋白周转以及碳水化合物的转运和代谢。差异表达基因KEGG富集分析表明,差异基因主要富集于代谢、植物激素和信号转导、淀粉和蔗糖代谢等代谢通路。通过代谢组学检测,筛选到差异显著代谢物135个,主要包括脂类、氨基酸及衍生物和黄酮类等,KEGG富集分析表明,差异代谢物主要富集于甘油磷脂代谢、糖基磷脂酰肌醇(GPI)-锚定生物合成和次生代谢产物的生物合成-未分类等过程,其中,缬氨酸,亮氨酸和异亮氨酸的降解代谢通路同时出现在了基因和代谢物KEGG代谢通路富集Top 20中。     

基于RNA-seq的番茄响应根结线虫MiPDCD6蛋白的SA信号通路基因鉴定与表达分析

为了揭示南方根结线虫MiPDCD6蛋白抑制番茄PTI免疫的分子机理，以番茄品种新金丰1号MiPDCD6超表达苗为试验材料，以番茄品种新金丰1号组培苗为对照，通过转录组测序技术分别对番茄MiPDCD6超表达苗和对照苗进行转录组测序。以番茄栽培品种Heinz 1706基因组作为参考基因组进行对比，利用FPKM法计算基因表达量，设定参数(|log₂ FC|>1且P<0.05)筛选差异表达基因。利用Gene Ontology(GO)数据库进行差异表达基因GO功能富集分析，统计每个GO term中的差异表达基因数目，计算出基因富集的显著性，找出富集显著的功能条目。利用KEGG数据库对差异表达基因进行Pathway富集分析，超几何分布检验方法计算每个Pathway中差异表达基因富集的显著性。以FDR和基因数量来衡量KEGG的富集程度。结合基因差异表达分析和功能富集分析，研究MiPDCD6蛋白对番茄PTI免疫相关通路基因的影响。结果发现：在超表达MiPDCD6番茄植株中，与野生型番茄相比，MiPDCD6过表达番茄植株中有2 366个差异表达基因(DEGs),其中1 ...     

芦笋两性花与雄花的碳氮比及氮化合物相关基因可变剪接的差异分析

为探明芦笋两性花、雄花发育与碳氮比及氮化合物的关系,测定了芦笋两性花与雄花性别分化的前期、中期、后期3个时期花蕾的碳、氮化合物含量,并采用RNA-seq技术分析了6个样品中相关氮化合物基因的可变剪接及表达量的差异。结果发现,在性别分化的前期、中期,两性花蕾中含碳量低于雄花蕾,但含氮量却显著高于雄花蕾;在性别分化后期,两性花蕾中含碳量与雄花蕾相当,但含氮量仍然显著高于雄花蕾。可见,较低的碳氮比有利于两性花蕾中雌蕊的发育。RNA-seq测序及GO富集分析表明,有86个氮化合物相关基因只在两性花蕾中发生可变剪接,这些基因主要编码核糖体蛋白、MADS-box转录因子和一些酶类物质等,并且在性别分化的前期、中期,分别有52,47个基因表达量高于雄花蕾。这与两性花蕾中氮化合物含量高于雄花蕾的差异一致。综上所述,氮化合物代谢相关基因的可变剪接及较高的表达量和较低的碳氮比有利于两性花雌蕊的发育。     

DNA甲基化参与调控马铃薯响应干旱胁迫的关键基因挖掘

植物受到干旱胁迫时,会通过DNA甲基化做出快速反应以帮助其应对胁迫。为探究在干旱胁迫下,DNA甲基化是如何影响基因转录表达,本研究对甘露醇模拟干旱和5-azadC(去甲基化)处理下,抗旱性不同的2个马铃薯品种(抗旱型,青薯9号;干旱敏感型,大西洋)进行转录组学分析,以Fold-change>2和校正后P<0.01进行差异表达基因(DEG)的筛选。GO富集分析发现,2种处理都共同显著富集到氧化应激和碳水化合物代谢过程相关的GO term。说明不同耐旱性马铃薯在响应干旱胁迫时,与这些GO term相关的基因也受DNA去甲基化调控。对既响应干旱又响应DNA去甲基化的1345个DEG进行KEGG功能富集发现,与植物抗旱相关的通路有植物MAPK信号途径、植物激素信号转导途径、植物谷胱甘肽代谢通路、糖酵解与糖异生和磷酸肌醇代谢通路。说明这些通路相关基因在大西洋和青薯9号2个抗旱性不同的马铃薯品种中,响应干旱的敏感性受DNA甲基化调控。接着对DEG上游1500 bp启动子区域进行顺式作用原件和甲基化CpG岛分析发现,干旱胁迫下参与植物谷胱甘肽代谢的GST基因通过DNA去甲基化来降低启动子区ABRE和CAAT-box作用元件的甲基化水平,进而激活该基因的表达以应对干旱胁迫。因此,利用比较转录组学分析干旱和DNA去甲基化处理下的差异基因,可挖掘到DNA甲基化参与调控马铃薯响应干旱胁迫的相关基因,为研究马铃薯干旱胁迫响应的表观遗传学机理提供新的研究思路。     

甜玉米乳熟期果皮厚度相关基因的加权基因共表达网络分析

为深入探究甜玉米果皮厚度发育的分子机制,利用华南农业大学甜玉米课题组选育的果皮厚度差异较大的甜玉米自交系M03和M08为试验材料,对其授粉后15,19,23 d的果皮进行转录组测序。联合差异表达基因分析与加权基因共表达网络(WGCNA)分析,鉴定甜玉米乳熟期果皮厚度相关的共表达基因模块,根据模块内基因的连接度鉴定核心基因,利用qRT-PCR验证核心基因表达量的可靠性。比较M03和M08不同时期果皮转录组数据共筛选出4 748个差异表达基因(DEGs),DEGs主要富集到碳水化合物代谢过程、植物-病原体相互作用、植物MAPK信号通路。WGCNA分析共构建了18个共表达模块,通过筛选得到4个具有生物学意义且与果皮厚度显著相关(相关系数的绝对值>0.60)的特异性共表达模块,分别是Turquoise、Yellow、Magenta和Pink模块。GO和KEGG功能富集分析结果均表明,特异性模块富集到的天冬氨酸、丙氨酸和谷氨酸代谢、半胱氨酸和蛋氨酸代谢以及植物MAPK信号通路在甜玉米果皮厚度发育中具有重要的作用。选取连接度最高的前20个基因进行核心基因筛选,通过基因功能注释最终筛选到包括M...     

太平洋牡蛎对海洋酸化响应的分子机制——基于转录组的分析

旨在通过对太平洋牡蛎海洋酸化条件下转录组的分析,探讨其响应海洋酸化的分子机制。利用NCBI中太平洋牡蛎海洋酸化胁迫下的RNA-seq数据（包括4个组：pH 7.8、7.4、7.0、6.6）,通过Fastp、Hisat2、Samtools以及R语言分析,对其进行了差异基因、GO和KEGG等分析。结果显示：pH 7.4 vs pH 7.8组、pH 7.0 vs pH 7.8组及pH 6.6 vs p H 7.8组分别有61个、93个和943个显著差异表达基因。通过对pH 6.6 vs p H 7.8组显著差异表达基因进行GO分析发现这些基因显著富集于代谢、遗传信息处理和人类疾病等分子功能上,KEGG富集分析发现这些基因显著富集在糖代谢、脂类代谢、免疫机能上。表达分析表明,随着pH的下降,与能量代谢和免疫相关的基因表达量显著下降,如糖代谢中的糖原磷酸化酶的表达量从41.5几乎下降到0;氨基酸代谢中的精氨酸激酶1的表达量从1438几乎下降到0,精氨酸激酶2的表达量从27下降到3;脂肪代谢中的脂肪酸合成酶的表达量从28下降到2.5。这些基因的表达显著抑制可能使牡蛎机体能量代谢发生紊乱;此外,与免...     

转录组鉴定‘三白’石榴品种籽粒发育相关差异表达基因

为了研究石榴籽粒发育遗传机制及其在育种中的应用,以‘三白’石榴品种开花后60 d和120 d籽粒为试材,进行转录组测序分析。对表达差异基因进行GO注释及KEGG富集分析。结果分析表明,与‘三白’石榴开花后60 d籽粒相比,在开花后120 d籽粒中总共检测到8 409个显著差异表达基因;KEGG富集分析表明差异表达基因显著富集在激素信号转导及其它代谢产物的合成和代谢等18条代谢通路上;5类可变剪接类型中鉴定出了76个靶基因,其中17个为差异表达基因。定位到的差异表达基因可以为解析石榴籽粒发育以及基因克隆提供理论指导。     

氧化石墨烯促进玉米根系生长的转录组分析

为了研究氧化石墨烯(graphene oxide, GO)的生物效应,从转录水平解析低浓度GO促进植物生长的分子机制,本研究以玉米为材料,分析了土壤中施加不同浓度(0, 25, 50, 100, 200 mg/L) GO溶液14 d后的生长表型,筛选到GO对玉米生长的最佳的浓度(50 mg/L),进而利用最佳浓度GO处理48 h后的根系mRNA进行转录组测序分析。转录组测序共获得548个差异表达基因,其中194个上调表达,354个下调表达。Gene Ontology富集分析表明GO处理导致大量胁迫刺激响应基因的表达变化。在预测与根系生长相关的转录因子中,BRN2 (BEARSKIN2)、NAC2 (NAC domain-containing protein 2)、MYB93 (MYB domain protein 93)和PIF3 (Phytochrome-interacting factor 3)可能是GO的响应基因。GO处理使大量乙烯相关转录因子表达下调,说明GO处理可能抑制了体内的乙烯信号通路。在代谢相关基因中,一组木葡聚糖糖基转移酶/水解酶基因(XTHs)可能与根系生长相关。本研究为理解GO调控植物生长的分子机制提供线索,并为GO在农林业领域的应用提供参考。     

基于转录组测序枣胚发育过程差异表达基因分析

为了探索枣胚发育的分子机理,本研究使用Illumina高通量测序平台对中度可育枣品种‘冷白玉’四个发育时期(原胚,小球形胚,心形胚,鱼雷形胚)的幼胚进行转录组测序,共获得135 743 080条高质量Reads,Q3092.88%。样品间共筛选出2 333个差异表达基因,其中有19个差异表达基因是所有样品组合共有的。原胚与小球形胚的幼胚(L2 vs L3)间的差异表达基因(differentially expressed gene, DEG)基因数(867)明显多于其他组合,小球形胚与鱼雷形胚(L3 vs L5)间的DEG基因数(503)明显多于原胚与心形胚(L2 vs L4)间的DEG基因数(498)。GO富集分析显示,四个发育时期的DEG主要共同富集到分子功能、细胞组分、生物过程等GO条目。KEGG富集分析发现,四个时期的DEGs主要共同富集到次生代谢合成、代谢过程、植物信号转导通路,主要参与了苯丙氨酸代谢、α-亚麻酸代谢、氰氨基酸代谢、光合作用-天线蛋白质等代谢途径。通过转录组测序分析了不同发育阶段的‘冷白玉’枣幼胚差异基因表达状况,为寻找枣胚发育的相关基因提供理论参考。