棉花β-tubulin基因家族的鉴定及其在纤维发育中的表达

【目的】β-微管蛋白是棉纤维细胞形态建成的基本结构单位，在纤维发育过程中具有重要作用。通过鉴定棉花β-tubulin基因家族成员，并进行生物信息学和表达模式分析，为深入探究β-tubulin基因在棉花纤维发育中的作用提供理论依据。【方法】利用BLAST方法，在4个棉种基因组中鉴定β-tubulin基因家族成员，并结合ProtParam tool分析理化性质、MEGA7.0构建系统进化树、Mapchart2.2绘制染色体定位图、MEME分析保守基序、PlantCARE分析启动子顺式作用元件；根据39个材料发育纤维转录组数据分析陆地棉β-tubulin基因家族的表达水平，并利用Spearman相关性分析鉴定影响纤维品质性状形成的β-tubulin基因。【结果】在陆地棉（Gossypium hirsutum,AD1）、海岛棉（Gossypium barbadense,AD2）、亚洲棉（Gossypium arboretum,A2）和雷蒙德氏棉（Gossypium raimondii,D5）基因组中分别鉴定到36、37、19和18个β-tubulin基因，且在四倍体棉种中的数目约是二倍体棉种数...     

亚洲棉Argonaute蛋白家族的生物信息学预测与分析

【目的】在全基因组水平上鉴定亚洲棉(Gossypium arboretum L.)AGO蛋白家族成员,分析其进化关系,为亚洲棉遗传改良提供基因资源。【方法】用AGO蛋白保守域在亚洲棉基因组数据库中进行同源比对,并对AGO蛋白家族的进化关系、基因结构、结构域、染色体定位、三级结构及表达模式等进行生物信息学分析。【结果】从亚洲棉蛋白数据库中共鉴定出11个AGO蛋白家族成员,其编码区长度为2 949~9 807bp,定位在6条染色体上,其中第9号染色体上分布有4个基因;11个AGO家族成员分为5个亚族,亚族内各成员间具有很高的保守性;AGO家族成员均具有典型的α-螺旋跨膜结构,属于碱性蛋白;从表达模式看,AGO家族成员在纤维发育不同时期和不同组织中的表达量存在较大差异,其中5个基因在棉纤维次生壁后期高表达,在其他时期表达水平较低,存在组织特异性。【结论】成功鉴定了亚洲棉AGO蛋白家族,明确了其基本结构的进化关系及基因表达特点。     

棉花蔗糖转化酶基因家族的生物信息学分析

为更好地了解蔗糖转化酶在棉花基因组中的数量和分布情况等,本研究利用拟南芥中17个蔗糖转化酶家族的基因为参考序列,在海岛棉、陆地棉、雷蒙德氏棉和亚洲棉基因组数据库中鉴定到的蔗糖转化酶基因家族成员分别为65、46、26和25个,结果表明,海岛棉所含有的蔗糖转化酶基因明显多于陆地棉,而雷蒙德氏棉比亚洲棉仅多了1个。通过基因结构分析发现,酸性蔗糖转化酶和中性/碱性蔗糖转化酶α亚族大都包含6或7个外显子,而β亚族大多含有4或5个外显子。此外,部分海岛棉与雷蒙德氏棉在基因结构上存在着较长的UTR区,而陆地棉与亚洲棉可能没有;酸性蔗糖转化酶家族基因所编码的氨基酸序列大都包含3个保守基序,分别为DNPNG、FRDP和WECPD,表明不同种属中同一基因家族的进化保守性;由聚类分析推测,VIN亚族基因功能可能与棉花纤维发育相关。     

芥菜TCP转录因子家族全基因组鉴定及表达分析

【目的】对芥菜基因组中TCP转录因子家族基因成员进行鉴定和序列分析,并分析TCP家族基因在茎发育时期的表达水平,为研究芥菜TCP基因功能和茎膨大机理提供基因源和理论依据。【方法】从芸薹属基因组数据库下载芥菜基因组数据,结合转录组数据,对芥菜TCP转录因子家族成员进行鉴定,并利用生物信息学方法对该家族成员进行理化性质、亚细胞定位、系统发育关系、基因结构、染色体定位及茎不同发育时期基因表达谱等方面的分析。【结果】芥菜TCP家族共有53个成员;系统进化树分析可分为两个大的亚组;TCP基因大都位于细胞核上,但理化性质差异很大;大部分基因没有内含子,基因结构保守性较高;在染色体上分布不均;表达模式分析显示TCP家族基因可能在芥菜茎膨大中起重要作用。【结论】鉴定了53个芥菜TCP家族基因成员,并发现PCF亚组成员在芥菜茎发育时期特异性表达,可能在芥菜茎的生长发育过程中起重要作用,该研究为芥菜及其他植物中TCP家族的鉴定和功能研究提供了理论基础。     

