Preprint
Article

This version is not peer-reviewed.

Effects of Sevoflurane on the Development of a Human Brain Microphysiological System

A peer-reviewed version of this preprint was published in:
International Journal of Molecular Sciences 2026, 27(7), 3322. https://doi.org/10.3390/ijms27073322

Submitted:

25 March 2026

Posted:

27 March 2026

You are already at the latest version

Abstract
Animal studies have shown that early life exposure to general anesthetics may impair brain development. However, the implications of this phenomenon in human patients remain unclear. In this study, we use an induced pluripotent stem cell (iPSC)-derived human brain microphysiological system (bMPS) to investigate the effects of early sevoflurane (SEV) exposure on human brain development. Human iPSCs were cultured and differentiated into neural progenitor cells (NPCs) and then into bMPS. At week 8, bMPSs were exposed to 2.4% SEV for 4h. Four weeks after exposure, immunofluorescence, Western blotting, and quantitative real time polymerase chain reaction (qPCR) were conducted to evaluate the alteration of nerve cells in bMPS. After SEV exposure, number of apoptotic cells increases and the level of neural differentiation marker decreases. The ratio of mature neurons over NPCs and mature oligodendrocytes over oligodendrocyte progenitor cells (OPCs) are reduced which leads to reduction of myelination. SEV also impedes the development of astrocytes and synaptogenesis, especially the formation of excitatory synapses. Meanwhile, SEV increases the expression of molecules in mammalian target of rapamycin (mTOR) signal pathway. In conclusion, early SEV exposure substantially disrupts the development of human brain tissue. The mTOR signal pathway is involved in this alteration.
Keywords: 
;  ;  ;  ;  ;  

1. Introduction

Over two hundred million general anesthesia are performed each year to allow patients to comfortably tolerate surgeries and procedures which would otherwise be impossible, [1] and it is widely accepted that the anesthetic itself has only transient effects on the patient. However, there is concern that some vulnerable categories of patients, especially young children and individuals with underlying brain disorders, might be at risk of lasting neurologic impairments due to toxicity effects caused by exposure to general anesthetic (GA) drugs on vulnerable elements on the brain during critical periods of brain vulnerability [2,3,4]. There is evidence from early human studies in this area that surgery and exposure to GA from late postnatal through early childhood periods may lead to cognitive impairment as measured by an increased incidence of learning disorders and worsening performance on school assessments [5,6]. More recent work has focused on the possibility that early life surgery and GA exposure may correlate with a variety of adverse behavior outcomes [7,8,9]. Based on these studies, the U.S. Food and Drug Administration issued a warning that lengthy or repeated exposure to general anesthetics and sedative drugs from the third trimester of prenatal development through the first 3 years of life may cause disruptions of brain development [10]. Human studies invariably have substantial and unavoidable confounds, including surgery and co-morbid disease, and thus much of the relevant literature that focuses specifically on the contribution of anesthetics has been conducted in animal models.
In rodent and non-human primate animal studies, there is little question that early life exposure to GA agents substantially interferes with brain development and causes lasting neurocognitive and neurobehavioral sequelae [11,12,13,14,15,16,17,18,19]. The preponderance of this literature, and nearly all of our knowledge of putative mechanisms of developmental anesthetic neurotoxicity, is in rodent models. It has been challenging to ascertain how these findings would be expected to translate to human patients, given the substantial interspecies differences in brain structure and developmental timeline. In this study, we aim to begin bridging this gap by using human brain cerebral microphysiological system (bMPS), a versatile and novel model for investigating the toxic effects of anesthetics on the developing brain [20]. The bMPS is a three-dimensional (3-D) culture system derived from induced human pluripotent stem cells (iPSCs). It faithfully replicates the key human brain cell types and generates the interactions between cells that are critical for brain development in living organisms (such as active synapses and axon myelination) [21,22]. The human bMPS represents a tremendous step forward in developmental neurotoxicology testing, as arguably it has a high degree of applicability to patients given the use of human-derived iPSC and what is presently available in the basic science armamentarium [23,24]. For this study, we have employed human bMPS which mimics the cellular microenvironment of the human cerebrum and is comprised of neurons, astrocytes, and oligodendrocytes in physiologically relevant ratios [25]. The critical elements of CNS structures and functions, such as neuron-neuron interaction (e.g., spontaneous electric field potentials), neuron-glia connectivity (e.g., myelination), and synaptogenesis, are recapitulated. At week 8, the bMPS were exposed to sevoflurane (SEV) for 4 hours. It is not possible to make an equivalence between human chronological age and bMPS days in vitro, given the vast heterogeneity in developmental timelines of the human brain in vivo. We selected 8 weeks in vitro for SEV exposure because, by this time, the brain microenvironment is established and large populations of neurons and glia have fully differentiated, while neurogenesis and glial development remain ongoing, and the putatively vulnerable processes of synaptogenesis and myelination are still in early stages [26,27]. After 4 additional weeks of growth, the bMPS were harvested and exanimated at week 12, representing a more mature stage of neural development (Figure 1). This study addresses an important translational gap by determining whether developmental effects of GA previously described in rodent models are also observed in human-derived three-dimensional brain models. Although the work does not primarily aim to define a new mechanism, it provides evidence that these effects extend to a human multicellular organoid system, which is a meaningful advance for evaluating human developmental neurotoxicity. Because cross-species concordance cannot be assumed for neurodevelopment processes, especially in the context of human brain development, confirming these effects in a human-relevant model represents more than simple replication of animal literature. Thus, the value of the present study lies in establishing human relevance and strengthening the translational foundation of this field. We tested the central hypothesis that developmental exposure to GA inhibits neural development in a human brain model. The bMPS used in this study is a valuable tool for studying neurodevelopmental neurotoxicity because it reproduces multiple key aspects of human brain development withing a single integrated platform. Specifically, it captures major neurodevelopmental processes, including proliferation, apoptosis, differentiation, migration, neurite growth, synaptogenesis, myelination, and neural network formation.
In this study, we focused on SEV, the most widely used primary anesthetic agent in early-life pediatric patients. Using quantitative immunofluorescence (IF) microscopy, Western blotting (WB), and quantitative polymerase chain reaction (qPCR), we assessed neuronal growth and differentiation, oligodendrocyte differentiation, myelination, astrocyte production, and synapse formation. We also examined the involvement of mammalian target of rapamycin (mTOR) signal pathway in this GA caused alteration in human brain.

2. Results

2.1. Effect of SEV Exposure on Gross Development

First, we investigated whether SEV exposure affects neural tissue development by measuring bMPS size. To accomplish this, we used light microscopy at the end of week 12 (4 weeks after exposure) to measure the maximal cross-sectional area for each human bMPS and to determine size. In the control (CON) group, the area is measured at 0.20 ± 0.03 mm2 (Figure 2A), and it decreases to 0.17 ± 0.03 mm2 in the SEV exposure group (Figure 2B). This represents a statistically significant decrease (n = 12; p < 0.05) (Figure 2C), and it is consistent with a lasting and substantial neurotoxic effect of early developmental anesthetic exposure to human brain tissue.

2.2. Effect of SEV Exposure on Developmental Apoptosis

Given the extensive literature suggesting that early developmental anesthetic exposure can cause an increase in apoptosis in a range of different brain cell types, [28,29] we counted the cells labeled with cleaved caspase-3, a death protease that is required for progression of apoptotic processes, in SEV-treated bMPS and CON. A mean density of 662.17 ± 159.03 per mm2 caspase-3 labeled cells is observed in the CON, and this number is significantly increased to 1,008.5 ± 164.84 per mm2 in SEV group (n = 6; p < 0.01) (Figure 2D–F). These findings suggest that SEV exposure results in a substantial upregulation of apoptotic pathways in developing human brain cells.

2.3. Effect of SEV Exposure on Neurogenesis

Numerous reports on rodent models using both in vivo and in vitro systems have shown evidence that early exposure to anesthetics can cause alterations in neurogenesis, which in turn result in structural changes and functional impairments [30]. The effects of SEV exposure on neurogenesis were evaluated by measuring the density of neural progenitors and mature neurons, and by examining their relative proportions as a measure of developmental progression. IF for the early neural progenitor marker Nestin and the mature neuron marker microtubule-associated protein 2 (MAP2) was conducted (Figure 3A,B). In general, no Nestin+/MAP2+ double-labeled cells are observed, which indicates that Nestin and MAP2 signals represent distinct cell populations as expected. Overall, the ratio of MAP2+ cells over Nestin+ cells in CON is 164.80 ± 55.87% as compared to 55.82 ± 22.67% in SEV exposure group (p < 0.01) (Figure 3C). Western blotting tests indicate the ratio of immuno-intensity of sex determining region Y-box 2 (SOX2), the neural differentiation marker, to β-actin in control (75.92 ± 12.71%) as compared to 53.93 ± 16.06% in SEV exposure condition (p < 0.05) (Figure 3D,E). The intensity of cortical neuronal molecule T-box brain1 (TBR1) in the SEV group is partially lower than that of the control (23.58 ± 7.04% vs. 31.48 ± 6.77%; p > 0.05) (Figure 3D,F). While this is not the only possible interpretation, it is suggestive of an impairment in the differentiation of neuronal precursors into mature neurons.

2.4. Effect of SEV Exposure on Oligodendrogenesis and Myelin Development

