1、Noncoding配列探索比較片山俊明 東大醫(yī)科研解析分野情報(bào)利用2005 10/11-12実戦系割C+O+A+8割大人事情內(nèi)容一部割愛(ài)FANTOM3 The Transcriptional Landscape of the Mammalian Genome Science 309:1559-1563 (2 Sep 2005)：哺乳動(dòng)物総合的解析哺乳動(dòng)物総合的解析RNARNA新大陸発見(jiàn)新大陸発見(jiàn)作出RNA（非RNA）、予想超23,000個(gè)以上発見(jiàn)。二重鎖DNA両方鎖転寫RNA（RNA）考數(shù)31,422個(gè)発見(jiàn)。 割転寫Noncoding重要話変重要次構(gòu)造 Noncoding RNA 次構(gòu)造機(jī)能

2、色中次構(gòu)造取探索方法流用分類Tandemly repeated DNA - 反復(fù)型Simple sequence repetitions - ,Blocks of tandemly repeated segments - 長(zhǎng)Dispersed repetitive DNA - 分散型Segmental duplications - 內(nèi)配列重複Transposable elements RNA intermediate (Class I) - Retrotransposon - LTR持、+env 仲間 Retroposon - LTR持、LINE, SINE DNA intermediate

3、(Class II) - IS (insertion sequence) - transposase (integrase) 持 Transposon - IS + 持 不明 萌要素 MITE (miniature inverted repeat trasposable element)分類不明 = 萌（） 動(dòng)物？植物？ 菌類？動(dòng)物？ 粘菌 (南方熊楠 - ChemRuby) 節(jié)足動(dòng)物？環(huán)形動(dòng)物？ 緩歩動(dòng)物門 ( - BioRuby)大人事情內(nèi)容一部割愛(ài) 苔棲、歩 東大駐車場(chǎng)南極 乾 tun 狀態(tài)変形水戻蘇生変形 耐乾燥：何年間乾眠 耐溫度：272151環(huán)境耐 耐線：57萬(wàn)OK （致死線量50

4、0） 耐圧力：微生物死絶3000気圧倍6000気圧OK 耐真空：真空曬OK詳 http:/ Miniature inverted transposable elementTransposon / ISMITEtransposaseTIRDRTIR DR 配列 書 ATGC 並CCTAGGCGAACCTTTAGCAGTAGCGACAAAAGCTACCTAGGCG 回CTAGGCGA 回TAGGCGAA 回 (例：8bp) 中出現(xiàn)頻度一定、？BioRubyCCTAGGCGAACCTTTAGCAGTAGCGACAAAAGCTA#!/usr/bin/env rubyrequire biocount =

5、 Hash.new(0)genome = Bio:Sequence:NA.new( )genome.window_search(8) do |subseq| countsubseq += 1endsorted = count.sort_by |subseq, num| num sorted.each do |subseq, num| puts #subseqt#numendcaagcccc 1ttccaaac 1 :agcga tcg 19cgatcgct 20ggcgatcg 21cgatcgcc 27gcgatcgc 51実際偏BioRuby Ruby言語(yǔ) 指向 言語(yǔ) 日本発（作）言語(yǔ) R

