基于简化基因组开发青冈和滇青冈微卫星引物

doi:10.12302/j.issn.1000-2006.202212002

国家林草科技领军期刊
中国精品科技期刊
中国高校百佳科技期刊
江苏省新闻出版政府奖期刊奖
RCCSE林学权威期刊（A+）
CSCD核心期刊
Scopus数据库收录期刊
中文核心期刊
SCD核心期刊

作者加群：102861116

微信公众号：南京林业大学学报

高级检索

南京林业大学学报（自然科学版） ›› 2024, Vol. 48 ›› Issue (6): 62-70.doi: 10.12302/j.issn.1000-2006.202212002

基于简化基因组开发青冈和滇青冈微卫星引物

欧阳泽怡¹^,²(), 李志辉¹, 牟虹霖¹, 姜小龙¹, 程勇³, 吴际友³^,^*()

1.中南林业科技大学林学院,湖南长沙 410004
2.湖南省植物园,湖南长沙 410116
3.湖南省林业科学院,湖南长沙 410004

收稿日期:2022-12-01 修回日期:2023-01-10 出版日期:2024-11-30 发布日期:2024-12-10
通讯作者: *吴际友(805610788@qq.com),研究员。
作者简介:
欧阳泽怡(455281453@qq.com),博士,工程师。
基金资助:
国家科技支撑计划(2012BAD21B03);湖南省自然科学基金青年基金项目(2021JJ41069);湖南省林业局林业科技攻关项目(XLK201706)

Designing of microsatellite primers for Quercus glauca and Q. glaucoides (Fagaceae) based on RAD-seq data

OUYANG Zeyi¹^,²(), LI Zhihui¹, MOU Honglin¹, JIANG Xiaolong¹, CHENG Yong³, WU Jiyou³^,^*()

1. College of Forestry, Central South University of Forestry and Technology, Changsha 410116, China
2. Hunan Botanical Garden, Changsha 410116, China
3. Hunan Academy of Forestry, Changsha 410004, China

Received:2022-12-01 Revised:2023-01-10 Online:2024-11-30 Published:2024-12-10

摘要/Abstract

摘要：

【目的】青冈(Quercus glauca)和滇青冈(Q. glaucoides)为珍贵树种,是东亚亚热带常绿阔叶林的重要建群种以及典型地理替代种,具有很高的生态和经济价值。开发其微卫星(simple sequence repeats,SSR)引物有助于分析其群体遗传格局及遗传多样性,为2种青冈林的管理和资源开发利用提供参考,也可为跨物种间的微卫星标记开发提供借鉴。【方法】分别基于3株青冈和3株滇青冈植株的简化基因组测序数据开发SSR引物。先后采用Stacks 2.0b软件process_radtags模型和SciRoKo 3.4软件对测序数据进行过滤和提取。利用Primer premier 6.0软件设计SSR引物。【结果】使用pyRAD 3.0.66软件对序列进行聚类,共鉴定出217个SSR位点,其中35%(76个)的SSR位点在青冈和滇青冈中均具有多态性。设计并筛选出28对SSR引物,分别在两个青冈群体和两个滇青冈群体(共48个个体)中进行巢式聚合酶链式反应(PCR)扩增。开发的28对SSR引物分布在青冈10条染色体上,SSR引物均在青冈和滇青冈个体中成功扩增,扩增率达到了90.7%。SSR分型分析结果表明,从遗传多样性方面看,所开发的28对引物中,共有176个等位基因;引物的等位基因数为3~13,平均为6.29个;期望杂合度为0.223~0.886,观察杂合度为 0.159~0.830。【结论】利用简化基因组数据可以快速、高效、经济地开发青冈和滇青冈通用微卫星引物,为后续2种青冈群体遗传学研究提供基础。也证明利用简化基因组测序数据,可开发近缘物种的通用微卫星引物。

关键词: 青冈, 滇青冈, 微卫星, 简化基因组

Abstract:

