JP7277438B2 - ゲノムデータ解析に近縁性を活用するシステム及び方法 - Google Patents

ゲノムデータ解析に近縁性を活用するシステム及び方法 Download PDF

Info

Publication number
JP7277438B2
JP7277438B2 JP2020513709A JP2020513709A JP7277438B2 JP 7277438 B2 JP7277438 B2 JP 7277438B2 JP 2020513709 A JP2020513709 A JP 2020513709A JP 2020513709 A JP2020513709 A JP 2020513709A JP 7277438 B2 JP7277438 B2 JP 7277438B2
Authority
JP
Japan
Prior art keywords
variants
sample
less
putative
novo
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2020513709A
Other languages
English (en)
Other versions
JP2020532989A (ja
JP2020532989A5 (ja
Inventor
ステープルズ、ジェフリー
ハベガー、ルーカス
リード、ジェフリー
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Regeneron Pharmaceuticals Inc
Original Assignee
Regeneron Pharmaceuticals Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Regeneron Pharmaceuticals Inc filed Critical Regeneron Pharmaceuticals Inc
Publication of JP2020532989A publication Critical patent/JP2020532989A/ja
Publication of JP2020532989A5 publication Critical patent/JP2020532989A5/ja
Priority to JP2023076708A priority Critical patent/JP2023103325A/ja
Application granted granted Critical
Publication of JP7277438B2 publication Critical patent/JP7277438B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B20/00ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B20/00ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
    • G16B20/20Allele or variant detection, e.g. single nucleotide polymorphism [SNP] detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/28Databases characterised by their database models, e.g. relational or object models
    • G06F16/284Relational databases
    • G06F16/288Entity relationship models
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B30/00ICT specially adapted for sequence analysis involving nucleotides or amino acids
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B40/00ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B5/00ICT specially adapted for modelling or simulations in systems biology, e.g. gene-regulatory networks, protein interaction networks or metabolic networks
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H10/00ICT specially adapted for the handling or processing of patient-related medical or healthcare data
    • G16H10/60ICT specially adapted for the handling or processing of patient-related medical or healthcare data for patient-specific data, e.g. for electronic patient records
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H70/00ICT specially adapted for the handling or processing of medical references
    • G16H70/60ICT specially adapted for the handling or processing of medical references relating to pathologies

Description

関連出願の相互参照
本出願は、2017年9月7日提出の米国仮特許出願第62/555,597号の利益を主張するものであり、この米国仮特許出願の内容は、その全体が、参照することにより本明細書に援用される。さらに、2018年9月7日提出の「System and Method for Predicting Relatedness in a Human Population」と題する同時係属出願もまた、その全体が、参照することにより援用される。
分野
本開示は、一般に、ゲノムデータの分析のための方法及びシステム、ならびに大集団コホートにおける近縁性を利用して、稀少遺伝的バリエーションを疾患及び疾患感受性に結び付けることに関する。より詳細には、本開示は、同祖性を確立し、遺伝的バリアントを複合ヘテロ接合性変異またはde novo変異としてフェージングするためのシステム及び方法に関する。
ヒトの病状は、環境要因によってのみならず、遺伝的要因によっても引き起こされ影響される。したがって、ヒト集団の遺伝的バリエーションについて理解することが、ヒト疾患の病因及び進行の理解、ならびにこれらの疾患を治療するための新規薬剤ターゲットの同定にとって重要である。
保健医療集団の遺伝的研究は、この点で、広範な健康管理データが使用できることのために特に有効である。これにより、遺伝的バリアントがヒトの病状にどのように寄与するかの研究が簡略化される。過去には、通常、そのような研究は、疾患遺伝子座をマッピングするための全ゲノム遺伝的連鎖解析に基づいており、これらは同定されると、同定後に、分子レベルでさらに詳細に解析することができた。
ここ数年にわたって、ハイスループットなDNAシークエンシングテクノロジが広く一般に使用できるようになったことで、数十万人のヒトのゲノムの並列シークエンシングが可能になっている。理論上、これらのデータは、ヒト疾患の遺伝的基礎を解読するために使用できる強力な情報源を意味する。しかしながら、これらの増え続けているデータセットは、このような非常に大きなデータセットの効率的な処理を継続するために、バイオインフォマティクスのツール及び解析パイプラインの持続的な革新が必要とされている。さらに、このような大きなデータセットにおける近縁性及び家族構造の有用性、ならびにこれをバリアントの同定及び特性評価に活用できる範囲に関しては、完全には認識されておらず、利用されていない。
大規模なゲノムデータの解析のための、改良されたバイオインフォマティクスのツールが必要とされている。本開示は、このニーズに応えるものである。
一態様では、本開示は、集団における遺伝的バリアントを、集団の近縁性を活用することによってフェージングする方法であって、複数のヒト被験者から得られた核酸配列サンプルのデータセットから、低クオリティの配列バリアントを除去すること、サンプルの1つ以上のそれぞれについて祖先スーパークラスの指定を確立すること、低クオリティサンプルをデータセットから除去すること、祖先スーパークラス内での被験者の第1の同祖性推定を生成すること、被験者の祖先スーパークラスから独立した被験者の第2の同祖性推定を生成すること、第2の同祖性推定の1つ以上に基づいて被験者を第1の第一度家族ネットワークにクラスタリングすること、第1の第一度家族ネットワーク内での被験者の第3の同祖性推定を生成すること、第1の同祖性推定と第3の同祖性推定とを統合して、統合同祖性推定を取得すること、統合同祖性推定に基づいて被験者の第2の第一度家族ネットワークを構築すること、ならびにバリアントを、統合同祖性推定及び第2の第一度家族ネットワークに応じて、複合ヘテロ接合性変異(CHM)であるかもしくはCHMではないものとしてフェージングし、またはバリアントを、統合同祖性推定及び第2の第一度家族ネットワークに応じて、de novo変異(DNM)として同定することを含む方法を提供する。
いくつかの例示的な実施形態では、第1の同祖性推定と第3の同祖性推定とを統合することは、第3の同祖性推定に特有のペアワイズ同祖性推定を補って、第1の同祖性推定を増強することを含む。
いくつかの例示的な実施形態では、バリアントを複合ヘテロ接合性変異(CHM)としてフェージングすることは、(1)集団アレル頻度に応じてバリアントをフェージングすること、(2)ハーディ・ワインベルグ平衡(HWE)から外れたバリアント、もしくは同じサンプル内の別のバリアントの10塩基対内のバリアント、またはその両方を除去すること、及びクオリティバイデプス(quality by depth)(QD)が約2以下、またはリードデプス(DP)が約5未満、または代替アレルバランス(AB)が約10%以下、またはこれらの組み合わせの一塩基多型(SNP)を除去すること、及びQDが約2以下、またはDPが約5未満、またはABが約10%以下、またはこれらの組み合わせの挿入または欠失(INDEL)を除去すること、(3)同じサンプル内及び同じ遺伝子内に1つ以上のバリアントの対が存在する潜在複合ヘテロ接合性変異(pCHM)として残りのバリアントを選択すること、ならびに(4)pCHMをシスpCHMまたはトランスpCHMとしてフェージングし、次いでトランスpCHMとしてフェージングされたpCHMをCHMとして分類することを含む。
いくつかの例示的な実施形態では、バリアントを複合ヘテロ接合性変異としてフェージングすることは、ハーディ・ワインベルグ平衡(HWE)から外れたバリアント、もしくは同じサンプル内の別のバリアントの10塩基対内のバリアント、またはその両方を除去すること、及びクオリティバイデプス(quality by depth)(QD)が約3以下、またはリードデプス(DP)が約7未満、または代替アレルバランス(AB)が約15%以下、またはこれらの組み合わせのSNPを除去すること、及びQDが約5以下、またはDPが約10未満、またはABが約20%以下、またはこれらの組み合わせの挿入または欠失(INDEL)を除去することを含む。
いくつかの例示的な実施形態では、本方法は、さらに、(1)機能的効果優先度に応じてCHMにスコアを付与すること、及び(2)サンプル当たりの遺伝子ごとに最も高い機能的効果優先度スコアを有するCHMを選択し、ヒトが同一遺伝子内に複数のCHMを有する場合に、タンパク質機能抑制をもたらす可能性が最も高いCHMが同定されるようにすることをさらに含む。
いくつかの例示的な実施形態では、バリアントをde novo変異としてフェージングすることは、(1)第2の第一度家族ネットワーク及びそのトリオのサンプルの中のバリアントを同定すること、(2)遺伝子型尤度スコアを、トリオ内の親サンプル及び対応する子サンプルのバリアントに割り当て、バリアントがde novo変異である確率を計算し、計算した確率が統計的に有意な場合に、そのバリアントを推定de novo変異として同定すること、(3)トリオの子サンプルのバリアントを同定し、トリオのいずれの親サンプルにもバリアントが存在しない場合に、そのバリアントを推定de novo変異として同定すること、(4)子サンプルに約35未満の遺伝子型クオリティ(GQ)アノテーションを有する推定de novo変異、または複数のヒト被験者由来のサンプル全体にわたって10以上の代替アレルカウント(AC)を有する推定de novo変異、または子サンプルに、約7未満のリードデプス(DP)と約4未満の代替DPとを有する推定de novo変異、またはどちらか一方の親サンプルに約2%を超えるアレルバランス(AB)を有する推定de novo変異、または子サンプルに約15%未満のアレルバランス(AB)を有する推定de novo変異、または子サンプルに約90%を超えるABを有する推定de novo変異、またはどちらか一方の親サンプルに代替アレルホモ接合性を有する推定de novo変異、またはこれらの組み合わせの推定de novo変異を除去することによって同定される推定de novo変異をフィルタ処理すること、ならびに(5)フィルタ処理済みの同定された推定de novo変異を組み合わせて、それによって推定de novo変異データセットを形成することを含む。
いくつかの例示的な実施形態では、本方法は、推定de novo変異データセット内の推定de novo変異を、推定de novo変異のアレルバランスが、子サンプルで約0.15以上、かつ各親サンプルで約0.02以下の場合、及び推定de novo変異のマッピングクオリティが約40未満ではない場合、及び推定de novo変異のクオリティバイデプス(quality by depth)(QD)の値が約2未満ではない場合、及び推定de novo変異のMACがサンプル全体にわたって約20未満の場合、及び推定de novo変異のソフトクリップリードが、推定de novo変異のキャリアにおけるバリアントサイトで約3以下の場合、及び推定de novo変異が、約4を超えるモノポリマーランを持つINDELではない場合に、中信頼度de novo変異として分類することをさらに含む。
いくつかの例示的な実施形態では、本方法は、中信頼度de novo変異を、中信頼度de novo変異が、親サンプルにおいて約90以上の遺伝子型クオリティアノテーションを有する場合、及び中信頼度de novo変異が、各親サンプルにおいて約10以上のリードデプスを有する場合、及び中信頼度de novo変異が、子サンプルにおいて約7以上の代替リードデプスを有する場合、及び中信頼度de novo変異が、SNPについて約3より大きいQDを有する場合、及び中信頼度de novo変異が、INDELについて約5より大きいQDを有する場合に、高信頼度de novo変異として分類することをさらに含む。
一態様では、本開示は、集団における複合ヘテロ接合性変異(CHM)を同定する方法であって、複数のヒト被験者由来のDNA配列サンプルの中のバリアントを同定すること、同定されたバリアントに基づいて被験者の祖先スーパークラスの指定を確立すること、祖先スーパークラス内での被験者の第1の同祖性推定を生成すること、被験者の祖先スーパークラスから独立した被験者の第2の同祖性推定を生成すること、第2の同祖性推定の1つ以上に基づいて被験者を第1の第一度家族ネットワークにクラスタリングすること、第1の第一度家族ネットワーク内での被験者の第3の同祖性推定を生成すること、第1の同祖性推定と第3の同祖性推定とを統合して、統合同祖性推定を取得すること、統合同祖性推定に基づいて第2の第一度家族ネットワークを構築すること、集団アレル頻度に応じて、サンプルの中のバリアントをフェージングすること、フェージングされたバリアントを、同じ被験者及び遺伝子における2つ以上のバリアントの存在に基づいて潜在CHMとして分類すること、ならびに潜在CHMを、同じ被験者及び遺伝子における別バリアントとシスであるまたはトランスであるとしてフェージングし、次いでトランスであるとしてフェージングされた潜在CHMをCHMとして分類することを含む方法を提供する。
いくつかの例示的な実施形態では、本方法は、被験者の祖先スーパークラスの指定が確立される前に、同定されたバリアントをフィルタ処理することをさらに含む。
いくつかの例示的な実施形態では、本方法は、被験者の第1の同祖性推定及び第2の同祖性推定が生成される前に、同定されたバリアントをフィルタ処理することをさらに含む。
いくつかの例示的な実施形態では、バリアントをフィルタ処理することは、複数のヒト被験者由来のサンプル全体にわたって約10%を超える代替アレル頻度を有するバリアントを除去すること、またはp値が約10-6超である、ハーディ・ワインベルグ平衡(HWE)に反するバリアントを除去すること、または複数のヒト被験者由来のサンプルの約5%超でコールが欠損するバリアントを除去すること、またはこれらの組み合わせを含む。
いくつかの例示的な実施形態では、本方法は、同定されたバリアントがフィルタ処理された後に、低クオリティサンプルを除去することを含む。
いくつかの例示的な実施形態では、低クオリティサンプルは、D-statが0.12を超えるか、もしくは20xのリードカバレッジが75%未満であるか、またはその両方であるサンプルである。
いくつかの例示的な実施形態では、第1の同祖性推定と第3の同祖性推定とを統合することは、第3の同祖性推定に特有のペアワイズ同祖性推定を補って、第1の同祖性推定を増強することを含む。
いくつかの例示的な実施形態では、同祖性推定は、サンプルペアの間のゲノム全域にわたる0、1、及び2のIBD値の計算を含む。
いくつかの例示的な実施形態では、本方法は、バリアントが集団アレル頻度に応じてフェージングされた後に、バリアントをフィルタ処理することをさらに含む。
いくつかの例示的な実施形態では、集団アレル頻度に応じてフェージングされたバリアントをフィルタ処理することは、ハーディ・ワインベルグ平衡(HWE)から外れたバリアント、もしくは同じサンプル内の別のバリアントの10塩基対内のバリアント、またはその両方を除去すること、及びクオリティバイデプス(quality by depth)(QD)が約2以下、またはリードデプス(DP)が約5未満、または代替アレルバランス(AB)が約10%以下、またはこれらの組み合わせのSNPを除去すること、及びQDが約2以下、またはDPが約5未満、またはABが約10%以下、またはこれらの組み合わせの挿入または欠失(INDEL)を除去することを含む。
いくつかの例示的な実施形態では、集団アレル頻度に応じてバリアントをフェージングすることは、ヒト被験者のDNA配列サンプルを、遺伝子間領域において、近似的に等しいサイズであり、実質的なセグメントオーバーラップ及び切断点を有するゲノムセグメントに分割することを含む。
いくつかの例示的な実施形態では、潜在CHMは、トリオのデータ、もしくは親子のデータ、もしくは完全同胞のデータ、もしくは遠縁の血縁者のデータ、もしくはそれらの組み合わせに基づいてフェージングされる、またはマイナーアレルカウント(MAC)に基づいてフェージングされる、または集団アレル頻度に基づいてフェージングされる、またはそれらの組み合わせである。
いくつかの例示的な実施形態では、本方法は、機能的効果優先度に応じてCHMにスコアを付与すること、及びサンプル当たりの遺伝子ごとに最も高い機能的効果優先度スコアを有するCHMを選択し、それによって医学的に重要な変異の集合を取得することをさらに含む。
いくつかの例示的な実施形態では、DNA配列サンプルはエクソーム配列を含む。
いくつかの例示的な実施形態では、複数のヒト被験者は、10Kを超える被験者を含む。
いくつかの例示的な実施形態では、低クオリティサンプルのD-statは、Kolmogorov-Smirnov(KS)検定を使用して、サンプルの実際のアレルバランスの分布を、予想されるアレルバランスの分布と比較することによって判定される。
いくつかの例示的な実施形態では、集団アレル頻度に応じてフェージングされたバリアントをフィルタ処理することは、ハーディ・ワインベルグ平衡(HWE)から外れたバリアント、もしくは同じサンプル内の別のバリアントの10塩基対内のバリアント、またはその両方を除去すること、及びクオリティバイデプス(quality by depth)(QD)が約3以下、またはリードデプス(DP)が約7未満、または代替アレルバランス(AB)が約15%以下、またはこれらの組み合わせのSNPを除去すること、及びQDが約5以下、またはDPが約10未満、またはABが約20%以下、またはこれらの組み合わせの挿入または欠失(INDEL)を除去することを含む。
別の態様では、本開示は、集団における複合ヘテロ接合性変異(CHM)を同定する非一時的なコンピュータ実行方法を提供する。概して、本非一時的なコンピュータ実行方法は、複数のヒト被験者由来のDNA配列サンプルの中のバリアントを同定するために、コンピューティングデバイスのデータプロセッサを使用すること、同定されたバリアントに基づいて被験者の祖先スーパークラスの指定を確立するために、データプロセッサを使用すること、祖先スーパークラス内での被験者の第1の同祖性推定を生成するために、データプロセッサを使用すること、被験者の祖先スーパークラスから独立した被験者の第2の同祖性推定を生成するために、データプロセッサを使用すること、第2の同祖性推定の1つ以上に基づいて被験者を第1の第一度家族ネットワークにクラスタリングするために、データプロセッサを使用すること、第1の第一度家族ネットワーク内での被験者の第3の同祖性推定を生成するために、データプロセッサを使用すること、第1の同祖性推定と第3の同祖性推定とを統合して、統合同祖性推定を取得するために、データプロセッサを使用すること、統合同祖性推定に基づいて第2の第一度家族ネットワークを構築するために、データプロセッサを使用すること、集団アレル頻度に応じて、サンプルの中のバリアントをフェージングするために、データプロセッサを使用すること、フェージングされたバリアントを、同じ被験者及び遺伝子における2つ以上のバリアントの存在に基づいて潜在CHMとして分類するために、データプロセッサを使用すること、ならびに潜在CHMを、同じ被験者及び遺伝子における別バリアントとシスであるまたはトランスであるとしてフェージングし、次いでトランスであるとしてフェージングされた潜在CHMをCHMとして分類するために、データプロセッサを使用することを含む。
いくつかの例示的な実施形態では、本非一時的なコンピュータ実行方法は、被験者の祖先スーパークラスの指定が確立される前に、同定されたバリアントをフィルタ処理するために、データプロセッサを使用することをさらに含む。
いくつかの例示的な実施形態では、本非一時的なコンピュータ実行方法は、被験者の第2の同祖性推定が生成される前に、同定されたバリアントをフィルタ処理するために、データプロセッサを使用することをさらに含む。
いくつかの例示的な実施形態では、バリアントをフィルタ処理することは、複数のヒト被験者由来のサンプル全体にわたって約10%を超える代替アレル頻度を有するバリアントを除去すること、またはp値が約10-6超である、ハーディ・ワインベルグ平衡(HWE)に反するバリアントを除去すること、または複数のヒト被験者由来のサンプルの約5%超でコールが欠損するバリアントを除去すること、またはこれらの組み合わせを含む。
いくつかの例示的な実施形態では、本非一時的なコンピュータ実行方法は、同定されたバリアントがフィルタ処理された後に、低クオリティサンプルを除去するために、データプロセッサを使用することをさらに含む。
いくつかの例示的な実施形態では、低クオリティサンプルは、D-statが0.12を超えるか、もしくは20xのリードカバレッジが75%未満であるか、またはその両方であるサンプルである。
いくつかの例示的な実施形態では、第1の同祖性推定と第3の同祖性推定とを統合することは、第3の同祖性推定に特有のペアワイズ同祖性推定を補って、第1の同祖性推定を増強することを含む。
いくつかの例示的な実施形態では、同祖性推定は、サンプルペアの間のゲノム全域にわたる0、1、及び2のIBD値の計算を含む。
いくつかの例示的な実施形態では、本非一時的なコンピュータ実行方法は、バリアントが集団アレル頻度に応じてフェージングされた後に、バリアントをフィルタ処理するために、データプロセッサを使用することをさらに含む。
いくつかの例示的な実施形態では、集団アレル頻度に応じてフェージングされたバリアントをフィルタ処理することは、ハーディ・ワインベルグ平衡から外れたバリアント、もしくは同じサンプル内の別のバリアントの10塩基対内のバリアント、またはその両方を除去すること、及びクオリティバイデプス(quality by depth)(QD)が約2以下、またはリードデプス(DP)が約5未満、または代替アレルバランス(AB)が約10%以下、またはこれらの組み合わせのSNPを除去すること、及びQDが約2以下、またはDPが約5未満、またはABが約10%以下、またはこれらの組み合わせの挿入または欠失(INDEL)を除去することを含む。
いくつかの例示的な実施形態では、集団アレル頻度に応じてバリアントをフェージングすることは、ヒト被験者のDNA配列サンプルを、遺伝子間領域において、近似的に等しいサイズであり、実質的なセグメントオーバーラップ及び切断点を有するゲノムセグメントに分割することを含む。
いくつかの例示的な実施形態では、潜在CHMは、トリオのデータ、もしくは親子のデータ、もしくは完全同胞のデータ、もしくは遠縁の血縁者のデータ、もしくはそれらの組み合わせに基づいてフェージングされる、またはマイナーアレルカウント(MAC)に基づいてフェージングされる、または集団アレル頻度に基づいてフェージングされる、またはそれらの組み合わせである。
いくつかの例示的な実施形態では、本非一時的なコンピュータ実行方法は、機能的効果優先度に応じてCHMにスコアを付与すること、及びサンプル当たりの遺伝子ごとに最も高い機能的効果優先度スコアを有するCHMを選択し、それによって医学的に重要な変異の集合を取得することのために、データプロセッサを使用することをさらに含む。
いくつかの例示的な実施形態では、DNA配列サンプルはエクソーム配列を含む。
いくつかの例示的な実施形態では、複数のヒト被験者は、10Kを超える被験者を含む。
いくつかの例示的な実施形態では、低クオリティサンプルのD-statは、KS検定を使用して、サンプルの実際のアレルバランスの分布を、予想されるアレルバランスの分布と比較することによって判定される。
いくつかの例示的な実施形態では、集団アレル頻度に応じてフェージングされたバリアントをフィルタ処理することは、ハーディ・ワインベルグ平衡(HWE)から外れたバリアント、もしくは同じサンプル内の別のバリアントの10塩基対内のバリアント、またはその両方を除去すること、及びクオリティバイデプス(quality by depth)(QD)が約3以下、またはリードデプス(DP)が約7未満、または代替アレルバランス(AB)が約15%以下、またはこれらの組み合わせのSNPを除去すること、及びQDが約5以下、またはDPが約10未満、またはABが約20%以下、またはこれらの組み合わせの挿入または欠失(INDEL)を除去することを含む。
別の態様では、本開示は、本方法及び非一時的なコンピュータ実行方法を実行するためのシステムを提供する。本システムは、概して、データプロセッサ、データプロセッサに結合されたメモリ、及びメモリに格納されたプログラムを含み、プログラムが、複数のヒト被験者由来のDNA配列サンプルの中のバリアントを同定する、同定されたバリアントに基づいて被験者の祖先スーパークラスの指定を確立する、祖先スーパークラス内での被験者の第1の同祖性推定を生成する、被験者の祖先スーパークラスから独立した被験者の第2の同祖性推定を生成する、第2の同祖性推定の1つ以上に基づいて被験者を第1の第一度家族ネットワークにクラスタリングする、第1の第一度家族ネットワーク内での被験者の第3の同祖性推定を生成する、第1の同祖性推定と第3の同祖性推定とを統合して、統合同祖性推定を取得する、統合同祖性推定に基づいて第2の第一度家族ネットワークを構築する、集団アレル頻度に応じて、サンプルの中のバリアントをフェージングする、フェージングされたバリアントを、同じ被験者及び遺伝子における2つ以上のバリアントの存在に基づいて潜在CHMとして分類する、ならびに潜在CHMを、同じ被験者及び遺伝子における別バリアントとシスであるまたはトランスであるとしてフェージングし、次いでトランスであるとしてフェージングされた潜在CHMをCHMとして分類するための命令を含む。
いくつかの例示的な実施形態では、プログラムは、被験者の祖先スーパークラスの指定が確立される前に、同定されたバリアントをフィルタ処理するための命令を含む。
いくつかの例示的な実施形態では、プログラムは、被験者の第1の同祖性推定及び第2の同祖性推定が生成される前に、同定されたバリアントをフィルタ処理するための命令を含む。
いくつかの例示的な実施形態では、バリアントをフィルタ処理することは、複数のヒト被験者由来のサンプル全体にわたって約10%を超える代替アレル頻度を有するバリアントを除去すること、またはp値が約10-6超である、ハーディ・ワインベルグ平衡(HWE)に反するバリアントを除去すること、または複数のヒト被験者由来のサンプルの約5%超でコールが欠損するバリアントを除去すること、またはこれらの組み合わせを含む。
いくつかの例示的な実施形態では、プログラムは、同定されたバリアントがフィルタ処理された後に、低クオリティサンプルを除去するための命令を含む。
いくつかの例示的な実施形態では、低クオリティサンプルは、D-statが0.12を超えるか、もしくは20xのリードカバレッジが75%未満であるか、またはその両方であるサンプルである。
いくつかの例示的な実施形態では、第1の同祖性推定と第3の同祖性推定とを統合することは、第3の同祖性推定に特有のペアワイズ同祖性推定を補って、第1の同祖性推定を増強することを含む。
いくつかの例示的な実施形態では、同祖性推定は、サンプルペアの間のゲノム全域にわたる0、1、及び2のIBD値の計算を含む。
いくつかの例示的な実施形態では、プログラムは、バリアントが集団アレル頻度に応じてフェージングされた後に、バリアントをフィルタ処理するための命令を含む。
いくつかの例示的な実施形態では、集団アレル頻度に応じてフェージングされたバリアントをフィルタ処理することは、ハーディ・ワインベルグ平衡(HWE)から外れたバリアント、もしくは同じサンプル内の別のバリアントの10塩基対内のバリアント、またはその両方を除去すること、及びクオリティバイデプス(quality by depth)(QD)が約2以下、またはリードデプス(DP)が約5未満、または代替アレルバランス(AB)が約10%以下、またはこれらの組み合わせのSNPを除去すること、及びQDが約2以下、またはDPが約5未満、またはABが約10%以下、またはこれらの組み合わせの挿入または欠失(INDEL)を除去することを含む。
いくつかの例示的な実施形態では、集団アレル頻度に応じてバリアントをフェージングすることは、ヒト被験者のDNA配列サンプルを、遺伝子間領域において、近似的に等しいサイズであり、実質的なセグメントオーバーラップ及び切断点を有するゲノムセグメントに分割することを含む。
いくつかの例示的な実施形態では、潜在CHMは、トリオのデータ、もしくは親子のデータ、もしくは完全同胞のデータ、もしくは遠縁の血縁者のデータ、もしくはそれらの組み合わせに基づいてフェージングされる、またはマイナーアレルカウント(MAC)に基づいてフェージングされる、または集団アレル頻度に基づいてフェージングされる、またはそれらの組み合わせである。
いくつかの例示的な実施形態では、プログラムは、機能的効果優先度に応じてCHMにスコアを付与する、及びサンプル当たりの遺伝子ごとに最も高い機能的効果優先度スコアを有するCHMを選択し、それによって医学的に重要な変異の集合を取得するための命令を含む。
いくつかの例示的な実施形態では、DNA配列サンプルはエクソーム配列を含む。
いくつかの例示的な実施形態では、複数のヒト被験者は、10Kを超える被験者を含む。
いくつかの例示的な実施形態では、低クオリティサンプルのD-statは、KS検定を使用して、サンプルの実際のアレルバランスの分布を、予想されるアレルバランスの分布と比較することによって判定される。
いくつかの例示的な実施形態では、集団アレル頻度に応じてフェージングされたバリアントをフィルタ処理することは、ハーディ・ワインベルグ平衡(HWE)から外れたバリアント、もしくは同じサンプル内の別のバリアントの10塩基対内のバリアント、またはその両方を除去すること、及びクオリティバイデプス(quality by depth)(QD)が約3以下、またはリードデプス(DP)が約7未満、または代替アレルバランス(AB)が約15%以下、またはこれらの組み合わせのSNPを除去すること、及びQDが約5以下、またはDPが約10未満、またはABが約20%以下、またはこれらの組み合わせの挿入または欠失(INDEL)を除去することを含む。
別の態様では、本開示は、集団におけるde novo変異(DNM)を同定する方法を提供する。概して、本方法は、複数のヒト被験者由来のDNA配列サンプルの中のバリアントを同定すること、同定されたバリアントに基づいて被験者の祖先スーパークラスの指定を確立すること、祖先スーパークラス内での被験者の第1の同祖性推定を生成すること、被験者の祖先スーパークラスから独立した被験者の第2の同祖性推定を生成すること、第2の同祖性推定の1つ以上に基づいて被験者を第1の第一度家族ネットワークにクラスタリングすること、第1の第一度家族ネットワーク内での被験者の第3の同祖性推定を生成すること、第1の同祖性推定と第3の同祖性推定とを統合して、統合同祖性推定を取得すること、統合同祖性推定に基づいて核家族を構築すること、核家族におけるバリアントを同定すること、構築された核家族のトリオの両親及び子に由来するサンプルのバリアントに遺伝子型尤度スコアを付与して、バリアントがde novo変異である確率を計算するとともに、独立して、トリオのどちらの親サンプルにも存在しない子サンプルのバリアントを単純に同定して、そのバリアントがde novo変異である確率を計算し、次いで両方の確率を組み合わせ、それによって推定de novo変異のデータセットを形成することを含む。
いくつかの例示的な実施形態では、本方法は、被験者の祖先スーパークラスの指定が確立される前に、同定されたバリアントをフィルタ処理することをさらに含む。
いくつかの例示的な実施形態では、本方法は、被験者の第2の同祖性推定が生成される前に、同定されたバリアントをフィルタ処理することをさらに含む。
いくつかの例示的な実施形態では、バリアントをフィルタ処理することは、複数のヒト被験者由来のサンプル全体にわたって約10%を超える代替アレル頻度を有するバリアントを除去すること、またはp値が約10-6超である、ハーディ・ワインベルグ平衡(HWE)に反するバリアントを除去すること、または複数のヒト被験者由来のサンプルの約5%超でコールが欠損するバリアントを除去すること、またはこれらの組み合わせを含む。
いくつかの例示的な実施形態では、本方法は、同定されたバリアントがフィルタ処理された後に、低クオリティサンプルを除去することを含む。
いくつかの例示的な実施形態では、低クオリティサンプルは、D-statが0.12を超えるか、もしくは20xのリードカバレッジが75%未満であるか、またはその両方であるサンプルである。
いくつかの例示的な実施形態では、第1の同祖性推定と第3の同祖性推定とを統合することは、第3の同祖性推定に特有のペアワイズ同祖性推定を補って、第1の同祖性推定を増強することを含む。
いくつかの例示的な実施形態では、同祖性推定は、サンプルペアの間のゲノム全域にわたる0、1、及び2のIBD値の計算を含む。
いくつかの例示的な実施形態では、遺伝子型尤度スコアは、複数の核家族における複数のヒト被験者由来のDNA配列サンプルに基づいている。
いくつかの例示的な実施形態では、本方法は、遺伝子型尤度スコアに基づいて、バリアントがde novo変異である確率が計算された後に、バリアントをフィルタ処理することをさらに含む。
いくつかの例示的な実施形態では、本方法は、どちらの親サンプルにも存在しない子サンプルのバリアントを単純に同定することに基づいて、バリアントがde novo変異である確率が計算された後に、バリアントをフィルタ処理することをさらに含む。
いくつかの例示的な実施形態では、バリアントをフィルタ処理することは、子サンプルに約35未満の遺伝子型クオリティ(GQ)アノテーションを有するバリアント、またはサンプル中で10以上の代替アレルカウント(AC)を有するバリアント、または子サンプルに、約7未満のリードデプス(DP)と約4未満の代替DPとを有するバリアント、またはどちらか一方の親サンプルに約2%を超えるアレルバランス(AB)を有するバリアント、または子サンプルに約15%未満のアレルバランス(AB)を有するバリアント、または子サンプルに約90%を超えるABを有するバリアント、またはどちらか一方の親サンプルに代替アレルホモ接合性を有するバリアント、またはこれらの組み合わせのバリアントを除去することを含む。
いくつかの例示的な実施形態では、本方法は、バリアントに品質管理メトリクスでアノテーションを付けることをさらに含む。
いくつかの例示的な実施形態では、本方法は、どちらの親サンプルにも存在しない子サンプルのバリアントを単純に同定することに基づいて、推定de novo変異が同定された後に、サンプルBAMファイルデータに基づいてバリアントをフィルタ処理することをさらに含む。
いくつかの例示的な実施形態では、本方法は、推定de novo変異を、推定de novo変異のアレルバランスが、子サンプルで約0.15以上の場合に、中信頼度de novo変異として分類することをさらに含む。
いくつかの例示的な実施形態では、本方法は、推定de novo変異を、推定de novo変異のアレルバランスが、各親サンプルで約0.02以下の場合に、中信頼度de novo変異として分類することをさらに含む。
いくつかの例示的な実施形態では、本方法は、推定de novo変異を、推定de novo変異のマッピングクオリティが約40未満ではない場合に、中信頼度de novo変異として分類することをさらに含む。
いくつかの例示的な実施形態では、本方法は、推定de novo変異を、推定de novo変異のクオリティバイデプス(quality by depth)(QD)の値が約2未満ではない場合に、中信頼度de novo変異として分類することをさらに含む。
いくつかの例示的な実施形態では、本方法は、推定de novo変異を、推定de novo変異のMACがサンプル全体にわたって約20未満の場合に、中信頼度de novo変異として分類することをさらに含む。
いくつかの例示的な実施形態では、本方法は、推定de novo変異を、推定de novo変異のソフトクリップリードが、推定de novo変異のキャリアにおけるバリアントサイトで約3以下の場合に、中信頼度de novo変異として分類することをさらに含む。
いくつかの例示的な実施形態では、本方法は、推定de novo変異を、推定de novo変異が、約4を超えるモノポリマーランを持つINDELではない場合に、中信頼度de novo変異として分類することをさらに含む。
いくつかの例示的な実施形態では、本方法は、推定de novo変異を、推定de novo変異のアレルバランス(AB)が、子サンプルで約0.15以上、かつ各親サンプルで約0.02以下の場合、及び推定de novo変異のマッピングクオリティ(MQ)が約40未満ではない場合、及び推定de novo変異のクオリティバイデプス(quality by depth)(QD)の値が約2未満ではない場合、及び推定de novo変異のマイナーアレルカウント(MAC)がサンプル全体にわたって約20未満の場合、及び推定de novo変異のソフトクリップリードが、推定de novo変異のキャリアにおけるバリアントサイトで約3以下の場合、及び推定de novo変異が、約4を超えるモノポリマーランを持つINDELではない場合に、中信頼度de novo変異として分類することをさらに含む。
いくつかの例示的な実施形態では、本方法は、中信頼度de novo変異を、中信頼度de novo変異が、親サンプルにおいて約90以上の遺伝子型クオリティ(GQ)アノテーションを有する場合、及び中信頼度de novo変異が、各親サンプルにおいて約10以上のリードデプス(DP)を有する場合、及び中信頼度de novo変異が、子サンプルにおいて約7以上の代替DPを有する場合、及び中信頼度de novo変異が、約3より大きいQDのSNPを有する場合、及び中信頼度de novo変異が、約5より大きいQDのINDELを有する場合に、高信頼度de novo変異として分類することをさらに含む。
いくつかの例示的な実施形態では、DNA配列サンプルはエクソーム配列を含む。
いくつかの例示的な実施形態では、複数のヒト被験者は、10Kを超える被験者を含む。
いくつかの例示的な実施形態では、低クオリティサンプルのD-statは、KS検定を使用して、サンプルの実際のアレルバランスの分布を、予想されるアレルバランスの分布と比較することによって判定される。
別の態様では、本開示は、集団におけるde novo変異(DNM)を同定する非一時的なコンピュータ実行方法を提供する。概して、本非一時的なコンピュータ実行方法は、複数のヒト被験者由来のDNA配列サンプルの中のバリアントを同定するために、コンピューティングデバイスのデータプロセッサを使用すること、同定されたバリアントに基づいて被験者の祖先スーパークラスの指定を確立するために、データプロセッサを使用すること、祖先スーパークラス内での被験者の第1の同祖性推定を生成するために、データプロセッサを使用すること、被験者の祖先スーパークラスから独立した被験者の第2の同祖性推定を生成するために、データプロセッサを使用すること、第2の同祖性推定の1つ以上に基づいて被験者を第1の第一度家族ネットワークにクラスタリングするために、データプロセッサを使用すること、第1の第一度家族ネットワーク内での被験者の第3の同祖性推定を生成するために、データプロセッサを使用すること、第1の同祖性推定と第3の同祖性推定とを統合して、統合同祖性推定を取得するために、データプロセッサを使用すること、統合同祖性推定に基づいて核家族を構築するために、データプロセッサを使用すること、核家族におけるバリアントを同定するために、データプロセッサを使用すること、構築された核家族のトリオの両親及び子に由来するサンプルのバリアントに遺伝子型尤度スコアを付与して、バリアントがde novo変異である確率を計算するとともに、独立して、トリオのどちらの親サンプルにも存在しない子サンプルのバリアントを単純に同定して、そのバリアントがde novo変異である確率を計算し、次いで両方の確率を組み合わせ、それによって推定de novo変異のデータセットを形成するために、データプロセッサを使用することを含む。
いくつかの例示的な実施形態では、本非一時的なコンピュータ実行方法は、被験者の祖先スーパークラスの指定が確立される前に、同定されたバリアントをフィルタ処理するために、データプロセッサを使用することをさらに含む。
いくつかの例示的な実施形態では、本非一時的なコンピュータ実行方法は、被験者の第2の同祖性推定が生成される前に、同定されたバリアントをフィルタ処理するために、データプロセッサを使用することをさらに含む。
いくつかの例示的な実施形態では、バリアントをフィルタ処理することは、複数のヒト被験者由来のサンプル全体にわたって約10%を超える代替アレル頻度を有するバリアントを除去すること、またはp値が約10-6超である、ハーディ・ワインベルグ平衡に反するバリアントを除去すること、または複数のヒト被験者由来のサンプルの約5%超でコールが欠損するバリアントを除去すること、またはこれらの組み合わせを含む。
いくつかの例示的な実施形態では、本非一時的なコンピュータ実行方法は、同定されたバリアントがフィルタ処理された後に、低クオリティサンプルを除去するために、データプロセッサを使用することをさらに含む。
いくつかの例示的な実施形態では、低クオリティサンプルは、D-statが0.12を超えるか、もしくは20xのリードカバレッジが75%未満であるか、またはその両方であるサンプルである。
いくつかの例示的な実施形態では、第1の同祖性推定と第3の同祖性推定とを統合することは、第3の同祖性推定に特有のペアワイズ同祖性推定を補って、第1の同祖性推定を増強することを含む。
いくつかの例示的な実施形態では、同祖性推定は、サンプルペアの間のゲノム全域にわたる0、1、及び2のIBD値の計算を含む。
いくつかの例示的な実施形態では、遺伝子型尤度スコアは、複数の核家族における複数のヒト被験者由来のDNA配列サンプルに基づいている。
いくつかの例示的な実施形態では、本非一時的なコンピュータ実行方法は、遺伝子型尤度スコアに基づいて、バリアントがde novo変異である確率が計算された後に、バリアントをフィルタ処理するために、データプロセッサを使用することをさらに含む。
いくつかの例示的な実施形態では、本非一時的なコンピュータ実行方法は、どちらの親サンプルにも存在しない子サンプルのバリアントを単純に同定することに基づいて、バリアントがde novo変異である確率が計算された後に、バリアントをフィルタ処理するために、データプロセッサを使用することをさらに含む。
いくつかの例示的な実施形態では、バリアントをフィルタ処理することは、子サンプルに約35未満の遺伝子型クオリティ(GQ)アノテーションを有するバリアント、またはサンプル全体にわたって10以上の代替アレルカウント(AC)を有するバリアント、または子サンプルに、約7未満のリードデプス(DP)と約4未満の代替DPとを有するバリアント、またはどちらか一方の親サンプルに約2%を超えるアレルバランス(AB)を有するバリアント、または子サンプルに約15%未満のアレルバランス(AB)を有するバリアント、または子サンプルに約90%を超えるABを有するバリアント、またはどちらか一方の親サンプルに代替アレルホモ接合性を有するバリアント、またはこれらの組み合わせのバリアントを除去することを含む。
いくつかの例示的な実施形態では、本非一時的なコンピュータ実行方法は、バリアントに品質管理メトリクスでアノテーションを付けるために、データプロセッサを使用することをさらに含む。
いくつかの例示的な実施形態では、本非一時的なコンピュータ実行方法は、どちらの親サンプルにも存在しない子サンプルのバリアントを単純に同定することに基づいて、推定de novo変異が同定された後に、サンプルBAMファイルデータに基づいてバリアントをフィルタ処理するために、データプロセッサを使用することをさらに含む。
いくつかの例示的な実施形態では、本非一時的なコンピュータ実行方法は、推定de novo変異を、推定de novo変異のアレルバランスが、子サンプルで約0.15以上の場合に、中信頼度de novo変異として分類するために、データプロセッサを使用することをさらに含む。
いくつかの例示的な実施形態では、本非一時的なコンピュータ実行方法は、推定de novo変異を、推定de novo変異のアレルバランスが、各親サンプルで約0.02以下の場合に、中信頼度de novo変異として分類するために、データプロセッサを使用することをさらに含む。
いくつかの例示的な実施形態では、本非一時的なコンピュータ実行方法は、推定de novo変異を、推定de novo変異のマッピングクオリティが約40未満ではない場合に、中信頼度de novo変異として分類するために、データプロセッサを使用することをさらに含む。
いくつかの例示的な実施形態では、本非一時的なコンピュータ実行方法は、推定de novo変異を、推定de novo変異のクオリティバイデプス(quality by depth)(QD)の値が約2未満ではない場合に、中信頼度de novo変異として分類するために、データプロセッサを使用することをさらに含む。
いくつかの例示的な実施形態では、本非一時的なコンピュータ実行方法は、推定de novo変異を、推定de novo変異のMACがサンプル全体にわたって約20未満の場合に、中信頼度de novo変異として分類するために、データプロセッサを使用することをさらに含む。
いくつかの例示的な実施形態では、本非一時的なコンピュータ実行方法は、推定de novo変異を、推定de novo変異のソフトクリップリードが、推定de novo変異のキャリアにおけるバリアントサイトで約3以下の場合に、中信頼度de novo変異として分類するために、データプロセッサを使用することをさらに含む。
いくつかの例示的な実施形態では、本非一時的なコンピュータ実行方法は、推定de novo変異を、推定de novo変異が、約4を超えるモノポリマーランを持つINDELではない場合に、中信頼度de novo変異として分類するために、データプロセッサを使用することをさらに含む。
いくつかの例示的な実施形態では、本非一時的なコンピュータ実行方法は、推定de novo変異を、推定de novo変異のアレルバランス(AB)が、子サンプルで約0.15以上、かつ各親サンプルで約0.02以下の場合、及び推定de novo変異のマッピングクオリティ(MQ)が約40未満ではない場合、及び推定de novo変異のクオリティバイデプス(quality by depth)(QD)の値が約2未満ではない場合、及び推定de novo変異のマイナーアレルカウント(MAC)がサンプル全体にわたって約20未満の場合、及び推定de novo変異のソフトクリップリードが、推定de novo変異のキャリアにおけるバリアントサイトで約3以下の場合、及び推定de novo変異が、約4を超えるモノポリマーランを持つINDELではない場合に、中信頼度de novo変異として分類するために、データプロセッサを使用することをさらに含む。
いくつかの例示的な実施形態では、本非一時的なコンピュータ実行方法は、中信頼度de novo変異を、中信頼度de novo変異が、親サンプルにおいて約90以上の遺伝子型クオリティ(GQ)アノテーションを有する場合、及び中信頼度de novo変異が、各親サンプルにおいて約10以上のリードデプス(DP)を有する場合、及び中信頼度de novo変異が、子サンプルにおいて約7以上の代替DPを有する場合、及び中信頼度de novo変異が、約3より大きいQDのSNPを有する場合、及び中信頼度de novo変異が、約5より大きいQDのINDELを有する場合に、高信頼度de novo変異として分類するために、データプロセッサを使用することをさらに含む。
いくつかの例示的な実施形態では、DNA配列サンプルはエクソーム配列を含む。
いくつかの例示的な実施形態では、複数のヒト被験者は、10Kを超える被験者を含む。
いくつかの例示的な実施形態では、低クオリティサンプルのD-statは、KS検定を使用して、サンプルの実際のアレルバランスの分布を、予想されるアレルバランスの分布と比較することによって判定される。
別の態様では、本開示は、システムを提供する。本システムは、例えば、本方法及び非一時的なコンピュータ実行方法を実行するために使用されてもよい。本システムは、概して、データプロセッサ、データプロセッサに結合されたメモリ、及びメモリに格納されたプログラムを含み、プログラムが、複数のヒト被験者由来のDNA配列サンプルの中のバリアントを同定する、同定されたバリアントに基づいて被験者の祖先スーパークラスの指定を確立する、祖先スーパークラス内での被験者の第1の同祖性推定を生成すること、被験者の祖先スーパークラスから独立した被験者の第2の同祖性推定を生成する、第2の同祖性推定の1つ以上に基づいて被験者を第1の第一度家族ネットワークにクラスタリングする、第1の第一度家族ネットワーク内での被験者の第3の同祖性推定を生成する、第1の同祖性推定と第3の同祖性推定とを統合して、統合同祖性推定を取得する、統合同祖性推定に基づいて核家族を構築する、核家族におけるバリアントを同定する、構築された核家族のトリオの両親及び子に由来するサンプルのバリアントに遺伝子型尤度スコアを付与して、バリアントがde novo変異である確率を計算するとともに、独立して、トリオのどちらの親サンプルにも存在しない子サンプルのバリアントを単純に同定して、そのバリアントがde novo変異である確率を計算し、次いで両方の確率を組み合わせ、それによって推定de novo変異のデータセットを形成するための命令を含む。
いくつかの例示的な実施形態では、プログラムは、被験者の祖先スーパークラスの指定が確立される前に、同定されたバリアントをフィルタ処理するための命令を含む。
いくつかの例示的な実施形態では、プログラムは、被験者の第2の同祖性推定が生成される前に、同定されたバリアントをフィルタ処理するための命令を含む。
いくつかの例示的な実施形態では、バリアントをフィルタ処理することは、複数のヒト被験者由来のサンプル全体にわたって約10%を超える代替アレル頻度を有するバリアントを除去すること、またはp値が約10-6超である、ハーディ・ワインベルグ平衡(HWE)に反するバリアントを除去すること、または複数のヒト被験者由来のサンプルの約5%超でコールが欠損するバリアントを除去すること、またはこれらの組み合わせを含む。
いくつかの例示的な実施形態では、プログラムは、同定されたバリアントがフィルタ処理された後に、低クオリティサンプルを除去するための命令を含む。
いくつかの例示的な実施形態では、低クオリティサンプルは、D-statが0.12を超えるか、もしくは20xのリードカバレッジが75%未満であるか、またはその両方であるサンプルである。
いくつかの例示的な実施形態では、第1の同祖性推定と第3の同祖性推定とを統合することは、第3の同祖性推定に特有のペアワイズ同祖性推定を補って、第1の同祖性推定を増強することを含む。
いくつかの例示的な実施形態では、同祖性推定は、サンプルペアの間のゲノム全域にわたる0、1、及び2のIBD値の計算を含む。
いくつかの例示的な実施形態では、遺伝子型尤度スコアは、複数の核家族における複数のヒト被験者由来のDNA配列サンプルに基づいている。
いくつかの例示的な実施形態では、プログラムは、遺伝子型尤度スコアに基づいて、バリアントがde novo変異である確率が計算された後に、バリアントをフィルタ処理するための命令を含む。
いくつかの例示的な実施形態では、プログラムは、どちらの親サンプルにも存在しない子サンプルのバリアントを単純に同定することに基づいて、バリアントがde novo変異である確率が計算された後に、バリアントをフィルタ処理するための命令を含む。
いくつかの例示的な実施形態では、バリアントをフィルタ処理することは、子サンプルに約35未満の遺伝子型クオリティ(GQ)アノテーションを有するバリアント、またはサンプル全体にわたって10以上の代替アレルカウント(AC)を有するバリアント、または子サンプルに、約7未満のリードデプス(DP)と約4未満の代替DPとを有するバリアント、またはどちらか一方の親サンプルに約2%を超えるアレルバランス(AB)を有するバリアント、または子サンプルに約15%未満のアレルバランス(AB)を有するバリアント、または子サンプルに約90%を超えるABを有するバリアント、またはどちらか一方の親サンプルに代替アレルホモ接合性を有するバリアント、またはこれらの組み合わせのバリアントを除去することを含む。
いくつかの例示的な実施形態では、プログラムは、バリアントに品質管理メトリクスでアノテーションを付けるための命令を含む。
いくつかの例示的な実施形態では、プログラムは、どちらの親サンプルにも存在しない子サンプルのバリアントを単純に同定することに基づいて、推定de novo変異が同定された後に、サンプルBAMファイルデータに基づいてバリアントをフィルタ処理するための命令を含む。
いくつかの例示的な実施形態では、プログラムは、推定de novo変異を、推定de novo変異のアレルバランスが、子サンプルで約0.15以上の場合に、中信頼度de novo変異として分類するための命令を含む。
いくつかの例示的な実施形態では、プログラムは、推定de novo変異を、推定de novo変異のアレルバランスが、各親サンプルで約0.02以下の場合に、中信頼度de novo変異として分類するための命令を含む。
いくつかの例示的な実施形態では、プログラムは、推定de novo変異を、推定de novo変異のマッピングクオリティが約40未満ではない場合に、中信頼度de novo変異として分類するための命令を含む。
いくつかの例示的な実施形態では、プログラムは、推定de novo変異を、推定de novo変異のクオリティバイデプス(quality by depth)(QD)の値が約2未満ではない場合に、中信頼度de novo変異として分類するための命令を含む。
いくつかの例示的な実施形態では、プログラムは、推定de novo変異を、推定de novo変異のMACがサンプル全体にわたって約20未満の場合に、中信頼度de novo変異として分類するための命令を含む。
いくつかの例示的な実施形態では、プログラムは、推定de novo変異を、推定de novo変異のソフトクリップリードが、推定de novo変異のキャリアにおけるバリアントサイトで約3以下の場合に、中信頼度de novo変異として分類するための命令を含む。
いくつかの例示的な実施形態では、プログラムは、推定de novo変異を、推定de novo変異が、約4を超えるモノポリマーランを持つINDELではない場合に、中信頼度de novo変異として分類するための命令を含む。
いくつかの例示的な実施形態では、プログラムは、推定de novo変異を、推定de novo変異のアレルバランス(AB)が、子サンプルで約15%以上、かつ各親サンプルで約2%以下の場合、及び推定de novo変異のマッピングクオリティ(MQ)が約40未満ではない場合、及び推定de novo変異のクオリティバイデプス(quality by depth)(QD)の値が約2未満ではない場合、及び推定de novo変異のマイナーアレルカウント(MAC)がサンプル全体にわたって約20未満の場合、及び推定de novo変異のソフトクリップリードが、推定de novo変異のキャリアにおけるバリアントサイトで約3以下の場合、及び推定de novo変異が、約4を超えるモノポリマーランを持つINDELではない場合に、中信頼度de novo変異として分類するための命令を含む。
いくつかの例示的な実施形態では、プログラムは、中信頼度de novo変異を、中信頼度de novo変異が、親サンプルにおいて約90以上の遺伝子型クオリティ(GQ)アノテーションを有する場合、及び中信頼度de novo変異が、各親サンプルにおいて約10以上のリードデプス(DP)を有する場合、及び中信頼度de novo変異が、子サンプルにおいて約7以上の代替DPを有する場合、及び中信頼度de novo変異が、約3より大きいQDのSNPを有する場合、及び中信頼度de novo変異が、約5より大きいQDのINDELを有する場合に、高信頼度de novo変異として分類するための命令を含む。
いくつかの例示的な実施形態では、DNA配列サンプルはエクソーム配列を含む。
いくつかの例示的な実施形態では、複数のヒト被験者は、10Kを超える被験者を含む。
いくつかの例示的な実施形態では、低クオリティサンプルのD-statは、KS検定を使用して、サンプルの実際のアレルバランスの分布を、予想されるアレルバランスの分布と比較することによって判定される。
いくつかの例示的な実施形態では、本方法、非一時的なコンピュータ実行方法、またはシステムは、構築された核家族のトリオの両親及び子に由来するサンプルのバリアントに遺伝子型尤度スコアを付与して、バリアントがde novo変異である確率を計算し、バリアントがde novo変異である確率が有意に高いバリアントを選択するとともに、独立して、トリオのどちらの親サンプルでもコールされない子サンプルでコールされたバリアントを単純に同定し、次いでde novo変異の2つのセットを組み合わせ、それによって推定de novo変異のデータセットを形成することを含む。
別の態様では、本開示は、ヒト集団における近縁性の予測モデルを提供する。予測モデルは、第1の集団データセットを確立すること、120年のバーンインフェーズを実行して、第2の集団データセットを確立すること、及び以下のステップを実施することによって、第2の集団データセットを修正することを含むプロセスによって準備し得る。(a)第2の集団データセット内の個人を、個人の年齢に応じて年齢プールに移動し、(b)第2の集団データセット内の独身男性及び独身女性から、いとこよりも遠縁である独身男性及び独身女性のペアをランダムに選択し、年齢パラメータによって指定される結婚によって結婚させ、結婚率パラメータによって指定された数に達するまでペアを選択する。(c)指定された離婚率で夫婦を離婚させる。夫婦は第2の集団データセットから無作為に選ばれ、離婚時に独身とマークされる。(d)指定された比率で2番目の集団データセットから無作為に独身男性及び独身女性または既婚カップルのペアを選択し、成功した受胎の目標数に達するまで、指定された出生率に従ってそれらを再現できるようにする。両親は、第1のいとこよりも遠縁であることに制限されており、第2の集団データセットの全ての個人は1年に1人の子供を持つことに制限される。(e)第2の集団データセットの個人が、特定の死亡率及び年齢パラメータにより、特定の死亡率で死亡することを許可する。(f)個人が第2の集団データセットに出入りすることを許可し、これにより、集団の年齢と性別の分布、及び第2の集団データセット内の結婚している生殖可能な年齢の者の割合が維持される。(g)個人が第2の集団データセット内を移動できるようにし、これにより、サブ集団から個人がランダムに選択され、サブ集団間の指定された移動率が達成されるまで、存在する場合は、無作為に別の集団に割り当てられる。ステップ(a)~(g)を1年間隔で所定の年数繰り返し、ステップは前の繰り返しから得られた集団データセットに適用される。
いくつかの例示的な実施形態では、第1の集団のデータセットを確立することは、いくつかのサブ集団及びサイズを指定することをさらに含む。
いくつかの例示的な実施形態では、第1の集団のデータセットを確立することは、ゼロと最大受胎年齢との間の第1の集団データセット内の個人に年齢を割り当てることをさらに含む。
いくつかの例示的な実施形態では、最大受胎年齢は49歳である。
いくつかの例示的な実施形態では、バーンインフェーズを実行することは、第2の集団データセット内の個人の出生数及び死亡数を等しく保ち、個人の正味の移動率をゼロに保つことをさらに含む。
いくつかの例示的な実施形態では、バーンインフェーズの実行はさらに、個人が最低出生可能年齢を超える年齢になると、個人の第2の集団のデータセットを青年プールから交配プールに移動すること、及び個人が最大受胎可能年齢を超える年齢になると、個人を交配プールから高齢者プールに移動すること、ならびに個人が移住または亡くなった場合、全ての年齢層から個人を除去することを含む。
いくつかの例示的な実施形態では、出生の最低年齢は15歳、最高受胎年齢は49歳である。
別の態様では、本開示は、予測モデルを使用する方法を提供し、個人の確認はランダムに行われる。
別の態様では、本開示は、予測モデルを使用する方法を提供し、個人を確認することは、クラスタ化された方法で実行される。
いくつかの例示的な実施形態では、個人を確認することは、確認された個人間の第一度または第二度の関係、またはその両方を含む、確認された個人に関する近縁性データ及び関連統計を収集することをさらに含む。
集団の近縁性を活用することによって、集団における遺伝的バリアントのフェージング/同定が行われる例示的な実施形態のフローチャートである。 集団における複合ヘテロ接合性変異(CHM)が同定される例示的な実施形態のフローチャートである。 集団におけるde novo変異(DNM)が同定される例示的な実施形態のフローチャートである。 例示的な実施形態による、ヒト集団における近縁性の予測モデルを作成する方法のフローチャートである。 同祖性が判定される例示的な実施形態のフローチャートを表す。 同祖性が判定される例示的な実施形態のフローチャートを表す。 同祖性が判定される例示的な実施形態のフローチャートを表す。 同祖性が判定される例示的な実施形態のフローチャートを表す。 集団において複合ヘテロ接合性変異(CHM)が同定/フェージングされる例示的な実施形態のフローチャートを表す。 集団において複合ヘテロ接合性変異(CHM)が同定/フェージングされる例示的な実施形態のフローチャートを表す。 集団において複合ヘテロ接合性変異(CHM)が同定/フェージングされる例示的な実施形態のフローチャートを表す。 集団においてde novo変異(DNM)が同定される例示的な実施形態のフローチャートを表す。DNMのコール、フィルタ処理、及び信頼度ランキングのワークフロー。GQ=遺伝子型クオリティである。MACはDiscovEHRのマイナーアレルカウントである。DP=DNMサイトのリードデプスである。AD=代替アレルデプスである。AB=代替アレルバランスである。MQ=マッピングクオリティである。QD=DiscovEHRpVCFと呼ばれる接合のクオリティバイデプス(クオリティバイデプス)である。ホモポリマーINDELは、同じヌクレオチドの4つ以上の連続した塩基対を持つINDELである。ブラックリストに登録された遺伝子には、PDE4DIP、PRAMEF1、PABPC3、NBPF10、NBPF14、嗅覚遺伝子(OR)、MUC遺伝子(MUC)、及びHLA遺伝子(HLA-)が含まれる。 集団においてde novo変異(DNM)が同定される例示的な実施形態のフローチャートを表す。DNMのコール、フィルタ処理、及び信頼度ランキングのワークフロー。GQ=遺伝子型クオリティである。MACはDiscovEHRのマイナーアレルカウントである。DP=DNMサイトのリードデプスである。AD=代替アレルデプスである。AB=代替アレルバランスである。MQ=マッピングクオリティである。QD=DiscovEHRpVCFと呼ばれる接合のクオリティバイデプス(クオリティバイデプス)である。ホモポリマーINDELは、同じヌクレオチドの4つ以上の連続した塩基対を持つINDELである。ブラックリストに登録された遺伝子には、PDE4DIP、PRAMEF1、PABPC3、NBPF10、NBPF14、嗅覚遺伝子(OR)、MUC遺伝子(MUC)、及びHLA遺伝子(HLA-)が含まれる。 A、B、C及びDは、様々なタイプの集団ベースのゲノム研究、及び対応するサンプリング方法の概要を提供するスキームを表し、厳密な確認が家族構造を増加させ、使用する必要がある統計分析アプローチに影響を与えることを示す。パネルAは、(1)従来の集団ベースのゲノム研究(灰色のボックス)、(2)医療集団ベースのゲノム(HPG)研究(緑色のボックス)、及び(3)家族ベースのゲノム研究(黄色のボックス)の概略図を示す。パネルBは、前述の3つの確認アプローチにおける家族構成の線グラフを示す。パネルCは、前述の3つの確認アプローチにおける家族構造の散布図を示す(線は、前述の3つの確認アプローチから確認された第一度及び第二度のペアワイズ関係を示す)。パネルDは、家族構成のレベルに基づいて4つのカテゴリに分類された統計分析アプローチを示す。 分析されたデータセット(DiscovEHRデータセット;実施例を参照)の中の潜在複合ヘテロ接合性変異(pCHM)のフェーズを決定するために行われたカスケード分析の概要を示す例示的な実施形態のフローチャートである。 例示的な動作環境である。 開示された方法を実行するために構成された複数のシステム構成要素を示す。 A、B、C、及びDは、例示的な実施形態によるDiscovEHRコホートからの第1の61Kのシークエンシングされた個人に見られる近縁性を示す。パネルAは、IBD0対IBD1のプロットを示す。 例示的な実施形態によるDiscovEHRコホートからの第1の61Kのシークエンシングされた個人に見られる近縁性を示す。パネルBは、分析されたコホート内の第一度家族ネットワークのサイズ分布をプロットしたヒストグラムを示す。 例示的な実施形態によるDiscovEHRコホートからの第1の61Kのシークエンシングされた個人に見られる近縁性を示す。パネルCは、ペアワイズIBD推定から再構築された、シークエンシングされた25の個人を含む第一度家族ネットワークの家系図を示す。 例示的な実施形態によるDiscovEHRコホートからの第1の61Kのシークエンシングされた個人に見られる近縁性を示す。パネルDは、7,084の個人の最大の第二度家族ネットワークを表現するスキームを示す。 例示的な実施形態によって確認されたデータセット内の確認された個人の数の関数として、DiscovEHRコホート内の近縁性の累増を示す。 例示的な実施形態によって確認されたデータセット内の確認された個人の数の関数として、DiscovEHRコホート内の近縁性の累増を示す。 例示的な実施形態による、61KのDiscovEHR参加者の中にいる第一度血縁者の確認と、シミュレートされた集団のランダムな確認との比較を示す。パネルAは、第一度血縁者ペアの確認を示す。 例示的な実施形態による、61KのDiscovEHR参加者の中にいる第一度血縁者の確認と、シミュレートされた集団のランダムな確認との比較を示す。パネルBは、複数の第一度血縁者を持つ個人の数の確認を示す。 例示的な実施形態による、シミュレートされた集団及び確認されたDiscoverEHRコホート内の第一度の近縁性の累増に適合する確認を示す。パネルAは、第一度血縁者のペアの蓄積を示す。 例示的な実施形態による、シミュレートされた集団及び確認されたDiscoverEHRコホート内の第一度の近縁性の累増に適合する確認を示す。パネルBは、1人以上の第一度血縁者がいる確定した参加者の割合を示す。 例示的な実施形態による、シミュレートされた集団及び確認されたDiscoverEHRコホート内の第一度の近縁性の累増に適合する確認を示す。パネルCは、第一度の関係の数の上限と下限とともにシミュレーションされた確認予測を示す。 例示的な実施形態による、シミュレートされた集団及び確認されたDiscoverEHRコホート内の第一度の近縁性の累増に適合する確認を示す。パネルDは、第一度血縁者が1人以上いる、確認された参加者の割合の上限と下限とともにシミュレーションされた予測を示す。 例示的な実施形態による、確認され拡張されたDiscovEHRコホートからの最初の92Kのシークエンシングされた個人を示す。パネルAは、IBD0対IBD1のプロットを示す。 例示的な実施形態による、確認され拡張されたDiscovEHRコホートからの最初の92Kのシークエンシングされた個人を示す。パネルBは、分析されたコホート内の第一度家族ネットワークのサイズ分布をプロットしたヒストグラムを示す。 例示的な実施形態による、確認され拡張されたDiscovEHRコホートからの最初の92Kのシークエンシングされた個人を示す。パネルCは、ペアワイズIBD推定から再構築された、シークエンシングされた25の個人を含む第一度家族ネットワークの家系図を示す。 例示的な実施形態による、確認され拡張されたDiscovEHRコホートからの最初の92Kのシークエンシングされた個人を示す。パネルDは、7,084の個人の最大の第二度家族ネットワークを表現するスキームを示す。 例示的な実施形態による、シミュレートされた集団のランダムな確認と比較した、92Kの拡張DiscovEHR参加者の中にいる第一度血縁者の確認の間の比較を示す。パネルAは、第一度血縁者ペアの確認を示す。 例示的な実施形態による、シミュレートされた集団のランダムな確認と比較した、92Kの拡張DiscovEHR参加者の中にいる第一度血縁者の確認の間の比較を示す。パネルBは、複数の第一度血縁者を持つ個人の数の確認を示す。 例示的な実施形態による、拡張されたDiscovEHRコホートにおける第一度近縁性の蓄積をシミュレートした集団及び確認の適合を示す。パネルAは、第一度血縁者のペアの蓄積を示す。 例示的な実施形態による、拡張されたDiscovEHRコホートにおける第一度近縁性の蓄積をシミュレートした集団及び確認の適合を示す。パネルBは、1人以上の第一度血縁者がいる確定した参加者の割合を示す。 例示的な実施形態による、拡張されたDiscovEHRコホートにおける第一度近縁性の蓄積をシミュレートした集団及び確認の適合を示す。パネルCは、第一度の関係の数の上限と下限とともにシミュレーションされた確認予測を示す。 例示的な実施形態による、拡張されたDiscovEHRコホートにおける第一度近縁性の蓄積をシミュレートした集団及び確認の適合を示す。パネルDは、第一度血縁者が1人以上いる、確認された参加者の割合の上限と下限とともにシミュレーションされた予測を示す。 例示的な実施形態による、確認されたDiscovEHRコホート内の第一度及び第二度の近縁性の蓄積をシミュレートした集団及び確認の適合を示す。パネルAは、第一度及び第二度の血縁者のペアの蓄積を示す。 例示的な実施形態による、確認されたDiscovEHRコホート内の第一度及び第二度の近縁性の蓄積をシミュレートした集団及び確認の適合を示す。パネルBは、1人以上の第一度及び第二度の血縁者がいる確定した参加者の割合を示す。 例示的な実施形態による、確認されたDiscovEHRコホート内の第一度及び第二度の近縁性の蓄積をシミュレートした集団及び確認の適合を示す。パネルCは、第一度及び第二度の関係の数の上限と下限とともにシミュレーションされた確認予測を示す。 例示的な実施形態による、確認されたDiscovEHRコホート内の第一度及び第二度の近縁性の蓄積をシミュレートした集団及び確認の適合を示す。パネルDは、第一度または第二度の血縁者が1人以上いる、確認された参加者の割合の上限と下限とともにシミュレーションされた予測を示す。 例示的な実施形態による、確認された拡張DiscovEHRコホート内の第一度及び第二度の近縁性の蓄積をシミュレートした集団及び確認の適合を示す。パネルAは、第一度及び第二度の血縁者のペアの蓄積を示す。 例示的な実施形態による、確認された拡張DiscovEHRコホート内の第一度及び第二度の近縁性の蓄積をシミュレートした集団及び確認の適合を示す。パネルBは、1人以上の第一度及び第二度の血縁者がいる確定した参加者の割合を示す。 例示的な実施形態による、確認された拡張DiscovEHRコホート内の第一度及び第二度の近縁性の蓄積をシミュレートした集団及び確認の適合を示す。パネルCは、第一度及び第二度の関係の数の上限と下限とともにシミュレーションされた確認予測を示す。 例示的な実施形態による、確認された拡張DiscovEHRコホート内の第一度及び第二度の近縁性の蓄積をシミュレートした集団及び確認の適合を示す。パネルDは、第一度または第二度の血縁者が1人以上いる、確認された参加者の割合の上限と下限とともにシミュレーションされた予測を示す。 例示的な実施形態による、DiscovEHRコホートにおける個人ごと同定された複合ヘテロ接合性変異(CHM)の数を示す。パネルAは、DiscovEHRコホートにおける個人あたりのCHMの数を示す。 例示的な実施形態による、DiscovEHRコホートにおける遺伝子ごとに同定された複合ヘテロ接合性変異(CHM)の数を示す。パネルBは、DiscovEHRコホートにおける遺伝子あたりのCHMの数を示す。 例示的な実施形態による、DiscovEHRコホートにおける個人ごとに同定されたde novo変異(DNM)の数を示す。パネルCは、DiscovEHRコホートのトリオの子供間のエクソン高信頼度DNMの数の分布を示す。 例示的な実施形態による、DiscovEHRコホートにおける遺伝子ごとに同定されたde novo変異(DNM)の数を示す。パネルDは、遺伝子ごとのいくつかの非同義のDNMを示す。 例示的な実施形態による、DiscovEHRデータセットについて同定されたフェージング済みの複合ヘテロ接合性変異体(CHM)バリアント間のゲノム距離の範囲を示すチャートである。 (A)大動脈瘤、(B)QT延長症候群、及び(C)甲状腺癌のバリアントを含む、既知の疾患の原因となるバリアントの分離を示すDiscovEHRコホートから再構成された家系図を示す。 LDLRにおける縦列重複の22/29のキャリアを含む、シークエンシングされたDiscovEHRから、およびシークエンシングされたコホートからの影響を受けていない10人の血縁関係(第一度または第二度)のある個人から再構成された家系図である。 92Kディスカバー参加者間の潜在複合ヘテロ接合性変異(pCHM)のフェーズを決定するための例示的な実施形態の決定カスケードである。 例示的な実施形態による、同定された複合ヘテロ接合性変異(CHM)及びde novo変異(DNM)についての結果となる拡張DiscovEHRコホートを示す。パネルAは、DiscovEHRコホートにおける個人あたりのCHM数の分布を示す。 例示的な実施形態による、同定された複合ヘテロ接合性変異(CHM)についての結果となる拡張DiscovEHRコホートを示す。パネルBは、遺伝子あたりのCHM数の分布を示す。 例示的な実施形態による、de novo変異(DNM)についての結果となる拡張DiscovEHRコホートを示す。パネルCは、DiscovEHRコホートのトリオの子供たちの間の3,415個のエクソン高信頼度及び中信頼度DNMの分布を示す。 例示的な実施形態による、de novo変異(DNM)についての結果となる拡張DiscovEHRコホートを示す。パネルDは、1つ以上の2,802遺伝子にわたる非同義DNMの分布を示す。 例示的な実施形態による、拡張DiscovEHRについて同定されたフェージング済みの複合ヘテロ接合性変異体(CHM)バリアント間のゲノム距離の範囲を示すチャートである。 例示的な実施形態による、拡張DiscovEHRデータセットで識別された、家族トリオの数、親の年齢を伴う家族のトリオ、1以上のエクソンDNMを伴う発端者、エクソンDNM、中/高信頼度DNM、単一ヌクレオチドDNM、中/高信頼度バリアント、及びランダムなバリアントを示すコホートプロファイルである。 例示的な実施形態による、拡張DiscovEHRコホートにおける信頼度レベルごと及びヒトごとに同定されたDNMの数を示す。パネルAは、拡張DiscovEHRコホートにおける信頼度レベルごとのDNM数の分布を示す。パネルBは、例示的な実施形態による、DNMが同定された拡張DiscovEHRコホートにおける個人あたりのDNM数分布を示す。 例示的な実施形態による、拡張DiscovEHRコホートにおける機能的効果DNMごとに同定されたDNM数分布をプロットするヒストグラムである。 例示的な実施形態による、拡張DiscovEHRコホートにおけるDNMのタイプ(トランジション、トランスバージョン、及びインデル)ごとに同定されたDNM数分布をプロットするヒストグラムである。 例示的な実施形態による、拡張DiscovEHRコホートにおける単一ヌクレオチドDNMのタイプ(→)ごとに同定されたDNM数分布をプロットするヒストグラムである。 例示的な実施形態による、拡張DiscovEHRコホートにおける染色体あたり10Mのエクソン塩基対ごとに同定されたDNM数分布をプロットするヒストグラムである。 例示的な実施形態による、拡張DiscovEHRコホートにおいてCGジヌクレオチド(慣習的にCpGと記され、「p」は2つの塩基間のリン酸塩を表す)が富化されていることがわかっているゲノムの領域で発生するDNMまたはランダムに選択されたバリアントのパーセンテージの分布をプロットする棒グラフである。 LDLRにおける新規FHを引き起こす縦列重複の25/37のキャリア、及び20の非キャリア、血縁関係(第一度または第二度)個人を含む、シークエンシングされた拡張DiscovEHRコホートから再構成された家系予測の画像を示す。 例示的な実施形態による、DiscovEHRコホートにおける出生時の父方及び母方の年齢と、子供において同定されたエクソンDNMの数との関係を示す。 例示的な実施形態による、DiscovEHRコホートにおける出生時の父方及び母方の年齢と、子供において同定されたエクソンDNMの数との関係を示す。 例示的な実施形態による、DiscovEHRコホートにおける子供の出産時の母親及び父親の年齢と、子供において同定されたDNMとの相関を示すチャートである。 例示的な実施形態による、拡張DiscovEHRコホートにおいて同定されたDNM及びランダムバリアントの病原性予測をプロットするヒストグラムである。
「a」という用語は、「少なくとも1つ」を意味すると解されるべきである。「約」及び「およそ」という用語は、当業者によって理解されるであろう標準的な変動を許すと解されるべきである。範囲が提供される場合は、端点が含まれる。
これまでの大規模なヒトゲノム研究では、通常は、多数の様々な地理的地域及び/またはヘルスケアシステムにわたってヒトサンプルを収集し、それらを組み合わせて分析用のコホートを生成した。このようなコホートでは、サンプリングされる個人の総数はしばしば大きいが、これらのコホートにおける近縁性及び家族構造の範囲は比較的小さい傾向があった。ゲノム解析との関連で一般に用いられる多くの統計的手法では、関連解析及び主成分分析を含めて、全てのサンプルが血縁関係にないことを求められる。そうでない場合は、これらの検定の統計出力にバイアスがかかり、結果的にp値をつり上げ、偽陽性所見をもたらすことになる(図8)(Kang et al.(2010),Nature Publishing Group 42,348-354、Sun and Dimitromanolakis(2012),Methods Mol.Biol.850,47-57、Devlin and Roeder(1999),Biometrics 55,997-104、及びVoight and Pritchard(2005),PLoS Genet 1,e32-10)。
近縁関係にあるサンプルがデータセットに少数しかない場合、データセットから家族構造を除去することは実行可能な選択肢である(Lek,et al.(2016),Nature Publishing Group 536,285-291、Fuchsberger et al.(2016),Nature Publishing Group 536,41-47、Locke et al.(2015),Nature 518,197-206、及びSurendran et al.(2016)Nat Genet 48,1151-1161)。また、血縁関係にないデータのサブセットが、主成分(PC)を算出し、次いで残りのサンプルをこれらのPCに投影するなどの統計分析のために十分である場合には、家族構造の除去はやはり実現可能な選択肢となる(Dewey et al.(2016),Science 354,aaf6814-aaf6814)。血縁関係にない個人の最大サイズのセットを研究者が保持するのに役立ついくつかの方法が存在する(Staples at al.(2013),Genet.Epidemiol.37,136-141、Chang at al.(2015),Gigascience 4,7)。残念ながら、血縁個人を除去するとサンプルサイズが小さくなるだけでなく、有用な関係の情報が破棄される。実際のところ、データセットが中程度の家族構造を有する場合であっても、そのような情報の損失は多くの分析では受け入れられない。
本開示は、複数対象者のゲノムサンプルのデータセットのうち家族構造及び家系構造ならびに近縁性に関する情報が、とりわけ、稀少遺伝的バリエーション(例えば、複合ヘテロ接合性変異及び/またはde novo変異)と病気との間のつながりを調査することを可能にする多くの分析の機会を与えることから、その情報は有用であるという認識に、少なくとも部分的に基づいている。
また、本開示は、ゲノムワイドな同祖性(IBD)の推定が、複数対象者のゲノムサンプルのデータセット内の近縁性、及び2組の個人間の近縁性のレベルを定量化するための優れたメトリックであるという認識にも、少なくとも部分的に基づいている。
正確なペアワイズ関係をモデル化するいくつかの統計的手法が開発されている。例えば、混合モデルを用いるゲノムワイド関連解析は、より高性能であり、交絡近縁性をモデル化しない方法より性能が優れているが(Kang et al.(2010),Nature Publishing Group 42,348-354、Zhang et al.(2010),Nat Genet 42,355-360、Yang et al.(2014),Nat Genet 46,100-106、Kirkpatrick and Bouchard-Cote(2016),arXiv q-bio.QM)、混合モデルは、家族構造の中に含まれる情報を完全には活用せず、数十万のサンプルと数百から数千の表現型とを含むデータセットに合わせて、実際のところ、拡大することはできない。ペアワイズ関係を家系を用いないQTL連鎖解析に使用する場合もある(Day-Williams et al.(2011),Genet.Epidemiol.35,360-370)。ペアワイズ関係推定(PCrelate)(Conomos et al.(2016),Am.J.Hum.Genet.98,127-148)及び主成分分析(PC-AiR)(Conomos et al.(2015),Genet.Epidemiol.39,276-293)のために、集団構造及び家族構造をモデル化する更なるソフトウェアパッケージがある。
従来のゲノムワイド関連解析とは対照的に、最近及び将来の大規模ゲノム研究、例えば、本開示に具体化される研究では、個別の地理的地域から、数万~数十万人もの参加者をサンプリングする。結果的に、これらの研究は、同じ地理的地域から遥かに大きな割合の人々、したがって、サンプルデータセット内の家族構造及び家系構造を確認して、従来の集団全体の関連解析では正しく評価されていない家族において分離する稀少バリアントを同定する。
いくつかの理由によって、そのような大規模なゲノム研究のデータには、家族構造及び遠縁の隠れた近縁性が豊富に含まれる。第一に、これらの研究は、特定の地理的地域から、例えばヘルスケアシステム集団を通じて大量にサンプリングし、判明される血縁関係にある個人の対の数は、単一集団からより多くのサンプルが確認されるにつれて、組み合わせ的に増加する(図8A)。第二に、同じ地理的地域に住んでいる家族は、共通の保険適用及び利便性のために、同じヘルスケアシステムにおいて同じ医師から医療を受けているものと考えられる。第三に、共通の遺伝的要因及び環境的要因は、特定の家族を対象にして医療連携の頻度を増加させる場合がある。移住率の低い集団では、家族構造及び遠縁の隠れた近縁性の両方とも、さらに一層顕著である(Henn et al.(2012),PLoS ONE 7,e34267)。家族構造の影響力は、連鎖、家系に基づく分析、IBDモデル化、及び非血縁関係者の分析に対するサンプリング方法の効果を通して観察され得る(図8、パネルD)。「連鎖」とは、情報価値のある1つ以上の家系を用いる従来の連鎖解析のことをいう。「家系に基づく分析」とは、血縁関係がない個人を含むより大きなコホート内の家系構造を用いる連鎖の範囲を超えた統計的手法のことをいう。「IBDモデル化」とは、家系構造全体を用いずに個人間のペアワイズ関係をモデル化する分析のことをいう。「非血縁関係者の分析」とは、コホート内の全ての個人が非血縁関係にあると仮定する分析のことをいう。
本開示では、家族構造に焦点を合わせ、実データ及び模擬データの両方を使用して、高レベルの家族構造を明示する。本開示の改良点の1つは、複合ヘテロ接合性変異(CHM)及び/またはde novo変異(DNM)を、従来型のアプローチよりも正確かつ確実に同定し及び/またはフェージングすることである(実施例の節に開示されているデータ参照)。
したがって、本開示は、確認された集団における遺伝的バリアントを、集団の近縁性を活用することによってフェージングする方法を提供する。例示的なフェージング方法を概説するフローが図1に与えられる。
本方法は、様々な集団における各種遺伝的バリアントに適用し得る。評価し得る遺伝的バリアントの種類の非限定的な例には、点突然変異、挿入、欠失、逆位、重複、及び多量体化が含まれる。集団の種類の非限定的な例には、単一医療ネットワーク集団、複数医療ネットワーク集団、人種的、文化的、もしくは社会的に同種もしくは異種の集団、混合年齢集団もしくは年齢に関して均質な集団、地理的に集中した集団もしくは地理的に分散した集団、またはこれらの組み合わせが含まれる。遺伝的バリアントを取得し得る手段の非限定的な例は、次のステップを含む。
- サンプルの調製及びシークエンシング(Dewey et al.(2016),Science 354,aaf6814-1~aaf6814-10)。
- シークエンシングの完了時に、各シークエンシングランからの生データがローカルバッファストレージに収集され、自動解析のために、DNAnexusプラットフォーム(Reid et al.(2014);BMC Bioinformatics 15,30)にアップロードされる。
- サンプルレベルのリードファイルが、CASAVA(Illumina Inc.,San Diego,CA)で生成され、BWA-memでGRCh38(Li及びDurbin(2009);Bioinformatics 25,1754-176;Li(2013);arXiv q-bio.GN)にアライメントされる。
- 結果として得られたBAMファイルが、ソート、重複のマーク付け、及び推定インデル周辺のリードの局所的なアライメント修正の実行のために、GATK(McKenna et al.(2010);Genome Res.20,1297-1303)及びPicardを使用して処理される。
- シークエンシングされたバリアントは、転写産物及び遺伝子に対する機能的影響を明らかにするために、Ensembl85遺伝子定義を使用して、snpEFF(Cingolani et al.(2012);Fly(Austin)6,80-92)でアノテーションを付けられる。
本方法は、上記のどのステップにも限定されず、配列バリアントの取得を、任意好適な手段によって行い得ると解される。
図1は、集団の近縁性を活用することによって、集団における遺伝的バリアントのフェージング/同定が行われる例示的な実施形態のフローチャートである。複数のヒト被験者から得られた核酸配列サンプルのデータセットから、低クオリティの配列バリアントを、任意好適な手段によってステップ1で除去し得る。そのような手段の非限定的な例としては、PLINK(Chang et al.(2015);Gigascience 4,7)、及び本実施例に開示されているものが挙げられる。
サンプルの1つ以上のそれぞれについて祖先スーパークラスの指定を、任意好適な手段によって、ステップ2で確立し得る。そのような手段の非限定的な例としては、PLINK(Chang et al.(2015);Gigascience 4,7)及び本実施例に開示されているものが挙げられる。
低クオリティサンプルを、任意好適な手段によって、ステップ3でデータセットから除去し得る。そのような手段の非限定的な例としては、Dewey et al.(2016),Science 354,aaf6814-1~aaf6814-10に開示されているもの、及び本実施例に開示されているものが挙げられる。
祖先スーパークラス内での被験者の第1の同祖性推定を、任意好適な手段によって、ステップ4で生成し得る。そのような手段の非限定的な例としては、PLINK(Chang et al.(2015);Gigascience 4,7)、及び本実施例に開示されているものが挙げられる。
被験者の祖先スーパークラスから独立した被験者の第2の同祖性推定を、ステップ5及びステップ6で生成し得、第2の同祖性推定の1つ以上に基づいて、任意好適な手段によって、被験者を第1の第一度家族ネットワークにクラスタリングし得る。そのような手段の非限定的な例としては、PLINK(Chang et al.(2015);Gigascience 4,7)、及び本実施例に開示されているものが挙げられる。
第1の第一度家族ネットワーク内での被験者の第3の同祖性推定を、任意好適な手段によって、ステップ7で生成し得る。そのような手段の非限定的な例としては、PLINK(Chang et al.(2015);Gigascience 4,7)、及び本実施例に開示されているものが挙げられる。
第1の同祖性推定と第3の同祖性推定とを、任意好適な手段によってステップ8で統合して、統合同祖性推定を取得し得る。そのような手段の非限定的な例としては、PLINK(Chang et al.(2015);Gigascience 4,7)、及び本実施例に開示されているものが挙げられる。
統合同祖性推定に基づいて被験者の第2の第一度家族ネットワークを、任意好適な手段によって、ステップ9で構築し得る。そのような手段の非限定的な例としては、PLINK(Chang et al.(2015);Gigascience 4,7)、及び本実施例に開示されているものが挙げられる。
ステップ10で、バリアントを、任意好適な手段によって、統合同祖性推定及び第2の第一度家族ネットワークに応じて、複合ヘテロ接合性変異(CHM)であるかもしくはCHMではないものとしてフェージングし得、またはバリアントを、任意好適な手段によって、統合同祖性推定及び第2の第一度家族ネットワークに応じて、de novo変異(DNM)として同定し得る。そのような手段の非限定的な例としては、図6及び図7ならびに本実施例に開示されているものが挙げられる。
図5~図7は、同祖性(IBD)推定を生成し、そのIBD推定を使用して、複合ヘテロ接合性変異(CHM)もしくは潜在複合ヘテロ接合性変異(pCHM)、またはde novo変異(DNM)として遺伝子バリアントをフェージングする手法を、限定するためではなく説明のために、基本的な動作ロジックを示す。ロジック中で識別されるプログラム(例えば、EAGLE、PLINKなど)は、それらが識別されるステップの例となるものであるが、そのようなプログラムは、そのようなステップを実行する唯一の方法ではないことが理解される。
バリアントを複合ヘテロ接合性変異(CHM)としてフェージングすることは、(1)集団アレル頻度に応じてバリアントをフェージングすること、(2)ハーディ・ワインベルグ平衡(HWE)から外れたバリアント、もしくは同じサンプル内の別のバリアントの10塩基対内のバリアント、またはその両方を除去すること、及びクオリティバイデプス(quality by depth)(QD)が約2以下、またはリードデプス(DP)が約5未満、または代替アレルバランス(AB)が約10%以下、またはこれらの組み合わせのSNPを除去すること、及びQDが約2以下、またはDPが約5未満、またはABが約10%以下、またはこれらの組み合わせの挿入または欠失(INDEL)を除去すること、(3)同じサンプル内及び同じ遺伝子内に1つ以上のバリアントの対が存在する潜在複合ヘテロ接合性変異(pCHM)として残りのバリアントを選択すること、ならびに(4)pCHMをシスpCHMまたはトランスpCHMとしてフェージングし、次いでトランスpCHMとしてフェージングされたpCHMをCHMとして分類することを含み得る。集団アレル頻度に応じてバリアントをフェージングすることは、EAGLE(Loh et al.(2016),Nat Genet 48,1443-1448)を含むがこれに限定されない任意好適な手段によって促進し得る。ある一定の選択基準を満たさないバリアントは除去され得、残りのバリアントが潜在複合ヘテロ接合性変異として選択され、潜在複合ヘテロ接合性変異は、本実施例に記載されたものを含む任意好適な手段によってフェージングされる。また、これらの例示的な実施形態は図6にも示される。
バリアントを複合ヘテロ接合性変異としてフェージングすることは、ハーディ・ワインベルグ平衡(HWE)から外れたバリアント、もしくは同じサンプル内の別のバリアントの10塩基対内のバリアント、またはその両方を除去すること、及びクオリティバイデプス(quality by depth)(QD)が約3以下、またはリードデプス(DP)が約7未満、または代替アレルバランス(AB)が約15%以下、またはこれらの組み合わせのSNPを除去すること、及びQDが約5以下、またはDPが約10未満、またはABが約20%以下、またはこれらの組み合わせの挿入または欠失(INDEL)を除去することを含み得る。これらのステップは、除外パラメータがより厳しいレベルに設定されることを除いて、本明細書の他の箇所に記載されるように実施してもよい。
いくつかの例示的な実施形態では、本方法は、さらに、(1)機能的効果優先度に応じてCHMにスコアを付与すること、及び(2)サンプル当たりの遺伝子ごとに最も高い機能的効果優先度スコアを有するCHMを選択することを含み、ヒトが同一遺伝子内に複数のCHMを有する場合に、タンパク質機能抑制をもたらす可能性が最も高いCHMが同定されるようにするものとする。これらのステップは、SIFT(Loh et al.(2016);Nat Genet 48,1443-1448)(損傷を与える(damaging))、PolyPhen2 HDIV45(損傷を与える(damaging)、及び損傷を与えるおそれがある(possibly damaging))、PolyPhen2 HVAR(損傷を与える(damaging)、及び損傷を与えるおそれがある(possibly damaging))、LRT46(有害(deleterious))、及びMutationTaster(Schwarz et al.(2014);Nat.Methods 11、361-362)(必然的に疾患を引き起こす(disease causing automatic)、及び疾患を引き起こす(disease causing))を含むが、これらに限定されない任意好適な手段によって行い得る。
バリアントをde novo変異としてフェージングすることは、(1)第2の第一度家族ネットワーク及びそのトリオのサンプルの中のバリアントを同定すること、(2)遺伝子型尤度スコアを、トリオ内の親サンプル及び対応する子サンプルのバリアントに割り当て、バリアントがde novo変異である確率を計算し、計算した確率が統計的に有意な場合に、そのバリアントを推定de novo変異として同定すること、(3)トリオの子サンプルのバリアントを同定し、トリオのいずれのの親サンプルにもバリアントが存在しない場合に、そのバリアントを推定de novo変異として同定すること、(4)子サンプルに約35未満の遺伝子型クオリティ(GQ)アノテーションを有する推定de novo変異、またはサンプル全体にわたって10以上の代替アレルカウント(AC)を有する推定de novo変異、または子サンプルに、約7未満のリードデプス(DP)と約4未満の代替DPとを有する推定de novo変異、またはどちらか一方の親サンプルに約2%を超えるアレルバランス(AB)を有する推定de novo変異、または子サンプルに約15%未満のアレルバランス(AB)を有する推定de novo変異、または子サンプルに約90%を超えるABを有する推定de novo変異、またはどちらか一方の親サンプルに代替アレルホモ接合性を有する推定de novo変異、またはこれらの組み合わせの推定de novo変異を除去することによって同定される推定de novo変異をフィルタ処理すること、ならびに(5)フィルタ処理済みの同定された推定de novo変異を組み合わせて、それによって推定de novo変異データセットを形成することを含み得る。これらのステップは、本実施例に記載されたものを含む任意好適な手段によって行い得る。また、これらの例示的な実施形態は図7にも示される。
いくつかの例示的な実施形態では、本方法は、さらに、推定de novo変異のアレルバランスが、子サンプルで約0.15以上、かつ各親サンプルで約0.02以下の場合、及び推定de novo変異のマッピングクオリティが約40未満ではない場合、及び推定de novo変異のクオリティバイデプス(quality by depth)(QD)の値が約2未満ではない場合、及び推定de novo変異のMACがサンプル全体にわたって約20未満の場合、及び推定de novo変異のソフトクリップリードが、推定de novo変異のキャリアにおけるバリアントサイトで約3以下の場合、及び推定de novo変異が、約4を超えるモノポリマーランを持つINDELではない場合、推定de novo変異データセット内の推定de novo変異を中信頼度de novo変異として分類することを含む。いくつかの例示的な実施形態では、本方法は、さらに、中信頼度de novo変異が、親サンプルにおいて約90以上の遺伝子型クオリティアノテーションを有する場合、及び中信頼度de novo変異が、各親サンプルにおいて約10以上のリードデプスを有する場合、及び中信頼度de novo変異が、子サンプルにおいて約7以上の代替リードデプスを有する場合、及び中信頼度de novo変異が、約3より大きいQDのSNPを有する場合、及び中信頼度de novo変異が、約5より大きいQDのINDELを有する場合に、中信頼度de novo変異を、高信頼度de novo変異として分類することを含む。どちらの例示的な実施形態も、本実施例に開示されたものを含むが、これらに限定されない、いずれかの方法で実施してもよい。
また、本開示は、集団における複合ヘテロ接合性変異(CHM)を同定する方法を提供する。CHMを同定する方法の例を示すフローチャートが図2に与えられる。
本方法は、任意手段によって得られる任意の種類のヒト被験者に由来するあらゆる種類のDNA配列サンプルに適用し得る。バリアントの非限定的な例には、点突然変異、挿入、欠失、逆位、重複、及び多量体化が含まれる。ヒト被験者の種類の非限定的な例には、単一医療ネットワーク集団、複数医療ネットワーク集団、人種的、文化的、もしくは社会的に同種もしくは異種の集団、混合年齢集団もしくは年齢に関して均質な集団、地理的に集中した集団もしくは地理的に分散した集団、またはこれらの組み合わせからのヒト被験者が含まれる。DNA配列サンプルは、Dewey et al.(2016),Science 354,aaf6814-1~aaf6814-10に開示されているものを含むが、これに限定されない多くの方法のいずれかで取得し得る。
いくつかの例示的な実施形態では、DNA配列サンプルはエクソーム配列を含む。エクソームDNAは、一般に使用される方法のいずれかによって、またはDewey et al.(2016),Science 354,aaf6814-1~aaf6814-10に記載の通りに単離し得る。
複数のヒト被験者由来のDNA配列サンプルの中のバリアントを、任意好適な手段によって、ステップ11で同定し得る。バリアントを同定し得る手段の非限定的な例は、次のステップを含む。
- シークエンシングの完了時に、各シークエンシングランからの生データがローカルバッファストレージに収集され、自動解析のために、DNAnexusプラットフォーム(Reid et al.(2014);BMC Bioinformatics 15,30)にアップロードされる。
- サンプルレベルのリードファイルが、CASAVAソフトウェア(Illumina Inc.,San Diego,CA)で生成され、BWA-memでGRCh38(Li及びDurbin(2009);Bioinformatics 25,1754-176;Li(2013);arXiv q-bio.GN)にアライメントされる。
- 結果として得られたBAMファイルが、ソート、重複のマーク付け、及び推定インデル周辺のリードの局所的なアライメント修正の実行のために、GATK(McKenna et al.(2010);Genome Res.20,1297-1303)及びPicardを使用して処理される。
- シークエンシングされたバリアントは、転写産物及び遺伝子に対する機能的影響を明らかにするために、Ensembl85遺伝子定義を使用して、snpEFF(Cingolani et al.(2012);Fly(Austin)6,80-92)でアノテーションを付けられる。
本方法は、上記のどのステップにも限定されず、配列バリアントの取得を、任意好適な手段によって行い得ると解される。
図2は、集団における複合ヘテロ接合性変異(CHM)が同定される例示的な実施形態のフローチャートである。同定されたバリアントに基づく被験者の祖先スーパークラスの指定を、任意好適な手段によって、ステップ12で確立し得る。祖先スーパークラス内での被験者の第1の同祖性推定を、任意好適な手段によって、ステップ13で生成し得る。被験者の祖先スーパークラスから独立した被験者の第2の同祖性推定を、任意好適な手段によって、ステップ14で生成し得る。ステップ15で、第2の同祖性推定の1つ以上に基づいて、任意好適な手段によって、被験者を第1の第一度家族ネットワークにクラスタリングし得る。第1の第一度家族ネットワーク内での被験者の第3の同祖性推定を、任意好適な手段によって、ステップ16で生成し得る。第1の同祖性推定と第3の同祖性推定とを、任意好適な手段によって、ステップ17で統合して、統合同祖性推定を取得し得る。ステップ18で、統合同祖性推定に基づいて、任意好適な手段によって、第2の第一度家族ネットワークを構築し得る。そのような手段の非限定的な例としては、PLINK(Chang et al.(2015);Gigascience 4,7)、及び本実施例に開示されているものが挙げられる。いくつかの例示的な実施形態では、同祖性推定は、サンプルペアの間のゲノム全域にわたる0、1、及び2のIBD値の計算を含む。
EAGLE(Loh et al.(2016),Nat Genet 48,1443-1448)を含むがこれに限定されない任意好適な手段によって、集団アレル頻度に応じてサンプルの中のバリアントを、ステップ19でフェージングし得る。
フェージングされたバリアントの対を、同じ被験者及び遺伝子における2つのそのバリアントの存在に基づいて、ステップ20で潜在CHMとして分類し得る。このことは、同じ人物の遺伝子内のへテロ接合性pLoF及び/または有害ミスセンスバリアントの全てのあり得る組み合わせを検査することによって確認された。
潜在CHMを、ステップ21でシスまたはトランスとしてフェージングし得、トランスとしてフェージングされた潜在CHMをCHMとして分類し得る。潜在CHMは、好適な手段のいずれかによってフェージングし得る。非限定的な例では、EAGLEを用いる集団アレル頻度ベースのフェージングと家系/関係ベースのフェージングとの組み合わせが使用されて、潜在CHMの相が、シスかトランスかを判定する(この例示的なプロセスは、図9にも示される)。
いくつかの例では、本方法は、被験者の祖先スーパークラスの指定が確立される前に、同定されたバリアントをフィルタ処理することをさらに含み、いくつかの例示的な実施形態では、本方法は、被験者の第2の同祖性推定が生成される前に、同定されたバリアントをフィルタ処理することをさらに含む。バリアントを、任意好適な手段によってフィルタ処理してもよい。そのような手段の非限定的な例としては、PLINK(Chang et al.(2015);Gigascience 4,7)、及び本実施例に開示されているものが挙げられる。
いくつかの例示的な実施形態では、バリアントをフィルタ処理することは、複数のヒト被験者由来のサンプル全体にわたって約10%を超える代替アレル頻度を有するバリアントを除去すること、またはp値が約10-6超である、ハーディ・ワインベルグ平衡(HWE)に反するバリアントを除去すること、または複数のヒト被験者由来のサンプルの約5%超でコールが欠損するバリアントを除去すること、またはこれらの組み合わせを含む。ある一定の選択基準を満たさないバリアントは除去され得、残りのバリアントが潜在複合ヘテロ接合性変異として選択され、潜在複合ヘテロ接合性変異は、本実施例に記載されたものを含む任意好適な手段によってフェージングされる。また、これらの例示的な実施形態は図6にも示される。
いくつかの例示的な実施形態では、本方法は、同定されたバリアントがフィルタ処理された後に、低クオリティサンプルを除去することを含む。低クオリティサンプルを、任意好適な手段によって除去し得る。そのような手段の非限定的な例としては、一般的に知られているDewey et al.(2016),Science 354,aaf6814-1~aaf6814-10に開示されているもの、及び本実施例に開示されているものが挙げられる。いくつかの例示的な実施形態では、D-statが0.12を超えるか、もしくは20xのリードカバレッジが75%未満であるか、またはその両方であるサンプルが、除去される低クオリティサンプルであるように、パラメータが調節される。
第1の同祖性推定と第3の同祖性推定とを統合することは、第3の同祖性推定に特有のペアワイズ同祖性推定を補って、第1の同祖性推定を増強することを含み得、これは、例えば限定されないが、PLINK(Chang et al.(2015);Gigascience4,7)、及び本実施例に開示される手段によって促進し得る。
いくつかの例示的な実施形態では、本方法は、バリアントが集団アレル頻度に応じてフェージングされた後に、バリアントをフィルタ処理することをさらに含み、このフェージングは、いくつかの例示的な実施形態では、ヒト被験者のDNA配列サンプルを、遺伝子間領域において、近似的に等しいサイズであり、実質的なセグメントオーバーラップ及び切断点を有するゲノムセグメントに分割することを含み得る。集団アレル頻度に応じてバリアントをフェージングすることは、EAGLE(Loh et al.(2016),Nat Genet 48,1443-1448)を含むがこれに限定されない任意好適な手段によって促進し得る。集団アレル頻度に応じてフェージングされたバリアントをフィルタ処理することは、ハーディ・ワインベルグ平衡(HWE)から外れたバリアント、もしくは同じサンプル内の別のバリアントの10塩基対内のバリアント、またはその両方を除去すること、及びクオリティバイデプス(quality by depth)(QD)が約2以下、またはリードデプス(DP)が約5未満、または代替アレルバランス(AB)が約10%以下、またはこれらの組み合わせのSNPを除去すること、及びQDが約2以下、またはDPが約5未満、またはABが約10%以下、またはこれらの組み合わせの挿入または欠失(INDEL)を除去することを含んでもよい。集団アレル頻度に応じてフェージングされたバリアントをフィルタ処理することは、ハーディ・ワインベルグ平衡(HWE)から外れたバリアント、もしくは同じサンプル内の別のバリアントの10塩基対内のバリアント、またはその両方を除去すること、及びクオリティバイデプス(quality by depth)(QD)が約3以下、またはリードデプス(DP)が約7未満、または代替アレルバランス(AB)が約15%以下、またはこれらの組み合わせのSNPを除去すること、及びQDが約5以下、またはDPが約10未満、またはABが約20%以下、またはこれらの組み合わせの挿入または欠失(INDEL)を除去することを含んでもよい。ある一定の選択基準を満たさないバリアントは除去され得、残りのバリアントが潜在複合ヘテロ接合性変異として選択され、潜在複合ヘテロ接合性変異は、本実施例に記載されたものを含む任意好適な手段によってフェージングされる。また、これらの例示的な実施形態は図6にも示される。
潜在CHMは、トリオのデータ、もしくは親子のデータ、もしくは完全同胞のデータ、もしくは遠縁の血縁者のデータ、もしくはそれらの組み合わせに基づいてフェージングされ得、またはマイナーアレルカウント(MAC)に基づいてフェージングされる、または集団アレル頻度に基づいてフェージングされる、またはそれらの組み合わせである。フェージングは、当技術分野で一般的に用いられる任意好適な方法によって促進し得る。非限定的な例では、EAGLEを用いる集団アレル頻度ベースのフェージングと家系/関係ベースのフェージングとの組み合わせが使用されて、潜在CHMをフェージングする。また、この例示的なプロセスは図9にも示される。
いくつかの例示的な実施形態では、本方法は、機能的効果優先度に応じてCHMにスコアを付与すること、及びサンプル当たりの遺伝子ごとに最も高い機能的効果優先度スコアを有するCHMを選択し、それによって医学的に重要な変異の集合を取得することをさらに含む。これらのステップは、SIFT(Loh et al.(2016);Nat Genet 48,1443-1448)(損傷を与える(damaging))、PolyPhen2 HDIV(損傷を与える(damaging)、及び損傷を与えるおそれがある(possibly damaging))、PolyPhen2 HVAR(損傷を与える(damaging)、及び損傷を与えるおそれがある(possibly damaging))、LRT(有害(deleterious))、及びMutationTaster(Schwarz et al.(2014);Nat.Methods 11、361-362)(必然的に疾患を引き起こす(disease causing automatic)、及び疾患を引き起こす(disease causing))を含むが、これらに限定されない任意好適な手段によって行い得る。
いくつかの例示的な実施形態では、低クオリティサンプルのD-statは、KS検定を使用して、サンプルの実際のアレルバランスの分布を、予想されるアレルバランスの分布と比較することによって判定される。
また、本開示は、集団におけるde novo変異(DNM)を同定する方法をも提供する。DNMを同定する方法の例を示すフローチャートが図3に与えられる。
本方法は、任意手段によって得られる任意の種類のヒト被験者に由来するあらゆる種類のDNA配列サンプルに適用し得る。バリアントの非限定的な例には、点突然変異、挿入、欠失、逆位、重複、及び多量体化が含まれる。ヒト被験者の種類の非限定的な例には、単一医療ネットワーク集団、複数医療ネットワーク集団、人種的、文化的、もしくは社会的に同種もしくは異種の集団、混合年齢集団もしくは年齢に関して均質な集団、地理的に集中した集団もしくは地理的に分散した集団、またはこれらの組み合わせからのヒト被験者が含まれる。DNA配列サンプルは、Dewey et al.(2016),Science 354,aaf6814-1~aaf6814-10に開示されているものを含むが、これに限定されない多くの方法のいずれかで取得し得る。
DNA配列サンプルは、エクソーム配列を含むか、またはエクソーム配列である。エクソームDNAは、一般に使用される方法のいずれかによって、またはDewey et al.(2016),Science 354,aaf6814-1~aaf6814-10に記載の通りに単離し得る。
複数のヒト被験者由来のDNA配列サンプルの中のバリアントを、任意好適な手段によって、ステップ22で同定し得る。バリアントを同定し得る手段の非限定的な例は、次のステップを含む。
- シークエンシングの完了時に、各シークエンシングランからの生データがローカルバッファストレージに収集され、自動解析のために、DNAnexusプラットフォーム(Reid et al.(2014);BMC Bioinformatics 15,30)にアップロードされる。
- サンプルレベルのリードファイルが、CASAVA(Illumina Inc.,San Diego,CA)で生成され、BWA-memでGRCh38(Li及びDurbin(2009);Bioinformatics 25,1754-176;Li(2013);arXiv q-bio.GN)にアライメントされる。
- 結果として得られたBAMファイルが、ソート、重複のマーク付け、及び推定インデル周辺のリードの局所的なアライメント修正の実行のために、GATK(McKenna et al.(2010);Genome Res.20,1297-1303)及びPicardを使用して処理される。
- シークエンシングされたバリアントは、転写産物及び遺伝子に対する機能的影響を明らかにするために、Ensembl85遺伝子定義を使用して、snpEFF(Cingolani et al.(2012);Fly(Austin)6,80-92)でアノテーションを付けられる。
本開示は、上記のどのステップにも限定されず、配列バリアントの取得を、任意好適な手段によって行い得ると解される。
図3は、集団におけるde novo変異(DNM)が同定される例示的な実施形態のフローチャートである。同定されたバリアントに基づく被験者の祖先スーパークラスの指定を、任意好適な手段によって、ステップ23で確立し得る。祖先スーパークラス内での被験者の第1の同祖性推定を、任意好適な手段によって、ステップ24で生成し得る。被験者の祖先スーパークラスから独立した被験者の第2の同祖性推定を、任意好適な手段によって、ステップ25で生成し得る。ステップ26で、第2の同祖性推定の1つ以上に基づいて、任意好適な手段によって、被験者を第1の第一度家族ネットワークにクラスタリングし得る。第1の第一度家族ネットワーク内での被験者の第3の同祖性推定を、任意好適な手段によって、ステップ27で生成し得る。第1の同祖性推定と第3の同祖性推定とを、任意好適な手段によって、ステップ28で統合して、統合同祖性推定を取得し得る。そのような手段の非限定的な例としては、PLINK(Chang et al.(2015);Gigascience 4,7)、及び本実施例に開示されているものが挙げられる。同祖性推定は、サンプルペアの間のゲノム全域にわたる0、1、及び2のIBD値の計算を含み得る。
さらに、統合同祖性推定に基づいて、ステップ29で核家族を構築し得る。核家族におけるバリアントを、ステップ30で同定し得る。ステップ31で、構築された核家族のトリオの両親及び子に由来するサンプルのバリアントに、遺伝子型尤度スコアを付与し得て、そのバリアントがde novo変異である確率を計算し得るとともに、独立して、トリオのどちらの親サンプルにも存在しない子サンプルのバリアントを単純に同定し得て、そのバリアントがde novo変異である確率を計算し得、次いで推定de novo変異の両方のセットを組み合わせ得、それによって、推定de novo変異のデータセットを形成する。上記のステップを実施する手段の非限定的な例としては、本実施例に開示されているものが挙げられる。
いくつかの例示的な実施形態では、本方法は、被験者の祖先スーパークラスの指定が確立される前に、同定されたバリアントをフィルタ処理することをさらに含み、いくつかの例示的な実施形態では、本方法は、被験者の第2の同祖性推定が生成される前に、同定されたバリアントをフィルタ処理することをさらに含む。バリアントを、任意好適な手段によってフィルタ処理してもよい。そのような手段の非限定的な例としては、PLINK(Chang et al.(2015);Gigascience 4,7)、及び本実施例に開示されているものが挙げられる。
バリアントをフィルタ処理することは、複数のヒト被験者由来のサンプル全体にわたって約10%を超える代替アレル頻度を有するバリアントを除去すること、またはp値が約10-6超である、ハーディ・ワインベルグ平衡(HWE)に反するバリアントを除去すること、または複数のヒト被験者由来のサンプルの約5%超でコールが欠損するバリアントを除去すること、またはこれらの組み合わせを含み得る。ある一定の選択基準を満たさないバリアントは除去され得、残りのバリアントが潜在複合ヘテロ接合性変異として選択され、潜在複合ヘテロ接合性変異は、本実施例に記載されたものを含む任意好適な手段によってフェージングされる。
いくつかの例示的な実施形態では、本方法は、同定されたバリアントがフィルタ処理された後に、低クオリティサンプルを除去することを含む。低クオリティサンプルを、任意好適な手段によって除去し得る。そのような手段の非限定的な例としては、一般的に知られていて、したがって本明細書ではこれ以上詳述されないDewey et al.(2016),Science 354,aaf6814-1~aaf6814-10に開示されているもの、及び本実施例に開示されているものが挙げられる。
いくつかの例示的な実施形態では、D-statが0.12を超えるか、もしくは20xのリードカバレッジが75%未満であるか、またはその両方であるサンプルが、除去される低クオリティサンプルであるように、パラメータが調節される。いくつかの例示的な実施形態では、低クオリティサンプルのD-statは、KS検定を使用して、サンプルの実際のアレルバランスの分布を、予想されるアレルバランスの分布と比較することによって判定される。
第1の同祖性推定と第3の同祖性推定とを統合することは、第3の同祖性推定に特有のペアワイズ同祖性推定を補って、第1の同祖性推定を増強することを含み得、これは、例えば限定されないが、PLINK(Changetal.(2015);Gigascience4,7)、及び本実施例に開示される手段によって促進され得る。
バリアントをフィルタ処理することは、子サンプルに約35未満の遺伝子型クオリティ(GQ)アノテーションを有するバリアント、またはサンプル全体にわたって10以上の代替アレルカウント(AC)を有するバリアント、または子サンプルに、約7未満のリードデプス(DP)と約4未満の代替DPとを有するバリアント、またはどちらか一方の親サンプルに約2%を超えるアレルバランス(AB)を有するバリアント、または子サンプルに約15%未満のアレルバランス(AB)を有するバリアント、または子サンプルに約90%を超えるABを有するバリアント、またはどちらか一方の親サンプルに代替アレルホモ接合性を有するバリアント、またはこれらの組み合わせのバリアントを除去することを含み得る。推定de novo変異を中信頼度de novo変異として分類することは、推定de novo変異のアレルバランス(AB)が、子サンプルで約15%以上、かつ各親サンプルで約2%以下の場合、及び推定de novo変異のマッピングクオリティ(MQ)が約40未満ではない場合、及び推定de novo変異のクオリティバイデプス(quality by depth)(QD)の値が約2未満ではない場合、及び推定de novo変異のマイナーアレルカウント(MAC)がサンプル全体にわたって約20未満の場合、及び推定de novo変異のソフトクリップリードが、推定de novo変異のキャリアにおけるバリアントサイトで約3以下の場合、及び推定de novo変異が、約4を超えるモノポリマーランを持つINDELではない場合に行われ得る。中信頼度de novo変異を高信頼度de novo変異として分類することは、中信頼度de novo変異が、親サンプルにおいて約90以上の遺伝子型クオリティ(GQ)アノテーションを有する場合、及び中信頼度de novo変異が、各親サンプルにおいて約10以上のリードデプス(DP)を有する場合、及び中信頼度de novo変異が、子サンプルにおいて約7以上の代替DPを有する場合、及び中信頼度de novo変異が、約3より大きいQDのSNPを有する場合、及び中信頼度de novo変異が、約5より大きいQDのINDELを有する場合に行われ得る。これらのステップは、任意好適な手段によって行い得る。
本実施例に記載されたものを含む。また、これらの例示的な実施形態は図7にも示される。
本明細書で使用される「D-stat」という用語は、低クオリティサンプルを同定するのに生成し使用し得るQCメトリクスを指す。サンプルのクオリティが低いのは汚染に起因する場合があり、下流の分析で問題を引き起こす場合がある。サンプルのD-statは、例えば、サンプルの実際のアレルバランスの分布を参照アレルバランスの分布(例えば、予想されるアレルバランスの分布)と比較することによって計算してもよい。この参照分布は、例えば、分析すべきサンプルをクエリするのに使用されるものと同じプラットフォームを使用してキャプチャされ配列決定される汚染の形跡が全くない複数のサンプルから計算してもよい。本明細書で使用されるD-statのQCメトリクス値は、p値の計算に先立ってK-S(Kolmogorov-Smirnov)検定から生成されるD統計値と同等である。D-statには単位がない。K-S検定のD統計値は、0と1との間の値になり、1は、参照分布及びサンプル分布の累積分布間の極大差を意味する。いくつかの例示的な実施形態では、サンプルの実際のアレルバランスの分布を、K-S検定に従って計算されたアレルバランスの予想される分布/参照分布と比較することによって、低クオリティサンプルが同定される。いくつかの例示的な実施形態では、特定のD-stat値を有すると判定されたサンプルは、低クオリティサンプルとみなされ、以降の分析から除外される。いくつかの例示的な実施形態では、低クオリティであるとみなされ、除外すべきサンプルのD-stat値は、0.01、0.02、0.03、0.04、0.05、0.06、0.07、0.08、0.09、0.1、0.11、または0.12よりも大きい。好ましい実施形態では、低クオリティであるとみなされ、除外すべきサンプルのD-stat値は、約0.12よりも大きい。さらに一層好ましい実施形態では、低クオリティであるとみなされ、除外すべきサンプルのD-stat値は、0.12よりも大きい。
記載されており、または例示されている方法のいずれも、非一時的なコンピュータ実行方法として、及び/またはシステムとして実施し得る。このために、当業者に知られている任意好適なコンピュータシステムを使用してよい。
図10は、本方法及びシステムが動作し得る例示的な環境201の様々な態様を示す。本方法は、デジタル及びアナログの両方の機器を使用する各種のネットワーク及びシステムで使用し得る。本明細書では、機能説明が提供され、それぞれの機能が、ソフトウェア、ハードウェア、またはソフトウェアとハードウェアとの組み合わせによって実行し得ることが示されている。
環境201は、ローカルデータ/処理センタ210を含み得る。ローカルデータ/処理センタ210は、1つ以上のコンピューティングデバイスの間の通信を助けるために、ローカルエリアネットワークなどの1つ以上のネットワークを含んでもよい。1つ以上のコンピューティングデバイスを用いて、生物学的データの保存、処理、分析、出力、及び/または視覚化を行い得る。環境201は、任意選択で、医療データプロバイダ220を含んでもよい。医療データプロバイダ220は、生物学的データの1つ以上の情報源を含み得る。例えば、医療データプロバイダ220は、1人以上の患者の医療情報が入手可能な1つ以上の医療システムを含み得る。医療情報は、例えば、病歴、医療専門家の見解及び注釈、検査室報告、診断、医師命令、処方箋、バイタルサイン、体液バランス、呼吸機能、血液パラメータ、心電図、X線写真、CTスキャン、MRIデータ、臨床試験結果、診断、予後、評価、入院及び退院の記録、ならびに患者登録情報を含み得る。医療データプロバイダ220は、1つ以上のコンピューティングデバイスの間の通信を促進するために、ローカルエリアネットワークなどの1つ以上のネットワークを含んでもよい。1つ以上のコンピューティングデバイスを用いて、医療情報の保存、処理、分析、出力、及び/または視覚化を行い得る。医療データプロバイダ220は、医療情報を匿名化し、この匿名化医療情報をローカルデータ/処理センタ210に提供し得る。匿名化医療情報は、医療情報を匿名化された状態に保ちながらも、ある患者の医療情報と別の患者の医療情報を区別するために、患者ごとに一意の識別子を含み得る。匿名化医療情報により、患者の個人情報が患者の特定の医療情報に関連付けされないようにする。ローカルデータ/処理センタ210は、匿名化医療情報を分析して、(例えば、国際疾病分類「ICD」コード及び/または医師診療行為用語「CPT」コードを付与することによって)各患者に1つ以上の表現型を付与してもよい。
環境201は、NGSシークエンシング施設230を含み得る。NGSシークエンシング設備230は、1つ以上のシークエンサ(例えば、Illumina HiSeq 2500、Pacific Biosciences PacBio RS IIなど)を含み得る。1つ以上のシークエンサは、エクソームシークエンシング、全エクソームシークエンシング、RNA-seq、全ゲノムシークエンシング、ターゲットシークエンシングなどのために構成され得る。例示的な態様では、医療データプロバイダ220は、匿名化医療情報に関連している患者由来の生体サンプルを提供してもよい。一意の識別子が用いられて、生体サンプルと、その生体サンプルに対応する匿名化医療情報との間の関連が維持され得る。NGSシークエンシング施設230は、生体サンプルを基に各患者のエクソームを配列決定し得る。シークエンシングに先立って生体サンプルを保管するために、NGSシークエンシング施設230は、バイオバンク(例えば、Liconic Instruments製)を備えてもよい。チューブ(患者に割り当てられた各チューブ)中に生体サンプルを収容することができ、各チューブは、サンプルをローカルデータ/処理センタ210に自動的に記録するためにスキャンされ得るバーコード(または他の識別子)を備えてもよい。NGSシークエンシング施設230は、均質のデータと効率的な無停止運用とを保証するために、シークエンシングの1つ以上のフェーズで用いるための1つ以上のロボットを備えてもよい。このようにして、NGSシークエンシング施設230は、年間数万のエクソームを配列決定することができる。一態様では、NGSシークエンシング施設230は、月当たり少なくとも1000、2000、3000、4000、5000、6000、7000、8000、9000、10,000、11,000、または12,000の全エクソームを配列決定する機能的能力を有する。
NGSシークエンシング施設230によって生成された生物学的データ(例えば、未加工のシークエンシングデータ)は、ローカルデータ/処理センタ210に転送され得、次いで、ローカルデータ/処理センタ210は、生物学的データをリモートデータ/処理センタ240に転送してもよい。リモートデータ/処理センタ240は、1つ以上のコンピューティングデバイスを含むクラウドベースのデータストレージ及び処理センタを備え得る。ローカルデータ/処理センタ210及びNGSシークエンシング施設230は、1つ以上の大容量ファイバ線を介して、リモートデータ/処理センタ240との間でデータ通信をすることができるが、他のデータ通信システム(例えば、インターネット)も企図されている。例示的な態様では、リモートデータ/処理センタ240は、例えばアマゾンウェブサービス(DNAnexus)などのサードパーティシステムを含んでもよい。リモートデータ/処理センタ240は、分析ステップの自動化を容易にすることができ、安全な方法で1人以上の協働者250とデータを共有することを可能にする。ローカルデータ/処理センタ210から生物学的データを受け取ると、リモートデータ/処理センタ240は、バイオインフォマティクスツールを使用して、一次及び二次のデータ分析のための自動化された一連のパイプラインステップを実行し、結果としてサンプルごとにアノテーション付きバリアントファイルを生成し得る。そのようなデータ分析(例えば、遺伝子型)の結果は、ローカルデータ/処理センタ210に返送され得て、例えば、ラボラトリ情報管理システム(LIMS)に組み込まれ得、各生体サンプルの状態を維持するように構成され得る。
その後、ローカルデータ/処理センタ210は、NGSシークエンシング施設230及びリモートデータ/処理センタ240によって得られた生物学的データ(例えば、遺伝子型)を、匿名化医療情報(同定された表現型を含む)と組み合わせて利用して、遺伝子型と表現型との間の関連を同定してもよい。例えば、ローカルデータ/処理センタ210は、表現型優先のアプローチを適用してもよく、このアプローチでは、特定の疾患領域、例えば心臓血管疾患に対する極端な血中脂質において治療可能性を有し得る表現型が定義される。もう一つの例は、典型的な併存疾患の範囲から保護されていると思われる個人を同定するための肥満患者の研究である。別のアプローチは、遺伝子型と、例えば、遺伝子Xが疾患Yの発症か、または疾患Yからの保護に関与しているという仮説とに端を発するものである。
例示的な態様では、1人または複数人の協働者250が、インターネット260などのネットワークを介して、生物学的データ及び/または匿名化医療情報の一部または全てにアクセスし得る。
図11に示す例示的な態様では、ローカルデータ/処理センタ210及び/またはリモートデータ/処理センタ240の1つ以上は、遺伝データ構成要素300、表現型データ構成要素310、遺伝的バリアント-表現型関連データ構成要素320、及び/またはデータ分析構成要素330のうちの1つ以上を含む1つ以上のコンピューティングデバイスを備え得る。遺伝データ構成要素300、表現型データ構成要素310、及び/または遺伝的バリアント-表現型関連データ構成要素320は、配列データの品質評価、参照ゲノムへのリードアライメント、バリアントの同定、バリアントのアノテーション、表現型の同定、バリアント-表現型関連の同定、データの可視化、それらの組み合わせなどのうちの1つ以上のために構成され得る。
例示的な態様では、構成要素の1つ以上は、全面的にハードウェアの実施形態、全面的にソフトウェアの実施形態、またはソフトウェア態様とハードウェア態様とを組み合わせた実施形態の形を取り得る。さらに、方法及びシステムは、記憶媒体中に取り入れられたコンピュータ可読プログラム命令(例えば、非一時的なコンピュータソフトウェア)を有するコンピュータ可読記憶媒体上のコンピュータプログラム製品の形を取ってもよい。より具体的に、本方法及びシステムは、ウェブ上で実行されるコンピュータソフトウェアの形を取ってもよい。ハードディスク、CD-ROM、光学記憶装置、または磁気記憶装置を含む任意好適なコンピュータ可読記憶媒体を利用してよい。
例示的な態様では、遺伝データ構成要素300は、1つ以上の遺伝的バリアントに機能的にアノテーションするように構成され得る。また、遺伝データ構成要素300は、1つ以上の遺伝的バリアントの保管、分析、受領などを行うように構成される場合もある。1つ以上の遺伝的バリアントは、1人以上の患者(被験者)から取得された配列データ(例えば、未加工の配列データ)から判断してアノテーションされ得る。例えば、1つ以上の遺伝的バリアントは、少なくとも100,000、200,000、300,000、400,000、または500,000人の被験者のそれぞれからアノテーションされ得る。1つ以上の遺伝的バリアントを機能的にアノテーションした結果、遺伝的バリアントデータが生成される。例として、遺伝的バリアントデータは、1つ以上のバリアントコールフォーマット(VCF)ファイルを含む場合がある。VCFファイルは、SNP、インデル、及び/または構造バリエーションのコールを構成するためのテキストファイル形式である。転写産物/遺伝子に対する機能的影響についてバリアントが評価され、潜在的な機能喪失(pLoF)候補が同定される。Ensembl75遺伝子定義を使用して、バリアントがsnpEffでアノテーションを付けられ、次いで機能的アノテーションが、バリアント(及び遺伝子)ごとにさらに処理される。
本明細書で提供される方法ステップの数字及び/または文字による連続したラベル付けは、方法またはそのいかなる例示的な実施形態をも特定の示された順序に限定することを意味しない。
特許、特許出願、公開された特許出願、アクセッション番号、技術論文、学術論文など、種々の出版物が本明細書全体を通して引用されている。これらの引用文献のそれぞれは、参照によりその全体があらゆる目的で本文書に援用される。
本開示は、本開示をより詳細に説明するために提供される以下の実施例を参照することによって、より完全に理解されるであろう。これらは例示を意図したものであり、本開示の範囲を限定するものと解釈されるべきではない。
実施例1.1
61Kのヒトエクソームのコホートにおける関係の推定と近縁性の内容
61Kのヒトエクソームのコホートを分析した。このコホートは、2014年に開始されたRegeneron Genetics Center(RGC)及びGeisinger Health System(GHS)による研究を起源としている(Dewey et al.(2016),Science 354,aaf6814-aaf6814)。このDiscovEHR研究では、移住率の低い集団を扱う単一のヘルスケアシステムに属する患者を密にサンプリングした。この61Kのヒトエクソームコホートは、本明細書ではDiscovEHRデータセットと呼ぶ。膨大な量の家族構造をDiscovEHRデータセット内で同定し、本明細書に開示されたシミュレーションでは、本研究で250Kのターゲットを確認したときに、データセット内の個人の70%~80%が、第一度または第二度の血縁者を持つであろうと見積もられた。
同祖性(IBD)推定を用いて、データセット内の様々な種類の家族関係を同定し、PRIMUS(Staples et al.(2014),Am.J.Hum.Genet.95,553-564)を用いて、ペアワイズ関係を種々の家族クラスに分類するとともに、家系を再構築した(実施例8でさらに説明する)。全エクソームシークエンシング(WES)データから遠縁の血縁者のIBD比率を正確に推定することに制限があるため、DiscovEHRデータセットサンプルの間に推定した第一度、第二度、及び信頼度の高い第三度の関係のみを含めた。
データセット内に、合計で、20の一卵性双生児、8,802の親子関係、6,122の完全同胞関係、及び約20,000の第二度関係を同定した(図12A)。第二度関係及び第三度関係のIBD共有分布は互いに重なり合うため、この研究では、2つの期待平均値の間にハードカットオフ(hard cutoff)を設定した。第三度関係(図12Aではアスタリスクで印が付けられている)は、エクソームデータの技術的な制限、ならびにより遠い関係クラス(例えば、第四度及び第五度)の期待される平均IBD比率の周囲の拡張変動及び重複変動のために、正確に推定することは困難である。次に、無向グラフを生成するために、個人をノードとし、関係をエッジとして扱った。第一度関係のみを用いて、第一度家族ネットワークと呼ばれる7,684の連結成分を同定した。図12Bは、第一度家族ネットワークのサイズの分布を示しており、これは2から25までのシークエンシングされた個人の範囲に広がる。同様に、7,136の第二度家族ネットワークが見つかり、最大のものは7,123の個人を含んでいた(データセット全体の約12%、図12D)。図12Dでは、第二度家族ネットワーク内の第一度家族ネットワークは、ネットワーク(図12Cに示す第一度家族ネットワークの家系を含む)内の個人の数に比例したサイズの赤いボックスとして表す。単一の個人は、青いエッジとして描いた第二度の関係で結ばれる黒いノードとして表す。
第二度家族ネットワーク内に、約4,500の第三度の関係を同定することもできた。表示される祖先グループ内のIBD推定のための最小IBDカットオフを緩和することにより、DiscovEHRデータセット内で50Kを大幅に超す第三度関係を同定した。ヨーロッパ系の個人はDiscovEHRデータセットの96.5%を構成するに過ぎないが(以下の表1a参照)、データセットに見つかるペアワイズ関係の圧倒的多数(>99%)はヨーロッパ系の個人を含む(以下の表1b参照)。いずれにしても、ヨーロッパ系の人々の間、非ヨーロッパ系の人々の間、及び種々の祖先を持つ個人の間には、多くの関係が存在する。例えば、子供が厳密に参照集団に適合しなかったために、子供に不明の祖先が割り当てられたヨーロッパ人の父親と東アジア人の母親とを含むトリオを、DiscovEHRデータセット中に見出した。
Figure 0007277438000001
Figure 0007277438000002
Figure 0007277438000003
血縁者の累増率は、実験的にサンプルを確認した割合を遥かに超え(図13A)、シミュレーションを通じて(図14A)、血縁者の累増率は、サンプルを確認する割合を遥かに超えていた。さらなるペアワイズ関係の累増は、より多くの個人がこれらの関係に関わる結果となった。現在、61Kの個人の50.4%は、DiscovEHRデータセットに第一度または第二度の血縁者を1人以上有する(図13B)。
実施例1.2
92Kのヒトエクソームのコホートにおける関係の推定と近縁性の内容
より大きな92,455のヒトエクソームの臨床コホートを分析した。このコホートは、2014年に開始されたRegeneron Genetics Center(RGC)及びGeisinger Health System(GHS)による進行中の研究を起源としている(Staples et al.(2018),Am.J.Hum.Genet.102(5):874-889)。同様に、この拡張されたDiscoverEHRコホートは、中央ペンシルベニアの移住率が低い大部分が農村の集団を扱う単一のヘルスケアシステムからの参加者を密にサンプリングしたものである。
準備しシークエンシングした第1の61Kサンプル(実施例1.1)を含むセットを、「VCRomeセット」と称した。残りの31Kサンプルのセットは、NimbleGenプローブによるキャプチャの代わりに、わずかに改変したバージョンのIDTのxGenプローブを用い、NimbleGenのVCRomeキャプチャ試薬によってカバーされるが、標準のxGenプローブによってはカバーされにくいゲノムの領域をキャプチャするために、補足プローブを用いたことを除いて、同じプロセスによって調製した。キャプチャした断片はストレプトアビジン結合ビーズに結合させ、非特異的DNA断片を、製造元(IDT)の推奨プロトコルに従って、一連のストリンジェント洗浄によって除去した。この第2のサンプルセットを、「xGenセット」と称した。バリアントコールを、GATKで生成した。GATKは、推定インデル周辺の各サンプルの、アライメントし、重複にマークをつけたリードを、局所的にアライメント修正するために用いた。サンプルがゲノムバリアントコールフォーマット(gVCF)のゲノム参照と異なる全てのエクソン位置を同定するために、GATKのHaplotypeCallerを使用して、INDELのアライメント修正した重複マーク付きのリードを処理した。各サンプル上のGATKのGenotypeGVCFと、単一ヌクレオチドバリアント(SNV)及びインデルの両方を参照と比較して同定する単一サンプルのバリアントコールフォーマット(VCF)ファイルを出力するランダムに選ばれた50個のサンプルのトレーニングセットとを用いて、遺伝子型決定を達成した。次に、単一サンプルのVCFファイルを使用して、両方のセットの単一サンプルVCFファイルからの全ての変異性サイトを含む疑似サンプルを作成した。さらに、200の単一サンプルgVCFファイルを疑似サンプルと共に共同でコールすることによって、VCRomeセット用の独立したpVCFファイルは、2つのキャプチャセットにわたる全ての変異性サイトで、各サンプルのコールまたは無コールを強制する。200サンプルのpVCFファイル全てを組み合わせて、VCRomeのpVCFファイルを作成し、このプロセスを繰り返してxGenのpVCFファイルを作成した。VCRomeのpVCFファイルとxGenのpVCFファイルとを組み合わせて、結合pVCFファイルを作成した。GRCh38に配列リードをアライメントし、Ensembl85遺伝子定義を用いることによって、バリアントにアノテーションを付けた。遺伝子定義は、アノテーション付きの開始及び停止を伴うタンパク質コードである19,467の遺伝子に対応する54,214の転写産物に制限した。サンプルのQCプロセスの後に、92,455のエクソームを分析のために残した。
92,455の個人の拡張したDiscovEHRデータセットから、43の一卵性双生児、16,476の親子関係、10,479の完全同胞関係、及び39,000の第二度関係を同定した(図16、パネルA)。無向グラフを生成するために、個人をノードとし、関係をエッジとして扱った。第一度関係のみを用いて、第一度家族ネットワークと呼ばれる12,594の連結成分を同定した。図16、パネルBは、第一度家族ネットワークのサイズの分布を示しており、これは2から25までのシークエンシングされた個人の範囲に広がる。同様に、最大のもので19,968(データセット全体の22%、図16、パネルC)の個人を含む、10,173の第二度家族ネットワークを同定した。また、第二度家族ネットワーク内に、約5,300の第三度の関係も同定した。第二度家族ネットワークを考慮することなく、祖先グループ内のIBD推定に低いIBDカットオフ(PIJ_HAT>0.09875)を使用することで、拡張したDiscovEHRコホート内に、100,000を超える第三度の関係を同定した。拡張したDiscovEHR個人の95.9%がヨーロッパ系であったことを考えると(表2a)、見つかったペアワイズ関係の大多数(98.6%)がヨーロッパ系の2つの個人間であったことは驚くに値しない(表2b)。それにもかかわらず、ヨーロッパ系の人々の間、同じ非ヨーロッパ系の祖先を持つ人々の間、及び異なる祖先を持つ個人の間に、多くの関係を同定した。例えば、ヨーロッパ人の親1名、東アジア人の親1名、ならびにゲノムの種類が混在しているために祖先が超集団に割り当てられていない子供を含むトリオがいくつかあった。重要なことには、実験的に(図17A)かつシミュレーションを通じて(図18A)、血縁者の累増率は、サンプルを確認する割合を遥かに超えることが判明した。このことは、データセット内で可能なペアワイズ関係の数が、データセットのサイズが大きくなるにつれて組み合わせ的に増加し、また、データセット内の以前は無関係であった個人が、新しく同定された関係に関与してくる可能性も増加することを考えると、予想されたことである。現在、拡張したDiscovEHRコホートのうちの39%の個人は、データセット内に少なくとも1人の第一度血縁者を有する可能性があり、参加者の56%は、データセット内に1人以上の第一度または第二度の血縁者を有する可能性がある(図17、パネルB)。
Figure 0007277438000004
Figure 0007277438000005
Figure 0007277438000006
Figure 0007277438000007
実施例2
SimProgenyを用いたシミュレーションと近縁性予想
DiscovEHR及び拡張DiscovEHRのデータセット内の関係ネットワークの成長のモデル化、理解、及び予測をしようとする試みで、シミュレーションフレームワーク(以下「SimProgeny」)を開発した。これは、複数のサブ集団にまたがって分散した数百年にわたる数百万人の血統をシミュレートすることが可能である。これらのシミュレートされた集団から、SimProgenyは、様々なサンプリングアプローチをモデル化し、研究者が、所与の集団及びサンプリングパラメータのセットに対して見出すことを期待するだろう近縁性の量を推定することができる(実施例17参照)。
SimProgenyを用いて、DiscovEHR及び拡張DiscovEHRの集団と、それぞれ、それらからの第1の61K参加者及び第1の92K参加者の確認とをシミュレートした。シミュレーションは、DiscovEHR及び拡張DiscovEHRの参加者を、集団からランダムにサンプリングしたのではなく、近親血縁者に対してデータセットを濃縮したことを示す。図14A及び図14Bに示すように、薄い線で結ばれた句読点で表示される周期的な「フリーズ」で、実データを計算した。また、61K人のフリーズにおいて同定したサンプル及び関係を取得し、次いで、確認順序をシャッフルして、61KのDiscovEHR参加者の前半は、後半に比べて、第一度の関係が濃縮されたことを示した。DiscovEHRを確認した実際の集団と同様のパラメータを使用して、様々なサイズの集団をシミュレートした。次に、これらの集団のそれぞれからランダムに確認を行って、どの集団のサイズが実データに最もよく適合しているかを確認した。重要な留意点は、これらの集団のサイズはどれも実データに適合せず、ランダムな確認アプローチは適合性が低いということである。第一度血縁者を濃縮する別の確認アプローチは、ランダムな確認と比べて、より良好な適合をもたらし得る。図14Aは、サイズ270Kの有効サンプリング集団における第一度血縁者ペアの確認が、実データのシャッフルしたバージョンによく適合するが、確認した61Kの参加者より少ない血縁者ペアの数を低く見積もり過ぎ、61Kの参加者を超える血縁者ペアの数を著しく高く見積もり過ぎることを示す。図14Bは、270Kの集団が、1人以上の第一度血縁者を持つ個人の数については、シャッフルされた実データに最もよく適合するが、実データにはあまり適合しないことを示す。
拡張したDiscovEHRデータセットを用いて、同様の結果が観察された(図17A及び図17B)。次に、92K人のフリーズにおいて同定したサンプル及び関係をシャッフルして、92Kの拡張DiscovEHR参加者の前半は、後半に比べて、第一度の関係が濃縮されたことを示した。次に、これらの集団のそれぞれからランダムに確認を行って、どの集団のサイズが実データに最もよく適合しているかを確認した。図17Aは、サイズ403Kの有効サンプリング集団における第一度血縁者ペアの確認が、実データのシャッフルしたバージョンによく適合するが、確認した92Kの参加者より少ない血縁者ペアの数を低く見積もり過ぎ、92Kの参加者を超える血縁者ペアの数を著しく高く見積もり過ぎることを示す。図17Bは、403Kの集団が、1人以上の第一度血縁者を持つ個人の数については、シャッフルされた実データに最もよく適合するが、実データにはあまり適合しないことを示し、拡張されたDiscovEHRの参加者がランダムに確認されなかったことを示唆する。
近親血縁者の濃縮は、DiscovEHR(図15A及び図15B)及び拡張DiscovEHR(図18A及び図18B)の実データに、より良好に適合するシミュレーションを生成するクラスタリングされた確認アプローチ(実施例17参照)を使用することによってモデル化した。図15及び図18では、共に、薄い線で結ばれた句読点で表示される周期的な「フリーズ」で、実データを計算した。ほとんどのシミュレーションパラメータは、実際の集団の人口統計及びDiscovEHR確認アプローチに関する情報に基づいて設定した。ただし、2つのパラメータが不明であり、実データへの適合を踏まえて選択した。すなわち、1)サンプルが確認された有効な集団サイズ、及び2)「クラスタリングされた確認」と呼ばれる、以前に確認された第一度血縁者を考慮して誰かが確認される可能性の増加である。全てのパネルは、推定される有効な集団サイズにまたがる、シミュレーションされた同じ3つの集団のサイズを示す。クラスタリングされた確認は、ポアソン分布の第一度血縁者の乱数と連動して、個人をランダムに確認することによってシミュレートした(ポアソン分布のラムダ値は凡例に示される)。これらのシミュレーション結果は、有効なサンプリング集団のサイズが約475Kの個人であり、ラムダ0.2のポアソン分布が第一度血縁者の濃縮と最もよく一致することを示唆した。このことは、現在の参加者の大多数が、GHS管轄区域全体(>250万人)に均等に分散するのではなく、この例ではPA州ダンビル(約500K人)などの特定の地域の地理的領域に住んでいるという理解と一致した。
実データに適度に適合するシミュレーションパラメータを同定した後、SimProgenyを使用して、DiscovEHR及び拡張DiscovEHRの研究が250Kの参加者という目標まで拡張したときに期待されるべき第一度関係の総計予測を得た。結果は、参加者の確認を同じようにして継続するならば、DiscovEHR(図15C)及び拡張DiscovEHR(図18C)では、DiscovEHR参加者の約60%が関与して(図15D)、かつ拡張DiscovEHR参加者の約60%が関与して(図18D)、約150Kの第一度関係が期待されるはずであることを示した。
次に、シミュレーション分析を、第二度関係を含めるように拡張し、シミュレーション結果から、250Kの参加者では、DiscovEHR(図19)及び拡張DiscovEHR(図20)の個人の70%超が関与する第一度及び第二度の関係を合わせた200Kを優に超えることが期待されることが示唆された。この分析では、図中の薄い線で結ばれた句読点で表示される周期的な「フリーズ」で、実データを計算した。ほとんどのシミュレーションパラメータは、実際の集団の人口統計及びDiscovEHR確認アプローチに関する情報に基づいて設定した。全てのパネルは、推定される有効な集団サイズにまたがる、シミュレーションされた同じ3つの集団のサイズを示す。クラスタリングされた確認は、ポアソン分布の第一度血縁者の乱数及び第二度血縁者の別個の乱数と連動して、個人をランダムに確認することによってシミュレートした(どちらのポアソン分布も、図の凡例に示されているラムダを有する)。
シミュレーション結果は、DiscovEHRのHPG研究において近縁性の明確な濃縮を示し、さらなる参加者の確認を続けたときに見られることが予想される膨大な量の近縁性に対する重要な洞察を提供した。
実施例3.1
近縁性をDicoverEHRデータセットの確率変数のように扱うことの代わりに、近縁性を活用すること
家系再構築ツールPRIMUS(Staples et al.(2014),Am.J.Hum.Genet.95,553-564.)を使用して、DiscovEHRデータセットにおいて、7,684の全ての第一度家族ネットワークを再構築し、IBD推定と報告された年齢とを考慮するとき、これらの家系の98.9%が一意に再構築されることが分かった。これらの家系には、1,081の核家族(925のトリオ、134のカルテット、19のクインテット、及び3つのセクステット)が含まれていた。以下の表3は、トリオの祖先別の内訳を示す。1,081の核家族を、それらの個々のトリオ構成要素に分割した。例えば、カルテットは、同じ親を持つ2つの別々のトリオに分けることになる。DiscovEHRコホートは、大半はヨーロッパ人であったため、トリオの大多数は、ヨーロッパ系の個人を含んでいた。不明の祖先を持つ個人は、概して、異なる系統を持つ親の子供であった。例えば、東アジア人-ヨーロッパ人-不明のトリオの3つ全ては、ヨーロッパ人の父親と東アジア人の母親とを含んでおり、混合された子供をもたらす。これらのヨーロッパ人-東アジア人の混合した個人に厳密に一致する参照集団がなかったため、これらは不明の祖先という結果となった。
Figure 0007277438000008
図12Cは、25のシークエンシングされた個人を含むDiscovEHRデータセットにおいて同定された最大の第一度家系を示す。これらの関係及び家系は、以下を含むいくつかの方法で使用した。
複合ヘテロ接合性変異
人類遺伝学の主な目的は、ヒトゲノム中の全ての遺伝子の機能をより深く理解することである。ホモ接合性機能喪失型変異(LoF)は、これらの「ヒトノックアウト」(KO)の表現効果を分析することにより、遺伝子機能を洞察する力を得るための強力なツールである。最近の大規模シークエンシング研究では、稀少(MAF<1%)ホモ接合性LoFが特に注目されており、多くの遺伝子-表現型相互作用の同定において重要になっている(Lek et al.(2016),Nature Publishing Group 536,285-291、Dewey et al.(2016),Science 354,aaf6814-aaf6814、Saleheen et al.(2017),Nature Publishing Group 544,235-239、及びNarasimhan et al.(2016),Science 352,474-477)。2つのヘテロ接合性LoFの稀少複合ヘテロ接合性変異(CHM)は、機能的には稀少ホモ接合性KOと同等であるが、これらは、このような大規模シークエンシング研究ではめったに調査されない(Lek et al.(2016),Nature Publishing Group 536,285-291、Dewey et al.(2016),Science 354,aaf6814-aaf6814、及びSaleheen et al.(2017),Nature Publishing Group 544,235-239)。LoFの稀少CHMの正確な同定は、(1)稀少CHMがヒト遺伝子KO数を実質的に増加させ、統計的検出力を改善するため、(2)稀少CHMのKOが、ホモ接合性キャリアを欠き得る極めて稀少なヘテロ接合性変異を含む可能性があるため、及び(3)稀少CHMが、「ヒトKOプロジェクト」のためのより完全なKOセットを提供するため、価値がある(Saleheen et al.(2017),Nature Publishing Group 544,235-239、Perdigoto(2017),Nat.Rev.Genet.18,328-329)。
DiscovEHRデータセット内の稀少CHMの調査を行った。第一に、有害であるという強力な証拠を持つ、推定LoF(pLoF、すなわち、ナンセンス変異、フレームシフト変異、またはスプライスサイト変異)またはミスセンスバリアントのいずれかである稀少ヘテロ接合性バリアントのペアで構成された、39,459の高クオリティの潜在CHM(pCHM)を同定した(実施例10参照)。第二に、pCHMは、EAGLEを使用するアレル頻度ベースのフェージングと、再構築された家系及び関係のデータを使用する家系ベースのフェージングとの組み合わせを使用してフェージングした(図9)。EAGLEは、トリオ検証に基づいて、91.4%の精度でpCHMをフェージングした(以下の表4参照)。ただし、このコホート内には広範囲にわたる家系及び関係のデータがあったため、pCHMのほぼ3分の1を、これらのデータに基づいて約100%の精度でフェージングすることができ(以下の表4を参照)、不正確なフェージングを約31%削減することが可能であった。フェージングされたpCHMは、シングルトンから1%MAFまでの全範囲に及んだ(以下の表5参照)。
Figure 0007277438000009
トリオの子に発生したMAF<1%かつMAC>1の全てのpCHMを、再構築したトリオを用いてフェージングし、「真」であると見なした。子において、寄与バリアントの1つ以上がde novoであると判定されたpCHMは除外した。次に、フェージング済みトリオのpCHMを使用して、他のフェージング方法を評価した。EAGLEの精度は、再構築した各核家族から1人の子供の第一度血縁者を全て除去し、次いで残りのデータセット中の全てのバリアントをフェージングすることによって評価した。EAGLEのフェージング済みpCHMを、トリオのフェージング済みpCHMと比較した。
Figure 0007277438000010
pCHMの精度は極めて稀少なバリアントでは低下する傾向があったため、2つのpCHMバリアントのうちのより稀少なMAFを用いて、pCHMをそれぞれの頻度ビンに入れた。MACが1のpCHMは、関係データを用いてフェージングし、pCHMキャリア中のde novo変異ではないと見なした。pCHMの不明な相は、片方または両方のpCHMバリアントがEAGLE(MAC=1または欠損>10%)によってフィルタ処理された結果であり、フェージング用の関係データがないために発生した。
処理後、pCHMの39%をトランスにてフェージングし、61Kの個人の11,375の間に分布する、13,335の稀少で有害なCHMの高信頼度セットを産出した(平均=0.22;最大=6;図21、パネルA)。シスでのpCHMバリアント間のゲノム距離の中央値(5,308bps)は、トランスでのバリアント間の中央値距離(11,201bps;図22)の半分より少し少なかった。CHMのほぼ3分の1は少なくとも1つのpLoFを含み、CHMの9.8%は2つのpLoFバリアントで構成されていた(以下の表6を参照)。19,467の標的化遺伝子のうち3,385を超える遺伝子が1つ以上のCHMキャリアを含み(以下の表7を参照)、そのうちの1,555(46%)は複数のキャリアを有していた。85を超えるCHMキャリアを持つ11の遺伝子は、ExACのpLIスコアに基づいて、ゲノムで最もLoF耐性があると推定された(Lek et al.(2016),Nature Publishing Group 536,285-291)(以下の表8を参照)。
Figure 0007277438000011
この表は、稀少(<1%MAF)pLoFとミスセンスバリアントとで構成されているCHMの内訳を提供する。また、これらのCHMのいくつがインデル-インデル、インデル-SNP、及びSNP-SNPペアリングを構成しているかも示す。
Figure 0007277438000012
Figure 0007277438000013
CHMが最も多い11遺伝子については、ExACによって報告されたそれらのpLIスコアを表示している。また、全ての遺伝子をpLIスコアでランク付けし、報告pLiスコアを持つ遺伝子の総数で割ることによって計算した、LoF許容値に対する各遺伝子のパーセンタイルも表示している。
SSPOのpLIスコアは、ExACによって報告されていない。
ヒトノックアウト遺伝子のよりロバストなセットを取得し、CHMの付加価値を実証するために、CHMを、61KのDiscovEHR参加者の間で見出された3,915のホモ接合性pLoFと組み合わせた。pLoF-pLoFのCHMは、推定KOを有する≧1及び≧10の個人の遺伝子数を、それぞれ15%及び54%だけ増加させた(上記の表6を参照)。KO分析にCHMを含めることの利点は、タンパク質の機能を混乱させると予測されるミスセンスバリアントを考慮した場合、さらに一層重要となった。つまり、CHMは、キャリアが1つ以上の遺伝子を28%増加させ、キャリアが10以上の遺伝子を246%増加させたが、この場合、遺伝子の両方のコピーは完全にノックアウトするかまたは攪乱されることが予測される。
トリオ検証の結果は、家族関係ベースのフェージングは100%で正確(750/750pCHM)であり、EAGLEフェージングは91.4%(459/502pCHM;上記表3参照)で、それほど正確ではないことを示した。また、Illuminaのリードデータの視覚的検証を、190のpCHM(115のシス及び79のトランス;126がEAGLEでフェージングし、74が家系/関係でフェージング)に対して実施した。視覚的検証は、家系/関係及びEAGLEのフェージングに対して、全体的な精度がそれぞれ95.8%及び89.9%であった(以下の表9を参照)。Illuminaのリードベースの検証結果は、トリオ検証の結果と一致していたが、Illuminaのリードベースの検証精度の結果は、トリオを用いるフェージングによって決定されたフェージング精度よりも低かったことに留意すべきである。この差は、シークエンシング及びバリアントコールのエラーを起こしやすいエクソンの小さな問題を抱えた領域での偽陽性pCHMの濃縮による可能性が高いと考えられる。
Figure 0007277438000014
61KのDiscovEHR参加者の中から、両方のバリアントが互いに75塩基対以内にある200のpCHMをランダムに選択し、次いで2つのバリアントにまたがるリードスタックを見ることによって、相を視覚的に検証した。10個(5%)については、両方のバリアントにオーバーラップするリードが無いか、またはリードが矛盾する結果(すなわち、一部のリードがシスを示し、別の一部のリードがトランスを示した)をもたらすかいずれかのために、リードスタックを用いて確信的にフェージングすることができなかった。
de novo変異
de novo変異(DNM)は、純化淘汰の低減に起因して、ヒトにおいて極端な表現型を生成する可能性が高い稀少バリエーションのクラスである。最近のシークエンシング研究は、DNMがヒトの遺伝的疾患の大きな要因であることを示しており(de Ligt et al.(2012),N.Engl.J.Med.367,1921-1929、Deciphering Developmental Disorders Study(2017).Prevalence and architecture ofde novo mutations in developmental disorders.Nature Publishing Group 542,433-438、及びFromer et al.(2014),Publishing Group 506,179-184)、DNMが遺伝子機能をより深く理解するための有益なツールであることを示している。
DiscovEHRデータセットから再構築された核家族を使用して、トリオで1,262の有効な子供のうちの887に分布する1,800の中信頼度及び高信頼度エクソンDNMを確信的にコールした(実施例12参照)。個人当たりのDNMの平均数は1.43で、最大は49であった(図21C)。PolyPhen2は、DNMの28.2%(N=507)を「ほぼ確実に損傷を与える(probably damaging)」と予測し、さらに8.6%(N=154)を「損傷を与えるおそれがある(possibly damaging)」と予測した。DNMは、1,597の遺伝子にわたって分布し(図21D)、1つの遺伝子のみが5超のDNMを収容していた。最も一般的なタイプのDNMは、非同義SNV(57.17%)であり、これに同義SNV(25.56%)が続いた。以下の表10は、DNMタイプの全内訳を提供し、様々な機能クラスに分類されるDNMの割合が、発育障害のある子供のDNMについての最近の研究で見出された割合とよく一致していることを示す。
Figure 0007277438000015
発達障害解読研究(DDD)(Deciphering Developmental Disorders Study(2017).Prevalence and architecture of de novo mutations in developmental disorders.Nature Publishing Group 542,433-438)。DDD論文はまた、我々の分析にもこの表にも含まれていない他のクラスの57のDNMを報告している。パーセンテージは、それに応じて調節した。
全ての機能クラスにまたがる、23の高信頼度DNM、30の中信頼度DNM、及び47の低信頼度DNMの視覚的検証を行おうと試みた。8つの中信頼度バリアント及び2つの低信頼度バリアントを、真陽性または偽陽性のDNMとして確信的にコールすることができなかった。残りの23/23(100%)の高信頼度DNM、19/22(86%)の中信頼度DNM、及び12/43(28%)の低信頼度DNMは、真陽性であると確証できた。また、視覚的検証では、10を超えるDNMを持つ個人の潜在的なDNMの大多数(40/49)が、偽陽性のコールである可能性が高いことも確認した。
家系におけるバリアントと表現型との分離
DiscovEHRデータセットの中から再構築された家系データを使用して、新規/稀少な集団のバリエーションと家族性バリアントとを区別し、この家系データを活用して、集団全体の関連解析で正当に評価されていない、家族において分離する高浸透性疾患のバリアントを同定した。これは、これらの家系を介して伝わる全ての既知のメンデル遺伝病の原因となるバリエーションの調査を意図したものではないが、家族性大動脈瘤(図23、パネルA)、QT延長症候群(図23、パネルB)、甲状腺癌(図23、パネルC)、家族性高コレステロール血症(FH;図24)(Maxwell,E.K.,et al.(2017).Profiling copy number variation and disease associations from 50,726 DiscovEHR Study exomes)など、いくつかの説明に役立つ実例を確認した。FHの例は、LDLRの新規家族性高コレステロール血症を引き起こす縦列重複の27/29キャリアが再構築されたことから、特に興味深いものであった。5つのさらなるキャリア(描画せず)もこの家系に含まれていた。LDL及び総コレステロールの上昇、ならびに冠動脈疾患及び早期発症虚血性心疾患(「虚血性心疾患年齢」は男性では55歳未満、女性では65歳未満)の有病率の増加は、重複キャリアと分離する。彼らの共通の祖先の履歴は、彼ら全員が、この重複事象を、約6世代前の共通の祖先から継承したことを示している。残りの2つのサンプルは互いに第一度血縁者であったが、それらは正常に遺伝子型決定されず、より大きな家系に接続することができなかった。
シークエンシング研究は、絶えず増え続ける大きさのヒト集団を捕集し、配列決定し続け、ヒトの近縁性の非常に複雑で絡み合った性質を明らかにしつつある。DiscovEHRデータセットでは、約35Kの第一度及び第二度の関係を同定し、7,684家系を再構築し、7,000人を超える参加者の第二度家族ネットワークを明らかにした。創設者集団の研究は、関係の複雑さを前々から強調しており(オールド・オーダー・アーミッシュ(McKusick,V.A.,HOSTETLER,J.A.,and EGELAND,J.A.(1964).GENETIC STUDIES OF THE AMISH,BACKGROUND AND POTENTIALITIES.Bull Johns Hopkins Hosp 115,203-222)、ハッタライト(Ober et al.(2001),The American Journal of Human Genetics 69,1068-1079),及びアシュケナジム系ユダヤ人 (Gusev et al.(2012),Mol.Biol.Evol.29,473-486)、非創始者集団の最近の研究は、広範なレベルの近縁性(UK Biobank(Bycroft et al.(2017).Genome-wide genetic data on ~500,000 UK Biobank participant),NHAMES(Malinowski et al.(2015),Front Genet 6,317)と、AncestryDNA(Han et al.(2017),Nat Commun 8,14238.)とを報告している。かっては大規模なシークエンシングコホートに少数の個人のみが関与していたが、密接した関係は、大規模な医療集団ベースのゲノム(HPG)研究で、大多数ではないにしても、大きな割合の個人を関わらせる可能性が高い。本明細書では、シミュレーションと実データとを通じて、多数の密接した家族関係、核家族、及び情報価値のある家系を取得できることを示した。この観察は、家族が同じヘルスケアシステムを訪問する傾向があり、遺伝的及び環境的疾患のリスクが類似しているため、HPG研究のために収集されたデータセットでは、より顕著である可能性が高かった。コホート全体のごく一部に過ぎないことを知っていて、近縁関係にある個人のペアを関連研究から単純に除外することは、もはや不可能であることが明らかになりつつある。最大サイズの非血縁関係のセットを取得する従来のアプローチは、HPGコホートのサイズを劇的に減少させ、これらのタイプのコホートで実施される多くの重要な疾患表現型分析には適さない。代わりに、本明細書に概説した近縁性情報を活用するために新しい方法が必要である。
本研究では、近縁性情報を活用する方法について、いくつかの方法を実証している。第一に、稀少複合ヘテロ接合性変異(CHM)のフェージング精度を改善した。EAGLEを用いてCHMの正確なフェージングが得られたが、家系及び関係ベースのフェージングは遥かに正確であり、pCHMのフェージングエラーは推定で31%減少した。ペアワイズ関係を使用したフェージングでは、2つのバリアントが2人の血縁者に一緒に現れる場合、それらはシスにあり、共通の祖先から一緒に分離されていると想定しているため、pCHMの関係に基づくフェージングの精度は、MAFが1%を超えるバリアントが含まれるとき、わずかに低下する可能性があった。独立して分離された2つの共通のバリアントが、複数の人々に同時に出現する可能性は非常に高く、その結果、アルゴリズムによって誤ってシスとしてフェージングされる。一般的なバリアントでは、集団アレル頻度を使用したフェージングの方が、関係ベースのフェージングよりも適切な場合がある。
第二に、HPG研究で同定した関係の家系再構築によって、いくつかの方法で使用できる有益なトリオと情報価値のある家系とを提供した。1,262の再構築したトリオを、1,800のDNMを見出すのに使用し、拡張した家系を通じて、既知の疾患原因変異を追跡することが可能であった。情報価値のある家系の数とサイズは、集団のより大きい部分をシークエンシングするにつれて増加し続け、より豊富な家系データセットを提供するようになる。家系を介した稀少バリアントの伝播は、それが実在し、より伝統的なメンデル遺伝学的アプローチの使用を可能にするという強力な証拠を提供するので、家系及び関係は、非常に稀少なバリアントに対して特に有用である。家系は、フォローアップ調査で小規模から中規模の家系を増補するために、DiscovEHRの患者に再連絡し、追加の家族を採用する能力と組み合わせると、特に有用であることが判明した。
近縁性を対処する必要のある厄介な問題として見るのではなく、貴重で未開発の遺伝学的洞察の源を利用する機会として捉えるべきである。ゲノムに基づく精密医療の時代が始まり、ますます成長するシークエンシングコホートに含まれる家族構造と遠縁の近縁性とを効果的にマイニングすることができる革新的な方法とツールとが必要とされている。
実施例3.2
近縁性を拡張DicoverEHRデータセットの確率変数のように扱うことの代わりに、近縁性を活用すること
拡張DiscovEHRデータセット内の12,574の第一度家族ネットワークの家系構造を、家系再構築ツールPRIMUSを使用して再構築した。LBD推定及び報告された参加者の年齢を考慮すると、これらの家系の98.9%が単一の家系構造に明確に再構築されることが分かった。これらの家系には、2,192の核家族(1,841のトリオ、297のカルテット、50のクインテット、3つのセクステット、及び1つのセプテット)が含まれる。表11は、トリオの祖先別の内訳を示す。図14、パネルCは、34のシークエンシングされた個人を含む最大の第一度家系を示す。
Figure 0007277438000016
複合ヘテロ接合性変異
稀少ヘテロ接合変異のペアで構成される57,355の高クオリティpCHMは、推定LoF(pLoF、すなわち、ナンセンス変異、フレームシフト変異、またはスプライスサイト変異)、または有害である強い証拠を有するミスセンスバリアントであることを認めた。次に、EAGLEによるアレル頻度ベースのフェージングと、再構築された家系及び関係のデータによる家系ベースのフェージングとの組み合わせを使用することによるpCHMのフェージングを実行した(図25)。トリオ検証は、EAGLEが、平均89.1%の精度でpCHMをフェージングしたことを示した(以下の表12)。ただし、このコホートでは、広範な家系及び関係のデータのため、pCHMの25.2%をフェージングし、非常に正確なトリオ及び関係のフェージングデータ(R 98.0%;表12)を有するトランスCHMの33.8%が、トランスCHMの不正確なフェージングを約3分の1だけ減らした。フェージングされたpCHMは、シングルトンから1%MAFまでの全範囲に及んだ(以下の表13参照)。
Figure 0007277438000017
トリオの子に発生したMAF<1%かつMAC>1の全てのpCHMを、再構築したトリオを用いてフェージングし、「真」であると見なした。子において、寄与バリアントの1つ以上がde novoであると判定されたpCHMは除外した。次に、フェージング済みトリオのpCHMを使用して、他方のフェージング方法を評価した。EAGLEの精度は、再構築した各核家族から1人の子供の第一度血縁者を全て除去し、次いで残りのデータセット中の全てのバリアントをフェージングすることによって評価した。EAGLEのフェージング済みpCHMを、トリオのフェージング済みpCHMと比較した。
Figure 0007277438000018
処理後、pCHMの40.3%をトランスにてフェージングし、92Kの個人の17,533の間に分布する、20,947の稀少で有害なCHMの高信頼度セットを産出した(1人当たり平均約0.23;1人当たり最大約10;図26、パネルA)。シスのpCHMバリアント間のゲノム距離の中央値(5,955bp)は、トランスでのpCHM間の中央値距離(11,600bp;図27)の半分より少し多かった。CHMのほぼ3分の1は少なくとも1つのpLoFを含み、CHMの8.9%は2つのpLoFバリアントで構成されていた(以下の表14を参照)。19,467の標的化遺伝子のうち4,216を超える遺伝子が1つ以上のCHMキャリアを含み(以下の表15を参照)、そのうちの2,468は複数のキャリアを有している(図26、パネルB)。ExACのpLIスコアは、CHMキャリアが125を超える10の遺伝子が、ゲノム内で最もLoF耐性がある可能性が高いことを示す。(以下の表16を参照)。
Figure 0007277438000019
この表は、稀少(<1%MAF)pLoFとミスセンスバリアントとで構成されているCHMの内訳を提供する。また、これらのCHMのいくつがインデル-インデル、インデル-SNP、及びSNP-SNPペアリングを構成しているかも示す。
Figure 0007277438000020
Figure 0007277438000021
CHMが最も多い10の遺伝子については、ExAC3によって報告されたそれらのpLIスコアを表示している。また、全ての遺伝子をpLIスコアでランク付けし、報告pLiスコアを持つ遺伝子の総数で割ることによって計算した、LoF許容値に対する各遺伝子のパーセンタイルも表示している。
SSPOのpLIスコアは、ExACによって報告されていない。
遺伝子の両方のコピーが同じ個人でノックアウトまたは攪乱される、よりロバストな遺伝子セットを得るとともに、CHMの付加価値を実証するために、CHMを、92KのDiscovEHR参加者の間から見出した6,560の稀少(MAF<1%)ホモ接合性pLoFと組み合わせた。pLoF-pLoFのCHMは、R1及びR20の個人でノックアウトされた遺伝子の数を、それぞれ15%及び61%だけ増加させた(以下の表16を参照)。CHMをKO分析に含めることの利点は、タンパク質機能を攪乱させると予測されるミスセンスバリアントを考慮した場合、さらに一層重要になる。組み合わされた20,364の稀少ホモ接合性pLOFと、有害ミスセンスバリアントとが、92Kの参加者に見つかった。ホモ接合性pLoFまたは予測される有害ミスセンスバリアントのキャリアは、完全にノックアウトまたは攪乱されると予測される多数の遺伝子を提供した。ただし、CHMのキャリアを含めると、R1個人でノックアウトまたは攪乱される遺伝子が26%多く、R20個人でノックアウトまたは攪乱される遺伝子が397%多かった(表15)。
de novo変異
92Kの拡張DiscovEHR参加者から再構築された核家族は、トリオで2,602の有効な子供のうちの1,783に分布する3,415の中信頼度及び高信頼度エクソンDNMを内密にコールすることができた(平均約1.31;最大約48;図26、パネルC)。PolyPhen2は、DNMの29.1%(n~995)を「ほぼ確実に損傷を与える(probably damaging)」と予測し、さらに9.2%(n~316)を、損傷を与えるおそれがある(possibly damaging)と予測している。DNMは2,802の遺伝子にわたって分布しており(図26、パネルD)、TTNが最も多く(9)収容する。DNMの最も一般的なタイプは、非同義SNV(58.5%)であり、これに同義SNV(24.3%)が続く。表17は、DNMタイプの全内訳を提供し、様々な機能クラスに分類されるDNMの割合が、発育障害のある子供のDNMについての最近の研究で見出された割合とほぼ一致することを示す。図7に記載の通り、DNMのコール、フィルタリング、及び信頼度順位付けのワークフローに従った。シークエンシングされた92,455のGHSエクソームのコホートから、2,602トリオを同定した(図28)。低信頼度、中信頼度、及び高信頼度のDNMに基づいてソートしたトリオから、6,645のエクソンDNMを同定した。拡張DiscovEHRデータセットから再構築した家族を使用して、2,602の家族トリオから、3,409の中信頼度及び高信頼度のエクソンDNM、ならびに3,045の単一ヌクレオチドDNMを確信的にコールした(図29、パネルA及びB)。コホートのほとんどの個人は、DNMは5未満であった。さらに、シークエンシングした92,455のGHSエクソームのコホートから、2,602のトリオを同定し、これらを低信頼度、中信頼度、高信頼度のバリアントに基づいてソートして、10,000のランダムバリアントを生成する73、192の中信頼度/高信頼度バリアントを提供した。
最も一般的なタイプのDNMは、非同義SNVであり、これに同義SNVが続いた。ストップロスSNVは、最も一少ないDNMであった。この結果は、61Kのエクソームシークエンシングデータを含むDiscovEHRコホートで得られた結果と類似していた(以下の表17を参照)。図30は、拡張DiscovEHRコホートで見出される中信頼度及び高信頼度のエクソンDNM(n=3409)のタイプの全内訳を提供し、異なる機能的効果クラスに分類されるDNMの割合を示す。
Figure 0007277438000022
発達障害解読研究(DDD)(Deciphering Developmental Disorders Study(2017).Prevalence and architecture of de novo mutations in developmental disorders.Nature Publishing Group 542,433-438)。DDD論文はまた、我々の分析にもこの表にも含まれていない他のクラスの57のDNMを報告している。パーセンテージは、それに応じて調節した。
図31は、拡張DiscovEHRコホートで見出した中信頼度及び高信頼度のエクソンDNM(n=3409)のタイプの内訳を提供し、トランジション、トランスバージョン、及びインデルによって引き起こされるDNMの比率を示す。拡張DiscovEHRで見出した中信頼度及び高信頼度のエクソンDNM(n=3409)のうち、トランジションによる変異の数は2038であり、トランスバージョンによる変異の数は1007であり、インデルによる変異の数は364であった。したがって、トランジションのトランスバージョンに対する比(Ti:Tv)は2:1であり、これは、他の研究から得られるトランジションのトランスバージョンに対する比と類似していた。単一ヌクレオチドDNM(n=3045)の中では、システインからチミン、グアニンからアデニンが最も一般的な変異であった(図32)。
中信頼度及び高信頼度のDNMは、常染色体の全体にわたって均一に分布していた。一方向のカイ二乗検定(χ検定)は、10Mのエクソン塩基対当たりのDNMが、ランダム分布(p=0.045)から大きく逸脱しないことを示した(図33)。
CGジヌクレオチド(慣習的にCpGと記されており、「p」は2つの塩基間のリン酸塩を表す)の変異は、ヒトにおける疾患の原因となる生殖細胞系変異の3分の1の原因である(Cooper and Krawczak(1990);Hum.Genet.85:55-74)。中信頼度/高信頼度DNM(n=3,409)のうち、約13%のDNMがCpGアイランドでのDNMによるものと説明された。ランダムなバリアント(n=10,000)のうち、約10%のDNMがCpGアイランドにおけるDNMに起因するものであることが説明された。DNMは、ランダムなバリアントよりもCpGアイランドで発生する可能性が高かった(χ=32.3661、df値=1、p値=1.28E-08)(図34)。このことは、これらのサイトの高可変性のためにより予想される。
全ての機能クラスにまたがる、23の高信頼度DNM、30の中信頼度DNM、及び47の低信頼度DNMの視覚的検証を行う試みがなされた。8つの中信頼度バリアント及び2つの低信頼度バリアントを、真陽性または偽陽性のDNMとして確信的にコールすることができなかった。残りの23/23(100%)の高信頼度DNM、19/22(86%)の中信頼度DNM、及び12/43(28%)の低信頼度DNMは、真陽性であると確証した。また、視覚的検証では、10を超えるDNMを持つ個人の潜在的なDNMの大多数(40/49)が、偽陽性のコールである可能性が最も高いことも確認した。
家系におけるバリアントと表現型との分離
92Kの拡張DiscovEHR参加者の中から再構築された家系データを使用して、稀少集団バリエーションと家族性バリアントとを区別し、これを活用して、家族において分離する高浸透性疾患のバリアントを同定した。これは、DiscoverEHRデータセットと同様に、これらの家系を介して伝わる全ての既知のメンデル遺伝病の原因となるバリエーションの調査を意図したものではないが、家族性大動脈瘤、QT延長症候群、甲状腺癌、家族性高コレステロール血症(FH[MLM:143890];図35)を同定している。CNVコールの更新時に、92Kのエクソーム間でFHを引き起こす縦列重複の37のキャリアを見つけている。これに基づいて、37のキャリアのうち30を単一の拡張家系に再構築した。キャリアの共通の祖先の履歴は、彼ら全員が、この重複事象を、約6世代前の共通の祖先から継承したことを示した。残りの7人のキャリアのうち2人は互いに第二度血縁者であるが、残りの7人のキャリアも図36の他のキャリアと遠縁であることを確認するために、遺伝子型決定の配列データを入手することはできなかった。図36に記載された家系では、キャリア及び非キャリアの状態は、各個人からのエクソームデータから決定され、最大LDLレベル(記号の下の値)の上昇、ならびに冠動脈疾患(CAD、赤い塗りつぶし)の有病率の増加及び重複キャリアで分離された純粋な高コレステロール血症(ICD 272.0;青)が見られた。5人のさらなるキャリア(描かれず)は、この家系の個人の遠縁の血縁者(第七度から第九度の血縁者)であることも判明した(図36)。
実施例4
患者とサンプル
予測モデルを、コホート、(A)61,720人の匿名化された患者のエクソームを持つDicovEHRコホート、及び(B)92,455人の匿名化された患者のエクソームを持つ拡張DicovEHRコホートに適用することにより、2セットのデータを収集した。
Geisinger Health System(GHS)から取得した両方のコホートの匿名化された患者参加者全員のシークエンシングを行った。全ての参加者が、MyCode(登録商標)CommunityHealthInitiative(Carey et al.(2016),Genet.Med.18,906-913)に参加することに同意し、Regeneron-GHS DiscovEHR研究(Dewey et al.(2016),Science 354,aaf6814-aaf6814)のゲノム解析に、DNAサンプルを提供した。全ての患者のエクソームは、対応する匿名化電子医療記録(EHR)にリンクさせた。シークエンシングされた最初の50,726人の詳細については、以前に公開されている(Dewey et al.(2016),Science 354,aaf6814-aaf6814、Abul-Husn et al.(2016),Science 354,aaf7000-aaf7000)。
この研究は、この研究に参加する家族を具体的に対象としたものではなく、ヘルスケアシステムと頻繁にやり取りする慢性的な健康上の問題を抱えた成人、ならびに冠動脈カテーテル検査室及び肥満治療室からの参加者を対象としたものであった。
実施例5
サンプル準備、シークエンシング、バリアントコール、及びサンプルQC
サンプルの準備とシークエンシングは、Deweyら(Dewey et al.(2016),Science 354,aaf6814-aaf6814)で以前に記載されている。
シークエンシングの完了時に、各Illumina Hiseq 2500ランからの生データがローカルバッファストレージに収集され、自動解析のために、DNAnexusプラットフォーム(Reid et al.(2014)15,30)にアップロードされた。サンプルレベルのリードファイルが、CASAVA(Illumina Inc.,San Diego,CA)で生成され、BWA-memでGRCh38(Li及びDurbin(2009);Bioinformatics 25,1754-1760、Li,H.(2013);arXiv q-bio.GN)にアライメントされた。結果のBAMファイルを、GATK及びPicardを使用して処理して、重複のソート、マーク付け、及び推定インデル周辺のリードの局所的なアライメント修正を行った。シークエンシングされたバリアントは、転写産物及び遺伝子に対する機能的影響を明らかにするために、Ensembl85遺伝子定義を使用して、snpEFF(Cingolani et al.(2012);Fly(Austin)6,80-92)でアノテーションを付けられた。遺伝子定義は、19,467の遺伝子に対応する、アノテーション付きの開始及び停止を伴うタンパク質コードである54,214の転写産物に制限した。
高頻度のホモ接合性、低配列データカバレッジ、または実際の一卵性双生児であることを確認できない遺伝的に同定された重複により示される低クオリティのDNA配列データを有する個人を除外した。61,019のエクソームが分析のために残った。サンプル準備、シークエンシング、バリアントコール、及びバリアントアノテーションに関する追加情報は、Dewey et al.(2016),Science 354,aaf6814-1~aaf6814-10に報告されている。
実施例6
主成分と祖先推定
PLINKv1.9を使用してデータセットをHapMap3(International HapMap 3 Consortium,Altshuler et al.(2010)、Nature Publishing Group 467,52-58)とマージし、両方のデータセットにあるSNPのみを保持した。また、次のPLINKフィルタを適用した。--maf 0.1 --geno 0.05 --snps-only --hwe 0.00001。HapMap3サンプルの主成分(PC)分析が計算され、PLINKを使用して各サンプルがデータセット内のそれらのPCに投入された。HapMap3サンプルのPCを使用して、5つの祖先スーパークラス、すなわち、アフリカ人(AFR)、混合アメリカ人(AMR)、東アジア人(EAS)、ヨーロッパ人(EUR)、南アジア人(SAS)のそれぞれのカーネル密度推定(KDE)を訓練した。KDEを使用して、各サンプルが各スーパークラスに属する尤度を計算した。各サンプルについて、祖先スーパークラスを尤度に基づいて割り当てた。サンプルが0.3を超える尤度を持つ2つの祖先グループを有する場合、AFRをEURに優先して割り当て、AMRをEURに優先して割り当て、AMRをEASに優先して割り当て、SASをEURに優先して割り当て、AMRをAFRに優先して割り当てた。それ以外の場合は「不明」になる(これは、EURとEASの集団の厳密な推定と、データセット内の混合の多い集団の包括的な推定とを提供するために行われた)。ゼロまたは3つ以上の祖先グループが十分に高い尤度を持っていた場合、サンプルに祖先の「不明」を割り当てた。祖先が不明なサンプルは、祖先ベースの同祖性(IBD)計算から除外した。
実施例7
IBD推定
次のフラグを使用して全データセットに対してPLINKを実行することにより、高クオリティで一般的なバリアントをフィルタ処理した。--maf 0.1 --geno 0.05 --snps-only --hwe 0.00001。次に、エクソームデータから正確なIBD推定を取得するために、2つのアプローチを取った。最初に、個人間のIBD推定を、祖先分析から決定された祖先スーパークラス(例えば、AMR、AFR、EAS、EUR、及びSASなど)内で計算した。次のPLINKフラグを、第二度関係までのIBD推定を取得するために使用した。--genome --min 0.1875。全てのサンプルが同様の祖先アレルを共有したため、これにより、より正確な関係推定が可能になった。ただし、このアプローチでは、異なる系統の個人間の関係、例えば、ヨーロッパ人の父親とアジア人の母親との子供を予測できなかった。
次に、祖先の異なる個人間の第一度の関係を把握するために、--min 0.3 PLINKオプションを使用して、全ての個人間でIBD推定を計算した。次に、個人を第一度家族ネットワークにグループ化し、この場合、ネットワークノードが個人であり、エッジが第一度第一度関係である。各第一度家族ネットワークを、prePRIMUSパイプライン(Staples et al.(2014);Am.J.Hum.Genet.95,553-564)を通じて実行し、サンプルの祖先を適切な祖先マイナーアレル頻度に一致させて、IBD推定を改善した。このプロセスは、各家族ネットワーク内の個人間の第一度及び第二度の関係を正確に推定した(0.15の最小PI_HAT)。
最後に、上記の2つのアプローチからのIBD推定を、家族ネットワーク由来のIBD推定から欠落している関係を、祖先ベースのIBD推定に加えることによって結合した。このアプローチは、類似した祖先の全てのサンプル間の第二度の関係と、全てのサンプル間の第一度の関係とに、正確なIBD推定をもたらした。
第三度血縁者のIBD比率は、分析が人為的に誇張されたIBD推定に起因して予測される第三度関係の過剰数をもたらすことが多いため、多様な系統を持つ大規模なエクソームシークエンシングデータセットから正確に推定することは困難である。祖先特異的IBD分析の間中に、--min 0.09875カットオフを使用して、DiscovEHR及び拡張DiscovEHRコホートに存在し得る第三度の関係の数を把握したが、これらはいずれのフェージングまたは家系ベースの分析にも使用しなかった。むしろ、本明細書で開示された関係に基づく分析では、第一度及び第二度の家族ネットワーク内で同定した信頼度の高い第三度の関係のみを使用した。
実施例8
家系再構築
DiscovEHR及び拡張DiscovEHRコホート内で同定した全ての第一度家族ネットワークをPRIMUSv1.9.0で再構築した。結合IBD推定は、遺伝的に得られた性別及びEHR報告年齢の範囲内でPRIMUSに提供された。PI_HAT>0.375の近縁性カットオフを、再構築が第一度家族ネットワークに制限されるように指定し、0.1875のカットオフを指定して、第二度ネットワークを定義した。
実施例9
アレル頻度ベースのフェージング
EAGLEv2.3(Loh et al.(2016);Nat Genet 48,1443-1448)を用いて、61,019エクソーム由来の全ての両アレルバリアントをフェージングした。DNAnexus内の解析を並列化するために、ゲノムは、500のバリアントと250Kの塩基対との最小オーバーラップを有する約40Kのバリアントのオーバーラップセグメントに分割した。目標は、遺伝子内の推定複合ヘテロ接合性変異をフェージングすることであったため、セグメント切断点が遺伝子間領域に生じるように注意を払った。
EAGLEが提供するgenetic_map_hg19.txt.gzファイルのhg19からGRCh38へのリフトオーバーを実行し、染色体を交換し、または染色体内の相対的な順序を変更して、染色体の位置とcMの位置とがどちらもソート時に昇順にならないようにする全てのバリアントを除去した。ほとんどの場合、このQCステップは、セントロメア周辺の逆位を除去することであった。代替染色体にマッピングしたSNPも除去した。合計で、330万のSNPのうち、2,783のみを遺伝子マップファイルから除去した。各セグメントのデータは、PLINK形式のファイルとしてEAGLEに提供され、次のEAGLEコマンドラインパラメータを使用して、DNAnexus上で実行した。
--geneticMapFile=genetic_map_hg19_withX.txt.GRCh38_liftover.txt.gz
--maxMissingPerIndiv 1
--genoErrProb 0.01
--numThreads=16
実施例10
複合ヘテロ接合性コール
目的は、機能喪失(pLoF)と推定されるバリアントの高信頼度複合ヘテロ接合性変異(CHM)コールを取得して、ノックアウトまたは攪乱された可能性のある遺伝子の両方のコピーを有するヒトを同定することであった。フレームシフト、終止コドン獲得、終止コドン喪失、開始コドン獲得、開始コドン喪失、またはスプライシングアクセプタもしくはドナー変更バリアントが発生した場合に、バリアントをpLoFとして分類した。2番目の、潜在的に有害なバリアントの拡張されたセットを作成し、これには、pLOF及びおそらく攪乱ミスセンスバリアントが含まれ、以下の5つの方法の全てによって有害であると予測されることによって定義された。SIFT(Loh et al.(2016);Nat Genet 48,1443-1448)(損傷を与える(damaging))、PolyPhen2 HDIV(損傷を与える(damaging)、及び損傷を与えるおそれがある(possibly damaging))、PolyPhen2 HVAR(損傷を与える(damaging)、及び損傷を与えるおそれがある(possibly damaging))、LRT(有害(deleterious))、及びMutationTaster(Schwarz et al.(2014);Nat.Methods 11、361-362)(必然的に疾患を引き起こす(disease causing automatic)、及び疾患を引き起こす(disease causing))。
稀少(代替アレルカウント<1%)潜在複合ヘテロ接合性変異(pCHM)は、同じ人物の遺伝子内のヘテロ接合性pLoF及び/または有害ミスセンスバリアントの全ての可能な組み合わせをテストすることによって同定した。ハーディ・ワインベルグ平衡(HWE)から外れており(PLINKで計算されたp値<10-15(Chang et al.(2015);Gigascience4,7.))、61Kサンプル全体で10%の欠損を超えているか、または同じ個人の10塩基対内に別のバリアントがある全てのバリアントを除外した。また、QD<3、AB<15%、及びリードデプス<7であるSNPと、QD<5、AB<20%、及びリードデプス<10であるINDELとを除外した。フィルタ処理後、25,031の個人に分布し、pCHMバリアントがトランスでフェージングされた場合に、ヒトの遺伝子の両方のコピーの機能をノックアウトまたは攪乱する可能性のある、39,459の高クオリティpCHMを得た。
次のステップは、pCHMをフェージングすることであった。EAGLEを用いる集団アレル頻度ベースのフェージングと家系/関係ベースのフェージングとの組み合わせを使用して、pCHMがシスかトランスかを判定した。図9は、DiscovEHRデータセット内の各pCHMの最も正確なフェージングを取得するために使用されたpCHMフェージングワークフローを示す。図2は、拡張DiscovEHRデータセット内の各pCHMの最も正確なフェージングを取得するために使用されたpCHMフェージングワークフローを示す。家系及び関係のフェージングは、EAGLEのフェージングよりも正確であることが証明されたため、家系及び関係のデータをフェージングに使用した。以下の表18は、様々なタイプの家族関係のpCHMの相を決定するために使用されるロジックを示す。残りの全てのpCHMでは、上記のEAGLEフェーズデータを使用した。シングルトンでのEAGLEのフェージング精度は、ランダムな推測と有意差がなかったため、EAGLEでフェージングした、バリアントのいずれかまたは両方がシングルトンであるpCHMを除外した(DiscovEHRデータセットについては以下の表19を、拡張DiscovEHRデータセットについては表20を参照)。DiscovEHRデータセットでは、pCHMの2つのバリアントが、100未満の同じマイナーアレルカウント(MAC)を有する場合、それらはほとんど常にシス(我々のトリオの間で37中の36)であり、EAGLEのpCHMフェージングの精度を超えていた。拡張DiscovEHRデータセットでは、pCHMの2つのバリアントが、100未満の同じマイナーアレルカウント(MAC)を有する場合、それらはほとんど常にシス(トリオの間で22中の22)であり、EAGLEのpCHMフェージングの精度を超えていた。
Figure 0007277438000023
Figure 0007277438000024
これらのルールは、各関係の最上位のルールから順に適用した。「?」結果は、pCHMをフェージングできなかったことを意味する。「NA」の結果は、その結果が発生してはならなかったことを示しており、シークエンシングエラーまたは他の非メンデリアンのバリアントの伝播の結果である可能性が高い。PC_relは、親子関係の非pCHMキャリアを指す。「稀少」とは、1%未満のMAFを指し、これには、本明細書で用いられる全てのバリアントが含まれる。
Figure 0007277438000025
MAFが1%未満の全てのpCHMは、pCHMを構成する2つのバリアントのうち頻度の低いものによってビニングした。pCHMのEAGLEフェージングとトリオで決定されたフェージングとを比較することにより、正しいコールと精度とを決定した。トリオ内でシスまたはトランスであると判定された誤ったEAGLEフェージングのpCHMの数も提供される。トリオの子で一方または両方のバリアントがde novoであると判断されたpCHMを除外した。MAC>6のpCHMは全て90年代半ばから上位にかけて同様の精度であったが、MACが2~6の場合に精度の低下が見られた。EAGLEのシングルトンフェージングは、ランダムな推測よりもパフォーマンスが大幅に向上していなかったため、EAGLEでフェージングされたシングルトンは、EAGLEでフェージングされたpCHMの全体的な精度を測定するときと同様に、フェージングされたpCHMの結果から除外した。
Figure 0007277438000026
MAFが1%未満の全てのpCHMは、pCHMを構成する2つのバリアントのうち頻度の低いものによってビニングした。pCHMのEAGLEフェージングとトリオで決定されたフェージングとを比較することにより、正しいコールと精度とを決定した。また、トリオ内でシスまたはトランスであると判定された誤ったEAGLEフェージングのpCHMの数も提供した。トリオの子供で一方または両方のバリアントがde novoであると判定されたpCHMは除外した。MAC>9のpCHMは全て同様の精度90を有するが、2~9のMACを有するEAGLEのpCHMフェージング精度が低下することが分かる。EAGLEのシングルトンを含むpCHMのフェージングは、ランダムな推測よりもパフォーマンスが大幅に低下しないため、EAGLEでフェージングされたシングルトンは、EAGLEでフェージングされたpCHMの全体的な精度を測定する場合と同様に、フェージングされたpCHMの結果から除外した。シングルトンバリアントのフェーズ精度が低いため、シングルトンバリアントを含む2,838のpCHMを除去した。したがって、残りの401個のシングルトンバリアントは、トリオ及び関係のデータのみでフェージングした。
データセット全体にわたるEAGLEのpCHMフェージングの正確な精度を測定するために、フェージング前に各核家族の1人の子供の全ての第一度の血縁者を除外して、データセット全体にわたってEAGLEを実行した。親のハプロタイプを含めると、データセットに親のないサンプルと比較して、トリオの子のフェージングの精度が向上するため、この取り除きが必要であった。
最後に、個人の同じ遺伝子内に複数のpCHMがあった場合、最も有害なプロファイルを持つpCHMのみを保持した(以下の表21を参照)。全てのpCHMの>99%をフェージングし、13,335の稀少複合ヘテロ接合性変異(CHM)を同定することが可能であった。
Figure 0007277438000027
人が同じ遺伝子に2つ以上のトランスpCHMを持っている場合に、この表の値を、最も有害なpCHMを同定して保持するために使用した。効果スコアは、2つのバリアントの機能的効果スコアを加算し、それらが全ての遺伝子転写産物に影響しなかった場合にペアにペナルティを課すことによって計算した。スコアのより低いpCHMは、最も有害であると予測され、保持された。
実施例11.1
DiscovEHRデータセットの複合ヘテロ接合性変異検証
フェージングの精度は、フェージング予測をトリオ及びIlluminaのリードを用いて行ったフェージングと比較することによって評価した。まず、トリオのフェーズ済みのpCHMを真として使用して、pCHMのフェージング精度を評価した。各家族関係のフェージングアプローチは、トリオのフェージングとは独立して行われたため、pCHMキャリアがトリオの子供である限り、各関係クラスのフェージングの精度を適切に測定することができた。上記の表4及び表12は、これらの稀少pCHMについて、家族関係ベースのフェージングの精度が100%正確であったことを示している。EAGLEのフェージングは、DiscovEHR及び拡張DiscovEHRデータセットで、それぞれ91.4%及び89.1%と正確性が劣っていた。DiscovEHRデータセットの場合、フェージングpCHMでのEAGLEの精度が、異なるマイナーアレル頻度の範囲で評価され、MACが6より大きい場合に95%以上、MACが2~6の場合、約77%の精度が一貫して得られた(上記の表19参照)。EAGLEフェージングは、予期したように、シングルトンでは機能しなかった。
次に、統合ゲノミクスビューア(IGV)(Robinson et al.(2011);Nat.Biotechnol.29,24-26)のリードスタックを見て、2つのバリアントが同じリードで発生するか、独立して発生するかを確認することにより、短いIlluminaリード(約75bp)で200pCHMを検証することを試みた。この検証プロセス中に、最初の欠失の終わりが2番目の欠失の10bps内にある2つの欠失で構成されるpCHMは、実際には1つの大きな欠失が2つの個別の欠失として誤ってコールされていることが分かった(39,459のpCHMのうちN=1,109)。15のみがトランスとしてフェージングされたため(CHMデータセット全体の約0.1%)、これらのpCHMは分析全体から除外されなかったが、検証のために200のpCHMを選択したときに除外した。短いリードを使用して、ランダムに選択された200のpCHMのうち、190を決定的にフェージングするためにリードを使用することができた。残りの10例は、シス及びトランスの両方のフェージングのリード証拠を示したが、これは、バリアントの一方または両方が偽陽性コールである可能性が最も高い。
実施例11.2
拡張DiscovEHRデータセットの複合ヘテロ接合性変異検証
DiscovEHRデータセットの場合、上記の表12は、稀少pCHMの家族ベースのフェージングの精度が99.6%(1060/1064pCHM)であったことを示している。EAGLEのフェージングは、89.1%(766/860pCHM;上記の表12)であり、それほど正確ではなかった。マイナーアレル頻度の様々な範囲でのEAGLEのpCHM位相精度を評価して、MACが9を超える場合にEAGLEが一貫して90%を超える精度を達成し、MACが2から9の場合に約77%の精度を達成することを確認した(上記の表20を参照)。EAGLEフェージングは、シングルトンでは不十分に機能した。
次に、統合ゲノミクスビューア(IGV)(Robinson et al.(2011);Nat.Biotechnol.29,24-26)のリードスタックを見て、2つのバリアントが同じリードで発生するか、独立して発生するかを確認することにより、短い(975bp)Illuminaリードで200pCHMを検証することを試みた。短いリードを使用して、190(115のシス及び79のトランス、EAGLEフェーズ済みの126及び家系または関係でフェーズ済みの74)のpCHMを選択した。残りの10例は、シス及びトランスのフェージングの両方のリード証拠を示したが、これは、バリアントの一方または両方が偽陽性のコールであったためという可能性が最も高い。視覚的検証では、家系及び関係のフェージング及びEAGLEのフェージングに対して、全体的な精度がそれぞれ95.8%及び89.9%であった(表22を参照)。Illuminaのリードベースの検証結果は、トリオの検証結果と一致していたが、Illuminaのリードベースの検証精度の結果は、トリオのフェージングの精度よりも低かった。この差は、シークエンシングエラーとバリアントコールエラーが発生しやすい問題を抱えた小さなエクソン領域での偽陽性pCHMの濃縮による可能性が高い。
Figure 0007277438000028
200のpCHMは、92Kの拡張DiscovEHR参加者の中からランダムに選び、両方のバリアントは互いに75塩基対以内であり、2つのバリアントにまたがるリードスタックを調べることによって、相を視覚的に検証した。10個(5%)については、両方のバリアントにオーバーラップするリードが無いか、またはリードが矛盾する結果(すなわち、一部のリードがシスを示し、別の一部のリードがトランスを示した)をもたらすかいずれかのために、リードスタックを用いて確信的にフェージングすることができなかった。
実施例12
de novo変異(DNM)検出
DNMを検出するための2つの異なるアプローチの結果を統合した。最初の方法はTrioDeNovo(Wei et al.(2015);Bioinformatics 31,1375-1381)であり、これは、子の変異性サイトごとに、子及び親の遺伝子型の尤度を読み取る。これらの尤度をベイジアンフレームワークに入力し、子供のバリアントがDNMであるという事後尤度を計算した。第2のプログラムはDeNovoCheck(https://sourceforge.net/projects/denovocheck)であり、これは、de Ligt,et al.の補足方法で説明されている(de Ligt et al.(2012);N.Engl.J.Med.367,1921-1929)。DeNovoCheckは、どちらの親にもなく、子でコールされていると識別された候補DNMのセットを取り込む。次に、BAMファイルを調べることによって、子にバリアントが存在するかどうか、及び両方の親にバリアントが存在しないかどうかを確認する。これらの潜在的DNMをフィルタリングし、結合セットにおける各DNMの信頼度レベルを種々のQCメトリクスを用いて評価した。図7は、このDNMコールプロセスを示し、適用したバリアントフィルタを示し、各DNMを低信頼度、中信頼度、または高信頼度のいずれかに分類するために使用した基準を提供する。
実施例13.1
DiscovEHRデータセットでの受胎時の親の年齢と子供のDNMの数との相関関係のテスト
この分析では、10以上のDNMを持つサンプルを外れ値(N=6の除外されたサンプル)として除外した。それらは、技術上の人為結果、または体細胞変異を示している可能性が高い。母方と父方の年齢には高い相関がある(rho=0.78、p=1.2x10^-262);共同でモデル化した場合、共線性のためにどちらも有意ではなかった(0.0053母体DNM/年、p=0.48;0.0076父親DNM/年、p=0.26;ポアソン回帰)(図36A及び36B)。次に、両親の年齢差(父方と母方の年齢)を、出生時の母方または父方の年齢のいずれかと共にテストした。父方と母方の両方の年齢は、DNM数を等しく予測することが判明した(つまり、年齢差は、母方または父方の年齢を所与として、DNMの数と有意に関連しなかった)。
エクソンDNM数の増加も、他の報告と一致して、出生時の母性(0.012DNM/年、p=0.011;ポアソン回帰;図37)と父性(0.011DNM/年;p=0.007)の両方の年齢に関して観察された(Deciphering Developmental Disorders Study(2017).Nature 542,433-438;Kong et al.(2012)Nature 542,433-438;Rahbari et al.(2016)Nat.Genet.48,126-133;及びWong et al.(2016)Nat.Commun.7,10486)。特に、出生時の母方と父方の年齢はデータセットにおいて高度に相関し(rho=0.78、p=1.2x10^-262;図38)、従って、比率は付加的ではなく、駆動因子として区別するための有意差は同定されなかった。
実施例13.2
拡張DiscovEHRデータセットでの受胎時の親の年齢と子供のDNMの数との相関関係のテスト
拡張DiscoverEHRコホートは、受胎時の親の年齢と子供のDNMの数との相関関係をテストした結果、DiscovEHRコホートと同様の結果を示した。他の報告と一致して、出生時の母性(0.011DNM/年、p=7.3x10-4;ポアソン回帰;図37)及び父性(0.010DNM/年;p=5.6x10-4)の年齢の両方に関して、エクソンDNMの数の増加が観察された。特に、出産時の母方及び父方の年齢はデータセット(r=0.79;図39)で高度に相関して、したがって、比率は付加的ではなく、駆動因子としても区別される有意差はない。
父方の年齢は、ポアソン分布を使用して、1人あたりのDNM数と相関していた(n=2587、係数=0.010、p=5.67E-4)。同様に、母性年齢は、ポアソン分布を使用して、1人あたりのDNM数と相関した(n=2587、係数=0.011、p=7.35E-4)。さらに、父方と母方の年齢も互いに相関している(R=0.79;p<10E-308)。
機能的予測アルゴリズム、すなわちSIFT(損傷を与える(damaging))、PolyPhen2HDIV(損傷を与える(damaging)、及び損傷を与えるおそれがある(possibly damaging))、PolyPhen2HVAR(損傷を与える(damaging)、及び損傷を与えるおそれがある(possibly damaging))、LRT(有害(deleterious))、及びMutationTaster(Schwarz et al.(2014);Nat.Methods 11、361-362)(必然的に疾患を引き起こす(disease causing automatic)、及び疾患を引き起こす(disease causing))を使用して、DNMの病原性が予測された。DNMの病原性予測は、ランダムなバリアント分布の予測とは大きく異なる(図38)。より高い割合のDNMも、非病原性の一致した予測を示した。DNMは、5/5アルゴリズムによって病原性であると予測される可能性が1.8倍高くなる。ランダムバリアントは、病原性の予測が一致しない可能性が1.5倍高い。
実施例14
LDLR縦列重複遠縁家系推定
本コホートでは匿名化された個人の真の家族歴を知ることは不可能であるが、PRIMUS再構成家系(Staples et al.(2014);Am.J.Hum.Genet.95,553-564.)、ERSA遠縁関係推定(Huff et al.(2011);GenomeRes.21,768-774.)、及び家系を接続するPADREの能力を用いて(Staples et al.(2016);The American Journal of Human Genetics 99,154-162)、LDLRにおける新規縦列重複の突然変異キャリアの最良の家系表現を同定した(Maxwell et al.(2017).Profiling copy number variation and disease associations from 50,726 DiscovEHR Study exomes)。HumanOmniExpress配列データは、以前はより遠い関係を推定するために使用されていた。
実施例15
SimProgeny
SimProgenyは、1つ以上のサブ集団に分散した何百万もの人々の集団をシミュレートし、数百年にわたって彼らの子孫を追跡できる。単純性と現実性の間の良好なバランスを見つけるために、ユーザが調整できるいくつかの主要な集団レベルパラメータを選択した(下の表23を参照)。これらのパラメータは、シミュレーションツールを比較的単純に保ちながら、実際の集団と家族の家系構造とを適切に近似するために選択された。デフォルト値は、米国の人口統計に基づいている。デフォルト値は両方のコホートで機能するように設定されており、SimProgenyコード(Webリソース)に含まれている構成ファイルを変更することで、これらのパラメータを簡単にカスタマイズして、様々な集団をモデル化できる。集団シミュレーションプロセスの詳細については、例17を参照されたい。
Figure 0007277438000029
DiscovEHRコホート用に開発されたフレームワークの場合、生殖の終了は49歳であり、拡張DiscovEHRコホート用に開発されたフレームワークの場合、生殖の終了は50歳であった。
集団のモデリングに加えて、SimProgenyは、遺伝的研究のために集団から個人を選択することをモデル化する2つの確認アプローチ、すなわちランダム確認及びクラスタ化サンプリングをシミュレーションする。ランダム確認は、集団中の各個人が置換なしに確認される等しい機会を与える。クラスタ化サンプリングは、近親血縁者を濃縮するアプローチであり、第一度及び第二度の血縁者の多くとともにランダムに個人を選択することによって行われる。第一度血縁者の数は、ユーザ指定の第一度の確認ラムダ(デフォルトは0.2)で分布されたポアソンから値をサンプリングすることによって決定される。第二度血縁者の数も同じ方法で決定され、デフォルトの第二度確認ラムダは0.03である。SimProgenyの確認オプションの詳細については、実施例17を参照されたい。
実施例16
基礎となるDiscovEHR集団のシミュレーションとその確認
シミュレーションモデルが過度に複雑にならないように、シミュレーションには、開始サイズが200K、300K、400K、450K、500K、550K、600K、及び1000Kの個々の集団が含まれていた。SimProgenyパラメータ(上記の表23を参照)は、公的に入手可能な国、州、及び郡レベルのデータと、GHSを通じて個人が確認された方法に関する我々自身の理解とに合わせて調整された。選択したパラメータの入手源は、補足ファイルSimulation_parameters.xlsで見ることができる。GHSは主に農村地域を対象としており、都市部より移住率が低い傾向があるため、ペンシルベニア州(PA)の平均的な移住率を減少させた。シミュレーションは120年のバーンイン期間で実行し、その後101年間継続した。シミュレートされた集団は約15%増加したが、これは、20世紀半ば以降のPAの成長と類似している。
ランダムの確認とクラスタ化された確認との両方を実行した。両確認アプローチに対して、集団の最初の5%の確認順序(ordered_sampling_proportanceパラメータで指定)を、共同研究の開始時にGHSバイオバンクの個人のランダムシークエンシング順序をモデル化するためにシャッフルした。このパラメータの選択は、ランダム確認には影響せず、クラスタリングされた確認でのペアワイズ関係の累増に対する影響は無視できたが、クラスタ化サンプリングで確認されたデータセットにおいて、変曲点(より高いラムダ値でより顕著である)を生じることにより、1つ以上の血縁者を持つ個人の割合には影響した。この変曲点は、実際のデータのフリーズプロセスをモデル化し、またバイオバンクからのシークエンシングサンプルと新しく確認された個人との間のよりスムーズな移行をモデル化したりする場合には、それほど顕著ではないであろう。
実施例17
SimProgeny集団及び確認シミュレーションプロセス
シミュレーションは、ユーザが指定した数のサブ集団及びサイズを初期化することによって開始させた。年齢は最初に0と最大繁殖年齢(デフォルトは49)との間に割り当てた。集団中の個人は、青年、繁殖力のある個人、高齢者という3つの年齢別プールのいずれかに存在していた。個人は、繁殖年齢(15のデフォルト)未満であればサブ集団の青年プールに割り当てられ、繁殖年齢の範囲内(デフォルトは15~49)であればサブ集団の交配プールに割り当てられた。個人は最低生殖可能年齢以上になると幼体プールから交配プールに移動した。同様に、繁殖可能年齢を超えた時点で、交配プールから高齢者の集団に移した。移住または死亡した個人は、全ての年齢層から除外した。初期集団を確立した後、シミュレーションは120年のバーンイン段階を実施して、出生数と死亡数とが等しく、正味移動率がゼロであることを必要としながら、入力パラメータとより密接に一致する家族関係と年齢分布とを確立した。バーンイン後、シミュレーションを、与えられた個人群成長と移動速度とで特定の年数で実行した。シミュレーションは1年ずつ増加し、特に明記しない限り、各サブ集団内で実施された以下のステップを毎年実施した。
1.年齢-年齢プールから年齢を超えた個人を次の年齢プールに移動する。
2.交際-一夫一婦の結婚に参加する独身男性と独身女性とをシミュレートする。このプロセスは、完全同胞関係の現実的な数を取得するために重要である。男性と女性のペアは、単一の生殖年齢の男女プールからランダムに選択され、彼らは、男性及び女性の「年齢別の結婚」パラメータで指定された年齢で結婚する可能性に基づいて正常に結婚する。ペアは、結婚率で定義された成功した結婚の数に達するまで引き出される。夫婦は、従兄弟よりも遠縁に限定される。バーンインフェーズ中は、ユーザが指定した初期の結婚率に達するまで、結婚率は2倍になる(デフォルトは、結婚している交配プールの66%である)。
3.分割-指定された離婚率で婚姻関係を解消する男女をシミュレートする。カップルはランダムに選択され、両方の個人は独身としてマークされる。
4.交わり-1年間に集団内で発生する可能性のある全ての生殖をシミュレートする。母親と父親のペアは、単一の出生年齢プールまたは既婚プールのいずれかから無作為に完全同胞率で定義された比率で選択される(デフォルトは、全ての出産の88%が既婚カップルである)。ペアが引き出され、目標とする受胎成功数に到達するまで、生殖の試みが行われる(デフォルトの出生率は1人あたり0.0219出生である)。受胎が成功する可能性は、母親の年齢と対応する出生率とに基づいている。両親は従兄弟よりも遠縁に制限されており、全ての個人は年に1人の子供を持つことが制限されている。
5.淘汰-亡くなった個人をシミュレートする。死亡率(デフォルトは1人あたり0.0095人の死亡)を使用して、特定の年に集団内で予想される死亡数を決定する。年齢別の男性と女性の死亡率パラメータを使用して、ランダムに選択された個人が亡くなる可能性に重みを付ける。0から1の間の乱数が、その人の年齢で死亡する確率を超える場合、その個人は保持され、別の個人が選択される。不運な個人は、死者のプールに追加され、生者の他のプールから除去される。120歳を超える全ての個人は、自動的に死者のプールに追加され、年間の目標死亡数にカウントされる。
6.移住-集団間の移住をシミュレートする。移住は、交配プールから無作為に1個人を選び、個人が結婚していて繁殖可能な年齢であれば、配偶者とともに集団から外すことによって行われる。青年と高齢者が去る割合は、繁殖可能な年齢の夫婦の数とともに記録される。移民は、年齢分布と繁殖可能な年齢の夫婦の数を維持する方法で行われる。まず、既存の集団から少年をランダムに選択し、同じ性別と年齢の新しい個人を青年プールに追加し、適切な割合の少年が追加されるまでこのプロセスを繰り返す。同じプロセスが高齢者に対しても繰り返される。次に、2人の繁殖力のある年齢の者が既存の人口から選択され、2人の新しい個人が対応する年齢とともに追加される。一方は男性、もう一方は女性に割り当てられ、2人の移民は結婚する。このステップは、夫婦の数が補充されるまで繰り返される。最後に、繁殖力のある年齢の者が、新しい青年を追加するのと同じプロセスで追加され、移民の目標数が達成されるまで繰り返される。このプロセスは、人口の年齢と性別の分布、及び結婚している繁殖可能な年齢の者の割合を維持するのに役立つ。
7.移植-サブ集団内で移動する人々をシミュレートする。サブ集団間の遺伝的隔離がないことをシミュレートするために、個人は集団全体の中でサブ集団間を移動することができる。集団全体にわたって単一の移動率が使用される。サブ集団からの個人はランダムに選択され、移植の所望の数が達成されるまで、他のサブ集団の1つにランダムに割り当てられる。この手順は、サブ集団が1つしかない場合、または移植率が0の場合は発生しない(デフォルトは、全集団移植の1%である)。
このシミュレーションは、各創始者とその子孫を追跡しながら、指定された時間の間進行する。
ランダムの確認とクラスタ化された確認との両方を実行した。両確認アプローチに対して、集団の最初の5%の確認順序(ordered_sampling_proportanceパラメータで指定)を、共同研究の開始時にGHSバイオバンクの個人のランダムシークエンシング順序をモデル化するためにシャッフルした。このパラメータの選択は、ランダム確認には影響せず、クラスタリングされた確認でのペアワイズ関係の累増に対する影響は無視できたが、クラスタ化サンプリングで確認されたデータセットにおいて、変曲点(より高いラムダ値でより顕著である)を生じることにより、1つ以上の血縁者を持つ個人の割合には影響した。この変曲点は、実際のデータのフリーズプロセスをモデル化し、またバイオバンクからのシークエンシングサンプルと新しく確認された個人との間のよりスムーズな移行をモデル化したりする場合には、それほど顕著ではないであろう。ユーザは、次のサブ集団のセットに移動する前に、1つ以上のサブ集団から確認することをシミュレートしたい場合に、サブ集団確認順序を指定することができる。デフォルトでは、最初に全てのサブ個人群をグループ化し、それらのサブ個人群をあたかも単一の個人群であるかのように確認する。ユーザは、他のサブ集団または全集団に移る前に確認された集団の最初の割合を指定することもできる。このプログラムは、集団全体の出力をpedファイル形式で作成し、確認されたサンプルのリストを確認された順に作成し、有用な集団及び確認統計を要約したいくつかの結果ファイルを作成した。
実施例18
所与の表現型の遺伝的原因の同定を補助するために家系構造を用いる方法は、典型的には、関連マッピング、連鎖解析、またはその両方に関する革新的なバリエーションを含む。このような方法には、MORGAN31、pVAAST15、FBAT(www.hsph.harvard.edu/fbat/fbat.htm)、QTDT(csg.sph.umich.edu/abecasis/qtdt/)、ROADTRIPS、rareIBD、及びRV-GDTが含まれる。使用する適切な方法は、表現型、遺伝様式、祖先、家系構造/サイズ、家系の数、及び非血縁関係のデータセットのサイズに依存する。関連性と家系を用いて遺伝子-表現型の関連性を直接調べることに加えて、それらは、追加のまたは改善されたデータを生成するための多くの他の方法、すなわち家系を認識した帰属、家系を認識したフェージング、メンデル式エラーチェック、複合ヘテロ接合性ノックアウト検出及びde novo変異コール、ならびにバリアントコール検証で使用することもできる。
本開示は、上記で説明し例示した例示的な実施形態に限定されるものではなく、添付の特許請求の範囲内で変形及び修正が可能である。

Claims (14)

  1. 集団における遺伝的バリアントを、前記集団の近縁性を活用することによってフェージングする方法であって、
    サンプルの1つ以上のそれぞれについて祖先スーパークラスの指定を確立すること、
    祖先スーパークラス内での被験者の第1の同祖性推定を生成すること、
    被験者の祖先スーパークラスから独立した被験者の第2の同祖性推定を生成すること、
    前記第2の同祖性推定の1つ以上に基づいて被験者を第1の第一度家族ネットワークにクラスタリングすること、
    第1の第一度家族ネットワーク内での被験者の第3の同祖性推定を生成すること、
    第1の同祖性推定と第3の同祖性推定とを統合して、統合同祖性推定を取得すること、
    統合同祖性推定に基づいて被験者の第2の第一度家族ネットワークを構築すること、ならびに
    バリアントを、統合同祖性推定及び第2の第一度家族ネットワークに応じて、複合ヘテロ接合性変異(CHM)であるかもしくはCHMではないものとしてフェージングする、またはバリアントを、統合同祖性推定及び第2の第一度家族ネットワークに応じて、de novo変異(DNM)として同定することを含む、前記方法。
  2. 第1の同祖性推定と第3の同祖性推定とを統合することは、前記第3の同祖性推定に特有のペアワイズ同祖性推定を補って、前記第1の同祖性推定を増強することを含む、請求項1に記載の方法。
  3. バリアントを複合ヘテロ接合性変異(CHM)としてフェージングすることは、
    集団アレル頻度に応じてバリアントをフェージングすること、
    ハーディ・ワインベルグ平衡(HWE)から外れたバリアント、もしくは同じサンプル内の別のバリアントの10塩基対内のバリアント、またはその両方を除去すること、及びクオリティバイデプス(quality by depth)(QD)が2以下、またはリードデプス(DP)が5未満、または代替アレルバランス(AB)が10%以下、またはこれらの組み合わせのSNPを除去すること、及びQDが2以下、またはDPが5未満、またはABが10%以下、またはこれらの組み合わせの挿入または欠失(INDEL)を除去すること、
    同じサンプル内及び同じ遺伝子内に1つ以上のバリアントの対が存在する潜在複合ヘテロ接合性変異(pCHM)として残りのバリアントを選択すること、ならびに
    pCHMをシスpCHMまたはトランスpCHMとしてフェージングし、次いでトランスpCHMとしてフェージングされた前記pCHMをCHMとして分類することを含む、請求項1または請求項2に記載の方法。
  4. バリアントを複合ヘテロ接合性変異としてフェージングすることは、
    ハーディ・ワインベルグ平衡(HWE)から外れたバリアント、もしくは同じサンプル内の別のバリアントの10塩基対内のバリアント、またはその両方を除去すること、及びクオリティバイデプス(quality by depth)(QD)が3以下、またはリードデプス(DP)が7未満、または代替アレルバランス(AB)が15%以下、またはこれらの組み合わせのSNPを除去すること、及びQDが5以下、またはDPが10未満、またはABが20%以下、またはこれらの組み合わせの挿入または欠失(INDEL)を除去することを含む、請求項3に記載の方法。
  5. 表現型の重症度に応じてCHMにスコアを付与すること、及び
    サンプル当たりの遺伝子ごとに最も高い表現型重症度スコアを有するCHMを選択し、前記サンプルのデータセット内の個人が同一遺伝子内に複数のCHMを有する場合に、タンパク質機能抑制をもたらす可能性が最も高い前記CHMが同定されるようにすることをさらに含む、請求項3または請求項4に記載の方法。
  6. バリアントをde novo変異としてフェージングすることは、
    第2の第一度家族ネットワーク及びその両親及び対応する子に由来するトリオのサンプルの中のバリアントを同定すること、
    遺伝子型尤度スコアを、トリオ内の親サンプル及び対応する子サンプルのバリアントに割り当て、前記バリアントがde novo変異である確率を計算し、計算した前記確率が統計的に有意な場合に、前記バリアントを推定de novo変異として同定すること、
    トリオの子サンプルのバリアントを同定し、前記トリオのどちらの親サンプルにも前記バリアントが存在しない場合に、前記バリアントを推定de novo変異として同定すること、
    以下の(1)~(6)の条件:(1)前記子サンプルに35未満の遺伝子型クオリティ(GQ)アノテーションを有する、(2)前記サンプル全体にわたって10以上の代替アレルカウント(AC)を有する、(3)前記子サンプルに、7未満のリードデプス(DP)と4未満の代替DPとを有する、またはどちらか一方の前記親サンプルに2%を超えるアレルバランス(AB)を有する、(4)前記子サンプルに15%未満のアレルバランス(AB)を有する、(5)前記子サンプルに90%を超えるABを有する、または(6)どちらか一方の前記親サンプルに代替アレルホモ接合性を有する、のうちの少なくとも1つあるいはそれらの組み合わせを満たす、推定de novo変異として同定された変異を除去することによって同定される推定de novo変異をフィルタ処理すること、ならびに
    フィルタ処理済みの同定された推定de novo変異を組み合わせて、それによって推定de novo変異データセットを形成することを含む、請求項1~5のいずれか1項に記載の方法。
  7. 前記推定de novo変異データセット内の推定de novo変異を、前記推定de novo変異のアレルバランスが、前記子サンプルで0.15以上、かつ各前記親サンプルで0.02以下の場合、及び前記推定de novo変異のマッピングクオリティが40未満ではない場合、及び前記推定de novo変異のクオリティバイデプス(quality by depth)(QD)の値が2未満ではない場合、及び前記推定de novo変異のMACが前記サンプル全体にわたって20未満の場合、及び前記推定de novo変異のソフトクリップリードが、前記推定de novo変異のキャリアにおける前記変異が生ずるゲノムの遺伝子座で3以下の場合、及び前記推定de novo変異が、4を超えるモノポリマーランを持つINDELではない場合に、中信頼度de novo変異として分類することをさらに含む、請求項6に記載の方法。
  8. 中信頼度de novo変異を、前記中信頼度de novo変異が、前記親サンプルにおいて90以上の遺伝子型クオリティアノテーションを有する場合、及び前記中信頼度de novo変異が、各前記親サンプルにおいて10以上のリードデプスを有する場合、及び前記中信頼度de novo変異が、前記子サンプルにおいて7以上の代替リードデプスを有する場合、及び前記中信頼度de novo変異が、3より大きいQDのSNPを有する場合、及び前記中信頼度de novo変異が、5より大きいQDのINDELを有する場合に、高信頼度de novo変異として分類することをさらに含む、請求項7に記載の方法。
  9. 低クオリティサンプルを前記サンプルのデータセットから除去することをさらに含み、前記低クオリティサンプルは、D-statが0.12を超えるか、もしくは20xのリードカバレッジが75%未満であるか、またはその両方であるサンプルである、請求項1~8のいずれか1項に記載の方法。
  10. データプロセッサ、前記データプロセッサに結合されたメモリ、及び前記メモリに格納されたプログラムを備え、前記プログラムが、
    複数のヒト被験者由来のDNA配列サンプルの中のバリアントを同定すること、
    同定されたバリアントに基づいて被験者の祖先スーパークラスの指定を確立すること、
    祖先スーパークラス内での被験者の第1の同祖性推定を生成すること、
    被験者の祖先スーパークラスから独立した被験者の第2の同祖性推定を生成すること、
    前記第2の同祖性推定の1つ以上に基づいて被験者を第1の第一度家族ネットワークにクラスタリングすること、
    第1の第一度家族ネットワーク内での被験者の第3の同祖性推定を生成すること、
    第1の同祖性推定と第3の同祖性推定とを統合して、統合同祖性推定を取得すること、
    統合同祖性推定に基づいて第2の第一度家族ネットワークを構築すること、
    集団アレル頻度に応じて、サンプルの中のバリアントをフェージングすること、
    フェージングされたバリアントを、同じ被験者及び遺伝子における2つ以上のバリアントの存在に基づいて潜在CHMとして分類すること、ならびに
    潜在CHMを、同じ被験者及び遺伝子における別バリアントとシスであるまたはトランスであるとしてフェージングし、次いでトランスであるとしてフェージングされた前記潜在CHMをCHMとして分類するための命令を含む、システム。
  11. 前記プログラムは、被験者の祖先スーパークラスの指定が確立される前に、同定されたバリアントをフィルタ処理するための命令を含み、
    前記バリアントをフィルタ処理することは、複数の前記ヒト被験者由来の前記サンプル全体にわたって10%を超える代替アレル頻度を有するバリアントを除去すること、またはp値が10 -6 超である、ハーディ・ワインベルグ平衡(HWE)に反するバリアントを除去すること、または複数の前記ヒト被験者由来の前記サンプルの5%超でコールが欠損するバリアントを除去すること、またはこれらの組み合わせを含む、請求項10に記載のシステム。
  12. 前記プログラムは、被験者の第2の同祖性推定が生成される前に、同定されたバリアントをフィルタ処理するための命令を含む、請求項10または請求項11に記載のシステム。
  13. 前記プログラムは、同定されたバリアントがフィルタ処理された後に、低クオリティサンプルを除去するための命令を含む、請求項11~1のいずれか1項に記載のシステム。
  14. 低クオリティサンプルは、D-statが0.12を超えるか、もしくは20xのリードカバレッジが75%未満であるか、またはその両方であるサンプルである、請求項13に記載のシステム。
JP2020513709A 2017-09-07 2018-09-07 ゲノムデータ解析に近縁性を活用するシステム及び方法 Active JP7277438B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2023076708A JP2023103325A (ja) 2017-09-07 2023-05-08 ゲノムデータ解析に近縁性を活用するシステム及び方法

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201762555597P 2017-09-07 2017-09-07
US62/555,597 2017-09-07
PCT/US2018/049967 WO2019051238A1 (en) 2017-09-07 2018-09-07 SYSTEMS AND METHODS OF OPERATING PARENTAL IN GENOMIC DATA ANALYSIS

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2023076708A Division JP2023103325A (ja) 2017-09-07 2023-05-08 ゲノムデータ解析に近縁性を活用するシステム及び方法

Publications (3)

Publication Number Publication Date
JP2020532989A JP2020532989A (ja) 2020-11-19
JP2020532989A5 JP2020532989A5 (ja) 2021-10-14
JP7277438B2 true JP7277438B2 (ja) 2023-05-19

Family

ID=63684599

Family Applications (3)

Application Number Title Priority Date Filing Date
JP2020513709A Active JP7277438B2 (ja) 2017-09-07 2018-09-07 ゲノムデータ解析に近縁性を活用するシステム及び方法
JP2020513751A Pending JP2020533679A (ja) 2017-09-07 2018-09-07 ヒト集団における関連性を予測するシステム及び方法
JP2023076708A Pending JP2023103325A (ja) 2017-09-07 2023-05-08 ゲノムデータ解析に近縁性を活用するシステム及び方法

Family Applications After (2)

Application Number Title Priority Date Filing Date
JP2020513751A Pending JP2020533679A (ja) 2017-09-07 2018-09-07 ヒト集団における関連性を予測するシステム及び方法
JP2023076708A Pending JP2023103325A (ja) 2017-09-07 2023-05-08 ゲノムデータ解析に近縁性を活用するシステム及び方法

Country Status (11)

Country Link
US (3) US11605444B2 (ja)
EP (2) EP3679575A1 (ja)
JP (3) JP7277438B2 (ja)
KR (3) KR20200050992A (ja)
CN (3) CN111247599A (ja)
AU (3) AU2018327331A1 (ja)
CA (2) CA3075266A1 (ja)
IL (2) IL272927A (ja)
MX (2) MX2020002665A (ja)
SG (2) SG11202001715YA (ja)
WO (2) WO2019051233A1 (ja)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
PL3207481T3 (pl) * 2014-10-14 2020-05-18 Ancestry.Com Dna, Llc Redukowanie błędu w przewidywanych powiązaniach genetycznych
JP2021536635A (ja) 2018-09-07 2021-12-27 リジェネロン・ファーマシューティカルズ・インコーポレイテッドRegeneron Pharmaceuticals, Inc. 家系エンリッチメント及び家系内の家族ベース解析のための方法及びシステム
CN110060737B (zh) * 2019-04-30 2023-04-18 上海诚明融鑫科技有限公司 一种基于最大频率虚拟个体的str快速比对方法及系统
CN110974203B (zh) * 2019-12-10 2020-12-04 电子科技大学 一种基于等值排列网络的心率非平衡性分析方法
CN111883207B (zh) * 2020-07-31 2022-08-16 武汉蓝沙医学检验实验室有限公司 一种生物学亲缘关系的鉴定方法
CN116209777A (zh) * 2020-10-27 2023-06-02 深圳华大基因股份有限公司 基于无创产前基因检测数据的亲缘关系判定方法和装置
WO2022109267A2 (en) 2020-11-19 2022-05-27 Regeneron Pharmaceuticals, Inc. Genotyping by sequencing
CN112863601B (zh) * 2021-01-15 2023-03-10 广州微远基因科技有限公司 病原微生物耐药基因归属模型及其建立方法和应用
CN113380328B (zh) * 2021-04-23 2023-06-20 中国人民大学 一种基于mRNA碱基的生物亲缘识别方法和系统
GB2621782A (en) * 2021-05-05 2024-02-21 Univ Leland Stanford Junior Methods and systems for analyzing nucleic acid molecules
CN115101130A (zh) * 2022-06-30 2022-09-23 山东大学 一种基于网络对抗学习的单细胞数据插补方法及系统

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013067001A1 (en) 2011-10-31 2013-05-10 The Scripps Research Institute Systems and methods for genomic annotation and distributed variant interpretation

Family Cites Families (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1288835A1 (en) * 2001-08-31 2003-03-05 Hitachi, Ltd. A method for defining gene hunting studies
WO2006093879A2 (en) * 2005-02-26 2006-09-08 Coco Communications Corporation Naming system layer
US20080055049A1 (en) * 2006-07-28 2008-03-06 Weill Lawrence R Searching methods
US20090177496A1 (en) * 2005-09-28 2009-07-09 Tuck Edward F Relationship prediction system
CA2630070A1 (en) * 2005-11-17 2007-05-31 Motif Biosciences, Inc. Systems and methods for the biometric analysis of index founder populations
WO2008025093A1 (en) * 2006-09-01 2008-03-06 Innovative Dairy Products Pty Ltd Whole genome based genetic evaluation and selection process
JP4778876B2 (ja) * 2006-11-02 2011-09-21 株式会社エー・アンド・デイ エンジン計測装置
US20130080365A1 (en) * 2011-04-13 2013-03-28 The Board Of Trustees Of The Leland Stanford Junior University Phased Whole Genome Genetic Risk In A Family Quartet
WO2013112948A1 (en) * 2012-01-26 2013-08-01 Nodality, Inc. Benchmarks for normal cell identification
US9177098B2 (en) * 2012-10-17 2015-11-03 Celmatix Inc. Systems and methods for determining the probability of a pregnancy at a selected point in time
US20140222349A1 (en) * 2013-01-16 2014-08-07 Assurerx Health, Inc. System and Methods for Pharmacogenomic Classification
EP2949204B2 (en) * 2013-06-14 2020-06-03 Keygene N.V. Directed strategies for improving phenotypic traits
AU2014284180B2 (en) * 2013-06-21 2020-03-19 Sequenom, Inc. Methods and processes for non-invasive assessment of genetic variations
US20150088541A1 (en) * 2013-09-26 2015-03-26 Univfy Inc. System and method of using personalized outcome probabilities to support the consumer in comparing costs and efficacy of medical treatments and matching medical provider with consumer
JP6399511B2 (ja) * 2013-09-30 2018-10-03 学校法人近畿大学 家系図作成支援装置及びコンピュータプログラム
US10468121B2 (en) 2013-10-01 2019-11-05 Complete Genomics, Inc. Phasing and linking processes to identify variations in a genome
US20150112706A1 (en) * 2013-10-17 2015-04-23 Ovuline, Inc. System and Methods for Personal health Analytics Technical Field
WO2015058093A1 (en) * 2013-10-18 2015-04-23 Seven Bridges Genomics Inc. Methods and systems for genotyping genetic samples
CA2937502A1 (en) * 2014-01-27 2015-07-30 Celmatix, Inc. Methods for assessing whether a genetic region is associated with infertility
US9670530B2 (en) 2014-01-30 2017-06-06 Illumina, Inc. Haplotype resolved genome sequencing
US9779460B2 (en) * 2014-02-18 2017-10-03 Marineh Tchakerian Systems, methods and non-transitory computer readable storage media for tracking and evaluating predictions regarding relationships
US20160042282A1 (en) * 2014-08-11 2016-02-11 Rashied Baradaran Amini Relationship evaluator
WO2016049878A1 (zh) * 2014-09-30 2016-04-07 深圳华大基因科技有限公司 一种基于snp分型的亲子鉴定方法及应用
EP3216002A4 (en) * 2014-11-06 2018-04-18 Ancestryhealth.com, LLC Predicting health outcomes
CN106169034B (zh) * 2016-05-26 2019-03-26 中国农业科学院作物科学研究所 基因组信息辅助育种方法i-基于snp聚类信息和pav变异信息的育种亲本选择

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013067001A1 (en) 2011-10-31 2013-05-10 The Scripps Research Institute Systems and methods for genomic annotation and distributed variant interpretation

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Chang C C et al.,Second-generation PLINK: rising to the challenge of larger and richer datasets,Gigascience,2015年,Vol. 4:7
Staples J. et al.,PRIMUS: improving pedigree reconstruction using mitochondrial and Y haplotypes,Bioinformatics,2016年,Vol. 32,pp. 596-598
Tewhey R. et al.,The importance of phase information for human genomics,Nat Rev Genet,2011年,Vol. 12,pp. 215-223

Also Published As

Publication number Publication date
KR20200065000A (ko) 2020-06-08
WO2019051238A1 (en) 2019-03-14
IL272927A (en) 2020-04-30
AU2018327326A1 (en) 2020-03-19
MX2020002665A (es) 2020-10-05
US20190205502A1 (en) 2019-07-04
SG11202001747XA (en) 2020-03-30
CN111279420A (zh) 2020-06-12
JP2020533679A (ja) 2020-11-19
AU2023282218A1 (en) 2024-01-04
CN111247599A (zh) 2020-06-05
JP2020532989A (ja) 2020-11-19
US20190074092A1 (en) 2019-03-07
JP2023103325A (ja) 2023-07-26
EP3679575A1 (en) 2020-07-15
CN117746978A (zh) 2024-03-22
CN111279420B (zh) 2024-01-16
KR102648634B1 (ko) 2024-03-19
MX2020002664A (es) 2020-09-25
IL272963A (en) 2020-04-30
WO2019051233A1 (en) 2019-03-14
KR20240038142A (ko) 2024-03-22
KR20200050992A (ko) 2020-05-12
CA3075266A1 (en) 2019-03-14
US11605444B2 (en) 2023-03-14
EP3679576A1 (en) 2020-07-15
AU2018327331A1 (en) 2020-03-12
CA3075182A1 (en) 2019-03-14
US20220336045A1 (en) 2022-10-20
SG11202001715YA (en) 2020-03-30

Similar Documents

Publication Publication Date Title
JP7277438B2 (ja) ゲノムデータ解析に近縁性を活用するシステム及び方法
Staples et al. Profiling and leveraging relatedness in a precision medicine cohort of 92,455 exomes
Wu et al. A comparison of humans and baboons suggests germline mutation rates do not track cell divisions
Hollenbach et al. Analytical methods for disease association studies with immunogenetic data
CN108277267A (zh) 检测基因突变的装置和用于对孕妇和胎儿的基因型进行分型的试剂盒
Mavengere et al. Paracoccidioides genomes reflect high levels of species divergence and little interspecific gene flow
Tian et al. Estimating the genome-wide mutation rate from thousands of unrelated individuals
Niehus et al. PopDel identifies medium-size deletions jointly in tens of thousands of genomes
Zhang et al. NyuWa Genome Resource: Deep Whole Genome Sequencing Based Chinese Population Variation Profile and Reference Panel
Löytynoja Thousands of human mutation clusters are explained by short-range template switching
Zeng et al. A pipeline for classifying relationships using dense SNP/SNV data and putative pedigree information
Greene Methods for Determining the Genetic Causes of Rare Diseases
Bergström Genomic insights into the human population history of Australia and New Guinea
이선호 New Methods for SNV/InDel Calling and Haplotyping from Next Generation Sequencing Data
Löytynoja Drastic mutations can be hidden in short-read mapping: thousands of mutation clusters in human genome are explained by short-range template switching
Kessler et al. De novo mutations across 1,465 diverse genomes reveal novel mutational insights and reductions in the Amish founder population
Hård et al. Conbase: a software for discovery of clonal somatic mutations in single cells through read phasing
Richter Noncoding Genetic Variants in Congenital Heart Disease
Colucci Next-generation kinship, ancestry and phenotypic deduction for forensic and genealogical analysis
Al-Jumaan et al. Low-pass Whole Genome Imputation Enables the Characterization of Polygenic Breast Cancer Risk in the Indigenous Arab Population
Williams Defining shared/distinct genetic mechanisms for respiratory infections and subsequent respiratory disease
Zeng A pipeline for classifying close family relationships with dense SNP data and putative pedigree information
KEPLER Analyzing NGS Data with Machine Learning-From IBD Segments to Copy Number Variations
Mukhopadhyay Computational methods for calculating meiotic recombination from nuclear pedigrees
Mollon Statistical Analysis of Genome-Genome Interaction with Reference to Kidney Transplant Outcome

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210903

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210903

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20221025

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20230125

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20230404

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20230508

R150 Certificate of patent or registration of utility model

Ref document number: 7277438

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150