Association of cerebrovascular dysfunction with the development of Alzheimer’s disease-like pathology in OXYS rats
BMC Genomics volume 19, Article number: 75 (2018)
Cerebrovascular dysfunction plays a critical role in the pathogenesis of Alzheimer’s disease (AD): the most common cause of dementia in the elderly. The involvement of neurovasculature disorders in the progression of AD is now increasingly appreciated, but whether they represent initial factors or late-stage pathological changes during the disease is unclear. Using senescence-accelerated OXYS rats, which simulate key characteristics of sporadic AD, we evaluated contributions of cerebrovascular alterations to the disease development. At preclinical, early, and advanced stages of AD-like pathology, in the hippocampus of OXYS and Wistar (control) rats, we evaluated (i) the blood vessel state by histological and electron-microscopic analyses; (ii) differences in gene expression according to RNA sequencing (RNA-Seq) to identify the metabolic processes and pathways associated with blood vessel function; (iii) the amount of vascular endothelial growth factor (VEGF) by western blot and immunohistochemical analysis.
We observed a loss of hippocampal blood vessel density and ultrastructural changes of those blood vessels in OXYS rats at the early stage of AD-like pathology. There were significant alterations in the vessels and downregulation of VEGF with an increased amount of amyloid β1–42 there at the advanced stage of the disease. According to RNA-Seq data analysis, major alterations in cerebrovascular processes of OXYS rats were associated with blood vessel development, circulatory system processes, the VEGF signaling pathway, and vascular smooth muscle contraction. At preclinical and early stages of the AD-like pathology, these processes were upregulated and then downregulated with age. At the advanced stage in OXYS rats, differentially expressed genes (DEGs) were associated with downregulation of cerebrovascular function as compared to Wistar rats. Among the 46 DEGs at the preclinical stage of the disease, 28 DEGs at the early stage, and among 85 DEGs at the advanced stage, using functional analysis and gene network construction, we identified genes (Nos1, P2rx4, Pla2g6, and Bdkrb2) probably playing a significant role in the development of cerebrovascular dysfunction in OXYS rats.
Changes in expression of the genes functionally associated with cerebrovascular processes already in the early period of life may contribute to the development of AD-like pathology in OXYS rats.
Vascular dysfunction is a universal feature of aging that contributes to the cognitive decline and neurodegeneration observed in Alzheimer’s disease (AD): the most common cause of dementia in the elderly [1, 2]. AD is characterized by a gradual decline in memory and cognition that correlate with synaptic dysfunction and neuronal loss. It was suggested over 30 years ago that vascular defects present in AD may be important for the development of the disease . The role of a neurovasculature disorder in the progression of AD is now becoming more fully appreciated. It is clear that microvascular deficits diminish cerebral blood flow and, consequently, the brain’s supply of oxygen, energy, substrates, and nutrients. Moreover, such deficits impair the clearance of neurotoxic molecules, such as amyloid β (Aβ), which are released into the brain interstitial fluid and accumulate in the parenchyma and vessel walls. Patients with AD frequently show focal changes in brain microcirculation. These changes include alterations in density and morphology of cerebral microvasculature, increased endothelial pinocytosis, a decrease in mitochondrial content, accumulation of collagen and perlecans in the basement membrane, loss of tight junctions and/or adherens junctions, and a blood-brain barrier breakdown with leakage of blood-borne molecules [4, 5].
Multiple epidemiological studies have shown a remarkable overlap among risk factors for a cerebrovascular disorder and sporadic, late-onset AD , but whether they represent initial factors or late-stage pathological changes during the disease is still a matter of debate. Growing data from brain imaging studies in humans and animal models suggest that cerebrovascular dysfunction may precede a cognitive decline and the onset of neurodegenerative changes in AD and animal models, yet the molecular and cellular mechanisms and targets that underlie the contribution of vascular disease to AD are not well understood .
Because in humans, the research on the pathogenesis of AD, especially early asymptomatic stages, is problematic, there is a need for animal models that closely mirror the human pathology. There is evidence that a suitable model of sporadic form of AD is senescence-accelerated OXYS rats. These animals spontaneously develop all key signs of AD such as structural neurodegenerative alterations, neuronal loss, synaptic damage, disturbances of the neurotrophic supply, Aβ1–42 peptide accumulation and tau hyperphosphorylation in the hippocampus, and impairment of the learning ability and memory [7,8,9,10].
Using MRI, we showed that at the age of 12 months, when AD was progressing, OXYS rats had structural and functional alterations in the cerebral blood flow typical of chronic ischemia such as reduced cerebral blood flow and a decline of cerebrovascular reactivity in response to a vasodilatory challenge, also defined as the cerebrovascular reserve [11, 12], as in patients with AD. In addition, we found that OXYS rats show altered hemorheological properties of blood resulting from abnormal red blood cell deformability and aggregation that are regarded as one of the mechanisms underlying the complex etiology of AD , impairing the oxygen transport efficiency of blood in AD.
In this study, we assessed the mechanisms of cerebrovascular dysfunction and its contributions to the development of AD-like pathology in OXYS rats. We compared these rats with the control Wistar rats at different stages of the disease in OXYS rats, including the preclinical stage. We used histological and electron-microscopic analyses to evaluate cerebrovascular alterations during development of the signs of AD in OXYS rats, a bioinformatic analysis of RNA-Seq data to identify the metabolic processes and pathways involved in the development of cerebrovascular dysfunction, and quantification of vascular endothelial growth factor (VEGF), which is a neuroprotective cytokine promoting neurogenesis and angiogenesis in the brain.
Structural cerebrovascular alterations in the hippocampus of OXYS rats
To evaluate cerebrovascular alterations during the development of AD-like pathology in OXYS rats, we conducted comparative analysis of the density and ultrastructural state of blood vessels in the hippocampal СА1 and СА3 regions and dentate gyrus in OXYS and Wistar rats at the age of 20 days and 5 and 18 months. We found that total blood vessel density in the hippocampus was lower in OXYS rats (F1,54 = 49.2, p < 0.0001). The maximal density of blood vessels in the hippocampus of OXYS and Wistar rats was observed at the age of 20 days, and this density decreased with age in both rat strains (F2,54 = 370.7, p < 0.0001; Fig. 1a, b). When we analyzed the vessel density in the hippocampal СА1 and СА3 regions and dentate gyrus, we found that in OXYS rats, a significant decrease of this parameter occurred at age 5 months in the dentate gyrus and at the age of 18 months in all hippocampal regions as compared to Wistar rats (p < 0.05; Fig. 1a).
According to electron-microscopic analysis (Fig. 1e–j), in contrast to Wistar rats, 5-month-old OXYS rats showed the first changes in structural organization of the vessels, such as an increase in electron density of the basal membrane of endotheliocytes (Fig. 1f) and hyperchromatic nuclei of endotheliocytes. By the age of 18 months, the degenerative changes in hippocampal vessels of OXYS rats had progressed and were characterized by thickening of the wall and a hyperchromatic endothelial nucleus, the cytoplasm of endotheliocytes contained vacuoles, and the membrane formed folds (Fig. 1h, g). Besides, in the cytoplasm of endotheliocytes, large vesicles formed, probably due to destructive changes in the organelles (mitochondria and the endoplasmic reticulum; Fig. 1j). In addition, the pericytes contained lipid aggregates in the cytoplasm, a large number of lysosomes, and low content of organelles. In the mitochondria of pericytes, we often observed destruction of cristae. Thus, these data indicate a hippocampal blood vessel density loss and ultrastructural changes in those blood vessels in OXYS rats in the period of active development of AD-like pathology, and significant alterations of these vessels at the progressive stage of the disease.
Differential expression of genes in the hippocampus of OXYS rats
To identify the pathways and biological functions involved in the development of cerebrovascular dysfunction in OXYS rats, we used RNA-Seq to obtain the lists of differentially expressed genes (DEGs) in the hippocampus of 20-day-old and 5- and 18-month-old OXYS rats compared to age-matched Wistar rats. Differential expression of genes was evaluated using the DESeq software . We identified 841 protein-coding mRNAs (at p < 0.01) that are differentially expressed in 20-day-old OXYS rats compared with age-matched Wistar rats; among those genes, 281 were downregulated, and 560 were upregulated. There were 489 genes that are differentially expressed in 5-month-old OXYS rats compared with age-matched Wistar rats; among these genes, 209 were downregulated, and 280 were upregulated. At 18 months of age, there were 1324 DEGs; among them, 750 were downregulated, and 574 were upregulated in OXYS rats.
Functional annotation of DEGs associated with cerebrovascular dysfunction in OXYS rats
Differentially expressed protein-coding genes were functionally annotated by means of Gene Ontology (GO) biological process categories and pathways as described in Methods. According to DAVID, at the age of 20 days in OXYS rats, seven GO terms (p < 0.05) were found to be associated with blood vessel processes; among these genes, four were upregulated (Additional file 1). In particular, GO terms were enriched for functional categories explicitly related to the blood circulation and blood pressure. At 5 months of age, there were two upregulated GO terms related to the blood vessel development. At 18 months of age, there were nine GO terms; among them, five were downregulated and were enriched for genes belonging to processes involved in blood vessel development and hemopoiesis (Additional file 1). In summary, there are 19 DEGs involved in biological processes associated with hippocampal blood vessels in 20-day-old OXYS rats; among these genes, 3 were downregulated, and 16 were upregulated. At the age of 5 months, there were 8 upregulated DEGs in OXYS rats associated with blood vessel processes. At the age of 18 months, among the 65 DEGs associated with blood vessel processes, mRNA expression of 26 genes was increased in OXYS rats and decreased for the other 39 genes, in comparison with Wistar rats (Additional file 2). Interstrain differences in the expression levels of genes Cxcl12 and Itga7 were detected at 5 and 18 months, and for the 7 common genes, interstrain differences were observed between ages 20 days and 18 months (Agt, Bdkrb2, Gaa, Nos1, P2rx4, Pcsk5, and Ptgs1), the expression of all genes increased, except for gene Bdkrb2.
In KEGG pathway analysis, we found 3 pathways (p < 0.05) associated with blood vessel processes in OXYS rats: the vascular smooth muscle contraction pathway was upregulated in OXYS rats at the age of 20 days and 5 months; the VEGF signaling pathway was upregulated at the age of 20 days and downregulated at the age of 18 months; the hematopoietic cell lineage pathway was downregulated at the age of 5 and 18 months (Additional file 2). Note that interstrain differences in the high expression level of gene Pla2g6, which is involved in the VEGF signaling pathway and vascular smooth muscle contraction, were observed at all 3 ages.
Age-associated transcriptome changes that are involved in cerebrovascular function in OXYS and Wistar rats
We found that in both rat strains, the genes with significantly changed expression with age were associated with a sufficiently large number of GO terms (p < 0.05) that are involved in blood vessel processes according to DAVID (Additional file 1). Among these genes, in OXYS and Wistar rats, the genes significantly downregulated between the ages of 20 days and 5 months as well as between the ages of 5 and 18 months were associated with blood vessel development, blood circulation, and its regulation (Fig. 2a). As for the KEGG pathway analysis, we found that in Wistar rats with age, the DEGs were associated with a change in vascular smooth muscle contraction, in the VEGF signaling pathway, and in the hematopoietic cell lineage pathway (Additional file 1). Whereas in OXYS rats, the vascular smooth muscle contraction pathway was significantly downregulated between the ages of 5 and 18 months, and the VEGF signaling pathway was downregulated in both age periods (Fig. 2a; Additional file 1). For a combined analysis, we obtained a total list of DEGs involved in GO terms according to DAVID and in pathways according to KEGG in OXYS and Wistar rats in both age periods (Additional file 1). We found, that in OXYS rats, 135 genes significantly changed expression between the ages of 20 days and 5 months (of these, 84 were downregulated, and 51 were upregulated), and 197 genes changed expression between the ages of 5 and 18 months (of these, 147 were downregulated, and 50 were upregulated; Fig. 2b). In Wistar rats, 150 genes significantly changed expression between the ages of 20 days and 5 months (of these, 92 were downregulated, and 58 were upregulated), and 190 genes changed expression levels between the ages of 5 and 18 months (of these, 114 were downregulated, and 76 were upregulated; Fig. 2b). Figure 2c shows a Venn diagram visualizing the overlaps among the lists of DEGs associated with blood vessel processes in OXYS and Wistar rats in the two periods. In OXYS rats, between the ages of 20 days and 5 months, 15 genes exclusively changing expression were mainly related to blood circulation and regulation of blood pressure (Npy, Stat1, Cdc42, Tac1, Apln, and Cyb11b2; Fig. 2c; Additional file 2). In Wistar rats during the same period, there were 17 genes that exclusively changed expression; they were mainly related to blood vessel morphogenesis and regulation of angiogenesis (Mylk, Prkx, Epha2, Adrb2, Serpine1, Notch1, Dll4, and Vegfc). In OXYS rats, between the ages of 5 and 18 months, 28 genes exclusively changing expression levels were mainly related to blood circulation and regulation of blood vessel size (Agt, Htr2a, Smad5, Gclm, P2rx4, Npy1r, Cacna1b, and Angpt1). By contrast, in Wistar rats during the same period, there were 32 genes that exclusively changed expression; they were mainly related to the regulation of lymphocyte proliferation and of immune responses (Il2ra, Il15, Cd38, Cd74, and Serpine1). Only in OXYS rats, in both age groups, changed expression level were observed for 5 genes, 3 of them are associated with blood vessel morphogenesis (Cited1, Plxnd, and Itgav). In Wistar rats, among 8 DEGs in both age groups, the change of expression of 4 genes was also associated with blood vessel morphogenesis (Gja5, Ctgf, Apold1, and Itga7; Fig. 2c; Additional file 2).
Major altered cerebrovascular processes in OXYS rats and construction of a gene network
Using the DAVID and KEGG analyses, we found that in 20-day-old and 5-month-old OXYS rats, cerebrovascular processes and pathways are upregulated (Fig. 2a). The genes that are differentially expressed during aging were downregulated, and the results revealed that at the age of 18 months in OXYS rats, DEGs were associated with downregulation of cerebrovascular function as compared with Wistar rats. Among these processes, the most significant alterations are associated with blood vessel development, circulatory system processes, and pathways associated with VEGF signaling, vascular smooth muscle contraction, and hematopoietic cell lineages (Fig. 2a). We next compiled a total list of DEGs involved in GO processes and pathways in OXYS rats in 3 age groups (Additional file 1). Already at the age of 20 days in OXYS rats, there are 46 DEGs (16 downregulated and 30 upregulated genes). At the age of 5 months, there were 28 DEGs (10 downregulated and 18 upregulated genes) in OXYS rats associated with blood vessel processes. At the age of 18 months, of the 85 DEGs associated with blood vessel processes, mRNA expression of 56 genes was decreased in OXYS rats and increased for the other 29 genes, in comparison with Wistar rats (Fig. 3a; Additional file 1). Interstrain differences in the expression levels of 13 genes were detected at the ages of 20 days and 5 and 18 months. It should be noted that mRNA levels of Cxcl12, Nhej1, Nos1, P2rx4, Pcsk5, Pla2g6, and Ptpn22 increased, and the expression of Cd4, Cd7, Hhip, Lig4, and Spta1 decreased in OXYS rats at all 3 ages. Only expression of the Bdkrb2 gene decreased in 20-day-old and 18-month-old OXYS rats and increased in 5-month-old rats. In addition, interstrain differences in the expression levels of upregulated genes Itgam and Pla2g2d were observed at 20 days and 5 months, and for the upregulated gene Itga7 and downregulated genes C3 and Zcchc2 interstrain differences in the expression levels were detected at ages 5 and 18 months. For the 7 common genes that were found between ages 20 days and 18 months, the expression of Fech, Gaa, Igfbp5, and Ptgs1 increased and the expression of Cav2, Dclre1c, and Itgav diminished.
We next used GeneMANIA to construct the interaction networks for the analysis of transcriptomic changes associated with functional categories of cerebrovascular dysfunction during the AD-like pathology development (Fig. 3c–e). It is evident that most of the genes are directly or indirectly related to one another. The GeneMANIA results on enrichment with GO terms for the network’s members are presented in Additional file 3. Among the most enriched GO terms were the regulation of blood pressure (p < 4.85E–10), blood circulation (p < 1.48E–9), vasodilation (p < 4.67E–7), and regulation of blood vessel size (p < 1.43E–6) in the gene set of 20-day-old rats (Fig. 3c); blood circulation (p < 2.85E–4), the regulation of blood pressure (p < 2.85E–4), muscle contraction (p < 1.08E–5), and blood vessel morphogenesis (p < 5.98E–4) in the gene set of 5-month-old rats (Fig. 3d); and blood vessel morphogenesis (p < 2.22E–17), blood circulation (p < 2.16E–11), muscle contraction (p < 9.49E–13), regulation of blood vessel size (p < 1.24E–11), and a vascular process in circulatory system (p < 4.06E–12) in the gene set of 18-month-old OXYS rats (Fig. 3e). Analysis of the gene networks in 20-day-old and 5- and 18-month-old animals revealed that the common genes probably play a significant role in the regulation of blood vessel processes in the hippocampus of OXYS rats during development of AD-like pathology: Nos1, P2rx4, Pla2g6, and Bdkrb2. Taken together, these results indicate that in OXYS rats, the development of AD signs is preceded by changes in the expression of hippocampal genes involved in cerebrovascular processes; these changes increased with the progression of disease pathology.
An age-dependent decrease of the VEGF level is accompanied by accumulation of Aβ1–42 in the hippocampal blood vessels of OXYS rats
On the basis of the comparative analysis of the functional annotation and pathways, we next examined the concentration of VEGF in the hippocampus of OXYS and Wistar rats at the ages of 20 days and 5 and 18 months. According to western blot analysis (Fig. 4a), the VEGF level in the hippocampus was lower in OXYS rats (F1,36 = 5.98, p < 0.02). The maximal level of VEGF in the hippocampus of OXYS and Wistar rats was observed at the age of 20 days, and this concentration decreased with age in both rat strains (F2,36 = 26.6, p < 0.0001), but in OXYS rats, its level at the age of 18 months was lower than that in Wistar rats (p < 0.04). Immunolabeling of tissue samples with an anti-VEGF antibody confirmed this result (Fig. 4b). Moreover, immunolabeling for Aβ1–42 revealed accumulation of Aβ1–42 with age in the hippocampal vessels of OXYS rats and increased colocalization of VEGF with Aβ1–42 at the age of 18 months (Fig. 4b).
There is growing evidence of microvascular pathophysiological alterations having a causal role in the development of AD . Our results indicate that cerebrovascular alterations also make a significant contribution to the development of AD-like pathology in OXYS rats. A decline of the hippocampal blood vessel density and ultrastructural abnormalities of these vessels in OXYS rats were revealed in the period of active manifestation of the signs of AD-like pathology, as were significant alterations of these vessels at the progressive stage of the disease, which imply reduced blood flow. Such alteration results in a mismatch between the supply and demand of oxygen and metabolic substrates during the functioning of cerebral tissue; this situation leads to neuronal dysfunction . It should be noted that the signs pointing to a decrease in blood flow were also present in the hippocampus of 5-month-old OXYS rats prior to Aβ deposition, which appears at the age of approximately a year  against the background of chronic ischemia and of a decline of cerebrovascular reactivity [11, 12]. A gradual decline in oxygen as well as glucose supply to the brain during aging or hypoxia manifested themselves as contributing factors to hypometabolism. Progression of cerebrovascular alterations impairs vascular clearance of Aβ from the brain and contributes directly to the cognitive symptoms and thereby accelerates the progression of AD . Here we demonstrated that significant cerebrovascular disturbances were observed in the hippocampus of OXYS rats at the age of 18 months, when the cognitive impairments and neurodegenerative processes in the brain were well pronounced and accompanied by significant accumulation of toxic forms of Аβ, formation of amyloid plaques, hyperphosphorylation of the tau protein, mitochondrial dysfunction, and a deficit of synapses [7,8,9].
Next, we focused on bioinformatic analyses of the data from whole-transcriptome profiling of the hippocampus to identify the molecular mechanism, the metabolic pathways, and biological functions involved in the development of the cerebrovascular dysfunction in OXYS rats. Aging-associated changes occurring at the molecular level contribute to cellular dysfunction and disease. It is important to determine the molecular background of the disease at preclinical stages, which cannot be studied in humans. Our study includes a comparative RNA-Seq analysis of the hippocampal transcriptome between OXYS and Wistar rats at the age of 20 days, which is the period of completion of postnatal development of the brain. By monitoring hippocampal changes of transcript levels at several time points throughout the development of AD-like pathology in OXYS rats, we revealed that already at the age of 20 days, OXYS rats have substantial changes in the expression of 841 protein-coding genes. These early changes before the first clinical signs of AD can serve as the foundation for the development and progression of AD at a later age. It is important that the overwhelming majority of these genes (560) was downregulated in OXYS rats in comparison with age-matched Wistar rats. This trend changed at the age of 5 months when the number of DEGs decreased to 489, and 280 of them were upregulated and were also detected at 18 months of age, when among 1324 DEGs, 750 genes were downregulated, and 574 genes were upregulated in OXYS rats.
In accordance with the purpose of the present study, we analyzed a pool of DEGs that encode gene products that may contribute to the cerebrovascular dysfunction in OXYS rats. In the hippocampus of OXYS and Wistar rats, we identified aging-related DEGs associated with a sufficiently large number of GO terms that are involved in blood vessel processes including blood vessel development, blood circulation, and its regulation. Moreover, the combined analysis of the total list of these DEGs involved in GO terms according to DAVID and in pathways according to KEGG in OXYS and Wistar rats identified similar numbers of DEGs associated with vessel processes. Nonetheless, these processes were represented by different sets of cerebrovascular functions and genes, thereby highlighting a qualitative difference in the age-related alterations of the cerebral vessels between healthy Wistar rats and OXYS rats developing the AD-like pathology. Although between the ages of 20 days and 5 months, genes that exclusively changed expression in OXYS rats were mainly related to blood circulation and regulation of blood pressure, in Wistar rats, such DEGs were related to blood vessel morphogenesis and regulation of angiogenesis. Between the ages of 5 and 18 months, genes exclusively changing expression levels in OXYS rats were mainly related to blood circulation and regulation of blood vessel size, whereas in Wistar rats during the same period, genes that exclusively changed expression were mainly related to the regulation of lymphocyte proliferation and regulation of immune responses. It is noteworthy that 22 of the 32 DEGs associated with the immune system were upregulated. Aging is characterized by chronic low-grade inflammation. Activation of inflammatory reactions could be both beneficial and detrimental to the brain and depends on the magnitude of neurodegenerative changes. In contrast, the imbalance of inflammatory responses in the central nervous system may serve as an initiating factor in many neurodegenerative diseases, including AD .
The main aim of this study was to identify the metabolic processes and pathways involved in the development of cerebrovascular dysfunction in OXYS rats. We found that at the preclinical stage of development of the AD-like pathology (age 20 days) as well in the period of its active manifestation (age 5 months) in the hippocampus of OXYS rats, cerebrovascular processes and pathways are upregulated in comparison with Wistar rats; this phenomenon can be considered a compensatory reaction that prevents the development of neurodegenerative processes. These results are consistent with our earlier findings that the neuronal populations in all regions of the hippocampus of OXYS rats increase at the age of 5 months as compared with age-matched Wistar rats owing to the enhanced neurogenesis, which we consider an endogenous compensatory mechanism of repair . During the cerebrovascular dysfunction, elevated neurogenesis associated with increased production of immature neurons failed to integrate the existing networks. Numerous studies have suggested that the dysfunctional neurogenesis in the brain leads subsequently to subtle early manifestations of a disease, which could, in turn, render neurons more vulnerable to AD and contribute to memory impairment .
In the hippocampus of OXYS rats, the genes that are differentially expressed during aging were downregulated, and as a result, at the age of 18 months, DEGs were associated with downregulation of cerebrovascular function as compared with Wistar rats. Among these processes, the most significant alterations are associated with blood vessel development, circulatory system processes, and VEGF pathways. VEGF, originally described as a key angiogenic factor, has recently been shown to perform multiple functions in a developing and adult nervous system by acting both on blood vessels and neurons . VEGF has also been shown to serve as a neuroprotective factor in AD. Relative to controls, patients with AD have lower levels of serum VEGF and cerebral capillary VEGF expression in the hippocampus .
It is noteworthy that the VEGF signaling pathway in OXYS rats was upregulated at the age of 20 days during completion of postnatal brain development and was downregulated at the age of 18 months: in the period of active progression of the AD-like pathology. The maximal concentration of VEGF in the hippocampus of both OXYS and Wistar rats was detected at the age of 20 days, and this concentration declined with age, becoming lower in OXYS rats than in Wistar rats at the age of 18 months. VEGF has high specificity for endothelial cells and is profoundly accumulated in (and colocalized with) amyloid plaques in the brain of AD patients . We observed increased colocalization of VEGF with Aβ1–42 in the hippocampal blood vessels of OXYS rats at age 18 months.
Using network analysis, we identified a core set of 4 interrelated DEGs from functional categories of cerebrovascular dysfunction that were common to 20-day-old and 5- and 18-month-old OXYS rats and probably play a significant role in the development of AD-like pathology in OXYS rats. Three of them (Nos1, P2rx4, Pla2g6) were upregulated at all 3 ages and 1 (Bdkrb2) downregulated at the age of 20 days and 18 months, and upregulated at the age of 5 months.
The neuronal nitric oxide synthase NOS1 family, a group of NO synthases that synthesize NO from L-arginine, is a pivotal mediator of neurotransmission, spatial learning, and cognition and has been repeatedly implicated in the neurotoxicity associated with neurodegenerative diseases . A similar NO synthase (endothelial NOS; eNOS), NOS1 plays an obligatory role in the regulation of cerebral blood flow and cell viability and in the protection of nerve cells or fibers from pathogenic factors associated with AD . All three NOS isoforms show aberrant patterns of expression in AD, altering intracellular signaling and routing oxidative stress in directions that are uncompounded ; however, data indicating changes in the NOS1 level and their involvement in amplifying or reducing neuronal damage are inconsistent, and to some degree, contradictory. Some studies on different models show overexpression and/or increased activity of nNOS during AD progression [27, 28], whereas others describe reduced nNOS expression and/or activity [29, 30].
The product of the P2rx4 gene belongs to the family of purinoceptors for ATP. Accumulating evidence reveals major participation of P2 receptors in AD pathology, including Aβ production and elimination, neuroinflammation, neuronal function, and cerebral blood flow . Receptor P2RX4 functions as a ligand-gated ion channel with high calcium permeability and performs critical functions in synaptic plasticity and neuronal survival. Aβ1–42 promotes accumulation of calcium-permeable purinergic receptor P2X4 in neurons and thus induces synaptic dysfunction and neuronal death . Experimental evidence from human AD tissue and studies involving selective purinergic receptor agonists and antagonists in vitro and in AD animal models have revealed that purinergic receptors represent novel therapeutic targets in AD .
The protein encoded by the Pla2g6 gene is an A2 phospholipase, a class of enzymes that is known to play multiple roles in the membrane phospholipid homeostasis and in the production of a variety of lipid mediators . Pla2g6 is the causative gene for a familial form of juvenile-onset dystonia parkinsonism and is associated with the pathogenesis of idiopathic Parkinson’s disease . Reduced phospholipase A2 (PLA2) activity has been reported in the brain and platelets of patients with AD . In the brain, inhibition of PLA2 suppresses the physiological secretion of APP: a mechanism that increases Аβ formation .
Bdkrb2 encodes a receptor for bradykinin (B2R) that elicits many responses including vasodilation, edema, smooth muscle spasm, and pain fiber stimulation . B2R activation may be intimately involved in Aβ-induced neuroinflammation: a key contributor to AD progression. Blockage of B2R can protect against the memory deficits induced by Aβ in mice, pointing to a neuroprotective role of B2R . Here, Bdkrb2 expression was found to be reduced at preclinical and advanced stages, and was increased at early stage of the AD-like pathology in OXYS rats.
AD is a multifactorial disease and involves several etiopathogenic mechanisms including cerebrovascular dysfunction. Understanding how neurovascular impairments can lead to AD-related neurodegeneration is necessary to develop effective drugs that prevent or reverse these changes. On the basis of our present and recent results, we can conclude that the development of cerebrovascular impairments in OXYS rats with age along with the age-related decrease in the density of cerebral blood vessels and their structural alterations including accumulation of Aβ, impairment of cerebral blood flow, and increased neuronal degeneration and susceptibility to hypoxia and ischemia. Moreover, we can conclude that changes in the cerebrovascular processes already in the early period of life may contribute in the future to the accelerated cognitive decline in the elderly, and become a factor risk of AD.
Animals and treatments
All the experimental procedures were in compliance with the European Communities Council Directive of 24 November 1986 (86/609/EEC). The protocol of the animal study was approved by the Commission on Bioethics of the Institute of Cytology and Genetics, Novosibirsk, Russia.
To assess the age-dependent transcriptomic changes, male senescence-accelerated OXYS rats (n = 12) and age-matched male Wistar rats (n = 12) at the age of 20 days and 5 and 18 months were obtained from the Center for Genetic Resources of Laboratory Animals at the Institute of Cytology and Genetics, the Siberian Branch of the Russian Academy of Sciences (RFMEFI61914X0005 and RFMEFI61914X0010). The OXYS strain was derived from the Wistar strain of rats at the Institute of Cytology and Genetics as described earlier  and was registered in the Rat Genome Database (http://rgd.mcw.edu/). At this point, we have the 109th generation of OXYS rats, with spontaneously developing cataract and accelerated senescence syndrome inherited in a linked manner. At the age of 4 weeks, the pups were weaned, housed in groups of 5 animals per cage (57 × 36 × 20 cm), and kept under standard laboratory conditions (22 °C ± 2 °C, 60% relative humidity, and 12-h light/12-h dark cycle; lights on at 9 a.m.). The animals were provided with standard rodent feed (PK-120-1; Laboratorsnab, Ltd., Moscow, Russia) and water ad libitum.
For RNA-Seq, 12 male OXYS rats and 12 age-matched male Wistar rats at ages 20 days (n = 3), 5 months (n = 3), and 18 months (n = 3) were euthanized by CO2 inhalation. After decapitation, the hippocampus was excised rapidly, placed in RNAlater (Ambion, cat. # AM7020), frozen, and stored at − 20 °C prior to analysis. Frozen rat tissues were lysed with the TRIzol Reagent (Invitrogen, Cat. #15596–018), and total RNA was isolated according to the manufacturer’s protocol. RNA quality and quantity were assessed using an Agilent Bioanalyser (Agilent).
More than 40 million single-end reads 50 bp long were obtained for each sample of hippocampal RNA, using Illumina nonstranded sequencing on an Illumina GAIIx instrument at the Genoanalitika Lab, Moscow  in accordance with standard Illumina protocols (mRNA-Seq Sample Prep Kit, cat. # 1004816). Briefly, polyadenylated mRNA was purified from total RNA using Sera-Mag Magnetic Oligo (dT) beads and then broken into small fragments by means of divalent cations and heating. Using a reverse transcriptase and random primers, we synthesized first- and second-strand cDNAs. The cDNA was processed in an end repair reaction with T4 DNA polymerase and Klenow DNA polymerase to blunt the termini. An “A” base was then added to the 3′ end of the blunt phosphorylated DNA fragments, and an Illumina adaptor with a single T overhang at its 3′ end was then ligated to the end of the DNA fragment, for hybridization in a single-read flow cell. After that, a size range of cDNA templates was selected, and these fragments were amplified on a cluster station using the Single-Read Cluster Generation Kit v2. Sequencing-by-synthesis (SBS) at 50-nucleotide length was performed by means of SBS v4 reagents on a Genome Analyzer IIx running the SCS2.8 software (Illumina, cat. #FC-940–4001).
Gene expression analysis
The sequencing data were preprocessed using the Cutadapt tool  to remove adapters and low-quality sequences. The resulting reads were mapped onto the Rnor_5.0 reference genome assembly in the TopHat2 software . The data were then converted into gene count tables using ENSEMBL and RefSeq gene annotations data. The resulting tables were subjected to the analysis of differential gene expression in the DESeq software . The genes with p value < 0.01 were selected as differentially expressed.
To identify the GO terms associated with blood vessel processes overrepresented in a DEG list, the detected DEGs were subjected to functional enrichment analyses by means of the DAVID (Database for Annotation, Visualization and Integrated Discovery, https://david-d.ncifcrf.gov/tools.jsp) tool. Pathway analysis of the DEGs associated with blood vessel processes was conducted using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways (http://www.genome.jp/kegg/).
Construction and analysis of gene interaction networks
The gene interaction networks associated with blood vessel processes were identified on the GeneMANIA web server (http://www.genemania.org/) with default parameters. To construct gene interaction networks associated with blood vessel processes from lists of DEGs (raw p value < 0.01), we selected genes that were in the following GO categories: blood vessel development, vasculature development, blood vessel morphogenesis, patterning of blood vessels, angiogenesis, blood vessel remodeling, endothelial cell development, vasculogenesis, artery morphogenesis, hemopoiesis, hemopoietic or lymphoid organ development, blood circulation, circulatory system process, vascular process in circulatory system, vasoconstriction, smooth muscle contraction, regulation of blood pressure, positive regulation of smooth muscle cell proliferation, regulation of angiogenesis, regulation of endothelial cell proliferation, regulation of vasoconstriction, positive regulation of endothelial cell proliferation, vascular endothelial growth factor receptor signaling pathway, regulation of smooth muscle cell proliferation according to DAVID and pathways: vascular smooth muscle contraction, VEGF signaling pathway, and hematopoietic cell lineage according to KEGG.
For this analysis, the hippocampus of 20-day-old (n = 6 to 8), 5-month-old (n = 6 to 8), and 18-month-old OXYS and Wistar rats (n = 12 to 16) was separated from the brain, placed in a microcentrifuge tube for protein isolation, and frozen in liquid nitrogen. Immunoblotting was carried out as previously described . Antibodies and dilutions used in this study were as follows: anti-VEGF (1:1000; # ab14078, Abcam) and anti-β-actin (1:1000; # ab1801, Abcam) and respective secondary antibodies (1:5000; ## ab6721, ab97135, Abcam). Quantitative densitometric analyses were performed on digitized images of immunoblots in the ImageJ software (NIH, USA).
The tissue preparation was conducted as described elsewhere . To estimate numeric density of the blood vessels in the middle molecular layer of the dentate gyrus, and in the CA1 and CA3 pyramidal layers of hippocampus, a set of 4–6 serial sections from each animal was used. A 40× objective lens (Axioplan 2, Zeiss) was used to determine the numeric density of capillaries per visual field.
This procedure was performed by a standard indirect method as described previously . Primary antibodies and dilutions were as follows: anti-VEGF (1:250; # ab14078, Abcam) and anti-amyloid-β1–42 (1:250; # ab10148, Abcam). After incubation with the respective secondary antibodies conjugated with Alexa Fluor 555 or Alexa Fluor 488 (## ab150170, ab150073, Abcam) diluted 1:250, the tissue slices were coverslipped with the Fluoro-shield mounting medium containing 4′,6-diamidino-2-phenylindole (DAPI; # ab104139, Abcam) and examined under an Axioplan 2 microscope (Zeiss).
For this purpose, the hippocampus samples (n = 4) were fixed with 2.5% glutaraldehyde in sodium cacodylate buffer and were prepared as described elsewhere . For evaluation of morphological features, capillaries in the pyramidal layer of the CA1 region were identified in the electron micrographs (30 photos per animal). Next, all the morphological changes that are located in the vessels were painted using software. The photos were processed in Adobe Photoshop. Ultrathin sections were studied and photographed by means of an electron microscope JEM-7A.
The data were subjected to ANOVA (Statistica 8.0 software). Two-way ANOVA was used to evaluate the age-dependent effects (age × genotype [strain]). The Newman–Keuls test was applied to significant main effects and interactions to assess the differences between some sets of means. The data were presented as mean ± SEM. The differences were considered statistically significant at p < 0.05.
Database for annotation, visualization, and integrated discovery
Differentially expressed gene
Kyoto encyclopedia of genes and genomes
Log2 fold change
Magnetic resonance imaging
high-throughput RNA sequencing
Vascular endothelial growth factor
Charidimou A, Boulouis G, Gurol ME, Ayata C, Bacskai BJ, Frosch MP, Viswanathan A, Greenberg SM. Emerging concepts in sporadic cerebral amyloid angiopathy. Brain. 2017; https://doi.org/10.1093/brain/awx047.
Toth P, Tarantini S, Csiszar A, Ungvari Z. Functional vascular contributions to cognitive impairment and dementia: mechanisms and consequences of cerebral autoregulatory dysfunction, endothelial impairment, and neurovascular uncoupling in aging. Am J Physiol Heart Circ Physiol. 2017;312(1):H1–H20. https://doi.org/10.1152/ajpheart.00581.2016.
Scheibel AB, Duong TH, Jacobs R. Alzheimer’s disease as a capillary dementia. Ann Med. 1989;21:103–7.
Zlokovic BV. Neurovascular pathways to neurodegeneration in Alzheimer's disease and other disorders. Nat Rev Neurosci. 2011;12(12):723–38. https://doi.org/10.1038/nrn3114.
Tarantini S, Tran CHT, Gordon GR, Ungvari Z, Csiszar A. Impaired neurovascular coupling in aging and Alzheimer's disease: contribution of astrocyte dysfunction and endothelial impairment to cognitive decline. Exp Gerontol. 2017;94:52–8. https://doi.org/10.1016/j.exger.2016.11.004.
de la Torre JC. Vascular risk factor detection and control may prevent Alzheimer's disease. Ageing Res Rev. 2010;9(3):218–25. https://doi.org/10.1016/j.arr.2010.04.002.
Stefanova NA, Kozhevnikova OS, Vitovtov AO, Maksimova KY, Logvinov SV, Rudnitskaya EA, Korbolina EE, Muraleva NA, Kolosova NG. Senescence-accelerated OXYS rats: a model of age-related cognitive decline with relevance to abnormalities in Alzheimer disease. Cell Cycle. 2014;13:898–09.
Stefanova NA, Muraleva NA, Korbolina EE, Kiseleva E, Maksimova KY, Kolosova NG. Amyloid accumulation is a late event in sporadic Alzheimer's disease-like pathology in nontransgenic rats. Oncotarget. 2015;6:1396–13.
Stefanova NA, Muraleva NA, Maksimova KY, Rudnitskaya EA, Kiseleva E, Telegina DV, Kolosova NG. An antioxidant specifically targeting mitochondria delays progression of Alzheimer's disease-like pathology. Aging (Albany NY). 2016;8(11):2713–33. https://doi.org/10.18632/aging.101054.
Rudnitskaya EA, Kolosova NG, Stefanova NA. Impact of changes in Neurotrophic supplementation on development of Alzheimer's disease-like pathology in OXYS rats. Biochemistry (Mosc). 2017;82(3):318–29. https://doi.org/10.1134/S0006297917030105.
Agafonova IG, Kolosova NG, Mishchenko NP, Chaikina EL, Stonik VA. Effect of histochrome on brain vessels and research and exploratory activity of senescence-accelerated OXYS rats. Bull Exp Biol Med. 2007;143:467–71.
Agafonova IG, Kotelnikov VN, Kolosova NG, Stonik VA. Comparative study on hypertension-induced cerebral vascular alterations in two rat lines. Appl Magn Reson. 2012;42:487–97.
Chang CY, Liang HJ, Chow SY, Chen SM, Liu DZ. Hemorheological mechanisms in Alzheimer's disease. Microcirculation. 2007;14(6):627–34.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol 2010;11:R106; PMID:20979621; https://doi.org/10.1186/gb-2010-11-10-r106
Snyder HM, Corriveau RA, Craft S, Faber JE, Greenberg SM, Knopman D, Lamb BT, Montine TJ, Nedergaard M, Schaffer CB, Schneider JA, Wellington C, Wilcock DM, Zipfel GJ, Zlokovic B, Bain LJ, Bosetti F, Galis ZS, Koroshetz W, Carrillo MC. Vascular contributions to cognitive impairment and dementia including Alzheimer's disease. Alzheimers Dement. 2015;11:710–7.
Jessen SB, Brazhe A, Lind BL, Mathiesen C, Thomsen K, Jensen K, Lauritzen M. GABAA receptor-mediated bidirectional control of synaptic activity, intracellular Ca2+, cerebral blood flow, and oxygen consumption in mouse Somatosensory cortex in vivo. Cereb Cortex. 2015;25(9):2594–609.
Stefanova NA, Muraleva NA, Skulachev VP, Kolosova NG. Alzheimer's disease-like pathology in senescence-accelerated OXYS rats can be partially retarded with mitochondria-targeted antioxidant SkQ1. J Alzheimers Dis. 2014;38(3):681–94.
Sagare AP, Bell RD, Zlokovic BV. Neurovascular dysfunction and faulty amyloid β-peptide clearance in Alzheimer disease. Cold Spring Harb Perspect Med. 2012;2(10) https://doi.org/10.1101/cshperspect.a011452.
Sochocka M, Diniz BS, Leszek J. Inflammatory response in the CNS: friend or foe? Mol Neurobiol. 2016. [Epub ahead of print]. PubMed PMID:27889895.
Lazarov O, Hollands C. Hippocampal neurogenesis: learning to remember. Prog Neurobiol. 2016;138-140:1–18. doi: https://doi.org/10.1016/j.pneurobio.2015.12.006.
Mackenzie F, Ruhrberg C. Diverse roles for VEGF-A in the nervous system. Development. 2012;139(8):1371–80.
Hohman TJ, Bell SP. Jefferson AL; Alzheimer’s Disease Neuroimaging Initiative. The role of vascular endothelial growth factor in neurodegeneration and cognitive decline: exploring interactions with biomarkers of Alzheimer disease. JAMA Neurol. 2015;72(5):520–9.
Yang SP, Bae DG, Kang HJ, Gwag BJ, Gho YS, Chae CB. Co-accumulation of vascular endothelial growth factor with beta-amyloid in the brain of patients with Alzheimer's disease. Neurobiol Aging. 2004;25(3):283–90.
Hardingham N, Dachtler J, Fox K. The role of nitric oxide in pre-synaptic plasticity and homeostasis. Front Cell Neurosci. 2013;7:190. Review.. Asiimwe N, Yeo SG, Kim MS, Jung J, Jeong NY. Nitric Oxide: Exploring the Contextual Link with Alzheimer's Disease. Oxid Med Cell Longev. 2016;2016:7205747
Toda N, Ayajiki K, Okamura T. Cerebral blood flow regulation by nitric oxide in neurological disorders. Can J Physiol Pharmacol. 2009;87(8):581–94.
Hannibal L. Nitric oxide homeostasis in neurodegenerative diseases. Curr Alzheimer Res. 2016;13(2):135–49.
Lüth HJ, Holzer M, Gertz HJ, Arendt T. Aberrant expression of nNOS in pyramidal neurons in Alzheimer's disease is highly co-localized with p21ras and p16INK4a. Brain Res. 2000;852(1):45–55.
Limón ID, Díaz A, Mendieta L, Chamorro G, Espinosa B, Zenteno E, Guevara J. Amyloid-beta (25-35) impairs memory and increases NO in the temporal cortex of rats. Neurosci Res. 2009;63(2):129–37.
Norris PJ, Faull RL, Emson PC. Neuronal nitric oxide synthase (nNOS) mRNA expression and NADPH-diaphorase staining in the frontal cortex, visual cortex and hippocampus of control and Alzheimer's disease brains. Brain Res Mol Brain Res. 1996;41(1–2):36–49.
Hartlage-Rübsamen M, Apelt J, Schliebs R. Fibrillary beta-amyloid deposits are closely associated with atrophic nitric oxide synthase (NOS)-expressing neurons but do not upregulate the inducible NOS in transgenic Tg2576 mouse brain with Alzheimer pathology. Neurosci Lett. 2001;302(2–3):73–6.
Erb L, Cao C, Ajit D, Weisman GA. P2Y receptors in Alzheimer's disease. Biol Cell. 2015;107(1):1–21.
Varma R, Chai Y, Troncoso J, Gu J, Xing H, Stojilkovic SS, Mattson MP, Haughey NJ. Amyloid-beta induces a caspase-mediated cleavage of P2X4 to promote purinotoxicity. NeuroMolecular Med. 2009;11(2):63–75.
Woods LT, Ajit D, Camden JM, Erb L, Weisman GA. Purinergic receptors as potential therapeutic targets in Alzheimer's disease. Neuropharmacology. 2016;104:169–79.
Tomiyama H, Yoshino H, Hattori N. Analysis of Pla2g6 in patients with frontotemporal type of dementia. Parkinsonism Relat Disord. 2011;17(6):493–4.
Miki Y, Tanji K, Mori F, Kakita A, Takahashi H, Wakabayashi K. PLA2G6 accumulates in Lewy bodies in PARK14 and idiopathic Parkinson's disease. Neurosci Lett. 2017;645:40–5.
Gattaz WF, Talib LL, Schaeffer EL, Diniz BS, Forlenza OV. Low platelet iPLA2 activity predicts conversion from mild cognitive impairment to Alzheimer's disease: a 4-year follow-up study. J Neural Transm (Vienna). 2014;121(2):193–200.
Gattaz WF, Forlenza OV, Talib LL, Barbosa NR, Bottino CM. Platelet phospholipase a(2) activity in Alzheimer's disease and mild cognitive impairment. J Neural Transm (Vienna). 2004;111(5):591–601.
Couture R, Blaes N, Girolami JP. Kinin receptors in vascular biology and pathology. Curr Vasc Pharmacol. 2014r;12(2):223–48.
Bicca MA, Costa R, Loch-Neckel G, Figueiredo CP, Medeiros R, Calixto JB. B2 receptor blockage prevents Aβ-induced cognitive impairment by neuroinflammation inhibition. Behav Brain Res. 2015;278:482–91.
Stefanova NA, Fursova A, Kolosova NG. Behavioral effects induced by mitochondria-targeted antioxidant SkQ1 in Wistar and senescence-accelerated OXYS rats. J Alzheimers Dis. 2010;21:479–91.
Genoanalitika Lab, Moscow. [http://www.genoanalytica.ru/]. Accessed 20 Dec 2017.
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMB Net J. 2011;17:10–2. https://doi.org/10.14806/ej.17.1.200
Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14(4):R36.
Stefanova NA, Maksimova KY, Kiseleva E, Rudnitskaya EA, Muraleva NA, Kolosova NG. Melatonin attenuates impairments of structural hippocampal neuroplasticity in OXYS rats during active progression of Alzheimer’s disease-like pathology. J Pineal Res. 2015;59:163–77.
The authors would like to thank Nikita Ershov, who was strongly involved in the analysis of sequencing data.
The RNA-Seq experiments were supported by a Russian Scientific Foundation grant (16–15-10005), whereas the histological examination and electron-microscopic analysis were supported by grants from the Russian Foundation for Basic Research (15–04-01938 and 15–04-06066) and by grant from the government of the Russian Federation, Resolution No. 220 (14.B25.31.0033). The publication costs for this article were funded by a Russian Scientific Foundation grant (16–15-10005).
Availability of data and materials
About this supplement
This article has been published as part of BMC Genomics Volume 19 Supplement 3, 2018: Selected articles from Belyaev Conference 2017: genomics. The full contents of the supplement are available online at https://bmcgenomics.biomedcentral.com/articles/supplements/volume-19-supplement-3.
Ethics approval and consent to participate
The study protocol was approved by the Commission on Bioethics of the Institute of Cytology and Genetics, Novosibirsk, Russia. All the animal experiments were in compliance with the European Communities Council Directive of 24 November 1986 (86/609/EEC). Every effort was made to minimize animal suffering, and rats’ hippocampi were collected after euthanasia by CO2 inhalation.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The lists of significantly enriched GO terms according to DAVID and pathways according to KEGG results, and the list of significantly differentially expressed protein-coding genes associated with cerebrovascular function identified by the RNA-Seq platform. (XLSX 97 kb)
The file contains a Venn diagram for significantly differentially expressed protein-coding genes associated with cerebrovascular function. (XLSX 217 kb)
The file contains significant terms associated with blood vessel processes according to GeneMANIA. (XLSX 47 kb)
About this article
Cite this article
Stefanova, N.A., Maksimova, K.Y., Rudnitskaya, E.A. et al. Association of cerebrovascular dysfunction with the development of Alzheimer’s disease-like pathology in OXYS rats. BMC Genomics 19 (Suppl 3), 75 (2018). https://doi.org/10.1186/s12864-018-4480-9