High-throughput dual-colour precision imaging for brain-wide connectome with cytoarchitectonic landmarks at the cellular level

Abstract

The precise annotation and accurate identification of neural structures are prerequisites for studying mammalian brain function. The orientation of neurons and neural circuits is usually determined by mapping brain images to coarse axial-sampling planar reference atlases. However, individual differences at the cellular level likely lead to position errors and an inability to orient neural projections at single-cell resolution. Here, we present a high-throughput precision imaging method that can acquire a co-localized brain-wide data set of both fluorescent-labelled neurons and counterstained cell bodies at a voxel size of 0.32 × 0.32 × 2.0 μm in 3 days for a single mouse brain. We acquire mouse whole-brain imaging data sets of multiple types of neurons and projections with anatomical annotation at single-neuron resolution. The results show that the simultaneous acquisition of labelled neural structures and cytoarchitecture reference in the same brain greatly facilitates precise tracing of long-range projections and accurate locating of nuclei.

ARTICLE

Received 25 Aug 2015 | Accepted 6 Jun 2016 | Published 4 Jul 2016

Hui Gong1,2,*, Dongli Xu1,2,*, Jing Yuan1,2, Xiangning Li1,2, Congdi Guo1,2, Jie Peng1,2, Yuxin Li1,2,Lindsay A. Schwarz3, Anan Li1,2, Bihe Hu1,2, Benyi Xiong1,2, Qingtao Sun1,2, Yalun Zhang1,2, Jiepeng Liu1,2, Qiuyuan Zhong1,2, Tonghui Xu1,2, Shaoqun Zeng1,2 & Qingming Luo1,2

The precise annotation and accurate identication of neural structures are prerequisites for studying mammalian brain function. The orientation of neurons and neural circuits is usually determined by mapping brain images to coarse axial-sampling planar reference atlases. However, individual differences at the cellular level likely lead to position errors and an inability to orient neural projections at single-cell resolution. Here, we present a high-throughput precision imaging method that can acquire a co-localized brain-wide data set of both uorescent-labelled neurons and counterstained cell bodies at a voxel size of0.32 0.32 2.0 mm in 3 days for a single mouse brain. We acquire mouse whole-brain

imaging data sets of multiple types of neurons and projections with anatomical annotation at single-neuron resolution. The results show that the simultaneous acquisition of labelled neural structures and cytoarchitecture reference in the same brain greatly facilitates precise tracing of long-range projections and accurate locating of nuclei.

1 Britton Chance Center for Biomedical Photonics, Wuhan National Laboratory for Optoelectronics-Huazhong University of Science and Technology, 1037 Luoyu Road, Wuhan 430074, China. 2 Department of Biomedical Engineering, Key Laboratory for Biomedical Photonics of Ministry of Education, Huazhong University of Science and Technology, 1037 Luoyu Road, Wuhan 430074, China. 3 Howard Hughes Medical Institute and Department of Biology, Stanford University, 385 Serra Mall, Stanford, California 94305, USA. * These authors contributed equally to this work. Correspondence and requests for materials should be addressed to Q.L. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 7:12142 | DOI: 10.1038/ncomms12142 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

DOI: 10.1038/ncomms12142 OPEN

High-throughput dual-colour precision imaging for brain-wide connectome with cytoarchitectonic landmarks at the cellular level

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12142

At the cellular level, the connectome precisely annotates a comprehensive map of neural connections in the brain and signicantly increases current understanding of how

functional brain states emerge from underlying structural substrates. Mapping a cellular mouse connectome requires centimetre-scale imaging with axon resolution. Combined with physical sectioning, Li et al.1, Ragan et al.2 and Gong et al.3 pioneered whole-brain optical microscopic imaging. Neurons with different functions have different sizes, shapes and locations, and even neighbouring neurons of the same cell type differ in their morphologies and projection patterns4,5. Thus, orienting neural projections at single-cell resolution is required for both the precise annotation and accurate identication of the spatial organization of neural structures. To locate neurons and neural circuits, researchers usually map brain images to coarse axial-sampling planar reference atlases. However, individual differences at the cellular level likely lead to position errors, making the precise annotation and accurate identication of projections at single-cell resolution difcult.

Here, we develop an automatic microscopy method called the brain-wide positioning system (BPS) for dissecting and locating neural structures with cytoarchitectonic landmarks at a single-cell resolution. BPS involves a whole-brain real-time counterstain protocol for simultaneously staining cytoarchitectonic landmarks during imaging and a high-throughput multi-channel brain-wide imaging system, termed wide-eld large-volume tomography (WVT), for accelerating imaging acquisition at single-cell resolution. We obtain ve full-volume, dual-colour mouse brain data sets of cell-type-specic uorescent protein-expressing neurons and neuronal projections and subsequently reconstruct the three-dimensional (3D) ne structure of these neurons at the axon level using cellular landmarks. This standardized precision imaging approach facilitates quantitative brain-wide studies of 3D neuronal morphology, dendritic arbor and bouton distributions, axonal projection densities and distances and neuronal cell types, based on cytoarchitectonic landmarks at the cellular level. We propose that this method represents a routine tool for highly accurate and precise analysis of the cellular connectome.

ResultsBrain-wide positioning system. As a proof of principle, we revealed the spatial distribution with anatomical annotation of an 8-week-old Thy1-GFP M-line transgenic mouse6 using BPS (Fig. 1 and Supplementary Fig. 1). The resin-embedded whole-brain sample was then xed in the anterior-posterior (A-P) direction on a 3D translation stage. The WVT system employed wide-eld two-channel fast structured illumination microscopy (SIM) to accelerate imaging acquisition. SIM, which was rst demonstrated by Neil et al.7, has optical sectioning capability comparable to confocal microscopy. The system automatically performed brain-wide data acquisition. The 3D translation stage moved the sample to the microtome, and subsequently the microtome removed the supercial layer, revealing a smooth surface. The sample was then moved to fast SIM, and images were acquired through mosaic scanning in each layer. To avoid the compromised effect of sectioning marks on the machined surface, we set the imaging plane slightly below the surface, generally at 12 mm. The sectioning-imaging cycles were repeated until entire brain acquisition was completed (Fig. 1a). To achieve whole-brain counterstaining, we immersed the sample in a uorescent nuclear staining solution during whole-brain imaging. After sectioning, the dye immediately penetrated the fresh surface, combining with nucleic acid inside cell bodies and staining the soma, proximal dendrites and axon hillocks (Fig. 1b and Supplementary Figs 2, 3). To improve the detection of weak uorescence signals in thick

tissue through SIM, we modied the sample processing protocol. We added a light absorber, Sudan black B (SBB)8, to the embedded resin to inhibit background uorescence and enhance the signal-to-noise ratio of the SIM imaging (Supplementary Fig. 4). We also used chemical reactivation of the quenched uorescent protein9 on the surface of the sample to enhance the in-focus GFP signal.

Whole-brain neuron-mapping with cytoarchitecture. The whole-brain data set from a Thy1-GFP M-line transgenic mouse, including GFP-labelled neurons and propidium iodide (PI)-stained cytoarchitecture, was acquired at a voxel resolution of 0.32 0.32 1 mm in a coronal plane. The images of the

