如何將GFF轉換為GTF文件

這篇文章主要介紹了如何將GFF轉換為GTF文件，具有一定借鑒價值，感興趣的朋友可以參考下，希望大家閱讀完這篇文章之后大有收獲，下面讓小編帶著大家一起了解一下。

思禮網(wǎng)站制作公司哪家好，找成都創(chuàng)新互聯(lián)！從網(wǎng)頁設計、網(wǎng)站建設、微信開發(fā)、APP開發(fā)、成都響應式網(wǎng)站建設等網(wǎng)站項目制作，到程序開發(fā)，運營維護。成都創(chuàng)新互聯(lián)從2013年創(chuàng)立到現(xiàn)在10年的時間，我們擁有了豐富的建站經(jīng)驗和運維經(jīng)驗，來保證我們的工作的順利進行。專注于網(wǎng)站建設就選成都創(chuàng)新互聯(lián)。

存儲基因和轉錄本的結構信息，gtf和gff3兩種格式都可以。在實際分析時，會需要轉換兩種格式。比如，NCBI 只提供了GFF格式的下載文件，我們需要轉換成GTF文件之后再使用。

完成這一任務，可以自己編寫腳本，也可以借助現(xiàn)成的工具。接下來看下每種工具的使用方法和特點。使用NCBI的GFF文件進行測試，鏈接如下

ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/000/001/405/GCF_000001405.38_GRCh48.p12/GCF_000001405.38_GRCh48.p12_genomic.gff.gz

1. gffread

gffread是由cufflinks的開發(fā)團隊提供的一款讀取gff文件的工具，可以實現(xiàn)gff文件轉換為gtf文件，用法如下

gffread -T GCF_000001405.38_GRCh48.p12_genomic.gff  -o hg38.gtf

生成的gtf示例如下

NC_000001.11    BestRefSeq      exon    11874   12227   .       +       .       transcript_id "rna0"; gene_id "gene0"; gene_name "DDX11L1";

gffread 生成的gtf文件中只提供了exon和CDS這兩種類型的結構信息，第九列的屬性也只有transcript_id, gene_id, gene_name這3種屬性，最關鍵的是， 這些ID沒有任何含義，我們更想要的是基因的Entrez ID和轉錄本的RefSeq 編號。

2. UCSC

UCSC采用GenePred 格式存儲基因和轉錄本的結構信息， 通過UCSC的小工具，我們可以借助GenePred格式來實現(xiàn)GFF轉換成GTF。用法如下

gff3ToGenePred  GCF_000001405.38_GRCh48.p12_genomic.gff  hg38.GenePred
genePredToGtf database hg38.GenePred hg38.gtf

首先用gff3ToGenePred轉換為GenePred格式，然后在利用genePredToGtf轉換為GTF格式。

生成的gtf示例如下

NC_012920.1    hg38.GenePred    transcript    15956    16023    .    -    .    gene_id "gene60958"; transcript_id "rna171196";  gene_name "gene60958";

相比gffreads, UCSC提供的類型更多，包括以下幾種

  1. exon
  2. CDS
  3. start_codon
  4. stop_codon
  5. transcript

雖然區(qū)間的類型變多了，但是屬性的問題和gffread是一樣的，除了gene_name屬性有價值外，其他屬性都沒有價值。

當然，網(wǎng)上也有各種別人編寫的腳本，但是或多或少多都存在問題。最佳的方案還是自己來編寫，首先需要搞清楚GTF文件中，我們需要哪些類型的區(qū)間信息。

對于GTF的實際使用而言，只需要exon區(qū)間信息就可以區(qū)分不同的轉錄本了，而且在定量的過程中，也是只需要參考exon的位置信息。所以，自己編寫轉換腳本，只需要保留eoxn信息。

另外一個問題就是，第九列提供哪些屬性。根據(jù)我的經(jīng)驗，只需要以下6個屬性

1. gene_id

2. gene_name

3. transcript_id

4. transcript_name

5. gene_type

6. transcript_type

gene_id可以用來保存不同數(shù)據(jù)庫中的基因ID,比如NCBI Entrez Id, Ensembl gene Id, 當然也可以和gene_name屬性保持一致；gene_name屬性用來保存gene  symbol,  相比id, symbol在文章中的使用頻率更高。

transcript_id和transcript_name表征轉錄本的id和名稱，可以是RefSeq ID，也可以是Ensembl transcript id, 用于區(qū)分不同的轉錄本。

gene_type和transcript_type表示基因和轉錄本的類型，比如是protein_coding, 還是lncRNA, rRNA等。在分析時，我們通常會根據(jù)類型選擇其中的部分轉錄本來分析，比如只分析蛋白編碼的轉錄本。

以上6種屬性就能夠滿足幾乎100%的場景，對于不同數(shù)據(jù)庫中的文件，只需要自己寫腳本提取這些信息，就可以了。

感謝你能夠認真閱讀完這篇文章，希望小編分享的“如何將GFF轉換為GTF文件”這篇文章對大家有幫助，同時也希望大家多多支持創(chuàng)新互聯(lián)，關注創(chuàng)新互聯(lián)行業(yè)資訊頻道，更多相關知識等著你來學習!

網(wǎng)頁標題：如何將GFF轉換為GTF文件
網(wǎng)頁路徑：http://bm7419.com/article8/igohip.html

成都網(wǎng)站建設公司_創(chuàng)新互聯(lián)，為您提供做網(wǎng)站、手機網(wǎng)站建設、App設計、品牌網(wǎng)站制作、網(wǎng)站設計、網(wǎng)頁設計公司

廣告

聲明：本網(wǎng)站發(fā)布的內容（圖片、視頻和文字）以用戶投稿、用戶轉載內容為主，如果涉及侵權請盡快告知，我們將會在第一時間刪除。文章觀點不代表本網(wǎng)站立場，如需處理請聯(lián)系客服。電話：028-86922220；郵箱：631063699@qq.com。內容未經(jīng)允許不得轉載，或轉載時需注明來源： 創(chuàng)新互聯(lián)