Extended field-of-view ultrathin microendoscopes for high-resolution two-photon imaging with minimal invasiveness

  1. Andrea Antonini
  2. Andrea Sattin
  3. Monica Moroni
  4. Serena Bovetti
  5. Claudio Moretti
  6. Francesca Succol
  7. Angelo Forli
  8. Dania Vecchia
  9. Vijayakumar P Rajamanickam
  10. Andrea Bertoncini
  11. Stefano Panzeri
  12. Carlo Liberale
  13. Tommaso Fellin  Is a corresponding author
  1. Optical Approaches to Brain Function Laboratory, Istituto Italiano di Tecnologia, Italy
  2. Nanostructures Department, Istituto Italiano di Tecnologia, Italy
  3. University of Genova, Italy
  4. Neural Coding Laboratory, Istituto Italiano di Tecnologia, Italy
  5. Neural Computation Laboratory, Center for Neuroscience and Cognitive Systems @UniTn, Istituto Italiano di Tecnologia, Italy
  6. Center for Mind and Brain Sciences (CIMeC), University of Trento, Italy
  7. Biological and Environmental Science and Engineering Division (BESE), King Abdullah University of Science and Technology (KAUST), Saudi Arabia

Abstract

Imaging neuronal activity with high and homogeneous spatial resolution across the field-of-view (FOV) and limited invasiveness in deep brain regions is fundamental for the progress of neuroscience, yet is a major technical challenge. We achieved this goal by correcting optical aberrations in gradient index lens-based ultrathin (≤500 µm) microendoscopes using aspheric microlenses generated through 3D-microprinting. Corrected microendoscopes had extended FOV (eFOV) with homogeneous spatial resolution for two-photon fluorescence imaging and required no modification of the optical set-up. Synthetic calcium imaging data showed that, compared to uncorrected endoscopes, eFOV-microendoscopes led to improved signal-to-noise ratio and more precise evaluation of correlated neuronal activity. We experimentally validated these predictions in awake head-fixed mice. Moreover, using eFOV-microendoscopes we demonstrated cell-specific encoding of behavioral state-dependent information in distributed functional subnetworks in a primary somatosensory thalamic nucleus. eFOV-microendoscopes are, therefore, small-cross-section ready-to-use tools for deep two-photon functional imaging with unprecedentedly high and homogeneous spatial resolution.

Introduction

The amount of information carried by neural ensembles and the impact that ensemble activity has on signal propagation across the nervous system and on behavior critically depend on both the information and tuning properties of each individual neuron and on the structure of correlated activity, either at the level of correlations between each pair of neurons or at the whole network level (Ni et al., 2018; Salinas and Sejnowski, 2001; Shahidi et al., 2019; Shamir and Sompolinsky, 2006; Panzeri et al., 1999). To study neuronal population coding, it is thus essential to be able to measure with high-precision, large signal-to-noise-ratio (SNR), and without distortions the activity of individual neurons and the relationship between them (Aharoni et al., 2019). Two-photon imaging makes it possible to record the activity of many hundreds of individual neurons simultaneously and provides reliable measures of correlated neuronal events (Yang et al., 2016; Kazemipour et al., 2019; Villette et al., 2019; Runyan et al., 2017; Rumyantsev et al., 2020). Light scattering within the brain, however, strongly affects the propagation of excitation and emission photons, making effective imaging increasingly difficult with tissue depth (Ji et al., 2010; Wang et al., 2014; Helmchen and Denk, 2005). Various strategies have been developed to improve imaging depth in multi-photon fluorescence microscopy (Kobat et al., 2009; Kobat et al., 2011; Theer et al., 2003; Horton et al., 2013; Ji et al., 2012; Mittmann et al., 2011; Tischbirek et al., 2015; Birkner et al., 2017), allowing the visualization of regions 1–1.6 mm below the brain surface. However, deeper imaging requires the use of implantable microendoscopic probes or chronic windows, which allow optical investigation of neural circuits in brain regions that would otherwise remain inaccessible (Jung and Schnitzer, 2003; Flusberg et al., 2008; Barretto et al., 2009; Resendez et al., 2016; Bocarsly et al., 2015).

A critical barrier to progress is the lack of availability of microendoscopic devices with small cross-sections that maintain cellular resolution across a large FOV, to allow high-resolution and high SNR two-photon population imaging on a large number of neurons while minimizing tissue damage. Preserving short- and long-range connectivity is fundamental to study physiological network dynamics in brain circuits. For example, severing thalamocortical fibers alters low-frequency spontaneous oscillations in the thalamocortical loop (Lemieux et al., 2014; Lemieux et al., 2015). Current microendoscopes for deep imaging are frequently based on gradient index (GRIN) rod lenses, which typically have diameter between 0.35–1.5 mm and are characterized by intrinsic optical aberrations (Barretto et al., 2009). These aberrations are detrimental in two-photon imaging because they decrease the spatial resolution and lower the excitation efficiency, especially in the peripheral parts of the FOV. This leads to uneven SNR and spatial resolution across the FOV and therefore restricted effective FOV (Wang and Ji, 2013; Ji et al., 2016). This is a serious issue when measuring correlated neuronal activity, because signals sampled at the periphery of the FOV will be more contaminated by neuropil or by activity of neighboring cells compared to signals sampled near the optical axis. Optical aberrations in GRIN microendoscopes can be corrected with adaptive optics which, however, requires significant modification of the optical path (Wang and Ji, 2013; Wang and Ji, 2012; Bortoletto et al., 2011) and may limit the temporal resolution of functional imaging over large FOVs (Wang and Ji, 2013). Alternatively, the combination of GRIN lenses of specific design with plano-convex lenses within the same microendoscopic probe has been used to increase the Numerical Aperture (NA) and to correct for aberrations on the optical axis (Barretto et al., 2009). However, technical limitations in manufacturing high-precision free-form optics with small lateral dimensions have so far prevented improvements in the performances of GRIN microendoscopes with lateral diameter <1 mm using corrective optical microelements (Matz et al., 2015).

Here, we report the design, development, and characterization of a new approach to correct aberrations and extend the FOV in ultrathin GRIN-based endoscopes using aspheric lenses microfabricated with 3D micro-printing based on two-photon lithography (TPL) (Liberale et al., 2010). These new endoscopic probes are ready-to-use and they require no modification of the optical set-up. We applied eFOV-microendoscopes to study the encoding of behavioral state-dependent information in a primary somatosensory thalamic nucleus of awake mice in combination with genetically-encoded calcium indicators.

Results

Optical simulation of eFOV-microendoscopes

Four types (type I-IV) of eFOV-microendoscopes of various length and cross section were designed, all composed of a GRIN rod, a glass coverslip and a microfabricated corrective aspheric lens (Figure 1). One end of the GRIN rod was directly in contact with one side of the glass coverslip while the corrective aspheric lens was placed on the other side of the coverslip. GRIN rods and corrective lenses were different in each of the four types of eFOV-microendoscopes (lateral diameter, 0.35–0.5 mm; length, 1.1–4.1 mm; all 0.5 NA). The glass coverslip was 100 µm thick for type I, III, IV eFOV-microendoscopes and 200 µm thick for type II eFOV-microendoscopes. This design did not require additional cannulas or clamping devices (Kim et al., 2012; Jung et al., 2004) that would increase the lateral size of the microendoscope assembly or reduce its usable length, respectively.

Optical design of eFOV-microendoscopes.

(a-d) Ray-trace simulations for the four different eFOV-microendoscopes (type I-IV). The insets show the profiles of corrective polymeric lenses used in the different eFOV-microendoscopes. For each eFOV-microendoscope, it is specified the thickness of the coverslip, the length, the diameter of the GRIN rod, the pitch of the GRIN rod, and the working distance, in air or in tissue, at which the simulated correction of aberrations was performed. See also Supplementary file 1 - Table 1.

For each type of GRIN rod used in the eFOV-microendoscopes, ray-trace simulations were performed at λ = 920 nm and determined the profile of the aspheric lens (Figure 1) that corrected optical aberrations and maximized the FOV (Figure 2). In the representative case of type I eFOV-microendoscopes, the simulated corrective lens had a diameter of 0.5 mm, height <40 µm, and the coefficients used in Equation (1) (see Materials and methods) to define the lens profile are reported in Supplementary file 1 - Table 1. For this type of eFOV-microendoscope, the simulated point-spread-function (PSF) at incremental radial distances from the optical axis (up to 200 μm) showed that the Strehl ratio of the system was >80% (i.e. a diffraction-limited condition is achieved according to the Maréchal criterion [Smith, 2008]) at a distance up to ~165 µm from the optical axis with the corrective lens, while only up to ~70 µm for the same optical system without the corrective lens (Figure 2a). This improvement led to a ~5 times increase in the area of the diffraction-limited FOV. Figure 2b–d reports the Strehl ratio for simulated uncorrected and corrected type II-IV eFOV-microendoscopes. Simulations showed that the area of the FOV was ~2–9 times larger for these other types of eFOV-microendoscopes, compared to microendoscopes without the corrective lens. We evaluated the simulated excitation PSF for the four types of microendoscopes (Figure 2—figure supplements 12). We found that the axial dimension of the PSF in lateral portions of the FOV remained smaller and more similar to the axial dimension of the PSF in the center of the FOV in corrected compared to uncorrected endoscopes.

Figure 2 with 4 supplements see all
Corrective lenses improve the simulated optical performance of ultrathin microendoscopes.

(a) Simulated diffraction PSFs to assess the Strehl ratio of the designed microendoscope (type I microendoscopes) without the corrective lens (uncorrected, left) and with the corrective lens (corrected, right). PSFs are shown color-coded according to their radial distance from the optical axis. The black dotted line represents the diffraction-limited condition, which was set at 80% (Maréchal criterion). Distances written in red indicate the radial positions at which the maximal normalized irradiance of the corresponding PSF was >80%. (b–d) Same as in (a) for type II (b), type III (c), and type IV (d) microendoscopes. See also Figure 2—figure supplements 14.

The simulated focal length in the absence and presence of the corrective microlens for the four different types of microendoscopes is reported in Supplementary file 1 - Table 2. All simulations reported above were performed maximizing aberration correction only in the focal plane of the microendoscopes. We explored the effect of aberration correction outside the focal plane. In corrected endoscopes, we found that the Strehl ratio was >0.8 (Maréchal criterion) in a 1.7–3.7 times larger volume compared to uncorrected endoscopes (Figure 2—figure supplement 3). We then investigated the effect of changing wavelength on the Strehl ratio. We found that the Strehl ratio remained >0.8 within at least ±15 nm from λ = 920 nm (Figure 2—figure supplement 4), which covered the limited bandwidth of our femtosecond laser.

Fabrication of eFOV-microendoscopes

Corrective lenses were experimentally fabricated using TPL (Liberale et al., 2010; Figure 3—figure supplement 1a,b) and plastic molding replication (Schaap and Bellouard, 2013) directly onto the glass coverslip (see Materials and methods). Experiments and optical characterization were performed using lens replica only. Fabricated lenses had profile largely overlapping with the simulated one (Figure 3—figure supplement 1c). The corrective lens was aligned to the appropriate GRIN rod using a customized optomechanical set-up (Figure 3—figure supplement 2a,b) to generate eFOV-microendoscopes. To experimentally validate the optical performance of fabricated eFOV-microendoscopes, we first coupled them with a standard two-photon laser-scanning system using a customized mount (Figure 3a,b and Figure 3—figure supplement 2c,d). We initially measured the on-axis spatial resolution by imaging subresolved fluorescence beads (diameter: 100 nm) at 920 nm. We found that eFOV-microendoscopes had similar on-axis axial resolution compared to uncorrected probes (Figure 3c,e,g,i left and Supplementary file 1 - Table 3). Given that most aberrations contribute to decrease optical performance off-axis, we repeated the measurements described above and measured the axial and lateral resolution at different radial distances from the center of the FOV. As expected by ray-trace simulations (Figure 2), eFOV-microendoscopes displayed higher and more homogeneous axial resolution in a larger portion of the FOV (Figure 3c,e,g,i left and Figure 3—figure supplement 3). Defining the effective FOV as that in which FWHMz <10 µm (dotted line in Figure 3c,e,g,i left), we found ~3.2–9.4 folds larger effective FOV in corrected microendoscopes compared to uncorrected probes (Supplementary file 1 - Table 3).

Figure 3 with 7 supplements see all
Optical characterization shows enlarged effective FOV in corrected ultrathin microendoscopes.

(a) Schematic of the eFOV-microendoscope mount for head implant. The GRIN rod is glued to one side of the glass coverslip, the microfabricated polymer lens to the other side of the coverslip. The coverslip is glued on a circular metal ring that facilitates fixation on the animal's skull. (b) Zoom in of the portion highlighted with the red dotted line in (a). (c) Axial (left) and lateral (right) spatial resolution (see Materials and methods for definition) evaluated using subresolved fluorescent beads (diameter: 100 nm) as a function of the radial distance from the center of the FOV for type I uncorrected (dashed line) and corrected (solid line) microendoscopes. Points represent values obtained by averaging at least eight measurements from three different probes (see Supplementary file 1 - Table 1), while error bar represents standard deviation (sd). Lines are quartic functions fitting the data. The red dashed line indicates a threshold value (10 µm) to define the limit of the effective FOV. (d) x,z projections (x, horizontal direction; z, vertical direction) of a z-stack of two-photon laser-scanning images of a subresolved fluorescent layer (thickness: 300 nm) obtained using a type I eFOV-microendoscope, without (uncorrected, top) and with (corrected, bottom) the microfabricated corrective lens. λexc = 920 nm. (e,f) Same as in (c,d) for type II eFOV-microendoscopes. (g,h) Same as in (c,d) for type III eFOV-microendoscopes. (i,j) Same as in (c,d) for type IV eFOV-microendoscopes. See also Figure 3—figure supplements 17 and Supplementary file 1 - Table 3.

Figure 3—source data 1

Experimental values of the axial and lateral spatial resolution as a function of radial distance from the center of the FOV in uncorrected and corrected microendoscopes.

https://cdn.elifesciences.org/articles/58882/elife-58882-fig3-data1-v3.xlsx

To visualize the profile of fluorescence intensity across the whole diameter of the FOV for both uncorrected and corrected probes, we used a subresolved thin fluorescent layer (thickness: 300 nm) as detailed in Antonini et al., 2014. Figure 3d,f,h,j shows the x, z projections of the z-stack of the subresolved fluorescent layer for uncorrected and corrected type I-IV microendoscopes. In agreement with the measurements of spatial resolution using subresolved fluorescent beads (Figure 3c,e,g,i), eFOV-microendoscopes displayed higher intensity and smaller FWHMz in peripheral portions of the FOV compared to uncorrected probes (Figure 3d,f,h,j). eFOV-microendoscopes were characterized by a curved FOV and this distortion was evaluated using a fluorescent ruler (Figure 3—figure supplement 4) and corrected for in all measurements of spatial resolution (Figure 3 and Supplementary file 1 - Table 3). The ability of eFOV-microendoscopes to image effectively larger FOV compared to uncorrected probes was also confirmed in biological tissue by imaging neurons expressing the green fluorescence protein (GFP) in fixed brain slices (Figure 3—figure supplement 5).

Validation of eFOV-microendoscopes for functional imaging in subcortical regions

To validate eFOV-microendoscopes performance for functional measurements in vivo, we first expressed the genetically-encoded calcium indicator GCaMP6s in the mouse hippocampal region in anesthetized mice (Figure 3—figure supplement 6) and in the ventral posteromedial nucleus of the thalamus (VPM), a primary sensory thalamic nucleus, in awake head-restrained mice (Figure 3—figure supplement 7). To this aim, we injected either adenoassociated viruses (AAVs) carrying a flex GCaMP6s construct together with AAVs carrying a Cre-recombinase construct under the CamKII promoter in the hippocampus of wild type mice (Figure 3—figure supplement 6a-a2) or AAV carrying a floxed GCaMP6s in the VPM (Figure 3—figure supplement 7a,b) of Scnn1a-Cre mice (Madisen et al., 2010), a mouse line which expresses Cre in a large subpopulation of VPM neurons. These experimental protocols resulted in the labeling of CA1 hippocampal and VPM neurons, respectively (Figure 3—figure supplements 67). We imaged spontaneous activities in the CA1 hippocampal region with type I eFOV-microendoscopes (Figure 3—figure supplement 6b-c1) and type III eFOV-microendoscopes (Figure 3—figure supplement 6d-e1) or spontaneous activities in the VPM with the longer type II eFOV-microendoscopes (Figure 3—figure supplement 7b–d). Tens to hundreds of active ROIs per single FOV were identified and could be imaged using eFOV-microendoscopes (Figure 3—figure supplements 67).

Higher SNR and more precise evaluation of pairwise correlation in eFOV-microendoscopes

To establish a quantitative relationship between the improved optical properties of eFOV-microendoscopes and their potential advantages for precisely detecting neuronal activity, we generated two-photon imaging t-series using synthetic GCaMP data. This approach allowed us to compare results of the simulation of calcium data for both uncorrected and corrected endoscopic probes with the known ground truth of neuronal activity. Simulated neuronal activity within a volumetric distribution of cells was generated according to known anatomical and functional parameters of the imaged region (the VPM in this case) and established biophysical properties of the indicator (Chen et al., 2013; Dana et al., 2019) (see Materials and methods). t-series were generated by sampling simulated neuronal activity across an imaging focal surface resembling the experimental data obtained using the representative case of type II GRIN lenses for both eFOV-microendoscopes and uncorrected probes (Figure 4a,b). To scan the imaging focal surface, we used an excitation volume which resembled the aberrated and aberration-corrected PSFs experimentally measured for uncorrected and eFOV-microendoscopes, respectively (Figure 3, see Materials and methods). Fluorescence traces were extracted from artificial t-series and compared between the uncorrected and corrected case. On average, a slightly larger number of ROIs could be identified in corrected probes (Figure 4c). Crucially, we observed a nonlinear interaction between the type of probe, the number of detected ROIs, and the SNR of the calcium traces (Figure 4d). Using corrected probes did not always allow the identification of a higher number of ROIs (p=0.096, two-way ANOVA with respect to probe type). Rather, the use of corrected probes allowed to segment more ROIs with high SNR, shifting the distribution of SNR across ROIs to higher mean SNR values (Figure 4d, p=1E-5 for ANOVA with respect to the interaction, and Figure 4e). Corrected endoscopes allowed to segment a smaller number of ROIs with low SNR likely because: (i) in the segmentation method we implemented, which was based on the ground truth distribution of the neurons in the simulated sample, at least two pixels belonging to a ground truth neuron were defined as a ROI; (ii) in uncorrected microendoscopes, the axial PSF largely increased as a function of the radial distance. The enlarged axial PSF in the lateral portions of the FOV augmented the probability of sampling voxels belonging to multiple neurons located at different z positions. Once projected in the 2D plane, the contribution of multiple neurons located at different z positions increased the probability of having pixels belonging to ROIs. An increased axial PSF thus led to an increased number of detected ROIs; (iii) corrected endoscopes had smaller axial PSF compared to uncorrected ones in lateral portions of the FOV potentially leading to smaller number of detected ROIs.

