1. Developmental Biology and Stem Cells
  2. Microbiology and Infectious Disease
Download icon

A conserved bacterial protein induces pancreatic beta cell expansion during zebrafish development

  1. Jennifer Hampton Hill
  2. Eric A Franzosa
  3. Curtis Huttenhower
  4. Karen Guillemin Is a corresponding author
  1. University of Oregon, United States
  2. Harvard T. H. Chan School of Public Health, United States
  3. The Broad Institute, United States
  4. Canadian Institute for Advanced Research, Canada
Research Article
Cited
3
Views
2,582
Comments
0
Cite as: eLife 2016;5:e20145 doi: 10.7554/eLife.20145

Abstract

Resident microbes play important roles in the development of the gastrointestinal tract, but their influence on other digestive organs is less well explored. Using the gnotobiotic zebrafish, we discovered that the normal expansion of the pancreatic β cell population during early larval development requires the intestinal microbiota and that specific bacterial members can restore normal β cell numbers. These bacteria share a gene that encodes a previously undescribed protein, named herein BefA (β Cell Expansion Factor A), which is sufficient to induce β cell proliferation in developing zebrafish larvae. Homologs of BefA are present in several human-associated bacterial species, and we show that they have conserved capacity to stimulate β cell proliferation in larval zebrafish. Our findings highlight a role for the microbiota in early pancreatic β cell development and suggest a possible basis for the association between low diversity childhood fecal microbiota and increased diabetes risk.

https://doi.org/10.7554/eLife.20145.001

eLife digest

For a long time, genes and carefully orchestrated chemical signals from the mother’s body have been known to shape the earliest phases of an embryo’s development. More recently, researchers have started to appreciate that environmental cues – such as signals from nearby bacteria – also shape animal development.

The body is home to a teeming community of bacteria and other microbes, called the microbiota. Because the gut houses more microbes than any other part of the body, the role that the microbiota plays in gut development has been investigated. Less is known about how the microbiota affects the development of the other organs involved in digestion, such as the pancreas. In developing fish and mammals, the pancreas grows its population of insulin-producing beta cells during the same period of development in which the microbial population of the gut becomes established.

Now, Hill et al. show that certain gut bacteria are necessary for the pancreas to populate itself with a robust number of beta cells during development. Normally, the number of beta cells in zebrafish larvae increases steadily in the first few days after hatching. However, developing zebrafish that were reared in a microbe-free environment maintained the same number of beta cells as they had before hatching. Exposing the microbe-free fish to certain bacteria restored their beta cell populations to normal levels. Further investigation revealed that these bacteria release a protein called BefA that causes the beta cells to multiply.

Some bacteria in humans produce proteins that are similar to BefA. Hill et al. performed experiments that showed that these proteins also stimulate beta cell development in microbe-free fish. Future studies are now needed to investigate the mechanism by which the proteins affect beta cell development, and to find out whether they have the same effect in humans and other animals. This will help us to understand whether a lack of gut microbes could contribute to the development of diseases, such as diabetes, that are characterized by insufficient numbers of beta cells.

https://doi.org/10.7554/eLife.20145.002

Introduction

Host-associated microbes play important roles in the development of animal digestive tracts (Bates et al., 2006; Semova et al., 2012; Sommer and Bäckhed, 2013). Using the gnotobiotic zebrafish model, our group has shown previously that resident microbes promote host processes in the developing intestine such as epithelial differentiation (Bates et al., 2006) and proliferation (Cheesman et al., 2011). The role of microbes in the development of other digestive organs remains underexplored, despite the fact that many diseases in peripheral digestive organs are correlated with microbial dysbiosis (Chang and Lin, 2016; Gülden et al., 2015). The ability to manipulate resident microbes in the larval zebrafish (Milligan-Myhre et al., 2011), combined with the optical transparency and sophisticated genetic tools of the zebrafish model, make it a powerful platform to investigate this question. Here, we use gnotobiotic zebrafish to demonstrate a role for resident microbes in promoting pancreatic β cell development.

The zebrafish has a well-characterized program of β cell development, which is highly conserved with that of mammals (Kinkel and Prince, 2009). In the zebrafish embryo, initial β cells arise from precursors within the dorsal and ventral pancreatic buds (Biemar et al., 2001; Field et al., 2003; Wang et al., 2011). The two buds fuse by 52 hr post fertilization (hpf), and give rise to the fully fated pancreas with only a single islet of hormone-secreting endocrine cells, by 3 days post fertilization (dpf) (Biemar et al., 2001; Field et al., 2003; Kumar, 2003). Coinciding with the approximate time of larval emergence from the chorion by 3 dpf, these newly fated β cells begin to expand (Chung et al., 2010; Dong et al., 2007; Hesselson et al., 2009; Kimmel et al., 2011; Moro et al., 2009). β cells derived from the dorsal bud become quiescent, while ventral bud derived β cells begin to undergo expansion via mechanisms of both proliferation and neogenesis (Hesselson et al., 2009). Between 3 and 6 dpf, the number of β cells within the primary islet will almost double (Moro et al., 2009). Intestinal colonization with microbes occurs concurrently with this early larval period of β cell expansion. Following development of the gut tube within the sterile embryo, the intestine of the emergent larva becomes open to the environment at both the mouth and the vent by 3.5 dpf, allowing for inoculation by environmental microbes (Bates et al., 2006). Within the larval gut, bacteria proliferate rapidly, such that a single species in mono-association can reach the luminal carrying capacity within several hours (Jemielita et al., 2014).

Human post-natal β cell expansion also occurs concurrently with intestinal tract colonization by commensal microbes. In utero, β cells are produced via differentiation from progenitors (Georgia et al., 2006; Stanger et al., 2007) and at birth this newly fated cell population begins to expand by self-proliferation (Georgia and Bhushan, 2004; Gregg et al., 2012; Kassem et al., 2000; Teta et al., 2007). β cell proliferation rates peak at 2 years of age and then steadily decline (Gregg et al., 2012). By 5 years of age, most of the β cell mass has become slow cycling and will not expand significantly again unless stimulated by elevated metabolic demands, such as obesity or pregnancy. At birth, infants are exposed to their mothers’ vaginal, fecal and skin associated microbes, which immediately begin to colonize the neonatal intestine (Biasucci et al., 2010; Dominguez-Bello et al., 2010; Palmer et al., 2007). By 3 years of age, the composition and complexity of the microbiota typically resembles that of an adult associated community (Murgas Torrazza and Neu, 2011; Palmer et al., 2007; Yatsunenko et al., 2012). However, factors such as diet, birth mode and antibiotic exposure can result in reduced microbial taxonomic diversity during these early years of life (Mueller et al., 2015). Notably, factors that reduce microbiota diversity are also associated with increased risk for diabetes mellitus (Knip et al., 2005). Loss of β cell function through autoimmunity results in abnormal glucose homeostasis and is the cause of type 1 diabetes (T1D) in humans. Recent studies have shown that decreased taxonomic diversity of the intestinal microbiota is correlated with T1D (Brown et al., 2011; Giongo et al., 2011). Indeed, loss of bacterial diversity precedes the onset of T1D in children, and may play a causative role in disease (Kostic et al., 2015).

To our knowledge, no one has yet investigated a role for the gut microbiota in the development of pancreatic β cells. Communication between the intestine and the pancreas is critical for overall homeostasis. The two organs are therefore connected physically, metabolically, and developmentally in order to carry out their essential functions. We propose that this established and important connection might also mediate the influence of resident microbes on developmental processes in the pancreas. Here we examine the effects of microbial colonization on initial expansion of zebrafish primary islet β cells. We find that β cell mass expansion, up to at least 6 dpf, is promoted by the presence of the microbiota. Using a culture collection of zebrafish intestinal bacteria, we show that certain strains can restore β cell expansion in germ free (GF) fish. We report the discovery of a secreted protein, shared among these strains and named herein β cell expansion factor A (BefA) that is sufficient to recapitulate this effect. Homologs of the befA gene are present in the genomes of a subset of human intestinal bacteria, and we show that two of the corresponding proteins share BefA’s capacity to induce β cell expansion in zebrafish.

Results

The microbiota is required for normal expansion of the larval β cell mass

