Austin Leikvoll
Prakash Kara- Department of Neuroscience, University of Minnesota, Minneapolis, MN, United States
Two-photon imaging of genetically-encoded calcium indicators (GECIs) has traditionally relied on intracranial injections of adeno-associated virus (AAV) or transgenic animals to achieve expression. Intracranial injections require an invasive surgery and result in a relatively small volume of tissue labeling. Transgenic animals, although they can have brain-wide GECI expression, often express GECIs in only a small subset of neurons, may have abnormal behavioral phenotypes, and are currently limited to older generations of GECIs. Inspired by recent developments in the synthesis of AAVs that readily cross the blood brain barrier, we tested whether an alternative strategy of intravenously injecting AAV-PHP.eB is suitable for two-photon calcium imaging of neurons over many months after injection. We injected C57BL/6 J mice with AAV-PHP.eB-Synapsin-jGCaMP7s via the retro-orbital sinus. After allowing 5 to 34 weeks for expression, we performed conventional and widefield two-photon imaging of layers 2/3, 4 and 5 of the primary visual cortex. We found reproducible trial-by-trial neural responses and tuning properties consistent with known feature selectivity in the visual cortex. Thus, intravenous injection of AAV-PHP.eB does not interfere with the normal processing in neural circuits. In vivo and histological images show no nuclear expression of jGCaMP7s for at least 34 weeks post-injection.
1. Introduction
Two-photon calcium imaging is a broadly used technique for recording neuronal population activity with cellular resolution in vivo (Grewe and Helmchen, 2009; Grienberger and Konnerth, 2012; Grienberger et al., 2022). Conventional two-photon microscopes record from relatively small fields of view (<1 × 1 mm). The very first in vivo two-photon calcium imaging study from the neocortex imaged 100–200 neurons in small patches of neocortex (<300 × 300 μm regions; Stosiek et al., 2003). Imaging such small regions of neural tissue has yielded many insights into local microcircuit processing in various sensory systems including visual, olfactory, auditory, taste and somatosensory cortices (Ohki et al., 2005; Li et al., 2008; Kara and Boyd, 2009; Stettler and Axel, 2009; Bandyopadhyay et al., 2010; Chen et al., 2011a; Clancy et al., 2015). Incremental improvements in laser power, objective lens technology, scanning speed, and detector sensitivity have all contributed to an explosive growth in two-photon calcium imaging (Litvina et al., 2019; Abdelfattah et al., 2022). Moreover, the BRAIN Initiative has been a major driving force in creating a new class of optical instruments that can measure computation across large populations of neurons (with cellular resolution) over multiple cortical regions (Jorgenson et al., 2015). Large field of view two-photon mesoscopes can now image regions up to 5 × 5 mm with sub-cellular resolution (Sofroniew et al., 2016; Demas et al., 2021; Yu et al., 2021). Exceeding these advancements in instrumentation, innovations regarding the mode of introducing calcium sensors into neurons and the sensitivity of these sensors have been rather remarkable over the past 20 years.
The first calcium imaging studies in vivo relied on synthetic chemical indicator dyes like Oregon Green BAPTA-1 AM (OGB-1 AM) (Stosiek et al., 2003; Ohki et al., 2005). Then genetically encoded calcium indicators (GECIs), like GCaMP, were progressively refined for high signal-to-noise changes in fluorescence, rapid response kinetics and long-term labeling of neurons (Nakai et al., 2001; Palmer and Tsien, 2006; Tian et al., 2009; Chen et al., 2013). Initially, GECIs were introduced to the brain with intracranial injections of adeno-associated virus (AAV) carrying a GECI transgene. Then, numerous transgenic mouse lines for jGCaMP6 (Thy1, Emx1, Ai94, Ai162) and GCaMP7 (G7NG817, TRE-G) were created (Dana et al., 2014; Sato et al., 2015; Monai et al., 2016; Daigle et al., 2018).
Intracranial AAV injections are not ideal for large field of view calcium imaging. Most notably, intracranial injections using micropipettes label a relatively small volume of tissue, typically a sphere ~1–2 mm in diameter, so labeling large volumes of the cortex with this method requires AAV injections be made at multiple sites. Intracranial injections also create a gradient effect in labeling–neurons closer to the center of the injection site express more GECI and are thus brighter than neurons at the periphery (Tian et al., 2009; Chen et al., 2013; Bedbrook et al., 2018). Because local administration of AAV delivers a high copy number of viral genomes to each transduced cell, this strategy can result in excessive levels of GECI expression, causing GECIs to accumulate in cell nuclei and potentially leading to neuronal dysfunction (Tian et al., 2009; Chen et al., 2013; Resendez et al., 2016; Yang et al., 2018; Dana et al., 2019). Intracranial injections also have practical disadvantages—they are relatively time-consuming and invasive procedures, and they require extensive practice to perform without damaging the underlying neural tissue.
GECI-expressing transgenic animals have the obvious advantage of brain-wide labeling (Dana et al., 2014, 2018). However, there are potential caveats to this approach, as some GCaMP lines have abnormal behavioral phenotypes likely arising from aberrant neocortical activity (Steinmetz et al., 2017). Moreover, two-photon laser powers used for in vivo imaging in transgenic lines are rarely provided, but when cited, they appear to be much higher than the average powers used in experiments where intracranial injections are performed. For example, the Thy1 line needed 145 mW for imaging layer 2/3 neurons (100–250 μm depth) (Dana et al., 2014). The Emx1 and Ai94 required up to 70 mW for imaging layer 2/3 (Huang et al., 2021). Finally, producing and validating new transgenic animal lines is a relatively slow process, limiting the speed at which new versions of GECIs can be incorporated into transgenic animals.
Virally-mediated, brain-wide expression of GECIs is now possible due to the development of AAV capsid variants capable of crossing the blood–brain barrier following intravenous injection (Deverman et al., 2016; Chan et al., 2017; Goertsen et al., 2022). The AAV-PHP.B capsid (Deverman et al., 2016) has already been used to achieve brain-wide GECI expression for functional imaging, but these studies required a high dose of virus (>1 × 1012 viral genomes per mouse) to achieve adequate expression. The newer AAV-PHP.eB capsid (Chan et al., 2017), an enhanced variant of AAV-PHP.B, is particularly promising for achieving brain-wide GECI expression because it appears to transduce the majority of neurons across cortical layers even with a substantially lower viral titer than its predecessor (Brown et al., 2021).
None of the previous studies using AAV-PHP.B or AAV-PHP.eB determined the reliability of neural responses over repeated trials (Allen et al., 2017; Hillier et al., 2017). In fact, one of the studies included “unreliable” cells that responded to only one stimulus repetition (Hillier et al., 2017). While behavioral modulation can influence trial-by-trial reliability of sensory responses (Stringer et al., 2019), the animals in this previous AAV-PHP.B study were anesthetized (Hillier et al., 2017). In an earlier study using synthetic calcium dyes, sensory-evoked trial-by-trial responses were found to be less reliable with Fluo-4 than OGB-1 (Bandyopadhyay et al., 2010).
Here, we show that intravenous injections of AAV-PHP.eB-Synapsin-jGCaMP7s is suitable for two-photon calcium imaging of sensory-evoked neural responses that have large amplitudes and extremely high trial-by-trial reliability. Moreover, while previous work assessed the viability of expression up to 10 weeks post injection, we demonstrate that this method provides stable jGCaMP7s expression for at least 34 weeks, with maintenance of robust calcium responses and no visible nuclear expression.
2. Materials and methods
2.1. Animals and AAV injections
All animal procedures were approved by the Institutional Animal Care and Use Committee at the University of Minnesota. Male and female C57BL/6 J mice of different ages were injected with AAV (n = 11 mice of postnatal age 23–31 days, n = 2 mice of postnatal age 48 days, n = 2 mouse postnatal age 80 days, mass 11–27 g; Jackson # 000664). The mice were anesthetized with isoflurane (4%–5% induction, 1–2% maintenance) and intravenously injected with 60–100 μl of sterile saline containing 4 × 1011 viral genomes of AAV-PHP.eB via the retro-orbital sinus of the left eye (see Supplementary note). Each AAV-PHP.eB capsid carried a plasmid for jGCaMP7s driven by the synapsin promoter and WPRE regulatory element (AAV-PHP.eB-Synapsin-jGCaMP7s-WPRE; Addgene # 104487-PHPeB). Virus injections took less than 10 min from the beginning of anesthesia induction to full recovery. Mice were imaged with two-photon microscopy 5 to 34 weeks after AAV injection. Note that the creators of AAV-PHP.eB found that even though expression is present within days after transduction, it does not reach a steady-state level until at least 3–4 weeks after transduction (Challis et al., 2019).
2.2. Surgical procedures
Cranial windows were installed over the left primary visual cortex prior to two-photon imaging. Mice were initially anesthetized with a bolus injection of fentanyl citrate (0.05 mg kg−1), midazolam (5 mg kg−1), and dexmedetomidine (0.25 mg kg−1). During imaging, continuous intraperitoneal infusion with a lower concentration mixture (fentanyl citrate: 0.002–0.03 mg kg−1 h−1, midazolam: 0.2–3.0 mg kg−1 h−1, and dexmedetomidine: 0.010–0.15 mg kg−1 h−1) was administered using a catheter connected to a syringe pump (O’Herron et al., 2020a). Heart and respiration rates, rectal temperature, and SpO2 levels were monitored throughout surgery and imaging. Lidocaine (0.02 ml, 2%) was injected subcutaneously before making incisions. A titanium headplate was attached to the skull using C&B Metabond dental cement (Parkell # S380). A craniectomy (3–4 mm in diameter) was made over the left primary visual cortex (V1), centered ~2.5 mm lateral to lambda and ~ 1.5 mm anterior to the transverse sinus. After performing a durotomy over V1, the cranial windows were sealed with agarose (1.5% in artificial cerebrospinal fluid) and a glass coverslip (5 mm diameter, 0.15 mm thickness; Warner # D263).
2.3. Imaging
Two-photon imaging of the visual cortex was performed with microscope from Bruker coupled to an InSight X3 laser set at 940 nm (Spectra-Physics), as we have used previously (Liu et al., 2020; Cho et al., 2022). Imaging was performed with either a 10X objective lens (NA 0.5, ThorLabs # TL10X–2P) or a 25X objective lens (NA 1.05, Olympus # XLPLN25XWMP2). On both objectives, the correction collar was set to 0.15 mm to compensate for the aberration from the coverslip. Images acquired with the 10X lens were 1.6 mm × 1.6 mm to 1.9 mm × 1.9 mm in size and collected at 1024 × 1,024 pixels to better resolve individual neurons. When imaging with the 10X lens, the integration time per pixel was set to 0.8 μs, and the frame period was 1.68 s. Images collected with the 25X objective were 500 μm × 500 μm to 535 μm × 535 μm in size and collected at 512 × 512 pixels. When imaging with the 25X lens, the integration time per pixel was set to 3.2–3.6 μs, and the frame periods ranged from 1.10–1.35 s.
2.4. Visual stimulation
Drifting square-wave grating stimuli (100% contrast, 1.5 Hz temporal frequency, and 0.033 cycles/degree spatial frequency) were presented on a 17-inch LCD monitor placed 15 cm from the mouse’s right eye. Blackout cloth was wrapped around the objective lens to prevent light from the visual stimulation monitor from entering the microscope’s photomultiplier tubes (O’Herron et al., 2012). Two types of visual stimuli were used, one to assess receptive field location and another to measure orientation selectivity. The retinotopic mapping stimulus was presented in a 3 × 3 grid (Shen et al., 2012; O’Herron et al., 2020b). The retinotopic stimulus consisted of vertically-oriented gratings drifting in one direction of motion within a circle occupying approximately one-ninth of the stimulus monitor. Drifting gratings were displayed in a sequence, starting with position 1 in the upper left of the screen. For determining orientation selectivity, the visual stimuli consisted of drifting gratings presented with 8 directions of motion in 45° steps. During functional imaging runs, visual stimuli were typically displayed for 5 imaging frames interleaved with 10 frames of blank (equiluminant uniform gray across the whole screen). Thus, grating stimuli were displayed for approximately 6 s preceded by approximately 12 s of blank. Each stimulus condition was repeated 3 times.
2.5. Histology
After imaging, mice were euthanized by intraperitoneal injection of pentobarbital (>180 mg kg−1) and transcardially perfused with phosphate-buffered saline (PBS) followed by 4% paraformaldehyde (PFA) in PBS. Brains were removed, placed in 4% PFA solution overnight at 4°C, transferred to PBS, and stored at 4°C. Coronal brain sections (50 μm thick) were cut using a vibratome (PELCO easiSlicer), mounted onto slides, dried overnight, and coverslipped with mountant (Invitrogen # P36980).
2.6. Data analysis
All images were analyzed with code written in MATLAB (MathWorks). Neuronal cell bodies were masked manually as follows. First the entire time series of TIFF images was averaged and then contrast optimized using adaptive histogram equalization. Donut shaped neurons across the entire field of view were then readily visible and could be manually selected, equivalent to using the ImageJ magic wand tool (Walker, 2014). Fluorescence time courses of each neuron F(t) were calculated by averaging over the pixels within each cell body mask. Only raw data were used to calculate time course and subsequent response statistics–no denoising was performed. The response ΔF/F were computed as (F1 − F0)/F0, where F1 was the average fluorescence across the entire stimulus window and F0 was the average fluorescence during the last approximately 6 s of the blank interval. When using retinotopic visual stimuli, significantly responsive neurons were defined by ANOVA across baseline and the nine retinotopic visual stimuli across three trials (p < 0.01). When using orientation and direction specific visual stimuli, significantly responsive neurons were defined by ANOVA across baseline and eight directions over three trials (p < 0.01). Trial-by-trial variability was measured as coefficient of variation (standard deviation/mean) for visually responsive neurons. For each neuron, the coefficient of variation was calculated for the response to its preferred stimulus. Thus, during retinotopic mapping runs, the single stimulus location (of the nine presented) which gave the largest mean response over three trials was used to calculate trial-by-trial variability. Similarly, for orientation/direction mapping runs, the single stimulus direction (of the eight presented) which produced the largest mean response was used to determine trial-by-trial variability.
3. Results
3.1. Expression pattern across cortical layers following intravenous injection of AAV GCaMP7
Three to 34 weeks after intravenous injection, we examined jGCaMP7s expression patterns with in vivo two-photon microscopy or postmortem histological samples in 15 mice e.g., (see Figures 1A–D; Supplementary Figure 1). There are very few bona fide neurons in layer 1 and thus, as expected, little GCaMP7 expression was seen in layer 1. In contrast, robust expression was present in cortical layer 2/3. In previous work, we showed that intracranial injection of AAV.hSyn.GCaMP6s produces weak expression in cortical layer 4 (Liu et al., 2020) and thus we had resorted to OGB-1 AM dye for layer 4 calcium imaging (O’Herron et al., 2020a). Here we show that with AAV-PHP.eB-Synapsin-jGCaMP7s, there is no “gap” in expression in layer 4 when imaging in vivo (Figure 1B; Supplementary Video 1) and when examining coronal histological sections (Figures 1C,D). In vivo, we exceeded the 450 μm two-photon depth limit for functional imaging of fluorophores when they are approximately homogeneously labeled across cortical layers (Takasaki et al., 2020). Thus, our in vivo imaging included layer 5 pyramidal cell bodies at 550 μm (see Figure 1). We could not image layer 6 neurons in vivo with two-photon microscopy, but layer 6 neurons are clearly labeled with this method as evident in histology (Figures 1C,D). In our histological samples, we also saw jGCaMP7s expression in multiple brain regions below the neocortex e.g., (see Figures 1C,D), but here we have focused on the expression pattern of AAV GCaMP7 specifically in V1, from several weeks post-injection of the AAV through more than 8 months post-injection. We saw no nuclear expression even at 16–34 weeks after injection of AAV-PHP.eB (Figure 1D; Supplementary Video 2).
Figure 1. Experimental approach and representative cortical expression pattern following intravenous injection of AAV GCaMP7. (A) Timeline and overview of experimental methods. Mice were intravenously injected with AAV-PHP.eB-Synapsin-jGCaMP7s via the retro-orbital sinus. Five to thirty-four weeks later, sensory-evoked neural responses in the visual cortex were recorded with two-photon microscopy and the tissue was thereafter preserved for postmortem histology. (B) Side view of a three-dimensional volume of GCaMP7 expressing neurons obtained in vivo, 16 weeks after intravenous injection of the AAV-PHP.eB. The x-y-z dimensions of the volume are 750 μm × 750 μm × 700 μm and the first z plane is 50 μm above the pia. The grid line interval is also 50 μm. Yellow arrows indicate the planes of imaging where visually-evoked responses were recorded. (C) Representative histological coronal section of the mouse visual cortex 8 weeks after intravenous injection of AAV-PHP.eB shows neural expression across cortical and sub-cortical layers. The mouse age at AAV injection was 31 days postnatal. (D) Histological coronal section taken 34 weeks after intravenous injection of AAV-PHP.eB. The mouse age at AAV injection was 23 days postnatal.
We next examined the fidelity of trial-by-trial responses in jGCaMP7s-expressing neurons in vivo using two different types of visual stimuli (retinotopy and orientation) in a subset of animals.
3.2. Reliable trial-by-trial neural responses to retinotopic visual stimuli after intravenous injection of AAV GCaMP7
Varying the location of the visual stimulus activated spatially distinct clusters of neurons in mouse V1 (n = 7 mice). Specific neurons responded optimally to a single location over repeated stimulus presentations (see Figures 2A–D). The shift in retinotopic location across the field of view (Figure 2B) matched the known organization of the coarse retinotopic map present in mouse V1 (Schuett et al., 2002; Zhuang et al., 2017). Neurons located centrally in V1 (and thus in the middle of the field of view of imaging) responded to medium elevation visual stimuli, e.g., in Figure 2, neurons 4 and 6 respond to stimulus position 4, which is directly between stimulus positions 1 and 7 (see Figure 2B). Most importantly, the trial-by-trial variability of GCaMP7 neural responses to retinotopic visual stimuli was low. In the large region imaged shown in Figure 2, the median coefficient of variation was 0.46 across all responsive neurons (see Figure 2D). In mouse V1 studies, typically full-field visual stimuli that vary in direction and orientation are presented. Therefore, we next examined the trial-by-trial reliability of neural responses to stimuli varying in direction and orientation.
Figure 2. Two-photon imaging of neural responses to retinotopic (position dependent) visual stimuli after intravenous injection of AAV GCaMP7. (A) In vivo image from layer 2/3 (imaging plane 206 μm below pia) of a large region (1,740 × 1,740 μm) from mouse visual cortex. The average laser power used was 50 mW with a ThorLabs 10× objective lens. The eight circled neurons (numbered 1–8) were randomly selected to show their visually-evoked responses (see panel B). The compass in the bottom right corner (shown in magenta color) represents the medial, lateral, anterior and posterior coordinates of the imaged region. (B) Top panel, schematic of the visual stimulus sequence used for retinotopic mapping. Small patches of drifting gratings were sequentially presented in nine locations within a 3 × 3 grid, starting with position 1 and ending with position 9. Bottom panel, time courses for the eight example neurons (numbered 1–8) shown in panel A. Gray shaded areas in the time course plots correspond to periods of visual stimulation. The bolded numbers above each shaded area represent the stimulus position. Fluorescence traces from three trials are superimposed for each neuron. X and Y axes scales are identical for all eight examples shown. (C) Heat map of calcium responses from all 699 neurons that were masked in the imaged region shown in panel A. Each horizontal line represents a single neuron trace (response averaged across three trials). (D) Box plot showing the visually-evoked trial-by-trial variability (measured as coefficient of variation) for the 551 neurons that were significantly responsive to visual stimuli (p < 0.01 ANOVA across baseline and the 9 retinotopic visual stimuli across three trials). For the box plots in this and subsequent figures, the middle red line represents the median, the box represents the 25 and 75% quartiles, and the whiskers represent the 1 and 99 percentiles.
3.3. Reliable trial-by-trial neural responses to orientation selective visual stimuli after intravenous injection of AAV GCaMP7
In all field of views examined, varying the direction and orientation of visual stimuli produced one population of highly tuned orientation-selective neurons and another population of untuned GCaMP7 neurons (n = 10 mice). Tuned and untuned neurons were always interspersed in the same imaging field of view. For example, compare neurons 1 and 2 in Figures 3A–C; neurons 3 and 4 in Figures 3E–G; neurons 5 and 6 in Figures 3I–K; and neurons 1 and 2 in Supplementary Figure 2. Previous studies using other types of calcium sensors, e.g., OGB-1 AM, established two populations of neurons, tuned excitatory neurons and untuned inhibitory neurons (Sohya et al., 2007; Kerlin et al., 2010). Since we wanted to examine the fidelity of intravenous injection of AAV GCaMP7, without potential interference from other reagents, we did not attempt to tag the various classes of inhibitory neurons. AAV-PHP.eB has already been shown to transduce excitatory neurons and the three major classes of inhibitory neurons (Brown et al., 2021). Most critically, orientation and direction specific stimuli evoked GCaMP7 neural responses with very low trial-by-trial variability. In the imaged regions shown in Figure 3, the median coefficient of variation was 0.59 in layer 2/3 (Figure 3D), 0.55 in layer 4 (Figure 3H) and 0.48 in layer 5 neurons (Figure 3L). Intracranial injections of AAV typically result in viable expression and functional imaging for 1–2 months post injection. Thereafter nuclear expression and attenuation of sensory-evoked responses have been observed (Tian et al., 2009; Chen et al., 2013). Therefore, we next determined if intravenous AAV injections can maintain normal expression and high trial-by-trial reliability of responses for several months post injection (see Figure 4; Supplementary Figure 3).
Figure 3. Two-photon imaging of neural responses to orientation and direction selective visual stimuli after intravenous injection of AAV GCaMP7. (A–D) Data from layer 2/3 (imaging plane 273 μm below the pia). The average laser power used was 30 mW. For the data shown in this and subsequent figures, a high NA 25× Olympus objective lens was used. The two circled neurons in panel A (numbered 1–2) were selected to show tuned and untuned visual responses, respectively (see time courses in panels B,C). For all time courses shown, three trials are overlaid. Panel D shows a box plot of the coefficient of variation of calcium responses for neurons that were significantly responsive to visual stimuli (p < 0.01, ANOVA across baseline and eight directions over three trials). (E–H) Data from layer 4 (imaging plane 400 μm below the pia). The average laser power used was 48 mW. Conventions for time courses and box plots are the same as shown in previous panels. (I–L) Data from layer 5 (imaging plane 550 μm below the pia). Average laser power during imaging was 56 mW. In some of the time courses, e.g., see panels C, G, and K, in the first blank period around time zero, there appears to be a “response.” This is simply residual transient from the 8th visual stimulus because of the robust and sustained GCaMP7s response. If a neuron does not respond to the 8th visual stimulus (315° direction drifting grating in this figure and in Figure 4) then there will be no residual transient in the first blank period. Therefore, the first trial (magenta) has a flat baseline in the first blank period but the second and third trials (blue and green) have the residual response from the 8th visual stimulus from the previous trials.
Figure 4. Robust visual responses persist even after 16 weeks of intravenous injection of AAV GCaMP7. (A–C) Data from layer 2/3. The average laser power during imaging was 26 mW. (D–F) Data from layer 4. The average power used was 32 mW. Conventions for two-photon imaging planes, time courses and box plots, as in previous figures.
3.4. Robust visual responses persist for many months after intravenous injection of GCaMP7
The distribution of animals imaged in vivo across post-injection periods was as follows. Five weeks after intravenous injection n = 1 mouse, 6 weeks post-injection n = 3 mice, 7 weeks post-injection n = 1 mouse, 8 weeks post-injection n = 2 mice, 9 weeks post-injection n = 1 mouse, 10 weeks post-injection n = 1 mouse, 11 weeks post-injection n = 1 mouse, 16 weeks post-injection n = 1 mouse and 34 weeks post-injection n = 1 mouse. Response amplitudes were typically >300% across all periods and all cortical layers examined (see Figures 2–4; Supplementary Figures 2,3 and Supplementary Video 2). The trial-by-trial variability of responses to visual stimuli (measured as the coefficient of variation) was consistently below 0.6 across all imaging periods (5–34 weeks post injection) and in all cortical layers examined (layers 2/3, layer 4 and layer 5; see Figures 2–4; Supplementary Figures 2,3).
4. Discussion
4.1. Reliable trial-by-trial feature selective responses with intravenous injection of AAV GCaMP7
We have shown that intravenous injection of AAV-PHP.eB-Synapsin-jGCaMP7s provides robust and stable expression for high fidelity two-photon imaging of sensory-evoked neural responses in mouse visual cortex. Previously, a comprehensive anatomical study that used advanced quantitative immunocytochemistry showed that intravenous AAV-PHP.eB labels approximately 50% of neocortical neurons (Challis et al., 2019). Here we show that in such a population of neurons across cortical layers, the widespread expression is sufficiently strong for conventional small field of view two-photon imaging and wide-field two-photon imaging. Moreover, feature selectivity such as retinotopy and orientation-selectivity was preserved with intravenous injection of AAV GCaMP7.
Most importantly, trial-by-trial sensory-evoked responses were highly reproducible in these intravenous AAV GCaMP7 experiments in mice, reminiscent of prior work using single-unit recordings from the cat visual cortex (Kara et al., 2000). While animals engaged in behavioral tasks may evoke eye movements and attention modulation that correlate with trial-by-trial variability (Stringer et al., 2019), our animals were anesthetized. Prior work in macaque monkeys has shown that when subjects maintain fixation, trial-by-trial variability of visually-evoked responses is markedly lower during the awake state compared to during anesthesia (Festa et al., 2021). Thus, it is not a foregone conclusion that anesthesia necessitates reliable sensory-evoked responses over repeated trials. For example, in the auditory and visual cortex, there have been reports of highly variable trial-by-trial responses to sensory stimulation during two-photon calcium imaging in anesthetized mice (Bandyopadhyay et al., 2010; Jia et al., 2010; Chen et al., 2011b). The type of anesthesia, the experimental preparation quality and sensor performance can all heavily impact the fidelity of responses (Bandyopadhyay et al., 2010; Grienberger et al., 2022). Even in visual cortical brain slices, monosynaptic activation of spiny stellate neurons in layer 4 can evoke responses with a wide range of trial-by-trial variability (Stratford et al., 1996). For example, the coefficient of variation of responses in layer 4 neurons is nearly an order of magnitude lower when stimulating LGN afferents compared to when stimulating layer 5/6 neurons (Stratford et al., 1996).
4.2. Prior use of intravenous injection of GCaMP6, GCaMP7 and GCaMP8
Previous studies have used AAV-PHP.eB-Syn-GCaMP6s for widefield one-photon imaging in neocortex (Michelson et al., 2019; Xiao et al., 2021). Additional work used AAV-PHP.eB-CAG-FLEX-jGCaMP7s in Cre-driver mice to achieve GECI expression in specific cortical layers and then used 2-photon imaging to measure the effects of anesthesia on spontaneous neuronal activity (Bharioke et al., 2022). However, to our knowledge no previous studies have demonstrated retinotopic and orientation-selective responses at the single-neuron level or the high-fidelity trial-to-trial responses achievable with intravenous injections of AAV-PHP.eB-Synapsin-GCaMP7s. Note that a recent AAV-PHP.eB-GCaMP8 study examined visual cortical responses to orientation and direction selective stimuli by intracranial injection of the AAV, not intravenous injection (Grodem et al., 2023).
4.3. Closing the layer 4 gap with intravenous injection of GCaMP7
Unlike previous work showing that intracranial injection of AAV2/9.hSyn.GCaMP6s.WPRE.SV40 produces very weak expression in cortical layer 4 (Liu et al., 2020), here we show that with intravenous injection of AAV-PHP.eB-Synapsin-jGCaMP7s, there is no “gap” in layer 4 expression.
4.4. Viable expression and high-fidelity sensory responses persist for many months after intravenous injection of GCaMP7
In all past studies examining AAV-PHP.B- and AAV-PHP.eB-mediated GCaMP expression, the post-injection period examined was 2–10 weeks (Allen et al., 2017; Hillier et al., 2017; Michelson et al., 2019; Grodem et al., 2023). We have extended this timeframe by demonstrating that expression remains stable and free of nuclear labeling for at least 34 weeks after transduction. Moreover, at this age, trial-by-trial responses had low variability, just as we found at earlier time points.
4.5. Future directions
Intravenous injection of AAV GCaMP7 appears to be suitable for long-term chronic two-photon imaging across multiple brain regions in projects that involve learning and memory, recovery from experience-dependent plasticity, aging and Alzheimer’s disease. The 34 week post-injection example we show (Supplementary Figure 3) is from a 37-week-old mouse, which is considered “middle-age” (Flurkey et al., 2007). It is possible that quality expression and reliable responses will be present even in “old age” mice (>72 weeks).
For any of the long-term imaging paradigms mentioned above, a key advantage of AAV-mediated GECI expression over transgenic animals is flexibility. To change the transgene carried by an AAV capsid is straightforward, whereas creating new transgenic animal lines is a difficult and lengthy process. This flexibility is especially advantageous for calcium imaging because new generations of GECIs with enhanced performance and cell-type specific promoters are regularly engineered (Grodem et al., 2023). Additionally, intravenous AAV injections could be used to label neurons in other transgenic lines, e.g., disease models (as noted above) without the need for crossbreeding with a transgenic GECI line.
Finally, it is important to note that the AAV-PHP.eB capsid can efficiently cross the blood brain barrier only in a subset of mouse strains, such as C57BL/6 J, which have particular isoforms of the Ly6a receptor (Hordeaux et al., 2019; Huang et al., 2019; Batista et al., 2020). However, several new AAV variants cross the blood brain barrier in all mouse strains and even in non-rodent species, including non-human primates (Ravindra Kumar et al., 2020; Goertsen et al., 2022). It is expected that refinement of these and other newly developed tools, such as soma-and axon-targeted expression of GCaMP (Grodem et al., 2023), will make two-photon calcium imaging in all mammalian species completely ‘plug-and-play’ in the near future.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The animal study was reviewed and approved by IACUC University of Minnesota.
Author contributions
AL and PK designed experiments. AL collected data. AL and PK performed data analysis and wrote the manuscript. Both authors contributed to the article and approved the submitted version.
Funding
This work was supported by NIH R01 MH111447, U01 NS115585, and NSF 1707287.
Acknowledgments
We thank D. Farinella and H. Jayakumar for optics support in the Kara laboratory, Viviana Gradinaru at Caltech for developing the PHP.eB technology and donating it to Addgene. We also thank Aaron Kerlin, Madhu Kannan, and Zane Crabtree for comments on the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2023.1181828/full#supplementary-material
References
Abdelfattah, A. S., Ahuja, S., Akkin, T., Allu, S. R., Brake, J., Boas, D. A., et al. (2022). Neurophotonic tools for microscopic measurements and manipulation: status report. Neurophotonics 9:013001. doi: 10.1117/1.NPh.9.S1.013001
Allen, W. E., Kauvar, I. V., Chen, M. Z., Richman, E. B., Yang, S. J., Chan, K., et al. (2017). Global representations of goal-directed behavior in distinct cell types of mouse neocortex. Neuron 94, 891–907.e6. doi: 10.1016/j.neuron.2017.04.017
Bandyopadhyay, S., Shamma, S. A., and Kanold, P. O. (2010). Dichotomy of functional organization in the mouse auditory cortex. Nat. Neurosci. 13, 361–368. doi: 10.1038/nn.2490
Batista, A. R., King, O. D., Reardon, C. P., Davis, C., Shankaracharya,, Philip, V., et al. (2020). Ly6a differential expression in blood–brain barrier is responsible for strain specific central nervous system transduction profile of AAV-PHP.B. Hum. Gene Ther. 31, 90–102. doi: 10.1089/hum.2019.186
Bedbrook, C. N., Deverman, B. E., and Gradinaru, V. (2018). Viral strategies for targeting the central and peripheral nervous systems. Annu. Rev. Neurosci. 41, 323–348. doi: 10.1146/annurev-neuro-080317-062048
Bharioke, A., Munz, M., Brignall, A., Kosche, G., Eizinger, M. F., Ledergerber, N., et al. (2022). General anesthesia globally synchronizes activity selectively in layer 5 cortical pyramidal neurons. Neuron 110, 2024–2040.e10. doi: 10.1016/j.neuron.2022.03.032
Brown, D., Altermatt, M., Dobreva, T., Chen, S., Wang, A., Thomson, M., et al. (2021). Deep parallel characterization of AAV tropism and AAV-mediated transcriptional changes via single-cell RNA sequencing. Front. Immunol. 12:730825. doi: 10.3389/fimmu.2021.730825
Challis, R. C., Ravindra Kumar, S., Chan, K. Y., Challis, C., Beadle, K., Jang, M. J., et al. (2019). Systemic AAV vectors for widespread and targeted gene delivery in rodents. Nat. Protoc. 14, 379–414. doi: 10.1038/s41596-018-0097-3
Chan, K. Y., Jang, M. J., Yoo, B. B., Greenbaum, A., Ravi, N., Wu, W.-L., et al. (2017). Engineered AAVs for efficient noninvasive gene delivery to the central and peripheral nervous systems. Nat. Neurosci. 20, 1172–1179. doi: 10.1038/nn.4593
Chen, X., Gabitto, M., Peng, Y., Ryba, N. J., and Zuker, C. S. (2011a). A gustotopic map of taste qualities in the mammalian brain. Science 333, 1262–1266. doi: 10.1126/science.1204076
Chen, X., Leischner, U., Rochefort, N. L., Nelken, I., and Konnerth, A. (2011b). Functional mapping of single spines in cortical neurons in vivo. Nature 475, 501–505. doi: 10.1038/nature10193
Chen, T. W., Wardill, T. J., Sun, Y., Pulver, S. R., Renninger, S. L., Baohan, A., et al. (2013). Ultrasensitive fluorescent proteins for imaging neuronal activity. Nature 499, 295–300. doi: 10.1038/nature12354
Cho, S., Roy, A., Liu, C. J., Idiyatullin, D., Zhu, W., Zhang, Y., et al. (2022). Cortical layer-specific differences in stimulus selectivity revealed with high-field fMRI and single-vessel resolution optical imaging of the primary visual cortex. NeuroImage 251:118978. doi: 10.1016/j.neuroimage.2022.118978
Clancy, K. B., Schnepel, P., Rao, A. T., and Feldman, D. E. (2015). Structure of a single whisker representation in layer 2 of mouse somatosensory cortex. J. Neurosci. 35, 3946–3958. doi: 10.1523/JNEUROSCI.3887-14.2015
Daigle, T. L., Madisen, L., Hage, T. A., Valley, M. T., Knoblich, U., Larsen, R. S., et al. (2018). A suite of transgenic driver and reporter mouse lines with enhanced brain-cell-type targeting and functionality. Cells 174:e422, 465–480.e22. doi: 10.1016/j.cell.2018.06.035
Dana, H., Chen, T.-W., Hu, A., Shields, B. C., Guo, C., Looger, L. L., et al. (2014). Thy1-GCaMP6 transgenic mice for neuronal population imaging in vivo. PLoS One 9:e108697. doi: 10.1371/journal.pone.0108697
Dana, H., Novak, O., Guardado-Montesino, M., Fransen, J. W., Hu, A., Borghuis, B. G., et al. (2018). Thy1 transgenic mice expressing the red fluorescent calcium indicator jRGECO1a for neuronal population imaging in vivo. PLoS One 13:e0205444. doi: 10.1371/journal.pone.0205444
Dana, H., Sun, Y., Mohar, B., Hulse, B. K., Kerlin, A. M., Hasseman, J. P., et al. (2019). High-performance calcium sensors for imaging activity in neuronal populations and microcompartments. Nat. Methods 16, 649–657. doi: 10.1038/s41592-019-0435-6
Demas, J., Manley, J., Tejera, F., Barber, K., Kim, H., Traub, F. M., et al. (2021). High-speed, cortex-wide volumetric recording of neuroactivity at cellular resolution using light beads microscopy. Nat. Methods 18, 1103–1111. doi: 10.1038/s41592-021-01239-8
Deverman, B. E., Pravdo, P. L., Simpson, B. P., Kumar, S. R., Chan, K. Y., Banerjee, A., et al. (2016). Cre-dependent selection yields AAV variants for widespread gene transfer to the adult brain. Nat. Biotechnol. 34, 204–209. doi: 10.1038/nbt.3440
Festa, D., Aschner, A., Davila, A., Kohn, A., and Coen-Cagli, R. (2021). Neuronal variability reflects probabilistic inference tuned to natural image statistics. Nat. Commun. 12:3635. doi: 10.1038/s41467-021-23838-x
Flurkey, K., Currer, J. M., and Harrison, D. E. (2007). The mouse in aging research in The mouse in biomedical research, (2nd edn.), eds. J. G. Fox, M. T. Davisson, F. W. Quimby, S. W. Barthold, C. E. Newcomer, and A. L. Smith. Second ed Burlington, MA: American College of Laboratory Animal Medicine-Elsevier, 637–672.
Goertsen, D., Flytzanis, N. C., Goeden, N., Chuapoco, M. R., Cummins, A., Chen, Y., et al. (2022). AAV capsid variants with brain-wide transgene expression and decreased liver targeting after intravenous delivery in mouse and marmoset. Nat. Neurosci. 25, 106–115. doi: 10.1038/s41593-021-00969-4
Grewe, B. F., and Helmchen, F. (2009). Optical probing of neuronal ensemble activity. Curr. Opin. Neurobiol. 19, 520–529. doi: 10.1016/j.conb.2009.09.003
Grienberger, C., Giovannucci, A., Zeiger, W., and Portera-Cailliau, C. (2022). Two-photon calcium imaging of neuronal activity. Nat. Rev. Methods Primers 2:67. doi: 10.1038/s43586-022-00147-1
Grienberger, C., and Konnerth, A. (2012). Imaging calcium in neurons. Neuron 73, 862–885. doi: 10.1016/j.neuron.2012.02.011
Grodem, S., Nymoen, I., Vatne, G. H., Rogge, F. S., Bjornsdottir, V., Lensjo, K. K., et al. (2023). An updated suite of viral vectors for in vivo calcium imaging using intracerebral and retro-orbital injections in male mice. Nat. Commun. 14:608. doi: 10.1038/s41467-023-36324-3
Hillier, D., Fiscella, M., Drinnenberg, A., Trenholm, S., Rompani, S. B., Raics, Z., et al. (2017). Causal evidence for retina-dependent and-independent visual motion computations in mouse cortex. Nat. Neurosci. 20, 960–968. doi: 10.1038/nn.4566
Hordeaux, J., Yuan, Y., Clark, P. M., Wang, Q., Martino, R. A., Sims, J. J., et al. (2019). The GPI-linked protein LY6A drives AAV-PHP.B transport across the blood-brain barrier. Mol. Ther. 27, 912–921. doi: 10.1016/j.ymthe.2019.02.013
Huang, Q., Chan, K. Y., Tobey, I. G., Chan, Y. A., Poterba, T., Boutros, C. L., et al. (2019). Delivering genes across the blood-brain barrier: LY6A, a novel cellular receptor for AAV-PHP.B capsids. PLoS One 14:e0225206. doi: 10.1371/journal.pone.0225206
Huang, L., Ledochowitsch, P., Knoblich, U., Lecoq, J., Murphy, G. J., Reid, R. C., et al. (2021). Relationship between simultaneously recorded spiking activity and fluorescence signal in GCaMP6 transgenic mice. elife 10:e51675. doi: 10.7554/eLife.51675
Jia, H., Rochefort, N. L., Chen, X., and Konnerth, A. (2010). Dendritic organization of sensory input to cortical neurons in vivo. Nature 464, 1307–1312. doi: 10.1038/nature08947
Jorgenson, L. A., Newsome, W. T., Anderson, D. J., Bargmann, C. I., Brown, E. N., Deisseroth, K., et al. (2015). The BRAIN initiative: developing technology to catalyse neuroscience discovery. Philos. Trans. R. Soc. Lond. Ser. B Biol. Sci. 370:20140164. doi: 10.1098/rstb.2014.0164
Kara, P., and Boyd, J. D. (2009). A micro-architecture for binocular disparity and ocular dominance in visual cortex. Nature 458, 627–631. doi: 10.1038/nature07721
Kara, P., Reinagel, P., and Reid, R. C. (2000). Low response variability in simultaneously recorded retinal, thalamic, and cortical neurons. Neuron 27, 635–646. doi: 10.1016/S0896-6273(00)00072-6
Kerlin, A. M., Andermann, M. L., Berezovskii, V. K., and Reid, R. C. (2010). Broadly tuned response properties of diverse inhibitory neuron subtypes in mouse visual cortex. Neuron 67, 858–871. doi: 10.1016/j.neuron.2010.08.002
Li, Y., Van Hooser, S. D., Mazurek, M., White, L. E., and Fitzpatrick, D. (2008). Experience with moving visual stimuli drives the early development of cortical direction selectivity. Nature 456, 952–956. doi: 10.1038/nature07417
Litvina, E., Adams, A., Barth, A., Bruchez, M., Carson, J., Chung, J. E., et al. (2019). BRAIN initiative: cutting-edge tools and resources for the community. J. Neurosci. 39, 8275–8284. doi: 10.1523/JNEUROSCI.1169-19.2019
Liu, C. J., Roy, A., Simons, A. A., Farinella, D. M., and Kara, P. (2020). Three-photon imaging of synthetic dyes in deep layers of the neocortex. Sci. Rep. 10:16351. doi: 10.1038/s41598-020-73438-w
Michelson, N. J., Vanni, M. P., and Murphy, T. H. (2019). Comparison between transgenic and AAV-PHP.eB-mediated expression of GCaMP6s using in vivo wide-field functional imaging of brain activity. Neurophotonics 6:1. doi: 10.1117/1.NPh.6.2.025014
Monai, H., Ohkura, M., Tanaka, M., Oe, Y., Konno, A., Hirai, H., et al. (2016). Calcium imaging reveals glial involvement in transcranial direct current stimulation-induced plasticity in mouse brain. Nat. Commun. 7:11100. doi: 10.1038/ncomms11100
Nakai, J., Ohkura, M., and Imoto, K. (2001). A high signal-to-noise Ca2+ probe composed of a single green fluorescent protein. Nat. Biotechnol. 19, 137–141. doi: 10.1038/84397
O’Herron, P., Levy, M., Woodward, J. J., and Kara, P. (2020a). An unexpected dependence of cortical depth in shaping neural responsiveness and selectivity in mouse visual cortex. eNeuro 7:7. doi: 10.1523/ENEURO.0497-19.2020
O’Herron, P., Shen, Z., Lu, Z., Schramm, A. E., Levy, M., and Kara, P. (2012). Targeted labeling of neurons in a specific functional micro-domain of the neocortex by combining intrinsic signal and two-photon imaging. J. Vis. Exp. 70:e50025. doi: 10.3791/50025
O’Herron, P., Summers, P. M., Shih, A. Y., Kara, P., and Woodward, J. J. (2020b). In vivo two-photon imaging of neuronal and brain vascular responses in mice chronically exposed to ethanol. Alcohol 85, 41–47. doi: 10.1016/j.alcohol.2019.12.001
Ohki, K., Chung, S., Ch’ng, Y. H., Kara, P., and Reid, R. C. (2005). Functional imaging with cellular resolution reveals precise micro-architecture in visual cortex. Nature 433, 597–603. doi: 10.1038/nature03274
Palmer, A. E., and Tsien, R. Y. (2006). Measuring calcium signaling using genetically targetable fluorescent indicators. Nat. Protoc. 1, 1057–1065. doi: 10.1038/nprot.2006.172
Ravindra Kumar, S., Miles, T. F., Chen, X., Brown, D., Dobreva, T., Huang, Q., et al. (2020). Multiplexed Cre-dependent selection yields systemic AAVs for targeting distinct brain cell types. Nat. Methods 17, 541–550. doi: 10.1038/s41592-020-0799-7
Resendez, S. L., Jennings, J. H., Ung, R. L., Namboodiri, V. M. K., Zhou, Z. C., Otis, J. M., et al. (2016). Visualization of cortical, subcortical and deep brain neural circuit dynamics during naturalistic mammalian behavior with head-mounted microscopes and chronically implanted lenses. Nat. Protoc. 11, 566–597. doi: 10.1038/nprot.2016.021
Sato, M., Kawano, M., Ohkura, M., Gengyo-Ando, K., Nakai, J., and Hayashi, Y. (2015). Generation and imaging of transgenic mice that express G-CaMP7 under a tetracycline response element. PLoS One 10:e0125354. doi: 10.1371/journal.pone.0125354
Schuett, S., Bonhoeffer, T., and Hubener, M. (2002). Mapping retinotopic structure in mouse visual cortex with optical imaging. J. Neurosci. 22, 6549–6559. doi: 10.1523/JNEUROSCI.22-15-06549.2002
Shen, Z., Lu, Z., Chhatbar, P. Y., O’Herron, P., and Kara, P. (2012). An artery-specific fluorescent dye for studying neurovascular coupling. Nat. Methods 9, 273–276. doi: 10.1038/nmeth.1857
Sofroniew, N. J., Flickinger, D., King, J., and Svoboda, K. (2016). A large field of view two-photon mesoscope with subcellular resolution for in vivo imaging. elife 5:e14472. doi: 10.7554/eLife.14472
Sohya, K., Kameyama, K., Yanagawa, Y., Obata, K., and Tsumoto, T. (2007). GABAergic neurons are less selective to stimulus orientation than excitatory neurons in layer II/III of visual cortex, as revealed by in vivo functional Ca2+ imaging in transgenic mice. J. Neurosci. 27, 2145–2149. doi: 10.1523/jneurosci.4641-06.2007
Steinmetz, N. A., Buetfering, C., Lecoq, J., Lee, C. R., Peters, A. J., Jacobs, E. A. K., et al. (2017). Aberrant cortical activity in multiple GCaMP6-expressing transgenic mouse lines. eneuro 4:ENEURO.0207-0217.2017. doi: 10.1523/ENEURO.0207-17.2017
Stettler, D. D., and Axel, R. (2009). Representations of odor in the piriform cortex. Neuron 63, 854–864. doi: 10.1016/j.neuron.2009.09.005
Stosiek, C., Garaschuk, O., Holthoff, K., and Konnerth, A. (2003). In vivo two-photon calcium imaging of neuronal networks. Proc. Natl. Acad. Sci. U. S. A. 100, 7319–7324. doi: 10.1073/pnas.1232232100
Stratford, K. J., Tarczy-Hornoch, K., Martin, K. A., Bannister, N. J., and Jack, J. J. (1996). Excitatory synaptic inputs to spiny stellate cells in cat visual cortex. Nature 382, 258–261. doi: 10.1038/382258a0
Stringer, C., Pachitariu, M., Steinmetz, N., Reddy, C. B., Carandini, M., and Harris, K. D. (2019). Spontaneous behaviors drive multidimensional, brainwide activity. Science 364:255. doi: 10.1126/science.aav7893
Takasaki, K., Abbasi-Asl, R., and Waters, J. (2020). Superficial bound of the depth limit of two-photon imaging in mouse brain. eNeuro 7:ENEURO.0255-0219.2019. doi: 10.1523/ENEURO.0255-19.2019
Tian, L., Hires, S. A., Mao, T., Huber, D., Chiappe, M. E., Chalasani, S. H., et al. (2009). Imaging neural activity in worms, flies and mice with improved GCaMP calcium indicators. Nat. Methods 6, 875–881. doi: 10.1038/nmeth.1398
Walker, T. (2014). Cell magic wand tool. Available at: https://github.com/fitzlab/CellMagicWand.
Xiao, D., Forys, B. J., Vanni, M. P., and Murphy, T. H. (2021). MesoNet allows automated scaling and segmentation of mouse mesoscale cortical maps using machine learning. Nat. Commun. 12:5992. doi: 10.1038/s41467-021-26255-2
Yang, Y., Liu, N., He, Y., Liu, Y., Ge, L., Zou, L., et al. (2018). Improved calcium sensor GCaMP-X overcomes the calcium channel perturbations induced by the calmodulin in GCaMP. Nat. Commun. 9:1504. doi: 10.1038/s41467-018-03719-6
Yu, C. H., Stirman, J. N., Yu, Y., Hira, R., and Smith, S. L. (2021). Diesel2p mesoscope with dual independent scan engines for flexible capture of dynamics in distributed neural circuitry. Nat. Commun. 12:6639. doi: 10.1038/s41467-021-26736-4
Keywords: brain-wide calcium imaging, two-photon microscopy, AAV-PHP GCaMP7, visual cortex (V1), blood brain barrier, trial-by-trial variability, long-term imaging
Citation: Leikvoll A and Kara P (2023) High fidelity sensory-evoked responses in neocortex after intravenous injection of genetically encoded calcium sensors. Front. Neurosci. 17:1181828. doi: 10.3389/fnins.2023.1181828
Edited by:
Fan Lam, University of Illinois at Urbana-Champaign, United StatesReviewed by:
Daniel Llano, University of Illinois at Urbana-Champaign, United StatesChen Ran, Harvard Medical School, United States
Copyright © 2023 Leikvoll and Kara. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Prakash Kara, cGthcmFAdW1uLmVkdQ==