棉籽大小数量性状核苷酸定位及候选基因初步筛选

棉籽是陆地棉(Gossypium hirsutum L.)生物产量的重要决定因素之一,棉籽大小影响种子的生长发育,最终影响籽棉产量。为发掘控制棉籽大小的遗传位点和候选基因,以中国农业科学院棉花研究所培育的陆地棉ZR014121和EZ60为亲本构建了200个重组自交系(recombinant inbred line,RIL)的群体,在2年2个地点(河南省安阳县、河北省威县)对棉籽粒长、粒宽、籽粒面积、籽粒周长、籽粒圆度和籽指这6个棉籽大小相关性状进行表型鉴定,结合液相芯片基因分型RIL群体得到的高质量单核苷酸多态性(single nucleotide polymorphism,SNP)数据进行全基因组关联分析(genome-wide association study,GWAS),挖掘控制相关性状的数量性状核苷酸(quantitative trait nucleotide,QTN)位点并筛选候选基因。结果发现:在4个环境下,6个棉籽大小性状的表型数据均符合正态分布,除籽粒圆度外,其他性状均呈显著正相关,且在RIL群体中存在丰富的表型变异,变异系数范围为3.46%~10.82%。GWAS共检测到47个与6个棉籽性状关联的SNP位点,分布于21条染色体上,其中有8个SNP位点在2个或2个以上环境中能被稳定检测。通过分析置信区间内基因的功能和转录组表达数据,初步推测Gh_D05G1018和Gh_A10G1731为控制棉籽生长发育的候选基因。这些研究结果为棉籽大小相关性状的遗传改良奠定基础。

红蓝光促进陆地棉愈伤组织诱导和增殖

【目的】陆地棉(Gossypium hirsutum)的遗传转化技术存在转化周期长、转化效率低等问题,其中,愈伤组织增殖速度慢是导致棉花转化周期长的关键因素。探究陆地棉愈伤组织增殖的最适光照条件,为缩短棉花转化周期提供技术基础。【方法】以陆地棉受体品种WC的下胚轴为外植体诱导愈伤组织,设置红光、蓝光、红蓝光(1﹕1)和白光(CK)4种不同光照处理,研究不同光质下愈伤组织的增殖速度、形态特征,以及相关基因表达水平和相关酶活性的差异,解析光照对愈伤组织增殖的影响并获得最优光照条件。【结果】不同光质对诱导棉花外植体出愈和愈伤组织生长具有明显影响。在4种不同光照处理下,红蓝光相较其他光质具有更显著地促进外植体出愈和愈伤组织增殖的效果,表现为愈伤组织分化时间更短,7—15 d,即与其他处理存在明显表型差异,7 d愈伤组织鲜重红蓝光(0.39 g)>蓝光(0.34 g)>白光(0.24 g)>红光(0.23 g),15 d时,红蓝光(1.15 g)>蓝光(0.98 g)>白光(0.69 g)>红光(0.51 g),红蓝光照射下的愈伤组织鲜重是白光对照的1.65倍,比单一蓝光或红光照射下的愈伤组织鲜重增加16.5%或125.5%;与此相一致的是,红蓝光处理下,第二周的愈伤增殖率也最高,达14.67%,比红光处理(7.17%)增加1倍多。4种光质下,愈伤增殖及体胚再生相关基因的表达量,红蓝光处理均显示最高。此外,红蓝光处理下愈伤组织的过氧化氢酶(CAT)活性显著提高,活性氧(ROS)含量低于白光对照。【结论】不同光质条件下,棉花外植体的愈伤组织增殖速度不同,红蓝光(1﹕1)促进外植体出愈和愈伤增殖表现最优,其次为蓝光,再次为白光,单独的红光不利于棉花出愈和愈伤增殖。红蓝光(1﹕1)处理诱导光受体及其互作因子、愈伤增殖等相关基因的表达,增强愈伤组织内过氧化氢酶的活性,降低活性氧含量,促进细胞分裂,进而促进愈伤增殖。

棉花抗旱性综合评价及其相关性分析

