Specification of positional identity in forebrain organoids.

Human brain organoids generated with current technologies recapitulate histological features of the human brain, but they lack a reproducible topographic organization. During development, spatial topography is determined by gradients of signaling molecules released from discrete signaling centers. We hypothesized that introduction of a signaling center into forebrain organoids would specify the positional identity of neural tissue in a distance-dependent manner. Here, we present a system to trigger a Sonic Hedgehog (SHH) protein gradient in developing forebrain organoids that enables ordered self-organization along dorso-ventral and antero-posterior positional axes. SHH-patterned forebrain organoids establish major forebrain subdivisions that are positioned with in vivo-like topography. Consistent with its behavior in vivo, SHH exhibits long-range signaling activity in organoids. Finally, we use SHH-patterned cerebral organoids as a tool to study the role of cholesterol metabolism in SHH signaling. Together, this work identifies inductive signaling as an effective organizing strategy to recapitulate in vivo-like topography in human brain organoids.

INTRODUCTION

Human pluripotent stem cells (hPSCs) have the intrinsic capacity to self-organize into multicellular organ-like structures called organoids[1,2]. Brain organoids recapitulate the cellular diversity and micro-architectural features characteristic of discrete brain regions, providing unprecedented opportunities to model human brain development and disease[3-5]. However, in these organoids, individual brain regions are ordered randomly, non-reproducibly and lack the characteristic in vivo anterior-posterior, dorso-ventral, and medio-lateral positioning that supports the emergence of complex brain structure and function[1,4,6]. The absence of a defined topography is a crucial shortcoming of current brain organoid technologies.

During brain development, topographic maps are generated by gradients of signaling activity that allow cells to acquire discrete regional identities as a function of their position[7,8]. Region-specific organoids attempt to overcome positional heterogeneity by restricting cellular identities to a single brain area such as the optic cup[9], adenohypophysis[10], forebrain[11], midbrain[12], or hindbrain[13], using bath-application of patterning factors[14]. More recently, several groups have created a dorso-ventral axis by fusing dorsal and ventral forebrain organoids[15-17]. While these strategies recapitulate some processes associated with long-range tissue interactions, such as cell migration, they fail to capture the full continuum of regional diversity encompassed by the human forebrain.

RESULTS

Establishing a SHH protein gradient in 3D hPSC culture

Sonic Hedgehog (SHH) is a signaling factor whose graded expression specifies the spatial organization of discrete progenitor domains across the neuraxis[18,19]. To test whether introduction of a SHH signaling gradient into developing forebrain organoids could specify positional domains, we engineered an inducible SHH-expressing hPSC line (iSHH) that could be embedded at one pole of a developing organoid (Fig. 1a). The iSHH line was generated by TALEN-mediated gene targeting into the AAVS1 locus of an hPSC line with constitutive RFP expression[20,21]. A constitutively active reverse tetracycline transactivator is expressed from one allele of the AAVS1 locus, while full-length human SHH is expressed from the other allele, under the control of a tetracycline response element (TRE) (Fig. 1a). This strategy resulted in titratable SHH expression (Fig. 1b) and commensurate post-translational palmitoylation (Supplementary Figure 1a). The iSHH line induced robust expression of the SHH target gene FOXA2 during neural differentiation in the presence of doxycycline, indicating appropriate biological activity (Fig. 1c).

Figure 1.

Establishment of dorso-ventral and antero-posterior developmental axes in SHH organoids. (a) Strategy to specify positional identity in forebrain organoids. SHH exhibits a ventral-high to dorsal-low gradient in the developing forebrain. A doxycycline-inducible SHH (iSHH) hPSC line was engineered using TALEN-mediated homologous recombination. The iSHH line could be positioned at one pole of developing forebrain organoids. (b) qRT-PCR analysis shows robust and titratable induction of SHH transcript expression in iSHH but not H9 hPSC line. Individual replicates plotted on bar graphs, mean ± S.D. One-way ANOVA with Dunnett Test. *** P=0.0003, **** P=0.0001. N=3 for each iSHH condition, N=4 for each H9 condition. (c) qRT-PCR analysis of SHH and FOXA2, a SHH target gene, during neural differentiation shows that iSHH expresses biologically active protein, comparable to levels from hPSC-derived floorplate. Individual replicates plotted on bar graphs, mean ± S.D. One-way ANOVA with Dunnett Test. ***P<0.001, **P<0.01. N=3 differentiations for all conditions. (d) Schematic of the method for generating forebrain organoids with a SHH protein gradient. Timeline and small molecules used for neural induction are shown. The specific timing of matrigel embedding can exhibit slight variation among differentiations. (e) Visualization of SHH protein gradient using tyramide signal amplification (TSA). A line scan was used to quantify SHH protein signal. Graph shows mean ± 95% confidence interval. Intensity of SHH signal (y-axis, green) normalized to maximum iSHH fluorescence. Origin (X-axis) set to position of maximum iSHH fluorescence. (f) Quantification of the distance at which the iSHH-organizer signal (ORG) and SHH protein signal (SHH) reach 25% of maximum intensity demonstrates that the SHH protein gradient extends beyond the limit of the organizer cells. Student two-tailed t-test. N=6 spheroids, 2 batches. (g) Schematic illustrating that PAX6 and NKX2.1 define dorsal and ventral forebrain, respectively. H9 organoids (H9) or SHH-organoids grown without doxycycline (−DOX) largely express PAX6 without detectable NKX2.1 expression, indicating dorsal identity. SHH-organoids grown in the presence of doxycycline (400 ng/ml, +DOX) suppress PAX6 and induce NKX2.1 near the iSHH organizer, indicating induction of ventral identity. Representative images are shown. H9, N=7 organoids; −DOX, N=8 organoids; SHH+Dox, N=11 organoids. Samples are from at least 2 batches. (h) Schematic illustrating that FOXG1 and OTX2 define anterior and posterior forebrain domains. H9 organoids (H9) or SHH-organoids grown without doxycycline (−DOX) largely express FOXG1 with only small amounts of OTX2 expression, indicating predominant telencephalic (anterior forebrain) identity. SHH-organoids grown in the presence of doxycycline (400 ng/ml, +DOX) suppress FOXG1 and induce OTX2 near the iSHH organizer, indicating induction of diencephalic (posterior forebrain) identity. Representative images are shown. H9, N=8 organoids; −DOX, N=8 organoids; +DOX, N=12 organoids. Samples are from at least 2 batches. (i) Schematic illustrating diencephalic expression of TCF7L2, SIX3, and DKK1. Organoids grown without doxycycline (−DOX) do not express TCF7L2 or SIX3. Organoids grown in the presence of doxycycline (400 ng/ml, +DOX) induce TCF7L2 and SIX3 at day 6. −DOX, N=6 organoids; +DOX, N=6 organoids. Samples are from at least 2 batches. (j) Dose-dependent activation of DKK1 by SHH. N=6-8 for each dose, organoids from two batches. Individual replicates plotted on bar graphs, mean ± S.D. One-way ANOVA, P < 0.0001, with Dunnett Test. **** P < 0.0001, ** P < 0.01. Scale Bars: 200 μm.

