Genetic mechanisms of dynamic functional connectivity density in diabetic retinopathy brains: a combined transcriptomic and resting-state functional magnetic resonance imaging study

Abstract

Background:

Diabetic retinopathy (DR) is a condition characterized by fundus lesions resulting from retinal microvascular leakage and obstruction linked to chronic progressive diabetes mellitus. Previous neuroimaging research has revealed both structural and functional changes in the brains of DR patients. Nevertheless, the variations in dynamic functional connectivity density (dFCD) within the brains of DR patients, along with the underlying molecular mechanisms connected to these changes, have yet to be fully understood.

Methods:

Forty-seven diabetic retinopathy (DR) patients and 46 healthy controls (HCs) matched for sex, age, and education were recruited for this study from the Department of Ophthalmology at the Jiangxi Provincial People’s Hospital. All subjects underwent resting-state functional magnetic resonance imaging scans to analyze the differences in dFCD between the two groups. Utilizing the Allen Human Brain Atlas, we conducted spatial correlation analyses integrating transcriptomic and neuroimaging data to pinpoint genes showing correlated expression levels with dFCD alterations in DR patients. Subsequently, we carried out gene enrichment, specific expression, and protein-protein interaction analyses.

Results:

In comparison to the HC group, the DR group exhibited significantly reduced dFCD variability in the left anterior cingulum, left superior occipital gyrus, and right postcentral gyrus. The abnormal dFCD variability is linked to 1,318 positively and 1,318 negatively associated genes, primarily enriched for biological functions such as ion channels, synapses, and cellular junctions. Specific expression analysis revealed that these genes were distinctly expressed in Purkinje neurons, cortex, and striatum brain regions. Furthermore, protein-protein interaction (PPI) analyses indicated that these positive and negative genes could organize PPI networks with the support of respective hub genes.

Conclusion:

our study identified altered dFCD variability in brain regions linked to visual and cognitive functions in DR patients. Moreover, transcriptome-neuroimaging correlation analyses revealed a spatial association between these dFCD changes and the genes with unique functional profiles. These genes were enriched in biologically significant functions and pathways, specific to certain cells and brain areas. Our research offers novel understandings of the genetic mechanisms influencing dFCD alterations in DR.

Introduction

Diabetes mellitus (DM) is a prevalent metabolic disorder associated with various microvascular complications (Xu, 2015; Bahtiyar et al., 2016; Bansal et al., 2023), among these, DR is a significant microvascular complication affecting the retina in diabetic individuals and stands as a leading cause of global blindness (Wong et al., 2016). A survey conducted in 2020 revealed a worldwide DR prevalence of 22.27% in diabetic populations, impacting approximately 103.12 million individuals (Teo et al., 2021). Risk factors for DR include hypertension, dyslipidemia, obesity, pregnancy, and prolonged diabetes duration (Lopes de Faria et al., 2011; Ting et al., 2016; Wondmeneh and Mohammed, 2024). The primary pathological manifestations of DR encompass capillary non-perfusion, vascular leakage, and degeneration, progressing to proliferative retinal detachment and eventual vision loss. Notably, retinal vessels share similarities in anatomy, physiology, and embryology with cerebral vasculature (Kwa et al., 2002; van de Kreeke et al., 2018). Notably, retinal vessels share similarities in anatomy, physiology, and embryology with cerebral vasculature (Hägg et al., 2013; Sanahuja et al., 2016). Individuals with DR face heightened susceptibility to dementia and Alzheimer’s disease (Ciudin et al., 2017; Biessels and Despa, 2018), and potential cognitive impairment (Crosby-Nwaobi et al., 2013; Chai et al., 2024). While evidence indicates brain alterations and abnormal central nervous system (CNS) activity in DR patients, the underlying neurophysiological mechanisms necessitate further exploration.

Resting-state functional magnetic resonance imaging (rs-fMRI), a non-invasive imaging technique, enables the examination of spontaneous intrinsic brain activity by measuring fluctuations in blood oxygen level-dependent (BOLD) signals without specific tasks or stimuli (Fox and Raichle, 2007; Barkhof et al., 2014). Studies using resting-state functional MRI have revealed diverse patterns of spontaneous brain activity alterations in DR patients, including amplitude of low-frequency fluctuations (ALFF), dynamic ALFF (dALFF), regional homogeneity (ReHo), and functional network connectivity (FNC), predominantly observed in brain regions like the middle occipital gyrus, calcarine, anterior cingulate gyrus, posterior cerebellar lobe, and networks associated with visual and cognitive functions (e.g., visual network, default mode network, executive control network, and salience network). These findings shed light on the neural mechanisms underlying visual and central nervous system abnormalities in DR (Huang et al., 2019, 2020, 2021a; Liao X.-L et al., 2019; Qi et al., 2020; Shi et al., 2022). Functional connectivity density (FCD) analysis, an extension of the FC method, assesses the importance of network nodes in functional integration by measuring their total connectivity strength with all other nodes in the brain network (Lv et al., 2018; Tomasi and Volkow, 2010). This approach, without requiring prior assumptions, is suited for exploratory analyses (Tomasi et al., 2016; Pang et al., 2017; Zhang et al., 2017). Moreover, the human brain is a complex dynamic system with non-stationary neural activity and constantly changing neural interactions (Hutchison et al., 2013; Calhoun et al., 2014; Liao W. et al., 2019), implying the importance of investigating dynamic neural communication across different brain regions in DR patients utilizing the dFCD approach. However, the specific alterations in dFCD among DR patients remain unclear.

DR has a hereditary component, with estimated heritability ranging from 6 to 33% (Kim et al., 2022). Multiple studies have demonstrated familial clustering of DR (Zhang et al., 2010; Cho and Sobrin, 2014), with notably high concordance rates among family members; 65% in twins with type 1 diabetes and 95% in twins with type 2 diabetes (Leslie and Pyke, 1982). These findings imply the genetic regulation of DR onset and severity. Genome-wide association studies (GWAS) have the potential to identify genetic risk loci associated with DR. A comprehensive meta-analysis has revealed several gene loci linked to DR that are crucial for cell survival, insulin signaling, and angiogenesis (Grassi et al., 2011). Moreover, various studies have highlighted the influence of genetic factors on resting-state brain function (Fornito et al., 2011). Nevertheless, the genetic underpinnings of neuroimaging changes in DR patients are not effectively captured by GWAS. This challenge is compounded by the fact that most loci identified by GWAS are situated in regulatory, rather than coding, regions of the genome, posing significant challenges for their biological interpretation (Edwards et al., 2013; Tam et al., 2019). Spatial correlation analyses integrating transcriptomics and neuroimaging data from the Allen Human Brain Atlas (AHBA), known for its case-control imaging discrepancies, were conducted to explore the genetic basis of macroscopic neuroimaging alterations (Arnatkeviciute et al., 2022; Xue et al., 2022). This method is widely used to unveil gene expression patterns associated with structural and functional variations in patients’ brains, although investigations on the genetic mechanisms underlying dFCD changes in DR are lacking.