为综合评价棉花种质资源的抗旱性,筛选优异抗旱品种(品系),并建立成株期产量抗旱预测模型。对50份棉花种质资源在苗期进行干旱处理,调查其表型性状并测定其生理指标,根据各性状/指标差异的显著性筛选抗旱指标,并计算抗旱系数和抗旱指数,利用主成分分析法和隶属函数法计算综合抗旱系数、综合抗旱指数、综合抗旱隶属函数值、加权抗旱性系数和抗旱性综合度量值,对棉花抗旱性进行综合评价和分级。进一步将苗期评价分级结果与成株期干旱胁迫下产量和纤维品质性状进行比较,以分析苗期抗旱鉴定模型对成株期预测的精度。苗期生理性状在提高抗旱鉴定模型精度上具有重要作用,基于苗期生理指标建立的抗旱鉴定模型能够很好地预测成株期的抗旱性。

基于重测序鉴定GbTMEM214基因在陆地棉基因组的插入位点

【目的】利用农杆菌介导法获得了转GbTMEM214基因陆地棉品系,旨在明确转基因棉花株系中T-DNA插入位点的序列特征及检测方法,为其生物安全评价提供依据。【方法】基于基因组重测序技术,利用BLASTn将测序数据与陆地棉TM-1基因组序列进行比对分析,设计特异性引物通过聚合酶链式反应(polymerase chain reaction,PCR)验证插入位点。【结果】携带目标基因GbTMEM214的T-DNA整合在陆地棉基因组D13号染色体57019068~57019106 bp,T-DNA插入导致受体基因组37 bp的片段缺失;通过PCR扩增获得携带GbTMEM214基因的T-DNA插入位点的侧翼序列,由此建立了转GbTMEM214基因棉花的特异性检测方法。【结论】基于基因组重测序技术获得了转GbTMEM214基因棉花的T-DNA插入位点及侧翼序列信息,可为转基因棉花的生物安全评价提供技术参考。

棉花独脚金内酯D27基因表达与功能分析

为探究D27基因在棉花纤维及植株发育过程中的生物学功能,本研究通过生物信息学方法鉴定得到棉属中18个D27基因成员,分析其理化性质、结构及进化关系,利用实时荧光定量PCR(qRT-PCR)技术探究了海岛棉纤维中6个GbD27基因的表达特异性,通过同源重组的方法获得了6个陆地棉GhD27基因,构建GhD27基因的病毒诱导的基因沉默(VIGS)载体并转化至棉花,使用qRT-PCR技术检测其沉默效率及独脚金内酯(SLs)处理下的表达模式,并借助酶联免疫吸附(ELISA)技术测定了沉默植株的内源独脚金内酯含量。结果显示,棉花D27基因家族的开放阅读框(ORF)长度为726~816 bp,编码的D27蛋白含241~271个氨基酸;D27蛋白主要由α-螺旋与无规卷曲构成;进化分析结果显示D27基因高度同源。海岛棉材料新海21的GbD27基因在0~35 d的棉纤维中均存在表达,可能对棉纤维发育存在重要调节作用。陆地棉材料新农大棉4号的GhD27基因在叶片和根部均得到有效沉默,表明已获得GhD27基因沉默植株。在陆地棉GhD27基因的沉默植株中,GhD27基因对独脚金内酯处理均有不同程度的响应,内源性独脚金内酯含量与GhD27基因的表达水平有一定关联。本研究为进一步了解棉花D27基因的分子生物学功能奠定了基础。

海岛棉种质资源遗传多样性及对黄萎病的抗性鉴定

【目的】通过田间抗病性与室内苗期抗病性联合鉴定和分析,掌握海岛棉种质资源的抗病性分布状况,为筛选出黄萎病的高抗资源提供基础。【方法】以本课题组的238份海岛棉种质资源为鉴定对象,于2023年在库尔勒和阿瓦提通过自然病田小区试验,并联合室内苗期接菌的方法,对海岛棉资源进行黄萎病抗病性鉴定和评价。【结果】对238份海岛棉种质资源的15个指标遗传多样性分析,变异系数范围为2.00%~36.97%,均值为16.05%。于不同环境进行抗性鉴定发现变异系数范围为8.40%~71.06%。结合4个不同环境分析未鉴定到对黄萎病完全免疫及高抗的种质资源材料,鉴定出30份抗病材料,194份耐病材料,14份感病材料。【结论】海岛棉种质资源对大部分黄萎病都具有抗耐性,可从海岛棉种质资源中鉴定和选育对黄萎病抗性强的种质资源。