We next developed a method to embed iSHH cells at one pole of an hPSC spheroid, mimicking a developmental organizer. First, 1,000 iSHH cells were seeded in low-attachment round bottom microwells and allowed to aggregate for 24 hours (day −2 to −1). The next day 10,000 wild-type H9 hPSC were seeded on top of the iSHH cells (day −1 to 0) and allowed to aggregate for another 24 hours, resulting in a reproducible chimeric iSHH-H9 spheroid in which a small cluster of iSHH cells are embedded within a larger spheroid (Supplementary Figure 1b, >90% efficiency). We refer to these chimeric 3D hPSC cultures as SHH-spheroids. For the specific induction of forebrain identity, we devised a strategy that combines aspects of several previous organoid induction protocols. In brief, SHH-spheroids were cultured in the presence of inhibitors of TGFβ, BMP, and WNT (3-inhibitor protocol) for 6-8 days to promote specification of anterior forebrain identity[11]. After 6-8 days of differentiation, spheroids were embedded in matrigel droplets to promote neuroepithelial organization and moved to an orbital shaker upon outgrowth of neuroepithelial buds, which was typically 4-6 days[4] (Supplementary Figure 1c). Doxycycline (400ng/ml) was added to the differentiation medium starting at day 0 to induce an asymmetric SHH signal (Fig. 1d). Under those conditions, the iSHH-organizer typically remained positioned at one end of the developing organoid, though we observed small clusters of organizer cells separated from the main organizer in 25% of instances (Supplementary Figure 1d). To assess whether this strategy generates a gradient of SHH protein, we performed immunocytochemistry at day 4 of the differentiation. In the absence of doxycycline, no SHH protein was detectable. In the presence of doxycycline, high levels of SHH expression were detected in the region of the iSHH cells, and the abundance of SHH protein was reduced as a function of distance from the SHH-organizer (Fig. 1e,f).

Specification of developmental axes in SHH-organoids

Functional subdivisions in the mammalian forebrain are located at discrete positions within a Cartesian coordinate system, whose axes are arranged in the dorso-ventral, medio-lateral, and antero-posterior directions. Each of the dorso-ventral and antero-posterior domains can be identified by the expression of characteristic transcription factors. PAX6 is expressed in the dorsal forebrain and NKX2.1 in the ventral forebrain (Fig. 1g). FOXG1 is expressed in the anterior forebrain, or telencephalon, while it is absent from the posterior forebrain, or diencephalon. In contrast, OTX2 is expressed in the diencephalon (Fig. 1h). When we differentiated H9 spheroids or SHH-spheroids in the absence of doxycycline, they predominantly expressed PAX6 and FOXG1 within self-organized neuroepithelium, in agreement with previous findings[11] (Fig. 1g,h). This indicates that the default identity of our 3D cultures is dorsal-anterior forebrain, which gives rise to the neocortex. SHH-spheroids grown in the presence of doxycycline (400 ng/ml) induced NKX2.1 near the organizer, while PAX6 expression was suppressed (Fig. 1g).

Surprisingly, FOXG1 expression was also suppressed near the organizer, while OTX2 expression was maintained (Fig. 1h), suggesting that SHH can lead to posterior forebrain specification. We therefore examined the expression of TCF7L2, which is expressed in the diencephalon[22,23] and SIX3[24], which is expressed in optic recess and hypothalamus, in SHH organoids. Both TCF7L2 and SIX3 were induced near organizer tissue, supporting a role for SHH in anterior-posterior patterning, as well as dorsoventral patterning. SHH exhibited dose-dependent induction of DKK1 (Fig. 1j), a secreted antagonist that shapes the anterior-posterior WNT gradient in vivo, that is commonly induced in regions of high WNT/β–catenin activity and that is a critical regulator of anterior-posterior patterning[25-27]. Thus, SHH might exert anterior-posterior patterning activity via regulation of WNT signaling[28].

Topographic organization of SHH-organoids

