陆地棉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启动子区包含了许多与植物激素反应、生长发育和逆境应答密切相关的顺式调控因子。通过对转录组数据分析以及荧光定量...     

棉花β-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基因，且在四倍体棉种中的数目约是二倍体棉种数...     

大白菜COBRA基因家族鉴定及其低温胁迫表达分析

【目的】鉴定大白菜COBRA基因家族成员,分析基因在低温胁迫下的表达并预测基因功能。【方法】利用生物信息学方法对大白菜的COBRA基因家族进行鉴定,并对基因结构、系统进化、启动子元件和表达模式进行分析。【结果】大白菜COBRA基因家族包含23个成员,命名为BrCOBL1～BrCOBL23,分布于大白菜基因组的8条染色体上,根据基因序列特点将其分为2个亚族。大白菜与同属于十字花科的拟南芥、甘蓝型油菜、甘蓝和芥菜基因组中共鉴定出147个COBRA基因,系统进化树分析表明以上COBRA基因家族成员可聚为8个亚组。基因表达分析显示：低温胁迫下,大白菜BrCOBL1、BrCOBL2、BrCOBL15和BrCOBL16的表达量显著升高。【结论】COBRA蛋白序列含有共同的保守基序,且某些基序特异存在于不同的亚族,表明这些基序可能与各亚族的功能分化有关。根据启动子分析结果推测BrCOBL基因可能参与大白菜的非生物胁迫响应及激素应答调控。     

小麦超氧化物歧化酶基因家族鉴定及盐胁迫下响应锌钾的表达分析

超氧化物歧化酶(SOD)是抗氧化系统的关键酶，在保护植物免受各种生物和非生物胁迫方面发挥着至关重要的作用。然而关于小麦(Triticum aestivum L.)SOD基因家族对盐胁迫的表达模式及响应锌(Zn)和钾(K)的功能知之甚少。本研究以‘西农979’为试验材料，基于全基因组分析，对小麦SOD家族成员进行鉴定和生物信息学分析，并分析相关基因在盐胁迫下响应Zn、K的表达情况。结果表明，在小麦中共鉴定出10个SOD基因(TaSDs),包括4个Cu/Zn-SOD基因(TaCSD1～TaCSD4)、4个Fe-SOD基因(TaFSD1～TaFSD4)和2个Mn-SOD基因(TaMSD1、TaMSD2)。系统发育分析表明，TaSDs可分为Cu/Zn-SODs、Fe/Mn-SODs两大类，共7个亚组。基因结构、序列基序及启动子分析表明，所有TaSDs都具有1～7个内含子和2～7个外显子，不同TaSDs具有不同的氨基酸序列组成和高度保守的活性位点残基，且每个TaSDs的启动子中存在许多与应激反应和植物激素相关的顺式作用元件。转录组表达谱分析表明，不同TaSDs对盐分表现出不同的表达模式，TaCS...     

黄瓜光敏色素互作因子的鉴定及表达分析

【目的】对黄瓜光敏色素互作因子（Phytochrome-interacting factor,PIF）家族基因进行系统分析和表达特征分析,为探究其在黄瓜生长发育与胁迫响应中的作用机理提供参考。【方法】基于黄瓜全基因组数据,通过生物信息学方法鉴定黄瓜PIF家族基因,并对家族基因的外显子—内含子结构、染色体定位、基因复制事件、系统进化树、保守基序及组织和胁迫相关表达模式等进行系统分析。【结果】从黄瓜基因组中鉴定出6个PIF基因,其不均等地分布于黄瓜的6条染色体上。基因复制事件分析结果表明,黄瓜PIF基因存在一个片段重复事件（CsPIF1a和CsPIF1b）。系统发育分析表明,来自黄瓜、拟南芥、辣椒、玉米和水稻的PIF蛋白可分为4个不同的组（PIF1、PIF2/3/6、PIF4/5、PIF7/8）。保守基序分析揭示了许多保守基序的存在,与PIF蛋白的分类一致。基因结构分析表明,CsPIF基因有5或6个内含子,与其他植物相比,内含子数量相对稳定。CsPIF基因在黄瓜各组织均有一定的表达,但在不同组织中表达量差异明显,CsPIF1a、CsPIF1b和CsPIF7分别在雌花、叶和卷须中表达量最高,...     

陆地棉膜结合脂肪酸去饱和酶基因家族的全基因组鉴定及表达分析

