Efforts to model OUD in vitro began with the generation of iPSC-derived 2D neuronal cultures from individuals either with opioid dependence or carrying genetic variants linked with increased opioid addiction risk. These studies were a response to the lack of patient- and gene-specific research into processes that bring about vulnerability to opioid dependence. The neuronal cultures themselves were patterned to represent cell types relevant to the neuropathology of OUD (Table 1).

In the first such investigation of its kind, Sheng et al. (2016a) generated iPSC-derived midbrain dopaminergic (DA) neurons from opioid-dependent subjects, motivated by the DA system’s association with reward and addiction. Expanding the use of this culture system in a parallel study, they also derived DA neurons from opioid-dependent individuals carrying variable number tandem repeat (VNTR) polymorphisms in the human dopamine transporter (hDAT) gene (Sheng et al., 2016b) associated with substance misuse (Heinz and Goldman, 2000). Relative to non-dependent controls, DA neurons from opioid-dependent subjects in both studies exhibited reduced expression of the dopamine D₂ receptor (Drd2). In addition, Sheng et al. (2016b) identified that increased VNTR length corresponded to lower DAT transcript levels, implying a role for this polymorphism in the regulation of hDAT gene expression. Interestingly, both Drd2 and hDAT expression levels were rescued by treatment with valproic acid (VPA), an anti-epileptic drug implicated in relapse prevention (Romão et al., 2022). Overall, the fidelity of these results to known DA pathway disruptions in OUD (Koob and Volkow, 2016; Burns et al., 2019) and prior neuroimaging studies of OUD patients (Wang et al., 1997; Volkow et al., 2004), reinforced the utility of opioid-dependent subject derived neurons for further studies of opioid dependence and treatment. These initial experiments by Sheng et al. (2016a,b) also highlighted the genetic tractability of iPSC-derived neuronal systems for the study of OUD, opening avenues for further examinations of underlying molecular dynamics through the modification or correction of disease-causing mutations (e.g., via CRISPR).

The use of the paradigm established by Sheng et al. (2016a,b) was next expanded by Halikere et al. (2020) to understand how molecular disruptions upstream of the DA system might enhance susceptibility to OUD. Evidence that DA neurons are excited through the suppression of inhibitory neurons following μ-opioid receptor (MOR) activation, prompted the team to explore the cellular repercussions of disrupting this pathway. To this end, Halikere et al. (2020) derived inhibitory neurons (iN) from individuals carrying the addiction-risk associated A118G single nucleotide polymorphism (SNP) in MOR. In these iNs, MOR activation by the μ-opioids DAMGO or morphine resulted in heightened inhibition, which manifested as reductions in synaptic release. This suppression of iN activity, in turn, implied increases in downstream DA activation and release, as occurs during acute opioid intoxication (Koob and Volkow, 2016; Uhl et al., 2019). Together, these findings constituted novel conceptual progress in understanding the etiology of opioid dependence at a cellular level and provided further rationale for generating OUD-relevant cell types in vitro.

While Sheng et al. (2016a,b) and Halikere et al. (2020) focused on understanding the role of gene variants in the pathogenesis of OUD, recent in vitro studies have shifted to modeling phases of the opioid addiction cycle (i.e., binge, withdrawal, and anticipation) and/or outcomes like overdose. In a cross-sectional study of heroin-dependent patients undergoing drug detoxification, Chen et al. (2022) used iPSC-derived neurons to ectopically express exosomal miRNAs identified in their blood during various stages of opioid withdrawal. The circulating miRNAs served as biomarkers for OUD progression and also influenced transcriptional programs associated with neurotransmitter dynamics, neurite outgrowth, and neural growth at a cellular level (Chen et al., 2022). The following year, Guo et al. (2023), developed a model of opioid overdose by generating iPSC-derived neurons representing the preBötzinger Complex (preBötC), a brainstem structure necessary for inspiratory rhythm generation, which is suppressed by opioids. These neurons exhibited dose-dependent cessations in activity due to four μ-opioids (fentanyl, codeine, DAMGO, and methadone) and recovery upon naloxone administration (Guo et al., 2023). Although the cells in both studies were not derived from OUD patients, they helped to demonstrate the utility of 2D neuronal cultures for parsing the cellular and molecular changes associated with specific phases of the addiction cycle, as well as identifying valuable biomarkers and therapeutic targets for its sequelae.