typical structure of the hippocampus, cortex and cellular layers of the cortex (Fig. 2a,b) demonstrate the method used in the present study to identify classical regions through PI-stained landmarks. Neuron arbors, axon boutons and dendritic spines were distinguished in the reconstruction of neurons (Fig. 2ce). The overlapping of soma with GFP-labelled neurons and PI-stained cellular nuclei (Fig. 2fh) illustrated single-cell co-location accuracy in the same brain using this BPS method. The sagittal reconstruction of continuous projections demonstrated that the high resolution and self-registration of BPS guarantee the data integrity of the whole-brain data set (Supplementary Fig. 5). Precision imaging through BPS not only reveals the complexity of bre orientations, even in a small tissue volume (400 400 400 mm), but also distinguishes individual axons in

the dense bre bundle (Fig. 3 and Supplementary Movie 1).

Reconstruction and localization of single pyramidal neurons. Furthermore, to illustrate the accurate co-localization of the same neuron types with cellular landmarks in the same region, we examined classical function columns and determined the complete 3D morphology of barrel cortex layer V/VI pyramidal neurons. The cellular organization of each cortical barrel column is whisker-specic10. The function, intrinsic connectivity and cellular composition of the barrel column were characterized in detail. We acquired another dual-colour whole-brain data set from Thy1-GFP M-line transgenic mouse at a voxel size of0.32 0.32 2 mm in 77 h (Supplementary Movie 2). The raw

data set was larger than 10.9 terabytes, including 4,834 coronal sections for each channel. We clearly identied regularly arranged barreloids in layer IV of the primary somatosensory cortex (S1) using perspective projections of the PI-channel image stack (Fig. 4a). The image of the PI-stained cortex illustrates a typical six-layer structure (Fig. 4b). The reconstruction of a GFP-expressing projection neuron located in layer V of the barrel eld of S1 shows the typical morphology of cortical pyramidal neurons (Fig. 4b and Supplementary Movie 3). The apical dendrites ascend through the barrel in layer IV, generating the tuft observed in layer I. Moreover, we identied and produced 3D reconstructions of all barrel columns in the mouse barrel cortex according to the PI signal and then located the reconstructed 50 layer V/VI pyramidal neurons in the barrel eld (Fig. 4c). The results showed various morphological distributions of these neurons among different barrels. The ability of BPS to continuously obtain high-quality images facilitates the tracing of individual pyramidal neurons (Supplementary Movies 4 and 5). We classied the reconstructed neurons as local (n 14) and long-range (n 36) projection neurons according

to the extent of the axonal arborizations. For long-range projection neurons, the axons extended toward the striatum or thalamus (n 16), midbrain (n 6) and pons or medulla

(n 14), and most of these neurons had other collaterals that

projected to distinct brain regions (Fig. 4d). In addition, we

2 NATURE COMMUNICATIONS | 7:12142 | DOI: 10.1038/ncomms12142 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12142 ARTICLE

a Reconstructed z-stack

Structured illuminate

b

Modulated RAW images

Record

Real-time counterstaining

Nuclei staining solution

Soma

Traverse

Block-face

Imageable depth

Uncombined dye mol.

Water mol.

Permeation

Axial

z-scan

Compute

Machined surface

Brain tissue

Dendrite

Sub-block

Axon

Spine

Resin block

Mosaic

xy-scan

Combined dye mol.

Nucleic acid mol.

IMAGING SECTIONING

Resin block

Figure 1 | Principle of BPS. (a) BPS Data acquisition overview. The resin-embedded sample moves between the SIM and the microtome. The SIM acquires three raw images to reconstruct an optical section and axially scans for Z-stack imaging over a certain range. The entire coronal section is imaged in a mosaic manner. (b) Schematic representation of the real-time cytoarchitecture staining. The dye molecules penetrate into the fresh resin-embedded sample surface and immediately combine with the nucleic acids inside the cell body. The soma is then visualized.

a

c d

e

b

f g h

Figure 2 | Brain-wide dual-colour imaging of a Thy1-GFP M-line mouse brain. (a,b) Maximum intensity projections of the coronal sections. The projections were 400 (a) and 1 mm (b) thick. (c) Merged images indicated by the white boxes in a,b. (d) An enlarged view of the area indicated by the upper white square in c, demonstrating the visualization of axonal boutons, indicated using arrowheads. (e) An enlarged view of the area indicated by the white rectangle in c demonstrating visualization of dendritic spines, indicated using arrowheads. (f,g) Raw signals of the enlarged views of the area indicated by the lower white square in c. (h) Merged image, demonstrating accurate co-localization at single soma resolution using PI staining. (c,h) The projection thicknesses of GFP and PI were 400 mm and 1 mm, respectively. Scale bars, (a,b) 1 mm; (c) 50 mm; (d) 40 mm; (e) 5 mm; and (fh) 10 mm.

classied these 50 pyramidal neurons as long-tufted (n 24) and

short-sparse (n 26) neurons (Fig. 4e). The apical dendrites of

the long-tufted neurons ascended to layer I and had abundant arbors, whereas those of the short-sparse neurons did not reach layer I and exhibited only rare branches. All local projection neurons were short-sparse neurons. Details on the localization and 3D reconstruction of these 50 neurons throughout the entire

brain are shown in Supplementary Movie 6. These results showed that the spatial relationships between pyramidal neurons and barrel columns in the mouse barrel cortex are diverse, indicating that these individual neurons might play different roles in whisker movement. Compared with traditional single-cell tracing, BPS method has two advantages. The morphology reconstructions of neurons are consecutive and complete,

NATURE COMMUNICATIONS | 7:12142 | DOI: 10.1038/ncomms12142 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12142

a b

Lateral

400 m

Ventral

400 m

400 m

Anterior

c

d

Figure 3 | Continuous neural projections with high resolution. (a) 3D reconstruction of a GFP-labelled image block without preprocessing. The image block is selected from the same whole-brain data set of Fig. 2 and located in the thalamus. (bd) Maximum intensity projections on transverse, sagittal and coronal sections, accordingly. Each projection depth is 400 mm. Scale bar, 50 mm.

including local and long-range axons. In addition, simultaneously imaging neural structure and cytoarchitecture during whole-brain imaging enables to locate the reconstructed structures according to their own cytoarchitecture.

Quantitative morphology analysis of pyramidal neurons. We also orientated and traced long-range projection neurons in the barrel cortex in the rst Thy-1 GFP data set. Consistent with the results from the second mouse brain, the axons of long-range projection neurons primarily extended toward the striatum or thalamus (n 5), midbrain (n 5), and pons or medulla (n 3).

Using these two data sets, we quantitatively compared the ne morphology of neurons projecting to the thalamus (thalamus group (THG), n 21) with those projecting to the midbrain

(midbrain group (MBG), n 11) (Fig. 5a and Supplementary

Table 1). The axonal lengths of the MBG were more than double those of the THG (Fig. 5c), whereas the numbers of axonal branches were similar between the two groups (Fig. 5d). The dendritic complexity of the neurons was signicantly different between the two groups. Both the dendrite lengths (Fig. 5e) and branch numbers (Fig. 5f) of the neurons in the MBG were much larger than those in the THG. Furthermore, we analysed and compared the morphology of the apical and basal dendrites in the two groups. The two groups had obvious differences in both apical dendrite lengths (Fig. 5g) and branch numbers (Fig. 5h).

For the basal dendrites, the lengths of the neurons in the two groups were obviously different (Fig. 5i), although there was no obvious difference in the branch numbers between the two groups (Fig. 5j). These results showed that the neurons projecting to the midbrain have longer axons and more complicated dendrite structures than the neurons projecting to the thalamus (Fig. 5a,b). Furthermore, these results indicate that in the barrel cortex, pyramidal neurons with longer axonal projections have a more complicated dendritic structure.

Mapping long-range projections with nuclei locations. The co-localization of axonal long-range projections with cellular landmarks is a prerequisite for the precise analysis of the cellular connectome. To illustrate the accurate identication of the nuclei that the axons pass through and precise tracing of single axons, we performed two-channel whole-brain imaging at a voxel size of0.32 0.32 2 mm on an 8-week-old C57BL/6J mouse injected

with adeno-associated virus vector plasmid (AAV, expressing GFP)11 in the cingulate cortex (Cg). We reconstructed GFP-labelled axons and conducted anterograde tracing of the brain-wide projections from the Cg (Fig. 6a and Supplementary Movie 7). We determined the destination nuclei according to the PI signal, and the results highlighted projection connectivities between the Cg region and the accumbens nucleus (Acb), posterior pretectal nucleus (PPT), ventral tegmental area (VTA)

4 NATURE COMMUNICATIONS | 7:12142 | DOI: 10.1038/ncomms12142 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12142 ARTICLE

a b

E1 2 3 4

5

D3

B1

C1

D1

D4

I

II/III

IV

V

VI

A1

c

d

Midbrain

C1

B1

Thalamus

P

Medulla

L

V

e

I

II/III

IV

V

VI

A1 C1 D3 D5 E3 A4 B1 B3 D5 E5

Figure 4 | Localization and reconstruction of pyramidal neurons in the Thy1-GFP M-line mouse barrel eld. (a) PI-channel volume rendering and the 220-mm-thick perspective projection of the barrel eld of the right hemisphere, showing layer IV cytoarchitecture. Scale bar, 200 mm. (b) A reconstructed pyramidal neuron and associated cytoarchitecture reference. (left) Neuron and the local layer IV cytoarchitecture of the D3 and D4 columns(410 340 80 mm image stack). The apical dendrites of this neuron are located in the barrel hollow of D3. The blue arrow indicates the location where

the neuron passes through the upper boundary of the image stack. (right) The cytoarchitecture from layer I to the corpus callosum surrounding the same neuron as shown on the left. The blue arrowheads indicate barrel walls. Image stack volume size: 760 850 150 mm. Scale bar, 100 mm. (c) 3D

reconstruction of the barrel walls and 50 representative neurons located in the barrel eld. The barrel walls were reconstructed according to PI-stained images and represented as gold loops (height: 80 mm). (d) The distribution of all long-range projection neurons among the 50 representative neurons. The thalamus, midbrain and medulla are the major projection regions. L, left; P, posterior; and V, ventral. (e) Ten reconstructed neurons (c) and corresponding barrel columns. The ve neurons shown on the left are long-tufted pyramidal neurons, and the remaining ve neurons are short-sparse neurons. The orange neuron located in D3 is the same neuron shown in b. Blue, green and orange colours represent neurons with axons extending the furthest to the striatum or thalamus, the midbrain and the pons or medulla, respectively. Pink denotes local projection neurons.

and other brain regions and such as the contralateral cortex and subcortical brain areas, consistent with previous studies12,13. Most notably, the high resolution of BPS revealed the distribution and location of the axon arbors and cells in any brain region. We identied 15 nuclei through which the projections passed (the inset of Fig. 6a). These regions received different projection patterns from the Cg, and the morphology and density of the cells varied among different areas (Fig. 6b). In particular, the BPS method obtained ne and accurate whole-brain imaging at a high spatial resolution (Fig. 6c), thereby facilitating not only the

distinction of dense GFP-labelled neurons but also the tracing of single-axon pathways. We extracted a 541 667 180 mm data

cube at the injection site and counted 746 AAV-labelled neurons (Supplementary Fig. 6 and Supplementary Movie 8). We traced an axon from the Cg to the deep white layer of the superior colliculus (DpWh) (the inset in Fig. 6c,d, Supplementary Fig. 7 and Supplementary Movies 9 and 10). This axon projected successively through eight brain regions, including corpus callosum (cc), lateral septal nucleus (LS), the nucleus of the vertical limb of the diagonal band (VDB), the horizontal limb of

NATURE COMMUNICATIONS | 7:12142 | DOI: 10.1038/ncomms12142 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12142

a b

I

II/III

IV

V

VI

Midbrain Thalamus

M

A

V

c e g i

12

Axon length (mm)

Total numberof axon branches

Dendrite length (mm)

10

8

6

4

2

Apical dendrite length (mm)

6

4

3

2

1

0

Basal dendrite length (mm)

4

3

2

1

0

***

10

***

***

5

*

8

6

4

2

0 THG MBG

0 THG MBG THG MBG THG MBG THG MBG

0 THG MBG THG MBG THG MBG

d f h j

8

6

4

2

Total number

of dendrite branches

140

100

80 60 40 20

0

***

Total number

of apical dendrite branches

100

80

60

40

20

0

of basal dendrite branches

50

40

30

20

10

0

*

120

***

Total number

Figure 5 | Quantitative comparison of the ne morphology of neurons. (a) 3D reconstruction of the barrel walls and two representative long-range projection neurons located in the barrel cortex. The barrel walls were reconstructed according to PI-stained images and represented as gold loops (height: 80 mm). The blue neuron was localized to the D5 barrel column and projected to the thalamus. The green neuron was localized in the E2 barrel column and projected to the midbrain. A, anterior; M, middle; and V, ventral. (b) The enlarged views of the dendrites of the two neurons in (a) showing the distinct dendrite complexities of the projection neurons in the two groups. (c,d) Histograms of the axonal lengths (Po0.001) and branch numbers (P 0.520) of

the long-range neurons, respectively. THG (n 15) and MBG (n 11) represent the neurons projecting to the thalamus and midbrain, respectively.

(e,f) Histograms of the dendrite lengths (Po0.001) and branch numbers (Po0.001) of long-range neurons, respectively. (g,h) Histograms of the apical dendrite lengths (Po0.001) and branch numbers (Po0.001) of long-range neurons, respectively. (i,j) Histograms of the basal dendrite lengths (P 0.017)

and branch numbers (P 0.045) of long-range neurons, respectively. We performed Students t-test on two groups by SPSS 22. Error bars were dened as

s.e.m. Condence level was set to 0.05 (P value). * represents Po0.05, and *** represents Po0.001.

the diagonal band (HDB), substantia innominata (SI), inferior colliculus (ic), lateral posterior thalamic nucleus, mediorostral part (LPMR) and the intermediate grey layer of the superior colliculus (InG), indicating that the cingulate has direct connections with brain regions involved in vision, pain and attention regulation. These results demonstrated that the BPS method has an advantage in tracing and co-localizing brain-wide dense projections.

Notably, different from previous method for axon tracing14, the high-spatial resolution of BPS enables us to perform the

microanatomical quantitative analysis of subcellular structure. To determine whether the terminal axon distribution could be revealed through BPS, we analysed the terminal axon arbors and boutons from a long-range projection neuron of the hypothalamus in the cortex (Fig. 2f). In a 128 128 400 mm

cube of layer II/III, the length of the axon arbor is 3,221.6 mm, and the number of boutons is 121. These quantitative results demonstrate that BPS not only facilitates the determination of connections among regions but also evaluates connection efciencies.

6 NATURE COMMUNICATIONS | 7:12142 | DOI: 10.1038/ncomms12142 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12142 ARTICLE

a

b

1. cc

2. CPu

3. M1/M2

4. ANT 5. MTN

6. VM/ZI

3

4

9

8

10 11

7 13

14

15

1

2 5

6

12

11. ml/VM

7. V1/V2 8. IC 9. MD 10. PAG

12. Zl/cp

13. IC 15. Ent

14. SC

c

d

Cg

cc

LPMR InG

SI

LS

DpWh

VDB

ic

HDB

Figure 6 | Anterograde tracing and localization of brain-wide neural projections. (a) Anterograde projections of AAV-GFP-labelled Cg neurons in the entire brain. The inset shows the main pattern of the Cg projecting to the whole brain, and the 15 passing regions are indicated in the global observation at low resolution. The grey sphere indicates the injection site. (b) Merged local maximum intensity projections of the sagittal planes of the corresponding positions indicated in a. Passing regions: (b1) corpus callosum, (b2) caudate putamen, (b3) primary/secondary motor cortex, (b4) anteroventral thalamic nucleus and anterodorsal thalamic nucleus, (b5) anteromedial thalamic nucleus and paratenial thalamic nucleus, (b6) ventromedial thalamic nucleus, zona incerta dorsal part (ZID), zona incerta ventral part (ZIV), and reticular thalamic nucleus, (b7) primary/secondary visual cortex, (b8) inferior colliculus (ic), (b9) mediodorsal thalamic nucleus central part, mediodorsal thalamic nucleus dorsal part and mediodorsal thalamic nucleus lateral part, (b10) periaqueductal gray, (b11) medial lemniscus and ZID, (b12) ZID/ZIV and cerebral peduncle, (b13) ic, (b14) superior colliculus and (b15) entorhinal cortex. The size of each panel is 300 320 mm. Green represents the maximum intensity projections of GFP-labelled axons. The projection thickness shown in

images 3, 5, 6, 9, 13 and 14 is 32 mm, and the projection thickness of the other images is 64 mm. Magenta represents PI-stained cells (2 mm thickness). (c) Projection pattern from Cg neurons on the coronal plane indicated with an arrow in (a). Green represents the maximum intensity projection of GFP-labelled projections (960 mm total thickness). Magenta represents PI-stained cells (2 mm thickness). The inset is an amplied image of the block shown in (c) containing sparse axons. (d) The reconstructed and localized axon indicated with arrowheads in (c). The cell body is located in the Cg region, and the axon connects to DpWh. Additionally, the axon bifurcates in the ic and LPMR brain regions. Purple represents the PI-counterstained cytoarchitecture background. Scale bar, 50 mm (b); 500 mm (c); 50 mm (the inset in c); and 100 mm (d).

In addition, we injected 90 nl of AAV-FLEX (expressing GFP)15 into the posterior region of the anterior hypothalamic area of a SOM-Cre (ref. 16) mouse and imaged the entire brain using dual-colour channels with BPS at a voxel size of0.32 0.32 2 mm. The results showed the whole-brain

projection patterns of the specic neuron type in the injected brain area (Fig. 7) and suggested that BPS in combination with different labelling methods facilitates the analysis of the locations and functions of specic neurons and neural circuits.

Moreover, we labelled brainstem neurons with long-range projections to olfactory bulb of a 7-week-old wild-type mouse. We imaged the entire brain using dual-colour channels through BPS at a voxel size of 0.32 0.32 2 mm (Fig. 8). Sparse

labelling in combination with BPS facilitates the visualization and orientation of local and long-range projections of brainstem neurons with cytoarchitecture in the whole brain (Fig. 8a and Supplementary Fig. 8). We identied the axon terminals of several neurons (Fig. 8bg). Thus, these results explicitly show that these long-range neurons project to the bilateral paraflocculus (PFI) as well, and possess large axon terminals (Fig. 8b,g). The labelled neurons were so sparse that the

detailed branches and axon terminals from a single axon could be identied (Fig. 8d).

DiscussionHere, we report a high-throughput method for the precise reconstruction and accurate localization of specic labelled neurons and projections in the whole brain. We observed the integrated cytoarchitecture of barrel columns in the barrel cortex using real-time counterstaining during whole-brain imaging and traced layer V/VI pyramidal neurons. In addition, we precisely traced a representative single-neuron projection on these PI-counterstained landmarks and identied the nuclei through which this projection passed. The diverse and complex spatial relationships between these neurons and nuclei suggest the necessity of reconstructing and analysing neuronal morphology and neural connections with accurate anatomical annotation. In this regard, BPS provides two major advantages.

First, the acquisition of cytoarchitecture landmarks in whole-brain imaging is simple and feasible. In traditional neurobiological experiments, the registration of neural images to a planar

NATURE COMMUNICATIONS | 7:12142 | DOI: 10.1038/ncomms12142 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12142

a

b

c d

a

b c d e

e

g

f

g

h

f

i

h

i

Figure 7 | Whole-brain imaging of GFP-labelled neurons in a SOM Cre-line mouse brain. (ae) Maximum intensity projections of different coronal sections showing the projection patterns with cytoarchitectonic references. The inset shows locations of images shown in ae. The projection thicknesses of GFP-signal were 60 mm (a,b), 120 mm (c) 200 mm (d) and 100 mm (e). Magenta represents PI-stained cells (2 mm thickness). (fi) Representative raw images showing the ne structure, such as axon boutons, of the long-range projections shown in (e). Scale bar, 1 mm (ae) and 15 mm (fi).

