基于的ssr 遗传多样性分析

NOTE 基于种群聚类,特有等位基因

基本概念

1 多态性条带

个人理解就是ssr引物1 在位点 扩增出的 a b 2个条带(2倍体),ab长度可能相同也可能不相同。如果a条带 在所有样本DNA中的长度一样,就不是多态性条带。如果有不相同就属于多态性条带?还是在学习中慢慢理解吧,不能单独去理解这一个问题花费很多时间。多态性条带应该是属于显性标记中的内容,文献中说ssr也有多态性片段概念

2 ISSR

下为ISSR 检测流程,ISSR 引物不分F R, 引物即是F 又是 R。ISSR引物可在dna 两条链中结合。需要某段序列两段具有 反向且互补的重复单位,这样单条引物可扩增到片段。这样需要反向且互补的重复单位间的这段序列尽量在样本间不显示差异。所以才说ISSR是显性,单个ISSR引物在样本中只能扩增出一条带(理论情况),其实不是一条带,只要存在上述情况就可扩增


what-is-issr.png

3 TP-M13-SSR

参考文献为
Schuelke, Markus. "An economic method for the fluorescent labeling of PCR fragments." Nature biotechnology 18.2 (2000): 233-234.

主要流程如下,但个人感激此图做的并不准确,带有FAM荧光基团M13序列并不与E中的带有F向引物的扩增产物序列退火(红色方框中)。而是与其互补序列退火.
image.png

对PCR原理还不是很清晰,画图了解以下,红色为M13序列,蓝色为F R引物, 橙色阴影部分表示引物与序列退火,QQQQ为目标序列,NNNN为非目标序列

3.1 首先 带有M13接头序列的正向引物和反向引物,对目标序列退火
M13F R
生成以下扩增产物
image.png

3.2 反向引物扩增带有M13接头序列的正向引物的序列获得
M13接头序列的正向引物和反向引物扩增产物

3.3 带有荧光基团的M13序列作为正向引物和剩余反向引物 扩增上步骤得到初步的扩增产物,最后由带有荧光基团的M13正向引物扩增得到的产物可由仪器可检测到
image.png

4 聚丙烯酰胺 电泳图

看文献,聚丙烯酰胺 电泳图中SSR引物有很多条带,怎么判断哪些是 产物条带?,下面这个文献中的电泳图是属于与荧光检测保持一致的,即纯合位点为一条带,杂合位点2条带。
聚丙烯酰胺 电泳图,参考:利用荧光SSR 分析中国糜子的遗传多样性和群体遗传结构

5 遗传多样性参数

5.1 观测杂合度 (Ho); Observed Heterozygosity (Codom Data)

杂合个体数与样本数的比值
Ho

5.2 期望杂合度(He) Expected heterozygosity (Codom Data)

1减去所有等位基因频率平方和,但不明白此参数能表明什么,
He

5.3 Fixation Index (F)(Codom Data),

F 趋近0 说明位点under 随机杂交,符合哈温平衡?大量正值表明 近亲杂交(inbreeding)

或 没有检测到等位 基因,负值说明
GenAlex中的F

5.4 有效等位基因数

理想种群下的

Ne

5.5 基因流 Nm

Nm

流程

1 引物获取 筛选

从转录组或基因组数据中获得批量SSR引物,也可从以往论文中选取引物,进行引物筛选。引物设计使用 oligo7 或 primer 等。oligo是免费的,primer 需付费使用。

筛选引物可通过琼脂糖,聚丙烯酰胺,毛细管电泳检测引物扩增效果,最后筛选到特异性较强的引物进行基因组扩增。
20200901 数据全部完成

2 数据分析

本次实验采用 荧光 毛细管电泳。以软件的使用作为实验内容。软件的具体用法见文集中的软件使用说明,其实手册已经写的很详细了。
最终数据结果如下,会得到具体片段长度的大小

毛细管电泳结果

2.1 分析遗传多样性参数

共24个位置无数据
使用GenAlex分析遗传多样性参数,但是不太明白这些参数的含义,还担心在计算的过程中会遗漏数据,就跑2次,看下结果一样不
分不同地区比较不同种源的遗传多样性指数差异。目前以省划分最低分类

杂合度低 与 自花授粉有关?

2.2 计算PIC

GenAlex无法计算PIC,使用Powermarker

2.3 进化树分析

Powermarker计算遗传距离,UPGMA或NJ 法建树
DARwin 也有文章 使用这个软件 ,很多文章使用 DARwin吗,Powermarker 比较简单
得到树文件,使用itol美化。统计每个亚群下的种群个数。

2.4 PCA 分析

GenAlEx 中有此功能

2.5 AMOVA 分析

variance components of the populations, 即种群遗传分化

可根据地区,聚类,structure 的种群分类来做此分析, 并比较不同分类下的遗传差异,种群内外变异差异。
(1)GenAlex中AMOVA分析需要确定以下参数

AMOVA参数
对于共显性数据有以下三个选项,其中Codom - Microsat适合SSR数据,不知道选择哪个选项。
image.png

(2)其他参数
也不太清楚每个参数的含义,就先选择默认参数
image.png

(3) 结果
会得到类似以下的结果

文献中amova结果

聚类 与 structure 之间的某些种群用venn 图表示共有样本
GenAlEx 可用,

2.4 种群structure 分析

每个种群中个体根据 probability score 划分为 pure 或admixture,单独划为一个subgroup?
根据第一次划分的种群,可以继续划分sub-sub group
(1) 使用Structure Harvester获得最佳k
(2)Clumpp 重复抽样从多个run获得结果
(3) distruct 画图

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 203,456评论 5 477
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 85,370评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 150,337评论 0 337
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,583评论 1 273
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,596评论 5 365
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,572评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,936评论 3 395
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,595评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,850评论 1 297
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,601评论 2 321
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,685评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,371评论 4 318
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,951评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,934评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,167评论 1 259
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 43,636评论 2 349
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,411评论 2 342