Figure 4 with 2 supplements see all
eFOV-microendoscopes allow higher SNR and more accurate evaluation of pairwise correlation.

(a) Schematic of the procedure for in silico simulation of imaging data. Neuronal activity was simulated within spheres located in a 3D volume, integrated over an elliptical PSF (blue) that was scanned on a curved FOV, projected on a 2D plane to generate artificial frames, and modulated through an intensity mask. Only voxels falling within the PSF contributed to the pixel signal (black trace). (b) Segmentation of in silico data for uncorrected (left, red lines indicate identified ROIs) and corrected (right, blue lines indicate identified ROIs) endoscopes. (c) Number of ROIs segmented in simulated FOVs for uncorrected (Unc.) and corrected (Cor.) microendoscopes. n = 54 segmentations from nine simulated FOVs, Wilcoxon signed-rank test, p=0.037. In this as well as in other figures, values from individual experiments are shown in gray, the average of all experiments in black and error bars indicate sem, unless otherwise stated. In this as well as in other figures: *p<0.05; **p<0.01; ***p<0.001; NS, not significant. (d) Number of segmented ROIs as a function of the peak-SNR threshold in artificial data from n = 9 simulated experiments. A two-way ANOVA with interactions showed a significant effect of peak-SNR threshold (p=1E-50) and of the interaction between peak-SNR threshold and probe type (p=1E-5) on the number of segmented ROIs, while the effect of probe type was not significant (p=0.096). (e) Average SNR of calcium signals under the different experimental conditions (peak-SNR threshold = 20 for the segmentation). n = 987 and 1603 ROIs for nine simulated experiments with uncorrected and corrected microendoscopes, respectively. Mann-Whitney test, p=5E-52. (f) Schematic representation of how radial distance (distr) and pairwise distance (distp) were calculated. (g), Pairwise correlation as a function of radial distance for simulated experiments with uncorrected (left) and corrected (right) microendoscopes. In this as well as other figures, lines show the linear regression of data. Shaded areas represent 95% confidence intervals. n = 738 and n = 869 pairwise correlations for uncorrected and corrected microendoscopes, respectively, from the n = 9 simulated experiments shown in (e). Linear regression fit of data: slope = 0.002, permutation test, p=0 and slope = −2E-4, permutation test, p=0.27 for uncorrected and corrected microendoscopes, respectively. (h) Distribution of the correlation of calcium signals with the first (left) or second (right) component of the ground truth for experiments with uncorrected and corrected microendoscopes. First component: n = 987 and 1603 ROIs for uncorrected and corrected microendoscopes, respectively, from n = 9 simulated experiments. Second component: n = 62 and 85 ROIs for uncorrected and corrected microendoscopes, respectively, from n = 9 experiments. (i) Schematic of the experimental configuration in awake animals. (j) Two-photon images of VPM neurons expressing GCaMP7f showing manually identified ROIs for uncorrected (left, red lines) and corrected (right, blue lines) type II microendoscopes. (k) Spatial density of ROIs identified in in vivo experiments. n = 9 FOVs and 6 FOVs from three animals with uncorrected and corrected microendoscopes, respectively. Student’s t-test, p=0.92. (l) SNR of segmented ROIs in in vivo recordings. n = 557 ROIs from 9 FOVs for uncorrected microendoscopes; n = 306 from 6 FOVs for corrected microendoscopes. Mann-Whitney test, p=0.0011. (m) Pairwise correlation as a function of radial distance for in vivo experiments. Number of pairwise correlations: n = 168 from 9 FOVs, n = 36 from 6 FOVs, and n = 92 from 24 FOVs for uncorrected, corrected (dataset 1), and corrected (dataset 2), respectively. Dataset two was obtained from experimental sessions performed during behavioral state monitoring as in Figure 6. Linear regression fit of data: slope = 0.0002, permutation test p=0.006 for uncorrected; slope = −0.0005, permutation test p=0.05 for dataset 1; slope = −0.0007, permutation test p=0.05 for dataset 2. See also Figure 4—figure supplements 12.

Figure 4—source data 1

Results of manual segmentation: # of ROIs, SNR, and pairwise correlations for simulated and experimental data.

Comparison between data obtained with uncorrected and corrected microendoscopes in silico and in vivo.

https://cdn.elifesciences.org/articles/58882/elife-58882-fig4-data1-v3.xlsx

Pairwise correlation between nearby neurons (distance between the center of neurons <20 µm) should not vary with the radial distance because in our simulations this value was constant across neurons. However, we found an artefactual increase of correlation strength with the radial distance of neuronal pairs in uncorrected endoscopes due to the cross-contamination of activity at different points generated by the larger and aberrated PSF without the corrective lens. In contrast, correlation strength remained constant in eFOV-microendoscopes (Figure 4f,g), suggesting that the PSF of the corrected probes was small enough across the FOV to decrease contamination of activity across neurons. This result is thus in agreement with the decreased spatial resolution observed in more distal parts of the FOV in uncorrected probes and the improved and homogeneous resolution across the FOV that is instead found in corrected microendoscopes (Figure 3). Overall, these findings suggest that signal corresponding to individual neurons could be more accurately extracted from ROIs across the FOV in corrected microendoscopes. We quantified this in synthetic data by evaluating, for each ROI, the correlation of the extracted calcium trace with the ground truth fluorescence signal generated by the simulated neuronal activity contained in that ROI (Figure 4h left). For those ROIs whose fluorescence dynamics were determined by more than one neuron, the correlation with the second most relevant cell was also calculated (Figure 4h right). We found that ROIs segmented from eFOV-microendoscopes displayed larger correlation with the ground truth signal of an individual neuron (Figure 4h left) and lower correlation with the ground truth signal of other nearby neurons (Figure 4h right) compared to uncorrected probes, suggesting that aberration correction allowed to collect more precisely from single cellular emitters and with decreased cross-contamination between neurons.

We experimentally validated the results of the simulations performing functional imaging in the VPM using type II eFOV-microendoscopes before and after the removal of the corrective microlens in awake head-restrained mice (Figure 4i,j). The number of ROIs detected under the two conditions was not significantly different (Figure 4k). However, we found increased average SNR of calcium signals in corrected compared to uncorrected probes (Figure 4l), confirming also in experimental data (as it happens in simulations, Figure 4e) that the use of an eFOV-microendoscope shifts the distribution of SNR across ROIs toward having a higher proportion of ROIs with large SNR value. Moreover, the linear fit of pairwise correlations as a function of radial distance for uncorrected endoscopes had a significantly positive slope (Figure 4m, left panel), indicating higher pairwise correlations in lateral compared to more central portions of the FOV. For corrected endoscopes, the slope of the linear fit was not significantly different from zero (Figure 4m, right panel), in agreement with the analysis of the artificial calcium t-series (Figure 4g). Overall, results of simulations and experiments demonstrate that correcting optical aberrations in eFOV-microendoscopes enabled higher SNR and more precise evaluation of pairwise correlations compared to uncorrected probes.

To evaluate if the segmentation method could affect these results, we compared the quality of the manual segmentation method used in previous experiments with that of a standard automated algorithm (e.g. CaImAn, [Giovannucci et al., 2019]) by computing precision, recall, and F1 score in simulated data (Figure 4—figure supplement 1). We found that the automated method had recall values, on average across SNR values,<0.4 (Figure 4—figure supplement 1a), leading to the detection of only a minority of the total number of neurons. In contrast, the manual method tended to have higher recall across SNR threshold values. Moreover, for low SNR threshold values the automated segmentation had precision values < 0.8 in both uncorrected and corrected endoscopes (Figure 4—figure supplement 1b), leading to identification of ROIs which did not correspond to cells in the ground truth. In contrast, the manual segmentation method tended to have larger values of precision across SNR threshold levels. Overall, F1 scores tended to be higher for the manual segmentation method compared to the automated one for both uncorrected and corrected endoscopes (Figure 4—figure supplement 1c). We extended the comparison between the manual and automated segmentation methods to the real data shown in Figure 4i–m. We observed that in uncorrected endoscopes the automated method identified smaller number of ROIs compared to manual segmentation (Figure 4—figure supplement 1d). In contrast, the number of ROIs identified with the automated approach and the manual method in t-series acquired with the corrected endoscope were not significantly different (Figure 4—figure supplement 1e). One potential explanation of this finding is that the automated segmentation method more efficiently segments ROIs with high SNR compared to the manual one. Since aberration correction significantly increases SNR of fluorescent signals, the automated segmentation performed as the manual segmentation method in corrected endoscopes.

We evaluated the effect of aberration correction on the output of the analysis in the simulated and experimental data shown in Figure 4 using an automated segmentation method (e.g. CaImAn). In simulated data, we found that using CaImAn the number of ROIs segmented in corrected endoscopes was consistently higher than in the uncorrected case across SNR thresholds (Figure 4—figure supplement 2a), similarly to what observed with manual segmentation (Figure 4d). Using CaImAn, SNR values of fluorescence events were significantly higher in corrected compared to uncorrected endoscopes (Figure 4—figure supplement 2b), similarly to what observed with manual segmentation (Figure 4e). Moreover, the slope of the linear fit of pairwise correlations as a function of radial position for uncorrected endoscopes had a significantly positive slope (Figure 4—figure supplement 2c left panel), indicating higher pairwise correlations in lateral compared to more central portions of the FOV. For corrected endoscopes, the slope of the linear fit was not significantly different from zero (Figure 4—figure supplement 2c right panel). This result is also in line with what previously observed with manual segmentation (Figure 4g). In experimental data, we found that using CaImAn SNR values of fluorescence events tended to be higher in corrected compared to uncorrected endoscopes (Figure 4—figure supplement 2d), a trend that was in line with what observed with manual segmentation (Figure 4l). The slope of the linear fit of pairwise correlations as a function of radial position for uncorrected endoscopes was significantly positive (Figure 4—figure supplement 2e), indicating higher pairwise correlations in lateral compared to more central portions of the FOV. For corrected endoscopes, the slope of the linear fit was not significantly different from zero (Figure 4—figure supplement 2e). Both results are in line with what previously observed with manual segmentation (Figure 4m). Overall, the results of the comparison between the manual and automated segmentation methods showed that improvements introduced by aberration correction in endoscopes were observed with both the automated and manual segmentation methods.

Spatial mapping of behavior state-dependent information in sensory thalamic nuclei in awake mice

We then focused our attention on the VPM, a key region which relays somatosensory (whisker) information to the barrel field of the primary somatosensory cortex (S1bf) through excitatory thalamocortical fibers (Feldmeyer et al., 2013). VPM also receives strong cortical feedback from corticothalamic axons of deep cortical layers. Cortical inputs to VPM has been proposed to strongly modulate thalamic activity. Thus, to study VPM physiology it is fundamental to preserve corticothalamic and thalamocortical connectivity. Electrophysiological recordings showed that VPM networks are modulated by whisking and behavioral state (Urbain et al., 2015; Moore et al., 2015; Poulet et al., 2012). However, how information about whisking and other behavioral state-dependent processes (e.g arousal, locomotion) are spatially mapped in VPM circuits at the cellular level is largely unknown. We used eFOV-microendoscopes to address this question.

As an important control experiment, we first confirmed that the ultrathin GRIN lenses that we used in our study (diameter ≤500 µm) did not significantly damage anatomical thalamocortical and corticothalamic connectivity, a difficult task to achieve with larger cross-section GRIN lenses or with chronic optical windows (Figure 3—figure supplement 7e). To this aim, we performed local co-injections in the VPM of Scnn1a-Cre mice of red retrobeads to stain corticothalamic projecting neurons with axons targeting the VPM and of an adenoassociated virus carrying a floxed GFP construct to stain thalamocortical fibers (Figure 5a). We evaluated the amount of thalamocortical and corticothalamic connectivity looking at the percentage of pixels displaying green and red signal in the S1bf region in endoscope-implanted vs. non-implanted mice (Figure 5a–c). In accordance with the known anatomy of the thalamocortical system (Feldmeyer, 2012), we found that the green signal was mostly localized in layer IV barrels and in layer V/VI while the red signal was largely restricted to layer VI (Figure 5d,e). Importantly, we found no difference in the percentage of pixels displaying green and red signals in implanted vs. non-implanted mice (Figure 5f).

Ultrathin microendoscope implantation preserves thalamo-cortical and cortico-thalamic connectivity between S1bf and VPM.

(a) Local injection of red retrobeads and AAVs transducing floxed eGFP was performed in the VPM of Scnn1a-Cre mice. (b) Confocal image showing a coronal slice from an injected control animal. Scale bar: 500 μm. (c) Same as (b) but for a mouse implanted with a type II eFOV-microendoscope (probe diameter: 500 μm). (d,e) Zoom in of the S1bf region highlighted in (b,c). Scale bar: 100 μm. (f) Percentage of labeled S1bf area with eGFP (green) and retrobeads (red) in control and implanted mice. Points indicate the value of fluorescence from three mice (counted three coronal slices from each animal), column bars indicate average ± sd. One-tailed Mann-Whitney, p=0.20 for eGFP and p=0.50 for red retrobeads, respectively.

Figure 5—source data 1

Percentage of labeled S1bf area with eGFP and retrobeads in control and implanted mice.

https://cdn.elifesciences.org/articles/58882/elife-58882-fig5-data1-v3.xlsx

We then used eFOV-microendoscopes to address the question of how information about motor behavior (e.g. locomotion and whisking) and internal states (e.g. arousal state) are mapped on VPM circuits at the cellular level. To this aim, we used eFOV-microendoscopes to perform GCaMP6s imaging in VPM circuits in awake head-restrained mice while monitoring locomotion, whisker mean angle, and pupil diameter (Figure 6a–d, see Materials and methods). We identified quiet (Q) periods, time intervals that were characterized by the absence of locomotion and whisker movements, and active (A) periods, intervals with locomotor activity, dilated pupils, and whisker movements. Active periods were further subdivided into whisking (W), whisking and locomotion (WL), and locomotion with no whisking (L). Figure 6e shows a histogram representing the amount of time spent in the different behavioral states. Mice whisk when they move, therefore L periods were rare. We found that Q periods showed calcium events that were sparsely distributed both across time and neurons (Figure 6c,d). In contrast, active periods displayed an increase in both frequency and amplitude of calcium signals across VPM neurons compared to Q periods (average frequency: fQ = 1.95 ± 0.02 Hz, fA = 2.22 ± 0.02 Hz, Student’s t-test p=2E-74, n = 24 t-series from four mice; average amplitude: AQ = 0.137 ± 0.005 ΔF/F0, AA = 0.245 ± 0.008 ΔF/F0, Student’s t -test p=2E-128, n = 24 t-series from four mice). This resulted in a significant increase in the average fluorescence across neurons during the active W and WL periods compared to Q periods (Figure 6f). The increase in the frequency of WL also correlated with pupil size (Figure 6g and Supplementary file 1 - Table 4), which reflects the arousal level of the animal (Busse et al., 2017).

High-resolution population dynamics in the VPM of awake mice during locomotion and free whisking.

(a) Schematic of the experimental set-up for the recording of locomotion, whisker mean angle, and pupil size in awake head-fixed mice during VPM imaging using type II eFOV-microendoscopes. (b) Two-photon image of GCaMP6s labeled VPM neurons in vivo. (c) Representative traces of locomotion (top), whisker mean angle (middle), and pupil diameter (bottom). Red and blue shades indicate periods of locomotion and whishing, respectively (see Materials and methods for definition). (d), ΔF/F0 over time for all different ROIs in the experiment shown in (c). (e) Percentage of frames spent by the animal in quite wakefulness (Q), whisking (W), locomotion (L), and whisking + locomotion (WL). n = 24 time series from four animals. One-way ANOVA with Bonferroni post-hoc correction, p=3E-23. (f) Average ΔF/F0 across ROIs under the different experimental conditions. n = 24 time series from four animals. One-way ANOVA with Bonferroni post-hoc correction, p=2E-16. (g) Distribution of the Q, W, and WL states as a function of pupil diameter. Kolmogorov-Smirnov test for comparison of distributions of Q, W and WL states: p=0.07 for Q vs. W states, p=2E-10 for Q vs. WL states, p=1E-5 for W vs. WL states. For the statistical comparison of Q, W and WL states in each range of pupil diameter a two-way ANOVA with Tukey-Kramer post-hoc correction was performed; see Supplementary file 1 - Table 4. The error bars represent s.e.m across different t-series.

Figure 6—source data 1

Percentage of frames, ΔF/F0, and distribution of pupil diameter across behavioral states.

https://cdn.elifesciences.org/articles/58882/elife-58882-fig6-data1-v3.xlsx

Cell-specific encoding of whisking-dependent information in distributed VPM subnetworks

We investigated how neuronal activity was modulated by an important behavioral variable: whether the mouse was whisking or not. We considered neuronal activity both at the single-cell and population level. We quantified the content of mutual information about whisking state (whether or not the mouse was whisking; shortened to whisking information hereafter) based on the fluorescence signals extracted from individual neurons (Figure 7a). We found that many neurons were informative about whisking, but only a fraction were particularly informative (Figure 7b). Highly-informative neurons were sparse and could be surrounded by low-information-containing neurons (Figure 7a). This indicates that while informative neurons are distributed across the FOV, information is strongly localized and highly cell-specific. Whisking information in individual cells measured with the eFOV-microendoscope was constant over the radial distance (Figure 7c), with no significant difference in the amount of information in ROI positioned at low radial distance (Dist. low in Figure 7d) and at high radial distance (Dist. high in Figure 7d). Moreover, adding neurons recorded at higher radial distances to a population decoder improved the amount of extracted whisking information compared to considering only neurons in the center of the FOV (Figure 7e). This did not imply that ROIs away from the center carry higher information (information in individual neurons was independent on the cell position within the FOV, Figure 7c), but only that information of cells in lateral portions of the FOV summed to the information of more central ROIs in a way that was not purely redundant. Together, these results suggest that the corrected endoscope was effective at capturing information also at the FOV borders. Information of single cells correlated positively with the SNR of their calcium signals (Figure 7f), with higher amount of information carried, on average, by cells displaying higher SNR (Figure 7g). This suggests that the benefit of eFOV-microendoscopes in providing higher SNR, demonstrated earlier, also translates in an ability to extract more information about the circuit’s encoding capabilities (Figure 7e).

Figure 7 with 1 supplement see all
Cell-specific encoding of behaviorally-dependent information in distributed VPM subnetworks.

