BioMed Research International
Volume 2016 (2016), Article ID 8758460, 13 pages
http://dx.doi.org/10.1155/2016/8758460

Research Article

Ultrastructural Mapping of the Zebrafish Gastrointestinal System as a Basis for Experimental Drug Studies

Delfine Cheng,¹ Gerald J. Shami,¹ Marco Morsch,² Roger S. Chung,² and Filip Braet1,3,4

¹School of Medical Sciences (Discipline of Anatomy and Histology), The Bosch Institute, The University of Sydney, Sydney, NSW 2006, Australia
²Faculty of Medicine and Health Sciences, Macquarie University, Sydney, NSW 2109, Australia
³Australian Centre for Microscopy & Microanalysis (ACMM), The University of Sydney, Sydney, NSW 2006, Australia
⁴Charles Perkins Centre, The University of Sydney, Sydney, NSW 2006, Australia

Received 17 March 2016; Accepted 4 May 2016

Academic Editor: Minjun Chen

Copyright © 2016 Delfine Cheng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Research in the field of gastroenterology is increasingly focused on the use of alternative nonrodent model organisms to provide new experimental tools to study chronic diseases. The zebrafish is a particularly valuable experimental platform to explore organ and cell structure-function relationships under relevant biological and pathobiological settings. This is due to its optical transparency and its close-to-human genetic makeup. To-date, the structure-function properties of the GIS of the zebrafish are relatively unexplored and limited to histology and fluorescent microscopy. Occasionally those studies include EM of a given subcellular process but lack the required full histological picture. In this work, we employed a novel combined biomolecular imaging approach in order to cross-correlate 3D ultrastructure over different length scales (optical-, X-ray micro-CT, and high-resolution EM). Our correlated imaging studies and subsequent data modelling provide to our knowledge the first detailed 3D picture of the zebrafish larvae GIS. Our results provide unequivocally a limit of confidence for studying various digestive disorders and drug delivery pathways in the zebrafish.

1. Introduction

Zebrafish (ZF, Danio rerio) exhibit a high degree of resemblance in their genetic profile (69% of the their genes have at least one human ortholog [1]), molecular mechanisms, cell development, and organ physiology to humans [2]. Their small size and translucent nature make them easy to manipulate and observe as a whole animal [3], which contribute to their attractiveness as model organisms for biological experimentation. Accordingly, ZF have been successfully employed to study various physiological and pathophysiological processes in embryogenesis, organogenesis, genomics, and cancerogenesis, as well as research in toxicology and drug studies (for reviews, see [4–6]).

Although ZF have been shown to be a beneficial complementary model to rodent in many research fields, they lack some of the typical mammalian organs such as lung, prostate, skin, and mammary glands. On the other hand, as vertebrates, ZF possess organs and tissues, such as heart, kidney, liver, pancreas, intestinal tract, and brain that display analogous structures and functions to those found in humans [7]. Furthermore, with the advances in genetics, any human tumour type has also been successfully genetically developed in the ZF with similar morphology, gene expression, and signalling pathways [8, 9]. Consequently, despite their infancy in research fields, ZF have already gained distinctive popularity in the fields of genetics and development [10, 11]: they have become a complementary platform to rodent experimental models in preclinical screening studies in the field of translational drug research, in particular, in the assessment of drug compounds delivered to the gastrointestinal system (GIS). These attributes further extended the use of ZF to include the fields of metabolic organ diseases, including cancer [4, 12–14].

The ZF digestive system organogenesis [15–18] and morphogenesis [19–21] have been described by a few, and histological information is available on FishNet [22], the Zebrafish Atlas [23], and The Zebrafish Atlas of Macroscopic and Microscopic Anatomy [24]. However, to-date, the literatures available regarding comparative imaging studies that explore the ZF GIS with different imaging modalities and fine ultrastructural studies of the ZF microarchitecture are particularly scarce, mainly limited to advanced light-laser optical microscopy approaches. Hence, we present here the first full image sets of a 12 dpf ZF larvae from the macro- to the nanometre scales, including detailed information on cellular and subcellular features of the digestive system organs and a comprehensive comparative analysis of ZF GIS, weighted against published literature from rodents.

2. Materials and Methods

Firstly, a dedicated sample preparation and imaging workflow was prepared after assessing different experimental approaches, allowing for the subsequent swift imaging of one sample across different microscopy platforms, including X-ray microcomputed tomography (micro-CT), light microscopy (LM), and electron microscopy (EM) (Figure 1, top panel). Secondly, following this workflow, we designed the subsequent sample manipulation processes so the entire sample was retained for whole-mount investigation, allowing for the ability to image multiple areas, multiple times, and across different beam-lines (Figure 1, lower panel).

Figure 1: Sample preparation and imaging workflow used for the observation and ultrastructural data correlation of a single zebrafish sample compatible with X-ray micro-CT, LM, and EM imaging modalities. LM imaging modality includes the array tomography technique, whereby serial sections are collected onto a glass slide and imaged using LM and back-scattered EM (BSEM). EM includes TEM, transmission electron tomography (TET), BSEM, and SBF-SEM. This sample preparation protocol not only allows for the sample to be compatible with all the different microscopy platforms but also provides superior ultrastructural preservation of the zebrafish larvae, compared to conventional protocols used for EM.

2.1. Zebrafish Animal Model

Zebrafish (wild-type, Danio rerio) were maintained at 28°C in a 13 h light and 11 h dark cycle. Embryos were collected by natural spawning and raised at 28.5°C in E3 solution according to standard protocols [25]. Note that, from previous organogenesis studies of the ZF digestive system, it was determined that from 6 days after fertilisation (dpf), when the yolk is exhausted, the digestive functions are comparable to those of an adult fish [15, 19, 21]. Therefore, we studied ZF larvae aged 12 dpf to map the ultrastructure of the GIS. This timeframe is of particular importance as we can take the advantages including the optical translucent properties of the animal and the small size of the fish to facilitate the different sample preparation steps needed, inherent to correlated biomolecular microscopy (i.e., staining and fixation). The sample can be processes as a whole animal (i.e., whole mount), which excludes dissection artefacts and is large enough to easily differentiate the internal organs in LM.

