Preliminary exploration of the co-regulation of Alzheimer’s disease pathogenic genes by microRNAs and transcription factors

Background: Alzheimer’s disease (AD) is the most common form of age-related neurodegenerative disease. Unfortunately, due to the complexity of pathological types and clinical heterogeneity of AD, there is a lack of satisfactory treatment for AD. Previous studies have shown that microRNAs and transcription factors can modulate genes associated with AD, but the underlying pathophysiology remains unclear.

Methods: The datasets GSE1297 and GSE5281 were downloaded from the gene expression omnibus (GEO) database and analyzed to obtain the differentially expressed genes (DEGs) through the “R” language “limma” package. The GSE1297 dataset was analyzed by weighted correlation network analysis (WGCNA), and the key gene modules were selected. Next, gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis for the key gene modules were performed. Then, the protein-protein interaction (PPI) network was constructed and the hub genes were identified using the STRING database and Cytoscape software. Finally, for the GSE150693 dataset, the “R” package “survivation” was used to integrate the data of survival time, AD transformation status and 35 characteristics, and the key microRNAs (miRNAs) were selected by Cox method. We also performed regression analysis using least absolute shrinkage and selection operator (Lasso)-Cox to construct and validate prognostic features associated with the four key genes using different databases. We also tried to find drugs targeting key genes through DrugBank database.

Results: GO and KEGG enrichment analysis showed that DEGs were mainly enriched in pathways regulating chemical synaptic transmission, glutamatergic synapses and Huntington’s disease. In addition, 10 hub genes were selected from the PPI network by using the algorithm Between Centrality. Then, four core genes (TBP, CDK7, GRM5, and GRIA1) were selected by correlation with clinical information, and the established model had very good prognosis in different databases. Finally, hsa-miR-425-5p and hsa-miR-186-5p were determined by COX regression, AD transformation status and aberrant miRNAs.

Conclusion: In conclusion, we tried to construct a network in which miRNAs and transcription factors jointly regulate pathogenic genes, and described the process that abnormal miRNAs and abnormal transcription factors TBP and CDK7 jointly regulate the transcription of AD central genes GRM5 and GRIA1. The insights gained from this study offer the potential AD biomarkers, which may be of assistance to the diagnose and therapy of AD.