Identification of Four Potential Biomarkers Associated With Coronary Artery Disease in Non-diabetic Patients by Gene Co-expression Network Analysis

Potential role of soluble CD40 receptor in chronic inflammatory diseases

The CD40 receptor and its ligand CD154 are widely expressed in various immune-competent cells. Interaction of CD154 with CD40 is essential for B-cell growth, differentiation, and immunoglobulin class switching. Many other immune-competent cells involved in innate and adaptive immunity communicate through this co-stimulatory ligand-receptor dyad. CD40-CD154 interaction is involved in the pathogenesis of numerous inflammatory and autoimmune diseases. While CD40 and CD154 are membrane-bound proteins, their soluble counterparts are generated by proteolytic cleavage or alternative splicing. This review summarises current knowledge about the impact of single nucleotide polymorphisms in the human CD40 gene and compensatory changes in the plasma level of the soluble CD40 receptor (sCD40) isoform in related pro-inflammatory diseases. It discusses regulation patterns of the disintegrin metalloprotease ADAM17 function leading to ectodomain shedding of transmembrane proteins, such as pro-inflammatory adhesion molecules or CD40. The role of sCD40 as a potential biomarker for chronic inflammatory diseases will also be discussed.

Identification and Validation of Immune-Related Genes Diagnostic for Progression of Atherosclerosis and Diabetes

Background

Atherosclerosis and type 2 diabetes mellitus contribute to a large part of cardiovascular events, but the underlying mechanism remains unclear. In this study, we focused on identifying the linking genes of the diagnostic biomarkers and effective therapeutic targets associated with these two diseases.

Methods

The transcriptomic datasets of atherosclerosis and type 2 diabetes mellitus were obtained from the GEO database. Differentially expressed genes analysis was performed by R studio software, and differential analysis including functional enrichment, therapeutic small molecular agents prediction, and protein-protein interaction analysis were applied to the common shared differentially expressed genes. Hub genes were identified and further validated using an independent dataset and clinical samples. Furthermore, we measured the expression correlations, immune cell infiltration, and diagnostic capability of the three key genes.

Results

We screened out 28 up-regulated and six down-regulated common shared differentially expressed genes. Functional enrichment analysis showed that cytokines and immune activation were involved in the development of these two diseases. Six small molecules with the highest absolute enrichment value were identified. Three critical genes (CD4, PLEK, and THY1) were further validated both in validation sets and clinical samples. The gene correlation analysis showed that CD4 was strongly positively correlated with PLEK, and ROC curves confirmed the good discriminatory capacity of CD4 and PLEK in two diseases.

We have established the co-expression network between atherosclerosis lesions progressions and type 2 diabetes mellitus, and identified CD4 and PLEK as key genes in the two diseases, which may facilitate both development of diagnosis and therapeutic strategies.

Identification and Validation of Key Gene Modules and Pathways in Coronary Artery Disease Development and Progression

The development and progression of atherosclerosis represent a chronic process involving complex molecular interactions. Therefore, identifying the potential hub genes and pathways contributing to coronary artery disease (CAD) development is essential for understanding its underlying molecular mechanisms. To this end, we performed transcriptome analysis of peripheral venous blood collected from 100 patients who were divided into four groups according to disease severity, including 27 patients in the atherosclerosis group, 22 patients in the stable angina group, 35 patients in the acute myocardial infarction group, and 16 controls. Weighted gene co-expression network analysis was performed using R programming. Significant module-trait correlations were identified according to module membership and genetic significance. Metascape was used for the functional enrichment of differentially expressed genes between groups, and the hub genes were identified via protein-protein interaction network analysis. The hub genes were further validated by analyzing Gene Expression Omnibus (GSE48060 and GSE141512) datasets. A total of 9,633 messenger ribonucleic acids were detected in three modules, among which the blue module was highly correlated with the Gensini score. The hub genes were significantly enriched in the myeloid leukocyte activation pathway, suggesting its important role in the progression of atherosclerosis. Among these genes, the Mediterranean fever gene (MEFV) may play a key role in the progression of atherosclerosis and CAD severity.

Transcriptomic Profiling Unravels Novel Deregulated Gene Signatures Associated with Acute Myocardial Infarction: A Bioinformatics Approach