To investigate a possible role for the microbiota in pancreas development and specifically in β cell expansion, we quantified total β cells in GF and conventionally reared (CV) Tg(-1.0insulin:eGFP) fish (diIorio et al., 2002) at 3, 4, 5 and 6 dpf (Figure 1A, Figure 1—source data 1). The number of β cells in CV fish increased steadily from 3 to 6 dpf (Figure 1A). However, the average number of β cells in GF fish remained static over this time (Figure 1A). Furthermore, at 6 dpf, the overall structure of β cells within the primary islet also appeared much less densely packed in GF than in CV fish (Figure 1B). This effect is not likely to be due to changes in initial differentiation of the β cell population since the total number of β cells is not different between GF and CV fish at 3 dpf (Figure 1A), a time at which exposure to bacteria is also limited.

The microbiota are required for normal expansion of the larval β cell mass.

(A) Total number of β cells per larva in GF (white box plots) and CV (grey box plots) fish at 3, 4, 5 and 6 dpf. In this, and in all subsequent figures, CV data are shown in grey box plots, and GF data, or statistically similar treatment groups, are shown in white box plots. In all relevant panels and remaining figures, box plot whiskers represent the 95% confidence interval of the data set. Single factor ANOVA indicates that gnotobiology of the fish was significant in determining the number of β cells present (F7=9.01, p=1.45e−8). Labels a, ab and b indicate the results of post hoc means testing (Tukey). The difference between GF and CV cell counts became significant at 6 dpf (t=−5.91, p<0.001). (B) Representative 2D slices from confocal scans through the primary islets of 6 dpf CV and GF Tg(-1.0insulin:eGFP) larvae. Each slice is taken from the approximate center of the islet structure. Insulin promoter expressing β cells are in green and nuclei are blue. Scale bar = 40 μM. (C) The average amount of glucose (pmol) per larva aged 6 dpf (* t17=−3.65, p<0.01).

https://doi.org/10.7554/eLife.20145.003

Because insulin from β cells functions to reduce levels of circulating glucose, we tested whether the β cell deficiency in GF larvae at 6 dpf affected the metabolic function of the fish by measuring free glucose levels. The amount of glucose detected in GF fish was significantly higher than in CV fish (Figure 1C, Figure 1—source data 2). These data suggest that GF fish, with a paucity of β cells, are less efficient at importing and processing glucose from the blood due to lower levels of circulating insulin. This is consistent with previous studies showing free glucose levels in zebrafish larvae to be correlated with β cell numbers (Andersson et al., 2012).

Only specific bacterial members of the zebrafish microbiota are sufficient to rescue normal expansion of the GF β cell mass

We developed an experimental timeline, depicted in Figure 2A, to test the capacity of individual zebrafish bacterial isolates to induce β cell expansion. We derived embryos GF at 0 dpf and allowed them to develop in this environment until after hatching. At 4 dpf, when the GF larvae have a patent gut tube, we inoculated them with defined microbes and/or microbial derived products by adding these directly to the embryo media. The fish were incubated with the treatment of interest for 48 hr before analysis of the β cell mass at 6 dpf.

Specific bacterial members of the zebrafish microbiota are sufficient to rescue normal expansion of the GF β cell mass.

(A) Experimental timeline for all subsequent zebrafish experiments, unless stated otherwise.Experimental manipulations are denoted by red text. Important zebrafish developmental events are denoted by black text. (B) Quantification of β cells in CV, GF and GF larvae treated at 4 dpf with either non-sterile tank water (XGF) or mono-associated with a specific bacterial strain. Bacterial mono-associations are labeled by genus. Different Aeromonas sp are labeled with a number (1, 2 or 3). *p<0.05, **p<0.01, ***p<0.001: Denotes treatment that is significantly different than GF by Tukey analysis. Additionally, here and in all subsequent figures, significant data sets (p<0.05 when compared to GF) are also highlighted as green box plots. (C) Bacterial isolates of the zebrafish gut and related strains are capable of forming mono-associations with larvae from 4 to 6 dpf. Quantification of the colony forming units (CFUs) per gut for each bacterial strain, assayed after 48-hr exposure to GF larvae. Dashed line denotes the limit of detection.

https://doi.org/10.7554/eLife.20145.006

We found that we could rescue β cell numbers to CV levels by the addition of non-sterile, normal fish tank water to GF larvae at 4 dpf (Figure 2B, Figure 2—source data 1), suggesting that development of the normal number β cells is dependent upon microbes or microbial-derived products present in the water. We next inoculated 4 dpf GF larvae with a selection of bacterial isolates from the zebrafish intestine (Stephens et al., 2015) as well as one other related strain (Bomar et al., 2013). We prioritized bacterial strains that were capable of forming robust mono-associations with larvae between 4 and 6 dpf, as measured by the number of bacteria found within the gut at 6 dpf (Figure 2C). We found that the mono-associations with three different species of the genus Aeromonas and one species of the genus Shewanella was sufficient to rescue GF β cell numbers to levels observed in CV fish (Figure 2B, Figure 2—source data 1). Importantly, other isolates such as Vibrio sp. and Delftia sp. were not sufficient to rescue this phenotype (Figure 2B, Figure 2—source data 1), indicating that only specific members of the microbiota are capable of inducing expansion of the β cell mass.

Aeromonas secretes a factor that rescues normal expansion of the GF β cell mass

Bacterial interactions with host organisms often involve secreted molecules. To test whether a secreted bacterial factor(s) could influence β cell expansion, we harvested cell free supernatant (CFS) from overnight cultures of each Aeromonas strain shown to rescue β cell expansion (Figure 2B) and added these to GF larvae at 4 dpf. For each of the three strains of Aeromonas tested, the CFS alone was able to restore β cell numbers in GF fish (Figure 3A, Figure 3—source data 1), indicating that a secreted factor (or factors) produced by these bacteria is (are) sufficient to induce β cell expansion. As a control, we also treated GF fish with CFS from a Vibrio sp. isolate, which colonized the zebrafish gut (Figure 2C, *), but did not induce β cell expansion (Figure 2B, *). We found the number of β cells in fish receiving Vibrio CFS was not significantly different from that of GF fish (Figure 3A, Figure 3—source data 1). Furthermore, the capacity to induce increased β cell numbers was lost when the Aeromonas 1 (A. veronii) CFS sample was treated with proteinase K (Figure 3A, Figure 3—source data 1), indicating that our secreted factor(s) of interest was likely to be a protein. Because of existing genetic reagents available for the A. veronii strain (Bomar et al., 2013), and its capacity to modulate traits of gnotobiotic zebrafish and other hosts (Bates et al., 2006; Cheesman et al., 2011; Graf, 1999; Rolig et al., 2015), we focused on this strain for the remainder of our analysis.

Figure 3 with 1 supplement see all
Aeromonas secretes a factor that rescues normal expansion of the GF β cell mass.

(A) Total β cell numbers in GF, CV and GF fish treated at 4 dpf with different cell free supernatant (CFS) samples. 'Aero.' refers to bacteria of the genus Aeromonas, with each number (1, 2, 3) denoting a separate species. '+ PK' indicates proteinase K addition to the CFS sample prior to treatment. *p<0.05, **p<0.01, ***p<0.001: Denotes treatment that is significantly different than GF by Tukey analysis. (B) Total β cell numbers in CV, GF and GF fish treated at 4 dpf with separate ammonium sulfate fractions (% AS) prepared from the Aeromonas 1ΔT2SS CFS. Note that the 60–80% ammonium sulfate fraction resulted in the greatest increase in β cell numbers. (C) Total β cells in GF, CV and GF fish treated with purified protein. 10165 represents purified protein from the M001_10165 locus. (D) Representative 2D slices from confocal scans through the primary islets of GF, CV and 10165 protein treated Tg(-1.0insulin:eGFP) 6 dpf larvae. Insulin promoter expressing β cells are shown in green and nuclei are blue. Scale bar = 40 μM.

https://doi.org/10.7554/eLife.20145.008

To narrow down the list of candidate proteins secreted by A. veronii, we tested whether the activity was present in the CFS of an A. veroniiΔT2SS mutant strain (Maltz and Graf, 2011) lacking a functional type 2 secretion system (T2SS), one of the major protein secretion pathways of Gram-negative bacteria. Despite the fact that it has a reduced secretome, CFS harvested from the A. veroniiΔT2SS strain was sufficient to rescue GF β cell numbers (Figure 3A, Figure 3—source data 1). Conveniently, this finding significantly decreased the number of candidate secreted A. veronii proteins with β cell expansion capacity. This result also suggested that our protein(s) of interest was secreted through an alternative mechanism.