【Objective】Quercus glauca and Q. glaucoides are valuable and dominant species in the subtropical evergreen broadleaf forests of East Asia. They represent typical geographical vicarious species with significant ecological and economic importance. Therefore, the development of SSR primers for these two species can facilitate the analysis of genetic patterns and genetic diversity for the management and resource development of evergreen broadleaf forests as well as provide a reference for the development of microsatellite markers across species.【Method】This study developed SSR primers based on RADseq data from three Q. glauca and three Q. glaucoides individuals, respectively. The sequencing data were filtered and extracted using the process_radtags model in Stacks 2.0b software and SciRoKo 3.4 software sequentially. SSR primers were designed using Primer premier V6.0 software.【Result】The sequences were clustered using pyRAD 3.0.66, identifying a total of 217 SSR loci, 35% (76) of which were polymorphic in both species. Twenty-eight SSR primer pairs were designed and validated in two Q. glauca populations and two Q. glaucoides populations (48 individuals in total) through nested polymerase chain reaction (PCR) amplification. The 28 SSR primer pairs are distributed across ten chromosomes of the Q. glauca genome and successfully amplified in Q. glauca and Q. glaucoides individuals, with an amplification rate of 90.7%. The SSR genotyping analysis detected a total of 176 alleles, with the number of alleles per primer ranging from 3 to 13, and an average of 6.29. The expected and observed heterozygosity of the primers ranged from 0.223 to 0.886 and from 0.159 to 0.830, respectively.【Conclusion】The universal microsatellite primers developed in this study using RADseq data for Q. glauca and Q. glaucoides provide a basis for further population genetics studies of these species. In addition, this study demonstrates that RADseq data can be employed to rapidly, efficiently, and be used to cost-effectively develop universal microsatellite primers for closely related species.

Key words: Quercus glauca, Quercus glaucoides, microsatellite, RAD-seq

中图分类号:

S795

欧阳泽怡,李志辉,牟虹霖,等. 基于简化基因组开发青冈和滇青冈微卫星引物[J]. 南京林业大学学报（自然科学版）, 2024, 48(6): 62-70.

OUYANG Zeyi, LI Zhihui, MOU Honglin, JIANG Xiaolong, CHENG Yong, WU Jiyou. Designing of microsatellite primers for Quercus glauca and Q. glaucoides (Fagaceae) based on RAD-seq data[J].Journal of Nanjing Forestry University (Natural Science Edition), 2024, 48(6): 62-70.DOI: 10.12302/j.issn.1000-2006.202212002.

图/表 9

表1

表2

表3

图1

图2

图3

图4

表4

图5

参考文献 25

[1]	姜小龙. 福建青冈和岭南青冈系统发育关系及居群遗传结构[D]. 长沙: 中南林业科技大学, 2020.
	JIANG X L. Phylogenetic relationship and population genetic structure of Quercus chungii and Q.championii[D]. Changsha: Central South University of Forestry & Technology, 2020.DOI:10.27662/d.cnki.gznlc.2020.000507.
[2]	DENK T, GRIMM G W, MANOS P S, et al. An updated infrageneric classification of the oaks:review of previous taxonomic schemes and synthesis of evolutionary patterns in Oaks Physiological Ecology.Exploring the Functional Diversity of Genus Quercus L.[M]. Cham:Springer, 2017,13-38. DOI: 10.1007/978-3-319-69099-5.
[3]	MANOS P S, DOYLE J J, NIXON K C. Phylogeny,biogeography,and processes of molecular differentiation in Quercus Subgenus Quercus (Fagaceae)[J]. Mol Phylogenet Evol, 1999, 12(3):333-349.DOI: 10.1006/mpev.1999.0614.
[4]	罗艳, 周浙昆. 青冈亚属植物的地理分布[J]. 云南植物研究, 2001, 23(1):1-16,28.
	LUO Y, ZHOU Z K. Phytogeography of Quercus subg. Cyclobalanopsis[J]. Acta Bot Yunnanica, 2001, 23(1):1-16,28.
[5]	中国科学院中国植物志编辑委员会. 中国植物志-第四十二卷,第一分册[M]. 北京: 科学出版社,1993.
[6]	郭双兴. 云南临沧晚中新世邦卖组植物群[J]. 古生物学报, 2011, 50(3):353-408.
	GUO S X. The late Miocene Bangmai flora from Lincang County of Yunnan,southwestern China[J]. Acta Palaeontol Sin, 2011, 50(3):353-408.DOI: 10.19800/j.cnki.aps.2011.03.008.
