Microarray expression profile of mRNAs and long noncoding RNAs and the potential role of PFK-1 in infantile hemangioma

Infantile hemangioma (IH) is the most common benign tumor in children. Long noncoding RNAs (lncRNAs) play a critical role in tumorigenesis. However, the expression levels and biological functions of lncRNAs in IH have not been well-studied. This study aimed to analyze the expression profile of lncRNAs and mRNAs in proliferating and involuting IHs. The expression profiles of lncRNAs and mRNAs in proliferating and involuting IHs were identified by microarray analysis. Subsequently, detailed bioinformatics analyses were performed. Finally, quantitative real-time polymerase chain reaction (qRT-PCR) and immunohistochemistry (IHC) analyses were conducted to validate the microarray results. In total, 146 differentially expressed (DE) lncRNAs and 374 DE mRNAs were identified. The DE mRNAs were enriched mostly in angiogenesis-related biological processes (BPs) and pathways by bioinformatics analysis. In addition, metabolism-related BPs (e.g., "glycogen biosynthetic process" and "metabolic process") and pathways (e.g., "oxidative phosphorylation") were identified. A lncRNA-mRNA co-expression network was constructed from 42 DE lncRNAs and 217 DE mRNAs. Twelve lncRNAs were predicted to have cis-regulated target genes. The microarray analysis results were validated by qRT-PCR using 5 randomly selected lncRNAs and 13 mRNAs. The IHC results revealed that both LOXL2 and FPK-1 exhibited higher protein expression levels in proliferating IH than in involuting IH. Moreover, inhibition of PFK-1 could suppress hemangioma-derived endothelial cell proliferation and migration, induce cell arrest, and reduce glucose uptake and lactate and ATP production. The findings suggest that the identified DE lncRNAs and mRNAs may be associated with the pathogenesis of IH. The data presented herein can improve our understanding of IH development and provide direction for further studies investigating the mechanism underlying IH.