关注公众号

关注公众号

手机扫码查看

手机查看

喜欢作者

打赏方式

微信支付微信支付
支付宝支付支付宝支付
×

水稻RGA基因组规模的鉴定与分析

2019.4.20

实验概要

本实验对粳稻蛋白质数据库进行BLASTP,获得所有粳稻的RGA序列,对RGA进行分类,分别统计了各类RGA在不同染色体上的分布情况,对两个亚种间RGA多态性进行了鉴定和开发。

实验步骤

1. 水稻RGA的搜索

搜集已报道的45个R基因,然后用它们对粳稻蛋白质数据库进行BLASTP (参数取E<10-10,最小长度为该基因的80%),获得所有粳稻的RGA序列。去除粳稻数据库中存在的重复克隆,建立一个粳稻RGA蛋白数据库。同时,根据每个BLASTP搜索中匹配最好的结果,得到这些粳稻RGA的核普酸序列。此外,去除粳稻数据库中存在由于的克隆重复造成的R基因的重复。

2. 水稻RGA的结构分类及其在染色体上的分布

利用Hmmer程序中的hmmsearch部分,采用前面建立的功能域列表,在新建立的数据库中搜索基因序列所包含的功能域。利用pepcoil程序分析序列中CC结构的可能性,数值大于90%的认为具备该结构。运用TM-HMMer分析TM转膜功能域。依据功能域分布的情况,对RGA进行分类,分别统计各类RGA在不同染色体上的分布情况。

3. 两个亚种间RGA多态性的鉴定和开发

将所有粳稻的RGA序列分别与粕稻基因组数据库进行TBLASTN联配,以确定釉稻中对应的等位基因。为消除非等位联配,在TBLASTN搜索中采用了严格的判别标准,将E值设为10-20。对初筛到的粕、粳稻RGA等位基因,进一步用sim4程序进行联配分析,去除匹配率≤85%且两条序列同时间断200bp以上的结果。接着运用diffseq程序分析SNP和InDel在RGA的基因组序列中的分布情况,将存在InDel的作为候选的RGA标记。以粳稻的基因组序列为模板,在INDEL位置的两侧各取100bp的序列,利用ePrimer3程序在取得的序列的两端进行引物设计。ePrimer3程序一般给出五对候选引物,我们取其设计最合理的一条(考虑到目标片段最大长度不大于1000bp),并进行电子PCR(e-PCR)验证。对得到的水稻RGA标记进行命名,规则为以OSR开头,后跟4个数字,例如:OSR0255。上述步骤主要通过编写pert脚本程序来实现。


推荐
热点排行
一周推荐
关闭