膜结合脂肪酸去饱和酶(fatty acid desaturase, FAD)是生物合成不饱和脂肪酸的关键酶。本研究以陆地棉基因组数据为基础，利用生物信息学方法对棉花FAD基因家族进行全基因组鉴定和进化分析。结果表明，陆地棉基因组中共含有37个FAD基因，进化分析发现这些基因分为4个亚家族，各亚家族成员数量不一，但相同亚家族成员含有相似的基因结构和保守基序。共线性分析发现28对复制基因全为片段复制基因，且都经历了严格的纯化选择作用。此外，在陆地棉FAD基因的启动子区域，发现了丰富的响应植物激素信号和逆境胁迫的顺式作用元件。转录组分析表明，陆地棉FAD基因家族响应干旱和盐胁迫，定量PCR分析表明施加外源褪黑素显著影响了FAD基因的表达。本研究结果为进一步解析FAD家族基因的功能奠定了基础。     

粉蕉PAL家族基因的鉴定及其在逆境胁迫下的表达分析

苯丙氨酸解氨酶(PAL)在植物生长发育和对生物/非生物胁迫的响应中起着重要作用，因此有必要对香蕉中的PAL家族基因进行详细的研究。本研究从粉蕉基因组中鉴定出8个MaPAL基因，分别命名为MaPAL1～MaPAL8,在系统发育树中MaPAL6基因单独被分成一支。根据基因系统进化关系、蛋白质结构域和基因结构的分析，8个MaPAL基因属于PAL家族基因。转录组分析结果表明，MaPAL基因参与了粉蕉的发育、成熟以及对生物/非生物胁迫的响应。MaPAL4基因在所有果实发育和成熟阶段都表现出高表达。MaPAL1和MaPAL7基因可能对非生物胁迫有响应，MaPAL2和MaPAL8基因的表达在香蕉枯萎病菌4号生理小种(Foc4)侵染时被明显抑制，可能响应Foc4的侵染。和巴西蕉相比，MaPAL基因在粉蕉中被低温和渗透胁迫诱导的程度要高。本研究结果为进一步研究香蕉PAL家族基因的功能提供了基础，也为香蕉遗传改良提供了潜在基因资源。     

小麦HP基因家族鉴定和分析

组氨酸磷酸转运蛋白HP(histidine phosphotransfer proteins)在植物的生长发育调控及逆境胁迫应答中发挥重要作用。为理解HP基因在小麦基因组中的进化特征和功能,研究通过生物信息学分析鉴定普通小麦的HP基因家族成员,对其理化性质、进化特征、基因结构、顺式作用元件以及在逆境胁迫条件下的表达模式进行了分析。结果表明,在小麦基因组鉴定到31个HP基因,编码蛋白质序列包括116～200个氨基酸。通过和其他植物的HP蛋白比较以及对蛋白结构、基因结构和基序的分析,显示小麦HP家族基因序列具有保守性。顺式作用元件预测表明,HP基因具有光、植物激素、干旱、低温等非生物胁迫响应相关的启动子调控元件。热图分析表明,HP基因在应答非生物胁迫中表达模式存在多样性;qRT-PCR分析表明,TaHP5-6B基因在低磷胁迫下受到强烈的诱导表达。综上所述,小麦全基因组鉴定的HP家族基因是非生物胁迫响应的重要基因资源。     

闽楠bZIP基因家族的鉴定与表达分析

【目的】基于闽楠全基因组数据对bZIP(碱性亮氨酸拉链)基因家族进行鉴定,并研究其在闽楠根系水分胁迫下的作用,确定关键调控基因。为进一步探究闽楠根系抵御逆境胁迫的分子机制提供了参考。【方法】利用生物信息学方法进行系统进化、保守基序、基因结构、启动子顺式势作用元件以及蛋白互作分析,并通过qRT-PCR验证基因在水分胁迫下的表达模式。【结果】在闽楠基因组中共鉴定出52个bZIP转录因子,划分为10个亚家族;顺式作用元件分析表明,PbbZIP可能响应光照、生物与非生物胁迫以及生长发育等生物学过程;转录组数据分析显示,PbbZIP15在干旱/水涝胁迫下的表达量较高,PbbZIP40在干旱胁迫下的表达量较高;qRT-PCR验证结果与转录组数据基本一致。【结论】PbbZIP15和PbbZIP40是调控闽楠根系水分胁迫的2个关键基因。     

陆地棉BGLU基因家族成员的全基因组鉴定与表达分析

植物BGLU基因在生物和非生物胁迫防御中起着重要作用，利用生物信息学手段对陆地棉（Gossypium hirsutum） BGLU家族成员进行全基因组鉴定，对序列特征、保守域结构、染色体定位、启动子元件等家族特征进行了分析，并进一步结合转录组数据和qRT-PCR分析了BGLU基因家族的病菌胁迫响应模式。结果表明，共鉴定出53个陆地棉BGLU基因，根据系统发育进化关系分为6个亚族，部分基因在病菌胁迫下特异性表达。荧光定量和转录组测序结果表明，GhBGLU5、GhBGLU14、GhBGLU18、GhBGLU26、GhBGLU34在抗病棉花中的表达水平显著高于感病品种，推测这些基因正向调控棉花黄萎病抗性。上述结果为进一步分析棉花BGLU家族基因的功能以及分子机制提供一定的理论基础。     

油橄榄Ⅲ类过氧化物酶家族的进化和表达分析

【目的】Ⅲ类过氧化物酶（PRX）是一类植物特有的氧化还原酶,在植物生长发育和胁迫响应方面具有十分重要的作用。本研究拟探讨油橄榄Ⅲ类PRX基因家族的进化和表达模式,旨在为油橄榄分子育种提供参考。【方法】利用生物信息学方法鉴定并分析了油橄榄Ⅲ类PRX基因家族的系统发育关系、分组、染色体分布、复制基因、基序分布、基因结构图、顺式作用元件分布和在不同组织和不同非生物胁迫下的表达,并对部分成员基因进行了实时荧光定量PCR验证。【结果】（1）鉴定得到了106个OePRX基因,根据其与拟南芥和毛果杨PRX蛋白序列的系统进化关系分为了14个组。（2）OePRX基因不均匀地分布于23条染色体上;片段复制是该基因家族扩张的主要动力,且复制基因在进化过程中受到了较强的纯化选择;与拟南芥相比,油橄榄PRX基因与毛果杨PRX基因的亲缘关系更近。（3）同一组中OePRX蛋白的等电点、分子量、基序分布、基因结构、信号肽分布、在不同组织中的表达模式均比较保守;OePRX基因启动子中含有较多的生长发育元件和激素应答元件;在油橄榄干热胁迫和水涝胁迫下,38%的OePRX基因显著差异表达。【结论】油橄榄OePRX基因家族明...     

光皮桦AP2/ERF基因家族鉴定与表达分析

  目的  深入研究AP2/ERF基因家族在光皮桦Betula luminifera生长发育及环境胁迫响应中的生物学功能。  方法  利用光皮桦基因组数据,通过生物信息学方法开展AP2/ERF基因家族鉴定、基因特征、系统进化、基因结构、保守基序、顺式作用元件、蛋白互作和表达分析。  结果  在光皮桦基因组中共鉴定到77个AP2/ERF基因,其编码的蛋白理化性质存在差异,大多数蛋白(60个)的理论等电点小于7.0。系统进化分析显示:这77个AP2/ERF转录因子属于5个亚家族,其中ERF亚家族最大,包含34个成员。光皮桦AP2/ERF各亚家族间的基因结构存在较大差异,其中AP2亚家族成员均具有6~9个内含子,而DREB亚家族基因则没有内含子;但AP2/ERF各亚家族内的不同成员具有相似的保守基序类型和分布。同时,AP2/ERF基因启动子上都存在大量与激素、调节、胁迫响应及生长发育相关的顺式作用元件。此外,互作网络分析预测不同的光皮桦AP2/ERF亚家族蛋白间存在广泛的互作关系。进一步的表达分析显示:绝大多数光皮桦AP2/ERF基因(71个)的表达存在较强组织特异性,且在高温胁迫下,多数ERF或DREB基因的表达发生显著变化,表明ERF和DREB基因在高温胁迫应答中可能发挥着重要作用。  结论  通过生物信息学分析获得光皮桦77个AP2/ERF基因,分属于5个亚家族。不同亚家族基因具有相似的基因结构、保守基序等特征。基因启动子区含激素、胁迫响应等相关的作用元件。基因表达具有较强的组织特异性,且多数ERF或DREB基因对高温胁迫有明显响应。图6表1参39     

高粱GRF基因家族鉴定及在非生物胁迫下的表达分析

生长调控因子（growth regulation factor,GRF）是植物中特有的一类转录因子,在植物激素调控、生长发育以及胁迫应答等多个领域扮演重要角色。对高粱GRF基因家族成员进行筛选鉴定,得到8个家族基因,分析了这些基因的理化性质、保守蛋白序列、基因功能及其共线性,分析了这些基因在非生物胁迫下转录组表达差异。结果表明,GRF基因家族结构较为保守,启动子元件含有干旱、低温、茉莉酸等响应元件,推测该家族基因参与生物及非生物胁迫应答响应。通过盐胁迫和烯效唑处理下的转录组数据分析,得到一个与耐盐相关的基因（SORBI_3002G297800）和一个与激素相关的基因（SORBI_3006G203400）,为进一步探究这些基因在高粱应对非生物胁迫中的作用提供了重要参考。     

大豆LIM转录因子家族鉴定及盐胁迫响应分析

为探究大豆LIM转录因子家族基因在应答盐胁迫中的功能,利用生物信息学方法,根据大豆全基因组数据,共鉴定出21个LIM基因,命名为GmLIM01～GmLIM21,分析大豆LIM基因家族的染色体位置、蛋白质的理化性质、进化关系、保守基序以及对非生物逆境胁迫的响应。结果表明,这些GmLIM基因在大豆20条染色体中的14条染色体上分布不均,片段重复是大豆LIM家族扩大的主要原因。根据系统发育进化分析,大豆21个LIM基因可分为5个亚家族,与其他9个物种的LIM基因一起构建进化树,也可分为5个亚家族。大豆的21个GmLIM基因共包含10个基序,大部分LIM蛋白含有motif 1、motif 2和motif 4这3个保守基序。此外,顺式调控元件的分析表明,在GmLIM基因的启动子序列中,与光响应、生长发育、植物激素和非生物逆境胁迫应答相关的调控元件非常丰富。定量PCR结果表明,在盐胁迫处理后大豆叶和根中GmLIM07、GmLIM17基因的表达量分别在3和6 h达到峰值,后随着处理时间的加长表达量下降。综上GmLIM基因在大豆的盐胁迫应答过程中具有重要的调控作用。     

陆地棉trihelix转录因子基因响应非生物胁迫表达谱分析

【目的】对陆地棉trihelix转录因子基因在多种胁迫下的应答表达模式进行研究分析,为棉花抗逆相关基因克隆及阐明棉花抗逆分子机制奠定基础。【方法】根据已知24个trihelix转录因子EST序列设计引物,通过RT-PCR和qRT-PCR的方法分析棉花（Gossypium hirsutum L.）GhGT在高盐（200 mmol.L-1NaCl）、干旱、低温（4℃）等非生物胁迫和外源激素ABA（100 μmol.L-1）处理下的表达模式。【结果】12个基因响应高盐和干旱胁迫应答,7个基因响应冷胁迫应答,应答ABA处理的基因9个。此外,有3个基因（GhGT2、GhGT18和GhGT23）在4种处理下均有应答。【结论】GhGTs在陆地棉的非生物胁迫适应过程中可能具有重要的作用。     

甘薯全基因组WRKY转录因子的基因鉴定与逆境胁迫表达分析

【目的】对甘薯(Ipomoea batatas)全基因组WRKY转录因子进行基因鉴定与逆境胁迫表达分析,为甘薯WRKY基因的应用提供参考。【方法】在甘薯全基因组序列的基础上,应用生物信息学方法对WRKY基因进行挖掘,分析基因潜在复制关系、保守结构域、亚家族系统进化树、保守基序和抗逆表达模式,并对甘薯SPF1基因与IbWRKY基因进行了比较和分析。【结果】甘薯全基因组上共有82个WRKY基因,可分为Ⅰ、Ⅱ和Ⅲ3个类型,其中类型Ⅰ和类型Ⅲ的成员数量分别为18和5,类型Ⅱ可进一步分为Ⅱ-a、Ⅱ-b、Ⅱ-c、Ⅱ-d和Ⅱ-e 5个亚类,其成员数量分别为4,14,20,8和13。甘薯WRKY基因不均匀分布于15条染色体上,其中25对基因存在潜在复制关系。保守结构域分析发现,大部分WRKY的结构域高度保守,但个别基因存在保守结构域缺失现象。拟南芥与甘薯的WRKY转录因子在系统进化树上呈现按物种少量聚集的现象。在甘薯WRKY转录因子家族中共发现10个保守基序,包括含有WRKY七肽结构域的基序1和基序7、含有锌指结构域的基序2和基序3,其中基序3为Ⅰ类型WRKY转录因子N端WRKY结构域所特有。对逆境胁迫下的甘薯转录组数据进行分析发现,甘薯苗期在蔓割病菌(Fusarium oxysporum f. sp.batatas)侵染后共有19个WRKY基因差异表达,甘薯块根贮藏期在低温胁迫下有34个WRKY基因差异表达。蛋白序列比较和分析结果表明,甘薯SPF1的氨基酸序列与IbWRKY32的氨基酸序列一致度最高,为85.4%。【结论】甘薯全基因组中鉴定得到82个WRKY基因,其结构域较为保守,在蔓割病胁迫与低温胁迫条件下均存在差异表达。     

荔枝GRF基因家族的全基因组鉴定及表达分析

为揭示荔枝（Litchi chinensis）GFR基因的功能,对荔枝生长调控因子（LcGRF）基因家族进行了全基因组鉴定和分析,并研究其在荔枝中的表达模式。基于荔枝基因组数据库,使用生物信息学软件对LcGRFs家族成员进行全基因组鉴定,分析基本理化性质、染色体定位、基因结构、进化关系、蛋白保守基序、顺式作用元件和时空表达情况进行系统分析。共获得12个GRF基因,不均匀地分布在10条染色体上,内含子2～5个。蛋白保守基序分析发现,荔枝GRF蛋白均含有保守的motif1(WRC)和motif2(QLQ),进化分析LcGRF划分为5个亚家族,LcGRFs启动子上存在大量的光、植物激素、非生物胁迫响应以及生长发育相关的顺式作用元件。不同转录组表达模式结果显示,各个组织中呈现出多样化的表达特征,表明不同成员可能在荔枝不同生长发育过程中发挥调控作用,参与调节荔枝的生长发育。研究显示,荔枝GRF家族成员有12个,划分为5个亚家族。     

冬瓜 WRKY 基因家族鉴定及表达分析

【目的】WRKY转录因子家族在植物生长发育、抵御胁迫等过程起着至关重要的作用。挖掘参与冬瓜（Benincasa hispida Cogn.）非生物胁迫的WRKY基因家族成员,探究其生物学功能。【方法】基于冬瓜基因组数据库鉴定冬瓜WRKY成员,利用生物信息学的方法系统分析其WRKY基因家族成员的蛋白理化性质、系统发育、保守基序、顺式作用元件,通过qPCR分析冬瓜WRKY的表达情况。【结果】冬瓜WRKY基因家族有57个成员,氨基酸数目在126～650之间,相对分子质量在13.92～71.68 kD之间;蛋白等电点在4.61～9.69之间。根据系统进化树将该家族分为3个类群,第二类群又分为5个亚组。冬瓜WRKY外显子有2～6个。染色体定位分析发现,有56个基因定位在12条染色体上,1个基因未定位在染色体上。保守基序分析显示,Motif 1和Motif 3存在于56个基因,且同一类群具有类似的保守基序结构。同时,冬瓜WRKY基因的启动子上均含有应激反应相关元件、激素响应元件。WRKY在冬瓜不同组织和果实发育时期特异表达。RT-PCR结果表明非生物胁迫和激素处理冬瓜可诱导部分WRKY基因的表达。...     

CoCl2胁迫下青稞苯丙氨酸解氨酶基因的克隆与表达分析

为阐明青稞幼苗对钴胁迫抗逆基因的响应模式,以青稞‘昆仑15号’为材料,克隆青稞苯丙氨酸解氨酶（PAL）基因,采用电子同源法克隆青稞苯丙氨酸解氨酶基因cDNA序列,并分析相关抗逆基因的表达变化情况。结果表明,该基因开放阅读框（ORF）序列有2 142bp,编码713个氨基酸,命名为HvPAL,其GenBank登录号为MK695675。经分析HvPAL编码的氨基酸序列有3个高度保守结构域,包含典型的植物苯丙氨酸解氨酶活性中心的保守氨基酸序列,与禾本科植物大麦和小麦的氨基酸序列有较高的相似性（90%）。PAL氨基酸的分子系统进化树显示青稞属于禾本科植物。采用实时荧光定量PCR法分析,与对照相比,70mg/L Co~(2+)处理下青稞幼苗SOD基因表达为最高（P0.05),150mg/L Co~(2+)处理下PAL、GSH和P5CS基因的表达量达到最大值（P0.05),200mg/L Co~(2+)处理下POD基因表达量最大（P0.05）。综上,青稞幼苗中POD、PAL、GSH、P5CS和SOD基因协同表达来抵抗和适应Co~(2+)胁迫。     

陆地棉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个成员,具有保守的基因结构和功能结构域,多样化的组织表达特征,大多数基因受激素和逆境调控。