基于盐胁迫转录组信息的蚕豆F-box基因家族分析

【目的】通过生物信息学方法分析蚕豆F-box基因家族成员的分布、结构及进化,研究家族成员在不同处理时间条件下的表达模式及对盐胁迫的响应,为该类基因生物学功能和盐胁迫机制的研究提供参考。【方法】基于蚕豆盐胁迫转录组测序(RNA-seq)数据,利用NR、Swiss-prot和PFAM 3个数据库和NCBI网站,对蚕豆F-box基因进行筛选注释;利用Web Logo 3、Prot Comp 9.0、MEGA-X和MEME等软件进行保守结构域、亚细胞定位、系统进化树和Motif等生物信息学分析。基于盐胁迫转录组数据分析蚕豆(yz17134耐盐和yz17078不耐盐)F-box基因家族在盐胁迫下的差异表达模式,并采用实时荧光定量PCR技术(qRT-PCR)检测部分家族成员在16和24 h的表达情况。【结果】基于盐胁迫转录组测序数据,注释得到161个蚕豆F-box基因,均含有F-box保守结构域。根据C端结构域的不同,将其分成11个亚族:FBX、FBXFBA、FBXLRR、FBXPP2、FBXKelch、FBXTUB、FBXFBD、FBXDUF、FBXACTIN、FBXWD40和FBO。保守结构域分析表明,F-box保守基序中包含1个极度保守的色氨酸残基。比较分析蚕豆F-box家族和拟南芥F-box家族共同构建的进化树,发现同一C端结构域的基因大多聚集在一起。亚细胞定位预测结果显示,124个F-box基因定位于细胞外,37个定位于细胞核中。基因结构分析表明,蚕豆F-box家族基因的DNA序列中均无内含子,且均由UTR区和CDS区组成。基于盐胁迫转录组数据的F-box差异表达模式分析表明,蚕豆F-box基因在2个不同处理时间点上的表达各不相同,在盐处理16 h的表达较为明显。q RT-PCR分析结果表明,在F-box家族成员中,共存在5个差异基因。其中Vf056266.1、Vf062764.1和Vf024236.1在盐处理16 h的表达量均上调,Vf060904.1和Vf045761.1在盐处理16 h的表达量均下调。【结论】蚕豆F-box基因家族注释得到161个蚕豆F-box基因,分为11个亚族。其中5个重要的F-box基因在不同盐处理时间的表达量存在差异。     

大白菜ACA基因家族的全基因组鉴定与表达分析

【目的】通过对大白菜ACA（Ca²⁺-ATPase）基因家族鉴定与表达分析,研究其家族基因间的共性与特性,为进一步揭示ACA家族进化关系提供数据支撑,为深入解析BraACAs在低温胁迫、盐胁迫以及自交不亲和方面的功能研究奠定基础。【方法】根据已报道的拟南芥ACA基因家族,同源比对出大白菜ACA基因家族,利用在线软件Expasy预测其分子量、理论等电点等理化性质;采用MEGA 5.0软件构建系统进化树;运用在线软件GSDS 2.0绘制基因结构图谱;TBtools对其染色体定位;McscanX软件进行拟南芥与大白菜ACA家族基因共线性分析;利用在线软件PlantCARE预测大白菜ACA基因家族启动子元件;通过在线工具Pfam和MEME进行蛋白保守结构域分析;利用qRT-PCR技术检测BraACAs在不同组织、非生物胁迫和自交异交授粉后的表达量。【结果】大白菜ACA基因家族有18个基因成员,分布在大白菜10条染色体上;根据进化树关系分成4组,分别包含3、4、4和7个成员;蛋白结构域分析显示,有13个成员包含N端自抑结构域。qRT-PCR结果表明,BraACAs主要在花与果荚中高表达;低温胁迫下,Bra002762与Bra035649表达量总体上调;盐胁迫下,Bra031701表达量显著上调;自交和异交授粉中,Bra003276和Bra024117差异性表达。对Bra002762、Bra035649、Bra031701、Bra003276和Bra024117亚细胞定位分析,发现这些基因均定位在细胞质膜上。【结论】大白菜ACA家族基因蛋白结构均含有4个ACA基因特有的高度保守结构域。该家族在大白菜不同组织中表达模式不同,5个ACA家族基因成员编码蛋白定位于细胞膜上,其中Bra002762、Bra035649、Bra031701与低温和盐胁迫响应有关;Bra003276和Bra024117与自交不亲和性相关。     

