Expression Prole of Peripheral Immune Cells-derived Coding and Long Non-coding RNAs in Patients With Proliferative Vitreoretinopathy

Introduction: Peripheral immune response has been revealed to play a critical role in proliferative vitreoretinopathy (PVR). However, the reliable immune-related factors that are acting as prognostic indicators or therapeutic targets for PVR remain to explore further. Methods: In the current study, we applied whole-transcriptome sequencing to prole peripheral blood mononuclear cells (PBMCs) from PVR patients and also analyzed lncRNA-mRNA interactions in peripheral immune cells to explore the pathways that might mediate immunopathology and resultant retinal damage in PVR. Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses and Ingenuity Pathway Analysis (IPA) were employed to classify the function of these differentially expressed genes (DEGs). Results: Compared to the controls, there were 319 genes upregulated, and 191 genes downregulated in PVR patients. GO, and KEGG enrichment analyses as well as IPA showed that these upregulated genes were signicantly enriched in immune-related and infection-relate terms. Immune-related gene NFKBIA, CXCL2, and CXCL8 were detected as hub-genes in the co-expression network, while lncRNAs such as AC007032.1, AC037198.2, AL929472.2, and SLED1 were highly co-expressed with them. lncRNA-mRNA interactions analysis also showed that putative targeted genes of these differentially expressed lncRNAs were also signicantly enriched in immune-related or infection-relate pathways. Conclusion: Our study highlights the transformation of immune-related genes/pathways in PVR by comparing controls, and validates several critical genes and lncRNAs, which are serving as potential diagnostic markers for PVR patients. closely related to PVR change. We found immunoglobulin and its receptors as well as antigen binding were all up regulated in PBMC from PVR patient. Further, RNAs related to biological process of adaptive immune response and lymphocytes activation were also up regulated indicating microbial infection might play a role in the disease process. In IPA analysis, we found in enriched canonical pathways, 23 while SLED1 was found up regulated in peripheral blood cells of systemic lupus erythematosus patients. 24 Additionally, within IPA analysis, virus infection relate pathways and immune relate networks were highlighted in PVR patients. Further, the most obvious changed LncRNA AC037198.2 and ZNF433-AS1 were selected to verify and were proved their actual change in PBMC from PVR patients. LncRNA AC037198.2 is associate with THBS1(thrombospondin 1) gene, which encoded a secreted protein to mediate cell-to-cell and cell-to-matrix interactions. As for ZNF433-AS1, this LncRNA can suppress ZNF433, which belongs to transcriptional factors with the zinc nger motif, and was that play an important role in multiple sclerosis, which is an autoimmune disease.


Introduction
Proliferative vitreoretinopathy (PVR) is a critically blinding complication that occurs during rhegmatogenous retinal detachment (RRD), before or after surgery. It is caused by the proliferation of glial cells or RPE cells to form a brous membrane at the neural retinal surface or even in the retina. 1 Immune cells involved in the pathogenesis of PVR, such as monocyte/macrophage in ltration and activity contributed to the progress of PVR. 2 And, some in ammatory-related genes are identi ed that can predict the susceptibility of PVR in populations. 3 Recently, the interaction of long non-coding RNAs (lncRNAs) and protein-coding RNAs (mRNAs) has been revealed playing a signi cant role in diseases related to in ammation and immunity. 4,5 However, despite some lncRNAs, such as MALAT1 that had been unveiled by some publications, 6,7 there are no reliable lncRNAs currently implicating in clinical practice acting as prognostic indicators or therapeutic targets for PVR patients. Thus, subsequent studies to identify more critical lncRNAs associated with PVR are warranted.
In the present study, we hypothesized that transcriptomes of peripheral immune cells in PVR might be altered, which are the potential cause for the initiation or progression of PVR. Peripheral blood samples were taken from patients undergoing standard three-port pars plana vitrectomy for indications of PVR secondary to RRD and the entire transcriptome sequencing was performed. Patients with idiopathic epiretinal membrane (iERM) underwent pars plana vitrectomy were used as controls. 8,9 Patients were excluded from analysis if they were with systemic diseases (e.g., diabetes, immunological diseases, infections, etc.) that could in uence systemic in ammation. With the sequencing data, we applied Gene Ontology (GO) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis and Ingenuity Pathway Analysis (IPA), as well as gene-lncRNA co-expression analysis to see the differences in transcriptome level between PVR patients and controls. We also validated some selected differentially expressed lncRNAs by quantitative real-time polymerase chain reaction (qPCR) assay. The results of the current study provide novel insights into PVR pathogenesis and treatment therapeutic targets.

