<?xml version="1.0" encoding="utf-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Artif. Intell.</journal-id>
<journal-title>Frontiers in Artificial Intelligence</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Artif. Intell.</abbrev-journal-title>
<issn pub-type="epub">2624-8212</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/frai.2024.1376042</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Artificial Intelligence</subject>
<subj-group>
<subject>Technology and Code</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Open and remotely accessible Neuroplatform for research in wetware computing</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name><surname>Jordan</surname> <given-names>Fred D.</given-names></name>
<xref ref-type="corresp" rid="c001"><sup>&#x002A;</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/2519502/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Kutter</surname> <given-names>Martin</given-names></name>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Comby</surname> <given-names>Jean-Marc</given-names></name>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Brozzi</surname> <given-names>Flora</given-names></name>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Kurtys</surname> <given-names>Ewelina</given-names></name>
<uri xlink:href="https://loop.frontiersin.org/people/2638198/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
</contrib-group>
<aff><institution>FinalSpark, Rue du Clos 12</institution>, <addr-line>Vevey</addr-line>, <country>Switzerland</country></aff>
<author-notes>
<fn fn-type="edited-by" id="fn0001">
<p>Edited by: Michael Taynnan Barros, University of Essex, United Kingdom</p>
</fn>
<fn fn-type="edited-by" id="fn0002">
<p>Reviewed by: Chiara Magliaro, University of Pisa, Italy</p>
<p>Qi Huang, Johns Hopkins University, United States</p>
</fn>
<corresp id="c001">&#x002A;Correspondence: Fred D. Jordan, <email>fred.jordan@finalspark.com</email></corresp>
</author-notes>
<pub-date pub-type="epub">
<day>02</day>
<month>05</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>7</volume>
<elocation-id>1376042</elocation-id>
<history>
<date date-type="received">
<day>24</day>
<month>01</month>
<year>2024</year>
</date>
<date date-type="accepted">
<day>11</day>
<month>03</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2024 Jordan, Kutter, Comby, Brozzi and Kurtys.</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Jordan, Kutter, Comby, Brozzi and Kurtys</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>Wetware computing and organoid intelligence is an emerging research field at the intersection of electrophysiology and artificial intelligence. The core concept involves using living neurons to perform computations, similar to how Artificial Neural Networks (ANNs) are used today. However, unlike ANNs, where updating digital tensors (weights) can instantly modify network responses, entirely new methods must be developed for neural networks using biological neurons. Discovering these methods is challenging and requires a system capable of conducting numerous experiments, ideally accessible to researchers worldwide. For this reason, we developed a hardware and software system that allows for electrophysiological experiments on an unmatched scale. The Neuroplatform enables researchers to run experiments on neural organoids with a lifetime of even more than 100 days. To do so, we streamlined the experimental process to quickly produce new organoids, monitor action potentials 24/7, and provide electrical stimulations. We also designed a microfluidic system that allows for fully automated medium flow and change, thus reducing the disruptions by physical interventions in the incubator and ensuring stable environmental conditions. Over the past three years, the Neuroplatform was utilized with over 1,000 brain organoids, enabling the collection of more than 18 terabytes of data. A dedicated Application Programming Interface (API) has been developed to conduct remote research directly via our Python library or using interactive compute such as Jupyter Notebooks. In addition to electrophysiological operations, our API also controls pumps, digital cameras and UV lights for molecule uncaging. This allows for the execution of complex 24/7 experiments, including closed-loop strategies and processing using the latest deep learning or reinforcement learning libraries. Furthermore, the infrastructure supports entirely remote use. Currently in 2024, the system is freely available for research purposes, and numerous research groups have begun using it for their experiments. This article outlines the system&#x2019;s architecture and provides specific examples of experiments and results.</p>
</abstract>
<kwd-group>
<kwd>wetware computing</kwd>
<kwd>organoid intelligence</kwd>
<kwd>biocomputing</kwd>
<kwd>synthetic biology</kwd>
<kwd>AI</kwd>
<kwd>biological neural network</kwd>
<kwd>hybrot</kwd>
</kwd-group>
<counts>
<fig-count count="7"/>
<table-count count="0"/>
<equation-count count="3"/>
<ref-count count="30"/>
<page-count count="16"/>
<word-count count="8370"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Organoid Intelligence</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="sec1">
<label>1</label>
<title>Introduction</title>
<p>The recent rise in wetware computing and consequently, artificial biological neural networks (BNNs), comes at a time when Artificial Neural Networks (ANNs) are more sophisticated than ever.</p>
<p>The latest generation of Large Language Models (LLMs), such as Meta&#x2019;s Llama 2 or OpenAI&#x2019;s GPT-4, fundamentally rely on ANNs.</p>
<p>The recent acceleration of ANN use in everyday life, such as in tools like ChatGPT or Perplexity combined with the explosion in complexity in the underlying ANN&#x2019;s architectures, has had a significant impact on energy consumption. For instance, training a single LLM like GPT-3, a precursor to GPT-4, approximately required 10 GWh, which is about 6,000 times the energy a European citizen uses per year. According to a recent publication the energy consumption projected may increase faster than linearly (<xref ref-type="bibr" rid="ref9">De Vries, 2023</xref>). At the same time, the human brain operates with approximately 86 billion neurons while consuming only 20&#x2009;W of power (<xref ref-type="bibr" rid="ref8">Clark and Sokoloff, 1999</xref>). Given these conditions, the prospect of replacing ANNs running on digital computers with real BNNs is enticing (<xref ref-type="bibr" rid="ref23">Smirnova et al., 2023</xref>). In addition to the substantial energy demands associated with training LLMs, the inference costs present a similarly pressing concern. Recent disclosures reveal that platforms like OpenAI generate over 100 billion words daily through services such as ChatGPT as reported by Sam Altman, the CEO of OpenAI. When we break down these figures, assuming an average of 1.5 tokens per word&#x2014;a conservative estimate based on OpenAI&#x2019;s own tokenizer data&#x2014;the energy footprint becomes staggering. Preliminary calculations, using the LLaMA 65B model (precursor to Llama 2) as a reference point, suggest energy expenditures ranging from 450 to 600 billion Joules per day for word generation alone (<xref ref-type="bibr" rid="ref22">Samsi et al., 2023</xref>). While necessary for providing AI-driven insights and interactions to millions of users worldwide, this magnitude of energy use underscores the urgency for more energy-efficient computing paradigms.</p>
<p>Connecting probes to BNNs is not a new idea. In fact, the field of multi-unit electrophysiology has an established state of the art spanning easily over the past 40&#x2009;years. As a result, there are already well-documented hardware and methods for performing functional electrical interfacing and micro-fluidics needed for nutrient delivery (<xref ref-type="bibr" rid="ref13">Gross et al., 1977</xref>; <xref ref-type="bibr" rid="ref19">Pine, 1980</xref>; <xref ref-type="bibr" rid="ref25">Wagenaar et al., 2005a</xref>; <xref ref-type="bibr" rid="ref16">Newman et al., 2013</xref>). Some systems are also specifically designed for brain organoids (<xref ref-type="bibr" rid="ref30">Yang et al., 2024</xref>). However, their research is mostly focused on exploring brain biology for biomedical applications (e.g., mechanisms and potential treatments of neurodegenerative diseases). The possibility of using these methods for making new computing hardware has not been extensively explored.</p>
<p>For this reason, there is comparatively less literature on methods that can be used to reliably program those BNNs in order to perform specific input&#x2013;output functions (as this is essential for wetware computing, not for biomedical applications). To understand what we need for programming of BNNs, it is helpful to look at analogous problem for ANNs.</p>
<p>For ANNs, the programming task involves finding the network parameters, globally denoted as <inline-formula>
<mml:math id="M1">
<mml:mi>S</mml:mi>
</mml:math>
</inline-formula> below, that minimize the difference <inline-formula>
<mml:math id="M2">
<mml:mi>L</mml:mi>
</mml:math>
</inline-formula> computed between expected output <inline-formula>
<mml:math id="M3">
<mml:mi>E</mml:mi>
</mml:math>
</inline-formula> and actual output <inline-formula>
<mml:math id="M4">
<mml:mi>O</mml:mi>
</mml:math>
</inline-formula>, for given inputs <inline-formula>
<mml:math id="M5">
<mml:mi>I</mml:mi>
</mml:math>
</inline-formula>, given the transfer function <inline-formula>
<mml:math id="M6">
<mml:mi>T</mml:mi>
</mml:math>
</inline-formula> of the ANN. This can be written as:</p>
<p><inline-formula>
<mml:math id="M7">
<mml:mi>L</mml:mi>
<mml:mo>=</mml:mo>
<mml:mi>f</mml:mi>
<mml:mfenced open="(" close=")" separators=",">
<mml:mi>O</mml:mi>
<mml:mi>E</mml:mi>
</mml:mfenced>
</mml:math>
</inline-formula>, with <inline-formula>
<mml:math id="M8">
<mml:mi>O</mml:mi>
<mml:mo>=</mml:mo>
<mml:mi>T</mml:mi>
<mml:mfenced open="(" close=")" separators=",">
<mml:mi>I</mml:mi>
<mml:mi>S</mml:mi>
</mml:mfenced>
</mml:math>
</inline-formula></p>
<p>where <inline-formula>
<mml:math id="M9">
<mml:mi>f</mml:mi>
</mml:math>
</inline-formula> is typically a function that equals 0 when <inline-formula>
<mml:math id="M10">
<mml:mi>O</mml:mi>
<mml:mo>=</mml:mo>
<mml:mi>E</mml:mi>
</mml:math>
</inline-formula>.</p>
<p>The same equation applies to BNNs. However, the key differences compared to ANNs include the fact that the network parameters <inline-formula>
<mml:math id="M11">
<mml:mi>S</mml:mi>
</mml:math>
</inline-formula> cannot be individually adjusted in the case of BNNs, and the transfer function <inline-formula>
<mml:math id="M12">
<mml:mi>T</mml:mi>
</mml:math>
</inline-formula> is both unknown and non-stationary. Therefore, alternative heuristics must be developed, for instance based on spatiotemporal stimulation patterns (<xref ref-type="bibr" rid="ref2">Bakkum et al., 2008</xref>; <xref ref-type="bibr" rid="ref14">Kagan et al., 2022</xref>; <xref ref-type="bibr" rid="ref4">Cai et al., 2023a,b</xref>). Such developments necessitate numerous electrophysiological experiments, including, for instance, complex closed-loop algorithms where stimulation is a function of the network&#x2019;s prior responses. These experiments can sometimes span days or months.</p>
<p>To facilitate long-term experiments involving a global network of research groups, we designed an open innovation platform. This platform enables researchers to remotely perform experiments on a server interfaced with our hardware. For example, our Neuroplatform enhances the chances of discovering the abovementioned stimulation heuristics. It should be noted that, outside of the field of neuroplasticity, similar open platforms were already proposed in 2023 (<xref ref-type="bibr" rid="ref17">O&#x2019;Leary et al., 2022</xref>; <xref ref-type="bibr" rid="ref1">Armer et al., 2023</xref>; <xref ref-type="bibr" rid="ref10">Elliott et al., 2023</xref>; <xref ref-type="bibr" rid="ref31">Zhang et al., 2023</xref>). However, to our knowledge, there are no platforms specifically dedicated to research related to biocomputing.</p>
</sec>
<sec id="sec2">
<label>2</label>
<title>Biological setup</title>
<p>The biological material used in our platform is made of brain spheroids [also called minibrains (<xref ref-type="bibr" rid="ref12">Govindan et al., 2021</xref>), brain organoids (<xref ref-type="bibr" rid="ref20">Qian et al., 2019</xref>), or neurospheres (<xref ref-type="bibr" rid="ref3">Brewer and Torricelli, 2007</xref>)] developed from Human iPSC-derived Neural Stem Cells (NSCs), following the protocol of Prof. Roux Lab (<xref ref-type="bibr" rid="ref12">Govindan et al., 2021</xref>). Based on the recent guidelines to clarify the nomenclature for defining 3D cellular models of the nervous system (<xref ref-type="bibr" rid="ref18">Pa&#x015F;ca et al., 2022</xref>), we can call those brain spheroids &#x201C;forebrain organoids&#x201D; (FOs). Generation of brain organoids from NSCs has been already described for both mouse (<xref ref-type="bibr" rid="ref7">Ciarpella et al., 2023</xref>), and human models (<xref ref-type="bibr" rid="ref15">Lee et al., 2020</xref>). Our protocol is based on the following steps: expansion phase of the NSCs, induction of the 3D structure, differentiation steps (using GDNF and BDNF), and maturation phase (<xref ref-type="fig" rid="fig1">Figures 1A</xref>,<xref ref-type="fig" rid="fig1">B</xref>). The <xref ref-type="fig" rid="fig1">Figure 1C</xref> is an image of the FO obtained using electronic microscope, it shows that it is a compact spheroid. The average shape of FOs obtained with this protocol are spheroids of a diameter around 500&#x2009;&#x03BC;m (<xref ref-type="bibr" rid="ref12">Govindan et al., 2021</xref>). Our experiments show that the FOs obtained can be kept alive in an orbital shaker for years, as previously demonstrated (<xref ref-type="bibr" rid="ref12">Govindan et al., 2021</xref>).</p>
<fig position="float" id="fig1">
<label>Figure 1</label>
<caption>
<p>FO generation and MEA setup. <bold>(A)</bold> Protocol used for the generation of forebrain organoids (FO). Neural progenitors are first thawed, plated and expanded in T25 flasks. They are then differentiated in P6 dishes on orbital shakers, and finally manually placed on the MEA. <bold>(B)</bold> Representative images showing various stages of FO formation and differentiation, taken at different time points. The scale bar represents 250&#x2009;&#x03BC;m. <bold>(C)</bold> Image of a whole FO taken with scanning electron microscope. The scale bar represents 100&#x2009;&#x03BC;m. <bold>(D)</bold> Microscope view of the FO (in white) sitting on the electrodes of the MEA, and the membrane. The hole in the membrane is not visible on the picture since it is hidden by the FO. The scale bar represents 500&#x2009;&#x03BC;m <bold>(E)</bold> Overview of the MEA, where the 32 electrodes are visible as 4 sets of 8 electrodes each. An FO is placed atop of each set of 8 electrodes, visible as a darker area. For each FO, the 2 circles correspond to a 2.5&#x2009;mm circular membrane with a central hole. The scale bar represents 1&#x2009;mm. <bold>(F)</bold> Cross-sectional view of the MEA setup, illustrating the air-liquid interface. The medium covering the FO is supplied from the medium chamber through the porous membranes.</p>
</caption>
<graphic xlink:href="frai-07-1376042-g001.tif"/>
</fig>
<p>Gene expression analysis of mature FOs vs. NSCs showed a marked upregulation of genes characteristic to neurons, oligodendrocytes and astrocytes in FOs compare to NSCs. More precisely, FOs expressed genes typically enriched in the forebrain, such as striatum, sub pallium, and layer 6 of motor cortex (<xref ref-type="bibr" rid="ref12">Govindan et al., 2021</xref>). Pathway enrichment analysis of FOs vs. NSCs demonstrated activation of biological processes like synaptic activity, neuron differentiation and neurotransmitter release (<xref ref-type="bibr" rid="ref12">Govindan et al., 2021</xref>).</p>
<p>At the age of 12&#x2009;weeks, FOs contain a high number of ramified neurons (<xref ref-type="bibr" rid="ref12">Govindan et al., 2021</xref>), and they are mature enough to be transferred to the electrophysiological measurement system (<xref ref-type="fig" rid="fig1">Figure 1A</xref>). In this setup, they have a life expectancy of several months, even with 24/7 experiments that include hours of electrical stimulations. This setup has a quick turnaround with occasional downtime &#x2013; about 1&#x2009;h &#x2013; during organoid replacements. Therefore, the platform maintains a high availability for experiments.</p>
</sec>
<sec id="sec3">
<label>3</label>
<title>Hardware architecture</title>
<sec id="sec4">
<label>3.1</label>
<title>Introduction</title>
<p>The remotely accessible hardware includes all the systems which are required to preserve homeostasis, monitor environmental parameters and perform electrophysiological experiments. These systems can be controlled interactively using our custom Graphical User Interface (GUI) or via Python scripts. All data is stored in a time-series database (InfluxDB), which can be accessed either via a GUI or via Python scripts. The users typically connect to the system using the Remote Desktop Protocol (RDP).</p>
<p>The platform is composed of several sub-systems, which can be accessed remotely via API calls over the internet, typically through Python.</p>
</sec>
<sec id="sec5">
<label>3.2</label>
<title>Multi-Electrode Array (MEA)</title>
<p>Our current platform features 4 MEAs. The MEAs were designed by Prof. Roux&#x2019;s Lab form Haute Ecole du Paysage, d&#x2019;Ing&#x00E9;nierie et d&#x2019;Architecture (HEPIA) and are described in <xref ref-type="bibr" rid="ref28">Wertenbroek et al. (2021)</xref>. Each MEA can accommodate 4 organoids, with 8 electrodes per organoid (<xref ref-type="fig" rid="fig1">Figure 1E</xref>).</p>
<p>The MEA setup utilizes an Air-Liquid-Interface (ALI) approach (<xref ref-type="bibr" rid="ref24">Stoppini et al., 1991</xref>), in which the organoids are directly placed on electrodes located atop of a permeable membrane (<xref ref-type="fig" rid="fig1">Figure 1D</xref>), with the medium flowing beneath this membrane in a 170&#x2009;&#x03BC;L chamber. As a result, a thin layer of medium, created by surface tension, separates the upper side of the organoids from the humidified incubator air. This arrangement is further protected by a lid partially covering the MEA (<xref ref-type="fig" rid="fig1">Figure 1F</xref>). This ALI method enables a higher throughput and higher stability compared to submerged approaches, since no dedicated coating is required, and it is less prone to have the organoids detaching from the electrodes.</p>
</sec>
<sec id="sec6">
<label>3.3</label>
<title>Electrophysiological stimulation and recording system</title>
<p>The electrodes in our system enable both stimulation and recording. The respective digital-to-analog and analog-to-digital conversions are performed by Intan RHS 32 headstages. Stimulations are executed using a current controller that ranges from 10&#x2009;nA to 2.5&#x2009;mA, and recordings are obtained by measuring the voltage on each electrode at a 30&#x2009;kHz sampling frequency with a 16 bits resolution giving an accuracy of 0.15&#x2009;&#x03BC;V. The headstages are connected to an Intan RHS controller, which in turn is connected to a computer via a USB port. The <xref ref-type="fig" rid="fig2">Figure 2A</xref> shows the electrical activity recorded for each of the 32 electrodes. It can be noticed that the recorded activity is different between each electrode. This difference comes from the facts that each set of 8 electrodes records a different FO and that for a given FO, electrodes record at a different location. This display is refreshed in real-time and also available 24/7 on our website at the URL <ext-link xlink:href="https://finalspark.com/live/" ext-link-type="uri">https://finalspark.com/live/</ext-link>. We compared the recording characteristics of this ALI setup to MCS MEA (60MEA200/30iR-Ti) monitoring a submerged FO, using the exact same Intan system for voltage conversion. The overlays of an action potential recorded, respectively, with the ALI and submerged versions are shown in <xref ref-type="fig" rid="fig2">Figures 2C</xref>,<xref ref-type="fig" rid="fig2">D</xref> and show similar signal characteristics.</p>
<fig position="float" id="fig2">
<label>Figure 2</label>
<caption>
<p>Recording system and user interface. <bold>(A)</bold> Electrical activity measured in &#x03BC;V over one second for each of the 32 electrodes. Each set of 8 electrodes records a different FO. <bold>(B)</bold> Graphical User Interface for manually controlling each of the microfluidic pumps. <bold>(C)</bold> Overlays of FO action potential recorded by the ALI system of the Neuroplatform. <bold>(D)</bold> Overlays of FO action potentials recorded with an MCS system. <bold>(E)</bold> Fluctuations of the flowrate of the medium within the microfluidic system, illustrating the cyclic variations induced by the peristaltic pump operating at 1 round per minute with 10 cams. <bold>(F)</bold> Temporal variations of the red component of the medium color, triggered by a sudden change in medium acidity, resulting in phenol red color change.</p>
</caption>
<graphic xlink:href="frai-07-1376042-g002.tif"/>
</fig>
</sec>
<sec id="sec7">
<label>3.4</label>
<title>Micro-fluidics</title>
<p>To sustain the life of the organoids on the MEA, Neuronal Medium (NM) needs to be constantly supplied. Our Neuroplatform is equipped with a closed-loop microfluidic system that allows for a 24/7 medium supply. The medium is circulating at a rate of 15uL/min. The medium flow rate is controlled by a BT-100 2&#x2009;J peristaltic pump and is continuously adjusted according to needs, for instance during experimental runs. The peristaltic pump is connected to the PC-control software using an RS485 interface, for programmed (i.e., in Python) or manual operations (<xref ref-type="fig" rid="fig2">Figure 2B</xref>). Additionally, <xref ref-type="fig" rid="fig3">Figure 3A</xref> depicts this microfluidic closed-loop circuit.</p>
<fig position="float" id="fig3">
<label>Figure 3</label>
<caption>
<p>Microfluidics. <bold>(A)</bold> Microfluidic system illustrating the continuously operating primary system, which ensures constant flow in the medium chamber, and the secondary system responsible for medium replacing every 48&#x2009;h. <bold>(B)</bold> Side view of the assembly, featuring the camera and the MEA. The entire assembly is enclosed with aluminum foil to ensure the lowest possible noise level. <bold>(C)</bold> Front view of the assembly, showing the intake and outtake of the microfluidic system, as well as the LED used during image capture.</p>
</caption>
<graphic xlink:href="frai-07-1376042-g003.tif"/>
</fig>
<p>The microfluidic circuit is made of 0.8&#x2009;mm (inside diameter, ID) tubing. Continuous monitoring of the microfluidic circuit and flow rate is achieved by using Fluigent flow-rate sensors, which connect to the Neuroplatform control center via USB. Data related to medium flow rate is stored in a database for later access. <xref ref-type="fig" rid="fig2">Figure 2E</xref> shows the cyclic variations in flow induced by the cams of the peristaltic pump.</p>
<p>A secondary microfluidic system is used to replace the medium in the closed-loop with fresh medium every 24&#x2009;h, a process illustrated in <xref ref-type="fig" rid="fig3">Figure 3A</xref>. This replacement is fully automated through a Python script and performed in the following consecutive steps:</p>
<list list-type="order">
<list-item>
<p>Set the rotary valve to select the path from the reservoir F50 to the syringe pump</p>
</list-item>
<list-item>
<p>Pump 2&#x2009;mL of old medium using the syringe pump</p>
</list-item>
<list-item>
<p>Set the rotary valve to select the path from the syringe pump to the waste F50</p>
</list-item>
<list-item>
<p>Push 2&#x2009;mL of old medium to the waste using the syringe pump</p>
</list-item>
<list-item>
<p>Set the rotary valve to select the path from the new medium in the F50 in the fridge to the syringe pump</p>
</list-item>
<list-item>
<p>Pump 2&#x2009;mL of fresh medium using the syringe pump</p>
</list-item>
<list-item>
<p>Set the rotary valve to select the path from the syringe pump to the reservoir F50</p>
</list-item>
<list-item>
<p>Push 2&#x2009;mL of fresh medium using the syringe pump</p>
</list-item>
</list>
</sec>
<sec id="sec8">
<label>3.5</label>
<title>Cameras</title>
<p>Each MEA is equipped with a 12.3-megapixel camera that can be controlled interactively or programmatically (i.e., through a Raspberry Pi) for still image capture or video recording. The camera is positioned below the MEA, while illumination is provided by a remotely controlled LED situated above the MEA. <xref ref-type="fig" rid="fig3">Figures 3B</xref>,<xref ref-type="fig" rid="fig3">C</xref> illustrate this assembly (the aluminum wrapping is used in order to minimize the noise). This setup is particularly useful for detecting various changes, such as cell necrosis, possible organoid displacement caused by microfluidics, variations in medium acidity (using color analysis since our medium contains Phenol red), contamination, neuromelanin production (which can happen when uncaging dopamine), overflows (where the medium inadvertently fills the chamber above the membrane), or bubbles in the medium. For the latter two events, dedicated algorithms automatically detect these issues and trigger an alert to the on-site operator.</p>
<p>Changes of acidity, for example, can be detected by measuring the average color over a pre-defined window. <xref ref-type="fig" rid="fig2">Figure 2F</xref> shows the evolution of the medium&#x2019;s red color component, with data points recorded hourly. The noticeable sudden drop is attributed to the pumping of medium with a slightly different acidity. This change in acidity results in a color alteration of the phenol red present in the medium.</p>
</sec>
<sec id="sec9">
<label>3.6</label>
<title>UV light controlled uncaging</title>
<p>It is also possible to release molecules at specific timings using a process called uncaging. In this method, a specific wavelength of light is employed to break open a molecular &#x201C;cage&#x201D; that contains a neuroactive molecule, such as Glutamate, NMDA or Dopamine. A fiber optic of 1,500&#x2009;&#x03BC;m core diameter and a numerical aperture of 0.5 is used to direct light in the medium within the MEA chamber. The current system, Prizmatix Silver-LED, operates at 365&#x2009;nm with an optical power of 260&#x2009;mW. The uncaging system is fully integrated into the Neuroplatform and can be programmatically controlled during experiment runs via our API (see section 5.3).</p>
</sec>
<sec id="sec10">
<label>3.7</label>
<title>Environmental measurements</title>
<p>The environmental conditions are monitored within two incubators. In both incubators, the following parameters are recorded: CO2, O2 concentrations, humidity, atmospheric pressure and temperature. Door-opening events are also logged since they have a major impact on measurements. The primary purpose of this monitoring is to ensure that experiments are performed in stable and reproducible environmental conditions.</p>
<p>All these parameters are displayed in real-time in a graphic interface showing both instant values as well as variations versus time of noise and flowrates (<xref ref-type="fig" rid="fig4">Figure 4A</xref>).</p>
<fig position="float" id="fig4">
<label>Figure 4</label>
<caption>
<p>Graphic user interface to monitor critical parameters in the incubators. <bold>(A)</bold> Graphical User Interface displaying critical environmental conditions for the incubator 1, where electrophysiological experiments are performed, as well as the incubator 2, where FO are maintained on an orbital shaker. <bold>(B)</bold> The display shows environmental data for incubator 1 for specific time periods, extracted from the database, with door opening events displayed as dashed line. Noise, Temperature, humidity and pressure are indicated by different colored lines. The units of each measurement are normalized between 0 and 1 for the selected time interval.</p>
</caption>
<graphic xlink:href="frai-07-1376042-g004.tif"/>
</fig>
<p>Incubator 1 houses the MEAs and the organoids used for electrophysiological experiments. In addition to the mentioned parameters, flowmeters are also utilized to report the actual flow rate of the microfluidic for each MEA, as depicted in the graph labelled &#x201C;Pump&#x201D; in <xref ref-type="fig" rid="fig4">Figure 4A</xref>. The system&#x2019;s state is indirectly monitored through the noise level of each MEA, as shown in the graph labelled &#x201C;Noise Intan&#x201D; in <xref ref-type="fig" rid="fig4">Figure 4A</xref>. The noise level is calculated based on the standard deviation of the electrical signals recorded by the electrodes over a 30&#x2009;ms period.</p>
<p>Incubator 2 houses the organoids which are kept in orbital shakers. Piezoelectric gyroscopes are used to measure the actual rotation speed of the orbital shakers.</p>
<p>Since all the data is logged in the database, it is also possible to access the historical measurements through a dedicated GUI (<xref ref-type="fig" rid="fig4">Figure 4B</xref>).</p>
</sec>
</sec>
<sec id="sec11">
<label>4</label>
<title>Software</title>
<sec id="sec12">
<label>4.1</label>
<title>General architecture</title>
<p>The core of the system relies on a computational notebook which provides access to 3 resources (<xref ref-type="fig" rid="fig5">Figure 5A</xref>):</p>
<list list-type="simple">
<list-item>
<p>1. A database where all the information regarding the Neuroplatform system is stored</p>
</list-item>
<list-item>
<p>2. The Intan software running on a dedicated PC, which is used for:</p>
</list-item>
</list>
<list list-type="bullet">
<list-item>
<p>Recording the number of detected spikes in a 200&#x2009;ms time window</p>
</list-item>
<list-item>
<p>Setting stimulation parameters</p>
</list-item>
</list>
<list list-type="simple">
<list-item>
<p>3. A Raspberry Pi for triggering current stimulation according to stimulation parameters</p>
</list-item>
</list>
<fig position="float" id="fig5">
<label>Figure 5</label>
<caption>
<p>Software setup and electrical stimulation. <bold>(A)</bold> General architecture of the Neuroplatform. The Jupyter Notebook serves as the main controller, enabling initiation and reading of spikes, configuration stimulation signals and access to database via, e.g., Python <bold>(B)</bold> Parameters of the stimulation current: settings optimally these parameters can elicit spikes. Through the Python API, parameters that can be adjusted for the bi-phasic stimulation signals include the duration (D1) and amplitude (A1) of the positive current phase, and, respectively, D2 and A2 for the negative current phase. Additionally, the polarity of the biphasic signal can be reversed to start with a negative current.</p>
</caption>
<graphic xlink:href="frai-07-1376042-g005.tif"/>
</fig>
</sec>
<sec id="sec13">
<label>4.2</label>
<title>Database</title>
<p>The Neuroplatform records monitored data 24/7 using InfluxDB, a database designed for time series. Other options are also available.</p>
<p>This database contains all the data coming from the hardware listed in Section 3.</p>
<p>The electrical activity of the neurons is also recorded 24/7 at a sampling rate of 30&#x2009;kHz. To minimize the volume of stored data, we designed a dedicated process that focuses on significant events, such as threshold crossings that are likely to be due to action potentials (spikes). The following pseudo code illustrates the implemented approach:</p>
<list list-type="simple">
<list-item>
<p>- Each 1-min write buffer to database</p>
</list-item>
<list-item>
<p>- Each 33&#x2009;&#x03BC;s</p>
</list-item>
<list-item>
<p>- For each electrode</p>
</list-item>
<list-item>
<p>- If, at time <italic>t</italic>, the voltage exceeds a threshold <italic>T</italic></p>
</list-item>
<list-item>
<p>- Store (in buffer) 3&#x2009;ms of data [<italic>t</italic>-1&#x2009;ms, <italic>t</italic>&#x2009;+&#x2009;2&#x2009;ms]</p>
</list-item>
<list-item>
<p>- Each 3&#x2009;s update T</p>
</list-item>
</list>
<p>Additionally, a timestamp corresponding to each detected event is also stored in the database, along with the maximum value of voltage during the 3&#x2009;ms spike waveform recording.</p>
<p>The threshold <italic>T</italic> is computed directly from voltage values sampled each 33&#x2009;&#x03BC;s, according to the following formula:</p>
<disp-formula id="E1">
<mml:math id="M13">
<mml:mi>T</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>6</mml:mn>
<mml:mo>&#x2217;</mml:mo>
<mml:mi>M</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>n</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mfenced open="{" close="}">
<mml:msub>
<mml:mi>&#x03C3;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mfenced>
</mml:mfenced>
</mml:math>
</disp-formula>
<p>Where <inline-formula>
<mml:math id="M14">
<mml:msub>
<mml:mi>&#x03C3;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:math>
</inline-formula> is the standard deviation computed over a set <italic>i</italic> of 30&#x2009;ms consecutive voltage values, and <inline-formula>
<mml:math id="M15">
<mml:mi>M</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>n</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow/>
</mml:mfenced>
</mml:math>
</inline-formula> represents the median function computed over 101 consecutive <inline-formula>
<mml:math id="M16">
<mml:mfenced open="{" close="}">
<mml:msub>
<mml:mi>&#x03C3;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mfenced>
</mml:math>
</inline-formula> values. The use of the median reduces the sensitivity to outliers, which is typically caused by action potentials. In our current setup, a multiplier of 6 on the median has proven to be a good compromise for achieving reliable spike detection.</p>
<p>Besides electric tension data, the number spikes recorded per minute is also computed and stored in the database every minute by a batch process.</p>
</sec>
<sec id="sec14">
<label>4.3</label>
<title>Recording electrical activity</title>
<p>As previously discussed, the communication among neurons is captured by the MEA and converted into a voltage signal sampled at 30&#x2009;kHz. The Neuroplatform offers two basic access modes to the recorded neural activity:</p>
<list list-type="order">
<list-item>
<p>Raw: raw sampling values.</p>
</list-item>
<list-item>
<p>Optimized: waveforms of the raw signal near neuronal spikes, available directly from the database.</p>
</list-item>
</list>
<p>In addition to the aforementioned features, the Neuroplatform offers even more advanced methods. For instance, it includes counting spikes over a fixed time period of 200&#x2009;ms following stimulation, with a 10&#x2009;ms delay suppressing the stimulation artifact.</p>
<p>From a technical perspective, accessing the number of spikes can be accomplished in two different ways:</p>
<list list-type="simple">
<list-item>
<p>- Retrieving the number of spikes per minute from the database</p>
</list-item>
<list-item>
<p>- Through direct communication with the PC managing the Intan controller for spike count</p>
</list-item>
</list>
<p>The second approach is required when the stimulation protocol demands real-time responsiveness. This is typically the case for certain closed-loop strategies. For instance, closed-loop stimulation strategies have been deployed in primary cortical cultures for effective burst control (<xref ref-type="bibr" rid="ref25">Wagenaar et al., 2005a</xref>,<xref ref-type="bibr" rid="ref26">b</xref>) and for goal-directed learning (<xref ref-type="bibr" rid="ref22">Samsi et al., 2023</xref>).</p>
</sec>
<sec id="sec15">
<label>4.4</label>
<title>Syntax for stimulations</title>
<p>Programmatically stimulating the FO on the Neuroplatform is accomplished by sending an electrical current to the MEA electrodes. The electrical current profile can be parameterized in a variety of ways, which is partly shown in <xref ref-type="fig" rid="fig5">Figure 5B</xref>. These parameters and controls include:</p>
<list list-type="simple">
<list-item>
<p>- Basic shape of stimulation signal:</p>
</list-item>
</list>
<list list-type="simple">
<list-item>
<p>o Bi-phasic</p>
</list-item>
<list-item>
<p>o Bi-phasic with interphase delay</p>
</list-item>
<list-item>
<p>o Tri-phasic</p>
</list-item>
</list>
<list list-type="simple">
<list-item>
<p>- Stimulation duration and intensity:</p>
</list-item>
</list>
<list list-type="simple">
<list-item>
<p>o Positive (A1) and negative (A2) electrical current intensity (typical 1uA, ranging from 0.1uA to 20uA)</p>
</list-item>
<list-item>
<p>o Duration of positive (D1) and negative (D2) stimulation currents</p>
</list-item>
</list>
<list list-type="simple">
<list-item>
<p>- Stimulation triggers</p>
</list-item>
</list>
<list list-type="simple">
<list-item>
<p>o Single start</p>
</list-item>
<list-item>
<p>o Table with collection of start triggers</p>
</list-item>
<list-item>
<p>o Pulse trains</p>
</list-item>
</list>
<list list-type="simple">
<list-item>
<p>- MEA electrodes</p>
</list-item>
</list>
<p>send_stim_param (electrodes, params)</p>
</sec>
</sec>
<sec>
<label>5</label>
<title>Examples of electrophysiological experiments</title>
<p>To demonstrate the effectiveness of the Neuroplatform, the following sections will provide an overview of several experiments conducted on the Neuroplatform at FinalSpark&#x2019;s Laboratories in Vevey, Switzerland.</p>
<sec>
<label>5.1</label>
<title>Modification of spontaneous activity</title>
<p>The spontaneous electrical activity of the FO can be represented by the concept of &#x201C;Center of Activity&#x201D; (CA) (<xref ref-type="bibr" rid="ref2">Bakkum et al., 2008</xref>) which is defined as a virtual position <inline-formula>
<mml:math id="M17">
<mml:mi>C</mml:mi>
</mml:math>
</inline-formula> on the MEA described by:</p>
<disp-formula id="E2">
<mml:math id="M18">
<mml:mi>C</mml:mi>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msubsup>
<mml:mstyle displaystyle="true">
<mml:mo stretchy="true">&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mn>8</mml:mn>
</mml:msubsup>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:mo>&#x22C5;</mml:mo>
<mml:mfenced open="(" close=")" separators=",">
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>Y</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:msubsup>
<mml:mstyle displaystyle="true">
<mml:mo stretchy="true">&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mn>8</mml:mn>
</mml:msubsup>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:math>
</disp-formula>
<p>Where <inline-formula>
<mml:math id="M19">
<mml:mfenced open="(" close=")" separators=",">
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>Y</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mfenced>
</mml:math>
</inline-formula> define the spatial position of the 8 electrodes and <inline-formula>
<mml:math id="M20">
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:math>
</inline-formula>is the number of spontaneous spikes detected. The interest of the concept of CA is that its position provides statistical information about the average location of the activity over the surface of the FO. The ability to change the position of the CA is interesting because it also shows the ability to memorize information in the state of the FO.</p>
<p>The coordinates of the CA can be modified using a high frequency stimulation. In the following experiment we use the following protocol:</p>
<list list-type="simple">
<list-item>
<p>1) Compute the CA using the number of detected spikes over 500&#x2009;ms</p>
</list-item>
<list-item>
<p>2) Goto 1,100x</p>
</list-item>
<list-item>
<p>3) Perform a 20&#x2009;Hz stimulation during 500&#x2009;ms using a bi-phasic current (negative first) of 2&#x2009;&#x03BC;A of 200 &#x03BC;S, for both phases, on one electrode</p>
</list-item>
<list-item>
<p>4) Wait 1&#x2009;s</p>
</list-item>
<list-item>
<p>5) Goto 5,100x</p>
</list-item>
</list>
<p><xref ref-type="fig" rid="fig6">Figure 6A</xref> displays the 100 measured positions of the CA corresponding to the spontaneous activity before the 20&#x2009;Hz stimulation in blue, and after the high-frequency stimulation in red (the average position is indicated by a cross). A close-up is shown in <xref ref-type="fig" rid="fig6">Figure 6B</xref>. The timestamps of the spontaneous activity, before and after stimulation, are presented in <xref ref-type="fig" rid="fig6">Figures 6C</xref>,<xref ref-type="fig" rid="fig6">D</xref>, respectively. Each graph shows one example of the 100 records of 500&#x2009;ms used to compute the CA location (showing a decrease of spontaneous firing activity of electrodes 3, 4 and 6). A noticeable shift in the average position (shown by a cross) of the CA can be observed before and after the high-frequency stimulation (as seen in <xref ref-type="fig" rid="fig6">Figure 6A</xref>), indicating a change of state of the biological network. A classifier based on a simple logistic regression is employed to predict if the network has received the 20&#x2009;Hz stimulation. In this particular experiment, the classification accuracy, computed from the confusion matrix, is 95.5%.</p>
<fig position="float" id="fig6">
<label>Figure 6</label>
<caption>
<p>Center of activity modification. <bold>(A)</bold> Graph showing the 2D layout of the 8 electrodes, the X and Y axis are normalized units showing the spatial coordinates of the electrodes. All electrodes can be used for both stimulation and reading. A 20&#x2009;Hz stimulation signal is applied to electrode 6. The 100 blue circles represent the positions of the Center of Activity (CA) before 20&#x2009;Hz stimulation, while the 100 red circles indicate the positions after the stimulation. The cross mark the average position. <bold>(B)</bold> A closer look at the two groups of CA. <bold>(C)</bold> Timestamps depicting the spontaneous activity over 500&#x2009;ms for each of the 8 electrodes before the high-frequency stimulation. <bold>(D)</bold> Spontaneous activity observed after the high-frequency stimulation, showing a lower activity of electrodes 6, 4 and 3, compared to <bold>(C)</bold>.</p>
</caption>
<graphic xlink:href="frai-07-1376042-g006.tif"/>
</fig>
<p>The Neuroplatform allows users to perform both the experimental part (including stimulation and reading operations) and the visualization of the CA displacement within the same Python source code. The 500&#x2009;ms 20&#x2009;Hz signal is generated directly by the Python source code shown below. The first trigger.send instruction sends the trigger for the stimulation on a specific electrode and time.sleep pauses the execution for 50&#x2009;ms.<inline-graphic xlink:href="frai-07-1376042-i001.tif"/></p>
<p>Despite the common perception of Python as being less than ideal for real-time signal processing due to its inherent latency, our empirical data reveals a time accuracy of under 1&#x2009;ms (on an Intel Xeon CPU E5-2690 v2 @ 3.00GHz), a level of precision that is satisfactory for the generation of tetanic signals.</p>
</sec>
<sec id="sec18">
<label>5.2</label>
<title>Optimization of stimulation parameters</title>
<p>In this example, the objective is to identify the set of stimulation parameters that can elicit the maximum number of action potentials within 200&#x2009;ms after a stimulation.</p>
<p>Depending on the FOs, their composition, and maturity, only specific combinations of electrodes and parameters can elicit spikes. In our experiment, we use an 8-electrode MEA and cycle through several stimulation signal parameters as shown in <xref ref-type="fig" rid="fig7">Figure 7A</xref>. Consequently, we need to test a total of 342 different parameter-electrode combinations. The following pseudo code illustrates the Python script used in this experiment.</p>
<list list-type="simple">
<list-item>
<p>1) For each set of stimulation parameters</p>
</list-item>
<list-item>
<p>2) For each stimulation electrode</p>
</list-item>
<list-item>
<p>3) For each recording electrode</p>
</list-item>
<list-item>
<p>4) During 15&#x2009;s, every 250&#x2009;ms</p>
</list-item>
<list-item>
<p>5) Decide between stimulating, or recording spontaneous activity, with a 50% probability</p>
</list-item>
<list-item>
<p>6) Record number of spikes during 200&#x2009;ms</p>
</list-item>
</list>
<fig position="float" id="fig7">
<label>Figure 7</label>
<caption>
<p>Neural activity stimulation and dopamine uncaging. <bold>(A)</bold> Graph depicting the number of spikes recorded over 250&#x2009;ms. The spike counts in orange were measured following a stimulation, while those in blue were measured during periods without stimulation. For clarity in visualization, a small bar is displayed even when no spikes are detected. <bold>(B)</bold> Diagram illustrating the different steps involved in the closed-loop uncaging process of dopamine, which is repeated 240 times. <bold>(C)</bold> Timestamps of action potentials from the 8 electrodes before and after stimulation (shown as red line), showcasing the elicited spikes. <bold>(D)</bold> Graph displaying the number of elicited spikes over the 240 steps of the closed-loop (in blue) alongside the activation events of the UV light source (red).</p>
</caption>
<graphic xlink:href="frai-07-1376042-g007.tif"/>
</fig>
<p>The aim of probabilistic stimulation and no stimulation in step 5 is to evaluate the difference between elicited and spontaneous spikes in a way that ensures there is no bias.</p>
<p>The bar chart in <xref ref-type="fig" rid="fig7">Figure 7A</xref> displays a segment of the experimental results. It shows a 15-s recording from a single electrode, corresponding to one execution of step 4 in the pseudo code above. Each bar represents the spike count during a 200&#x2009;ms period, repeated every 250&#x2009;ms. The orange bars in this plot are the result of the parameters selected in step 1 of the pseudo code. The blue bars represent no-stimulation periods, thus corresponding to the spontaneous activity of the neurons.</p>
<p>From <xref ref-type="fig" rid="fig7">Figure 7A</xref>, we can see that this particular combination of electrode and parameters reliably elicits responses.</p>
<p>In practice, the Python script can also be used to automatically display the 342 graphs similar to <xref ref-type="fig" rid="fig7">Figure 7A</xref>, allowing the operator to select the optimal set of parameters. Additionally, it can compute a scalar metric to characterize the &#x201C;efficiency&#x201D; of the parameters, and automatically identify the optimal parameters.</p>
<p>An example of a parameter maximization metric is given in the equation below. Let us denote <inline-formula>
<mml:math id="M21">
<mml:msub>
<mml:mi>&#x03BC;</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
</mml:math>
</inline-formula> and <inline-formula>
<mml:math id="M22">
<mml:msub>
<mml:mi>&#x03BC;</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
</mml:math>
</inline-formula> the average number of spikes recorded spontaneously or after a stimulation, respectively, and <inline-formula>
<mml:math id="M23">
<mml:msub>
<mml:mi>&#x03C3;</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
</mml:math>
</inline-formula> and <inline-formula>
<mml:math id="M24">
<mml:msub>
<mml:mi>&#x03C3;</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
</mml:math>
</inline-formula> as their standard deviations. The following metric is used:</p>
<disp-formula id="E3">
<mml:math id="M25">
<mml:mi>m</mml:mi>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mfenced open="|" close="|">
<mml:mrow>
<mml:msub>
<mml:mi>&#x03BC;</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x03BC;</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
<mml:mrow>
<mml:mo>max</mml:mo>
<mml:mfenced open="(" close=")" separators=",">
<mml:msub>
<mml:mi>&#x03C3;</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>&#x03C3;</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
</mml:mfenced>
</mml:mrow>
</mml:mfrac>
</mml:math>
</disp-formula>
<p>The set of parameters that maximize this metric can then be utilized to perform other experiments requiring elicited spikes, such as investigating the effect of pharmacological agents on a biological network&#x2019;s ability to react quickly to stimulation.</p>
</sec>
<sec id="sec19">
<label>5.3</label>
<title>UV light-induced uncaging of molecules</title>
<p>&#x2018;Uncaging&#x2019; is a pivotal technique in cellular biology, enabling the precise control of molecular interactions within cells (<xref ref-type="bibr" rid="ref11">Gienger et al., 2020</xref>). It involves the use of photolabile caged compounds that are activated by specific light wavelengths, releasing bioactive molecules in a targeted and timely manner. This method is particularly valuable for studying dynamic processes in neural networks and intracellular signaling, offering real-time insights into complex biological mechanisms.</p>
<p>Our Neuroplatform is equipped with all necessary components to perform uncaging. In this example, we investigate closed-loop stimulation, where dopamine is used to reward the network when more spikes are elicited by the same stimulation. The release of the dopamine is achieved through the uncaging of CNV-dopamine using the UV system described in section 3.6.</p>
<p><xref ref-type="fig" rid="fig7">Figure 7B</xref> shows the flow chart of the closed-loop uncaging process. The optimal stimulation parameters are first found using the technique shown in 5.2 (in this case, a current of 4uA, biphasic with 100uS per phase), which is sent successively to each of the 8 electrodes with a delay of 10&#x2009;ms between each electrode.</p>
<p><xref ref-type="fig" rid="fig7">Figure 7C</xref> shows the response timestamps of the 8 electrodes for a period of 1,200&#x2009;ms, 600&#x2009;ms before and after the stimulation. The stimulation event is indicated by the vertical red line. It is interesting to observe that in this particular case, most of the elicited spikes originate from 2 electrodes, specifically electrode 112 and electrode 119.</p>
<p>The Python source code implementing the closed-loop process illustrated in <xref ref-type="fig" rid="fig7">Figure 7B</xref> is provided below. We would like to highlight here how concise the code is. With only 13 lines of code, the entire closed-loop process has been implemented.<inline-graphic xlink:href="frai-07-1376042-i002.tif"/></p>
<p>The graph in <xref ref-type="fig" rid="fig7">Figure 7D</xref> shows the variation in the number of spikes elicited during the execution of the script above across 5&#x2009;h. A general increase in the number of elicited spikes can be observed. However, it is obviously not possible to establish causality between the closed-loop strategy and the observed increase with this single experiment alone. The primary purpose of this closed-loop experiment is to demonstrate the flexibility offered by the Neuroplatform.</p>
</sec>
</sec>
<sec id="sec20">
<label>6</label>
<title>External users of the Neuroplatform</title>
<p>Access to the Neuroplatform is freely available for research purposes. For researchers lacking lab infrastructure, the Neuroplatform provides the capability to conduct real-time experiments on biological networks. Additionally, it allows others to replicate results obtained in their own lab. The database is shared between all research groups, however the Python scripts and Jupyter Notebooks are in private sections.</p>
<p>In 2023, 36 academic groups proposed research projects, of which 8 were selected. At the time of writing, 4 of these have already yielded some results:</p>
<list list-type="bullet">
<list-item>
<p>University C&#x00F4;te d&#x2019;Azure, CNRS, NeuroMod Institute and Laboratoire JA Dieudonn&#x00E9;: investigates the functional connectivity of FO and how electrical stimulation can modify it.</p>
</list-item>
<list-item>
<p>University of Michigan, investigates stimulation protocols that induce global changes in electrical activity of a FO.</p>
</list-item>
<list-item>
<p>Free University of Berlin, investigates stimulation protocols that induce changes in the electrical activity of a FO. Additionally, this research employs machine learning tools to extract information from neural firing patterns and to develop well-conditioned responses. Moreover, it utilizes both shallow and deep reinforcement learning techniques to identify optimal training strategies, aiming to elicit reproducible behaviors in the FO.</p>
</list-item>
<list-item>
<p>University of Exeter, Department of Mathematics and Statistics, Living Systems Institute, investigates storing and retrieving of spatiotemporal spiking patterns, using closed-loop experiments that combine mathematical models of synaptic communication with the Neuroplatform.</p>
</list-item>
<list-item>
<p>Lancaster University Leipzig and University of York: characterizes computational properties of FOs under the reservoir computing model, with a view to building low-power environmental sensors.</p>
</list-item>
<list-item>
<p>Oxford Brookes University, School of Engineering, Computing and Mathematics: investigating the properties of emerging dynamics and criticality within neural organizations using the FOs.</p>
</list-item>
<list-item>
<p>University of Bath, ART-AI, IAH: using the free energy principle and active inference to study the learning capabilities of neurons, embodied in a virtual environment.</p>
</list-item>
<list-item>
<p>University of Bristol: stimulating of FOs based on data gathered from an artificial tactile sensor. Use machine learning techniques to interpret the FO&#x2019;s output, investigating their ability to process real-world data.</p>
</list-item>
</list>
</sec>
<sec id="sec21">
<label>7</label>
<title>Discussion and conclusion</title>
<p>The Neuroplatform has now been operational 24/7 for the past 4&#x2009;years. During this time, the organoids on the MEA have been replaced over 250 times. Considering that we place at least 4 organoids per MEA, and change all the organoids simultaneously, this amounts to testing over 1,000 organoids. Initially, their lifetime was only a few hours, but various improvements, especially related to the microfluidics setup, have extended this to up to 100&#x2009;days in best cases. It is important to note that the spontaneous activity of the organoids can vary over their lifetime, a factor that must be taken into consideration when conducting experiments (<xref ref-type="bibr" rid="ref27">Wagenaar et al., 2006</xref>). Additionally, we observed that the minimum current required to elicit spikes, computed using the method described in section 5.2, is increasing over the lifetime of the organoid. This phenomenon may be linked to an impedance increase caused by glial encapsulation (<xref ref-type="bibr" rid="ref21">Salatino et al., 2017</xref>).</p>
<p>The 24/7 recording strategy as described in section 4.2, results in the constant growth of the database. As of this writing, its size has reached 18 terabytes. This volume encompasses the recording of over 20 billion individual action potentials, each sampled at a 30&#x2009;kHz resolution for 3&#x2009;ms. This extensive dataset is significant not only due to its size but also because it was all recorded in a similar <italic>in-vitro</italic> environment, as described in section 3.2. We are eager to share this data with any interested research group.</p>
</sec>
<sec id="sec22">
<label>8</label>
<title>Future extensions</title>
<p>In the future, we plan to extend the capabilities of our platform to manage a broader range of experimental protocols relevant to wetware computing. For example, we aim to enable a remote control over the injection of specific molecules into the medium, facilitating remote experiments that involve pharmacological manipulation of neuronal activity. This expansion will provide additional degrees of freedom for the automatic optimization of parameters influencing neuroplasticity.</p>
<p>Currently, as detailed in Chapter 2, only one differentiation protocol is used for generating organoids. We plan to introduce additional types of organoid generation protocols soon, with the aim of exploring a broader range of possibilities.</p>
<p>Although 32 research groups requested to access to the Neuroplatform, our current infrastructure only allows us to accommodate 7 groups, considering our own research needs as well. We are in the process of scaling-up the AC/DC hardware system to support more users simultaneously. Additionally, we are currently limited to executing close-loop algorithms for neuroplasticity on one single FO, as these algorithms require sending in real-time adapted simulation signals to each FO. Our software is being updated to run closed-loops in parallel on up to 32 FO.</p>
</sec>
<sec sec-type="methods" id="sec23">
<label>9</label>
<title>Methods</title>
<sec id="sec24">
<label>9.1</label>
<title>Brain organoid generation</title>
<p>Human forebrain organoids were originated as described in <xref ref-type="bibr" rid="ref12">Govindan et al. (2021)</xref>. Briefly, Human Neural Stem Cells derived from the human induced pluripotent stem (hiPS) cell line (ThermoFisher), were plated in flasks coated with CellStart (Fisher Scientific) and amplified in Stempro NSC SFM kit (ThermoFischer) complete medium: KnockOut D-MEM/F12, 2&#x2009;mM of GlutaMAX, 2% of StemPro Neural supplement, 20&#x2009;ng/mL of Human FGF-basic (FGF-2/bFGF) Recombinant Protein, and 20&#x2009;ng/mL of EGF Recombinant Human Protein (Fisher Scientific). Cells were then detached with StemPro<sup>&#x2122;</sup> Accutase (Gibco) and plated in p6 at the concentration of 250,000 cells/well. The plates were sealed with breathable adhesive paper and leads, placed on an orbital shaker at 80&#x2009;rpm, and culture for 7&#x2009;days at 37&#x00B0;C 5% CO2. After one week the newly formed spheroids were put in differentiation medium I (Diff I), containing DMEM/F-12, GlutaMAX<sup>&#x2122;</sup> supplement (Gibco), 2% BSA, 1X of Stempro&#x00AE; hESC Supplement, 20&#x2009;ng/mL of BDNF Recombinant Human Protein (Invitrogen), 20&#x2009;ng/mL of GDNF Recombinant Human Protein (Gibco), 100&#x2009;mM of N6,2&#x2032;-O-Dibutyryladenosine 3&#x2032;,5&#x2032;-cyclic monophosphate sodium salt, and 20&#x2009;mM of 2-Phospho-L-ascorbic acid trisodium salt. After one week, brain spheroids were put in differentiation medium II (Diff II) made of 50% of Diff I and 50% of Neurobasal Plus (Invitrogen). After 3&#x2009;weeks of culture in Diff II, brain organoids were plated in Neurobasal Plus and kept in the orbital shaker until the transfer on the MEA. Medium was change once per week.</p>
</sec>
<sec id="sec25">
<label>9.2</label>
<title>Electron microscopy analysis of FOs</title>
<p>Mature FOs were fixed in 2.5% Glutaraldehyde in 0.1&#x2009;M phosphate buffer pH 7.4, at RT. After 24&#x2009;h the samples were processed as described in <xref ref-type="bibr" rid="ref6">Cakir et al. (2019)</xref> at the Electron Microscopy Facility of University of Lausanne. The whole FO images were acquired with Quanta FEG 250 Scanning Electron Microscope.</p>
</sec>
<sec id="sec26">
<label>9.3</label>
<title>Transfer of FOs on MEA</title>
<p>MEA connected with the microfluid system was moved from the incubator and placed on a 12.3-megapixel camera system (with an optical lens of 16&#x2009;mm of focal, giving a magnification power of 21x) inside the cell culture hood. The lid was removed to access the top of the liquid/air interface. Sterile Hydrophilic PTFE MEMBRANE Hole &#x2018;confetti&#x2019; (diameter 2.5&#x2009;mm, diameter of the hole 0.7&#x2009;mm) (HEPIA) were positioned on top of each electrode and left there 2&#x2009;min to absorb the medium. FOs were collected from the plate using wide bore pipette tips (Axygen) and placed in the middle of confetti, in a 10&#x2009;&#x03BC;L drop of medium. The position of the organoids was adjusted with the help of sterile forceps. After all the organoids were put on place, the chamber was covered with the plate sealer Greiner Bio-One<sup>&#x2122;</sup> BREATHseal<sup>&#x2122;</sup> Sealer (Fisher Scientific), and with the MEA lid. MEA containing the organoids were placed immediately back in the cell incubator and were ready to be used for recording and stimulation. A similar procedure was used for the positioning of organoids on MCS MEA (60MEA200/30iR-Ti). In this case the Hydrophilic PTFE MEMBRANE was not used and organoids were directly laid on the electrodes in a 30&#x2009;&#x03BC;L drop of medium. Recording of organoid activity was performed immediately afterwards.</p>
</sec>
<sec id="sec27">
<label>9.4</label>
<title>System design and assembly</title>
<p>Cell culture media was stored in a 50&#x2009;mL Falcon tube with a multi-port delivery cap (ElveFlow) and stored at 4&#x00B0;C. Each reservoir delivery cap contained a single 0.8&#x2009;mm ID &#x00D7; 1.6&#x2009;mm OD PTFE tubing (Darwin Microfluidics), sealed by a two-piece PFA Fittings and ferrule threaded adapter (IDEX), extending from the bottom of the reservoir to an inlet port on the 4-port valve head of the RVM Rotary Valve (Advance Microfluidics SA). Sterile air is permitted to refill the reservoir through a 0.22-&#x03BC;m filter (Milian) fixed to the cap to compensate for syringe pump medium withdrawal. A similar PTFE tubing and PFA Fittings and adapters were used to connect the syringe pump to the 4-port valve head of the RVM Rotary Valve (Advance Microfluidics SA). Each PTFE tubing coming from the distribution valve connects with a 50&#x2009;mL falcon tube inside the cell culture incubator (Binder) and to a borosilicate glass bottle (Milian) to collect discarded cell culture medium.</p>
<p>A secondary microfluid system made of 0.8&#x2009;mm ID &#x00D7; 1.6&#x2009;mm OD PTFE tubing, were used to connect each 50&#x2009;mL falcon tube inside the cell culture incubator with its own MEA (HEPIA). The connection was through a precise peristaltic pump BT100-2&#x2009;J (Darwin Microfluidics) containing 10 rollers. A compute module (Raspberry Pi 4) controlled the peristaltic pump and the Rotary Valve, through a custom application program interface (API), using RS485 interface and RS-232 interface, respectively. A Fluigent flow-rate sensor connected via USB to the Raspberry Pi 4 allowed the monitoring of the flow rate inside the microfluidic system between the peristaltic pump and the MEA. Python was used to develop the software required to carry out automation protocols.</p>
</sec>
<sec id="sec28">
<label>9.5</label>
<title>Uncaging of dopamine</title>
<p>Carboxynitroveratryl (CNV)-caged dopamine (Tocris Bioscience) was dissolved in Neurobasal Plus at the concentration of 1&#x2009;mM, and injected in the fluidic system. After 3&#x2009;h from the injection, the uncaging experiment started as described in paragraph 5.3. UV Silver-LED fiber-coupled LED (Prizmatix) was used to uncage the dopamine at the wavelength of 365&#x2009;nm for 800&#x2009;ms each time.</p>
</sec>
</sec>
<sec sec-type="data-availability" id="sec29">
<title>Data availability statement</title>
<p>The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.</p>
</sec>
<sec sec-type="ethics-statement" id="sec30">
<title>Ethics statement</title>
<p>Ethical approval was not required for the studies on humans in accordance with the local legislation and institutional requirements because only commercially available established cell lines were used.</p>
</sec>
<sec sec-type="author-contributions" id="sec31">
<title>Author contributions</title>
<p>FJ: Writing &#x2013; original draft, Writing &#x2013; review &#x0026; editing. MK: Writing &#x2013; original draft, Writing &#x2013; review &#x0026; editing. J-MC: Writing &#x2013; original draft, Writing &#x2013; review &#x0026; editing. FB: Writing &#x2013; original draft, Writing &#x2013; review &#x0026; editing. EK: Writing &#x2013; original draft, Writing &#x2013; review &#x0026; editing.</p>
</sec>
</body>
<back>
<sec sec-type="funding-information" id="sec32">
<title>Funding</title>
<p>The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.</p>
</sec>
<ack>
<p>We thank Steve M. Potter and Daniel Burger for their multiple advices and editing, as well as Mathias Reusser for the figures.</p>
</ack>
<sec sec-type="COI-statement" id="sec33">
<title>Conflict of interest</title>
<p>FJ, MK, J-MC, FB, and EK are employed at FinalSpark, Switzerland.</p>
</sec>
<sec id="sec100" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="ref1"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Armer</surname> <given-names>C.</given-names></name> <name><surname>Letronne</surname> <given-names>F.</given-names></name> <name><surname>DeBenedictis</surname> <given-names>E.</given-names></name></person-group> (<year>2023</year>). <article-title>Support academic access to automated cloud labs to improve reproducibility</article-title>. <source>PLoS Biol.</source> <volume>21</volume>:<fpage>e3001919</fpage>. doi: <pub-id pub-id-type="doi">10.1371/journal.pbio.3001919</pub-id>, PMID: <pub-id pub-id-type="pmid">36595506</pub-id></citation></ref>
<ref id="ref2"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bakkum</surname> <given-names>D. J.</given-names></name> <name><surname>Chao</surname> <given-names>Z. C.</given-names></name> <name><surname>Potter</surname> <given-names>S. M.</given-names></name></person-group> (<year>2008</year>). <article-title>Spatio-temporal electrical stimuli shape behavior of an embodied cortical network in a goal-directed learning task</article-title>. <source>J. Neural Eng.</source> <volume>5</volume>, <fpage>310</fpage>&#x2013;<lpage>323</lpage>. doi: <pub-id pub-id-type="doi">10.1088/1741-2560/5/3/004</pub-id>, PMID: <pub-id pub-id-type="pmid">18714127</pub-id></citation></ref>
<ref id="ref3"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Brewer</surname> <given-names>G. J.</given-names></name> <name><surname>Torricelli</surname> <given-names>J. R.</given-names></name></person-group> (<year>2007</year>). <article-title>Isolation and culture of adult neurons and neurospheres</article-title>. <source>Nat. Protoc.</source> <volume>2</volume>, <fpage>1490</fpage>&#x2013;<lpage>1498</lpage>. doi: <pub-id pub-id-type="doi">10.1038/nprot.2007.207</pub-id></citation></ref>
<ref id="ref4"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cai</surname> <given-names>H.</given-names></name> <name><surname>Ao</surname> <given-names>Z.</given-names></name> <name><surname>Tian</surname> <given-names>C.</given-names></name> <name><surname>Wu</surname> <given-names>Z.</given-names></name> <name><surname>Liu</surname> <given-names>H.</given-names></name> <name><surname>Tchieu</surname> <given-names>Z.</given-names></name> <etal/></person-group>. (<year>2023a</year>). <article-title>Brain organoid reservoir computing for artificial intelligence</article-title>. <source>Nat Electron</source> <volume>6</volume>, <fpage>1032</fpage>&#x2013;<lpage>1039</lpage>. doi: <pub-id pub-id-type="doi">10.1038/s41928-023-01069-w</pub-id></citation></ref>
<ref id="ref5"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cai</surname> <given-names>H.</given-names></name> <name><surname>Ao</surname> <given-names>Z.</given-names></name> <name><surname>Tian</surname> <given-names>C.</given-names></name> <name><surname>Wu</surname> <given-names>Z.</given-names></name> <name><surname>Liu</surname> <given-names>H.</given-names></name> <etal/></person-group>. (<year>2023b</year>). <article-title>Brain organoid computing for artificial intelligence</article-title>. <source>bioRxiv</source> <comment>[Preprint]</comment>. doi: <pub-id pub-id-type="doi">10.1101/2023.02.28.530502</pub-id></citation></ref>
<ref id="ref6"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cakir</surname> <given-names>B.</given-names></name> <name><surname>Xiang</surname> <given-names>Y.</given-names></name> <name><surname>Tanaka</surname> <given-names>Y.</given-names></name> <name><surname>Kural</surname> <given-names>M. H.</given-names></name> <name><surname>Parent</surname> <given-names>M.</given-names></name> <name><surname>Kang</surname> <given-names>Y. J.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>Engineering of human brain organoids with a functional vascular-like system</article-title>. <source>Nat. Methods</source> <volume>16</volume>, <fpage>1169</fpage>&#x2013;<lpage>1175</lpage>. doi: <pub-id pub-id-type="doi">10.1038/s41592-019-0586-5</pub-id>, PMID: <pub-id pub-id-type="pmid">31591580</pub-id></citation></ref>
<ref id="ref7"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ciarpella</surname> <given-names>F.</given-names></name> <name><surname>Zamfir</surname> <given-names>R. G.</given-names></name> <name><surname>Campanelli</surname> <given-names>A.</given-names></name> <name><surname>Pedrotti</surname> <given-names>G.</given-names></name> <name><surname>Di Chio</surname> <given-names>M.</given-names></name> <name><surname>Bottani</surname> <given-names>E.</given-names></name> <etal/></person-group>. (<year>2023</year>). <article-title>Generation of mouse hippocampal brain organoids from primary embryonic neural stem cells</article-title>. <source>STAR Protoc.</source> <volume>4</volume>:<fpage>102413</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.xpro.2023.102413</pub-id>, PMID: <pub-id pub-id-type="pmid">37454299</pub-id></citation></ref>
<ref id="ref8"><citation citation-type="book"><person-group person-group-type="author"><name><surname>Clark</surname> <given-names>D. D.</given-names></name> <name><surname>Sokoloff</surname> <given-names>L.</given-names></name></person-group> (<year>1999</year>). &#x201C;<article-title>Circulation and energy metabolism of the brain</article-title>&#x201D; in <source>Basic neurochemistry: Molecular, cellular and medical aspects</source>. <edition>6th</edition> ed (<publisher-loc>New York</publisher-loc>: <publisher-name>Philadelphia Lippincott-Raven Publishers</publisher-name>), <fpage>637</fpage>&#x2013;<lpage>669</lpage>.</citation></ref>
<ref id="ref9"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>De Vries</surname> <given-names>A.</given-names></name></person-group> (<year>2023</year>). <article-title>The growing energy footprint of artificial intelligence</article-title>. <source>Joule</source> <volume>7</volume>, <fpage>2191</fpage>&#x2013;<lpage>2194</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.joule.2023.09.004</pub-id></citation></ref>
<ref id="ref10"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Elliott</surname> <given-names>M. A. T.</given-names></name> <name><surname>Schweiger</surname> <given-names>H. E.</given-names></name> <name><surname>Robbins</surname> <given-names>A.</given-names></name> <name><surname>Vera-Choqqueccota</surname> <given-names>S.</given-names></name> <name><surname>Ehrlich</surname> <given-names>D.</given-names></name> <name><surname>Hernandez</surname> <given-names>S.</given-names></name> <etal/></person-group>. (<year>2023</year>). <article-title>Internet-connected cortical organoids for project-based stem cell and neuroscience education</article-title>. <source>eNeuro</source> <volume>10</volume>. doi: <pub-id pub-id-type="doi">10.1523/ENEURO.0308-23.2023</pub-id>, PMID: <pub-id pub-id-type="pmid">38016807</pub-id></citation></ref>
<ref id="ref11"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gienger</surname> <given-names>M.</given-names></name> <name><surname>H&#x00FC;bner</surname> <given-names>H.</given-names></name> <name><surname>L&#x00F6;ber</surname> <given-names>S.</given-names></name> <name><surname>K&#x00F6;nig</surname> <given-names>B.</given-names></name> <name><surname>Gmeiner</surname> <given-names>P.</given-names></name></person-group> (<year>2020</year>). <article-title>Structure-based development of caged dopamine D2/D3 receptor antagonists</article-title>. <source>Sci. Rep.</source> <volume>10</volume>:<fpage>829</fpage>. doi: <pub-id pub-id-type="doi">10.1038/s41598-020-57770-9</pub-id>, PMID: <pub-id pub-id-type="pmid">31965029</pub-id></citation></ref>
<ref id="ref12"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Govindan</surname> <given-names>S.</given-names></name> <name><surname>Batti</surname> <given-names>L.</given-names></name> <name><surname>Osterop</surname> <given-names>S. F.</given-names></name> <name><surname>Stoppini</surname> <given-names>L.</given-names></name> <name><surname>Roux</surname> <given-names>A.</given-names></name></person-group> (<year>2021</year>). <article-title>Mass generation, neuron labeling, and 3D imaging of Minibrains</article-title>. <source>Front. Bioeng. Biotechnol.</source> <volume>8</volume>:<fpage>582650</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fbioe.2020.582650</pub-id>, PMID: <pub-id pub-id-type="pmid">33598450</pub-id></citation></ref>
<ref id="ref13"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gross</surname> <given-names>G. W.</given-names></name> <name><surname>Rieske</surname> <given-names>E.</given-names></name> <name><surname>Kreutzberg</surname> <given-names>G. W.</given-names></name> <name><surname>Meyer</surname> <given-names>A.</given-names></name></person-group> (<year>1977</year>). <article-title>A new fixed-array multimicroelectrode system designed for long-term monitoring of extracellular single unit neuronal activity in vitro</article-title>. <source>Neurosci. Lett.</source> <volume>6</volume>, <fpage>101</fpage>&#x2013;<lpage>105</lpage>. doi: <pub-id pub-id-type="doi">10.1016/0304-3940(77)90003-9</pub-id></citation></ref>
<ref id="ref14"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kagan</surname> <given-names>B. J.</given-names></name> <name><surname>Kitchen</surname> <given-names>A. C.</given-names></name> <name><surname>Tran</surname> <given-names>N. T.</given-names></name> <name><surname>Habibollahi</surname> <given-names>F.</given-names></name> <name><surname>Khajehnejad</surname> <given-names>M.</given-names></name> <name><surname>Parker</surname> <given-names>B. J.</given-names></name> <etal/></person-group>. (<year>2022</year>). <article-title>In vitro neurons learn and exhibit sentience when embodied in a simulated game-world</article-title>. <source>Neuron</source> <volume>110</volume>, <fpage>3952</fpage>&#x2013;<lpage>3969.e8</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.neuron.2022.09.001</pub-id>, PMID: <pub-id pub-id-type="pmid">36228614</pub-id></citation></ref>
<ref id="ref15"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lee</surname> <given-names>S. E.</given-names></name> <name><surname>Shin</surname> <given-names>N.</given-names></name> <name><surname>Kook</surname> <given-names>G.</given-names></name> <name><surname>Kong</surname> <given-names>D.</given-names></name> <name><surname>Kim</surname> <given-names>N. G.</given-names></name> <name><surname>Choi</surname> <given-names>S. W.</given-names></name> <etal/></person-group>. (<year>2020</year>). <article-title>Human iNSC-derived brain organoid model of lysosomal storage disorder in Niemann&#x2013;pick disease type C</article-title>. <source>Cell Death Dis.</source> <volume>11</volume>:<fpage>1059</fpage>. doi: <pub-id pub-id-type="doi">10.1038/s41419-020-03262-7</pub-id></citation></ref>
<ref id="ref16"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Newman</surname> <given-names>J. P.</given-names></name> <name><surname>Zeller-Townson</surname> <given-names>R.</given-names></name> <name><surname>Fong</surname> <given-names>M.</given-names></name> <name><surname>Arcot Desai</surname> <given-names>S.</given-names></name> <name><surname>Gross</surname> <given-names>R. E.</given-names></name> <name><surname>Potter</surname> <given-names>S. M.</given-names></name></person-group> (<year>2013</year>). <article-title>Closed-loop, multichannel experimentation using the open-source NeuroRighter electrophysiology platform</article-title>. <source>Front. Neur. Circ.</source> <volume>6</volume>:<fpage>98</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fncir.2012.00098</pub-id>, PMID: <pub-id pub-id-type="pmid">23346047</pub-id></citation></ref>
<ref id="ref17"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>O&#x2019;Leary</surname> <given-names>G.</given-names></name> <name><surname>Khramtsov</surname> <given-names>I.</given-names></name> <name><surname>Ramesh</surname> <given-names>R.</given-names></name> <name><surname>Perez-Ignacio</surname> <given-names>A.</given-names></name> <name><surname>Shah</surname> <given-names>P.</given-names></name> <name><surname>Chameh</surname> <given-names>H. M.</given-names></name> <etal/></person-group>. (<year>2022</year>). <article-title>OpenMEA: open-source microelectrode Array platform for bioelectronic interfacing</article-title>. <source>BioRxiv</source>. doi: <pub-id pub-id-type="doi">10.1101/2022.11.11.516234</pub-id></citation></ref>
<ref id="ref18"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pa&#x015F;ca</surname> <given-names>S. P.</given-names></name> <name><surname>Arlotta</surname> <given-names>P.</given-names></name> <name><surname>Bateup</surname> <given-names>H. S.</given-names></name> <name><surname>Camp</surname> <given-names>G.</given-names></name> <name><surname>Cappello</surname> <given-names>S.</given-names></name> <name><surname>Gage</surname> <given-names>F. H.</given-names></name> <etal/></person-group>. (<year>2022</year>). <article-title>A nomenclature consensus for nervous system organoids and assembloids</article-title>. <source>Nature</source> <volume>609</volume>, <fpage>907</fpage>&#x2013;<lpage>910</lpage>. doi: <pub-id pub-id-type="doi">10.1038/s41586-022-05219-6</pub-id>, PMID: <pub-id pub-id-type="pmid">36171373</pub-id></citation></ref>
<ref id="ref19"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pine</surname> <given-names>J.</given-names></name></person-group> (<year>1980</year>). <article-title>Recording action potentials from cultured neurons with extracellular microcircuit electrodes</article-title>. <source>J Neurosci. Meth</source> <volume>2</volume>, <fpage>19</fpage>&#x2013;<lpage>31</lpage>. doi: <pub-id pub-id-type="doi">10.1016/0165-0270(80)90042-4</pub-id>, PMID: <pub-id pub-id-type="pmid">7329089</pub-id></citation></ref>
<ref id="ref20"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Qian</surname> <given-names>X.</given-names></name> <name><surname>Song</surname> <given-names>H.</given-names></name> <name><surname>Ming</surname> <given-names>G. L.</given-names></name></person-group> (<year>2019</year>). <article-title>Brain organoids: advances, applications and challenges</article-title>. <source>Development</source> <volume>146</volume>:<fpage>dev166074</fpage>. doi: <pub-id pub-id-type="doi">10.1242/dev.166074</pub-id>, PMID: <pub-id pub-id-type="pmid">30992274</pub-id></citation></ref>
<ref id="ref21"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Salatino</surname> <given-names>J. W.</given-names></name> <name><surname>Ludwig</surname> <given-names>K. A.</given-names></name> <name><surname>Kozai</surname> <given-names>T. D. Y.</given-names></name> <name><surname>Purcell</surname> <given-names>E. K.</given-names></name></person-group> (<year>2017</year>). <article-title>Glial responses to implanted electrodes in the brain</article-title>. <source>Nat. Biomed. Eng.</source> <volume>1</volume>, <fpage>862</fpage>&#x2013;<lpage>877</lpage>. doi: <pub-id pub-id-type="doi">10.1038/s41551-017-0154-1</pub-id>, PMID: <pub-id pub-id-type="pmid">30505625</pub-id></citation></ref>
<ref id="ref22"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Samsi</surname> <given-names>S.</given-names></name> <name><surname>Zhao</surname> <given-names>D.</given-names></name> <name><surname>McDonald</surname> <given-names>J.</given-names></name> <name><surname>Li</surname> <given-names>B.</given-names></name> <name><surname>Michalea</surname> <given-names>A.</given-names></name> <name><surname>Jones</surname> <given-names>M.</given-names></name> <etal/></person-group>. (<year>2023</year>). <article-title>From words to Watts: benchmarking the energy costs of large language model inference</article-title>. <source>arXiv</source>:2310.03003 [cs.CL] doi: <pub-id pub-id-type="doi">10.48550/arXiv.2310.03003</pub-id></citation></ref>
<ref id="ref23"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Smirnova</surname> <given-names>L.</given-names></name> <name><surname>Caffo</surname> <given-names>B. S.</given-names></name> <name><surname>Gracias</surname> <given-names>D. H.</given-names></name> <name><surname>Huang</surname> <given-names>Q.</given-names></name> <name><surname>Morales Pantoja</surname> <given-names>I. E.</given-names></name> <name><surname>Tang</surname> <given-names>B.</given-names></name> <etal/></person-group>. (<year>2023</year>). <article-title>Organoid intelligence (OI): the new frontier in biocomputing and intelligence-in-a-dish</article-title>. <source>Front. Sci.</source> <volume>1</volume>:<fpage>1017235</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fsci.2023.1017235</pub-id></citation></ref>
<ref id="ref24"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Stoppini</surname> <given-names>L.</given-names></name> <name><surname>Buchs</surname> <given-names>P.-A.</given-names></name> <name><surname>Muller</surname> <given-names>D.</given-names></name></person-group> (<year>1991</year>). <article-title>A simple method for organotypic cultures of nervous tissue</article-title>. <source>J. Neurosc. Methods</source> <volume>37</volume>, <fpage>173</fpage>&#x2013;<lpage>182</lpage>. doi: <pub-id pub-id-type="doi">10.1016/0165-0270(91)90128-M</pub-id></citation></ref>
<ref id="ref25"><citation citation-type="other"><person-group person-group-type="author"><name><surname>Wagenaar</surname> <given-names>D. A.</given-names></name> <name><surname>DeMarse</surname> <given-names>T. B.</given-names></name> <name><surname>Potter</surname> <given-names>S. M.</given-names></name></person-group>, "MeaBench: a toolset for multi-electrode data acquisition and on-line analysis.," In 2nd Intl. IEEE EMBS conference on neural Engigneering, Arlington, VA, USA, (<year>2005a</year>).</citation></ref>
<ref id="ref26"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wagenaar</surname> <given-names>D. A.</given-names></name> <name><surname>Madhavan</surname> <given-names>R.</given-names></name> <name><surname>Pine</surname> <given-names>J.</given-names></name> <name><surname>Potter</surname> <given-names>S. M.</given-names></name></person-group> (<year>2005b</year>). <article-title>Controlling bursting in cortical cultures with closed-loop multi-electrode stimulation</article-title>. <source>J. Neuroscience</source> <volume>25</volume></citation></ref>
<ref id="ref27"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wagenaar</surname> <given-names>D. A.</given-names></name> <name><surname>Pine</surname> <given-names>J.</given-names></name> <name><surname>Potter</surname> <given-names>S. M.</given-names></name></person-group> (<year>2006</year>). <article-title>An extremely rich repertoire of bursting patterns during the development of cortical cultures</article-title>. <source>BMC Neurosci.</source> <volume>7</volume>:<fpage>11</fpage>. doi: <pub-id pub-id-type="doi">10.1186/1471-2202-7-11</pub-id></citation></ref>
<ref id="ref28"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wertenbroek</surname> <given-names>R.</given-names></name> <name><surname>Thoma</surname> <given-names>Y.</given-names></name> <name><surname>Mor</surname> <given-names>F. M.</given-names></name> <name><surname>Grassi</surname> <given-names>S.</given-names></name> <name><surname>Heuschkel</surname> <given-names>M. O.</given-names></name> <name><surname>Roux</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2021</year>). <article-title>SpikeOnChip: a custom embedded platform for neuronal activity recording and analysis</article-title>. <source>IEEE Trans. Biomed. Circuits Syst.</source> <volume>15</volume>, <fpage>743</fpage>&#x2013;<lpage>755</lpage>. doi: <pub-id pub-id-type="doi">10.1109/TBCAS.2021.3097833</pub-id>, PMID: <pub-id pub-id-type="pmid">34280107</pub-id></citation></ref>
<ref id="ref30"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yang</surname> <given-names>X.</given-names></name> <name><surname>Forr&#x00F3;</surname> <given-names>C.</given-names></name> <name><surname>Li</surname> <given-names>L. T.</given-names></name> <name><surname>Miura</surname> <given-names>L.</given-names></name> <name><surname>Zaluska</surname> <given-names>T. J.</given-names></name> <name><surname>Tsai</surname> <given-names>C. T.</given-names></name> <etal/></person-group>. (<year>2024</year>). <article-title>Kirigami electronics for long-term electrophysiological recording of human neural organoids and assembloids</article-title>. <source>Nat. Biotechnol.</source> doi: <pub-id pub-id-type="doi">10.1038/s41587-023-02081-3</pub-id>, PMID: <pub-id pub-id-type="pmid">38253880</pub-id></citation></ref>
<ref id="ref31"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>X.</given-names></name> <name><surname>Dou</surname> <given-names>Z.</given-names></name> <name><surname>Kim</surname> <given-names>S. H.</given-names></name> <name><surname>Upadhyay</surname> <given-names>G.</given-names></name> <name><surname>Havert</surname> <given-names>D.</given-names></name> <name><surname>Kang</surname> <given-names>S.</given-names></name> <etal/></person-group>. (<year>2023</year>). <article-title>Mind in vitro platforms: versatile, scalable, robust, and open solutions to interfacing with living neurons</article-title>. <source>Adv. Sci. (Weinh)</source>:<fpage>e2306826</fpage>. doi: <pub-id pub-id-type="doi">10.1002/advs.202306826</pub-id>. <comment>Online ahead of print</comment>.</citation></ref>
</ref-list>
</back>
</article>