Analyses of circRNA and mRNA Profiles in Vogt–Koyanagi–Harada Disease

Abstract

Recent studies revealed that circular RNAs (circRNAs) are important in numerous biological process and involved in autoimmune diseases. However, their role in Vogt–Koyanagi–Harada (VKH) disease, a classical autoimmune disease, is not yet known. This research aimed to study the expression profile of mRNAs, microRNAs (miRNAs) and circRNAs and investigate the influence of circRNAs on the pathogenesis of VKH disease. We identified circRNAs, miRNAs, and mRNAs expression profiles in CD4+ T cells between 4 VKH patients and 3 healthy controls using the whole-transcriptome sequencing (RNA-seq) technique. We discovered that a total of 5088 mRNAs, 451 circRNAs and 433 miRNAs were differently expressed. The GO and KEGG pathway enrichment analyses were performed for significantly differentially expressed circRNAs and mRNAs. GSEA was conducted for all mRNAs. The functional enrichment suggested that the inflammatory response, the adaptive immune response, NF-kappa B signaling pathway, Th17 cell differentiation, Th1 and Th2 cell differentiation and T cell receptor signaling pathway were associated with VKH disease. In addition, based on the immune-related genes we screened, the circRNA-miRNA-mRNA ceRNA network was analyzed and constructed. Ten differently expressed mRNAs (LAT, ZAP70, ITK, ICOS, RASGRP1, PAG1, PLCG1, PRKCQ, LCK, CARD11) and 5 differently expressed circRNAs (hsa_circ_0033144, hsa_circ_0000233, hsa_circ_0000396, hsa_circ_0001924, hsa_circ_0001320) were selected to be validated by Real-time qPCR (RT-qPCR). The results of RT-qPCR turned out to be consistent with RNA-seq data. Further analysis showed that hsa_circ_0001320 and hsa_circ_0001924 may serve as crucial candidate marker genes of VKH disease. These results reveal that circRNAs may have a crucial immunomodulatory function in the pathophysiological process of VKH disease.