Acute myocardial infarction (AMI) is a severe disease with elevated morbidity and mortality rate worldwide. This is attributed to great losses of cardiomyocytes, which can trigger the alteration of gene expression patterns. Although several attempts have been made to assess the AMI biomarkers, to date their role in rescuing myocardial injury remains unclear. Therefore, the current study investigated three independent microarray-based gene expression datasets from AMI patients (n = 85) and their age-sex-matched healthy controls (n = 70), to identify novel gene signatures that might be involved in cardioprotection. The differentially expressed genes (DEGs) were analyzed using 'GEO2R', and weighted gene correlation network analysis (WGCNA) was performed to identify biomarkers/modules. We found 91 DEGs, of which the number of upregulated and downregulated genes were 22 and 5, respectively. Specifically, we found that the deregulated genes such as ADOR-A3, BMP6, VPS8, and GPx3, may be associated with AMI. WGCNA revealed four highly preserved modules among all datasets. The 'Enrichr' unveiled the presence of miR-660 and STAT1, which is known to affect AMI severity. Conclusively, these genes and miRNA might play a crucial role the rescue of cardiomyocytes from severe damage, which could be helpful in developing appropriate therapeutic strategies for the management of AMI.

Relationship between the Soluble F11 Receptor and Annexin A5 in African Americans Patients with Type-2 Diabetes Mellitus

Type 2 diabetes mellitus (T2DM) is characterized by endothelial dysfunction, increased thrombogenicity, and inflammation. The soluble human F11 receptor (sF11R) and annexin A5 (ANXA5) play crucial roles in inflammatory thrombosis and atherosclerosis. We examined the relationship between circulating sF11R and ANXA5 and their impact on endothelial function. The study included 125 patients with T2DM. Plasma levels of sF11R and ANXA5 were quantified by ELISA. Microvascular function was assessed using the vascular reactivity index (VRI). Large artery stiffness was assessed by carotid-femoral pulse wave velocity (PWV). Carotid intima-media thickness (CIMT) was assessed by B-mode ultrasound imaging. The mean age of patients in the study was 59.7 ± 7.8 years, 78% had hypertension, 76% had dyslipidemia, and 12% had CKD. sF11R correlated positively with ANXA5 levels (β = 0.250, p = 0.005), and correlated inversely with VRI and total nitic oxide (NO), (β = −0.201, p = 0.024; β = −0.357, p = 0.0001, respectively). Multivariate regression analysis revealed that sF11R was independently associated with ANXA5 in the total population and in patients with HbA1c > 6.5% (β = 0.366, p = 0.007; β = 0.425, p = 0.0001, respectively). sF11R and ANXA5 were not associated with vascular outcome, suggesting that they may not be reliable markers of vascular dysfunction in diabetes. The clinical significance of sF11R/ANXA5 association in diabetes warrants further investigation in a larger population.

LncRNA LYPLAL1-DT screening from type 2 diabetes with macrovascular complication contributes protective effects on human umbilical vein endothelial cells via regulating the miR-204-5p/SIRT1 axis

Long noncoding RNAs (lncRNAs) are involved in diabetes related diseases. However, the role of lncRNAs in the pathogenesis of type 2 diabetes with macrovascular complication (DMC) has seldomly been recognized. This study screened lncRNA profiles of leukocytes from DMC patients and explored protective role of lncRNA LYPLAL1-DT in endothelial cells (EC) under high glucose (HG) and inflammatory conditions (IS). Between DMC and healthy controls, 477 differential expression lncRNAs (DE-lncRNAs) were identified. The enrichment and pathway analysis showed that most of the DE-lncRNAs belonged to inflammatory, metabolic, and vascular diseases. A total of 12 lncRNAs was validated as significant DE-lncRNAs in expanding cohorts. Furthermore, these DE-lncRNAs were shown to be significantly related to hypoxia, HG, and IS in EC, especially lncRNA LYPLAL1-DT. LYPLAL1-DT overexpression results in the promotion of the proliferation, and migration of EC, as well as an elevation of autophagy. Overexpressed LYPLAL1-DT reduces the adhesion of monocytes to EC, boosts anti-inflammation, and suppresses inflammatory molecules secreted in the medium. Mechanistically, LYPLAL1-DT acts as competing endogenous RNA (ceRNA) by downregulating miR-204-5p, therefore enhancing SIRT1 and protecting EC autophagy function; thus, alleviating apoptosis. Finally, exosome sequencing revealed LYPLAL1-DT expression was 4 times lower in DMC cells than in healthy samples. In general, we identified LYPLAL1-DT having protective effects on EC as ceRNA mediated through the miR-204-5p/SIRT1 pathway. Therefore, it inhibits the autophagy of EC as well as modulating systemic inflammation. This approach could be regarded as a new potential therapeutic target in DMC.

Weighted Co-Expression Network Analysis Identifies RNF181 as a Causal Gene of Coronary Artery Disease

Background: Coronary artery disease (CAD) exerts a global challenge to public health. Genetic heritability is one of the most vital contributing factors in the pathophysiology of CAD. Co-expression network analysis is an applicable and robust method for the interpretation of biological interaction from microarray data. Previous CAD studies have focused on peripheral blood samples since the processes of CAD may vary from tissue to blood. It is therefore necessary to find biomarkers for CAD in heart tissues; their association also requires further illustration.