Given that a SHH protein gradient could specify distinct dorso-ventral and antero-posterior positional domains, we next sought to determine whether the resulting forebrain subdivisions were aligned according to an anatomically appropriate topography. The position of the presumptive forebrain subdivisions in vivo can be defined by their distance from the SHH source in the hypothalamus. The neocortex (PAX6+/FOXG1+) represents the anterior and dorsal subdivision of the forebrain which is located most distal from the SHH source. Just ventral to the neocortex is the lateral ganglionic eminence (LGE) characterized by the co-expression expression of GSH2+/FOXG1+, followed by the medial ganglionic eminence (MGE) characterized by co-expression of NKX2.1+/FOXG1+, the antero-dorsal hypothalamus characterized by expression of NKX2.2+, and finally the ventro-posterior hypothalamus characterized by the expression of NKX2.1+/FOXG1-, representing the brain region located most ventral and proximal to the SHH source (Fig. 2a).

Figure 2.

In vivo-like topographic organization of major forebrain subdivisions in SHH-organoids. (a) In situ hybridization images depicting forebrain topography. The neocortex co-expresses FOXG1 (blue) and PAX6 (green), and is located distal to the SHH source (red). The lateral ganglionic eminence (LGE) co-expresses FOXG1 and GSH2 (orange), and is just ventral to the neocortex. The medial ganglionic eminence (MGE) co-expresses FOXG1 and NKX2.1 (light purple), and is just ventral to the LGE. NKX2.1 is also expressed in the ventro-posterior hypothalamus (dark purple), and is distinguished from the MGE by exclusion of FOXG1 expression. The anterior hypothalamus expresses NKX2.2 (cyan) and resides in the gap between MGE and Hypothalamic NKX2.1 expression. SHH (red) is strongly expressed in the ventro-posterior hypothalamus, and is also weakly expressed in the ventral MGE. Images are from © 2008 Allen Institute for Brain Science. Allen Developing Mouse Brain Atlas. Available from developingmouse.brain-map.org. (b) Quantification of position of regional domains in SHH organoids (400 ng/ml doxycycline). Immunofluorescent signals of each regional marker are plotted as a function of distance from the iSHH organizer (see methods and Supplementary Figure 2 for details of quantification). Mean ± S.E.M. N=30 organoids from 6 batches. (c) Sections from a single organoid showing that SHH-organoids develop with in vivo-like topography. A neocortex-like domain that co-expresses FOXG1 and PAX6 is distal to the SHH organizer (red). A GSH2 domain is observed partially overlapping with and adjacent to the PAX6 domain. A FOXG1/NKX2.1 co-expression domain occupies an intermediate position between the neocortex-like domain and the SHH organizer. Finally, NKX2.2 and NKX2.1 are expressed in the FOXG1 negative territory, with NKX2.1 most proximal to the SHH organizer. Colors of arrowheads correspond to regional identities in panel 2a and delineate presumptive forebrain domains. (d) SHH-organoids grown without doxycycline treatment mainly express FOXG1 (N=6 organoids, 2 batches) and PAX6 (N=6 organoids, 2 batches), with some GSH2 expression (N=9 organoids, 3 batches). No NKX2.1 or NKX2.2 is observed (N=10 organoids, 3 batches). (e) The length of ZO-1 apical neuroepithelia differs between PAX6+ and NKX2.1+ regions. PAX6+ regions contain circular, rosette-like neuroepithelia, while NKX2.1+ regions contain highly-extended neuroepithelia. High-magnification images correspond to insets. Quantification shows average neuroepithelium length in PAX6 versus NKX2.1 domains. Dots represent individual organoids. Student two-tailed t-test. N=7 organoids, 2 batches for each condition. Scale bars: 200 μm (Low magnification), 50 μm (High magnification).

When differentiated in the presence of doxycycline (400 ng/ml), SHH-spheroids formed at least five topographically distinct presumptive forebrain domains. Quantification of the distance of forebrain domains from the SHH-organizer (Supplementary Figure 2), identified by protein expression, revealed that the positioning of these domains mimics the pattern observed in vivo (Fig. 2b). Topographic organization could be triggered in a similar manner in at least five additional hPSC or iPSC lines using the same H9-based iSHH organizer, though the specific growth rates and size of regional domains could differ (Supplementary Figure 3). We next identified presumptive neocortical (FOXG1+/PAX6+), LGE (GSH2+), MGE (FOXG1+/NKX2.1+), anterior hypothalamic (NKX2.2+) and ventro-posterior hypothalamic (NKX2.1+/FOXG1-) regions by immunocytochemistry (Fig. 2c). FOXG1 expression was absent from the organizer tissue itself, but robust expression of NKX2.1 was observed in at least a subset of the organizer cells, suggesting hypothalamic identity[29] (Supplementary Figure 4a). The resulting topographic patterning is in contrast to the patterning exerted by bath application of SHH agonist to forebrain organoids, which results in subpallium-restricted identities[15,16]. Topographic patterning could also be achieved in the absence of matrigel (Supplementary Figure 4b). SHH-spheroids differentiated in the absence of doxycycline expressed markers of dorsal-anterior forebrain, including PAX6 and FOXG1, with a small proportion of GSH2 expressing cells. These 3D cultures in the absence of doxycycline did not express the ventral identity genes NKX2.1 and NKX2.2 (Fig. 2d).

Interestingly, tissue cytoarchitecture differed between areas that are proximal (hypothalamic-like) and distal (telencephalic-like) to the SHH source. At day 20, PAX6+ distal neuroepithelia acquired circular, rosette-like morphologies (Fig. 2e, inset), consistent with the self-organizing radial structures described for forebrain organoids[4,11]. In contrast, the NKX2.1+ proximal region contained many thin and highly extended neuroepithelia (Fig. 2e, inset) as observed previously for 3D structures of early hypothalamic lineages[10,30]. This appears to be restricted to diencephalic NKX2.1+ regions, as presumptive MGE domains (FOXG1+/NKX2.1+) typically acquired circular, rosette-like morphologies (Fig. 2c Supplementary Figure 4c). Our data indicate that an asymmetric SHH cue enables the ordered patterning and topographical organization of brain regions without interfering in the self-organization process, which leads to the establishment of region-specific tissue microarchitectures. We term these 3D cultures SHH-organoids.