[7]	郭双兴. 四川西部高原上新世植物群[J]. 古生物学报, 1978, 17(3):343-350,373.
	GUO S X. Pliocene floras of western Sichuan[J]. Acta Palaeontol Sin, 1978, 17(3):343-350,373.DOI: 10.19800/j.cnki.aps.1978.03.007.
[8]	罗冉, 吴委林, 张旸, 等. SSR分子标记在作物遗传育种中的应用[J]. 基因组学与应用生物学, 2010, 29(1):137-143.
	LUO R, WU W L, ZHANG Y, et al. SSR marker and its application to crop genetics and breeding[J]. Genom Appl Biol, 2010, 29(1):137-143.DOI: 10.3969/gab.029.000137.
[9]	SELKOE K A, TOONEN R J. Microsatellites for ecologists:a practical guide to using and evaluating microsatellite markers[J]. Ecol Lett, 2006, 9(5):615-629.DOI: 10.1111/j.1461-0248.2006.00889.x.
[10]	DAVEY J W, HOHENLOHE P A, ETTER P D, et al. Genome-wide genetic marker discovery and genotyping using next-generation sequencing[J]. Nat Rev Genet, 2011, 12(7):499-510.DOI: 10.1038/nrg3012.
[11]	胡亚亚, 刘兰服, 冀红柳, 等. 简化基因组测序技术研究进展[J]. 江苏师范大学学报(自然科学版), 2018, 36(4):63-68.
	HU Y Y, LIU L F, JI H L, et al. Research progress on the reduced-representation genome sequencing technique[J]. J Jiangsu Norm Univ (Nat Sci Ed), 2018, 36(4):63-68.DOI: 10.3969/j.issn.2095-4298.2018.04.012.
[12]	宁馨, 姜小龙, 邓敏, 等. 基于简化基因组数据开发岭南青冈微卫星引物[J]. 植物研究, 2020, 40(4):629-634.
	NING X, JIANG X L, DENG M, et al. Development of microsatellite primers of Quercus championii with RAD-seq data[J]. Bull Bot Res, 2020, 40(4):629-634.DOI: 10.7525/j.issn.1673-5102.2020.04.018.
[13]	李慧, 刘东超, 徐瑞瑞, 等. 基于RAD-seq技术的金银花SSR标记开发及鉴定[J]. 北京林业大学学报, 2021, 43(6):108-117.
	LI H, LIU D C, XU R R, et al. Development and identification of SSR markers based on RAD-seq of Lonicera japonica[J]. J Beijing For Univ, 2021, 43(6):108-117.DOI: 10.12171/j.1000-1522.20200337.
[14]	GAO Y, YIN S, LIU C, et al. A rapid approach for SSR development in Amorphophallus paeoniifolius using RAD-seq[J]. Taiwania, 2018, (63): 281-285.
[15]	王久利, 陈世龙, 邢睿, 等. 椭圆叶花锚简化基因组的SSR信息分析及SSR引物开发[J]. 植物研究, 2018, 38(2):292-297.
	WANG J L, CHEN S L, XING R, et al. Simplified genome SSR information and development of SSR primers of Halenia ellipitica (Gentianaceae)[J]. Bull Bot Res, 2018, 38(2):292-297.DOI: 10.7525/j.issn.1673-5102.2018.02.018.