Anesthetic exposure has also been shown to disrupt the formation of myelin through the disruption of oligodendrocyte (OL) differentiation and OL precursor cell (OPC) maturation [15]. To test for effects of anesthetic on OL development, bMPSs were co-immunolabeled with the OPC marker neural/glial antigen2 (NG2) and the mature OL marker adenomatous polyposis coli (APC) (Figure 4A,B). There are numerous NG2+ and APC+ cells visualized in the bMPSs, but no double-labeled cells are observed, confirming the specificity of the markers. The ratio of APC+ cell number over NG2+ cells in the SEV exposure group (229.85 ± 106.83%) is significantly lower than in the CON conditions (486.75 ± 152.4%) (p < 0.05) (Figure 4C). These data suggest that SEV exposure results in a lasting impairment of differentiation of OPCs to mature OL phenotype.
Next, we investigated the effects of sevoflurane exposure on myelin development using quantitative IF measurement of myelin basic protein (MBP) levels (Figure 4D–I). Compared to CON (100 ± 25.86%), there is a reduced intensity of MBP immunoreactivity in the SEV group (67.82 ± 9.48%; p < 0.05) (Figure 4J). Additionally, we measured the sum length of myelinated axons, identified as processes that co-immunolabeled for axon marker neurofilament (NF) 200 and MBP (Figure 4F,I). Early exposure to SEV decreases the total length of NF200+/MBP+ processes (6.39 ± 1.8 mm vs. 9.1 ± 2.2 mm; p < 0.05) (Figure 4K). Taken together, these findings suggest that early developmental exposure to sevoflurane results in a lasting reduction in myelination in human bMPS. This is potentially, at least partially, resultant from a reduced number of mature OLs available to provide myelination.

2.5. Effect of SEV Exposure on Astrocyte Production

Previous work in rodent models has shown that astrocytes are also potentially susceptible to the toxic effects of exposure to anesthetics during development [31]. We tested cell loss of astrocytes by measuring the density of cells positive for the astrocyte marker glial fibrillary acidic protein (GFAP) in our human bMPS model. Compared to CON (Figure 5A–C), a reduction in the density of GFAP+ cells is seen in SEV exposed bMPSs (Figure 5D–F). We observed that GFAP immunoreactivity partially decreases from 100 arbitrary units ± 32.1% in control conditions to 69.4 arbitrary units ± 23.4% in SEV exposure conditions (p > 0.05) (Figure 5G) and the number of GFAP+ astrocytes in SEV group is lower than control (206.5 ± 56.2/mm2 vs. 125.3 ± 31.7/mm2 p < 0.01) (Figure 5H). These findings suggest that early developmental exposure to sevoflurane substantially disrupts astrocyte development in human brain tissue.

2.6. Effect of SEV Exposure on Synaptogenesis

Numerous in vitro and in vivo animal model studies have shown lasting aberrations in synaptic development, including synapse loss, which result from early anesthetic exposure [13,28]. To investigate the effects of SEV exposure on synaptogenesis, we employed IF to evaluate the intensity of immunolabeled puncta for synaptophysin (SYP), a presynaptic marker, in bMPSs (Figure 6A,B). SEV exposure significantly reduces the number of synaptophysin+ terminals as compared with control (6,679.0/mm2 vs. 4,298.2/mm2; p < 0.05) (Figure 6C). The vesicular glutamate transportor-1 (vGLUt1) is a representative excitatory presynaptic marker. The mRNA of vGlut1 was examined with quantitative real time polymerase chain reaction (qPCR). SEV dramatically reduces the relative expression of vGlut1 compared to CON (76.6 ± 8.7% vs. 100 ± 4.1%; p < 0.001) (Figure 6D). The Western blotting data shows the ratio of intensity of glutamate NMDA receptor-2B (GluN2B), a post-synaptic component molecule, over β-actin is partially declined in SEV than CON (106.5 ± 19.6% vs. 133.3 ± 21.1%; p > 0.05) (Figure 6E,F). These findings indicate that early developmental exposure to anesthetics causes a lasting and substantial reduction of presynaptic puncta density and synaptic marker abundance in human brain tissue.

2.7. Effect of SEV Exposure on Mammalian Target of Rapamycin (mTOR)

The expression of molecules related to the mTOR signaling pathway was evaluated by IF and Western blotting (WB). Our previous data showed that GA exposure pathologically upregulates the expression of the mTOR signal pathway in the rodent model. The phospho-S6 ribosomal protein (pS6), a critical downstream molecule, is used as a marker to evaluate the mTOR level in this study. The pS6 was predominantly distributed in the outer layer of bMPSs and located in both neurons and glial cells (Figure 7A,B). Compared to control (100 ± 30.13%), intensity of pS6 reactivity in SEV exposed bMPSs increases to 153.28 ± 34.58% with a statistical difference (p < 0.01) (Figure 7C). WB data show that SEV exposure also alters the expression and phosphorylation of mTOR-related signal molecule ribosomal protein 70 S6 kinase (p70S6K) (Figure 7D). In SEV conditions, the ratio of intensity of total p70S6K (t-p70S6K) over β-actin was partially higher than CON (89.15 ± 14.67% vs. 68.53 ± 11.53%; p > 0.05) (Figure 7E) and the ratio of phosphorylated p70S6K (p-p70S6K) over t-p70S6K in SEV group are significantly elevated than CON (76.84 ± 9.17% vs. 55.16 ± 8.91%; p < 0.01) (Figure 7F). These data indicate that SEV exposure predominantly upregulates phosphorylation of mTOR signal molecules.

3. Discussion