We next used ammonium sulfate precipitation to further separate proteins within the A. veroniiΔT2SS CFS. Each of the fractions was able to increase β cells in GF fish (Figure 3B, Figure 3—source data 2), suggesting that either A. veroniiΔT2SS produces multiple proteins with this activity, or that the effector was present to some extent within each fraction. Since the 60–80% fraction was able to induce the greatest increase in β cell numbers (Figure 3B), we used mass spectrometry to analyze the content of this fraction, which led to the identification of 163 proteins (Supplementary file 1). To identify promising candidates from this list, we took advantage of the fact that our zebrafish-associated bacterial isolates, for which we have drafted genome sequences (Stephens et al., 2015), differed in their capacity to induce β cells (Figure 2B). Using basic local alignment search tool (BLAST) we identified those proteins from our candidate list that were, first, predicted to be encoded by the genomes of the four bacterial strains with β cell expansion capacity, and second, absent from the strains lacking this capacity. Our analysis identified one single candidate gene, denoted by the locus tag, M001_10165 (10165), predicted to encode a putative protein of 261 amino acids. Consistent with the candidate protein being found in the CFS, the putative protein contained a predicted N-terminal secretion sequence.

To test whether 10165 encoded the secreted protein responsible for inducing β cell expansion, we cloned the gene into an inducible expression vector in E. coli strain BL21, which contains no 10165 homologues in its genome. We expressed and purified the 10165 protein to homogeneity, as confirmed by SDS-page gel electrophoresis (Figure 3—figure supplement 1). Purified protein was added to flasks of 4 dpf GF zebrafish larvae. This treatment was sufficient to rescue β cell numbers to CV levels by 6 dpf (Figure 3C, Figure 3—source data 3). The islets of larvae treated with the purified protein were visibly expanded compared to those of GF animals (Figure 3D). Therefore, we have named this protein β cell expansion factor A (BefA) after its observed activity in zebrafish.

BefA is required for Aeromonas to induce GF β cell expansion

To determine whether the befA (10165) locus is necessary for A. veronii to induce an increase in β cell numbers, we generated an A. veroniiΔbefA mutant strain by replacing the coding region of befA with a chloramphenicol resistance gene. To ensure that the loss of the befA gene would not affect the ability of A. veronii to form mono-associations with larvae, we performed growth and colonization assays and saw no deficiency in either the in vitro growth rate (Figure 4—figure supplement 1A) or the ability of A. veroniiΔbefA to colonize the GF intestine compared to the wild-type (WT) strain (Figure 4A). However, when inoculated in a 1:1 ratio together with A. veroniiWT, the A. veroniiΔbefA strain showed a small yet reproducible fitness disadvantage as measured by colonization level and competition index after 48 hr (Figure 4—figure supplement 1B,C). This result indicates that BefA confers some colonization benefit for A. veronii within the larval gut.

Figure 4 with 1 supplement see all
BefA is required for Aeromonas to induce GF β cell expansion.

(A) Quantification of the colony forming units (CFUs) per gut in GF fish mono-associated (MA) with either wild type (WT) or mutant (△befA) A. veronii strains for 48 hr. Dashed line denotes the limit of detection (B) Total β cells in GF fish that have been mono-associated with △befA, treated with CFS from either WT or △befA, treated with purified BefA, or have been inoculated with a combination of these. **p<0.01, ***p<0.001: Denotes treatment that is significantly different than GF by Tukey analysis.

https://doi.org/10.7554/eLife.20145.013

GF fish were mono-associated with the A. veroniiΔbefA strain, or treated with its CFS from 4 to 6 dpf. Neither treatment was sufficient to rescue β cell numbers to CV levels (Figure 4B, Figure 4—source data 1). However, mono-associations of A. veroniiΔbefA could be complemented in trans with either CFS from A. veroniiWT or purified BefA protein, which resulted in the restoration of the β cell population (Figure 4B, Figure 4—source data 1). Taken together, these data demonstrate that the BefA protein is necessary in an A. veronii mono-association for early β cell expansion and suggests that A. veronii only produces a single effector of host β cell expansion.

BefA facilitates β cell expansion by inducing proliferation

Proliferation is the primary mode of human neonatal β cell expansion (Gregg et al., 2012; Kassem et al., 2000; Teta et al., 2007). In 4–6 dpf zebrafish larvae, proliferation also contributes to β cell expansion (Field et al., 2003; Hesselson et al., 2009; Moro et al., 2009). Therefore, we investigated whether CV larvae had higher levels of β cell proliferation than GF larvae. 4 dpf larvae were treated with the thymadine analog, 5-ethynyl-2’-deoxyuridine (EdU) for 48 hr to mark cells that underwent proliferation during this time window. We found that, by 6 dpf, CV larvae had significantly more EdU labeled insulin-expressing cells than GF larvae (Figure 5A,B, Figure 5—source data 1). Next we asked whether treatment of GF larvae with BefA was sufficient to restore β cell proliferation to CV levels. We found that BefA-treated GF larvae had EdU incorporation similar to CV fish and significantly greater than GF (Figure 5A,B, Figure 5—source data 1). CFS from our A. veroniiΔbefA strain was not sufficient to increase proliferation rates in GF fish (Figure 5B, Figure 5—source data 1). Our results show that BefA is sufficient to increase cell proliferation that gives rise to an expanded β cell population during early larval development. Furthermore, BefA seems to be the only product of the A. veronii CFS that is capable of inducing this cell proliferation.

Figure 5 with 1 supplement see all
BefA facilitates β cell mass expansion through proliferation.

(A, D & E) Representative 2D slices from confocal scans through the primary islets of GF, CV and BefA (10165) protein treated 6 dpf larvae. Scale bars = 40 μM. (A) Insulin promoter expressing β cells are shown in green, all nuclei are blue, and EdU containing nuclei are magenta. Left hand panels are a merge of all three markers. For ease of resolving cells that are double positive for both insulin and EdU, the right hand panels show the location of insulin outlined by white dashed lines. (B) Percentage of EdU positive β cells in CV, GF or GF treated with either purified BefA or CFS from A. veroniiΔbefA cultures (△befA CFS). ***p<0.001: Denotes treatment that is significantly different than GF by Tukey analysis. (C) Total EdU positive exocrine cells quantified from the approximate central longitudinal plane of the pancreas in each fish. (D) Insulin promoter expressing β cells are shown in green, all nuclei are blue, and α cells, stained with anti-glucagon antibody are magenta. (E) Somatostatin promoter expressing δ cells are shown in white, all nuclei are blue, and β cells stained with anti-insulin antibody are outlined in green. (F) Total α cells in GF, CV and GF fish treated with BefA. (G) Total δ cells in GF, CV and GF fish treated with BefA.

https://doi.org/10.7554/eLife.20145.016

In zebrafish larvae, both the proliferation of existing β cells as well as the proliferation of progenitors contribute to the expansion of β cells that occurs between 4 and 6 dpf (Dong et al., 2007; Field et al., 2003). Because our 48-hr EdU pulse labeled β cells born from both events, our experiment did not distinguish the exact cell population undergoing proliferation in response to BefA. Due to their low rates of proliferation, dividing β cells are difficult to detect without pulse labeling. Neogenesis of β cells from progenitors is also rare, but can be detected as the appearance of insulin positive cells in the extra-pancreatic duct (EPD) (Dong et al., 2007; Hesselson et al., 2009). We quantified insulin expressing cells in the EPD in 6 dpf CV and GF larvae. In a survey of over 500 Tg(-1.0insulin:eGFP) larvae, we found a slight but significant increase in EPD-localized insulin expressing cells in CV versus GF fish (Figure 5—figure supplement 1), suggesting that the microbiota increases endocrine progenitor proliferation. Whether the microbiota also promote proliferation of mature β cells in the islet and whether BefA promotes the proliferation of one or both of these cell populations remains to be determined.

