Transcriptome analysis identified a novel 3-LncRNA regulatory network of transthyretin attenuating glucose induced hRECs dysfunction in diabetic retinopathy

BACKGROUND: Diabetic retinopathy (DR) is the leading cause of blindness in the working age population. Transthyretin (TTR) showed a significantly decreased concentration in DR patients and exerted a visual protective effect by repressing neovascularization. This work intended to identify long non coding RNAs (lncRNAs) and explore their potential mechanism underlying the protective role of TTR. METHODS: Transcriptome of human retinal endothelial cells (hRECs) treated with low glucose (LG), high glucose (HG) or high glucose with 4 μM TTR (HG + TTR) was conducted. Differentially expressed lncRNAs, mRNAs and TTR related lncRNAs and mRNA were acquired. Functional annotation and Gene Set Enrichment Analysis were applied to analyse TTR affected pathways and processes. Weighted gene co-expression network analysis (WGCNA) was implemented to obtain hub modules and genes. LncRNA-mRNA regulatory networks were constructed based on cis, trans and competing endogenous RNAs acting mode. QRT-PCR was conducted to validate the expression of lncRNAs in aqueous humor and serum samples from 30 DR patients and 10 normal controls. RESULTS: RNA-sequencing of hRECs treated with low glucose (LG), high glucose (HG) or high glucose with 4 μM TTR (HG + TTR) was conducted. 146,783 protein-coding transcripts, 12,403 known lncRNA transcripts and 1184 novel non-coding transcripts were characterized. A total of 11,407 differentially expressed mRNAs (DE-mRNAs), 679 differentially expressed lncRNAs (DE-lncRNAs) in HG group versus LG group, 6206 DE-mRNAs and 194 DE-lncRNAs in HG + TTR versus HG group were obtained, respectively. 853 TTR-mRNAs and 48 TTR-lncRNAs were acquired, and functionally involved in cell cycle, apoptosis, inflammation signalling pathway, response to oxidative stress, neovascularization and autophagy. The WGCNA analysis identified a hub module of 133 genes, with the core function of oxidative stress response, angiogenesis, MAPK pathway, cell proliferation and apoptosis. After qRT-PCR validation, a 3-lncRNA regulatory network was proposed. At last, lncRNAs MSTRG.15047.3 and AC008403.3 showed significantly relative higher expression levels in both aqueous humor and serum samples, compared with normal controls, and FRMD6-AS2 was significantly down-regulated. CONCLUSIONS: TTR regulated mRNAs and biological processes including oxidative stress, inflammation signalling and autophagy. A 3-lncRNA regulatory network was characterized underlying TTR repressing neovascularization, and showed potential diagnostic performance in DR.