基于甘薯耐盐转录组测序的SSR和SNP特征分析

doi:10.11923/j.issn.2095-4050.cjas2021-0124

摘要/Abstract

摘要：

旨在对甘薯耐盐转录组中潜在的SSR和SNP位点进行挖掘及特征分析，以完善甘薯分子标记。试验以不同耐盐甘薯品种的转录组序列为基础，利用MISA、GATK软件分析了转录本中的SSR和SNP信息。结果显示，研究共获得SSR位点33192个，分布在24323条Unigenes中，出现频率为21.11%。其中，6271条Unigenes含有超过1个以上的SSR位点。SSR类型以单核苷酸重复SSR和双核苷酸重复SSR为主。在重复基元中，A/T和AG/CT所占的比例最高，为优势基元。在157252条Unigenes中共挖掘到SNP位点7691906个，分布密度为0.08个/bp。其中转换类型有4729922个，颠换类型有2961984个。在所有突变类型中，C/T和A/G含量最高，分别为占总数的32.34%和29.15%。综上所述，甘薯转录组具有较丰富的SSR和SNP标记，多态性较高，可为甘薯抗逆品种筛选、种质培育等研究奠定理论基础。

关键词: 甘薯, 盐胁迫, 转录组, SSR, SNP

Abstract:

The aims are to mine and analyze the potential SSR and SNP loci in sweet potato transcriptome under salt stress, and enrich the molecular markers of sweet potato. Based on the transcriptome sequences of different salt tolerant sweet potato varieties, the SSR and SNP were analyzed by MISA and GATK software. The results showed that, a total of 33192 SSR loci were identified in 24323 Unigenes with a frequency of 21.11%. 6271 Unigenes contained more than one SSR locus. Mono-nucleotide repeat SSR and di-nucleotide repeat SSR were the main types. The two repeat motifs, A/T and AG/CT, had the highest proportion and were the dominant motifs of sweet potato SSR. After searching SNP loci in 157252 Unigenes, 7691906 SNP loci were identified with a distribution density of 0.08 SNP per bp. The predicted SNPs included 4729922 transitions and 2961984 transversions. C/T and A/G were the main types, accounting for 32.34% and 29.15% of all the SNPs. Taken together, the transcriptome of sweet potato has abundant SSR and SNP markers with high polymorphism, and the study could lay a theoretical basis for the selection of stress resistant varieties and germplasm cultivation of sweet potato.

Key words: sweet potato, salt stress, transcriptome, SSR, SNP

中图分类号:

S531

张小红, 彭琼, 鄢铮. 基于甘薯耐盐转录组测序的SSR和SNP特征分析[J]. 农学学报, 2022, 12(11): 55-61.

ZHANG Xiaohong, PENG Qiong, YAN Zheng. SSR and SNP in Salt Tolerant Sweet Potato: Characteristics Analysis Based on Transcriptome Sequencing[J]. Journal of Agriculture, 2022, 12(11): 55-61.

图/表 6

参考文献 43