brain reference atlas is typically required. Individual differences among subjects, unavoidable deformation during tissue preparation, and interval-sampling planar reference atlases inevitably generate location errors in these forced registrations. Thus, staining the cytoarchitecture landmarks and labelling the neural structures in the same brain are necessary to avoid these errors. However, low-permeability nucleic acid dyes, such as PI and DAPI, are traditionally used in slide staining to label cytoarchitecture landmarks, and these stains are difcult to use when staining the whole mouse brain. Here, we propose a new concept of real-time counterstaining during whole-brain imaging, rather than whole-brain counterstaining. The low permeability of nucleic acid dyes results in supercial staining during whole-brain imaging without background interference from deep layers in the same channel. Moreover, the sectioning-staining-imaging cycle guarantees consistent PI staining in all coronal planes. This real-time staining approach avoids additional sample preparation. The stained cytoarchitecture and labelled neurons in the same eld of view (FOV) are simultaneously imaged. Thus, neither additional data acquisition time for anatomical landmarks nor extra two-channel registration is needed. This concept can be

applied to other mechanical section-based whole-brain imaging techniques, such as serial two-photon (STP) tomography and uorescence micro-optical sectioning tomography (fMOST), for acquiring anatomical references in the same brain.

Second, WVT efciently shortens the time required for data acquisition during full-volume whole-brain imaging at single-neuron resolution. The time required for data acquisition is one of the decisive factors for the use of whole-brain imaging as a routine approach in neuroanatomical studies. Currently, STP is considered the fastest approach; however, this approach sacrices the continuity of the whole-brain imaging data set. fMOST and 2p-fMOST (ref. 17) distinguish the detailed neural structures in full-volume whole-brain imaging, requiring more than 1 week for completion. As all these techniques employ a point-scanning approach, it is challenging to improve imaging throughput. Taking advantage of the high throughput of the wide-eld and the background inhibition of SIM, BPS achieves full-volume imaging of the mouse brain with high resolution within 3 days. The high continuity and self-registration of the whole-brain imaging data set through BPS facilitates the tracing of the detailed morphology of neurons, including axons, dendritic arbors,

8 NATURE COMMUNICATIONS | 7:12142 | DOI: 10.1038/ncomms12142 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12142 ARTICLE

a

d

c

d

f

b

g

f

e

b c e g

Figure 8 | Reconstruction of sparse-labelled brainstem neurons with long-range projections. (a) Maximum intensity projections of the coronal sections. The projection thickness was 900 mm. (bg) Enlarged views of the single axon arbor and axon terminals of the neurons shown in a. The reconstructed boutons are indicated as white arrowheads. The large axon terminals were indicated as orange arrowheads. Scale bar, 1 mm (a); 100 mm (d); 50 mm (f); and 25 mm (b,c,e,g).

boutons and spines in the whole brain. With the rapid development of an objective with a larger FOV, a camera with a larger sensor size, and a 3D stage with longer travel, we expect that the throughput of whole-brain imaging will further improve. Compared with point-scanning whole-brain imaging, the BPS approach has signicant advantages, with potential for imaging larger tissues, such as whole primate brains.

Recently, the CUBIC technique18 has been a notable breakthrough in counterstaining mouse brains. However, it remains unknown whether the chemically treated samples are sufciently hard to facilitate mechanical sectioning, such as the samples used in STP or fMOST. Thus far, CUBIC is employed in light-sheet microscopy19 to image counterstained brains and has not shown potential for use in detailed whole-brain labelling and counterstaining images at axonal resolution.

Compared with current methods, BPS technology represents the technical advance of providing co-located cytoarchitecture for neural structure in the same whole brain. In addition, this technology also has the advantages of automation, high throughput, high resolution and robustness to accelerate the acquisition of high-resolution whole-brain data sets. BPS enables us to fast acquire both the neural structures and their own anatomic reference simultaneously. It helps to accurately identify and precisely annotate the brain-wide neural structures, which is difcult to be achieved by previous methods. Potentially, BPS can be routinely used to examine the intersubject variability of neurons. Consistent with previous studies20,21, the results of the present study demonstrate that in different brains of mice that are the same strain, even the same types of neurons in the same brain areas differ in terms of morphology and distribution. Revealing the intersubject variability of an axonal projection is a fundamental requirement in neuroanatomy. Axonal projections of individual neurons correlate with sublaminar location, dendritic morphology, intrinsic and synaptic properties and local connectivity patterns of those neurons, respectively2224. The accurate tracing of axonal projections relies on the precise localization of brain regions or nuclei. Benetting from the cytoarchitecture of the brain, the brain-wide images obtained

using the BPS system can facilitate a comparison of the ne intersubject differences between axonal projections. Neuro-scientists have made signicant advances in the genetic labelling of specic neural circuits25. The combination of the genetic dissection of neural circuits and BPS technology will accelerate studies of the intersubject variability of neural circuits.

Full-volume imaging at high resolution inevitably generates large data sets, resulting in a crucial challenge to store, transfer, compute, manage and analyse these data. The most fundamental and essential challenge to mine large neural data sets is to recognize and reconstruct the complete morphology of single neurons from the raw data. This challenge is recognized worldwide as an open question and bottleneck. As a breakthrough, the BigNeuron project of the Allen Institute for Brain Science and the Blue Brain Project in Europe have focused on computational efforts. Recently, Markram et al.26 reported a simulation method of partial cortex circuits based on vast amounts of experimental and model data.

In summary, the BPS method described in the present study facilitate the acquisition of a more detailed morphology of neurons and accurately identify the brain regions or nuclei in which these are located, connected with and passing through. We propose that this method has many potential applications for neuroscience and will promote future cellular connectome studies using a traditional interval-sampling 2D atlas to precisely generate personalized 3D maps for each brain at cellular resolution. This method could also contribute to cell type, projectome and connectome studies.

Methods

Animals. All the animal experiments followed procedures approved by the Institutional Animal Ethics Committee of Huazhong University of Science and Technology or the Administrative Panel on Laboratory Animal Care (APLAC) at Stanford University. An 8-week-old C57BL/6J male mouse, two 8-week-old Thy1-GFP M-line transgenic mice (Jackson Laboratory, Bar Harbor, ME, USA),a 12-week-old SOM-Cre transgenic mouse (Huang Lab, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA) and a 7-week-old wild-type mouse (mixed genetic strain) were used.

NATURE COMMUNICATIONS | 7:12142 | DOI: 10.1038/ncomms12142 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12142

Tissue preparation. All histological procedures have been previously described3,17,27. Briey, the mice were anaesthetized using a 1% solution of sodium pentobarbital and subsequently intracardially perfused with 0.01 M PBS (Sigma-Aldrich Inc., St Louis, MO, USA), followed by 4% paraformaldehyde (Sigma-Aldrich Inc., St Louis, MO, USA) and 2.5% sucrose in 0.01 M PBS. The brains were excised and post-xed in 4% paraformaldehyde at 4 C for 24 h. After xation, each intact brain was rinsed overnight at 4 C in a 0.01 M PBS solution and subsequently dehydrated in a graded ethanol series (50, 70 and 95% ethanol, changing from one concentration to the next every 1 h at 4 C). We modied the previous resin-embedding approach3,10,19 to inhibit background uorescence (Supplementary Fig. 4). Briey, after dehydration, the brains were immersed in a graded glycol methacrylate (GMA) series (Ted Pella Inc., Redding, CA, USA), including0.2% SBB (70%, 85%, and 100% GMA for 2 h each and 100% GMA overnight at 4 C). Subsequently, the samples were impregnated in a prepolymerization GMA solution for 3 days at 4 C and embedded in a vacuum oven at 48 C for 24 h. The 100% GMA solution comprised 67 g of A solution, 2.8 g of deionized water, 29.4 g of B solution, 0.2 g of SBB and 0.6 g of AIBN as an initiator. The 70% and 85% GMA solutions (wt wt 1) were prepared from 95% ethanol and 100% GMA.

