生物資訊學
生物醫學組學背景知識及測序原理
[編輯]1.生物大分子與中心法則
2.生物資訊學簡介與測序技術的發展歷程、技術原理、測序策略與基本流程:一代、二代、三代測序、單細胞測序
3.測學原理:基因組、轉錄組、表觀組、宏基因組、宏轉錄組、蛋白質組、代謝組
生物資訊資料格式與資料庫
[編輯]1.生物資訊格式:fastq/fasta、SAM/BAM、vcf、gtf/gff、bed、maf、測序資料的格式
2.生物資訊資料庫:序列資料庫(NCBI、UCSC、Ensembl)、綜合資料庫(ENCODE)、癌症資料庫(TCGA)
3.常用專業功能資料庫使用方法:GO、COG、UniProt、KEGG、PubMed、OMIM
生物資訊分析基本技能
[編輯]學習資源:https://zh.wikiversity.org/wiki/School:生物資訊學
生物資訊分析方法與流程及常用軟體
[編輯]- 測序資料的預處理原理、方法及流程
- 生物資訊分析流程的結構與功能概述
高通量測序資料的品質控制
[編輯]高通量測序技術主要指第二代測序技術(Next Generation Sequencing Technology)和第三代測序技術。第二代測序包含了 454測序、Illumina測序和SOLiD測序這3種測序技術;第三代測序技術一般指PacBio單分子即時測序技術。8年以前,主要的測序資料來自 454測序現今;而後,Illumina測序成為主流,到如今不斷發展,成為主要的測序平台;近年來,PacBio測序也不斷發展,成為了主流的第三代測序,特別在微生物基因組測序方面發展迅速。因此,本章重點講述Illumina資料的品質控制。
序列比對
[編輯]高通量測序資料比對,就是將測序得到的reads定位到基因組序列上。由於測序的資料量比較大,因此比對軟體需要能快速將reads比對到參考序列上,並且能並列化執行。對lllumina測序或454測序得到的short reads進行比對的常用軟體有Bowtie、BWA、(基因組、原核轉錄組) HISAT和Tophat(真核轉錄組)。對PacBio測序得到的long reads進行比對的常用軟體是Blasr。
- 使用Bowtie和Bowtie2進行短序列比對
- 使用BWA進行READS比對
- 使用TOPHAT對轉錄組資料進行比對
- 使用HISAT2對轉錄組資料進行比對
- 使用Subread對高通量資料進行比對(未完成)
- 使用STAR對轉錄組資料進行比對:STARsolo用於單細胞轉錄組分析(未完成)
- 使用BLASR對PACBIO資料進行比對(未完成)
有參考基因組的轉錄組分析
[編輯]- 使用featureCounts進行alignment-based的定量
- 使用salmon進行alignment-free的定量:salmon的Alevin用於單細胞轉錄組分析
- 使用Kallisto進行alignment-free的定量:使用bustools的bus命令分析單細胞轉錄組資料
甲基化分析方法與流程
[編輯]全基因組甲基化測序 (WGBS)分析流程
[編輯]WGBS(Whole-genome bisulfite sequencing)被視為甲基化測序的「金標準」,其原理是用 Bisulfite 處理,將基因組中未發生甲基化的 C 鹼基轉換成 U,進行PCR擴增後變成T,與原本具有甲基化修飾的 C 鹼基區分開來,再結合高通量測序技術,與參考序列比對,即可判斷 CpG/CHG/CHH 位點是否發生甲基化,特別適用於繪製單鹼基解析度的全基因組 DNA 甲基化圖譜。
Methyl RAD
[編輯]Methyl RAD技術(Wang et al,Open Biol,2015)使用甲基化修飾依賴性內切酶如FspEI、MspJL、LpnPI、AspBHI等,此類內切聘辨識DNA上發生甲基化的胞曉淀在辨識位置的下游隔一定距離切割雙鏈,若DNA雙鏈具有中心對稱甲基化狀態
則可以切割產生一個固定長度的雙鏈DNA片段,對酶切產生的標籤建庫測序,即可進行甲基化位點的定性和相對定量分析。
m6A 甲基化測序(針對轉錄mRNA水平)分析流程
[編輯]m6A(N6-methyladenosine,6-甲基腺嘌呤)是真核生物mRNA最常見的一種轉錄後修飾占到RNA甲基化修飾的80%。已知絕大部分真核生物中mRNA在5' Cap處存在的甲基化修飾,作用包括維持mRNA穩定性、mRNA前體剪下、多腺苷酸化、mRNA運輸與翻譯起始等,而3' polyA發生的甲基化修飾有助於出核轉運翻譯起始以及與polyA結合蛋白一起維持mRNA的結構穩定。
目前,在全轉錄組範圍檢測m6A修飾的主流技術是MeRIP-seq(m6A-seq),該技術使用N6-甲基腺嘌呤抗體富集高甲基化的RNA片段,然後結合高通量測序,在全轉錄組範圍檢測m6A修飾。
單細胞資料分析方法與流程
[編輯]單細胞轉錄組資料分析
[編輯]
單細胞測序是指在單個細胞水平上進行測序。單細胞轉錄組測序(single cell RNA Seq,scRNA-seq)是指對於單個細胞水平上 將mRNA反轉錄擴增後進行高通量測序的技術。
單細胞測序通過在單個細胞水平上進行測序,解決了用組織樣本無法獲得不同細胞間的異質性資訊或樣本量太少無法進行常規測序的難題,為科學家研究單個細胞的行為、機制等提供了新的方向。 單細胞基因組測序主要包括四個步驟:單細胞分離→擴增→高 通量測序→資料分析。
目前常規的測序主要是數百萬甚至更多細胞的混合DNA樣本。這種方法能夠得到基因表達資訊,但是對其進行研究得到的結果只是一群細胞中訊號的平均值,或者只代表其中占優勢數量的細胞資訊,單個細胞獨有的特性被忽視。 單細胞RNA-seq能夠獨立地提供每個細胞的RNA表達譜,並鑑定異質細胞群中的稀有細胞。儘管腫瘤異質性可歸因於累積突變,但即使是遺傳上相同的細胞在相同環境下也可能表現出基因和蛋白表達水平的差異,單細胞RNA-seq 就能夠發現這些稀有個體。比如在腫瘤組織中,腫塊中心的細胞,腫塊周圍的細胞,淋巴轉移灶的細胞,以及遠端轉移的細胞,其基因組和轉錄組等遺傳資訊,是存在差異的。
- 單細胞轉錄組上游分析(從原始資料到表達矩陣)
- 單細胞轉錄組下游分析(基於表達矩陣)
- 單細胞轉錄組線上分析工具及pipline教程匯總
- 單細胞轉錄組資料庫
單細胞scATACseq分析
[編輯]- Cicero:scATACseq分析(結合Monocle),cis-調控網路構建和但細胞染色體發育軌跡
單細胞甲基化分析
[編輯]- MethBank:單細胞甲基化資料庫 http://bigd.big.ac.cn/methbank
單細胞分析工具匯總
[編輯]- 單細胞分析工具匯總:https://www.scrna-tools.org
- 單細胞資料庫及工具:https://www.singlecell.de/index.php/resources/databases/
生物序列比對及典型軟體的使用方法
[編輯]bwa、muscle、MUMmer,核酸、蛋白資料庫介紹與blast搜尋
4.基因組denovo組裝分析流程:常用軟體使用方法:SOAPdenovo、Canu等
5.基因組組裝及基因組注釋分析流程:重複序列、基因結構與功能、ncRNA、細胞器基因組注釋
6.轉錄組組裝常用軟體使用方法:Trinity、Cufflinks、Tophat
7.轉錄組表現組的差異表達分析流程及差異表達分析軟體使用方法
8.高通量資料視覺化常用軟體使用方法
9.分子進化分析常用軟體使用方法
10.重測序變異檢測與注釋分析常用軟體使用方法