Therefore, In our study, we initially assessed and compared dFCD differences between DR patients and controls. Subsequently, utilizing gene expression data from the AHBA database integrated with dFCD maps, we conducted spatial correlation analyses to identify genes whose expression levels correlated with dFCD changes in DR. Finally, gene enrichment, specific expression, and protein-protein interaction analyses were carried out on the identified genes to elucidate their functions and understand the molecular mechanisms behind dFCD alterations in DR.

Materials and methods

Participants

Forty-seven diabetic retinopathy (DR) patients and 46 healthy controls (HCs) matched for sex, age, and education were recruited for this study from the Department of Ophthalmology at the Jiangxi Provincial People’s Hospital.

The inclusion criteria for diabetic retinopathy (DR) were as follows: (1) Fasting blood glucose levels of equal to or greater than 7.0 mmol/L, random blood glucose levels of equal to or greater than 11.1 mmol/L, or 2-h glucose levels of equal to or greater than 11.1 mmol/L; (2) All DR patients were diagnosed with non-proliferative DR without macular edema. The DR classification followed the Early Treatment Diabetic Retinopathy Study grading scheme, encompassing mild and moderate non-proliferative DR, severe non-proliferative retinopathy, and non-high-risk and high-risk proliferative DR.

The exclusion criteria for diabetic retinopathy (DR) were as follows: (1) Patients with proliferative DR complicated by retinal detachment and vitreous hemorrhage; (2) Presence of additional ocular complications (such as cataracts, glaucoma, high myopia, or optic neuropathy); (3) Patients diagnosed with diabetic nephropathy (urine albumin/creatinine ratio > 30 mg/g persisting for more than 3 months); (4) Diabetic neuropathy.

The inclusion and exclusion criteria for healthy controls (HC) were as follows: (1) Fasting blood glucose levels < 7.0 mmol/L, random blood glucose levels < 11.1 mmol/L, glycosylated hemoglobin levels < 6.5%; (2) Absence of ocular diseases (such as myopia, cataracts, glaucoma, optic neuritis, or retinal degeneration); (3) Binocular visual acuity ≥ 1.0; (4) No history of ocular surgery; (5) Absence of psychiatric disorders. All participants met the eligibility criteria for MRI scans.

The study protocol adhered to the Declaration of Helsinki and received approval from the Research Ethics Committee of Jiangxi Provincial People’s Hospital. Before participating in the study, all participants were informed about the study’s objectives, methods, and potential risks, and provided written informed consent to participate.

MRI acquisition

Functional images were acquired using a gradient-echo-planar imaging sequence. Participants were instructed to maintain a relaxed state with closed eyes, refrain from focusing on specific thoughts, and avoid falling asleep. Whole-brain T1-weighted images were obtained through a three-dimensional brain volume imaging (3D-BRAVO) MRI protocol with the following parameters: repetition time (TR)/echo time (TE) = 8.5/3.3 ms, slice thickness = 1.0 mm, no intersection gap, acquisition matrix = 256 × 256, field of view = 240 × 240 mm², and flip angle = 12°. Functional images were acquired using a gradient-echo-planar imaging sequence with the following parameters: TR/TE = 2,000 ms/25 ms, slice thickness = 3.0 mm, gap = 1.2 mm, acquisition matrix = 64 × 64, flip angle = 90°, field of view = 240 × 240 mm², voxel size = 3.6 × 3.6 × 3.6 mm³, and 35 axial slices. All participants were instructed to relax with closed eyes, avoiding specific thoughts or falling asleep.

Data preprocessing

The preprocessing steps utilized the Data Processing and Analysis of Brain Imaging (DPABI) toolbox (DPABI)¹ (Yan et al., 2016), a tool based on Statistical Parametric Mapping (SPM12) incorporated in MATLAB 2013a (MathWorks, Natick, MA, United States). The pre-processing protocol consisted of the following steps: (1) Removal of the initial 10 volumes, (2) Correction for slice timing effects and motion, excluding participants with head motion exceeding 2 mm or rotation exceeding 1.5° during scanning (Van Dijk et al., 2012). (3) Normalization of data to Montreal Neurological Institute (MNI) 152 space at a resolution of 3 mm × 3 mm × 3 mm, and (4) Spatial smoothing through convolution with a 6 mm × 6 mm × 6 mm isotropic Gaussian kernel at full width at half maximum.

Dynamic functional connectivity density analysis

The analysis of dFCD was processed with the DynamicBC software. dFCD analysis was conducted using the sliding-window approach, with a window size of 50 TR and a sliding step of 10 TR. We obtained a global FCD map in each window by computing Pearson’s correlations between the truncated time course of all pairs of voxels within the automated anatomical labeling-90 (AAL-90) atlas (comprising 45 cortical and subcortical regions in each hemisphere) (Tzourio-Mazoyer et al., 2002), yielding a set of sliding-window FCD maps for each subject. We used r = 0.2 as the correlation coefficient threshold to define the connectivity between two voxels. If their correlation coefficient was larger than 0.2, then connectivity was present between them. The threshold was selected to eliminate the weak correlations induced by noise (Li et al., 2019). Subsequently, the temporal variability was estimated by computing the SD of FCD across sliding windows. In consideration that the global signal regression may induce controversial negative correlations (Fox et al., 2009; Murphy et al., 2009), all our analyses were performed based on positive correlation above a threshold of 0.2.

Brain gene expression data processing

Brain gene expression data were obtained from the AHBA dataset², which was derived from six human postmortem donors. The original expression data of more than 20,000 genes at 3,702 spatially distinct brain tissue samples were processed using a newly proposed pipeline (Arnatkeviciute et al., 2019). The dataset from AHBA was processed following recommendations from previous studies to ensure consistency and reproducibility of results. Partial Least Squares (PLS) regression was employed to investigate the regional dFCD maps in relation to the transcriptional activity of all 10,027 genes.