Cell-type diversity within SHH-organoids

To explore the degree to which topography is maintained over time, we cultured SHH-organoids for up to 70 days. FOXG1+ and OTX2+ domains remained largely discrete at this stage, typically with the OTX2+ domain located proximal to the organizer, while the FOXG1+ domain was located distal to the organizer (Supplementary Figure 5a). Radially organized PAX6, TBR2, TBR1 cerebral cortex-like tissue and DARPP32+ striatum-like tissue emerged in distal domains (Supplementary Figure 5b-d), while hypothalamic-like tissue, expressing LHX6, OTP, POMC, and TH, was found in the immediate vicinity of the organizer cells (Supplementary Figure 5d-f). In some instances organoid topography was more difficult to discern because the organizer seemed to have dispersed throughout the tissue. In addition, tissue necrosis in the center of those larger organoids may degrade organization over time.

The identity of LHX6+ cells including the presence of cortical interneurons was further defined by differentiating SHH-organoids using an LHX6-citrine hPSC line. As expected, LHX6+, putative hypothalamic precursors, emerge near the organizer, with some cells co-expressing NKX2.1 (Supplementary Figure 6a). A subset of LHX6+ cells co-expresses FOXG1+ (Supplementary Figure 6b), while a subset is negative for FOXG1, suggesting the presence of both MGE and hypothalamic LHX6 lineages. Some FOXG1+/LHX6+ cells exhibit a leading-process morphology characteristic of migrating interneurons (Supplementary Figure 6c). Cells expressing SST, PV, or CR were observed in SHH-organoids (Supplementary Figure 6d-f), suggesting the emergence of diverse interneuron populations. However, PV+ cells did not co-express LHX6 indicating that those PV+ cells many not represent cortical interneurons.

Mechanism of long-range SHH signaling

We next sought insight into the mechanisms by which a SHH signaling gradient establishes forebrain topography. Using the suppression of PAX6 as a read-out of SHH signaling activity, we compared the radial extent of PAX6 suppression in 2D versus 3D cultures. High-dox conditions (2mg/ml) were used in assessing SHH activity. In 3D cultures, we observe the suppression of PAX6 expression at an average distance of 335 ± 54 μm from the organizer tissue (Fig. 3a). Interestingly, this range is similar to the distances reported for SHH signaling in vivo based on studies across multiple tissues of the developing mouse or chick embryos[31-33]. However, in 2D cultures in which a similar 1:10 ratio of iSHH to wildtype cells is maintained, we saw suppression of PAX6 at much shorter average distances (75 ± 31 μm) from the SHH source tissue (Fig. 3a).

Figure 3.

Characterization of long-range SHH signaling mechanism. (a) Comparison of distance of SHH-signaling activity after 12 days of differentiation in 3D and 2D culture using 2ug/ml doxycycline. In 3D culture, SHH acts over 335 ± 54 μm from the iSHH-organizer, assessed by suppression of PAX6 expression. In 2D culture, SHH acts 75 ± 31 μm from the iSHH-organizer cells. Insets show high-magnification images. Quantification shows mean ± S.D. Individual replicates are plotted on graphs. Student two-tailed t-test. **** P < 0.0001. N=5 organoids, 2 batches (3D), n = 4 differentiations (2D). (b) Depiction of experiment to distinguish between two theoretical mechanisms for long-range SHH activity in 3D culture. A shell of GFP-expressing hPSCs is embedded around the iSHH organizer (see methods). A temporal model is supported if the all of SHH activity remains encapsulated within the GFP shell. A concentration-dependent model is supported if SHH activity extends past the GFP shell. (c) After 20 days of differentiation in 2μg/ml doxycycline, SHH activity is observed to extend past the GFP shell. The green dotted circle indicates the boundary of the GFP shell and the purple dotted circle indicates the boundary of the SHH activity, assessed by suppression of PAX6 activity. Quantifications show that the GFP shell extends 268 ± 24 μm from the iSHH-organizer, while SHH activity extends 508 ± 86 μm from the iSHH-organizer. Quantification shows mean ± S.D. Student two-tailed t-test. * P = 0.028. Individual replicates are plotted on graphs. N=5 organoids, 2 batches. Scale bars: 200 μm (low magnification), 50 μm (high magnification).

The observation that the range of SHH signaling activity is restricted in 2D raises the question how SHH establishes long-range gradients in 3D. We generated 2 theoretical models based on known mechanisms of SHH signaling in vivo (Fig. 3b). The first model is a concentration-dependent model. In this model, SHH signaling activity is determined by the local concentration of SHH morphogen, and is dependent on mechanisms that can transport SHH protein over long distances, such as diffusion, facilitated transport, or cell-to-cell relay[31,33-37]. The second model is a temporal model. In this model, SHH signaling is restricted to tissue immediately adjacent to the source, which is relieved of SHH signaling as tissue grows away from the source[38]. Tissue that grows away from a SHH source will record a brief duration of signaling activity, while tissue that remains near the SHH source will record a long duration of signaling activity. This temporal model does not rely on transport of SHH over long distances.

