枳漆酶基因家族鉴定及其响应盐胁迫的表达分析

漆酶(LAC)是植物木质素生物合成中催化木质素单体聚合的关键酶，在调控植物生长发育和胁迫响应中发挥着重要作用。对枳漆酶基因家族进行鉴定和分析，并探究其在盐胁迫下的表达模式，可为进一步研究枳漆酶基因功能提供重要的参考信息。采用生物信息学手段鉴定枳漆酶基因家族成员，对其家族成员的理化性质、基因结构、系统进化关系、启动子顺式作用元件等进行分析，并通过qRT-PCR方法分析其盐胁迫表达模式。结果表明，枳基因组中共有20个LAC基因家族成员，其中18个分布在5条已知染色体上，2个分布在未知染色体上，被预测定位到细胞膜和细胞核中；共有6～14个外显子，5～13个内含子，9～11个motif;系统进化树分析结果显示，20个枳LAC基因家族成员与17个拟南芥LAC基因家族成员共分成7个亚组，20个枳LAC基因家族成员分布在其中的6组；20个枳LAC基因家族成员与拟南芥LAC基因间存在19对共线性关系；启动子区域含有24种顺式作用元件，其中厌氧诱导元件、干旱响应元件和茉莉酸甲酯响应元件数量最多；16个枳LAC基因在盐胁迫下显著上调表达，推测枳漆酶基因参与了盐胁迫响应。     

玉米DREB转录因子家族的全基因组鉴定与分析

借助生物信息学手段,从玉米DREB家族基因的鉴定、基因结构、顺式作用元件等多角度进行解析,以期了解其在抗逆途径中的潜在功能。结果表明：玉米中共有87个DREB家族成员,分为6个亚组,理化性质分析结果表明其编码的氨基酸序列长度为125～452个氨基酸,均编码亲水性蛋白;在玉米10条染色体上均有分布,在4号、5号染色体长臂处较为集中;仅有少数基因含有内含子,大多数成员只以CDS的形式存在且在各亚组之间保守;启动子区域除含有MBS等响应干旱胁迫的元件外,还含有响应赤霉素、ABA、MeJA等相关调控途径的元件;物种共线性分析结果表明,高粱与玉米DREB家族成员之间存在共线性特征,且部分基因呈现倍性关系;基于转录组数据的表达模式分析显示,在干旱、盐、干旱和盐双重胁迫处理下,仅有几个基因的表达水平显著升高,多数DREB家族基因以下调的表达方式参与胁迫响应。     

番茄脱落酸受体PYL基因家族全基因组鉴定及表达分析

脱落酸受体家族包含许多功能基因，其中PYR/PYL/RCARs(PYL)在植物生长发育和防御反应等方面发挥重要生物学功能。PYL基因表达可参与ABA通路调节，对非生物胁迫作出响应。PYL基因家族成员在其他植物中已被鉴定和研究，但在番茄中研究较少。研究首次在番茄中鉴定出20个PYL基因家族成员，分别命名为SlPYL1～SlPYL20。利用生物信息学方法，分析其理化性质、基因结构、系统发育关系、染色体定位、共线性关系和顺式作用元件等，通过顺式作用元件分析发现其含有胁迫响应元件、光响应相关元件和植物生长发育相关元件，预示PYL基因家族功能多样性。转录组数据分析表明，干旱和高温处理前后SlPYL基因表达模式不同，q RTPCR分析发现6个SlPYLs基因全部响应干旱和外源ABA处理，部分响应冷胁迫和盐胁迫。研究为提高番茄逆境抵抗能力研究提供新的候选基因。     

毛果杨ZHD家族全基因组水平鉴定及在干旱胁迫下的表达分析

  目的  对毛果杨Populus trichocarpa ZHD (PtrZHD)家族进行生物信息学以及干旱胁迫下表达特性分析,为研究PtrZHD在干旱胁迫中的功能提供参考。  方法  利用生物信息学方法从全基因组水平鉴定出毛果杨ZHD家族全部成员,并对其进化、理化性质、基因结构、保守基序、启动子顺式作用元件和表达特性进行分析。  结果  毛果杨ZHD家族包括21个成员,可分为7个亚家族;有8对同源基因,且非同义替换率(K_a)/同义替换率(K_s)值远小于1。该家族成员理化性质存在差异,但其结构较为保守,均含有Motif 1;启动子区含有数量不等的激素和非生物胁迫响应元件,不同基因之间响应元件的种类存在差异。在毛果杨PtrZHDs中,分别有1、7和13个基因在根、茎和叶组织中具有偏好性表达特征;PtrZHD家族成员对干旱胁迫的响应具有组织和时间表达特异性,在根、茎和叶部组织中各成员的表达量不同,但随着干旱胁迫时间的增加均呈先上升后下降的趋势。  结论  PtrZHD家族基因对干旱胁迫有不同程度的响应,可调控毛果杨对干旱胁迫的应答。图6表2参27     

