Transcriptome expression profiles associated with diabetic nephropathy development

Abstract

The objective of this study was to identify different transcriptome expression profiles involved in the pathogenesis of diabetic nephropathy (DN) and to illustrate the diagnostic and therapeutic potential of mRNAs, long noncoding RNAs (lncRNAs), and circular RNAs (circRNAs) in DN progression. The participants were divided into four groups: normoalbuminuria (group DM), microalbuminuria (group A2), macroalbuminuria (group A3) and healthy controls (group N). There were three individuals in each group for sequencing. Transcriptome sequencing analysis was performed on the peripheral blood of all the participants to identify the differential expression of mRNAs, lncRNAs, and circRNAs between intervention groups and controls. The functional enrichment analysis, the short time-series expression miner (STEM) program, and the miRNA-circRNA-mRNA network were further conducted. To verify the reproducibility of transcriptome sequencing, 10 and 30 blood samples were collected from the control and diseased groups, respectively. Four candidate biomarkers were selected from differentially expressed circRNAs (circ_0005379, circ_0002024, and circ_0000567, and circ_0001017) and their concentrations in the blood were measured using quantitative PCR (qPCR). In the comparison of A2 with N, 549 mRNAs, 1259 lncRNAs, and 12 circRNAs were screened. In the comparison of A3 with N, 1217 mRNAs, 1613 lncRNAs, and 24 circRNAs were screened. Moreover, in the comparison of diabetes mellitus (DM) with N, 948 mRNAs, 1495 lncRNAs, and 25 circRNAs were screened. Functional enrichment analysis showed that differentially expressed mRNAs were related to insulin secretion, insulin resistance, and inflammation, while differentially expressed lncRNAs were mainly associated with crossover junction endodeoxyribonuclease activity. In STEM analysis, a total of 481 mRNAs and 152 differential expression circRNAs showed a significant tendency. The key relationships in the miRNA-circRNA-mRNA network were identified, such as hsa-miR-103a-3p-circ_0005379-PTEN, hsa-miR-497-5p-circ_0002024-IGF1R and hsa-miR-1269a-circ_0000567-SOX6. In addition, qPCR showed consistent results with RNA sequencing. We found that differentially expressed mRNAs, lncRNAs, and circRNAs participated in DN development. Circ_0005379, circ_0002024, and circ_0000567 could be adopted as potential biomarkers for DN.

Keywords: Diabetic nephropathy; Transcriptome sequencing; circRNAs; lncRNAs; mRNAs.