删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

甘蓝型油菜BnGS3和BnGhd7的同源克隆及其与油菜产量相关性状的关系

本站小编 Free考研考试/2021-12-26

薛志飞, 王夏, 李付鹏, 马朝芝*,
华中农业大学植物科学技术学院 / 作物遗传改良国家重点实验室 / 国家油菜工程技术研究中心, 湖北武汉 430070;

Homologous Cloning of BnGS3 and BnGhd7 Genes in Brassica napus and Their Relationship with Rapeseed Yield-related Traits

XUEZhi-Fei, WANGXia, LIFu-Peng, MAChao-Zhi*,
College of Plant Science and Technology of Huazhong Agricultural University / National Key Laboratory of Crop Genetic Improvement / National Research Center of Rapeseed Engineering and Technology, Wuhan 430070, Hubei, China;
通讯作者:马朝芝, E-mail: yuanbeauty@mail.hzau.edu.cn
收稿日期:2017-04-19
接受日期:2017-11-21
网络出版日期:2018-02-12
版权声明:2018作物学报编辑部作物学报编辑部
基金资助:本研究由国家重点研发计划专项(2016YFD0100803)资助
作者简介:
-->zfxue@webmail.hzau.edu.cn, Tel: 027-87281807



展开

摘要
甘蓝型油菜为异源四倍体, 基因组结构复杂, 而水稻基因组与油菜基因组具有一定的共线性。本研究利用水稻产量相关基因GS3Ghd7序列信息, 通过同源克隆方法, 获得了甘蓝型油菜的同源基因BnGS3BnGhd7BnGS3有6个外显子, ORF全长666 bp, 编码222个氨基酸。BnGS3蛋白具有水稻GS3四个保守结构域中的VWF结构, 属于A型。BnGhd7含有1个外显子, ORF全长1014 bp, 编码337个氨基酸。BnGhd7蛋白具有N端的B-Box和C端的CCT两个重要的结构域。BnGS3BnGhd7分别位于A2和A10连锁群, 比较测序得到BnGS3的多态性标记brgs-16及BnGhd7的多态性标记brghd-3和ghd7-7, 其中brghd-3与千粒重(P < 0.05)、ghd7-7与株高(P < 0.01)正相关, ghd7-7与开花期负相关(P < 0.05)。上述结果表明, 利用水稻功能基因序列信息克隆甘蓝型油菜的同源基因是可行的, 为油菜功能基因研究提供了一种有效途径。

关键词:甘蓝型油菜;同源克隆;产量;比较测序;相关性分析
Abstract
Brassica napus is an allotetraploid with complex genomic structure, but Brassica genome has a collinearity with that of rice. In this study, the homologous genes BnGS3 and BnGhd7 of Brassica napus were obtained by homologous cloning of the yield related genes GS3 and Ghd7 of rice. BnGS3 has six exons with a 666 bp of ORF and encodes 222 amino acids. BnGS3 protein has VWF structure, one of the four conserved domains of rice GS3, and belongs to type A. BnGhd7 contains one exon having ORF with a full-length of 1014 bp and encoding 337 amino acids. BnGhd7 protein has two important domains, the N-terminal B-Box and the C-terminal CCT. BnGS3 and BnGhd7 were located in linkage groups A2 and A10, respectively. Polymorphism markers brgs-16 of BnGS3 and polymorphism markers brghd-3 and ghd7-7 of BnGhd7 were obtained by comparative sequencing, among which brghd-3 with thousand kernel weight (P < 0.05) and ghd7-7 with plant height (P <0.01) were positively correlated, and ghd7-7 was negatively correlated with flowering stage (P < 0.05). The results indicate that it is feasible to clone the homologous genes of Brassica napus using rice functional gene sequence information, which provides an effective way in Brassica napus functional gene research.

Keywords:Brassica napus;homologous cloning;yield;comparative sequencing;correlation analysis

-->0
PDF (1668KB)元数据多维度评价相关文章收藏文章
本文引用格式导出EndNoteRisBibtex收藏本文-->
薛志飞, 王夏, 李付鹏, 马朝芝. 甘蓝型油菜BnGS3BnGhd7的同源克隆及其与油菜产量相关性状的关系[J]. 作物学报, 2018, 44(02): 297-305 https://doi.org/10.3724/SP.J.1006.2018.00297
XUE Zhi-Fei, WANG Xia, LI Fu-Peng, MA Chao-Zhi. Homologous Cloning of BnGS3 and BnGhd7 Genes in Brassica napus and Their Relationship with Rapeseed Yield-related Traits[J]. Acta Agronomica Sinica, 2018, 44(02): 297-305 https://doi.org/10.3724/SP.J.1006.2018.00297
油菜是世界上最主要的油料作物之一, 提高产量始终是油菜育种的一个重要目标。研究产量相关性状、发展分子标记进行辅助选择将是提高产量的有效途径。现已完成甘蓝型冬油菜品种Tapidor全基因组测序, 预测有101 040个基因, 在测序植物中基因数目最多[1]。甘蓝型油菜(Brassica napus, AACC=38)是由白菜(Brassica rapa, AA=20)和甘蓝(Brassica oleracea, CC=18)自然杂交、加倍而成的异源四倍体[2], A/C基因组在进化过程中发生了大量的重组和突变, 增加了基因组结构的复杂性, 导致甘蓝型油菜功能基因组研究难度加大。油菜与水稻(Oryza sativa)基因组虽然差异大, 但仍然具有共线性[3]。作为高等模式生物, 水稻的大量功能基因已被克隆, 利用水稻基因序列克隆油菜同源基因, 很可能是研究油菜基因功能的有效手段。
无论水稻还是油菜, 生育期、千粒重和籽粒大小都是重要的产量相关性状, 控制这些性状的很多基因或QTL已在水稻中被克隆。其中, GS3是控制粒长与千粒重的主效QTL, 同时也是粒宽和粒厚相关的微效QTL[4]; Ghd7同时影响株高、穗粒数和抽穗期等多个性状, 是一个控制水稻产量的多效基因[5]GS3全长5363 bp, 其编码蛋白包含4个保守结构域, 分别是磷脂酞乙醇胺结合蛋白区域(phosphatidylethanolamine-binding protein, PEBP-like)、跨膜区域、生长因子受体家族富半胱氨酸同源区域(tumour necrosis factor receptor, TNFR或nerve growth factor receptor, NGFR)和von Willebrand因子C型结构域(von Willebrand factor type C, VWFC)[6]Ghd7编码257个氨基酸残基, 其蛋白C端具有一个明显的CCT (constans, constans-like, TOCl)结构域, 对控制植物花期的重要基因CO发挥核心功能, Ghd7缺少位于CO基因N端的另一个重要的B-BOX结构域[7]。甘蓝型油菜是否有GS3Ghd7的同源基因, 其同源基因是否与生育期和千粒重等产量性状有关, 目前还未见相关的研究报道。
本研究根据水稻GS3Ghd7基因序列, 通过比较基因组学方法, 克隆甘蓝型油菜的同源基因, 并对候选基因进行比较测序, 分析其与产量性状的相关性, 目的在于探究利用模式植物水稻功能基因信息研究油菜相关基因的可行性。

1 材料与方法

1.1 植物材料

选用58份甘蓝型油菜自交系和由Tapidor × Ningyou 7创建而来具有202个DH系(double haploid)的TN群体[8]。58份材料中, 自交不亲和系SI-1300用于同源基因克隆, 57份自交亲和系用于表型性状的相关分析, TN群体用于连锁图谱分析和基因定位。

1.2 表型调查与分析

2009年10月上旬, 将57份自交亲和系分别在湖北武汉和黄冈两地种植, 均采用随机区组设计, 3行区, 行距0.27 m, 行长3.50 m, 小区面积2.84 m2, 定苗后每行20株、每小区60株。记录初花期(小区内50%植株开第1朵花)和终花期日期, 分别计算播种-初花期、初花期-终花期的天数。成熟后随机连续取10株, 考察株高、一次有效分枝数、单株角果总数、角果粒数、单株产量、千粒重和小区产量。
利用SAS 9.4软件进行方差分析(ANOVA)以及表型与基因型的简单相关分析(simple correlation analysis)。

1.3 引物设计及同源基因克隆

使用天根公司植物基因组DNA提取试剂盒(OSR- M301)提取油菜苗期幼嫩叶片总DNA。
使用软件Oligo v. 7设计同源克隆引物。根据水稻GS3的DNA序列(GenBank登录号为AB743898.1)设计引物bngs-1-1、bngs-1-2和bngs-1-3 (见附表1), 但是该3对引物在甘蓝型油菜中没有特异扩增, 因此利用GS3蛋白序列重新进行同源搜索, 获得同源的拟南芥NP_680175.2蛋白序列, 利用其对应的基因序列信息在NCBI中再次进行同源比对, 获得白菜型油菜的同源序列AC189411.2。根据拟南芥NP_680175.2基因序列设计引物gs3-1和gs3-4, 根据白菜型油菜AC189411.2的基因序列设计了3对引物, 分别是gs3-2、gs3-3和gs3-5 (见附表1)。
Supplementary table 1
附表1
附表1同源克隆BnGS3BnGhd7的引物
Supplementary table 1Primers used in homologous clonging of BnGS3 and BnGhd7
引物
Primer
正向序列
Forward sequence (5'-3')
反向序列
Reverse sequence (5'-3')
退火温度
Tm (°C)
bngs-1-1TACTCATCCTCCTCCTCCCGCTCATCATCGAGAACT55
bngs-1-2GCGTAATGCCTCAGCTACGACCCATCATCATCGAGA55
bngs-1-3GCATGATGCTTTCACCACCGATCATCATCGAGAACT55
gs3-1TGTTAGGGTTTCTGTTGGTGGCCAGATGCTGCAAAGAGTAAG59
gs3-2GGGTTTCTGTTGGTGGGCTTCCCATGAGGTATGTCAGCAT59
gs3-3TTCTAGACAATGAATGGCAGGACGAAAGAGTCGACGACACTG58
gs3-4TTCTAGACAATGAATGGCAGTGCAACTGCAAGATCAAATG59
gs3-5CCGCAGCGGAAAAGGTATGTGTCCAGATGCTGCAAAGAGT59
ghd7-2-4CCAAAAGCCAACGTCACCATGCGATCAGCGACCATTAAAG56
ghd7-2-5CCAAAAGCCAACGTCACCATGGAGCGACCGAAAACTACAT54
ghd7-2-6TCCCTCACCAACAACAAACCATCCTTGGTCTTTTCTCTGC55
ghd7-1CAAAAGCCAACGTCACCATCCGACCATTAAAGAACAGGCT58
ghd7-2TTCGGATTCGGTTCTGGTTCGCTCACATGATTGACAGACT59
ghd7-3AGCAGAGGCGGCTTCTTGGTGCGATCAGCGACCATTAAAG59
ghd7-4CCAAAAGCCAACGTCACCATGCGATCAGCGACCATTAAAG58
ghd7-5CCAAAAGCCAACGTCACCATGGAGCGACCGAAAACTACAT59
ghd7-6TCCCTCACCAACAACAAACCATCCTTGGTCTTTTCTCTGC59
ghd7-7GGGTTGTTCCACTTCAGGTTACCCATGGAAAGTGGTAGAT59
ghd7-8CAAAAGCCAACGTCACCATCCATCCGATATTTTTGTCTCC60
ghd7-9GATCCATCACGGGCCATAACGGGTTTGAAACTGTTGTCTC59