陆地棉GhRTE基因家族的鉴定与耐盐性分析

RTE(reversion-to-ethylene)是一种植物乙烯响应的负调节因子,参与植物生长发育并在非生物胁迫的逆境中起着重要作用。通过对陆地棉(Gossypium hirsutum)RTE家族成员的鉴定及其在盐胁迫下的表达模式进行解析,为进一步探索棉花RTE家族成员的功能提供一定的研究指导。通过棉花全基因组数据,鉴定棉花RTE家族成员,以陆地棉为主要研究对象,结合海岛棉(Gossypium barbadense)、亚洲棉(Gossypium arboreum)、雷蒙德氏棉(Gossypium raimondii)等植物的RTE基因家族成员进行物种内的进化分析。同时,利用转录组学方法,研究陆地棉RTE家族基因的表达量及其对盐胁迫的应答规律。结果表明,棉花RTE蛋白家族共鉴定出24个RTE基因,其中来自陆地棉和海岛棉的各有8个,有4个来自雷蒙德氏棉,4个来自亚洲棉。系统发育树分析将24个RTE基因分成3个亚族,GhRTE基因同一亚族成员间具有相似的基因结构和保守基序。GhRTE启动子区包含了许多与植物激素反应、生长发育和逆境应答密切相关的顺式调控因子。通过对转录组数据分析以及荧光定量...     

陆地棉TCP家族基因鉴定及组织表达分析

【目的】鉴定、分析陆地棉TCP家族基因,为进一步研究TCP基因在棉花植株内的生物学功能奠定基础。【方法】基于2019年最新组装的陆地棉TM-1参考基因组,通过Pfam网站获取TCP家族HMM模型文件,使用HMMER网站鉴定陆地棉TCP基因,CottonFGD网站获取陆地棉TCP基因组蛋白序列。利用MapInspect进行染色体定位、MEGA 7.0进行多序列比对聚类分析、MEME网站motif预测、TBtools软件鉴定基因结构以及陆地棉TCP基因组织特异性表达分析。【结果】共鉴定出63个陆地棉TCP基因,其中39个Class I亚族,24个ClassⅡ亚族;ClassⅡ亚族中包括17个CIN亚类,7个CYC/TB1亚类。染色体定位发现该63个TCP基因在陆地棉22条染色体上均有分布,其中A组染色体上分布有33个基因、D组上分布有30个基因。所有TCP蛋白均包含TCP结构域。TCP基因外显子、内含子结构及长度在同一亚家族内具有相似性。TCP家族基因中有12个基因在陆地棉纤维中优势表达,32个基因在花器官中优势表达,16个基因在营养器官：根、茎和叶中优势表达。【结论】获得了与陆地棉生长发...     

玉米ZmEXO70s基因家族鉴定及苗期耐热性关联分析