90份转BT基因抗虫棉品种(系)在新疆早熟棉区的适应性评价

【目的】评价转BT基因抗虫棉品种(系)的适应性与应用价值,为棉花新品种的选育提供适合的亲本。【方法】以90份转BT基因抗虫棉为材料,研究其农艺性状和品质性状的变异情况及相关性,并以主要农艺性状和品质性状对90份转BT基因抗虫棉进行主成分和聚类分析。【结果】供试90份转BT基因抗虫棉品种(系)纤维品质、农艺性状间差异显著,其中44份转BT基因抗虫棉为早熟材料,46份转BT基因抗虫棉为中早熟材料;前5个主成分特征值均大于1,累计贡献率达71.62%,第1主成分与产量有关,第2主成分与纤维品质有关,第3主成分与植株性状有关。90份转BT基因抗虫棉划分为4个类群,第Ⅰ类群籽指性状表现好;第Ⅱ类群断裂比强度、伸长率较好;第Ⅲ类群单株铃数、铃重、衣分较高;第Ⅳ类群纤维品质较好。【结论】筛选出了G21-2、G21-14、G21-1、G21-77等综合性状优异、铃重均在5.00 g以上、衣分在42%以上的转BT基因抗虫棉品种(系)18份。

283份陆地棉机采相关农艺性状的遗传多样性分析

为了探究283份陆地棉种质资源株型和产量相关性状的遗传多样性,筛选陆地棉机采株型优异品种,对283份来自国内外的陆地棉种质资源的11个农艺性状进行相关性分析、主成分分析和聚类分析,评价各性状的遗传变异情况,利用转录组数据和KEGG通路富集对果枝长度QTL定位区间基因进行表达量分析和基因筛选。结果表明,在陆地棉11个农艺性状中,遗传多样性指数以始节高最高,为2.01;变异系数以下部果枝长度最高,且各性状之间均存在遗传多样性。聚类分析表明,283份陆地棉资源分为五类,基于棉花株型、棉花产量等表型筛选出了高产材料1份,短果枝材料14份,始节高大于25 cm的材料4份,紧凑株型材料5份,同时在果枝长度QTL区间中筛选出2个可能影响果枝长度的候选基因。

磷胁迫对不同品种棉花苗期根系的影响

为了探究低磷胁迫下不同品种的棉花对磷的反应机理,本试验采用水培方法,对21个不同基因型的棉花品种进行高磷(KH_(2)PO_(4)1×10^(-3)mol/L)处理和低磷(KH_(2)PO_(4)1×10^(-5)mol/L)处理,通过主成分分析等统计方法,探讨了不同磷含量对棉花幼苗期根系形态的影响。结果表明,低磷处理后的棉花的根体积、根表面积均明显增加;随着磷含量的降低,棉花根系生长也会呈下降趋势。筛选得到1个磷高效棉花品种、2个磷中等棉花品种和18个磷低效棉花品种。棉苗总根长度、根面积、根表面积、根体积等指标在不同磷处理和各品种之间有明显差别,可以用来评价棉花品种的耐磷能力。

4个异源四倍体棉种自然群体遗传多样性解析

利用4个四倍体棉种(陆地棉野生种系、海岛棉、毛棉、达尔文氏棉)共116份材料的9个表型性状数据进行遗传多样性和聚类分析发现:毛棉的遗传多样性高于其他3个四倍体棉种;气孔导度与蒸腾速率、细胞间隙CO_(2)浓度,叶蜜腺到叶柄距离与叶蜜腺数、叶绿素含量、叶面积,叶蜜腺数与叶绿素含量、叶面积呈极其显著正相关关系;主茎茸毛数与叶面积、叶蜜腺到叶柄距离、叶绿素含量、叶蜜腺数呈极其显著负相关关系;通过对4个四倍体棉种表型性状进行主成分分析和聚类分析,得出本研究中的陆地棉野生种系材料与海岛棉材料之间亲缘关系最近,但与毛棉材料亲缘关系最远。本研究为有效利用陆地棉野生种系和野生棉种质资源提供了理论依据。

