Abstract
C57BL/6 is the most commonly used mouse strain in neurobehavioral research, serving as a background for multiple transgenic lines. However, C57BL/6 exhibit behavioral and sensorimotor disadvantages that worsen with age. We bred FVB/NJ females and C57BL/6J males to generate first-generation hybrid offspring (FVB/NJ x C57BL/6J)F1. The hybrid mice exhibit reduced anxiety-like behavior, improved learning, and enhanced long-term spatial memory. In contrast to both progenitors, hybrids maintain sensorimotor performance upon aging and exhibit improved long-term memory. The hybrids are larger than C57BL/6J, exhibiting enhanced running behavior on a linear track during freely-moving electrophysiological recordings. Hybrids exhibit typical rate and phase coding of space by CA1 pyramidal cells. Hybrids generated by crossing FVB/NJ females with transgenic males of a C57BL/6 background support optogenetic neuronal control in neocortex and hippocampus. The hybrid mice provide an improved model for neurobehavioral studies combining complex behavior, electrophysiology, and genetic tools readily available in C57BL/6 mice.
Significance Statement
Because of genetic tools, mice are increasingly used for neuroscience experiments that were traditionally performed using rats. However, mice are smaller and “less intelligent” than rats, limiting the size of neural implants and the complexity of performed behaviors. Here, we show that these problems are exacerbated by the widespread use of inbred C57BL/6J mice, which are small animals with behavioral and sensorimotor deficits that worsen with age. In line with the established principle of hybrid vigor, we found that first-generation offspring of C57BL/6J and FVB/NJ mice exhibit improved learning and memory, maintain sensorimotor performance on aging, are larger, and exhibit improved behavior while carrying a chronic implant. The hybrid mice allow genetic control while supporting complex behavior over prolonged durations.
Introduction
Many advances in behavioral and biomedical research rely on lab animals. In choosing an animal model, there are always two conflicting considerations. On the one hand, an organism as similar as possible to humans is desired. On the other hand, the least advanced organism for answering the research question is preferred because of ethical considerations. The lab mouse (Mus musculus) has emerged as a key model in balancing these demands, sharing multiple physiological systems and genes with humans while exhibiting simple mating. Over the years, numerous inbred mouse strains have been developed, allowing genetic modifications, labeling, and manipulation of specific proteins, cells, and organs (Beck et al., 2000). However, there are phenotypic differences between inbred strains used as background for transgenic mice (Upchurch and Wehner, 1988; Võikar et al., 2001; Wahlsten et al., 2003; Brown and Wong, 2007; O’Leary et al., 2011, 2013; Kafkafi et al., 2017). Hence, mouse strain selection is essential when designing a scientific project.
The mouse strain used most often in neurobehavioral research is C57BL/6. The C57BL/6 Jackson Laboratory substrain, C57BL/6J, provided the first extensively sequenced mouse genome (Mouse Genome Sequencing Consortium, 2002) and is among the most widely used inbred strains (Altman and Katz, 1979; Mekada et al., 2009). An advantage of inbreeding is that differences between individuals of the same strain are minimal. The genetic similarity of inbred animals is especially useful for knock-out studies, which may require homozygous animals (Silva et al., 1997). Indeed, a large variety of transgenic mice on the C57BL/6 background is available. However, inbreeding exposes undesired phenotypical traits because of homozygous recessive alleles and does not guarantee stability over generations because of genetic drift (Brekke et al., 2018). C57BL/6-derived mice exhibit known phenotypic disadvantages, including sensitivity to pain, addiction, impaired balance, age-dependent hearing loss, and increased anxiety (Crawley, 1996; Mogil et al., 1999; Ouagazzal et al., 2006). Thus, the behavioral tasks studied using C57BL/6 are limited.
One way to balance the requirements of convenient genetic control and complex behavior is to generate hybrids of C57BL/6J and another strain. Hybrids inherit one allele from the C57BL/6J parent, maintaining transgenic properties in a heterozygous manner. Previous work has shown that offspring of C57BL/6J (C57) and FVB/NJ (FVB) mice consume more alcohol than either progenitor, serving as a preferable model for alcohol consumption (Blednov et al., 2005, 2010). FVB is an inbred albino strain which carries a recessive allele causing retinal degeneration (Pittler and Baehr, 1991; Taketo et al., 1991). The C57 and FVB strains have distinct genealogies (Beck et al., 2000; Mekada et al., 2009), increasing genetic heterogeneity of the hybrid offspring. Thus, hybrids may serve as a potentially favorable animal model.
Here, we bred C57 males with FVB females to determine whether the first-generation hybrids (HYB) are a preferable model for neurobehavioral studies. We tested mice of the three strains (C57, FVB, and HYB) using a battery of standard phenotyping assays (Crawley, 2008; Fuchs et al., 2011). Compared with either inbred strain, HYB exhibited similar sensorimotor performance, reduced anxiety-like behavior, faster learning, and improved memory, which were maintained or further improved on aging. HYB mice were physically larger, and during electrophysiological recordings exhibited enhanced running on a linear track compared with C57, while CA1 neurons had similar place coding properties. Furthermore, transgenic HYB supported optogenetic neuronal control in neocortex and CA1. Together, the results suggest that the hybrid mice constitute a preferable model for systems neuroscience studies combining behavioral tests, electrophysiology, and genetic targeting.
Materials and Methods
Experimental model and subject details
A total of 176 freely-moving adult mice were used in this study. A total of 170 male mice were used for phenotyping, of which 61 were C57BL/6J (C57; JAX #000664, The Jackson Laboratory); 52 were FVB/NJ (FVB; JAX #001800); and 57 were hybrid (HYB; Table 1), offspring of an FVB female and a C57 male. Six transgenic mice were used for electrophysiological recordings (Extended Data Fig. 5-1). One was single-transgenic and hybrid, generated by crossing an FVB female with a parvalbumin (PV)-Cre male (#008069). Two were dual-transgenic, generated by crossing CaMKII-Cre females (JAX #005359) with Ai32 males (#012569). One was dual-transgenic, generated by crossing a PV-Cre female with an Ai32 male. Two were dual-transgenic and hybrid, generated by crossing FVB females with second-generation (CaMKII-Cre x Ai32) males. All mice were bred in-house. After separation from the parents, animals were housed in groups of same-litter siblings. Animals were held on a reverse dark/light cycle (dark phase, from 8 A.M. until 8 P.M.). Data recorded from CA1 during linear track behavior were used in a previous report (Sloin et al., 2022). All animal handling procedures were in accordance with Directive 2010/63/EU of the European Parliament, complied with Israeli Animal Welfare Law (1994), and approved by the Tel Aviv University Institutional Animal Care and Use Committee (IACUC #01-16-051, #01-19-017, and #01-21-051).
Mice used for phenotyping assays
Phenotyping study design
Mice of each strain were divided into two age groups, three and nine months old (Table 1). Each age group was further divided into two subgroups: one subgroup was tested on test Battery A, and a second subgroup was used for test Battery B. Animals assigned to Battery A were subjected to assays in the following order: elevated plus maze (1 d); open field (1 d); rotarod (5 d); treadmill (3 d); optomotor drum (1 d); and forced swim test (1 d). Animals assigned to Battery B were tested on the catwalk (1 d) and on the Morris water maze (MWM; 5 d). Mice were used in the phenotyping cages at least a week postbattery and were counterbalanced from Batteries A and B (nine-month-old mice only). The division of assays into batteries was done to limit the number of assays that each mouse would perform. The order of assays in each battery was designed so that each battery began with voluntary, low stress eliciting assays. Between assays, mice were given a minimum of 24 h without human interaction. All tests were initiated at beginning of the dark phase (8 A.M.) and were administered to a single mouse at a time. All equipment was thoroughly cleaned with Virusolve before and between trials. Except for the rotarod and treadmill tests, behavior during all tests was recorded using GigE cameras (ac1300-60gm mono, Basler; frame rate, 25 Hz). Commercial software (Ethovision 15XT; Noldus Information Technology; Noldus et al., 2001) was used to analyze all video files.
Elevated plus maze
The elevated plus maze is used to evaluate anxiety-like behavior based on the natural uneasiness of rodents toward open, elevated fields (Rodgers and Dalvi, 1997; Crawley, 2000; Komada et al., 2008). The apparatus consists of a four-armed platform resembling a “+” shape, positioned 40 cm above the floor. Two arms are confined by walls (“closed”; L × W × H: 35 × 5 × 15 cm), whereas two other arms are not walled (“open”; 35 × 5 cm). Similar arms face one another. In the test, the mouse was initially placed at the center of the maze facing one of the closed arms and then allowed to move freely for 7 min. The time spent exploring the open arms (“open”) and the time spent exploring the closed arms (“closed”) were used to derive an index (open – closed)/(open + closed) that served as a contraindication to anxiety-like behavior (“contra-anxiety index”).
Forced swim test
The forced swim test is used to evaluate depressive-like behavior based on induced “behavioral despair” (Porsolt et al., 1977; Castagné et al., 2011; Can et al., 2012). The apparatus consists of a clear Plexiglas cylinder (24 cm in height, 19 cm in diameter) filled with 16 cm of water at 22°C. In the test, the mouse was placed in the cylinder for 7 min and then moved to a heated cage until the fur dried completely. Typically, the mouse gradually stopped swimming before being removed. Freeze (immobility) duration, defined as the time the mouse remained floating motionless in the water, was measured during the last 5 min as a manifestation of behavioral despair.
Gait analysis
The gait analysis apparatus (CatWalk XT, Noldus Information Technology) enables the assessment of voluntary gait and locomotion in mice (Ängeby-Möller et al., 2000; Crowley et al., 2018). The apparatus consists of a hardware system with a glass walkway (L × W: 130 × 20 cm), a GigE video camera, and a software package for the quantitative assessment of animal footprints. The walkway is illuminated from above by red light and from the side by green light. The green light is internally reflected within the glass, except at touched points. The walkway is connected to a dark goal box and enclosed by an adjustable tunnel at one end. In the test, mice were placed at the walkway entrance and allowed to run freely. A successful run was defined when the animal traversed the track without pausing. Every animal was tested until three successful runs were completed. For every parameter, the average of the three runs was used for analysis.
Morris water maze
The MWM is used to evaluate learning and spatial long-term memory based on the natural tendency of rodents to attempt to escape a body of water (Morris, 1981, 1984; Vorhees and Williams, 2006; Bromley-Brits et al., 2011). The apparatus is a circular pool (150 cm diameter, 10 cm depth) filled with water and maintained at a temperature of 23°C. A transparent platform (5 × 5 cm) was fixed 1 cm below the water surface at a constant location. The pool was situated in a room containing distal visual cues, and proximal cues were placed on the inner walls of the pool. The mouse was introduced to the pool at different starting points and allowed 60 s to find the platform and additional 15 s to stay on the platform. If the mouse failed to find the platform within 60 s, the animal was guided to the platform and allowed to stay there for 15 s. The animal swam three times every day during four consecutive “learning days.” During the learning days, the time before finding the platform (“latency to platform”) and the total distance traveled were measured. Changes between days were used to evaluate learning performance.
On the fifth day of the assay (“probe day”), the platform was removed, and the mouse was allowed to swim for 60 s. The time spent in every quadrant was recorded. The fraction of time spent in the “target” quadrant (i.e., close to the location of the missing platform), the number of visits to the target quadrant, and the latency to platform were used to evaluate spatial long-term memory performance.
Open field
The open field test provides a way to systematically assess motivational behavior by quantifying exploration of a novel environment and general locomotor activity (Christmas and Maxwell, 1970; Prut and Belzung, 2003). In the test, the mouse was initially placed in one of the corners of an open-top Plexiglas box (L × W: 50 × 50 cm, raised 40 cm above the floor), and behavior was recorded for 15 min. Activity, defined as the fraction of pixels that changed in the entire arena between consecutive frames (frame rate, 25 Hz), was measured throughout the run duration to evaluate exploratory behavior.
Optomotor visual test
The optomotor drum is based on an apparatus developed for immobile mice (Mitchiner et al., 1976) and allows assessing visual behavior in freely-moving mice (Abdeljalil et al., 2005). The apparatus consists of an elevated stationary platform (20 cm) surrounded by a drum (internal diameter, 39 cm) with vertically oriented black and white stripes on the inside, each spanning 10°. In the test, the mouse was initially habituated to the platform for 2 min, with the drum stationary. The drum then rotated at 2 rpm counterclockwise for 2 min, stopped for 30 s, and rotated clockwise for 2 min. The number of head turns (15° movements at the drum speed) and the cumulative duration of head turns were measured.
Phenotyping cages
The PhenoTyper (model 3000, Noldus Information Technology) is an instrumented home cage in which rodent behavior is automatically monitored through a video-based and event-based system (De Visser et al., 2006; Maroteaux et al., 2012; Grieco et al., 2021). The cage (L × W × H: 30 × 30 × 35 cm) is made of transparent Plexiglas walls with an opaque Plexiglas floor. The cage is equipped with a watering station, a feeding station, a running wheel, and a shelter in one corner (H × D: 10 × 9 cm, transparent material). The feeding and watering stations are equipped with beam-breaking devices, allowing automatic recording of feeding behavior and water intake (number of feeds and licks). The lid of the cage lid is equipped with an infrared-sensitive video camera (768 × 576 pixels) and several infrared LEDs, allowing continuous recording of animal position. Bedding covered the floor, and food and water were provided ad libitum. The mouse was introduced to the cage at the beginning of the dark phase (8 A.M.). Testing lasted 3 d, during which no human interference took place.
Rotarod test
Mice were subjected to a five-lane accelerating rotarod (Ugo Basile) to evaluate motor learning and balance (Jones and Roberts, 1968; Crawley, 2000, 2003; Shiotsuki et al., 2010). The apparatus consists of a 3.2 cm (diameter) horizontal rod elevated 10 cm from the ground. The paradigm was comprised of five consecutive days. On every day, mice were subjected to five trials, of which the duration of the three longest trials were averaged. A trial began with the rod rotating at 4 rpm and gradually accelerating for 5 min, up to a maximum of 50 rpm. The duration until the animal fell from the rod (“latency to fall”) was measured.
Treadmill
The treadmill apparatus (Panlab, Harvard Apparatus) is used to assess maximal endurance in mice (Marques-Aleixo et al., 2015). The apparatus consists of a five-lane motorized treadmill with an electric shock zone at one end of the track. Each lane is 38 × 7 × 7 cm (L × W × H). A current shock (0.2 mA) was used to encourage running; cumulative shock duration of 2 s defined the maximal ability of mice to run. The protocol included two training days and one testing day. On the first training day, treadmill speed was fixed at 5 cm/s. The second training day consisted of two parts: (1) walking on the treadmill for 300 s at a constant speed (5 cm/s); (2) subjecting the mice to 330 s of locomotion at a variable speed (increasing from 5 to 21 cm/s by 1 cm/s every 20 s). On the third (testing) day, speed was initially set to 5 cm/s and was gradually increased by 1 cm/s every 20 s. The trial ended when the cumulative shock duration was 2 s. Trial duration and distance to failure were measured.
Indices for behavioral assays
An “aging index” was defined for every parameter as the difference between the median value of that parameter for nine-month-old mice and the median for three-month-old mice of the same strain, divided by the sum. We used a bootstrap procedure to estimate index dispersion and determine statistical significance. The null hypothesis of the index being equal to zero is equivalent to no consistent changes between the two age groups. In the procedure, we resampled the data (e.g., 12 values from three-month-old mice and 15 values from nine-month-old mice) with replacement many (10,000) times. For each resampling iteration, we computed the aging index. The reported index is the mean over all iterations, and the SEM is the SD over all iterations. For a one-sided alternative hypothesis (increase/decrease), p-values are the fraction of indices below/above zero.
A “C57 index” was defined in an equivalent manner, where the index was defined as the difference between the median value of the parameter of interest for the tested strain (FVB or HYB) and the median value of the same parameter for C57 mice of the same age group, divided by the sum. Mean, SEM, and p-values were estimated using bootstrapping as for the aging indices.
When parameters are strictly positive, the difference divided by the sum, (a – b)/(a + b), provides a valid estimate. However, some parameters of interest can take negative values, for instance, the “contra-anxiety index” ranges from −1 to 1. Then, an increase from one negative value to another will yield a negative “aging index” (or “C57 index”). Thus, when the parameter of interest was itself an index, we replaced the ratio (a – b)/(a + b) with the a/b ratio.
Probes and surgery
Every animal used in electrophysiological experiments (Extended Data Fig. 5-1) was implanted with a multi-shank silicon probe attached to a movable microdrive, equipped with optical fibers following previously described procedures (Stark et al., 2012; Noked et al., 2021). The probes used were Stark64 (Diagnostic Biochips), Buzaski32 (NeuroNexus), and Dual-sided64 (Diagnostic Biochips). The Stark64 probe consists of six shanks, spaced horizontally 200 μm apart, with each shank consisting of 10–11 recording sites spaced vertically 15 μm apart. The Buzaski32 probe consists of four shanks, spaced horizontally 200 μm apart, with each shank consisting of eight recording sites spaced vertically 20 μm apart. The Dual-sided64 probe consists of two dual-sided shanks, spaced horizontally 250 μm apart, with each shank consisting of 16 channels on each side (front and back), spaced vertically 20 μm apart.
Before probe implantation, the single-transgenic hybrid was injected with a DIO-hChR2 viral vector (rAAV5/EF1a-DIO-hChR2(H134R)-eYFP; 3.2 × 1012 IU/ml; University of North Carolina viral core facility; courtesy of K. Deisseroth). The solution was injected stereotactically (Kopf) into the neocortex and hippocampus at 8 different depths (AP −1.6, ML 1.1, DV 0.4–1.8 at 0.2 mm increments; 25 nl/site; Nanoject III, Drummond).
Probes were implanted in the parietal neocortex above the right hippocampus (AP/ML, −1.6/1.1 mm; 45° angle to the midline) under isoflurane (1%) anesthesia. Following recovery from anesthesia, linear-track animals were placed on a water-restriction schedule that guaranteed at least 40 ml/kg of water (corresponding to 1 ml/25 g mouse) on every recording day. Recordings were conducted 5 d/week, and animals received free water on the sixth day. After every one to five recording sessions, the probe was translated vertically downwards by up to 70 μm.
Histology
After recordings have ended, the implanted mice were deeply anesthetized with pentobarbital (100 mg/kg) and perfused with 0.1 m phosphate buffered saline (PBS; pH 7.4) and 4% paraformaldehyde (PFA). The brains were removed and postfixed overnight in PFA. Coronal sections (70 μm) were cut on a vibratome (VT1000S, Leica) and collected in PBS. Sections were mounted in Fluoromount with DAPI (F6057-20ML, Sigma) and imaged with a wide-field fluorescence microscope (Axio Scope A1, Zeiss).
Linear track sessions
For the linear track analyses, neuronal activity was recorded in 4.5 [1.9 10.1]-h sessions (median [interquartile interval, IQR]). At the beginning of every session, neural activity was recorded while the animal was in the home cage. The animal was then placed on a 150 cm linear track that extended between two 10 × 10 cm square platforms. Each platform included a water delivery port. Mice were under water restriction and were trained to repeatedly traverse the track for a water reward of 3–10 μl. Over all sessions, mice ran 167 [132 200] one-direction trials over about 1 h (Extended Data Fig. 5-2). Trials with a mean running speed below 10 cm/s were excluded from analyses. Animals were equipped with a three-axis accelerometer (ADXL-335, Analog Devices) for monitoring head movements. Head position and orientation were tracked in real-time using two head-mounted LEDs, a machine vision camera (ace 1300-1200uc, Basler), and a dedicated system (“Spotter,” Gaspar et al., 2019).
Spike detection and sorting
Neural activity was filtered, amplified, multiplexed, and digitized on the headstage (0.1–7500 Hz, x192; 16 bits, 20 kHz; RHD2132 or RHD2164, Intan Technologies) and then recorded by an RHD2000 evaluation board (Intan Technologies). Offline, spikes were automatically detected and sorted into single units using KlustaKwik3 (Kadir et al., 2014; Rossant et al., 2016) for shanks with up to 11 sites/shank or KiloSort2 (Pachitariu et al., 2016) for 16 channel shanks. Automatic spike sorting was followed by manual adjustment of the clusters. Only well-isolated units were used for further analyses [amplitude >40 μV; L-ratio <0.05 (Schmitzer-Torbert et al., 2005); ISI index <0.2 (Fee et al., 1996)]. Units were classified into putative pyramidal cells (PYR) or PV-like interneurons (INT) using a Gaussian mixture model (Stark et al., 2013).
Place field and phase precession analysis
Based on the linear track data, spatial information, place fields, and phase precession were determined for every PYR (Sloin et al., 2022). Briefly, place fields were defined as regions spanning 15–100 cm in which the firing rate increased compared with the on-track spontaneous firing rate (p < 0.05, Bonferroni-corrected Poisson test). Theta phase precession was quantified for each place field using a circular-linear analysis (Schmidt et al., 2009; Kempter et al., 2012). The circular-linear model yielded the precession slope, a, and the resultant length of the residuals, R, indicating model fit; statistical significance was determined by a permutation test (Sloin et al., 2022). Precession effect size was quantified as the ratio between the fit of spikes to the circular-linear model, R, divided by the median of 300 model fits to randomly permuted phase/position pairs.
Optogenetic stimulation
To determine the effect of optogenetic stimulation on spiking, 50-ms blue-light pulses were administered. In every session, illumination was conducted for every shank separately. In CA1, pulses were given in a total of 131 shanks during 55 sessions, using light power of 2.43 [0.96 3.98] μW. In a given stimulation experiment, pulses were applied 137 [51 209] times. In the neocortex (Extended Data Fig. 6-3), pulses were given in a total of 63 shanks during 21 sessions, using light power of 11.07 [5.25 21.71] μW. In a given stimulation experiment, pulses were applied 150 [75 300] times. Light-induced firing rate gain was defined as the mean firing rate during illumination, divided by the mean firing rate during baseline (in the lack of illumination on any shank). Units were determined as light-activated if the Poisson probability of seeing the observed number of spikes (or more) during illumination was < 0.05, based on the baseline firing.
Statistical analyses
In all statistical tests used in this study, a significance threshold of α = 0.05 was used. All descriptive statistics (n, median, IQR, mean, SEM) can be found in the results, the figure legends, and the tables. Since most tested parameters did not follow a standard normal distribution (Kolmogorov–Smirnov normality test), nonparametric statistical tests were used for all analyses. All statistical analyses were conducted in MATLAB (MathWorks). Differences between medians of two groups were tested with Mann–Whitney’s U test. Differences between medians of three groups or more were tested with Kruskal–Wallis nonparametric analysis of variance and corrected for multiple comparisons using Tukey’s procedure. Differences between medians measured along two dimensions (e.g., strain and day) were tested with a two-way Kruskal–Wallis analysis of variance. Wilcoxon’s signed-rank test was employed to determine whether a group median is distinct from zero. To estimate whether a given fraction was smaller or larger than expected by chance, an exact binomial test was used. Differences between the proportions of observations of two categorical variables were tested with a likelihood ratio (G-) test of independence. Bonferroni’s correction was employed in cases of G-test multiple comparisons.
Results
Hybrid mice are larger than C57 mice and exhibit similar sensorimotor performance
We decided to use offspring of C57 and FVB because the two inbred strains are of different genealogies (Beck et al., 2000) and exhibit distinct phenotypic disadvantages (Crawley, 2008). FVB dams produce large litters (Silver, 1995) and were therefore selected as the maternal strain. To maximize the genetic similarity between individual mice, we focused on (FVB/NJ x C57BL/6J)F1, the first-generation offspring (HYB). We found that sedentary HYB were larger than either parental strain. At the age of three months, the median [IQR] weight of HYB was 32 [29.6 34.9] g, compared with 26.9 [25.1 28.1] g for C57 (p = 9.6 × 10−10; Kruskal–Wallis test, corrected for multiple comparisons; Fig. 1; Extended Data Fig. 1-1A). Higher weights are beneficial for electrophysiological experiments in freely-moving animals, in which animal size limits the weight of the implanted apparatus. Examination of spontaneous home cage behavior (Extended Data Fig. 1-1B) and gait analysis (Extended Data Fig. 1-1C) revealed consistent differences between the three strains. Thus, the HYB constitute a unique, larger strain.
Hybrid mice weigh more than both progenitors. Mice of the three strains were weighted periodically at ages between 75 and 105 d (∼3-month-old; A) or between 255 and 285 d (∼9-month-old; B), and weights were averaged over measurements for every mouse. Every box plot shows median and interquartile range (IQR), whiskers extend for 1.5 times the IQR in every direction, and a plus indicates an outlier. P-value indicated by text is for a one-way Kruskal–Wallis test; ***p < 0.001, Kruskal–Wallis test, corrected for multiple comparisons. See also Extended Data Figures 1-1, 1-2, and 1-3.
Extended Data Figure 1-1
Hybrids are larger and exhibit distinct patterns of daily conduct in the home cage. A, HYB weigh more than both progenitors. Top, Mice of the three strains were weighed periodically over six months, from three to nine months old (n = 20 C57, 17 FVB, and 15 HYB). Each line connects measurements of a single animal. P-values indicated by text are for a strain by day two-way Kruskal–Wallis test. Dashed lines, age ranges used for statistical analyses (center and bottom panels). Every box plot shows median and interquartile range (IQR), whiskers extend for 1.5 times the IQR in every direction, and a plus indicates an outlier. For every mouse, weights were averaged over measurements made in the age of 75–105 d (“3 months”; center) and 255–285 d (“9 months”; bottom); ***p < 0.001, Kruskal–Wallis test, corrected for multiple comparisons. B, Hybrids are least active and spend more time in the shelter during the dark phase. Nine-month-old mice (n = 16 C57, 15 FVB, and 11 HYB) were placed in individual phenotyping cages and the distance covered, time in the shelter, feeding, and licking were sampled every hour over 3 d. a, Top, During both dark (8 A.M. until 8 P.M.) and light phases, HYB cover less distance. Bands show SEM. Second row, During the dark phase, HYB spend more time in the shelter than C57 or FVB. Third row, During both dark and light phases, HYB feed less than either C57 or FVB. Bottom, The number of licks is not consistently different between the strains. b, Mean metric (distance covered, sheltering, feeding/licking rate) during the dark (top half) and light (bottom) phases. Error bars, SEM; **p < 0.01, ***p < 0.001, Kruskal–Wallis test, corrected for multiple comparisons. c, Diurnal indices, defined as the difference between the metric (e.g., distance covered) during the light phase (L) and the dark phase (D), divided by the sum; *p < 0.05, ***p < 0.001, bootstrap test. Mice of all strains exhibit diurnal changes in all metrics tested. d, Pie charts show the fraction of time spent by mice of every strain carrying three specific behaviors: sheltering (being in the shelter), feeding (breaking the photobeam of the feeding station), and wheeling (being on the wheel). “Other” corresponds to home cage roaming, including licking. The time mice dedicated to each behavior differs between strains (***p < 0.001, strain, behavior, and interaction effects, two-way Kruskal–Wallis test). Compared to the dark phase, mice of all strains spend less time on the wheel and feeding during the light phase (***p < 0.001, Kruskal–Wallis test). C, Gait patterns of the three strains. a, Gait regularity does not differ consistently between strains or between age groups. Here and in all other panels, p-values indicated in text in the top and center are for Kruskal–Wallis test (e.g., 0.93). Bottom, Aging indices; *p < 0.05, ***p < 0.001 bootstrap test. b, The front paws base of support (BOS) of HYB increases with aging. c, The hind paws BOS is distinct for the three strains; *p < 0.05, **p < 0.01, ***p < 0.001, Kruskal–Wallis test, corrected for multiple comparisons. d–f, Track crossing duration, maximum variation, and cadence do not differ consistently between strains in either age group. Download Figure 1-1, EPS file.
Extended Data Figure 1-2
Results of phenotyping assays for every age group and strain. aValues indicate median [IQR] for every age group, strain, and metric. For learning on the MWM, values indicate mean ± SEM. Download Figure 1-2, DOCX file.
Extended Data Figure 1-3
Hybrid sensorimotor performance is maintained upon aging. A, Reaction to visual stimuli on the optomotor drum of HYB and C57 is similar and remains intact at an older age. In contrast, FVB do not react to rotating visual stimuli at either age group. Top row, Three-month-old mice. Center, Nine-month-old mice. Bottom, Aging indices (a) FVB mice do not make any head turns, indicating poor visual abilities. Here and in other panels, P-values indicated in text are for a Kruskal–Wallis test; **p < 0.01, ***p < 0.001, Kruskal–Wallis test, corrected for multiple comparisons. b, Directional bias in the optomotor drum. Every trial included 2 min of 2 rpm counterclockwise (CCW) and 2 min of 2 rpm clockwise (CW) drum rotations. Directionality indices (DI) were defined as the number of head turns during CCW versus CW drum rotations. c, Time spent making head movements. FVB mice do not move their heads as the drum turns, whereas the time spent moving the heads does not differ consistently between C57 and HYB in either age group. B, Motor endurance deteriorates for older C57 and FVB mice. Center, At the age of nine months, HYB cover longer distances than FVB; **p < 0.01, Kruskal–Wallis test, corrected for multiple comparisons. Bottom, In contrast to C57 and FVB, HYB maintain performance between the two age groups. Bars, mean aging index; error bars, SEM; **p < 0.01, ***p < 0.001, bootstrap test compared to a no-change null. C, Balance performance deteriorates upon aging for C57 and FVB. In contrast, HYB balance performance is maintained upon aging. Top and center, Mean (SEM) time to fall on five consecutive training days; *** next to the text indicate p < 0.001 for a two-way, strain by day, Kruskal–Wallis test. */*** at right: p < 0.05/p < 0.001, one-way Kruskal–Wallis test, measuring within-strain, across-day effect. In the older age group, HYB remained on the apparatus for the longest duration before falling. Bottom, Aging indices were computed for each strain every day; *p < 0.05, ***p < 0.001, bootstrap test. Inset, Aging indices, averaged over days. Download Figure 1-3, EPS file.
To assess the sensorimotor capabilities of HYB, we subjected three-month-old mice to three standard assays (n = 12 C57, 12 FVB, and 11 HYB; Table 1; Extended Data Fig. 1-2). Visual behavior was quantified in an optomotor drum test (Abdeljalil et al., 2005), in which mice follow vertical visual stimuli by making directional head movements. C57 and HYB made similar numbers of head turns (C57: 7.5 [3 11]; HYB: 7 [5 10], p = 0.98, Kruskal–Wallis test; Extended Data Fig. 1-3A). In contrast, FVB did not make any head turns (0 [0 0]; C57 vs FVB: p = 0.0011; HYB vs FVB: p = 0.003), consistent with blindness because of retinal degeneration (Pittler and Baehr, 1991). Motor endurance was tested by running the mice on a treadmill at gradually increasing speed, and measuring the distance to failure (Marques-Aleixo et al., 2015). HYB stayed on the treadmill for 222 [30 262] s, not consistently different from C57 (227 [148 242] s; p = 0.21; Kruskal–Wallis test) or from FVB (120 [73 197] s; p = 0.65; Extended Data Fig. 1-3B). Balance was assessed by placing the mice on an accelerating rotarod for five consecutive days, and measuring the latency to fall (Jones and Roberts, 1968). C57 and FVB exhibited increased latency to fall over days, indicative of balance learning (C57: p = 0.03, FVB: p = 0.03, within-strain Kruskal–Wallis test; Extended Data Fig. 1-3C). HYB latency to fall did not increase over days (p = 0.12, Kruskal–Wallis test). To summarize, at the age of three months, visual behavior, and motor endurance of HYB are similar to those of C57.
Hybrid mice exhibit reduced anxiety and improved learning and memory
To characterize neuropsychiatric phenotypes and cognitive performance, we tested three-month-old mice using five assays (Table 1). Anxiety-like behavior was studied using the elevated plus maze (Rodgers and Dalvi, 1997) and quantified by a “contra-anxiety index” that ranges from −1 to 1, defined as time spent in (open – closed)/(open + closed) arms. Indices closer to one indicate that more time was spent in the two “open” arms, corresponding to reduced anxiety-like behavior. C57 spent consistently more time in the closed arms compared with the open arms (−0.59 [−0.66 −0.55]; p = 0.016, Wilcoxon’s signed-rank test), as did HYB (−0.34 [−0.55 −0.22], p = 0.016; Fig. 2A; Extended Data Fig. 2-1A). However, HYB spent more time in the open arms than C57, indicating reduced anxiety-like behavior (p = 0.034; Kruskal–Wallis test; Fig. 2A; Extended Data Fig. 1-2).
Young hybrids exhibit reduced anxiety. A, At the age of three months, HYB exhibit reduced anxiety-like behavior compared with C57. Anxiety-like behavior was tested in an elevated plus maze and quantified using a “contra-anxiety index,” defined as the time spent in the open minus the time spent in the closed arms, divided by the sum. Lined *p < 0.05, Kruskal–Wallis test. Here, and in C, *p < 0.05, **p < 0.01, Wilcoxon’s signed-rank test, comparing to chance level (horizontal dashed line). B, HYB exhibit improved learning in the MWM task compared with C57 and FVB, as measured by the latency to platform (a) and by the distance traveled (b) over days. Bands show mean and SEM; *** next to the top text indicate p < 0.001 for a two-way, strain by day, Kruskal–Wallis test; */**/*** next to strain names indicate p < 0.05/p < 0.01/p < 0.001 for a within-strain, across-day one-way Kruskal–Wallis test. C, Spatial long-term memory behavior on the MWM. a, When performance is quantified by the fraction of time spent in the target quadrant, there is no consistent difference between the strains. b, When performance is quantified by the fraction of visits to the target quadrant, HYB exhibit improved performance compared with C57 and to FVB. All other conventions are the same as in A. D, Compared with C57, three-month-old HYB exhibit improved anxiety-like behavior, improved exploration, and decreased balance performance, as quantified by the “C57 index.” Positive Indices indicate improved performance compared with C57. Circles (and lines) at vertices show mean (and SEM) indices; filled circles indicate significant changes compared with C57 (p < 0.05, bootstrap test). See also Extended Data Figure 2-1.
Extended Data Figure 2-1
Hybrids exhibit reduced anxiety-like behavior, improved learning, and improved memory. A, HYB explore the open arms of an elevated plus maze more than C57 mice, particularly at an older age, indicating reduced anxiety-like behavior. Here and in subsequent panels, Top and center, Every box plot shows median and interquartile range (IQR), whiskers extend for 1.5 times the IQR in every direction, and a “+” sign indicates an outlier. Here and in Eb, Ec, P-values indicated by text are for a one-way Kruskal–Wallis test; lined *p < 0.05, **p < 0.01, ***p < 0.001, Kruskal–Wallis test, corrected for multiple comparisons; *p < 0.05, **p < 0.01, ***p < 0.001, Wilcoxon’s signed-rank test, comparing indices to chance level. Bottom, Aging index. Bars and errors show mean and SEM, respectively; ***p < 0.001, bootstrap test, compared to a no-change null (index of zero). B, Depressive-like behavior is not consistently different among HYB and C57 mice, and is reduced among older FVB mice. C, Open field exploration increases at an older age for all strains. Top, In the three-month-old age group, HYB were more active in the open field than FVB; **p < 0.01, Kruskal–Wallis test, corrected for multiple comparisons. Center, In the nine-month-old group, C57 were more active than FVB. Bottom, Mice of all three strains exhibited increased activity upon aging, evident as above-zero aging indices; **p < 0.01, ***p < 0.001, bootstrap test. D, Learning on the MWM task is improved for HYB compared to C57, and is maintained at older age. a, Mean (and SEM) latency to finding the platform (limited to 60 s). Here and in c, **/*** next to the text indicate p < 0.01/p < 0.001 in a two-way, strain by day, Kruskal–Wallis test; */**/*** at right indicate p < 0.05/p < 0.01/p < 0.001, one-way Kruskal–Wallis test, measuring the within-strain, across-day effect. In both age groups. only HYB showed consistent across-day learning. Bottom, Aging indices, averaged over all learning days. Learning performance of C57 deteriorates at an older age. b, Day-by-day “C57 indices” (Fig. 1D), measuring the change in latency to the platform between FVB (blue) or HYB (green) and C57 mice. Bands show mean and SEM; *p < 0.05, **p < 0.01, ***p < 0.001, bootstrap test. Insets, Mean (and SEM) C57 indices, averaged over all learning days; p-values, geometric mean over all learning days. c, C57 and HYB travel progressively shorter distances over training days, indicating learning. All conventions are the same as in a. E, On the probe day of the MWM task, HYB spend the longest time in and the shortest latency to the target quadrant. a, Time spent in every quadrant; *p < 0.05, **p < 0.01, ***p < 0.001, Wilcoxon’s signed-rank test, comparing time in each quadrant to chance level (0.25 of the total time). b, Time spent in the target quadrant, for three-month-old (top) and nine-month-old (center) mice. Here and in c, Horizontal dashed lines indicate chance level of 0.25, expected for random exploration; other conventions are the same as in A. Compared to C57 or FVB, nine-month-old HYB spend a larger fraction of time in the target quadrant. Bottom, Aging indices. Compared to younger HYB, older HYB spend more time in the target quadrant. c, In both age groups, the fraction of visits to the target quadrant is highest for HYB. d, Latency to the platform location, measured during the probe day, is the shortest for HYB. Download Figure 2-1, EPS file.
Depressive-like behavior was tested using a forced swim test, in which longer freeze duration indicates behavioral despair (Porsolt et al., 1977). Mice of all three strains exhibited similar freeze durations (C57: 124 [94 135] s; FVB: 209 [130 261] s; HYB: 201 [85 290] s; p = 0.23, Kruskal–Wallis test; Extended Data Fig. 2-1B). Motivational behavior was quantified by activity levels in an open field (Christmas and Maxwell, 1970). We did not observe consistent differences between C57 and HYB in overall activity in the field (C57: 0.23 [0.17 0.27] %; HYB: 0.28 [0.2 0.31] %; p = 0.36, Kruskal–Wallis test; Extended Data Fig. 2-1C). Thus, compared with C57, HYB exhibit similar depressive-like behavior and similar motivational behavior.
To assess learning and memory behavior, we trained three-month-old mice on the MWM task (Morris, 1981; n = 10 C57, 9 FVB, and 11 HYB; Table 1). During four consecutive days, mice were placed in a pool with an invisible submerged platform. Over days, only HYB exhibited a decrease in the latency to platform, indicative of learning (C57: p = 0.67, FVB: p = 0.08, HYB: p = 0.03, intrastrain Kruskal–Wallis tests; Fig. 2B). Pooled over days, the mean ± SEM latency to the platform was shorter for HYB (19.3 ± 3.9 s), compared with C57 (21.7 ± 4.3 s) and to FVB (41.1 ± 3.8 s; p = 1.7 × 10−13, Kruskal–Wallis test; Fig. 2Ba). Similar results were obtained when examining the distance traveled (Terry, 2009) instead of the latency to platform (Fig. 2Bb). Hence, compared with both progenitors, three-month-old HYB exhibit improved learning performance.
To assess spatial long-term memory, the mice were tested in the MWM during a “probe” day, without a platform. Mice of the three strains spent similar fractions of time in the target quadrant (p = 0.13, Kruskal–Wallis test; Fig. 2Ca). However, when spatial long-term memory was quantified by the fraction of visits to the target quadrant, HYB exhibited the largest fraction of visits to the target quadrant, indicative of spatial long-term memory (C57: 0.28 [0.26 0.27]; FVB: 0.21 [0.15 0.3]; HYB: 0.36 [0.33 0.44]; p = 0.0014, Kruskal–Wallis test; Fig. 2Cb). Thus, in contrast to both progenitors, three-month-old HYB exhibit spatial long-term memory behavior.
To quantify differences between HYB (or FVB) and the C57, we treated the C57 strain as a reference and computed “C57 indices” for the eight assays (Fig. 2D). For a given metric, the C57 index is defined as the difference between the median value of HYB (or FVB) and the median value of C57, divided by the sum. Therefore, a C57 index above zero indicates higher values of the metric for the tested strain, compared with C57 performance. HYB performance was similar to C57 in assays testing visual behavior (mean ± SEM C57 index: −0.042 ± 0.19, p = 0.36, bootstrap test), motor endurance (−0.06 ± 0.2, p = 0.47), depressive-like behavior (0.21 ± 0.23, p = 0.19), learning (day-averaged latency to platform; 0.05 ± 0.2, p = 0.25), and spatial long-term memory (time in target quadrant; −0.011 ± 0.08, p = 0.47). In contrast, HYB exhibited reduced anxiety-like behavior (0.31 ± 0.1, p = 0.011) and increased exploratory behavior (0.089 ± 0.05, p = 0.031). Thus, at the age of three months, HYB motor endurance and visual behavior is similar to C57, whereas the HYB exhibit reduced anxiety-like behavior.
Older hybrids maintain sensorimotor performance and exhibit improved memory
To examine the suitability of the HYB for long-term studies that require stable behavioral performance, we applied all assays also to nine-month-old mice (Table 1). Akin to the three-month-old HYB, nine-month-old HYB exhibited reduced anxiety-like behavior compared with C57 (contra-anxiety index: C57: −0.56 [−0.75 −0.46]; HYB: −0.15 [−0.27 0.15]; p = 4.1 × 10−6; Kruskal–Wallis test; Fig. 3A; Extended Data Fig. 2-1A). The nine-month-old HYB exhibited improved learning compared with C57, as indicated by shorter latency to platform (C57: 33.3 ± 2.9 s; HYB: 17.8 ± 3 s; p = 1.7 × 10−13, Kruskal–Wallis test; Fig. 3B; Extended Data Fig. 2-1D). Furthermore, nine-month-old HYB displayed improved spatial long-term memory compared with C57, as indicated by a larger fraction of time spent in the target quadrant (C57: 0.3 [0.24 0.4]; HYB: 0.55 [0.36 0.62]; p = 0.002, Kruskal–Wallis test) and by a larger fraction of visits to the target quadrant (C57: 0.35 [0.27 0.4]; HYB: 0.45 [0.36 0.5]; p = 0.013; Kruskal–Wallis test; Fig. 3C; Extended Data Fig. 2-1E). Correspondingly, the C57 indices of HYB indicated reduced anxiety-like behavior (mean ± SEM C57 index: 0.49 ± 0.1, p = 0.00,001), improved learning (day-averaged latency to platform; 0.35 ± 0.14, p = 0.002), and improved spatial long-term memory (time in target quadrant; 0.3 ± 0.06, p = 0.0003; Fig. 3D). The performance of nine-month-old HYB was similar to C57 in assays testing visual behavior (−0.031 ± 0.1, p = 0.34, bootstrap test), motor endurance (0.14 ± 0.12, p = 0.09), balance (0.04 ± 0.05, p = 0.14), exploration (−0.01 ± 0.023, p = 0.35), and depressive-like behavior (0.04 ± 0.08, p = 0.32). Thus, nine-month-old HYB exhibit sensorimotor properties which are similar to nine-month-old C57, while exhibiting reduced anxiety-like behavior and improved learning and memory performance.
Older hybrids exhibit improved learning and long-term memory. A, At the age of nine months, HYB exhibit reduced anxiety-like behavior compared with C57. All conventions are the same as in Figure 2A. B, Nine-month-old HYB exhibit improved learning in the MWM task compared with C57 and FVB, indicated by shorter latency to platform (a) and shorter distance travelled (b) over days. All conventions are the same as in Figure 2B. C, Compared with C57 and FVB, nine-month-old HYB exhibit improved spatial long-term memory, indicated by larger fraction of time spent in the target quadrant (a) and by a larger fraction of visits to the target quadrant (b) during the probe day. All conventions are the same as in Figure 2C. D, Nine-month-old HYB exhibit improved anxiety-like behavior, improved learning, and improved spatial long-term memory behavior compared with C57. All conventions are the same as in Figure 2D.
To quantify performance differences between three- and nine-month-old mice of the same strain, we defined an “aging index” (Fig. 4). For a given performance metric, the aging index is the difference between the median value obtained for nine- and three-month-old mice, divided by the sum. Thus, an aging index of zero indicates no age-related changes in performance, whereas positive indices indicate age-dependent increases. Visual behavior was maintained across age groups for all strains (mean ± SEM aging indices: C57, 0.03 ± 0.17, p = 0.41; FVB, not available; HYB, 0.04 ± 0.13, p = 0.28, bootstrap test; Fig. 4Aa). While C57 and FVB exhibited reduced motor endurance on aging (C57: −0.41 ± 0.1, p = 0.0001; FVB: −0.5 ± 0.13, p = 0.0014), HYB endurance was maintained (−0.23 ± 0.22, p = 0.11; Fig. 4Ab). Finally, while C57 and FVB balance performance declined between age groups (C57: −0.2 ± 0.06, p < 0.05; FVB: −0.31 ± 0.09, p < 0.05), HYB balance was maintained (0.1 ± 0.13, p = 0.75; Fig. 4Ac). Thus, while the motor capabilities of C57 and FVB deteriorate upon aging, HYB motor performance remains stable across the two age groups studied.
Hybrids maintain sensorimotor performance and exhibit improved long-term memory on aging. A, Sensorimotor behavior remains intact at an older age for HYB but not for C57 or FVB mice. Positive aging indices indicate improvement across three- and nine-month-old mice of the same strain. a, Visual behavior of HYB and C57 remains intact at an older age, whereas FVB do not react to rotating visual stimuli at any age. Here and in B, bars indicate mean; error bars, SEM; *p < 0.05, **p < 0.01, ***p < 0.001, bootstrap test compared with a no-change null. b, Motor endurance is reduced for older C57 and FVB and maintained for older HYB. c, Balance performance deteriorates for C57 and FVB on aging and is maintained for older HYB. B, Psychiatric-like and cognitive performance across age groups. All conventions are the same as in A. a, HYB, but not C57 or FVB, exhibit reduced anxiety-like behavior upon aging. b, Learning performance of C57, but not FVB or HYB, deteriorates upon aging. c, Spatial long-term memory performance of C57 and FVB is maintained upon aging, whereas HYB memory performance improves. C, Comparison of performance changes in eight behavioral assays using the “aging index.” All conventions are the same as in Figure 2D.
Compared with C57, three-month-old HYB exhibited reduced anxiety-like behavior, improved learning, and improved long-term memory (Fig. 2). We tested whether performance on the three assays was maintained across age groups in every strain. Compared with three-month-old mice of the same strain, anxiety-like behavior was reduced for nine-month-old HYB (mean ± SEM aging index: 0.26 ± 0.11, p = 0.0005, bootstrap test), but not for C57 (0.04 ± 0.12, p = 0.37) or FVB (0.22 ± 0.17, p = 0.07; Fig. 4Ba). Learning performance was reduced in older C57 mice (−0.28 ± 0.14, p = 0.03), but maintained for older HYB (0.04 ± 0.21, p = 0.43; Fig. 4Bb). Both C57 and FVB maintained spatial long-term memory performance between age groups (C57: 0.008 ± 0.13, p = 0.28; FVB: −0.07 ± 0.16, p = 0.16). In contrast, older HYB exhibited improved spatial long-term memory (0.31 ± 0.07; p = 0.0003; Fig. 4Bc). Thus, anxiety-like behavior is reduced and memory performance improves on aging in HYB (Fig. 4C).
Hybrid mice exhibit enhanced linear track running and typical CA1 place coding
The results reported in Figures 1-4 show that HYB exhibit certain advantages over C57. However, improved performance on behavioral assays spanning several minutes does not necessarily predict enhanced behavior over daily recording sessions that span hours and include the additional weight of an implanted apparatus. To assess the suitability of freely-moving HYB for electrophysiological experiments, we performed high-density recordings from two HYB and two C57 mice running back and forth on a 150-cm linear track (Extended Data Fig. 5-1). Behavior during the first 14 sessions on the track was assessed to quantify performance at similar experience levels. HYB conducted a median [IQR] of 164 [139 215] one-direction trials per session, while C57 ran 136 [113 178] trials (p = 0.015, Mann–Whitney test; Fig. 5A; Extended Data Fig. 5-2). Furthermore, HYB exhibited higher running speeds, 31 [22 44] cm/s compared with 23 [18 30] cm/s for C57 (p = 0.007, Mann–Whitney test). Thus, freely-moving HYB exhibit enhanced motor abilities and exploration during prolonged electrophysiological recordings.
Hybrid mice exhibit enhanced linear track running and typical CA1 place coding. A, Compared with C57 mice, HYB carry out more trials and run faster on a 150-cm-long linear track. Only the first 14 sessions from every mouse were employed (two C57 and two HYB mice); *p < 0.05, **p < 0.01, Mann–Whitney test. Dashed vertical lines indicate group medians. a, The number of one-direction trials. b, Mean running speed over trials. B, Location of pyramidal cell (PYR) somata relative to the center of the CA1 pyramidal cell layer (dashed horizontal line) is not consistently different between strains. The estimated depths of PYR recorded by high-density silicon probes were partitioned into 5-μm bins. Positive numbers correspond to PYR closer to str. oriens. n.s.: p > 0.05, Mann–Whitney test. C, Units recorded from both strains exhibit increased firing rates within specific regions of the linear track. Each row represents a unit; firing rates on right (R) to left (L) runs are concatenated with L to R runs and scaled to the 0–1 (white-black) range for presentation purposes. Bar graphs (right) show peak on-track firing rates. D, A PYR recorded from a HYB exhibits a typical place field and theta phase precession. Top, Firing rate as a function of position (mean ± SEM over 82 right-to-left trials). Bottom, Theta phase and animal position at the time of every spike; phase of zero corresponds to theta peak. Running direction is presented from left to right, and vertical dashed lines indicate place field limits. a represents phase precession slope, and R represents the goodness of fit of spikes to the precession model; **p < 0.01, permutation test. Left, Wide-band (0.1–7500 Hz) spike waveforms (mean ± SD) recorded on eight consecutive sites, vertically spaced by 20 μm. Bottom left, Autocorrelation histogram. E, HYB and C57 exhibit similar CA1 spatial rate coding. a, Spatial information rate. n.s.: p > 0.05, Mann–Whitney test. Here and in c, vertical dashed lines indicate group medians. b, Fraction of CA1 PYR with one or more place fields out of all units active and stable on the track; n.s.: p > 0.05, likelihood ratio test. Here and in Fa, ***p < 0.001, exact binomial test, compared with chance level (0.05; horizontal dashed line). Error bars, SEM. c, Place field size; n.s.: p > 0.05, Mann–Whitney test. F, HYB and C57 exhibit similar CA1 spatial phase coding. a, Fraction of fields exhibiting theta phase precession; n.s.: p > 0.05, likelihood ratio test. b, Precession slope size. c, Precession effect size; n.s.: p > 0.05, Mann–Whitney test. See also Extended Data Figures 5-1 and 5-2.
Extended Data Figure 5-1
Mice used for electrophysiological recordings. aFor C57 animals, the male parent was the reporter, and the female parent was the driver. bAt time of implantation. Download Figure 5-1, DOCX file.
Extended Data Figure 5-2
Linear track behavior and place coding in individual mice. aNumber of one direction trials, median [IQR] over sessions. bWell-isolated PYR, active and stable on the linear track. cSI, spatial information. Download Figure 5-2, DOCX file.
To examine whether CA1 pyramidal neurons (PYR) exhibit similar place coding in HYB and C57 mice tested under the same conditions, we analyzed the activity of CA1 pyramidal layer PYR active on the track (C57: n = 779, HYB: n = 960; Fig. 5B; Extended Data Fig. 5-2). In both strains, most active PYR exhibited increased firing rates in specific parts of the track (Fig. 5C). An example PYR recorded in HYB exhibited typical spatial rate coding (place field; O’Keefe and Dostrovsky, 1971; Moser et al., 2008) and spatial phase coding (theta phase precession; O’Keefe and Recce, 1993; Fig. 5D). At the group level, there was no consistent difference in the spatial information carried by PYR, with 1.28 [0.47 3.18] bits/s for HYB PYR and 1.37 [0.63 2.87] bits/s for C57 PYR (p = 0.56, Mann–Whitney test; Fig. 5Ea). Place fields were detected in 443/779 (57%) of C57 PYR and in 508/960 (53%) of HYB PYR (Fig. 5Eb). Both fractions were above expected by chance (p < 1.11 × 10−16, exact binomial test) and were not consistently different from each other (p = 0.37, likelihood ratio test). Furthermore, place field sizes were not consistently different, spanning 40 [27.5 57.5] cm for C57 PYR and 40 [25 57.5] cm for HYB PYR (603 C57 fields and 649 HYB fields; p = 1, Mann–Whitney test; Fig. 5Ec). Thus, CA1 PYR recorded in HYB and in C57 mice exhibit similar rate coding of space.
To determine whether HYB CA1 PYR exhibit typical phase coding of space, we examined theta phase precession within PYR place fields. Of all fields, 308/603 (51%) C57 PYR and 353/649 (54%) HYB PYR exhibited consistent spatial theta phase precession (p < 1.11 × 10−16, exact binomial test; Fig. 5Fa). The fraction of place fields with precession was not consistently different between the two strains (p = 0.51, likelihood ratio test). Precession slope size was not consistently different between strains, being 0.02 [0.01 0.03] cyc/cm for C57 PYR place fields and 0.02 [0.02 0.03] cyc/cm for HYB PYR place fields (p = 0.071, Mann–Whitney test; Fig. 5Fb). Finally, precession effect size was not consistently different between strains, with 1.6 [1.28 2.13] for C57 phase precessing fields and 1.54 [1.24 2.06] for HYB fields (p = 0.29, Mann–Whitney test; Fig. 5Fc). Thus, CA1 PYR recorded in HYB and in C57 exhibit similar phase coding of space.
Transgenic hybrid mice enable optogenetic control of individual neurons
A major advantage of C57 mice for neurobehavioral research is the myriad transgenic strains available on the same background. To combine genetic control available for C57 with the enhanced behavior of the HYB, we first generated dual-transgenic mice on a C57 background by breeding CaMKII-Cre females with Ai32 males. We then crossed the male offspring with FVB females. Opsin expression in the resulting dual-transgenic and hybrid mice was verified using histology (Extended Data Fig. 6-1A). The dual-transgenic HYB are susceptible to blue-light optogenetic activation of projection neurons, including neocortical and CA1 PYR. Optogenetic PYR activation was conducted in transgenic HYB implanted with multi-site optoelectronic probes (Extended Data Fig. 5-1). A total of 5568 PYR and 1066 putative interneurons (INT) were recorded from the CA1 pyramidal layer during n = 55 sessions (Extended Data Fig. 6-2). In every recording session, we tested the response of the recorded units to 50-ms blue-light pulses. During illumination, PYR exhibited increased firing rates (Fig. 6A, purple). Of all recorded PYR, 1213 (22%) exhibited a consistent increase in firing rate during illumination (p < 1.11 × 10−16, exact binomial test).
Transgenic hybrid mice enable optogenetic control of individual CA1 neurons. A, Directly activated PYR and an indirectly activated interneuron (INT) recorded in hippocampal region CA1 of a HYB dual-transgenic (CaMKII::ChR2) mouse. a, The units exhibit a cross-correlation histogram (center; no light condition) consistent with monosynaptic excitation (p < 0.001, Poisson test). b, Local field and spiking responses to a single 50-ms blue-light pulse (0.67 μW). The PYR and the INT both spike during the light. c, Raster plots (top) and PSTHs (bottom) during blue light pulses for the PYR (left) and the INT (right). For visualization purposes, raster plots include 200 pulses. d, Top parts, Mean PSTHs of 22 PYR (purple) and four INT (green) recorded simultaneously on the illuminated shank (S1). Bottom parts, Greyscale, PSTHs of all 48 PYR and eight INT recorded simultaneously during the session. Each row shows the PSTH of one unit, scaled to the 0–1 (white-black) range for visualization purposes. Arrows indicate the example PYR and INT. B, Activation probability of CA1 units depends on the distance from the illuminated shank. Dataset includes a total of 5568 PYR and 1066 INT. a, Bars show the fraction of units with increased firing rate (p < 0.05, Poisson test) during single-shank illumination. Error bars, SEM; ***p < 0.001, exact binomial test, comparing the fraction of light-activated units to chance level (horizontal dashed line). b, Median firing rate gain of activated PYR and INT. Error bars, IQR. Gain is shown only for distances with an above-chance number of optically activated units. P-values indicated in text are for a one-way Kruskal–Wallis test; ***p < 0.001, Kruskal–Wallis test. See also Extended Data Figures 6-1, 6-2, and 6-3.
Extended Data Figure 6-1
Opsin expression in transgenic and injected hybrids. A, Opsin expression in neocortex and hippocampus of mA234, a double transgenic and hybrid mouse (FVB/NJ x CaMKII-Cre x Ai32) used in linear track (Fig. 5) and optogenetic (Fig. 6) experiments. EYFP fluorescence (green) is used as a reporter for ChR2 expression. Left, Low-magnification (2.5×) image of neocortex and hippocampus of a transgenic mouse. Right, High-magnification (20×) image, showing EYFP in neocortical somata. B, Opsin expression in neocortex and hippocampus of a HYB injected with a viral vector (rAVV5/CaMKII-hChR2(H134R)-mCherry) in the right hemisphere. mCherry fluorescence (red) is used as a reporter for ChR2 expression. All other conventions are the same as in A. Download Figure 6-1, EPS file.
Extended Data Figure 6-2
Optogenetic activation of individual neurons in individual mice. aWell-isolated PYR with increased firing rate during optogenetic stimulation (p < 0.05, Poisson test), out of all PYR recorded in the same sessions. bLight-activated PYR recorded on the illuminated shank. Download Figure 6-2, DOCX file.
Extended Data Figure 6-3
Transgenic hybrids support robust optogenetic activation of individual neocortical neurons. A, B, Optogenetic activation in the neocortex of hybrid freely-moving dual-transgenic CaMKII::ChR2 mice. A, Example directly activated PYR and indirectly activated INT recorded in HYB neocortex. a, The units exhibit cross-correlation histogram (center; no light condition) consistent with monosynaptic excitation (p < 0.001, Poisson test). Bottom, Auto-correlation histograms (ACHs). Wide-band (0.1–7500 Hz) spike waveforms (mean ± SD) recorded on the same shank; intersite vertical spacing is 20 μm. b, Local field and spiking responses of the units to a single 50-ms blue-light pulse (1.71 μW). Both units spike during the light. c, Raster plots (top) and PSTHs (bottom) during 230 blue light pulses for the PYR (left) and indirectly activated INT (right). The raster plot shows 200 representative repetitions. d, Top, Mean PSTHs of 26 PYR (purple) and seven INT (green) recorded simultaneously on the illuminated shank during light pulses. Bottom, Greyscale, PSTH of all unit recorded during the session (49 PYR and 21 INT). S1, S2, S3, and S4 denote different shanks. Each row corresponds to the PSTH of one unit, scaled to the 0–1 (white-black) range for visualization purposes. Arrows indicate the example units. B, Activation probability and firing rate gain of neocortical units depend on the distance from the illuminated shank. Dataset includes a total of 1546 PYR and 517 INT. a, Bars show the fraction of units with increased firing rate (p < 0.05, Poisson test) during single-shank illumination as a function of the horizontal distance from the illuminated shank. Error bars, SEM; ***p < 0.001, exact binomial test, compared to chance level (0.05, horizontal dashed line). Neocortical PYR and INT recorded up to 600 μm away from the illuminated shank are activated. b, Median firing rate gain of activated PYR and INT as a function of distance from the illuminated shank. Error bars, IQR. Gain is considered only for distances with an above-chance number of activated units. P-values indicated in text are for a one-way Kruskal–Wallis test; ***p < 0.001, Kruskal–Wallis test, corrected for multiple comparisons. Gain differs between activated PYR recorded on the illuminated shank and activated PYR recorded 400 μm away. C, Parvalbumin-immunoreactive (PV) cell activation in a hybrid freely-moving single-transgenic PV-Cre mouse, injected with a Cre-dependent viral vector for ChR2. a, Subnetwork of a PYR (PYR1; purple) and a PV cell (green), recorded simultaneously using a high-density optoelectronic probe. ACHs and cross-correlations (no light condition), consistent with reciprocal monosynaptic excitation and inhibition (p < 0.001, Poisson test). b, Wide-band (0.1–7500 Hz) traces and spikes of the PV cell and another unit, PYR2. Top three traces are neuronal recordings; fourth trace indicates instantaneous light power emitted by the fiber attached to the recording shank; and the bottom trace indicates head acceleration. During blue light pulses (62 μW), the PV cell spikes, whereas PYR2 spiking is suppressed. c, Raster plot (top) and PSTH (bottom) during blue light pulses for the light-activated PV cell (left) and the suppressed PYR (right). Arrows indicate the example units. d, PSTHs of six units recorded by the illuminated shank during blue light pulses. Bottom, Each row corresponds to the PSTH of one unit, scaled to the 0–1 (white-black) range for visualization purposes. Top, Mean PSTHs. Download Figure 6-3, EPS file.
As previously reported in rat hippocampus (Stark et al., 2012), activation probability depended on the horizontal distance from the light source. Of the PYR recorded on the illuminated shank, 772/1527 (47%) exhibited optical activation (p < 1.11 × 10−16, exact binomial test), compared with 340/1984 (17%) units recorded on the adjacent shank (horizontal spacing, 200 μm; p = 2.5 × 10−13, exact binomial test; Fig. 6Ba, left). The fraction of optically-activated PYR was lower on shanks farther away (0 vs 400 μm away: p < 1.1 × 10−16, Bonferroni-corrected likelihood ratio test). Moreover, the firing rate gain of the PYR that responded to focal illumination depended on the distance from the illuminated shank, being 6.5 [3.59 13.97] for same-shank PYR, compared with 4.1 [2.49 7.35] for PYR recorded on the adjacent shank (p = 9.6 × 10−10, Kruskal–Wallis test; Fig. 6Bb, left). Therefore, CA1 PYR in dual-transgenic HYB exhibit robust distance-dependent responses to focal illumination.
Optogenetic stimulation led to the indirect activation of INT via monosynaptic inputs from excitatory neurons (Fig. 6A, green). A total of 241/1066 (23%) of CA1 INT exhibited light-induced activation (p < 1.11 × 10−16, exact binomial test; Fig. 6Ba, right). The fraction of light activated INT was 34% (92/273) on the illuminated shank and lower on shanks farther away (400 μm: 40/258, 16%, p = 1.3 × 10−4, Bonferroni-corrected likelihood ratio test). In contrast to PYR, the firing rate gain of light-activated INT did not differ consistently between shanks (p = 0.38, Kruskal–Wallis test; Fig. 6Bb, right). Thus, consistent with previous observations in C57 mice (Stark et al., 2013), indirect INT activation in CA1 leads to more widespread activation than the local activation achieved in PYR.
To examine HYB optogenetic activation in a brain region with a distinct architecture, we performed optical stimulation in the neocortex of the same mice (n = 21 sessions; Extended Data Fig. 6-3). Similar to CA1 units, neocortical PYR and INT exhibited induced spiking (Extended Data Fig. 6-3A), with activation probability that decreased with distance from the illuminated shank (Extended Data Fig. 6-3Ba). A total of 274/409 (67%) neocortical PYR recorded on the illuminated shank exhibited optical activation, compared with 722/1527 (47%) CA1 PYR (p = 1.2 × 10−4, likelihood ratio test). Furthermore, 331/591 (56%) neocortical PYR exhibited optical activation on the adjacent shank (200 μm away), compared with 340/1984 (17%) CA1 PYR (p < 3.5 × 10−11, likelihood ratio test). Compared with CA1, focal illumination in the neocortex activated PYR located farther away from the illuminated shank. Light-induced firing rate gain did not differ between the illuminated shank (7.6 [4.54 16.03]) and the adjacent shank (7.6 [4.08 13.94]; p = 0.79, Kruskal Wallis test; Extended Data Fig. 6-3Bb, left). The difference between brain regions is consistent with different network topologies, since in contrast to the neocortex, CA1 PYR do not exhibit abundant recurrent excitation (Thomson and Radpour, 1991; Buhl and Whittington, 2007). Neocortical units on distant shanks may be activated indirectly via excitatory synaptic connections from directly activated PYR closer to the illumination source.
An alternative approach for establishing optogenetic control in HYB is to cross FVB females with single-transgenic driver male mice and to inject the offspring with a viral vector that expresses a reporter gene. To demonstrate feasibility, we generated a PV-Cre x FVB mouse and injected the animal with a viral vector that allows Cre-dependent expression of ChR2. Focal illumination in the transgenic HYB induced direct PV activation and indirect PYR silencing (Extended Data Fig. 6-3C). In summary, crossing transgenic C57 with FVB yields transgenic HYB offspring, suitable for optogenetic experiments.
Discussion
Compared with the parental C57BL/6J progenitor strain, young first-generation hybrid offspring (FVB/NJ x C57BL/6J)F1 exhibited reduced anxiety-like behavior. Upon aging, HYB anxiety-like behavior was further reduced, and learning and long-term memory performance improved. In contrast, learning and memory performance of the parental C57BL/6J and the maternal FVB/NJ progenitor strains did not improve at older age. Furthermore, motor performance of the progenitor strains deteriorated upon aging, whereas HYB sensorimotor abilities were maintained. HYB were larger, ran faster, and performed more trials during electrophysiological experiments on the linear track. HYB and C57 exhibited similar CA1 place cell rate and phase coding. Finally, optogenetic manipulations were readily achieved in transgenic HYB.
Mice have emerged as a dominant model for biomedical research, mainly because of their small size, high reproductive rate, and genetic similarities to humans. The same properties facilitated the engineering of multiple transgenic lines. The maintenance of transgenic lines, particularly the generation of double-transgenic or multitransgenic mice, is greatly facilitated by inbred mice (Silva et al., 1997). However, inbred mice are suboptimal for neurobehavioral studies because of idiosyncratic recessive traits. Furthermore, C57 exhibit age-dependent deterioration, limiting the duration of prolonged experiments. C57 mice are also smaller than other strains, which is a disadvantage for electrophysiological studies involving implanted headgear. Since headgear weight is limited by body size, the range of studied behaviors is limited by animal weight. A single well-defined outbreeding step produces larger offspring that exhibit enhanced behavior even during electrophysiology experiments, maintaining the ability of genetic targeting.
One of the challenges in mouse-based research is adjusting the strain to the specific study question. Strain selection can potentially affect research outcome, since strain characteristics may interact with genetic manipulations (Sittig et al., 2016). Furthermore, strain characteristics could interact with experimental manipulations. For example, mice that suffer from retinal degeneration cannot fully use visual information and exhibit impaired behavior in hippocampal-based navigation tasks that rely on visual cues (Nguyen et al., 2000). A second challenge is to ensure the stability of the tested phenotype under the null condition throughout the study period. In many cases, the effects of a given manipulation are studied over weeks to months (Prevot et al., 2019; Namdar et al., 2020; Rahn et al., 2021). Deterioration of control phenotypic behavior could confound potential manipulation-based intergroup differences. The present work showed that at the age of three months, HYB were larger than C57 and exhibited reduced anxiety-like behavior, but sensorimotor performance did not differ. Upon aging, C57 sensorimotor capabilities and learning performance deteriorated. In contrast, HYB sensorimotor capabilities did not deteriorate on aging, and the older HYB exhibited reduced anxiety-like behavior and improved long-term memory performance. Thus, HYB exhibit particularly stable performance during aging, which is potentially useful for long-term studies.
To provide a quantitative description of the HYB, we supplemented the array of standard assays with gait assessment and home cage behavior. The gait of all strains was inconsistent with abnormalities seen in motor disorders (Preisig et al., 2016; Rahn et al., 2021), yet every strain exhibited distinct gait patterns that were largely preserved upon aging. In addition, we observed apparent interstrain differences of daily conduct in the home cage, including distinct time division between sheltering, wheel running, and feeding. Together, the weight, daily conduct, and gait differences indicate that the spontaneous behavior of the (FVB/NJ x C57BL/6J)F1 hybrids is distinct from that of both progenitor strains.
There are some situations in which hybrid mice are not advantageous. First, compared with the C57 progenitors, HYB exhibited similar or improved performance in every assay used, with one exception: balance performance in the three-month-old age group. HYB impaired balance may result from weight differences between the strains since HYB weigh more, and weight affects rotarod performance (Mao et al., 2015). Although HYB balance performance improved upon aging, hybrids may be unsuitable for studies involving balance. Second, while the reduced anxiety exhibited by HYB is beneficial for studies involving cognitive and motor tasks, reduced anxiety-like behavior may be suboptimal for other purposes. For instance, the overall reduced anxiety-like behavior and the age-related changes make the HYB less suitable for long-term studies of anxiety. Third, some genes expressed by progenitors may be dominant; for instance, FVB contain a Disc1 mutation which has been shown to cause cognitive impairment even when heterozygous (Koike et al., 2006; Ritchie and Clapcote, 2013). Nevertheless, we found that the potentially deleterious effects of dominant genes on HYB are smaller than the adverse effects of inbreeding. Third, because HYB mice are more robust, genetic knock-out manipulations may be less likely to show a consistent phenotypic change (Sittig et al., 2016). Fourth, it is not straightforward to use hybrids if a homozygous allele is required or if knock-out of both copies of a given gene is required (although see Silva et al., 1997).
For studies that do not involve genetic manipulations, HYB can be used as-is, providing improved performance while maintaining minimal intersubject differences. Specifically, HYB can be implanted with high-density electrode arrays for neocortical and hippocampal recordings over multiple months. To combine transgenic control with enhanced behavior, we used two distinct strategies. One simple strategy involves using a homozygous parent of a C57-background, e.g., a driver line, yielding heterozygous HYB offspring that present one transgenic allele. For conditional expression of a reporter gene, HYB can be injected with a viral vector (as in Extended Data Figs. 6-1B, 6-3C). A second possibility is for one of the parents to be multi-transgenic, e.g., derived by crossing C57-based driver and reporter lines (as in Fig. 6; Extended Data Figs. 6-1A, 6-3A,B). Yet a third strategy is to use a C57-based driver line and backcross the reporter onto an FVB/NJ background. Because of the large body of transgenic mice developed on the C57 background, the number of possible combinations is very large, allowing to harness the behavioral potential of the hybrids for a wide range of studies.
We focused on first generation (F1) hybrids to minimize genetic differences between individuals. While further breeding is expected to increase intersubject genetic variability, F2 hybrids will necessarily express recessive traits not expressed by the F1 generation. For instance, 25% of the F2 hybrids are expected to be blind because of a functional copy of the Pde6b gene (Pittler and Baehr, 1991) maintained in every F1 hybrid. Thus, unless careful backcrossing is performed, F2 hybrids should not be used. Second, we employed HYB which are offspring of FVB females and C57 males. Full characterization is required before the use of reciprocal hybrids, offspring of a C57 female and an FVB male. However, because of the reduced litter sizes of female C57 compared with female FVB (Taketo et al., 1991; Silver, 1995), the usage of the reciprocal hybrids is not recommended. Finally, since interstrain differences are also sex-dependent (Brown and Wong, 2007), female HYB may exhibit distinct phenotypes.
To conclude, we showed that a single breeding step produces a hybrid mouse strain that exhibits improved behavioral performance while maintaining genetic control. The hybrid vigor and the enhanced behavioral capabilities may yield more trials in every experimental session, shorter learning periods, and higher accuracy in complex tasks. By increasing yield and accuracy, the usage of HYB may allow uncovering unexplored neuronal mechanisms. Enhanced behavior also has ethical advantages, since faster learning and higher yield per animal translate to fewer experimental animals or sessions required to answer a given question. The combination of genetic tools and the enhanced behavioral capabilities of the hybrid mice offers a unique opportunity for studying the neuronal basis of behavior.
Acknowledgments
Acknowledgments: We thank Ronni Cohen for constructive comments.
Footnotes
The authors declare no competing financial interests.
This work was supported by the United States-Israel Binational Science Foundation (BSF) Grant 2015577; the European Research Council Grant 679253; the Israel Science Foundation Grant 638/16; the Israel Science Foundation FIRST Program Grant 1871/17; the Rosetrees Trust Grant A1576; the Canadian Institutes of Health Research (CIHR), the International Development Research Centre (IDRC), the Israel Science Foundation (ISF), and the Azrieli Foundation Grant 2558/18; and the Zimin Institute.
This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.