−論文−●グリッドコンピューティングによるヒト染色体DNA塩基配列の周期性探索
新製品開発研究所 久米 勝嘉

ヒトのDNAに潜むであろう長大な周期性を解析し、遺伝子病との関連を明らかにするNTTデータ・CellComputing大規模実証実験・遺伝子病治療研究プロジェクト（2002年12月から2003年4月末まで行われた㈱NTTデータによるCellComputing大規模実証実験研究テーマの一つとして東亞合成㈱ 新製品開発研究所

吉田徹彦のプロジェクトテーマが採用1））の成果の一つとして、ヒト染色体DNA塩基配列中に10,000塩基以上の長大な周期性を有するタンデムリピート領域を見出したのでここに報告する。

1 緒

Homeがその応用例として有名である5）。

言

今回、配列の一部分が非常に長い周期で繰り返し出現するタンデムリピートの探索をおこなった。
ヒトゲノム中に数多く存在するこ

2003年4月にヒトゲノム配列の完全解読（全ゲノム配列の99%の2）部分に対して精度99.99％）が終了した 。 今後はこのゲノム配とが知られている繰り返し配列には、一定長の配列が連続して繰列情報を解読し、その機能を明らかにしてゆくことで、最終的にはり返すタンデムリピートと呼ばれるタイプのものと、ゲノム中に数多く医療分野、医薬品開発などに利用することが期待されている。ヒの類似配列が存在している分散型反復配列と呼ばれるタイプのもトゲノム塩基配列に対する研究の多くは遺伝子に着目したものでのが存在する。このうちタンデムリピートについては、先に述べたよあり、ゲノム配列中に存在する遺伝子の発見、遺伝子機能の解

うに周期の短いものについては、各種の遺伝子病と関係している明を行うものである。
しかし、ヒトゲノム配列中の遺伝子領域の占ものがあることが知られているが、周期の長いものについてはまだめる割合は25%程度であり、さらにこの遺伝子領域中でタンパク質検討が行なわれていない。また、塩基配列中に含まれるタンデムリに翻訳されるエクソン領域は約1%と全体に占める比率は小さい。

ピート配列を探索する作業を実施するためのソフトウェアが多数提

また、ヒトゲノム配列中には塩基配列比率で35%以上にもなる様々案されているが、既存の探索ソフトウェアには必要な計算量の問

3）な繰り返し配列が含まれていることが知られている 。この繰り返題から検出可能な繰り返し周期の長さや、取り扱える塩基の全長し配列はハンチントン病の原因となる3塩基配列（CAG）の繰り返等に制約が存在している。たとえばBenson6）らのTandem Repeat

しのように疾病との関係が明らかにされたものも存在するが、そのFinderは、Fig.1に示すようなタンデムリピートを配列と周期につい多くは機能的な意味が判明していないものが大半である。て一定の基準内での変動を許容しながら検出するプログラムであ

このようにヒトゲノム配列には、まだ十分に解析が進んでいないるが検出可能な周期は2,000以下に限られる。

領域や機能の判明していない塩基配列が多数存在しており、こ
こには未知のまだ解読できていない情報やルールが何らかの形でゲノム配列中に記述されていることが予想される。
しかし、ヒトゲ
ノム配列のデータ量は膨大であり、新規なルールの発見のための解析を実用的な期間内で実施することは困難である。
今回、㈱NTTデータによるGridComputing技術のデモンストレーションであるCellComputing大規模実証実験4）に参画することで、Fig.1 Tandem repeat

数万台規模のコンピュータを利用することにより得られる膨大な計算能力を用いて大規模な解析を実施することが可能となった。

