Transmembrane signaling molecules play a key role in the pathogenesis of IgA nephropathy: a weighted gene co-expression network analysis study

Transcriptome meta-analysis and validation to discovery of hub genes and pathways in focal and segmental glomerulosclerosis

BACKGROUND

Focal segmental glomerulosclerosis (FSGS), a histologic pattern of injury in the glomerulus, is one of the leading glomerular causes of end-stage renal disease (ESRD) worldwide. Despite extensive research, the underlying biological alterations causing FSGS remain poorly understood. Studying variations in gene expression profiles offers a promising approach to gaining a comprehensive understanding of FSGS molecular pathogenicity and identifying key elements as potential therapeutic targets. This work is a meta-analysis of gene expression profiles from glomerular samples of FSGS patients. The main aims of this study are to establish a consensus list of differentially expressed genes in FSGS, validate these findings, understand the disease's pathogenicity, and identify novel therapeutic targets.

METHODS

After a thorough search in the GEO database and subsequent quality control assessments, seven gene expression datasets were selected for the meta-analysis: GSE47183 (GPL14663), GSE47183 (GPL11670), GSE99340, GSE108109, GSE121233, GSE129973, and GSE104948. The random effect size method was applied to identify differentially expressed genes (meta-DEGs), which were then used to construct a regulatory network (STRING, MiRTarBase, and TRRUST) and perform various pathway enrichment analyses. The expression levels of several meta-DEGs, specifically ADAMTS1, PF4, EGR1, and EGF, known as angiogenesis regulators, were analyzed using quantitative reverse transcription polymerase chain reaction (RT-qPCR).

RESULTS

The identified 2,898 meta-DEGs, including 665 downregulated and 669 upregulated genes, were subjected to various analyses. A co-regulatory network comprising 2,859 DEGs, 2,688 microRNAs (miRNAs), and 374 transcription factors (TFs) was constructed, and the top molecules in the network were identified based on degree centrality. Part of the pathway enrichment analysis revealed significant disruption in the angiogenesis regulatory pathways in the FSGS kidney. The RT-qPCR results confirmed an imbalance in angiogenesis pathways by demonstrating the differential expression levels of ADAMTS1 and EGR1, two key angiogenesis regulators, in the FSGS condition.

CONCLUSION

In addition to presenting a consensus list of differentially expressed genes in FSGS, this meta-analysis identified significant distortions in angiogenesis-related pathways and factors in the FSGS kidney. Targeting these factors may offer a viable strategy to impede the progression of FSGS.

Metabolome panels as potential noninvasive biomarkers for primary glomerulonephritis sub-types: meta-analysis of profiling metabolomics studies

Primary glomerulonephritis diseases (PGDs) are known as the top causes of chronic kidney disease worldwide. Renal biopsy, an invasive method, is the main approach to diagnose PGDs. Studying the metabolome profiles of kidney diseases is an inclusive approach to identify the disease's underlying pathways and discover novel non-invasive biomarkers. So far, different experiments have explored the metabolome profiles in different PGDs, but the inconsistencies might hinder their clinical translations. The main goal of this meta-analysis study was to achieve consensus panels of dysregulated metabolites in PGD sub-types. The PGDs-related metabolome profiles from urine samples in humans were selected in a comprehensive search. Amanida package in R software was utilized for performing the meta-analysis. Through sub-type analyses, the consensus list of metabolites in each category was obtained. To identify the most affected pathways, functional enrichment analysis was performed. Also, a gene-metabolite network was constructed to identify the key metabolites and their connected proteins. After a vigorous search, among the 11 selected studies (15 metabolite profiles), 270 dysregulated metabolites were recognized in urine of 1154 PGDs and control samples. Through sub-type analyses by Amanida package, the consensus list of metabolites in each category was obtained. Top dysregulated metabolites (vote score of ≥ 4 or ≤ - 4) in PGDs urines were selected as main panel of meta-metabolites including glucose, leucine, choline, betaine, dimethylamine, fumaric acid, citric acid, 3-hydroxyisovaleric acid, pyruvic acid, isobutyric acid, and hippuric acid. The enrichment analyses results revealed the involvement of different biological pathways such as the TCA cycle and amino acid metabolisms in the pathogenesis of PGDs. The constructed metabolite-gene interaction network revealed the high centralities of several metabolites, including pyruvic acid, leucine, and choline. The identified metabolite panels could shed a light on the underlying pathological pathways and be considered as non-invasive biomarkers for the diagnosis of PGD sub-types.

Clear Cell Renal Cell Carcinoma: A Comprehensive in silico Study in Searching for Therapeutic Targets

INTRODUCTION

Clear cell renal cell carcinoma (ccRCC) is recognized as one of the leading causes of illness and death worldwide. Understanding the molecular mechanisms in ccRCC pathogenesis is crucial for discovering novel therapeutic targets and developing efficient drugs. With the application of a comprehensive in silico analysis of the ccRCC-related array sets, the main objective of this study was to discover the top molecules and pathways in the pathogenesis of this cancer.

METHODS

ccRCC microarray datasets were downloaded from the Gene Expression Omnibus database, and after quality checking, normalization, and analysis using the Limma algorithm, differentially expressed genes (DEGs) were identified, considering the adjusted p value <0.049. The intensity values of the identified DEGs were introduced to the Weighted Gene Co-Expression Network Analysis (WGCNA) algorithm to construct co-expression modules. Functional enrichment analyses were performed using the DEGs in the disease-correlated module, and hub genes were identified among the top genes in a protein-protein interaction network and the disease most correlated module. The expression analysis of hub genes was done by utilizing GEPIA, and the GSCA server was used to compare the expression patterns of hub genes in ccRCC and other cancers. DGIdb database was utilized to identify the hub gene-related drugs.

RESULTS

Three datasets, including GSE11151, GSE12606, and GSE36897, were retrieved, merged, normalized, and analyzed. Using WGCNA, the DEGs were clustered into eight different modules. Translocation of ZAP-70 to immunological synapse, endosomal/vacuolar pathway, cell surface interactions at the vascular wall, and immune-related pathways were the topmost enriched terms for the ccRCC-correlated DEGs. Twelve genes including PTPRC, ITGAM, TLR2, CD86, PLEK, TYROBP, ITGB2, RAC2, CSF1R, CCR5, CCL5, and LCP2 were introduced as hub genes. All the 12 hub genes were upregulated in ccRCC samples and showed a positive correlation with the infiltration of different immune cells. According to the DGIdb database, 127 drugs, including tyrosine kinase inhibitors, glucocorticoids, and chemotaxis targeting molecules, were identified to interact with the hub genes.

CONCLUSION

By utilizing an integrative bioinformatics approach, this experiment shed light on the underlying pathways in the pathogenesis of ccRCC and introduced several potential therapeutic targets for repurposing or developing novel drugs for an efficient treatment of this cancer. Our next step would be to assess the gene expression profiles of the identified hubs in different cell populations in the tumor microenvironment.

Integrative analyses of biomarkers and pathways for diabetic nephropathy

Background: Diabetic nephropathy (DN) is a widespread diabetic complication and a major cause of terminal kidney disease. There is no doubt that DN is a chronic disease that imposes substantial health and economic burdens on the world's populations. By now, several important and exciting advances have been made in research on etiopathogenesis. Therefore, the genetic mechanisms underlying these effects remain unknown. Methods: The GSE30122, GSE30528, and GSE30529 microarray datasets were downloaded from the Gene Expression Omnibus database (GEO). Analyses of differentially expressed genes (DEGs), enrichment of gene ontology (GO), the Kyoto Encyclopedia of Genes and Genomes (KEGG), and gene set enrichment analysis (GSEA) were performed. Protein-protein interaction (PPI) network construction was completed by the STRING database. Hub genes were identified by Cytoscape software, and common hub genes were identified by taking intersection sets. The diagnostic value of common hub genes was then predicted in the GSE30529 and GSE30528 datasets. Further analysis was carried out on the modules to identify transcription factors and miRNA networks. As well, a comparative toxicogenomics database was used to assess interactions between potential key genes and diseases associated upstream of DN. Results: Samples from 19 DNs and 50 normal controls were identified in the GSE30122 dataset. 86 upregulated genes and 34 downregulated genes (a total of 120 DEGs). GO analysis showed significant enrichment in humoral immune response, protein activation cascade, complement activation, extracellular matrix, glycosaminoglycan binding, and antigen binding. KEGG analysis showed significant enrichment in complement and coagulation cascades, phagosomes, the Rap1 signaling pathway, the PI3K-Akt signaling pathway, and infection. GSEA was mainly enriched in the TYROBP causal network, the inflammatory response pathway, chemokine receptor binding, the interferon signaling pathway, ECM receptor interaction, and the integrin 1 pathway. Meanwhile, mRNA-miRNA and mRNA-TF networks were constructed for common hub genes. Nine pivotal genes were identified by taking the intersection. After validating the expression differences and diagnostic values of the GSE30528 and GSE30529 datasets, eight pivotal genes (TYROBP, ITGB2, CD53, IL10RA, LAPTM5, CD48, C1QA, and IRF8) were finally identified as having diagnostic values. Conclusion: Pathway enrichment analysis scores provide insight into the genetic phenotype and may propose molecular mechanisms of DN. The target genes TYROBP, ITGB2, CD53, IL10RA, LAPTM5, CD48, C1QA, and IRF8 are promising new targets for DN. SPI1, HIF1A, STAT1, KLF5, RUNX1, MBD1, SP1, and WT1 may be involved in the regulatory mechanisms of DN development. Our study may provide a potential biomarker or therapeutic locus for the study of DN.

Comparative co-expression analysis of RNA-Seq transcriptome revealing key genes, miRNA and transcription factor in distinct metabolic pathways in diabetic nerve, eye, and kidney disease

Diabetes and its related complications are associated with long term damage and failure of various organ systems. The microvascular complications of diabetes considered in this study are diabetic retinopathy, diabetic neuropathy, and diabetic nephropathy. The aim is to identify the weighted co-expressed and differentially expressed genes (DEGs), major pathways, and their miRNA, transcription factors (TFs) and drugs interacting in all the three conditions. The primary goal is to identify vital DEGs in all the three conditions. The overlapped five genes (AKT1, NFKB1, MAPK3, PDPK1, and TNF) from the DEGs and the co-expressed genes were defined as key genes, which differentially expressed in all the three cases. Then the protein-protein interaction network and gene set linkage analysis (GSLA) of key genes was performed. GSLA, gene ontology, and pathway enrichment analysis of the key genes elucidates nine major pathways in diabetes. Subsequently, we constructed the miRNA-gene and transcription factor-gene regulatory network of the five gene of interest in the nine major pathways were studied. hsa-mir-34a-5p, a major miRNA that interacted with all the five genes. RELA, FOXO3, PDX1, and SREBF1 were the TFs interacting with the major five gene of interest. Finally, drug-gene interaction network elucidates five potential drugs to treat the genes of interest. This research reveals biomarker genes, miRNA, TFs, and therapeutic drugs in the key signaling pathways, which may help us, understand the processes of all three secondary microvascular problems and aid in disease detection and management.

A Bioinformatic Approach Based on Systems Biology to Determine the Effects of SARS-CoV-2 Infection in Patients with Hypertrophic Cardiomyopathy

Recently, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of coronavirus disease 2019 (COVID-19), has infected millions of individuals worldwide. While COVID-19 generally affects the lungs, it also damages other organs, including those of the cardiovascular system. Hypertrophic cardiomyopathy (HCM) is a common genetic cardiovascular disorder. Studies have shown that HCM patients with COVID-19 have a higher mortality rate; however, the reason for this phenomenon is not yet elucidated. Herein, we conducted transcriptomic analyses to identify shared biomarkers between HCM and COVID-19 to bridge this knowledge gap. Differentially expressed genes (DEGs) were obtained using the Gene Expression Omnibus ribonucleic acid (RNA) sequencing datasets, GSE147507 and GSE89714, to identify shared pathways and potential drug candidates. We discovered 30 DEGs that were common between these two datasets. Using a combination of statistical and biological tools, protein-protein interactions were constructed in response to these findings to support hub genes and modules. We discovered that HCM is linked to COVID-19 progression based on a functional analysis under ontology terms. Based on the DEGs identified from the datasets, a coregulatory network of transcription factors, genes, proteins, and microRNAs was also discovered. Lastly, our research suggests that the potential drugs we identified might be helpful for COVID-19 therapy.

Screening of the Key Genes and Signalling Pathways for Diabetic Nephropathy Using Bioinformatics Analysis

BACKGROUND

This study aimed to identify biological markers for diabetic nephropathy (DN) and explore their underlying mechanisms.

METHODS

Four datasets, GSE30528, GSE47183, GSE104948, and GSE96804, were downloaded from the Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) were identified using the "limma" package, and the "RobustRankAggreg" package was used to screen the overlapping DEGs. The hub genes were identified using cytoHubba of Cytoscape. Logistic regression analysis was used to further analyse the hub genes, followed by receiver operating characteristic (ROC) curve analysis to predict the diagnostic effectiveness of the hub genes. Correlation analysis and enrichment analysis of the hub genes were performed to identify the potential functions of the hub genes involved in DN.

RESULTS

In total, 55 DEGs, including 38 upregulated and 17 downregulated genes, were identified from the three datasets. Four hub genes (FN1, CD44, C1QB, and C1QA) were screened out by the "UpSetR" package, and FN1 was identified as a key gene for DN by logistic regression analysis. Correlation analysis and enrichment analysis showed that FN1 was positively correlated with four genes (COL6A3, COL1A2, THBS2, and CD44) and with the development of DN through the extracellular matrix (ECM)-receptor interaction pathway.

CONCLUSIONS

We identified four candidate genes: FN1, C1QA, C1QB, and CD44. On further investigating the biological functions of FN1, we showed that FN1 was positively correlated with THBS2, COL1A2, COL6A3, and CD44 and involved in the development of DN through the ECM-receptor interaction pathway. THBS2, COL1A2, COL6A3, and CD44 may be novel biomarkers and target therapeutic candidates for DN.