Next Article in Journal
A Molecular Networking Strategy: High-Throughput Screening and Chemical Analysis of Brazilian Cerrado Plant Extracts against Cancer Cells
Next Article in Special Issue
Fenamates as Potential Therapeutics for Neurodegenerative Disorders
Previous Article in Journal
CALINCA—A Novel Pipeline for the Identification of lncRNAs in Podocyte Disease
Previous Article in Special Issue
Invertebrate Models Untangle the Mechanism of Neurodegeneration in Parkinson’s Disease
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cells
Volume 10
Issue 3
10.3390/cells10030693
cells-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Voluntary Wheel Running Did Not Alter Gene Expression in 5xfad Mice, but in Wild-Type Animals Exclusively after One-Day of Physical Activity

by
Anna Wierczeiko
1,2,†,
Lena Gammel
3,†,
Konstantin Radyushkin
2,
Vu Thu Thuy Nguyen
3,
Hristo Todorov
1,
Susanne Gerber
1,* and
Kristina Endres
3,*
1
Working Group Computational Systems Genetics (CSG), Institute of Human Genetics, University Medical Center, Johannes Gutenberg University, 55131 Mainz, Germany
2
Working Group Mouse Behavioral Unit (MBU), Leibniz Institute for Resilience Research (LIR), 55122 Mainz, Germany
3
Working Group Healthy Aging and Neurodegeneration, Department of Psychiatry and Psychotherapy, University Medical Center, Johannes Gutenberg University, 55131 Mainz, Germany
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Cells 2021, 10(3), 693; https://doi.org/10.3390/cells10030693
Submission received: 10 February 2021 / Revised: 1 March 2021 / Accepted: 17 March 2021 / Published: 20 March 2021
(This article belongs to the Collection Advances in Neurodegenerative Disease)
Browse Figures
Versions Notes

Abstract

:
Physical activity is considered a promising preventive intervention to reduce the risk of developing Alzheimer’s disease (AD). However, the positive effect of therapeutic administration of physical activity has not been proven conclusively yet, likely due to confounding factors such as varying activity regimens and life or disease stages. To examine the impact of different routines of physical activity in the early disease stages, we subjected young 5xFAD and wild-type mice to 1-day (acute) and 30-day (chronic) voluntary wheel running and compared them with age-matched sedentary controls. We observed a significant increase in brain lactate levels in acutely trained 5xFAD mice relative to all other experimental groups. Subsequent brain RNA-seq analysis did not reveal major differences in transcriptomic regulation between training durations in 5xFAD mice. In contrast, acute training yielded substantial gene expression changes in wild-type animals relative to their chronically trained and sedentary counterparts. The comparison of 5xFAD and wild-type mice showed the highest transcriptional differences in the chronic and sedentary groups, whereas acute training was associated with much fewer differentially expressed genes. In conclusion, our results suggest that different training durations did not affect the global transcriptome of 3-month-old 5xFAD mice, whereas acute running seemed to induce a similar transcriptional stress state in wild-type animals as already known for 5xFAD mice.

1. Introduction

As the population ages, the number of people with degenerative diseases is simultaneously increasing. Dementia is the fifth most frequent cause of death worldwide, with Alzheimer’s disease (AD) as the most prevalent type of dementia in the elderly [1]. AD is characterized by a progressive decline in memory functions, behavioral impairments, and the loss of social skills, which means that quality of life begins to deteriorate well before death [2,3]. Key neuropathological features of this devastating disorder include the accumulation of extracellular β-amyloid (Aβ) plaques, intracellular neurofibrillary tangles (NFT) composed of hyperphosphorylated tau proteins, and region-specific synaptic as well as cellular degeneration, together with an increase in inflammation and oxidative stress [4,5,6,7]. In addition, metabolic dysfunctions such as abnormal glucose uptake and brain insulin resistance, likely increasing degeneration and cognitive impairment, have also been described in AD patients [8]. Besides the known pathological hallmarks, the leading cause of progressive brain atrophy is still unknown, explaining the lack of successful AD treatments. It became, however, generally accepted that the underlying mechanisms are polyfactorial and depend on the complex interplay of multiple (partly unknown) genetic and non-genetic variables [9,10,11,12]. Furthermore, it is widely recognized that AD already seems to develop decades before clinical symptoms occur [13,14,15,16]. Since prophylactic pharmacological therapies in advance of a possible disease outbreak might be ethically not reasonable or hampered by lack of a long-term compliance, lifestyle-modifying interventions, and preventions such as diet and physical activity have become increasingly attractive in the neurodegenerative research field.
Physical activity (e.g., aerobic exercise or resistance training) has already been associated with several beneficial effects that reduce the risk of developing AD in both, humans and animal models [17,18,19,20]. The duration of physical activity being studied can be divided into acute, one single bout of physical activity, and chronic, consistent activity over an extended time period [21]. Several studies suggested that regular training can improve reaction time, increase the size of the hippocampus (the brain region mainly responsible for learning and memory), and delay age-related memory impairments as well as the overall cognitive decline in older individuals with and without mild-cognitive impairment (MCI) or early AD [18,22,23,24]. Even after just one exercise session, short- and long-term memory and immune cell function were presumably stimulated in normal aging and MCI cohorts [25,26,27]. In transgenic rodent AD models, chronically trained animals exhibited a reduction in Aβ deposition and tau phosphorylation, anti-inflammatory modifications, improved brain oxygenation, and general cognitive improvements as compared to sedentary littermates [13,28,29,30,31,32,33]. Furthermore, another study showed that only three days of voluntary wheel running was associated with enhanced neurogenesis in rats [34].
However, there is still insufficient evidence that physical activity provides an effective method of fighting AD, as there are also several studies that have not observed a positive link between physical activity and brain function improvement [35,36,37,38,39]. Recently, Hansson et al. conducted a large longitudinal study over 20 years that indicated no reduction in the frequency of developing AD by comparing exercising participants (n = 197,685) with sedentary individuals (n = 197,684) [37]. The same results were obtained in a subsequent study on chronically trained AD model mice (5xFAD) [37]. The evaluation of neuroinflammation and non-cognitive parameters in voluntarily running young 5xFAD mice, in which no pathological features had yet developed, also showed no changes compared to sedentary mice. Moreover, signs of disease acceleration could even be identified within the trained transgenic group [38]. The inconsistency of physical activity effects on cognitive function and AD could likely be explained by different starting points such as age or disease state, duration (e.g., acute or chronic), and activity types [18].
On the neurochemical level, specific neurotransmitters, neuromodulators, and lactate concentrations were shown to increase after short and long-term exercise [40]. Furthermore, some genes and proteins that are predominantly involved in neurogenesis were suggested to be differentially expressed after both, acute and regular physical activity, in humans and in animal models, with brain-derived neurotrophic factor (BDNF) being the most prominent [34,41,42]. The global transcriptional response to different types of physical activity in AD has rarely been investigated. To better understand the molecular biology behind physical activity in AD, we examined in this study the blood and brain lactate levels as well as the transcriptomes of 3-month-old transgenic 5xFAD and wild-type mice after one day (acute) or four weeks (chronic) voluntary wheel running. Comparing the gene expression of 5xFAD mice at an age at which the pathology is not severely manifested with age-matched controls could provide further insights into fundamental molecular mechanisms of short- and long-term exercise in both AD and healthy individuals.

2. Materials and Methods

2.1. Animals

Male B6SJL-Tg(APPSwFlLon, PSEN1*M146L*L286V)6799Vas/Mmjax (5xFAD) mice (Jackson Lab, Bar Harbor, Maine, USA) were crossbred with female C57BL/6J mice and maintained as heterozygous transgenics on the C57BL/6J background. The animals were single-caged in type II cages with a 12-h day/night cycle. Food and water were available ad libitum. The non-transgenic littermates were used as the control. All procedures were performed in accordance with the European Communities Council Directive regarding care and use of animals for experimental procedures and were approved by local authorities (Landesuntersuchungsamt Rheinland-Pfalz; approval number G 14-1-087). Mice were randomly assigned to three groups: untrained (sedentary), access to a saucer wheel (15-cm diameter; med associates Inc.) for 30 days, or access to a saucer wheel for 18 h before sacrifice. The saucer wheel contained a WLAN connection to a computer for measuring running counts. Food consumption and body weight were assessed regularly.

2.2. Tissue Preparation

Mice were sacrificed after isoflurane anesthesia, and truncal blood was collected. Blood glucose was measured immediately from two independent droplets with an AccuCheck mobile device (Roche). Heart, muscle, abdominal fat, and the left brain-hemisphere were dissected and shock-frozen in liquid nitrogen with subsequent storage at −80 °C. The right hemisphere was cut into small cubicles and incubated at Room Temperature for 20 min in RNAlater (Qiagen, Hilden, Germany) before being stored at −80 °C. Serum was obtained from truncal blood by two-times centrifugation after clotting and stored at −80 °C.

2.3. Lactate Measurement

