Abstract
The medial prefrontal cortex (mPFC) plays a critical role in the control of cognition and emotion. Reciprocal circuits between the mPFC and basolateral amygdala (BLA) are particularly important for emotional control. However, the neurons and synapses that link these brain regions remain largely unknown. Here we examine long-range connections between the mouse mPFC and BLA, using whole-cell recordings, optogenetics, and two-photon microscopy. We first identify two non-overlapping populations of layer 2 pyramidal neurons that directly project to either the BLA or contralateral mPFC. We then show that pyramidal neurons projecting to the BLA receive much stronger excitatory inputs from this same brain region. We next assess the contributions of both presynaptic and postsynaptic mechanisms to this cell-type and input-specific connectivity. We use two-photon mapping to reveal differences in both the synaptic density and subcellular targeting of BLA inputs. Finally, we simulate and experimentally validate how the number, volume, and location of active spines all contribute to preferential synaptic drive. Together, our findings reveal a novel and strong reciprocal circuit that is likely to be important for how the mPFC controls cognition and emotion.
Introduction
The medial prefrontal cortex (mPFC) helps control cognition and emotion (Goldman-Rakic, 1995; Miller, 2000; Fuster, 2001), as highlighted by its dysfunction in neuropsychiatric disorders (Egan and Weinberger, 1997). This ability depends on reciprocal circuits between the mPFC and other cortical and subcortical regions (LeDoux, 2000; Sotres-Bayon et al., 2012). Long-range interactions between the mPFC and basolateral amygdala (BLA) are particularly important for emotional control (Herry et al., 2008; Burgos-Robles et al., 2009). In principle, these interactions could occur through indirect networks involving several brain regions or cortical layers, or via direct communication between neurons in the mPFC and BLA. However, little is known about the projection neurons and synaptic connections that might constitute these reciprocal circuits.
Pyramidal neurons in the mPFC are capable of sending diverse long-range outputs to other brain regions. The targets of layer 5 (L5) pyramidal neurons include the contralateral mPFC (cmPFC), striatum, and pons (Morishima and Kawaguchi, 2006; Dembrow et al., 2010; Morishima et al., 2011; Otsuka and Kawaguchi, 2011). As in other cortical areas, these projection neurons are well known to have distinct properties, including dendritic morphology, intrinsic physiology, and synaptic connectivity (Hattox and Nelson, 2007; Le Bé et al., 2007; Brown and Hestrin, 2009). In the mPFC, L2 pyramidal neurons also send long-range outputs, including to both the BLA and cmPFC (Gabbott et al., 2005; Hirai et al., 2012). These projection neurons thus appear well positioned to help govern the cognitive and emotional functions of the mPFC. However, much less is known about the cellular and synaptic properties of superficial pyramidal neurons compared with those found in deeper layers.
The mPFC also receives long-range excitatory inputs from an assortment of other brain regions (Hoover and Vertes, 2007). These inputs can arrive in layers that are not directly sampled by the dendrites of projection neurons, as in other agranular cortices (Mao et al., 2011; Hooks et al., 2013). To participate in reciprocal circuits, these inputs and neurons must first be linked by interlaminar connections (Weiler et al., 2008; Anderson et al., 2010). In contrast, excitatory inputs from the BLA and cmPFC densely arborize within superficial layers of the mPFC (Sesack et al., 1989; Bacon et al., 1996). These inputs are poised to directly synapse onto L2 pyramidal neurons whose dendrites sample from these layers (Little and Carter, 2012). This unique anatomy could enable bidirectional communication between these pyramidal neurons and the BLA and cmPFC.
Here we examine reciprocal circuits linking L2 pyramidal neurons in the mPFC with both the BLA and cmPFC. We first define two non-overlapping populations of intermingled L2 pyramidal neurons that project to these brain regions. We then establish that neurons projecting to the BLA receive strong excitatory inputs from this same region. We show that this preferential connectivity is attributable to increased numbers of synapses targeting large spines in proximal dendrites. Together, our findings reveal a novel circuit that allows direct communication between the mPFC and BLA. Given the role of long-range interactions in cognitive and emotional control, our findings have important implications for both prefrontal function and disease.
Materials and Methods
Preparation.
We studied L2 pyramidal neurons in acute slices from the prelimbic mPFC of P21 to P28 Swiss Webster mice of either sex. All procedures followed guidelines established by the New York University animal welfare committee. Mice were anesthetized with an intraperitoneal injection of a lethal dose of ketamine/xylazine. After anesthesia, mice were perfused intracardially with an ice-cold solution containing the following (in mm): 65 sucrose, 75 NaCl, 25 NaHCO3, 1.25 NaH2PO4, 25 glucose, 2.5 KCl, 1 CaCl2, 5 MgCl2, 0.4 Na-ascorbate, and 2 Na-pyruvate (bubbled with 95% O2/5% CO2). Coronal sections (300 μm thick) were cut in this solution and transferred to artificial CSF (ACSF) containing the following (in mm): 119 NaCl, 25 NaHCO3, 1.4 NaH2PO4, 25 glucose, 2.5 KCl, 2 CaCl2, 1 MgCl2, 0.4 Na-ascorbate, and 2 Na-pyruvate (bubbled with 95% O2/5% CO2). Slices were recovered for 30 min at 34°C and then stored for at least 30 min at 24°C. All experiments were conducted at 30–34°C, except the two-photon mapping experiments, which were at 22–24°C. For all experiments, 10 μm gabazine and 10 μm d-serine were included to block GABAA receptors and prevent NMDA receptor (NMDAR) desensitization, respectively. For most voltage-clamp experiments, 1 μm TTX and 0.1 mm 4-AP were included to block action potentials (APs) and restore presynaptic glutamate release, respectively. For paired-pulse stimulation experiments, these drugs were omitted to preserve physiological presynaptic APs. For strontium (Sr) experiments, 2 mm Ca in the ACSF was replaced with 2 mm Sr to evoke asynchronous release. In some experiments, 10 μm NBQX or 10 μm 3-((R)-2-carboxypiperazine-4-yl)-propyl-1-phosphonic acid (CPP) were added to block AMPA receptors (AMPARs) and NMDARs, respectively. All chemicals were from Sigma or Tocris Bioscience.
Stereotaxic injections.
To target different neurons and inputs, we performed stereotaxic injections on P9–P12 mice, as described previously (Little and Carter, 2012). Injection site coordinates were relative to bregma (mediolateral axis, dorsoventral axis, and rostrocaudal axis: cmPFC, +0.3 mm, −2.3 to −1.3 mm, and +1.7 mm; BLA, +2.8 mm, −4.6 mm, and −0.7 mm). For anterograde tracing, we injected 100–200 nl of virus [adeno-associated virus 2/9 (AAV2/9)–CAG–mCherry or AAV2/9–CAG–EGFP; University of Pennsylvania Vector Core]. For retrograde labeling, we injected a 0.2% dilution of Alexa Fluor-conjugated cholera toxin subunit B (CTB–Alexa 488 or CTB–Alexa 647; Invitrogen) or undiluted retrobeads (green or red; Lumafluor). For combined retrograde labeling and optogenetics, we injected a 1:1 mixture of 100–200 nl of virus (AAV2/9–CAG–hChR2–mCherry; University of Pennsylvania Vector Core) and either red or green retrobeads. After injections, animals were returned to their cages for 2–3 weeks before being used for experiments.
Electrophysiology.
We targeted L2 pyramidal neurons located in a dense band (approximately two to three neurons or ∼30 μm thick) at the L1/L2 border, ∼200 μm from the pial surface, using infrared-differential interference contrast. Corticoamygdalar (CA) and corticocortical (CC) neurons were identified based on the presence of red or green retrobeads. For voltage-clamp experiments, borosilicate pipettes (2–5 MΩ) were filled with the following (in mm): 135 Cs-gluconate, 10 HEPES, 10 Na-phosphocreatine, 4 Mg2-ATP, 0.4 NaGTP, 0.6 D600, 10 TEA-Cl, 2 QX-314-Cl, and 10 EGTA, pH 7.35 with CsOH (290–295 mOsm). For current-clamp recordings, we used the following (in mm): 135 K-gluconate, 7 KCl, 10 HEPES, 10 Na-phosphocreatine, 4 Mg2-ATP, 0.4 NaGTP, and 0.5 EGTA, pH 7.35 with KOH (290–295 mOsm). For imaging experiments, 40 μm Alexa Fluor 594 and 1 mm Fluo-4FF (Invitrogen) were included in the internal solution. We waited at least 15 min after patching before beginning mapping experiments, to allow for complete dye diffusion throughout the cell. Physiology data were collected with a Multiclamp 700B amplifier. Signals were filtered at 5 kHz for current-clamp recordings and at 2 kHz for voltage-clamp recordings and sampled at 10 kHz. Series resistance in voltage-clamp recordings was <25 MΩ and not compensated.
Optogenetics.
Glutamate release was triggered by activating channelrhodopsin-2 (ChR2) in the presynaptic terminals of different long-range excitatory inputs to the mPFC, as described previously (Little and Carter, 2012). In all optogenetic experiments, nearby CA and CC pairs were recorded in random order within the same slice. For wide-field illumination experiments, a single 0.5–8 ms pulse of 473 nm light from a light-emitting diode (LED) triggered release. When assessing responses in current clamp, data were only included when at least one neuron in the pair fired an AP in response to stimulation. For two-photon mapping and paired-pulse ratio (PPR) experiments, the back focal plane of the objective was filled with collimated light from a 473 nm laser (Laserglow Technologies) using a fiber-optic cable (Thorlabs), planoconvex lenses (Edmund Optics), and a dichroic mirror (Chroma Technologies), as described previously (Chalifoux and Carter, 2011). Beam diameter at the slice surface was ∼30 μm, allowing for focal ChR2 activation. Laser power was ∼4 mW. For two-photon mapping experiments, a fast shutter (Uniblitz) was placed in the detection path to protect the photomultiplier tubes. Laser power and duration were chosen to elicit reliable EPSCs and limit photodamage. For PPR experiments, the laser beam was focused ∼200–300 μm below the soma in L5 to prevent direct depolarization of presynaptic terminals in L1 and L2.
Two-photon microscopy.
Two-photon imaging was performed on a custom microscope, as described previously (Chalifoux and Carter, 2010). Briefly, a titanium:sapphire laser (Coherent) tuned to 810 nm was used to excite Alexa Fluor 594 and Fluo-4FF to image morphology and detect Ca signals, respectively. Ca signals were quantified as the change in Fluo-4FF fluorescence [green (G)] normalized to the Alexa Fluor 594 fluorescence [red (R)], giving units of ΔG/R. These signals were then normalized to the G/R value measured with a saturating concentration of Ca added to the internal solution in a thin-walled pipette, giving final measurements in ΔG/Gsat. Imaging used a 63×, 1.0 numerical aperture (NA) objective (Olympus).
Two-photon mapping.
Active spines were mapped in the dendrites using a variant of a previously described technique (Little and Carter, 2012; MacAskill et al., 2012). As in other experiments, pairs of CA and CC neurons were sequentially patched in random order. After allowing dyes to equilibrate, a systematic search was performed to detect active spines in dendritic segments. Sampling began at dendrites terminating in either L1 or L2 and at segments either proximal or distal to the soma. For each segment, the x, y, and z positions relative to soma were noted for offline analysis of spine locations. A high-resolution two-photon image stack (0.1 × 0.1 × 0.2 μm) was also acquired for analysis of spine morphology. Identical numbers of segments were imaged at similar dendritic locations for both neurons in each pair to ensure unbiased sampling. The synapse density was computed as the number of detected spines in a neuron, divided by the number of sampled segments. At the end of each mapping experiment, a two-photon image stack of the cell was acquired (1.0 × 1.0 × 1.0 or 0.5 × 0.5 × 1.0 μm) for analysis of morphology.
Histology and confocal microscopy.
Mice were anesthetized and perfused intracardially, as described above. Brains were blocked and transferred to 4% paraformaldehyde in 0.01 m PBS for 16–24 h at 4°C before being transferred to 0.01 m PBS. Slices were cut on a VT-1000S vibratome (Leica) at 50 μm thickness and placed on gel-coated glass slides. ProLong Gold antifade reagent with DAPI (Invitrogen) was applied with 0.01 m PBS to the surface of the slices, which were then covered with a glass coverslip. All confocal images were taken on a Leica TCX SP5. Images of anterograde anatomy and retrograde anatomy used for CA and CC localization were taken with a 40×, 1.5 NA oil-immersion objective at 1024 × 1024 pixels (0.34 × 0.34 μm). To ensure that channel bleed-through did not influence our colocalization analysis, fluorophores were chosen with widely separated excitation and emission spectra and were sequentially acquired as dual-fluorescence images.
Simulations.
To predict the probability that CA and CC neurons generate APs in response to BLA input, we modified a previously described model of L2 pyramidal neurons (Little and Carter, 2012). We first sampled the volumes and locations of spines contacted by BLA inputs at CA and CC neurons. We then predicted the size of the somatic unitary EPSP (uEPSP) at each of these spines, based on two-photon uncaging experiments (Little and Carter, 2012). We linearly summed individual uEPSPs to reach an experimentally determined threshold of −48 mV. For each neuron, we simulated 1–100 synapses and computed the spike probability over 1000 runs. To independently assess the impact of spine volume or dendritic location, we also generated predictions in which (1) the CA distribution of spine volume was replaced by the CC distribution or (2) the CA distribution of dendritic location was replaced by the CC distribution. To simulate the effect of the synapse density of BLA inputs, we also ran simulations in which CC neurons received an experimentally determined 50% fewer inputs. Spike probability curves were plotted as a function of the fraction of total number of BLA inputs rather than the absolute numbers because the latter are unknown.
Data analyses.
Imaging and physiology data were acquired using National Instruments boards and custom software written in MATLAB (MathWorks). Online spine Ca signal detection was performed in MATLAB. Image preprocessing was performed in NIH ImageJ and MATLAB. Morphological analysis was conducted in NeuronStudio (Dumitriu et al., 2011). Ca imaging and physiology analysis was performed in Igor Pro (Wavemetrics). In current-clamp experiments, output gain was quantified as the slope of the frequency/current (F/I) curve. In Sr experiments, quantal EPSCs (qEPSCs) were detected in 600 ms windows after the stimulus, using threshold-based, event-detection functions in NeuroMatic, with the threshold set at 4 standard deviations (SDs) above noise. Spontaneous EPSCs (sEPSCs) were detected in a 600 ms window before the stimulus.
For analysis of two-photon image stacks, all visible dendrites of each mapped neuron were reconstructed in NeuronStudio, as described previously (Little and Carter, 2012). To quantify the distributions of axons from our confocal images and determine the colocalization of CA and CC neurons, background subtraction was performed on each slice, followed by manual alignment. For axon analysis, intensity profiles were averaged across multiple slices (for each input, n = 3 animals, 7 slices), as described previously (Little and Carter, 2012). For colocalization, dual-fluorescent pixels were detected and overlaid on a dual-color image of retrogradely labeled CA and CC neurons. Somata with colocalized pixels, as well as large DAPI-positive nuclei, were manually counted in regions 40 × 100 μm centered on L2 of mPFC.
EPSC and qEPSC amplitudes were calculated as the average value over a 1 ms window around the peak. Standard errors (SEs) were calculated via bootstrap resampling (10,000 surrogates). For experiments with sequential recordings, summary comparisons are expressed as both medians ± SE of raw data and medians ± SE of the ratios of these data. Some data are presented in box-plot form, showing the median, interquartile range, and 10–90% range (whiskers). Average EPSCs, qEPSCs, and EPSPs are shown in figures as individual trials or arithmetic mean ± SE. Statistical comparisons were performed using the two-tailed Mann–Whitney rank-sum test, except in Figure 6, in which a bootstrap permutation test on data medians was used. Correlations are reported as the linear correlation coefficient (r), and significance of correlations was tested with a linear correlation test.
Results
Distinct populations of L2 projection neurons
We examined how pyramidal neurons in superficial layers of the prelimbic mPFC communicate with the BLA and cmPFC. We first used anterograde labeling to establish where inputs from these two brain regions are distributed in the mPFC. We injected AAV–EGFP into the BLA and AAV–mCherry into the cmPFC and examined the profiles of fluorescently labeled axons. We found that BLA axons arborized in a narrow band at the border of L1 and L2, whereas cmPFC axons were more distributed throughout the superficial layers (Fig. 1A). These stereotyped profiles suggest that BLA and cmPFC inputs could make unique synaptic connections onto the dendrites of L2 pyramidal neurons in the mPFC.
We next used retrograde labeling to identify which populations of pyramidal neurons project to the BLA and cmPFC. We coinjected fluorescently tagged CTB into the BLA (CTB–Alexa 488) and cmPFC (CTB–Alexa 647) and detected fluorescently labeled corticoamygdalar (CA) and corticocortical (CC) neurons in the mPFC. CA neurons were concentrated near the L1/L2 border, whereas CC neurons were distributed throughout the superficial layers (Fig. 1B). Labeled neurons were a minority of the total L2 population (CA = 9.8 ± 0.7%; CC = 4.9 ± 1.1%; n = 3 animals; Fig. 1C). However, CA neurons constituted the majority of labeled neurons (CA = 62 ± 7%; CC = 31 ± 6%), and these two cell types did not overlap (0.0 ± 0.1%). These results indicate that CA and CC neurons are two distinct populations of pyramidal neurons in the superficial layers of the mPFC.
Morphology and physiology of CA and CC neurons
In deeper layers of the mPFC, different populations of L5 pyramidal neurons have unique dendritic arborizations (Morishima and Kawaguchi, 2006; Dembrow et al., 2010; Morishima et al., 2011; Otsuka and Kawaguchi, 2011). We next sought to establish whether CA and CC neurons also had distinct morphological properties. For these experiments, we injected green retrobeads into the BLA and red retrobeads into the cmPFC, which are taken up by axons and retrogradely transported to the mPFC (Hattox and Nelson, 2007). We then used whole-cell recordings to selectively record from labeled L2 pyramidal neurons (Fig. 2A). We filled CA and CC neurons with fluorescent dyes and used two-photon microscopy to characterize their dendritic arbors (Little and Carter, 2012). We found that the morphology of these neurons was largely similar, as assessed by their total dendritic lengths (CA = 2.8 ± 0.2 mm, n = 19 cells; CC = 3.1 ± 0.3 mm, n = 17 cells; p = 0.79; Fig. 2B). There were only subtle morphological differences, with CC neurons having more dendrites at distal locations (Fig. 2C). These findings suggest that the dendrites of these two neurons may sample excitatory inputs arriving in the superficial layers from the BLA and cmPFC.
Depending on their targets, different populations of deep-layer pyramidal neurons also exhibit a range of firing properties (Morishima and Kawaguchi, 2006; Dembrow et al., 2010; Morishima et al., 2011; Otsuka and Kawaguchi, 2011). We next used current-clamp recordings to characterize the physiological properties of CA and CC neurons. Current steps evoked repetitive APs in both populations (Fig. 2D) and more firing with increasing current (Fig. 2E). Over many recordings, we found that these neurons had indistinguishable resting potentials (CA = −77.7 ± 2.8 mV, n = 9 cells; CC = −80.5 ± 4.3 mV, n = 7 cells; p = 0.72; Fig. 2F), input resistances (CA = 176 ± 16 MΩ; CC = 197 ± 30 MΩ; p = 0.44; Fig. 2G), threshold potentials (CA = −48.3 ± 2.0 mV; CC = −47.9 ± 3.4 mV; p = 0.46; Fig. 2H), and output gains (CA = 2.6 ± 0.34; CC = 2.4 ± 0.32; p = 0.86; Fig. 2I). These results indicate that CA and CC neurons have equivalent physiological properties and may similarly transform synaptic inputs into APs.
Cell- and input-specific synaptic connectivity
Although their anatomical and physiological properties appear similar, CA and CC neurons may participate in distinct circuits. We next used a combination of optogenetics and targeted recordings to assess whether BLA and mPFC inputs make selective contacts onto these two types of pyramidal neurons. To selectively activate inputs from the BLA and cmPFC, we injected AAV–ChR2 into one of these two regions. To identify projection cells, we also injected green retrobeads into the BLA and red retrobeads into the cmPFC of the same animals. After allowing for viral expression and bead transport, we then made sequential whole-cell recordings from pairs of nearby CA and CC neurons while stimulating excitatory inputs with a wide-field LED (Fig. 3A). These experiments were performed in the presence of TTX (1 μm) to block polysynaptic activity and 4-AP (0.1 mm) to restore glutamate release (Petreanu et al., 2009; Cruikshank et al., 2010; Little and Carter, 2012). These targeted recordings of nearby CA and CC neurons enabled paired comparisons of synaptic strength for different inputs and minimized response variability caused by infection efficiency or slice orientation.
In voltage-clamp recordings at −70 mV, with CPP (10 μm) to block NMDARs and gabazine (10 μm) to block GABAARs, we found that both BLA and cmPFC inputs generated AMPAR EPSCs at neighboring CA and CC neurons. However, BLA inputs evoked larger responses at CA neurons (8 ms pulse: CA = 88 ± 20 pA; CC = 31 ± 18 pA; n = 10 pairs; p = 0.007; Fig. 3A). In contrast, cmPFC inputs evoked equivalent responses at the two cell types (CA = 280 ± 133 pA; CC = 288 ± 172 pA; n = 8 pairs; p = 0.25; Fig. 3B). Similar results were observed using a range of light durations, which recruit varying numbers of presynaptic inputs (Fig. 3C,D). On average, BLA inputs generated larger AMPAR EPSCs at CA neurons, whereas cmPFC evoked very similar responses at the two cell types (CC/CA ratio: BLA = 0.25 ± 0.17; cmPFC = 1.07 ± 0.12; p = 0.012; Fig. 3E). We also observed equivalent targeting with NMDAR EPSCs recorded at +40 mV in the presence of NBQX (10 μm) to block AMPARs (CC/CA ratio: BLA = 0.28 ± 0.15, n = 6 pairs; cmPFC = 1.09 ± 0.18, n = 9 pairs; p = 0.036; Fig. 3F). These results indicate that the mPFC participates in direct reciprocal circuits with both the BLA and cmPFC. Moreover, the mPFC establishes a preferentially strong reciprocal circuit with the BLA, because inputs from the BLA show biased synaptic strength onto their associated projection neurons.
Different postsynaptic properties of BLA inputs
In principle, a variety of presynaptic and postsynaptic parameters could contribute to the differential strength of BLA inputs at CA and CC neurons (Gil et al., 1999). One possibility is that the amplitude of quantal responses to BLA inputs differs at these L2 pyramidal neurons. To test this, we substituted extracellular calcium for Sr (2 mm), which enhances asynchronous release and reveals input-specific quantal EPSCs (qEPSCs) (Goda and Stevens, 1994; Hull et al., 2009; MacAskill et al., 2012). Under these conditions, wide-field illumination of BLA inputs generated qEPSCs at both CA and CC neurons (Fig. 4A). However, these responses were significantly larger at CA neurons (CA = 15.7 ± 0.4 pA; CC = 11.0 ± 1.0 pA; CC/CA ratio = 0.75 ± 0.06; n = 7 pairs; p = 0.016; Fig. 4B,C). Moreover, qEPSCs evoked by BLA inputs were significantly larger than spontaneous EPSCs (sEPSCs) before the light pulse at CA neurons (12.6 ± 0.6 pA; p = 0.015) but not CC neurons (12.2 ± 1.6 pA; p = 0.58; Fig. 4B). In contrast, qEPSCs evoked by cmPFC inputs were similar between the two cell types (CA = 10.3 ± 0.4 pA; CC = 10.7 ± 1.9 pA; CC/CA ratio = 1.2 ± 0.2; n = 6 pairs; p = 0.55; Fig. 4C). These results indicate that the individual connections formed by BLA inputs onto CA neurons have enhanced synaptic strength relative to CC neurons. However, although these responses were larger in CA neurons, the CC/CA ratios for qEPSCs were smaller than for compound EPSCs, suggesting that other factors also contribute to differences in overall synaptic strength.
Similar presynaptic properties of BLA inputs
Another possibility is that presynaptic factors, including the probability of release, could differ at BLA inputs onto the two projection cells. We next examined the paired-pulse ratio (PPR) of BLA synapses onto CA and CC neurons, which gives an indication of release probability. For these experiments, we activated BLA inputs with a 473 nm laser spot (∼30 μm diameter) positioned within L5 (∼200–300 μm from L2) (Petreanu et al., 2007). This focal stimulation avoided direct stimulation of presynaptic terminals, which could otherwise disrupt physiological glutamate release (Zhang and Oertner, 2007). In voltage-clamp recordings at −70 mV, we found that paired-pulse stimulation of BLA inputs (interstimulus interval, 50 ms) evoked AMPAR EPSCs at both CA and CC neurons (Fig. 5A). These responses were blocked by TTX (1 μm) and not rescued by subsequent 4-AP (0.1 mm), confirming that they depend on presynaptic activity (n = 6 neurons; Fig. 5B). Consistent with our initial experiments, the first and second EPSCs were larger at CA neurons (EPSC1: CA = 109.4 ± 20.0 pA; CC = 26.0 ± 30.6 pA; EPSC2: CA = 92.4 ± 33.1 pA; CC = 24.6 ± 22.0 pA; n = 13 pairs; Fig. 5C). However, PPR was similar at CA and CC neurons (EPSC2/EPSC1 ratio: CA = 0.98 ± 0.11; CC = 0.90 ± 0.10; p = 0.74; Fig. 5D). Accordingly, there was a high correlation between PPR at CA and CC neurons within each pair (r = 0.69; p = 0.01). These results indicate that the release probability of BLA inputs is similar at these neurons, suggesting that presynaptic factors do not play a major role in producing differences in synaptic strength at these two cell types.
Subcellular synaptic connectivity of BLA inputs
The presynaptic and postsynaptic properties we examined cannot fully explain the large differences in synaptic strength of BLA inputs onto CA and CC neurons. One additional possibility is that BLA inputs also form different numbers of connections onto these two cell types. We next used two-photon microscopy to directly map functional connections at the level of individual spines (Little and Carter, 2012; MacAskill et al., 2012). For these experiments, we sequentially patched neighboring CA and CC neurons and filled them with both a morphological indicator (Alexa Fluor 594) and Ca indicator (Fluo-4FF). We then focused a laser over dendritic segments to activate ChR2 and trigger presynaptic release from BLA inputs. In voltage-clamp recordings at +10 mV, we imaged evoked NMDAR Ca signals to detect spines responding to glutamate release. This approach allowed us to detect functional BLA synapses onto spines located throughout the dendritic arbor of both populations of projection cells.
We first assessed the density of BLA connections onto equivalent dendritic segments in neighboring CA and CC neurons (Fig. 6A). We restricted our two-photon mapping to segments within 150 μm of the L1 border, in which our anatomy indicated that BLA axons are concentrated. We defined synapse density as the number of active spines divided by the number of sampled dendritic segments, with equal sampling of segments across each pair. By mapping multiple segments in each neuron, we determined that synapse density was greater at CA neurons (CA = 0.25 ± 0.04, n = 43 spines; CC = 0.17 ± 0.04, n = 26 spines; CC/CA ratio = 0.50 ± 0.06, n = 21 pairs; p = 0.004; Fig. 6B). These findings indicate that BLA inputs are twice as likely to form connections onto CA neurons.
Additional analysis of our mapping experiments revealed two additional ways in which the subcellular connections made by BLA inputs differ between these neurons. First, we found that the head volumes of BLA-contacted spines were larger at CA neurons (CA = 0.68 ± 0.13 μm3; CC = 0.41 ± 0.13 μm3; p = 0.0052; Fig. 6C). Spines contacted by BLA inputs were also larger than adjacent spines at CA neurons (0.41 ± 0.08 μm3; p = 0.006) but not CC neurons (0.44 ± 0.05 μm3; p = 0.16). These findings are consistent with the results from our Sr experiments and indicate that BLA inputs preferentially contact larger spines on CA neurons but innervate spines on CC neurons more randomly.
Surprisingly, we found that BLA inputs were also more likely to synapse on spines at more proximal locations on CA neurons (CA = 44.3 ± 6.3 μm; CC = 66.9 ± 11.4 μm; p = 0.028; Fig. 6D). Importantly, the sampling distances in these mapping experiments were equivalent at both cell types (CA = 47.3 ± 3.0 μm; CC = 46.2 ± 3.4 μm; p = 0.96). This dendritic targeting was also observed when plotting median spine probability as a function of vertical displacement toward the pial surface (Fig. 6E) or horizontal displacement along the L1/L2 border (Fig. 6F). These findings indicate that BLA inputs target spines near the soma of CA neurons. Together, these two-photon mapping experiments reveal how synapse density and subcellular targeting differ at the two cell types.
Impact of synaptic connectivity on functional output
Our findings revealed synaptic mechanisms that could allow BLA inputs to be stronger at CA neurons than nearby CC neurons. We next determined the relative contributions of these mechanisms to differences in overall synaptic drive at these two populations. We used a linear model of EPSP summation to assess the effects of spine volume, dendritic location, and synapse density on AP firing. We first simulated the effects of subcellular targeting independent of any differences in synapse density (Fig. 7A). We ran the model while sampling the distributions of spine volume and dendritic location measured from CA and CC neurons. We found that fewer BLA inputs were needed to fire CA neurons (53 ± 0.8% of total; n = 20 simulations) than CC neurons (68 ± 0.4%; p = 2 × 10−6). Intermediate values were found when independently swapping either spine volume (61 ± 0.4%; p = 2 × 10−6) or dendritic location (58 ± 0.5%; p = 2 × 10−6) between CC and CA neurons, indicating that both of these parameters contribute to the ability to fire these neurons. Finally, to assess the impact of the synapse density, we scaled the total number of BLA inputs onto CC neurons by 50%, consistent with our mapping experiments (Fig. 7B). We found that reducing synapse density in this way effectively prevented CC neurons from responding by firing APs. Together, these simulations predict that subcellular connectivity enables BLA inputs to form preferentially strong synapses at CA neurons.
To validate these modeling predictions, we next performed current-clamp recordings with wide-field stimulation of excitatory inputs. These experiments were performed in the absence of TTX and 4-AP, which would otherwise perturb synaptic responses and AP firing. To avoid complications from polysynaptic activity, we also restricted our analysis to a short time window after the stimulus pulse (20 ms). We found that BLA inputs reliably triggered APs from resting potentials at CA neurons (8 ms pulse: spike probability = 0.7 ± 0.20; n = 7 pairs) but not nearby CC neurons (0.0 ± 0.20; p = 0.031; Fig. 7C). In contrast, cmPFC inputs reliably triggered APs in both cell types with equal likelihood (CA = 1.0 ± 0.19; CC = 1.0 ± 0.14; p = 0.5; Fig. 7D). Overall, BLA inputs exhibit a strong bias in driving CA neurons (CC/CA ratio = 0.0 ± 0.16), which was not found for cmPFC inputs (CC/CA ratio = 1.0 ± 0.09; Fig. 7E). These results highlight the importance of differences in synaptic strength and show that BLA inputs can preferentially drive feedback interactions in the mPFC.
Discussion
We examined reciprocal circuits that connect the mPFC and two of its key target brain regions. We first identified distinct L2 pyramidal neurons that project to either the BLA (CA neurons) or cmPFC (CC neurons). We then showed that these neurons receive direct excitatory inputs from their own projection targets. Inputs from the BLA are much stronger onto CA neurons than neighboring CC neurons, indicating preferential synaptic connectivity. Using two-photon mapping, we demonstrated that synapse density, spine volume, and dendritic location all differ at CA and CC neurons. These properties combine to enable BLA inputs to drive AP output most effectively in CA neurons. These results highlight the synaptic organization of a direct reciprocal circuit between the mPFC and BLA.
Reciprocal circuits involving L2 pyramidal neurons
In many cortical areas, functional interactions with other brain regions are organized through polysynaptic circuits. In granular cortex, long-range excitatory inputs often arrive in L4, whereas equivalent outputs arise from deeper layers (Douglas and Martin, 2004). However, the mPFC is agranular and lacks a defined L4, meaning that this canonical circuit does not apply. Instead, both deep and superficial layers receive long-range inputs from cortical and subcortical regions (Sesack et al., 1989; McDonald, 1991; Hoover and Vertes, 2007). Pyramidal neurons in different layers also send long-range outputs to a variety of targets (McDonald et al., 1996; Gabbott et al., 2005). Local processing of different inputs undoubtedly occurs within and between individual cortical layers (Kampa et al., 2006; Le Bé et al., 2007; Weiler et al., 2008; Brown and Hestrin, 2009; Anderson et al., 2010; Hooks et al., 2011; Kiritani et al., 2012). However, our anatomy suggested that superficial neurons participate in direct reciprocal circuits with cmPFC and BLA. These neurons are thus uniquely positioned to directly link these distant regions without additional interlaminar signaling.
Many studies have highlighted the diversity of pyramidal neurons in deeper layers of agranular frontal and motor cortices. Intermingled L5 pyramidal neurons in the mPFC that project to distinct cortical and subcortical targets often have different morphological and physiological properties (Morishima and Kawaguchi, 2006; Dembrow et al., 2010; Morishima et al., 2011; Otsuka and Kawaguchi, 2011). We found that L2 pyramidal neurons in the mPFC also have distinct targets, including the BLA and cmPFC, consistent with recent anatomy (Hirai et al., 2012). However, we observed minimal differences in dendritic morphology and intrinsic physiology of CA and CC neurons. These results imply that these neurons may be similarly capable of receiving and responding to long-range excitatory inputs, unlike those in deeper layers. They also suggest that these neurons would be difficult to distinguish based on anatomy or intrinsic physiology alone. Because of these similarities and their scarcity throughout most cortical layers, these neurons may be challenging to isolate with blind in vivo recordings.
Our initial anatomy suggested that both CA and CC neurons could receive long-range excitatory inputs from the BLA and cmPFC. Using optogenetics allowed us to demonstrate functional connections, indicating that both neurons participate in reciprocal circuits. By blocking firing with TTX and restoring release with 4-AP, our recording conditions ensured direct, monosynaptic connections (Petreanu et al., 2009; Cruikshank et al., 2010; Little and Carter, 2012). The amplitudes of light-evoked synaptic responses reflect a variety of experimental factors and are difficult to compare across inputs. However, taking the ratio of the responses at CA and CC neurons normalizes for this variability and allows these comparisons (MacAskill et al., 2012). Using this approach, we found that inputs from the BLA preferentially contact CA neurons over CC neurons. In contrast, inputs from the cmPFC show no preferential synaptic connectivity onto these two populations of neurons. These findings indicate that the type of projection neuron helps define the strength of its long-range connections. These results complement previous findings in other agranular cortices of long-range (Mao et al., 2011; Hooks et al., 2013) and local (Weiler et al., 2008; Anderson et al., 2010; Kiritani et al., 2012) connections onto distinct subpopulations of projection neurons. Together, these studies highlight how synaptic connectivity throughout the cortex is exquisitely organized to enable specific interactions between different cell types.
Mechanisms underlying differences in synaptic strength
One of our key findings was that BLA inputs make preferentially stronger synapses onto CA neurons. In principle, a variety of presynaptic and postsynaptic factors could contribute to stronger synaptic strength, including quantal response amplitude, probability of release, and the number of connections (Gil et al., 1999). Our Sr experiments allowed us to isolate quantal responses (MacAskill et al., 2012) and highlighted both input and cell specificity. The larger quantal responses at CA neurons are likely attributable to the number of postsynaptic AMPARs (Matsuzaki et al., 2001), which are correlated with the size of dendritic spines (Arellano et al., 2007). These differences in quantal responses could also result from differential positioning of synapses in the dendrites (Williams and Mitchell, 2008). These experiments highlight how postsynaptic factors play an important role in determining the strength of long-range inputs.
In contrast, our paired-pulse experiments indicated minimal differences in the probability of release from BLA inputs at different neurons. Importantly, these recordings were made under conditions that preserve presynaptic APs, because direct activation of ChR2 can perturb release from the presynaptic terminal (Zhang and Oertner, 2007; Schoenenberger et al., 2011). One concern is that these monosynaptic EPSCs could be contaminated by polysynaptic activation. However, this is unlikely in our experiments, because pulse durations were kept short (under 2 ms) and the differences in EPSC amplitudes are similar to wide-field stimulation. Although these PPR experiments are only a first step in analyzing presynaptic differences at these synapses, our findings suggest that other factors determine the strength of BLA inputs onto CA and CC neurons.
The observed differences in quantal response amplitude could not fully account for the large disparity in compound response amplitude. An additional factor that could influence synaptic strength is the number of synaptic connections (Cruikshank et al., 2007; Bagnall et al., 2011; Franks et al., 2011). To assess this possibility, we used two-photon mapping to directly measure synapse density at spines in dendritic segments. This technique provides a readout of functional connectivity at the level of individual spines (Little and Carter, 2012; MacAskill et al., 2012). It is conceptually similar to mapping experiments used to identify functional inputs onto neurons in vivo (Chen et al., 2011). Our two-photon mapping experiments revealed that BLA inputs are twice as likely to make connections onto CA neurons. These experiments also showed differences in subcellular targeting, with BLA inputs contacting larger and more proximal spines at CA neurons than nearby CC neurons. These results demonstrate how multiple synaptic factors distinguish BLA connections onto different projection neurons.
In other neurons, synaptic strength is influenced by dendritic location (Williams and Stuart, 2002), spine volume (Matsuzaki et al., 2001; Arellano et al., 2007), or their combination (Richardson et al., 2009; MacAskill et al., 2012). These subcellular properties can also strongly influence the strength of excitatory inputs onto L2 pyramidal neurons (Little and Carter, 2012). Using a model that incorporated our previous two-photon uncaging results, we found that spine and dendrite targeting causes BLA inputs to be stronger at CA neurons. However, synapse density also proved critical in determining overall synaptic strength and drive at these neurons. Thus, a combination of factors allows a common afferent to differentially excite nearby cortical projection neurons that are otherwise similar in their dendritic anatomy and intrinsic physiology. In the future, it will be interesting to determine how these factors are changed in fear conditioning and other emotional behaviors that can alter interactions between mPFC and BLA.
Functional roles of synaptic connections between mPFC and BLA
The mPFC governs a range of higher-level behaviors, including working memory and attention (Goldman-Rakic, 1995; Miller, 2000; Fuster, 2001). One of the best-characterized functions of the mPFC is the control of emotional behaviors, which are thought to critically depend on long-range interactions with the BLA (LeDoux, 2000; Maren and Quirk, 2004). Activity within the prelimbic mPFC correlates with the expression of conditioned fear (Burgos-Robles et al., 2009). Moreover, pharmacological inactivation of the prelimbic mPFC prevents the expression of this behavior (Corcoran and Quirk, 2007). Distinct populations of neurons in the BLA project to the mPFC and play unique functional roles during the learning and extinction of conditioned fear (Herry et al., 2008). Moreover, activating the BLA may enable sustained firing of mPFC neurons during fear conditioning (Sotres-Bayon and Quirk, 2010; Sotres-Bayon et al., 2012). Although indirect networks could mediate these interactions, our findings highlight the presence of direct connections between the mPFC and BLA. Equivalent closed-loop neural architectures are proposed to enhance the gain of other local and long-range circuits (Douglas et al., 1995; Mongillo et al., 2008). In the future, it will be interesting to examine how these long-range interactions are influenced by local connections within and between different layers of the mPFC.
Footnotes
This work was supported by National Institutes of Health Grant R01-MH085974 (A.G.C.). We thank members of the Carter laboratory for helpful discussions and comments on this manuscript.
The authors declare no competing financial interests.
- Correspondence should be addressed to Adam G. Carter, Center for Neural Science, New York University, 4 Washington Place, New York, NY 10003. adam.carter{at}nyu.edu