To experimentally test these models, we embedded the iSHH organizer within a circumscribed ring of GFP-expressing hPSCs. This was achieved by sequential plating of 1,000 iSHH cells, followed by 1,000 GFP cells 8 hours later, followed by 10,000 H9 cells the next day. We reasoned that if the activity of SHH extends past the GFP boundary, this would support a long-range concentration-dependent mechanism. On the other hand, if the GFP boundary expands coincident with the extent of SHH activity, this would support a temporal model in which tissue growth rather than protein transport establishes long-range SHH activity (Fig. 3b). After growing SHH-GFP-organoids for 20 days in doxycycline (2mg/ml), we observed that the GFP-expressing cells stayed within 285 ± 44 μm to the SHH organizer cells. The activity of SHH, determined as suppression of PAX6 expression, extended to 430 ± 96 μm, past the green boundary (Fig. 3c). These data provide evidence for a concentration-dependent signaling mechanism in which SHH protein is transported long distances away from producing cells. Still, because the long-range SHH signal degrades in 2D culture while a short-range signal persists, our data do not exclude the possibility of a concurrent short-range signaling mechanism in which the activity SHH is restricted to the immediate vicinity of producing cells. In fact, during limb patterning, SHH is known to act by both short- and long-range signaling mechanisms, and such an interplay of short- and long-range signaling may also be critical during forebrain development[29,39-41].

SHH-organoids as a model of neurodevelopment

The reproducible topography of SHH-organoids might make them suitable as a tool to study signaling mechanisms during early forebrain patterning. Cholesterol is integrally linked to the SHH pathway in many ways. Cholesterol is an agonist of the smoothened receptor[42,43]. SHH long-range signaling is dependent on the cholesterol post-translational modification of SHH[34,35,44], cholesterol rich lipid rafts[36,45,46], and can be facilitated by cholesterol-rich macromolecular soluble complexes[47,48].

To explore the multifaceted relationship of cholesterol to SHH signaling, we treated SHH-organoids (400ng/ml doxycycline) from day 0-10 with small molecules that inhibit cholesterol synthesis at different points along the biosynthetic pathway (Fig. 4a). AY9944 is a specific inhibitor of 7-dehydrocholesterol reductase (DHCR7), which facilitates the terminal step of cholesterol synthesis. Lovastatin is an inhibitor of HMG-CoA reductase, an early step in cholesterol biosynthesis. Cyclopamine, a Smoothened antagonist that blocks all SHH activity, was used as a control. At day 20 in the no drug condition, SHH-organoids established distinct positional domains marked by the expression of PAX6, GSH2, NKX2.2, and NKX2.1 at discrete distances from the SHH-organizer (Fig. 4b). As expected, SHH-organoids grown in the presence of cyclopamine (1 μM[49]) were nearly uniform for PAX6 expression, with some GSH2 expression, consistent with a near complete inhibition of SHH signaling activity (Fig. 4b,c). AY9944 (1.25 μM[50]) treatment reduced the efficacy and range of SHH signaling activity, shifting the organoid to a more dorsal identity. NKX2.2 induction was largely blocked, while cells just adjacent to the organizer expressed GSH2 and PAX6 (Fig. 4b,c). Lovastatin (5 μM[50]) treatment caused similar, but milder effects, than AY9944. There was a partial block in NKX2.2 induction, with concomitant increase in GSH2 expression (Fig. 4b,c), consistent with a dorsal shift and reduction in the efficacy and range of SHH signaling activity. Notably, the organizer cells retained diencephalic NKX2.1 expression in the presence of AY9944 or lovastatin but not cyclopamine (Supplementary Figure 7a), suggesting that inhibition of cholesterol synthesis has minimal impact on local SHH activity.

Figure 4.

Statins perturb tissue growth and SHH signaling. (a) Experimental design to test effect of statins on SHH-organoid development. (b) Four distinct positional domains are identified in day 20 SHH-organoids (400 ng/ml doxycycline). SHH-organoids treated with cyclopamine (1 μM) largely express PAX6, with some GSH2 immunoreactivity, suggesting a near complete inhibition of SHH activity. SHH-organoids treated by AY9944 (1.25μM) fail to induce NKX2.2 and only retain NKX2.1 expression in the immediate vicinity of the organizer, suggesting a strong reduction the range of SHH-signaling activity. SHH organoids treated with Lovastatin (5 μM) exhibit a moderate reduction in NKX2.2 induction, with a concomitant increase in the relative area of the GSH2 expression domain. Summary of drug phenotypes and potential impact on SHH signaling activity is depicted at bottom of panel. (c) Quantification of area of each positional domain, relative to total area. Dots represent individual organoids, Error bars are S.D. One-way ANOVA with Dunnett Test. No drug N=8-14; Cyclopamine N=5-14; AY9944 N=8-16; Lovastatin N=7-13. * P < 0.05; *** P < 0.001, **** P < 0.0001. (d) Lovastatin produces a dose-dependent reduction in the distance at which NKX2.2 is induced from the organizer (red dotted boundary). Insets show NKX2.2 cells in boxed regions. Frequency histograms plot the relative distribution of organizer and NKX2.2 positive cells as a function of distance from the center of the organizer. Dotted lines represent S.E.M. Bar graphs quantify the average distance between the center of NKX2.2 and organizer distributions, illustrated by horizontal grey bars in frequency histograms. Dots represent individual organoids. Error bars are S.D. One-way ANOVA with Dunnett Test. No drug, N=23; Lovastatin (Lova) 2 μM N=18; Lova 5 μM N=24; Lova 20 μM N=23; Simvastatin (Simv) N=12; Atorvastatin (Atorv) N=10 Scale bars: 200 μm (Day 20): 50 μm (Day 6).

In addition to the effects on patterning, all drugs affected organoid growth. Overall organoid size was reduced in all treatment conditions (Fig 4b,c). However, AY9944 and Lovastatin did not reduce the size of the organizer itself (Supplementary Figure 7b), consistent with the interpretation that cholesterol-synthesis inhibition has differential effects on local and long-range SHH activity. To dissociate growth and patterning phenotypes, we analyzed organoids at day 6, during the early stages of SHH patterning, and prior to the major phase of tissue growth (which occurs after matrigel embedding). We used NKX2.2, a direct target of GLI[51], as a read-out of SHH signaling activity. By measuring the distance at which NKX2.2 is expressed from the SHH producing cells, we could infer the distance of SHH signaling activity. Lovastatin reduced the distance at which NKX2.2 is induced from the organizer in a dose-dependent manner (Fig. 4d). AY9944 strongly inhibited NKX2.2 induction at all concentrations tested (Supplementary Figure 7c). Neither AY9944 nor lovastatin grossly perturbed processing of SHH peptide length (Supplementary Figure 7d). Together, these data support a model in which inhibition of cholesterol synthesis by Lovastatin or AY9944 result in a complex developmental phenotype characterized by impaired growth and a reduction in the range of SHH signaling activity, though perhaps in distinct ways.

DISCUSSION

This study reveals that an asymmetric morphogenetic cue provides positional information from which in vivo-like topography can emerge in brain organoids. Furthermore, we present evidence that SHH uses a long-range signaling mechanism during forebrain organoid development, recapitulating its in vivo behavior[40]. The emergence of long-range activity in 3D that is not present in 2D cultures highlights the importance of tissue geometry on the spatial dynamics of signaling activity.

The strategy of generating topographically patterned organoids allowed us to investigate how alterations in cholesterol metabolism perturb SHH-dependent patterning. Our data suggest that while both AY9944 and Lovastatin impair SHH-dependent patterning, they may do so in different ways. AY9944 has previously been shown to block Smoothened activation[50], consistent with the reduction in NKX2.2 induction at day 6 (Supplementary Figure 7c). On the other hand, Lovastatin reduced the range of SHH signaling at day 6 without reducing the intensity of signaling. While future studies are required to determine mechanism of action, there are multiple possible explanations for this observation. First, statin treatment could disrupt lipid raft integrity[52-56], which is required for long-range hedgehog signaling[46] and perhaps SHH oligomerization at the membrane[57,58]. Second, statins could lead to removal of the C-terminal SHH cholesterol modification, which is required for long-range signaling[34,35,44]. While it is unlikely that statins block the initial cholesterylation of SHH during auto-cleavage[59], SHH can subsequently lose its C-terminal cholesterol through proteolytic cleavage, or ‘shedding’, at the cell membrane[60]. Notably, disruption of lipid rafts by sterol-reduction, either with statins or methyl-β-cyclodextrin, increases shedding of multiple membrane-associated proteins, including SHH [52,53,60]. Third, statin treatment could reduce the formation of cholesterol-rich multimeric SHH complexes involved in long-range transport[47,48]. Finally, it is possible that statins perturb SHH patterning through isoprenoid pathway inhibition or off-target effects given the high concentration of drug used. Statins have been reported to be associated with fetal maldevelopment, including holoprosencephaly, in a small number of cases[61]. Our data suggest that some of these cases may be due to impaired SHH signaling, although caution must be exercised given the high concentration of statin used in our experiments (5μM) compared to the typical plasma concentration of statin (15-20 nM)[62].

A key feature of the patterning strategy presented here is the establishment of well-ordered positional domains during early organoid development. In future studies it will be important to determine the extent to which the topography of these domains can be retained during long-term in vitro culture for the study of region-specific neuronal diversity as well as connectivity and function. It is possible that additional tissue engineering approaches will be required to facilitate the maintenance of topography over time and to prevent increasing cell mixing or inappropriate migration, processes that could blur the boundaries between positional domains over time. For example, manipulating soluble ECM components has been used to improve the maintenance of radially organized structures within organoids over time[63].

For disease modeling, current organoid technologies have been successful at modeling cortical progenitor dynamics that ultimately result in changes in brain size. Pre-patterning and fusion of pallial- and subpallial-specific spheroids have been used to study specific phenotypes related to cortical interneuron migration and connectivity[16]. The establishment of forebrain topography in organoids shown here opens the possibility of studying a wide range of phenotypes in a single organoid system. In particular, it would be interesting to study complex neurodevelopmental diseases that have been hypothesized to alter or derive from regional specification during forebrain patterning, including autism[64], epilepsy[65],or suprasellar pediatric gliomas[66]. In the future it will important to develop strategies, such as bio-engineering approaches to reduce the variability of topography in long-term cultures. Finally, it will be exciting to test whether our strategy could be used as a general approach to establish topographies across all regions of the CNS beyond the forebrain.

ONLINE METHODS

hPSC maintenance and monolayer differentiation