(a) Spatial map of neurons encoding whisking information. The pseudocolor scale shows significantly informative neurons (see Materials and methods). Data are pooled from 24 t-series from four animals. See also Figure 7—figure supplement 1. (b) Distribution of information content in individual ROIs. The individual ROIs information content distribution was fitted using a double exponential function (R2 = 0.99). n = 808 ROIs from 24 time series. (c) Information content of individual neurons vs. their radial distance from the center of the FOV. n = 842 neurons from 24 time series. Linear regression fit: slope = −3E-5, Student’s t-test p=0.13. Pearson correlation coefficient: −0.053, Student’s t-test p=0.13. (d) Information content for ROIs located at low radial distance (distr low) and in lateral portions of the FOV (distr high). n = 425 and 417 ROIs from 24 t-series for distr low and distr high, respectively. Mann-Whitney test, p=0.42. (e) Information content of neural populations as a function of the distance from the FOV center (% radius) used to include ROIs in the population. One-way ANOVA repeated measurements with Bonferroni post hoc, p=1E-18. Data pooled from 24 time series. (f) Information content of individual neurons vs. SNR. n = 842 neurons from 24 time series. Linear regression fit: slope = 9E-4, Student’s t-test p=2.7E-4. Pearson correlation coefficient: 0.125, p=2.7E-4. (g) Information content for ROIs with low SNR (SNR low) and high SNR (SNR high). n = 421 and 421 ROIs from 24 t-series for distr low and distr high, respectively. Mann-Whitney test, p=5E-5. (h,i) Synergistic (gray) or redundant (dark green) information within pairs of neurons is shown as a function of pairwise correlation (h) and as a function of radial distance (i). n = 61 and 31 pairs of neurons for synergistic and redundant information content, respectively. Data from 24 time series. Linear regression fit in (h): slope = 0.008, permutation test p=0.45 and slope = −0.02, permutation test p=0 for synergistic and redundant information, respectively. Linear regression fit in (i): slope = −0.0001, permutation test p=0.007 and slope = 2E-5, permutation test p=0.38 for synergistic and redundant information, respectively. (j) Top: representative grayscale matrix showing the activation coefficient across time for 4 NMF modules emerging in a FOV containing 37 neurons. Periods of whisking are shown in red bars and shades. Bottom: ROIs belonging to the four different modules shown in the top panel. Each colored circle represents a ROI belonging to the specified module and its radius is proportional to the ROI weight within that module. The corresponding activation coefficients are presented in the upper panel. (k) Module sparseness as a function of the whisking modulation index (WMI). n = 213 modules from 24 time series. Pearson correlation coefficient: −0.164, p=0.016. (l) Normalized number of modules as a function of their spatial spread for the experimental data (gray) and for a null distribution obtained by randomly shuffling the spatial position of ROIs in the FOV (white). Number of modules: n = 213 and n = 2400 for the 24 t-series of the in vivo dataset and for the null distribution, respectively. (m) Left: Normalized number of modules as a function of their Jaccard index for pairs of modules identified within the same FOV for experimental data (gray) and for a null distribution (white) obtained by randomly shuffling the spatial position of ROIs within the FOV. Number of modules: n = 1704 and n = 34080 for the 24 t-series of the experimental data and for the null distribution, respectively. Right: same as in the left panel for Cosine similarity coefficients.

Figure 7—source data 1

Information theoretical analysis and non-negative matrix factorization results.

Values of cell-specific and population information analysis.

https://cdn.elifesciences.org/articles/58882/elife-58882-fig7-data1-v3.xlsx

We also considered the redundancy and synergy of whisking information carried by pairs of simultaneously recorded nearby (distance between neurons < 20 µm) neurons. This analysis is important because how pairwise correlations shape the redundancy and synergy of information representation is fundamental to the understanding of population codes (Runyan et al., 2017; Rumyantsev et al., 2020; Averbeck et al., 2006). We found that redundancy and synergy of pairs recorded with the eFOV-microendoscopes had, on average, more negative and more positive values in pairs with stronger correlations, respectively (Figure 7h). Importantly, and as a consequence of the fact that eFOV-microendoscopes avoid the artificial increase of pairwise correlations close to the FOV borders, we did not observe an increase of synergy or redundancy between nearby neurons close to the FOV border (Figure 7i). This shows that aberration correction helps avoiding the generation of artificially biased estimates of synergy and redundancy near the FOV border.

We finally turned to analyzing the properties of firing at the level of the whole population recorded in the FOV. We applied non-negative matrix factorization (NMF) (Lee and Seung, 1999) to identify subpopulations of neurons (modules, Figure 7j) characterized by correlated activity. Detected modules were differentially activated in time and were sparsely distributed in space (Figure 7j). Moreover, we found that modules could be oppositely modulated by whisking, with the activity of some modules being enhanced and the activity of some other modules being depressed by whisking (Figure 7j). We computed the whisking modulation index (WMI, see Materials and methods for definition) and found that the large majority of modules was positively modulated by whisking (WMI >0 for 89.6 ± 0.4% of total modules number, n = 24 FOVs), while the activity of a minority of modules was suppressed during whisking periods (negatively modulated, WMI <0 for 10.4 ± 0.4% of total modules number, n = 24 FOVs). Sparseness of modules appeared to be negatively correlated with the WMI, suggesting that those few modules that were negatively modulated by whisking were also characterized by few, but highly-informative neurons (within the ensemble) (Figure 7k). In contrast, modules with high WMI values were less sparse, suggesting similar activity (and information) across most of the neurons belonging to these ensembles. Single modules covered distances of hundreds of μm spanning the whole eFOV (Figure 7l). We showed that the spatial distances covered by functional modules were higher than distances obtained by chance, using a permutation test (Figure 7l). This suggests that corrected probes allow unveiling functional relationship between groups of neurons spanning the whole eFOV. Neurons could belong to different ensembles as quantified by the distribution of the values of the Jaccard index among pairs of modules (Figure 7m, left). This distribution showed a peak toward the value 1 (the two modules were composed by the same ROIs). However, when ROIs belonged to more than one module they tended to have module-specific weight (i.e. different weights for different modules). In fact, the distribution of values for the cosine similarity, an index which considers the weight of ROIs within a module (see Materials and methods), was shifted toward smaller values compared to the distribution of the Jaccard index values (Figure 7m).

Discussion

Improving optical performances in ultrathin (diameter ≤0.5 mm) microendoscopes with built-in optical elements is a major technical challenge. Since the insertion of the probe irreversibly damages the tissue above the target area, reducing the size of the probe and consequently its invasiveness is of utmost importance when imaging deep brain regions. In this study, we designed, developed, characterized, and successfully validated a new approach to correct aberrations in ultrathin GRIN-based endoscopes using aspheric lenses microfabricated with 3D micro-printing based on TPL (Liberale et al., 2010).

Corrective lenses were fabricated on glass coverslips, which were aligned and assembled with the GRIN rod to form an aberration-corrected microendoscope. This optical design resulted in improved axial resolution and extended effective FOV (Figure 3), in good agreement with the predictions of the optical simulations (Figure 2). The absolute values of the axial PSF computed from optical simulations (Figure 2—figure supplements 12) were generally smaller in the optical simulations compared to real measurements (Figure 3 and Figure 2—figure supplements 12). This may be due to multiple reasons. First, high-order aberrations were not included in the simulations. Second, in simulations, although the intensity of the excitation PSF was small in lateral portion of the FOV (Figure 2), a Gaussian function could still well fit the dim intensity distribution and provide a clear quantification of the PSF dimension. In experimental measurements of fluorescence emitted by subresolved beads, the more degraded PSF in the lateral portions of the FOV would result in low efficacy of the excitation beam in stimulating fluorescence which would result in low SNR fluorescence signals and introduce large variability in the fit. Third, variabilities in some of the experimental parameters (e.g. the distance between the GRIN back end and the focusing objective) were not considered in the simulations.

Aberration correction in GRIN microendoscopes can be achieved using adaptive optics (AO) (Wang and Ji, 2013; Wang and Ji, 2012; Bortoletto et al., 2011; Lee and Yun, 2011). For example, using pupil-segmentation methods for AO, diffraction-limited performance across an enlarged FOV was obtained in GRIN-based endoscopes with diameter of 1.4 mm (Wang and Ji, 2013; Wang and Ji, 2012) and, in principle, this approach could be extended to probes with smaller diameter. Nevertheless, AO through pupil segmentation requires significant modification of the optical setup and the use of an active wavefront modulation system (e.g. a liquid crystal spatial light modulator) which needs the development of ad-hoc software control. Moreover, AO through pupil segmentation may limit the temporal resolution of the system, since multiple AO corrective patterns must be applied to obtain an aberration-corrected extended FOV (Wang and Ji, 2013). Compared to AO approaches, the technique developed in this study does not necessitate modification of the optical path nor the development of ad-hoc computational approaches. Moreover, it is easily coupled to standard two-photon set-ups, and does not introduce limitations in the temporal resolution of the imaging system. A potential alternative to the approach described in this study would be to place a macroscopic optical element of the desired profile in a plane optically conjugated to the objective back aperture along the optical path. This solution could have the advantage of being manufactured using a more standard techniques. However, it would require significant change in the optical set-up in contrast to the built-in correction method that we describe in the present study. Moreover, this macroscopic optical element would have to be changed according to the type of microendoscope used.

Using synthetic calcium data, we demonstrated that the improved optical properties of eFOV-microendoscopes directly translate into important advantages for measuring neural population parameters in vivo. Namely, they achieve a higher SNR of calcium signals and a more precise evaluation of pairwise correlations compared to uncorrected GRIN lenses, predictions that were all confirmed experimentally in awake mice (Figure 4). Importantly, synthetic calcium data also allowed us to evaluate the impact of correcting optical aberrations on the accuracy in extracting neuronal activity and population codes from calcium imaging data. We found larger correlation of extracted calcium traces with the known ground truth of neuronal spiking activity in eFOV-microendoscopes compared to uncorrected probes. Traces extracted from eFOV-microendoscopes were less contaminated by neighboring neurons compared to uncorrected probes, in agreement with the higher and more homogeneous spatial resolution of eFOV-microendoscopes (Figure 4). All these achievements were obtained without increasing the lateral size of the probe, thus minimizing tissue damage in biological applications.

Studying neuronal population codes requires the measurement of neuronal population activity with high-precision, large SNR, and without introducing artificial bias on the activity of individual neurons and the measures of relationship between them, such as pairwise correlations. In particular, pairwise correlations are thought to be fundamental for information coding, signal propagation, and behavior (Ni et al., 2018; Salinas and Sejnowski, 2001; Shahidi et al., 2019; Panzeri et al., 1999; Runyan et al., 2017). Here, we demonstrate that the homogeneous spatial resolution, which characterized eFOV-microendoscopes (Figure 3), allowed an unbiased computation of pairwise correlations, a higher correlation of extracted calcium traces with the ground truth neuronal activity, and a smaller contamination of extracted signals by neighboring cells (Figure 4). Several studies suggested that even small biases in measuring single cell and pairwise properties, either, for example, in incorrectly measuring the average amount of correlations or the heterogeneity of single cell tuning and of correlation values, may lead to large biases in determining how these populations encode information (Shamir and Sompolinsky, 2006; Panzeri et al., 1999; Ecker et al., 2011). eFOV-microendoscopes allow to remove the artifacts and biases introduced by uncorrected GRIN endoscopes in measuring both individual cell information properties and correlation between each and every pair of neurons. The advantages introduced by eFOV-microendoscopes are therefore essential for unraveling the true nature of population codes in deep brain structure.

Corrected endoscopes were characterized by a curved FOV. In the case of type II corrected endoscopes, the change in the z coordinate in the focal plane was up to ~75 µm (Figure 3). This z value was smaller for all other corrected endoscope types (Figure 3). The observed field curvature of corrected endoscopes may impact imaging in brain regions characterized by strong axially organized anatomy (e.g. the pyramidal layer of the hippocampus), but would not significantly affect imaging in regions with homogeneous cell density within the z range described above (<75 µm for type II corrected microendoscopes).

We used the unique features of the eFOV-microendoscopes to study how highly correlated activity is mapped in the VPM thalamic nucleus of awake mice with unprecedented combination of high spatial resolution across the FOV and minimal invasiveness (Figures 57). The VPM is a primary somatosensory thalamic nucleus which relays sensory information to the S1bf (Sherman, 2012). However, VPM also receives strong cortical innervation which deeply affects VPM activity (Crandall et al., 2015; Mease et al., 2014; Temereanca and Simons, 2004). We first showed that the small cross-section of the eFOV-microendoscopes developed here preserved thalamocortical and corticothalamic connectivity (Figure 5), a fundamental prerequisite for VPM physiology and a hardly achievable task with larger cross-section GRIN lenses or with chronic windows (Figure 3—figure supplement 7e). We then imaged GCaMP6s-expressing VPM neurons while monitoring locomotion, whisker movement, and pupil diameter (Figure 6). We found cell-specific encoding of whisking information in distributed functional VPM subnetworks. Most individual neurons encoded significant amount of whisking information generating distributed networks of informative neurons in the VPM (Figure 7). However, the amount of encoded information was highly cell-specific, with high-information-containing neurons being sparsely distributed in space and surrounded by low-information-containing cells. Sparse distribution of information content has been similarly observed in other brain areas (Runyan et al., 2017; Ince et al., 2013). At the population level, we observed the presence of whisking modulated functional ensembles of neurons, which were oppositely modulated by whisking. Some ensembles displayed enhanced activity upon whisking, while some other ensembles showed suppressed activity by whisking. Interestingly, single neurons could belong to multiple functional ensembles, but their weight within one ensemble was ensemble-specific (Figure 7). Overall, the application of eFOV-microendoscopes revealed the complexity and cellular specificity of the encoding of correlated behavioral state-dependent information in a primary thalamic sensory nucleus.

One important area of future development for eFOV-microendoscopes will be to determine whether the approach we described in this study could be used to correct optical aberrations in GRIN lenses different from the ones described here (e.g. longer GRIN rods). A second area of interest will be to develop corrective lenses or compound corrective lenses for 3D imaging in larger volumes. This should be possible given that, in the present study, we designed the corrective lenses to maximize aberration correction only in the focal plane of the endoscope and in optical simulations we found that the Strehl ratio was >0.8 in a 1.7–3.7 times larger volume in corrected compared to uncorrected endoscopes (Figure 2—figure supplement 3). Finally, new fabrication materials (Weber et al., 2020) may allow to develop solutions for effective chromatic aberrations (Figure 2—figure supplement 4) in two-photon multi-wavelength applications.

In summary, we developed a new methodology to correct for aberrations in ultrathin microendoscopes using miniaturized aspheric lenses fabricated with 3D printing based on TPL. This method is flexible and can be applied to the GRIN rods of different diameters and lengths that are required to access the numerous deep regions of the mammalian brain. Corrected endoscopes showed improved axial resolution and up to nine folds extended effective FOV, allowing high- resolution population imaging with minimal invasiveness. Importantly, we demonstrated that eFOV-microendoscopes enable more precise extraction of population codes from two-photon imaging recordings. Although eFOV-microendoscopes have been primarily applied for functional imaging in this study, we expect that their use can be extended to other applications. For example, eFOV-microendoscopes could be combined with optical systems for two-photon holographic optogenetic manipulations (Packer et al., 2012; Rickgauer and Tank, 2009; Papagiakoumou et al., 2010) and for simultaneous functional imaging and optogenetic perturbation (Packer et al., 2015; Rickgauer et al., 2014; Carrillo-Reid et al., 2016; Forli et al., 2018; Marshel et al., 2019) using a diffractive optical element to provide patterned illumination of neurons but also to correct for z-defocus with an appropriate lens function. Moreover, besides its applications in the neuroscience field, eFOV-microendoscopes could be used in a large variety of optical applications requiring minimally invasive probes, ranging from cellular imaging (Kim et al., 2012; Ghosh et al., 2011) to tissue diagnostic (Huland et al., 2012; Huland et al., 2014). Importantly, applications of ultrathin eFOV-microendoscopes to other fields of research will be greatly facilitated by the built-in aberration correction method that we developed. This provides a unique degree of flexibility that allows using ready-to-use miniaturized endoscopic probes in a large variety of existing optical systems with no modification of the optical path.

Materials and methods

Key resources table
Reagent type
(species)
or resource
DesignationSource or referenceIdentifiersAdditional information
Strain, strain
background
(M. musculus)
C57BL/6JCharles RiverRRID:IMSR_JAX:000664
Genetic reagent
(M. musculus)
B6;C3-Tg(Scnn1a-
cre)3Aibs/J
The Jackson LaboratoryRRID:IMSR_JAX:009613
Recombinant
DNA reagent
pAAV.Syn.Flex.
GCaMP6s.WPRE.SV40
Penn Vector CoreRRID:Addgene_100845;
Addgene viral prep
# 100845-AAV1
Chen et al., 2013
Recombinant
DNA reagent
pGP-AAV-syn-
FLEX-jGCaMP7f-WPRE
AddgeneRRID:Addgene_104492;
Addgene viral prep
# 104492-AAV1
Dana et al., 2016
Recombinant
DNA reagent
AAV pCAG-
FLEX-EGFP-WPRE
Penn Vector CoreRRID:Addgene_51502;
Addgene viral prep
# 51502-AAV1
Oh et al., 2014
Recombinant
DNA reagent
AAV.CaMKII0.4.
Cre.SV40
Penn Vector CoreRRID:Addgene_105558;
Addgene viral prep
# 105558-AAV1
Commercial
assay or kit
Kwik-CastWorld Precision InstrumentsCat# KWIK-CAST
Commercial
assay or kit
Sylgard Silicone
Elastomer
Dow IncCat# Sylgard 164
Commercial
assay or kit
Norland Optical
Adhesive 63
NorlandCat# NOA 63
Commercial
assay or kit
GRIN lensGrintechCat# NEM-050-
25-10-860-S
Commercial
assay or kit
GRIN lensGrintechCat# NEM-050-
43-00-810-S-1.0p
Commercial
assay or kit
GRIN lensGrintechCat# GT-IFRL-035-
cus-50-NC
Commercial
assay or kit
GRIN lensGrintechCat# NEM-035-
16air-10–810 S-1.0p
Chemical compound,
drug
bisBenzimide H 33342
trihydrochloride (Hoechst)
Sigma-AldrichCat# B2261; CAS: 23491-52-3
Chemical compound,
drug
Red RetrobeadsLumaFluor IncRed Retrobeads
Software, algorithmZemax OpticStudio 15Zemaxhttps://www.zemax.com/products/opticstudio
Software, algorithmMATLAB R2017aMathworksRRID:SCR_001622;
https://it.mathworks.com/products/matlab.html
Software, algorithmGraphPad PRISMGraphPad PRISMRRID:SCR_002798;
https://www.graphpad.com/
Software, algorithmImageJ/FijiFijiRRID:SCR_002285;
http://fiji.sc/
Software, algorithmNoRMCorrePnevmatikakis and Giovannucci, 2017https://github.com/flatironinstitute/NoRMCorre
Software, algorithmCaImAnGiovannucci et al., 2019https://github.com/flatironinstitute/CaImAn-MATLAB
Software, algorithmPopulation Spike Train
Factorization Toolbox
for Matlab Version 1.0
Onken et al., 2016https://stommac.eu/index.php/code
Software, algorithmLIBSVMChang and Lin, 2011https://www.csie.ntu.edu.tw/~cjlin/libsvm/
Software, algorithmInformation
Breakdown ToolBox
Magri et al., 2009N/A
Software, algorithmSoftware used in this
paper for generation of
artificial time series
https://github.com/moni90/eFOV_microendoscopes_simFigure 4a–h,
Figure 4—figure supplement 1a–c
and
Figure 4—figure supplement 2a–c
Software, algorithmSoftware to compute
recall, precision,
and F1 score
Soltanian-Zadeh et al., 2019https://github.com/soltanianzadeh/STNeuroNet
OtherBasler ace cameraBasler AGCat# acA800-510um
OtherOptical encoderBroadcomAEDB-9140-A13
OtherZortrax M200 3D printerZortraxM200
OtherZ-ULTRAT 3D
printer filament
ZortraxZ-ULTRAT
OtherArduino UnoArduinoArduino Uno

