A Framework for Detecting Noncoding Rare-variant Associations of Large-scale Whole-genome Sequencing Studies

Large-scale whole-genome sequencing studies have enabled analysis of noncoding rare-variant (RV) associations with complex human diseases and traits. Variant-set analysis is a powerful approach to study RV association. However, existing methods have limited ability in analyzing the noncoding genome. We propose a computationally efficient and robust noncoding RV association detection framework, STAARpipeline, to automatically annotate a whole-genome sequencing study and perform flexible noncoding RV association analysis, including gene-centric analysis and fixed window-based and dynamic window-based non-gene-centric analysis by incorporating variant functional annotations. In gene-centric analysis, STAARpipeline uses STAAR to group noncoding variants based on functional categories of genes and incorporate multiple functional annotations. In non-gene-centric analysis, STAARpipeline uses SCANG-STAAR to incorporate dynamic window sizes and multiple functional annotations. We apply STAARpipeline to identify noncoding RV sets associated with four lipid traits in 21,015 discovery samples from the Trans-Omics for Precision Medicine (TOPMed) program and replicate several of them in an additional 9,123 TOPMed samples. We also analyze five non-lipid TOPMed traits.

Citing Articles

Streamlining Large-Scale Genomic Data Management: Insights from the UK Biobank Whole-Genome Sequencing Data.

Li X, Wood A, Yuan Y, Zhang M, Huang Y, Hawkes G medRxiv. 2025; .

PMID: 39974066 PMC: 11838927. DOI: 10.1101/2025.01.27.25321225.

Assessment of the functionality and usability of open-source rare variant analysis pipelines.

Riccio C, Jansen M, Thalen F, Koliopanos G, Link V, Ziegler A Brief Bioinform. 2025; 26(1).

PMID: 39907318 PMC: 11795309. DOI: 10.1093/bib/bbaf044.

RetroFun-RVS: A Retrospective Family-Based Framework for Rare Variant Analysis Incorporating Functional Annotations.

Mangnier L, Ruczinski I, Ricard J, Moreau C, Girard S, Maziade M Genet Epidemiol. 2025; 49(2):e70001.

PMID: 39876583 PMC: 11775437. DOI: 10.1002/gepi.70001.

Zim4rv: an R package to modeling zero-inflated count phenotype on regional-based rare variants.

Liu X, Li Y, Fan Q BMC Bioinformatics. 2025; 26(1):18.

PMID: 39819419 PMC: 11740424. DOI: 10.1186/s12859-024-06029-5.

Review on GPU accelerated methods for genome-wide SNP-SNP interactions.

Ren W, Liang Z Mol Genet Genomics. 2024; 300(1):10.

PMID: 39738695 DOI: 10.1007/s00438-024-02214-6.

References

Manolio T, Collins F, Cox N, Goldstein D, Hindorff L, Hunter D . Finding the missing heritability of complex diseases. Nature. 2009; 461(7265):747-53. PMC: 2831613. DOI: 10.1038/nature08494. View

Wainschtein P, Jain D, Zheng Z, Cupples L, Shadyab A, McKnight B . Assessing the contribution of rare variants to complex trait heritability from whole-genome sequence data. Nat Genet. 2022; 54(3):263-273. PMC: 9119698. DOI: 10.1038/s41588-021-00997-7. View

Hernandez R, Uricchio L, Hartman K, Ye C, Dahl A, Zaitlen N . Ultrarare variants drive substantial cis heritability of human gene expression. Nat Genet. 2019; 51(9):1349-1355. PMC: 6730564. DOI: 10.1038/s41588-019-0487-7. View

Taliun D, Harris D, Kessler M, Carlson J, Szpiech Z, Torres R . Sequencing of 53,831 diverse genomes from the NHLBI TOPMed Program. Nature. 2021; 590(7845):290-299. PMC: 7875770. DOI: 10.1038/s41586-021-03205-y. View

Flannick J, Mercader J, Fuchsberger C, Udler M, Mahajan A, Wessel J . Exome sequencing of 20,791 cases of type 2 diabetes and 24,440 controls. Nature. 2019; 570(7759):71-76. PMC: 6699738. DOI: 10.1038/s41586-019-1231-2. View

Van Hout C, Tachmazidou I, Backman J, Hoffman J, Liu D, Pandey A . Exome sequencing and characterization of 49,960 individuals in the UK Biobank. Nature. 2020; 586(7831):749-756. PMC: 7759458. DOI: 10.1038/s41586-020-2853-0. View

Zhang F, Lupski J . Non-coding genetic variants in human disease. Hum Mol Genet. 2015; 24(R1):R102-10. PMC: 4572001. DOI: 10.1093/hmg/ddv259. View

Khurana E, Fu Y, Chakravarty D, Demichelis F, Rubin M, Gerstein M . Role of non-coding sequence variants in cancer. Nat Rev Genet. 2016; 17(2):93-108. DOI: 10.1038/nrg.2015.17. View

Lee P, Lee C, Li X, Wee B, Dwivedi T, Daly M . Principles and methods of in-silico prioritization of non-coding regulatory variants. Hum Genet. 2017; 137(1):15-30. PMC: 5892192. DOI: 10.1007/s00439-017-1861-0. View

10.

Moore J, Purcaro M, Pratt H, Epstein C, Shoresh N, Adrian J . Expanded encyclopaedias of DNA elements in the human and mouse genomes. Nature. 2020; 583(7818):699-710. PMC: 7410828. DOI: 10.1038/s41586-020-2493-4. View

11.

Bansal V, Libiger O, Torkamani A, Schork N . Statistical analysis strategies for association studies involving rare variants. Nat Rev Genet. 2010; 11(11):773-85. PMC: 3743540. DOI: 10.1038/nrg2867. View

12.

Lee S, Abecasis G, Boehnke M, Lin X . Rare-variant association analysis: study designs and statistical tests. Am J Hum Genet. 2014; 95(1):5-23. PMC: 4085641. DOI: 10.1016/j.ajhg.2014.06.009. View

13.

Kiezun A, Garimella K, Do R, Stitziel N, Neale B, McLaren P . Exome sequencing and the genetic basis of complex traits. Nat Genet. 2012; 44(6):623-30. PMC: 3727622. DOI: 10.1038/ng.2303. View

14.

Li B, Leal S . Methods for detecting associations with rare variants for common diseases: application to analysis of sequence data. Am J Hum Genet. 2008; 83(3):311-21. PMC: 2842185. DOI: 10.1016/j.ajhg.2008.06.024. View

15.

Wu M, Lee S, Cai T, Li Y, Boehnke M, Lin X . Rare-variant association testing for sequencing data with the sequence kernel association test. Am J Hum Genet. 2011; 89(1):82-93. PMC: 3135811. DOI: 10.1016/j.ajhg.2011.05.029. View

16.

Li X, Li Z, Zhou H, Gaynor S, Liu Y, Chen H . Dynamic incorporation of multiple in silico functional annotations empowers rare variant association analysis of large whole-genome sequencing studies at scale. Nat Genet. 2020; 52(9):969-983. PMC: 7483769. DOI: 10.1038/s41588-020-0676-4. View

17.

Morrison A, Huang Z, Yu B, Metcalf G, Liu X, Ballantyne C . Practical Approaches for Whole-Genome Sequence Analysis of Heart- and Blood-Related Traits. Am J Hum Genet. 2017; 100(2):205-215. PMC: 5294677. DOI: 10.1016/j.ajhg.2016.12.009. View

18.

Li Z, Li X, Liu Y, Shen J, Chen H, Zhou H . Dynamic Scan Procedure for Detecting Rare-Variant Association Regions in Whole-Genome Sequencing Studies. Am J Hum Genet. 2019; 104(5):802-814. PMC: 6507043. DOI: 10.1016/j.ajhg.2019.03.002. View

19.

He Z, Xu B, Buxbaum J, Ionita-Laza I . A genome-wide scan statistic framework for whole-genome sequence data analysis. Nat Commun. 2019; 10(1):3018. PMC: 6616627. DOI: 10.1038/s41467-019-11023-0. View

20.

Natarajan P, Peloso G, Zekavat S, Montasser M, Ganna A, Chaffin M . Deep-coverage whole genome sequences and blood lipids among 16,324 individuals. Nat Commun. 2018; 9(1):3391. PMC: 6107638. DOI: 10.1038/s41467-018-05747-8. View