APOE is featured as a major apolipoprotein that mediates the cholesterol and lipid metabolism upon binding the (very) low-density lipoprotein receptors, including low density lipoprotein receptor (LDLR), very low density lipoprotein receptor (VLDR), and some of the low density lipoprotein receptor-related protein (LRP1-12) genes family . It binds lipids and thus forms intermediate density lipoproteins further cleared from the plasma into the liver by receptor-mediated endocytosis. The liver and brain are major APOE secretion organs in the vertebrates . Liver originated plasma APOE cannot pass the blood brain barrier (BBB), making brain APOE synthesis and metabolism independent .
Notably, ApoE expression in the rodent brain is two orders of magnitude higher than any other apolipoprotein gene as observed from contemporary RNA-Seq data. While basal ApoE gene expression can be observed in virtually all cell types, the cells maintaining the highest ApoE expression, other than hepatocytes in the liver, are brain astrocytes and microglia , , .
The fact that astrocytes synthesize lipids and at the same time highly express ApoE – which mediates cholesterol as well as amyloid and lipid metabolism – underlines the ultimate significance of ApoE functioning specifically in the brain. Notably, astrocytes shed “mini energy stations,” large membrane vesicles comprised of lipids, mitochondria, and ATPs . It was also shown that ApoE eps4 allele carriers exhibit altered energy metabolism pathways in the brain . The role of brain APOE in cholesterol and amyloid beta homeostasis and its implication in Alzheimer’s disease (AD) was presented as a new paradigm in a recent review by . Similarly, the “Triad” of ultimate factors such as ApoE, mitochondrial haplotype and sexual dimorphism as major risk factors of AD etiology was pointed out in . All of these observations underscore vital APOE neuroenergetic components in neural diseases such as AD. Along with ApoE involvement in phagocytosis of cell corpses  and astrocyte mediated synapse elimination in the brain , ApoE stands as a vital gene in brain catabolism.
While cholesterol rich lipoproteins are required in neurons to support synaptogenesis and maintenance of synaptic connections, another function of APOE is delivering Amyloid beta (APP/ABeta)-APOE complexes to cellular catabolism via neuronal APOE receptors . With many APOE/ABeta complexes entering the cytoplasm via neuronal LDLR receptors , another mediator is the pro-neurotrophin receptor, sortilin (SORT1), which also proved to be one of the major endocytic pathways for clearance of APOE/ABeta complexes in neurons , . It was shown that a lack of Sort1 receptor expression in mice results in accumulation of APOE and ABeta in the brain and in aggravated plaque burden .
It is also necessary to mention the signalling role of APOE protein within the brain. It was shown that APOE suppresses food intake via a mechanism enhancing proopiomelanocortin (POMC) signalling in the hypothalamus . ApoE takes part in RhoA  and NMDA mediated  pathways. It also plays a distinct role in microvasculature maintenance . Previous studies revealed that CNS neurons in mice express ApoE in response to excitotoxic injury , implying brain defense functions. ApoE gene is also expressed in primary cultured human CNS astrocytes  and in many human neuronal cell lines, including SY-5Y, Kelly, and NT2 cells .
The human ApoE gene E4 (ApoE4) allele is the major known genetic risk factor for AD. In most clinical studies, ApoE4 carriers account for 65–80 % of all AD cases. ApoE mouse models are widely used in neurobiology and in the pathogenesis of AD . Thus, ApoE locus and closely located neighboring genes are the subject of close attention both in mice and in humans . In addition to ApoE, particular attention has been drawn to its neighbor, Tomm40 . Tomm40 has been reported to be involved in the predisposition to AD and non-pathologic cognitive decline based on “mitochondrial cascade hypothesis,” which proposes a genetic contribution toward mitochondrial durability and function in the brain in particular , , . Due to high linkage disequilibrium between ApoE and Tomm40 in humans, however , it is hard to elucidate an independent impact of Tomm40 on AD development.
The high linkage disequilibrium rate in ApoE extended locus region and its high evolutionary conservation across mammalian species  imply that there might be some impact on gene expression coordination within the locus. In the current study, we analyzed the joint gene expression rate in the vicinity of ApoE in mice for 8 genes.
A specific ApoE implication in stress response metabolism was first reported in 1996 . Further studies indicate a significant interaction between ApoE expression status and stressful life events , . Stress is a factor of enhanced risk for AD . In a series of animal knockout studies ,  it has been postulated that the effects of stress are mediated by glucocorticoids and depend on the presence of APOE. The effects of ApoE expression rate variation may alter susceptibility to environmental factors such as stress. The threshold at which the stress can result in damage to neurons may differ depending on ApoE isoforms proportions . CNS injury is supposed to induce neuronal expression of ApoE to participate in neuronal repair or to protect neurons from injury , .
The present research is designed to study the interference of gene expression variation in the ApoE locus at stressful loads in brain regions of male mice. We used the model of chronic social conflicts in daily agonistic interactions, featuring strong stress induction in the mice. The gene expression in the ApoE locus was primarily studied in the hypothalamus, the main brain region of stress regulation, as well as in other brain regions – the hippocampus, ventral tegmental area, striatum, and midbrain raphe nuclei.
We pursued the issue of the consistency of ApoE expression assessment and its coordination with the genomic context by considering the neighboring genes. Additionally, we targeted the small subset of genes which are the partners of ApoE in metabolic processes and signal transduction to assess their possible interaction in brain regions. Finally, as APOE metabolism involves lipoprotein receptors, we considered a list of lipoprotein receptors to elucidate their performance specifics in various brain regions and their relation to ApoE expression rate.
2 Materials and Methods
Adult male mice of the C57BL/6J strain from Animal Breeding Facility, Branch of Institute of Bioorganic Chemistry of the RAS (ABF BIBCh, RAS) (Pushchino, Moscow region) were used. The animals were housed under standard conditions (12:12 h light/dark regime, switch-on at 8:00 a.m.; food (pellets) and water available (ad libitum). Mice were weaned at one month of age and housed in groups of 8 to 10 in plastic cages (36 × 23 × 15 cm). Experiments were performed on mice 10 to 12 weeks of age. All procedures were in compliance with the European Communities Council Directive 2010/63/EU on September 22, 2010. The study was approved by Scientific Council N 9 of the Institute of Cytology and Genetics SD RAS of March 24, 2010, N 613.
2.2 Generation of Stressful Loads Under Chronic Agonistic Interactions in Male Mice
Stressful loads were induced by daily agonistic interactions in male mice in situations of chronic social conflicts, which were accompanied by positive and negative social experience (wins and defeats) , . Pairs of weight-matched animals were each placed in a steel cage (28 × 14 × 10 cm), bisected by a perforated, transparent partition that allowed the animals to see, hear, and smell each other, but prevented physical contact. The animals were left undisturbed for 2 or 3 days to adapt to new housing conditions and sensory contact before they were exposed to encounters. Every afternoon (14:00–17:00 p.m. local time), the cage lid was replaced with a transparent one. Five minutes later, (the period necessary for individuals’ activation), the partition was removed for 10 min to encourage agonistic interactions. The superiority of one of the mice was firmly established within two or three encounters with the same opponent. The superior mouse behaved by attacking, biting, and chasing another mouse, who in response displayed only defensive behaviors (sideways postures, upright postures, withdrawal, lying on the back, or freezing). If aggressive interaction between males lasted for more than 3 min, or less in some cases, the attack was discontinued by lowering the partition to prevent damage to the loser. Each defeated mouse (loser) was exposed to the same winner for three days, while afterwards, each loser was placed, once a day after the fight, in an unfamiliar cage with an unfamiliar winner behind the partition. Each victorious mouse (winners, aggressive mice) remained in its original cage. This procedure was performed once a day for 20 days and yielded an equal number of winners and losers.
In the behavioral study, the following three groups of animals were used: (a) defeated mice – chronically defeated mice after 20 days of agonistic interactions; (b) aggressive mice – chronically aggressive mice that exhibited daily aggression during twenty days of agonistic interactions; (c) controls – mice without a consecutive experience of agonistic interactions.
Aggressive and defeated mice with the most prominent behavioral phenotypes were selected for the transcriptome analysis. All mice were simultaneously decapitated, 24 h after the last agonistic interaction and the control animals. The brain regions were dissected by one experimenter according to the map presented in the Allen Mouse Brain Atlas (http://mouse.brain-map.org/static/atlas). All biological samples were placed in RNAlater solution (Life Technologies, USA) and were stored at −70 °C until sequencing. In the current study we did not pursue the task of elucidating the differential background between the aggressive and defeated mice, but treated the samples as replicas of affected mice, gaining the insight on their common features.
The brain regions were selected for the analysis based on their functions and localization of neurons of neurotransmitter systems. These areas are as follows: the midbrain raphe nuclei, which comprises a multifunctional brain region and contains the majority of serotonergic neuronal bodies; the ventral tegmental area (VTA), which contains the bodies of dopaminergic neurons, is widely implicated in the natural reward circuitry of the brain, and is important in cognition, motivation, drug addiction, and emotions related to several psychiatric disorders; the striatum, which is responsible for the regulation of motor activity and stereotypical behaviors and is also potentially involved in various processes; the hippocampus, which belongs to the limbic system and is essential for memory consolidation and storage, and plays an important role in neurogenesis, and emotional mechanisms; and the hypothalamus, which regulates the stress reaction and many other physiological functions.
Thus, we maintained 45 samples comprised of 3 groups (control, aggressive mice, defeated mice) and 5 brain regions (hypothalamus, striatum, ventral tegmental area, hippocampus, midbrain raphe nuclei). Each group was represented with 3 replicas.
Cuffdiff software was used  to elucidate differentially expressed genes in the region.
The collected samples were sequenced at JSC Genoanalytica (http://genoanalytica.ru/, Moscow, Russia), and the mRNA was extracted using a Dynabeads mRNA Purification Kit (Ambion, Thermo Fisher Scientific, Waltham, MA, USA). cDNA libraries were constructed using the NEBNext mRNA Library PrepReagent Set for Illumina (New England Biolabs, Ipswich, MA, USA) following the manufacturer’s protocol and were subjected to Illumina sequencing. More than 20 million reads were obtained for each sample. The resulting “fastq” format files were used to align all reads to the GRCm38.p3 reference genome using the TopHat aligner . DAVID Bioinformatics Resources 6.7 (http://david.ncifcrf.gov) was used for the description of differentially expressed gene ontology. The Cufflinks program was used to estimate the gene expression levels in FPKM (fragments per kilobase of transcript per million mapped reads) and subsequently identify the differentially expressed genes in the analyzed and control groups. Each brain area was considered separately in nine animals total. Only annotated gene sequences were used in the following analysis. Genes were considered differentially expressed at p ≤ 0.05. q < 0.05 was also taken into consideration.
We have previously conducted studies of gene expression in males in similar experiments using the RT-PCR method with larger samples for each compared experimental group, i.e. aggressive and defeated mice (>10 animals). The direction and extent of changes in expressions of the Tph2, Slc6a4, Bdnf, Creb1, and Gapdh genes in the midbrain raphe nuclei of defeated or aggressive males compared with the control produced by two methods, including RT-PCR ,  and RNA-Seq, are generally consistent. This finding suggests that the transcriptome analyses of the data provided by the Company Genoanalitika have been verified, and this method reflects the actual processes that occur in the brain under our experimental paradigm.
2.4 Previously Published Expression Data in Mouse Brain
In order to cross-validate the results obtained, we employed the unique resource from Barres Lab, Stanford University, USA , . The brain cell type specific RNA-Seq data proved to be highly concordant with our RNA-Seq data pool when we averaged them to use as a “cell mix,” represented in our brain region samples. In particular, when we compared average FPKM values of 32 apolipoprotein (Apo*) genes expression, we achieved Pearson correlation coefficient to be r = 0.9943; df = 31; p < 10−5 , which underlines significant consistency in expression variation. It was also consistent for other genes as will be shown in the results section.
2.4.1 Mouse RNA-Seq Data of Barres Lab in Seven Cortex Cell Types
For mouse cortex, 7 cell types were retrieved with original methods , : astrocytes, neurons, oligodendrocyte precursor cells, newly formed oligodendrocytes, myelinating oligodendrocytes, microglia, and endothelial cells.
2.5 Genes of Interest in the ApoE Locus
Based on the Hi-C data on chromatin domains  and linkage disequilibrium in humans, we selected 7 genes in the vicinity of ApoE, that are Bcam, Pvrl2, Tomm40, ApoE, Apoc1, Apoc2, Apoc4, Clptm1 ordered by DNA strand. The short description of the genes is presented in Table 1. All are located on the same DNA strand. As was reported previously , the gene order within the 1mb of the locus is highly conservative throughout mammals, which implies evolutionary constrains on the locus .
We also targeted low-density lipoprotein receptors expression rates to assess the correlation with ApoE (Table 2). Additionally, we included Pomc, which is a hypothalamus-specific gene involved in neuroendocrine stress response and ApoE-mediated signal transduction, as reported previously . The large molecule of POMC is the source of several important, biologically active substances which may be connected with regulation of lipid metabolism, stress, and psychoemotional state. All of them changed expression rate in aggressive and defeated mice relative to the controls.
2.6 Statistical Analysis
Pearson correlation coefficient was used in all instances of correlation analysis. The significance rates for Pearson correlation values were derived based on Student t-test with corresponding degrees of freedom.
Principal component (PC) analysis has been carried out using XLStat software (www.xlstat.com). Pearson correlation coefficient has been used as a similarity metric during pairwise similarity matrix construction in PC analysis.
The multi-dimensional scaling (MDS) plots were built with XLStat software as well on the basis of Pearson pairwise similarity matrix. The target functional for minimization was Kruskal’s stress statistic. Initial configuration was set to random. Number of repetitions was set to 5. Stop conditions were: convergence step less that 1E−5; maximal number of iterations: 500.
ANOVA analysis has been carried out by XLStat software (www.xlstat.com). Fisher test was used for elucidation of variation significance in ANOVA.
3 Results and Discussion
3.1 ApoE Expression Locus Rate Evaluation and Data Assessment
The 7 adjacent genes at ApoE region may be grouped according to their location relative to cell compartments (Table 1). Bcam and Clptm1 genes are cell membrane receptors, Pvrl2 is a cell membrane protein, Tomm40 is a component of mitochondrial import translocase complex providing import of proteins into mitochondria, and ApoE, Apoc1, Apoc2, and Apoc4 are secreted proteins.
To evaluate the consistency of ApoE gene expression rate in mouse brain regions, we selected 7 consequent genes within ApoE location vicinity and assessed their joint average expression rate from our data (Supplementary Table 1), as well as from Barres Lab (, ; Supplementary Table 2). The compilation of averaged FPKM values across correspondent samples, as well as their standard deviation in mouse brain, for 2 projects are presented in Table 3/Figure 1. As it can be seen (Table 3; Figure 1), while genes from two projects maintain the same order of FPKM values, the ApoE expression rate is more than an order of magnitude higher than that of the neighboring genes. As the data was drawn from various brain regions/cell types, we present the overall compiled data set in Supplementary Table 1.
As we witnessed tight coordination within apolipoprotein cluster (Table 4), and at the same time negligible expression values of Apoc2 and Apoc4 in data considered (Table 3), we omitted from further consideration Apoc2, Apoc4 genes. Thus, we further considered Bcam, Pvrl2, Tomm40, ApoE, Apoc1 and Clptm1 genes.
3.2 ApoE Expression Rate in the Control, Aggressive and Defeated Mice Based on RNA-Seq Sequencing
The differential gene expression list obtained with CuffDiff analysis and corresponding FPKM values are presented in Supplementary Table 3. We assessed ApoE expression rate in the control mice, and the mice with the chronic agonistic interactions experience which was accompanied by stress, pooled into 3 different groups (see Materials and Methods; Supplementary Tables 1 and 3). We observed the consistent significant elevation of ApoE expression rate between the control and affected groups (Figure 2; Supplementary Table 3). The pairwise expression correlation values are presented in Table 4.
3.3 Proopiomelanocortin (Pomc) and App Genes are Positively Correlated with ApoE Expression in Mouse Brain, While Sort1 is Antagonistic to ApoE
As mentioned above, ApoE is grossly associated with AD etiology. Based on literature mining of other AD-associated genes, we selected 3 genes that were reported to associate with AD, and also maintain highly non-random expression covariation with ApoE in our data (Table 5). Their descriptions are presented below.
Pomc gene is the source of up to 10 enzymatically derived hormone related peptides involved in pain and energy homeostasis, melanocyte stimulation, and immune modulation and is involved in neuroendocrine stress response . Pomc function in the hypothalamus neurons differs from that of the pituitary gland, another major Pomc source . In particular, Pomc gene is known for the regulation of food intake in rats in the presence of ApoE . These facts underline the involvement of the ApoE-Pomc tandem in specific energy metabolism signalling cascades in the brain.
APP, the main constituent of Amyloid plagues, is targeted by APOE . It was elucidated recently that APOE epsilon 2, 3, 4 isoforms differentially stimulate APP expression in neurons ; Supplementary Figure 1), APOE eps4 isoform being the highest stimulator.
Sortilin 1 (SORT1) is also considered the APOE- related receptor responsible to major uptake of APOE-APP complexes  and is specifically expressed in astrocytes ,  (Supplementary Figure 1). One of the closest homologues of Sort1 involved in the process is sorting protein-related receptor (SORL1) which is associated with AD pathology in human  and is preferentially expressed in microglia (compared to Sort1) ,  (Supplementary Figure 1).
We found that both Pomc and App expression rates are significantly correlated with ApoE expression (r = 0.618; p < 3E−6 and r = 0.659; p < 4E−7, correspondingly; df = 44; Table 5). Notably, though the correlation values of Pomc and App with ApoE are similar, they constitute significant differences in gene pathways. In particular, we observed Pomc expression solely in the hypothalamus of the mice studied, consistent with previous research studies , . In contrast, non-basal App, as well as ApoE expression, was observed throughout all 45 samples considered. On the neural cells level, ApoE and App maintain antagonistic expression throughout 7 cell types  (Supplementary Figure 1). Pearson correlation rate across the cells is r = −0.91; df = 6; p < 6E−6). The highest ApoE expression is observed in microglia and astrocytes with App maintaining minimal expression within them. So the ApoE-App expression correlation observed across brain regions reflects coordinated inter-cellular interaction between asrocytes and neurons reported in . Thus, it could be interpreted that ApoE-App covariation represents the metabolic interaction background observed in all brain regions, while Pomc pairing with Apoe apparently implies a specific (signalling) pathway, particularly in the hypothalamus.
Consistent with , Sort1 is negatively correlated with all other genes considered (Table 5): ApoE (r = −0.45; p < 1E−3), App (r = −0.42; p < 2E−3) and Pomc (r = −0.578; p < 1.5E−5). Considering Sort1/Sorl1 genes interaction, their expression was highly correlated. It was observed in all samples, and both genes’ expression was significantly down in the hypothalamus featuring the highest ApoE/App expression. The significant negative association of Sorl1 with ApoE (r = −0.49; df = 44; p < 3E−4; not shown) in the brain regions may mediate Sorl1 highly confident AD association reported previously ,  (Table 5). We did not manage to provide any additional interpretation on the Sort1/Sorl1 negative correlation with Pomc, ApoE, App observed in this instance besides partly presented in .
3.4 Lipoprotein Receptors and ApoE Expression Covariation in Brain Regions
We assessed expression rate of available (with average FPKM > 0.2) 10 Lrp* receptors in ours and 10 Lrp* receptors in  projects, and also added to the list LDLR, VLDLR, and Sort1 receptors (Figure 3). As can be seen from Figure 3 based on  data, 4 receptors (Lrp1, Lrp10, Ldlr and Sort1) manifest major receptor expression profile in neural cell types by expression rate.
In our data we observed similar major receptors’ expression pattern (Supplementary Table 1). Also, we found that Sort1 expression rate is significantly elevated in the hippocampus and downregulated in the hypothalamus, being contrary to ApoE. Lrp11 also maintains the hippocampus-specific expression, while Lrp1 can be considered a hypothalamus-specific receptor (Supplementary Table 1).
Overall expression rate of lipoprotein receptors considered was 170 ± 7 FPKM units in four brain regions, and in the hippocampus it is 200.4. Thus, overall lipoprotein receptors’ expression homeostasis proved to be rather conserved across brain regions, except for in the hippocampus due to the Sort1 expression outburst (Supplementary Table 1).
To elaborate on the relations between ApoE and lipoprotein receptors expressed in mouse brain, we built up multi-dimensional scaling (MDS) plots based on pairwise Pearson correlations of 17 Ldr* receptors and the ApoE gene (Figure 4A, B). We observed 3 statistically significant clusters, which consisted of one ApoE-synergistic and 2 ApoE-antagonistic ones.
Lrp1 expression rate is highly correlated with ApoE expression profile across 45 samples (r = 0.858; df = 44; p < 7.6E−7; Figure 4A), next is Lrp10 (r = 0.767; df = 44; p < 1.2E−5; Figure 4A). It was also proved that Lrp1 is differentially expressed in aggressive and defeated mice, compared to the control in the hypothalamus manifesting the same regression pattern (Table 3). Thus, Lrp1 is confirmed as the major partner of ApoE in the brain, as reported previously .
Data from  confirms the high correlation rate elucidated through our data (Figure 4A) between ApoE and Lrp1 (r = 0.925; df = 11; p < 1.2E−5), and Lrp5 (r = 0.65; df = 11; p < 0.011) across various brain cell types (Figure 4B), though other clusters are distinct from brain regions’ ones.
We also performed an ANOVA analysis of brain region-specific variation of ApoE and lipoprotein receptors to assess non-randomness of lipoprotein expression rate. R2 residuals’ values correspond to non-randomness of the brain-specific expression. As can be seen from the resulting Table 6, most lipoprotein receptors maintain highly non-random brain region-specific expression. Its specificity ultimately bifurcates between the hypothalamus and hippocampus/striatum clusters (Table 6). Based on the results, we can assume that each brain region is characterized with specific lipoprotein receptors’ expression profile, some of which (Lrp1, Lrp10, Lrp5) respond to ApoE expression rate alteration.
4 General Discussion
Bearing in mind the distinct, confirmed involvement of the ApoE locus in neuropathologic etiology, including AD, we aimed at an analysis of the ApoE gene and its genomic and pathway contexts using the expression rate co-variation assessment across the sample of 45 brain region RNA-Seq sets. The contemporary RNA-Seq data ,  have unambiguously confirmed that the astrocytes and microglia are the major sources of ApoE (Supplementary Table 2), achieving as high as 32540.3 FPKM units for certain exons’ peak instances in microglia ,  (Supplementary Table 4). Thus, one overlooked issue concerning the ApoE phenomenon is its high expression rate in the mouse brain. This, nevertheless, obviously corroborates the fact that the brain is the “fattiest” organ in the human body and is composed of roughly 60 % fat, while APOE is needed for lipid metabolism in the cells. According to the human protein atlas resource, APOE protein expression maintains the highest score in brain tissues (http://www.proteinatlas.org/ENSG00000130203-APOE/tissue), and it is even higher than in the liver for protein expression. This fact underscores the intense transcription of ApoE which belongs to the top 20 highest expressed genes in mouse brain regions.
The overall expression rate of ApoE gene in all mouse brain regions considered proved to be extraordinarily high, and, in our case, it averaged 1100 FPKM units (Table 7) across brain regions. To get a more elaborated view and underline ApoE high expression rate, we compiled the list of genes which maintain average expression rate more than 700 in FPKM values in our data. The results are presented in Table 7.
The expression rate of the ApoE gene enters the list of the top 21 highest, expressed both in our sample (Table 7) and one from top 17 genes from mouse brain cell types data  (Table 8). Discarding ribosomal and other non-specific genes, we may state that ApoE expression rate in mouse is uniquely high.
The preliminary analysis of highly expressed gene ratio revealed that some of the genes represent non-specific, ubiquitous proteins such as ribosomal proteins and others. We replicated the results for previously published data ,  (Table 8), and found that 5 genes are replicated in both lists (Table 7 and Table 8; marked in bold) underscoring highly non-random replication of highly expressed genes in the brain for distinct transcriptome projects (p < 1.5E−11 for random null hypothesis of 5 genes coincidence based on binomial test for 20 top genes pool with p = 0.001).
A class of highly expressed genes is reported to maintain specific means for their regulation and features . In particular, ribosomal proteins encoding genes maintain specific transcription subunit TATA-box-binding protein-related factor 2 (TRF2) as a core subunit for their transcription , . A significant part of highly expressed genes maintain miRNA and other noncoding RNA within their vicinity . Also, CTCF binding sites and remote enhancers are common for their expression regulation . We also observed a range of CTCF binding sites along the ApoE region.
Interestingly, while ApoE expression rate differs significantly between astrocytes and cortical neurons (FPKM values of 3006.2 vs. 200.4, correspondingly; Supplementary Table 2; Figure 1) , the amount of highly expressed genes proves to be the closest in these two cell types compared to other 5 cell types , revealing that they are structurally similar by that metric (highly expressed genes number; HEGN) . While close co-variation between ApoE and its neighbors has been observed throughout all brain regions considered (Table 4), the highest ApoE elevation rate of affected groups related to the controls was observed in the hypothalamus (Figure 5).
The genes coordination rate has been the highest therein, as well as across brain regions (Supplementary Table 3). Notably, one of the highest expression covariation rates with ApoE manifests Tomm40 (r = 0.964 df = 8; p < 1.3E−13; Figure 5; Supplementary Table 3). Tomm40 expression modulation might affect the mitochondrial metabolism homeostasis reported previously .
The issue of gene co-expression and the tendency for similar expressed genes to cluster in eukaryotes have been studied throughout the genome research era . Still, no implicit evidence for cystronic-like organization for the majority of eukaryotic genes exists and such cases are majorly confined to closely located paralogous gene clusters .
We should herein distinguish the functional and non-specific co-expression coordination. While the co-regulation issue discussed in the previous paragraph refers to the former instance, in the case of ApoE, we observed non-specific co-regulation based on disproportional expression rates of genes located densely at the chromosome  while maintaining focal highly expressed ApoE gene. Non-specific co-expression may thus arise from an abundance of transcription factors and chromatin accessibility within the region due to high transcription intensity at ApoE locus.
Based on the covariation analysis, we also found that ApoE turnover by lipoprotein receptors manifests highly complicated, brain region-specific homeostasis which may affect its proper function in a case of misbalance. While stated earlier , this is the first report of highly significant Apoe-LRP1 covariation based on RNA-Seq data, elucidating the major receptor counterpart for ApoE across brain regions.
Considering the mouse model of chronic agonistic interactions, we observed the highly statistically significant elevation (FDR < 0.007 in aggressive mice and FDR < 0.001 in defeated mice) of ApoE gene expression in the hypothalamus of affected mice, but not in other brain regions. The observation was additionally supported by a concordant expression shift of surrounding genes. The high correlation between neighboring genes expression in the vicinity of the ApoE gene was observed in all brain regions (Table 4). The increased expression of the ApoE gene might represent a metabolism rate increase featuring the behaviour affection. In particular, the stabilization of ApoE expression in the course of aging by means of aerobic exercise was proved to prevent astrocyte dysfunction in mice . Notably, the absence of ApoE (ApoE deficient mice) did not reveal that relation . High elevation of expression of the ApoE gene in the hypothalamus, both in the defeated and the aggressive mice, directly indicates the effects of social stress and anxiety on both participants of agonistic interactions, as has been shown previously , .
Overall, the causes and consequences of ApoE expression elevation we observed in hypothalamus RNA-Seq data can be interpreted as it is depicted in Figure 6.
The major cause of ApoE elevation in the hypothalamus can be the neuroendocrine stress response, which induces Pomc expression and adrenocorticotropic hormone in particular . We observed high Pomc elevation in stressed mice compared to mice in the control group (Supplementary Table 1; p < 2E−5 for stress groups merged).
The interaction of Pomc and ApoE has been reported previously , , and we observe their covariation with high confidence (r = 0.618; p < 3E−6). It may be considered both ways, since intracerebroventricular administration of APOE was reported to increase Pomc expression .
To summarize the causes and consequences of ApoE expression elevation (Figure 6), we would like to stress the issues we observed and reported in the results:
We report significant covariation of ApoE and App expression (r = 0.659; p < 2E−7). It was elucidated in a recent paper  that APOE induction rate of APP is isoform dependent, with APOE eps4 inducing App expression the most, which might explain the APOE allelic effect on AD etiology.
Previously we observed a metabolic increase in the hypothalamus of stress-induced mice based on Rpl/Rps cytoplasmic ribosomal genes expression elevation . The overall index of covariation of ApoE with Rps, Rpl genes was equal to 0.65 (p < 1E−6).
The Apoe expression elevation in mice under social stress (Supplementary Table 1) impacts the neighboring genes apparently in a non-specific manner. The covariation of Bcam, Pvrl2, Tomm40, Apoc1, and Clptm1, listed in Table 4, results in the following corresponding p-values: 1.37E−08, 7.88E−09, 2.96E−07, 4.85E−06, 8.63E−06. In the hypothalamus only, the corresponding values are: 3.14E−08, 4.29E−14, 1.35E−13, 0.226, 0.003.
Due to the sufficient independent replications number in RNA-Seq data, we were able to derive highly confident covariation assessments. Coordination with brain specific cell type data ,  accommodated cross-validation of our data.
This study was supported by the Russian Science Foundation (Nо. 14-15-00063).
Liu M, Kuhel DG, Shen L, Hui DY, Woods SC. Apolipoprotein E does not cross the blood-cerebrospinal fluid barrier, as revealed by an improved technique for sampling CSF from mice. Am J Physiol Regul Integr Comp Physiol. 2012;303:R903–8. CrossrefPubMedGoogle Scholar
Zhang Y, Chen K, Sloan SA, Bennett ML, Scholze AR, O’Keeffe S, et al. An RNA-sequencing transcriptome and splicing database of glia, neurons, and vascular cells of the cerebral cortex. J Neurosci. 2014;34:11929–47. CrossrefPubMedGoogle Scholar
Zhang Y, Sloan SA, Clarke LE, Caneda C, Plaza CA, Blumenthal PD, et al. Purification and characterization of progenitor and mature human astrocytes reveals transcriptional and functional differences with mouse. Neuron. 2016;89:37–53. CrossrefPubMedGoogle Scholar
Falchi AM, Sogos V, Saba F, Piras M, Congiu T, Piludu M. Astrocytes shed large membrane vesicles that contain mitochondria, lipid droplets and ATP. Histochem Cell Biol. 2013;139:221–31. CrossrefPubMedGoogle Scholar
Perkins M, Wolf AB, Chavira B, Shonebarger D, Meckel JP, Leung L, et al. Altered energy metabolism pathways in the posterior cingulate in young adult apolipoprotein E epsilon 4 carriers. J Alzheimers Dis. 2016;53:95–106. CrossrefGoogle Scholar
Grainger DJ, Reckless J, McKilligin E. Apolipoprotein E modulates clearance of apoptotic bodies in vitro and in vivo, resulting in a systemic proinflammatory state in apolipoprotein E-deficient mice. J Immunol. 2004;173:6366–75. CrossrefGoogle Scholar
Chung WS, Verghese PB, Chakraborty C, Joung J, Hyman BT, Ulrich JD, et al. Novel allele-dependent role for APOE in controlling the rate of synapse pruning by astrocytes. Proc Natl Acad Sci USA. 2016;113:10186–91. CrossrefGoogle Scholar
Castellano JM, Deane R, Gottesdiener AJ, Verghese PB, Stewart FR, West T, et al. Low-density lipoprotein receptor overexpression enhances the rate of brain-to-blood Abeta clearance in a mouse model of beta-amyloidosis. Proc Natl Acad Sci USA. 2012;109:15502–7. CrossrefGoogle Scholar
Carlo AS, Gustafsen C, Mastrobuoni G, Nielsen MS, Burgert T, Hartl D, et al. The pro-neurotrophin receptor sortilin is a major neuronal apolipoprotein E receptor for catabolism of amyloid-beta peptide in the brain. J Neurosci. 2013;33:358–70. CrossrefPubMedGoogle Scholar
Casey CS, Atagi Y, Yamazaki Y, Shinohara M, Tachibana M, Fu Y, et al. Apolipoprotein E inhibits cerebrovascular pericyte mobility through a RhoA protein-mediated pathway. J Biol Chem. 2015;290:14208–17. CrossrefPubMedGoogle Scholar
May P, Rohlmann A, Bock HH, Zurhove K, Marth JD, Schomburg ED, et al. Neuronal LRP1 functionally associates with postsynaptic proteins and is required for normal motor function in mice. Mol Cell Biol. 2004;24:8872–83. CrossrefPubMedGoogle Scholar
Xu Q, Bernardo A, Walker D, Kanegawa T, Mahley RW, Huang Y. Profile and regulation of apolipoprotein E (ApoE) expression in the CNS in mice with targeting of green fluorescent protein gene to the ApoE locus. J Neurosci. 2006;26:4985–94. PubMedCrossrefGoogle Scholar
Dupont-Wallois L, Soulie C, Sergeant N, Wavrant-de Wrieze N, Chartier-Harlin MC, Delacourte A, et al. ApoE synthesis in human neuroblastoma cells. Neurobiol Dis. 1997;4:356–64. CrossrefPubMedGoogle Scholar
Payton A, Sindrewicz P, Pessoa V, Platt H, Horan M, Ollier W, et al. A TOMM40 poly-T variant modulates gene expression and is associated with vocabulary ability and decline in nonpathologic aging. Neurobiol Aging. 2016;39:217e1–7. CrossrefGoogle Scholar
Roses A, Sundseth S, Saunders A, Gottschalk W, Burns D, Lutz M. Understanding the genetics of APOE and TOMM40 and role of mitochondrial structure and function in clinical pharmacology of Alzheimer’s disease. Alzheimers Dement. 2016;12:687–94. PubMedCrossrefGoogle Scholar
Sheffler J, Moxley J, Sachs-Ericsson N. Stress, race, and APOE: understanding the interplay of risk factors for changes in cognitive functioning. Aging Ment Health. 2014;18:784–91. PubMedCrossrefGoogle Scholar
de Kloet ER, Grootendorst J, Karssen AM, Oitzl MS. Gene x environment interaction and cognitive performance: animal studies on the role of corticosterone. Neurobiol Learn Mem. 2002;78:570–7. PubMedCrossrefGoogle Scholar
Grootendorst J, Kempes MM, Lucassen PJ, Dalm S, de Kloet ER, Oitzl MS. Differential effect of corticosterone on spatial learning abilities in apolipoprotein E knockout and C57BL/6J mice. Brain Res. 2002;953:281–5. CrossrefPubMedGoogle Scholar
Peavy GM, Lange KL, Salmon DP, Patterson TL, Goldman S, Gamst AC, et al. The effects of prolonged stress and APOE genotype on memory and cortisol in older adults. Biol Psychiatry. 2007;62:472–8. CrossrefPubMedGoogle Scholar
Boyarskikh UA, Bondar NP, Filipenko ML, Kudryavtseva NN. Downregulation of serotonergic gene expression in the Raphe nuclei of the midbrain under chronic social defeat stress in male mice. Mol Neurobiol. 2013;48:13–21. CrossrefPubMedGoogle Scholar
Smagin DA, Boyarskikh UA, Bondar NP, Filipenko ML, Kudriavtseva NN. Reduction of serotonergic gene expression in the raphe nuclei of the midbrain under positive fighting experience in male mice. Adv Biosci Biotechnol. 2013;4:36–44. CrossrefGoogle Scholar
Rao SS, Huntley MH, Durand NC, Stamenova EK, Bochkov ID, Robinson JT, et al. A 3D map of the human genome at kilobase resolution reveals principles of chromatin looping. Cell. 2014;159:1665–80. CrossrefPubMedGoogle Scholar
Jenks BG. Regulation of proopiomelanocortin gene expression an overview of the signaling cascades, transcription factors, and responsive elements involved. Ann NY Acad Sci. 2009;1163:17–30. CrossrefGoogle Scholar
Huang YW, Zhou B, Wernig M, Sudhof TC. ApoE2, ApoE3, and ApoE4 differentially stimulate APP transcription and A beta secretion. Cell. 2017;168:427.e21–441. Google Scholar
Trakhtenberg EF, Pho N, Holton KM, Chittenden TW, Goldberg JL, Dong L. Cell types differ in global coordination of splicing and proportion of highly expressed genes. Sci Rep. 2016;6:32249. PubMedCrossrefGoogle Scholar
Hurst LD, Ghanbarian AT, Forrest AR, Huminiecki L, Consortium F. The constrained maximal expression level owing to haploidy shapes gene content on the mammalian X chromosome. PLoS Biol. 2015;e100231513. PubMedGoogle Scholar
Engreitz JM, Haines JE, Perez EM, Munson G, Chen J, Kane M, McDonel PE, Guttman M, Lander ES. Local regulation of gene expression by lncRNA promoters, transcription and splicing. Nature. 2016;539:452–5. CrossrefPubMedGoogle Scholar
Eldholm V, Haugen A, Zienolddiny S. CTCF mediates the TERT enhancer- promoter interactions in lung cancer cells: identification of a novel enhancer region involved in the regulation of TERT gene. Int J Cancer. 2014;134:2305–13. CrossrefPubMedGoogle Scholar
Soto I, Graham LC, Richter HJ, Simeone SN, Radell JE, Grabowska W, et al. APOE stabilization by exercise prevents aging neurovascular dysfunction and complement induction. PLoS Biol. 2015;e100227913. PubMedGoogle Scholar
Hoeve ES, Kelly G, Luz S, Ghanshani S, Bhatnagar S. Short-term and long-term effects of repeated social defeat during adolescence or adulthood in female rats. Neuroscience. 2013;249:63–73. CrossrefPubMedGoogle Scholar
Raber J, Akana SF, Bhatnagar S, Dallman MF, Wong D, Mucke L. Hypothalamic-pituitary-adrenal dysfunction in Apoe(-/-) mice: possible role in behavioral and metabolic alterations. J Neurosci. 2000;20:2064–71. Google Scholar
Smagin DA, Kovalenko IL, Galyamina AG, Orlov YL, Babenko VN, Kudryavtseva NN. Heterogeneity of brain ribosomal genes expression following repeated experience of aggression in male mice as revealed by RNA-Seq. Mol Neurobiol 2017. (Epub ahead of print). Google Scholar
The online version of this article (https://doi.org/10.1515/jib-2017-0024) offers supplementary material, available to authorized users.
About the article
Published Online: 2017-09-13
Conflict of interest statement: Authors state no conflict of interest. All authors have read the journal’s Publication ethics and publication malpractice statement available at the journal’s website and hereby confirm that they comply with all its parts applicable to the present scientific work
Citation Information: Journal of Integrative Bioinformatics, Volume 14, Issue 3, 20170024, ISSN (Online) 1613-4516, DOI: https://doi.org/10.1515/jib-2017-0024.
©2017, Vladimir N. Babenko et al., published by De Gruyter, Berlin/Boston. This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 3.0 License. BY-NC-ND 3.0