蚕豆自噬基因鉴定及响应干旱胁迫分析

【目的】根据转录组测序结果鉴定蚕豆自噬(Augophagy, ATG)基因,分析其对干旱胁迫的响应,为深入研究该类基因在干旱胁迫下的作用机制提供理论参考。【方法】利用蚕豆干旱胁迫转录组测序(RNA-seq)数据,根据swiss-port、NCBI、Ensemble等数据库的注释,筛选和鉴定蚕豆ATG基因;通过Smart、MEGA-X、MEME等软件进行保守结构域、亚细胞定位预测、系统发育树和motif等进行生物信息学分析。基于转录组数据,分析ATG基因在蚕豆CDAS105(耐旱)和鄂蚕1号(敏感)品种中干旱胁迫下的表达差异。【结果】筛选注释得到12个蚕豆ATG基因,均含有Autophagy结构域。根据C端结构域的不同将其分成5个亚族:ATG3、ATG6、ATG8、ATG18和ATG26。蚕豆ATG基因与拟南芥、玉米、水稻、酵母共同构建的进化树,发现同一C端结构域的基因聚集在一起。亚细胞定位预测结果显示:11个蚕豆ATG基因定位于细胞外,1个ATG蛋白定位于细胞质中。基因结构分析表明,蚕豆ATG家族基因的DNA序列中均无内含子,且均由UTR区和CDS区组成。保守基序分析鉴定出这些蛋白含有20个mortif。基于干旱胁迫转录组数据的ATG基因表达模式分析结果显示,ATG18b、ATG18c和ATG8a在干旱处理16或64 h后在抗旱品种中较敏感品种表达上调,且ATG8a在抗旱品种中随着干旱处理时间的增加表达进一步增强。【结论】从蚕豆中鉴定得到了12个ATG基因,分布于5个亚家族,ATG8a、ATG18b和ATG18c响应干旱胁迫,可以作为潜在的基因资源改良蚕豆种子萌发耐旱能力。     

玉米钙依赖蛋白激酶38(ZmCDPK38)的生物信息学及表达特性分析

【目的】 研究ZmCDPK38与其它物种中同源基因之间的遗传进化关系,通过与亲缘关系较近且功能已知的同源基因比较,对ZmCDPK38的功能进行理论预判。确定ZmCDPK38基因在不同组织和不同胁迫条件下的表达特性。【方法】 生物信息学分析使用在线分析工具及MEGAX软件;基因表达分析采用Real-time PCR。【结果】 ZmCDPK38与ZmCDPK28及高粱等6个其它植物物种的CDPK的氨基酸序列具有较高的一致性,遗传进化关系上较近。氨基酸序列包含4个结构域,具有CDPKs家族氨基酸的典型特征,Ser/Thr蛋白激酶结构域相对保守。ZmCDPK38基因在叶片中的表达量最高,其次是茎部、根、雄穗及雌穗;受盐胁迫及干旱胁迫诱导,ZmCDPK38基因的表达水平显著上升,在低温和高温胁迫条件下基因表达未见明显变化。【结论】 ZmCDPK38基因作为CDPK家族成员,在蛋白序列和结构上具备典型的CDPK家族成员特征;ZmCDPK38基因响应干旱胁迫与盐胁迫,可能在干旱应答与盐胁迫应答反应中发挥一定的功能。     

番茄DIR基因家族鉴定及其对非生物胁迫响应的分析

