Identification of key genes, pathways and potential therapeutic agents for IgA nephropathy using an integrated bioinformatics analysis

Mendelian randomization analysis reveals causal relationships between circulating cell traits and renal disorders

Purpose

The aim of this study was to investigate the causal relationships between circulating cell traits and risk of renal disorders.

Methods

We applied a comprehensive two-sample Mendelian randomization (MR) analysis. Single nucleotide polymorphisms (SNPs) from publicly available genome-wide association studies (GWAS) databases were utilized. Genetically predicted instrumental variables of human blood cell traits were extracted from Blood Cell Consortium (BCX) while data on renal diseases was obtained from Finngen consortium. The primary MR analysis was conducted using the inverse variance weighted (IVW) method, with the weighted median (WM) and MR-Egger models used as additional methods. Sensitivity analyses, including MR-PRESSO, radial regression and MR-Egger intercept were conducted to detect outliers and assess horizontal pleiotropy. We further utilized the leave-one-out analysis to assess the robustness of the results. Causal associations were considered significant based on false rate correction (FDR), specifically when the IVW method provided a pFDR < 0.05.

Results

Our results demonstrated that both white blood cell (WBC) count (OR = 1.50, 95% CI = 1.10-2.06, pFDR = 0.033, pIVW = 0.011) and lymphocyte count (OR = 1.50, 95% CI = 1.13-1.98, pFDR = 0.027, pIVW = 0.005) were causally associated with a higher risk of IgA nephropathy. Furthermore, WBC count was identified as a significant genetic risk factor for renal malignant neoplasms (OR = 1.23, 95% CI = 1.06-1.43, pFDR = 0.041, pIVW = 0.007). Additionally, an increased level of genetically predicted eosinophils was found to be causally associated with a higher risk of diabetic nephropathy (OR = 1.21, 95% CI = 1.08-1.36, pFDR = 0.007, pIVW = 0.001). No evidence of pleiotropy was determined.

Conclusion

Our findings provide evidence of causal associations of circulating WBC count, lymphocyte count and IgA nephropathy, WBC count and renal malignant neoplasms, and eosinophil count and diabetic nephropathy. These results have the potential to contribute to the development of novel diagnostic options and therapeutic strategies for renal disorders.

Comparative transcriptomic analysis revealed dynamic changes of distinct classes of genes during development of the Manila clam (Ruditapes philippinarum)

Background

The Manila clam Ruditapesphilippinarum is one of the most economically important marine shellfish. However, the molecular mechanisms of early development in Manila clams are largely unknown. In this study, we collected samples from 13 stages of early development in Manila clam and compared the mRNA expression pattern between samples by RNA-seq techniques.

Results

We applied RNA-seq technology to 13 embryonic and larval stages of the Manila clam to identify critical genes and pathways involved in their development and biological characteristics. Important genes associated with different morphologies during the early fertilized egg, cell division, cell differentiation, hatching, and metamorphosis stages were identified. We detected the highest number of differentially expressed genes in the comparison of the pediveliger and single pipe juvenile stages, which is a time when biological characteristics greatly change during metamorphosis. Gene Ontology (GO) enrichment analysis showed that expression levels of microtubule protein-related molecules and Rho genes were upregulated and that GO terms such as ribosome, translation, and organelle were enriched in the early development stages of the Manila clam. Kyoto Encyclopedia of Genes and Genomes pathway analysis showed that the foxo, wnt, and transforming growth factor-beta pathways were significantly enriched during early development. These results provide insights into the molecular mechanisms at work during different periods of early development of Manila clams.

Conclusion

These transcriptomic data provide clues to the molecular mechanisms underlying the development of Manila clam larvae. These results will help to improve Manila clam reproduction and development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08813-0.

Identification of key candidate genes for IgA nephropathy using machine learning and statistics based bioinformatics models

Immunoglobulin-A-nephropathy (IgAN) is a kidney disease caused by the accumulation of IgAN deposits in the kidneys, which causes inflammation and damage to the kidney tissues. Various bioinformatics analysis-based approaches are widely used to predict novel candidate genes and pathways associated with IgAN. However, there is still some scope to clearly explore the molecular mechanisms and causes of IgAN development and progression. Therefore, the present study aimed to identify key candidate genes for IgAN using machine learning (ML) and statistics-based bioinformatics models. First, differentially expressed genes (DEGs) were identified using limma, and then enrichment analysis was performed on DEGs using DAVID. Protein-protein interaction (PPI) was constructed using STRING and Cytoscape was used to determine hub genes based on connectivity and hub modules based on MCODE scores and their associated genes from DEGs. Furthermore, ML-based algorithms, namely support vector machine (SVM), least absolute shrinkage and selection operator (LASSO), and partial least square discriminant analysis (PLS-DA) were applied to identify the discriminative genes of IgAN from DEGs. Finally, the key candidate genes (FOS, JUN, EGR1, FOSB, and DUSP1) were identified as overlapping genes among the selected hub genes, hub module genes, and discriminative genes from SVM, LASSO, and PLS-DA, respectively which can be used for the diagnosis and treatment of IgAN.

Identification of Hub Genes and Therapeutic Agents for IgA Nephropathy Through Bioinformatics Analysis and Experimental Validation

Background

IgA nephropathy (IgAN) is the most common primary glomerular disease and the leading cause of the end-stage renal disease in the world. The pathogenesis of IgAN has not been well elucidated, and yet treatment is limited. High-throughput microarray has been applied for elucidating molecular biomarkers and potential mechanisms involved in IgAN. This study aimed to identify the potential key genes and therapeutics associated with IgAN using integrative bioinformatics and transcriptome-based computational drug repurposing approach.

Methods

Three datasets of mRNA expression profile were obtained from the gene expression omnibus database and differentially expressed genes (DEGs) between IgAN glomeruli and normal tissue were identified by integrated analysis. Gene ontology and pathway enrichment analyses of the DEGs were performed by R software, and protein-protein interaction networks were constructed using the STRING online search tool. External dataset and immunohistochemical assessment of kidney biopsy specimens were used for hub gene validation. Potential compounds for IgAN therapy were obtained by Connectivity Map (CMap) analysis and preliminarily verified in vitro. Stimulated human mesangial cells were collected for cell proliferation and cell cycle analysis using cell counting kit 8 and flow cytometry, respectively.

Results

134 DEGs genes were differentially expressed across kidney transcriptomic data from IgAN patients and healthy living donors. Enrichment analysis showed that the glomerular compartments underwent a wide range of interesting pathological changes during kidney injury, focused on anion transmembrane transporter activity and protein digestion and absorption mostly. Hub genes (ITGB2, FCER1G, CSF1R) were identified and verified to be significantly upregulated in IgAN patients, and associated with severity of renal lesions. Computational drug repurposing with the CMap identified tetrandrine as a candidate treatment to reverse IgAN hub gene expression. Tetrandrine administration significantly reversed mesangial cell proliferation and cell cycle transition.

Conclusion

The identification of DEGs and related therapeutic strategies of IgAN through this integrated bioinformatics analysis provides a valuable resource of therapeutic targets and agents of IgAN. Especially, our findings suggest that tetrandrine might be beneficial for IgAN, which deserves future research.

Identification of key biomarkers and signaling pathways and analysis of their association with immune cells in immunoglobulin A nephropathy

Introduction

Immunoglobulin A nephropathy (IgAN) is the most common glomerular disease worldwide, with a poor prognosis. The aim of our study was to identify key biomarkers and their associations with immune cells to aid in the study of IgAN pathology and immunotherapy.

Material and methods

The data of IgAN were downloaded from a public database. The metaMA package and limma package were used to identify differentially expressed mRNAs (DEmRNAs) and differentially expressed miRNAs (DEmiRNAs), respectively. Biological functions of the DEmRNAs were analyzed. Machine learning was used to screen the mRNA biomarkers of IgAN. Pearson's correlation coefficient was used to analyze the correlation between mRNA biomarkers, immune cells and signaling pathways. Moreover, we constructed a miRNAs-mRNAs targeted regulatory network. Finally, we performed in vitro validation of the identified miRNAs and mRNAs.

Results

1205 DEmRNAs and 125 DEmiRNAs were identified. In gene set enrichment analysis (GSEA), tumor necrosis factor α (TNF-α) signaling via nuclear factor κB (NF-κB), apoptosis and MTORC-1 signaling were inhibited in IgAN. 8 mRNA biomarkers were screened by machine learning. In addition, the distribution of 8 immune cell types was found to be significantly different between normal controls and IgAN by difference analysis. Pearson correlation coefficient analysis demonstrated that AKAP8L was significantly negatively correlated with CD4⁺ memory T-cells. AKAP8L was also significantly negatively correlated with TNF-α signaling via NF-κB, apoptosis, and MTORC-1 signaling. Subsequently, 5 mRNA biomarkers predicted corresponding negative regulatory miRNAs.

Conclusions

The identification of 8 important biomarkers and their correlation with immune cells and biological signaling pathways provides new ideas for further study of IgAN.

Comprehensive analysis of IgA nephropathy expression profiles: identification of potential biomarkers and therapeutic agents

Background

IgA nephropathy (IgAN) is a kidney disease recognized by the presence of IgA antibody depositions in kidneys. The underlying mechanisms of this complicated disease are remained to be explored and still, there is an urgent need for the discovery of noninvasive biomarkers for its diagnosis. In this investigation, an integrative approach was applied to mRNA and miRNA expression profiles in PBMCs to discover a gene signature and novel potential targets/biomarkers in IgAN.

Methods

Datasets were selected from gene expression omnibus database. After quality control checking, two datasets were analyzed by Limma to identify differentially expressed genes/miRNAs (DEGs and DEmiRs). Following identification of DEmiR-target genes and data integration, intersecting mRNAs were subjected to different bioinformatic analyses. The intersecting mRNAs, DEmiRs, related transcription factors (from TRRUST database), and long-non coding RNAs (from LncTarD database) were used for the construction of a multilayer regulatory network via Cytoscape.

Result

“GSE25590” (miRNA) and “GSE73953” (mRNA) datasets were analyzed and after integration, 628 intersecting mRNAs were identified. The mRNAs were mainly associated with “Innate immune system”, “Apoptosis”, as well as “NGF signaling” pathways. A multilayer regulatory network was constructed and several hub-DEGs (Tp53, STAT3, Jun, etc.), DEmiRs (miR-124, let-7b, etc.), TFs (NF-kB, etc.), and lncRNAs (HOTAIR, etc.) were introduced as potential factors in the pathogenesis of IgAN.

Conclusion

Integration of two different expression datasets and construction of a multilayer regulatory network not only provided a deeper insight into the pathogenesis of IgAN, but also introduced several key molecules as potential therapeutic target/non-invasive biomarkers.