Animal models

Request a detailed protocol

Experimental procedures involving animals have been approved by the Istituto Italiano di Tecnologia Animal Health Regulatory Committee, by the National Council on Animal Care of the Italian Ministry of Health (authorization # 1134/2015-PR, # 689/2018-PR) and carried out according to the National legislation (D.Lgs. 26/2014) and to the legislation of the European Communities Council Directive (European Directive 2010/63/EU). Experiments were performed on adult (8–14 week old) mice. C57BL/6J mice (otherwise called C57, Charles River #000664, Calco, IT) were used in Figure 3—figure supplement 6. Data reported in Figures 47 and Figure 3—figure supplement 7 were obtained from B6;C3-Tg(Scnn1a-cre)3Aibs/J (JAX #009613, Jackson Laboratory, Bar Harbor, USA) mice crossed with C57 mice (otherwise called Scnn1a-Cre). Both male and female animals were used in this study. Animals were housed in individually ventilated cages under a 12 hr light:dark cycle. Access to food and water was ad libitum. The number of animals used for each experimental set of data is specified in the text or in the corresponding Figure legend.

Methods details

Design and simulation of corrective lenses and of eFOV-microendoscopes

Request a detailed protocol

Simulations were run with OpticStudio15 (Zemax, Kirkland, WA) to define the profile of the aspheric corrective lens to be integrated in the aberration-corrected microendoscopes, with the aim to achieve: (i) a full-width half maximum (FWHM) lateral resolution <1 µm at the center of the FOV; (ii) a FWHM axial resolution below <10 µm; (iii) a working distance between 150 µm and 220 µm into living brain tissue. The wavelength used for simulations was λ = 920 nm. The surface profile of corrective aspheric lenses was described in Optic Studio Manual, 2017:

(1) Z(r)=cr21+1-1+kc2r2+nαnr2n

Since GRIN lenses have intrinsic spherical aberration, the optimization for the shape of the corrective lenses started with the profile of a Schmidt corrector plate (Born and Wolf, 1999) as initial guess; the parameters c, k, αn (with n = 1–8) in Equation (1) were then automatically varied in order to maximize the Strehl ratio (Dorband et al., 2012) over the largest possible area of the FOV (Supplementary file 1 - Table 1). A fine manual tuning of the parameters was performed for final optimization. Ray-trace simulations were performed considering the material used in lens replica (i.e. NOA63). Simulated two-photon PSFs (Figure 2—figure supplements 12) were determined by 3D sampling the squared calculated Strehl ratio. Intensity profiles derived from sections in the x, y, z directions of the simulated two-photon PSFs were fitted with Gaussian curves and analyzed as for experimental PSFs (see below).

To evaluate the effect of corrective lenses on the 3D image space (Figure 2—figure supplement 3), the Strehl ratio was calculated for different defocused working distances, which were simulated by changing the spacing between the microscope objective and the microendoscope back end. Different radial distances were simulated on the same fields used in Figure 2. To evaluate the axial chromatic aberrations (Figure 2—figure supplement 4), the Strehl ratio was calculated for different wavelengths. For those wavelengths, the defocused working distance which maximizes the Strehl ratio, was also determined.

Corrective lens manufacturing and microendoscope assembly

Request a detailed protocol

The optimized aspheric lens structure obtained with simulations was exported into a 3D mesh processing software (MeshLab, ISTI-CNR, Pisa, IT) and converted into a point cloud dataset fitting the lens surface (with ~300 nm distance among first neighborhood points). Two-photon polymerization with a custom set-up (Liberale et al., 2010) including a dry semi-apochromatic microscope objective (LUCPlanFLN 60x, NA 0.7, Olympus Corp., Tokyo, JP) and a near infrared pulsed laser beam (duration: 100 fs; repetition rate: 80 MHz; wavelength: 780 nm; FemtoFiber pro NIR, Toptica Photonics, Graefelfing, DE) was used for the fabrication of the corrective lenses. A drop of resin (4,4’-Bis(diethylamino)benzophenone photoinitiator mixed with a diacrylate monomer), sealed between two coverslips, was moved by a piezo-controlled stage (model P-563.3CD, PI GmbH, Karlsruhe, DE) with respect to the fixed laser beam focus, according to the 3D coordinates of the previously determined point cloud with precision of 20 nm. Output laser power was ~15 mW at the sample. Once the surface was polymerized, the lens was dipped for ~2 min in methanol followed by ~1 min immersion in isopropyl alcohol and finally exposed to UV light (λ = 365 nm; 3 Joule / cm2[Salinas and Sejnowski, 2001]) to fully polymerize the bulk of the structure.

For fast generation of multiple lens replicas, a molding (Schaap and Bellouard, 2013) technique was used. To this end, polydimethylsiloxane (PDMS, Sylgard 164, 10:1 A:B, Dow Corning, Auburn, MI) was casted onto the lens and hardened by heat cure in a circulating oven at 80°C for approximately 30 min. The resulting bulked structure of solid PDMS was then used as negative mold. A drop of a UV-curable optically-clear adhesive with low fluorescent emissivity (NOA63, Norland Products Inc, Cranbury, NJ) was deposited on the negative mold, pressured against a coverslip (diameter: 5 mm) of appropriate thickness (thickness: 100 or 200 µm depending on the eFOV-microendoscope type, Figure 1) and hardened by UV exposure. One side of the UV-curable adhesive was in contact with the mold, the other side was instead attached to the coverslip. After UV curing, by gently pulling the glass coverslip away, the lens made of UV-curable adhesive detached easily from the PDMS mold, while remaining firmly attached to the coverslip. The coverglass with the lens attached was then glued onto a metal ring. The yield for 3D printed lenses and lens replica was ~100%. One end of the appropriate GRIN rod (NEM-050-25-10-860-S, type I; NEM-050-43-00-810-S-1.0p, type II; GT-IFRL-035-cus-50-NC, type III; NEM-035-16air-10–810 S-1.0p, type IV, Grintech GmbH, Jena, DE) was attached perpendicularly to the other surface of the coverslip using NOA63. Alignment of the corrective lens and the GRIN rod was performed under visual guidance using an opto-mechanical stage, custom-built using the following components (Figure 3—figure supplement 2): camera (DCC1645C, Thorlabs, Newton, NJ), fine z control (SM1Z, Thorlabs, Newton, NJ), coarse z control (L200/M, Thorlabs, Newton, NJ), xyz control (MAX313D/M, Thorlabs, Newton, NJ), high power UV LED (M375L3, Thorlabs, Newton, NJ), long pass dichroic mirror (FF409-Di02, Semrock, Rochester, NY), tube lens (AC254-150-A, Thorlabs, Newton, NJ), objective (UPlanFLN 4×, 0.13NA, Olympus, Milan, IT), xy control (CXY1, Thorlabs, Newton, NJ), custom GRIN rod holder, and fiber optic holder (HCS004, Thorlabs, Newton, NJ). An additional and removable coverglass or a silicone cap (Kwik-Cast Sealant, World Precision Instruments, Friedberg, DE) was glued on the top of every support ring to keep the polymeric corrective lens clean and to protect it from mechanical damage.

Optical characterization of eFOV-microendoscopes

Request a detailed protocol

Optical characterization of eFOV-microendoscopes was carried out with a two-photon laser-scanning microscope equipped with a wavelength-tunable, ultrashort-pulsed, mode-locked Ti:Sapphire laser source (Ultra II Chameleon, pulse duration: 160 fs; repetition rate: 80 MHz; wavelength: 920 nm; Coherent Inc, Santa Clara, CA) and a commercial Prairie Ultima IV scanhead (Bruker Corporation, Milan, IT). For all measurements, the wavelength was set at 920 nm. The optomechanical assembly used for the eFOV-microendoscope characterization is shown in Figure 3—figure supplement 2c. The coupling objective was EC Epiplan-Neofluar 20x, 0.5NA (Zeiss, Oberkochen, DE). The z control (SM1Z) and xy control (CXY2) were purchased from Thorlabs (Newton, NJ). Spatial resolution of each microendoscope was evaluated using subresolved spherical fluorescent beads (diameter: 100 nm, Polyscience, Warrington, PA), following a previous spatial calibration using a custom fluorescent ruler (Motic, Xiamen, CN). The same ruler was used to evaluate the distortion of the FOV. To visualize the curvature of the imaging field, thin (thickness: 300 nm) fluorescent slices (Antonini et al., 2014) were used. Fluorescent samples were deposited on a microscope slide and imaged through the endoscope assembly aligned to the microscope objective, with or without the corrective microlens above the coverslip, using the coupling apparatus described in Figure 3—figure supplement 2. Imaging was performed with the distal end of the GRIN rod immersed in a droplet of water placed on the slide. We observed no appreciable damage on the lens over imaging sessions. In considering this, please note that the distance between the focal plane of the microscope objective and the endoscope assembly was ~100 µm and it was fixed for all measurements. Given the imaging field curvature of endoscopes, for both the ruler and the thin fluorescent slices (planar samples), the acquisition of z-series of images (512 pixels x 512 pixels, with 1 µm axial step) was performed.

Viral injections and microendoscope implantation

Request a detailed protocol

Adeno-associated viruses (AAVs) AAV1.Syn.flex.GCaMP6s.WPRE.SV40, AAV1.CAG.Flex.eGFP.WPRE.bGH, AAV1.CaMKII0.4.Cre.SV40 were purchased from the University of Pennsylvania Viral Vector Core. AAV1.Syn.flex.GCaMP7f.WPRE.SV40 was purchased from Addgene (Teddington, UK) Animals were anesthetized with isoflurane (2% in 1 L/min O2), placed into a stereotaxic apparatus (Stoelting Co, Wood Dale, IL) and maintained on a warm platform at 37°C. The depth of anesthesia was assessed by monitoring respiration rate, heartbeat, eyelid reflex, vibrissae movements, reactions to tail and toe pinching. 2% lidocaine solution was injected under the skin before surgical incision. A small hole was drilled through the skull and 0.5–1 µl (30–50 nl/min, UltraMicroPump UMP3, WPI, Sarasota, FL) of AAVs containing solution was injected at stereotaxic coordinates: 1.4 mm posterior to bregma (P), 1 mm lateral to the sagittal sinus (L), and 1 mm deep (D) to target the hippocampal CA1 region; 1.7 mm P, 1.6 mm L, and 3 mm D to target the VPM. Co-injection of AAV1.Syn.flex.GCaMP6s.WPRE.SV40 and AAV1.CaMKII0.4.Cre.SV40 (1:1) was performed to express GCaMP6s in hippocampus CA1 pyramidal cells of C57 mice (Figure 3—figure supplement 6). Injection of AAV1.Syn.flex.GCaMP6s.WPRE.SV40 (1:4 in saline solution) (Figures 67 and Figure 3—figure supplement 7) or AAV1.Syn.flex.GCaMP7f.WPRE.SV40 (1:4 in saline solution) (Figure 4) in the Scnn1a-Cre mice was performed to express GCaMP6/7 in the VPM. Following virus injection a craniotomy (~600×600 µm2 or ~400×400 µm2 depending on the endoscope size) was performed over the neocortex at stereotaxic coordinates: 1.8 mm P and 1.5 mm L to image the hippocampus; 2.3 mm P and 2 mm L to reach the VPM. A thin column of tissue was suctioned with a glass cannula (ID, 300 µm and OD, 500 µm; Vitrotubs, Vitrocom Inc, Mounting Lakes, NJ) and the microendoscope was slowly inserted in the cannula track, using a custom holder, down to the depth of interest and secured by acrylic adhesive and dental cement to the skull. If necessary, metal spacers (thickness:~100 µm) were glued on the flat coverslip surface to obtain the desired protrusion distance of the GRIN rod. For experiments in awake animals (Figures 4, 6 and 7 and Figure 3—figure supplement 7), a custom metal head plate was sealed on the skull using dental cement to assure stable head fixation during two-photon imaging. An intraperitoneal injection of antibiotic (BAYTRIL,Bayer, DE) and dexamethasone (MSD Animal Health, Milan, IT) was performed to prevent infection and inflammation. Animals were then positioned under a heat lamp and monitored until recovery.

To evaluate thalamocortical (TC) and corticothalamic (CT) anatomical connectivity (Figure 5) after implantation of microendoscopes, Scnn1a-Cre mice were injected, as described above, at 100 nl/min (UltraMicroPump UMP3, WPI, Sarasota, FL) with AAV1.CAG.Flex.eGFP.WPRE.bGH virus (1:2 in saline solution) and red retrobeads (1:8 in saline solution, Fluorescent Latex Microspheres, LumaFluor Inc, Durham, NC) in VPM after tissue suctioning and in the absence of tissue suctioning. The total injected volume was 250 nl.

Functional imaging with eFOV-microendoscopes in vivo

Request a detailed protocol

For experiments in anesthetized conditions (Figure 3—figure supplement 6), three to five weeks after injection, mice were anesthetized with urethane (16.5%, 1.65 g*kg−1) and placed into a stereotaxic apparatus to proceed with imaging. Body temperature was measured with a rectal probe and kept at 37°C with a heating pad. Depth of anesthesia was assured by monitoring respiration rate, eyelid reflex, vibrissae movements, and reactions to pinching the tail and toe. In some experiments, oxygen saturation was controlled by a pulseoxymeter (MouseOx, Starr Life Sciences Corp., Oakmont, PA).

For experiments in behaving mice (Figures 4, 6 and 7 and Figure 3—figure supplement 7), imaging was performed two to four weeks after the endoscope implant, following 7–10 days of habituation, in which mice were placed daily on the set-up, each day for a longer time duration, up to 45 min (Gentet et al., 2010). Mice were allowed to run spontaneously on the wheel. During experiments, recording sessions were up to five 5 min long (frame rate typically ~3 Hz) and they were interleaved by 5 min in which no imaging was performed. We did not collect fluorescent signals at frame rates higher than 4 Hz. This was because we aimed at imaging the largest possible FOV and our experimental set-up was equipped with regular galvanometric mirrors. For scanning imaging of GCaMP6/7-expressing neurons, the same microscope set-up used for the optical characterization of eFOV-microendoscopes was used and GCaMP6/7 fluorescence was excited at 920 nm (laser power: 28–90 mW).

Measurement of whisker angle, pupil size, and locomotion

Request a detailed protocol

In Figures 67, whisker movements in the whisker pad contralateral to the recording site were imaged with a high-speed,Basler acA800 (Basler, Ahrensburg, DE; acquisition rate: 150 Hz) through a 45° tilted mirror placed below the whiskers. Illumination was provided by an array of infrared LEDs (emission wavelength: 800 nm) fixed to the microscope objective and aligned to the whiskers and the mirror. Imaging of the contralateral eye was performed during each experimental session with a Basler acA800 (Basler, Ahrensburg, DE), coupled with a long pass filter (HP 900, Thorlabs, Newton, NJ) at 80 Hz acquisition rate. Illumination light was provided by the pulsed laser sourced which was used to perform two-photon microendoscopic imaging (λ = 920 nm). Locomotor activity was measured with an optical encoder (AEDB-9140-A13, Broadcom, San Jose, CA) mounted under the wheel.

Immunohistochemistry

Request a detailed protocol

Deeply anesthetized animals were transcardially perfused with 0.01 M PBS (pH 7.4) followed by 4% paraformaldehyde. Brains were post-fixed for 6 hr, cryoprotected with 30% sucrose solution in 0.1 M PBS, and serially cut in coronal sections (thickness: 40–50 µm) using a HM 450 Sliding Microtome (Thermo Fisher). Sections were counterstained with Hoechst (1:300, Sigma–Aldrich, Milan, IT), mounted, and coverslipped with a DABCO [1,4-diazobicyclo-(2,2,2)octane]-based antifade mounting medium. Fluorescence images were acquired with a Leica SP5 inverted confocal microscope (Leica Microsystems, Milan, IT).

For the evaluation of the anatomical connections between VPM and S1bf, mice were perfused after 10 days from the injection and GRIN lens implantation. 50 µm thick coronal brain slices were cut, counterstained with Hoechst (1:300, Sigma–Aldrich, Milan, IT), and mounted with an Antifade Mounting Medium (Vectashield, Burlingame, CA). Confocal images were acquired with a Nikon Eclipse scope (Nikon, Milan, IT).

Simulations

Request a detailed protocol
Geometrical considerations
Request a detailed protocol

In Figure 4a–h, neurons were simulated as spheres with an average radius r_mean which was estimated from recorded data (r_mean = 7.95 ± 2.33 µm, mean ± sd). Some variability was introduced in the neurons size, sampling it from a normal distribution with mean r_mean and standard deviation r_sigma = 1.31 (within the measured sd). At the center of each neuron a nuclear region was added, such that the spherical shell surrounding the nucleus had a width of variable size (r_shell randomly sampled from a normal distribution with mean = 4 µm and sd = 1 µm). The nuclear region did not express GCaMP6s and the fluorescence signal could be collected only from the spherical shell surrounding the nucleus. Simulated neurons were randomly placed in a volume of size 500 × 500×80 µm3 with no overlap between neurons up to the point that the volume was filled with cells or that neural density reached the value 83,100 ± 7,900 cells/mm3 (Keller et al., 2018). The resolution of the spatial volume was 0.5 µm/pixel in the x and y direction, 1 µm/pixel in the z direction.

Neural activity
Request a detailed protocol

Neural spiking activity was simulated as the sum of Poisson processes. Each neuron was assigned with a mean spiking rate (rho = 0.4, arbitrary selected), following a binary synchronicity matrix with value one for neurons with common inputs (common inputs probability = 0.8, arbitrary selected). The activity of each neuron was the sum of an independent Poisson process and as many common Poisson processes as the neurons with shared variability due to common inputs. The spiking rate of the summed Poisson processes was the mean spiking rate for that neuron. We assigned to randomly selected groups of neurons a shared common input generated as a Poisson process. This resulted in average pairwise correlations in the ground truth calcium activity larger than zero (average correlation 0.084 ± 0.001, n = 5833) and independent from the radial distance. Calcium activity and fluorescence traces were then generated using the equations in [Friedrich et al., 2017]. An autoregressive model of order 1 (with parameter γ = 0.7) was selected to convolve spike trains into calcium activity. A model with supralinearities and Hill saturation was used to convert calcium activity into fluorescence intensity (model parameters: baseline = 0; single event amplitude = 1500; Hill saturation coefficient = 1; dissociation constant = 1; noise variance = 0.05).

Generation of fluorescence time series
Request a detailed protocol

The size and the resolution of the simulated FOV were set to 500 × 500 µm2 and 2.5 µm/pixel, respectively. The resolution was adjusted according to the changes in the magnification factor (estimated from experimental data Figure 3—figure supplement 4), obtaining a non-uniform resolution in the FOV.

To generate the synthetic t-series, we used the measurements experimentally obtained from corrected and uncorrected microendoscopes (Figure 3). For corrected microendoscopes, the synthetic imaging focal surface was a spherical shell with curvature radius estimated from the measurements of fluorescent films (curvature radius: 400 μm) and the excitation volume was an ellipsoid resembling the aberration-corrected and experimentally measured PSF. For uncorrected microendoscopes, the synthetic imaging surfaces were two spherical shells with curvature radius estimated from the data (curvature radius: 265 μm and 2000 μm, respectively) and the excitation volume was an ellipsoid resembling the aberration-uncorrected and experimentally measured PSF.

Excitation volumes were scanned along the imaging focal surface (or surfaces for uncorrected microendoscopes), such that their axial direction was always orthogonal to the imaging focal surface(s). All the voxels falling within the excitation volumes contributed to the signal of the corresponding pixel in the FOV, resulting in one of the following possible three conditions:

  • If the pixel was in the edge of the FOV (radial distance >250 µm), its signal was randomly sampled from a normal distribution, with mean and standard deviation estimated from experimental data. For this and the following conditions, we selected the best fitting distribution for the signal mean through log-likelihood maximization across four alternative models: a normal, a gamma, a log-normal distribution, and a Gaussian mixture model. Dark noise mean was best fitted by a Gaussian mixture model (component 1: proportion = 0.37; mean = 137.48; sd = 48.96; component 2: proportion = 0.63; mean = 126.83; sd = 5.02, Log-likelihoods: Normal = −3.3E5, Gamma = −3.0E5, Log-normal = −2.9E5, Gaussian mixture model = −2.6E5), while the standard deviation of the dark noise depended on the dark noise mean in a linear way (p0 = −175.39, p1 = 1.57). The simulated dark noise was generated with the mean randomly sampled from the Gaussian Mixture Modeling (GMM) distribution and the standard deviation linearly dependent from the mean.

  • If the pixel was in the central part of the FOV (radial distance ≤250 µm) but no neurons were within the excitation volume, the pixel signal was randomly sampled from a normal distribution with mean and standard deviation estimated from experimental data. The mean intensity of pixels that were neither in the edges nor belonging to ROIs were fitted using a lognormal distribution (mean = 5.43, sd = 0.36. Log-likelihoods: Normal = −1.46E6, Gamma = −1.39E6, Log-normal = −1.37E6, Gaussian mixture model = −1.39E6) and the best linear fit between the squared root of the mean intensity and the intensity sd was computed (p0 = −162.55, p1 = 18.28). Simulated noise in the FOV was generated as Gaussian noise with mean randomly sampled from the lognormal distribution and sd linearly dependent from the squared root of the mean.

  • If the pixel was in the central part of the FOV (radial distance ≤250 µm) and at least one neuron was in the excitation volume(s), each voxel in the excitation volume(s) was assigned either Gaussian noise (estimated as in the previous condition. Log-likelihoods: Normal = −2.70E5, Gamma = −2.60E5, Log-normal = −2.58E5, Gaussian mixture model = −2.59E5) in case no neurons were in that voxel, or the fluorescence intensity of the neuron sampled by that voxel. In case a neuron was contained in a voxel, Gaussian noise was also added to the neuron signal. The mean of the added Gaussian noise was zero, while the sd was proportional to the square root of the mean intensity of the voxel, with the coefficients estimated from a linear fit between the square root of the mean intensity and the intensity standard deviation of pixels assigned to ROIs in experimental data (p0 = −132.44, p1 = 16.94). The activity of all the voxels falling within the excitation volume(s) was then averaged to obtain the pixel’s fluorescence intensity. The intensity of each pixel signal was finally modulated as a function of the radial position within the FOV, accordingly to the optical characterization of corrected and uncorrected microendoscopes using the radial intensity obtained imaging the subresolved fluorescent layer (Figure 3).

In simulations, the imaging rate of the t-series was set to 5 Hz.

Segmentation of simulated time series
Request a detailed protocol

We segmented simulated time series using two approaches: an automated procedure that we developed to resemble manual segmentation, and a standard automated approach (CaImAn [Giovannucci et al., 2019]). The automated procedure that we developed for the segmentation of synthetic t-series was based on the ground truth spatial distribution of neurons in the FOV. We first associated with each imaged neuron its spatial footprint, which consists of all the pixels collecting signal from that neuron. This segmentation would be the ideal segmentation. However, for experimental data, the ground truth is not available neither to users segmenting the FOV nor to automated segmentation algorithms. For this reason, we modified the ideal segmentation to obtain a more realistic situation. We reasoned that very small ROIs are not likely to be detected and therefore we removed all the ROIs with few pixels. We tried different thresholds for the minimum number of pixels composing a ROI (n = 5,10,15 pixels). We observed that this parameter did not have an effect on the results of the comparison between corrected and uncorrected microendoscopes and therefore we set it to 5. We then reasoned that overlapping ROIs could be distinguished only if their overlap was not perfect and therefore we merged ROIs with high overlap. The overlap between two ROIs was defined as the fraction between overlapping pixels and total number of pixels of the smallest ROI. We merged ROIs with overlap larger than 70, 80% and 90%. This parameter had no effect on the results of the comparison between corrected and uncorrected microendoscopes and we set it to 80%. Ultimately, we considered that ROIs whose fluorescence signal had low SNR could not be discriminated from noise. We reasoned that a single event was sufficient to segment a ROI and we defined a peak-SNR as the ratio between the maximum peak in the fluorescence trace and the baseline noise (defined as the sd of those parts of the fluorescence trace with signal lower than the 25% of the fluorescence distribution). We considered in the final segmentation only ROIs whose signal peak-SNR was higher than a threshold of 5, 10, 15, 20, 25, and 30 (Figure 4).

For the CaImAn segmentation, we tested different values of the SNR parameter (0.25, 0.5, 1, 1.5, 2), that represented a lower SNR threshold for a ROI to be kept in the final segmentation.

Statistics and analysis

Statistics

Values are expressed as mean ±sem, unless otherwise stated; the number of samples (n) and p values are reported in the Figure legends or in the text. No statistical methods were used to pre-determine sample size. All recordings with no technical issues were included in the analysis and blinding was not used in this study. Statistical analysis was performed with MATLAB software (Mathworks, Natick, MA) and GraphPad Prism software (GraphPad Software, San Diego, CA). A Kolmogorov-Smirnov test was run on each experimental sample to test for normality and to test the equality of the distributions in Figures 6g and 7m. The significance threshold was always set at 0.05. When comparing two paired populations of data, paired Student’s t-test or paired Wilcoxon signed-rank test (Figure 4c, Figure 4—figure supplement 1d,e) were used to calculate statistical significance in case of normal and non-normal distribution, respectively. Unpaired Student’s t-test (Figure 4k, Figure 4—figure supplement 2b,d) and Mann-Whitney test (Figures 4e,l and 7d,g) were used for unpaired comparisons of normally and non-normally distributed data, respectively. One-way analysis of variance (ANOVA) with Bonferroni post-hoc correction was used to compare the dependence of multiple classes from a single factor (Figures 6e,f and 7e). Two-way ANOVA with Bonferroni post-hoc correction was used to compare the dependence of multiple classes from two factors (Figure 6g). Two-way ANOVA with the interaction factor and with Bonferroni post-hoc correction was used in Figure 4d and Figure 4—figure supplement 2a. To fit linear data, linear regression was used. The significance of linear regression coefficients being different from zero was assessed using a permutation test, where we built a null distribution of linear regression coefficients by first destroying the relationship between the variable of interest (n = 5000 repetitions) and then fitting a linear regression model to each repetition of the permuted data (Figures 4g,m and 7h,i, and Figure 4—figure supplement 2c,e). Pearson correlation coefficients were used to test the dependence between variables (Figure 7c,f,k). The significance of the Pearson correlation coefficients was assessed using a Student’s t-test. All tests were two-sided, unless otherwise stated. Information theoretical analyses, NMF module identification, module characterization, and SVM classification were performed on MATLAB software (Mathworks, Natick, MA) using available toolboxes (Magri et al., 2009) or custom written codes.

Analysis of confocal images

Request a detailed protocol

Three mice were unilaterally injected with AAV-eGFP and red retrobeads after tissue aspiration and then implanted with the endoscope. Three mice were injected with AAV-GFP and red retrobeads without tissue aspiration and they were not implanted (controls). Three confocal images were acquired from fixed slices for each hemisphere at different focal planes (minimal distance between planes: 20 µm). Images in the red and green acquisition channels were blurred with a Gaussian filter (sigma = 2 µm) and binarized with a triangle thresholding method. S1bf was manually identified using anatomical cues from Hoechst labeling in each sample. To quantify the amount of preserved TC and CT connections within a given area of S1bf, we computed the fraction of pixels showing suprathreshold pixel intensity out of the total of pixel of the chosen area. A single value for each sample, obtained by averaging between confocal images of different FOVs, was used to run the one-tailed Mann-Whitney test for the different acquisition channels (Figure 5).

Analysis of field distortion and calibration of the pixel size

Request a detailed protocol

A regular fluorescent grid spanning the FOV was imaged in order to evaluate the distortion in the FOV. The number of pixels necessary to span 10 µm in the x and y direction was measured as a function of the distance from the FOV center. A magnification factor which varied along the radial directions was evaluated by computing the ratio between the measured number of pixels in the distorted (microscope objective coupled with GRIN-lens-based microendoscope) and undistorted (microscope objective alone) conditions. The estimated magnification factor (from x and y directions) was fitted using a quadratic curve (corrected: p0 = 0.76, p1 = −6.24E-04, p2 = 1.95E-05, norm of residuals = 1.78; Uncorrected: p0 = 0.73, p1 = −2.91E-04, p2 = 8.11e-06, norm of residuals = 0.24). The magnification factor was used to correctly calibrate experimental measurements in Figures 3, 4, 6 and 7, Figure 3—figure supplement 3, Figure 4—figure supplements 12, and Figure 7—figure supplement 1.

To measure the PSF as a function of the radial position within the FOV, z-stacks of subresolved fluorescent beads (diameter: 100 nm) were taken at different distances from the optical axis. Intensity profiles obtained from sections in the x, y, z directions of the PSFs were fitted with Gaussian curves and their FWHM was defined as x, y, and z resolution, respectively. Lateral resolution was calculated as the average of x and y resolution. Axial resolution coincided with the z resolution. When, due to aberrations in the lateral portion of the FOV, the intensity profile in the z direction was better fitted with a sum of two Gaussian curves instead of a single one, the axial resolution was defined as the axial distance between the two peaks of the best fitting curves. For each group of measurements at a specific distance from the optical axis, outliers were identified using the Rout method (GraphPad Software, San Diego, CA) and excluded from the data. Mean and standard deviation of resolutions were plotted against radial distance (Figure 3). Data were fitted with a symmetric quartic function to respect the cylindrical geometry of the optical system and the maximal FOV radial extent was determined as the radial distance at which the axial resolution fitting curve crossed a 10 µm threshold. PSF measurements were conducted using at least three eFOV-microendoscope for each type.

Analysis of fluorescence t-series

Request a detailed protocol

Experiments in VPM were analyzed with a customized graphical user interface (GUI) in MATLAB (version R2017a; Mathworks, Natick, MA). The GUI enabled the loading of data saved from the microscope acquisition software, included the motion correction algorithm (NoRMCorre) described in Pnevmatikakis and Giovannucci, 2017, facilitated the manual segmentation of the FOV, and allowed the deconvolution of neural activity from the recorded fluorescent dynamics. ROIs were drawn by visualizing single frames or temporal projections of the FOV and, within the selected regions, only those pixels necessary to maximize the peak-signal-to-noise ratio (SNR) of the mean intensity were selected. Peak-SNR was defined as:

peakSNR= maxt[tδt,t+δt]f(t)avg(fbaseline)std(fbaseline)

where fbaseline is the portion of intensity trace lower than the 25th percentile of the intensity distribution and t' is the peak time instant. After segmentation, the GUI performed the deconvolution of normalized calcium activity from the fluorescence extracted from ROIs, by using the algorithm provided in Giovannucci et al., 2019; Pnevmatikakis et al., 2016. The algorithm was based on the fit of the fluorescence activity with an autoregressive model. We used models of order 1 if the acquisition rate was low (<2 Hz), otherwise order 2. At the end of the pre-processing, a structure containing all the extracted information (acquisition parameters, ROIs spatial footprints, ROIs fluorescence activity, deconvolved activity, and normalized calcium activity) was saved and used for the subsequent analyses.

Analysis of synthetic and experimental time series

Request a detailed protocol
Quantification of ROIs number
Request a detailed protocol

For synthetic data in Figure 4 and Figure 4—figure supplement 2, the number of segmented ROIs was computed as a function of the probe type (uncorrected or corrected microendoscopes), of the peak-SNR value used in the segmentation procedure, and of the interaction between these two factor.

For experimental data in Figures 4, 6 and 7, background intensity was not always uniform in the endoscopes FOV and in some regions no ROIs could be detected. In order to discount this factor, the count of the ROIs was normalized to the brighter part of the FOV, obtaining a measurement of the ROIs density (number of ROIs divided by the total bright area). To detect the dark background regions, the edges of the FOV (where mostly dark noise was collected) were used as intensity threshold. All the parts of the FOV with mean intensity lower than 85th percentile of the threshold were discarded for the normalization of the ROIs count. The remaining part of the FOV was considered as bright area and used for the ROIs density analyses. For the benchmarking of manual and automatic segmentation in Figure 4—figure supplement 1, we computed recall, precision, and F1 score using the code provided in Soltanian-Zadeh et al., 2019.

SNR of calcium activity
Request a detailed protocol

For both synthetic and experimental data, the SNR was defined as:

SNR=power(Cdf)power(FrawCdf)=i(Cdf(ti))2i(Fraw(ti)Cdf(ti))2,

where Fraw and Cdf are the z-scored raw fluorescence intensity and deconvoluted calcium activity, respectively. SNR was evaluated to measure the quality of the extracted ROIs signal.

Correlation with ground truth activity
Request a detailed protocol

For synthetic data, we computed the correlation between the calcium activity of each segmented ROI and the ground truth calcium activity of neurons contributing to that ROI. In case more neurons were merged during the automated segmentation that we developed, we sorted the merged neurons for decreasing correlation with the corresponding ROI. We defined ‘source neuron’ the neuron with highest correlation with the ROI (Figure 4h left) and considered the correlation with the other merged neurons (only the second highest correlation is shown) as a measure of signals’ contamination between nearby neurons (Figure 4h right).

Pairwise correlation
Request a detailed protocol

For each pair of nearby ROIs (distance between ROIs centers < 20 μm), we computed pairwise correlation of the extracted calcium activity as a function of the radial distance of the ROIs pair. We defined the radial distance of each pair as the distance between the FOV’s center and the center of the segment connecting the two ROIs’ centers. To measure changes in pairwise correlations exclusively caused by changes in imaging resolution, we normalized the pairwise correlations of nearby ROIs by subtracting the value of pairwise correlations between distant ROIs (distance >60 μm) placed at the same radial distance. We then fitted the normalized pairwise correlation as a function of radial distance using linear regression.

Analysis of behavioral parameters

Request a detailed protocol

Videos of whisker movements were binarized with the Ridge Detection plugin of ImageJ/Fiji in order to individuate pixels corresponding to whiskers. Videos were then processed in MATLAB (Mathworks, Natick, MA) to extract the whisker mean angle. To this aim, all whiskers were fitted with straight lines and for each frame the mean angle of all the lines was calculated with respect to the horizontal direction of the FOV. Once the mean angle of the imaged whiskers was calculated for each frame, this signal was processed with a moving standard deviation over a 400 ms window and a Gaussian filter over a 50 ms window. Whisking and no whisking periods were identified by binarizing the mean whisker angle with a temporal and amplitude threshold. While the temporal threshold was fixed at 200 ms, the amplitude threshold was extracted by manually identifying whisking periods in ~1/10th of the full-length video and using this manual classification to find the best amplitude threshold with a ROC analysis. Temporal gaps between whisking periods shorter than 0.5 s were considered whisking periods, and linear interpolation was used to obtain the whiskers mean angle in frames in which less than four whiskers were detected.

For the analysis of pupil diameter, movies were analyzed with MATLAB (Mathworks, Natick, MA). Each frame was thresholded with the Otsu’s method and the region corresponding to the pupil was approximated with an ellipse. The length of the major axis of the ellipse was considered the pupil diameter. Linear interpolation was used for frames in which the pupil was not properly detected.

Detection of locomotion periods was performed using a threshold criterion on the wheel speed (Pakan et al., 2016). The wheel speed signal was downsampled at 40 Hz and an instant was considered to be part of a locomotion epoch if it met the following conditions: (i) instantaneous speed >1 cm/s; (ii) low-pass filtered speed (short pass filter at 0.25 Hz) >1 cm/s; (iii) average speed over 2 s windows >0.1 cm/s. Temporal gaps between locomotion periods shorter than 0.5 s were considered periods of locomotion.

Four behavioral states were defined in Figure 6: (i) quiet (Q) when neither locomotion nor whisking was observed, (ii) whisking (W), when whisking but no locomotion was observed; (iii) locomotion (L), when locomotion but not whisking was detected; (iv) whisking and locomotion (WL), when both locomotion and whisking were detected. L epochs were extremely rare (1.45 ± 0.75% of total acquisition time, mean ± sem) and were not considered in the analysis. For the SVM analysis and for the NMF analysis (Figures 6 and 7), we considered just two states: quiet (Q) and active (A), with A being the union of W and WL.

Analysis of calcium signals across states

Request a detailed protocol

To compare the amplitude of calcium activity across behavioral states, the deconvolved activity of each ROI was averaged in each of the three states (Q, W and WL).

To measure whether calcium activity was further modulated by arousal, we discretized the pupil size in ten bins and measured the distribution of the average calcium activity for each behavioral state, separately.

Information theoretic analysis

Request a detailed protocol

For information theoretical analyses, we used the MATLAB toolbox provided in Magri et al., 2009. To compute whisking information encoded in the calcium signal of single ROIs, we computed the amount of mutual information (Quian Quiroga and Panzeri, 2009) that the calcium signals carried about whether the animal was in a quiet (Q) or active (A) state. For each state, nT time points, where nT is the number of time points spent in the less frequent state, were randomly sampled without replacement. For each ROI then, the calcium activity in the selected time points was used to compute the information carried by that ROI regarding the state of the animal (the details were as follows: direct method, quadratic extrapolation for bias correction, number of response bins = 2). The amount of information was considered significantly different from 0 only if the real information value was larger than the 95th percentile of the distribution of n = 500 shuffled information values (obtained after randomly shuffling the whisking state label). This procedure was repeated n = 100 times, randomly sampling different time points and the reported information values were computed as the average information encoded across different iterations.

We sorted ROIs according to their information content and fitted the distribution of information content of individual ROIs using a double exponential function, where the information carried by the i-th ROI was given by:

IROIi=a*exp-b*i+c*exp(-d*i)

We used the R (Salinas and Sejnowski, 2001) coefficient to assess the goodness of the fit. We computed the Pearson correlation coefficient to check for dependence between the information carried by individual ROIs and the ROIs radial distance (or the ROIs SNR). To test whether ROIs at different radial distances (or with different SNR) carried a different amount of information, we split the ROIs in two groups. ROIs with low radial distance (or SNR) where those ROIs whose radial distance (or SNR) was lower than the median of the radial distances (or SNR) distribution. ROIs with high radial distance (or SNR) where those ROIs whose radial distance (or SNR) was higher (or equal) than the median of the radial distances (or SNR) distribution.

To compute information carried by pairs of neurons, we considered only pairs of nearby neurons (distance between neurons < 20 μm). We computed the amount of synergistic information carried by the neuronal pair as a function of the pairwise correlation between neurons and as a function of the pair’s radial distance. Data in Figure 7h were permutation subtracted. Specifically, information values were computed for n = 500 random permutations of the whisking variable label of the data, which destroyed the relationship between the whisking variable and the neuronal response. We then subtracted from the value of information obtained for real data the mean of the information obtained across permutations. We used the same permutations to compute a null distribution for synergy and redundancy values. n = 43 out of n = 61 pairs (70.4%) had significant synergy. n = 22 out of n = 31 pairs (70.1%) had significant redundancy. We fitted the synergistic information using a linear function for synergistic and redundant pairs, separately. We computed the linear fit using both all pairs (shown in Figure 7h) and only pairs showing significant synergy/redundancy (pairwise correlations: Synergy, slope = 0.009, permutation test p=0.2; Redundancy, slope = 0.05, permutation test p=0. Radial distance: Synergy, slope = −1E-4, permutation test p=0.002; Redundancy, slope = 7.1E-6, permutation test p=0.46). Results were not affected by this choice.

To compute information about the whisking state from a large population of neurons, we first decoded the whisking state from the single-trial population activity, and then we expressed the decoding performance as mutual information contained in the confusion matrix as in Equation (11) of Quian Quiroga and Panzeri, 2009. We performed single-trial population decoding as follows. A Gaussian kernel SVM was trained to classify the animal state (Q or A) by observing population activity (Chang and Lin, 2011). For each state, nT time points were randomly sampled (without replacement) and split into two equal and balanced sets. One set was used as training set for the SVM, and on this data, a ten-fold cross validation was performed over a fixed grid for the SVM hyperparameters. The performance of the SVM was then tested using the test set. This procedure was repeated n = 100 times by randomly sampling different time points. The reported classification accuracy was the average information encoded across different iterations. To check whether an increase in the number of imaged ROIs led to a better classification of whisking state (Figure 7e), we computed information for neuronal population of gradually increasing size. At first, we considered only the ROIs in the central portion of the FOV (distance from FOV center < ¼ of FOV radius) and then the other ROIs (distance steps = ¼ of FOV radius) were incrementally added for the training and testing phase of the SVM.

Non-negative matrix factorization (NMF)

Request a detailed protocol

To compute NMF (Figure 7), two states (Q and A) were considered. For each state, nT time points were randomly sampled (without replacement) and split into two equal and balanced sets (a training set and a test set). The number of NMF modules was selected based on the ability of a linear discriminant analysis (LDA) classifier trained on the reduced data to predict the presence or absence of whisking (Onken et al., 2016). The dimensionality of the training set and the test set were at first reduced to k, with k ranging from one to the number of ROIs of the dataset. Then, for each factorization, the LDA classifier was trained on the training set to predict the behavioral state variable, and its performance was tested on the test set. The final dimension for the NMF was selected as the number of modules at which the first elbow in the performance plot (performance increase <0.4%) was found. Then, the dimensionality of the entire dataset was reduced by computing the NMF with the selected number of modules. For each module in the obtained factorization, the following quantities were computed:

  • Sparseness: sparseness=jwji2jwj(i), where wj(i) denotes the contribution of the jth ROI to the ith module. Sparseness values close to 1 indicate that few ROIs contribute heavily to the module, while sparseness values close to 0 indicate that the contribution of ROIs to the modules is more homogeneous.

  • Whisking modulation index (WMI): WMI=meanactA-mean(actQ)meanactA+mean(actQ), where actA and actQ denote the activation coefficients in each behavioral state. WMI > 0 indicates that the module’s activity is increased during A state, while WMI < 0 indicates that the module’s activity is reduced during A state.

  • Spatial spread: we defined as spatial spread the shortest path that connected the ten ROIs with highest weights (or all the ROIs in case a module was composed by fewer ROIs).

To compare pairs of modules, we computed the following similarity measures (Figure 7m):

  • The Jaccard index is the fraction between the number of ROIs belonging to both modules and the total number of ROIs composing the two modules (without repetitions). It ranges between 0 and 1 and assumes value 0 if two modules do not share common ROIs, value one if two modules are composed by exactly the same ROIs. The Jaccard index does not take into account the weights of ROIs in the modules.

  • Cosine similarity = wi*viwi2*vi2, where wi and vi represent the weight of the i-th ROI in each of the two modules. As for the Jaccard index, cosine similarity ranges between 0 and 1 and it takes value 0 if two modules do not share common ROIs. Contrarily to the Jaccard index, cosine similarity takes into account also the weights of the ROIs in the modules and it is 1 only if two ROIs are composed by exactly the same ROIs with equal weights.

For both the Jaccard index and Cosine similarity, we computed null distributions. Specifically, for each NMF factorization we reassigned ROIs randomly within each module by shuffling their weights. We did not consider modules composed by single ROIs.

Data and software availability

Request a detailed protocol

The datasets shown in Figures 4, 6 and 7 and corresponding figure supplements are available at: https://doi.org/10.17632/wm6c5wzs4c.1.

The software used in this paper to generate artificial t-series is available at: https://github.com/moni90/eFOV_microendoscopes_sim (Antonini, 2020; copy archived at https://github.com/elifesciences-publications/eFOV_microendoscopes_sim).

Data availability

The datasets shown in Figures 1-7 and corresponding figure supplements are available at: https://doi.org/10.17632/wm6c5wzs4c.2. The software used in this paper to generate artificial t-series is available at: https://github.com/moni90/eFOV_microendoscopes_sim (copy archived at https://github.com/elifesciences-publications/eFOV_microendoscopes_sim) Numerical data for graphs represented in figures 3-7, figure 2-figure supplement 2-4, figure 3-figure supplement 4, figure 4-figure supplement 1-2 are provided as source data.

The following data sets were generated
    1. Antonini A
    2. Sattin A
    3. Moroni M
    4. Bovetti S
    5. Moretti C
    6. Succol F
    7. Forli A
    8. Vecchia D
    9. Rajamanickam VP
    10. Bertoncini A
    11. Panzeri S
    12. Liberale C
    13. Fellin T
    (2020) Mendeley Data
    Dataset of Extended field-of-view ultrathin microendoscopes for high-resolution two-photon imaging with minimal invasiveness.
    https://doi.org/10.17632/wm6c5wzs4c.2

References

  1. Book
    1. Smith WJ
    (2008)
    Modern Optical Engineering
    McGraw-Hill, NY.

Decision letter

  1. Ryohei Yasuda
    Reviewing Editor; Max Planck Florida Institute for Neuroscience, United States
  2. Andrew J King
    Senior Editor; University of Oxford, United Kingdom
  3. Darcy S Peterka
    Reviewer; Columbia University, United States
  4. Kaspar Podgorski
    Reviewer; Janelia Research Campus, United Kingdom

In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.

Acceptance summary:

In this revised manuscript, the authors demonstrated an extended field-of-view (eFOV), thin (<500μm) microendoscopes by adding aberration correcting microlenses fabricated by TPL. The performance of the microlens corrected eFOV endoscopes was evaluated with simulations on synthetic calcium data and in vivo activity imaging. In addition, the eFOV microendoscopes were used to investigate VPM activity correlated with locomotion, whisker movement, and pupil diameter. Cell-specific encoding in VPM during these behaviors is further scrutinized with statistical and machine learning methods. The presentation of the methods and the results are clear. The methods are useful and practical for increasing the FOV of ultrathin microendoscopes, removing one of the current limitations of small GRIN lens applications. The authors significantly improved the manuscript based on the reviewer's comments.

Decision letter after peer review:

Thank you for submitting your article "Extended field-of-view ultrathin microendoscopes for high-resolution two-photon imaging with minimal invasiveness" for consideration by eLife. Your article has been reviewed by three peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Andrew King as the Senior Editor. The following individuals involved in review of your submission have agreed to reveal their identity: Darcy S Peterka (Reviewer #2); Kaspar Podgorski (Reviewer #3).

The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.

We would like to draw your attention to changes in our revision policy that we have made in response to COVID-19 (https://elifesciences.org/articles/57162). Specifically, when editors judge that a submitted work as a whole belongs in eLife but that some conclusions require a modest amount of additional new data, as they do with your paper, we are asking that the manuscript be revised to either limit claims to those supported by data in hand, or to explicitly state that the relevant conclusions require additional supporting data.

Our expectation is that the authors will eventually carry out the additional experiments and report on how they affect the relevant conclusions either in a preprint on bioRxiv or medRxiv, or if appropriate, as a Research Advance in eLife, either of which would be linked to the original paper.

Your manuscript describes the development and use of microfabricated corrector lenses for improving GRIN lens imaging. These correction lenses increase the usable field-of-view (FOV) of a given GRIN lens, and improves the resolution. The approach was validated in a variety of samples. While the reviewers agree that the work will be potentially high impact, they feel that revisions with additional analyses would significantly improve the paper.

Reviewer #1:

In this manuscript, the authors demonstrated an extended field-of-view (eFOV), thin (<500μm) microendoscopes by adding aberration correcting microlenses fabricated bt TPL. Performance of the microlens corrected eFOV endoscopes was demonstrated with simulations on synthetic calcium data and in vivo activity imaging. In addition, the eFOV microendoscopes were used to investigate VPM activity correlated with locomotion, whisker movement, and pupil diameter. Cell-specific encoding in VPM during these behaviors are further scrutinized with statistical and machine learning methods. The presentation of the methods and the results are clear. The methods are useful and reasonably practical for increasing the FOV of ultrathin microendoscopes, which is a major limitation for small GRIN lenses. Overall, I recommend publication of the paper in eLife. Listed below are a number of concerns that we hope the authors can address in their revision to further improve the manuscript.

1) For the design and fabrication of the eFOV-microendoscope:

a) The authors presented four types of GRIN lenses (Type I-IV). Four different models of the GRIN lens are shown in subsection “Corrective lens manufacturing and microendoscope assembly”. It will probably help if the authors can explicitly indicate the Types with the model number.

b) It is claimed that the experimental result is similar to ray-trace simulation as shown in Figure 2. However, there is no comparison to support the claim. The high-order coefficient in the gradient profile in GRIN lens is usually not included in the optical model used for ray-trace. High-order aberration will thus not be reflected in the simulation. It would be important to know the difference between the simulation result and the experimental result. For example, the authors can add experimental data in Figure 2 or add simulation data in Figure 3.

