3D visualization of macromolecule synthesis

  1. Timothy J Duerr
  2. Ester Comellas
  3. Eun Kyung Jeon
  4. Johanna E Farkas
  5. Marylou Joetzjer
  6. Julien Garnier
  7. Sandra J Shefelbine
  8. James R Monaghan  Is a corresponding author
  1. Department of Biology, Northeastern University, United States
  2. Department of Mathematics, Laboratori de Càlcul Numeric (LaCàN), Universitat Politècnica de Catalunya (UPC), Spain
  3. Department of Mechanical and Industrial Engineering, Northeastern University, United States
  4. University of Technology of Compiègne, France
  5. Department of Bioengineering, Northeastern University, United States
  6. Institute for Chemical Imaging of Living Systems, Northeastern University, United States
6 figures, 3 videos, 2 tables and 4 additional files

Figures

Outline of staining/analysis pipeline and exemplary images.

(A) Overview of entire sample preparation, imaging, and data analysis pipeline. (B–D) Once macromolecules are labeled in vivo, synthesis can be visualized throughout the injected animal. Here we show DNA synthesis (EdU) in the torso of a stage 52 larvae (B) and both DNA synthesis and protein glycosylation (GlcNAz) in the hand (C) and finger (D). For panels B-D, animals were pulsed with the corresponding macromolecule analog(s) for 3 hr. Images from panels B and C were uncleared and imaged at 5× magnification. Image from panel D was also uncleared and imaged at 10× magnification. Scale bars for panels B and C = 600 µm for each axis. Scale bars for panel D = 200 µm for each axis.

Figure 2 with 5 supplements
Dual staining of macromolecule synthesis in whole-mount imaging.

(A–D) Stitched and fused 3D reconstruction of 13 dpa blastemas stained for multiple macromolecules obtained by LSFM. (E–H) Single Z-plane from A-D that represents the entirety of the blastema. (I–L) Tissue section from identically treated limbs as A-H showing similar macromolecule staining patterns, indicating that the whole-mount staining method does not alter macromolecule synthesis staining patterns. Scale bars for panels A-D = 600 µm for each axis. Scale bars for panels E-L = 200 µm.

Figure 2—figure supplement 1
Color blind friendly images from Figure 2.

(A–I) Color blind friendly images of DAPI, EdU, and AHA stained limbs. (J–R) Color blind friendly images of DAPI, 5-EU, and AHA stained limbs. (S–Aa) Color blind friendly images of DAPI, EdU, and GlcNAz stained limbs. (Ab-Aj) Color blind friendly images of DAPI, 5-EU, and GlcNAz stained limbs. A-C = 800 µm for each axis. Scale bars for panels J-L, S-U, Ab-Ad = 400 µm for each axis. Scale bars for panels D-I, M-R, V-Aa, Ae-Aj = 200 µm.

Figure 2—figure supplement 2
Single staining of macromolecule synthesis in whole-mount imaging.

(A–C) Stitched and fused 3D reconstruction of 13 dpa blastemas stained for one macromolecule obtained by LSFM. (D–F) Single Z-plane from A-C that represents the entirety of the blastema. (G–I) Tissue section from identically treated limbs as A-F. (J–Aa) Grayscale images of A-I. Scale bars for panels A-C = 600 µm in each axis. Scale bars for panels D-I = 200 µm. Scale bars for grayscale images are identical to those in panels A-I.

Figure 2—figure supplement 3
Subcellular resolution obtained with LSFM.

(A) Cropped panel from Figure 2E. (B) Zoom in on panel A magnified 12.7×. (C) Zoom in on panel B magnified 6×. Scale bars = 100 µm.

Figure 2—figure supplement 4
Specificity of GlcNAz staining.

(A–D) Tissue section of a regenerating axolotl limb where the click-it cocktail for GlcNAz was applied before staining with GlcNAc antibodies (Ab). (E–H) Tissue section of a regenerating axolotl limb where GlcNAc antibodies were applied before treatment with the click-it cocktail for GlcNAz, demonstrating the specificity of the in vivo GlcNAz labeling. Scale bars = 100 µm.

