A CRISPR toolbox for generating intersectional genetic mouse models for functional, molecular, and anatomical circuit mapping

The functional understanding of genetic interaction networks and cellular mechanisms governing health and disease requires the dissection, and multifaceted study, of discrete cell subtypes in developing and adult animal models. Recombinase-driven expression of transgenic effector alleles represents a significant and powerful approach to delineate cell populations for functional, molecular, and anatomical studies. In addition to single recombinase systems, the expression of two recombinases in distinct, but partially overlapping, populations allows for more defined target expression. Although the application of this method is becoming increasingly popular, its experimental implementation has been broadly restricted to manipulations of a limited set of common alleles that are often commercially produced at great expense, with costs and technical challenges associated with production of intersectional mouse lines hindering customized approaches to many researchers. Here, we present a simplified CRISPR toolkit for rapid, inexpensive, and facile intersectional allele production. Briefly, we produced 7 intersectional mouse lines using a dual recombinase system, one mouse line with a single recombinase system, and three embryonic stem (ES) cell lines that are designed to study the way functional, molecular, and anatomical features relate to each other in building circuits that underlie physiology and behavior. As a proof-of-principle, we applied three of these lines to different neuronal populations for anatomical mapping and functional in vivo investigation of respiratory control. We also generated a mouse line with a single recombinase-responsive allele that controls the expression of the calcium sensor Twitch-2B. This mouse line was applied globally to study the effects of follicle-stimulating hormone (FSH) and luteinizing hormone (LH) on calcium release in the ovarian follicle. The lines presented here are representative examples of outcomes possible with the successful application of our genetic toolkit for the facile development of diverse, modifiable animal models. This toolkit will allow labs to create single or dual recombinase effector lines easily for any cell population or subpopulation of interest when paired with the appropriate Cre and FLP recombinase mouse lines or viral vectors. We have made our tools and derivative intersectional mouse and ES cell lines openly available for non-commercial use through publicly curated repositories for plasmid DNA, ES cells, and transgenic mouse lines.


