Janina Feige-Diller1,2*†
Marisol Herrera-Rivero3*†
Anika Witten3,4
Monika Stoll2,3
Sylvia Kaiser1
S. Helene Richter1,2
Norbert Sachser1,2- 1Department of Behavioral Biology, University of Münster, Münster, Germany
- 2DFG RTG EvoPAD, WWU Münster, Münster, Germany
- 3Department of Genetic Epidemiology, Institute of Human Genetics, University of Münster, Münster, Germany
- 4Core Facility Genomics, Medical Faculty, University of Münster, Münster, Germany
Background: During early phases of life, such as prenatal or early postnatal development and adolescence, an organism's phenotype can be shaped by the environmental conditions it experiences. According to the Match-Mismatch hypothesis (MMH), changes to this environment during later life stages can result in a mismatch between the individual's adaptations and the prevailing environmental conditions. Thus, negative consequences in welfare and health can occur. We aimed to test the MMH in the context of food availability, assuming adolescence as a sensitive period of adaptation.
Methods: We have previously reported a study of the physiological and behavioral effects of match and mismatch conditions of high (ad libitum) and low (90% of ad libitum intake) food availability from adolescence to early adulthood in female C57BL/6J mice (n = 62). Here, we performed RNA-sequencing of the livers of a subset of these animals (n = 16) to test the effects of match and mismatch feeding conditions on the liver transcriptome.
Results: In general, we found no effect of the match-mismatch situations. Contrarily, the amount of food available during early adulthood (low vs. high) drove the differences we observed in final body weight and gene expression in the liver, regardless of the amount of food available to the animals during adolescence. Many of the differentially expressed genes and the corresponding biological processes found to be overrepresented overlapped, implicating common changes in various domains. These included metabolism, homeostasis, cellular responses to diverse stimuli, transport of bile acids and other molecules, cell differentiation, major urinary proteins, and immunity and inflammation.
Conclusions: Our previous and present observations found no support for the MMH in the context of low vs high food availability from adolescence to early adulthood in female C57BL/6J mice. However, even small differences of approximately 10% in food availability during early adulthood resulted in physiological and molecular changes with potential beneficial implications for metabolic diseases.
Background
Phenotypic plasticity allows the shaping of an individual's phenotype by the environment it experiences during early phases of life (1). Therefore, the developmental course of an individual might be adjusted according to predictions of the future environmental conditions based on information gained throughout these early phases (1, 2). As stated by the Match-Mismatch hypothesis (MMH), an individual thus effectively matched to the prevailing and future environmental conditions will be less prone to illnesses later in life. In contrast, if the environment drastically changes during the mature phase, a mismatch of the individual's phenotype and the environmental conditions can occur, possibly causing an increased disease susceptibility (3–7). The consequences of such a mismatch might be especially severe if fundamental needs, such as nutrition, are affected. Evidence in humans indicates a mismatch effect caused by discrepancies in food availability between prenatal or early postnatal life and adulthood. Data from historical events, e.g., the Dutch potato famine, regarding a mismatching situation of low food availability during early life and high food availability during later life indicated an increased risk of various cardiovascular and metabolic diseases (8–12). Similarly, individuals shaped by high food availability during early life may have an increased risk to suffer from malnutrition and comorbidities in situations of scarce nutrition, such as prison camps or famines (1, 13). Contrarily, individuals who experienced a matched low food availability were found to have smaller stature and a metabolism that favors the laying down of fat (3, 6), causing a reduced risk to suffer when exposed to poor nutrition (14). Furthermore, in the opposite situation of matched high food availability, a decreased risk of obesity and comorbidities during adulthood has been reported (15–17).
So far, match-mismatch effects have been mainly reported for discrepancies between the prenatal or neonatal phase and later life (9, 14–16, 18, 19). However, adolescence also represents a sensitive developmental period (20, 21) as it is, for example, associated with a modulation of the hypothalamic-pituitary-adrenal (HPA) axis (22), which plays an important role in energy regulation (23). Therefore, we hypothesized that adolescence represents a sensitive period for the long-term shaping of an individual's metabolism. To investigate this, our group applied a match-mismatch design on female mice exposed to low or high food availability in adolescence and early adulthood. Recently, we reported that the match-mismatch situation showed no significant effects on the welfare and health of the animals, as assessed by physiological and behavioral measurements. However, various effects of low vs. high food availability were observed, including lower relative liver weights and glucocorticoid secretion in a situation of plentiful nutrition. Taken together, the effects suggested a short-term adjustment of the metabolism and immune system to the prevailing situation of food availability (24). Here, we used RNA-sequencing to investigate changes in the liver transcriptome of a subset of these animals exposed to matching or mismatching food availability in adolescence and early adulthood. The liver is fundamentally important for lipid and glucose metabolisms [reviewed in Rui (25)], but also strongly involved in immune function (26–28). Consistent with our previous findings, we found no match-mismatch effects. Transcriptional changes in the liver were rather induced by low vs. high food availability during early adulthood, and implicated a wide range of biological processes related to metabolism as well as, to some extent, to immunity and inflammation.
Methods
Experimental Animals and Collection of Tissue Samples
All mice used for the present study were part of our previously published work, and detailed experimental protocols have been described (24). Briefly, female C57BL/6J mice were assigned to matching or mismatching conditions of food availability during adolescence and early adulthood to test the Match-Mismatch hypothesis (Figure 1A). High food availability was simulated by ad libitum access to a standard laboratory food (Altromin 1324, Altromin GmbH, Lage, Germany). In contrast, low food availability was characterized by feeding the animals once per day with a mildly restricted diet (10% reduction of the ad libitum intake). In phase 1, two groups of adolescent animals (postnatal day-PND- 28–70 ± 1) were exposed to low (L) food availability, while another two groups were exposed to high (H) food availability for a period of 6 weeks. During this phase, animals were not considered fully-grown. L and H cages were paired according to matching body weights. The amount of food in the L cage was then kept to approximately 90% of its respective H cage. In this manner, the continued growth of animals in the L cages, despite of food restriction, was warranted (29, 30). In phase two, when animals were considered to have reached early adulthood and full growth (PND 71 ± 2), they were kept either in the same (match) or the opposite (mismatch) situation of food availability for another 6 weeks. During this phase, animals exposed to low food availability were either kept to 90% of their own maximum body weight (group HL) or to that of the mice in the respective paired cage (group LL) during the last week of the previous phase. The groups studied thus consisted of the match or mismatch conditions HH (n = 15), HL (n = 15), LH (n = 16) and LL (n = 16). Subsequently, all animals underwent a series of behavioral tests over a period of 5 weeks (PND 111 – 142 ± 1), while remaining in the same situation as after the match/mismatch. At least 5 days after the last behavioral test, 16 randomly selected animals were anesthetized using 2.5% isoflurane in oxygen and decapitated. After dissection and removal of excessive fat tissue, the weights of the hearts, kidneys, livers, spleens and adrenal glands were measured. Livers were snap-frozen in liquid nitrogen and stored at −70°C.
Figure 1. (A) Experimental design. Mice were kept at high (H) or low (L) food availability conditions during 6 weeks of adolescence. Upon reaching early adulthood, food availability either remained the same (match) or was changed (mismatch) for another 11 weeks. During the last 5 weeks, behavioral tests were conducted for a previous publication (24). Subsequently, final body weights and body length were assessed; liver samples were weighed and collected for RNA-sequencing. (B) Progressive increases in body weight were assessed at weeks 1, 2, 6, 7, 8, 12 and at the end of the experiment. Afterwards, liver weight was recorded. Significant differences (Tukey HSD adj. p < 0.05) between groups are denoted by (*). L, low food availability; H, high food availability.
All experimental procedures complied with the regulations covering animal experimentation within Germany (Animal Welfare Act) and the EU (European Communities Council DIRECTIVE 2010/63/EU). The study was approved by the corresponding local (Gesundheits-und Veterinäramt Münster, Nordrhein-Westfalen) and federal authorities (Landesamt für Natur, Umwelt und Verbraucherschutz Nordrhein-Westfalen ‘LANUV NRW,” reference number 84–02.04.2018.A067).
Statistical Testing of Descriptive Phenotypes
Differences between the match and mismatch groups in terms of body weight at weeks 1, 2, 6, 7, 8, and 12, final body weight, final body length, weight of the liver and relative weight of the liver were tested through analysis of variance (ANOVA), followed by the Tukey's honest significant difference (HSD) method applied with a family-wise confidence level of 95%. Significance was set at p < 0.05.
RNA Extraction, Library Preparation and Sequencing
Liver tissue samples from four animals in each match or mismatch group (HH, LL, HL and LH) were destined for analysis of their transcriptomes via RNA-sequencing (RNA-seq). Total RNA was isolated using the Direct-zol RNA Microprep Kit (Zymo Research), followed by a DNase digestion step. Library preparation was carried out upon mRNA enrichment with the NEBNext Poly(A) mRNA Magnetic Isolation Module, using the NEBNext Ultra II RNA Directional Library Prep Kit for Illumina (New England BioLabs). Single read sequencing took place on a NextSeq 500 System (Illumina), using the corresponding NextSeq 500 High Output Kit v2.5, with a read length of 75 base pairs. The integrity of the RNA and quality of the library were assessed using a TapeStation 4200 (Agilent).
Data Pre-processing
Using a molecular barcode, the data was automatically demultiplexed using the Illumina bcl2fastq2 Conversion Software v2.20. FastQ files underwent two rounds of quality control, pre-trimming and post-alignment, using FastQC v0.11.7 (31). Removal of Illumina adapters and low-quality sequences was performed with Trimmomatic v0.38 (32). Reads of length<15 bases, as well as leading and/or trailing bases with quality<3 or no base call, and bases with average quality<15 in a 4-base sliding window were removed. Alignment was performed with HISAT2 v2.1.0 (33) using the mouse genome assembly mm10 (Mus musculus, GRCm38). Mapped reads (primary alignments) were sorted by read name using SAMtools v1.8 (34), and read counts were calculated with HTSeq v0.11.2 (35).
Differential Expression and Functional Enrichment Analyses
Differential expression between all possible combinations of match or mismatch groups (HH vs. HL, HH vs. LH, HL vs. LH, LL vs. HH, LL vs. HL and LL vs. LH) was assessed using DESeq2 (36). Raw read counts were filtered to remove genes with <10 counts prior to analysis. DESeq2 tests for the statistical significance of coefficients in a negative binomial generalized linear model using the Wald test, and corrects the p-values for multiple comparisons according to the Benjamini-Hochberg method. These models were adjusted for experimental batch of the animals. Genes were considered differentially expressed when adjusted-p (p.adj) <0.05 (a log2 fold change (lfc) threshold was not applied). Furthermore, to provide a biological context to these findings, each list of differentially expressed genes (DEGs) was subjected to functional enrichment analysis for Gene Ontology biological processes (GO_BP) in Mus musculus using BiNGO v3.0.3 (37) for Cytoscape v3.7.1 (38). GO_BP terms were considered significantly enriched following a hypergeometric test for overrepresentation, corrected for multiple comparisons according to the Benjamini-Hochberg method (p.adj<0.05).
Results
A basic description of the subsample used for transcriptomic profiling of the livers of mice placed in match or mismatch combinations of low and high food availability from our previous report (24) is shown in Table 1. In this subsample of 16 mice, we observed differences in the animals' body weights from week 8, when mice in the LH and HH groups started to show significant weight gains (week 8 p = 0.0165, week 12 p = 0.0077, final measurement p = 0.0017) in comparison to animals subjected to the LL and HL conditions. Liver weights differed only between mismatch groups, as the higher and lower values were found for these mice. Livers of animals in the LH group were significantly heavier than those in the HL group (Figure 1B). From RNA-seq, alignment rates > 97% and over 25 million reads were obtained for all samples. Visualization of the first two principal components (PCs) of the normalized read counts showed clustering of the liver transcriptomes according to food availability during phase two, which could explain up to 46% of the phenotypic variance between all samples (Figure 2A).
Table 1. Basic description of the sample used for transcriptomic profiling by RNA-sequencing.
Figure 2. (A) Plot of the first 2 principal components (PCs) of the normalized read counts. (B) Overlaps of differentially expressed genes (DEGs) and (C) the corresponding enriched Gene Ontology biological processes (GO_BPs) among all comparisons. No DEGs and/or GO_BP terms overlapping among all comparisons performed were identified. L, low food availability; H, high food availability.
The results obtained from the series of differential expression (Figure 3, Supplementary Tables 1–6) and corresponding functional enrichment analyses performed (Supplementary Table 7) are summarized in Table 2. In general, judging by the numbers of DEGs found for each comparison and, similar to what was observed for the final body weight, the liver transcriptomes changed significantly between match or mismatch groups with opposing feeding conditions during early adulthood. This suggested that the observed molecular changes were induced by the later (early adulthood) rather than by the earlier (adolescence) diet scheme, and that there was no significant effect of the match-mismatch situation in our study. Moreover, many overlaps were observed among DEGs (Figure 2B) and GO_BP terms (Figure 2C) between comparisons. When food availability was the same during early adulthood, changes in the liver transcriptomes between animals were limited. Nevertheless, a higher proportion of the GO_BP terms enriched in the latter case were specific for these comparisons. Table 3 provides a summary of the comparison-specific GO_BP terms found for each group comparison performed.
Figure 3. Overview of differential expression results for all comparisons. MA plots show the statistically significant (adjusted p < 0.05) dysregulations obtained in each comparison marked in red.
Table 2. Summary of results from the differential expression and functional enrichment analyses.
Table 3. Summary of biological processes enriched for each statistical comparison in a specific fashion.
In more detail, differential expression analysis between the match-mismatch groups LL and LH showed the largest number of DEGs (808; Supplementary Table 6), with similar numbers of these being up- (397) or down-regulated (411). The top 5 DEGs were: sult2a5 (p.adj = 2.76 x 10−39, lfc = −3.01), plin4 (p.adj = 1.44 x 10−34, lfc = −2.96), mup3 (p.adj = 1.91 x 10−33, lfc = 1.94), sult2a3 (p.adj = 6.91 x 10−31, lfc = −3.53) and gm10804 (p.adj = 2.07 x 10−22, lfc = −4.75). A great proportion of DEGs and enriched GO_BP terms were not specific for this comparison, where a particularly high overlap with the comparison between match groups (LL vs. HH; about 40% in DEGs and 56% in GO_BPs) was found (Figures 2B,C). Specific functional enrichments for the comparison LL vs. LH accounted for about 20% of the enriched terms (86/426) and included GO_BP terms related to the cardiovascular and respiratory systems, as well as to metabolism, such as carbohydrate homeostasis and responses to nutrient levels, and to cellular localization, cell adhesion and glycosylation, among others (Supplementary Table 7). In contrast, the analyses comparing HH vs. LH and LL vs. HL found only three (rap1gap: p.adj = 0.011, lfc = −0.88; cyp26b1: p.adj = 0.017, lfc = 1.64; and mt1: p.adj = 0.018, lfc = −0.98) and 19 DEGs, respectively (Supplementary Tables 2, 5). Between LL and HL, the top 5 DEGs identified were: sult2a7 (p.adj = 3.88 x 10−5, lfc = −2.88), gadd45g (p.adj = 5.13 x 10−4, lfc = 1.34), 1810008I18Rik (p.adj = 1.17 x 10−3, lfc = −0.83), srebf1 (p.adj = 1.27 x 10−3, lfc = −0.79) and lrtm1 (p.adj = 2.28x10−3, lfc = −1.16). Because in these cases we allowed functional enrichments to take place with only one overlapping gene, however, there were large numbers of significantly overrepresented GO_BP terms (40 and 104, respectively). From these, 47.5% were specific for the HH vs. LH comparison, while 41.3% were specific for the LL vs. HL comparison (Supplementary Table 7). Nevertheless, these latter enrichment results should be interpreted with caution.
Moreover, analysis of both match groups (LL vs. HH) resulted in 452 DEGs and 291 overrepresented GO_BP terms (Supplementary Tables 4, 7). The top 5 DEGs were: pde6c (p.adj = 2.24 x 10−27, lfc = −2.22), mup17 (p.adj = 2.03 x 10−25, lfc = 3.18), mup3 (p.adj = 5.67 x 10−25, lfc = 1.85), mup6 (p.adj = 2.61 x 10−23, lfc = 2.17) and plin4 (p.adj = 2.33 x 10−21, lfc = −2.92). From the GO_BP terms, 37 (12.7%) were comparison-specific and included terms such as androgen metabolism, embryonic development, heart contraction, and lipoprotein particle remodeling, among others. On the other hand, comparing both mismatch groups (HL vs. LH) resulted in 691 DEGs and 467 GO_BP terms (Supplementary Tables 3, 7), showing a particularly important overlap with the analysis of HH vs. HL (45.5% in DEGs and 61.4% in GO_BPs; Figures 2B,C). This overlap was even found within the top 5 DEGs in these comparisons with the gene mup7 (HL vs. LH: p.adj = 2.04 x 10−29, lfc = 5.03; HH vs. HL: p.adj = 4.98 x 10−19, lfc = −4.86). For the HL vs. LH comparison, 86 (18.4%) of the overrepresented GO_BP terms were specific and included terms related to metabolism, cell signaling and cell death, immunity, beta-oxidation and cholesterol storage, among others (Supplementary Table 7).
As mentioned above, many overlaps among comparisons were identified for both, DEGs and enriched GO_BP terms (Figures 2B,C). In the case of DEGs, at least four distinctive functional categories common to all comparisons were recognized among a maximum of 20 of the most significant genes obtained from each differential expression analysis. These categories were metabolism, immunity, transport of bile acids and major urinary proteins (Mups). Together, top DEGs resulting from our comparisons highlighted the importance of 64 genes in the transcriptional responses to patterns of food availability in the mouse liver. These included genes coding for various sulfotransferases (e.g. sult2a5, sult3a1, sult5a1), cytochrome P450 proteins (e.g. cyp4a14, cyp2j9), solute carriers (e.g. slc22a7, slc10a2, slco1a1), Mups (e.g. mup9, mup14, mup21), and regulatory molecules, including metabolic (e.g. mid1ip1, ppp1r3c, rarres1), immune (e.g. dntt, egr1, gadd45g) and circadian (e.g. dbp) regulators (Figure 4A). Furthermore, in the case of overrepresented GO_BPs, 156 terms were found commonly enriched in a maximum of four to five of the performed comparisons (Supplementary Table 7). These terms could be grouped into the following broad functional categories: metabolism, homeostasis, immunity, differentiation, cellular responses, localization, molecular transport, oxidation-reduction and apoptosis/cell death. Some of the terms in these categories have been manually selected to represent these top shared GO_BPs and are shown in Figure 4B. No DEGs and/or GO_BP terms overlapping among all comparisons performed were identified.
Figure 4. (A) Top genes. Selection was based on the gene being differentially expressed within the top 20 most significant hits in more than one comparison, and belonging to a recognizable cluster of functional categories. TBA: transport of bile acids. (B) Representative top shared Gene Ontology biological processes (GO_BPs). Selected terms are representative of all different GO_BP terms enriched in at least four of the comparisons. Only significant (adjusted p < 0.05) genes and functional terms are shown.
Discussion
According to the MMH, negative effects on health can occur if the environment in which an individual was shaped differs from the environment it later lives in (3–7). Noticeably, match-mismatch studies so far focused on the long-term shaping of the phenotype during prenatal and early postnatal development. Contrarily, adolescence has only recently gained attention as a time period during which shaping of the phenotype in response to the prevailing environmental conditions can occur (20, 21).
To test the MMH assuming adolescence as the sensitive period of developmental phenotypic plasticity, we used female C57BL/6J mice exposed to matching or mismatching situations of low and high food availability in adolescence and early adulthood. Here, we reported our study of differential gene expression in the liver using RNA-seq. In accordance with our previous assessments of physiological and behavioral parameters (24), our results suggested no match-mismatch effects, indicating that no long-term shaping occurred in response to food availability during adolescence. However, relatively large changes in the liver transcriptome were observed according to the amount of food available to the animals during early adulthood. These changes involved mainly a plethora of metabolic pathways and, in a lower extent, immune processes. Similarly, in our previous publication (24), various effects of low vs. high food availability were observed regarding the activation of the HPA axis (39), as well as of organ weights, including liver and spleen. Interestingly, we previously observed increased liver weights in mice of the LL group, compared to mice in the HH group, when the measurements were based on a subsample of 10 animals per group (24) whereas, in the sequencing subsample reported herein, comprised of 4 mice per group, the differences in liver weights were found between the match-mismatch LH and HL conditions, which aligned better with our expectations. Taken together, our results suggested a short-term adjustment of the metabolism and, possibly, the immune system to the prevailing situation of food availability.
In accordance with our findings, an impact on gene expression levels of the currently prevailing food availability has been previously reported. For example, a study showed that prolonged feeding of a high-fat diet to adult mice induced significant differences in the expression of hepatic genes associated with metabolism, immune system and other biological processes. In particular, several genes involved in defense responses against a lipotoxic environment, oxidative stress and inflammatory processes were upregulated (40). Furthermore, it was also shown that feeding a restricted diet (60% of ad libitum intake) decreased the expression levels in heart, liver and hypothalamus of genes associated with inflammation (41). This is in line with a chronic elevation of corticosterone release observed previously in mice and rats fed a restrictive diet (24, 42), which is hypothesized to enhance protective mechanisms against stressors (42–44) and is likely to have a profound anti-inflammatory effect (45). Correspondingly, there is compelling evidence for health improvements in rats fed with even mildly restrictive (down to 10% of ad libitum intake) diets (46–48). Taken together, the evidence suggests that high food intake is associated with inflammation and activation of the immune system [reviewed in Reilly & Saltiel (49) and Berg & Scherer (50)]. This may lead to conditions linked to chronic inflammation, such as obesity and its comorbidities, and increase the risk of cardiovascular disease [reviewed in Karczewskiet al. (51), Libby (52)]. Indeed, low-grade chronic inflammation is suggested as an underlying cause of the well-established negative health consequences of ad libitum feeding in rodents (53–57). Therefore, our study agrees with previous evidence and supports the notion that a high food intake due to unrestricted food availability might act as a risk factor for metabolic disease partially through the activation of inflammatory and immune processes while, contrarily, mild to moderate food restriction might have anti-inflammatory and immunosuppressive effects that favor healthier metabolic outcomes. However, we found no compelling evidence to support long-term metabolic effects of exposure to restricted or unrestricted diets in adolescence nor to the match-mismatch situation, suggesting that the molecular changes we observed herein were driven by acute responses to the current feeding situation.
A more detailed inspection of the top 20 DEGs from all comparisons performed in our study revealed that, beyond metabolism and immunity, other top genes belonged to the Mups and bile acid transporters. Mups are mainly expressed in the liver and secreted into the urine (58). Recently, they have been linked to metabolic diseases, such as obesity and type 2 diabetes (59). In accordance with the present study, decreasing the caloric intake of mice was shown to cause a reduction of expression levels of various Mup genes in the liver (60, 61). Such a decrease of expression levels in a situation of nutritional scarcity might reduce energy expenditure by decreasing protein loss via the urine (58). In contrast, increased expression levels when food is available in abundance might help to alleviate the negative consequences of increased food intake by utilizing unneeded energy for Mup production, as well as increasing locomotor activity and glucose tolerance (62). Moreover, Mups have been associated with social communication via pheromones in the urine [reviewed in Charkoftakiet al. (59)]. While the literature focuses on the effects of male Mups on females [e.g., (63, 64)], the role of female Mup expression on conspecifics and within the producing animal is less clear. Taken together with our previous observation of a change to social interest in an unfamiliar female caused by differences in food availability, the differential expression of Mups possibly indicates a link of the internal Mup expression pattern and the social behavior of females. Moreover, bile acid transporters have recently come into focus in the context of metabolic diseases (65). Bile acids (BA) are synthesized from cholesterol in the liver and secreted into the small intestine, where they are crucial for the absorption of dietary fatty acids, cholesterol and fat-soluble vitamins (66). In addition, they also act as signaling molecules, affecting for example bile acid metabolism, glucose homeostasis and lipid metabolism (67, 68). A disruption of the cholesterol and BA homeostasis is associated with metabolic diseases (67, 68). Bile acid transporters are likely to play an important role in maintaining BA homeostasis, as they are indispensable for the reabsorption of 95% of BA from the intestines (69). Indeed, there is evidence that a disruption of the well characterized ileal apical sodium bile acid cotransporter (Slc10a2) can cause an increased fecal loss of BA resulting in increased BA synthesis (69). A high expression level of transporters under high food availability, as observed herein, might lead to increased reuptake of BA from the intestines, in turn inhibiting the synthesis of new BA. This would lead to less absorption of cholesterol and fats, possibly alleviating the negative effects of high food intake (70).
We acknowledge some limitations of our study. Negative effects associated with a mismatch in humans, such as metabolic diseases, occur predominantly in later life stages (7). Therefore, effects on health and welfare caused by a mismatch situation in mice might not yet have manifested during early adulthood. Nonetheless, in case of a lasting effect of food availability during adolescence on late adulthood, first indications could be expected during earlier life stages. Instead, only effects of the currently prevailing food availability during early adulthood were found. Moreover, most evidence for the MMH has been gathered from situations of very low food availability, such as during a famine, compared to food in abundance (8–12, 18, 71). The discrepancies of approximately 10% between high and low food availability in our study might have been too mild to cause any long-term changes to the phenotype during adolescence.
Conclusions
We hypothesized that adolescence might represent a sensitive period for the long-term metabolic shaping. As such, we expected to observe physiological, behavioral as well as molecular differences in animals that were exposed to different schemes of mildly restrictive (10% for 11 weeks)/unrestrictive food availability during adolescence compared to early adulthood. Nevertheless, results from either our previous physiological and behavioral studies or the current molecular study support a match-mismatch effect between adolescence and early adulthood when a mildly restrictive diet is implemented. Taken together, physiological, behavioral and molecular differences observed in the liver transcriptome between our study groups mainly pointed toward short-term changes in metabolic and immune pathways. These results might have implications for the development of obesity and other metabolic diseases. Further work should explore the possibility of a match-mismatch effect between adolescence and (early) adulthood at higher dietary restrictions and systematically assess whether dietary restriction in adolescence results beneficial for the health status later in life.
Data Availability Statement
Raw and processed RNA-seq data has been deposited into NCBI's Gene Expression Omnibus (GEO) database and is accessible with the ID GSE200731.
Ethics Statement
The animal study was reviewed and approved by Gesundheits-und Veterinäramt Münster, Nordrhein-Westfalen and Landesamt für Natur, Umwelt und Verbraucherschutz Nordrhein-Westfalen.
Author Contributions
MS, SK, SR, and NS conceived the study and designed the experiments. JF-D carried out the experimental procedures and sample collection, and prepared the manuscript. MH-R processed and analyzed the data and prepared the manuscript. AW supervised sequencing. All authors reviewed and approved the final version of the manuscript.
Funding
This work was supported by a grant from the German Research Foundation (DFG) to NS (281125614/GRK2220, Project B5).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors would like to thank Elvira Barg and Marianne Jansen-Rust from the Core Facility Genomics of the Medical Faculty at the University of Münster, and Kimberly Kubski and Hanna Sicking from the Department of Behavioral Biology at the University of Münster, for their excellent technical assistance.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2022.910762/full#supplementary-material
Abbreviations
BA, bile acids; DEGs, differentially expressed genes; GO_BP, Gene Ontology biological processes; H, high food availability; HH, high/high food availability match group; HL, high/low food availability mismatch group; L, low food availability; LH, low/high food availability mismatch group; LL, low/low food availability match group; MMH, Match-Mismatch hypothesis; Mups, major urinary proteins; PCs, principal components; PND, postnatal day; RNA, ribonucleic acid; RNA-seq, RNA-sequencing; SD, standard deviation.
References
1. Bateson P. Fetal experience and good adult design. Int J Epidemiol. (2001) 30:928–34. doi: 10.1093/ije/30.5.928
2. Gluckman PD, Hanson MA. Developmental origins of disease paradigm: a mechanistic and evolutionary perspective. Pediatr Res. (2004) 56:311–7. doi: 10.1203/01.PDR.0000135998.08025.FB
3. Bateson P, Gluckman P, Hanson M. The biology of developmental plasticity and the Predictive Adaptive Response hypothesis. J Physiol. (2014) 592:2357–68. doi: 10.1113/jphysiol.2014.271460
4. Bateson P, Barker D, Clutton-Brock T, Deb D, D'Udine B, Foley RA, et al. Developmental plasticity and human health. Nature. (2004) 430:419–21. doi: 10.1038/nature02725
5. Gluckman PD, Cutfield W, Hofman P, Hanson MA. The fetal, neonatal, and infant environments—the long-term consequences for disease risk. Early Hum Dev. (2005) 81:51–9. doi: 10.1016/j.earlhumdev.2004.10.003
6. Gluckman PD, Hanson MA, Pinal C. The developmental origins of adult disease. Matern Child Nutr. (2005) 1:130–41. doi: 10.1111/j.1740-8709.2005.00020.x
7. Gluckman PD, Hanson MA, Beedle AS. Early life events and their consequences for later disease: a life history and evolutionary perspective. Am J Hum Biol. (2007) 19:1–19. doi: 10.1002/ajhb.20590
8. Lindeboom M, Portrait F, van den Berg GJ. Long-run effects on longevity of a nutritional shock early in life: The Dutch Potato famine of 1846–1847. J Health Econ. (2010) 29:617–29. doi: 10.1016/j.jhealeco.2010.06.001
9. Roseboom TJ. Prenatal exposure to famine and aging. In: AV, editor. Early Life Origins of Aging and Longevity Healthy Aging and Longevity. Healthy Aging and Longevity. Cham;: Springer (2019). p. 233–44. doi: 10.1007/978-3-030-24958-8_12
10. Roseboom TJ, Painter RC, van Abeelen AFM, Veenendaal MVE, de Rooij SR. Hungry in the womb: what are the consequences? Lessons from the Dutch famine. Maturitas. (2011) 70:141–5. doi: 10.1016/j.maturitas.2011.06.017
11. Roseboom T, de Rooij S, Painter R. The Dutch famine and its long-term consequences for adult health. Early Hum Dev. (2006) 82:485–91. doi: 10.1016/j.earlhumdev.2006.07.001
12. Roseboom TJ, Van Der Meulen JHP, Ravelli ACJ, Osmond C, Barker DJP, Bleker OP. Effects of prenatal exposure to the Dutch famine on adult disease in later life: an overview. Twin Res. (2001) 4:293–8. doi: 10.1375/twin.4.5.293
13. Chali D, Enquselassie F, Gesese MA. case-control study on determinants of rickets. Ethiop Med J. (1998) 36:227–34.
14. Forrester TE, Badaloo AV, Boyne MS, Osmond C, Thompson D, Green C, et al. Prenatal factors contribute to the emergence of kwashiorkor or marasmus in severe undernutrition: evidence for the predictive adaptation model. PLoS ONE. (2012) 7:e35907. doi: 10.1371/journal.pone.0035907
15. Barker DJP, Forsén T, Uutela A, Osmond C, Eriksson JG. Size at birth and resilience to effects of poor living conditions in adult life: longitudinal study. Br Med J. (2001) 323:1273. doi: 10.1136/bmj.323.7324.1273
16. Eriksson JG, Forsen T, Tuomilehto J, Osmond C, Barker DJ. Early growth and coronary heart disease in later life: longitudinal study. Br Med J. (2001) 322:949–53. doi: 10.1136/bmj.322.7292.949
17. Gluckman PD, Hanson MA, Spencer HG, Bateson P. Environmental influences during development and their later consequences for health and disease: implications for the interpretation of empirical studies. Proc R Soc B Biol Sci. (2005) 272:671–7. doi: 10.1098/rspb.2004.3001
18. Heijmans BT, Tobi EW, Stein AD, Putter H, Blauw GJ, Susser ES, et al. Persistent epigenetic differences associated with prenatal exposure to famine in humans. Proc Nat Acad Sci. (2008) 105:17046–9. doi: 10.1073/pnas.0806560105
19. Stein AD, Zybert PA, van de Bor M, Lumey LH. Intrauterine famine exposure and body proportions at birth: the Dutch Hunger Winter. Int J Epidemiol. (2004) 33:831–6. doi: 10.1093/ije/dyh083
20. Sachser N, Kaiser S, Hennessy MB. Behavioral profiles are shaped by social experience: when, how and why. Philoso Transac R Soc Lond B Biol Sci. (2013) 368:20120344. doi: 10.1098/rstb.2012.0344
21. Sachser N, Hennessy MB, Kaiser S. The adaptive shaping of social behavioral phenotypes during adolescence. Biol Lett. (2018) 14:20180536. doi: 10.1098/rsbl.2018.0536
22. Romeo RD. Pubertal maturation and programming of hypothalamic–pituitary–adrenal reactivity. Front Neuroendocrinol. (2010) 31:232–40. doi: 10.1016/j.yfrne.2010.02.004
23. Busch DS, Hayward LS. Stress in a conservation context: a discussion of glucocorticoid actions and how levels change with conservation-relevant variables. Biol Conserv. (2009) 142:2844–53. doi: 10.1016/j.biocon.2009.08.013
24. Feige-Diller J, Palme R, Kaiser S, Sachser N, Richter SH. The impact of varying food availability on health and welfare in mice: Testing the Match-Mismatch hypothesis. Physiol Behav. (2021) 228:113193. doi: 10.1016/j.physbeh.2020.113193
25. Rui L. Energy metabolism in the liver. Compr Physiol. (2014) 4:177–97. doi: 10.1002/cphy.c130024
26. Robinson MW, Harmon C, O'Farrelly C. Liver immunology and its role in inflammation and homeostasis. Cell Mol Immunol. (2016) 13:267–76. doi: 10.1038/cmi.2016.3
27. Crispe IN. The liver as a lymphoid organ. Annu Rev Immunol. (2009) 27:147–63. doi: 10.1146/annurev.immunol.021908.132629
28. Nemeth E, Baird AW, O'Farrelly C, editors. Microanatomy of the liver immune system. Sem Immunopatho. (2009) 33:333–43. doi: 10.1007/s00281-009-0173-4
29. Hashimoto T, Watanabe S. Chronic food restriction enhances memory in mice–analysis with matched drive levels. Neuroreport. (2005) 16:1129–33. doi: 10.1097/00001756-200507130-00019
30. Kliewer KL, Ke JY, Lee HY, Stout MB, Cole RM, Samuel VT, et al. Short-term food restriction followed by controlled refeeding promotes gorging behavior, enhances fat deposition, and diminishes insulin sensitivity in mice. J Nutr Biochem. (2015) 26:721–8. doi: 10.1016/j.jnutbio.2015.01.010
31. Andrews S. Babraham Bioinformatics-Fastqc: A Quality Control Tool for High Throughput Sequence Data. (2010). Available online at: http://www.bioinformatics.babraham.ac.uk/projects/fastqc/.
32. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. (2014) 30:2114–20. doi: 10.1093/bioinformatics/btu170
33. Kim D, Paggi JM, Park C, Bennett C, Salzberg SL. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol. (2019) 37:907–15. doi: 10.1038/s41587-019-0201-4
34. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAM tools. Bioinformatics. (2009) 25:2078–9. doi: 10.1093/bioinformatics/btp352
35. Anders S, Pyl PT, Huber W. HTSeq—a Python framework to work with high-throughput sequencing data. Bioinformatics. (2015) 31:166–9. doi: 10.1093/bioinformatics/btu638
36. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. (2014) 15:550. doi: 10.1186/s13059-014-0550-8
37. Maere S, Heymans K, Kuiper M. BiNGO a Cytoscape plugin to assess overrepresentation of Gene Ontology categories in Biological Networks. Bioinformatics. (2005) 21:3448–9. doi: 10.1093/bioinformatics/bti551
38. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. (2003) 13:2498–504. doi: 10.1101/gr.1239303
39. Palme R. Non-invasive measurement of glucocorticoids: advances and problems. Physiol Behav. (2019) 199:229–43. doi: 10.1016/j.physbeh.2018.11.021
40. Kim S, Sohn I, Ahn J-I, Lee K-H, Lee YS, Lee YS. Hepatic gene expression profiles in a long-term high-fat diet-induced obesity mouse model. Gene. (2004) 340:99–109. doi: 10.1016/j.gene.2004.06.015
41. Fu C, Hickey M, Morrison M, McCarter R, Han E-S. Tissue specific and non-specific changes in gene expression by aging and by early stage CR. Mech Aging Dev. (2006) 127:905–16. doi: 10.1016/j.mad.2006.09.006
42. Levay EA, Tammer AH, Penman J, Kent S, Paolini AG. Calorie restriction at increasing levels leads to augmented concentrations of corticosterone and decreasing concentrations of testosterone in rats. Nutr Res. (2010) 30:366–73. doi: 10.1016/j.nutres.2010.05.001
43. Yu BP, Chung HY. Stress resistance by caloric restriction for longevity. Ann N Y Acad Sci. (2001) 928:39–47. doi: 10.1111/j.1749-6632.2001.tb05633.x
44. Mattson MP, Chan SL, Duan W. Modification of brain aging and neurodegenerative disorders by genes, diet, and behavior. Physiol Rev. (2002) 82:637–72. doi: 10.1152/physrev.00004.2002
45. Masoro EJ. Overview of caloric restriction and aging. Mech Aging Dev. (2005) 126:913–22. doi: 10.1016/j.mad.2005.03.012
46. Duffy PH, Seng JE, Lewis SM, Mayhugh MA, Aidoo A, Hattan DG, et al. The effects of different levels of dietary restriction on aging and survival in the Sprague-Dawley rat: implications for chronic studies. Aging Clin Exp Res. (2001) 13:263–72. doi: 10.1007/BF03353422
47. Duffy PH, Lewis SM, Mayhugh MA, Trotter RW, Latendresse JR, Thorn BT, et al. The effects of different levels of dietary restriction on neoplastic pathology in the male Sprague-Dawley rat. Aging Clin Exp Res. (2004) 16:448–56. doi: 10.1007/BF03327400
48. Duffy PH, Lewis SM, Mayhugh MA, Trotter RW, Thorn BT, Feuers RJ, et al. The effects of different levels of dietary restriction on non-neoplastic diseases in male Sprague-Dawley rats. Aging Clin Exp Res. (2004) 6:68–78. doi: 10.1007/BF03324535
49. Reilly SM, Saltiel AR. Adapting to obesity with adipose tissue inflammation. Nature Rev Endocrinol. (2017) 13:633. doi: 10.1038/nrendo.2017.90
50. Berg AH, Scherer PE. Adipose tissue, inflammation, and cardiovascular disease. Circ Res. (2005) 96:939–49. doi: 10.1161/01.RES.0000163635.62927.34
51. Karczewski J, Sledzińska E, Baturo A, Jończyk I, Maleszko A, Samborski P, et al. Obesity and inflammation. Eur Cytokine Netw. (2018) 29:83–94. doi: 10.1684/ecn.2018.0415
52. Libby P. Inflammation and cardiovascular disease mechanisms. Am J Clin Nutr. (2006) 83:456S–60S. doi: 10.1093/ajcn/83.2.456S
53. Keenan KP, Laroque P, Ballam GC, Soper KA, Dixit R, Mattson BA, et al. The effects of diet, ad libitum overfeeding, and moderate dietary restriction on the rodent bioassay: the uncontrolled variable in safety assessment. Toxicol Pathol. (1996) 24:757–68. doi: 10.1177/019262339602400620
54. Keenan KP C, Ballam G, Haught DG, Laroque P. Nutrition In: Krinke GJ, editor. The Laboratory Rat Academic Press: Elsevier (2000). doi: 10.1016/B978-012426400-7.50044-3
55. Keenan KP, Hoe C-M, Mixson L, Mccoy CL, Coleman JB, Mattson BA, et al. Diabesity: a polygenic model of dietary-induced obesity from ad libitum overfeeding of Sprague–Dawley rats and its modulation by moderate and marked dietary restriction. Toxicol Pathol. (2005) 33:650–74. doi: 10.1080/01926230500311222
56. Roe FJC, Lee PN, Conybeare G, Kelly D, Matter B, Prentice D, et al. The Biosure Study: influence of composition of diet and food consumption on longevity, degenerative diseases and neoplasia in Wistar rats studied for up to 30 months post weaning. Food Chem Toxicol. (1995) 33:S1–S100. doi: 10.1016/0278-6915(95)80200-2
57. Roe FJC. 1200-Rat Biosure Study: Desing and Overview of Results. In: Fishbein L, editor. Biological Effects of Dietary Restriction. Heidelberg: Springer (1991). doi: 10.1007/978-3-642-58181-6_26
58. Beynon RJ, Hurst JL, Gaskell SJ, Hubbard SJ, Humphries RE, Malone N, et al. Mice, MUPs and myths: structure-function relationships of the major urinary proteins. Chemical Signals. (2001) 9:149–56. doi: 10.1007/978-1-4615-0671-3_19
59. Charkoftaki G, Wang Y, McAndrews M, Bruford EA, Thompson DC, Vasiliou V, et al. Update on the human and mouse lipocalin (LCN) gene family, including evidence the mouse Mup cluster is result of an “evolutionary bloom”. Hum Genomics. (2019) 13:1–14. doi: 10.1186/s40246-019-0191-9
60. Giller K, Huebbe P, Doering F, Pallauf K, Rimbach G. Major urinary protein 5, a scent communication protein, is regulated by dietary restriction and subsequent re-feeding in mice. Proc R Soc B Biol Sci. (2013) 280:20130101. doi: 10.1098/rspb.2013.0101
61. Dhahbi JM, Kim H-J, Mote PL, Beaver RJ, Spindler SR. Temporal linkage between the phenotypic and genomic responses to caloric restriction. Proc Nat Acad Sci. (2004) 101:5524–9. doi: 10.1073/pnas.0305300101
62. Hui X, Zhu W, Wang Y, Lam KS, Zhang J, Wu D, et al. Major urinary protein-1 increases energy expenditure and improves glucose intolerance through enhancing mitochondrial function in skeletal muscle of diabetic mice. J Biol Chem. (2009) 284:14050–7. doi: 10.1074/jbc.M109.001107
63. Mucignat-Caretta C, Caretta A, Cavaggioni A. Acceleration of puberty onset in female mice by male urinary proteins. J Physiol. (1995) 486:517–22. doi: 10.1113/jphysiol.1995.sp020830
64. Marchlewska-Koj A, Cavaggioni A, Mucignat-Caretta C, Olejniczak P. Stimulation of estrus in female mice by male urinary proteins. J Chem Ecol. (2000) 26:2355–66.
65. Yang N, Dong Y-Q, Jia G-X, Fan S-M, Li S-Z, Yang S-S, et al. ASBT(SLC10A2): A promising target for treatment of diseases and drug discovery. Biomed Pharmacother. (2020) 132:110835. doi: 10.1016/j.biopha.2020.110835
66. Lundåsen T, Andersson E-M, Snaith M, Lindmark H, Lundberg J, Östlund-Lindqvist A-M, et al. Inhibition of Intestinal Bile Acid Transporter Slc10a2 improves Triglyceride metabolism and normalizes elevated plasma glucose levels in mice. PLoS One. (2012) 7:e37787. doi: 10.1371/journal.pone.0037787
67. Li T, Chiang JYL. Bile acid signaling in liver metabolism and diseases. J Lipids. (2012) 2012:754067. doi: 10.1155/2012/754067
68. Li T, Chiang JYL. Bile acid signaling in metabolic disease and drug therapy. Pharmacol Rev. (2014) 66:948–83. doi: 10.1124/pr.113.008201
69. Dawson PA, Haywood J, Craddock AL, Wilson M, Tietjen M, Kluckman K, et al. Targeted deletion of the ileal bile acid transporter eliminates enterohepatic cycling of bile acids in mice. J Biol Chem. (2003) 278:33920–7. doi: 10.1074/jbc.M306370200
70. La Frano MR, Hernandez-Carretero A, Weber N, Borkowski K, Pedersen TL, Osborn O, et al. Diet-induced obesity and weight loss alter bile acid concentrations and bile acid–sensitive gene expression in insulin target tissues of C57BL/6J mice. Nutr Res. (2017) 46:11–21. doi: 10.1016/j.nutres.2017.07.006
Keywords: Match-Mismatch hypothesis, RNA-sequencing (RNA-seq), food availability changes, differential expression (DE) analysis, liver transcriptome, metabolism, immunity and inflammation
Citation: Feige-Diller J, Herrera-Rivero M, Witten A, Stoll M, Kaiser S, Richter SH and Sachser N (2022) The Impact of Varying Food Availability on Gene Expression in the Liver: Testing the Match-Mismatch Hypothesis. Front. Nutr. 9:910762. doi: 10.3389/fnut.2022.910762
Received: 01 April 2022; Accepted: 02 June 2022;
Published: 01 July 2022.
Edited by:
Safarina G. Malik, Mochtar Riady Institute for Nanotechnology, IndonesiaReviewed by:
Sung Wook Kang, Baylor University, United StatesYuqian Ye, National Cancer Center, China
Copyright © 2022 Feige-Diller, Herrera-Rivero, Witten, Stoll, Kaiser, Richter and Sachser. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Marisol Herrera-Rivero, Marisol.Herrera@uni-muenster.de; Marisol.Herrerarivero@ukmuenster.de; Janina Feige-Diller, j.feige-diller@uni-muenster.de
†These authors have contributed equally to this work and share first authorship