Abstract
Several neurodevelopmental disorders are associated with increased mTOR activity that results in pathogenic neuronal dysmorphogenesis (i.e., soma and dendrite overgrowth), leading to circuit alterations associated with epilepsy and neurologic disabilities. Although an mTOR analog is approved for the treatment of epilepsy in one of these disorders, it has limited efficacy and is associated with a wide range of side effects. There is a need to develop novel agents for the treatment of mTOR-pathway related disorders. Here, we developed a medium-throughput phenotypic assay to test drug efficacy on neurite morphogenesis of mouse neurons in a hyperactive mTOR condition. Our assay involved in utero electroporation (IUE) of a selective population of cortical pyramidal neurons with a plasmid encoding the constitutively active mTOR activator, Rheb, and tdTomato. Labeled neurons from the somatosensory cortex (SSC) were cultured onto 96-well plates and fixed at various days in vitro or following Torin 1 treatment. Automated systems were used for image acquisition and neuron morphologic measurements. We validated our automated approach using traditional manual methods of neuron morphologic assessment. Both automated and manual analyses showed increased neurite length and complexity over time, and decreased neurite overgrowth and soma size with Torin 1. These data validate the accuracy of our automated approach that takes hours compared with weeks when using traditional manual methods. Taken together, this assay can be scaled to screen 32 compounds simultaneously in two weeks, highlighting its robustness and efficiency for medium-throughput screening of candidate therapeutics on a defined population of wild-type or diseased neurons.
Significance Statement
Preclinical studies and screens often rely on cell lines and traditional techniques that are time consuming and introduce human bias during analysis compared with automated methods. Some of these techniques include the manual tracing of cells to quantify morphologic changes in response to various treatments. We developed an assay that allows medium-throughput morphologic analysis in a selective population of neurons in vitro while expediting data collection and analysis through automation. This assay is modifiable and is applicable for a wide range of disease conditions.
Introduction
Several neurologic disorders are associated with abnormalities of neuronal development. One such group of neurodevelopmental disorders are mTORopathies that arise from gene variants leading to increased mTOR activity and epilepsy (Crino et al., 2002). mTOR hyperactivity occurs in neocortical pyramidal neurons and results in increased neuronal soma size, dendritic complexity, and alterations in connectivity, all of which contribute to epilepsy (Feliciano et al., 2013). Although one mTOR blocker, everolimus, is approved for the treatment of epilepsy in one of these disorders, tuberous sclerosis complex (Franz et al., 2021), this drug has severe side effects and seizures remain in most patients. While research has focused on developing novel drugs and treatment strategies, screening their efficacy in vivo has remained challenging considering that monitoring seizure activity is difficult and labor intensive. An alternative is to use a medium-throughput to high-throughput phenotypic assay to assess the effect of novel drugs or knock-down of specific molecules on the morphologic development of neurons expressing hyperactive mTOR.
The pharmaceutical industry as well as academic research laboratories have developed phenotypic assays to efficiently analyze the effect of many compounds on cells. These assays are critically important to investigate novel treatments in a wide range of disorders. For neurite outgrowth-based assays, various neuronal systems have been implemented. This includes the use of neuronal cell lines, such as SH-SY5Y cells, which have been differentiated into neurons and immunostained for neurite markers (Dravid et al., 2021; Schikora et al., 2021). Others have used nucleofected cortical neurons in 96-well plates (Blackmore et al., 2010), sparsely transfected primary neurons using automated imaging and analysis (Sharma et al., 2012), and human iPSC-derived neurons (primarily interneurons, with some Layer V cortical pyramidal-like neurons) immunostained for neurite markers (Sirenko et al., 2014; Sherman and Bang, 2018). In addition to using high-throughput automated imaging, studies have incorporated autonomous computational assays for quantitative analysis (Wu et al., 2010; Sharma et al., 2012; Schikora et al., 2021). Although some studies have achieved sparse labeling of neurons, which improves image analysis, none have labeled specific neuronal populations.
Genetic manipulation via in utero electroporation (IUE) offers a reliable method to selectively label distinct populations of cells. IUE is a technique involving the injection of DNA plasmid into the ventricles of embryos and applying an electrical current that creates transient pores in the cell membrane. This allows the DNA to enter the targeted neural progenitors lining the ventricle. For example, IUE at embryonic day (E)15 will specifically label Layer II/III cortical pyramidal neurons (Molyneaux et al., 2007).
Here, we provide proof of principle for the use of specific sequential strategies using IUE to induce a disease phenotype in a subset of murine cortical neurons, with the purpose of testing numerous compounds on diseased neurons in a medium-throughput manner using an IN Cell Analyzer 2200 Imaging System and a customized CellProfiler pipeline (Stirling et al., 2021).
Materials and Methods
Animals
Experiments were performed according to guidelines set forth by the Yale University Institutional Animal Care and Use Committee and National Institutes of Health Guide for the Care and Use of Laboratory Animals. E14 pregnant CD-1 mice were obtained from Charles River Laboratories. Mice were housed in littermate groups under pathogen-free conditions with a 12/12 h light/dark cycle.
In Utero electroporation
IUE was performed by aseptic survival surgery on a total of three E15 pregnant CD1 mice, 1 d postacclimation to their cage. Weight and approximate age were documented before surgery. Mice were injected with 10 mg/kg buprenorphine 30 min before anesthetization. Mice were placed into an anesthesia induction chamber and 3% isoflurane was applied. Following appropriate anesthesia, mice were placed onto a predisinfected surgical table with a heat source and attached to a nose cone with 2–3% isoflurane flow. Mouse incision regions were shaved, and the skin prepared for aseptic surgery. A midline ventral laparotomy was performed. Uterine horns were slowly removed with ringed forceps. DNA solution (∼1–1.5 μg/μl) containing a dual expression plasmid encoding a constitutively active Rheb, mutant Y35L with tdTomato (T2A) under the chicken β-actin promoter with CMV enhancer (CAG) promoter diluted in PBS with 0.1% fast green was prepared before surgery. DNA solution (1.5 μl) was injected into the lateral ventricle using a sterile pulled glass pipette. Electrodes (model 520; BTX) were soaked in 0.9% saline solution and placed on the heads of the embryos, and five, 40-V square pulses of 50-ms duration with 950-ms intervals were applied using a pulse generator (ECM830; BTX). After injection and electroporation, the incision was closed with absorbable sutures and skin closed with 9-mm autoclips. Mice were monitored for the following 2 d. The date of pup birth was recorded, and mice killed at the indicated time point neonatally.
Microdissection and cell culture
Hibernate E (ThermoFisher Scientific; A1247601) with 2% B27 (ThermoFisher Scientific; 17 504–044) and 2 mm GlutaMAX (ThermoFisher Scientific; 35 050–061) was first placed on ice, and papain digestion solution (Worthington; LK003176) was warmed to 37°C for 30 min before use. Mice were screened following IUE, and at postnatal day (P)1–P2 the brains of IUE-positive mice were collected and placed in a 60 mm-diameter Petri dish with cold Hibernate E solution. The positive somatosensory cortex (SSC) region of three to four brains were microdissected and cortices placed in 5 ml of cold Hibernate E solution on ice in a 15-ml conical tube. Dissected cortices were kept in Hibernate E solution on ice for the duration of all microdissections. Next, the Hibernate E solution was carefully removed, and 5 ml of prewarmed and activated papain digestion solution was added to enzymatically digest the cortices. The cap of the tube was sealed with Parafilm, and then incubated in a 37°C water bath for 15 min. At this step MEM solution (ThermoFisher Scientific; 11095-080) with 0.6% wt/vol of D-Glucose (Millipore Sigma; G8769) and 5% fetal bovine serum (ThermoFisher Scientific; 10082147) was prewarmed. After 15 min, the tube was removed from the water bath, sprayed with 70% ethanol, and the papain digestion solution was removed. A total of 10 ml of prewarmed Hibernate E solution was added. Once the cortices settled to the bottom of the tube, this process was repeated to wash out the remaining papain digestion solution. Finally, the supernatant was removed and replaced with 5 ml of prewarmed MEM solution. Cortices were gently triturated four times using a pipette bulb with a large bore sterile fire-polished glass Pasteur pipette (∼1 mm in diameter) and this was repeated with a slightly smaller bore sterile fire-polished glass Pasteur pipette. The tubes were then centrifuged at 400 × g for 2 min, the supernatant removed, and 3–8 ml of MEM solution added (volume depending on size of tissue microdissected and the number of brains used) followed by gentle trituration to resuspend cells. The MEM solution containing suspended cells was then transferred into a new 15-ml tube though a 40-μm cell strainer (ThermoFisher Scientific; 22363547).
The cells were counted using a hemocytometer and cultured in a glass bottom 96-well plate (Cellvis; P96-1.5H-N) pretreated with 50 μg/ml Poly-D-Lysine (ThermoFisher Scientific; A389890-01) and 30 μg/ml laminin (Millipore Sigma; 11243217001) at a density of 30 000 cells per well. Two hours later, half the media was replaced with neurobasal media (ThermoFisher Scientific; 21103-049) containing 2 mm GlutaMAX and 2% B27. For neuron maintenance, half of the media was replaced with an equal amount of fresh neurobasal solution daily. Cells were fixed with 4% paraformaldehyde (PFA) for 30 min at the conclusion of the experiment.
Immunocytochemistry
Cells were treated in blocking solution composed of 1× TBS (ThermoFisher Scientific; J60764) with 0.1% Tween 20 (Millipore Sigma; P7949) and 2% BSA (Millipore Sigma; A7906) for 1 h and incubated overnight in the primary antibody against RFP (1:500; Rockland; 600-401-379). Cells were again washed three times for 10 min each in wash buffer composed of 1× TBS with Tween 20, and again incubated overnight in anti-rabbit 555 secondary antibody (1:500; ThermoFisher Scientific; A32732).
Drug treatment
A total of 100 nm Torin 1 (Millipore Sigma; 475991) or an equimolar amount of DMSO (0.01%, controls) was added to neurons 2 h after seeding and replaced daily along with half the media. Drug treatments were done in four wells per condition. Cells were collected at 3 d in vitro (DIV) for drug treatments.
Automated imaging
Using a General Electric IN Cell Analyzer 2200 Imaging System, images were acquired in an automated fashion (Fig. 1). Fluorescent (CY3) images were acquired at a 10× magnification with a 450-ms exposure time in a grid like pattern for each well to avoid overlap.
Automated neurite and soma measurements
For automated neurite measurements, a custom pipeline in CellProfiler was developed to measure and export raw data containing both the total neurite length per cell and the number of terminal branches per cell. The settings are as follows: neurites suppressed for identification of soma as primary object (diameter, 12–45 pixels), neurites re-enhanced, and neurites identified as secondary objects based on primary object (soma) location. To convert the image to binary and skeletonize the neurites, we used an adaptive threshold strategy with a smoothing scale of 0.5, a correction factor of 2, a lower and upper bound on threshold of 0.05–1.0, and 10 pixels as the adaptive window size. Finally, the skeletonized neurites were automatically measured using the “measure object skeleton” and “measure object size shape” modules, and readouts such as total neurite length per cell and the number of terminal branches were exported to excel for statistical analysis by the “export to spreadsheet” module.
For automated soma size measurements, another CellProfiler pipeline was generated, where neurites were suppressed using the enhance or suppress features module, and primary objects (soma) were identified. Identification parameters include a 10- to 50-pixel diameter, an adaptive minimum cross entropy threshold with a smoothing scale of 2, a threshold correction factor of 1, a lower and upper bound on threshold of 0.35–1.0, and the shape method to distinguish clumped objects. Finally, the “measure object size shape” module was used to automatically measure the pixel area of primary objects (soma) and the results were exported to excel for statistical analysis. Pipeline files can be found on the CellProfiler website for download.
Because of the high number of cells analyzed and the automated nature of the measurements, the comparison between groups (particularly involving DIV9 neurons, which had very extensive neurite arbors) are more easily visualized if the y-axis is truncated to omit one to eight upper limit cells from the graph. These were not omitted from analysis or statistics, and only represent a small fraction of the 375 (Fig. 2) or 176 (Fig. 3) cells analyzed. Number of visually omitted individual points per graph for Figure 2 by y-axis truncation: (DIV3 = 0 for all graphs), Figure 2D, DIV6 = 1, DIV9 = 7; Figure 2E, DIV6 = 2, DIV9 = 6; Figure 2F, DIV6 = 0, DIV9 = 3; Figure 2G, DIV6 = 0, DIV9 = 6. Number of visually omitted individual points per graph for Figure 3 by y-axis truncation: (Torin 1 = 0 for all graphs), Figure 3C, DMSO = 8; Figure 3D, DMSO = 4; Figure 3E, DMSO = 2; Figure 3F, DMSO = 6; Figure 3G, DMSO = 3; Figure 3H, DMSO = 0.
Manual neurite and soma measurements
Images were uploaded to FIJI (ImageJ 1.53q) and neurites were traced using the Simple Neurite Tracer plug-in. Total length and number of terminal branchpoints per cell were extracted from FIJI based on the manual tracing. The freehand tracing tool in FIJI allowed for quantification of cell size by tracing the soma of tdTomato+ neurons and extracting the area.
Statistical analysis
Statistical analyses using unpaired two-tailed Student’s t test and ordinary one-way ANOVA with Tukey’s post hoc test were performed using GraphPad Prism version 9.3.1.
Results
Figure 1A illustrates the experimental paradigm from neuron labeling to analysis that constitute our new phenotypic assay. We performed IUE of one plasmid encoding both a constitutively active Y35L mutant Rheb and tdTomato into the lateral ventricle of E15 mouse embryos. We targeted E15 neural stem cells (radial glia) of the somatosensory cortex (SSC) that generate Layer II/III pyramidal neurons (Molyneaux et al., 2007; Lin et al., 2016; Sokolov et al., 2018). RhebY35L is known to increase mTOR activity, resulting in neuron dysmorphogenesis including increased soma size and dendrite overgrowth (Zhao et al., 2019). Four days later, the SSC containing tdTomato-positive neurons were microdissected at postnatal day 0–1 and cultured onto a glass bottom 96-well plate at a density of 30,000 cells per well. This approach led to sparse tdTomato-positive neuron plating (Fig. 2A–C). On the indicated days in vitro (DIV), neurons were fixed, and immunocytochemistry was performed for tdTomato to amplify neurite labeling. Next, neurons were imaged in an automated fashion using an IN Cell Analyzer 2200 Imaging System. Following image acquisition, morphologic measurements were performed by feeding the IN Cell Analyzer 2200 images into a pipeline of customized modules in CellProfiler. These modules include various postprocessing steps such as identifying structures and adaptive thresholding (detailed in Materials and Methods; Fig. 1B–E). To validate this phenotypic assay, we examined neurite outgrowth over time and tested the impact of the mTOR blocker, Torin 1, on neurite outgrowth and soma size in a hyperactive mTOR condition. This experimental paradigm was completed in under two weeks including analysis.
Automatic versus manual assessment of neurite outgrowth in vitro over time
To measure baseline morphologic changes in RhebY35L neurons over time, we compared neurite outgrowth at DIV3, DIV6, and DIV9. Using the protocol outlined in Figure 1, significant increases in neurite length and complexity were identified over time. Total neurite length per cell increased by 99% from DIV3 to DIV6 (p < 0.01) and 74% from DIV6 to DIV9 (p < 0.0001; N = 3 wells for each time point for all experiments; Fig. 2A–D). Additionally, the number of terminal branches per cell, an indicator of neurite complexity, increased by 79% from DIV3 to DIV6 (p < 0.05) and 93% from DIV6 to DIV9 (p < 0.0001; Fig. 2E). To validate our automated analysis, we manually traced and measured the neurites of labeled neurons using FIJI software. This confirmed the changes in neurite morphology found using the automated pipeline, with an increase in the total neurite length per cell (DIV3-to-DIV6: 127%, p < 0.001; and DIV6-to-DIV9: 79%, p < 0.0001) and the number of terminal branches per cell (DIV3-to-DIV6: 57%, p < 0.05; and DIV6-to-DIV9: 71%, p < 0.0001; Fig. 2F,G). These data suggest that the automated analysis accurately recapitulates manual analysis and is sensitive and robust enough to detect statistically significant differences in neurite length and complexity over time.
Automatic and manual assessment of Torin 1-induced morphologic changes in RhebY35L neurons
Considering mTOR’s well-documented role in dendrite outgrowth, we chose the mTOR inhibitor Torin 1 to assess the efficacy of using our pipeline to monitor drug-induced morphologic changes. We repeated the protocol outlined in Figure 1 and treated cultured RhebY35L neurons daily with 100 nm Torin 1 or vehicle until DIV3 (equimolar DMSO, 0.01%). We chose to treat for 3 d considering neurites are not overly complex at this time point, allowing less neurite overlap and more accurate automated detection. We found that Torin 1 reduced total neurite length per cell by 66% (p < 0.001; Fig. 3A–C). Torin 1 also decreased the number of terminal branches per cell by 54% (p < 0.0001; Fig. 3D). As mTOR activity influences cell size, we used a second modified CellProfiler pipeline to automatically measure soma size and found a 29% reduction following Torin 1 treatment (p < 0.01; Fig. 3E). To confirm that these changes were accurately measured, we manually traced and measured neurites using FIJI software. The manual analysis reflected the results of the automated analysis, showing Torin 1 reduced the total neurite length per cell (71% decrease; p < 0.0001), the number of terminal branches per cell (62% decrease; p < 0.001), and soma size (30% decrease; p < 0.01; Fig. 3F–H). Thus, our automated phenotypic assay is a viable alternative to manual analysis to efficiently analyze the morphologic effect of drugs on sparsely labeled neurons.
Discussion
Here, we describe a versatile novel medium-throughput assay to investigate morphologic changes in primary mouse cortical neurons in an automated fashion. As opposed to cell lines, primary neuronal cultures offer a more clinically relevant model for studying neurologic disorders. We also used IUE, which allows us to express a plasmid of our choice that encodes constitutively active Rheb to increase mTOR activity, as reported in several neurodevelopmental disorders, and assess the impact of drug treatment on neuron morphology.
We validated our experimental paradigm by successfully measuring and quantifying changes in neurite outgrowth over time and following Torin 1 treatment using the CellProfiler pipeline. As expected, our assay detected increased neurite length and complexity over time and Torin 1 reduced neurite length, complexity, and soma size compared with vehicle treated neurons. To further confirm the accuracy of the automated measurements, we performed traditional manual tracings using FIJI Simple Neurite Tracer (Schindelin et al., 2012). The manual analysis gave similar results to those with the automated analysis. However, while manual tracing took over a week, our automated approach took hours. It is nevertheless important to note the discrepancy in mean values recorded using the automated versus the manual method. The automated method underestimated the total length and complexity compared with the manual tracing method. This discrepancy is because of the thresholding applied using automated analysis. A threshold is applied to avoid false neurite labeling (e.g., because of neurite crossing), resulting in the omission of fainter distal portions of the neurites. However, the changes between experimental conditions remain consistent with automated and manual analyses. Collectively, our automated paradigm is a viable alternative to the traditional manual analysis for quantifying neuron morphology in a more efficient manner.
To label cells for automated neurite imaging and tracing, we used IUE. Other labeling techniques can be combined with our approach such as sparse transfection, nucleofection, lentiviral infection of cultured rodent neurons or human iPSC-derived neurons, and AAV infection in vivo. It can also be combined with the use of transgenic mice or rats. All these techniques allow sparse labeling and the targeting of neurons or glia instead of using cell lines differentiated into a neuron-like state that do not express all neuronal or glial properties (Dravid et al., 2021; Schikora et al., 2021). Drug screening may thus be more clinically translatable. The major advantage of IUE over plasmid transfection, nucleofection, and lentivirus infection in vitro, is the specificity of the labeled cell type and region targeted without the need for using specific promoters. Here, we targeted E15 cortical radial glia that generate Layer II/III pyramidal neurons (Molyneaux et al., 2007; Lin et al., 2016; Sokolov et al., 2018). IUE at an earlier time point would allow targeting of deep layer pyramidal neurons, and IUE at later timepoints (E18) would preferentially label astrocytes (Molyneaux et al., 2007). Targeting the ganglionic eminence would label interneurons (Anderson et al., 1997; Kepecs and Fishell, 2014). Another advantage of IUE over the above mentioned in vitro labeling approaches is that cell viability is increased considering IUE occurs days before culturing the neurons (Blackmore et al., 2010; Sharma et al., 2012). Compared with transgenic mice, IUE is simply more versatile in expressing plasmids that encode any protein of interest or knock-down systems (shRNA or CRISPR/Cas9) without the need to generate double or triple transgenic mice. Finally, IUE is comparable to in vivo AAV considering that the cost of AAV production has decreased. Specific promoters can allow AAV to be expressed in excitatory versus inhibitory neurons (interneurons) versus different types of glia. The promoters may not be yet specific enough to achieve labeling of Layer II/III versus IV or V pyramidal neurons, but most screens may not need this level of specificity.
Our assay, however, has a couple of minor limitations. The CellProfiler algorithm does not efficiently differentiate between axons and dendrites. This did not affect data from our studies since mTOR affects both axon and dendrite growth (Kumar et al., 2005; Gong et al., 2015). For future studies examining the effect of drugs on either axon or dendrite growth, it would be important to perform immunocytochemistry for specific neurite markers, such as MAP2 for dendrites and tau for axons, and optimize the CellProfiler modules to only measure MAP2 or tau positive objects. Another limitation is that IUE requires technical expertise in handling and manipulating embryonic tissues as well as precise injection capability. A level of consistency with the IUE and microdissection procedure is essential for optimally labeling cultured neurons.
In conclusion, we have developed a new phenotypic assay that can screen the morphologic effects of compounds on specific neuronal populations in a robust and selective fashion in vitro. The newly designed pipelines in CellProfiler can be adapted to analyze the morphology of different populations of neurons and glia as well as non-neuronal cells in vitro. Finally, the assay is medium-throughput and can provide clinically relevant data on ∼32 compounds in under two weeks.
Acknowledgments
Acknowledgments: We thank the Yale Center for Molecular Discovery for the helpful conversations and expertise, in particular Laura Abriola, for operating the IN Cell Analyzer for image acquisition.
Footnotes
The authors declare no competing financial interests.
This work was supported by the National Institutes of Health/National Institute of Neurological Disorders and Stroke Grant R01 NS093704 and by the Swebilius Foundation (A.B.).
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.