(57)【要約】  本開示は、変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得することと、前記少なくとも１つの遺伝子シーケンシング断片の属性情報に基づいて、前記変異遺伝子座候補の、前記遺伝子座の位置に関連する配列特徴および非配列特徴を特定することと、前記配列特徴および前記非配列特徴に基づいて、前記変異遺伝子座候補の遺伝子変異を認識することと、を含む遺伝子変異認識方法、装置および記憶媒体に関する。本開示の実施例は、遺伝子の配列特徴と非配列特徴を組み合わせて、変異遺伝子座の特徴をより全面的に解析し、遺伝子変異認識の正確性を高めることができる。【選択図】図１
2022500773.tif 000002
【特許請求の範囲】
【請求項１】
変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得することと、
前記少なくとも１つの遺伝子シーケンシング断片の属性情報に基づいて、前記変異遺伝子座候補の、遺伝子座の位置に関連する配列特徴および非配列特徴を特定することと、
前記配列特徴および前記非配列特徴に基づいて、前記変異遺伝子座候補の遺伝子変異を認識することと、を含むことを特徴とする、遺伝子変異認識方法。
【請求項２】
前記属性情報は、配列属性情報を含み、
前記少なくとも１つの遺伝子シーケンシング断片の属性情報に基づいて、前記変異遺伝子座候補の配列特徴を特定することは、
前記変異遺伝子座候補の遺伝子位置情報に基づいて、前記変異遺伝子座候補が位置する所定の遺伝子座区間を特定することと、
前記少なくとも１つの遺伝子シーケンシング断片の、前記所定の遺伝子座区間内に位置する各遺伝子座の、遺伝子座の位置に関連しかつ遺伝子属性を表す配列属性情報を取得することと、
前記所定の遺伝子座区間に位置する各遺伝子座の配列属性情報に基づいて、前記変異遺伝子座候補の配列特徴を生成することと、を含むことを特徴とする、請求項１に記載の方法。
【請求項３】
前記少なくとも１つの遺伝子シーケンシング断片の、前記所定の遺伝子座区間内に位置する各遺伝子座の配列属性情報を取得することは、
前記各遺伝子座において前記少なくとも１つの遺伝子シーケンシング断片の遺伝子型を特定することと、
前記各遺伝子座において遺伝子を各遺伝子型毎にカウントすることと、を含むことを特徴とする、請求項２に記載の方法。
【請求項４】
前記少なくとも１つの遺伝子シーケンシング断片の、前記所定の遺伝子座区間内に位置する各遺伝子座の配列属性情報を取得することは、
各遺伝子シーケンシング断片の遺伝子配列と参照ゲノム配列との比較結果に基づいて、前記各遺伝子座において各遺伝子シーケンシング断片の欠失遺伝子の遺伝子型を特定することと、
前記各遺伝子座において前記少なくとも１つの遺伝子シーケンシング断片の欠失遺伝子を各遺伝子型毎にカウントすることと、を含むことを特徴とする、請求項２に記載の方法。
【請求項５】
前記少なくとも１つの遺伝子シーケンシング断片の、前記所定の遺伝子座区間内に位置する各遺伝子座の配列属性情報を取得することは、
各遺伝子シーケンシング断片の遺伝子配列と参照ゲノム配列との比較結果に基づいて、前記各遺伝子座において各遺伝子シーケンシング断片の挿入遺伝子の遺伝子型を特定することと、
前記各遺伝子座において前記少なくとも１つの遺伝子シーケンシング断片の挿入遺伝子を各遺伝子型毎にカウントすることと、を含むことを特徴とする、請求項２に記載の方法。
【請求項６】
前記配列属性情報は、
参照遺伝子の遺伝子型、各遺伝子型毎の遺伝子数、各遺伝子型毎の欠失遺伝子数、各遺伝子型毎の挿入遺伝子数の少なくとも１つを含むことを特徴とする、請求項１から５のいずれか一項に記載の方法。
【請求項７】
前記属性情報は、非配列属性情報を含み、
前記少なくとも１つの遺伝子シーケンシング断片の属性情報に基づいて、前記変異遺伝子座候補の非配列特徴を特定することは、
前記少なくとも１つの遺伝子シーケンシング断片の、遺伝子座の位置に関連しなくかつ遺伝子属性を表す非配列属性情報を取得することと、
前記少なくとも１つの遺伝子シーケンシング断片の非配列属性情報に基づいて、前記変異遺伝子座候補の非配列特徴を特定することと、を含むことを特徴とする、請求項１から６のいずれか一項に記載の方法。
【請求項８】
前記非配列情報は、
比較品質、プラス鎖およびマイナス鎖の偏り、遺伝子シーケンシング断片の長さ、エッジの偏りの少なくとも１つを含むことを特徴とする、請求項７に記載の方法。
【請求項９】
前記少なくとも１つの遺伝子シーケンシング断片の非配列属性情報に基づいて、前記変異遺伝子座候補の非配列特徴を特定することは、
各遺伝子シーケンシング断片の各遺伝子座の比較品質に基づいて、各遺伝子シーケンシング断片の比較品質を特定することであって、前記比較品質は、遺伝子シーケンシング断片の各遺伝子配列毎のシーケンシングの正確性を表すことと、
各遺伝子シーケンシング断片の比較品質に基づいて、前記変異遺伝子座候補に対応する非配列特徴を特定することと、を含むことを特徴とする、請求項８に記載の方法。
【請求項１０】
前記少なくとも１つの遺伝子シーケンシング断片の非配列属性情報に基づいて、前記変異遺伝子座候補の非配列特徴を特定することは、
各遺伝子シーケンシング断片の属する遺伝子鎖がプラス鎖またはマイナス鎖のいずれであるかの情報に基づいて、前記少なくとも１つの遺伝子シーケンシング断片の、プラス鎖である遺伝子鎖とマイナス鎖である遺伝子鎖との割合を特定することと、
前記プラス鎖である遺伝子鎖とマイナス鎖である遺伝子鎖との割合に基づいて、前記変異遺伝子座候補に対応する非配列特徴を特定することと、を含むことを特徴とする、請求項８に記載の方法。
【請求項１１】
前記配列特徴および前記非配列特徴に基づいて、前記変異遺伝子座候補の遺伝子変異を認識することは、
前記配列特徴および前記非配列特徴を統合し、前記変異遺伝子座候補の統合特徴を得ることと、
前記変異遺伝子座候補の統合特徴に基づいて、前記変異遺伝子座候補の遺伝子変異を認識することと、を含むことを特徴とする、請求項１から１０のいずれか一項に記載の方法。
【請求項１２】
前記変異遺伝子座候補の統合特徴に基づいて、前記変異遺伝子座候補の遺伝子変異を認識することは、
前記変異遺伝子座候補の統合特徴に基づいて、前記変異遺伝子座候補の遺伝子の変異性を示す変異値を得ることと、
前記変異値が所定の閾値以上である場合、前記変異遺伝子座候補の遺伝子に変異があると判定することと、を含むことを特徴とする、請求項１１に記載の方法。
【請求項１３】
前記変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得することは、
体細胞遺伝子による遺伝子シーケンシングで得られた遺伝子シーケンシング断片を取得することと、
前記遺伝子シーケンシング断片の遺伝子配列を参照ゲノム配列と比較し、比較結果を得ることと、
前記比較結果に基づいて、前記体細胞遺伝子に異常がある変異遺伝子座候補を特定することと、
前記変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得することと、を含むことを特徴とする、請求項１から１２のいずれか一項に記載の方法。
【請求項１４】
変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得するための取得モジュールと、
前記少なくとも１つの遺伝子シーケンシング断片の属性情報に基づいて、前記変異遺伝子座候補の、遺伝子座の位置に関連する配列特徴および非配列特徴を特定するための特定モジュールと、
前記配列特徴および前記非配列特徴に基づいて、前記変異遺伝子座候補の遺伝子変異を認識するための認識モジュールと、を含むことを特徴とする、遺伝子変異認識装置。
【請求項１５】
前記属性情報は、配列属性情報を含み、
前記特定モジュールは、
前記変異遺伝子座候補の遺伝子位置情報に基づいて、前記変異遺伝子座候補が位置する所定の遺伝子座区間を特定する第１の特定サブモジュールと、
前記少なくとも１つの遺伝子シーケンシング断片の、前記所定の遺伝子座区間内に位置する各遺伝子座の、遺伝子座の位置に関連しかつ遺伝子属性を表す配列属性情報を取得する第１の取得サブモジュールと、
前記所定の遺伝子座区間に位置する各遺伝子座の配列属性情報に基づいて、前記変異遺伝子座候補の配列特徴を生成する第１の生成サブモジュールと、を含むことを特徴とする、請求項１４に記載の装置。
【請求項１６】
前記第１の取得サブモジュールは、具体的に、
前記各遺伝子座において前記少なくとも１つの遺伝子シーケンシング断片の遺伝子型を特定することと、
前記各遺伝子座において遺伝子を各遺伝子型毎にカウントすることと、に用いられることを特徴とする、請求項１５に記載の装置。
【請求項１７】
前記第１の取得サブモジュールは、具体的に、
各遺伝子シーケンシング断片の遺伝子配列と参照ゲノム配列との比較結果に基づいて、前記各遺伝子座において各遺伝子シーケンシング断片の欠失遺伝子の遺伝子型を特定することと、
前記各遺伝子座において前記少なくとも１つの遺伝子シーケンシング断片の欠失遺伝子を各遺伝子型毎にカウントすることと、に用いられることを特徴とする、請求項１５に記載の装置。
【請求項１８】
前記第１の取得サブモジュールは、具体的に、
各遺伝子シーケンシング断片の遺伝子配列と参照ゲノム配列との比較結果に基づいて、前記各遺伝子座において各遺伝子シーケンシング断片の挿入遺伝子の遺伝子型を特定することと、
前記各遺伝子座において前記少なくとも１つの遺伝子シーケンシング断片の挿入遺伝子を各遺伝子型毎にカウントすることと、に用いられることを特徴とする、請求項１５に記載の装置。
【請求項１９】
前記配列属性情報は、
参照遺伝子の遺伝子型、各遺伝子型毎の遺伝子数、各遺伝子型毎の欠失遺伝子数、各遺伝子型毎の挿入遺伝子数の少なくとも１つを含むことを特徴とする、請求項１４から１８のいずれか一項に記載の装置。
【請求項２０】
前記属性情報は非配列属性情報を含み、
前記特定モジュールは、
前記少なくとも１つの遺伝子シーケンシング断片の、遺伝子座の位置に関連しなくかつ遺伝子属性を表す非配列属性情報を取得する第２の取得サブモジュールと、
前記少なくとも１つの遺伝子シーケンシング断片の非配列属性情報に基づいて、前記変異遺伝子座候補の非配列特徴を特定する第２の特定サブモジュールと、を含むことを特徴とする、請求項１４から１９のいずれか一項に記載の装置。
【請求項２１】
前記非配列情報は、
比較品質、プラス鎖およびマイナス鎖の偏り、遺伝子シーケンシング断片の長さ、エッジの偏りの少なくとも１つを含むことを特徴とする、請求項２０に記載の装置。
【請求項２２】
前記第２の特定サブモジュールは、具体的に、
各遺伝子シーケンシング断片の各遺伝子座の比較品質に基づいて、各遺伝子シーケンシング断片の比較品質を特定することであって、前記比較品質は、遺伝子シーケンシング断片の各遺伝子配列毎のシーケンシングの正確性を表すことと、
各遺伝子シーケンシング断片の比較品質に基づいて、前記変異遺伝子座候補に対応する非配列特徴を特定することと、に用いられることを特徴とする、請求項２１に記載の装置。
【請求項２３】
前記第２の特定サブモジュールは、具体的に、
各遺伝子シーケンシング断片の属する遺伝子鎖がプラス鎖またはマイナス鎖のいずれであるかの情報に基づいて、前記少なくとも１つの遺伝子シーケンシング断片の、プラス鎖である遺伝子鎖とマイナス鎖である遺伝子鎖との割合を特定することと、
前記プラス鎖である遺伝子鎖とマイナス鎖である遺伝子鎖との割合に基づいて、前記変異遺伝子座候補に対応する非配列特徴を特定することと、に用いられることを特徴とする、請求項２１に記載の装置。
【請求項２４】
前記認識モジュールは、
具体的に前記配列特徴および前記非配列特徴を統合し、前記変異遺伝子座候補の統合特徴を得る統合サブモジュールと、
前記変異遺伝子座候補の統合特徴に基づいて、前記変異遺伝子座候補の遺伝子変異を認識する認識サブモジュールと、を含むことを特徴とする、請求項１４から２３のいずれか一項に記載の装置。
【請求項２５】
前記認識サブモジュールは、具体的に、
前記変異遺伝子座候補の統合特徴に基づいて、前記変異遺伝子座候補の遺伝子の変異性を示す変異値を得ることと、
前記変異値が所定の閾値以上である場合、前記変異遺伝子座候補の遺伝子に変異があると判定することと、に用いられることを特徴とする、請求項２４に記載の装置。
【請求項２６】
前記取得モジュールは、具体的に、
体細胞遺伝子による遺伝子シーケンシングで得られた遺伝子シーケンシング断片を取得することと、
前記遺伝子シーケンシング断片の遺伝子配列を参照ゲノム配列と比較し、比較結果を得ることと、
前記比較結果に基づいて、前記体細胞遺伝子に異常がある変異遺伝子座候補を特定することと、
前記変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得することと、に用いられることを特徴とする、請求項１４から２５のいずれか一項に記載の装置。
【請求項２７】
プロセッサと、
プロセッサが実行可能な命令を記憶するメモリと、を含み、
前記プロセッサは、前記実行可能な命令を呼び出すことによって請求項１から１３のいずれか一項に記載の方法を実現することを特徴とする遺伝子変異認識装置。
【請求項２８】
コンピュータプログラム命令が記憶されている不揮発性コンピュータ可読記憶媒体であって、前記コンピュータプログラム命令がプロセッサにより実行されると、請求項１から１３のいずれか一項に記載の方法を実現することを特徴とする、不揮発性コンピュータ可読記憶媒体。
【発明の詳細な説明】
【関連出願の相互参照】
【０００１】
本開示は、２０１９年３月２９日に中国特許局に提出された、出願番号２０１９１０２５１８９１．０、出願の名称「遺伝子変異認識方法、装置および記憶媒体」の中国特許出願の優先権を主張し、その開示の全てが参照によって本開示に組み込まれる。
【技術分野】
【０００２】
本開示はコンピュータ技術の分野に関し、特に遺伝子変異認識方法、装置および記憶媒体に関する。
【背景技術】
【０００３】
バイオ技術の発展に伴い、ヒト遺伝子の配列は遺伝子シーケンシング技術によって測定できるようになり、遺伝子配列の解析は、遺伝子のさらなる研究および組み換えの基盤となり得る。現在、第２世代の遺伝子シーケンシング技術は、第１世代の遺伝子シーケンシング技術に比べて、遺伝子シーケンシングの効率が大幅に改善され、遺伝子シーケンシングのコストが低減されるとともに、遺伝子シーケンシングの正確性が保たれている。第１代のシーケンシング技術の場合、１つのヒトゲノムのシーケンシングを完了するまで３年の時間を要することがあるが、第２世代のシーケンシング技術を使用すると、時間を僅か１週間に短縮できる。
【発明の概要】
【０００４】
以上に鑑み、本開示は、遺伝子変異認識の解決手段を提供する。
【０００５】
本開示の一態様によれば、
変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得することと、
前記少なくとも１つの遺伝子シーケンシング断片の属性情報に基づいて、前記変異遺伝子座候補の、遺伝子座の位置に関連する配列特徴および非配列特徴を特定することと、
前記配列特徴および前記非配列特徴に基づいて、前記変異遺伝子座候補の遺伝子変異を認識することと、を含む遺伝子変異認識方法が提供される。
【０００６】
１つの可能な実施形態では、前記属性情報は、配列属性情報を含み、
前記少なくとも１つの遺伝子シーケンシング断片の属性情報に基づいて、前記変異遺伝子座候補の配列特徴を特定することは、
前記変異遺伝子座候補の遺伝子位置情報に基づいて、前記変異遺伝子座候補が位置する所定の遺伝子座区間を特定することと、
前記少なくとも１つの遺伝子シーケンシング断片の、前記所定の遺伝子座区間内に位置する各遺伝子座の、遺伝子座の位置に関連しかつ遺伝子属性を表す配列属性情報を取得することと、
前記所定の遺伝子座区間に位置する各遺伝子座の配列属性情報に基づいて、前記変異遺伝子座候補の配列特徴を生成することと、を含む。
【０００７】
１つの可能な実施形態では、
前記少なくとも１つの遺伝子シーケンシング断片の、前記所定の遺伝子座区間内に位置する各遺伝子座の配列属性情報を取得することは、
前記各遺伝子座において前記少なくとも１つの遺伝子シーケンシング断片の遺伝子型を特定することと、
前記各遺伝子座において遺伝子を各遺伝子型毎にカウントすることと、を含む。
【０００８】
１つの可能な実施形態では、
前記少なくとも１つの遺伝子シーケンシング断片の、前記所定の遺伝子座区間内に位置する各遺伝子座の配列属性情報を取得することは、
各遺伝子シーケンシング断片の遺伝子配列と参照ゲノム配列との比較結果に基づいて、前記各遺伝子座において各遺伝子シーケンシング断片の欠失遺伝子の遺伝子型を特定することと、
前記各遺伝子座において前記少なくとも１つの遺伝子シーケンシング断片の欠失遺伝子を各遺伝子型毎にカウントすることと、を含む。
【０００９】
１つの可能な実施形態では、
前記少なくとも１つの遺伝子シーケンシング断片の、前記所定の遺伝子座区間内に位置する各遺伝子座の配列属性情報を取得することは、
各遺伝子シーケンシング断片の遺伝子配列と参照ゲノム配列との比較結果に基づいて、前記各遺伝子座において各遺伝子シーケンシング断片の挿入遺伝子の遺伝子型を特定することと、
前記各遺伝子座において前記少なくとも１つの遺伝子シーケンシング断片の挿入遺伝子を各遺伝子型毎にカウントすることと、を含む。
【００１０】
１つの可能な実施形態では、
前記配列属性情報は、参照遺伝子の遺伝子型、各遺伝子型毎の遺伝子数、各遺伝子型毎の欠失遺伝子数、各遺伝子型毎の挿入遺伝子数の少なくとも１つを含む。
【００１１】
１つの可能な実施形態では、
前記属性情報は非配列属性情報を含み、
前記少なくとも１つの遺伝子シーケンシング断片の属性情報に基づいて、前記変異遺伝子座候補の非配列特徴を特定することは、
前記少なくとも１つの遺伝子シーケンシング断片の、遺伝子座の位置に関連しなくかつ遺伝子属性を表す非配列属性情報を取得することと、
前記少なくとも１つの遺伝子シーケンシング断片の非配列属性情報に基づいて、前記変異遺伝子座候補の非配列特徴を特定することと、を含む。
【００１２】
１つの可能な実施形態では、
前記非配列情報は、比較品質、プラス鎖およびマイナス鎖の偏り、遺伝子シーケンシング断片の長さ、エッジの偏りの少なくとも１つを含む。
【００１３】
１つの可能な実施形態では、
前記少なくとも１つの遺伝子シーケンシング断片の非配列属性情報に基づいて、前記変異遺伝子座候補の非配列特徴を特定することは、
各遺伝子シーケンシング断片の各遺伝子座の比較品質に基づいて、各遺伝子シーケンシング断片の比較品質を特定することであって、前記比較品質は、遺伝子シーケンシング断片の各遺伝子配列毎のシーケンシングの正確性を表すことと、
各遺伝子シーケンシング断片の比較品質に基づいて、前記変異遺伝子座候補に対応する非配列特徴を特定することと、を含む。
【００１４】
１つの可能な実施形態では、
前記少なくとも１つの遺伝子シーケンシング断片の非配列属性情報に基づいて、前記変異遺伝子座候補の非配列特徴を特定することは、
各遺伝子シーケンシング断片の属する遺伝子鎖がプラス鎖またはマイナス鎖のいずれであるかの情報に基づいて、前記少なくとも１つの遺伝子シーケンシング断片の、プラス鎖である遺伝子鎖とマイナス鎖である遺伝子鎖との割合を特定することと、
前記プラス鎖である遺伝子鎖とマイナス鎖である遺伝子鎖との割合に基づいて、前記変異遺伝子座候補に対応する非配列特徴を特定することと、を含む。
【００１５】
１つの可能な実施形態では、
前記配列特徴および前記非配列特徴に基づいて、前記変異遺伝子座候補の遺伝子変異を認識することは、
前記配列特徴および前記非配列特徴を統合し、前記変異遺伝子座候補の統合特徴を得ることと、
前記変異遺伝子座候補の統合特徴に基づいて、前記変異遺伝子座候補の遺伝子変異を認識することと、を含む。
【００１６】
１つの可能な実施形態では、
前記変異遺伝子座候補の統合特徴に基づいて、前記変異遺伝子座候補の遺伝子変異を認識することは、
前記変異遺伝子座候補の統合特徴に基づいて、前記変異遺伝子座候補の遺伝子の変異性を示す変異値を得ることと、
前記変異値が所定の閾値以上である場合、前記変異遺伝子座候補の遺伝子に変異があると判定することと、を含む。
【００１７】
１つの可能な実施形態では、
前記変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得することは、
体細胞遺伝子による遺伝子シーケンシングで得られた遺伝子シーケンシング断片を取得することと、
前記遺伝子シーケンシング断片の遺伝子配列を参照ゲノム配列と比較し、比較結果を得ることと、
前記比較結果に基づいて、前記体細胞遺伝子に異常がある変異遺伝子座候補を特定することと、
前記変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得することと、を含む。
【００１８】
本開示の別の態様によれば、
変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得するための取得モジュールと、
前記少なくとも１つの遺伝子シーケンシング断片の属性情報に基づいて、前記変異遺伝子座候補の、遺伝子座の位置に関連する配列特徴および非配列特徴を特定するための特定モジュールと、
前記配列特徴および前記非配列特徴に基づいて、前記変異遺伝子座候補の遺伝子変異を認識するための認識モジュールと、を含む遺伝子変異認識装置が提供される。
【００１９】
１つの可能な実施形態では、
前記属性情報は、配列属性情報を含み、
前記特定モジュールは、
前記変異遺伝子座候補の遺伝子位置情報に基づいて、前記変異遺伝子座候補が位置する所定の遺伝子座区間を特定する第１の特定サブモジュールと、
前記少なくとも１つの遺伝子シーケンシング断片の、前記所定の遺伝子座区間内に位置する各遺伝子座の、遺伝子座の位置に関連しかつ遺伝子属性を表す配列属性情報を取得する第１の取得サブモジュールと、
前記所定の遺伝子座区間に位置する各遺伝子座の配列属性情報に基づいて、前記変異遺伝子座候補の配列特徴を生成する第１の生成サブモジュールと、を含む。
【００２０】
１つの可能な実施形態では、
前記第１の取得サブモジュールは、具体的に、
前記各遺伝子座において前記少なくとも１つの遺伝子シーケンシング断片の遺伝子型を特定することと、
前記各遺伝子座において遺伝子を各遺伝子型毎にカウントすることと、に用いられる。
【００２１】
１つの可能な実施形態では、
前記第１の取得サブモジュールは、具体的に、
各遺伝子シーケンシング断片の遺伝子配列と参照ゲノム配列との比較結果に基づいて、前記各遺伝子座において各遺伝子シーケンシング断片の欠失遺伝子の遺伝子型を特定することと、
前記各遺伝子座において前記少なくとも１つの遺伝子シーケンシング断片の欠失遺伝子を各遺伝子型毎にカウントすることと、に用いられる。
【００２２】
１つの可能な実施形態では、
前記第１の取得サブモジュールは、具体的に、
各遺伝子シーケンシング断片の遺伝子配列と参照ゲノム配列との比較結果に基づいて、前記各遺伝子座において各遺伝子シーケンシング断片の挿入遺伝子の遺伝子型を特定することと、
前記各遺伝子座において前記少なくとも１つの遺伝子シーケンシング断片の挿入遺伝子を各遺伝子型毎にカウントすることと、に用いられる。
【００２３】
１つの可能な実施形態では、
前記配列属性情報は、参照遺伝子の遺伝子型、各遺伝子型毎の遺伝子数、各遺伝子型毎の欠失遺伝子数、各遺伝子型毎の挿入遺伝子数の少なくとも１つを含む。
【００２４】
１つの可能な実施形態では、
前記属性情報は非配列属性情報を含み、
前記特定モジュールは、
前記少なくとも１つの遺伝子シーケンシング断片の、遺伝子座の位置に関連しなくかつ遺伝子属性を表す非配列属性情報を取得する第２の取得サブモジュールと、
前記少なくとも１つの遺伝子シーケンシング断片の非配列属性情報に基づいて、前記変異遺伝子座候補の非配列特徴を特定する第２の特定サブモジュールと、を含む。
【００２５】
１つの可能な実施形態では、
前記非配列情報は、比較品質、プラス鎖およびマイナス鎖の偏り、遺伝子シーケンシング断片の長さ、エッジの偏りの少なくとも１つを含む。
【００２６】
１つの可能な実施形態では、
前記第２の特定サブモジュールは、具体的に、
各遺伝子シーケンシング断片の各遺伝子座の比較品質に基づいて、各遺伝子シーケンシング断片の比較品質を特定することであって、前記比較品質は、遺伝子シーケンシング断片の各遺伝子配列毎のシーケンシングの正確性を表すことと、
各遺伝子シーケンシング断片の比較品質に基づいて、前記変異遺伝子座候補に対応する非配列特徴を特定することとに用いられることを特徴とする。
【００２７】
１つの可能な実施形態では、
前記第２の特定サブモジュールは、具体的に、
各遺伝子シーケンシング断片の属する遺伝子鎖がプラス鎖またはマイナス鎖のいずれであるかの情報に基づいて、前記少なくとも１つの遺伝子シーケンシング断片の、プラス鎖である遺伝子鎖とマイナス鎖である遺伝子鎖との割合を特定することと、
前記プラス鎖である遺伝子鎖とマイナス鎖である遺伝子鎖との割合に基づいて、前記変異遺伝子座候補に対応する非配列特徴を特定することとに用いられることを特徴とする。
【００２８】
１つの可能な実施形態では、
前記認識モジュールは、
具体的に前記配列特徴および前記非配列特徴を統合し、前記変異遺伝子座候補の統合特徴を得る統合サブモジュールと、
前記変異遺伝子座候補の統合特徴に基づいて、前記変異遺伝子座候補の遺伝子変異を認識する認識サブモジュールと、を含む。
【００２９】
１つの可能な実施形態では、
前記認識サブモジュールは、具体的に、
前記変異遺伝子座候補の統合特徴に基づいて、前記変異遺伝子座候補の遺伝子の変異性を示す変異値を得ることと、
前記変異値が所定の閾値以上である場合、前記変異遺伝子座候補の遺伝子に変異があると判定することと、に用いられる。
【００３０】
１つの可能な実施形態では、
前記取得モジュールは、具体的に、
体細胞遺伝子による遺伝子シーケンシングで得られた遺伝子シーケンシング断片を取得することと、
前記遺伝子シーケンシング断片の遺伝子配列を参照ゲノム配列と比較し、比較結果を得ることと、
前記比較結果に基づいて、前記体細胞遺伝子に異常がある変異遺伝子座候補を特定することと、
前記変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得することと、に用いられる。
【００３１】
本開示の別の態様によれば、プロセッサと、プロセッサが実行可能な命令を記憶するメモリとを含み、前記プロセッサは上記方法を実行するように構成されている遺伝子変異認識装置が提供される。
【００３２】
本開示の別の態様によれば、コンピュータプログラム命令が記憶されている不揮発性コンピュータ可読記憶媒体であって、前記コンピュータプログラム命令がプロセッサにより実行されると、上記方法を実現する不揮発性コンピュータ可読記憶媒体が提供される。
【００３３】
本開示の実施例は、変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得し、少なくとも１つの遺伝子シーケンシング断片の属性情報に基づいて、変異遺伝子座候補の配列特徴および非配列特徴を特定することができ、これにより、特定された配列特徴および非配列特徴に基づいて、変異遺伝子座候補の遺伝子変異を認識することができる。ここで、配列特徴は、遺伝子座の位置に関連する特徴であってよく、非配列特徴は、遺伝子座の位置に関連しない特徴であってよく、これにより、遺伝子変異認識のプロセスにおいて、遺伝子の配列特徴および非配列特徴を組み合わせて、変異遺伝子座の特徴をより全面的に解析し、生殖細胞系遺伝子の変異、およびノイズやエラーによる干渉を除去し、遺伝子変異をより良好に認識し、遺伝子変異認識の正確性を高めることができる。
【００３４】
以下、図面を参照しながら例示的な実施例について詳細に説明することにより、本開示の他の特徴および態様は明瞭になる。
【図面の簡単な説明】
【００３５】
明細書に含まれ且つ明細書の一部を構成する図面は明細書と共に本開示の例示的な実施例、特徴および態様を示し、さらに本開示の原理を解釈するものである。
【図１】本開示の一実施例に係る遺伝子変異認識方法のフローチヤ一トを示す。
【図２】本開示の一実施例に係る変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得するプロセスのフローチヤ一トを示す。
【図３】本開示の一実施例に係る変異遺伝子座候補の配列特徴を特定するプロセスのフローチヤ一トを示す。
【図４】本開示の一実施例に係る変異遺伝子座候補の非配列特徴を特定するプロセスのフローチヤ一トを示す。
【図５】本開示の一実施例に係る変異遺伝子座候補の遺伝子変異を認識するプロセスのフローチヤ一トを示す。
【図６】本開示の一実施例に係るニューラルネットワークモデルのブロック図を示す。
【図７】本開示の一実施例に係る遺伝子変異認識装置のブロック図を示す。
【図８】本開示の一例示的実施例により示される遺伝子変異認識のための装置のブロック図を示す。
【発明を実施するための形態】
【００３６】
以下に図面を参照しながら本開示の様々な例示的実施例、特徴および態様を詳細に説明する。図面における同じ符号は同じまたは類似する機能の要素を表す。図面において実施例の様々な態様を示したが、特に断らない限り、比例で図面を作る必要がない。
【００３７】
ここの用語「例示的」とは、「例、実施例として用いられることまたは説明ためのなもの」を意味する。ここで「例示的」に説明するいかなる実施例も他の実施例より優れたものと理解すべきではない。
【００３８】
また、本開示をより効果的に説明するために、以下の具体的な実施形態において様々な具体的詳細を示す。当業者であれば、何らかの具体的詳細がなくなるにも関わらず、本開示は同様に実施できるということを理解すべきである。いくつかの実施例では、本開示の趣旨を強調するよう、当業者に既知の方法、手段、要素および回路に対する詳細な説明を省略する。
【００３９】
本開示の実施例により提供される遺伝子変異認識の解決手段は、変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得することができ、それにより少なくとも１つの遺伝子シーケンシング断片に基づいて、変異遺伝子座候補の遺伝子変異を認識することができる。遺伝子変異認識のプロセスにおいて、少なくとも１つの遺伝子シーケンシング断片の配列属性情報に基づいて配列特徴を生成し、少なくとも１つの遺伝子シーケンシング断片の非配列属性情報に基づいて非配列特徴を生成し、その後、配列特徴および非配列特徴によって変異遺伝子座候補の遺伝子変異を認識することができ、これにより、少なくとも１つの遺伝子シーケンシング断片の配列属性情報および非配列属性情報を統合し、遺伝子シーケンシング断片の配列属性情報をより全面的に利用することができる。
【００４０】
関連技術では、サポートベクターマシーンや、ランダムフォレストなど従来の機械学習方法を用いて遺伝子変異認識を行うことが一般的である。このような方法は、簡単に実現できるが、変異遺伝子座候補の近傍の遺伝子配列の配列属性情報の利用が困難であり、遺伝子のデータ量が一定の程度まで増加すると、遺伝子変異認識の効果がボトルネックに入る。また、ディープラーニングの方法を用いて、ニューラルネットワークによって遺伝子変異を認識する関連技術もある。しかし、ニューラルネットワークは、遺伝子配列の非配列情報を統合し難く、遺伝子データをより全面的に解析することはできない。本開示の実施例では、遺伝子変異認識において、マルチモーダル情報が統合されたニューラルネットワークモデルを用いて変異遺伝子座候補の配列特徴および非配列特徴を抽出することができ、これにより、遺伝子配列の配列属性情報および非配列属性情報を統合し、遺伝子データをより全面的に解析し、生殖細胞系遺伝子の変異、およびノイズやエラーによる干渉を除去し、遺伝子変異をより良好に認識することができる。以下、実施例により遺伝子変異の認識プロセスを詳細に説明する。
【００４１】
図１は本開示の一実施例に係る遺伝子変異認識方法のフローチヤ一トを示す。該遺伝子変異の認識方法は、遺伝子変異認識装置または他の処理機器により実行される。ここで、遺伝子変異認識装置は、ユーザ機器（Ｕｓｅｒ  Ｅｑｕｉｐｍｅｎｔ、ＵＥ）、モバイル機器、ユーザ端末、端末、セルラフォン、コードレス電話機、ＰＤＡ（Ｐｅｒｓｏｎａｌ  Ｄｉｇｉｔａｌ  Ａｓｓｉｓｔａｎｔ）、ハンドヘルド機器、コンピュータ機器、車載機器、ウェアラブル機器などでもよい。あるいは、遺伝子変異認識装置は、サーバでもよい。いくつかの可能な実施形態では、該遺伝子変異の認識方法は、プロセッサによってメモリに記憶されているコンピュータ読み取り可能な命令を呼び出すことにより実現される。
【００４２】
図１に示すように、該遺伝子変異の認識方法は以下のステップを含む。
【００４３】
ステップ１１、変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得する。
【００４４】
本開示の実施例では、遺伝子変異認識装置は、遺伝子シーケンシングにより得られた遺伝子シーケンシング断片を取得し、その後、遺伝子シーケンシングにより得られた遺伝子シーケンシング断片から、変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得することができる。ここで、遺伝子シーケンシング断片は、遺伝子シーケンシングを経て遺伝子型がラベリングされている遺伝子配列と解釈することができ、各遺伝子シーケンシング断片の長さは同一であってもよいし、異なっていてもよい。長さが異なる場合、各遺伝子シーケンシング断片の長さを所定の長さ範囲内とすることができ、これにより、各遺伝子シーケンシング断片の長さが比較的近似していることを保証できる。遺伝子型は、塩基タイプと解釈することができ、遺伝子型は、シトシン（Ｃ）、グアニン（Ｇ）、アデニン（Ａ）、チミジン（Ｔ）を含んでよく、これにより、遺伝子シーケンシング断片は、ＡＧＣＴを含む遺伝子配列であってもよい。ここで、変異遺伝子座候補は、遺伝子配列に異常がある遺伝子座であってもよい。遺伝子配列の遺伝子座は、遺伝子配列の位置を表すことができ、各遺伝子座には、少なくとも１つの遺伝子シーケンシング断片が存在し、すなわち、同じ遺伝子座には、遺伝子シーケンシングにより得られた少なくとも１つの遺伝子シーケンシング断片が存在してもよい。これに応じて、変異遺伝子座候補は少なくとも１つの遺伝子シーケンシング断片に対応し、ここで、この少なくとも１つの遺伝子シーケンシング断片は全てこの遺伝子座で異常が生じる。変異遺伝子座候補は少なくとも１つであってよく、各変異遺伝子座候補は、少なくとも１つの遺伝子シーケンシング断片に対応することが可能である。理解の便宜上、本開示の実施例は、１つの変異遺伝子座候補として説明する。
【００４５】
ステップ１２、前記少なくとも１つの遺伝子シーケンシング断片の属性情報に基づいて、前記変異遺伝子座候補の、遺伝子座の位置に関連する配列特徴および非配列特徴を特定する。
【００４６】
本開示の実施例では、変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得した後、該変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片の属性情報を抽出し、抽出された属性情報に基づいて、該変異遺伝子座候補の配列特徴および非配列特徴を生成することができる。属性情報は、配列属性情報および非配列属性情報を含んでもよい。配列属性情報は、遺伝子座位置に関連する、遺伝子シーケンシング断片の遺伝子属性を表す情報であってもよい。非配列属性情報は、遺伝子座の位置による制限を受けることがなく且つ遺伝子属性を表す情報であってもよい。属性情報を抽出する際に、該遺伝子座候補に対応する複数の遺伝子シーケンシング断片をランダムに選択し、ランダムに選択された複数の遺伝子シーケンシング断片の属性情報を抽出するようにしてもよく、また、該遺伝子座候補に対応する各遺伝子シーケンシング断片の属性情報を抽出してもよい。
【００４７】
ここで、配列属性情報を抽出する際に、該変異遺伝子座候補において少なくとも１つの遺伝子シーケンシング断片の配列属性情報を抽出してもよいし、該変異遺伝子座候補の近傍遺伝子座において少なくとも１つの遺伝子シーケンシング断片の配列属性情報を抽出してもよい。ここで、変異遺伝子座候補の配列特徴を特定する時、畳み込み層およびプーリング層を備えるニューラルネットワークモデルを利用して、変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片に対して、変異遺伝子座候補の配列特徴を抽出することができる。該ニューラルネットワークモデルは、２つの分岐構造を含み、１つの分岐は、畳み込み層およびプーリング層を含み、遺伝子シーケンシング断片の配列特徴を抽出するものである。もう１つの分岐は、遺伝子シーケンシング断片の非配列特徴を抽出するものである。これにより、該ニューラルネットワークモデルは複数種類のモーダル情報（配列属性情報および非配列属性情報）を統合し、変異遺伝子座候補の遺伝子変異を認識することができる。変異遺伝子座候補の非配列特徴を特定する時、上記ニューラルネットワークモデルを用いて、該ニューラルネットワークモデルのもう１つの分岐により少なくとも１つの遺伝子シーケンシング断片の非配列特徴を抽出することができ、該分岐構造は全結合層を含んでよく、全結合層は、位置による制限を受けることがない非配列特徴の抽出に使用され得る。
【００４８】
ステップ１３、前記配列特徴および前記非配列特徴に基づいて、前記変異遺伝子座候補の遺伝子変異を認識する。
【００４９】
本開示の実施形態では、変異遺伝子座候補の配列特徴および非配列特徴を特定した後、配列特徴および非配列特徴を融合し、該変異遺伝子座候補の遺伝子変異を認識することができ、例えば、上記ニューラルネットワークモデルを用いて、該変異遺伝子座候補の遺伝子変異の有無、または該変異遺伝子座候補の遺伝子の遺伝子配列異常はノイズなどに起因するかを判断することができる。
【００５０】
本開示の実施例では、変異遺伝子座候補の配列特徴および非配列特徴に基づいて変異遺伝子座候補の遺伝子変異を認識することができ、これにより、遺伝子シーケンシングデータをより全面的に解析することができる。変異遺伝子座候補の遺伝子変異を認識する時、まず変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得する必要がある。本開示の実施例では、変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得するプロセスをさらに提供する。
【００５１】
図２は、本開示の一実施例に係る変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得するプロセスのフローチヤ一トを示す。１つの可能な実施形態では、変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得することは、以下のステップを含み得る。
【００５２】
ステップ１１１、体細胞遺伝子による遺伝子シーケンシングで得られた遺伝子シーケンシング断片を取得する。
【００５３】
ここで、体細胞遺伝子による遺伝子シーケンシングにより、少なくとも１つの遺伝子シーケンシング断片を得ることができ、遺伝子シーケンシング断片は、体細胞遺伝子に遺伝子型がラベリングされた配列であってもよい。体細胞遺伝子を遺伝子シーケンシングして、遺伝子シーケンシング断片における各遺伝子の遺伝子型を得ることができるだけでなく、さらに遺伝子シーケンシング断片における各遺伝子が位置する遺伝子座の遺伝子位置情報を得ることもできる。同じ遺伝子座は少なくとも１つの遺伝子シーケンシング断片に対応し得る。
【００５４】
１つの可能な実施形態では、体細胞遺伝子に対して遺伝子シーケンシングを行うことで少なくとも１つの遺伝子シーケンシング断片を取得することができ、遺伝子シーケンシングで得られた遺伝子シーケンシング断片を予備処理することができ、ここの予備処理方法は、交差汚染のスクリーニング、シーケンシング品質のスクリーニング、比較品質のスクリーニング、リード長さ異常のスクリーニングなどを含み得る。予備処理によって、交差汚染がある遺伝子シーケンシング断片、およびシーケンシング品質や比較品質が低く、リード長さに異常がある遺伝子シーケンシング断片を除外することができる。
【００５５】
ステップ１１２、前記遺伝子シーケンシング断片の遺伝子配列を参照ゲノム配列と比較し、比較結果を得る。
【００５６】
本開示の実施例では、体細胞遺伝子による遺伝子シーケンシングで得られた遺伝子シーケンシング断片を取得した後、取得された遺伝子シーケンシング断片の遺伝子配列を同じ遺伝子座の参照ゲノム配列と比較し、比較結果を得ることができる。例えば、遺伝子シーケンシングを行って得られた各遺伝子シーケンシング断片を同じ遺伝子座の参照ゲノム配列と比較し、遺伝子シーケンシング断片の遺伝子配列が参照ゲノム配列と異なる遺伝子座を特定することができる。同じ遺伝子座にある少なくとも１つの遺伝子シーケンシング断片を同じ遺伝子座の参照ゲノム配列と比較し、少なくとも１つの遺伝子シーケンシング断片の遺伝子配列が参照ゲノム配列と異なる遺伝子座を特定することもできる。
【００５７】
ステップ１１３、前記比較結果に基づいて、前記体細胞遺伝子に異常がある変異遺伝子座候補を特定する。
【００５８】
本開示の実施例では、比較結果に基づいて、遺伝子シーケンシング断片の遺伝子配列が参照ゲノム配列と異なる遺伝子座を特定することができ、該遺伝子座に対応する少なくとも１つの遺伝子シーケンシング断片に対して、該遺伝子座で変異が生じる遺伝子シーケンシング断片の占める割合が所定の割合より大きい場合、該遺伝子座は変異遺伝子座候補であると判断することができ、そうでない場合、該遺伝子座は変異遺伝子座候補ではないと考えてよい。該遺伝子座における遺伝子シーケンシング断片の遺伝子配列と参照ゲノム配列との違いは、シーケンシングエラーによるものである可能性があるため、このようにすれば、遺伝子シーケンシングの誤りによる遺伝子配列異常の現象を低減することができる。
【００５９】
ステップ１１４、前記変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得する。
【００６０】
本開示の実施例では、変異遺伝子座候補を特定した後、変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得することができる。ここで、各変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片における該変異遺伝子座候補の遺伝子配列は、同じ遺伝子座の参照ゲノム配列と異なってもよい。ここで、変異遺伝子座候補は少なくとも１つであってもよい。
【００６１】
上記の変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得するプロセスによって、変異遺伝子座候補を正確に特定することができるだけでなく、遺伝子シーケンシングで得られた遺伝子シーケンシング断片から、変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を特定することもできる。
【００６２】
本開示の実施例では、変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片の配列属性情報に基づいて、該変異遺伝子座候補の配列特徴を特定することができ、これにより、変異遺伝子座候補の遺伝子変異を認識する時、変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片の配列属性を考慮することができる。以下、一例により変異遺伝子座候補の配列特徴を特定するプロセスを詳細に説明する。
【００６３】
図３は、本開示の一実施例に係る変異遺伝子座候補の配列特徴を特定するプロセスのフローチヤ一トを示す。図３に示すように、上記ステップ１２は、
前記変異遺伝子座候補の遺伝子位置情報に基づいて、前記変異遺伝子座候補が位置する所定の遺伝子座区間を特定するステップ１２１ａと、
前記少なくとも１つの遺伝子シーケンシング断片の、前記所定の遺伝子座区間内に位置する各遺伝子座の、遺伝子座の位置に関連しかつ遺伝子属性を表す配列属性情報を取得するステップ１２２ａと、
前記所定の遺伝子座区間に位置する各遺伝子座の配列属性情報に基づいて、前記変異遺伝子座候補の配列特徴を生成するステップ１２３ａと、を含んでもよい。
【００６４】
本開示の実施例の一例では、各変異遺伝子座候補について、少なくとも１つの遺伝子シーケンシング断片が存在してもよい。遺伝子変異認識の正確性を高めるには、該変異遺伝子座候補の配列属性情報に加えて、さらに該変異遺伝子座候補の近傍の遺伝子座の配列属性情報も考慮することができる。変異遺伝子座候補の配列特徴を特定する時、変異遺伝子座候補の遺伝子位置情報に基づいて、該変異遺伝子座候補が位置する所定の遺伝子座区間を特定することができ、例えば、変異遺伝子座候補の前後の１５０個の塩基対の区間を変異遺伝子座候補が位置する所定の遺伝子座区間とする。その後、該所定の遺伝子座区間に位置する各遺伝子座について、該遺伝子座において少なくとも１つの遺伝子シーケンシング断片の配列属性情報を取得し、該遺伝子座の配列属性情報から該遺伝子座に対応する配列特徴を生成することが可能である。配列特徴は、配列特徴ベクトルで表すことができる。変異遺伝子座候補が位置する所定の遺伝子座区間内の少なくとも１つの遺伝子座に対応する少なくとも１つの配列特徴ベクトルから、変異遺伝子座候補の配列特徴行列を生成することができる。例えば、変異遺伝子座候補が位置する所定の遺伝子座区間は３つの遺伝子座ｂ１、ｂ２、ｂ３を含み、３つの遺伝子座に対応する配列特徴ベクトルはそれぞれａ１、ａ２、ａ３である場合、変異遺伝子座候補の配列特徴行列は［ａ１  ａ２  ａ３］となり、ここで、ａ１、ａ２、ａ３の配列特徴は、ｂ１、ｂ２、ｂ３の配列属性情報にそれぞれ対応する。
【００６５】
ここで、配列属性情報は、参照ゲノムの遺伝子型、各遺伝子型毎の遺伝子数、各遺伝子型毎の欠失遺伝子数、各遺伝子型毎の挿入遺伝子数を含んでもよいが、これらに限定されない。参照ゲノムの遺伝子型は、変異遺伝子座候補における参照ゲノムの遺伝子型であってもよい。各遺伝子型毎の遺伝子数は、該変異遺伝子座候補における少なくとも１つの遺伝子シーケンシング断片の各遺伝子型毎の遺伝子数であってもよい。例えば、該変異遺伝子座候補は５つの遺伝子シーケンシング断片に対応し、該変異遺伝子座候補における各遺伝子シーケンシング断片の遺伝子型はそれぞれＡ、Ｃ、Ｃ、Ｇ、Ｇである場合、各遺伝子型毎の遺伝子数はそれぞれ、Ａが１つ、Ｃが２つ、Ｇが２つである。各遺伝子型毎の欠失遺伝子数は、該変異遺伝子座候補における少なくとも１つの遺伝子シーケンシング断片の各遺伝子型毎の欠失遺伝子数であってもよい。例えば、該変異遺伝子座候補における各遺伝子シーケンシング断片の欠失の遺伝子型はそれぞれＡ、Ｃ、Ｃ、Ｇ、Ｇである場合、各遺伝子型毎の欠失遺伝子数はそれぞれ、Ａが１つ、Ｃが２つ、Ｇが２つである。各遺伝子型毎の挿入遺伝子数は、該変異遺伝子座候補における少なくとも１つの遺伝子シーケンシング断片の各遺伝子型毎の挿入遺伝子数であってよもい。例えば、該変異遺伝子座候補における各遺伝子シーケンシング断片の挿入の遺伝子型はそれぞれＡ、Ｃ、Ｃ、Ｇ、Ｇである場合、各遺伝子型毎の挿入遺伝子数はそれぞれＡが１つ、Ｃが２つ、Ｇが２つである。
【００６６】
１つの可能な実施形態では、所定の遺伝子座区間内の各遺伝子座において少なくとも１つの遺伝子シーケンシング断片の配列属性情報を取得する時、該所定の遺伝子座区間内の各遺伝子座において、少なくとも１つの遺伝子シーケンシング断片の遺伝子型を特定し、該遺伝子座に対応する遺伝子を各遺伝子型毎にカウントすることができ、これにより、変異遺伝子座候補において、該遺伝子座に対応する少なくとも１つの遺伝子シーケンシング断片の各遺伝子型毎の遺伝子数を特定することができる。
【００６７】
１つの可能な実施形態では、所定の遺伝子座区間内の各遺伝子座において少なくとも１つの遺伝子シーケンシング断片の配列属性情報を取得する時、各遺伝子シーケンシング断片の遺伝子配列と参照ゲノム配列との比較結果に基づき、該所定の遺伝子座区間内の各遺伝子座において、各遺伝子シーケンシング断片の欠失遺伝子の遺伝子型を特定し、該遺伝子座において少なくとも１つの遺伝子シーケンシング断片の欠失遺伝子を各遺伝子型毎にカウントし、これにより、変異遺伝子座候補において、該遺伝子座に対応する少なくとも１つの遺伝子シーケンシング断片の各遺伝子型毎の欠失遺伝子の数を特定することができる。
【００６８】
１つの可能な実施形態では、所定の遺伝子座区間内の各遺伝子座において少なくとも１つの遺伝子シーケンシング断片の配列属性情報を取得する時、各遺伝子シーケンシング断片の遺伝子配列と参照ゲノム配列との比較結果に基づき、該所定の遺伝子座区間内の各遺伝子座において各遺伝子シーケンシング断片の欠失遺伝子の遺伝子型を特定し、該遺伝子座において、少なくとも１つの遺伝子シーケンシング断片のの挿入遺伝子を各遺伝子型毎にカウントし、これにより、変異遺伝子座候補において、該遺伝子座に対応する少なくとも１つの遺伝子シーケンシング断片の各遺伝子型毎の挿入遺伝子数を特定することができる。
【００６９】
例えば、配列属性情報は参照ゲノムの遺伝子型、各遺伝子型毎の遺伝子数、各遺伝子型毎の欠失遺伝子数、各遺伝子型毎の挿入遺伝子数を含むと仮定する場合、変異遺伝子座候補の配列特徴を特定する時、変異遺伝子座候補が位置する所定の遺伝子座区間内の各遺伝子座において、変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片の上記４つの情報を抽出することができ、例えば、変異遺伝子座候補に対応する５つの遺伝子シーケンシング断片の場合、所定の遺伝子座区間内のある遺伝子座について、参照ゲノムの遺伝子型、５つの遺伝子シーケンシング断片の各遺伝子型毎の遺伝子数、５つの遺伝子シーケンシング断片の各遺伝子型毎の欠失遺伝子数、および５つの遺伝子シーケンシング断片の各遺伝子型毎の挿入遺伝子数をそれぞれ特定することができる。その後、該遺伝子座に対応する少なくとも１つの配列属性情報を統合することで、該遺伝子座の配列特徴を得ることができる。変異遺伝子座候補の配列特徴は、所定の遺伝子座区間内の各遺伝子座の配列特徴を含んでもよい。
【００７０】
本開示の実施例の一例では、変異遺伝子座候補の遺伝子変異を認識する時、変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片の配列属性だけでなく、さらに少なくとも１つの遺伝子シーケンシング断片の非配列属性も考慮した。以下、一例により変異遺伝子座候補の非配列特徴を特定するプロセスを詳細に説明する。
【００７１】
図４は本開示の一実施例に係る変異遺伝子座候補の非配列特徴を特定するプロセスのフローチヤ一トを示す。図４に示すように、上記ステップ１２は、
前記少なくとも１つの遺伝子シーケンシング断片の、遺伝子座の位置に関連しなくかつ遺伝子属性を表す非配列属性情報を取得するステップ１２１ｂと、
前記少なくとも１つの遺伝子シーケンシング断片の非配列属性情報に基づいて、前記変異遺伝子座候補の非配列特徴を生成するステップ１２２ｂと、を含んでもよい。
【００７２】
本開示の実施例の一例では、遺伝子変異認識の正確性を高めるには、少なくとも１つの遺伝子シーケンシング断片の配列属性情報に加えて、さらに少なくとも１つの遺伝子シーケンシング断片の非配列属性情報も考慮することができる。ここで、非配列情報は、比較品質、プラス鎖およびマイナス鎖の偏り、遺伝子シーケンシング断片の長さ、エッジの偏りの少なくとも１つを含んでもよい。変異遺伝子座候補の非配列特徴を特定する時、少なくとも１つの遺伝子属性配列リードの非配列属性情報を取得し、その後、取得された非配列属性情報から変異遺伝子座候補の非配列特徴を生成することができる。
【００７３】
１つの可能な実施形態では、前記少なくとも１つの遺伝子シーケンシング断片の非配列属性情報に基づいて、前記変異遺伝子座候補の非配列特徴を特定する時、各遺伝子シーケンシング断片内の各遺伝子座の比較品質に基づいて、各遺伝子シーケンシング断片の比較品質を特定し、その後、各遺伝子シーケンシング断片の比較品質に基づいて、前記変異遺伝子座候補に対応する非配列特徴を特定することができる。ここで、比較品質は、遺伝子シーケンシング断片の各遺伝子配列の遺伝子シーケンシングの正確性を表すことができる。ある遺伝子配列の比較品質が所定の値より低い場合、該遺伝子配列の遺伝子シーケンシングによる遺伝子型が正確でないと考えられ、これにより、比較品質は変異遺伝子座候補の遺伝子変異の有無を判断するための１つの参考要素となり得る。例えば、変異遺伝子座候補は少なくとも１つの遺伝子シーケンシング断片に対応する場合、各遺伝子配列の比較品質に基づいて、各遺伝子シーケンシング断片の比較品質を特定することができ、１つの遺伝子シーケンシング断片を例として、該遺伝子シーケンシング断片に含まれる遺伝子配列の比較品質の平均値または中間値を該遺伝子シーケンシング断片の比較品質としてもよく、該遺伝子シーケンシング断片から少なくとも１つの遺伝子配列をランダムに選択し、選択された少なくとも１つの遺伝子配列の比較品質の平均値または中間値を該遺伝子シーケンシング断片の比較品質としてもよい。その後、各遺伝子シーケンシング断片の比較品質から該変異遺伝子座候補に対応する比較品質を取得し、例えば、該変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片の比較品質の平均値または中間値を計算し、該変異遺伝子座候補に対応する比較品質を取得し、これにより、該変異遺伝子座候補に対応する比較品質に基づいて変異遺伝子座候補に対応する非配列特徴を特定することができる。
【００７４】
１つの可能な実施形態では、少なくとも１つの遺伝子シーケンシング断片の非配列属性情報に基づいて、変異遺伝子座候補の非配列特徴を特定する時、各遺伝子シーケンシング断片の属する遺伝子鎖がプラス鎖またはマイナス鎖のいずれであるかの情報に基づいて、少なくとも１つの遺伝子シーケンシング断片の、プラス鎖である遺伝子鎖とマイナス鎖である遺伝子鎖との割合を特定し、その後、特定されたプラス鎖である遺伝子鎖とマイナス鎖である遺伝子鎖との割合に基づいて、変異遺伝子座候補に対応する非配列特徴を特定することができる。ここで、プラス鎖およびマイナス鎖の偏りは、遺伝子シーケンシング断片の属する遺伝子鎖においてプラス鎖とマイナス鎖との割合であってよく、遺伝子鎖は、プラス鎖およびマイナス鎖を含み、ここで、プラス鎖は、リボ核酸（ＲＮＡ）の塩基配列と同じであるデオキシリボ核酸（ＤＮＡ）単一鎖であってよく、マイナス鎖は、リボ核酸（ＲＮＡ）の塩基配列に相補的なデオキシリボ核酸（ＤＮＡ）単一鎖であってよい。例えば、変異遺伝子座候補は５つの遺伝子シーケンシング断片に対応し、そのうち、３つの遺伝子シーケンシング断片は遺伝子鎖のプラス鎖に対応し、２つの遺伝子シーケンシング断片は遺伝子鎖のマイナス鎖に対応する場合、プラス鎖およびマイナス鎖の偏りは３：２である。
【００７５】
１つの可能な実施形態では、少なくとも１つの遺伝子シーケンシング断片の非配列属性情報に基づいて、変異遺伝子座候補の非配列特徴を特定する時、各遺伝子シーケンシング断片の遺伝子シーケンシング断片長さに基づいて、変異遺伝子座候補の非配列特徴を特定することができる。遺伝子シーケンシング断片の長さは、各遺伝子シーケンシング断片が有する塩基配列の長さであってもよく、例えば、１つの遺伝子シーケンシング断片に４つの塩基配列が含まれる場合、該遺伝子シーケンシング断片の長さは４であり、各遺伝子シーケンシング断片の長さから変異遺伝子座候補の非配列特徴を特定してもよく、少なくとも１つの遺伝子シーケンシング断片の長さの中間値または平均値から変異遺伝子座候補の非配列特徴を特定してもよい。
【００７６】
１つの可能な実施形態では、少なくとも１つの遺伝子シーケンシング断片の非配列属性情報に基づいて、変異遺伝子座候補の非配列特徴を特定する時、各遺伝子シーケンシング断片のエッジの偏りに基づいて、変異遺伝子座候補の非配列特徴を特定することができる。ここで、エッジの偏りは、ある遺伝子座について、当該遺伝子座がエッジ位置に位置する遺伝子シーケンシング断片と当該遺伝子座が中間位置に位置する遺伝子シーケンシング断片との割合であってよい。例えば、遺伝子シーケンシング断片を３つのセグメントに均一に分けて、そのうち、遺伝子シーケンシング断片両端の２セグメントをエッジ位置とし、遺伝子シーケンシング断片中央の１セグメントを中間位置とし、変異遺伝子座候補が５つの遺伝子シーケンシング断片に対応する場合、変異遺伝子座候補は、３つの遺伝子シーケンシング断片においてエッジ位置に位置し、２つの遺伝子シーケンシング断片において中間位置に位置すれば、該変異遺伝子座候補のエッジの偏りは３：２となる。これに応じて、変異遺伝子座候補の各遺伝子シーケンシング断片でのエッジの偏りから、変異遺伝子座候補の非配列特徴を特定してもよく、少なくとも１つの遺伝子シーケンシング断片に対応するエッジの偏りの中間値または平均値から、変異遺伝子座候補の非配列特徴を特定してもよい。
【００７７】
以上のようにして、変異遺伝子座候補における少なくとも１つの遺伝子シーケンシング断片の非配列属性情報から変異遺伝子座候補の非配列特徴を生成することができ、遺伝子変異認識時に変異遺伝子座候補の非配列特徴を考慮することができる。これにより、遺伝子変異認識がより正確になる。非配列特徴を特定する時、非配列属性情報のいずれか１つ以上の情報の組み合わせから、少なくとも１つの遺伝子シーケンシング断片の非配列特徴を生成することができる。
【００７８】
以下、一例により変異遺伝子座候補の遺伝子変異を認識するプロセスを説明する。
【００７９】
図５は、本開示の一実施例に係る変異遺伝子座候補の遺伝子変異を認識するプロセスのフローチヤ一トを示す。図５に示すように、上記ステップ１３は、
前記配列特徴および前記非配列特徴を統合し、前記変異遺伝子座候補の統合特徴を得るステップ１３１と、
前記変異遺伝子座候補の統合特徴に基づいて、前記変異遺伝子座候補の遺伝子変異を認識するステップ１３２と、を含み得る。
【００８０】
本開示の実施例では、変異遺伝子座候補の配列特徴および非配列次元特徴を特定した後、ニューラルネットワークモデルを用いて、配列特徴および非配列特徴を統合し、配列特徴から形成された配列特徴行列と、非配列特徴から形成された非配列特徴行列とを１つの特徴行列として合成し、特徴が統合された統合特徴行列を取得し、その後、ニューラルネットワークモデルを用いて、該統合特徴行列に基づいて変異候補遺伝子座での遺伝子変異を認識することができる。このようにして、ニューラルネットワークモデルを用いて、変異遺伝子座候補に対応する配列属性情報および非配列属性情報を統合することができ、これにより、遺伝子シーケンシングデータをより全面的に解析することができ、遺伝子変異の認識はより正確になる。トレーニング中において、トレーニングサンプルとして、一塩基多型（Ｓｉｎｇｌｅ  Ｎｕｃｌｅｏｔｉｄｅ  Ｐｏｌｙｍｏｒｐｈｉｓｍ、ＳＮＰ）のある遺伝子シーケンシング断片、挿入／欠失（Ｉｎｓｅｒｔｉｏｎ／Ｄｅｌｅｔｉｏｎ、ＩｎＤｅｌ）のある遺伝子シーケンシング断片を選択することができ、これにより、トレーニング後に得られた遺伝子変異認識モデルはＳＮＰ、ＩｎＤｅｌの遺伝子変異を効果的に認識することができる。
【００８１】
１つの可能な実施形態では、前記変異遺伝子座候補の統合特徴に基づいて、前記変異遺伝子座候補の遺伝子変異を認識することは、前記変異遺伝子座候補の統合特徴に基づいて、前記変異遺伝子座候補の遺伝子の変異性を示す変異値を得ることと、前記変異値が所定の閾値以上である場合、前記変異遺伝子座候補の遺伝子に変異があると判定することと、を含んでもよい。ここで、遺伝子の変異性を示す変異値は、該変異遺伝子座候補で変異が生じる可能性を表すものであってもよい。例えば、変異値が大きいほど、該変異遺伝子座候補で変異が生じる可能性は大きい。上記ニューラルネットワークを用いて２次元特徴を処理して変異値を得て、変異値に基づいて変異遺伝子座候補の遺伝子に変異があるか否かを判断することができる。１つの可能な実施形態では、変異値を０から１の範囲内とすることができる。所定の閾値は、適用場面に応じて、例えば０．３、０．５に設定することができ、変異値が所定の閾値より大きい場合、該変異遺伝子座候補の遺伝子に変異が生じると考えてよく、そうでない場合、該変異遺伝子座候補の遺伝子に変異がないと考えてよい。
【００８２】
本開示の実施例では、ニューラルネットワークモデルを用いて変異遺伝子座候補の遺伝子変異を認識することができ、該ニューラルネットワークモデルは、変異遺伝子座候補の配列特徴および非配列特徴を抽出することができる。本開示の実施例では、ニューラルネットワークモデルの構造をさらに提供する。
【００８３】
図６は、本開示の一実施例に係るニューラルネットワークモデルのブロック図を示す。図６に示すように、ニューラルネットワークモデルは第１の分岐および第２の分岐の２つの分岐構造を含んでよもい。第１の分岐は、畳み込み層およびプーリング層を含み、変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片の配列特徴の抽出に使用される。第２の分岐は、全結合層を含み、変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片の非配列特徴の抽出に使用される。ニューラルネットワークモデルは、変異遺伝子座候補の配列特徴および非配列特徴を抽出した後、配列特徴および非配列特徴を統合し、例えば、配列特徴の配列特徴行列と非配列特徴の非配列特徴行列とを結合させ、特徴が統合された統合特徴行列を得ることができ、その後、全結合層により変異遺伝子座候補の変異値を得ることができる。
【００８４】
本開示の実施例では、変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片の配列属性情報および非配列属性情報を抽出し、配列属性情報および非配列属性情報を統合することで得られた統合特徴を用いて、遺伝子変異を認識し、これにより、変異遺伝子座候補に対応する配列属性情報および非配列属性情報を総合的に考慮し、遺伝子シーケンシング情報をより全面的に解析し、遺伝子候補遺伝子座での遺伝子変異をより良好に認識し、生殖細胞系遺伝子の変異、およびノイズやエラーによる干渉を除去し、遺伝子変異認識の正確率は高まる。
【００８５】
具体的な実施形態の上記方法において、各ステップの記述順序は厳しい実行順序であるというわけではなく、実施プロセスの何の制限にもならず、各ステップの具体的な実行順序はその機能と可能な内在的論理に依存することが当業者に理解される。
【００８６】
図７は、本開示の実施例に係る遺伝子変異認識装置のブロック図を示し、図７に示すように、前記遺伝子変異認識装置は、
変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得する取得モジュール７１と、
前記少なくとも１つの遺伝子シーケンシング断片の属性情報に基づいて、前記変異遺伝子座候補の、遺伝子座の位置に関連する配列特徴および非配列特徴を特定する特定モジュール７２と、
前記配列特徴および前記非配列特徴に基づいて、前記変異遺伝子座候補の遺伝子変異を認識するための認識モジュール７３と、を含む。
【００８７】
１つの可能な実施形態では、
前記属性情報は、配列属性情報を含み、
前記特定モジュール７２は、
前記変異遺伝子座候補の遺伝子位置情報に基づいて、前記変異遺伝子座候補が位置する所定の遺伝子座区間を特定する第１の特定サブモジュールと、
前記少なくとも１つの遺伝子シーケンシング断片の、前記所定の遺伝子座区間内に位置する各遺伝子座の、遺伝子座の位置に関連しかつ遺伝子属性を表す配列属性情報を取得する第１の取得サブモジュールと、
前記所定の遺伝子座区間に位置する各遺伝子座の配列属性情報に基づいて、前記変異遺伝子座候補の配列特徴を生成する第１の生成サブモジュールと、を含む。
【００８８】
１つの可能な実施形態では、
前記第１の取得サブモジュールは、具体的に、
前記各遺伝子座において前記少なくとも１つの遺伝子シーケンシング断片の遺伝子型を特定することと、
前記各遺伝子座において遺伝子を各遺伝子型毎にカウントすることと、に用いられる。
【００８９】
１つの可能な実施形態では、
前記第１の取得サブモジュールは、具体的に、
各遺伝子シーケンシング断片の遺伝子配列と参照ゲノム配列との比較結果に基づいて、前記各遺伝子座において各遺伝子シーケンシング断片の欠失遺伝子の遺伝子型を特定することと、
前記各遺伝子座において前記少なくとも１つの遺伝子シーケンシング断片の欠失遺伝子を各遺伝子型毎にカウントすることと、に用いられる。
【００９０】
１つの可能な実施形態では、
前記第１の取得サブモジュールは、具体的に、
各遺伝子シーケンシング断片の遺伝子配列と参照ゲノム配列との比較結果に基づいて、前記各遺伝子座において各遺伝子シーケンシング断片の挿入遺伝子の遺伝子型を特定することと、
前記各遺伝子座において前記少なくとも１つの遺伝子シーケンシング断片の挿入遺伝子を各遺伝子型毎にカウントすることと、に用いられる。
【００９１】
１つの可能な実施形態では、
前記配列属性情報は、参照遺伝子の遺伝子型、各遺伝子型毎の遺伝子数、各遺伝子型毎の欠失遺伝子数、各遺伝子型毎の挿入遺伝子数の少なくとも１つを含む。
【００９２】
１つの可能な実施形態では、
前記属性情報は、非配列属性情報を含み、
前記特定モジュールは、
前記少なくとも１つの遺伝子シーケンシング断片の、遺伝子座の位置に関連しなくかつ遺伝子属性を表す非配列属性情報を取得する第２の取得サブモジュールと、
前記少なくとも１つの遺伝子シーケンシング断片の非配列属性情報に基づいて、前記変異遺伝子座候補の非配列特徴を特定する第２の特定サブモジュールと、を含む。
【００９３】
１つの可能な実施形態では、
前記非配列情報は、比較品質、プラス鎖およびマイナス鎖の偏り、遺伝子シーケンシング断片の長さ、エッジの偏りの少なくとも１つを含む。
【００９４】
１つの可能な実施形態では、
前記第２の特定サブモジュールは、具体的に、
各遺伝子シーケンシング断片の各遺伝子座の比較品質に基づいて、各遺伝子シーケンシング断片の比較品質を特定することであって、前記比較品質は、遺伝子シーケンシング断片の各遺伝子配列毎のシーケンシングの正確性を表すことと、
各遺伝子シーケンシング断片の比較品質に基づいて、前記変異遺伝子座候補に対応する非配列特徴を特定することとに用いられる。
【００９５】
１つの可能な実施形態では、
前記第２の特定サブモジュールは、具体的に、
各遺伝子シーケンシング断片の属する遺伝子鎖がプラス鎖またはマイナス鎖のいずれであるかの情報に基づいて、前記少なくとも１つの遺伝子シーケンシング断片の、プラス鎖である遺伝子鎖とマイナス鎖である遺伝子鎖との割合を特定することと、
前記プラス鎖である遺伝子鎖とマイナス鎖である遺伝子鎖との割合に基づいて、前記変異遺伝子座候補に対応する非配列特徴を特定することとに用いられる。
【００９６】
１つの可能な実施形態では、
前記認識モジュール７３は、
具体的に前記配列特徴および前記非配列特徴を統合し、前記変異遺伝子座候補の統合特徴を得る統合サブモジュールと、
前記変異遺伝子座候補の統合特徴に基づいて、前記変異遺伝子座候補の遺伝子変異を認識する認識サブモジュールと、を含む。
【００９７】
１つの可能な実施形態では、
前記認識サブモジュールは、具体的に、
前記変異遺伝子座候補の統合特徴に基づいて、前記変異遺伝子座候補の遺伝子の変異性を示す変異値を得ることと、
前記変異値が所定の閾値以上である場合、前記変異遺伝子座候補の遺伝子に変異があると判定することと、に用いられる。
【００９８】
１つの可能な実施形態では、
前記取得モジュール７１は、具体的に、
体細胞遺伝子による遺伝子シーケンシングで得られた遺伝子シーケンシング断片を取得することと、
前記遺伝子シーケンシング断片の遺伝子配列を参照ゲノム配列と比較し、比較結果を得ることと、
前記比較結果に基づいて、前記体細胞遺伝子に異常がある変異遺伝子座候補を特定することと、
前記変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得することと、に用いられる。
【００９９】
いくつかの実施例では、本開示の実施例で提供された装置に備えた機能またはモジュールは、上記方法実施例に記載の方法を実行するものであって、その具現化については上記方法実施例の説明を参照してもよく、簡単化するために、ここで重複説明は割愛する。
【０１００】
図８は一例示的実施例により示される遺伝子変異認識のための装置１９００のブロック図である。例えば、装置１９００はサーバとして提供できる。図８を参照して、装置１９００は、さらに一つ又は複数のプロセッサを含む処理コンポーネント１９２２、および、処理コンポーネント１９２２が実行可能な命令、例えばアプリケーションプログラムを記憶する、メモリ１９３２を代表とするメモリ資源を含む。メモリ１９３２に記憶されたアプリケーションプログラムはそれぞれが１グループの命令に対応する一つ以上のモジュールを含むことができる。また、処理コンポーネント１９２２は命令を実行することで上記方法を実行するように構成される。
【０１０１】
装置１９００はさらに、装置１９００の電源管理を実行するように構成された一つの電源コンポーネント１９２６、装置１９００をネットワークに接続するように構成された一つの有線または無線ネットワークインタフェース１９５０、および一つの入出力（Ｉ／Ｏ）インタフェース１９５８を含むことができる。装置１９００はメモリ１９３２に記憶されたオペレーティングシステム、例えばＷｉｎｄｏｗｓ  ＳｅｒｖｅｒＴＭ、Ｍａｃ  ＯＳ  ＸＴＭ、ＵｎｉｘＴＭ、ＬｉｎｕｘＴＭ、ＦｒｅｅＢＳＤＴＭまたは類似するものに基づいて動作できる。
【０１０２】
例示的な実施例では、さらに、非揮発性のコンピュータが読み取り可能な記憶媒体、例えば、装置１９００の処理コンポーネント１９２２によって実行されることで上記方法を実行するためのコンピュータプログラム命令が記憶されているメモリ１９３２が提供される。
【０１０３】
本開示はシステム、方法および／またはコンピュータプログラム製品であってもよい。コンピュータプログラム製品は、プロセッサに本開示の各態様を実現させるためのコンピュータ可読プログラム命令がロードされているコンピュータが読み取り可能な記憶媒体を含むことができる。
【０１０４】
コンピュータが読み取り可能な記憶媒体は、命令実行機器により使用される命令を保持および記憶可能な有形機器であってもよい。コンピュータが読み取り可能な記憶媒体は、例えば、電気記憶装置、磁気記憶装置、光記憶装置、電磁記憶装置、半導体記憶装置または上記の任意の適当な組み合わせであってもよいが、これらに限定されない。コンピュータが読み取り可能な記憶媒体のさらに具体的な例（非包括的リスト）としては、携帯型コンピュータディスク、ハードディスク、ランダムアクセスメモリ（ＲＡＭ）、読み取り専用メモリ（ＲＯＭ）、消去可能プログラマブル読み取り専用メモリ（ＥＰＲＯＭまたはフラッシュメモリ）、静的ランダムアクセスメモリ（ＳＲＡＭ）、携帯型コンパクトディスク読み取り専用メモリ（ＣＤ−ＲＯＭ）、デジタル多用途ディスク（ＤＶＤ）、メモリスティック、フロッピーディスク、機械的符号化装置、例えば命令が記憶されているせん孔カードまたはスロット内突起構造、および上記の任意の適当な組み合わせを含む。ここで使用されるコンピュータが読み取り可能な記憶媒体は、例えば無線電波または他の自由に伝播される電磁波、導波路または他の伝送媒体を経由して伝播される電磁波（例えば、光ファイバーケーブルを通過するパルス光）、または電線を経由して伝送される電気信号のような、瞬時信号自体であると解釈されるものではない。
【０１０５】
ここで記述したコンピュータ可読プログラム命令は、コンピュータが読み取り可能な記憶媒体から各計算／処理機器にダウンロードされてもよいし、またはネットワーク、例えばインターネット、ローカルエリアネットワーク、広域ネットワークおよび／または無線ネットワークを介して外部のコンピュータまたは外部記憶装置にダウンロードされてもよい。ネットワークは、銅伝送ケーブル、光ファイバー伝送、無線伝送、ルーター、ファイアウォール、交換機、ゲートウェイコンピュータおよび／またはエッジサーバを含むことができる。各計算／処理機器内のネットワークアダプタカードまたはネットワークインタフェースはネットワークからコンピュータ可読プログラム命令を受信し、該コンピュータ可読プログラム命令を転送し、各計算／処理機器内のコンピュータが読み取り可能な記憶媒体に記憶する。
【０１０６】
本開示の動作を実行するためのコンピュータプログラム命令はアセンブラ命令、命令セットアーキテクチャ（ＩＳＡ）命令、機械語命令、機械依存命令、マイクロコード、ファームウェア命令、状態設定データ、または一つまたは複数種類のプログラミング言語の任意の組み合わせで書かれたソースコードまたは目標コードであってもよい。前記プログラミング言語は、Ｓｍａｌｌｔａｌｋ、Ｃ＋＋などのオブジェクト指向プログラミング言語、および「Ｃ」言語または類似するプログラミング言語などの一般的な手続き型プログラミング言語を含む。コンピュータ可読プログラム命令は、完全にユーザのコンピュータにおいて実行されてもよく、部分的にユーザのコンピュータにおいて実行されてもよく、スタンドアロンソフトウェアパッケージとして実行されてもよく、部分的にユーザのコンピュータにおいてかつ部分的にリモートコンピュータにおいて実行されてもよく、または完全にリモートコンピュータもしくはサーバにおいて実行されてもよい。リモートコンピュータに関与する場合、リモートコンピュータは、ローカルエリアネットワーク（ＬＡＮ）または広域ネットワーク（ＷＡＮ）を含む任意の種類のネットワークを経由してユーザのコンピュータに接続されてもよく、または、（例えばインターネットサービスプロバイダを利用してインターネットを経由して）外部コンピュータに接続されてもよい。いくつかの実施例では、コンピュータ可読プログラム命令の状態情報を利用することで、例えばプログラマブル論理回路、フィールドプログラマブルゲートアレイ（ＦＰＧＡ）またはプログラマブル論理アレイ（ＰＬＡ）などの、コンピュータ可読プログラム命令を実行して本開示の各態様を実現できる電子回路をパーソナライズする。
【０１０７】
なお、ここで本開示の実施例に係る方法、装置（システム）およびコンピュータプログラム製品のフローチャートおよび／またはブロック図を参照しながら本開示の各態様を説明したが、フローチャートおよび／またはブロック図の各ブロック、およびフローチャートおよび／またはブロック図の各ブロックの組み合わせは、いずれもコンピュータ可読プログラム命令によって実現できることを理解すべきである。
【０１０８】
これらのコンピュータ可読プログラム命令は、汎用コンピュータ、専用コンピュータまたは他のプログラマブルデータ処理装置のプロセッサへ提供されてもよい。それにより、これらの命令はコンピュータまたは他のプログラマブルデータ処理装置のプロセッサによって実行されると、フローチャートおよび／またはブロック図の一つまたは複数のブロックにおいて指定された機能／動作を実現する手段を創出するような機械が生じる。また、これらのコンピュータ可読プログラム命令は、コンピュータが読み取り可能な記憶媒体に記憶し、それによってコンピュータ、プログラマブルデータ処理装置および／または他の機器を特定の方式で動作させるようにしてもよい。それにより、命令が記憶されるコンピュータが読み取り可能な記憶媒体は、フローチャートおよび／またはブロック図の一つまたは複数のブロックにおいて指定された機能／動作の各態様を実現する命令を含む製品を含む。
【０１０９】
コンピュータ可読プログラムはコンピュータ、他のプログラマブルデータ処理装置、または他の機器にロードしてもよい。これにより、コンピュータ、他のプログラマブルデータ処理装置または他の機器において一連の動作ステップを実行させて、コンピュータ実施プロセスを生成する。それにより、コンピュータ、他のプログラマブルデータ処理装置、または他の機器において実行される命令はフローチャートおよび／またはブロック図の一つまたは複数のブロックにおいて指定された機能／動作を実現する。
【０１１０】
図面のうちフローチャートおよびブロック図は、本開示の複数の実施例に係るシステム、方法およびコンピュータプログラム製品の実現可能なシステムアーキテクチャ、機能および動作を示す。この点では、フローチャートまたはブロック図における各ブロックは一つのモジュール、プログラムセグメントまたは命令の一部分を代表することができ、前記モジュール、プログラムセグメントまたは命令の一部分は、指定された論理機能を実現するための一つまたは複数の実行可能な命令を含む。いくつかの置換としての実現形態では、ブロックに表記される機能は、図面に付したものと異なる順序で実現してもよい。例えば、二つの連続的なブロックは実質的に同時に実行してもよく、また、係る機能によって、それらは逆な順序で実行してもよい場合がある。なお、ブロック図および／またはフローチャートにおける各ブロック、およびブロック図および／またはフローチャートにおけるブロックの組み合わせは、指定される機能または動作を実行するハードウェアに基づく専用システムによって実現してもよいし、または専用ハードウェアとコンピュータ命令との組み合わせによって実現してもよいことに注意すべきである。
【０１１１】
以上、本開示の各実施例を記述したが、上記説明は例示的なものに過ぎず、網羅的なものではなく、かつ披露された各実施例に限定されるものでもない。当業者にとって、説明された各実施例の範囲および精神から逸脱することなく、様々な修正および変更が明らかである。本明細書に選ばれた用語は、各実施例の原理、実際の適用またはマーケティングにおける技術への技術的改善を好適に解釈するか、または当該技術的分野の一般的な当業者に本文に披露された各実施例を理解させるためのものである。
【図１】
2022500773.tif 000003
【図２】
2022500773.tif 000004
【図３】
2022500773.tif 000005
【図４】
2022500773.tif 000006
【図５】
2022500773.tif 000007
【図６】
2022500773.tif 000008
【図７】
2022500773.tif 000009
【図８】
2022500773.tif 000010
【手続補正書】
【提出日】令和3年3月16日(2021.3.16)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得することと、
前記少なくとも１つの遺伝子シーケンシング断片の属性情報に基づいて、前記変異遺伝子座候補の、遺伝子座の位置に関連する配列特徴および非配列特徴を特定することと、
前記配列特徴および前記非配列特徴に基づいて、前記変異遺伝子座候補の遺伝子変異を認識することと、を含むことを特徴とする、遺伝子変異認識方法。
【請求項２】
前記属性情報は、配列属性情報を含み、
前記少なくとも１つの遺伝子シーケンシング断片の属性情報に基づいて、前記変異遺伝子座候補の配列特徴を特定することは、
前記変異遺伝子座候補の遺伝子位置情報に基づいて、前記変異遺伝子座候補が位置する所定の遺伝子座区間を特定することと、
前記少なくとも１つの遺伝子シーケンシング断片の、前記所定の遺伝子座区間内に位置する各遺伝子座の、遺伝子座の位置に関連しかつ遺伝子属性を表す配列属性情報を取得することと、
前記所定の遺伝子座区間に位置する各遺伝子座の配列属性情報に基づいて、前記変異遺伝子座候補の配列特徴を生成することと、を含むことを特徴とする、請求項１に記載の方法。
【請求項３】
前記少なくとも１つの遺伝子シーケンシング断片の、前記所定の遺伝子座区間内に位置する各遺伝子座の配列属性情報を取得することは、
前記各遺伝子座において前記少なくとも１つの遺伝子シーケンシング断片の遺伝子型を特定することと、
前記各遺伝子座において遺伝子を各遺伝子型毎にカウントすることと、を含むことを特徴とする、請求項２に記載の方法。
【請求項４】
前記少なくとも１つの遺伝子シーケンシング断片の、前記所定の遺伝子座区間内に位置する各遺伝子座の配列属性情報を取得することは、
各遺伝子シーケンシング断片の遺伝子配列と参照ゲノム配列との比較結果に基づいて、前記各遺伝子座において各遺伝子シーケンシング断片の欠失遺伝子の遺伝子型を特定することと、
前記各遺伝子座において前記少なくとも１つの遺伝子シーケンシング断片の欠失遺伝子を各遺伝子型毎にカウントすることと、を含むことを特徴とする、請求項２に記載の方法。
【請求項５】
前記少なくとも１つの遺伝子シーケンシング断片の、前記所定の遺伝子座区間内に位置する各遺伝子座の配列属性情報を取得することは、
各遺伝子シーケンシング断片の遺伝子配列と参照ゲノム配列との比較結果に基づいて、前記各遺伝子座において各遺伝子シーケンシング断片の挿入遺伝子の遺伝子型を特定することと、
前記各遺伝子座において前記少なくとも１つの遺伝子シーケンシング断片の挿入遺伝子を各遺伝子型毎にカウントすることと、を含むことを特徴とする、請求項２に記載の方法。
【請求項６】
前記配列属性情報は、
参照遺伝子の遺伝子型、各遺伝子型毎の遺伝子数、各遺伝子型毎の欠失遺伝子数、各遺伝子型毎の挿入遺伝子数の少なくとも１つを含むことを特徴とする、請求項１から５のいずれか一項に記載の方法。
【請求項７】
前記属性情報は、非配列属性情報を含み、
前記少なくとも１つの遺伝子シーケンシング断片の属性情報に基づいて、前記変異遺伝子座候補の非配列特徴を特定することは、
前記少なくとも１つの遺伝子シーケンシング断片の、遺伝子座の位置に関連しなくかつ遺伝子属性を表す非配列属性情報を取得することと、
前記少なくとも１つの遺伝子シーケンシング断片の非配列属性情報に基づいて、前記変異遺伝子座候補の非配列特徴を特定することと、を含むことを特徴とする、請求項１から６のいずれか一項に記載の方法。
【請求項８】
前記非配列情報は、
比較品質、プラス鎖およびマイナス鎖の偏り、遺伝子シーケンシング断片の長さ、エッジの偏りの少なくとも１つを含むことを特徴とする、請求項７に記載の方法。
【請求項９】
前記少なくとも１つの遺伝子シーケンシング断片の非配列属性情報に基づいて、前記変異遺伝子座候補の非配列特徴を特定することは、
各遺伝子シーケンシング断片の各遺伝子座の比較品質に基づいて、各遺伝子シーケンシング断片の比較品質を特定することであって、前記比較品質は、遺伝子シーケンシング断片の各遺伝子配列毎のシーケンシングの正確性を表すことと、
各遺伝子シーケンシング断片の比較品質に基づいて、前記変異遺伝子座候補に対応する非配列特徴を特定することと、を含むことを特徴とする、請求項８に記載の方法。
【請求項１０】
前記少なくとも１つの遺伝子シーケンシング断片の非配列属性情報に基づいて、前記変異遺伝子座候補の非配列特徴を特定することは、
各遺伝子シーケンシング断片の属する遺伝子鎖がプラス鎖またはマイナス鎖のいずれであるかの情報に基づいて、前記少なくとも１つの遺伝子シーケンシング断片の、プラス鎖である遺伝子鎖とマイナス鎖である遺伝子鎖との割合を特定することと、
前記プラス鎖である遺伝子鎖とマイナス鎖である遺伝子鎖との割合に基づいて、前記変異遺伝子座候補に対応する非配列特徴を特定することと、を含むことを特徴とする、請求項８に記載の方法。
【請求項１１】
前記配列特徴および前記非配列特徴に基づいて、前記変異遺伝子座候補の遺伝子変異を認識することは、
前記配列特徴および前記非配列特徴を統合し、前記変異遺伝子座候補の統合特徴を得ることと、
前記変異遺伝子座候補の統合特徴に基づいて、前記変異遺伝子座候補の遺伝子の変異性を示す変異値を得ることと、
前記変異値が所定の閾値以上である場合、前記変異遺伝子座候補の遺伝子に変異があると判定することと、を含むことを特徴とする、請求項１１に記載の方法。
【請求項１２】
前記変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得することは、
体細胞遺伝子による遺伝子シーケンシングで得られた遺伝子シーケンシング断片を取得することと、
前記遺伝子シーケンシング断片の遺伝子配列を参照ゲノム配列と比較し、比較結果を得ることと、
前記比較結果に基づいて、前記体細胞遺伝子に異常がある変異遺伝子座候補を特定することと、
前記変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得することと、を含むことを特徴とする、請求項１から１１のいずれか一項に記載の方法。
【請求項１３】
変異遺伝子座候補に対応する少なくとも１つの遺伝子シーケンシング断片を取得するための取得モジュールと、
前記少なくとも１つの遺伝子シーケンシング断片の属性情報に基づいて、前記変異遺伝子座候補の、遺伝子座の位置に関連する配列特徴および非配列特徴を特定するための特定モジュールと、
前記配列特徴および前記非配列特徴に基づいて、前記変異遺伝子座候補の遺伝子変異を認識するための認識モジュールと、を含むことを特徴とする、遺伝子変異認識装置。
【請求項１４】
プロセッサと、
プロセッサが実行可能な命令を記憶するメモリと、を含み、
前記プロセッサは、前記実行可能な命令を呼び出すことによって請求項１から１２のいずれか一項に記載の方法を実現することを特徴とする遺伝子変異認識装置。
【請求項１５】
コンピュータプログラム命令が記憶されている不揮発性コンピュータ可読記憶媒体であって、前記コンピュータプログラム命令がプロセッサにより実行されると、請求項１から１２のいずれか一項に記載の方法を実現することを特徴とする、不揮発性コンピュータ可読記憶媒体。
【国際調査報告】
2022500773.tif 000011
2022500773.tif 000012
2022500773.tif 000013
2022500773.tif 000014