玉米脱水素基因家族的鉴定与分析

脱水素(dehydrin,DHN)属于LEA蛋白第二家族成员,是一种植物中广泛存在的亲水性蛋白,在干旱、低温和高盐等非生物胁迫的环境中起到重要作用。通过生物信息学方法对玉米全基因组DHN家族成员进行了鉴定,并进一步对其系统发育、基因结构、染色体定位和基因复制以及表达模式进行了系统分析。结果显示,在玉米DHN基因家族中共有5个家族成员,多物种系统发育进化树分析、基因结构以及基序分析都表明DHN家族成员在进化上具有高度的保守性。基因复制和物种间微共线性分析表明,在5个玉米DHN基因中存在着1对片段复制基因(ZmDHN1-ZmDHN2),玉米、高粱和水稻3个物种间存在2对直系同源基因(ZmDHN2-Sb04g032250.1,ZmDHN2-Os02g44870.1)。通过转录组表达数据分析表明,玉米DHN家族基因在不同发育时期具有不同的组织表达模式;同时,诱导表达模式分析表明ZmDHN基因的表达受到盐和干旱胁迫的显著诱导。该研究结果将为进一步鉴定玉米DHN家族重要的基因成员并对其开展功能分析奠定基础。     

玉米GRAS基因家族的全基因组鉴定及生物信息学分析

GRAS基因家族是植物中广泛存在的一类转录因子,在植物生长发育、生物和非生物逆境胁迫、光信号和激素信号应答等多个过程中发挥重要作用。对玉米GRAS基因家族成员的理化性质、染色体定位、系统发育、顺式作用元件等特征进行了分析。结果表明,在玉米全基因组中共鉴定出49个ZmGRAS基因,不均匀地分布于1～10号染色体上,编码蛋白质理化性质差异较大,可能在不同的微环境下发挥作用。系统进化分析将GRAS蛋白分为8个亚家族,可能在调节自身生长发育、逆境应答等过程中具有重要作用。玉米GRAS基因家族的启动子区含有激素应答、光响应、胁迫应答等多种顺式作用元件,推测其可能响应激素、胁迫等多种信号。共线性分析显示,具有共线性关系的基因可能是染色体片段复制的结果,且属于同一亚家族,具有相似的结构和功能。研究结果为进一步研究玉米GRAS基因的功能和逆境胁迫响应机制提供了依据。     

水稻OsRIP基因家族的全基因组鉴定及其表达分析

水稻(Oryza sativa)中RIP基因家族编码典型的E3泛素连接酶，在生长发育、逆境响应过程中发挥重要作用。利用水稻基因组数据，采用生物信息学方法鉴定了水稻中的RIP基因家族成员，鉴定到6个RIP家族基因，命名为OsRIP1～6，其编码区长度在906～1086 bp之间，分布在5条染色体上，内含子数量为2～3个。水稻中RIP基因被分为3个亚家族，每个亚家族的内含子数量、基因结构、motif基本一致。蛋白结构域分析结果表明，在水稻中，RIP家族的所有成员均含有SINA和RING结构域。共线性分析显示，水稻与玉米存在最多的共线对，与拟南芥不存在共线对。对启动子顺式作用元件分析发现，RIP基因启动子区有多个非生物胁迫响应元件、光响应元件及激素响应元件。另外，还分析RIP基因家族在不同组织器官、不同发育时期及昼夜节律的表达模式，结果表明在根、茎及花器官中表达水平较高；OsRIP6在播种后83 d表达水平最高，而其他成员在播种后48～69 d的表达处于较高水平；昼夜表达的分析结果显示，在光照后10 h该基因家族成员的表达水平较高。在ABA、JA激素和干旱胁迫、盐胁迫处理下，OsRIP1～5...     

大豆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基因在大豆的盐胁迫应答过程中具有重要的调控作用。     

玉米CPK基因家族鉴定及干旱胁迫下表达分析

干旱严重影响玉米的产量和品质。钙依赖性蛋白激酶(Calcium-dependent protein kinase, CPK)是一种Ca2+结合蛋白，在植物非生物胁迫应答方面发挥重要作用。本研究对玉米37个ZmCPKs的染色体定位、系统进化、基因结构、保守结构域和理化性质进行了分析。其中，36个ZmCPK基因不均匀的分布在9条染色体上；通过系统进化分析将其分为4个亚组，同一亚组成员具有相似的基因结构和保守结构域；ZmCPKs启动子分析表明，ZmCPKs启动子区含有大量的逆境响应元件，包括MBS、DRE和ABRE等；qTeller数据库和转录组数据分析表明，ZmCPKs基因家族成员在营养生长和生殖生长期均有表达，其中，26个基因在干旱胁迫下表达量发生变化。综合启动子元件及转录组数据，共筛选出9个基因作为干旱响应候选基因。进一步qRT-PCR验证显示，这9个基因均受PEG诱导表达。本研究可为玉米耐旱性的遗传改良提供优异的基因资源。     

花生KNOX基因家族的全基因组鉴定及组织表达分析

为研究花生KNOX基因家族在生长发育以及非生物胁迫响应中的功能，本研究采用生物信息学手段在栽培种花生基因组水平上对KNOX基因家族成员进行鉴定，并对其理化性质、基因结构、蛋白质保守结构域、系统进化、启动子顺式作用元件以及组织与逆境下的表达模式进行分析。结果显示：花生基因组中共有26个KNOX基因家族成员，分布在17条染色体上，绝大多数编码酸性蛋白，均为亲水性蛋白，并全部定位于细胞核。系统发育分析可将花生KNOX家族分为ClassⅠ和ClassⅡ两个亚家族，并进一步分为4个亚组。启动子分析发现一些与生长发育、激素和胁迫响应相关的顺式作用调控元件。全生育期组织表达分析结果显示，ClassⅠ类基因表达组织特异性明显，主要在花、营养茎尖、生殖茎尖、果针等组织中高表达；ClassⅡ类基因时空表达更广泛。在7种非生物胁迫处理下，部分基因表达量变化较大，其中5个基因响应特定胁迫，RT-qPCR验证了2个特异性响应PEG胁迫的基因。本研究为进一步探究花生KNOX基因在生长发育、逆境响应中的功能奠定了基础。     

豌豆CBF基因家族的生物信息学及低温胁迫下的表达分析

【目的】为了更加全面地探究CBF基因家族在豌豆中的进化及生物学特征。【方法】通过生物信息学方法从基因家族鉴定、结构分析、蛋白基本理化性质、亚细胞定位、染色体定位、顺式作用元件、蛋白质相互作用与表达分析等方面，对豌豆PsCBF基因家族进行系统鉴定和分子特征分析。【结果】在豌豆中共鉴定出21个PsCBF基因，染色体定位显示18个PsCBF基因分布在1、2、4、6与7号染色体上，其余3个未组装到染色体上(PsCBF19～PsCBF21)。21个PsCBF基因编码蛋白的氨基酸数在123(PsCBF09)～349(PsCBF11)，且均属于疏水性蛋白。亚细胞定位预测显示18个PsCBF基因(85.71%)定位在细胞核上，二级结构主要以无规则延伸为主。蛋白基序显示同一类群的基因都具相似的保守基序特征且该家族成员包含1～7个保守基序；同一类群的基因也具有相似的基因结构，且该家族基因含1～2个外显子。PsCBF基因启动子区分析表明其含有多种顺势作用元件，在激素反应，组织特异表达以及逆境响应中起作用。蛋白聚类网络图显示其蛋白可能在豌豆应对非生物胁迫中起调节作用。此外，本研究发现低温胁迫后K_2<...     

小麦TaBTLs基因家族全基因组鉴定与表达模式分析

【目的】鉴定和分析小麦TaBTLs基因家族成员。【方法】通过生物信息学方法，分析小麦TaBTLs家族成员的数量、生化性质、基因结构、共线性关系、系统进化关系和启动子顺式作用元件；利用小麦RNA-Seq数据，分析TaBTL基因在不同组织器官和发育时期的表达模式以及对非生物胁迫的响应水平。【结果】小麦共有35个TaBTLs基因家族成员，所有成员均含RING-H2和BZF保守结构域，相对分子量介于29.64～45.16 ku之间。TaBTLs基因分布在除2A、2B和2D染色体以外的所有染色体上，且各成员之间存在大量重复事件。不同物种间BTLs基因分为7个亚族。TaBTLs基因在启动子区域存在大量的植物生长发育和逆境响应顺式作用元件；该基因家族成员在不同的组织和器官中广泛表达，并且受盐和干旱胁迫诱导。【结论】研究结果为阐明TaBTLs基因家族的进化关系和进一步研究其家族成员的功能提供理论依据和前期工作基础。     

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

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

甜瓜GRAS家族全基因组的鉴定与表达分析

为了研究甜瓜(Cucumis melo L.)GRAS基因家族的结构特征、进化模式和生长发育过程中的表达模式，对甜瓜GRAS全基因组的基本信息、染色体位置、共线性、系统发育树、基因结构、表达模式、顺式元件进行了全面分析。结果表明，从甜瓜中鉴定出37个GRAS基因；这些基因在甜瓜的11条染色体上呈不均匀分布，并在6、7、8号染色体上出现串联重复；甜瓜GRAS家族可分为HAM、LS、SCR、SCL3、DELLA、LISCL、SHR、PAT1 8个亚族；基因结构分析显示，含有1～3个内含子的基因有8个，没有内含子的基因有29个；启动子分析，鉴定到生长素、脱落酸、赤霉素、低温、干旱、光信号响应作用元件以及参与分生组织和胚乳表达的元件；表达模式分析结果表明，HAM、DELLA亚族的基因几乎在所有组织中都有表达，LS亚族在根和茎中表达含量较高，PAT1亚族在甜瓜的整个生长发育过程中发挥着重要作用。得出结论，甜瓜全基因组共获得GRAS基因37个，共可分为8个亚族，不同的亚族具有不同的基因结构、蛋白质结构和表达模式；不同亚族在甜瓜生长发育的不同阶段具有差异性表达，它们共同参与调控甜瓜的发育过程。     

番茄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相互作用,共同参与调控番茄植株的逆境胁迫。     

基于PEG胁迫表达谱发掘调控谷子干旱胁迫的泛素连接酶E3基因

[目的]探析谷子调控干旱胁迫的E3基因。[方法]本文基于二代测序技术获得勾勾母鸡咀(GG)和晋汾16(JF)在PEG胁迫下的转录表达谱数据,运用生物信息学检索筛选出11个干旱相关的E3基因,对其理化性质、功能域、启动子元件、表达模式以及与其它物种E3基因的系统进化分析。[结果]蛋白等电点介于4.78~8.88之间,功能域主要为ZF_UBR型、ZF_RING和U_BOX型三类,启动子序列包含干旱诱导响应元件MBS或脱落酸响应元件ABRE等,其中Si001665m和Si022170m与PUB22/23亲缘关系较近,Si029090m与PUB19亲缘关系较近。[结论]初步确定这些是调控谷子干旱胁迫的E3候选基因,可能受干旱或脱落酸诱导表达,这为进一步研究谷子干旱响应相关的泛素连接酶E3基因提供理论依据。     

甘蓝型油菜WOX基因家族的鉴定与表达分析

对甘蓝型油菜WOX基因家族进行研究，利用生物信息学方法鉴定到58个家族成员，对其理化性质、系统发育树、基因结构、保守结构域、染色体位置、顺式作用元件和表达模式等进行分析，结果表明：氨基酸长度为121~397，分子质量为13 962.77~44 157.46 u，等电点为5.23~9.63，亲水性均小于0，不稳定系数为42.12~75.63，亚细胞定位在细胞核和叶绿体；系统进化树将WOX基因分为3大类和9个亚支；相同亚支的成员motif和基因结构相似或相同；除 BnWOX1、 BnWOX3、 BnWOX4、 BnWOX15、 BnWOX32、 BnWOX42、 BnWOX43、 BnWOX51、 BnWOX53 基因外，其余 49 个BnWOX基因不均匀分布在19条染色体中的18条染色体上(ChrC06 除外)，有38对串联重复基因；基因上游序列中包含转录因子结合位点、逆境响应、光响应、激素响应、分生组织调控等顺式作用元件等；各组织器官表达分析表明，BnWOX基因在根、茎尖、种子、角果中表达量较高；对WOX基因进行定量分析发现，在冷胁迫和干旱胁迫下基因表达量有显著变化，表明BnWOX基因可能在油菜响应冷/旱胁迫的调控机制中发挥着重要作用。     

