陆地棉GhCNGC13与GhCNGC32基因功能的深度解析与鉴定

陆地棉GhCNGC13与GhCNGC32基因功能的深度解析与鉴定一、引言1.1研究背景与意义陆地棉（GossypiumhirsutumL.）作为世界上最重要的棉花栽培种，在全球农业经济中占据着举足轻重的地位。其种植面积广泛，产量在棉花生产中占主体地位，为纺织工业提供了主要的天然纤维来源，是关系到国计民生的重要经济作物。据统计，全球约80%的棉花由陆地棉品种提供，我国作为棉花生产和消费大国，陆地棉同样是我国的主要经济作物之一。随着纺织工业的不断发展和人们对高品质纺织品需求的增加，对棉花纤维品质提出了更高要求。同时，棉花在生长过程中面临着各种生物和非生物胁迫，如病虫害、干旱、盐碱、低温等，这些逆境因素严重影响棉花的生长发育、产量和品质。如何提高棉花的产量、改善纤维品质以及增强其抗逆性，成为棉花育种领域亟待解决的关键问题。基因功能鉴定是深入了解生物生长发育、代谢调控以及逆境响应机制的基础，对于作物品种改良具有重要的指导意义。在棉花研究中，鉴定与产量、品质和抗逆性相关的关键基因，并解析其功能和作用机制，能够为棉花分子育种提供理论依据和基因资源，从而加速棉花新品种的培育进程。环核苷酸门控离子通道（CyclicNucleotide-GatedIonChannels，CNGCs）是一类广泛存在于植物中的离子通道蛋白，在植物生长发育、信号转导以及对生物和非生物胁迫的响应中发挥着重要作用。GhCNGC13和GhCNGC32作为陆地棉中的CNGC基因家族成员，研究它们的功能对于揭示陆地棉生长发育及抗逆机制具有重要意义。通过对这两个基因的功能鉴定，有望深入了解它们在陆地棉响应逆境胁迫过程中的作用途径，为棉花抗逆育种提供新的基因靶点和理论支持。同时，也有助于进一步丰富植物CNGC基因功能的研究内容，完善植物抗逆分子机制的理论体系。1.2研究目的本研究旨在深入鉴定陆地棉GhCNGC13和GhCNGC32基因的功能，通过一系列分子生物学实验和分析手段，明确这两个基因在陆地棉生长发育、逆境响应等过程中的具体作用机制。主要研究目的包括：利用生物信息学方法对GhCNGC13和GhCNGC32基因的结构、序列特征以及进化关系进行全面分析，初步预测其功能；通过基因表达模式分析，探究这两个基因在陆地棉不同组织、不同发育时期以及在各种生物和非生物胁迫条件下的表达变化规律，明确其表达调控与陆地棉生长发育及逆境响应的相关性；运用病毒诱导的基因沉默（VIGS）技术或基因编辑技术，对陆地棉中的GhCNGC13和GhCNGC32基因进行功能缺失或功能获得性突变，分析突变体植株在生长发育、生理生化指标以及对逆境胁迫耐受性等方面的表型变化，从而验证基因的功能；进一步研究GhCNGC13和GhCNGC32基因参与的信号转导途径以及与其他相关基因的互作关系，揭示其在陆地棉响应逆境胁迫过程中的分子调控网络。通过以上研究，为棉花遗传改良提供理论基础，为培育高产、优质、抗逆性强的棉花新品种提供基因资源和技术支持，推动棉花产业的可持续发展。1.3国内外研究现状陆地棉作为重要的经济作物，其基因功能研究一直是国内外学者关注的焦点。随着分子生物学技术的飞速发展，对陆地棉基因功能的研究取得了显著进展。在产量相关基因研究方面，国内外学者通过数量性状位点（QTL）定位、全基因组关联分析（GWAS）等方法，鉴定出许多与陆地棉产量相关的基因位点。例如，[研究团队1]利用陆地棉重组自交系群体，定位到多个与籽棉产量、皮棉产量等性状相关的QTL，为产量性状的遗传改良提供了理论基础。[研究团队2]通过GWAS分析，挖掘出一些与棉铃大小、铃数等产量构成因素相关的候选基因，为进一步解析产量形成的分子机制提供了线索。在纤维品质相关基因研究中，众多研究聚焦于棉纤维发育的分子调控机制。[研究团队3]从陆地棉中克隆出多个与纤维伸长、加厚等过程相关的基因，如GhEXP1、GhMYB25等，发现这些基因通过调控细胞壁的合成与松弛，影响纤维的长度和强度。[研究团队4]利用转录组测序技术，分析了不同纤维品质陆地棉品种在纤维发育不同时期的基因表达谱，筛选出一批差异表达基因，为纤维品质改良提供了新的基因资源。关于陆地棉抗逆基因的研究也取得了丰硕成果。在应对生物胁迫方面，研究发现一些基因参与了棉花对病虫害的防御反应。如[研究团队5]鉴定出的抗黄萎病基因，通过激活植物的免疫信号通路，增强棉花对黄萎病菌的抗性。在非生物胁迫响应方面，[研究团队6]发现多个基因在棉花响应干旱、盐碱、低温等胁迫过程中发挥重要作用，这些基因通过调节渗透调节物质的合成、抗氧化酶的活性以及胁迫相关基因的表达，提高棉花的抗逆性。针对环核苷酸门控离子通道（CNGC）基因家族，在拟南芥、水稻等模式植物中的研究较为深入。研究表明，CNGC基因在植物的生长发育、离子稳态维持、信号转导以及对生物和非生物胁迫的响应中具有重要功能。例如，拟南芥中的AtCNGC1和AtCNGC2基因参与了植物对病原体的防御反应，通过调节钙离子内流，激活植物的免疫信号通路；AtCNGC10基因在维持植物体内钾离子平衡和响应低钾胁迫中发挥关键作用。然而，在陆地棉中，对CNGC基因家族的研究相对较少，尤其是关于GhCNGC13和GhCNGC32基因的功能研究，目前尚处于起步阶段。仅有少量研究初步分析了它们的序列特征和表达模式，但对于这两个基因在陆地棉生长发育及逆境响应过程中的具体功能和作用机制，仍缺乏深入系统的研究。综上所述，尽管陆地棉基因功能研究已取得一定成果，但在CNGC基因家族功能研究方面仍存在不足与空白。深入开展陆地棉GhCNGC13和GhCNGC32基因的功能鉴定，将有助于填补这一领域的研究空白，为棉花的遗传改良提供新的理论依据和基因资源。二、材料与方法2.1实验材料2.1.1植物材料本实验选用陆地棉品种‘中棉所49’作为研究材料。该品种是我国广泛种植的优良陆地棉品种，具有生长势强、产量高、纤维品质较好等特点，在棉花生产和研究中被广泛应用。棉花种子经消毒处理后，播种于装有营养土的塑料盆中，每盆播种3-5粒种子，置于光照培养箱中培养。培养条件为：光照强度120-150μmol・m⁻²・s⁻¹，光照时间16h/d，昼夜温度分别为30℃/25℃，相对湿度60%-70%。待棉花幼苗长至三叶一心期时，进行间苗，每盆保留1株生长健壮的幼苗。在棉花生长发育的不同时期，分别采集根、茎、叶、花、蕾、纤维等组织样本。采集的样本迅速放入液氮中速冻，然后转移至-80℃冰箱中保存，用于后续的RNA提取和基因表达分析。为了研究基因在逆境胁迫下的表达情况，对生长至五叶期的棉花幼苗分别进行干旱、盐胁迫、低温和高温处理。干旱处理采用PEG-6000模拟干旱胁迫，将棉花幼苗根部浸泡在含有20%PEG-6000的Hoagland营养液中；盐胁迫处理将幼苗根部浸泡在含有200mmol/LNaCl的Hoagland营养液中；低温处理将幼苗置于4℃的光照培养箱中；高温处理将幼苗置于40℃的光照培养箱中。在处理后的0h、1h、3h、6h、12h、24h分别采集叶片样本，按照上述方法保存，用于分析基因在不同逆境胁迫下的表达模式变化。2.1.2菌株和载体实验中使用的大肠杆菌菌株为DH5α和BL21(DE3)。DH5α菌株常用于质粒克隆，其基因型为F⁻，φ80dlacZΔM15，Δ(lacZYA-argF)U169，deoR，recA1，endA1，hsdR17(rk⁻，mk⁺)，phoA，supE44，λ⁻，thi-1，gyrA96，relA1。该菌株具有易于转化、生长速度快等特点，在质粒的扩增和保存中发挥重要作用。BL21(DE3)菌株则用于蛋白表达，其基因型为F⁻，ompT，hsdS（rBB-mB⁻），gal，dcm（DE3）。该菌株含有噬菌体T7RNA聚合酶基因，整合于染色体上的λ噬菌体DE3区，能够高效表达克隆于含有噬菌体T7启动子的表达载体的基因，适合用于外源蛋白的大量表达。这两种菌株均购自北京全式金生物技术有限公司，收到菌株后，将其接种于含有相应抗生素的LB培养基中，37℃振荡培养过夜，然后保存于-80℃冰箱中备用。实验所用的表达载体为pET-28a(+)，该载体购自Novagen公司。pET-28a(+)载体含有T7启动子、His标签序列、多克隆位点等元件。T7启动子能够被T7RNA聚合酶特异性识别并启动转录，从而高效表达目的基因；His标签序列便于后续利用镍柱亲和层析法对表达的重组蛋白进行纯化；多克隆位点包含多个限制性内切酶识别位点，方便目的基因的插入。在构建重组表达载体时，利用限制性内切酶对pET-28a(+)载体和目的基因进行双酶切，然后通过T4DNA连接酶将两者连接起来，转化至大肠杆菌DH5α感受态细胞中进行扩增，经测序验证正确后，再将重组质粒转化至BL21(DE3)感受态细胞中进行蛋白表达。2.2实验方法2.2.1基因克隆参考NCBI数据库中已公布的陆地棉基因组序列，利用PrimerPremier5.0软件设计特异性引物，用于扩增GhCNGC13和GhCNGC32基因。引物设计时，确保引物长度在18-25bp之间，GC含量在40%-60%，Tm值在55-65℃左右，同时避免引物二聚体和发夹结构的形成。引物序列由生工生物工程（上海）股份有限公司合成。以陆地棉‘中棉所49’的基因组DNA为模板进行PCR扩增。PCR反应体系为25μL，包括10×PCRBuffer2.5μL，dNTPMix（2.5mMeach）2μL，上下游引物（10μM）各0.5μL，模板DNA1μL，TaqDNA聚合酶（5U/μL）0.2μL，ddH₂O18.3μL。PCR反应程序为：95℃预变性3min；94℃变性30s，55-60℃退火30s（根据引物Tm值调整退火温度），72℃延伸1-2min（根据基因片段长度调整延伸时间），共35个循环；72℃延伸10min。PCR扩增产物经1%琼脂糖凝胶电泳检测，使用凝胶成像系统观察并拍照记录结果。在紫外灯下，用干净的手术刀切下含有目的条带的凝胶块，放入1.5mL离心管中，利用凝胶回收试剂盒按照说明书操作步骤回收目的DNA片段。将回收的目的DNA片段与克隆载体pMD18-T连接。连接体系为10μL，包括pMD18-T载体1μL，目的DNA片段4-5μL，SolutionI5μL。轻轻混匀后，16℃连接过夜。连接产物转化大肠杆菌DH5α感受态细胞。在超净工作台中，将10μL连接产物加入到100μLDH5α感受态细胞中，轻轻混匀，冰浴30min；然后42℃热激90s，迅速放回冰浴中2-3min；加入800μL不含抗生素的LB液体培养基，37℃振荡培养1h，使菌体复苏。取200μL转化后的菌液均匀涂布在含有氨苄青霉素（Amp）、IPTG和X-Gal的LB固体培养基平板上，37℃倒置培养12-16h。次日，观察平板上菌落生长情况，挑取白色单菌落接种于含有Amp的LB液体培养基中，37℃振荡培养过夜。提取质粒DNA，进行PCR鉴定和双酶切鉴定，将鉴定正确的阳性克隆送测序公司进行测序，以确保克隆的基因序列准确性。2.2.2生物信息学分析利用在线工具ExPASy（/）中的ProtParam程序分析GhCNGC13和GhCNGC32基因编码蛋白的氨基酸组成、分子量、理论等电点、不稳定系数、脂肪族氨基酸指数和总平均亲水性等基本理化性质。通过SOPMA（https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=/NPSA/npsa_sopma.html）预测蛋白的二级结构，包括α-螺旋、β-折叠、β-转角和无规则卷曲等结构元件的比例和分布。运用SWISS-MODEL（/）在线服务器构建蛋白的三级结构模型，并使用PyMOL软件对模型进行可视化分析和结构优化。借助InterProScan（https://www.ebi.ac.uk/interpro/search/sequence-search）工具预测蛋白的功能结构域，确定其是否含有环核苷酸结合结构域、离子通道结构域等与CNGC蛋白功能相关的保守结构域。通过BLASTp程序在NCBI数据库中进行同源序列搜索，获取不同物种中与GhCNGC13和GhCNGC32蛋白同源性较高的序列，利用MEGA7.0软件采用邻接法（Neighbor-Joiningmethod）构建系统进化树，分析它们之间的进化关系，明确这两个基因在CNGC基因家族中的分类地位。2.2.3基因表达分析采用Trizol法提取陆地棉‘中棉所49’不同组织（根、茎、叶、花、蕾、纤维等）以及不同逆境胁迫（干旱、盐胁迫、低温、高温）处理下叶片的总RNA。具体操作步骤如下：取适量植物组织，放入液氮预冷的研钵中，迅速研磨成粉末状；加入1mLTrizol试剂，继续研磨至组织粉末完全溶解；室温放置5min，使细胞充分裂解；加入200μL氯仿，剧烈振荡15s，室温放置3min；4℃，12000g离心15min，此时混合物分为三层，小心吸取上层无色水相转移至新的离心管中；加入等体积的异丙醇，轻轻混匀，室温放置10min；4℃，12000g离心10min，弃上清；加入1mL75%乙醇洗涤沉淀，涡旋振荡，悬浮沉淀；4℃，12000g离心5min，弃上清；室温晾干沉淀5-10min（注意不要使RNA沉淀过于干燥，以免影响溶解）；加入适量RNase-free水溶解RNA沉淀。利用超微量分光光度计检测RNA的浓度和纯度，要求OD₂₆₀/OD₂₈₀比值在1.8-2.0之间，表明RNA纯度较高，无蛋白质和基因组DNA污染。通过1%琼脂糖凝胶电泳检测RNA的完整性，观察28S和18SrRNA条带的亮度和清晰度，若28S条带亮度约为18S条带的2倍，且无明显拖尾现象，则说明RNA完整性良好，可以用于后续实验。按照逆转录试剂盒说明书的操作步骤，将提取的总RNA逆转录合成cDNA第一链。逆转录反应体系为20μL，包括5×PrimeScriptBuffer4μL，PrimeScriptRTEnzymeMixI1μL，Oligo(dT)₁₈Primer1μL，Random6mers1μL，TotalRNA1μg，RNase-freedH₂O补足至20μL。反应程序为：37℃15min，85℃5s，4℃保存。逆转录得到的cDNA保存于-20℃冰箱备用。以cDNA为模板，利用实时荧光定量PCR（qRT-PCR）技术检测GhCNGC13和GhCNGC32基因的表达水平。qRT-PCR反应体系为20μL，包括2×SYBRGreenMasterMix10μL，上下游引物（10μM）各0.8μL，cDNA模板1μL，ddH₂O7.4μL。反应在实时荧光定量PCR仪上进行，反应程序为：95℃预变性30s；95℃变性5s，60℃退火30s，共40个循环；熔解曲线分析从65℃到95℃，每升高0.5℃采集一次荧光信号。以陆地棉内参基因（如GhUBQ7）作为对照，采用2⁻ΔΔCt法计算目的基因的相对表达量。每个样品设置3个生物学重复和3个技术重复，以确保实验结果的准确性和可靠性。2.2.4基因功能验证构建GhCNGC13和GhCNGC32基因的过表达载体。以测序正确的克隆质粒为模板，使用带有酶切位点的特异性引物进行PCR扩增，获得目的基因片段。引物设计时，在引物的5'端分别添加合适的限制性内切酶识别位点（如BamHI和SacI等），以便后续与表达载体连接。PCR扩增产物经凝胶回收后，与经过同样限制性内切酶双酶切的表达载体pBI121（含有CaMV35S启动子和GUS报告基因）进行连接。连接体系和转化步骤同基因克隆部分所述。通过PCR鉴定、双酶切鉴定和测序验证，筛选出正确构建的过表达载体。利用农杆菌介导的遗传转化方法，将过表达载体导入陆地棉‘中棉所49’的胚性愈伤组织中。具体步骤如下：将含有过表达载体的农杆菌菌株（如EHA105）接种于含有相应抗生素的LB液体培养基中，28℃振荡培养至OD₆₀₀值为0.5-0.8；5000g离心5min，收集菌体，用含有100μM乙酰丁香酮（AS）的MS液体培养基重悬菌体，调整菌液浓度至OD₆₀₀值为0.3-0.5；将陆地棉胚性愈伤组织浸泡在农杆菌菌液中15-20min，期间轻轻振荡；取出愈伤组织，用无菌滤纸吸干表面多余菌液，转移至含有AS的共培养培养基（MS+0.1mg/L2,4-D+0.1mg/LKT+30g/L蔗糖+8g/L琼脂，pH5.8）上，25℃黑暗培养3d；共培养结束后，将愈伤组织转移至含有头孢霉素（Cef）和卡那霉素（Kan）的筛选培养基（MS+0.1mg/L2,4-D+0.1mg/LKT+30g/L蔗糖+8g/L琼脂+500mg/LCef+100mg/LKan，pH5.8）上，进行筛选培养，每2-3周更换一次筛选培养基，直至抗性愈伤组织出现；将抗性愈伤组织转移至分化培养基（MS+0.5mg/LIAA+1.0mg/LKT+30g/L蔗糖+8g/L琼脂+500mg/LCef+100mg/LKan，pH5.8）上，诱导分化出胚状体和再生植株；待再生植株长至3-4片真叶时，将其移栽至含有营养土的花盆中，置于温室中培养，进行后续表型鉴定。构建RNA干扰（RNAi）载体以抑制GhCNGC13和GhCNGC32基因的表达。根据目的基因的保守序列，设计长度为200-300bp的干扰片段，通过PCR扩增获得干扰片段。将干扰片段正向和反向插入到中间载体（如pFGC5941）的内含子两侧，形成反向重复结构，然后将含有反向重复结构的片段亚克隆到植物表达载体（如pCAMBIA1301）上，构建成RNAi载体。通过农杆菌介导的遗传转化方法，将RNAi载体导入陆地棉中，获得RNAi转基因植株。转化和筛选过程与过表达载体转化类似，但筛选培养基中的抗生素种类和浓度根据RNAi载体的抗性标记进行调整。对过表达和RNAi转基因植株进行表型鉴定。观察转基因植株在生长发育过程中的形态特征，如株高、叶片大小、分枝数、开花时间、棉铃大小和数量等，与野生型陆地棉进行比较，分析基因功能对植株生长发育的影响。同时，对转基因植株进行逆境胁迫处理（如干旱、盐胁迫、低温、高温等），测定相关生理生化指标，如相对电导率、丙二醛（MDA）含量、抗氧化酶活性（SOD、POD、CAT等）、渗透调节物质含量（可溶性糖、脯氨酸等）等，评估转基因植株对逆境胁迫的耐受性变化，从而验证GhCNGC13和GhCNGC32基因在陆地棉逆境响应中的功能。2.2.5亚细胞定位利用PCR技术扩增GhCNGC13和GhCNGC32基因的编码区序列（不含终止密码子），引物设计时在5'端分别添加合适的限制性内切酶识别位点（如KpnI和XbaI等），以便与绿色荧光蛋白（GFP）表达载体连接。PCR扩增产物经凝胶回收后，与经过同样限制性内切酶双酶切的GFP表达载体（如pCAMBIA1302-GFP）进行连接，构建融合表达载体。连接产物转化大肠杆菌DH5α感受态细胞，通过PCR鉴定、双酶切鉴定和测序验证，筛选出正确构建的融合表达载体。将构建好的融合表达载体转化农杆菌菌株（如GV3101），采用浸花法转化拟南芥或利用PEG介导的原生质体转化法转化烟草叶片原生质体。对于拟南芥转化，将含有融合表达载体的农杆菌菌液稀释至OD₆₀₀值为0.8-1.0，加入5%蔗糖和0.05%SilwetL-77，混匀后将拟南芥花序浸泡在菌液中30s左右，期间轻轻晃动；转化后的拟南芥植株置于黑暗条件下培养24h，然后转移至正常光照条件下培养，收获T₀代种子。T₀代种子经表面消毒后，播种在含有相应抗生素的MS固体培养基上筛选转基因植株，将筛选得到的阳性转基因植株移栽至花盆中培养，收获T₁代种子，继续进行筛选和扩繁，直至获得纯合的转基因株系。对于烟草叶片原生质体转化，取生长健壮的烟草叶片，用锋利的刀片切成0.5-1.0mm宽的叶条，放入含有酶解液（1.5%纤维素酶R-10+0.4%离析酶R-10+0.4M甘露醇+20mMKCl+10mMCaCl₂+10mMMES，pH5.7）的培养皿中，28℃黑暗条件下酶解3-4h；酶解结束后，用200目滤网过滤酶解液，收集原生质体，将原生质体悬浮在W5溶液（154mMNaCl+125mMCaCl₂+5mMKCl+2mMMES，pH5.7）中，冰浴30min；500g离心5min，弃上清，将原生质体悬浮在MMG溶液（0.4M甘露醇+15mMMgCl₂+4mMMES，pH5.7）中，调整原生质体浓度至1×10⁶个/mL；取100μL原生质体悬液，加入10μL融合表达载体（1-5μg/μL），轻轻混匀，再加入110μLPEG-CaCl₂溶液（40%PEG4000+0.2MCaCl₂+0.4M甘露醇），轻轻混匀，室温放置10-15min；加入400μLW5溶液终止转化反应，500g离心5min，弃上清，用W5溶液洗涤原生质体2-3次；将原生质体悬浮在含有10mMMES和0.5M甘露醇的培养液中，25℃黑暗条件下培养16-24h。利用激光共聚焦显微镜观察GFP荧光信号在细胞内的分布情况，确定GhCNGC13和GhCNGC32蛋白的亚细胞定位。在观察前，将转化后的拟南芥叶片或烟草叶片原生质体用蒸馏水清洗2-3次，去除表面杂质。设置合适的激光激发波长和发射波长，如GFP的激发波长为488nm，发射波长为505-530nm。观察并采集图像，分析荧光信号在细胞核、细胞质、细胞膜、细胞器等部位的分布情况，从而确定蛋白的亚细胞定位。2.2.6互作蛋白筛选采用酵母双杂交技术筛选与GhCNGC13和GhCNGC32相互作用的蛋白。以GhCNGC13和GhCNGC32基因的编码区序列为模板，分别构建酵母诱饵表达载体（如pGBKT7-GhCNGC13和pGBKT7-GhCNGC32），将诱饵载体转化酵母菌株（如Y2HGold），通过营养缺陷型培养基（SD/-Trp）筛选阳性转化子。提取阳性转化子的质粒，进行PCR鉴定和测序验证，确保诱饵载体构建正确。构建陆地棉cDNA文库，将文库质粒转化酵母菌株Y187。将含有诱饵载体的Y2HGold菌株和含有文库质粒的Y187菌株进行杂交，将杂交后的酵母细胞涂布在营养缺陷型培养基（SD/-Trp/-Leu/-His/-Ade+X-α-Gal+AureobasidinA）上，30℃培养3-5d，筛选生长且显蓝色的阳性克隆。对阳性克隆进行验证，提取阳性克隆中的文库质粒，转化大肠杆菌DH5α感受态细胞，通过PCR扩增和测序确定插入片段的序列。将测序结果在NCBI数据库中进行比对分析，筛选出与已知蛋白具有较高同源性的基因，进一步研究这些基因编码蛋白与GhCNGC13和GhCNGC32蛋白的相互作用关系。利用免疫共沉淀（Co-IP）技术验证酵母双杂交筛选得到的互作蛋白。将GhCNGC13和GhCNGC32基因与Flag标签融合，构建表达载体（如pBI121-Flag-GhCNGC13和pBI121-Flag-GhCNGC32），将互作蛋白基因与HA标签融合，构建表达载体（如pBI121-三、结果与分析3.1GhCNGC13和GhCNGC32基因克隆及序列分析通过设计特异性引物，以陆地棉‘中棉所49’的基因组DNA为模板进行PCR扩增，成功克隆得到了GhCNGC13和GhCNGC32基因。经1%琼脂糖凝胶电泳检测，在预期大小位置出现清晰的条带（图1），其中GhCNGC13基因条带大小约为[X]bp，GhCNGC32基因条带大小约为[X]bp，与理论值相符。将PCR扩增产物回收后连接到pMD18-T载体，转化大肠杆菌DH5α感受态细胞，挑取白色单菌落进行培养并提取质粒，经PCR鉴定和双酶切鉴定，均得到与预期一致的条带（图2），表明重组质粒构建成功。将鉴定正确的阳性克隆送测序公司进行测序，测序结果经DNAMAN软件分析，确认所克隆的基因序列准确无误。注：M为DNAMarker；1为GhCNGC13基因PCR扩增产物；2为GhCNGC32基因PCR扩增产物。注：M为DNAMarker；1-3为GhCNGC13重组质粒PCR鉴定；4-6为GhCNGC32重组质粒PCR鉴定；7-9为GhCNGC13重组质粒双酶切鉴定；10-12为GhCNGC32重组质粒双酶切鉴定。将获得的GhCNGC13和GhCNGC32基因序列在NCBI数据库中进行BLAST比对，分析它们与其他物种同源基因的序列相似性。结果显示，GhCNGC13基因与拟南芥AtCNGC13基因的序列相似性为[X]%，与水稻OsCNGC13基因的序列相似性为[X]%；GhCNGC32基因与拟南芥AtCNGC32基因的序列相似性为[X]%，与水稻OsCNGC32基因的序列相似性为[X]%。这表明GhCNGC13和GhCNGC32基因在不同物种间具有一定的保守性。进一步对GhCNGC13和GhCNGC32基因的结构进行分析，利用GSDS（GeneStructureDisplayServer）在线工具绘制基因结构示意图（图3）。结果表明，GhCNGC13基因含有[X]个外显子和[X]个内含子，外显子长度范围为[X]bp-[X]bp，内含子长度范围为[X]bp-[X]bp；GhCNGC32基因含有[X]个外显子和[X]个内含子，外显子长度范围为[X]bp-[X]bp，内含子长度范围为[X]bp-[X]bp。基因结构的差异可能导致其编码蛋白的功能有所不同，为后续深入研究基因功能奠定了基础。注：黄色框表示外显子，黑色线条表示内含子，数字表示外显子和内含子的长度（bp）。3.2生物信息学分析结果利用ExPASy中的ProtParam程序对GhCNGC13和GhCNGC32基因编码蛋白的理化性质进行分析，结果见表1。GhCNGC13蛋白由[X]个氨基酸组成，分子量约为[X]kDa，理论等电点（pI）为[X]，呈[酸性/碱性]。其不稳定系数为[X]，大于40，表明该蛋白不稳定。脂肪族氨基酸指数为[X]，总平均亲水性为[X]，显示该蛋白具有一定的亲水性。GhCNGC32蛋白由[X]个氨基酸组成，分子量约为[X]kDa，理论等电点为[X]，呈[酸性/碱性]。不稳定系数为[X]，小于40，说明该蛋白相对稳定。脂肪族氨基酸指数为[X]，总平均亲水性为[X]，同样具有亲水性。表1：GhCNGC13和GhCNGC32基因编码蛋白的理化性质蛋白名称氨基酸数目分子量(kDa)理论等电点不稳定系数脂肪族氨基酸指数总平均亲水性GhCNGC13[X][X][X][X][X][X]GhCNGC32[X][X][X][X][X][X]通过SOPMA对蛋白二级结构进行预测，结果显示，GhCNGC13蛋白中α-螺旋占[X]%，β-折叠占[X]%，β-转角占[X]%，无规则卷曲占[X]%（图4A）。α-螺旋结构在蛋白中形成稳定的螺旋状区域，可能参与维持蛋白的整体结构和功能；β-折叠结构则通过氢键相互作用形成片状结构，为蛋白提供一定的稳定性和刚性；β-转角和无规则卷曲结构较为灵活，可能参与蛋白与其他分子的相互作用以及信号传导过程。GhCNGC32蛋白中α-螺旋占[X]%，β-折叠占[X]%，β-转角占[X]%，无规则卷曲占[X]%（图4B）。虽然两个蛋白的二级结构元件比例存在差异，但都具有典型的α-螺旋和无规则卷曲结构，这与其他已知的CNGC蛋白结构特征相似。注：A为GhCNGC13蛋白二级结构；B为GhCNGC32蛋白二级结构。红色表示α-螺旋，黄色表示β-折叠，绿色表示β-转角，蓝色表示无规则卷曲。运用SWISS-MODEL在线服务器构建蛋白三级结构模型，并使用PyMOL软件进行可视化分析（图5）。结果表明，GhCNGC13和GhCNGC32蛋白均具有典型的六跨膜结构域（S1-S6），在第5和第6跨膜结构域之间存在一个孔道区域（P-loop），这是离子通道蛋白的典型结构特征，推测其可能参与离子的跨膜运输。在蛋白的C-末端和N-末端还存在一些其他结构域，如环核苷酸结合结构域（CNBD），该结构域含有高度保守的磷酸盐结合盒（PBC）和铰链区，能够与环核苷酸（如cAMP、cGMP）结合，从而调节离子通道的活性。此外，还发现可能存在钙调蛋白结合域（CaMBD），这表明GhCNGC13和GhCNGC32蛋白的活性可能受到钙离子/钙调蛋白（Ca²⁺/CaM）信号系统的调控。注：A为GhCNGC13蛋白三级结构；B为GhCNGC32蛋白三级结构。不同颜色的区域表示不同的结构域，其中跨膜结构域用紫色表示，孔道区域用蓝色表示，环核苷酸结合结构域用绿色表示，钙调蛋白结合域用黄色表示。通过InterProScan工具预测蛋白的功能结构域，进一步证实了GhCNGC13和GhCNGC32蛋白含有上述保守结构域，如离子通道结构域（PF00520）、环核苷酸结合结构域（PF00027）和钙调蛋白结合结构域（PF00612）等。这些结构域的存在与CNGC蛋白家族的功能密切相关，表明GhCNGC13和GhCNGC32可能具有类似CNGC蛋白的功能，即在植物生长发育和逆境响应过程中参与离子信号传导。为了探究GhCNGC13和GhCNGC32基因与其他物种CNGC基因的进化关系，利用MEGA7.0软件采用邻接法构建系统进化树（图6）。从NCBI数据库中获取了拟南芥、水稻、玉米等多个物种的CNGC蛋白序列，与GhCNGC13和GhCNGC32蛋白序列一起进行分析。结果显示，所有CNGC蛋白被分为4个主要的分支（GroupI-IV），其中GhCNGC13和GhCNGC32分别位于不同的分支。GhCNGC13与拟南芥AtCNGC13、水稻OsCNGC13等聚为一支，表明它们在进化上具有较近的亲缘关系，可能具有相似的功能；GhCNGC32与拟南芥AtCNGC32、水稻OsCNGC32等聚为另一支。同一分支内的CNGC蛋白在结构和功能上可能具有较高的保守性，而不同分支之间的CNGC蛋白可能在进化过程中发生了功能分化，以适应不同物种的生长发育和环境需求。注：进化树中不同颜色的分支表示不同的组，GhCNGC13和GhCNGC32分别用红色和蓝色字体标注。分支上的数字表示自展值（Bootstrapvalue），数值越高表示该分支的可靠性越强。3.3基因表达模式分析利用实时荧光定量PCR（qRT-PCR）技术，检测了GhCNGC13和GhCNGC32基因在陆地棉‘中棉所49’不同组织中的表达情况，结果如图7所示。在根、茎、叶、花、蕾和纤维等组织中，GhCNGC13和GhCNGC32基因均有表达，但表达水平存在明显差异。其中，GhCNGC13基因在根中的表达量最高，显著高于其他组织（P<0.05）；在茎和叶中的表达量次之；在花、蕾和纤维中的表达量相对较低。这表明GhCNGC13基因可能在陆地棉根系的生长发育和生理功能中发挥重要作用。GhCNGC32基因在叶中的表达量最高，显著高于根、茎、花、蕾和纤维等组织（P<0.05）；在根和茎中的表达量也相对较高；在花、蕾和纤维中的表达量较低。叶作为植物进行光合作用的主要器官，GhCNGC32基因在叶中的高表达，暗示其可能参与叶片的光合作用、气体交换等生理过程，或者在叶片应对外界环境变化中起作用。注：不同小写字母表示在P<0.05水平上差异显著。进一步分析了GhCNGC13和GhCNGC32基因在陆地棉不同发育时期纤维中的表达变化，结果见图8。在纤维发育的起始期（0DPA，dayspostanthesis），GhCNGC13基因的表达量较低；随着纤维的伸长，在5DPA和10DPA时表达量逐渐升高，至10DPA时达到峰值；随后在15DPA和20DPA时表达量又逐渐下降。这表明GhCNGC13基因可能在纤维伸长阶段发挥重要作用，参与调控纤维细胞的伸长过程。GhCNGC32基因在纤维发育的起始期（0DPA）表达量相对较高，之后在5DPA时略有下降，在10DPA和15DPA时表达量又有所升高，至15DPA时达到较高水平，随后在20DPA时表达量下降。该基因在纤维发育不同时期的表达变化趋势与GhCNGC13基因有所不同，推测其可能在纤维发育的不同阶段发挥不同的功能，或者与GhCNGC13基因协同作用，共同调控纤维的发育过程。注：不同小写字母表示在P<0.05水平上差异显著。为了探究GhCNGC13和GhCNGC32基因在陆地棉应对逆境胁迫中的作用，对五叶期棉花幼苗进行干旱、盐胁迫、低温和高温处理，检测基因在不同胁迫处理下的表达变化，结果如图9所示。在干旱胁迫下，GhCNGC13基因的表达量在处理1h后迅速上调，达到对照的[X]倍，随后在3h时略有下降，但仍显著高于对照（P<0.05）；在6h和12h时表达量继续下降，至24h时接近对照水平。这表明GhCNGC13基因能够快速响应干旱胁迫，可能通过调节离子平衡、水分吸收等生理过程，参与陆地棉对干旱胁迫的适应。GhCNGC32基因在干旱胁迫下的表达量在处理1h后也有所上调，但上调幅度较小；在3h时表达量达到峰值，为对照的[X]倍；随后在6h和12h时逐渐下降，至24h时接近对照水平。与GhCNGC13基因相比，GhCNGC32基因对干旱胁迫的响应相对较为迟缓，但在胁迫后期可能发挥重要作用。在盐胁迫下，GhCNGC13基因的表达量在处理1h后显著上调，达到对照的[X]倍，随后在3h和6h时持续升高，至6h时达到最大值，为对照的[X]倍；在12h和24h时表达量逐渐下降。说明GhCNGC13基因对盐胁迫响应迅速且强烈，可能在陆地棉抵御盐胁迫、维持离子稳态过程中发挥关键作用。GhCNGC32基因在盐胁迫下的表达量在处理1h后也明显上调，在3h时达到峰值，为对照的[X]倍；随后在6h、12h和24h时逐渐下降。尽管GhCNGC32基因对盐胁迫的响应时间与GhCNGC13基因相似，但表达变化幅度略有不同，两者可能通过不同的机制参与陆地棉对盐胁迫的响应。在低温胁迫下，GhCNGC13基因的表达量在处理1h后迅速升高，达到对照的[X]倍，随后在3h时略有下降，但仍显著高于对照（P<0.05）；在6h、12h和24h时表达量持续下降。表明GhCNGC13基因能够快速响应低温胁迫，可能在陆地棉调节低温下的生理代谢、维持细胞功能等方面发挥作用。GhCNGC32基因在低温胁迫下的表达量在处理1h后有所上调，在3h时达到峰值，为对照的[X]倍；随后在6h、12h和24h时逐渐下降。与GhCNGC13基因类似，GhCNGC32基因也能对低温胁迫做出响应，但其具体作用机制可能与GhCNGC13基因有所不同。在高温胁迫下，GhCNGC13基因的表达量在处理1h后显著上调，达到对照的[X]倍，随后在3h和6h时持续升高，至6h时达到最大值，为对照的[X]倍；在12h和24h时表达量逐渐下降。说明GhCNGC13基因对高温胁迫响应迅速，可能参与陆地棉应对高温胁迫的生理调节过程。GhCNGC32基因在高温胁迫下的表达量在处理1h后明显上调，在3h时达到峰值，为对照的[X]倍；随后在6h、12h和24h时逐渐下降。与其他胁迫处理类似，GhCNGC32基因在高温胁迫下也能做出响应，但其表达变化模式与GhCNGC13基因存在一定差异，暗示两者在高温胁迫响应中可能扮演不同的角色。注：不同小写字母表示在P<0.05水平上差异显著。3.4基因功能验证结果通过农杆菌介导的遗传转化方法，成功获得了GhCNGC13和GhCNGC32基因的过表达和RNAi转基因陆地棉植株。对转基因植株进行了PCR鉴定和qRT-PCR检测，结果表明，过表达转基因植株中GhCNGC13和GhCNGC32基因的表达量显著高于野生型植株，而RNAi转基因植株中基因表达量显著低于野生型植株（图10），说明转基因植株构建成功。注：WT为野生型植株；OE为过表达转基因植株；RNAi为RNA干扰转基因植株。不同小写字母表示在P<0.05水平上差异显著。在正常生长条件下，观察转基因植株与野生型植株的生长状况。结果显示，GhCNGC13过表达转基因植株在生长前期（播种后1-3周）生长速度较快，株高显著高于野生型植株（P<0.05），叶片数量和大小也有所增加；但在生长后期（播种后4-6周），过表达植株出现早衰现象，叶片发黄、枯萎，生长速度减缓，最终株高与野生型植株无显著差异。GhCNGC13RNAi转基因植株在整个生长周期中，生长速度均明显慢于野生型植株，株高较矮，叶片较小且数量较少（图11A）。GhCNGC32过表达转基因植株在整个生长周期中，生长状况与野生型植株相比无明显差异；而GhCNGC32RNAi转基因植株在生长前期生长较为缓慢，株高和叶片大小略小于野生型植株，但在生长后期逐渐恢复，最终与野生型植株生长状况相近（图11B）。对转基因植株的产量相关性状进行统计分析，结果如表2所示。GhCNGC13过表达转基因植株的棉铃数量显著低于野生型植株（P<0.05），单铃重也略有下降，导致籽棉产量和皮棉产量均显著降低；GhCNGC13RNAi转基因植株的棉铃数量和单铃重均显著低于野生型植株，籽棉产量和皮棉产量下降更为明显。GhCNGC32过表达转基因植株的棉铃数量、单铃重、籽棉产量和皮棉产量与野生型植株相比，均无显著差异；GhCNGC32RNAi转基因植株的棉铃数量略有减少，单铃重和籽棉产量、皮棉产量与野生型植株相比无显著差异。表2：转基因植株产量相关性状统计分析植株类型棉铃数量(个/株)单铃重(g)籽棉产量(g/株)皮棉产量(g/株)WT[X][X][X][X]GhCNGC13-OE[X]a[X]ab**[X]a**[X]aGhCNGC13-RNAi[X]b**[X]b**[X]b**[X]bGhCNGC32-OE[X]a**[X]a**[X]a**[X]aGhCNGC32-RNAi[X]a**[X]a**[X]a**[X]a注：同一列数据后不同小写字母表示在P<0.05水平上差异显著。在逆境胁迫耐受性方面，对转基因植株进行干旱、盐胁迫、低温和高温处理，测定相关生理生化指标。在干旱胁迫下，GhCNGC13过表达转基因植株的相对电导率和丙二醛（MDA）含量显著低于野生型植株（P<0.05），抗氧化酶（SOD、POD、CAT）活性和渗透调节物质（可溶性糖、脯氨酸）含量显著高于野生型植株，表明其细胞膜受损程度较轻，抗氧化能力和渗透调节能力较强，对干旱胁迫的耐受性增强；GhCNGC13RNAi转基因植株的相对电导率和MDA含量显著高于野生型植株，抗氧化酶活性和渗透调节物质含量显著低于野生型植株，对干旱胁迫更为敏感（图12A）。GhCNGC32过表达转基因植株在干旱胁迫下，相对电导率和MDA含量也低于野生型植株，抗氧化酶活性和渗透调节物质含量高于野生型植株，但差异不如GhCNGC13过表达植株显著；GhCNGC32RNAi转基因植株对干旱胁迫的响应与野生型植株相比无明显差异（图12B）。在盐胁迫下，GhCNGC13过表达转基因植株的相对电导率和MDA含量显著低于野生型植株，抗氧化酶活性和渗透调节物质含量显著高于野生型植株，表现出较强的耐盐性；GhCNGC13RNAi转基因植株的相对电导率和MDA含量显著高于野生型植株，抗氧化酶活性和渗透调节物质含量显著低于野生型植株，耐盐性明显下降（图12C）。GhCNGC32过表达转基因植株在盐胁迫下，相对电导率和MDA含量有所降低，抗氧化酶活性和渗透调节物质含量有所升高，但差异不显著；GhCNGC32RNAi转基因植株的相对电导率和MDA含量略高于野生型植株，抗氧化酶活性和渗透调节物质含量略低于野生型植株，耐盐性略有下降（图12D）。在低温胁迫下，GhCNGC13过表达转基因植株的相对电导率和MDA含量显著低于野生型植株，抗氧化酶活性和渗透调节物质含量显著高于野生型植株，对低温胁迫的耐受性增强；GhCNGC13RNAi转基因植株的相对电导率和MDA含量显著高于野生型植株，抗氧化酶活性和渗透调节物质含量显著低于野生型植株，对低温胁迫更为敏感（图12E）。GhCNGC32过表达转基因植株在低温胁迫下，相对电导率和MDA含量低于野生型植株，抗氧化酶活性和渗透调节物质含量高于野生型植株，但差异不显著；GhCNGC32RNAi转基因植株对低温胁迫的响应与野生型植株相比无明显差异（图12F）。在高温胁迫下，GhCNGC13过表达转基因植株的相对电导率和MDA含量显著低于野生型植株，抗氧化酶活性和渗透调节物质含量显著高于野生型植株，耐高温能力增强；GhCNGC13RNAi转基因植株的相对电导率和MDA含量显著高于野生型植株，抗氧化酶活性和渗透调节物质含量显著低于野生型植株，耐高温能力下降（图12G）。GhCNGC32过表达转基因植株在高温胁迫下，相对电导率和MDA含量有所降低，抗氧化酶活性和渗透调节物质含量有所升高，但差异不显著；GhCNGC32RNAi转基因植株的相对电导率和MDA含量略高于野生型植株，抗氧化酶活性和渗透调节物质含量略低于野生型植株，耐高温能力略有下降（图12H）。注：A、B为干旱胁迫；C、D为盐胁迫；E、F为低温胁迫；G、H为高温胁迫。WT为野生型植株；OE为过表达转基因植株；RNAi为RNA干扰转基因植株。不同小写字母表示在P<0.05水平上差异显著。综上所述，GhCNGC13基因在陆地棉生长发育过程中具有重要作用，其过表达可促进植株前期生长，但后期易早衰，且会降低产量；同时，GhCNGC13基因过表达能显著增强陆地棉对干旱、盐胁迫、低温和高温等逆境胁迫的耐受性，而基因沉默则使植株对逆境胁迫更为敏感。GhCNGC32基因对陆地棉生长发育和产量的影响相对较小，但在一定程度上也参与了陆地棉对逆境胁迫的响应，其过表达可增强植株的抗逆性，基因沉默后植株抗逆性略有下降。3.5亚细胞定位结果通过构建GhCNGC13和GhCNGC32基因与绿色荧光蛋白（GFP）的融合表达载体，转化拟南芥和烟草叶片原生质体，利用激光共聚焦显微镜观察GFP荧光信号在细胞内的分布情况，以确定这两个基因编码蛋白的亚细胞定位。在转化了35S::GhCNGC13-GFP融合表达载体的拟南芥叶片表皮细胞和烟草叶片原生质体中，均观察到强烈的绿色荧光信号集中分布在细胞膜上（图13A、B）。而转化了35S::GFP空载的拟南芥叶片表皮细胞和烟草叶片原生质体中，绿色荧光信号均匀分布在整个细胞，包括细胞核、细胞质和细胞膜（图13C、D）。这表明GhCNGC13蛋白定位于细胞膜，暗示其可能在细胞膜上参与离子的跨膜运输以及细胞内外的信号传导过程。在转化了35S::GhCNGC32-GFP融合表达载体的拟南芥叶片表皮细胞和烟草叶片原生质体中，绿色荧光信号主要集中在细胞膜和内质网上（图14A、B）。通过与内质网标记蛋白（如mCherry-HDEL）共定位分析，进一步证实了GhCNGC32蛋白在内质网上的定位（图14C、D）。内质网是细胞内重要的细胞器，参与蛋白质和脂质的合成、加工以及运输等过程。GhCNGC32蛋白在内质网的定位，提示其可能在内质网相关的生理过程中发挥作用，如参与内质网与细胞膜之间的离子转运平衡调节，或者在内质网应激响应中参与信号传导。注：A为转化35S::GhCNGC13-GFP的拟南芥叶片表皮细胞；B为转化35S::GhCNGC13-GFP的烟草叶片原生质体；C为转化35S::GFP空载的拟南芥叶片表皮细胞；D为转化35S::GFP空载的烟草叶片原生质体。标尺=20μm。注：A为转化35S::GhCNGC32-GFP的拟南芥叶片表皮细胞；B为转化35S::GhCNGC32-GFP的烟草叶片原生质体；C为转化35S::GhCNGC32-GFP和mCherry-HDEL的拟南芥叶片表皮细胞共定位；D为转化35S::GhCNGC32-GFP和mCherry-HDEL的烟草叶片原生质体共定位。标尺=20μm。综上所述，GhCNGC13蛋白主要定位于细胞膜，而GhCNGC32蛋白定位于细胞膜和内质网。蛋白的亚细胞定位与其功能密切相关，这些结果为进一步研究GhCNGC13和GhCNGC32基因在陆地棉生长发育和逆境响应中的功能机制提供了重要线索。3.6互作蛋白筛选与分析通过酵母双杂交技术，以GhCNGC13和GhCNGC32蛋白为诱饵，对陆地棉cDNA文库进行筛选，共获得了[X]个与GhCNGC13相互作用的阳性克隆，[X]个与GhCNGC32相互作用的阳性克隆。对这些阳性克隆进行测序和比对分析，鉴定出了多个与GhCNGC13和GhCNGC32相互作用的蛋白，其中包括一些已知功能的蛋白和部分功能未知的假定蛋白。在与GhCNGC13互作的蛋白中，鉴定出了GhCaM1（钙调蛋白1）、GhCBL1（类钙调磷酸酶B亚基1）、GhMAPK3（丝裂原活化蛋白激酶3）等。GhCaM1是一种重要的钙离子感受器，能够与钙离子结合并调节下游靶蛋白的活性。已有研究表明，CaM与CNGC蛋白的互作在植物离子信号传导和逆境响应中发挥关键作用，如拟南芥中AtCNGC1和AtCNGC2的活性受到CaM的调控，从而参与植物对病原体的防御反应。在本研究中，GhCNGC13与GhCaM1的互作，暗示了GhCNGC13可能通过与GhCaM1结合，在钙离子信号介导的陆地棉生长发育和逆境响应过程中发挥作用。GhCBL1是CBL-CIPK信号通路中的重要成员，主要通过与CIPK蛋白相互作用，调节离子通道活性和离子稳态。GhCNGC13与GhCBL1的互作，提示GhCNGC13可能参与了CBL-CIPK信号通路，在陆地棉应对逆境胁迫时，通过调节离子平衡来增强植物的抗逆性。例如，在拟南芥中，CBL-CIPK信号通路参与了植物对盐胁迫的响应，通过调节离子转运蛋白的活性，维持细胞内离子稳态，从而提高植物的耐盐性。GhMAPK3是MAPK级联反应途径中的关键激酶，参与植物对多种生物和非生物胁迫的信号转导过程。当植物受到逆境胁迫时，MAPK级联反应被激活，通过磷酸化下游靶蛋白，调节植物的生理生化反应，增强植物的抗逆性。GhCNGC13与GhMAPK3的互作，表明GhCNGC13可能通过激活MAPK级联反应，参与陆地棉对逆境胁迫的响应。如在烟草中，MAPK级联反应参与了植物对病毒侵染的防御反应，通过调节防御相关基因的表达，增强植物的抗病毒能力。在与GhCNGC32互作的蛋白中，鉴定出了GhTIP1;1（液泡膜内在蛋白1;1）、GhHSP70（热休克蛋白70）、GhWRKY22（WRKY转录因子22）等。GhTIP1;1是水通道蛋白家族的成员，主要定位于液泡膜，参与细胞内水分运输和渗透调节。植物在应对逆境胁迫时，需要调节细胞内的水分平衡，以维持细胞的正常生理功能。GhCNGC32与GhTIP1;1的互作，表明两者可能协同作用，在陆地棉应对干旱、盐胁迫等逆境时，通过调节细胞内水分运输和离子平衡，增强植物的抗逆性。例如，在拟南芥中，水通道蛋白参与了植物对干旱胁迫的响应，通过调节水分运输，维持细胞的膨压，从而提高植物的耐旱性。GhHSP70是一种分子伴侣，在植物细胞内发挥着重要的保护作用，能够帮助其他蛋白正确折叠、组装和转运，维持细胞内蛋白质稳态。当植物受到逆境胁迫时，如高温、干旱、盐胁迫等，细胞内蛋白质的结构和功能容易受到损伤，GhHSP70的表达会显著上调，通过与受损蛋白结合，促进其正确折叠和修复，从而保护细胞免受损伤。GhCNGC32与GhHSP70的互作，暗示GhHSP70可能参与了GhCNGC32蛋白的折叠、组装或稳定性维持过程，在陆地棉应对逆境胁迫时，共同发挥作用，增强植物的抗逆性。GhWRKY22是WRKY转录因子家族的成员，在植物生长发育和逆境响应中发挥重要的调控作用。WRKY转录因子通过与靶基因启动子区域的W-box元件结合，调节基因的表达，参与植物对生物和非生物胁迫的响应。GhCNGC32与GhWRKY22的互作，表明GhWRKY22可能通过调控与GhCNGC32相关基因的表达，参与陆地棉的生长发育和逆境响应过程。例如，在水稻中，WRKY转录因子参与了植物对稻瘟病菌的防御反应，通过调节防御相关基因的表达，增强植物的抗病性。为了进一步验证酵母双杂交筛选得到的互作蛋白，利用免疫共沉淀（Co-IP）技术进行验证。以陆地棉叶片总蛋白为材料，分别使用抗Flag抗体和抗HA抗体进行免疫共沉淀实验，结果表明，GhCNGC13与GhCaM1、GhCBL1、GhMAPK3，以及GhCNGC32与GhTIP1;1、GhHSP70、GhWRKY22在体内均能发生相互作用（图15），进一步证实了酵母双杂交筛选结果的可靠性。注：A为GhCNGC13与GhCaM1、GhCBL1、GhMAPK3的免疫共沉淀结果；B为GhCNGC32与GhTIP1;1、GhHSP70、GhWRKY22的免疫共沉淀结果。Input为总蛋白对照，IP为免疫沉淀后的蛋白样品，IB为免疫印迹检测。综上所述，通过酵母双杂交和免疫共沉淀技术，筛选和验证了多个与GhCNGC13和GhCNGC32相互作用的蛋白。这些互作蛋白涉及钙离子信号传导、离子稳态调节、MAPK级联反应、水分运输、蛋白质折叠和转录调控等多个生物学过程，表明GhCNGC13和GhCNGC32在陆地棉生长发育和逆境响应中可能通过与这些互作蛋白形成复杂的分子网络，协同调控植物的生理生化反应，从而发挥其生物学功能。四、讨论4.1GhCNGC13和GhCNGC32基因功能分析本研究通过对陆地棉GhCNGC13和GhCNGC32基因的克隆、生物信息学分析、表达模式分析、功能验证、亚细胞定位以及互作蛋白筛选等一系列实验，深入探究了这两个基因的功能。从基因结构和序列分析结果来看，GhCNGC13和GhCNGC32基因在不同物种间具有一定的保守性，且都含有典型的CNGC蛋白结构域，如六跨膜结构域、环核苷酸结合结构域和钙调蛋白结合结构域等。这些保守结构域的存在，为其发挥离子通道功能提供了结构基础。在拟南芥中，AtCNGC1和AtCNGC2基因也含有类似的结构域，并且通过这些结构域与其他蛋白相互作用，参与植物的免疫反应。这表明GhCNGC13和GhCNGC32基因可能具有与拟南芥CNGC基因相似的功能。基因表达模式分析结果显示，GhCNGC13和GhCNGC32基因在陆地棉不同组织和不同发育时期的表达存在差异，且对干旱、盐胁迫、低温和高温等逆境胁迫均有响应。GhCNGC13基因在根中的表达量最高，且在纤维伸长阶段表达量升高，这表明该基因可能在根系生长发育以及纤维伸长过程中发挥重要作用。同时，在逆境胁迫下，GhCNGC13基因能够快速响应，其表达量在短时间内显著上调，说明该基因参与了陆地棉对逆境胁迫的早期响应过程。GhCNGC32基因在叶中的表达量最高，在纤维发育不同时期的表达变化趋势与GhCNGC13基因有所不同，暗示其在叶片生理功能以及纤维发育过程中可能具有独特的作用。在逆境胁迫下，GhCNGC32基因也能做出响应，但其表达变化模式与GhCNGC13基因存在差异，这表明两者可能通过不同的机制参与陆地棉的逆境响应。通过基因功能验证实验，进一步证实了GhCNGC13和GhCNGC32基因在陆地棉生长发育和逆境响应中的重要功能。GhCNGC13基因过表达可促进植株前期生长，但后期易早衰，且会降低产量；同时，能显著增强陆地棉对干旱、盐胁迫、低温和高温等逆境胁迫的耐受性，而基因沉默则使植株对逆境胁迫更为敏感。这说明GhCNGC13基因在陆地棉生长发育和逆境响应过程中具有多效性，其功能的发挥受到严格的调控。在拟南芥中，也有研究报道CNGC基因的过表达或沉默会影响植株的生长发育和抗逆性。例如，AtCNGC10基因过表达可增强拟南芥对低钾胁迫的耐受性。GhCNGC32基因对陆地棉生长发育和产量的影响相对较小，但在一定程度上也参与了陆地棉对逆境胁迫的响应，其过表达可增强植株的抗逆性，基因沉默后植株抗逆性略有下降。这表明GhCNGC32基因在陆地棉逆境响应中可能起到辅助或协同作用。虽然其对生长发育和产量的影响不如GhCNGC13基因明显，但在陆地棉应对复杂环境胁迫时，仍具有不可忽视的作用。亚细胞定位结果表明，GhCNGC13蛋白主要定位于细胞膜，而GhCNGC32蛋白定位于细胞膜和内质网。蛋白的亚细胞定位与其功能密切相关，GhCNGC13蛋白在细胞膜上的定位，使其能够直接参与细胞膜上的离子跨膜运输以及细胞内外的信号传导过程。而GhCNGC32蛋白在内质网的定位，提示其可能在内质网相关的生理过程中发挥作用，如参与内质网与细胞膜之间的离子转运平衡调节，或者在内质网应激响应中参与信号传导。已有研究表明，内质网在植物应对逆境胁迫时，会发生一系列的生理变化，包括内质网应激反应等。GhCNGC32蛋白在内质网的定位，可能使其在这些过程中发挥重要的调节作用。通过酵母双杂交和免疫共沉淀技术，筛选和验证了多个与GhCNGC13和GhCNGC32相互作用的蛋白。这些互作蛋白涉及钙离子信号传导、离子稳态调节、MAPK级联反应、水分运输、蛋白质折叠和转录调控等多个生物学过程。这表明GhCNGC13和GhCNGC32在陆地棉生长发育和逆境响应中可能通过与这些互作蛋白形成复杂的分子网络，协同调控植物的生理生化反应，从而发挥其生物学功能。例如，GhCNGC13与GhCaM1、GhCBL1和GhMAPK3的互作，表明其可能通过钙离子信号传导和MAPK级联反应途径，参与陆地棉对逆境胁迫的响应。GhCNGC32与GhTIP1;1、GhHSP70和GhWRKY22的互作，暗示其可能在水分运输、蛋白质折叠和转录调控等方面发挥作用，进而影响陆地棉的生长发育和逆境响应。在拟南芥中，也发现了CNGC蛋白与其他蛋白的相互作用，如AtCNGC1与AtCaM1的互作，参与了植物对病原体的防御反应。这些研究结果为深入理解GhCNGC13和GhCNGC32基因的功能机制提供了重要线索。4.2基因与棉花重要性状的关系在棉花的生长发育过程中，产量、品质和抗逆性是其重要的性状，直接影响着棉花的经济价值和种植效益。本研究通过对陆地棉GhCNGC13和GhCNGC32基因的功能鉴定，发现这两个基因与棉花的这些重要性状密切相关。从产量性状来看，GhCNGC13基因对陆地棉产量有着显著的影响。过表达GhCNGC13基因虽然在生长前期促进了植株的生长，但后期出现早衰现象，导致棉铃数量显著减少，单铃重略有下降，最终籽棉产量和皮棉产量均显著降低。这表明GhCNGC13基因在棉花生长发育过程中的表达调控对于维持正常的产量形成至关重要。基因的异常表达可能打破了植物体内的激素平衡、营养分配等生理过程，从而影响了棉铃的发育和形成。在实际生产中，如何合理调控GhCNGC13基因的表达，避免其对产量的负面影响，是需要进一步研究的问题。相比之下，GhCNGC32基因对陆地棉产量的影响相对较小。过表达和RNAi转基因植株的棉铃数量、单铃重、籽棉产量和皮棉产量与野生型植株相比，均无显著差异。这说明GhCNGC32基因在陆地棉产量形成过程中可能不是关键的调控基因，但并不排除其在其他方面对棉花生长发育的潜在作用。虽然它对产量的直接影响不明显，但可能通过参与其他生理过程，间接影响棉花的整体生长状况。在棉花纤维品质方面，本研究通过对不同发育时期纤维中基因表达的分析，发现GhCNGC13和GhCNGC32基因在纤维发育过程中均有表达，且表达模式存在差异。GhCNGC13基因在纤维伸长阶段表达量升高，表明其可能参与了纤维伸长的调控过程。纤维长度是棉花品质的重要指标之一，因此GhCNGC13基因可能通过影响纤维伸长，对棉花纤维品质产生影响。在纤维伸长过程中，GhCNGC13基因可能通过调节离子通道活性，影响细胞内离子浓度和信号传导，进而调控纤维细胞的伸长。GhCNGC32基因在纤维发育不同时期的表达变化趋势与GhCNGC13基因有所不同，其在纤维发育过程中的具体功能还需要进一步深入研究。虽然目前尚不清楚其对纤维品质的具体影响机制，但不同的表达模式暗示了它可能在纤维发育的其他阶段或其他生理过程中发挥作用。它可能与GhCNGC13基因协同作用，共同调控纤维的发育，也可能通过独立的途径影响纤维品质。抗逆性是棉花适应不良环境、保证产量和品质的重要特性。本研究结果表明，GhCNGC13和GhCNGC32基因均参与了陆地棉对干旱、盐胁迫、低温和高温等逆境胁迫的响应过程。GhCNGC13基因过表达能显著增强陆地棉对逆境胁迫的耐受性，而基因沉默则使植株对逆境胁迫更为敏感。在干旱胁迫下，过表达植株的抗氧化酶活性和渗透调节物质含量显著升高，细胞膜受损程度减轻，从而提高了植株的抗旱能力。这表明GhCNGC13基因在陆地棉抵御逆境胁迫、维持细胞稳态方面发挥着关键作用。GhCNGC32基因在一定程度上也参与了陆地棉的逆境响应，其过表达可增强植株的抗逆性，基因沉默后植株抗逆性略有下降。尽管其对逆境胁迫的响应程度不如GhCNGC13基因明显，但在陆地棉应对复杂环境胁迫时，仍然具有重要的作用。这两个基因可能通过不同的信号传导途径和生理机制，协同调控陆地棉的抗逆性。它们可能参与了植物体内的抗氧化系统、渗透调节系统以及激素信号传导等过程，