Materials and Methods: To filter for causal genes, an analysis of microarray expression profiles, GSE12504 and GSE22253, was performed with weighted gene co-expression network analysis (WGCNA). Co-expression modules were constructed after batch effect removal and data normalization. The results showed that 7 co-expression modules with 8,525 genes and 1,210 differentially expressed genes (DEGs) were identified. Furthermore, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were conducted. Four major pathways in CAD tissue and hub genes were addressed in the Hybrid Mouse Diversity Panel (HMDP) and Human Protein Atlas (HPA), and isoproterenol (ISO)/doxycycline (DOX)-induced heart toxicity models were used to validate the hub genes. Lastly, the hub genes and risk variants were verified in the CAD cohort and in genome-wide association studies (GWAS).

Results: The results showed that RNF181 and eight other hub genes are perturbed during CAD in heart tissues. Additionally, the expression of RNF181 was validated using RT-PCR and immunohistochemistry (IHC) staining in two cardiotoxicity mouse models. The association was further verified in the CAD patient cohort and in GWAS.

Conclusion: Our findings illustrated for the first time that the E3 ubiquitination ligase protein RNF181 may serve as a potential biomarker in CAD, but further in vivo validation is warranted.

Multi-Omic Approaches to Identify Genetic Factors in Metabolic Syndrome

Metabolic syndrome (MetS) is a highly heritable disease and a major public health burden worldwide. MetS diagnosis criteria are met by the simultaneous presence of any three of the following: high triglycerides, low HDL/high LDL cholesterol, insulin resistance, hypertension, and central obesity. These diseases act synergistically in people suffering from MetS and dramatically increase risk of morbidity and mortality due to stroke and cardiovascular disease, as well as certain cancers. Each of these component features is itself a complex disease, as is MetS. As a genetically complex disease, genetic risk factors for MetS are numerous, but not very powerful individually, often requiring specific environmental stressors for the disease to manifest. When taken together, all sequence variants that contribute to MetS disease risk explain only a fraction of the heritable variance, suggesting additional, novel loci have yet to be discovered. In this article, we will give a brief overview on the genetic concepts needed to interpret genome-wide association studies (GWAS) and quantitative trait locus (QTL) data, summarize the state of the field of MetS physiological genomics, and to introduce tools and resources that can be used by the physiologist to integrate genomics into their own research on MetS and any of its component features. There is a wealth of phenotypic and molecular data in animal models and humans that can be leveraged as outlined in this article. Integrating these multi-omic QTL data for complex diseases such as MetS provides a means to unravel the pathways and mechanisms leading to complex disease and promise for novel treatments. © 2022 American Physiological Society. Compr Physiol 12:1-40, 2022.

Ischemic heart injury leads to HIF1-dependent differential splicing of CaMK2γ

Ischemic heart disease is a leading cause of heart failure and hypoxia inducible factor 1 (HIF1) is a key transcription factor in the response to hypoxic injury. Our lab has developed a mouse model in which a mutated, oxygen-stable form of HIF1α (HIF-PPN) can be inducibly expressed in cardiomyocytes. We observed rapid cardiac dilation and loss of contractility in these mice due to lower expression of excitation-contraction coupling genes and reduced calcium flux. As alternative splicing plays an underappreciated role in transcriptional regulation, we used RNA sequencing to search for splicing changes in calcium-handling genes of HIF-PPN hearts and compared them to previous sequencing data from a model of myocardial infarction (MI) to select for transcripts that are modified in a pathological setting. We found overlap between genes differentially expressed in HIF-PPN and post-MI mice (54/131 genes upregulated in HIF-PPN hearts at 1 day and/or 3 days post-MI, and 45/78 downregulated), as well as changes in alternative splicing. Interestingly, calcium/calmodulin dependent protein kinase II, gamma (CAMK2G) was alternatively spliced in both settings, with variant 1 (v1) substantially decreased compared to variants 2 (v2) and 3 (v3). These findings were also replicated in vitro when cells were transfected with HIF-PPN or exposed to hypoxia. Further analysis of CAMK2γ protein abundance revealed only v1 was detectable and substantially decreased up to 7 days post-MI. Rbfox1, a splicing factor of CAMK2G, was also decreased in HIF-PPN and post-MI hearts. Subcellular fractionation showed CAMK2γ v1 was found in the nuclear and cytoplasmic fractions, and abundance decreased in both fractions post-MI. Chromatin immunoprecipitation analysis of HIF1 in post-MI hearts also demonstrated direct HIF1 binding to CAMK2G. CaMK2 is a key transducer of calcium signals in both physiological and pathological settings. The predominantly expressed isoform in the heart, CaMK2δ, has been extensively studied in cardiac injury, but the specific role of CaMK2γ is not well defined. Our data suggest that loss of CaMK2γ after MI is HIF1-dependent and may play an important role in the heart's calcium signaling and transcriptional response to hypoxia.

