Literature DB >> 27672559

Behavioral and omics analyses study on potential involvement of dipeptide balenine through supplementation in diet of senescence-accelerated mouse prone 8.

Nobuhiro Wada1, Satoru Yamanaka2, Junko Shibato1, Randeep Rakwal3, Satoshi Hirako4, Yuzuru Iizuka5, Hyounju Kim5, Akiyo Matsumoto5, Ai Kimura4, Fumiko Takenoya6, Genta Yasunaga7, Seiji Shioda1.   

Abstract

This study investigates effects of dipeptide balenine, as a major component of whale meat extract (hereafter, WME), supplementation on senescence-accelerated mouse prone 8 (SAMP8), an Alzheimer's disease (AD) model at level of learning and memory formation and brain expression profiles genome-wide in brain. Mice fed experimental balenine (+ WME) supplemented diet for 26 weeks were subjected to four behavioral tests - open field, Y-maze, new object recognition, and water-filled multiple T-maze - to examine effects on learning and memory. Brain transcriptome of SAMP8 mice-fed the WME diet over control low-safflower oil (LSO) diet-fed mice was delineated on a 4 × 44 K mouse whole genome DNA microarray chip. Results revealed the WME diet not only induced improvements in the learning and memory formation but also positively modulated changes in the brain of the SAMP8 mouse; the gene inventories are publically available for analysis by the scientific community. Interestingly, the SAMP8 mouse model presented many genetic characteristics of AD, and numerous novel molecules (Slc2a5, Treh, Fbp1, Aldob, Ppp1r1a, DNase1, Agxt2l1, Cyp2e1, Acsm1, Acsm2, and Pah) were revealed over the SAMR1 (senescence-accelerated mouse resistant 1) mouse, to be oppositely regulated/recovered under the balenine (+ WME) supplemented diet regime by DNA microarray and bioinformatics analyses. Our present study demonstrates an experimental strategy to understand the effects of dipeptide balenine, prominetly contained in meat diet, on SAMP8, providing new insight into whole brain transcriptome changes genome-wide. The gene expression data has been deposited into the Gene Expression Omnibus (GEO): GSE76459. The data will be a valuable resource in examining the effects of natural products, and which could also serve as a human model for further functional analysis and investigation.

Entities:  

Keywords:  Alzheimer's disease; Balenine; Behavior; Mouse; SAMP8; Whole genome DNA microarray

Year:  2016        PMID: 27672559      PMCID: PMC5030327          DOI: 10.1016/j.gdata.2016.09.004

Source DB:  PubMed          Journal:  Genom Data        ISSN: 2213-5960


Introduction

The dipeptide, balenine (beta-alanyl-N tau-methyl histidine), has been identified long back in the extracts of muscles from 9 species of mammal, including man and chicken by comparative chromatography in combination with amino acid sequencing techniques [1]. As mentioned by Harris and Milne, 1987, “the dipeptide contained equimolar amounts of beta-alanine and N tau-methyl histidine with beta-alanine as the N-terminus [1]. Further, it was proposed that balenine, like carnosine and anserine is to be regarded as a normal constituent of muscle” [1]. To briefly give a background to the current study, our interest in looking at novel molecules for their potential beneficial effects on the human body was developed while investigating the nucleoprotein (NP) from salmon soft roe [2], especially how NP diet affects the brain and neurodegenerative disorders/diseases therein [3]. Among the major human diseases, and of great concern, are the neurodegenerative disorders, such as Parkinson's disease and Alzheimer's disease (AD), and which seriously limit the health and well-being of the world's aging population. Thus, the brain (and especially neurodegenerative disorders/diseases) becomes a hot topic for current researchers in this field, including our own group. With this backdrop, our present research examines the imidazole dipeptide, balenine – a major component of the meat, especially whale meat that is known to contain balenine in high concentration [4], [5], [6], in context of it being an important anti-fatigue food. This assumption is based on previous research, which shows the anti-fatigue effects with a ddY [7] swimming mouse model [8], and in humans under physical load and by daily activities [9]. These two studies [8], [9] suggested that whale meat through diet, in mice and as a supplement for human subjects, was most probably involved in “high efficiency for muscle endurance in mice” [8] and “attenuating high workload induced fatigue in humans” [9], indicating possible utilization of the whale meat as an anti-fatigue food, and maybe more. As balenine is the major constituent, and the objective of our research, we used the whale meat extract (WME) as a test sample to conduct the present research. The availability of the whale meat at Kyodo Senpaku Co. Ltd. (Tokyo, Japan) also helped to utilize it for our research as a candidate meat sample. It should be re-emphasized that we were primarily interested in the imidazole dipeptide, balenine, contained in the meat as an essential molecule/compound [4], [5], [6] behind these observed benefits for human health [8], [9], and not about whaling per se. To note, research into human wellness includes studies into the major life-style diseases and the benefits of not only searching for a suitable and healthy lifestyle but also and more importantly utilizing the potential of naturally occurring and protective compounds in our environment (plant or animal). In the present study, we report the results of our investigation into the effects of balenine through WME-supplemented diet in a mouse model of AD, the senescence-accelerated mouse prone 8 (SAMP8), using behavior and molecular genetic approaches to better understand the central nervous system dysfunction [10]. The reason for using SAMP8 is that it is a widely used model for AD and numerous studies have been conducted in relation to cognitive impairment (learning and memory deficits) and pathology, including oxidative stress and gene and protein expression analyses (for review see [11], [12]). Our present study builds up on two sets of experiments investigating the experimental balenine (+ WME)-supplemented diet over the low-safflower oil (hereafter, LSO) diet as a control diet, on SAMP8 and SAMR1 (senescence-accelerated mouse resistant 1) learning and memory deficits (experiment 1) and whole genome DNA microarray-based transcriptomics profiling in conjunction with Ingenuity Pathway Analysis (IPA; experiment 2). Results revealed that balenine (+ WME) supplementation in the diet of SAMP8 mouse resulted in an increase in the level of learning and memory formation and positive changes in the transcriptome of the brain in regards to cellular energy metabolism, suggested through the observation of recovery of gene expressions in the SAMP8 model over the not supplemented (control diet) mouse, and the potential involvement of the dipeptide balenine, as a major source of the beta-amino acid, beta-alanine, and other peptides to the cells.

Results

The experimental design and strategy for investigating effects of balenine (+ WME)-supplemented diet on SAMP8 mice

With an aim to investigate the potential beneficial effects of balenine (+ WME) we have here utilized a mouse model of SAMP8, and which also serves as an AD model. In this study, we linked behavioral analysis with omics approach via whole genome DNA microarray chip experiment and IPA analysis to reveal changes in the SAMP8 mice brain with or without WME supplemented diet and over the SAMR1 mice brain. As schematically presented in Fig. 1, the SAMP8 and SAMR1 mice were fed WME or control LSO diets for 26 weeks, followed by tests to determine improvements in learning and memory and thereafter extraction of total RNA for unraveling the brain transcriptome to first identify changes genome-wide and then attempt to delineate potential molecular factors behind the suggested improvements in learning and memory. Previous studies have also reported the WME-induced anti-fatigue effect in mouse and human [8], [9], and that was one reason for carrying out the present research. The other reason for why the meat component was used on the SAMP8 mice model, we provide a preliminary data on the radical (superoxide anion, O2 −) scavenging activity between different edible parts of the whale oil and meat (see Online Resource 1[13]). These data revealed that meat strongly removed O2– compared to other groups in vitro, suggesting its potential as an anti-oxidant or containing anti-oxidative properties, which was useful in context of the SAMP8 mice model. Next, we present below the results of each of the four behavioral tests and the transcriptome profiles of the SAMP8 mouse brains under balenine (+ WME) supplemented diet.
Fig. 1

The experimental strategy for investigating the effects of balenine (+ WME) supplementation diet in the senescence-accelerated mouse prone 8 (SAMP8) mouse. 24-weeks-old SAMP8 and senescence-accelerated mouse resistant 1 (SAMR1) mice were housed in individual cages till 50-weeks-old. Following the chow (CE-2) for 24 weeks, mice were fed low-safflower oil (LSO) diet as control diet or balenine through whale meat extract (WME) diet powders until sacrificed for 26 weeks. Two analyses were performed, four behavior tests for learning and memory followed by whole-genome DNA microarray experiment using whole brain in conjunction with bioinformatics analysis. Details are in the text.

Behavioral tests on the SAMP8 and SAMR1 mice fed balenine (+ WME)-supplemented diet

To investigate potential balenine (+ WME) diet-induced mechanisms underlying effects for learning and memory in SAMP8 as AD model mouse, we performed four tests. Firstly, mice were evaluated by open field test as measurement of locomotor activity. SAMR1 mice had significantly reduced locomotor activity (SAMR1 LSO; 852.8 ± 60.5). No differences between SAMP8 LSO and SAMP8 WME groups were observed (SAMP8 LSO; 1520.5 ± 93.9, SAMP8 WME; 1405.7 ± 76.1) (Fig. 2a). Secondly, the Y-maze test results showed that in the SAMP8 WME group degradation of spatial memory formation significantly improved compared to SAMP8 LSO group (SAMR1 LSO; 49.6 ± 4.2, SAMP8 LSO; 41.5 ± 4.0, SAMP8 WME; 56.6 ± 1.8%; p < 0.05) (Fig. 2b). Thirdly, for the NOR test, there were no differences among any of the groups (SAMR1 LSO; 65.0 ± 4.9, SAMP8 LSO; 44.9 ± 2.3, SAMP8 WME; 54.4 ± 4.5%) (Fig. 2c). Fourthly, the WMTM test revealed that the SAMP8 WME diet group had reduced escape time from the water in a time-dependent manner, and the observed escape time was significantly different at the 3rd day compared to SAMP8 LSO group (SAMP8 LSO; 135.9 ± 31.3, SAMP8 WME; 52.5 ± 5.8 2nd, p < 0.02), but the number of inserted blind array were no different among any of the groups (Fig. 2d).
Fig. 2

The four behavioral tests for examining the balenine (+ WME)-supplemented diet or not (LSO) effects in the SAMP8 mouse. We performed open field test (a), Y-maze test (b), novel object recognition test (NOR; c), and water-filled multiple T-maze test (WMTM; d). (a) Locomotor activity of SAMR1 mice were significantly decreased compared to SAMP8 LSO and WME group mice. However, no significant differences were observed between SAMP8 LSO group and balenine (+ WME) group. (b) In the evaluation of learning impairment using Y-maze test, the SAMP8 WME group showed a significant alteration in behavior compared to the SAMP8 LSO group. No differences were seen between SAMP8 (+ WME) and SAMR1 LSO groups (one-way comparison post-hoc test; *p < 0.05). (c) In the training for novel object recognition, first during the training period two objects (A and A′) were presented to a mouse. Second, during the testing phase, one of the conditioned objects was replaced by a differently colored novel object (A and B). During training all mice groups were found to equally touch the objects (A and A′). During testing, SAMR1 LSO and SAMP8 (+ WME) groups increased the percentage of touching novel objects, however there were no differences among groups. (d) In the escape latency from multiple T-maze in water maze, the SAMP8 (+ WME) group showed a significantly faster escape latency on day 3 compared to SAMP8 LSO group (repeated measures one-way ANOVA, Bonferroni's multiple comparison tests; #p < 0.02). LSO: low-safflower oil.

Global gene expression (transcriptomics) analysis of SAMP8 mice fed with balenine (+ WME)-supplemented diet over SAMP8 LSO (− WME) fed diet

DNA microarray analysis revealed 274 and 224 up (≧ 1.5-fold)- and down (≦ 0.75-fold)-regulated genes in the balenine (+ WME)-supplemented diet SAMP8 mouse whole brain, interestingly over the much higher number of genome-wide changes in gene expressions in the − WME (i.e. LSO diet only) SAMP8 mouse, where 904 and 613 up (≧ 1.5-fold)- and down (≦ 0.75-fold)-regulated genes were being altered (Fig. 3). It can also be said that the gene abnormalities (i.e. changed gene expressions) seen in the SAMP8 model mice brains are being strongly affected by the balenine (+ WME) supplementation. As to what is the nature of these gene expressions was the next question, we answer below. These genes are listed in Online Resource 2 (+ WME, up- and down-regulated gene lists) and Online Resource 3 (− WME (LSO), up- and down-regulated gene lists). These results provide not only the first inventory of WME supplemented diet-influenced genes (see Online Resource 2) in the animal (mouse) brain but also a comprehensive list of gene expressions in a mouse model of AD, i.e. SAMP8 (see Online Resource 3). To note, log transformed fold (log2 fold) change is used to represent the expression levels [14]. The reason for using ≧/≦ 1.5/0.75-fold as the cut-off value in the present experiment was to consider even smaller fold-change value genes in the experiment for downstream bioinformatics analyses. Moreover, as our experiment utilizes the two-color dye-swap approach (see below, Section 4.5) a highly confident (severe selection) differentially expressed gene list was generated in the analysis.
Fig. 3

Differentially expressed genes in the whole brain of SAMP8 mouse with or without balenine (WME) supplementation. The histograms show the up (red)- and down (blue)-regulated gene (≧/≦ 1.5/0.75-fold; log transformed fold (log2 fold) change) numbers in each sample condition. Balenine (+ WME) indicates SAMP8 mice fed with the experimental diet containing balenine as described in the Materials and methods section. − WME refers to the SAMP8 mice fed the control diet (LSO, low-safflower oil).

Comparative analysis of the genes oppositely regulated in the SAMP8 mice fed the balenine (+ WME)-supplemented diet and their role in cellular energy metabolism

The above DNA microarray analysis data also revealed 92 genes expressions which were recovered in the SAMP8 (+ WME) mouse brain over the same genes expression status in the SAMP8 LSO (− WME) mouse brain. The color codes (red, up-regulated; green, down-regulated) mark these genes in Table 1. Each gene recovery, either up-expressed or down-expressed, is indicated by their ratios. Looking at the genes expressions oppositely regulated by the diets (+ WME or LSO), and their functions in major metabolic pathways in the cell, some major trends were noticed in possible modulations of the energy metabolism via the tricarboxylic acid (TCA) cycle. Some of these genes (Slc2a5, Treh, Fbp1, Aldob, Ppp1r1a, DNase1, Agxt2l1, Cyp2e1, Acsm1 and Acsm2, and Pah) expressions are mentioned in a schematic pathway revolving around the said metabolic pathways depicting their up- or down-regulations status, in Fig. 4.
Table 1

Genes oppositely regulated in the SAMP8 mouse fed the balenine (+ WME)-supplemented diet. Color codes (red, up-expressed; green, down-expressed) mark the genes expressions recovering in the SAMP8 balenine (+ WME) mouse brain over the same genes expression state in the SAMP8 LSO (− WME) mouse brain. The genes differentially expressed in the SAMP8 mice brain (LSO diet) were determined using SAMR1 mice brain(LSO diet) gene expression data as a reference. LSO, Low-safflower oil.

Fig. 4

The genes of the major cellular metabolic pathways influenced in the whole brain of SAMP8 mouse with or without balenine (WME) supplementation. The red and blue bold arrows indicate up-regulated and down-regulated expression of genes (italicized) located intermediate to the respective substrates and products. Balenine (+ WME) indicates SAMP8 mice fed with the experimental diet containing balenine as described in the Materials and methods section. − WME refers to the SAMP8 mice fed the control diet (LSO, low-safflower oil).

Functional categorization and IPA analyses of differentially expressed brain genes of SAMP8 mice fed with balenine (+ WME)-supplemented diet and SAMP8 LSO (− WME) fed diet

The functionally categorized genes post-analysis using the SABiosciences PCR array list (Qiagen) are presented as color coded up-regulated (red) and down-regulated (green) genes numbers, for each category in both the SAMP8 (+ WME) and SAMP8 LSO (− WME) mouse brain (see Online Resource 4). At a glance the number of genes and gene categories highly expressing in the SAMP8 mouse brain under control (LSO) diet were reversed in the balenine (+ WME)-fed SAMP8 mouse brain (see also, Table 1). For a wider bioinformatics analysis, a biological function and network analysis was performed using the IPA tool generating numerous networks and pathways revealing the balenine (+ WME)-supplemented diet affecting the SAMP8 mouse brain in complex ways (see ESM 5, ESM 6). Each sample condition IPA analysis also provided us with a list of top molecules that are gene candidates for the potential mechanism/s underlying the balenine (+ WME)-supplemented diet effects in the SAMP8 model mouse (Fig. 5).
Fig. 5

The top molecules identified by IPA analysis in the whole brain of SAMP8 mouse with or without balenine (WME) supplementation. Gene names are in blue, and the up- and down-regulated genes along with the fold changes are indicated by red and green bold arrows for each sample condition. Balenine (+ WME) indicates SAMP8 mice fed with the experimental diet containing balenine as described in the Materials and methods section. − WME refers to the SAMP8 mice fed the control diet (LSO, low-safflower oil).

Discussion

Our study indicates that the balenine (+ WME)-supplemented diet has first, brought about a change in the behavior of the SAMP8 mice, namely, the improvement in the learning and memory formation, and secondly, caused numerous positive genome-wide changes in the brain of SAMP8 mouse over the LSO diet (control)-fed mouse. Inadvertently, these genomic data also revealed that the SAMP8 (AD model) mouse-fed LSO (− WME) had more numbers of changed gene expressions (1517 genes) than the SAMP8 mouse-fed (+ WME) diet (498 genes) (Fig. 3). This can be expected as the SAMP8 is an AD model mice [11], [12]. Interestingly many genes known to be affected in AD were also found to be present among these differentially expressed genes and when comparing the SAMP8 brains to the SAMR1 brains (seealso Table 1 and Fig. 5). We discuss some of these genes below in context of the recovery of gene expressions in the brain following the feeding of mouse with the balenine (+ WME)-experimental diet over the control diet (LSO, − WME). The present study thus provides genes and their potential roles/functions in context of the utilization of the balenine as a possible supplement for human wellness, especially in regards to the AD model used in this study. In the behavioral part of this study, SAMP8 balenine (+ WME)-diet group showed significantly differences compared to SAMP8 LSO diet (control) group in Y-maze and WMTM tests. Two previous reports have shown that WME-induced the anti-fatigue effect in mouse and human [8], [9], concluding that balenine as imidazole dipeptide [4], [5], [6] in the WME might affect the muscle endurance. Searching the literature we found that another imidazole dipeptide, carnosine, has been reported to have protective effects from ischemic vascular dementia [15], and that dietary supplementation of carnosine suppressed the accumulation amyloid-β in hippocampal region in 3 × Tg-AD mice as an AD model [16]. The chicken breast meat is known to abundantly contain carnosine [17]. That balenine is abundantly and specifically present in whale meat has also been reported [18], but there is no study to our knowledge on its (balenine) phygiological role. However, based on previous researches and our preliminary data (see Online Resource 1), we have considered the present study whereby, the results suggest an improvement in learning and memory function (Fig. 2) in SAMP8 balenine (+ WME)-fed supplemented diet mice. Our results also lead us to speculate that balenine may exert a similar function as the imidazole dipeptide carnosine [19]. This may include possible roles as activators of the enzyme adenosine triphosphate, formation of energy-rich phosphorus compounds, buffers (protective action on the mitochondrial structures), as have been discussed for anserine and carnosine, mostly in the muscle [19]. This however needs to thoroughly examined in the mouse brain and proved in the future studies. As will be discussed below, the balenine (+ WME)-supplemented diet is seen to cause, in the brain of the SAMP8 mouse, a recovery in the energy requirements or balance of the brain power via providing energy substrates and critical amino acids, including the non-essential amono acid beta-alanine (Fig. 4), that could be directly derived from balenine, in the WME degradation. Table 1 lists the genes (92) oppositely regulated/recovered by the balenine (+ WME)-supplemented diet. These genes present insight into the positive effects of the balenine-rich (WME) diet on the SAMP8 mouse, which we find to be impaired in metabolic functions as exemplified by the plotting of these genes on the energy metabolism pathway/s shown in Fig. 4. Though speculative at this stage, this DNA microarray-based gene expression data and analysis points twoards the fact that the balenine (WME) diet might exert a positive effect on the glucose and gluconeogenesis, ethanol fermenation, fatty acid degradation (a major route is the mitochondrial beta-oxidation pathway), aromatic amino acids, and finally the TCA cycle. In other words, it implies that balenine (+ WME)-supplemented diet might provide the necessary energy resources impaired in the SAMP8 mouse to provide the brain with necessary metabolites and amino acids for its proper functioning. Here, in relation to the above mentioned balenine (+ WME) positively affecting learning and memory functions in the SAMP8 mice and the glucose metabolism pathway genes being affected genome-wide, we would like to mention some previous literature pointing to a link between learning memory and glucose metabolism in aging brain and in causes of disorders causing dementia [20], [21], [22], [23], [24], [25]. For example, a study by Mayer and co-workers (1990) inducing irreversible hyperglycemia and depression of cerebral glucose metabolism via the intracerebroventricular (i.c.v.) or intraperitoneal (i.p.) administration of streptozotocin (STZ), which is known to damage insulin-producing cells, examined the behavioral functions (locomotor activity, learning and memory) in rat [24]. Results revealed that the STZ-i.c.v. injection rat was impaired in learning tasks compared with control groups and the STZ-i.p. injection group, suggesting impairment of passive avoidance learning might reflect a poorer acquisition and/or retention of memory in these rats [24]. This lead those authors to further suggest that poor cerebral regulation of glucose would be a primary metabolic event involved in impairments of learning and memory [24]. In 2003, Hoyer summarized studies specifically focused on the role of cerebral glucose metabolism and the function of insulin in the brain regarding memory capacity under normal and pathological conditions, including AD [25], [26]. An interesting result to mention from that paper is the finding that Ginkgo biloba extract (EGb 761®) exerts a beneficial effect on STZ-damaged rats returning them to completely normal values of cerebral energy metabolism, including partial compensation of the deficits in learning, memory and cognitive ability [25]. Indirectly, these data provide support to the beneficial effects of balenine (+ WME) diet on abnormalities in behavior and possible link to the recovery of the energy metabolism in the SAMP8 mouse model. In the following paragraphs, we discuss some of these genes presented in Table 1, first taking an example of the glucose metabolism/transport-related Slc2a5 gene. Slc2a5, which is a solute carrier family 2 (facilitated glucose transporter) member 5 (Glut5), is down-expressed in balenine(+ WME) over the LSO (− WME) diet. Although several of the GLUTs have been identified in the brain, it is GLUT1, which is the key transporter of glucose through the blood brain barrier [27]. The GLUT 5 (Glut5/Slc2a5), known as a fructose transporter having very low affinity for glucose, has been shown to be expressed in human and rat microglia cells but its functional significance remains unknown [27], [28]. Although glucose metabolism alterations have been suggested to play a role in AD, the mechanism remains unclear as also the link to AD [27]. Nevertheless, altered brain glucose levels and metabolism are suggested to be impaired in neurodegenerative diseases, and this might be the case with the SAMP8 model. Therein, we come across the second gene, namely Treh [29], which conversely to Slc2a5 is highly expressed in the balenine(+ WME) condition. This implies the TREHALASE, catalyzing the hydrolysis of trehalose into glucose providing an alternative source of energy [30], [31], is available in the balenine (+ WME) SAMP8 mouse and might be utilized for the production of glucose, based on the Treh gene expression level. It is also interesting to note that this implies an additional availability of trehalose, and which if might be proven in later studies, would be important factor in the improved learning and memory we find in the SAMP8 (+ WME) mice. The reason may lie in the fact that this disaccharide is also known as an antioxidant protecting cells from injury, suppressing protein aggregation including having cellular and behavioral beneficial effects in animal models of neurodegenerative disorders [32]. These results lead us to speculate on the availability of glucose in the SAMP8 mice fed the balenine (+ WME)-supplemented diet, which might help in recovering the said glucose deficiency in these mice. Similarly, we see the increased up-expression of the Fbp1 (fructose-1,6-bisphosphate 1) and Aldob (aldolase B, fructose-bisphosphate) genes in the (+ WME) SAMP8 mouse, again pointing towards activated glycolysis pathway, and leading to the TCA cycle. The other significant genes positively regulated by the balenine (+ WME) supplemented diet in SAMP8 mouse were, alanine-glyoxylate aminotransferase 2-like 1 (Agxt2l1; glutamate catabolism, promoting the conversion of l-alanine + glyoxylate to pyruvate + glycine, metabolism of l-glutamate), cytochrome P450, family 2, subfamily e, polypeptide 1 (Cyp2e1; conversion of ethanol to acetaldehyde), acyl-CoA synthetase medium-chain family member 1 and 2 (Acsm1, Acsm2; priming or activation of the fatty acids by forming fatty acyl-CoA thioesters), phenylalanine hydroxylase (Pah; tetrahydrobiopterin (BH4)-dependent enzyme phenylalanine 4-hydroxylase converts phenylalanine to tyrosine), and deoxyribonuclease 1 (DNase1; nucleates to ribose-5-P) (Fig. 4). Further, than their generally known functions, these genes are suggested here, in this study, to act also for the production of substrates to be utilized for the energy metabolism in the mitochondria as shown in Fig. 4. Other than the possibility of the above genes presented in Fig. 4, towards improved energy metabolism in the (+ WME) SAMP8 mouse brain, we would also like to discuss some unique characteristics of these gene products in relation to disease, specifically neurodegenerative disorders. The case of the novel gene up-expressed, Agxt2l1, a class II aminotransferase [33] is presented. First, the AGXT2L1 is a brain enzyme which might be involved in the pathogenesis of severe neuropsychiatric disturbances, and second, despite it being a transaminase was recently shown to have a stronger lyase activity (phospholipid metabolism) than an aminotransferase in human [34]. Taking the case of the human study, AGXT2L1 acts on O-phosphoethanolamine (PEA) irreversibly degrading it to acetaldehyde, phosphate and ammonia [34]. As described by Schiroli and co-workers (2013) [34], and we quote “PEA is both a byproduct of sphingolipid degradation and a precursor of glycerophospholipid biosynthesis and it can accumulate when membrane remodeling or degradation increase”; such an accumulation can lead to mitochondrial dysfunction [35]. Thus, AGXT2L1 might act as a protector of the toxic metabolite increase in cells or in the case of disease might serve as a compensatory mechanism to alleviate the effects of excessive membrane degradation [34]; is this the case in the SAMP8 mouse brain, though we do not know at present. The Agxt2l1 gene was also recently shown to be one of the three strongly up-regulated brain genes in the mouse model of psychiatric disorders/neurological diseases, such as compound C and lithium therapy, which have been used to study the effects on chromatin-modifying enzymes such as the histone deacetylase [36]. These data suggest an involvement of the Agxt2l1 gene in AD, though it should be noted that this gene function needs further investigation and validation as a novel target of lithium and neuroleptics, thus providing greater relevance to its decreases in depression and suicide and increases in bipolar disorder and schizophrenia [37]. Another exciting finding of an up-expressed gene is Pah, which, interestingly, causes a reduction in the amino acid phenylalanine turnover in chronic inflammation conditions that is known to affect dopamine, adrenaline and noradrenaline synthesis, as tyrosine is the precursor of these neurotransmitters [38]. Moreover, a cross-regulation between the immune response and severe metabolic changes is known and a connection therein with depression and cognitive impairment, which is in turn linked with changes in amino acid metabolism and neurotransmitter (serotonin, L-DOPA, dopamine, noradrenaline and adrenaline) synthesis [38], [39], [40], [41]. As recently reported by Gostner and co-workers (2015), the low-inflammatory stress-associated dysregulation of amino acid metabolism may be part of the pathophysiology of common HIV-associated neuropsychiatric conditions [38]. These reports led us to speculate that SAMP8 mice model of AD brains might have an amino acid imbalance and that is being possibly recovered in part by the up-expression of genes, such as Pah, with the balenine (+ WME)-supplemented diet. Herein, again the role of balenine (in the WME) may be called to question and further investigation. From the functional gene categories, the gene numbers in antifungal response, NFκB signaling targets, cancer inflammation and immunity crosstalk, inflammatory responses and autoimmunity related categories were highly dominant in the SAMP8 mouse (LSO diet) brain. Contrastingly, these functional categeories were considerably reduced in gene numbers in the SAMP8 mouse-fed balenine (+ WME) diet, brain (see Online Resource 3). Looking carefully at the literature it was found that brain infections have recently become an important topic in relation to AD [42], [43], [44]. A most recent publication showed that different brain regions (central nervous system) of AD patients are infected with fungal pathogens [45]. However, as those authors mention, it is not clear as to whether the fungal agents were the cause of AD or the patients, who were elderly, had a poor adaptive immune response resulting in them having fungal infections [45]. As also recently reported, researchers have considered the idea that AD might be an infectious disease or that infectious agents, include bacteria are risk factors for AD [46], [47]. The gene category for antibacterial response was also dominant in the SAMP8 (LSO diet) mouse brain (see Online Resource 3). Neuroinflammation or the inflammatory processes also are thought to play an important part in the pathogenesis of AD [48], [49], [50], and these gene categories are also well represented in the SAMP8 (LSO diet) mouse brain (see Online Resource 3). Combining the powerful approach of DNA microarray chip technology with the IPA bioinformatic tool, the possible brain gene interactions and pathways/networks could be delineated. These networks, one through six, for both SAMP8 mouse-fed balenine (WME) and LSO diets are shown in ESM 5, ESM 6. The top molecules presented in Fig. 5, of which some can also be seen in the Table 1 (oppositely regulated expressions in the balenine + WME mouse-fed diet), suggest their importance as key genes in the said positive effects of balenine (+ WME) supplementation. It would be interesting to analyze these genes (Kap, Cyp2j5, S100G, Keg1, UGT2B28, Igkv6-14, HLA-A, CCL19, etc.) individually and further bioinformatics analysis would shed light on their possible role/function and modulation by balenine (+ WME)-supplemented diet in the SAMP8 mouse. For example, the most highly up-regulated and top molecule (Fig. 5), namely, Kap gene that was also the most down-regulated gene in SAMP8 mouse-fed LSO diet showed a dramatic recovery in its expression with the balenine (+ WME) diet (Table 1). Kap encodes for the kidney androgen-regulated protein precursor [51], and its functions in the kidney are being widely investigated but nothing is known about it in the brain. Kap expression and presence raises new questions about its possible role in the brain of the SAMP8 mouse. In summary, this research is a first such study describing the effects of balenine (in WME) as a dietry supplement in a mice model (SAMP8) using both behavioral tests and high-throught DNA microarray-based omics approach, which provides new insight into the effects of the meat supplemented diet on the AD model mouse brain. To note, a previous study has also investigated the aging-related gene expression profiles in the frontal cortext of SAMP8 [52]. In particular, the energy metabolism pathways discussed above and amino acids provision to the brain point towards the beneficial positive effects of balenine (+ WME) diet on the SAMP8 mouse, which are interesting but require further experimental evidence. Our study helps set up a model system to investigate further the effects of specific whale meat compoments, such as the imidazole dipeptide balenine, and most probably other natural meats/products for possible beneficial effects on the human body. In addition, the obtained gene information resource in this study might serve well in investigating and designing new therapies that target these genes/proteins/enzymes of the said metabolic pathways for AD. The upcoming studies in our group will not only look directly at the effects of balenine in animal models of AD (SAMP8)/neurological diseases, but also utlize other omics-based approaches, proteomics [53] and metabolomics to try link the genes from this microarray analysis to their products for a better understanding of how the balenine (+ WME) supplemented diet affects the observed positive effects on brain recovery and function in the AD model mouse. We would argue at the end that this study was conducted only as a test to try get an insight into the positive effects of the dipeptide balenine (though dietry supplementation) a major component of the whale meat that formed/forms the diet of not only ancient, and to some extent current, Japanese people, but also in other whale eating traditions in the world.

Materials and methods

Animals and husbandary, supplemented diet, and sampling of whole brain

Six-weeks-old SAMP8 and SAMR1 mice (CLEA, Tokyo, Japan) were housed at the Animal Institution Facility in Showa University, and maintained in individual cages in a ventilated animal room with controlled temperature and relative humidity under a 12 h light: 12 h dark regime (8:00 AM, lights turned on). Mice were fed chow (CE-2, CLEA Japan) and tap water ad libitum until 24-weeks-old. Then, 24-weeks-old SAMs mice were given experimental diets, balenine (+ WME)-supplemented diet and LSO as a control diet, and both in a powder form. The WME was prepared from red meat of Antarctic minke whale Balaenoptera bonaerensis taken from the Japanese Whale Research Program under Special Permit in the Antarctic-Phase II in 2009/2010 by heat, enzyme and drying treatments, and as provided by Kyodo Senpaku Co. Ltd., Tokyo (Japan). The quality standard of whale meat was measured by Marugei Co. Ltd., (Hyogo, Japan), and results of the analysis are shown in Online Resource 7. The LSO diet components are shown in Online Resource 8. SAMP8 mice were randomly given balenine (in WME) or LSO diet, respectively. SAMR1 mice were given only the LSO diet, until 50-weeks-old. The behavioral tests, as detailed below, were performed at the timing of 49-weeks-old for eight days. At the end of the experiment (50-weeks of age) following the behavioral analysis and the last day of the feeding, the mice were removed their cages, decapitated and their brains carefully removed on ice. The whole brains were quickly frozen in liquid nitrogen in sterile freeze tubes, and stored at − 80 °C, till extraction of total RNA (Online Resource 9) for downstream DNA microarray analysis (see also 4.4, 4.5). All animal studies were conducted in accordance with the “Standards Relating to the Care Management of Experimental Animals” (Notice No. 6 of the Office of Prime Minister dated March 27, 1980) and with approval from the Animal Use Committee of Showa University (Approval Number: #04093).

Behavioral tests

To evaluate the effect of WME on the mouse, behavioral tests were carried out as follows: open field test, Y-maze test, new object recognition test (NOR) and water-filled multiple T-maze (WMTM) [54]. The number of animals in each test was as follows: SAMP8 LSO, n = 8; SAMP8 WME, n = 10; SAMR1 LSO, n = 10. Behavioral testing was performed between 10:00 AM and 6:00 PM. To habituate the animals to the testing room conditions, mice were transferred to the behavior testing room 1 h (9:00 AM to 10:00 AM) prior to beginning the first trial. The test room was not illuminated by the overhead lights but by using indirect lighting that pointed upwards to the roof. All apparatus were purchased from Bio-Research Co. Ltd. (Tokyo, Japan), and analysis of behavior test were analyzed by the SMART software 2.0 (Bio-Research Co. Ltd.). Each test was performed consequently for 8 days (open field test, 1 day; Y-maze test, 1 day; NOR, 2 days; WMTM, 4 days). We describe below, in detail, each test, and the images of each of the four apparatus set-up are schematically depicted in Online Resource 10.

Open field test

Mice were moved to the test room in their home cages and were handled by the base of their tails at all times. Each mouse was placed into the one of the four corners of the open field (45 cm × 45 cm × 40 cm) and allowed to explore the apparatus for 5 min (see Online Resource 10a). After 5 min, the mouse were returned to its home cage, and the open field was cleaned with oxidative water and permitted to dry between tests, followed by the second mouse.

Y-maze test

The Y-maze apparatus consisted of three arms made of gray plastic joined in the middle to form a “Y” shape (see Online Resource 10b). The walls of the arms were combinated, and the one arm component spec is 300 × 60 × 150 mm. The mouse was placed into one of the arms of the maze (start arm) and allowed to explore freely all three arms of the maze for 8 min. The number and the sequence of arms entered were recorded. The dependent variables were activity, defined as the number of arms entered, and percent alternation, calculated as the number of alternations (entries into three different arms consecutively) divided by the total possible alternations.

Novel object recognition test

The NOR test was performed using the open field box (see, open field test above) (see Online Resource 10c). Mice were individually habituated to the apparatus for 10 min at the 1st day (training day). The mouse can freely explore the open field area and touch the novel objects (A and A′) for 10 min. The animal behavior was recorded and we did a touch count, i.e., the number of times the mouse touched the objects was counted. On the next day, mice were individually let into the open field area again. We observed the mouse noticing change in the color of the object (A and B). If the mouse remembered the color of the object from the previous day, the tendency of the mouse moving towards the novel object is greater. This behavior alteration was recorded and the number of touches to the novel object was counted. Only 5 min were given to each mouse on the 2nd day.

Water-filled multiple T-maze test

The WMTM test is known as the multiple learning and memory test [54]. WMTM consists of seven T-maze combinations (20 cm wall height × 9 cm width), which are filled with water (22 ± 1 °C) to a height of 16 cm (see Online Resource 10d). During training (day 1), mice were habituated to the water by swimming three times a day. Mice remembered the goal position in order to escape from the water. During testing (days 2, 3 and 4), mice were placed at the start position in the water and allowed to swim freely for a maximum of 180 s. We measured when the mouse entered into the blind alley or not, and the time taken to arrive at the goal position. The interval time of this test is 2 h, and we performed this test three times in a day for three consecutive days.

Statistical analysis

All statistical analyses were performed with the Sigma Plot 12 (Hulinks, Tokyo, Japan) software package. Differences between three groups were assessed using the one-way analysis of variance (ANOVA) followed by Bonferroni's post-hoc test. Values were reported as the mean ± SE. Statistical significance was defined as p < 0.05.

Total RNA extraction

The deep-frozen whole mouse brain in a freeze tube was transferred from the deep freezer to a container with liquid nitrogen. The tissue was then carefully transferred to a pre-chilled (in liquid nitrogen) mortar and ground with a pestle to a very fine powder with liquid nitrogen [55], [56], [57]. Sample powders were immediately used for the next step of total RNA extraction by the QIAGEN RNeasy Mini Kit (QIAGEN, Maryland, USA) or stored at − 80 °C. Prior to DNA microarray analysis, the total RNA, quantity and quality measured spectrophotometrically with NanoDrop (Thermo Scientific, Wilmington, DE, USA) and re-confirmed using formaldehyde-agarose gel electrophoresis [55], [56], [57]. The obtained total RNA quality measurements revealed good quality as demonstrated by the A260/280 values > 1.8, A260/230 > 1.8; in our experiments, aiming for a higher value of > 2.0 to 2.3 [55], [56], [57]. Further, we obtained total RNA in optimum quantity of > 300 ng/μl. The total RNA was visualized for confirming pre-experimental integrity of the total RNA subunits by gel electrophoresis as shown in Online Resource 9.

DNA microarray chip-based whole genome brain transcriptome profiling

The total RNA that was extracted from the brain tissues for each animal was pooled in each group (SAMP8 LSO; SAMP8 WME; SAMR1 LSO) prior to DNA microarray analysis. As mentioned in Online Resource 9, for total RNA extraction individual biological replicates for SAMP8 LSO (n = 5), SAMP8 + WME (balenine) (n = 4), and SAMR1 LSO ()− WME) (n = 5) were used. These RNAs were pooled for application to the DNA microarray chip for one experiment. In addition, two technical replicates were performed by using two DNA microarray chips following the two-color dye-swap approach (detailed in Fig. 1). The reason for pooling the total RNA and not individual tissues is first, due to the usually larger variability at the tissue level than at the RNA level for estimating the gene expression, and second, in the use of lesser number of chips that was an economical/cost-effective option for the highly expensive microarray analysis [55], [58], [59]. All the steps followed in sequence provided confidence that we had the best quality total RNA for use in the DNA microarray chip. As a last step, the Agilent mouse whole genome 4 × 44 K (G4122F) DNA slide (composed of four chips on one slide) was used for the microarray analysis, which was performed according to the Agilent instructions (Agilent Technologies, Santa Clara, CA, USA) and detailed in our publications [55], [56], [57], [60]. Our established DNA microarray analysis experiment is also based on the two-color dye-swap approach providing a more-stringent selection condition for changed gene expression profiling than a simple single/two-color approach [55]. Briefly, total RNA (200 ng) was labeled with either Cy3 or Cy5 dye using an Agilent Low-RNA-Input Fluorescent Linear Amplification Kit (Agilent Technologies). Fluorescently labeled targets of control (SAMP8 −/SAMR1 −) as well as treatment (SAMP8 +/SAMP8 −) indicated in Fig. 1 were hybridized to the same microarray slide with 60-mer probes. A flip-labeling (dye-swap or reverse labeling with Cy3 and Cy5 dyes) procedure was followed to nullify dye bias associated with unequal incorporation of the two Cy dyes into cDNA. Differentially expressed genes were considered genes that were up-regulated in chip-1 but down-regulated in chip-2 for each group. Hybridization and wash processes were performed according to the manufacturer's instructions, and the hybridized microarrays were scanned using an Agilent Microarray scanner G2505C. To detect significantly differentially expressed genes between the control and treated samples, each slide image was processed using Agilent Feature Extraction software (version 11.0.1.1). This program measures the Cy3 and Cy5 signal intensities of whole probes. Dye bias tends to be signal-intensity-dependent; therefore, the software selected probes using a set-by-rank consistency filter for dye normalization. Said normalization was performed by LOWESS (locally weighted linear regression), which calculates the log ratio of dye-normalized Cy3- and Cy5 signals and the final error of the log ratio. Statistical significance of the log ratio values obtained for individual microarray probes was calculated using the “most conservative error model” default function of the Feature Extraction software, by which both a propagated and a universal error model are evaluated, and the higher (more conservative estimate of error) p-value of two error models was reported (Agilent Technologies [55]). For the analyses, threshold of significant differentially expressed genes is < 0.01 (for the confidence that the feature was not differentially expressed). Additionally, erroneous data generated due to artifacts were eliminated before data analysis by the software.

Functional categorization of differentially expressed genes

Pathway- and disease state-focused gene classification of the differentially expressed genes in the brain were classified based on the available categories of > 100 biological pathways or specific disease states in the SABiosciences PCR array list (Qiagen; www.sabiosciences.com) for Mus musculus.

IPA analysis

The biological function and network analyses were performed using IPA (Ingenuity® Systems, www.ingenuity.com). The dataset from the microarray [balenine diet (WME) supplementation in SAMP8 (+) or in SAMP8 (−)], which includes the differentially expressed (≧/≦ 1.5/0.75-fold compared with control) genes and their corresponding fold change values, was uploaded as an Excel spread sheet into the IPA tool. To create gene networks, the genes were overlaid onto a global molecular network that was developed from information that was contained in the ingenuity knowledge base. The functional analysis identified the biological functions that were most significant to the dataset (p-value < 0.05) according to a right-tailed Fisher's exact test.

Access to gene array data

The outputs of microarray analysis used in this study are available under the series number GSE 76459 (whole genome - http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE76459; at the NCBI Gene Expression Omnibus (GEO) public functional genomics data repository (http://www.ncbi.nlm.nih.gov/geo/info/linking.html). The following are the supplementary data related to this article.

ESM 1

Radical scavenging activity between different edible parts of the whale (oil and meat). For the determination of superoxide anion (O2 −), produced by HX-XOD reaction (total volume: 200 μl) 2 mM HX in 0.1 M phosphate buffer (pH 7.4), 50 ml PB, 0.5 mM DETAPAC 20 ml, 8% DMPO 30 ml, sample (in DMPO) 40 ml, H2O 30 ml, and XOD (0.5 U/ml in PB) 30 ml [13]. Each sample was measured with time constant and scanning time after 1min, respectively. These data showed that Mink whale meat strongly removed O2 − compared to other groups in vitro.

ESM 2

SAMP8 mouse fed balenine (+ WME) supplemented diet influenced up- and down-regulated genes in the whole brain. The low-safflower (LSO) diet fed mice were the control.

ESM 3

SAMP8 mouse brain up- and down-regulated genes. The mice were fed the control, low-safflower (LSO) diet.

ESM 4

Pathway- and disease states-focused gene classifications of genes in the whole brain of SAMP8 mouse with or without balenine (+ WME) supplementation. For each category, the number of genes are indicated in both the SAMP8 (+ WME) and SAMP8 LSO (− WME) mouse brain, where the red and green color codes indicate the up-regulated and down-regulated genes. All up- and down-regulated genes (≧/≦ 1.5/0.75-fold) for each sample microarray data (see Fig. 3) were used for the analysis using an Excel program developed by our group at Showa University. + WME indicates SAMP8 mice fed with the experimental diet containing the whale meat extract powder as described in the Materials and methods section. Low-safflower oil (LSO) refers to the SAMP8 mice fed the control diet.

ESM 5

The top six networks (1, 2, 3. 4, 5 and 6) by IPA analysis for the whole brain of SAMP8 mouse fed the balenine (+ WME) supplemented diet. The genes (up-regulated – red; down-regulated – green) are networked based on the evidence in the IPA. For further details, see the Materials and methods section.

ESM 6

The top six networks (1, 2, 3. 4, 5 and 6) by IPA analysis for the whole brain of SAMP8 mouse fed the control low-safflower oil (LSO) diet. The genes (up-regulated – red; down-regulated – green) are networked based on the evidence in the IPA. For further details, see the Materials and methods section.

ESM 7

Results of quality standard analysis of the whale meat used in this study.

ESM 8

Component of the experimental diet, with or without whale meat extract.

ESM 9

Total RNA extraction. To verify the quality of this RNA, the yield and purity were determined spectrophotometrically with NanoDrop (Thermo Scientific, Wilmington, DE, USA) and confirmed using formaldehyde-agarose gel electrophoresis. SAMP8 LSO, n = 5; SAMP8 WME (balenine), n = 4; SAMR1 LSO, n = 5; the remaining mice brains were stored for conducting immunohistochemistry analysis. LSO, low-safflower oil.

ESM 10

The apparatus for the behavior tests. (a) Open field box, 45 cm (w) × 45cm (d) × 40 cm (h) square box. (b) Y-maze, 30 cm (w) × 60 cm (d) × 15 cm (h); Y-shape construction. (c) NOR, styrol cylinder of 5 cm diameter set on the open field box. (d) WMTM, 82 cm (w) × 82 cm (d) × 20 cm (h) multiple T-maze apparatus; straight was 20 cm, maze pathway was 9 cm wide, and water was filled to a height of 16 cm at 22 ± 1 °C.

Conflict of interest

The authors declare no conflict of interests.
  52 in total

1.  ISOLATION OF 3-METHYL HISTIDINE FROM WHALEMEAT EXTRACT AND THE PREPARATION OF SOME DERIVATIVES.

Authors:  D H COCKS; P O DENNIS; T H NELSON
Journal:  Nature       Date:  1964-04-11       Impact factor: 49.962

Review 2.  Immune attack: the role of inflammation in Alzheimer disease.

Authors:  Frank L Heppner; Richard M Ransohoff; Burkhard Becher
Journal:  Nat Rev Neurosci       Date:  2015-06       Impact factor: 34.870

3.  A new diabetic mouse model derived from the ddY strain.

Authors:  Ichiro Noge; Yoshiyuki Kagawa; Toshio Maeda
Journal:  Biol Pharm Bull       Date:  2010       Impact factor: 2.233

Review 4.  The molecular evolution of pyridoxal-5'-phosphate-dependent enzymes.

Authors:  P K Mehta; P Christen
Journal:  Adv Enzymol Relat Areas Mol Biol       Date:  2000

5.  Computed tomography and positron emission transaxial tomography evaluations of normal aging and Alzheimer's disease.

Authors:  M J de Leon; S H Ferris; A E George; B Reisberg; D R Christman; I I Kricheff; A P Wolf
Journal:  J Cereb Blood Flow Metab       Date:  1983-09       Impact factor: 6.200

Review 6.  The senescence-accelerated prone mouse (SAMP8): a model of age-related cognitive decline with relevance to alterations of the gene expression and protein abnormalities in Alzheimer's disease.

Authors:  D Allan Butterfield; H Fai Poon
Journal:  Exp Gerontol       Date:  2005-07-18       Impact factor: 4.032

Review 7.  Neuroimmunomodulation in the pathogenesis of Alzheimer's disease.

Authors:  Inelia Morales; Gonzalo Farías; Ricardo B Maccioni
Journal:  Neuroimmunomodulation       Date:  2010-02-04       Impact factor: 2.492

8.  Ethanolamine and phosphoethanolamine inhibit mitochondrial function in vitro: implications for mitochondrial dysfunction hypothesis in depression and bipolar disorder.

Authors:  Josephine S Modica-Napolitano; Perry F Renshaw
Journal:  Biol Psychiatry       Date:  2004-02-01       Impact factor: 13.382

Review 9.  Senescence-accelerated mouse (SAM) with special references to neurodegeneration models, SAMP8 and SAMP10 mice.

Authors:  Toshio Takeda
Journal:  Neurochem Res       Date:  2009-02-07       Impact factor: 3.996

Review 10.  Disturbed Amino Acid Metabolism in HIV: Association with Neuropsychiatric Symptoms.

Authors:  Johanna M Gostner; Kathrin Becker; Katharina Kurz; Dietmar Fuchs
Journal:  Front Psychiatry       Date:  2015-07-14       Impact factor: 4.157

View more
  2 in total

1.  The Body-wide Transcriptome Landscape of Disease Models.

Authors:  Satoshi Kozawa; Ryosuke Ueda; Kyoji Urayama; Fumihiko Sagawa; Satsuki Endo; Kazuhiro Shiizaki; Hiroshi Kurosu; Glicia Maria de Almeida; Sharif M Hasan; Kiyokazu Nakazato; Shinji Ozaki; Yoshinori Yamashita; Makoto Kuro-O; Thomas N Sato
Journal:  iScience       Date:  2018-03-29

2.  Detection of Balenine in Mouse Plasma after Administration of Opah-Derived Balenine by HPLC with PITC Pre-Column Derivatization.

Authors:  Yasutaka Shigemura; Yu Iwasaki; Yoshio Sato; Tomomi Kato; Takuya Seko; Kenji Ishihara
Journal:  Foods       Date:  2022-02-18
  2 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.