【目的】鉴定玉米EXO70s（ZmEXO70s）基因家族成员,分析其遗传变异与玉米耐热性的关联性,为揭示其在玉米耐热方面的功能和分子机制打下基础。【方法】以拟南芥EXO70s（AtEXO70s）基因家族成员序列为参考,利用MaizeGDB和NCBI数据库从玉米全基因组中鉴定出ZmEXO70s基因家族成员,对其进行基因结构及系统进化分析,根据其在系统发育进化树上的位置对其进行系统命名,并基于转录组测序（RNA-Seq）数据分析ZmEXO70s基因家族成员在不同组织及非生物胁迫下的表达模式。最后,鉴定玉米自交群体AM368中ZmEXO70s基因SNP位点,结合玉米AM368群体苗期耐热存活率进行基因家族关联分析。【结果】从玉米全基因组序列中共鉴定出36个ZmEXO70s基因,分布在10条染色体上,长度为228～3726 bp,编码75～1241个氨基酸残基,蛋白分子量为8.6～136.6 kD,理论等电点（pI）介于4.5～9.9,大部分蛋白pI小于7.0。拟南芥、水稻和玉米的EXO70s蛋白被分为9组（A组～I组）,其中ZmEXO70s蛋白在各组中均有分布。大多数ZmEXO70s基因不含内含子,只有少数基因含有内含子,且内含子数量不同。36个ZmEXO70s蛋白共含8种基序（Motif 1～Motif 8）,主要集中在C端,说明这些蛋白端序列较保守,且Motif 1～Motif 8组成Exo70保守结构域。ZmEXO70s基因在不同组织中的表达水平均存在明显差异,其中ZmEXO70B3基因在雄穗、花药、叶片和花丝中高表达,ZmEXO70C1a基因只在雄穗和花药中高表达,ZmEXO70D2b基因在胚乳和萌发种子中低表达,ZmEXO70G1c基因在花丝和萌发种子中低表达。有23个ZmEXO70s基因至少可响应一种非生物胁迫,说明这些ZmEXO70s基因参与非生物胁迫响应过程。ZmEXO70D2a和ZmEXO70E1基因的遗传变异与玉米苗期耐热性具有显著相关性（P≤0.01）。【结论】ZmEXO70s基因家族成员在系统发育进化上较保守,其基因组复制事件可能发生在禾本科植物分化后,且大多数ZmEXO70s基因被保留,仅部分基因被丢失。ZmEXO70s基因可能在非生物胁迫响应过程中发挥重要作用,ZmEXO70D2a和ZmEXO70E1基因可作为调控玉米苗期耐热性重要候选基因。     

黄瓜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调控下游基因表达的功能。     

海岛棉GH9基因家族成员鉴定及分析

植物GH9基因在细胞壁的生物合成和重塑中起着重要作用,对海岛棉GH9基因进行全基因组鉴定与分析,挖掘其在棉纤维发育中的作用。从海岛棉全基因组中鉴定到53个GH9基因,分析其理化性质、基因结构、染色体分布、进化历程、棉纤维发育过程中表达模式及转录调控。结果表明,53个GbGH9s分为A、B、C三组,定位在22条染色体上;多倍化事件是海岛棉GH9基因家族扩张的主要驱动力,基因复制后经历严格的选择约束;A组多个成员在棉纤维发育整个时期尤其是次生壁加厚期高表达,C组多个成员在棉纤维起始和伸长期优势表达,B组具有相对复杂的表达模式,生长素和乙烯可能在GbGH9s的转录调控中发挥重要作用。此外,通过回交将GbGH9B6导入陆地棉可提高棉纤维强度。海岛棉GH9基因家族成员的鉴定及分析有助于明确其进化和功能,为后续研究提供重要参考。     

海岛棉与陆地棉不同纤维发育时期苯丙烷代谢相关基因表达差异比较

[目的]比较研究在海岛棉与陆地棉纤维发育过程中苯丙烷代谢相关基因(PAL和4CL)表达模式的差异.[方法]利用荧光定量PCR,检测不同品种来源纤维在不同发育时期目的基因表达量.[结果]发现PAL1及4CL1基因的表达量随纤维发育天数的增加而增加,在陆地棉军棉1号中,基因表达量在15 DPA达到最高,而在海岛棉新海20号中PAL1及4CL1基因的表达高峰期则为20DPA.[结论]发现PAL1基因与4CL1基因的表达模式在海岛棉与陆地棉纤维发育过程中出现差异,海岛棉中两个基因的表达高峰期均晚于陆地棉,为逐步揭示海岛棉和陆地棉纤维品质差异的分子机理奠定了一定实验基础.     

棉花FLA基因家族的全基因组鉴定及Gh FLA05在棉纤维发育中的功能分析

