Identifying Significant Genes and Functionally Enriched Pathways in Familial Hypercholesterolemia Using Integrated Gene Co-Expression Network Analysis

Gene expression profiling in elderly patients with familial hypercholesterolemia with and without coronary heart disease

BACKGROUND AND AIMS

Elderly familial hypercholesterolemia (FH) patients are at high risk of coronary heart disease (CHD) due to high cholesterol burden and late onset of effective cholesterol-lowering therapies. A subset of these individuals remains free from any CHD event, indicating the potential presence of protective factors. Identifying possible cardioprotective gene expression profiles could contribute to our understanding of CHD prevention and future preventive treatment. Therefore, this study aimed to investigate gene expression profiles in elderly event-free FH patients.

METHODS

Expression of 773 genes was analysed using the Nanostring Metabolic Pathways Panel, in peripheral blood mononuclear cells (PBMCs) from FH patients ≥65 years without CHD (FH event-free, n = 44) and with CHD (FH CHD, n = 39), and from healthy controls ≥70 years (n = 39).

RESULTS

None of the genes were differentially expressed between FH patients with and without CHD after adjusting for multiple testing. However, at nominal p < 0.05, we found 36 (5%) differentially expressed genes (DEGs) between the two FH groups, mainly related to lipid metabolism (e.g. higher expression of ABCA1 and ABCG1 in FH event-free) and immune responses (e.g. lower expression of STAT1 and STAT3 in FH event-free). When comparing FH patients to controls, the event-free group had fewer DEGs than the CHD group; 147 (19%) and 219 (28%) DEGs, respectively.

CONCLUSIONS

Elderly event-free FH patients displayed a different PBMC gene expression profile compared to FH patients with CHD. Differences in gene expression compared to healthy controls were more pronounced in the CHD group, indicating a less atherogenic gene expression profile in event-free individuals. Overall, identification of cardioprotective factors could lead to future therapeutic targets.

Identification of Plasma Exosomes hsa_circ_0001360 and hsa_circ_0000038 as Key Biomarkers of Coronary Heart Disease

Background

Coronary heart disease (CHD) is the leading cause of death and disability worldwide. Accumulating evidence reveals that atherosclerosis (AS), characterized by systemic, chronic, and multifocal disease, and is the primary pathological basis of cardiovascular diseases, including CHD. However, the molecular underpinnings of CHD are still far from well understood. Our study attempted to identify aberrant plasma exosome-derived circRNAs and key exosomal circRNA biomarkers for CHD.

Methods

The expression profiles of mRNAs, circRNAs, and lncRNAs in the blood exosomes of CHD patients and healthy controls were obtained from the exoRBase database. The corresponding miRNAs of the differentially expressed mRNAs, circRNAs, and lncRNAs were predicted via ENCORI and the miRcode database. LncRNAs/circRNAs and mRNAs with the cotargeted miRNAs were selected to construct an interaction network. Multiple machine learning algorithms have been used to explore potential biomarkers, followed by verification in patients with CHD using real-time quantitative reverse transcription-polymerase chain reaction (RT-qPCR).

Results

Based on the cutoff criterion of P < 0.05, we identified 85 differentially expressed circRNAs (4 upregulated and 81 downregulated), 43 differentially expressed lncRNAs (24 upregulated and 19 downregulated), and 312 differentially expressed mRNAs (55 upregulated and 257 downregulated). Functional enrichment analysis revealed that the differentially expressed mRNAs were involved mainly in neutrophil extracellular trap (NET) formation and the nucleotide-binding oligomerization domain- (NOD-) like receptor signaling pathway. Further analysis revealed that the DEGs in the circRNA/lncRNA-miRNA-mRNA interaction network were closely related to lipid and atherosclerotic signaling pathways. Hsa_circ_0001360 and hsa_circ_0000038 were identified as potential biomarkers for CHD based on three machine learning algorithms. The relative expression levels of hsa_circ_0001360 and hsa_circ_0000038 were significantly altered in plasma exosomes from patients with CHD. ROC curve analysis revealed that the areas under the curve (AUCs) were 0.860, 0.870, and 0.940 for hsa_circ_0001360, hsa_circ_0000038, and the two-gene combination, respectively.

Conclusion

The circRNA/lncRNA-miRNA-mRNA interaction network might help to elucidate the pathogenesis of CHD. Hsa_circ_0001360 combined with hsa_circ_0000038 might be an important diagnostic biomarker.

Network pharmacology based anti-diabetic attributes of bioactive compounds from Ocimum gratissimum L. through computational approach

The present research was framed to determine the key compounds present in the plant Ocimum gratissimum L. targeting protein molecules of Diabetes Mellitus (DM) by employing In-silico approaches. Phytochemicals previously reported to be present in this herb were collated through literature survey and public phytochemical databases, and their probable targets were anticipated using BindingDB (p ≥ 0.7). STRING and KEGG pathway databases were employed for pathway enrichment analysis. Homology modelling was executed to elucidate the structures of therapeutic targets. Further, Phytocompounds from O. gratissimum were subjected for docking with four therapeutic targets of DM by using AutoDock vina through POAP pipeline implementation. 30 compounds were predicted to target 136 protein molecules including aldose reductase, DPP4, alpha-amylase, and alpha-glucosidase. Neuroactive ligand-receptor interaction, MAPK, PI3K-Akt, starch and insulin resistance were predicted to have potentially modulation by phytocompounds. Based on the phytocompound's binding score with the four targets of DM, Rutin scored the lowest binding energy (-11 kcal/mol) with Aldose reductase by forming 17 intermolecular interactions. In conclusion, based on the network and binding score, phytocompounds from O. gratissimum have a synergistic and considerable effect in the management of DM via multi-compound, multi-target, and multi-pathway mechanisms.

Rare variant burden analysis from exomes of three consanguineous families reveals LILRB1 and PRSS3 as potential key proteins in inflammatory bowel disease pathogenesis

Background

Inflammatory bowel disease (IBD) is a chronic autoimmune disorder characterized by severe inflammation and mucosal destruction of the intestine. The specific, complex molecular processes underlying IBD pathogenesis are not well understood. Therefore, this study is aimed at identifying and uncovering the role of key genetic factors in IBD.

Method

The whole exome sequences (WESs) of three consanguineous Saudi families having many siblings with IBD were analyzed to discover the causal genetic defect. Then, we used a combination of artificial intelligence approaches, such as functional enrichment analysis using immune pathways and a set of computational functional validation tools for gene expression, immune cell expression analyses, phenotype aggregation, and the system biology of innate immunity, to highlight potential IBD genes that play an important role in its pathobiology.

Results

Our findings have shown a causal group of extremely rare variants in the LILRB1 (Q53L, Y99N, W351G, D365A, and Q376H) and PRSS3 (F4L and V25I) genes in IBD-affected siblings. Findings from amino acids in conserved domains, tertiary-level structural deviations, and stability analysis have confirmed that these variants have a negative impact on structural features in the corresponding proteins. Intensive computational structural analysis shows that both genes have very high expression in the gastrointestinal tract and immune organs and are involved in a variety of innate immune system pathways. Since the innate immune system detects microbial infections, any defect in this system could lead to immune functional impairment contributing to IBD.

Conclusion

The present study proposes a novel strategy for unraveling the complex genetic architecture of IBD by integrating WES data of familial cases, with computational analysis.

Potential Biomarkers for Parkinson Disease from Functional Enrichment and Bioinformatic Analysis of Global Gene Expression Patterns of Blood and Substantia Nigra Tissues

The Parkinson disease (PD) is the second most common neurodegenerative disorder affecting the central nervous system and motor functions. The biological complexity of PD is yet to reveal potential targets for intervention or to slow the disease severity. Therefore, this study aimed to compare the fidelity of blood to substantia nigra (SN) tissue gene expression from PD patients to provide a systematic approach to predict role of the key genes of PD pathobiology. Differentially expressed genes (DEGs) from multiple microarray data sets of PD blood and SN tissue from GEO database are identified. Using the theoretical network approach and variety of bioinformatic tools, we prioritized the key genes from DEGs. A total of 540 and 1024 DEGs were identified in blood and SN tissue samples, respectively. Functional pathways closely related to PD such as ERK1 and ERK2 cascades, mitogen-activated protein kinase (MAPK) signaling, Wnt, nuclear factor-κB (NF-κB), and PI3K-Akt signaling were observed by enrichment analysis. Expression patterns of 13 DEGs were similar in both blood and SN tissues. Comprehensive network topological analysis and gene regulatory networks identified additional 10 DEGs functionally connected with molecular mechanisms of PD through the mammalian target of rapamycin (mTOR), autophagy, and AMP-activated protein kinase (AMPK) signaling pathways. Potential drug molecules were identified by chemical-protein network and drug prediction analysis. These potential candidates can be further validated in vitro/in vivo to be used as biomarkers and/or novel drug targets for the PD pathology and/or to arrest or delay the neurodegeneration over the years, respectively.

Exome Sequencing Identifies the Extremely Rare ITGAV and FN1 Variants in Early Onset Inflammatory Bowel Disease Patients

BACKGROUND

Molecular diagnosis of early onset inflammatory bowel disease (IBD) is very important for adopting suitable treatment strategies. Owing to the sparse data available, this study aims to identify the molecular basis of early onset IBD in Arab patients.

METHODS

A consanguineous Arab family with monozygotic twins presenting early onset IBD was screened by whole exome sequencing (WES). The variants functional characterization was performed by a series of computational biology methods. The IBD variants were further screened in in-house whole exome data of 100 Saudi cohorts ensure their rare prevalence in the population.

RESULTS

Genetic screening has identified the digenic autosomal recessive mode of inheritance of ITGAV (G58V) and FN1 (G313V) variants in IBD twins with early onset IBD. Findings from pathogenicity predictions, stability and molecular dynamics have confirmed the deleterious nature of both variants on structural features of the corresponding proteins. Functional biology data suggested that both genes show abundant expression in gastrointestinal tract and immune organs, involved in immune cell restriction, regulation of different immune related pathways. Data from knockout mouse models for ITGAV gene has revealed that the dysregulated expression of this gene impacts intestinal immune homeostasis. The defective ITGAV and FN1 involved in integrin pathway, are likely to induce intestinal inflammation by disturbing immune homeostasis.

CONCLUSIONS

Our findings provide novel insights into the molecular etiology of pediatric onset IBD and may likely pave way in developing genomic medicine.