【目的】鉴定番茄DIR基因家族所有成员,并对其基因结构、编码蛋白的理化性质、系统进化、染色体定位、共线性、启动子元件、表达模式、互作转录因子及内源竞争RNA预测等进行了解析,为探究DIR在番茄生长发育和逆境胁迫中的作用提供参考。【方法】采用生物信息学方法,基于全基因组数据对番茄DIR基因家族成员进行鉴定,运用Phytozome、MEME、PlantCARE、opsRNATarget和plantcircnet等在线网站获取染色体位置、保守基序、顺式作用元件、存在互作的转录因子、miRNA和circRNA等信息。利用Mapteace、TBtools、Cytoscape、OmicShare等作图软件及工具绘制染色体定位图、进化树、DIR与转录因子关系图、ceRNA网络等。采用NCBI的基因数据库结合转录组测序及qRT-PCR试验研究DIR基因家族在逆境下的表达情况。【结果】共鉴定到番茄中27个DIR基因,将其命名为SlDIR1—SlDIR27,分别位于12条染色体上,且大部分基因位于染色体末端,其基因结构、基序和结构域相对保守,其中22个SlDIR具有一个外显子的经典结构。番茄DIR基因家族同拟南芥的共线性关系远高于水稻和大豆。基于系统进化关系将27个番茄DIR成员分为3个不同的亚家族。转录组数据表明大部分DIR基因在番茄根部具有较高的表达量。此外,DIR基因的启动子区域含有多个与干旱、低温等非生物胁迫,以及MeJA、ABA、SA等激素诱导相关的顺式作用元件,且预测到与激素、生长发育、非生物胁迫相关的ERF、E2F/DP、MYB等互作转录因子。结合转录组数据分析,分别有5、10、10和13个SlDIR在农药、干旱、盐和冷胁迫后显著上调表达,其中,SlDIR23受到以上4种胁迫的诱导表达,而SlDIR8、SlDIR13和SlDIR20特异性响应冷胁迫的诱导,SlDIR17特异性响应盐胁迫。番茄DIR的ceRNA调控表明,miR-156与靶基因SlDIR8可能共同作用调控番茄的逆境胁迫。【结论】共鉴定出番茄DIR家族基因27个,不均匀地分布在12条染色体上,在根部有较高的表达量。SlDIR1、SlDIR13和SlDIR14等具有MeJA、ABA、SA等激素响应元件,其中,SlDIR6仅具有MeJA元件,SlDIR27仅具有SA响应元件。另外,SlDIR2、SlDIR14、SlDIR23等参与干旱、盐、低温等多种逆境胁迫,其中SlDIR23在不同胁迫处理下均可被激活。此外,DIR基因和转录因子、非编码RNA相互作用,共同参与调控番茄植株的逆境胁迫。     

葡萄SAP家族的鉴定与表达分析

