Skip to content

Aggregated Variant Calls (AggV2)

AggV2 is an aggregate of a large number of gVCFs. It was made from release 10 of the 100kGP project and comprises 78,195 germline genomes. These were all the consented germline genomes that were aligned to GRCh38 and passed quality control available in release 10.

As part of AggV2, we provide:

AggV2 contains information on participants who have since withdrawn consent from research. You cannot use them in any new analyses. It is extremely important to remove these samples from your analyses and only use samples included in the latest data release.

The list of samples for the consented participants can be found in the aggregate_gvcf_sample_stats table in LabKey, for the latest data release, or in the current samples file, located in /gel_data_resources/main_programme/aggregation/aggregate_gVCF_strelka/aggV2/docs/.

To filter the aggregate to these samples, all bcftools commands should include the flag -S /gel_data_resources/main_programme/aggregation/aggregate_gVCF_strelka/aggV2/docs/<nameofcurrentfile>.

Submit a ticket to the Genomics England Service desk if you are unsure of how to filter the dataset for any other use.

Description

We have aggregated 78,195 germline gVCFs (genomic VCFs) from the 100,000 Genomes Project which we made available as a multi-sample VCF dataset (aggV2). aggV2 comprises over 722 million annotated single nucleotide variants and small indels (≤50 bp) from quality controlled rare disease and cancer germline whole genomes.

All samples in the dataset were sequenced with 150bp paired-end reads in a single lane of an Illumina HiSeq X instrument and uniformly processed on the Illumina North Star Version 4 Whole Genome Sequencing Workflow (NSV4, version 2.6.53.23); which comprises the iSAAC Aligner (version 03.16.02.19) and Starling Small Variant Caller(version 2.4.7). Samples were aligned to the Homo Sapiens NCBI GRCh38 assembly with decoys. 

The dataset was constructed from the aggregation of single-sample gVCFs using the Illumina software gVCF genotyper (version: 2019.02.26). Variant normalisation and decomposition was implemented by vt (version 0.57721). Genomic annotation and calculation of allele statistics (count, frequency etc.) was performed using Ensembl VEP and bcftools respectively. 

The multi-sample VCF is split into 1,371 roughly equal 'chunks' across the genome for faster processing. Each chunk contains the full set of samples and is in the VCF.gz file format with accompanying tabix index files (.tbi). Chromosomes 1-22, X, Y, and M are included. 

Extended details

Each step of the pipeline to generate aggV2 is documented in the sections below: 

FAQs and code book

An FAQs section regarding all aggV2 queries can be found here: aggV2 FAQ

A code book of popular queries to help you use aggV2 is found here: aggV2 Code Book

Manifest

The aggV2 dataset comprises four main parts: 

  1. A multi-sample VCF file for each chunk containing the genotypes and per variant quality metrics and filter flags
  2. A corresponding VCF file for each chunk containing the functional (genomic) annotation and allele statistics for all variants
  3. The aggregate_gvcf_sample_stats table in the current 100kGP LabKey folder which contains all sample quality metrics and accompanying meta-data
  4. Associated files, such as Principal Components across all included samples, information on sample relatedness, and assignment of predicted super-population to each sample. Some of this information is also provided in LabKey, in the aggregate_gvcf_sample_stats table mentioned above.

Location

All aggV2 outputs can be found in the following folder within the Genomics England Research Environment: 

/gel_data_resources/main_programme/aggregation/aggregate_gVCF_strelka/aggV2/

This folder is accessible from the Desktop Environment and from the HPC as shown below:

Desktop:

HPC:

Phased data for aggV2

As of 100kGP Data Release v15 we have made phased data available on behalf of Sinan Shi from the University of Oxford. Documentation and details can be found here: Aggv2 Phased Data (Provided by University of Oxford)

Masked pgen files for aggV2

Further to the multi-sample VCF files for aggV2, we also provide masked pgen files (per chunk and per chromosome) in the following directory:

/gel_data_resources/main_programme/aggregation/aggregate_gVCF_strelka/aggV2/genomic_data/pgen_allvariants

For these files, we masked low-quality genotypes by setting them to missing using the bcftools setGT module: For autosomes, we masked genotypes having DP<10 or genotype quality (GQ)<20 or heterozygote genotypes failing an ABratio binomial test with P-value < 10−3. For chrX, we masked females as for autosomes. We masked male genotypes having DP<5 or GQ<20.

Help and support

Please reach out via the Genomics England Service Desk for any issues related to the aggV2 aggregation or companion datasets, including "aggV2" in the title/description of your inquiry.