<title>Introduction</title> <p>The ease of access to the embryo and subsequent manipulability has made the chick a reliable and powerful system for developmental biology. The power of this system has been enhanced recently by the availability of genomic data from both whole genome sequencing ##REF##15592404##(2004)## and large-scale EST projects (<xref rid="bib6 bib10 bib19 bib27" ref-type="bibr">Boardman et al., 2002; Carre et al., 2006; Hubbard et al., 2005; Kim et al., 2006</xref>). This opens up new opportunities for identifying all the genes that mediate the development of the embryo and its constituent parts and then using high throughput methods to test their function (##REF##12560806##Brown et al., 2003##). One of the first steps in this process is to document where and when specific genes are expressed in the embryo and a large amount of gene expression data is being generated. A repository of chick embryo gene expression data, GEISHA, has already been pioneered by Antin and colleagues (##REF##14991723##Bell et al., 2004##). This consists of a collection of photographs of embryos in which gene expression has been assayed mainly by using whole mount <italic>in-situ</italic> hybridisation, in some cases accompanied by sectioned material. The database that we describe here is different in that gene expression is visualised in 3D using OPT and expression patterns are mapped onto 3D digitised images. Here we describe how we have started to establish such a database of 3D gene expression patterns for the developing chick wing and investigated some of the practicalities involved.</p> <p>OPT was developed at the MRC Human Genetics Unit in Edinburgh and is one of a number of new microscopy techniques that have been developed in the last few years that allow capture of 3D image data. OPT has already been used to study the development of human, mouse, fly and plant embryos (<xref rid="bib13 bib26 bib28 bib34 bib40" ref-type="bibr">DeLaurier et al., 2006; Kerwin et al., 2004; Lee et al., 2006; McGurk et al., 2007; Sharpe et al., 2002</xref>) and one of its advantages is that it captures the three dimensional distribution of gene expression in an intact embryo.</p> <p>In order to compare large numbers of gene expression patterns, a number of recent atlas projects have taken the approach of mapping gene expression data to digital reference models. For example, the Edinburgh Mouse Atlas of Gene Expression (EMAGE) deals with 2D data in this way (##REF##15043218##Baldock et al., 2003##), section <italic>in-situ</italic> derived gene expression data in the mouse brain have been mapped to 3D models generated by Micromagnetic resonance imaging (##REF##17151600##Lein et al., 2007##) and a Zebrafish 3D anatomical Atlas has been produced based on sectioned material for the projection of gene expression data (##UREF##3##Verbeek et al., 1999##). Projects such as EMAGE (##REF##15043218##Baldock et al., 2003##) and GENEPAINT (##REF##14681479##Visel et al., 2004##) have begun to build large queryable databases containing both whole mount and section <italic>in-situ</italic> data for mouse embryos. Since the system has already been set up for mouse embryos, the establishment of a parallel database for the chick should allow direct comparisons between gene expression patterns in the two organisms. We have adopted the Bookstein thin plate spine algorithm for mapping our 3D data, which has previously been used extensively in morphometric analysis (<xref rid="bib1 bib9 bib17 bib45" ref-type="bibr">Albertson and Kocher, 2001; Bruner et al., 2004; Harmon, 2007; Yeh, 2002</xref>). Once gene expression data are assembled in a digital atlas, powerful modern data mining techniques can be used to examine relationships potentially leading to unexpected discoveries.</p> <p>We have focussed on the developing chick wing bud. The wing bud is a good test system for investigating the power of a 3D database because it is a structure with no significant morphological features at early stages. Many insights into the mechanisms that pattern the vertebrate limb have come from studies on chick embryos (##REF##15296968##Tickle, 2004##). In the long term, the ability to compare multiple patterns of gene expression should enable us to identify synexpression domains, complementary patterns and possibly also discrete boundaries of gene expression. We have optimised protocols to maximise consistency of initial whole mount <italic>in-situ</italic> hybridisation, OPT capture and mapping data to a reference model. In this paper, we have studied a number of previously described genes in terms of their expression patterns. This has allowed us to perform a pairwise analysis of overlap of expression for an initial set of genes and to identify some previously unappreciated features of expression with respect to dorso-ventral distribution. We have also used computational techniques based on microarray analysis to look for specific regions of the limb where genes are co-expressed.</p>
<title>Materials and methods</title> <title>Embryo preparation</title> <p>White Leghorn chick eggs were incubated in a humidified incubator at 38 °C for the appropriate time for the desired developmental stage as determined by the Hamburger and Hamilton stages (##REF##24539719##Hamburger and Hamilton, 1951##). Eggs were then windowed, embryos removed to ice-cold Phosphate Buffered Saline (PBS) (0.02 M phosphate, 0.15 M NaCl) and cleaned of extra-embryonic membranes. Eyes and forebrain were punctured with a tungsten needle to reduce trapping, and embryos were transferred to 4% (w/v) ice-cold paraformaldehyde (PFA) overnight. The embryos were then put through a graded methanol series at 4 °C; ending in 2, 100% methanol washes and stored at − 20 °C until use.</p> <title>Plasmid preparation and probe synthesis</title> <p>The plasmids used for the different genes were <italic>Shh</italic> (##REF##7916661##Echelard et al., 1993##), <italic>Fgf8</italic> (##REF##8548816##Crossley et al., 1996##), <italic>Msx1</italic> (##REF##2565278##Hill et al., 1989##), and <italic>Tbx3</italic> (##REF##9550719##Isaac et al., 1998##). All plasmids were linearised and transcribed according to their sources. EST clones acquired from ARK genomics were used as probes for <italic>Wnt3a</italic> (EST clone 603102629F1), <italic>Wnt5a</italic> (EST clone 603799237F1), <italic>Lmx1</italic> (EST clone 603127966F1) and <italic>HoxD13</italic> (EST clone 603499362F1). All EST clones were in pBluescript II KS+, which was linearised with Not1 (NEB) and transcribed with T3 RNA polymerase (Roche) to produce antisense probes. Plasmids were grown up using standard protocols and purified using Qiagen plasmid mini kits and individual clones were sequenced to check their identity.</p> <p>The RNA probes were synthesised accordance with standard protocols (<xref rid="bib33 bib37" ref-type="bibr">Maniatis et al., 1982; Nieto et al., 1996</xref>) and purified using the ProbeQuant G-50 spin column system (Amersham Biosciences). In some cases probe purification was performed using phenol chloroform extraction and Lithium Chloride precipitation as detailed in ##REF##8722478##Nieto et al. (1996)##.</p> <title><italic>In-situ</italic> hybridisation</title> <p>The <italic>in-situ</italic> protocol used was a modified version of that of ##REF##8722478##Nieto et al. (1996)##, full details of the modified protocol are in <xref rid="app1" ref-type="sec">supplementary materials</xref>.</p> <p>Before scanning under UV, embryos require some further washes to remove excess NBT–BCIP. Embryos were washed twice for 10 min in PBS at RT. Embryos were then moved to 10× TBST and allowed to equilibrate at RT, this should take between 10 and 20 min depending on the size of the embryo. Embryos were then washed 3 times for 20 min in 1× TBST and left to wash overnight in fresh 1× TBST at 4 °C. Embryos were washed 3 times for 5 min in PBT at RT and then fixed overnight in 4%PFA–DEPC–PBS at 4 °C. Embryos were washed briefly a further 2 times in PBS and then refixed in formal saline.</p> <title>Section <italic>in-situ</italic> hybridisation</title> <p>Section <italic>in-situs</italic> were performed on HH22 embryonic limbs according to the method of ##REF##11118473##Moorman et al. (2001)##. The <italic>Wnt5a</italic> probe was the same as above.</p> <title>Tyramide Signal Amplification (TSA)</title> <p>To modify the <italic>in-situ</italic> protocol for the fluorescent TSA colour reaction, glutaraldehyde was removed in the fixation step to minimise autofluorescence. Secondary detection was performed using a peroxidase (POD) linked anti-DIG antibody. The colour reaction was performed according to the Alexa Fluor 568 kit manufacturer’s instructions (Invitrogen — T-20914) with volumes increased to accommodate whole embryos.</p> <title>OPT sample preparation and scanning</title> <p>Standard reference embryos were fixed in 4% PFA/0.2% glutaraldehyde mix, which produces a stronger autofluorescent signal than PFA alone. The addition of 0.2% glutaraldehyde to the fix was not necessary for embryos that had been <italic>in-situ</italic> hydbridised, due to the presence of glutaraldehyde in the fixative steps of that protocol. Reference embryos were stored in 100% methanol until scanning, at which point they were taken back through a methanol series to PBS and briefly to water. Embryos having undergone <italic>in-situ</italic> hybridisation were washed 3 times for 20 min in PBS to remove storage fixative. In order to remove excess salts, embryos were washed twice for 10 min in distilled water and subsequently left overnight in distilled water followed by 1 wash of 10 min in fresh distilled water. OPT scanning was carried out following the protocol set out in ##REF##11964482##Sharpe et al. (2002)##, for more detailed protocols see <xref rid="app1" ref-type="sec">supplementary materials</xref>. The resulting output was in the Wlz file format which the Edinburgh Mouse atlas Project (EMAP) and EMAGE projects have utilised and which can store 3D grey scale information. The MRC HGU has produced a set of software tools for manipulating data in this format and these tools were used to convert the data into a format that could be imported in to the AMIRA software package.</p> <title>3D mapping</title> <p>The mapping of the 3D gene expression data to the reference models was performed using the AMIRA 4.1 software from Mercury Computer Systems. The data to be mapped were first roughly aligned with the reference model. Two corresponding sets of landmarks were then set up between an isosurface for the reference embryo and an isosurface for the fluorescent/anatomical data from the scan to be mapped. The landmarks were based on prominent morphological landmarks such as the AER, the region where the limb attaches to the flank and to proportional distances along the main axes of the limb. The fluorescent/anatomical data was warped, using a Bookstein thin plate spline method (##UREF##0##Bookstein, 1989##) provided by the AMIRA software and based on the previously defined landmark sets. Provided the resulting warped fluorescent/anatomical data seemed consistent with the reference limb’s morphology the same warp was then applied to the brightfield channel data. For full details see methods in <xref rid="app1" ref-type="sec">supplementary materials</xref>.</p> <title>Real-time PCR</title> <p>Chick embryos (incubated for 4 days at 38 °C) <italic>i.e.</italic>, approx. stage 22–23, were harvested in ice-cold PBS, the limb buds removed and the distal third cut off with tungsten needles. These pieces were immediately transferred to RNALater (Qiagen). A similar procedure was carried out on the proximal and median thirds. RNA preps were made of pooled limb sections with 20 sections for each region using the Qiagen RNA Easy micro kit, and checked using an Agilent bioanalyzer, using their RNA 6000 nanochip. The integrity values of all RNA samples were between 9.9 and 10 and the 28S and 18S ratio between 1.9 and 2.1. These values indicate little degradation or contamination.</p> <p>Real-time PCR was carried out using an Applied BioSystems HT-7900 machine in a manner similar to ##REF##17148754##Jesmin et al. (2007)## using the FastStart Taqman Probe master (Rox) standard reaction mix (Roche). Primers were selected from the Roche Universal Library using their online software. For the <italic>Wnt5a</italic> reaction, probe 52 was used, and for <italic>β-actin</italic>, probe 43 was used. Here chick probes are not automatically checked against other possible hybridisation targets, so this was carried out manually by Blasting the candidate sequences against the chick genome in Ensembl. The primers used for the <italic>Wnt5a</italic> reaction were: forward 5′catgatgaacctacacaatga 3′; reverse 5′ ccacgtcagccaggttgta 3′. And for the <italic>β-actin</italic> reaction were: F 5′ cacacaagtgcccatttacga 3′; R 5′ caagtccagacgcaggatg 3′. For full protocol see <xref rid="app1" ref-type="sec">supplementary materials</xref>.</p> <title>Computational analysis</title> <p>Simple arithmetical analyses were produced using either the AMIRA software package or Wlz based software tools. AMIRA’s arithmetic module was used to produce averages of multiple datasets and to produce masked datasets for domains of coexpression. Wlz based software developed by the MRC HGU was used to derive medians from multiple datasets and also to derive mean grey level intensity values for both discrete domains and serial sections through the limb. For fuller details of these image manipulations see <xref rid="app1" ref-type="sec">supplementary materials and methods</xref>.</p> <p>For more complex computational analyses each of the experimental gene-expression spatial distributions has been mapped into the standard coordinate framework defined by the model limb. To analyse the gene-expression patterns we first divided the limb into 560 non-overlapping sub-regions each of 10 × 10 × 10 voxels. Each of these was used to sample the experimental gene-expression patterns. For each experimental pattern, the mean gene-expression strength (integrated optical density divided by the volume) within each box was calculated. If the box was partially external to the limb then only the intersecting volume was considered. By this means a 2D matrix of mean expression strengths across the limb was calculated. Each row of the matrix for a given gene is a low-resolution representation of the pattern and each column for a given sample-region is the genetic “signature” for that spatial location.</p> <p>The resulting matrix of gene expression values was analysed using the TMEV4 package from TIGR. The data were analysed using a hierarchical clustering method (##REF##9843981##Eisen et al., 1998##) to produce a nested tree of gene expression pattern similarity based on a Euclidean distance metric. A nested tree was also produced of the similarity of the individual sample regions of the limb making up the 3D data model based on gene expression. The resulting tree was then used to identify clusters of similar expression made of small groups of regions at the terminus of long branches. These regions were used to produce larger 3D domains corresponding to the whole volume occupied by the regions comprising each cluster, which were subsequently visualised using the AMIRA software package.</p>
<title>Results and discussion</title> <title>Assessment of efficiency of whole mount <italic>in-situ</italic> hybridisation and scanning with OPT as a method of detecting gene expression domains</title> <p>An initial technical issue we encountered with OPT scanning of WISH (whole mount <italic>in-situ</italic> hybridisation) specimens was that strong <italic>in-situ</italic> colour reaction staining can block autofluorescence and prevent the capture of portions of the anatomical data required for subsequent mapping to a reference model. To obviate this problem, we identified a particular depth of staining with the NBT–BCIP substrate, suitable for OPT scanning, which captures an extensive range of the expression pattern and allows the visualisation of dynamic gradients without causing a substantial dropout of the anatomical data necessary for mapping (see <xref rid="app1" ref-type="sec">supplementary materials Fig. S</xref>1). To test our standard <italic>in-situ</italic> hybridisation protocol, monitor probe penetration of the embryonic limb, and assess the ability of the OPT system to identify graded patterns of expression within deep tissues, we focussed on <italic>Wnt5a</italic>. <italic>Wnt5a</italic> has been reported to have a proximo-distal gradient of expression based on both radioactive section <italic>in-situ</italic> hybridisation and northern blots of distinct portions of the limb (##REF##8297789##Dealy et al., 1993##) and is known to be expressed throughout mesodermal regions of the limb.</p> <p>We first assayed expression by WISH (method modified after ##REF##8722478##Nieto et al., 1996## see <xref rid="app1" ref-type="sec">supplementary data</xref>) ##FIG##0##Fig. 1##A). The <italic>Wnt5a</italic> whole mount was scanned using OPT and gene expression data mapped onto a reference limb (##FIG##1##Fig. 2##D), from which virtual sections were derived (##FIG##0##Fig. 1##B). These virtual sections were then compared with section <italic>in-situs</italic> (##FIG##0##Fig. 1##C) performed as in ##REF##11118473##Moorman et al. (2001)##. This comparison shows that the virtual section captures the extent and range of the <italic>Wnt5a</italic> expression pattern as accurately as the section <italic>in-situ</italic> with the exception of some apical ectodermal ridge (AER) expression (##FIG##0##Fig. 1##C arrowed). For a further illustration of the effective capture of expression patterns using OPT see <xref rid="app1" ref-type="sec">supplementary data</xref> (Fig. S3–5).</p> <p>We then measured the mean grey level signal intensity in all individual sections along the proximo-distal axis. The plot of these data (##FIG##0##Fig. 1##D) shows low levels of signal in the proximal region (Red, slice 1–28), either very low expression or background. The mean grey level intensity then climbs steeply in the medial region (Orange, slice 29–52) from ∼ 30 up to 100. In the distal region (Green, 53–75) the mean grey level intensity increases less steeply and then levels off at a mean intensity of around 160. In the final 5 slices the mean grey level intensity drops rapidly. This shows the capability of the OPT imaging method to allow a detailed analysis of graded patterns of expression.</p> <p>We also compared <italic>Wnt5a</italic> expression levels as measured from OPT scans of whole mount <italic>in-situs</italic> with real-time RT-PCR analysis (##FIG##0##Fig. 1##F). For both OPT and RT-PCR analyses the limb bud was divided into three regions of equal length along the proximo-distal axis designated proximal, medial and distal (##FIG##0##Fig. 1##E). For RT-PCR the sample tissue was dissected into the three regions of equal length and samples from 10 embryos were pooled. In the case of the OPT data this division was performed digitally using AMIRA’s segmentation software based upon guidance from the researcher who performed the initial microdissection. The real time RT-PCR data produced relative values for the expression of <italic>Wnt5a</italic> as follows; expression in the proximal region was taken as the reference expression level, the medial region had a 5 fold increase over the proximal region and the distal region a 19.7 fold increase. The OPT based analysis produced relative values for medial and distal regions of 4.5 fold and 6 fold increases over the value for the proximal region respectively. Therefore, WISH/OPT captures the graded nature of the expression along the proximo-distal axis, indeed the correspondence in the proximal and medial regions is striking, but not across the whole quantitative range captured using RT-PCR as the correspondence falls off dramatically at the higher levels seen in the distal region. The limitations in the captured range of expression may be due to limitations of the WISH detection method, <italic>i.e.</italic>, a nonlinear relationship between signal intensity and RNA quantity. This allows comparisons of the level of expression of a particular gene within a particular scan but not the comparison of absolute levels between different scans, although high and low regions of expression could be compared.</p> <p>The quality of data capture for the gradient of <italic>Wnt5A</italic>, both from selected domains and from serial virtual sections, suggests that WISH/OPT is suitable for examining complex patterns of graded expression in tissues within developing embryos, but not for quantification of signal with high accuracy and quantitative comparisons of mRNA levels between samples.</p> <title>Reference models for comparative analysis</title> <p>An initial requirement for meaningful comparison of gene expression patterns is a common spatial reference framework onto which different patterns can be mapped. We have produced a panel of such reference frameworks for several embryonic stages; whole embryos, isolated and fixed according to a standardised protocol, were collected at Hamburger Hamilton (HH) (##REF##24539719##Hamburger and Hamilton, 1951##) stages from HH18 through to HH25. Embryos were then scanned by OPT using autofluorescence stimulated by a UV lamp and reconstructed to produce 3D models. These models were rendered using the AMIRA software package to show morphology and gross anatomy of the embryos (##FIG##1##Figs. 2##A–H). This resulted in a clear 3D visualisation of the embryo and revealed details such as the AER – the thickened layer of epithelium that rims the distal limb buds – in the models from HH20 onwards (##FIG##1##Figs. 2##A–F, blue arrows indicate AER). This is best appreciated when the model is rotated (##FIG##1##Fig. 2##H″, blue arrows indicate AER) (to view models in 3D see movies in <xref rid="app1" ref-type="sec">supplementary data Fig. S6</xref>). Measurements of the length (<italic>L</italic>; along a line from anterior join of bud and trunk to posterior join) and width (<italic>W</italic>; a line from distal tip to trunk perpendicular to line <italic>L</italic>) of wing buds were made on the dorsal plan view (##FIG##1##Fig. 2##, Table 1) and the <italic>L</italic>/<italic>W</italic> ratio calculated in order to stage the reference embryo (<xref rid="app1" ref-type="sec">Fig. S2</xref>). The procedure was the same as one would use to stage a living chick embryo using the staging criteria of Hamburger and Hamilton (##FIG##1##Fig. 2##, Table 2) (<xref rid="app1" ref-type="sec">supplementary Fig. S2</xref>).</p> <p>Reference models for specific regions, such as the developing limb buds (##FIG##1##Figs. 2##A′–H′), can be extracted from whole embryo models and used for mapping of gene expression in these regions. Within the AMIRA program, such extracted models remain in register with the model of the whole embryo from which they are derived, thus allowing maintenance of a consistent positional system between all expression patterns mapped. These relative positions are also maintained when files are exported to the Wlz file format, which can be used to store 3 dimensional greyscale image data, using software tools developed by the MRC HGU. Subsequent mappings and analyses of gene expression reported here were performed on an extracted data set for the right wing bud of the late stage HH22 reference embryo (##FIG##1##Fig. 2##E′). Further reference sets for other stages can be easily generated.</p> <title>Reproducibility of 3D mapping of gene expression within and between labs</title> <p>To test the reliability of our <italic>in-situ</italic> protocol and 3D mapping, we focussed on the expression domain of <italic>Sonic hedgehog</italic> (<italic>Shh</italic>) in HH22 wing buds. <italic>Shh</italic> is expressed in the polarizing region at the posterior margin of the wing bud and <italic>Shh</italic> expression correlates with maps of polarizing activity (##REF##8269518##Riddle et al., 1993##). We compared data generated from wing buds in a single round of whole mount <italic>in-situ</italic> hybridisation experiments and from wing buds processed in three different labs.</p> <p>A round of <italic>in-situ</italic> hybridisations using standardised protocols for <italic>Shh</italic> expression were carried out on 4 embryos and the sense control probe was used on 2 embryos. All 4 embryos assayed using the Shh anti-sense probe were treated identically and detection was carried out in the same tube. Nevertheless there were differences in intensity of staining of <italic>Shh</italic> transcripts in the polarizing region (##FIG##2##Figs. 3##A–D) with wing buds of one embryo (##FIG##2##Fig. 3##C) showing very faint staining. Wing buds of control embryos (data not shown) had no visible background staining. <italic>In-situs</italic> of embryos from Edinburgh (##FIG##2##Fig. 3##E) and Dublin (##FIG##2##Fig. 3##F) showed similar localisation of <italic>Shh</italic> transcripts in the wing bud and one embryo from each site was then scanned together with the four embryos from the run carried out in Dundee.</p> <p>All six embryos were OPT scanned through different channels to capture a) autofluorescence to represent anatomy and b) the staining pattern under visible light. Having reconstructed 3D representations of each, we then digitally extracted the right wing bud and accompanying flank using the same spatial parameters. Co-visualisation of both anatomy (grey) and expression (orange) with volume rendering (##FIG##2##Figs. 3##A′–F′) shows the representation of the original <italic>in-situ</italic> data (##FIG##2##Figs. 3##A–F) following OPT scanning and reconstruction.</p> <p>The patterns of <italic>Shh</italic> expression in each of the six wing buds were then mapped in 3D to the HH22 reference wing (##FIG##1##Fig. 2##E). ##FIG##2##Figs. 3##A″–F″ shows heat maps of intensity of <italic>Shh</italic> expression in one section, taken across the antero-posterior/dorso-ventral (A-Po/Do-V) axes of the HH22 stage reference wing bud in a plane situated next to the AER (##FIG##2##Fig. 3##L), for the individual patterns of expression for the six source wing buds. Signal intensity of mapped gene expression data is represented by the heatmap in ##FIG##2##Fig. 3##A″ corresponding to grey scale values between 1 and 255, this measure of intensity is not suitable for precise quantitative comparisons between samples but allows visualisation of the differing levels of expression within a sample within the limits of the WISH methodology. Of the 4 wing buds from the Dundee laboratory one showed a weak signal (##FIG##2##Fig. 3##C″), with a maximal intensity of only 28. The sense controls were both very clean with no signal (data not shown). Scans for Edinburgh (##FIG##2##Fig. 3##E″) and Dublin (##FIG##2##Fig. 3##F″) had a localisation very similar to ##FIG##2##Fig. 3##A″ but an intensity of expression much closer to ##FIG##2##Fig. 3##D″.</p> <p>Both unique domains of expression and intersecting domains of expression can be derived for these data sets. Visualisation of the unique expression domains in 3 dimensions shows no unique domains for the specimens shown in A–D and only small peripheral regions for the higher intensity signal data sets from specimens in E and F (Data not shown). The intersect of the expression domains between all the scans (##FIG##2##Fig. 3##G) is restricted by the small domain of the outlying data set (##FIG##2##Fig. 3##C″), if we remove this outlying data from the calculation we have an intersect domain almost twice as large (##FIG##2##Fig. 3##H).</p> <p>Clearly, there is some variation between individual scans and an occasional extreme outlier, but our data suggest a clear common domain of gene expression is identifiable. To produce a reliable and robust domain which could compensate for the variations we see in individual <italic>in-situs</italic>, we incorporated the data from our multiple scans into one domain. We produced means of the data from the scans and corrected for background from the controls, using both raw and normalised data (##FIG##2##Figs. 3##I–J). The mean of the datasets (##FIG##2##Fig. 3##I) was heavily influenced by high intensity samples (compare ##FIG##2##Fig. 3##E″ and ##FIG##2##Fig. 3##I). To correct for the variation in signal intensity, whether as a result of differences in the <italic>in-situ</italic> itself or from the scanning steps, data sets were normalised by stretching their entire grey value range to cover the maximal range of 0–255. Such correction had almost no effect (##FIG##2##Fig. 3##J), although there was a small extension of the domain proximally. As an alternative to normalisation to account for variability and extreme outlying values we derived a median value for each voxel based on the grey levels of all of the scans (##FIG##2##Fig. 3##K). This median value seemed less dominated by outliers and a better representation of the range seen across the differing samples although it produced a more conservative domain than the simple mean since it removed areas where signal was not apparent in more than half of the samples.</p> <p>These results suggest that best practice for producing a reliable domain of expression for comparison is to perform several <italic>in-situs</italic> developed to suitable stain intensity and merge the resultant data. A minimum of four scans seems advisable to contribute to the merged data set and these may then be mapped to a common reference and a mean or median expression pattern derived. A mean of the patterns appears to better emphasise the extent of the expression domain while the median produces a more conservative domain less affected by outliers. As our analysis shows, more complicated treatment of the data seems to make little difference to the resulting domains although, for <italic>in-situs</italic> with persistently low signal, normalisation might help emphasise gradients of expression in some cases.</p> <p>This averaging or otherwise merging together of multiple samples is less necessary in the case of previously well characterised genes where the expected pattern of expression is already known and a representative sample can be confidently identified. This approach should be most valuable in the situation where the gene expression is either poorly characterised or unknown, as is likely to be the case in large scale screens.</p> <title>Comparative analysis of 3D Gene expression patterns</title> <p>OPT is a rapid method of capturing the 3D expression pattern and can allow data on multiple genes to be integrated into a common framework, therefore we used 3D warping of OPT data to our reference models to produce such an integrated data set. Expression patterns for the genes <italic>Shh</italic>, <italic>HoxD13</italic>, <italic>Fgf8</italic>, <italic>Msx1</italic>, <italic>Lmx1</italic>, <italic>Wnt5a</italic> and <italic>Tbx3</italic> (for representative <italic>in-situs</italic> see <xref rid="app1" ref-type="sec">supplementary Fig. S6</xref>) were mapped to our HH22 reference model (##FIG##1##Fig. 2##E) using AMIRA’s 3D warping capabilities. This stage was chosen as it represents a well-developed limb bud but still consisting mainly of undifferentiated mesenchyme cells. These genes were chosen for particular characteristics of expression such as dorsal restriction, <italic>Lmx1</italic>, specific expression in the AER, <italic>Fgf8</italic>, specific expression in the mesoderm, <italic>Shh</italic>, or particular gradients of expression, such as the proximo-distal gradient of <italic>Wnt5a</italic>. Particular specimens for scanning were chosen based on <italic>in-situ</italic> quality in comparison to others processed with them, usually the best example from 5–6 <italic>in-situs</italic>. The resulting mappings were then visualised in 3D (see <xref rid="app1" ref-type="sec">supplementary materials Fig. S8</xref>), virtual sections were derived along specific planes (##FIG##3##Fig. 4##A), gene expression patterns and intensity were visualised on the anterior–posterior/dorsal–ventral plane (##FIG##3##Fig. 4##A.i). The expression domains of several pairs of genes were co-visualised on the antero-posterior/dorso-ventral (A-Po/Do-V) (##FIG##3##Fig. 4##C, section plane in ##FIG##3##Fig. 4##A.i), antero-posterior/proximo-distal (A-Po/Pr-Di) (##FIG##3##Fig. 4##D, section plane in ##FIG##3##Fig. 4##A.ii) and dorso-ventral/proximo-distal (Do-V/Pr-Di) (##FIG##3##Fig. 4##E, section plane in ##FIG##3##Fig. 4##A.iii) planes to allow some specific comparisons to be drawn (##FIG##3##Figs. 4##C–E; for 3D visualisations see <xref rid="app1" ref-type="sec">supplementary materials Fig. S11–18</xref>). These mappings show both expected features, such as the dorsal expression pattern of <italic>Lmx1</italic>, and novel features such as an apparent gradient of expression from ventral to dorsal in <italic>Wnt5a</italic>. Indeed several dorso-ventral asymmetries were by far the most striking features to emerge.</p> <p>One such previously unappreciated asymmetry was in <italic>Shh</italic> expression (##FIG##3##Fig. 4##B), which shows that <italic>Shh</italic> expression extends further anteriorly on the ventral side of the limb. A transverse section through the limb along plane E shows that the asymmetry is more complex and the domain of <italic>Shh</italic> expression is skewed with the more proximal regions of the limb showing a more dorsal expression of <italic>Shh</italic> (##FIG##3##Fig. 4##E.ii). Planes taken at more anterior levels lose this obvious skewing (see <xref rid="app1" ref-type="sec">supplementary Fig. S9</xref>) and elements of this distribution are confirmed by both normal <italic>Shh</italic> whole mount <italic>in-situ</italic> and a double <italic>in-situ</italic> for <italic>Shh</italic> and <italic>Fgf8</italic> (see <xref rid="app1" ref-type="sec">supplementary materials Fig. S10</xref>). This may explain the observations of ##REF##7697724##Yang and Niswander (1995)## who reported no apparent dorso-ventral asymmetry in <italic>Shh</italic> expression at HH24 based on sectioned whole mount <italic>in-situ</italic> hybridisation data. Alternatively this may be due to the different stages assayed. A more dorsal distribution in the more proximal regions of the limb may suggest that a dorsally localised signal, such as <italic>Wnt7a</italic>, might be maintaining dorsal expression proximally while more distal expression would be maintained evenly across the dorso-ventral axis by signals from the apical ectodermal ridge. Indeed it is already known that Wnt7a plays a role in maintaining <italic>Shh</italic> expression in the polarising region (<xref rid="bib38 bib44" ref-type="bibr">Parr and McMahon, 1995; Yang and Niswander, 1995</xref>). Furthermore ##REF##10445500##Kawakami et al. (1999)## reported that <italic>Frizzled10</italic> (a Wnt receptor) colocalises with <italic>Shh</italic> dorsally and suggested that <italic>Shh</italic> expression in this part of the polarising region might be regulated by Wnt7a (##REF##11142678##Kawakami et al., 2000##).</p> <p>The pairwise comparisons similarly produced both expected and unexpected results. <italic>HoxD13</italic> shows striking asymmetrical dorso-ventral expression (##FIG##3##Figs. 4##C.v and E.v) and the ventral margin of <italic>HoxD13</italic> expression appears to coincide with that of <italic>Lmx1</italic> (##FIG##3##Figs. 4##C.viii and E.viii; Movie in <xref rid="app1" ref-type="sec">supplementary Fig. S18</xref>), suggesting a possible regulatory relationship. The dorso-ventral asymmetry in <italic>Hoxd13</italic> had been previously noted by Duboule et al. from 3D reconstructions based on radioactive section <italic>in-situs</italic> (##UREF##2##Olivo et al., 1993##) but at this time <italic>Lmx1</italic> was not known. Recent lineage tracing studies in the developing mouse limb have shown the existence of a dorsoventral lineage restriction compartment further suggesting that there is still considerable complexity in the dorsoventral organisation of the limb to be discovered (##REF##17715176##Arques et al., 2007##).</p> <p>Comparison of <italic>Shh</italic> and <italic>Fgf8</italic> expression (##FIG##3##Figs. 4##D.iii and E.iii; movie in <xref rid="app1" ref-type="sec">supplementary Fig. S13</xref>) shows an unexpected overlap of expression in the mesoderm. This shows that there are important limits of spatial resolution to the current mapping procedure given the well-characterised expression of these genes in mesoderm and apical ectodermal ridge respectively. Since the AER is one of the principle morphological features of the limb, it is heavily used in the landmarking process, which is the first stage of mapping expression data to the reference limb. Strong <italic>in-situ</italic> staining, as seen with <italic>Fgf8</italic>, can block autofluorescence and prevent the capture of the anatomical data for the ridge. A comparison of <italic>Wnt5a</italic> and <italic>Fgf8</italic> (Movie in <xref rid="app1" ref-type="sec">supplementary Fig. S11</xref>) also shows considerable overlap except for the most anterior region of <italic>Fgf8</italic> expression (##FIG##3##Fig. 4##C.i) and a persistent ‘leading edge’ of <italic>Fgf8</italic> expression in the most distal portion of the limb (##FIG##3##Fig. 4##E.i). In this case, we would expect to see considerable overlap due to expression of <italic>Wnt5a</italic> in the ridge, although it is not clear whether the common domain of expression accurately represents just ridge expression.</p> <title>OPT using Tyramide Signal Amplification</title> <p>To address the problem of the loss of anatomical landmarks when looking at strongly expressed ectodermal genes, such as <italic>Fgf8</italic>, we used a Tyramide Signal Amplification (TSA) kit (Invitrogen) to enhance a fluorescent colour reaction and avoid the blocking effect seen with chromagenic substrates such as NBT–BCIP. We compared the expression of <italic>Tbx3</italic> and <italic>Fgf8</italic> at HH24, mapped to the HH24 reference limb (##FIG##1##Fig. 2##G′), from a normal NBT–BCIP based colour reaction for <italic>Tbx3</italic> and a TSA fluorescent colour reaction for <italic>Fgf8</italic> (##FIG##4##Figs. 5##A, B). This approach provided a much more accurate localisation of <italic>Fgf8</italic> expression to the apical ridge compared to that seen in ##FIG##3##Fig. 4##. When we carried out a pairwise comparison between <italic>Fgf8</italic> and <italic>Tbx3</italic> expression patterns, the overlap was very much reduced compared to that seen between <italic>Fgf8</italic> and <italic>Shh</italic> in ##FIG##3##Fig. 4##D.iii and compared with <italic>Fgf8</italic> and <italic>Tbx3</italic> from our original data set in which we had used an NBT–BCIP reaction for both genes (##FIG##4##Fig. 5##C). While these mapped localisations are not sufficient for accurately discriminating expression from expressing tissues close together or very thin tissue layers a database containing such mapped data would be linked to the original 3D scans allowing visualisation of the data.</p> <title>Computational comparisons of mapped gene expression patterns</title> <p>Despite the shortcomings discussed above, our overall mapping strategy provides a valuable tool for analysing spatial and temporal relationships of multiple complex patterns.</p> <p>To begin to apply computational methods to analyse these multiple data sets at once and to look for similarities in patterns of gene expression, we utilised software produced by the MRC HGU to manipulate 3D image data in the Wlz format, a format used by EMAGE. AMIRA files were converted to Wlz format and a set of sub-regions of 10 × 10 × 10 voxels were defined for the early HH22 reference limb producing a coarse sampling of the 3D model. The expression data for the previously analysed genes and an additional gene, <italic>Wnt3a</italic>, were then used to derive mean expression values from the grey levels of each OPT scan for the newly defined volumes. This produced a matrix of common positional IDs and expression levels for each gene.</p> <p>These data were then analysed using hierarchical clustering for both the genes and the regions defined by the coarse sampling. This analysis was performed and visualised using the TMEV package (<ext-link ext-link-type="uri" xlink:href="http://www.tm4.org/mev.html">http://www.tm4.org/mev.html</ext-link>) (##FIG##5##Fig. 6##A). The clustering of these genes fits our expectations with <italic>Fgf8</italic> and <italic>Wnt3a</italic>, both known to be expressed in the apical ridge (<xref rid="bib4 bib25 bib35" ref-type="bibr">Barrow et al., 2003; Kengaku et al., 1998; McQueeney et al., 2002</xref>), treeing out together. Visualisation of the clustered regions shows that they largely form contiguous spatial domains (##FIG##5##Figs. 6##B–G; <xref rid="app1" ref-type="sec">supplementary Fig. S19</xref>); some of these domains are associated with specific regions of the limb. Cluster B lies in a plane through the limb along the proximal–distal axis at the level of the anterior margin of <italic>Shh</italic> expression (##FIG##5##Fig. 6##B), cluster D is associated with the AER in the anterior of the limb (##FIG##5##Fig. 6##D) while cluster E is in the dorsal margin of the limb (##FIG##5##Fig. 6##E). None of the clusters visualised here corresponds simply to the expression pattern of one particular gene.</p> <p>This simple level of clustering computational analysis shows the potential for methods developed to study gene expression data from other sources, such as microarray data, to be applied to the study of 3D gene expression data. Not only can we identify similarly expressed genes using this method but it should also be possible to identify specific regions of the limb that may be important in the regulation of gene expression, such as novel signalling centres.</p>
<title>Results and discussion</title> <title>Assessment of efficiency of whole mount <italic>in-situ</italic> hybridisation and scanning with OPT as a method of detecting gene expression domains</title> <p>An initial technical issue we encountered with OPT scanning of WISH (whole mount <italic>in-situ</italic> hybridisation) specimens was that strong <italic>in-situ</italic> colour reaction staining can block autofluorescence and prevent the capture of portions of the anatomical data required for subsequent mapping to a reference model. To obviate this problem, we identified a particular depth of staining with the NBT–BCIP substrate, suitable for OPT scanning, which captures an extensive range of the expression pattern and allows the visualisation of dynamic gradients without causing a substantial dropout of the anatomical data necessary for mapping (see <xref rid="app1" ref-type="sec">supplementary materials Fig. S</xref>1). To test our standard <italic>in-situ</italic> hybridisation protocol, monitor probe penetration of the embryonic limb, and assess the ability of the OPT system to identify graded patterns of expression within deep tissues, we focussed on <italic>Wnt5a</italic>. <italic>Wnt5a</italic> has been reported to have a proximo-distal gradient of expression based on both radioactive section <italic>in-situ</italic> hybridisation and northern blots of distinct portions of the limb (##REF##8297789##Dealy et al., 1993##) and is known to be expressed throughout mesodermal regions of the limb.</p> <p>We first assayed expression by WISH (method modified after ##REF##8722478##Nieto et al., 1996## see <xref rid="app1" ref-type="sec">supplementary data</xref>) ##FIG##0##Fig. 1##A). The <italic>Wnt5a</italic> whole mount was scanned using OPT and gene expression data mapped onto a reference limb (##FIG##1##Fig. 2##D), from which virtual sections were derived (##FIG##0##Fig. 1##B). These virtual sections were then compared with section <italic>in-situs</italic> (##FIG##0##Fig. 1##C) performed as in ##REF##11118473##Moorman et al. (2001)##. This comparison shows that the virtual section captures the extent and range of the <italic>Wnt5a</italic> expression pattern as accurately as the section <italic>in-situ</italic> with the exception of some apical ectodermal ridge (AER) expression (##FIG##0##Fig. 1##C arrowed). For a further illustration of the effective capture of expression patterns using OPT see <xref rid="app1" ref-type="sec">supplementary data</xref> (Fig. S3–5).</p> <p>We then measured the mean grey level signal intensity in all individual sections along the proximo-distal axis. The plot of these data (##FIG##0##Fig. 1##D) shows low levels of signal in the proximal region (Red, slice 1–28), either very low expression or background. The mean grey level intensity then climbs steeply in the medial region (Orange, slice 29–52) from ∼ 30 up to 100. In the distal region (Green, 53–75) the mean grey level intensity increases less steeply and then levels off at a mean intensity of around 160. In the final 5 slices the mean grey level intensity drops rapidly. This shows the capability of the OPT imaging method to allow a detailed analysis of graded patterns of expression.</p> <p>We also compared <italic>Wnt5a</italic> expression levels as measured from OPT scans of whole mount <italic>in-situs</italic> with real-time RT-PCR analysis (##FIG##0##Fig. 1##F). For both OPT and RT-PCR analyses the limb bud was divided into three regions of equal length along the proximo-distal axis designated proximal, medial and distal (##FIG##0##Fig. 1##E). For RT-PCR the sample tissue was dissected into the three regions of equal length and samples from 10 embryos were pooled. In the case of the OPT data this division was performed digitally using AMIRA’s segmentation software based upon guidance from the researcher who performed the initial microdissection. The real time RT-PCR data produced relative values for the expression of <italic>Wnt5a</italic> as follows; expression in the proximal region was taken as the reference expression level, the medial region had a 5 fold increase over the proximal region and the distal region a 19.7 fold increase. The OPT based analysis produced relative values for medial and distal regions of 4.5 fold and 6 fold increases over the value for the proximal region respectively. Therefore, WISH/OPT captures the graded nature of the expression along the proximo-distal axis, indeed the correspondence in the proximal and medial regions is striking, but not across the whole quantitative range captured using RT-PCR as the correspondence falls off dramatically at the higher levels seen in the distal region. The limitations in the captured range of expression may be due to limitations of the WISH detection method, <italic>i.e.</italic>, a nonlinear relationship between signal intensity and RNA quantity. This allows comparisons of the level of expression of a particular gene within a particular scan but not the comparison of absolute levels between different scans, although high and low regions of expression could be compared.</p> <p>The quality of data capture for the gradient of <italic>Wnt5A</italic>, both from selected domains and from serial virtual sections, suggests that WISH/OPT is suitable for examining complex patterns of graded expression in tissues within developing embryos, but not for quantification of signal with high accuracy and quantitative comparisons of mRNA levels between samples.</p> <title>Reference models for comparative analysis</title> <p>An initial requirement for meaningful comparison of gene expression patterns is a common spatial reference framework onto which different patterns can be mapped. We have produced a panel of such reference frameworks for several embryonic stages; whole embryos, isolated and fixed according to a standardised protocol, were collected at Hamburger Hamilton (HH) (##REF##24539719##Hamburger and Hamilton, 1951##) stages from HH18 through to HH25. Embryos were then scanned by OPT using autofluorescence stimulated by a UV lamp and reconstructed to produce 3D models. These models were rendered using the AMIRA software package to show morphology and gross anatomy of the embryos (##FIG##1##Figs. 2##A–H). This resulted in a clear 3D visualisation of the embryo and revealed details such as the AER – the thickened layer of epithelium that rims the distal limb buds – in the models from HH20 onwards (##FIG##1##Figs. 2##A–F, blue arrows indicate AER). This is best appreciated when the model is rotated (##FIG##1##Fig. 2##H″, blue arrows indicate AER) (to view models in 3D see movies in <xref rid="app1" ref-type="sec">supplementary data Fig. S6</xref>). Measurements of the length (<italic>L</italic>; along a line from anterior join of bud and trunk to posterior join) and width (<italic>W</italic>; a line from distal tip to trunk perpendicular to line <italic>L</italic>) of wing buds were made on the dorsal plan view (##FIG##1##Fig. 2##, Table 1) and the <italic>L</italic>/<italic>W</italic> ratio calculated in order to stage the reference embryo (<xref rid="app1" ref-type="sec">Fig. S2</xref>). The procedure was the same as one would use to stage a living chick embryo using the staging criteria of Hamburger and Hamilton (##FIG##1##Fig. 2##, Table 2) (<xref rid="app1" ref-type="sec">supplementary Fig. S2</xref>).</p> <p>Reference models for specific regions, such as the developing limb buds (##FIG##1##Figs. 2##A′–H′), can be extracted from whole embryo models and used for mapping of gene expression in these regions. Within the AMIRA program, such extracted models remain in register with the model of the whole embryo from which they are derived, thus allowing maintenance of a consistent positional system between all expression patterns mapped. These relative positions are also maintained when files are exported to the Wlz file format, which can be used to store 3 dimensional greyscale image data, using software tools developed by the MRC HGU. Subsequent mappings and analyses of gene expression reported here were performed on an extracted data set for the right wing bud of the late stage HH22 reference embryo (##FIG##1##Fig. 2##E′). Further reference sets for other stages can be easily generated.</p> <title>Reproducibility of 3D mapping of gene expression within and between labs</title> <p>To test the reliability of our <italic>in-situ</italic> protocol and 3D mapping, we focussed on the expression domain of <italic>Sonic hedgehog</italic> (<italic>Shh</italic>) in HH22 wing buds. <italic>Shh</italic> is expressed in the polarizing region at the posterior margin of the wing bud and <italic>Shh</italic> expression correlates with maps of polarizing activity (##REF##8269518##Riddle et al., 1993##). We compared data generated from wing buds in a single round of whole mount <italic>in-situ</italic> hybridisation experiments and from wing buds processed in three different labs.</p> <p>A round of <italic>in-situ</italic> hybridisations using standardised protocols for <italic>Shh</italic> expression were carried out on 4 embryos and the sense control probe was used on 2 embryos. All 4 embryos assayed using the Shh anti-sense probe were treated identically and detection was carried out in the same tube. Nevertheless there were differences in intensity of staining of <italic>Shh</italic> transcripts in the polarizing region (##FIG##2##Figs. 3##A–D) with wing buds of one embryo (##FIG##2##Fig. 3##C) showing very faint staining. Wing buds of control embryos (data not shown) had no visible background staining. <italic>In-situs</italic> of embryos from Edinburgh (##FIG##2##Fig. 3##E) and Dublin (##FIG##2##Fig. 3##F) showed similar localisation of <italic>Shh</italic> transcripts in the wing bud and one embryo from each site was then scanned together with the four embryos from the run carried out in Dundee.</p> <p>All six embryos were OPT scanned through different channels to capture a) autofluorescence to represent anatomy and b) the staining pattern under visible light. Having reconstructed 3D representations of each, we then digitally extracted the right wing bud and accompanying flank using the same spatial parameters. Co-visualisation of both anatomy (grey) and expression (orange) with volume rendering (##FIG##2##Figs. 3##A′–F′) shows the representation of the original <italic>in-situ</italic> data (##FIG##2##Figs. 3##A–F) following OPT scanning and reconstruction.</p> <p>The patterns of <italic>Shh</italic> expression in each of the six wing buds were then mapped in 3D to the HH22 reference wing (##FIG##1##Fig. 2##E). ##FIG##2##Figs. 3##A″–F″ shows heat maps of intensity of <italic>Shh</italic> expression in one section, taken across the antero-posterior/dorso-ventral (A-Po/Do-V) axes of the HH22 stage reference wing bud in a plane situated next to the AER (##FIG##2##Fig. 3##L), for the individual patterns of expression for the six source wing buds. Signal intensity of mapped gene expression data is represented by the heatmap in ##FIG##2##Fig. 3##A″ corresponding to grey scale values between 1 and 255, this measure of intensity is not suitable for precise quantitative comparisons between samples but allows visualisation of the differing levels of expression within a sample within the limits of the WISH methodology. Of the 4 wing buds from the Dundee laboratory one showed a weak signal (##FIG##2##Fig. 3##C″), with a maximal intensity of only 28. The sense controls were both very clean with no signal (data not shown). Scans for Edinburgh (##FIG##2##Fig. 3##E″) and Dublin (##FIG##2##Fig. 3##F″) had a localisation very similar to ##FIG##2##Fig. 3##A″ but an intensity of expression much closer to ##FIG##2##Fig. 3##D″.</p> <p>Both unique domains of expression and intersecting domains of expression can be derived for these data sets. Visualisation of the unique expression domains in 3 dimensions shows no unique domains for the specimens shown in A–D and only small peripheral regions for the higher intensity signal data sets from specimens in E and F (Data not shown). The intersect of the expression domains between all the scans (##FIG##2##Fig. 3##G) is restricted by the small domain of the outlying data set (##FIG##2##Fig. 3##C″), if we remove this outlying data from the calculation we have an intersect domain almost twice as large (##FIG##2##Fig. 3##H).</p> <p>Clearly, there is some variation between individual scans and an occasional extreme outlier, but our data suggest a clear common domain of gene expression is identifiable. To produce a reliable and robust domain which could compensate for the variations we see in individual <italic>in-situs</italic>, we incorporated the data from our multiple scans into one domain. We produced means of the data from the scans and corrected for background from the controls, using both raw and normalised data (##FIG##2##Figs. 3##I–J). The mean of the datasets (##FIG##2##Fig. 3##I) was heavily influenced by high intensity samples (compare ##FIG##2##Fig. 3##E″ and ##FIG##2##Fig. 3##I). To correct for the variation in signal intensity, whether as a result of differences in the <italic>in-situ</italic> itself or from the scanning steps, data sets were normalised by stretching their entire grey value range to cover the maximal range of 0–255. Such correction had almost no effect (##FIG##2##Fig. 3##J), although there was a small extension of the domain proximally. As an alternative to normalisation to account for variability and extreme outlying values we derived a median value for each voxel based on the grey levels of all of the scans (##FIG##2##Fig. 3##K). This median value seemed less dominated by outliers and a better representation of the range seen across the differing samples although it produced a more conservative domain than the simple mean since it removed areas where signal was not apparent in more than half of the samples.</p> <p>These results suggest that best practice for producing a reliable domain of expression for comparison is to perform several <italic>in-situs</italic> developed to suitable stain intensity and merge the resultant data. A minimum of four scans seems advisable to contribute to the merged data set and these may then be mapped to a common reference and a mean or median expression pattern derived. A mean of the patterns appears to better emphasise the extent of the expression domain while the median produces a more conservative domain less affected by outliers. As our analysis shows, more complicated treatment of the data seems to make little difference to the resulting domains although, for <italic>in-situs</italic> with persistently low signal, normalisation might help emphasise gradients of expression in some cases.</p> <p>This averaging or otherwise merging together of multiple samples is less necessary in the case of previously well characterised genes where the expected pattern of expression is already known and a representative sample can be confidently identified. This approach should be most valuable in the situation where the gene expression is either poorly characterised or unknown, as is likely to be the case in large scale screens.</p> <title>Comparative analysis of 3D Gene expression patterns</title> <p>OPT is a rapid method of capturing the 3D expression pattern and can allow data on multiple genes to be integrated into a common framework, therefore we used 3D warping of OPT data to our reference models to produce such an integrated data set. Expression patterns for the genes <italic>Shh</italic>, <italic>HoxD13</italic>, <italic>Fgf8</italic>, <italic>Msx1</italic>, <italic>Lmx1</italic>, <italic>Wnt5a</italic> and <italic>Tbx3</italic> (for representative <italic>in-situs</italic> see <xref rid="app1" ref-type="sec">supplementary Fig. S6</xref>) were mapped to our HH22 reference model (##FIG##1##Fig. 2##E) using AMIRA’s 3D warping capabilities. This stage was chosen as it represents a well-developed limb bud but still consisting mainly of undifferentiated mesenchyme cells. These genes were chosen for particular characteristics of expression such as dorsal restriction, <italic>Lmx1</italic>, specific expression in the AER, <italic>Fgf8</italic>, specific expression in the mesoderm, <italic>Shh</italic>, or particular gradients of expression, such as the proximo-distal gradient of <italic>Wnt5a</italic>. Particular specimens for scanning were chosen based on <italic>in-situ</italic> quality in comparison to others processed with them, usually the best example from 5–6 <italic>in-situs</italic>. The resulting mappings were then visualised in 3D (see <xref rid="app1" ref-type="sec">supplementary materials Fig. S8</xref>), virtual sections were derived along specific planes (##FIG##3##Fig. 4##A), gene expression patterns and intensity were visualised on the anterior–posterior/dorsal–ventral plane (##FIG##3##Fig. 4##A.i). The expression domains of several pairs of genes were co-visualised on the antero-posterior/dorso-ventral (A-Po/Do-V) (##FIG##3##Fig. 4##C, section plane in ##FIG##3##Fig. 4##A.i), antero-posterior/proximo-distal (A-Po/Pr-Di) (##FIG##3##Fig. 4##D, section plane in ##FIG##3##Fig. 4##A.ii) and dorso-ventral/proximo-distal (Do-V/Pr-Di) (##FIG##3##Fig. 4##E, section plane in ##FIG##3##Fig. 4##A.iii) planes to allow some specific comparisons to be drawn (##FIG##3##Figs. 4##C–E; for 3D visualisations see <xref rid="app1" ref-type="sec">supplementary materials Fig. S11–18</xref>). These mappings show both expected features, such as the dorsal expression pattern of <italic>Lmx1</italic>, and novel features such as an apparent gradient of expression from ventral to dorsal in <italic>Wnt5a</italic>. Indeed several dorso-ventral asymmetries were by far the most striking features to emerge.</p> <p>One such previously unappreciated asymmetry was in <italic>Shh</italic> expression (##FIG##3##Fig. 4##B), which shows that <italic>Shh</italic> expression extends further anteriorly on the ventral side of the limb. A transverse section through the limb along plane E shows that the asymmetry is more complex and the domain of <italic>Shh</italic> expression is skewed with the more proximal regions of the limb showing a more dorsal expression of <italic>Shh</italic> (##FIG##3##Fig. 4##E.ii). Planes taken at more anterior levels lose this obvious skewing (see <xref rid="app1" ref-type="sec">supplementary Fig. S9</xref>) and elements of this distribution are confirmed by both normal <italic>Shh</italic> whole mount <italic>in-situ</italic> and a double <italic>in-situ</italic> for <italic>Shh</italic> and <italic>Fgf8</italic> (see <xref rid="app1" ref-type="sec">supplementary materials Fig. S10</xref>). This may explain the observations of ##REF##7697724##Yang and Niswander (1995)## who reported no apparent dorso-ventral asymmetry in <italic>Shh</italic> expression at HH24 based on sectioned whole mount <italic>in-situ</italic> hybridisation data. Alternatively this may be due to the different stages assayed. A more dorsal distribution in the more proximal regions of the limb may suggest that a dorsally localised signal, such as <italic>Wnt7a</italic>, might be maintaining dorsal expression proximally while more distal expression would be maintained evenly across the dorso-ventral axis by signals from the apical ectodermal ridge. Indeed it is already known that Wnt7a plays a role in maintaining <italic>Shh</italic> expression in the polarising region (<xref rid="bib38 bib44" ref-type="bibr">Parr and McMahon, 1995; Yang and Niswander, 1995</xref>). Furthermore ##REF##10445500##Kawakami et al. (1999)## reported that <italic>Frizzled10</italic> (a Wnt receptor) colocalises with <italic>Shh</italic> dorsally and suggested that <italic>Shh</italic> expression in this part of the polarising region might be regulated by Wnt7a (##REF##11142678##Kawakami et al., 2000##).</p> <p>The pairwise comparisons similarly produced both expected and unexpected results. <italic>HoxD13</italic> shows striking asymmetrical dorso-ventral expression (##FIG##3##Figs. 4##C.v and E.v) and the ventral margin of <italic>HoxD13</italic> expression appears to coincide with that of <italic>Lmx1</italic> (##FIG##3##Figs. 4##C.viii and E.viii; Movie in <xref rid="app1" ref-type="sec">supplementary Fig. S18</xref>), suggesting a possible regulatory relationship. The dorso-ventral asymmetry in <italic>Hoxd13</italic> had been previously noted by Duboule et al. from 3D reconstructions based on radioactive section <italic>in-situs</italic> (##UREF##2##Olivo et al., 1993##) but at this time <italic>Lmx1</italic> was not known. Recent lineage tracing studies in the developing mouse limb have shown the existence of a dorsoventral lineage restriction compartment further suggesting that there is still considerable complexity in the dorsoventral organisation of the limb to be discovered (##REF##17715176##Arques et al., 2007##).</p> <p>Comparison of <italic>Shh</italic> and <italic>Fgf8</italic> expression (##FIG##3##Figs. 4##D.iii and E.iii; movie in <xref rid="app1" ref-type="sec">supplementary Fig. S13</xref>) shows an unexpected overlap of expression in the mesoderm. This shows that there are important limits of spatial resolution to the current mapping procedure given the well-characterised expression of these genes in mesoderm and apical ectodermal ridge respectively. Since the AER is one of the principle morphological features of the limb, it is heavily used in the landmarking process, which is the first stage of mapping expression data to the reference limb. Strong <italic>in-situ</italic> staining, as seen with <italic>Fgf8</italic>, can block autofluorescence and prevent the capture of the anatomical data for the ridge. A comparison of <italic>Wnt5a</italic> and <italic>Fgf8</italic> (Movie in <xref rid="app1" ref-type="sec">supplementary Fig. S11</xref>) also shows considerable overlap except for the most anterior region of <italic>Fgf8</italic> expression (##FIG##3##Fig. 4##C.i) and a persistent ‘leading edge’ of <italic>Fgf8</italic> expression in the most distal portion of the limb (##FIG##3##Fig. 4##E.i). In this case, we would expect to see considerable overlap due to expression of <italic>Wnt5a</italic> in the ridge, although it is not clear whether the common domain of expression accurately represents just ridge expression.</p> <title>OPT using Tyramide Signal Amplification</title> <p>To address the problem of the loss of anatomical landmarks when looking at strongly expressed ectodermal genes, such as <italic>Fgf8</italic>, we used a Tyramide Signal Amplification (TSA) kit (Invitrogen) to enhance a fluorescent colour reaction and avoid the blocking effect seen with chromagenic substrates such as NBT–BCIP. We compared the expression of <italic>Tbx3</italic> and <italic>Fgf8</italic> at HH24, mapped to the HH24 reference limb (##FIG##1##Fig. 2##G′), from a normal NBT–BCIP based colour reaction for <italic>Tbx3</italic> and a TSA fluorescent colour reaction for <italic>Fgf8</italic> (##FIG##4##Figs. 5##A, B). This approach provided a much more accurate localisation of <italic>Fgf8</italic> expression to the apical ridge compared to that seen in ##FIG##3##Fig. 4##. When we carried out a pairwise comparison between <italic>Fgf8</italic> and <italic>Tbx3</italic> expression patterns, the overlap was very much reduced compared to that seen between <italic>Fgf8</italic> and <italic>Shh</italic> in ##FIG##3##Fig. 4##D.iii and compared with <italic>Fgf8</italic> and <italic>Tbx3</italic> from our original data set in which we had used an NBT–BCIP reaction for both genes (##FIG##4##Fig. 5##C). While these mapped localisations are not sufficient for accurately discriminating expression from expressing tissues close together or very thin tissue layers a database containing such mapped data would be linked to the original 3D scans allowing visualisation of the data.</p> <title>Computational comparisons of mapped gene expression patterns</title> <p>Despite the shortcomings discussed above, our overall mapping strategy provides a valuable tool for analysing spatial and temporal relationships of multiple complex patterns.</p> <p>To begin to apply computational methods to analyse these multiple data sets at once and to look for similarities in patterns of gene expression, we utilised software produced by the MRC HGU to manipulate 3D image data in the Wlz format, a format used by EMAGE. AMIRA files were converted to Wlz format and a set of sub-regions of 10 × 10 × 10 voxels were defined for the early HH22 reference limb producing a coarse sampling of the 3D model. The expression data for the previously analysed genes and an additional gene, <italic>Wnt3a</italic>, were then used to derive mean expression values from the grey levels of each OPT scan for the newly defined volumes. This produced a matrix of common positional IDs and expression levels for each gene.</p> <p>These data were then analysed using hierarchical clustering for both the genes and the regions defined by the coarse sampling. This analysis was performed and visualised using the TMEV package (<ext-link ext-link-type="uri" xlink:href="http://www.tm4.org/mev.html">http://www.tm4.org/mev.html</ext-link>) (##FIG##5##Fig. 6##A). The clustering of these genes fits our expectations with <italic>Fgf8</italic> and <italic>Wnt3a</italic>, both known to be expressed in the apical ridge (<xref rid="bib4 bib25 bib35" ref-type="bibr">Barrow et al., 2003; Kengaku et al., 1998; McQueeney et al., 2002</xref>), treeing out together. Visualisation of the clustered regions shows that they largely form contiguous spatial domains (##FIG##5##Figs. 6##B–G; <xref rid="app1" ref-type="sec">supplementary Fig. S19</xref>); some of these domains are associated with specific regions of the limb. Cluster B lies in a plane through the limb along the proximal–distal axis at the level of the anterior margin of <italic>Shh</italic> expression (##FIG##5##Fig. 6##B), cluster D is associated with the AER in the anterior of the limb (##FIG##5##Fig. 6##D) while cluster E is in the dorsal margin of the limb (##FIG##5##Fig. 6##E). None of the clusters visualised here corresponds simply to the expression pattern of one particular gene.</p> <p>This simple level of clustering computational analysis shows the potential for methods developed to study gene expression data from other sources, such as microarray data, to be applied to the study of 3D gene expression data. Not only can we identify similarly expressed genes using this method but it should also be possible to identify specific regions of the limb that may be important in the regulation of gene expression, such as novel signalling centres.</p>
<title>Conclusion</title> <p>We have developed improved technology for the production of 3D atlases of gene expression (##REF##15043218##Baldock et al., 2003##). Specifically, we have shown that OPT is a reliable and efficient way of visualising 3D patterns of gene expression in the chick limb and have been able to directly compare different patterns on reference models using a 3D warping technique. This technique allows visualisation both of samples too large for confocal microscopy (##REF##18377226##Welten et al., 2006##) and those too small for good resolution with microMRI (see ##REF##18045352##Li et al., 2007## for visualisation of chick wings at later stages). Furthermore, visualisation of gene expression is still rudimentary with microMRI (##REF##17234603##Liu et al., 2007a##).</p> <p>As more 3D patterns of gene expression are mapped, simple pairwise comparisons will no longer be sufficient to analyse more complex relationships between groups of genes and then the computational approaches we have used here will be greatly beneficial, indeed similar methods have been used to study the expression of around 20,000 genes in the mouse brain (<xref rid="bib29 bib32" ref-type="bibr">Lein et al., 2007; Liu et al., 2007b</xref>). Although we have only compared 7 genes known to be expressed in chick limb development we have already revealed some previously unappreciated asymmetries and relationships, particularly with respect to the dorso-ventral axis. These techniques will also be more generally applicable to different developing structures and organisms.</p>
<p>Currently at the UK MRC Human Genetics Unit, Western General Hospital, Crewe Road, Edinburgh, UK.</p> <p>Currently at the School of Biology, Bute Building, University of St. Andrews, St. Andrews, UK.</p> <p>Currently at the School of Biology and Biochemistry, Bath University, Bath, UK.</p> <p>Chick embryos are good models for vertebrate development due to their accessibility and manipulability. Recent large increases in available genomic data from both whole genome sequencing and EST projects provide opportunities for identifying many new developmentally important chicken genes. Traditional methods of documenting when and where specific genes are expressed in embryos using wholemount and section <italic>in-situ</italic> hybridisation do not readily allow appreciation of 3-dimensional (3D) patterns of expression, but this can be accomplished by the recently developed microscopy technique, Optical Projection Tomography (OPT). Here we show that OPT data on the developing chick wing from different labs can be reliably integrated into a common database, that OPT is efficient in capturing 3D gene expression domains and that such domains can be meaningfully compared. Novel protocols are used to compare 3D expression domains of 7 genes known to be involved in chick wing development. This reveals previously unappreciated relationships and demonstrates the potential, using modern genomic resources, for building a large scale 3D atlas of gene expression. Such an atlas could be extended to include other types of data, such as fate maps, and the approach is also more generally applicable to embryos, organs and tissues.</p> <title>Keywords</title>
<title>Supplementary data</title> <p></p> <p></p> <title>Acknowledgments</title> <p>We thank Calum Thomson from the Light Microscopy facility at the University of Dundee for help with sectioning. This work was supported in part by funding from the Medical Research Council and by funding from The Royal Society, Cheryll Tickle is a Royal Society Professor.</p> <p>We would also like to thank Harris Morrison for maintaining the OPT facility at HGU, and for his help and advice, Venkat Shanmugasundaram for his help with warping expression domains using AMIRA and Peter Hohenstein for assistance with real time RT-PCR.</p>
<fig id="fig1"><label>Fig. 1</label><caption><p>Comparison of WISH/OPT and other methods for detecting a gradient of expression. (A) Dorsal view of wholemount <italic>in-situ</italic> hybridisation of <italic>Wnt5a</italic> in a HH22 wing bud. (B) Virtual section of OPT data scanned from A mapped to a reference limb. (C) Section <italic>in-situ</italic> of Wnt5a from a HH22 wing bud, arrow indicates expression in the AER. (D) A plot of the mean signal intensity in virtual slices of the OPT data set taken along the proximo-distal axis of the limb with 0 representing the most proximal position and 75 the most distal. Colouring under the line represents the domain to which the slices belong, coloured as in panel E. (E) A surface rendering showing the early HH22 reference limb and the three assayed limb domains; proximal (red), medial (orange) and distal (green), arrows indicate the antero-posterior (A-Po) and proximo-distal (Pr-Di) axes of the limb. (F) Comparison of levels of expression in three domains assayed by OPT (purple) and real time RT-PCR (blue), error bars represent standard errors of ± 0.29, ± 1.4 and ± 1.6 for the RT-PCR measurements in the proximal medial and distal regions respectively. OPT values were based on the mean grey level intensity within the domain and standardised against the mean intensity value of the proximal domain to get a relative expression. Domain labels are coloured as in panel E.</p></caption></fig> <fig id="fig2"><label>Fig. 2</label><caption><p>3D reference models of whole embryos over several stages. (A–H) Volume renderings of OPT scans of whole embryos at reference stages HH19 (A), HH20 (B), HH21 (C), Early stage HH22 (D), Late Stage HH22 (E), HH23 (F), HH24 (G), HH25 (H), scale bars (orange) for panels A–H indicate 1 mm. Blue arrows indicate where the AER is visible on rendered embryos. (A′–H′) Plan views of digitally extracted right wing buds from corresponding whole embryo scans. Panel A′ shows arrows indicating the antero-posterior (A-Po) and proximo-distal (Pr-Di) axes of the limb. (H″) Distal view showing AER (blue arrow), white arrows indicate the antero-posterior (A-Po) and dorso-ventral (Do-V) axes of the limb. Scale bars (orange) indicate 300 μm in panels A′–H′ and H″. (Table 1) Measurements of wing bud length to width ratios (<italic>L</italic>/<italic>W</italic>) made on the plan views. (Table 2) Stages according to the original <italic>L</italic>/<italic>W</italic> measurements of ##REF##24539719##Hamburger and Hamilton (1951)##. The measurements for later stages are not included in Table 1 as they are not covered by Hamburger and Hamilton.</p></caption></fig> <fig id="fig3"><label>Fig. 3</label><caption><p>Reproducibility of 3D mapping of gene data. Photographs of right wing buds from four whole mount <italic>in-situ</italic> hybridisations for <italic>Shh</italic> from Dundee (A–D), one from Edinburgh (E) and one from Dublin (F). (A′–F′) OPT scans of the <italic>in-situs</italic> shown in panels A–F using both the fluorescence channel, for the anatomy (grey), and brightfield, for the signal (orange), visualised as volume renderings. (A″–F″) Gene expression data derived from the OPT scans in panels A′–F′ are displayed on a correspondingly labelled view of an A-Po/Pr-Di section (L) of the early HH22 wing bud (##FIG##1##Fig. 2##D), signal intensity is represented according to the heatmap in panel A″. (G) Intersect of domains of expression, in white, of panels A″–F″. (H) Intersect of domains of expression, in white, of panels A″, B″ and D″–F″, excluding the restrictive outlying data set C. (I) Mean of data in panels A″–F″. (J) Mean of data in panels A″–F″ after normalisation. (K) Median of data in panels A″–F″. The orange scale bar in all panels represents 300 μm, the scale shown in panel A″ is consistent for subsequent panels through to K. (L) A 3D model of the early HH22 wing bud showing the A-Po/Pr-Di plane of section, the nearby AER is indicated by blue arrows.</p></caption></fig> <fig id="fig4"><label>Fig. 4</label><caption><p>Mapping of gene expression for <italic>Shh</italic>, <italic>HoxD13</italic>, <italic>Fgf8</italic>, <italic>Msx1</italic>, <italic>Lmx</italic>, <italic>Wnt5a</italic> and <italic>Tbx3</italic>. (A) 3 views of the late stage HH22 reference model showing the position of planes of section used for subsequent analyses. (A.i) A dorsal view of the limb showing the position of the A-Po/Do-V plane seen in panels B and C. (A.ii) A distal view of the limb showing the position of the A-Po/Pr-Di plane seen in panel D. (A.iii) A dorsal view of the limb showing the position of the Do-V/Pr-Di plane seen in panel E. (B) A virtual section through the reference limb along the A-Po/Do-V plane shown in panel A.i with mapped gene expression patterns for <italic>Shh</italic>, <italic>HoxD13</italic>, <italic>Fgf8</italic>, <italic>Msx1</italic>, <italic>Lmx1</italic>, <italic>Wnt5a</italic> and <italic>Tbx3</italic>. Sections from the three planes shown in panel A are (C) an A-Po/Do-V plane, (D) an A-Po/Pr-Di plane, (E) a Do-V/Pr-Di plane. Pairwise comparisons were visualised on these planes for expression domains of i) <italic>Fgf8</italic> (green) and <italic>Wnt5a</italic> (red)<italic>,</italic> ii) <italic>Shh</italic> (red) and <italic>Tbx3</italic> (green), iii) <italic>Fgf8</italic> (green) and <italic>Shh</italic> (red), iv) <italic>Fgf8</italic> (green) and <italic>Msx1</italic> (red), v) <italic>Wnt5a</italic> (green) and <italic>HoxD13</italic> (red), vi) <italic>HoxD13</italic> (red) and <italic>Tbx3</italic> (green), vii) <italic>HoxD13</italic> (red) and <italic>Shh</italic> (green), viii) <italic>HoxD13</italic> (red) and <italic>Lmx1</italic> (green). Comparison numberings are consistent between panels C and E<italic>.</italic> Regions of overlap are in yellow. All scale bars (orange) represent 300 μm.</p></caption></fig> <fig id="fig5"><label>Fig. 5</label><caption><p>Improved localisation of <italic>Fgf8</italic> using fluorescent <italic>in-situ</italic> hybridisation. (A) An A-Po/Do-V plane through the HH24 reference limb (##FIG##1##Fig. 2##G′), (B) An A-Po/Pr-Di plane through the HH24 reference limb and (C) An A-Po/Pr-Di plane through the HH22 reference limb as in ##FIG##3##Fig. 4##D. All 3 panels show expression of both <italic>Fgf8</italic> (red) and <italic>Tbx3</italic> (blue) and the overlap of the domains of expression (yellow). Expression data for <italic>Tbx3</italic> are NBT–BCIP derived in all cases. <italic>Fgf8</italic> expression data in panels A and B are derived from Tyramide signal amplified FISH and NBT–BCIP derived in panel C. Note marked reduction in overlap in panels A and B compared with panel C.</p></caption></fig> <fig id="fig6"><label>Fig. 6</label><caption><p>Computational analysis of gene expression data. (A) A hierarchical clustering of the coarse sampled gene expression pattern data for both genes and limb regions. Each column represents a gene expression pattern and each cell within a column of the matrix grid represents a discrete 10 × 10 × 10 spatial volume in the reference model. Both the gene and the spatial volume data have been hierarchically clustered and trees derived showing similarity relationships for genes at the top and for spatial volumes to the left. Cells in the matrix grid are coloured according to the accompanying heatmap for signal intensity and represent the mean signal intensity for a particular gene within a particular spatial volume. This matrix does not show data for all volumes in the reference model. To the right of the matrix are blocks of colour corresponding to clusters located at the end of particularly long branches of the hierarchical clustering tree to the left, the lettered clusters correspond to the visualisations in the subsequent panels. (B–G) Visualisations of the lettered clusters from panel A in a volume rendered view of the limb, anatomy is rendered in grey and the clusters are rendered in the colours corresponding to their labelling on the right hand side of the matrix (B) purple, (C) blue, (D) green, (E) Cyan,(F) pink, (G) red. (A movie showing these clusters is available in the supplemental materials Fig. S19). Scale bars in orange correspond to 300 μm; arrows to orient the 3 principal axes, A-Po axis (blue), Pr-Di axis (red) and Do-V axis (green). Intersect of arrows represents posterior, proximal and dorsal ends of the axes.</p></caption></fig>
<supplementary-material content-type="local-data" id="d32e1381"><caption><title>Fig. S1</title><p>Suggested standard depth of colour development for OPT scanning. HH22 <italic>Wnt5a in-situ</italic> hybridisation showing suggested stain development. B) A colourstrip with the suggested endpoint stain depth highlighted in red.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="d32e1392"><caption><title>Fig. S2</title><p>This graph compares the measurements of wing bud length/width ratios from reference embryos (##FIG##1##Figs. 2##A′–F′) to the original values stipulated by Hamburger and Hamilton. The reference model measurements are shown in blue and the error bar represents the variation between the left and right wing buds of the reference embryos. The Hamilton–Hamburger measurements are shown in red and the error bars represent the range of values Hamburger and Hamilton associate with specific stages. Embryos with <italic>L</italic>/<italic>W</italic> values outside the range covered by the Hamburger Hamilton series are not shown on this graph.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="d32e1409"><caption><title>Fig. S3</title><p>Efficiency of data capture using OPT. HH Stage 22 whole mount <italic>in-situ</italic> hybridisation showing <italic>Msx1</italic> gene expression. B) Raw image data of embryo during OPT data capture. C) Processed image data showing combined brightfield (Red) and fluorescence (Green) channels.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="d32e1423"><caption><title>Fig. S4</title><p>Raw movie of embryo rotating during OPT data capture.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="d32e1432"><caption><title>Fig. S5</title><p>Movie showing rotating combined data for the brightfield (Red) and fluorescence (Green) channels.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="d32e1440"><caption><title>Fig. S6</title><p>Movie showing the rotation of the HH25 stage embryo from ##FIG##1##Fig. 2##H, the AER is readily appreciable on all limbs.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="d32e1451"><caption><title>Fig. S7</title><p>Whole mount <italic>in-situs</italic> of HH22 embryos assayed for the expression of A) <italic>Msx1</italic>, B) <italic>Lmx1</italic>, C) <italic>Wnt5a</italic>, D) <italic>Wnt3a</italic>, E) <italic>Tbx3</italic>, F) <italic>HoxD13</italic>, G) <italic>Shh</italic> and H) <italic>Fgf8</italic>.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="d32e1487"><caption><title>Fig. S8</title><p>A movie showing an initial surface render of the anatomy of the limb bud which then fades to reveal a surface rendering of the domain of <italic>Shh</italic> (Red) expression is shown, this is then joined by the domains of <italic>HoxD13</italic> (Purple), <italic>Fgf8</italic> (Green), <italic>Tbx3</italic> (Yellow), <italic>Msx1</italic> (Pink), <italic>Wnt5a</italic> (Cyan) and <italic>Lmx1</italic> (Blue) in sequence.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="d32e1517"><caption><title>Fig. S9</title><p>A movie showing serial virtual sections along the A-Po/Do-V plane going from distal to proximal through the HH22 reference limb. The expression pattern of <italic>Shh</italic> is displayed on these sections according to the heatmap in ##FIG##2##Fig. 3##A″. The dorsal bias of <italic>Shh</italic> in more proximal sections is apparent.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="d32e1534"><caption><title>Fig. S10</title><p>Asymmetry of <italic>Shh</italic> expression. A) A distal view of a HH stage 23 <italic>in-situ</italic> showing the expression of <italic>Shh</italic>. The expression is predominantly in the dorsal region of the limb with respect to the AER, indicated by the blue arrow. B) A dorsal view of a HH Stage 23 double <italic>in-situ</italic> for <italic>Fgf8</italic> and <italic>Shh</italic>, both developed with NBT–BCIP. <italic>Shh</italic> expression in the ZPA indicated with red arrow. C) A posterior view of the same HH Stage 23 double <italic>in-situ</italic> for <italic>Fgf8</italic> and <italic>Shh</italic>. The domain of <italic>Shh</italic> expression, indicated by the red arrow, is clearly shifted to the dorsal side of the limb relative to the AER as marked by the expression of <italic>FGF8</italic>.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="d32e1580"><caption><title>Fig. S11</title><p>A) Movie showing rotation of HH22 reference limb displaying expression domains of both <italic>Wnt5a</italic> (Cyan) and <italic>Fgf8</italic> (Green).</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="d32e1594"><caption><title>Fig. S12</title><p>Movie showing rotation of HH22 reference limb displaying expression domains of both <italic>Tbx3</italic> (Yellow) and <italic>Shh</italic> (Red).</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="d32e1608"><caption><title>Fig. S13</title><p>Movie showing rotation of HH22 reference limb displaying expression domains of both <italic>Shh</italic> (Red) and <italic>Fgf8</italic> (Green).</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="d32e1622"><caption><title>Fig. S14</title><p>Movie showing rotation of HH22 reference limb displaying expression domains of both <italic>Msx1</italic> (Pink) and <italic>Fgf8</italic> (Green).</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="d32e1636"><caption><title>Fig. S15</title><p>Movie showing rotation of HH22 reference limb displaying expression domains of both <italic>Wnt5a</italic> (Cyan) and <italic>HoxD13</italic> (Purple).</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="d32e1651"><caption><title>Fig. S16</title><p>Movie showing rotation of HH22 reference limb displaying expression domains of both <italic>Tbx3</italic> (Yellow) and <italic>HoxD13</italic> (Purple).</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="d32e1665"><caption><title>Fig. S17</title><p>Movie showing rotation of HH22 reference limb displaying expression domains of both <italic>Shh</italic> (Red) and <italic>HoxD13</italic> (Purple).</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="d32e1679"><caption><title>Fig. S18</title><p>Movie showing rotation of HH22 reference limb displaying expression domains of both <italic>Lmx1</italic> (Blue) and <italic>HoxD13</italic> (Purple).</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="d32e1693"><caption><title>Fig. S19</title><p>Movie showing volume renderings of domains representing clusters of gene expression. The clusters appear in alphabetical order relating to their labelling in ##FIG##4##Figs. 5##B–G and are coloured as in ##FIG##4##Fig. 5##.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="d32e1707"><caption><title>Paper Materials and methods revised</title></caption></supplementary-material>
<fn-group><fn id="d32e1711" fn-type="supplementary-material"><label>Appendix A</label><p>Supplementary data associated with this article can be found, in the online version, at <ext-link ext-link-type="doi" xlink:href="10.1016/j.ydbio.2008.01.031">doi:10.1016/j.ydbio.2008.01.031</ext-link>.</p></fn></fn-group>
<graphic xlink:href="gr1"/> <graphic xlink:href="gr2"/> <graphic xlink:href="gr3"/> <graphic xlink:href="gr4"/> <graphic xlink:href="gr5"/> <graphic xlink:href="gr6"/>
<media xlink:href="mmc1.pdf"/> <media xlink:href="mmc2.pdf"/> <media xlink:href="mmc3.pdf"/> <media xlink:href="mmc4.mov"/> <media xlink:href="mmc5.avi"/> <media xlink:href="mmc6.mpg"/> <media xlink:href="mmc7.pdf"/> <media xlink:href="mmc8.mpg"/> <media xlink:href="mmc9.mpg"/> <media xlink:href="mmc10.pdf"/> <media xlink:href="mmc11.mpg"/> <media xlink:href="mmc12.mpg"/> <media xlink:href="mmc13.mpg"/> <media xlink:href="mmc14.mpg"/> <media xlink:href="mmc15.mpg"/> <media xlink:href="mmc16.mpg"/> <media xlink:href="mmc17.mpg"/> <media xlink:href="mmc18.mpg"/> <media xlink:href="mmc19.mpg"/> <media xlink:href="mmc20.doc"/>
[{"surname": ["Bookstein"], "given-names": ["F."], "article-title": ["Principal warps: thin-plate splines and the decomposition of deformations"], "source": ["IEEE Trans. Pattern Anal. Mach. Intell."], "volume": ["11"], "year": ["1989"], "fpage": ["567"], "lpage": ["585"]}, {"surname": ["Maniatis", "Fritsch", "Sambroook"], "given-names": ["T.", "E.", "J."], "chapter-title": ["Molecular Cloning: A Laboratory Manual"], "year": ["1982"], "publisher-name": ["Cold Spring Harbor Laboratory Press"], "publisher-loc": ["New York"]}, {"surname": ["Olivo", "Izpis\u00faa Belmonte", "Tickle", "Boulin", "Duboule"], "given-names": ["J.", "J.C.", "C.", "C.", "D."], "article-title": ["Reconstruction from serial sections: a tool for developmental biology. Application to Hox genes expression in chicken wing buds"], "source": ["BioImaging"], "volume": ["1"], "year": ["1993"], "fpage": ["115"], "lpage": ["158"]}, {"surname": ["Verbeek", "den Broeder", "Boon", "Doerry", "van Raaij", "Zivkovic"], "given-names": ["F.J.", "M.J.", "P.J.B.B.", "E.", "E.J.", "D."], "article-title": ["Standard 3D digital atlas of zebrafish embryonic development for projection of experimental data"], "source": ["Proc. SPIE"], "volume": ["3964"], "year": ["1999"], "fpage": ["242"], "lpage": ["252"]}, {"surname": ["Yeh"], "given-names": ["J."], "article-title": ["The effect of miniaturized body size on skeletal morphology in frogs"], "source": ["Evol. Int. J. Org. Evol."], "volume": ["56"], "year": ["2002"], "fpage": ["628"], "lpage": ["641"]}]
Dev Biol. 2008 May 1; 317(1):13-23