Bipolar disorder (BD), also known as psychiatric disorder, affects millions of people all over the world. The aim of this investigation was to screen and verify hub genes involved in BD as well as to explore potential molecular mechanisms. The next generation sequencing (NGS) dataset GSE124326 was downloaded from the Gene Expression Omnibus (GEO) database, which contained 480 samples, including 240 BD and 240 normal controls. Differentially expressed genes (DEGs) were filtered and subjected to gene ontology (GO) and pathway enrichment analyses. A protein–protein interaction (PPI) network and module analysis were used to identify hub genes. The miRNA-hub gene regulatory network and TF-hub gene regulatory network were also constructed. Receiver operating characteristic curve (ROC) analysis was used to validate hub genes. In this work, 957 DEGs, including 477 up regulated and 480 down regulated, were obtained from NGS data. The GO and pathway enrichment analyses of the DEGs showed that the up regulated genes were enriched in the neutrophil degranulation, immune system, transport, cytoplasm and enzyme regulator activity, and the down regulated genes were enriched in extracellular matrix organization, diseases of metabolism, multicellular organismal process, cell periphery and metal ion binding. Finally, through analyzing the PPI network, miRNA-hub gene regulatory network and TF-hub gene regulatory network, we screened hub genes UBB, UBE2D1, TUBA1A, RPL11, RPS24, NOTCH3, CAV1, CNBD2, CCNA1 and MYH11 by the Cytoscape software. The current investigation illustrates a characteristic NGS data in BD, which might contribute to the interpretation of the progression of BD and provide novel biomarkers and therapeutic targets for BD.