Name: 植物におけるマイクロRNA転写物を正確かつ効率的に分析するバイオインフォマティクスパイプライン (動画) | JoVE Journal
Uploaded: 2020-01-21T12:30:00.000+00:00
Duration: 6 min 34 s
Description: 学ぶ植物マイクロRNAのバイオインフォマティクス解析手法。miRDeep-P2パイプラインを用いて高精度かつ高速に転写物を注釈付けし、大規模ゲノムにおける偽陽性削減と処理時間短縮を実現する計算的手法です。

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

この研究は、北京農林科学アカデミー(KJCX201917、KJCX20180425、KJCX20180204)から中国国立自然科学財団(31621001)からLLに支援されています。

1. インストールとテスト
<ol><li>必要な依存関係をダウンロード: Bowtie222と RNAfold23.コンパイル済みパッケージをお勧めします。<ol>
<li>読み取りマッピングツールである Bowtie2 をホーム サイト (<a href="http://bowtie-bio.sourceforge.net/bowtie2/index.shtml" target="_blank">http://bowtie-bio.sourceforge.net/bowtie2/index.shtml</a>) からダウンロードします。</li>
		<li>RNA二次構造を予測するために使用されるウィーンパッケージのツールであるRNAfold<a href="http://www.tbi.univie.ac.at/~ivo/RNA/" target="_blank">をhttp://www.tbi.univie.ac.at/~ivo/RNA/</a>からダウンロードします。</li>
		<li>miRDP2 をインストールする前に、これら 2 つの依存関係が正しくインストールされていることを確認し、bash 環境ファイル (.bashrc など) をカスタマイズして、これら 2 つの依存関係に対して正しい PATH を設定します。 注: Bowtie24などの他のマッピング ツールも miRDP2 に適しています。Bowtie または Bowtie2 は、バージョン 1.1.3 以降で使用できます。</li>
	</ol></li>
	<li>miRDP2 パッケージをダウンロードするには<a href="https://sourceforge.net/projects/mirdp2/files/latest_version/" target="_blank">、https://sourceforge.net/projects/mirdp2/files/latest_version/</a>に移動し、tarball ファイルをフェッチします。</li>
	<li>miRDP2をインストールする前に、Perlがパス内にあることを確認してください。miRDP2 をインストールするには、ダウンロードした tarball ファイルのすべての内容を 1 つのフォルダ (1.4.2 のコマンド ライン) に抽出し、フォルダ パスを PATH に設定します。 注: miRDP2 を実行するには、8 GB 以上の RAM と 100 GB のストレージを備えたコンピュータまたはコンピューティング ノードをお勧めします。</li>
	<li>MiRDP2 パイプラインをテストします。<ol>
<li>miRDP2 が正しくインストールされているかどうかをテストするには、テスト データと https://sourceforge.net/projects/mirdp2/files/TestData/<a href="https://sourceforge.net/projects/mirdp2/files/TestData/" target="_blank">で見</a>つかった期待される出力を使用します。テストデータは、1つのフォーマットされたGSMシーケンシングファイルと1つのアラビドプシスタリアナゲノムファイルが含まれています。</li>
		<li>ダウンロードしたすべてのファイルを現在の作業ディレクトリに移動します。 mv miRDP2-v*.tar.gz テストデータ.tar.gz ncRNA_rfam.tar.gz <user_selected_folder> cd <user_selected_folder>
</li>
		<li>圧縮された tarball ファイルを抽出します。 tar –xvzf miRDP2-v*.tar.gz tar –xvzf テストデータ.tar.gz tar –xvzf ncRNA_rfam.tar.gz
</li>
		<li>アラビドプシスゲノム参照指数を構築します。 bowtie2-ビルド-f ./テストデータ/TAIR10_genome.fa./テストデータ/TAIR10_genome
</li>
		<li>ncRNA参照インデックスを構築します。 bowtie2-build -f ./ncRNA_rfam.fa ./1.1.3/スクリプト/インデックス/rfam_index
</li>
		<li>miRDP2 パイプラインを実行します。 bash ./1.1.3/miRDP2-v1.1.3_pipeline.bash –g ./テストデータ/TAIR10_genome.fa -i ./ テストデータ/TAIR 10_genome –f ./TestData/GSM2094927.fa – o . 注: 使用される Linux コマンドは太字フォントと斜体フォントで、コマンドラインオプションは斜体です。*は miRDP2 のバージョンを示します (現在のバージョンは 1.1.3 です)。bowtie2-build コマンドには約 10 分かかり、miRDP2 パイプラインは数分以内に終了する必要があります。</li>
	</ol></li>
	<li>テスト出力を確認します。<ol>
<li>'GSM2094927-15-0-10' という名前のフォルダは、すべての中間ファイルと結果を含む<user_selected_folder>に自動的に生成されることに注意してください。</li>
		<li>タブ区切り出力ファイル GSM2094927-15-0-10_filter_P_prediction、予測 miRNA の最終出力には、染色体 ID、ストランド方向、代表読み取り ID、前駆体 ID、成熟 miRNA 位置、前駆体を示す列が含まれていることを確認します。位置、成熟したシーケンス、および前駆体シーケンス。詳細な分析を容易にするために、このファイルから派生した追加のベッド ファイルに注意してください。</li>
		<li>完成した手順に関する情報を提供するファイル "progress_log" と、プログラムの出力と警告を含むファイル "script_log" と "script_err" を確認します。 注: 現在、クラスタ サーバー上の CentOS リリース 6.5、PC Windows システム上の Cygwin 2.6.0 など、2 つの Linux プラットフォームで miRDP2 をテストしており、perl をサポートする同様のシステムでも miRDP2 が動作するはずです。</li>
	</ol></li>
</ol>2. 新しい miRNA の識別
<ol><li>パイプラインを実行する前に、入力読み取りが正しい形式に前処理されていることを確認してください。 注: miRDP2 の新しいバージョン 1.1.3 では、元の FASTQ フォーマット・ファイルを入力として受け入れることができますが、読み取りのフォーマット処理は以前のバージョンと同様に行われます。<ol>
<li>まず、深いシーケンス読み取りの 5' および 3' の端からアダプターを取り外します (存在する場合)。</li>
		<li>次に、ディープ シーケンス読み取りを FASTA 形式に解析します。</li>
		<li>第 3 に、同一のシーケンスで読み取りが単一の一意の FASTA エントリで表されるように冗長性を削除します。</li>
		<li>最後に、すべての FASTA 識別子が一意であることを確認します。各シーケンス識別子は、深いシーケンス データセットで取得された正確なシーケンスのコピー番号を示す '_x' と整数で終わる必要があります。一意の FASTA 識別子を確保する 1 つの方法は、ID に実行中の番号を含めます。参照については、テスト データ<a href="https://sourceforge.net/projects/mirdp2/files/TestData/" target="_blank">(</a>https://sourceforge.net/projects/mirdp2/files/TestData/ のファイル GSM2094927.fa を参照してください。</li>
		<li>正しい形式の読み取りの例については、以下を参照してください。 >read0_x29909 TTGGATTGAAGGGAGCTA >read1_x36974 TTCCACAGCTTCTCTCTCTG >read2_x32635 TTCCACAGCTTCTCTCTCTCTCTTACTT</li>
	</ol>
</li>
	<li>参照インデックスを作成します。<ol>
<li>ゲノム参照のために、時間を節約するために、目的の種のゲノム配列がインデックス化されている場合は、iGenomesウェブサイト<a href="https://support.illumina.com/sequencing/sequencing_software/igenome.html" target="_blank">(https://support.illumina.com/sequencing/sequencing_software/igenome.html)</a>からBowtie2インデックスファイルをダウンロードします。それ以外の場合、ユーザーは参照シーケンスにインデックスを付け、ゲノムシーケンスの再インデックス化が必要になる可能性があるため、プロジェクトが終了するまでしばらくの間インデックスファイルを保持します。ゲノム参照のインデックス作成方法の詳細については、bowtie2 マニュアル (<a href="http://bowtie-bio.sourceforge.net/bowtie2/manual.shtml" target="_blank">http://bowtie-bio.sourceforge.net/bowtie2/manual.shtml</a>) に記載されています。</li>
		<li>別の非miRNA ncRNA指標は、他の非コードRNA断片から騒々しい配列を排除するためにも必要である。このファイルは、rRNA、tRNA、snRNA、およびsnoRNAを含むRfamからの主要なncRNA配列のコレクションである。このインデックスを作成するには、インデックスを正しく配置して名前を付ける必要があります(<miRDP2_version>/スクリプト/インデックス/rfam_index)、パート1.4を参照してください。</li>
	</ol>
</li>
	<li>miRDP2 を実行します。<ol>
<li>miRDP2 を使用してディープ シーケンス データから新しい miRNA を検出するには、パッケージで bash スクリプトを実行して分析パイプラインを開始します (手順 1.4 の例を参照してください)。 <path_to_miRDP2_folder>/miRDP2-v*.*_pipeline.bash –g <genome_file> -i <path_to_index/index_prefix> -f <seq_file > -lt;output_folder> * はパイプライン bash スクリプトのバージョンを示します。変更可能なパラメータには、1) 読み取りがマップできる異なる場所の数、2) bowtie2 を実行するための不一致の数、3) RPM のしきい値 (100 万あたり読み取り) の 3 つのパラメーターがあります。それぞれ-L、-M、および –R オプションを使用してこれらを変更します。詳細な説明はセクション3.1にあります。</li>
	</ol>
