Alzheimer's disease rewires gene coexpression networks coupling different brain regions.

Connectome studies have shown how Alzheimer's disease (AD) disrupts functional and structural connectivity among brain regions. But the molecular basis of such disruptions is less studied, with most genomic/transcriptomic studies performing within-brain-region analyses. To inspect how AD rewires the correlation structure among genes in different brain regions, we performed an Inter-brain-region Differential Correlation (Inter-DC) analysis of RNA-seq data from Mount Sinai Brain Bank on four brain regions (frontal pole, superior temporal gyrus, parahippocampal gyrus and inferior frontal gyrus, comprising 264 AD and 372 control human post-mortem samples). An Inter-DC network was assembled from all pairs of genes across two brain regions that gained (or lost) correlation strength in the AD group relative to controls at FDR 1%. The differentially correlated (DC) genes in this network complemented known differentially expressed genes in AD, and likely reflects cell-intrinsic changes since we adjusted for cell compositional effects. Each brain region used a distinctive set of DC genes when coupling with other regions, with parahippocampal gyrus showing the most rewiring, consistent with its known vulnerability to AD. The Inter-DC network revealed master dysregulation hubs in AD (at genes ZKSCAN1, SLC5A3, RCC1, IL17RB, PLK4, etc.), inter-region gene modules enriched for known AD pathways (synaptic signaling, endocytosis, etc.), and candidate signaling molecules that could mediate region-region communication. The Inter-DC network generated in this study is a valuable resource of gene pairs, pathways and signaling molecules whose inter-brain-region functional coupling is disrupted in AD, thereby offering a new perspective of AD etiology.

Analysis of drought and heat stress response genes in rice using co-expression network and differentially expressed gene analyses.

Studies on Oryza sativa (rice) are crucial for improving agricultural productivity and ensuring global sustenance security, especially considering the increasing drought and heat stress caused by extreme climate change. Currently, the genes and mechanisms underlying drought and heat resistance in rice are not fully understood, and the scope for enhancing the development of new strains remains considerable. To accurately identify the key genes related to drought and heat stress responses in rice, multiple datasets from the Gene Expression Omnibus (GEO) database were integrated in this study. A co-expression network was constructed using a Weighted Correlation Network Analysis (WGCNA) algorithm. We further distinguished the core network and intersected it with differentially expressed genes and multiple expression datasets for screening. Differences in gene expression levels were verified using quantitative real-time polymerase chain reaction (PCR). OsDjC53, MBF1C, BAG6, HSP23.2, and HSP21.9 were found to be associated with the heat stress response, and it is also possible that UGT83A1 and OsCPn60a1, although not directly related, are affected by drought stress. This study offers significant insights into the molecular mechanisms underlying stress responses in rice, which could promote the development of stress-tolerant rice breeds.

Decoding the principle of cell-fate determination for its reverse control.

Understanding and manipulating cell fate determination is pivotal in biology. Cell fate is determined by intricate and nonlinear interactions among molecules, making mathematical model-based quantitative analysis indispensable for its elucidation. Nevertheless, obtaining the essential dynamic experimental data for model development has been a significant obstacle. However, recent advancements in large-scale omics data technology are providing the necessary foundation for developing such models. Based on accumulated experimental evidence, we can postulate that cell fate is governed by a limited number of core regulatory circuits. Following this concept, we present a conceptual control framework that leverages single-cell RNA-seq data for dynamic molecular regulatory network modeling, aiming to identify and manipulate core regulatory circuits and their master regulators to drive desired cellular state transitions. We illustrate the proposed framework by applying it to the reversion of lung cancer cell states, although it is more broadly applicable to understanding and controlling a wide range of cell-fate determination processes.

Identification and verification of a novel signature that combines cuproptosis-related genes with ferroptosis-related genes in osteoarthritis using bioinformatics analysis and experimental validation.