6、uby on Rails ？世界中中 Mac OS X 最初入 BioRuby、用 配列解析 KEGG 今年 ChemRuby 共IPA未踏採(cǎi)択 、今日開(kāi)発者人來(lái)。藍(lán)藻 Anabaena sp. PCC7120 高頻度配列？BLAST數(shù)十保存配列次構(gòu)造手 Stockholm 次構(gòu)造表現(xiàn) Emacs (ralee-mode.el) 利用 http:/www.sanger.ac.uk/Users/sgj/ralee/; /.emacs 設(shè)定例; Emacs lisp /lib/lisp/ 以下掘置整理(let (default-directory /lib/lisp) (normal-top-le

7、vel-add-subdirs-to-load-path); /lib/lisp/ralee/ 以下 ralee-*.el ; .stk 開(kāi) ralee-mode 起動(dòng)設(shè)定(autoload ralee-mode ralee-mode Yay! RNA things t)(setq auto-mode-alist (cons (.stk$ . ralee-mode) auto-mode-alist)Emacs (ralee-mode.el)（作者）Stockholm # STOCKHOLM 1.0 seq1 catgcgaaacgtcaagctgggcatc seq2 cattcgaaaggt

8、catggtgcgcaat seq3 catgggaaaccacacagtggccatt #=GC SS_cons . /少例？(INFERNAL Userguide.pdf )ralee mode C-f, C-b, C-n, C-p速移動(dòng)(20文字,20行単位）文字移動(dòng)矢印C-c C-p, C-c C-o対応 (C-c C-o 方畫面分割). (), C-d配列挿入、削除C-c C-i, C-c C-d全配列対挿入、削除C-c C-fRNAfold 畫像 rna.ps 出力欠點(diǎn)？：pseudoknot 扱INFERNAL CM 作成 配列対 Rfam 作成検索使 http:/www.gen

9、/eddy/infernal/INFERNAL cmbuild% cmbuild mite.cm mite.stk Stokholm mite.stk CM mite.cm 生成INFERNAL cmbuild% cmsearch mite.cm genome.fa 配列対、次構(gòu)造考慮検索INFERNAL良點(diǎn)簡(jiǎn)単INFERNAL悪點(diǎn)遅make 大規(guī)模解析 () 遅並列計(jì)算 http:/www.sanger.ac.uk/Software/Rfam/help/scripts/search/rfam_scan.pl rfam_scan 配列対 Rfam 効率的検索 Pe

10、rl Rfam BLAST 絞込手順1: shotgun.rb KEGG 全生物種 1000bp 11000bp 毎切#!/usr/bin/env rubyrequire bioBio:FlatFile.auto(ARGF) do |ff| ff.each do |entry| name, seq = entry.entry_id, entry.naseq begin i = 0 seq.window_search(11000, 10000) do |subseq| puts subseq.to_fasta(#name:segment_#i*10000+1 (11000bp), 60) i +

11、= 1 end puts seq.subseq(i*10000+1).to_fasta(#name:segment_#i*10000+1, 60) rescue puts seq.to_fasta(#name:segment_all, 60) end endend手順2: Makefile 適當(dāng)數(shù)毎分割% make -j 200# 検索対象配列QUERIES := $(wildcard query/*/*)# 配列結(jié)果生成result/%: query/%.fa bin/rfam_scan.pl -d ./rfam -f tab $ $ echo $ finished.txt# 作結(jié)果名一覧生

12、成infernal: $QUERIES:query/%.fa=result/%結(jié)果月末始 266 生物種中現(xiàn)在 137 生物種orz.終 KEGG DAS 載http:/das.hgc.jp/明日、止。（計(jì)算中） make 大丈夫。再起動(dòng)後% make -j 200続、話 MITE 戻INFERNAL 釣 MITE BLAST, HMMER 高感度、高精度 取、崩端拾 近縁藍(lán)藻種検索 計(jì)150ana 84ava 66284466領(lǐng)域比較 MITE 挿入部位前後 500bp 種間比較MITE挿入遺伝子遺伝子anaavaEMBOSS EMBOSS ？ needle stretcher (needl

13、e 足場(chǎng)合？） water (blast, ssearch 、）needle 大文字小文字保存197313.19841 451 gtcaagttggtgagtcagcaatctgatcaatttcaccataatgccaatat 500 |.|.|.| 1250971.1251 451 gtcaagatggtgagtcagcaatctgatcaatttcaccacaatgtcaat- 498197313.19841 501 TAGGACTTACGCATCCAGGTTGTCTGTTAAGACTGGGTGTAAGGGGGTAA 550 1250971.1251 499 - 498197313.19

14、841 551 GGGTATAGCCCTGTACCCCTACACCCTTCTCCAAACCCTTGATTTTTCGT 600 1250971.1251 499 - 498197313.19841 601 TTTCATGCGTAAGTCCTAaatataagtaaggatgtcaacaaagatgcgag 650 .|.|1250971.1251 499 -gaaggatgtcaacaaagatggatg 522197313.19841 651 gttcttgggaaggcatgattcaatcttgctcaatctt- 687 |.|.| 1250971.1251 523 gttcttggga

15、aggtatgattcaatgttgctcaatctttgttgtctttgct 572興味領(lǐng)域大文字見(jiàn)両端合結(jié)果 ana-specific結(jié)果 ana-ava-common伝播時(shí)期？ MITE 作成 中崩 MITE 群 ClustalW 、內(nèi)部數(shù)変化進(jìn)化速度違部分削 手？XCED 宮田研（隈、加藤 et al.）統(tǒng)合環(huán)境 http:/www.biophys.kyoto-u.ac.jp/katoh/programs/align/xced/ XCED 內(nèi)蔵高精度（MAFFT）手xced 図（xced ）ClustalW XCED FASTA 吐出 ClustalW 読込作成% clustalw

16、mite.dnd R - ape /統(tǒng)計(jì) ape 系統(tǒng)樹(shù)描/src/contrib/Descriptions/ape.html% R install.packages(ape) example(plot.phylo) q()勝手、。問(wèn)題出場(chǎng)合 root R 起動(dòng) sudo R 。ape 系統(tǒng)樹(shù)描lab4ut unrooted tree 用向 axial horizontal 指定font 04 bold italic （）cex font 大試行錯(cuò)誤設(shè)定 tree = read.tree(tree.dnd)# 無(wú)根系統(tǒng)樹(shù)描 postscript(tree_unrooted.ps, horizontal=F, pointsize=5) plot(tree, u, lab4ut=axial, font = 4, cex = 0.5) dev.off()# 有根系統(tǒng)樹(shù)描 postscript(tree_rooted.p