Circular RNA expression profiles in human bronchial epithelial cells treated with beryllium sulfate.

Circular RNAs (circRNAs), is a novel type of endogenous non-coding RNAs (ncRNAs) participated in the pathogenesis of many diseases. Beryllium is one of the carcinogenesis elements. However, the mechanism and function of circRNAs in human bronchial epithelial cells (16HBE) induced by beryllium sulfate (BeSO) was rarely reported. Therefore, the high-throughput RNA sequencing analysis was performed to detect the circRNA profiles between control groups and BeSO-induced groups. Furthermore, circRNA-miRNA-mRNA network, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, and PPI network analysis were used for bioinformatics analysis. CircRNA sequencing analysis revealed that 36 circRNAs were up-regulated and 35 circRNAs were down-regulated in the BeSO-exposed groups. The selected circRNAs were verified by real-time fluorescent quantitative PCR (qRT-PCR). Hsa_circ_0004214 and hsa_circ_0003586 were validated to be up-regulated, hsa_circ_0047958, hsa_circ_0001944, and hsa_circ_0008982 were down-regulated. The circRNA-miRNA-mRNA network annotated the key signaling pathway including cellular senescence, TNF signaling pathway, NF-kappa B signaling pathway, HIF-1 signaling pathway, and Hippo signaling pathway. The PPI network indicated the most circRNAs might participate in the BeSO toxicity by acting as a sponge for the miR-663b through JAK-STAT signaling pathway. In summary, our study suggests that circRNAs may play roles in the mechanism of beryllium toxicity.