Genetic diversity and molecular genealogy of local silkworm varieties

doi:10.3724/SP.J.1003.2013.11035

Abstract

Abstract:

In order to explore the genetic diversity and systematic differentiation pattern among silkworm varieties, aiming to guide hybridization breeding, we sequenced a total of 72 Bmamy2gene fragments from local silkworm varieties. The analysis of nucleotide sequence diversity and systematic differentiation indicated that there was rich genovariation in the sequencing region of Bmamy2gene, and the base mutation rate is 5.6- 8.2%, the haplotype diversity is 0.8294, and the nucleotide diversity is 0.0236±0.00122, suggesting Bmamy2being a better marking gene with rich nucleotide sequence diversity, based on which the genetic diversity among different local silkworm varieties can be identified. The same heredity population structure is proclaimed by several analysis methods that every clade consisting of varieties from different geosystems and ecological types, while the varieties from the same geosystem and ecotype belong to different clades in the phylogeny. There is no population structure pattern that different varieties claded together according to geosystem or ecotype. It can be speculated that the silkworm origins from mixture of kinds of several voltinism mulberry silkworm, Bombyx mandarina, while the domestication events took place in several regions, from which the domesticated mulberry silkworms are all devoting to the domesticated silkworm population of today.

Key words: Bombyx mori, germplasm resources, genetic diversity, molecule genealogy

Zhouhe Du,Junfeng Liu,Binbin Liu,Yanchun Zuo,Jianmei Wu,Yi-an Chen,Jianfei Zhang,Cheng Lu. Genetic diversity and molecular genealogy of local silkworm varieties[J]. Biodiv Sci, 2013, 21(3): 315-325.

Figures/Tables 9

References 38

1	Arunkumar KP, Muralidhar M, Nagaraju J (2006) Molecular phylogeny of silkmoths reveals the origin of domesticated silkmoth, Bombyx mori from Chinese Bombyx mandarina and paternal inheritance of Antheraea proylei mitochondrial DNA. Molecular Phylogenetics and Evolution, 40, 419-427.
2	Bandelt HJ, Forster P, Röhl A (1999) Median-joining networks for inferring intraspecific phylogenies. Molecular Biology and Evolution, 16, 37-48.
3	Brown WM (1983) Evolution of animal mitochondrial DNA. In: Evolution of Genes and Proteins (eds Nei M, Koehned RK), pp. 62-88. Sinauer Associates, Sunderland MA.
4	Chen LY (陈丽媛), Zhao QL (赵巧玲), Shen XJ (沈兴家), Zhang ZF (张志芳), Tang SM (唐顺明), Xu AY (徐安英), Zhang GZ (张国政), Guo XJ (郭锡杰) (2007) Nucleotide sequences in the A+T rich regions of mitochondrial DNA from local races of Bombyx mori and their molecular evolution. Science of Sericulture(蚕业科学), 33, 5-13. (in Chinese with English abstract)
5	Cheng XY (成新跃), Zhou HZ (周红章), Zhang GX (张广学) (2000) Perspective of molecular biological techniques applied in insect systematics. Acta Zootaxonomica Sinica(动物分类学报), 25, 121-133. (in Chinese with English abstract)
6	Du ZH (杜周和) (2009) The Origin and Differentiation of Silkworm, Bombyx mori, Based on the Bmamy2 Gene (基于Bmamy2基因的家蚕起源与分化研究), PhD dissertation, Southwest University, Chongqing. (in Chinese)
7	Du ZH (杜周和), Dong ZP (董占鹏), Liu JF (刘俊凤), Liu BB (刘彬斌), Lu C (鲁成) (2008) Cloning, expression pattern and sequence analysis of Bmamy2 gene in the silkwom, Bombyx mori. Science of Sericulture(蚕业科学), 34, 642-649. (in Chinese with English abstract)
8	Du ZH (杜周和), Liu JF (刘俊凤), Liu BB (刘彬斌), Dong ZP (董占鹏), Yu QY (余泉友), Lu C (鲁成), Chen YA (陈义安) (2009) Genetic diversity of the wild mulberry silkworm, Bombyx mandarina, in China and its phylogenetic relationship with the domesticated silkworm, Bombyx mori, based on amy gene. Acta Entomologica Sinica(昆虫学报), 52, 1338-1348. (in Chinese with English abstract)
9	Excoffier L, Smouse PE, Quattor JM (1992) Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics, 131, 479-491.
10	Fang SM (房守敏), Yu QY (余泉友), Li B (李斌), Dai FY (代方银), Miao XX (苗雪霞), Huang YP (黄勇平), Lu C (鲁成) (2007) Comparison of SSR, AFLP and RAPD markers for genetic analysis of domestic silkworm. Scientia Agricultura Sinica(中国农业科学), 40, 2926-2930. (in Chinese with English abstract)
11	Fang SM (房守敏), Zhang L (张烈), Lu C (鲁成) (2010) Study and application in origin and evolution of mitochondrial DNA in the silkmoths. Chinese Bulletin of Entomology(昆虫知识), 47, 439-445. (in Chinese with English abstract)
12	Gyllensten U, Wharton D, Josefsson A, Wilson AC (1991) Paternal inheritance of mitochondrial DNA in mice. Nature, 352, 255-257.
13	Hu XL, Cao GL, Xue RY, Zheng XJ, Zhang X, Duan HR, Gong CL (2010) The complete mitogenome and phylogenetic analysis of Bombyx mandarina variety Qingzhou. Molecular Biology Report, 37, 2599-2608.
14	Jiang YL(蒋猷龙) (1982) The Origin of Silkworm, Bombyx mori, and Its Differentiation(家蚕的起源与分化). Jiangsu Science and Technology Press, Nanjing(in Chinese)
15	Jiang YL (蒋猷龙) (1989) Advance of the study on the origin of silkworm and its differentiation. Newsletter of Agriculture and Animal Husbandry(农牧情报研究), (7), 43-47. (in Chinese)
16	Kenji Y, Hideki S, Toshiki T, Eiichi K, Makoto K (2011) Nucleotide sequence variation in mitochondrial COI gene among 147 silkworm (Bombyx mori) varieties from Japanese, Chinese, European and moltinism classes. Genes & Genetic Systems, 86, 315-323.
17	Li AL, Zhao QL, Tang SM, Zhang ZF, Pan SY, Shen GF (2005) Molecular phylogeny of the domesticated silkworm, Bombyx mori, based on the sequences of mitochondrial cytochrome b genes. Journal of Genetics, 84, 137-142.
18	Li D, Guo YR, Shao HJ, Tellier LC, Wang J, Xiang ZH, Xia QY (2010) Genetic diversity, molecular phylogeny and selection evidence of the silkworm mitochondria implicated by complete resequencing of 41 genomes. BMC Evolutionary Biology, 10, 81.
19	Li HW (李焕文) (1985) Exploration of the origin and voltinism change of silkworm, Bombyx mori, and mulberry silkworm, Bombyx mandarina. Newsletter of Sericultural Science(蚕学通讯), (1), 41-45. (in Chinese)
20	Li MW (李木旺) (2005) The Application of Microsatellite Markers in Silkworm (Bombyx mori): Genetic Diversity, QTL Analysis, Gene Mapping and Marker Assisted Selection (微卫星遗传标记在家蚕中的应用研究——遗传多样性分析、QTL分析、基因定位与分子标记辅助育种). PhD dissertation, Institute of Plant Physiology & Ecology, Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, Shanghai. (in Chinese)
21	Lu C (鲁成), Yu HS (余红仕), Xiang ZH (向仲怀) (2002a) Molecular systematic studies on Chinese mandarina silkworm (Bombyx mandarina M.) and domestic silkworm (Bombyx mori L.). Scientia Agricultura Sinica(中国农业科学), 35, 94-101. (in Chinese with English abstract)
22	Lu C (鲁成), Zhao AC (赵爱春), Zhou ZY (周泽扬), Xiang ZH (向仲怀), Wan CL (万春玲) (2001) AFLP analysis of Chinese mulberry wild silkworm (Bombyx mandarina) and domestic silkworm (Bombyx mori). Acta Sericologica Sinica(蚕业科学), 27, 243-252. (in Chinese with English abstract)
23	Lu C (鲁成), Yu HS (余红仕), Xiang ZH (向仲怀) (2002b) The genetic diversity and phylogenetic relationship of Bombyx mandarina and Bombyx mori from China based on RAPD analysis. Acta Entomologica Sinica(昆虫学报), 45, 198-203. (in Chinese with English abstract)
24	Moritz C, Dowling TE, Brown WM (1987) Evolution of animal mitochondrial DNA relevance for population biology and systematic. Annual Review of Ecology and Systematics, 18, 269-289.
25	Rozas J, Sánchez-DelBarrio JC, Messeguer X, Roza R (2003) DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics, 19, 2496-2497.
26	Tamura K, Dudley J, Nei M, Kumar S (2007) MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Molecular Biology and Evolution, 24, 1596-1599.
27	Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG (1997) The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Research, 25, 4876-4882.
28	Wei S(卫斯) (1985) Primary investigation of the origin times of sericulture cultivation in China. In: Study of Agricultural History, Vol. 6(农史研究, 第六辑). China Agriculture Press, Beijing. (in Chinese)
29	Wei S (卫斯) (2005) The mark stela in our modern archeology history times: the course and importance of the first excavation in Xiyin Village historical site. .
30	Wu W (2006) Bioinformatics Analysis of the Comprise and Evolution of Animal Mitochondrial Genome (动物线粒体基因组组成与进化的生物信息学研究). PhD dissertation, Southwest University, Chongqing. (in Chinese)
31	Xia QY, Guo YR, Zhang Z, Li D, Xuan ZL, Li Z, Dai FY, Li YG, Cheng DJ, Li RQ, Cheng TC, Jiang T, Becquet C, Xu X, Liu C, Zha XF, Fan W, Lin Y, Shen YH, Jiang L, Jensen J, Hellmann I, Tang S, Zhao P, Xu HF, Yu C, Zhang GJ, Li J, Cao JJ, Liu SP, He NJ, Zhou Y, Liu H, Zhao J, Ye C, Du ZH, Pan GQ, Zhao AC, Shao HJ, Zeng W, Wu P, Li CF, Pan MH, Li JJ, Yin XY, Li DW, Wang J, Zheng HS, Wang W, Zhang XQ, Li SG, Yang HM, Lu C, Nielsen R, Zhou ZY, Wang J, Xiang ZH, Wang J (2009) Complete Resequencing of 40 genomes reveals domestication events and genes in silkworm (Bombyx). Science, 326, 433-436.
32	Xia QY (夏庆友), Zhou ZY (周泽扬), Lu C (鲁成), Xiang ZH (向仲怀) (1998) Molecular phylogenetic study on the racial differentiation of Bombyx mori by random amplified polymorphic DNA (RAPD) markers. Acta Entomologica Sinica(昆虫学报), 41, 32-40. (in Chinese with English abstract)
33	Xiang ZH(向仲怀) (1994) The Inheritance and Thremmatology of Silkworm, Bombyx mori(家蚕遗传育种学). China Agriculture Press, Beijing(in Chinese)
34	Yoshitake N (吉武成美), Jiang YL (蒋猷龙) (1987) The origin of silkworm, Bombyx mori, and its differentiation. Acta Sericologica Sinica(蚕业科学), 13, 182-184. (in Chinese)
35	Zhang H (张会), Du GS (杜桂森) (2002) The advances of molecular systematics on bryophytes. Journal of Capital Normal University (Natural Sciences Edition) (首都师范大学学报 (自然科学版) ), 23(2), 57-60. (in Chinese with English abstract)
36	Zhou KM(周匡明) (2006) The inspiration from the archeological spun silk pieces discovered at Qianshanyang site. In: Collectanea of Silk History (蚕业史论文选). Literature and History Press in China, Beijing. (in Chinese)
37	Zhou KM (周匡明) (1982) Study on the origin of sericulture cultivation. Agricultural Archaeology(农业考古), (1), 133-138. (in Chinese)
38	Zhou ZY (周泽扬), Li B (李斌), Ravikumar G, Quan GX (全国兴), Tamura T (田村俊树) (2004) Evolution studies on domestic silkworm (Bombyx mori L.) by AFLP markers. Southwest China Journal of Agricultural Sciences(西南农业学报), 17, 240-243. (in Chinese with English abstract)