Instrument. The WVT system is shown in Supplementary Fig. 1a. A mercury lamp (X-Cite exacte, Lumen Dynamics, Mississauga, Ontario, Canada) was used as the light source, providing high exibility of wavelength selection. The collimated excitation light was directed to a digital micro-mirror device (DMD, XD-ED01N, X-digit, Shanghai, China) to generate the illumination grid pattern. The DMD comprised 1,024 768 electronically controllable micro-mirrors, and the pitch size

of each micro-mirror was 13.68 mm. When the micro-mirror was set to on, the illumination light was reected into the light path. In contrast, in the off state, the micro-mirror blocked the illumination light. The modulated illumination light was transmitted through a lens (f 150 mm) and a water immersion objective

(1.0 NA, XLUMPLFLN 20XW, Olympus, Shinjuku, Tokyo, Japan) and focused on the sample. The uorescent light from the sample was collected through the objective and detected using two scientic complementary metal-oxide semiconductor cameras (ORCA-Flash 4.0, Hamamatsu Photonics K.K., Hamamatsu, Japan). The sensor array of the camera was 2,048 2,048 pixels with a 6.5 mm pixel

size. A piezoelectric translational stage (P-725 PIFOC Long-Travel Objective Scanner, E-753 Digital Piezo Controller, PI GmbH, Karlsruhe, Germany) moved the objective for axial scanning. The actual imaging format was set to 1,700 1,800

pixels to t the size of the modulated light eld of the DMD. The resin-embedded whole-brain sample was subsequently xed using cyanoacrylate in a 200

90 53 mm sample box. The sample box was screwed onto a high-precision 3D

translation stage (ABL20020-ANT130-AVL125, Aerotech Inc., Pittsburgh, PA, USA). The 3D translation stage moved the sample for mosaic scanning and sectioning. The travel ranges were 200, 60 and 25 mm along the x, y and z axes, respectively. The encoded resolution of the x, y and z axes was 10, 1 and 4.5 nm, respectively. The repeatability of the x, y and z axes was 200, 75 and300 nm, respectively.

Sectioning and imaging were split in the WVT system, different from micro-optical sectioning tomography (MOST) and fMOST. For WVT, the microtome was used to only remove the imaged surface. The microtome was based on a 45 diamond knife (Diatome AG, Nidau, Switzerland), similar to MOST and fMOST.

We controlled the 3D stage motion for sectioning and mosaic scanning using a computer. Running data acquisition and controlling other hardware in the WVT were achieved using customized C acquisition software in a workstation. The

acquisition and motion control software were communicated via a TCP/IP protocol. The acquired data were saved as TIFF les in a storage array (PowerVault MD1200 with a PERC H810 Host-RAID adapter, Dell Inc., Round Rock, Texas, USA).

We imaged 0.2 mm FluoSpheres Carboxylate-Modied Microspheres (Molecular Probes, Eugene, OR, USA) in the green and red channels to examine the point spread function of the WVT system. After setting the pattern period on the DMD plane to 109.44 mm (13.68 mm pixel 1 8 pixels), a 3D image stack of

the beads was acquired with a z step of 200 nm. We reconstructed the xz cross-section views of a selected bead in both channels, as shown in Supplementary Fig. 1b,c. The lateral and axial full-width at half-maximum values were 0.55 and2.20 mm in the green channel and 0.62 and 2.59 mm in the red channel, respectively (Supplementary Fig. 1dg).

Penetration of uorescent nuclear staining solution. We estimated and compared the penetration performance of PI and 40, 6-diamidino-2-phenylindole (DAPI) in the resin-embedded sample to obtain an optimized penetration protocol for the uorescent nuclear staining solution. First, we studied the PI staining effect in GMA resin-embedded brain tissue at different concentrations and selected2 mg ml 1 for use in whole-brain imaging. Then we acquired time-dependent images of the same mosaic at the depth of the imaging plane for this concentration (Supplementary Fig. 2) using the WVT. The system sectioned and imaged at a randomly selected single mosaic. The PI molecules penetrated the fresh surface immediately after sectioning. The time interval between sectioning and imaging of this mosaic was 15 s. Subsequently, we imaged this mosaic at 30 s intervals. The results indicated that the penetration of the PI solution was fast enough to facilitate staining in real time at the depth of the imaging plane. Thus, the data acquisition

time was only affected through the imaging and sectioning speed and coronal plane size, not the penetration time. We also examined the features of DAPI staining and observed that the penetration of DAPI solution was too slow to achieve real-time staining. Thus, PI was used as the cytoarchitecture dye in the present study, and the staining parameters of PI were further optimized.

Whole-brain imaging with real-time PI staining. The sample was immersed in a water bath containing PI. Whole-brain imaging was performed in the water bath. To avoid the compromising effects of sectioning marks on the machined surface, the imaging plane was set slightly below the surface, generally at 12 mm. Before data acquisition, we focused on the top surface of the sample and adjusted the objective to obtain a clear image. Subsequently, we moved the objective down to set the imaging plane below the surface of the sample block. After these preparations, the imaging parameters were set, and the WVT system automatically performed the sectioning and imaging to complete the brain-wide data acquisition. We also exibly adjusted the imaging parameters, such as range of interest, exposure time and so on. To enhance the in-focus GFP signal, we added Na2CO3 into the water bath9. Most of the GFP molecules were preserved in a nonuorescent state, rather than directly damaged, through chromophore protonation during the resin-embedding procedure. These uorescent signals were chemically recovered to the uorescent state using 0.05 M Na2CO3 during imaging.

Sectioning was achieved through a relative motion between the xed diamond knife and the 3D translation stage in the WVT. The x axis of the translation stage was the sectioning direction, and the sectioning width was 2 mm. The sectioning thickness was exible. The y and z axes provided the necessary additional movement to cover the sample surface for sectioning. Thus, a smooth fresh surface was exposed for imaging.

Subsequently, the sample was imaged. PI solution served as the immersion liquid for the objective lens. The liquid level of the solution was maintained at a level higher than the bottom surface of the objective lens during data acquisition. The GFP and PI molecules were simultaneously excited, and the emitted uorescence signals were divided using a dichroic mirror and detected using two cameras. Three phase-shifted raw images with a phase step of p/2 were required to obtain an optical section image for each imaging channel28. The axial background inhibition was achieved through optical sectioning using SIM. When necessary, axial scanning was subsequently executed using the piezoelectric translational stage. Subsequently, the sample was moved to the next mosaic FOV with a 10 mm overlap between the adjacent FOVs. The mosaic imaging process was repeated until the entire coronal section was acquired. In addition, we used a recirculating ltration device1 to maintain a attened section, remove the cutting chips, purify the PI solution and maintain a uniform PI concentration. The actual contrast of cytoarchitecture staining in the whole-brain data set showed that there was no obvious decay of the PI concentration during whole-brain imaging (Supplementary Fig. 3).

