Saivageethi Nuthikattu
email: snuthikattu@ucdavis.edu

Dragan Milenkovic
email: dragan.milenkovic@inra.fr

Jennifer E Norman
email: jenorman@ucdavis.edu

John Rutledge
email: jcrutledge@ucdavis.edu

Amparo Villablanca
email: avillablanca@ucdavis.edu

Inhibition of Soluble Epoxide Hydrolase Is Protective against the Multiomic Effects of a High Glycemic Diet on Brain Microvascular Inflammation and Cognitive Dysfunction

Keywords: multi-omics, microvascular, brain, dementia, high glycemic diet, soluble epoxide hydrolase inhibitor, maless

Diet is a modifiable risk factor for cardiovascular disease (CVD) and dementia, yet relatively little is known about the effect of a high glycemic diet (HGD) on the brain's microvasculature. The objective of our study was to determine the molecular effects of an HGD on hippocampal microvessels and cognitive function and determine if a soluble epoxide hydrolase (sEH) inhibitor (sEHI), known to be vasculoprotective and anti-inflammatory, modulates these effects. Wild type male mice were fed a low glycemic diet (LGD, 12% sucrose/weight) or an HGD (34% sucrose/weight) with/without the sEHI, trans-4-[4-(3-adamantan-1-yl-ureido)-cyclohexyloxy]-benzoic acid (t-AUCB), for 12 weeks. Brain hippocampal microvascular gene expression was assessed by microarray and data analyzed using a multi-omic approach for differential expression of protein and non-proteincoding genes, gene networks, functional pathways, and transcription factors. Global hippocampal microvascular gene expression was fundamentally different for mice fed the HGD vs. the LGD. The HGD response was characterized by differential expression of 608 genes involved in cell signaling, neurodegeneration, metabolism, and cell adhesion/inflammation/oxidation effects reversible by t-AUCB and hence sEH inhibitor correlated with protection against Alzheimer's dementia. Ours is the first study to demonstrate that high dietary glycemia contributes to brain hippocampal microvascular inflammation through sEH.

Introduction

Dementias are the seventh leading cause of death globally and contribute significantly to health care costs [START_REF] Duplantier | A Critical Review of the Study of Neuroprotective Diets to Reduce Cognitive Decline[END_REF]. Several studies suggest that a high-fat diet or Western diet (high fat and high glycemic content) can lead to reduced cognitive function [START_REF] Wie ¸ckowska-Gacek | Western Diet Induces Impairment of Liver-Brain Axis Accelerating Neuroinflammation and Amyloid Pathology in Alzheimer's Disease[END_REF][START_REF] Leigh | Inflammation and the gut microbiome: Mechanisms for obesity-associated cognitive impairment[END_REF][START_REF] Buie | Obesity-related cognitive impairment: The role of endothelial dysfunction[END_REF][START_REF] Rutkowsky | Reduced cognitive function, increased blood-brain-barrier transport and inflammatory responses, and altered brain metabolites in LDLr -/-and C57BL/6 mice fed a western diet[END_REF][START_REF] Nguyen | Obesity and cognitive decline: Role of inflammation and vascular changes[END_REF]. It has also become increasingly recognized that vasculature plays an important role in the development of dementias [START_REF] Gorelick | Vascular cognitive impairment and dementia[END_REF]. Our group has previously demonstrated the multi-omic and lipotoxic effect of a Western diet on the brain microvasculature and its negative consequences on cognitive function in male and female mice [START_REF] Aung | Lipotoxic brain microvascular injury is mediated by activating transcription factor 3-dependent inflammatory and oxidative stress pathways[END_REF][START_REF] Nuthikattu | The Western Diet Regulates Hippocampal Microvascular Gene Expression: An Integrated Genomic Analyses in Female Mice[END_REF][START_REF] Nuthikattu | Lipotoxic Injury Differentially Regulates Brain Microvascular Gene Expression in Male Mice[END_REF][START_REF] Nuthikattu | Sex-Dependent Molecular Mechanisms of Lipotoxic Injury in Brain Microvasculature: Implications for Dementia[END_REF]. While there have been numerous mechanistic studies focusing on cognitive function with high-fat diets, few studies have explored the impact of high glycemia in the absence of high levels of dietary fat.

There is compelling epidemiological data to suggest that the effects of a high glycemic diet (HGD) on cognitive function and the brain are an important area of study. High blood glucose and high dietary glycemic load were both found to be related to poorer performance in perceptual speed and spatial ability [START_REF] Seetharaman | Blood glucose, Diet-based glycemic load and cognitive aging among dementia-free older adults[END_REF]. An HGD has also been associated with a greater cerebral amyloid burden [START_REF] Taylor | A high-glycemic diet is associated with cerebral amyloid burden in cognitively normal older adults[END_REF]. Further, consumption of a high glycemic index afternoon snack was associated with cognitive decline in apolipoprotein E4 allele carriers [START_REF] Gentreau | High Glycemic Load Is Associated with Cognitive Decline in Apolipoprotein E ε4 Allele Carriers[END_REF]. However, there is relatively sparse animal data on the effects of an HGD on the brain and cognitive function. In rats, an HGD had detrimental effects on memory [START_REF] Saikrishna | Combined Administration of Monosodium Glutamate and High Sucrose Diet Accelerates the Induction of Type 2 Diabetes, Vascular Dysfunction, and Memory Impairment in Rats[END_REF][START_REF] Wong | High-sucrose diets in male rats disrupt aspects of decision making tasks, motivation and spatial memory, but not impulsivity measured by operant delay-discounting[END_REF], disrupted hypothalamic redox homeostasis [START_REF] Żebrowska | High-Sugar Diet Disrupts Hypothalamic but Not Cerebral Cortex Redox Homeostasis[END_REF], and increased hippocampal endoplasmic reticulum stress [START_REF] Pinto | Early and sustained exposure to high-sucrose diet triggers hippocampal ER stress in young rats[END_REF]. In a mouse model of Alzheimer's disease, an HGD was found to increase neuroinflammation and cortical levels of Amyloid-β [START_REF] Yeh | A high-sucrose diet aggravates Alzheimer's disease pathology, attenuates hypothalamic leptin signaling, and impairs food-anticipatory activity in APPswe/PS1dE9 mice[END_REF]. These studies suggest that an HGD is associated with cognitive impairment in animal models.

The impact of high glycemia on cognition may be reversible. One study in rats found that a high sucrose diet increased brain unesterified arachidonic acid and the activity of enzymes facilitating the release of arachidonic acid from phospholipids [START_REF] Taha | Upregulated expression of brain enzymatic markers of arachidonic and docosahexaenoic acid metabolism in a rat model of the metabolic syndrome[END_REF]. Arachidonic acid is metabolized by cytochrome P450 enzymes to epoxyeicosatrienoic acid isomers (EETs), which are short-lasting and locally active neuroprotective, vasodilatory, and antiinflammatory [START_REF] Wang | Neuroprotective effects of epoxyeicosatrienoic acids[END_REF] signaling molecules. Soluble epoxide hydrolase (sEH) is an enzyme that converts EETs into dihydroxyeicosatrienoic acids (DHETs), which have less biological activity [START_REF] Wang | Neuroprotective effects of epoxyeicosatrienoic acids[END_REF][START_REF] Zarriello | Humble beginnings with big goals: Small molecule soluble epoxide hydrolase inhibitors for treating CNS disorders[END_REF]. Inhibiting sEH activity thus increases the amount of beneficial EETs [START_REF] Zarriello | Humble beginnings with big goals: Small molecule soluble epoxide hydrolase inhibitors for treating CNS disorders[END_REF]. Studies have implicated soluble epoxide hydrolase (sEH) in many disorders of the central nervous system, including Parkinson's disease, white matter hyperintensities, vascular cognitive decline/impairment, and Alzheimer's disease [START_REF] Zarriello | Humble beginnings with big goals: Small molecule soluble epoxide hydrolase inhibitors for treating CNS disorders[END_REF][START_REF] Ren | Soluble epoxide hydrolase plays a key role in the pathogenesis of Parkinson's disease[END_REF][START_REF] Nelson | Role of soluble epoxide hydrolase in age-related vascular cognitive decline[END_REF][START_REF] Yu | Soluble Epoxide Hydrolase-Derived Linoleic Acid Oxylipins in Serum Are Associated with Periventricular White Matter Hyperintensities and Vascular Cognitive Impairment[END_REF][START_REF] Griñán-Ferré | Pharmacological Inhibition of Soluble Epoxide Hydrolase as a New Therapy for Alzheimer's Disease[END_REF]. Inhibitors of sEH (sEHI) have been shown to be protective in animal models of stroke [START_REF] Simpkins | Soluble epoxide inhibition is protective against cerebral ischemia via vascular and neural protection[END_REF][START_REF] Shaik | Soluble epoxide hydrolase inhibitor trans-4-[4-(3-adamantan-1-yl-ureido)-cyclohexyloxy]-benzoic acid is neuroprotective in rat model of ischemic stroke[END_REF]. Inhibition of sEH is of great clinical interest as it has also been shown to reduce neuroinflammation and cognitive impairment in animal models of cerebral hypoperfusion and type 1 and type 2 diabetes mellitus [START_REF] Chen | Soluble epoxide hydrolase inhibition Promotes White Matter Integrity and Long-Term Functional Recovery after chronic hypoperfusion in mice[END_REF][START_REF] Hao | Soluble epoxide hydrolase inhibition alleviated cognitive impairments via NRG1/ErbB4 signaling after chronic cerebral hypoperfusion induced by bilateral carotid artery stenosis in mice[END_REF][START_REF] Matin | Soluble epoxide hydrolase inhibition improves cognitive function and parenchymal artery dilation in a hypertensive model of chronic cerebral hypoperfusion[END_REF][START_REF] Minaz | An inhibitor of soluble epoxide hydrolase ameliorates diabetes-induced learning and memory impairment in rats[END_REF][START_REF] Wu | Inhibition of soluble epoxide hydrolase (sEH) protects hippocampal neurons and reduces cognitive decline in type 2 diabetic mice[END_REF]. Further, the use of sEHI and genetic knockout of the sEH gene reduces cognitive impairment in animal models of age-related cognitive decline and Alzheimer's disease [START_REF] Griñán-Ferré | Pharmacological Inhibition of Soluble Epoxide Hydrolase as a New Therapy for Alzheimer's Disease[END_REF][START_REF] Ghosh | An epoxide hydrolase inhibitor reduces neuroinflammation in a mouse model of Alzheimer's disease[END_REF][START_REF] Lee | Genetic deletion of soluble epoxide hydrolase delays the progression of Alzheimer's disease[END_REF].

The hippocampus is central to the formation of memory [START_REF] Knierim | The hippocampus[END_REF], and dysfunction of the microvasculature can contribute to the development of dementia [START_REF] Han | Cerebral microvascular dysfunction and neurodegeneration in dementia[END_REF]. To our knowledge, no studies have been published to date examining the effects of an HGD on brain hippocampal microvascular gene expression. The objectives of this study were to use a male murine model to comprehensively characterize the effect of an HGD on neurovascular function through hippocampal microvascular multi-omics and to assess the impact of an HGD on cognitive function. We hypothesized that an HGD would result in injurious differential gene expression changes characterized by brain hippocampal microvascular oxidation, inflammation, and blood-brain barrier disruption. Further, we aimed to determine whether the deleterious genomic effects of an HGD could be mitigated by inhibiting sEH.

Materials and Methods

Experimental Animals and Soluble Epoxide Hydrolase Inhibitor (sEHI) Treatment

19-week-old C57BL/6J wild type (WT) male mice (Jackson Laboratories, stock 000664) were fed a standard chow diet (catalog no 0915 from Envigo Teklad Diets, Madison, WI, USA) and allowed to acclimate for one week prior to beginning the study procedures. Mice were then fed for 12 weeks either a Low Glycemic Index diet (LGD, catalog no TD.08485 from Envigo Teklad Diets, Madison, WI, USA) composed of 13% fat, 19.1% protein, 67.9% carbohydrate, as percent kcal, containing 12% sucrose by weight, or a High Glycemic Index diet (HGD, catalog no. TD.05230, Envigo Teklad Diets, Madison, WI, USA) composed of 12.6% fat, 18.7% protein, 68.7% carbohydrate, as percent kcal, containing 34% sucrose by weight. Mice receiving each diet were given 10mg/L of soluble epoxide hydrolase inhibitor (sEHI), trans-4-[4-(3-adamantan-1-yl-ureido)-cyclohexyloxy]-benzoic acid (t-AUCB) (Cayman Chemical, Ann Arbor, MI, USA) containing 1% v/v polyethylene glycol 400 (PEG400) (Millipore, Burlington, MA, USA) in the drinking water for 12 weeks at which point mice were 32 weeks of age and sacrificed. Mice consumed approximately 7 to 7.5 mL of water each day, consistent with previously published work [START_REF] Gui | Soluble epoxide hydrolase inhibitors, t-AUCB, regulated microRNA-1 and its target genes in myocardial infarction mice[END_REF], and 2.5 to 3 mg of t-AUCB (sEHI) per kg per day. There were a total of four experimental treatment groups (n = 7 mice/gp): LGD alone, LGD with sEHI (LGD+sEHI), HGD alone, and HGD with sEHI (HGD+sEHI). Mice were randomly assigned to the dietary groups.

Animals were housed one mouse per cage in a temperature-and humidity-controlled environment with a 12 h light/dark cycle in the University of California, Davis Mouse Biology Program. Body weight was measured at baseline and at the completion of the dietary intervention period, and activity, water, and food intake were monitored daily by vivarium staff. The research was conducted in conformity with the Public Health Service Policy on Humane Care and Use of Laboratory Animals and all protocols approved by the Institutional Animal Care and Use Committee of the University of California, Davis.

Serum Lipid, Glucose, and Insulin Assays

Mice were fasted overnight for 8 hours, and blood obtained by submandibular nick blood draw for the pre diet samples, and by ventricular puncture at the time of sacrifice following completion of the dietary feeding period for the post diet samples. Blood samples were stored at -80 • C. Lipid, glucose and insulin levels were measured in fasted serum samples. Total cholesterol (TC), high-density lipoprotein cholesterol (HDL), and lowdensity lipoprotein cholesterol (LDL) were measured using enzymatic assays from Fisher Diagnostics (Middleton, VA, USA), and precipitation separation from AbCam (Cambridge, MA, USA) adapted to a microplate format. Glucose was measured using enzymatic assays from Fisher Diagnostics (Middleton, VA, USA), and insulin was determined by electrochemiluminescence from Meso Scale Discovery (Rockville, MD, USA) according to the manufacturer's instructions. All serum assays were performed by the UC Davis Mouse Metabolic Phenotyping Center (MMPC) on non-pooled serum samples.

Isolation and Cryosection of Murine Brain Hippocampus

Following completion of the 12 week dietary feeding period, mice were anesthetized by intraperitoneal xylazine/ketamine and euthanized by exsanguination during the light phase of their light/dark cycle. Intact brains were rapidly removed under RNAse free conditions, cut into regions including the temporal lobe segment containing the hippocampus, and embedded using HistoPrep Frozen Tissue Embedding Media (Fisher Scientific, Pittsburgh, PA, USA). To identify the hippocampus and hippocampal neurons, brain sections in the medial aspect of the temporal lobe were stained with hematoxylin and visualized with microscopy as previously described [START_REF] Aung | Lipotoxic brain microvascular injury is mediated by activating transcription factor 3-dependent inflammatory and oxidative stress pathways[END_REF]. The hippocampus was then coronally cryosectioned (8 µm, Leica Frigocut 2800n Cryostat, Leica Biosystems, Buffalo Grove, IL, USA) and placed on charged RNA-free PEN Membrane Glass slides, treated with RNAlater®-ICE (Life Technologies, Grand Island, NY, USA) to prevent RNA degradation, and stored at -80 • C until use. When ready for use, cryosections from the hippocampal segments were submerged in nuclease-free water and dehydrated in desiccant.

Laser Capture Microdissection (LCM) of Hippocampal Microvessels

For analysis of gene transcriptome of hippocampal brain microvessels, endothelial microvessels (<20um) were first identified in the hippocampal brain cryosections by alkaline phosphatase staining utilizing 5-bromo-4-chloro-3-indolyl phosphate/nitro blue tetrazolium chloride (BCIP/NBT) substrate as previously described [START_REF] Ball | Isolating vessels from the mouse brain for gene expression analysis using laser capture microdissection[END_REF]. Laser capture microdissection (LCM) was then used to isolate the microvascular endothelium in hippocampal cryosections by capture of the entire vessel wall under direct microscopic visualization using a Leica LMD6000 Laser Microdissection Microscope (Leica Microsystems, Wetzlar, Germany). Microvessels were not categorized by hippocampal region or subregion, although they primarily corresponded to endothelial enriched sections in hippocampus dorsal segments that would have included CA1 and CA3 regions.

RNA Extraction from Laser Captured Brain Microvessels

Total RNA was extracted from the laser-captured hippocampal brain microvessels (300 per mice, 3 mice per experimental group) using an Arcturus PicoPure™ RNA Isolation Kit (Thermo Fisher Scientific, Santa Clara, CA, USA) according to the manufacturer's instructions. The quality of the RNA from the LCM-derived vessels was assessed by Nanodrop. RNA quantification was performed according to Affymetrix RNA quantification kit with SYBR Green I and ROX™ Passive Reference Dye protocol (Affymetrix, Santa Clara, CA, USA).

Microarray Hybridization and Transcriptome Analysis

For transcriptomics analysis, we used Clariom D Mouse Array (one array per mouse), containing more than 7 million probes for protein-coding and protein non-coding genes such as micro RNAs (miRNAs), small nucleolar RNAs (snoRNAs), and long non-coding RNAs (LncRNAs) (Thermo Fisher, Santa Clara, CA, USA). RNA (122.3 pg) was used to prepare cRNA and sscDNA using GeneChip®WT Pico Kit (Thermo Fisher, Santa Clara, CA, USA). SscDNA (5.5 µg) was fragmented by uracil-DNA glycosylase (UDG) and apurinic/apyrimidinic endonuclease 1 (APE 1) and labeled by terminal deoxynucleotidyl transferase (TdT) using the DNA Labeling Reagent that is covalently linked to biotin. Fragmented and labeled sscDNA samples were then submitted to the UC Davis Genome Center shared resource core for hybridization, staining, and scanning using Thermo Fisher Scientific WT array hybridization protocol following the manufacturer's protocol. Hybridization of fragmented and labeled sscDNA samples was performed using GeneChip™Hybridization Oven 645, and samples were then washed and stained using GeneChip™ Fluidics Station 450. The arrays were scanned using GeneChip™ Scanner 3000 7G (Thermo Fisher Scientific, Santa Clara, CA, USA). Quality control of the microarrays and data analysis was performed using Thermo Fisher Scientific Transcriptome Analysis Console software version 4.0.2. We have deposited the microarray data in GEO, and the accession number is GSE185057.

Bioinformatic Analysis

Bioinformatics analysis of differentially expressed genes was performed by two of the study investigators (SN and DM) using multiple software tools. We compared the following study groups: (A) HGD to LGD (to determine the effect of the HGD diet) and (B) HGD+sEHI to HGD (to determine the effect of the inhibitor on the HGD diet) as shown in the flow chart in Figure 1.