Conclusions:
The lines presented here are representative examples of outcomes possible with the successful application of our genetic toolkit for the facile development of diverse, modifiable animal models. This toolkit will allow labs to create single or dual recombinase effector lines easily for any cell population or subpopulation of interest when paired with the appropriate Cre and FLP recombinase mouse lines or viral vectors. We have made our tools and derivative intersectional mouse and ES cell lines openly available for non-commercial use through publicly curated repositories for plasmid DNA, ES cells, and transgenic mouse lines.
Keywords: Intersectional genetics, Gene targeting, CRISPR/Cas9, Cre, Flp, Dre, DREADDs, Fluorescent reporter Background Targeted expression of effector molecules, like fluorescent markers, calcium reporters, optogenetic actuators, or exogenous ligand-responsive receptors (DREADDs) [1], are increasingly applied in a variety of fields for greater precision, quantitative expression level control, and reduced side effects compared to previous methods for labeling and manipulation. For example, the P1 bacteriophage cyclization recombinase gene (Cre)/ locus of crossover in P1 (LoxP) system is used to achieve permanent cell-type restriction by using a promoter or enhancer to express the site-specific recombinase Cre [2]. In such examples, Cre transgenes are paired with a constitutively active, but conditional, allele where expression of an effector molecule, such as enhanced green fluorescent protein (eGFP), is interrupted by a LoxP-Stop-LoxP cassette that is recombined out by Cre, which enables expression in targeted, Cre-expressing cells. The use of various effector molecules in this paradigm enables fluorescent marking, neuronal perturbation, molecular affinity pull-downs, activity tracking, and other studies [1,.
In many fields, it is becoming increasingly clear that recombinase expression based on a single gene does not offer the resolution needed for a variety of developmental or targeting applications [29][30][31]. Indeed, application of intersectional genetics has led to new progress in various fields including neural circuits, cell-type lineage, and embryonic development [32][33][34]. Intersectional genetics adds needed resolution by employing a dual recombinase system using both Cre recombinase as well as a second recombinase, flippase (FLP), to activate a conditional effector allele only in cells where both recombinases have been expressed in the same cell (though not necessarily concurrently) [20,[35][36][37][38][39][40][41][42]. With new methodologies being developed to use Cre and FLP not only as traditional genetic markers, but also as activity and connectivity markers, unique combinatorial cell-type definitions become possible [43][44][45][46][47]. Thus, a resource consisting of multiple dual recombinase intersectional alleles that each express different effector molecules would add significant value and needed resolution to our ability to deconstruct neural circuits on multiple levels. This technology could also be applied to a multitude of other fields where intersectionally defined subpopulations of target cells exist and may play different roles in the measured outcomes.
Although mouse intersectional technology provides relatively benign access to otherwise inaccessible populations of cells, few laboratories have generated single transgenic or intersectional genetic mouse lines in house for several reasons. The complexity and size of the final targeting vectors puts them beyond present (cost-effective) commercial DNA synthesis capabilities, thus requiring some level of recombinant DNA cloning and precluding straightforward production of large, highfidelity ssDNA donors that can facilitate pronuclear clustered regularly interspaced short palindromic repeats (CRISPR)-mediated targeting [48]. While intersectional targeting vectors are available from the Addgene plasmid repository, they are finished vectors that require significant reverse engineering or modifications for use in a new context or approach. To our knowledge, there are no modular intersectional targeting vectors that are publicly available for facile and rapid production of new targeting alleles for the generation of intersectional mouse models. Furthermore, vector stability and other in vitro difficulties combined with the expense and time associated with target vector insertion and mouse line production limit the number of intersectional mouse lines available for public use. Thus, the production of intersectional genetic mouse lines has been largely limited to a small number of pioneering labs or resource-rich institutions such as the Howard Hughes Medical Institute and the Allen Institute for Brain Science [41,49].
To address these pitfalls and make the production of intersectional genetic mouse models more widely feasible, we aimed to produce a freely available resource toolbox consisting of several intersectional and single recombinase-responsive Rosa26 targeting vectors for rapid, facile, and cost-effective generation of complex mouse lines using CRISPR/ CRISPR-associated protein 9 (Cas9)-mediated homologous recombination in mouse embryonic stem cells and early oocytes. In oocytes, genomic insertions/deletions (in/dels) and short targeted insertions are readily produced whereas large construct targeting in oocytes has shown more limited but growing success. Additionally, requisite equipment and facilities are difficult to access and out of reach for many investigators. ES cell approaches are well established and widely available, allowing for rapid screening and identification, and, if using morula aggregation, require much simpler methodologies and equipment for mouse generation.
Toward this, we produced 7 intersectional mouse lines using a dual recombinase system; one mouse line with a single recombinase system, and three additional ES cell lines to study how functional, molecular, and anatomical features relate to each other in building the circuits that underlie physiology and behavior. As a proof of principle, we applied three of these lines to different neuronal populations for anatomical mapping and functional in vivo characterization in respiratory control. Next, we globally applied the single recombinaseresponsive line, which controls the expression of the calcium sensor Twitch-2B. Twitch-2B was expressed globally in the generated mouse line to study the effects of follicle-stimulating hormone (FSH) and luteinizing hormone (LH) on calcium release in the ovarian follicle. The publication and availability of this technology will allow for the seamless production of a highly diverse group of mouse lines that can be used to generate animal models of human disease, label specific cell populations for developmental or connectivity studies, or modulate cellular activity in established disease model lines, among other possibilities. All reagents and vectors used or generated in this study are now openly available for not-for-profit research.

Vector design and optimization
For each of the targeting vectors generated, the intersectional or Dre-responsive cassette was knocked into a well-established site in the Rosa26 locus ( Fig. 1) [50,51]. For positive ES cell clone selection, the targeting vector was simplified and improved in several ways compared to earlier approaches [52]. A simplified neomycin resistance cassette was integrated into the intersectional cassette before the second LoxP-flanked stop cassette to utilize the CAG promoter and polyadenylation (pA) sequences of the FLP recombinase target (FRT) flanked stop cassette, eliminating prior use of an additional 3phosphoglycerate kinase (PGK) promoter and Bovine Growth Hormone polyA (BGHpA) signal. To use CRISPR/ Cas9, we cloned an sgRNA that was close to the 5′/3′ homology junction into the px330 vector [53], which expresses both Cas9 and the subcloned sgRNA (px330_ Rosa26_sgRNA). A Woodchuck hepatitis virus (WHV) posttranscriptional regulatory element (WPRE) and BGHpA were added at the end of the expression cassette to enhance effector molecule expression and limit reliance on the disrupted Rosa locus for transcript termination [52]. The homology arms were significantly shortened to 1 kb to remove repetitive genomic sequences and stabilize the vector for prokaryotic propagation. Lastly, the terminal nonhomology Diphtheria Toxin A (DTA) chain negative selection cassette was removed, allowing for the complete targeting vector to be functionally tested via cell culture or in utero electroporation, which was not possible with the presence of the terminal DTA selection cassette without an additional subcloning step. Redesigned vectors and their Addgene ID numbers are outlined in Table 1 for public distribution.

Embryonic stem cell electroporation
We optimized traditional but widespread ES cell targeting by applying CRISPR/Cas9-mediated homologydirected repair (HDR) [54]. Our initial electroporation (EP) experiments using an earlier, more complex version of the Rosa26 intersectional targeting vector (containing the longer 4.2 kb 3′ homology arm, DTA, and no CRISPR/Cas9) showed targeting rates ranging from 6 to 14% (over the course of four electroporations we saw targeting rates of 6/48, 3/48, 7/48, and 3/48 clones, mouse strains listed in Table 2). To determine the effect of CRISPR/Cas9 on targeting efficiency, we coelectroporated the px330_Rosa26_sgRNA vector expressing Cas9 and a Rosa26-specific sgRNA with the optimized RR5 targeting vector at four different molar ratios (0:1, 0.5:1, 1:1, and 10:1), proportionally decreasing the amount of donor vector to accommodate the increasing px330_Rosa26_sgRNA vector per EP (18-20 μg total DNA per EP for 1 × 10 7 ES cells) (Fig. 1 C). Under the 0:1 px330 vector:donor vector condition, we saw a 6% targeting rate similar to those seen in electroporations without Cas9 and with a much longer 3′ homology arm and negative selection, suggesting that the shortening did not have a large impact on targeting efficiency under our EP conditions. Under the 0.5:1, 1:1, and 10:1 conditions, we observed 33%, 27%, and 58% targeting, respectively, suggesting that higher ratios of the Cas9 vector resulted in increased targeting efficiency (Fig. 1 D).
Since the 10:1 ratio had the highest targeting efficiency, despite having significantly less donor DNA, we used this ratio in subsequent EPs. Because the targeting efficiency was notably high, we next attempted a single electroporation with multiple vectors targeting the same locus, but containing equimolar amounts of five different cassettes, while keeping the overall px330_Rosa26_ sgRNA vector to total donor vector ratio at 10:1. The five vectors consisted of: (1) a Cre/FLP-responsive mutant methionyl-tRNA synthetase (RR6) for selective labeling of newly synthesized proteins; (2) a Cre/FLP-responsive modified G-protein-coupled receptor (RR7); (3) a Cre/ FLP-responsive bicistronic reporter with H2B-TagBFP and sfGFP separated by a self-cleaving peptide derived from porcine teschovirus-1 targeting 2A (p2a) element, and tdTomato expressed in the subtractive population (cells that express FLP but not Cre) (RR9); (4) a Dreresponsive tricistronic reporter with histone 2B-tagged with blue fluorescent protein (H2B-TagBFP), superfolder GFP (sfGFP), and synaptophysin-tdTomato separated by p2a elements (RR11); and (5) a Dre-responsive bicistronic reporter with H2B-TagBFP and sfGFP separated by a p2a element (RR10). We saw a 52% targeting efficiency and successful targeting of all five cassettes at varying efficiencies (Table 3). Due to the lower targeting efficiency of RR7 (2%), we attempted another 10:1 electroporation using the RR7 donor vector alone and obtained a 63% targeting efficiency (genotyping results shown in Fig. 1). Thus, our results show that we are able to target and recover as many as five Rosa alleles in a single EP, significantly increasing efficiency and reducing costs toward intersectional mouse generation.

Oocyte targeting
Given the high efficiency of targeting, we also attempted to create a founder line (RR8) through direct oocyte Targeting knock-in was determined using PCR primers that spanned from outside the Rosa26 homology arms to either the CAG promoter (5′ end) or WPRE (3′ end). Amplification of a band indicates targeting. Shown are results from a targeting event with over 60% targeting efficiency (RR7). C Four conditions with different Cas9:targeting vector ratios were used in our initial study: a 0:1 with no Cas9, 0.5:1, 1:1, and 10:1. D Targeting efficiency results from the different ratios. The 10:1 Cas9:targeting vector ratio showed a 5-10-fold increase over an electroporation with no Cas9. Shown in red are results from previous electroporations using the traditional Rosa26 targeting vector with more commonly used longer homology arms injection of the optimized Rosa26 targeting vector with a Cre-responsive calcium indicator, Twitch [55] (total cassette size without homology arms, 6.5 kb), using pronuclear injection of Cas9 protein (30 ng/μl), sgRNA (20 ng/μl), and double-stranded DNA plasmid (2 ng/μl). A total of 1266 embryos were injected yielding 129 mice, of which 7 genotyped positive for the calcium indicator but only one targeted successfully (< 1% targeting efficiency of mice born), suggesting that direct oocyte injection is inefficient with this system under the specified parameters.

Off-target analysis
Genomic sequences that are similar to the sgRNA used for targeted double-stranded breaks may cause unintended gene mutations or editing at off-target sites. To account for this possibility, we predicted the off-target sites for each sgRNA using the crispr.mit.edu tool and selected the top 5 sites for follow-up. We PCR amplified and sequenced these loci from three correctly targeted ES cell clones and did not detect any genetic changes (Fig. 2). While we cannot rule out off-target effects in other loci, these data suggest that off-target effects are not prevalent in this setting.

Mouse line derivation
The goal of our efforts was to produce freely available genetic tools that delineate and access distinct populations for multifaceted circuit mapping or functional characterization studies. With this toolkit in hand, we collectively produced 11 intersectional or Dre-responsive alleles available as either mouse lines (8) or ES cells (3) ( Table 2). To perturb neuron function, we produced three lines that express metabotropic DREADD receptors, modified human M4 muscarinic receptor (hM4D) (RR1), modified human M3 muscarinic receptor (hM3D) (RR2), and RASSL (Rs)-EGFP (RR7), and a fourth line expressing the ionotropic pharmacologically selective actuator module (PSAM) receptor (RR4) [24,56,57]. Two lines enable molecular characterization: one expresses the EGFP-L10A fusion protein for ribosomal affinity purification and capture of translating mRNAs (RR3) [58], while the second line expresses methionyl-tRNA synthetase (Met-RS) that incorporates an artificial amino acid into nascent peptides (RR6) [13]. One line enables ratiometric calcium imaging via expression of Twitch2B (RR8). Our neuro-anatomical mapping lines described below also incorporate a tagged histone (H2B), offering the possibility of chromatin affinity isolations. We Table 1 Publicly available modular targeting vectors for rapid generation of recombinase-responsive ES cells and mouse lines created four alleles of differing Cre, FLP, and Dre recombinase-responsive configurations that fluorescently mark distinct cellular compartments for unambiguous cell counts, morphological characterizations, and projection mapping (RR5, RR9, RR10, RR11). The derivation of lines RR5, RR6, and RR7 demonstrated that ES pluripotency was maintained and that germline transmission was not diminished in our Cas9-mediated single and multiplexed electroporations after correctly targeted clones were selected and injected into blastocysts for chimera generation.

Select mouse line characterization
Upon dual recombinase expression, the RR5 tricistronic multi-color reporter allele expresses three spectrally separated and modified fluorescent proteins to highlight the nucleus in blue, fill the neuron in green, and emphasize presynaptic contacts in red. At the site of dual recombinase expression, targeted cells are brightly labeled with TagBFP, sfGFP, and tdTomato where TagBFP fluorescence is constrained to the nucleus by an H2B fusion, sfGFP is unmodified so that fluorescence is seen throughout the cell body, and tdTomato is fused to Table 2 Rosa26 knock-in alleles generated with or without CRISPR/Cas9 methods. Shown are targeted Rosa26 alleles, function, length of cassette, current status, and expected MGI name (to be determined) in grey rows synaptophysin [59] so that fluorescence is excluded from nuclear areas and primarily seen in projection areas. Individual nuclei can be resolved using blue fluorescence and co-localization with sfGFP-labeled somata.
We used RR5 to evaluate the functional activity and specificity in vivo of our CRISPR/Cas9 approach to generating mouse lines in three distinct contexts; (1) germline recombinase expression, (2) viral recombinase expression, and (3) combinatorial retrograde viral and germline recombinase expression to target single gene defined neurons by a specific projection field. First, to demonstrate genetically restricted, germline expression of recombinases in the intersectional RR5 line, we bred the RR5 line to a double dopamine beta-hydroxylase (DBH) p2aFLPo [54]; Bactin-Cre recombinase driver to express the tricistronic fluorescent cassette in DBH-positive noradrenergic (NA) neurons in the brainstem (  Second, we verified viral expression of Cre and FLP recombinase where adult RR5 mice were stereotaxically injected with equal titer amounts of adeno-associated viruses (AAVs) expressing Cre and FLP ( Clear expression of all three fluorescent proteins were resolved at the site of injection to the appropriate cellular locations. Third, defining properties of the subpopulations can be extended to include any combination of gene expression or anatomical location and projection target when injection of retrograde viral vectors is applied. To utilize this option, we genetically marked a neuronal subtype defined by the partial overlap of a genetic selector and projection target by injecting canine adenovirus 2 (CAV2)-Cre virus [62,63] into the basolateral amygdala of DBH p2aFLPo ; RR5 mice that express only FLP in all noradrenergic neurons. CAV2-Cre virus efficiently transduces axon terminals, thus genetically marking presynaptic neurons. In this context, only noradrenergic neurons expressing DBH that project to the amygdala express both Cre and FLP and the resulting tricistronic fluorescent cassette (Fig. 5A). We only observed recombination in the brainstem noradrenergic nuclei and in the locus coeruleus, primarily on the ipsilateral side ( Fig. 5B-E) with some sparse marking on the contralateral side ( Fig. 5F-I), in agreement with a prior study [64]. Overlapping red and green puncta (with no blue) could be seen in the injected amygdala, suggesting that the marked neurons project to the amygdala, as expected given the nature of the CAV2-Cre virus ( Fig. 5J-M). We also observed collateral projections to several additional areas in the mid and forebrain, including the dorsal raphé ( Fig. 5N-Q), reticulotegmental pontine Before in vivo characterization of the DREADD systems coded by RR1 and RR2, we first examined the ability of the RR1 and RR2 lines to modulate dopamine beta-hydroxylase (DBH)-defined noradrenergic neuron activity at the cellular level. We conducted whole-cell recordings from postnatal day (P)30-P60 locus coeruleus neurons expressing one of 3 DREADD receptors: hM3D (RR2), hM4D (RR1), or Di (previously published). Each cassette exists with FRT and LoxP bound stop cassettes, so each line was initially bred with Bactin_FLPe to remove the FLP-dependent stop cassette [65]. The RC::FrePe was also bred with Bactin_FLPe to remove the FLP-dependent stop cassette and then bred to  the cre-only dependent DREADD lines [66]. This compound allowed for fluorescent labeling of neurons where Cre recombinase, and thus DREADD receptors, are present. These compounds resulted in mice wherein GFP and a DREADD were expressed in TgDBH-Cre-defined cells for acute brain slice visualization (See Additional File 2: Fig. S2 for further genetic details). RC::PDi [52] was included as it represents a similarly constructed hM4D intersectional allele, but lacks the modifications and optimizations made in our targeting system (Fig. 1A). . These data confirmed functional expression of each DREADD system. All raw data from these experiments are available on Zenodo (10.5281/zenodo.5794197). We next performed whole-body plethysmography to assess changes in respiratory function caused by CNO-DREADD-mediated perturbation of targeted neurons in vivo. Significant changes in any of the calculated respiratory variables demonstrate a measurable effect of DBH-defined noradrenergic neuronal activity on respiratory control and/ or ventilatory response to CO 2 . If change was seen in room air, then basal respiratory control was affected by the change in activity of these neurons. If changes were seen in 5% CO 2 , then changes in the activity of these neurons affected the hypercapnic ventilatory response. DBH-Cre-defined NA neurons were inhibited by crossing mouse line RR1P to Tg(Dbhcre)KH212Gsat (TgDBH-Cre) mice. We observed no changes under room air conditions (Fig. 8A). However, under hypercapnic conditions, we saw a reduction in ventilation (Fig. 8B), with decreases in V T (tidal . Measured values include respiratory frequency (B), tidal volume (C), ventilation (D), oxygen consumption (E), and ventilation normalized to oxygen consumption (G). No difference was seen under room air conditions but DBH-Cre; RR1P animals showed a deficit in volume, ventilation, and ventilation normalized to oxygen consumption under hypercapnic conditions. G No difference in temperature was seen between DBH-Cre; RR1P animals and sibling controls. *p < 0.05; **p < 0.01; ***p < 0.001 volume) (− 26.6%, p = 0.0021, Fig. 8D),V E (ventilation) (− 26.8%, p = 0.026, Fig. 8E), andV E /V O2 (ventilatory equivalents for oxygen) (− 28.7%, p < 0.0001, Fig. 8G). No significant differences were observed in body temperature between experimental animals and sibling controls (Fig. 8H). Notably, sibling controls showed no difference in respiratory parameters pre-and post-CNO administration. All raw data from these experiments are available on Zenodo (10.5281/zenodo.5781763).
We then tested the applicability of the RR2P line in vivo. After CNO-DREADD-mediated stimulation of DBH-Cre-defined neurons, we saw significant increases under room air ventilation (Fig. 9 A) in f R (Respiratory frequency) (+ 28.1%, p = 0.0076, Fig. 9 C), V T (+ 43.3%, p = 0.00077, Fig. 9D),V E (+ 83.6%, p < 0.0001, Fig. 9E), andV O2 (Oxygen consumption) (+ 50.9%, p < 0.0001, Fig. 9F). The increase inV E was proportional to the increase inV O2 , however, soV E /V O2 did not change (Fig. 9G). Under hypercapnic conditions, as compared to pre-CNO values,V O2 was increased (+ 37.2%, p = 0.0015, Fig. 9F), resulting in a reducedV E /V O2 (− 24.2%, p = 0.04, Fig. 9G). We also observed a deficit in body temperature in DBH-Cre; RR2P animals 30 min after the end of the assay at room temperature (~1.5 h after CNO injection) as compared to sibling controls (Fig. 9H), but not immediately at exit from the respiratory chamber, which is kept at a thermoneutral temperature, 30-32°C. Sibling controls showed no difference in respiratory parameters pre-or post-CNO administration. The respiratory phenotypes seen here clearly demonstrate an acute and cell-autonomous involvement of the whole NA system in CO 2 chemosensitivity and baseline metabolism while avoiding the confounds of developmental compensatory events or off-target effects. All raw data from these experiments are available on Zenodo (10.5281/zenodo.5781775).
Finally, we characterized a fourth line, produced by oocyte CRISPR-mediated recombination to express the Twitch2B ratiometric calcium indicator upon Cre recombinase expression. Calcium signaling was recorded in outer granulosa cells surrounding an oocyte taken from a germline recombined (RR8; Bact_Cre) mouse line (Fig. 10A). Upon application of 10 μM ionomycin in the presence of 10 mM CaCl 2 , the yellow fluorescent protein (YFP)/ cyan fluorescent protein (CFP) ratio increased more than 3-fold (Fig. 10B). The data demonstrate how the targeting system here can be used to build conditional mouse lines that are useful beyond the nervous system and throughout the body [67]. All raw data from these experiments are available on Zenodo (10.5281/zenodo.5781661).
Collectively, these data highlight how our single recombinase and intersectional lines can be used to functionally assess dispersed and difficult to access populations, setting the stage for further characterization and dissection of a variety of systems in a given measured outcome through intersectional subdivision by genetic or viral methods.

Vector design and optimization
Our tricistronic mouse line vector (RR5) was a complex vector to design. Thus, the discussion of its construction applies to the other simpler vectors and best describes choices made to optimize the vectors. Our initial goal was to create an intersectional mouse strain expressing three fluorescent components to highlight the nucleus, whole neuron, and presynaptic contacts, enabling facile study of cell counts, cellular morphology, and axonal and projection targets (RR5 consisting of H2B-TagBFP, sfGFP, and synaptophysin-tdTomato separated by p2a elements) without the need for signal amplification using antibodies. However, likely due to the high number of repetitive elements in the Rosa26 3′ homology arm and the high complexity and repetitive components of the intersectional cassette (i.e., two stop cassettes consisting of multiple SV40 polyadenylation sequences, tdTomato protein dimer, double loxP and FRT sites, and two p2a self-cleaving peptides separating the elements), we had no success in stably subcloning the intersectional cassette into a low copy vector (p15a ori) containing the original Rosa26 homology arms (1.1 kb 5′ homology arm and 4.2 kb 3′ homology arm) and DTA negative selection gene [52]. Thus, we sought to make several improvements on the prior approaches used to build intersectional genetic mice. First, we sought to stabilize the vector by shortening the homology arms and eliminating repetitive genomic elements found at the Rosa26 locus. Second, we removed the DTA negative selection terminal non-homology to eliminate the need to validate intersectional cassettes separately before placing the cassette in the targeting vector. Third, we reduced the size of the selection cassette. Prior intersectional constructs carried an additional PGK promoter and BGHpA tail for neomycin expression. We eliminated both promoter and BGHpA as the functions of these elements can easily be served by the CAG promoter and SV40pA in the stop cassette. With smaller homology arms, it was possible that our targeting efficiency would fall. Therefore, we introduced a CRISPR/Cas9 targeting strategy to enhance homologous recombination in either ES cells or oocytes. Surprisingly, in ES cells, the shorter homology arms performed equally well as the full-length homology arms (see "Discussion" below). The resulting strategy offered a simplified, smaller, and highly stable vector that remains efficient at targeting the Rosa26 locus. Through this optimization process, we built a set of baseline vectors  79%N 2 ). B Representative breathing traces from a DBH-Cre; RR2P animal before and after CNO administration under hypercapnic conditions (5%CO 2 /21%O 2 /74%N 2 ). C-G Quantification of respiratory and metabolic parameters under room air and hypercapnic conditions in DBH-Cre; RR2P animals (n = 16) and sibling controls (n = 14). Measured values include respiratory frequency (C), tidal volume (D), ventilation (E), oxygen consumption (F), and ventilation normalized to oxygen consumption (G). After CNO administration, DBH-Cre; RR2P animals showed increased respiratory frequency, volume, ventilation, and oxygen consumption under room air conditions. Under hypercapnic conditions, DBH-Cre; RR2P animals showed increased oxygen consumption resulting in a reduced ventilation to oxygen consumption value. I DBH-Cre; RR2P animals showed a significant deficit in temperature 30 min after the end of the assay. *p < 0.05; **p < 0.01; ***p < 0.001 that can be used to build any intersectional genetic targeting allele rapidly and inexpensively by simply cloning in a cDNA of interest into one or both of the multiple cloning sites for both intersectional and subtractive expression.

Embryonic stem cell electroporation
While CRISPR strategies are widespread for knockout mutations and small deletions and insertions in mouse zygotes, consistent knock-in of large targeting cassettes (> 4 kb) still remains a challenge with fewer but increasing successes [48,68,69]. To address this limitation, our new multiplexed methodology enabled us to: (1) simplify our base targeting vector to less than 5 kb in total length; (2) see successful targeting of added cassettes up to 11 kb; (3) increase the rate of targeting by 5-10 folds over previous Rosa26 targeting attempts (that used a significantly longer 3′ homology arm) under our ES cell strain and conditions; and (4) further reduce cost by coelectroporating five different targeting vectors in a single ES cell electroporation that were easily resolved to produce new mouse lines. While targeting rates for the Rosa26 locus can vary greatly in the literature, we compare only to our own experiments, as significant variability in targeting efficiency from lab to lab or facility to facility can arise from the several factors in ES cell electroporation ranging from electroporation conditions (e.g., buffer ionic strength, adjuvants, field strength, field duration, field shape, square wave vs exponential decay) to culture conditions (ES cell strain and genetic background, feeder type, media composition, inclusion of LIF). Notably, when we co-electroporated a CRISPR/ Cas9 plasmid to facilitate targeting, there was a proportional increase in targeting efficiency despite reduced amounts of targeting vector (total DNA in the EP was capped at 18-20 μg per EP). At the highest efficiencies, we could further subdivide the small amount of targeting vector across multiple plasmids, which allowed us to introduce and recover as many as five vectors in a single EP. This multiplexed approach significantly reduced the costs of mouse production via ES cells. With these goals met, multiplexed ES cell targeting enables the rapid and high-throughput production of intersectional alleles that can be readily distributed throughout the mouse research community for further studies.

Oocyte targeting
Oocyte injection was only attempted with one line, RR8. Less than 1% of the embryos contained a properly targeted vector, which in our experience suggests that direct oocyte injection is not efficient. Future optimization and the development of new technologies may increase the efficiency of targeting for more rapid generation of transgenic mice; for example, a study showed 10-20% targeting of 8-11 kb cassettes to the Rosa26 locus in oocytes using longer homology arms, while another showed successful targeting of a 5.5-kb cassette into the Rosa26 locus in rats with co-injection of two ssDNA oligos [69,70]. As noted above, there can be significant variability in targeting across labs based upon a variety of factors including background strain (vs strain from which the homology arms are derived), the use of Cas9 mRNA or protein, and site of injection (pro-nucleus vs cytoplasm, single vs 2-cell stage). Thus, our outcomes here may reflect our core facility capabilities as much as inherent limitations in consistent Cas9-mediated oocyte targeting of large constructs. The successful outcomes here and in a prior publication demonstrate the use of the system to develop lines capable of functional imaging [67].

Select mouse line characterization
Intersectional expression of fluorescent proteins (FPs) has been used to great effect for anatomical characterizations of targeted neuron populations [66,71]. Characterizations aimed at counting or highlighting single cells or cell populations made with a single FP can be made difficult, however, in cases where cell, axon, or dendrite density (neuropil) is high. Usually, one or two FPs are deployed to either fill the cell or highlight specific features, which necessitates the use of multiple reporter alleles to clearly observe cell number, morphology, and projections. To ameliorate these issues, we sought to create a mouse line (RR5) that highlights three specific features in an intersectionally defined neuron. Overall, this novel three-color mouse line enables intersectional and simultaneous labeling of three subcellular compartments. The polycistronic fluorescent labeling cassette gives unambiguous cell counts and clear visualization of cell morphology, axonal projections, and synaptic contacts. The strength of expression of all three colors is likely due to the use of p2a elements, rather than internal ribosome entry site (IRES) elements, which cause expression level reduction. Reduced expression from an IRES may stem from cryptic splice elements in the IRES that could interact with the intron in the CAG promoter or splice acceptors in some stop cassettes [72]. Our use of p2a appears to bypass those issues as all three fluorescent proteins could be readily visualized without additional enhancement as compared to other approaches [38].
We also show that combinations of germline, viral, and retrograde viral expression of recombinases in the RR5 line can all be used to clearly map the shared innervation of multiple regions by a small subset of genetic, projection, and/ or intersectionally defined neurons. In contrast to prior intersectional studies using two genetic delimiters to define a subpopulation with a single fluorescent protein (with or without an inert retrograde marker injected into separate sites in the brain), our approach enables facile collateral mapping across the whole nervous system with higher resolution that is growing in popularity [36,64,71]. These results demonstrate that despite the high complexity and repetitive nature of the RR5 cassette, we could efficiently target this allele, the allele remains stable through germline transmission, and it is responsive to Cre and FLP recombinase expression. The restriction of expression to anatomically defined NA neurons as a whole or as a projection subset further establishes the recombinase specificity of our targeting schematic. This approach should also translate equally well to further anatomical, functional, and molecular characterizations of NA and other subtypes defined by their projection patterns with the use of the additional lines reported here.
The outcomes of either inhibiting (RR1P; TgDBH_Cre) or stimulating (RR2P; TgDBH_Cre) NA neurons as defined by TgDBH_Cre strongly support the functionality of these two lines in an in vivo setting. NA neurons are strongly implicated in control of respiratory homeostasis and chemosensitivity [8,[73][74][75][76][77][78][79]. Previous studies focused on demonstrating the role of NA neurons in breathing bear methodological caveats such as lack of resolution due to overly broad lesions or injections [80], developmental and non-cell-autonomous compensation after gene mutations, restraint or anesthesia in vivo, and they typically lack concurrent metabolic measurements ( V E /V O2 ) [81]; thus, direct comparisons are difficult. Nonetheless, previous work has been reiterated in many ways with our results here in that focal pharmacological NA lesions decrease respiratory frequency and hypercapnic ventilatory response [75]. Here we show in conscious, unrestrained mice that inhibition of TgDBH_ Cre-defined noradrenergic neurons led to significantly reduced tidal volume, ventilation, andV E /V O2 following hypercapnic exposure. On the other hand, activation of TgDBH_Cre-defined noradrenergic neurons led to significantly increased respiratory frequency, tidal volume, ventilation,V O2 , andV E /V O2 in the room air breathing condition. These changes likely stem from a combined increase in central and peripheral (sympathetic) outflow increasing overall metabolic rate (V O2 ). Although breathing andV O2 are both increased, a small but significant mismatch occurs in hypercapnia suggesting that global excitation may destabilize hypercapnic reflexes as well. Notably, after removal from the thermoneutral recording chamber, the body temperature dropped significantly. This drop is consistent with prior studies that show central NA outflow to the hypothalamus negatively regulates body temperature [82,83]. The magnitude of changes seen for both stimulation and inhibition of TgDBH_Cre defined NA neurons in various breathing parameters, particularly for hypercapnia, is in line with expectations of a highly redundant chemosensory system in the respiratory network along with the fact that the NA system is neuromodulatory in nature [37,52,[84][85][86][87]. Altogether, the numerous breathing, metabolic, and temperature changes observed support the functionality of these lines for in vivo chemogenetic manipulations. These data represent modulation of the noradrenergic system that, methodologically, significantly deviates from previously published studies interrogating this system and provides an alternative to those investigating control of breathing that reduces the influence of unintended effects.
In conclusion, the methods presented here represent a significant improvement in feasibility and affordability for labs to generate their own intersectional mouse models for application in a variety of fields. Our vector design and application allow for large cassettes to be targeted to the Rosa26 locus with high throughput, increased efficiency, germline transmission, and preserved pluripotency. These protocols and techniques can be used to investigate most, if not all, neural circuits in the context of development and/or disease. Implementing these techniques also serves as an alternative to costprohibitive commercially available vectors and mouse lines. Together, we present a facile, cost-effective method for producing gene-specific intersectional mouse lines.

Conclusions
During development, many distinct cellular subtypes arise and intercalate to create the complex cells, organs, and organ systems that constitute our behavior and physiology. Key to understanding this organization and how it may be perturbed in diseases is the ability to identify and access discrete cell subtypes in the developing and adult mouse for multifaceted studies. It is now clear that even within narrowly defined cell types once thought as homogeneous, significant diversity is found at multiple levels including genetic and molecular signatures, activity patterns, synaptic connectivity, and projection patterns. Given this complexity, it is often difficult to define and access specific cellular populations for study, particularly during embryogenesis, where in utero development makes access problematic. In neuroscience, the power of studying combinatorially defined neuronal populations increases significantly when multiple distinct neuronal features, such as birthdate, collateral projection targets, molecular profiles, or a functional requirement in a given physiological process or behavior can be examined in parallel to reveal deeper insights into the mechanisms involved in the establishment and maintenance of neural circuit organization. However, the ability to carry out such studies is restricted by the limited number of publicly available intersectional mouse lines as well as the inherent difficulty in developing new intersectional alleles.
We present a novel CRISPR/Cas9-mediated system for consistently targeting large cassettes to the Rosa26 locus for high-throughput and parallel generation of intersectional and conditional mouse lines. These studies demonstrate that CRISPR/Cas9 can be used to increase targeting efficiency at the Rosa26 locus in mouse embryonic stem cells with shortened 1-1.2-kb homology arms while preserving pluripotency and the ability to transmit alleles via the germline. This technology allows for the opportunity to generate mouse lines that express a wide range of effector molecules in various cell types, enabling the developmental, anatomical, molecular, and functional characterization of cellular organization in behavior and physiology. To democratize the generation of intersectional and conditional targeted ES cells and mouse lines, we have generated a publicly available vector toolkit for single-step cloning into Rosa26-targeting vectors for double and single recombinase-responsive cassettes (Table 1). We have also generated several mouse lines available, without restriction, to the not-forprofit research community enabling high specificity intersectional access of neuronal and other populations for cross correlative mapping approaches using functional neuronal perturbation, molecular profiling, and anatomical characterizations for multifaceted studies (Table 2). Overall, these resources and the modular nature of the intersectional approach make rapid and low-cost production of large numbers of intersectional alleles possible that can be efficiently distributed as ES cells or mouse lines throughout the mouse research community. All vectors (Addgene; 97007-97012, 99142), mouse lines (MMRRC; 043513-043519), and ES cells (BCM ES Cell Core Facility) have been made publicly and unconditionally available to the academic research community.

Construction of targeting and Cas9 vectors
The simplified Rosa26 targeting vector was derived using standard cloning procedures from a targeting vector used in previous studies [52] consisting of a p15a origin of replication, kanamycin selection cassette, DTA negative selection cassette, and Rosa26 homology (1081 bp 5′ homology arm and 4342 bp 3′ homology arm). The total vector size was 9084 bps. The DTA was removed and the 5′ homology arm was shortened to 1025 bp fragment (deleting the sgRNA target sequence (chr6:113,076,062-113,077,087 GRCm38/mm10)) while the 3′ homology arm was shortened to 1231 bps (chr6:113,074,801-113,076,032 GRCm38/mm10), maintaining the insertion site within the first intron and reducing the total vector length to 4694 bps and eliminating a number of repetitive regions that often resulted in vector instability when constructing intersectional targeting vectors.
For the Cas9 expressing vector, we selected an sgRNA target sequence that was close to the 5′ and 3′ junction of the Rosa26 gene locus (ACTGGAGTTGCAGATC ACGA with PAM motif GGG; chr6:113,076,040-113,076,05 GRCm38/mm10). The selected sgRNA was cloned into the BbsI sites of the px330 vector that expresses Cas9 (final vector is called px330_Rosa26_ sgRNA).
Knock-in cassettes were assembled in either (1) a previously used intersectional template plasmid (RR1-4, RR6-7 OL963770-OL963773; OL963775-OL963776); (2) a newly constructed intersectional template plasmid with neomycin promoter and pA elements removed (RR5 and RR9, OL963774; OL963778); (3) a Cre-only derivative of the newly constructed intersectional template plasmid with neomycin promoter and pA elements removed (RR8, OL963777); or (4) a constructed Dre-responsive template plasmid (RR10 and RR11 OL963779; OL963780). The intersectional plasmid consisted of a ubiquitous CAG promoter sequence; an FRT-flanked stop cassette consisting of a PGK-neo sequence (for positive selection of targeted ES cells) and three SV40 pA sequences; a LoxP-flanked stop cassette containing mCherry and a PBS302 stop cassette; and a cloning site for insertion followed by a WPRE and BGHpA sequence. The new intersectional template consisted of a ubiquitous CAG promoter sequence; an FRT-flanked stop cassette consisting of a neomycin sequence and His3_SV40 pA; a LoxP-flanked subtractive cloning site stop cassette consisting of a His3_SV40 pA; and a cloning site for insertion followed by WPRE and BGHpA sequences. The Dre-responsive template consisted of a Rox-flanked stop cassette consisting of a neomycin sequence and PBS302 stop cassette, and a cloning site followed by a WPRE and BGHpA sequence.
For assembly of complete targeting vectors, the cDNAs of interest were PCR amplified and cloned into the corresponding vector. The intersectional or Dre-responsive cassettes were then cut out with PacI or PacI/AscI and cloned into the shortened Rosa26 targeting vector.
To facilitate single-step cloning of complete targeting vectors, we generated a vector toolkit (Table 1) by constructing and moving empty CAG cassettes in the Rosa26 targeting vector. We generated four template vectors with neomycin resistance: (1) a Cre/FLP-responsive targeting vector; (2) Cre only; (3) FLP only; and (4) Dre-responsive targeting vector. The intersectional template vector consists of the CAG promoter, an FRTflanked neomycin resistance gene and stop cassette, a subtractive cloning site (EcoRV) and stop cassette flanked by LoxP sites, an intersectional cloning site (PmeI), and WPRE and BHGpA elements. Each single template vector consists of the CAG promoter, stop cassette(s) and neomycin resistance cassette flanked by recombinase sites, cloning site (SwaI), and WPRE and bgh polyA elements. Additionally, for possible direct oocyte injections, we also generated a Cre/FLP-responsive targeting vector that does not contain a neomycin resistance gene, with a subtractive cloning site (EcoRV) and intersectional cloning site (SwaI). All cloning sites are blunt restriction enzyme sites for cloning facilitation. Sequences of the template vectors and px330 vector and plasmids are available through Addgene (ID #s 97007-97012, 99142).

Generation of knock-in mice
Embryonic stem (ES) cells (AB2.2) were electroporated with 15-20 μg of varying ratios of the px330_Rosa26_sgRNA vector to the targeting vector. Neomycin selected clones were screened for homologous recombination using PCR genotyping. Targeted clones were identified using PCR genotyping for 5′ and 3′ targeting from outside the homology arm and were considered to be successful targeting events if positive for both 5′ and 3′ genotyping. We used pairs 5′CGCCTAAAGAAGAGGCTGTG (Rosa26-F) and 5′GGGCGTACTTGGCATATGAT (CAG-R), producing a 1450 bp band; and 5′AATCAACCTCTGGATTA-CAAAATTT (WPRE-F) and 5′TGGCTCCTCTGTCCA-CAGTT (Rosa26-R), producing a 2472-bp band. Select targeted clones were microinjected into C57B1/J6 blastocysts and chimeric males were bred to wildtype C57B1/J6 females to achieve germline transmission.

Pronuclear injection
We constructed a targeting vector by cloning a calcium indicator, Twitch-2B [55], into the Cre-responsive Rosa26 targeting vector. The Baylor College of Medicine Embryonic Stem (ES) Cell Core generated the Rosa26specific sgRNA and Cas9 protein for injection. Pronuclear injections were performed by the Baylor College of Medicine Genetically Engineered Mouse (GEM) Core using the following parameters: 30 ng/μl Cas9 protein, 20 ng/μl sgRNA, and 2 ng/μl dsDNA plasmid. Potential founders were screened for targeting as described above.
Established mouse lines reside at the Mutant Mouse Regional Resource Center (MMRRC) and will be unconditionally available to the academic research community (MMRRC 043513-043519). Targeted ES cells have been archived and are available upon request under the condition that any new mouse lines are deposited in a public repository.
All animal experiments were performed with the approval of IACUC.

Off-target analysis
Off-target sequences were identified using the Optimized CRISPR Design tool (crispr.mit.edu). The top 5 sequences for each sgRNA were amplified from selected targeted ES cells and sequenced to determine if any mutations or changes occurred in off-target sites.

Histology
For RR5 expression, animals were sacrificed and transcardially perfused with 0.1 M phosphate-buffered saline (PBS) then with 4% paraformaldehyde (PFA) in PBS. Brains were dissected out and drop fixed for 2 h in 4% PFA before a PBS rinse and equilibration in 20% sucrose in PBS. Brains were sectioned into 30-40-μm coronal sections and mounted on slides. Images were collected on a Zeiss confocal LSM780 microscope or Leica TCS SPE confocal microscope.

Slice preparation
Slice preparation from the mouse LC follows an N-Methyl-D-glucamine (NMDG) slicing protocol [88,89]. Briefly, animals were deeply anesthetized using 3% isoflurane. After decapitation, the brain was removed and placed into cold (0-4°C) oxygenated NMDG solution containing 93 mM NMDG, 93 mM HCl, 2.5 mM KCl, 1.2 mM NaH 2 PO 4 , 30 mM NaHCO 3 , 20 mM HEPES, 25 mM glucose, 5 mM sodium ascorbate, 2 mM Thiourea, 3 mM sodium pyruvate, 10 mM MgSO 4 , and 0.5 mM CaCl 2 , pH 7.35 (all from SIGMA-ALDRICH). Horizontal slices were prepared using a vibratome (200 μm thick) using zirconia blades. The brain slices were kept at 37.0 ± 0.5°C in oxygenated NMDG solution for 10 min. They were then transferred to an artificial cerebrospinal fluid (ACSF) containing 125 mM NaCl, 2.5 mM KCl, 1.25 mM NaH 2 PO 4 , 25 mM NaHCO 3 , 1 mM MgCl 2 , 25 mM glucose, and 2 mM CaCl 2 (pH 7.4) for at least 1 h prior to the beginning of recordings. During the recording sessions, the slices were submerged in a commercially available chamber (Luig Neumann, Order No. 200-100 500 0150-M) and were stabilized with a fine nylon net attached to a custom-designed platinum ring. This recording chamber was continuously perfused with oxygenated physiological solution throughout the recording session.

Recordings
Whole-cell recordings were performed as described previously [89][90][91]. Briefly, patch pipettes (2-7 MΩ) were filled with an internal solution containing 120 mM potassium gluconate, 10 mM HEPES, 4 mM KCl, 4 mM MgATP, 0.3 mM Na 3 GTP, 10 mM sodium phosphocreatine, and 0.5% biocytin (pH 7.25). Whole-cell recordings from up to 8 LC neurons were performed using two Quadro EPC 10 amplifiers (HEKA Electronic, Germany). PatchMaster (HEKA) and custom-written Matlab-based programs (Mathworks) were used to operate the recording system and perform online and offline data analysis. In current-clamp recordings, neurons were first current clamped at~− 40 pA to prevent spontaneous firing. Prior to investigating the effect of drugs, we calculated spike thresholds and recorded firing patterns in response to sustained depolarizing currents by injecting increasing current steps (+10pA). Continuous recordings were obtained from LC neurons current clamped at − 40 to 0 pA during drug wash-on experiments. We also calculated other intrinsic electrophysiological parameters, such as the input resistance, membrane time constant, spike amplitude, and after-hyperpolarization (AHP) [89][90][91].

Plethysmography
Plethysmography on conscious, unrestrained mice was carried out as described on 6-12-week-old adult animals with at least 12 animals in both the experimental and control group for a given experiment [52,87]. Mice were subjected to a 5-day habituation protocol with each day consisting of several minutes of handling, temperature taken by rectal probe, intraperitoneal saline injection, and 30 min in the plethysmography chamber. Mice were then tested within 1 week of the last day of conditioning.
On the day of testing, mice were taken from their home cage and weighed, and rectal temperature was taken. Animals were then placed into an airtight, temperaturecontrolled (~32°C) plethysmography chamber and allowed to acclimate for at least 20 min in room air (21% O 2 /79% N 2 ) conditions. After acclimation and measurement under room air, the chamber gas was switched to a hypercapnic mixture of 5% CO 2 /21% O 2 /74% N 2 for 20 min. Chamber gas was then switched back to room air for 20 min. The mice were briefly removed for rectal temperature measurement and intraperitoneal injection of clozapine-N-oxide (CNO, National Institute of Mental Health Chemical Synthesis and Drug Supply Program) dissolved in saline (0.1 mg/mL) for an effective concentration 1 mg/kg. The animal was returned to the chamber for another 20 min of room air, 20 min of hypercapnia, and 20 min of room air. The animal was then removed from the chamber and rectal temperature was taken immediately afterwards and again 30 min after the termination of the experiment. The animal was placed in its own cage during these 30 min at the ambient room temperature (~23°C).

Plethysmography data collection and analysis
Plethysmography pressure changes were measured using a Validyne DP45 differential pressure transducer and CD15 carrier demodulator in comparison to a reference chamber and recorded with LabChart Pro in real time. Waveforms were analyzed offline using LabChart Pro and custom-written MATLAB code (hosted on Zenodo 10.5281/zenodo.5768587) to determine respiratory frequency (f R ), tidal volume (V T ) [52], ventilation (V E ), oxygen consumption (V O2 ), ventilatory equivalents for oxygen (V E /V O2 ), and pattern analysis. Respiratory waveforms were collected offline during periods when the animal was at rest and readings were free from movement artifacts. A minimum of 1-min cumulative data compiled from traces at least 10 s long from the last 10 min of a given experimental condition were analyzed. O 2 consumption was determined by comparing the gas composition between calibration in an empty chamber and live breathing using an AEI oxygen sensor and analyzer. Chamber temperature was constantly monitored and recorded using a ThermoWorks MicroThermo 2 with probe and was recorded with LabChart Pro in real time.

Plethysmography statistics
Results (f R , V T ,V E ,V O2 ,V E /V O2 ) for room air and hypercapnic data were compared between all cohorts using a linear mixed-effects regression model with animal type (experimental or control) and CNO administration (preor post-injection) as fixed effects and animal ID as a random effect. Temperature data was compared using a linear mixed-effects regression model with animal type (experimental or control) as a fixed effect. Residuals were independent and identically distributed as a normal distribution, which matches our model assumptions (Additional Files 3-5: Fig. S3-5). The residual plot for ventilatory equivalents of oxygen in room air for hM4D (Fig. 8) shows a slightly different distribution pattern (Additional File 4: Fig. S4). However, this is because the random effect coefficients are unusually small. This data still follows the normal and independent residuals assumption. Random effect coefficient distribution is not as easily assessed with so few datapoints and is not critical to statistical outcome if violated [92]. The p values reported correspond to statistical significance of the conditional interaction between animal type and CNO administration. A p value of < 0.05 was used to indicate statistical significance and standard error of the mean is shown on all charts.

Imaging of Twitch 2 in mouse preovulatory follicles
Preovulatory follicles were isolated and imaged as previously described [93]. Briefly, antral follicles were dissected from 23-to 26-day-old RR8; Bact_Cre mice. Follicles were cultured for 24-30 h on organotypic membranes (Millipore; cat. No. PICMORG50) in the presence of follicle-stimulating hormone. The follicle was held in a perfusion slide consisting of a plastic slide (ibidi) and a glass coverslip and assembled using silicon grease. The slide was constructed such that medium containing ovine LH (National Hormone and Peptide Program; 10 μg/mL) could be perfused through a 200-μm-deep channel holding the follicle. Temperature was maintained at 30-34°C, by use of a warm air blower (Nevtek). Preovulatory Follicles were imaged using a Zeiss Pascal confocal microscope with a × 40/1.2 NA objective. Images were collected every 10 s. Measurements were corrected for autofluorescence and for spectral bleed-through of CFP into the YFP channel. Ratios were calculated by dividing the mean CFP intensity in each region of interest by the mean YFP intensity. Data analysis was done using ImageJ and Excel software. Data is representative of 4 follicles.
PJH completed experiments, analyzed data, and prepared the manuscript. PGF completed experiments, analyzed data and prepared the manuscript. JP completed experiments, analyzed data, and prepared the manuscript. AC completed statistical analyses and interpretations and prepared the manuscript. JJS completed experiments, analyzed data, and prepared the manuscript. VKM assisted in the completion of experiments. PJZ completed experiments, analyzed data and prepared the manuscript. JRE completed experiments, analyzed data, and prepared the manuscript. GA statistical analyses and interpretations and prepared the manuscript. XJ completed experiments, analyzed data, and prepared the manuscript. BRA completed experiments, analyzed data, and prepared the manuscript. AST completed experiments, analyzed data and prepared the manuscript. M C-M assisted with experimental design and prepared the manuscript. RR conceptualized the study, performed experiments, analyzed data, and prepared the manuscript. All authors read and approved the final manuscript. Availability of data and materials All datasets, animals, and materials generated and/ or used in this study are publicly available (Addgene and GenBank, the Mutant Mouse Resource and Research Centers Supported by NIH (MMRRC)), or available from the corresponding author on reasonable request. All data generated or analyzed during this study are included in this published article and supplementary information files. Specifically, raw data can be found on Zenodo for electrophysiology experiments resulting in Fig. 7

Declarations
Ethics approval and consent to participate All experiments reported herein were conducted with explicit approval of and oversight by both Baylor College of Medicine Institutional Animal Care and Use Committee (IACUC) and University of Connecticut IACUC and abide by all state and national regulations regarding animal research for each work site.

Consent for publication
Not applicable.