棉花GR基因家族的全基因组鉴定及分析

【目的】谷胱甘肽还原酶基因(Glutathione reductase gene,GR)家族参与植物生长发育和非生物胁迫响应等生物进程,但其在棉花中的特性及功能尚不清楚。本研究通过在异源四倍体陆地棉(Gossypium hirsutum)、海岛棉(G. barbadense)及其可能的二倍体祖先种亚洲棉(G. arboreum)和雷蒙德氏棉(G. raimondii)中对GR基因家族全基因组鉴定与特性分析,分析棉种分化及异源四倍体棉花形成过程中GR基因的进化历程,探讨其在非生物胁迫响应中的作用,为后续相关研究提供理论基础。【方法】利用生物信息学方法鉴定陆地棉、海岛棉、雷蒙德氏棉、亚洲棉的GR基因家族成员;解析GR基因家族成员的理化性质、序列特征、染色体位置、系统发育及表达模式。【结果】共鉴定到18个GR基因家族成员,陆地棉、海岛棉、雷蒙德氏棉和亚洲棉中GR基因数目分别为6、6、3和3个。系统发育分析发现GR基因分为2个亚组,同一亚组基因具有相似的外显子数目和基因结构。对同源基因的非同义突变率(Ka)及同义突变率(Ks)分析发现Ka/Ks值均小于1,表明在进化过程中GR基因经历了较强的纯化选择作用。陆地棉GR基因的表达模式分析表明,所有的GR基因均积极响应胁迫环境,但在不同的非生物胁迫下,基因的表达模式有明显差别。【结论】本研究探讨了GR基因家族在亚洲棉、雷蒙德氏棉、海岛棉和陆地棉基因组中的进化及功能,可为棉花GR基因后续研究提供理论基础。     

雷蒙德氏棉MAPKKK基因家族全基因组筛选及其同源基因在陆地棉中表达分析

【目的】丝裂原活化蛋白质激酶激酶激酶(Mitogen-activated protein kinase kinase kinase,MAPKKK)家族在植物的胁迫反应和发育过程中起重要调控作用。本研究旨在筛选雷蒙德氏棉MAPKKK基因并分析其功能。【方法】以已鉴定的拟南芥MAPKKK蛋白序列为种子序列,在已发表的雷蒙德氏棉全基因组数据库中,通过本地BLAST以及Pfam和SMART鉴定雷蒙德氏棉MAPKKK基因家族成员;采用MEGA5、GSDS在线工具以及Mapchart进行进化树、基因结构及染色体定位分析;利用已有的陆地棉芯片数据进行响应逆境胁迫和纤维不同发育时期的表达谱分析。【结果】系统鉴定了114个雷蒙德氏棉MAPKKK家族基因,根据基因结构及进化树分析分为Raf、ZIK和MEKK三个亚家族。染色体定位表明,该基因家族广泛分布于13条染色体上,并存在基因复制。与最近公布的78个雷蒙德氏棉MAPKKK家族基因相比对,获得序列完全相同的基因47个。【结论】上述研究结果有助于了解雷蒙德氏棉MAPKKK基因家族的进化与功能,为后续研究棉花乃至棉属MAPKKK基因的功能奠定基础。     

棉花CLE多肽家族的全基因组鉴定与生物信息学分析

【目的】CLE(CLAVATA3/Embryo surrounding region-related)多肽是一类小分子分泌蛋白,在植物生长发育过程中,对于维持分生组织细胞增殖与分化间的平衡起到重要的信号转导和调控作用。CLE家族是迄今为止报道的最大的植物多肽分子家族,也是近十年来研究最为热门的植物多肽激素。尽管CLE基因家族在拟南芥和水稻等模式植物中取得较多的研究进展,但在棉花中有关CLE多肽基因家族的研究尚属空白。通过对CLE多肽家族的鉴定及演化研究,可以为进一步挖掘棉花中的功能多肽及其信号转导研究提供一定的理论基础。【方法】本文以亚洲棉(Gossypium arboreum L.)、雷蒙德氏棉(G. raimondii Ulbrich)、陆地棉(G. hirsutum L.)和海岛棉(G. barbadense L.)为研究对象,利用生物信息学方法和已鉴定的拟南芥CLE基因序列,在棉花基因组中进行同源比对,并对4个棉种中的CLE家族成员进行全基因组鉴定和分析。【结果】共得到148个棉花CLE候选基因。陆地棉和海岛棉中CLE基因数量都分别是亚洲棉和雷蒙德氏棉的2倍。聚类分析将所有CLE基因分为5个亚组。Ka/Ks分析表明大部分基因都经历了负选择作用,有11对基因经历了显著的正选择作用。转录组数据分析发现,除海岛棉中的Gb CLE39、GbCLE13、GbCLE43,陆地棉中Gh CLE34、Gh CLE9以及亚组棉中的GaCLE4外,大部分CLE基因在棉花组织中的表达量较低,这与多肽在植物体内含量较低相一致。CLE核心保守基序分析发现了12个棉花特有的多肽,其中有2个多肽来源于受到强烈正选择作用的基因,因此在后续研究中可以对其进行进一步深入分析。【结论】本文利用生物信息学、比较基因组学等方法,首次对棉花的CLE基因家族进行了鉴定和分析,对于进一步挖掘棉花中的功能多肽及其信号调控网络研究提供了一定的理论基础。     

