Transcriptomic clock predicts vascular changes of prodromal diabetic retinopathy

Scientific Reports volume 13, Article number: 12968 (2023) Cite this article

494 Accesses

2 Altmetric

Metrics details

Diabetic retinopathy is a common complication of long-term diabetes and that could lead to vision loss. Unfortunately, early diabetic retinopathy remains poorly understood. There is no effective way to prevent or treat early diabetic retinopathy until patients develop later stages of diabetic retinopathy. Elevated acellular capillary density is considered a reliable quantitative trait present in the early development of retinopathy. Hence, in this study, we interrogated whole retinal vascular transcriptomic changes via a Nile rat model to better understand the early pathogenesis of diabetic retinopathy. We uncovered the complexity of associations between acellular capillary density and the joint factors of blood glucose, diet, and sex, which was modeled through a Bayesian network. Using segmented regressions, we have identified different gene expression patterns and enriched Gene Ontology (GO) terms associated with acellular capillary density increasing. We developed a random forest regression model based on expression patterns of 14 genes to predict the acellular capillary density. Since acellular capillary density is a reliable quantitative trait in early diabetic retinopathy, and thus our model can be used as a transcriptomic clock to measure the severity of the progression of early retinopathy. We also identified NVP-TAE684, geldanamycin, and NVP-AUY922 as the top three potential drugs which can potentially attenuate the early DR. Although we need more in vivo studies in the future to support our re-purposed drugs, we have provided a data-driven approach to drug discovery.

Diabetic retinopathy is a common diabetes complication of long-term diabetes that results in damage to the retinal blood vessels and can progress towards vision loss. In several developed countries, including the United States, diabetic retinopathy is the leading cause of new cases of blindness among adults aged 20 to 741. People living with diabetes, regardless of clinical subtypes, have roughly one-third of a chance to develop diabetic retinopathy2, similarly affecting both males and females in general3,4.

At the early stages, non-proliferative diabetic retinopathy (DR) can be detected by screening when vision is not yet affected. However, there is no effective way to slow down or stop the progression of the disease to prevent vision loss, where patients are frequently left untreated until later stages. Moreover, diabetic retinopathy is not reversible and will continue to progress. Eventually, some cases will develop into proliferative diabetic retinopathy, where fibrosis develops on the retinal surface, and abnormal hemorrhage-prone new vessels lead to vision loss5. To prevent vision loss from diabetic retinopathy, we need to understand the early pathogenesis of asymptomatic retinopathy. Unfortunately, these early processes remain poorly understood partly due to insufficient histological material from human donors and limited animal models which can recapitulate the natural progression of human diabetic retinopathy6.

Nile rat (Arvicanthis niloticus), native to Northern Africa, is a diurnal rodent with a general diet of grass stems7. However, when fed a conventional rodent chow in a laboratory environment, the Nile rat quickly develops diet-induced diabetes in both sexes. To clarify, conventional rodent chow, though meant to sustain healthy common laboratory rodents, is hypercaloric, relative to the native diet, for the Nile rat species and will be thereafter described as a high-caloric diet. To prevent diet-induced diabetes in a laboratory environment, the Nile rat can be fed a high-fiber diet8. Thus, diet-induced diabetes in the Nile rat appear to follow the natural progression of type 2 diabetes in humans. Our previous study demonstrated that diabetic Nile rats could develop retinopathy with advanced retinal lesions similar to clinical features linked to vision loss in humans, such as macular edema, capillary non-perfusion, and proliferative diabetic retinopathy9. Specifically, Nile rat lesions included retinal leakage, retinal thickening, capillary non-perfusion, intraretinal microvascular abnormalities, and apparent proliferative diabetic retinopathy, based on observed neovascularization beyond the retinal surface9. Typically, other rodent models with diet-induced diabetic retinopathy do not exhibit the exact microvascular damage seen in humans and often lack proliferative neovascularization. However, neovascularization was observed in relatively older Nile rats and age might play a role in neovascularization. In that study, we also demonstrated that elevated acellular capillary density is a reliable quantitative trait present in the initial development of retinopathy, preceding additional vascular dysfunction in pericytes and endothelial cells9. The use of retinal acellular capillary density as one of the earliest markers of retinopathy is consistent in both Nile rat9 and in humans10. Therefore, we reasoned that retinal vasculature change is an important causal factor in initiating and driving the pathogenesis of diabetic retinopathy. Since retinal acellular capillary density is one of the initial detectable changes relevant to retinal vascular dysfunction, we could use acellular capillary density to mark early disease progression in diabetic retinopathy.

In this study, we look for gene signatures that correspond to the progression of increasing acellular capillary density in the vascular network of the whole retina. We performed RNA-Seq on freshly isolated retinal vascular tissue from 28 male and female Nile rats on either a high-caloric or high-fiber diet. We knew from our previous study that the acellular capillary density is a bilateral phenotype where left and right eyes are highly correlated within individual Nile rats9. Taking advantage of the bilateral nature of this phenotype, we could use one eye to determine the acellular capillary density that required fixed tissue and the contralateral eye to harvest fresh tissue for RNA sequencing. We identified gene signatures associated with increasing acellular capillary density, a quantitative marker of early retinopathy. Using a random forest machine learning model, we identified a combination of fourteen genes that can quantitatively predict acellular capillary density that can serve as a vascular-specific transcriptomic clock denoting the progression of retinopathy.

Leveraging a large-scale chemical compound/drug perturbated gene expression signature database (LINCS L1000), we identified NVP-TAE684, geldanamycin, and NVP-AUY922 as the top three potential drugs which can reverse the acellular capillary density correlated genes signatures, thereby providing putative drug candidates to address the early stages of retinopathy. Interestingly, both geldanamycin and NVP-AUY922 are Hsp90 inhibitors11,12. Studies have shown that Hsp90 inhibitors can prevent retinal degeneration in models of retinitis pigmentosa and age-related macular degeneration13. Although our re-purposed candidate drugs require further validation in vivo, we have provided a novel data-driven approach to prioritize compounds that can potentially block the progression of retinopathy prior to permanent vision impairment. Our study also reveals retinal vascular transcriptomics associated with the early pathogenesis of retinopathy.

All animal experiments were approved by the University of California Santa Barbara, Institutional Animal Care and Use Committee, and conducted in accordance with the NIH Guide for the Care and Use of Laboratory Animals; study protocol 893. UCSB founder Nile rats were derived from the Brandeis University colony of the KC Hayes Laboratory. Nile rats in UCSB are housed in a 12-h light cycle (10 a.m. to 10 p.m.) room at 21–26 °C in a conventional facility with individually ventilated cages and are provided autoclaved Sanichips as bedding material, with additional cage enrichment, including wooden blocks and plastic bones. Insulin was not administered to the diabetic rats. To avoid serious adverse effects from diabetes, we euthanized all animals with over 5% weight loss or a large drop of over 100 mg/dL in RBG. To identify these animals, we do a weekly check of the weight and RBG of all Nile rats after they reached RBG ≥ 500 mg/dL.