c) For all the experiments, is it performed with the 2-photon polymerized lens with resin or with lens replica using a molded UV-curable adhesive? Is the aberration correction performance different from these two fabrication methods? Which material is the ray-tracing simulation based on?

d) What is the yield of the lens fabricated? How is the cured UV-curable adhesive detached from the PDMS mold? During imaging sessions, was there any damage on the lens?

e) How does the correction microlens affect the focal length of the GRIN lens. In other words, what is the allowable change in working distance before and after the correction?

f) There appears to be a much larger field curvature with the corrected GRIN lens. It would be helpful to discuss the impact of this on functional imaging.

g) Instead of AO with an active element (e.g., SLM), one could also put a fixed lens with the desired curvature in the place of the SLM. Compared to the microlens approach, one needs to modify the excitation path somewhat, but the lens can be fabricated using standard techniques and may be more accessible to most research labs. It might be informative to the readers to comment on this alternative approach.

2) For the data taken with the eFOV-microendoscope, there is ~2x improvement in FOV for type III lens (Figure 3G). However, in Figure 3—figure supplement 5C, the images before and after correction for type III lens is almost the same.

3) For the data simulation and analysis:

a) Analysis was shown to correlate whisking, locomotion and cell-specific encoding. Figure 6C and 6D show that locomotion is usually correlated with a higher ΔF/F. Could this be a result of the motion artifact during locomotion from the head-constraint mouse?

