一个棉花液泡转化酶基因的克隆与功能分析
徐文亭, 王诚, 徐晓洋, 牛二利, 蔡彩平, 郭旺珍*
南京农业大学作物遗传与种质创新国家重点实验室 / 教育部杂交棉创制工程研究中心, 江苏南京210095
* 通讯作者(Corresponding author): 郭旺珍, E-mail:moelab@njau.edu.cn
摘要

棉纤维正常发育需要大量的蔗糖供应。转化酶是生物体内蔗糖代谢途径中的关键酶之一, 对转化酶基因的结构和功能研究将有助于揭示复杂的棉纤维发育分子机制, 并为纤维品质改良提供优良的基因资源。本研究以棉纤维突变体im和遗传标准系TM-1纤维发育中显著差异表达的EST序列(GenBank登录号为EY196825)为探针, 通过电子克隆方法, 结合cDNA及基因组全长基因PCR扩增、测序验证, 克隆到一个棉花液泡转化酶基因GhVacInv2a(GenBank登录号为KF305322)。该基因ORF全长1857 bp, 编码618个氨基酸, 与已登录的陆地棉GhVacInv2基因(GenBank登录号为FJ864677)序列一致性为99% (E=0)。基因组序列分析表明,GhVacInv2a在二倍体棉种非洲棉和雷蒙德氏棉中含1个拷贝, 在四倍体陆地棉和海岛棉中存在2个拷贝, 其中GhVacInv2aGhVacInv2分属于A、D亚组同源基因。该同源基因含7个外显子, 6个内含子。Q-PCR分析表明, 该基因在花药中表达量最高, 在快速伸长的纤维组织中优势表达。在13~19 DPA的纤维中, 其在TM-1中的表达水平极显著地高于im突变体。开发SNP标记, 将GhVacInv2a基因定位在异源四倍体棉花第3染色体上。标记与性状的关联分析显示,GhVacInv2a与纤维强度存在显著相关(P=0.0087), 表明GhVacInv2a在棉花纤维品质形成中可能发挥重要功能。

关键词: 棉花; 纤维发育; 液泡转化酶2; 克隆; 功能分析
Cloning and Functional Analysis ofGhVacInc2a Encoding Vacuolar Invertase in Cotton
XU Wen-Ting, WANG Cheng, XU Xiao-Yang, NIU Er-Li, CAI Cai-Ping, GUO Wang-Zhen*
State Key Laboratory of Crop Genetics and Germplasm Enhancement / Hybrid Cotton R&D Engineering Research Center, Ministry of Education, Nanjing Agricultural University, Nanjing 210095, China
Abstract

Cotton fiber development depends on a large supply of sugar. Invertase plays a central role in sucrose metabolism. Structural and functional analysis of gene encoding invertase will conductively reveal the complex molecular mechanism of fiber development, and also provide elite gene resources for improvement of cotton fiber quality. In the study, a novel gene encoding vacuolar invertase was obtained by sillico cloning, based on a differentially expressed EST (GenBank accession number EY196825) betweenG. hirsutum acc. TM-1 andim mutant, combined with reconfirmation of ORF information in transcriptional and genomic level. BlastN showed that the gene shared 99% identity ofGhVacInv2reported previously (GenBank accession number FJ864677), and designated asGhVacInv2a (GenBank accession number KF305322).GhVacInv2a contained an open reading frame of 1857 bp that encoded a polypeptide of 618 amino acids, with seven exons and six introns at genome sequence.GhVacInv2a comprised one copy in diploid cotton speciesG. herbaceum andG. raimondii, and two copies in tetraploid cotton speciesG. hirsutum acc. TM-1 andG. barbadense cv. Hai 7124, withGhVacInv2ain A-subgenome andGhVacInv2 in D-subgenome, respectively. Q-PCR expression analysis showed that the set of homoelogous genes performed the highest expressive abundance in anther in comparison with other tissues. In developing fibers, it was predominantly expressed from fiber rapid elongation stage, and there existed significant difference between TM-1 andim mutant during 13 to 19 days post anthesis (DPA) in fiber tissues. Further,GhVacInv2a was located on chromosome 3 by developing subgenome-specific SNP marker, and the association analysis showed that there was a significant correction (P=0.0087) betweenGhVacInv2a and fiber strength, suggesting thatGhVacInv2aprobably plays a key role in fiber qualities formation.

Keyword: Cotton; Fiber development; Vacuolar invertase 2; Cloning; Functional analysis

在大多数植物中, 光合作用固定的碳源常以蔗糖形式转运, 利用蔗糖从源叶中输送能量和碳到植物发育的各个组织中。转化酶在植物碳源利用上起着关键调节作用, 它将蔗糖分解为葡萄糖和果糖, 随后这2种己糖参与多个生物代谢途径[ 1]。根据最适pH值, 转化酶可被分为3种类型。中性和碱性转化酶的最适pH为7.0~7.8, 在多种组织中表达。酸性转化酶的最适pH值在4.5~5.0之间, 可进一步细分为细胞壁型(CIN)和液泡型(VIN)[ 1, 2]。葡萄[ 3]、马铃薯[ 4]、番茄[ 5]等物种中研究显示, 液泡转化酶(VIN)在决定细胞碳源分配中起着重要的调节作用。目前普遍认为转化酶通过调节细胞内的渗透势来参与细胞发育。因为VIN水解1分子蔗糖变成2分子的己糖, 增加了渗透势调节物质, 胞内膨压加大, 有利于驱动细胞扩张[ 6]。但也有报道认为, 在糖类物质并非是主要渗透调节物的组织细胞中, VIN可能是通过作用于细胞壁的延伸来促进细胞扩张[ 7]。其他研究也表明VIN还参与诸如植物生长素、逆境刺激等一系列应答反应[ 6, 8, 9]

棉纤维是从胚珠外表皮发育而来的单细胞无分支结构。纤维的发育可分为起始、伸长、次生壁加厚及脱水成熟4个相互重叠的阶段[ 10]。因在发育过程中表现出的优良特性, 棉纤维被视为研究细胞发育的理想模式系统[ 11]。目前, 棉花上一共报道了 GhVIN1[ 12] GhVacInv1 GhVacInv2[ 13] 3个VIN基因。其中, GhVIN1是与棉纤维伸长相关的VIN基因, 它的表达水平与纤维伸长速率相吻合。抑制和过表达 GhVIN1分别会减缓和促进纤维伸长。将 GhVIN1转入拟南芥中, 该基因能补偿拟南芥 VIN插入突变体的表型, 并能促进野生型拟南芥根伸长。而 GhVacInv1 GhVacInv2的不同组织器官表达特征表明, 这2个基因均在纤维伸长的细胞中优势表达, 可能参与纤维细胞的扩展、伸长。

本研究基于一个在棉纤维突变体 im和遗传标准系TM-1纤维发育中显著差异表达的EST序列(GenBank登录号为EY196825)为探针, 利用电子克隆, 结合cDNA及基因组全长基因PCR扩增、测序验证, 获得一个与棉花纤维发育相关的液泡转化酶基因 GhVacInv2a。进一步对其结构、功能及分子进化特征系统研究, 以期为深入阐明棉纤维发育的分子机制, 并用于纤维品质改良提供基因资源。

1 材料与方法
1.1 试验材料

2010—2012年于南京农业大学江浦试验站种植陆地棉遗传标准系TM-1 ( G. hirsutum acc. TM-1)和 im纤维不成熟突变体 ( G. hirsutum acc. im, 次生壁发育缺陷的纤维突变体) 材料。在植株生长盛花期, 分别采集2个材料开花后10、13、16、19、22、25 DPA (days post anthesis)的纤维, 于液氮中速冻, -70℃冰箱保存。根、茎和叶样品取自有五、六片真叶的棉花幼苗。其他二倍体A染色体组的阿非利加棉 ( G. herbaceumvar. africanum)、D染色体组的雷蒙德氏棉( G. raimondii)和异源四倍体栽培棉种海岛棉海7124 ( G. barbadense cv. Hai7124)均种植于南京农业大学牌楼实验基地。

1.2 基因克隆及氨基酸序列分析

基于遗传标准系TM-1和棉纤维突变体 im竞争杂交芯片结果, 获得一个在TM-1中优势表达的EST序列(GenBank登录号为EY196825)。以该EST序列为探针, 通过BlastN比对程序, 搜索陆地棉EST数据库(http://www.ncbi.nlm.nih.gov/)。将获得的EST序列用CAP3程序(http://pbil.univ-lyon1.fr/ cap3.php)拼接后形成Contig, 然后以Contig为探针, 再次使用BlastN比对程序检索陆地棉EST数据库, 直至检出的Contig不能再延伸。对电子拼接获得的序列进行ORF预测, 参考BlastX结果完成其命名。利用电子拼接的序列信息设计引物, 通过RT-PCR和PCR方法分别扩增转录组cDNA和基因组DNA, 进行目标基因序列验证及不同棉种同源基因克隆分析。表1为本研究中涉及的所有引物信息。由南京金斯瑞生物科技有限公司完成所有引物合成和DNA序列测定。

将获得的cDNA全长序列运用NCBI的CD搜索工具(http://www.ncbi.nlm.nih.gov/Structure/cdd/ wrpsb.cgi)进行编码产物保守结构域分析。利用ExPASy的Compute p I/Mw tool (http://us.expasy.org/ tools/pi_tool.html)预测编码氨基酸的等电点和分子量大小, 运用SignalP程序对由 GhVacInv2a基因编码的氨基酸序列进行N端信号肽结构的预测。

使用HMMER3.0程序[ 14]根据 GhVIN2基因的保守结构域搜索DOE Joint Genome Institute释放的棉花雷蒙德氏棉( Gossypium raimondii)蛋白数据库(http://www.phytozome.net/cotton.php), 获得雷蒙德氏棉D基因组中与 GhVacInv2a基因功能类似的家族

表1 本研究所用引物 Table 1 Primer pairs used in this study

基因。用ClustalX 1.8软件对本研究克隆的 GhVacInv2a基因与棉花、拟南芥、杨树、葡萄等物种的转化酶基因编码产物进行多重序列比对, 用MEGA5 (http://www.megasoftware.net/)绘制序列的进化树。

1.3 总DNA和RNA的提取

利用改进的热硼酸法提取13、16、19、22和25 DPA的纤维组织RNA[ 15], 采用CTAB酸酚法[ 16]提取5~6片真叶棉花的根、茎和叶总RNA; 采用CTAB法[ 17]从新鲜幼嫩的叶片提取DNA。

1.4 实时荧光定量PCR

利用实时荧光定量PCR分析 GhVacInv2a同源基因在TM-1和 im突变体不同组织器官及13~25 DPA纤维中的基因表达水平。共设3次生物学重复。以棉花 histon 3基因(GenBank登录号为AF024716)为内参照对表达数据进行标准化。使用罗氏公司Light Cycler Fast Start DNA Master SYBR Green I Kit (Roche), 应用生物公司AB7500 Real-time PCR检测系统(Applied Biosystems)进行PCR扩增并记录荧光信号值。定量PCR程序为95℃ 10 min; 95℃ 15 s, 60℃ 1 min, 40个循环; 以0.2℃ s-1的速度从65℃逐渐升至95℃来制作溶解曲线。为了确保定量PCR结果的可信度, 对 GhVacInv2a基因与内参照 histon 3基因均制作标准曲线来确定基因的扩增效率, 确保标准曲线 R2>0.98, 扩增效率在90%~105%之间。使用Pfaffl方法分析基因在不同材料和发育时期的相对表达水平[ 18]

其中, E: 扩增效率; target: 目标基因; reference: 内参基因; control: 对照; sample: 分析样品。

1.5 基因定位及标记性状关联分析

根据该基因在遗传标准系TM-1和海岛棉海7124中基因组序列的差异, 开发SNP标记, 利用本实验构建的包含138个单株的[(TM-1×海7124)× TM-1] BC1作图群体进行基因初步定位[ 19]。根据[(TM-1×海7124)×TM-1] BC1作图群体 GhVacInv2a基因的定位结果, 利用以TM-1为背景渐渗海7124染色体片段的单片段置换系CSIL028与 im突变体构建的F2作图群体[ 20], 对 GhVacInv2a基因进一步精细定位。并结合F2:3家系纤维品质及产量性状进行基因与性状的关联分析。

2 结果与分析
2.1 GhVacInv2a基因cDNA全长序列的获得及其分析

在陆地棉遗传标准系TM-1 ( G. hirsutum acc. TM-1)和不成熟纤维突变体( im)纤维发育竞争杂交芯片上, 选取一个TM-1材料优势表达的EST序列(GenBank登录号为EY196825), 以该EST序列为种子序列在NCBI棉花数据库中进行比对。通过电子克隆方法拼接获得一个最长为2113 bp的序列。根据电子延伸的序列设计PCR引物(表1), 以TM-1开花后10 d的纤维组织cDNA为模板, 扩增出长为2003 bp的cDNA序列, 测序验证和原电子拼接序列一致。生物信息学分析显示, 该序列与已报道的陆地棉液泡转化酶 GhVacInv2基因(GenBank登录号为FJ864677)在核苷酸和蛋白水平上均有99%的一致性。进一步序列比对发现该基因与蓖麻( Ricinus communis)的转化酶(β-果糖呋喃苷酶)、木薯( Manihot esculenta)的液泡转化酶分别具有74%和73%的较高同源性。因此, 参照其在棉花中同源基因的命名将该基因命名为 GhVacInv2a( Gossypium hirsutum VacInv2a; GenBank登录号为KF305322)。 GhVacInv2a基因ORF长为1857 bp, 编码618个氨基酸。用Expasy p I/Mw程序对该基因编码氨基酸一级结构分析发现, 该基因等电点为5.25, 分子量为69.2 kD。根据NCBI CDS (http://www.ncbi.nlm.nih.gov/Structure/cdd/)分析, 其编码的氨基酸残基含有GH32_B_ Fructosidase-motif, 在100~420氨基酸之间有属于GH43_ 62_32_68superfamily保守域。运用SignalP程序对该基因的氨基酸序列分析, 发现N端有信号肽序列, 预测在49~50位氨基酸之间裂解。推测该基因是一个蛋白前体, 它在被合成之后将运输到液泡中。

2.2 GhVacInv2a同源基因在不同棉种基因组序列分析

利用表1中设计2套巢式PCR引物, 分别扩增TM-1、海7124以及与其祖先种亲缘关系最近的二倍体棉种阿非利加棉和雷蒙德氏棉基因组DNA, 将获得的PCR产物回收、连接、转化, 随机挑选多个阳性克隆测序并序列比对分析, 获得该基因在二倍体和四倍体棉种中的基因序列。运用在线软件 GSDS (Gene Structure Display Server, http://gsds.cbi. pku.edu.cn/chinese.php)对基因结构进行分析, 该基因包含7个外显子, 6个内含子。系统进化分析表明(图1), 该基因在四倍体棉种中形成A亚组和D亚组两大分支, 可明确追溯到其祖先来源, 说明在四倍体棉种中来源于A亚组与D亚组重复基因是独立进化的。序列比对发现, 本文克隆的 GhVacInv2a基因属于A亚组, 前人报道的陆地棉转化酶的 GhVacInv2属于D亚组。

2.3 GhVacInv2a基因的分子进化分析

为进一步明确 GhVacInv2a与已报道的棉花和其他植物中亲缘关系相近的基因间进化关系, 使用HMMER3.0软件在雷蒙德氏棉( Gossypium raimondii) 基因组(http://www.phytozome.net/cotton.php)中发现

图1 GhVacInv2a类同源基因在不同棉种中的系统进化A: 阿非利加棉; D: 雷蒙德氏棉; TM-1_At: TM-1的A亚组; TM-1_Dt: TM-1的D亚组; H7124_At: 海7124的A亚组; H7124_Dt: 海7124的D亚组。GhVIN2代表 GhVacInv2aFig. 1 Phylogenetic analysis for GhVacInv2aorthologs in four cotton speciesA: G. herbaceum L. var. africanum; D: G. raimondii; TM-1_At: A subgenome of G. hirsutumacc.TM-1; TM-1_Dt: D subgenome of G. hirsutumacc.TM-1; H7124_At: A subgenome of G. barbadense cv. Hai7124; H7124_Dt: D subgenome of G. barbadense cv. Hai7124. GhVIN2: abbrevation for GhVacInv2a.

14个与 GhVacInv2a功能类似的转化酶基因。同时从NCBI网站中搜索了陆地棉( Gossypium hirsutum)、葡萄( Vitis vinifera)、蓖麻( Ricinus communis)、杨树( Populus trichocarpa)等植物中与 GhVacInv2a相似性较高的基因, 与本研究克隆的 GhVacInv2a进行氨基酸的多重序列比较和进化树分析。结果表明该进化树可分为两大类, 第1类为细胞壁转化酶和水解酶家族基因; 第2类为液泡转化酶; 本实验克隆的 GhVacInv2a属于第2类, 与已报道的陆地棉液泡转化酶2 ( GhVacInv2)同源性最高, 其次与蓖麻的转化酶(β-果糖呋喃苷酶)、葡萄液泡转化酶1也有较高同源性。

2.4 GhVacInv2a的表达分析

分别提取TM-1中根、茎、叶、花瓣、花药、不同发育时期纤维样品的RNA, 将其反转为cDNA, 利用实时荧光PCR (Q-PCR)扩增, 完成该基因在不同组织器官的表达分析。结果表明, GhVacInv2a在花药中的表达量最高, 其次为花瓣, 在根和茎中表达水平非常低。在纤维快速伸长期10DPA优势表达, 在纤维伸长与次生壁加厚重叠期表达水平也较高, 到纤维发育19DPA后表达水平显著降低(图3)。说明该基因主要在纤维发育的伸长及次生壁加厚早期发挥作用。

进一步分析(图4)显示, GhVacInv2a在TM-1和 im突变体中的基因表达模式类似, 但表达丰度有显著差异。尤其在13、16、19 DPA, GhVacInv2a在TM-1中的表达水平分别比 im突变体中高出28.8、13.4和8.5倍。暗示在纤维发育过程中, GhVacInv2a im

图2 来源不同植物 GhVacInv2a同源基因系统进化分析Fig. 2 Phylogenetic analysis of GhVacInv2a orthologs in different plant speciesVv: Vitis vinifera; Pt: Populus trichocarpa; Rc: Ricinus communis; Gorai: Gossypium raimondii; Gh: Gossypium hirsutum; GhVIN2: abbrevation for GhVacInv2a.

图3 GhVacInv2a基因在TM-1不同组织器官中表达水平R、S、L、P、A分别代表TM-1中根、茎、叶、花瓣、花药组织, 0代表0 DPA的胚珠, 10、13、16、19、22、25分别代表10、13、16、19、22、25 DPA的纤维组织。Fig. 3 Relative expression level of GhVacInv2a in different tissues and organs in TM-1R, S, L, P, A indicate root, stem, leaf, petal and anther tissues in TM-1, respectively. 0 indicates ovules at 0 DPA, and 10, 13, 16, 19, 22, 25 indicate fibers at 10, 13, 16, 19, 22, 25 DPA, respectively.

图4 GhVacInv2a在TM-1和 im不同纤维发育期的表达水平13、16、19、22、25分别代表13、16、19、22、25 DPA的 纤维组织。a和A分别表示0.05和0.01 显著水平上的差异。Fig. 4 Relative expression level of GhVacInv2a in different fiber developmental stages of TM-1 and im13, 16, 19, 22, 25 indicate fibers at 13, 16, 19, 22, 25 DPA, respectively. “a” and “A” indicate singnificant difference at 0.05 and 0.01 levels, respectively.

突变体中的极显著低水平表达, 可能与 im突变体的次生壁发育缺陷密切相关。

2.5 GhVacInv2a基因的染色体定位与标记性状关联分析

根据 GhVacInv2a基因在TM-1和海7124棉种间的SNP位点差异, 开发SNP引物扩增作图群体[(TM-1×H7124)×TM-1]BC1, 根据群体多态性[ 19], 初步将 GhVacInv2a定位到棉花异源四倍体第3染色体上。进一步利用以TM-1为背景渐渗海岛棉海7124第3染色体片段的单片段置换系CSIL028与 im突变体构建的F2群体, 根据群体多态性精细定位, 将 GhVacInv2a定位于标记 GhPEL与NAU8256之间, 分别相距1.3 cM和0.3 cM。

为了进一步明确 GhVacInv2a对纤维品质性状的影响, 结合CSIL028 × imF2:3家系纤维长度(FL)、纤维强度(FS)、马克隆值(Mic)、伸长率(FE)、整齐度(FU)、短纤维指数(SFI)、成熟度(Mat)和衣分(LP) 7个纤维品质性状和1个产量性状进行标记与性状关联分析[ 20]。由于CSIL028 × im F2群体纤维品质性状和衣分表型显著受 im突变基因影响, 先将分离群体按纤维表型分为正常与突变两类, 然后再分别进行单标记关联分析。在群体表型为正常类型的243个单株中, 海岛棉纯合/杂合纤维强度基因型均值为32.54 cN tex-1, 而陆地棉纯合基因型均值为 31.30 cN tex-1, 该基因与纤维强度存在极显著相关( P=0.0087), 说明 GhVacInv2a在棉花纤维强度的形成中可能发挥重要功能。

3 讨论

液泡转化酶(VIN)将液泡中的蔗糖水解为葡萄糖和果糖来调节细胞内蔗糖的浓度。许多研究报道 VIN基因在快速扩张的组织细胞, 如马铃薯块茎[ 21]、胡萝卜根[ 22]及发育中的玉米籽粒[ 23]中有较高的表达丰度或酶活性。 VIN参与调控果实及贮藏器官细胞内糖组分比例[ 4, 24, 25, 26]。通过调节渗透势 VIN参与调控细胞伸长与组织扩增[ 12, 22, 27]。研究显示 VIN还响应逆境和激素信号[ 21, 28, 29, 30, 31, 32]

本研究从棉纤维中克隆了一个编码液泡转化酶的cDNA全长基因 GhVacInv2a, 它与Taliercio等[ 13]报道的陆地棉液泡转化酶基因 GhVacInv2有很高的相似性。基因表达分析显示 GhVacInv2a除了在花药、花瓣中有很高的表达量外, 在快速伸长的纤维中具有较高的表达丰度, 在纤维伸长与次生壁加厚重叠期表达水平也较高, 而在起始和次生壁加厚期的纤维中则维持极低的表达水平。和对照TM-1相比, im突变体13~25 DPA的纤维中 GhVacInv2a始终维持一个低的表达水平。特别是在13~19 DPA的纤维中, im突变体中 GhVacInv2a的表达量比TM-1中的低8.5~28.8倍, 暗示2个材料在纤维伸长阶段可能存在不同的蔗糖利用水平。有研究报道 VIN也会受到胞内糖浓度的反馈调控[ 31, 33]。我们推测在 im突变体纤维中存在糖组分比例失衡, 细胞中没有多余的蔗糖可被转运至液泡中临时存储, 受液泡内糖浓度的反馈调控, im突变体纤维细胞中仅需维持一个低的VIN活性水平就可满足液泡中蔗糖代谢。

异源四倍体中 GhVacInv2a定位于第3染色体上, 与果胶裂解酶基因 GhPEL相距1.3 cM, 与 im基因相距约10 cM。利用基于群体的QTL分析, Wang等[ 20]在该区间内检测到几个纤维品质显著相关QTL。因此, GhVacInv2a, GhPEL im基因之间是否存在某种调控作用有必要进一步的研究。

GhVacInv2a与棉纤维强度存在极显著相关, 暗示 GhVacInv2a可能参与棉纤维强度的形成。纤维强度是棉花重要的纤维品质性状, 属于数量性状, 其形成机理较复杂。目前, 在棉花上尚未有鉴定分离控制纤维强度基因的报道。 GhVacInv2a可能参与纤维次生壁发育过程中的蔗糖转化, 与纤维强度及纤维品质密切相关。然而, GhVacInv2a在纤维发育中的确切作用, 特别是在次生壁发育缺陷的 im突变体中 GhVacInv2a在纤维发育过程中扮演的角色, 以及该基因与 im基因在纤维发育通路中的关系值得进一步深入研究。

4 结论

从遗传标准系TM-1中克隆了一个编码液泡转化酶的 GhVacInv2a。该基因包含7个外显子和6个内含子, ORF全长1857 bp, 编码618个氨基酸。在二倍体棉种阿非利加棉和雷蒙德氏棉基因组中各含1个拷贝, 在四倍体陆地棉、海岛棉基因组中存在2个拷贝, 本文克隆的 GhVacInv2a和前人报道的 GhVacInv2分属于A、D亚组同源基因。利用SNP标记将 GhVacInv2a定位于第3染色体上。 GhVacInv2a在花药中表达量最高, 在快速伸长阶段的纤维细胞和纤维伸长与次生壁加厚重叠期优势表达。13~19 DPA的纤维中, GhVacInv2a在TM-1中表达水平极显著高于 im不成熟纤维突变体。 GhVacInv2a与纤维强度存在极显著相关。推测该基因功能与纤维次生壁发育过程中的蔗糖转化有关。

The authors have declared that no competing interests exist.

作者已声明无竞争性利益关系。The authors have declared that no competing interests exist.

参考文献
[1] Sturm A, Tang G Q. The sucrose-cleaving enzymes of plants and crucial for development, growth and carbon partitioning. Trends Plant Sci, 1999, 4: 401-407 [本文引用:2] [JCR: 11.808]
[2] Roitsch T, González M C. Function and regulation of plant invertases: sweet sensations. Trends Plant Sci, 2004, 9: 606-613 [本文引用:1] [JCR: 11.808]
[3] Davies C, Robinson S P. Sugar accumulation in grape berries: cloning of two putative vacuolar invertase cDNAs and their expression in grapevine tissues. Plant Physiol, 1996, 111: 275-283 [本文引用:1] [JCR: 6.555]
[4] Greiner S, Rausch T, Sonnewald U, Herbers K. Ectopic expression of a tobacco invertase inhibitor homolog prevents cold-induced sweetening of potato tubers. Nat Biotechnol, 1999, 17: 708-711 [本文引用:2] [JCR: 32.438]
[5] Klann E M, Chetelat R T, Bennett A B. Expression of acid invertase gene controls sugar composition in tomato (Lycopersicon) fruit. Plant Physiol, 1993, 103: 863-870 [本文引用:1] [JCR: 6.555]
[6] Ruan Y L, Jin Y, Yang Y J, Li G J, Boyer J S. Sugar input, metabolism, and signaling mediated by invertase: roles in development, yield potential, and response to drought and heat. Mol Plant, 2010, 3: 942-955 [本文引用:2] [JCR: 6.126] [CJCR: 0.9487]
[7] Kohorn B D, Kobayashi M, Johansen S, Riese J, Huang L F, Koch K, Fu S, Dotson A, Byers N. An Arabidopsis cell wall-associated kinase required for invertase activity and cell growth. Plant J, 2006, 46: 307-316 [本文引用:1] [JCR: 6.582]
[8] Roitsch T, González M C. Function and regulation of plant invertases: sweet sensations. Trends Plant Sci, 2004, 9: 606-613 [本文引用:1] [JCR: 11.808]
[9] Wind J, Smeekens S, Hanson J. Sucrose: metabolite and signaling molecule. Phytochemistry, 2010, 71: 1610-1614 [本文引用:1] [JCR: 3.05]
[10] Basra A S, Malik C. Development of the cotton fiber. Int Rev Cytol, 1984, 89: 65-113 [本文引用:1]
[11] Kim H J, Triplett B A. Cotton fiber growth in planta and in vitro. Models for plant cell elongation and cell wall biogenesis. Plant Physiol, 2001, 127: 1361-1366 [本文引用:1] [JCR: 6.555]
[12] Wang L, Li X R, Lian H, Ni D A, He Y K, Chen X Y, Ruan Y L. Evidence that high activity of vacuolar invertase is required for cotton fiber and Arabidopsis root elongation through osmotic dependent and independent pathways, respectively. Plant Physiol, 2010, 154: 744-756 [本文引用:2] [JCR: 6.555]
[13] Taliercio E, Scheffler J, Scheffler B. Characterization of two cotton (Gossypium hirsutum L. ) invertase genes. Mol Biol Rep, 2010, 37: 3915-3920 [本文引用:2] [JCR: 2.506]
[14] Eddy S R. Accelerated profile HMM searches. PLoS Comput Biol, 2011, 7: e1002195 [本文引用:1] [JCR: 4.867]
[15] 武耀廷, 刘进元. 一种高效提取棉花不同组织总RNA的热硼酸改良法. 棉花学报, 2004, 16(2): 67-71
Wu Y T, Liu J Y. A modified hot borate method for efficient isolation of total RNA from different cotton tissues. Cotton Sci, 2004, 16(2): 67-71 (in Chinese with English abstract) [本文引用:1] [CJCR: 1.462]
[16] 蒋建雄, 张天真. 利用CTAB/酸酚法提取棉花组织总RNA. 棉花学报, 2003, 15: 166-167
Jiang J X, Zhang T Z. Extraction of total RNA in cotton tissues with CTAB-acidic phenolic method. Cotton Sci, 2003, 15: 166-167 (in Chinese with English abstract) [本文引用:1] [CJCR: 1.462]
[17] Paterson A H, Brubaker C L, Jonathan F. Rapid method for extraction of cotton (Gossypium spp. ) genomic DNA suitable for RFLP or PCR analysis. Plant Mol Biol Rep, 1993, 11: 122-127 [本文引用:1] [JCR: 5.319]
[18] Pfaffl M W. A new mathematical model for relative quantification in real-time RT-PCR. Nucl Acids Res, 2001, 29: e45 [本文引用:1]
[19] Zhao L, Lv Y D, Cai C P, Tong X C, Chen X D, Zhang W, Du H, Guo X H, Guo W Z. Toward allotetraploid cotton genome assembly: integration of a high-density molecular genetic linkage map with DNA sequence information. BMC Genomics, 2012, 13: 539 [本文引用:2] [JCR: 4.397]
[20] Wang C, Zhang T Z, Guo W Z. The im mutant gene negatively affects many aspects of fiber quality traits and lint percentage in cotton. Crop Sci, 2013, 53: 27-37 [本文引用:3] [JCR: 1.513]
[21] Ross H A, Davies H V, Burch L R, Viola R, McRae D. Developmental changes in carbohydrate content and sucrose degrading enzymes in tuberising stolons of potato (Solanum tuberosum). Physiol Plant, 1994, 90: 748-756 [本文引用:2] [JCR: 6.555]
[22] Tang G Q, Luscher M, Sturm A. Antisense repression and vacuolar and cell wall invertase in transgenic carrot alters early plant development and sucrose partitioning. Plant Cell, 1999, 11: 177-189 [本文引用:2] [JCR: 9.251]
[23] Andersen M N, Asch F, Wu Y, Jensen C R, Naested H, Mogensen V O, Koch K E. Soluble invertase expression is an early target of drought stress during the critical, abortion-sensitive phase of young ovary development in maize. Plant Physiol, 2002, 130: 591-604 [本文引用:1] [JCR: 6.555]
[24] Klann E M, Hall B, Bennett A B. Antisense acid invertase (TIV1) gene alters soluble sugar composition and size in transgenic tomato fruit. Plant Physiol, 1996, 112: 1321-1330 [本文引用:1] [JCR: 6.555]
[25] Davies C, Robinson S P. Sugar accumulation in grape berries: cloning of two putative vacuolar invertase cDNAs and their expression in grapevine tissues. Plant Physiol, 1996, 111: 275-283 [本文引用:1] [JCR: 6.555]
[26] Yau Y, Simon P. A 2. 5-kb insert eliminates acid soluble invertase isozyme II transcript in carrot (Daucus carota L. ) roots, causing high sucrose accumulation. Plant Mol Biol, 2003, 53: 151-162 [本文引用:1] [JCR: 3.518]
[27] Sergeeva L I, Keurentjes J J B, Bentsink L, Vonk J, van der Plas L H W, Koornneef M, Vreugdenhil D. Vacuolar invertase regulates elongation of Arabidopsis thaliana roots as revealed by QTL and mutant analysis. Proc Nat Acad Sci USA, 2006, 103: 2994-2999 [本文引用:1]
[28] Long J C, Zhao W, Rashotte A M, Muday G K, Huber S C. Gravity-stimulated changes in auxin and invertase gene expression in maize pulvinal cells. Plant Physiol, 2002, 128: 591-602 [本文引用:1] [JCR: 6.555]
[29] Trouverie J, Thévenot C, Rocher J P, Sotta B, Prioul J L. The role of abscisic acid in the response of a specific vacuolar invertase to water stress in the adult maize leaf. J Exp Bot, 2003, 54: 2177-2186 [本文引用:1] [JCR: 5.242]
[30] Mitsuhashi W, Sasaki S, Kanazawa A, Yang Y Y, Kamiya Y, Toyomasu T. Differential expression of acid invertase genes during seed germination in Arabidopsis thaliana. Biosci Biotechnol Biochem, 2004, 68: 602-608 [本文引用:1]
[31] Trouverie J, Chateau-Joubert S, Thévenot C, Jacquemot M P, Prioul J L. Regulation of vacuolar invertase by abscisic acid or glucose in leaves and roots from maize plantlets. Planta, 2004, 219: 894-905 [本文引用:2] [JCR: 3.347]
[32] Li Z, Palmer W M, Martin A P, Wang R, Rainsford F, Jin Y, Patrick JW, Yang Y, Ruan Y L. High invertase activity in tomato reproductive organs correlates with enhanced sucrose import into, and heat tolerance of, young fruit. J Exp Bot, 2012, 63: 1155-1166 [本文引用:1] [JCR: 5.242]
[33] Xu J, Avigne W T, McCarty D R, Koch K E. A similar dichotomy of sugar modulation and developmental expression affects both paths of sucrose metabolism: evidence from a maize invertase gene family. Plant Cell, 1996, 8: 1209-1220 [本文引用:1] [JCR: 9.251]