【目的】通过生物信息学分析明确SAP（Stress Associated Protein）家族在葡萄基因组中的数量、结构,利用qRT-PCR技术分析SAP家族的生物节律和在激素、非生物胁迫处理下的表达模式,为进一步探究SAP家族在葡萄中的作用奠定基础。【方法】在前期从山葡萄‘双优’（Vitis amurensis cv. Shuangyou）和欧洲葡萄‘红地球’（Vitis vinifera cv. Red Globe）0℃寒胁迫下转录组数据库中筛选出表达量明显上调的基因SAP15基础上,利用NCBI和葡萄基因组数据库中的BLAST功能,根据SAP家族的保守结构域,鉴定葡萄基因组中的SAP。利用生物信息学软件DNAMAN5.0、MEME、GSDS2.0、ExPASy和MEGA6等对VvSAPs序列、基因结构、蛋白质结构、理化性质、染色体定位和进化树等进行分析。采用qRT-PCR方法检测VvSAP家族的生物节律和在激素、逆境处理下的表达特征。【结果】从葡萄（Vitis vinifera）全基因组中鉴定得到15个SAP家族成员,均含有AN1保守结构域,大多含有A20保守结构域。按照其保守结构域和在染色体上的位置,依次命名为VvSAP1—VvSAP15,可分为Class Ⅰ、Class Ⅱ和Class Ⅲ三类。理化性质分析表明,VvSAP家族编码氨基酸数目介于109—293,理论等电点分布在7.99—9.68。染色体定位分析发现,该基因家族的15个基因分布于葡萄的9条染色体上,其中第8条染色体上分布最多,有3个VvSAP。亚细胞定位分析表明,VvSAP家族主要在细胞核、叶绿体和细胞骨架中进行表达。VvSAP家族的二级结构主要以无规则卷曲和α-螺旋为主。基因结构分析表明,VvSAP1—VvSAP12的DNA序列中均无内含子,VvSAP13—VvSAP15的DNA序列中有1个内含子,基因结构高度保守。qRT-PCR分析结果表明,在VvSAP家族成员中,VvSAP1和VvSAP9在无处理和各处理下（激素和非生物胁迫）均表达极低或不表达,初步鉴定为假基因。除VvSAP10在400 mmol·L ^-1 NaCl处理下呈下调表达,其余基因均在50 μmol·L ^-1 ABA、100 μmol·L ^-1SA和400 mmol·L ^-1NaCl处理下呈上调表达,其中,VvSAP10—VvSAP14显著响应50 μmol·L ^-1 ABA胁迫,处理24 h后相对表达量分别上调为0 h的37.19、36.63、21.69、58.34和267.35倍;VvSAP8和VvSAP11在NaCl处理4 h后相对表达量最高,分别为0 h的13.16和12.42倍;在4℃低温胁迫下,相对表达量上调最高的是VvSAP15,处理后8 h相对表达量为0 h的35.90倍。 【结论】从葡萄基因组中共鉴定出15个SAP家族成员,初步鉴定得出VvSAP1和VvSAP9为假基因,15个基因分别分布于9条染色体上,并且结构进化高度保守。所有成员均相响应逆境,且均有昼夜节律变化,但该基因家族成员在不同逆境胁迫下的响应程度和在不同时间点的表达模式存在一定的差异。     

芥菜全基因组中Hsf基因家族鉴定与分析

【目的】通过鉴定与分析芥菜基因组中Hsf转录因子家族成员,为芥菜Hsf基因功能研究与遗传改良提高抗逆性提供理论依据。【方法】利用生物信息学方法分析芥菜Hsf热激转录因子家族成员的功能结构、保守基序、启动子顺式作用元件、系统进化、共线性及采用RNA-seq验证Hsf低温胁迫的基因表达。【结果】芥菜基因组中鉴定出71个Hsf基因家族成员,分布于18条染色体上,归类为3个亚家族,蛋白均含有DBD和HR-A/B结构域。BjuHsf启动子区域包含与逆境胁迫、激素、生长发育等相关顺式作用元件。系统进化与共线性分析表明,芥菜Hsf家族成员与大白菜Hsf家族成员具有更近的亲缘关系。在低温胁迫下,BjuHsf基因表达分析表明8个BjuHsf基因显著上调表达。【结论】这些显著差异表达BjuHsf基因与芥菜抗寒性极其相关,可作为芥菜耐寒性遗传改良的候选基因。     

黄瓜ERF基因家族鉴定及其在雌花芽分化中的表达分析

【目的】通过以黄瓜9930_V2版本基因组为参照进行生物信息学分析,对ERF基因家族在基因组中的数量、结构以及表达特征进行分析,为研究ERF转录因子在黄瓜雌花分化与发育中的作用提供数据支持。【方法】根据已报道的拟南芥ERF,利用黄瓜基因组数据库中9930_V2版本基因组进行BLAST比对,通过MEGA、MEME、TBtools、ExPASy等工具进行生物信息学分析。采用qRT-PCR方法检测不同性型黄瓜材料、雌花发育初期不同阶段中ERF基因家族成员的表达水平。采用酵母单杂交方法验证家族成员与乙烯响应元件GCC-box的互作。【结果】从黄瓜材料9930基因组中鉴定得到138个ERF基因家族成员,共分为10个亚族,编码氨基酸长度介于126—745。按照基因家族成员在染色体上的位置分布,将其命名为CsERF1-CsERF138。多序列比对和motif分析结果表明,黄瓜ERF基因家族均具有AP2/ERF结构域,其中4个成员具有B3结构域。表达分析结果显示,在不同性型材料中共有19个ERF家族成员差异表达,其中9个在FFMMAA基因型中高表达,10个在ffMMAA基因型中高表达。通过雌花芽发育初期ERF家族成员的表达趋势分析,发现31个ERF随子房发育表达上调,30个表达下调。初步证明CsERF9和CsERF31具有结合GCC-box元件的功能。【结论】从黄瓜基因组中鉴定出138个ERF基因家族成员,均拥有1个或多个AP2/ERF结构域;其中部分成员在不同性型材料中差异表达,并可能参与雌花分化初期的基因表达调控;部分成员具有结合保守元件GCC-box调控下游基因表达的功能。     

盐胁迫下盐穗木DNA聚合酶λ基因的克隆和表达分析

【目的】开展盐胁迫下DNA损伤修复基因的表达研究,有助于揭示DNA损伤修复与盐生植物耐盐的相关性。【方法】根据盐胁迫下盐穗木转录组测序结果,利用RACE技术克隆获得了盐穗木DNA损伤修复基因HcDNA聚合酶λ(HcDNA polλ)基因,开放阅读框1 335 bp,编码444个氨基酸。【结果】保守结构域分析显示,HcDNA polλ基因属于DNA聚合酶X家族成员,系统进化分析显示HcDNA polλ为独立的分支,亚细胞定位于细胞核,无信号肽且不含跨膜区域的亲水蛋白。实时荧光定量PCR分析表明,随着NaCl胁迫浓度的增加,同化枝的HcDNA polλ基因的表达逐渐上调,在300 mmol/L NaCl处理7和14 d时,HcDNA polλ的表达分别增加3和5倍。最后在700 mmol/L NaCl处理14 d时,HcDNA polλ的表达增加20倍,达到峰值。【结论】HcDNA polλ基因的表达受盐胁迫的诱导。     

西伯利亚杏全基因组CBF基因家族鉴定及低温胁迫下的表达模式

【目的】利用西伯利亚杏全基因组数据库鉴定CBF家族成员，并分析CBF转录因子保守域特点与功能及其在低温胁迫下的表达情况。【方法】根据西伯利亚杏全基因组数据库，利用BLASTP和HMM Search搜索西伯利亚杏CBF基因，通过CDD和MEME验证其保守结构域。使用Prot-Param、WOLF PSORT、MAGA X、Phyre2、PlantCARE等生物信息学工具分别分析蛋白理化性质、亚细胞定位、系统发育、蛋白三级结构、染色体定位及顺式调控元件。利用实时荧光定量PCR分析西伯利亚杏CBF基因在低温胁迫后的表达模式。【结果】本研究共鉴定出4个西伯利亚杏CBF基因，它们均具有AP2基因家族保守结构域PKKPAGR和DSAWR序列。西伯利亚杏CBF蛋白包含206～246个氨基酸，蛋白理论分子量平均值为25.80 KDa,平均等电点为5.92。PsCBF3和PsCBF4在基因结构上较为相似，亲缘关系更近。4个PsCBFs亚细胞定位均位于细胞核内。启动子顺式作用元件分析结果显示CBF家族成员可参与植物激素和非生物胁迫响应。西伯利亚杏CBF成员均属于酸性和亲水性蛋白。低温胁迫后，CBF基因主要...     

盐穗木miR167d和预测靶基因ARF8在盐胁迫下的表达与相关性分析