Altogether, such rapid developments in the generation of iPSC-derived neuronal models of OUD provoked questions regarding their fidelity to in vivo signatures of the disorder. To address this, Mendez et al. (2023), engineered novel iPSC-derived cortical neurons from skin fibroblasts of individuals who had died of an opioid overdose. Following chronic treatment with morphine, these neurons remarkably showed transcriptional alterations paralleling those observed in the postmortem, ex vivo frontal cortex tissue of individuals with OUD (Mendez et al., 2023). These included developmental and synaptic genes associated with substance use disorders (Gallo et al., 2018; Seney et al., 2021), as well G-protein-coupled receptor (GPCR) pathways, of interest given that opioid receptors are GPCRs themselves. While there are caveats related to extrapolating disease signatures from postmortem samples (e.g., cellular and molecular deterioration), these findings helped configure an informed, preliminary picture of how efficacious iPSC-derived neuronal cultures can be for recapitulating key molecular features of OUD.

Although Mendez et al. (2023) assembled a strong case for the utility of iPSC-derived neuronal cultures in modeling and studying OUD, their study also highlighted pitfalls associated with their simplicity. Due to their two-dimensional growth patterns, these neurons lack the requisite interactions between heterogenous cell types, multi-dimensional cell–cell contact and communication, and nutrient/oxygen diffusion that confer relevance to in vivo neurophysiology. This lack of tissue complexity and organization impacts cellular growth, development, and survival, which complicates interpretations of disease mechanisms (Centeno et al., 2018; Mendez et al., 2023; Mendez and Walss-Bass, 2024). It was these technical gaps in 2D neuronal culture that spurred attempts to recapitulate OUD via 3D brain organoids or spheroids (Table 1).

The inaugural steps toward this objective were taken by studies testing the efficacy of 3D neural spheroids as a high-throughput screening (HTS) platform for compounds intended to model, diagnose, or treat OUD. Despite their limited structural organization compared to organoids, spheroids (self-assembled spheres of different neural cell-types) were chosen for HTS due to their shorter incubation times and relatively higher homogeneity. Boutin et al. (2022) used cortical spheroids to test a library of neuroactive compounds targeting opioid receptors or psychoactive compounds linked with depression, anxiety, and analgesia, which are sequelae of long-term opioid misuse. Following drug exposure, activity changes in these spheroids, represented by fluctuations in calcium fluorescence, were measured using a fluorescent imaging plate reader (FLIPR). MOR agonists were found to have an inhibitory effect, reducing the count and increasing the spacing of calcium activity peaks (Boutin et al., 2022). The consistency of this response with the MOR-activation-induced suppression of cortical neuron activity and synaptic loss in animal studies, demonstrated the potential of this culture system for modeling OUD (Chang et al., 1997; Robinson and Kolb, 1999).

This prospect enabled Strong et al. (2023) to expand the utility of this methodology beyond drug screening to disease modeling. Their team generated novel iPSC-derived neural spheroids mimicking the prefrontal cortex (PFC) and ventral tegmental area (VTA), key regions involved in opioid addiction. Importantly, these spheroids reproducibly retained cell-type compositions that conferred both physiological relevance and region specificity. Considering the vital role of neuronal-glial interactions in neural communication, all spheroids were generated using 90% neurons and 10% astrocytes. Neuronal subtypes in PFC and VTA spheroids were also included in ratios that corresponded to postmortem examinations of the human brain and resulted in unique calcium activity phenotypes. Using this system, Strong et al. (2023) were able to model regional responses to both the intoxication and withdrawal phases of OUD via chronic treatment with and deprivation of the MOR agonist DAMGO. During chronic treatment, PFC-like spheroids experienced reductions in calcium activity peak counts, while treatment and withdrawal both increased peak count in VTA-like spheroids. Although the PFC deficits were rescued by naloxone, the same was not true for the VTA spheroids, indicating fundamental differences in recovery from opioid exposure between brain regions (Strong et al., 2023). This study introduced the first intentional iPSC-derived 3D model of OUD in vitro; its value reinforced by the mechanistic insights it provided into region-specific responses to chronic opioid exposure. Moreover, Strong et al. (2023) contributed technical advancements that will prove valuable for future studies of OUD in vitro. These included the successful incorporation of genetically encoded biosensors for continuous neuronal activity monitoring in spheroids, and the fusion of VTA- and PFC-like spheroids into assembloids with functional neural circuitry that can be altered by designer drugs.