In this study, we report that sevoflurane exposure in human induced pluripotent stem cell (iPSC) derived bMPS causes a substantial disruption in the normal brain cell development. The human bMPS consists of a variety of cell types found in the human brain, including neurons, oligodendrocytes, and astrocytes, precursor cells associated with these cells, such as neural progenitors and OPCs, cell-type specific specializations including myelin, and synapses that are critical for neural transmission [24,25,26,27]. We realize that, in our study, the use of a single iPSC line is an important limitation We chose a 4-hour exposure to 2.4% sevoflurane using a previously established in vitro exposure system [28,29,30,31,32]. In our preliminary experiments, we tested a series of sevoflurane doses (1%; 2%; 2.4%; 4%; 8%). The 2.4% sevoflurane concentration chosen in this study is the lowest dose that causes significant development alterations in human bMPS. In addition, this dose corresponds to similar experiments from other researchers in which minimum alveolar concentration (MAC) of mice was tested and comparison of water/gas partition coefficient of sevoflurane in vitro and blood/gas partition coefficient in vivo was performed [33,34]. We did not conduct the anesthetic concentration stability over time in this study. However, according to our previous research using same exposure technique, anesthetic concentration stability over time should be ensured [35]. We realize that analysis of anesthetic concentration in the culture medium tested with HPLC is necessary in our further study.
Given that the one MAC of sevoflurane varies between 2% and 3.3% during birth and 3 years of life, [36] this represents our intent to model the effects of a single, relatively high dose, lengthy general anesthetic free from the confounds of physiological perturbation. We found evidence of gross perturbation in development, enhanced chronic apoptosis, impaired progression of neuron and oligodendrocyte development, reduced numbers of mature neurons, oligodendrocytes, and astrocytes, and a loss of synapses. While these findings do not prove that anesthetics are toxic to humans in vivo, nor do they speak to the scope of putative toxicity, they do establish that the phenomenon of developmental anesthetic neurotoxicity is not limited to animal model systems and does occur in developing human brain tissue.
Animal studies, which have been conducted on a limited basis in non-human primates, but the predominance of which are in rats and mice, have shown that general anesthetic exposure in early development causes impaired neurocognitive performance and pathological alteration of CNS structure [7,8,9,10,11,12,13,14,15,37,38]. However, while animal models recapitulate human brain development in many respects, there are many key differences, both at the cellular and systems levels, including divergent developmental timelines, distinct genetic and epigenetic properties, and substantially dissimilar degrees of structural complexity. To date, conclusive evidence of neuropathology resulting from developmental anesthetic neurotoxicity has not been established in human studies, and it is unlikely that this will occur in the near future, given that it would depend on non-invasive and harmless measurement techniques and also would require rigorous controls that are practically challenging to establish. Recent advances in stem cell culture models, such as the bMPS model we employed, provide the best currently available alternative. In this system, iPSCs simulate the process of human CNS development, including early-stage neural tube formation, neuroepithelium differentiation and regional specification, in a three-dimensional structure that establishes spontaneous neural transmission and activity. This model has been used to evaluate neurodevelopmental disorders (NDDs) as diverse as microcephaly, autism, focal epilepsy, and other neurological diseases [39,40]. Cultured cerebral organoids from patient-specific iPSCs or gene interference display a variety of microcephaly characteristics, including size reduction, decrease of neuron number, and enhanced cell death [20,41,42]. In another study of Zika virus-infected organoids, pathological alterations resembling microcephaly, including smaller bMPS size, increased cell death, reduced neural progenitor proliferation, and decreased neuronal cell-layer volume, were observed [43,44,45]. Our data showed similar pathological alterations to other NDDs tested in this model system at a gross level--the size of bMPS that were exposed to sevoflurane were smaller than controls. We found increased levels of apoptosis, reduced progression of progenitor cells, and ultimately a loss of mature neurons, oligodendrocytes, and astrocytes, suggesting that these factors account for the observed gross loss of bMPS size. Obviously, this phenomenon is not observed in humans who were exposed to general anesthesia during brain development—if microcephaly were a consequence of anesthesia exposure, this would surely have been discovered. However, unlike many other NDDs that have been studied with bMPS, anesthetic exposure is a transient toxic stimulus that presumably alters the developmental trajectory of a cohort of vulnerable cells, rather than a continuously present insult that interferes with most or all cells in the brain, and hence a subtler phenotype would be expected in vivo.
Although most investigations of developmental anesthetic neurotoxicity have focused on neurons, [5,12] brain function and development dependent heavily on myelin-forming oligodendrocytes, which undergo critical developmental events during the putative window of vulnerability to anesthetic toxicity. Myelination is essential in establishing connectivity in the growing brain by facilitating rapid and synchronized information transfer across the nervous system. Myelin, rather than simply providing insulation to facilitate signal transduction in the developed brain, is critically important for the development of brain circuitry [46]. Myelination involves proliferation and differentiation of oligodendrocyte progenitor cells (OPCs), maturation of oligodendrocytes (OLs), and ensheathment of axons. Damage or loss of myelin function is associated with many neurologic and psychiatric disorders, including multiple sclerosis, Alzheimer’s disease, Parkinson’s disease, and CNS injury [47,48,49,50]. Currently, most of our understanding of neurotoxicity related to oligodendrocyte development and myelination is based on animal studies, [15,51] but recent work in brain organoids has begun to extend this work to human tissue [26,52,53,54]. Proliferation and differentiation of OLs and the formation of myelination are known to be susceptible to numerous neurotoxic insults, including alcohol, [55,56] cuprizone, [57] phenylalanine, [58,59] and lead [60,61]. Our work in this area showed that GA exposure in a rodent model resulted in chronic impairments of oligodendrocyte proliferation and differentiation in the fimbria of the hippocampus, which were linked to cognitive dysfunction [15]. The present study raises the prospect that a similar phenomenon of impaired OPC function as a cause of reduced myelin generation may be at work in human tissue. We calculate the ratio of mature OLs over OPCs to evaluate OL differentiation, but this does not represent the OL proliferation. In addition, there are other possibilities, such as selective cell loss, differential marker expression, and alteration of OL survival. The underlying molecular mechanisms have not been fully elucidated, however, the potential for doing so in human OPCs and oligodendrocytes now exists.
The other major functional finding in our study is a possible reduction in synaptic number, and it is the most significant, given the critical role that synapses play in neural computation and transmission. However, to confirm the loss of synapses, further electrical physiological evidence is necessary. This validation for the synaptic effects was not performed in the present study and should be acknowledged as a limitation, even though the model has been shown to be electrical active in our previous publications. Our analysis does not indicate how much of this may be due to putative neuronal loss versus loss of synapse density per neuron, and we would speculate that both factors are likely at work based on previous work in this field. Numerous studies in animal models have demonstrated synaptic loss because of developmental anesthetic neurotoxicity and multiple mechanisms have been suggested. To our best knowledge, this is the first demonstration of synaptic loss related to early anesthesia exposure to human tissue. While we did not pursue molecular mechanisms, numerous explanations related to cellular function are inherent in our results. It is inherently obvious that neuronal loss either via apoptosis or via failure of neural progenitor cell progression, has the potential to result in synapse loss. Early work in hippocampal and retinal primary neuron culture made it clear that astrocytes or astrocyte- conditioned media promote synaptic development, and more sophisticated work has since confirmed that astrocytes play a key role in the formation and maintenance of synaptic connections during development [62]. Interestingly, this may to some extent be a bidirectional relationship in which astrocytes are dependent on synaptic activity to engender their further growth and development [63]. The mechanisms by which astrocytes and synapses interact include numerous trophic and other signaling factors [64]. While less well-established, it is also likely that synapse formation and myelin development are engaged in a mutually supportive bi-directional developmental process through activity-dependent modulation of plasticity [65]. We speculate that the effects that are noted in our bMPS model, in which neurons, oligodendrocytes, astrocytes, myelin, and synapses all show developmental disruption, are interrelated and self-reinforcing.
Our previous work in rodent models has indicated that exposure to general anesthetics (GA) disrupts neural development in the early postnatal period by inappropriately increasing activity in the mTOR pathway, and both behavioral and histological changes could be reversed by mTOR inhibition [13,15]. In the present study, as in animal studies, we observe a lasting alteration in the tone of mTOR signaling downstream molecules in both neurons and glial cells in bMPS after SEV exposure, which is probably related to the developmental disruption. The mTOR pathway is an intracellular signaling pathway that regulates cellular activities including proliferation, differentiation, apoptosis, metabolism, transmitter release, and other biological processes [66]. In the past decade, many studies have implicated mTOR signaling in neural developmental disorders (NDDs) [67]. We propose that the effect of anesthesia neurotoxicity on neural development in bMPS model should be further studied. Anesthetic-caused genetic alterations of different cellular components, such as neurons, oligodendrocytes, astrocytes, microglial cells, myelin, and synapses, should be detected. The reagents having potential functions to mitigate impairment of the neurotoxicity, such as mTOR inhibitors, should be identified. We hope further investigations in this field with the bMPS model will be conducted.
Our findings are only the beginning of what we hope will be substantial further work in the cultured human brain tissue models to elucidate the nature of developmental anesthetic neurotoxicity and the mechanisms behind observed phenomena in human tissue. This model system does have limitations, chief among them that development is driven by genetics and non-patterned activity, rather than by external stimulation. Furthermore, the system we chose represents the human forebrain in terms of cell composition and interactions but lacks the organization of any distinct human brain area that occurs in vivo and thus does not allow for investigation of the development of specific brain circuitry that might be deemed particularly vulnerable to anesthetic toxicity. Despite the bMPS does not recapitulate all in vivo determinants including systemic physiology, vasculature, and sensory-driven patterned activity, this model remains valuable for developmental neurotoxicity studies because it recapitulates major neurodevelopmental processes, including proliferation, apoptosis, differentiation, migration, neurite outgrowth, synaptogenesis, myelination, and neural network formation. In addition, although it lacks regional organization, it supports complex neuron–glia interactions in a multicellular environment in which neurons, astrocytes, and oligodendrocytes arise through coordinated developmental waves. Terminology in the field varies, and the word “brain organoid” is often used to refer generically to any human stem cell-derived three-dimensional culture. We have chosen to refer to our system as a “brain microphysiological system” (bMPS) in contrast to approaches that result in the generation of six distinct cerebral cortical lamina or a hippocampal structure with four distinct cornu ammonis subfields, which we would deem an organoid, and this approach could be employed in future work. While stem-cell-based human brain tissue cultures cannot perfectly model the human brain, we believe they hold great promise, particularly for the investigation of genetic and epigenetic mechanisms of anesthetic neurotoxicity that may be unique to human cells. They are also ideally suited for exploring selective vulnerabilities that may arise from human neurodevelopmental genetic diseases that can be modeled by deriving stem cells directly from patients with the disease background. Furthermore, it will be of great value to discern whether previously established mechanisms of anesthetic toxicity in rodent models translate as a whole or in part in culture systems comprised of human brain tissue. Our hope is that these approaches can help address the unanswered questions that remain about how and to what degree the human brain is adversely affected by developmental anesthetic exposure.

4. Materials and Methods

4.1. Generation of bMPS