[1]	BRUCKNER A W. AFLP-based genetic diversity assessment of global sweetpotato (Ipomoea batatas (L.) Lam.) germplasm resources:Progress toward the development of a sweetpotato core collection[D]. North Carolina: North Carolina State University, 2004.
[2]	聂世现. 燃料乙醇专用型甘薯品种筛选及乙醇发酵技术研究[D]. 合肥: 安徽大学, 2010.
[3]	陆建珍, 汪翔, 秦建军, 等. 中国甘薯种植业时空布局分析及产业发展建议[J]. 天津农业科学, 2020, 26(3):53-62.
[4]	陆国权. 中国甘薯育成品种系谱[J]. 中国甘薯, 1990(4):26-28.
[5]	李强, 刘庆昌, 马代夫, 等. 中国甘薯主要育成品种的遗传多样性及遗传趋势[J]. 江苏农业学报, 2009, 25(2):253-259.
[6]	赵冬兰, 郑立涛, 唐君, 等. 甘薯种质资源遗传稳定性及遗传多样性SSR分析[J]. 植物遗传资源学报, 2011, 12(3):389-395.
[7]	张林甦, 刘承贵, 蒙慧彤, 等. 基于蛇足石杉转录组的SNP和SSR多态性分析[J]. 黔南民族医专学报, 2018, 31(4):235-238.
[8]	周军永, 陆丽娟, 刘茂, 等. 基于李府贡枣转录组测序的SSR和SNP特征分析[J]. 江苏农业科学, 2019, 47(4):51-54.
[9]	严佳文, 肖图舰, 袁启凤, 等. 基于转录组序列的火龙果SSR和SNP多态性分析[J]. 热带作物学报, 2018, 39(7):85-90.
[10]	许家磊. 基于甘薯徐781和徐薯18转录组测序的SNP标记开发[D]. 北京: 中国农业科学院, 2015.
[11]	JIANG D, YE Q L, WANG F S, et al. The mining of citrus EST-SNP and its application in cultivar discrimination[J]. Agricultural sciences in China, 2010, 9(2):179-190. doi: 10.1016/S1671-2927(09)60082-1 URL
[12]	陆丹, 牛楠, 李玥莹. SSR标记技术在植物基因组研究上的应用[J]. 沈阳师范大学学报(自然科学版), 2010(1):83-85.
[13]	MCCOUCH S R, ZHAO K, WRIGHT M, et al. Development of genome-wide SNP assays for rice[J]. Breeding science, 2010, 60(5):524-535. doi: 10.1270/jsbbs.60.524 URL
[14]	肖炳光, 邱杰, 曹培健, 等. 利用基因组简约法开发烟草SNP标记及遗传作图[J]. 作物学报, 2014, 40:397-404.
[15]	何虎翼, 谭冠宁, 何新民, 等. 几种基于PCR的分子标记在甘薯遗传多样性研究中的应用[J]. 生物技术进展, 2014, 4(4):245-250.
[16]	张明飞, 于卓, 于肖夏, 等. 四倍体马铃薯SRAP分子遗传连锁图谱的构建[J]. 草业学报, 2019, 28(9):190-199.
[17]	WANG Z, LI J, LUO Z, et al. Characterization and development of EST-derived SSR markers in cultivated sweetpotato (Ipomoea batatas)[J]. BMC plant biology, 2011, 11(1):139-139. doi: 10.1186/1471-2229-11-139 URL
[18]	张超凡, 黄艳岚, 周虹, 等. 湖南甘薯品种遗传多样性的SSR分析[J]. 湖南农业科学, 2011(21):1-2.
[19]	XIE F L, BURKLEW C E, YANG Y F, et al. De novo sequencing and a comprehensive analysis of purple sweet potato (Impomoea batatas L.) transcriptome[J]. Planta, 2012, 236(1):101-113. doi: 10.1007/s00425-012-1591-4 URL
[20]	ZHAO H Y, ZHANG S S, WANG F B, et al. Comparative transcriptome analysis of purple-fleshed sweet potato provides insights into the molecular mechanism of anthocyanin biosynthesis[J]. Frontiers of agricultural science & engineering, 2018, 5(2):214-225.
[21]	GRABHERR M G, HAAS B J, YASSOUR M, et al. Trinity: reconstructing a full-length transcriptome without a genome from RNA-Seq data[J]. Nature biotechnology, 2011, 29(7):644-652. doi: 10.1038/nbt.1883 URL
[22]	UNTERGASSER A, CUTCUTACHE I, KORESSAAR T, et al. Primer3--new capabilities and interfaces[J]. Nucleic acids research, 2012, 40(15):e115. doi: 10.1093/nar/gks596 URL
[23]	SCHULER G D. Sequence mapping by electronic PCR[J]. Genome research, 1997, 7(5):541-550. pmid: 9149949
[24]	DOBIN A, DAVIS C A, SCHLESINGER F, et al. STAR: ultrafast universal RNA-seq aligner[J]. Bioinformatics, 2013, 29(1):15-21. doi: 10.1093/bioinformatics/bts635 pmid: 23104886
[25]	MCKENNA A, HANNA M, BANKS E, et al. The genome analysis toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data[J]. Genome research, 2010, 20:1297-1303. doi: 10.1101/gr.107524.110 pmid: 20644199
[26]	魏明明, 陈钰辉, 刘富中, 等. 基于转录组测序的茄子SSR标记开发[J]. 植物遗传资源学报, 2016, l7(6):1082-1091.
[27]	杨新笋. 基于SSR, SNP和形态学标记的甘薯种质资源遗传多样性研究[D]. 北京: 中国农业大学, 2016.
[28]	王学文. SSR和SNP遗传分子标记与作物分子育种[C]//. 遗传学进步促进粮食安全与人口健康高峰论坛论文集. 郑州: 中国遗传学会, 2012:20.
[29]	杨帅, 杨喜翠, 何方, 等. 贵紫麦1号转录组EST-SSR标记发掘与应用[J]. 河南农业科学, 2020, 49(5):15-30.
[30]	刘强, 赵丽, 郭虹霞, 等. 玉米野生近缘种大刍草幼苗叶片的转录组SSR分析[J]. 玉米科学, 2020, 28(2):73-79.
[31]	王洋洋, 单文琪, 徐文龙, 等. 印度南瓜转录组SSR信息分析及其多态性研究[J]. 园艺学报, 2016, 43(3):578-586.
[32]	刘峰, 王运生, 田雪亮, 等. 辣椒转录组SSR挖掘及其多态性分析[J]. 园艺学报, 2012, 39(1):168-174.
[33]	WANG Z Y, FANG B P, CHEN J Y, et al. De novo assembly and characterization of root transcriptome using Illumina paired-end sequencing and development of cSSR markers in sweet potato (Ipomoea batatas)[J]. BMC genomics, 2010, 11(1):726. doi: 10.1186/1471-2164-11-726 URL
[34]	LI R, ZHAI H, CHEN K, et al. De Novo transcriptome sequencing of the orange-fleshed sweet potato and analysis of differentially expressed genes related to carotenoid biosynthesis[J]. International journal of genomics, 2015, 2015:1-10.
[35]	ZHU H, ZHOU Y Y, ZHAI H, et al. Transcriptome profiling reveals insights into the molecular mechanism of drought tolerance in sweetpotato[J]. Journal of integrative agriculture, 2019, 18(1):13-27.
[36]	林艳. 木荷基因组SSR与SNP位点挖掘及种质资源遗传多样性分析[D]. 广州: 华南农业大学, 2018.
[37]	郑燕, 张耿, 吴为人. 禾本科植物微卫星序列的特征分析和比较[J]. 基因组学与应用生物学, 2011, 30(5):9-16.
[38]	蒋超, 袁媛, 刘贵明, 等. 基于EST-SSR的金银花分子鉴别方法研究[J]. 药学学报, 2012, 47(6):803-810.
[39]	魏秀清, 许玲, 章希娟, 等. 莲雾转录组SSR信息分析及其分子标记开发[J]. 园艺学报, 2018, 45(3):541-551.
[40]	朱林. 基于140种植物全基因组序列的SSRs比较分析与长SSRs潜在功能的初探[D]. 泰安: 山东农业大学, 2020.
[41]	TEMNYKH S, DECLERCK G, LUKASHOVA A, et al. Computational and experimental analysis of microsatellites in rice (Oryza sativa L.):frequency, length variation, transposon associations, and genetic marker potential[J]. Genome research, 2001, 11(8):1441-1452. doi: 10.1101/gr.184001 URL
[42]	罗怀容, 施鹏, 张亚平. 单核苷酸多态性的研究技术[J]. 遗传, 2001, 23(5):471-476.
[43]	BROOKES A J. The essence of SNPs[J]. Gene, 1999, 234(2):177-186. doi: 10.1016/s0378-1119(99)00219-x pmid: 10395891

项目	数量
评估的序列数目	157252
评估的序列总碱基数 (bp)	90649057
识别的SSR总数	33192
包含SSR的序列数目	24323
包含1个以上SSR的序列数目	6271
以复合物形式存在的SSR数目	4471
单核苷酸重复SSR数目	18718
双核苷酸重复SSR数目	8121
三核苷酸重复SSR数目	5565
四核苷酸重复SSR数目	601
五核苷酸重复SSR数目	129
六核苷酸重复SSR数目	58

项目	数量
评估的序列数目	157252
评估的序列总碱基数 (bp)	90649057
识别的SSR总数	33192
包含SSR的序列数目	24323
包含1个以上SSR的序列数目	6271
以复合物形式存在的SSR数目	4471
单核苷酸重复SSR数目	18718
双核苷酸重复SSR数目	8121
三核苷酸重复SSR数目	5565
四核苷酸重复SSR数目	601
五核苷酸重复SSR数目	129
六核苷酸重复SSR数目	58

重复类型	重复基元类型总数	重复基元	数量	占比/%
单核苷酸重复	2	A/T	18391	55.41
单核苷酸重复	2	C/G	327	0.99
双核苷酸重复	4	AG/CT	3706	11.17
		AT/AT	3185	9.60
		AC/GT	1220	3.68
		CG/CG	10	0.03