序号 No.	品种 Variety	代号 Code	生态类型 Ecotype	序号 No.	品种 Variety	代号 Code	生态类型 Ecotype
1	小石丸 Xiaoshiwan	XSW	CV1	37	萨尼斯 Sanisi	SNS	EV1
2	二毛 Ermao	EM	CV1	38	罗尼7 Luoni7	LN7	EV1
3	土三虎斑 Tusanhuban	TSHB	CV1	39	法408 Fa408	F408	EV1
4	中土 Zhongtu	ZT	CV1	40	欧18 Ou18	O18	EV1
5	达县土种 Daxiantuzhong	DXTZ	CV1	41	意16 Yi16	Y16	EV1
6	土一肉黄圆 Tuyirouhuangyuan	TYRHY	CV1	42	乌D WuD	WD	EV2
7	土八 Tuba	Tba	CV1	43	乌G WuG	WG	EV2
8	鲁1 Lu1	L1	CV1	44	乌B WuB	WB	EV2
9	Re黑 Rehei	RH	CV1	45	乌2 Wu2	W2	EV2
10	鲁10 Lu10	L10	CV1	46	日110 Ri110	R110	T
11	土种01 Tuzhong01	TZ01	CV1	47	大造 Dazao	DZ	T
12	诸桂 Zhugui	ZG	CV1	48	高花 Gaohua	GH	T
13	温州种 Wenzhouzhong	WZZ	CV1	49	第五白卵 Diwubailuan	DWBL	T
14	芙蓉灰卵 Furonghuiluan	FRHL	CV1	50	演金黄白 Yanjinhuangbai	YJHB	T
15	灰色卵 Huiseluan	HSL	CV1	51	15-100	15-100	T
16	盱眙种 Xuyizhong	XYZ	CV1	52	07-310	07-310	T
17	余杭24 Yuhang24	YH24	CV1	53	高白 Gaobai	GB	T
18	沛县种 Peixianzhong	PXZ	CV1	54	琼山海南 Qiongshanhainan	QSHN	T
19	邯郸种 Handanzhong	HDZ	CV1	55	655	655	T
20	中4010 Zhong4010	Z4010	CV2	56	防四 Fangsi	FS	T
21	泰 Tai	T	CV2	57	印度种 Yinduzhong	YDZ	T
22	玉溪2号 Yuxi2hao	YX	CV2	58	加秋 Jiaqiu	JQ	T
23	平 Ping	P	CV2	59	白夏B BaixiaB	BXB	T
24	C108	C108	CV2	60	天龙青白 Tianlongqingbai	TLQB	JV1
25	余杭2 Yuhang2	YH2	CV2	61	赤熟 Chishu	Cshu	JV1
26	土白 Tubai	Tbai	CV2	62	日9 Ri9	R9	JV1
27	松花形吴 Songhuaxingwu	SHXW	CV2	63	姬蚕 Jican	JC	JV1
28	保黄 Baohuang	BH	CVd	64	黑子 Heizi	HZ	JV1
29	MCH101	MCH101	CVd	65	日10 Ri10	R10	JV1
30	BT924	BT924	CVd	66	绵蚕 Miancan	MC	JV1
31	BH863	BH863	CVd	67	日限2 Rixian2	RX2	JV2
表1 (续) Table 1 (continued)
序号 No.	品种 Variety	代号 Code	生态类型 Ecotype	序号 No.	品种 Variety	代号 Code	生态类型 Ecotype
32	琼十 Qiongshi	QS	CVd	68	赢纹形 Yingwenxing	YWX	JV2
33	四海 Sihai	SH	CVd	69	J115	J115	JV2
34	兰十 Lanshi	LanS	CVd	70	830	830	JV2
35	罗萨 Luosa	LuoS	EV1	71	春四 Chunsi	CS	JV2
36	阿利可斯 Alikesi	ALKS	EV1	72	T8	T8	JV2

