The orchestration of coding and non-coding genes at the CDH13 locus in coronary artery disease

Abstract

Abstract Background The CDH13 locus was associated with coronary artery disease (CAD) by genome-wide association studies (GWAS). However, the causal gene(s) and underlying disease mechanisms at this locus remain undetermined. Purpose We intended to pinpoint and functionally validate the causal coding and non-coding genes at the CDH13 locus and delineate the underlying mechanisms for CAD. Methods and results We conducted the transcriptome-wide association analysis (TWAS) using nine types of CAD-relevant tissues from Stockholm-Tartu Atherosclerosis Reverse Network Engineering Task (STARNET) and Genotyped Tissue Expression (GTEx) projects and identified CDH13 (T-cadherin) and four long non-coding RNA (lncRNA) genes as candidate causal genes at The CDH13 locus. One of the antisense(AS) lncRNAs, CDH13-AS2, was predicted to bind on the 3’ untranslated region (UTR) of CDH13 using LncRRIsearch and RNAup servers. The binding interaction was validated by dCas13-medicated RNA immunoprecipitation. Weighted correlation network analysis (WGCNA) using transcriptomic data from CAD-relevant tissues of ~600 human subjects suggested a positive expression correlation between CDH13 and CDH13-AS2 in artery tissues. CRISPR/Cas9-based knockout (KO) of CDH13 or CDH13-AS2 in human umbilical vein endothelial cells (HUVECs) indicated synergistically atherogenic phenotypes, including reduced proliferation and migration, and increased apoptosis and monocyte adhesion. The data matched with the downregulation of CDH13 in artery samples of atherosclerosis patients. The cellular phenotypes were also in line with our in vivo data of Cdh13 and Apoe double-KO (Cdh13-/-/Apoe-/-) mice, which showed increased aortic lesion areas when fed on a high-fat diet compared to Apoe-/- mice. Furthermore, prediction using TargetScan, miRWalk, and scanMiRApp databases identified the potential binding of several EC-expression microRNAs (mir) on 3’UTR of CDH13. Several disease-associated variants were mapped to the predicted binding sites of a few microRNAs, including mir-let7, mir-30, and mir-125. The dual luciferase-based RNA interference (RNAi) assay confirmed the binding of mir-let7 on 3’UTR of CDH13. Transcriptional activation of CDH13-AS2 by the enzymatically inactivated Cas9 (dCas9) system diminished the binding of mir-let7 and increased CDH13 expression. Conclusions At the CDH13 locus, CDH13 and CDH13-AS2 are two causal genes, whose downregulation likely contributes to CAD. CDH13-AS2 bound and stabilized CDH13 partially by preventing mir-let7 mediated degradation of CDH13 mRNA, suggesting therapeutic potentials.