b) To simulate synthetic calcium data, different models (Gaussian Mixture, lognormal, etc.) are used to obtain the pixel intensity value as described in subsection “Generation of fluorescence time series”. It is mentioned that the model is obtained by fitting the experimental data. What criteria (e.g. comparing likelihood) is used to choose those specific models? What other distribution is tested to choose the best-fit model?

c) As mentioned in subsection “Analysis of field distortion and calibration of the pixel size”, the author used a gaussian profile to define axial resolution. It is probably better to fit with a Lorentzian function I(z)=〖(1/(1+〖(z/zR)〗^2 ))〗^2 for 2-photon excitation, where zR is the Rayleigh range of the gaussian focus. Not sure how large the difference will be between these two fits.

d) Figure 4D. It is not clear to me why the ROI number is smaller for the corrected lens than the uncorrected when the peak SNR value is small (e.g., less than 10). Can the authors add some explanations?

e) The authors claim that "pairwise correlation increased as a function of the radial distance of the pair in uncorrected probes compared to corrected ones." Figure 4M does not seem to show noticeable increase as a function of the radial distance for the uncorrected lens (left panel). Rather, it seems to show a noticeable decrease as a function of the radial distance for the corrected lens (right panel). Therefore, the statement does not seem to be supported by Figure 4M. This is confusing, and some explanation is needed.

4) There is a typo in the first sentence of subsection “Fabrication of eFOV-microendoscopes”. I believe it is (Figure 3—figure supplement 1C) instead of (Figure 1C).

Reviewer #2:

In this paper, the authors describe the creation and use of microfabricated corrector lenses for improving GRIN lens imaging. GRIN lenses are one of a few technologies (and the only one that is "mature") that allow for high resolution optical imaging in brains at depths beyond ~1.5 mm. These correction lenses a) increase the usable field-of-view (FOV) of a given GRIN lens, and b) improve the fidelity of signals extracted from regions-of-interest (ROIs), by decreasing the point spread function (PSF) volume. Both are important for high fidelity recording of neuronal populations. The paper is well written, and I find their demonstrations compelling. I think the work is of broad interest, and suitable for publication in eLife, after addressing a few concerns and questions.

My comments below are mainly critical, but I note that I think the paper overall is very good, and thorough, and I have high enthusiasm.

Given that the surface profiles begin their optimization from an analytical expression, It would be nice if the authors could comment on the likely maximal possible correction possible – that is on a GRIN lens with a GRIN relay of significant pitch, do they think microlenses would be similarly manufacturable or effective?

In the simulated imaging data, especially the figure, and SNR conclusions drawn from that, I think it is important to more clearly identify the results as "simulations". The authors do label and describe the data as such, however, because of the look of the figure, and in the Discussion, it is very easy to misinterpret as measured experimental data, rather than an in silico sim – perhaps panel shading on a-h would help?

Regarding the simulations, it was not clear what the "imaging" rate was set to be. I am curious why processes where not included in the in silico simulation. I would expect the corrected FOVs to have even better performance than uncorrected, if processes were included, and it would have been nice to explore. Further, I am not quite sure why the authors did not use a standard method to extract the putative ROIs, and compare that with their ground truth – it seems they considered the flaws of standard practice, and biased their sampling of ROIs directly – but I would think using a standard pipeline instead would have been more informative.

Along with this, in the real experiments, different methods were used to pick the ROIs – I am wondering if the authors could comment on the rational for the choices? I say this also, because while one would prefer to record more faithful signals directly, it isn't clear to me that some of the additional "mixing" that would occur between the ROIs would not be removed through source factorization.

The slow imaging rate also leads to higher correlations. Did the authors also collect signals at high frame rates? The basic optical properties would not change at higher rates, but 2-3 Hz is far from typical functional imaging rates. The gains in SNR may be more material at these higher frame rates too. Was it simply a microscope limitation that prohibited faster imaging?

For Figure 7, It is not clear to me why there should be a bias toward information away from the center assuming (semi) random placement of the GRIN- do the authors think that this is real, or simply a measurement artifact – the center of the FOV is nominally furthest from completely healthy tissue, so the slightly lower information reflects network damage? That would also affect the magnitude (but likely not the fundamental interpretation) of the spatial distribution of cells in the modules.

Reviewer #3:

In this paper, the authors describe a novel approach for improving two-photon imaging through thin GRIN-lens endoscopes by using a nanofabricated corrective lens. They validate the approach with in vitro measurements and simulations, and apply the technique to imaging a variety of samples, most notably somatosensory thalamus in awake, behaving mice.

The technique seems easy to apply and could advance the state of the art in several labs. Although I haven't done such experiments myself, I know from other groups that experiments are often sorely limited by their field of view. This paper presents an exciting solution to this problem.

I am enthusiastic about the work, but point out the following issues:

1) With uncorrected GRIN endoscopes, 3D imaging is possible. Do the corrective lenses affect the axial range that can be scanned this way, and if they do, how does this impact the total of neurons that can be recorded with a corrected vs uncorrected endoscope? Zemax simulations would be sufficient to address this.

2) I think the paper could do a better job of motivating the need for thinner endoscopes, particularly in the Introduction. Figure 5 leaves me wondering whether connections would really be disrupted by a larger endoscope. Figure 3—figure supplement 7 makes it clear that larger endoscopes and cannulae require removing a lot of structures, but the paper leaves it to the reader to imagine the consequences. Have there been published accounts of failure rates for different endoscope sizes? Perhaps your own experiences, or studies of e.g. amygdala or hypothalamus?

3) The methods used to segment/analyze calcium imaging data (subsection “Segmentation of simulated time series”) come off as somewhat straw-man. There are a variety of activity-based tools that can separate mixed signals from overlapping cells (PCA/ICA, NMF, CNMF, etc.), and these have been popular for (1P) endoscopic imaging. I wonder if such methods narrow the gap between uncorrected and corrected SNRs and correlations. The statements at the end of paragraph five of the Discussion seem a bit strong if this possibility has not been explored.

4) How chromatic are these corrections? The authors propose simultaneous functional imaging and optogenetic perturbations with these corrective lenses, but this would require the correction be useful across a large wavelength range (perhaps ~900nm and ~1040nm), which needs more evidence. Chromatic issues can easily affect 2P imaging even over the bandwidth of a femtosecond laser.

https://doi.org/10.7554/eLife.58882.sa1

Author response

Reviewer #1:

In this manuscript, the authors demonstrated an extended field-of-view (eFOV), thin (<500μm) microendoscopes by adding aberration correcting microlenses fabricated bt TPL. Performance of the microlens corrected eFOV endoscopes was demonstrated with simulations on synthetic calcium data and in vivo activity imaging. In addition, the eFOV microendoscopes were used to investigate VPM activity correlated with locomotion, whisker movement, and pupil diameter. Cell-specific encoding in VPM during these behaviors are further scrutinized with statistical and machine learning methods. The presentation of the methods and the results are clear. The methods are useful and reasonably practical for increasing the FOV of ultrathin microendoscopes, which is a major limitation for small GRIN lenses. Overall, I recommend publication of the paper in eLife. Listed below are a number of concerns that we hope the authors can address in their revision to further improve the manuscript.

1) For the design and fabrication of the eFOV-microendoscope:

a) The authors presented four types of GRIN lenses (Type I-IV). Four different models of the GRIN lens are shown in subsection “Corrective lens manufacturing and microendoscope assembly”. It will probably help if the authors can explicitly indicate the Types with the model number.

The correspondence between the type and the model number of the GRIN lenses used in the study is now indicated in paragraph two of subsection “Corrective lens manufacturing and microendoscope assembly”.

b) It is claimed that the experimental result is similar to ray-trace simulation as shown in Figure 2. However, there is no comparison to support the claim. The high-order coefficient in the gradient profile in GRIN lens is usually not included in the optical model used for ray-trace. High-order aberration will thus not be reflected in the simulation. It would be important to know the difference between the simulation result and the experimental result. For example, the authors can add experimental data in Figure 2 or add simulation data in Figure 3.

Following the reviewer’s comment, we ran simulations (now shown in new Figure 2—figure supplements 1-2) to quantify the excitation PSF for the four types of endoscopes and to compare these simulation results with the experimental findings displayed in Figure 3. Similarly to what observed in real experiments, in optical simulations we found that the axial dimension of the PSF in lateral portions of the FOV remained smaller and more similar to the axial dimension of the PSF in the center of the FOV in corrected compared to uncorrected endoscopes. Thus, these new optical simulations predicted enlarged FOV in corrected microendoscopes, a result that we experimentally validated (Figure 3). We edited the text in paragraph two of the Results and figures (new Figure 2—figure supplements 1-2) to include these new findings.

It must be noted that the absolute values of the axial PSF were generally smaller in the optical simulations compared to real measurements. This may be due to multiple reasons. First, as suggested by the reviewer, high-order aberrations were not included in the simulations. Second, in simulations, although the intensity of the excitation PSF was small in lateral portion of the FOV (Figure 2), a Gaussian function could still reliably fit the dim intensity distribution and provide a clear quantification of the PSF dimension. In experimental measurements of fluorescence emitted by subresolved beads, the more degraded PSF in the lateral portions of the FOV would result in low efficacy of the excitation beam in stimulating fluorescence, which would result in low SNR fluorescence signals. This would introduce large variability in the fit. Third, small variability in some of the experimental parameters (e.g., the distance between the GRIN back end and the focusing objective) were not reflected in the simulations.

c) For all the experiments, is it performed with the 2-photon polymerized lens with resin or with lens replica using a molded UV-curable adhesive? Is the aberration correction performance different from these two fabrication methods? Which material is the ray-tracing simulation based on?

Experiments and optical characterization were performed using lens replica only. Ray-trace simulation were performed considering the material used in lens replica (i.e., NOA63). We edited the text to make this clearer.

d) What is the yield of the lens fabricated? How is the cured UV-curable adhesive detached from the PDMS mold? During imaging sessions, was there any damage on the lens?

We thank the reviewer for their comments and we inserted the requested information in the text. Specifically:

“The yield for 3D printed lenses and lens replica was ~ 100 %.”

During molding, one drop of UV-curable adhesive was first put on a coverslip, which was then pressed against the mold. One side of the UV-curable adhesive was in contact with the mold, the other side was instead attached to the coverslip. After UV curing, by gently pulling the glass coverslip away the lens made of UV curable adhesive detached easily from the PDMS mold, while remaining firmly attached to the coverslip. This is now more clearly described in subsection “Corrective lens manufacturing and microendoscope assembly”.

We observed no appreciable damage on the lens over imaging sessions (now reported in subsection “Optical characterization of eFOV-microendoscopes”). In considering this, please note that laser beam is not directly focused on the polymer lens, but rather about 100 µm above it (Figure 1).

e) How does the correction microlens affect the focal length of the GRIN lens. In other words, what is the allowable change in working distance before and after the correction?

The simulated focal length in the absence and presence of the corrective microlens for the four different types of endoscopes is now reported in a novel table (new Supplementary file 1-table 2). The difference in focal length between uncorrected and corrected endoscopes is in the range 2-23 µm.

f) There appears to be a much larger field curvature with the corrected GRIN lens. It would be helpful to discuss the impact of this on functional imaging.

A short paragraph describing the impact of the larger field curvature of corrected endoscopes on functional imaging has been added in the Discussion. It reads: “Corrected endoscopes are characterized by a curved FOV. In the case of type II corrected endoscopes, the change in the z coordinate in the focal plane can be up to 75 µm (Figure 3). This z value is smaller for all other corrected endoscope types (Figure 3). The observed field curvature of corrected endoscopes may impact imaging in brain regions characterized by strong axially organized anatomy (e.g., the pyramidal layer of the hippocampus), but would not significantly affect imaging in regions with homogeneous cell density within the z range described above (< 75 µm for type II corrected microendoscopes).”

g) Instead of AO with an active element (e.g., SLM), one could also put a fixed lens with the desired curvature in the place of the SLM. Compared to the microlens approach, one needs to modify the excitation path somewhat, but the lens can be fabricated using standard techniques and may be more accessible to most research labs. It might be informative to the readers to comment on this alternative approach.

Following the reviewer’s comment, we modified the text to discuss this possibility. The new text reads: “A potential alternative to the approach describe in this study would be to place a macroscopic optical element of the desired profile in a plane optically conjugated to the objective back aperture along the optical path. This solution could have the advantage of being manufactured using more standard techniques. However, it would require significant change in the optical set-up, in contrast to the built-in correction method that we describe in the present study. Moreover, this macroscopic optical element would have to be changed according to the type of microendoscope used”.

2) For the data taken with the eFOV-microendoscope, there is ~2x improvement in FOV for type III lens (Figure 3g). However, in Figure 3—figure supplement 5C, the images before and after correction for type III lens is almost the same.

We substituted the images in Figure 3—figure supplement 5C with new ones that are more representative of the results shown in Figure 3G.

3) For the data simulation and analysis:

a) Analysis was shown to correlate whisking, locomotion and cell-specific encoding. Figure 6C and 6D show that locomotion is usually correlated with a higher ΔF/F. Could this be a result of the motion artifact during locomotion from the head-constraint mouse?

We think it is unlikely that the observed increase in ΔF/F is due to motion artifacts during locomotion because: i) we implemented a motion correction procedure in our data analysis; ii) the effect of potential motion artefacts on ΔF/F values would likely be that of increasing dispersion rather than biasing the mean towards higher values, as observed in Figure 6F. This is because a change in focus or a x,y translation due to a motion artefact would likely have similar probability of increasing or decreasing ΔF/F values.

b) To simulate synthetic calcium data, different models (Gaussian Mixture, lognormal, etc.) are used to obtain the pixel intensity value as described in subsection “Generation of fluorescence time series”. It is mentioned that the model is obtained by fitting the experimental data. What criteria (e.g. comparing likelihood) is used to choose those specific models? What other distribution is tested to choose the best-fit model?

We fitted the experimental data using four alternative models: a normal, a γ, a log-normal distribution, and a Gaussian mixture model. We selected the best-fit model by maximizing the log-likelihood.

In the following tables, the log-likelihood values of the used distributions pixels in the center and lateral portion of the FOV are reported for the reviewer’s inspection (highlighted in green the one we used under each condition):

Central pixels (without ROIs)

Central pixels (with ROIs)

Edges” pixels (“dark noise”):

c) As mentioned in subsection “Analysis of field distortion and calibration of the pixel size”, the author used a gaussian profile to define axial resolution. It is probably better to fit with a Lorentzian function I(z)=〖(1/(1+〖(z/zR)〗^2 ))〗^2 for 2-photon excitation, where zR is the Rayleigh range of the gaussian focus. Not sure how large the difference will be between these two fits.

Following the reviewer’s comment, we used the Lorentzian function to fit the fluorescence intensity distribution for type II microendoscopes. Axial resolution values computed using the fit with the Lorentzian function were not significantly different compared to those obtained fitting with a Gaussian function (Wilcoxon matched-pairs signed rank test, p = 0,93, n = 132 PSF measurements). We therefore decided to leave the results obtained with the Gaussian function in the manuscript.

d) Figure 4D. It is not clear to me why the ROI number is smaller for the corrected lens than the uncorrected when the peak SNR value is small (e.g., less than 10). Can the authors add some explanations?

We thank the reviewer for raising this point. The smaller number of ROIs when the SNR is low in corrected endoscopes is due to:

1) In the segmentation method we implemented, which is based on the ground truth distribution of the neurons in the simulated sample, at least two pixels belonging to a ground truth neuron are defined as a ROI.

2) In uncorrected endoscopes, the axial PSF largely increases as a function of the radial distance. The enlarged axial PSF in the lateral portions of the FOV augments the probability of sampling voxels belonging to multiple neurons located at different z positions. Once projected in the 2D plane, the contribution of multiple neurons located at different z positions increases the probability of having pixels belonging to ROIs. An increased axial PSF thus leads to an increased number of detected ROIs.

3) Corrected endoscopes have smaller axial PSF compared to uncorrected ones and thus smaller number of detected ROIs.

We added a short text to better explain this result in subsection “Higher SNR and more precise evaluation of pairwise correlation in eFOV-microendoscopes”.

e) The authors claim that "pairwise correlation increased as a function of the radial distance of the pair in uncorrected probes compared to corrected ones." Figure 4M does not seem to show noticeable increase as a function of the radial distance for the uncorrected lens (left panel). Rather, it seems to show a noticeable decrease as a function of the radial distance for the corrected lens (right panel). Therefore, the statement does not seem to be supported by Figure 4M. This is confusing, and some explanation is needed.

We thank the reviewer for raising this point and we apologize with them if the indicated sentence was imprecise and generated confusion. What we meant was that the linear fit of pairwise correlations as a function of radial position for uncorrected endoscopes had a significantly positive slope (Figure 4M, left panel, slope = 0.0002, permutation test p = 0.006), indicating higher pairwise correlations in lateral compared to more central portions of the FOV. For corrected endoscopes, the slope of the linear fit was not significantly different from zero (Figure 4M, right panel, slope = -0.0005, permutation test p = 0.05 for dataset 1; slope = -0.0007, permutation test p = 0.05 for dataset 2). We modified the text to correct the previous imprecision and clarify this point.

4) There is a typo in the first sentence of subsection “Fabrication of eFOV-microendoscopes”. I believe it is (Figure 3—figure supplement 1C) instead of (Figure 1C).

Thank you for spotting this error. It has been fixed.

Reviewer #2:

In this paper, the authors describe the creation and use of microfabricated corrector lenses for improving GRIN lens imaging. GRIN lenses are one of a few technologies (and the only one that is "mature") that allow for high resolution optical imaging in brains at depths beyond ~1.5 mm. These correction lenses a) increase the usable field-of-view (FOV) of a given GRIN lens, and b) improve the fidelity of signals extracted from regions-of-interest (ROIs), by decreasing the point spread function (PSF) volume. Both are important for high fidelity recording of neuronal populations. The paper is well written, and I find their demonstrations compelling. I think the work is of broad interest, and suitable for publication in eLife, after addressing a few concerns and questions.

My comments below are mainly critical, but I note that I think the paper overall is very good, and thorough, and I have high enthusiasm.

Given that the surface profiles begin their optimization from an analytical expression, It would be nice if the authors could comment on the likely maximal possible correction possible – that is on a GRIN lens with a GRIN relay of significant pitch, do they think microlenses would be similarly manufacturable or effective?

We thank the reviewer for this interesting question. It is expected that optical aberrations of GRIN-rod lenses become more pronounced as the GRIN-rod length increases. This raises the possibility that the needed correcting lens profile could significantly differ from our current initial guess (the Schmidt corrector plate profile). To properly answer the reviewer’s question, there are thus a number of issues that need to be evaluated altogether, namely: a) up to what length of the GRIN-rod, for a given NA, the Schmidt corrector plate profile is an effective starting guess for the lens optimization process; b) once such an optimized correcting profile is found, to what extent the introduced correction is able to increase the effective FOV; c) how accurately the complex correcting lens profile could be reproduced by 3D micro-printing. Because our present approach satisfactory works for the GRIN-rod lengths considered in our study (1.1 mm – 4.07 mm), and because of the amount of work needed to address the manifold and interrelated aspects of the asked question, we deemed it to be the subject of a whole new study. As consequence, we currently cannot provide a final answer to this question, but we aim at addressing the matter in the continuation of our work. Please note that in the paper we do not make any claim that the method is applicable to GRIN lenses of any length, and we made sure in revision that this limitation in scope is acknowledged.

In the simulated imaging data, especially the figure, and SNR conclusions drawn from that, I think it is important to more clearly identify the results as "simulations". The authors do label and describe the data as such, however, because of the look of the figure, and in the Discussion, it is very easy to misinterpret as measured experimental data, rather than an in silico sim – perhaps panel shading on a-h would help?

We thank the reviewer for their suggestion. We followed the indication of the reviewer and modified Figure 4 to make clear, using boxes, which data come from simulations and which from real experiments.

Regarding the simulations, it was not clear what the "imaging" rate was set to be.

In the simulations the “imaging” rate was set to 5 Hz. This is indicated in subsection “Generation of fluorescence time series”.

I am curious why processes where not included in the in silico simulation. I would expect the corrected FOVs to have even better performance than uncorrected, if processes were included, and it would have been nice to explore.

We did not address the impact of including processes in the simulations, because the axial PSF of the GRIN microendoscopes was in the range 8-10 µm and the labeling that we achieved in our preparation was rather dense. Under these experimental conditions, resolving individual neuronal processes was complicated and we limited our goal to image fluorescent signals from neuronal cell bodies. We agree with the reviewer this is an interesting point, which can be addressed in future investigations.

Further, I am not quite sure why the authors did not use a standard method to extract the putative ROIs, and compare that with their ground truth – it seems they considered the flaws of standard practice, and biased their sampling of ROIs directly – but I would think using a standard pipeline instead would have been more informative.

Following the reviewer’s comment, we first compared the quality of the manual segmentation described in this study with that of a standard automated algorithm (e.g., CaImAn, (Giovannucci et al., 2019) by computing precision, recall, and F1 score in simulated data (new Figure 4—figure supplement 1). We found that the automated method was characterized by recall values which were, on average across SNR values, < 0.4 (new Figure 4—figure supplement 1A), leading to the detection of only a minority of the total number of neurons. In contrast, the manual method led to higher recall across SNR threshold values. Notably, both manual and automated methods yielded larger recall in corrected endoscopes compared to uncorrected ones. Moreover, for low SNR threshold values the automated segmentation had precision values < 0.8 in both uncorrected and corrected endoscopes (new Figure 4—figure supplement 1B), leading to identification of ROIs which did not correspond to cells in the ground truth. In contrast, the manual segmentation method had much larger values of precision across SNR threshold levels. Overall, F1 scores were higher for the manual segmentation method compared to the automated one for both uncorrected and corrected endoscopes (new Figure 4—figure supplement 1C). Because of these results showing better performance of the manual segmentations in simulated data (an observation now included in the revised manuscript), we took the decision to present in the main figures of the paper results based on manual segmentation.

In the revised manuscript, we extended the comparison between the manual and automated segmentation methods to real data. We observed that in uncorrected endoscopes CaImAn identified smaller number of ROIs compared to manual segmentation (new Figure 4—figure supplement 1D). In contrast, the number of ROIs identified with CaImAn and the manual method in t-series acquired with the corrected endoscope were not significantly different (new Figure 4—figure supplement 1E). One potential explanation of this finding is that the automated segmentation method more efficiently segments ROIs with high SNR compared to the manual one. Since aberration correction significantly increases SNR of fluorescent signals, the automated segmentation performed as the manual segmentation method in corrected endoscopes.

Overall, we believe that, in the context of our study, the use of the manual segmentation method rather than an automated segmentation algorithm is still preferable, because it allowed a more precise evaluation of ROIs in simulated data and fairer comparison of corrected and uncorrected endoscopes in real data. Most importantly, improvements introduced by corrected endoscopes could be observed with both the manual and the automated segmentation methods (see response to the next point). We believe that the comparisons suggested by reviewer are of interest and we feel that the new results in Figure 4—figure supplement 1 improve our manuscript. We are grateful to the reviewer for this suggestion.

Along with this, in the real experiments, different methods were used to pick the ROIs – I am wondering if the authors could comment on the rational for the choices? I say this also, because while one would prefer to record more faithful signals directly, it isn't clear to me that some of the additional "mixing" that would occur between the ROIs would not be removed through source factorization.

In all experiments presented in the main paper, we used our customized manual pipeline for the analysis of fluorescence t-series. In Figure 3—figure supplement 6 we used a different method, the PCA/ICA algorithm described in Mukamel et al., Neuron 2009, to segment neurons and extract calcium signals. This was originally done to show that data recorded with corrected endoscopes could be analyzed with standard automated methods. Following the comments of this reviewer and that of reviewer #3 asking for a systematic evaluation of how much the improvements introduced by aberration correction depend on the analysis pipeline, we removed the PCA/ICA analysis in Figure 3—figure supplement 6 and we systematically evaluated the effect of aberration correction on the output of the analysis in the simulated and experimental data shown in Figure 4 using both the manual segmentation method and an automated segmentation method (e.g., CaImAn).

In simulated data, we found that with CaImAn the number of ROIs segmented in corrected endoscopes was consistently higher than in uncorrected endoscopes across SNR thresholds (new Figure 4—figure supplement 2A), similarly to what observed with manual segmentation (Figure 4D). Using CaImAn, SNR values of fluorescence events were significantly higher in corrected compared to uncorrected endoscopes (new Figure 4—figure supplement 2B), similarly to what observed with manual segmentation (Figure 4E). Moreover, the linear fit of pairwise correlations as a function of radial position for uncorrected endoscopes had a significantly positive slope (new Figure 4—figure supplement 2C left panel), indicating higher pairwise correlations in lateral compared to more central portions of the FOV. For corrected endoscopes, the slope of the linear fit was not significantly different from zero (new Figure 4—figure supplement 2C right panel). This result is also in line with what we previously observed with manual segmentation (Figure 4G). The description of these novel results is added in the text.

In experimental data, we found that SNR values of fluorescence events tended to be higher in corrected compared to uncorrected endoscopes (new Figure 4—figure supplement 2D), a trend that was in line with what observed with manual segmentation (Figure 4L). The slope of the linear fit of pairwise correlations as a function of radial position for uncorrected endoscopes was significantly positive (new Figure 4—figure supplement 2E), indicating higher pairwise correlations in lateral compared to more central portions of the FOV. For corrected endoscopes, the slope of the linear fit was not significantly different from zero (new Figure 4—figure supplement 2E). Both results are in line with what previously observed with manual segmentation (Figure 4M). The description of these novel results is added in the text.

Overall, the results of the systematic comparison between the manual and automated segmentation methods show that improvements introduced by aberration correction in endoscopes could be observed with both the manual and the automated segmentation methods. We thank the reviewer for raising this point.

The slow imaging rate also leads to higher correlations. Did the authors also collect signals at high frame rates? The basic optical properties would not change at higher rates, but 2-3 Hz is far from typical functional imaging rates. The gains in SNR may be more material at these higher frame rates too. Was it simply a microscope limitation that prohibited faster imaging?

We did not collect fluorescent signals at frame rates higher than 4 Hz. This was because we aimed at imaging the largest possible FOV and our experimental set up was equipped with regular galvanometric mirrors and not with resonant mirrors. We added a sentence in the text to highlight this limitation of our study.

For Figure 7, It is not clear to me why there should be a bias toward information away from the center assuming (semi) random placement of the GRIN- do the authors think that this is real, or simply a measurement artifact – the center of the FOV is nominally furthest from completely healthy tissue, so the slightly lower information reflects network damage? That would also affect the magnitude (but likely not the fundamental interpretation) of the spatial distribution of cells in the modules.

At the single cell level, information was evenly distributed across the FOV (Figure 7C). We believe the reviewer’s comment refers to data shown in Figure 7E, which contains data pooled from a population of neurons. In Figure 7E, we observed an increase in the information carried by the population of neurons when we incrementally considered larger populations by adding ROIs at larger radial distances. This does not imply that ROIs away from the center carry higher information, but only that their information sums to the information of more central ROIs in a way that is not purely redundant. We apologize if the description of Figure 7E was confusing, and we edited the text to make it clearer.

Reviewer #3:

In this paper, the authors describe a novel approach for improving two-photon imaging through thin GRIN-lens endoscopes by using a nanofabricated corrective lens. They validate the approach with in vitro measurements and simulations, and apply the technique to imaging a variety of samples, most notably somatosensory thalamus in awake, behaving mice.

The technique seems easy to apply and could advance the state of the art in several labs. Although I haven't done such experiments myself, I know from other groups that experiments are often sorely limited by their field of view. This paper presents an exciting solution to this problem.

I am enthusiastic about the work, but point out the following issues:

1) With uncorrected GRIN endoscopes, 3D imaging is possible. Do the corrective lenses affect the axial range that can be scanned this way, and if they do, how does this impact the total of neurons that can be recorded with a corrected vs uncorrected endoscope? Zemax simulations would be sufficient to address this.