Contemporary with the development of the neural spheroid model of OUD, Ho et al. (2022) generated the first 3D organoid model of this disorder. Given the PFC’s role in drug reward, withdrawal and relapse during addiction, the group generated iPSC-derived forebrain organoids from individuals with OUD. Subsequently, they used this model to examine mechanisms of differential drug action in opioid-dependent subjects at a single-cell level. Focusing specifically on oxycodone and buprenorphine, two of the most prescribed opioids in the United States, Ho et al. (2022) conducted single-nucleus RNA-sequencing (snRNA-seq) and found that both drugs alter the expression of distinct genes and molecular pathways. While buprenorphine selectively influenced transcriptional regulation in glia, oxycodone activated immune-response associated signaling (STAT1 and type 1 interferon) across several neural cell types in OUD-derived forebrain organoids (Ho et al., 2022). Not only did this research establish a brain organoid model of OUD pathophysiology, but it also established a preliminary repository of drug- and cell-type specific molecular changes associated with opioid-exposure in dependent subjects that may be used for further mechanistic probing or therapeutic development.

Although limited, the technical and conceptual progress made with regard to modeling OUD via 3D organoid and spheroid cultures has been promising. These initial studies have established a strong framework upon which further innovations and mechanistic investigations in the field of addiction research may be conducted in vitro. This is especially vital when it comes to maternal OUD, which has not yet been modeled or explored using iPSC-derived neurons, organoids, or spheroids. Prior evidence suggests that drug metabolism and pharmacokinetics are significantly altered by pregnancy and can vary considerably between individuals (Farid et al., 2008; Costantine, 2014; Feghali et al., 2015), which makes the development of patient-, tissue-, gene-, and cell type-specific in vitro models even more critical. As these methodologies continue to evolve, their application toward understanding cellular and molecular mechanisms of opioid dependence and addiction during pregnancy will be of prime importance for the improvement of maternal and fetal health, and the identification of novel clinical interventions.

In contrast to in vitro models of OUD, iPSC-derived 2D neuronal cultures used to study the effects of opioid exposure were developed contemporaneously with organoid models of POE. Given this, they cannot strictly be considered foundational for the development of more complex 3D culture systems in this field. Additionally, the usual aim of these studies was not to model POE, but to engineer neural cell types relevant for screening therapeutics that confer neuroprotection or non-opioid based analgesia (Table 2). Therefore, we will only briefly touch upon their findings, focusing instead on the technical advancements that make these cultures pertinent to investigations of POE in vitro.

The first neuronal cultures relevant to the study of POE were developed by Ju et al. (2021), who used iPSCs to generate neurons expressing μ-(MOR) and κ-opioid receptors (KOR). Opioids bind to three major opioid receptors throughout the central and peripheral nervous systems: mu-(MOR), kappa (KOR), and delta (DOR). Early preclinical studies of nervous system opioid pharmacodynamics revealed the broad distribution these receptors throughout the CNS, and highlighted differences in receptor expression between fetal/neonatal and adult brains (Barg and Simantov, 1989; Rius et al., 1991; Wittert et al., 1996; Zhu et al., 1998). MORs and KORs are the first opioid-receptors to appear in the fetal brain, while DORs appear postnatally (Rius et al., 1991; Zhu et al., 1998). This finding highlights the possibility that opioids may have distinct effects depending on developmental stage, making the MOR and KOR expressing neurons generated by Ju et al. (2021) an exceedingly relevant model system for POE. In addition, because these neurons originate from exfoliated renal epithelial cells in urine, they afford a unique level of scalability based on the ready supply of source material that can be clinically and non-invasively obtained from maternal OUD subjects or opioid-exposed neonates (Ju et al., 2021).