</li>
	<li>miRDP2 出力を確認します。<ol>
<li>出力フォルダは <output_folder> の下に自動的に生成され、'<seq_file_name>-15-0-10' という名前になります。最後の 3 つの数値は、それぞれパラメーター 1、2、および 3 の値 (この場合はデフォルト) を示します。ファイル <seq_file_name>_filter_P_prediction には、新しく更新されたプラント miRNA アノテーション基準を満たす最終的な予測 miRNA の情報が含まれています。出力ファイルの形式の詳細については、パート 1.4 で説明します。</li>
	</ol>
</li>
</ol>3. miRDP2 を使用した変更と注意
<ol><li>変更可能なパラメータ<ol>
<li>'-L' オプションを使用して、読み取りをマップできる場所の数の制限を設定します (パラメーター 1)。読み取りマッピングは、繰り返しシーケンスに関連付けられている可能性があり、miRNA に関連付けられている可能性があります。既定の設定は 15 です。特定の種の場合、多数のメンバーを有するmiRNAファミリーがある場合、ゲノムの風景に適応するために最初のパラメータを手動で増加させてもよい。</li>
		<li>'-M' オプションを使用して、ボウタイに許可される不一致を設定します (パラメーター 2)。既定の設定は 0 です。</li>
		<li>'-R' オプションを使用して、成熟 miRNA に対応する可能性のある読み取りのしきい値を設定します (パラメーター 3)。時間の消費と誤検知を減らすには、RPM で読み取りをフィルター処理します。特定の RPM しきい値を超える読み取りのみが、バックグラウンド ノイズではなく miRNA の成熟したシーケンスを表し、さらに分析するために保持されます。既定の設定は 10 RPM です。</li>
		<li>これらのパラメーターを変更すると、パフォーマンスと時間の消費に影響を与える可能性があることに注意してください。一般に、パラメーター 1 と 2 の増加とパラメーター 3 の減少は、より厳しい結果を生成し、実行時間を長くし、その逆も同様です。</li>
	</ol>