hPSC (H9 and derivatives (EF1α::GFP (GFP), EF1 α::RFP; TRE-SHH (iSHH), MEL1, HUES6, HUES8, LHX6-citrine[1]) and iPSC (J1 and 348) were maintained with Essential 8 medium or Essential 8 flex (E8, Thermo, A15117001 or A28558501) in feeder-free conditions on vitronectin (VTN-N) substrate (Thermo, A14700). hPSCs were passaged as clumps with EDTA solution (0.5 M EDTA/PBS). For neural monolayer differentiation, hPSCs were dissociated to single cells and plated on matrigel substrate (BD Biosciences, 354234) in E8 at a density of 250,000 cells/cm2 in the presence of ROCK inhibitor (Y-27632, 10 μM, Tocris 1254) (Day −1). From Day 0 to 7, cells were cultured in Essential 6 medium (E6, Thermo, A1516401) in the presence of TGFβ and BMP inhibitors (LDN193189, 100 nM, Stem Cell Technologies, 72142; SB431542, 10 μM, Tocris, 1614). From day 7 to 12, cells were cultured in E6 alone. Floorplate was specified by addition of SHH (200 ng/ml, R&D Systems 464-SH) and CHIR99021 (0.7 μM, Tocris 4423). Media was changed every day during the differentiation.

SHH inducible line

The iSHH hPSC was generated according to a previously described gene targeting strategy[2]. Briefly, two donor constructs were targeted to the first intron of the AAVS1 locus using TALEN-facilitated homologous recombination. One construct contained a constitutively expressed reverse tetracycline transactivator (M2rtTA) and the second construct contained a tetracycline inducible (TRE) SHH cassette. The SHH cassette was generated using full length human SHH cDNA (Genecopoeia T1004). Targeting was performed in a constitutive RFP expressing line (WA09 derived, EF1α::RFP) so iSHH cells could be visualized[20].

SHH-spheroid aggregation

1,000 iSHH cells were dissociated to single cell and aggregated in low-attachment round bottom microwells for 24 hours in the presence of ROCK inhibitor (Y-27632, 10 μM) and WNT inhibitor (XAV939, 5μM). For figure 3c, the GFP boundary was established by first plating iSHH cells and allowing them to aggregate for 8 hours. 1,000 GFP cells were then plated on top and allowed to aggregate overnight. The next day, 10,000 wildtype H9 hPSCs were dissociated to single cell and allowed to aggregate on top of the iSHH cells, in E8 medium with ROCK inhibitor (Y-27632, 10 μM) and WNT inhibitor (XAV939, 5μM, Tocris 3748).

SHH organoid differentiation

SHH-spheroids were cultured in E6 with 3-inhibitors (LDN193189, 100 nM; SB431542, 10 μM; XAV939, 5 μM) plus doxycycline (400 ng/mL) until spheroid tissue began to brighten and have smooth edges, typically 6-8 days. SHH-spheroids were then embedded in matrigel drops and cultured in a N2/B27 based organoid media as previously described[3], except NeuroBrew B21 (Miltenyi Biotec) was substituted for B27. SHH-spheroids were moved an orbital shaker once formation of translucent, neuroepithelial buds was observed, typically 4 - 7 days after embedding. Organoids that failed to form translucent, neuroepithelial buds were discarded. Half media was changed every other day.

Histology and immunocytochemistry

Day 20 organoids were fixed in 4% PFA overnight at 4°C, and washed three times with PBS the next day. After fixation, tissue was cryoprotected in 30% sucrose/PBS and sectioned at 30μm on a cryostat (Leica 3050 S). Sections were blocked for 30 minutes in 10%FBS, 1%BSA, 0.3%triton in PBS, and incubated as floating section in primary antibody overnight. The next day, sections were washed with PBS and incubated as floating section in secondary antibody for 3 hours at room temperature. Day 6 organoids were fixed in 4% PFA for 2 hours at 4°C, and washed three times with PBS. After fixation, tissue was cryoprotected in 30% sucrose/PBS and sectioned at 20μm and mounted directly onto slides prior to immunofluorescence assays.

Quantification of organoid patterning and SHH activity

Organoid topography was quantified by first dividing each image into a grid of 20μm × 20μm squares (day 20 analysis) or 5μm × 5μm squares (day 6 analysis) that form a Cartesian coordinate grid. Each square is a discrete region of interest (ROI) that is associated with an X and Y coordinate. The origin of the grid was calculated for each image as the “center of mass” of the iSHH-organizer cells: CMx = Σ mixi / M and CMy = Σ miyi / M, where m = grey value intensity of iSHH cells for individual ROIs, × or y = coordinate, and M = total intensity for all ROIs, and the sum includes all x or y coordinates in the image. Next, ROI’s that were positive for expression of a given protein (e.g. PAX6) were defined by thresholding staining intensity and determining ROIs with positive immunoreactivity. The linear distance from the origin to each ROI (sqrt((ROIx-CMx)2 + (ROIy-CMy)2)) was calculated. We attempted to analyze every sixth section for each regional marker per organoid (2-5 sections typically), thus sampling the organoids in 150μm intervals using uniform random sampling.

In figure 2c and 4d ROI frequencies are plotted as a function of distance from the center of mass. In figure 3c the GFP shell distance is defined as the radial distance from the center of mass at which the ratio of GFP ROIs / DAPI ROIs drops below 35%. The SHH-activity distance is the radial distance from the center of mass to the edge of the PAX6-negative territory, defined as the distance at which the ratio of PAX6-negative ROIs / DAPI ROIs dropped below 95% (PAX6-negative ROIs were obtained by subtracting PAX6 thresholded pixels from DAPI thresholded pixels). The 35% and 95% parameters were determined prior to analysis as accurately cut off values by comparing with manual distance measurements. Area measurements in figure 4b are determined by calculating the number of positive ROIs positive for a given gene (e.g. PAX6) divided by the ROIs positive for DAPI. Exclusion criteria for analysis: 1) organoid batches that failed to induce FOXG1 expression (Supplementary Figure 2b), 2) organoids containing a split organizer, in which we observe two independent patterning centers (Supplementary Figure 2c), 3) Sections that lacked organizer cells, a center of mass could not be defined, and 4) sections that had large necrotic centers as this could confound distance analysis.

The distance of SHH patterning activity in figure 3a was obtained by measuring the distance of lines from the edge of PAX6 territory to the nearest organizer cells. At least 5 lines were measured per section. All quantifications were performed using ImageJ. Microscopy was performed using a standard inverted epifluorescence microscope (Olympus IX71 or Zeiss Axio Observer). Images were acquired using Cell Sens (Olympus) or Zen Pro (Zeiss) software. Min, max and gamma (midtone) adjustments were applied uniformly to images during processing with Adobe Photoshop Creative Cloud.

RNA extraction and RT-PCR

RNA was extracted using TRIzol (Invitrogen, 15596026) followed by chloroform extraction. RNA was precipitated in isopropanol and resuspended in nuclease free ddH2O. cDNA synthesis was performed using 1 μg of RNA (iScript, Bio-Rad, 1708840). RT-PCR was performed with EvaGreen Supermix (Bio-Rad, 1725202). Three to four individual samples were combined into a single replicate when collecting RNA from spheroids.

