Abstract
A fundamental strategy in sensory coding is parallel processing, whereby unique, distinct features of sensation are computed and projected to the central target in the form of submodal maps. It remains unclear, however, whether such parallel processing strategy is employed in the main olfactory system, which codes the complex hierarchical odor and behavioral scenes. A potential scheme is that distinct subsets of projection neurons in the olfactory bulb (OB) form parallel projections to the targets. Taking advantage of the observation that the distinct projection neurons develop at different times, we developed a Cre-loxP-based method that allows for birthdate-specific labeling of cell bodies and their axon projections in mice. This birthdate tag analysis revealed that the mitral cells (MCs) born in an early developmental stage and the external tufted cells (TCs) born a few days later form segregated parallel projections. Specifically, the latter subset converges the axons onto only two small specific targets, one of which, located at the anterolateral edge of the olfactory tubercle (OT), excludes widespread MC projections. This target is made up of neurons that express dopamine D1 but not D2 receptor and corresponds to the most anterolateral isolation of the CAP compartments (aiCAP) that were defined previously. This finding of segregated projections suggests that olfactory sensing does indeed involve parallel processing of functionally distinct submodalities. Importantly, the birthdate tag method used here may pave the way for deciphering the functional meaning of these individual projection pathways in the future.
Significance Statement
Odorant receptors form an ordered odorant map on the main olfactory bulb (MOB). This spatial representation disappears in most subsequent targets by “diffuse divergence and random convergence” of OB projections. We revisited these projections using a novel method that can genetically dissect distinct subsets of OB projection neurons based on their neuronal birthdates. Our birthdate tag analysis exposed parallel segregated projections formed by early-born mitral cells (MCs) and late-born external tufted cells (TCs) in otherwise apparently random olfactory networks. The results show that the birthdate tag approach is effective for uncovering a hidden circuit structure formed by intermingling neuronal subsets in the olfactory as well as other nervous systems.
Introduction
Parallel processing is a fundamental strategy in sensory nervous systems (Young, 1998). In the visual system, for example, distinct submodalities such as color, orientation, or direction, are extracted from the visual scene by unique parallel computations and represented as individual higher-order maps in the central target (Nassi and Callaway, 2009). The olfactory system consists of two complete parallel pathways; the physically separated nasal and vomeronasal organs detect psychochemically distinct stimuli, which are in turn conveyed through segregated relay stations, the main olfactory bulb (MOB) and accessory olfactory bulb (AOB), to distinct sets of central targets (Dulac and Wagner, 2006). Within each pathway, a few subpathways are found peripherally between sensory neurons and the OB (Mori and Yoshihara, 1995; Brechbühl et al., 2008; Dewan et al., 2013; Greer et al., 2016). It remains unknown whether these subpathways extend parallel processing beyond the OB.
The MOB processes information of the complex odor scene through a large repertoire of odorant receptors. The nasal sensory neurons that express the same single odorant receptor converge their axons onto a few fixed glomeruli in the MOB, thereby representing the stereotypical odor map (Mori and Sakano, 2011). The projection neurons in the MOB then relay the output of single glomerulus to central target areas. Anterograde axon tracing of projection neurons assigned to the identified glomerulus revealed what appeared to be random, divergent projections of their axon terminals covering all of approximately ten target areas (Ghosh et al., 2011; Sosulski et al., 2011). Retrograde trans-synaptic tracing and postsynaptic electrophysiological recordings indicated that each central target neuron receives a combinatorial integration of outputs from multiple random glomeruli (Davison and Ehlers, 2011; Miyamichi et al., 2011). As a result, the MOB topography is undetectable other than in exceptional areas such as the pars externa (pE), a subdivision of the anterior olfactory nucleus (AON; Schoenfeld and Macrides, 1984; Yan et al., 2008) and parts of the amygdala, in which dorsal or ventral glomerular output is overrepresented (Miyamichi et al., 2011; Inokuchi et al., 2017). Instead, olfactory information processing is generally characterized by diffuse divergence from a single glomerulus to all the target areas and random convergence from multiple glomeruli to each target neuron (Wilson and Sullivan, 2011), making it challenging to deduce the logic underlying the information processing.
A critical question in understanding the olfactory logic relates to the connection specificity of projection neuron subtypes. The projection neurons are classified into mitral cells (MCs) and tufted cells (TCs) based on the location of their cell bodies in MOB layers. The information of each glomerulus is relayed by ∼20 MCs and ∼50 TCs to the central targets (Nagayama et al., 2014). TCs are a heterogeneous group, and further subclassified into internal, middle TCs (imTCs) and external TCs (eTCs) subtypes by their locations within the external plexiform layer (EPL) and the glomerular layer (GLL). These projection neuron subtypes have distinct intrabulbar connections (Mori et al., 1983; Orona et al., 1984) and electrophysiological properties (Nagayama et al., 2004; Griff et al., 2008a; Fukunaga et al., 2012), suggesting that they could carry different kinds of information. A few studies have reported that TCs project only to the anterior olfactory targets involving the AON, olfactory tubercle (OT), and anterior piriform cortex (PC; Haberly and Price, 1977; Scott, 1981; Nagayama et al., 2010). In particular, single cell dye labeling has revealed that the terminal branches of several imTCs accumulate around the boundary between the OT and anterior PC (Igarashi et al., 2012). However, since there are no good molecular markers to distinguish the MOB projection neurons, it remains unclear whether these heterogeneous populations simply make a convergent map or segregated maps in the olfactory target areas.
A potential approach to dissect the OB projection neurons is to make use of differences in their birthdates. Neuronal birthdating using nucleotide analogs is a well-established standard for the classification of neuronal subtypes (Bayer and Altman, 1987). It has led to many important discoveries while undergoing very little technological improvement. The obvious limitation is that it remains an entirely histologic method, capable only of visualizing cell nuclei that have incorporated the analogs in fixed samples.
This study developed a novel birthdating method by which Cre-loxP recombination is induced in a neuronal birthdate-dependent manner. Using this method, we demonstrate the segregation of axon trajectories of OB projection neuron subsets into parallel projections. Specifically, a subset of eTCs birthdate tagged at embryonic day (E)15.5 projected to only two small targets deviated from other projections. Their selective target in the OT was then analyzed in more details.
Materials and Methods
Mouse breeding and experiments
All animal procedures were performed in accordance with the animal care committee's regulations at the National Institute of Genetics and RIKEN Center for Biosystems Dynamics Research (BDR). Although both male and female mice were used for generation of experimental materials, datasets were obtained from only male mice after their sexes were externally apparent around one week of age. The sexes of younger mice were undetermined, but it is highly likely that only male mice were used for sampling because the transgene of the birthdate tag driver Nerog2CreER(G2A) seems to be located in the Y chromosome.
The day on which a vaginal plug was detected was designated as E0.5, and the day of birth was designated as postnatal day (P)0. Tamoxifen (TM) treatment was performed by intraperitoneally injecting 250 μl of corn oil (Sigma-Aldrich catalog #C8267) containing 9 mM TM (Sigma-Aldrich catalog #T5648) and 5 mM progesterone (Wako catalog #161-14531) into a staged pregnant mouse. Because TM often delays delivery, when pups were not born by E19.5, they were collected by caesarian delivery and given to ICR foster mothers.
Mouse lines
The mouse lines used in this study are listed in Table 1. The birthdate tag driver line Neurog2CreER(G2A) (accession number CDB0512T−1: http://www2.clst.riken.jp/arg/TG%20mutant%20mice%20list.html) was established using a genomic BAC clone RP23-360O22 encoding the mouse neurogenin (Neurog)2 gene obtained from the BACPAC Resource Center (Children’s Hospital Oakland Research Institute, Oakland, CA). The entire coding sequence in exon2 was replaced by a TM-inducible Cre recombinase, CreERT2 (Feil et al., 1997; a generous gift from Dr. Pierre Chambon, University of Strasbourg) through homologous recombination as described elsewhere (Oginuma et al., 2008). The BAC recombinant was injected into fertilized eggs with a C57BL/6 background, and the resulting mice were assayed for integration of the transgene by PCR of the genomic DNA. The three mouse lines that independently received the transgene integration were then assayed for the Cre recombinase activities by crossing with ROSA26R Cre reporter mice (Soriano, 1999). For the assay, the pregnant mice were injected with TM solution at E12.5 or E14.5, and brains of the embryos were dissected out at E18.5–E19.5 and whole-mount stained with X-gal (5-bromo-4-chloro-3-indoyl-β-D-galactopyranoside) as described elsewhere (Saga et al., 1992). Through this screen, the birthdate tag line Neurog2CreER(G2A) was selected by the highest recombination rate in OB neurons. Transgene containment was routinely determined with PCR using internal Cre recombinase primers, 5′-TAAAGATATCTCACGTACTGACGGTG-3′ and 5′-TCTCTGACCAGAGTCATCCTTAGC-3′, resulting in amplification of 300-bp fragments.
The Cdhr1tTA mouse line (accession number CDB0535T: http://www2.clst.riken.jp/arg/TG%20mutant%20mice%20list.html) was established using a mouse genomic BAC clone RP23-55P20 (BACPAC Resource Center), which encodes cadherin related family member (Cdhr)1 gene that exhibits OB-specific expression (Nagai et al., 2005). At 25 bp upstream of the translation initiation site of this gene, a loxP-LacZ-loxP-tTA cassette was inserted through homologous recombination. In this cassette, the LacZ gene followed by a transcription termination site was flanked by loxP sequences and used as a STOP sequence, which followed by the tetracycline-controlled transactivator (tTA) gene in an optimized form for mammalian expression (Inamura et al., 2012; a generous gift from Dr. Kazuhiro Ikenaka, National Institute for Physiologic Sciences). The original loxP site that had existed in the vector backbone (pBACe3.6) was deleted by replacement using p23loxZeo (a generous gift from Dr. Junji Takeda, Osaka University). The constructed recombinant BAC was then injected into fertilized eggs with a C57BL/6 background. Among the six independent transgenic lines obtained, the Cdhr1tTA was characterized by the OB-specific strongest expression of the internal LacZ gene. The Cdhr1tTA transgene integration was routinely determined with PCR using internal LacZ primers, 5′-CATTGGCGTAAGTGAAGCGAC-3′ and 5′-ATCCCAGCGGTCAAAACAGG-3′, resulting in amplification of 381-bp fragments
The ROSA26-TREmGFP (accession number CDB0300K: http://www2.clst.riken.jp/arg/TG%20mutant%20mice%20list.html) was developed as reported previously (Abe et al., 2011). The tetracycline response element (TRE) cloned from pTRE tight2 (Addgene plasmid # 19407, a generous gift from Dr. Markus Ralser, University of Cambridge) was connected to a GFP gene containing a palmitoylation sequence of GAP43 N terminally, and the construct was flanked by 2xcHS4 insulators (Stewart et al., 2009). To generate a targeting vector, the 2xcHS4 TRE-mGFP with an frt-flanked neomycin resistant gene was inserted into the pMC1-DT-A-ROSA26 vector (Abe et al., 2011). The targeting vector was introduced into the HK3i embryonic stem cells derived from C57BL/6 (Kiyonari et al., 2010). After screening by PCR and Southern blotting, the targeted ES clones obtained were injected into ICR 8-cell stage embryos to generate chimeric mice, after which the chimeric mice were mated with Tg(ACTFLPe)9205Dym transgenic mice (Rodríguez et al., 2000) to remove the neomycin cassette. The offspring were genotyped with PCR using the primers, 5′-GCGATCACATGGTCCTGC-3′ and 5′-AATCGGCCGCTCTAGAAC-3′, for 305-bp fragments from the knock-in mutant allele.
The TaumGFP-nLacZ mice (Hippenmeyer et al., 2005; JAX stock #021162, The Jackson Laboratory) with a CD-1 mixed background were provided by Dr. Silvia Arber (Friedrich Miescher Institute for Biomedical Research), and backcrossed with C57BL/6 wild-type mice for at least four generations before use in this study. The TREtdTomato-sypGFP mice (Li et al., 2010; JAX stock #12345) with the FVB background were obtained from the Jackson laboratory and backcrossed with C57BL/6 mice at least three generations before use in this study.
Histochemistry
Animals were anesthetized and transcardially perfused with 4% paraformaldehyde (PFA)/PBS. The brains were dissected and further fixed with 4% PFA/PBS for 12–24 h, immersed in 30% sucrose/PBS, and frozen in the mixed solution of OCT-compound (Sakura Finetek) and 30% sucrose/PBS at a ratio of 2:1. Coronal sections with thickness depending on mouse age (16 μm for E19.5–P0, 20 μm for P2–P16, 30 μm for P20–P23, unless otherwise indicated) were cut on a cryostat, placed on MAS-coated glass slides (Matsunami Glass), and immunostained as described previously (Mita et al., 2015). The primary antibodies used were as follows: rabbit anti-β-gal (1:5000, Invitrogen catalog #A11132, an old lot); chicken anti-β-gal (1:2000, Abcam catalog #ab9361, RRID: AB_307210); rabbit anti-GFP (1:1000, #598, MBL International catalog #598, RRID: AB_591816); rabbit anti-TBR2 (1:1000, Abcam catalog #ab23345, RRID: AB_778267); rabbit anti-cholecystokinin 26-33 (1:500, Sigma-Aldrich catalog #C2581, RRID: AB_258806); rabbit anti-calretinin (1:5000, Millipore catalog #AB5054, RRID: AB_2068506), rabbit anti-tyrosine hydroxylase (1:1000, Millipore catalog #AB152, RRID: AB_390204); rabbit anti-calbindin D-28K (1:1000, Millipore catalog #AB1778, RRID: AB_2068336); rabbit anti-GFP (1:1000, #598, MBL International catalog #598, RRID: AB_591816); rabbit anti-DARPP-32 (1:1000, Abcam catalog #ab40801, RRID: AB_731843); goat anti-dopamine receptor 1 (1:1000, Frontier Institute catalog #D1R-Go-Af1000, RRID: AB_2571594); rabbit anti-dopamine receptor 2 (1:1000, Frontier Institute catalog #D2R-Rb-Af960, RRID: AB_2571596); and mouse anti-GAD67 (1:300, Millipore catalog #MAB5406, RRID: AB_2278725) antibodies. Chicken instead of rabbit anti-β-gal antibody was used in double staining double staining with with EdU (Fig. 1C) and with rabbit primary antibodies (Fig. 2C–F). In the immunostaining shown in Figures 2C–F, 11D, sections were autoclaved at 105°C for 2 min in an antigen retrieval solution (HistoVT One, Nakalai Tesque) before the primary antibody reactions. The secondary antibodies used were donkey Alexa Fluor 488-conjugated anti-rabbit IgG (1:1000, Life Technologies catalog #A-21206, RRID: AB_141708), donkey Cy3-conjugated anti-rabbit IgG (1:1000, Jackson ImmunoResearch catalog #711-165-152, RRID: AB_2307443), donkey Alexa Fluor 647-conjugated anti-rabbit IgG (1:1000, Jackson ImmunoResearch catalog #711-605-152, RRID: AB_2492288), donkey Cy3-conjugated anti-mouse IgG (1:1000, Jackson ImmunoResearch catalog #715-165-150, RRID: AB_2340813), donkey Alexa Fluor 488-conjugated anti-chicken IgY (1:1000, Jackson ImmunoResearch catalog #703-545-155, RRID: AB_2340375), donkey Alexa Fluor 488-conjugated anti-goat IgG (1:1000, Life Technologies catalog #A-11055, RRID: AB_2534102), and donkey Alexa Fluor 647-conjugated anti-goat IgG (1:1000, Jackson ImmunoResearch catalog #705-605-147, RRID: AB_2340437) antibodies. Some sections were counterstained with DAPI (4',6-diamidino-2-phenylindole; Wako catalog #045-30361) and FluoroMyelin RED (Invitrogen catalog #F34652) according to the manufacturer's instructions.
Extended Data Figure 1-1
Number of neurons counted in Fig.1. Download Figure 1-1, TIF file.
Image acquisition and analysis
The fluorescent images were captured with an inverted confocal microscope (Olympus IX81 FV1000) through 10×, 20×, and 40× objectives (UPlanSAPO) using Fluoview (FV-ASW ver4.02) software and automatically pseudo-colored. Wide images were automatically constructed via multi-area time-lapse with the built-in mosaic-imaging tool. Only the fluorescent images in Figure 1E were taken with a CCD camera (Olympus DP71) on a conventional fluorescent microscope (Zeiss Axioplan2) through a 20× objective lens (FLUAR) using Cell Sens Standard software. The obtained images were rotated and cropped, and the brightness and contrast were non-linearly adjusted equally across of the entire image using Photoshop CS5 software (Adobe).
EdU incorporation assay
A thymidine analog, 5-ethynyl-2'-deoxyuridine (EdU; Tokyo Chemical Industry catalog #E1057) was intraperitoneally injected into a pregnant mouse (50 mg/kg body weight) at the indicated time point either before or after the TM injection that was given at the fixed E12.5. Brains were dissected from the embryos at E19.5 and cryosectioned as explained above. The sections were antigen-retrieved as described above and immunostained with chicken anti-β-gal and the secondary antibody first. Subsequently, the incorporated EdU was detected in 0.1 M Tris (pH 7.6), 2 mM CuSO4, 3 μM Alexa Fluor 555 azide (Thermo Fisher Scientific catalog #A20012), and 10 mM ascorbic acid for 40 min at room temperature. The numbers of neurons labeled for β-gal and those doubly labeled for β-gal and EdU were manually counted in eight independent visual fields (containing 127-307 β-gal-positive neurons) in MOB sections for each mouse with a 40× objective lens (Plan-Apochromat) under a fluorescent microscope (Zeiss Axioplan2) using filter sets 38 HE (ex470/40, em525/50) and 15 (ex549/12, em590). Because the plug-based staging was adopted in this study, the exact time of conception was uncertain, resulting in interlitter variation of development of up to 12 h (Heimer et al., 1987). Therefore, to reduce the litter bias, at least four different litters were used for each quantification; –6 h: 10 mice from four litters; 0 h: 11 mice from five litters; 6 h: 11 mice from four litters; 12 h: 11 mice from four litters; 18 h: 12 mice from four litters; 24 h: 10 mice from four litters.
Quantification of birthdate-tagged OB neurons
P14 or P15 OBs were coronally cut into serial 20 μm-sections. A series of every ten sections were immunostained with rabbit anti-β-gal antibody after the antigen retrieval for TaumGFP-nLacZ reporter (Figs. 1D,G, 2A,B) and with rabbit anti-GFP antibody for ROSA26-TREmGFP reporter (Fig. 3C,D). For the TREtdTomato-sypGFP reporter (Fig. 3A,B), tdTomato-positive cells were counted without antibody staining. MCs were defined as the cells whose cell bodies were positioned in the MC layer (MCL); imTCs in the deep 2/3rd of the EPL; eTCs from the upper 1/3rd of the EPL to the lower 1/3rd of the GLL. The number of the reporter-tagged neuron subtypes within a 500 μm-width of OB layer was manually counted under a fluorescent microscope (Axioplan2. Zeiss) with a 20× objective lens (FLUAR). For each single mouse, 42–74, 37–81, and 37–81 visual fields were counted bilaterally from the entire OBs of TaumGFP-nLacZ (Fig. 1D,G), TREtdTomato-sypGFP (Fig. 3A,B) and ROSA26-TREmGFP (Fig. 3C,D) reporter lines, respectively. The numbers of neuron subtypes were then summed up for each mouse and normalized by the dimension of the examined layers, leading to the value per 500-μm width of layer (Figs. 1D, 3A,C). The proportion of OB neuron subtypes was calculated without normalization by the tissue dimension using the datasets of three mice (Figs. 1G, 3B,D), so that it represents the total cell composition in the whole OB. The exact numbers of neurons counted are shown in Extended Data Figures 1-1, 3-1.
Extended Data Figure 3-1
Number of neurons counted in Fig.3. Download Figure 3-1, TIF file.
For distribution analyses of birthdate-tagged neurons (Fig. 2A,B), the number of β-gal-positive MCs and eTCs were counted from the rostral, dorsal, lateral, medial, and ventral sectors of MOBs. The rostral counts were from coronal sections positioned 200–600 μm from the rostral end of the OB. The other counts were from the most dorsal, lateral, medial and ventral visual fields placed in each OB coronal section at the level of 1000–2800 μm from the rostral end.
The proportion of birthdate-tagged neurons that expressed neuron subtype markers (Fig. 2C–F) was calculated by manually counting the neurons in MOB sections of three P14 mice with a 40× objective lens (Plan-Apochromat) under a fluorescent microscope (Zeiss Axioplan2).
Fluorescent dextran labeling
Whole brains were dissected in Hanks’ balanced salt solution from P14.5 Neurog2CreER(G2A); Cdhr1tTA; TREtdTomato-sypGFP triple heterozygotes that had been given TM at E15.5. A small crystal of tetramethylrhodamine or fluorescein-conjugated fixable dextran (Thermo Fisher catalog #D3308, D3306) was prepared following the described methods (Glover, 1995; Sato et al., 1998) and inserted into the anterior most position of the OT that was accumulated with tdTomato-expressing axons under a fluorescence dissection microscopy (MZFLII, Leica). Until the dextran was transported through the axons, the brains were roll-bottle cultured for 3–4 h at 37°C in Neurobasal medium (Thermo Fisher) containing B27 supplements (Thermo Fisher) in an atmosphere of 95% O2/5% CO2 using a whole embryo culture apparatus (10-310, Ikemoto Rica Kogyo). The brains were then fixed with 4% PFA/PBS and subjected for histochemical analyses as described above. The dextran injections were performed on the both sides of 12 brains, and a confined labeling of the targeted cell cluster in the most anterolateral OT was successful in nine brains (at least one side). In all these successful cases, a small bundle of dextran-labeled axons coursing through the polymorph layer of the OT was observed as shown in the Figure 12C,D.
Virus labeling
AAV2-TRETurboFP635 (3.0 × 1012 vg/ml) and VSVG (Lenti)-TREGFP (4.0 × 1012 vg/ml) were prepared as described previously (Kato et al., 2007, 2018). Neurog2CreER(G2A); Cdhr1tTA double heterozygotes mice (two to five months of age) that had been given TM at E15.5 were anesthetized with medetomidine hydrochloride (0.3 mg/kg) and sodium pentobarbital (20 mg/kg), and head-fixed in an apparatus (EDMS11-054 SR-5M-S, Narishige) under a dissection microscopy. A small area of the bone above the OB was carefully cut out using a microgrinder (UC210, Urawa) with a diamond point (ED101). A pulled borosilicate glass pipette (3B100-75-10, Sutter Instrument) filled with the virus solution was placed in a holder (HI-7, Narishige) on a manual course manipulator (Narishige) and inserted into the medial or lateral side of the MOB according to the atlas of the mouse brain (Paxinos and Franklin, 2008). More specifically, the pipette was pierced through the midline of each OB positioned between 1/3 and 2/3 distance along the anterior-posterior axis, avoiding blood vessels, and penetrated medially for 75–100 μm or laterally for 125–150 μm at ∼60° oblique angle to the surface plane. The virus solution (1 μl) was injected slowly over 10 min using an oil syringe injector (IM-6, Narishige). The window over the OB was then covered with the removed bone and sealed with Vetbond Tissue Adhesive (3M Science #1469). At three to four weeks after the virus injection, the brains were dissected out for histochemical examination. The double virus injections were performed unilaterally to the OB of 13 mice, among which four were successful for focal labeling of both of the viruses and exhibited qualitatively same results as shown in Figure 13.
Statistical analysis
The exact numbers of animals used in individual experiments were described in the figure legends or methods. Numerical values were calculated using Microsoft Excel (Mac 2011, Microsoft). No data were excluded for the quantification. Due to the limitation in the number of mice with the appropriate combination of the recombinant genes, sample sizes of this study were inevitably small, and therefore, all the individual data points were represented in graphs for accuracy (Figs. 1C,D, 2A,B, 3A,C). Because this study does not need any dichotomous judgment, null hypothesis significance testing involving p value was not performed. When the sample size is sufficiently large (Fig. 1C), bootstrap confidence was inferred by the estimation statistics (https://www.estimationstats.com/; Ho et al., 2019).
Availability of data and materials
All data and materials generated in this study are available from the corresponding author on reasonable request. The mouse lines generated in this study are deposited at RIKEN for BDR (http://www2.clst.riken.jp/arg/TG%20mutant%20mice%20list.html) and provided under materials transfer agreements.
Results
Generation of birthdate tag driver mouse
The birthdate tagging takes advantage of neuronal differentiation genes, which are expressed within a short time window immediately after neurons exit their final cell cycle (Fig. 1A). For example, the expression of a basic helix-loop-helix transcription factor, Neurog2, is sharply up-regulated in neuronally committed cells, sustained during the neuronal differentiation phase, and then down-regulated in the maturing phase (Shimojo et al., 2008; Ochiai et al., 2009; Telley et al., 2016). Using enhancers of such neuronal differentiation genes, the TM-inducible Cre recombinase, or CreER gene (Feil et al., 1997) has been expressed transiently in newly born neurons in several genetically modified mice (Fig. 1A; Zirlinger et al., 2002; Battiste et al., 2007; Kim et al., 2011; Winpenny et al., 2011). Although these mice have been primarily used for cell lineage analyses related to the manipulated genes (Zirlinger et al., 2002; Ma and Wang, 2006; Battiste et al., 2007; Kim et al., 2011; Florio et al., 2012), they are theoretically applicable for neuronal birthdating (Sudarov et al., 2011; Toma et al., 2014), because a single administration of TM at a certain developmental stage is expected to induce recombination of loxP sequences only in the CreER-expressing neurons that are in the postmitotic differentiation phase (Fig. 1A).
Because the lineage of OB projection neurons expresses Neurog2 (Schuurmans and Guillemot, 2002), we first tested an existing mouse line in which the endogenous Neurog2 gene is replaced by the CreER gene (Zirlinger et al., 2002). However, even with a sufficient dose of TM, the Cre-loxP recombination efficiency was too low for birthdating analyses of OB neurons. Thus, we decided to generate a more efficient mouse line aimed at neuronal birthdating by a brute-force transgenic approach using a potential enhancer of Neurog2. We generated many transgenic mice and successfully isolated a birthdate tag driver line, Neurog2CreER(G2A) (Fig. 1B; Table 1), that robustly induces loxP recombination in a TM-dependent manner.
Birthdate tagging of OB projection neurons
We first tested whether the Neurog2CreER(G2A) line would work for birthdate labeling of OB neurons using a global neuronal Cre reporter, TaumGFP-nLacZ (Hippenmeyer et al., 2005), which expresses dual nucleus- and membrane-anchored reporters in widespread neurons on loxP recombination (Fig. 1B; Table 1). When the Neurog2CreER(G2A) mouse was crossed with the TaumGFP-nLacZ reporter, a TM injection during gestation stages was successful in labeling different types of OB neurons depending on the injection stage (Fig. 1D–G). More specifically, projection neurons of the AOB were selectively labeled from the very earliest TM injection at E10.5–E11.5 (hereafter called TM10.5-11.5; Fig. 1D). The labeling of MCs in the MOB peaked at TM12.5. Neurons labeled by subsequent injections gradually shifted from deeper to more superficial surfaces of the EPL (Fig. 1E). Judging from the location of their cell bodies (Fig. 1E) and dendrites (Fig. 1F; Mori et al., 1983; Orona et al., 1984), imTCs and eTCs were most effectively labeled at TM13.5 and TM15.5, respectively (Fig. 1D). Figure 1G shows the estimated proportion of OB neuron subtypes tagged at individual TM stages in the whole OB.
To characterize the cell-cycle state of OB neurons that underwent TM-induced recombination, we conducted a double injection of TM at the fixed E12.5 stage and a thymidine-analog EdU at a certain time point either before or after the TM injection (Fig. 1C). The emergence of EdU and reporter double-positive neurons indicated that OB neurons were most susceptible to the TM-induced recombination 6–12 h after final DNA synthesis (Fig. 1C). Taking this time lag into account, the time course of the TM-induced OB labeling (Fig. 1D,G) matches well with the results of previous neuronal birthdating studies using 3H-thymidine in OB neurons (Hinds, 1968; Bayer, 1983; Grafe, 1983).
We next examined the spatial distribution of the birthdate-tagged neurons along the tangential plane of different MOB sectors (Fig. 2A,B). Unlike a previous report describing enrichment of late-born MCs in the ventral MOB (Imamura et al., 2011), the reporter-tagged MCs (Fig. 2A) and eTCs (Fig. 2B) were distributed in a more or less random manner across all MOB sectors regardless of the time of TM injection.
Lastly, neuron marker analyses confirmed that the TM-tagged neurons indeed belong to the MOB projection neuron subtypes. Most TM12.5-, TM13.5-, and TM15.5-tagged neurons expressed the excitatory projection neuron marker TBR2 (Mizuguchi et al., 2012; Stahl et al., 2016; TM12.5: 99.2 ± 0.4%, TM13.5: 99.6 ± 0.3%, TM15.5: 100.0 ± 0.0%, mean ± SEM; Fig. 2C–E). The majority were also positive for cholecystokinin, a widely used but less clear-cut immunohistochemical marker for both MCs and TCs (Seroogy et al., 1985; Cheetham et al., 2015; TM12.5: 64.8 ± 5.7%, TM13.5: 71.6 ± 6.5%, TM15.5: 66.5 ± 5.9%; Fig. 2C–E). None of these TM-tagged neurons expressed interneuron markers (Kosaka et al., 1998; Allen et al., 2007; Fig. 2C–E). The TM17.5-tagged neurons (Fig. 2F) were also positive for TBR2 (100.0 ± 0.0%) but less positive for cholecystokinin (16.1 ± 2.3%); many of the reporter-tagged cells were distributed over a surface area beyond the cholecystokinin-positive cell band at the boundary between the EPL and GLL (Fig. 2F). However, because none expressed the OB interneuron markers (Fig. 2F), it seems generally acceptable to classify these cells as a subset of eTCs based on their cell body positions (Macrides and Schneider, 1982; Hayar et al., 2004b; Antal et al., 2006). Taken together, these results indicate that birthdate tagging indeed subclassifies MOB projection neurons.
Mouse lines to visualize axon trajectories of birthdate-tagged OB neurons
Our next focus was axon trajectories of the birthdate-tagged OB projection neurons completed in adult mice. Because labeling of neurons outside the OB interferes with the observation of labeled axons projected out of the OB, we introduced an OB-specific transgene Cdhr1tTA, by which the tetracycline-controlled transactivator tTA is expressed in an OB neuron-specific manner on excision of the loxP-flanked sequence (Fig. 3; Table 1). The expression of tTA was then reported with amplification using its responsive TRE reporters, TREtdTomato-sypGFP (Li et al., 2010; Fig. 3A) or ROSA26-TREmGFP (Fig. 3C). Thus, three genes including the birthdate tag driver Neurog2CreER(G2A) were combined (Fig. 3A,C) to specifically visualize birthdate-tagged OB axons.
This genetic design contained two caveats. First, the Cdhr1 enhancer is not completely OB-specific (Nagai et al., 2005); a small fraction of neurons in olfactory target areas were occasionally labeled intensely with the reporter depending on TM injection stages. Second, the two TRE reporters marginally differ in labeling of OB neurons even in the same combination with the other two genes (Fig. 3). Specifically, using the TREtdTomato-sypGFP reporter, MCs of the MOB were only faintly labeled and their axons were hardly detected, whereas AOB neurons and eTCs were intensely labeled (Fig. 3A,B). This may be due to cell-type specific epigenetic silencing of the TRE in transgenes as suggested by previous reports (Zhu et al., 2007; Oyer et al., 2009). The ROSA26-TREmGFP reporter, in contrast, nicely visualized MCs of the MOB, but counted fewer numbers of AOB neurons and eTCs (Fig. 3C,D). We assume that these numbers are likely to be underestimated in the quantification, because this membrane-anchored GFP reporter only weakly labels cell bodies while intensely labeling fibers, thereby greatly obscuring cell bodies of AOB neurons and eTCs that are buried in fiber-dense regions. Each reporter thus has its drawback and advantage, but, more importantly, the labeling of OB subtypes was strictly TM-dependent and reproducible in each reporter line. Furthermore, the overall time course of the subtype labeling was consistent with that using the global TaumGFP-nLacZ reporter (Fig. 1D,G). Therefore, by using both reporters and carefully comparing results, we conducted trajectory analyses of birthdate-tagged OB neurons.
Axon trajectories of TM10.5-tagged OB neurons
Using the TREtdTomato-sypGFP reporter, the birthdate tagging at TM10.5 visualized the specific AOB trajectories (Fig. 4). On this early TM injection, the projection neurons in the AOB were selectively labeled (Fig. 4A,A”), and their axons projected almost exclusively to the known AOB target areas (von Campenhausen and Mori, 2000; Mohedano-Moriano et al., 2007). More specifically, the TM10.5 axons clearly coursed through the deep dorsal part of the lateral olfactory tract (LOT; Fig. 4B,C, arrowheads), as reported before (Inaki et al., 2004), and reached all the AOB-specific targets (BAOT, MeA, and PMCo in Fig. 4D–F), but no other areas. When the ROSA26-TREmGFP reporter was instead used for visualization, axons were found to project to both of the AOB and MOB targets, displaying a similar pattern to that of TM11.5 axons with the same reporter that will be described next.
Axon trajectories of TM11.5-tagged OB neurons
With TM11.5 tagging, the TREtdTomato-sypGFP reporter still visualized the AOB-specific projections as those of TM10.5 (Fig. 4). The ROSA26-TREmGFP reporter instead visualized the widespread mixed projections from both the AOB and MOB (Fig. 5). Specifically, in the OB, AOB neurons and MCs of the MOB were both intensely labeled (Fig. 5A). The labeled basal dendrites in the MOB were confined to the deep part of the EPL (Fig. 5A’,A”), confirming the labeling of MCs (Mori et al., 1983; Orona et al., 1984). Their axons covered all the widespread olfactory areas (Fig. 5B–I) including both of the AOB and MOB targets (Schwob and Price, 1984; Hintiryan et al., 2012). There was, however, one exceptional vacant space that excluded these axons in the anterolateral edge of the OT abutting the LOT (Fig. 5D,D’,E,E’, asterisks), which we will focus on later.
Axon trajectories of TM12.5-tagged OB neurons
The TM12.5 tagging using the ROSA26-TREmGFP reporter brought a more specific labeling of MCs of the MOB, leaving AOB neurons mostly unlabeled (Fig. 6A). The basal dendrites in the MOB were more densely labeled, compared with those at TM11.5, but still confined to the deep part of the EPL (Fig. 6A’,A”), supporting the selective labeling of MCs.
At this TM injection stage, labeling of AOB targets was extremely faint (Fig. 6G–I), displaying a complementary image to that of TM10.5 axons (Fig. 4D–F). The TM12.5 axons instead projected to all the known MOB target areas (Fig. 6B–I; Schwob and Price, 1984; Walz et al., 2006; Hintiryan et al., 2012). Their wide diffuse projections were likely due to extensive branching of these axons as shown previously (Ghosh et al., 2011; Sosulski et al., 2011; Igarashi et al., 2012), and indeed recapitulated the general image of divergent central olfactory projections. However, as observed for the TM11.5 axons, we found a small oval space that excluded the TM12.5 axons and remained unlabeled, in the anterolateral edge of the OT abutting the LOT (Fig. 6D,D’,E,E’, asterisks). The wide remaining part of the OT was uniformly and densely filled with TM12.5 axons (Fig. 6E,F).
Axon trajectories of TM13.5-tagged OB neurons
A notable characteristic of TM 13.5 tagging is the labeling of imTCs, although their proportion among all labeled neurons was still below 50% using either of the reporters (Fig. 3B,D). Cell bodies of the labeled imTCs were scattered in varying positions across the EPL (Fig. 7A–A”,I–I”). The labeled basal dendrites, projected from the whole population tagged at this TM stage, were confined to the upper part of the EPL (Fig. 1A’,A”,I’,I”), indicating that imTCs and eTCs are the main sources of the dendrites (Mori et al., 1983; Orona et al., 1984). The internal plexiform layer at the most superficial part of the granule cell layer was also labeled (Fig. 1A,A’,I,I’, arrows), suggesting that at least some of the TM13.5-tagged population have intrabulbar association fibers, bifurcated branches of OB axons coursing through this layer to connect to the other side of the OB (Schoenfeld et al., 1985; Belluscio et al., 2002).
The axon trajectories visualized with the two reporters appeared more or less similar (Fig. 7B–H,J–P). Compared with earlier TM injection stages, the projection area was drastically reduced, and covered only the anterior part of the olfactory target areas, namely all the subdivisions of the AON (Fig. 7B,C,J,K), the lateral half of the OT (Fig. 7D–F,L–N), and the anterior part of the PC (Fig. 7D,E,L,M). Within the LOT, the reporter-labeled axons, bundling together, coursed through the most ventral part down to the caudal end (Fig. 7G,O, arrowheads), at which the axons appeared to simply stall without leaving for any target (Fig. 7H,P). At this TM injection stage, ectopic neurons outside the OB were frequently labeled in olfactory target areas (Fig. 7E,F or L–P). This undesirable labeling interfered with the detailed mapping of OB axons, but the small space in the anterolateral edge of the OT appeared to be weakly accessible to the TM13.5 axons (Fig. 7D,D’,E,E’,L,L’,M,M’, asterisks). At this point, the two reporters seemed to make a slight difference: in the ROSA26-TREmGFP reporter, the axons were dense in the surrounding area (Fig. 7D,D’,E,E’), whereas in the TREtdTomato-sypGFP reporter, the axons more densely projected to the core space in the anterolateral edge of the OT (Fig. 7L,L’,M,M’).
Axon trajectories of TM15.5-tagged OB neurons
On TM15.5 tagging, eTCs at the boundary between the EPL and GLL were prominently labeled in both reporters (Fig. 8A,A’,G,G’). Their cell bodies clung along their basal dendrites that projected horizontally in the upper part of the EPL (Fig. 8A’,A”,G’,G”). Although Individual neurons exhibited substantially different morphologies, e.g., horizontally-oriented cells with only one horizontal basal dendrite or small pyramidal ones with multiple basal dendrites (Fig. 8A”,G”), all these characteristics appeared to belong to a subclass of eTCs (Macrides and Schneider, 1982; Orona et al., 1984; Hayar et al., 2004a; Antal et al., 2006; Takahashi et al., 2016). The internal plexiform layer that contains intrabulbar association fibers was intensely labeled at this TM injection stage (Fig, 8A,A’,G,G”, arrows).
The TM15.5 OB axons showed remarkable convergence onto only two small domains within the wide olfactory target areas, regardless of which reporter was used (Fig. 8B–F,H–L). One is that small space in the anterolateral edge of the OT abutting the LOT (Fig. 8D,D’,E,E’,J,J’,K,K’), which excluded TM11.5 and TM12.5 axons. At closer inspection, the 15.5 axons arose from the entire MOB (Fig. 8A,G), ran for a short distance along the surface of the LOT (Fig. 9A), shifted ventrally from the other LOT axons (Fig. 9B), and finally reached this space (Fig. 9C). Here, the dense axon terminals embraced a small, discrete cluster of cells (Figs. 8D’,J’, dotted circles, 9C, arrowheads). This cell cluster was evident in DAPI staining and identifiable based on its unique location abutting the LOT, but has not been described so far. In anterior coronal sections, this target of TM15 axons appeared as one or two isolated cell clusters medially adjacent to the LOT underneath the arched cell layer of the PC (Figs. 8D,J, 9C). Posteriorly, it was positioned just laterally next to the lateral end of the OT cell layer that bent dorsally like a hook, forming a sidewall of the LOT (Fig. 8E,E’,K,K’), and eventually fused with this hook of the OT cell layer at a more posterior level (Fig. 9D, arrowheads). The wide remaining part of the OT did not seem to receive the TM15.5 axon projections (Fig. 8E,F,K,L).
In addition to this new target cell cluster, the TM15.5 axons projected onto a second target, the pE in the AON (Fig. 8B,H). The AON is a major olfactory target consisting of five subdivisions of unknown function. The four subdivisions other than the pE did not seem to be innervated by the TM15.5 axons (Fig. 8C,I), whereas all five subdivisions seemed to be innervated by TM11.5-13.5 axons (Figs. 5–7); it is however not possible to precisely distinguish the axons that terminate on the pE from those that simply penetrate the pE to terminate on the juxtaposed pars lateralis (Scott et al., 1985).
Axon trajectories of TM17.5-tagged OB neurons
The TM17.5 tagging labeled eTCs that were more superficially positioned than those labeled at TM15.5 (Fig. 10A,A”). Most of these neurons did not have basal dendrites (Fig. 10A”) and resembled to another subset of eTCs described previously (Hayar et al., 2004b; Antal et al., 2006; Takahashi et al., 2016). The labeling of intrabulbar association fibers was not evident at this TM injection stage (Fig. 10A,A’). The TM17.5 axons projected out of the OB only in the anterior part of the LOT (Fig. 10B,C). These axons passed close to the pE (Fig. 10B), and a fraction of them even approached the cell cluster at the anterolateral edge of the OT (Fig. 10D,D’), but did not take off the tract and penetrate into any of the target areas.
Characterization of the target cell cluster of TM15.5 axons
We next characterized the newly identified target of TM15.5 axons in the anterolateral OT. Although typical olfactory targets consist of excitatory principal neurons, the OT is a part of the basal ganglia populated with GABAergic spiny neurons that receive dopamine signals (Millhouse and Heimer, 1984; Murata et al., 2015). The cell cluster targeted by TM15.5 axons strongly expressed dopamine-regulated neuronal phosphoprotein-32 (DARPP), a dopamine-responsive neuron marker (Fig. 11A,E), confirming that it is indeed a component of the OT area. This cell cluster also strongly expressed dopamine receptor D1, but lacked D2 receptor expression (Fig. 11B,D). Dopaminergic axons presumably from the ventral tegmental area (Fallon et al., 1978; Haberly and Price, 1978a; Fallon, 1983; Mooney et al., 1987) projected over this cell cluster as well as the other remaining parts of the OT (Fig. 11C). All these observations suggest that this target cell cluster receives dual inputs, one from the dopaminergic system and the other from the glutamatergic olfactory system through a subset of eTCs, and integrates the information in a similar manner to the striatal direct pathway.
Anatomically, the OT exhibits a chimeric architecture consisting of the cell-layered cortical part, the CAP compartments, and the Calleja islands (Fallon et al., 1978; Meyer and Wahle, 1986; Murata et al., 2015). The D2 receptor-negative characteristic of the target cluster of TM15.5 indicates that it belongs to the CAP compartments (Murata et al., 2015). CAP compartments are defined as about a dozen of periodic grooves toward which the OT cell layer protrudes like ripples approaching and touching the brain surface (Hosoya and Hirata, 1974; Meyer and Wahle, 1986; Murata et al., 2015). Accordingly, each CAP is observed as a brain surface patch that often accompanies a single Calleja island on its bottom. Our observations, taken together, indicate that the target cell cluster of TM15.5 axons represents the most anterolateral isolation of the CAP compartments (aiCAP). Thus, we hereafter refer to this target as aiCAP, although it remains unclear how it was classified in previous studies, due to its atypical isolated configuration. Indeed, at higher magnification, neurons in the aiCAP were simply packed (Fig. 11F) and were not associated with a Calleja island or in direct contact with the brain surface.
To assess where the aiCAP projects, we focally injected fluorescent dextran, an anterograde axon tracer, into this cell cluster (Fig. 12A,B), visually guided by the dense accumulation of TM15.5-tagged axons. The afferent axons of the aiCAP projected posteriorly, creating a small fascicle within the GAD67-expressing plexus in the polymorph layer of the OT (Fig. 12C), the area suggested to be analogous to the globus pallidus in the basal ganglia (Young et al., 1984; Heimer et al., 1987). While the dextran-labeled axons projected posteriorly for 300–600 μm in distance, they densely surrounded neurons that were sparsely embedded in the fascicle (Fig. 12D). These local neurons may be the target of afferent axons of the aiCAP.
Target map representations by TM15.5 axons
The pE of the AON is the only target in which the odor map is topographically represented (Schoenfeld and Macrides, 1984; Yan et al., 2008). We therefore tested whether the TM15.5 axons were responsible for this topographic projection. In mice that had been birthdate tagged at TM15.5, two reporter viruses, each encoding a different TRE-driven fluorescent protein (Fig. 13), were focally injected into the medial and lateral parts of the MOB, respectively. This procedure brightly and specifically labeled eTCs positioned on the medial and the lateral sides in different colors (Fig. 13A). Their axons showed a clear topographic segregation in the pE (Fig. 13B), underlining the contribution of TM15.5 axons to the topographic connection from the MOB to the pE. Such a topography was not observed in the aiCAP, within which axon terminals of TM15.5 axons from different MOB regions intermingled (Fig. 13C). Therefore, this target cell cluster in the anterolateral OT appeared to contain a whole representation of the MOB map.
Development of axon projections to the aiCAP
Lastly, we examined the development of axon projections to the aiCAP. By P4, this target of TM15.5 axon already begun to express DARPP and was identifiable (Fig. 14A). However, the 15.5 axons only penetrated this target from P10 through P14 (Fig. 14B–D). The TM12.5 axons, in contrast, began their projections earlier, penetrating the wide olfactory target areas by P4 (Fig. 14E). Interestingly, however, even at initial projection stages, these axons were clearly excluded from the aiCAP (Fig. 14E–H), which appeared to be saved for later innervation by the TM15.5 axons.
Discussion
The present study investigated axon trajectories of OB neurons born at different times using a newly developed birthdate tag method. Through comparisons of axon projections from differently mixed populations of OB projection neuron subtypes, a new circuit structure emerged from apparently random central olfactory networks; a subset of eTCs represented by the TM15.5 population form a segregated projection only to two small targets, the pE in the AON and the aiCAP in the OT. The pE appeared to be a convergent target of different projection neuron subtypes, whereas the aiCAP appeared to be a selective, specific target of the eTC subset and to contain a complete representation of the MOB odor map. Although it remains unclear how accessible this target is to other subtypes, the majority of MCs represented by TM11.5 and TM12.5 populations simply avoided projecting to the aiCAP, while the other wide olfactory target areas were closely governed by these neurons. These results indicate that the main olfactory system contains at least two segregated projection pathways branching at the glomeruli.
Potential information features represented by the eTC pathway
The branched parallel pathway involving multiple representations of the sensory map is a general organizing principle of sensory systems. Owing to this organization, sensory systems can extract unique features of information from a common input (Young, 1998; Nassi and Callaway, 2009). Since the main olfactory system follows this principle, it can also extract qualitatively different information from the common odor map represented in the MOB.
What kind of features can be represented in the eTC channel? The underdeveloped basal dendrites of eTCs have suggested that they receive weaker lateral inhibitions (Macrides and Schneider, 1982; Orona et al., 1984; Antal et al., 2006; Griff et al., 2008a; Economo et al., 2016). Electrophysiological analyses show that eTCs are activated at lower odor concentrations and an earlier time window than MCs (Schneider and Scott, 1983; Ezeh et al., 1993; Fukunaga et al., 2012; Igarashi et al., 2012). Thus, following odor stimulation, eTCs may play a more active role in early processing of odor information, employing their specific targets. Their functions can be significantly different from those of MCs that are activated later after odor stimulation and convey signals to the other wide target areas that can integrate top-down signals, reflecting various behavioral scenes (Shiotani et al., 2018).
Potential functions of the eTC channel can also be inferred from their specific targets. The neurochemical characteristics of the aiCAP suggest that this pathway belongs to the dopamine reward system (Fallon et al., 1978; Haberly and Price, 1978b; Young et al., 1984; Mooney et al., 1987; Wesson and Wilson, 2011; Gadziola et al., 2015). Considering that the connections develop after the juvenile stage (Fig. 14), the related behavioral output may be based on motivation-related olfactory learning, rather than innate functions, as suggested before (Murata et al., 2015). The pE exclusively projects to the MOB on the contralateral side (Haberly and Price, 1978a; Schoenfeld and Macrides, 1984). Through its topographic connections from the ipsilateral MOB, this subnucleus has been suggested to mediate interactions between the two bilateral glomeruli that receive the same odorant receptor input (Yan et al., 2008; Kikuta et al., 2010).
Fortunately, these entire sets of distinct possibilities are now directly testable with the help of the birthdate tag method. In other words, this study provides the basis for future studies at the functional level that would confirm the anatomic observations reported in this study.
Heterogeneous properties of eTC subtypes
Some previous studies have characterized eTCs as interneurons that do not project out of the OB (Macrides et al., 1985; Schoenfeld et al., 1985; Liu and Shipley, 1994; Nagayama et al., 2014). Based on the present results, we can assume retrospectively that it was extremely difficult to incidentally hit their small targets with a retrograde axon tracer or antidromic electrode stimulation, and that the afferent projections of eTCs may have simply remained undetected in some previous experiments (Haberly and Price, 1977; Skeen and Hall, 1977; Schneider and Scott, 1983).
The present study, at the same time, shows that eTCs are a subdividable population. One subset tagged at TM15.5 contained the projection neurons that send extremely confined afferents to the pE and aiCAP. The E15.5 population had basal dendrites elongating in the upper part of the EPL and dense intrabulbar association fibers (Fig. 8), and more than half were positive for cholecystokinin (Fig. 2E). Another subset tagged at TM17.5 seemed to consist of mostly interneurons, although some had axons projecting into the LOT without any ultimate targets. The TM17.5 population barely had basal dendrites or intrabulbar association fibers (Fig. 10), and the majority were cholecystokinin-negative (Fig. 2F). Interestingly, these characteristics of eTCs subsets based on their birthdates agree well with those of eTC subgroups classified by other criteria (Macrides and Schneider, 1982; Schoenfeld et al., 1985; Liu and Shipley, 1994; Hayar et al., 2004b; Antal et al., 2006; Takahashi et al., 2016).
Uncertainties related to imTCs
Projections of imTCs have remained enigmatic in this study, since we are so far unable to separate them from MCs or eTCs using the birthdate tag approach. The imTCs have been considered representative of TCs, and, accordingly, much of our knowledge about TCs is actually derived from imTCs (Mori et al., 1983; Scott et al., 1985; Ezeh et al., 1993; Nagayama et al., 2004; Griff et al., 2008b; Ghosh et al., 2011). A previous study has visualized complete axon trajectories of single individual imTCs associated with the electrically identified glomeruli (Igarashi et al., 2012). Although only a limited number of cells were analyzed, these imTCs show a characteristic distribution of bushy terminals in the pE and pars posterior of the AON as well as clustered terminals along the border between the OT and the medial edge of the anterior PC. Indeed, these terminal domains resemble those of the TM13.5 population that substantially included imTCs (Fig. 7).
Traditionally imTCs have been considered abundant members in MOB projection neurons (Haberly and Price, 1977; Nagayama et al., 2014), but in the present study, imTCs represented only a minor fraction of the birthdate-tagged OB neurons. Before discussing this seeming discrepancy, however, we should re-estimate the number of imTCs in the mouse MOB. We noticed that the widely-used markers for MOB projection neurons, such as TBR2, cholecystokinin, or Cdhr1, all label only sparse cells in the imTC-containing EPL zone, but much more abundant cells in the two flanking bands that contain MCs and eTCs (Fig. 2; Stahl et al., 2016), mirroring the distribution of birthdate-tagged OB subtypes in this study. Moreover, a recent study shows that a large fraction of neurons in the imTC-containing EPL zone are actually parvalbumin-containing inhibitory interneurons (Kato et al., 2013).
Axon guidance of OB neurons with different birthdates
OB axons project through the LOT in a chronotopic manner (Inaki et al., 2004; Yamatani et al., 2004). Namely, within the tract, early growing axons occupy deeper regions while later growing axons are simply added to the ventral surface side of the preexisting axons. The present study confirmed that this temporal order of axon growth, in fact, reflects the birthdate order of the projection sources; the TM10.5-tagged neurons mainly consisting of AOB neurons projected the axons deep into the LOT (Fig. 4), whereas a more ventral surface space of the LOT were occupied by axons of late-born neurons including TM13.5-17.5 populations (Figs. 7–10).
This chronotopic axon organization within the LOT has been discussed as a potential consequence of physical competition among growing axons over open space on a first-come-first-served basis (Yamatani et al., 2004; Hirata and Iwai, 2019). A different form of guidance seemed to operate for the late-stage targeting of the axon terminals onto the aiCAP. The developmental processes observed here (Fig. 14) suggest that axons from the differently aged neurons do not compete for this target but are selectively guided to their specific targets, from the beginning, by distinct molecular signals.
Neuronal birthdating based on the Neurog2 expression
The traditional birthdating method using thymidine analogs determines the timing of the final S phase in neuronal progenitors (Bayer and Altman, 1987), whereas the present birthdate tag makes use of the expression timing of the Neurog2 gene. Neurog2 is expressed in an oscillatory, fluctuating manner in neural progenitors and in a more sustained manner in neuronally committed cells (Shimojo et al., 2008; Imayoshi and Kageyama, 2014). The latter strong expression acts as the master determinant that directs neural progenitors into the neuronal lineage (Schuurmans and Guillemot, 2002; Imayoshi and Kageyama, 2014), and is likely to be exploited for Cre-loxP recombination in our birthdate tag approach. Indeed, the reported time course of Neurog2 expression in the G1/G0 phase after the final mitosis in cerebellar Purkinje cells (Florio et al., 2012) matches very well with the recombination time course observed in the present study (Fig. 1C).
Neuronal differentiation has been related to the last cell cycle (McConnell, 1989; McConnell and Kaznowski, 1991). More recent studies show that these events can actually be uncoupled experimentally, but are interconnected in vivo at various levels (Hardwick and Philpott, 2014). The key element of the relationship is that up-regulated neuronal differentiation genes such as Neurog2 put an end to the cell cycle, by regulating cell cycle components (Lacomme et al., 2012). There is however a clear deviation from this relationship. Basal progenitors in the cortex (Ochiai et al., 2009; Telley et al., 2016) are neuronally committed cells, but undergo one or two additional rounds of cell cycles to generate solely neurons. They express Neurog2, and in fact, this expression is supposed to be the critical inducer of these fate-restricted progenitors (Miyata et al., 2004; Britz et al., 2006). Thus, Neurog2 expression may be a metric more directly related to the timing of the neuronal commitment than the cell cycle exit.
Recently, new birthdating methods have revolutionized our understanding of neural development (Tabata and Nakajima, 2001; Saito, 2006; Govindan et al., 2018). For example, a new metric based on M phase-dependent cell positioning within the ventricular zone has yielded unexpected results about corticogenesis (Telley et al., 2016; Oberst et al., 2019). Likewise, different metrics of neuronal birthdating may highlight distinct biological consequences under certain situations.
Future perspectives
Recent studies have revealed many more maps in the higher-order brain centers than previously expected, resulting from branched parallel pathways (Young, 1998; Laramée and Boire, 2014; Glasser et al., 2016; Tsukano et al., 2017). We thus wonder whether the detection of an olfactory map may be especially difficult, due to its uniqueness; the map does not represent real anatomic space, but rather repertoires of receptors, whose lack of spatial continuity makes it difficult to detect repetitive individual maps. Using an alternative approach, this study as well as a previous Drosophila study (Jefferis et al., 2001) demonstrate the effectiveness of anatomic dissection of projection neuron subtypes by birthdate-based classification. It is possible that this approach, in the future, will lead to the further subclassification of projection neuron subtypes and reveal additional maps that are buried in the seemingly random olfactory networks. Furthermore, the birthdate tag method will be useful for dissecting functional neuronal networks in other nervous systems as well.
Acknowledgments
Acknowledgements: We thank Yash Hiromi, Shinichi Aizawa, Lena Iwai, and Yan Zhu for comments on this manuscript; Markus Ralser for pTREtight2; Pierre Chambon for CreERT2; Junji Takeda for p23loxZeo; Silvia Arber for TaumGFP-nLacZ mice; and Kazuhiro Ikenaka for the optimized tTA gene.
Footnotes
The authors declare no competing financial interests.
This work was supported by Ministry of Education, Culture, Sports, Science and Technology (MEXT); Japan Society for the Promotion of Science (JSPS); KAKENHI Grants 17H05776, 17H05587, and 16H04659 and by a Research Organization of Information and Systems (ROIS) Challenging Exploratory Research Projects for the Future grant (T.H.).
This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.