亚洲棉和雷蒙德氏棉MATE基因家族生物信息学及其同源基因在陆地棉中的表达分析

多药和有毒化合物排出(MATE)蛋白家族是1个次级转运蛋白家族。为了更好地了解亚洲棉和雷蒙德氏棉中MATE蛋白的种类和数量,利用2种二倍体棉花基因组的氨基酸和c DNA数据库对MATE基因家族进行筛选,并分析亚洲棉和雷蒙德氏棉全基因组中MATE基因的种类和进化关系。结果显示,在亚洲棉基因组中初步鉴定了34个MATE基因,而在雷蒙德氏棉中筛选出42个MATE基因。基因结构和系统进化的比较分析表明,进化关系近的MATE基因在结构上基本是一致的。同时对陆地棉中克隆的Gh TT12及与其同一族的部分MATE基因进行荧光定量分析,推测其所在同一族的MATE基因功能可能与转运原花青素有关。这些结果为进一步深入分析棉花MATE基因的功能以及在原花青素转运中的作用提供了理论基础。     

棉花中HD-Zip Ⅳ家族基因的全基因组鉴定及特征分析

同源异型亮氨酸拉链蛋白(homeodomain-leucine zipper,HD-Zip)是高等植物特有的一类转录因子,分成Ⅰ~Ⅳ4个亚家族。HD-ZipⅣ转录因子主要参与表皮结构的发育调控,在其他作物中已经广泛研究,但是棉花中还未进行系统的分析。本研究中,我们分别对陆地棉、亚洲棉、雷蒙德氏棉进行全基因组的鉴定,分别得到15个、12个、14个HD-ZipⅣ基因。除一个雷蒙德氏棉基因发生丢失外,陆地棉中几乎保留了所有的亚洲棉、雷蒙德氏棉中的HD-ZipⅣ基因。分析表明这些基因编码的蛋白都有HD-ZipⅣ转录因子特征性的结构域,且进化树上关系近的成员拥有相似的基因结构和基序组成。大部分陆地棉HD-ZipⅣ基因经历了纯化选择,且有一些经历正向选择。同时,进一步分析得出结论:片段复制在陆地棉HD-ZipⅣ基因的演化中起重要作用。组织表达分析表明大部分HD-ZipⅣ基因主要在发育着的胚珠和纤维中表达,暗示该家族与棉花中种子形成、纤维发育有关。总之,本研究为棉花中HD-ZipⅣ基因的功能研究及应用提供了参考依据。     

棉花CML基因家族成员鉴定与功能分析

【目的】类钙调素(Calmodulin-like,CML)蛋白是1种重要的Ca2+传感器,在调节植物应激反应机制中起重要作用。对CML基因家族进行全基因组学分析,以便深入研究CML基因在棉花逆境胁迫下的作用。【方法】利用生物信息学的方法进行了棉花CML家族成员的鉴定。利用转录组数据和反转录聚合酶链式反应分析陆地棉(Gossypium hirsutum L.)CML基因在盐胁迫下的表达模式。利用病毒诱导基因沉默技术对Gh CML44-2基因进行沉默并进行功能验证。【结果】在陆地棉(Gossypium hirsutum L.)、雷蒙德氏棉(G. raimondii Ulbrich)、亚洲棉(G. arboreum L.)中分别获得了154个、74个、78个CML蛋白。进化树和保守基序分析表明,GhCML蛋白分为8个亚类,都含有保守的EF-hand结构域;在盐胁迫下,107个Gh CML基因在叶和根中有显著的差异表达,且在根和叶中的差异表达模式不同。启动子分析表明,这些基因启动子中存在多种不同的胁迫响应元件;利用病毒诱导基因沉默技术成功沉默Gh CML44-2的棉株相较于对照更加不耐盐。【结论】该研究结果有助于了解棉花CML基因家族的进化与功能,为后续研究其功能提供了一定的理论依据。     

雷蒙德氏棉和亚洲棉Alfin-like PHD finger家族的生物信息学分析

Alfin-like PHD finger是植物中特有的一类转录调控因子,在调控植物生长发育、非生物胁迫响应及抗病反应等过程中起重要作用。为全面了解Alfin-like PHD finger基因家族在雷蒙德氏棉和亚洲棉基因组中的数目、分布情况及家族各成员间的关系,进而研究和揭示Alfin-like PHD finger在棉花抵抗逆境胁迫中的作用。运用生物信息学的方法在雷蒙德氏棉(D5)和亚洲棉(A2)全基因组中分别鉴定出12个和10个Alfin-like PHD finger基因家族成员,并对家族成员的基因结构、染色体定位、保守结构域及进化关系等方面进行分析。结果表明,雷蒙德氏棉和亚洲棉中Alfin-like PHD finger转录因子具有高度保守的DUF3594结构域和PHD结构域,两个二倍体棉花中该家族基因成员之间为一对一关系;雷蒙德氏棉中的12个成员基因分布在5条染色体上,亚洲棉中的10个成员基因也分布在5条染色体上;亚细胞定位预测发现,该家族的基因分别定位于细胞核、叶绿体基质、过氧化物酶体和细胞质基质;根据结构域差异和系统发育分析结果,将Alfin-like PHD finger转录因子家族分为3个进化分支。上述结果为分离鉴定棉花中新的逆境胁迫响应基因提供参考。     

雷蒙德氏棉LBD基因家族的鉴定及进化表达分析

家族转录因子是后基因组时代的研究热点。植物中的家族转录因子参与调控了许多重要的生物学过程,包括形态建成、信号转导、环境应激反应。植物特有的(lateral organ boundaries domain,LBD)/ASL(ASYMMETRIC LEAVES2-LIKE)家族转录因子包含保守的类似锌指结构CX2CX6CX3C基序,在调控拟南芥、水稻、玉米等模式植物生长发育过程中起到至关重要的作用。然而,对于棉花LBD基因的功能目前我们还知之甚少。本研究从雷蒙德氏棉基因组中鉴定出66个LBD基因,不均匀的分布在13条染色体上。这些基因可分为ClassⅠ和ClassⅡ两大类,细分为Ⅰa、Ⅰb、Ⅰc、Ⅰd、Ⅱa和Ⅱb等6个亚类。该家族基因结构简单,外显子数目不超过3个。从雷蒙德氏棉LBD家族成员中共找出15段保守基序,并对(lateral organ boundaries,LOB)结构域的保守性进行了分析。基因表达模式分析揭示出雷蒙德氏棉LBD基因家族的表达具有一定时空特异性。本研究为今后雷蒙德氏棉LBD基因功能的验证奠定了基础。     

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

【目的】磷是植物三大必需矿质元素之一,对植物的生长发育至关重要。Pht1家族的磷酸盐转运蛋白在植物磷吸收和转运方面具有重要作用,然而关于该基因的系统研究工作尚很少开展。本研究旨在进行Pht1家族成员全基因组鉴定和表达模式分析。【方法】通过生物信息学的方法对陆地棉Pht1家族成员的基因结构、编码蛋白质的结构、染色体定位、基因复制和表达情况等进行全面分析。【结果】(1)在陆地棉的基因组中共鉴定到17个磷酸盐转运蛋白基因(GhPT),其中A亚组包含8个GhPT,D亚组包含9个GhPT;(2)棉花GhPT蛋白之间序列相似性很高,并均具有12个疏水的跨膜区域;(3)进化分析表明这些GhPT蛋白主要聚为2大组(GroupⅠ和GroupⅡ),位于同一组的GhPT的编码基因大部分具有相似的内含子/外显子分布模式;(4)17个GhPT基因不均匀分布在A亚组和D亚组中的5条染色体上,串联复制和片段复制可能导致GhPT基因在陆地棉中的扩增;(5)表达模式分析表明,GhPT6和GhPT14在根中表达量最高,且同时响应低磷和低钾胁迫诱导并上调表达。【结论】这些结果将有助于深入了解Pht1家族基因的功能及离子信号途径相互作用的分子机制。     

