Identification and analysis of epithelial‐mesenchymal transition‐related key long non‐coding RNAs in hypospadias

Abstract EMT dysfunction is a dominant mechanisms of hypospadias. Thus, identification of EMT‐related lncRNAs based on transcriptome sequencing data of hypospadias might provide novel molecular markers and therapeutic targets for hypospadias. First, the microarray data related to hypospadias were downloaded from Gene Expression Omnibus (GEO). Besides, the differentially expressed lncRNAs and messenger RNAs (mRNAs) related to EMT were screened to construct lncRNA‐mRNA co‐expression interaction pairs. In addition, the microRNA (miRNA) prediction analysis was performed through bioinformatics methods to construct a ceRNA network. Moreover, function prediction and function enrichment and pathway analyses were also performed. Finally, the core EMT‐related lncRNAs were verified based on mRNA expression changes and cell functions. A total of 6 EMT‐related lncRNAs were identified and 123 mRNA‐lncRNA co‐expression interaction pairs were screened in this study. Additionally, a ceRNA regulatory network comprising 17 mRNAs, 4 lncRNAs, and 28 miRNAs was constructed based on the prediction of hypospadias‐related miRNAs. The validation results of the dataset GSE121712 revealed that only BEX1 was positively correlated with the expression of the lncRNA GNAS‐AS1 (r = 0.874, P < 0.01), both of which had high expression. The cell experiment results demonstrated that interfering with the expression of GNAS‐AS1 significantly promoted the proliferation, migration, and EMT of cells. Importantly, it was confirmed that GNAS‐AS1 can serve as a ceRNA and play an important role in the EMT of hypospadias. Hence, it may be considered as a potential target in the treatment of this disease.