WGS

材料と対象者

2003年から2017年度にバイオバンク・ジャパンに登録された症例: 1,964名 心筋梗塞(ICD10: I21): 1,765症例 認知症(ICD10: F00-03): 199症例

実験方法

WGS

プラットフォーム

Illumina [HiSeq X Five]

試料説明

末梢血から抽出したgDNA

リード重複除去

Picard 2.10.6

再アライメント・塩基品質補正

GATK 3.7

ライブラリ構築

TruSeq DNA PCR-Free Prep kit

断片化

超音波断片化

リードタイプ

Paired-end

リード長

151 bp

解析方法

GATK 3.7 HaplotypeCaller

フィルタリング

リードのbase qualityが全体的に悪い検体、リード毎の%GCの結果にて異常を示した検体を除去。 Alignment後、Low mapping rate検体、Insert sizeがおかしい検体、メタデータの性別情報とalingment結果より推定される性別情報が不一致な検体、性染色体異常疑いの検体を除去。 Genotyping時に、VQSR、DP/GP filter (DP < 5, GQ < 20, DP > 60 && GQ < 95を除去)、heterozygosity filter (F>=0.05 を除去)、HWE filter (p < 10-6を除去)、Repeat & Low Complexity filterを実施。 1000 genomes projectと合わせたPCAを実施し、日本人クラスタから大きく外れる検体を除外。 その後、Genome-In-A-Bottleプロジェクトから公開されているHighConfidenceRegionリストに記載のある領域のバリアントにフラグを付与。

マッピング

BWA mem 0.7.12

マッピング品質

GATK 3.7 HaplotypeCallerで変異コール時にMAPQ<20のリードを除外

リファレンス配列

GRCh37/hg19(hs37d5)

カバレッジ

心筋梗塞: 15.0x、認知症: 30.0x

SNV数

76,768,387(常染色体) 2,898,518(X 染色体)

INDEL数

10,202,908(常染色体) 410,435(X 染色体)

JGAデータセットAccession

JGAD000496(bam、vcf)

総データ量

188.4 TB(fastq、bam、vcf)

利用ポリシー

NBDC policy