Human iPSCs NIBSC8 (passage 10-13, female origin) were kindly provided by The National Institute for Biological Standards and Control, NIBSC, UK. The iPSCs were cultured in mTeSRTM Plus medium on vitronectin-coated plates at 37 °C, 5% CO2, and 5% O2. After 80 - 90% confluency, 2 x 105 iPSC were seeded as small colonies per well in a Matrigel-coated 6–well plate in mTeSR Plus medium supplemented with 10 μM Y-27632. The next day, the medium was changed to neural induction medium (Neurobasal medium and 1X neural induction supplement), and the medium was changed every second day. During the induction period, cell morphology was checked under an optical microscope. After seven days of neural induction, iPSCs differentiated into neural progenitor cells (NPCs) and expanded into the Neural Expansion medium (50% Neurobasal, 50% Advanced DMEM/F12 media, 1X Neural induction supplement, and 5 μM Y-27632) for five passages before freezing the NPC stock. Y-27632 was removed the next day after passing. After passage five, NPCs were cultured without Y-27632. Cultures were transferred to normoxia (~20% O2, 37 °C, and 5% CO2 and the medium was changed every second day. Passages 6-15 were used to generate bMPS. NPCs at 90-100% confluence were treated with Gentle Cell Dissociation Reagent (GCDR) for 5 min at RT, after which the reagent was removed. Fresh Neural Expansion medium was added, and the NPCs were detached and counted. A total of 2 x 106 cells per well were plated in uncoated 6-well plates. After 2 days, the Neural Expansion medium was replaced with differentiation medium (B-27TM Plus Neuronal Culture System, 1% Glutamax, 0.01 μg/ml human recombinant GDNF, 0.01 μg/ml human recombinant BDNF, 1% Pen/Strep/Glutamine). Cultures were maintained at 37 °C, 5% CO2 and 20% O2 under gyratory shaking (88 rpm,19 mm orbit) for 8 weeks, with medium changes three times per week (Figure 1A). The procedures were adapted from previously established protocols [25,26,27]. In this study, total ten independent batches of bMPSs were used for all experiments including immunofluorescence (IF) staining, Western blotting (WB) and quantitative real-time PCR (qPCR). Four batches were for IF staining, and six for WB and qPCR analysis.

4.2. Sevoflurane Exposure

At the end of week 8, sevoflurane exposure in bMPS was performed in a sterile environment. 6-well plates with cultures were removed from the incubator (5% CO2, 37 °C), placed in a fume hood, and randomly divided into CON and SEV groups. The plates were placed in identical air-tight, humidified chambers (Billups-Rothenberg, Del Mar, CA, USA) as previously described [28,32]. Before the exposure day, chambers were thoroughly cleaned with 75% ethanol and UV irradiated for 24 h. After the plates were placed on the mesh of the chamber (lids were placed beside), the chambers were sealed. Both gas inlet and outlet remained open at this time. The chambers for SEV group were connected to an agent-specific tabletop portable anesthetic vaporizer (Supera-Vet, Vaporizer Sales and Services Inc., Rockmart, GA, USA) that delivered 2.4% sevoflurane mixed with 5% CO2 and 95% air carrier gas. A calibrated flowmeter was used to deliver carrier gas at a flow rate of 12 L/min. Cultures in CON group were exposed to carrier gas alone. After a 15-min equilibration, the gas inlet and outlet of the chamber were tightly clipped and the chambers with cultured bMPS were placed on the gyratory shaker inside the incubator to maintain the temperature at 37 °C (Figure 1B). According to our previous research, anesthetic concentration stability over time in this exposure technique should be ensured [35]. After 4 h of exposure, the chambers were removed from the incubator. The chambers were unsealed and the sevoflurane mix was released. The cultures were placed back in the incubator for 4 weeks and maintained under gyratory shaking, with medium changed three times per week.

4.3. bMPS Size Measurement

At the end of week 12, cultured bMPSs from both control and sevoflurane groups were observed in a phase-contrast microscope (Nikon, Japan), and photos were taken at 10× magnification. Twelve bMPSs from each group were randomly selected from two independent batches of cultures, six from each, for area (mm2) measurement using ImageJ software. Statistics (t-test) were performed with the GraphPad Prism 8 software (Figure 2A-C). Two researchers participated in the experiments. The person performing image analysis was blinded to the experimental grouping in all cases.

4.4. Immunofluorescence (IF) Staining

After size measurement at the end of week 12, bMPSs were harvested. After twice washing with 0.1M phosphate buffer saline (PBS), bMPSs were fixed with 4% paraformaldehyde (PFA) in PBS at 4 °C for 60 minutes. After 3 x 60 minutes washing with PBS in 1% BSA, bMPSs were transferred into PBS with 0.1% Triton (PBST) in a 24-well plate. The bMPSs were permeabilized in PBST at 4 °C for 30 minutes and blocked in 300 μl 100% BlockAid per well onto a shaker at 4ºC for 60 min. The bMPSs were incubated with the following primary antibodies in 10% BlockAid/Triton for 48 hours at 4ºC: (1) Rabbit anti Caspase 3 (1:100; Abcam; USA) for apoptosis. (2) Rabbit anti Nestin for neural progenitor cell (NPC) (1:250; BioVision, USA) mixed with mouse anti microtubule-associated protein 2 (MAP2) for mature neuron (1:250; Millipore Sigma, USA); (3) Rabbit anti–neural/glial antigen 2 (NG2) for oligodendrocyte progenitor cell – OPC (1:200; Millipore, USA) mixed with mouse anti–adenomatous polyposis coli-CC1 (APC-CC1) for mature oligodendrocytes (1:300; Millipore); (4) mouse anti myelin basic protein (MBP) for myelin (1:500; Santa Cruz Biotechnology, USA) mixed with rabbit anti neurofilament 200 (NF200) for axon (1:100; Sigma-Aldrich, USA); (5) Mouse anti glial fibrillary acidic protein (GFAP) for astrocyte (1:400; Millipore Sigma; USA); (6) Mouse anti synaptophysin (SYP) for synapses (1:250; Abcam, USA) mixed with rabbit anti neuronal nuclei (NeuN) for neuron (1:300; EDM Millipore, USA). and (7) rabbit anti phospho-S6 (1:1,000; Cell Signaling, USA). After 3 x 60-minute washes with 800 μl organoid washing buffer (OWB: PBS with 0.5% BSA and 0.1% Triton), bMPSs were incubated in following secondary antibodies in 10% BlockAid/Triton for 24 hours at 4ºC on a shaker: Alexa 488–conjugated goat anti-rabbit IgG (1:300; Invitrogen, USA) mixed with Cy3-conjugated goat anti-mouse IgG (1:600; Jackson ImmunoResearch Labs, USA), or Alexa 488–goat anti-mouse IgG (1:300; Invitrogen) mixed with Cy3-conjugated goat anti-rabbit IgG (1:600; Jackson ImmunoResearch Labs). After 3 x 60 mins washing with OWB and a 60-minute wash with 1% BSA in PBS, bMPSs were moved onto glass slides and cover-slipped with Immu-Mount (Epredia, UK) [15,68,69].

4.5. Quantitative Analysis of IF Labeling

Total four independent batches of cultured bMPSs were used for IF staining. The IF-stained bMPSs were observed and imaged at 10× magnification using a Leica 4000 confocal microscope (Leica, Wetzlar, Germany). In each staining, six or eight (pS6) bMPSs with immunoreactivity were randomly selected from two independent batches of cultures (2-4 from each) were selected for further quantitative analysis. Identical photo exposure was set for both groups. Based on the photos taken from these images, data analysis was quantitatively conducted with ImageJ (National Institutes of Health) [15,68,69].

4.5.1. Cell and Synaptic Puncta Counting

In IF-staining cases for Caspase3, MAP2, Nestin, NG2, APC and GFAP labeled cells, and synaptophysin positive synaptic puncta were counted. The criteria for caspase-3 positive cells were labeling of dead cell nuclei. For other neural cells, immunoreactivity was located in the cytoplasm and proximal processes, and the profiles of cells were identified. For terminal synapses, synaptophysin-positive puncta, either contacting neuronal cytoplasm/proximal dendrites or not, were counted. Labeled cells/synapses were identified with green (Alexa 488) or red (Cy3) channel. Images (10x) were opened and initialized in ImageJ. First, the region of each single bMPS was outlined using the ‘‘Freehand’’ tool. “Analyze” to “Measure” tools were selected, and the areas (mm2) were measured. Then, “Plugins,” “Analysis,” and ‘‘Cell Counter’’ tools were selected, and each labeled cell inside bMPSs was clicked, with which each counted cell was marked, preventing the same cell from being counted twice. The numbers of counted cells in each bMPS were automatically recorded after selecting “Analyze” and “Measure” tools. The number of labeled cells per mm2 was calculated according to the area of bMPSs. We then used the ratio of MAP2-positive cells over Nestin-positive cells, and the ratio of APC+ over NG2+ cells to evaluate the maturation of neuronal and oligodendrocyte lineage cells. For the study of synaptophysin+ puncta, 40x high-power images from five fields of each bMPS (up, down, left, right, center) were randomly selected. After puncta counting, average numbers of five fields were calculated.

4.5.2. Measurement of Intensity of Immunoreactivity

The labeling for MBP, GFAP, and pS6 were quantitatively evaluated by measuring intensity of immunoreactivity. The images were opened with ImageJ, and an outline of bMPSs was drawn with “Freehand” tool. The “set measurements” was selected from the analyze menu, and “integrated density” was activated. A region outside the bMPSs was selected as background. The final intensity of MBP, Synaptophysin, GFAP, and pS6 equals the measured density minus background.

4.5.3. Measurement of Myelinated Axons

The summed length of MBP+/NF200+ double-labeled axons was measured with ImageJ. In the double-labeling images, the “Freehand line” tool was selected. All MBP+/NF200+ labeled axons (yellow color) in the images were drawn and digitized by clicking “Measure” from “Analyze”. The summed length of these axons (mm) in each bMPS was calculated. For all analysis, two researchers participated in the experiments. The experimenter performing image analysis was blinded to the experimental grouping in all cases.

4.6. Western Blotting (WB)

The bMPSs were lysed and homogenized. Samples were running on Bis-Tris Protein Gels for 1 h. The proteins were transferred to nitrocellulose-blotting membranes. Blots were probed with the following antibodies: (1) sex determining region Y-box 2 (SOX2) for neural differentiation (1:1,000; Cell Signaling); (2) cortical neuronal marker T-box brain1 (TBR1) for cortical neurons (1:1,000; Cell Signaling); (3) NMDA receptor 2B (GluN2B) for postsynaptic component (1:1,000; Cell Signaling); (4) total (t)-p70S6 kinase (1:1,000; Cell Signaling) for mTOR signaling pathway; (5) phosphor (p)-p70S6K for phosphorylation of mTOR molecule (1:1,000; Cell Signaling); (6) β-actin as a standard marker (1:2,000; Cell Signaling). Blots were visualized using the ECL substrate kit. Images were acquired using the ChemiDoc imaging system and quantified with ImageJ (NIH) software. The ratios of band density of SOX2, TBR1, GluN2B, and t-p70S6K over β-actin; and p-p70S6K over t-p70S6K were calculated. Quantitative analysis was performed with the ImageJ program [15,67,68]. Total six independent batches of bMPSs were used for WB study. For each batch, all bMPSs from CON or SEV group were lysed with lysis buffer and the supernatant from each group of each batch was collected into 10 small tubes (15 ul, 20-30 mg/ml protein) for gel running. For each tested molecule, six blots with immunoreactivity from two batches of six batches were selected for further quantitative analysis.

4.7. Quantitative Real-Time PCR (qPCR)

The bMPSs were harvested and immediately frozen with dry ice and then placed into −80° C for storage. The total RNA samples were prepared with “RNeasy Plus Micro Kit” (Qiagen; Cat# 74034). The cDNA synthesis reaction was performed with 1 µg of total RNA plus “qScript cDNA SuperMix” (Quanta Biosciences; Cat# 95048). The specific human vGlut1 primers (forward: GAGCGCAAGTACATCGAGGA; reverse: CTTGCTGATCTCGAAGCCGA) and h-GAPDH primers (forward: CAGAAGACTGTGGATGGCCC; reverse: CCACCTGGTGCTCAGTGTAG) were applied. The qPCR was performed using Bio-Rad “Sso-Advanced Universal SYBR Green Supermix” (Cat# 172-5271) in Bio-Rad CFX96 Real-Time System (C1000 Tauch Thermal Cycler). The annealing temperature was set at 53 °C and the number of cycles was fixed to 40. Gene expression was analyzed using the ∆∆CT method with CT as the threshold cycle. The relative levels of target genes, which were normalized to naive CT value, were reported as 2−∆∆CT [69,70]. Total six independent batches of bMPSs were used for qPCR study. Because large quantity of bMPS is needed for this study, the tissue for six repeated tests were collected from six batches.

4.8. Statistical Analysis

Statistical analyses were performed using GraphPad Prism 8 software. Sample size was determined based on prior experience with this experimental design. Samples for each individual experiment (IF, WB, or qPCR for each marker) were taken from two independent bMPS cultures (two independent batches), and sample sizes are reported on a per-group basis. All comparisons between the CON and SEV groups were analyzed using the two-tailed Mann–Whitney U test (nonparametric test). Statistical significance was defined as p < 0.05. Data are presented as mean ± standard deviation (SD). The sample size (n) represents the number of bMPS analyzed per group.

5. Conclusions

In summary, we conclude that sevoflurane, a frequently used general anesthetic reagent, substantially inhibits the brain development in a human brain microphysiological system (bMPS) and has lasting disruptive effects on the development of neurons, oligodendrocytes, and astrocytes, resulting in loss of synapses and reduced myelination. The aberrant over-expression of mTOR signal molecules is probably involved in this alteration.

Author Contributions

Conceptualization, C.D.M. and Q.L.; methodology, C.D.M. and Q.L.; formal analysis, Q.L.; investigation, Q.L., L.D., E.M.P., N.M.; resources, C.D.M. and Q.L.; data curation, Q.L.; writing—original draft preparation, Q.L.; writing—review and editing, C.D.M.; supervision, C.D.M. and L.S.; project administration, C.D.M. and Q.L.; funding acquisition, C.D.M. and Q.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Institute of General Medicine to C.D.M. (1R01GM120519-01), ACCM StAAR Award to Q.L. (80067671), and Maryland Stem Cell Research Fund (MSCRF) to Q.L (2025-R2-MSCRFL-00025).

Institutional Review Board Statement

The human bMPS protocol was approved by Institutional Review Board (IRB) and Stem Cell Research Committee (SCRO) of Johns Hopkins University (Protocol number: IRB00507275; Approval date: June 10, 2025) for studies involving human stem cells.

Data Availability Statement

The datasets generated during the current study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

iPSC induced pluripotent stem cell
bMPS brain microphysiological system
NPC neural progenitor cell
SEV sevoflurane
CON control
IF immunofluorescence
WB Western blotting
qPCR quantitative real time polymerase chain reaction
MAP2 microtubule-associated protein 2
SOX2 sex determining region Y-box 2
TBR1 T-box brain1
GA General Anesthetic
mTOR Mammalian target of rapamycin
pS6 Phospho-S6 ribosomal protein
p-p70S6K phosphorylated ribosomal protein 70 S6 kinase
t-p70S6K total ribosomal protein 70 S6 kinase
GFAP Glial fibrillary acidic protein
APC adenomatous polyposis coli
OL oligodendrocyte
NG2 neural/glial antigen 2
OPC oligodendrocyte progenitor cell
MBP Myelin basic protein
NF neurofilament
vGLUt1 vesicular glutamate transportor-1
GluN2B glutamate NMDA receptor-2B
mRNA Message ribonucleic acid
NDD neurodevelopmental disorder

References

  1. Weiser, T.G.; Regenbogen, S.E.; Thompson, K.D.; Haynes, A.B.; Lipsitz, S.R.; Berry, W.R.; Gawande, A.A. An estimation of the global volume of surgery: A modelling strategy based on available data. Lancet 2008, 372, 139–144. [Google Scholar] [CrossRef] [PubMed]
  2. Eckenhoff, J.E. Relationship of anesthesia to postoperative personality changes in children. AMA Am J Dis Child 1953, 86, 587–591. [Google Scholar] [CrossRef]
  3. Davidson, A.J.; Disma, N.; de Graaff, J.C.; Withington, D.E.; Dorris, L.; Bell, G.; Stargatt, R.; Bellinger, D.C.; Schuster, T.; Arnup, S.J.; et al. Neurodevelopmental outcome at 2 years of age after general anaesthesia and awake-regional anaesthesia in infancy (GAS): An international multicentre, randomised controlled trial. Lancet 2016, 387, 239–250. [Google Scholar] [CrossRef]
  4. Sun, L.S.; Li, G.; Miller, T.L.; Salorio, C.; Byrne, M.W.; Bellinger, D.C.; Ing, C.; Park, R.; Radcliffe, J.; Hays, S.R.; et al. Association between a single general anesthesia exposure before age 36 months and neurocognitive outcomes in later childhood. JAMA 2016, 315, 2312–2320. [Google Scholar] [CrossRef]
  5. Vutskits, L.; Xie, Z. Lasting impact of general anaesthesia on the brain: Mechanisms and relevance. Nat Rev Neurosci 2016, 17, 705–717. [Google Scholar] [CrossRef]
  6. Kachmar, A.G.; Irving, S.Y.; Connolly, C.A.; Curley, M.A. A systematic review of risk factors associated with cognitive impairment after pediatric critical illness. Pediatr Crit Care Med 2018, 19, e164–e171. [Google Scholar] [CrossRef]
  7. Ing, C.; Brambrink, A.M. Mayo Anesthesia Safety in Kids continued: two new studies and a potential redirection of the field. Br J Anaesth 2019, 122, 716–719. [Google Scholar] [CrossRef] [PubMed]
  8. Warner, D.O.; Zaccariello, M.J.; Katusic, S.K.; Schroeder, D.R.; Hanson, A.C.; Schulte, P.J.; Buenvenida, S.L.; Gleich, S.J.; Wilder, R.T.; Sprung, J.; et al. Neuropsychological and Behavioral Outcomes after Exposure of Young Children to Procedures Requiring General Anesthesia: The Mayo Anesthesia Safety in Kids (MASK) Study. Anesthesiology 2018, 129, 89–105. [Google Scholar] [CrossRef] [PubMed]
  9. Ing, C.; Vutskits, L. Unanswered questions of anesthesia neurotoxicity in the developing brain. Curr Opin Anaesthesiol 2023, 36, 510–515. [Google Scholar] [CrossRef]
  10. Center for Drug Evaluation and Research. Drug safety and availability – FDA drug safety communication: FDA review results in new warnings about using general anesthetics and sedation drugs in young children and pregnant women. U.S. Food and Drug Administration. Available online: https://www.fda.gov/Drugs/DrugSafety/ (accessed on 19 July 2019).
  11. Clausen, N.G.; Hansen, T.G.; Disma, N. Anesthesia Neurotoxicity in the Developing Brain: Basic Studies Relevant for Neonatal or Perinatal Medicine. Clin Perinatol 2019, 46, 647–656. [Google Scholar] [CrossRef] [PubMed]
  12. Useinovic, N.; Jevtovic-Todorovic, V. Controversies in anesthesia-induced developmental neurotoxicity. Best Pract Res Clin Anaesthesiol 2023, 37, 28–29. [Google Scholar] [CrossRef] [PubMed]
  13. Kang, E.; Jiang, D.; Ryu, Y.; Lim, S.; Kwak, M.; Gray, C.D.; Xu, M.; Choi, J.H.; Junn, S.; Jieun Kim, J.; et al. Early postnatal exposure to isoflurane causes cognitive deficits and disrupts development of newborn hippocampal neurons via activation of the mTOR pathway. PLoS Biol 2017, 15, e2001246. [Google Scholar] [CrossRef]
  14. Neudecker, V.; Xu, J.; Thomas, M.A.; Penberthy, K.K.; Kang, E.; Berg, D.A.; O’Meara, A.M.I.; Brambrink, A.M.; Mintz, C.D. An Update on Preclinical Research in Anesthetic-Induced Developmental Neurotoxicity in Nonhuman Primate and Rodent Models. J Neurosurg Anesthesiol 2023, 35, 104–113. [Google Scholar] [CrossRef] [PubMed]
  15. Li, Q.; Mathena, R.P.; Xu, J.; Eregha, O.N.; Wen, J.; Mintz, C.D. Early postnatal exposure to isoflurane disrupts oligodendrocyte development and myelin formation in the mouse hippocampus. Anesthesiology 2019, 131, 1077–1091. [Google Scholar] [CrossRef]
  16. Wagner, M.; Ryu, Y.K.; Smith, S.C.; Patel, P.; Mintz, C.D. Review: effects of anesthetics on brain circuit formation. J Neurosurg Anesthesiol 2014, 26, 358–362. [Google Scholar] [CrossRef]
  17. Borzage, M.T.; Peterson, B.S. A Scoping Review of the Mechanisms Underlying Developmental Anesthetic Neurotoxicity. Anesth Analg 2025, 140, 409–426. [Google Scholar] [CrossRef]
  18. Bleeser, T.; Brenders, A.; Hubble, T.R.; Van de Velde, M.; Deprest, J.; Rex, S.; Devroe, S. Preclinical evidence for anaesthesia-induced neurotoxicity. Best Pract Res Clin Anaesthesiol 2023, 37, 16–27. [Google Scholar] [CrossRef]
  19. Jevtovic-Todorovic, V. General Anesthetics and Neurotoxicity: How Much Do We Know? Anesthesiol Clin 2016, 34, 439–451. [Google Scholar] [CrossRef]
  20. Lancaster, M.A.; Renner, M.; Martin, C.A.; Wenzel, D.; Bicknell, L.S.; Hurles, M.E.; Homfray, T.; Penninger, J.M.; Jackson, A.P.; Knoblich, J.A. Cerebral organoids model human brain development and microcephaly. Nature 2013, 501, 373–379. [Google Scholar] [CrossRef] [PubMed]
  21. Hogberg, H.T.; Smirnova, L. The future of 3D brain cultures in developmental neurotoxicity testing. Front Toxicol 2022, 4, 808620. [Google Scholar] [CrossRef]
  22. Qian, X.; Nguyen, H.N.; Song, M.M.; Hadiono, C.; Ogden, S.C.; Hammack, C.; Yao, B.; Hamersky, G.R.; Jacob, F.; Zhong, C.; et al. Brain-region-specific organoids using mini-bioreactors for modeling ZIKV exposure. Cell 2016, 165, 1238–1254. [Google Scholar] [CrossRef]
  23. Mulder, L.A.; Depla, J.A.; Sridhar, A.; Wolthers, K.; Pajkrt, D.; Vieira de Sá, R. A beginner’s guide on the use of brain organoids for neuroscientists: a systematic review. Stem Cell Res Ther 2023, 14, 87. [Google Scholar] [CrossRef]
  24. Setia, H.; Muotri, A.R. Brain organoids as a model system for human neurodevelopment and disease. Semin Cell Dev Biol 2019, 95, 93–97. [Google Scholar] [CrossRef]
  25. Pamies, D.; Barreras, P.; Block, K.; Makri, G.; Kumar, A.; Wiersma, D.; Smirnova, L.; Zang, C.; Bressler, J.; Christian, K.M.; et al. A human brain microphysiological system derived from induced pluripotent stem cells to study neurological diseases and toxicity. ALTEX 2017, 34, 362–376. [Google Scholar] [CrossRef]
  26. Morales Pantoja, I.E.; Ding, L.; Leite, P.E.C.; Marques, S.A.; Romero, J.C.; Alam El Din, D.M.; Zack, D.; Chamling, X.; Smirnova, L. A novel approach to increase glial cell populations in brain microphysiological systems. Adv Biol 2023, 8, e2300198. [Google Scholar] [CrossRef]
  27. Romero, JC.; Berlinicke, C.; Chow, S.; Duan, Y.; Wang, Y.; Chamling, X.; Smirnova, L. Oligodendrogenesis and myelination tracing in a CRISPR/Cas9-engineered brain microphysiological system. Front Cell Neurosci 2023, 16, 1094291. [Google Scholar] [CrossRef]
  28. Xu, J.; Mathena, R.P.; Xu, M.; Wang, Y.; Chang, C.; Fang, Y.; Zhang, P.; Mintz, C.D. Early Developmental Exposure to General Anesthetic Agents in Primary Neuron Culture Disrupts Synapse Formation via Actions on the mTOR Pathway. Int J Mol Sci 2018, 19, 2183. [Google Scholar] [CrossRef] [PubMed]
  29. Yang, M.; Wei, H. Anesthetic neurotoxicity: Apoptosis and autophagic cell death mediated by calcium dysregulation. Neurotoxicol Teratol 2017, 60, 59–62. [Google Scholar] [CrossRef]
  30. Kang, E.; Berg, D.A.; Furmanski, O.; Jackson, W.M.; Ryu, Y.K.; Gray, C.D.; Mintz, C.D. Neurogenesis and developmental anesthetic neurotoxicity. Neurotoxicol Teratol 2017, 60, 33–39. [Google Scholar] [CrossRef] [PubMed]
  31. Ryu, Y.K.; Khan, S.; Smith, S.C.; Mintz, C.D. Isoflurane impairs the capacity of astrocytes to support neuronal development in a mouse dissociated coculture model. J Neurosurg Anesthesiol 2014, 26, 363–368. [Google Scholar] [CrossRef] [PubMed]
  32. Mintz, C.D.; Barrett, K.M.; Smith, S.C.; Benson, D.L.; Harrison, N.L. Anesthetics interfere with axon guidance in developing mouse neocortical neurons in vitro via a gamma-aminobutyric acid type A receptor mechanism. Anesthesiology 2013, 118, 825–833. [Google Scholar] [CrossRef]
  33. Jiang, M.; Tang, T.; Liang, X.; Li, J.; Qiu, Y.; Liu, S.; Bian, S.; Xie, Y.; Fang, F.; Cang, J. Maternal sevoflurane exposure induces temporary defects in interkinetic nuclear migration of radial glial progenitors in the fetal cerebral cortex through the Notch signalling pathway. Cell Prolif 2021, 54, e13042. [Google Scholar] [CrossRef]
  34. Shang, J.; Li, B.; Fan, H.; Liu, P.; Zhao, W.; Chen, T.; Chen, P.; Yang, L. Sevoflurane promotes premature differentiation of dopaminergic neurons in hiPSC-derived midbrain organoids. Front Cell Dev Biol 2022, 10, 941984. [Google Scholar] [CrossRef]
  35. Mintz, C.D.; Smith, S.C.; Barrett, K.M.; Benson, D.L. Anesthetics interfere with the polarization of developing cortical neurons. J Neurosurg Anesthesiol. 2012, 24, 368–375. [Google Scholar] [CrossRef]
  36. Lerman, J.; Sikich, N.; Kleinman, S.; Yentis, S. The pharmacology of sevoflurane in infants and children. Anesthesiology 1994, 80, 814–824. [Google Scholar] [CrossRef]
  37. Xu, J.; Mathena, R.P.; Singh, S.; Kim, J.; Long, J.J.; Li, Q.; Junn, S.; Blaize, E.; Mintz, C.D. Early developmental exposure to repetitive long duration of midazolam sedation causes behavioral and synaptic alterations in a rodent model of neurodevelopment. J Neurosurg Anesthesiol 2019, 31, 151–162. [Google Scholar] [CrossRef] [PubMed]
  38. Jevtovic-Todorovic, V.; Hartman, R.E.; Izumi, Y.; Benshoff, N.D.; Dikranian, K.; Zorumski, C.F.; Olney, J.W.; Wozniak, D.F. Early exposure to common anesthetic agents causes widespread neurodegeneration in the developing rat brain and persistent learning deficits. J Neurosci 2003, 23, 876–882. [Google Scholar] [CrossRef]
  39. Wang, L.; Owusu-Hammond, C.; Sievert, D.; Gleeson, J.G. Stem cell-based organoid models of neurodevelopmental disorders. Biol Psychiatry 2023, 93, 622–631. [Google Scholar] [CrossRef] [PubMed]
  40. Trujillo, C.A.; Muotri, A.R. Brain organoids and the study of neurodevelopment. Trends Mol Med 2018, 24, 982–990. [Google Scholar] [CrossRef] [PubMed]
  41. Trimborn, M.; Bell, S.M.; Felix, C.; Rashid, Y.; Jafri, H.; Griffiths, P.D.; Neumann, L.M.; Krebs, A.; Reis, A.; Sperling, K.; et al. Mutations in microcephalin cause aberrant regulation of chromosome condensation. Am J Hum Genet 2004, 75, 261–266. [Google Scholar] [CrossRef]
  42. Thomas, C.A.; Tejwani, L.; Trujillo, C.A.; Negraes, P.D.; Herai, R.H.; Mesci, P.; Macia, A.; Crow, Y.J.; Muotri, A.R. Modeling of TREX1-dependent autoimmune disease using human stem cells highlights L1 accumulation as a source of neuroinflammation. Cell Stem Cell 2017, 21, 319–331. [Google Scholar] [CrossRef]
  43. Krenn, V.; Bosone, C.; Burkard, T.R.; Spanier, J.; Kalinke, U.; Calistri, A.; Salata, C.; Rilo Christoff, R.; Pestana Garcez, P.; Mirazimi, A.; et al. Organoid modeling of Zika and herpes simplex virus 1 infections reveals virus-specific responses leading to microcephaly. Cell Stem Cell 2021, 28, 1362–1379. [Google Scholar] [CrossRef]
  44. Majolo, F.; Marinowic, D.R.; Moura, A.Á.; Machado, D.C.; da Costa, J.C. Use of induced pluripotent stem cells (iPSCs) and cerebral organoids in modeling the congenital infection and neuropathogenesis induced by Zika virus. J Med Virol 2019, 91, 525–532. [Google Scholar] [CrossRef]
  45. Nickel, M.; Gu, C. Regulation of central nervous system myelination in higher brain functions. Neural Plast 2018, 2018, 6436453. [Google Scholar] [CrossRef]
  46. Williamson, J.M.; Lyons, D.A. Myelin Dynamics Throughout Life: An Ever-Changing Landscape. Front Cell Neurosci 2018, 12, 424. [Google Scholar] [CrossRef]
  47. Chiaravalloti, N.D.; DeLuca, J. Cognitive impairment in multiple sclerosis. Lancet Neurol 2008, 7, 1139–1151. [Google Scholar] [CrossRef] [PubMed]
  48. Bourbon-Teles, J.; Bells, S.; Jones, D.K.; Coulthard, E.; Rosser, A.; Metzler-Baddeley, C. Myelin breakdown in human Huntington’s disease: multi-modal evidence from diffusion MRI and quantitative magnetization transfer. Neurosci 2019, 403, 79–92. [Google Scholar] [CrossRef] [PubMed]
  49. Li, Q.; Houdayer, T.; Liu, S.; Belegu, V. Induced neural activity promotes an oligodendroglia regenerative response in the injured spinal cord and improves motor function after spinal cord injury. J Neurotrauma 2017, 34, 3351–3361. [Google Scholar] [CrossRef]
  50. Brambrink, A.M.; Back, S.A.; Riddle, A.; Gong, X.; Moravec, M.D.; Dissen, G.A.; Creeley, C.E.; Dikranian, K.T.; Olney, J.W. Isoflurane-induced apoptosis of oligodendrocytes in the neonatal primate brain. Ann Neurol 2012, 72, 525–535. [Google Scholar] [CrossRef]
  51. Jiang, D.; Lim, S.; Kwak, M.; Ryu, Y.K.; Mintz, C.D. The changes of oligodendrocytes induced by anesthesia during brain development. Neural Regen Res 2015, 10, 1386–1387. [Google Scholar] [PubMed]
  52. Marton, R.M.; Miura, Y.; Sloan, S.A.; Li, Q.; Revah, O.; Levy, R.J.; Huguenard, J.R.; Pașca, S.P. Differentiation and maturation of oligodendrocytes in human three-dimensional neural cultures. Nat Neurosci 2019, 22, 484–491. [Google Scholar] [CrossRef]
  53. Chesnut, M.; Paschoud, H.; Repond, C.; Smirnova, L.; Hartung, T.; Zurich, M.G.; Hogberg, H.T.; Pamies, D. Human iPSC-derived model to study myelin disruption. Int J Mol Sci 2021, 22, 9473. [Google Scholar] [CrossRef]
  54. Chesnut, M.; Hartung, T.; Hogberg, H.; Pamies, D. Human oligodendrocytes and myelin In vitro to evaluate developmental neurotoxicity. Int J Mol Sci 2021, 22, 7929. [Google Scholar] [CrossRef]
  55. Okamoto, H.; Miki, T.; Lee, K.Y.; Yokoyama, T.; Kuma, H.; Wang, Z.Y.; Gu, H.; Li, H.P.; Matsumoto, Y.; Irawan, S.; et al. Oligodendrocyte myelin glycoprotein (OMgp) in rat hippocampus is depleted by chronic ethanol consumption. Neurosci Lett 2006, 406, 76–80. [Google Scholar] [CrossRef] [PubMed]
  56. Arzua, T.; Yan, Y.; Jiang, C.; Logan, S.; Allison, R.L.; Wells, C.; Kumar, S.N.; Schäfer, R.; Bai, X. Modeling alcohol-induced neurotoxicity using human induced pluripotent stem cell-derived three-dimensional cerebral organoids. Transl Psychiatry 2020, 10, 347. [Google Scholar] [CrossRef]
  57. Zirngibl, M.; Assinck, P.; Sizov, A.; Caprariello, A.V.; Plemel, J.R. Oligodendrocyte death and myelin loss in the cuprizone model: an updated overview of the intrinsic and extrinsic causes of cuprizone demyelination. Mol Neurodegener 2022, 17, 34. [Google Scholar] [CrossRef]
  58. Schoemans, R.; Aigrot, M.S.; Wu, C.; Marée, R.; Hong, P.; Belachew, S.; Josse, C.; Lubetzki, C.; Bours, V. Oligodendrocyte development and myelinogenesis are not impaired by high concentrations of phenylalanine or its metabolites. Inherit Metab Dis 2010, 33, 113–120. [Google Scholar] [CrossRef] [PubMed]
  59. Kim, J.; Lee, S.; Lee, J.; Park, J.C.; Kim, K.H.; Ko, J.M.; Park, S.H.; Kim, S.K.; Mook-Jung, I.; Lee, J.Y. Neurotoxicity of phenylalanine on human iPSC-derived cerebral organoids. Mol Genet Metab 2022, 136, 132–144. [Google Scholar] [CrossRef] [PubMed]
  60. Rothenberg, S.J.; Poblano, A.; Garza-Morales, S. Prenatal and perinatal low level lead exposure alters brainstem auditory evoked responses in infants. NeuroToxicology 1994, 15, 695–699. [Google Scholar]
  61. Tiffany-Castiglioni, E. Cell culture models for lead toxicity in neuronal and glial cells. NeuroToxicology 1993, 14, 513–536. [Google Scholar]
  62. Fossati, G.; Matteoli, M.; Menna, E. Astrocytic Factors Controlling Synaptogenesis: A Team Play. Cells 2020, 9, 2173. [Google Scholar] [CrossRef]
  63. Stipursky, J.; Romão, L.; Tortelli, V.; Neto, V.M.; Gomes, F.C. Neuron-glia signaling: Implications for astrocyte differentiation and synapse formation. Life Sci 2011, 89, 524–531. [Google Scholar] [CrossRef] [PubMed]
  64. Chung, W.S.; Allen, N.J.; Eroglu, C. Astrocytes Control Synapse Formation, Function, and Elimination. Cold Spring Harb Perspect Biol 2015, 7, a020370. [Google Scholar] [CrossRef] [PubMed]
  65. Fields, R.D. Myelination: an overlooked mechanism of synaptic plasticity. Neuroscientist 2005, 11, 528–531. [Google Scholar] [CrossRef]
  66. Laplante, M.; Sabatini, D.M. mTOR signaling in growth control and disease. Cell 2012, 149, 274–293. [Google Scholar] [CrossRef] [PubMed]
  67. Costa-Mattioli, M.; Monteggia, L.M. mTOR complexes in neurodevelopmental and neuropsychiatric disorders. Nat Neurosci 2013, 16, 1537–1543. [Google Scholar] [CrossRef]
  68. Li, Q.; Mathena, R.P.; Eregha, O.N.; Mintz, C.D. Effects of Early Exposure of Isoflurane on Chronic Pain via the Mammalian Target of Rapamycin Signal Pathway. Int J Mol Sci. 2019, 20, 5102. [Google Scholar] [CrossRef]
  69. Li, Q.; Mathena, R.P.; Li, F.; Dong, X.; Guan, Y.; Mintz, C.D. Effects of Early Exposure to Isoflurane on Susceptibility to Chronic Pain Are Mediated by Increased Neural Activity Due to Actions of the Mammalian Target of the Rapamycin Pathway. Int J Mol Sci. 2023, 24, 13760. [Google Scholar] [CrossRef]
  70. Lee, J.A.; Bae, D.H.; Choi, W.H.; Cho, C.H.; Bang, Y.S.; Yoo, J. Effects of Sevoflurane Exposure on Fetal Brain Development Using Cerebral Organoids. J Mol Neurosci. 2022, 72, 2440–2450. [Google Scholar] [CrossRef]
Figure 1. Experimental Model System (A). Timeline for human induced pluripotent stem cell (iPSC) derived brain microphysiological system (bMPS) generation. The iPSCs were induced into neural progenitor cells (NPCs) in a 7-day neural induction. NPCs were cultured in uncoated 6-well plates under gyratory shaking (88 rpm) for eight weeks to produce bMPSs. (B). Sevoflurane exposure at end of week 8. The bMPSs were placed in an air-tight plastic chamber, which was connected to anesthetic vaporizer that delivers 2.4% sevoflurane mixed with carrier gas (5% CO2 and 95% air) for 15 min equilibration. Cultures exposed to carrier gas alone were used as control. Then, the chamber was placed on gyratory shaker in culture incubator (5% CO2) for 4 hours. After sevoflurane exposure, brain bMPSs were removed from the chamber and continuously cultured until the end of week 12.
Figure 1. Experimental Model System (A). Timeline for human induced pluripotent stem cell (iPSC) derived brain microphysiological system (bMPS) generation. The iPSCs were induced into neural progenitor cells (NPCs) in a 7-day neural induction. NPCs were cultured in uncoated 6-well plates under gyratory shaking (88 rpm) for eight weeks to produce bMPSs. (B). Sevoflurane exposure at end of week 8. The bMPSs were placed in an air-tight plastic chamber, which was connected to anesthetic vaporizer that delivers 2.4% sevoflurane mixed with carrier gas (5% CO2 and 95% air) for 15 min equilibration. Cultures exposed to carrier gas alone were used as control. Then, the chamber was placed on gyratory shaker in culture incubator (5% CO2) for 4 hours. After sevoflurane exposure, brain bMPSs were removed from the chamber and continuously cultured until the end of week 12.
Preprints 205133 g001
Figure 2. Effect of sevoflurane exposure on bMPS size and apoptosis. (A-B). At end of week 12 (4 weeks after exposure), bMPSs were observed in light microscope and photos were taken with 10 x magnification. The surface area of bMPSs from control (CON) and sevoflurane exposure (SEV) groups were measured with ImageJ. Scale bar = 0.5 mm. (C). Quantitative comparison reveals a smaller bMPS size with sevoflurane exposure. n = 12 per group; *: p < 0.05. Error bars: SD. (D-E). The bMPSs were immunostained with apoptotic marker cleaved Caspase-3 (green). Scale bar = 0.2 mm. (F). Statistics. SEV significantly increases the number of caspase-3 positive cells. n = 6; *: p < 0.05; **: p < 0.01.
Figure 2. Effect of sevoflurane exposure on bMPS size and apoptosis. (A-B). At end of week 12 (4 weeks after exposure), bMPSs were observed in light microscope and photos were taken with 10 x magnification. The surface area of bMPSs from control (CON) and sevoflurane exposure (SEV) groups were measured with ImageJ. Scale bar = 0.5 mm. (C). Quantitative comparison reveals a smaller bMPS size with sevoflurane exposure. n = 12 per group; *: p < 0.05. Error bars: SD. (D-E). The bMPSs were immunostained with apoptotic marker cleaved Caspase-3 (green). Scale bar = 0.2 mm. (F). Statistics. SEV significantly increases the number of caspase-3 positive cells. n = 6; *: p < 0.05; **: p < 0.01.
Preprints 205133 g002
Figure 3. Effects of SEV exposure on neural and neuronal development. (A, B) Representative images of bMPSs immunostained with early NPC marker Nestin (green) and mature neuron marker MAP2 (red). Boxes in photos A and B represent the locations of high-power images at lower right corners respectively. Scale bar = 0.2 mm. The number of Nestin labeled cells is increased in the SEV group as compared to the CON group. By contrast MAP2 positive cell number in SEV is lower than in CON. Consequently, the ratio of MAP2+ cells over Nestin+ NPCs is reduced by SEV exposure. (C). Quantitative data for A and B. (D). Western blotting images for neural differentiation marker SOX2 and cortical neuronal marker TBR1 in bMPS. SEV decreases the ratio of SOX2 intensity over β-actin, which indicates the neural differentiation is reduced. Reduction of the level of TBR1 means neuronal development is partially impeded. Quantitative results are shown in (E) and (F). n = 6. ns: no significance. *: p < 0.05; **: p < 0.01.
Figure 3. Effects of SEV exposure on neural and neuronal development. (A, B) Representative images of bMPSs immunostained with early NPC marker Nestin (green) and mature neuron marker MAP2 (red). Boxes in photos A and B represent the locations of high-power images at lower right corners respectively. Scale bar = 0.2 mm. The number of Nestin labeled cells is increased in the SEV group as compared to the CON group. By contrast MAP2 positive cell number in SEV is lower than in CON. Consequently, the ratio of MAP2+ cells over Nestin+ NPCs is reduced by SEV exposure. (C). Quantitative data for A and B. (D). Western blotting images for neural differentiation marker SOX2 and cortical neuronal marker TBR1 in bMPS. SEV decreases the ratio of SOX2 intensity over β-actin, which indicates the neural differentiation is reduced. Reduction of the level of TBR1 means neuronal development is partially impeded. Quantitative results are shown in (E) and (F). n = 6. ns: no significance. *: p < 0.05; **: p < 0.01.
Preprints 205133 g003
Figure 4. Effect of sevoflurane exposure on oligodendrocyte differentiation and myelination. (A, B). Representative images of bMPSs labeled with OPC marker NG2 (green) and the mature oligodendrocyte marker APC (red). Boxes in (A) and (B) are shown at higher magnification respectively. (C). The ratio of mature oligodendrocytes (APC+) over OPCs (NG2+) is significantly decreased in SEV than CON groups, which means SEV exposure inhibits oligodendrocyte differentiation. (D-I). Images of bMPSs from CON (D, E, F) and SEV exposure (G, H, I) groups. bMPSs were double stained with the myelin specific marker MBP (D, G) and axonal marker NF200 (E, H). The intensity of MBP immunoreactivity (D, G) and the summed length of MBP+/NF200+ double-labeled processes, which represent myelinated axons, were measured in the merged images (F, I). There is significant reduction in MBP+ immunoreactivity (J) and myelinated axon length (K) in the SEV group as compared to the CON group. Scale bars = 0.2 mm. n = 6; *: p < 0.05.
Figure 4. Effect of sevoflurane exposure on oligodendrocyte differentiation and myelination. (A, B). Representative images of bMPSs labeled with OPC marker NG2 (green) and the mature oligodendrocyte marker APC (red). Boxes in (A) and (B) are shown at higher magnification respectively. (C). The ratio of mature oligodendrocytes (APC+) over OPCs (NG2+) is significantly decreased in SEV than CON groups, which means SEV exposure inhibits oligodendrocyte differentiation. (D-I). Images of bMPSs from CON (D, E, F) and SEV exposure (G, H, I) groups. bMPSs were double stained with the myelin specific marker MBP (D, G) and axonal marker NF200 (E, H). The intensity of MBP immunoreactivity (D, G) and the summed length of MBP+/NF200+ double-labeled processes, which represent myelinated axons, were measured in the merged images (F, I). There is significant reduction in MBP+ immunoreactivity (J) and myelinated axon length (K) in the SEV group as compared to the CON group. Scale bars = 0.2 mm. n = 6; *: p < 0.05.
Preprints 205133 g004
Figure 5. Sevoflurane reduces astrocyte production. (A-F). Images of bMPSs from CON (A-C) and SEV exposure (D-F) groups. bMPSs were double stained with GFAP for astrocytes (red) and TOTO for nuclei (blue). Compared to CON, both GFAP reactivity (G) and the number of GFAP+ cells (H) are reduced in SEV exposure group. Scale bar = 0.2 mm. n = 6; ns: no significance **: p < 0.01.
Figure 5. Sevoflurane reduces astrocyte production. (A-F). Images of bMPSs from CON (A-C) and SEV exposure (D-F) groups. bMPSs were double stained with GFAP for astrocytes (red) and TOTO for nuclei (blue). Compared to CON, both GFAP reactivity (G) and the number of GFAP+ cells (H) are reduced in SEV exposure group. Scale bar = 0.2 mm. n = 6; ns: no significance **: p < 0.01.
Preprints 205133 g005
Figure 6. Sevoflurane exposure impairs synaptogenesis in human bMPS. (A, B) Representative images of bMPSs double labeled with the presynaptic marker, synaptophysin (SYP; red), and the neuronal marker NeuN (green), with inset boxes representing the location of high magnification images respectively. The number of synaptophysin positive puncta is reduced in SEV group as compared to CON. (C). Statistics for synaptophysin positive puncta. (D). The mRNA of excitatory presynaptic molecule vGlut1 was examined with qPCR. SEV exposure dramatically attenuates the level of vGlut1 than control. (E). Western blotting data show SEV exposure reduces the expression of excitatory postsynaptic molecule GluN2B in human bMPS. (F). Statistics for E. Scale bar = 0.2 mm. n = 6; ns: no significance; *: p < 0.05; ***: p < 0.001.
Figure 6. Sevoflurane exposure impairs synaptogenesis in human bMPS. (A, B) Representative images of bMPSs double labeled with the presynaptic marker, synaptophysin (SYP; red), and the neuronal marker NeuN (green), with inset boxes representing the location of high magnification images respectively. The number of synaptophysin positive puncta is reduced in SEV group as compared to CON. (C). Statistics for synaptophysin positive puncta. (D). The mRNA of excitatory presynaptic molecule vGlut1 was examined with qPCR. SEV exposure dramatically attenuates the level of vGlut1 than control. (E). Western blotting data show SEV exposure reduces the expression of excitatory postsynaptic molecule GluN2B in human bMPS. (F). Statistics for E. Scale bar = 0.2 mm. n = 6; ns: no significance; *: p < 0.05; ***: p < 0.001.
Preprints 205133 g006
Figure 7. Effect of sevoflurane exposure on mTOR signaling pathway. (A, B). The bMPSs were immunostained for representative mTOR pathway marker pS6. Immunoreactivity is observed in all cellular components including neurons, glial cells, and axons. (C). The intensity of reactivity of pS6 in SEV cases is evidently higher than control, which means sevoflurane exposure aberrantly increases the expression of mTOR pathway. Scale bar = 0.2 mm; n = 8. (D). Western blotting analysis shows expression and phosphorylation of a specific downstream effector molecule p70S6 kinase. (E, F). SEV exposure partially enhances the level of total (t-) p70S6K and significantly increases phosphorylated (p-) p70S6K. n = 6; ns: no significance; **: p < 0.01.
Figure 7. Effect of sevoflurane exposure on mTOR signaling pathway. (A, B). The bMPSs were immunostained for representative mTOR pathway marker pS6. Immunoreactivity is observed in all cellular components including neurons, glial cells, and axons. (C). The intensity of reactivity of pS6 in SEV cases is evidently higher than control, which means sevoflurane exposure aberrantly increases the expression of mTOR pathway. Scale bar = 0.2 mm; n = 8. (D). Western blotting analysis shows expression and phosphorylation of a specific downstream effector molecule p70S6 kinase. (E, F). SEV exposure partially enhances the level of total (t-) p70S6K and significantly increases phosphorylated (p-) p70S6K. n = 6; ns: no significance; **: p < 0.01.
Preprints 205133 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

Disclaimer

Terms of Use

Privacy Policy

Privacy Settings

© 2026 MDPI (Basel, Switzerland) unless otherwise stated