(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2024-09-10
(45)【発行日】2024-09-19
(54)【発明の名称】体細胞バリアント検出のための方法および組成物
(51)【国際特許分類】
G16B 20/00 20190101AFI20240911BHJP
【ＦＩ】
G16B20/00
【請求項の数】
61
(21)【出願番号】P 2020572675
(86)(22)【出願日】2019-10-30
(65)【公表番号】
(43)【公表日】2022-01-31
(86)【国際出願番号】 US2019058895
(87)【国際公開番号】W WO2020092591
(87)【国際公開日】2020-05-07
【審査請求日】2022-10-24
(31)【優先権主張番号】62/754,094
(32)【優先日】2018-11-01
(33)【優先権主張国・地域又は機関】US
(73)【特許権者】
【識別番号】500358711
【氏名又は名称】イルミナ インコーポレイテッド
(74)【代理人】
【識別番号】100107489
【弁理士】
【氏名又は名称】大塩 竹志
(72)【発明者】
【氏名】ジュ， ジン ヒュン
【審査官】山崎 誠也
(56)【参考文献】
【文献】特表２０１９－５３５０８１（ＪＰ，Ａ）
【文献】特表２０１８－５００６２５（ＪＰ，Ａ）
【文献】米国特許出願公開第２０１６／０３０００１４（ＵＳ，Ａ１）
【文献】ComprehensiveEvaluation of Illumina’s TruSight? Tumor 170 Panel to Estimate Tumor MutationalBurden，[online]，2017年，[2023年12月26日検索], インターネット<URL:https://jp.illumina.com/content/dam/illumina-marketing/documents/products/documents/tst-170-panel-to-estimate-tumor-mutational-burden-aacr-2017-poster-5358.pdf>
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ｇ１６Ｂ ５／００－９９／００
(57)【特許請求の範囲】
【請求項1】
複数のバリアントにおける体細胞バリアントを識別するためのコンピューター実装方法であって、（ａ）前記コンピューターによって、体細胞バリアントおよび生殖系列バリアントを含む前記複数のバリアントにデータベースフィルターを適用するステップであって、
前記複数のバリアントにおける第１の生殖系列バリアントを決定する工程を含み、前記第１の生殖系列バリアントが各々、閾値対立遺伝子カウントよりも大きいまたはそれと等しい、バリアントの第１の参照セットにおける対立遺伝子カウントを有する、ステップ；
（ｂ）前記コンピューターによって、前記複数のバリアントに近接フィルターを適用するステップであって、
（ｉ）前記複数のバリアントのバリアントを複数のビン中にビニングする工程であって、ゲノムの同じ領域中に位置するバリアントが、同じビン中にビニングされる、工程、
（ｉｉ）前記複数のバリアントにおけるデータベースバリアントを決定する工程であって、データベースバリアントが、バリアントの第２の参照セット中に存在する、工程、および
（ｉｉｉ）前記複数のバリアントにおける第２の生殖系列バリアントを決定する工程であって、前記第２の生殖系列バリアントが各々、前記第２の生殖系列バリアントと同じビン中の少なくとも１つのデータベースバリアントの対立遺伝子頻度の近似範囲内の対立遺伝子頻度を有し、前記近似範囲が、第２の生殖系列バリアントの対立遺伝子頻度から０．０５の最大値および最小値を有する範囲であるか、または前記近似範囲が、第２の生殖系列バリアントの対立遺伝子頻度を中心として、第２の生殖系列バリアントの対立遺伝子頻度の二項分布から２標準偏差の最大値および最小値を有する範囲である、工程
を含む、ステップ；ならびに
（ｃ）前記コンピューターによって、前記複数のバリアントから識別された第１および第２の生殖系列バリアントを除去することによって、前記複数のバリアントにおける体細胞バリアントを決定するステップ
を含む、方法。
【請求項2】
前記コンピューターによって、（ａ）および（ｂ）が連続的に実施される、請求項１に記載の方法。
【請求項3】
前記コンピューターによって、（ｂ）が（ａ）の前に実施される、請求項１に記載の方法。
【請求項4】
前記閾値対立遺伝子カウントが５である、請求項１から３のいずれか一項に記載の方法。
【請求項5】
前記閾値対立遺伝子カウントが１０である、請求項４に記載の方法。
【請求項6】
バリアントの前記第１および第２の参照セットが、同じ参照セットである、請求項１から５のいずれか一項に記載の方法。
【請求項7】
バリアントの前記第１または第２の参照セットが、複数の個体についてのバリアントのデータベースを含む、請求項１から６のいずれか一項に記載の方法。
【請求項8】
バリアントの前記第１または第２の参照セットが、ゲノム集約データベース（ｇｎｏｍＡＤ）および１０００ゲノムデータベースから選択される少なくとも１つのデータベースを含む、請求項１から７のいずれか一項に記載の方法。
【請求項9】
ゲノムの前記同じ領域が、同じ染色体内にある、請求項１から８のいずれか一項に記載の方法。
【請求項10】
ゲノムの前記同じ領域が、同じ染色体腕内にある、請求項１から９のいずれか一項に記載の方法。
【請求項11】
ゲノムの前記同じ領域が、同じ染色体サイトバンド内にある、請求項１から１０のいずれか一項に記載の方法。
【請求項12】
ゲノムの前記同じ領域が、１０Ｍｂの領域内にある、請求項１から１１のいずれか一項に記載の方法。
【請求項13】
近接フィルターを適用する前記ステップが、前記コンピューターによって、０．９よりも大きいまたはそれと等しい対立遺伝子頻度を有する第２の生殖系列バリアントを識別する工程をさらに含む、請求項１から１２のいずれか一項に記載の方法。
【請求項14】
近接フィルターを適用する前記ステップが、前記コンピューターによって、前記複数のバリアントにおける第２の生殖系列バリアントを識別する工程をさらに含み、前記第２の生殖系列バリアントが、バリアントの前記第２の参照セット中に存在するデータベースバリアントである、請求項１から１３のいずれか一項に記載の方法。
【請求項15】
前記近似範囲が、第２の生殖系列バリアントの前記対立遺伝子頻度から０．０５の最大値および最小値を有する範囲である、請求項１から１４のいずれか一項に記載の方法。
【請求項16】
前記近似範囲が、第２の生殖系列バリアントの対立遺伝子頻度を中心として、第２の生殖系列バリアントの前記対立遺伝子頻度の二項分布から２標準偏差の最大値および最小値を有する範囲である、請求項１から１５のいずれか一項に記載の方法。
【請求項17】
前記第２の生殖系列バリアントが、前記第２の生殖系列バリアントと同じビン中の少なくとも５つのデータベースバリアントの対立遺伝子頻度に近接した閾値内の対立遺伝子頻度を有する、請求項１から１６のいずれか一項に記載の方法。
【請求項18】
前記第２の生殖系列バリアントが、前記第２の生殖系列バリアントと同じビン中の少なくとも１０個のデータベースバリアントの対立遺伝子頻度に近接した閾値内の対立遺伝子頻度を有する、請求項１から１７のいずれか一項に記載の方法。
【請求項19】
前記複数のバリアントが、前記コンピューターによって、腫瘍細胞を含む生体試料から配列データを得ることによって識別される、請求項１から１８のいずれか一項に記載の方法。
【請求項20】
前記複数のバリアントが、前記コンピューターによって、前記配列データを参照配列とアラインさせること、および前記配列データ中のバリアントを識別することによって識別される、請求項１９に記載の方法。
【請求項21】
腫瘍細胞を含む前記生体試料が、血清試料、大便試料、血液試料、腫瘍試料から選択される、請求項１９または２０に記載の方法。
【請求項22】
前記腫瘍試料が固定される、請求項２１に記載の方法。
【請求項23】
腫瘍の腫瘍突然変異量を決定するコンピューター実装方法であって、前記コンピューターによって、腫瘍細胞を含む生体試料から配列データを得るステップ；
前記コンピューターによって、前記配列データから複数のバリアントを決定するステップ；および
請求項１から２２のいずれか一項に記載の方法に従って、前記コンピューターによって、複数のバリアントにおける体細胞バリアントの数を決定するステップであって、体細胞バリアントの前記数が、前記腫瘍の前記腫瘍突然変異量である、ステップ
を含む、方法。
【請求項24】
遺伝的変異データを分析するための電子システムであって、プロセッサー上で実行される、腫瘍細胞を含む生体試料由来の配列データから複数のバリアントを識別するように適合されたインフォマティクスモジュールであって、前記複数のバリアントが、体細胞バリアントおよび生殖系列バリアントを含む、インフォマティクスモジュール；
前記複数のバリアントから第１の生殖系列バリアントを除去するように適合されたデータベースフィルターモジュールであって、前記第１の生殖系列バリアントが各々、閾値対立遺伝子カウントよりも大きいまたはそれと等しい、バリアントの第１の参照セットにおける対立遺伝子カウントを有する、データベースフィルターモジュール；
前記複数のバリアントから第２の生殖系列バリアントを除去するように適合された近接フィルターモジュールであって、前記近接フィルターモジュールが、
複数のビンを返すように適合されたビニングサブモジュールであって、各ビンが、ゲノムの同じ領域中に位置する前記複数のバリアントのバリアントを含む、ビニングサブモジュール、
前記複数のバリアントにおけるデータベースバリアントを返すように適合された識別サブモジュールであって、データベースバリアントが、バリアントの第２の参照セット中に存在する、識別サブモジュール、および
前記複数のバリアントから第２の生殖系列バリアントを除去するように適合された除去サブモジュールであって、前記第２の生殖系列バリアントが各々、前記第２の生殖系列バリアントと同じビン中の少なくとも１つのデータベースバリアントの対立遺伝子頻度の近似範囲内の対立遺伝子頻度を有し、前記近似範囲が、第２の生殖系列バリアントの対立遺伝子頻度から０．０５の最大値および最小値を有する範囲であるか、または前記近似範囲が、第２の生殖系列バリアントの対立遺伝子頻度を中心として、第２の生殖系列バリアントの対立遺伝子頻度の二項分布から２標準偏差の最大値および最小値を有する範囲である、除去サブモジュール
を含む、近接フィルターモジュール；ならびに
前記複数のバリアントから除去されなかったバリアントを返すように適合されたディスプレイモジュール
を含む、電子システム。
【請求項25】
インフォマティクスモジュールが、バリアントアノテーションツールを含む、請求項２４に記載のシステム。
【請求項26】
前記閾値対立遺伝子カウントが５である、請求項２４または２５に記載のシステム。
【請求項27】
前記閾値対立遺伝子カウントが１０である、請求項２６に記載のシステム。
【請求項28】
バリアントの前記第１および第２の参照セットが、同じ参照セットである、請求項２４から２７のいずれか一項に記載のシステム。
【請求項29】
バリアントの前記第１または第２の参照セットが、複数の個体についてのバリアントのデータベースを含む、請求項２４から２８のいずれか一項に記載のシステム。
【請求項30】
バリアントの前記第１または第２の参照セットが、ゲノム集約データベース（ｇｎｏｍＡＤ）および１０００ゲノムデータベースから選択される少なくとも１つのデータベースを含む、請求項２４から２９のいずれか一項に記載のシステム。
【請求項31】
ゲノムの前記同じ領域が、同じ染色体内にある、請求項２４から３０のいずれか一項に記載のシステム。
【請求項32】
ゲノムの前記同じ領域が、同じ染色体腕内にある、請求項２４から３１のいずれか一項に記載のシステム。
【請求項33】
ゲノムの前記同じ領域が、同じ染色体サイトバンド内にある、請求項２４から３２のいずれか一項に記載のシステム。
【請求項34】
ゲノムの前記同じ領域が、１０Ｍｂの領域内にある、請求項２４から３３のいずれか一項に記載のシステム。
【請求項35】
前記除去サブモジュールが、前記複数のバリアントから、０．９よりも大きいまたはそれと等しい対立遺伝子頻度を有するバリアントを除去するように適合される、請求項２４から３４のいずれか一項に記載のシステム。
【請求項36】
前記除去サブモジュールが、前記複数のバリアントから、バリアントの前記第２の参照セット中に存在するデータベースバリアントを除去するように適合される、請求項２４から３５のいずれか一項に記載のシステム。
【請求項37】
前記近似範囲が、第２の生殖系列バリアントの前記対立遺伝子頻度から０．０５の最大値および最小値を有する範囲である、請求項２４から３６のいずれか一項に記載のシステム。
【請求項38】
前記近似範囲が、第２の生殖系列バリアントの対立遺伝子頻度を中心として、第２の生殖系列バリアントの前記対立遺伝子頻度の二項分布から２標準偏差の最大値および最小値を有する範囲である、請求項２４から３７のいずれか一項に記載のシステム。
【請求項39】
前記第２の生殖系列バリアントが、前記第２の生殖系列バリアントと同じビン中の少なくとも５つのデータベースバリアントの対立遺伝子頻度に近接した閾値内の対立遺伝子頻度を有する、請求項２４から３８のいずれか一項に記載のシステム。
【請求項40】
前記第２の生殖系列バリアントが、前記第２の生殖系列バリアントと同じビン中の少なくとも１０個のデータベースバリアントの対立遺伝子頻度に近接した閾値内の対立遺伝子頻度を有する、請求項２４から３９のいずれか一項に記載のシステム。
【請求項41】
腫瘍細胞を含む前記生体試料が、血清試料、大便試料、血液試料、腫瘍試料から選択される、請求項２４から４０のいずれか一項に記載のシステム。
【請求項42】
前記腫瘍試料が固定される、請求項４１に記載のシステム。
【請求項43】
複数のバリアントにおける体細胞バリアントを識別するためのコンピューター実装方法であって、前記コンピューターによって、請求項１から２２のいずれか一項に記載の方法を実施するステップを含む、方法。
【請求項44】
複数のバリアントにおける体細胞バリアントを識別するためのコンピューター実装方法であって、（ａ）前記コンピューターによって、腫瘍細胞を含む生体試料由来の配列データから複数のバリアントを受け取るステップであって、前記複数のバリアントが、体細胞バリアントおよび生殖系列バリアントを含む、ステップ；
（ｂ）前記コンピューターによって、前記複数のバリアントにデータベースフィルターを適用するステップであって、
前記複数のバリアントについてドキュメントのインデックスを作成する工程、
前記インデックスを用いてバリアントの第１の参照セットを検索して、前記インデックスにおける第１の生殖系列バリアントを識別する工程であって、前記第１の生殖系列バリアントが各々、閾値対立遺伝子カウントよりも大きいまたはそれと等しい、バリアントの前記第１の参照セットにおける対立遺伝子カウントを有する、工程、および
前記インデックスから識別された第１の生殖系列バリアントを除去して、第１のフィルタリングされたバリアントのインデックスを作成する工程
を含む、ステップ；
（ｃ）前記コンピューターによって、第１のフィルタリングされたバリアントの前記インデックスに近接フィルターを適用するステップであって、
（ｉ）ゲノムの異なる領域について複数のビンを作成する工程、
（ｉｉ）第１のフィルタリングされたバリアントの前記インデックスのバリアントをビニングする工程であって、ゲノムの同じ領域中に位置するバリアントが、同じビン中にビニングされる、工程、
（ｉｉｉ）第１のフィルタリングされたバリアントの前記インデックスを用いてバリアントの第２の参照セットを検索して、第１のフィルタリングされたバリアントの前記インデックスにおけるデータベースバリアントを識別する工程、
（ｉｖ）第２の生殖系列バリアントを識別することによって、第１のフィルタリングされたバリアントの前記インデックスから第２の生殖系列バリアントのインデックスを生成する工程であって、前記第２の生殖系列バリアントが各々、前記第２の生殖系列バリアントと同じビン中の少なくとも１つのデータベースバリアントの対立遺伝子頻度の近似範囲内の対立遺伝子頻度を有し、前記近似範囲が、第２の生殖系列バリアントの対立遺伝子頻度から０．０５の最大値および最小値を有する範囲であるか、または前記近似範囲が、第２の生殖系列バリアントの対立遺伝子頻度を中心として、第２の生殖系列バリアントの対立遺伝子頻度の二項分布から２標準偏差の最大値および最小値を有する範囲である、工程、および
（ｖ）第１のフィルタリングされたバリアントの前記インデックスから識別された第２の生殖系列バリアントを除去して、体細胞バリアントのインデックスを作成し、それによって、前記複数のバリアントにおける体細胞バリアントを識別する工程
を含む、ステップ
を含む、方法。
【請求項45】
前記閾値対立遺伝子カウントが５である、請求項４４に記載の方法。
【請求項46】
前記閾値対立遺伝子カウントが１０である、請求項４５に記載の方法。
【請求項47】
バリアントの前記第１および第２の参照セットが、同じ参照セットである、請求項４４から４６のいずれか一項に記載の方法。
【請求項48】
バリアントの前記第１または第２の参照セットが、複数の個体についてのバリアントのデータベースを含む、請求項４４から４７のいずれか一項に記載の方法。
【請求項49】
バリアントの前記第１または第２の参照セットが、ゲノム集約データベース（ｇｎｏｍＡＤ）および１０００ゲノムデータベースから選択される少なくとも１つのデータベースを含む、請求項４４から４８のいずれか一項に記載の方法。
【請求項50】
ゲノムの前記同じ領域が、同じ染色体内にある、請求項４４から４９のいずれか一項に記載の方法。
【請求項51】
ゲノムの前記同じ領域が、同じ染色体腕内にある、請求項４４から５０のいずれか一項に記載の方法。
【請求項52】
ゲノムの前記同じ領域が、同じ染色体サイトバンド内にある、請求項４４から５１のいずれか一項に記載の方法。
【請求項53】
ゲノムの前記同じ領域が、１０Ｍｂの領域内にある、請求項４４から５２のいずれか一項に記載の方法。
【請求項54】
第２のフィルタリングされたバリアントのインデックスを生成する前記工程が、前記コンピューターによって、０．９よりも大きいまたはそれと等しい対立遺伝子頻度を有する第２の生殖系列バリアントを識別することをさらに含む、請求項４４から５３のいずれか一項に記載の方法。
【請求項55】
第２のフィルタリングされたバリアントのインデックスを生成する前記工程が、前記コンピューターによって、前記複数のバリアントにおける第２の生殖系列バリアントを識別することをさらに含み、前記第２の生殖系列バリアントが、バリアントの前記第２の参照セット中に存在するデータベースバリアントである、請求項４４から５４のいずれか一項に記載の方法。
【請求項56】
前記近似範囲が、第２の生殖系列バリアントの前記対立遺伝子頻度から０．０５の最大値および最小値を有する範囲である、請求項４４から５５のいずれか一項に記載の方法。
【請求項57】
前記近似範囲が、第２の生殖系列バリアントの対立遺伝子頻度を中心として、第２の生殖系列バリアントの前記対立遺伝子頻度の二項分布から２標準偏差の最大値および最小値を有する範囲である、請求項４４から５６のいずれか一項に記載の方法。
【請求項58】
前記第２の生殖系列バリアントが、前記第２の生殖系列バリアントと同じビン中の少なくとも５つのデータベースバリアントの対立遺伝子頻度に近接した閾値内の対立遺伝子頻度を有する、請求項４４から５７のいずれか一項に記載の方法。
【請求項59】
前記第２の生殖系列バリアントが、前記第２の生殖系列バリアントと同じビン中の少なくとも１０個のデータベースバリアントの対立遺伝子頻度に近接した閾値内の対立遺伝子頻度を有する、請求項４４から５８のいずれか一項に記載の方法。
【請求項60】
腫瘍細胞を含む前記生体試料が、血清試料、大便試料、血液試料、腫瘍試料から選択される、請求項４４から５８のいずれか一項に記載の方法。
【請求項61】
前記腫瘍試料が固定される、請求項６０に記載の方法。
【発明の詳細な説明】
【技術分野】
【0001】
発明の分野
本明細書で提供される方法およびシステムの一部の実施形態は、単一の試料から得られた配列データからのバリアント呼び出しに関する。一部の実施形態では、体細胞バリアントは、試料中のバリアントのバリアント対立遺伝子頻度およびゲノム中のその位置に基づいて、生殖系列バリアントから識別され得る。
本明細書で提供される方法およびシステムの一部の実施形態は、単一の試料から得られた配列データからのバリアント呼び出しに関する。一部の実施形態では、体細胞バリアントは、試料中のバリアントのバリアント対立遺伝子頻度およびゲノム中のその位置に基づいて、生殖系列バリアントから識別され得る。
【背景技術】
【0002】
発明の背景
ＤＮＡ突然変異は、がんの原因であり、がんの研究および処置の焦点である。次世代シーケンシング（ＮＧＳ）は、現代のシーケンサーが生成できる膨大な数の読み取りに起因して、ｄｅ ｎｏｖｏ突然変異検出のための有望なテクノロジーである。理論的には、ゲノム試料中の全ての突然変異またはバリアントを、バリアント対立遺伝子頻度（ＶＡＦ）またはゲノム領域とは関係なく、十分な読み取り深度を前提として観察することができる。しかし、自信を持ってバリアントを呼び出すことは、読み取りにおけるノイズに起因して、些末なことではない。いくつかのバイオインフォマティクスツールが、シーケンシング読み取りからバリアントを見出すために開発されており、かかる手順は、典型的には、以下の３つの構成要素からなる：読み取り処理、マッピングおよびアラインメント、ならびにバリアント呼び出し。
ＤＮＡ突然変異は、がんの原因であり、がんの研究および処置の焦点である。次世代シーケンシング（ＮＧＳ）は、現代のシーケンサーが生成できる膨大な数の読み取りに起因して、ｄｅ ｎｏｖｏ突然変異検出のための有望なテクノロジーである。理論的には、ゲノム試料中の全ての突然変異またはバリアントを、バリアント対立遺伝子頻度（ＶＡＦ）またはゲノム領域とは関係なく、十分な読み取り深度を前提として観察することができる。しかし、自信を持ってバリアントを呼び出すことは、読み取りにおけるノイズに起因して、些末なことではない。いくつかのバイオインフォマティクスツールが、シーケンシング読み取りからバリアントを見出すために開発されており、かかる手順は、典型的には、以下の３つの構成要素からなる：読み取り処理、マッピングおよびアラインメント、ならびにバリアント呼び出し。
【0003】
読み取り処理のために、通常は読み取りの３’末端近傍の低品質の塩基、およびシーケンシングアダプタなどの外因性配列は、ＤＮＡ試料読み取り処理ツールからトリミングされる。第２に、クリーンにされた読み取りは、バリアントが参照ゲノムに由来し得る場所を決定するために、マッピングおよびアラインメントツールを使用してマッピングされ、次いで、塩基ごとにアラインされる。第３のステップである、バリアント呼び出しのプロセスは、ライブラリー調製、試料富化、シーケンシングおよびマッピング／アラインメントから生じるアーチファクトから、実際のバリアントを分離するために使用される。配列データからのバリアント呼び出しの改善された方法が、引き続き必要とされている。
【発明の概要】
【課題を解決するための手段】
【0004】
発明の要旨
一部の実施形態は、複数のバリアントにおける体細胞バリアントを識別するための方法であって、（ａ）体細胞バリアントおよび生殖系列バリアントを含む複数のバリアントを得るステップ；（ｂ）複数のバリアントにデータベースフィルターを適用するステップであって、複数のバリアントにおける第１の生殖系列バリアントを決定する工程を含み、第１の生殖系列バリアントが各々、閾値対立遺伝子カウントよりも大きいまたはそれと等しい、バリアントの第１の参照セットにおける対立遺伝子カウントを有する、ステップ；（ｃ）複数のバリアントに近接フィルターを適用するステップであって、（ｉ）複数のバリアントのバリアントを複数のビン中にビニングする工程であって、ゲノムの同じ領域中に位置するバリアントが、同じビン中にビニングされる、工程、（ｉｉ）複数のバリアントにおけるデータベースバリアントを決定する工程であって、データベースバリアントが、バリアントの第２の参照セット中に存在する、工程、および（ｉｉｉ）複数のバリアントにおける第２の生殖系列バリアントを決定する工程であって、第２の生殖系列バリアントが各々、第２の生殖系列バリアントと同じビン中の少なくとも１つのデータベースバリアントの対立遺伝子頻度の近似範囲内の対立遺伝子頻度を有する、工程を含む、ステップ；ならびに（ｄ）複数のバリアントから識別された第１および第２の生殖系列バリアントを除去することによって、複数のバリアントにおける体細胞バリアントを決定するステップを含む、方法を含む。
一部の実施形態は、複数のバリアントにおける体細胞バリアントを識別するための方法であって、（ａ）体細胞バリアントおよび生殖系列バリアントを含む複数のバリアントを得るステップ；（ｂ）複数のバリアントにデータベースフィルターを適用するステップであって、複数のバリアントにおける第１の生殖系列バリアントを決定する工程を含み、第１の生殖系列バリアントが各々、閾値対立遺伝子カウントよりも大きいまたはそれと等しい、バリアントの第１の参照セットにおける対立遺伝子カウントを有する、ステップ；（ｃ）複数のバリアントに近接フィルターを適用するステップであって、（ｉ）複数のバリアントのバリアントを複数のビン中にビニングする工程であって、ゲノムの同じ領域中に位置するバリアントが、同じビン中にビニングされる、工程、（ｉｉ）複数のバリアントにおけるデータベースバリアントを決定する工程であって、データベースバリアントが、バリアントの第２の参照セット中に存在する、工程、および（ｉｉｉ）複数のバリアントにおける第２の生殖系列バリアントを決定する工程であって、第２の生殖系列バリアントが各々、第２の生殖系列バリアントと同じビン中の少なくとも１つのデータベースバリアントの対立遺伝子頻度の近似範囲内の対立遺伝子頻度を有する、工程を含む、ステップ；ならびに（ｄ）複数のバリアントから識別された第１および第２の生殖系列バリアントを除去することによって、複数のバリアントにおける体細胞バリアントを決定するステップを含む、方法を含む。
【0005】
一部の実施形態では、（ｂ）および（ｃ）は、連続的に実施される。
【0006】
一部の実施形態では、（ｃ）は、（ｂ）の前に実施される。
【0007】
一部の実施形態では、閾値対立遺伝子カウントは、５である。一部の実施形態では、閾値対立遺伝子カウントは、１０である。
【0008】
一部の実施形態では、バリアントの第１および第２の参照セットは、同じ参照セットである。
【0009】
一部の実施形態では、バリアントの第１または第２の参照セットは、複数の個体についてのバリアントのデータベースを含む。一部の実施形態では、バリアントの第１または第２の参照セットは、ゲノム集約データベース（ｇｎｏｍＡＤ）および１０００ゲノムデータベースから選択される少なくとも１つのデータベースを含む。
【0010】
一部の実施形態では、ゲノムの同じ領域は、同じ染色体内にある。一部の実施形態では、ゲノムの同じ領域は、同じ染色体腕内にある。一部の実施形態では、ゲノムの同じ領域は、同じ染色体サイトバンド（cytoband）内にある。一部の実施形態では、ゲノムの同じ領域は、１０Ｍｂの領域内にある。
【0011】
一部の実施形態では、近接フィルターを適用するステップは、０．９よりも大きいまたはそれと等しい対立遺伝子頻度を有する第２の生殖系列バリアントを識別する工程をさらに含む。
【0012】
一部の実施形態では、近接フィルターを適用するステップは、複数のバリアントにおける第２の生殖系列バリアントを識別する工程をさらに含み、第２の生殖系列バリアントは、バリアントの第２の参照セット中に存在するデータベースバリアントである。
【0013】
一部の実施形態では、近似範囲は、第２の生殖系列バリアントの対立遺伝子頻度から０．０５の最大値および最小値を有する範囲である。
【0014】
一部の実施形態では、近似範囲は、第２の生殖系列バリアントの対立遺伝子頻度を中心として、第２の生殖系列バリアントの対立遺伝子頻度の二項分布から２標準偏差の最大値および最小値を有する範囲である。
【0015】
一部の実施形態では、第２の生殖系列バリアントは、第２の生殖系列バリアントと同じビン中の少なくとも５つのデータベースバリアントの対立遺伝子頻度に近接した閾値内の対立遺伝子頻度を有する。一部の実施形態では、第２の生殖系列バリアントは、第２の生殖系列バリアントと同じビン中の少なくとも１０個のデータベースバリアントの対立遺伝子頻度に近接した閾値内の対立遺伝子頻度を有する。
【0016】
一部の実施形態では、（ａ）は、腫瘍細胞を含む生体試料から配列データを得る工程を含む。一部の実施形態は、配列データを参照配列とアラインさせる工程、および配列データ中のバリアントを識別する工程もまた含む。
【0017】
一部の実施形態では、腫瘍細胞を含む生体試料は、血清試料、大便試料、血液試料、腫瘍試料から選択される。一部の実施形態では、腫瘍試料は、固定される。
【0018】
一部の実施形態は、腫瘍の腫瘍突然変異量を決定する方法であって、腫瘍細胞を含む生体試料から配列データを得るステップ；配列データから複数のバリアントを決定するステップ；および上述の実施形態のいずれか１つの方法に従って、複数のバリアントにおける体細胞バリアントの数を決定するステップであって、体細胞バリアントの数が、腫瘍の腫瘍突然変異量である、ステップを含む、方法を含む。
【0019】
一部の実施形態は、腫瘍を処置する方法であって、腫瘍の腫瘍突然変異量を決定する方法に従って、１０体細胞バリアントよりも大きいまたはそれと等しい腫瘍突然変異量を有する腫瘍を決定するステップ；および有効量のチェックポイント阻害剤を投与することによって、腫瘍を処置するステップを含む、方法を含む。
【0020】
一部の実施形態では、腫瘍は、結腸直腸腫瘍、肺腫瘍、子宮内膜腫瘍、子宮腫瘍、胃腫瘍、黒色腫、乳房腫瘍、膵腫瘍、腎臓腫瘍、膀胱腫瘍および脳腫瘍からなる群から選択される。
【0021】
一部の実施形態では、チェックポイント阻害剤は、ＣＴＬＡ－４阻害剤、ＰＤ－１阻害剤およびＰＤ－Ｌ１阻害剤からなる群から選択される。一部の実施形態では、チェックポイント阻害剤は、イピリムマブ、ニボルマブ、ペムブロリズマブ、スパルタリズマブ、アテゾリズマブ、アベルマブおよびデュルバルマブからなる群から選択される。
【0022】
一部の実施形態は、遺伝的変異データを分析するための電子システムであって、プロセッサー上で実行され、腫瘍細胞を含む生体試料由来の配列データから複数のバリアントを識別するように適合されたインフォマティクスモジュールであって、複数のバリアントが、体細胞バリアントおよび生殖系列バリアントを含む、インフォマティクスモジュール；複数のバリアントから第１の生殖系列バリアントを除去するように適合されたデータベースフィルターモジュールであって、第１の生殖系列バリアントが各々、閾値対立遺伝子カウントよりも大きいまたはそれと等しい、バリアントの第１の参照セットにおける対立遺伝子カウントを有する、データベースフィルターモジュール；複数のバリアントから第２の生殖系列バリアントを除去するように適合された近接フィルターモジュールであって、近接フィルターモジュールが、複数のビンを返すように適合されたビニングサブモジュールであって、各ビンが、ゲノムの同じ領域中に位置する複数のバリアントのバリアントを含む、ビニングサブモジュール、複数のバリアントにおけるデータベースバリアントを返すように適合された識別サブモジュールであって、データベースバリアントが、バリアントの第２の参照セット中に存在する、識別サブモジュール、および複数のバリアントから第２の生殖系列バリアントを除去するように適合された除去サブモジュールであって、第２の生殖系列バリアントが各々、第２の生殖系列バリアントと同じビン中の少なくとも１つのデータベースバリアントの対立遺伝子頻度の近似範囲内の対立遺伝子頻度を有する、除去サブモジュールを含む、近接フィルターモジュール；ならびに複数のバリアントから除去されなかったバリアントを返すように適合されたディスプレイモジュールを含む、電子システムを含む。
【0023】
一部の実施形態では、インフォマティクスモジュールは、バリアントアノテーションツールを含む。
【0024】
一部の実施形態では、閾値対立遺伝子カウントは、５である。一部の実施形態では、閾値対立遺伝子カウントは、１０である。
【0025】
一部の実施形態では、バリアントの第１および第２の参照セットは、同じ参照セットである。
【0026】
一部の実施形態では、バリアントの第１または第２の参照セットは、複数の個体についてのバリアントのデータベースを含む。一部の実施形態では、バリアントの第１または第２の参照セットは、ゲノム集約データベース（ｇｎｏｍＡＤ）および１０００ゲノムデータベースから選択される少なくとも１つのデータベースを含む。
【0027】
一部の実施形態では、ゲノムの同じ領域は、同じ染色体内にある。一部の実施形態では、ゲノムの同じ領域は、同じ染色体腕内にある。一部の実施形態では、ゲノムの同じ領域は、同じ染色体サイトバンド内にある。一部の実施形態では、ゲノムの同じ領域は、１０Ｍｂの領域内にある。
【0028】
一部の実施形態では、除去サブモジュールは、複数のバリアントから、０．９よりも大きいまたはそれと等しい対立遺伝子頻度を有するバリアントを除去するように適合される。
【0029】
一部の実施形態では、除去サブモジュールは、複数のバリアントから、バリアントの第２の参照セット中に存在するデータベースバリアントを除去するように適合される。
【0030】
一部の実施形態では、近似範囲は、第２の生殖系列バリアントの対立遺伝子頻度から０．０５の最大値および最小値を有する範囲である。
【0031】
一部の実施形態では、近似範囲は、第２の生殖系列バリアントの対立遺伝子頻度を中心として、第２の生殖系列バリアントの対立遺伝子頻度の二項分布から２標準偏差の最大値および最小値を有する範囲である。
【0032】
一部の実施形態では、第２の生殖系列バリアントは、第２の生殖系列バリアントと同じビン中の少なくとも５つのデータベースバリアントの対立遺伝子頻度に近接した閾値内の対立遺伝子頻度を有する。一部の実施形態では、第２の生殖系列バリアントは、第２の生殖系列バリアントと同じビン中の少なくとも１０個のデータベースバリアントの対立遺伝子頻度に近接した閾値内の対立遺伝子頻度を有する。
【0033】
一部の実施形態では、腫瘍細胞を含む生体試料は、血清試料、大便試料、血液試料、腫瘍試料から選択される。一部の実施形態では、腫瘍試料は、固定される。
【0034】
一部の実施形態は、複数のバリアントにおける体細胞バリアントを識別するためのコンピューター実装方法であって、上述の方法のいずれか１つの方法を実施するステップを含む、方法を含む。
【0035】
一部の実施形態は、複数のバリアントにおける体細胞バリアントを識別するためのコンピューター実装方法であって、（ａ）腫瘍細胞を含む生体試料由来の配列データから複数のバリアントを受け取るステップであって、複数のバリアントが、体細胞バリアントおよび生殖系列バリアントを含む、ステップ；（ｂ）複数のバリアントにデータベースフィルターを適用するステップであって、複数のバリアントについてドキュメントのインデックスを作成する工程、インデックスを用いてバリアントの第１の参照セットを検索して、インデックスにおける第１の生殖系列バリアントを識別する工程であって、第１の生殖系列バリアントが各々、閾値対立遺伝子カウントよりも大きいまたはそれと等しい、バリアントの第１の参照セットにおける対立遺伝子カウントを有する、工程、およびインデックスから識別された第１の生殖系列バリアントを除去して、第１のフィルタリングされたバリアントのインデックスを作成する工程を含む、ステップ；（ｃ）第１のフィルタリングされたバリアントのインデックスに近接フィルターを適用するステップであって、（ｉ）ゲノムの異なる領域について複数のビンを作成する工程、（ｉｉ）第１のフィルタリングされたバリアントのインデックスのバリアントをビニングする工程であって、ゲノムの同じ領域中に位置するバリアントが、同じビン中にビニングされる、工程、（ｉｉｉ）第１のフィルタリングされたバリアントのインデックスを用いてバリアントの第２の参照セットを検索して、第１のフィルタリングされたバリアントのインデックスにおけるデータベースバリアントを識別する工程、（ｉｖ）第２の生殖系列バリアントを識別することによって、第１のフィルタリングされたバリアントのインデックスから第２の生殖系列バリアントのインデックスを生成する工程であって、第２の生殖系列バリアントが各々、第２の生殖系列バリアントと同じビン中の少なくとも１つのデータベースバリアントの対立遺伝子頻度の近似範囲内の対立遺伝子頻度を有する、工程、および（ｖ）第１のフィルタリングされたバリアントのインデックスから識別された第２の生殖系列バリアントを除去して、体細胞バリアントのインデックスを作成し、それによって、複数のバリアントにおける体細胞バリアントを識別する工程を含む、ステップを含む、方法を含む。
【0036】
一部の実施形態では、閾値対立遺伝子カウントは、５である。一部の実施形態では、閾値対立遺伝子カウントは、１０である。
【0037】
一部の実施形態では、バリアントの第１および第２の参照セットは、同じ参照セットである。
【0038】
一部の実施形態では、バリアントの第１または第２の参照セットは、複数の個体についてのバリアントのデータベースを含む。一部の実施形態では、バリアントの第１または第２の参照セットは、ゲノム集約データベース（ｇｎｏｍＡＤ）および１０００ゲノムデータベースから選択される少なくとも１つのデータベースを含む。
【0039】
一部の実施形態では、ゲノムの同じ領域は、同じ染色体内にある。一部の実施形態では、ゲノムの同じ領域は、同じ染色体腕内にある。一部の実施形態では、ゲノムの同じ領域は、同じ染色体サイトバンド内にある。一部の実施形態では、ゲノムの同じ領域は、１０Ｍｂの領域内にある。
【0040】
一部の実施形態では、第２のフィルタリングされたバリアントのインデックスを生成する工程は、０．９よりも大きいまたはそれと等しい対立遺伝子頻度を有する第２の生殖系列バリアントを識別することをさらに含む。
【0041】
一部の実施形態では、第２のフィルタリングされたバリアントのインデックスを生成する工程は、複数のバリアントにおける第２の生殖系列バリアントを識別することをさらに含み、第２の生殖系列バリアントは、バリアントの第２の参照セット中に存在するデータベースバリアントである。
【0042】
一部の実施形態では、近似範囲は、第２の生殖系列バリアントの対立遺伝子頻度から０．０５の最大値および最小値を有する範囲である。
【0043】
一部の実施形態では、近似範囲は、第２の生殖系列バリアントの対立遺伝子頻度を中心として、第２の生殖系列バリアントの対立遺伝子頻度の二項分布から２標準偏差の最大値および最小値を有する範囲である。
【0044】
一部の実施形態では、第２の生殖系列バリアントは、第２の生殖系列バリアントと同じビン中の少なくとも５つのデータベースバリアントの対立遺伝子頻度に近接した閾値内の対立遺伝子頻度を有する。一部の実施形態では、第２の生殖系列バリアントは、第２の生殖系列バリアントと同じビン中の少なくとも１０個のデータベースバリアントの対立遺伝子頻度に近接した閾値内の対立遺伝子頻度を有する。
【0045】
一部の実施形態では、腫瘍細胞を含む生体試料は、血清試料、大便試料、血液試料、腫瘍試料から選択される。一部の実施形態では、腫瘍試料は、固定される。
【図面の簡単な説明】
【0046】
【0047】
【0048】
図２Ｂは、フィルター決定された体細胞バリアント（黒色で塗りつぶされた丸）およびフィルター決定された生殖系列バリアント（灰色で塗りつぶされた丸）による、各バリアントの染色体位置に従う種々のバリアントについてのＶＡＦを示すグラフである。
【0049】
【0050】
【0051】
図４Ｂは、データベースフィルターのみおよび近接フィルターでフィルタリングした、フィルター決定された体細胞バリアント（黒色で塗りつぶされた丸）およびフィルター決定された生殖系列バリアント（灰色で塗りつぶされた丸）による、各バリアントの染色体位置に従う種々のバリアントについてのＶＡＦを示すグラフである。
【0052】
【0053】
【0054】
【0055】
【0056】
図８Ｂは、腫瘍のみアッセイとＷＥＳ腫瘍－正常アッセイとの間の腫瘍突然変異量（ＴＭＢ）の比較を示すグラフである。
【発明を実施するための形態】
【0057】
詳細な説明
本明細書で提供される方法およびシステムの一部の実施形態は、単一の試料から得られた配列データからのバリアント呼び出しに関する。一部の実施形態では、体細胞バリアントは、試料中のバリアントの対立遺伝子頻度およびゲノム中のバリアントの位置に基づいて、生殖系列バリアントから識別され得る。本明細書で使用される場合、「バリアント」には、核酸分子内の多型が含まれ得る。多型には、挿入、欠失、可変長タンデムリピート、単一ヌクレオチド突然変異、および構造的バリアント、例えば、転座、コピー数変異、またはそれらの組合せが含まれ得る。本明細書で使用される場合、「生殖系列バリアント」には、個体の胚細胞および全ての細胞中に存在するバリアントが含まれ得る。本明細書で使用される場合、「体細胞バリアント」には、個体の腫瘍細胞中に存在するが、他の細胞中には存在しないバリアントが含まれ得る。
本明細書で提供される方法およびシステムの一部の実施形態は、単一の試料から得られた配列データからのバリアント呼び出しに関する。一部の実施形態では、体細胞バリアントは、試料中のバリアントの対立遺伝子頻度およびゲノム中のバリアントの位置に基づいて、生殖系列バリアントから識別され得る。本明細書で使用される場合、「バリアント」には、核酸分子内の多型が含まれ得る。多型には、挿入、欠失、可変長タンデムリピート、単一ヌクレオチド突然変異、および構造的バリアント、例えば、転座、コピー数変異、またはそれらの組合せが含まれ得る。本明細書で使用される場合、「生殖系列バリアント」には、個体の胚細胞および全ての細胞中に存在するバリアントが含まれ得る。本明細書で使用される場合、「体細胞バリアント」には、個体の腫瘍細胞中に存在するが、他の細胞中には存在しないバリアントが含まれ得る。
【0058】
伝統的に、体細胞バリアントと生殖系列バリアントとの間のバリアント呼び出しは、腫瘍試料から得られたデータと、対応する正常試料から得られたデータとの間の比較に依存してきた。しかし、伝統的なバリアント呼び出しは、それに見合った試料が入手可能であることおよび２セットのデータが得られることを要求する。本明細書で提供される実施形態は、個体由来の単一の試料から取得した配列データからのバリアント呼び出しに関する。単一の試料を使用することは、それに見合った試料の必要性を低減させ得、腫瘍試料およびそれに見合った正常試料の両方について配列データを得るために必要なコストを低減させ得る。
【0059】
一部の実施形態は、試料、例えば、腫瘍細胞を含む個体由来の試料から配列データを得ること、配列データを参照と比較して、配列データにおける複数のバリアントを識別すること、ならびに１つまたは複数のフィルターをバリアントに適用して、生殖系列バリアントおよび体細胞バリアントを識別することに関する。一部の実施形態では、フィルターには、近接フィルターが含まれ得る。一部の実施形態では、近接フィルターは、ゲノム中のバリアントの位置に従って、複数のバリアントを複数のビン中にビニングすることを含む。ビニングされたバリアントの一部は、バリアントの１つまたは複数の参照セットにおける対応するバリアントの存在によって、生殖系列バリアントとして識別され得る。特徴付けられていないビニングされたバリアントが、特徴付けられていないバリアントと同じビン中の１つまたは複数の識別された生殖系列バリアントの対立遺伝子頻度と類似の対立遺伝子頻度を有する場合、特徴付けられていないビニングされたバリアントは、生殖系列バリアントであると決定され得る。一部の実施形態は、データベースフィルターを適用して、生殖系列バリアントを識別することも含む。データベースフィルターは、バリアントの１つまたは複数の参照セット中の対応するバリアントの対立遺伝子カウントに従って、生殖系列バリアントを識別し得る。一部の実施形態では、データベースフィルターおよび近接フィルターは、生殖系列バリアントを識別するために、複数のバリアントに適用され得る。一部の実施形態では、体細胞バリアントは、生殖系列バリアントとして識別されるバリアントである。体細胞バリアントの数は、腫瘍の腫瘍突然変異量を示し得る。
【0060】
腫瘍突然変異量は、最近の研究が腫瘍突然変異量とチェックポイント阻害剤免疫療法の有効性との間の相関を示した後で、がん治療選択のための重要なバイオマーカーとして出現してきた。腫瘍突然変異量を計算する際には、生殖系列バリアントを識別し、それをフィルタリングして除くことが有用である。生殖系列バリアントには、個体がそれを持って生まれた（または腫瘍と正常細胞との間で共有される）が、参照ゲノムと比較してバリアントとして検出されるバリアントが含まれ得る。これらのバリアントは、腫瘍細胞を正常細胞から識別することに寄与せず、したがって、正確にフィルタリングされて除かれない場合、腫瘍突然変異量の過大評価をもたらし得る。実施形態は、腫瘍について腫瘍突然変異量を決定するステップ、腫瘍突然変異量に従って腫瘍のための処置を選択するステップ、およびそれを必要とする対象に処置を投与するステップを含む。
ある特定の方法
ある特定の方法
【0061】
本明細書で提供される方法およびシステムの一部の実施形態は、体細胞バリアントおよび生殖系列バリアントを含む複数のバリアントにおける体細胞バリアントを識別するための方法に関する。一部の実施形態では、１つまたは複数のフィルターを使用して、複数のバリアントから生殖系列バリアントをフィルタリングすることができる。かかるフィルターの例には、データベースフィルターおよび近接フィルターが含まれる。
【0062】
一部の実施形態では、複数のバリアントにデータベースフィルターを適用することができる。データベースフィルターを使用して、バリアントを生殖系列バリアントとして識別し、複数のバリアントからバリアントを除去することができる。データベースフィルターは、複数のバリアントの特定のバリアントについての、データベースにおける対応するバリアントの対立遺伝子カウントに関連し得る。
【0063】
複数の各バリアントについて、参照データベースは、データベースにおける対応するバリアントについて検索され得る。参照データベースは、複数の個体についてのバリアントのデータベースを含み得る。本明細書で提供される実施形態で有用なデータベースの例には、ｇｎｏｍＡＤエクソームおよびｇｎｏｍＡＤゲノムデータベースを含むゲノム集約データベース（ｇｎｏｍＡＤ）、ならびに１０００ゲノムデータベース（Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｇｅｎｏｍｅ Ｓａｍｐｌｅ Ｒｅｓｏｕｒｃｅ）が含まれる。例えば、その全体が参照によって組み込まれるLek, M., et al., (2016) Nature 536:285-292を参照のこと。総対立遺伝子カウントは、１つまたは複数の参照データベースにおける対応するバリアントについて決定され得る。対立遺伝子カウントは、バリアントが観察されるデータベース内の観察の総数を示し得る。例えば、対応するバリアントについてのデータベースにおける１０の対立遺伝子カウントは、ホモ接合性バリアントについて少なくとも５つの試料、またはヘテロ接合性バリアントについて最大１０個の試料において、対応するバリアントが観察されたことを示している。一部の実施形態では、対立遺伝子カウントは、１つよりも多いデータベースにおいて観察された最も高い対立遺伝子カウントであり得る。ある特定の閾値対立遺伝子カウントよりも大きいまたはそれと等しい対立遺伝子カウントを有する対応するバリアントを有するバリアントは、生殖系列バリアントとして識別され得る。一部の実施形態では、閾値対立遺伝子カウントは、２、３、４、５、６、７、８、９、１０、１１、１２、１３、１４、１５、１６、１７、１８、１９および２０よりも大きいまたはそれと等しい場合がある。
【0064】
一部の実施形態では、複数のバリアントに近接フィルターを適用することができる。データベースフィルターを使用して、バリアントを生殖系列バリアントとして識別し、複数のバリアントからバリアントを除去することができる。近接フィルターは、複数のバリアントのある特定のバリアントの対立遺伝子頻度、ゲノムの領域中のバリアントの位置、およびバリアントの対立遺伝子頻度の、ゲノムの同じ領域における識別された生殖系列バリアントの対立遺伝子頻度との近接に関連し得る。一部の実施形態では、複数のバリアントのバリアントは、ゲノムの同じ領域中に位置するバリアントが、同じビン中に選別またはビニングされるように、複数のビン中に選別またはビニングされ得る。一部の実施形態では、ゲノムの同じ領域は、同じ染色体内、染色体の同じ腕内、同じ染色体サイトバンド内にあり得る。一部の実施形態では、ゲノムの同じ領域は、同じ連続する１００Ｍｂ、５０Ｍｂ、４０Ｍｂ、３０Ｍｂ、２０Ｍｂ、１０Ｍｂ、５Ｍｂ、１Ｍｂ内、または上述の数のうち任意の２つの間の任意の範囲内であり得る。
【0065】
一部の実施形態では、近接フィルターは、どのビニングされたバリアントが生殖系列バリアントとして容易に識別可能かを決定することも含む。例えば、ビニングされたバリアントは、１つまたは複数の参照データベース中に存在する対応するバリアントを有し得、生殖系列バリアントとして識別され得る。
【0066】
一部の実施形態では、近接フィルターは、試料中の閾値頻度よりも大きいまたはそれと等しい対立遺伝子頻度を有するバリアントが生殖系列バリアントであることを決定することを含む。一部のかかる実施形態では、０．７、０．８、０．９または１．０よりも大きいまたはそれと等しい対立遺伝子頻度を有するバリアントは、生殖系列バリアントとして識別され得る。
【0067】
一部の実施形態では、近接フィルターは、生殖系列バリアントとして識別されていないバリアントについての対立遺伝子頻度の近似範囲を決定することを含む。バリアントについての対立遺伝子頻度の近似範囲は、バリアントの対立遺伝子頻度の上および下の対立遺伝子頻度の範囲を含み得る。一部の実施形態では、近似範囲は、バリアントの対立遺伝子頻度から０．０１、０．０２、０．０３、０．０４、０．０５、０．０６、０．０７、０．０８、０．０９、または上述の数のうち任意の２つの間の範囲内の任意の数の最大値および最小値を有する範囲である。例えば、０．２の対立遺伝子頻度および０．０５の近似範囲を有するバリアントについて、近似範囲の最小値および最大値は、それぞれ、０．１５および０．２５の対立遺伝子頻度である。
【0068】
一部の実施形態では、近似範囲は、所与のバリアントについての支持的証拠が二項プロセスによって生成されると仮定して、二項分布の２（ｎ）標準偏差の値によって決定される。例えば、対立遺伝子頻度（ｘ）、カバレッジ（ｙ）を有するバリアントについて、近似範囲（ｚ）は、
ｚ＝ｎ＊ｓｑｒｔ（ｙ＊ｘ＊（１－ｘ））／ｙ
であり得る。
ｚ＝ｎ＊ｓｑｒｔ（ｙ＊ｘ＊（１－ｘ））／ｙ
であり得る。
【0069】
例えば、０．２の対立遺伝子頻度、１００のシーケンシングのカバレッジ／深度を有するバリアントについて、近似範囲は０．０８であり、近似範囲の最小値および最大値は、それぞれ、０．１２および０．２８の対立遺伝子頻度である。一部の実施形態では、近似範囲は、バリアントの対立遺伝子頻度の上および下、０．０５、またはバリアントの対立遺伝子頻度の二項分布から２（ｎ）標準偏差のいずれか高い方である。
【0070】
一部の実施形態では、バリアントが、バリアントと同じビン中の１つまたは複数の識別された生殖系列バリアントの近似範囲内の対立遺伝子頻度を有する場合、バリアントは、生殖系列バリアントとして識別され得る。一部の実施形態では、バリアントが、バリアントと同じビン中の１、２、３、４、５、６、７、８、９または１０個よりも多い識別された生殖系列バリアントの近似範囲内の対立遺伝子頻度を有する場合、バリアントは、生殖系列バリアントとして識別され得る。一部の実施形態では、バリアントが、バリアントと同じビン中の５つよりも多い識別された生殖系列バリアントの近似範囲内の対立遺伝子頻度を有する場合、バリアントは、生殖系列バリアントとして識別され得る。例えば、バリアントが、バリアントと同じビン中の５つよりも多い識別された生殖系列バリアントの近似範囲内の対立遺伝子頻度を有する場合に、バリアントが生殖系列バリアントとして識別される実施形態では、０．２の対立遺伝子頻度、０．０５の近似範囲を有し、したがって、０．１５の範囲最小値および０．２５の範囲最大値を有し、第７染色体を示すビン中にビニングされたバリアントは、５つよりも多い識別された生殖系列バリアントがバリアントの近似範囲中の対立遺伝子頻度を有し、第７染色体を示すビン中にビニングされた場合、生殖系列バリアントとして識別される。
【0071】
一部の実施形態では、近接フィルターは、生殖系列バリアントとして識別されていないバリアントである体細胞バリアントを識別する。一部の実施形態では、腫瘍由来のシーケンシングデータから得られた体細胞バリアントの数は、腫瘍の腫瘍突然変異量である。
【0072】
一部の実施形態では、複数のバリアントにデータベースフィルターまたは近接フィルターを適用して、生殖系列バリアントを識別し、複数のバリアントからそれを除去することができる。一部の実施形態では、データベースフィルターおよび近接フィルターは、連続的に適用され得る。例えば、かかるデータベースフィルターの出力物は、近接フィルターの入力物として使用され得る。逆に、近接フィルターの出力物は、データベースフィルターの入力物として使用され得る。
ある特定の電子システムおよびコンピューター実装方法
ある特定の電子システムおよびコンピューター実装方法
【0073】
本明細書で提供される方法およびシステムの一部の実施形態は、遺伝的変異データを分析するための電子システムを含む。一部のかかる実施形態では、本明細書に記載されるデータベースフィルターおよび／または本明細書に記載される近接フィルターは、生殖系列バリアントを識別するために、遺伝的変異データに適用され得る。
【0074】
一部の実施形態は、プロセッサー上で実行される、腫瘍細胞を含む生体試料由来の配列データから複数のバリアントを識別するように適合されたインフォマティクスモジュールであって、複数のバリアントが、体細胞バリアントおよび生殖系列バリアントを含む、インフォマティクスモジュールを含み得る。
【0075】
一部の実施形態は、複数のバリアントから生殖系列バリアントを除去するように適合されたデータベースフィルターモジュールであって、生殖系列バリアントが各々、閾値対立遺伝子カウントよりも大きいまたはそれと等しい、バリアントの参照セットにおける対立遺伝子カウントを有する、データベースフィルターモジュールを含む。一部の実施形態では、閾値対立遺伝子カウントは、２、３、４、５、６、７、８、９、１０、１１、１２、１３、１４、１５、１６、１７、１８、１９および２０よりも大きいまたはそれと等しい場合がある。
【0076】
一部の実施形態は、複数のバリアントから生殖系列バリアントを除去するように適合された近接フィルターモジュールを含む。一部の実施形態では、近接フィルターモジュールは、複数のビンを返すように適合されたビニングサブモジュールであって、各ビンが、ゲノムの同じ領域中に位置する複数のバリアントのバリアントを含む、ビニングサブモジュールを含み得る。一部の実施形態では、複数のバリアントのバリアントは、ゲノムの同じ領域中に位置するバリアントが、同じビン中に選別またはビニングされるように、複数のビン中に選別またはビニングされ得る。一部の実施形態では、ゲノムの同じ領域は、同じ染色体内、染色体の同じ腕内、同じ染色体サイトバンド内にあり得る。一部の実施形態では、ゲノムの同じ領域は、同じ連続する１００Ｍｂ、５０Ｍｂ、４０Ｍｂ、３０Ｍｂ、２０Ｍｂ、１０Ｍｂ、５Ｍｂ、１Ｍｂ内、または上述の数のうち任意の２つの間の任意の範囲内であり得る。
【0077】
一部の実施形態では、近接フィルターモジュールは、複数のバリアントにおけるデータベースバリアントを返すように適合された識別サブモジュールであって、データベースバリアントが、バリアントの参照セット中に存在する、識別サブモジュールを含み得る。
【0078】
一部の実施形態では、近接フィルターモジュールは、複数のバリアントから生殖系列バリアントを除去するように適合された除去サブモジュールであって、生殖系列バリアントが各々、生殖系列バリアントと同じビン中の少なくとも１つのデータベースバリアントの対立遺伝子頻度の近似範囲内の対立遺伝子頻度を有する、除去サブモジュールを含み得る。一部の実施形態では、近接フィルターは、生殖系列バリアントとして識別されていないバリアントについての対立遺伝子頻度の近似範囲を決定することを含む。一部の実施形態では、近似範囲は、バリアントの対立遺伝子頻度から０．０１、０．０２、０．０３、０．０４、０．０５、０．０６、０．０７、０．０８、０．０９、または上述の数のうち任意の２つの間の範囲内の任意の数の最大値および最小値を有する範囲である。一部の実施形態では、近似範囲は、バリアントの対立遺伝子頻度の二項分布から２標準偏差の最大値および最小値を有する範囲である。一部の実施形態では、近似範囲は、バリアントの対立遺伝子頻度の上および下、０．０５、またはバリアントの対立遺伝子頻度の二項分布から２（ｎ）標準偏差、のうち高い方である。
【0079】
一部の実施形態では、バリアントが、バリアントと同じビン中の１つまたは複数の識別された生殖系列バリアントの近似範囲内の対立遺伝子頻度を有する場合、バリアントは、生殖系列バリアントとして識別され得る。一部の実施形態では、バリアントが、バリアントと同じビン中の１、２、３、４、５、６、７、８、９または１０個よりも多い識別された生殖系列バリアントの近似範囲内の対立遺伝子頻度を有する場合、バリアントは、生殖系列バリアントとして識別され得る。一部の実施形態では、除去サブモジュールは、閾値頻度よりも大きいまたはそれと等しい対立遺伝子頻度を有するバリアントを除去するように適合される。一部のかかる実施形態では、０．７、０．８、０．９または１．０よりも大きいまたはそれと等しい対立遺伝子頻度を有するバリアントは、生殖系列バリアントとして識別され得る。一部の実施形態では、除去サブモジュールは、複数のバリアントから、バリアントの参照セット中に存在するデータベースバリアントを除去するように適合される。
【0080】
本明細書で提供される一部の実施形態は、複数のバリアントにおける体細胞バリアントを識別するためのコンピューター実装方法を含む。一部のかかる実施形態は、腫瘍細胞を含む生体試料由来の配列データから複数のバリアントを受け取るステップであって、複数のバリアントが、体細胞バリアントおよび生殖系列バリアントを含み得る、ステップを含み得る。
【0081】
一部の実施形態は、複数のバリアントにデータベースフィルターを適用するステップを含む。一部のかかる実施形態は、複数のバリアントについてドキュメントのインデックスを作成する工程、インデックスを用いてバリアントの参照セットを検索して、インデックスにおける生殖系列バリアントを識別する工程を含む。一部の実施形態では、生殖系列バリアントは各々、閾値対立遺伝子カウントよりも大きいまたはそれと等しい、バリアントの参照セットにおける対立遺伝子カウントを有する。一部の実施形態では、閾値対立遺伝子カウントは、２、３、４、５、６、７、８、９、１０、１１、１２、１３、１４、１５、１６、１７、１８、１９および２０よりも大きいまたはそれと等しい場合がある。一部の実施形態は、インデックスから識別された生殖系列バリアントを除去して、第１のフィルタリングされたバリアントのインデックスを作成する工程も含む。
【0082】
一部の実施形態は、第１のフィルタリングされたバリアントのインデックスに近接フィルターを適用するステップを含む。一部のかかる実施形態は、ゲノムの異なる領域について複数のビンを作成する工程を含む。一部の実施形態は、第１のフィルタリングされたバリアントのインデックスのバリアントをビニングする工程であって、ゲノムの同じ領域中に位置するバリアントが、同じビン中にビニングされる、工程を含む。一部の実施形態では、ゲノムの同じ領域は、同じ染色体内、染色体の同じ腕内、同じ染色体サイトバンド内にあり得る。一部の実施形態では、ゲノムの同じ領域は、同じ連続する１００Ｍｂ、５０Ｍｂ、４０Ｍｂ、３０Ｍｂ、２０Ｍｂ、１０Ｍｂ、５Ｍｂ、１Ｍｂ内、または上述の数のうち任意の２つの間の任意の範囲内であり得る。
【0083】
一部の実施形態は、第１のフィルタリングされたバリアントのインデックスを用いてバリアントの参照セットを検索して、第１のフィルタリングされたバリアントのインデックスにおけるデータベースバリアントを識別する工程を含む。
【0084】
一部の実施形態は、生殖系列バリアントを識別することによって、第１のフィルタリングされたバリアントのインデックスから生殖系列バリアントのインデックスを生成する工程を含む。一部の実施形態では、生殖系列バリアントは各々、第２の生殖系列バリアントと同じビン中の少なくとも１つのデータベースバリアントの対立遺伝子頻度の近似範囲内の対立遺伝子頻度を有する。一部の実施形態では、近似範囲は、バリアントの対立遺伝子頻度から０．０１、０．０２、０．０３、０．０４、０．０５、０．０６、０．０７、０．０８、０．０９、または上述の数のうち任意の２つの間の範囲内の任意の数の最大値および最小値を有する範囲である。一部の実施形態では、近似範囲は、バリアントの対立遺伝子頻度の二項分布から２標準偏差の最大値および最小値を有する範囲である。一部の実施形態では、近似範囲は、バリアントの対立遺伝子頻度の上および下、０．０５、またはバリアントの対立遺伝子頻度の二項分布から２（ｎ）標準偏差のうち高い方である。
【0085】
一部の実施形態では、バリアントが、バリアントと同じビン中の１つまたは複数の識別された生殖系列バリアントの近似範囲内の対立遺伝子頻度を有する場合、バリアントは、生殖系列バリアントとして識別され得る。一部の実施形態では、バリアントが、バリアントと同じビン中の１、２、３、４、５、６、７、８、９または１０個よりも多い識別された生殖系列バリアントの近似範囲内の対立遺伝子頻度を有する場合、バリアントは、生殖系列バリアントとして識別され得る。一部の実施形態では、生殖系列バリアントは、閾値頻度よりも大きいまたはそれと等しい対立遺伝子頻度を有するバリアントとして識別され得る。一部のかかる実施形態では、０．７、０．８、０．９または１．０よりも大きいまたはそれと等しい対立遺伝子頻度を有するバリアントは、生殖系列バリアントとして識別され得る。
【0086】
一部の実施形態は、第１のフィルタリングされたバリアントのインデックスから識別された生殖系列バリアントを除去して、体細胞バリアントのインデックスを作成し、それによって、複数のバリアントにおける体細胞バリアントを識別する工程を含む。一部の実施形態では、腫瘍由来のシーケンシングデータから得られた体細胞バリアントの数は、腫瘍の腫瘍突然変異量である。
処置の方法
処置の方法
【0087】
方法およびシステムの一部の実施形態は、腫瘍を処置する方法を含む。一部のかかる実施形態では、腫瘍中に存在する体細胞バリアントの数は、本明細書で提供される方法およびシステムによって決定され得る。例えば、配列データが腫瘍から得られ得、複数のバリアントが配列データから識別され得、生殖系列バリアントが識別され得、複数のバリアントから除去され得、それによって、複数のバリアントにおける体細胞バリアントを識別し得る。一部の実施形態では、生殖系列バリアントは、データベースフィルターおよび／または近接フィルターのうち１つまたは複数を適用することによって、識別され得、複数のバリアントから除去され得、それによって、フィルターのうち１つまたは複数を適用することによって除去されなかった体細胞バリアントを識別し得る。一部の実施形態では、腫瘍由来のシーケンシングデータから得られた体細胞バリアントの数は、腫瘍の腫瘍突然変異量である。一部の実施形態では、腫瘍突然変異量は、ゲノム領域当たりの体細胞バリアントの平均数、例えば、５０ｋｂ、１００ｋｂ、１Ｍｂ、１０Ｍｂ、１００Ｍｂなど当たりの突然変異として計算される。腫瘍突然変異量は、ゲノム全体またはその一部分をシーケンシングすることによって、サンプリングされ得る。例えば、ゲノムの一部分は、１つまたは複数の目的のゲノム領域、例えば、腫瘍遺伝子パネル、完全エクソーム、部分エクソームなどについて富化することによって、シーケンシングされ得る。
【0088】
腫瘍を処置する一部の実施形態は、腫瘍突然変異量閾値よりも大きいまたはそれと等しい腫瘍突然変異量を有する腫瘍を決定するステップ、および腫瘍を有効量の治療剤と接触させるステップを含み得る。一部の実施形態は、腫瘍を有する対象を処置するステップを含み、ＴＭＢ閾値よりも大きいまたはそれと等しい腫瘍突然変異量を有する腫瘍を決定するステップ、および対象に有効量の治療剤を投与するステップを含み得る。一部の実施形態では、腫瘍突然変異量閾値は、２、３、４、５、６、７、８、９、１０、１５、２０、３０、４０、５０、６０、７０、８０、９０、１００、２００、３００、４００、５００、６００、７００、８００、９００、１０００または上述の数のうち任意の２つの間の範囲中の任意の数であり得る。治療剤の例には、化学療法剤が含まれる。一部の実施形態では、治療剤には、チェックポイント阻害剤が含まれ得る。チェックポイント阻害剤の例には、ＣＴＬＡ－４阻害剤、ＰＤ－１阻害剤およびＰＤ－Ｌ１阻害剤が含まれる。一部の実施形態では、チェックポイント阻害剤には、イピリムマブ、ニボルマブ、ペムブロリズマブ、スパルタリズマブ、アテゾリズマブ、アベルマブおよびデュルバルマブが含まれ得る。腫瘍の例には、結腸直腸腫瘍、肺腫瘍、子宮内膜腫瘍、子宮腫瘍、胃腫瘍、黒色腫、乳房腫瘍、膵腫瘍、腎臓腫瘍、膀胱腫瘍および脳腫瘍が含まれる。本明細書に含まれる方法およびシステムで処置され得るがんの追加の例は、その全体が本明細書に参照によって明示的に組み込まれる、米国特許出願公開第２０１８０２１８７８９号に列挙されている。
試料
試料
【0089】
一部の実施形態は、生体試料から配列データを得るステップを含む。一部の実施形態では、生体試料は、腫瘍細胞を含み得る。一部の実施形態では、生体試料には、血清試料、大便試料、血液試料および腫瘍試料が含まれ得る。一部の実施形態では、生体試料は、固定される。
【0090】
一部の実施形態では、対象は、生体試料を提供し得る。生体試料は、対象によって産生される任意の物質であり得る。一般に、生体試料は、対象から採取される任意の組織または対象によって産生される任意の物質である。生体試料の例には、血液、血漿、唾液、脳脊髄液（ＣＳＦ）、頬組織、尿、糞便、皮膚、毛髪、臓器組織が含まれ得る。一部の実施形態では、生体試料は、固形腫瘍または固形腫瘍の生検である。一部の実施形態では、生体試料は、ホルマリン固定パラフィン包埋（ＦＦＰＥ）組織試料である。生体試料は、核酸を含む任意の生体試料であり得る。生体試料は、対象に由来し得る。対象は、哺乳動物、爬虫類、両生類、鳥類または魚類であり得る。哺乳動物の例には、ヒト、類人猿、オランウータン、サル、チンパンジー、ウシ、ブタ、ウマ、げっ歯類、トリ、爬虫類、イヌ、ネコ、イルカまたは他の動物が含まれる。爬虫類の例には、トカゲ、ヘビ、アリゲーター、カメ（turtle）、クロコダイル、イグアナおよびカメ（tortoise）が含まれる。両生類の例には、ヒキガエル、カエル、イモリおよびサンショウウオが含まれる。鳥類の例には、ニワトリ、カモ、ガチョウ、ペンギン、ダチョウ、ツノメドリおよびフクロウが含まれる。魚類の例には、ナマズ、ウナギ、サメ、金魚およびメカジキが含まれる。一部の実施形態では、対象はヒトである。
ある特定のシステムおよび方法
ある特定のシステムおよび方法
【0091】
一部の実施形態は、本明細書に記載される方法を実施するための、コンピューターベースのシステムおよびコンピューター実装方法を含む。一部の実施形態では、システムは、試料中のバリアント、例えば、生殖系列バリアントおよび／または体細胞バリアントの存在または非存在を決定し報告するために利用され得る。システムは、１つまたは複数のクライアントコンポーネントを含み得る。１つまたは複数のクライアントコンポーネントは、ユーザーインターフェースを含み得る。システムは、１つまたは複数のサーバーコンポーネントを含み得る。サーバーコンポーネントは、１つまたは複数のメモリロケーションを含み得る。１つまたは複数のメモリロケーションは、データ入力物を受け取るように構成され得る。データ入力物は、シーケンシングデータを含み得る。シーケンシングデータは、対象由来の核酸試料から生成され得る。システムは、１つまたは複数のコンピュータープロセッサーをさらに含み得る。１つまたは複数のコンピュータープロセッサーは、１つまたは複数のメモリロケーションに操作可能に連結され得る。１つまたは複数のコンピュータープロセッサーは、シーケンシングデータを参照配列にマッピングするようにプログラミングされ得る。１つまたは複数のコンピュータープロセッサーは、シーケンシングデータから複数のバリアントの存在または非存在を決定するようにさらにプログラミングされ得る。１つまたは複数のコンピュータープロセッサーは、少なくとも１つのフィルターを遺伝的変異体に適用して、生殖系列バリアントを識別するように、さらにプログラミングされ得る。フィルターの例には、データベースフィルターおよび近接フィルターが含まれる。１つまたは複数のコンピュータープロセッサーは、識別されたバリアントのインデックスから識別された生殖系列バリアントを除去するようにさらにプログラミングされ得る。１つまたは複数のコンピュータープロセッサーは、スクリーン上でのディスプレイのための出力物を生成するようにさらにプログラミングされ得る。出力物は、複数のバリアントにおける生殖系列バリアントおよび／または体細胞バリアントを識別する１つまたは複数の報告を含み得る。
【0092】
方法およびシステムの一部の実施形態は、１つまたは複数のクライアントコンポーネントを含み得る。１つまたは複数のクライアントコンポーネントは、１つもしくは複数のソフトウェアコンポーネント、１つもしくは複数のハードウェアコンポーネント、またはそれらの組合せを含み得る。１つまたは複数のクライアントコンポーネントは、１つまたは複数のサーバーコンポーネントを介して１つまたは複数のサービスにアクセスできる。１つまたは複数のサービスは、ネットワークを介して、１つまたは複数のクライアントコンポーネントによってアクセスされ得る。「サービス」は、任意の製品、方法、機能、またはシステムの使用を指すために、本明細書で使用される。例えば、ユーザーは、遺伝子検査を注文することができる。注文は、システムの１つまたは複数のクライアントコンポーネントを介して行われ得、要求は、ネットワークを介して、システムの１つまたは複数のサーバーコンポーネントに送信され得る。ネットワークは、インターネット、インターネットおよび／もしくはエクストラネット、またはインターネットと通信するイントラネットおよび／もしくはエクストラネットであり得る。ネットワークは、一部の例では、電気通信および／またはデータネットワークである。ネットワークは、分散型コンピューティング、例えば、クラウドコンピューティングを可能にできる、１つまたは複数のコンピューターサーバーを含み得る。一部の例では、コンピューターシステムの助けを借りたネットワークは、コンピューターシステムに連結されたデバイスがクライアントまたはサーバーとして挙動することを可能にし得るピアツーピアネットワークを実装し得る。
【0093】
システムの一部の実施形態は、１つまたは複数のメモリロケーション、例えば、ランダムアクセスメモリ、リードオンリーメモリ、フラッシュメモリ；電子記憶ユニット、例えば、ハードディスク；通信インターフェース、例えば、１つまたは複数の他のシステムとの通信のためのネットワークアダプタ、ならびに周辺デバイス、例えば、キャッシュ、他のメモリ、データ記憶および／または電子ディスプレイアダプタを含み得る。メモリ、記憶ユニット、インターフェースおよび周辺デバイスは、通信バスを介して、マザーボードなどのＣＰＵと通信する。記憶ユニットは、データを記憶するためのデータ記憶ユニットまたはデータリポジトリであり得る。一例では、１つまたは複数のメモリロケーションは、受け取ったシーケンシングデータを記憶できる。
【0094】
方法およびシステムの一部の実施形態は、１つまたは複数のコンピュータープロセッサーを含み得る。１つまたは複数のコンピュータープロセッサーは、例えば、記憶されたシーケンシングデータにアクセスするために、１つまたは複数のメモリロケーションに操作可能に連結され得る。１つまたは複数のコンピュータープロセッサーは、本明細書に記載される方法を実施するために、機械実行可能コードを実装できる。例えば、１つまたは複数のコンピュータープロセッサーは、シーケンシングデータ入力物を参照配列にマッピングするためならびに／または生殖系列バリアントおよび／もしくは体細胞バリアントを識別するために、機械可読コードを実行できる。
【0095】
本明細書で提供される方法およびシステムの一部の実施形態は、機械実行可能コードまたは機械可読コードを含み得る。一部のかかる実施形態では、機械実行可能コードまたは機械可読コードは、ソフトウェアの形態で提供され得る。使用の間に、コードは、プロセッサーによって実行され得る。一部の場合には、コードは、記憶ユニットから検索され得、プロセッサーによる即座のアクセスのために、メモリ上に記憶され得る。一部の実施形態では、電子記憶ユニットは、除外され得、機械実行可能命令は、メモリ上に記憶される。コードは、コードを実行するように適合されたプロセッサーを有する機械との使用のためにプリコンパイルおよび構成され得、ランタイムの間にコンパイルされ得、またはランタイムの間に解読され得る。コードは、プリコンパイルされた、アズコンパイルされた（as-compiled）または解読された様式でコードが実行されるのを可能にするように選択され得るプログラミング言語で提供され得る。
【0096】
本明細書で提供されるシステムおよび方法の一部の実施形態、例えば、コンピューターシステムは、プログラミングで具体化され得る。テクノロジーの種々の態様は、典型的には、ある型の機械可読媒体上に保持されるまたはある型の機械可読媒体で具体化される機械（またはプロセッサー）実行可能コードおよび／または関連データの形態での、「製品」または「製造品」であると考えられ得る。機械実行可能コードは、電子記憶ユニット、例えば、メモリまたはハードディスク上に記憶され得る。「記憶」型媒体は、ソフトウェアプログラミングのためにいつでも一時的でない（non-transitory）記憶を提供し得る、コンピューターの有形メモリ、プロセッサーなどのうちいずれかもしくは全て、またはその関連モジュール、例えば、種々の半導体メモリ、テープドライブ、ディスクドライブなどを含み得る。ソフトウェアの全てまたは部分は、時折、インターネットまたは種々の他の電気通信ネットワークを介して通信され得る。例えば、かかる通信は、１つのコンピューターまたはプロセッサーから別のコンピューターまたはプロセッサーへの、例えば、管理サーバーまたはホストコンピューターからアプリケーションサーバーのコンピュータープラットフォームへの、ソフトウェアのロードを可能にし得る。したがって、ソフトウェア要素を有し得る別の型の媒体には、例えば、ローカルデバイス間の物理インターフェースを横断して、有線および光学地上通信線ネットワークを介して、ならびに種々のエアリンク（air-link）を通じて使用される、光波、電波および電磁波が含まれる。かかる波を伝える物理要素、例えば、有線または無線リンク、光リンクなどもまた、ソフトウェアを有する媒体とみなされ得る。本明細書で使用される場合、一時的でない有形「記憶」媒体に限定されない限り、コンピューターまたは機械「可読媒体」などの用語は、実行のために命令をプロセッサーに提供することに関与する任意の媒体を指す。
【0097】
本明細書に開示される方法およびシステムの一部の実施形態は、１つもしくは複数の電子ディスプレイを含み得る、またはそれと通信し得る。電子ディスプレイは、コンピューターシステムの一部であり得る、または直接、もしくはネットワークを介して、コンピューターシステムに連結され得る。コンピューターシステムは、本明細書に開示される種々の特性および機能を提供するためのユーザーインターフェース（ＵＩ）を含み得る。ＵＩの例には、限定なしに、グラフィカルユーザーインターフェース（ＧＵＩ）およびウェブベースのユーザーインターフェースが含まれる。ＵＩは、ユーザーが本明細書に記載される方法およびシステムを利用できるインタラクティブツールを提供できる。例として、本明細書で想定されるＵＩは、医療関係者が遺伝子検査を注文でき、検査すべき遺伝的変異体のリストをカスタマイズでき、生物医学的報告を受け取りそれを見ることができる、ウェブベースのツールであり得る。
【0098】
本明細書に開示される方法およびシステムの一部の実施形態は、１つもしくは複数のデータベース、１つもしくは複数のアッセイ、１つもしくは複数のデータもしくは結果、１つもしくは複数のアッセイに基づくまたはそれに由来する１つもしくは複数の出力物、１つもしくは複数のデータもしくは結果に基づくまたはそれに由来する１つもしくは複数の出力物、あるいはそれらの組合せからのデータおよび／または情報に基づく、生物医学的データベース、ゲノムデータベース、生物医学的報告、疾患報告、症例対照分析およびレアバリアント発見分析を含み得る。
【実施例】
【0099】
（実施例１）
試料比較によって体細胞バリアントを識別する
配列データを、個体由来の腫瘍試料および正常試料について得た。バリアントを、配列データにおいて識別した。腫瘍試料中の生殖系列バリアントを、腫瘍試料中に存在するバリアントおよび正常試料中に存在しないバリアントを比較することによって識別した。図２Ａは、体細胞バリアント（黒色で塗りつぶされた丸）および生殖系列バリアント（灰色で塗りつぶされた丸）による、各バリアントの染色体位置に従う異なるバリアントについてのバリアント対立遺伝子頻度（ＶＡＦ）を示すグラフである。この方法は、個体由来の２つの試料を必要とした。
（実施例２）
バリアントのデータベースフィルタリング
試料比較によって体細胞バリアントを識別する
配列データを、個体由来の腫瘍試料および正常試料について得た。バリアントを、配列データにおいて識別した。腫瘍試料中の生殖系列バリアントを、腫瘍試料中に存在するバリアントおよび正常試料中に存在しないバリアントを比較することによって識別した。図２Ａは、体細胞バリアント（黒色で塗りつぶされた丸）および生殖系列バリアント（灰色で塗りつぶされた丸）による、各バリアントの染色体位置に従う異なるバリアントについてのバリアント対立遺伝子頻度（ＶＡＦ）を示すグラフである。この方法は、個体由来の２つの試料を必要とした。
（実施例２）
バリアントのデータベースフィルタリング
【0100】
配列データを、実施例１の腫瘍試料のみについて得た。バリアントを、配列データにおいて識別した。簡潔に述べると、バリアント呼び出しパイプラインから呼び出されたバリアントを、アノテーションツールＮｉｒｖａｎａ（Ｉｌｌｕｍｉｎａ、Ｓａｎ Ｄｉｅｇｏ）を使用してアノテーションした。Ｎｉｒｖａｎａは、ゲノムバリアント、例えば、単一ヌクレオチドバリアント、マルチヌクレオチドバリアント、挿入、欠失、コピー数バリアントの臨床グレードのアノテーションを提供した。Ｎｉｒｖａｎａへの入力物は、バリアントコールフォーマット（ＶＣＦ）であり、出力物は、全てのアノテーションおよび試料情報の構造化ＪＳＯＮ提示であった。
【0101】
識別されたバリアントについて、総対立遺伝子カウントを、バリアント対立遺伝子頻度およびカバレッジと共に、ゲノム集約データベース（ｇｎｏｍＡＤ）エクソーム、ｇｎｏｍＡＤゲノムおよび１０００ゲノムデータベース中の所与のバリアントについて解析した。これらの総対立遺伝子カウントは、異なる下位集団を横断するデータベース内の観察の総数を示した。各バリアントについて、３つ全てのデータベースにおいて観察された最大対立遺伝子カウントを、ゲノムデータベースと比較したそのより大きい試料サイズを利用しながら、エクソームデータベース中のカバーされなかった領域を考慮に入れるために、取り込んだ。フィルタリング戦略は、潜在的な生殖系列バリアントとして、１０よりも大きいまたはそれと等しい最大対立遺伝子カウントを有するバリアントにマークした。所与のバリアントについて、データベースにおける１０の対立遺伝子カウントは、それらが全てホモ接合性であった場合には少なくとも５つの試料において、またはそれらが全てヘテロ接合性であった場合には最大１０個の試料において、それが観察されたことを意味している。図２Ｂは、フィルター決定された体細胞バリアント（黒色で塗りつぶされた丸）およびフィルター決定された生殖系列バリアント（灰色で塗りつぶされた丸）による、各バリアントの染色体位置に従う種々のバリアントについてのバリアント対立遺伝子頻度（ＶＡＦ）を示すグラフである。これは、データベースフィルタリングのみが、バリアントを誤って呼び出したことを実証した。
（実施例３）
バリアントの近接フィルタリング
（実施例３）
バリアントの近接フィルタリング
【0102】
配列データを、個体由来の腫瘍試料のみについて得た。バリアントを、配列データにおいて識別した。バリアントに実施例２のデータベースフィルターを適用した。近接フィルターを使用して、データベース中に見出されなかったバリアントをさらにフィルタリングして除いた。
【0103】
近接フィルターは、位置的に密に近接したデータベースフィルタリングされたバリアントの情報を使用した。データベース中に見出されず、０．９よりも低い対立遺伝子頻度を有した所与のバリアントについて、同じ染色体上のバリアントを、フィルタリングされていないバリアントのバリアント対立遺伝子頻度の所与の範囲内で検索した。９０％よりも大きい対立遺伝子頻度を有するバリアントは、いずれのさらなる処理もなしに、生殖系列としてマークした。範囲を、所与のバリアントについての支持的証拠が二項プロセスによって生成されると仮定して、０．０５および二項分布の２標準偏差の最大値として決定した。例えば、フィルタリングされていないバリアントが、１００のカバレッジと共に０．２の対立遺伝子頻度を有した場合、範囲は、０．０５と２＊ｓｑｒｔ（１００＊０．２＊（１－０．２））／１００＝０．０８との間の最大値であり、それは、０．０８であった。これは、両方の方向での０．０８の範囲へと変換でき、０．１２と０．２８との間の対立遺伝子頻度を有する全てのバリアントを、同じ染色体から検索した。引き続いて、固定された閾値を上回る検索されたバリアントの数をチェックし、これを５に設定した。必要とされるバリアントの数に達した時点で、本発明者らは、バリアントの０．９５に設定したかなりの割合がデータベースフィルターによってフィルタリングされたかどうかをチェックした。両方の条件を満たした場合に、バリアントを近接フィルターによってマークした。言い換えれば、バリアントが、データベース中に見出された、対立遺伝子頻度空間中の十分な数のバリアントによって取り囲まれた場合、これも生殖系列バリアントとみなした。このフィルターは、およそ５０％または１００％の期待されたバリアント対立遺伝子頻度を有する正常領域中の、および対立遺伝子頻度分布がシフトし得るコピー数バリアント領域中の、生殖系列バリアントを除去した。
【0104】
図３（左パネル）は、データベースフィルターのみでフィルタリングした、フィルター決定された体細胞バリアント（黒色で塗りつぶされた丸）およびフィルター決定された生殖系列バリアント（灰色で塗りつぶされた丸）による、各バリアントについて第１～７染色体についての染色体位置に従う種々のバリアントについてのバリアント対立遺伝子頻度（ＶＡＦ）を示すグラフである。図３（右パネル）は、特定のフィルター決定された体細胞バリアント（黒色の丸）が選択され、いくつかのフィルター決定された生殖系列バリアント（灰色の丸）を包含するバリアントから範囲が引き出された、第７染色体上に位置するバリアントについての拡大図である。選択されたフィルター決定された体細胞バリアント（黒色の丸）が生殖系列バリアントとして呼び出されるべきであるという決定は、選択されたバリアントの対立遺伝子頻度の、ある特定の数のすでに識別された生殖系列バリアントの対立遺伝子頻度への近接に基づいてなされ得る。
【0105】
図４Ａは、データベースフィルターのみでフィルタリングした、フィルター決定された体細胞バリアント（黒色で塗りつぶされた丸）およびフィルター決定された生殖系列バリアント（灰色で塗りつぶされた丸）による、各バリアントの染色体位置に従う種々のバリアントについてのバリアント対立遺伝子頻度（ＶＡＦ）を示すグラフである。図４Ｂは、データベースフィルターのみおよび近接フィルターでフィルタリングした、フィルター決定された体細胞バリアント（黒色で塗りつぶされた丸）およびフィルター決定された生殖系列バリアント（灰色で塗りつぶされた丸）による、各バリアントの染色体位置に従う種々のバリアントについてのバリアント対立遺伝子頻度（ＶＡＦ）を示すグラフである。図４Ｂは、図４Ａにおいて体細胞バリアントとして示されたある特定の推定偽陽性が、図４Ｂにおいて生殖系列バリアントとして識別されたことを示している。例えば、約０．４および０．３の対立遺伝子頻度を有する、第７染色体上に位置する識別された体細胞バリアント（図４Ａ）は、近接フィルターを適用した場合、生殖系列バリアントとして識別された（図４Ｂ）。
（実施例４）
標的化シーケンシングを用いて腫瘍突然変異量を測定する
（実施例４）
標的化シーケンシングを用いて腫瘍突然変異量を測定する
【0106】
この実施例は、ホルマリン固定パラフィン包埋（ＦＦＰＥ）腫瘍試料中の腫瘍突然変異量（ＴＭＢ）を測定するための標的化次世代シーケンシングアッセイに関する。図５は、アッセイのための、例となるワークフローを示す。配列データを、１．３３Ｍｂのエクソンサイズを有する１．９４Ｍｂのパネルサイズ中の５２３個の遺伝子について、腫瘍試料から得た。シーケンシングを、固有分子識別子（unique molecular identifier）（ＵＭＩ）を用い、Ｉｌｌｕｍｉｎａ ＮｅｘｔＳｅｑ（商標）５００／５５０プラットフォームを使用して実施した。データ分析を、５％のバリアント対立遺伝子頻度（ＶＡＦ）においてバリアントを検出するためのパイプラインを使用して実施した。技術的ノイズ除去のために、ＵＭＩからの情報、および試料特異的エラープロファイルを利用するバリアント呼び出しアルゴリズムを使用して、異なるＦＦＰＥ品質の試料にわたる均一なバリアント呼び出しパフォーマンスを確実にした。ＴＭＢ計算から生殖系列バリアントを正確に除去するために、大規模な公開データベースからの情報を、各バリアントの測定されたカバレッジおよびバリアント対立遺伝子頻度と統合する、上述の実施例のデータベースフィルターおよび近接フィルターと実質的に類似したハイブリッド戦略を使用した。
【0107】
簡潔に述べると、配列データを得、参照とアラインさせ、バリアントを識別した。データベースフィルターおよび近接フィルターを使用して、識別されたバリアントから生殖系列バリアントをフィルタリングし、ＴＭＢを、図６に示されるパイプラインと実質的に類似したワークフローで計算した。合計１７０対の腫瘍－正常試料を分析して、生殖系列フィルタリングおよびＴＭＢパフォーマンスを評価した（表１）。１０８試料対のサブセットは、全エクソームシーケンシング（ＷＥＳ）でも分析した。
表1
表1
【表1】
【0108】
技術的ノイズ除去のために、正常ＦＦＰＥ試料のコレクション（Ｎ＝１７６）中の偽陽性バリアントの数を評価した。試料品質とは無関係に、試料１つ当たり平均０．６３の偽陽性が観察され（Ｒ２＝０．００１）、９２．６％の試料が、２以下の偽陽性バリアント（ＶＡＦ＜２０％）を含んだ。さらに、５％に近いバリアントを有し、９８．７％の感度を達成した、ＦＦＰＥおよび細胞系混合試料のセットを検査した。
【0109】
生殖系列フィルタリングパフォーマンスを、表１に記載される１７０の腫瘍／正常試料対を使用して評価した。小さいバリアント（ＳＮＶ、挿入／欠失）生殖系列バリアントフィルタリングにおいて、試料１つ当たり平均して１．３よりも少ない生殖系列バリアントを残す、９９．７％を上回る全体的フィルタリング率に達した。近接フィルタリングの追加は、体細胞突然変異に対する最小限の影響を有する一方で、偽陽性の数を顕著に低減させた。図７は、データベースのみ（約３生殖系列残留／Ｍｂにグラフのピーク）およびハイブリッド戦略（約０生殖系列残留／Ｍｂにグラフのピーク）によるフィルタリング後の残存する生殖系列バリアントカウントの分布を示す。
【0110】
ＴＭＢ再現性を、３つの演算子を横断して、４つの細胞系および４つのＦＦＰＥ試料を含む８つの異なる試料において評価した。各試料の平均および標準偏差（ＳＤ）を計算した。表２は、各々１２の複製にわたる、４つの細胞系および４つのＦＦＰＥ試料において評価したＴＭＢ再現性を列挙する。
表2
表2
【表2】
【0111】
合わせると、腫瘍のみアッセイによって生成されたＴＭＢ測定値は、腫瘍／正常アッセイ対試料から生成された推定値と高度に相関した（Ｒ２＝０．９９３、Ｎ＝１６９、ＴＭＢ＜２００の試料のみ）。腫瘍のみアッセイのＴＭＢ推定値は、全エクソームシーケンシングを通じて得られたＴＭＢ値とも、高い相関を示した（Ｒ２＝０．９３１、Ｎ＝１０５、ＷＥＳ ＴＭＢ＜１００の試料のみ）。図８Ａは、腫瘍のみアッセイと腫瘍／正常アッセイとの間のＴＭＢ比較を示す。図８Ｂは、腫瘍のみアッセイとＷＥＳ腫瘍－正常アッセイとの間のＴＭＢ比較を示す。
【0112】
最後に、１０のＴＭＢ閾値を用いて、９４．７４％の陽性一致率（ＰＰＡ）および９６．０８％の陰性一致率（ＮＰＡ）が実証された。全体的分類一致は、ＴＭＢ高試料およびＴＭＢ低試料の識別において、９５．３７％であった。表３は、ＴＭＢ分類パフォーマンスを列挙する。
表3
表3
【表3】
【0113】
上述の結果は、データベースフィルターおよび近接フィルターを用いた腫瘍のみについてのアッセイが、ＦＦＰＥ試料中のＴＭＢをロバストに測定する能力を実証した。さらに、ＴＭＢ推定値は、高い分類一致で、ＷＥＳベースの測定値との高レベルの相関を示した。
【0114】
用語「含む（comprising）」は、本明細書で使用される場合、「含む（including）」、「含む（containing）」または「～によって特徴付けられる」と同義であり、包含的またはオープンエンドであり、さらなる未列挙の要素または方法ステップを排除しない。
【0115】
上記記載は、本発明のいくつかの方法および材料を開示している。本発明は、方法および材料における改変、ならびに製造方法および装置における変更を受けることができる。かかる改変は、本開示の検討および本明細書に開示される発明の実施から、当業者に明らかとなる。結果として、本発明が本明細書に開示される特定の実施形態に限定されることは意図しないが、それが本発明の真の範囲および精神内に入る全ての改変および変更をカバーすることを意図する。
【0116】
公開および未公開の出願、特許ならびに文献参照が含まれるがこれらに限定されない、本明細書で引用される全ての参考文献は、それらの全体が参照によって本明細書に組み込まれ、これにより、本明細書の一部を構成する。参照によって組み込まれる刊行物および特許または特許出願が、本明細書中に含まれる開示と矛盾する限り、本明細書が、任意のかかる矛盾する資料に優先され、かつ／またはそれよりも優位であるものとする。
本発明は、例えば、以下の項目を提供する。
（項目１）
複数のバリアントにおける体細胞バリアントを識別するための方法であって、
（ａ）体細胞バリアントおよび生殖系列バリアントを含む複数のバリアントを得るステップ；
（ｂ）前記複数のバリアントにデータベースフィルターを適用するステップであって、
前記複数のバリアントにおける第１の生殖系列バリアントを決定する工程を含み、前記第１の生殖系列バリアントが各々、閾値対立遺伝子カウントよりも大きいまたはそれと等しい、バリアントの第１の参照セットにおける対立遺伝子カウントを有する、ステップ；
（ｃ）前記複数のバリアントに近接フィルターを適用するステップであって、
（ｉ）前記複数のバリアントのバリアントを複数のビン中にビニングする工程であって、ゲノムの同じ領域中に位置するバリアントが、同じビン中にビニングされる、工程、
（ｉｉ）前記複数のバリアントにおけるデータベースバリアントを決定する工程であって、データベースバリアントが、バリアントの第２の参照セット中に存在する、工程、および
（ｉｉｉ）前記複数のバリアントにおける第２の生殖系列バリアントを決定する工程であって、前記第２の生殖系列バリアントが各々、前記第２の生殖系列バリアントと同じビン中の少なくとも１つのデータベースバリアントの対立遺伝子頻度の近似範囲内の対立遺伝子頻度を有する、工程
を含む、ステップ；ならびに
（ｄ）前記複数のバリアントから識別された第１および第２の生殖系列バリアントを除去することによって、前記複数のバリアントにおける体細胞バリアントを決定するステップ
を含む、方法。
（項目２）
（ｂ）および（ｃ）が連続的に実施される、項目１に記載の方法。
（項目３）
（ｃ）が（ｂ）の前に実施される、項目１に記載の方法。
（項目４）
前記閾値対立遺伝子カウントが５である、項目１から３のいずれか一項に記載の方法。
（項目５）
前記閾値対立遺伝子カウントが１０である、項目４に記載の方法。
（項目６）
バリアントの前記第１および第２の参照セットが、同じ参照セットである、項目１から５のいずれか一項に記載の方法。
（項目７）
バリアントの前記第１または第２の参照セットが、複数の個体についてのバリアントのデータベースを含む、項目１から６のいずれか一項に記載の方法。
（項目８）
バリアントの前記第１または第２の参照セットが、ゲノム集約データベース（ｇｎｏｍＡＤ）および１０００ゲノムデータベースから選択される少なくとも１つのデータベースを含む、項目１から７のいずれか一項に記載の方法。
（項目９）
ゲノムの前記同じ領域が、同じ染色体内にある、項目１から８のいずれか一項に記載の方法。
（項目１０）
ゲノムの前記同じ領域が、同じ染色体腕内にある、項目１から９のいずれか一項に記載の方法。
（項目１１）
ゲノムの前記同じ領域が、同じ染色体サイトバンド内にある、項目１から１０のいずれか一項に記載の方法。
（項目１２）
ゲノムの前記同じ領域が、１０Ｍｂの領域内にある、項目１から１１のいずれか一項に記載の方法。
（項目１３）
近接フィルターを適用する前記ステップが、０．９よりも大きいまたはそれと等しい対立遺伝子頻度を有する第２の生殖系列バリアントを識別する工程をさらに含む、項目１から１２のいずれか一項に記載の方法。
（項目１４）
近接フィルターを適用する前記ステップが、前記複数のバリアントにおける第２の生殖系列バリアントを識別する工程をさらに含み、前記第２の生殖系列バリアントが、バリアントの前記第２の参照セット中に存在するデータベースバリアントである、項目１から１３のいずれか一項に記載の方法。
（項目１５）
前記近似範囲が、第２の生殖系列バリアントの前記対立遺伝子頻度から０．０５の最大値および最小値を有する範囲である、項目１から１４のいずれか一項に記載の方法。
（項目１６）
前記近似範囲が、第２の生殖系列バリアントの対立遺伝子頻度を中心として、第２の生殖系列バリアントの前記対立遺伝子頻度の二項分布から２標準偏差の最大値および最小値を有する範囲である、項目１から１５のいずれか一項に記載の方法。
（項目１７）
前記第２の生殖系列バリアントが、前記第２の生殖系列バリアントと同じビン中の少なくとも５つのデータベースバリアントの対立遺伝子頻度に近接した閾値内の対立遺伝子頻度を有する、項目１から１６のいずれか一項に記載の方法。
（項目１８）
前記第２の生殖系列バリアントが、前記第２の生殖系列バリアントと同じビン中の少なくとも１０個のデータベースバリアントの対立遺伝子頻度に近接した閾値内の対立遺伝子頻度を有する、項目１から１７のいずれか一項に記載の方法。
（項目１９）
（ａ）が、腫瘍細胞を含む生体試料から配列データを得る工程を含む、項目１から１８のいずれか一項に記載の方法。
（項目２０）
前記配列データを参照配列とアラインさせる工程、および前記配列データ中のバリアントを識別する工程をさらに含む、項目１９に記載の方法。
（項目２１）
腫瘍細胞を含む前記生体試料が、血清試料、大便試料、血液試料、腫瘍試料から選択される、項目１９または２０に記載の方法。
（項目２２）
前記腫瘍試料が固定される、項目２１に記載の方法。
（項目２３）
腫瘍の腫瘍突然変異量を決定する方法であって、
腫瘍細胞を含む生体試料から配列データを得るステップ；
前記配列データから複数のバリアントを決定するステップ；および
項目１から２２のいずれか一項に記載の方法に従って、複数のバリアントにおける体細胞バリアントの数を決定するステップであって、体細胞バリアントの前記数が、前記腫瘍の前記腫瘍突然変異量である、ステップ
を含む、方法。
（項目２４）
腫瘍を処置する方法であって、
項目２３に記載の方法に従って、１０体細胞バリアントよりも大きいまたはそれと等しい腫瘍突然変異量を有する腫瘍を決定するステップ、
有効量のチェックポイント阻害剤を投与することによって、前記腫瘍を処置するステップ
を含む、方法。
（項目２５）
前記腫瘍が、結腸直腸腫瘍、肺腫瘍、子宮内膜腫瘍、子宮腫瘍、胃腫瘍、黒色腫、乳房腫瘍、膵腫瘍、腎臓腫瘍、膀胱腫瘍および脳腫瘍からなる群から選択される、項目２４に記載の方法。
（項目２６）
前記チェックポイント阻害剤が、ＣＴＬＡ－４阻害剤、ＰＤ－１阻害剤およびＰＤ－Ｌ１阻害剤からなる群から選択される、項目２４または２５に記載の方法。
（項目２７）
前記チェックポイント阻害剤が、イピリムマブ、ニボルマブ、ペムブロリズマブ、スパルタリズマブ、アテゾリズマブ、アベルマブおよびデュルバルマブからなる群から選択される、項目２４から２６のいずれか一項に記載の方法。
（項目２８）
遺伝的変異データを分析するための電子システムであって、
プロセッサー上で実行される、腫瘍細胞を含む生体試料由来の配列データから複数のバリアントを識別するように適合されたインフォマティクスモジュールであって、前記複数のバリアントが、体細胞バリアントおよび生殖系列バリアントを含む、インフォマティクスモジュール；
前記複数のバリアントから第１の生殖系列バリアントを除去するように適合されたデータベースフィルターモジュールであって、前記第１の生殖系列バリアントが各々、閾値対立遺伝子カウントよりも大きいまたはそれと等しい、バリアントの第１の参照セットにおける対立遺伝子カウントを有する、データベースフィルターモジュール；
前記複数のバリアントから第２の生殖系列バリアントを除去するように適合された近接フィルターモジュールであって、前記近接フィルターモジュールが、
複数のビンを返すように適合されたビニングサブモジュールであって、各ビンが、ゲノムの同じ領域中に位置する前記複数のバリアントのバリアントを含む、ビニングサブモジュール、
前記複数のバリアントにおけるデータベースバリアントを返すように適合された識別サブモジュールであって、データベースバリアントが、バリアントの第２の参照セット中に存在する、識別サブモジュール、および
前記複数のバリアントから第２の生殖系列バリアントを除去するように適合された除去サブモジュールであって、前記第２の生殖系列バリアントが各々、前記第２の生殖系列バリアントと同じビン中の少なくとも１つのデータベースバリアントの対立遺伝子頻度の近似範囲内の対立遺伝子頻度を有する、除去サブモジュール
を含む、近接フィルターモジュール；ならびに
前記複数のバリアントから除去されなかったバリアントを返すように適合されたディスプレイモジュール
を含む、電子システム。
（項目２９）
インフォマティクスモジュールが、バリアントアノテーションツールを含む、項目２８に記載のシステム。
（項目３０）
前記閾値対立遺伝子カウントが５である、項目２８または２９に記載のシステム。
（項目３１）
前記閾値対立遺伝子カウントが１０である、項目３０に記載のシステム。
（項目３２）
バリアントの前記第１および第２の参照セットが、同じ参照セットである、項目２８から３１のいずれか一項に記載のシステム。
（項目３３）
バリアントの前記第１または第２の参照セットが、複数の個体についてのバリアントのデータベースを含む、項目２８から３２のいずれか一項に記載のシステム。
（項目３４）
バリアントの前記第１または第２の参照セットが、ゲノム集約データベース（ｇｎｏｍＡＤ）および１０００ゲノムデータベースから選択される少なくとも１つのデータベースを含む、項目２８から３３のいずれか一項に記載のシステム。
（項目３５）
ゲノムの前記同じ領域が、同じ染色体内にある、項目２８から３４のいずれか一項に記載のシステム。
（項目３６）
ゲノムの前記同じ領域が、同じ染色体腕内にある、項目２８から３５のいずれか一項に記載のシステム。
（項目３７）
ゲノムの前記同じ領域が、同じ染色体サイトバンド内にある、項目２８から３６のいずれか一項に記載のシステム。
（項目３８）
ゲノムの前記同じ領域が、１０Ｍｂの領域内にある、項目２８から３７のいずれか一項に記載のシステム。
（項目３９）
前記除去サブモジュールが、前記複数のバリアントから、０．９よりも大きいまたはそれと等しい対立遺伝子頻度を有するバリアントを除去するように適合される、項目２８から３８のいずれか一項に記載のシステム。
（項目４０）
前記除去サブモジュールが、前記複数のバリアントから、バリアントの前記第２の参照セット中に存在するデータベースバリアントを除去するように適合される、項目２８から３９のいずれか一項に記載のシステム。
（項目４１）
前記近似範囲が、第２の生殖系列バリアントの前記対立遺伝子頻度から０．０５の最大値および最小値を有する範囲である、項目２８から４０のいずれか一項に記載のシステム。
（項目４２）
前記近似範囲が、第２の生殖系列バリアントの対立遺伝子頻度を中心として、第２の生殖系列バリアントの前記対立遺伝子頻度の二項分布から２標準偏差の最大値および最小値を有する範囲である、項目２８から４１のいずれか一項に記載のシステム。
（項目４３）
前記第２の生殖系列バリアントが、前記第２の生殖系列バリアントと同じビン中の少なくとも５つのデータベースバリアントの対立遺伝子頻度に近接した閾値内の対立遺伝子頻度を有する、項目２８から４２のいずれか一項に記載のシステム。
（項目４４）
前記第２の生殖系列バリアントが、前記第２の生殖系列バリアントと同じビン中の少なくとも１０個のデータベースバリアントの対立遺伝子頻度に近接した閾値内の対立遺伝子頻度を有する、項目２８から４３のいずれか一項に記載のシステム。
（項目４５）
腫瘍細胞を含む前記生体試料が、血清試料、大便試料、血液試料、腫瘍試料から選択される、項目２８から４４のいずれか一項に記載のシステム。
（項目４６）
前記腫瘍試料が固定される、項目４５に記載のシステム。
（項目４７）
複数のバリアントにおける体細胞バリアントを識別するためのコンピューター実装方法であって、項目１から２２のいずれか一項に記載の方法を実施するステップを含む、方法。
（項目４８）
複数のバリアントにおける体細胞バリアントを識別するためのコンピューター実装方法であって、
（ａ）腫瘍細胞を含む生体試料由来の配列データから複数のバリアントを受け取るステップであって、前記複数のバリアントが、体細胞バリアントおよび生殖系列バリアントを含む、ステップ；
（ｂ）前記複数のバリアントにデータベースフィルターを適用するステップであって、
前記複数のバリアントについてドキュメントのインデックスを作成する工程、
前記インデックスを用いてバリアントの第１の参照セットを検索して、前記インデックスにおける第１の生殖系列バリアントを識別する工程であって、前記第１の生殖系列バリアントが各々、閾値対立遺伝子カウントよりも大きいまたはそれと等しい、バリアントの前記第１の参照セットにおける対立遺伝子カウントを有する、工程、および
前記インデックスから識別された第１の生殖系列バリアントを除去して、第１のフィルタリングされたバリアントのインデックスを作成する工程
を含む、ステップ；
（ｃ）第１のフィルタリングされたバリアントの前記インデックスに近接フィルターを適用するステップであって、
（ｉ）ゲノムの異なる領域について複数のビンを作成する工程、
（ｉｉ）第１のフィルタリングされたバリアントの前記インデックスのバリアントをビニングする工程であって、ゲノムの同じ領域中に位置するバリアントが、同じビン中にビニングされる、工程、
（ｉｉｉ）第１のフィルタリングされたバリアントの前記インデックスを用いてバリアントの第２の参照セットを検索して、第１のフィルタリングされたバリアントの前記インデックスにおけるデータベースバリアントを識別する工程、
（ｉｖ）第２の生殖系列バリアントを識別することによって、第１のフィルタリングされたバリアントの前記インデックスから第２の生殖系列バリアントのインデックスを生成する工程であって、前記第２の生殖系列バリアントが各々、前記第２の生殖系列バリアントと同じビン中の少なくとも１つのデータベースバリアントの対立遺伝子頻度の近似範囲内の対立遺伝子頻度を有する、工程、および
（ｖ）第１のフィルタリングされたバリアントの前記インデックスから識別された第２の生殖系列バリアントを除去して、体細胞バリアントのインデックスを作成し、それによって、前記複数のバリアントにおける体細胞バリアントを識別する工程
を含む、ステップ
を含む、方法。
（項目４９）
前記閾値対立遺伝子カウントが５である、項目４８に記載の方法。
（項目５０）
前記閾値対立遺伝子カウントが１０である、項目４９に記載の方法。
（項目５１）
バリアントの前記第１および第２の参照セットが、同じ参照セットである、項目４８から５０のいずれか一項に記載の方法。
（項目５２）
バリアントの前記第１または第２の参照セットが、複数の個体についてのバリアントのデータベースを含む、項目４８から５１のいずれか一項に記載の方法。
（項目５３）
バリアントの前記第１または第２の参照セットが、ゲノム集約データベース（ｇｎｏｍＡＤ）および１０００ゲノムデータベースから選択される少なくとも１つのデータベースを含む、項目４８から５２のいずれか一項に記載の方法。
（項目５４）
ゲノムの前記同じ領域が、同じ染色体内にある、項目４８から５３のいずれか一項に記載の方法。
（項目５５）
ゲノムの前記同じ領域が、同じ染色体腕内にある、項目４８から５４のいずれか一項に記載の方法。
（項目５６）
ゲノムの前記同じ領域が、同じ染色体サイトバンド内にある、項目４８から５５のいずれか一項に記載の方法。
（項目５７）
ゲノムの前記同じ領域が、１０Ｍｂの領域内にある、項目４８から５６のいずれか一項に記載の方法。
（項目５８）
第２のフィルタリングされたバリアントのインデックスを生成する前記工程が、０．９よりも大きいまたはそれと等しい対立遺伝子頻度を有する第２の生殖系列バリアントを識別することをさらに含む、項目４８から５７のいずれか一項に記載の方法。
（項目５９）
第２のフィルタリングされたバリアントのインデックスを生成する前記工程が、前記複数のバリアントにおける第２の生殖系列バリアントを識別することをさらに含み、前記第２の生殖系列バリアントが、バリアントの前記第２の参照セット中に存在するデータベースバリアントである、項目４８から５８のいずれか一項に記載の方法。
（項目６０）
前記近似範囲が、第２の生殖系列バリアントの前記対立遺伝子頻度から０．０５の最大値および最小値を有する範囲である、項目４８から５９のいずれか一項に記載の方法。
（項目６１）
前記近似範囲が、第２の生殖系列バリアントの対立遺伝子頻度を中心として、第２の生殖系列バリアントの前記対立遺伝子頻度の二項分布から２標準偏差の最大値および最小値を有する範囲である、項目４８から６０のいずれか一項に記載の方法。
（項目６２）
前記第２の生殖系列バリアントが、前記第２の生殖系列バリアントと同じビン中の少なくとも５つのデータベースバリアントの対立遺伝子頻度に近接した閾値内の対立遺伝子頻度を有する、項目４８から６１のいずれか一項に記載の方法。
（項目６３）
前記第２の生殖系列バリアントが、前記第２の生殖系列バリアントと同じビン中の少なくとも１０個のデータベースバリアントの対立遺伝子頻度に近接した閾値内の対立遺伝子頻度を有する、項目４８から６２のいずれか一項に記載の方法。
（項目６４）
腫瘍細胞を含む前記生体試料が、血清試料、大便試料、血液試料、腫瘍試料から選択される、項目４８から６２のいずれか一項に記載の方法。
（項目６５）
前記腫瘍試料が固定される、項目６４に記載の方法。
本発明は、例えば、以下の項目を提供する。
（項目１）
複数のバリアントにおける体細胞バリアントを識別するための方法であって、
（ａ）体細胞バリアントおよび生殖系列バリアントを含む複数のバリアントを得るステップ；
（ｂ）前記複数のバリアントにデータベースフィルターを適用するステップであって、
前記複数のバリアントにおける第１の生殖系列バリアントを決定する工程を含み、前記第１の生殖系列バリアントが各々、閾値対立遺伝子カウントよりも大きいまたはそれと等しい、バリアントの第１の参照セットにおける対立遺伝子カウントを有する、ステップ；
（ｃ）前記複数のバリアントに近接フィルターを適用するステップであって、
（ｉ）前記複数のバリアントのバリアントを複数のビン中にビニングする工程であって、ゲノムの同じ領域中に位置するバリアントが、同じビン中にビニングされる、工程、
（ｉｉ）前記複数のバリアントにおけるデータベースバリアントを決定する工程であって、データベースバリアントが、バリアントの第２の参照セット中に存在する、工程、および
（ｉｉｉ）前記複数のバリアントにおける第２の生殖系列バリアントを決定する工程であって、前記第２の生殖系列バリアントが各々、前記第２の生殖系列バリアントと同じビン中の少なくとも１つのデータベースバリアントの対立遺伝子頻度の近似範囲内の対立遺伝子頻度を有する、工程
を含む、ステップ；ならびに
（ｄ）前記複数のバリアントから識別された第１および第２の生殖系列バリアントを除去することによって、前記複数のバリアントにおける体細胞バリアントを決定するステップ
を含む、方法。
（項目２）
（ｂ）および（ｃ）が連続的に実施される、項目１に記載の方法。
（項目３）
（ｃ）が（ｂ）の前に実施される、項目１に記載の方法。
（項目４）
前記閾値対立遺伝子カウントが５である、項目１から３のいずれか一項に記載の方法。
（項目５）
前記閾値対立遺伝子カウントが１０である、項目４に記載の方法。
（項目６）
バリアントの前記第１および第２の参照セットが、同じ参照セットである、項目１から５のいずれか一項に記載の方法。
（項目７）
バリアントの前記第１または第２の参照セットが、複数の個体についてのバリアントのデータベースを含む、項目１から６のいずれか一項に記載の方法。
（項目８）
バリアントの前記第１または第２の参照セットが、ゲノム集約データベース（ｇｎｏｍＡＤ）および１０００ゲノムデータベースから選択される少なくとも１つのデータベースを含む、項目１から７のいずれか一項に記載の方法。
（項目９）
ゲノムの前記同じ領域が、同じ染色体内にある、項目１から８のいずれか一項に記載の方法。
（項目１０）
ゲノムの前記同じ領域が、同じ染色体腕内にある、項目１から９のいずれか一項に記載の方法。
（項目１１）
ゲノムの前記同じ領域が、同じ染色体サイトバンド内にある、項目１から１０のいずれか一項に記載の方法。
（項目１２）
ゲノムの前記同じ領域が、１０Ｍｂの領域内にある、項目１から１１のいずれか一項に記載の方法。
（項目１３）
近接フィルターを適用する前記ステップが、０．９よりも大きいまたはそれと等しい対立遺伝子頻度を有する第２の生殖系列バリアントを識別する工程をさらに含む、項目１から１２のいずれか一項に記載の方法。
（項目１４）
近接フィルターを適用する前記ステップが、前記複数のバリアントにおける第２の生殖系列バリアントを識別する工程をさらに含み、前記第２の生殖系列バリアントが、バリアントの前記第２の参照セット中に存在するデータベースバリアントである、項目１から１３のいずれか一項に記載の方法。
（項目１５）
前記近似範囲が、第２の生殖系列バリアントの前記対立遺伝子頻度から０．０５の最大値および最小値を有する範囲である、項目１から１４のいずれか一項に記載の方法。
（項目１６）
前記近似範囲が、第２の生殖系列バリアントの対立遺伝子頻度を中心として、第２の生殖系列バリアントの前記対立遺伝子頻度の二項分布から２標準偏差の最大値および最小値を有する範囲である、項目１から１５のいずれか一項に記載の方法。
（項目１７）
前記第２の生殖系列バリアントが、前記第２の生殖系列バリアントと同じビン中の少なくとも５つのデータベースバリアントの対立遺伝子頻度に近接した閾値内の対立遺伝子頻度を有する、項目１から１６のいずれか一項に記載の方法。
（項目１８）
前記第２の生殖系列バリアントが、前記第２の生殖系列バリアントと同じビン中の少なくとも１０個のデータベースバリアントの対立遺伝子頻度に近接した閾値内の対立遺伝子頻度を有する、項目１から１７のいずれか一項に記載の方法。
（項目１９）
（ａ）が、腫瘍細胞を含む生体試料から配列データを得る工程を含む、項目１から１８のいずれか一項に記載の方法。
（項目２０）
前記配列データを参照配列とアラインさせる工程、および前記配列データ中のバリアントを識別する工程をさらに含む、項目１９に記載の方法。
（項目２１）
腫瘍細胞を含む前記生体試料が、血清試料、大便試料、血液試料、腫瘍試料から選択される、項目１９または２０に記載の方法。
（項目２２）
前記腫瘍試料が固定される、項目２１に記載の方法。
（項目２３）
腫瘍の腫瘍突然変異量を決定する方法であって、
腫瘍細胞を含む生体試料から配列データを得るステップ；
前記配列データから複数のバリアントを決定するステップ；および
項目１から２２のいずれか一項に記載の方法に従って、複数のバリアントにおける体細胞バリアントの数を決定するステップであって、体細胞バリアントの前記数が、前記腫瘍の前記腫瘍突然変異量である、ステップ
を含む、方法。
（項目２４）
腫瘍を処置する方法であって、
項目２３に記載の方法に従って、１０体細胞バリアントよりも大きいまたはそれと等しい腫瘍突然変異量を有する腫瘍を決定するステップ、
有効量のチェックポイント阻害剤を投与することによって、前記腫瘍を処置するステップ
を含む、方法。
（項目２５）
前記腫瘍が、結腸直腸腫瘍、肺腫瘍、子宮内膜腫瘍、子宮腫瘍、胃腫瘍、黒色腫、乳房腫瘍、膵腫瘍、腎臓腫瘍、膀胱腫瘍および脳腫瘍からなる群から選択される、項目２４に記載の方法。
（項目２６）
前記チェックポイント阻害剤が、ＣＴＬＡ－４阻害剤、ＰＤ－１阻害剤およびＰＤ－Ｌ１阻害剤からなる群から選択される、項目２４または２５に記載の方法。
（項目２７）
前記チェックポイント阻害剤が、イピリムマブ、ニボルマブ、ペムブロリズマブ、スパルタリズマブ、アテゾリズマブ、アベルマブおよびデュルバルマブからなる群から選択される、項目２４から２６のいずれか一項に記載の方法。
（項目２８）
遺伝的変異データを分析するための電子システムであって、
プロセッサー上で実行される、腫瘍細胞を含む生体試料由来の配列データから複数のバリアントを識別するように適合されたインフォマティクスモジュールであって、前記複数のバリアントが、体細胞バリアントおよび生殖系列バリアントを含む、インフォマティクスモジュール；
前記複数のバリアントから第１の生殖系列バリアントを除去するように適合されたデータベースフィルターモジュールであって、前記第１の生殖系列バリアントが各々、閾値対立遺伝子カウントよりも大きいまたはそれと等しい、バリアントの第１の参照セットにおける対立遺伝子カウントを有する、データベースフィルターモジュール；
前記複数のバリアントから第２の生殖系列バリアントを除去するように適合された近接フィルターモジュールであって、前記近接フィルターモジュールが、
複数のビンを返すように適合されたビニングサブモジュールであって、各ビンが、ゲノムの同じ領域中に位置する前記複数のバリアントのバリアントを含む、ビニングサブモジュール、
前記複数のバリアントにおけるデータベースバリアントを返すように適合された識別サブモジュールであって、データベースバリアントが、バリアントの第２の参照セット中に存在する、識別サブモジュール、および
前記複数のバリアントから第２の生殖系列バリアントを除去するように適合された除去サブモジュールであって、前記第２の生殖系列バリアントが各々、前記第２の生殖系列バリアントと同じビン中の少なくとも１つのデータベースバリアントの対立遺伝子頻度の近似範囲内の対立遺伝子頻度を有する、除去サブモジュール
を含む、近接フィルターモジュール；ならびに
前記複数のバリアントから除去されなかったバリアントを返すように適合されたディスプレイモジュール
を含む、電子システム。
（項目２９）
インフォマティクスモジュールが、バリアントアノテーションツールを含む、項目２８に記載のシステム。
（項目３０）
前記閾値対立遺伝子カウントが５である、項目２８または２９に記載のシステム。
（項目３１）
前記閾値対立遺伝子カウントが１０である、項目３０に記載のシステム。
（項目３２）
バリアントの前記第１および第２の参照セットが、同じ参照セットである、項目２８から３１のいずれか一項に記載のシステム。
（項目３３）
バリアントの前記第１または第２の参照セットが、複数の個体についてのバリアントのデータベースを含む、項目２８から３２のいずれか一項に記載のシステム。
（項目３４）
バリアントの前記第１または第２の参照セットが、ゲノム集約データベース（ｇｎｏｍＡＤ）および１０００ゲノムデータベースから選択される少なくとも１つのデータベースを含む、項目２８から３３のいずれか一項に記載のシステム。
（項目３５）
ゲノムの前記同じ領域が、同じ染色体内にある、項目２８から３４のいずれか一項に記載のシステム。
（項目３６）
ゲノムの前記同じ領域が、同じ染色体腕内にある、項目２８から３５のいずれか一項に記載のシステム。
（項目３７）
ゲノムの前記同じ領域が、同じ染色体サイトバンド内にある、項目２８から３６のいずれか一項に記載のシステム。
（項目３８）
ゲノムの前記同じ領域が、１０Ｍｂの領域内にある、項目２８から３７のいずれか一項に記載のシステム。
（項目３９）
前記除去サブモジュールが、前記複数のバリアントから、０．９よりも大きいまたはそれと等しい対立遺伝子頻度を有するバリアントを除去するように適合される、項目２８から３８のいずれか一項に記載のシステム。
（項目４０）
前記除去サブモジュールが、前記複数のバリアントから、バリアントの前記第２の参照セット中に存在するデータベースバリアントを除去するように適合される、項目２８から３９のいずれか一項に記載のシステム。
（項目４１）
前記近似範囲が、第２の生殖系列バリアントの前記対立遺伝子頻度から０．０５の最大値および最小値を有する範囲である、項目２８から４０のいずれか一項に記載のシステム。
（項目４２）
前記近似範囲が、第２の生殖系列バリアントの対立遺伝子頻度を中心として、第２の生殖系列バリアントの前記対立遺伝子頻度の二項分布から２標準偏差の最大値および最小値を有する範囲である、項目２８から４１のいずれか一項に記載のシステム。
（項目４３）
前記第２の生殖系列バリアントが、前記第２の生殖系列バリアントと同じビン中の少なくとも５つのデータベースバリアントの対立遺伝子頻度に近接した閾値内の対立遺伝子頻度を有する、項目２８から４２のいずれか一項に記載のシステム。
（項目４４）
前記第２の生殖系列バリアントが、前記第２の生殖系列バリアントと同じビン中の少なくとも１０個のデータベースバリアントの対立遺伝子頻度に近接した閾値内の対立遺伝子頻度を有する、項目２８から４３のいずれか一項に記載のシステム。
（項目４５）
腫瘍細胞を含む前記生体試料が、血清試料、大便試料、血液試料、腫瘍試料から選択される、項目２８から４４のいずれか一項に記載のシステム。
（項目４６）
前記腫瘍試料が固定される、項目４５に記載のシステム。
（項目４７）
複数のバリアントにおける体細胞バリアントを識別するためのコンピューター実装方法であって、項目１から２２のいずれか一項に記載の方法を実施するステップを含む、方法。
（項目４８）
複数のバリアントにおける体細胞バリアントを識別するためのコンピューター実装方法であって、
（ａ）腫瘍細胞を含む生体試料由来の配列データから複数のバリアントを受け取るステップであって、前記複数のバリアントが、体細胞バリアントおよび生殖系列バリアントを含む、ステップ；
（ｂ）前記複数のバリアントにデータベースフィルターを適用するステップであって、
前記複数のバリアントについてドキュメントのインデックスを作成する工程、
前記インデックスを用いてバリアントの第１の参照セットを検索して、前記インデックスにおける第１の生殖系列バリアントを識別する工程であって、前記第１の生殖系列バリアントが各々、閾値対立遺伝子カウントよりも大きいまたはそれと等しい、バリアントの前記第１の参照セットにおける対立遺伝子カウントを有する、工程、および
前記インデックスから識別された第１の生殖系列バリアントを除去して、第１のフィルタリングされたバリアントのインデックスを作成する工程
を含む、ステップ；
（ｃ）第１のフィルタリングされたバリアントの前記インデックスに近接フィルターを適用するステップであって、
（ｉ）ゲノムの異なる領域について複数のビンを作成する工程、
（ｉｉ）第１のフィルタリングされたバリアントの前記インデックスのバリアントをビニングする工程であって、ゲノムの同じ領域中に位置するバリアントが、同じビン中にビニングされる、工程、
（ｉｉｉ）第１のフィルタリングされたバリアントの前記インデックスを用いてバリアントの第２の参照セットを検索して、第１のフィルタリングされたバリアントの前記インデックスにおけるデータベースバリアントを識別する工程、
（ｉｖ）第２の生殖系列バリアントを識別することによって、第１のフィルタリングされたバリアントの前記インデックスから第２の生殖系列バリアントのインデックスを生成する工程であって、前記第２の生殖系列バリアントが各々、前記第２の生殖系列バリアントと同じビン中の少なくとも１つのデータベースバリアントの対立遺伝子頻度の近似範囲内の対立遺伝子頻度を有する、工程、および
（ｖ）第１のフィルタリングされたバリアントの前記インデックスから識別された第２の生殖系列バリアントを除去して、体細胞バリアントのインデックスを作成し、それによって、前記複数のバリアントにおける体細胞バリアントを識別する工程
を含む、ステップ
を含む、方法。
（項目４９）
前記閾値対立遺伝子カウントが５である、項目４８に記載の方法。
（項目５０）
前記閾値対立遺伝子カウントが１０である、項目４９に記載の方法。
（項目５１）
バリアントの前記第１および第２の参照セットが、同じ参照セットである、項目４８から５０のいずれか一項に記載の方法。
（項目５２）
バリアントの前記第１または第２の参照セットが、複数の個体についてのバリアントのデータベースを含む、項目４８から５１のいずれか一項に記載の方法。
（項目５３）
バリアントの前記第１または第２の参照セットが、ゲノム集約データベース（ｇｎｏｍＡＤ）および１０００ゲノムデータベースから選択される少なくとも１つのデータベースを含む、項目４８から５２のいずれか一項に記載の方法。
（項目５４）
ゲノムの前記同じ領域が、同じ染色体内にある、項目４８から５３のいずれか一項に記載の方法。
（項目５５）
ゲノムの前記同じ領域が、同じ染色体腕内にある、項目４８から５４のいずれか一項に記載の方法。
（項目５６）
ゲノムの前記同じ領域が、同じ染色体サイトバンド内にある、項目４８から５５のいずれか一項に記載の方法。
（項目５７）
ゲノムの前記同じ領域が、１０Ｍｂの領域内にある、項目４８から５６のいずれか一項に記載の方法。
（項目５８）
第２のフィルタリングされたバリアントのインデックスを生成する前記工程が、０．９よりも大きいまたはそれと等しい対立遺伝子頻度を有する第２の生殖系列バリアントを識別することをさらに含む、項目４８から５７のいずれか一項に記載の方法。
（項目５９）
第２のフィルタリングされたバリアントのインデックスを生成する前記工程が、前記複数のバリアントにおける第２の生殖系列バリアントを識別することをさらに含み、前記第２の生殖系列バリアントが、バリアントの前記第２の参照セット中に存在するデータベースバリアントである、項目４８から５８のいずれか一項に記載の方法。
（項目６０）
前記近似範囲が、第２の生殖系列バリアントの前記対立遺伝子頻度から０．０５の最大値および最小値を有する範囲である、項目４８から５９のいずれか一項に記載の方法。
（項目６１）
前記近似範囲が、第２の生殖系列バリアントの対立遺伝子頻度を中心として、第２の生殖系列バリアントの前記対立遺伝子頻度の二項分布から２標準偏差の最大値および最小値を有する範囲である、項目４８から６０のいずれか一項に記載の方法。
（項目６２）
前記第２の生殖系列バリアントが、前記第２の生殖系列バリアントと同じビン中の少なくとも５つのデータベースバリアントの対立遺伝子頻度に近接した閾値内の対立遺伝子頻度を有する、項目４８から６１のいずれか一項に記載の方法。
（項目６３）
前記第２の生殖系列バリアントが、前記第２の生殖系列バリアントと同じビン中の少なくとも１０個のデータベースバリアントの対立遺伝子頻度に近接した閾値内の対立遺伝子頻度を有する、項目４８から６２のいずれか一項に記載の方法。
（項目６４）
腫瘍細胞を含む前記生体試料が、血清試料、大便試料、血液試料、腫瘍試料から選択される、項目４８から６２のいずれか一項に記載の方法。
（項目６５）
前記腫瘍試料が固定される、項目６４に記載の方法。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