不同水分胁迫处理对陆地棉核心种质资源的影响

为研究不同水分胁迫下陆地棉对干旱的响应能力并筛选棉花抗旱关键指标和优异抗旱种质资源,选取30份陆地棉核心种质为试验材料,以全生育期正常灌水处理为对照,设置播种后和花铃期各浇水1次和花铃期断水2次2种胁迫处理,在蕾期、花铃期和吐絮期测定相应指标,通过描述性统计分析、差异分析、主成分分析和相关性分析,采用抗旱综合度量值(D)进行各材料的抗旱性评价。结果表明,干旱胁迫对不同棉花材料的生长发育均有不同程度的影响,播种后和花铃期各浇一水处理的D值离散较大(0.285~0.774),能更清晰地区分不同品种的抗旱性。利用D值可将30份材料分为4类:第Ⅰ类为抗旱材料,包括‘中棉所41’‘新陆早7号’等6个品种;第Ⅱ类为中抗材料,包括‘晋棉46’‘新陆早31号’等11个品种;第Ⅲ类为敏感材料,包括‘中棉所17’‘鲁1138’等5个品种;第Ⅳ类为高敏材料,包括‘新陆中8号’‘酒棉8号’等8个品种。果枝始节数、果枝数、株高、皮面产量和衣分5个指标对干旱较为敏感,可作为棉花抗旱评价的关键指标。以上研究结果可为棉花抗旱育种提供参考。

陆地棉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启动子区包含了许多与植物激素反应、生长发育和逆境应答密切相关的顺式调控因子。通过对转录组数据分析以及荧光定量分析,初步验证该家族GhRTE6基因可能在棉花生长发育以及盐胁迫响应中发挥作用,可为阐明GhRTE基因耐盐的分子机制提供理论依据。

转基因棉花耐除草剂鉴定及遗传稳定性分析

以农杆菌介导喷花法将海岛棉中的纤维发育优异的目的基因GbHCT与耐草铵膦除草剂基因Bar导入到陆地棉中,经过体积分数梯度筛选,确定了非抗性植株在800μL/L体积分数下死亡,由于除草剂筛选可能会有假阳性结果,对除草剂抗性植株还需目的基因及标记基因的进一步鉴定。PCR鉴定结果显示,抗除草剂基因Bar与目的基因GbHCT连续两代稳定遗传,分析其农艺性状,2个转基因棉花的株高与对照相比均有所升高,有效铃数与对照组相比也有显著差异,初步判断这些农艺性状的改变与目的基因GbHCT的插入有关。试验为后续研究GbHCT基因在纤维发育的功能研究进一步提供了参考。

陆地棉TRM基因家族的鉴定及纤维品质相关优异单倍型分析

【目的】鉴定陆地棉TRM基因家族序列及理化性质,分析与棉纤维品质相关的优异单倍型基因差异。【方法】利用生物信息学分析陆地棉TRM基因家族进化关系、理化性质和聚类表达;利用基因单倍型效应分析筛选调控纤维品质性状(纤维长度、比强度、马克隆值)的候选基因。【结果】陆地棉TRM基因家族编码的氨基酸为376~1093,等电点为4.64~9.56。亚细胞定位预测发现多达60个陆地棉TRM基因家族定位于细胞核中。TRM基因家族含有较多的光响应、激素响应、胁迫响应和生长发育相关的元件。60个TRM基因家族在纤维发育时期优势表达,调控棉花纤维发育。每个基因的单倍型个数为1~8,并筛选到纤维长度、强度以及马克隆值相关的优异单倍型TRM基因分别有14、18和15个,其中有11个基因同时具有长度、强度、马克隆值3种改良纤维品质的优异单倍型。GH_D09G0775的Hap_4单倍型和GH_D03G1434的Hap_3单倍型在增加纤维长度、强度的同时降低了马克隆值。【结论】在陆地棉(Gossypium hirsutum L.)中鉴定出75个TRM基因家族成员分布在24条染色体上,系统进化将其分为ClusterⅠ~Ⅲ3个亚族。

基于无人机多光谱影像的棉花SPAD值及叶片含水量估测

利用多光谱快速准确获取棉花花铃期旱胁迫性状对棉花后期产量和纤维品质具有重要帮助。为实现棉花关键抗旱指标的快速检测,该研究以253份棉花品种为材料,设置干旱胁迫和正常灌溉2个处理,在花铃期通过大疆精灵4多光谱无人机获得图像,分析花铃期干旱胁迫和正常处理棉花的光谱反射率,结合地面调查叶绿素相对值(S_(PAD))和叶片含水量(L_(WC))数据,结合神经网络算法中的径向基函数进行模型的预测。结果显示:棉花干旱胁迫前后除红光反射值略上升,其余4个光谱均小于正常处理;通过径向基函数模型估测S_(PAD)和L_(WC)的最佳估测模型都为二次函数,决定系数R^(2)分别为0.8488和0.9366,均方根误差R_(MSE)分别为2.005和0.930,相对误差R_(E)分别为0.004和0.011;将2个模型应用于试验区影像,对S_(PAD)及L_(WC)预测值进行聚类分析,根据S_(PAD)聚类结果,试验棉花旱情等级划分为特旱、重旱、中旱、轻旱、无旱5类,根据LWC聚类结果,试验棉花旱情等级划分为特旱、重旱、中旱、干旱、轻旱、无旱6类。以上模型及分类效果良好且符合棉花干旱胁迫的特点,研究结果可为快速评价棉花抗旱性、挖掘棉花抗旱基因提供技术支持。

基于主成分分析和隶属函数的棉花种子萌发期耐涝性鉴定评价

参照棉花耐渍涝性鉴定技术标准,在棉花种子萌发期模拟淹水胁迫,采用主成分分析、隶属函数法及聚类分析方法,分析了涝渍胁迫处理对6份陆地棉材料出苗率、株高、主根长、子叶叶绿素含量(SPAD值)、幼苗鲜物质质量、幼苗干物质质量的影响。结果表明:幼苗鲜物质质量、株高、主根长、出苗率、幼苗干物质质量和SPAD值可以作为棉花萌发期耐涝性鉴定的指标;6份参试材料的耐涝性强弱顺序依次为402、晋棉36、冀丰908、宛棉3号、苏棉9号、TM-1,其中402和晋棉36为耐涝性较强的材料。相关研究方法可为批量鉴定和评价棉花种质资源的耐涝性提供依据。

235份陆地棉表型性状遗传多样性分析

本研究利用235份陆地棉种质资源材料,对棉花株型、生育期、产量性状展开研究。通过对16个表型性状进行遗传多样性分析、相关性分析、主成分分析、聚类分析及隶属函数计算得到综合评价值,对陆地棉种质资源进行综合评价。结果表明,235份陆地棉群体具有丰富的遗传多样性,性状遗传多样性指数为1.969~2.107。不同农艺性状指标通过主成分分析得到6个独立的因子,主效因子贡献率为6.20%~23.80%,累计贡献率为81.573%。依据特点分为产量相关因子、棉花生长结构因子、早熟相关性状因子等因子。聚类分析按照整体性状分为Ⅳ类,第Ⅰ类群占比15%,株高较高,株型松散;第Ⅱ类群占比27%,株型较为紧凑,生育期较短,产量相对较高,适合新疆地区的机械采收;第Ⅲ类群占比24%,综合表现最差;第Ⅳ类群占比34%,特点为生育期较长,产量一般。通过计算综合评价值并与表型性状进行相关性分析得出,综合评价值与株高、单铃质量、衣分、单铃皮棉质量和有效铃数等指标呈显著正相关。最后,利用线性回归方程逐步回归得到株高、有效铃数、中部果枝长度、落叶等级、单铃质量、有效果枝数、第一果枝高度、全生育期和中部果枝夹角等9个表型性状,这9个表型性状能够综合评价陆地棉的遗传多样性。本研究结果可为陆地棉种质资源的合理利用及亲本选育提供理论依据。

利用InDel标记划分玉米自交系的杂种优势群

利用北京市农林科学院玉米研究中心构建的多重扩增靶向捕获测序基因分型InDel标记组合,对102份玉米自交系进行InDel标记的基因分型。利用UPGMA方法将102份玉米自交系划分为6个杂种优势群包括瑞德群、塘四平头、兰卡斯特群、旅大红骨群、P群及热带和亚热带混血群。类群分析结果与现有的杂种群和已知谱系高度一致。

陆地棉种质资源材料的耐热性评价及指标筛选

【目的】评价陆地棉种质资源材料的耐热性及筛选指标,为选育耐热性棉花品种及研究耐热性机理提供参考依据。【方法】基于2年试验,选用24份棉花种质资源,随机区组设计,测定相关指标,综合评价24份陆地棉种质资源材料的耐热性。【结果】24份陆地棉资源材料可分为3类。其中第一类为耐热材料,第二类为中间型材料,第三类为敏热型材料,建立陆地棉耐热性的评价模型:D=-0.111+0.175X 2+0.470X 3+1.211X 4(R^(2)=0.9148),建立的最优回归方程可以预测陆地棉材料的耐热性。籽棉产量、单铃重和皮棉产量可作为棉花大田耐热的评价指标,而株高和衣分与耐热性无关。【结论】筛选出4份耐热陆地棉资源材料,建立了以籽棉产量、单铃重和皮棉产量作为陆地棉大田耐热性的评价指标,以过氧化物酶、丙二醛和叶绿素含量作为室内耐热性的评价指标。

棉花β-tubulin基因家族的鉴定及其在纤维发育中的表达

【目的】β-微管蛋白是棉纤维细胞形态建成的基本结构单位,在纤维发育过程中具有重要作用。通过鉴定棉花β-tubulin基因家族成员,并进行生物信息学和表达模式分析,为深入探究β-tubulin基因在棉花纤维发育中的作用提供理论依据。【方法】利用BLAST方法,在4个棉种基因组中鉴定β-tubulin基因家族成员,并结合ProtParam tool分析理化性质、MEGA7.0构建系统进化树、Mapchart2.2绘制染色体定位图、MEME分析保守基序、PlantCARE分析启动子顺式作用元件;根据39个材料发育纤维转录组数据分析陆地棉β-tubulin基因家族的表达水平,并利用Spearman相关性分析鉴定影响纤维品质性状形成的β-tubulin基因。【结果】在陆地棉(Gossypiumhirsutum,AD1)、海岛棉(Gossypiumbarbadense,AD2)、亚洲棉(Gossypiumarboretum,A2)和雷蒙德氏棉(Gossypiumraimondii,D5)基因组中分别鉴定到36、37、19和18个β-tubulin基因,且在四倍体棉种中的数目约是二倍体棉种数目的二倍,系统进化将其分为ClusterⅠ—ClusterⅤ共5个亚组。系统进化与共线性分析发现,与陆地棉β-tubulin基因家族相比,海岛棉与二倍体亚洲棉和雷蒙德氏棉亲缘关系更近。保守结构域均具有典型的Tubulin和Tubulin-C。理化性质分析表明,该家族基因编码的氨基酸数目为421—508,等电点为4.68—5.09。启动子顺式作用元件分析获得生长发育响应相关元件、激素响应相关元件和胁迫响应相关元件等,表明β-tubulin基因参与细胞的生长调节。利用TM-1转录组数据对36个陆地棉β-tubulin基因表达模式进行聚类分析,发现有42%基因在纤维中优势表达;此外,利用海陆群体动态纤维转录组数据筛选到1、6和11个β-tubulin基因分别与纤维马克隆值、比强度和纤维长度显著相关,其中4个基因同时影响纤维长度和比强度性状。【结论】4个棉种中共鉴定出110个β-tubulin基因家族成员,氨基酸理化性质和序列高度保守而启动子序列调控元件多样;筛选出陆地棉纤维优势表达的β-tubulin基因家族成员,并发掘潜在调控纤维发育的候选基因。