 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 5502 | 回復(fù): 6 | ||||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||||
QQ894064647新蟲 (初入文壇)
|
[交流]
QIIME--16s測序生物信息培訓(xùn)課程(原創(chuàng)) 已有6人參與
|
|||||
|
QIIME--16s測序生物信息培訓(xùn)課程(原創(chuàng)) 名詞解釋 16S rRNA 16S rRNA 基因是編碼原核生物核糖體小亞基的基因，長度約為1542bp，其分子大小適中，突變率小，是細(xì)菌系統(tǒng)分類學(xué)研究中最常用和最有用的標(biāo)志。其序列包含可變區(qū)（Variable region）和保守區(qū)（constant region）。可變區(qū)序列因細(xì)菌而異，且變異程度與細(xì)菌的系統(tǒng)發(fā)育密切相關(guān)，通過提取環(huán)境樣品的DNA，并擴(kuò)增和檢測其中的一段或幾段區(qū)域，可以了解環(huán)境樣品(土壤、腸道、海洋)中細(xì)菌的種類和豐度。 OTU 如果序列之間，比如不同的 16S rRNA序列的相似性高于97%就可以把它定義為一個(gè)OTU，每個(gè)OTU對應(yīng)于一個(gè)不同的16S rRNA序列，也就是每個(gè)OTU對應(yīng)于一個(gè)不同的細(xì)菌（微生物）種。通過OTU分析，就可以知道樣品中的微生物多樣性和不同微生物的豐度。 測序區(qū)域 由于16s rDNA較長（1.5kb），考慮到測序平臺(tái)的讀長限制，我們只能對可變區(qū)進(jìn)行測序。16s rDNA包含有9個(gè)可變區(qū)，分別是v1-v9。一般我們對v3，v3-v4和v6可變區(qū)域進(jìn)行擴(kuò)增和測序。 主要內(nèi)容 QIIME安裝 測序數(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可能沒有開通VT(virtualization technology)，趕緊進(jìn)入BIOS開通吧。 2.Linux常見命令 大家如果沒有用過linux的話，請學(xué)會(huì)一下幾個(gè)常用的linux命令吧。如果大家想對我們的課程有詳細(xì)的詳解的話，大家可以在掏寶搜索“16s測序分析”就能找到我們的課程了。 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)該是唯一的。 第三列必須為LinkerPrimerSequence，為擴(kuò)增樣品的引物。 后面的列可以根據(jù)樣品的特點(diǎn)加以描述，但是每一列必須包含至少兩個(gè)值，就是指定樣品信息，就可以在這個(gè)位置設(shè)置了。 最后一列必須是Description。每一個(gè)樣品不一樣的地方，必須不一樣額。 5.join_paired_ends.py (fastq-join) 6.質(zhì)量過濾 (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文件對序列分類。 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é)中常用來估計(jì)物種總數(shù)，由Chao (1984) 最早提出。Chao1值越大代表物種總數(shù)越多。 計(jì)算公式：Schao1=Sobs+n1(n1-1)/2(n2+1) Shannon-Wiener 曲線，是利用shannon指數(shù)來進(jìn)行繪制的，反映樣品中微生物多樣性的指數(shù)，利用各樣品的測序量在不同測序深度時(shí)的微生物多樣性指數(shù)構(gòu)建曲線，以此反映各樣本在不同測序數(shù)量時(shí)的微生物多樣性。 當(dāng)曲線趨向平坦時(shí)，說明測序數(shù)據(jù)量足夠大，可以反映樣品中絕大多數(shù)的微生物物種信息。 Simpson：用來估算樣品中微生物的多樣性指數(shù)之一，在生態(tài)學(xué)中常用來定量的描述一個(gè)區(qū)域的生物多樣性。Simpson 指數(shù)值越大，說明群落多樣性越高。 Rank-Abundance曲線：同時(shí)解釋樣品所含物種的豐富程度和均勻程度。橫坐標(biāo)代表物種排序的數(shù)量；縱坐標(biāo)代表觀測到的相對豐度。物種的豐富程度由曲線在橫軸上的長度來反映，曲線越寬，表示物種的組成越豐富；物種組成的均勻程度由曲線的形狀來反映，曲線越平坦，表示物種組成的均勻程度越高。而曲線快速陡然下降表明樣本中的優(yōu)勢菌群所占比例很高，多樣性較低。 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)樹。 它用來表示系統(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分析 對找到的OTU，我們會(huì)分析其在各個(gè)樣品中的豐度，分析每個(gè)OTU在對照組和實(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)載，請與作者聯(lián)系。大家如果對該課程該興趣，可以查看作者的用戶名，也可以在掏寶搜索“16s測序分析視頻”或者“生物信息視頻”。 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é) |
» 猜你喜歡
284求調(diào)劑
已經(jīng)有10人回復(fù)
-
一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑
已經(jīng)有4人回復(fù)
-
07化學(xué)280分求調(diào)劑
已經(jīng)有4人回復(fù)
-
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有12人回復(fù)
-
求材料，環(huán)境專業(yè)調(diào)劑
已經(jīng)有3人回復(fù)
-
335求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿吉大化學(xué)322求調(diào)劑
已經(jīng)有4人回復(fù)
-
環(huán)境學(xué)碩288求調(diào)劑
已經(jīng)有8人回復(fù)
-
341求調(diào)劑(一志愿湖南大學(xué)070300)
已經(jīng)有6人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
whynot2015
木蟲 (小有名氣)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 2862.8
- 紅花: 2
- 帖子: 79
- 在線: 129.6小時(shí)
- 蟲號: 2731604
- 注冊: 2013-10-17
- 性別: GG
- 專業(yè): 微生物學(xué)
3樓2016-03-17 11:17:16
mbin89
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1080.7
- 帖子: 42
- 在線: 14.7小時(shí)
- 蟲號: 4354340
- 注冊: 2016-01-13
- 專業(yè): 作物栽培與耕作學(xué)
4樓2016-09-20 23:07:58
17864859282
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 88.5
- 帖子: 4
- 在線: 7.8小時(shí)
- 蟲號: 5583487
- 注冊: 2017-02-15
- 專業(yè): 土壤生態(tài)學(xué)
★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
|
您好請教您個(gè)問題：如果對細(xì)菌16s測序，輸出結(jié)果中能不能得到處理與對照相比哪個(gè)屬或種有差異？謝謝 發(fā)自小木蟲Android客戶端 |
6樓2017-05-13 01:14:42
7樓2018-12-07 16:11:00
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開開心心沒煩惱 2026-03-23 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 335求調(diào)劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|
|
[考研] 求調(diào)劑 +7 | 十三加油 2026-03-21 | 7/350 |
|
|
|
[考研] 一志愿吉大化學(xué)322求調(diào)劑 +3 | 17501029541 2026-03-23 | 4/200 |
|
|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7 | jhybd 2026-03-23 | 11/550 |
|
|
|
[考研] 一志愿211 初試270分 求調(diào)劑 +4 | 谷雨上岸 2026-03-23 | 5/250 |
|
|
|
[考研] 316求調(diào)劑 +7 | 梁茜雯 2026-03-19 | 7/350 |
|
|
|
[考研] 276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷 +8 | 材料學(xué)257求調(diào)劑 2026-03-23 | 9/450 |
|
|
|
[考研] 317求調(diào)劑 +12 | 申子申申 2026-03-19 | 18/900 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3 | 晨昏線與星海 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿華中科技大學(xué)071000，求調(diào)劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|
|
[考研] 311求調(diào)劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
|
[考研] 22 350 本科985求調(diào)劑，求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-18 | 6/300 |
|
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
|
[考研] 生物學(xué)調(diào)劑招人！�。� +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|