序号 No.	品种 Variety	代号 Code	生态类型 Ecotype	序号 No.	品种 Variety	代号 Code	生态类型 Ecotype
1	小石丸 Xiaoshiwan	XSW	CV1	37	萨尼斯 Sanisi	SNS	EV1
2	二毛 Ermao	EM	CV1	38	罗尼7 Luoni7	LN7	EV1
3	土三虎斑 Tusanhuban	TSHB	CV1	39	法408 Fa408	F408	EV1
4	中土 Zhongtu	ZT	CV1	40	欧18 Ou18	O18	EV1
5	达县土种 Daxiantuzhong	DXTZ	CV1	41	意16 Yi16	Y16	EV1
6	土一肉黄圆 Tuyirouhuangyuan	TYRHY	CV1	42	乌D WuD	WD	EV2
7	土八 Tuba	Tba	CV1	43	乌G WuG	WG	EV2
8	鲁1 Lu1	L1	CV1	44	乌B WuB	WB	EV2
9	Re黑 Rehei	RH	CV1	45	乌2 Wu2	W2	EV2
10	鲁10 Lu10	L10	CV1	46	日110 Ri110	R110	T
11	土种01 Tuzhong01	TZ01	CV1	47	大造 Dazao	DZ	T
12	诸桂 Zhugui	ZG	CV1	48	高花 Gaohua	GH	T
13	温州种 Wenzhouzhong	WZZ	CV1	49	第五白卵 Diwubailuan	DWBL	T
14	芙蓉灰卵 Furonghuiluan	FRHL	CV1	50	演金黄白 Yanjinhuangbai	YJHB	T
15	灰色卵 Huiseluan	HSL	CV1	51	15-100	15-100	T
16	盱眙种 Xuyizhong	XYZ	CV1	52	07-310	07-310	T
17	余杭24 Yuhang24	YH24	CV1	53	高白 Gaobai	GB	T
18	沛县种 Peixianzhong	PXZ	CV1	54	琼山海南 Qiongshanhainan	QSHN	T
19	邯郸种 Handanzhong	HDZ	CV1	55	655	655	T
20	中4010 Zhong4010	Z4010	CV2	56	防四 Fangsi	FS	T
21	泰 Tai	T	CV2	57	印度种 Yinduzhong	YDZ	T
22	玉溪2号 Yuxi2hao	YX	CV2	58	加秋 Jiaqiu	JQ	T
23	平 Ping	P	CV2	59	白夏B BaixiaB	BXB	T
24	C108	C108	CV2	60	天龙青白 Tianlongqingbai	TLQB	JV1
25	余杭2 Yuhang2	YH2	CV2	61	赤熟 Chishu	Cshu	JV1
26	土白 Tubai	Tbai	CV2	62	日9 Ri9	R9	JV1
27	松花形吴 Songhuaxingwu	SHXW	CV2	63	姬蚕 Jican	JC	JV1
28	保黄 Baohuang	BH	CVd	64	黑子 Heizi	HZ	JV1
29	MCH101	MCH101	CVd	65	日10 Ri10	R10	JV1
30	BT924	BT924	CVd	66	绵蚕 Miancan	MC	JV1
31	BH863	BH863	CVd	67	日限2 Rixian2	RX2	JV2
表1 (续) Table 1 (continued)
序号 No.	品种 Variety	代号 Code	生态类型 Ecotype	序号 No.	品种 Variety	代号 Code	生态类型 Ecotype
32	琼十 Qiongshi	QS	CVd	68	赢纹形 Yingwenxing	YWX	JV2
33	四海 Sihai	SH	CVd	69	J115	J115	JV2
34	兰十 Lanshi	LanS	CVd	70	830	830	JV2
35	罗萨 Luosa	LuoS	EV1	71	春四 Chunsi	CS	JV2
36	阿利可斯 Alikesi	ALKS	EV1	72	T8	T8	JV2

	序列长度 Sequence length (bp)	碱基插入/缺失位点数 Insert/delete (%)	变异位点数 Variable site (%)	单个突变位点 Singleton variable site (%)	简约信息位点 Parsimony variable site (%)