Circular RNA circitga7 accelerates glioma progression via miR-34a-5p/VEGFA axis

Circular RNAs (circRNAs) are a group of noncoding RNAs derived from back-splicing events. CircRNA is reported to be involved in various tumor progressions, including glioma. Although there are a few reports of circular RNAs participating in gliomas, it is still unclear whether circular RNAs regulate the occurrence of gliomas. In our research, we found that the expression of circITGA7 in glioma tissues and glioma cells increased significantly. Knocking down circITGA7 can significantly inhibit the proliferation of glioma cells and reduce cell metastasis. Through analysis and dual-luciferase report assay, we found that circITGA7 acts as a sponge for miR-34a-5p targeting VEGFA in glioma. Our study showed that circITGA7 regulates the proliferation and metastasis of glioma cell lines (SW1783&U373) by regulating the miR-34a-5p/VEGFA pathway. In conclusion, our study revealed a regulatory loop for the circITGA7/miR-34a-5p/VEGFA axis to regulate glioma development.

Exome sequencing in patient-parent trios suggests new candidate genes for early-onset primary sclerosing cholangitis

BACKGROUND & AIMS

Primary sclerosing cholangitis (PSC) is a rare bile duct disease strongly associated with inflammatory bowel disease (IBD). Whole-exome sequencing (WES) has contributed to understanding the molecular basis of very early-onset IBD, but rare protein-altering genetic variants have not been identified for early-onset PSC. We performed WES in patients diagnosed with PSC ≤ 12 years to investigate the contribution of rare genetic variants to early-onset PSC.

METHODS

In this multicentre study, WES was performed on 87 DNA samples from 29 patient-parent trios with early-onset PSC. We selected rare (minor allele frequency < 2%) coding and splice-site variants that matched recessive (homozygous and compound heterozygous variants) and dominant (de novo) inheritance in the index patients. Variant pathogenicity was predicted by an in-house developed algorithm (GAVIN), and PSC-relevant variants were selected using gene expression data and gene function.

RESULTS

In 22 of 29 trios we identified at least 1 possibly pathogenic variant. We prioritized 36 genes, harbouring a total of 54 variants with predicted pathogenic effects. In 18 genes, we identified 36 compound heterozygous variants, whereas in the other 18 genes we identified 18 de novo variants. Twelve of 36 candidate risk genes are known to play a role in transmembrane transport, adaptive and innate immunity, and epithelial barrier function.

CONCLUSIONS

The 36 candidate genes for early-onset PSC need further verification in other patient cohorts and evaluation of gene function before a causal role can be attributed to its variants.

Identifying RBM47, HCK, CD53, TYROBP, and HAVCR2 as Hub Genes in Advanced Atherosclerotic Plaques by Network-Based Analysis and Validation

Background: Atherosclerotic cardiovascular diseases accounted for a quarter of global deaths. Most of these fatal diseases like coronary atherosclerotic disease (CAD) and stroke occur in the advanced stage of atherosclerosis, during which candidate therapeutic targets have not been fully established. This study aims to identify hub genes and possible regulatory targets involved in treatment of advanced atherosclerotic plaques.

Material/Methods: Microarray dataset GSE43292 and GSE28829, both containing advanced atherosclerotic plaques group and early lesions group, were obtained from the Gene Expression Omnibus database. Weighted gene co-expression network analysis (WGCNA) was conducted to identify advanced plaque-related modules. Module conservation analysis was applied to assess the similarity of advanced plaque-related modules between GSE43292 and GSE28829. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of these modules were performed by Metascape. Differentially expressed genes (DEGs) were mapped into advanced plaque-related modules and module membership values of DEGs in each module were calculated to identify hub genes. Hub genes were further validated for expression in atherosclerotic samples, for distinguishing capacity of CAD and for potential functions in advanced atherosclerosis.

Results: The lightgreen module (MElightgreen) in GSE43292 and the brown module (MEbrown) in GSE28829 were identified as advanced plaque-related modules. Conservation analysis of these two modules showed high similarity. GO and KEGG enrichment analysis revealed that genes in both MElightgreen and MEbrown were enriched in immune cell activation, secretory granules, cytokine activity, and immunoinflammatory signaling. RBM47, HCK, CD53, TYROBP, and HAVCR2 were identified as common hub genes, which were validated to be upregulated in advanced atherosclerotic plaques, to well distinguish CAD patients from non-CAD people and to regulate immune cell function-related mechanisms in advanced atherosclerosis.

Conclusions: We have identified RBM47, HCK, CD53, TYROBP, and HAVCR2 as immune-responsive hub genes related to advanced plaques, which may provide potential intervention targets to treat advanced atherosclerotic plaques.