[16]	CATCHEN J M, AMORES A, HOHENLOHE P, et al. Stacks:building and genotyping loci De novo from short-read sequences[J]. G3 Genes\|genomes\|genetics, 2011, 1(3):171-182.DOI: 10.1534/g3.111.000240.
[17]	KOFLER R, SCHLÖTTERER C, LELLEY T. SciRoKo:a new tool for whole genome microsatellite search and investigation[J]. Bioinformatics, 2007, 23(13):1683-1685.DOI: 10.1093/bioinformatics/btm157.
[18]	EATON D A R. PyRAD:Assembly of de novo RADseq loci for phylogenetic analyses[J]. Bioinformatics, 2014, 30(13):1844-1849.DOI: 10.1093/bioinformatics/btu121.
[19]	SINGH V K, MANGALAM A K, DWIVEDI S, et al. Primer premier:program for design of degenerate primers from a protein sequence[J]. BioTechniques, 1998, 24(2):318-319.DOI: 10.2144/98242pf02.
[20]	SCHUELKE M. An economic method for the fluorescent labeling of PCR fragments[J]. Nat Biotechnol, 2000, 18(2):233-234.DOI: 10.1038/72708.
[21]	HOLLAND M M, PARSON W. GeneMarker© HID:a reliable software tool for the analysis of forensic STR data[J]. J Forensic Sci, 2011, 56(1):29-35.DOI: 10.1111/j.1556-4029.2010.01565.x.
[22]	PEAKALL R, SMOUSE P E. Genalex 6:genetic analysis in Excel.Population genetic software for teaching and research[J]. Mol Ecol Notes, 2006, 6(1):288-295.DOI: 10.1111/j.1471-8286.2005.01155.x.
[23]	APARICIO J M, ORTEGO J, CORDERO P J. What should we weigh to estimate heterozygosity,alleles or loci?[J]. Mol Ecol, 2006, 15(14):4659-4665.DOI: 10.1111/j.1365-294X.2006.03111.x.
[24]	兰进茂, 覃瑞, 夏婧. 华蟹甲的简化基因组测序及SSR引物开发[J]. 中南民族大学学报(自然科学版), 2021, 40(3):258-263.
	LAN J M, QIN R, XIA J. Simplified genome SSR information and the development of SSR primers in Sinacalia tangutica (Asteraceae)[J]. J South Central Univ Natl (Nat Sci Ed), 2021, 40(3):258-263.DOI: 10.12130/znmdzk.20210307.
[25]	LV S Z, CHENG S, WANG Z Y, et al. Draft genome of the famous ornamental plant Paeonia suffruticosa[J]. Ecol Evol, 2020, 10(11):4518-4530.DOI: 10.1002/ece3.5965.

类型 type	样本编号 sample ID	采集地 location of collection	东经/(°) longitude (E)	北纬/(°) latitude (N)	海拔/m elevation
	gDM16207_5	浙江省遂昌县	118.85	28.38	718
	gDM18375	江西省吉安市	114.20	26.52	406
简化基因组测序RADseq	gDM18581	西藏自治区察隅县	97.08	28.56	1 645
	sDM11817	云南省晋宁县	102.77	24.70	2 004
	sDM17689	云南省巧家县	102.92	27.28	2 189
	sDM17825	云南省广南县	104.97	23.77	1 715
	gDM20391-gDM20402	江西省铅山县	117.75	27.89	658
	gDM21167-gDM21178	日本埼玉县	139.42	35.98	35
扩增验证PCR	sDM12688-sDM12693 sDM13035-sDM13040	云南省元谋县	101.95	25.72	2 093
	sDM10272-sDM 10275 sDM10280-sDM10287	广西壮族自治区隆林县	105.19	24.65	1 497

样本编号 sample ID	序列数 number of sequences	过滤后序列数 number of filter sequences	SSR聚类位点数 number of SSR Clustering Loci	一致性位点数 number of consistent loci	杂合度 heterozygosity	测序错误率/% sequencing error rate
gDM16207_5	9 391 215	8 073 849	20 415	1 808	0.006 656	0.070 44
gDM18375	9 095 105	8 382 814	23 195	1 808	0.007 887	0.080 43
gDM18581	9 140 314	8 397 947	26 498	1 808	0.006 120	0.067 51
sDM11817	9 515 559	8 166 681	28 418	1 808	0.007 409	0.069 65
sDM17689	9 364 847	8 376 899	30 229	1 808	0.006 690	0.079 68
sDM17825	8 466 497	7 689 333	33 170	1 808	0.007 648	0.075 63

引物号 primer code	正向序列(5'—3') forward sequence	反向引物序列(5'—3') forward sequence	NCBI序列号 accession number
GS2	TTGATGAACTGCGTACTCC	AAGCGAACCCAATTCCAC	MZ170721
GS3	ATGATTTACCGGATCTGGAG	CATGATTCAGTGCCATCTG	MZ170722
GS6	CTTGTGCCGACATTCTTG	ACCACCAGCCTGAAATTC	MZ170723
GS15	ATTTCTCTTCTGGGTCGGA	GCTGACCTTTCTTGACTTTG	MZ170724
GS19	GGTGAGTAGACGGAGGTT	CCTTCACTTCTTCTTCATCCT	MZ170725
GS28	CACCAAGTGTGGCCTTTC	TGAAGAGCAGAAATGGATGG	MZ170726
GS32	TCTGGTATAGAATCCAACGA	GCCCTTCATACTTCTTTGTT	MZ170727
GS33	GGTTCTGGTTCTGGTTCTG	CGCTACAAATTCCTTGCTTT	MZ170728
GS35	CTCATCTCCTACCTCATTCTT	TCCCTTTGCTGACTTTGT	MZ170729
GS74	TTCCTACTTCCTAGTGAGAG	ATGAGGTCAGGCTCCTTA	MZ170730
GS83	CACAATCACAATCACCAACA	TCCTTGTCTTCCTCCTCTT	MZ170731
GS84	TGGCTTCTTCTTCCTATATGAG	GCGTATGCTGGTAAGTGAT	MZ170732
GS88	TCGGCACGGATCGTAATG	CAACCCATGAAGCTGAACT	MZ170733
GS105	CGCTTTCGCTTTGAAATTG	CTACACAAATTGCGTTTCAG	MZ170734
GS117	AGGATCTGGATCAGTAAGC	AACGTAGCACAGTAACATG	MZ170735
GS119	ATCCGGTCATGTACTTCATT	GCAACCAAACGGAGAGTA	MZ170736
GS122	CGCCATTTGGTTGTAATCT	GCTGTTAGATGATCCGTTAG	MZ170737
GS150	GGGATGAAACAAATTAGGAC	CCAATCGTTATTTAGCAAGG	MZ170738
GS160	CCTGTCTTCTTCTTCTTCTTC	ACTAGGAGATGATGATATACCC	MZ170739
GS165	CGGGTTTGGTCCAACAAC	CGCACTGAGCTTCCCTTT	MZ170740
GS169	GGCTGGCTTCCAAATCTT	AATATGGAATGGACCACTCA	MZ170741
GS171	TTCCTCTAGTCTCTTCAGTCT	TGAGATACATACTTGGCTTACC	MZ170742
GS172	CTCTCAGCTACAACAACAA	AGACCCATTCTCACTTCAA	MZ170743
GS180	CCAAAGAACTGGAGCTTTATG	CTTTCACCATGTGTCCCTC	MZ170744
GS185	TAGGTGGCGGTCTAGGTT	AGCTTTGCACAGAGTCTC	MZ170745
GS186	TCTGTGTTATGGACTCTGAG	ACAACGTAAGACGGTAAGA	MZ170746
GS209	TTCAGTGAAATATCACCTGC	CCAAGAAATCATCCTCTTCTC	MZ170747
GS216	AAGTGAGACTCGGTTCAA	GGCTTCCCTCTAAGTTTCT	MZ170721

引物号 primer code	重复基序 repeat motifs	序列长度/bp sequences length	个体数 sample size	N_Aa	N_Ab	期望杂合度 expected heterozygosity	观测杂合度 observed heterozygosity
GS2	TGT	138~156	47	5	7	0.739	0.319
GS3	GTT	137~155	48	2	5	0.471	0.333
GS6	CT	134~160	42	4	13	0.862	0.714
GS15	CCA	125~134	44	3	4	0.553	0.341
GS19	TGGTTT	114~126	47	3	3	0.594	0.617
GS28	CT	133~153	46	5	11	0.841	0.783
GS32	AACAAA	109~133	39	2	5	0.642	0.410
GS33	GGTTCT	121~133	48	3	3	0.484	0.750
GS35	AG	121~143	47	4	12	0.873	0.830
GS74	GA	124~128	33	3	3	0.432	0.212
GS83	AG	108~130	46	5	12	0.886	0.739
GS84	CA	129~133	44	2	3	0.223	0.159
GS88	TGT	121~127	40	3	3	0.405	0.250
GS105	AG	115~129	41	5	8	0.808	0.463
GS117	AACAC	121~141	33	2	5	0.730	0.636
GS119	TTTGC	117~127	43	2	3	0.577	0.326
GS122	AG	130~136	42	2	4	0.584	0.500
GS150	ATAAA	139~154	46	2	4	0.664	0.283
GS160	TTC	119~138	43	4	8	0.810	0.465
GS165	GAA	120~129	40	3	4	0.669	0.475
GS169	GTGA	145~153	46	3	3	0.591	0.326
GS171	CT	105~121	45	3	10	0.836	0.733
GS172	AAG	134~152	34	3	7	0.755	0.382
GS180	AGA	108~132	43	5	9	0.740	0.442
GS185	TC	146~148	48	2	6	0.736	0.396
GS186	GA	146~168	48	5	9	0.795	0.792
GS209	GA	124~144	48	4	7	0.678	0.667
GS216	AGT	124~136	47	4	5	0.701	0.404

基于简化基因组开发青冈和滇青冈微卫星引物

Designing of microsatellite primers for Quercus glauca and Q. glaucoides (Fagaceae) based on RAD-seq data

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

图/表 9

参考文献 25

相关文章 15

编辑推荐

Metrics

本文评价

作者

读者

审稿

[1]	罗芊芊, 李峰卿, 肖德卿, 邓章文, 王建华, 周志春. 两个南方红豆杉天然居群的交配系统分析[J]. 南京林业大学学报（自然科学版）, 2023, 47(5): 80-86.
[2]	邓平, 赵英, 王霞, 陈秋佑, 吴敏. 水杨酸对NaHCO₃胁迫下桂西北喀斯特地区青冈栎种子萌发的影响[J]. 南京林业大学学报（自然科学版）, 2021, 45(4): 114-122.
[3]	李峰卿, 王秀花, 楚秀丽, 张东北, 吴小林, 周生财, 叶明. 缓释肥N/P比及加载量对5种珍贵树种1年生苗生长和养分库构建的影响[J]. 南京林业大学学报（自然科学版）, 2020, 44(1): 72-80.
[4]	李伟强, 戴晓港, 李小平, 李淑娴. 美洲黑杨种质材料倍性鉴定[J]. 南京林业大学学报（自然科学版）, 2019, 43(5): 51-58.
[5]	张晓晨,赵洋,熊中人,伊贤贵,陈昕. 宝华山青冈种群年龄结构及点格局分析[J]. 南京林业大学学报（自然科学版）, 2018, 42(06): 77-83.
[6]	刘延惠,丁访军,舒德远,崔迎春,侯贻菊,赵文君. 茂兰喀斯特原生林细叶青冈树干液流环境响应特征[J]. 南京林业大学学报（自然科学版）, 2017, 41(03): 77-85.
[7]	马秋月,廖卓毅,张得芳,戴晓港,陈赢男,李淑娴. 碧桃花瓣转录组微卫星特征分析[J]. 南京林业大学学报（自然科学版）, 2015, 39(03): 34-38.
[8]	孔凡明,王小龙,陈赢男,李淑娴,徐进,尹佟明. 松属编码区微卫星特征和相应基因功能分析[J]. 南京林业大学学报（自然科学版）, 2014, 38(02): 47-51.
[9]	董维,苏舒,游录鹏,黄守程,戚金亮,陆桂华,黄应华,杨永华. 高粱F₆代群体分蘖数的QTL定位[J]. 南京林业大学学报（自然科学版）, 2013, 37(02): 55-58.
[10]	史洁,尹佟明,管宏伟,戴晓港,陈金慧,施季森. 油茶基因组微卫星特征分析[J]. 南京林业大学学报（自然科学版）, 2012, 36(02): 47-51.
[11]	刘菁菁，戴晓港，王洁，李淑娴，尹佟明. 杨树微卫星序列对基因表达频率的影响及表达序列中微卫星特征的分析[J]. 南京林业大学学报（自然科学版）, 2011, 35(01): 11-14.
[12]	栾启福，姜景民，张守攻. 松属GAs相关ESTSSR标记与火炬松×洪都拉斯加勒比松苗高相关性分析[J]. 南京林业大学学报（自然科学版）, 2011, 35(01): 6-10.
[13]	刘映良，薛建辉. 鼠类及种子特征对水青冈种子命运的影响[J]. 南京林业大学学报（自然科学版）, 2010, 34(05): 12-16.
[14]	曹建新,姜远标,张劲峰,苏文华,张光飞. 不同土壤条件下滇青冈和元江栲幼苗生长特征的比较[J]. 南京林业大学学报（自然科学版）, 2009, 33(01): 79-82.
[15]	杨礼旦,王安文,李朝志. 水青冈群落物种多样性及乔木种群分布格局[J]. 南京林业大学学报（自然科学版）, 2005, 29(03): 107-110.