One problem, however, was that the neurons generated by Ju et al. (2021) did not possess any regional or subtype identity. Consequently, recent efforts have focused on generating cell types more specific to mechanisms of opioid action in vivo. Given the broad use of opioids for pain management, Deng et al. (2023), Nimbalkar et al. (2023), and Röderer et al. (2023) worked on generating and modifying iPSC-derived sensory nociceptive neurons as an experimental platform for screening alternative analgesics. Nociception is the process of communicating electrical impulses generated by noxious stimuli. It is important to note that while fetal nociceptive pathways are thought to develop by as early as 7–10 weeks of gestation, the inception of pain sensation or perception remains a controversial topic (Thill, 2022). Therefore, conservatively, these cultures provide a unique opportunity to explore the onset and mechanisms of fetal nociceptive responsivity to opioids. Specifically, the utility of such iPSC-derived sensory nociceptors is conferred by their expression of opioid-receptors (MOR, KOR, DOR, as well as the nociception opioid peptide receptor, NOP) and activity suppression upon μ-opioid exposure (Nimbalkar et al., 2023; Röderer et al., 2023). However, the timing of opioid receptor expression and opioid responsivity varied between protocols. In Nimbalkar et al. (2023), only MOR and KOR were expressed by day 21 in culture, but not DOR. Meanwhile, all opioid receptors were expressed after 21 days of differentiation in the Röderer et al. (2023) study, although signaling through these receptors was not noted until day 70. These differences highlight the caveat of variability that may arise with the use of iPSC-derived in vitro models for the study of POE (Volpato and Webber, 2020; Beekhuis-Hoekstra et al., 2021; Brunner et al., 2023; Nath et al., 2023).

Nevertheless, the integration of iPSC-derived nociceptors with multi-electrode arrays (MEA) (Nimbalkar et al., 2023) and the generation of peptidergic and non-peptidergic sensory neurons (Deng et al., 2023) in these studies have expanded the utility of this model for studying opioid effects on electrophysiology and cellular subtypes in the fetal brain. The longitudinal investigation of iPSC-derived nociceptor maturity by Röderer et al. (2023) also supplied helpful information regarding windows of opioid-responsivity in these cultures, with fentanyl only inhibiting the activity of protein kinase A-II (required for pain sensitization) after 70 days of differentiation. Therefore, even though such 2D neuronal cultures have not yet provided insights into the neurodevelopmental effects of POE, they remain advantageous, well-studied platforms upon which future studies may be established.

With respect to the development of 3D organoid models of POE, progress has primarily centered around studying the neurodevelopmental impact of opioid-based pharmacotherapies clinically recommended for the Medication-Assisted Treatment (MAT) of maternal OUD (Substance Abuse and Mental Health Services Administration, 2016; The American College of Obstetricians and Gynecologists, 2017) (Table 2). This is due to the rising rates of pregnant women seeking treatment for OUD (Martin et al., 2015; Krans et al., 2019), motivated by the need to improve their own health and prevent neonatal opioid withdrawal (Cleveland and Bonugli, 2014; Frazer et al., 2019; Macfie et al., 2020).

Since the 1970s, MAT using methadone, a synthetic opioid analgesic and full μ-opioid receptor agonist, has been primary standard of care for opioid-addiction during pregnancy (Payte, 1991; Center for Substance Abuse Treatment, 2005; Farid et al., 2008; Krans et al., 2019). However, evidence of methadone’s ability to readily cross the placenta and accumulate in animal and human fetal tissues (Farid et al., 2008; de Castro et al., 2011; Kongstorp et al., 2019; Badhan and Gittins, 2021), in addition to its association with long-term neurocognitive deficits (Wong et al., 2014; Bier et al., 2015; Chen et al., 2015; Hauser and Knapp, 2018; Grecco et al., 2021; Levine et al., 2021; Lum et al., 2021), have led to concerns regarding its effects on fetal neural development in utero. These apprehensions have been compounded by methadone’s tendency to cause Neonatal Abstinence Syndrome (NAS), a collection of symptoms associated with withdrawal from POE, which yields CNS hyperirritability and autonomic nervous system dysfunction (Jones et al., 2010; Gaalema et al., 2012).

