Screening of potential biomarkers in uterine leiomyomas disease via gene expression profiling analysis

The present study aimed to screen potential biomarkers for uterine leiomyomas disease, particularly target genes associated with the mediator of RNA polymerase II transcription subunit 12 (MED12) mutation. The microarray data of GSE30673, including 10 MED12 wild-type myometrium, 8 MED12 mutation leiomyoma and 2 MED12 wild-type leiomyoma samples, were downloaded from the Gene Expression Omnibus database. Compared with myometrium samples, differently-expressed genes (DEGs) in the MED12 mutation and wild-type leiomyoma samples were identified using the Limma package. The two sets of DEGs obtained were intersected to screen common DEGs. The DEGs in the MED12 mutation and wild-type leiomyoma samples, and common DEGs were defined as group A, B and C. Gene Ontology (GO) and pathway enrichment analyses were performed using the Database for Annotation, Visualization and Integrated Discovery online tool. Based on the Kyoto Encyclopedia of Genes and Genomes database, pathway relation networks were constructed. DEGs in GO terms and pathways were intersected to screen important DEGs. Subsequently, a gene co-expression network was constructed and visualized using Cytoscape software. Reverse transcription-quantitative polymerase chain reaction was used to detect the expression levels of important DEGs. A total of 1,258 DEGs in group A were screened, and enriched for extracellular matrix (ECM) organization and ECM-receptor interaction. In addition, a total of 1,571 DEGs in group B were enriched for cell adhesion. Furthermore, 391 DEGs were involved in extracellular matrix organization. Pathway relation networks of group A, B and C were constructed with nodes of 48, 39, and 28, respectively. Finally, 135 important DEGs were obtained, including Acyl-CoA synthetase medium-chain family member 3, protein S (α) (PROS1) and F11 receptor. A gene co-expression network with 68 nodes was constructed. The expression of caspase 1 (CASP1) and aldehyde dehydrogenase 1 family member A1 (ALDH1A1) was significant higher in SK-UT-1 compared with that in PHM1-31 cells, while the expression of PROS1 was significant lower in SK-UT-1 cells. These results that CASP1, ALDH1A1 and PROS1 may be potential biomarkers for uterine leiomyomas. Furthermore, hematopoietic prostaglandin D synthase and carbonyl reductase 3 (CBR3) may be particular genes associated with the MED12 mutation in this disease.