枳Rboh基因家族全基因组鉴定及表达分析

植物呼吸爆发氧化酶(Rboh)可调控活性氧的产生,在植物生长发育和胁迫响应中起着重要作用。基于已公布的枳全基因组数据,利用生物信息学方法对枳Rboh基因家族成员进行鉴定,并分析其理化特征、基因结构、保守域信息、系统进化关系、启动子顺式作用元件和组织表达模式。结果表明：枳基因组中共有7个Rboh家族成员,分布在4条染色体和未知染色体上,被预测定位到细胞膜上;共有12～15个外显子,4～5个保守结构域,8～10个motif;系统进化分析分成5个亚组,与拟南芥Rboh基因之间存在更多的共线性关系;枳Rboh基因启动子区域共有38种顺式作用元件,包括生长发育元件、激素响应元件和逆境响应元件3大类型;枳Rboh基因表达具有组织特异性。这一研究为进一步研究枳Rboh基因功能提供了重要参考信息。     

小麦IQM基因家族鉴定及非生物胁迫下表达分析

IQM基因是钙调素结合蛋白家族中的重要分支,在植物生长发育和应激反应中发挥重要作用。本研究利用生物信息学方法在小麦全基因组中鉴定出23个IQM基因家族成员,对其染色体位置、理化性质、系统进化关系、基因结构、蛋白保守结构域、启动子顺式作用元件和基因表达特性等进行了系统分析。结果表明,TaIQM基因家族成员随机分布在小麦18条染色体上,亚细胞定位结果显示所有基因均位于细胞核中,系统进化分析将其分为3个亚类,同一亚类间基因结构、蛋白保守结构域相似度较高,推测其功能相似;顺式作用元件分析表明,TaIQM基因家族成员启动子中含有多种与逆境胁迫及生长发育相关顺式作用元件;转录组数据分析显示,TaIQM基因家族成员在不同时期的根、茎、叶、穗和籽粒中表达量存在显著性差异;qRT-PCR分析显示,在小麦苗期地上部和地下部,TaIQM基因响应干旱、低温、高温、NaCl、ABA等多种胁迫,显示上调或下调表达。初步推断这些基因可能通过Ca²⁺信号通路参与非生物胁迫调控,研究结果为全面解析TaIQM基因结构与生物学功能、非生物胁迫响应分子机制提供依据。     

毛竹磷转运蛋白Ⅰ家族基因鉴定及表达模式

  目的  鉴定毛竹Phyllostachys edulis磷转运蛋白Ⅰ (phosphate transporter 1, PHTⅠ)家族基因,分析其表达模式。  方法  利用生物信息学方法,鉴定毛竹PHTⅠ家族成员,分析基因启动子调控元件、编码蛋白的理化性质、基因结构、氨基酸保守基序、基因在染色体的位置、组织表达特异性、基因适应性进化及系统进化等。  结果  毛竹中共鉴定出20个PHTⅠ家族基因 (PePHTs),分布在10条染色体上,均定位于细胞膜。每个基因都含有1~2个内含子,PePHTs启动子序列中包含干旱、低温等非生物胁迫以及赤霉素等激素类响应元件。毛竹PHTⅠ大部分为碱性蛋白质,分子量为48.61~76.37 kDa,理论等电点为6.84~9.30,疏水性值均大于0,都属于疏水蛋白。基因适应性进化分析显示:大多数PePHTs的选择压力值小于0,说明多数基因受到负选择压力。转录组表达图谱表明:PePHTs在不同组织中的表达存在差异性,说明该基因家族在毛竹生长发育过程中发挥着不同的作用。系统进化树表明:PePHTs都聚类在第Ⅰ亚家族,并且优先和水稻Oryza sativa聚类在同一支上。  结论  PHTⅠ家族在植物吸收和转运磷的过程中扮演了重要作用。本研究结果为深入研究毛竹PHTⅠ家族基因的功能奠定了理论基础。图5表1参45