【目的】盐穗木(Halostachys caspica)是一种藜科多年生盐生灌木,对盐分适应性极强。miR167是作用于生长素信号通路上的内源非编码小RNA分子,它通过靶向生长素响应因子ARFs(auxin response factors)来调控生长素响应基因的表达。研究从高盐(600 mmol/L NaCl)处理48 h盐穗木根的小RNA文库中筛选到差异表达的miR167d,并从该物种的转录组数据中预测到其靶基因ARF8。该文开展了高盐胁迫下盐穗木miR167d和预测靶基因ARF8的表达模式和相关性分析。【方法】通过荧光定量PCR试验检测和分析高盐胁迫下盐穗木miR167d和预测靶基因的表达模式和相关性;利用烟草瞬时表达试验间接地鉴定盐穗木miR167d对预测靶基因ARF8的靶向关系;采用同源克隆方法结合RACE技术获得盐穗木miR167d的靶基因ARF8的全长序列,并进行生物信息学分析。【结果】(1) 600 mmol/L NaCl处理盐穗木,其同化枝中miR167d和ARF8的表达均受到诱导,胁迫48 h时,miR167d的表达最高并与靶基因的表达呈现显著的负相关性,HcARF8基因的表达随胁迫处理时间的延长呈现先增加后下降的趋势,推测盐穗木miR167d和ARF8两基因可能在响应盐胁迫过程中发挥作用。(2) 构建融合GFP的含预测靶基因HcARF8的植物表达载体并转化农杆菌,烟草瞬时表达试验的结果表明,拟南芥miR167d对盐穗木ARF8基因具有剪切作用,间接证明盐穗木miR167d对预测的靶基因HcARF8具有靶向切割作用。(3) 克隆获得的盐穗木ARF8基因全长2 861 bp,ORF 2 442 bp,编码813个氨基酸,在编码区与盐穗木miR167d高度互补匹配。生物信息学分析该基因编码的蛋白高度保守,与同科植物甜菜ARF8同源性达到87%,都具有能与生长素相关元件(B3 DNA绑定元件、生长素响应元件,生长素诱导转录IAA超家族的元件)结合的功能域。【结论】盐穗木miR167d和HcARF8基因具有靶向关系,高盐胁迫处理盐穗木,其同化枝中两基因的表达都受到诱导并呈现显著的负相关性。这些结果为后续阐明盐穗木miR167d与其作用的靶基因ARF8的生物学功能奠定基础。     

柑橘Rboh家族鉴定及其对激素和柑橘溃疡病的响应

【目的】活性氧（reactive oxygen species,ROS）是植物抗病反应中重要的信号分子,而Rboh是参与活性氧产生的关键酶之一。本研究通过鉴定和分析柑橘全基因组中Rboh基因家族生物信息学特性、生物胁迫相关植物激素和溃疡病菌（Xanthomonas citri subsp. citri,Xcc）侵染诱导下表达模式,为研究Rboh基因家族与柑橘抗溃疡病过程的相关性,进一步研究和利用柑橘Rboh基因打下基础。【方法】利用过氧化物酶数据库RedOxiBase和柑橘（甜橙）CAP数据库获得柑橘Rboh家族序列信息;利用生物信息学软件对CsRboh进行理化性质、系统进化关系、染色体定位、基因结构、蛋白功能结构域、保守基序和启动子元件系统分析。通过实时荧光定量PCR（qRT-PCR）分析植物激素水杨酸（salicylic acid,SA）、茉莉酸（jasmonic acid,JA）、脱落酸（abscisic acid,ABA）和溃疡病菌处理下CsRboh表达模式。【结果】鉴定得到7个柑橘Rboh家族成员（CsRboh01—CsRboh07）,其编码蛋白包含784—946个不等的氨基酸残基,等电点（PI）分布在8.67—9.40,主要定位于细胞膜结构和细胞质。根据系统进化树可将柑橘Rboh家族划分为I—IV 4个亚家族。CsRboh家族基因不均匀分布在甜橙的5条染色体,均含有典型的呼吸爆发NADPH氧化酶结构域（NADPH_Ox）、铁还原酶跨膜组件（Ferric_reduct）、FAD结合结构域（FAD_binding）和铁还原酶NAD结合结构域（NAD_binding）4个功能结构域。CsRboh家族各基因的启动子区域含不同激素响应元件,数目和类型各有差异。qRT-PCR结果显示CsRboh基因家族各成员可响应激素和溃疡病菌诱导,在不同激素和溃疡病菌处理下感病品种晚锦橙和抗病品种四季橘中表达模式具有差异。结合系统进化树聚类分析、同源基因功能研究及其启动子区域顺式作用元件分析发现,CsRboh02、CsRboh04和CsRboh06在抗、感两个品种中受柑橘溃疡病菌诱导表现出明显不同的表达趋势和表达量。【结论】CsRboh02、CsRboh04和CsRboh06可能与柑橘品种抗、感性密切相关,是3个有潜力的抗溃疡病分子育种候选基因,初步确定CsRboh家族基因在柑橘响应溃疡病过程中发挥关键作用。     

