A fully automated high-throughput workflow for 3D-based chemical screening in human midbrain organoids
Abstract
Three-dimensional (3D) culture systems have fueled hopes to bring about the next generation of more physiologically relevant high-throughput screens (HTS). However, current protocols yield either complex but highly heterogeneous aggregates (‘organoids’) or 3D structures with less physiological relevance (‘spheroids’). Here, we present a scalable, HTS-compatible workflow for the automated generation, maintenance, and optical analysis of human midbrain organoids in standard 96-well-plates. The resulting organoids possess a highly homogeneous morphology, size, global gene expression, cellular composition, and structure. They present significant features of the human midbrain and display spontaneous aggregate-wide synchronized neural activity. By automating the entire workflow from generation to analysis, we enhance the intra- and inter-batch reproducibility as demonstrated via RNA sequencing and quantitative whole mount high-content imaging. This allows assessing drug effects at the single-cell level within a complex 3D cell environment in a fully automated HTS workflow.
eLife digest
In 1907, the American zoologist Ross Granville Harrison developed the first technique to artificially grow animal cells outside the body in a liquid medium. Cells are still grown in much the same way in modern laboratories: a single layer of cells is placed in a warm incubator with nutrient-rich broth. These cell layers are often used to test new drugs, but they cannot recapitulate the complexity of a real organ made from multiple cell types within a living, breathing human body.
Growing three-dimensional miniature organs or 'organoids' that behave in a similar way to real organs is the next step towards creating better platforms for drug screening, but there are several difficulties inherent to this process. For one thing, it is hard to recreate the multitude of cell types that make up an organ. For another, the cells that do grow often fail to connect and communicate with each other in biologically realistic ways. It is also tough to grow a large number of organoids that all behave in the same way, making it hard to know whether a particular drug works or whether it is just being tested on a 'good' organoid.
Renner et al. have been able to overcome these issues by using robotic technology to create thousands of identical, mid-brain organoids from human cells in the lab. The robots perform a series of precisely controlled tasks – including dispensing the initial cells into wells, feeding organoids as they grow and testing them at different stages of development. These mini-brains, which are the size of the head of a pin, mimic the part of the brain where Parkinson's disease first manifests. They can be used to test new drugs for Parkinson's, and to better understand the biology of the brain. Perhaps more importantly, other types of organoids can be created using the same technique to model diseases that affect other areas of the brain, or other organs altogether. For example, Renner et al. also generated forebrain organoids using an automated approach for both generation and analysis.
This research, which shows that organoids can be grown and tested in a fully automated, reproducible and scalable way, creates a platform to quickly, cheaply and easily test thousands of drugs for Parkinson's and other difficult-to-treat diseases in a human setting. This approach has the potential to reduce research waste by increasing the chances that a drug that works in the lab will also ultimately work in a patient; and reduce animal experiments, as drugs that do not work in human tissues will not proceed to animal testing.
Introduction
A number of uniquely human diseases, including Parkinson’s disease, would greatly benefit from a comprehensive human cellular in vitro model that recapitulates key characteristics of midbrain tissues in a high-throughput-compatible format. Three-dimensional (3D) cell culture in general and the ability to generate organ-like aggregates (‘organoids’) in particular have found a rapid following over the past few years (Sato et al., 2009; Eiraku et al., 2011; Nakano et al., 2012; Lancaster et al., 2013; Quadrato et al., 2017; Paşca et al., 2015; Iefremova et al., 2017; Takasato et al., 2015; Dye et al., 2015; Takebe et al., 2013; McCracken et al., 2014; Pașca, 2018) due to their potential to mimic cellular niches more closely than 2D cell cultures. These approaches promise to develop next-generation high-throughput screens (HTS) that can provide more relevant predictions of drug efficacy and toxicity (Fatehullah et al., 2016; Ranga et al., 2014; Fang and Eglen, 2017; Dutta et al., 2017; Ho et al., 2018; Chen et al., 2018; Friese et al., 2019) as they may allow better modeling of pathologies with complex interactions of several cell types in specific cellular niches (Qian et al., 2016; Mariani et al., 2015; Ogawa et al., 2015; Verissimo et al., 2016; Vlachogiannis et al., 2018; Czerniecki et al., 2018).
3D culture in the form of spheroids has long been established, especially in the field of cancer biology (Sutherland et al., 1971), and used in various HTS applications (Kelm et al., 2003; Senkowski et al., 2015; Wenzel et al., 2014; Kenny et al., 2015). While these model systems are already more complex and potentially more physiologically relevant than 2D culture (Pickl and Ries, 2009), they display a much simpler, less organ-like 3D cell organization (Fang and Eglen, 2017) and do not mimic functional features of the organ as broadly and as closely as organoid tissues. Moreover, many of the 3D-based screens performed so far depend on whole-aggregate-based readouts such as size, morphology, and cell viability (Vlachogiannis et al., 2018; Ivanov et al., 2014; Friedrich et al., 2009; Ivanov et al., 2015; Hou et al., 2018; Kang et al., 2015), which make it challenging to gain mechanistic insights into cells or sub-populations of cells in the context of their niches.
In contrast, complex organoids have emerged as a promising research tool due to their unique resemblance to human tissues, defined by an organ-like architecture composed of different tissue-specific cell types and the capability to mimic organ functions (Lancaster and Knoblich, 2014). One structure of particular interest in the context of disease modeling and drug development is the midbrain due to its role in the highly prevalent Parkinson’s disease (affecting 2–3% of the population aged 65 and above Poewe et al., 2017) and many other developmental disorders (Barkovich, 2012; Doherty et al., 2013). While the current state of the art allows the generation of different types of midbrain organoids (Qian et al., 2016; Jo et al., 2016; Monzel et al., 2017), the rigorous standardized organoid production and quantification methods needed for high-throughput screening have been elusive. The established protocols tend not to focus on scalable, homogeneous organoids with quantitatively predictable morphology, cellular composition, and local cell organization. Obligatory extensive manual handling often including cumbersome matrix embedding steps render them challenging for scale-up (Tong et al., 2018). Furthermore, common analysis methods (e.g. sectioning and immunostaining, RNA sequencing) do not scale well for HTS applications.
Here, we present a fully automated workflow in a standard 96-well format that combines generation, maintenance, whole mount immunostaining, tissue clearing, and high-content imaging of automated midbrain organoids (AMOs) (see Figure 1a). The resulting AMOs are similar to published midbrain organoids with regard to their expression of midbrain-specific markers and cell populations, yet maintain a reproducible homogeneous phenotype. They mimic relevant organ function in the form of spontaneous, highly synchronized neural activity indicating functional cellular coupling across the entire AMO. Their high homogeneity, reproducibility, culture format, and fast development of approximately one month render them ideal for high-throughput screening applications. Moreover, our combined whole mount immunostaining and clearing workflow abolishes the need for labor-intensive tissue sections and allows for quantitative whole mount high-content analysis of entire organoids with single-cell resolution. Our automation of the entire workflow from seeding to analysis in standard plates allows for easy scale-up and implementation into existing screening facilities.
Results
Automation enables high-throughput-compatible production of homogenous midbrain organoids
Screening applications require biological systems that operate within predictable physiological parameters. In order to limit cellular heterogeneity during differentiation, we produced human AMOs starting from small molecule neural precursor cells (smNPCs) (Reinhardt et al., 2013a), which in turn originate from pluripotent stem cells (PSCs). The neural-restricted developmental potential of smNPCs still allows the self-organization required for the formation of a 3D architecture (Monzel et al., 2017; Di Lullo and Kriegstein, 2017) during differentiation toward a midbrain fate. To further reduce batch-to-batch variability, we also omitted matrigel embedding and standardized mechanical stresses by using an automated liquid handling system (ALHS). Starting from seeding the organoids, all following steps including maintenance, fixation, whole mount staining, and clearing are performed in a fully scalable automated fashion using a 96-channel-pipetting head in a robotic ALHS (see Figure 1a). The resulting AMOs show little intra- and inter-batch variability in size distribution (see Figure 1b, average coefficient of variation (CV) within one batch 3.56%; min 2.2%, max 5.6%), morphology (see Figure 1c), and cellular composition and organization (see Figure 2), making them ideal for HTS-approaches. Furthermore, our workflow generates one aggregate per well, maintained independently from others, thus minimizing batch effects due to paracrine signaling observed in bioreactor-based strategies (Quadrato et al., 2017). If paracrine signaling is desired, our workflow can easily accommodate several aggregates per well. The fully automated workflow operates with very high efficiency, retaining 99.7% (standard deviation 0.7%) of samples for automated seeding, aggregation, and maturation steps over 30 days and 96.5% (standard deviation 3.1%) of samples for fixation, whole mount staining, clearing, and transfer to flat bottom imaging plates over 12 days. Lastly, 6.1% (standard deviation 1.3%) of these samples are rejected during high-content imaging for presence of dust, damage, or fibers (see Figure 1d and Supplementary file 1 for source data).
Automated midbrain organoids express typical neural and midbrain markers and show structural organization
In order to characterize protein localization in our AMOs (>600 μm diameter) and assess the efficiency of their neural/midbrain differentiation at a cellular resolution and in a HTS-compatible manner, we adapted an extended 3D-staining protocol (Lee et al., 2016) and combined it with benzyl alcohol and benzyl benzoate (BABB)-based tissue clearing (Dent et al., 1989). BABB-based clearing proved to be both the fastest and most efficient method in a comparison of different clearing protocols (see Figure 1—figure supplement 1). The combination of whole mount staining and clearing allows the 3D reconstruction of entire AMOs via confocal imaging and enables further detailed 3D quantification and analysis, for example tracing of neurites throughout the whole aggregate, which cannot be performed using typical tissue sectioning procedures (see Video 1).
The immunostaining results are depicted as either single confocal optical slices (see Figure 2a–d,f/g) or maximum intensity projections (MIP, see Figure 2e). Already at day 25, the AMOs contained large numbers of neurons as indicated by the expression of Map2 (Shafit-Zagardo and Kalcheva, 1998) (also see Video 1), β-tubulin III (TUBB3) (Leandro-García et al., 2010) (see Figure 2e), and doublecortin (Gleeson et al., 1999) (DCX, see Figure 2c/d). Presence of high levels of tyrosine hydroxylase (TH, Figure 2a/b, also see Figure 2—figure supplements 1 and 2), the rate-limiting enzyme in dopamine synthesis (Nagatsu, 1995), as well as the expression of the transcription factors Foxa2, Lmx1a, Nurr1, and Pitx3 (Hegarty et al., 2013) (see Figure 2—figure supplements 1 and 2), are consistent with differentiation toward a dopaminergic midbrain fate. While other neuronal subtypes, specifically GABAergic (vGAT) and glutamatergic (vGLUT1) neurons, are present in AMOs (see Figure 2—figure supplement 1), their abundance is low compared to dopaminergic neurons. As commonly seen in all 3D neural cultures, AMOs retain a population of neural precursors identified by the expression of Sox2 (Ellis et al., 2004) (see Video 1, Figure 2a/b, and Figure 2—figure supplement 2), Brn2 (Dominguez et al., 2013; Figure 2c/d and Figure 2—figure supplement 2), and the more general neural marker nestin (Hendrickson et al., 2011) (see Figure 2a/b and Figure 2—figure supplement 2).
Over time, AMOs matured further. Expression of synapsin (Thiel, 1993) (see Figure 2—figure supplements 1 and 2) as well as the presynaptic marker synaptophysin and postsynaptic marker homer (Tadokoro et al., 1999) frequently colocalized with each other on Map2-positive neurites (see Figure 2f) and indicated the presence of synapses. Since gliogenesis follows neurogenesis in vivo (Miller and Gauthier, 2007), we expected the emergence of astrocytes after the initial formation of neurons. Consistently, AMOs contained GFAP and S100b double-positive astrocytes (Götz et al., 2015) at later stages (see Figure 2g).
In the cortex, neurons form cortical layers with distinct markers. Cortical 3D models recapitulate this layer organization to a degree (Lancaster et al., 2013; Paşca et al., 2015; Qian et al., 2016; Mariani et al., 2015; Bhaduri et al., 2020). In contrast, the midbrain does not possess the typical layer organization of the cortex, hence published midbrain organoids are devoid of cortex-like layers (Qian et al., 2016; Monzel et al., 2017). Due to their self-organizing nature, typical published midbrain organoids initially form random local subdomains of organized tissue within the bulk of the organoid, often in the form of rosettes (Jo et al., 2016), a hallmark of the very early stages of neural development (Perrier et al., 2004; Elkabetz et al., 2008). This makes their morphology harder to predict within batches and leads to heterogeneity rendering screening strategies more challenging. We have optimized our AMOs to not form distinct random local subdomains; rather, the different cell types within the AMOs (i.e. neurons, astrocytes, and neural progenitors) self-organize into different concentric zones with distinct cellular orientations spanning the entire organoid (see Figure 2a–c). The outermost layer of the AMOs contains few nuclei with a dense, circumferentially oriented layer of TH+/nestin+/DCX+ cell processes. Cellular orientation changes in the underlying zone closer to the core, with TH+ dopaminergic and DCX+ neurons showing a clear radial alignment (see Figure 2b/c). The next zone, separating this region of radially organized neurons and the core, contains circumferentially oriented DCX+ neurons and few Brn2+ neural precursors (see Figure 2b/c). The core itself includes mostly neural precursors and few neurons. Within each zone, the different cell types are homogeneously distributed. In the context of HTS-compatibility, this homogeneity is advantageous, as this trait allows us to extrapolate data even from single medial confocal slices, drastically shortening acquisition times (see Figure 8—figure supplement 1). This inter-sample similarity can provide a uniform baseline for compound testing and render self-organized 3D human neural cell niches amenable to HTS strategies.
Ultrastructural analysis of AMOs (see Figure 2—figure supplement 3) supported the immunofluorescence data, revealing a dense 3D cell architecture consistent with neuronal cell bodies surrounded by nerve fibers. Analyzing the nerve fibers at a higher magnification revealed regular-spaced neurofilaments and microtubules. Moreover, vesicles with the characteristic size and localization of synaptic vesicles are frequently found within these nerve fibers.
Further quantitative real-time PCR (qPCR) analysis demonstrated increasing expression levels of various neural (DCX, Map2, NEFL, NeuN, TBR2, TUBB3, Syt1), midbrain (TH, NURR1, NKX6-1, EN1, GIRK2, AADC), and glia-specific (GLAST, MBP, S100b) markers at different developmental stages with concomitant decreases in neural precursor markers (Brn2, nestin, Pax6, Sox1, Sox2), confirming neural maturation toward a midbrain identity over time (Figure 3). Finally, we replicated these immunostaining and gene expression results with a second independent cell line (‘AMO line 2’, see Figure 2—figure supplement 2 and Figure 3—figure supplement 1), demonstrating the applicability of our workflow to cells with different genetic backgrounds and origin.
Calcium imaging and electrophysiological analysis reveal spontaneous and synchronized activity throughout entire organoids
To assess functional coupling of individual cells within the AMOs we first performed Fluo-4 acetoxymethyl ester (AM)-based calcium imaging, which can be used as a readout for spiking activity of neurons (Grienberger and Konnerth, 2012). In addition to frequent spontaneous activity of individual cells, we observed aggregate-wide synchronous and periodic calcium spikes (see Video 2) in all analyzed AMOs. To characterize this behavior further, we defined different regions of interest (ROIs) and assessed the change in fluorescence intensity over time in each region (see Figure 4). Measuring the entire AMO reveals two consecutive spikes in Fluo-4 brightness, with a period of approximately 30 s (see Figure 4a). When we subdivided the measured area into four quadrants, we observed synchronized spiking activity in all four resulting ROIs (see Figure 4b). This parallel activity pattern could be found at many structural levels of the AMO, even for single cells (see Figure 4c/e). Changing the time scale revealed additional levels of synchronicity between selected single cells, in addition to aggregate-wide spikes (see Figure 4d/e). Considering the calcium-imaging analyses along with the existence of synaptic vesicles on the ultrastructural level (see Figure 2—figure supplement 3), the verification of synapses via immunostaining (see Figure 2f and Figure 2—figure supplements 1 and 2) as well as synaptotagmin 1 (Syt1) via qPCR (see Figure 3 and Figure 3—figure supplement 1), our results support the presence of functionally coupled and spontaneously active neurons within the AMOs. The synchronous spiking patterns suggest that not only a small number of neurons but, in fact, the entire aggregate may be functionally connected. Large-scale synchronous bursting behavior can also be observed in several developing brain regions in vivo (Ben-Ari, 2001) and in brain slices in vitro (Silva et al., 1991).
Given the high reproducibility of synchronized calcium activity across all tested samples, we decided to evaluate the feasibility of using it as a functional readout in screening settings. Thus, we performed Fluo-4 AM-based calcium imaging on younger (day 35) AMOs and measured the resulting fluorescence signals on a standard plate reader. While they showed a shorter periodicity than the older samples, all tested AMOs displayed distinct peaks in fluorescence intensity resembling the synchronous activity patterns seen during spinning disk microscopy analysis (see Figure 4f), also in a second independent AMO line (see Figure 2—figure supplement 2). Treatment with the known calcium channel blocker cobalt(II) chloride completely abolished these peaks (see Figure 4g). Since the synchronous calcium activity of the AMOs and its modulation by inhibitors can be measured easily via HTS-friendly standard plate readers or specialized FLIPR Ca imagers (Sirenko et al., 2019), AMOs may be a promising 3D model of human neural activity that allows directly assessing midbrain related organ functions in HTS.
Multielectrode array (MEA) measurements revealed spontaneous electrical activity in 35-day-old AMOs (see Figure 4h/i). The field potentials of several electrodes in proximity to a 35-day-old AMO (see Figure 4i) oscillated in synchrony over time with two concurrent main frequencies at 1 Hz and 14 Hz (see Figure 4j, upper and lower panels). Figure 4k shows a Fast Fourier Transformation (FFT) based on the data shown in Figure 4j. Uncoordinated single-cell activity can hardly account for such robust and spatially long-ranging electric field oscillations. Rather, the concurrent and covariant signals at disparate electrodes support a widespread, synchronized electrical activity encompassing the entire AMO. Taken together, this data further supports the functional coupling of entire AMOs indicated by the calcium imaging experiments.
Finally, we characterized the electrophysiological properties of single cells from the AMOs using voltage patch-clamping. A stepwise increase of the membrane holding potential from −70 to +60 mV with 10 mV increments elicited transmembrane currents that consisted of a fast-activating, fast-inactivating inward current followed by a slower activating, slowly deactivating outward current ranging from a few hundred pA to several nA (see Figure 4l). The I–V curves of both currents are typical for sodium inward and potassium outward currents through voltage-gated channels (see Figure 4m; Reinhardt et al., 2013a; Simard et al., 1993; Cummins et al., 1994; Reinhardt et al., 2013b). Furthermore, the current-clamp recordings demonstrated that these cells generated action potentials (APs) in response to current injections (see Figure 4n). The average membrane potential of the recorded cells was −41.9 ± 15.2 mV (n = 29). These typical excitable, neuron-like electrophysiological properties could be detected as early as day 25 (see Figure 4—figure supplement 1) and in 29 of the 62 recorded cells. The rest of the cells possessed only outward currents of a few hundred pA by stepping to +60 mV (see Figure 4—figure supplement 1) and were unable to generate APs in response to current injections. These may represent other cell types present in AMOs like astrocytes and neural precursors.
RNA sequencing supports differentiation toward a human midbrain-like fate and homogeneous, predictable gene expression of automated midbrain organoids
To characterize AMOs on the level of global gene expression, we performed RNA sequencing of single organoids from three independent batches of AMOs (i.e. cells were separately thawed, seeded, and cultured) and compared the results with published RNA sequencing data sets of primary human tissues (Roost et al., 2015) and established midbrain organoids (Jo et al., 2016). Consistent with successful neural differentiation, AMOs were most similar to the brain and spinal cord in a panel with data from 21 human fetal tissues (see Figure 5a). Moreover, AMOs also correlated well with published data sets from different midbrain(-like) samples including primary human tissue (see Figure 5b). On a global gene expression level, AMOs more closely resembled the primary human midbrain samples than published midbrain organoids (AMOs = correlation 0.78, published midbrain organoids = 0.72, see Figure 5b). We also included publicly available RNAseq data from three prenatal human cortex samples (Jaffe et al., 2015) in the comparison as non-midbrain controls. As expected, the cortical samples showed high correlation with each other but less with the midbrain samples (see Figure 5b; for further comparisons between AMOs and cerebral organoids via the protocol by Lancaster et al., 2013. see Figure 5—figure supplement 1). Taken together, AMOs resemble published midbrain organoids as well as primary human midbrain tissue at the level of global gene expression.
Since homogeneity and reproducibility are crucial for screening applications, we next examined the variance of AMOs on the gene expression level and compared it to that of published midbrain organoids (Jo et al., 2016). This revealed that AMOs were consistently more reproducible within and between different batches than current midbrain protocols, as illustrated by the principal component analysis (PCA) plotted in Figure 5c. The AMOs from three independently thawed and cultured batches (n = 64 separately processed single organoids in total) clustered much more closely together than the published midbrain organoids (n = 6). This further underlines the utility of AMOs as a 3D cellular platform for HTS strategies. In screening settings, the wells at the edges of plates often display different readouts than those located toward the center of the plate (‘edge-effects’) (Malo et al., 2006). Therefore, we sequenced half of a 96-well plate for one AMO batch and tested for differences resulting from well location within the plate (group ‘one inside’ = center of the plate vs. ‘one outside’ = edge in Figure 5a–c). Importantly, in the PCA plot the AMOs clustered independently of their position on the plate (Figure 5c) and the groups also showed no apparent differences in any of the other analyses (Figure 5a/b), indicating that AMOs exhibit no measurable edge effects at the global gene expression level and further substantiating the high reproducibility of our protocol.
To further investigate the differences between AMOs and established midbrain organoid protocols, we performed gene ontology (GO) (Ashburner et al., 2000; Supek et al., 2011) analysis of the genes significantly upregulated (padj. < 0.05) in AMOs compared to the previously used published midbrain organoids (Jo et al., 2016). This analysis yielded almost exclusively GO terms connected to neuronal and synaptic activity (Figure 5d; for a complete list of GO terms see Supplementary file 2). Consistent with the previously described synchronous activity patterns (see Figure 4), this further illustrates the physiological relevance and efficient neural differentiation of AMOs.
Automated whole mount immunostaining is highly quantitative and reveals homogeneity of automated midbrain organoids
While immunofluorescence-based screening-compatible techniques of whole 3D aggregates have been reported, they can only detect cells in the outer layers of large organoids (Vergara et al., 2017), or they use small aggregates of approximately 100 μm diameter (Verissimo et al., 2016) or cystic organoids (Czerniecki et al., 2018), both of which can be penetrated by antibodies and fluorescence illumination more easily. In contrast, our workflow is custom-tailored for automation and allows the quantification of entire dense, large-scale aggregates (>800 μm diameter) with single-cell resolution and high sensitivity, as highlighted by a dose-response assay for 3D cellular detection (see Figure 6a). We mixed cells labeled with CellTracker deep red dye with unlabeled cells at known proportions, aggregated them to sizes similar or exceeding that of AMOs (750 μm and 950 μm, see Figure 6a), cleared them, and then analyzed them on a confocal high-content imaging system. The resulting relationship between the amount of tracked cells and measured brightness was highly linear (R2 >0.99), illustrating the quantitative nature of our optical HTS 3D whole mount analysis workflow.
Next, we demonstrated the homogeneity of AMOs at the protein level. A fully automated 96-well based whole mount optical analysis (see Figure 6b left) illustrated the ability to detect both abundant filamentous structures (neural marker Map2) and nuclear markers (Sox2) in a HTS-compatible manner (see Figure 6b right, single slice from one aggregate). Using nuclear markers like Sox2, our technique allowed quantification at single-cell resolution by identifying, counting, and summing the brightness of Sox2+ nuclei for each imaged confocal plane (see Figure 6c/d/f/h). Filamentous, abundant signals like Map2 could be quantified throughout 3D aggregates by summing the overall mean brightness for each confocal plane (see Figure 6e/g). The comparison of three 96-well plates from independent batches revealed the uniform cellular composition of AMOs within and between independently thawed and cultured batches (see Figure 6d–g) (Average CVSox2 = 5%, CVMap2 = 9%).
Positional analysis detected effects of plate position (edge effects) for Map2 levels but not Sox2 levels with about 10% reduced Map2 brightness of samples in the center of the plate (Figure 6—figure supplement 1) compared with the wells at the edge. Considered together with the absence of edge effects in the RNA sequencing results, this may indicate that only a specific subset of proteins is altered by edge conditions, while the vast majority of cellular processes is uniform throughout the plate.
smNPC-derived AMOs are morphologically, structurally, and functionally more homogeneous than automated hiPSC-derived organoids
Since differentiation outcomes and kinetics are known to vary considerably between cell lines, we decided to benchmark our AMOs against a protocol that can be implemented using the same starting cell line as our AMOs and that is adaptable to the same automation and analysis workflow that we established for our midbrain model. In this comparison, smNPC-derived AMOs (line 2) and the hiPSC organoids share the same cell line of origin. We compared our smNPC-derived cultures to hiPSC-derived 3D neural organoids based on a core protocol by Paşca et al., 2015 (also described in more detail by Sloan et al., 2018) with modifications, as they share a number of key traits with our AMOs (for an overview of our automated protocol and the differences to the published original, see Figure 7—figure supplement 1). They are self-aggregated and self-organized, and they do not require the addition of an external matrix for proper development. Furthermore, to eliminate any potential bias due to manual handling, we adapted the cortical protocol to our automation pipeline. As a result, any remaining variability did not originate from handling but from stochastic biological processes. For detailed characterization, hiPSC cortical organoids underwent our established automated whole mount staining and clearing procedure (see Figure 7—figure supplement 2) as well as qPCR (see Figure 3—figure supplement 1), confirming differentiation toward their correct cortical fate in our workflow. Out of twelve full 96-well plates we were, with our workflow, able to generate and maintain cortical organoids in all but one well, where the organoid got lost during the 30 days of automated culture. This further underlines the adaptability and efficiency of our midbrain protocol for other organoid types.
Compared to smNPC-derived AMOs, morphology and zonal arrangement of neural subpopulations in hiPSC organoids varied to a larger degree (see Figure 7a–d). Three independently cultured batches of automated hiPSC organoids showed up to a 5-fold higher coefficient of variation in cellular viability and up to a 10-fold higher coefficient of variation for organoid size than AMOs (see Figure 7g/h, for individual organoid size and viability data see Figure 7e/f).
We performed high-content analyses at the protein level analogous to the data in Figure 6 and found that the variation of Sox2 and Map2 content of the automated hiPSC-organoids was larger than for AMOs, even when we normalized for strongly variable sizes (see Figure 7i/j). The acquisition of high-content data at the same hardware settings for both types of organoids (from the same cell line of origin) also allowed for a direct comparison of cell-type-specific signals. hiPSC-organoids contained a distinctly lower amount of Map2 per area, indicating less efficient/delayed neuronal maturation.
In conclusion, AMOs are more homogeneous with regard to their morphology, cellular structure, size, viability, and protein expression than the automated hiPSC organoids generated from the same cell line and under strictly standardized conditions.
Automated midbrain organoids possess functional characteristics of midbrain tissue and allow assessment of neural subpopulations for high-throughput screening
A clinically relevant midbrain model requires dopaminergic activity. To further confirm midbrain-specific function of AMOs, we generated electrophysiological data from multielectrode arrays together with specific molecular agonists and antagonists for dopaminergic, GABAergic, and glutamatergic pathways. Functional responses were consistent with a midbrain identity of AMOs (see Figure 8e/f/g/k). Electric field oscillations of AMOs responded strongly to the addition of modulators of the dopamine pathway, but yielded more limited responses when targeting the GABA pathway (see Figure 8e/f/k). This correlates well with the lower abundance of GABAergic neurons in AMOs as detected by immunostaining (see Figure 2—figure supplement 1 for a vGAT staining). The electrical activity increases following the stimulation with glutamate agonists and subsequently decreases after addition of the antagonist (Figure 8g/k), despite the low abundance of glutamatergic cells (see vGLUT1 staining in Figure 2—figure supplement 1). However, this is consistent with previous studies reporting an increase in firing rate and burst-firing for dopaminergic neurons treated with glutamate agonists which can be counteracted by NMDA receptor antagonists (Meltzer et al., 1997; Wang and French, 1993). Automated hiPSC-derived cortical organoids served as negative controls to validate midbrain-specific MEA data. As expected following cortical differentiation, control hiPSC organoids had a weak response to dopaminergic modulation and reacted most strongly to glutamatergic modulation (see Figure 8h–j,k). The overall higher electric field amplitude in AMOs compared to automated hiPSC-organoids further illustrated the accelerated differentiation and faster maturation achieved by our workflow over hiPSC-based organoid protocols of the same age.
To further evaluate cellular maturity as one key factor in physiologically relevant screens, we determined the percentage of Sox2-negative mature cell types within a large number of AMOs via high-content imaging. Testing AMOs from two separate lines with three independently cultured batches each, we detected between 70% and 80% mature cell types after only 30 days of differentiation with very little variation within and between batches and cell lines (see Figure 8d).
The maturity and dopaminergic midbrain identity of AMOs was further confirmed by the spontaneous and unstimulated secretion of dopamine into the culture medium at a similar level as found in adult human cerebrospinal fluid (CSF, see Figure 8c; Goldstein et al., 2012). Sequestered dopamine may constitute an additional readout that is functionally relevant in models of the human midbrain.
To further assess the physiological relevance of our workflow as well as its ability to quantify drug effects in 3D cell cultures, we treated AMOs from two different cell lines with increasing concentrations of toxins specific for dopaminergic neurons, namely 6-Hydroxy Dopamine (6OHD) and 1-Methyl-4-phenylpyridinium (MPP+) (Reinhardt et al., 2013b; Blum et al., 2001; Simola et al., 2007; Meredith and Rademacher, 2011). High-content imaging allowed us to discern distinct effects of these nigral toxins on TH-positive neurons compared to general MAP2-positive neurons in AMOs. While both compounds reduced the number of all neurons in a dose-dependent manner, levels of TH-positive dopaminergic neurons decreased more strongly and at lower concentrations (see Figure 8a/b), with 6OHD having the more specific effect. Overall, our HTS workflow detected subpopulation-specific effects within dense, whole-mount-stained organoids with very little variation between replicates and between AMOs from separate cell lines.
We next evaluated the use of cleaved caspase 3 (cCasp3) as a general marker for cell death in AMOs for HTS toxicity studies (see Figure 8—figure supplement 1). If used judiciously, cCasp3 staining can yield quantitative dose-response curves, for example for the addition of known cell toxicants such as G418 (see Figure 8—figure supplement 1a, b, d). When used in co-staining scenarios, this makes cCasp3 potentially attractive to quantify the viability of a wide range of cellular subpopulations, such as Sox2+ precursor cells (see Figure 8—figure supplement 1c, e).
However, cCasp3 only labels cells that are currently undergoing apoptosis, and this signal disappears when dead cells are cleared from the tissue. Consequently, we observed that the cCasp3 signal dropped at high concentrations of toxicants (see Figure 8—figure supplement 1g), necessitating a careful optimization of assay timing and assay windows. This makes cCasp3 challenging to use for primary screens, where compounds with a wide range of toxicities and kinetics may need to be evaluated in parallel. Other viability assays, such as the ubiquitous CellTiter-Glo 3D assay, provide a much broader dynamic range and are less sensitive to timing, as they quantify living cells instead of a transient signal of dying cells (see Figure 8—figure supplement 1i). Unfortunately, CellTiter-Glo cannot distinguish survival of cellular subpopulations, since it relies on a non-discriminatory lysis protocol. Thus, CellTiter-Glo may be most suitable for primary screens and cCasp3 staining as a follow up to probe effects on specific cellular subpopulations. Another common readout in 3D cell culture screens is aggregate size (Yao et al., 2020; Mittler et al., 2017; Kim et al., 2020). Interestingly, the largest cross-sectional area of organoids did not correlate with G418 doses or CellTiter-Glo survival data in the conditions and time scales tested here, but remained largely constant despite escalating toxin levels (see Figure 8—figure supplement 1h). Our automated workflow is compatible with all three types of assay and the results demonstrate that AMOs show lower variance in cell death/survival in toxicity studies than iPSC-organoids derived from the same cell line (see Figure 8—figure supplement 1g, i).
Overall, this data highlights the suitability of AMOs as an ideal tool for midbrain-specific drug and toxicity screening with various high-throughput-compatible and functional readouts that can be directly applied to large-scale screening campaigns as well as more detailed secondary analyses to gain mechanistic insight into the behaviors of distinct cellular subpopulations in the context of their niche.
Discussion
In this report, we present a fully scalable, automated workflow for the generation, maintenance, fixation, immunostaining, clearing, and optical high-content analysis of human midbrain organoids. In designing our workflow, we used a standard 96-well format along with optimized protocols providing full compatibility with standard liquid handlers in most screening facilities. By omitting cumbersome matrix embedding steps, we eradicate variation originating from positional differences of cells within extracellular matrix (ECM) droplets during embedding. Importantly, our workflow is also compatible with manual pipetting, and thus can benefit labs without access to liquid handlers. By starting with neural precursor cells rather than PSCs, we can accelerate and streamline the neural differentiation process resulting in highly reproducible, homogeneous, functionally coupled, and electrophysiologically active human 3D microtissues. In contrast to the highly organized layers of the cortex, which can be partially recapitulated in cerebral/cortical organoids (Lancaster et al., 2013; Qian et al., 2016; Mariani et al., 2015; Bhaduri et al., 2020), different tissues of the midbrain have a much less distinctive morphological organization at the cellular scale. While our AMOs do not represent the complexity of the actual human midbrain, we demonstrate the presence of key cell populations with a degree of structure and organization similar to that of other published midbrain organoids (Qian et al., 2016; Jo et al., 2016; Monzel et al., 2017). AMOs’ uniformity in key parameters such as size, cellular organization, gene expression, and protein levels enables large-scale 3D-based HTS strategies. Moreover, the presence of aggregate-wide, reproducible, synchronous calcium and electrical activity in human neural 3D structures may enable neural activity to be used as a HTS-compatible, simple readout for phenotypic screens including toxicology studies and drug screening for disorders with altered brain activity.
While other groups have previously reported the use of human neural precursor cell (hNPC)-derived neurospheres for toxicity testing (Moors et al., 2009), as well as a hiPSC-derived high-throughput-compatible spheroid model including the use of calcium oscillations as a readout to evaluate neurotoxicity (Sirenko et al., 2019), these aggregates do not display structural tissue self-organization to the same extent as our organoids. Importantly, none of the available high-throughput-compatible protocols generates midbrain-specific organoids, but rather they focus on a cortical fate (Sirenko et al., 2019). Our model of the human midbrain opens up the potential to perform 3D organoid-based HTS of midbrain-specific disorders including the highly prevalent Parkinson’s disease. One of the most widely used models for Parkinson’s disease is the ablation of dopaminergic neurons with toxins such as 6-hydroxy dopamine or MPP+ (Reinhardt et al., 2013b; Blum et al., 2001; Simola et al., 2007; Meredith and Rademacher, 2011). Using these compounds as a dose-response model, our workflow demonstrates the unbiased, quantitative, cell-type-specific assessment of dopaminergic neurons in a physiologically relevant and complex human 3D model system with full scalability and high-throughput compatibility.
Our automated whole mount optical analysis can also be directly applied to other established 3D protocols and may help uncover phenotypes that manifest themselves primarily in 3D or in distinct subpopulations of cells. The neural precursor-based proof-of-principle demonstrated here may be suitable, for example, in screens for diseases affecting primarily neural precursors including Zika virus infection. While this condition has been intensively studied using organoids, platform technologies for 3D high-throughput screens are still lacking and in high demand (Qian et al., 2017).
Furthermore, other non-midbrain organoid protocols may benefit from our workflow; for example, we have adapted cortical organoids by Paşca et al., 2015/Sloan et al., 2018 to be scaled-up and analyzed in our system as automated hiPSC-based control organoids for our midbrain model. Automation of this existing manual protocol necessitated some modifications (summarized in Figure 7—figure supplement 1). For example, in order to create a workflow that is as scalable and as repeatable as possible, we decided to use accutase-digested single cell suspensions originating from feeder-free iPSC-cultures rather than dispase-treated colonies lifted off partially or whole from feeder-based iPSC-cultures. The use of colonies as a starting point makes colony/aggregate sizes difficult to control, especially in an automated setting, resulting in a variable starting point for each organoid generated in this fashion. We also aggregated and maintained the resulting organoids in SBS-compatible U-bottom plates that facilitate automated handling of all organoids as separate biological replicates, each in its own well, rather than in 100 mm dishes as in the original publication. These modifications to the original protocol may have changed the outcome compared to the published original; however, we demonstrate proper forebrain differentiation at the RNA and protein levels, and provide detailed single organoid-based high-content data from three independent batches for select parameters. We also believe that these modifications provide the most rigorous comparison between iPSC- and precursor-based organoid approaches by starting both protocols with well-defined and easily controllable cell numbers as single cell suspensions. By utilizing colonies as starting point, we would have biased reproducibility against the automated iPSC-protocol, thus weakening our comparison.
Despite our efforts to standardize iPSC-based organoids in a modified automated workflow, AMOs displayed a distinctly lower variance in a broad set of parameters, including size, cellular subpopulations, and survival in toxicity studies. This may indicate that mechanical standardization alone cannot compensate for the innate variability of PSCs during their self-organization.
Although the use of precursors cells as a starting population provides benefits, the comparison to hiPSC-based organoids also highlights certain limitations of this approach. Our choice of generating AMOs by seeding more committed neural precursor cells instead of hiPSCs improves the predictability of the differentiation outcome while maintaining the ability to form self-organized tissues. However, it also sacrifices some of the hiPSCs’ broader cell fate potential and complexity resulting in simpler structures compared to hiPSC-based organoids. In addition, smNPCs can only give rise to cells of the ectoderm, which excludes, for example, the formation of microglia. Although smNPCs are capable of efficiently generating oligodendrocytes (Reinhardt et al., 2013a; Ehrlich et al., 2017), they generally arise very late in neural development (Goldman and Kuypers, 2015), and other organoid protocols also demonstrate the presence of oligodendrocytes after 100 days of differentiation or more (Marton et al., 2019). These long time scales are not very attractive for the design of high-throughput screening strategies. Instead, we favor approaches that ectopically add a known number of oligodendrocytes or microglia to our AMOs and thus reach a screenable state much earlier. Similar strategies have recently been described for the introduction of microglia into hiPSC-based neural organoids (Lin et al., 2018; Abud et al., 2017; Muffat et al., 2018), albeit without the ability to perform HTS. AMOs fill the gap between heterogeneous PSC-derived organoids and established but less physiological spheroid- and 2D-based screening formats.
The emerging science of 3D cell culture promises to probe the effects of drugs on single cell types (here: dopaminergic neurons) in an intact niche as part of a complex 3D in vitro culture. Understanding the effects in context will be an essential factor to better understand biology at the tissue level. Currently, 3D cell culture science is in a phase of transition, where traditional manual low throughput protocols need to be adapted to unbiased higher throughput workflows. This is a prerequisite to mine the information necessary to merge organoid science with the promises of big data and machine learning to better tackle the complexities of understanding 3D biology. Therefore, we see our workflow as a contribution to help others translate successful strategies such as cell painting (Bray et al., 2016) and other big-data-generating high-content strategies (Friese et al., 2019; Scheeder et al., 2018) into biologically relevant 3D models. If successful, this will help to open up single-cell-based phenotypic discovery to the third dimension by taking 3D-based HTS approaches beyond bulk techniques such as cell survival or gross morphology. Taken together, we hope that our automation approach can contribute to establishing a next generation of cellular 3D in vitro disease models that allow unbiased, quantitative, high-throughput access to human tissue-surrogates in a dish.
Materials and methods
smNPC culture
View detailed protocolAll cell lines used in this study tested negative for mycoplasma contamination in PCR- and sequencing-based analyses. Unless otherwise noted, all cells and 3D aggregates were maintained at 37°C and 5% CO2. The human small molecule precursor cells (smNPCs) were generated and characterized during a previous study (Reinhardt et al., 2013a) and cultured as described there, with minor modifications. We grew smNPCs in 1% (v/v) Matrigel (BD)-coated 6-well plates (Sarstedt) in N2B27 medium supplemented with the small molecules smoothened agonist (SAG) (0.5 μM, Cayman Chemical) and CHIR 99021 (3 μM, Axon MedChem). N2B27 consisted of DMEM-F12 (Thermo Fisher) and Neurobasal Medium (Thermo Fisher) at a 1:1 ratio, enriched with 1:400 diluted N2 supplement (Thermo Fisher), and 1:200 diluted B27 supplement without vitamin A (Thermo Fisher), 1% penicillin/streptomycin/glutamine (Thermo Fisher), and 200 μM ascorbic acid (Sigma-Aldrich). Typically, we exchanged medium every other day. The cells were split every 5–7 days at a splitting ratio of 1:10 to 1:20 via accutase (Sigma-Aldrich) treatment for approximately 15 min at 37°C, yielding a single-cell solution. To stop the digestion, the cells were diluted in DMEM-F12 with 0.1% bovine serum albumin (BSA) (Thermo Fisher) and centrifuged at 220 g for 2 min. The cell pellet was resuspended in fresh smNPC medium (N2B27 with 0.5 μM SAG and3 μM CHIR) and plated on Matrigel-coated 6-well plates.
Definition of batch
Request a detailed protocolEach batch consisted of independently frozen and thawed cells of the same passage and cell line, passaged separately and seeded, maintained, fixed, stained, and cleared separately from the other batches.
AMO generation
View detailed protocolAll liquid handling steps (seeding, maintenance, and fixation of organoids) were performed by a Beckman Coulter Biomek FXP liquid handling station equipped with a 96-channel-pipetting head and an attached Cytomat incubator (Thermo Fisher). After digestion by accutase, we seeded 9000 smNPCs in each well of a conical 96-well plate (Thermo Fisher) in smNPC medium and allowed them to aggregate for 2 days. To increase inter-cell adhesion, we added 0.4% (w/v) polyvinyl alcohol (PVA, Sigma-Aldrich). Starting at day 2, cells undergo ventral patterning over 4 days in two feedings by removal of CHIR99021 in the continued presence of 1 μM SAG and the addition of 1 ng/mL brain derived neurotrophic factor (BDNF, PeproTech) and 1 ng/mL glial cell line-derived neurotrophic factor (GDNF, PeproTech). After ventralization, we removed SAG on day 6, further supported midbrain differentiation and maturation by the addition of 1 ng/mL transforming growth factor beta 3 (TGFβ−3, PeproTech), and 100 μM dibutyryl cyclic adenosine monophosphate (dbcAMP, Sigma-Aldrich). To boost maturation and cell survival during the rest of the neural maturation, we increased the concentration of BDNF and GDNF to 2 ng/mL each starting at day 6. A single dose of 5 ng/mL Activin A (eBioscience) was added on day 6 only. Depending on the desired degree of maturity, the duration of the maturation phase can be prolonged to 100 days and longer. Organoids were fed every other day for the entire duration of culture.
Size measurement of organoids
Request a detailed protocolFor size measurements, we took brightfield images of randomly selected AMOs or hiPSC organoids using a stereo microscope (Leica MZ10 F, camera: Leica DFC425 C). Images were processed with ImageJ/Fiji (Schindelin et al., 2012) using a custom-tailored standardized workflow. The auto threshold function was used to discriminate samples from the background followed by a measurement of their area with the analyze particles function. The measured area corresponds to the largest cross-section of the organoid. Data were outputted to Microsoft Excel and GraphPad Prism v8.4.2 (Graphpad Software, Inc) for further analysis. Coefficients of variation (CVs) were calculated via CV = standard deviation/mean.
Whole mount staining and clearing
View detailed protocolIn order to analyze protein expression in 3D in a HTS-compatible manner, we adapted a whole mount staining protocol based on Lee et al., 2016. for large-scale 3D aggregates and optimized it for use in an automated liquid handling system (Beckman Coulter Biomek FXP, 96-channel-pipetting head). After fixation with 4% PFA (VWR) for 10–15 min, we stained the aggregates with primary and secondary antibodies (Alexa Fluor secondary antibodies, Thermo Fisher) for 6 days each. A list of all primary antibodies and concentrations can be found in Supplementary file 3. We diluted the antibodies in a blocking and permeabilization solution (6% BSA, 0.5% Triton-X 100 (Roth), 0.1% (w/v) sodium azide (Sigma-Aldrich) in PBS (Sigma-Aldrich)) and renewed it every 2 days. Between primary and secondary antibody incubation as well as after the staining procedure we washed the samples 5 times for 1 hr with 0.1% Triton X-100 in PBS. This extremely long staining procedure allows the antibodies to fully penetrate the aggregates despite their large size and high density. To enable full penetration by microscope illumination, the whole mount staining procedure is followed by BABB-based tissue clearing (Dent et al., 1989). First, the samples were dehydrated stepwise through a methanol (Roth) series (25%, 50%, 70%, 90%, 100%, 15 min each). Next, they were transferred to an organic solvent-resistant cyclo-olefin 96-well plate (‘Screenstar’, Greiner Bio-One). The samples were incubated for 30 min in 1:1 methanol/BABB (benzyl benzoate (Sigma-Aldrich) and benzyl alcohol (Sigma-Aldrich) 1:1) and subsequently kept in BABB for imaging. We used Imaris v9.1.2 (Bitplane, Oxford Instruments) for 3D rendering of confocal slices to produce Video 1.
Quantitative real-time PCR
Request a detailed protocolWe performed RNA isolation for quantitative real-time PCR (qPCR) analysis using the NucleoSpin RNA XS kit (Macherey-Nagel) according to the manufacturer’s instructions. Depending on the age and protocol used, we pooled organoids from one batch in order to yield enough RNA for downstream analysis; AMOs: 32 (d6), 24 (d16), or 18 (d22 and d30), hiPSC organoids: 32 (d6), 14 (d22), or 10 (d30). We determined RNA concentration and purity using a NanoDrop 8000 spectrophotometer (Thermo Fisher) and performed reverse transcription according to standard protocols using 1000 ng RNA per reaction. qPCR was done in triplicates on a Quantstudio 5 Real-Time PCR System (Applied Biosystems) with iTaq Universal SYBR Green Supermix (Bio-Rad) and 3.2 ng RNA equivalents per 10 μL reaction. Cycling conditions were 2 min at 95°C followed by 40 cycles of 15 s at 95°C and 60 s at 60°C. We calculated the relative expression using the ΔΔct method and normalized to undifferentiated smNPCs (AMO line 1) collected before aggregation as well as using GAPDH as a housekeeping gene. Alternatively, gene expression was quantified with the Biomark 48.48 integrated fluidic circuit (IFC) Delta Gene assay (Fluidigm) according to the manufacturer’s instructions. Briefly, following 14 cycles of preamplification, the samples were subjected to an exonuclease I (New England Biolabs) treatment (37°C for 30 min and 80°C for 15 min) and diluted twentyfold with DNA Suspension buffer (TEKnova). The samples (in duplicates) and assay mixtures were loaded onto a 48.48 microfluidic IFC chip and run on the BioMark real-time PCR reader (Fluidigm) where they were amplified and measured according to manufacturer’s instructions. Here, data analysis was performed using the BioMark real-time PCR analysis software 4.3.1 (Fluidigm) with standard settings. Again, undifferentiated smNPCs (AMO line 1) were used as a reference and GAPDH served as housekeeping gene. All data was transferred to Microsoft Excel for further processing and GraphPad Prism v8.4.2 for plotting. A list of all used primers can be found in Supplementary file 4.
Calcium imaging
Request a detailed protocolFor calcium imaging, we added 10 μM cell-permeant Fluo-4 AM (Thermo Fisher) diluted in AMO medium to the aggregates and incubated for 60 min at 37°C. For inhibitor studies cobalt(II) chloride (Sigma-Aldrich) was added to the medium at a concentration of 2 mM together with the Fluo-4 AM. Imaging was performed using a Dragonfly spinning disc confocal microscope (Andor, Oxford Instruments) at a frequency of 10 Hz for 4 min. Data analysis was performed using ImageJ/Fiji (Schindelin et al., 2012). First, different ROIs were defined as depicted in Figure 4. Then, the mean fluorescence intensity in those ROIs was measured over time and plotted using GraphPad Prism v8.0.2. The video was assembled via ImageJ/Fiji (Schindelin et al., 2012) and the frame rate accelerated to compress 4 min real time at 10 Hz into 20 s running time. Alternatively, we measured fluorescence intensity on a Synergy Mx plate reader (BioTek), acquired data with the Gen5 software (BioTek) and outputted it to Microsoft Excel and GraphPad Prism v8.0.2 for further analysis and plotting.
Electrophysiological analysis of AMOs by microelectrode array
Request a detailed protocolElectrophysiological measurements on microelectrode arrays (USB-MEA256system, Multichannel Systems) were performed on electrode areas of 9-well MEAs as previously described (Piccini et al., 2017). The MEAs were plasma-cleaned and coated with 1:75 diluted Matrigel (Corning) in KO-DMEM (Invitrogen) overnight and additionally for 2 hr before seeding with a 0.1% gelatin (Sigma-Aldrich) in PBS solution at room temperature. The MEAs were pre-warmed to 37°C, organoids were transferred to the electrode area of the MEAs, and allowed to attach for approximately 28 hr. Subsequently, MEA chambers with attached samples were used for electrophysiological recordings at 37°C. To study the effects of different pharmacological modulators, organoids were first measured under basal conditions (i.e. no addition of compounds) to record a reference signal. We then added specific pathway activators to each sample chamber of a MEA and recorded the electric field potential after a brief period of equilibration. After recording the signal for the activators, we next added the inhibitors for the respective pathways to the same well and repeated the measurement procedure as before. This guaranteed that we measured the change of electric field potential of each sample and the basal activity of each sample could serve as an internal control. The compounds used were dopamine hydrochloride (10 μm, Sigma-Aldrich), risperidone (10 nM, Sigma-Aldrich), GABA (10 μm, Sigma-Aldrich), bicuculline (1 μm, Sigma-Aldrich), glutamate/glycine (100 μm each, Sigma-Aldrich), and ketamine (10 μm, Sigma-Aldrich). It is possible to culture organoids on MEAs for longer and to achieve a more widespread contact with MEA electrodes. However, organoids flatten and change morphology and possibly cellular composition upon prolonged attachment, and we sought to measure the organoids’ electrical activity as close to their usual spherical state as possible. 28 hr allowed sufficient time for attachment of the aggregates to record electric fields from still spherical samples while providing enough mechanical cell-substrate connection to transport MEA substrates from the incubator to the recording rig. After the recordings, it was possible to gently remove the organoids without causing observable damage. Although not optimized (as we only performed measurements at one time point per organoid), this opens up the possibility to transfer the organoids back to standard culture conditions after MEA measurements and thus perform multiple measurements over time. Datasets were recorded with Cardio2D software (Multichannel Systems). Analyses were performed using the software Cadio2D+ (Multichannel Systems) and Origin v9.0 (OriginLab) on exported data. Discrete fast Fourier analyses in Origin (Blackman window) was used to assess frequencies of autonomous activity of the organoids. To compare neural sample activity via electric field oscillations (Figure 8k), we pre-processed the raw MEA data with a Savitzky-Golay filter (curves were smoothed with a window size of 50 to preserve peak data while removing noise) in Origin v9.0 and formed the total sum of the absolute values of each oscillation over 15 s. The results were outputted to Microsoft Excel, reformatted and then transferred to GraphPad Prism v8.4.2 for plotting.
Electrophysiological analysis of single cells by patch-clamping
Request a detailed protocolDue to the morphology of AMOs (high optical density and the fact that most cell bodies are located in a depth of at least 10–20 µm), it was technically impossible to perform the patch-clamp measurements on intact aggregates. Therefore, the organoids were treated with 1 mg/ml trypsin and then mechanically dispersed to obtain single cells. These were seeded on PDL-coated coverslips and cultured for 1–3 days in AMO medium (we stated the age of AMOs at the time of dissociation). The transmembrane currents were recorded from isolated cells using the whole-cell configuration of the patch-clamp technique (Hamill et al., 1981). The patch pipettes were fabricated from borosilicate glass on a Sutter P1000 (Sutter Instrument company) pipette puller. When filled with pipette solution, they had a tip resistance of 4–6 MΩ. Recordings were done using an EPC-10 amplifier (HEKA Elektronik) and Patchmaster acquisition software (HEKA Elektronik). Series resistance, liquid junction potential, pipette and whole-cell capacitance were cancelled electronically. Bath solution contained (mM): NaCl 140, KCl 2.4, MgCl2 1.2, CaCl2 2.5, HEPES 10, D-glucose 10, pH 7.4 and the pipette solution contained (mM): K-aspartate 125, NaCl 10, EGTA 1, MgATP 4, HEPES 10, D-glucose 10, pH 7.4 (KOH). We performed all experiments at room temperature. Recordings of current-voltage relationship (I-V curves) were done in voltage-clamp mode at a holding potential of −70 mV. Recordings of evoked action potentials were performed in current-clamp mode. Data were analyzed using Patcher’s Power Tool routine for IgorPro (WaveMetrics), SciDAVis (http://scidavis.sourceforge.net/) and Origin Pro 2019 (Origin Lab). To reveal the shape of I-V curves, single traces were normalized to the peak amplitude and then averaged.
RNA sequencing
Request a detailed protocolTo isolate RNA of single AMOs and organoids we used the Direct-zol-96 RNA kit (Zymo Research) according to the manufacturer’s instructions. We assessed RNA concentration and purity using a NanoDrop 8000 spectrophotometer and RNA integrity with a Bioanalyzer (Agilent Technologies) per standard protocols. Next, mRNA was enriched using the NEBNext Poly(A) Magnetic Isolation Module (NEB) followed by strand-specific cDNA NGS library preparation (NEBNext Ultra II Directional RNA Library Prep Kit for Illumina, NEB). The size of the resulting library was controlled by use of a D1000 ScreenTape (Agilent 2200 TapeStation) und quantified using the NEBNext Library Quant Kit for Illumina (NEB). Equimolar pooled libraries were sequenced in a single read mode (75 cycles) on the NextSeq 500 System (Illumina) using v2 chemistry yielding in an average QScore distribution of 95% >= Q30 score and subsequent demultiplexed and converted to FASTQ files by means of bcl2fastq v2.20 Conversion software (Illumina).
RNA sequencing analysis
Request a detailed protocolWe aligned the RNA sequencing reads to the human genome hg19 with TopHat2 aligner (v2.1.1) (Kim et al., 2013), using default input parameters. Gene annotation from Ensembl (version GRCh37.87) were used in the mapping process. The number of reads that were mapped to each gene was counted using the Python package HTSeq (v0.7.2) (Anders et al., 2015) with ‘htseq-count – mode union – stranded no’. For the correlation analysis, sequencing data of 21 different human fetal organs and midbrain(-like) samples were obtained from GSE66302 and E-MTAB-4868, respectively. Human prenatal cortex (at 24 weeks post-conception) RNA-seq datasets were obtained from www.nature.com/neuro/journal/v18/n1/extref/nn.3898-S9.zip. Reads were mapped to the human genome as described above. RPKM values (Reads Per Kilobase of transcript per Million mapped reads) for each gene were computed by Cufflinks (v2.2.1). We selected genes with high expression (log RPKM > 1) for further analysis. Based on the expression of the selected genes, Pearson correlations were calculated. We averaged the correlation coefficients for biological replicates. PCA and differential expression analysis were performed with raw counts using the R package DESeq2 (v1.18.1). Genes were considered as deregulated if |log2FC| > 2 and FDR < 0.05 using Benjamini-Hochberge multiple test adjustment (Benjamini and Hochberg, 1995). Gene Ontology (GO) term enrichment was analyzed with the bioinformatics web server Gorilla (Eden et al., 2009) and visualized with REViGO (Supek et al., 2011).
For the comparison with cerebral organoids (Figure 5—figure supplement 1), the dispersion within groups (Figure 5—figure supplement 1b) was calculated using the average distance between data points and centroids on the PCA plot (Figure 5—figure supplement 1a).
All RNA sequencing data generated by us was deposited to the NCBI GEO database (GSE119060) and can be accessed at https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE119060.
Quantification of whole mount staining and clearing
View detailed protocolTo assess how quantitative our imaging workflow is, we performed a dilution experiment. We mixed unlabeled smNPCs with different percentages (1.25%, 2.5%, 5%, 10%, 20%, 40%) of CellTracker deep red dye (Life Technologies)-labeled cells (labeling according to standard protocols, dye concentration 1:20,000) and aggregated them in smNPC maintenance medium with 0.4% PVA. To explore the effects of overall aggregate size on quantitation, we generated aggregates with 100,000 as well as 200,000 cells in total. After 1 day of aggregation, the aggregates were fixed with 4% PFA, subjected to BABB-based tissue clearing, imaged, and analyzed as described below.
General workflow for high-content imaging and analysis
After staining and clearing, we achieved uniform aggregate positioning within the wells by tilting the plates off the horizontal plane at 60 degrees for 1 min. Image acquisition was carried out in an Operetta high content imager (Perkin Elmer) and images were analyzed in Harmony 4.1 software and/or Columbus version 2.6.0.
Each high-content imaging experiment requires a customized set of parameters adjusted for the size of the aggregates, marker-wavelengths, their morphology and distribution within the sample, marker combinations, and signal and background intensities. As the details of each workflow are particular to our hardware and software setup, we first provide a general description of our workflow that focuses on explaining the principal steps in a manner that allows platform-independent reproduction. We then provide the detailed analysis pipeline for the data shown in Figure 6 (quantification of Sox2 and Map2 in AMO line 1) with all its unique parameters and steps in the Materials and methods section titled ‘Detailed high-content imaging analysis example for the data shown in Figure 6'.
Acquisition
View detailed protocolAs our version of high content imager did not include automatic sample-finding capability, we started by acquiring well overview images at 2x magnification. Although uniform, locations of aggregates were not identical from well to well, necessitating acquisition of tiled 3D stacks of adjacent fields of view (FOV) with 10x magnification covering all possible locations of aggregates within a small area of given wells. The 2x overview images guided the choice of location for the 10x FOVs. This resulted in tiled non-overlapping stacks of images, some of which contained aggregates, and some of which did not. These 10x images entered further image analysis.
Analysis
In general, after hardware-based background correction, image analysis was performed separately for each FOV and confocal plane, and results were either reported as mean or sum for all FOVs and planes per well as indicated in the figures and figure legends. Further analysis consisted of three principal stages: (1) sample identification, (2) background correction, (3) segmentation, and (4) data reporting.
To identify aggregates and distinguish them from occasional contaminations, we summed all acquired channels (or in case of high background noise in particular channels used a single one with homogeneous signal) and smoothed the resulting image with either a gaussian or a median kernel to remove local outlier pixels. Simple thresholding identified organoids as regions of interest (ROI), which were further refined via morphology (e.g. area and roundness) and intensity to exclude incorrectly identified objects such as dust or fibers.
As large spherical aggregates possessed a non-uniform background signal, we performed local background correction before image segmentation. Subtracting a blurred version of the image from the raw image removed most of the non-uniform haze in each plane. To boost fine structures, sliding-parabola algorithms provided additional local contrast when needed.
Structures of interest within each bona fide aggregate identified in (1), e.g. nuclei or cells, were segmented by either simple thresholds or by dedicated proprietary algorithms (often called ‘find nuclei’ or ‘find cells’) and further refined by morphology and intensity if needed.
Finally, the parameters of interest (often intensity, area, number of objects) were measured from the raw image channel based on the previously defined and selected ROIs. If object segmentation was not possible (e.g. due to insufficient resolving power at 10x for highly filamentous markers or low signal to noise ratios), raw/total integrated signal intensity within the defined sample region was reported. In order to generate results for the entire 3D aggregate, which may span several FOVs, datasets for all planes and fields were combined. Often, it is also necessary to normalize the data to the aggregate area to account for sample variations across planes and wells.
Detailed high-content imaging analysis example for the data shown in Figure 6
View detailed protocolAfter staining and clearing, we achieved uniform aggregate positioning within the wells by tilting the plates off the horizontal plane at 60° for 1 min. Image acquisition was carried out in an Operetta high content imager (Perkin Elmer) and images were analyzed in Harmony 4.1 software. We acquired a total of 16 confocal planes in three channels (DAPI, Sox2 488 nm, and MAP2 647 nm) with an inter-plane spacing of 36.6 μm for a total stack of 549 μm, covering the entire aggregate height (clearing and dehydration steps tend to shrink and flatten aggregates slightly). To define the aggregate region on each image plane, all three channels were summed, filtered with a median filter to remove small localized features, and bright areas were identified via the 'find image region' function. After cleaning the edge of the aggregate region by dilation and erosion steps of 10 and 3 pixels, respectively, we identified bona fide AMOs by selecting for regions with a minimum of 300 arbitrary brightness units (abu) and 4000 μm2 size. In order to better isolate Sox2+ nuclei from the general background, we ran a sliding parabola algorithm with a curvature setting of 2 across each image plane in the 488 nm channel. Nuclei were then identified within each aggregate region via the ‘find nuclei’ function, algorithm 'M' and further selected to be Sox2+ if they were larger than 10 µm2 and brighter than 1200 abu. We excluded image artifacts, small dust particles, and overlapping nuclei by omitting nuclei brighter than 6000 abu and larger than 70 μm2 from further quantification. For final output, the number and total brightness of nuclei in 488 nm and of aggregate regions in 647 nm were summed for all planes and all fields of view in each well and transferred to Microsoft Excel and TIBCO Spotfire for further annotation, analysis, and plotting. We omitted data from wells that contained dust particles within the same FOV as organoids, incompletely imaged aggregates due to improper positioning, or AMOs that have been damaged or lost during culture or downstream processing. Plate 1, 2, and 3 represent independent batches of separately thawed and cultured cells of the same frozen batch.
Human-induced pluripotent stem cell culture
View detailed protocolFor cortical organoid generation (see method section below), human induced Pluripotent Stem Cells (hiPSCs) (parental line of AMO line 2 smNPCs) were generated and characterized during a previous study (Reinhardt et al., 2013b). Cells were cultured in Vitronectin (Thermo Fisher)-coated 6-well plates in mTeSR Plus medium (Stemcell technologies) supplemented with 1% penicillin/streptomycin. The medium was changed every other day and cells were passaged at a ratio of 1:10 to 1:15, using accutase when they reached 80–90% confluency. After splitting, the medium was supplemented with the ROCK inhibitor Y-27632 (10 μM, tebu-bio) until the first media exchange.
For the generation of iPSC-based cerebral organoids (see method section below and Figure 5—figure supplement 1) hiPSC culture was performed feeder-free using modified FTDA medium (Frank et al., 2012) in 1% (v/v) Matrigel-coated 6-well plates with iPSCs from the same line as AMO line 2. FTDA medium consisted of DMEM-F12 supplemented with 1% human serum albumin (Biological Industries), 1% Chemically Defined Lipid Concentrate (Life Technologies), 0.1% Insulin-Transferrin-Selenium (BD), 1% penicillin/streptomycin/glutamine. We fed the iPSCs daily and added 10 ng/mL FGF2 (PeproTech GmbH), 0.2 ng/mL TGFβ3 (PeproTech GmbH), 50 nM Dorsomorphin (Enzo Life Sciences), 5 ng/mL Activin A (eBioscience), 20 nM C59 (Tocris) before each media exchange. We split the iPSCs as single cells every 3–5 days using accutase for approximately 10 min at 37°C. We transferred 600,000 cells per well of a 6-well plate to be seeded to DMEM-F12 with 0.1% BSA and centrifuged at 220 g for 2 min. We resuspended the cell pellet in fresh FTDA medium supplemented with 1:2000 ROCK inhibitor Y-27632 (tebu-bio) and plated the iPSCs on Matrigel-coated 6-well plates.
Generation of automated cortical organoids
Request a detailed protocolAs a control for our AMOs, we generated cortical hiPSC organoids from the same hiPSC line used to derive smNPCs for AMO line 2 (Reinhardt et al., 2013b). After manual 2D culture of hiPCs, all steps were fully automated using our liquid handling system (Beckman Coulter) with attached incubator (Thermo Fisher). Generally, we followed the protocol previously published by Paşca et al., 2015. (and also described in more detail by Sloan et al., 2018), with adaptations for our automation pipeline (see Figure 7—figure supplement 1). Starting with 90–100% confluent cultures, we detached hiPSCs with accutase and seeded 10,000 cells per well in ultra-low attachment U-bottom plates (Corning). Cortical organoid medium consisted of DMEM F-12, 20% Knock-out Serum replacement (GIBCO), 1% penicillin/streptomycin/glutamine, 1% Non-essential amino acids (Sigma-Aldrich), and 0.2% 2-Mercaptoethanol (Thermo Fisher). For the first 6 days, we supplemented the cortical organoid medium with 5 μM dorsomorphin (Enzo Life Sciences) and 10 μM SB-431542 (Biomol). During seeding only, we also added 10 μM ROCK inhibitor Y-27632. Aggregates were fed every 3 days using an automated liquid handling system. From day 6 to 24, culture medium supplements were exchanged to EGF and FGF2 (both 20 ng/ml, PeproTech) and afterwards BDNF and NT3 (metabion) (both 20 ng/ml).
Organoid viability assay
View detailed protocolTo measure the viability of individual organoids we used the CellTiter-Glo 3D Cell Viability Assay (Promega) according to the manufacturer’s instructions. The entire procedure was performed using an automated liquid handling system (Beckman Coulter) and is thus fully scalable and HTS-compatible. In short, the reagent and the AMOs were brought to room temperature in their 96-well culture plates for 30 min and the media volume of each 96-well was adjusted to 55 µl. We added an equal volume (55 µl) of the CellTiter Glo 3D reagent and let it shake on a Thermomixer (Eppendorf) at 900 rpm for 5 min before incubating the samples protected from light at room temperature for 25 min. To prevent cross-talk between wells when measuring the luminescence, we next transferred the contents from the clear 96-well culture plates to opaque white 384-well Lumitrac plates (Greiner) with two technical replicates per sample. Luminescence was recorded immediately after transfer with a Synergy Mx plate reader (BioTek). The results were outputted to Microsoft Excel, reformatted and then transferred to GraphPad Prism v8.4.2 for plotting. Coefficients of variation (CVs) were calculated via CV = standard deviation/mean.
Dopaminergic neuron-specific toxicity testing and quantification
Request a detailed protocolOn day 47 of differentiation, AMOs received Tox medium (TM) overnight to remove antioxidants present in the B27 media supplement. TM consisted of DMEM-F12 supplemented with 1% N2 and 1% penicillin/streptomycin/glutamine. The next day, AMOs were treated with different concentrations (0, 50, 100, 250, 500 μm) of either 6-Hydroxydopamine hydrochloride (Sigma-Aldrich) or 1-Methyl-4-phenylpyridinium iodide (MPP+, Sigma-Aldrich) in TM. To refresh compounds, AMOs received TM including toxins at the same concentrations as before 24 hr after the first application. After a total incubation of 48 hr, the medium was changed back to standard AMO medium, and the samples were cultured for 6 more days to allow cell death to occur and the dead cells to be cleared from the organoids. After 6 days, the samples were fixed with 4% PFA, whole mount stained for Map2 and TH, BABB-cleared, and subjected to high-content confocal imaging as described above. We performed image analysis following the steps as outlined in the high-content analysis section with slight modifications to accommodate the individual brightness, morphology, and background characteristics of the staining. Briefly, after identifying AMOs, we ran a Gaussian smoothing algorithm across the TH 647 nm (10 px width) and Map2 488 nm (five px width) channels and subtracted the smoothed images from the raw images to better isolate TH+ and Map2+ cells from the background. In the TH 647 nm channel, cells were then identified with the ‘find cells’ function, algorithm ‘C’ and further selected to be TH+ if they were brighter than 200 abu and larger than 25 μm2. For Map2, an additional sliding parabola algorithm with a curvature of 10 was run across the subtracted image to further reduce background noise. Map2+ cells were then also identified using the ‘find cells’ function algorithm ‘C’. As final output, the total intensity of the identified TH+ cells in the raw 647 nm channel and the Map2+ cells in the raw 488 nm channel were summed for all fields of view and confocal planes per well. The data was then transferred to Microsoft Excel for further analysis, including normalization to the organoid area, and plotted using GraphPad Prism v8.4.2.
Dopamine secretion
Request a detailed protocolWe collected the cell culture supernatant from 35 days old organoids 40 hr after feeding and measured its dopamine content using a Dopamine ELISA Kit (Abnova) according to the manufacturer’s instructions. Measurements were performed in duplicates and the sample concentrations were calculated based on a standard curve.
Determining the proportion of mature cells
Request a detailed protocolSince there is no single marker that identifies all different mature cell populations within the AMOs, but all immature smNPCs express the precursor marker Sox2, which is downregulated upon maturation, we defined Sox2-negative cells as mature for the purpose of this analysis. Using confocal high-content imaging analysis as described above, we quantified the number of Sox2-positive cells as well as the total cell number (based on DAPI-stained nuclei) within the entire organoids. We then calculated the ratio of Sox2-negative cells as 1 - (number of Sox2-positive cells) / (total cell number). Calculations were performed in Microsoft Excel and data was transferred to GraphPad Prism v8.4.2 for plotting.
ScaleSQ clearing
Request a detailed protocolScalesSQ tissue clearing was performed as previously described by Hama et al., 2015. Briefly, we incubated organoids in ScaleSQ (22.5% (w/v)) D-sorbitol (Sigma) and 9.1 M urea (Sigma) for 2 hr at 37°C. Next, we exchanged the solution to ScaleS4 (40% (w/v)) D-sorbitol, 10% (w/v) glycerine (Roth), and 15–25% (v/v) DMSO (Sigma) also at 37°C. ScalseS4 was renewed after 2 hr and the organoids were maintained at 37°C until analysis. This is critical, as the clearing effect decreases visibly after 1 hr at room temperature.
X-Clarity clearing
Request a detailed protocolWe performed X-Clarity tissue clearing using an X-CLARITY Polymerization System (Logos biosystems) per manufacturer’s instructions. First, organoids were fixed with 4% PFA overnight at 4°C. The next day, the organoids were embedded in a hydrogel monomer solution and incubated 12–24 hr, followed by a 3-hr polymerization step at −90 bar and 37°C in the X-CLARITY Polymerization System. The embedded organoids were transferred to a tissue container and cleared in the electrophoretic tissue chamber for 2 hr at 1.2 A and 37°C.
ClearT clearing
Request a detailed protocolFor ClearT tissue clearing we followed a protocol by Kuwajima et al., 2013. In short, fixed organoids were incubated for 30 min at room temperature in increasing formamide (AppliChem) concentrations: 20%, 40%, 80%, and 95% (AppliChem). Samples were directly used for further applications or stored at 4°C.
CUBIC clearing
Request a detailed protocolCUBIC tissue clearing was modified after Susaki et al., 2015. Briefly, we incubated the organoids for 6 hr at room temperature in CUBIC-1 ½ (one part water and one part CUBIC-1: 25% (w/w) quadrol (Sigma), 25% (w/w) urea, and 15% (w/w) triton X-100 in water). Afterwards, samples were incubated for 24 hr in CUBIC-1 at 4°C and 24 hr in CUBIC-2 (25% (w/w) urea, 50% (w/w) sucrose, and 10% (w/w) triethanolamine (Sigma)).
DAPI detection depth for evaluation of clearing efficiency
Request a detailed protocolWe stained the organoids with 0.5 μg/mL DAPI (Sigma) in PBS for 24 hr and subsequently subjected them to the different tissue clearing protocols. Afterwards, a LSM 700 scanning confocal microscope (Zeiss) was used to acquire z-stacks of the stained and cleared organoids. Three XZ and three YZ cross-sections per aggregate (n = 10 aggregates per clearing method) were used to quantify the maximum depth at which the DAPI signal could still be detected at a given brightness threshold. The depth of each cross section was measured manually at n = 5 different positions for each slice and n = 10 organoids per clearing protocol using ImageJ/Fiji. The data was exported to Microsoft Excel for further processing and GraphPad Prism v7.0 for plotting and statistical analysis.
Light transmission analysis for evaluation of clearing efficiency
Request a detailed protocolFollowing clearing, we took brightfield images of the organoids (n = 10 per protocol, n = 6 for X-Clarity) using an upright stereomicroscope (Leica MZ10 F, camera: Leica DFC425 C) under standardized brightfield transmission conditions. The mean brightness of the organoid area, measured by ImageJ/Fiji, served as measure for the amount of light transmitted through the aggregate. The data was exported to Microsoft Excel for further processing and GraphPad Prism v7.0 for plotting and statistical analysis.
Statistical analysis for clearing quantification
Request a detailed protocolFor statistical analysis of the clearing protocol comparisons, we used GraphPad Prism v7.0 and performed unpaired, two-tailed t-tests with α = 0.05 as normal distribution and equal variances can be assumed for the analyzed data.
Electron microscopy
Request a detailed protocolFor electron microscopy analysis, we used n = 3 AMOs from one batch at day 32. Samples were initially fixed for 3 hr with 2% glutaraldehyde, 2% paraformaldehyde in 0.2 M cacodylate buffer, pH 7.2. Afterwards, the specimen was dissected into smaller pieces and post-fixed with 1% osmium tetroxide containing 1.5% potassium hexacyanoferrate. Samples were dehydrated stepwise, including an over-night uranyl-block staining step in 70% ethanol. The specimen was orientated and flat embedded in epon. In total, three samples sectioned under different angles were analyzed at the electron microscope to visualize different aspects of the AMOs.
Generation of iPSC-based cerebral organoids
Request a detailed protocolFor iPSC-derived cerebral organoid generation, we followed the protocol by Lancaster et al., 2013. with minor modifications. Briefly, we dissociated iPSCs to single cells by accutase treatment and plated 9000 cells per well in a conical 96-well plate in low FGF stem cell medium (DMEM-F12 with knockout serum replacement (KOSR, Thermo Fisher) 1:5, fetal bovine serum (Biochrom) 1:33.3, 1% penicillin/streptomycin/glutamine, 1% non-essential amino acids (NEAA, Sigma-Aldrich), 2-mercaptoethanol (Thermo Fisher) 1:143, 4 ng/μL FGF2, 50 μm ROCK inhibitor Y-27632, and 0.4% PVA on seeding day only to facilitate aggregation). We exchanged the medium every other day, FGF2 and Y-27632 were withdrawn on day 6. Neural induction was started on day 8 (neural induction medium: DMEM-F12 with KOSR 1:5, 1%penicillin/streptomycin/glutamine, 1% non-essential amino acids, N2 supplement 1:100, and Heparin (Sigma-Aldrich) 1 μg/mL) and continued for 6 days with media changes every other day. On day 13, we embedded the aggregates into 30 μL matrigel droplets and transferred them to 6 cm2 suspension tissue culture dishes (Sarstedt) in cerebral organoid differentiation medium (DMEM-F12 and Neurobasal 1:1 with 1% penicillin/streptomycin/glutamine, 1% NEAA, N2 upplement 1:200, B27 supplement without vitamin A 1:100, Insulin (Sigma-Aldrich) 1:4000, and 2-mercaptoethanol 1:285714). We placed the culture dishes on a shaker at 37°C and 5% CO2 and fed the organoids every other day. On day 20 the B27 supplement was replaced by B27 with Vitamin A (Thermo Fisher) and organoids were cultured until day 30 or 45.
Initial cleaved caspase 3-based detection of G418 toxicity (Figure 8—figure supplement 1a–f)
Request a detailed protocolAt day 50, we treated AMOs with increasing concentrations (0, 5, 50, 100, 250, 500, 1000 μg/mL) of G418 (Sigma-Aldrich) added directly to the culture medium. After 2 days, we renewed the medium (including identical toxin concentrations) and fixed the aggregates after a total of 4 days of treatment. Fixation, whole mount immunostaining for cCasp3 and Sox2 as well as BABB-based clearing was performed as outlined above. Image analysis followed the steps as outlined in the high-content analysis section with slight modifications to accommodate the individual brightness, morphology, and background characteristics of the cCasp3 staining. Briefly, after identifying AMOs and Sox2+ cells as described previously, the cCasp3 channel was background corrected by running a sliding parabola algorithm with a curvature setting of 10 across each confocal slice of the AMO. We identified apoptotic cells via the ‘find nuclei’ function in the 647 nm channel, algorithm 'M' and further selected them to be cCasp3+ if they were larger than 11 μm2, smaller than 100 μm2, and brighter than 2700 abu. We considered cells to be Sox2/cCasp3 double-positive if they fulfilled the criteria for both filters at the same time. The results were outputted to Microsoft Excel, reformatted and then transferred to GraphPad Prism v8.0.2 for plotting, data analysis, and curve fitting.
Additional G418 toxicity testing with a broader range of concentrations (Figure 8—figure supplement 1g–i)
Request a detailed protocolOn day 30 of differentiation, AMOs received Tox medium (TM) overnight to remove antioxidants present in the B27 media supplement. TM consisted of DMEM-F12 supplemented with 1% N2 and 1% penicillin/streptomycin/glutamine. On day 31, we added different concentrations of G418 (ranging from 0 to 10,000 μg/mL) dissolved in TM and renewed the medium (including the different compound concentrations) once two days later. After a total treatment time of 96 hr, the samples were analyzed by their size, the CellTiter-Glo 3D Cell Viability Assay, and cCasp3 staining as described before. In the case of cCasp3 staining, the image analysis parameters were used with minor modifications to account for the individual characteristics of the staining, for example signal intensity, aggregate size, background. The calculation of the organoid size/area plotted in Figure 8—figure supplement 1h was performed on the same fluorescence images as the cCasp3 quantification in subfigure g) according to the methods outlined in the description of the high content analysis. Data was transferred to Microsoft Excel for further analysis and GraphPad Prism v8.4.2 for plotting and curve fitting.
Data availability
All RNA sequencing data generated by us was deposited to the NCBI GEO database (GSE119060).
-
NCBI Gene Expression OmnibusID GSE119060. A fully automated high throughput-workflow for human neural organoids.
-
NCBI Gene Expression OmnibusID GSE66302. A human fetal transcriptional atlas.
-
ArrayExpressID E-MTAB-4868. Transcriptome profiling of DA neurons, human midbrain-like organoids and prenatal midbrain.
-
NCBI BioProjectID PRJNA245228. RNAseq data of 36 samples across human brain development by age group from LIBD.
References
-
Developmental disorders of the midbrain and hindbrainFrontiers in Neuroanatomy 6:7.https://doi.org/10.3389/fnana.2012.00007
-
Developing networks play a similar melodyTrends in Neurosciences 24:353–360.https://doi.org/10.1016/S0166-2236(00)01813-0
-
Controlling the false discovery rate: a practical and powerful approach to multiple testingJournal of the Royal Statistical Society: Series B 57:289–300.https://doi.org/10.1111/j.2517-6161.1995.tb02031.x
-
Pluripotent stem cell platforms for drug discoveryTrends in Molecular Medicine 24:805–820.https://doi.org/10.1016/j.molmed.2018.06.009
-
Functional properties of rat and human neocortical voltage-sensitive sodium currentsJournal of Neurophysiology 71:1052–1064.https://doi.org/10.1152/jn.1994.71.3.1052
-
A whole-mount immunocytochemical analysis of the expression of the intermediate filament protein vimentin in xenopusDevelopment 105:61–74.
-
The use of brain organoids to investigate neural development and diseaseNature Reviews Neuroscience 18:573–584.https://doi.org/10.1038/nrn.2017.107
-
Disease modeling in stem Cell-Derived 3D organoid systemsTrends in Molecular Medicine 23:393–410.https://doi.org/10.1016/j.molmed.2017.02.007
-
Three-Dimensional cell cultures in drug discovery and developmentSLAS Discovery : Advancing Life Sciences R & D 22:456–472.https://doi.org/10.1177/1087057117696795
-
Organoids as an in vitro model of human development and diseaseNature Cell Biology 18:246–254.https://doi.org/10.1038/ncb3312
-
Spheroid-based drug screen: considerations and practical approachNature Protocols 4:309–324.https://doi.org/10.1038/nprot.2008.226
-
The convergence of stem cell technologies and phenotypic drug discoveryCell Chemical Biology 26:1050–1066.https://doi.org/10.1016/j.chembiol.2019.05.007
-
Cerebrospinal fluid biomarkers of central catecholamine deficiency in Parkinson's disease and other synucleinopathiesBrain : A Journal of Neurology 135:1900–1913.https://doi.org/10.1093/brain/aws055
-
ScaleS: an optical clearing palette for biological imagingNature Neuroscience 18:1518–1529.https://doi.org/10.1038/nn.4107
-
Improved patch-clamp techniques for high-resolution current recording from cells and cell-free membrane patchesPflügers Archiv - European Journal of Physiology 391:85–100.https://doi.org/10.1007/BF00656997
-
Disease modeling using 3D organoids derived from human induced pluripotent stem cellsInternational Journal of Molecular Sciences 19:936.https://doi.org/10.3390/ijms19040936
-
Concave microwell array-mediated three-dimensional tumor model for screening anticancer drug-loaded nanoparticlesNanomedicine: Nanotechnology, Biology and Medicine 11:1153–1161.https://doi.org/10.1016/j.nano.2015.02.009
-
Method for generation of homogeneous multicellular tumor spheroids applicable to a wide variety of cell typesBiotechnology and Bioengineering 83:173–180.https://doi.org/10.1002/bit.10655
-
Comparison of cell and Organoid-Level analysis of Patient-Derived 3D organoids to evaluate tumor cell growth dynamics and drug responseSLAS Discovery : Advancing Life Sciences R & D 25:744–754.https://doi.org/10.1177/2472555220915827
-
Statistical practice in high-throughput screening data analysisNature Biotechnology 24:167–175.https://doi.org/10.1038/nbt1186
-
Modulation of dopamine neuronal activity by glutamate receptor subtypesNeuroscience & Biobehavioral Reviews 21:511–518.https://doi.org/10.1016/S0149-7634(96)00030-9
-
MPTP mouse models of Parkinson's disease: an updateJournal of Parkinson's Disease 1:19–33.https://doi.org/10.3233/JPD-2011-11023
-
High-Content monitoring of drug effects in a 3D spheroid modelFrontiers in Oncology 7:293.https://doi.org/10.3389/fonc.2017.00293
-
Human neurospheres as three-dimensional cellular systems for developmental neurotoxicity testingEnvironmental Health Perspectives 117:1131–1138.https://doi.org/10.1289/ehp.0800207
-
Tyrosine hydroxylase: human isoforms, structure and regulation in physiology and pathologyEssays in Biochemistry 30:15–35.
-
Directed differentiation of cholangiocytes from human pluripotent stem cellsNature Biotechnology 33:853–861.https://doi.org/10.1038/nbt.3294
-
Drug discovery through stem cell-based organoid modelsAdvanced Drug Delivery Reviews 69-70:19–28.https://doi.org/10.1016/j.addr.2014.02.006
-
Machine learning and image-based profiling in drug discoveryCurrent Opinion in Systems Biology 10:43–52.https://doi.org/10.1016/j.coisb.2018.05.004
-
Fiji: an open-source platform for biological-image analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019
-
Making sense of the multiple MAP-2 transcripts and their role in the neuronMolecular Neurobiology 16:149–162.https://doi.org/10.1007/BF02740642
-
Ionic channel currents in cultured neurons from human cortexJournal of Neuroscience Research 34:170–178.https://doi.org/10.1002/jnr.490340204
-
The 6-hydroxydopamine model of Parkinson's diseaseNeurotoxicity Research 11:151–167.https://doi.org/10.1007/BF03033565
-
Advanced CUBIC protocols for whole-brain and whole-body clearing and imagingNature Protocols 10:1709–1727.https://doi.org/10.1038/nprot.2015.085
-
Growth of multicell spheroids in tissue culture as a model of nodular carcinomasJournal of the National Cancer Institute 46:113–120.https://doi.org/10.1093/jnci/46.1.113
Article and author information
Author details
Funding
H2020 European Research Council (grant agreement No [669168])
- Henrik Renner
- Martha Grabos
- Mandy Otto
- Dagmar Zeuschner
- Hans R Schöler
- Jan M Bruder
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Britta Trappmann and Martin Weiss for their help with the spinning disc confocal imaging for the calcium activity experiments, as well as Laura Gonzales Cano, Yotam Menuchin-Lasowski, Michele Boiani, and Eva Kutejova for helpful discussion of the manuscript.
Copyright
© 2020, Renner et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 17,712
- views
-
- 1,983
- downloads
-
- 134
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Neuroscience
- Stem Cells and Regenerative Medicine
Fibro-adipogenic progenitors (FAPs) are muscle-resident mesenchymal progenitors that can contribute to muscle tissue homeostasis and regeneration, as well as postnatal maturation and lifelong maintenance of the neuromuscular system. Recently, traumatic injury to the peripheral nerve was shown to activate FAPs, suggesting that FAPs can respond to nerve injury. However, questions of how FAPs can sense the anatomically distant peripheral nerve injury and whether FAPs can directly contribute to nerve regeneration remained unanswered. Here, utilizing single-cell transcriptomics and mouse models, we discovered that a subset of FAPs expressing GDNF receptors Ret and Gfra1 can respond to peripheral nerve injury by sensing GDNF secreted by Schwann cells. Upon GDNF sensing, this subset becomes activated and expresses Bdnf. FAP-specific inactivation of Bdnf (Prrx1Cre; Bdnffl/fl) resulted in delayed nerve regeneration owing to defective remyelination, indicating that GDNF-sensing FAPs play an important role in the remyelination process during peripheral nerve regeneration. In aged mice, significantly reduced Bdnf expression in FAPs was observed upon nerve injury, suggesting the clinical relevance of FAP-derived BDNF in the age-related delays in nerve regeneration. Collectively, our study revealed the previously unidentified role of FAPs in peripheral nerve regeneration, and the molecular mechanism behind FAPs’ response to peripheral nerve injury.
-
- Cell Biology
- Stem Cells and Regenerative Medicine
Our understanding of the transitions of human embryonic stem cells between distinct stages of pluripotency relies predominantly on regulation by transcriptional and epigenetic programs with limited insight on the role of established morphological changes. We report remodeling of the actin cytoskeleton of human embryonic stem cells (hESCs) as they transition from primed to naïve pluripotency which includes assembly of a ring of contractile actin filaments encapsulating colonies of naïve hESCs. Activity of the Arp2/3 complex is required for the actin ring, to establish uniform cell mechanics within naïve colonies, promote nuclear translocation of the Hippo pathway effectors YAP and TAZ, and effective transition to naïve pluripotency. RNA-sequencing analysis confirms that Arp2/3 complex activity regulates Hippo signaling in hESCs, and impaired naïve pluripotency with inhibited Arp2/3 complex activity is rescued by expressing a constitutively active, nuclear-localized YAP-S127A. Moreover, expression of YAP-S127A partially restores the actin filament fence with Arp2/3 complex inhibition, suggesting that actin filament remodeling is both upstream and downstream of YAP activity. These new findings on the cell biology of hESCs reveal a mechanism for cytoskeletal dynamics coordinating cell mechanics to regulate gene expression and facilitate transitions between pluripotency states.