文章信息
- 武玉珍, 冯睿芝, 张峰
- WU Yuzhen, FENG Ruizhi, ZHANG Feng
- 基于ISSR分子标记的褐马鸡亲缘关系分析
- Analysis on genetic relationship of Brown-eared pheasant (Crossoptilon mantchuricum) based on ISSR molecular marker
- 生态学报, 2015, 35(4): 1059-1067
- Acta Ecologica Sinica, 2015, 35(4): 1059-1067
- http://dx.doi.org/10.5846/stxb201408051563
-
文章历史
- 收稿日期:2014-08-05
- 修订日期:2014-12-10
2. 复旦大学生物医学研究院, 上海 200032;
3. 山西大学生命科学院, 太原 030006
2. Institute of Biomedical Sciences, Fudan University, Shanghai 200032, China;
3. School of Life Science, Shanxi University, Taiyuan 030006, China
褐马鸡(Crossoptilon mantchuricum)为中国特有珍稀鸟类,国家一级保护动物,世界濒危物种。目前褐马鸡只分布在中国华北的部分地区,由于太行山和黄河的阻隔,褐马鸡的栖息地被分割成三块孤立的自然地理区域,形成了三个独立的褐马鸡种群,即山西吕梁山的中部种群、河北小五台山与北京东灵山的东部种群、陕西黄龙山的西部种群。山西省是褐马鸡分布中心和多度中心。长期以来,对于褐马鸡的研究主要集中在生物学、生态学、资源调查及保护等方面[1, 2, 3],近年来以分子生物学技术进行褐马鸡遗传多样性方面的研究[4, 5, 6, 7]有少量报道,表明对褐马鸡的研究开始从宏观进入微观分子水平。但采用ISSR分子标记对褐马鸡亲缘关系的研究目前尚未见报道。ISSR (Inter-Simple Sequence Repeats)是Zietkiewicz等[8]于1994年创建的一种新的分子标记技术,它是根据基因组中广泛存在的微卫星序列,设计通用引物对基因组DNA进行PCR扩增,从而检测基因组中微卫星序列的变异。ISSR结合了SSR和 RAPD 的优点,操作简单、稳定性好、检测方便,近年来广泛应用在种群遗传多样性[9, 10]、种质资源鉴定[11, 12]、亲缘关系[13, 14]和系统发育[15]等方面的研究。
本文采用ISSR分子标记技术,对山西庞泉沟国家自然保护区和太原动物园的两个褐马鸡种群的35个个体进行亲缘关系分析,旨在了解这两个分别代表山西野生种群和笼养种群褐马鸡的遗传背景和种质潜质,为保护褐马鸡种质资源和褐马鸡的野外放归与再引入提供科学的理论依据。
1 材料与方法 1.1 材料实验材料:庞泉沟褐马鸡样品18份,其中肌肉样品3份,采自庞泉沟保护区繁育场死亡幼雏;血液样品15份,采自保护区繁育场混养的成年褐马鸡。太原动物园血液样品17份,均采自动物园笼养种鸡,实验样品共35份。
1.2 方法 1.2.1 DNA提取分别取 0.34g肌肉、10—70 μL的血液,采用常规的SDS/蛋白酶K裂解,酚/氯仿法提取基因组总DNA,详细操作参照武玉珍[16]文献。
1.2.2 ISSR引物的筛选依据加拿大不列颠哥伦比亚大学(UBC)公布的第9套ISSR引物序列,参照Abbot[17]与白秀娟[18]的研究结果,选择其中的20条引物(表 1),均由北京奥科公司合成。通过4个随机样本的选择,从20条ISSR引物中筛选出10条扩增产物稳定、条带清晰、重复性好的引物,再对全部材料进行PCR扩增。
引物编号 Primer No. | 引物序列(5′-3′) Primer sequence | 退火温度/℃ Annealing temperature | 引物编号 Primer No. | 引物序列(5′-3′) Primer sequence | 退火温度/℃ Annealing temperature |
802 | ATATATATATATATATT | 36 | 815 | CTCTCTCTCTCTCTCTG | 54 |
803 | ATATATATATATATATC | 38 | 816 | CACACACACACACACAT | 52 |
805 | TATATATATATATATAC | 38 | 817 | CACACACACACACACAA | 52 |
806 | TATATATATATATATAG | 38 | 818 | CACACACACACACACAG | 54 |
807 | AGAGAGAGAGAGAGAGT | 52 | 820 | GTGTGTGTGTGTGTGTC | 54 |
808 | AGAGAGAGAGAGAGAGC | 54 | 822 | TCTCTCTCTCTCTCTCA | 52 |
809 | AGAGAGAGAGAGAGAGG | 54 | 824 | TCTCTCTCTCTCTCTCG | 54 |
811 | GAGAGAGAGAGAGAGAC | 54 | 827 | ACACACACACACACAG | 54 |
812 | GAGAGAGAGAGAGAGAA | 52 | 829 | TGTGTGTGTGTGTGTGC | 54 |
814 | CTCTCTCTCTCTCTCTA | 52 | 830 | TGTGTGTGTGTGTGTGG | 54 |
在25μL反应体系中含有10×PCR Buffer 2.5μL、25mmol/L Mg2+ 2μL、1.25 mmol/L dNTP 4μL、5.0 μmol/L引物5.0μL、2 U/μL Taq酶0.5μL、10 ng—150 ng/μL DNA模板1 μL,余量由ddH2O补足,在PTC-150 PCR仪上进行扩增。PCR反应条件为:94 ℃预变性5 min,94 ℃变性45 s,退火温度38—54 ℃ 45 s,72 ℃延伸1 min,40个循环,最后72 ℃延伸7min。每条引物的PCR反应退火温度见表 1。
1.2.4 电泳检测PCR产物用1.5%的琼脂糖凝胶进行电泳,DNA marker为DL2000 (购自大连宝生物公司),采用凝胶成像系统观察拍照,并记录电泳结果。
1.2.5 谱带分析ISSR扩增谱带以 0、1数据统计,在同一引物,同一位点,根据扩增产物的有(1)无(0)得到二元资料,形成ISSR分析的原始数据矩阵(0,1)矩阵。利用PhylTools6.0软件计算两个种群各个体之间的Nei′s无偏差遗传距离和遗传相似度,应用组平均法对褐马鸡35个个体的(0,1)原始矩阵进行聚类分析,得到聚类结果。
2 研究结果 2.1 ISSR扩增图谱本实验中共试验了20 条ISSR 引物,其中10条可扩增出清晰条带,将其用于供试35个样本PCR扩增,引物编号分别为806、807、808、809、811、812、815、816、817和818,除806外,引物的退火温度在50℃以上,保证了PCR 反应较高的特异性。10条引物共扩增出65条清晰稳定的条带,扩增片段大小在300—2000bp之间,每个引物扩增的位点从5到9不等,平均为6.5,65个扩增位点中有50个位点具有多态性,多态性位点比率为76.9%。引物对部分样品的扩增结果见图 1, 图 2, 图 3。
2.2 遗传距离与遗传相似系数利用PhylTools软件计算出了褐马鸡两个种群各个体之间的Nei′s遗传距离(表 2和表 3)、遗传相似系数(表 4和表 5)。
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | 15 | 16 | 17 |
1—17 分别为:dwy 1—17 号个体 | ||||||||||||||||
1 | ||||||||||||||||
2 | 0.2000 | |||||||||||||||
3 | 0.3333 | 0.3846 | ||||||||||||||
4 | 0.0909 | 0.1111 | 0.4286 | |||||||||||||
5 | 0.0909 | 0.1111 | 0.2857 | 0.2000 | ||||||||||||
6 | 0.3333 | 0.3846 | 0.1111 | 0.2857 | 0.4286 | |||||||||||
7 | 0.2000 | 0.0000 | 0.3846 | 0.1111 | 0.1111 | 0.3846 | ||||||||||
8 | 0.2000 | 0.2500 | 0.3846 | 0.3333 | 0.1111 | 0.5384 | 0.2500 | |||||||||
9 | 0.0909 | 0.1111 | 0.2857 | 0.2000 | 0.0000 | 0.4286 | 0.1111 | 0.1111 | ||||||||
10 | 0.5000 | 0.4000 | 0.2000 | 0.4545 | 0.4545 | 0.2000 | 0.4000 | 0.4000 | 0.4545 | |||||||
11 | 0.2308 | 0.2727 | 0.1250 | 0.3333 | 0.1667 | 0.2500 | 0.2727 | 0.2727 | 0.1667 | 0.2308 | ||||||
12 | 0.0909 | 0.1111 | 0.2857 | 0.2000 | 0.0000 | 0.4286 | 0.1111 | 0.1111 | 0.0000 | 0.4545 | 0.1667 | |||||
13 | 0.3333 | 0.3846 | 0.0000 | 0.4286 | 0.2857 | 0.1111 | 0.3846 | 0.3846 | 0.2857 | 0.2000 | 0.1250 | 0.2857 | ||||
14 | 0.3333 | 0.3846 | 0.0000 | 0.4286 | 0.2857 | 0.1111 | 0.3846 | 0.3846 | 0.2857 | 0.2000 | 0.1250 | 0.2857 | 0.0000 | |||
15 | 0.2857 | 0.3333 | 0.0588 | 0.3846 | 0.2308 | 0.1765 | 0.3333 | 0.3333 | 0.2308 | 0.1429 | 0.0667 | 0.2308 | 0.0588 | 0.0588 | ||
16 | 0.2857 | 0.3333 | 0.0588 | 0.3846 | 0.2308 | 0.1765 | 0.3333 | 0.3333 | 0.2308 | 0.1429 | 0.0667 | 0.2308 | 0.0588 | 0.0588 | 0.0000 | |
17 | 0.0909 | 0.1111 | 0.4286 | 0.0000 | 0.2000 | 0.2857 | 0.1111 | 0.3333 | 0.2000 | 0.4546 | 0.3333 | 0.2000 | 0.4286 | 0.4286 | 0.3846 | 0.3846 |
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | 15 | 16 | 17 | 18 |
1—18 分别为:pqg 1—18 号个体 | |||||||||||||||||
1 | |||||||||||||||||
2 | 0.5556 | ||||||||||||||||
3 | 0.4545 | 0.6667 | |||||||||||||||
4 | 0.5556 | 0.0000 | 0.6667 | ||||||||||||||
5 | 0.5000 | 1.0000 | 1.0000 | 1.0000 | |||||||||||||
6 | 0.4286 | 0.5556 | 0.2727 | 0.5556 | 0.6667 | ||||||||||||
7 | 0.5000 | 0.7143 | 0.5556 | 0.7143 | 0.6000 | 0.1667 | |||||||||||
8 | 0.5556 | 1.0000 | 0.3333 | 1.0000 | 1.0000 | 0.5556 | 0.4286 | ||||||||||
9 | 0.3846 | 0.5000 | 0.4000 | 0.5000 | 0.6364 | 0.0769 | 0.0909 | 0.5000 | |||||||||
10 | 0.2500 | 0.8182 | 0.3846 | 0.8182 | 0.2857 | 0.2500 | 0.4286 | 0.6364 | 0.3333 | ||||||||
11 | 0.3333 | 0.7143 | 0.5556 | 0.7143 | 0.4000 | 0.6667 | 1.0000 | 1.0000 | 0.8182 | 0.2857 | |||||||
12 | 0.4545 | 0.6667 | 0.0000 | 0.6667 | 1.0000 | 0.2727 | 0.5556 | 0.3333 | 0.4000 | 0.3846 | 0.5556 | ||||||
13 | 0.2000 | 0.8000 | 0.3333 | 0.8000 | 0.3846 | 0.3333 | 0.5385 | 0.6000 | 0.4286 | 0.0588 | 0.2308 | 0.3333 | |||||
14 | 0.2857 | 1.0000 | 0.4545 | 1.0000 | 0.3333 | 0.4286 | 0.5000 | 0.5556 | 0.5385 | 0.1250 | 0.3333 | 0.4545 | 0.0667 | ||||
15 | 0.2857 | 1.0000 | 0.4545 | 1.0000 | 0.3333 | 0.4286 | 0.5000 | 0.5556 | 0.5385 | 0.1250 | 0.3333 | 0.4545 | 0.0667 | 0.0000 | |||
16 | 0.6364 | 1.0000 | 0.5000 | 1.0000 | 0.5556 | 0.2727 | 0.1111 | 0.3333 | 0.2000 | 0.3846 | 1.0000 | 0.5000 | 0.5000 | 0.4545 | 0.4545 | ||
17 | 0.3333 | 1.0000 | 0.5000 | 1.0000 | 0.2308 | 0.3333 | 0.3846 | 0.6000 | 0.4286 | 0.0588 | 0.3846 | 0.5000 | 0.1250 | 0.0667 | 0.0667 | 0.3333 | |
18 | 0.1667 | 1.0000 | 0.5556 | 1.0000 | 0.4000 | 0.6667 | 0.6000 | 0.4286 | 0.6364 | 0.2857 | 0.4000 | 0.5556 | 0.2308 | 0.1667 | 0.1667 | 0.5556 | 0.2308 |
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | 15 | 16 | 17 | |
1—17 分别为:dwy 1—17 号个体 | |||||||||||||||||
1 | |||||||||||||||||
2 | 0.8000 | ||||||||||||||||
3 | 0.6667 | 0.6154 | |||||||||||||||
4 | 0.9091 | 0.8889 | 0.5714 | ||||||||||||||
5 | 0.9091 | 0.8889 | 0.7143 | 0.8000 | |||||||||||||
6 | 0.6667 | 0.6154 | 0.8889 | 0.7143 | 0.5714 | ||||||||||||
7 | 0.8000 | 1.0000 | 0.6154 | 0.8889 | 0.8889 | 0.6154 | |||||||||||
8 | 0.8000 | 0.7500 | 0.6154 | 0.6667 | 0.8889 | 0.4615 | 0.7500 | ||||||||||
9 | 0.9091 | 0.8889 | 0.7143 | 0.8000 | 1.0000 | 0.5714 | 0.8889 | 0.8889 | |||||||||
10 | 0.5000 | 0.6000 | 0.8000 | 0.5455 | 0.5455 | 0.8000 | 0.6000 | 0.6000 | 0.5455 | ||||||||
11 | 0.7692 | 0.7273 | 0.8750 | 0.6667 | 0.8333 | 0.7500 | 0.7273 | 0.7273 | 0.8333 | 0.7692 | |||||||
12 | 0.9091 | 0.8889 | 0.7143 | 0.8000 | 1.0000 | 0.5714 | 0.8889 | 0.8889 | 1.0000 | 0.5455 | 0.8333 | ||||||
13 | 0.6667 | 0.6154 | 1.0000 | 0.5714 | 0.7143 | 0.8889 | 0.6154 | 0.6154 | 0.7143 | 0.8000 | 0.8750 | 0.7143 | |||||
14 | 1.0000 | 0.6667 | 0.6154 | 1.0000 | 0.5714 | 0.7143 | 0.8889 | 0.6154 | 0.6154 | 0.7143 | 0.8000 | 0.8750 | 0.7143 | ||||
15 | 0.7143 | 0.6667 | 0.9412 | 0.6154 | 0.7692 | 0.8235 | 0.6667 | 0.6667 | 0.7692 | 0.8571 | 0.9333 | 0.7692 | 0.9412 | 0.9412 | |||
16 | 0.7143 | 0.6667 | 0.9412 | 0.6154 | 0.7692 | 0.8235 | 0.6667 | 0.6667 | 0.7692 | 0.8571 | 0.9333 | 0.7692 | 0.9412 | 0.9412 | 1.0000 | ||
17 | 0.9091 | 0.8889 | 0.5714 | 1.0000 | 0.8000 | 0.7143 | 0.8889 | 0.6667 | 0.8000 | 0.5455 | 0.6667 | 0.8000 | 0.5714 | 0.5714 | 0.6154 | 0.6154 |
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | 15 | 16 | 17 | 18 |
1—18 分别为:pqg 1—18号个体 | |||||||||||||||||
1 | |||||||||||||||||
2 | 0.4444 | ||||||||||||||||
3 | 0.5455 | 0.3333 | |||||||||||||||
4 | 0.4444 | 1.0000 | 0.3333 | ||||||||||||||
5 | 0.5000 | 0.0000 | 0.0000 | 0.0000 | |||||||||||||
6 | 0.5714 | 0.4444 | 0.7273 | 0.4444 | 0.3333 | ||||||||||||
7 | 0.5000 | 0.2857 | 0.4444 | 0.2857 | 0.4000 | 0.8333 | |||||||||||
8 | 0.4444 | 0.0000 | 0.6667 | 0.0000 | 0.0000 | 0.4444 | 0.5714 | ||||||||||
9 | 0.6154 | 0.5000 | 0.6000 | 0.5000 | 0.3636 | 0.9231 | 0.9091 | 0.5000 | |||||||||
10 | 0.7500 | 0.1818 | 0.6154 | 0.1818 | 0.7143 | 0.7500 | 0.5714 | 0.3636 | 0.6667 | ||||||||
11 | 0.6667 | 0.2857 | 0.4444 | 0.2857 | 0.6000 | 0.3333 | 0.0000 | 0.0000 | 0.1818 | 0.7143 | |||||||
12 | 0.5455 | 0.3333 | 1.0000 | 0.3333 | 0.0000 | 0.7273 | 0.4444 | 0.6667 | 0.6000 | 0.6154 | 0.4444 | ||||||
13 | 0.8000 | 0.2000 | 0.6667 | 0.2000 | 0.6154 | 0.6667 | 0.4615 | 0.4000 | 0.5714 | 0.9412 | 0.7692 | 0.6667 | |||||
14 | 0.7143 | 0.0000 | 0.5455 | 0.0000 | 0.6667 | 0.5714 | 0.5000 | 0.4444 | 0.4615 | 0.8750 | 0.6667 | 0.5455 | 0.9333 | ||||
15 | 0.7143 | 0.0000 | 0.5455 | 0.0000 | 0.6667 | 0.5714 | 0.5000 | 0.4444 | 0.4615 | 0.8750 | 0.6667 | 0.5455 | 0.9333 | 1.0000 | |||
16 | 0.3636 | 0.0000 | 0.5000 | 0.0000 | 0.4444 | 0.7273 | 0.8889 | 0.6667 | 0.8000 | 0.6154 | 0.0000 | 0.5000 | 0.5000 | 0.5455 | 0.5455 | ||
17 | 0.6667 | 0.0000 | 0.5000 | 0.0000 | 0.7692 | 0.6667 | 0.6154 | 0.4000 | 0.5714 | 0.9412 | 0.6154 | 0.5000 | 0.8750 | 0.9333 | 0.9333 | 0.6667 | |
18 | 0.8333 | 0.0000 | 0.4444 | 0.0000 | 0.6000 | 0.3333 | 0.4000 | 0.5714 | 0.3636 | 0.7143 | 0.6000 | 0.4444 | 0.7692 | 0.8333 | 0.8333 | 0.4444 | 0.7692 |
从表 2可以看出:太原动物园17个个体间的遗传距离为0.0000—0.5348,平均值为0.2409。其中,dwy8与dwy6之间遗传距离最大为0.5348,其次是dwy1与 dwy10间为0.5000,而dwy7与dwy2,dwy9与dwy5,dwy12分别与dwy5、dwy9,dwy13与dwy3,dwy14分别与dwy3、dwy13,dwy16与dwy15,dwy17与dwy4间遗传距离为0.0000。
从表 3可以看出:庞泉沟保护区18个个体间的遗传距离在0.000—1.000之间,平均值为0.4939。其中,pqg2与pqg5、pqg 8、pqg14—18,pqg3与pqg5,pqg4与pqg5、pqg 8、pqg14—18,pqg5与pqg8、pqg12,pqg7与pqg11,pqg8与pqg11,pqg11与pqg16遗传距离最远为1.0000;其次是pqg10与pqg2、pqg4以及pqg9与pqg11间遗传距离为0.8182;而pqg4与pqg2、pqg12与pqg3、pqg15与pqg14遗传距离最近为0.0000。
从表 4可以看出:太原动物园种群中,dwy7与dwy2,dwy9与dwy5,dwy12分别与dwy5、dwy9,dwy13与dwy3,dwy14分别与dwy1、dwy4,dwy16与dwy15,dwy17与dwy4之间遗传相似度最高为1.0000,遗传距离最近;而dwy8与dwy6之间遗传相似度最低为0.4615,遗传距离最远。
从表 5可以看出:庞泉沟自然保护区种群中,pqg4与pqg2,pqg12与pqg3,pqg15与pqg14遗传相似度最高为1.0000,遗传距离最近;而pqg2、pqg4分别与pqg5、pqg8、pqg14、pqg15、pqg16、pqg17、pqg18,pqg5分别与pqg3、pqg8、 pqg12,pqg11与pqp8、pqg7,遗传一致度最低为0.0000,遗传距离最远。
2.3 亲缘关系聚类图应用组平均法对褐马鸡35个个体的(0,1)原始矩阵进行聚类分析,结果见图 4。从图 4看出:35个个体聚为3类,其中24(pqg7)为单独一类;太原动物园除dwy7外,其它个体聚为一类;庞泉沟自然保护区的个体聚为一类。大多数同一种群的不同个体首先聚在一起,显示出明确的生态地理条件相似的种群优先相聚,也就是说生态地理条件相似的种群亲缘关系较近。pqg7显示出与两个种群的其他个体亲缘关系相对较远;dwy6与动物园种群的其他个体亲缘关系相对较远;dwy7显示出与庞泉沟种群具有较近的亲缘关系。
3 分析与讨论 3.1 ISSR实验的稳定性要获得一个物种不同种群的遗传背景信息,必需检测大量的样本,检测方法的稳定性也是非常重要的。ISSR引物一般是由17bp的碱基组成,对基因组DNA具有快速、高效、灵敏的检测特性,而且可重复性好,操作简便、成本较低,实验的稳定性较高[19]。由于ISSR的稳定性受Mg2+浓度、引物浓度、退火温度等因素的影响,因此,进行ISSR大量样本的扩增前,先对ISSR反应的条件进行了优化,采用统一的PCR反应条件。对20条ISSR引物进行了筛选,其中有10条引物能扩增出条带清晰且具多态性的带型。在这10条引物中(AG)n与(CA)n各有3条,(GA)n有2条,说明褐马鸡基因组中存在大量的(CA)、(AG)二核苷酸重复序列。在引物筛选的过程中,退火温度较高的引物扩增效果均较退火温度低的好,筛选出所选用的10条引物只有1条退火温度较低,其它9条引物退火温度均较高。在ISSR图谱中出现某种群的特征条带,引物815的扩增图谱只有庞泉沟种群的个体出现大于1000 bp 的片段,引物811、812、809对庞泉沟自然保护区种群的扩增有特异性。引物811可以将这35个个体区分开。因此,只要引物选择合适,扩增出种群的特征性谱带,有助于种群的遗传结构分析,为进一步克隆、测序,寻找种群特有的分子标记提供了可能。
应用mtDNA标记和ISSR标记两种方法,对这两个褐马鸡种群进行了遗传多样性的研究[5, 7],二者的结果基本一致,这表明无论用线粒体DNA标记(mtDNA),还是细胞核DNA标记(ISSR),研究的结果基本相同,可以互为印证。这说明ISSR标记是一种研究种群遗传多样性、遗传结构、种质资源及亲缘关系的稳定有效工具。
3.2 亲缘关系分析山西庞泉沟国家自然保护区与太原动物园两个褐马鸡种群个体间遗传相似性的平均值分别为0.5061与0.7591,遗传距离的平均值分别为:0.4939与0.2409,35只褐马鸡个体总的Nei′s遗传距离和相似性系数平均值分别为:0.3723和0.6310。分析表明两个种群的遗传变异有差异,但差异未达到统计学显著性水平,两个种群个体之间亲缘关系较近。聚类分析结果表明:遗传距离较近的个体聚在一起,然后与遗传距离较远的个体再聚。遗传距离较近的个体亲缘关系近,而遗传距离较远的个体,亲缘关系较远。大多数同一种群的不同个体首先聚在一起,显示出生态地理条件相似的种群优先相聚,但也存在不同种群聚在一起的情况。如图 4中pqg7单独聚为一类,显示出与两个种群亲缘关系较远,推测可能是来自保护区分布较远的不同种群的后代。图 4中dwy7与庞泉沟种群聚为一类,显示出与庞泉沟种群有较近的亲缘关系,推测这两个褐马鸡种群可能具有共同的亲源。这一结论与太原动物园饲养褐马鸡的历史相符合。太原动物园1959年由芦芽山(宁武县)提供褐马鸡种源,开始试验人工饲养,曾经4次较大数量补充野生种源,1979—1983年从芦芽山补充60余只,1992年和2014年分别从庞泉沟自然保护区补充野生种源8只、6只,2011年从芦芽山补充野生卵孵化9只。到2000年建成了全国最大的褐马鸡人工饲养繁育基地,为全国各地动物园提供褐马鸡种源和观赏鸡。
两个种群遗传相似性和遗传距离的差异(表 2—表 5),反应出两个种群饲养方式的不同。庞泉沟自然保护区的个体间遗传相似性的平均值低于太原动物园,遗传距离的平均值高于太原动物园,说明庞泉沟种群个体间的亲缘关系较远,动物园种群个体之间的亲缘关系较近。这种状况与两个种群不同的饲养方式有关。庞泉沟自然保护区繁育场采取仿野生的散养、混养方式,并且每年从保护区野外捡拾褐马鸡卵进行孵化补充种群,繁育场的种群虽然有可能近亲繁殖,但因为不断有新种源补充,群体的亲缘关系仍较远。太原动物园从1996年起选配12对成年褐马鸡实行“一夫一妻”的分笼配对繁殖饲养,这种方式有效避免了近交繁殖,但是种源有限,更新也较慢,群体的亲缘关系仍比较近。因为饲养的用途不同,所采取的饲养方式也不同。
本试验对褐马鸡两个种群35个个体进行遗传距离聚类及亲缘关系分析,可作为优质种源和再引入亲本选择时的理论依据。作为亲本应尽量选择遗传距离大于0.5的组合,同时避免遗传距离小于 0.25的亲本组合。可重点选择庞泉沟自然保护区野生种群的pqg2、pqg4分别与pqg5、pqg8、pqg14、pqg15、pqg16、pqg17、pqg18,pqg5与pqg3、pqg8、 pqg12,pqg11与pqp7、pqg8等,这些个体之间遗传相似性最低(0.0000),遗传距离最远(1.0000),亲缘关系最远,可以作为优质亲本种源。
3.3 ISSR标记在野生动物研究中的应用前景ISSR标记常用作一种植物(农作物、药材等)多个品种、品系之间亲缘关系的区分和鉴定。但野生动物个体之间外部形态差异不大,亲缘关系主要与生态地理环境有关,特别是像褐马鸡这种飞行能力较差的地栖鸟类。应用ISSR标记技术可以有效分辨动物个体之间的亲缘关系,同时可以筛选出亲缘关系较远的个体,如独自聚为一类的pqg7个体(图 4),作为人工饲养繁育的优质种源和褐马鸡再引入的原始种源。栖息地的破碎化致使我国的褐马鸡被分割为3个孤立种群(山西、河北、陕西),各种群之间因地理阻隔难有基因交流,遗传多样性较低[5, 7]。以ISSR标记研究褐马鸡亲缘关系、寻找发现优质种源,对褐马鸡的种群发展和科学保护具有重要的学术价值和现实意义。
ISSR标记作为以PCR为基础的分子标记,除具有操作简单、成本低、快速、灵敏的优点外,还有所需 DNA模板量少、稳定和多态条带丰富等优点,是研究物种遗传背景的有效工具。ISSR标记已广泛应用于植物的遗传连锁图谱构建、种质资源鉴定、基因定位、植物分类、进化和遗传多样性等研究[20],而在动物研究中应用较少[21],仅有少数家禽与水产鱼虾蟹的研究。ISSR标记用于濒危野生动物的研究则更少,仅见东北虎[18] (Panthera tigris altaica)、金钱豹[22] (Panthera pardus)和北海狮[23] (Eumetopias jubatus)、獐(Hydropotes inermis) [24]等研究,而用于珍稀野生鸟类的研究鲜见报道。本研究将ISSR标记技术用于中国特有的濒危鸟类——褐马鸡的亲缘关系研究,扩展了鸟类遗传背景和系统进化等方面的研究手段,是有意义的尝试和探索,对其他野生动物的遗传分析也不失为一种便捷稳定的方法。由于采样数量及引物数量较少,检测到的位点有限,对褐马鸡圈养和野生种群亲缘关系的全面评价,有待于进一步深入研究。
致谢: 山西省自然保护区管理站张龙胜教授、山西庞泉沟国家自然保护区管理站武建勇高级工程师及太原动物园卫泽珍副主任、孟庆珍兽医师等对样品采集中给予帮助,特此致谢。[1] | 李宏群, 廉振民, 刘晓莉. 中国褐马鸡的研究现状及其保护措施. 延安大学学报: 自然科学版, 2009, 28(2): 92-96. |
[2] | 张国钢, 郑光美, 张正旺, 郭建荣, 王建平, 宫树龙. 山西芦芽山褐马鸡越冬栖息地选择的多尺度研究. 生态学报, 2005, 25(5): 952-957. |
[3] | 武玉珍, 冯睿芝. 褐马鸡的濒危原因及保护措施. 中国家禽, 2013, 35(8): 49-50. |
[4] | 冯宁, 杨淑娟, 徐振武, 高学斌, 王万云, 尹祚华. 褐马鸡地理种群的遗传多样性研究. 陕西师范大学学报: 自然科学版, 2007, 35(S1): 77-80. |
[5] | 武玉珍, 王孟本, 张峰. 褐马鸡圈养种群的mtDNA控制区多态性. 生态学报, 2010, 30(11): 2958-2964. |
[6] | 武玉珍, 冯睿芝, 张峰. 珍禽褐马鸡线粒体DNA控制区结构和亲缘关系. 生态学杂志, 2013, 32(12): 3243-3249. |
[7] | 武玉珍, 冯睿芝. 褐马鸡圈养种群ISSR标记的遗传多样性. 林业科学, 2013, 49(11): 103-108. |
[8] | Zietkiewicz E, Rafalske A, Labuda D. Genome fingerprinting by simple sequence repeat (SSR)-anchored polymerase chain reaction amplification. Genomics, 1994, 20(2): 178-183. |
[9] | 吴则焰, 刘金福, 洪伟, 潘东明, 郑世群, 何中声. 孑遗植物水松不同年龄级种群遗传多样性的ISSR分析. 生态学杂志, 2012, 31(8): 1911-1916. |
[10] | 张杰, 王佳, 李浩宇, 张慧荣, 王迎春. 濒危植物蒙古扁桃不同地理种群遗传多样性的ISSR分析. 生态学报, 2012, 32(14): 4443-4452. |
[11] | 孙利娜, 席梦利, 徐进, 吴祝华, 施季森. 百合辐射突变体的ISSR鉴定技术研究. 现代农业科技, 2011, (1): 219-221. |
[12] | 范付华, 杨勇胜, 李庆宏, 邓勇, 文晓鹏. 利用ISSR分析技术鉴定贵州枇杷新种质的研究. 北京林业大学学报, 2012, 34(4): 52-57. |
[13] | 刘本英, 李友勇, 唐一春, 王丽鸳, 成浩, 王平盛. 云南茶树资源遗传多样性与亲缘关系的ISSR分析. 作物学报, 2010, 36(3): 391-400. |
[14] | 严华, 张冬梅, 罗玉兰, 鲁琳. 38种国兰亲缘关系的ISSR分析. 分子植物育种, 2010, 8(4): 736-741. |
[15] | 权俊萍, 何树兰, 彭峰, 郑玉红, 夏冰. 部分百里香属植物分子系统学研究. 西北植物学报, 2008, 28(8): 1566-1572. |
[16] | 武玉珍. 濒危鸟类褐马鸡遗传多样性及保护研究 [D]. 太原: 山西大学, 2008: 12-18. |
[17] | Abbot P. Individual and population variation in invertebrates revealed by inter-simple sequence repeats (ISSR). The Journal of Insect Science, 2001, 1: 8. |
[18] | 白秀娟. 圈养东北虎ISSR指纹分析初报. 兽类学报, 2004, 24(1): 90-92. |
[19] | 郑芳, 吕秀玲, 孙红英, 周开亚, 赵强, 高伟, 徐信荣. ISSR标记在河蟹种质检测中的应用. 中国水产科学, 2007, 14(1): 46-51. |
[20] | 周延清. DNA分子标记技术在植物研究中的应用. 北京: 化学工业出版社, 2005: 143-161. |
[21] | 林杰君, 鲍毅新, 刘军, 张旭. ISSR分子标记及其在动物遗传结构研究中的应用. 生态学杂志, 2012, 31(5): 1319-1326. |
[22] | 朱军, 王立斌, 郭东龙, 马恩波, 任竹梅. 基于ISSR标记技术的金钱豹遗传多样性分析. 山西大学学报: 自然科学版, 2008, 31(2): 244-247. |
[23] | Koyama S, Fujita S, Hirota T, Satoh T, Obara Y, Hoshino H, Wada A, Burkanov V N, Wada K. Genetic structure of steller sea lion (Eumetopias jubatus) rookeries in the sea of Okhotsk. Zoological Studies, 2008, 47(6): 781-787. |
[24] | 林杰君, 鲍毅新, 刘军, 王艳妮, 张旭. 舟山群岛四个养殖獐种群遗传多样性和遗传结构. 生态学报, 2013, 33(11): 3460-3469. |