Visualization of SHH protein gradient

SHH-spheroids were fixed in 4% PFA for 6 hours at 4C, then washed three times in PBS. Standard immunocytochemistry was performed on wholemount SHH-spheroids using anti-SHH primary antibody (1:100, clone 5E1, DSHB). Secondary antibody was amplified using a tyramide signal amplification kit (Thermo, B40941). Fluorescence labeled spheroids were visualized using an inverted epifluorescence microscope. Gradients were quantified using line scans that transected the organizer tissue. Quantifications were performed using ImageJ.

Cell-based SHH Palmitoylation Assay

Synthesis of [I125]iodopalmitate was carried out as previously described[4-6]. iSHH hPSC were dissociated to single cell and seeded as a confluent monolayer on matrigel substrate at a density of 250,000 cells / cm2 and differentiated for 9 days in E6 with LDN193189, (100 nM) and SB431542 (10 μM) at increasing concentrations of doxycycline (0, 125, 250, 500, 1000, 2000 ng/ml). On day 9, iSHH cells were incubated for 1 h at 37 °C with DMEM containing 2% dialyzed FBS and then labeled for 4 h with 10μCi of [I125]iodopalmitate. The cells were lysed and subjected to immunoprecipitation with an anti-SHH antibody, and the immunoprecipitates were diluted in 100 μl of 2× SDS-PAGE sample buffer containing 10 mM DTT. The samples were analyzed on a 12.5% SDS-PAGE gel, followed by phosphorimaging on a Typhoon FLA-7000 phosphorimager. An aliquot of each immunoprecipitate was analyzed for SHH protein expression by Western blotting.

SHH Western Blot

Cells were harvested and washed with PBS. Cell pellets were lysed in RIPA buffer (Boston Bioproducts Inc. NC9517624) supplemented with protease inhibitors (Thermo 78425) and ran through a syringe several times. Lysate was incubated for 30 minutes on ice and then centrifuged to isolate the supernatant. For quantification of protein concentration, 3ul of lysate was added to 300ul of Precision Red (Cytoskeleton, Inc. ADV02-A) and absorbance (600nm) was read on a plate reader and the concentration was calculated by: (Abssample – Absbackground) × 12.5. For Western blots, 10μg of protein was resuspended in 2X Laemmli buffer and denatured at 100 C for 5 minutes. Samples were loaded and run on a 4–12% pre-cast Bis-Tris gel (Thermo NP0322BOX) and transferred overnight onto Nitrocellulose membrane (Thermo 88018). Membranes were blocked in 5% milk and primary antibodies were incubated overnight at 4C. Blots were visualized using ECL (Perkin Elmer NEL104001EA).

Statistical Analysis

All reported measurements are from distinct samples. When comparing samples across treatment groups, samples were collected, processed, and analyzed in parallel. In the figure legends, “n” refers to an independent biological replicate (e.g. single organoid or single monolayer culture). “Organoid” refers to an individual, distinct organoid from 1 well of culture dish. “Batch” refers to a group of organoids that were differentiated in the same dish on the same day. See Life Sciences Reporting Summary for more details. Statistical tests and biological replicates for each experiment are listed in the figure legends and below.

Data Availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Figure	Test	Sample	N
1b	One-way ANOVA with Dunnett Test	H9	4 cultures / condition
		iSHH	3 cultures / condition
1c	One-way ANOVA with Dunnett Test	H9	3 cultures / condition
		iSHH	3 cultures / condition
		floorplate	3 cultures / condition
1e/1f	Two-tailed t-test	iSHH+Dox	6 spheroids
1g		H9	7 organoids
		−Dox	8 organoids
		+Dox	11 organoids
1h		H9	8 organoids
		−Dox	8 organoids
		+Dox	12 organoids
1i		−Dox	6 organoids
		+Dox	6 organoids
1j	One-way ANOVA with Dunnett Test	0 ng/ml dox	8 organoids / condition
1j	One-way ANOVA with Dunnett Test	80 ng/ml dox	8 organoids / condition
		400 ng/ml dox	7 organoids / condition
		2000 ng/ml dox	6 organoids / condition
2b,c			30 organoids
2d		FOXG1	6 organoids
		PAX6	6 organoids
		GSH2	9 organoids
		NKX2.1	10 organoids
		NKX2.2	10 organoids
2e	Two-tailed t-test	PAX6/ZO1	7 organoids
		NKX2.1/ZO1	7 organoids
3a	Two-tailed t-test	3D culture	5 organoids
		2D culture	4 cultures
3c	Two-tailed t-test	Day 20	5 organoids
4b	One-way ANOVA with Dunnett Test	No Drug/PAX6	8 organoids
		No Drug/GSH2	14 organoids
		No Drug/NKX2.1	14 organoids
		No Drug/NKX2.2	14 organoids
		Cyclopamine/PAX6	5 organoids
		Cyclopamine/GSH2	12 organoids
		Cyclopamine/NKX2.1	14 organoids
		Cyclopamine/NKX2.2	14 organoids
		AY9944/PAX6	8 organoids
		AY9944/GSH2	15 organoids
		AY9944/NKX2.1	16 organoids
		AY9944/NKX2.2	12 organoids
		Lova5/PAX6	7 organoids
		Lova5/GSH2	13 organoids
		Lova5/NKX2.1	13 organoids
		Lova5/NKX2.2	13 organoids
Fig 4d	One-way ANOVA with Dunnett Test	No drug	23 organoids
		Lova2	18 organoids
		Lova 5	24 organoids
		Lova20	23 organoids
		Simv	12 organoids
		Atorv	10 organoids