長周期のタンデムリピートの探索が困難であるのは、探索対象内のGridComputingとは、ネットワーク上の多数のコンピュータを利用すべての可能性のある繰返し周期について検討する必要があるするための技術の総称であり、今後のコンピュータの利用方法をため計算量が周期の大きさとデータ量に比例して大きくなることと、大きく変える技術として高い関心が寄せられている。その一つのさらに繰り返す塩基配列が同一ではない場合を考慮したり、繰り返応用方法として多数のコンピュータに計算作業を分散させることしの周期が完全に一致しない場合を考慮すると急速に計算量がで超高速のコンピュータとして膨大な計算処理を行うことが期待増加することにある。今回、長周期のタンデムリピートを探索するにされており、地 球 外 知 的 生 命 体 の探 査を目的としたS E T I＠

あたり、1周期の長さが1万塩基以上のものを対象とすることとした。

東亞合成研究年報

2

TREND 2004 第7号

列の組み合わせの位置情報をすべてリストアップする。入力塩基一方、繰り返しの周期が変動することは今回は考慮しないこととした。
また、繰り返し周期が長い場合、配列全体が一様に類似している配列データはfasta形式の塩基配列データを使用する。結果出力フ可能性は非常に低いと考えられるため、Fig.2に示したようなタンデ

ァイルは、比較元配列の開始位置、比較元配列と比較先配列の塩

ムリピート繰り返し単位の配列のうち、その一部分のみが類似して基配列距離、2配列の一致確率が記述される。
また、比較する配列いるタイプの繰り返し配列を探索することとした。これは一定長のの長さ、出力する一致確率の下限値、探索する2配列の間隔の範

塩基配列Seqorgに類似した塩基配列Seqsimが一定の周期間隔D囲のパラメータが制御ファイルで指定可能であり、このパラメータをで出現するタイプの繰り返し領域と表現することもできる。このよう調節することでCellComputing大規模実証実験における各端末な部分配列が一定の周期で繰り返す領域を部分タンデムリピートでの計算分担量の調整をおこなった。

2．
1．
2 使用データと呼ぶこととし、繰り返しの回数nに応じてn回繰り返し部分タンデムリピートと呼ぶことにする。例えばFig.2のケースは4回繰り返し部

解析対象となる塩基配列データは、CellComputing大規模実証分タンデムリピートである。

実験の開始時点で塩基配列の配列解析が終了し、配列データが公開されていた第13,14,20,21及び22番染色体の塩基配列データ
をNCBIのwebサイ
ト内8）のHumanGenomeResourcesより、その時点（2002年10月）での最新の配列データをダウンロードして利用した。
入手した塩基配列データは、一部にまだ塩基配列が決定されていないギャップ部分が含まれているため、ギャップ部分に関してはギャップ長さの分の記号Nとして表現し、配列決定されている部分と結Fig.2 Partly repeated tandem repeat合した。類似配列の探索計算はCellComputing実験参加者のPCで実施されるが、Lateenプログラムは解析対象の塩基配列をすべ

今回の探索では部分タンデムリピートの周期Dが1万以上、部分タてメモリ上に展開するためPCのメモリが少ない場合、動作に不具ンデムリピート中の繰り返す部分配列の長さを100塩基、部分配列合が生じる可能性がある。今回、実験参加者の平均的なメモリ搭の類似度が70%以上のものを探索することとした。

載量を考慮し、いったん結合した染色体配列データを2,500万塩基単位で分割した。このとき分割による探索漏れをなくすため1,2002 実

験

万塩基分をオーバーラップさせて分割をおこなった。このため探索する類似配列の間隔の上限は1,200万塩基とした。

2．
1 GridComputing用プログラムLateen

2．
1．
3 動作パラメータ

2．
1．
1 プログラム

Lateenの動作を制御するパラメータのうち、比較する塩基配列㈱NTTデータのCellComputingは、GridComputingを構成の単位は100塩基、類似配列として採用する配列の一致確率のするための技術としてすでに実績のあるUnited Devices社で下限は0.7とした。並列化のための計算対象の分割は、計算対象開発されたMetaProcessorPlatform技術を採用している。

データ自体の分割、及び分割データ内の計算担当部分の分割によ

CellComputing用のアプリケーションはプログラムの配布先のオり実施した。計算対象データとして、第13,14,20,21及び22番の染TMのみ）

色体配列をそれぞれ2,500万塩基単位で分割した21種の塩基配で動作するプログラムであり、かつMetaProcessorPlatformの列データを利用した。この分割された各塩基配列データについて、ペレーションシステム環境（今回はMicrosoft Windows機能を利用するためのUD関数を追加するためにC言語のソー検討を行う部分タンデムリピートの周期の範囲を109種類に分割した。

スコードが用意できるものである必要がある7）。このため既存のプ

これにより合計1,537種類に計算処理作業を分割した。

ログラムの流用は断念し、探索プログラムの基本部分は新規に作成することとした。この基本部分に㈱NTTデータによりUD関2．
2 CellComputing大規模実証実験

数の付加によるセキュリティ機能、動作状況の把握機能の追加、CellComputing大規模実証実験は、2002年12月20日から2003

動作時に表示されるスクリーンセーバープログラムの追加作業が年 4月3 0日まで実 施された。㈱ N T Tデータ社 内に設 置された

行われ、CellComputing大規模実証実験で利用するプログラム

CellComputingサーバーより、Lateenプログラムと分割された染色Lateen（Large-sequence Analytical Tool Energized by an

体塩基配列データ及び計算範囲を指定したパラメータファイルがExtensive Network）を作成した。

CellComputing実験参加者のPCに配信される。Lateenプログラム

Lateenは、まず探索対象の染色体配列から一定長の塩基配列は配信されたCellComputing参加者の各PCで動作し類似領域のSeqorgを取り出し,Seqorgに一定以上の類似度で一致する同じ長さ探索を実施する。なお、動作中はFig.3のようなスクリーンセーバーの配列Seqsimの位置を探索する。この作業をSeqorgの位置を変化画面が表示される。

させて実施し、染色体内の一定の距離範囲の類似する2つの配東亞合成研究年報

3TREND 2004 第7号

3．
2 長周期リピートの探索結果
Lateenの実行により得られた類似配列情報データを解析し、繰り返し回数が多く、周期の長い部分タンデムリピートの探索を行なった。繰り返す部分配列の一致確率が70％以上で類似配列の繰り返し回数が3回以上、かつ一定間隔で繰り返す周期が1万以上という条件で探索した結果、発見された部分タンデムリピートの総数は11,826個であった。繰り返しの回数については4回繰り返すものが最大であり5回以上繰り返すものは発見できなかった。各染色体毎の3回、4回繰り返しの長周期部分タンデムリピートの数についてはそれぞれTable1,2に示した。条件に適合する長周期のn回繰り返しの部分タンデムリピートの数は、各染色体とも繰り返しの回
Fig.3 CellComputing screen saver

数、配列の一致確率の設定により大きく変化した。

各PCでの計算終了後、結果出力である類似配列の組み合わTable1 Number of 3times repeat region

せデータリストは暗号化され、インターネット経由でCellComputingサーバーに返送される。データの返送を確認した後、再度新しい配列データと計算範囲の指定ファイルが送信され異なるエリアの探索が開始される。
また、今回の実験では各PCで正常に計算が行なわれたかどうかを検証するために同じ設定のパラメータファイルの組み合わせを複数のPCに送付しその結果が一致するかどうかで結果の信頼性を確認した。返送されたデータについては結果の検証後、Table2 Number of 4times repeat region

一括してオフラインで受け取り解析作業の実施をおこなった。

2．
3 データ処理
Lateenにより収集された類似配列の組み合わせデータから、長周期部分タンデムリピートの探索を実施した。探索作業は比較元
配列Seqorgに対する類似配列Seqsimの先頭塩基位置とSeqsimの先頭塩基の位置の差で表される類似配列間の距離についてリス

部分タンデムリピートの周期について最大のものは、3回繰り返しトを作成し、このリストから配列間距離が整数比となるSeq orgとの場合、探索範囲の上限である600万であり、4回繰り返しの場合

Seqsimの組み合わせを探索することで実施した。

は約193万であった。3回繰り返しの部分タンデムリピートの周期について各染色体別にその分布状況をFig.4に示した。部分タンデ

3 結果と考察

ムリピートの数はいずれの染色体においても20万以下の周期のものが多い傾向が見られた。

3．
1 CellComputingでの計算処理結果
第13,14,20,21及び22番染色体の塩基配列データ中の類似領域の探索作業を計算対象を1,537分割して実施した。分割した各
job 1つあたりの実行時間は、参加者の標準的なPC（PentiumIII
1 G H z 相 当 ）で 計 算を 行った 場 合 、約 2 4 時 間 であった 。
CellComputing大規模実証実験に参加したPCの台数は最終的に12,206台であり、これらのPCで分担して冗長度10（ 同じ計算を10台で独立に実行）で計算を実施した。対象範囲の検討に必要な計算は4ヶ月間のCellComputimng大規模実証実験期間中にすべて完了した。Lateenプログラムにより収集された類似配列の組み合わせは、総計293,472,764組であり、出力データサイズはgzip形式で圧縮した状態で15.8GBであった。
Fig.4 Distribution of periodic distance東亞合成研究年報

4TREND 2004 第7号

さらに今回探索された部分タンデムリピートの染色体上での分布をみるために開始位置と周期の分布についてプロットした結果をFig.5-9に示す。染色体塩基配列中の部分タンデムリピートの開始位置の出現場所の分布には偏りがあり特定の領域に集中する傾向が見られた。
しかし、部分タンデムリピートが多く出現する場所は今回検討した染色体毎に異なっており共通する傾向は見られなかった。
また、部分タンデムリピートの周期と開始位置の間の関係についても各染色体間で共通の傾向は見られなかった。

Fig.7 3times repeat distribution in Chromosome20

Fig.5 3times repeat distribution in Chromosome13

Fig.8 3times repeat distribution in Chromosome21

Fig.6 3times repeat distribution in Chromosome14

Fig.9 3times repeat distribution in Chromosome22

東亞合成研究年報

5TREND 2004 第7号

この4回繰り返している100塩基長の配列の由来をTable5に示す。

今回検出された長周期部分タンデムリピートのうち、4回繰返しの部分タンデムリピートは偶然に出現する可能性は非常に低く、なんらかのルールを反映した領域が検出されていると考え、その領域のTable5 Source of repeat region

検討を行なった。今回検出された4回繰返しの部分タンデムリピートの染色体塩基配列分割データ上の位置、周期に関するデータをTable3に示した。染色体塩基配列分割データは1,200万塩基を重複させ先頭から2,500万塩基長単位で分割したものであり、例えば第13番染色体で発見された7箇所の繰り返し部分配列は、いずれChr13̲7は13番染色体の7番めの分割データを示している。

も34塩基を一単位とする繰り返し配列であった。4回繰り返しの部
分タンデムリピートが発見された領域は34塩基単位の類似した配Table3 Start position of each repeated similar sequence

列が少々の変異を伴いながら約50,000塩基ほど繰り返すミニサテ
ライ
ト領域に存在し、発見された4回繰り返し部分タンデムリピートは、この長く続くミニサテライト領域の規則性を反映したものであることが判明した。このミニサテライ
ト領域について、同じくAを赤、Tを黄、Cを緑、Gを青として表現し、幅を170塩基（34塩基×5）とした2D色彩法で表現した結果をFig.13に示す。
また、4回繰り返す部分配列をTable4に、この配列情報を2D色彩法9）で表示したものをFig.10,11,12に示す。
ここで各塩基はAを赤、Tを黄、Cを緑、Gを青として表現した。
また、4回繰り返している部分
配列は図中の白線に挟まれた部分である。2D画像で表現することにより、検出された領域は類似性の高い部分配列が同じ周期で4回繰り返して出現したものであることが確認された。

Table4 Sequences of repeated region

Fig.10 4times repeated region of Chr13̲7Fig.11 4times repeated region of Chr20̲4Fig.13 Minisatellite region of Chrmosome 13