苹果果实大小相关的ARF-Aux/IAA互作组合筛选

【目的】通过转录组学和生物信息学,在苹果全基因组中对可能互作的MdARF和MdIAA进行鉴定,为明确相关基因功能和解析生长素调控苹果果实大小的分子机理奠定基础。【方法】对野生大果型‘皇家嘎啦’和miRNA172p过表达的转基因小果型‘皇家嘎啦’进行不同发育时期和不同组织材料的转录组测序,对测序结果进行基因的功能注释及差异表达分析。利用转基因小果和野生型大果的转录组数据,筛选在果实发育中表达的苹果MdARFs和MdIAAs基因家族成员,通过逐一在两个家族间计算基因时空表达的相关系数,筛选可能互作的MdARF-MdIAA组合,将从拟南芥基因组中下载的23个ARFs和34个Aux/IAAs、番茄基因组中下载的21个ARFs和25个Aux/IAAs,分别与互作候选MdARFs和MdIAAs进行比对,并进一步构建系统发育树。使用MEME和TBtools对苹果候选互作对中的MdARFs和MdIAAs蛋白进行Motif分析。利用STRING蛋白互作预测数据库进行同源映射,构建苹果中的蛋白-蛋白互作网络,进一步的确认候选互作对,最终得到苹果中通过互作参与果实发育可能性最高的MdARF-MdIAA组合。【结果】分别对野生型‘皇家嘎啦’和miR172OX转基因‘皇家嘎啦’盛花期后两周的全果和盛花期后4周的果皮、果肉和果核进行转录组测序,共生成178.19 Gb的数据量,各项指标均表明,3个生物学重复在所有组织类型上均具有高度一致性。在转录组数据中,共鉴定到38个MdARFs和27个MdIAAs在至少一个文库中的FPKM值大于2,在苹果果实发育时期表达。通过计算Pearson相关系数对表达的MdARFs和MdIAAs两两进行相关性分析,其中8对MdARF-MdIAA的相关系数大于0.9或小于-0.9,作为初步筛选的候选互作组合。将8对组合中的MdARFs和MdIAAs分别与拟南芥和番茄中的ARFs和IAAs进行序列比对并构建系统进化树后发现,MdARF6和MdARF19与起转录激活作用的AtARFs同属一个分支。而MdARF2、MdARF4和MdARF9则与起转录抑制作用的AtARFs具有较近的亲缘关系。Motif分析结果显示,候选MdARF、MdIAA蛋白中均包含Motif 2和Motif 5。Motif 2和Motif 5分别对应IAA蛋白中的保守结构域Motif IV和Motif III。互作蛋白在拟南芥中进行同源映射校验后,最终得到两对MdARF-MdIAA组合可用于进一步的功能验证。【结论】苹果MdARF和MdIAA家族成员,在果实发育时期,有8对组合在表达量上存在显著的相关性,进一步同源映射确认互作后,最终确定MdARF4-MdIAA17和MdARF4-MdIAA19两对互作组合,极有可能通过互作传递生长素信号参与调控苹果果实发育。     

盐胁迫下盐穗木miR393b与预测靶基因HcTIR1的表达及相关性

【目的】研究盐胁迫下盐穗木生长素通路中miR393b和预测的靶基因TIR1的表达模式和相关性。【方法】采用生物信息学方法预测盐穗木miR393b的靶基因;通过qRT-PCR方法检测盐胁迫下盐穗木miR393b及预测靶基因HcTIR1(transport inhibitor response1)的相对表达模式,并分析其相关性。【结果】miR393b成熟体在不同植物种中高度保守;盐穗木miR393b预测的靶基因为TIR1;盐胁迫处理,盐穗木同化枝中两基因响应高盐胁迫,并具有一定的负相关性。【结论】盐穗木miR393b和预测的靶基因响应高盐胁迫,二者具有一定的负相关。该结果为进一步探讨盐穗木miR393b与预测靶基因TIR1的生物学功能奠定基础。