Nile rats are either fed a high-fiber diet (Lab Diet 5L3M; Newco Speciality, Rancho Cucamonga, CA, USA) or a high-caloric diet (Formulab Diet 5008; Newco Speciality, Rancho Cucamonga, CA, USA)14. The percentage of crude fiber is 23% for the high-fiber diet and only 4% for the high-caloric diet. The ratio for the percentage of calories provided by carbohydrate, fat, and protein were 67:10:23 for the high-fiber diet and 56:17:27 for the high-caloric diet. Random blood glucose (RBG) levels are measured using a glucometer (Contour) every four weeks (between 2 and 5 p.m.) starting at weaning age (4 weeks old) with blood from a tail prick. There were 32 Nile rats in total involved in this study. Four out of thirty-two Nile rat samples failed at the RNA-seq library step (whole retinal vasculature). Thus, we have 28 Nile rats in total for RNA-seq data analysis. The detailed Nile rat information, such as acellular capillary density, average random blood glucose, duration of diabetes, sex, age, diet, and whether it is included in the RNA-seq data analysis, is shown in Supplemental Table 1. All animals are monitored for signs of development of diabetic complications that would be considered distressful to the animal, which include but are not limited to dehydration (severe polydipsia and polyuria, an indication of kidney failure), or pronounced and sustained weakness or lethargy, when we will apply the humane endpoint to these distressed animals. Euthanasia is carried out by carbon dioxide inhalation. This study is reported in accordance with ARRIVE guidelines (https://arriveguidelines.org).

From each animal, we did a trypsin digest on one retina and collected the retinal vasculature for RNA-Seq from the other retina. The retinal vasculature was loosened from other retinal cells using osmotic lysis, where the retina was placed in RNase-free water for 1 h at 4 °C on a shaker. Next, water was aspirated, and the remaining retinal tissue was incubated in a 10% DNase solution at 37 °C for 15 min. Then, the retina was transferred to a petri-dish with RNase-free water. Working under a dissecting microscope, we used a P200 pipet to squirt water on the retina until the inner limiting membrane was detached and neuronal cells were washed away from the retinal vasculature. Finally, the resulting retinal vasculature is preserved in RNAlater overnight at 4 °C and stored at − 80 °C. For the trypsin digest, we used our previously published method9.

Total RNAs from 32 Nile rat vascular samples were isolated using trizol (ThermoFisher #15596018) and chloroform phase separations followed by the RNeasy mini protocol (Qiagen #74106) with optional on-column DNase digestion (Qiagen #79254). One hundred nanograms of total RNA was used to prepare sequencing libraries using the LM-Seq (Ligation Mediated Sequencing) protocol15. RNAs were selected using the NEB Next Poly A + Isolation Kit (NEB #E7490S/L). Poly A + fractions were eluted, primed, and fragmented for 7 min at 85 °C. First stand cDNA synthesis was performed using SmartScribe Reverse Transcriptase (Takara Bio USA #639538), and RNA was removed. cDNA fragments were purified with AMpure XP beads (Beckman Coulter #A63881). The 5′ adapter was ligated, and 18 cycles of amplification were performed. These final indexed cDNA libraries were quantified, normalized, multiplexed, and run as single-end reads on the HiSeq 3000 sequencer. Four out of thirty-two Nile rat RNA-seq samples were failed during the library preparation step, and thus were excluded in the RNA-seq data analysis. Hence, there are 28 Nile rat RNA-seq samples in the downstream data analysis.

We used our lab-developed CRSP tool16 to estimate the gene expression values of Nile rat retina vascular tissues. Specifically, the retinal vasculature RNA-seq reads were mapped to our lab pre-built Nile rat transcriptomic contigs, assembled via Trinity tool17 based on ultra-deep sequencing (2.3 billion × 2 paired-end reads) of a pooled 22 major Nile rat organs16. The transcriptomic contigs were annotated based on mouse proteins via NCBI BLASTX tool18 with E-value < 10–5. A contig level abundance estimation is performed via RSEM19. For each Nile grass rat retina vascular RNA-seq sample, gene-level RNA-seq mapped read counts are obtained by summing the counts of the contigs that are mapped to proteins with the same gene symbol. The relative abundances, in terms of transcripts per million (TPM), for genes were computed by first normalizing each gene’s fragment count by the sum of the “effective lengths” (length less the read length) of the contigs mapped to that gene and then scaling the resulting values such that they sum to one million over all genes. The TPM values were further normalized via the median-of-ratio method by EBSeq20 R package.

Micrographs for quantification were taken with a Canon Rebel XSi digital camera (Canon, Tokyo, Japan) attached to an Olympus CKX41 microscope via an LM scope C-mount. For acellular capillary counts, 12 micrographs were taken from 3 randomly selected areas in each of the 4 retinal quadrants and analyzed at a final magnification of ×200. We categorized abnormal capillary features into three groups: the first group is definitive acellular capillaries; the second group is pericyte processes (described in our paper9) and the third group is ambiguous abnormal features that were very short or sometimes partially obscured by other vessels. The total area used for quantification corresponded to approximately 8.5% of the whole retinal area. The counts were quantified using FIJI computer software without any additional plugins by two masked graders.

To detect gene expression patterns associated with acellular capillary density levels, we performed a segmented regression analysis. Segmented regression is a method in which the independent variable is partitioned into intervals and a separate linear regression to fit each interval. For each gene, the number of breakpoints was determined by the lowest Bayesian Information Criterion (BIC)21 via enumerating all K (K <  = N) a possible number of breakpoints, where N was the number of time points. An optimal model will be found for every gene. We further defined the goodness of fit of the optimal model as the adjusted correlation coefficient R2:

where is the optimally chosen K for a specific gene and N is the number of time points. To further estimate the false discovery rate (FDR), we permutated the data and assessed the background distribution of R2. We only consider genes with FDR < 5% as a good fitting. The segmented regression was implemented by Trendy22 R package.

Gene ontology (GO) enrichment analysis was performed using the R package (“allez”)23. The enrichment p-values were further adjusted by Benjamini–Hochberg (BH) multiple test correction.

We used the random forest regression (RFR) model to predict the retina acellular capillary density from gene expression patterns. RFR is a supervised machine learning algorithm that uses an ensemble learning method for regression. The ensemble learning method is a technique that combines predictions from multiple machine learning models to make a more accurate prediction than a single model. RFR constructs a multitude of decision trees at training time. It takes advantage of two machine-learning techniques: bagging24 and random feature selection. In bagging, each tree is trained on a bootstrap sample of the training data, and aggregating prediction results make predictions from trees. RFR randomly selects a subset of features to split at each node when growing a tree. To assess the prediction performance of the random forest algorithm, RFR performs a type of cross-validation in parallel with the training step by using the so-called out-of-bag (OOB) samples. Specifically, each tree is grown using a particular bootstrap sample in training. Since bootstrapping is sampling with replacement from the training data, some sequences will be ‘left out’ of the sample, while others will be repeated in the sample. The ‘left out’ sequences constitute the OOB sample. On average, each tree is grown using about 1 − e−1 ≅ 2/3 of the training sequences, leaving e−1 ≅ 1/3 as OOB. Because OOB sequences have not been used in the tree construction, one can use them to estimate the prediction performance25. The RF algorithm was implemented by the “randomForest” R package26. We first ranked the gene importance based on the permutation test to achieve the best prediction performance via a minimum number of genes. Specifically, we permutated each gene's expression patterns and calculated the prediction error changes for each permutation. Intuitively, if a gene is essential, the permutation will result in a significant decrease in the prediction performance. Hence, we can rank the gene importance based on the prediction error changes for each permutation. Then, we used a forward feature selection method (i.e., adding each gene at a time, starting with the top 3 most essential genes) to select a minimum number of genes that can achieve the best performance. The prediction performance was calculated by Pearson correlation coefficient (R) between observed and predicted acellular capillary densities. We identified that combining the top 14 genes achieved the best Pearson correlation coefficient (R). Hence, we used the gene expression patterns of the top 14 ranked genes as features to build an RFR model to predict acellular capillary density.

Increased acellular capillary density is one of the earliest detectable signs of diabetic retinopathy9,10. Acellular capillaries are capillary segments with no nuclei along their length, as shown in Fig. 1 (Nile rat retina acellular capillaries; indicated as arrows). Hyperglycemia is a well-known risk factor for increased acellular capillaries9,10. However, it is unclear how acellular capillary density is related to three relevant factors of type 2 diabetes: blood glucose levels, diet, and sex. First, we performed a Spearman's rank correlation analysis (N = 32; Nile rats; See Supplementary Table 1) to investigate the correlation between random blood glucose (RBG) level and acellular capillary density. As shown in Fig. 2A, the RBG is moderately associated with acellular capillary density (Spearman’s correlation coefficient Rho = 0.31, P-value = 0.09), indicating a positive but complex relationship between RBG and acellular capillary density, where the trend is more obvious at the pre-diabetic range of RBG. We also observed that some Nile rats with high acellular capillary density do not have hyperglycemia. Next, we investigated the role of diet. Nile rats on a high-caloric diet have higher RBG compared to Nile rats on a high-fiber diet (Fig. 2B, Wilcoxon rank sum exact test, P-value = 1.38e-05). In our study cohort, 61.9% of Nile rats on the high-caloric diet had not yet developed hyperglycemia, defined as RBG > 100 mg/dL9. Hence, we could investigate if diet could affect acellular capillary density independent of RBG. Additionally, since retinopathy can have differences between males and females, we also included sex in our acellular capillary density analysis. To understand individual factors contributing to increased acellular capillary density, we developed a Bayesian network (BN) to calculate the conditional probability of high acellular capillary density based on RBG, diet, and sex. As shown in Fig. 2C, the likelihood of a high acellular capillary density is progressively increased when factors of blood glucose level, diet, and sex are combined. For example, a Nile rat on a high-caloric diet will have a 43.07% chance of developing acellular capillary density greater than 10 counts per mm2. This probability increases to 62.0% if this Nile rat also has an RBG > 100 mg/dL. Furthermore, if this Nile rat is male, the possibility will further increase to 80.56%. Hence, a male Nile rat on a high-caloric diet with a high blood glucose level will have a much higher chance of presenting high acellular capillary density, a quantitative marker for early retinopathy.

Demonstration of Nile rat retina acellular capillaries and the overall experimental design for RNA-seq. Arrows indicate acellular capillaries. Acellular capillaries are capillary-sized vessel tubes without nuclei which are considered a key symptom of the early stage of retinopathy.

Factors associated with increased acellular capillaries density (ACD). (A) Random blood glucose (RBG) is moderately associated with ACD (Spearman's rank correlation coefficient Rho = 0.31, P-value = 0.09); (B) High-caloric diet is significantly associated with high random blood glucose (RBG) (Wilcoxon rank sum exact test, P-value = 1.38e-05); (C) Bayesian network modeling the conditional probability of relative higher ACD (> 10 counts per mm2) suggests that combined factors of high RBG, high-caloric diet, and the male will result in a twofold increase the likelihood of high ACD than any factors alone.

We hypothesized that the retinal vasculature contains molecular signals driving the early stages of retinopathy pathogenesis. To test our hypothesis, we performed RNA-seq on these whole retinal vascular tissues (See methods) of Nile rats. We have 32 Nile rat in total with a large acellular capillary density dynamic range from 2.9 to 17.9 (median: ~ 10) counts per mm2. However, four Nile rat samples failed to pass quality control (QC) during the RNA-seq library step and thus were excluded from the RNA-seq data analysis (Supplementary Table 1). Hence, there are 28 Nile rat retinal vascular RNA-seq data in total. The gene expression values of retinal vasculatures were quantified as normalized TPM (transcript per million) values. We performed a segmented regression analysis to identify gene signatures associated with elevated acellular capillary density (see Methods). As shown in Fig. 3, we identified four major gene expression pattern categories associated with elevated ACD: Up (gene expression up-regulated); Down (gene expression Down-regulated); Up-Down (gene expression first goes up and then followed by down-regulation); Down-Up (gene expression first goes down and then followed by up-regulation). Examples of each pattern can be seen in Fig. 3B. Interestingly, the gene ontology (GO) term (Astrocyte development) is the top one enriched in genes with the Up pattern. Studies have shown that glucose tightly controls the morphological functions of astrocytes27, suggesting a direct or indirect impact of hyperglycemia on retinal vasculature transcriptomic changes. The down-regulated genes (negatively correlated with acellular capillary density) are enriched in DNA and cellular repair/response-related functions, such as inter-strand cross-link repair. It has been known that macular (a light-sensitive layer in the center of the retina) degeneration is tightly associated with impaired retina DNA damage and repair mechanism28,29, which is characterized by the sensitivity of retinal pigment epithelial (RPE) to blue and UV lights exposure. Interestingly, the impaired DNA damage repair mechanism is observed in both macular degeneration and early retinopathy. The Up-Down and Down-Up patterns are associated with a wide range of biological functions (See Supplementary Table 2), such as electron transport chain (Up-Down pattern) and regulation of neuron maturation (Down-Up) patterns, indicating the complexity of dynamic transcriptomic changes in response to elevated acellular capillary density. The top enriched gene ontology terms for each pattern are shown in Fig. 3C.

Different gene expression patterns are associated with elevated acellular capillaries density (ACD). (A) Segmented regression analysis identified four patterns (Up, Down, Up-Down and Down-Up) associated with elevated ACD with a False Discovery Rate (FDR) < 5%; (B) Examples of each pattern; (C) Top 3 enriched gene ontology (GO) terms for each pattern.

We hypothesize that the pathological characteristics of a dysfunctional retina are encoded in their transcriptomic patterns. Thus, subtle changes associated with the early progression of diabetic retinopathy, as quantified by acellular capillary density, can be inferred from a combination of gene expression patterns. To identify a minimum number of genes that can be used to impute acellular capillary density, we develop a random forest regression (RFR) model. We ranked the gene importance based on the prediction performance changing by permutation test (i.e., permutating gene expression patterns for each gene). Then we used a forward feature selection method (i.e., adding each gene at a time) to select a minimum number of genes that can achieve the best performance (See Methods). As shown in Fig. 4A, we identified a combination of 14 genes that can predict acellular capillary density with high accuracy (Pearson correlation coefficient R = 0.91; Fig. 4B). It suggests that acellular capillary density is uniquely encoded in gene expression patterns of these 14 genes, and a machine-learning regression model can decode the transcriptomic signals associated with this early retinopathy feature.

Random forest regression model to predict acellular capillaries density. (A) Random forest regression model identified a combination of 14 genes that can quantitively predict acellular capillaries density. (B) The predicted acellular capillaries density is highly correlated with the observed acellular capillaries density.

Currently, there are no effective ways to treat or reverse early retinopathy. Emerging evidence suggests that the reversal of disease gene expression patterns correlates well with drug efficacy and reveals therapeutic targets30. Specifically, if a small molecule or compound can reverse or partially reverse disease marker gene expression patterns, it can be a potential disease treatment drug. Given that elevated acellular capillary density is one of the earliest quantitative features of diabetic retinopathy, we hypothesize that if a compound can reverse or partially reverse the acellular capillary density-associated transcriptome, it could be a putative drug candidate to alter the progression toward vision loss. To this end, we identified genes that are positively or negatively correlated with acellular capillary density levels via Spearman’s rank correlation analysis (Absolute Correlation coefficient (Rho) > 0.4; Fig. 5). Absolute Rho values greater than 0.4 are generally considered to have correlated patterns31. We then searched LINCS L1000 chemical compound perturbated gene expression database32 to find drugs that could potentially reverse this pattern (i.e., a compound that can down-regulate the gene expression for acellular capillary density positively correlated genes while up-regulating the gene expression for acellular capillary density negatively correlated genes). We found NVP-TAE684, geldanamycin, and NVP-AUY922 to be the top three compounds identified to potentially reverse the retinal vascular phenotype associated with high acellular capillary density (Fig. 5). NVP-TAE684 is an anaplastic lymphoma kinase (ALK) inhibitor. ALK is a receptor tyrosine kinase belonging to the insulin receptor superfamily sharing a high degree of homology with leukocyte tyrosine kinase (LTK)33. Both geldanamycin and NVP-AUY922 are Hsp90 inhibitors11,12. Interestingly, studies have shown that a short-time treatment of Hsp90 inhibitor can stimulate a retina stress response with molecular chaperone expression, preventing retinal degeneration in models of retinitis pigmentosa and age-related macular degeneration13. Although no study shows that Hsp90 inhibitors can be used to treat early retinopathy, there is evidence that Hsp90 inhibitors, such as AUY922, can protect and repair microvascular endothelial cells related to endothelial barrier dysfunction in a different biological context (e.g., human lung34). Retina endothelial dysfunction is a significant complication that drives diabetic retinopathy35,36. Hence, it is likely that Hsp90 inhibitors (geldanamycin and NVP-AUY922) can potentially treat or protect early retinopathy via modulating microvascular endothelial cells.

Identification of putative candidates that can reduce acellular capillary density via reverse transcriptome. The heatmap shows genes with |Rho|> 0.4 between gene expression values and acellular capillaries density (ACD). Searching LINCS L1000 chemical compound perturbated gene expression database identified NVP-TAE684, geldanamycin, and NVP-AUY922 as the top three candidates.

To further evaluate our predicted drugs, we performed a text-mining approach. We listed a set of retinopathy-associated anatomical features: Microaneurysm, Pericyte loss, Dot haemorrhages, Intraretinal haemorrhages, Hard exudates, Cotton wool spots, Neovascularization, Venous beading, Ischemia. We then used a text mining tool (KinderMiner: https://www.kinderminer.org/login) to calculate whether the predicted compounds (NVP-TAE684, geldanamycin and NVP-AUY922) are statistically enriched in any of the retinopathy-associated anatomical features as listed above via scanning 34,006,217 PubMed articles). We found geldanamycin is statistically enriched in “Neovascularization term” (P-value = 0.001) and “Ischemia” term (P-value = 8.7E-10). Since both geldanamycin and NVP-AUY922 are Hsp90 inhibitors, we further calculated the enrichment P-values of the term “Hsp90 inhibitors” in retinopathy-associated anatomical features based on KinderMiner tool. The term “Hsp90 inhibitors” is statistically enriched in “Neovascularization term” (P-value = 0.02) and “Ischemia” term (P-value = 6.6E-5). This indirectly supported that our predicted Hsp90 inhibitors (geldanamycin and NVP-AUY922) could potentially attenuate or mitigate prodromal diabetic retinopathy.

Diabetic retinopathy is a progressive disease that will eventually lead to vision loss. Although current technological advances have led to easy detection of early retinopathy prior to any vision impairment, strategies to address early retinopathy and prevent vision loss are lacking. Human tissues presenting early retinopathy are scarcely available. Hence animal models are needed. Diabetic Nile rats develop retinal lesions similar to human diabetic retinopathy associated with vision loss9, which may arise from shared molecular mechanisms.

In our previous study, we found that acellular capillary density is a quantitative feature associated with risk factors of diabetic retinopathy at least in the Nile rat model. The acellular capillary density is apparent earlier with a large dynamic range compared to other vascular characteristic changes such as pericyte loss (small dynamic range) and neuronal lesions (apparent later) in the Nile rat. We hypothesize that molecular signatures initiating early retinopathy are encoded in the retinal vasculature and therefore are also associated with acellular capillary density. We developed a machine learning model via gene expression patterns of fourteen genes that can quantitively predict the acellular capillary density. Since acellular capillary density is a quantitative feature associated with early diabetic retinopathy, our model can be used as a retinopathy early progression “transcriptomic clock” measuring the severity or the progression of early retinopathy. However, we acknowledge that this model was based on Nile rat diabetic retinopathy data. Comparisons with other animal models and humans are needed to validate the gene expression pattern and if this prediction model is consistent across species.

Emerging evidence suggests that reversal of disease gene expression patterns correlates well with drug efficacy30. Hence, we hypothesize that if a small molecule or a drug can reverse the majority or a significant fraction of the DR-driven gene expression changes, it can be a putative DR treatment drug candidate. We used acellular capillary density-associated up- and down-regulated gene signatures to identify drugs that can potentially reverse or partially reverse such patterns via searching LINCS L1000 chemical compound perturbated gene expression database32. The fundamental basis for this approach assumes that targeting early vascular changes is the key to blocking early retinopathy. This assumption is based on our prior study that the elevated acellular capillary density was a reliable quantitative trait present in the initial development of retinopathy preceding additional vascular dysfunction in pericytes and endothelial cells9. Among the top three predicted drugs, geldanamycin, and NVP-AUY922 are Hsp90 inhibitors11,12. Some studies have shown that Hsp90 inhibitors can protect and repair microvascular endothelial cells related to endothelial barrier dysfunction but in a different biological context: human lungs34. Although we need more in vivo studies in the future to support our findings, we provide a data-driven approach to re-purpose treatment drugs, which can largely shorten the traditional drug discovery path. However, we also acknowledge that in vivo and in vitro studies are needed to validate these re-purposed drugs further.

Our study not only gives transcriptomic insights into the pathogenesis of early-stage retinopathy but also provides a novel data-driven approach to prioritize and re-purpose treatment drugs.

We have submitted the RNA-seq data and the gene expression data (TPMs and Mapping counts) to GEO (accession number: GSE220672).

Diabetic retinopathy

Bayesian network

Random blood glucose

Acellular capillary density

Gene ontology

Random forest regression

Out-of-bag

Retinal pigment epithelial

Klein, R. & Klein, B. E. Vision disorders in diabetes. Diabetes Am. 1, 293 (1995).

Google Scholar

Zhang, X. et al. Prevalence of diabetic retinopathy in the United States, 2005–2008. JAMA 304, 649–656. https://doi.org/10.1001/jama.2010.1111 (2010).

Article CAS PubMed PubMed Central Google Scholar

Nakayama, Y. et al. Retrospective exploratory analyses on gender differences in determinants for incidence and progression of diabetic retinopathy in Japanese patients with type 2 diabetes mellitus. Endocr. J. 68, 655–669. https://doi.org/10.1507/endocrj.EJ20-0630 (2021).

Article CAS PubMed Google Scholar

de Ritter, R. et al. Sex differences in the association of prediabetes and type 2 diabetes with microvascular complications and function: The Maastricht study. Cardiovasc. Diabetol. 20, 102. https://doi.org/10.1186/s12933-021-01290-x (2021).

Article CAS PubMed PubMed Central Google Scholar

Davis, M. D. et al. Risk factors for high-risk proliferative diabetic retinopathy and severe visual loss: Early Treatment Diabetic Retinopathy Study Report #18. Invest. Ophthalmol. Vis. Sci. 39, 233–252 (1998).

CAS PubMed Google Scholar

Olivares, A. M. et al. Animal models of diabetic retinopathy. Curr. Diab. Rep. 17, 93. https://doi.org/10.1007/s11892-017-0913-0 (2017).

Article PubMed PubMed Central Google Scholar

Senzota, R. A case of rodent: Ungulate resource partitioning. J. Mammal. 64, 326–329 (1983).

Article Google Scholar

Toh, H., Thomson, J. A. & Jiang, P. Maternal high-fiber diet protects offspring against type 2 diabetes. Nutrients 13, 94. https://doi.org/10.3390/nu13010094 (2020).

Article CAS PubMed PubMed Central Google Scholar

Toh, H. et al. Vascular changes in diabetic retinopathy-a longitudinal study in the Nile rat. Lab. Invest. https://doi.org/10.1038/s41374-019-0264-3 (2019).

Article PubMed PubMed Central Google Scholar

Tan, F. et al. Associated risk factors in the early stage of diabetic retinopathy. Eye Vis. (Lond.) 6, 23. https://doi.org/10.1186/s40662-019-0148-z (2019).

Article PubMed Google Scholar

Eccles, S. A. et al. NVP-AUY922: A novel heat shock protein 90 inhibitor active against xenograft tumor growth, angiogenesis, and metastasis. Cancer Res. 68, 2850–2860. https://doi.org/10.1158/0008-5472.CAN-07-5256 (2008).

Article CAS PubMed Google Scholar

Ochel, H. J., Eichhorn, K. & Gademann, G. Geldanamycin: the prototype of a class of antitumor drugs targeting the heat shock protein 90 family of molecular chaperones. Cell Stress Chaperones 6, 105–112. https://doi.org/10.1379/1466-1268(2001)006%3c0105:gtpoac%3e2.0.co;2 (2001).

2.0.co;2" data-track-action="article reference" href="https://doi.org/10.1379%2F1466-1268%282001%29006%3C0105%3Agtpoac%3E2.0.co%3B2" aria-label="Article reference 12" data-doi="10.1379/1466-1268(2001)0062.0.co;2">Article CAS PubMed PubMed Central Google Scholar

Aguila, M. & Cheetham, M. E. Hsp90 as a Potential therapeutic target in retinal disease. Adv. Exp. Med. Biol. 854, 161–167. https://doi.org/10.1007/978-3-319-17121-0_22 (2016).

Article CAS PubMed PubMed Central Google Scholar

Chaabo, F., Pronczuk, A., Maslova, E. & Hayes, K. Nutritional correlates and dynamics of diabetes in the Nile rat (Arvicanthis niloticus): A novel model for diet-induced type 2 diabetes and the metabolic syndrome. Nutr. Metab. (Lond.) 7, 29. https://doi.org/10.1186/1743-7075-7-29 (2010).

Article CAS PubMed Google Scholar

Hou, Z. et al. A cost-effective RNA sequencing protocol for large-scale gene expression studies. Sci. Rep. 5, 9570. https://doi.org/10.1038/srep09570 (2015).

Article CAS PubMed PubMed Central Google Scholar

Toh, H. et al. A Nile rat transcriptomic landscape across 22 organs by ultra-deep sequencing and comparative RNA-seq pipeline (CRSP). Comput. Biol. Chem. 102, 107795. https://doi.org/10.1016/j.compbiolchem.2022.107795 (2022).

Article CAS PubMed Google Scholar

Grabherr, M. G. et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29, 644–652. https://doi.org/10.1038/nbt.1883 (2011).

Article CAS PubMed PubMed Central Google Scholar

Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. Basic local alignment search tool. J. Mol. Biol. 215, 403–410. https://doi.org/10.1016/S0022-2836(05)80360-2 (1990).

Article CAS PubMed Google Scholar

Li, B. & Dewey, C. N. RSEM: Accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 12, 323. https://doi.org/10.1186/1471-2105-12-323 (2011).

Article CAS Google Scholar

Leng, N. et al. EBSeq: An empirical Bayes hierarchical model for inference in RNA-seq experiments. Bioinformatics 29, 1035–1043. https://doi.org/10.1093/bioinformatics/btt087 (2013).

Article CAS PubMed PubMed Central Google Scholar

Schwarz, G. Estimating the dimension of a model. In The Annals of Statistics, 461–464 (1978).

Bacher, R. et al. Trendy: Segmented regression analysis of expression dynamics in high-throughput ordered profiling experiments. BMC Bioinf. 19, 380. https://doi.org/10.1186/s12859-018-2405-x (2018).

Article CAS Google Scholar

Newton, M. A., Quintana, F. A., Den Boon, J. A., Sengupta, S. & Ahlquist, P. Random-set methods identify distinct aspects of the enrichment signal in gene-set analysis. Ann. Appl. Stat. 1, 85–106 (2007).

Article MathSciNet MATH Google Scholar

Breiman, L. Bagging predictors. Mach. Learn. 24, 123–140 (1996).

Article MATH Google Scholar

Svetnik, V. et al. Random forest: A classification and regression tool for compound classification and QSAR modeling. J. Chem. Inf. Comput. Sci. 43, 1947–1958. https://doi.org/10.1021/ci034160g (2003).

Article CAS PubMed Google Scholar

Liaw, A. & Wiener, M. Classification and regression by randomForest. R. News 2, 18–22 (2002).

Google Scholar

Lee, C. Y., Dallerac, G., Ezan, P., Anderova, M. & Rouach, N. Glucose tightly controls morphological and functional properties of astrocytes. Front. Aging Neurosci. 8, 82. https://doi.org/10.3389/fnagi.2016.00082 (2016).

Article CAS PubMed PubMed Central Google Scholar

Blasiak, J. & Szaflik, J. P. DNA damage and repair in age-related macular degeneration. Front. Biosci. (Landmark Ed) 16, 1291–1301. https://doi.org/10.2741/3789 (2011).

Article CAS PubMed Google Scholar

Szaflik, J. P. et al. DNA damage and repair in age-related macular degeneration. Mutat. Res. 669, 169–176. https://doi.org/10.1016/j.mrfmmm.2009.06.008 (2009).

Article CAS PubMed Google Scholar

Chen, B. et al. Reversal of cancer gene expression correlates with drug efficacy and reveals therapeutic targets. Nat. Commun. 8, 16022. https://doi.org/10.1038/ncomms16022 (2017).

Article ADS CAS PubMed PubMed Central Google Scholar

Akoglu, H. User’s guide to correlation coefficients. Turk. J. Emerg. Med. 18, 91–93. https://doi.org/10.1016/j.tjem.2018.08.001 (2018).

Article PubMed PubMed Central Google Scholar

Subramanian, A. et al. A next generation connectivity map: L1000 platform and the first 1,000,000 profiles. Cell 171, 1437-1452.e1417. https://doi.org/10.1016/j.cell.2017.10.049 (2017).

Article CAS PubMed PubMed Central Google Scholar

Lemmon, M. A. & Schlessinger, J. Cell signaling by receptor tyrosine kinases. Cell 141, 1117–1134. https://doi.org/10.1016/j.cell.2010.06.011 (2010).

Article CAS PubMed PubMed Central Google Scholar

ColungaBiancatelli, R. M. L., Solopov, P., Gregory, B. & Catravas, J. D. The HSP90 inhibitor, AUY-922, protects and repairs human lung microvascular endothelial cells from hydrochloric acid-induced endothelial barrier dysfunction. Cells 10, 1489. https://doi.org/10.3390/cells10061489 (2021).

Article CAS Google Scholar

Gui, F., You, Z., Fu, S., Wu, H. & Zhang, Y. Endothelial dysfunction in diabetic retinopathy. Front. Endocrinol. (Lausanne) 11, 591. https://doi.org/10.3389/fendo.2020.00591 (2020).

Article PubMed Google Scholar

Siasos, G. et al. Role of endothelial dysfunction and arterial stiffness in the development of diabetic retinopathy. Diabetes Care 38, e9–e10. https://doi.org/10.2337/dc14-1741 (2015).

Article PubMed Google Scholar

Download references

We thank Jennifer Bolin and Jessica Antosiewicz-Bourget for their technical support (RNA-seq). We also thank Xiujie Li for the discussion and suggestions for this project. We thank the support from the Center for Gene Regulation in Health and Disease (GRHD) at Cleveland State University and Prof. Anton Komar.

This study was supported by the Garland Initiative for Vision, funded by the William K. Bowes Jr. Foundation. Peng Jiang is also funded by the Defense Advanced Research Projects Agency (DARPA) (D20AC00002).

Neuroscience Research Institute, University of California Santa Barbara, Santa Barbara, CA, USA

Huishi Toh, Alexander Smolentsev, Ryan Sadjadi, Dennis Clegg & James A. Thomson

Department of Molecular, Cellular and Developmental Biology, University of California Santa Barbara, Santa Barbara, CA, USA

Dennis Clegg & James A. Thomson

Morgridge Institute For Research, Madison, WI, 53706, USA

Ron Stewart & James A. Thomson

Department of Biological, Geological and Environmental Sciences, Cleveland State University, Cleveland, OH, 44115, USA

Jingqi Yan & Peng Jiang

Center for Gene Regulation in Health and Disease, Cleveland State University, Cleveland, OH, 44115, USA

Jingqi Yan & Peng Jiang

Center for RNA Science and Therapeutics, School of Medicine, Case Western Reserve University, 10900 Euclid Avenue, Cleveland, OH, 44106, USA

Peng Jiang

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

Nile rat experiment: H.T., A.S. & R.S. (Ryan); Data analysis & modeling: P.J.; Technical & lab support: D.C., J.Y., R.S. (Ron); Manuscript writing: H.T., P.J.; Supervision: J.A.T. & P.J.

Correspondence to Peng Jiang.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Toh, H., Smolentsev, A., Sadjadi, R. et al. Transcriptomic clock predicts vascular changes of prodromal diabetic retinopathy. Sci Rep 13, 12968 (2023). https://doi.org/10.1038/s41598-023-40328-w

Download citation

Received: 31 March 2023

Accepted: 08 August 2023

Published: 10 August 2023

DOI: https://doi.org/10.1038/s41598-023-40328-w

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.