Fig.12 4times repeated region of Chr21̲2

東亞合成研究年報

6TREND 2004 第7号

3．
3 まとめ

また、第20,21番染色体での繰り返し部分配列はいずれもAluファミリーに属する配列であった。Aluはヒトゲノム中に大量に存在する分ヒトDNA塩基配列について、従来検討されたことのない長周期散型反復配列であり、30億塩基のヒトゲノム全体に100万個以上がの繰り返し構造（長周期部分タンデムリピート）の探索をおこなった。

存在するとされている。このため平均では約3,000塩基あたり1回

具体的には、第13,14,20,21及び22番染色体について、100塩基長出現することになる。Aluがランダムに分布すると仮定した場合の3塩基配列が類似度70%以上で10,000塩基以上の一定の周期で3回繰返しの部分タンデムリピート総数の推定値と測定値の比較を回以上繰り返す条件を満たす領域を探索した。5回以上一定の行い、その結果をTable6に示した。

長周期で部分配列が繰り返す領域は存在しなかったが、4回繰り返す領域を7箇所、3回繰り返す領域を11,826箇所で発見した。これら部分タンデムリピートの染色体中での出現位置は一様ではなくTable6 Calculated .vs Observed Repeat region number

偏りが存在し、ゲノム配列中の何らかの特徴、規則性を反映していると考えられる。

謝Calc. =（Number of total possible sequence pair）/（Alu distance）

辞

CellComputing大規模実証実験に参加する機会を与えて頂き、プログラム開発、計算の実施にあたって多大なご支援をいただい実際に発見された3回繰返し、4回繰返しの部分タンデムリピートのた ㈱ N T T データ技 術 開 発 本 部 の 鑓 水リーダー、副 田 様 、

総数はAlu配列がランダムに分布していた場合に予想されるリピ

CellComputingグループメンバーの皆様に深く感謝いたします。
まート領域数に比較して数万分の1と少ない。これはFig.4-9に見らたCellComputing大規模実証実験に参加頂いた参加者の皆様れる長周期部分タンデムリピートの存在位置の分布の大きな偏りに深く感謝いたします。

を反映しているものと考えられる。この偏りの原因はAluの存在密度が染色体中の場所により異なること10）が考えられるが、他にも引用文献等部分タンデムリピートを構成する制約条件が存在している可能性も考えられる。この制約条件の探索する試みの一つとして4回繰

1） CellComputing「遺伝子病治療研究プロジェクト」, http：//

返しの部分タンデムリピート中の繰り返し出現している部分領域をwww.cellcomputing.jp/project/index̲b.html

含む遺伝子を調査した。その結果をTable7に示す。Aluに関係

2） International Consortium Completes Human Genome Project

した部分配列が繰り返している20,21番染色体中の領域の場合、press release, http：//www.ornl.gov/TechResources/

計8箇所の繰り返し出現する部分配列のうち4箇所が何らかの遺Human̲Genome/project/50yr/press4̲2003.htm

伝子中に含まれていた。ヒトゲノム中の遺伝子領域はゲノム中の3） J.C.Venter, et al Science 291,1304（2001）.

25%程度であることから遺伝子が高密度に存在する領域であるこ

4） CellComputing 大規模実証実験,とが長周期の部分タンデムリピートが存在する条件である可能性

http：//www.cellcomputing.jp/test/index2.html

がある。

5） SETI@home, http：//setiathome.ssl.berkeley.edu
6） G.Benson, Nucleic Acids Res., 27, 573（1999）.
7） UNITED DEVICESTM MetaProcessor Platform Version2.2

Table7 Genes related with repeat region

Application Developer s Guide（2002）.
8） NCBI Human Genome Resources,
http：//www.ncbi.nlm.nih.gov
9） T.Yoshida,N.Obata,K.Oosawa, J.Mol.Biol, 298,343（2000）.
10） The Genome International Sequencing Consortium,
Nature 409, 860（ 2001）.

東亞合成研究年報

7

TREND 2004 第7号