The mature cortex consists of billions of neurons organized into a six-layered (L 1-6) structure. These neurons interact via short and long-range connections to form the complex circuitry which results in the emergent property of cognition. This cortical layout is highly conserved across mammals1,2. The neurons of the cortex occupy two major classes, GABAergic (inhibitory; iN) and glutamatergic (excitatory), and from those two major classes, dozens to hundreds of subtypes form, identified through distinct neuronal processes, circuit membership and gene expression patterns3. Glutamatergic excitatory neurons possess axonal projections which synapse at various brain regions, and projection properties dominate cell type grouping in single-cell RNA data3. These neurons form circuits with their inputs serving to promote action potential firing and downstream neural activity. Typography and topography (i.e. cell type and location) inform neuron function. For instance, glutamatergic pyramidal tract neurons are most common in deep layer L5 (ventrally located) and are associated with executing voluntary movements and planning. These project to subcortical targets like the striatum, thalamus, tectum and pons4. In contrast, glutamatergic intratelenecephalic trajectories connect excitatory neurons between cortical layers or cortical brain regions, and span most layers3. GABAergic inhibitory neurons have two major subclasses which reflect their point of origin5. Adenosine Deaminase RNA Specific B2 (ADARB2+) expressing inhibitory neurons are formed in the caudal ganglionic eminence, whereas LIM Homeobox 6 (LHX6+) neurons form in the medial ganglionic eminence (CGE and MGE, respectively). Inhibitory neurons produce the small molecule GABA, and modulate neuronal circuitry through dampening neuronal firing. Glia, which have previously been considered to be connective cells, have been found to serve critical roles in maintaining synapse integrity and cortical function. Astrocytes mediate blood brain barrier and have been shown to mediate synapse formation, elimination, and plasticity. Oligodendrocytes enable salutatory conduction of actionable potentials throughout the brain. Microglia are resident immune cells, participate in phagocytosis and inflammation response. Neuronal circuitry is an active field of study – which benefits from understanding the unique states and types of neurons and glia6.
The generation of cortical circuitry is critical. Dysfunction during corticogenesis has been implicated in multiple neurodevelopmental disorders7–9, and evolutionary changes between humans and other apes have been linked to the rapid expansion of the human cortex10. In early embryonic development (GW4 in humans; E10 in mice; where “E’ is embryonic day post conception and “GW” is gestational week), the ectodermal neural tube expands and compartmentalizes into the prosencephalon (forebrain), the mesencephalon (midbrain) and the rhombencephalon (hindbrain)11. The maturing forebrain subdivides along the dorsal-ventral axis into the pallium and subpallium, respectively. The pallium generates the bulk of the cerebral cortex and the subpallium forms the MGE and CGE (Figure 3)12. As the pallium develops, neuroepithelial cells differentiate to RG, named as such for their radial projection from the ventricular zone (VZ) towards the dorsal surface of the pallium, and for their combined marker set of neuroepithelial and astroglial expression patterns13. RG divide asymmetrically, both producing newborn RG to replenish the pool of stem cells, as well as forming intermediate progenitors (IPs) and a subset newborn Cajal-Retzius neurons (CR) directly. However most CR neurons are born exterior to the developing pallium and migrate tagentially into the marginal zone (MZ). IPs move dorsally to populate the subventricular zone (SVZ) while CRs migrate further to develop in the cortical plate (CP). RG continue to mitotically cycle while their nuclei rhythmically move dorsally up the VZ during G2/M phase (basal RG or bRG), and ventrally for/during/in S-phase along cellular projections in a process known as interkinetic nuclear migration (ventricular RG, or vRG; IKNM)13. This process continues through cortical development, with the self-replenishing pool of RG generate IPs and expand the VZ. The process of self-replenishing symmetric divisions and asymmetric neuron generation is partially regulated through the balance of key epigenetic regulators of cells, transcription factors PAX6 and EMX2, respectively14,15. IP cells, not anchored to the apical VZ, populate an outer area of the SVZ, split by an inner fiber layer (IFL), forming outer RG cells (oRG). IPs continue to divide and differentiate forming the cortex in an inside out manner, generating deep layer neurons, then the more superficial layers16. Notably, RG and IPs are known to express messenger RNA (mRNA) associated with deep and superficial layer neurons markers prior to differentiation, though they don’t express the resultant proteins. This is regulated through post-transcriptional repression mechanism and suggests a priming of RG/IPs throughout maturation17. A mature subset of IP, oRG cells form non-neuronal glial cells, such as oligodendrocytes, and astrocytes which permeate across the cortical layers18. In humans, these oRG cells are abundant and self-renew, a characteristic that has been postulated to lead to the human specific cortical expansion19–21. In recent work, differential gene expression of the transcription factor FOXO3 and genes a part of the mTOR pathway have been implicated in oRG formation and self-renewal22,23. After cortical layer formation, RG eventually self-consume into pairs of neurons.
To fully dissect the epigenomic dynamics of corticogenesis, a robust model system is needed. A model system must be both faithful to the subject of study, as well as mutable. Major considerations persist in our ability to understand cortical development both in terms of what may go awry in neurodevelopmental disorders, and what leads to the human-specific expansion of the cortex24. However, the necessary reductionist study to uncover the epigenomic landscape responsible for cortical layer stratification faces major hurdles. First and foremost is sample rarity; human and non-human primate fetal tissue is difficult to obtain. Mouse models lack several key cortical sub-regions and cell types, including the more elaborate organization of progenitors — namely the OSVZ and the oRG found within18. In addition, the developing human cortical plate and subplate, containing CR cells, have distinct cell subtypes missing in mouse1,2. Regions of accelerated mutation since human divergence from chimpanzees reveal the importance of non-coding regions. 92% of human accelerated regions (HARs; 663/721 HARs) fall outside of transcribed sites, and are enriched for enhancer-like activity or transcription-factor binding motifs25. Further, these sites are seen to be active in early embryonic forebrain development25,26. Secondly there is also the need for genetic manipulation. Necessity and sufficiency are largely determined through gene knockout and rescue experiments — corticogenesis is no different. An emerging model system must allow for both genetic manipulation and recapitulate human-specific aspects of development.
Cortical organoids are self-organizing three-dimensional cultures that model features of the developing human cerebral cortex27. They are an adaptation of a 2D cortical “rosette” method that modelled early polarization of neuroepithelial cells and neural tube formation. Induction of human embryonic stem cells (hESC), or induced pluripotent stem cells (iPSCs) to the ectodermal lineage generates cellular aggregates called embryoid bodies (EBs). Neuroectodermal lineage priming is done through in vitro differentiation of stem cells in decreased basic fibroblast growth factor (bFGF) and a high dose of ROCK inhibitors to limit cell death28,29. From here the protocol deviates from the 2D cortical rosette method to allow for three-dimensional cortical layering. EBs aggregate and are cultured in suspension in a Neurobasal medium with additives to support neural progenitors and their progeny. Shortly thereafter, EBs are embedded into matrigel, an artificial extracellular matrix, which acts as a scaffold for cell migration. EBs expand in the matrigel to form organoids containing fluid-filled cavities reminiscent of brain ventricles, and buds of neuroepithelium that replicate early to mid-gestation of cortical development. Cortical organoids do not form blood vessels and thus as they expand to up to 4 mm in diameter, the diffusion of oxygen and nutrients to the core decreases. This leads to necrotic centers if grown in culture for multiple months. To mitigate necrosis, cortical organoid protocols all feature a form of agitation to facilitate movement of nutrient rich media through the organoids. To achieve this agitation, organoids are cultured in spinning bioreactors30, or on orbital shakers, and previous groups have reportedly maintained organoids in culture for excess of 18 months31. The original organoid protocol did not include cortical region specification and was largely undirected, showing arealization of both the forebrain (FOXG1+), mid brain (OTX1/2+), ventral forebrain (NKX2.1+) and even retinal tissue. Developments in cortical organoids differentiation have revealed that they can be selectively induced to form different brain regions based on small molecule addition to the culture media. For instance, SMAD inhibitors such as dorsomorphin and SB-4321542 induce rapid neural differentiation to the dorsal forebrain state, while retinoic acid presence in early organoid induction is caudalizing30.
Cortical organoids develop in a shorter time frame than native corticogenesis occurs, yet they follow the same cell differentiation progression. Exact times vary by protocol, however a generalized timing is as follows. Within 15-20 days in vitro (DIV15-20, where DIV0 is the original induction of stem cells to ectodermal lineage), cells form continuous neuroepithelia, surrounding fluid-filled cavities (similar to neural rosettes). Pluripotency markers OCT4 and NANOG begin to diminish, while neural identity markers SOX1 and PAX6 increase32,33. By DIV30, a radially organized CP begins to form. This region is positive to pre-plate marker TBR1, and contains RELN expressing CR cells32. Bulk RNA analysis shows at this point that organoids closely resemble prefrontal cortex at GW (gestational week) 8-930. Around DIV60-75, organoids exhibit rudimentary separation between early-born deep layer corticofugal neurons (BCL11B+) and late-born superficial layer (SATB2+, POU3F2+), depending on protocol27,30,31,34. Additionally, it is around this time that the human-specific oRG cells (SOX2+,HOPX+) begin to populate27. Around DIV90-100, organoids become more closely correlated to fetal prefrontal cortex at GW17-2530. Progeny of oRG begin to form astrocytes (GFAP+)30. As organoids age further we begin to see the formation of dopaminergic neurons (TH+) and mature astrocytes (DIV180)35. Organoids in directed protocols that lack a ventral (NKX2.1+ region) did not exhibit interneuron formation27,30,31. However those with the ventral marker showed late formation of interneurons. This is expected, given that interneurons are known to migrate from the ventrally located LGE/MGE during corticogenesis. When cerebral organoids are generated from mouse embryonic stem cells, they lack both the IFL and oRG31. Further supporting evidence that organoid models can recapitulate RG behavior, is a study in which GFP was electroporated, followed by a pulse of BrdU to track lineage divisions in proliferating cells. The authors reported that daughter cells formed after the BrdU pulse chase included both RG and IPs, suggesting this cerebral organoids capture the asymmetric division potential of RG 18.
Organoid cortical models are not without limitations. In RNA comparisons between organoids and primary fetal tissue samples, organoids consistently enrich for genes associated with cellular stress, glycolysis, and electron transport pathways22,24,27. However, it has been demonstrated that this can be alleviated with culturing alterations and is likely induced in the early stages of ectodermal lineage priming of pluripotent stem cells. Organoids transplanted into a mouse cortex appropriated mouse oligodendrocytes and astrocytes, which led to decreased cellular stress signals27. New protocols have introduced vascularization processes to address glycolysis concerns as well36. Xenografted organoids show a higher correlation of radial glia maturation to primary sample age over organoid age22. Further the directed differentiation of organoids is not perfect. Mesodermal linage cells have been uncovered, despite early patterning to the neuroectodermal fate35 and organoids are not homogenous in forebrain cortical area, with many showing both primary visual cortex (V1-like) and prefrontal cortex (PFC-like) signatures22,27. This is partially to be expected given the belief that thalamic input helps define areal signature37,38. Organoids tend to lack the diversity of cell subtypes that form over time in the human cortex27. Despite nuances in organoid differentiation when compared to native human corticogenesis, this model system remains extremely promising for a battery of previously untestable hypotheses and closely resembles early corticogenesis.