SSR基因分型在遗传多样性研究中的应用进展

发稿时间：2018-08-07来源：天昊生物

 

SSR（简单重复序列，或称微卫星）作为一类遗传分子标记，在种质资源的遗传多样性分析中扮演着十分重要的角色。近期发表在BMC Plant Biology上的一篇文章，就利用SSR基因分型结果系统解析了地中海盆地和中亚地区栽培葡萄和野生葡萄的遗传多样性情况。

 

 

发表期刊：BMC Plant Biology 发表时间：2018-6-27 影响因子：3.93

 

摘要一览

背景

高加索和中国之间的山区被认为是葡萄驯化的中心。尽管中亚在葡萄种植历史上非常重要，但与地中海盆地周围的野生和栽培葡萄相比，关于该地区葡萄遗传变异程度和分布的信息有限。本研究的主要目的是调查高加索、中亚和地中海盆地野生和栽培葡萄种质的遗传多样性和相互关系，了解基因流动、可能的驯化事件和适应性渗入。

结果

利用20个SSR标记对地中海盆地周围和中亚收集的1378种野生和栽培葡萄进行了基因分型。遗传数据经过分析(聚类分析、主坐标分析和群体结构分析)确定了群体分组情况，结果还通过Nei’s遗传距离、成对F_ST分析和分配检验得到验证。所有这些分析确定了三个遗传群体: G1，来自克罗地亚、法国、意大利和西班牙的野生资源；G2，来自亚美尼亚、阿塞拜疆和格鲁吉亚的野生资源；G3，来自西班牙、法国、意大利、格鲁吉亚、伊朗、巴基斯坦和土库曼斯坦的品种，其中包括来自格鲁吉亚和克罗地亚的少量野生资源。乔治亚的野生资源与同一地区种植的葡萄聚集在一起，但也与西欧聚集在一起，支持乔治亚成为葡萄驯化的古老中心。此外，聚类分析表明，西欧野生葡萄与来自同一地区的栽培葡萄聚在一起，表明栽培的西方proles occidentalis区比来自东欧的野生葡萄对葡萄酒葡萄的早期发展贡献更大。

结论

通过对基因型遗传关系的分析，为地中海盆地和中亚野生和栽培品种间遗传关系提供的新的证据。遗传结构表明有相当多的基因流动，这限制了这两个亚种之间的分化。结果还表明，在野生葡萄被驯化的地区，有着混合祖先葡萄的出现。

 

 

实验方法

用ABI 3500进行SSR基因分型，GeneMapper 4.0对等位基因条带判读。

 

数据分析方法

遗传多样性

对群体间和所有群体的遗传多样性进行估计。用标准化SSR基因型数据确定不同等位基因数(Na)、有效等位基因数(Ne)、Shannon’s指数(I)、观测杂合度(Ho)和预期杂合度(He)。这些参数由GenAlEx软件完成。每个基因座F-统计(F_IS，F_IT，F_ST) 以及每个群体的F_IS值由FSTAT和Arlequin软件完成。每个群体的等位基因丰富度(AR)和特有等位基因丰富度(PAR)是用rarefaction估计，该方法弥补了HP-Rear1.1软件对群体间样本量(即稀有等位基因丰富度)的差异。使用GENEPOP 3.4软件的特有等位基因法，估计了12个葡萄群体以及两个亚种之间的世代迁移有效数(Nm)。

遗传关系与分化

使用R语言Poppr包构建NJ系统进化树，距离矩阵为Nei’s距离，无根树由R语言ape包实现。为了检验树中的层次结构所代表实际距离的程度, 使用R语言计算获得共表性相关系数(CCC)。Hclust R函数用来校正NJ层次聚类。为了阐明各地理群体内部和之间的遗传关系, 对多位点微卫星数据进行了主坐标分析 (PCoA), 然后利用R语言adegenet包排列到地理分组中。利用GenAlEx软件中的成对Nei’s遗传距离和成对 F_ST 进行聚类验证和多元分析。最后,用Arlequin软件对种群和遗传分组中遗传变异的划分进行AMOVA分析。

群体结构分析

利用STRUCTURE软件对SSR分型数据进行基于贝叶斯模型的聚类分析, 确定遗传分组的最佳数目。STRUCTURE根据基因型数据将个体分配到多个群中(K), 从而最大限度地减少与哈迪-温伯格和连锁不平衡的偏差。该程序使用马尔可夫链估计P (X|K)，用以判断符合K簇假设的后验概率。分析根据隶属度系数(Q值)将个体分配给每个K簇, 得到统一性最好的K值。STRUCTURE设置为忽略群体信息, 并使用相关等位基因频率的admixture模型, 因为它被认为是微小种群结构的最佳选择。K变化从2到10, 并进行了20次复制运行, 以量化变化的可能性，在整个运行中最大可能性 (Ln P(D)) 的 K 值通常被推断为最可能的簇数。然而，对K的解释应该小心对待, 因为它只是提供一个临时近似，有时真正的和精细的群体结构可能被STRUCTURE软件错过。因此, 本研究还使用了一个特殊的ΔK来选择统计的最佳簇数 (K)。

 

 

研究结果

表1、根据地理来源对1378种栽培和野生资葡萄的分类列表

 

  表2、为1378种不同葡萄基因型多样性指数计算结果 

 

 

表3、野生和栽培葡萄群体遗传多样性估计

 

 

图1、1378种栽培葡萄和野生葡萄NJ聚类图

 

 

 

图2、1378种栽培葡萄和野生葡萄PCoA图

 

 

 

图3、1378种栽培葡萄和野生葡萄群体结构（K = 3）分析

 

   

表4、野生和栽培葡萄群体内成对Nei’s遗传距离(对角线以下)和F_ST值(对角线以上)估计

 

 

 

 

 

结论

 1、本研究对葡萄sativa和sylvestris两个亚种进行了系统分析，SSR基因分型结果表现出高度的多态性和杂合性，说明葡萄亚种内部和亚种之间都存在显著的多样性(表2和3)。

• 来自佐治亚地区的sativa和sylvestris具有最高的有效等位基因数量和等位基因丰富度，这表明该地区是酿酒葡萄多样性的中心。
•  对Ho、He值及F_IS的比较结果表明，来自同一野生种群的个体之间存在高度的遗传关系，特别是在西班牙和格鲁吉亚的种群中，这种交配会影响个体和群体的动态，增加近亲交配几率。然而，在随机交配群体中，一些野生群体的F_IS值接近于零，这些相反的结果可以用样本大小导致的亚群中等位基因的随机遗传漂移来解释。
• 与其他研究类似，本研究也发现sylvestris样本中观察到的多样性水平的降低。由于森林砍伐和城市化，世界许多地方的sylvestris品种被认为濒临灭绝，支离破碎。人为的和自然的地理障碍也可能导致野生群体在他们的自然栖息地被隔离，并可能导致显著的近亲繁殖，减少不同地理群体内部和之间的基因流动，从而降低杂合度。
•  除了意大利品种外，栽培品种的F_IS值接近于零，表明为随机交配。意大利群体的负F_IS值表明杂合子过多，但在统计学上没有显著性(表3)。大多数栽培群体中纯合子的缺失表明，它们是由具有不同群体统计学(奠基者效应、瓶颈、扩散)和选择历史的种质组成的。种质收集通常是基因型的混合物。因此，这些集合中的地理群体表现出相对较高的分化水平，导致杂合度高于预期水平。

 

2、聚类分析和PCoA分析表明，亚种内部和亚种之间的遗传结构和分化具有显著差异(图1和图2)。

• 两种分析都发现西欧野生葡萄和高加索地区野生样品之间有明显的区别，而法国和西班牙野生葡萄关系密切，有着密切的遗传关系。 
• 西班牙野生葡萄表现出等级分化，表明相邻种群间的基因流动导致了种群结构的垫脚石模型。或者，等级差异可能是不同地理区域气候差异的结果。
• 克罗地亚的sylvestris与欧洲的sylvestris个体相关，形成了一个基础姐妹群，表明它们有着共同的基因库。包括亚美尼亚、阿塞拜疆和格鲁吉亚在内的外高加索野生葡萄藤组成了一个独特的亚群。
• 在sativa中，观察到来自格鲁吉亚的两组不同的品种，一组是意大利、法国和西班牙品种的姐妹品种(图1)，另一组与意大利sativa和格鲁吉亚sylvestris亚组密切相关。这一结果可能表明，中亚和高加索地区的第一批驯化品种在西欧资源中留下了遗传足迹。这种遗传亲缘关系也可能是地中海地区早期育种计划的反映，在那里，具有优良果实的sylvestris或杂交野生藤本植物被用于与驯化品系的杂交。
• PCoA描绘的总体分化模式非常类似于NJ聚类分析。来自格鲁吉亚和亚美尼亚的sylvestris品种中的集群重叠，并与来自格鲁吉亚、巴基斯坦和土库曼斯坦的栽培形式密切相关。

3、STRUCTURE分析支持主要群体之间的差异，而一些群体之间的细微差异，尤其是那些混合群体之间的差异并不明显(图3)。

• 遗传结构的贝叶斯推断表明，这两个亚种之间有相当大的基因流动，并有适度的分化。这些结果表明酿酒葡萄的种植和酿酒促进了野生葡萄的驯化、新品种的创造和葡萄历史早期生长技术的进步。
• 本研究结果表明西欧sylvestris向一些当前西欧品种存在着渗入现象。
• 对高度濒危的栽培葡萄树祖先作为抗病育种基因资源的潜力提供了的支持。鉴于欧亚、北非和亚洲野生葡萄种质和种质在很大程度上还未被开发，它们的鉴定、保存以及生物和非生物抗性和浆果品质的特征对于葡萄酒和葡萄产业的未来非常重要。

关于天昊：

天昊生物拥有多种SSR检测平台及SSRseq^TM等专利技术，可以根据客户项目需求，提供不同数量样本和位点的高性价比微卫星检测服务，欢迎联系我们具体咨询！邮箱：techsupport@geneskies.com 电话：400-065-6886