|
|
||||||||
Review
Experimental and Molecular Cardiology Group, Academic Medical Center, 1105 AZ, Amsterdam, The Netherlands
| ABSTRACT |
|---|
|
|
|---|
mouse; embryo; heart; in situ hybridization
| INTRODUCTION |
|---|
|
|
|---|
Three-dimensional (3D) reconstructions of embryonic development were originally based on computer-aided manual tracing of the organs of interest (11, 24), and the reconstructed organs were illustrated by medical artists, e.g., Ref. 27. With the advent of digital cameras, reconstruction methods based on digitized images have become commonplace. A full review of past reconstruction methods would do no justice to those researchers who had to implement their reconstruction protocols with the hardware and software available at that time. Recently, methods were published based on episcopic image capturing (6, 26). These episcopic methods, which acquire an image just before sectioning, enable one to accurately obtain aligned high-resolution digital stacks, based on fluorescence. Another new technique, optical projection tomography (OPT), allows the reconstruction of fluorescent and nonfluorescent whole mount specimens from transmission images acquired from different directions (20). Micro-computer tomography (µCT) and micro-magnetic resonance imaging (µMRI) can be used for 3D visualization of tissue density and water, respectively, at the microscopic level (2, 9). Still missing, however, is a protocol that is adequate for 3D reconstruction of gene-specific staining patterns in embryonic organs. Of the above mentioned methods, the OPT method is suitable for the visualization of such signals. However, the whole mount specimens themselves prevent reliable staining to occur, due to penetration problems of the staining agents, especially at later embryonic stages (20). Visualization of whole mount stained specimens with confocal microscopy (8) suffers from similar penetration problems and from decay of intensities with increasing depth (18). To accurately reconstruct patterns of gene expression, tissue sectioning remains required. With sections, the staining agent only needs to penetrate a few micrometers, which allows complete staining of the material. Another benefit of choosing sectioned material is that it enables the correlation of different labels on single or serial sections and the reconstruction of local quantitative differences in cell density and cell division (21).
The relative contributions of variance at the different levels in research show that most of the variation originates from the biological variation between specimens (7). Therefore, the expression "do more, less well" has been coined as the maxim for morphometric research. In accord with this maxim, this paper describes a 3D reconstruction protocol for patterns of gene expression. The 3D quality of the resulting reconstructions, compared with some other methods, can be considered "low," but it is sufficient for most biological goals. For the same cost, this protocol can be used to generate multiple reconstructions instead of one single "high"-quality reconstruction. Such an approach yields the necessary information on biological variation between specimens, within the time and financial constraints of the research project. In this paper the use of this protocol is illustrated by using it to map and reconstruct the patterns of a myocardium-specific in situ hybridization (ISH) staining of a series of embryonic mouse hearts ranging from embryonic day 8.5 to 14.5 (ED 8.5 to 14.5). A detailed morphological account will be presented elsewhere (Soufan et al., unpublished results).
| MATERIALS AND METHODS |
|---|
|
|
|---|
-MHC, ß-MHC, cTnI, SERCA, and MLC-2v (Table 1). For each embryonic day, the ISH mix composition was optimized to yield a homogeneous staining of all myocardial parts of the heart.
|
Image acquisition.
The stained sections were digitized with a peltier-cooled digital monochrome CCD camera (Photometrics; 1,317 x 1,035 pixels, 12 bit, linear over its whole dynamic range) attached to a Zeiss Axiophot microscope, with normal bright-field illumination. Color images were obtained by using a liquid crystal tunable filter (Micro*Color, RGB-MS; Cambridge Research & Instrumentation, http://www.cri-inc.com), placed as a slider in the microscope. By capturing the same image using a red, green, and blue exclusive filter, color images are obtained (RGB).
The images are immediately processed by a custom-written macro in a camera control and image-processing software program (PMIS Image Processing Software; GKR Computer Consulting, http://www.gkrcc.com), which applies a flat-field correction and reduces each channel from 12-bit to 8-bit, effectively storing a 24-bit color image (Fig. 1A) for each section (Fig. 1, "image acquisition"). Additionally, a 2 x 2 binning procedure is applied, which resizes the 1,317 x 1,035 images to 658 x 517 pixels. This size reduction also speeds up the entire process as each image is reduced from 4 MB to 1 MB. The resizing of the images only leads to a minor loss of quality as the in-plane resolution of the image, after reduction, corresponds to the resolution of the microscope objectives. Throughout this paper, a set of images obtained from sections belonging to the same embryo will be referred to as a "stack."
|
Thresholding.
Each individual image of a stack was then processed in Image-Pro Plus by a custom-written "dynamic threshold" macro. This dynamic threshold macro converts the 24-bit color images into binary images containing only the signal of interest (Fig. 1B), regardless of the between-section variation in background and signal intensities (Fig. 1, "automatic thresholding"). The dynamic threshold applied by this macro is based on an automatic analysis of particle and area profiles within the image. A detailed description and validation of this dynamic threshold procedure will be published separately (A. T. Soufan, J. Hagoort, and J. M. Ruijter, unpublished results). In short, the image is split into several classes based on the behavior of a parameter that is extracted from the image using the selected area and particle number at different threshold values. This results in a stack of aligned binary images displaying the myocardium, which is examined to evaluate the specific signal isolated by the macro. In sections in which the macro was not successful, the signal was isolated manually (Fig. 1, "manual thresholding"). If less than 80% of the sections were thresholded as intended, then the current stack was discarded and a new set of sections was prepared and digitized (Fig. 1, "x% success?").
Lumen labeling.
To label the lumen of the heart, the stack of binary images was again processed in Image-Pro Plus, using the custom-written "label lumen" macro (Fig. 1, "label lumen"). This macro consists of an automatic and an interactive part. First, enclosed structures are automatically labeled in a separate color, thus effectively labeling the cardiac lumen. Second, the macro enables the manual closure of "spaces" that were not fully enclosed by myocardium, but by nonmyocardial tissue of the heart. The user places caps to close these structures in each section. The lumen macro converts the binary images into 8-bit gray-scale images with four gray values in which a value of 0 (black) stands for the background, 127 for lumen, 150 for user-drawn caps, and 255 (white) for myocardium (Fig. 1C, inverted colors).
It should be noted that the labeling of the lumen is an optional step of the reconstruction protocol (Fig. 1, "fill lumen?"). Initially, this step was included because of our interest in the volumetrics of the developing embryonic mouse heart. In retrospect, the separate visualization of the lumen turned out to be of great importance for the understanding of the structural complexity.
Manual inspection and correction in 2D.
The stack of labeled images was then loaded into Amira and used to visually inspect and manually repair the images, e.g., damaged sections, dust particles, using Amiras limited, but sufficient, interactive editing tools. At this stage, sections that were beyond repair were deleted from the stack and replaced by (combinations of) adjacent sections. Effectively, this repair procedure is performed in 2D (Fig. 1, "2D correction"), but, using Amira, 3D information of adjacent images can be used to guide the process.
Labeling.
To visualize the stack in 3D using Amira, it is necessary to convert the 8-bit gray-scale stacks into another format, the so-called LabelField file, which is closely related to uniform scalar fields. This conversion was done using the LabelVoxel module of Amira. In this LabelVoxel module, three thresholds were placed (at 100, 130, and 180), effectively separating the four gray values defined in the lumen labeling section (Fig. 1, "labeling").
After the conversion, the labels are enhanced using tools contained in the segmentation editor. Besides manual editing, a smooth step in 3D is implemented to compensate for the stair-stepping and surface breakage present in the z-direction, which are inherent to sectioning. The result is a LabelField file containing the boundaries of each individual gray value in 3D.
Surface conversion.
The LabelField file containing the structures of interest was then converted into a 3D scalar field, using the generalized marching cubes (GMC) module of Amira (Fig. 1, "surface conversion"). The module computes a triangular mesh of polygons of the interfaces between different structures in a LabelField file. This is done by means of a GMC algorithm. The resulting surface file can be graphically drawn, e.g., rendered (Fig. 1D). This render is already adequate for judging the shape of the heart (Fig. 1F) or measuring the volume of the individual structures. The surface render can also be manually rotated and scaled. Upon inspection of the result in 3D, holes may be visible in the structures (Fig. 1, "holes?"). These are often the result of misalignment between sections in the z-direction or by warped sections. Holes that appear to be artifacts are manually closed by making the necessary alterations in the 3D LabelField file (Fig. 1, "3D correction").
Volume conversion.
To enable the user to virtually cut the reconstruction and obtain virtual cross sections, it is necessary to convert the surface file to a volume file, as a structure in the surface render is simply a hollow shell (Fig. 1D). The conversion of a surface file to a volume file is done using the TetraGen module, which creates a volumetric grid of tetrahedrons (Fig. 1, "volume conversion"). The resulting volume render can now be virtually sectioned correctly in any direction (Fig. 1E).
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
|
Sectioning vs. whole mount methods.
In literature, several methods for reconstructing embryos have been described (2, 6, 8, 9, 20, 22, 26, 28). These reconstruction methods can be separated into two groups: invasive (the tissue is serially sectioned and stained) and noninvasive (the embryos are stained whole mount). Comparing the invasive methods against the noninvasive methods, one has to conclude that the main advantage of the invasive method is that gene expression patterns can be determined unambiguously. Sectioning allows reliable staining of the material because the staining agent only needs to penetrate a couple of micrometers. Undoubtedly, whole mount methods using very small tissue samples also do not suffer from penetration problems, and indeed none are described to occur in embryos of ED 10.0 (8). However, the noninvasive methods suffer from penetration problems at later stages of development, thus creating an uncertainty in the visualization of the signal, as reported for OPT on an ED 11.5 mouse embryo (20).
With regard to the 3D quality, the noninvasive methods are undeniably superior to invasive methods. The main difference between the two groups of methods lies in the resolution and quality of a stack in the z-direction. The z-resolution in the noninvasive methods is, most often, identical to the xy-resolution (isotropic sets). In the invasive method, the z-resolution is directly linked to the section thickness. Thus the z-quality is impaired by the section thickness and reduced even further by the imperfect alignment and morphological deformation introduced by sectioning. The latter is not entirely true for the episcopic image-capturing methods (6, 26), where image acquisition is performed before sectioning of the material. However, because the staining occurs before cutting, a tradeoff is made between the high z-quality and the poor penetration of the stain. On the other hand, the xy-resolution of the invasive methods is still higher than most of the whole mount methods. It is true that this does not compensate for the low z-resolution (thus 3D quality), but the high xy-resolution is a requirement for many molecular and embryological studies, e.g., studies that aim to assess the relationship between cell density and division (21) and patterns of gene expression during development.
Required resources.
Most papers do not discuss all the required resources, i.e., hardware, software, and time, in much detail. The hardware requirements for the protocol presented in this study are very lenient. No custom-made hardware has been used: one only needs to be able to digitize stained sections and have access to a decent desktop computer. Other methods, such as µMRI, µCT, and OPT (2, 9, 20, 28) require highly specialized and dedicated equipment and software. The software required for the current reconstruction protocol is also not very "heavy." The two software programs used were designed for end users (e.g., they do not require programming skills) and can be purchased at the sites described in the MATERIALS AND METHODS. The custom-written macros to automate the "thresholding" and "label lumen" steps of the protocol are available on request.
Software.
The initial protocol made use of five programs, which was considered unwieldy as it required several conversions of image formats and made the digital reconstruction unnecessarily time-consuming. To alleviate this, a reduction in the number of programs was required. A variety of programs was tested to obtain a slim, yet versatile, package of software, that still contained all the required visualization and image-processing tools. The final package consists of only two programs: Amira and Image-Pro Plus. Both programs have a solid base of developers, and updates are being released regularly. To further decrease the required amount of time, we automated as many steps as possible. This was successfully done for two steps: the thresholding step and the labeling of the lumen. The thresholding macro, used for converting color images (24-bit) to binary images (1-bit) by placing a dynamic threshold (Fig. 1, "automatic thresholding"), operated as intended for most stacks. Although the macro was written to compensate for biological and histological variation between sections (images) and between embryos (stacks), an average of 1 in 20 images still needed intervention by the user. A plethora of image "traits" could result in faulty automatic thresholding. For instance, the atria could display an overall lower signal than the ventricles, which resulted in the signal in the atria being identified as background and omitted from the binary image. The "label lumen" macro, included as an optional step in the protocol, also worked for most images. Enclosed structures were automatically identified and given a separate color to discriminate them from background and myocardium. Falsely labeled areas did arise in the outflow region in between the large arteries and in those areas of the heart were the "outside" is enclosed by myocardium, as might occur in the atrioventricular sulcus region. These false lumen areas had to be removed manually.
Accessibility of the protocol.
Initially, a reconstruction took about 7 days to complete. The implementation of a standard protocol, custom macros, and the revision of the software package reduced the time required for a reconstruction to 23 days. In both cases the histological steps (also 23 days) are not taken into account. However, reconstruction is riddled with pitfalls caused by the amalgamation of biology with computer technology. To smoothly deal with these pitfalls, experts from both fields should be close at hand during the reconstruction process.
After more than 20 reconstructions, a 3D reconstruction protocol was developed that demands few resources and is suitable for the visualization of gene expression in embryonic organs. All that the reconstruction method described herein requires is a stack of stained sections as input, an image-capturing device, a personal computer, two software programs, and two downloadable macros.
Reconstruction quality and biological variation.
As mentioned above, the resolution in the z-direction of the reconstructions obtained by using the protocol is "low" compared with some other methods. This is mainly caused by morphological deformation, which is inherent to sectioning. Correction of deformation using external marker-based automatic congruencing (22) has not been attempted, as it is almost impossible to place external markers close enough to the embryonic heart. Although it is possible to manually correct for deformation and to obtain a technically "perfect" reconstruction (25), this strategy has not been implemented in the current protocol, because this would not result in a biologically better reconstruction. Rather than investing the available amount of time and money in a single technically perfect reconstruction, we opted for multiple reconstructions. This chosen scheme would still give us the required information with regard to the morphology, but, most importantly, also insight into the biological variation present in the developing heart. A technically perfect reconstruction, however, would give us no insight into the variation present between different biological specimens. This rationale is in accordance with the recommendation "do more, less well" that is given for morphometric measurements (7). The consequence of this choice is that definite statements with regard to cellular details cannot be made based on the resulting reconstructions. The quality of the reconstructions (Fig. 2), however, is still sufficient to fulfill most biological goals. When one wishes to discuss details, one has to refer to the original sections. The reconstruction then allows the morphologist to place these details in their proper spatial context. Also, the reconstructions can be used to locate the proper position of individual sections from other experiments (Fig. 3). To this purpose, a complete series of reconstructed hearts (ED 8.5 to ED 16.5; 23 reconstructions per time point) will be made available. For details, E-mail to bioinfo{at}amc.uva.nl, subject: 3Dmouseheart.
|
|
|
|
The correlation of gene expression data and quantitative morphological changes, both in their proper morphological context, will then serve to gain insight in the functional genomics of heart development.
| ACKNOWLEDGMENTS |
|---|
We were supported by Netherlands Heart Foundation Grant M96.002.
| FOOTNOTES |
|---|
Article published online before print. See web site for date of publication (http://physiolgenomics.physiology.org).
Address for reprint requests and other correspondence: J. M. Ruijter, Dept. of Anatomy and Embryology, Academic Medical Center, Meibergdreef 15, Rm. K2-283, 1105 AZ Amsterdam, the Netherlands (E-mail: j.m.ruijter{at}amc.uva.nl).
10.1152/physiolgenomics.00182.2002.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
B. van Wijk, G. van den Berg, R. Abu-Issa, P. Barnett, S. van der Velden, M. Schmidt, J. M. Ruijter, M. L. Kirby, A. F.M. Moorman, and M. J.B. van den Hoff Epicardium and Myocardium Separate From a Common Precursor Pool by Crosstalk Between Bone Morphogenetic Protein- and Fibroblast Growth Factor-Signaling Pathways Circ. Res., August 28, 2009; 105(5): 431 - 441. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Schlueter and T. Brand A right-sided pathway involving FGF8/Snai1 controls asymmetric development of the proepicardium in the chick embryo PNAS, May 5, 2009; 106(18): 7485 - 7490. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. E. Mangoni and J. Nargeot Genesis and Regulation of the Heart Automaticity Physiol Rev, July 1, 2008; 88(3): 919 - 982. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. T.M. Mommersteeg, W. M.H. Hoogaars, O. W.J. Prall, C. de Gier-de Vries, C. Wiese, D. E.W. Clout, V. E. Papaioannou, N. A. Brown, R. P. Harvey, A. F.M. Moorman, et al. Molecular Pathway for the Localized Formation of the Sinoatrial Node Circ. Res., February 16, 2007; 100(3): 354 - 362. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. T. Soufan, G. van den Berg, J. M. Ruijter, P. A.J. de Boer, M. J.B. van den Hoff, and A. F.M. Moorman Regionalized Sequence of Myocardial Cell Growth and Proliferation Characterizes Early Chamber Formation Circ. Res., September 1, 2006; 99(5): 545 - 552. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. T.M. Mommersteeg, A. T. Soufan, F. J. de Lange, M. J.B. van den Hoff, R. H. Anderson, V. M. Christoffels, and A. F.M. Moorman Two Distinct Pools of Mesenchyme Contribute to the Development of the Atrial Septum Circ. Res., August 18, 2006; 99(4): 351 - 353. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. M. Christoffels, M. T.M. Mommersteeg, M.-O. Trowe, O. W.J. Prall, C. de Gier-de Vries, A. T. Soufan, M. Bussen, K. Schuster-Gossler, R. P. Harvey, A. F.M. Moorman, et al. Formation of the Venous Pole of the Heart From an Nkx2-5-Negative Precursor Population Requires Tbx18 Circ. Res., June 23, 2006; 98(12): 1555 - 1563. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. T. Soufan, M. J.B. van den Hoff, J. M. Ruijter, P. A.J. de Boer, J. Hagoort, S. Webb, R. H. Anderson, and A. F.M. Moorman Reconstruction of the Patterns of Gene Expression in the Developing Mouse Heart Reveals an Architectural Arrangement That Facilitates the Understanding of Atrial Malformations and Arrhythmias Circ. Res., December 10, 2004; 95(12): 1207 - 1215. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. J. de Lange, A. F.M. Moorman, R. H. Anderson, J. Manner, A. T. Soufan, C. d. G.-d. Vries, M. D. Schneider, S. Webb, M. J.B. van den Hoff, and V. M. Christoffels Lineage and Morphogenetic Analysis of the Cardiac Valves Circ. Res., September 17, 2004; 95(6): 645 - 654. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. M.H Hoogaars, A. Tessari, A. F.M Moorman, P. A.J de Boer, J. Hagoort, A. T Soufan, M. Campione, and V. M Christoffels The transcriptional repressor Tbx3 delineates the developing central conduction system of the heart Cardiovasc Res, June 1, 2004; 62(3): 489 - 499. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Glueck Note from the Deputy Editor: reviews in mouse CV phenotyping Physiol Genomics, May 13, 2003; 13(3): 183 - 183. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |