Enteroendocrine cells (EECs) are rare cells found along the intestinal mucosa and represent the most extensive endocrine system in the human body. EECs sense nutrients in the gut lumen and, in response, secrete various hormones that act locally or at a distance.
1Function and mechanisms of enteroendocrine cells and gut hormones in metabolism.
Enteroendocrine cells release numerous peptides including glucagon-like peptides 1 and 2 (GLP-1, GLP-2), peptide YY (PYY), gastric inhibitory polypeptide (GIP), secretin (SCT), cholecystokinin (CCK), ghrelin (GHRL), neurotensin (NTS), or somatostatin (SST), and amines such as serotonin (5-HT). Various dietary metabolites such as fatty acids generated by the processing of triglycerides or short-chain fatty acids produced by the fermentation of dietary fibers by the gut microbiota are detected by EECs, thereby modulating their activity. Enteroendocrine hormones emerged as essential regulators of energy homeostasis via their control of intestinal absorption, food intake and insulin secretion.
1Function and mechanisms of enteroendocrine cells and gut hormones in metabolism.
, 2Enteroendocrine cells and lipid absorption.
, 3Enteroendocrine regulation of nutrient absorption.
Importantly, GLP-1 receptor agonists are already used in the treatment of diabetes and obesity. Thus, a better understanding of the differentiation and function of the EEC system could lead to innovative strategies to manipulate them and develop novel treatments for metabolic diseases.
A constitutive lack of EECs induces lipid malabsorption. It demonstrates the critical role of enteroendocrine hormones in energy uptake and growth of neonate mice, which frequently die when born without EECs.
13- Mellitzer G.
- Beucher A.
- Lobstein V.
- et al.
Loss of enteroendocrine cells in mice alters lipid absorption and glucose homeostasis and impairs postnatal survival.
Similarly, human newborns with mutations in NEUROG3 lack all EECs in a condition termed “enteric anendocrinosis” and present with severe malabsorptive diarrhea requiring parenteral nutrition.
14- Wang J.
- Cortina G.
- Wu S.V.
- et al.
Mutant neurogenin-3 in congenital malabsorptive diarrhea.
A recent study by McCauley and colleagues shed light on the role of EEC hormones in linking nutrient sensing and absorption.
15- McCauley H.A.
- Matthis A.L.
- Enriquez J.R.
- et al.
Enteroendocrine cells couple nutrient sensing to nutrient absorption by regulating ion transport.
Despite the role of EECs in lipid absorption, the transcriptomic program induced by lipids in EECs is poorly explored. Furthermore, in these studies,
13- Mellitzer G.
- Beucher A.
- Lobstein V.
- et al.
Loss of enteroendocrine cells in mice alters lipid absorption and glucose homeostasis and impairs postnatal survival.
,14- Wang J.
- Cortina G.
- Wu S.V.
- et al.
Mutant neurogenin-3 in congenital malabsorptive diarrhea.
humans and mice with mutations in NEUROG3 have impaired EEC development since birth. Therefore, the importance of the EEC system in regulating energy homeostasis could not be addressed thoroughly in healthy adult mice because most mice with a constitutive EEC deficiency die in the first weeks of life. Consequently, potential adaptation mechanism to sustained malabsorption in adults could also not be identified. Growing evidence points towards the importance of the gut microbiota in digestion and intestinal absorption. Yet, the role of EEC-microbiota interactions in these processes and the pathogenesis of metabolic diseases are understudied. Herein, by blocking the formation of EECs in adult males, we aimed to investigate the EEC system’s function in maintaining energy homeostasis and unravel intestinal transcriptomic response to the lack of enteroendocrine hormones and in interaction with the gut microbiota.
Discussion
The role of enteroendocrine hormones in energy homeostasis has been investigated using knockouts of individual hormone genes or receptors. However, the importance of the intestinal enteroendocrine system in healthy adult mice remains unknown. Our study describes physiological and intestinal transcriptional adaptative responses to a long-term, global, EEC depletion. We found impaired lipid malabsorption and food efficiency. Yet, EECs are dispensable for survival when mice are fed with a standard CD but are required when a lipid-enriched diet is given. We reported a rapid remodeling of the intestinal microbiota upon EECs removal, steady lipid malabsorption, and adaptative metabolic and transcriptomic responses. In addition, we provided novel insights into the role of EECs in lipid metabolism by unveiling EEC transcriptional response to an HFD, which revealed an increased expression of a subset of hormones transcripts and other regulators of lipid metabolism. All experiments have been performed in males; we thus cannot conclude whether EEC-deficient female mice have a similar phenotype.
After induction of EEC deficiency, we first noticed that mice lost weight. However, when fed with an HFD, the weight dropped severely, and experimentation needed to be stopped. This observation agrees with the failure to thrive and the death of most murine neonates with intestinal EEC-deficiency when fed with high fat-containing mother’s milk.
13- Mellitzer G.
- Beucher A.
- Lobstein V.
- et al.
Loss of enteroendocrine cells in mice alters lipid absorption and glucose homeostasis and impairs postnatal survival.
Interestingly, Rfx6
ΔAdInt mice,
12- Piccand J.
- Vagne C.
- Blot F.
- et al.
Rfx6 promotes the differentiation of peptide-secreting enteroendocrine cells while repressing genetic programs controlling serotonin production.
which lack all EECs but enterochromaffin cells, only transiently lost weight. Despite lipid malabsorption, these mice caught up rapidly with wild-types, suggesting that serotonin/5-HT prevented weight loss, by promoting stored energy release.
37- Sumara G.
- Sumara O.
- Kim J.K.
- et al.
Gut-derived serotonin is a multifunctional determinant to fasting adaptation.
This might explain the more pronounced phenotype in Neurog3
ΔAdint compared with Rfx6
ΔAdInt mice. Bomb calorimetry experiments revealed decreased food efficiency in EEC-deficient mice, likely resulting from lipid malabsorption. We found that lipid absorption was spatially delayed and increased in the ileum 1 month after EECs removal. A similar distal shift of lipid absorption was observed in another mouse model of extensive EEC deficiency.
38- Terry N.A.
- Ngaba L.V.
- Wilkins B.J.
- et al.
Lipid malabsorption from altered hormonal signaling changes early gut microbial responses.
It could be a compensatory mechanism for the duodenal malabsorption resulting from EEC deficiency.
To identify the mechanisms underlying lipid malabsorption and adaptative responses of EEC-deficient mice fed a CD, we profiled the transcriptome of their duodenum and ileum and scrutinized the expression of known regulators of the different steps of lipid absorption.
21- Abumrad N.A.
- Davidson N.O.
Role of the gut in lipid homeostasis.
Although the expression of important fatty acid and cholesterol absorption regulators were unchanged in EEC-deficient mice, we found that transcripts of proteins involved in chylomicrons formation (
Apob/Apob-48 or
Ces1g), ileal bile absorption (Fabp6) or cholesterol uptake (
Abca1) were decreased in EEC-deficient mice. Also, blood concentrations of lipase, whose secretion is stimulated by CCK and SCT enteroendocrine hormones, were decreased in EEC-depleted mice, suggesting triglyceride hydrolyses could be impaired. Thus, lipid malabsorption in EEC-depleted mice might result from combined alterations of fat digestion by lipase, fatty acid and cholesterol transport, and chylomicron lipoprotein assembly.
We next focused on upregulated genes to identify possible adaptation mechanisms implemented in EEC-deficient mice to improve energy homeostasis. In the duodenum, increased expression of
Scd1,
Pcsk9, and Ldlr suggested lipogenesis and cholesterol uptake to be stimulated. Upregulated genes were numerous in the ileum. Notably,
Pla2g4c, the top gene strongly enriched in the EEC-deficient ileum, encodes the phospholipase A2 gamma, recently shown in a human hepatocellular carcinoma cell line to promote lipid droplet formation.
25- Su X.
- Liu S.
- Zhang X.
- et al.
Requirement of cytosolic phospholipase A2 gamma in lipid droplet formation.
However, such a function was so far not reported in enterocytes. Another gene highly enriched is the cilia and flagella-associated protein 61,
Cfap61. Of note, cilia have been shown to drive adipogenesis through ciliary signaling in preadipocytes.
39- Hilgendorf K.I.
- Johnson C.T.
- Mezger A.
- et al.
Omega-3 fatty acids activate ciliary FFAR4 to control adipogenesis.
Stimulation of enterocytic lipid resynthesis could be another compensatory mechanism, owing to the ileal increased expression of the acyl-CoA cholesterol acyl transferase
Acat2, known to esterify sterols in enterocytes’ endoplasmic reticulum.
40Intestinal lipid absorption and lipoprotein formation.
Notably, transcriptome analysis also suggested that carbohydrate metabolism and intestinal gluconeogenesis is promoted in the ileum of EEC-deficient mice. Transcriptional adaptation in the ileum might result from the upregulation of the lipogenic and glycolytic transcription factor Mlixpl/ChREBP, thought to integrate metabolic signals to modulate its activity according to nutritional status.
41- Abdul-Wahed A.
- Guilmeau S.
- Postic C.
Sweet sixteenth for ChREBP: established roles and future goals.
Thus, intestinal transcriptomic alterations in EEC-deficient mice reflect the lack of enteroendocrine secretion and impaired lipid metabolism and adaptative responses promoting lipogenesis and carbohydrate metabolism to overcome lipid malabsorption failures and improve energy homeostasis. Importantly, the dramatic weight loss of EEC-deficient mice fed an HFD could result from their limited possibility to shift to carbohydrates as an alternative energy source. Indeed, the transcriptomic analysis revealed that the expression of genes controlling the uptake and processing of dietary carbohydrates are not enhanced unlike CD nutrition conditions.
The phenotyping of mice lacking EECs and corresponding hormones confirmed the role of enteroendocrine hormones in lipid metabolism. However, we profiled their transcriptome in wild-type mice to gain further insights into EEC’s response to an HFD. We found that transcripts for
Pyy,
Sct,
Cck,
Nts, Gip, and Vgf hormones were increased, suggesting the need to overproduce them to compensate for a high-fat intake. These findings suggest a previously unappreciated role of GIP and VGF in EEC response to lipids. Indeed, the function of GIP in lipid metabolism has not been extensively studied, but recent findings report that GIP receptor-deficient mice, fed an HFD, have reduced adiposity.
42- Boer G.A.
- Keenan S.N.
- Miotto P.M.
- et al.
GIP receptor deletion in mice confers resistance to high-fat diet-induced obesity via alterations in energy expenditure and adipose tissue lipid metabolism.
VGF is co-expressed with serotonin in intestinal enterochromafin cells.
12- Piccand J.
- Vagne C.
- Blot F.
- et al.
Rfx6 promotes the differentiation of peptide-secreting enteroendocrine cells while repressing genetic programs controlling serotonin production.
Hypothalamic VGF is important for appetite regulation and VGF-deficient mice have reduced stored lipid in white adipose tissue.
43- Benchoula K.
- Parhar I.S.
- Hwa W.E.
The molecular mechanism of VGF in appetite, lipids, and insulin regulation.
Nevertheless, the role of intestinal VGF and derivative peptides is unknown. Besides hormone transcripts, we also found that genes involved in lipid metabolism (
ApoA5,
Ptgs2/
Cox2,
MglII/Magl) were upregulated specifically in EECs in HFD conditions but not in other intestinal epithelial cells. In summary, we established the transcriptomic signature of EECs in response to an HFD, thereby generating novel hypothesis on their specific role in lipid metabolism.
The gut microbiota is an essential regulator of nutrient digestion, absorption, and energy harvest.
44- Bäckhed F.
- Ding H.
- Wang T.
- et al.
The gut microbiota as an environmental factor that regulates fat storage.
,45- Martinez-Guryn K.
- Hubert N.
- Frazier K.
- et al.
Small intestine microbiota regulate host digestive and absorptive adaptive responses to dietary lipids.
Alterations of gut microbiota have been shown in diverse metabolic diseases.
46Gut microbiota in human metabolic health and disease.
Various metabolites, including short chain fatty acids and products of dietary fiber fermentation, mediate its metabolic role. Short-chain fatty acids serve as an energy source but can also signal via fatty acid receptors at the surface of EECs modulating both their hormonal secretion
47- Psichas A.
- Sleeth M.L.
- Murphy K.G.
- et al.
The short chain fatty acid propionate stimulates GLP-1 and PYY secretion via free fatty acid receptor 2 in rodents.
,48- Tolhurst G.
- Heffron H.
- Lam Y.S.
- et al.
Short-chain fatty acids stimulate glucagon-like peptide-1 secretion via the G-protein–coupled receptor FFAR2.
and transcriptome.
49- Arora T.
- Akrami R.
- Pais R.
- et al.
Microbial regulation of the L cell transcriptome.
Whether, in turn, enteroendocrine hormones influence microbiota composition has not been studied directly. Analysis of the microbiota in the stools of mice lacking EECs and fed a regular chow diet revealed a severe drop in alpha- and beta-diversity, 4 weeks after the induction of EECs depletion and before weight loss.
Bacteroides and
Lactobacillus were found progressively enriched, whereas several Firmicutes genera including
Lachnospiraceae and
Muribaculaceae decreased. An increase in gram negative bacteria
Proteobacteria was also observed in EEC-deficient mice. Microbiota remodeling in EEC-deficient mice could be either causal to lipid malabsorption or secondary to gut dietary content changes due to residual unabsorbed lipids, which would favor the growth of subsets of bacteria. Indeed, the diet is essential for microbial composition as both starvation
50- Jawahar J.
- McCumber A.W.
- Lickwar C.R.
- et al.
Starvation causes changes in the intestinal transcriptome and microbiome that are reversed upon refeeding.
and HFD
51- David L.A.
- Maurice C.F.
- Carmody R.N.
- et al.
Diet rapidly and reproducibly alters the human gut microbiome.
, 52- Hildebrandt M.A.
- Hoffmann C.
- Sherrill-Mix S.A.
- et al.
High-fat diet determines the composition of the murine gut microbiome independently of obesity.
, 53- Turnbaugh P.J.
- Bäckhed F.
- Fulton L.
- et al.
Diet-induced obesity is linked to marked but reversible alterations in the mouse distal gut microbiome.
have been shown to impact microbiota communities. For instance, HFD decreased
Bacteroidetes, whereas
Firmicutes and
Proteobacteria increased.
52- Hildebrandt M.A.
- Hoffmann C.
- Sherrill-Mix S.A.
- et al.
High-fat diet determines the composition of the murine gut microbiome independently of obesity.
The
Firmicutes/Bacteroidetes ratio was, however, unaltered in mice lacking EECs and fed with a CD. Thus, despite possible increased lipids in the gut lumen, microbiota changes differ from those reported in an HFD. Like mice lacking all EECs, this increased
Firmicutes/Bacteroidetes ratio was also not observed in mice lacking neurotensin, which actually attenuated HFD-induced intestinal dysbiosis. Like EEC-deficient mice, neurotensin-deficient mice have reduced intestinal fat absorption.
54- Li J.
- Li X.
- Song J.
- et al.
Absence of neurotensin attenuates intestinal dysbiosis and inflammation by maintaining Mmp7/α-defensin axis in diet-induced obese mice.
,55- Li J.
- Song J.
- Zaytseva Y.Y.
- et al.
An obligatory role for neurotensin in high-fat-diet-induced obesity.
Changes in gut microbiota composition may result from changes in the host ecosystem such as transit time, pH. Studies suggest that
Firmicutes and
Proteobacteria are more tolerant to these factors.
56- Donaldson G.P.
- Lee S.M.
- Mazmanian S.K.
Gut biogeography of the bacterial microbiota.
Interestingly, the expression of genes encoding antimicrobial proteins such as
Reg3b,
Reg3g, or their upstream regulator
Retlnb57- Hogan S.P.
- Seidu L.
- Blanchard C.
- et al.
Resistin-like molecule β regulates innate colonic function: Barrier integrity and inflammation susceptibility.
was increased in the small intestine of EEC-deficient mice (both in CD and HFD). These antimicrobial peptides can be activated by intestinal bacteria to promote the spatial segregation of the microbiota from the epithelium.
58- Propheter D.C.
- Chara A.L.
- Harris T.A.
- et al.
Resistin-like molecule β is a bactericidal protein that promotes spatial segregation of the microbiota and the colonic epithelium.
Interestingly, alteration of the microbiota composition can induce increased expression of
Reg3g, which is required for microbiota beneficial effect on body weight and energy balance.
59- Shin J.H.
- Bozadjieva-Kramer N.
- Shao Y.
- et al.
The gut peptide Reg3g links the small intestine microbiome to the regulation of energy balance, glucose levels, and gut function.
Higher levels of antimicrobial peptides in the duodenum could also suggest an increase in bacterial load in the duodenum of EEC-deficient mice. Lipid malabsorption and an increased intestinal anti-microbial response were also reported in mice lacking Arx transcription factor, which is necessary for the specification of Sct-, Cck-, Glp1-, Glp2-, Pyy-, Nts-, and Gip-expressing cells.
38- Terry N.A.
- Ngaba L.V.
- Wilkins B.J.
- et al.
Lipid malabsorption from altered hormonal signaling changes early gut microbial responses.
Despite the severe alteration of gut microbiota diversity in EEC-deficient mice, fecal short chain fatty acids were unaltered. Recent studies suggest that the small intestinal microbiota adapts to dietary lipid variations and plays an important role in lipid digestion and absorption.
45- Martinez-Guryn K.
- Hubert N.
- Frazier K.
- et al.
Small intestine microbiota regulate host digestive and absorptive adaptive responses to dietary lipids.
Whether microbiota communities and their metabolites are modified in the small intestine of EEC-deficient mice remains to be explored.
We further investigated the metabolic alterations of EEC-deficient mice using PICRUSt2 approach to predict the functional capability of their microbiota. The loss of EEC induced a decreased ability to maintain bacterial homeostasis as indicated by a lower abundance of pathways involved in chemotaxis, bacterial motility, and biosynthesis of antimicrobial products (beta-lactams and ansamycins). This suggests that the gut microbiota of EEC-deficient mice may have a lower fitness that could predispose to later infection.
60- Akahoshi D.T.
- Bevins C.L.
Flagella at the host-microbe interface: key functions intersect with redundant responses.
Among different factors, the nutritional quality of the carbon source could participate in the alterations of these biological processes.
61- Colin R.
- Ni B.
- Laganenka L.
- et al.
Multiple functions of flagellar motility and chemotaxis in bacterial physiology.
The loss of EEC likely led the gut microbiota of EEC-deficient mice to shift towards increasing carbohydrate metabolism, including gluconeogenesis and lipoic acid metabolism, to find another source of carbon and energy. Interestingly, a metabolic shift towards carbohydrate metabolism was also suggested by the mouse transcriptomic data. Metabolic adaptation of the gut microbiota may participate in host adaptation following tamoxifen treatment. Also, inferred xenobiotic degradation pathways (aminobenzoate degradation, chlorocyclohexane and chlorobenzene degradation, nitrotoluene degradation) were over-represented in mutant mice. They may witness an increase in toxic compounds and optimization of energy source availability. Several predictive pathways (glutathione metabolism, selenocompound metabolism, and lipoic acid metabolism) increased in EEC-deficient mice were related to redox homeostasis balance and may participate in metabolic deregulation and oxidative stress.
62- Lennicke C.
- Rahn J.
- Kipp A.P.
- et al.
Individual effects of different selenocompounds on the hepatic proteome and energy metabolism of mice.
In particular, glutathione has been reported to affect bacterial pathogenesis and virulence.
63New roles for glutathione: modulators of bacterial virulence and pathogenesis.
In summary, the predictive functional analysis highlighted the alteration of several metabolic pathways that participate in both gut microbiota and host adaptation to the loss of EEC. However, it is still challenging to discriminate the pathways that compensate for the loss of energy source availability from pathogenic processes.
Overall, our studies show that in adult mice: (1) EECs are not essential for survival when mice are fed with a standard CD but are vital under high-fat nutrition, which triggers EEC-specific transcriptional changes in wild-type mice; (2) EECs are essential to optimize lipid absorption and food efficacy and maintain weight under a standard diet; (3) in the absence of EECs, intestinal lipid absorption is spatially delayed, and transcriptional changes suggest adaptation mechanisms to remedy lipid malabsorption and adjust energy utilization; (4) blocking EEC formation rapidly leads to a severe drop of gut microbiota diversity and remodeling. These studies using mouse models may contribute to shedding light on the pathophysiology of human anendocrinosis and adaptive mechanisms that could be targeted to improve intestinal malabsorption and limit the need for parenteral nutrition.
Material and Methods
Study Approval
All animal studies were approved by the French Ministry of “Enseignement supérieur de la Recherche et de l’Innovation” ethical committee, in compliance with the European legislation on care and use of laboratory animals (agreement no. APAFIS#22786-2019111511478935 v3).
Animals and Animal Handling
Mice had ad libitum access to water and 2 different types of diets from weaning: standard CD (Safe D04 diet – Safe Diets) or a rodent HFD with 60 kcal % fat (D12492 – Research Diets). Animals were maintained on a 12-hour light/dark cycle.
Inducible intestine-specific Neurog3-deficient mice, Neurog3
ΔAdInt, were generated by crossing Neurog3
flox/flox mice
13- Mellitzer G.
- Beucher A.
- Lobstein V.
- et al.
Loss of enteroendocrine cells in mice alters lipid absorption and glucose homeostasis and impairs postnatal survival.
with the Villin-CreER
T2 line
(gift of Dr S. Robine). Mice were maintained on C57BL/6 background. Genotyping was performed as previously described.
13- Mellitzer G.
- Beucher A.
- Lobstein V.
- et al.
Loss of enteroendocrine cells in mice alters lipid absorption and glucose homeostasis and impairs postnatal survival.
In all this study, Neurog3
flox/flox; Villin-Cre
ERT2 animals were referred to as mutants or Neurog3
ΔAdInt and Neurog3
flox/flox littermates as controls. To induce Neurog3 deletion, 8- to 13-week-old males were treated with tamoxifen (10 mg/30 g body weight mouse) (Sigma-Aldrich) by oral gavage once per day, every second day for 5 days. Mutants and controls were treated with tamoxifen concurrently. Studies were performed with adult males.
Metabolic Studies
Metabolic studies were performed by the Mouse Clinical Institute (MCI) metabolic platform (
http://www.ics-mci.fr). Body weights were recorded weekly all over the study. At 4 and 10 weeks after the tamoxifen treatment, body composition was determined by magnetic resonance imaging (Minispec+ analyzer, Bruker) and normalized to the weight. At 5 weeks after tamoxifen, food and water consumption as well as stool and urine production were measured over 24 hours in metabolic cages (Tecniplast) after 3 days of acclimatation. The animals were housed individually for 4 days with ad libitum access to food and water. The energy content of stools was evaluated in a bomb calorimeter (C503 control, IKA). The energy excreted was calculated as feces energy content (Cal/g) × feces weight (g). The energy ingested was calculated as food intake per day (g) × diet calorific value (3.34 kcal/g). Food efficiency is the ratio of energy ingested/energy excreted per day.
Blood Analyses
All blood analyses were performed by the MCI metabolic platform (
http://www.ics-mci.fr). At the end of the metabolic tests, a terminal intracardiac blood sample was taken. Different blood chemistry parameters and a panel of metabolic markers were measured on the plasma samples: total cholesterol, HDL and LDL cholesterol, triglycerides, free fatty acids, glycerol, biliary salts, and lipase.
Short-chain Fatty Acids Measurements
Short-chain fatty acids (acetate, propionate, butyrate, isobutyrate, and valerate) contents were determined in stool samples by gas chromatography-mass spectrometry (GC-MS). Stool samples were weighed (20–50 mg) and dissolved in 200 μL of sodium hydroxide solution at 0.15 M (Sigma-Aldrich). Short-chain fatty acids dosage was performed as described in Cuffaro et al (PMID: 33803291).
Sterol Measurements
Sterol (cholesterol, desmosterol, lathosterol, campesterol, lanosterol, and sitosterol) contents were determined in stool samples by GC-MS. For quantification, standard solutions of sterols (ranging from 0–500 μg/mL for cholesterol and 0-2 μg/mL for other sterols) were made by serial dilutions of 500 μg/mL and 2 mg/mL (for cholesterol and others sterols, respectively) solution in cyclohexane. Stool samples were weighed (20–50 mg) and diluted 1/10 in distilled water. Sterols were extracted from stool homogenates and standard solutions (100 μL) with 1.5 mL of a cyclohexane and 20 μl of Internal Standards Solution (D7-cholesterol [500 μg/mL] and D6-campesterol [2 μg/mL]). After centrifugation (10 minutes, 4100 rpm, 15 °C), supernatants were dried under a nitrogen stream at 60 °C. Dried samples were dissolved in 500 μL of ethanol 90% and 250 μL of NaOH 1N and heated at 100 °C for 3.5 hours for saponification. Distilled water was added (250 μL), and sterols were extracted with 1.75 mL cyclohexane. After centrifugation (5 minutes, 3000 rpm, 15 °C), supernatants were collected and dried under nitrogen at 60 °C. Dried samples were derivatized with 40 μL of N-methyl-N-trimethylsilyl-trifluoroacetamide (40 minutes, 70 °C, Sigma Aldrich). Samples were evaporated to dryness under a nitrogen stream (60 °C) and reconstituted with cyclohexane (250 μL). GC-MS analyses were performed using 6890 and 5975C inert systems (Agilent Technologies). Sterols and internal standards were eluted on a HP-1MS column (Agilent Technologies, 30 m × 0.25 mm ID × 0.25 μm film) with helium used as gas vector at 1.2 mL/min flow rate. Five microliters were injected into the analytical systems, the injector being set at 250 °C. Separation was achieved using a temperature gradient (70 °C/min from 150 °C to 285 °C, 20 °C/min to 300 °C, 300 °C for 7.5 minutes), and MS ran on electron impact mode for ionization. Single-ion monitoring mode was used at m/z ratios of 301, 308, 343, 458, 382, 393, 349, and 396 for cholesterol, D7-cholesterol, desmosterol, lathosterol, campesterol, lanosterol, D6-campesterol, and sitosterol, respectively.
Immuno-histochemistry and Histopathology
Mouse tissues were fixed in paraformaldehyde 4% at 4 °C overnight and processed as previously described.
12- Piccand J.
- Vagne C.
- Blot F.
- et al.
Rfx6 promotes the differentiation of peptide-secreting enteroendocrine cells while repressing genetic programs controlling serotonin production.
,13- Mellitzer G.
- Beucher A.
- Lobstein V.
- et al.
Loss of enteroendocrine cells in mice alters lipid absorption and glucose homeostasis and impairs postnatal survival.
Primary antibodies used for immunostainings: goat anti-ChromagraninA (ChrA (C20) sc-1488; Santa Cruz Biotechnology); goat anti-secretin ([S-21] sc-26630; Santa Cruz Biotechnology); rat anti-somatostatin (YC7, MAB354; Chemicon). Secondary antibodies are conjugated to Alexa Fluor 594, DyLight 488, DyLight 549, and DyLight 649 (Jackson ImmunoResearch). Standard histology techniques were performed in collaboration with the MCI (
http://www.ics-mci.fr).
Morphometric Analyzes
The length of villi (at least 69 units under a CD; at least 20 under an HFD) and colon glands (at least 63 units under a CD) were measured per genotype in well oriented crypt-villus units using the Fiji software (n = 3–8). Evaluation of Oil Red O staining in intestinal epithelium and lamina propria was performed on 2 sections per segment per animal (n = 5). Blind scoring was performed according to size and frequency as follows: absence (0), rare very small droplets (0.25), rare small droplets (0.5), frequent small droplets (1), and ubiquitous big droplets (2).
RNA isolation and Reverse Transcriptase-quantitative Polymerase Chain Reaction
The intestine from adult Neurog3ΔAdInt and control males was collected 18 weeks (males under CD) and 5 weeks (males under HFD) after tamoxifen administration. Total RNA from duodenum, jejunum, ileum, proximal colon, or distal colon was extracted using NucleoSpin RNA Plus, Mini Kit (Macherey-Nagel) according to the manufacturer’s instructions. Reverse transcription was performed using Transcriptor Reverse Transcriptase (Roche) and Random primers (Roche). Quantitative polymerase chain reactions were performed in a Light Cycler 480 (Roche) with Light Cycler 480 Probes Master Mix (Roche) using mouse-specific TaqMan primers and probes (ThermoFisher) recognizing Neurog3 (Mm00437606_s1), Cck (Mm00446170_m1), ChgA (Mm00514341_m1), Defa1 (Mm02524428_g1), Gcg (Mm00801712_m1), Ghrl (Mm00445450_m1), Gip (Mm00433601_m1), mKi67 (Mm01278617_m1), Lgr5 (Mm00458299_m1), Lyz1 (Mm00657323_m1), Muc2 (Mm00458299_m1), Nts (Mm00481140_m1), Olfm4 (Mm01320260_m1), Pcna (Mm00448100_g1), Pyy (Mm00520715_m1), Sct (Mm00441235_g1), Sst (Mm00436671_m1), Tph1 (Mm00493794_m1). Gene expression levels were normalized to Actb (4352933E) or Rplp0 (Mm01974474_gH). Relative changes in gene expression were determined by the 2-ΔΔCt method.
RNA Sequencing and Bioinformatic Analysis
Neurog3
eYFP/+ mice
66- Mellitzer G.
- Martín M.
- Sidhoum-Jenny M.
- et al.
Pancreatic islet progenitor cells in neurogenin 3-yellow fluorescent protein knock-add-on mice.
(C57BL/6N background), 8 weeks old, were fed with standard CD (Safe D04) or HFD (D12492 Research diet) for 4 weeks. The jejunum was minced and villi enriched by EDTA treatment. Cells were dissociated with dispase and eYFP+ and eYFP- cells sorted (FACS Aria II), RNA was extracted using RNeasy Qiagen columns. Double-stranded cDNAs were prepared from total RNA using Ovation RNA-Seq System V2 (NuGEN Part #7102-32). Libraries were then generated using Ovation SP Ultralow Library Systems (NuGeN, Part #8033) following manufacturer instructions and sequenced on an Illumina HiSeq2500 sequencer (single-end 50bp reads). Reads were mapped onto the mm10 assembly of mouse genome using Tophat v2.0.10
67- Kim D.
- Pertea G.
- Trapnell C.
- et al.
TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions.
with bowtie v2.1.0.
68- Langmead B.
- Salzberg S.L.
Fast gapped-read alignment with Bowtie 2.
Gene expression was quantified using HTSeq version 0.5.4p3 and gene annotations from Ensembl release 74. Statistical analysis was performed with R 3.0.1 and DESeq2 1.0.19
69- Love M.I.
- Huber W.
- Anders S.
Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2.
Bioconductor library with Benjamini and Hochberg (BH)
P-value adjustment method.
RNA sequencing was carried out on the duodenum and ileum from adult Neurog3
ΔAdInt and control males (n = 4/genotype) fed with standard CD or HFD. After tamoxifen administration, samples were collected upon 18 weeks (CD) or 5 weeks (HFD) after. Libraries were prepared using TruSeq Stranded Total RNA Library Prep Gold kit (Illumina) with either rRNA Removal Mix-Gold (Illumina) or riboPOOL kit (siTOOLs Biotech) at the step of ribodepletion and sequenced on an Illumina HiSeq 4000 or NextSeq 2000 (single-end 50bp reads). Reads were mapped onto the mm39 assembly of mouse genome using STAR
70- Dobin A.
- Davis C.A.
- Schlesinger F.
- et al.
STAR: ultrafast universal RNA-seq aligner.
version 2.5.3a. Gene expression was quantified using htseq-count
71- Anders S.
- Pyl P.T.
- Huber W.
HTSeq—a Python framework to work with high-throughput sequencing data.
version 0.6.1p1 and gene annotations from Ensembl release 104. Statistical analysis was performed using R 3.3.2 and DESeq2 1.16.1 Bioconductor library and independent hypothesis weighting
72- Ignatiadis N.
- Klaus B.
- Zaugg J.B.
- et al.
Data-driven hypothesis weighting increases detection power in genome-scale multiple testing.
P-value adjustment method. Gene Ontology enrichment analysis was performed using METASCAPE. Data were deposited in GEO repository (GSE218600, GSE224026, GSE221473).
High Throughput 16S Metagenomic Analysis
Samples and library preparations and preliminary bioinformatic analyses were performed by the MCI (
http://www.ics-mci.fr). The feces of each mouse were sampled before and 8, 15, 22, and 30 days after tamoxifen administration. DNA was extracted using Nucleospin DNA stool kit (Macherey-Nagel). The region V3-V4 of the 16S rRNA gene was amplified by polymerase chain reaction and paired-end sequenced (2 250bp) using MiSeq Illumina (University of Laval, Canada). FROGS (galaxy-supported pipeline) was used to filter the sets of 16S V3-V4 amplicon sequences described in Escudié et al.
73- Escudié F.
- Auer L.
- Bernard M.
- et al.
FROGS: Find, Rapidly, OTUs with Galaxy Solution.
FROGS was used to produce abundance tables of OTUs and their taxonomic affiliation.
73- Escudié F.
- Auer L.
- Bernard M.
- et al.
FROGS: Find, Rapidly, OTUs with Galaxy Solution.
SILVA was used as reference database for taxonomy affiliation.
74- Quast C.
- Pruesse E.
- Yilmaz P.
- et al.
The SILVA ribosomal RNA gene database project: improved data processing and web-based tools.
Alpha diversity was estimated by the observed OTUs and Shannon index metrics and visualized using R package dplyr v1.0.7, ggplot2 v3.3.5, and ggstatsplot v0.7.2. α-diversity was compared using the Wilcoxon rank-sum test. Beta-diversity ordinations were calculated on the Bray-Curtis dissimilarity matrix using Phyloseq v1.30.0
75phyloseq: An R package for reproducible interactive analysis and graphics of microbiome census data.
and ggplot2 packages. Difference between groups was tested using permutational multivariate analysis of variance (PERMANOVA) computed with the Adonis function of the vegan package 2.5.7 with default parameters. Differential abundance analyses were performed at the OTUs level with DESeq2 using Phyloseq v1.30.0 and DESeq2 v1.26.0. The size factor estimation was set to use “poscounts” as suggested for microbiota analysis. Only results displaying an adjusted
P-value < .05 after BH correction were reported and visualized using ggplot2.
Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt2, https://github.com/picrust/picrust2; v2.5.1) was used to predict and investigate the metagenome of the bacteria identified using 16S sequences. The analysis was performed using default parameters.
76- Douglas G.M.
- Maffei V.J.
- Zaneveld J.R.
- et al.
PICRUSt2 for prediction of metagenome functions.
Pathway abundance predictions were computed based on Kyoto Encyclopedia of Genes and Genomes database. ALDEX2 (
https://github.com/ggloor/ALDEx2)
77- Fernandes A.D.
- Macklaim J.M.
- Linn T.G.
- et al.
ANOVA-Like Differential Expression (ALDEx) Analysis for Mixed Population RNA-Seq.
Cliquez ou appuyez ici pour entrer du texte. statistical analysis was performed to compare pathway abundances of centered-log-ratio (CLR)-transformed data, considering Welch’s
t-test with BH correction. Differentially inferred pathways were visualized in MA plot. Only pathways with a BH adjusted
P-value < .05 were considered significant between the 2 conditions and reported in the bar plot.
Statistics
All data are expressed as mean ± standard deviation (SD). Statistical analysis was performed on GraphPad Prism 8.0 or R software with significance set at P-value < .05. Longitudinal data and grouped analyzes with 2 factors were analyzed using 2-way analysis of variance (ANOVA) and Holm-Sidak posthoc test. Single comparisons were made with Mann-Whitney tests. The statistical tests and sample sizes are indicated in the figure legends.
Acknowledgment
The authors thank Guillaume Pavlovic and Loic Lindner, Marie-France Champy and Erwan Grandgirard (Institut Clinique de la Souris) for help in analyzing the microbiota, mice phenotyping and imaging, Cédric Lemay and Filipe de Vadder for helpful discussions, and GG lab members for critical reading of the manuscript.
CRediT Authorship Contributions
Gerard Gradwohl, PhD (Conceptualization: Lead; Formal analysis: Lead; Funding acquisition: Lead; Supervision: Lead; Validation: Lead; Writing – original draft: Lead)
Florence Blot (Formal analysis: Lead; Methodology: Lead; Validation: Lead)
Justine Marchix (Data curation: Lead; Formal analysis: Lead; Software: Lead; Writing – original draft: Lead)
Miriam Ejarque (Formal analysis: Lead; Methodology: Lead; Validation: Lead)
Sara Jimenez Correa (Data curation: Lead; Formal analysis: Lead; Methodology: Lead; Validation: Lead; Visualization: Lead; Writing – original draft: Lead)
Aline Meunier (Methodology: Lead; Validation: Lead; Visualization: Lead)
Céline Keime (Data curation: Lead; Formal analysis: Lead; Methodology: Lead; Validation: Lead; Visualization: Lead)
Camille Trottier (Formal analysis: Equal)
Mikaël Croyal (Formal analysis: Lead; Methodology: Lead; Validation: Lead)
Céline Lapp (Methodology: Supporting)
Maxime M. Mahe (Funding acquisition: Lead; Methodology: Lead; Supervision: Lead; Writing – original draft: Lead)
Adèle De Arcangelis (Formal analysis: Lead; Methodology: Lead; Project administration: Lead; Supervision: Lead; Validation: Lead; Writing – original draft: Lead)