To test whether BefA activity was specific to endocrine tissue, or whether it acts as a nonspecific pro-proliferative stimulant in the pancreas, we analyzed its ability to induce proliferation in exocrine pancreatic tissue by treating Tg(ptf1a:eGFP) larvae (Thisse et al., 2004) with EdU and BefA from 4 to 6 dpf and quantifying proliferative eGFP positive cells. We found no difference in the level of exocrine cell proliferation across GF, CV and BefA treatments (Figure 5C). To test whether β cells were the only endocrine cell type in the islet to be responsive to BefA, we also quantified the total number of glucagon-expressing α (Figure 5D) and somatostatin-expressing δ (Figure 5E) cells in GF, CV and BefA treated fish. We again found no difference in the total numbers of these cells across treatments (Figure 5F,G). These results suggest that in the pancreas, β cells alone are responsive to the presence of BefA.

BefA homologs are produced by members of the human gut microbiota and have conserved function

We wondered if BefA-like proteins are produced by the human microbiota. Phylogenetic analysis of related sequences in bacterial genomes uncovered close homologs (at least 82% amino acid sequence identity) in many, but not all, species of the Aeromonas, Vibrio, and Photobacterium genera. We also found an example of a highly related sequence in the human-associated species Enterococcus gallinarum, which was likely acquired through a horizontal gene transfer event (Figure 6A). Widening the search to include more distant homologs identified potentially related genes in three additional human-associated genera: Enterobacter, Escherichia, and Klebsiella (Figure 6B).

Figure 6 with 1 supplement see all
Homologs of BefA encoded in the human microbiome have conserved function in zebrafish.

(A) Close homologs of BefA across microbial species. Each species is represented by its closest BefA homolog, with a minimum allowed amino acid sequence identity of 50% (relative to the query sequence). Notably, the Enterococcus gallinarum homolog clusters among homologs from the Aeromonas genus, which is evidence of a possible lateral gene transfer event. (B) A view of the BefA phylogeny including more distant homologs (sequence identity >20%) and grouped by genus. The portion of the tree represented in A is contained in the light gray box. In both panels, red numbers indicate branch support (values closer to 1 are better supported); branches with support values <0.5 have been collapsed. Blue clades indicate genera that were associated with humans in metagenomes produced during the Human Microbiome Project (HMP). Black arrowheads indicate genera tested for functional conservation in panel D. Scale bars indicate amino acid substitutions per amino acid site. (C) SDS-page gel: 1 = ladder, 2 = CFS from induction of E. coli BL21 carrying an empty vector, 3 = CFS from induction of E. coli BL21 carrying vector with an Enterococcus gallinarum homolog, estimated size of 29 kDa, lane 4 = CFS from induction of E. coli BL21 carrying vector with Enterobacter aerogenes homolog, estimated size of 21 kDa. White arrows indicate induced proteins. (D) Total β cells in CV, GF and GF fish that have been treated with either induced BL21 E. coli supernatant dominated by the homologous BefA protein encoded from Enterococcus gallinarum (E. gal. homolog) and Enterobacter aerogenes (E. aero. homolog), or induced supernatant from an empty vector control. *p<0.05, **p<0.01, ***p<0.001: Denotes treatment that is significantly different than GF by Tukey analysis.

https://doi.org/10.7554/eLife.20145.019

We tested whether representative BefA-like proteins from human-associated bacteria had the capacity to induce β cell expansion in our gnotobiotic zebrafish model. We cloned into BL21 E. coli two befA-like genes: the more closely related homologue from Enterococcus gallinarum and a more distantly related homologue from Enterobacter aerogenes. The amino acid sequence alignment of these two homologs against the Aeromonas BefA sequence is shown in Figure 6—figure supplement 1. Both the Aeromonas and Enterococcus sequences contain a short N-terminal hydrophobic secretion signal, which is not predicted in the more distant Enterobacter sequence. The most conserved region of these proteins is the C-terminal portion, which contains a putative SYLF domain of unknown function. Induction of expression of each gene in E. coli yielded CFS that were dominated by each of the respective homologous proteins, in contrast to the CFS from control E. coli expressing an empty vector (Figure 6C). Upon addition of these supernatants to GF larval zebrafish, we observed rescue of β cell numbers to the CV level with both the Enterococcus gallinarum and Enterobacter aerogenes proteins, but not the empty vector control (Figure 6D, Figure 6—source data 1). These results indicate that members of human-associated microbiota produce secreted proteins capable of inducing β cell expansion.

Discussion

Using a gnotobiotic zebrafish model, we have discovered a class of proteins produced by resident gut bacteria that have the capacity to increase the expansion of pancreatic β cells during early zebrafish development. BefA and related homologues are predicted to contain a C-terminal SYLF domain, which has been described in proteins from organisms in all kingdoms of life, including humans, but for which little is known functionally beyond a possible role in lipid binding (Hasegawa et al., 2011). Genes encoding BefA and related proteins are found in a small subset of all bacteria genera, with a predominance in genera of host-associated bacteria, but befA homologues are not ubiquitously present in any of these genera.

Our finding of a role for specific secreted bacterial proteins in β cell development raises the possibility of a new link between the resident microbiota and diseases of β cell paucity, such as diabetes mellitus. Type 1 diabetes (T1D), is caused by both genetic and environmental factors, as indicated by the 50% disease discordance among monozygotic twins (Akerblom et al., 2002). One environmental factor associated with T1D is microbiota composition (Gülden et al., 2015). Mechanistic models for the role of the microbiota in T1D etiology have focused on the capacity of the microbiota to modulate the development and function of the immune system, and thus influence the propensity of genetically susceptible individuals to develop autoimmunity to β cell antigens (Gülden et al., 2015). Multiple aspects of host immune cell development and function known to play a role in T1D are altered by the loss of microbes, including development of lymphoid tissue (Macpherson and Harris, 2004) and T cell differentiation and function (Alam et al., 2011; Farkas et al., 2015; Ivanov et al., 2008). We hypothesize an additional role for the early microbiota in establishing the β cell population size that would either buffer against, or render individuals susceptible to, β cell depletion by autoimmune destruction.

In humans, β cells undergo a period of postnatal expansion, before becoming quiescent around age two (Gregg et al., 2012). Differences in β cell growth during this time are thought to account for the wide variation in β cell mass observed in adults (Wang et al., 2015b). The idea that early life β cell census could influence diabetes risk is supported by studies in both rodents and humans, and has been theorized as an important risk factor for type 2 diabetes (Kaijser et al., 2009), a disease which is also influenced strongly by microbiota composition (Cox and Blaser, 2014). Compromised β cell development in rats results in an insufficient number of cells to adequately control glucose metabolism (Figliuzzi et al., 2010). In mice, perinatal β cell proliferation rates can be tuned via the modulation of Gi-GPCR signaling (Berger et al., 2015). Changes to early β cell proliferation capacity in these mice correlates directly with adult β cell mass, which subsequently impacts glucose regulation (Berger et al., 2015). Furthermore, meta analysis of human data has revealed a correlation between an early age of β cell loss and more rapid onset of T1D (Klinke, 2008), consistent with the model that failure to generate a reserve of β cells early in development increases disease risk.

We hypothesize that neonatal microbiomes with a low abundance of BefA equivalents would result in reduced β cell proliferation, lower β cell mass, and increased diabetes risk. We do not know how many different microbiota-derived molecules can stimulate β cell proliferation, but for the case of befA homologues, we know these to be sparsely distributed in bacterial genomes, such that microbiomes of low taxonomic diversity could lack these genes. The idea that microbiota-derived factors capable of protecting against diabetes are not widely conserved is consistent with human microbiota profiling data (Morgan et al., 2013), our own functional assays of bacteria in gnotobiotic zebrafish, and other gnotobiotic rodent experiments. For example, specific bacterial lineages have been shown to attenuate disease in diabetes models, including Segmented Filamentous Bacteria (SFB) in the non-obese diabetic (NOD) mouse (Kriegel et al., 2011; Yurkovetskiy et al., 2013) and Lactobacillus johnsonii in the Biobreeding rat model (Valladares et al., 2010). Furthermore, Wen and colleagues have shown that certain microbial assemblages, but not others, confer disease protection in neonatal NOD mice (Peng et al., 2014). Additional recent work by Wen and colleagues demonstrates early development as a critical window for microbiota modulation of disease risk in NOD mice (Hu et al., 2015). We have shown that BefA acts during early developmental stages in zebrafish, and we hypothesize that β cell expansion during this developmental window is important for disease prevention, and may be a critical period for clinical intervention for infants at risk for T1D development. Further work will be required to determine whether BefA is capable of inducing proliferation of adult β cells in zebrafish or other animals.

Why certain bacteria produce BefA is unclear. In the context of the zebrafish intestinal environment, BefA confers a slight colonization advantage to A. veronii, however this is unlikely to be related to its capacity to induce β cell mass, because the colonization requirement is only apparent in the context of co-colonization with wild type A. veronii that induce normal β cell numbers. It is possible that bacterial modulation of host β cell number serves a purpose for the bacteria not measured in our assay. Alternatively, bacteria may produce BefA for a purpose independent of β cell expansion and the host simply uses this bacterial molecule as a cue for its own developmental program. Learning the molecular basis for BefA sensing by the host, and whether it interacts directly or indirectly with β cells, will help shed light on the nature and evolutionary conservation of this interspecies signaling. It will also be important to understand the bacterial function of BefA in order to be able to manipulate its abundance for potential therapeutic purposes.

The incidence of autoimmune diseases such as T1D has been increasing markedly in developed nations over the past several decades. One theory to explain this phenomenon is the disappearing microbiota hypothesis, which proposes that over time, as our modern lifestyles have become increasingly sterile, we have lost ancestral microbial symbionts important for specific aspects of our health (Blaser and Falkow, 2009). Our discovery of a specific class of bacterial proteins that promote β cell expansion in early development is consistent with the hypothesis that loss of specific microbial taxa from gut microbiota could underlie increased diabetes risk. Specifically, we suggest that BefA-like proteins promote the establishment of a robust β cell population that is more resilient to subsequent β cell loss. Because befA is a relatively rare component of the microbiome, we cannot measure it directly from available metagenomic sequence data to test our hypothesis that befA abundance correlates with reduced diabetes risk. The low abundance of befA in metagenomes also highlights the challenge of discovering disease determinants from metagenomic data, and emphasizes the importance of functional screening approaches. The larval zebrafish has served as a valuable high-throughput vertebrate model for the identification of new compounds and pathways that can increase β cell numbers exogenously (Andersson et al., 2012; Wang et al., 2015a). We have employed the gnotobiotic zebrafish to explore how microbial cues modulate β cell development. Our discovery of BefA highlights the importance of the microbiota in shaping the development of an extra-intestinal tissue and influencing the overall metabolic state of the host. We postulate that resident bacteria are a rich and underexplored source of functionally conserved molecules that shape early host development in ways that impact disease risk in later life.

Material and methods

Gnotobiotic zebrafish

All zebrafish experiments were performed using protocols approved by the University of Oregon Institutional Care and Use Committee and followed standard protocols. Zebrafish embryos were derived germ-free (GF) as previously described (Bates et al., 2006). XGF and mono-associated larvae were also generated as previously described (Bates et al., 2006), except that all bacterial inoculate were added to GF flasks at 4 dpf at a final concentration of 106 CFUs/mL. In experiments quantifying the colonization levels of bacterial isolates, each strain was added to the embryo media (EM) and incubated with the larvae for 48 hr at 27°C. Larvae were sacrificed at 6 dpf, immediately before the gut was removed and homogenized in a small sample of sterile EM. Dilutions of this gut slurry were plated onto tryptic soy agar and allowed to incubate overnight at 30°C. Colonies from each gut were quantified. A minimum of 10 guts per mono-association or di-association were analyzed.

Free glucose assay

To measure β cell function in GF and CV zebrafish larvae, levels of free glucose were measured at 6 dpf using a free glucose assay kit (BioVision, Milpitas, CA) as described previously (Andersson et al., 2012; Gut et al., 2013) except that only 10 larvae were combined per tube. Three to five biological replicates (sets of 10 larvae) were completed for both GF and CV treatments each time the assay was conducted. Data shown here were combined from 3 separate experimental assays or technical replicates.

Cell free supernatant

GF fish were inoculated with secreted bacterial products at 4 dpf by adding cell free supernatant (CFS) at a final concentration of 500 ng/mL to the water of the sterile flasks. CFS was harvested from a 50 mL overnight culture of the specified bacterial strain. The cultures were centrifuged at 7000 g for 10 min at 4°C. The supernatant was then filtered through a 0.22-µm sterile tube top filter (Corning Inc., Corning, NY); sterile supernatant was concentrated at 4°C for 1 hr at 3000 g with a centrifugal device that has a 10 kDa weight cut off (Pall Life Sciences, Port Washington, NY).

For experiments utilizing proteinase K (Qiagen, Hilden, Germany), the enzyme was added to samples of CFS at a final concentration of 100 μg/mL and allowed to incubate at 55°C for 1 hr before inactivating the enzyme at 90°C for 10 min.

Ammonium sulfate fractionation

Ammonium sulfate fractionation was performed on un-concentrated, sterile CFS from a 50 mL overnight culture by slowly adding 100% ammonium sulfate until solutions of 20%, 40%, 60% and 80% ammonium sulfate were achieved. These solutions were prepared at 4°C. Precipitated proteins were collected from each fraction by centrifugation at 4°C and 14,000 g for 15 min. The proteins were resuspended in cold EM and dialyzed for 2–3 hr at 4°C before adding them to 4 dpf GF larvae at a final concentration of 500 ng/mL.

Mass spectrometry

The 60–80% ammonium sulfate fraction of the A. veroniiΔT2SS CFS was sent to the Proteomics Lab at Oregon Health and Science University in Portland, OR for protein identification (partial sequencing) analysis.

Protein expression and purification

The nucleotide sequence for the befA gene from was amplified from A. veronii using the following forward and reverse PCR primers respectively: 5’-GCCCATATGatgaacaagcgtaactggttgctg-3’ and 5’-GGCCTCGAGgcggctcgtttcagtcaagtc-3’. The nucleotide sequences for both the Enterococcus gallinarum and Enterobacter aerogenes befA gene homologs were obtained from NCBI and subsequently synthesized by GenScript, Piscataway, NJ. Each of these two genes was then cloned separately into the pET-21b plasmid (Novagen, Darmstadt, Germany), which contains an IPTG inducible promoter. A His•Tag was added to the C-terminal of the original BefA protein sequence for subsequent purification. As a control, a second version was also constructed lacking the tag. These vectors were then transformed into BL21 Escheria coli (RRID:WB_HT115(DE3)), treated with 0.5 – 1.0 mM IPTG during exponential growth phase (OD600 = 0.4–0.6) and allowed to grow for 3–4 more hours at 30°C. This resulted in both a CFS and cell lysate dominated by our proteins of interest, as confirmed via SDS-page gel electrophoresis by the presence of dark bands of the expected sizes for each protein. These bands were absent from BL21 cultures carrying an empty pET-21b vector. The CFS from these inductions was added to GF zebrafish at 4 dpf at a final concentration of 500 ng/mL.

For purification of BefA, IPTG induced BL21 cells were sonicated at 32,000 g in a 50 nM Tris, 150 mM NaCl buffer (buffer A). The supernatant was then added to a solution of nickel beads (Thermo Scientific HisPur Ni-NTA Resin, Waltham, MA) to capture the His•tag. The beads were washed several times in a 30 mM imidazol solution in buffer A and subsequently eluted in 300 mM imidazole solution in buffer A. The isolation of pure BefA was confirmed with SDS-page gel electrophoresis by the presence of a single band of about 29 kDa in size. Purified BefA was added to 4 dpf GF fish at a final concentration of 500 ng/mL.

Experimental bacterial strains

To create the A. veroniiΔbefA mutant strain, a vector containing a chloramphenicol resistance cassette was transformed into SM10 E. coli. Conjugation between wild-type Aeromonas veronii HM21 and the vector carrying SM10 E. coli strain was carried out, allowing the chloramphenicol resistance gene to replace the befA locus in A. veronii via allelic exchange. Candidate mutants were selected for loss of the plasmid and maintenance of chloramphenicol resistance. Insertion of the chloramphenicol cassette into the befA locus was verified in these candidates by PCR. Joerg Graf graciously provided us with the A. veroniiΔT2SS strain (Maltz and Graf, 2011).

Primary islet cell type quantifications and EdU staining in larvae

Tg(-1.0insulin:eGFP) (RRID:ZFIN_ZDB-GENO-100513-10, ZIRC, Eugene, OR) (diIorio et al., 2002) zebrafish embryos were used to visualize and quantify the total number of β cells in developing larvae. Tg(insulin:PhiYFP-2a-nsfB, sst2:mCherry) (RRID:ZFIN_ZDB-GENO-120217-6) (Wang et al., 2015a) were obtained from Jeff Mumm and were used to visualize and quantify δ cells. All experiments were analyzed at 6 dpf unless otherwise specified. At all time points in all experiments, larvae were fixed with 4% paraformaldehyde supplemented with 0.01% TritonX-100 (Thermo Fisher Scientific, Waltham, MA) at 4°C overnight, or at room temperature for 2–3 hr, and then washed with PBS. The following antibodies were used to distinguish α and β cells: guinea-pig anti-insulin (Dako Cat# A0564, RRID:AB_10013624, Carpinteria, CA), mouse anti-glucagon (Sigma-Aldrich Cat# G2654, RRID:AB_259852), St. Louis, MO), rabbit anti-GFP (Molecular Probes Cat# A-11122, RRID:AB_221569), mouse anti-mCherry (Abcam Cat# ab125096, RRID:AB_11133266, Cambridge, MA), Alexa Fluor 488 goat anti-rabbit (Thermo Fisher Scientific, Waltham, MA), anti-mouse Cy3 (Jackson ImmunoResearch Laboratories Inc., West Grove, PA), Alexa Fluor 488 goat anti-guinea-pig (Thermo Fisher Scientific, Waltham, MA), and TO-PRO-3-Iodide (642/661) (Thermo Fisher Scientific, Waltham, MA).

For experiments quantifying proliferation, EdU was added at 4 dpf directly to the EM at a final concentration of 0.1 mg/mL. The Click-iT EdU Imaging Kit (Thermo Fisher Scientific, Waltham, MA) was used to process the EdU label in whole fixed zebrafish prior to antibody staining, according to the manufacturer’s protocols. Whole, antibody-stained larvae were mounted for confocal microscopy (BioRad Radiance 2100) with their right side facing up against the cover slip, which was flattened sufficiently to spread out the cells within the islet for optimal quantification of individual cells. For quantification of β cells and other primary islet cells, the entire endocrine portion of the pancreas was scanned using a 60X objective (Nikon Eclipse E600FN), and Fiji (RRID:SCR_002285) (Schindelin et al., 2012) software was used to analyze each image stack. For quantification of pancreatic exocrine tissue proliferation, Tg(ptf1a:eGFP) (RRID:ZFIN_ZDB-GENO-080111-1, ZIRC, Eugene, OR) (Thisse et al., 2004) zebrafish were scanned through the entire pancreas with a 20X objective (Nikon Eclipse E600FN) and Fiji was used to analyze the percentage of proliferative cells in single sections from the center of the organ. Images were prepared for publication using the open source Inkscape software (RRID:SCR_014479).

For experiments quantifying insulin-expressing cells in the region of the EPD, zebrafish were processed as described above, and analyzed on a Leica fluorescent microscope using a 2x objective.

BefA phylogenetic analysis

We screened for BefA homologs across microbial species using a blastp-based (Altschul et al., 1997) search of the UniProt Knowledgebase (UniProt Consortium, 2015) (version 6/2015); default search parameters were changed to allow (i) a maximum E-value of 1.0 and (ii) an arbitrarily large number of database hits. We classified database hits as 'close homologs' if amino acid sequence identity exceeded 50% (relative to the query length) and 'distant homologs' if their percent identity exceeded 20%. For phylogenetic analysis at the species level, each species was represented by the hit of highest percent identity to BefA among isolates of that species (if any); an analogous procedure was used for genus-level analysis. Aligned portions of database sequences were isolated and multiply aligned with MUSCLE (RRID:SCR_011812) (Edgar, 2004). Phylogenetic trees were constructed from these multiple sequence alignments using PhyML (RRID:SCR_014629) (Guindon and Gascuel, 2003) and visualized within the Phylogeny.fr webserver (Dereeper et al., 2008). Microbial genera were classified as 'human-associated' if they occurred with relative abundance >0.01% in at least 5 metagenomes from the Human Microbiome Project (Huttenhower et al., 2012) as profiled by MetaPhlAn (RRID:SCR_004915) (Segata et al., 2012). Secretion signal peptides were predicted from amino acid sequences using SignalP (Petersen et al., 2011).

Statistical analysis

Appropriate sample sizes for all experiments were estimated a priori using a power of 80% and a significance level of 0.05. From preliminary experiments we estimated variance and effect. For larval β cell quantification, these parameters suggested using a sample size of 30 in order to detect significant changes between treatment groups. Therefore, each experiment contained about 10–15 biological replicates or individual fish per treatment group, although some larger experiments had fewer biological replicates due to limited material. Entire experiments or technical replicates were repeated multiple times, resulting in pooled data sets of about 20–50 biological replicates. These data are represented in the figures as box and whisker plots, which display the data median (line within the box), first and third quartiles (top and bottom of the box), and 95% confidence interval (whiskers). Any data point falling outside the 95% confidence interval is represented as a solid dot. These pooled data were analyzed through the statistical software RStudio. For experiments comparing just two differentially treated populations, a Student’s t-test with equal variance assumptions was used. For experiments measuring a single variable with multiple treatment groups, a single factor ANOVA with post hoc means testing (Tukey) was utilized. A p-value of less than 0.05 was required to reject the null hypothesis that no difference existed between groups of data.

References

  1. 1
  2. 2
  3. 3
  4. 4
  5. 5
  6. 6
  7. 7
  8. 8
  9. 9
  10. 10
  11. 11
  12. 12
    Dysbiosis in gastrointestinal disorders
    1. C Chang
    2. H Lin
    (2016)
    Best Practice & Research Clinical Gastroenterology 30:3–15.
    https://doi.org/10.1016/j.bpg.2016.02.001
  13. 13
  14. 14
  15. 15
  16. 16
  17. 17
  18. 18
  19. 19
  20. 20
  21. 21
  22. 22
  23. 23
  24. 24
  25. 25
  26. 26
  27. 27
    Symbiosis of Aeromonas veronii biovar sobria and Hirudo medicinalis, the medicinal leech: a novel model for digestive tract associations
    1. J Graf
    (1999)
    Infection and Immunity 67:1–7.
  28. 28
  29. 29
  30. 30
  31. 31
  32. 32
  33. 33
  34. 34
  35. 35
    Structure, function and diversity of the healthy human microbiome
    1. C Huttenhower
    2. D Gevers
    3. R Knight
    4. S Abubucker
    5. JH Badger
    6. AT Chinwalla
    7. HH Creasy
    8. AM Earl
    9. MG FitzGerald
    10. RS Fulton
    11. MG Giglio
    12. K Hallsworth-Pepin
    13. EA Lobos
    14. R Madupu
    15. V Magrini
    16. JC Martin
    17. M Mitreva
    18. DM Muzny
    19. EJ Sodergren
    20. J Versalovic
    21. AM Wollam
    22. KC Worley
    23. JR Wortman
    24. SK Young
    25. Q Zeng
    26. KM Aagaard
    27. OO Abolude
    28. E Allen-Vercoe
    29. EJ Alm
    30. L Alvarado
    31. GL Andersen
    32. S Anderson
    33. E Appelbaum
    34. HM Arachchi
    35. G Armitage
    36. CA Arze
    37. T Ayvaz
    38. CC Baker
    39. L Begg
    40. T Belachew
    41. V Bhonagiri
    42. M Bihan
    43. MJ Blaser
    44. T Bloom
    45. V Bonazzi
    46. J Paul Brooks
    47. GA Buck
    48. CJ Buhay
    49. DA Busam
    50. JL Campbell
    51. SR Canon
    52. BL Cantarel
    53. PSG Chain
    54. I-MA Chen
    55. L Chen
    56. S Chhibba
    57. K Chu
    58. DM Ciulla
    59. JC Clemente
    60. SW Clifton
    61. S Conlan
    62. J Crabtree
    63. MA Cutting
    64. NJ Davidovics
    65. CC Davis
    66. TZ DeSantis
    67. C Deal
    68. KD Delehaunty
    69. FE Dewhirst
    70. E Deych
    71. Y Ding
    72. DJ Dooling
    73. SP Dugan
    74. W Michael Dunne
    75. A Scott Durkin
    76. RC Edgar
    77. RL Erlich
    78. CN Farmer
    79. RM Farrell
    80. K Faust
    81. M Feldgarden
    82. VM Felix
    83. S Fisher
    84. AA Fodor
    85. LJ Forney
    86. L Foster
    87. V Di Francesco
    88. J Friedman
    89. DC Friedrich
    90. CC Fronick
    91. LL Fulton
    92. H Gao
    93. N Garcia
    94. G Giannoukos
    95. C Giblin
    96. MY Giovanni
    97. JM Goldberg
    98. J Goll
    99. A Gonzalez
    100. A Griggs
    101. S Gujja
    102. S Kinder Haake
    103. BJ Haas
    104. HA Hamilton
    105. EL Harris
    106. TA Hepburn
    107. B Herter
    108. DE Hoffmann
    109. ME Holder
    110. C Howarth
    111. KH Huang
    112. SM Huse
    113. J Izard
    114. JK Jansson
    115. H Jiang
    116. C Jordan
    117. V Joshi
    118. JA Katancik
    119. WA Keitel
    120. ST Kelley
    121. C Kells
    122. NB King
    123. D Knights
    124. HH Kong
    125. O Koren
    126. S Koren
    127. KC Kota
    128. CL Kovar
    129. NC Kyrpides
    130. PS La Rosa
    131. SL Lee
    132. KP Lemon
    133. N Lennon
    134. CM Lewis
    135. L Lewis
    136. RE Ley
    137. K Li
    138. K Liolios
    139. B Liu
    140. Y Liu
    141. C-C Lo
    142. CA Lozupone
    143. R Dwayne Lunsford
    144. T Madden
    145. AA Mahurkar
    146. PJ Mannon
    147. ER Mardis
    148. VM Markowitz
    149. K Mavromatis
    150. JM McCorrison
    151. D McDonald
    152. J McEwen
    153. AL McGuire
    154. P McInnes
    155. T Mehta
    156. KA Mihindukulasuriya
    157. JR Miller
    158. PJ Minx
    159. I Newsham
    160. C Nusbaum
    161. M O’Laughlin
    162. J Orvis
    163. I Pagani
    164. K Palaniappan
    165. SM Patel
    166. M Pearson
    167. J Peterson
    168. M Podar
    169. C Pohl
    170. KS Pollard
    171. M Pop
    172. ME Priest
    173. LM Proctor
    174. X Qin
    175. J Raes
    176. J Ravel
    177. JG Reid
    178. M Rho
    179. R Rhodes
    180. KP Riehle
    181. MC Rivera
    182. B Rodriguez-Mueller
    183. Y-H Rogers
    184. MC Ross
    185. C Russ
    186. RK Sanka
    187. P Sankar
    188. J Fah Sathirapongsasuti
    189. JA Schloss
    190. PD Schloss
    191. TM Schmidt
    192. M Scholz
    193. L Schriml
    194. AM Schubert
    195. N Segata
    196. JA Segre
    197. WD Shannon
    198. RR Sharp
    199. TJ Sharpton
    200. N Shenoy
    201. NU Sheth
    202. GA Simone
    203. I Singh
    204. CS Smillie
    205. JD Sobel
    206. DD Sommer
    207. P Spicer
    208. GG Sutton
    209. SM Sykes
    210. DG Tabbaa
    211. M Thiagarajan
    212. CM Tomlinson
    213. M Torralba
    214. TJ Treangen
    215. RM Truty
    216. TA Vishnivetskaya
    217. J Walker
    218. L Wang
    219. Z Wang
    220. DV Ward
    221. W Warren
    222. MA Watson
    223. C Wellington
    224. KA Wetterstrand
    225. JR White
    226. K Wilczek-Boney
    227. Y Wu
    228. KM Wylie
    229. T Wylie
    230. C Yandava
    231. L Ye
    232. Y Ye
    233. S Yooseph
    234. BP Youmans
    235. L Zhang
    236. Y Zhou
    237. Y Zhu
    238. L Zoloth
    239. JD Zucker
    240. BW Birren
    241. RA Gibbs
    242. SK Highlander
    243. BA Methé
    244. KE Nelson
    245. JF Petrosino
    246. GM Weinstock
    247. RK Wilson
    248. O White
    249. Human Microbiome Project Consortium
    (2012)
    Nature 486:207–214.
    https://doi.org/10.1038/nature11234
  36. 36
  37. 37
  38. 38
  39. 39
  40. 40
  41. 41
  42. 42
  43. 43
  44. 44
  45. 45
  46. 46
  47. 47
  48. 48
  49. 49
  50. 50
  51. 51
  52. 52
  53. 53
  54. 54
  55. 55
  56. 56
  57. 57
  58. 58
  59. 59
  60. 60
  61. 61
  62. 62
  63. 63
  64. 64
  65. 65
  66. 66
  67. 67
  68. 68
  69. 69
  70. 70
  71. 71
  72. 72

Decision letter

  1. Didier Y R Stainier
    Reviewing Editor; Max Planck Institute for Heart and Lung Research, Germany

In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.

Thank you for submitting your article "A conserved bacterial protein induces pancreatic beta cell expansion during zebrafish development" for consideration by eLife. Your article has been favorably evaluated by Marianne Bronner as the Senior Editor and three reviewers, including Daniel Hesselson (Reviewer #2) and a member of our Board of Reviewing Editors.

The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.

Summary:

In this elegant paper from the Guillemin lab, the authors report the surprising discovery of a bacterial protein that positively regulates the number of pancreatic β cells in zebrafish by inducing β cell proliferation. The data are in general of high quality (see below for some concerns) as is the writing, and the significance of the findings is very high.

Essential revisions:

1) The variability of the wild-type sample in Figure 1 (especially at 6 dpf) is really quite surprising (previous publications have not reported such variability) and could be due to a number of factors including the method used to count (transgene vs. Ab), the imaging platform, the counting method itself (by hand vs using a particular software), whether the fish are well synchronized and healthy (and thus with a similar amount of yolk left), and whether they are fed or not(the variability in control animals in the later figures appears to be reduced).

[Ideally, the combination of a membrane and a nuclear marker would be best.]

It is also possible that the variability in Figure 1 comes from non-random sampling of the cohort. If the authors started with a pool of animals at 3 dpf and chose the healthiest/most developed larvae each day then collected the remainder at 6 dpf that would explain the apparent loss of β cells and increased variability between 5 dpf and 6 dpf in the CV controls. Perhaps this experiment should be repeated with a defined sampling plan.

2) In the germ free animals (e.g., Figure 1B), what is the identity of the other cells in the islet? The authors should stain the animals for glucagon and somatostatin to determine A) whether CV vs GF conditions affect differentiation of the other endocrine cell types and B) if α/δ cell mass is affected whether this activity is also mediated by BefA. It is also possible that BefA affects hormone expression rather than cell differentiation which would also be an interesting observation.

