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

【目的】苯丙氨酸解氨酶是苯丙烷类代谢的关键酶,对调节植物生长和抗逆有重要作用。从全基因组层次了解棉花苯丙氨酸解氨酶基因家族的进化和功能,可为相关研究提供数据支撑。【方法】利用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胁迫下有不同的表达方式。【结论】本研究明确了苯丙氨酸解氨酶在棉花基因组中的数量、理化性质、进化和表达特征,为苯丙氨酸解氨酶的深入研究提供了参考。     

陆地棉PPO基因全基因组鉴定及对黄萎病菌的响应分析

【目的】多酚氧化酶(Polyphenol oxidases,PPO)广泛参与植物抵抗病虫害等生物逆境胁迫的过程。本研究旨在研究棉花中的PPO基因及其表达模式,验证多酚氧化酶与棉花抗黄萎病的关系。【方法】分析了黄萎病菌V991侵染下异源四倍体陆地棉TM-1根系中PPO酶活力变化,并利用TM-1的基因组数据库,鉴定相关基因,并进行生物信息学和表达分析。【结果】共筛选到13个候选GhPPO基因。这些基因都不含内含子,所编码的蛋白包含2个铜离子结合位点和PPO的保守序列(酪氨酸酶、DWL和KFDV保守结构域)。进化分析显示,陆地棉PPO基因的种内相似性大于种间相似性,并且相互之间形成了多对重复基因。GhPPO的表达量差异较大,少数基因具有组织特异表达的特性,多数基因在棉花不同组织中表达量都很低。实时荧光定量聚合酶链式反应分析结果表明:GhPPO6D在棉花根系中表达量较高,且在黄萎病菌V991侵染后上调。【结论】GhPPO6D可能参与了棉花与黄萎病菌的互作过程,与V991侵染后棉花根系PPO酶活力上升相关。     

陆地棉MIKC~C基因家族的全基因组分析

【目的】MIKC~C是1类保守的转录因子家族,参与调控植物的开花时间和花器官发育。通过对MIKC~C家族进行全基因组学分析,为深入研究MIKC~C在棉花开花及花器官发育的分子调控机理提供基础。【方法】利用HMMER 3.0及pfam种子文件鉴定棉花全基因组MIKC~C基因,结合其表达量进行偏向表达及聚类分析。【结果】在陆地棉基因组中共鉴定发现100个MIKC~C基因,分为11个亚类。表达聚类分析显示,棉花MIKC~C基因的表达模式大致可分为4个不同的类群,说明这类基因在棉花进化中出现了功能分化。100个MIKC~C基因中,有12个基因具有miRNA靶位点。【结论】研究结果显示MIKC~C家族基因在棉花纤维中存在功能分化,而且可能受到miRNA的调控,这为进一步研究该家族基因的功能提供信息参考。     

陆地棉抗坏血酸过氧化物酶基因家族全基因组生物信息学分析

【目的】抗坏血酸过氧化物酶(Ascorbate peroxidase,APX)是过氧化物酶家族中的1类成员,可以清除活性氧。APX家族的生物信息学分析有助于更好地了解该家族基因。【方法】利用生物信息学方法,从中国农业科学院棉花研究所提供的陆地棉序列中鉴定出14条APX基因家族序列,并对这14个成员的染色体定位、亚细胞定位、进化关系、理化性质、信号肽、保守基序以及其编码蛋白的结构特征进行了预测分析,并初步分析了部分基因的表达。【结果】APX基因分布广泛,分散在不同的染色体上;且大部分APX蛋白定位在细胞质,是亲水性脂溶蛋白;有3个保守基序,保守性较强;基因起源和进化比较复杂;二级结构的主要元件为α-螺旋和无规则卷曲。根据表达分析推测,第1亚组和第2亚组的基因在棉纤维发育的不同时期起作用。【结论】这些结果有望为APX家族成员的结构与功能挖掘提供参考。     

陆地棉叶片发育相关基因GhRH39克隆与功能分析