</li>
	<li>冗長性と miRNA*<ol>
<li>miRDP2 からの出力 miRNA は、既知の miRNA とは異なる場合があることに注意してください。これは主に、成熟したmiRNAの不均一性またはmiRNAとmiRNA*の相対的な豊富さの2つの理由の1つによるものであることがわかりました。これは、前駆体の最適な長さの選択および既知のmiRNA遺伝子のプロファイリングに影響を与えないことがわかった。</li>
	</ol>
</li>
</ol>

<ol><li>Ghildiyal, M., Zamore, P. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Small+silencing+RNAs%3A+an+expanding+universe%5BTitle%5D)+AND+%22Nature+Reviews+Genetics%22%5BJournal%5D)">Small silencing RNAs: an expanding universe.</a> Nature Reviews Genetics. 10 (2), 94-108 (2009).</li>
<li>Bartel, D. P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((MicroRNAs%3A+target+recognition+and+regulatory+functions%5BTitle%5D)+AND+%22Cell%22%5BJournal%5D)">MicroRNAs: target recognition and regulatory functions.</a> Cell. 136 (2), 215-233 (2009).</li>
<li>Moran, Y., Agron, M., Praher, D., Technau, U. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+evolutionary+origin+of+plant+and+animal+microRNAs%5BTitle%5D)+AND+%22Nature+Ecology+Evolution%22%5BJournal%5D)">The evolutionary origin of plant and animal microRNAs.</a> Nature Ecology Evolution. 1 (3), 27(2017).</li>
<li>Xie, Z., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Expression+of+Arabidopsis+MIRNA+genes%5BTitle%5D)+AND+%22Plant+Physiology%22%5BJournal%5D)">Expression of Arabidopsis MIRNA genes.</a> Plant Physiology. 138 (4), 2145-2154 (2005).</li>
<li>Zhao, X., Zhang, H., Li, L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Identification+and+analysis+of+the+proximal+promoters+of+microRNA+genes+in+Arabidopsis%5BTitle%5D)+AND+%22Genomics%22%5BJournal%5D)">Identification and analysis of the proximal promoters of microRNA genes in Arabidopsis.</a> Genomics. 101 (3), 187-194 (2013).</li>
<li>Bologna, N. G., Mateos, J. L., Bresso, E. G., Palatnik, J. F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+loop-to-base+processing+mechanism+underlies+the+biogenesis+of+plant+microRNAs+miR319+and+miR159%5BTitle%5D)+AND+%22EMBO+JOURNAL%22%5BJournal%5D)">A loop-to-base processing mechanism underlies the biogenesis of plant microRNAs miR319 and miR159.</a> EMBO JOURNAL. 28 (23), 3646-3656 (2009).</li>
<li>Rogers, K., Chen, X. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Biogenesis%2C+turnover%2C+and+mode+of+action+of+plant+microRNAs%5BTitle%5D)+AND+%22Plant+Cell%22%5BJournal%5D)">Biogenesis, turnover, and mode of action of plant microRNAs.</a> Plant Cell. 25 (7), 2383-2399 (2013).</li>
<li>Voinnet, O. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Origin%2C+biogenesis%2C+and+activity+of+plant+microRNAs%5BTitle%5D)+AND+%22Cell%22%5BJournal%5D)">Origin, biogenesis, and activity of plant microRNAs.</a> Cell. 136 (4), 669-687 (2009).</li>
<li>Iwakawa, H. O., Tomari, Y. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+Functions+of+MicroRNAs%3A+mRNA+Decay+and+Translational+Repression%5BTitle%5D)+AND+%22Trends+in+Cell+Biology%22%5BJournal%5D)">The Functions of MicroRNAs: mRNA Decay and Translational Repression.</a> Trends in Cell Biology. 25 (11), 651-665 (2015).</li>
<li>Lee, R. C., Feinbaum, R. L., Ambros, V. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+C.+elegans+heterochronic+gene+lin-4+encodes+small+RNAs+with+antisense+complementarity+to+lin-14%5BTitle%5D)+AND+%22Cell%22%5BJournal%5D)">The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14.</a> Cell. 75 (5), 843-854 (1993).</li>
<li>Wightman, B., Ha, I., Ruvkun, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Posttranscriptional+regulation+of+the+heterochronic+gene+lin-14+by+lin-4+mediates+temporal+pattern+formation+in+C.+elegans%5BTitle%5D)+AND+%22Cell%22%5BJournal%5D)">Posttranscriptional regulation of the heterochronic gene lin-14 by lin-4 mediates temporal pattern formation in C. elegans.</a> Cell. 75 (5), 855-862 (1993).</li>
<li>Friedlander, M. R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Discovering+microRNAs+from+deep+sequencing+data+using+miRDeep%5BTitle%5D)+AND+%22Nature+Biotechnology%22%5BJournal%5D)">Discovering microRNAs from deep sequencing data using miRDeep.</a> Nature Biotechnology. 26 (4), 407-415 (2008).</li>
<li>Yang, X., Li, L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((miRDeep-P%3A+a+computational+tool+for+analyzing+the+microRNA+transcriptome+in+plants%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">miRDeep-P: a computational tool for analyzing the microRNA transcriptome in plants.</a> Bioinformatics. 27 (18), 2614-2615 (2011).</li>
<li>Meyers, B. C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Criteria+for+annotation+of+plant+MicroRNAs%5BTitle%5D)+AND+%22Plant+Cell%22%5BJournal%5D)">Criteria for annotation of plant MicroRNAs.</a> Plant Cell. 20 (12), 3186-3190 (2008).</li>
<li>Yang, X., Zhang, H., Li, L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Global+analysis+of+gene-level+microRNA+expression+in+Arabidopsis+using+deep+sequencing+data%5BTitle%5D)+AND+%22Genomics%22%5BJournal%5D)">Global analysis of gene-level microRNA expression in Arabidopsis using deep sequencing data.</a> Genomics. 98 (1), 40-46 (2011).</li>
<li>Kuang, Z., Wang, Y., Li, L., Yang, X. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((miRDeep-P2%3A+accurate+and+fast+analysis+of+the+microRNA+transcriptome+in+plants%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">miRDeep-P2: accurate and fast analysis of the microRNA transcriptome in plants.</a> Bioinformatics. , (2018).</li>
<li>Kozomara, A., Birgaoanu, M., Griffiths-Jones, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((miRBase%3A+from+microRNA+sequences+to+function%5BTitle%5D)+AND+%22Nucleic+Acids+Research%22%5BJournal%5D)">miRBase: from microRNA sequences to function.</a> Nucleic Acids Research. 47 (1), 155-162 (2019).</li>
<li>Griffiths-Jones, S., Saini, H. K., van Dongen, S., Enright, A. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((miRBase%3A+tools+for+microRNA+genomics%5BTitle%5D)+AND+%22Nucleic+Acids+Research%22%5BJournal%5D)">miRBase: tools for microRNA genomics.</a> Nucleic Acids Research. 36, Database issue 154-158 (2008).</li>
<li>Axtell, M. J., Meyers, B. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Revisiting+Criteria+for+Plant+MicroRNA+Annotation+in+the+Era+of+Big+Data%5BTitle%5D)+AND+%22Plant+Cell%22%5BJournal%5D)">Revisiting Criteria for Plant MicroRNA Annotation in the Era of Big Data.</a> Plant Cell. 30 (2), 272-284 (2018).</li>
<li>Taylor, R. S., Tarver, J. E., Hiscock, S. J., Donoghue, P. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Evolutionary+history+of+plant+microRNAs%5BTitle%5D)+AND+%22Trends+in+Plant+Science%22%5BJournal%5D)">Evolutionary history of plant microRNAs.</a> Trends in Plant Science. 19 (3), 175-182 (2014).</li>
<li>Kozomara, A., Griffiths-Jones, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((miRBase%3A+annotating+high+confidence+microRNAs+using+deep+sequencing+data%5BTitle%5D)+AND+%22Nucleic+Acids+Research%22%5BJournal%5D)">miRBase: annotating high confidence microRNAs using deep sequencing data.</a> Nucleic Acids Research. 42, Database issue 68-73 (2014).</li>
<li>Langmead, B., Salzberg, S. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Fast+gapped-read+alignment+with+Bowtie+2%5BTitle%5D)+AND+%22Nature+Methods%22%5BJournal%5D)">Fast gapped-read alignment with Bowtie 2.</a> Nature Methods. 9 (4), 357-359 (2012).</li>
<li>Lorenz, R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((ViennaRNA+Package+2.0%5BTitle%5D)+AND+%22Algorithms+for+Molecular+Biology%22%5BJournal%5D)">ViennaRNA Package 2.0.</a> Algorithms for Molecular Biology. 6, 26(2011).</li>
<li>Langmead, B., Trapnell, C., Pop, M., Salzberg, S. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Ultrafast+and+memory-efficient+alignment+of+short+DNA+sequences+to+the+human+genome%5BTitle%5D)+AND+%22Genome+Biology%22%5BJournal%5D)">Ultrafast and memory-efficient alignment of short DNA sequences to the human genome.</a> Genome Biology. 10 (3), 25(2009).</li>
<li>An, J., Lai, J., Sajjanhar, A., Lehman, M. L., Nelson, C. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((miRPlant%3A+an+integrated+tool+for+identification+of+plant+miRNA+from+RNA+sequencing+data%5BTitle%5D)+AND+%22BMC+Bioinformatics%22%5BJournal%5D)">miRPlant: an integrated tool for identification of plant miRNA from RNA sequencing data.</a> BMC Bioinformatics. 15, 275(2014).</li>
<li>Lei, J., Sun, Y. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((miR-PREFeR%3A+an+accurate%2C+fast+and+easy-to-use+plant+miRNA+prediction+tool+using+small+RNA-Seq+data%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">miR-PREFeR: an accurate, fast and easy-to-use plant miRNA prediction tool using small RNA-Seq data.</a> Bioinformatics. 30 (19), 2837-2839 (2014).</li>
<li>Evers, M., Huttner, M., Dueck, A., Meister, G., Engelmann, J. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((miRA%3A+adaptable+novel+miRNA+identification+in+plants+using+small+RNA+sequencing+data%5BTitle%5D)+AND+%22BMC+Bioinformatics%22%5BJournal%5D)">miRA: adaptable novel miRNA identification in plants using small RNA sequencing data.</a> BMC Bioinformatics. 16, 370(2015).</li>
<li>Mathelier, A., Carbone, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((MIReNA%3A+finding+microRNAs+with+high+accuracy+and+no+learning+at+genome+scale+and+from+deep+sequencing+data%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">MIReNA: finding microRNAs with high accuracy and no learning at genome scale and from deep sequencing data.</a> Bioinformatics. 26 (18), 2226-2234 (2010).</li>
<li>Zhu, Q. H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+diverse+set+of+microRNAs+and+microRNA-like+small+RNAs+in+developing+rice+grains%5BTitle%5D)+AND+%22Genome+Research%22%5BJournal%5D)">A diverse set of microRNAs and microRNA-like small RNAs in developing rice grains.</a> Genome Research. 18 (9), 1456-1465 (2008).</li>
<li>Fahlgren, N., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((MicroRNA+gene+evolution+in+Arabidopsis+lyrata+and+Arabidopsis+thaliana%5BTitle%5D)+AND+%22Plant+Cell%22%5BJournal%5D)">MicroRNA gene evolution in Arabidopsis lyrata and Arabidopsis thaliana.</a> Plant Cell. 22 (4), 1074-1089 (2010).</li>
<li>Fromm, B., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+Uniform+System+for+the+Annotation+of+Vertebrate+microRNA+Genes+and+the+Evolution+of+the+Human+microRNAome%5BTitle%5D)+AND+%22Annual+Review+of+Genetics%22%5BJournal%5D)">A Uniform System for the Annotation of Vertebrate microRNA Genes and the Evolution of the Human microRNAome.</a> Annual Review of Genetics. 49, 213-242 (2015).</li>
<li>Blevins, T., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Identification+of+Pol+IV+and+RDR2-dependent+precursors+of+24+nt+siRNAs+guiding+de+novo+DNA+methylation+in+Arabidopsis%5BTitle%5D)+AND+%22Elife%22%5BJournal%5D)">Identification of Pol IV and RDR2-dependent precursors of 24 nt siRNAs guiding de novo DNA methylation in Arabidopsis.</a> Elife. 4, 09591(2015).</li>
<li>Zhai, J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+One+Precursor+One+siRNA+Model+for+Pol+IV-Dependent+siRNA+Biogenesis%5BTitle%5D)+AND+%22Cell%22%5BJournal%5D)">A One Precursor One siRNA Model for Pol IV-Dependent siRNA Biogenesis.</a> Cell. 163 (2), 445-455 (2015).</li>
<li>Werner, S., Wollmann, H., Schneeberger, K., Weigel, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Structure+determinants+for+accurate+processing+of+miR172a+in+Arabidopsis+thaliana%5BTitle%5D)+AND+%22Current+Biology%22%5BJournal%5D)">Structure determinants for accurate processing of miR172a in Arabidopsis thaliana.</a> Current Biology. 20 (1), 42-48 (2010).</li>
<li>Mateos, J. L., Bologna, N. G., Chorostecki, U., Palatnik, J. F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Identification+of+microRNA+processing+determinants+by+random+mutagenesis+of+Arabidopsis+MIR172a+precursor%5BTitle%5D)+AND+%22Current+Biology%22%5BJournal%5D)">Identification of microRNA processing determinants by random mutagenesis of Arabidopsis MIR172a precursor.</a> Current Biology. 20 (1), 49-54 (2010).</li>
<li>Vitsios, D. M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Mirnovo%3A+genome-free+prediction+of+microRNAs+from+small+RNA+sequencing+data+and+single-cells+using+decision+forests%5BTitle%5D)+AND+%22Nucleic+Acids+Research%22%5BJournal%5D)">Mirnovo: genome-free prediction of microRNAs from small RNA sequencing data and single-cells using decision forests.</a> Nucleic Acids Research. 45 (21), 177(2017).</li>
</ol>

著者たちは何も開示する必要はない。

NGSの出現により、多様な種29、30におけるsRNAシーケンシングデータの増加量から多数のmiRNA遺伝子座が同定されている。集中型コミュニティデータベースmiRBase21では、寄託されたmiRNAアイテムが過去10年間でほぼ100倍に増加しています。しかし、動物のmiRNAと比較して、植物miRNAは、識別/注釈をより複雑にする多くのユニークな機能を持っています13,14.
まず、植物miRNAの前駆体は、長さと構造がより可変である(補足ファイル2)16。動物miRNA前駆体の比較的均一な長さが70〜90 ntの周りのように、植物前駆体の長さは、いくつかのフォールドによって異なり、数百nts<sup class="x...

生物学における過去20年間で最もエキサイティングな発見の一つは、ゲノム1の多様な機能を調節するsRNA種の増殖的役割である。特に、miRNAは真核生物において20~24-nt sRNAの重要なクラスを構成し、主にライフサイクル開発段階を通じて著名な遺伝子調節因子として転写後のレベルで機能するとともに、刺激およびストレス応答2、3において機能する。植物において、miRNAは、一般に個々の転写単位4、5としてRNAポリメラーゼIIによって転写されるpri-miRNAと呼ばれる一次転写物から生じる。進化的に保存された細胞機械(動物のドロシャ・RNase III、植物におけるDICER様)によって処理され、pri-miRNAは、分子内幹ループ構造を形成する配列を含む即時miRNA前駆体、プレmiRNAに切除される。その後、プレミRNAは、機能鎖、成熟miRNA、およびあまり頻繁に機能しないパートナー、miRNA*2、8からなる二本鎖中間体、すなわちmiRNA二重鎖に処理される。RNA誘導サイレンシング複合体(RISC)にロードされた後、成熟したmiRNAは、配列相補性に基づいてmRNA標的を認識することができ、その結果、負の調節機能2、8をもたらす。miRNAは、ターゲットトランスクリプトを不安定にするか、ターゲット翻訳を防ぐことができますが、前者の方法は植物8、9で支配されています。
線虫カエノルハブディティス・エレガンス10、11における最初のmiRNAの発見以来、特にNGS法の入手後に、miRNA同定とその機能解析に多くの研究が行われてきた。NGS法の広範な応用は、前駆体の幹ループ構造や成熟miRNAおよびmiRNA*における配列読み取りの優先的蓄積など、miRNAのユニークな特徴を捕捉するように設計された計算ツールの利用を大幅に促進しました。その結果、研究者は多様な種のmiRNAを同定する上で顕著な成功を収めています。先に説明した確率モデル12に基づいて、NGSデータから植物miRNAを発見するための最初の計算ツールであるmiRDeep-P13を開発しました。miRDeep-Pは、より可変的な前駆体長と大きなパラロゴスファミリー13、14、15を特徴とするプラントmiRNAのデコードの課題を克服することを特に目的としました。そのリリース後、このプログラムは何千回もダウンロードされ、40以上の植物種16でmiRNA転写物に注釈を付けるために使用されています。miRDeep-PのようなNGSベースのツールによって推進され、パブリックmiRNAリポジトリmiRBase17に登録されたmiRNAの数が劇的に増加しており、2008年18年には約500 miRNA項目(リリース2.0)に比べて38,000以上のmiRNAアイテムがホストされています(リリース22.1)。
しかし、植物miRNAアノテーションから2つの新しい課題が生じている。第一に、偽陽性の高い比率は、次の理由で植物miRNAアノテーション16、19の品質に大きな影響を与えました:1)NGS sRNAライブラリからの内因性短絡干渉RNA(siRNA)の希薄さは、厳格なmiRNAアノテーション基準がないためにmiRNAとして誤って注釈付けされました。2)事前miRNA情報を持たない種については、NGSデータに基づいて予測される偽陽性を排除することは困難である。例としてmiRBaseを使用して、テイラーら20は、パブリックリポジトリ21(リリース21)の植物miRNAエントリの3分の1が説得力のある証拠を欠いて、植物miRNAファミリーの4分の3さえ疑わしいことを発見しました。第二に、大規模で複雑なゲノム16を有する植物miRNAを予測するための非常に時間のかかるプロセスとなる。これらの課題を克服するために、新しいフィルタリング戦略を追加し、スコアリングアルゴリズムを見直し、プラントmiRNAアノテーションの新しい基準を統合することでmiRDeep-Pを更新し、新しいバージョンのmiRDP2をリリースしました。さらに、徐々にゲノムサイズが増加するNGS sRNAデータセット(アラビドプシス、米、トマト、トウモロコシ、小麦)を用いてmiRDP2をテストしました。他の5つの広く使用されているツールとその古いバージョンと比較して、miRDP2はこれらのsRNAデータを解析し、精度と感度を向上させ、miRNA転写物をより速く分析しました。
miRDP2 パッケージの内容 miRDP2 パッケージは、準備された bash スクリプトによって順番に実行する必要がある 6 つの文書化された Perl スクリプトで構成されます。6 つのスクリプトのうち、3 つのスクリプト (convert_bowtie_to_blast.pl、 filter_alignments.pl、およびexcise_candidate.pl) は miRDeep-P から継承されます。その他のスクリプトは、元のバージョンから変更されます。6 つのスクリプトの機能については、以下で説明します。
preprocess_reads.plは、入力読み取りを含め、長すぎる読み取りまたは短すぎる読み取り (&lt;19 nt または &gt;25 nt)、Rfam ncRNA シーケンスと相関する読み取り、および RPM (100 万分の読み取り) による読み取りが 5 未満である。次に、スクリプトは既知の miRNA 成熟配列に関連付けられた読み取りを取得します。入力ファイルは、FASTA/FASTQ 形式のオリジナル読み取りであり、miRNA および ncRNA シーケンスへの読み取りマッピングの bowtie2 出力です。
RPM を計算する式は次のとおりです。
<img alt="Equation 1" src="/files/ftp_upload/59864/59864equ01.jpg">
convert_bowtie_to_blast.plは、ボウタイ形式を BLAST 解析形式に変更します。BLAST 解析形式は、標準の NCBI BLASToutput 形式から派生したカスタム表形式で区切られた形式です。
filter_alignments.plは、ディープ シーケンシング読み取りのアライメントをゲノムにフィルター処理します。部分的な位置合わせとマルチアライメント読み取り (ユーザー指定の周波数カットオフ) をフィルター処理します。基本入力は、BLAST 解析形式のファイルです。
excise_candidate.plは、配置された読み取りをガイドラインとして使用して、参照シーケンスから潜在的な前駆体シーケンスを切り取ります。基本入力は、BLAST 解析形式のファイルと FASTA ファイルです。出力は、FASTA 形式のすべての潜在的な前駆体シーケンスです。
mod-miRDP.plには、植物固有のパラメータを使用してスコアリングシステムを変更することにより、コアmiRDeep-Pアルゴリズムから変更される2つの入力ファイル、シグネチャファイルと構造ファイルが必要です。入力ファイルはドットブラケット前駆体構造ファイルであり、配布署名ファイルを読み取ります。
mod-rm_redundant_meet_plant.plには、chromosome_length、前駆体、mod-miRDP.plによって生成されたoriginal_predictionの 3 つの入力ファイルが必要です。これは、2つの出力ファイル、非冗長予測ファイル、および新たに更新されたプラントmiRNA基準によってフィルタリングされた予測ファイルを生成します。出力ファイルの形式の詳細については、セクション 1.4 で説明します。

<table><tbody><tr><td>コンピュータ/コンピューティング ノード</td><td>N/A</td> N<td>/A</td><td>Perl が必要。少なくとも 8 GB の RAM と 100 GB のストレージを推奨</td></tr></tbody></table>

植物におけるマイクロrna転写物を正確かつ効率的に分析するバイオインフォマティクスパイプライン

A Bioinformatics Pipeline to Accurately and Efficiently Analyze the MicroRNA Transcriptomes in Plants

マイクロRNA(miRNA)は、転写後のレベルで遺伝子発現を調節する上で強力な役割を果たす植物や動物に広く存在する20〜24ヌクレオチド(nt)内因性小さなRNA(sRNA)である。次世代シーケンシング(NGS)法によるsRNAライブラリのシーケンシングは、過去10年間にmiRNA転写物を同定および分析するために広く採用されており、miRNA発見の急速な増加をもたらしました。しかし、配列化されたsRNAライブラリーの深さの増加と植物ゲノムのサイズと複雑さにより、植物miRNAアノテーションには2つの大きな課題が生じます。第一に、他の多くのタイプのsRNA、特に、sRNAライブラリからの短い干渉RNA(siRNA)は、多くの計算ツールによってmiRNAとして誤ってアノトされています。第二に、大規模で複雑なゲノムを持つ植物種におけるmiRNA転写物を解析するための非常に時間のかかるプロセスになります。これらの課題を克服するために、最近、新しいフィルタリング戦略を採用し、スコアリングアルゴリズムを見直し、新たに更新された植物miRNAを組み込むことにより、miRDeep-P(miRNAトランスクリプトーム解析のための一般的なツール)をmiRDeep-P2(略してmiRDP2)にアップグレードしました。注釈の条件。アラビドプシス、米、トマト、トウモロコシ、小麦を含むゲノム複雑度の増加を伴う5つの代表的な植物における配列SRNA集団に対してmiRDP2を試験した。この結果は、miRDP2 がこれらのタスクを非常に高い効率で処理したことを示しています。さらに、miRDP2は感度と精度に関する他の予測ツールを上回っています。まとめると、植物miRNAトランスクリプトムを分析するための迅速かつ正確なツールとしてmiRDP2を実証し、植物のmiRNAに対する注釈を改善するコミュニティを支援する上で有用なツールです。

本明細書に記載されるmiRNAアノテーションパイプラインmiRDP2は、アラビドプシス・タリアナ、オリザ・サティバ(米)、ソラナム・ライコペルシカム(トマト)、ゼア・メイス(トウモロコシ)およびトリティカム・アエスティブム(図1A)を含む、徐々に増加したゲノム長を有する5つの植物種から10の公共sRNA-seqライブラリーに適用される。全体として、各種について、異なる組織からの2つの代表的なsRNAライブラリ(一意の読み取りに崩壊し、プロトコルセクションの詳細)およびそれらは、2つの入力として処理される(表1)。5つのmiRNA計算予測ツール(miRDeep-P13、miRPlant 25、miR-PREFeR26、miRA27、miReNA28...

バイオインフォマティクスパイプライン、すなわちmiRDeep-P2(略してmiRDP2)は、植物miRNA基準とオーバーホールされたアルゴリズムを備えており、特に複雑で大きなゲノムを持つ種に対して、植物のマイクロRNA転写物を正確かつ効率的に分析することができました。

植物におけるマイクロRNA転写物を正確かつ効率的に分析するバイオインフォマティクスパイプライン

MicroRNAs (miRNAs) are 20- to 24-nucleotide (nt) endogenous small RNAs (sRNAs) extensively existing in plants and animals that play potent roles in ...

miRDeep2は、植物開発において重要な転写調節因子であり、環境問題への対応に不可欠な植物マイクロRNAを正確に特定するために使用できます。miRDeep2は著しく短い実行時間を必要とし、特に大きなゲノムを有する植物のマイクロRNAを予測するために感度と精度の両方で優れた性能を発揮します。偽陽性と長い処理時間は、植物マイクロRNAアノテーションにおける課題です。新しいフィールディング戦略を追加し、スコアリングアルゴリズムを見直し、厳しい基準を統合することで、miRDeep2はこれらの問題を克服することができます。高信頼マイクロRNAアノテーションは、ゲノムの多様な機能を調節する上でのマイクロRNAの役割を発見するための基本です。ステップバイステップのガイダンスは、初めてマイクロRNAアオナタを使用するのに役立ちます。経験豊富な研究者にとって、この方法は他のツールよりもmiRDeep2を使用することの利点と利点を理解するのに役立ちます。miRDeep2 パッケージをインストールするには、miRDeep2 Web ページに移動して、tarball ファイルを取得します。次に、ダウンロードしたファイルのすべての内容を 1 つのフォルダーに展開し、フォルダーパスをパスに設定します。miRDeep2パイプラインをテストするには、テストデータと1つのフォーマットされたGSMシーケンシングファイルと1つのシロイヌナナゲノムファイルを含む予想出力をダウンロードし、ダウンロードしたすべてのファイルを現在の作業ディレクトリに移動します。圧縮されたタールボールファイルを抽出した後、シロイヌナズナのゲノム参照インデックスと非コードRNA参照インデックスを構築します。中間ファイルと結果をすべて含むuser_selected_folderに、フォルダが自動的に生成されます。miRDeep2 パイプラインは、テストデータを使用して実行できます。テスト出力を確認するには、タブ区切り出力ファイルを表示します。予測されたマイクロRNAの最終的な出力には、染色体ID、ストランド方向、代表的な読み取りID、前駆体ID、成熟miRNA位置、前駆体位置、成熟配列および前駆体配列を示すカラムが含まれます。次に、完了したステップに関する情報を提供するprogress_logファイルと、プログラムの出力と警告を含むscript_logファイルとscript_errファイルを確認します。パイプラインを実行する前に、入力読み取りが適切な形式に前処理されていることを確認するには、深いシーケンスの読み取りの 5 つの素端と 3 つの素端からアダプターを削除し、すべての FAST A 識別子が一意であることを確認します。各シーケンス識別子は、アンダースコア x と、深いシーケンス データセットで取得された正確なシーケンスのコピー番号を示す整数で終わる必要があります。一意のFAST A識別子を確保するために、参照インデックスを構築ID.Toにランニング番号を含め、対象の種のゲノム配列がインデックス化されている場合は、iGenomesウェブサイトからBowtie 2インデックスファイルをダウンロードしてください。次に、リボソームRNAを含むRNAファムからの主要な非コード配列を含む非マイクロRNA非コードRNAインデックスを構築し、RNA、小核RNA、小核核RNAを伝達し、他の非コードRNAフラグメントから騒々しい配列を除外する。miRDeep2 を使用して、深いシーケンスデータから新しいマイクロRNA を検出するには、パッケージ内で bash スクリプトを実行して分析パイプラインを開始します。読み取り先の異なる場所の数、Bowtie 2 を実行するための不一致数、および 100 万分の 1 読み取りのしきい値を必要に応じて変更できます。miRDeep2 出力を確認するには、自動生成されたoutput_folderのデータを表示します。この代表的な分析では、miRDeep2マイクロRNAアノテーションパイプラインを、示されたとおりに徐々に増加したゲノムサイズを有する5つの植物種から10のパブリックsRNA配列ライブラリに適用した。各種について、異なる組織からの2つの代表的な小RNAライブラリーとそのインデックスゲノム配列が2つの入力として処理される。以前の方法を使用すると、ゲノム処理には100時間以上かかるか、ゲノムの長さのために解析の途中で停止する場合があります。ただし、miRDeep2 は、これらの予測プロセスを、数分から時間まで著しく短い期間で終了します。このテストで使用された2つの配列付きシロイヌナズナシスの小さなRNAについては、miRDeep2は他のツールと比較して感度と精度の両方でより優れたパフォーマンスを発揮しました。プログラムの入力インデックスが正しいことを確認してください。たとえば、Bowtie インデックスでのみ Bowtie を使用し、大きなゲノムには大きなインデックス オプションを使用します。得られたマイクロRNAの標的は、マイクロRNA機能に関する洞察を提供できるシーケンシングデータを使用して予測することができる。miRDeep2は、特定の植物種におけるほとんどのマイクロRNAを正確かつ敏感に同定するために使用することができるので、マイクロRNA機能全体としての役割を検討することができる。

a-bioinformatics-pipeline-to-accurately-efficiently-analyze-microrna

Research

JoVE Journal

Genetics

遺伝学

このビデオでは、植物組織から抽出した全RNAからマイクロRNA(miRNA)を検出し定量するハイブリダイゼーション技術であるノーザンブロットを実演します。...

ノーザンブロットは、植物組織の全RNA抽出物中の特定のマイクロRNAを検出し、定量します。

RNA-Seq解析は、非モデル生物の適応形質の分子基盤を特定するためにますます重要になってきています。ここでは、蚊の飼育からRNAシーケンシング、RNA-Seqデータのバイオインフォマティクス解析まで、休眠時と非休眠型のAedes...

アジアのタイガーモスキートにおける日長休眠のRNA-seqの分析のための実験とバイオインフォマティクス·プロトコル、ヒトスジシマカ

循環マイクロRNAは、最近様々な癌および他の疾患のための有望で新規なバイオマーカーとして浮上している。この記事の目的は、マイクロRNAの循環の豊富さを定量化し、決定するために利用可能である三つの異なるプローブベースのリアルタイムPCRプラットフォームと方法を議論することです。...

プローブベースのリアルタイムPCRは、マイクロRNAの定量的測定のためのアプローチ

このプロトコルは、結合したリボソームの数に基づいて植物組織から転写物を抽出し、分画する簡単な方法を説明しています。これにより、翻訳活性の全体的な推定と、特定のmRNAの翻訳状態の決定が可能になります。...

工場ポリソームプロファイリングのための簡単​​な方法

このレポートでは、高速かつ信頼性の高い携帯電話の miRNAs の mRNA ターゲットの検証方法について説明します。合成ビオチン化核酸 LNA のロックを用いた miRNA を模倣をキャプチャする使用していた方法は、mRNA をターゲットします。プルダウンにストレプトアビジン コートした磁気ビーズを採用するその後、ターゲット mRNA qPCR...

ビオチン ベース プルダウン法検証 mRNA をターゲットの携帯電話の miRNAs

ここで肺式を評価する方法について述べる予測される miRNAs...

肺のマイクロ Rna 発情周期にわたってプロファイリングのオゾン暴露マウスの

ここでは、小型RNAを分離し、マイクロRNAを強化し、ハイスループットシーケンス用のサンプルを準備するためのステップバイステップの戦略について説明します。次に、オープン ソース ツールを使用して、シーケンスの読み取りを処理し、マイクロRNA...

MicroRNAを分離およびシーケンシングし、オープンソースツールを使用してそれらを分析するための完全なパイプライン

我々は、ポテトウイルスX(PVX)ベースのマイクロRNAサイレンシング(VbMS)システムの詳細なプロトコルを提示し、ポテト中の内因性マイクロRNA(miRNA)を機能的に特徴付ける。目的のmiRNAの標的模倣(TM)分子はPVXベクターに統合され、一過性にジャガイモで発現して標的miRNAまたはmiRNAファミリーを沈黙させる。...

ポテトウイルスXベースマイクロRNAサイレンシング(VbMS)ポテト中。

この方法は、全RNA抽出物からmiRNAを検出するためのノーザンハイブリダイゼーション技術の使用を実証する。...

植物マイクロRNAの検出と定量のためのRNAブロット分析

ここでは、1)既知および新規のmiRNAをより正確に同定でき、2)完全に自動化され、自由に利用できる新しい完全に自動化されたmiRNAパイプラインmirMachineを紹介します。ユーザーは、短い送信スクリプトを実行して、完全に自動化された mirMachine...

mirMachine:植物miRNAアノテーションのワンストップショップ

学ぶ植物マイクロRNAのバイオインフォマティクス解析手法。miRDeep-P2パイプラインを用いて高精度かつ高速に転写物を注釈付けし、大規模ゲノムにおける偽陽性削減と処理時間短縮を実現する計算的手法です。

植物におけるマイクロRNA転写物を正確かつ効率的に分析するバイオインフォマティクスパイプライン

さらに動画を探す