Gene enrichment analysis

Gene enrichment analysis was performed utilizing the DAVID Functional Annotation Bioinformatics Microarray Analysis platform, accessible at https://metascape. org/gp/index.html#/main/step1. Gene Ontology (GO) was employed to assess biological functions, covering biological processes (BP), cellular components (CC), and molecular functions (MF). Pathway databases, particularly the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, were utilized to evaluate pertinent biological pathways. Enrichment analysis plots were created with SRplot available at http://www.bioinformatics.com.cn/.

Specific expression analysis

Cell type-specific expression analysis, brain region-specific expression analysis, and temporal-specific expression analyses were performed at doughertytools.wustl.edu/CSEAtool.html to identify the distinct cell types, brain regions, and developmental stages showing an overrepresentation of these genes. A specificity index probability (pSI) was utilized with two thresholds (0.01 and 0.001) to evaluate the enrichment of genes in specific terms relative to others.

Protein-protein interaction analysis

Protein-protein interaction (PPI) analysis was conducted using STRING v11.0³ to explore whether genes linked to brain functional alterations could form a PPI network. PPI pairs with a confidence interaction score > 0.9 were extracted for analysis. Hub genes were identified as those within the top 15% of degree values, representing the number of edges connected to a gene. Furthermore, the Human Brain Transcriptome database⁴ was utilized to delineate the spatial-temporal expression trajectory of hub genes exhibiting the highest degree values.

An overview of the analysis pipeline is illustrated in Figure 1.

FIGURE 1

Statistical analysis

The demographic and clinical characteristics of the groups were analyzed using SPSS version 20.0 software (SPSS Inc., Chicago, IL, United States), employing the Chi-square test and independent-samples t-test, with a p-value of less than 0.05 considered indicative of a significant statistical difference. Intragroup patterns of dFCD maps were assessed through one-sample t-tests using SPM8 software.

Furthermore, two-sample t-tests were conducted to investigate group disparities in the dFCD maps (50TRs), employing the Gaussian random field (GRF) method for correcting multiple comparisons. Age and sex were regressed as covariates, and the analysis was performed with SPM8 software at a voxel level of p < 0.001 with GRF correction and a cluster level of p < 0.05.

Additionally, partial least squares (PLS) regression was utilized to explore the relationship between regional changes in dFCD and transcriptional activity for all 10,027 genes.

Results

Demographic and clinical characteristics

Significant differences were observed in best-corrected visual acuity (P < 0.001) between the two groups. However, there were no significant differences in sex, age, or body mass index between the two groups. Further details are provided in Table 1.

Differences in dFCD variability between DR and HC

The spatial distribution of dFCD variability in the DR group and the HC group is shown in Figures 2A,B. The results of the two-sample t-tests indicate that, compared to the HC group, the DR group exhibited significantly decreased dFCD values in the left anterior cingulum, left superior occipital gyrus and right postcentral gyrus (two tailed, voxel level of p < 0.001 with GRF correction and a cluster level of p < 0.05) (Figures 2C,D; Table 2).

FIGURE 2

Gene enrichment results

To elucidate the biological functions and pathways associated with brain functional alterations, the t-values of dFCD maps (Figure 3A) between two groups and gene (Figure 3B) was analyzed. we conducted functional enrichment analyses. For the positive genes related to dFCD, there were a total of 1,318 genes (Figure 3C). Gene Ontology (GO) analysis revealed 7 biological processes (BP), including gamma-aminobutyric acid signaling pathway, synaptic vesical lumen acidification, synaptic transmission GABAergic, inhibitory synapse assembly, regulation of postsynaptic membrane potential, neurological system process, and chloride transmembrane transport. Additionally, 7 cellular components (CC) were identified, such as GABA-A receptor complex, transmembrane transporter complex, chloride channel complex, GABA-ergic synapse, cytosol, membrane, and cytoplasm. Moreover, 7 molecular functions (MF) were annotated, including GABA-gated chloride ion channel activity, inhibitory extracellular ligand-gated ion channel activity, GABA-A receptor activity, extracellular ligand-gated ion channel activity, chloride channel activity, ATPase activity, and protein binding (Figure 3E). In terms of KEGG pathways, 15 pathways were enriched, including Biosynthesis of unsaturated fatty acids, Nicotine addiction, GABAergic synapse and Fatty acid metabolism, etc. (Figure 3F).

FIGURE 3

For the negative genes related to dFCD, there were also 1,318 genes identified (Figure 3D). GO analysis revealed 7 biological processes, such as cell-cell adhesion mediated by cadherin, adherens junction organization, cell-cell junction assembly, SCF-dependent proteasomal ubiquitin-dependent protein catabolic process, Rho protein signal transduction, regulation of translation, and cell morphogenesis. Additionally, 7 cellular components and 7 molecular functions were identified (Figure 3G), along with 15 enriched KEGG pathways, including Glycosaminoglycan biosynthesis-keratan sulfate, Various types of N-glycan biosynthesis, and Glycosphingolipid biosynthesis-lacto and neolacto series, etc. (Figure 3H).

Specific expression analysis

For the 1,318 positive genes associated with dFCD, cell-specific expression analysis revealed their distinct expression across various cell types, notably in Cb.Pcp2 (Figure 4A; Table 3). Moreover, brain-specific expression analysis highlighted their specific expression patterns in different brain regions, particularly in the Cortex and Striatum (Figure 4C; Table 4). Temporal-specific expression analysis further elucidated their differential expression across developmental stages and brain regions, however, no significant differential expression was found in this study (Figure 4E; Table 5).

FIGURE 4

In contrast, cell-specific expression analysis of the 1,318 negative genes linked to dFCD did not reveal significant expression (Figure 4B; Table 3). Brain-specific expression analysis showcased their distinct expression in different brain regions, notably in the Striatum (Figure 4D; Table 4). Temporal-specific expression analysis delineated their varied expression across developmental stages and brain regions, similar to the positive genes, we did not find significant differential expression (Figure 4E; Table 5).

PPI network and hub genes

PPI analysis uncovered that 1,318 positive genes linked to dFCD formed an interlinked PPI network (Figure 5A). Additionally, we elucidated the spatial-temporal expression patterns of three key genes with the highest degree values: RBX1 (Figure 5B), DDB2 (Figure 5C), and CBX2 (Figure 5D).

FIGURE 5

Similarly, PPI analysis revealed that 1,318 negative genes associated with dFCD constructed an interconnected PPI network (Figure 6A). We further outlined the spatial-temporal expression trajectories of three pivotal genes with the highest degree values: RBX1 (Figure 6B), DCAF13 (Figure 6C), and CUL4A (Figure 6D).

FIGURE 6

Discussion

This study represents an innovative approach integrating functional magnetic resonance imaging with transcriptome-neuroimaging spatial correlation to explore the genetic underpinnings of dFCD variability in diabetic retinopathy (DR). Our investigation revealed a significant reduction in dFCD variability within the left anterior cingulum, left superior occipital gyrus, and right postcentral gyrus in the DR group compared to the healthy control (HC) group. Transcriptome-neuroimaging correlation analyses identified two groups of DR-related genes (PLS1 + and PLS1-) that exhibited substantial enrichment in biological functions and pathways related to ion channels, synaptic function, neuronal systems, and cellular junctions. Specific expression analyses indicated the unique expression of these genes in Purkinje neurons, the cortex, and striatum brain regions, with no developmental period-specific expression observed. Additionally, genes linked to altered brain function were capable of forming a protein-protein interaction (PPI) network supported by hub genes. These insights suggest that the disrupted brain function in DR patients arises from intricate interactions among multiple genes with diverse functional characteristics, offering novel perspectives on the molecular mechanisms linked to dFCD alterations in DR patients.

Significantly reduced dFCD variability was observed in the left anterior cingulum, left superior occipital gyrus, and right postcentral gyrus in the diabetic retinopathy (DR) group compared to the healthy control (HC) group. The anterior cingulate gyrus, a component of the salience network, plays a role in higher functions like attention allocation (Pardo et al., 1990), reward anticipation, decision making, and emotion (Jackson et al., 2006). A study by Liao et al showed significantly reduced ReHo values in the right anterior cingulate gyrus in patients with DR and hypothesized that this reflected dysfunction of the anterior cingulate gyrus (Liao X.-L et al., 2019). Similar to their results, we found significantly reduced dFCD variability in the left anterior cingulate gyrus in DR patients, which may imply anterior cingulate gyrus dysfunction and lead to corresponding cognitive impairment in DR patients. The superior occipital gyrus is located in the occipital lobe and is mainly responsible for the processing of visual information. It has been found that functional connectivity from the bilateral lingual gyrus to the left superior occipital gyrus was reduced in the DR group compared to healthy controls (Huang et al., 2021b), which suggests disrupted connectivity of the visual network and the dorsal visual pathway, and may reflect impaired top-down modulation and visuospatial function in patients with DR. In another study, gray matter density (GMD) in the right occipital lobe of DR patients was significantly lower than that of healthy controls (Wessels et al., 2006). In addition, there is a significant decrease in intra-network functional connectivity (FC) of the visual network (VN) in DR patients (Huang et al., 2020). Combining these studies, we suggest that the reduced dFCD variability in the left superior occipital gyrus may reflect impaired visual function in DR patients. The postcentral gyrus is located in sensorimotor network (SMN) (Chenji et al., 2016), which plays an important role in motion control. One study found reduced local clustering in the postcentral gyrus in DR patients, which correlates with motor speed (van Duinkerken et al., 2016). Postcentral gyrus is closely associated with PVA-related (primary visual areas) spontaneous brain activity (Wang et al., 2008), and their connection affects the processing of spatial visual information (Ringach et al., 1996). We hypothesize that reduced dFCD variability in the postcentral gyrus of DR patients may cause abnormalities in motor control, such as abnormal eye movements.

Transcription-neuroimaging spatial association analysis is an effective method to study the genetic mechanisms underlying neuroimaging phenotypes in the patient’s brain (Zhang et al., 2022), and can identify genes associated with neuroimaging phenotypes (Fornito et al., 2019). Using transcriptional-neuroimaging spatial association analysis, we identified two groups of genes, PLS + and PLS-, which correlate with changes in dFCD in DR patients. Enrichment analyses showed that these genes were mainly enriched for biological functions and pathways such as ion channels, synapses and cellular junctions. It was shown that ion channels are closely related to neuronal growth and differentiation, membrane potential generation, signal transduction and neurotransmitter release (Moody and Bosma, 2005; Smith and Walsh, 2020; Kumar et al., 2016; Hull and Isom, 2018). Alterations in the expression levels of somatostatin (SRIF) or specific SRIF receptors affect the morphology of retinal cell types (i.e., optic rod cells) and may influence transmitter release in photoreceptors by modulating voltage-gated K + and Ca2 + currents, making SRIF a potentially promising drug SRIF for the treatment of DR (Casini et al., 2005). 11-cis-retinal intervention improves aberrant outer retinal ion channel closure, thereby alleviate diabetes-induced outer retinal layer dysfunction (Berkowitz et al., 2012). This implies that ion channels may be an important therapeutic target for DR. Inter-neuronal communication is largely dependent on synaptic function (Pereda, 2014), and in hyperglycemia, both neurochemical (synaptic) and electrical (gap junction) modes of communication between retinal cells are affected (Eleftheriou et al., 2020). Diabetes causes a cumulative loss of neural dendrites and synapses in the inner retina, which results in significant thinning of the inner plexiform layer (IPL) (Barber et al., 1998). Hyperphosphorylated tau proteins contribute to diabetic retinal neurodegeneration by disrupting synaptic and mitochondrial function in retinal ganglion cells (RGCs) (Zhu et al., 2018). This suggests the importance of synaptic dysfunction in the pathophysiological process of DR. In DR, the blood-retinal barrier (BRB) is disrupted (Xu et al., 2011). Inflammatory cytokine levels are elevated in the vitreous fluid of patients with DR (Daruich et al., 2018), whereas pro-inflammatory cytokines have been found to cause down-regulation of tight junctions in bovine retinal endothelial cells through the NF-κB pathway, and, increase retinal endothelial cell permeability (Aveleira et al., 2010). In the mouse retina, hyperglycemic memory (HGM) induces sustained oxidative stress, mitochondrial membrane potential collapse and fission, adhesion junction disassembly and subsequent vascular leakage after normalization of blood glucose (Lee et al., 2022). The association between disruption of the blood-retinal barrier (BRB) and the down-regulation or disintegration of cellular connections is evident. Conversely, the overexpression of endothelial CYP2J2 preserves the distribution of endothelial tight junctions and adherens junctions in an ANXA1-dependent manner, thus safeguarding the integrity of the BRB and shielding retinal ganglion cells from loss following ischemia-reperfusion injury (Zhao et al., 2022).

The PLS + genes exhibited specific expression mainly in the cortex and striatum, while the PLS- genes showed predominant expression in the striatum, indicating a focus of diabetic retinopathy (DR) effects on the cerebral cortex. Moreover, the PLS + gene demonstrated particular expression in Purkinje neurons (Cb.Pcp2) located in the cerebellar cortex (Hirano, 2018), pivotal for integrating cerebellar inputs and outputs. The regular firing activity of Purkinje neurons is maintained through precise ion channel control (Raman and Bean, 1999), and ion channel dysregulation leads to perturbations in the intrinsic membrane excitability of Purkinje neurons, causing abnormal firing in Purkinje neurons (Bushart and Shakkottai, 2019; Cook et al., 2021; Huang and Shakkottai, 2023), which corresponds to the enrichment of these genes in ion channels. We did not find significant specific expression of PLS- genes at the cellular level. Furthermore, time-specific expression analyses, using predefined thresholds, did not reveal specific expression of dFCD-related genes at distinct developmental stages.

PPI analyses have revealed that dFCD-related genes can form a PPI network supported by functionally significant hub genes crucial for various biological processes. For instance, the RBX1 gene encodes RBX1 (RING box protein 1), essential for developmental processes. RBX1 suppression triggers the DNA damage response (DDR), resulting in G2-M arrest, senescence, and apoptosis (Jia et al., 2011). Another example is DDB2 (damage-specific DNA-binding protein 2), a component of the group E gene associated with dyschromic dry dermatosis, participating in the initiation of molecular ubiquitination within the ubiquitin ligase complex involving DDB1 and CUL4A (cullin 4A) (Stoyanova et al., 2009; Yan et al., 2011). DDB2 directly interacts with progestin and adipoQ receptor 3 (PAQR3), facilitating PAQR3 ubiquitination, consequently enhancing PI3K/AKT axis signaling to ameliorate insulin resistance and mitigate retinopathy in diabetic mice (Xiao et al., 2022).

Several limitations are evident in this study. Firstly, the gene expression data originated from the brains of six postmortem donors, while the neuroimaging data were obtained from 93 living participants (Zeng et al., 2012; Hawrylycz et al., 2015). To address the differential sources, we focused on genes displaying more conserved expression profiles by setting a DS threshold, facilitating our transcriptome-neuroimaging spatial correlation analysis. Secondly, due to limited gene expression data in the right hemisphere and notable variations between cortical and subcortical regions, only tissue samples from the left cerebral cortex were included, potentially introducing bias. Thirdly, oversampling in spatial correlation analyses is a common concern, which may inflate the rate of false positives. Lastly, gene identification related to dynamic functional connectivity density (dFCD) was achieved indirectly through spatial correlation analysis of transcriptional neuroimaging, lacking direct experimental validation. Further exploration and validation through animal experiments are recommended for future studies.

Conclusion

Our study identified altered dFCD variability in brain regions linked to visual and cognitive functions in DR patients. This was manifested by significantly lower dFCD variability in the left anterior cingulate gyrus, left supraoccipital gyrus, and right postcentral gyrus in DR group. Moreover, transcriptome-neuroimaging correlation analyses revealed a spatial association between these dFCD changes and the genes with unique functional profiles. These genes were enriched in biologically significant functions and pathways, specific to certain cells and brain areas. Our research offers novel understandings of the genetic mechanisms influencing dFCD alterations in DR. The findings are also important for the diagnosis of DR and the discovery of potential therapeutic targets for DR.

References

• 1

  ArnatkeviciuteA.FulcherB. D.FornitoA. (2019). A practical guide to linking brain-wide gene expression and neuroimaging data.Neuroimage189353–367. 10.1016/j.neuroimage.2019.01.011

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 2

  ArnatkeviciuteA.FulcherB. D.BellgroveM. A.FornitoA. (2022). Imaging transcriptomics of brain disorders.Biol. Psychiatry Glob. Open Sci.2319–331. 10.1016/j.bpsgos.2021.10.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3

  AveleiraC. A.LinC.-M.AbcouwerS. F.AmbrósioA. F.AntonettiD. A. (2010). TNF-α signals through PKCζ/NF-κB to alter the tight junction complex and increase retinal endothelial cell permeability.Diabetes592872–2882. 10.2337/db09-1606

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4

  BahtiyarG.GuttermanD.LebovitzH. (2016). Heart failure: A major cardiovascular complication of diabetes mellitus.Curr. Diab. Rep.16:116. 10.1007/s11892-016-0809-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5

  BansalS.BurmanA.TripathiA. K. (2023). Advanced glycation end products: Key mediator and therapeutic target of cardiovascular complications in diabetes.World J. Diabetes141146–1162. 10.4239/wjd.v14.i8.1146

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 6

  BarberA. J.LiethE.KhinS. A.AntonettiD. A.BuchananA. G.GardnerT. W. (1998). Neural apoptosis in the retina during experimental and human diabetes.Early onset and effect of insulin. J. Clin. Invest.102783–791. 10.1172/JCI2425

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7

  BarkhofF.HallerS.RomboutsS. A. R. B. (2014). Resting-state functional MR imaging: A new window to the brain.Radiology27229–49. 10.1148/radiol.14132388

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 8

  BerkowitzB. A.BissigD.PatelP.BhatiaA.RobertsR. (2012). Acute systemic 11-cis-retinal intervention improves abnormal outer retinal ion channel closure in diabetic mice.Mol. Vis.18372–376.

  • Google Scholar

• 9

  BiesselsG. J.DespaF. (2018). Cognitive decline and dementia in diabetes mellitus: Mechanisms and clinical implications.Nat. Rev. Endocrinol.14591–604. 10.1038/s41574-018-0048-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10

  BushartD. D.ShakkottaiV. G. (2019). Ion channel dysfunction in cerebellar ataxia.Neurosci. Lett.68841–48. 10.1016/j.neulet.2018.02.005

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11

  CalhounV. D.MillerR.PearlsonG.AdalıT. (2014). The chronnectome: Time-varying connectivity networks as the next frontier in fMRI data discovery.Neuron84262–274. 10.1016/j.neuron.2014.10.015

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 12

  CasiniG.CatalaniE.Dal MonteM.BagnoliP. (2005). Functional aspects of the somatostatinergic system in the retina and the potential therapeutic role of somatostatin in retinal disease.Histol. Histopathol.20615–632. 10.14670/HH-20.615

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13

  ChaiY.-H.HanY.-P.ZhangJ.-Y.ZhouJ.-B. (2024). Diabetic retinopathy and brain structure, cognition function, and dementia: A bidirectional mendelian randomization study.J. Alzheimers Dis.971211–1221. 10.3233/JAD-231022

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 14

  ChenjiS.JhaS.LeeD.BrownM.SeresP.MahD.et al (2016). Investigating default mode and sensorimotor network connectivity in amyotrophic lateral sclerosis.PLoS One11:e0157443. 10.1371/journal.pone.0157443

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15

  ChoH.SobrinL. (2014). Genetics of diabetic retinopathy.Curr. Diab. Rep.14:515. 10.1007/s11892-014-0515-z

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 16

  CiudinA.EspinosaA.Simó-ServatO.RuizA.AlegretM.HernándezC.et al (2017). Type 2 diabetes is an independent risk factor for dementia conversion in patients with mild cognitive impairment.J. Diabetes Complications311272–1274. 10.1016/j.jdiacomp.2017.04.018

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17

  CookA. A.FieldsE.WattA. J. (2021). Losing the beat: Contribution of purkinje cell firing dysfunction to disease, and its reversal.Neuroscience462247–261. 10.1016/j.neuroscience.2020.06.008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18

  Crosby-NwaobiR. R.SivaprasadS.AmielS.ForbesA. (2013). The relationship between diabetic retinopathy and cognitive impairment.Diabetes Care363177–3186. 10.2337/dc12-2141

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 19

  DaruichA.MatetA.MoulinA.KowalczukL.NicolasM.SellamA.et al (2018). Mechanisms of macular edema: Beyond the surface.Prog. Retin. Eye Res.6320–68. 10.1016/j.preteyeres.2017.10.006

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20

  EdwardsS. L.BeesleyJ.FrenchJ. D.DunningA. M. (2013). Beyond GWASs: Illuminating the dark road from association to function.Am. J. Hum. Genet.93779–797. 10.1016/j.ajhg.2013.10.012

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 21

  EleftheriouC. G.IvanovaE.SagdullaevB. T. (2020). Of neurons and pericytes: The neuro-vascular approach to diabetic retinopathy.Vis. Neurosci.37:E005. 10.1017/S0952523820000048

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22

  FornitoA.ArnatkevičiūtėA.FulcherB. D. (2019). Bridging the gap between connectome and transcriptome.Trends Cogn. Sci.2334–50. 10.1016/j.tics.2018.10.005

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 23

  FornitoA.ZaleskyA.BassettD. S.MeunierD.Ellison-WrightI.YücelM.et al (2011). Genetic influences on cost-efficient organization of human cortical functional networks.J. Neurosci.313261–3270. 10.1523/JNEUROSCI.4858-10.2011

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 24

  FoxM. D.RaichleM. E. (2007). Spontaneous fluctuations in brain activity observed with functional magnetic resonance imaging.Nat. Rev. Neurosci.8700–711. 10.1038/nrn2201

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25

  FoxM. D.ZhangD.SnyderA. Z.RaichleM. E. (2009). The global signal and observed anticorrelated resting state brain networks.J. Neurophysiol.1013270–3283. 10.1152/jn.90777.2008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 26

  GrassiM. A.TikhomirovA.RamalingamS.BelowJ. E.CoxN. J.NicolaeD. L. (2011). Genome-wide meta-analysis for severe diabetic retinopathy.Hum. Mol. Genet.202472–2481. 10.1093/hmg/ddr121

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 27

  HäggS.ThornL. M.PutaalaJ.LiebkindR.HarjutsaloV.ForsblomC. M.et al (2013). Incidence of stroke according to presence of diabetic nephropathy and severe diabetic retinopathy in patients with type 1 diabetes.Diabetes Care364140–4146. 10.2337/dc13-0669

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28

  HawrylyczM.MillerJ. A.MenonV.FengD.DolbeareT.Guillozet-BongaartsA. L.et al (2015). Canonical genetic signatures of the adult human brain.Nat. Neurosci.181832–1844. 10.1038/nn.4171

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29

  HiranoT. (2018). Purkinje neurons: Development, morphology, and function.Cerebellum17699–700. 10.1007/s12311-018-0985-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30

  HuangH.ShakkottaiV. G. (2023). Targeting ion channels and purkinje neuron intrinsic membrane excitability as a therapeutic strategy for cerebellar ataxia.Life (Basel)13:1350. 10.3390/life13061350

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31

  HuangX.TongY.QiC.-X.DanH.-D.DengQ.-Q.ShenY. (2020). Large-scale neuronal network dysfunction in diabetic retinopathy.Neural Plast.2020:6872508. 10.1155/2020/6872508

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32

  HuangX.TongY.QiC.-X.XuY.-T.DanH.-D.ShenY. (2019). Disrupted topological organization of human brain connectome in diabetic retinopathy patients.Neuropsychiatr. Dis. Treat.152487–2502. 10.2147/NDT.S214325

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 33

  HuangX.WenZ.QiC.-X.TongY.ShenY. (2021a). Dynamic changes of amplitude of low-frequency fluctuations in patients with diabetic retinopathy.Front. Neurol.12:611702. 10.3389/fneur.2021.611702

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 34

  HuangX.XieB.-J.QiC.-X.TongY.ShenY. (2021b). Abnormal intrinsic functional network hubs in diabetic retinopathy patients.Neuroreport32498–506. 10.1097/WNR.0000000000001620

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35

  HullJ. M.IsomL. L. (2018). Voltage-gated sodium channel β subunits: The power outside the pore in brain development and disease.Neuropharmacology13243–57. 10.1016/j.neuropharm.2017.09.018

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36

  HutchisonR. M.WomelsdorfT.AllenE. A.BandettiniP. A.CalhounV. D.CorbettaM.et al (2013). Dynamic functional connectivity: Promise, issues, and interpretations.Neuroimage80360–378. 10.1016/j.neuroimage.2013.05.079

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37

  JacksonP. L.BrunetE.MeltzoffA. N.DecetyJ. (2006). Empathy examined through the neural mechanisms involved in imagining how I feel versus how you feel pain.Neuropsychologia44752–761. 10.1016/j.neuropsychologia.2005.07.015

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38

  JiaL.BickelJ. S.WuJ.MorganM. A.LiH.YangJ.et al (2011). RBX1 (RING box protein 1) E3 ubiquitin ligase is required for genomic integrity by modulating DNA replication licensing proteins.J. Biol. Chem.2863379–3386. 10.1074/jbc.M110.188425

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 39

  KimJ.JensenA.KoS.RaghavanS.PhillipsL. S.HungA.et al (2022). Systematic heritability and heritability enrichment analysis for diabetes complications in UK Biobank and ACCORD studies.Diabetes711137–1148. 10.2337/db21-0839

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 40

  KumarP.KumarD.JhaS. K.JhaN. K.AmbastaR. K. (2016). Ion Channels in neurological disorders.Adv. Protein Chem. Struct. Biol.10397–136. 10.1016/bs.apcsb.2015.10.006

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41

  KwaV. I. H.van der SandeJ. J.StamJ.TijmesN.VroolandJ. L.et alAmsterdam Vascular (2002). Retinal arterial changes correlate with cerebral small-vessel disease.Neurology591536–1540. 10.1212/01.wnl.0000033093.16450.5c

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 42

  LeeY.-J.JeonH.-Y.LeeA.-J.KimM.HaK.-S. (2022). Dopamine ameliorates hyperglycemic memory-induced microvascular dysfunction in diabetic retinopathy.FASEB J.36:e22643. 10.1096/fj.202200865R

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43

  LeslieR. D.PykeD. A. (1982). Diabetic retinopathy in identical twins.Diabetes3119–21. 10.2337/diab.31.1.19

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 44

  LiJ.DuanX.CuiQ.ChenH.LiaoW. (2019). More than just statics: Temporal dynamics of intrinsic brain activity predicts the suicidal ideation in depressed patients.Psychol. Med.49852–860. 10.1017/S0033291718001502

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 45

  LiaoW.LiJ.JiG.-J.WuG.-R.LongZ.XuQ.et al (2019). Endless fluctuations: Temporal dynamics of the amplitude of low frequency fluctuations.IEEE Trans. Med. Imaging382523–2532. 10.1109/TMI.2019.2904555

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46

  LiaoX.-L.YuanQ.ShiW.-Q.LiB.SuT.LinQ.et al (2019). Altered brain activity in patients with diabetic retinopathy using regional homogeneity: A resting-state fMRI study.Endocr. Pract.25320–327. 10.4158/EP-2018-0517

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 47

  Lopes de FariaJ. B.SilvaK. C.Lopesde FariaJ. M. (2011). The contribution of hypertension to diabetic nephropathy and retinopathy: The role of inflammation and oxidative stress.Hypertens Res.34413–422. 10.1038/hr.2010.263

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 48

  LvH.WangZ.TongE.WilliamsL. M.ZaharchukG.ZeinehM.et al (2018). Resting-state functional MRI: Everything that nonexperts have always wanted to know.AJNR Am. J. Neuroradiol.391390–1399. 10.3174/ajnr.A5527

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 49

  MoodyW. J.BosmaM. M. (2005). Ion channel development, spontaneous activity, and activity-dependent development in nerve and muscle cells.Physiol. Rev.85883–941. 10.1152/physrev.00017.2004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 50

  MurphyK.BirnR. M.HandwerkerD. A.JonesT. B.BandettiniP. A. (2009). The impact of global signal regression on resting state correlations: Are anti-correlated networks introduced?Neuroimage44893–905. 10.1016/j.neuroimage.2008.09.036

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 51

  PangY.CuiQ.DuanX.ChenH.ZengL.ZhangZ.et al (2017). Extraversion modulates functional connectivity hubs of resting-state brain networks.J. Neuropsychol.11347–361. 10.1111/jnp.12090

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 52

  PardoJ. V.PardoP. J.JanerK. W.RaichleM. E. (1990). The anterior cingulate cortex mediates processing selection in the Stroop attentional conflict paradigm.Proc. Natl. Acad. Sci. U S A.87256–259. 10.1073/pnas.87.1.256

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 53

  PeredaA. E. (2014). Electrical synapses and their functional interactions with chemical synapses.Nat. Rev. Neurosci.15250–263. 10.1038/nrn3708

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 54

  QiC.-X.HuangX.ShenY. (2020). Altered intrinsic brain activities in patients with diabetic retinopathy using amplitude of low-frequency fluctuation: A resting-state fMRI study.Diabetes Metab. Syndr. Obes.132833–2842. 10.2147/DMSO.S259476

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 55

  RamanI. M.BeanB. P. (1999). Ionic currents underlying spontaneous action potentials in isolated cerebellar Purkinje neurons.J. Neurosci.191663–1674. 10.1523/JNEUROSCI.19-05-01663.1999

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 56

  RingachD. L.HawkenM. J.ShapleyR. (1996). Binocular eye movements caused by the perception of three-dimensional structure from motion.Vision Res.361479–1492. 10.1016/0042-6989(95)00285-5

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 57

  SanahujaJ.AlonsoN.DiezJ.OrtegaE.RubinatE.TravesetA.et al (2016). Increased burden of cerebral small vessel disease in patients with type 2 diabetes and retinopathy.Diabetes Care391614–1620. 10.2337/dc15-2671

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 58

  ShiW.-Q.ZhangM.-X.TangL.-Y.YeL.ZhangY.-Q.LinQ.et al (2022). Altered spontaneous brain activity patterns in patients with diabetic retinopathy using amplitude of low-frequency fluctuation.World J. Diabetes1397–109. 10.4239/wjd.v13.i2.97

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 59

  SmithR. S.WalshC. A. (2020). Ion channel functions in early brain development.Trends Neurosci.43103–114. 10.1016/j.tins.2019.12.004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 60

  StoyanovaT.RoyN.KopanjaD.BagchiS.RaychaudhuriP. (2009). DDB2 decides cell fate following DNA damage.Proc. Natl. Acad. Sci. U S A.10610690–10695. 10.1073/pnas.0812254106

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 61

  TamV.PatelN.TurcotteM.BosséY.ParéG.MeyreD. (2019). Benefits and limitations of genome-wide association studies.Nat. Rev. Genet.20467–484. 10.1038/s41576-019-0127-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 62

  TeoZ. L.ThamY.-C.YuM.CheeM. L.RimT. H.CheungN.et al (2021). Global prevalence of diabetic retinopathy and projection of burden through 2045: Systematic review and meta-analysis.Ophthalmology1281580–1591. 10.1016/j.ophtha.2021.04.027

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 63

  TingD. S. W.CheungG. C. M.WongT. Y. (2016). Diabetic retinopathy: Global prevalence, major risk factors, screening practices and public health challenges: A review.Clin. Exp. Ophthalmol.44260–277. 10.1111/ceo.12696

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 64

  TomasiD.VolkowN. D. (2010). Functional connectivity density mapping.Proc. Natl. Acad. Sci. U S A.1079885–9890. 10.1073/pnas.1001414107

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 65

  TomasiD.Shokri-KojoriE.VolkowN. D. (2016). High-resolution functional connectivity density: Hub locations, sensitivity, specificity, reproducibility, and reliability.Cereb. Cortex263249–3259. 10.1093/cercor/bhv171

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 66

  Tzourio-MazoyerN.LandeauB.PapathanassiouD.CrivelloF.EtardO.DelcroixN.et al (2002). Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain.Neuroimage15273–289. 10.1006/nimg.2001.0978

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 67

  van de KreekeJ. A.NguyenH. T.KonijnenbergE.TomassenJ.den BraberA.Ten KateM.et al (2018). Retinal and cerebral microvasculopathy: Relationships and their genetic contributions.Invest. Ophthalmol. Vis. Sci.595025–5031. 10.1167/iovs.18-25341

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 68

  Van DijkK. R. A.SabuncuM. R.BucknerR. L. (2012). The influence of head motion on intrinsic functional connectivity MRI.Neuroimage59431–438. 10.1016/j.neuroimage.2011.07.044

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 69

  van DuinkerkenE.IjzermanR. G.KleinM.MollA. C.SnoekF. J.ScheltensP.et al (2016). Disrupted subject-specific gray matter network properties and cognitive dysfunction in type 1 diabetes patients with and without proliferative retinopathy.Hum. Brain Mapp.371194–1208. 10.1002/hbm.23096

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 70

  WangK.JiangT.YuC.TianL.LiJ.LiuY.et al (2008). Spontaneous activity associated with primary visual cortex: A resting-state FMRI study.Cereb. Cortex18697–704. 10.1093/cercor/bhm105

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 71

  WesselsA. M.SimsekS.RemijnseP. L.VeltmanD. J.BiesselsG. J.BarkhofF.et al (2006). Voxel-based morphometry demonstrates reduced grey matter density on brain MRI in patients with diabetic retinopathy.Diabetologia492474–2480. 10.1007/s00125-006-0283-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 72

  WondmenehT. G.MohammedJ. A. (2024). Prevalence of diabetic retinopathy and its associated risk factors among adults in Ethiopia: A systematic review and meta-analysis.Sci. Rep.14:28266. 10.1038/s41598-024-78596-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 73

  WongT. Y.CheungC. M. G.LarsenM.SharmaS.SimóR. (2016). Diabetic retinopathy.Nat. Rev. Dis. Primers2:16012. 10.1038/nrdp.2016.12

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 74

  XiaoH.SunX.LinZ.YangY.ZhangM.XuZ.et al (2022). Gentiopicroside targets PAQR3 to activate the PI3K/AKT signaling pathway and ameliorate disordered glucose and lipid metabolism.Acta Pharm. Sin. B122887–2904. 10.1016/j.apsb.2021.12.023

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 75

  XuH.-Z.SongZ.FuS.ZhuM.LeY.-Z. (2011). RPE barrier breakdown in diabetic retinopathy: Seeing is believing.J. Ocul. Biol. Dis. Infor.483–92. 10.1007/s12177-011-9068-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 76

  XuR.-S. (2015). Pathogenesis of diabetic cerebral vascular disease complication.World J. Diabetes654–66. 10.4239/wjd.v6.i1.54

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 77

  XueK.LiangS.YangB.ZhuD.XieY.QinW.et al (2022). Local dynamic spontaneous brain activity changes in first-episode, treatment-naïve patients with major depressive disorder and their associated gene expression profiles.Psychol. Med.522052–2061. 10.1017/S0033291720003876

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 78

  YanC.-G.WangX.-D.ZuoX.-N.ZangY.-F. (2016). DPABI: Data processing & analysis for (Resting-State) brain imaging.Neuroinformatics14339–351. 10.1007/s12021-016-9299-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 79

  YanY.ZhangX.LegerskiR. J. (2011). Artemis interacts with the Cul4A-DDB1DDB2 ubiquitin E3 ligase and regulates degradation of the CDK inhibitor p27.Cell Cycle104098–4109. 10.4161/cc.10.23.18227

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 80

  ZengH.ShenE. H.HohmannJ. G.OhS. W.BernardA.RoyallJ. J.et al (2012). Large-scale cellular-resolution gene profiling in human neocortex reveals species-specific molecular signatures.Cell149483–496. 10.1016/j.cell.2012.02.052

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 81

  ZhangC.CaiH.XuX.LiQ.LiX.ZhaoW.et al (2022). Genetic architecture underlying differential resting-state functional connectivity of subregions within the human visual cortex.Cereb. Cortex322063–2078. 10.1093/cercor/bhab335

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 82

  ZhangX.GaoY.ZhouZ.WangJ.ZhouQ.LiQ. (2010). Familial clustering of diabetic retinopathy in Chongqing, China, type 2 diabetic patients.Eur. J. Ophthalmol.20911–918. 10.1177/112067211002000516

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 83

  ZhangZ.LiaoM.YaoZ.HuB.XieY.ZhengW.et al (2017). Frequency-specific functional connectivity density as an effective biomarker for adolescent generalized anxiety disorder.Front. Hum. Neurosci.11:549. 10.3389/fnhum.2017.00549

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 84

  ZhaoB.HuangJ.LouX.YaoK.YeM.MouQ.et al (2022). Endothelial CYP2J2 overexpression restores the BRB via METTL3-mediated ANXA1 upregulation.FASEB J.36:e22619. 10.1096/fj.202201061RR

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 85

  ZhuH.ZhangW.ZhaoY.ShuX.WangW.WangD.et al (2018). GSK3β-mediated tau hyperphosphorylation triggers diabetic retinal neurodegeneration by disrupting synaptic and mitochondrial functions.Mol. Neurodegener.13:62. 10.1186/s13024-018-0295-z

  • Pubmed Abstract
  • CrossRef
  • Google Scholar