- Open Access
Starvation causes changes in the intestinal transcriptome and microbiome that are reversed upon refeeding
BMC Genomics volume 23, Article number: 225 (2022)
The ability of animals and their microbiomes to adapt to starvation and then restore homeostasis after refeeding is fundamental to their continued survival and symbiosis. The intestine is the primary site of nutrient absorption and microbiome interaction, however our understanding of intestinal adaptations to starvation and refeeding remains limited. Here we used RNA sequencing and 16S rRNA gene sequencing to uncover changes in the intestinal transcriptome and microbiome of zebrafish subjected to long-term starvation and refeeding compared to continuously fed controls.
Starvation over 21 days led to increased diversity and altered composition in the intestinal microbiome compared to fed controls, including relative increases in Vibrio and reductions in Plesiomonas bacteria. Starvation also led to significant alterations in host gene expression in the intestine, with distinct pathways affected at early and late stages of starvation. This included increases in the expression of ribosome biogenesis genes early in starvation, followed by decreased expression of genes involved in antiviral immunity and lipid transport at later stages. These effects of starvation on the host transcriptome and microbiome were almost completely restored within 3 days after refeeding. Comparison with published datasets identified host genes responsive to starvation as well as high-fat feeding or microbiome colonization, and predicted host transcription factors that may be involved in starvation response.
Long-term starvation induces progressive changes in microbiome composition and host gene expression in the zebrafish intestine, and these changes are rapidly reversed after refeeding. Our identification of bacterial taxa, host genes and host pathways involved in this response provides a framework for future investigation of the physiological and ecological mechanisms underlying intestinal adaptations to food restriction.
Starvation is a state of severe caloric restriction regularly experienced by many animal species and a significant portion of the human population. In humans, starvation can be the result of environmental or socioeconomic conditions including war, famine, and poverty . It can also occur alongside pathologies such as anorexia nervosa and cancer . In animals, periods of absolute or relative starvation can result from seasonal changes such as drought and severe cold, or from behaviors such as nesting, lactation, migration, and hibernation . This wide range of circumstances leading to starvation across the animal kingdom evokes a range of progressive physiologic adaptations to starvation across different species. Indeed, previous studies have reported similarity and divergence in starvation physiology across animal taxa such as humans, rodents, polar bears, penguins, reptiles, amphibians, fish, and insects . However, previous studies have largely focused on tissue histopathologies associated with starvation, whereas effects on the underlying physiological processes remain incompletely understood.
Across many animal species, starvation leads to a progressive decrease in metabolic rate . Increased blood glycerol, which serves as a gluconeogenic precursor, is also common in starved animals, as are fluctuations in free fatty acids [3, 4]. The overall depletion in energy stores leads to weight loss, which is generally greater in endotherms when compared to ectotherms . Starvation is also associated with a gradual reduction in mass in important organs such as the liver, skeletal muscle, and intestine [4, 7]. These effects necessitate a recovery from starvation to restore optimal function to these organs. Inherently linked to starvation, the return to homeostasis following starvation is facilitated by a refeeding response that gradually reverses starvation-induced adaptations and restores energy balance. Physiological responses to starvation and subsequent refeeding are dynamic and complex, involving coordination between major organ systems via nutritional and hormonal signals. The ultimate outcome of these starvation responses is often the preservation of lean body mass while favoring the depletion of energy stores such as glycogen and fat [8, 9]. However, despite these effects, starvation often results in lasting defects on bone density, pancreatic function, and mental development long after refeeding [10,11,12]. Thus, improved understanding of these dynamic physiological processes could lead to new approaches to reduce morbidities and mortalities associated with starvation in humans and other animals .
Previous studies on the effects of refeeding after starvation have largely focused on tissues such as liver, skeletal muscle, brain, and pancreatic islets [14,15,16,17,18]. We have a relatively poor understanding of the transcriptional starvation and refeeding responses in the intestine. The intestine is the major site of dietary nutrient sensing and absorption, and harbors complex communities of microorganisms (microbiome). Previous studies in humans and rodent models have shown that intestinal microbiome composition changes in response to starvation and diet composition with distinct contributions to the nutritional physiology of their hosts [19,20,21,22,23,24,25,26]. These findings informed more recent studies that have investigated microbiome-targeted therapeutics for alleviating starvation and its associated developmental defects [27,28,29]. However, gut microbial responses to starvation have been largely limited to mammals, and our understanding of intestinal physiological responses to starvation and feeding in any animal remains quite limited.
Animal models provide opportunities to study the processes that underlie starvation and refeeding responses in vertebrates, resulting in a general understanding that may be translated to humans . Poikilothermic vertebrates such as cyprinid fishes are particularly interesting due to their capacity to endure prolonged starvation periods. In response to prolonged starvation, cyprinids such as carp exhibit a reduction in intestinal thickness and weight, altered enterocyte morphology, and a decrease in body weight and liver size, similar to the starvation response in mice [17, 31,32,33,34,35]. Zebrafish (Danio rerio) survive up to 4 weeks of starvation as adults, and a suite of genomic and genetic resources facilitate the investigation of their physiology . Using in vivo imaging to monitor white adipose tissues as a measure of energy storage, we previously showed that prolonged starvation in adult zebrafish leads to progressive mobilization of fat stored in white adipose tissues, which is replenished in response to refeeding [37,38,39]. Because adipose tissues develop progressively during juvenile and adult stages, the duration of starvation required to completely mobilize adipose lipid increases with animal age (e.g., from 1 week in juveniles up to 3 weeks in adults) [37,38,39]. However, the impact of prolonged starvation and refeeding on the zebrafish intestine has not been explored.
The zebrafish intestine displays extensive cellular and physiological homology to that of mammals, and harbors a microbiome that varies in composition as a function of age and diet composition [40,41,42,43,44,45]. The presence and composition of the intestinal microbiome in zebrafish impacts the host by regulating dietary nutrient absorption, epithelial renewal, and inflammation [41, 43, 46,47,48,49,50]. By comparing patterns of gene expression and accessible chromatin in intestinal epithelial cells from zebrafish, mouse, and human, we recently discovered a conserved transcriptional regulatory network conserved across 420 million years of vertebrate evolution . Building upon this recent work, here we define the impact of prolonged starvation and refeeding on gene expression in the adult zebrafish intestine, and use these results to predict the physiological processes and transcription regulatory pathways that underlie the response to starvation and refeeding. We also show how the taxonomic composition of the adult zebrafish intestinal microbiome is altered during the same prolonged starvation and refeeding regimen.
Starvation is accompanied by significant changes in gut microbiome composition that are reversed during refeeding
To determine the influence of starvation on zebrafish microbiome composition and intestinal gene expression, zebrafish were reared under conventional conditions using a standard diet to early adulthood (60 days post fertilization or dpf). Animals were then moved into clean tanks and randomly assigned into one of two treatment groups: one group was starved for 21 days followed by 21 days refeeding, and a reference control group was consistently fed across the same 42-day time course (Fig. 1A). This 21 day starvation regimen was selected because it is sufficient in adult zebrafish to completely deplete stored lipid from adipose tissues and reduce body weight and liver size, whereas subsequent refeeding largely restores total adipose tissue lipid, body weight, and liver size within 14 days [37, 38, 52, 53]. We performed 16S rRNA gene sequencing on whole intestinal samples from zebrafish at 0, 1, 3, 7, and 21 days post-starvation (dpS), then at 1, 3, 7, and 21 days post-refeeding (dpR) with a standard diet (Fig. 1A). Age matched siblings fed the same standard diet on a daily basis served as reference controls and were sampled at the same time points.
During these early adult stages, zebrafish fed normally continued to display somatic growth as expected. Measurements of animal size as standard length (SL) and height at anterior anal fin (HAA) revealed that somatic growth in starved fish was largely arrested compared to control fish (Fig. 1B-C). Starved fish were significantly smaller than fed fish by 7 days post-starvation and this trend continued beyond the end of the starvation period. Starved animals resumed growth after refeeding, though they remained significantly smaller than fed fish throughout the duration of the experiment (p < 0.05, two-way ANOVA with Bonferroni correction) (Fig. 1B-C). We observed no mortality in any of these conditions consistent with previous studies [38, 53]. Starvation therefore caused a general arrest in somatic growth which was restored upon refeeding.
Analysis of 16S rRNA gene sequencing data from intestinal samples revealed the impacts of prolonged starvation and refeeding on intestinal microbiome composition. Overall the intestinal microbiomes of starved zebrafish maintained a higher Faith’s PD diversity compared to fed controls. Both conditions displayed an initial loss of diversity by 1dpS, perhaps due to stress caused by tank transfer at 0dpS when the experiment began. However starved communities maintained significant higher diversity from 1dpS through the end of starvation at 21dpS and again at 21dpR (p < 0.05, ANOVA and Tukey HSD) (Fig. 1D). Beta diversity analysis of community composition using weighted UniFrac distances showed that starved and fed communities began to differ by 1 dpS, with the centroid distances being greatest at 3dpS and 21 dpS (Fig. 2A). During refeeding, the starved fish samples quickly returned to a composition more similar to fed controls (Fig. 2A). PERMANOVA further confirmed that starvation, experimental time point, and the interaction are significant factors (p < 0.05, R2 = 0.017, 0.17 and 0.07, respectively) affecting gut microbiome composition. Thus, prolonged starvation induced detectable shifts in overall composition of gut bacterial communities that were reversed quickly after refeeding.
Starvation increases similarity between zebrafish gut microbiomes and their surrounding water environment
The microbiome residing in the zebrafish intestine exists in continuity with that of the surrounding water environment, however these communities typically display distinct compositions [42,43,44]. The ecological processes contributing to these differences remain unclear, but could include non-neutral processes such as host selection [54, 55] or the magnitude of dispersal between the intestine and the surrounding environment . To test if starvation and refeeding alter the relationship between the gut and environmental microbiomes, we compared weighted UniFrac distance between matched gut and environmental samples in starved/refed and control fish normalized by the day 0 values. The distance between gut and environmental samples increased between 0 dpS and 21 dpS (Fig. 1E), perhaps reflecting restoration of homeostasis after the stress of transfer into new tanks that occurred at 0 dpS. However, from 7 dpS to 7 dpR the distance between gut and environmental samples was greater for starved than fed controls. This suggests that the fed controls had an intestinal microbiome composition more similar to the environment compared to starved fish during those stages. At 21 dpR these distances were not statistically different between treatment groups. These results indicate that starvation increases similarity between gut and environmental communities but refeeding restores differences between these communities to levels achieved under constant feeding conditions.
Vibrio bacteria are significantly enriched in the intestine during starvation
We next sought to identify the specific bacterial taxa that were significantly affected by starvation and refeeding using LEfSe . LEfSe identified 120 genus-level taxa that reached a logarithmic LDA score of 2.0 (Table S1). This set of affected taxa included 12 abundant genera (median relative abundance > 0.1% across all samples; shown in Fig. 2B) including starvation-induced depletion of Plesiomonas and enrichment of Vibrio. Strikingly, Vibrio reached a maximum relative abundance of 87% (median 34%) in the intestines of starved zebrafish at 21dpS, which was markedly higher than that of fed controls (maximum 62%, median 4%) at the same time point (Fig. 2C). In contrast, starvation led to reduced relative abundance of Plesiomonas in the intestine by 1dpS continuing through 7dpR (Fig. 2E). These effects of starvation on Vibrio and Plesiomonas sp. in the starved guts were reflected in the environmental samples of the starved fish (Fig. 2D-F). Importantly, none of the phyla or orders that were significantly depleted or enriched were significantly correlated with SL after Bonferroni correction (see Table S2). This suggests that their depletion and enrichment are due to the dietary treatment, and not simply the growth arrest observed in starved animals (Fig. 1C-D). These results establish that prolonged starvation leads to significant alterations in intestinal microbiome composition including marked enrichment of Vibrio genus members, and that these alterations in microbiome composition are largely normalized within 1–3 days of refeeding [41, 42, 48].
Starvation and refeeding leads to distinct changes in intestinal gene expression that vary with the duration of starvation
Previous work in vertebrates has shown that starvation can significantly affect host gene expression in multiple organs [4, 7, 17, 33, 52]. Our analysis of intestinal microbiomes during starvation and refeeding suggested distinct stages - early starvation when microbiome effects are initially observed (i.e., 3 dpS), late starvation when microbiome alterations are greatest (i.e., 21 dpS), and early refeeding when microbiome composition is largely normalized (i.e., 3 dpR). We therefore dissected whole intestinal tracts from 3dpS, 21dpS, and 3dpR adult zebrafish and their fed age-matched controls for RNA-seq analysis (3–4 biological replicate samples/condition; Fig. 1A). Principal components analysis (PCA) of these data revealed similarities between biological replicates (Fig. 3A). Similar to the observed effects on the intestinal microbiome, PCA indicated the impact of starvation on intestinal gene expression was greater at 21dpS than 3dpS or 3dpR. We then used DEseq2 analysis to identify genes differentially expressed in starved/refed fish compared to their fed controls at each timepoint (Table S3A). In accord with our PCA analysis, the number of significant differentially expressed genes increased from 87 genes at 3dpS to 182 genes at 21dpS, and was reduced to 11 genes by 3dpR (Fig. 3B). This further supports that starvation has a progressive impact on the intestinal transcriptome through 3dpS and 21dpS which is largely normalized by 3dpR.
As a control and to estimate the influence of the developmental time covered during the experiment, we compared differential gene expression between the fed timepoints. There were 135 genes found to be significant in these comparisons, but only 30 of these were represented among significant genes in our Starved vs. Fed comparisons (18 at 3dpS, 15 at 21dpS, 2 at 3dpR; listed in Table S3B and shown in Fig. S1B) and they were removed from our subsequent analyses of starvation effects (see Fig. S2 for results with those genes retained in the analysis). Hierarchical clustering of the log2 fold changes in transcript abundance revealed distinct groups of genes that were upregulated or downregulated in response to starvation, including striking differences between the response to starvation at 3dpS and 21dpS (Fig. 3C). Though gene expression differences between starved/refed fish and fed controls was largely restored by 3dpR, there was a small set of 9 genes that continued to be differentially expressed even at 3dpR (Fig. 3B; Table S3A). Although this list of persistent genes was too small to permit functional enrichment analysis, it does suggest a potentially small group of intestinal functions that remain altered after restoration of feeding or that respond to both starvation as well as to refeeding. These persistently different genes are discussed below in context.
Among the genes known to be starvation responsive, we first examined transcript levels of elovl2, a fatty acid elongase previously shown to be downregulated in zebrafish during starvation [58,59,60]. Elovl2 has also been implicated in inducing insulin secretion in response to glucose in mice, and fatty acid elongases have been extensively studied in fish as they function in biosynthesis of long-chain polyunsaturated fatty acids, which are commercially important in fish aquaculture [58, 61, 62]. In accord, elovl2 was significantly downregulated by 3dpS, was one of the most significantly downregulated genes in starved fish at 21dpS, and was also consistently expressed across the control fed fish group (Fig. 3D-E, Fig. S1A). This downregulation suggested a reduction in intestinal fatty acid synthesis during starvation.
To understand which biological processes are impacted by starvation, we performed Gene Ontology (GO) term searches of four groups of genes from our dataset; genes significantly upregulated at 3dpS or 21dpS, and genes significantly downregulated at 3dpS or 21dpS (Figs. 3C, S1C-D, and S3). We first identified distinct, non-overlapping functions that were enriched early in starvation (i.e., at 3dpS) and late in starvation (i.e., at 21dpS). For example, functions enriched only among upregulated genes at 3dpS, and not 21dpS, included “ribosome” and “ribosome large subunit biogenesis” (Fig. S1E). These included the ribosome biogenesis factor nsa2 and gtpbp4 which is involved in biogenesis of the 60S ribosomal subunit, which were significantly upregulated at 3dpS but not 21dpS. However, the function “ribosomal large subunit assembly” was enriched among genes upregulated at 21dpS and not 3dpS. This included some genes that were only significantly increased at 21dpS such as ruvbl1 and srfbp1, and others that were significantly increased at both 3dpS and 21dpS such as rsl24d1, ptges3l, and gltscr2 (Table S3A). Overall, genes involved in ribosome biogenesis were induced more strongly at 3dpS compared to 21dpS, suggesting it is a relatively early response to starvation with aspects that continue through 21dpS (Fig. S3CD, Table S3A).
Also among the genes upregulated specifically at 3dpS was the heat shock protein hsp90ab1, a molecular chaperone previously shown to be upregulated in adult zebrafish liver in response to starvation . The most significantly upregulated gene in our dataset was the enteropeptidase/enterokinase tmprss15, that converts trypsinogen into active trypsin which in turn activates pancreatic enzymes and potentially also antimicrobial proteins in the intestinal lumen (Fig. 3D-F) [63, 64]. Intestinal expression of tmprss15 was not affected by starvation at 3dpS, but was upregulated 6-fold by 21dpS (Fig. 3D-F). Notably, a deficiency in TMPRSS15 has been shown to confer a lean, starvation-like phenotype in humans, consistent with its known key role in nutrient digestion and absorption . Upregulation of tmprss15 in the starved zebrafish intestine suggests potential adaptive increases in nutrient digestion programs to salvage nutrients from the intestinal lumen, or in antimicrobial defense against an altered and potentially pro-inflammatory microbiome.
Similarly, functions enriched specifically among downregulated genes at 3dpS and not at 21dpS included “metabolism of lipids”, “regulation of cell proliferation”, and “ubiquitin-dependent protein catabolic process” (Fig. S1E). This included downregulation of the acyltransferase lclat1 and fatty acid binding protein fabp1b.1 at 3dpS but not 21dpS. However, related functions “glycerophospholipid metabolism”, “lipid transport”, and “lipid metabolic process” were enriched among downregulated genes at 21dpS but not 3dpS. These included the phospholipase pla2g12b, the fatty acid desaturase fads2, the lipid transfer protein scp2a, and multiple apolipoproteins including apoa1a, apoa4b.1, and apobb.1 (Figs. 3D, S1E, S3AB) . Notably, apolipoprotein genes have been shown to be downregulated in starved rainbow trout livers . Although these genes were significantly downregulated only at 21dpS, most began trending towards downregulation at 3dpS. Yet other genes involved in these functions including elovl2 and fabp1b.1 were significantly downregulated at both 3dpS and 21dpS with a larger difference at 21dpS (Figs. 3D and S1A, Table S3A). Notably, the transporter slc31a1/ctr1 involved in dietary copper uptake was also significantly downregulated at both timepoints. Thus, while shorter durations of starvation such as 3 days lead to a downregulation of some metabolic functions, most of the genes involved in lipid metabolism are not significantly downregulated until 21 days of starvation. Prolonged starvation therefore leads to reduced expression of genes involved in lipid biosynthesis and transport, perhaps representing an adaptation to the prolonged absence of dietary fats and other nutrients. However, this contrasts with shorter periods of starvation, such as 48 h, where other zebrafish studies have observed an increase in lipid catabolism, potentially to increase available energy and improve resistance to cold .
The genes significantly downregulated by starvation were also enriched for host immune functions. For example, the signal transducer stat1b, which is required for inflammatory responses in the intestine and for myeloid development in zebrafish [69, 70], and the interferon responsive gene ifit8 are downregulated by 3dpS and continuing through 21dpS. By 21dpS, ifit15 and the antiviral protein rsad2 are also significantly downregulated. In accord, downregulated genes at 21dpS were enriched for functions involved in “defense response to virus”. Finally, the carboxypeptidase cpa5 which is considered a marker for mast cells in zebrafish  was also significantly downregulated at 21dpS, suggesting a potential reduction in mast cell number or activity in the intestinal tissue. Although several immune-related genes were downregulated in starved fish, complement proteins c6 and c7 were both significantly upregulated in starved fish (Table S3). Our analysis of genes downregulated during starvation therefore suggests a reduction or impairment in immune function and inflammatory tone in the intestine during starvation, along with significant reductions in lipid metabolism and lipoprotein production. Reduced immune function during starvation may represent a mechanism contributing to the microbial community alterations observed at those timepoints.
While there were too few significant genes after refeeding at 3dpR to permit analysis of functional enrichment, several of these genes were suggestive of potential intestinal functions. This included increased expression at 3dpR of genes encoding the tandem-duplicated trypsin-like serine proteases prss59.1 and prss59.2. This small set of genes also included three mitochondrial enzymes beta carotene dioxygenase-like gene bco2l, involved in cleavage of dietary carotenoids into retinoids towards Vitamin A synthesis; and dimethylglycine dehydrogenase dmgdh, involved in glycine synthesis and production of sarcosine in the choline oxidation pathway. Notably, Dmgdh was previously shown to be induced in mouse livers upon fasting, and reduced in the livers of ground squirrels preparing for hibernation [72, 73]. Of the 11 genes differentially expressed in 3dpR refed fish compared to fed controls, 6 were also differentially expressed at 21dpS including prss59.1, prss59.2, and bco2l. These may represent starvation adaptations that remain altered after restoration of feeding or that respond to both starvation as well as to refeeding.
Although we had already removed from this analysis any genes that were differentially expressed between fed control timepoints (Table S3B, Fig. S2), we wanted to further evaluate whether there were broader biological processes that may have been differential between those fed control samples that could affect our comparisons with starved/refed animals. We therefore performed GO term analysis of genes identified as significantly different between our fed control timepoints. The GO term “lipid metabolic process” was significantly enriched among genes that were significantly downregulated in 21dpS fed relative to 3dpS fed fish. Conversely, the GO terms “lipid localization” and “response to lipid” were significantly enriched among genes that were significantly upregulated in 21dpS fed relative to 3dpS fed fish (Table S3D). Importantly, the GO term “lipid metabolic process” was also enriched in genes that were significantly downregulated in 21dpS starved relative to 21dpS fed fish, even after genes that were significant in our control analysis were removed (Fig. S1E). This raised the possibility that our observed impacts of starvation on lipid metabolism genes here may be driven in part by unusually low expression of certain lipid metabolic genes in 21dpS fed fish, whereas other related lipid metabolic functions may be unusually high in 21dpS fed fish relative to the other fed timepoints. We therefore evaluated the log2 fold changes of genes from this control analysis alongside genes that were significantly different between starved and fed fish to discern if some of these differences may be driven by the control 21dpS fed fish (Fig. S3). We found that genes under the GO terms “ribosome” and “ribosome large subunit biogenesis” do not have differential expression in starved fish that is affected by unusual gene expression in the fed fish (Fig. S3C-D). In contrast, a subset of genes such as pdk3b, syt1b, apoa1a, apoa4b.1, and fads2 which are significantly downregulated in starved fish at 21dpS relative to 21dpS fed, may be due in part to unusually high expression of these genes in 21dpS fed fish (Fig. S3A-B). However, most genes emphasized here, such as elovl2, pla2g12b, slc31a1, and many others, are not affected by abnormalities within the fed fish cohort and are likely true biological effects of the starvation treatment. To confirm this, we generated an independent cohort of adult zebrafish subjected to 21 days of starvation or fed normally and used quantitative RT-PCR to evaluate intestinal expression of lipid metabolism genes. In accord with our RNA-seq data from the original experimental cohort (Table S3), starved fish at 21dpS displayed significantly reduced expression of elovl2, pla2g12b, and apoa4b.2 compared to 21dpS fed controls (Fig. S4).
Several genes responsive to starvation are also responsive to high fat feeding
To interpret which starvation-responsive genes from our dataset responded transcriptionally across a broad range of nutrient availability, and which ones may constitute a starvation-specific response, we referenced our intestinal RNA-seq results against previously published RNA-seq data comparing digestive tracts from zebrafish larvae that were either unfed or fed a high-fat meal (chicken egg yolk) . This revealed a large overlap in significantly differentially-expressed genes (Fig. 4A). Particularly, genes involved in lipid transport and metabolism such as fabp1b.1 and pla2g12b that were downregulated during starvation were upregulated during high fat feeding in zebrafish, underscoring the ability of these genes to respond to nutrients in zebrafish. Several genes involved in immune function such as rsad2, stat1b, and ifit15 were downregulated during starvation, and were upregulated after high fat feeding. Also among the overlapping genes was the enteropeptidase tmprss15 which was upregulated during starvation but downregulated by high fat feeding, and gamma butyrobetaine hydroxylase bbox1 which which was upregulated by both starvation and high fat feeding.
While there was an overlap between genes in the above datasets that implicated them in the intestinal nutrient response, several genes that were significantly affected by starvation were not significantly affected by high fat feeding. These genes included the complement factor c6, the fatty acid elongase elovl2, and the phospholipase pla2g4c. These findings suggest that some classes of genes involved in lipid transport or inflammation may be differentially regulated by factors uniquely associated with starvation and not nutrient excess inherent to high fat feeding. Alternatively, these differences could be ascribed to transcriptional responses unique to zebrafish life stages (adult vs larvae) or organs (intestine vs complete digestive tract including intestine, liver, pancreas, and swim bladder), or to indirect effects of high fat egg yolk feeding that are unrelated to nutrition.
A small subset of genes responsive to starvation are also responsive to microbial colonization
We and others have shown that intestinal gene expression is regulated in part by the presence and composition of the intestinal microbiome [55, 74,75,76]. Our 16S rRNA sequence data revealed that starvation induced marked and reversible alterations to gut microbiome composition including enrichment of Vibrio sp., members of which have been shown to be pro-inflammatory, and a decrease in similarity in microbiome composition between starved fish and their environmental samples . Although these results suggest altered gut microbial ecology during starvation, our study design did not permit us to causally link our observed changes in intestinal transcriptome and microbiome. Therefore, in order to identify transcriptional responses to starvation that may also be sensitive to microbiome, we compared our RNA-seq data to a previous study investigating the effect of microbial colonization on larval zebrafish digestive tracts . We found no correlation between the two datasets, implying that there may not be extensive overlaps between transcriptional responses to microbial colonization and starvation at these timepoints (Fig. 4B). This modest overlap may be due to transcriptional responses unique to zebrafish life stages (adult vs larvae) or organs (intestine vs complete digestive tract including intestine, liver, pancreas, and swim bladder). However, we did identify several overlapping genes that were significant in both datasets. Two complement factors, c6 and c7, which were upregulated in colonized zebrafish, were also upregulated during starvation. In addition, two genes involved in the processing of major histocompatibility complex, mad2l and erap1, were also among the overlapping genes. A small set of genes involved in lipid metabolism and intracellular cholesterol transport, such as elovl2, fads2, fabp1b.1 and scp2a were significantly downregulated by both colonization and starvation. The enteropeptidase tmprss15, which was differentially expressed during both starvation and high fat feeding, was also significantly upregulated by microbiome colonization. Overall, this comparison identified candidate genes that respond to both starvation and microbiome colonization.
Starvation-responsive genes may also be controlled by the transcription factor hepatocyte nuclear factor 4 alpha (HNF4A)
We next sought to identify transcription factors putatively linked to the regulation of the starvation response. Using HOMER , we queried the genomic regions near all genes identified as significantly upregulated or downregulated by starvation. We restricted our search to regions within the gene body plus the flanking 10 kb upstream and downstream that we previously identified as accessible chromatin in the zebrafish intestine . This revealed vertebrate transcription factor motifs significantly enriched near starvation responsive genes at either 3dpS or 21dpS (Fig. 5A; Table S3C). PAX5 motifs were significantly enriched near downregulated genes, while FOXA1 motifs were significantly enriched near upregulated genes, both at 3dpS and 21dpS. Both PAX5 and FOXA1 have been implicated in intestinal development in mice [78, 79]. Further, HNF4A motifs were significantly enriched near downregulated genes at both 3dpS and 21dpS. We previously showed that the nuclear receptor HNF4A mediates host transcriptional responses to microbial colonization in zebrafish . In addition, HNF4A in other animals is required for intestine-specific gene regulation and has conserved roles in glucose homeostasis, gluconeogenesis, and lipid metabolism, indicating that starvation-linked genes may be under the control of HNF4A [80,81,82,83]. This observation suggested that HNF4 transcription factors might facilitate both responses to starvation and changes to host microbiome.
Based on these previous findings, a comparison of our starvation dataset to an RNA-Seq dataset from hnf4a mutant zebrafish  demonstrated that genes putatively controlled by hnf4a were significantly downregulated in our dataset (Fig. 5B). GO Term analysis revealed that “lipid metabolic process”, “viral response”, and a variety of other metabolic functions were significantly enriched among these overlapping genes. Specifically, the genes pla2g12b and elovl2, and other genes involved in lipid metabolism that were downregulated significantly in starved fish, were all downregulated in hnf4a−/− fish, implying that lipid metabolic responses to starvation might be positively regulated by Hnf4a. Similarly, several of the immune response genes differentially regulated in starved fish such as ifit15, ifit8, c6, and erap1a were also differentially regulated in hnf4a −/− fish, suggesting that the immune response to starvation may also be partly influenced by Hnf4a function.
Interaction between the microbiome and host metabolism is known to occur in diverse pathophysiological contexts including starvation and malnutrition. However, few previous studies have simultaneously explored changes in host gene expression and microbiome composition as a function of starvation . We focused here on the intestine as the animal’s primary interface with the gut microbiome and dietary nutrients. Our RNA-Seq data suggests that cytoplasmic translation, ribosomal genes, and ribosomal synthesis genes are upregulated in the zebrafish intestine early in starvation, whereas DNA repair, and vitamin and cofactor metabolism genes become upregulated at 21dpS. Similarly, some pathways significantly downregulated at 21dpS were distinct from those downregulated at 3dpS, with 21dpS including genes involved in antiviral response, arginine and proline metabolism, and glycerophospholipid metabolism, among others. The distinct functions encoded at 3dpS and 21dpS suggest different stages of starvation, as previously described in zebrafish liver . In sharp contrast to previous studies in other organs, only a handful of genes were differentially expressed after refeeding in the zebrafish intestine, suggesting that the adaptive physiology displayed by the intestine during prolonged starvation is rapidly reversible after refeeding. In starved and refed zebrafish livers, upregulated genes are enriched for functions such as the TCA cycle, proteasome assembly, oxidative phosphorylation, and DNA replication and repair . Similar compensatory mechanisms have been observed to accompany refeeding in cattle livers, as well as in salmon and trout muscle [85,86,87]. These results suggest that the intestine may be particularly plastic in its adaptation to starvation and refeeding compared to other organs such as muscle and liver.
To explore gene regulatory mechanisms underlying the intestinal response to starvation, we provide evidence that hnf4a may regulate a substantial number of these starvation-associated changes, expanding the already large number of physiologic functions associated with this gene. Considering that hnf4a activity is suppressed by the microbiome in zebrafish and mice , hnf4a may link alterations in the host microbiome and transcriptome during starvation. Future studies could test the impact of starvation on Hnf4a occupancy using ChIP-Seq, or on chromatin accessibility or histone modifications in the intestinal epithelium to identify cis-regulatory regions involved in coordinating the starvation response. Our data also provide numerous candidate genes that can be used in future experiments to explore the specificity, regionality, and regulation of the starvation responses in the zebrafish intestine.
Although our RNA-Seq data suggests many commonalities between the starvation response in zebrafish and other vertebrates, it also highlights unique ways in which the zebrafish intestine may adapt to long-term starvation. For example, we observed an induction of complement proteins, ribosomal proteins, and a downregulation of the antiviral response during starvation. This is in contrast to rainbow trout liver where starvation was reported to reduce expression of ribosomal proteins . Meanwhile, genes significantly downregulated at 21dpS included several involved in the antiviral response. These pathways have not been previously reported in other animals in the context of starvation and thus could represent adaptive mechanisms unique to the zebrafish.
Starved zebrafish exhibited significantly reduced growth that was not fully recovered during a 21-day refeeding timeline. This suggests that full somatic recovery from starvation may require more time, or that there are permanent somatic changes associated with starvation. In contrast, we find that the changes that starvation induces in the zebrafish intestinal transcriptome and microbiome are rapidly normalized after refeeding. Whereas starvation significantly affected the expression of over 200 unique genes in the intestine compared to fed controls, refeeding for just 3 days restored normal levels of expression for all but 9 genes (Fig. 3B-C, Fig. S2D). Similarly, intestinal microbial communities subjected to starvation displayed significantly increased diversity (Fig. 1D) and altered composition (Fig. 2A) compared to fed controls, yet those differences were largely normalized within 1 to 3 days of refeeding. By comparison, another animal that undergoes prolonged starvation, the hibernating ground squirrel, maintained baseline levels of intestinal microbiome diversity during early stages of winter hibernation, reduced diversity later in the winter, and then increased diversity upon refeeding in the spring . That boost in diversity upon refeeding was attributed to new bacterial taxa associated with the introduced food. For this study, no samples of food-associated bacterial taxa were taken, so we cannot distinguish between these two possible explanations. Regardless, the distinct effects of starvation and refeeding on intestinal microbiome diversity in zebrafish and ground squirrels underscores the importance of studying the ecology and physiology of prolonged starvation and refeeding in diverse animal hosts.
Our results also provide insight into the specific bacterial lineages that are most sensitive to starvation and refeeding in the zebrafish intestine. We previously demonstrated that Vibrio and Plesiomonas genera are part of a core gut microbiome of zebrafish . We speculate that the opposing changes in relative abundance of these two taxa likely reflect differing abilities to survive in the altered environment of the starved gut (Fig. 2B, C, E). Vibrio sp. are common members of the intestinal microbiome in zebrafish, and their relative abundance correlates positively with intestinal inflammation [41, 42, 48]. This suggests that starvation-induced alterations in the relative abundance of Vibrio sp. and other bacteria in the zebrafish intestine might be linked with alterations in intestinal gene expression. During starvation in chickens, intestinal mucus is known to increase in abundance and thickness, possibly creating a competitive advantage for mucin-degrading bacteria [90, 91]. Vibrio spp. can degrade intestinal mucus, which may be why there is an observed increase in Vibrio during starvation [92, 93]. Conversely, Plesiomonas may be less suited for survival during prolonged starvation periods within the gut. It remains unclear if these changes in relative abundance were accompanied by alterations in microbial community density, which could be explored in future studies. It is striking that these and other starvation-induced perturbations to gut microbiome composition, similar to host gene expression in the gut, were largely restored within 1 to 3 days after refeeding. This underscores remarkable plasticity in intestinal physiology and microbial ecology in response to starvation and refeeding.
Unless otherwise stated, all fish were maintained on a 14-h light cycle at 28 °C. All zebrafish used for these experiments were born on the same day from 1 (RNA-Seq) or 3 (16S rRNA gene amplicon sequencing and qRT-PCR assays) breeding pairs from a single sibship.
Animals used for RNA-Seq and 16S rRNA gene amplicon sequencing were generated at Duke University as follows. Fertilized embryos were transferred into Petri dishes containing egg water (6 g sea salt, 1.5 g calcium sulfate, 0.75 g sodium bicarbonate, 10–12 drops methylene blue, 10 L water) at a density of 50 embryos/dish incubated at 28.5 °C. At 1-day post-fertilization (dpf), embryos were transferred to 3 L tanks containing 500 mL water from a recirculating zebrafish aquaculture system (system water). Each tank contained 10 (16S rRNA gene sequencing) or 30 (RNA-Seq) embryos. Fish were then maintained under standard husbandry on a recirculating zebrafish aquaculture system until the start of the experiment at 60dpf. Zebrafish were then randomly transferred into four (RNA-Seq) or eight (16S rRNA gene sequencing) clean 10 L tanks at a density of 44 (RNA-Seq) or 67 (16S rRNA gene sequencing) fish per tank, with half the tanks receiving no food for the following 21 days (Fig. 1A). Following the 21 days of starvation, feedings for all tanks were allowed to occur as per standard husbandry: two feedings of live Artemia per day interspersed with two feedings of Gemma 300 (Skretting). Over the 21 days of starvation and 21 days of refeeding, all fish remained on the recirculating zebrafish aquaculture system and we observed no mortality in any condition or experiment.
Animals used for qRT-PCR validation were generated at the Carnegie Institution as follows. Fertilized embyros were transferred at 5dpf to 10 L tanks, with each tank containing approximately 70 larvae. Fish were fed a regimen of Gemma Micro (Skretting) with gradually increasing pellet sizes. For the first 2 weeks of feeding, larvae were fed exclusively Gemma Micro 75, and then transitioned to a diet of Gemma Micro 150 supplemented with live Artemia (www.artemia-international.com) for the subsequent 3 weeks, at which point fish were transitioned to anadult diet of Gemma Micro 300 with only occasional live brine supplementation. Adults raised in this manner were transferred to one of 12 experimental 3 L tanks that were balanced for density and gender ratio, and half of those tanks received no food for the following 21 days.
All fish to be sampled on a particular day were collected prior to the first daily feeding in the fish facility. Samples for 16S rRNA gene amplicon sequencing were collected at 0 days post-starvation (0dpS), 1 dpS, 3dpS, 7dpS, 21dpS, 1 day post-re-feed (dpR), 3dpR, 7dpR, and 21dpR (Fig. 1A) with six randomly selected fish at each time point per tank were euthanized by tricaine overdose (0.83 mg/ml tricaine). Fish were imaged on a dissecting scope to facilitate subsequent standard length (SL) and height at anterior of anal fin (HAA) measurements . Intestinal tracts were then dissected from each fish and placed individually in lysis buffer (20 mM Tris-HCl (pH 8.0), 2 mM EDTA (pH 8.0), 1% Triton X-100, flash-frozen in a dry-ice/ethanol bath, and stored at − 80 °C until DNA extraction.
Samples for RNA-Seq were collected at 3dpS, 21dpS, and 3dpR. At each time point, three randomly selected fish per tank were euthanized by tricaine overdose (0.83 mg/ml tricaine). Fish were imaged on a dissecting scope to facilitate subsequent standard length (SL) . Intestinal tracts were then dissected from each fish and placed individually in 2 mL cryovials filled with TRIzol reagent (Thermo Fisher, 15,596,026), flash-frozen in a dry ice-ethanol bath, and stored at − 80 °C until RNA extraction.
Samples for qRT-PCR were collected at 21dpS. Fish were euthanized by tricaine overdose, and their intestines were dissected and stored individually in 0.5 mL of RNAlater.
16S rRNA gene sequencing
Genomic DNA was extracted from individual zebrafish intestinal tracts using Qiagen DNeasy Blood and Tissue Kits (Qiagen, modified with bead-beating as previously described) . Genomic DNA was subsequently used as template for PCR amplification of the v4 region of 16S rRNA gene and 150 bp/cycle paired-end sequencing was performed on an Illumina HiSeq 2000 Sequencing System (see Table S4 for primers) at the University of Oregon Genomics and Cell Characterization Core Facility.
16S rRNA gene sequence bioinformatic and statistical analysis
FASTQ files were demultiplexed and split by sample ID using QIIME (v1.9.1). Within RStudio version 3.4.1, the files were then quality filtered, trimmed, denoised, merged, checked for chimeras, and assigned taxonomy using DADA2. Taxonomic assignments were made using the Silva v132 database. Data analysis used the R packages vegan and phyloseq. LEFsE was accessed through the Huttenhower Galaxy website: https://huttenhower.sph.harvard.edu/galaxy/. The following packages were used for analyses and the creation of figures: ape, picante, vegan, ggtree, SEPP, plotyly, heatmaply, and pairwiseadonis [95,96,97,98,99,100,101,102].
RNA extraction and sequencing
Frozen whole intestinal samples stored at − 80 °C were homogenized using Zirconium oxide beads (Biospec, 11,079,107) and a Vortex Genie2 (Scientific Industries, 1311-V) fitted with a Vortex Adapter (Scientific Industries, 13,000-V1–24) in three 45-s intervals. Samples were put on ice in-between homogenization to prevent overheating. Following homogenization, a phase separation was performed by adding 200ul of chloroform to each sample and mixing by vigorous inversion 15 times. Samples were then incubated at room temperature for 3 min and centrifuged at 12000rcf for 15 min at 4 °C. 500ul of the aqueous upper phase from each sample was then transferred to a new Eppendorf tube, to which 500uL 70% Ethanol in DEPC water was added and vortexed. Following phase separation, samples were DNase treated and total RNA was extracted via column purification using the PureLink DNase Set (Thermo Fisher, 12,185,010) and the PureLink RNA Mini kit (Thermo Fisher, 12,183,025) according to the manufacturer’s instructions. Final sample quality and concentration were assessed via spectrophotometry and samples were stored at -80 °C until submission to the Duke Sequencing and Genomic Technologies Core. RNA-seq libraries were prepared and sequenced by Duke Sequencing and Genomic Technologies Core on an Illumina HiSeq 2500 to generate 50 bp single-end reads (SR).
All raw zebrafish RNA-seq data was processed on the Galaxy server . Raw fastq files were trimmed using Trim Galore . Trimmed fastq files were then mapped to the zebrafish genome (GRCz10) using STAR using default settings to generate BAM files, which were converted to counts using HTSeq (v.0.9.1). BAM files were converted to bigWig files using the bamCompare tool (v2.5.0) with RPKM normalization before visualization on the UCSC Genome Browser [105,106,107].
TPM expression values were obtained for transcripts via Salmon . Pairwise differential gene expression tests were carried out with DESeq2 using counts files generated by HTSeq [108, 109]. For comparisons between starved and fed fish, the default significance threshold of adjusted p-value 0.05 was used for each comparison. For comparisons across fed fish controls, the significance threshold was defined as the gene either having an absolute log2 fold change greater than 1.0 or a p-value less than 0.05 (See Figs. S1, S3, and Table S3B). Versions of the analysis that include the genes removed by this filtering step are available in Fig. S2.
Hierarchical clustering of log2 fold change values for genes was performed using Cluster 3.0, and heat maps were generated using Java Treeview [110, 111].
HOMER software (http://homer.ucsd.edu/homer/motif/) analysis was performed on significantly upregulated and downregulated genes at both 3dpS and 21dpS (3dpS starved/3dpS fed and 21dpS starved/21dpS fed, respectively), using regions within the gene body plus the flanking 10 kb upstream and downstream that we previously identified as accessible chromatin in the zebrafish intestine  using the findMotifs.pl command. A motif score was obtained by taking the -log10 values of the p-values assigned by HOMER. A motif was then deemed ‘enriched’ amongst either upregulated or downregulated genes at each timepoint (3dpS or 21dpS) based on whether it had a higher motif score among the upregulated or downregulated gene sets.
For comparisons with the larval zebrafish egg yolk feeding dataset, log2 fold changes in 21dpS fed fish relative to 21dpS starved were compared to log2 fold changes in larval zebrafish digestive tracts 4 h after egg yolk feeding (i.e. “HF 4h logFC”) obtained from Supplementary Table 1 in , the raw data for which is available at accession GSE87704. For comparisons with hnf4a mutant and microbially colonization datasets, data was obtained from Supplemental Table 2 in , using log2 fold changes comparing digestive tracts from hnf4a homozygous mutant and wild-type 6dpf zebrafish larvae raised under conventionalized ex-germ-free conditions (“MutCV/WTCV”) and from wild-type 6dpf zebrafish larvae reared under germ-free of ex-germ-free conventionalized conditions (“WTGF/WTCV”), respectively. Raw data from  is available at accession GSE90462.
Quantitative RT-PCR assays
Dissected intestines from 16 starved and 15 fed 5–6-month-old fish were stored in RNAlater and shipped to the Rawls lab. Tissue was homogenized with bead beater tubes (Bertin Corp, P000912-LYSK0-A) and RNA was extracted with chloroform, precipitated with isopropanol, washed with ethanol, and treated with DNase (Invitrogen, AM1906). RNA concentration was assessed using Qubit RNA BR Assay Kit (ThermoFisher, Q10210). 800 ng of RNA were then reverse transcribed using the iScript cDNA Synthesis kit (Bio-Rad, 1708891). Quantitative real-time PCR was performed on QuantBio Studio6 (Applied Biosystems) using gene-specific primers for pla2g12b (F 5′-CAGTACCGCTGCAGATATGGT-3′, R 5′- ATTCGGTACCTGGAAGCCAAG-3′), elovl2 (F 5′- ACAGTTTTCAGCTGTCCCGT-3′, R 5′- CATCCTCTCACGCGGGTATC-3′), apoa4b.2 (F 5′-TTGTGGTCTTTGCACTTGCT-3′, R 5′- TCATCTTGACGGTTTCCTCTG-3′) and ef1a (F 5′- CTTCTCAGGCTGACTGTGC-3′, R 5′- CCGCTAGGATTACCCTCC-3′). 6 samples (3 male, 3 female) in each group were randomly chosen and analyzed in duplicate 25 μl reactions using 2X SYBR Green SuperMix (PerfeCTa, Hi Rox, Quanta Biosciences, 95055) and normalized to the expression of ef1a as a house-keeping gene. Expression profile and associated statistical parameters were determined using the ∆∆CT method and graphed in Prism (GraphPad).
Availability of data and materials
All quality filtering parameters for generating the sequence variants, ASV table and figures for the 16S rRNA analysis can be found at: https://github.com/alexmccumber/fishguts. The raw 16S rRNA gene amplicon FASTQ files can be accessed from the European Nucleotide Archive under project access number PRJEB31503. Raw and processed RNA-Seq data is available on NCBI GEO at the accession GSE140821.
Watts, M. J. & Bohle, H. G. The space of vulnerability: the causal structure of hunger and famine: Progress in Human Geography (2016) doi:https://doi.org/10.1177/030913259301700103.
Fearon K, et al. Definition and classification of cancer cachexia: an international consensus. Lancet Oncol. 2011;12:489–95.
Nagy K. Food requirements of wild animals: predictive equations for free-living mammals, reptiles, and birds; 2018.
McCue MD. Starvation physiology: reviewing the different strategies animals use to survive a common challenge. Comp Biochem Physiol A Mol Integr Physiol. 2010;156:1–18.
Storey KB, Storey JM. Metabolic Rate Depression and Biochemical Adaptation in Anaerobiosis. Hibernation and Estivation Quarterly Review Biol. 1990;65:145–74.
Nagy KA, Girard IA, Brown TK. energetics of free-ranging mammals, reptiles, and birds. Ann Review Nutri Palo Alto. 1999;19:247–77.
Thomas DR. Loss of skeletal muscle mass in aging: examining the relationship of starvation, sarcopenia and cachexia. Clin Nutr. 2007;26:389–99.
Love AH. Metabolic response to malnutrition: its relevance to enteral feeding. Gut. 1986;27:9–13.
Moore FD. Energy and the maintenance of the body cell mass. J Parenter Enter Nutr. 1980;4:228–60.
McFarlane H, et al. Immunity, transferrin, and survival in kwashiorkor. Br Med J. 1970;4:268–70.
Champakam S, Srikantia SG, Gopalan C. Kwashiorkor and mental development. Am J Clin Nutr. 1968;21:844–52.
Acheson RM, Macintyre MN. The effects of acute infection and acute starvation on skeletal development. Br J Exp Pathol. 1958;39:37–45.
Palesty JA, Dudrick SJ. The goldilocks paradigm of starvation and Refeeding. Nutr Clin Pract. 2006;21:147–54.
Cahill GF. Starvation in man. N Engl J Med. 1970;282:668–75.
Felig P, Owen OE, Wahren J, Cahill GF. Amino acid metabolism during prolonged starvation. J Clin Invest. 1969;48:584–94.
Désert C, et al. Transcriptome profiling of the feeding-to-fasting transition in chicken liver. BMC Genomics. 2008;9:611.
Drew RE, et al. Effect of starvation on transcriptomes of brain and liver in adult female zebrafish (Danio rerio). Physiol Genomics. 2008;35:283–95.
Hakvoort, T. B. M. et al. Interorgan coordination of the murine adaptive response to fasting. J. Biol. Chem. jbc.M110.216986 (2011) doi:https://doi.org/10.1074/jbc.M110.216986.
Crawford PA, et al. Regulation of myocardial ketone body metabolism by the gut microbiota during nutrient deprivation. PNAS. 2009;106:11276–81.
David LA, et al. Diet rapidly and reproducibly alters the human gut microbiome. Nature. 2014;505:559–63.
Hildebrandt, M. A. et al. High-fat diet determines the composition of the murine gut microbiome independently of obesity. Gastroenterology 137, 1716–1724.e1–2 (2009).
Turnbaugh PJ, Bäckhed F, Fulton L, Gordon JI. Diet-induced obesity is linked to marked but reversible alterations in the mouse distal gut microbiome. Cell Host Microbe. 2008;3:213–23.
Basolo, A. et al. Effects of underfeeding and oral vancomycin on gut microbiome and nutrient absorption in humans. Nature Medicine 1–10 (2020) doi:https://doi.org/10.1038/s41591-020-0801-z.
Backhed, F. et al. gut microbiota as an environmental factor that regulates fat storage. Proceedings of the National Academy of Sciences of the United States of America 2004.
Smith MI, et al. Gut microbiomes of Malawian twin pairs discordant for kwashiorkor. Science. 2013;339:548–54.
Ridaura VK, et al. Gut microbiota from twins discordant for obesity modulate metabolism in mice. Science. 2013;341:1241214.
Cowardin CA, et al. Mechanisms by which sialylated milk oligosaccharides impact bone biology in a gnotobiotic mouse model of infant undernutrition. PNAS. 2019;116:11988–96.
Hibberd, M. C. et al. The effects of micronutrient deficiencies on bacterial species from the human gut microbiota. Sci Transl Med 2017 9:eaal4069.
Gehrig, J. L. et al. Effects of microbiota-directed foods in gnotobiotic animals and undernourished children. Science 2019 365, eaau4732.
Leulier F, et al. Integrative physiology: at the crossroads of nutrition, microbiota, animal physiology, and human health. Cell Metab. 2017;25:522–34.
Gas N, Noailliac-Depeyre J. Studies on intestinal epithelium involution during prolonged fasting. J Ultrastruct Res. 1976;56:137–51.
Segner, H. & Braunbeck, T. Hepatocellular adaptation to extreme nutritional conditions in ide,Leuciscus idus melanotus L. (Cyprinidae). A morphofunctional analysis. Fish Physiol Biochem 1988 5, 79–97.
Lenaerts K, et al. Starvation induces phase-specific changes in the proteome of mouse small intestine. J Proteome Res. 2006;5:2113–22.
Saotome I, Curto M, McClatchey AI. Ezrin is essential for epithelial organization and villus morphogenesis in the developing intestine. Dev Cell. 2004;6:855–64.
Habold C, Foltzer-Jourdainne C, Le Maho Y, Lignot J-H. Intestinal apoptotic changes linked to metabolic status in fasted and refed rats. Pflugers Arch - Eur J Physiol. 2006;451:749–59.
Goldsmith MI, Iovine MK, O’Reilly-Pol T, Johnson SL. A developmental transition in growth control during zebrafish caudal fin development. Dev Biol. 2006;296:450–7.
Flynn EJ, Trent CM, Rawls JF. Ontogeny and nutritional control of adipogenesis in zebrafish (Danio rerio). J Lipid Res. 2009;50:1641–52.
McMenamin SK, Minchin JEN, Gordon TN, Rawls JF, Parichy DM. Dwarfism and increased adiposity in the gh1 mutant zebrafish vizzini. Endocrinology. 2013;154:1476–87.
Minchin JEN, Rawls JF. A classification system for zebrafish adipose tissues. Dis Model Mech. 2017;10:797–809.
Ng AN et al. Formation of the digestive system in zebrafish: III. Intestinal epithelium morphogenesis. Developmental Biol 2005 286:114–135.
Rawls JF, Samuel BS, Gordon JI. Gnotobiotic zebrafish reveal evolutionarily conserved responses to the gut microbiota. PNAS. 2004;101:4596–601.
Stephens WZ, et al. The composition of the zebrafish intestinal microbial community varies across development. ISME J. 2016;10:644–54.
Semova I, et al. Microbiota regulate intestinal absorption and metabolism of fatty acids in the Zebrafish. Cell Host Microbe. 2012;12:277–88.
Wong S. et al. Ontogenetic Differences in Dietary Fat Influence Microbiota Assembly in the Zebrafish Gut. mBio 2015;6:e00687–e00615.
Wallace KN, Pack M. Unique and conserved aspects of gut development in zebrafish. Dev Biol. 2003;255:12–29.
Cheesman SE, Neal JT, Mittge E, Seredick BM, Guillemin K. Epithelial cell proliferation in the developing zebrafish intestine is regulated by the Wnt pathway and microbial signaling via Myd88. Proc Natl Acad Sci U S A. 2011;108(Suppl 1):4570–7.
Bates JM, et al. Distinct signals from the microbiota promote different aspects of zebrafish gut differentiation. Dev Biol. 2006;297:374–86.
Rolig AS, Parthasarathy R, Burns AR, Bohannan BJM, Guillemin K. Individual members of the microbiota disproportionately modulate host innate immune responses. Cell Host Microbe. 2015;18:613–20.
Murdoch CC, et al. Intestinal serum amyloid a suppresses systemic neutrophil activation and bactericidal activity in response to microbiota colonization. PLoS Pathog. 2019;15:e1007381.
Kanther M, et al. Microbial colonization induces dynamic temporal and spatial patterns of NF-κB activation in the Zebrafish digestive tract. Gastroenterology. 2011;141:197–207.
Lickwar CR, et al. Genomic dissection of conserved transcriptional regulation in intestinal epithelial cells. PLoS Biol. 2017;15:e2002054.
Jia J, et al. Microarray and metabolome analysis of hepatic response to fasting and subsequent refeeding in zebrafish (Danio rerio). BMC Genomics. 2019;20:919.
Minchin JEN, Scahill CM, Staudt N, Busch-Nentwich EM, Rawls JF. Deep phenotyping in zebrafish reveals genetic and diet-induced adiposity changes that may inform disease risk. J Lipid Res. 2018;59:1536–45.
Burns AR, et al. Contribution of neutral processes to the assembly of gut microbial communities in the zebrafish over host development. ISME J. 2016;10:655–64.
Rawls JF, Mahowald MA, Ley RE, Gordon JI. Reciprocal gut microbiota transplants from Zebrafish and mice to germ-free recipients reveal host habitat selection. Cell. 2006;127:423–33.
Ar, B. et al. Interhost dispersal alters microbiome assembly and can overwhelm host innate immunity in an experimental zebrafish model. Proceedings of the National Academy of Sciences of the United States of America vol. 114 http://pubmed.ncbi.nlm.nih.gov/28973938/ (2017).
Segata N, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12:R60.
Morais S, Monroig O, Zheng X, Leaver MJ, Tocher DR. Highly unsaturated fatty acid synthesis in Atlantic Salmon: characterization of ELOVL5- and ELOVL2-like Elongases. Mar Biotechnol. 2009;11:627–39.
Zeituni, E. M. et al. Endoplasmic reticulum lipid flux influences enterocyte nuclear morphology and lipid-dependent transcriptional responses. J. Biol. Chem. jbc.M116.749358 (2016) doi:https://doi.org/10.1074/jbc.M116.749358.
Wang Y, et al. Tissue-specific, nutritional, and developmental regulation of rat fatty acid elongases. J Lipid Res. 2005;46:706–15.
Cruciani-Guglielmacci C, et al. Molecular phenotyping of multiple mouse strains under metabolic challenge uncovers a role for Elovl2 in glucose-induced insulin secretion. Molecular Metabolism. 2017;6:340–51.
Alimuddin, Kiron, V., Satoh, S., Takeuchi, T. & Yoshizaki, G. Cloning and over-expression of a masu salmon (Oncorhynchus masou) fatty acid elongase-like gene in zebrafish. Aquaculture 2008; 282, 13–18.
Antalis TM, Shea-Donohue T, Vogel SN, Sears C, Fasano A. Mechanisms of disease: protease functions in intestinal mucosal pathobiology. Nat Clin Pract Gastroenterol Hepatol. 2007;4:393–402.
Mukherjee S, Hooper LV. Antimicrobial defense of the intestine. Immunity. 2015;42:28–39.
Holzinger A, et al. Mutations in the Proenteropeptidase gene are the molecular cause of congenital Enteropeptidase deficiency. Am J Hum Genet. 2002;70:20–5.
Thierer JH, Ekker SC, Farber SA. The LipoGlo reporter system for sensitive and specific monitoring of atherogenic lipoproteins. Nat Commun. 2019;10:3426.
Salem M, Silverstein J, Rexroad CE, Yao J. Effect of starvation on global gene expression and proteolysis in rainbow trout (Oncorhynchus mykiss). BMC Genomics. 2007;8:328.
Lu D-L, et al. Fasting enhances cold resistance in fish through stimulating lipid catabolism and autophagy. J Physiol. 2019;597:1585–603.
Song H, Yan Y, Titus T, He X, Postlethwait JH. The role of stat1b in zebrafish hematopoiesis. Mech Dev. 2011;128:442–56.
Richmond CA, et al. JAK/STAT-1 signaling is required for reserve intestinal stem cell activation during intestinal regeneration following acute inflammation. Stem Cell Reports. 2017;10:17–26.
Dobson JT, et al. Carboxypeptidase A5 identifies a novel mast cell lineage in the zebrafish providing new insight into mast cell fate determination. Blood. 2008;112:2969–72.
Hindle AG, Grabek KR, Epperson LE, Karimpour-Fard A, Martin SL. Metabolic changes associated with the long winter fast dominate the liver proteome in 13-lined ground squirrels. Physiol Genomics. 2014;46:348–61.
Kamata S, et al. 2D DIGE proteomic analysis reveals fasting-induced protein remodeling through organ-specific transcription factor(s) in mice. FEBS Open Bio. 2018;8:1524–43.
Davison JM, et al. Microbiota regulate intestinal epithelial gene expression by suppressing the transcription factor hepatocyte nuclear factor 4 alpha. Genome Res. 2017;27:1195–206.
Camp, J. G. et al. Microbiota modulate transcription in the intestinal epithelium without remodeling the accessible chromatin landscape. Genome Res. gr.165845.113 (2014) doi:https://doi.org/10.1101/gr.165845.113.
Larsson E, et al. Analysis of gut microbial regulation of host gene expression along the length of the gut and regulation of gut microbial ecology through MyD88. Gut. 2012;61:1124–31.
Heinz S, et al. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol Cell. 2010;38:576–89.
Xu H, et al. Transcriptional inhibition of intestinal NHE8 expression by glucocorticoids involves Pax5. Am J Physiol-Gastrointestinal Liver Physiol. 2010;299:G921–7.
Kaestner KH. Making the liver what it is: the many targets of the transcriptional regulator HNF4α. Hepatology. 2010;51:376–7.
Roman AKS, Aronson BE, Krasinski SD, Shivdasani RA, Verzi MP. Transcription factors GATA4 and HNF4A control distinct aspects of intestinal homeostasis in conjunction with transcription factor CDX2. J Biol Chem. 2015;290:1850–60.
Palanker L, Tennessen JM, Lam G, Thummel CS. Drosophila HNF4 regulates lipid mobilization and β-oxidation. Cell Metab. 2009;9:228–39.
Frochot V, et al. The transcription factor HNF-4α: a key factor of the intestinal uptake of fatty acids in mouse. Am J Physiol-Gastrointestinal Liver Physiol. 2012;302:G1253–63.
Barry WE, Thummel CS. The Drosophila HNF4 nuclear receptor promotes glucose-stimulated insulin secretion and mitochondrial function in adults. Elife. 2016;5.
Xia JH, et al. The intestinal microbiome of fish under starvation. BMC Genomics. 2014;15:266.
Rescan P-Y, et al. Dynamic gene expression in fish muscle during recovery growth induced by a fasting-refeeding schedule. BMC Genomics. 2007;8:438.
Connor EE, et al. Enhanced mitochondrial complex gene function and reduced liver size may mediate improved feed efficiency of beef cattle during compensatory growth. Funct Integr Genomics. 2010;10:39–51.
Rescan P-Y, Le Cam A, Rallière C, Montfort J. Global gene expression in muscle from fasted/refed trout reveals up-regulation of genes promoting myofibre hypertrophy but not myofibre production. BMC Genomics. 2017;18:447.
Carey HV, Walters WA, Knight R. Seasonal restructuring of the ground squirrel gut microbiota over the annual hibernation cycle. Am J Physiol Regul Integr Comp Physiol. 2013;304:R33–42.
Roeselers G, et al. Evidence for a core gut microbiota in the zebrafish. ISME J. 2011;5:1595–608.
Schroeder BO. Fight them or feed them: how the intestinal mucus layer manages the gut microbiota. Gastroenterol Rep (Oxf). 2019;7:3–12.
Smirnov A, Sklan D, Uni Z. Mucin dynamics in the chick small intestine are altered by starvation. J Nutr. 2004;134:736–42.
Szabady RL, Yanta JH, Halladin DK, Schofield MJ, Welch RA. TagA is a secreted protease of Vibrio cholerae that specifically cleaves mucin glycoproteins. Microbiology. 2011;157:516–25.
Bhowmick R, et al. Intestinal adherence of Vibrio cholerae involves a coordinated interaction between colonization factor GbpA and Mucin. Infect Immun. 2008;76:4968–77.
Parichy DM, Elizondo MR, Mills MG, Gordon TN, Engeszer RE. Normal table of postembryonic zebrafish development: staging by externally visible anatomy of the living fish. Dev Dyn. 2009;238:2975–3015.
Paradis E, Claude J, Strimmer K. APE: analyses of Phylogenetics and evolution in R language. Bioinformatics. 2004;20:289–90.
Sievert, C. et al. plotly: Create Interactive Web Graphics via ‘plotly. js’. R package version 2017;4, 110.
Galili T, O’Callaghan A. Sidi, J. & Sievert, C. heatmaply: an R package for creating interactive cluster heatmaps for online publishing. Bioinformatics. 2018;34:1600–2.
Kembel SW, et al. Picante: R tools for integrating phylogenies and ecology. Bioinformatics. 2010;26:1463–4.
Oksanen J, et al. The vegan package. Community ecology package. 2007;10:631–7.
Yu G, Smith DK, Zhu H, Guan Y, Lam TT-Y. ggtree: an r package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods Ecol Evol. 2017;8:28–36.
Janssen, S. et al. Phylogenetic Placement of Exact Amplicon Sequences Improves Associations with Clinical Information. mSystems 2018;3:e00021–e00018.
McMurdie P. J. & Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PloS one 2013;8:e61217.
Afgan E, et al. The galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2018 update. Nucleic Acids Res. 2018;46:W537–44.
Krueger, F. Trim galore. A wrapper tool around Cutadapt and FastQC to consistently apply quality and adapter trimming to FastQ files 2015.
Casper J, et al. The UCSC genome browser database: 2018 update. Nucleic Acids Res. 2018;46:D762–9.
Dobin A, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21.
Ramírez F, et al. deepTools2: a next generation web server for deep-sequencing data analysis. Nucleic Acids Res. 2016;44:W160–5.
Patro R, Duggal G, Love MI, Irizarry RA, Kingsford C. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods. 2017;14:417–9.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.
Saldanha A, Java J. Treeview—extensible visualization of microarray data. Bioinformatics. 2004;20:3246–8.
de Hoon MJL, Imoto S, Nolan J, Miyano S. Open source clustering software. Bioinformatics. 2004;20:1453–4.
We are grateful to Dr. Joshua Granek who provided instruction as a part of IBIEM and that without we would have been unable to complete this work, and to Dr. Doug Turnbull for helpful advice.
This work was supported by grants to B.J.M.B., K.G., and J.F.R from the NIH (R01-GM095385) and the Gordon and Betty Moore Foundation; and to J.F.R from the NIH (R01-DK093399, R01-DK111857, and R01-DK081426). J.J. was supported by an American Heart Association Predoctoral Fellowship (825770). A.W.M. and C.A. were supported by the NSF Research Traineeship Program “Integrative Bioinformatics for Investigating and Engineering Microbiomes” or IBIEM (1545220). A.W.M. was also supported by the Center for Biomolecular and Tissue Engineering (NIH T32-GM008555). The work was also supported by funding to S.A.F. from the Carnegie Institution for Science Endowment and the G. Harold and Leila Y. Mathers Charitable Foundation.
Ethics approval and consent to participate
All zebrafish experiments were conducted in conformity with the Public Health Service Policy on Humane Care and Use of Laboratory Animals using protocols A165–13-06 and A115–16-05 approved by the Institutional Animal Care and Use Committee (IACUC) of Duke University and protocol #139 approved by IACUC of the Carnegie Institution for Science. This study was completed in compliance with ARRIVE guidelines.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1.
Additional file 2.
Additional file 3.
Additional file 4.
Additional file 5.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Jawahar, J., McCumber, A.W., Lickwar, C.R. et al. Starvation causes changes in the intestinal transcriptome and microbiome that are reversed upon refeeding. BMC Genomics 23, 225 (2022). https://doi.org/10.1186/s12864-022-08447-2