【目的】研究陆地棉DEAD-box RNA解旋酶基因GhRH39在叶片发育中的功能。【方法】借助生物信息学方法分析GhRH39的基因结构和进化关系,利用实时定量聚合酶链反应PCR (Quantitative real-time polymerase chain reaction, qRT-PCR)分析该基因在棉花不同组织和叶片不同发育时期的表达情况,利用病毒诱导的基因沉默(Virus-induced gene silencing, VIGS)技术对该基因进行沉默。观察阳性植株表型,检测其光合色素含量变化,对阳性植株中叶绿体发育和光合色素合成相关基因进行表达量检测。【结果】Gh RH39编码620个氨基酸,且序列较为保守。qRT-PCR分析表明,该基因在根、茎、叶片、顶芽、花瓣、纤维中均有表达,在叶片中表达量较高,且在叶片中的表达量随着叶片发育进程而变化。利用VIGS技术成功降低了GhRH39的表达水平,基因沉默的阳性棉株出现失绿表型,叶绿素a、叶绿素b、类胡萝卜素3种光合色素含量均降低。阳性沉默棉株的叶绿体发育和光合色素合成相关基因表达量出现一定程度的下降。【结论】GhRH39基因影响棉花叶绿素和类胡萝卜素合成,同时影响叶绿体发育过程。     

棉花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基因家族的进化与功能,为下一步深入研究该基因家族在棉花生长发育和响应红光信号通路中的分子调控机理提供基础。     

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

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

棉花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基因后续研究提供理论基础。     

四倍体棉花GAE基因家族的鉴定及其在棉纤维发育中的表达分析

【目的】通过对陆地棉及海岛棉GAE基因家族进行全基因组分析,为深入研究GAE基因家族在棉纤维发育中的功能提供理论依据。【方法】利用生物信息学分析软件,对GAE基因家族进行全基因组鉴定,根据转录组数据分析在纤维发育过程中的表达规律。【结果】陆地棉GhGAE家族包含21个成员,海岛棉GbGAE家族有22个成员,分为3个亚组,多数Gh GAEs和Gb GAEs基因没有内含子,共分布在12条染色体上。GAE氨基酸序列保守基序有4个,保守性较强,且GAE蛋白均定位于高尔基体膜。根据GAEs在陆地棉、海岛棉纤维发育不同时期的表达变化,将其分为纤维起始期高表达、纤维伸长期高表达、次生壁增厚期高表达、全时期低表达4类模式。其中GhGAE01、GhGAE02、Gh GAE11、Gh GAE12在陆地棉中高水平表达,GbGAE01、Gb GAE02、GbGAE11、Gb GAE12在海岛棉中高水平表达,推测这些基因可能在纤维发育过程中发挥着重要作用。【结论】上述结果为研究GAE基因家族在棉纤维发育中的功能提供了参考依据。     

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

【目的】类甜蛋白(Thaumatin-like proteins,TLP)参与多种植物病原真菌侵染后的防御反应。本研究通过在全基因组分析TLP基因,为深入探究其在介导棉花抗黄萎病通路中的分子机理提供基础。【方法】分别利用生物信息学方法和荧光定量聚合酶链式反应技术对陆地棉中的TLP基因进行全基因组鉴定和响应棉花黄萎病菌侵染表达分析。【结果】共鉴定出88个TLP基因;通过系统发育分析和序列特征分析可将棉花TLP家族分为10类。TLP基因家族外显子数量介于1～5之间,内含子介于0～4之间,同组成员具有相似的基因结构。大部分TLP含有5个保守的motif,具有相同的基序组织模式,均为Motif 5_4_2_3_1。染色体定位显示87个TLP基因定位于陆地棉20条染色体上,其中A亚组和D亚组各分布有42和45个,其在染色体上的分布存在串联重复现象。响应病原菌侵染表达分析显示,黄萎病菌侵染可以诱导6个候选基因的表达,不同抗性品种相对表达量达到峰值的时间不同,且耐病品种(GZ-1)上述基因表达量比感病品种(86-1)有增高的趋势。【结论】该研究结果有助于了解棉花TLP基因家族的进化与功能。     

陆地棉苗期低温响应基因GhZAT10的克隆及功能研究

【目的】挖掘棉花耐冷相关基因,为培育棉花耐冷品种奠定基础。【方法】克隆陆地棉基因GhZAT10(Zinc finger of Arabidopsis thaliana 10),通过实时荧光定量聚合酶链反应(Quantitative real-time polymerase chain reaction,qRT-PCR)分析冷处理前后该基因在根、茎、叶的表达;通过生物信息学方法分析Gh ZAT10基因结构及其编码的蛋白质性质及进化关系;通过Gateway技术构建蛋白表达载体35S::GhZAT10-GFP进行亚细胞定位;利用病毒诱导的基因沉默技术研究低温胁迫下GhZAT10在棉花耐冷反应中的功能。【结果】GhZAT10基因开放阅读框长度为813 bp(base pair,碱基对),编码270个氨基酸;GhZAT10在根中表达量高于茎和叶,低温胁迫下在3种组织中均上调表达。GhZAT10蛋白无信号肽,不包含跨膜螺旋,且其活性与其磷酸化调控关系密切。陆地棉GhZAT10与可可、拟南芥、甜橙中的ZAT10蛋白亲缘关系较近;GhZAT10蛋白定位在细胞核;沉默GhZAT10基因使棉花的耐冷能力减弱。【结论】GhZAT10蛋白属于C_2H_2型锌指蛋白,在陆地棉冷胁迫响应信号途径中具有积极作用。     

棉花核酸外切酶基因GhWRN的克隆及功能验证

【目的】本研究旨在探索棉花GhWRN基因在生长发育中的功能。【方法】用生物信息学方法分析GhWRN基因的结构特征、进化关系及上游1 500 bp的启动子片段并对其功能进行初步分析,利用实时荧光定量-聚合酶链式反应(Quantitative real-time polymerase chain reaction,qRT-PCR)技术分析该基因在不同熟性棉花品种、不同组织及不同激素处理下的表达特性。构建过表达载体转化拟南芥,观察转基因拟南芥的表型。利用qRT-PCR分析拟南芥开花信号途径中标记基因的表达变化。【结果】以陆地棉cDNA为模板克隆得到GhWRN基因,生物信息学分析表明该基因编码产物含有1个保守的WRN-exo结构域,无内含子,不具有跨膜结构,为非分泌蛋白。qRT-PCR表明该基因在早熟棉花品种的二叶期开始上调表达,在植株各组织中均有表达,在雄蕊和苞片表达量较高;外源激素ABA、IAA处理后0.5 h该基因极显著上调表达。与野生型相比转基因拟南芥开花时间提前,莲座叶数量减少,qRT-PCR分析表明花发育相关关键基因AtSOC1、AtFT和AtFUL上调表达。【结论】GhWRN基因与促进植物开花相关,为进一步探究棉花的开花机制提供新线索。     

陆地棉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基因,对于揭示棉花纤维品质等重要性状的遗传调控机制及分子育种具有一定的理论意义和应用价值。     

基于表达谱分析陆地棉DUF642基因家族抗逆功能

【目的】DUF642(Domain of unknown function 642 genes)是1个未知功能基因家族,在植物逆境胁迫响应中发挥重要作用。本研究旨在分析其在棉花抗逆反应中的功能。【方法】基于陆地棉基因组数据,系统地鉴定DUF642基因家族。利用生物信息学方法分析基因结构、理化性质、亚细胞定位、亲缘关系、启动子顺式作用元件等。应用转录组数据和实时定量聚合酶链式反应技术分析该家族基因在不同组织和多种逆境(冷、热、干旱、盐和黄萎病菌)胁迫下的表达规律。【结果】陆地棉基因组中包含23个DUF642基因,分布于14条染色体和1条Scaffold;编码的蛋白含有1～2个保守的DUF642结构域,多定位在细胞膜。DUF642家族基因可分成4个亚组,在进化上相对保守。DUF642基因具有较为广泛的组织表达类型,其中根和叶中表达较高。通过启动子和胁迫表达分析,推测Gh DUF642-09参与棉花抗盐胁迫;GhDUF642-08和GhDUF642-19参与棉花抗黄萎病;GhDUF642-02、GhDUF642-14和Gh DUF642-17参与棉花生长发育和抗逆胁迫。【结论】本研究结果为进一步研究DUF642基因家族功能和棉花抗逆分子机制提供了重要的理论依据。     

陆地棉Dof基因家族的全基因组鉴定及分析

【目的】单锌指DNA结合蛋白(DNA binding with one finger, Dof)是1类植物特有的转录因子,在植物生长发育与非生物胁迫响应中发挥非常重要的作用。本研究旨在对陆地棉Dof转录因子家族进行全基因组分析,为深入研究Dof基因参与植物生长发育的调控机制提供参考。【方法】基于最新发布的陆地棉基因组数据,利用生物信息学手段对陆地棉Dof基因家族进行全基因组鉴定,并系统分析Dof基因家族成员的理化性质、序列特征、基因复制、系统进化和表达谱。【结果】陆地棉中一共鉴定出118个Dof基因。经聚类分析将其分成9个亚家族,同一亚家族的陆地棉Dof基因有相似的内含子/外显子和基序分布模式。基因复制分析表明,全基因组复制可能导致陆地棉Dof基因的扩增。陆地棉Dof基因启动子区域存在不同植物激素和逆境胁迫响应元件。转录组分析表明,陆地棉Dof基因在不同的组织、发育时期和逆境胁迫下聚为3个表达模式不同的类群,暗示了其功能多样性。【结论】陆地棉Dof基因家族的全基因组鉴定和分析,有助于了解其进化和功能,为后续研究提供重要参考。     

