Expression profile analysis of long non-coding RNAs in condyloma acuminatum in a Chinese female population by RNA-sequencing and bioinformatics [RNA-seq]

This study aimed to explore the differential lncRNAs in CA and their probable function on CA development.We used RNA-sequencing to assess the genome-wide expression levels of lncRNAs in CA and paired adjacent normal tissues. We further validated the candidate lncRNAs in larger-size of CA specimen using RT-qPCR. Furthermore, the functional enrichment analysis for these candidate lncRNAs and differential mRNAs were analyzed using GO terms, KEGG and STRING database. The coexpressed mRNAs for the candidate lncRNAs, calculated by Pearson?s correlation coefficient, were also analyzed using GO analysis.A total of 657 lncRNAs and 1486 mRNAs were found to dysregulated in CA compared to adjacent mucosal tissues. The microarray data of candidate lncRNAs, including HOXC13-AS, HOTAIR, PICSAR, FGF14-AS1, ADGRD1-AS1, TMEM108-AS1, TUSC7 and FGF10-AS1, were validated with a larger-size of specimen using RT-qPCR. The functional GO term and pathways analysis of DEGs are associated with the pathogenesis and development of CA, including cell proliferation and development, cellular adhesion, immune defense, cell migration and chemotaxis, angiogenesis, and maintaining skin hemostasis.LncRNAs that were differentially expressed between CA and normal tissues may be helpful to explore effective recurrence risk markers and potential intervention targets for CA.