使用FusionMap检测融合基因

欢迎关注微信公众号《生信修炼手册》!

融合基因是指两个或者多个基因联合起来,共同转录形成一个转录本,融合基因可以作为某些疾病的特异分子标记,常见的有以下几种

  1. bcr/abl融合基因存在于95%以上的慢性粒细胞白血病患者中;
  2. AML1/ETO融合基因主要见于急性粒细胞白血病部分分化型患者中;
  3. CBFβ/MYH11融合基因是M4Eo型白血病的分子标志;
  4. PML/RARα融合基因是急性早幼粒细胞白血病(APL)的分子标志

利用RNA_seq的数据,我们可以检测融合转录本,从而识别融合基因。 fusionMap就是这样的一款工具,官网如下

http://www.arrayserver.com/wiki/index.php?title=FusionMap

通过两种方式来检测融合基因

  • 对于没有比对上基因组的序列,即unmapped reads, 通过识别Fusion junction-spanning reads 来识别融合基因;这部分reads覆盖到了融合基因的连接点,在连接点两侧的序列分别比对到构成融合的多个基因上
  • 对于比对上基因组的reads, 通过识别 Inter-transcript read pairs来识别融合基因,这部分reads虽然没有直接覆盖连接点,但是其R1端和R2端分别比对到不同的基因
    示意图如下
  • 4
    点赞
  • 7
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
可以使用 MATLAB 中的线性选择方法来融合两种运动目标检测的结果,以下是一些具体步骤: 1. 准备两种不同的运动目标检测算法的结果。每种算法都应该生成一幅显著图像,用于表示每个像素的显著性程度。 2. 将两种显著图像进行归一化,使它们的像素值在[0,1]范围内。可以使用 MATLAB 中的函数“mat2gray”来进行归一化。 3. 对两种显著图像进行线性组合,得到融合后的显著图像。可以使用以下公式进行线性组合: `fusionMap = a * saliencyMap1 + (1 - a) * saliencyMap2` 其中,`a`是一个权重因子,表示对第一种显著图像的重要程度。可以根据实际情况调整`a`的值。 4. 对融合后的显著图像进行阈值处理,得到二值化的显著图像。可以使用 MATLAB 中的函数“imbinarize”来进行阈值处理。 以下是一个简单的示例代码,演示如何使用线性选择方法融合两种运动目标检测的结果: ```matlab % 加载两种显著图像 saliencyMap1 = imread('saliencyMap1.jpg'); saliencyMap2 = imread('saliencyMap2.jpg'); % 归一化显著图像 saliencyMap1 = mat2gray(saliencyMap1); saliencyMap2 = mat2gray(saliencyMap2); % 线性组合显著图像 a = 0.8; fusionMap = a * saliencyMap1 + (1 - a) * saliencyMap2; % 阈值处理融合后的显著图像 threshold = graythresh(fusionMap); binaryMap = imbinarize(fusionMap, threshold); % 显示结果 figure(1); imshow(saliencyMap1); title('Saliency map 1'); figure(2); imshow(saliencyMap2); title('Saliency map 2'); figure(3); imshow(fusionMap); title(sprintf('Fusion map (a=%.2f)', a)); figure(4); imshow(binaryMap); title('Binary map'); ``` 需要注意的是,线性选择方法只是一种简单的融合方法,可能无法充分利用两种检测算法的优势。可以考虑使用更复杂的融合方法,例如基于深度学习的融合方法,来提高检测性能。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值