The Principal Component Analysis (PCA) plot of identified differentially expressed genes (DEG) was obtained through ClustVis [START_REF] Metsalu | ClustVis: A web tool for visualizing clustering of multivariate data using Principal Component Analysis and heatmap[END_REF]. miRNA targets of DEG were identified using Mienturnet [START_REF] Licursi | MIENTURNET: An interactive web tool for microRNA-target enrichment and networkbased analysis[END_REF]. LncRNAs of DEG were identified using LncRRIsearch [START_REF] Fukunaga | LncRRIsearch: A Web Server for lncRNA-RNA Interaction Prediction Integrated With Tissue-Specific Expression and Subcellular Localization Data[END_REF] and Rtools CBRC [43]. Canonical pathway analysis was conducted using GeneTrial2 online database [44,[START_REF] Stockel | Multi-omics enrichment analysis using the GeneTrail2 web service[END_REF]. Networks were constructed and visualized using Cytoscape software (version 3.7.1) [46][START_REF]Cytoscape Network Analyzer Application[END_REF][START_REF] Shannon | Cytoscape: A software environment for integrated models of biomolecular interaction networks[END_REF][START_REF] Su | Biological network exploration with Cytoscape 3[END_REF]. Data preparation was performed with the use of several R packages including splitstackshape [START_REF]Splitstackshape[END_REF], data.table [START_REF] Rdata | [END_REF], dplyr [START_REF] Tidyverse | [END_REF][START_REF] Tidyverse | [END_REF] and string [START_REF] Tidyverse | [END_REF][START_REF]String Tidyverse Github[END_REF]. Pathway networks were built for pathways enriched from a global pathway analysis, considering all omic layers components together.

Transcription factor analyses were performed using Enrichr [START_REF] Chen | Enrichr: Interactive and collaborative HTML5 gene list enrichment analysis tool[END_REF][START_REF] Kuleshov | Enrichr: A comprehensive gene set enrichment analysis web server 2016 update[END_REF][START_REF] Xie | Gene Set Knowledge Discovery with Enrichr[END_REF]. Hierarchical clustering and heat map representations of differentially expressed genes (DEG) were performed using PermutMatrix software [59,[START_REF] Caraux | PermutMatrix: A graphical environment to arrange gene expression profiles in optimal linear order[END_REF].

We performed Pearson's correlation analysis between the genes differentially expressed by the HGD vs. LGD and the HGD+sEHI vs. HGD. We used the ggpubr package in R [START_REF]ggpubr R package[END_REF] to obtain the correlation coefficient and the significance level, as well as the scatter plot with regression line and confidence interval. We then performed prediction of neurodegenerative disease trait by performing correlation analysis between the obtained changes in the expression of genes following the HGD vs. LGD or the HGD+sEHI vs. HGD and the gene expression profiles observed in patients with Alzheimer's disease or cognitive disorders. The human genome data were extracted from the gene expression profile database, which was deposited in the GEO (gene expression omnibus) database. The identification of differentially expressed genes between patients and healthy volunteers was performed using GEO2R [START_REF][END_REF], an NCBI web tool that allows comparisons between two or more groups of samples in the GEO series to identify differentially expressed genes across the experimental conditions. The differentially expressed genes were screened according to p-values < 0.05. Pearson's correlation analysis between the genes identified as differentially expressed, following the HGD vs. LGD or the HGD + sEHI vs. HGD, and patients with neurogenerative diseases was performed using ggpubr package in R [START_REF]ggpubr R package[END_REF]. The Principal Component Analysis (PCA) plot of identified differentially expressed genes (DEG) was obtained through ClustVis [START_REF] Metsalu | ClustVis: A web tool for visualizing clustering of multivariate data using Principal Component Analysis and heatmap[END_REF]. miRNA targets of DEG were identified using Mienturnet [START_REF] Licursi | MIENTURNET: An interactive web tool for microRNA-target enrichment and networkbased analysis[END_REF]. LncRNAs of DEG were identified using LncRRIsearch [START_REF] Fukunaga | LncRRIsearch: A Web Server for lncRNA-RNA Interaction Prediction Integrated With Tissue-Specific Expression and Subcellular Localization Data[END_REF] and Rtools CBRC [43]. Canonical pathway analysis was conducted using GeneTrial2 online database [44,[START_REF] Stockel | Multi-omics enrichment analysis using the GeneTrail2 web service[END_REF]. Networks were constructed and visualized using Cytoscape software (version 3.7.1) [46][START_REF]Cytoscape Network Analyzer Application[END_REF][START_REF] Shannon | Cytoscape: A software environment for integrated models of biomolecular interaction networks[END_REF][START_REF] Su | Biological network exploration with Cytoscape 3[END_REF]. Data preparation was performed with the use of several R packages including splitstackshape [START_REF]Splitstackshape[END_REF], data.table [START_REF] Rdata | [END_REF], dplyr [START_REF] Tidyverse | [END_REF][START_REF] Tidyverse | [END_REF] and string [START_REF] Tidyverse | [END_REF][START_REF]String Tidyverse Github[END_REF]. Pathway networks were built for pathways enriched from a global pathway analysis, considering all omic layers components together.

Transcription factor analyses were performed using Enrichr [START_REF] Chen | Enrichr: Interactive and collaborative HTML5 gene list enrichment analysis tool[END_REF][START_REF] Kuleshov | Enrichr: A comprehensive gene set enrichment analysis web server 2016 update[END_REF][START_REF] Xie | Gene Set Knowledge Discovery with Enrichr[END_REF]. Hierarchical clustering and heat map representations of differentially expressed genes (DEG) were performed using PermutMatrix software [59,[START_REF] Caraux | PermutMatrix: A graphical environment to arrange gene expression profiles in optimal linear order[END_REF].

We performed Pearson's correlation analysis between the genes differentially expressed by the HGD vs. LGD and the HGD+sEHI vs. HGD. We used the ggpubr package in R [START_REF]ggpubr R package[END_REF] to obtain the correlation coefficient and the significance level, as well as the scatter plot with regression line and confidence interval. We then performed prediction of neurodegenerative disease trait by performing correlation analysis between the ob- Interaction of diet and inhibitor effects was performed using Thermo Fisher Scientific Transcriptome Analysis Console software version 4.0.2. Interaction was defined by a ≥2 or ≤-2 delta fold change and a p-value < 0.05 when comparing HGD+sEHI vs. LGD+sEHI to HGD vs. LGD.

Statistical Methods

For microarray, ANOVA ebayes (Thermo Fisher Scientific Transcriptome Analysis Console software, Santa Clara, CA, USA) were used for statistical analysis of microvessel transcriptomes. All genes from the microarray with p < 0.05 and ±2.0-fold change were considered as differentially expressed. Mean body weight and plasma lipid levels were expressed as means ± standard error of the mean (SEM), and significance was determined at p ≤ 0.05 using unpaired Student's t-tests (GraphPad software, La Jolla, CA, USA).

Results

The dietary treatment resulted in the expected weight gain in the study mice as follows:

The mean weight for male mice at 20 weeks of age, on the chow diet, prior to initiation of the study diets, was 30 g and increased significantly (p < 0.05) after 12 weeks in both diet groups (LGD mean 36 g, HGD mean 34g). The soluble epoxide hydrolase inhibitor (sEHI) had no effect on body weight.

Total cholesterol (TC) levels at the end of the feeding period increased significantly (p < 0.05) in all the groups when compared to the baseline measurement at 20 weeks of age and did not statistically differ between the LGD (154.9 mg/dL) and HGD (146.2 mg/dL). The sEHI had no significant effect on total cholesterol levels.

Glucose and insulin levels also increased significantly (p < 0.05) in both the LGD (416.1 mg/dL) and HGD (374.7 mg/dL) treatment groups but did not statistically differ between them. The sEHI had no significant effect on glucose and insulin levels.

. Global Gene Expression and Hierarchical Clustering

To define the molecular mechanisms in brain hippocampal microvessels in response to the HGD, we began by assessing global gene expression using principal component analysis (PCA), a genetic distance visualization tool that shows relatedness between populations. PCA plot analysis showed that the global gene expression profiles of mice on the HDG and LGD were distinctly different from each other (Figure 2A). Using loading plot analysis, we further defined that the genes relevant to the separation of the PCA had opposite expression patterns with the HGD vs the LGD diet. Plots for a few genes (Higd2a, Gm24400, Snora30, Cox5b, Mir5125, Acta2, and Map3k7cl) that contributed to separation of the two dietary groups are provided as examples in Figure 2B.

We then performed hierarchical clustering of global gene expression profiles. Hierarchical clustering groups similar data points together and then organizes the clusters into a hierarchy. Using this strategy, we further confirmed that genes with higher levels of expression with the HGD had lower levels of expression with the LGD, and vice versa (Figure 3). Thus, the effect of the HGD on global hippocampal microvascular gene expression was fundamentally opposite compared to the LGD.

Differential Gene Expression

To study the effect of the HGD on gene expression, we compared the HGD to the LGD (flow chart Figure 1A). Statistical analysis of microarray data revealed 608 differentially expressed genes (DEGs) in hippocampal microvessels following the HGD, with the majority of the DEGs being up-regulated by the HGD (468 genes up-regulated vs. 140 genes downregulated) when compared to the LGD, Supplemental Figure S1. The fold-change varied from 2.0 to 25.94 for up-regulated genes and from -16.65 to -2.01 for down-regulated genes (see Supplemental Table S1 for a complete listing of the DEGs). To our knowledge, we also show for the first time that the HGD regulates the expression of protein-coding genes (221) and non-coding genes [START_REF] Di Francesco | Global changes in DNA methylation in Alzheimer's disease peripheral blood mononuclear cells[END_REF] in male murine hippocampal microvessels (Supplemental Figure S2), including 32 long non-coding RNAs (lncRNAs), 25 microRNAs (miRNAs), and 21 small nucleolar RNAs (snoRNAs). The remaining 309 DEGs were other genes (pseudogenes, ribosomal RNAs (rRNAs), unassigned genes, and multiple-complex genes). Principal Component Analysis (PCA) and loading plot of genes expressed in hippocampal microvessels of the high glycemic diet (HGD) and the low glycemic diet (LGD). A) PCA scatter plot of the microarray data shows the trends of the expression profiles of the hippocampal microvasculature in the high glycemic diet (HGD, red circles) and low glycemic diet (LGD, green circles), respectively. The PCA plot captures the variance in a dataset in terms of principal components and displays the most significant of these on the x and y axes. The percentages of the total variation that are accounted for by the 1st and 2nd principal components are shown on the x-and y-axes labels. The data are shown for three biological replicates for each dietary group. B) Loading plot of genes (violet circles) relevant to the separation of PCA. Blue arrows show box plots with expression levels of a few genes (Higd2a, Gm24400, Snora30, Cox5b, Mir5125, Acta2, Map3k7cl, Snora30 and Map3k7cl) that contribute to separation of the two dietary groups (HGD, LGD). In the box plots, the black dots represent gene expression levels, the notch indicates 95% confidence interval around the median of each group, and yellow diamond indicates the average gene expression of each group. We then performed hierarchical clustering of global gene expression profiles. Hierarchical clustering groups similar data points together and then organizes the clusters into a hierarchy. Using this strategy, we further confirmed that genes with higher levels of expression with the HGD had lower levels of expression with the LGD, and vice versa (Figure 3). Thus, the effect of the HGD on global hippocampal microvascular gene expression was fundamentally opposite compared to the LGD.

Differential gene expression.

To study the effect of the HGD on gene expression, we compared the HGD to the LGD (flow chart Figure 1, Section A). Statistical analysis of microarray data revealed 608 differentially expressed genes (DEGs) in hippocampal microvessels following the HGD, with the majority of the DEGs being up-regulated by the HGD (468 genes up-regulated vs. 140 genes down-regulated) when compared to the LGD, Supplemental Figure S1. The fold-change varied from 2.0 to 25.94 for up-regulated genes and from -16.65 to -2.01 for down-regulated genes (see Supplemental Table S1 for a complete listing of the DEGs). To our knowledge, we also show for the first time that the HGD regulates the expression of protein-coding genes (221) and non-coding genes [START_REF] Di Francesco | Global changes in DNA methylation in Alzheimer's disease peripheral blood mononuclear cells[END_REF] in male murine hippocampal microvessels (Supplemental Figure S2), including 32 long non-coding RNAs (lncRNAs), 25 microRNAs (miRNAs), and 21 small nucleolar RNAs (snoRNAs). The remaining 309 DEGs were other genes (pseudogenes, ribosomal RNAs (rRNAs), unassigned genes, and multiple-complex genes).

Pathways and networks for coding and non-coding differentially expressed genes.

Next, we performed bioinformatic analysis to find cellular pathways involving differentially expressed (DE) protein-coding genes. Among the 56 cellular pathways we identified were those involved in the regulation of neurodegenerative diseases (e.g.,

HGD LGD

Pathways and Networks for Coding and Non-Coding Differentially Expressed Genes

Next, we performed bioinformatic analysis to find cellular pathways involving differentially expressed (DE) protein-coding genes. Among the 56 cellular pathways we identified were those involved in the regulation of neurodegenerative diseases (e.g., Alzheimer's disease), pathways involved in the regulation of cellular energy pathways (e.g., oxidative phosphorylation), and cellular metabolism (e.g., fatty acid metabolism). We also observed several cellular signaling pathways, cell adhesion pathways, as well as other pathways such as cell cycle or protein processing in endoplasmic reticulum (Supplemental Figure S3). Therefore, the HGD led to genomic modification in brain microvasculature pathways primarily by modulating the expression of genes involved in neurodegeneration, metabolism, and cell signaling. We also performed bioinformatics analyses of DEGs to identify potential transcription factors (TFs) whose activity could be modulated by the HGD and result in the observed genomic effects. Enrichr database analysis of the top 25 TFs and the DEGs regulated by them is shown in Supplemental Figure S4. The most statistically significant transcription factors were STAT3 (Signal Transducer and Activator of Transcription 3) involved in focal adhesion, DNMT1 (DNA Methyl Transferase 1), and PPARA (Peroxisome Proliferator-Activated Receptor Alpha) involved in Alzheimer's disease (Supplemental Table S2).

Microarray analysis also revealed that the HGD could induce differential expression of non-coding RNAs (miRNAs, lncRNAs, and snoRNAs). Of the 25 DE miRNAs in HGD vs. LGD (Supplemental Table S3), 12 were down-regulated (fold change -5.81 to -2.01) and 13 up-regulated (fold change 2.01 to 25.94). Using MIENTURNET software and database interrogation, we identified 442 potential target genes for 17 of the 25 miRNAs. The network of interactions between these miRNAs and their target genes is presented in Supplemental Figure S5. While most genes were the target of a single miRNA, we showed redundancy in that some genes were targets of two or three different miRNAs. Pathway analyses of miRNA target genes revealed they were involved in pathways regulating cell transduction, cell-cell adhesion, permeability, and neurofunction (Supplemental Figure S6). Eleven of the miRNA target genes pathways were in common with the protein-coding DEG pathways, such as Hypoxia-inducible factor 1 (HIF-1) signaling involved in endothelial cell function. Pathways unique to miRNA targets were primarily related to inflammation and cell signaling, whereas protein-coding DEGs specific pathways were primarily involved in cellular metabolism.

Together with miRNAs, our analysis also revealed 32 DE lncRNAs (Supplemental Table S4) following the HGD. Among these, 4 were down-regulated (fold-change -16.65 to -2.02), and 28 were up-regulated (fold-change 2.02 to 5.34). Using LncRRIsearch and Rtools CBRC databases, we were able to identify 458 potential target genes for 5 of the 32 lncRNAs (Supplemental Figure S7). Pathway analysis of these target genes showed that they were involved in pathways such as nitric oxide signaling, N-cadherin that regulate vascular endothelial function, as well as Alzheimer's disease (Supplemental Figure S8), which was also one of the 7 pathways in common with the miRNA targets pathways.

We also studied the expression of snoRNAs in the HGD compared to LGD. Among the 21 DE snoRNAs, 9 were up-regulated (fold change 2.36 to 19.44), and 12 were downregulated (fold change -10.37 to -2.07) (Supplemental Table S5). A literature review did not identify any known target genes for the DE snoRNAs.

Integrated Analysis of Differentially Expressed Genes, Key Pathways and Networks

Following the individual omic analysis, we performed integrated analysis of all the identified DEGs including mRNAs, miRNAs and their targets, lncRNAs and their targets, and the identified potential transcription factors. This analysis allowed us to obtain networks of interactions (Figure 4A), and showed that in comparison to the LGD, the HGD significantly impacted the expression levels of different RNA types which interact and form a large molecular network. This molecular network can have a significant effect on cellular functions. Therefore, in order to determine if there was a pattern of functional coordination in the molecular differential expression pattern of the HGD on brain microvessels, we performed integrated pathway and network analysis using DE mRNAs, miRNAs and lncRNAs targets for the HGD compared to LGD. This analysis indeed revealed differential regulation of 5 key cellular pathways including for neurodegenerative diseases (such as Alzheimer's disease), cell signaling pathways (such as PPAR signaling and phosphoinositide-3-kinase-protein kinase B (PI3K-Akt) signaling), cell adhesion and mobility (including focal adhesion), cellular metabolism (including oxidative phosphorylation and electron transport chain), and other cellular pathways (such as mRNA processing and oxidative damage) (Figure 4B and4C). Integrated pathways for DEGs, transcription factors, miRNAs and their targets, and lncRNAs and their targets are shown for the focal adhesion pathway (Figure 4D) and the Alzheimer's disease pathway (Figure 4E). Pathways are discussed in further detail in the Discussion section.

Effect of the Soluble Epoxide Hydrolase Inhibitor (sEHI) on the Hippocampal Microvascular Genome of Mice fed the High Glycemic Diet

In order to determine whether the soluble epoxide hydrolase inhibitor (sEHI) could inhibit the seemingly deleterious molecular effects of the HGD on hippocampal microvessels (upregulation of genes in pathways such as PPAR signaling, PI3K-Akt signaling that play an important role in oxidative stress, inflammation and Alzheimer's disease), we again performed hierarchical clustering of global gene expression profiles for the LGD and the HGD in the presence of the inhibitor (Figure 5). Since the gene expression profile of the LGD group did not substantially differ from the LGD with the inhibitor (Figure 5), the effect of the inhibitor was primarily analyzed in reference to the HGD in our analytic comparisons. Interestingly, in the presence of the sEHI, the gene expression profile of the HGD on hippocampal microvessels was nearly completely reversed and similar to that on the LGD.

Using PCA analysis, we were further able to show that the gene expression profile of hippocampal microvessels following the HGD was distinctly different from that of the HGD+sEHI (Figure 6A). We used loading plot to identify genes important to the separation of PCA (Figure 6B). Genes such as Gm24400, Resp18, Rheb, Brms1l, Ndufa13, and Ighv5-12-4 showed opposite expression in the HGD with inhibitor when compared to the HGD alone.

Statistical analysis of microarray data showed that there were a larger number of DEGs (1701) for the comparison of the HGD+sEHI vs. HGD (Supplemental Table S6) than there were for the HGD vs. LGD comparison. The sEHI primarily down-regulated both proteincoding and non-coding DEGs in the HGD (Figure 7A, and Supplemental Figures S9 andS10). We then performed correlation analysis between the DEGs of the HGD+sEHI vs. HGD and the HGD vs. LGD and identified a highly significant negative correlation with the sEHI (Figure 7B). This suggests that the inhibitor counteracts the effects of the HGD on differential gene expression. disease), we again performed hierarchical clustering of global gene expression profiles for the LGD and the HGD in the presence of the inhibitor (Figure 5). Since the gene expression profile of the LGD group did not substantially differ from the LGD with the inhibitor (Figure 5), the effect of the inhibitor was primarily analyzed in reference to the HGD in our analytic comparisons. Interestingly, in the presence of the sEHI, the gene expression profile of the HGD on hippocampal microvessels was nearly completely reversed and similar to that on the LGD. Most of the protein-coding genes down-regulated by the sEHI were involved in similar pathways (such as Alzheimer's disease, oxidative phosphorylation, and fatty acid metabolism) activated by the HGD alone (Supplemental Figure S11), suggesting that the inhibitor may offset the effect of the HGD on differential gene expression by targeting similar pathways. We also found 17 TFs in common between the HGD with and without inhibitors such as NEUROG3 (Neurogenin 3), MECP2 (Methyl-CpG Binding Protein 2), KLF9 (Kruppel Like Factor 9), and PPARA (Supplemental Table S2). Target genes of these TFs were up-regulated by the HGD but down-regulated by the inhibitor (Supplemental Figure S12). These target genes were involved in neurological disease pathways (e.g., Alzheimer's disease), oxidative phosphorylation (e.g., electron transport chain), apoptosis-related pathways (e.g., proteasome degradation), and cell signaling pathways (e.g., HF-1 signaling). Principal Component Analysis (PCA) and loading plot of genes expressed in hippocampal microvessels for the high glycemic diet (HGD) with and without soluble epoxide hydrolase inhibitor (sEHI). (A) PCA scatter plot of the microarray data shows the trends of the expression profiles of the hippocampal microvasculature of the high glycemic diet with (HGD+sEHI, green circles) and without (HGD, red circles) the soluble epoxide hydrolase inhibitor compared to the low glycemic diet (LGD, blue circles). The PCA plot captures the variance in a dataset in terms of principal components and displays the most significant of these on the x and y axes. The percentages of the total variation that are accounted for by the 1st and 2nd principal components are shown on the x-and y-axes labels. The data are shown for three biological replicates for each dietary group. (B) Loading plot of genes (violet circles) relevant to the separation of PCA. Blue arrows show box plots with the expression levels of a few genes (Gm24400, Resp18, Rheb, Brms1l, Ndufa13, Ighv5-12-4) that contribute to separation of the three dietary groups (HGD, HGD+sEHI, LGD). In the box plots, the black dots represent gene expression levels, the notch indicates 95% confidence interval around the median of each group, and the yellow diamond indicates the average gene expression of each group.

LGD HGD+sEHI HGD

A)

PCA Plot Most of the protein-coding genes down-regulated by the sEHI were involved in similar pathways (such as Alzheimer's disease, oxidative phosphorylation, and fatty acid metabolism) activated by the HGD alone (Supplemental Figure S11), suggesting that the inhibitor may offset the effect of the HGD on differential gene expression by targeting similar pathways. We also found 17 TFs in common between the HGD with and without inhibitors such as NEUROG3 (Neurogenin 3), MECP2 (Methyl-CpG Binding Protein 2), KLF9 (Kruppel Like Factor 9), and PPARA (Supplemental Table S2). Target genes of these TFs were up-regulated by the HGD but down-regulated by the inhibitor (Supplemental Figure S12). These target genes were involved in neurological disease pathways (e.g., Alzheimer's disease), oxidative phosphorylation (e.g., electron transport chain), apoptosis-related pathways (e.g., proteasome degradation), and cell signaling pathways (e.g., HF-1 signaling).

A higher number of DE non-coding RNAs (171 miRNAs, 127 snoRNAs, and 80 lncRNAs) were modulated by the sEHI (Supplemental Tables S3, S4, and S5, respectively) than the HGD alone. Most of the DE miRNAs that were up-regulated by the HGD were down-regulated by the inhibitor by targeting 300 genes (Supplemental Figure S13) A higher number of DE non-coding RNAs (171 miRNAs, 127 snoRNAs, and 80 lncRNAs) were modulated by the sEHI (Supplemental Tables S3,S4, and S5, respectively) than the HGD alone. Most of the DE miRNAs that were up-regulated by the HGD were downregulated by the inhibitor by targeting 300 genes (Supplemental Figure S13) that were involved in apoptosis-related pathways (such as PTEN dependent cell cycle arrest and Ubiquitin mediated proteolysis) (Supplemental Figure S14) and pathways related to cell degradation and death. The majority of DE lncRNAs that were up-regulated by the HGD were down-regulated by the inhibitor and again targeted a greater number of genes (785) compared to the HGD alone (Supplemental Figure S15) and were mainly involved in cell signaling (Supplemental Figure S16).

We then performed an integrated analysis of all the identified DEGs, including mR-NAs, miRNAs and their targets, lncRNAs and their targets, and the potential TFs. We obtained a large molecular network of interactions of the different types of RNAs modulated by the inhibitor (Figure 8A). Integrated pathway analysis of this molecular network (Figure 8B) revealed that the inhibitor regulated similar pathways as the HGD, namely neuro-related pathways (such as Alzheimer's disease), cell signaling, cell adhesion (such as focal adhesion), and cellular metabolism, but once again, in opposite directions to the HGD. In toto, these data suggest that the inhibitor may counteract the effect of an HGD by robustly reversing the deleterious pattern of differential gene expression on cellular function and on the operative TFs, particularly by down-regulating the protein-coding and non-protein-coding DEGs that are up-regulated by the HGD.

Correlation of Genomics with Human Alzheimer's Disease Data

To evaluate the genomic modifications induced by the HGD on the development of neurodegenerative diseases, we analyzed the correlations between the genes modulated by the HGD and the genes associated with Alzheimer's disease. We analyzed the available GEO datasets with GEO2R and differentially expressed mRNAs were correlated with those extracted from the analyzed publications using the R program. These analyses showed that the gene expression profile in the brain of patients with Alzheimer's disease (GSE118553) was correlated, with high statistical significance, to the gene expression profile identified in our mice after the consumption of the HGD (Figure 9A). The gene expression profile of mice on the HGD in our study was also correlated with the gene expression profile obtained in another study (GSE132903) of brains of patients presenting with Alzheimer's disease (Figure 9B). This analysis suggests that consumption of the HGD induces changes in the expression of genes in murine hippocampus that correlates with the development of Alzheimer's disease in humans. Using the same approach, we also tested correlations between the genes modulated in our study by the HGD in the presence of the sEHI and the gene expression profiles obtained in both of the above noted genomic studies using the brains of patients with Alzheimer's disease. Interestingly, we observed a negative correlation between the human AD data and the gene expression profile of our brain microvasculature for animals on the HGD with the sEHI (Figure 9C,D). These results suggest that the inhibitor can prevent, at least partially, changes in the expression of genes induced by the HGD that are associated with the development of Alzheimer's disease in humans and can therefore exert a potentially protective effect against this neurological disorder. that were involved in apoptosis-related pathways (such as PTEN dependent cell cycle arrest and Ubiquitin mediated proteolysis) (Supplemental Figure S14) and pathways related to cell degradation and death. The majority of DE lncRNAs that were up-regulated by the HGD were down-regulated by the inhibitor and again targeted a greater number of genes (785) compared to the HGD alone (Supplemental Figure S15) and were mainly involved in cell signaling (Supplemental Figure S16).

We then performed an integrated analysis of all the identified DEGs, including mRNAs, miRNAs and their targets, lncRNAs and their targets, and the potential TFs. We obtained a large molecular network of interactions of the different types of RNAs modulated by the inhibitor (Figure 8A). Integrated pathway analysis of this molecular network (Figure 8B) revealed that the inhibitor regulated similar pathways as the HGD, namely neuro-related pathways (such as Alzheimer's disease), cell signaling, cell adhesion (such as focal adhesion), and cellular metabolism, but once again, in opposite directions to the HGD. In toto, these data suggest that the inhibitor may counteract the effect of an HGD by robustly reversing the deleterious pattern of differential gene expression on cellular function and on the operative TFs, particularly by down-regulating the protein-coding and non-protein-coding DEGs that are up-regulated by the HGD. genomic studies using the brains of patients with Alzheimer's disease. Interestingly, we observed a negative correlation between the human AD data and the gene expression profile of our brain microvasculature for animals on the HGD with the sEHI (Figures 9C and9D). These results suggest that the inhibitor can prevent, at least partially, changes in the expression of genes induced by the HGD that are associated with the development of Alzheimer's disease in humans and can therefore exert a potentially protective effect against this neurological disorder.

A)

Interaction of Effects of the Glycemic Diet and the Soluble Epoxide Hydrolase Inhibitor

In order to more clearly define how the sEHI modulated the response of hippocampal microvessels to an HGD, we used two-factor analysis to identify the DEGs with a significant interaction between the effects of the diet and the effects of the inhibitor. For this analysis, the LGD+sEHI group was included as it was necessary to use this additional control group for the two-factor analysis and to determine the delta fold change [(HGD+sEHI vs.

LGD+sEHI) vs. (HGD vs. LGD)]. A positive delta fold change indicates a gene that is "more up-regulated" or "less down-regulated" by HGD in sEHI treated mice (compared to those with no inhibitor). A negative delta fold change indicates a gene that is "more downregulated" or "less up-regulated" by an HGD in sEHI treated mice (compared to those with no inhibitor). We identified a total of 1543 DEGs (Supplemental Table S7) that had a significant interaction between diet and inhibitor, the majority of which had a negative delta fold change. In many cases, the direction of the regulation of the DEGs by the HGD was opposite in direction in the presence of the inhibitor when compared to the effects of the HGD on that same gene in the absence of the inhibitor. The DEGs with a negative delta fold change were diversified between RNA types, including coding, non-coding, multiple complex, and pseudogenes. The 583 DEGs with a positive delta fold change were primarily non-coding genes. Pathway analyses using the KEGG and WikiPathway databases for all of the coding DEGs (positive and negative delta fold change) demonstrated that the sEHI modulated the effects of the HGD for pathways involved in the electron transport chain, oxidative phosphorylation, focal adhesion, HIF-1 signaling pathway, and pathways active in neurological diseases (Alzheimer's pathway) (Figure 10A).

Discussion

In this study, we characterized the effect of a high glycemic diet (HGD) on murine neurovascular function and cognition using hippocampal microvascular multiomics of mice exposed to high or low glycemic index diets. Utilizing similar techniques of large-scale transcriptome gene profiling, we have previously shown the effect of lipid stress and a Western diet (WD) on hippocampal microvascular injury and cognitive decline in male and female mice ApoE-/-and LDL-R-/-mice [START_REF] Aung | Lipotoxic brain microvascular injury is mediated by activating transcription factor 3-dependent inflammatory and oxidative stress pathways[END_REF][START_REF] Nuthikattu | The Western Diet Regulates Hippocampal Microvascular Gene Expression: An Integrated Genomic Analyses in Female Mice[END_REF][START_REF] Nuthikattu | Lipotoxic Injury Differentially Regulates Brain Microvascular Gene Expression in Male Mice[END_REF][START_REF] Nuthikattu | Sex-Dependent Molecular Mechanisms of Lipotoxic Injury in Brain Microvasculature: Implications for Dementia[END_REF]. However, the molecular footprint of the glycemic component of the WD on the brain microvasculature, and whether and how the HGD may contribute to microvascular injury, was up to now unknown. A HGD is a risk factor for dementia [START_REF] Gentreau | High Glycemic Load Is Associated with Cognitive Decline in Apolipoprotein E ε4 Allele Carriers[END_REF] and is associated with poor cognitive performance [START_REF] Saikrishna | Combined Administration of Monosodium Glutamate and High Sucrose Diet Accelerates the Induction of Type 2 Diabetes, Vascular Dysfunction, and Memory Impairment in Rats[END_REF][START_REF] Wong | High-sucrose diets in male rats disrupt aspects of decision making tasks, motivation and spatial memory, but not impulsivity measured by operant delay-discounting[END_REF], but whether the effect can be reversed has also been unknown. Given that in murine models of diabetes, soluble epoxide hydrolase inhibitors (sEHI) reduce neuroinflammation and cognitive decline [START_REF] Minaz | An inhibitor of soluble epoxide hydrolase ameliorates diabetes-induced learning and memory impairment in rats[END_REF][START_REF] Wu | Inhibition of soluble epoxide hydrolase (sEH) protects hippocampal neurons and reduces cognitive decline in type 2 diabetic mice[END_REF][START_REF] Ghosh | An epoxide hydrolase inhibitor reduces neuroinflammation in a mouse model of Alzheimer's disease[END_REF], we also aimed to examine their role in our model.

We used a relevant model of glycemic index and demonstrated the expected significant differences in weight, cholesterol, glucose, and insulin levels consistent with what has been published previously for these experimental models [START_REF] Yeh | A high-sucrose diet aggravates Alzheimer's disease pathology, attenuates hypothalamic leptin signaling, and impairs food-anticipatory activity in APPswe/PS1dE9 mice[END_REF][START_REF] Van Schothorst | Effects of a high-fat, low-versus high-glycemic index diet: Retardation of insulin resistance involves adipose tissue modulation[END_REF][START_REF] Zhu | The Short-Term Effect of a High-Glycemic Diet on Mouse Obesity and Intestinal Microbiota Composition[END_REF][START_REF] Sousa | The effects of diet composition and chronic obesity on muscle growth and function[END_REF]. To our knowledge, the present study shows for the first time that the HGD has distinct differential effects on gene expression that are different from the LGD and characterized by up-regulation of differentially expressed protein-coding genes, transcription factors, as well as non-protein-coding genes (miRNAs, snoRNAs, and lncRNAs), that are involved in five major cellular pathways: neurodegenerative diseases (Alzheimer's disease), endothelial cell function (focal adhesion), cell signaling (PPAR signaling, PI3K-Akt signal- An additional finding of the two-factor analysis was that when the individual DEGs of the relevant pathways shown in Figure 10A (77 DEGs) were examined, it again became apparent that the sEHI reverted the gene expression profile of brain hippocampal microvascular endothelium of mice on the HGD to a profile similar to that of the LGD. Examples of three genes (Rock1, Ndufa2, and Cox5b) are shown in Figure 10B. These DEGs were selected because they were each represented in multiple relevant pathways including focal adhesion (Rock1), electron transport chain (Ndufa2 and Cox5b), and oxidative phosphorylation (Ndufa2 and Cox5b).

Discussion

In this study, we characterized the effect of a high glycemic diet (HGD) on murine neurovascular function and cognition using hippocampal microvascular multiomics of mice exposed to high or low glycemic index diets. Utilizing similar techniques of largescale transcriptome gene profiling, we have previously shown the effect of lipid stress and a Western diet (WD) on hippocampal microvascular injury and cognitive decline in male and female mice ApoE-/-and LDL-R-/-mice [START_REF] Aung | Lipotoxic brain microvascular injury is mediated by activating transcription factor 3-dependent inflammatory and oxidative stress pathways[END_REF][START_REF] Nuthikattu | The Western Diet Regulates Hippocampal Microvascular Gene Expression: An Integrated Genomic Analyses in Female Mice[END_REF][START_REF] Nuthikattu | Lipotoxic Injury Differentially Regulates Brain Microvascular Gene Expression in Male Mice[END_REF][START_REF] Nuthikattu | Sex-Dependent Molecular Mechanisms of Lipotoxic Injury in Brain Microvasculature: Implications for Dementia[END_REF]. However, the molecular footprint of the glycemic component of the WD on the brain microvasculature, and whether and how the HGD may contribute to microvascular injury, was up to now unknown. A HGD is a risk factor for dementia [START_REF] Gentreau | High Glycemic Load Is Associated with Cognitive Decline in Apolipoprotein E ε4 Allele Carriers[END_REF] and is associated with poor cognitive performance [START_REF] Saikrishna | Combined Administration of Monosodium Glutamate and High Sucrose Diet Accelerates the Induction of Type 2 Diabetes, Vascular Dysfunction, and Memory Impairment in Rats[END_REF][START_REF] Wong | High-sucrose diets in male rats disrupt aspects of decision making tasks, motivation and spatial memory, but not impulsivity measured by operant delay-discounting[END_REF], but whether the effect can be reversed has also been unknown. Given that in murine models of diabetes, soluble epoxide hydrolase inhibitors (sEHI) reduce neuroinflammation and cognitive decline [START_REF] Minaz | An inhibitor of soluble epoxide hydrolase ameliorates diabetes-induced learning and memory impairment in rats[END_REF][START_REF] Wu | Inhibition of soluble epoxide hydrolase (sEH) protects hippocampal neurons and reduces cognitive decline in type 2 diabetic mice[END_REF][START_REF] Ghosh | An epoxide hydrolase inhibitor reduces neuroinflammation in a mouse model of Alzheimer's disease[END_REF], we also aimed to examine their role in our model.

We used a relevant model of glycemic index and demonstrated the expected significant differences in weight, cholesterol, glucose, and insulin levels consistent with what has been published previously for these experimental models [START_REF] Yeh | A high-sucrose diet aggravates Alzheimer's disease pathology, attenuates hypothalamic leptin signaling, and impairs food-anticipatory activity in APPswe/PS1dE9 mice[END_REF][START_REF] Van Schothorst | Effects of a high-fat, low-versus high-glycemic index diet: Retardation of insulin resistance involves adipose tissue modulation[END_REF][START_REF] Zhu | The Short-Term Effect of a High-Glycemic Diet on Mouse Obesity and Intestinal Microbiota Composition[END_REF][START_REF] Sousa | The effects of diet composition and chronic obesity on muscle growth and function[END_REF]. To our knowledge, the present study shows for the first time that the HGD has distinct differential effects on gene expression that are different from the LGD and characterized by up-regulation of differentially expressed protein-coding genes, transcription factors, as well as non-proteincoding genes (miRNAs, snoRNAs, and lncRNAs), that are involved in five major cellular pathways: neurodegenerative diseases (Alzheimer's disease), endothelial cell function (focal adhesion), cell signaling (PPAR signaling, PI3K-Akt signaling), cellular metabolism (oxidative phosphorylation, electron transport chain), and other cellular pathways (mRNA processing, oxidative damage). Below we discuss our multi-omics integrative analysis for the focal adhesion and Alzheimer's disease pathways to illustrate the complexity and characterize the molecular regulation in response to the HGD.

Integrative analyses allowed us to identify the protein-coding genes, TFs, and the non-coding genes modulated by the HGD in the focal adhesion pathway. Cdc42 (Cell Division Control Protein Homolog 42) is a protein-coding DEG up-regulated by the HGD and targeted by STAT3 (Signal Transducer and Activator of Transcription 3) transcription factor. Cdc42 is up-regulated in aging endothelial cells via the focal adhesion pathway [START_REF] Cho | Morphological adjustment of senescent cells by modulating caveolin-1 status[END_REF][START_REF] Flentje | Small GTPases and Their Role in Vascular Disease[END_REF] and mediates the activation of proinflammatory genes [START_REF] Ito | A crucial role for CDC42 in senescence-associated inflammation and atherosclerosis[END_REF]. STAT3 transcription factor is activated by oxidative stress [START_REF] Zouein | STAT3 and Endothelial Cell-Cardiomyocyte Dialog in Cardiac Remodeling[END_REF][START_REF] Bourgeais | Oxidative metabolism in cancer: A STAT affair?[END_REF] and is associated with neurovascular inflammation and Alzheimer's disease [START_REF] Millot | STAT3 inhibition protects against neuroinflammation and BACE1 upregulation induced by systemic inflammation[END_REF]. The HGD also up-regulated differentially expressed miRNA miR-1902 that targets ITGB (Integrin β). Activation of ITGB induces a reduction in endothelial barrier function by destabilizing intercellular junctions and abnormal remodeling of extracellular matrix (ECM) [START_REF] Faurobert | CCM1-ICAP-1 complex controls β1 integrin-dependent endothelial contractility and fibronectin remodeling[END_REF]. Another example of differentially expressed miRNAs up-regulated by the HGD is miR-5125 and miR-692 which target PKC (Protein Kinase C). PKC is associated with endothelial dysfunction, insulin resistance [START_REF] Naruse | Activation of vascular protein kinase C-beta inhibits Akt-dependent endothelial nitric oxide synthase function in obesity-associated insulin resistance[END_REF], and neuroinflammation [START_REF] Tang | Protein kinase C-delta inhibition protects blood-brain barrier from sepsis-induced vascular damage[END_REF]. The HGD also up-regulated differential expression of LncRNA Gm6117 that targets MLCK (Myosin Light Chain Kinase). Activation of MLCK induces endothelial vascular hyperpermeability, microvascular barrier dysfunction, and inflammation [START_REF] Shen | Myosin light chain kinase in microvascular endothelial barrier function[END_REF]. Another interesting pathway activated by the HGD was the vascular endothelial growth factor (VEGF) signaling pathway. VEGF is associated with an increase in microvascular permeability by activating PKC [START_REF] Wu | VEGF induces NO-dependent hyperpermeability in coronary venules[END_REF], which then signals the activation of MLCK [START_REF] Shen | Myosin light chain kinase in microvascular endothelial barrier function[END_REF]. Taken together, this molecular cascade suggests that an HGD may promote endothelial cell dysfunction by microvascular hyper permeability and activating proinflammatory TFs, protein, and non-protein-coding genes.

We also performed an integrated analysis of the Alzheimer's disease pathway and identified TFs, protein-coding, and non-coding RNAs modulated by the HGD. Mitochondrial complex IV protein-coding DEG was up-regulated by the HGD and targeted by DNMT1 (DNA Methyl Transferase 1) transcription factor. Increased expression of complex IV has been observed in the hippocampus of AD patients [START_REF] Manczak | Differential expression of oxidative phosphorylation genes in patients with Alzheimer's disease: Implications for early mitochondrial dysfunction and oxidative damage[END_REF], and DNMT1 expression is up-regulated in late-onset AD [START_REF] Di Francesco | Global changes in DNA methylation in Alzheimer's disease peripheral blood mononuclear cells[END_REF]. Alterations in mitochondrial function are also associated with diabetes [START_REF] Potenza | Diabetes and Alzheimer's Disease: Might Mitochondrial Dysfunction Help Deciphering the Common Path?[END_REF]. Another example of a DEG up-regulated by the HGD was LRP1 (LDLR-related protein). Expression of LRP1 is increased in the hippocampus of AD patients [START_REF] Akram | Association of ApoE and LRP mRNA levels with dementia and AD neuropathology[END_REF], particularly in the microvasculature close to amyloid plaques in AD brains [START_REF] Ruzali | LRP1 expression in cerebral cortex, choroid plexus and meningeal blood vessels: Relationship to cerebral amyloid angiopathy and APOE status[END_REF][START_REF] Donahue | LRP-1, and amyloid-beta protein in Alzheimer's disease[END_REF]. DE miRNA, miR-5125 was activated by the HGD and is known to target DKK1 (Dickkopf-1). Increased expression of DKK1 has been detected in the plasma and brain of AD patients and AD transgenic mice [START_REF] Caricasole | Induction of Dickkopf-1, a negative modulator of the Wnt pathway, is associated with neuronal degeneration in Alzheimer's brain[END_REF][START_REF] Rosi | Increased Dickkopf-1 expression in transgenic mouse models of neurodegenerative disease[END_REF]. Another target of miR-5125 was LPL (Lipoprotein lipase). LPL accumulates in senile plaques [START_REF] Rebeck | Diverse senile plaque-associated proteins are ligands of an apolipoprotein E receptor, the alpha 2-macroglobulin receptor/low-density-lipoprotein receptor-related protein[END_REF], and its expression has been shown to be increased in the hippocampus of AD mice [START_REF] Blain | Could lipoprotein lipase play a role in Alzheimer's disease?[END_REF]. LncRNAs Gm6117 and Gm16339 were also up-regulated by the HGD and target NMDAR (N-methyl-D-aspartate receptor). An increase in NMDAR activity induces apoptosis and neurodegeneration in AD [START_REF] Wang | Role of Glutamate and NMDA Receptors in Alzheimer's Disease[END_REF]. Yet another interesting pathway modulated by the HGD was the PI3K-Akt signaling pathway. Activation of the PI3K/Akt pathway promotes oxidative stress-mediated cell death [START_REF] Nogueira | Akt determines replicative senescence and oxidative or oncogenic premature senescence and sensitizes cells to oxidative apoptosis[END_REF]. Thus, our findings suggest that an HGD may be involved in previously defined Alzheimer's disease-related pathways via deleterious effects on the brain microvasculature by promoting mitochondrial dysfunction, apoptosis, and neurodegeneration, and similar to the focal adhesion pathway, via up-regulation of pathway-associated TFs, protein-coding, and non-coding genes.

Our work with the soluble epoxide hydrolase inhibitor (sEHI) trans-4-[4-(3-adamantan-1-yl-ureido)-cyclohexyloxy]-benzoic acid (t-AUCB) showed several novel findings. sEH enzyme degrades anti-inflammatory and neuroprotective fatty acid epoxides (EETs) to biologically less active diols (DHETs) which are known to be involved in vascular cognitive impairment and Alzheimer's disease [START_REF] Nelson | Role of soluble epoxide hydrolase in age-related vascular cognitive decline[END_REF][START_REF] Yu | Soluble Epoxide Hydrolase-Derived Linoleic Acid Oxylipins in Serum Are Associated with Periventricular White Matter Hyperintensities and Vascular Cognitive Impairment[END_REF]. In animal models of stroke, sEH inhibitors have been shown to be neuroprotective [START_REF] Simpkins | Soluble epoxide inhibition is protective against cerebral ischemia via vascular and neural protection[END_REF][START_REF] Shaik | Soluble epoxide hydrolase inhibitor trans-4-[4-(3-adamantan-1-yl-ureido)-cyclohexyloxy]-benzoic acid is neuroprotective in rat model of ischemic stroke[END_REF]. In our study, we show for the first time that the sEHI reversed the deleterious phenotype of the HGD by down-regulating the differentially expressed genes up-regulated by the HGD and their TFs, and also by differential expression of non-coding RNAs (miRNA, snoRNA, and lncRNA) and that this proceeds by targeting the same pathways up-regulated by the HGD. When looking at specific DEGs targeted by the inhibitor, we identified interesting DEGs. For example, the presence of the sEHI down-regulated the expression of Cdc42, and the deletion of Cdc42 is known to have protective effects against chronic inflammation in endothelial cells [START_REF] Flentje | Small GTPases and Their Role in Vascular Disease[END_REF]. The sEHI also down-regulated the expression of the genes targeted by the PPAR (Peroxisome Proliferator-Activated Receptor) transcription factor, which is known to play a protective role against neuroinflammation in AD [START_REF] Heneka | Impact and Therapeutic Potential of PPARs in Alzheimer's Disease[END_REF]. Furthermore, the sEHI down-regulated DE miRNAs, miR-5125, and miR-692, and their targets, DKK1 and PKC, respectively. Inhibition of DKK1 has been shown to improve cognitive impairment [START_REF] Ortiz-Matamoros | Chronic infusion of Wnt7a, Wnt5a and Dkk-1 in the adult hippocampus induces structural synaptic changes and modifies anxiety and memory performance[END_REF] and protect against neurotoxicity [START_REF] Fortress | Hippocampal Wnt Signaling: Memory Regulation and Hormone Interactions[END_REF]. PKC inhibition also protects against vascular injury of the blood-brain barrier [START_REF] Tang | Protein kinase C-delta inhibition protects blood-brain barrier from sepsis-induced vascular damage[END_REF]. Additionally, the sEHI also down-regulated DE lncRNA Gm6117 and its target, MLCK. MLCK inhibitors play a protective role against vascular hyperpermeability and microvascular dysfunction [START_REF] Shen | Myosin light chain kinase in microvascular endothelial barrier function[END_REF]. These findings indicate that the sEHI served to protect against the neurotoxicity, endothelial-and neuro-inflammation, and blood-brain barrier permeability injury associated with the HGD. Furthermore, our correlation analysis of the DEGs modulated by the HGD and inhibitor with the gene expression profiles obtained from the brains of Alzheimer's disease patients showed that the HGD is positively correlated, and the inhibitor is negatively correlated with the genes associated with Alzheimer's disease.

Our 2-factor interactive analysis served to further identify DEGs having a significant interaction between diet and the soluble epoxide hydrolase inhibitor. These DEGs were involved in Alzheimer's disease, electron transport chain, oxidative phosphorylation, and focal adhesion. The sEHI also was shown to restore the gene expression profile of hippocampal microvessels of mice on the HGD to a profile similar to that of the LGD, including for genes such as Rock1 (Rho-associated protein kinase 1), which was increased by an HGD and decreased by the sEHI. ROCK1 is activated in Alzheimer's disease, and reduction of ROCK1 protects against AD by depleting amyloid-β levels in the brain [START_REF] Henderson | Rho-associated protein kinase 1 (ROCK1) is increased in Alzheimer's disease and ROCK1 depletion reduces amyloid-β levels in brain[END_REF]. In vitro exposure of brain microvascular endothelial cells to EETs, which generally increase after inhibition of sEH, has been shown to inhibit rho kinase (ROCK) activity [START_REF] Gupta | Soluble epoxide hydrolase: Sex differences and role in endothelial cell survival[END_REF]. Taken together, our molecular results suggest that the soluble epoxide hydrolase inhibitor may offset the effect of an HGD by robustly reversing the deleterious pattern of differential gene expression to that similar to an LGD and likely thereby protecting against vascular injury, neuroinflammation, and neurodegeneration. Furthermore, the sEHI was associated with improved cognition by down-regulating the TFs, protein-coding, and non-protein-coding DEGs that were up-regulated by the HGD.

Conclusions and Clinical Implications

In alignment with our hypothesis, our results show that an HGD has deleterious effects on differential gene expression of brain microvessels in memory centers in the brain by up-regulating the protein-coding and non-coding genes involved in alterations of mitochondrial function, oxidation, inflammation, and microvascular dysfunction. These processes play an important role in the pathophysiology of dementia. We also showed that inhibition of sEH protects against neuroinflammation, apoptosis, vascular injury, and cognitive decline by down-regulating the DEGs up-regulated by an HGD. The clinical role of sEHIs has been reviewed [START_REF] Zarriello | Humble beginnings with big goals: Small molecule soluble epoxide hydrolase inhibitors for treating CNS disorders[END_REF]. In clinical trials, sEHIs play a protective role in hypertension [START_REF] Chen | Pharmacokinetics and pharmacodynamics of AR9281, an inhibitor of soluble epoxide hydrolase, in single-and multiple-dose studies in healthy human subjects[END_REF], and in animal models, they demonstrate neuroprotection in stroke [START_REF] Simpkins | Soluble epoxide inhibition is protective against cerebral ischemia via vascular and neural protection[END_REF][START_REF] Shaik | Soluble epoxide hydrolase inhibitor trans-4-[4-(3-adamantan-1-yl-ureido)-cyclohexyloxy]-benzoic acid is neuroprotective in rat model of ischemic stroke[END_REF], cerebral hypoperfusion, and in diabetes type 1 and type 2 [START_REF] Chen | Soluble epoxide hydrolase inhibition Promotes White Matter Integrity and Long-Term Functional Recovery after chronic hypoperfusion in mice[END_REF][START_REF] Hao | Soluble epoxide hydrolase inhibition alleviated cognitive impairments via NRG1/ErbB4 signaling after chronic cerebral hypoperfusion induced by bilateral carotid artery stenosis in mice[END_REF][START_REF] Matin | Soluble epoxide hydrolase inhibition improves cognitive function and parenchymal artery dilation in a hypertensive model of chronic cerebral hypoperfusion[END_REF][START_REF] Minaz | An inhibitor of soluble epoxide hydrolase ameliorates diabetes-induced learning and memory impairment in rats[END_REF][START_REF] Wu | Inhibition of soluble epoxide hydrolase (sEH) protects hippocampal neurons and reduces cognitive decline in type 2 diabetic mice[END_REF]. Our studies show that sEHIs may also be promising therapeutic targets in the microvascular endothelial dysfunction that accompanies Alzheimer's and vascular dementias.

Supplementary Materials:

The following are available online at https://www.mdpi.com/article/10 .3390/nu13113913/s1, Figure S1: Volcano plot of differential gene expression changes in hippocampal microvessels for the high glycemic diet (HGD) compared to the low glycemic diet (LGD), Figure S2: Distribution of differentially expressed RNAs in hippocampal microvessels for the high glycemic diet (HGD) when compared to the low glycemic diet (LGD), Figure S3: Histogram of differentially expressed protein coding genes pathways in hippocampal microvessels for the high glycemic diet (HGD) when compared to the low glycemic diet (LGD), Figure S4. Target gene networks of differentially expressed transcription factors (TFs) in hippocampal microvessels for the high glycemic diet (HGD) compared to the low glycemic diet (LGD), Figure S5. Target gene networks of differentially expressed miRNAs in hippocampal microvessels with the high glycemic diet (HGD) compared to the low glycemic diet (LGD), Figure S6: Histogram of differentially expressed miRNA targets pathways in hippocampal microvessels with the high glycemic diet (HGD) when compared to the low glycemic diet (LGD), Figure S7. Target gene networks of differentially expressed LncRNAs in hippocampal microvessels with the high glycemic diet (HGD) compared to the low glycemic diet (LGD), Figure S8: Histogram of differentially expressed lncRNA targets pathways in hippocampal microvessels with the high glycemic diet (HGD) when compared to the low glycemic diet (LGD), Figure S9: Volcano plot of differential gene expression changes in hippocampal microvessels for the high glycemic diet (HGD) with soluble epoxide hydrolase inhibitor (sEHI) compared to without sEHI treatment, Figure S10: Distribution of differentially expressed RNAs in hippocampal microvessels from the high glycemic diet (HGD) with soluble epoxide hydrolase inhibitor (sEHI) compared to without sEHI treatment, Figure S11: Histogram of protein coding differentially expressed genes pathways in hippocampal microvessels for the high glycemic diet (HGD) with soluble epoxide hydrolase inhibitor (sEHI) compared to without sEHI treatment, Figure S12: Target gene networks of differentially expressed transcription factors (TFs) in hippocampal microvessels for the high glycemic diet (HGD) with and without soluble epoxide hydrolase inhibitor (sEHI), Figure S13: Target gene networks of differentially expressed miRNAs in hippocampal microvessels for the high glycemic diet (HGD) with and without soluble epoxide hydrolase inhibitor (sEHI), Figure S14: Histogram of differentially expressed miRNA targets pathways in hippocampal microvessels for the HGD with soluble epoxide hydrolase inhibitor (sEHI) compared to without sEHI treatment, Figure S15: Target gene networks of differentially expressed LncRNAs in hippocampal microvessels for the high glycemic diet (HGD) with and without soluble epoxide hydrolase inhibitor (sEHI), Figure S16: Histogram of differentially expressed lncRNA targets pathways in hippocampal microvessels for the high glycemic diet (HGD) with soluble epoxide hydrolase inhibitor (sEHI) compared to without the sEHI; Table S1: Differentially expressed genes for the high glycemic diet (HGD) compared to the low glycemic diet (LGD), Table S2: Top 30 transcription factors (TFs) differentially expressed by a high glycemic diet (HGD) with and without soluble epoxide hydrolase inhibitor (sEHI) in hippocampal microvascular endothelium, Table S3: Effect of the high glycemic diet (HGD) with and without soluble epoxide hydrolase inhibitor (sEHI) on the expression of microRNAs (miRNAs) in hippocampal microvessels, Table S4: Effect of the high glycemic diet (HGD) with and without soluble epoxide hydrolase inhibitor (sEHI) on the expression of long non-coding RNAs (lncRNAs) in hippocampal microvessels, Table S5: Effect of the high glycemic diet (HGD) with and without soluble epoxide hydrolase inhibitor (sEHI) on the expression of small nucleolar RNAs (snoRNAs) in hippocampal microvessels, Table S6: Differentially expressed genes for the high glycemic diet (HGD) with soluble epoxide hydrolase inhibitor (sEHI) compared to without sEHI, Table S7: Differentially expressed genes with a significant interaction of diet and inhibitor.

Author Contributions: Individual contributions of the authors were as follows. Conceptualization, A.V. and J.R.; methodology, S.N., A.V.; software, S.N., A.V., D.M., J.E.N.; validation, A.V., S.N., D.M. and J.E.N.; formal analysis, S.N., D.M., A.V., J.E.N.; investigation, S.N.; resources, A.V. and J.R.; data curation, A.V., D.M., S.N., J.E.N.; writing-original draft preparation, A.V. and S.N.; writing-review and editing, A.V., S.N., D.M., J.E.N.; visualization, A.V., S.N., D.M., J.E.N.; supervision, A.V., J.R., D.M.; project administration, S.N.; funding acquisition, A.V. and J.R. All authors have read and agreed to the published version of the manuscript.

Nutrients 2021 , 31 Figure 1 .

 2021311 Figure 1. Research methodology flow chart. The flow chart shows the steps in the genomic analysis for: A) the high glycemic diet compared to the low glycemic diet (HGD vs LGD), and B) the high glycemic diet with the soluble epoxide hydrolase inhibitor (sEHI) compared to without the sEHI (HGD + sEHI vs HGD).

Figure 1 .

 1 Figure 1. Research methodology flow chart. The flow chart shows the steps in the genomic analysis for: (A) the high glycemic diet compared to the low glycemic diet (HGD vs LGD), and (B) the high glycemic diet with the soluble epoxide hydrolase inhibitor (sEHI) compared to without the sEHI (HGD + sEHI vs HGD).

3. 1 .

 1 Effect of the High Glycemic Diet on the Hippocampal Microvascular Genome 3.1.1

Figure 2 .

 2 Figure 2. Principal Component Analysis (PCA) and loading plot of genes expressed in hippocampal microvessels of the high glycemic diet (HGD) and the low glycemic diet (LGD). A) PCA scatter plot of the microarray data shows the trends of the expression profiles of the hippocampal microvasculature in the high glycemic diet (HGD, red circles) and low glycemic diet (LGD, green circles), respectively. The PCA plot captures the variance in a dataset in terms of principal components and displays the most significant of these on the x and y axes. The percentages of the total variation that are accounted for by the 1st and 2nd principal components are shown on the x-and y-axes labels. The data are shown for three biological replicates for each dietary group. B) Loading plot of genes (violet circles) relevant to the separation of PCA. Blue arrows show box plots with expression levels of a few genes (Higd2a, Gm24400, Snora30, Cox5b, Mir5125, Acta2, Map3k7cl, Snora30 and Map3k7cl) that contribute to separation of the two dietary groups (HGD, LGD). In the box plots, the black dots represent gene expression levels, the notch indicates 95% confidence interval around the median of each group, and yellow diamond indicates the average gene expression of each group.

Figure 2 .

 2 Figure 2. Principal Component Analysis (PCA) and loading plot of genes expressed in hippocampal microvessels of the high glycemic diet (HGD) and the low glycemic diet (LGD). (A) PCA scatter plot of the microarray data shows the trends of the expression profiles of the hippocampal microvasculature in the high glycemic diet (HGD, red circles) and low glycemic diet (LGD, green circles), respectively. The PCA plot captures the variance in a dataset in terms of principal components and displays the most significant of these on the x and y axes. The percentages of the total variation that are accounted for by the 1st and 2nd principal components are shown on the x-and y-axes labels. The data are shown for three biological replicates for each dietary group. (B) Loading plot of genes (violet circles) relevant to the separation of PCA. Blue arrows show box plots with expression levels of a few genes (Higd2a, Gm24400, Snora30, Cox5b, Mir5125, Acta2, Map3k7cl, Snora30 and Map3k7cl) that contribute to separation of the two dietary groups (HGD, LGD). In the box plots, the black dots represent gene expression levels, the notch indicates 95% confidence interval around the median of each group, and yellow diamond indicates the average gene expression of each group.

Figure 3 .

 3 Figure 3. Hierarchical clustering of differentially expressed genes in hippocampal microvessels of the low glycemic diet (LGD) and the high glycemic diet (HGD). Upregulated genes are in red and downregulated genes are in blue. The data are shown for three biological replicates for each dietary group.

Figure 3 .

 3 Figure 3. Hierarchical clustering of differentially expressed genes in hippocampal microvessels of the low glycemic diet (LGD) and the high glycemic diet (HGD). Upregulated genes are in red and downregulated genes are in blue. The data are shown for three biological replicates for each dietary group.

Figure 4 Figure 4 .

 44 Figure 4. Cont.

Figure 5 .Figure 5 .

 55 Figure 5. Hierarchical clustering of differentially expressed genes in hippocampal microvessels for the low glycemic diet (LGD) and the high glycemic diet (HGD) with and without soluble epoxide hydrolase inhibitor (sEHI). Hierarchical clustering of differentially expressed gene profiles in hip-

Figure 6 .

 6 Figure 6. Principal Component Analysis (PCA) and loading plot of genes expressed in hippocampal microvessels for the high glycemic diet (HGD) with and without soluble epoxide hydrolase inhibitor (sEHI). (A) PCA scatter plot of the microarray data shows the trends of the expression profiles of the hippocampal microvasculature of the high glycemic diet with (HGD+sEHI, green circles) and without (HGD, red circles) the soluble epoxide hydrolase inhibitor compared to the low glycemic diet (LGD, blue circles). The PCA plot captures the variance in a dataset in terms of principal components and displays the most significant of these on the x and y axes. The percentages of the total variation that are accounted for by the 1st and 2nd principal components are shown on the x-and y-axes labels. The data are shown for three biological replicates for each dietary group. (B) Loading plot of genes (violet circles) relevant to the separation of PCA. Blue arrows show box plots with the expression levels of a few genes (Gm24400, Resp18, Rheb, Brms1l, Ndufa13, Ighv5-12-4) that contribute to separation of the three dietary groups (HGD, HGD+sEHI, LGD). In the box plots, the black dots represent gene expression levels, the notch indicates 95% confidence interval around the median of each group, and the yellow diamond indicates the average gene expression of each group.

Figure 6 .

 6 Figure 6. Principal Component Analysis (PCA) and loading plot of genes expressed in hippocampal microvessels for the high glycemic diet (HGD) with and without soluble epoxide hydrolase inhibitor (sEHI). (A) PCA scatter plot of the microarray data shows the trends of the expression profiles of the hippocampal microvasculature of the high glycemic diet with (HGD+sEHI, green circles) and without (HGD, red circles) the soluble epoxide hydrolase inhibitor compared to the low glycemic diet (LGD, blue circles). The PCA plot captures the variance in a dataset in terms of principal components and displays the most significant of these on the x and y axes. The percentages of the total variation that are accounted for by the 1st and 2nd principal components are shown on the x-and y-axes labels. The data are shown for three biological replicates for each dietary group. (B) Loading plot of genes (violet circles) relevant to the separation of PCA. Blue arrows show box plots with the expression levels of a few genes (Gm24400, Resp18, Rheb, Brms1l, Ndufa13, Ighv5-12-4) that contribute to separation of the three dietary groups (HGD, HGD+sEHI, LGD). In the box plots, the black dots represent gene expression levels, the notch indicates 95% confidence interval around the median of each group, and the yellow diamond indicates the average gene expression of each group.

Figure 7 .

 7 Figure 7. Cont.

Figure 7 .

 7 Figure 7. The effect of the high glycemic diet with and without soluble epoxide hydrolase inhibitor (sEHI) on differential gene expression in the hippocampal microvasculature. A) Heat map showing individual differentially expressed genes (DEGs) in rows, and the two different experimental dietary comparison groups in columns, as follows: column 1: high glycemic diet compared to low glycemic diet (HGD vs LGD); column 2: HGD with soluble epoxide hydrolase inhibitor compared to without inhibitor (HGD+sEHI vs HGD). The data are shown for three biological replicates for each dietary group. Up-regulated genes are shown in red, and down-regulated genes in green. The rectangular box insert shows the list of genes with opposite gene expression patterns in the two comparisons (HGD vs. LGD compared to HGD +sEHI vs. HGD). B) Pearson's correlation analysis of the DEGs between the comparisons, HGD vs. LGD and HGD+sEHI vs. HGD.