【背景】伴随着棉纺织工艺水平的提升和人们对高品质纺织品的追求，提升棉花纤维品质日益重要。类成束阿拉伯半乳糖蛋白（fasciclin-like arabinogalactan proteins,FLAs）在棉纤维起始发育、次生壁合成等过程中可能具有重要作用。【目的】通过对棉花FLA基因家族进行全面鉴定与分析，研究该家族成员的共性特征及特异性表达模式，为FLA在棉纤维发育中的功能研究提供参考。【方法】根据棉花全基因组数据，使用HMMER 3.0对棉花FLA基因家族成员进行鉴定，并通过Pfam、Smart等软件进一步确认。使用ExPASy、TMHMM分析蛋白理化性质及跨膜结构域，应用MEGA、MCScanX、GSDS、MEME、TBtools、Jalview等工具进行进化树构建、染色体定位、共线性分析和蛋白保守结构域序列比对等。通过转录组数据分析陆地棉FLA基因在不同组织中的表达情况。利用实时荧光定量聚合酶链式反应（quantitative real-time polymerase chain reaction,qRT-PCR）检测GhFLAs在不同纤维品质材料的胚珠及不同发育时期纤维中的表...     

水稻SBP基因家族的全基因组研究

【目的】水资源紧缺严重威胁我国的粮食安全,优化水稻的节水抗旱性状是推广旱作农业节水的最为经济有效的方法。Squamosa promoter binding protein(SBP)基因家族是一类具有重要应用价值的植物特异性转录因子,它们可以影响植物的株型、产量和抗逆性,为节水抗旱育种提供了重要资源。然而,SBP家族成员在水稻发育和胁迫响应中的功能目前还未全部鉴定,调控抗逆性研究也未进行。【方法】利用SBP结构域一致性序列搜索水稻基因组,再进行基因结构和系统发育分析、共线性分析及基因表达模式分析。【结果】鉴定出19个SBP-like(OsSPL)基因,分布于10条染色体上,编码蛋白序列长度为217～1 141 aa。基因结构和系统发育分析表明水稻OsSPL基因可分成4组。共线性分析表明,水稻基因组内OsSPL-4和OsSPL-5形成1个串联重复基因簇,11个OsSPL,形成7个跨染色体分布的重复基因对。水稻和拟南芥SPL形成5个同源基因对,与玉米形成45个同源基因对。Evorepro数据库中OsSPL基因表达模式分析表明水稻叶片冷胁迫下OsSPL-1和OsSPL-3表达下调,OsSPL-...     

黄瓜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调控下游基因表达的功能。     

海岛棉类受体胞质激酶RLCK VI家族全基因组鉴定与胁迫表达分析

【目的】对海岛棉（Gossypium barbadense）类受体胞质激酶RLCK VI(GbRLCK VI)家族基因进行全基因组分析，为深入研究RLCK VI家族基因参与棉花生长发育和抗逆的调控机制提供参考。【方法】基于最新发布的海岛棉基因组数据，利用生物信息学手段对GbRLCK VI家族基因进行全基因组鉴定，并系统分析该家族基因成员的理化性质、序列特征、基因复制、系统进化和表达特征。【结果】海岛棉中共鉴定出39个GbRLCK VI家族基因，经聚类分析将其分为A、B两组，其中A组22个，B组17个，均含有1个激酶结构域，分布于16条染色体，多数位于质膜。基因复制分析表明，该家族在进化过程中发生染色体片段复制事件；Ka/Ks分析显示，所有基因对Ka/Ks均小于1，表明GbRLCK VI家族基因在进化过程中可能经历了严格的纯化选择作用。转录组分析表明，GbRLCK VI家族基因在10个不同组织中的表达模式不同，11个基因在花器官中优势表达，9个基因在根茎叶中优势表达；逆境胁迫下的表达分析显示，8个基因在干旱、盐、黄萎病胁迫下优势表达，4个基因只在黄萎病胁迫下优势表达，说明GbRLCKVI...     

玉米CesA家族的系统鉴定及功能研究

【目的】系统鉴定玉米纤维素合成酶(CesA)家族成员,明确其进化关系及基因功能。【方法】采用系统生物学方法,在全基因组水平上对玉米CesA家族成员进行鉴定,进一步结合基因结构变化对CesA家族成员的进化关系进行系统阐述;通过RNA测序和核糖体图谱,系统分析玉米CesA家族成员响应干旱胁迫的表达模式,并从基因复制角度分析该家族基因的扩增方式,同时基于构建的转录调控网络,对玉米CesA基因的功能进行初步探索。【结果】鉴定出16个玉米CesA家族成员,分为7个子家族,其中有3个是新发现的成员;CesA家族成员在玉米发育过程中对干旱胁迫的响应有明显的组织特异性和发育时期特异性,干旱胁迫对玉米种子发育过程的影响较对植株生长过程的影响大;CesA家族成员在合成纤维素过程中相互协同作用,成员间存在功能冗余。【结论】系统鉴定出了16个玉米CesA家族成员,进一步明确了其彼此之间的功能协作及其在种子发育和植株生长过程中对干旱胁迫的响应程度。     