ghr-miR160负调控棉花叶片衰老机理研究

【目的】陆地棉(Gossypium hirsutum L.)叶片早衰会导致棉花减产降质。MicroRNA(miRNA)是1类长度约为21 nt、在转录后水平调节基因表达的内源非编码小分子RNA。研究ghr-miR160对叶片衰老的调控具有重要意义。【方法】对ghr-miR160进行生物信息学分析、棉花子叶表达模式分析及过表达载体转化拟南芥验证分析。【结果】对ghr-miR160及pre-miR160结构分析表明,miR160在进化过程中高度保守。qRT-PCR结果显示,ghr-miR160在生长35 d的早熟不早衰品种辽4086子叶中优势表达。农杆菌介导转化拟南芥发现,该基因可使拟南芥出现晚抽薹、晚开花和叶片衰老延迟现象。【结论】ghr-miR160可能对棉花叶片衰老起负调控作用,同时为进一步阐明其作用机制奠定了基础。     

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

【目的】REM (Reproductive meristem)基因家族编码一类转录因子在植物生长发育过程中发挥重要作用,但在棉花中未见报道。【方法】基于陆地棉基因组数据和公共数据库中的转录组数据,运用生物信息学方法对陆地棉REM基因家族成员进行系统鉴定,并对该基因家族的理化性质、基因结构、组织表达特异性以及胁迫条件下的表达规律等进行分析。【结果】陆地棉REM基因家族包含79个成员,分布在25条染色体上,可分为5个亚组。每个成员编码的蛋白质至少含有1个B3结构域,大部分定位在细胞核。基因上游2 kb序列中存在激素及胁迫响应等顺式作用元件。组织表达特性分析结果显示,大部分REM基因在胚珠或纤维中的表达量较其他组织更高。逆境胁迫下的表达分析显示,29个基因响应棉花冷、热、盐或干旱胁迫。对6个基因在纤维不同发育时期表达量进行分析,发现这些基因的表达与已发表的转录组数据基本一致。【结论】明确了REM基因家族在基因组中的分布特征、结构特征以及系统进化特征,根据转录组数据初步揭示了该家族基因在生长发育和抗逆中的功能。     

陆地棉蔗糖磷酸合成酶基因家族的鉴定及表达分析

【目的】蔗糖磷酸合成酶(Sucrose phosphate synthase,SPS)是调控植物蔗糖代谢合成途径的关键酶,在植物光合产物的积累与分配方面发挥着重要作用,然而棉花中SPS基因的系统研究尚很少开展。本研究旨在对陆地棉SPS基因进行全基因组鉴定,并对它们的表达特性进行系统分析。【方法】基于已公布的陆地棉基因组序列,利用生物信息学和荧光定量聚合酶链式反应(Quantitative real-time polymerase chain reaction,qRT-PCR)等方法对陆地棉SPS家族基因的蛋白结构、进化关系、基因结构特征、染色体定位、基因复制和表达特性进行分析。【结果】(1)在陆地棉基因组中,共鉴定到10个Gh SPS基因(Gh SPS1-Gh SPS10);(2)Gh SPS蛋白具有植物SPS家族特有的两个保守的蛋白结构域和3个相对保守的蛋白磷酸化位点;(3)进化分析表明,Gh SPS蛋白可聚为A、B和C共3个亚族,其中A亚族成员最多,包含6个GhSPS蛋白;(4)位于同一亚族的GhSPS基因具有相似的外显子-内含子分布模式,但是外显子/内含子数目在不同亚族间差异很大;(5)GhSPS基因均匀地分布在陆地棉A亚组和D亚组的5条染色体上,片段复制可能导致了GhSPS基因在陆地棉基因组中的扩增;(6)转录组分析表明,不同亚族GhSPS基因具有不同的组织表达模式,A亚族Gh SPS基因在被检测的各个组织均有较高的表达,B亚族GhSPS基因主要在叶片中高表达,C亚族Gh SPS基因主要在纤维、叶片和花瓣中高表达;(7)进一步荧光定量PCR分析表明,GhSPS4在叶片中表达量很高,GhSPS1在叶片和花瓣中表达量较高,Gh SPS7和Gh SPS10在被检测的各个组织均有较高的表达,该结果与转录组分析结果相对一致。【结论】陆地棉SPS基因家族包含10个成员,分布在5条染色体上,可分为3个亚族,不同亚族成员呈现出不同的表达模式,为后续深入解析陆地棉SPS家族基因的功能奠定理论基础。