Figure 7 .

 7 Figure 7. The effect of the high glycemic diet with and without soluble epoxide hydrolase inhibitor (sEHI) on differential gene expression in the hippocampal microvasculature. (A) Heat map showing individual differentially expressed genes (DEGs) in rows, and the two different experimental dietary comparison groups in columns, as follows: column 1: high glycemic diet compared to low glycemic diet (HGD vs LGD); column 2: HGD with soluble epoxide hydrolase inhibitor compared to without inhibitor (HGD+sEHI vs HGD). The data are shown for three biological replicates for each dietary group. Up-regulated genes are shown in red, and down-regulated genes in green. The rectangular box insert shows the list of genes with opposite gene expression patterns in the two comparisons (HGD vs. LGD compared to HGD+sEHI vs. HGD). (B) Pearson's correlation analysis of the DEGs between the comparisons, HGD vs. LGD and HGD+sEHI vs. HGD.

Figure 8 .

 8 Figure 8. Cont.

1, 13 , 31 Figure 8 .Figure 8 .

 133188 Figure 8. Effect of the high glycemic diet (HGD) with soluble epoxide hydrolase inhibitor (sEHI) on cellular pathways of differentially expressed protein-coding genes, miRNA targets, and LncRNA targets in the hippocampal microvasculature. A) Network of differentially expressed protein-coding genes (grey circles), transcription factors (red hexagons), miRNAs (green diamonds) and their targets (blue circles), LncRNAs (purple rectangles) and their targets (blue circles) in the hippocampal microvessels for the high glycemic diet (HGD) with soluble epoxide hydrolase inhibitor (sEHI) when compared to HGD without sEHI. B) Histogram of subset of relevant pathways of differentially expressed protein-coding genes, miRNA targets, and LncRNA targets in the hippocampal microvasculature for the HGD+sEHI vs HGD. The data are shown for three biological replicates for each dietary group. Statistically significant pathways (p < 0.05) were identified using

Figure 9 .Figure 9 .

 99 Figure 9. Correlation of the high glycemic diet (HGD) with and without inhibitor to Alzheimer's disease. Pearson correlation analysis between the genes modulated by the HGD vs LGD and the gene expression profiles observed in patients with Alzheimer's disease extracted from GEO database number A) GSE118553 and B) GSE132903. Correlation between genes modulated by the

Figure 10 .

 10 Figure 10. Analysis of DEGs with a significant interaction of diet and inhibitor. A) Relevant enriched WikiPathways and KEGG Pathways for DEGs with an interaction between diet and inhibitor effects. B) Signal plots of selected DEGs with an interaction between diet and inhibitor effects. Data shown are from three biological replicates (mean +/-standard deviation). The DEGs displayed are cytochrome c oxidase subunit Vb (Cox5b), NADH dehydrogenase (ubiquinone) 1 alpha subcomplex 2 (Ndufa2), and Rho-associated coiled-coil containing protein kinase 1 (Rock1).

Figure 10 .

 10 Figure 10. Analysis of DEGs with a significant interaction of diet and inhibitor. (A) Relevant enriched WikiPathways and KEGG Pathways for DEGs with an interaction between diet and inhibitor effects. (B) Signal plots of selected DEGs with an interaction between diet and inhibitor effects. Data shown are from three biological replicates (mean +/-standard deviation). The DEGs displayed are cytochrome c oxidase subunit Vb (Cox5b), NADH dehydrogenase (ubiquinone) 1 alpha subcomplex 2 (Ndufa2), and Rho-associated coiled-coil containing protein kinase 1 (Rock1).

Funding:

  This work was supported by an award from the Richard A. and Nora Eccles Foundation (A20-0111), the Richard A. and Nora Eccles Harrison Endowed Chair in Diabetes Research (J.C.R.), and the Frances Lazda Endowed Chair in Women's Cardiovascular Medicine (A.C.V). Laser Capture Microdissection and Cryo-Sectioning were conducted at the Cellular and Molecular Imaging core facility at the University of California Davis (UC Davis) Center for Health and the Environment. The laser capture microscope was funded by NIHS10RR-023555. Microarray Scanning and Hybridization was done by the Genomics Shared Resource at the UC Davis Medical Center. We are also grateful for the technical support and/or services provided to our research by the UC Davis Mouse Metabolic Phenotyping Center which assisted with the metabolic assays supported by U24 DK092993. Institutional Review Board Statement: The institutional review board of the University of California, Davis, the Institutional Animal Care and Use Committee (IACUC) approved this project protocol number 20943 on 18 April 2019. Informed Consent Statement: Not applicable.

Acknowledgments:

We thank Nikita Patel, Anthony Pham, Taarini Hariharan, Nejma Wais, Ryan Vinh, and Corey Buckley for technical assistance in this project.

Data Availability Statement:

We have deposited the microarray data in GEO and the accession number is GSE185057.

Conflicts of Interest:

The authors declare no conflict of interest.