Decoding the Variant-to-Function Relationship for<i>LIPA</i>, a Risk Locus for CAD
Fang Li,Elise Flynn,Jianting Shi,Xun Wu,Ziyi Wang,Chenyi Xue,Haoxiang Cheng,Yujiao Meng,Jian Cui,Yizhou Zhu,Annie Rozenblyum,Jeana Chun,Antonio Hernandez-Ono,Babak Razani,Marit Westerterp,Robert C Bauer,Yousin Suh,Ke Hao,Tuuli Lappalainen,Hanrui Zhang
DOI: https://doi.org/10.1101/2022.11.12.516293
2022-01-01
Abstract:ABSTRACTBackgroundGenome-wide association studies revealed a robust association between genetic variants in theLIPA(lysosomal acid lipase) gene and coronary artery diseases (CAD), but not lipid traits. QTL studies support that the risk alleles ofLIPACAD variants are associated with higherLIPAmRNA and enzyme activity in human monocytes. Yet the variant-to-function relationship and how this important locus impacts disease etiology has not been fully established. Herein, we aim to determine the causal variant(s), involved cell type, and the target gene, establish the causality of the variant-to-function relationship, and elucidate how increased myeloid LIPA impacts atherosclerosisin vivo.MethodsWe apply functional genomic datasets, post-GWAS prioritization pipelines, and molecular biology techniques, incuding eQTL, enzyme activity-QTL, high-resolution Tri-HiC, ChIP-seq, and site-directed mutagenesis and luciferase assay to connect functional variants to the candidate genes in the causal cell type. To establish how increased myeloidLIPAimpacts atherosclerosis, we generated myeloid-specificLipaoverexpression mice(LipaTg).ResultsPost-GWAS pipelines supportLIPAas the candidate causal gene at the locus. In human monocyte-derived macrophages,LIPAmRNA, protein and enzyme activity were higher in the risk allele carriers of CAD variants. High-resolution Tri-HiC and luciferase assay confirmed an intronic enhancer region showing strong interaction with theLIPApromoter. Within the enhancer region, the risk alleles of rs1412444/rs1412445 and rs1320496 demonstrate enhanced binding to PU. 1, and acted as the functional variants with risk alleles leading to increased enhancer activity. The risk allele of rs1320496 is predicted to create a motif binding site for PU.1. The functional genomic data together support thatLIPAis the candidate causal gene in the locus, and the risk alleles of CAD led to increased LIPA in a myeloid cell-specific manner. Consistently, mice with myeloid-specificLipaoverexpression on aLdlr-/-background showed significantly increased atherosclerotic lesion size and lesion macrophage area without affecting plasma cholesterol. ScRNA-seq analysis showed thatLipaTgled to reduced lipid-enriched yet increased inflammatory macrophage subsets, and activation chemokine signaling pathway. This was further confirmed by reduced neutral lipid accumulation in both plaque and peritoneal macrophages and significantly increased monocytes infiltration into the lesion inLipaTgmice.ConclusionsWe established thatLIPArisk alleles drive increased myeloid LIPA and aggravate atherosclerosis.CLINICAL PERSPECTIVEWhat is New?CAD GWAS variants at theLIPAlocus led to increased macrophage LIPA expression and enzyme activity.Myeloid-specific overexpression ofLipaexacerbates atherosclerosis.Our study connected the genetic variation to the involved cell type and the target gene, and the disease mechanism for this important locus.What are the Clinical Implications?GWAS and meta-analyses have identified over 200 loci for CAD. Establishing the candidate genes and their mechanistic studies inform novel biological mechanisms and therapeutic application.There is strong statistical evidence linkingLIPAwith CAD. By leveraging functional genomic studies and transgenic mice, our work established the direct causality thatLIPArisk alleles drive increased myeloid LIPA and aggravate atherosclerosis. Establishing the variant-to-function relationship for this locus informs that increasing myeloid LIPA may not be a therapeutic strategy for CAD, despite the essential role of LIPA in regulating lysosomal lipid metabolism.