It was these contraindications that prompted Wu et al. (2020) and Yao et al. (2020) to generate the first iPSC-derived organoid models of POE. These studies integrate human iPSC derived-cortical organoid (hCO) cultures (Trujillo et al., 2019) with immunofluorescence, MEA, or patch-clamp electrophysiology techniques to probe how methadone alters neural growth and function in the embryonic brain. Yao et al. (2020) observed that methadone dose- and timeline-dependently alters the growth of hCOs, while also having a significant effect on neuronal and neural network function. Methadone suppressed the firing of spontaneous action potentials by hCOs attached to MEA plates, which the group hypothesized was likely due to the drug’s concurrent reduction of synaptic transmission (i.e., diminished frequency and amplitude of spontaneous excitatory post-synaptic currents) and voltage-dependent sodium currents that support the initiation of action potential burst firing (Yao et al., 2020).

While Yao et al., 2020 examined the first 3-month of hCO culture, Wu et al. (2020) extended this timeline to track the electrophysiological consequences of methadone exposure in 3–6 month-old hCOs, a period corresponding to neuronal and network activity maturation in utero (Trujillo et al., 2018, 2019). They uncovered that 12-weeks of chronic exposure to methadone suppresses the maturation of neuronal membrane properties and excitability via the impairment of voltage-dependent ion channel functions (Wu et al., 2020). Combined with Yao et al. (2020)’s results, these findings provided strong evidence that prenatal methadone exposure causes delays in the onset and progression of neural maturation in the fetal cortex. A subsequent study conducted by Dwivedi et al. (2023) contributed further proof of this effect. Bulk mRNA-sequencing of 2-month-old hCOs that had been chronically treated with methadone for 50 days yielded a robust transcriptional response, pointing toward interrelated alterations in functional components of the synapse, underlying extracellular matrix (ECM), and cilia. Methadone’s impact on molecular processes of synaptic assembly and activity during synaptogenesis in hCOs reinforced the drug’s deleterious influence on neuronal communication and, therefore, maturation of cortical functions (Dwivedi et al., 2023).

Taken together, Wu et al. (2020) and Yao et al. (2020) constituted the first proof-of-concept studies for using brain organoids to study the neurodevelopmental effects of opioids. Alongside Dwivedi et al. (2023), the findings from these investigations provided valuable insights into the structural and functional impact of methadone on fetal cortico-genesis. More specifically, they also supplied the first cellular- and molecular-evidence that methadone impacts synaptogenesis and synapse biology in the human fetal brain. The results from all three papers have helped to initiate a broad picture of how prenatal methadone exposure may give rise to long-term neurologic deficits.

Alongside methadone, buprenorphine is another widely employed opioid-based pharmacotherapy for the treatment of maternal OUD (Substance Abuse and Mental Health Services Administration, 2016; The American College of Obstetricians and Gynecologists, 2017). The use of this drug during pregnancy has become progressively more common, in part due to its inherent pharmacology (Zedler et al., 2016; Krans et al., 2016a). Unlike methadone, buprenorphine’s nature as a partial MOR and NOP receptor agonist and KOR antagonist with low intrinsic activity means that it can be administered by outpatients with lower overdose risk and fewer drug interactions (Suarez et al., 2022). Moreover, several randomized controlled trials have demonstrated that buprenorphine yields better neonatal outcomes than methadone (Fischer et al., 2006; Kakko et al., 2008; Lacroix et al., 2011; Metz et al., 2011), including fewer signs of NAS and less time or morphine required to treat the syndrome (Jones et al., 2010). Studies in both animals and humans have also shown that prenatal exposure to buprenorphine generates superior neurocognitive outcomes, birth weights, head circumferences, and risks of preterm birth than methadone (Coyle et al., 2012; Zedler et al., 2016; Kongstorp et al., 2019, 2020; Kinsella et al., 2022; Suarez et al., 2022; Gottlieb et al., 2023). Despite these advantages, however, buprenorphine readily crosses the placental barrier (Nanovskaya et al., 2002) and has been linked to adverse postnatal behavioral sequelae (Hung et al., 2013; Sundelin Wahlsten and Sarman, 2013; Tobon et al., 2019), cellular-level alterations in neurogenesis (Pettit et al., 2012; Wu et al., 2014), and deficits in myelination (Sanchez et al., 2008; Eschenroeder et al., 2012). It is these contradictory consequences associated with prenatal buprenorphine exposure that have driven recent efforts to mechanistically dissect its neurodevelopmental effects using 3D organoids.

In 2022, Nieto-Estévez et al. took the first step in this regard by exposing iPSC-derived human cortical spheroids (hCS) and subpallial spheroids (hSS) (expressing markers of developing excitatory and inhibitory interneurons, respectively) to buprenorphine. Their study aimed to dissect the effects of buprenorphine on the crucial yet precarious excitation/inhibition balance that belies cortical network activity (Nieto-Estévez et al., 2022). It was also the first investigation of its kind to use assembloids (Bagley et al., 2017; Birey et al., 2017), fusions of region-specific organoids or spheroids, in order to investigate POE in vitro. Upon fusing the hCS and hSS, Nieto-Estévez et al. (2022) observed increased inhibitory interneuron migration from subpallial to cortical spheroids as well as an increase in network activity in response to chronic buprenorphine treatment. Although this latter result seems contradictory, evidence suggests that the inhibitory neurotransmitter GABA exerts an excitatory influence during embryonic development (Leinekugel et al., 1999) that may impact synapse formation and function (Wolf et al., 1986; Caillard et al., 1999a,b). Taken together, these findings suggests that buprenorphine influences both the development, spatial organization, and activity of inhibitory neurons in the cortex.

Interestingly and in contrast to the iPSC-derived neuronal and cortical organoid models of POE cited above, the hCS or hSS generated by Nieto-Estévez et al. (2022) did not express the major opioid receptor subtypes MOR, KOR, or DOR. Buprenorphine was instead found to bind and signal through the NOP receptor, an opioid G-protein coupled receptor expressed throughout the human fetal cortex that does not respond to opioids with known misuse liability (Neal et al., 2002; Zaveri, 2016). Dysregulation of NOP signaling has been linked to psychiatric disease, depression, and memory deficits (Wang et al., 2009; Post et al., 2016; Khan et al., 2018), all of which are sequelae associated with POE. While this feature enabled the team to study buprenorphine’s effects on the fetal brain via the NOP receptor, it limits the use of this model in future studies intending to investigate the drug’s action through canonical opioid receptors.

Unlike Nieto-Estévez et al. (2022), the iPSC-derived cerebral organoids (CeO) generated by Fernandes et al. (2022) were confirmed to express MOR, DOR, and KOR within 30 days of differentiation in both neurons and glial cells. The expression of opioid receptors on glial cells in this model was unique, given discrepancies in findings regarding the presence of these receptors on astrocytes in vivo (Stiene-Martin et al., 2001; Machelska and Celik, 2020). Moreover, this was first use of a region non-specific organoid to study POE in vitro (Fernandes et al., 2022). Using this model, Fernandes et al. (2022) found that modulating opioid receptor activity using buprenorphine increased apoptosis, astrogliogenesis, glial cell maturation, and dopamine release in CeO cultures, consequences of chronic opioid exposure that have been observed in prior studies. In parallel, the team also subjected their CeOs to the bone marrow stem cell secretome (BMSCSec), which has been explored as an antinociceptive treatment alternative to opioid-based analgesia (Brini et al., 2017; Gama et al., 2018; Khatab et al., 2018). Interestingly, the BMSCSec almost exactly mirrored buprenorphine’s effects, apart from increasing dopamine release. On top of highlighting the effects of buprenorphine on glia during development, the use of CeOs and the stem cell secretome in this study also constituted a technical advancement for POE research. While CeOs provide a novel platform to investigate opioid effects in neurons and glia across developing brain regions (Lancaster and Knoblich, 2014a), the secretome introduces a new method of modulating opioid receptor activity that can be used to individuate the prenatal effects of opioids.

