 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 7 | 1/1 | 返回列表 |
| 查看: 5493 | 回復(fù): 6 | ||||||
QQ894064647新蟲 (初入文壇)
|
[交流]
QIIME--16s測(cè)序生物信息培訓(xùn)課程(原創(chuàng)) 已有6人參與
|
|
QIIME--16s測(cè)序生物信息培訓(xùn)課程(原創(chuàng)) 名詞解釋 16S rRNA 16S rRNA 基因是編碼原核生物核糖體小亞基的基因，長(zhǎng)度約為1542bp，其分子大小適中，突變率小，是細(xì)菌系統(tǒng)分類學(xué)研究中最常用和最有用的標(biāo)志。其序列包含可變區(qū)（Variable region）和保守區(qū)（constant region）�？勺儏^(qū)序列因細(xì)菌而異，且變異程度與細(xì)菌的系統(tǒng)發(fā)育密切相關(guān)，通過(guò)提取環(huán)境樣品的DNA，并擴(kuò)增和檢測(cè)其中的一段或幾段區(qū)域，可以了解環(huán)境樣品(土壤、腸道、海洋)中細(xì)菌的種類和豐度。 OTU 如果序列之間，比如不同的 16S rRNA序列的相似性高于97%就可以把它定義為一個(gè)OTU，每個(gè)OTU對(duì)應(yīng)于一個(gè)不同的16S rRNA序列，也就是每個(gè)OTU對(duì)應(yīng)于一個(gè)不同的細(xì)菌（微生物）種。通過(guò)OTU分析，就可以知道樣品中的微生物多樣性和不同微生物的豐度。 測(cè)序區(qū)域 由于16s rDNA較長(zhǎng)（1.5kb），考慮到測(cè)序平臺(tái)的讀長(zhǎng)限制，我們只能對(duì)可變區(qū)進(jìn)行測(cè)序。16s rDNA包含有9個(gè)可變區(qū)，分別是v1-v9。一般我們對(duì)v3，v3-v4和v6可變區(qū)域進(jìn)行擴(kuò)增和測(cè)序。 主要內(nèi)容 QIIME安裝 測(cè)序數(shù)據(jù)處理 OTU分析與注釋 Venn圖繪制 Heatmap圖繪制 樣品復(fù)雜度分析 (Alpha diversity) 多樣品比較分析 (Beta diversity) 系統(tǒng)樹的構(gòu)建 其他分析 詳細(xì)內(nèi)容 1.Qiime的安裝(http://qiime.org/) 在http://qiime.org/官網(wǎng)下載VirtualBox和QIIME虛擬機(jī)，然后安裝VirtualBox，導(dǎo)入QIIME虛擬機(jī)，具體圖示如下。  我的電腦為什么安裝不了？那是因?yàn)槟碾娔X可能沒(méi)有開通VT(virtualization technology)，趕緊進(jìn)入BIOS開通吧。 2.Linux常見(jiàn)命令 大家如果沒(méi)有用過(guò)linux的話，請(qǐng)學(xué)會(huì)一下幾個(gè)常用的linux命令吧。如果大家想對(duì)我們的課程有詳細(xì)的詳解的話，大家可以在掏寶搜索“16s測(cè)序分析”就能找到我們的課程了。 cd 切換到目錄 pwd 顯示當(dāng)前目錄完整路徑 mkdir/touch 新建目錄/文件 ls 查看文件與目錄 cp 復(fù)制文件 mv 移動(dòng)文件、目錄或更名 rm 刪除文件或目錄 zip/tar 解壓縮文件 history 查看歷史命令 wget 命令行的下載工具 shutdown/reboot 關(guān)機(jī)/重啟 3.FASTQ數(shù)據(jù)下載 主要有以下兩個(gè)網(wǎng)站下載原始的FASTQ： EBI (http://www.ebi.ac.uk/) SRA (http://www.ncbi.nlm.nih.gov/sra/) 4.樣品信息文件(map.txt) #SampleID BarcodeSequence LinkerPrimerSequence SampleType Description SRR1370913 AGCTGACTAGTC GTGCCAGCMGCCGCGGTAA tumor SRR1370913 SRR1370914 AGCTGACTAGTG GTGCCAGCMGCCGCGGTAA tumor SRR1370914 SRR1370915 AGCTGACTAGTA GTGCCAGCMGCCGCGGTAA normal SRR1370915 SRR1370920 AGCTGACTAGTT GTGCCAGCMGCCGCGGTAA normal SRR1370920 第一列必須是#SampleID，為可區(qū)分的數(shù)字，字母或 。每一行的這個(gè)值應(yīng)該是唯一的。 第二列必須為BarcodeSequence。每一行的這個(gè)值應(yīng)該是唯一的。 第三列必須為L(zhǎng)inkerPrimerSequence，為擴(kuò)增樣品的引物。 后面的列可以根據(jù)樣品的特點(diǎn)加以描述，但是每一列必須包含至少兩個(gè)值，就是指定樣品信息，就可以在這個(gè)位置設(shè)置了。 最后一列必須是Description。每一個(gè)樣品不一樣的地方，必須不一樣額。 5.join_paired_ends.py (fastq-join) 6.質(zhì)量過(guò)濾 (split_libraries_fastq.py) convert to fasta Read too short after quality truncation N characters 7.OTU分析與注釋--pick_open_reference_otus.py pick_otus.py------從fasta文件中提取OTUs，提取方法有cd-hit，blast，Mothur，usearch。 make_otu_table.py-------將提取出的OTUs制作成out_table，即.biom文件。 assign_taxonomy.py------使用reference中的reference_seqs和 taxonomy文件對(duì)序列分類。 biom add-metadata-------將taxa數(shù)據(jù)加入到.biom文件。 8.Venn圖繪制 使用VennDiagram package可以輕松獲得和高水平文章上類似的Venn（維恩）圖，代碼如下： 9.Heatmap圖 10.core_diversity_analyses.py（Alpha和Beta分析） single_rarefaction.py------Rarify the OTU table beta_diversity.py-------beta多樣性 principal_coordinates.py------PCoA圖形化 alpha_diversity.py-------Alpha多樣性 make_rarefaction_plots.py-------稀釋曲線 summarize_taxa.py------Summarize Taxonomy plot_taxa_summary.py------Taxa summary plots Alpha多樣性是指一個(gè)特定區(qū)域或者生態(tài)系統(tǒng)內(nèi)的多樣性，常用的度量指標(biāo)有Chao1 豐富度估計(jì)量（Chao1 richness estimator） 、香農(nóng)多樣性指數(shù)（Shannon diversity index）、辛普森多樣性指數(shù)（Simpson diversity index）等。 計(jì)算菌群豐度：Chao； 計(jì)算菌群多樣性：Shannon、Simpson。 Chao1：是用chao1 算法估計(jì)群落中含OTU 數(shù)目的指數(shù)，chao1 在生態(tài)學(xué)中常用來(lái)估計(jì)物種總數(shù)，由Chao (1984) 最早提出。Chao1值越大代表物種總數(shù)越多。 計(jì)算公式：Schao1=Sobs+n1(n1-1)/2(n2+1) Shannon-Wiener 曲線，是利用shannon指數(shù)來(lái)進(jìn)行繪制的，反映樣品中微生物多樣性的指數(shù)，利用各樣品的測(cè)序量在不同測(cè)序深度時(shí)的微生物多樣性指數(shù)構(gòu)建曲線，以此反映各樣本在不同測(cè)序數(shù)量時(shí)的微生物多樣性。 當(dāng)曲線趨向平坦時(shí)，說(shuō)明測(cè)序數(shù)據(jù)量足夠大，可以反映樣品中絕大多數(shù)的微生物物種信息。 Simpson：用來(lái)估算樣品中微生物的多樣性指數(shù)之一，在生態(tài)學(xué)中常用來(lái)定量的描述一個(gè)區(qū)域的生物多樣性。Simpson 指數(shù)值越大，說(shuō)明群落多樣性越高。 Rank-Abundance曲線：同時(shí)解釋樣品所含物種的豐富程度和均勻程度。橫坐標(biāo)代表物種排序的數(shù)量；縱坐標(biāo)代表觀測(cè)到的相對(duì)豐度。物種的豐富程度由曲線在橫軸上的長(zhǎng)度來(lái)反映，曲線越寬，表示物種的組成越豐富；物種組成的均勻程度由曲線的形狀來(lái)反映，曲線越平坦，表示物種組成的均勻程度越高。而曲線快速陡然下降表明樣本中的優(yōu)勢(shì)菌群所占比例很高，多樣性較低。 Beta多樣性是生態(tài)系之間的種多樣性，它包含分類單位的比較生態(tài)系內(nèi)的每一獨(dú)特性。即衡量群落之間的差別。Beta多樣性不僅描述生境內(nèi)生物種類的數(shù)量，同時(shí)也考慮到這些種類的相同性及其彼此之間的位置。不同生境間或某一生態(tài)梯度上不同地段間生物種類的相似性越差，Beta生物多樣性越高。 PCoA分析，即主坐標(biāo)分析（PCoA，Principal Coordinate Analysis）可以從復(fù)雜的多維度數(shù)據(jù)中提取和可視化展示主要的高信息量突變,與PCA的主要思想類似。它將距離矩陣的信息轉(zhuǎn)變成在一系列直交軸中顯示，這樣某種變量的最大值體現(xiàn)在第一個(gè)成分坐標(biāo)中，而另一種變量的最大值體現(xiàn)在另一個(gè)成分坐標(biāo)上。距離越近的點(diǎn)代表樣本間相似度越高，理論上生物學(xué)重復(fù)的樣本應(yīng)該聚在一起。 11.系統(tǒng)樹 系統(tǒng)發(fā)生樹（英文：phylogenetic tree或evolutionary tree）是表明被認(rèn)為具有共同祖先的各物種相互間演化關(guān)系的樹，又被稱為系統(tǒng)發(fā)育樹、系統(tǒng)演化樹、系統(tǒng)進(jìn)化樹、種系發(fā)生樹、演化樹、進(jìn)化樹、系統(tǒng)樹。 它用來(lái)表示系統(tǒng)發(fā)生研究的結(jié)果，用它描述物種之間的進(jìn)化關(guān)系。 filter_otus_from_otu_table.py------挑選豐度>1‰的OTU filter_fasta.py------獲得豐度>1‰OTU的fasta文件 make_phylogeny.py-------構(gòu)建系統(tǒng)樹 ggtree.R------系統(tǒng)樹可視化 為了更好的展示系統(tǒng)樹，我們采用了最新的ggtree R包實(shí)現(xiàn)系統(tǒng)樹的可視化。我們主要展示兩個(gè)系統(tǒng)樹，即circos樹和熱圖樹。ggtree.R代碼如下： 12.差異OTU分析 對(duì)找到的OTU，我們會(huì)分析其在各個(gè)樣品中的豐度，分析每個(gè)OTU在對(duì)照組和實(shí)驗(yàn)組中的差異情況。 OTU baseMean log2FoldChange stat pvalue padj taxonomy 1787644 44.801 4.315 4.364 1.28E-05 0.012 k__Bacteria; p__Cyanobacteria; c__Chloroplast; o__Streptophyta; f__; g__; s__ 539866 68.372 4.198 4.058 4.95E-05 0.028 k__Bacteria; p__Proteobacteria; c__Gammaproteobacteria; o__Thiotrichales; f__Piscirickettsiaceae; g__; s__ 特別聲明，此文章是作者原版，如需要轉(zhuǎn)載，請(qǐng)與作者聯(lián)系。大家如果對(duì)該課程該興趣，可以查看作者的用戶名，也可以在掏寶搜索“16s測(cè)序分析視頻”或者“生物信息視頻”。 1.png 2.png 3.png 4.png 5.png 6.png 7.png 8.png 9.png 10.png 11.png 12.png 13.png 14.png 15.png 16.png 17.png 18.png 19.png 20.png 21.png 22.png 23.png 24.png 25.png 26.png 27.png 28.png |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 分子生物 | 統(tǒng)計(jì)軟件資源帖 | 生物信息學(xué) |
» 猜你喜歡
一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿武理材料305分求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料專碩英一數(shù)二306
已經(jīng)有3人回復(fù)
-
085600材料與化工調(diào)劑 324分
已經(jīng)有7人回復(fù)
-
311求調(diào)劑
已經(jīng)有3人回復(fù)
-
26調(diào)劑/材料/英一數(shù)二/總分289/已過(guò)A區(qū)線
已經(jīng)有7人回復(fù)
-
295求調(diào)劑
已經(jīng)有5人回復(fù)
-
354求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料專業(yè)求調(diào)劑
已經(jīng)有3人回復(fù)
-
求材料調(diào)劑
已經(jīng)有7人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
cynlicious
新蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 177.4
- 散金: 5
- 紅花: 1
- 帖子: 236
- 在線: 61.6小時(shí)
- 蟲號(hào): 2428938
- 注冊(cè): 2013-04-21
- 性別: MM
- 專業(yè): 環(huán)境微生物學(xué)
非常感謝 
2樓2015-10-13 10:36:12
whynot2015
木蟲 (小有名氣)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 2862.8
- 紅花: 2
- 帖子: 79
- 在線: 129.6小時(shí)
- 蟲號(hào): 2731604
- 注冊(cè): 2013-10-17
- 性別: GG
- 專業(yè): 微生物學(xué)
3樓2016-03-17 11:17:16
mbin89
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1080.7
- 帖子: 42
- 在線: 14.7小時(shí)
- 蟲號(hào): 4354340
- 注冊(cè): 2016-01-13
- 專業(yè): 作物栽培與耕作學(xué)
4樓2016-09-20 23:07:58
charles08
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 9 (幼兒園)
- 金幣: 5700.3
- 散金: 2103
- 紅花: 27
- 帖子: 1927
- 在線: 603小時(shí)
- 蟲號(hào): 823046
- 注冊(cè): 2009-08-06
- 性別: GG
- 專業(yè): 植物病理學(xué)
5樓2016-12-12 20:53:57
17864859282
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 88.5
- 帖子: 4
- 在線: 7.8小時(shí)
- 蟲號(hào): 5583487
- 注冊(cè): 2017-02-15
- 專業(yè): 土壤生態(tài)學(xué)
★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
|
您好請(qǐng)教您個(gè)問(wèn)題：如果對(duì)細(xì)菌16s測(cè)序，輸出結(jié)果中能不能得到處理與對(duì)照相比哪個(gè)屬或種有差異？謝謝 發(fā)自小木蟲Android客戶端 |
6樓2017-05-13 01:14:42
7樓2018-12-07 16:11:00
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 學(xué)員KgpCIM 的主題更新
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +7 | 手機(jī)用戶 2026-03-17 | 7/350 |
|
|---|---|---|---|---|
|
[考博] 環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
|
[基金申請(qǐng)] 被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 7/350 |
|
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢(mèng)微 2026-03-15 | 3/150 |
|
|
|
[考研]
 機(jī)械專碩調(diào)劑
+3
|
笨笨兔子 2026-03-12 | 3/150 |
|
|
|
[考研] 26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識(shí)可樂(lè) 2026-03-14 | 5/250 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼？~+123456 2026-03-14 | 3/150 |
|
|
|
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8 | dream…… 2026-03-12 | 8/400 |
|
|
|
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機(jī)器學(xué)習(xí) 一志愿哈工大 +4 | 林yaxin 2026-03-12 | 4/200 |
|
|
|
[考研] 求b區(qū)學(xué)校調(diào)劑 +3 | 周56 2026-03-11 | 3/150 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|