Methods And Design
Ethics statement and clinical sample collection This study was approved by the ethics committees of the Zhongshan Ophthalmic Center. The peripheral blood samples were taken in accordance with the Declaration of Helsinki and written consent was obtained from all participants. Patients diagnosed as primary rhegmatogenous retinal detachment (RRD) with serious PVR (≥ Grade C), 10 and were scheduled to have pars plana vitrectomy (PPV) from October 2018 to February 2019 were included (PVR group). Patients diagnosed as idiopathic epiretinal membrane (iERM) and scheduled to have PPV during the same period were included as negative controls (iERM group). The blood samples were collected from patients before surgery.

RNA extraction and cDNA library construction
Using Ficoll-Paque™ PREMIUM Media (GE Healthcare Life Sciences, Massachusetts, America) and SepMate™-50 (STEMCELL Technologies, Vancouver, Canada), approximately 5 mL of anticoagulated peripheral blood was centrifuged at 500 g for 30 min, and the thin white layer was collected as PBMC. RNeasy Mini Kit (Qiagen Corporation, Hilden, Germany) was used to extract the total RNA from PBMC. After using Agilent 4200 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) to determine the RNA integrity number and quantity, the rRNA was removed by Epicentre Ribo-zero rRNA Removal Kit (Epicentre, Madison, WI, USA). Subsequently, the DNA was cleaned by DNase and then captured and puri ed by magnet beads (Vazyme, Nanjing, China). The puri ed RNA was interrupted into short fragments by adding fragmentation buffer, then the rst-strand cDNA was synthesized and double-strand cDNA was obtained thereafter with VAHTS Universal V6 RNA-seq Library Prep Kit for Illumina (Vazyme, Nanjing, China). Adapters were then connected to the cDNAs and Agilent 4200 Bioanalyzer, as well as qPCR, were used to verify the fragment size and amplify the templates. The constructed library was then loaded on the Illumina HiSeq X Ten system for sequencing.

Sequence analysis and functional annotation
We used FastQC rst to check the quality of raw reads and then processed with Cutadapt software to generate clean reads by trim adaptor sequences and removed the low-quality sequences. 11 HISAT2 was used to align the clean reads to the hg19 human reference genome. 12 Mapped reads used the featureCount software to obtain gene expression pro les, [13][14][15] and then the differential expression analysis was executed via the DESeq2 package on R programmer. 16 |Fold change| > 1.5 and P-value < .05 were decided upon as the signi cance of differentially expressed genes. GO and KEGG enrichment analyses were then performed using the Clusterpro ler R package. 17 For further functional analysis, differentially expressed transcriptome between PVR and iERM patients, including gene symbols and expression values were uploaded into IPA software (Qiagen, Germany).
The canonical pathways, diseases and biofunctions as well as gene networks analysis were processed.

QPCR validation
Total RNA was isolated as mentioned above, and the cDNA was then synthesized with HiScript II Q Select RT SuperMix for qPCR (Vazyme, Nanjing, China). Roche lightcycler 96 was then used to perform qPCR. Beta-actin was used as the internal control. The melting curve was used to con rm reaction speci city and relative expression was calculated by the 2 −ΔΔCt method. Primer sequences and product length were listed in Table 1.

Characteristics of the subjects
There were six males and six females in the PVR group, and eight males and ve females in the iERM group. The average age of the PVR group was 53.3 ± 10.6 years old, and that of the iERM group was 60.0 ± 9.7 years old. The differences in age and gender between PVR and iERM groups were not statistically different (P > .05).
Identi cation of differentially expressed transcriptome Using 1.5-fold expression difference as a cutoff, 510 genes were found differentially expressed between PVR patients and iERM patients, among which 319 were upregulated and 191 were downregulated (Fig. 1A). As expected, these two kinds of patients could be clustered into separated groups using the differentially expressed genes (DEGs) (Fig. 1B), highlighting the apparent transcriptomic difference between PVR and iERM. Differentially expressed transcripts between PVR and iERM was also analyzed in this study.
We also performed an unsupervised cluster analysis with the DETs. Similar to the gene clustering analysis, samples were signi cantly separated in accordance with two patient groups.

