博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
如何提取一个转录本的3'UTR区域的序列
阅读量:7096 次
发布时间:2019-06-28

本文共 770 字,大约阅读时间需要 2 分钟。

在做microRNA 和 mRNA 相互作用预测的时候,大家都知道microRNA 作用的靶点是位于mRNA 的3'UTR取,所以只需要提取mRNA 对应的3'UTR 区的序列去做分析即可;

那么如何提取一个mRNA的3'UTR区呢?

在UCSC数据库中,提供了3'UTR区序列的下载,以人类hg19为例, 利用table browser 浏览器选择对应的序列

链接:http://genome.ucsc.edu/cgi-bin/hgTables

按下图所示进行选择

点击get output 按钮,在弹出的页面选择 genomic

点击 submit 按钮,在弹出的页面勾选需要的区域,这里我们只选择 3'UTR区域

然后点击下方的get sequence 按钮,在浏览器中保存文件即可。

UCSC为我们提供了自动化的下载转录本特定区域的功能,如果我们自己来完成这件事,又该如何去做?

其实只需要两步:

1)第一步,确定每个mRNA的3' UTR区在基因组上的位置;

2) 第二步,根据基因组上的位置,从基因组上提取对应的序列就可以了;

如何定义一个转录本的3’UTR区呢,我们看UCSC是如何定义的,

以转录本NM_033487 为例,从UCSC下载的该转录本的序列为

利用NCBI的 nucleotide 数据库检索该转录本序列,链接如下 https://www.ncbi.nlm.nih.gov/nuccore/NM_033487

在对应的页面可以看到 该转录本的ploy A 尾开始的位置为2824;

在对应的序列中,可以看出poly A 尾之前的5bp的序列为aggaa, 和 UCSC对应的3’UTR区是一致的

从UCSC下载的3’UTR序列的长度为523bp,对应的的基因组位置为 chr1:1570603-1571125;

 

转载地址:http://vrxql.baihongyu.com/

你可能感兴趣的文章
使用Jmeter+Maven+Jenkins实现接口自动化测试
查看>>
用Class ID做网络分流
查看>>
nginx配置详解
查看>>
shell脚本之循环结构&函数的使用(二)
查看>>
Java基础学习总结(23)——GUI编程
查看>>
navigationItem标题以及左右button的创建...
查看>>
RabbitMQ学习总结(7)——Spring整合RabbitMQ实例
查看>>
zookeeper初探三 java客户端连接
查看>>
ImportError: Settings cannot be imported 错误处理
查看>>
宝玉论“后生可畏”,存于此以自勉
查看>>
linux mkdir命令
查看>>
可变参数列表
查看>>
网络综合实训学习资料
查看>>
iOS9 App Thinning(应用瘦身)功能介绍
查看>>
LVS(DR) + Keepalive实现负载均衡
查看>>
JavaScript基础教程(一)
查看>>
关于容灾的那些事儿
查看>>
查询?还是计算?这不再是个问题!(二)
查看>>
SQL Server Alwayson读写分离配置
查看>>
【源资讯 第41期】前端项目一言不合就分叉;工程师一言不合就删库??
查看>>