Streamlining Large-Scale Genomic Data Management: Insights from the UK Biobank Whole-Genome Sequencing Data
Preprint
Abstract
Biobank-scale Whole-Genome Sequencing (WGS) studies are increasingly pivotal in unraveling the genetic bases of diverse health outcomes. However, managing and analyzing these datasets’ sheer volume and complexity presents significant challenges. We propose vcf2agds , an all-in-one toolkit that efficiently converts WGS data from Variant Call Format (VCF) format to the annotated Genomic Data Structure (aGDS) format, significantly reducing data size while supporting seamless genomic and functional data integration for comprehensive genetic analyses. The toolkit was applied to the UK Biobank 500k WGS data, resulting in twenty-three aGDS files, one for each chromosome, which collectively compressed 1,473.85 Tebibytes of pVCF data into 1.10 Tebibytes. Utilizing these aGDS files, we conducted a functionally informed rare variant association analysis of total cholesterol employing the STAARpipeline and detected 480 genome-wide significant coding and noncoding associations. Overall, vcf2agds offers a streamlined approach facilitating the efficient management and analysis of biobank-scale WGS data across hundreds of thousands of samples.Keywords
This publication has 15 references indexed in Scilit:
- Whole-genome sequencing in 333,100 individuals reveals rare non-coding single variant and aggregate associations with heightNature Communications, 2024
- Genomic data in the All of Us Research ProgramNature, 2024
- A framework for detecting noncoding rare-variant associations of large-scale whole-genome sequencing studiesNature Methods, 2022
- The sequences of 150,119 genomes in the UK BiobankNature, 2022
- Sequencing of 53,831 diverse genomes from the NHLBI TOPMed ProgramNature, 2021
- Dynamic incorporation of multiple in silico functional annotations empowers rare variant association analysis of large whole-genome sequencing studies at scaleNature Genetics, 2020
- Graphtyper enables population-scale genotyping using pangenome graphsNature Genetics, 2017
- A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing dataBioinformatics, 2011
- The variant call format and VCFtoolsBioinformatics, 2011
- PLINK: A Tool Set for Whole-Genome Association and Population-Based Linkage AnalysesThe American Journal of Human Genetics, 2007