Figure 2—figure supplement 5
Comparison of imaging in PBS and 67%TDE.

(A–B) Regenerating axolotl limb before (A) and after (B) clearing with 67% TDE. (C–D) Regenerating axolotl limb treated with 0.5% trypsin (C) and cleared in 67% TDE (D). (E–F) Single Z-plane of 13 dpa blastema imaged in PBS (E) or cleared and imaged in 67% TDE (F). Red indicates EdU staining whereas blue represents DAPI staining. Scale bars for panels E-F = 100 µm (G–H) Pixel intensity map of PBS imaged blastema (G) and 67% TDE imaged blastema (H). Scale bars are in units of microns.

Workflow for 3D, multiscale analysis of the regenerating axolotl humerus.

Multiscale analysis of a 35 dpa regenerating axolotl humerus, stained for AHA (red) and EdU (green). The humerus in the image stack was (A) aligned along the proximodistal (P–D) axis and (B) its morphology was segmented. The resulting mask was used to analyze the organ volume and shape by (D) reslicing it along the P-D axis and studying the cross-sections obtained. (C) The segmented morphology in B was used to mask the green channel for cellular- and molecular-level analyses. (E) Cells in the humerus were segmented and their spatial distribution was analyzed to obtain cellular number and density. (F) The masked image stack in C was resliced along the P-D axis and the pixel maps of the cross-sections were used to characterize the molecular intensity and distribution within the humerus. The vertical yellow line in A and C indicates the plane of amputation.

3D quantification across scales of a regenerating axolotl humerus.

(A) The cross-sections of the humerus in Figure 3D were analyzed with the Fiji plugin BoneJ to quantify humerus shape and volume. (A’) The cross-sectional area along the proximodistal (P–D) axis provides a measure of volume distribution along the humerus. (A”) The ratio of the maximum chord length from the minor axis (2 R1) with respect to the maximum chord length from the major axis (2 R2) provides a measure of cross-sectional circularity in the humerus. Values closer to 1.0 in the proximal side indicate a more circular cross-section in this zone. (B) The Fiji plugin Trainable Weka Segmentation 3D and 3D Objects Counter were used in the cellular analysis of proliferating chondrocytes illustrated in Figure 3E. (B’) The number of EdU+ cells within a 50 µm slice along the P-D axis was divided by the slice volume to obtain a density-like measure. (B”) The center of mass of each cell was plotted in 3D, with point size and color proportional to the segmented cell volume. (C) The molecular intensity and distribution were analyzed based on the resliced pixel intensity maps of the masked green channel in Figure 3F. (C’) Mean intensity of each slice perpendicular to the P-D axis. (C”) The histogram of the EdU staining in the humerus can provide a measure of the DNA synthesis rate. The vertical yellow line in the top row images indicates the plane of amputation.

3D visualization of DNA synthesis in innervated/denervated regenerating limbs.

(A) Schematic of experimental design used to obtain samples from B-O. (B–O) Time course of regeneration in innervated and 24 hr denervated limbs at 0, 6, 9, 12, 15, 18, and 25 dpa. Scale bars for panels B-O = 600 µm for each axis.

Figure 6 with 1 supplement
3D quantification of DNA synthesis in innervated/denervated regenerating limbs.

(A) A cube with sides of 175 µm was cropped along the proximodistal axis 250 µm from the distal tip of each blastema. P = proximal, D = distal. (B) Violin plots illustrate the pixel intensity of the innervated vs denervated blastema cubes. Comparison of mean intensity values (marked with a cross) of the same animal confirms that innervated blastemas have faster DNA synthesis rates than their denervated counterparts.

Figure 6—figure supplement 1
Pixel intensity histograms from innervated and denervated limbs.

Histograms depicting EdU pixel intensity from innervated and denervated limbs. Data are the same as in Figure 6B.

Videos

Video 1
Rotating axolotl hand stained for EdU (green) and GlcNAz (red).
Video 2
Rotating axolotl limb stained for EdU (green), AHA (red), and DAPI (blue).
Video 3
Scroll through of Z-stack from Video 2.

Tables

Table 1
Monomer analogs used to demonstrate the whole-mount visualization method.
NameMacromolecule analogBiological processClick-it modificationReference
5-ethynyl-2′-deoxyuridine (EdU)ThymidineDNA synthesisAlkyneSalic and Mitchison, 2008
5-Ethynyl Uridine (5-EU)UracilTranscriptionAlkyneJao and Salic, 2008
L-Azidohomoalanine (AHA)MethionineTranslationAzideWang et al., 2017
Azide-modified glucosamine (GlcNAz)GlucoseamineProtein glycosylationAzideLaughlin et al., 2006
Key resources table
Reagent type
(species) or resource
DesignationSource or referenceIdentifiersAdditional information
Strain, strain background (Ambystoma mexicanum)d/d axolotlAmbystoma genetic stock centerRRID:AGSC_101L
Commercial assay, kit5-Ethynyl-2’-deoxyuridine
(EdU)
clickchemistrytools.comCat# 1149Monomer analog
Commercial assay, kit5-Ethynyl Uridine (5-EU)clickchemistrytools.comCat# 1261Monomer analog
Commercial assay, kitL-Azidohomoalanine (AHA)clickchemistrytools.comCat# 1066Monomer analog
Commercial assay, kitN-azidoacetylglucosamine-tetraacylated (Ac4GlcNAz)clickchemistrytools.comCat# 1085Monomer analog
Commercial assay, kitAFDye 488 Azideclickchemistrytools.comCat# 1275Fluorescent azide
Commercial assay, kitAFDye 594 Azideclickchemistrytools.comCat# 1295Fluorescent azide
Commercial assay, kitAFDye 594 Alkyneclickchemistrytools.comCat# 1297Fluorescent alkyne
Chemical compound, drug2,2′-Thiodiethanol (TDE)SigmaCat# 166782Clearing agent
Chemical compound, drugL(+)-Ascorbic acid sodium saltThermoCat# AC352680050Click-it cocktail component
Chemical compound, drugCopper(II) sulfate pentahydrate, 98+%ThermoCat# AC423615000Click-it cocktail component
Chemical compound, drugSlowFade Gold Antifade MountantThermoCat# S36936Mountant

Additional files

Source code 1

Source code for multiscale analysis in Figures 3 and 4.

Source code used in the 3D quantification across scales of a regenerating axolotl humerus depicted in Figures 3 and 4. A tutorial describing the process in detail is provided together with the ImageJ and Matlab scripts. This material is also available on Zenodo at https://doi.org/10.5281/zenodo.3891878. The original image used as the starting point is available at Northeastern University’s Digital Repository and also upon request to the authors.

https://cdn.elifesciences.org/articles/60354/elife-60354-code1-v2.zip
Source code 2

Source code for DNA synthesis analysis in Figure 6.

Source code used in the 3D quantification of DNA synthesis in innervated/denervated regenerating limbs depicted in Figure 6. Annotated ImageJ and Matlab scripts are provided. The raw images are available at Northeastern University’s Digital Repository and also upon request to the authors.

https://cdn.elifesciences.org/articles/60354/elife-60354-code2-v2.zip
Supplementary file 1

Imaging parameters for LSFM and confocal microscopy.

Key imaging parameters for imaging whole tissues with a 5× objective in LSFM, or imaging 10 µm tissue sections with a 20× objective in confocal microscopy.

https://cdn.elifesciences.org/articles/60354/elife-60354-supp1-v2.xlsx
Transparent reporting form
https://cdn.elifesciences.org/articles/60354/elife-60354-transrepform-v2.pdf

Download links

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

Downloads (link to download the article as PDF)

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

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

  1. Timothy J Duerr
  2. Ester Comellas
  3. Eun Kyung Jeon
  4. Johanna E Farkas
  5. Marylou Joetzjer
  6. Julien Garnier
  7. Sandra J Shefelbine
  8. James R Monaghan
(2020)
3D visualization of macromolecule synthesis
eLife 9:e60354.
https://doi.org/10.7554/eLife.60354