新窗口打开
水稻Ghd7具有保守的CCT结构域, CCT结构域是CO基因家族编码蛋白的核心功能域。同源比对发现Ghd7与甘蓝型油菜AY280868.2序列和拟南芥CO基因具有相似的同源性。根据水稻Ghd7 (GenBank登录号为JF926542.1)、甘蓝型油菜AY280868.2的序列信息设计12对引物(见附表1)。
PCR体系使用1×扩增缓冲液配制, 加入2 mmol L-1的MgCl2、0.1 mmol L-1的dNTP、左右引物各0.2 μmol L-1、0.5 U Taq酶、50 ng模板DNA, ddH2O定容至20 μL。PCR反应程序为94°C运行2 min, 94°C运行1 min, 55°C运行30 s, 72°C运行45 s, 35个循环; 72°C运行5 min; 于4°C保存产物。

1.4 候选基因及其编码产物的序列分析

用SSCP (single-strand conformation polymorphism)技术可以有效检测InDel和SNP所造成的序列差异, 但PCR产物需要经高温变性后才能检测到单链折叠的构象差异, 已经被用来开发甘蓝型油菜分子标记[9]。考虑到SSCP对检测片段大小有严格要求, 在8%聚丙烯酰胺凝胶的基础上, 又配制了3.5%和6%非变性聚丙烯酰胺凝胶, 分别检测>800 bp和400~500 bp的片段。使用天根公司纯化回收试剂盒(DP214)回收PCR产物, 使用pMD18-T Vector (TaKaRa)进行T-A克隆, 挑选转化的大肠杆菌感受态(K12菌株制备)阳性单克隆送北京擎科新业生物技术公司测序, 使用软件SeqMan Pro (http://www.dnastar.com/)分析测序结果。
通过数据库NCBI (https://www.ncbi.nlm.nih.gov/)分别对GS3Ghd7的序列进行同源比对。用ClustalW完成多序列对位排列分析[10]; 利用在线软件Softberry (http:// www.softberry.com/) [11]和Genscan (http://genes.mit.edu/ GENSCAN.html)预测基因结构, 用在线软件InterProScan (http://www.ebi.ac.uk/interpro/search/sequence-search)[12,13]和Conserved Domain Database (https://www.ncbi.nlm. nih.gov/cdd/)预测同源蛋白结构域。

1.5 基因型分析

参照MapMaker/Exp 3.0[14,15]的方法记载分子标记的基因型。

2 结果与分析

2.1 GS3在甘蓝型油菜中的同源克隆

2.1.1 GS3在甘蓝型油菜中的同源克隆 根据水稻GS3基因序列设计的引物bngs-1-1、bngs-1-2和bngs-1-3, 在甘蓝型油菜SI-1300中为非特异扩增(图1-A), 说明直接利用水稻序列设计引物不能在油菜上得到有效扩增。而根据白菜型油菜AC189411.2序列设计的3对引物则可检测到特异扩增产物(图1-B)。扩增片段的测序结果与白菜型油菜AC189411.2序列一致, 拼接产生3040 bp的DNA片段。预测基因BnGS3位于拼接片段的1271~2965 bp处, 有6个外显子, 开放阅读框(ORF)全长666 bp, 编码222个氨基酸(图2)。
显示原图|下载原图ZIP|生成PPT
图1根据同源序列设计引物的扩增结果
A: 由水稻GS3序列设计的引物在水稻(1~3泳道)和甘蓝型油菜(4~6泳道)中的扩增结果。M: DL2000 marker; 1和4泳道: 引物bngs-1-1; 2和5泳道: 引物bngs-1-2; 3和6泳道: 引物bngs-1-3。B: 由白菜型油菜AC189411.2基因序列设计的引物在甘蓝型油菜SI-1300中的扩增结果。M: DL2000 marker; 1: 引物gs3-2; 2: 引物gs3-3; 3: 引物gs3-5。

-->Fig. 1Amplification results of primers designed according to homologous sequences
A: amplification profile in rice (lanes 1-3) and B. napus (lanes 4-6) with primers designed from rice GS3 gene sequence. M: DL2000 marker; Lanes 1 and 4: primer bngs-1-1; Lanes 2 and 5: primer bngs-1-2; Lanes 3 and 6: primer bngs-1-3. B: amplification profile in B. napus SI-1300 with primers designed from B. rapa AC189411.2 gene sequence. M: DL2000 marker; 1: primer gs3-2; 2: primer gs3-3; 3: primer gs3-5.

-->

显示原图|下载原图ZIP|生成PPT
图2BnGS3基因结构预测
TSS: 转录起始位点; CDSf: 第1个外显子翻译起始位置; CDSi: 具有多个外显子; CDSo: 只有1个外显子; CDSI: 最后1个外显子翻译终止位置; Poly A: 转录终止位置; 数字表示碱基的位置。

-->Fig. 2Prediction of BnGS3 gene construction by FGENSH 2.0
TSS: transcription start site; CDSf: start coding site of the first exon; CDSi: inter exons of coding sequence; CDSo: only one exon of coding sequence; CDSI: last exon of coding sequence; Poly A: the site of termination of transcription; numbers indicate the site of base pair.

-->

根据BnGS3序列信息进行同源搜索, 发现与甘蓝型油菜EST序列DY004631同源性较高, 而BnGS3蛋白与GS3蛋白的序列同源性较低。借助预测蛋白的保守序列进行高级结构分析, 发现BnGS3蛋白缺少跨膜区域、TNFR/NGFR和PEBP-like的功能结构域, 只具有属于VWF结构域中的VWFA型的Iso_dh super family结构, 而GS3为VWFC型。甘蓝型油菜BnGS3与水稻GS3有一定的同源性, 但序列差异大。
2.1.2 Ghd7在甘蓝型油菜中的同源克隆 以甘蓝型油菜AY280868.2为模板设计的引物ghd7-2、ghd7-5、ghd7-6、ghd7-7和ghd7-9 (见附表1), 在甘蓝型油菜SI-1300中均有一条清晰扩增带(图3)。
显示原图|下载原图ZIP|生成PPT
图3引物ghd7-2、ghd7-5、ghd7-6、ghd7-7和ghd7-9在甘蓝型油菜SI-1300中的扩增
泳道M: DNA ladder DL2000; 泳道1: ghd7-2; 泳道2: ghd7-5; 泳道3: ghd7-6; 泳道4: ghd7-7; 泳道5: ghd7-9。

-->Fig. 3Amplification of primers ghd7-2, ghd7-5, ghd7-6, ghd7-7, and ghd7-9 in B. napus SI-1300
M: DNA ladder DL2000; 1: ghd7-2; 2: ghd7-5; 3: ghd7-6; 4: ghd7-7; 5: ghd7-9.

-->

扩增片段长度为2725 bp, 测序结果与甘蓝型油菜AY280868.2一致。通过结构预测发现在714~2700 bp处存在一个候选基因BnGhd7, 含有1个外显子, ORF全长1014 bp, 编码337个氨基酸(图4)。
显示原图|下载原图ZIP|生成PPT
图4BnGhd7基因结构预测
TSS: 转录起始位点; CDSf: 第1个外显子翻译起始位置; CDSi: 具有多个外显子; CDSo: 只有1个外显子; CDSI: 最后1个外显子翻译终止位置; Poly A: 转录终止位置; 数字表示碱基的位置。

-->Fig. 4Prediction of BnGhd7 gene construction by FGENSH 2.0
TSS: transcription start site; CDSf: start coding site of the first exon; CDSi: inter exons of coding sequence; CDSo: only one exon of coding sequence; CDSI: last exon of coding sequence; Poly A: the site of termination of transcription; numbers indicate the site of base pair.

-->

根据BnGhd7的ORF序列进行同源搜索, 发现了多个高度同源的甘蓝型油菜EST序列。使用InterProScan软件对预测蛋白BnGhd7进行结构分析, 发现了2个重要的结构域, 其中N端存在1个B-Box结构域, C端存在1个CCT结构域, 与水稻Ghd7有高度同源性, 属于典型的CO基因家族。

2.2 甘蓝型油菜BnGS3BnGhd7基因的变异位点检测

2.2.1 BnGS3基因的变异位点检测 在TN群体中随机选取15份材料, 通过引物gs3-2、gs3-3和gs3-5发现疑似变异位点, 但甘蓝型油菜作为异源四倍体存在基因多拷贝现象, 直接测序不能有效发现变异位点。根据测序结果重新设计了18对引物(编号brgs-1至brgs-18), 扩增范围介于400~1500 bp, 与引物gs3-2、gs3-3和gs3-5共同覆盖BnGS3的全长并形成有效重叠(见附表2)。通过SSCP技术在57份材料中发现多个突变位点, 筛选突变频率大于5的位点, 绘制BnGS3的比较测序结果(图5)。
Supplementary table 2
附表2
附表2检测预测基因BnGS3BnGhd7变异位点所用引物
Supplementary table 2Primers used in detecting BnGS3 and BnGhd7 polymorphic alleles
引物
Primer
正向序列
Forward sequence (5'-3')
反向序列
Reverse sequence (5'-3')
退火温度
Tm (°C)
gs3-2GGGTTTCTGTTGGTGGGCTTCCCATGAGGTATGTCAGCAT59
gs3-5CCGCAGCGGAAAAGGTATGTGTCCAGATGCTGCAAAGAGT59
gs3-3TTCTAGACAATGAATGGCAGGACGAAAGAGTCGACGACACTG58
brgs-1TAGTTTGGTGCTCACTCCTGTCCTGCCATTCATTGTCTAG57
brgs-2TAGACAATGAATGGCAGGATCTACCTTCGACCACAACCAG58
brgs-3CGCCTCAACGAAGACTGGTTACGTCTGAGATGATCAAATG59
brgs-4ATGTTGCCGTTGACCACGATCTGCAAGATCAAATGGTCAT60
brgs-5TCGACTGCAACTGCGGAAACGGCGAGTTTCTGCTGTTGTT60
brgs-6CACGCTCAAACACGGGAAACCCCAATGATACCTGCGTAAG61
brgs-7TTGCCTCGCCGTCGGTTAGTCCTCGCTTCGAATCAACACG59
brgs-8CGAAGCCTCGTGGGTGTTAGGGATCCAACTTGGTGATCAG60
brgs-9GTTGGATCCGTACCGTACATAGCAATGATATCGGTTTGGT58
brgs-10TCTCAAATTTGATGGGAAGCCCCATGAGGTATGTCAGCAT59
brgs-11GCTGACATACCTCATGGGAGGTTGCTGTTGTTCTTGTCCG58
brgs-12CTTTCGTTGTCTGCCCGGATGGACCATATACATGTTCACC58
brgs-13TCTATGAAAGCAACATGACGGCACCGTTGTGTTATGTTTG60
引物
Primer
正向序列
Forward sequence (5'-3')
反向序列
Reverse sequence (5'-3')
退火温度
Tm (°C)
brgs-14TTCTAGACAATGAATGGCAGAGAACACAAAAACGTACGTC60
brgs-15TAAAGCTGAGAGCAGACTTGGTCGACCATCATGTCTTAAT58
brgs-16CAACAAAGTCGTAGCTTAGGATTCTGACCCAATGATACCT59
brgs-17CTTTCGTTGTCTGCCCGGATGGTGTTAGTTCGTGCATGTG59
brgs-18GTATATGGTCCTTTCAAACGTTGCCCACTCAAATTAATTG58
ghd7-2TTCGGATTCGGTTCTGGTTCGCTCACATGATTGACAGACT59
ghd7-5CCAAAAGCCAACGTCACCATGGAGCGACCGAAAACTACAT59
ghd7-6TCCCTCACCAACAACAAACCATCCTTGGTCTTTTCTCTGC59
ghd7-7GGGTTGTTCCACTTCAGGTTACCCATGGAAAGTGGTAGAT59
ghd7-9GATCCATCACGGGCCATAACGGGTTTGAAACTGTTGTCTC59
brghd-1TTATGTCGGGTTCGAATCGTATGGTTGGTCAAACGTGTAT60
brghd-2TTGATGCAGAGGAGGCAGACTTCACCAAGAAAGCAATCAC60
brghd-3TTATGTCGGGTTCGAATCGTATGGTTGGTCAAACGTGTAT58
brghd-7ACTTTATGTTGCCGTTGACCGCAACTGCAAGATCAAATGG60
brghd-8CATGATGGTCGACTGCAACTGCGAGTTTCTGCTGTTGTTG57


新窗口打开
显示原图|下载原图ZIP|生成PPT
图5BnGS3比较测序结果示意图
*: 点突变; -: 小片段缺失; ^: 点突变密集区。

-->Fig. 5Schematic diagram of BnGS3 comparative sequencing
*: point mutation; -: small deletion; ^: point mutations concentrated area.

-->

引物brgs-16扩增片段长851 bp, 具有突变频率最高的多态性位点, 位于BnGS3第4个外显子区域的1010 bp到1050 bp处, 属于缺失突变位点(图6)。brgs-16扩增片段完全包含在引物gs3-5的扩增区域内, 且测序结果一致, 证明结果可靠。
显示原图|下载原图ZIP|生成PPT
图6引物brgs-16扩增片段序列比较
-->Fig. 6Sequence alignment of amplified fragments in accessions by primer brgs-16
-->

2.2.2 BnGhd7基因的变异位点检测 设计可以有效覆盖目标区段的10对引物ghd7-2、ghd7-5~ghd7-7、ghd7-9、brghd-1~brghd-3、brghd-7和brghd-8 (见附表2), 在15份材料进行扩增。扩增片段测序, 仅有2个有效突变位点, 引物ghd7-7扩增片段在BnGhd7编码区下游有1个300 bp的缺失突变位点, 引物brghd-3扩增片段在BnGhd7起始位点上游700 bp处有1个突变位点(图7)。
显示原图|下载原图ZIP|生成PPT
图7BnGhd7比较测序结果示意图
*: 点突变; —: 长片段的缺失。

-->Fig. 7Schematic diagram of BnGhd7 comparative sequencing
*: point mutation; —: long fragment deletion.

-->

引物ghd7-7扩增片段长1578 bp, 该缺失位点介于BnGhd7下游非编码区的1488~1783 bp处。引物brghd-3位于BnGhd7起始位点上游700 bp处, 引物ghd7-2与brghd-3的重叠范围序列一致, 证明结果可靠(图8)。
显示原图|下载原图ZIP|生成PPT
图8不同材料中brghd-3测序结果的比较
-->Fig. 8Alignment of the sequencing results of brghd-3
-->

2.3 变异位点与产量相关性状的相关性分析

在武汉和黄冈两地, 分别对57份甘蓝型油菜自交系进行产量相关性状的表型鉴定, 除武汉点播种至初花期天数外, 各性状的变异系数都超过5% (表1)。材料之间差异明显, 能够用来进行相关性分析。
Table 1
表1
表157份甘蓝型油菜自交系9个产量相关性状的变异
Table 1Variations of nine yield-related traits in 57 B. napus inbred lines
性状
Trait
平均数
Average
标准差
SD
变异区间
Variation interval
变异系数
CV (%)
武汉点 Wuhan location
株高 Plant height (cm)146.6613.7093.53-169.609.34
一次有效分枝 Number of primary branch6.191.084.40-9.6717.37
单株角果数 Number of siliques per plant206.6045.76134.60-312.2022.15
角果粒数 Number of seeds per silique19.513.8910.30-28.1119.91
单株产量 Yield per plant (g)10.912.345.25-16.3621.42
小区产量 Plot yield (g)452.6994.12209.31-694.4820.79
千粒重 Thousand-seed weight (g)3.620.582.28-5.0716.08
播种-初花期天数 Days from sowing to initial bloom159.194.87151.00-168.003.06
初花期-终花期天数 Days from initial bloom to end flowering17.194.539.00-26.0026.37
黄冈点 Huanggang location
株高 Plant height (cm)162.2614.6089.73-191.539.00
一次有效分枝 Number of primary branch6.491.183.47-9.6718.18
单株角果数 Number of siliques per plant206.2539.89116.60-308.9319.34
角果粒数 Number of seeds per silique20.604.0012.81-32.4919.45
单株产量 Yield per plant (g)13.402.149.72-18.0215.95
小区产量 Plot yield (g)558.95102.46311.37-842.0518.33
千粒重 Thousand-seed weight (g)3.770.622.25-5.4316.46
播种-初花期天数 Days from sowing to initial bloom158.138.88135.00-171.675.61
初花期-终花期天数 Days from initial bloom to end flowering27.047.0917.00-49.0026.21


新窗口打开
分别将3个标记brgs-16、ghd7-7和brghd-3与黄冈测得的表型数据进行相关性分析。BnGS3的标记brgs-16与表型数据不存在相关性。BnGhd7的标记brghd-3与千粒重显著相关, ghd7-7与花期显著负相关, 与株高极显著相关(表2)。
Table 2
表2
表2标记与产量等相关性状之间的相关分析结果(黄冈)
Table 2Correlation coefficients between markers and traits (site: Huanggang)
性状 Traitbrgs-16brghd-3ghd7-7
株高 Plant height0.01-0.050.34**
一次有效分枝 Number of first branches0.100.080.17
单株角果数 Number of siliques per plant0.06-0.010.12
角果粒数 Number of seeds per silique-0.11-0.20-0.05
单株产量 Yield per plant0.00-0.13-0.06
小区产量 Plot yield0.18-0.050.21
千粒重 Thousand-seed weight0.120.30*0.14
播种-初花期天数 Days from sowing to initial bloom-0.06-0.130.25
初花期-终花期天数 Days from initial bloom to end flowering0.100.16-0.32*

* and ** indicate significant correlation at the 5% and 1% probability levels, respectively.***分别表示在5%和1%概率水平显著相关。
新窗口打开
以TN群体作为定位群体, 引物brgs-15和引物ghd7-7在两亲本间具有多态性扩增产物。brgs-15是1个SSCP标记, 定位于A2连锁群(见附图1); ghd7-7是1个显性标记, 定位于A10连锁群(见附图2)。因而BnGS3BnGhd7被分别定位于油菜A2和A10连锁群。
显示原图|下载原图ZIP|生成PPT
附图1brgs-15被定位在A2连锁群
-->Supplementary fig. 1brgs-15 located in linkage group A2
-->

显示原图|下载原图ZIP|生成PPT
附图2ghd7-7被定位于A10连锁群
-->Supplementary fig. 2ghd7-7 located in linkage group A10
-->

3 讨论

作为第1个基因组被完整测序的植物, 拟南芥成为遗传学、分子生物学等领域的典型研究对象[16]。对模式生物拟南芥的深入解析曾加速了水稻等一些更高等植物的研究[17]。近年来水稻作为高等模式生物成为新的研究热点并在基因组和功能基因的研究中取得了巨大成果[18]。禾本科作物水稻和十字花科芸薹属甘蓝型油菜在进化上亲缘关系较远, 在基因组结构和组成上存在很大差异, 缺乏利用水稻基因信息在甘蓝型油菜中克隆同源基因的报道, 但比较基因组学研究结果表明, 两种作物基因组仍然具有共线性。本研究利用水稻GS3Ghd7的序列设计引物扩增油菜基因组DNA, 没有扩增产物或无特异扩增。因而, 尝试将拟南芥作为同源克隆的媒介, 首先通过GS3的序列信息分别在拟南芥和白菜中发现了同源序列NP_ 680175.2AC189411.2, 分别设计引物, 但是, 只有依据白菜AC189411.2设计引物gs3-2、gs3-3和gs3-5在甘蓝型油菜SI-1300中能有效扩增。同样地, 根据与Ghd7同源的甘蓝型油菜AY280868.2序列设计引物, 在SI-1300有特异扩增带。说明利用远缘物种如水稻的功能基因进行油菜基因功能研究时, 借助拟南芥为媒介, 搜索油菜或其两个祖先种甘蓝和白菜序列的同源片段, 设计引物进行同源克隆, 是一种有效的途径。
本研究克隆到油菜基因BnGS3BnGhd7, 比较基因变异位点得到BnGS3的多态性标记brgs-16、Bghd7的多态性标记brghd-3和brgs-16。BnGS3有6个外显子, ORF全长666 bp, 编码222个氨基酸。水稻GS3全长5363 bp, 包含4个内含子和5个外显子, 阅读框长699 bp。GS3编码的蛋白由232个氨基酸组成, 包含4个保守的结构域, 分别为磷脂酞乙醇胺结合蛋白区域(PEBP-like)、跨膜区域、生长因子受体家族富半胱氨酸同源区域(TNFR)/ (NGFR)和von Willebrand因子C型结构域(VWFC)[19,20]。BnGS3虽然只具有水稻GS3四个保守结构域中的von Willebrand factor(VWF)结构, 且属于A型(水稻为C型), 但该结构域已被发现存在于许多胞外蛋白中, 参与调节生长因子信号传导[21,22]BnGS3的标记brgs-16与产量性状不存在相关性。因而甘蓝型油菜BnGS3与水稻GS3有一定的同源性, 但基因功能上的相似性需要进一步验证。
BnGhd7含有1个外显子, ORF全长1014 bp, 编码337个氨基酸。BnGhd7蛋白具有N端的B-Box和C端的CCT两个重要的结构域, 据此认为属于CO基因家族, 与水稻Ghd7有高度同源性, 植物中CO基因家族在调控开花期上发挥重要功能[23,24]。相关性分析发现BnGhd7不仅与花期相关, 也与株高和千粒重相关, BnGhd7很可能与水稻Ghd7具有相似的功能。
BnGS3BnGhd7被分别定位在油菜A2连锁群、A10连锁群。A2连锁群存在大量产量及相关性状的QTL, 包括花期、株高、全株角果数、角果粒数、千粒重、种子产量等[25]。由于brgs-15标记位于连锁群末端, 与A2连锁群包含的产量性状相关性低, 因而BnGS3的标记brgs-16与产量性状不存在相关性。ghd7-7被定位到A10连锁群, 该区域同样存在大量产量及相关性状QTL, 包括产量、花期、成熟期、株高、全株角果数、千粒重和种子产量等性状[26,27]。A10连锁群上与花期相关的QTL区间与拟南芥第5染色体顶端的区间具有共线性, 在该区域具有拟南芥开花基因COFYFLC等, 而ghd7-7也被定位于该区间内。
The authors have declared that no competing interests exist.
作者已声明无竞争性利益关系。

参考文献 原文顺序
文献年度倒序
文中引用次数倒序
被引期刊影响因子

[1]Chalhoub B, Denoeud F, Liu S, Parkin I A, Tang H, Wang X, Chiquet J, Belcram H, Tong C, Samans B, Corréa M.Early allopolyploid evolution in the post-neolithicBrassica napus oilseed genome
. Science, 2014, 345: 950-953
[本文引用: 1]
[2]Nagaharu U.Genome analysis in Brassica with special reference to the experimental formation ofB. napus and peculiar mode of fertilization
. Jpn J Bot, 1935, 7: 389-452
[本文引用: 1]
[3]Devos K M, Gale M D.Genome relationships: the grass model in current research
.Plant Cell, 2000, 12: 637-646
[本文引用: 1]
[4]Fan C, Xing Y, Mao H, Lu T, Han B, Xu C, Li X, Zhang Q.GS3, a major QTL for grain length and weight and minor QTL for grain width and thickness in rice, encodes a putative transmembrane protein
. Theor Appl Genet, 2006, 112: 1164-1171
[本文引用: 1]
[5]Xue W, Xing Y, Weng X, Zhao Y, Tang W, Wang L, Zhou H, Yu S, Xu C, Li X, Zhang Q.Natural variation inGhd7 is an important regulator of heading date and yield potential in rice
. Nat Genet, 2008, 40: 761-767
[本文引用: 1]
[6]Mao H, Sun S, Yao J, Wang C, Yu S, Xu C, Li X, Zhang Q.Linking differential domain functions of theGS3 protein to natural variation of grain size in rice
. Proc Natl Acad Sci USA, 2010, 107: 19579-19584
[本文引用: 1]
[7]Wang C, Chen S, Yu S.Functional markers developed from multiple loci inGS3 for fine marker-assisted selection of grain length in rice
. Theor Appl Genet, 2011, 122: 905-913
[本文引用: 1]
[8]Trick M, Long Y, Meng J, Bancroft I.Single nucleotide polymorphism (SNP) discovery in the polyploidBrassica napus using Solexa transcriptome sequencing
. Plant Biotechnol J, 2009, 7: 334-346
[本文引用: 1]
[9]李媛媛, 陈庆芳, 傅廷栋, 马朝芝. 利用SSCP技术分析甘蓝型油菜10个功能基因序列差异
. 作物学报, 2012, 38: 43-49
[本文引用: 1]

Li Y Y, Chen Q F, Fu T D, Ma C Z.Polymorphism analysis of ten functional genes inBrassica napus using SSCP method
. Acta Agron Sin, 2012, 38: 43-49 (in Chinese with English abstract)
[本文引用: 1]
[10]Thompson J D, Gibson T, Higgins D G.Multiple sequence alignment using ClustalW and ClustalX. In: Current Protocols in Bioinformatics, John Wiley and Sons, 2002. pp 2-3 [本文引用: 1]
[11]Shahmuradov I A, Gammerman A J, Hancock J M, Bramley P M, Solovyev V V.PlantProm: a database of plant promoter sequences
.Nucl Acids Res, 2003, 31: 114-117
[本文引用: 1]
[12]Zdobnov E M, Apweiler R.InterProScan: an integration platform for the signature-recognition methods in InterPro
.Bioinformatics, 2001, 17: 847-848
[本文引用: 1]
[13]Quevillon E, Silventoinen V, Pillai S, Harte N, Mulder N, Apweiler R, Lopez R.InterProScan: protein domains identifier
.Nucleic Acids Res, 2005, 33: W116-W120
[本文引用: 1]
[14]Lander E S, Green P, Abrahamson J, Barlow A, Daly M J, Lincoln S E, Newburg L.MAPMAKER: an interactive computer package for constructing primary genetic linkage maps of experimental and natural populations
.Genomics, 1987, 1: 174-181
[本文引用: 1]
[15]Lincoln S E, Daly M J, Lander E S.Constructing genetic linkage maps with MAPMAKER/EXP Version 3.0: a tutorial and reference manual
. In: A Whitehead Institute for Biomedical Research Technical Report, 1993. pp 78-79
[本文引用: 1]
[16]Arabidopsis Genome Initiative.Analysis of the genome sequence of the flowering plantArabidopsis thaliana
. Nature, 2000, 408:796
[本文引用: 1]
[17]Meyerowitz E M.Prehistory and history of Arabidopsis research
.Plant Physiol, 2001, 125: 15-19
[本文引用: 1]
[18]Zhang Q, Li J, Xue Y, Han B, Deng X W.Rice 2020: a call for an international coordinated effort in rice functional genomics
.Mol Plant, 2008, 1: 715-719
[本文引用: 1]
[19]Mao H, Sun S, Yao J, Wang C, Yu S, Xu C, Li X, Zhang Q.Linking differential domain functions of theGS3 protein to natural variation of grain size in rice
. Proc Natl Acad Sci USA, 2010, 107: 19579-19584
[本文引用: 1]
[20]Wang C, Chen S, Yu S.Functional markers developed from multiple loci inGS3 for fine marker-assisted selection of grain length in rice
. Theor Appl Genet, 2011, 122: 905-913
[本文引用: 1]
[21]Abreu J G, Coffinier C, Larraın J, Oelgeschläger M, De Robertis E M. Chordin-like CR domains and the regulation of evolutionarily conserved extracellular signaling systems
.Gene, 2002, 287: 39-47
[本文引用: 1]
[22]O’Leary J M, Hamilton J M, Deane C M, Valeyev N V, Sandell L J, Downing A K. Solution structure and dynamics of a prototypical chordin-like cysteine-rich repeat (von Willebrand Factor type C module) from collagen IIA
.J Biol Chem, 2004, 279: 53857-53866
[本文引用: 1]
[23]Nemoto Y, Kisaka M, Fuse T, Yano M, Ogihara Y.Characterization and functional analysis of three wheat genes with homology to the CONSTANS flowering time gene in transgenic rice
.Plant J, 2003, 36: 82-93
[本文引用: 1]
[24]Miller T A, Muslin E H, Dorweiler J E.A maize CONSTANS- like gene,conz1, exhibits distinct diurnal expression patterns in varied photoperiods
. Planta, 2008, 227: 1377-1388
[本文引用: 1]
[25]Li Y, Shen J, Wang T, Chen Q, Zhang X, Fu T, Meng J, Tu J, Ma C.QTL analysis of yield-related traits and their association with functional markers in Brassica napus L
. Aust J Agric Res, 2007, 58: 759-766
[本文引用: 1]
[26]Shi J, Li R, Qiu D, Jiang C, Long Y, Morgan C, Bancroft I, Zhao J, Meng J.Unraveling the complex trait of crop yield with quantitative trait loci mapping inBrassica napus
. Genetics, 2009, 182: 851-861
[本文引用: 1]
[27]Qiu D, Morgan C, Shi J, Long Y, Liu J, Li R, Zhuang X, Wang Y, Tan X, Dietrich E, Weihmann T.A comparative linkage map of oilseed rape and its use for QTL analysis of seed oil and erucic acid content
.Theor Appl Genet, 2006, 114: 67-80
[本文引用: 1]
相关话题/基因 序列 结构 设计 信息