BACKGROUND: Exploring the pathogenesis of osteoarthritis (OA) is important for its prevention, diagnosis, and treatment. Therefore, we aimed to construct novel signature genes (c-FRGs) combining cuproptosis-related genes (CRGs) with ferroptosis-related genes (FRGs) to explore the pathogenesis of OA and aid in its treatment. MATERIALS AND METHODS: Differentially expressed c-FRGs (c-FDEGs) were obtained using R software. Enrichment analysis was performed and a protein-protein interaction (PPI) network was constructed based on these c-FDEGs. Then, seven hub genes were screened. Three machine learning methods and verification experiments were used to identify four signature biomarkers from c-FDEGs, after which gene set enrichment analysis, gene set variation analysis, single-sample gene set enrichment analysis, immune function analysis, drug prediction, and ceRNA network analysis were performed based on these signature biomarkers. Subsequently, a disease model of OA was constructed using these biomarkers and validated on the GSE82107 dataset. Finally, we analyzed the distribution of the expression of these c-FDEGs in various cell populations. RESULTS: A total of 63 FRGs were found to be closely associated with 11 CRGs, and 40 c-FDEGs were identified. Bioenrichment analysis showed that they were mainly associated with inflammation, external cellular stimulation, and autophagy. CDKN1A, FZD7, GABARAPL2, and SLC39A14 were identified as OA signature biomarkers, and their corresponding miRNAs and lncRNAs were predicted. Finally, scRNA-seq data analysis showed that the differentially expressed c-FRGs had significantly different expression distributions across the cell populations. CONCLUSION: Four genes, namely CDKN1A, FZD7, GABARAPL2, and SLC39A14, are excellent biomarkers and prospective therapeutic targets for OA.

CircRNA-associated ceRNA regulatory networks as emerging mechanisms governing the development and biophysiopathology of epilepsy.

The etiology of epilepsy is ascribed to the synchronized aberrant neuronal activity within the brain. Circular RNAs (circRNAs), a class of non-coding RNAs characterized by their circular structures and covalent linkage, exert a substantial influence on this phenomenon. CircRNAs possess stereotyped replication, transience, repetitiveness, and paroxysm. Additionally, MicroRNA (miRNA) plays a crucial role in the regulation of diverse pathological processes, including epilepsy. CircRNA is of particular significance due to its ability to function as a competing endogenous RNA, thereby sequestering or inhibiting miRNA activity through binding to target mRNA. Our review primarily concentrates on elucidating the pathological and functional roles, as well as the underlying mechanisms, of circRNA-miRNA-mRNA networks in epilepsy. Additionally, it explores the potential utility of these networks for early detection and therapeutic intervention.

Synergistic epistasis among cancer drivers can rescue early tumors from the accumulation of deleterious passengers.

Epistasis among driver mutations is pervasive and explains relevant features of cancer, such as differential therapy response and convergence towards well-characterized molecular subtypes. Furthermore, a growing body of evidence suggests that tumor development could be hampered by the accumulation of slightly deleterious passenger mutations. In this work, we combined empirical epistasis networks, computer simulations, and mathematical models to explore how synergistic interactions among driver mutations affect cancer progression under the burden of slightly deleterious passengers. We found that epistasis plays a crucial role in tumor development by promoting the transformation of precancerous clones into rapidly growing tumors through a process that is analogous to evolutionary rescue. The triggering of epistasis-driven rescue is strongly dependent on the intensity of epistasis and could be a key rate-limiting step in many tumors, contributing to their unpredictability. As a result, central genes in cancer epistasis networks appear as key intervention targets for cancer therapy.

Automatic design of gene regulatory mechanisms for spatial pattern formation.

Gene regulatory mechanisms (GRMs) control the formation of spatial and temporal expression patterns that can serve as regulatory signals for the development of complex shapes. Synthetic developmental biology aims to engineer such genetic circuits for understanding and producing desired multicellular spatial patterns. However, designing synthetic GRMs for complex, multi-dimensional spatial patterns is a current challenge due to the nonlinear interactions and feedback loops in genetic circuits. Here we present a methodology to automatically design GRMs that can produce any given two-dimensional spatial pattern. The proposed approach uses two orthogonal morphogen gradients acting as positional information signals in a multicellular tissue area or culture, which constitutes a continuous field of engineered cells implementing the same designed GRM. To efficiently design both the circuit network and the interaction mechanisms-including the number of genes necessary for the formation of the target spatial pattern-we developed an automated algorithm based on high-performance evolutionary computation. The tolerance of the algorithm can be configured to design GRMs that are either simple to produce approximate patterns or complex to produce precise patterns. We demonstrate the approach by automatically designing GRMs that can produce a diverse set of synthetic spatial expression patterns by interpreting just two orthogonal morphogen gradients. The proposed framework offers a versatile approach to systematically design and discover complex genetic circuits producing spatial patterns.

Data-driven modeling of core gene regulatory network underlying leukemogenesis in IDH mutant AML.

Acute myeloid leukemia (AML) is characterized by uncontrolled proliferation of poorly differentiated myeloid cells, with a heterogenous mutational landscape. Mutations in IDH1 and IDH2 are found in 20% of the AML cases. Although much effort has been made to identify genes associated with leukemogenesis, the regulatory mechanism of AML state transition is still not fully understood. To alleviate this issue, here we develop a new computational approach that integrates genomic data from diverse sources, including gene expression and ATAC-seq datasets, curated gene regulatory interaction databases, and mathematical modeling to establish models of context-specific core gene regulatory networks (GRNs) for a mechanistic understanding of tumorigenesis of AML with IDH mutations. The approach adopts a new optimization procedure to identify the top network according to its accuracy in capturing gene expression states and its flexibility to allow sufficient control of state transitions. From GRN modeling, we identify key regulators associated with the function of IDH mutations, such as DNA methyltransferase DNMT1, and network destabilizers, such as E2F1. The constructed core regulatory network and outcomes of in-silico network perturbations are supported by survival data from AML patients. We expect that the combined bioinformatics and systems-biology modeling approach will be generally applicable to elucidate the gene regulation of disease progression.

Trans-eQTL mapping in gene sets identifies network effects of genetic variants.

Nearly all trait-associated variants identified in genome-wide association studies (GWASs) are noncoding. The cis regulatory effects of these variants have been extensively characterized, but how they affect gene regulation in trans has been the subject of fewer studies because of the difficulty in detecting trans-expression quantitative loci (eQTLs). We developed trans-PCO for detecting trans effects of genetic variants on gene networks. Our simulations demonstrate that trans-PCO substantially outperforms existing trans-eQTL mapping methods. We applied trans-PCO to two gene expression datasets from whole blood, DGN (N = 913) and eQTLGen (N = 31,684), and identified 14,985 high-quality trans-eSNP-module pairs associated with 197 co-expression gene modules and biological processes. We performed colocalization analyses between GWAS loci of 46 complex traits and the trans-eQTLs. We demonstrated that the identified trans effects can help us understand how trait-associated variants affect gene regulatory networks and biological pathways.

Plasticity-led and mutation-led evolutions are different modes of the same developmental gene regulatory network.

The standard theory of evolution proposes that mutations cause heritable variations, which are naturally selected, leading to evolution. However, this mutation-led evolution (MLE) is being questioned by an alternative theory called plasticity-led evolution (PLE). PLE suggests that an environmental change induces adaptive phenotypes, which are later genetically accommodated. According to PLE, developmental systems should be able to respond to environmental changes adaptively. However, developmental systems are known to be robust against environmental and mutational perturbations. Thus, we expect a transition from a robust state to a plastic one. To test this hypothesis, we constructed a gene regulatory network (GRN) model that integrates developmental processes, hierarchical regulation, and environmental cues. We then simulated its evolution over different magnitudes of environmental changes. Our findings indicate that this GRN model exhibits PLE under large environmental changes and MLE under small environmental changes. Furthermore, we observed that the GRN model is susceptible to environmental or genetic fluctuations under large environmental changes but is robust under small environmental changes. This indicates a breakdown of robustness due to large environmental changes. Before the breakdown of robustness, the distribution of phenotypes is biased and aligned to the environmental changes, which would facilitate rapid adaptation should a large environmental change occur. These observations suggest that the evolutionary transition from mutation-led to plasticity-led evolution is due to a developmental transition from robust to susceptible regimes over increasing magnitudes of environmental change. Thus, the GRN model can reconcile these conflicting theories of evolution.

CoVar: A generalizable machine learning approach to identify the coordinated regulators driving variational gene expression.

Network inference is used to model transcriptional, signaling, and metabolic interactions among genes, proteins, and metabolites that identify biological pathways influencing disease pathogenesis. Advances in machine learning (ML)-based inference models exhibit the predictive capabilities of capturing latent patterns in genomic data. Such models are emerging as an alternative to the statistical models identifying causative factors driving complex diseases. We present CoVar, an ML-based framework that builds upon the properties of existing inference models, to find the central genes driving perturbed gene expression across biological states. Unlike differentially expressed genes (DEGs) that capture changes in individual gene expression across conditions, CoVar focuses on identifying variational genes that undergo changes in their expression network interaction profiles, providing insights into changes in the regulatory dynamics, such as in disease pathogenesis. Subsequently, it finds core genes from among the nearest neighbors of these variational genes, which are central to the variational activity and influence the coordinated regulatory processes underlying the observed changes in gene expression. Through the analysis of simulated as well as yeast expression data perturbed by the deletion of the mitochondrial genome, we show that CoVar captures the intrinsic variationality and modularity in the expression data, identifying key driver genes not found through existing differential analysis methodologies.

An integrated approach to understand the regulatory role of miR-27 family in breast cancer metastasis.

One of the prime reasons of increasing breast cancer mortality is metastasizing cancer cells. Owing to the side effects of clinically available drugs to treat breast cancer metastasis, it is of utmost importance to understand the underlying biogenesis of breast cancer tumorigenesis. In-silico identification of potential RNAs might help in utilizing the miR-27 family as a therapeutic target in breast cancer. The experimentally verified common interacting mRNAs for miR27 family are retrieved from three publicly available databases- TargetScan, miRDB and miRTarBase. Finally on comparing the common genes with HCMDB and GEPIA data, four breast cancer-associated differentially expressed metastatic mRNAs (GATA3, ENAH, ITGA2 and SEMA4D) are obtained. Corresponding to the miR27 family and associated mRNAs, interacting drugs are retrieved from Sm2mir and CTDbase, respectively. The interaction network-based approach was utilized to obtain the hub RNAs and triad modules by employing the 'Cytohubba' and 'MClique' plugins, respectively in Cytoscape. Further, sample-, subclass- and promoter methylation-based expression analyses reveals GATA3 and ENAH to be the most significant mRNAs in breast cancer metastasis having >10% genetic alteration in both METABRIC Vs TCGA datasets as per their oncoprint analysis via cBioPortal. Additionally, survival analysis in Oncolnc reveals SEMA4D as survival biomarker. Interactions among the miR27 family, their target mRNAs and drugs interacting with miRNAs and mRNAs can be extensively explored in both in-vivo and in-vitro setups to assess their therapeutic potential in the diminution of breast cancer.

A domain swapping strategy to create modular transcriptional regulators for novel topology in genetic network.

Transcriptional regulators generate connections between biological signals and genetic outputs. They are used robustly for sensing input signals in building genetic circuits. However, each regulator can only generate a fixed connection, which generates constraints in linking multiple signals for more complex processes. Recent studies discovered that a domain swapping strategy can be applied to various regulator families to create modular regulators for new signal-output connections, significantly broadening possibilities in circuit design. Here we review the development of this emerging strategy, the use of resulting modular regulators for creating novel genetic response behaviors, and current limitations and solutions for further advancing the design of modular regulators.

Gene networks and systems biology in Alzheimer's disease: Insights from multi-omics approaches.

Despite numerous studies in the field of dementia and Alzheimer's disease (AD), a comprehensive understanding of this devastating disease remains elusive. Bulk transcriptomics have provided insights into the underlying genetic factors at a high level. Subsequent technological advancements have focused on single-cell omics, encompassing techniques such as single-cell RNA sequencing and epigenomics, enabling the capture of RNA transcripts and chromatin states at a single cell or nucleus resolution. Furthermore, the emergence of spatial omics has allowed the study of gene responses in the vicinity of amyloid beta plaques or across various brain regions. With the vast amount of data generated, utilizing gene regulatory networks to comprehensively study this disease has become essential. This review delves into some techniques employed in the field of AD, explores the discoveries made using these techniques, and provides insights into the future of the field.

Analysis of ceRNA Network and Identification of Potential Treatment Target and Biomarkers of Endothelial Cell Injury in Sepsis.

Background: Sepsis is a complex clinical syndrome caused by a dysregulated host immune response to infection. This study aimed to identify a competing endogenous RNA (ceRNA) network that can greatly contribute to understanding the pathophysiological process of sepsis and determining sepsis biomarkers. Methods: The GSE100159, GSE65682, GSE167363, and GSE94717 datasets were obtained from the Gene Expression Omnibus (GEO) database. Weighted gene coexpression network analysis was performed to find modules possibly involved in sepsis. A long noncoding RNA-microRNA-messenger RNA (lncRNA-miRNA-mRNA) network was constructed based on the findings. Single-cell analysis was performed. Human umbilical vein endothelial cells were treated with lipopolysaccharide (LPS) to create an in vitro model of sepsis for network verification. Reverse transcription-polymerase chain reaction, fluorescence in situ hybridization, and luciferase reporter genes were used to verify the bioinformatic analysis. Result: By integrating data from three GEO datasets, we successfully constructed a ceRNA network containing 18 lncRNAs, 7 miRNAs, and 94 mRNAs based on the ceRNA hypothesis. The lncRNA ZFAS1 was found to be highly expressed in LPS-stimulated endothelial cells and may thus play a role in endothelial cell injury. Univariate and multivariate Cox analyses showed that only SLC26A6 was an independent predictor of prognosis in sepsis. Overall, our findings indicated that the ZFAS1/hsa-miR-449c-5p/SLC26A6 ceRNA regulatory axis may play a role in the progression of sepsis. Conclusion: The sepsis ceRNA network, especially the ZFAS1/hsa-miR-449c-5p/SLC26A6 regulatory axis, is expected to reveal potential biomarkers and therapeutic targets for sepsis management.

Single-cell biclustering for cell-specific transcriptomic perturbation detection in AD progression.

The pathogenesis of Alzheimer disease (AD) involves complex gene regulatory changes across different cell types. To help decipher this complexity, we introduce single-cell Bayesian biclustering (scBC), a framework for identifying cell-specific gene network biomarkers in scRNA and snRNA-seq data. Through biclustering, scBC enables the analysis of perturbations in functional gene modules at the single-cell level. Applying the scBC framework to AD snRNA-seq data reveals the perturbations within gene modules across distinct cell groups and sheds light on gene-cell correlations during AD progression. Notably, our method helps to overcome common challenges in single-cell data analysis, including batch effects and dropout events. Incorporating prior knowledge further enables the framework to yield more biologically interpretable results. Comparative analyses on simulated and real-world datasets demonstrate the precision and robustness of our approach compared to other state-of-the-art biclustering methods. scBC holds potential for unraveling the mechanisms underlying polygenic diseases characterized by intricate gene coexpression patterns.

Identification of microRNA-mRNA Regulatory Networks with Therapeutic Values in Alzheimer's Disease by Bioinformatics Analysis.

Background: Alzheimer's disease (AD) is the most prevalent neurological disorder worldwide, affecting approximately 24 million individuals. Despite more than a century of research on AD, its pathophysiology is still not fully understood. Objective: Recently, genetic studies of AD have focused on analyzing the general expression profile by employing high-throughput genomic techniques such as microarrays. Current research has leveraged bioinformatics advancements in genetic science to build upon previous efforts. Methods: Data from the GSE118553 dataset used in this investigation, and the analyses carried out using programs such as Limma and BioBase. Differentially expressed genes (DEGs) and differentially expressed microRNAs (DEmiRs) associated with AD identified in the studied areas of the brain. Target genes of the DEmiRs identified using the MultiMiR package. Gene ontology (GO) completed using the Enrichr website, and the protein-protein interaction (PPI) network for these genes drawn using STRING and Cytoscape software. Results: The findings introduced DEGs including CTNNB1, PAK2, MAP2K1, PNPLA6, IGF1R, FOXL2, DKK3, LAMA4, PABPN1, and GDPD5, and DEmiRs linked to AD (miR-106A, miR-1826, miR-1253, miR-10B, miR-18B, miR-101-2, miR-761, miR-199A1, miR-379 and miR-668), (miR-720, miR-218-2, miR-25, miR-602, miR-1226, miR-548K, miR-H1, miR-410, miR-548F2, miR-181A2), (miR-1470, miR-651, miR-544, miR-1826, miR-195, miR-610, miR-599, miR-323, miR-587 and miR-340), and (miR-1282, miR-1914, miR-642, miR-1323, miR-373, miR-323, miR-1322, miR-612, miR-606 and miR-758) in cerebellum, frontal cortex, temporal cortex, and entorhinal cortex, respectively. Conclusions: The majority of the genes and miRNAs identified by our findings may be employed as biomarkers for prediction, diagnosis, or therapy response monitoring.

Dissecting gene expression networks in the developing hippocampus through the lens of NEIL3 depletion.

Gene regulation in the hippocampus is fundamental for its development, synaptic plasticity, memory formation, and adaptability. Comparisons of gene expression among different developmental stages, distinct cell types, and specific experimental conditions have identified differentially expressed genes contributing to the organization and functionality of hippocampal circuits. The NEIL3 DNA glycosylase, one of the DNA repair enzymes, plays an important role in hippocampal maturation and neuron functionality by shaping transcription. While differential gene expression (DGE) analysis has identified key genes involved, broader gene expression patterns crucial for high-order hippocampal functions remain uncharted. By utilizing the weighted gene co-expression network analysis (WGCNA), we mapped gene expression networks in immature (p8-neonatal) and mature (3 m-adult) hippocampal circuits in wild-type and NEIL3-deficient mice. Our study unveiled intricate gene network structures underlying hippocampal maturation, delineated modules of co-expressed genes, and pinpointed highly interconnected hub genes specific to the maturity of hippocampal subregions. We investigated variations within distinct gene network modules following NEIL3 depletion, uncovering NEIL3-targeted hub genes that influence module connectivity and specificity. By integrating WGCNA with DGE, we delve deeper into the NEIL3-dependent molecular intricacies of hippocampal maturation. This study provides a comprehensive systems-level analysis for assessing the potential correlation between gene connectivity and functional connectivity within the hippocampal network, thus shaping hippocampal function throughout development.

Quantifying Landscape-Flux via Single-Cell Transcriptomics Uncovers the Underlying Mechanism of Cell Cycle.

Recent developments in single-cell sequencing technology enable the acquisition of entire transcriptome data. Understanding the underlying mechanism and identifying the driving force of transcriptional regulation governing cell function directly from these data remains challenging. This study reconstructs a continuous vector field of the cell cycle based on discrete single-cell RNA velocity to quantify the single-cell global nonequilibrium dynamic landscape-flux. It reveals that large fluctuations disrupt the global landscape and genetic perturbations alter landscape-flux, thus identifying key genes in maintaining cell cycle dynamics and predicting associated functional effects. Additionally, it quantifies the fundamental energy cost of the cell cycle initiation and unveils that sustaining the cell cycle requires curl flux and dissipation to maintain the oscillatory phase coherence. This study enables the inference of the cell cycle gene regulatory networks directly from the single-cell transcriptomic data, including the feedback mechanisms and interaction intensity. This provides a golden opportunity to experimentally verify the landscape-flux theory and also obtain its associated quantifications. It also offers a unique framework for combining the landscape-flux theory and single-cell high-through sequencing experiments for understanding the underlying mechanisms of the cell cycle and can be extended to other nonequilibrium biological processes, such as differentiation development and disease pathogenesis.

Identification of metabolic pathways and key genes associated with atypical parkinsonism using a systems biology approach.

Atypical parkinsonism (AP) is a group of complex neurodegenerative disorders with marked clinical and pathophysiological heterogeneity. The use of systems biology tools may contribute to the characterization of hub-bottleneck genes, and the identification of its biological pathways to broaden the understanding of the bases of these disorders. A systematic search was performed on the DisGeNET database, which integrates data from expert curated repositories, GWAS catalogues, animal models and the scientific literature. The tools STRING 11.0 and Cytoscape 3.8.2 were used for analysis of protein-protein interaction (PPI) network. The PPI network topography analyses were performed using the CytoHubba 0.1 plugin for Cytoscape. The hub and bottleneck genes were inserted into 4 different sets on the InteractiveVenn. Additional functional enrichment analyses were performed to identify Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and gene ontology for a described set of genes. The systematic search in the DisGeNET database identified 485 genes involved with Atypical Parkinsonism. Superimposing these genes, we detected a total of 31 hub-bottleneck genes. Moreover, our functional enrichment analyses demonstrated the involvement of these hub-bottleneck genes in 3 major KEGG pathways. We identified 31 highly interconnected hub-bottleneck genes through a systems biology approach, which may play a key role in the pathogenesis of atypical parkinsonism. The functional enrichment analyses showed that these genes are involved in several biological processes and pathways, such as the glial cell development, glial cell activation and cognition, pathways were related to Alzheimer disease and Parkinson disease. As a hypothesis, we highlight as possible key genes for AP the MAPT (microtubule associated protein tau), APOE (apolipoprotein E), SNCA (synuclein alpha) and APP (amyloid beta precursor protein) genes.