 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 7 | 1/1 | 返回列表 |
| 查看: 4674 | 回復(fù): 6 | |||
yuehedou木蟲 (小有名氣)
|
[求助]
miRNA差異表達(dá)分析時(shí)reads數(shù)目的校準(zhǔn)（TPM）問題
|
|
在做miRNA的差異表達(dá)分析時(shí)，很多文章里都要使用TPM（Transcripts per million，公式為：單一miRNA reads數(shù)×10^6/總reads數(shù)）來校準(zhǔn)miRNA的表達(dá)量，這樣做是為了去除pcr擴(kuò)增文庫時(shí)的偏倚嗎？ 然后，如果使用TPM進(jìn)行校準(zhǔn)，那也會(huì)產(chǎn)生另外的差異啊我覺得，舉個(gè)極端的例子：樣品A中實(shí)際總共表達(dá)了10個(gè)reads，其中miR1表達(dá)了2個(gè)；樣品B總共表達(dá)了100個(gè)reads，其中miR1表達(dá)了20個(gè)，——實(shí)際上在兩個(gè)樣品中miR1是顯著差異表達(dá)的，但經(jīng)TPM校準(zhǔn)后，兩個(gè)樣品中其表達(dá)量就都成了20 000，就一樣了么，那么校準(zhǔn)豈不是產(chǎn)生了原來沒有的錯(cuò)誤？？ 我覺得miRNA是短序列，不存在轉(zhuǎn)錄組中序列片段長短的問題，那么直接使用reads的count作為其相應(yīng)miRNA的表達(dá)量豈不是更好？ 歡迎討論！ |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 miRNA測序數(shù)據(jù)分析 |
» 猜你喜歡
不合理蛙科研實(shí)驗(yàn)之小鼠實(shí)驗(yàn)：嚴(yán)謹(jǐn)設(shè)計(jì)，解析生命機(jī)制的重要載體
已經(jīng)有0人回復(fù)
-
不合理蛙科研實(shí)驗(yàn)之重金屬檢測：精準(zhǔn)篩查，守護(hù)健康與環(huán)境的防線
已經(jīng)有0人回復(fù)
-
化學(xué)工程及工業(yè)化學(xué)論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有269人回復(fù)
-
不合理蛙科研實(shí)驗(yàn)之“蛙測重金屬我背鍋三千”
已經(jīng)有0人回復(fù)
-
不合理蛙科研實(shí)驗(yàn)之“鼠逃三次我發(fā)三篇SCI”
已經(jīng)有0人回復(fù)
-
納米粒度分析
已經(jīng)有3人回復(fù)
-
林科院林化所招收材料與化工專業(yè)碩士，坐標(biāo)南京
已經(jīng)有0人回復(fù)
-
留學(xué)--博士招生
已經(jīng)有16人回復(fù)
-
化學(xué)工程，本211，求調(diào)劑，ab區(qū)不限
已經(jīng)有4人回復(fù)
-
邵陽學(xué)院食品與化學(xué)工程學(xué)院碩士調(diào)劑
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- Illumina高通量RNA測序數(shù)據(jù)分析中不同標(biāo)準(zhǔn)化方法的綜合比較
已經(jīng)有6人回復(fù)
myibm369
木蟲 (正式寫手)
~~~~~~~~~~~~~~~~~
- 應(yīng)助: 17 (小學(xué)生)
- 金幣: 3127.1
- 紅花: 2
- 帖子: 457
- 在線: 113.8小時(shí)
- 蟲號(hào): 86473
- 注冊(cè): 2005-08-15
- 性別: GG
- 專業(yè): 生物信息學(xué)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
yuehedou: 金幣+18, ★★★★★最佳答案, 非常給力！ 2014-05-06 15:10:14
感謝參與，應(yīng)助指數(shù) +1
yuehedou: 金幣+18, ★★★★★最佳答案, 非常給力！ 2014-05-06 15:10:14
|
這其實(shí)是個(gè)normalization的過程，因?yàn)榧偃缒阌卸鄠�(gè)樣品（》=2），測序得到的total mapping reads的數(shù)目是不一樣的，大家必須normalize到同一個(gè)水平才方便比較。 你說的實(shí)際應(yīng)該是指細(xì)胞內(nèi)真正表達(dá)的數(shù)量吧，但是實(shí)際上通過測序你是拿不到真正意義上的絕對(duì)表達(dá)量，只是一個(gè)相對(duì)的表達(dá)量。我們?cè)谶M(jìn)行2個(gè)或者多個(gè)樣品genes或者miRNA表達(dá)差異分析比較的時(shí)候，是默認(rèn)總體RNA表達(dá)的一樣多的。簡單來說就是你例子中sample A和B的RNA總量是一致的，不過A測序測到了10個(gè)reads，B測序測的比較多，拿到了100個(gè)reads，如果A也達(dá)到B的測序深度，其實(shí)miR1也就是20個(gè)reads，表達(dá)量是一致的。 不知道你明白否 |
2樓2014-05-06 09:24:15
yuehedou
木蟲 (小有名氣)
- 應(yīng)助: 26 (小學(xué)生)
- 金幣: 2171.4
- 紅花: 2
- 帖子: 215
- 在線: 190.5小時(shí)
- 蟲號(hào): 1996646
- 注冊(cè): 2012-09-12
- 性別: GG
- 專業(yè): 生物信息學(xué)
|
我覺得total mapping reads的數(shù)目不一樣才是正常的，因?yàn)槭遣煌瑯悠贰?fieldset> 你說的實(shí)際應(yīng)該是指細(xì)胞內(nèi)真正表達(dá)的數(shù)量吧…… 我說的是實(shí)際表達(dá)量。我認(rèn)為只要能控制建庫均勻度、測序深度，那得到的數(shù)據(jù)反映的應(yīng)該就是實(shí)際表達(dá)豐度，即使這種反映是將原來的表達(dá)量以一種集體相對(duì)的方式表現(xiàn)出來，例如將實(shí)際的表達(dá)量整體放大。因?yàn)閷?shí)際總體和單個(gè)的表達(dá)量就不一致，在測序深度同樣為10的情況下，樣品A中miR1表達(dá)了1，放大10倍就測得了10；B中miR1表達(dá)了10，放大10倍就成了100. 那我們后邊又進(jìn)行校準(zhǔn)，就人為地抹掉了這種差異——實(shí)際上它們卻是有顯著差異的。 |
3樓2014-05-06 11:53:03
myibm369
木蟲 (正式寫手)
~~~~~~~~~~~~~~~~~
- 應(yīng)助: 17 (小學(xué)生)
- 金幣: 3127.1
- 紅花: 2
- 帖子: 457
- 在線: 113.8小時(shí)
- 蟲號(hào): 86473
- 注冊(cè): 2005-08-15
- 性別: GG
- 專業(yè): 生物信息學(xué)
|
其實(shí)你假設(shè)下就知道了，上機(jī)測序的total miRNA量 A為10ng其中miR1有1ng，B為10ng其中miRNA1有0.1ng，那么你最后測序得到的結(jié)果比例還是這樣的，A中miR1為10%，B中miR1為1%，都normalization到100，那么A的miR1是10，B的是1。 <quote>在測序深度同樣為10的情況下，樣品A中miR1表達(dá)了1，放大10倍就測得了10；B中miR1表達(dá)了10，放大10倍就成了100. 那我們后邊又進(jìn)行校準(zhǔn)，就人為地抹掉了這種差異——實(shí)際上它們卻是有顯著差異的。</quote> 校準(zhǔn)還是10：100哦，是有差異的啊。 |
4樓2014-05-06 12:13:49
yuehedou
木蟲 (小有名氣)
- 應(yīng)助: 26 (小學(xué)生)
- 金幣: 2171.4
- 紅花: 2
- 帖子: 215
- 在線: 190.5小時(shí)
- 蟲號(hào): 1996646
- 注冊(cè): 2012-09-12
- 性別: GG
- 專業(yè): 生物信息學(xué)
5樓2014-05-06 15:09:26
wangjiying
鐵蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 970.1
- 紅花: 2
- 帖子: 68
- 在線: 66.5小時(shí)
- 蟲號(hào): 779955
- 注冊(cè): 2009-05-26
- 專業(yè): 畜牧學(xué)
6樓2015-08-19 00:44:37
yuehedou
木蟲 (小有名氣)
- 應(yīng)助: 26 (小學(xué)生)
- 金幣: 2171.4
- 紅花: 2
- 帖子: 215
- 在線: 190.5小時(shí)
- 蟲號(hào): 1996646
- 注冊(cè): 2012-09-12
- 性別: GG
- 專業(yè): 生物信息學(xué)
7樓2015-08-20 16:05:04
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 yuehedou 的主題更新
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +7 | 爭取九點(diǎn)睡 2026-03-28 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 071000生物學(xué)求調(diào)劑，初試成績343 +7 | 小小甜面團(tuán) 2026-03-25 | 7/350 |
|
|
|
[考研] 283求調(diào)劑 +3 | A child 2026-03-28 | 3/150 |
|
|
|
[考研] 085600，材料與化工321分求調(diào)劑 +9 | 大饞小子 2026-03-28 | 9/450 |
|
|
|
[考研] 322求調(diào)劑 +5 | 舊吢 2026-03-24 | 5/250 |
|
|
|
[考研] 328求調(diào)劑 +7 | 嗯滴的基本都 2026-03-27 | 7/350 |
|
|
|
[考研] 266求調(diào)劑 +11 | 陽陽哇塞 2026-03-27 | 12/600 |
|
|
|
[考研]
 一志愿 西北大學(xué) 總分282 英語一62 求調(diào)劑
+7
|
18419759900 2026-03-25 | 8/400 |
|
|
|
[考研] 一志愿吉大071010，316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
|
[考研] 0703化學(xué)338求調(diào)劑！ +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
|
[考研] 276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷 +10 | 材料學(xué)257求調(diào)劑 2026-03-23 | 11/550 |
|
|
|
[考研] 359求調(diào)劑 +4 | 王了個(gè)楠 2026-03-25 | 4/200 |
|
|
|
[考研] 284求調(diào)劑 +11 | junqihahaha 2026-03-26 | 12/600 |
|
|
|
[考研] 調(diào)劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
|
[考研] 352求調(diào)劑 +4 | 大米飯！ 2026-03-22 | 4/200 |
|
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7 | jhybd 2026-03-23 | 12/600 |
|
|
|
[考研] 070300，一志愿北航320求調(diào)劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|
|
|
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|