マイクロアレイで発現変動を示した遺伝子の中からPubMedから入手したエビデンスリストに含まれる遺伝子を抽出するコードを作成しました。ご紹介します。
解析の前準備
- Rの導入についてはこちらから。
- XML, XLConnectパッケージの導入
＊パッケージの一括導入は下記コードで可能です。
[code language=”R”]
install.packages(c(“XML”, “XLConnect”))
[/code]
- 発現変動を示した遺伝子名(Official Symbol)が入力されたエクセルファイル（試験結果ファイル）
＊遺伝子名(Official Symbol)の他にfold changeやflagの情報が入力されていても問題はありません。必須情報は遺伝子名です。
PubMedのエビデンスリストを入手
遺伝子網羅解析を実施する方はPubMedをご存知かと思います。参考までに、PubMedで遺伝子のエビデンス入手方法を紹介します。
- PubMedにアクセスする。http://www.ncbi.nlm.nih.gov/pubmed
- 検索入力枠左から検索対象のデータベースを「Gene」に変更してエビデンスリストを入手したいキーワードを入力する。検索するとキーワードに関連する遺伝子が検索結果に表示されます。
- 検索結果の右上に表示されるSend toをクリックしFile→FormatをxmlとしCreate Filを押す。そうすると、xmlファイルがダウンロードできます。
＊なお、キーワードによってはxmlファイルサイズが大変大きくなりますのでご注意を。
コードの実行
実行すると、試験結果ファイルからエビデンスリストに含まれる遺伝子を抽出し、最終行にPubMedの遺伝子機能のサマリーが追加された「結果.csv」が出力されます。詳細はコードにコメントを記述していますのでご確認ください。
[code language=”R”]
###ライブラリーの読み込み#####
library(tcltk)
library(XML)
library(XLConnect)
########
###解析遺伝子データで遺伝子名列番号を入力#####
AnaGeneName
少しでも、あなたの実験解析が楽になりますように！！
もし、不明点がありましたらこちらからお問い合わせください。
Amazon audibleの登録の紹介
プライム会員限定で2024年7月22日まで3か月無料体験キャンペーン開催中です。無料体験後は月額1,500円で聞き放題です。なお、聞き放題対象外の本はAudible会員であれば非会員価格の３０％引きで購入することが可能です。
Amazon audibleはプロのナレーターが朗読した本をアプリで聞くことができるサービスで、オフライン再生も可能です。通勤や作業のお供にAmazon audibleのご登録はいかがでしょうか。
・AmazonのAudible