Gene Ontology analysis and KEGG analysis
To explore the biology underlying the deferentially expressed gene further, we performed an overrepresentation analysis of GO terms and KEGG pathways using the R program with the ClusterPro ler package. It is very interesting to note that most of the upregulated genes were enriched in immune system-related terms, including immunoglobulin, complement, and immune response ( Fig. 2A); while the downregulated GO terms were irrelevant to the immune system ( Figure S3A; supplementary le).
Upregulated genes of PVR patients compared with iERM patients were enriched in many pathways from KEGG (Fig. 2B). To our surprise, many infectious pathways, including malaria, legionellosis, and Chagas disease, were involved with high signi cance. Besides, pathways affecting the immune system, including IL-17 signaling, TNF signaling, as well as rheumatoid arthritis, were enriched. Much different from the upregulated genes, the downregulated genes were poorly enriched in KEGG pathway enrichment ( Figure S3B; supplementary le).

Enriched Pathway, diseases, biofunctions and interaction network by IPA analysis
We used IPA to analysis the signi cant changed 752 transcriptomes (P < 0.05 and fold-change > 1.5) between patient PVR and iERM patients. With IPA, we revealed these changed genes are closely related to 26 canonical pathways. The most signi cant of those are shown in Table 2. 5 out of top ten pathways are related to innate and adaptive immune cells especially Th17 cells and its cytokines IL-17A and IL-17F. With IPA, we also found similar change as we revealed in GO and KEGG analysis: Pathways related to immune reaction were also enriched. For the top four scored networks (Table S1, supplementary le), 2 of them were immune function related. These two networks are shown in Fig. 3: Cell-To-Cell Signaling and Interaction, Cellular Movement, Immune Cell Tra cking (score 23, Fig. 3A); and Gastrointestinal Disease, In ammatory Disease, In ammatory Response (score 21, Fig. 3B). mRNA and lncRNA co-expression analysis As LncRNAs are also important participants in disease process. Based on the correlation of the signi cantly regulated protein-coding mRNA with lncRNA, we constructed co-expression networks to analyze their interaction and to nd out the potential therapeutic target. Pearson correlation coe cient was calculated for each pair and signi cantly correlated RNA pairs were chosen (r > 0.8, P < .05 as the threshold). As shown in Fig. 4, immune-related gene NFKBIA, and chemokines CXCL2 and CXCL8 were of high hubness in the co-expression network, highlighting their important biological role in the difference between PVR and iERM. Besides, we found several lncRNAs were highly co-expressed with these genes, including AC007032.1, AC037198.2, AL929472.2, SLED1, etc., indicating the key regulation values of them.
We searched further to look for the lncRNA-mRNA interactions, inspecting the chromosome location of regulated lncRNAs to see if they had up-or down-stream 10 kb or overlapping mRNAs. As a result, 630 lncRNAs that had nearby DEGs were found, including 290 upregulated and 340 downregulated lncRNAs. GO, and KEGG analyses and also IPA were conducted as above ( Fig. 5 and Table 3). Like that of DETs GO analysis, the upregulated lncRNAs had more enriched GO terms compared to the downregulated lncRNAs though few of which was immune related (Fig. 5A and Figure S4A). For the KEGG analysis, fewer pathways were enriched( Figure 5B and Figure S4B). results above might due to many functions of lncRNA that are still unclear. However, in IPA analysis, though none of the top ten enriched pathways were immune related pathways, we found Virus Entry via Endocytic Pathways were enriched. Further, the top 2 scored networks are relate to Ophthalmic Disease or Immunological Disease respectively (Table   S2; supplementary le).

Validation of differentially expressed lncRNAs
We used qPCR assay to validate the deferentially expressed lncRNAs, while AC037198.2, and ZNF433-AS1 were selected based on differential expression and co-expression analyses. As shown in Fig. 6, these lncRNAs were signi cantly different among PVR and iERM patients.

