Use of Genetically Encoded Calcium Indicators (GECIs) Combined with Advanced Motion Tracking Techniques to Examine the Behavior of Neurons and Glia in the Enteric Nervous System of the Intact Murine Colon

Introduction

The function of neural networks, and their role in the generation of complex patterns of motor behavior in the gut is inherently difficult to study because of the large number of interacting neurons that have emergent properties that can’t readily be predicted from recording from individual neurons (Bayguinov et al., 2010a,b,c).

The enteric nervous system (ENS) in the bowel consists of an intricate system of layered and interconnected ganglionated neural networks that regulate motility, secretion and blood flow (Lomax and Furness, 2000; Furness, 2012; Okamoto et al., 2014; Smith et al., 2015). The largest and best known layer of neurons, the myenteric plexus, located between the longitudinal and circular muscle layers, is a two-dimensional sheet of interconnected ganglia that generates and coordinates many aspects of motor behavior of the colon (Bayguinov et al., 2010a,b). In the large bowel, most myenteric ganglia contain many different functional classes of neurons: sensory neurons, several ascending and descending interneurons, as well as excitatory and inhibitory motor neurons (Lomax and Furness, 2000; Smith et al., 2007).

Myenteric neurons are surrounded and embedded in a network of glial cells and their processes, which were largely thought to be the “glue” that held ganglia together (Gulbransen and Sharkey, 2012). Glia cells respond to a variety of neurotransmitter agonists (e.g., ACh, ATP, 5-HT, NK2, substance P) by releasing intracellular stores of Ca²⁺, and appear to be innervated or in close proximity to nerve varicosities (e.g., including seroternergic and sympathetic nerves: Broadhead et al., 2012; Gulbransen and Sharkey, 2009, 2012; Boesmans et al., 2013a; Okamoto et al., 2014). Furthermore, glial cells can display prolonged Ca²⁺ responses compared to myenteric neurons, peaking some ~40 s following neural reflex stimulation (Broadhead et al., 2012). This is more prolonged than glial Ca²⁺ responses observed after addition of purinergic, 5-HT and ACh agonists (Gulbransen and Sharkey, 2009; Boesmans et al., 2013a; McClain et al., 2014). Ca²⁺ waves in glial processes, which surround neurons, appear to summate in their soma where they possibly release substances such as prostaglandins, nitric oxide or purines that either help terminate or enhance the excitability of specific neurons (MacEachern et al., 2011; Broadhead et al., 2012; Heredia et al., 2012; Smith et al., 2015). These observations support the notion that enteric glia, similar to those in the peripheral and central nervous systems, may not just sense, but also modulate neurotransmission (Robitaille, 1998).

Using genetically-encoded calcium indicators such as GCaMP3 (Tian et al., 2009; Wilms and Häusser, 2009; Akerboom et al., 2012; Yamada and Mikoshiba, 2012; Zariwala et al., 2012; Sun et al., 2013) expressed in various enteric neurons and glia (Boesmans et al., 2015; Foong et al., 2015), we show here that the behavior of these cells and their role in colonic motor patterns can be readily analyzed in the undissected isolated colon using low-power imaging of GCaMP3. The specificity and quality of labeling of neurons and glia with GCaMP3 allows a much better understanding of how motor behaviors emerge from the activity of individual neurons in the gut both in situ and in vivo.

Methods

Generation of Specific GCaMP3 Mice

All mice (C57BL/6 background of either sex) were housed in a transgenic animal facility under pathogen-free conditions on a 12 h light/dark cycle with food and water ad libitum. The following mice were generated for this study: Wnt1-Cre; Rosa26-Lox-STOP-Lox-GCaMP3 mice (GCaMP3 targeted to neural crest derivatives, including enteric neurons and glia via Wnt1-Cre transgenic mice; Jax #009107; heretofore referred to as Wnt1-GCaMP3), GFAP-Cre; Rosa26-Lox-STOP-Lox-GCaMP3 (GCaMP3 targeted to glia via mice expressing Cre from the human GFAP promoter; Jax #004600; heretofore referred to as GFAP-GCaMP3), ChAT-Cre; Rosa26-Lox-STOP-Lox-GCaMP3 (GCaMP3 targeted to cholinergic neurons via ChAT-Cre BAC transgenic mice; MMRRC #37336; heretofore referred to as ChAT-GCaMP3) and nNOS-CreER; Rosa26-Lox-STOP-Lox-GCaMP3 (GCaMP3 targeted to nitrergic neurons via nNOS-CreER transgenic mice; Jax #014541; heretofore referred to as nNOS-GCaMP3). Mice expressing a conditional (i.e., Cre-dependent) allele of GCaMP3 in the Rosa26 locus (Rosa26-Lox-STOP-Lox-GCaMP3) were described previously (Jax # 014538; Zariwala et al., 2012). nNOS-GCaMP3 mice were given intraperitoneal injections of tamoxifen (TM) between the ages of 4–6 weeks in order to activate CreER and hence GCaMP3 expression. TM was first dissolved in 100% ethanol to a concentration of 100 μg/μl by vortexing for 15 min. Then this TM stock solution was further dissolved 1:5 in safflower oil, vortexed for 20 min and sonicated for 30 min immediately prior to injection. 100 μl of this working TM solution (2 mg) was injected intraperitoneally for three consecutive days, and animals were sacrificed 5–10 days after the last injection.

These procedures were in accordance with National Institutes of Health guidelines for the care and use of laboratory animals and approved by the Animal Ethics Committee at the University of Nevada, Reno.

Tissue Preparation

Mice were euthanized by inhalation of a 5% concentration of isoflurane, followed by cervical dislocation. A ventral midline incision was made and the whole colon (proximal colon and distal colon) was carefully excised. The entire large intestine was removed from the transgenic mice and the mesentery trimmed and pinned to the bottom of a 30 mL Sylgard-lined (Dow Corning Corp., Midland, MI, USA) organ bath in 4°C oxygenated Krebs Ringer Buffer solution (KRB, concentrations in mM: NaCl: 120.35; KCl: 5.9; NaHCO₃: 15.5; MgCl₂: 1.2; NaH₂PO₄: 1.2; glucose: 11.5; CaCl₂: 2.5; pH 7.4) oxygenated with 97% O₂ and 3% CO₂. The ends of the colon (~10 length) were cut open along the mesenteric border and pinned flat, mucosal side uppermost, to allow for brush stimulation of the mucosa to evoke neural reflexes that elicit colonic migrating motor complexes (CMMCs; Bayguinov et al., 2010a,b; Heredia et al., 2009). In some preparations the ends of the colon were attached to a frog heart clip to allow the colon to be briefly elongated in a localized pinned region by up to 10% of its length in order to evoke polarized descending inhibitory and excitatory reflexes triggered by colonic elongation as shown in Figure 1. Oral elongation (longitudinal stretch at the oral end of the preparation) inhibits CMMCs, whereas anal elongation (longitudinal stretch at the anal end of the preparation) generates a single CMMC (Heredia et al., 2010). The preparation was allowed to equilibrate for ~1 h as the bath gradually reached the recording temperature of 36–37°C. Most of our Ca²⁺ imaging was mainly restricted to the mid colon region approximately in the middle of the whole colon (~40 long) so that we could compare activity in GCaMP3 cells with that obtained previously with Ca²⁺ dyes (Bayguinov et al., 2010a,b, 2012; Heredia et al., 2010; Broadhead et al., 2012).

FIGURE 1

Figure 1. Schematic of the experimental setup used to record Ca²⁺ activity in targeted cells in the isolated intact mouse colon. The ends of the preparation were slit and pinned open to allow for mucosal stimulation using a brush and/or longitudinal stretch stimulation using clips. In the middle of the preparation where recordings were taken, the colon was stabilized by pinning through the outer edges.

Image Acquisition

Ca²⁺ imaging recordings were performed on an Olympus BX51WI upright fluorescence microscope using 10 and 20× Fluor water immersion lenses (Olympus, Center Valley, PA, USA). The tissue was excited at 488 nm, using a X-Cite series 120Q (Lumen Dynamics, Ontario, Canada) and a modified GFP dichroic cube (excitation 488 nm; emission 543 nm; Chroma Technology Corp., Bellows Falls, VT, USA). Movies were captured with an Andor iXon +897 EMCCD camera (normal recording ~2000 frames @ 32.4 Hz) using Andor Solis 4.14 software (Andor Technology, Belfast, UK). Solis (16-bit, tif) files were analyzed on a MacPro desktop computer (Apple Inc., Cupertino, CA, USA) using in-house analysis software (Volumetry G8c, G. W. Hennig).

Drugs Used

Hexamethonium bromide was purchased from Sigma-Aldrich (St. Louis, MO, USA).

Results

Ca²⁺-induced fluorescence in enteric neurons and glia expressing GCaMP3 was significantly brighter with better signal-to-noise ratio compared to Ca²⁺ dyes (e.g., Fluo-4, see Bayguinov et al., 2010a,b; Okamoto et al., 2012) such that no dissection of the tissue was needed. Labeled cells were observed throughout the entire gut. By focusing through the serosa and longitudinal muscle, we could easily visualize the myenteric plexus without confocal imaging. This allowed, in tandem with better cameras, imaging of much greater fields of view while still being able to resolve individual cells and their Ca²⁺ transients (10×). At low magnifications, Ca²⁺ transients in cell bodies could be measured, however, activity in fibers and varicosities could not be resolved at low power but could be resolved at higher power (Bayguinov et al., 2012). Ca²⁺ transients in cells in the myenteric plexus or submucous plexus (after focusing even deeper though the underlying circular muscle) were easily resolved without the need to dissect of any of the outer layers of the colon. The results presented below are produced from a representative subset of the animals tested during the development of the transgenic mice variants, and illustrate the Ca²⁺ responses in each of the targeted cells.

Modified Image Analysis Suitable for GCaMP3 Expressing Cells

An obvious benefit of using GCaMP3 expressed in cells is that imaging can be performed without dissection, thus better preserving circumferential pathways and motor responses during a CMMC. However, without dissecting and stretching preparations flat (Bayguinov et al., 2010a), we were less physically able to restrain movements. Similarly, due to the absence of non-specific contamination of the surface of the preparation common with Ca²⁺ dyes, we were unable to utilize motion tracking routines that performed frame-by-frame comparisons using the pattern of static fluorescence. In this section, we demonstrate modified motion tracking routines, suitable for structures that display dynamic changes in intensity and demonstrate distortion correction rather than displacement correction to compensate for non-uniform contraction and relaxation.

Automated Nuclei Tracking

The dark nuclei in cells expressing GCaMP3 were used as reference points since they had a relatively uniform size and shape. Briefly, a quiescent period (100–200 ms) in the recording was averaged to reduce noise (Figure 2A), then a threshold was ramped through the entire intensity range and particles generated for those areas matching a criteria of size and circularity of a nucleus (area = 25–50 μm², radius_perimeter:radius_area = 0.95–1.4). These “nuclear” particles were saved as a particle mask used for tracking (Figure 2B) and a separate particle mask was derived from the “nuclear” particles by performing a dilate routine (2 pixels) and removing the original “nuclear” particles. These peri-nuclear “rings” were used to extract changes in Ca²⁺ intensity from the cytoplasm of cells at a consistent radius around the nucleus (Figure 2C). Any manual corrections were done on the particle masks. ROIs were automatically created using bounding rectangles of individual particles and the ROIs and the “nuclear” particles were transferred to the original recording. Instead of using the absolute intensity value of the nuclei in the movie which alters during activity in a cell, the differential intensity between the nucleus and cytoplasm (at points ±2 pixels either side of the edge of the bitmask) was calculated every 30° (12 radial sampling points per nuclei). The centroid point within a search range that returned the maximum sum of differential intensities was stored and represents the most centralized position of the nucleus. This tracking was extremely fast and could process 15 neurons in a 2000 frame movie (±5 pixel search grid) in less than a second. The trajectories of each nuclei were smoothed (±60 ms see Figure 2D).

FIGURE 2

Figure 2. Advanced motion tracking and analysis of cells labeled with GCaMP3 in myenteric ganglia. The lack of non-fluorescing landmarks in tissues prevented the use of tracking routines based on the shape and intensity pattern of a reference region. Instead, to stabilize motion associated with CMMCs, the relatively darker nuclei (vs. cytoplasm) of these cells were used to generate ovoid “nuclear” particles. (A) Fluorescence image of a myenteric ganglion from a Wnt1-Cre-GCaMP3 mouse and (B) the same image showing these particles. (C) Nuclear particles were dilated to generate rings that were used to sample Ca²⁺ intensity changes in the cytoplasm around the nucleus. (D) Tracking of the nuclei was achieved using a radial differential routine around the edge of the “nuclear” particle, which maximized the sum of the sampled intensity slopes calculated from cytoplasm to the nucleus. The tracked trajectories are overlaid. (E) The tracked nuclei were used to generate movement vector maps (red color of reference polygon represents a mostly longitudinal displacement) which, when reversed, could be applied to the original movie to stabilize distortion to the reference shape (F) (arrows show direction and magnitude of the distortion correction), or used to calculate longitudinal and circular micro-motions of the preparation. (G) Plots of Ca²⁺ intensity for each neuron were then used to identify Ca²⁺ transients based on the upstroke slope and duration of each Ca²⁺ transient to generate an event sequence (H) (white triangles). The distance and angle between neurons that fired sequentially could be represented as (I) a 3D event sequence path plot, or (J) the time and position of firing events around a particular reference cell could represented.

Distortion Vector Maps

The movement trajectories of tracked cell nuclei were saved and used to generate movement vector maps. Traditional motion-tracking routines calculate displacement and/or rotation of the whole field of view (FOV), but do not take into account differential rotation, contraction and elongation (distortion). To calculate distortion within an area (e.g., ganglia), a reference frame was chosen and a bounding polygon that encompassed the outer edges of multiple tracked points was generated (Figures 2E,F). Tessellation of the polygon into individual triangles using tracking points within the bounding polygon enabled better estimates of any compression/elongation within the bounding polygon. Bilinear interpolation was used to blend the relative movements between the 3 pairs of trajectories used to create each triangle, allowing either a color movement offset map (per pixel X [red] and Y [green] offsets, see Figure 2E) or vector map (per pixel angle and distance) to be constructed for every frame of the movie.

The movement vector maps were applied to the original movie allowing complete stabilization of displacement, rotation, compression and expansion (distortion) within the reference bounding polygon (Figure 2F). After distortion had been stabilized in an area, fixed ROIs or bitmasks could be applied to ganglia for calculation of spatio-temporal maps (ST maps) or other analyses without motion artifacts. As an added feature, the information stored in movement vector maps can be used to quantify differential circular and longitudinal displacement and contraction within the bounding polygon.

Non-Biased Measurement of Cytoplasmic Ca²⁺ Transients

After successful tracking of cells within ganglia using “nuclear” particles, the previously calculated “ring” masks were transferred to the tracked ROIs in the original movie, replacing the original “nuclear” particles masks. Ca²⁺-induced intensity changes were then measured in cells under the “ring” (Figures 2D,G). Average Ca²⁺-induced intensity changes within the “rings” are shown in (Figure 2G) although other intensity measurements are possible using this method, including the maximal intensity changes and the radial uniformity of Ca²⁺-induced fluorescence within the ring.

Ca²⁺ Transient Event Discrimination and Event Sequences

After traces of Ca²⁺-induced intensity changes were generated per cell, they were smoothed (±30 ms average) and a discriminator used to detect the upstroke velocity and amplitude of Ca²⁺ transients (3 s window: Figure 2H). The discriminated traces were smoothed (±0.5 s) and peaks identified above the half maximum amplitude level. The temporal position of each discriminated Ca²⁺ transient in each cell was stored as an event sequence that saved the ROI number, frame and classification structure of each cell (e.g., ganglia). This allowed the spatial and temporal characteristics of neuronal Ca²⁺ events to be analyzed as 3D representations to better illustrate the sequence of activation of cells during a CMMC. The event sequence path (Figure 2I) is a graphical representation of the order in which cells in a myenteric ganglion fired during a CMMC where time at which a cell fired during the recording is color-coded. Multiple stimulations using the same stimuli can be used to determine the stability in the order of Ca²⁺ events and associations between individual neurons during the sequence of activation. Relative event plots (Figure 2J) centers on a reference cell and generates space-time vectors to the position of other cells that fired before or after the reference cell.

Responses of Neurons and Glia to Different Stimuli

The lack of photobleaching of Ca²⁺-induced fluorescence in cells expressing GCaMP3 allowed experiments to be performed over many hours. This flexibility allowed a particular FOV to be imaged repeatedly during different stimulations or with drugs added. Advances in analysis allowed Ca²⁺ transients in individual cells within a ganglia could be extracted and compared. To illustrate the power of GCaMP3 in monitoring Ca²⁺ transients in neurons and glia we illustrate the power of the techniques by giving representative examples below in Figures 3–7. Anal elongation/mucosal stroking is usually a more effective stimulus to evoke a CMMC than oral elongation/mucosal stroking (Heredia et al., 2009, 2010, 2013; Bayguinov et al., 2010a). To ascertain the different responses in neurons and glia in a myenteric ganglia to oral and anal elongation, we used higher power imaging (20×) to track and record individual Ca²⁺ responses (Figure 3A). These evoked responses were compared to a spontaneous CMMC that occurred prior to the elongation stimuli. During a spontaneous CMMC (Figure 3B), a number of different patterns of Ca²⁺ transients were observed in individual cells: (i) complex Ca²⁺ transient activity (see Figure 3B: cells 1, 2 and 6) likely due to the variations in summation of ongoing Ca²⁺ transients; (ii) individual Ca²⁺ transients, likely associated with neuronal action potentials (see Figure 3B: cell 14); and (iii) long-duration Ca²⁺ transients (see Figure 3B: cells 4, 5 and 12) either representative of tetanic burst of action potentials in neurons or slower Ca²⁺ transients in glia. During the spontaneous CMMC, a cluster of activity in myenteric cells was observed ~10 s after the onset of contraction (see Figure 3B: bottom green and red traces representing longitudinal and circular displacement respectively) in eight cells, while three cells were quiescent throughout. During oral elongation, there was a ~5 s delay from the end of the stimulus to the start of contraction (see movement traces at bottom of Figure 3B). There was no apparent clustering of firing, although activity in three cells was increased (see Figure 3C: cells 2, 5 and 14). Two of the cells active during the spontaneous CMMC remained quiescent throughout the length of the recording after oral elongation (Figure 3C cells 7 and 12). In stark contrast to spontaneous or CMMCs evoked by oral elongation, anal elongation produced a considerable Ca²⁺ response in all measured cells clustered within 1–5 s of the cessation of the stimulus, and was coincident with a large circular and longitudinal contraction (Figure 3D: note the different scaling of traces. The Ca²⁺ transients next to the blue asterisks are equivalent in amplitude). This pronounced difference in response of myenteric neurons and glia to the polarity of elongation stimuli reveals the possible duality of the CMMC acting as both an ongoing “feed-forward” motor pattern that does not require local stimulation, as well as a response more similar to a peristaltic reflex with a strongly polarized, but long-lasting response to pronounced local stimulation.

FIGURE 3

Figure 3. Responses of myenteric ganglion cells to oral and anal elongation in Wnt1-GCaMP3 mice. (A) Myenteric ganglion cells that remained in the FOV (during the contraction responses) and had sufficient nuclear edge contrast to enable tracking (labeled from 1–15). A neuron resembling a Dogiel Type II neuron is indicated by the orange asterisk. (B) During a spontaneous CMMC, a number of neurons showed Ca²⁺ transients, particularly ~10 after the onset of contraction. (C) During oral elongation a similar overall pattern was observed, except a few cells remained quiescent throughout the recording period and the resulting contraction was larger (see green longitudinal displacement and red circular displacement traces at bottom). (D) In marked contrast, anal elongation produced a robust Ca²⁺ response in all measured neurons, all within 1–2 s of each other following the elongation stimulus. Note the difference in the intensity scale bar between (B,C) and (D). The Ca²⁺ transients indicated by the blue asterisks in (C,D) are of comparable amplitude.

FIGURE 4

Figure 4. Overall pattern of activation of all enteric neurons and glia in myenteric and submucous plexi using Wnt1-GCaMP3-labeled neural crest derived cells. (A) Pseudo-colored sequence of images depicting the Ca²⁺-induced fluorescence of myenteric neurons and glia together during a CMMC evoked by anal elongation. Note the staggered activation of patches of neurons/glia in multiple ganglia with the rightmost ganglia displaying the most number of active neurons/glia initially followed, somewhat haphazardly by the middle and leftmost ganglia later on during the CMMC. (B) Temporal delay map was constructed by overlaying the area of thresholded Ca²⁺ transients at different time points and illustrates this apparent spread graphically over 15 s. The structure near the blue asterisk is a submucosal ganglia. (C) ST Map of the apparent propagation of activity between ganglia allowing an overall velocity of spread to be estimated. (D) A trace showing the time course of the distortion associated with the CMMC. (E) Higher-power image of an underlying smaller submucosal ganglia and (F) a plot of Ca²⁺-induced fluorescence in a selected cell in this ganglia.

FIGURE 5

Figure 5. Activation sequence of neurons and glia in colonic myenteric ganglia using Wnt1-GCaMP3 and ChAT- GCaMP3 labeled cells. (A) Low- power image of multiple myenteric ganglia in which neurons and glia were both labeled with Wnt1-GCaMP3. (B) Overall measurement of Ca²⁺ activity from both cell types from the entire ganglia (blue trace) shows a build-up of Ca²⁺ over the first ~4 s of the CMMC followed by a biphasic decline in Ca²⁺ intensity over the following ~ 30 s. The CMMC was evoked by anal mucosal stimulation. Measuring selected regions in the ganglia that contained primarily neurons (green trace; see green dashed circle in C: 4 s) or glia (magenta trace: see magenta dashed circle in C: 26 s) based on their morphology illustrates how Ca²⁺ transients in neurons preceded glial activity by 1–2 s. This is better illustrated in the sequential images displayed in (C) where the mesh-like appearance of glia is visible in the 7 s panel, but not at 4 s panel. (D) To determine whether the number of neurons that fired affected the activation of glia, we added a submaximal dose of hexamethonium (5 μM) to compromise cholinergic nicotinic transmission between enteric neurons. Stimulation of the mucosa some distance away resulted in very few neurons becoming active at the site of recording and no observable glial response. Images were background subtracted to highlight those cells that fired after the stimulus. (E) To examine the neuronal response during the complex without concomitant glial activity, the activity in ChAT-GCaMP3 cholinergic neurons was analyzed. Images at 3 time points are displayed in pseudo-color and the time course of Ca²⁺ intensity in a small group of neurons (F) shows the characteristics of the increase of activity over the first 3–7 s, followed by a more rapid decline in the subsequent 7 s.

FIGURE 6

Figure 6. Ca²⁺ activity in enteric neurons expressing nNOS- GCaMP3 and ChAT- GCaMP3. (A_i) Image of nNOS-GCaMP3 neurons in colonic myenteric ganglia and (A_ii) an activity map of those neurons displaying spontaneous Ca²⁺ transients. A trace of spontaneous Ca²⁺ activity is shown at the bottom of the panel. (B) Image in the left panel shows another myenteric ganglia in which (right panel) the Ca²⁺ activity in selected neurons was plotted in relation to the distortion produced by a spontaneous CMMC (bottom trace). nNOS neurons displayed different responses during the onset of the CMMC: the top two traces show the activation of two formerly silent neurons, while the bottom two traces show the cessation/silencing of two formerly active neurons. These likely represent interneurons and inhibitory motor neurons respectively. (C_i) Low-power image of ChAT-GCaMP3-expressing neurons in multiple colonic myenteric ganglia from which the rightmost ganglia was magnified (C_ii) and an activity map created (C_iii). (D) Image in the left panel shows another ganglia in which the Ca²⁺ activity (right panel) in selected ChAT-GCaMP3 neurons was plotted. These neurons responded with excitatory responses during the initial phase of a CMMC evoked by mucosal stimulation (three strokes; see bottom trace) and likely represent a mixed population of IPANs, interneurons and excitatory motor neurons.

FIGURE 7

Figure 7. Ca²⁺ activity in enteric glia expressing GFAP-GCaMP3 during the colonic CMMC. (A) Fluorescence image of glial cell bodies and their extensive processes in a colonic myenteric ganglia. In addition to these glia, two neurons appear to express GCaMP3, as shown in the bottom left hand corner. (B) Colored temporal delay map illustrates that the activation of enteric glia and their processes during a CMMC was not uniform, but occurred over a period of 4.5 s in a rather “patchy” pattern. (C) Sequence of thresholded images used to generate the map in (B) revealed significant delays between the activation of different parts of the glial network in the same ganglia. (D) In some regions, the pattern of activation allowed the rate of spread to be calculated using ST maps (see dashed region in A). (E) The overall increase in GFAP-GCaMP3 Ca²⁺ intensity within a myenteric ganglia appears tightly correlated to the initiation of distortion (movement).

Propagation Patterns of Ca²⁺ Activity in Enteric Neurons and Glia During the CMMC in Wnt1-GCaMP3 Mice

To demonstrate how low-power (10×) imaging can resolve propagation of neural activity associated with a CMMC, the example in Figure 4 shows a FOV ~1 mm² in size that encompasses three rows of ganglia. Lower-power imaging (10×) allowed single myenteric ganglion cells to be resolved and reduced the pixel shift associated with movement. Mucosal stimulation with a brush or a brief anal elongation were both effective stimuli to evoke a CMMC (Heredia et al., 2009, 2010, 2013). In Figure 4, anal elongation evoked a complex pattern of neuron and glial Ca²⁺ transients. Immediately after the stimulation, occasional neuron cell bodies scattered throughout the FOV became active (Figure 4A). The cells consisted of a variety of neurons, some of which were often large neurons that were ovoid in shape, similar to Dogiel Type 2 neurons (Figure 4A, 0 s; Bayguinov et al., 2010a; Okamoto et al., 2014). Over the following 10 s, there appeared to be a slow propagating wave of activity in which large numbers of myenteric ganglion cells displayed robust, prolonged Ca²⁺ transients going from right (anal) to left (oral) (see Figure 4A panels 0–10 s). This propagation can be better appreciated in the temporal delay map shown in Figure 4B where the time of peak Ca²⁺-induced fluorescence intensity was color-coded (spectrum). Significant delays in the activation of myenteric ganglion cells can be appreciated between, but also within individual ganglia (see middle ganglion in Figure 4B). The 10× lens has a large focal depth that allows for the significant delay in activation of cells in a submucosal ganglion located on a different focal plane on the left hand side of the FOV to be examined (Figure 4B: see blue asterisk next to ganglia). The speed of apparent propagation of the wave of activation of myenteric ganglion cells was calculated to be ~130 μms⁻¹ from the ST map shown in Figure 4C. Distortion produced by muscular contraction began ~5 s after the initial activation of myenteric ganglion cells, with the most rapid phase of contraction occurring when the largest number of myenteric ganglion cells were active in the FOV. The distortion (contraction) continued to increase for some time (~5 s) after Ca²⁺ transients in these cells had dissipated (Figures 4C,D). These results demonstrate that an evoked CMMC involves a slow spreading but strong activation of most cells in myenteric ganglia, followed significantly later by the activation of cells in submucosal ganglia. Submucosal ganglion cells could be imaged in detail using a higher-power (20×) lens and displayed ongoing Ca²⁺ transients between complexes (Figures 4E,F).

Temporal Activation Pattern of Neurons and Glia During a CMMC

When we examined the time courses of Ca²⁺ activity from myenteric neurons and glia (Figures 5A,B), it was clear that Ca²⁺ transients in glia cells occurred after the initial responding neurons (compare Figure 5C at 4 s with that at 7 s). Neurons were distinguished based on their larger size (nucleus and cytoplasmic area). The Ca²⁺ intensity of cells in a region of the ganglia that did not contain active neurons showed a 2–3 s delay compared to a region composed of mainly neurons (magenta and green circles in Figure 5C respectively). The combined Ca²⁺ intensity from the whole ganglion (Figure 5B blue line), neurons (Figure 2B green line) and glia (Figure 5B magenta line) all showed a biphasic response with a fast increase in overall Ca²⁺-induced fluorescence peaking at ~5 s after the first responders, followed by a fast decrease of Ca²⁺-induced fluorescence lasting another ~5 s. Following the faster rise and fall in overall Ca²⁺-induced fluorescence, there was a prolonged “plateau” or slowly declining phase lasting 20–30 s. During this phase, patches of glia cells were observed to fire sporadically as illustrated by the small peak in the glial trace (Figure 5B, 26 s).

To begin to examine the role of neuron-to-neuron communication on the responses of glia during a CMMC, we used a low dose of hexamethonium (5–10 μM, n = 4) to compromise cholinergic nicotinic transmission between neurons. This abolished nearly all neural activity, although a few neurons per ganglia were observed to display Ca²⁺ transients in response to mucosa or stretch stimulation (Figure 5D, green circle). However, no Ca²⁺ transients in glia were observed, suggesting either a critical mass of active neurons is required to activate glia, and/or that nicotinic transmission affects glia directly (Broadhead et al., 2012; Boesmans et al., 2013a; Petrov et al., 2014). To examine the response of neurons alone in more detail, we measured overall intensity of Ca²⁺ in a ganglion from a ChAT-GCaMP3 mouse. Similar to responses obtained from regions in Wnt1-GCaMP3 mice containing all classes of neurons, the time course of Ca²⁺-induced intensity changes from a group of ChAT-GCaMP3 neurons showed a similar biphasic response to anal elongation without glial contamination (Figures 5E,F).

Ca²⁺ Activity in nNOS-GCaMP3 Neurons During the CMMC

Between CMMCs, a number nNOS-GCaMP3 neurons exhibited spontaneous Ca²⁺ transients (Figure 6: see arrowhead and trace in Figures 6Ai,ii,B and Figure 6B: neurons 3 and 4). During a spontaneous CMMC, a proportion of nNOS-GCaMP3 neurons reduced their Ca²⁺ transients (Figure 6B: neurons 3 and 4), consistent with these neurons being inhibitory motor neurons (Bayguinov et al., 2010a). Other nNOS-GCaMP3 neurons, which exhibited little spontaneous activity between CMMCs, displayed long, sustained trains of Ca²⁺ transients at the onset of the CMMC (Figure 6B: neurons 1 and 2). These neurons are likely to be a subclass of interneurons as has been reported previously using Ca²⁺ imaging via conventional Ca²⁺ dyes (Bayguinov et al., 2010a).

Ca²⁺ Activity in ChAT-GCaMP3 Neurons During the CMMC

ChAT-GCaMP3 neurons were more numerous than nNOS-GCaMP3 neurons (compare Figures 6Ai,ii,B with Figures 6Ci,ii,D), which is not surprising since cholinergic neurons (~60% of all myenteric neurons) comprise many functionally different classes of neurons than nNOS neurons (Lomax and Furness, 2000; Hao et al., 2013). Some ChAT-GCaMP3 neurons exhibited ongoing Ca²⁺ transients between CMMCs and are likely to be mainly interneurons and inhibitory motor neurons (Figure 6D: neurons 1 and 2) as tonic inhibition dominates during this period (Dickson et al., 2010a,b). However, many ChAT-GCaMP3 neurons increased their activity during a CMMC (Figure 6D: neuron 4) and are likely to be excitatory motor neurons (Hao et al., 2013). Furthermore, by focusing deeper into the tissue through the circular muscle apparently separate classes of active nNOS and ChAT neurons could be observed in the submucous plexus and Henle’s plexus without removing the mucosa as was previously necessary (data not shown: see Okamoto et al., 2012, 2014).

These preliminary results mirror previous studies describing the behavior of NOS and ChAT neurons using Ca²⁺ dyes (Bayguinov et al., 2010a), but with a number of advantages including: (i) neurons could be visualized in preparations that were undissected (fully intact tubes with mucosa); (ii) no post hoc staining was needed to determine the class of neuron, as only those neurons with nNOS or ChAT expressed the Ca²⁺ indicator; (iii) the signal to noise of Ca²⁺-induced fluorescence in these cells enabled activity in these neurons to be observed in multiple ganglia simultaneously; and (iv) the lack of photobleaching and continual replenishment of GCaMP3 due to ongoing expression allowed activity in enteric neurons to be monitored for hours, allowing the effect of different drugs or stimulations to be observed in the same FOV. Morphological and functional differences in ChAT and nNOS enteric neurons allowed some differentiation of whether particular neurons were sensory (IPAN), ascending or descending inter- or motor neurons, but for unequivocal classification, live or post hoc labeling for known markers of these different functional classes of neurons will be necessary (Bayguinov et al., 2010a, 2012).

Ca²⁺ Activity in GFAP-GCaMP3 Labeled Glia During the CMMC

We also studied the characteristic of activation of the myenteric glia network in GFAP-CaMP3 mice during the CMMC (Figure 7). On occasion, 1 or 2 large neurons, including Type II neurons, were labeled in GFAP-GCaMP3 mice (Figure 7A, see bottom two neurons). It has also been shown that a subpopulation of GFAP-positive progenitor cells can develop into neurons in the central nervous system (CNS; Casper and McCarthy, 2006) and in the ENS it has been shown that glial cells can gain a neurogenic potential in certain conditions (Joseph et al., 2011; Laranjeira et al., 2011). In any case, the labeling of these neurons was easy to exclude from analysis of glial response based on their size and background intensity. In GFAP-GCaMP3 mice, myenteric glial cell somas had a size of 14.0 ± 2.2 μm (long axis) and 6.9 ± 0.9 μm (short axis; 60 cells, n = 4) that is consistent with them being much smaller than neurons (Broadhead et al., 2012).

There was little activation of glial cells between CMMCs, although on occasions some random glial cell processes exhibited Ca²⁺ transients. However, the pattern of Ca²⁺ transients in the myenteric glial network following initiation of a CMMC was unique. Following anal elongation there was a slow spread of activity through the myenteric glia (Figures 7A–C) that can best be described as “stepped” or “patchy”. While we do not know the corresponding neural activity occurring within the myenteric ganglia, the propagation of activity was not smooth or coherent, nor could a wave front be resolved (see temporal delay map in Figure 7B). Whether there is bona fide propagation in the glial network (Figure 7D) or whether the activation of glia is dependent on adjacent neural activity remains to be determined. However, taking overall Ca²⁺ transient intensity throughout the ganglia shows a steady build-up of activity over 3–4 s during the initial phase of the CMMC (Figure 7E).

Discussion

Use of Ca²⁺ Indicators to Monitor Enteric Neurons and Glia

Calcium indicators such as Fluo3/4 have been used effectively to measure Ca²⁺ transients in enteric neurons, glia, pacemaker cells, smooth muscle and varicosities in order to monitor their activity during gut motor behaviors (Stevens et al., 1999; Spencer et al., 2006; Bayguinov et al., 2010a,b, 2012; Broadhead et al., 2012; Okamoto et al., 2012; Boesmans et al., 2013a,b; Smith et al., 2015). However, these indicators can be subject to photobleaching, dispersal from certain target cells, de-esterification issues and often need complex and time consuming dissection to target the dye to the cells of interest (Bayguinov et al., 2010a,b; Okamoto et al., 2012). While the kinetics of GCaMP3 are slower than Ca²⁺ indicator dyes, we chose this Ca²⁺ indicator as it has the brightest fluorescence (quantum yield) suitable for detecting Ca²⁺ transients when using low power imaging. The firing frequency of enteric neurons is usually less than 10 Hz (Mazzuoli and Schemann, 2012), and can likely resolve most Ca²⁺ transients generated by action potentials (Tian et al., 2009; Yamada and Mikoshiba, 2012).

In this study in which the GCaMP3 Ca²⁺ indicator was directly expressed in cells of interest, Ca²⁺-induced GCaMP3 fluorescence was bright enough that no dissection was required and visualization of either myenteric or submucosal neurons/glial cells could be readily made along the whole gut by focusing through the longitudinal or circular muscle layers respectively in an intact colon. Therefore, this method leaves the structure and connectivity of these plexuses and their relation to other cells intact, including the mucosa. Moreover, the shapes and relative sizes of different functional classes of neurons are well preserved (Lomax and Furness, 2000). We demonstrated that GCaMP3 appears to be good indicator of Ca²⁺ transient activity in enteric neurons and glia in the gut at physiological temperatures in Wnt1-GCaMP3 (all ENS neurons and glia), ChAT-GCaMP3 (enteric cholinergic neurons), nNOS-GCaMP3 (enteric nNOS neurons) and GFAP-GCaMP3 (enteric glia) mice both between and during the powerful, neurally-mediated CMMC (Bywater et al., 1989; Heredia et al., 2009; Bayguinov et al., 2010a).

Although there is some variation in the brightness of Ca²⁺ transients in different neurons and glia, this likely relates to variation in their Ca²⁺ handling ability. Importantly, myenteric neurons expressing GCaMP3 appear to respond normally to physiological stimuli such as mucosal stroking or elongation suggesting that the neural reflex circuitry is intact and not altered by the expression of the foreign protein or by the genetic manipulations necessary to generate these mice.

Also, by using ChAT-GCaMP3 and nNOS-GCaMP3 mice we can begin monitoring Ca²⁺ transient activity in broadly defined cholinergic and nitrergic neurons in myenteric and submucosal ganglia, without the necessity of time-consuming post-staining immunohistochemistry, as was previously necessary (see Bayguinov et al., 2010a,b, 2012; Broadhead et al., 2012; Okamoto et al., 2012). One caveat in utilizing such strategies, however, is that the expression of transgenes using promoters for cell-specific markers is not always matched to the endogenous expression pattern of the marker. For example, in a pair of studies using transgenic ChAT-Cre mice, it was found that transgene expression was not observed in all cells positive for ChAT-immunoreactivity (IR) and conversely that transgene expression was observed in a small percentage of cells negative for ChAT-IR (Hao et al., 2013; Erickson et al., 2014). The observation that expression of GFP in transgenic ChAT-GFP mice, as well as tdTomato expression in ChAT-Cre; Rosa26-Lox-STOP-Lox-tdTomato mice, are both found in ChAT-IR-negative cells (Erickson et al., 2014) supports the idea that this latter pattern represents ectopic transgene expression by ChAT-Cre mice rather than transient expression of endogenous ChAT, since transgenic ChAT-GFP mice do not lineage-trace cells the way ChAT-Cre; Rosa26-Lox-STOP-Lox-tdTomato mice do. In the future, these techniques will be able to be extended to identify the roles of subclasses of enteric neurons with different chemical markers.

Characteristics of the CMMC Revealed by GCaMP3 Expressed in Different Classes of Enteric Neurons and Glia

The CMMC is a propagating contraction that is essentially a rhythmic peristaltic wave involving most cells across the colonic wall including myenteric and submucosal pacemaker cells and mucosal 5-HT release from enterochromaffin cells (Bayguinov et al., 2010a,b; Smith and Gershon, 2015; Smith et al., 2015). The CMMC appears to be equivalent to the high amplitude propagating contractions (HAPCs) observed in human colon (Smith et al., 2015).

Several ideas have been previously proposed for the generation and propagation of neurally-mediated motor complexes along the small and large intestine (see Smith et al., 2015). These ideas usually emphasize particular elements of the ENS as being more or less important for propulsive motor behaviors, e.g., a turning off of tonic inhibition, recurrent connections between IPANs, peristaltic reflexes and local reflexes generated by a pellet (Christensen et al., 1974; Bywater et al., 1989, 1998; Furness et al., 1996; Thomas et al., 2004; Heredia et al., 2009, 2013; Dinning et al., 2014). However, motor complexes likely invariably involve two main events: (1) generating an imbalance between inhibition and excitation to the smooth muscle to favor contraction and (2) propagating this activity along the gut. Because of gut movement, few groups have observed the interaction of all neurons in many ganglia covering a respectably-sized region of gut to verify their individual hypotheses (Bayguinov et al., 2010a,b; Okamoto et al., 2012).

What we deduce from our current studies using GCaMP3 mice is that CMMC generation likely involves the activation and inhibition of many types of neurons such as IPANs, interneurons, cholinergic motor neurons, inhibitory motor neurons, as well as glial cells. Immediately following a stimulus to activate a CMMC, a number of myenteric neurons, including likely Dogiel Type II neurons (IPANs), were activated; their activity lasts through much of the CMMC but is most pronounced in the first 10 s. However, subsequent to the initial activation of neurons and glia, there was a long lasting “plateau” or slowly decaying phase of activity over 20–30 s that is similar (albeit somewhat shorter) than the overall contractile responses occurring during a CMMC contraction (Dickson et al., 2010a). Furthermore, the apparent slow propagation of activation of neurons in adjacent ganglia in the myenteric plexus (~130 μms⁻¹) may contribute to the slow propagation of the contraction associated with the CMMC (0.8s⁻¹). The CMMC does not appear to be exclusively a simple peristaltic reflex, as described for the small intestine (Furness et al., 1996) that conducts rapidly through the myenteric plexus (Heredia et al., 2013), but is much slower, more elaborate and long-lasting—likely involving both a simple reflex activated immediately after the stimulus, followed by longer lasting complex interactions between neurons and glia.

The most effective stimuli for generating a CMMC was anal mucosal stimulation or anal elongation of the colon (Heredia et al., 2009, 2010) since these stimuli elicited a synchronized, prolonged and robust burst of Ca²⁺ transients compared to the activation of the same myenteric neurons observed during a spontaneous CMMC. Presumably this is because descending inhibitory nerve pathways are cut at the anal end of the colon, whereas ascending neural pathways are intact, although the point of origin of a spontaneous CMMC is difficult to detect (Spencer et al., 2005; Heredia et al., 2009, 2013; Bayguinov et al., 2010a).

The complexity of the neural circuitry underlying motor complexes is suggested by the fact that while a small percentage of neurons could be activated by stimuli in the presence of submaximal doses of hexamethonium (nicotinic antagonist), there was no CMMC-like contraction or glial response, indicating a critical number of neural connections (involving nicotinic synapses) are needed to generate the CMMC. In future, other characteristics of the CMMC can be explored in detail using GECIs, such as the observations that 5-HT antagonists (acting at 5-HT1A, 5-HT3, 5-HT4 and 5-HT7 receptors) inhibit or block the CMMC contraction (Heredia et al., 2009, 2013; Dickson et al., 2010b; Smith and Gershon, 2015; Smith et al., 2015; Yu et al., 2015). Although the CMMC contraction increases in frequency following the inhibition of nitric oxide synthesis, it doesn’t propagate effectively down the colon (Heredia et al., 2010). The relative importance and contribution of cellular activity of cholinergic, serotonergic and nitric oxide neurons can be assessed more effectively in mice in which GCaMP3 is expressed in specific cells.

New Analyses to Extract Information from Intact Segments of Gut

A formidable challenge in monitoring gut neural activity is the ongoing distortion of the tissue generated by spontaneous movements of the colon. Movement compensation is especially critical to monitor in the gut since the ganglionic distortion may be a stimulus to activate certain classes of neurons (Spencer et al., 2003; Mazzuoli and Schemann, 2012). Others have attempted to circumvent the gut movement problem by a number of ways including: (i) recording Ca²⁺ activity from neurons or glial cells either in small preparations of myenteric ganglia from which the longitudinal or circular smooth muscle has been removed; (ii) using cultured neurons or explants of myenteric ganglia; or (iii) using lower non-physiological temperatures or using L-type Ca²⁺ channel blockers or drugs that interfere with actin/myosin interactions (e.g., wortmannin/ML7) to reduce contractile activity (Okamoto et al., 2012; Boesmans et al., 2013a,b, 2015; McClain et al., 2014). Each of these methods invariably introduces some degree of artifact into the responses of cells or tissue.

With the ability to image neurons and glia in undissected tissue with very little background staining, existing analysis methods were inadequate. Especially at low power, the distortions produced by regional contractions and relaxations were never uniform across the FOV. Similarly, the lack of granular background staining (common with Ca²⁺ dyes) meant that the only objects that were distinct enough to track, also had the most changes in brightness due to Ca²⁺ transients. Nuclear differential tracking techniques combined with vector maps appear to offer an elegant solution to extracting data from tissue labeled with GECIs. The event sequence plots allow for detailed examination of activation patterns in relation to the whole FOV, or individual neurons, and can be used to compare the effect of drugs or stimulation patterns in minute detail.

In conclusion, the use of GECIs to study neuronal and glial behavior allows gut preparations to be left fully intact, and allows examination of the complex interactions between both myenteric and submucosal ganglia. As GECIs continue to develop (Akerboom et al., 2012), their excellent signal to noise ratio in a near physiological setting will allow GECIs in in vivo experiments to be performed with equal ease. By incorporating the new forms of analysis, the full complexity of neuronal and glial Ca²⁺ events that generate the CMMC or other gut behaviors involving neurons and glia can begin to be elucidated both in vitro and in vivo.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This project was supported by grants from the National Institutes of Health: 5P20RR018751-09 (GWH and TKS); 8P20GM103554 (TWG); PO1 NIH/NIDDK DK41315-26 (SDK). TWG produced the GFAP-GCaMP3, Wnt1-GCaMP3, ChAT-GCaMP3 mice. DJH assisted TWG. SDK produced the nNOS-GCaMP3 mouse. TKS, SDK, RJC and MCS performed the imaging experiments and analyzed and selected recordings. GWH developed the analysis software, analyzed the data and constructed the figures. GWH, TWG and TKS wrote the initial draft and finalized the manuscript that was reviewed and approved by all investigators.

Abbreviations

ChAT, Choline acetyl transferase; CMMC, colonic migrating motor complex; ENS, Enteric nervous system; FOV, field of view; GCaMP3, Genetically encoded Ca²⁺ indicator is fused to the calmodulin (CaM) and the M13 domain of a myosin light chain kinase; GECI, Genetically-encoded Ca²⁺ indicator; GFAP, Glial fibrillary acidic protein; nNOS, neuronal nitric oxide synthase; ROI, region of interest; Wnt1, wingless-type MMTV integration site family.

References

Akerboom, J., Chen, T. W., Wardill, T. J., Tian, L., Marvin, J. S., Mutlu, S., et al. (2012). Optimization of a GCaMP calcium indicator for neural activity imaging. J. Neurosci. 32, 13819–13840. doi: 10.1523/JNEUROSCI.2601-12.2012

PubMed Abstract | CrossRef Full Text | Google Scholar

Bayguinov, P. O., Broadhead, M. J., Okamoto, T., Hennig, G. W., and Smith, T. K. (2012). Activity in varicosities within the myenteric plexus between and during the colonic migrating motor complex in the isolated murine large intestine. Neurogastroenterol. Motil. 24, e185–e201. doi: 10.1111/j.1365-2982.2012.01892.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Bayguinov, P. O., Hennig, G. W., and Smith, T. K. (2010a). Calcium activity in different classes of myenteric neurons underlying the migrating motor complex in the murine colon. J. Physiol. 588, 399–421. doi: 10.1113/jphysiol.2009.181172

PubMed Abstract | CrossRef Full Text | Google Scholar

Bayguinov, P. O., Hennig, G. W., and Smith, T. K. (2010b). Ca2+ imaging of activity in ICC-MY during local mucosal reflexes and the colonic migrating motor complex in the murine large intestine. J. Physiol. 588, 4453–4474. doi: 10.1113/jphysiol.2010.196824

PubMed Abstract | CrossRef Full Text | Google Scholar

Bayguinov, P. O., Hennig, G. W., and Smith, T. K. (2010c). Generation of complex neural behavior in a mammalian nervous system. Physiol. News 80, 29–32.

Boesmans, W., Cirillo, C., Van den Abbeel, V., Van den Haute, C., Depoortere, I., Tack, J., et al. (2013a). Neurotransmitters involved in fast excitatory neurotransmission directly activate enteric glial cells. Neurogastroenterol. Motil. 25, 151–160. doi: 10.1111/nmo.12065

PubMed Abstract | CrossRef Full Text | Google Scholar

Boesmans, W., Martens, M. A., Weltens, N., Hao, M. M., Tack, J., Cirillo, C., et al. (2013b). Imaging neuron-glia interactions in the enteric nervous system. Front. Cell. Neurosci. 7:183. doi: 10.3389/fncel.2013.00183

PubMed Abstract | CrossRef Full Text | Google Scholar

Boesmans, W., Hao, M. M., and Vanden Berghe, P. (2015). Optical tools to investigate cellular activity in the intestinal wall. J. Neurogastroenterol. Motil. 21, 337–351. doi: 10.5056/jnm15096

PubMed Abstract | CrossRef Full Text | Google Scholar

Broadhead, M. J., Bayguinov, P. O., Okamoto, T., Heredia, D. J., and Smith, T. K. (2012). Ca2+ transients in myenteric glial cells during the colonic migrating motor complex in the isolated murine large intestine. J. Physiol. 590, 335–350. doi: 10.1113/jphysiol.2011.219519

PubMed Abstract | CrossRef Full Text | Google Scholar

Bywater, R. A., Small, R. C., and Taylor, G. S. (1989). Neurogenic slow depolarizations and rapid oscillations in the membrane potential of circular muscle of mouse colon. J. Physiol. 413, 505–519. doi: 10.1113/jphysiol.1989.sp017666

PubMed Abstract | CrossRef Full Text | Google Scholar

Bywater, R. A., Spencer, N. J., Fida, R., and Taylor, G. S. (1998). Second-, minute- and hour-metronomes of intestinal pacemakers. Clin. Exp. Pharmacol. Physiol. 25, 857–861. doi: 10.1111/j.1440-1681.1998.tb02167.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Casper, K. B., and McCarthy, K. D. (2006). GFAP-positive progenitor cells produce neurons and oligodendrocytes throughout the CNS. Mol. Cell. Neurosci. 31, 676–684. doi: 10.1016/j.mcn.2005.12.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Christensen, J., Anuras, S., and Hauser, R. L. (1974). Migrating spike bursts and electrical slow waves in the cat colon: effect of sectioning. Gastroenterology 66, 240–247.

PubMed Abstract

Dickson, E. J., Heredia, D. J., McCann, C. J., Hennig, G. W., and Smith, T. K. (2010a). The mechanisms underlying the generation of the colonic migrating motor complex in both wild-type and nNOS knockout mice. Am. J. Physiol. Gastrointest. Liver Physiol. 298, G222–G232. doi: 10.1152/ajpgi.00399.2009

PubMed Abstract | CrossRef Full Text | Google Scholar

Dickson, E. J., Heredia, D. J., and Smith, T. K. (2010b). Critical role of 5-HT1A, 5-HT3 and 5-HT7 receptor subtypes in the initiation, generation and propagation of the murine colonic migrating motor complex. Am. J. Physiol. Gastrointest. Liver Physiol. 299, G144–G157. doi: 10.1152/ajpgi.00496.2009

PubMed Abstract | CrossRef Full Text | Google Scholar

Dinning, P. G., Wiklendt, L., Omari, T., Arkwright, J. W., Spencer, N. J., Brookes, S. J., et al. (2014). Neural mechanisms of peristalsis in the isolated rabbit distal colon: a neuromechanical loop hypothesis. Front. Neurosci. 8:75. doi: 10.3389/fnins.2014.00075

PubMed Abstract | CrossRef Full Text | Google Scholar

Erickson, C. S., Lee, S. J., Barlow-Anacker, A. J., Druckenbrod, N. R., Epstein, M. L., and Gosain, A. (2014). Appearance of cholinergic myenteric neurons during enteric nervous system development: comparison of different ChAT fluorescent mouse reporter lines. Neurogastroenterol. Motil. 26, 874–884. doi: 10.1111/nmo.12343

PubMed Abstract | CrossRef Full Text | Google Scholar

Foong, J. P., Hirst, C. S., Hao, M. M., McKeown, S. J., Boesmans, W., Young, H. M., et al. (2015). Changes in nicotinic neurotransmission during enteric nervous system development. J. Neurosci. 35, 7106–7115. doi: 10.1523/JNEUROSCI.4175-14.2015

PubMed Abstract | CrossRef Full Text | Google Scholar

Furness, J. B. (2012). The enteric nervous system and neurogastroenterology. Nat. Rev. Gastroenterol. Hepatol. 9, 286–294. doi: 10.1038/nrgastro.2012.32

PubMed Abstract | CrossRef Full Text | Google Scholar

Furness, J. B., Bornstein, J. C., Kunze, W. A., Bertrand, P. P., Kelly, H., and Thomas, E. A. (1996). Experimental basis for realistic large-scale computer simulation of the enteric nervous system. Clin. Exp. Pharmacol. Physiol. 23, 786–792. doi: 10.1111/j.1440-1681.1996.tb01180.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Gulbransen, B., and Sharkey, K. A. (2009). Purinergic neuron-glia signaling in the enteric nervous system. Gastroenterology 136, 1349–1358. doi: 10.1053/j.gastro.2008.12.058

PubMed Abstract | CrossRef Full Text | Google Scholar

Gulbransen, B., and Sharkey, K. A. (2012). Novel functional roles for enteric glia in the gastrointestinal tract. Nat. Rev. Gastroenterol. Hepatol. 9, 625–632. doi: 10.1038/nrgastro.2012.138

PubMed Abstract | CrossRef Full Text | Google Scholar

Hao, M. M., Bornstein, J. C., and Young, H. M. (2013). Development of myenteric cholinergic neurons in ChAT-Cre;R26R-YFP mice. J. Comp. Neurol. 521, 3358–3370. doi: 10.1002/cne.23354

PubMed Abstract | CrossRef Full Text | Google Scholar

Heredia, D. J., Dickson, E. J., Bayguinov, P. O., Hennig, G. W., and Smith, T. K. (2009). Localized release of serotonin (5-hydroxytryptamine) by a fecal pellet regulates migrating motor complexes in murine colon. Gastroenterology 136, 1328–1338. doi: 10.1053/j.gastro.2008.12.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Heredia, D. J., Dickson, E. J., Bayguinov, P. O., Hennig, G. W., and Smith, T. K. (2010). Colonic elongation inhibits pellet propulsion and migrating motor complexes in the murine large bowel. J. Physiol. 588, 2919–2934. doi: 10.1113/jphysiol.2010.191445

PubMed Abstract | CrossRef Full Text | Google Scholar

Heredia, D. J., Gershon, M. D., Koh, S. D., Corrigan, R. D., Okamoto, T., and Smith, T. K. (2013). Important role of mucosal serotonin in colonic propulsion and peristaltic reflexes: in vitro analyses in mice lacking tryptophan hydroxylase 1. J. Physiol. 591, 5939–5957. doi: 10.1113/jphysiol.2013.256230

PubMed Abstract | CrossRef Full Text | Google Scholar

Heredia, D. J., Grainger, N., McCann, C. J., and Smith, T. K. (2012). Insights from a novel model of slow-transit constipation generated by partial outlet obstruction in the murine large intestine. Am. J. Physiol. Gastrointest. Liver Physiol. 303, G1004–G1016. doi: 10.1152/ajpgi.00238.2012

PubMed Abstract | CrossRef Full Text | Google Scholar

Joseph, N. M., He, S., Quintana, E., Kim, Y. G., Núñez, G., and Morrison, S. J. (2011). Enteric glia are multipotent in culture but primarily form glia in the adult rodent gut. J. Clin. Invest. 121, 3398–3411. doi: 10.1172/JCI58186

PubMed Abstract | CrossRef Full Text | Google Scholar

Laranjeira, C., Sandgren, K., Kessaris, N., Richardson, W., Potocnik, A., Vanden Berghe, P., et al. (2011). Glial cells in the mouse enteric nervous system can undergo neurogenesis in response to injury. J. Clin. Invest. 121, 3412–3424. doi: 10.1172/JCI58200

PubMed Abstract | CrossRef Full Text | Google Scholar

Lomax, A. E., and Furness, J. B. (2000). Neurochemical classification of enteric neurons in the guinea-pig distal colon. Cell Tissue Res. 302, 59–72. doi: 10.1007/s004410000260

PubMed Abstract | CrossRef Full Text | Google Scholar

MacEachern, S. J., Patel, B. A., McKay, D. M., and Sharkey, K. A. (2011). Nitric oxide regulation of colonic epithelial ion transport: a novel role for enteric glia in the myenteric plexus. J. Physiol. 589, 3333–3348. doi: 10.1113/jphysiol.2011.207902

PubMed Abstract | CrossRef Full Text | Google Scholar

Mazzuoli, G., and Schemann, M. (2012). Mechanosensitive enteric neurons in the myenteric plexus of the mouse intestine. PLoS One 7:e39887. doi: 10.1371/journal.pone.0039887

PubMed Abstract | CrossRef Full Text | Google Scholar

McClain, J. L., Fried, D., Gomez-Suarez, R. A., Leinninger, G. M., Sévigny, J., Parpura, V., et al. (2014). Ca2+ responses in enteric glia are mediated by connexin-43 hemichannels and modulate colonic transit in mice. Gastroenterology 146, 497–507. doi: 10.1053/j.gastro.2013.10.061

PubMed Abstract | CrossRef Full Text | Google Scholar

Okamoto, T., Barton, M. J., Hennig, G. W., Birch, G. C., Grainger, N., Corrigan, R. D., et al. (2014). Extensive projections of myenteric serotonergic neurons suggest they comprise the central processing unit in the colon. Neurogastroenterol. Motil. 26, 556–570. doi: 10.1111/nmo.12302

PubMed Abstract | CrossRef Full Text | Google Scholar

Okamoto, T., Bayguinov, P. O., Broadhead, M. J., and Smith, T. K. (2012). Ca(2+) transients in submucous neurons during the colonic migrating motor complex in the isolated murine large intestine. Neurogastroenterol. Motil. 24, 769–778. doi: 10.1111/j.1365-2982.2012.01934.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Petrov, K. A., Girard, E., Nikitashina, A. D., Colasante, C., Bernard, V., Nurullin, L., et al. (2014). Schwann cells sense and control acetylcholine spillover at the neuromuscular junction by α7 nicotinic receptors and butyrylcholinesterase. J. Neurosci. 34, 11870–11883. doi: 10.1523/JNEUROSCI.0329-14.2014

PubMed Abstract | CrossRef Full Text | Google Scholar

Robitaille, R. (1998). Modulation of synaptic efficacy and synaptic depression by glial cells at the frog neuromuscular junction. Neuron 21, 847–855. doi: 10.1016/s0896-6273(00)80600-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, T. K., and Gershon, M. D. (2015). Cross talk proposal: 5-HT is necessary for peristalsis. J. Physiol. 593, 3225–3227. doi: 10.1113/jp270182

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, T. K., Park, K. J., and Hennig, G. W. (2015). Colonic migrating motor complexes, high amplitude propagating contractions, neural reflexes and the importance of neuronal and mucosal serotonin. J. Neurogastroenterol. Motil. 20, 423–446. doi: 10.5056/jnm14092

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, T. K., Spencer, N. J., Hennig, G. W., and Dickson, E. J. (2007). Recent advances in enteric neurobiology: mechanosensitive interneurons. Neurogastroenterol. Motil. 19, 869–878. doi: 10.1111/j.1365-2982.2007.01019.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Spencer, N. J., Hennig, G. W., Bayguinov, P., Park, K. J., Lee, H. T., Heredia, D. J., et al. (2006). Activation of neuronal circuitry and Ca2+ waves in longitudinal and circular muscle during colonic MMCs and the consequences of rectal aganglionosis in mice. Am. J. Physiol. Gastrointest. Liver Physiol. 292, 546–555. doi: 10.1152/ajpgi.00352.2006

PubMed Abstract | CrossRef Full Text

Spencer, N. J., Hennig, G. W., Dickson, E., and Smith, T. K. (2005). Synchronization of enteric neuronal firing during the murine colonic MMC. J. Physiol. 564, 829–847. doi: 10.1113/jphysiol.2005.083600

PubMed Abstract | CrossRef Full Text | Google Scholar

Spencer, N. J., Hennig, G. W., and Smith, T. K. (2003). Stretch-activated neuronal pathways to longitudinal and circular muscle in guinea pig distal colon. Am. J. Physiol. Gastrointest. Liver Physiol. 284, G231–G241. doi: 10.1152/ajpgi.00291.2002

PubMed Abstract | CrossRef Full Text | Google Scholar

Stevens, R. J., Publicover, N. G., and Smith, T. K. (1999). Induction and regulation of Ca2+ waves by enteric neural reflexes. Nature 399, 62–66. doi: 10.1038/19973

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, X. R., Badura, A., Pacheco, D. A., Lynch, L. A., Schneider, E. R., Taylor, M. P., et al. (2013). Fast GCaMPs for improved tracking of neuronal activity. Nat. Commun. 4:2170. doi: 10.1038/ncomms3170

PubMed Abstract | CrossRef Full Text | Google Scholar

Tian, L., Hires, S. A., Mao, T., Huber, D., Chiappe, M. E., Chalasani, S. H., et al. (2009). Imaging neural activity in worms, flies and mice with improved GCaMP calcium indicators. Nat. Methods 6, 875–881. doi: 10.1038/nmeth.1398

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomas, E. A., Sjövall, H., and Bornstein, J. C. (2004). Computational model of the migrating motor complex of the small intestine. Am. J. Physiol. Gastrointest. Liver Physiol. 286, G564–G572. doi: 10.1152/ajpgi.00369.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilms, C. D., and Häusser, M. (2009). Lighting up neural networks using a new generation of genetically encoded calcium sensors. Nat. Methods 6, 871–872. doi: 10.1038/nmeth1209-871

PubMed Abstract | CrossRef Full Text | Google Scholar

Yamada, Y., and Mikoshiba, K. (2012). Quantitative comparison of novel GCaMP-type genetically encoded Ca(2+) indicators in mammalian neurons. Front. Cell. Neurosci. 6:41. doi: 10.3389/fncel.2012.00041

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, Y., Chen, J. H., Li, H., Yang, Z., Du, X., Hong, L., et al. (2015). Involvement of 5-HT3 and 5-HT4 receptors in colonic motor patterns in rats. Neurogastroenterol. Motil. 27, 914–928. doi: 10.1111/nmo.12550

PubMed Abstract | CrossRef Full Text | Google Scholar

Zariwala, H. A., Borghuis, B. G., Hoogland, T. M., Madisen, L., Tian, L., De Zeeuw, C. I., et al. (2012). A Cre-dependent GCaMP3 reporter mouse for neuronal imaging in vivo. J. Neurosci. 32, 3131–3141. doi: 10.1523/JNEUROSCI.4469-11.2012

PubMed Abstract | CrossRef Full Text | Google Scholar