This aim of delineating the effects opioid-based pharmacotherapies was carried forward by Yao et al. (2023), who used iPSC-derived cortical organoids to identify how buprenorphine and methadone differentially affect cellular growth and neuronal activity in the developing cortex. In this study, buprenorphine was found to have a milder effect than methadone on neural growth and activity in hCOs. Although 5–10 times less buprenorphine is required to achieve withdrawal relief than methadone, even at equivalent concentrations buprenorphine did not suppress neural network action potential firing rates. In fact, pre-treatment of hCOs with buprenorphine consistently blocked the severe growth suppressive effects of methadone and the drug even enhanced growth at higher (10 μM) concentrations. Yao et al. (2023) proposed that these distinct effects of methadone and buprenorphine on growth and neural activity are due to their contrasting activity at κ-opioid and NMDA receptors, respectively. Buprenorphine’s antagonism of KOR activity, which has been implicated in cell proliferation, differentiation, and death, as well as its lack of NMDA receptor antagonism were proposed to be the underlying cause of its tempered influence on hCO growth and function (Yao et al., 2023). Cumulatively, these results bring to light novel mechanistic details that may help to explain buprenorphine’s long-supposed superiority when it comes to neonatal outcomes.

As the above-cited articles demonstrate, investigations into prenatal opioid exposure have predominantly been geared toward exploring the effects of opioids used to treat OUD. However, two recent studies exemplify a potential shift in focus toward other opioids as well (Table 2). In a broad exploration of how narcotic and neuropsychiatric-risk factors influence neurodevelopment, Notaras et al. (2021) exposed iPSC-derived forebrain organoids to a panel of “enviromimetic” chemicals and analyzed downstream alterations in transcription, proteomics, and metabolomics. Exposure to opioids was modeled using the endogenous MOR agonist endomorphin, which is central to nervous system pain relief and stress response pathways (Fichna et al., 2007). Interestingly, endomorphin elicited proteomic changes associated with axon guidance, cellular stress response, and RNA regulation, paralleling the effects of cannabinoids, nicotine, and ethanol. Metabolic analyses also revealed converging effects of all treatments on L-Phenylalanine and GTP expression, pointing toward increases in systemic stress (Fernandes et al., 2010) and disruptions in energetics (Leone et al., 2010; Montenegro-Venegas et al., 2010) during corticogenesis. Altogether, this study contributed to the growing body of knowledge surrounding the influence of MOR activation on normative cortical development. Importantly, it also opened a door for future comparative and/or synchronous in vitro explorations of opioids and other neurodevelopmentally noxious factors (Notaras et al., 2021).

Nevertheless, endomorphin is not an exogenous opioid, and there are no indications that the neuropeptide’s mimetics are misused during pregnancy. This gap was remedied in a contemporaneous study by Kim et al. (2021), who conducted single-cell RNA-sequencing of iPSC-derived midbrain organoids exposed to fentanyl. Fentanyl is a potent opioid analgesic that is prescribed for the management of severe pain in preterm neonates (Hall and Shbarou, 2009) and during pregnancy, most often during labor (Shoorab et al., 2013). However, fentanyl and its derivatives retain a high degree of misuse potential. As a stark reminder of this fact, fentanyl accounted for 39% of drug overdose deaths in the United States in 2017 (Hedegaard et al., 2019). Therefore, Kim et al. (2021)’s novel examination of fentanyl’s effects on the human fetal midbrain was particularly timely. Corroborating the opioid-induced dysregulation of midbrain dopamine reward pathways (Wei et al., 2018), the group found that acute fentanyl exposure increased dopamine release in the organoids. In contrast, chronic fentanyl treatment arrested the fate determination of neural progenitor cells and altered the expression of synaptic activity and neuronal projection pathways (Kim et al., 2021). These findings were analogous to the neurodevelopmental effects of methadone and buprenorphine reported in cortical or cerebral organoids (Wu et al., 2020; Yao et al., 2020, 2023; Fernandes et al., 2022; Nieto-Estévez et al., 2022; Dwivedi et al., 2023). All the same, this study by Kim et al. (2021) unfurled a list of novel possibilities when it came to studies of POE in organoids, especially regarding the opioid-types and brain-regions modeled.