玉米CesA家族的系统鉴定及功能研究

【目的】系统鉴定玉米纤维素合成酶（CesA）家族成员,明确其进化关系及基因功能。【方法】采用系统生物学方法,在全基因组水平上对玉米CesA家族成员进行鉴定,进一步结合基因结构变化对CesA家族成员的进化关系进行系统阐述;通过RNA测序和核糖体图谱,系统分析玉米CesA家族成员响应干旱胁迫的表达模式,并从基因复制角度分析该家族基因的扩增方式,同时基于构建的转录调控网络,对玉米CesA基因的功能进行初步探索。【结果】鉴定出16个玉米CesA家族成员,分为7个子家族,其中有3个是新发现的成员;CesA家族成员在玉米发育过程中对干旱胁迫的响应有明显的组织特异性和发育时期特异性,干旱胁迫对玉米种子发育过程的影响较对植株生长过程的影响大;CesA家族成员在合成纤维素过程中相互协同作用,成员间存在功能冗余。【结论】系统鉴定出了16个玉米CesA家族成员,进一步明确了其彼此之间的功能协作及其在种子发育和植株生长过程中对干旱胁迫的响应程度。     

海岛棉与陆地棉纤维C4H1基因克隆及表达分析

[目的]研究海岛棉与陆地棉纤维不同发育时期C4H1基因的表达和结构差异.[方法]分别从海岛棉和陆地棉纤维中克隆C4H1基因的cDNA全长,并利用DNAMAN、MEGA3.1等软件进行C4H1基因的生物信息学分析;采用半定量和实时荧光定量PCR方法,比较海岛棉和陆地棉不同发育时期的棉纤维中C4H1基因的表达水平.[结果]C4H1基因在海岛棉与陆地棉纤维发育不同时期表达具有差异,海岛棉中该基因的表达高峰期在20 DPA,而陆地棉中该基因的表达高峰期在15 DPA.该基因的核酸序列在两种棉花中存在差异,并导致其编码蛋白的氨基酸序列也存在差异,且海岛棉C4H1蛋白中Thr(苏氨酸)含量高于陆地棉,蛋白等电点也耍高于陆地棉.[结论]海岛棉中C4H1基因的表达高峰期晚于陆地棉,两个棉种C4H1蛋白氨基酸序列的差异可能导致其理化性质及生物学功能的改变,为发掘海岛棉纤维品质优异基因奠定了基础.     

大豆LEA基因家族全基因组鉴定、分类和表达

 【目的】鉴定大豆全基因组中的LEA家族基因,对其进行定位、分类以及组织表达分析,为植物LEA基因的功能研究与利用提供基础。【方法】利用大豆基因组数据库,通过生物信息学手段,鉴定并获得大豆LEA家族基因的全序列、定位和拷贝数信息;通过序列比对进行进化和分类分析;利用大豆发育表达芯片数据、NCBI中UniGene的EST表达数据进行组织表达谱分析。【结果】系统地分析鉴定了36个大豆的LEA家族基因,根据结构域的差异和系统发育分析将这些LEA基因分成8个亚家族。基因定位分析结果表明,这些基因分布于大豆的16条染色体上,启动子分析表明,几乎全部LEA基因的启动子区含有逆境反应顺式作用元件。各个发育阶段表达谱的分析结果表明,多数成员至少在一个组织中表达,10个差异表达的基因中有5个在种子发育时期优势表达,另外5个在其它部位优势表达。【结论】通过全基因组扫描,获得大豆基因组的36个LEA家族基因,它们分属于8个不同的亚家族,分布于16条大豆染色体上,启动子区含有逆境相关顺式作用元件,基因表达具有一定特异性。     

陆地棉CRK基因家族的鉴定及其表达分析

【目的】富含半胱氨酸类受体激酶（CRK）是植物中最大的类受体激酶家族之一,在植物生长发育、激素信号传导和抗逆境胁迫中发挥重要作用。从全基因组水平鉴定陆地棉CRK基因家族并进行生物信息学和表达模式分析,为研究和利用陆地棉CRK基因家族奠定基础。【方法】从Pfam数据库下载stress-antifung结构域氨基酸序列,应用BLASTp程序搜索棉花基因组数据库,鉴定棉花CRK基因家族;利用Compute pI/Mw tool、SignalP、TMHMM Server V2.0、WoLF POSRT等在线工具预测陆地棉CRK家族蛋白的分子量、信号肽、跨膜结构域和亚细胞定位等;用ClustalX1.8软件对棉花和拟南芥CRK蛋白质进行氨基酸序列比对,MEGA5.0分析棉花和拟南芥CRK蛋白的系统进化关系;使用TBtools制作陆地棉CRK基因家族的染色体定位、基因结构和蛋白质结构域示意图;应用植物顺式调控元件数据库PlantCARE分析棉花启动子序列;通过植物磷酸化位点数据库PlantPhos预测陆地棉CRK家族蛋白的磷酸化位点;从NCBI数据库下载RNA-Seq数据,利用转录组定量工具Kallisto计算TPM值,通过在线工具Morpheus绘制陆地棉CRK家族基因表达热图。【结果】陆地棉基因组中有70个CRK基因,分布于14条染色体,其中52个基因（74.3%）集中串联成簇分布于A6/D6、A9/D9、A10/D10染色体,且在A/D染色体组之间呈现高度共线性关系。编码302-901个氨基酸,58个蛋白质（82.9%）具有跨膜结构域,主要定位于叶绿体、质膜和胞外。磷酸化位点预测结果表明,陆地棉和拟南芥CRK有5个相同的磷酸化位点基序,包括3种丝氨酸磷酸化位点基序和2种苏氨酸磷酸化位点基序。65个陆地棉CRK基因的启动子区（92.9%）至少含有一种逆境激素响应元件,69个基因启动子区（98.6%）至少含有一种生物或非生物胁迫响应元件。根据RNA-Seq数据分析结果,陆地棉CRK基因可分为3种不同的组织表达特征类型;盐、干旱、冷、热胁迫以及接种大丽轮枝菌均可以导致部分陆地棉CRK基因表达水平的改变。GhCRK25在根、茎、叶和胚珠中优势表达,在纤维中几乎不表达,ABA、GA3、SA、PEG-6000、氯化钠和大丽轮枝菌Vd991处理均能刺激GhCRK25迅速上调表达。应用病毒诱导的基因沉默技术（VIGS）沉默GhCRK25可导致棉花对大丽轮枝菌Vd991更为敏感。【结论】陆地棉CRK基因家族有70个成员,具有保守的基因结构和功能结构域,多样化的组织表达特征,大多数基因受激素和逆境调控。     

番茄LBD基因家族的全基因组序列鉴定及其进化和表达分析

【目的】从番茄全基因组中鉴定LBD基因,并进行基因进化、基因结构、染色体定位以及组织表达和诱导表达分析,为番茄LBD基因的功能研究与利用奠定基础。【方法】利用番茄基因组数据库,通过生物信息学手段,鉴定番茄LBD家族成员;采用MEGA5软件进行系统进化树分析;通过perl程序、MapDRAW及GSDS工具进行基因结构及染色体定位分析;利用已有的番茄芯片数据进行组织表达谱和基因表达响应分析。【结果】系统分析鉴定了 46 个番茄LBD 家族基因,根据基因结构及系统进化分析将其分成class I与class II两类,细分为5个亚家族（Ia、Ib、Ic、Id与II）。基因定位表明,12 条染色体中的 10 条均有LBD 基因,该基因家族的分布具有广泛性。各个发育阶段表达谱的分析结果表明,LBD 家族基因具有不同的组织表达模式,在各个发育时期均有LBD 基因的表达。基因响应分析发现,不同LBD基因响应不同的外界信号。【结论】通过全基因组分析,番茄LBD家族基因包括 46 个成员,在进化上分为两大类,5 个亚家族,分布于 10 条染色体上,组织表达模式及基因响应具有多样性。这些信息为番茄LBD基因家族的功能分析奠定了基础。