Integration of multidimensional splicing data and GWAS summary statistics for risk gene discovery

• Published in: PLoS genetics Vol. 18; no. 6; p. e1009814
• Main Authors: Ji, Ying, Wei, Qiang, Chen, Rui, Wang, Quan, Tao, Ran, Li, Bingshan
• Format: Journal Article
• Language: English
• Published: United States Public Library of Science 01.06.2022; Public Library of Science (PLoS)
• Subjects: Alzheimer's disease; Biology and Life Sciences; Chi-square test; Cholesterol; Correlation analysis; Gene expression; Genetic diversity; Genome-wide association studies; Genome-Wide Association Study - methods; Genomes; Genotypes; Humans; Medicine and Health Sciences; Mental disorders; Neurodegenerative diseases; Phenotype; Polymorphism, Single Nucleotide; Power; Quantitative Trait Loci - genetics; Regularization methods; Research and Analysis Methods; Schizophrenia; Simulation; Sparsity; Splicing; Statistical analysis; Statistics
• ISSN: 1553-7404; 1553-7390; 1553-7404
• DOI: 10.1371/journal.pgen.1009814

A common strategy for the functional interpretation of genome-wide association study (GWAS) findings has been the integrative analysis of GWAS and expression data. Using this strategy, many association methods (e.g., PrediXcan and FUSION) have been successful in identifying trait-associated genes via mediating effects on RNA expression. However, these approaches often ignore the effects of splicing, which can carry as much disease risk as expression. Compared to expression data, one challenge to detect associations using splicing data is the large multiple testing burden due to multidimensional splicing events within genes. Here, we introduce a multidimensional splicing gene (MSG) approach, which consists of two stages: 1) we use sparse canonical correlation analysis (sCCA) to construct latent canonical vectors (CVs) by identifying sparse linear combinations of genetic variants and splicing events that are maximally correlated with each other; and 2) we test for the association between the genetically regulated splicing CVs and the trait of interest using GWAS summary statistics. Simulations show that MSG has proper type I error control and substantial power gains over existing multidimensional expression analysis methods (i.e., S-MultiXcan, UTMOST, and sCCA+ACAT) under diverse scenarios. When applied to the Genotype-Tissue Expression Project data and GWAS summary statistics of 14 complex human traits, MSG identified on average 83%, 115%, and 223% more significant genes than sCCA+ACAT, S-MultiXcan, and UTMOST, respectively. We highlight MSG's applications to Alzheimer's disease, low-density lipoprotein cholesterol, and schizophrenia, and found that the majority of MSG-identified genes would have been missed from expression-based analyses. Our results demonstrate that aggregating splicing data through MSG can improve power in identifying gene-trait associations and help better understand the genetic risk of complex traits.