雷蒙德氏棉和亚洲棉TCP基因家族的生物信息学分析

TCP基因家族是植物中特有的一类转录因子家族,涉及了植物的多种生长途径以及各种生理生化反应的信号传导。为了更好地了解TCP基因家族在雷蒙德氏棉和亚洲棉基因组中的数量和分布情况等,通过生物信息学方法,在雷蒙德氏棉(D5)和亚洲棉(A2)全基因组中分别鉴定出37个TCP基因家族成员,并对TCP家族成员进行基因结构、染色体定位、保守域、进化关系等分析。结果表明,雷蒙德氏棉和亚洲棉全基因组分别编码37个TCP转录因子成员。雷蒙德氏棉30个成员基因分布在10条染色体上,亚洲棉37个成员基因分布在13条染色体上;内含子/外显子结构分析表明,TCP基因结构较为简单,大部分不含内含子;功能结构域分析,发现所有TCP转录因子具有高度保守的TCP结构域;根据结构域差异和系统发育分析的结果,将TCP基因家族分成2个亚族,进一步划分为PCF、CIN、CYC/TB1三个亚类。以上结果为进一步研究棉花TCP基因家族各成员的功能奠定一定的理论基础。     

棉花FAR1/FHY3基因家族的全基因组分析

【目的】FAR1/FHY3是1类起源于转座子酶的转录因子,在光敏色素A(phyA)信号通路的启动中起到非常重要的作用。通过对FAR1/FHY3家族进行全基因组学分析,为深入研究FAR1/FHY3在棉花光信号通路中的分子调控机理提供基础。【方法】利用生物信息学方法对亚洲棉(Gossypium arboreum L.)、雷蒙德氏棉(Gossypium raimondii Ulbrich)和陆地棉(Gossypium hirsutum L.)中的FAR1/FHY3基因进行了全基因组鉴定和分析。【结果】共鉴定出88个FAR1/FHY3基因,其中陆地棉中数量最少。聚类分析将所有成员分为3个亚组。Ka/Ks值表明大部分基因都经历了负选择过程。转录组数据分析发现,大部分FAR1/FHY3基因在棉花叶片中表达,少数在茎、花瓣和花托中表达。陆地棉中的Gh FAR14在4种组织中均有较高水平的表达,且与拟南芥FHY3具有较高的同源性,推测其在陆地棉phyA信号通路的启动过程中可能起到关键作用。【结论】该研究结果有助于了解棉花FAR1/FHY3基因家族的进化与功能,为下一步深入研究该基因家族在棉花生长发育和响应红光信号通路中的分子调控机理提供基础。     

陆地棉MADS-box家族基因鉴定及组织特异性表达分析

【目的】MADS-box蛋白是1种重要的转录因子,参与调控棉花生长发育进程。鉴定分析陆地棉MADS-box家族基因,可为研究其生物学功能奠定基础。【方法】基于2019年最新组装的陆地棉TM-1参考基因组,通过HMMER 3.0软件鉴定陆地棉TypeⅠ型及MIKC型MADS-box基因。利用MapInspect、MEGA7.0、MEME和TBtools软件以及omicshare网站进行染色体定位、多序列比对聚类分析、Motif预测、基因结构鉴定以及组织特异性表达分析。【结果】从全基因组水平上鉴定出181个陆地棉MADS-box基因(66个TypeⅠ型和115个MIKC型),染色体定位发现该181个基因在陆地棉26条染色体上均有分布,其中A组染色体上分布有78个基因、D组上分布有103个基因。系统进化分析显示TypeⅠ型MADS-box蛋白有3个亚家族,MIKC型蛋白包括MIKC*型和含有10个亚家族的MIKCC型;motif预测结果显示所有MADS-box蛋白均含有MADS结构域,基因结构分析结果显示,外显子、内含子结构和长度在各亚家族内具有同一性;组织特异性表达分析发现MADS-box家族基因中有33个在纤维组织中优势表达,103个在花器官中优势表达,41个基因在根茎中优势表达。【结论】本研究通过生物信息学方法,鉴定并获得了与陆地棉开花调控和纤维发育相关的MADS-box基因,对于揭示棉花纤维品质等重要性状的遗传调控机制及分子育种具有一定的理论意义和应用价值。     

棉花UGPase基因鉴定与生物信息学分析

【目的】尿苷二磷酸葡萄糖焦磷酸化酶(Uridine diphosphate glucose pyrophosphorylase,UGPase)是糖合成代谢途径中的一种重要限速酶,在植物纤维细胞发育中起重要作用,但是关于棉花中UGPase基因的生物信息学研究却很少。本研究采用生物信息学手段鉴定和分析棉花UGPase基因家族特性,为其相关研究提供重要线索。【方法】以陆地棉、海岛棉、亚洲棉和雷蒙德氏棉4种棉花为主要研究对象,结合藻类、水稻、葡萄、可可、榴莲等19种植物的UGPase基因进行系统进化分析和陆地棉UGPase基因的表达特性分析。【结果】系统发育分析显示UGPase基因被分为UGPase-A和UGPase-B两大类。保守基序及适应性进化分析显示,UGPase-B类基因在进化中十分保守,UGPase-A类基因在进化中与UGPase-B类基因存在较大差异。棉花、葡萄、可可、榴莲UGPase基因同源性分析显示,棉花中大多数UGPase基因来源于异源四倍体棉花形成之前的基因组加倍事件。UGPase-A类在棉纤维发育起始和伸长期表达量较高,UGPase-B类在营养器官中表达量较高。【结论】本研究明确了UGPase两类基因的结构、功能与进化,为后续研究奠定基础。     

陆地棉Nudix基因家族的全基因组鉴定及表达分析

【目的】Nudix是一类能够催化各种核苷二磷酸衍生物水解的酶,具有维持遗传物质稳定,响应逆境胁迫等生物学功能。本研究旨在对陆地棉Nudix基因进行全基因组分析,为深入研究Nudix基因家族参与棉纤维生长发育调控机制提供参考。【方法】通过生物信息学方法对陆地棉基因组(ZJU_v2.1)的Nudix基因进行全基因组鉴定,并系统分析Nudix基因家族成员的理化性质、序列特征、基因复制、系统进化和表达情况。【结果】从陆地棉基因组中共鉴定到76个GhNudix基因,经聚类分析将其分为7个亚组(CladeⅠ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、Ⅵ和Ⅶ)。基因复制分析表明,片段重复是陆地棉GhNudix基因家族成员扩增的主要原因。转录组数据表明GhNudix基因的表达不但有组织特异性而且有时间特异性。其中在纤维发育的起始和伸长阶段高调表达的GhNudix基因主要分布在CladeⅠ和Ⅱ;且大部分CladeⅠ和Ⅱ亚组的GhNudix基因启动子含有响应生长素和赤霉素的顺式作用元件,由此推测CladeⅠ和Ⅱ亚组的GhNudix基因在陆地棉纤维的起始和伸长阶段发挥作用。【结论】陆地棉GhNudix基因家族的全基因组鉴定和分析,为后续研究其在棉纤维发育中的作用机制提供参考。     

棉属D基因组3个野生棉种Bet v1基因鉴定及功能分析

【目的】对棉花D基因组中Bet v 1基因家族进行分析,比较其在不同抗性棉种间的表达模式差异,为深入研究Bet v 1基因在棉花抗黄萎病中的作用提供理论依据。【方法】通过生物信息学方法对D基因组中Bet v 1基因进行鉴定。通过雷蒙德氏棉(Gossypium raimondii,D5)、三裂棉(G. trilobum,D8)和瑟伯氏棉(G.thurberi,D1)转录组和实时荧光定量聚合酶链式反应(qRT-PCR)分析Bet v 1基因在黄萎病菌处理下的表达模式。利用病毒诱导的基因沉默技术(VIGS)对Bet v 1基因进行功能鉴定。【结果】[棉花D基因组中包含59个成员,其中57个基因带有内含子,分布于8条染色体上,多数为亲水性蛋白并定位于细胞质。]黄萎病菌胁迫条件下3个野生棉种的Bet v 1基因表达量与其抗病水平一致。将不同表达水平的基因分为3组,其中第3组基因响应黄萎病菌侵染并在抗病棉种瑟伯氏棉中高表达,表明该组基因可能与黄萎病胁迫应答反应有关。从中筛选出高水平表达的Bet v 1基因,在陆地棉中沉默相应的同源基因,棉株感病加重,揭示该基因在棉花抵御黄萎病菌侵染过程中起正调控作用。【结论】响应黄萎病菌胁迫的Bet v 1基因在棉花抗黄萎病复杂的生物过程中至关重要。本研究结果为棉花Bet v 1家族基因的深入研究提供依据,为进一步解析棉花Bet v 1基因的功能奠定基础。     

陆地棉Nudix基因家族的全基因组鉴定及表达分析

【目的】Nudix是一类能够催化各种核苷二磷酸衍生物水解的酶,具有维持遗传物质稳定,响应逆境胁迫等生物学功能。本研究旨在对陆地棉Nudix基因进行全基因组分析,为深入研究Nudix基因家族参与棉纤维生长发育调控机制提供参考。【方法】通过生物信息学方法对陆地棉基因组(ZJU＿v2.1)的Nudix基因进行全基因组鉴定,并系统分析Nudix基因家族成员的理化性质、序列特征、基因复制、系统进化和表达情况。【结果】从陆地棉基因组中共鉴定到76个Gh Nudix基因,经聚类分析将其分为7个亚组(CladeⅠ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、Ⅵ和Ⅶ)。基因复制分析表明,片段重复是陆地棉Gh Nudix基因家族成员扩增的主要原因。转录组数据表明GhNudix基因的表达不但有组织特异性而且有时间特异性。其中在纤维发育的起始和伸长阶段高调表达的GhNudix基因主要分布在CladeⅠ和Ⅱ;且大部分CladeⅠ和Ⅱ亚组的GhNudix基因启动子含有响应生长素和赤霉素的顺式作用元件,由此推测CladeⅠ和Ⅱ亚组的GhNudix基因在陆地棉纤维的起始和伸长阶段发挥作用。【结论】陆地棉GhNudix基因家族的全基因组鉴定和分析,为后续研究其在棉纤维发育中的作用机制提供参考。     

棉花苯丙氨酸解氨酶基因家族的生物信息学分析

【目的】苯丙氨酸解氨酶是苯丙烷类代谢的关键酶,对调节植物生长和抗逆有重要作用。从全基因组层次了解棉花苯丙氨酸解氨酶基因家族的进化和功能,可为相关研究提供数据支撑。【方法】利用HMMER和BLASTP从棉花基因组筛选苯丙氨酸解氨酶基因家族,利用MEGA-X进行进化分析,利用MCScanX进行共线性分析,并利用NCBI SRA的数据进行表达分析。【结果】本研究分别从陆地棉、海岛棉、亚洲棉和雷蒙德氏棉基因组中筛选到12、13、7、8个苯丙氨酸解氨酶基因,依据进化分析和共线性关系数据,我们推测四种棉花的二倍体祖先起初应有6个苯丙氨酸解氨酶基因。顺式作用元件分析结果表明,苯丙氨酸解氨酶基因家族含有较多的光响应、激素响应、胁迫响应和生长发育相关的元件。苯丙氨酸解氨酶的蛋白序列具有保守的基序。陆地棉的12个苯丙氨酸解氨酶在低温、高温、盐和PEG胁迫下有不同的表达方式。【结论】本研究明确了苯丙氨酸解氨酶在棉花基因组中的数量、理化性质、进化和表达特征,为苯丙氨酸解氨酶的深入研究提供了参考。     

棉花苯丙氨酸解氨酶基因家族的生物信息学分析

【目的】苯丙氨酸解氨酶是苯丙烷类代谢的关键酶,对调节植物生长和抗逆有重要作用。从全基因组层次了解棉花苯丙氨酸解氨酶基因家族的进化和功能,可为相关研究提供数据支撑。【方法】利用HMMER和BLASTP从棉花基因组筛选苯丙氨酸解氨酶基因家族,利用MEGA-X进行进化分析,利用MCScanX进行共线性分析,并利用NCBI SRA的数据进行表达分析。【结果】本研究分别从陆地棉、海岛棉、亚洲棉和雷蒙德氏棉基因组中筛选到12、13、7、8个苯丙氨酸解氨酶基因,依据进化分析和共线性关系数据,我们推测四种棉花的二倍体祖先起初应有6个苯丙氨酸解氨酶基因。顺式作用元件分析结果表明,苯丙氨酸解氨酶基因家族含有较多的光响应、激素响应、胁迫响应和生长发育相关的元件。苯丙氨酸解氨酶的蛋白序列具有保守的基序。陆地棉的12个苯丙氨酸解氨酶在低温、高温、盐和PEG胁迫下有不同的表达方式。【结论】本研究明确了苯丙氨酸解氨酶在棉花基因组中的数量、理化性质、进化和表达特征,为苯丙氨酸解氨酶的深入研究提供了参考。