2.2. Sample Preparation

For all microscopy examination purposes, ZF larvae at 12 dpf were collected and fixed in 4% paraformaldehyde + 2.5% glutaraldehyde in cacodylate buffer (4% sucrose + 0.15 mM CaCl₂ in 0.1 M sodium cacodylate buffer), overnight at 4°C. The samples were prepared following a protocol modified from Deerinck et al. [26] and Tapia et al. [27], whereby whole animals were exposed to the following solutions of heavy metal stains and mordanting agents: 2% osmium tetroxide (OsO₄) + 1.5% potassium ferrocyanide (KFeCN) for 2 h at 4°C, 1% thiocarbohydrazide (TCH) for 20 min, 2% OsO₄ for 20 min in the dark, 1% aqueous uranyl acetate at 4°C overnight, and Walton’s lead aspartate for 30 min at 60°C. The samples were thoroughly washed with distilled water in between staining steps, then dehydrated through a series of ethanol, gradually infiltrated with a series of EPON (hard grade) dilutions over a period of 3 days, and incubated in pure resin over 2 days. Finally, each ZF was carefully positioned, head-down, at the bottom of a BEEM® capsule (bottleneck, size 00), filled with fresh resin, and polymerised at 60°C overnight.

2.3. X-Ray Microcomputed Tomography (Micro-CT)

Micro-CT imaging was performed on resin-embedded ZF samples using a Skyscan 1072 system (Bruker microCT, UK), with no filter, operating at 40 kV and over 180 deg rotation angle. The individual X-ray images were reconstructed using NRecon (Bruker microCT, UK) and volumetric data were processed and analysed using Avizo (FEI Software) and IMOD (Boulder, Colorado, USA).

2.4. Light Microscopy (LM) and Array Tomography