三核苷酸重复	10	AAT/ATT	1606	4.84
		AAG/CTT	1325	3.99
		ATC/ATG	603	1.82
		CCG/CGG	585	1.76
		AAC/GTT	459	1.38
		剩余三核苷酸重复基元	987	2.97
四核苷酸重复	28	AAAT/ATTT	231	0.70
		AAAG/CTTT	89	0.27
		ACAT/ATGT	65	0.20
		AAAC/GTTT	56	0.17
		剩余四核苷酸重复基元	160	0.48
五核苷酸重复	40	AAAAT/ATTTT	31	0.09
		AAAAC/GTTTT	28	0.08
		AAAAG/CTTTT	9	0.03
		AAGAG/CTCTT	9	0.03
		剩余五核苷酸重复基元	52	0.16
六核苷酸重复	36	AAAAAC/GTTTTT	6	0.02
		AAAAAT/ATTTTT	6	0.02
		AAAGAG/CTCTTT	3	0.01
		AAAGGC/CCTTTG	2	0.01
		剩余六核苷酸重复基元	41	0.12

重复类型	重复基元类型总数	重复基元	数量	占比/%
单核苷酸重复	2	A/T	18391	55.41
单核苷酸重复	2	C/G	327	0.99
双核苷酸重复	4	AG/CT	3706	11.17
		AT/AT	3185	9.60
		AC/GT	1220	3.68
		CG/CG	10	0.03
三核苷酸重复	10	AAT/ATT	1606	4.84
		AAG/CTT	1325	3.99
		ATC/ATG	603	1.82
		CCG/CGG	585	1.76
		AAC/GTT	459	1.38
		剩余三核苷酸重复基元	987	2.97
四核苷酸重复	28	AAAT/ATTT	231	0.70
		AAAG/CTTT	89	0.27
		ACAT/ATGT	65	0.20
		AAAC/GTTT	56	0.17
		剩余四核苷酸重复基元	160	0.48
五核苷酸重复	40	AAAAT/ATTTT	31	0.09
		AAAAC/GTTTT	28	0.08
		AAAAG/CTTTT	9	0.03
		AAGAG/CTCTT	9	0.03
		剩余五核苷酸重复基元	52	0.16
六核苷酸重复	36	AAAAAC/GTTTTT	6	0.02
		AAAAAT/ATTTTT	6	0.02
		AAAGAG/CTCTTT	3	0.01
		AAAGGC/CCTTTG	2	0.01
		剩余六核苷酸重复基元	41	0.12

重复类型	重复次数									总计
重复类型	5	6	7	8	9	10	11	12	>12	总计
单核苷酸	0	0	0	0	0	5667	3174	2122	7755	18718
双核苷酸	0	2679	1579	1104	711	512	536	590	410	8121
三核苷酸	2930	1283	591	265	59	71	58	52	256	5565
四核苷酸	415	146	12	11	4	6	2	3	2	601
五核苷酸	99	16	9	4	0	0	1	0	0	129
六核苷酸	26	15	10	4	1	0	1	0	1	58
总计	3470	4139	2201	1388	775	6256	3772	2767	8424	33192