The optical simulations we did to design the corrective lenses were performed maximizing aberration correction only in the focal plane of the endoscope. Following the reviewer’s comment, we explored the effect of aberration correction outside the focal plane using new Zemax simulations. In corrected endoscopes, we found that the Strehl ratio was > 0.8 (Maréchal criterion) in a larger volume compared to uncorrected endoscopes (new Figure 2—figure supplement 3). The range of volume with high Strehl ratio in corrected endoscopes was 1.7-3.7 larger than the volume with high Strehl ratio in uncorrected probes, leading to similarly larger number of neurons being imaged in the enlarged volume. This is now described on lines (paragraph three subsection “Optical simulation of eFOV-microendoscopes”).

2) I think the paper could do a better job of motivating the need for thinner endoscopes, particularly in the Introduction. Figure 5 leaves me wondering whether connections would really be disrupted by a larger endoscope. Figure 3—figure supplement 7 makes it clear that larger endoscopes and cannulae require removing a lot of structures, but the paper leaves it to the reader to imagine the consequences. Have there been published accounts of failure rates for different endoscope sizes? Perhaps your own experiences, or studies of e.g. amygdala or hypothalamus?

Following the reviewer’s comment, we modified the Introduction to more clearly motivate the need of thin endoscopic probes for deep brain functional imaging in paragraph two of the Introduction. We are not aware of a systematic description on the impact of different endoscope sizes on the functionality of the targeted brain region, but we have cited relevant literature showing that surgical lesion of reciprocal connectivity alters thalamocortical dynamics.

3) The methods used to segment/analyze calcium imaging data (subsection “Segmentation of simulated time series”) come off as somewhat straw-man. There are a variety of activity-based tools that can separate mixed signals from overlapping cells (PCA/ICA, NMF, CNMF, etc.), and these have been popular for (1P) endoscopic imaging. I wonder if such methods narrow the gap between uncorrected and corrected SNRs and correlations. The statements at the end of paragraph five of the Discussion seem a bit strong if this possibility has not been explored.

Following the reviewer’s comment we evaluated the effect of aberration correction on the output of the analysis in the simulated and experimental data shown in Figure 4 using both the manual segmentation method and an automated segmentation method (e.g., CaImAn).

In simulated data, we found that with CaImAn the number of ROIs segmented in corrected endoscopes was consistently higher than in uncorrected endoscopes across SNR thresholds (new Figure 4—figure supplement 2A), similarly to what observed with manual segmentation (Figure 4D). Using CaImAn, SNR values of fluorescence events were significantly higher in corrected compared to uncorrected endoscopes (new Figure 4—figure supplement 2B), similarly to what observed with manual segmentation (Figure 4E). Moreover, the linear fit of pairwise correlations as a function of radial position for uncorrected endoscopes had a significantly positive slope (new Figure 4—figure supplement 2C left panel), indicating higher pairwise correlations in lateral compared to more central portions of the FOV. For corrected endoscopes, the slope of the linear fit was not significantly different from zero (new Figure 4—figure supplement 2C right panel). This result is also in line with what previously observed with manual segmentation (Figure 4G). The description of these novel results is added in subsection “Higher SNR and more precise evaluation of pairwise correlation in eFOV-microendoscopes”.

In experimental data, we found that SNR values of fluorescence events tended to be higher in corrected compared to uncorrected endoscopes (new Figure 4—figure supplement 2D), a trend that was in line with what observed with manual segmentation (Figure 4L). The slope of the linear fit of pairwise correlations as a function of radial position for uncorrected endoscopes was significantly positive (new Figure 4—figure supplement 2E), indicating higher pairwise correlations in lateral compared to more central portions of the FOV. For corrected endoscopes, the slope of the linear fit was not significantly different from zero (new Figure 4—figure supplement 2E). Both results are in line with what previously observed with manual segmentation (Figure 4M).

The description of these novel results is added in subsection “Higher SNR and more precise evaluation of pairwise correlation in eFOV-microendoscopes”.

Overall, the results of the systematic comparison between the manual and automated segmentation methods show that improvement introduced by aberration correction in endoscopes could be observed with both the manual and automated segmentation methods. We thank the reviewer for raising this point.

4) How chromatic are these corrections? The authors propose simultaneous functional imaging and optogenetic perturbations with these corrective lenses (Discussion final paragraph), but this would require the correction be useful across a large wavelength range (perhaps ~900nm and ~1040nm), which needs more evidence. Chromatic issues can easily affect 2P imaging even over the bandwidth of a femtosecond laser.

Thanks for this comment. All optical simulations described in the first submission were performed at a fixed wavelength (λ = 920 nm). Following the reviewer’s request, we explored the effect of changing wavelength on the Strehl ratio using new Zemax simulations. We found that the Strehl ratio remains > 0.8 within ± 15 nm from λ = 920 nm (new Figure 2—figure supplement 4), which covers the limited bandwidth of our femtosecond laser but not that needed for all-optical imaging and manipulation experiments. To perform this latter type of experiments, the SLM that is used in the stimulation pathway should also be used to correct for the z-defocus with an appropriate lens function. Alternatively, future microfabrication work will be needed to design more complex corrective lens design, which efficiently compensates chromatic aberrations in the 900-1100 nm wavelength range. These possibilities are now more clearly discussed (Discussion paragraphs eight and nine).

https://doi.org/10.7554/eLife.58882.sa2

Article and author information

Author details

  1. Andrea Antonini

    1. Optical Approaches to Brain Function Laboratory, Istituto Italiano di Tecnologia, Genova, Italy
    2. Nanostructures Department, Istituto Italiano di Tecnologia, Genova, Italy
    Contribution
    Data curation, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing - original draft
    Contributed equally with
    Andrea Sattin
    Competing interests
    No competing interests declared
  2. Andrea Sattin

    1. Optical Approaches to Brain Function Laboratory, Istituto Italiano di Tecnologia, Genova, Italy
    2. University of Genova, Genova, Italy
    3. Neural Coding Laboratory, Istituto Italiano di Tecnologia, Genova and Rovereto, Italy
    Contribution
    Data curation, Software, Formal analysis, Validation, Investigation, Visualization, Writing - original draft
    Contributed equally with
    Andrea Antonini
    Competing interests
    No competing interests declared
  3. Monica Moroni

    1. Neural Coding Laboratory, Istituto Italiano di Tecnologia, Genova and Rovereto, Italy
    2. Neural Computation Laboratory, Center for Neuroscience and Cognitive Systems @UniTn, Istituto Italiano di Tecnologia, Rovereto, Italy
    3. Center for Mind and Brain Sciences (CIMeC), University of Trento, Rovereto, Italy
    Contribution
    Data curation, Software, Formal analysis, Validation, Investigation, Visualization, Writing - original draft
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-1852-7217
  4. Serena Bovetti

    Optical Approaches to Brain Function Laboratory, Istituto Italiano di Tecnologia, Genova, Italy
    Contribution
    Data curation, Software, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing - original draft
    Competing interests
    No competing interests declared
  5. Claudio Moretti

    1. Optical Approaches to Brain Function Laboratory, Istituto Italiano di Tecnologia, Genova, Italy
    2. University of Genova, Genova, Italy
    Contribution
    Data curation, Software, Formal analysis, Validation, Investigation, Visualization, Methodology
    Competing interests
    No competing interests declared
  6. Francesca Succol

    Optical Approaches to Brain Function Laboratory, Istituto Italiano di Tecnologia, Genova, Italy
    Contribution
    Validation, Investigation, Methodology
    Competing interests
    No competing interests declared
  7. Angelo Forli

    1. Optical Approaches to Brain Function Laboratory, Istituto Italiano di Tecnologia, Genova, Italy
    2. Neural Coding Laboratory, Istituto Italiano di Tecnologia, Genova and Rovereto, Italy
    Contribution
    Formal analysis, Supervision, Validation
    Competing interests
    No competing interests declared
  8. Dania Vecchia

    1. Optical Approaches to Brain Function Laboratory, Istituto Italiano di Tecnologia, Genova, Italy
    2. Neural Coding Laboratory, Istituto Italiano di Tecnologia, Genova and Rovereto, Italy
    Contribution
    Formal analysis, Visualization, Writing - original draft
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-6091-538X
  9. Vijayakumar P Rajamanickam

    1. Optical Approaches to Brain Function Laboratory, Istituto Italiano di Tecnologia, Genova, Italy
    2. Nanostructures Department, Istituto Italiano di Tecnologia, Genova, Italy
    3. Biological and Environmental Science and Engineering Division (BESE), King Abdullah University of Science and Technology (KAUST), Thuwal, Saudi Arabia
    Contribution
    Validation, Methodology
    Competing interests
    No competing interests declared
  10. Andrea Bertoncini

    Biological and Environmental Science and Engineering Division (BESE), King Abdullah University of Science and Technology (KAUST), Thuwal, Saudi Arabia
    Contribution
    Methodology
    Competing interests
    No competing interests declared
  11. Stefano Panzeri

    1. Neural Coding Laboratory, Istituto Italiano di Tecnologia, Genova and Rovereto, Italy
    2. Neural Computation Laboratory, Center for Neuroscience and Cognitive Systems @UniTn, Istituto Italiano di Tecnologia, Rovereto, Italy
    Contribution
    Formal analysis, Supervision, Methodology, Writing - original draft, Writing - review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-1700-8909
  12. Carlo Liberale

    1. Nanostructures Department, Istituto Italiano di Tecnologia, Genova, Italy
    2. Biological and Environmental Science and Engineering Division (BESE), King Abdullah University of Science and Technology (KAUST), Thuwal, Saudi Arabia
    Contribution
    Conceptualization, Supervision, Funding acquisition, Methodology, Writing - original draft, Project administration, Writing - review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-5653-199X
  13. Tommaso Fellin

    1. Optical Approaches to Brain Function Laboratory, Istituto Italiano di Tecnologia, Genova, Italy
    2. Neural Coding Laboratory, Istituto Italiano di Tecnologia, Genova and Rovereto, Italy
    Contribution
    Conceptualization, Resources, Supervision, Funding acquisition, Writing - original draft, Project administration, Writing - review and editing
    For correspondence
    tommaso.fellin@iit.it
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-2718-7533

Funding

European Research Council (NEURO-PATTERNS)

  • Tommaso Fellin

National Institutes of Health (BRAIN Initiative NS090576)

  • Tommaso Fellin

National Institutes of Health (BRAIN Initiative NS107464)

  • Stefano Panzeri
  • Tommaso Fellin

National Institutes of Health (BRAIN Initiative NS109961)

  • Stefano Panzeri

Seventh Framework Programme (DESIRE)

  • Tommaso Fellin

FIRB (RBAP11X42L)

  • Tommaso Fellin

Flag-Era JTC Human Brain Project (SLOW-DYN)

  • Stefano Panzeri
  • Tommaso Fellin

IIT interdisciplinary grant

  • Carlo Liberale
  • Tommaso Fellin

King Abdullah University of Science and Technology (BAS/1/1064-01-01)

  • Carlo Liberale

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank M Dal Maschio for discussion at an initial stage of the project, F Nespoli for preliminary analysis, and B Sabatini and A Begue for critical reading an early version of this manuscript. We thank DS Kim and the GENIE project for the constructs Addgene viral prep # 100845-AAV1 and #104492-AAV1, H Zeng for the construct Addgene viral prep # 51502-AAV1, and JM Wilson for the construct Addgene viral prep # 105558-AAV1. This work was supported by an IIT interdisciplinary grant and in part by ERC (NEURO-PATTERNS), NIH Brain Initiative (U01 NS090576, U19 NS107464, R01NS109961), FP7 (DESIRE), MIUR FIRB (RBAP11 × 42L), and Flag-Era JTC Human Brain Project (SLOW-DYN). CL acknowledges support from KAUST under baseline funding BAS/1/1064-01-01.

Ethics

Animal experimentation: Experimental procedures involving animals have been approved by the Istituto Italiano di Tecnologia Animal Health Regulatory Committee, by the National Council on Animal Care of the Italian Ministry of Health (authorization # 1134/2015-PR, # 689/2018-PR) and carried out according to the National legislation (D.Lgs. 26/2014) and to the legislation of the European Communities Council Directive (European Directive 2010/63/EU).

Senior Editor

  1. Andrew J King, University of Oxford, United Kingdom

Reviewing Editor

  1. Ryohei Yasuda, Max Planck Florida Institute for Neuroscience, United States

Reviewers

  1. Darcy S Peterka, Columbia University, United States
  2. Kaspar Podgorski, Janelia Research Campus, United Kingdom

Publication history

  1. Received: May 13, 2020
  2. Accepted: October 12, 2020
  3. Accepted Manuscript published: October 13, 2020 (version 1)
  4. Accepted Manuscript updated: October 15, 2020 (version 2)
  5. Version of Record published: November 24, 2020 (version 3)

Copyright

© 2020, Antonini et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 3,097
    Page views
  • 513
    Downloads
  • 18
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

Download links

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

Downloads (link to download the article as PDF)

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

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

  1. Andrea Antonini
  2. Andrea Sattin
  3. Monica Moroni
  4. Serena Bovetti
  5. Claudio Moretti
  6. Francesca Succol
  7. Angelo Forli
  8. Dania Vecchia
  9. Vijayakumar P Rajamanickam
  10. Andrea Bertoncini
  11. Stefano Panzeri
  12. Carlo Liberale
  13. Tommaso Fellin
(2020)
Extended field-of-view ultrathin microendoscopes for high-resolution two-photon imaging with minimal invasiveness
eLife 9:e58882.
https://doi.org/10.7554/eLife.58882

Further reading

    1. Medicine
    2. Neuroscience
    Sachin Sharma, Russell Littman ... Olujimi A Ajijola
    Research Article Updated

    The cell bodies of postganglionic sympathetic neurons innervating the heart primarily reside in the stellate ganglion (SG), alongside neurons innervating other organs and tissues. Whether cardiac-innervating stellate ganglionic neurons (SGNs) exhibit diversity and distinction from those innervating other tissues is not known. To identify and resolve the transcriptomic profiles of SGNs innervating the heart, we leveraged retrograde tracing techniques using adeno-associated virus (AAV) expressing fluorescent proteins (GFP or Td-tomato) with single cell RNA sequencing. We investigated electrophysiologic, morphologic, and physiologic roles for subsets of cardiac-specific neurons and found that three of five adrenergic SGN subtypes innervate the heart. These three subtypes stratify into two subpopulations; high (NA1a) and low (NA1b and NA1c) neuropeptide-Y (NPY) -expressing cells, exhibit distinct morphological, neurochemical, and electrophysiologic characteristics. In physiologic studies in transgenic mouse models modulating NPY signaling, we identified differential control of cardiac responses by these two subpopulations to high and low stress states. These findings provide novel insights into the unique properties of neurons responsible for cardiac sympathetic regulation, with implications for novel strategies to target specific neuronal subtypes for sympathetic blockade in cardiac disease.

    1. Neuroscience
    Hang Hu, Rachel E Hostetler, Ariel Agmon
    Research Article Updated

    Oscillations of extracellular voltage, reflecting synchronous, rhythmic activity in large populations of neurons, are a ubiquitous feature in the mammalian brain, and are thought to subserve important, if not fully understood roles in normal and abnormal brain function. Oscillations at different frequency bands are hallmarks of specific brain and behavioral states. At the higher end of the spectrum, 150-200 Hz ripples occur in the hippocampus during slow-wave sleep, and ultrafast (400-600 Hz) oscillations arise in the somatosensory cortices of humans and several other mammalian species in response to peripheral nerve stimulation or punctate sensory stimuli. Here we report that brief optogenetic activation of thalamocortical axons, in brain slices from mouse somatosensory (barrel) cortex, elicited in the thalamorecipient layer local field potential (LFP) oscillations which we dubbed “ripplets”. Ripplets originated in the postsynaptic cortical network and consisted of a precisely repeating sequence of 2‑5 negative transients, closely resembling hippocampal ripples but, at ~400 Hz, over twice as fast. Fast-spiking (FS) inhibitory interneurons fired highly synchronous 400 Hz spike bursts entrained to the LFP oscillation, while regular-spiking (RS), excitatory neurons typically fired only 1-2 spikes per ripplet, in antiphase to FS spikes, and received synchronous sequences of alternating excitatory and inhibitory inputs. We suggest that ripplets are an intrinsically generated cortical response to a strong, synchronous thalamocortical volley, and could provide increased bandwidth for encoding and transmitting sensory information. Importantly, optogenetically induced ripplets are a uniquely accessible model system for studying synaptic mechanisms of fast and ultrafast cortical and hippocampal oscillations.