To reduce the data volume, the reconstruction of the SIM image was executed online using acquisition software. Only the reconstructed images were saved. The data set, comprising 1,700 1,800 pixel-sized mosaics, was saved at a 16-bit depth

in the Lempel-Ziv-Welch (LZW) compression TIF format. After collection, the data set was sent toa PB-sized distributed storage via a standard 10-gigabit bre.

We acquired the data set in Fig. 3 after sectioning at a 1 mm thickness and imaging at a 0.32 0.32 1 mm voxel size. To conserve the data acquisition time,

we tried another acquisition scheme at a 0.32 0.32 2 mm voxel size for all other

data acquisitions. We acquired a z-stack of two images in each FOV after axially scanning at a step of 2 mm and subsequently sectioning at a 4 mm thickness. The results (Figs 4 and 68) demonstrated that this acquisition scheme shortened the data acquisition time without sacricing the axon resolution power. In addition, we also changed the mosaic scanning range according to the prole characteristics of the mouse brain to reduce useless images of the surrounding embedding medium. We achieved the rapid data acquisition of a dual-colour whole-brain data set in 77 h (Fig. 4). The data set included 561,440 mosaic images in 4,834 layers. The system performed sectioning 2,417 times and imaged 280,720 z-stacks. Sectioning lasted for a total of 24.4 s to section column by column. The imaging time of each z-stack was determined as the exposure time of the camera, phase-shift time of DMD, axial scanning time of the objective, settle time of the 3D stages and the time for online data processing. In this case, the imaging time for each z-stack was 780 ms on average. The data acquisition time was variable for different labelling technologies.

Virus injections. AAV-CAG and AAV-CAG-FLEX expressing GFP (Serotype: 9, UNC Gene Therapy Center Vector Core, Chapel Hill, NC, USA) were used as anterograde tracers. The stereotaxic coordinates for the target areas were based on the Mouse Brain in Stereotaxic Coordinates Atlas29. Using a pressure injector (Nanoject II; Drummond Scientic Co., Broomall, PA, USA), 90 nl of AAV-CAG was injected into the Cg of an 8-week-old C57BL/6J mouse (0.38 mm A-P, 0.3 mm medial-lateral (M-L) and 1.8 mm dorsal-ventral (D-V)) and 90 nl of AAV-CAG-FLEX was injected into the posterior region of the anterior hypothalamic area of a SOM-Cre mouse (0.94 mm A-P, 0.4 mm M-L and 5 mm D-V). 0.7 ml of CAV-Cre and 0.5 ml AAV-FLEx(loxP)-TVA-GFP were injected into the olfactory bulb (from bregma: 0.75 mm L, 4 mm A and 1 mm V) and the ipsilateral locus coeruleus

10 NATURE COMMUNICATIONS | 7:12142 | DOI: 10.1038/ncomms12142 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12142 ARTICLE

(from lambda: 0.8 mm L, 0.8 mm P and 3.2 mm V) of a 7-week-old wild-type mouse, respectively. After surgery, the animals were returned to standard living conditions for 21 days until they were sacriced for brain sample preparation.

Image preprocessing. Image preprocessing for mosaic stitching (Supplementary Fig. 9) and illumination correction30 were performed for both the GFP and PI channels. The mosaics of each coronal section were stitched to obtain an entire section based on accurate spatial orientation and neighbouring overlap. Lateral illumination correction was performed section by section. Correction coefcient along each direction was determined by calculating mean intensity along each direction and corresponding polynomial curves tting. Equalizing the brightness of the different coronal sections was performed for axial illumination correction by quantifying the average grey-scale values of the images. Image preprocessing was implemented in C and optimized in parallel using the Intel MPI Library

(v3.2.2.006, Intel, Santa Clara, CA, USA). Image preprocesses for mouse brain data set at the voxel resolution of 0.32 0.32 2 mm and 0.32 0.32 1 mm were

executed on a computing server (72 cores, 2 GHz per core) within 6 and 12 h, respectively. All full coronal sections were saved at an 8-bit depth in LZW compression TIFF format after image preprocessing.

Visualization and reconstruction. We visualized the data set using Amira software (v 5.2.2, FEI, Mrignac Cedex, France) to generate the gures and Movies. The preprocessed data set was imported into Amira software using a desktop graphical workstation (T7600 with two Intel E5-2687w CPUs, 256 GB memory and an Nvidia K6000 graphics card, Dell Inc., Round Rock, Texas, USA). To process the TB-sized data on a single workstation, we transformed the data format from TIFF to the native LDA type using Amira. The visualization process included extracting the data in range of interest, sampling or interpolation, reslicing the images, identifying the maximum intensity projection, volume and surface rendering, and generating the Movies using the main module of Amira. The segmentation editor module of Amira was utilized for the manual outline segmentation of whole mouse brain.

We applied the lament editor module of Amira to brain-wide tracing of long-range axons in 3D by human-machine interaction. From soma, we continuously loaded data blocks along axons and dendrites into Amira. We assigned initial and terminal points of bres in the loaded block, and then Amira automatically calculated the path that the bre walked along between these two points. This procedure was repeated until the whole neural morphology was reconstructed. The tracing results with original position information were saved in SWC format.

We determined the barrel columns according to the counterstained images. The original PI-labelled data were sampled to 2 2 2 mm. According to the Mouse

Brain in Stereotaxic Coordinates Atlas29, we approximately localized the barrel cortex and resliced the sections parallel to the barrel cortex. We obtained a series of minimum intensity projections and subsequently identied the barrel columns in layer IV. According to the position information in the SWC le, we identied the corresponding coronal sections and obtained the maximum intensity projection sequences.

Besides, we used the NeuroGPS31 software to locate the soma in a chosen data block, as shown in Supplementary Movie 8.

Because recognition using the human eye is the golden standard of image segmenting in image processing3235, three persons back-to-back checked all results in this manuscript to guarantee the quality of the segmented and traced neurons.

Statistics. We compared the lengths and branch numbers of the axons, apical dendrites and basal dendrites of neurons projecting to the thalamus (THG) and midbrain (MBG). All parameters were measured using Neurolucida Explorer (MBF Bioscience, Williston, VT, USA). We performed Students t-test on two groups using SPSS software (v 22, IBM, New York, USA). The condence level was set to0.05 (P value), and all results are presented as the meanss.e.m.

Data availability. The authors declare that the data supporting the ndings of this study are available within the article and its Supplementary Information les.

References

1. Li, A. et al. Micro-optical sectioning tomography to obtain a high-resolution atlas of the mouse brain. Science 330, 14041408 (2010).

2. Ragan, T. et al. Serial two-photon tomography for automated ex vivo mouse brain imaging. Nat. Methods 9, 255258 (2012).

3. Gong, H. et al. Continuously tracing brain-wide long-distance axonal projections in mice at a one-micron voxel resolution. Neuroimage 74, 8798 (2013).

