Abstract
Deciphering expression and variants in cardiovascular disease genes among heart failure population for precision medicine Cardiovascular disease (CVD) is the leading cause of death in the United States and around the globe. 1 Despite significant advancements in CVD diagnostics, prevention, and treatment , approximately half of the affected patients reportedly die within 5 years of receiving a diagnosis. 2 The risk factors contributing to the development of CVD and response to therapy in an individual patient are highly variable. Evidence from the Framingham Heart Study suggests that CVD has a complex multifactorial aetiology including a genetic component. 3 Genomics information, including high-quality sequenced DNA and RNA sequencing (RNA-seq) of transcribed genes, informs us of a CVD patient's inherent genetic makeup with the most comprehensive view of the genome. 4 DNA-based gene variant detection when combined with RNA-seq-driven gene expression analysis has the potential to reveal novel and sensitive biomarkers and stratify CVD patient populations based on their disease risk. 5,6 The genetic variants predisposing to CVD span from rare and deleterious mutations that may be responsible for familial aggregation. 6 Investigating differentially expressed genes (DEGs) and disease-causing variants can support finding the root cause of uncertainties in patient care. 7 Understanding of the genetic basis of complex CVD can hamper genetic risk scoring which can now outperform traditional risk factors in risk prediction. 7 Various genomics studies have been conducted recently to discover underlying genetic aetiology in CVD patients, especially those suspected of HF disease. 7 Several CVD-related genes have been reported with significant mutations and expression differences among CVD patients. 7 Application of intelligent and integrative data analysis approaches involving genomics and transcriptomics data will not only help understand the pathophysiology of CVD but also reduce heteroge-neity in disease subtypes. To improve the deciphering of CVD mechanisms, here, we report our investigation of expression and variants among known genes that are responsible for the development of CVD. 5,6 Supporting this study, we developed an open-source, cross-platform, interactive, and user-friendly bioinformatics pipeline (GVViZ) for RNA-seq data pre-processing, and gene expression data analysis, annotation with relevant diseases, and heatmap visualization without requiring a strong computational background from the user. 8 In addition, we have developed a new bioinformatics pipeline (JWES), for variant discovery and interpretation, and big data modelling and visualization. 9 In this study, we conducted analysis of gene expression, disease-causing gene-variants, and associated phenotypes among CVD populations, with a focus on high-risk Heart Failure (HF). 5,6 We built a cohort of CVD patients, 40 male and 21 female individuals (n = 61), aged between 45 to 92, with self-described race (42 Whites, 7 Blacks or African Americans, 1 Asian, and 11 of unknown race). We collected blood samples , performed RNA-seq and gene expression analysis to generate transcriptomic profiles (Figure 1A). We performed in-depth gene expression analysis and annotation of RNA-seq data using GVViZ revealed regulation of genes known for HF (Figure 1B) and other CVDs (Figure 1C). Subsequent analyses were performed based on race and gender. 5 Our analysis identified altered expression pathways of genes with gender differences in middle-aged to frail CVD patients. 5 Next, we processed whole-genome sequencing (WGS) data using JWES and identified mutations among annotated genes for HF and other CVD patients (Figure 1D). We annotated these mutations to identify functional and nonfunctional mutations potentially associated with HF and other CVDs (Figure 1D). Mutation percentage in CVD genes was notably higher in HF patients compared to other CVD phenotypes (Figure 2A, B). In total, we detected 1 039 750 single nucleotide variants (SNV), insertion and deletion events. The most common mutation types in HF and CVD genes were intronic, flanking (5′ and 3′ flank) mutations. 6 Mutations in these genes have been linked to aberrant expression in CVD. WGS allowed us to do in-depth analysis of CVD genes as RNA-seq cannot detect any of the variants located in noncoding DNA regions. 6 We identified mutations among four genes with altered expression and having significant mutations associated with HF and other CVDs (FLNA, CST3, LGALS3, and HBA1). Notables are the missense mutations in FLNA, CT3, and LGALS3. Next,