3) Does the microbiota/BefA act solely by stimulating proliferation of mature β cells or does it also have effects on endocrine progenitors? In mammals these processes are temporally separated with the vast majority of neogenesis occurring before birth under sterile conditions followed by diet (and now possibly microbiota) induced proliferation postnatally. In zebrafish both of these processes occur concurrently during early larval development and contribute to the expansion of the β cell mass from 3 to 6 dpf. The 48 h EdU pulse used in Figure 5 would have labeled proliferating endocrine progenitor cells that subsequently differentiated into β cells in addition to proliferating β cells. Phospho-histone H3 staining should be used to clarify whether BefA stimulates proliferation of mature β cells.

https://doi.org/10.7554/eLife.20145.023

Author response

[…] Essential revisions:

1) The variability of the wild-type sample in Figure 1 (especially at 6 dpf) is really quite surprising (previous publications have not reported such variability) and could be due to a number of factors including the method used to count (transgene vs. Ab), the imaging platform, the counting method itself (by hand vs using a particular software), whether the fish are well synchronized and healthy (and thus with a similar amount of yolk left), and whether they are fed or not (the variability in control animals in the later figures appears to be reduced).

[Ideally, the combination of a membrane and a nuclear marker would be best.]

It is also possible that the variability in Figure 1 comes from non-random sampling of the cohort. If the authors started with a pool of animals at 3 dpf and chose the healthiest/most developed larvae each day then collected the remainder at 6 dpf that would explain the apparent loss of β cells and increased variability between 5 dpf and 6 dpf in the CV controls. Perhaps this experiment should be repeated with a defined sampling plan.

We are confident that the fish used for Figure 1A were controlled with respect to health, developmental synchronization, feeding (all fish are starved for GF derivation protocols), and experimental method. The data were also randomly sampled, as we raised a separate pool of fish for each day samples were taken.

The data in our manuscript is displayed using box plots, with whiskers denoting 5-95% of the values of the entire data set, whereas the data from other manuscripts quantifying β cell number is most often displayed using a histogram, with error bars representing either the standard error (SE) or the standard error of the mean (SEM) (Andersson et al. Cell Met. 2012, Hesselson et al. PNAS. 2009, Wang et al. eLife. 2015, Ye et al. Development. 2015). We chose to use box plots to represent our data to convey maximal information and transparency, as discussed by Krzywinski and Altman in Nature Methods in 2014. It is difficult to get an indication of the total variation of a data set from SE and SEM. These measurements are usually much smaller than the total variation we have shown in our box plots. Therefore, it is understandable that our figures appear to have larger than normal variation. For comparison, we transformed Figure 1A into a histogram with error bars representing SEM, and placed it next to a similar graph (Figure 3, panel G) from Andersson et al. Cell Met. 2012, which also used SEM. Of note, the error bars from both works are similar in size. The exact numerical values for the standard deviation as well as the SEM for all of our data can also be found in the corresponding source data file. Descriptions of our graphical representations are also indicated in our manuscript within the figure legends and methods.

2) In the germ free animals (e.g., Figure 1B), what is the identity of the other cells in the islet? The authors should stain the animals for glucagon and somatostatin to determine A) whether CV vs GF conditions affect differentiation of the other endocrine cell types and B) if α/δ cell mass is affected whether this activity is also mediated by BefA. It is also possible that BefA affects hormone expression rather than cell differentiation which would also be an interesting observation.

We have added quantifications of both glucagon-expressing α and somatostatin-expressing δ cells in GF, CV and BefA treated larvae to Figure 5. Text descriptions of these results start in the last paragraph of the subsection “BefA facilitates β cell expansion by inducing proliferation”. We found no change in the total numbers of either of these cell types across treatments.

3) Does the microbiota/BefA act solely by stimulating proliferation of mature β cells or does it also have effects on endocrine progenitors? In mammals these processes are temporally separated with the vast majority of neogenesis occurring before birth under sterile conditions followed by diet (and now possibly microbiota) induced proliferation postnatally. In zebrafish both of these processes occur concurrently during early larval development and contribute to the expansion of the β cell mass from 3 to 6 dpf. The 48 h EdU pulse used in Figure 5 would have labeled proliferating endocrine progenitor cells that subsequently differentiated into β cells in addition to proliferating β cells. Phospho-histone H3 staining should be used to clarify whether BefA stimulates proliferation of mature β cells.

As suggested, we performed anti-PH3 antibody staining on 6 dpf larvae to determine whether BefA or other microbes could increase proliferation of mature β cells. We found PH3 marked β cells to be very rare in all treatments. Based on the initial results from our experiment of about 50 fish per treatment, we calculated the number of samples needed to adequately power a study that wouldbe able to detect a difference in PH3 staining between our treatments. We approximated that 1300 fish per treatment were required to be confident that no difference actually exists in PH3 staining. We are unable to perform an experiment of this scale and therefore we cannot support or refute the possibility that BefA acts by inducing proliferation of existing β cells.

In parallel, we measured the effect of the microbiota on neogenesis. Since differentiation of β cells from 4-6 dpf occurs within the extra-pancreatic duct (EPD) (Dong et al., 2007; Hesselson et al., 2009), we quantified the number of insulin positive cells located in the region of the EPD in 6 dpf GF and CV larvae. Again, detection of β cells in this region at a single moment in time is a relatively rare event, but we were able to use larvae from past experiments to gather sufficient sample sizes to power this study. We found that there were significantly more EPD localized insulin+ cells in CV fish than in GF fish. We have included a supplemental figure (Figure 5—figure supplement 1) with these data. We also quantified EPD localized β cells in BefA treated fish, however we did not have a sufficient sample size (which we calculated to be about 200 fish) to detect a statistically significant difference, although the result is trending in the direction of CV larvae.

Unfortunately, we cannot conclusively determine the mechanism by which BefA increases EdU incorporation in β cells. We have changed the text to reflect the possibilities discussed above, in the subsection “BefA facilitates β cell expansion by inducing proliferation”. We believe that this question is important, and should be more extensively explored in a future publication focused on how the host detects and interprets BefA.

https://doi.org/10.7554/eLife.20145.024

Article and author information

Author details

  1. Jennifer Hampton Hill

    1. Institute of Molecular Biology, University of Oregon, Eugene, United States
    Contribution
    JHH, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  2. Eric A Franzosa

    1. Biostatistics Department, Harvard T. H. Chan School of Public Health, Boston, United States
    2. The Broad Institute, Cambridge, United States
    Contribution
    EAF, Acquisition of data, Analysis and interpretation of data
    Competing interests
    The authors declare that no competing interests exist.
  3. Curtis Huttenhower

    1. Biostatistics Department, Harvard T. H. Chan School of Public Health, Boston, United States
    2. The Broad Institute, Cambridge, United States
    Contribution
    CH, Acquisition of data, Analysis and interpretation of data
    Competing interests
    The authors declare that no competing interests exist.
  4. Karen Guillemin

    1. Institute of Molecular Biology, University of Oregon, Eugene, United States
    2. Humans and the Microbiome Program, Canadian Institute for Advanced Research, Toronto, Canada
    Contribution
    KG, Conception and design, Analysis and interpretation of data, Drafting or revising the article
    For correspondence
    1. kguillem@uoregon.edu
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon 0000-0001-6004-9955

Funding

National Institutes of Health (P50GM098911)

  • Karen Guillemin

National Institutes of Health (T32 GM007413-37)

  • Jennifer Hampton Hill

National Institutes of Health (P01HD22486)

  • Karen Guillemin

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank J Graf for providing us with the Aeromonas veronii HM21strains, T Wiles for guidance in the construction of the Aeromonas veronii HM21ΔbefA strain, A Rolig and A Banse for technical advice with CFS preparation, E Goers Sweeny for technical advice with protein purification, P De Verteuil for assistance with colonization experiments, R Sockol and the University of Oregon Zebrafish Facility for maintenance of zebrafish lines, and JPostlethwait, J Eisen, A Powell, A Rolig and T Jones for critical reading of the manuscript. Research reported in this publication was supported by the National Institute of General Medical Sciences of the National Institutes of Health under award numbers P50GM098911 and T32 GM007413-37. Grant P01HD22486 provided support for the Oregon Zebrafish Facility. The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH.

Ethics

Animal experimentation: This study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. All zebrafish experiments were performed using protocols approved by the University of Oregon Institutional Care and Use Committee under protocol number 15-83 and institutional Animal Welfare Assurance number A-3009-01.

Reviewing Editor

  1. Didier Y R Stainier, Reviewing Editor, Max Planck Institute for Heart and Lung Research, Germany

Publication history

  1. Received: July 28, 2016
  2. Accepted: November 3, 2016
  3. Version of Record published: December 13, 2016 (version 1)

Copyright

© 2016, Hill et al

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 2,582
    Page views
  • 630
    Downloads
  • 3
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, Scopus, PubMed Central.

Comments

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Download citations (links to download the citations from this article in formats compatible with various reference manager tools)

Open citations (links to open the citations from this article in various online reference manager services)

  1. Further reading

Further reading

    1. Developmental Biology and Stem Cells
    Tang Hai et al
    Research Article
    1. Biochemistry
    2. Neuroscience
    Hansen Wang et al.
    Research Article