Identification of DEGs in ACC cells after overexpression of Ptch1
To determine the effects of Ptch1 on ACC, we analyzed the RNA-seq data from the ACC cells with the overexpression of Ptch1. A total of 1160 genes were identified with the threshold of P < 0.01. The top up- and down-regulated genes were indicated by the heatmap and volcano plot (Figure 1). The top ten DEGs were listed in Table 1.
Enrichment analysis of DEGs after overexpression of Ptch1 in ACC cells
To further figure out the potential biological processes in ACC cells with the overexpression of Ptch1, we performed the KEGG and GO analyses (Figure 2). We identified the top ten KEGG signaling pathways, including “Calcium signaling pathway”, “Pathogenic Escherichia coli infection”, “Phospholipase D signaling pathway”, “Oxytocin signaling pathway”, “ECM−receptor interaction”, “Platelet activation”, “Complement and coagulation cascades”, “PPAR signaling pathway”, “Arrhythmogenic right ventricular cardiomyopathy”, and “Adipocytokine signaling pathway”. We identified the top ten BP items of GO enrichment, including “axon development”, “axonogenesis”, “cell junction assembly”, “regulation of nervous system development”, “synapse organization”, “axon guidance”, “neuron projection guidance”, “cell−cell adhesion via plasma−membrane adhesion molecules”, “regulation of synapse organization”, and “regulation of synapse structure or activity”. We identified the top ten CC items of GO enrichment, including “synaptic membrane”, “collagen−containing extracellular matrix”, “membrane raft”, “membrane microdomain”, “neuron to neuron synapse”, “postsynaptic density”, “asymmetric synapse”, “postsynaptic membrane”, “intrinsic component of synaptic membrane”, and “intrinsic component of postsynaptic membrane”. We also identified the top ten MF items of GO enrichment, including “DNA−binding transcription activator activity, RNA polymerase II−specific”, “DNA−binding transcription activator activity”, “peptide receptor activity”, “cell adhesion mediator activity”, “calcium−dependent protein binding”, “cell−cell adhesion mediator activity”, “scaffold protein binding”, “extracellular matrix structural constituent conferring tensile strength”, “ionotropic glutamate receptor activity”, and “platelet−derived growth factor binding”.
PPI network analysis
To determine the relationship among the DEGs, we create the PPI network by using 867 nodes and 2328 edges (Combined score > 0.2 as a cutoff, Cytoscope software). Table 2 showed the top ten genes with the highest scores. The top two significant clusters were indicated in Figure 3. We further analyzed the PPI and DEGs with Reactome map (Figure 4) and identified the top ten functional processes including "Response of EIF2AK1 (HRI) to heme deficiency", "Negative regulation of activity of TFAP2 (AP-2) family transcription factors", "Acyl chain remodelling of PG", "Ligand-receptor interactions", "Dissolution of Fibrin Clot", "Defective CHST3 causes SEDCJD", "Collagen degradation", "RUNX3 regulates RUNX1-mediated transcription", "CHL1 interactions", and "Acyl chain remodelling of PS" (Supplemental Table S1).