Brain samples were homogenized with a Tissue Lyzer, and stainless-steel beads (Qiagen, Hilden, Germany) in 600 µL 0.5M KH2PO4-buffer supplemented with protease inhibitor (cOmplete, Roche) pH7 for 5min at 50 Hz. 20 µL of the homogenate and 2.5 µL of serum were subjected to the enzymatic assay.
Lactate was measured following the protocol from Lin and colleagues [43]. In brief, the assay uses coupling of lactate-oxidase and peroxidase for the conversion of 2,2′-azino-bis (3-ethylbenzothiazoline-6-sulfonic acid) (ABTS, Sigma Aldrich, Steinheim, Germany). As a standard, lithium-lactate was used, and the absorbance of the chromogenic product was measured at 405 nm. The protein content of all samples was assessed with Roti-Nanoquant (Roth) and used for normalization.

2.4. RNA-Seq Library Preparation and Sequencing

Next-generation sequencing library preparation was performed with Illumina’s TruSeq stranded mRNA LT Sample Prep Kit following Illumina’s standard protocol (Part # 15,031,047 Rev, E). The libraries were prepared with a starting amount of 302 ng and amplified in 12 PCR cycles, profiled in a DNA 1000 Chip on a 2100 Bioanalyzer (Agilent Technologies) and quantified using the Qubit dsDNA HS Assay Kit, in a Qubit 2,0 Fluorometer (Life technologies). All 24 samples were pooled in equimolar ratio and sequenced on 1 NextSeq 500 Highoutput FC, PE for 2 × 42 cycles plus 7 cycles for the index read.

2.5. RNA-Seq Analysis

After merging five technical replicates and down-sampling one sequencing file, the reads were trimmed using BBDuk (version 38.06) [44]. For the differential expression analysis, the trimmed reads were mapped to the UCSC reference genome of Mus musculus obtained through Illumina iGenomes (version mm10, latest release in May 2012) [45], using the splice-aware mapper STAR (version 2.7.0d, default options of 2-pass mapping) [46]. The aligned reads were then counted per gene and sample (based on the mm10 UCSC annotation file from iGenomes) using FeatureCounts provided by SubRead (version 1.6.2, default options for single-end reads) [47].
The differential expression analysis (DEA) was done in R (version 4.0.3) and RStudio (version 1.3.1073) using the R package DESeq2 (version 1.30.0) [48,49,50]. In order to reduce the impact of extreme expressions, three genes, Psen1, App, and Thy1, which are known to be highly differentially expressed due to the genotype of 5xFAD mice [51], were removed from the count tables. The normalized gene counts of these marker genes can be found in Supplementary Figure S1. After normalizing the filtered gene counts by the size factors provided by DESeq2, a principal component analysis was performed for all mice and each genotype separately [52]. Additionally, PCA was conducted based on the separate training groups and is shown in the Supplementary Figure S2. The p-values of the DEA were adjusted using the Benjamini–Hochberg method, and the cut-off for significant differential regulation was set at adjusted p-value < 0.05 [53].
Gene Ontology (GO) term enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed for all subsets of differentially regulated genes using the R package clusterProfiler (version 3.18.0) [54]. Furthermore, the gene subsets were examined for brain cell-type enrichment using the Fisher’s Exact Test implemented in R [55]. The cell-type-associated gene list consisted of five brain cell types, including neurons, astrocytes, oligodendrocytes, microglia, and endothelial cells, and is based on a previously published brain transcription analysis [56,57]. The gene universe (i.e., the gene’s background set for the statistical analysis) was set to the genes considered for the differential expression analysis for all tests. After adjusting the p-values of the overrepresentation test using the Benjamini–Hochberg method, the terms with an adjusted p-value < 0.05 were considered as significantly enriched.
All plots were created using the R packages ggplot2 (version 3.3.2) and venn (version 1.9) [58,59]. The raw fastq files and the count table were uploaded to Gene Expression Omnibus (GEO) under the accession number GSE164798 (see Data Availability Statement).

3. Results

3.1. Voluntary Running Wheel Usage in 5xFAD Mice

Physical activity has often been discussed as beneficial in AD or animal models of the disorder. To analyze the potential impact of acute and chronic physical voluntary training, we investigated male mice with the 5xFAD genotype and their wild-type littermates, starting at the age of 2 months. At this age, no behavioral deficits have been reported in the mice; however, deposition of Aβ is observable already at 1.5 months of age [51]. Mice were single-caged and divided into three groups: an untrained group with no access to a saucer wheel, a chronically trained group with full access over the 30 days of the experiment, and the third group with access only within the last night before sacrifice (see Figure 1A). Running activity was assessed by a wireless reporting system: within the first night after adding the saucer wheel, animals of both genotypes quickly adjusted to the new enrichment during the first hours after being exposed to the saucer wheel. This is indicated by a significant count increase in both genotypes starting at 22 p.m. (comparison vs. starting point at 10 a.m., Figure 1B). No apparent differences occurred in the starting phase in both genotypes despite slight but non-significantly higher arousal in 5xFAD mice. Approximately 150 counts for wheel turning per hour were reached at 12 a.m. With prolonged access to the saucer wheel, counts as high as 2800 per hour were reached on average (Figure 1C) with a clear peak in the early dark phase (between 6 and 10 p.m). 5xFAD mice showed a statistically significant elevation of wheel usage at several points of time, resulting in an increased mean running distance (5xFAD: 8.6 ± 0.32 km per night, wild-type: 7.5 ± 0.41 km per night) but no change in the circadian rhythm of usage.

3.2. Physical Parameters of 5xFAD Mice under Acute and Chronic Saucer Wheel Usage

Physical activity might affect general parameters such as body weight or feeding behavior (caloric intake) that must be considered. Therefore, body weight was measured twice weekly. As depicted in Figure 2A and reported previously [60], male 5xFAD mice, in general, had a lower body weight as compared to wild-type controls at the age of 2 months but also at the end of the experiment. However, body weight did not change due to the respective activity group. Food consumption was elevated in both chronic training groups independently of the genotype when comparing with the starting point and the end of the experiment (Figure 2B). Interestingly, the increase in food intake was smaller in chronically trained 5xFAD males than in their wild-type littermates. Within the last night, the acute training phase led to a slight non-significant increase in food consumption in both groups.
Next, we aimed at quantifying parameters that characterize the training effect in the three groups of mice. For this purpose, blood glucose levels were measured from truncal blood upon sacrifice and the weight of the Gastrocnemius muscle, the heart, and abdominal fat mass was determined (Figure 3). Animals of all three groups were sacrificed in an alternating order to prevent circadian effects. As expected, we observed a reduction in fat mass with regular training (as seen in Figure 3D) and an increase in heart mass (Figure 3B). These parameters were comparable for both genotypes even though the heart mass elevation was not statistically significant in transgenic mice.
More astonishing was the effect of acute saucer wheel training on 5xFAD mice: in all examined parameters, the acutely trained group of 5xFAD males differed from their wild-type littermates. For example, the increase in blood glucose levels that was only subtly observable in wild-type mice reached statistical significance in transgenic animals (Figure 3A). The heart mass was reduced in the acutely trained transgenic animals compared to controls (Figure 3B). Moreover, we measured significantly lower abdominal fat mass amounts relative to the wild-type animals after acute training (Figure 3D). All these findings lead to the assumption that 5xFAD mice differ metabolically from their wild-type littermates.
An important metabolite for brain homeostasis that has only gained attention within the last few years is lactate. Both genotypes were indistinguishable in terms of serum lactate (Figure 4A). However, the lactate level of brain homogenates was quite striking in the AD mouse model: while chronic activity displayed no effect at all as compared to untrained mice, the acute physical activity for 18 h resulted in doubled lactate levels (Figure 4B), indicating that the physical activity had a direct effect not only on muscle or heart weight but also on brain metabolism.

3.3. RNA-Seq Analysis of 5xFAD and Wild-Type Mice with Chronic and Acute Physical Activity

To unravel if the pronounced impact of acute training on 5xFAD mice also corresponded to transcriptomic changes, we employed high-throughput sequencing of RNA extracted from the preserved brain hemispheres of the transgenic and wild-type mice.
After normalization and filtering of the aligned RNA-seq reads by DESeq2, 20,132 of 24,418 genes returned a non-zero total read count. First, we performed a principal component analysis (PCA) [52] to determine whether the highest variances within the gene counts were related to genotype and training differences. In the PCA, including all samples, principal component 1 (PC1) captured 70% of the variance in the data (Figure 5A). Along PC2, which accounted for 9% of the variance, wild-type mice were partially separated from the 5xFAD mice, indicating genotype-associated differences. Within the 5xFAD mice, the training groups did not cluster together based on the first and second PC and the first PC already captured 90% of the variance (Figure 5B). This indicated that the effect of training on the brain transcriptomic profile was relatively small in comparison to general biological differences in 5xFAD mice. We observed the same trend for untrained and chronically trained wild-type mice, whereas all acutely trained wild-type mice were clearly separated from animals exposed to chronic or no training.

3.3.1. Differential Expression Analysis between Different Physical Activity Levels in 5xFAD and Wild-Type Mice

In order to investigate the impact of different levels of physical activity on RNA expression regulation, we compared the transcriptomes of the different activity groups in either 5xFAD or wild-type mice. Interestingly, no transcriptional differences could be detected among the pairwise combinations of sedentary, acutely and chronically trained 5xFAD mice using an adjusted p-value (p-adj.) < 0.05 except for one gene—Secretogranin-1 (Chgb) (Supplementary Figure S3). Chgb is a neuroendocrine secretory granule protein that was significantly downregulated in 5xFAD mice with one day of wheel running compared to those with no training (p-adj. = 0.02, log2 fold change (LFC) = −0.275). Furthermore, the difference was almost statistically significant comparing acutely trained vs. untrained mice (p-adj. = 0.093, LFC = −0.258). Among the wild-type mice, the acutely trained animals displayed an altered expression pattern compared to both sedentary and chronically trained mice with 835 and 1617 differentially expressed genes, respectively (DEGs) (Figure 6A–D; Supplementary Table S1). Out of these genes, 595 were commonly up- or downregulated in 4-week trained and untrained mice compared to wild-type animals with one day of wheel running (Figure 6D). These results indicated that acute or persistent physical activity did not alter the basic gene regulation of transgenic mice. In contrast, acute activity was associated with a pronounced effect on gene regulation in wild-type mice.

3.3.2. Functional Annotation of the Altered Gene Expression of Acutely Trained Wild-Type Mice

To uncover the biological background of the transcriptional differences in wild-type mice possibly evoked by one day of voluntary wheel running, we applied the Gene Ontology (GO) term, KEGG pathway, and cell-type enrichment analyses to the significantly up- and downregulated genes that we found in either acute vs. sedentary or acute vs. chronic physical activity. All terms or cells with p-adj. < 0.05 were considered as significant based on the statistical overrepresentation test.
The genes that were upregulated in acutely trained wild-type mice were significantly enriched for GO terms mainly involved in Biological Processes of DNA metabolism and transcriptional regulation like “DNA repair” (p-adj. = 1.42 × 104), “histone modification” (p-adj. = 2.16 × 102), and “mRNA splicing, via spliceosome” (p-adj. = 1.22 × 104) (Figure 6E; Supplementary Table S2). On the other hand, many enriched terms of the downregulated genes were involved in general neuronal mechanisms including “axonogenesis” (p-adj. = 3.46 × 108), “post-synapse assembly” (p-adj. = 9.6 × 106) and “regulation of neurotransmitter receptor activity” (p-adj. = 8.42 × 106) (Figure 6E; Supplementary Table S2). Interestingly, the KEGG pathway analysis of all acute-specific genes resulted in a significant enrichment of terms related to neurodegenerative diseases with a p-adj. < 0.05. Additionally, the “Spliceosome” pathway was also overrepresented in the acute-specific genes (p-adj. = 4.71 × 102) (Supplementary Table S2). The cell-type enrichment analysis supported the GO term, and KEGG pathway results as the acute-specific downregulated genes were significantly enriched for neuronal cells (p-adj. = 1.32 × 104) (Figure 6F, Supplementary Table S2). This finding indicated that the short-term wheel running of wild-type mice might induce a neuronal stress state in the brain, altering the homeostasis of neuron functions and possibly overregulating basic DNA and RNA mechanistic processes in parallel.

3.3.3. Differential Expression Analysis between 5xFAD and Wild-Type Mice after Different Physical Activity Levels

Next, we investigated the differences in gene expression between 5xFAD and wild-type animals at each training duration level. The DEA between untrained 5xFAD and wild-type mice revealed 141 DEGs (p-adj. < 0.05) (Figure 7A; Supplementary Table S3). Since environmental factors did not influence these mice, the detected gene expression differences were likely to be disease-related. Regarding the transcriptional patterns following chronic training, the comparison between the genotypes yielded the most considerable difference in expression with 477 DEGs. Notably, more than 80% of these genes were associated with a significantly lower expression in 5xFAD mice (Figure 7B, Supplementary Table S3). Eighty genes were commonly regulated, comparing 5xFAD with wild-type mice after no training and 30 days of voluntary physical activity (Figure 7D).
Acute training resulted in 14 DEGs between 5xFAD and wild-type mice, of which 13 were upregulated (Figure 7C; Supplementary Table S3). All upregulated genes found after acute activity showed a significant upregulation in 5xFAD mice after chronic and/or no training compared to wild-type animals (Figure 7D). Moreover, nine of these 13 upregulated genes (Cst7, Ccl6, Itgax, Gfap, Clec7a, Ccp4, Cd68, Ccl3, and Trem2) are already known to be differentially expressed in the 5xFAD mouse model [61]. After one day of wheel running, the only downregulated gene was again Chgb, the same gene we previously found to be downregulated when comparing acutely trained with chronically trained 5xFAD mice.

3.3.4. Functional Annotation of Genotype-Specific Genes Found after Chronic and No Training

The genes that were upregulated in the 5xFAD model without any training as compared to wild-type mice were mainly enriched for GO terms of immune-related processes (Supplementary Table S4). The same was true for the upregulated genes found in chronically trained 5xFAD compared to wild-type animals as well as for the shared upregulated genes of the sedentary and chronically trained transgenic vs. wild-type mice (Supplementary Table S4). Moreover, the microglia cell type was significantly enriched in the upregulated genes in 5xFAD compared to wild-type mice after chronic training and especially without training (Figure 7E, Supplementary Table S4).

4. Discussion

Our study demonstrated that brain lactate levels increased dramatically after acute voluntary wheel running (1 day) only in 5xFAD mice but not in wild-type animals. In contrast, plasma lactate levels remained unchanged in all animals after physical activity. Interestingly, the whole brain’s subsequent molecular analysis revealed that the transcriptional patterns of 5xFAD mice after acute, chronic, and no training were very similar. In contrast, the acutely trained wild-type mice showed a significantly altered transcriptome compared with both groups either receiving no access to the wheel or 30 days of voluntary wheel running.

4.1. Brain Lactate Was Drastically Increased in 5xFAD Mice after Acute Physical Activity

Lactate is a metabolite released from muscle cells during and after intense physical activity [62]. Its conversion from pyruvate by lactate dehydrogenase is facilitated under conditions of reduced oxygen availability, so-called aerobic glycolysis [63,64,65]. As a result, after physical exercise, there is an acute time-restricted increase in lactate levels in circulating blood, which depends on the intensity of physical activity and individual cardiovascular fitness [66,67,68]. In addition to its role in replenishing pyruvate in the liver to provide substrates for gluconeogenesis, lactate is preferentially consumed and metabolized by neurons in an activity-dependent manner [69,70]. Several studies on rodent models have shown that neuronal lactate uptake plays an essential role in long-term memory functions and neurogenesis [71,72]. As an additional fuel of the brain, lactate is proposed as one of the main actors mediating the beneficial effects of physical activity on cognition.
In contrast to our findings, El Hayek and colleagues reported elevated lactate concentrations in healthy mice’s hippocampus after 30 days of voluntary wheel running, which was associated with cognitive improvements [73]. Even though we used the same training regime and the same mouse strain, a possible reason for the non-rising lactate levels after chronic training in our study could be that we measured the lactate concentration in whole-brain hemisphere tissue rather than in a specific region. Moreover, the mice in El Hayek’s study were one month younger than the mice in our analysis. On the one hand, this may indicate that regular physical activity induces region-specific increases in lactate levels, while the brain’s overall concentration remains unchanged. On the other hand, the animals’ age may also play an essential role in lactate enrichment due to physical effort.
Furthermore, we observed increased lactate brain levels in acutely trained 5xFAD animals but not in their wild-type counterparts. Therefore, this finding could indicate an overreaction to a possible neuronal activation and increased energy demands during physical activity in the AD model mice in order to compensate for a disease-driven energy deficit, which is common in AD [74]. The overreaction could be further enhanced by stress or arousal caused by the new motor stimulation of one-day trained animals. Previous studies have shown that rodents exposed to a new environment or behavioral testing exhibit increased aerobic glycolysis leading to higher hippocampal lactate levels and improved memory functions and gene expression changes [64]. Whether this reaction is enhanced under pathological conditions has not been investigated yet.
Furthermore, our results suggest that the lactate increase observed in acutely trained 5xFAD mice was mainly attributable to synthesis within the brain rather than transport from the blood. One possible lactate production mechanism within the brain is associated with the astrocyte-neuron-lactate shuttle (ANLS). This model hypothesizes that neurons release glutamate during neuronal activation, transport it to astrocytes that generate lactate, which is then transported back to neurons and used as an energy source for neuronal functions [64]. In line with this, Zuend and colleagues recently reported that lactate is released by astrocytes and taken up by neurons due to acute environmental arousal in mice [75]. In contrast, Díaz-García and colleagues provided evidence that acutely excited neurons perform aerobic glycolysis by themselves rather than taking up lactate produced by astrocytes [76,77]. Unfortunately, our experimental design does not allow us to ascertain if the lactate source in our study was attributable to astrocyte or direct neuronal production.

4.2. One Day of Voluntary Wheel Running Altered Transcription Only in Wild-Type but Not in 5xFAD Mice

Multiple studies investigating acute or chronic training were able to show beneficial effects on cognition, Aβ plaques, brain size, immune functions, or expression of neurotrophins in different transgenic mouse models of Alzheimer’s disease [18,22,23,24]. However, other studies could not confirm the positive influence of physical exercise on AD [35,36,37,38,39]. Recently, the 5xFAD mouse model was investigated after chronic voluntary wheel training, and neither cognitive improvements nor decreased Aβ levels and no reduction in AD-driven neuroinflammation were detected [37,38]. This is in accordance with our findings at the transcriptome level. One month of voluntary wheel training did not change the global gene expression in 5xFAD mice compared to sedentary animals. However, as we did not consider behavior or immune-related functions, we cannot conclude that regular activity cannot improve the pathology of the 5xFAD mouse model. Furthermore, other biological layers like the metabolome or proteome could be positively influenced by acute and chronic activity, while the transcriptome might play a minor role.
Additionally, the secretory-related gene Chgb was significantly downregulated in acutely trained 5xFAD mice compared to chronically trained and sedentary transgenic mice. The deposition of Aβ has previously been reported to be related to the downregulation of genes involved in the secretory pathway, including Chgb [78]. The slight reduction in Chgb expression potentially indicates that the short-term activity might enhance the AD-like decline by downregulating important secretory pathway factors. However, this finding contradicts the potentially beneficial impact of the lactate increases on the neuronal function within the acutely trained 5xFAD group. As physical activity is suggested to be more of a preventive strategy against cognitive decline, it is possible that the disease phenotype of the aggressive 5xFAD mouse model is already too advanced at two months of age to be alleviated by an external intervention such as physical activity.
Regarding the physical activity duration in wild-type mice, only one day of voluntary wheel running altered the transcriptional pattern as compared to sedentary mice. In humans, one single bout of exercise is suggested to positively influence cognition as well as mood and emotional state to a small extent [34]. Several studies investigating neurochemical measurements after one exercise session have shown that specific neurotransmitters, neurotrophins, and other essential metabolites related to neurogenesis and brain plasticity are enhanced after this short time [34]. In our study, DEA revealed that axons’ growth and development within acutely trained wild-type mice were likely inhibited. Many under-expressed genes were associated with axonogenesis assuming that these genes are positively correlated to this process. However, human studies have also reported negative effects after acute exercise [34,79]. For instance, acute stress and acute exercise share several mechanisms like the activation of the hypothalamic–pituitary–adrenal (HPA) axis and the sympathetic nervous system, both related to immunity.
In contrast to regular physical activity, acute exercise has been shown to reduce immunity and increase susceptibility to infections [80]. Furthermore, elevated oxidative stress levels have been correlated with acute exercise, which is involved in functional alterations, especially in the central nervous system [79]. Thus, changes in the transcriptome of the acutely trained wild-type mice in our study might correspond to a transient response to acute stress. In contrast, regular physical activity leads to an adaptation and a return of the transcriptome to a homeostatic state.

4.3. Immune-Related Genes Were Upregulated in Untrained and Chronically Trained 5xFAD Mice Compared to Their Wild-Type Counterparts

Furthermore, we investigated activity-genotype specific transcriptional regulation by comparing 5xFAD with wild-type mice after one day, 30 days, or no voluntary wheel training. Importantly, sedentary and chronically trained animals revealed several disease-driven gene expression differences. Here, we could demonstrate that the over-expressed genes in chronically trained or untrained 5xFAD mice compared to wild-type littermates were mainly involved in immune response processes. Neuroinflammation is one of the significant hallmarks of AD and has recently been suggested to be the main driver of the early progression of AD [81]. Within the 5xFAD model, immune-system changes are likely to occur prior to plaque deposition development in 2-month-old 5xFAD mice [82]. In addition, we found that the genes upregulated in 5xFAD mice were significantly associated with microglial cells. Microglia are responsible for activating the brain’s innate immune system, and their accumulation contributes to neurodegenerative pathologies [40,83,84]. Interestingly, the comparison of the genotypes after regular training yielded a much larger transcriptional gap between 5xFAD and wild-type animals than in sedentary mice, and the majority of genes were downregulated. However, no biological function was overrepresented for these under-expressed genes.
In contrast, acutely trained wild-type mice, which differed from the other wild-type groups, had a similar transcriptomic profile to their 5xFAD counterparts. In agreement with our findings, a human study showed that the human ortholog of Chgb, the only downregulated gene in acutely trained 5xFAD compared to wild-type mice, was significantly downregulated in AD patients compared to healthy controls [85]. Notably, many immune-response-associated genes were no longer differentially regulated between the transgenic and wild-type mice, as was the case for the chronically trained or sedentary animals. These findings corroborate the assumption that acute physical activity triggered a stress response in wild-type animals, which shifted their transcriptome profile in the direction of their transgenic counterparts.
Taken together, our results suggest that the brain’s global transcriptome might not be the primary mediator of the beneficial effects found in healthy and AD individuals due to regular physical activity. Interestingly, acute training was associated with drastic changes in both 5xFAD and wild-type mice, however, in entirely different ways. Brain lactate concentrations and other physiological parameters were altered in 5xFAD mice. In contrast, the brain transcriptome of wild-type animals changed its expression patterns after acute training. These genotype-specific differences should be further investigated in order to elucidate the real link between acute physical activity and AD.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4409/10/3/693/s1, Figure S1: Normalized gene counts of mutated and mutation-correlated genes that produce the 5xFAD mouse model, Figure S2: Principal Component Analysis of normalized counts from training-specific groups, Figure S3: Volcano plots of all pairwise contrasts of training groups in 5xFAD samples. Table S1: Differential expression results from the pairwise comparisons between acute and chronic training and between acute and no training of wild-type mice. Table S2: Functional annotation results of differential expressed genes observed from the comparisons of different activity regimens in wild-type mice, Table S3: Differential expression results from the pairwise comparisons between 5xFAD and wild-type mice after no training, chronic training and acute training, Table S4: Functional annotation results of differential expressed genes observed for the comparisons between 5xFAD and wild-type mice after each training duration.

Author Contributions

Data acquisition and curation: L.G.; data analysis: A.W., H.T., L.G. and V.T.T.N.; writing—original draft preparation, A.W.; writing—review and editing: H.T., S.G. and K.E.; interpretation: A.W., K.E., S.G. and H.T.; assistance with technical equipment: K.R.; supervision: S.G. and K.E. All authors have read and agreed to the published version of the manuscript.

Funding

This work was largely funded by the NMFZ-grant “Transcriptional and Post-Transcriptional Causes of Physiological Aging”. H.T. and S.G. acknowledge funding from the ReALity initiative. A.W. acknowledges funding from the Leibniz Institute for Resilience Research (LIR) gGmbH and the IDSAIR initiative. The funders had no role in the research design, data collection and analysis, writing of the manuscript and the decision to publish.

Institutional Review Board Statement

All procedures were performed in accordance with the European Communities Council Directive regarding care and use of animals for experimental procedures and were approved by local authorities (Landesuntersuchungsamt Rheinland-Pfalz; approval number G 14-1-087, date of approval: 28 July 2015).

Data Availability Statement

The raw RNA-seq data and the corresponding feature count table discussed in this publication have been deposited in NCBI’s Gene Expression Omnibus [86] and are accessible through GEO Series accession number GSE164798 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE164798).

Acknowledgments

Support by the IMB Genomics Core Facility for the RNA sequencing is gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Roth, G.A.; Abate, D.; Abate, K.H.; Abay, S.M.; Abbafati, C.; Abbasi, N.; Abbastabar, H.; Abd-Allah, F.; Abdela, J.; Abdelalim, A.; et al. Global, Regional, and National Age-Sex-Specific Mortality for 282 Causes of Death in 195 Countries and Territories, 1980–2017: A Systematic Analysis for the Global Burden of Disease Study 2017. Lancet 2018, 392, 1736–1788. [Google Scholar] [CrossRef] [Green Version]
  2. Conde-Sala, J.L.; Turró-Garriga, O.; Piñán-Hernández, S.; Portellano-Ortiz, C.; Viñas-Diez, V.; Gascón-Bayarri, J.; Reñé-Ramírez, R. Effects of Anosognosia and Neuropsychiatric Symptoms on the Quality of Life of Patients with Alzheimer’s Disease: A 24-Month Follow-up Study. Int. J. Geriatr. Psychiatry 2016, 31, 109–119. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. McKhann, G.; Drachman, D.; Folstein, M.; Katzman, R.; Price, D.; Stadlan, E.M. Clinical Diagnosis of Alzheimer’s Disease: Report of the NINCDS-ADRDA Work Group under the Auspices of Department of Health and Human Services Task Force on Alzheimer’s Disease. Neurology 1984, 34, 939–944. [Google Scholar] [CrossRef] [Green Version]
  4. Luque-Contreras, D.; Carvajal, K.; Toral-Rios, D.; Franco-Bocanegra, D.; Campos-Peña, V. Oxidative Stress and Metabolic Syndrome: Cause or Consequence of Alzheimer’s Disease? Oxidative Med. Cell. Longev. 2014, 2014, 1–11. [Google Scholar] [CrossRef]
  5. Calsolaro, V.; Edison, P. Neuroinflammation in Alzheimer’s Disease: Current Evidence and Future Directions. Alzheimer’s Dement. 2016, 12, 719–732. [Google Scholar] [CrossRef]
  6. Hardy, J.A.; Higgins, G.A. Alzheimer’s Disease: The Amyloid Cascade Hypothesis. Science 1992, 256, 184–185. [Google Scholar] [CrossRef] [PubMed]
  7. Arnold, S.E.; Hyman, B.T.; Flory, J.; Damasio, A.R.; Van Hoesen, G.W. The Topographical and Neuroanatomical Distribution of Neurofibrillary Tangles and Neuritic Plaques in the Cerebral Cortex of Patients with Alzheimer’s Disease. Cereb. Cortex 1991, 1, 103–116. [Google Scholar] [CrossRef]
  8. de la Monte, S.M. Insulin Resistance and Neurodegeneration: Progress Towards the Development of New Therapeutics for Alzheimer’s Disease. Drugs 2017, 77, 47–65. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Hinz, F.I.; Geschwind, D.H. Molecular Genetics of Neurodegenerative Dementias. Cold Spring Harb. Perspect. Biol. 2017, 9, a023705. [Google Scholar] [CrossRef] [Green Version]
  10. Bellenguez, C.; Grenier-Boley, B.; Lambert, J.C. Genetics of Alzheimer’s Disease: Where We Are, and Where We Are Going. Curr. Opin. Neurobiol. 2020, 61, 40–48. [Google Scholar] [CrossRef]
  11. Ruffini, N.; Klingenberg, S.; Schweiger, S.; Gerber, S. Common Factors in Neurodegeneration: A Meta-Study Revealing Shared Patterns on a Multi-Omics Scale. Cells 2020, 9, 2642. [Google Scholar] [CrossRef]
  12. Hewel, C.; Kaiser, J.; Wierczeiko, A.; Linke, J.; Reinhardt, C.; Endres, K.; Gerber, S. Common MiRNA Patterns of Alzheimer’s Disease and Parkinson’s Disease and Their Putative Impact on Commensal Gut Microbiota. Front. Neurosci. 2019, 13, 113. [Google Scholar] [CrossRef] [Green Version]
  13. Zhang, L.; Tang, W.; Chao, F.; Zhou, C.; Jiang, L.; Zhang, Y.; Liang, X.; Tang, J.; Qi, Y.; Yang, H.; et al. Four-Month Treadmill Exercise Prevents the Decline in Spatial Learning and Memory Abilities and the Loss of Spinophilin-Immunoreactive Puncta in the Hippocampus of APP/PS1 Transgenic Mice. Neurobiol. Dis. 2020, 136, 104723. [Google Scholar] [CrossRef] [PubMed]
  14. Fagan, A.M.; Head, D.; Shah, A.R.; Marcus, D.; Mintun, M.; Morris, J.C.; Holtzman, D.M. Decreased Cerebrospinal Fluid Aβ42 Correlates with Brain Atrophy in Cognitively Normal Elderly. Ann. Neurol. 2009, 65, 176–183. [Google Scholar] [CrossRef] [Green Version]
  15. Deming, Y.; Li, Z.; Kapoor, M.; Harari, O.; Del-Aguila, J.L.; Black, K.; Carrell, D.; Cai, Y.; Fernandez, M.V.; Budde, J.; et al. Genome-Wide Association Study Identifies Four Novel Loci Associated with Alzheimer’s Endophenotypes and Disease Modifiers. Acta Neuropathol. 2017, 133, 839–856. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Morris, J.C.; Price, J.L. Pathologic Correlates of Nondemented Aging, Mild Cognitive Impairment, and Early-Stage Alzheimer’s Disease. J. Mol. Neurosci. 2001, 17, 101–118. [Google Scholar] [CrossRef]
  17. Jia, R.; Liang, J.; Xu, Y.; Wang, Y. Effects of Physical Activity and Exercise on the Cognitive Function of Patients with Alzheimer Disease: A Meta-Analysis. BMC Geriatr. 2019, 19, 181. [Google Scholar] [CrossRef]
  18. De la Rosa, A.; Olaso-Gonzalez, G.; Arc-Chagnaud, C.; Millan, F.; Salvador-Pascual, A.; García-Lucerga, C.; Blasco-Lafarga, C.; Garcia-Dominguez, E.; Carretero, A.; Correas, A.G.; et al. Physical Exercise in the Prevention and Treatment of Alzheimer’s Disease. J. Sport Health Sci. 2020, 9, 394–404. [Google Scholar] [CrossRef] [PubMed]
  19. Orgeta, V.; Mukadam, N.; Sommerlad, A.; Livingston, G. The Lancet Commission on Dementia Prevention, Intervention, and Care: A Call for Action. Ir. J. Psychol. Med. 2019, 36, 85–88. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Rashid, M.H.; Zahid, M.F.; Zain, S.; Kabir, A.; Hassan, S.U. The Neuroprotective Effects of Exercise on Cognitive Decline: A Preventive Approach to Alzheimer Disease. Cureus 2020, 12, e6958. [Google Scholar] [CrossRef] [Green Version]
  21. Roig, M.; Nordbrandt, S.; Geertsen, S.S.; Nielsen, J.B. The Effects of Cardiovascular Exercise on Human Memory: A Review with Meta-Analysis. Neurosci. Biobehav. Rev. 2013, 37, 1645–1666. [Google Scholar] [CrossRef]
  22. Erickson, K.I.; Hillman, C.; Stillman, C.M.; Ballard, R.M.; Bloodgood, B.; Conroy, D.E.; Macko, R.; Marquez, D.X.; Petruzzello, S.J.; Powell, K.E. Physical Activity, Cognition, and Brain Outcomes: A Review of the 2018 Physical Activity Guidelines. Med. Sci. Sports Exerc. 2019, 51, 1242–1251. [Google Scholar] [CrossRef]
  23. Panza, G.A.; Taylor, B.A.; MacDonald, H.V.; Johnson, B.T.; Zaleski, A.L.; Livingston, J.; Thompson, P.D.; Pescatello, L.S. Can Exercise Improve Cognitive Symptoms of Alzheimer’s Disease? J. Am. Geriatr. Soc. 2018, 66, 487–495. [Google Scholar] [CrossRef] [Green Version]
  24. Brown, B.M.; Peiffer, J.J.; Martins, R.N. Multiple Effects of Physical Activity on Molecular and Cognitive Signs of Brain Aging: Can Exercise Slow Neurodegeneration and Delay Alzheimer’s Disease? Mol. Psychiatry 2013, 18, 864–874. [Google Scholar] [CrossRef] [Green Version]
  25. Sellami, M.; Gasmi, M.; Denham, J.; Hayes, L.D.; Stratton, D.; Padulo, J.; Bragazzi, N. Effects of Acute and Chronic Exercise on Immunological Parameters in the Elderly Aged: Can Physical Activity Counteract the Effects of Aging? Front. Immunol. 2018, 9, 2187. [Google Scholar] [CrossRef] [Green Version]
  26. Segal, S.K.; Cotman, C.W.; Cahill, L.F. Exercise-Induced Noradrenergic Activation Enhances Memory Consolidation in Both Normal Aging and Patients with Amnestic Mild Cognitive Impairment. J. Alzheimer’s Dis. 2012, 32, 1011–1018. [Google Scholar] [CrossRef] [Green Version]
  27. Von Ah Morano, A.E.; Dorneles, G.P.; Peres, A.; Lira, F.S. The Role of Glucose Homeostasis on Immune Function in Response to Exercise: The Impact of Low or Higher Energetic Conditions. J. Cell. Physiol. 2020, 235, 3169–3188. [Google Scholar] [CrossRef]
  28. Adlard, P.A.; Perreau, V.M.; Pop, V.; Cotman, C.W. Voluntary Exercise Decreases Amyloid Load in a Transgenic Model of Alzheimer’s Disease. J. Neurosci. 2005, 25, 4217–4221. [Google Scholar] [CrossRef]
  29. Trigiani, L.J.; Lacalle-Aurioles, M.; Bourourou, M.; Li, L.; Greenhalgh, A.D.; Zarruk, J.G.; David, S.; Fehlings, M.G.; Hamel, E. Benefits of Physical Exercise on Cognition and Glial White Matter Pathology in a Mouse Model of Vascular Cognitive Impairment and Dementia. Glia 2020, 68, 1925–1940. [Google Scholar] [CrossRef]
  30. Liu, Y.; Chu, J.M.T.; Yan, T.; Zhang, Y.; Chen, Y.; Chang, R.C.C.; Wong, G.T.C. Short-Term Resistance Exercise Inhibits Neuroinflammation and Attenuates Neuropathological Changes in 3xTg Alzheimer’s Disease Mice. J. Neuroinflammation 2020, 17, 4. [Google Scholar] [CrossRef] [Green Version]
  31. Belarbi, K.; Burnouf, S.; Fernandez-Gomez, F.J.; Laurent, C.; Lestavel, S.; Figeac, M.; Sultan, A.; Troquier, L.; Leboucher, A.; Caillierez, R.; et al. Beneficial Effects of Exercise in a Transgenic Mouse Model of Alzheimer’s Disease-like Tau Pathology. Neurobiol. Dis. 2011, 43, 486–494. [Google Scholar] [CrossRef]
  32. Ryan, S.M.; Kelly, Á.M. Exercise as a Pro-Cognitive, pro-Neurogenic and Anti-Inflammatory Intervention in Transgenic Mouse Models of Alzheimer’s Disease. Ageing Res. Rev. 2016, 27, 77–92. [Google Scholar] [CrossRef] [Green Version]
  33. Lu, X.; Moeini, M.; Li, B.; de Montgolfier, O.; Lu, Y.; Bélanger, S.; Thorin, É.; Lesage, F. Voluntary Exercise Increases Brain Tissue Oxygenation and Spatially Homogenizes Oxygen Delivery in a Mouse Model of Alzheimer’s Disease. Neurobiol. Aging 2020, 88, 11–23. [Google Scholar] [CrossRef]
  34. Basso, J.C.; Suzuki, W.A. The Effects of Acute Exercise on Mood, Cognition, Neurophysiology, and Neurochemical Pathways: A Review. Brain Plast. 2017, 2, 127–152. [Google Scholar] [CrossRef] [Green Version]
  35. Madden, D.J.; Blumenthal, J.A.; Allen, P.A.; Emery, C.F. Improving Aerobic Capacity in Healthy Older Adults Does Not Necessarily Lead to Improved Cognitive Performance. Psychol. Aging 1989, 4, 307–320. [Google Scholar] [CrossRef]
  36. Toots, A.; Littbrand, H.; Boström, G.; Hornsten, C.; Holmberg, H.; Lillemor, L.O.; Lindelöf, N.; Nordstrom, P.; Gustafson, Y.; Rosendahl, E. Effects of Exercise on Cognitive Function in Older People with Dementia: A Randomized Controlled Trial. J. Alzheimer’s Dis. 2017, 60, 323–332. [Google Scholar] [CrossRef] [Green Version]
  37. Hansson, O.; Svensson, M.; Gustavsson, A.M.; Andersson, E.; Yang, Y.; Nägga, K.; Hållmarker, U.; James, S.; Deierborg, T. Midlife Physical Activity Is Associated with Lower Incidence of Vascular Dementia but Not Alzheimer’s Disease. Alzheimer’s Res. Ther. 2019, 11, 87. [Google Scholar] [CrossRef]
  38. Svensson, M.; Andersson, E.; Manouchehrian, O.; Yang, Y.; Deierborg, T. Voluntary Running Does Not Reduce Neuroinflammation or Improve Non-Cognitive Behavior in the 5xFAD Mouse Model of Alzheimer’s Disease. Sci. Rep. 2020, 10, 1–10. [Google Scholar] [CrossRef] [Green Version]
  39. Frederiksen, K.S.; Gjerum, L.; Waldemar, G.; Hasselbalch, S.G. Effects of Physical Exercise on Alzheimer’s Disease Biomarkers: A Systematic Review of Intervention Studies. J. Alzheimer’s Dis. 2018, 61, 359–372. [Google Scholar] [CrossRef]
  40. Valenzuela, P.L.; Castillo-García, A.; Morales, J.S.; de la Villa, P.; Hampel, H.; Emanuele, E.; Lista, S.; Lucia, A. Exercise Benefits on Alzheimer’s Disease: State-of-the-Science. Ageing Res. Rev. 2020, 62, 101108. [Google Scholar] [CrossRef]
  41. Lista, I.; Sorrentino, G. Biological Mechanisms of Physical Activity in Preventing Cognitive Decline. Cell. Mol. Neurobiol. 2010, 30, 493–503. [Google Scholar] [CrossRef]
  42. Archer, T. Physical Exercise Alleviates Debilities of Normal Aging and Alzheimer’s Disease. Acta Neurol. Scand. 2011, 123, 221–238. [Google Scholar] [CrossRef]
  43. Lin, C.-Y.; Chen, S.-H.; Kou, G.-H.; Kuo, C.-M. An Enzymatic Microassay for Lactate Concentration in Blood and Hemolymph. Acta Zool. Taiwanica 1999, 10, 91–101. [Google Scholar]
  44. Zerbino, D.R.; Achuthan, P.; Akanni, W.; Amode, M.R.; Barrell, D.; Bhai, J.; Billis, K.; Cummins, C.; Gall, A.; Girón, C.G.; et al. Ensembl 2018. Nucleic Acids Res. 2018, 46, D754–D761. [Google Scholar] [CrossRef] [PubMed]
  45. Illumina, I. IGenomes. Available online: https://support.illumina.com/sequencing/sequencing_software/igenome.html (accessed on 21 October 2020).
  46. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast Universal RNA-Seq Aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef] [PubMed]
  47. Liao, Y.; Smyth, G.K.; Shi, W. FeatureCounts: An Efficient General Purpose Program for Assigning Sequence Reads to Genomic Features. Bioinformatics 2014, 30, 923–930. [Google Scholar] [CrossRef] [Green Version]
  48. R Studio Team. R Studio: Integrated Development Environment for R; R Studio, Inc.: Boston, MA, USA, 2020. [Google Scholar]
  49. R Core Team. R: A Language and Environment for Statistical Computing; R Core Team: Vienna, Austria, 2020. [Google Scholar]
  50. Love, M.I.; Huber, W.; Anders, S. Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [Green Version]
  51. Oakley, H.; Cole, S.L.; Logan, S.; Maus, E.; Shao, P.; Craft, J.; Guillozet-Bongaarts, A.; Ohno, M.; Disterhoft, J.; Van Eldik, L.; et al. Intraneuronal β-Amyloid Aggregates, Neurodegeneration, and Neuron Loss in Transgenic Mice with Five Familial Alzheimer’s Disease Mutations: Potential Factors in Amyloid Plaque Formation. J. Neurosci. 2006, 26, 10129–10140. [Google Scholar] [CrossRef]
  52. Todorov, H.; Fournier, D.; Gerber, S. Principal Components Analysis: Theory and Application to Gene Expression Data Analysis. Genom. Comput. Biol. 2018, 4, e100041. [Google Scholar] [CrossRef] [Green Version]
  53. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B (Methodol.) 1995, 57, 289–300. [Google Scholar] [CrossRef]
  54. Yu, G.; Wang, L.-G.; Han, Y.; He, Q.-Y. ClusterProfiler: An R Package for Comparing Biological Themes among Gene Clusters. Omics: A J. Integr. Biol. 2012, 16, 284–287. [Google Scholar] [CrossRef]
  55. Fisher, R.A. Statistical Methods for Research Workers; Springer: New York, NY, USA, 1992; pp. 66–70. [Google Scholar]
  56. Bagot, R.C.C.; Cates, H.M.M.; Purushothaman, I.; Lorsch, Z.S.S.; Walker, D.M.M.; Wang, J.; Huang, X.; Schlüter, O.M.M.; Maze, I.; Peña, C.J.J.; et al. Circuit-Wide Transcriptional Profiling Reveals Brain Region-Specific Gene Networks Regulating Depression Susceptibility. Neuron 2016, 90, 969–983. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Zhang, Y.; Chen, K.; Sloan, S.A.; Bennett, M.L.; Scholze, A.R.; O’Keeffe, S.; Phatnani, H.P.; Guarnieri, P.; Caneda, C.; Ruderisch, N.; et al. An RNA-Sequencing Transcriptome and Splicing Database of Glia, Neurons, and Vascular Cells of the Cerebral Cortex. J. Neurosci. 2014, 34, 11929–11947. [Google Scholar] [CrossRef] [PubMed]
  58. Wickham, H. Ggplot2: Elegant Graphics for Data Analysis; Springer: New York, NY, USA, 2016. [Google Scholar]
  59. Dusa, A. Venn: Draw Venn Diagrams. 2018. Available online: https://CRAN.R-project.org/package=venn (accessed on 9 February 2021).
  60. Brandscheid, C.; Schuck, F.; Reinhardt, S.; Schäfer, K.-H.; Pietrzik, C.U.; Grimm, M.; Hartmann, T.; Schwiertz, A.; Endres, K. Altered Gut Microbiome Composition and Tryptic Activity of the 5xFAD Alzheimer’s Mouse Model. J. Alzheimer’s Dis. 2017, 56, 775–788. [Google Scholar] [CrossRef] [PubMed]
  61. Landel, V.; Baranger, K.; Virard, I.; Loriod, B.; Khrestchatisky, M.; Rivera, S.; Benech, P.; Féron, F. Temporal Gene Profiling of the 5XFAD Transgenic Mouse Model Highlights the Importance of Microglial Activation in Alzheimer’s Disease. Mol. Neurodegener. 2014, 9, 33. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Proia, P.; di Liegro, C.M.; Schiera, G.; Fricano, A.; Di Liegro, I. Lactate as a Metabolite and a Regulator in the Central Nervous System. Int. J. Mol. Sci. 2016, 17, 1450. [Google Scholar] [CrossRef] [Green Version]
  63. Valvona, C.J.; Fillmore, H.L.; Nunn, P.B.; Pilkington, G.J. The Regulation and Function of Lactate Dehydrogenase A: Therapeutic Potential in Brain Tumor. Brain Pathol. 2016, 26, 3–17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Dienel, G.A.; Cruz, N.F. Aerobic Glycolysis during Brain Activation: Adrenergic Regulation and Influence of Norepinephrine on Astrocytic Metabolism. J. Neurochem. 2016, 138, 14–52. [Google Scholar] [CrossRef] [Green Version]
  65. Marti, H.H.; Jung, H.H.; Pfeilschifter, J.; Bauer, C. Hypoxia and Cobalt Stimulate Lactate Dehydrogenase (LDH) Activity in Vascular Smooth Muscle Cells. Pflügers Arch. Eur. J. Physiol. 1994, 429, 216–222. [Google Scholar] [CrossRef]
  66. Lev-Vachnish, Y.; Cadury, S.; Rotter-Maskowitz, A.; Feldman, N.; Roichman, A.; Illouz, T.; Varvak, A.; Nicola, R.; Madar, R.; Okun, E. L-Lactate Promotes Adult Hippocampal Neurogenesis. Front. Neurosci. 2019, 13, 403. [Google Scholar] [CrossRef] [PubMed]
  67. van Hall, G. Lactate Kinetics in Human Tissues at Rest and during Exercise. Acta Physiol. 2010, 199, 499–508. [Google Scholar] [CrossRef]
  68. Hering, G.O.; Hennig, E.M.; Riehle, H.J.; Stepan, J. A Lactate Kinetics Method for Assessing the Maximal Lactate Steady State Workload. Front. Physiol. 2018, 9, 310. [Google Scholar] [CrossRef] [Green Version]
  69. Larrabee, M.G. Lactate Metabolism and Its Effects on Glucose Metabolism in an Excised Neural Tissue. J. Neurochem. 1995, 64, 1734–1741. [Google Scholar] [CrossRef]
  70. Smith, D.; Pernet, A.; Hallett, W.A.; Bingham, E.; Marsden, P.K.; Amiel, S.A. Lactate: A Preferred Fuel for Human Brain Metabolism in Vivo. J. Cereb. Blood Flow Metab. 2003, 23, 658–664. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Dienel, G.A. The Metabolic Trinity, Glucose–Glycogen–Lactate, Links Astrocytes and Neurons in Brain Energetics, Signaling, Memory, and Gene Expression. Neurosci. Lett. 2017, 637, 18–25. [Google Scholar] [CrossRef]
  72. Suzuki, A.; Stern, S.A.; Bozdagi, O.; Huntley, G.W.; Walker, R.H.; Magistretti, P.J.; Alberini, C.M. Astrocyte-Neuron Lactate Transport Is Required for Long-Term Memory Formation. Cell 2011, 144, 810–823. [Google Scholar] [CrossRef] [Green Version]
  73. El Hayek, L.; Khalifeh, M.; Zibara, V.; Abi Assaad, R.; Emmanuel, N.; Karnib, N.; El-Ghandour, R.; Nasrallah, P.; Bilen, M.; Ibrahim, P.; et al. Lactate Mediates the Effects of Exercise on Learning and Memory through Sirt1-Dependent Activation of Hippocampal Brain-Derived Neurotrophic Factor (BDNF). J. Neurosci. 2019, 39, 2369–2382. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Cunnane, S.; Nugent, S.; Roy, M.; Courchesne-Loyer, A.; Croteau, E.; Tremblay, S.; Castellano, A.; Pifferi, F.; Bocti, C.; Paquet, N.; et al. Brain Fuel Metabolism, Aging, and Alzheimer’s Disease. Nutrition 2011, 27, 3–20. [Google Scholar] [CrossRef] [Green Version]
  75. Zuend, M.; Saab, A.S.; Wyss, M.T.; Ferrari, K.D.; Hösli, L.; Looser, Z.J.; Stobart, J.L.; Duran, J.; Guinovart, J.J.; Barros, L.F.; et al. Arousal-Induced Cortical Activity Triggers Lactate Release from Astrocytes. Nat. Metab. 2020, 2, 179–191. [Google Scholar] [CrossRef]
  76. Magistretti, P.J.; Allaman, I. Lactate in the Brain: From Metabolic End-Product to Signalling Molecule. Nat. Rev. Neurosci. 2018, 19, 235–249. [Google Scholar] [CrossRef]
  77. Díaz-García, C.M.; Mongeon, R.; Lahmann, C.; Koveal, D.; Zucker, H.; Yellen, G. Neuronal Stimulation Triggers Neuronal Glycolysis and Not Lactate Uptake. Cell Metab. 2017, 26, 361–374.e4. [Google Scholar] [CrossRef] [PubMed]
  78. Kuo, C.-C.; Chiang, A.W.; Lewis, N.E. Dysregulation of the Secretory Pathway Connects Alzheimer’s Disease Genetics to Aggregate Formation. bioRxiv 2020. [Google Scholar] [CrossRef]
  79. Hendrix, J.; Nijs, J.; Ickmans, K.; Godderis, L.; Ghosh, M.; Polli, A. The Interplay between Oxidative Stress, Exercise, and Pain in Health and Disease: Potential Role of Autonomic Regulation and Epigenetic Mechanisms. Antioxidants 2020, 9, 1166. [Google Scholar] [CrossRef]
  80. Kakanis, M.W.; Peake, J.; Brenu, E.W.; Simmonds, M.; Gray, B.; Hooper, S.L.; Marshall-Gradisnik, S.M. The Open Window of Susceptibility to Infection after Acute Exercise in Healthy Young Male Elite Athletes. Exerc. Immunol. Rev. 2010, 16, 119–137. [Google Scholar] [CrossRef] [Green Version]
  81. Guerriero, F.; Sgarlata, C.; Francis, M.; Maurizi, N.; Faragli, A.; Perna, S.; Rondanelli, M.; Rollone, M.; Ricevuti, G. Neuroinflammation, Immune System and Alzheimer Disease: Searching for the Missing Link. Aging Clin. Exp. Res 2017, 29, 821–831. [Google Scholar] [CrossRef]
  82. Boza-Serrano, A.; Yang, Y.; Paulus, A.; Deierborg, T. Innate Immune Alterations Are Elicited in Microglial Cells before Plaque Deposition in the Alzheimer’s Disease Mouse Model 5xFAD. Sci. Rep. 2018, 8, 1550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Passamonti, X.L.; Tsvetanov, X.K.; Jones, P.; Bevan-Jones, W.; Arnold, R.; Borchert, R.; Mak, E.; Su, L.; Rowe, X.J. Neurobiology of Disease Neuroinflammation and Functional Connectivity in Alzheimer’s Disease: Interactive Influences on Cognitive Performance. J. Neurosci. 2019, 39, 7218–7226. [Google Scholar] [CrossRef] [Green Version]
  84. Heppner, F.L.; Ransohoff, R.M.; Becher, B. Immune Attack: The Role of Inflammation in Alzheimer Disease. Nat. Rev. Neurosci. 2015, 16, 358–372. [Google Scholar] [CrossRef]
  85. Chi, L.M.; Wang, X.; Nan, G.X. In Silico Analyses for Molecular Genetic Mechanism and Candidate Genes in Patients with Alzheimer’s Disease. Acta Neurol. Belg. 2016, 116, 543–547. [Google Scholar] [CrossRef]
  86. Edgar, R.; Domrachev, M.; Lash, A.E. Gene Expression Omnibus: NCBI Gene Expression and Hybridization Array Data Repository. Nucleic Acids Res. 2002, 30, 207–210. [Google Scholar] [CrossRef] [Green Version]
Cells 10 00693 g001 550
Figure 1. Training groups and the comparison of first night and chronic saucer wheel usage. (A) Schematic of the three groups used for the investigation. All animals were single-caged and kept in type II macrolon cages for 30 days. The untrained group had no access to a saucer wheel, the chronically trained group kept the wheel over the whole time period, and the acutely trained mice had access only for the last night. Wheel turning counts per hour were measured for new access to the saucer wheel (B) and the chronic usage (C). Data are presented as mean of the performance over the 30 days of the chronically trained mice (n = 10 for wild-type, n = 8 for 5xFAD mice). Error bars are not visualized for clarity of the graph. Statistical analysis: multiple unpaired t-tests (* p < 0.05; *** p < 0.001).
Figure 1. Training groups and the comparison of first night and chronic saucer wheel usage. (A) Schematic of the three groups used for the investigation. All animals were single-caged and kept in type II macrolon cages for 30 days. The untrained group had no access to a saucer wheel, the chronically trained group kept the wheel over the whole time period, and the acutely trained mice had access only for the last night. Wheel turning counts per hour were measured for new access to the saucer wheel (B) and the chronic usage (C). Data are presented as mean of the performance over the 30 days of the chronically trained mice (n = 10 for wild-type, n = 8 for 5xFAD mice). Error bars are not visualized for clarity of the graph. Statistical analysis: multiple unpaired t-tests (* p < 0.05; *** p < 0.001).
Cells 10 00693 g001
Cells 10 00693 g002 550
Figure 2. Bodyweight development and food intake of 5xFAD and wild-type mice depending on physical training. (A) Body weight was measured at the same time of the day two times weekly. Here, the start weight and the weight at the end of the experiment are indicated. (B) Food consumption was measured weekly and is calculated per day. Data are presented as mean + standard error of the mean (SEM) (n = 6 for untrained and acutely trained, n = 7 for chronically trained mice). Statistical analysis: ordinary one-way ANOVA with the Fisher’s least significant difference (LSD) test (* p < 0.05; *** p < 0.001).
Figure 2. Bodyweight development and food intake of 5xFAD and wild-type mice depending on physical training. (A) Body weight was measured at the same time of the day two times weekly. Here, the start weight and the weight at the end of the experiment are indicated. (B) Food consumption was measured weekly and is calculated per day. Data are presented as mean + standard error of the mean (SEM) (n = 6 for untrained and acutely trained, n = 7 for chronically trained mice). Statistical analysis: ordinary one-way ANOVA with the Fisher’s least significant difference (LSD) test (* p < 0.05; *** p < 0.001).
Cells 10 00693 g002
Cells 10 00693 g003 550
Figure 3. Physical parameters in 5xFAD and wild-type mice after different activity regimens. Upon sacrifice, truncal blood glucose was measured (A) and Gastrocnemius of both hind legs (B), heart (C), and abdominal fat (D) were dissected and weighed. Data are presented as mean + SEM (n = 6 for untrained and acutely trained, n = 7 for chronically trained mice). Statistical analysis: ordinary one-way ANOVA with the Fisher LSD test (* p < 0.05; ** p < 0.01; *** p < 0.001).
Figure 3. Physical parameters in 5xFAD and wild-type mice after different activity regimens. Upon sacrifice, truncal blood glucose was measured (A) and Gastrocnemius of both hind legs (B), heart (C), and abdominal fat (D) were dissected and weighed. Data are presented as mean + SEM (n = 6 for untrained and acutely trained, n = 7 for chronically trained mice). Statistical analysis: ordinary one-way ANOVA with the Fisher LSD test (* p < 0.05; ** p < 0.01; *** p < 0.001).
Cells 10 00693 g003
Cells 10 00693 g004 550
Figure 4. Lactate levels in serum (A) and brain tissue (B) homogenates of 5xFAD and wild-type mice. Serum and brain homogenate lactate were assessed by an enzymatic assay. Values were normalized to protein content and are presented as mean + SEM (n = 6 for untrained and acutely trained, n = 7 for chronically trained mice). Statistical analysis: ordinary one-way ANOVA with the Fisher LSD test (*** p < 0.001).
Figure 4. Lactate levels in serum (A) and brain tissue (B) homogenates of 5xFAD and wild-type mice. Serum and brain homogenate lactate were assessed by an enzymatic assay. Values were normalized to protein content and are presented as mean + SEM (n = 6 for untrained and acutely trained, n = 7 for chronically trained mice). Statistical analysis: ordinary one-way ANOVA with the Fisher LSD test (*** p < 0.001).
Cells 10 00693 g004
Cells 10 00693 g005 550
Figure 5. Principal component analysis of normalized RNA counts from 5xFAD and wild-type mice after different training regimens. Using principal component analysis, the highest variances within the gene counts corresponding to principal component 1 (PC1) and principal component 2 (PC2) were calculated and plotted for (A) all (n = 24), (B) 5xFAD (n = 12), and (C) wild-type (n = 12) mice. In (A), dark green dots represent wild-type samples, while light green dots correspond to 5xFAD mice. The shapes of the symbols and the colors in (B,C) show the different training groups (red + circle = acute, orange + triangle = chronic, and yellow + squares = sedentary). The corresponding legend is located on the bottom.
Figure 5. Principal component analysis of normalized RNA counts from 5xFAD and wild-type mice after different training regimens. Using principal component analysis, the highest variances within the gene counts corresponding to principal component 1 (PC1) and principal component 2 (PC2) were calculated and plotted for (A) all (n = 24), (B) 5xFAD (n = 12), and (C) wild-type (n = 12) mice. In (A), dark green dots represent wild-type samples, while light green dots correspond to 5xFAD mice. The shapes of the symbols and the colors in (B,C) show the different training groups (red + circle = acute, orange + triangle = chronic, and yellow + squares = sedentary). The corresponding legend is located on the bottom.
Cells 10 00693 g005
Cells 10 00693 g006 550
Figure 6. Differential expression analysis results of different physical activity levels among wild-type mice. The volcano plots in (A–C) visualize the negative log10 of the adjusted p-value and the corresponding log2-fold change of the differential expression analyses between (A) chronic training and no training (B) acute training and no training (C) acute training and chronic training of wild-type mice. Each group consisted of four samples. Positive and negative log2-fold change values correspond to up- and downregulation in chronic (A) and acute training (B) compared to no training and acute training (C) compared to chronic training. Differentially expressed genes (DEGs) with an adjusted p-value < 0.05 are colored in red or blue corresponding to up- or downregulation, respectively. The five most significant genes are labeled by their gene symbols. All non-significant genes are colored gray. (D) shows the intersection of the DEGs found in each comparison. Empty fields represent a zero overlap. In (E), the number of genes for significantly enriched Gene Ontology (GO) terms of the acute-specific up- and downregulated genes from the differential expression analysis (DEA) of acute vs. sedentary and acute vs. chronic groups is visualized using bar plots. The color represents the significance of the adjusted p-value. (F) shows the cell type enrichment results for the acute-training-specific genes. The right and left bars correspond to upregulated and downregulated genes represented by the red and blue arrows’ direction. Green bars correspond to significantly enriched cell-types, while gray represents no significance.
Figure 6. Differential expression analysis results of different physical activity levels among wild-type mice. The volcano plots in (A–C) visualize the negative log10 of the adjusted p-value and the corresponding log2-fold change of the differential expression analyses between (A) chronic training and no training (B) acute training and no training (C) acute training and chronic training of wild-type mice. Each group consisted of four samples. Positive and negative log2-fold change values correspond to up- and downregulation in chronic (A) and acute training (B) compared to no training and acute training (C) compared to chronic training. Differentially expressed genes (DEGs) with an adjusted p-value < 0.05 are colored in red or blue corresponding to up- or downregulation, respectively. The five most significant genes are labeled by their gene symbols. All non-significant genes are colored gray. (D) shows the intersection of the DEGs found in each comparison. Empty fields represent a zero overlap. In (E), the number of genes for significantly enriched Gene Ontology (GO) terms of the acute-specific up- and downregulated genes from the differential expression analysis (DEA) of acute vs. sedentary and acute vs. chronic groups is visualized using bar plots. The color represents the significance of the adjusted p-value. (F) shows the cell type enrichment results for the acute-training-specific genes. The right and left bars correspond to upregulated and downregulated genes represented by the red and blue arrows’ direction. Green bars correspond to significantly enriched cell-types, while gray represents no significance.
Cells 10 00693 g006
Cells 10 00693 g007 550
Figure 7. Differential expression analysis results of 5xFAD compared to wild-type mice at different levels of physical activity. The volcano plots in A–C visualize the negative log10 of the adjusted p-value and the corresponding log2-fold change of the differential expression analyses between (A) sedentary 5xFAD vs. wild-type mice, (B) chronically trained 5xFAD vs. wild-type mice, and (C) acutely trained 5xFAD vs. wild-type mice. Each group consisted of four samples. Positive and negative log2-fold change values correspond to up- and downregulated genes in 5xFAD mice at the respective level of physical activity. Differentially expressed genes with an adjusted p-value < 0.05 are colored in red or blue, corresponding to up- or downregulation. The five most significant genes are labeled by their gene symbols. All non-significant genes are colored gray. (D) shows the intersection of the DEGs found in each comparison. Empty fields represent a zero overlap. (E) shows the cell type enrichment results for the sedentary-specific (left) and chronic-specific DEGs (right) while the right and left bars correspond to upregulated and downregulated genes represented by the direction of the red and blue arrows. Green bars visualize significantly enriched cell-types, while gray represents no significance.
Figure 7. Differential expression analysis results of 5xFAD compared to wild-type mice at different levels of physical activity. The volcano plots in A–C visualize the negative log10 of the adjusted p-value and the corresponding log2-fold change of the differential expression analyses between (A) sedentary 5xFAD vs. wild-type mice, (B) chronically trained 5xFAD vs. wild-type mice, and (C) acutely trained 5xFAD vs. wild-type mice. Each group consisted of four samples. Positive and negative log2-fold change values correspond to up- and downregulated genes in 5xFAD mice at the respective level of physical activity. Differentially expressed genes with an adjusted p-value < 0.05 are colored in red or blue, corresponding to up- or downregulation. The five most significant genes are labeled by their gene symbols. All non-significant genes are colored gray. (D) shows the intersection of the DEGs found in each comparison. Empty fields represent a zero overlap. (E) shows the cell type enrichment results for the sedentary-specific (left) and chronic-specific DEGs (right) while the right and left bars correspond to upregulated and downregulated genes represented by the direction of the red and blue arrows. Green bars visualize significantly enriched cell-types, while gray represents no significance.
Cells 10 00693 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wierczeiko, A.; Gammel, L.; Radyushkin, K.; Nguyen, V.T.T.; Todorov, H.; Gerber, S.; Endres, K. Voluntary Wheel Running Did Not Alter Gene Expression in 5xfad Mice, but in Wild-Type Animals Exclusively after One-Day of Physical Activity. Cells 2021, 10, 693. https://doi.org/10.3390/cells10030693

AMA Style

Wierczeiko A, Gammel L, Radyushkin K, Nguyen VTT, Todorov H, Gerber S, Endres K. Voluntary Wheel Running Did Not Alter Gene Expression in 5xfad Mice, but in Wild-Type Animals Exclusively after One-Day of Physical Activity. Cells. 2021; 10(3):693. https://doi.org/10.3390/cells10030693

Chicago/Turabian Style

Wierczeiko, Anna, Lena Gammel, Konstantin Radyushkin, Vu Thu Thuy Nguyen, Hristo Todorov, Susanne Gerber, and Kristina Endres. 2021. "Voluntary Wheel Running Did Not Alter Gene Expression in 5xfad Mice, but in Wild-Type Animals Exclusively after One-Day of Physical Activity" Cells 10, no. 3: 693. https://doi.org/10.3390/cells10030693

APA Style

Wierczeiko, A., Gammel, L., Radyushkin, K., Nguyen, V. T. T., Todorov, H., Gerber, S., & Endres, K. (2021). Voluntary Wheel Running Did Not Alter Gene Expression in 5xfad Mice, but in Wild-Type Animals Exclusively after One-Day of Physical Activity. Cells, 10(3), 693. https://doi.org/10.3390/cells10030693

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop