CN101974629A

CN101974629A - 一种虚拟合成物种调查异源多倍体物种起源的方法

Info

Publication number: CN101974629A
Application number: CN 201010519822
Authority: CN
Inventors: 钱伟; 梅家琴; 李加纳; 李勤菲; 钱论文; 傅鹰
Original assignee: Southwest University
Current assignee: Southwest University
Priority date: 2010-10-26
Filing date: 2010-10-26
Publication date: 2011-02-16

Abstract

本发明涉及一种虚拟合成物种调查异源多倍体物种起源的方法，其包括：1)采用PCR技术对自然异源多倍体及其亲本物种进行全基因组水平的分子标记分析；2)根据合成异源多倍体过程中，亲本具有的分子特征能在合成后代中得到体现的原则，利用亲本物种的分子标记数据，模拟合成种的分子标记数据；3)利用遗传关系分析软件和数据统计分析软件，分析由不同亲本种类型模拟合成的虚拟异源多倍体与自然异源多倍体的亲缘关系，根据遗传差异大小，推导异源多倍体中亚基因组的来源。本发明避免了人工合成的繁琐过程，可同时对待测亲本物种的所有变异类型进行分析，亦能扣除另一(其它)亲本种基因组对分析结果的影响，方法简单，时间周期短，耗资少，节省人力物力。

Description

一种虚拟合成物种调查异源多倍体物种起源的方法

技术领域

本发明属于生物技术领域，具体涉及一种虚拟合成物种调查异源多倍体物种起源的方法。适用于调查亲本物种变异类型较多的各种异源多倍体物种的起源。

背景技术

多倍体物种广泛存在于高等植物中,根据染色体组是否来源于相同物种,可分为同源多倍体和异源多倍体。很多重要作物比如甘蓝型油菜、陆地棉和普通栽培小麦都是异源多倍体物种，其染色体组来自于至少两个不同的亲本物种。

尽管很多异源多倍体物种的染色体组组成已被揭示，但是由于亲本种的变异类型众多，遗传变异广泛，人们对于异源多倍体物种各亚基因组的供体亲本类型并不是很清楚，由此也引起了对其起源和进化上的争论。

调查异源多倍体亲缘关系和起源可通过种间杂交、细胞学和分子生物学方法等方法进行, 其中一种重要的方法是：利用亲本物种合成并重演异源多倍体物种的形成。该方法主要包括以下过程：利用亲本物种进行种间杂交，人工合成异源多倍体，再与现有的自然异源多倍体进行分子生物学等方面的比较，最后根据来源于不同类型的亲本物种的人工合成种与自然异源多倍体的相似程度来推测该异源多倍体的原始祖先种类型。该方法能较大程度地模拟异源多倍体遗传物质的起源，是一种较为直接的方法。但该法存在以下主要缺点：人工合成异源多倍体过程繁琐，涉及远缘杂交，胚挽救，染色体加倍等程序，工作量大，且难以成功获得源于所有亲本物种类型的人工合成种。

另一种方法是，对亲本物种不同类型和自然异源多倍体进行比较，根据亲本种各变异类型与自然异源多倍体遗传差异程度，推测该异源多倍体的亲本种类型。该法虽然简单易行，但异源多倍体包含了来源于两个或更多物种的不同基因组，在分析某一种亲本物种与异源多倍体关系的过程中，异源多倍体中源于另一（其它）亲本物种的基因组会对分析造成干扰。

现代分子生物学的快速发展为研究异源多倍体的起源提供新的途径，如通过测序比较亲本物种和自然异源多倍体某些基因的DNA序列差异，通过细胞质DNA的分子标记，比较亲本物种和自然异源多倍体细胞质遗传物质的差异等，但以上该方法均需事前获得某特异基因的序列和开放出特异的引物，并且其揭示的信息量在基因组水平上非常有限。

发明内容

本发明的目的在于提供了一种虚拟合成物种调查异源多倍体物种起源的方法，该方法结合现代分子生物学方法模拟人工合成异源多倍体物种的分子特征，为研究异源多倍体的起源提供了新的途径，本发明的主要特点是：采用了人工合成异源多倍体的原理，但避免了人工合成的繁琐过程，可同时对待测亲本物种的所有变异类型进行分析，亦能扣除另一（其它）亲本种基因组对分析结果的影响，方法简单，时间周期短，耗资少，节省人力物力。

为了达到上述目的，本发明采用以下技术措施：

本发明的原理是：很多重要作物都是异源多倍体物种，其起源和进化上一直存在着争论。研究异源多倍体起源物种起源的多种方法中，利用亲本种合成并重演物种的形成是研究其起源最直接和重要的方法，但是该法过程复杂，且难以将亲本种的所有变异类型都囊括在内。基于该法，本发明阐述了一种利用虚拟合成物种研究异源多倍体起源的方法，其基本原理是：在合成异源多倍体物种的过程中，亲本物种的遗传物质传递给后代，亲本物种具有的分子特征能在合成后代中得到体现，因此对待测亲本不同变异类型和另一亲本物种进行全基因组的分子扫描，继而模拟出不同来源的异源多倍体的分子特征（即虚拟合成异源多倍体），并与自然异源多倍体材料的分子特征进行比较，根据遗传差异的大小推导自然异源多倍体物种的亲本类型。基于虚拟合成物种的策略研究异源多倍体物种起源的优点是，该法能扣除另一（其它）亲本物种基因组对分析结果的影响；能对亲本种的任何变异类型进行异源多倍体的虚拟合成，可在尽可能大的范围内寻找异源多倍体原始亲本供体类型；无需经历杂交、胚挽救、染色体加倍等过程，省时省力。

本发明提出的虚拟合成物种调查异源多倍体物种起源的方法步骤如下：

1）将亲本物种及自然异源多倍体全基因组分子扫描：采用SSR、AFLP等常规PCR技术对自然异源多倍体及其亲本物种进行全基因组水平分子标记特征扫描；

2）虚拟合成异源多倍体：根据合成异源多倍体过程中，亲本具有的分子特征能在合成后代中得到体现的原则，双亲或亲本之一中标记为“1”的则在虚拟合成异源多倍体中记为“1”，在双亲中都为“0”的则在虚拟合成异源多倍体中记为“0”，据此利用双亲的分子标记数据模拟出合成种的分子标记数据；

3）亲缘关系分析：利用NTSYS软件计算各亲本与自然异源多倍体，以及各虚拟合成异源多倍体与自然异源多倍体间的遗传距离，并对分子标记数据进行UPGMA聚类分析和主成分分析（PCA）；用STRUCTURE软件分析二倍体亲本的不同变异类型与自然异源多倍体之间的血缘关系；用ARLEQUIN软件，按亲本物种进行分子标记方差分析(AMOVA)，检验各亲本物种对虚拟合成异源多倍体遗传差异的贡献率；用SAS软件对来源于同一亲本物种内不同个体的虚拟合成异源多倍体之间的遗传距离进行方差分析（ANOVA），并对亲本与自然异源多倍体、虚拟合成异源多倍体与自然异源多倍体间的遗传距离进行方差分析及多重比较。通过以上分析揭示不同亲本种类型模拟合成的虚拟异源多倍体与自然异源多倍体相互间遗传差异的大小，根据遗传差异越小，该亲本类型与异源多倍体越为接近，则该类型越可能是异源多倍体中亚基因组的来源的原则，推导异源多倍体中亚基因组的来源。

本发明的特征在于上述步骤的2）-3）。步骤2）根据合成异源多倍体过程中，亲本具有的分子特征能在合成后代中得到体现的原则，利用双亲的分子标记数据模拟出合成种的分子标记数据。该步骤避免了传统的人工合成异源多倍体的繁琐过程，能模拟出所有亲本物种变异类型所产生的异源多倍体的标记特征，从而做到在更广的范围里寻找异源多倍体种亚基因组的来源；克服了直接利用亲本种与自然异源多倍体进行遗传差异分析的不足，在分析待测亲本物种与异源多倍体的关系时，通过虚拟合成多倍体物种将另一亲本种基因组的分子特征纳入其中，再经步骤3）与自然异源多倍体进行比较，从而扣除另一亲本种基因组对遗传差异的影响。

本发明针对的多倍体为异源多倍体，即染色体来源于至少两种不同物种的多倍体。本发明适用于亲本物种变异类型较多的异源多倍体物种的祖先种类型推测，是一种适用范围广，简单高效的调查异源多倍体起源的方法。

本发明涉及的分子标记可以是基于PCR技术的全基因组水平分子标记种的任意一种或多种，如AFLP, SSR, RAPD, ISSR等；所采用的数据分析软件可以是SAS、SPSS等软件及其他具有类似功能的软件；所采用的遗传关系分析软件可以选择NT-SYS、STRUCTURE等软件，其它具有相应功能的软件均可采用。本发明所涉及的实验方法如基于PCR的全基因组水平分子标记，遗传关系的分析均为分子生物学常规技术手段，为已有技术；各类分析软件可通过相应网站获得；各类试剂均可从商业途径获得。

本发明优点在于：可同时对待测亲本物种的所有变异类型进行分析，能扣除另一亲本种基因组对遗传差异的影响，简单高效，省时省力。

附图说明

图1 为利用虚拟合成物种研究异源多倍体物种起源的原理示意图。通过比较虚拟合成多倍体和自然异源多倍体的全基因组水平分特征，揭示两者间的遗传差异，从而推导该自然异源多倍体的亲本种类型。双箭头表示两者之间的标记特征可相互推导，直线表示两者之间进行遗传差异比较。

图2 为虚拟异源多倍体的合成原理示意图。箭头左边部分为合成异源多倍体物种中双亲染色体的传递行为：不同物种合成异源多倍体时，双亲的遗传物质均传递给后代，亲本种具有的分子特征能在合成后代中体现。其中，P₁和P₂分别代表亲本1和亲本2，竖立长方条代表染色体，两侧的阿拉伯数字表示标记位点，相同的数字代表位点相同，长方条中出现横线代表某位点上的标记类型为“1”（如亲本1中的位点1和位点2，以及亲本2中的位点1和位点3），没有横线则代表该位点标记类型为“0”（如亲本1和亲本2中的位点4）。箭头右边部分为根据双亲标记类型模拟合成异源多倍体标记类型的过程：对于虚拟合成异源多倍体，在双亲中同时存在（如位点1）或在双亲之一中存在（如位点2和位点3）的标记类型记为“1”，在双亲中都没有（如位点4）的标记类型记为“0”。

图3为甘蓝、白菜型油菜和自然甘蓝型油菜基于SSR和AFLP标记的聚类关系图。三个物种材料分别各自成群，其中甘蓝的十七种变异类型可划分为三个主要类群：甘蓝类群I，甘蓝类群II和甘蓝类群III。

图4为甘蓝、白菜型油菜和自然甘蓝型油菜间的血缘关系示意图。横坐标的阿拉伯数字为材料编号（与图3相同），横线下的罗马数字I ~ V分别代表甘蓝类群I、甘蓝类群II、甘蓝类群III、白菜型油菜和自然甘蓝型油菜，纵坐标代表血缘所占比例；图中不同的颜色代表不同的血缘。结果显示，按照血缘所分类群与聚类分析中所划分类群（图3）相符，三个甘蓝类群间存在三种截然不同的血缘组成，自然甘蓝型油菜具有部分与甘蓝类群I中甘蓝相同的血缘（浅灰色）。

图5为甘蓝、白菜型油菜、虚拟甘蓝型油菜和自然甘蓝型油菜基于SSR和AFLP标记的主成分分析图（PCA）。横坐标表示第一主成分，纵坐标表示第二主成分；正方形代表白菜型油菜，四角星代表自然甘蓝型油菜，三角形代表甘蓝（其中白色、浅灰色和深灰色三角形分别代表甘蓝类群I、II 和 III中的甘蓝），圆形代表虚拟合成甘蓝型油菜（其中白色、浅灰色和深灰色圆分别代表由甘蓝类群I、II 和 III中的甘蓝与白菜型油菜虚拟合成的甘蓝型油菜）。PCA结果表明，第一主成分解释的变异度为48.3%，第二主成分解释的变异度为27.0%；虚拟合成甘蓝型油菜比其亲本物种更接近于自然甘蓝型油菜，其中来源于甘蓝类群I的虚拟合成甘蓝型油菜与自然甘蓝型油菜的遗传特征最接近。

具体实施方式

本发明突出的实质性特点可从下述实施例得以体现，但并不对本发明作任何限制。

本实施例选择的异源多倍体物种为甘蓝型油菜（基因组为AACC），其亲本物种为甘蓝（基因组为CC）和白菜型油菜（基因组为AA），待测亲本物种为甘蓝。针对涉及十七种不同变异类型的39份甘蓝材料，利用4份不同遗传背景的白菜型油菜与其合成了156份虚拟甘蓝型油菜，通过与6份不同遗传背景的甘蓝型油菜进行遗传差异分析，发现甘蓝型油菜中C亚基因组很可能来源于栽培甘蓝及部分与其遗传关系较近的野生甘蓝，与前人研究结果相符。该实例证明了本发明方法（一种虚拟合成物种调查异源多倍体物种起源的方法）具有可行性。

实施例1：

虚拟合成甘蓝型油菜调查甘蓝型油菜中C亚基因组来源的方法，其步骤如下：

甘蓝和白菜型油菜（亲本物种）及甘蓝型油菜全基因组水平分子特征扫描：采用SSR和AFLP标记技术，对甘蓝的十七种变异类型、不同遗传背景的白菜型油菜，以及不同遗传背景的自然甘蓝型油菜进行全基因组水平分子标记分析；

虚拟合成甘蓝型油菜：根据甘蓝和白菜型油菜的标记类型，按照图2所示原则模拟出合成甘蓝型油菜的分子标记数据；

亲缘关系分析：利用NTSYS-PC (http://www.exetersoftware.com/cat/ntsyspc/ntsyspc.html下载)、STRUCTURE (http:// pritch. bsd. uchicago. edu/ structure. Html下载)、ARLEQUIN (http://cmpg.unibe.ch/software/arlequin3下载)等遗传关系分析软件和SAS数据统计分析软件(http://www.sas.com/下载)，分析来源于甘蓝不同变异类型的虚拟合成甘蓝型油菜与自然甘蓝型油菜的遗传差异，推导甘蓝型油菜中C亚基因组可能的来源。

本实施例中所用甘蓝材料是从IPK（the Institute of Plant Genetics and Crop Plant Research in Germany）、 UPM（the Universidad Politécnica de Madrid in Spain）和CGN（the Centre for Genetic Resources in the Netherlands）征集而来，白菜型油菜来源于西南大学种质资源库，甘蓝型油菜中的冬性材料和春性材料来自德国NPZ公司（Norddeutsche Pflanzenzucht Hans-Georg Lembke KG），半冬性材料来源于西南大学种质资源库。分子标记和遗传差异分析是本领域中所常采用的方法；总DNA提取、基于PCR的分子标记都有成熟的技术体系，其所涉及的试剂均可通过商业途径获得；使用的分析软件均是网站公开供免费使用的。

实施例2：亲本物种及自然甘蓝型油菜全基因组水平分子特征扫描

1）材料准备：

收集待测亲本的多种变异类型，具有不同遗传背景的另一亲本材料数份，以及具有不同遗传背景的自然异源多倍体数份。本实施例中的待测亲本物种为甘蓝，共39份，涉及十七种不同的变异类型；另一亲本物种为白菜型油菜，共4份，其中欧洲白菜和中国白菜各2份；自然甘蓝型油菜6份，其中冬性、春性和半冬性材料各2份。

2）提取材料全基因组DNA，具体步骤如下：

将2-3g幼叶用液氮磨成粉末，迅速转入15 mL离心管中；向管中加入4 mL 65 ℃预热的CTAB抽提液（2% CTAB，1.4M NaCl，0.02M EDTA，0.1M Tris-HCl）及80μL 的β-巯基乙醇，盖严后迅速摇匀，于65 ℃水浴60 min，每隔15min轻摇一次；室温下冷却后，向管中加入4 mL氯仿/异戊醇（v：v = 24：1），轻柔颠倒混匀，乳化20 min；之后12000 rpm离心10 min；转移上清至另一干净的离心管中，加入3 mL异丙醇，颠倒混匀，约1至2分钟后出现白色絮状沉淀，然后置于冰上静置20 min； 10000 rpm离心1min之后，倒掉管中液体，留下管底的DNA；向其中加入1 mL 75%乙醇，并全部转入1.5 mL离心管中；在75 %乙醇中漂洗沉淀DNA 5 min，然后于5000 rpm离心1min，倒弃乙醇，再重复本步骤一次；无水乙醇漂洗一次，5 min，弃乙醇，沉淀于37 ℃恒温箱中干燥2 h；之后用500 μL TE溶解DNA，并向其中加入10 μL RNAase，轻柔混匀后于37 ℃水浴锅中水浴1 h；之后向管中加入500 μL酚、氯仿、异戊醇的混合液（V_酚：V_氯仿：V_异戊醇 = 25：24：1），轻轻颠倒混匀使乳化20 min；12000 rpm离心20 min，吸取上清至一新的1.5mL离心管中，向上清中加入300μL已在-20 ℃下预冷20 min的异丙醇，冰浴沉淀30 min；12000 rpm、4 ℃条件下离心10 min后，倒弃液体；沉淀用500μL 的75％乙醇漂洗二次，每次5分钟；无水乙醇漂洗5分钟，弃乙醇，沉淀于37 ℃恒温箱中干燥1 h；之后向管中加入200 μL灭菌重蒸水，溶解DNA；完全溶解后检测浓度，并将浓度调至100 ng/μL。

3） AFLP标记：

将DNA调整到100 ng/μL后，首先用Mse I（Trul I）和PST I对总DNA进行酶切，按Vos等（1995）的报道设计接头和引物，由上海生工合成。接头退火后与酶切后的DNA进行连接，10倍稀释连接产物后，利用预扩增引物进行预扩增，10倍稀释预扩产物（Pre-A），作为选择性扩增的底物。选择性扩增的反应体系如下：

成分	体积（μL）
		ddH₂O	10.2
10 × Taq Buffer	2.0
		MgCl₂（25 mM）	1.5
P-primer (50 ng/μL)	0.6
		M-primer (50 ng/μL)	0.1
dNTPs (10 mM)	0.4
		Taqase (5 U/μL)	0.2
DNA (Pre-A)	5.0
		总体积	20.0

PCR程序为:

94 ℃ 30 s，[94 ℃ 30 s，65 ℃（-0.7 ℃/cycle）30 s，72 ℃ 1 min] ×13 cycles，[94 ℃ 30 s，56 ℃ 30 s，72 ℃ 1 min] ×23 cycles，72 ℃ 10 min

反应结束后加入10μL反应终止液（98%去离子甲酰胺，10mM EDTA（pH8.0），0.025%二甲苯氰，0.025%溴酚兰）于94 ℃变性3min后迅速于0℃冰水混合物中冷击，之后即用于电泳检测。电泳结束后进行数据统计，同一位点上，有带记为“1”，无带记为“0”。

4）SSR标记：

PCR反应体系如下：

成分	体积（μL）
		ddH₂O	6.38
10 × Taq Buffer	1.2
		MgCl₂（25 Mm）	1.0
dNTPs（10 Mm）	0.3
		F-Primer（50 ng/μL）	0.5
P-Primer（50 ng/μL）	0.5
		Taqase（5 U/μL）	0.12
DNA templete (50-100 ng/μL)	2.0
		总体积	12.0

PCR反应程序为：94 ℃ 5 min，[94 ℃ 45 s，55 ℃ 45 s，72 ℃ 1 min] ×35 cycles，72 ℃ 10 min

PCR反应结束后进行非变性聚丙烯酰胺凝胶电泳检测，数据统计方式为：同一多态性位点上，有带记为“1”，无带记为“0”。

结果：本实施例中11对AFLP选择性引物在49份材料中共检测出355个多态性位点，83对SSR特异引物扩增出464条多态性带，共计检测到819个清晰可读的多态性位点。

实施例3：虚拟合成甘蓝型油菜

根据合成异源多倍体过程中，亲本种具有的分子特征能在合成后代中得到体现的原则（图1），利用亲本物种的分子标记数据模拟出合成种的分子标记数据，具体模拟方式如下（图2）：同一位点上，双亲或亲本之一中为“1”的则在虚拟合成甘蓝型油菜中记为“1”，在双亲中都为“0”的则在虚拟合成甘蓝型油菜中记为“0”。本实施例中用39份甘蓝和4份白菜型油菜共模拟出39 × 4 = 156份虚拟合成甘蓝型油菜，共有819个多态性标记位点。

实施例4：亲缘关系分析

1）利用NTSYS-PC软件估算各亲本与自然甘蓝型油菜间，以及各虚拟合成甘蓝型油菜与自然甘蓝型油菜间的遗传距离。

2）利用多种数据分析系统，验证是否成功扣除白菜型油菜基因组对遗传差异的干扰，具体实施方式如下：

采用ARLEQUIN软件，按亲本物种进行分子标记方差分析(AMOVA)，检验各亲本物种（甘蓝和白菜型油菜）对虚拟合成甘蓝型油菜之间遗传差异的贡献率。结果：甘蓝亲本对虚拟合成甘蓝型油菜间遗传差异的贡献率为32.14%，白菜型油菜亲本的贡献率为16.36%，表明本实施例中虚拟合成甘蓝型油菜间的遗传差异主要来源于甘蓝亲本。

用SAS软件，对来源于同一亲本物种内不同个体的虚拟合成甘蓝型油菜之间的遗传距离进行方差分析（ANOVA）。本实施例结果：来源于4份不同的白菜型油菜的虚拟合成甘蓝型油菜间的遗传距离没有显著差异（F_3, 932= 1.10, P = 0.35），而来源于39份不同的甘蓝的虚拟合成甘蓝型油菜间的遗传距离存在显著差异（F_2, 933 = 567.8, P < 0.001）。

上述结果表明，本实施例中虚拟合成甘蓝型油菜间的差异主要来源于不同的甘蓝，而不同白菜型油菜分子特征的导入对其影响很小。因此，通过对自然甘蓝型油菜与虚拟合成甘蓝型油菜进行遗传差异比较，自然甘蓝型油菜中因白菜型油菜遗传成分（即A基因组）造成的差异可以被虚拟合成甘蓝型油菜中导入的白菜型油菜成分所扣除。

3）采用多种分析软件分析甘蓝型油菜中C亚基因组的来源，具体实施方式如下：

利用NTSYS-PC软件对亲本、虚拟合成甘蓝型油菜、自然甘蓝型油菜进行UPGMA聚类分析和主成分分析（PCA），调查甘蓝亲本的多种变异类型之间、甘蓝亲本与自然甘蓝型油菜之间、以及来源于不同甘蓝变异类型的虚拟合成甘蓝型油菜与自然甘蓝型油菜之间的遗传关系；利用STRUCTURE软件，分析甘蓝亲本的不同变异类型与自然甘蓝型油菜之间的血缘关系。结果：甘蓝亲本多种变异类型之间、亲本与自然甘蓝型油菜之间的遗传关系如图3和图4所示，涉及17种变异类型的39份甘蓝被划分为三个主要类群（甘蓝类群I，甘蓝类群II和甘蓝类群III），各类群间具有不同的血缘；而白菜型油菜和自然甘蓝型油菜则分别各自成群，其中甘蓝型油菜具有白菜型油菜以及部分与甘蓝类群I相同的血缘。亲本物种与自然甘蓝型油菜之间、虚拟合成甘蓝型油菜与自然甘蓝型油菜之间的遗传关系如图5，虚拟合成甘蓝型油菜比其亲本物种更接近于自然甘蓝型油菜，其中来源于甘蓝类群I的虚拟合成甘蓝型油菜与自然甘蓝型油菜的遗传特征最接近。

以甘蓝的不同变异类型为单位，用SAS软件对甘蓝与自然甘蓝型油菜、虚拟合成甘蓝型油菜与自然甘蓝型油菜间的遗传距离进行方差分析及多重比较。本实施例结果如表1，虚拟合成甘蓝型油菜与自然甘蓝型油菜间的遗传距离显著小于相应甘蓝亲本与自然甘蓝型油菜的遗传距离；来源于甘蓝类群I的虚拟甘蓝型油菜与自然甘蓝型油菜的遗传距离显著小于来源于甘蓝类群II和III的虚拟甘蓝型油菜与自然甘蓝型油菜间的遗传距离；来源于甘蓝类群I内不同变异类型的虚拟甘蓝型油菜与自然甘蓝型油菜间的遗传距离差异很小。

表1

Figure 2010105198222100002DEST_PATH_IMAGE001

本实施例中多种分析结果一致表明，甘蓝型油菜与甘蓝类群I（栽培甘蓝及部分野生甘蓝）中的甘蓝有很近的亲缘关系，甘蓝型油菜中的C亚基因组很可能来源于栽培甘蓝及部分与其遗传关系较近的野生甘蓝。

Claims

1.一种虚拟合成物种调查异源多倍体物种起源的方法，其步骤如下：

1）亲本物种及自然异源多倍体全基因组水平分子扫描：采用PCR技术对自然异源多倍体及其亲本物种进行全基因组水平的分子标记分析；

2）虚拟合成异源多倍体：根据合成异源多倍体过程中，亲本具有的分子特征能在合成后代中得到体现的原则，利用亲本物种的分子标记数据，模拟合成种的分子标记数据；

3)亲缘关系分析：利用NTSYS软件计算各亲本与自然异源多倍体，以及各虚拟合成异源多倍体与自然异源多倍体间的遗传距离，并对分子标记数据进行UPGMA聚类分析和主成分分析（PCA）；用STRUCTURE软件分析二倍体亲本的不同变异类型与自然异源多倍体之间的血缘关系；用ARLEQUIN软件，按亲本物种进行分子标记方差分析(AMOVA)，检验各亲本物种对虚拟合成异源多倍体遗传差异的贡献率；用SAS软件对来源于同一亲本物种内不同个体的虚拟合成异源多倍体之间的遗传距离进行方差分析（ANOVA），并对亲本与自然异源多倍体、虚拟合成异源多倍体与自然异源多倍体间的遗传距离进行方差分析及多重比较；通过以上分析揭示不同亲本种类型模拟合成的虚拟异源多倍体与自然异源多倍体相互间遗传差异的大小，根据遗传差异越小，该亲本类型与异源多倍体越为接近，则该类型越可能是异源多倍体中亚基因组的来源的原则，推导异源多倍体中亚基因组的来源。

2.根据权利要求1所述的虚拟合成物种调查异源多倍体物种起源的方法，其特征在于：所述其中步骤2）中，根据双亲或亲本之一中标记为“1”的则在虚拟合成异源多倍体中记为“1”，在双亲中都为“0”的则在虚拟合成异源多倍体中记为“0”的原则，据此利用双亲的分子标记数据模拟出合成种的分子标记数据。

3.根据权利要求1或2所述的虚拟合成物种调查异源多倍体物种起源的方法，其特征在于：步骤1）用于分子标记的PCR技术采用全基因组水平分子标记AFLP, SSR, RAPD, ISSR中的任意一种或多种。

4.根据权利要求1或2所述的虚拟合成物种调查异源多倍体物种起源的方法，其特征在于：步骤3）采用的数据分析软件是SAS、SPSS软件；所采用的遗传关系分析软件选择NT-SYS、STRUCTURE软件。