Discussion
PVR is still the leading cause of vitreoretinal surgery failure, mainly through retina re-detachment and even intraretinal brosis. The incidence of PVR is estimated to be 5-10%. 18 However, as the mechanisms of PVR are still not very clear, it is di cult to predict or treat the condition despite the many efforts that are still being made. In the present study, our results identi ed differentially expressed mRNA and lncRNA in PBMCs of PVR patients and revealed that gene expression pro le and molecular signature of PVR patients.
The proliferation of cells, mainly RPE and glial cells, is the essential point of PVR development. Clinicians have tried for more than four decades to inhibit the proliferation of cells in the vitreous to prevent PVR but have not had impressive progress. 18 This situation has raised the question of whether there is anything abnormal out of the eye in this disease. The answer to this might come from the immune system.
In the review by Pastor et al., the authors proposed that ischemia, blood-retina barrier breakdown, and in ammation lead to the nal PVR based on the collaborative genetic study named "Retina 4 Project", which found 30 in ammatory-related genes were responsible for PVR. 4 Further, a single nucleotide polymorphism (SNP) analysis in peripheral blood from PVR patients shown that TNF locus which encompasses the gene of TNFα contributes to the development of PVR. 19 It is suggested that PVR might not only be a "local in ammatory condition'', but also could be affected by system regulation.
Compare with these studies, with RNA from peripheral blood mononuclear cells and with GO, KEGG and IPA analysis, we revealed that the immune related pathways or components are closely related to PVR change. We found immunoglobulin and its receptors as well as antigen binding were all up regulated in PBMC from PVR patient. Further, RNAs related to biological process of adaptive immune response and lymphocytes activation were also up regulated indicating microbial infection might play a role in the disease process. In IPA analysis, we found in enriched canonical pathways, CCL1, CCL3, CXCL8 and especially Th17 cells and its related cytokines IL-17A and IL-17F were emphasized. In previous publication, Th17 and its cytokines are in ammatory mediator to RPE cells, 20 and RPE cells in in ammatory condition were thought to be a key point to PVR formation. 18 Infection had long been considered as a trigger to immune disease. In the KEGG analysis, to our surprise, rather than immune pathways, infectious related pathways related to amoebiasis, malaria or chagas disease were enriched. As in uveitis, Forrester et al. agreed that infection may directly or indirectly related to noninfectious uveitis in the eye. 21 In this consideration, remind us that infection might be a potential cause of PVR.
LncRNAs are widely expressed in monocytes, macrophages, neutrophils and implicated in the process of in ammation and immunity. Various molecular functions have been ascribed to lncRNAs, including gene regulation in cis, regulation of mRNA stability, and modulation of protein function. In PVR patients, Zhou et al. not only demonstrated that the expression of MALAT1 was signi cantly upregulated in the proliferating membrane, also found MALAT1 was signi cantly up-regulated in the peripheral blood. 6 MALAT1 can inhibit the DNA binding activity of NF-κB, reduce the production of in ammatory cytokines, and down-regulate the autoimmune in ammatory response. The knockdown of MALAT1 can increase lipopolysaccharide (LPS)-induced expression of TNFα and IL-6. 22 However, to our surprise, we did not nd MALAT1 were upregulated in PBMC in PVR patient.
In our research, by analysis the mRNA and lncRNA co-expression in PBMC, we found immune-related gene NFKBIA, and chemokines CXCL2 and CXCL8 and their associate LncRNA AC007032.1, AC037198.2, AL929472.2, SLED1 were highly associated with PVR. AC007032.1 is associated with immunomodulatory cytokine Nampt, 23 while SLED1 was found up regulated in peripheral blood cells of systemic lupus erythematosus patients. 24 Additionally, within IPA analysis, virus infection relate pathways and immune relate networks were highlighted in PVR patients.
Further, the most obvious changed LncRNA AC037198.2 and ZNF433-AS1 were selected to verify and were proved their actual change in PBMC from PVR patients. LncRNA AC037198.2 is associate with THBS1(thrombospondin 1) gene, which encoded a secreted protein to mediate cell-to-cell and cell-to-matrix interactions. As for ZNF433-AS1, this LncRNA can suppress ZNF433, which belongs to transcriptional factors with the zinc nger motif, and was found that play an important role in multiple sclerosis, which is an autoimmune disease. 25

Conclusions
In summary, we provide a landscape of differential expression pro le of mRNAs and lncRNAs between PVR and controls and construct an mRNA-lncRNA co-expression network based on the DETs. Pathway enrichment analyses offer novel insights into the pathogenesis of PVR, indicating that PVR might be related with abnormal immune system or previous infection. More importantly, some deferentially expressed lncRNAs, like LncRNA-AC037198.2 and ZNF433-AS1were appeared in our study, which might be potential molecular signatures for PVR. These results will provide hence our understanding of this disease and provide novel therapeutic targets for PVR patients.