4. Wickersham, I. R. et al. Monosynaptic restriction of transsynaptic tracing from single, genetically targeted neurons. Neuron 53, 639647 (2007).

5. Hezel, M., Ebrahimi, F., Koch, M. & Dehghani, F. Propidium iodide staining: a new application in uorescence microscopy for analysis of cytoarchitecture in adult and developing rodent brain. Micron 43, 10311038 (2012).

6. Feng, G. et al. Imaging neuronal subsets in transgenic mice expressing multiple spectral variants of GFP. Neuron 28, 4151 (2000).

7. Neil, M. A. A., Juskaitis, R. & Wilson, T. Method of obtaining optical sectioning by using structured light in a conventional microscope. Opt. Lett. 22, 19051907 (1997).

8. Sun, Y. et al. Sudan black b reduces autouorescence in murine renal tissue. Arch. Pathol. Lab. Med. 135, 13351342 (2011).

9. Xiong, H. et al. Chemical reactivation of quenched uorescent protein molecules enables resin-embedded uorescence microimaging. Nat. Commun. 5, 3992 (2014).

10. Woolsey, Thomas A. & Loos, H. V. d. The structural organization of layer IV in the somatosensory region (SI) of mouse cerebral cortex: the description of a cortical eld composed of discrete cytoarchitectonic units. Brain Res. Bull. 17, 205242 (1970).

11. Chamberlin, N. L., Du, B., de Lacalle, S. & Saper, C. B. Recombinant adeno-associated virus vector: use for transgene expression and anterograde tract tracing in the CNS. Brain Res. 793, 169175 (1998).

12. Oh, S. W. et al. A mesoscale connectome of the mouse brain. Nature 508,

207214 (2014).

13. Allen Institute for Brain Science. Allen Mouse Brain Connectivity Atlas. Available at: http://connectivity.brain-map.org

Web End =http://connectivity.brain-map.org (2015).

14. Kita, T. & Kita, H. The subthalamic nucleus is one of multiple innervation sites for long-range corticofugal axons: a single-axon tracing study in the rat.J. Neurosci. 32, 59905999 (2012).15. Schnutgen, F. et al. A directional strategy for monitoring Cre-mediated recombination at the cellular level in the mouse. Nat. Biotechnol. 21, 562565 (2003).

16. Taniguchi, H. et al. A resource of Cre driver lines for genetic targeting of GABAergic neurons in cerebral cortex. Neuron 71, 9951013 (2011).

17. Zheng, T. et al. Visualization of brain circuits using two-photon uorescence micro-optical sectioning tomography. Opt. Express 21, 98399850 (2013).

18. Susaki, E. A. et al. Whole-brain imaging with single-cell resolution using chemical cocktails and computational analysis. Cell 157, 726739 (2014).

19. Dodt, H.-U. et al. Ultramicroscopy: three-dimensional visualization of neuronal networks in the whole mouse brain. Nat. Methods 4, 331336 (2007).

20. Jiang, X. et al. Principles of connectivity among morphologically dened cell types in adult neocortex. Science 350, aac9462 (2015).

21. Tasic, B. et al. Adult mouse cortical cell taxonomy revealed by single cell transcriptomics. Nat. Neurosci. 19, 335346 (2016).

22. Brown, S. P. & Hestrin, S. Intracortical circuits of pyramidal neurons reect their long-range axonal targets. Nature 457, 11331136 (2009).

23. Groh, A. et al. Cell-type specic properties of pyramidal neurons in neocortex underlying a layout that is modiable depending on the cortical area. Cereb. Cortex 20, 826836 (2010).

24. Oberlaender, M. et al. Three-dimensional axon morphologies of individual layer 5 neurons indicate cell type-specic intracortical pathways for whisker motion and touch. Proc. Natl Acad. Sci. USA 108, 41884193 (2011).

25. Huang, Z. J. Toward a genetic dissection of cortical circuits in the mouse. Neuron 83, 12841302 (2014).

26. Markram, H. et al. Reconstruction and simulation of neocortical microcircuitry. Cell 163, 456492 (2015).

27. Yang, Z., Hu, B., Zhang, Y., Luo, Q. & Gong, H. Development of a plastic embedding method for large-volume and uorescent-protein-expressing tissues. PLoS ONE 8, e60877 (2013).

28. Fukano, T. & Miyawaki, A. Whole-eld uorescence microscope with digital micromirror device: imaging of biological samples. Appl. Opt. 42, 41194124 (2003).

29. Paxinos, G. & Franklin, K. B. The Mouse Brain in Stereotaxic Coordinates (Academic Press, 2001).

30. Ding, W. et al. Automatic macroscopic density artefact removal in a Nissl-stained microscopic atlas of whole mouse brain. J. Microsc. 251, 168177 (2013).

31. Quan, T. et al. NeuroGPS: automated localization of neurons for brain circuits using L1 minimization model. Sci. Rep. 3, 1414 (2013).

32. Peng, H. et al. BigNeuron: large-scale 3D neuron reconstruction from optical microscopy images. Neuron 87, 252256 (2015).

33. Meijering, E. Neuron tracing in perspective. Cytometry A 77A, 693704 (2010).34. Gillette, T. A., Brown, K. M. & Ascoli, G. A. The DIADEM metric: comparing multiple reconstructions of the same neuron. Neuroinform 9, 233245 (2011).

35. Ferreira, T. A. et al. Neuronal morphometry directly from bitmap images. Nat. Methods 11, 982984 (2014).

Acknowledgements

We thank the MOST group members of Britton Chance Center for Biomedical Photonics for assistance with experiments and comments on the manuscript. We are also grateful toZ. Josh Huang and Haohong Li for providing the SOM-Cre mouse strain, Xiaojun Wang

NATURE COMMUNICATIONS | 7:12142 | DOI: 10.1038/ncomms12142 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 11

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12142

for valuable assistance with the injection experiments and Liqun Luo for constructive comments. This work was nancially supported through The 973 project no. 2015CB755602, NSFC projects nos 61421064, 91432105, 91432116, 91232000 and the director fund of the WNLO.

Author contributions

Q.L. and H.G. conceived of and designed the study and wrote the manuscript. D.X., J.Y., A.L., H.G. and Q.L. constructed the sectioning and microscope. D.X., J.Y. and Q.Z. performed data acquisition. H.G., X.L., L.A.S., B.H., Q.S. and S.Z. designed the tissue preparation and built the virus tracing samples. H.G., C.G., J.P., Y.L., A.L., B.X., Y.Z. and J.L. analysed the data, performed image processing and visualization. J.Y., T.X. and H.G. modied the manuscript. All authors discussed the conceptual and practical implications of the methods.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Gong, H. et al. High-throughput dual-colour precision imaging for brain-wide connectome with cytoarchitectonic landmarks at the cellular level.

Nat. Commun. 7:12142 doi: 10.1038/ncomms12142 (2016).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

r The Author(s) 2016

12 NATURE COMMUNICATIONS | 7:12142 | DOI: 10.1038/ncomms12142 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications