棉花COL家族基因的鉴定、表达与进化分析

1、棉花COL家族基因的鉴定、表达与进化分析棉花是全球性的重要经济作物之一,植物学分类属于被子植物门(Angiospermae)、双子叶植物纲(Dicotyledoneae)、锦葵目(Malvales)、锦葵科(Malvaceae)、棉族(Gossypieae)、棉属(Gossypium)。棉属共50个种,包括45个二倍体及5个四倍体种。四倍体棉种包括3个野生棉,即毛棉、黄褐棉及达尔文氏棉,此外陆地棉和海岛棉是目前世界上栽培面积最广的两个栽培种。栽培种由野生种驯化而成,由于长期栽培、驯化及选择作用,野生棉伴随着对光周期敏感的多年生短日照生长习性转变为光周期不敏感的一年生生长习性。发掘棉花光周期相关

2、基因,研究光周期开花途径的作用机制,对于棉花适应多样性的环境及种植面积的扩展十分有益。CO(CONSTAN)COL(CONSTANS-LIKE)基因是调控植物光周期开花的一个重要转录因子。CCS上游生物钟相关基因(GIGANTEA)G的影响，并激活下游的成花素FT(FLOWERINGLOCUST滨达.对于特定基因或位点的选择作用，会导致栽培种群体内核酸多样性水平降低、基因频率的降低或升高而偏离中性期望及连锁不平衡的增加。目前水稻中已经证明COL因是选择的靶标基因。鉴于CO瑶族基因在光周期开花途径中的保守功能，棉花中CO瑶族基因信息的局限性以及驯化选择中作用未知。本研究根据拟南芥、水稻的COLB

3、族基因信息从全基因组水平分离棉花的CO瑶族基因，对其进行分类、结构特征、定位、组织表达、昼夜节律表达及分子进化等系统分析,主要研究结果如下:1.利用比较基因组学，根据拟南芥、水稻的COL家族基因信息，比较雷蒙德氏棉数据库，鉴定了23个棉花COLa族基因根据系统进化分析，该家族基因分为3类，其中I类基因包括COL1-8,均含有两个外显子和一个内含子。AA分析显示有2个B-box及1个CC琢才域（COL8除外）；II类基因包括COL9-11，也包含两个外显子和一个内含子，有1个B-box及1个CCT域；m类基因包括COL12-23，该类基因包含一个完整的B-box,一个分化的锌指结构及一个CC锚构

4、域。该类基因结构比较多变,包括2个或3个内含子。2.利用实时荧光定量PCR技术，研究了20个COLa族基因在棉花TM-1根、茎、叶、花药、花瓣等不同组织及0DPA10DPAR20DPA千维发育不同阶段的表达特性。结果显示，I类中的八个基因中，COL1-5,COL8等六个基因在叶中高表达,COL6和COL7t花中优势表达。此外，COL2-5在纤维中表达也较高。II类中三个基因均在叶中优势表达。田类中九个基因，除COL14COL2吸COL2汾别在茎、叶及纤维中优势表达外，其余六个基因均是组成型表达。3 .在LD或SD条件下,分别处理陆地棉TM-1及海岛棉H7124，连续48h取三叶一心期叶片进行C

5、O瑶族基因表达特征分析。结果显示已分析的20个基因中有18个基因的表达表现为昼夜节律性，表明棉花COIB族基因参与棉花光周期调控的开花途径,在棉花的光周期开花途径中起重要作用。4 .选择包括野生棉、陆地棉及海岛棉等25个棉材料及外类群长梗肖槿为研究材料，对CO/Hd1所在的I类中8个基因进行分子进化相关分析。通过各基因基因组与cDNA行多序歹I比对分析，结果表明这8个COLS因非常保守，均含有两个外显子和一个内含子。8个COL!因其DN睑长变异范围为1,030bp-1,611bp,内含子长度的变异范围为77bp-680bp。通过比较不同棉种间同一亚组其同源基因的变异发现,COL2,COL6及C

6、OL8ft供试材料间的内含子或外显子H区发生插入/缺失,导致长度变异。而同一材料其A-及D亚组部分同源基因的结构变异表明，COL4和COL琳外显子I或外显子H发生插入/缺失，导致At及Dt间长度变异。5.系统进化分析表明,A组供体阿非利加棉与异源四倍体棉A亚组聚在一起,D组供体与异源四倍体棉D亚组聚为一类。184个成对比较中,其Avs.D、Atvs.Dt成对比较的进化速率高于Avs.At及Dvs.Dt进化速率占到98.37%。进一步分析证明Avs.D与Atvs.Dt进化速率之间表现极高的正相关,最小的相关系数r值为0.797。说明异源四倍体中,8个COLB族基因的A-及D-亚组间是独立进化的。

7、6.通过成对比较8个COLa族基因组合，以及每个基因在所有四倍体棉种及AD基因组间的核酸多样性(冗)。结果表明8个基因组合中,DvsDt核酸多样性(0.01051),显著大于AvsAt的核酸多样性(0.00586)(P=4.9E-21)。单个基因分析表明，COL2-COL5,COL及COL8a6个基因在Dt中的进化速率显著快于At,COL6是在At的进化速率显著快于Dt,而COL仰因的进化速率在At和Dt之间无显著差异。综上所述,8个COLB族基因在At及Dt间的进化速率不同，大部分基因其Dt的进化速率快于Ato7.分别分析野生棉、陆地棉及海岛棉的8个COL1因位点其核酸多样性冗total,结

8、果显示野生棉的冗total,显著高于陆地棉(0.00369vs0.00139)(P=0.001)及海岛棉(0.00369vs0.00035)(P=7.3E-7),表明陆地棉及海岛棉在驯化过程中存在瓶颈效应。进一步比较B-box,Var及CCT&三个结构域在不同棉种中的核酸多样性，发现B-box(0.00284vs0.00119,P=0.028)及Var(0.00243vs0.00119,P=0.014)结构域的核酸多样性显著高于CCTM才域,B-box及Var结构域间的核酸多样性无显著差异。说明B-box及Var结构域在进化过程中发生较大变异,而CC锚构域在不同棉种间高度保守,功能上也比较保守。8.通过对8个CO瑶族基因进行TajimasD,FuandLisD和F中性检验,表明陆地棉在COL1MCOL型存在大量的低频率位点，这与陆地棉最近发生的正选择作用