外显子 Exon	493	2(0.4)*	10(2.0)	2(0.4)	8(1.6)
内含子 Intron	670	139(20.7)	55(8.2)	7(1.0)	48(7.2)
全序列 Entire	1,163	141(12.1)	65(5.6)	9(0.8)	56(4.8)

	序列长度 Sequence length (bp)	碱基插入/缺失位点数 Insert/delete (%)	变异位点数 Variable site (%)	单个突变位点 Singleton variable site (%)	简约信息位点 Parsimony variable site (%)
外显子 Exon	493	2(0.4)*	10(2.0)	2(0.4)	8(1.6)
内含子 Intron	670	139(20.7)	55(8.2)	7(1.0)	48(7.2)
全序列 Entire	1,163	141(12.1)	65(5.6)	9(0.8)	56(4.8)

	CV2	EV1	T	CV1M3	JV2	CV1M4	JV1	CVd
CV2
EV1	0.025
T	0.027	0.024
CV1M3	0.026	0.015	0.023
JV2	0.026	0.012	0.024	0.016
CV1M4	0.026	0.013	0.024	0.017	0.014
JV1	0.027	0.025	0.024	0.025	0.026	0.025
CVd	0.025	0.025	0.024	0.025	0.026	0.025	0.024
EV2	0.025	0.044	0.029	0.039	0.044	0.042	0.026	0.026