ZF larvae were orientated sagittal to minimise the block face surface. The head was trimmed off to just behind the eyes and the block was carefully faced up to just behind the inner ears or to the start of the pharyngeal pad, just prior to exposing the oesophagus and the liver. Excess resin around the tissue is trimmed away with a razor blade and a mix of Welwood® glue and xylene (1 : 2) was applied to the sides of the block, following the methods developed by Micheva and Smith [28] and Blumer et al. [29] for array tomography sectioning. We found that applying glue only to the bottom of the trapezoid was sufficient and preferable, as the dried layer of glue from the top and the sides frequently does not cut well and builds up and interferes with the sectioning. This method successfully generated long, uninterrupted serial sections with minimal section loss. A ultramicrotome (Ultracut 7, Leica Microsystems, Heerbrugg, Switzerland) and a histojumbo diamond knife (Diatome, USA) were used to create ribbons of 50 serial sections 0.5 μm thick, which were placed, in order, onto a glass slide previously placed inside the knife boat (Supplementary information 1 in Supplementary Material available online at http://dx.doi.org/10.1155/2016/8758460).

After every ribbon, the slide was dried, stained with 0.5% toluidine blue, and observed with LM for orientation purposes. Once the entire GIS had been sectioned, each section containing the GIS was sequentially imaged using a light microscope (DM6000, Leica, Germany). The ImageJ plugin StackReg (Biomedical Imaging Group, EPFL, Lausanne, Switzerland [30]) or the IMOD plugin midas (Boulder Lab for 3D Electron Microscopy of Cells, Colorado, USA) was used to either automatically or manually align the images relative to a reference image chosen within the stack (Supplementary information 2). The aligned stack was then used to create 3D models of the GIS organs using 3dMod (IMOD plugin), whereby volumetric and morphometric analysis can be generated. Counting sections from the array tomogram allows for precise localisation (within 0.5 μm) of the different organs, their size, and relative positions with each other.

2.5. Back-Scattered Scanning Electron Microscopy (BSEM)

Following LM imaging, the glass slides were carbon-coated and mounted on a stub in preparation for Scanning Electron Microscopy (SEM) imaging. Two lines of silver paint were applied from the top surface of the slides to the stubs to increase conductivity. Sections were then imaged by detection of back-scattered electrons using a SEM (Sigma VP FEG SEM, ZEISS, Germany) operating at 3.8 kV. Consecutive sections were imaged, aligned, and modelled following the same method described previously for array tomography with LM. In this way, dozens or even hundreds of consecutive sections can be imaged over multiple areas, multiple times.

2.6. Transmission Electron Microscopy (TEM)

At any time during sectioning for array tomography, ultrathin sections of 70 nm can be collected on 200 mesh copper grids for TEM observation. Typically, after every 50 sections, the ribbons of sections were checked. When an area of interest is apparent on the LM sections, an ultrathin section can be collected at this point which is adjacent to the last section from the last ribbon. Next, retrieved sections on grids were poststained with 2% aqueous uranyl acetate and Reynold’s lead citrate solutions for 10 min each and imaged with a TEM (1400 TEM JEOL, Tokyo, Japan), operating at 120 kV.

3. Results and Discussion

Currently, ZF have gained exponential momentum as an experimental animal model in biomedical research fields [31–33]. Although the ZF model has shown to be a viable addition to other animal models (e.g., pigs and dogs), and even rodents, they lack some of the typical mammalian organs as outlined earlier. Providentially, ZF possess a GIS that displays an analogous function to humans, including such organs as a liver, pancreas, and gallbladder, with the notable addition of a swim bladder [7]. This unique aspect of the ZF, together with their relative small size and translucent properties—directly benefiting sample processing and imaging—make them very appealing model in the fields of gastroenterology and hepatology. The importance of these features have been convincingly demonstrated throughout the literature, such as for liver development [34] and regeneration [35], hepatocellular carcinoma [36], hepatic nanoparticle-targeting [37], hypertriglyceridemia-mediated pancreatic organ abnormalities [38], drug-induced liver injury [39], intestinal inflammation [40], and gut-associated nutritional programming [41].

Despite the fervent application of ZF in the investigations of various functional aspects of liver and pancreas- and gut-associated diseases, their fine micro- and nanoanatomical structures are relatively unexamined. Further, there is no detailed comparison of ZF GIS to its mammalian and human counterparts. While many teams only apply state-of-the-art live-cell biomolecular optical imaging techniques, limited studies have subsequently verified their findings at the nanoscale, throughout multiple dimensions (i.e., , , and ). Indeed, high-resolution microscopy is the only way to provide important complementary histological information for many subcellular processes observed using optical microscopy [42–44].

In this contribution, we outline an alternative multimodal imaging route to fill in the missing histological pieces of the ZF GIS puzzle (Figure 1). By doing so, we firstly illustrated the unique overall ultrastructural resemblance of the ZF GIS to the mammalian GIS as examined by standard TEM (Figure 2). Next, we expanded those results by employing contemporary three-dimensional (3D) microscopy techniques (Figures 3–5) and finally reviewed our fine structure findings, including our morphometric data, against the existing literature (Table 1).

Table 1: Comparison between the zebrafish and the rat digestive systems, including the GIT and the accessory organs. The information on the ZF digestive system presented in the table below mainly result from the data analysis obtained in this study, complemented by information providing various literature sources [20, 24, 40, 54–56, 60–63]. The information on the rat ZF are mainly compiled from a review of the literature, including the following references [64–70].

(a)

Zebrafish liver
Zebrafish pancreas
Zebrafish gut

Size
12 dpf: volume = 5 061 840 µm³. Transverse diameter = 349 µm. Anterior-posterior diameter = 520 µm
Adult: volume = 0.535 mm³. Transverse diameter = 2.08 mm. Anterior-posterior diameter = 1.64 mm. (In males, ∼2.10% of total body weight. In females, ∼4.51% of total body weight)
Size
12 dpf: Volume = 1 738 670 µm³. Transverse diameter = 290 µm. Anterior-posterior diameter = 500 µm
Adult: volume = 0.3 mm³. Transverse diameter = 1.9 mm. Anterior-posterior diameter = 2.5 mm
Size
Long tube of about 400 µm. The interior is also folded and lined with a layer of simple columnar epithelium

Location
Posterior to inner ears and pharyngeal pad, anterior to intestinal bulb
Location
Between the intestine and the swim bladder. The head starts right posterior to the liver
Location
From esophagus to anus

Shape
3 lobes, boomerang shape. The left lobe, the largest, crosses the midline, under the swim bladder and the right lobe extends ventrally towards the head of the pancreas
Organisation
The pancreas is very diffused and acinar cells are scattered
Oesophagus
Located under the trachea. Short and muscular

Organisation of hepatocytes
Tubular
Endocrine tissue
One principal islet (or Brockman body) in the head (50 μm diameter). 2-3 secondary islets surround the principal islet
Composition
Mostly goblet cells and a few enterocytes

Composition
Hepatocytes, endothelial cells, bile duct epithelial cells. Kupffer cells seem to be absent
Beta cells
Insulin granules (160 nm)
Role
First place of enzymatic digestion

Hepatocytes
65% of total liver volume
Alpha cells
Glucagon granules (182 nm)
Intestinal bulb
Dilated, up to 80 µm in diameter

Size
Polygonal, 14–17 µm in size. Hepatocytes form plates, lined with sinusoids and biliary ducts
Other-cells
Delta-cells produce somatostatin
Composition
Enterocytes with very long microvilli at the beginning (up to 7.5 µm long and 115 nm in diameter) and then shortened to 2-3 µm long

Organelles
Nucleus, rER, mitochondria, Golgi apparatus, glycogen, lipid droplets, lysosomes
Vascular system
Rich vasculature (12% of the islet)
Role
Fat absorption

Biliary system
18% of liver volume
Exocrine tissue
Production of pancreatic digestive enzymes
Midintestine
Narrow, folded 3 times on itself and measures 120–140 µm in diameter

Bile canaliculi
Presence of microvilli on the surface. Diameter = 2.22 µm
Acinar cell
Polyhedral, 8–10 µm in size, surrounding a central acinar duct
Composition
More goblet cells than in the intestinal bulb. Enterocytes have shorter microvilli (2-3 µm long, 115 nm in diameter)

Bile ducts
Intrahepatic, and extrahepatic ducts
Organelles
Zymogen granules (2.30 µm), nucleus, ER, mitochondria
Role
Proteins absorption

Vascular system
17% of liver volume
Pancreatic system
Main pancreatic ducts through the middle of the pancreas head
Posterior intestine
Narrow and measures 80–90 µm in diameter

Endothelial cells
Fenestrated, average diameter = 130 nm.

Composition
Rare or no microvilli

Role
Osmoregulation

(b)

Rat liver
Rat pancreas
Rat gut

Size
Transverse diameter = 7.5–8 cm, superior-inferior diameter = 2.8–4.2 cm, anterior-posterior diameter = 2.2–2.5 cm. 5% of total body weight, mean weight = 13.6 g
Size
Weight = 804 mg
Size
The long tube measures up to 2 m in length

Location
Upper right portion of the abdomen, beneath the diaphragm and above the stomach. A small portion extends into the upper left quadrant
Location
In the cranial abdominal cavity, between the stomach and the small intestine
Location
From mouth to anus

Shape
4 lobes: left, middle, right, and caudate
Organisation
Diffused, lobulated in appearance and divided into 3 parts: biliary, duodenal, and gastrosplenic
Oesophagus
Long tube from mouth to stomach, approx. 20 cm long

Organisation of hepatocytes
Lobular
Endocrine tissue
2–3.5% of the pancreas. 5000 islets of Langerhans scattered within the organ
Composition
Stratified squamous epithelium

Composition
Hepatocytes, endothelial cells, bile duct epithelial cells, Kupffer cells
Beta cells
65–80% of the islet. Production of insulin granules (207 nm)
Role
Transport food from mouth to stomach

Hepatocytes
80.6% of total liver volume
Alpha cells
15–20% of the islet. Production of glucagon granules (172 nm)
Stomach
Functions are comparable to the ZF intestinal bulb, approx. 6 cm long

Size
Polyhedral, 27 µm in diameter. Overlapping plate-like sheets (trabeculae) form the three-dimensional structure of the liver lobule
Other-cells
Delta-cells (3–10% of the islet) produce somatostatin. PP-cells (3–5% of the islet) produce pancreatic polypeptides
Composition
Enterocytes are the main cell types (90% of the villi surface)

Organelles
Nucleus (25% are binucleate), mitochondria, rER, sER, lysosomes, Golgi apparatus, peroxisomes, lipid droplets, free ribosomes, lipoproteins, glycogen, polyosomes
Vascular system
Rich vasculature (10% of the islet)
Role
Storage of food, start of enzymatic digestion

Biliary system
0.2% of liver
Exocrine tissue
95% of the pancreas
Small intestine
Functions are comparable to the midintestine, approx. 1–1.5 m long

Bile canaliculi
Presence of microvilli on the surface. Diameter = 1.5 µm
Acinar cell
Polyhedral, 10 µm in diameter, surrounding a central acinar duct.
Composition
Divided into 3 parts: duodenum, jejunum, and ileum. Enterocytes are the main cell types

Bile ducts
Intrahepatic, interlobular, and extrahepatic bile ducts
Organelles
Zymogen granules (500–800 nm), nucleus, ER, mitochondria
Role
Absorption of nutrients

Vascular system
19.2% of liver volume
Pancreatic system
Anterior pancreatic duct (main duct) occupies 50–60% of the pancreas. Draining the pancreatic juice to the stomach and small intestine
Large intestine
Functions are comparable to the posterior intestine, approx. 22–26 cm long

Endothelial cells
Fenestrated, average diameter = 6.5 µm

Composition
Enterocytes

Role
Absorption of water and left over digested nutrients

Figure 2: Parasagittal section of a 12 dpf ZF larvae stained with toluidine blue and imaged with light microscopy, showing the different components of the digestive system (c). Corresponding EM images of the different regions include the oesophageal area, rich in goblet cells (Gob) (a), the pancreas, with a pancreatic duct (PD) in the middle and surrounded by acinar cells (Aci) rich in zymogen granules (b). The liver and its hepatocytes (Hep) surrounded by sinusoids (Sin) and its network of bile ducts are shown in (d), as well as the intestine lined with enterocytes (Ent) rich in villi forming the intestinal brush border (BB) in (e). (SB) is the swim bladder. Scale bar = 20 μm (LM) and 5 μm (EM).

Figure 3: Zebrafish larvae (12 dpf) digestive system imaged using X-ray, LM, and EM (BSEM and TEM). At any positions (here, sections 341, 376, and 431 are shown as examples), micro-CT images and model can be viewed as cross-sections. Corresponding LM images of toluidine blue stained sections (500 nm) can be retrieved by mean of measuring distances from recognisable organs in the X-ray data. Back-scattered SEM images are generated from the same sections as the LM sections. TEM images are generated from adjacent sections from the LM ones. Colour code for micro-CT model: swim bladder (yellow), pancreas (green), intestine (pink), and liver (blue). Colour code for EM images: hepatocyte (blue), islet of Langerhans (green), and intestinal brush border (pink). Scale bars = 100 μm (micro-CT) and 2 μm (TEM).

(a)

(b)

(c)

(d)

Figure 4: Zebrafish larvae (12 dpf) model of liver (blue) and its vasculature (red), generated by serial LM imaging of 416 consecutive sections of 500 nm. (a) Dorsal view of the liver. (b) Vasculature of the liver (17% of total volume) represented in (a). (c) Combined liver model and its vasculature. (d) Same as (c), viewed from a different angle and clipped opened to visualise the internal vasculature. Liver vasculature was modelled by thresholding the grey values corresponding to the vessels and sinusoids from individual LM images. For full animation, see Supplementary information. Scale bar = 50 μm.

Figure 5: Adult zebrafish digestive system reconstruction by X-ray micro-CT, showing the GIT (pink), liver (blue), pancreas (green), swim bladder (yellow), and oocytes (red). Cross-sections are shown on the bottom line for different positions (here, positions 285, 340, and 473 are used as examples). Scale bar = 20 mm.

The sample preparation protocol used in this study was modified from Deerinck et al. [26] and Tapia et al. [27] and was originally developed for the preparation of biological samples for serial-block face sectioning and back-scattered SEM imaging. Although the sample processing time is significantly longer (6-7 days) than that of a conventional sample preparation protocol for EM (3-4 days), we have still opted to base our protocol on the two protocols mentioned previously: the aforementioned series of heavy metals stains and mordanting agents have shown great contrast improvement and conductivity of biological tissues, critical for high-resolution EM imaging, without interfering with conventional LM staining. The use of reduced osmium (RO) with potassium ferrocyanide not only improved membrane preservation but also contrasted glycogen granules in hepatocytes [45, 46]. The addition of sucrose and CaCl₂ to the buffer further improved the ultrastructural preservation by stabilising the osmolarity [47] throughout the processing. X-ray imaging was used here as a relative quick way to image the GIS in its entirety. Due to the small size of the animal (in average 5 mm × 700 μm), a full scan of the GIS usually takes 2 hours while it would take a few days to scan an entire rat at a comparable resolution. Gross organ sizes and their relative positions within the fish body can be observed. The small sample size and the enhanced contrasting protocol used in this work also made the segmentation and modelling of the internal organs relatively easy. In fact, ZF samples processed with the protocol herein displayed superior contrast relative to samples processed using a standard protocol involving staining with osmium tetroxide only. In fact, previous studies have shown that soft tissue contrast can be enhanced for X-ray micro-CT imaging by heavy metal staining [48, 49]. Imaging of the sample embedded in resin also exhibits a higher signal-to-noise ratio compared to sample immersed in ethanol (Supplementary information 3). This may be related to the enhanced stability of the sample—increased signal—as well as the embedding media—decreased noise—during the scan acquisition. By imaging the entire intact ZF in the micro-CT, information such as 3D connectivity is retained and significant internal tissue damage resulting from sample preparation can be assessed. Since micro-CT is a nondestructive imaging technique, the same sample block can be retrieved and sectioned, and subsequent data can be correlated with histological (i.e., LM) and/or subcellular ultrastructural studies (i.e., EM) (Figure 3), as the sample was processed using a protocol compatible to all the mentioned imaging modalities. Sectioning for array tomography is the most time consuming step: it will take an experienced microtomist 6–8 hours to section through the entire GIS sagittal and collect the ribbons of sections on glass slides, following which the samples are conserved for future LM and EM imaging and analysis.

Micro-CT analysis allowed us to determine that the GIT of a 12 dpf ZF was approximately 2.5 mm long, from the oesophagus to the anus.

The liver, the largest of the digestive glands, plays a pivotal role in the maintenance of metabolic function and excretion [50]. Our micro-CT and array tomography reconstructions showed that the liver is a trilobe organ with a boomerang-like shape and lies ventrally and anterior to the swim bladder, surrounding the intestine (Figure 3). These observations are supported by previous histological observations [20]. In a 12 dpf larvae, we measured the liver to be approximately 25 μm (anterior-posterior diameter) and occupying an average volume of 5.10⁶ μm³ (approximately 2% of a 12 dpf ZF total body volume) within which 17% is occupied by blood vessels (Figure 4). Compared to a rat liver, a ZF liver volume is hence calculated to be 2 400 000 times smaller in volume and 2.2% less vascularised. From previous histological observations, the ZF liver is similar to that of other teleosts [51]: there are no portal triads nor hepatocytes arranged in plates but rather tubules of hepatocytes among portal veins, hepatic arteries, large biliary ducts, and sinusoids, which are distributed stochastically within the parenchyma [52]. Morphologically, the hepatocytes we analysed were polygonal, measured 17 μm in diameter (Figure 3), and formed tubules of small bile ducts, derived from the bile canaliculi, which was intrahepatic and had an average diameter of 2.2 μm. EM data showed vessels and sinusoids also lined with a monolayer of endothelial cells, exhibiting fenestrations varying from 80 to 230 nm in diameter, which means they can be up to double the size of the liver endothelial cells fenestrations reported in rats [53]. Notable was the apparent absence of Kupffer cells throughout the hepatic sinusoids, a distinct point of difference to mammalian livers.

Like in other animals who possess a pancreas, the ZF pancreas is a dichotomic organ and the site of glucohomeostasis [21, 54]. In the ZF, the pancreas is diffused and located around the liver and intestine: the exocrine part is mainly composed of pancreatic ducts and acinar cells. From our correlative imaging studies, the acinar cells were polyhedral, measured 8–10 μm in diameter, and were populated with relatively large zymogen granules (2–2.5 μm, compared to 500–800 nm found in rats) which discharge their contents into the pancreatic ducts to form the pancreatic juice. A principal islet of Langerhans (or Brockman body), identified in the pancreas head and measuring 50 μm in diameter, along with scattered secondary islets forms the endocrine component. The islets composition is as found in mammals, whereby a rich network of blood vessels surrounds neuroendocrine cells, namely, alpha-cells, beta cells, and delta-cells, which, respectively, produce the hormones glucagon, insulin, and somatostatin. Gamma and polypeptide producing- (PP-) cells appear to be absent (Figure 3) [24]. From our light microscopy analysis, the vasculature within the principal islet occupied approximately 12% of the islet. Glucagon and insulin granules were smaller than those found in the rat, measuring, respectively, 182 nm and 160 nm compared to 172 nm and 207 nm.

The gastrointestinal tract (GIT) is a long 3-fold tube, starting from the oesophagus to the anus, with morphological similarities to mammals: a simple columnar epithelium with folds (no real villi) including enterocytes, goblet cells, and endocrine cells. On the apical side, microvilli are present and the enterocytes are joined by tight junctions, which we determined by TEM to be 200–400 nm across. In the ZF, the oesophagus is short and muscular and is mainly composed of goblet cells and taste buds as it is the first site of enzymatic digestion. The intestine follows the oesophagus and histologists have divided it into 3 parts: the intestinal bulb, the midintestine, and the posterior intestine which can be differentiated by their shape, density of goblet cells, and the length of microvilli forming the brush border [18, 55]. Paneth cells, crypts, and organised lymphoid structures are absent across the entire intestine [55]. The intestinal bulb plays the role of food storage and fat absorption [24]. From our X-ray and LM observations, the “bulb” was dilated up to 80 μm wide; however, it could probably extend more following a large meal. Enterocytes were the main cell types and our analysis revealed that the enterocytes at the very beginning of the intestine have the longest microvilli (up to 7.5 μm long) after which they shorten to about 2-3 μm within the intestinal bulb (Figure 3). The midintestine is narrow and believed to be the site of protein absorption [19, 56]. Enterocytes in this area presented microvilli of similar length to those found in the intestinal bulb and were surrounded by more goblet cells than in the intestinal bulb. In the posterior intestine, goblet cells were still present but enterocytes presented short or few microvilli. In fact, most of the digestion process has already taken place and this part is mainly involved in osmoregulation [18].

Finally, we demonstrated that our whole-mount multimodal imaging approach could also be applicable to adult ZF (Figure 5). The versatile imaging workflow presented in this study validates the anatomical fine structures of the adult ZF GIS and compares the larvae ZF GIS systematically to the rodent experimental models (Table 1).

4. Conclusion

The sample preparation protocol presented here offered not only great versatility—compatible with X-ray, light, and electron microscopy—but also optimum ultrastructural preservation of ZF larvae as well as improved sample contrast and conductivity, required for high-resolution X-ray and EM imaging. The workflow can be easily adapted to incorporate in vivo fluorescence or fluorescently labelled structures by adding a fluorescent live imaging step prior to sample processing. The concept of Correlative Light and Electron Microscopy (CLEM) can then be applied to complement dynamic functional information with high-resolution ultrastructural details, on the same sample. The use of tracers or fiducial markers such as laser etching techniques [57] or carbon-coating of pattern on glass slides [58, 59] can facilitate the relocalisation of the area, cell,or subcellular structure of interest between the macro- and nanometre scales imaging. This workflow presents a real advantage in the fields of research aiming at exploring drug transport and xenobiotic metabolism within the digestive glands (e.g., liver and pancreas) or studying and modifying malignant cell behaviour via novel anticancer therapeutic approaches (e.g., colon, pancreas, and liver cancer) (Figure 6). In fact, the low cost, small size, and relative speed at which drugs can be tested in ZF have already made it a popular model for the aforementioned studies. By following the entire workflow, gross differences between experimental and control fish can be rapidly determined by X-ray micro-CT, while detailed analysis of the drug treatment upon pathology can be subsequently evaluated, on the same sample. While all imaging modalities might not be available within a certain research laboratory or institute, this workflow can be adapted to only incorporate the imaging instrumentation available. Since the sample has been prepared to accommodate all imaging techniques aforementioned, one can skip imaging modalities which are unavailable without affecting the next. Also, depending on the study, one can also limit the imaging techniques to only include the ones needed to answer the biological question.

Figure 6: Illustration showing the different routes of administration possible in zebrafish larvae to study the uptake, transport, metabolism, and efficacy of therapeutic drug- and/or cell-based approaches. (I) Indirect administration of complexes dissolved in water or administered directly via the oral route, mixed with food pellets. (II-III) Local and targeted microinjection of fluorescent macromolecular complexes at the site of interest or the use of microcapillary needles to deposit genetically modified cells within the digestive glands (II) or intestines (III). Note that those three administering routes are typically employed in routine preclinical screening studies in rodent models and human studies as well underpinning the relevance of the zebrafish model to investigate the pharmacology, toxicology, and effectiveness of new therapeutic interventions. Taking advantage of the optical translucent properties of the larvae, subsequent whole-mount live-cell imaging allows systematic monitoring of the treatment regimes using fluorescent navigation. The results can be combined with correlated electron microscopy techniques as depicted under Figures 1–4. Colour legend for the zebrafish: swim bladder (yellow); stomach and intestines (purple); liver (blue); pancreas (green).

Our comparative microscopy maps and concomitant delivery of 3D imaging models of the ZF digestive system provide a comprehensive overview of the ZF GIS and establish the much needed foundation as well as a limit of confidence in the use of ZF for future research in gastrodigestive related illnesses.

Competing Interests

The authors declare that they have no competing interests.

Authors’ Contributions

Delfine Cheng and Gerald J. Shami contributed equally.

Acknowledgments

The authors are grateful for the facilities and technical assistance from the Australian Microscopy & Microanalysis Research Facilities (AMMRF) at the Australian Centre for Microscopy & Microanalysis (ACMM) and the Advance Imaging Facility of the Bosch Institute, at The University of Sydney, and the Zebrafish facility at Macquarie University.

References

K. Howe, M. D. Clark, C. F. Torroja et al., “The zebrafish reference genome sequence and its relationship to the human genome,” Nature, vol. 496, no. 7446, pp. 498–503, 2013.

M. L. Kent, C. Harper, and J. C. Wolf, “Documented and potential research impacts of subclinical diseases in zebrafish,” ILAR Journal, vol. 53, no. 2, pp. 126–134, 2012.

C. Wittmann, M. Reischl, A. H. Shah, R. Mikut, U. Liebel, and C. Grabher, “Facilitating drug discovery: an automated high-content inflammation assay in zebrafish,” Journal of Visualized Experiments, no. 65, article e4203, 2012.

M. Gemberling, T. J. Bailey, D. R. Hyde, and K. D. Poss, “The zebrafish as a model for complex tissue regeneration,” Trends in Genetics, vol. 29, no. 11, pp. 611–620, 2013.

A. T. Nguyen, V. Koh, J. M. Spitsbergen, and Z. Gong, “Development of a conditional liver tumor model by mifepristone-inducible Cre recombination to control oncogenic kras V12 expression in transgenic zebrafish,” Scientific Reports, vol. 6, article 19559, 2016.

Y. Nishimura, A. Inoue, S. Sasagawa et al., “Using zebrafish in systems toxicology for developmental toxicity testing,” Congenital Anomalies, vol. 56, no. 1, pp. 18–27, 2016.

J. R. Goldsmith and C. Jobin, “Think small: zebrafish as a model system of human pathology,” Journal of Biomedicine and Biotechnology, vol. 2012, Article ID 817341, 12 pages, 2012.

J. F. Amatruda, J. L. Shepard, H. M. Stern, and L. I. Zon, “Zebrafish as a cancer model system,” Cancer Cell, vol. 1, no. 3, pp. 229–231, 2002.

S. Liu and S. D. Leach, “Zebrafish models for cancer,” Annual Review of Pathology: Mechanisms of Disease, vol. 6, no. 1, pp. 71–93, 2011.

G. J. Lieschke and P. D. Currie, “Animal models of human disease: zebrafish swim into view,” Nature Reviews Genetics, vol. 8, no. 5, pp. 353–367, 2007.

R. T. Peterson and C. A. MacRae, “Systematic approaches to toxicology in the zebrafish,” Annual Review of Pharmacology and Toxicology, vol. 52, no. 1, pp. 433–453, 2012.

I. A. Drummond, “Kidney development and disease in the zebrafish,” Journal of the American Society of Nephrology, vol. 16, no. 2, pp. 299–304, 2005.

A. V. Gore, K. Monzo, Y. R. Cha, W. Pan, and B. M. Weinstein, “Vascular development in the zebrafish,” Cold Spring Harbor Perspectives in Medicine, vol. 2, no. 5, Article ID a006684, 2012.

J. Wang, R. Karra, A. L. Dickson, and K. D. Poss, “Fibronectin is deposited by injury-activated epicardial cells and is necessary for zebrafish heart regeneration,” Developmental Biology, vol. 382, no. 2, pp. 427–435, 2013.

T. Tao and J. Peng, “Liver development in zebrafish (Danio rerio),” Journal of Genetics and Genomics, vol. 36, no. 6, pp. 325–334, 2009.

N. Tiso, E. Moro, and F. Argenton, “Zebrafish pancreas development,” Molecular and Cellular Endocrinology, vol. 312, no. 1-2, pp. 24–30, 2009.

E. A. Ober, H. A. Field, and D. Y. R. Stainier, “From endoderm formation to liver and pancreas development in zebrafish,” Mechanisms of Development, vol. 120, no. 1, pp. 5–18, 2003.

K. N. Wallace, S. Akhter, E. M. Smith, K. Lorent, and M. Pack, “Intestinal growth and differentiation in zebrafish,” Mechanisms of Development, vol. 122, no. 2, pp. 157–173, 2005.

A. N. Y. Ng, T. A. De Jong-Curtain, D. J. Mawdsley et al., “Formation of the digestive system in zebrafish: III. Intestinal epithelium morphogenesis,” Developmental Biology, vol. 286, no. 1, pp. 114–135, 2005.

H. A. Field, E. A. Ober, T. Roeser, and D. Y. R. Stainier, “Formation of the digestive system in zebrafish. I. Liver morphogenesis,” Developmental Biology, vol. 253, no. 2, pp. 279–290, 2003.

H. A. Field, P. D. Si Dong, D. Beis, and D. Y. R. Stainier, “Formation of the digestive system in zebrafish. II. Pancreas morphogenesis,” Developmental Biology, vol. 261, no. 1, pp. 197–208, 2003.

R. J. Bryson-Richardson, S. Berger, T. F. Schilling et al., “FishNet: an online database of zebrafish anatomy,” BMC Biology, vol. 5, article 34, 2007.

Jake Gittlen Cancer Research Foundation, Zebrafish Atlas, The Pennsylvania State University, 2013.

J. Schwendinger-Schreck, “The zebrafish: atlas of macroscopic and microscopic anatomy,” Yale Journal of Biology and Medicine, vol. 86, no. 3, p. 435, 2013.

M. Westerfield, The Zebrafish Book: A Guide for the Laboratory Use of Zebrafish Danio∗ (Brachydanio) Rerio, ZFIN, Eugene, Ore, USA, 2000, http://zfin.org/zf_info/zfbook/zfbk.html.

T. J. Deerinck, E. A. Bushong, A. Thor, and M. H. Ellisman, Protocol, University of California, San Diego, Calif, USA, 2010.

J. C. Tapia, N. Kasthuri, K. J. Hayworth et al., “High-contrast en bloc staining of neuronal tissue for field emission scanning electron microscopy,” Nature Protocols, vol. 7, no. 2, pp. 193–206, 2012.

K. D. Micheva and S. J. Smith, “Array tomography: a new tool for imaging the molecular architecture and ultrastructure of neural circuits,” Neuron, vol. 55, no. 1, pp. 25–36, 2007.

M. J. F. Blumer, P. Gahleitner, T. Narzt, C. Handl, and B. Ruthensteiner, “Ribbons of semithin sections: an advanced method with a new type of diamond knife,” Journal of Neuroscience Methods, vol. 120, no. 1, pp. 11–16, 2002.

P. Thévenaz, U. E. Ruttimann, and M. Unser, “A pyramid approach to subpixel registration based on intensity,” IEEE Transactions on Image Processing, vol. 7, no. 1, pp. 27–41, 1998.

C. A. MacRae and R. T. Peterson, “Zebrafish as tools for drug discovery,” Nature Reviews Drug Discovery, vol. 14, no. 10, pp. 721–731, 2015.

J. Barriuso, R. Nagaraju, and A. Hurlstone, “Zebrafish: a new companion for translational research in oncology,” Clinical Cancer Research, vol. 21, no. 5, pp. 969–975, 2015.

A. Schlegel and P. Gut, “Metabolic insights from zebrafish genetics, physiology, and chemical biology,” Cellular and Molecular Life Sciences, vol. 72, no. 12, pp. 2249–2260, 2015.

J. Chu and K. C. Sadler, “New school in liver development: lessons from zebrafish,” Hepatology, vol. 50, no. 5, pp. 1656–1663, 2009.

K. C. Sadler, K. N. Krahn, N. A. Gaur, and C. Ukomadu, “Liver growth in the embryo and during liver regeneration in zebrafish requires the cell cycle regulator, uhrf1,” Proceedings of the National Academy of Sciences of the United States of America, vol. 104, no. 5, pp. 1570–1575, 2007.

J.-W. Lu, “Zebrafish as a disease model for studying human hepatocellular carcinoma,” World Journal of Gastroenterology, vol. 21, no. 42, pp. 12042–12058, 2015.

L. Evensen, P. L. Johansen, G. Koster et al., “Zebrafish as a model system for characterization of nanoparticles against cancer,” Nanoscale, vol. 8, no. 2, pp. 862–877, 2016.

C. Liu, K. P. Gates, L. Fang et al., “Apoc2 loss-of-function zebrafish mutant as a genetic model of hyperlipidemia,” Disease Models and Mechanisms, vol. 8, no. 8, pp. 989–998, 2015.

P. Zhang, X. Jiang, X. Nie et al., “A two-photon fluorescent sensor revealing drug-induced liver injury via tracking γ-glutamyltranspeptidase (GGT) level in vivo,” Biomaterials, vol. 80, pp. 46–56, 2016.

S. Brugman, “The zebrafish as a model to study intestinal inflammation,” Developmental & Comparative Immunology, 2016.

E. Perera and M. Yúfera, “Soybean meal and soy protein concentrate in early diet elicit different nutritional programming effects on juvenile zebrafish,” Zebrafish, vol. 13, no. 1, pp. 61–69, 2016.

F. G. A. Faas, M. C. Avramut, B. M. van den Berg, A. Mieke Mommaas, A. J. Koster, and R. B. G. Ravelli, “Virtual nanoscopy: generation of ultra-large high resolution electron microscopy maps,” Journal of Cell Biology, vol. 198, no. 3, pp. 457–469, 2012.

S. J. Nixon, R. I. Webb, M. Floetenmeyer, N. Schieber, H. P. Lo, and R. G. Parton, “A single method for cryofixation and correlative light, electron microscopy and tomography of zebrafish embryos,” Traffic, vol. 10, no. 2, pp. 131–136, 2009.

J. R. Anderson, B. W. Jones, J.-H. Yang et al., “A computational framework for ultrastructural mapping of neural circuitry,” PLoS Biology, vol. 7, no. 3, Article ID e1000074, 2009.

D. L. White, J. E. Mazurkiewicz, and R. J. Barrnett, “A chemical mechanism for tissue staining by osmium tetroxide-ferrocyanide mixtures,” Journal of Histochemistry and Cytochemistry, vol. 27, no. 7, pp. 1084–1091, 1979.

E. Wisse, F. Braet, H. Duimel et al., “Fixation methods for electron microscopy of human and other liver,” World Journal of Gastroenterology, vol. 16, no. 23, pp. 2851–2866, 2010.

M. A. Hayat, “Methods of fixation,” in Fixation for Electron Microscopy, M. A. Hayat, Ed., chapter 7, pp. 200–208, Academic Press, New York, NY, USA, 1981.

C. C. Wong, I. S. Curthoys, and S. J. O'Leary, “Heavy metal staining, a comparative assessment of gadolinium chloride and osmium tetroxide for inner ear labyrinthine contrast enhancement using X-ray microtomography,” Acta Oto-Laryngologica, vol. 133, no. 1, pp. 22–27, 2012.

R. Mizutani and Y. Suzuki, “X-ray microtomography in biology,” Micron, vol. 43, no. 2-3, pp. 104–115, 2012.

S. R. Z. Abdel-Misih and M. Bloomston, “Liver anatomy,” Surgical Clinics of North America, vol. 90, no. 4, pp. 643–653, 2010.

R. C. Hardman, D. C. Volz, S. W. Kullman, and D. E. Hinton, “An in vivo look at vertebrate liver architecture: three-dimensional reconstructions from medaka (Oryzias latipes),” The Anatomical Record, vol. 290, no. 7, pp. 770–782, 2007.

A. L. Menke, J. M. Spitsbergen, A. P. M. Wolterbeek, and R. A. Woutersen, “Normal anatomy and histology of the adult zebrafish,” Toxicologic Pathology, vol. 39, no. 5, pp. 759–775, 2011.

F. Braet and E. Wisse, “Structural and functional aspects of liver sinusoidal endothelial cell fenestrae: a review,” Comparative Hepatology, vol. 1, article 1, 2002.

L. Gnügge, D. Meyer, and W. Driever, “Pancreas development in zebrafish,” Methods in Cell Biology, vol. 76, pp. 531–551, 2004.

K. N. Wallace and M. Pack, “Unique and conserved aspects of gut development in zebrafish,” Developmental Biology, vol. 255, no. 1, pp. 12–29, 2003.

A. J. Trotter, A. C. Parslow, and J. K. Heath, “Morphological analysis of the zebrafish digestive system,” in Zebrafish: Methods and Protocols, vol. 546 of Methods in Molecular Biology, pp. 289–315, Humana Press, New York, NY, USA, 2009.

I. Kolotuev, D. J. Bumbarger, M. Labouesse, and Y. Schwab, “Targeted ultramicrotomy: a valuable tool for correlated light and electron microscopy of small model organisms,” in Methods in Cell Biology, T. Müller-Reichert and P. Verkade, Eds., pp. 203–222, Academic Press, New York, NY, USA, 2012.

K.-I. Kobayashi, D. Cheng, M. Huynh, K. R. Ratinac, P. Thordarson, and F. Braet, “Imaging fluorescently labeled complexes by means of multidimensional correlative light and transmission electron microscopy: practical considerations,” in Correlative Light and Electron Microscopy, T. Müller-Reichert and P. Verkade, Eds., pp. 1–20, Elsevier, 2012.

P. Verkade, “Moving EM: the rapid transfer system as a new tool for correlative light and electron microscopy and high throughput for high-pressure freezing,” Journal of Microscopy, vol. 230, no. 2, pp. 317–328, 2008.

B. I. Kuperman and V. V. Kuz'mina, “The ultrastructure of the intestinal epithelium in fishes with different types of feeding,” Journal of Fish Biology, vol. 44, no. 2, pp. 181–193, 1994.

Y. Uchiyama and M. Watanabe, “Morphometric and fine-structural studies of rat pancreatic acinar cells during early postnatal life,” Cell and Tissue Research, vol. 237, no. 1, pp. 123–129, 1984.

M. Pack, L. Solnica-Krezel, J. Malicki et al., “Mutations affecting development of zebrafish digestive organs,” Development, vol. 123, no. 1, pp. 321–328, 1996.

N. G. Kan, D. Junghans, and J. C. I. Belmonte, “Compensatory growth mechanisms regulated by BMP and FGF signaling mediate liver regeneration in zebrafish after partial hepatectomy,” The FASEB Journal, vol. 23, no. 10, pp. 3516–3525, 2009.

A. A. Elayat, M. M. el-Naggar, and M. Tahir, “An immunocytochemical and morphometric study of the rat pancreatic islets,” Journal of Anatomy, vol. 186, part 3, pp. 629–637, 1995.

S. Nussey and S. Whitehead, “The endocrine pancreas,” in Endocrinology: An Integrated Approach, chapter 2, BIOS Scientific, Oxford, UK, 2001.

J. R. Henderson and M. C. Moss, “A morphometric study of the endocrine and exocrine capillaries of the pancreas,” Quarterly Journal of Experimental Physiology, vol. 70, no. 3, pp. 347–356, 1985.

P. M. Motta, G. Macchiarelli, S. A. Nottola, and S. Correr, “Histology of the exocrine pancreas,” Microscopy Research and Technique, vol. 37, no. 5-6, pp. 384–398, 1997.

J. Jo, M. Y. Choi, and D.-S. Koh, “Size distribution of mouse Langerhans islets,” Biophysical Journal, vol. 93, no. 8, pp. 2655–2666, 2007.

T. T. Kararli, “Comparison of the gastrointestinal anatomy, physiology, and biochemistry of humans and commonly used laboratory animals,” Biopharmaceutics and Drug Disposition, vol. 16, no. 5, pp. 351–380, 1995.

I. Nagy, Á. Pap, and V. Varró, “Time-course of changes in pancreatic size and enzyme composition in rats during starvation,” International Journal of Pancreatology, vol. 5, no. 1, pp. 35–45, 1989.