Publications

Environmental perturbations during the first years of life are a major factor in psychiatric diseases. Phencyclidine (PCP), a drug of abuse, has psychomimetic effects, and neonatal subchronic administration of PCP in rodents leads to long-term behavioral changes relevant for schizophrenia. The cerebellum is increasingly recognized for its role in diverse cognitive functions. However, little is known about potential cerebellar changes in models of schizophrenia. Here, we analyzed the characteristics of the cerebellum in the neonatal subchronic PCP model. We found that, while the global cerebellar cytoarchitecture and Purkinje cell spontaneous spiking properties are unchanged, climbing fiber/Purkinje cell synaptic connectivity is increased in juvenile mice. Neonatal subchronic administration of PCP is accompanied by increased cFos expression, a marker of neuronal activity, and transient modification of the neuronal surfaceome in the cerebellum. The largest change observed is the overexpression of Ctgf, a gene previously suggested as a biomarker for schizophrenia. This neonatal increase in Ctgf can be reproduced by increasing neuronal activity in the cerebellum during the second postnatal week using chemogenetics. However, it does not lead to increased climbing fiber/Purkinje cell connectivity in juvenile mice, showing the complexity of PCP action. Overall, our study shows that administration of the drug of abuse PCP during the developmental period of intense cerebellar synaptogenesis and circuit remodeling has long-term and specific effects on Purkinje cell connectivity and warrants the search for this type of synaptic changes in psychiatric diseases.

Computing the spike-triggered average (STA) is a simple method to estimate linear receptive fields (RFs) in sensory neurons. For random, uncorrelated stimuli, the STA provides an unbiased RF estimate, but in practice, white noise at high resolution is not an optimal stimulus choice as it usually evokes only weak responses. Therefore, for a visual stimulus, images of randomly modulated blocks of pixels are often used. This solution naturally limits the resolution at which an RF can be measured. Here, we present a simple super-resolution technique that can overcome these limitations. We define a novel stimulus type, the shifted white noise (SWN), by introducing random spatial shifts in the usual stimulus to increase the resolution of the measurements. In simulated data, we show that the average error using the SWN was 1.7 times smaller than when using the classical stimulus, with successful mapping of 2.3 times more neurons, covering a broader range of RF sizes. Moreover, successful RF mapping was achieved with brief recordings of light responses, lasting only about 1 min of activity, which is more than 10 times more efficient than the classical white noise stimulus. In recordings from mouse retinal ganglion cells with large scale multielectrode arrays, we successfully mapped 21 times more RFs than when using the traditional white noise stimuli. In summary, randomly shifting the usual white noise stimulus significantly improves RFs estimation, and requires only short recordings.

In the study of biological networks, one of the major challenges is to understand the relationships between network structure and dynamics. In this paper, we model in vitro cortical neuronal cultures as stochastic dynamical systems and apply a method that reconstructs directed networks from dynamics [Ching and Tam, Phys. Rev. E 95, 010301(R) (2017)] to reveal directed effective connectivity, namely, the directed links and synaptic weights, of the neuronal cultures from voltage measurements recorded by a multielectrode array. The effective connectivity so obtained reproduces several features of cortical regions in rats and monkeys and has similar network properties as the synaptic network of the nematode Caenorhabditis elegans, whose entire nervous system has been mapped out. The distribution of the incoming degree is bimodal and the distributions of the average incoming and outgoing synaptic strength are non-Gaussian with long tails. The effective connectivity captures different information from the commonly studied functional connectivity, estimated using statistical correlation between spiking activities. The average synaptic strengths of excitatory incoming and outgoing links are found to increase with the spiking activity in the estimated effective connectivity but not in the functional connectivity estimated using the same sets of voltage measurements. These results thus demonstrate that the reconstructed effective connectivity can capture the general properties of synaptic connections and better reveal relationships between network structure and dynamics.

Neuropathic pain such as that seen in diabetes mellitus, results in part from central sensitisation in the dorsal horn. However, the mechanisms responsible for such sensitisation remain unclear. There is evidence that disturbances in the integrity of the spinal vascular network can be causative factors in the development of neuropathic pain. Here we show that reduced blood flow and vascularity of the dorsal horn leads to the onset of neuropathic pain. Using rodent models (type 1 diabetes and an inducible endothelial specific vascular endothelial growth factor receptor 2 knockout mouse) that result in degeneration of the endothelium in the dorsal horn we show that spinal cord vasculopathy results in nociceptive behavioural hypersensitivity. This also results in increased hypoxia in dorsal horn neurons, depicted by increased expression of hypoxia markers hypoxia inducible factor 1𝛼, glucose transporter 3 and carbonic anhydrase 7. Furthermore, inducing hypoxia via intrathecal delivery of dimethyloxalylglycine leads to the activation of dorsal horn neurons as well as mechanical and thermal hypersensitivity. This shows that hypoxic signalling induced by reduced vascularity results in increased hypersensitivity and pain. Inhibition of carbonic anhydrase activity, through intraperitoneal injection of acetazolamide, inhibited hypoxia induced pain behaviours. This investigation demonstrates that induction of a hypoxic microenvironment in the dorsal horn, as occurs in diabetes, is an integral process by which neurons are activated to initiate neuropathic pain states. This leads to the conjecture that reversing hypoxia by improving spinal cord microvascular blood flow could reverse or prevent neuropathic pain.

The cerebellar cortex microcircuit is characterized by a highly ordered neuronal architecture having a relatively simple and stereotyped connectivity pattern. For a long time, this structural simplicity has incorrectly led to the idea that anatomical considerations would be sufficient to understand the dynamics of the underlying circuitry. However, recent experimental evidence indicates that cerebellar operations are much more complex than solely predicted by anatomy, due to the crucial role played by neuronal and synaptic properties. To be able to explore neuronal and microcircuit dynamics, advanced imaging, electrophysiological techniques and computational models have been combined, allowing us to investigate neuronal ensembles activity and to connect microscale to mesoscale phenomena. Here, we review what is known about cerebellar network organization, neural dynamics and synaptic plasticity and point out what is still missing and would require experimental assessments. We consider the available experimental techniques that allow a comprehensive assessment of circuit dynamics, including voltage and calcium imaging and extracellular electrophysiological recordings with multi-electrode arrays (MEAs). These techniques are proving essential to investigate the spatiotemporal pattern of activity and plasticity in the cerebellar network, providing new clues on how circuit dynamics contribute to motor control and higher cognitive functions.

Retinal neurons are remarkedly diverse based on structure, function and genetic identity. Classifying these cells is a challenging task, requiring multimodal methodology. Here, we introduce a novel approach for retinal ganglion cell (RGC) classification, based on pharmacogenetics combined with immunohistochemistry and large-scale retinal electrophysiology. Our novel strategy allows grouping of cells sharing gene expression and understanding how these cell classes respond to basic and complex visual scenes. Our approach consists of several consecutive steps. First, the spike firing frequency is increased in RGCs co-expressing a certain gene (Scnn1a or Grik4) using excitatory DREADDs (designer receptors exclusively activated by designer drugs) in order to single out activity originating specifically from these cells. Their spike location is then combined with post hoc immunostaining, to unequivocally characterize their anatomical and functional features. We grouped these isolated RGCs into multiple clusters based on spike train similarities. Using this novel approach, we were able to extend the pre-existing list of Grik4-expressing RGC types to a total of eight and, for the first time, we provide a phenotypical description of 13 Scnn1a-expressing RGCs. The insights and methods gained here can guide not only RGC classification but neuronal classification challenges in other brain regions as well.

‍

How does the dynamics of neurons in a network respond to changes in synaptic weights? Answer to this question would be important for a full understanding of synaptic plasticity. In this article, we report our numerical study of the effects of changes in inhibitory synaptic weights on the spontaneous activity of networks of spiking neurons with conductance-based synapses. Networks with biologically realistic features, which were reconstructed from multi-electrode array recordings taken in a cortical neuronal culture, and their modifications were used in the simulations. The magnitudes of the synaptic weights of all the inhibitory connections are decreased by a uniform amount subjecting to the condition that inhibitory connections would not be turned into excitatory ones. Our simulation results reveal that the responses of the neurons are heterogeneous: while the firing rate of some neurons increases as expected, the firing rate of other neurons decreases or remains unchanged. The same results show that heterogeneous responses also occur for an enhancement of inhibition. This heterogeneity in the responses of neurons to changes in inhibitory synaptic strength suggests that activity-induced modification of synaptic strength does not necessarily generate a positive feedback loop on the dynamics of neurons connected in a network. Our results could be used to understand the effects of bicuculline on spiking and bursting activities of neuronal cultures. Using reconstructed networks with biologically realistic features enables us to identify a long-tailed distribution of average synaptic weights for outgoing links as a crucial feature in giving rise to bursting in neuronal networks and in determining the overall response of the whole network to changes in synaptic strength. For networks whose average synaptic weights for outgoing links have a long-tailed distribution, bursting is observed and the average firing rate of the whole network increases upon inhibition suppression or decreases upon inhibition enhancement. For networks whose average synaptic weights for outgoing links are approximately normally distributed, bursting is not found and the average firing rate of the whole network remains approximately constant upon changes in inhibitory synaptic strength.

Retinal drug toxicity screening is essential for the development of safe treatment strategies for a large number of diseases. To this end, retinal organoids derived from human pluripotent stem cells (hPSCs) provide a suitable screening platform due to their similarity to the human retina and the ease of generation in large-scale formats. In this study, two hPSC cell lines were differentiated to retinal organoids, which comprised all key retinal cell types in multiple nuclear and synaptic layers. Single-cell RNA-Seq of retinal organoids indicated the maintenance of retinal ganglion cells and development of bipolar cells: both cell types segregated into several subtypes. Ketorolac, digoxin, thioridazine, sildenafil, ethanol, and methanol were selected as key compounds to screen on retinal organoids because of their well-known retinal toxicity profile described in the literature. Exposure of the hPSC-derived retinal organoids to digoxin, thioridazine, and sildenafil resulted in photoreceptor cell death, while digoxin and thioridazine additionally affected all other cell types, including Müller glia cells. All drug treatments caused activation of astrocytes, indicated by dendrites sprouting into neuroepithelium. The ability to respond to light was preserved in organoids although the number of responsive retinal ganglion cells decreased after drug exposure. These data indicate similar drug effects in organoids to those reported in in vivo models and/or in humans, thus providing the first robust experimental evidence of their suitability for toxicological studies.

Large-scale multi-site biosensors are essential to probe the olfactory bulb (OB) circuitry for understanding the spatiotemporal dynamics of simultaneous discharge patterns. Current ex-vivo biosensing techniques are limited to recording a small set of neurons and cannot provide an adequate resolution, which hinders revealing the fast dynamic underlying the information coding mechanisms in the OB circuit. Here, we demonstrate a novel biohybrid OB-CMOS biosensing platform to decipher the cross-scale dynamics of the OB electrogenesis and quantify the distinct neuronal coding properties. The approach with 4096-microelectrodes offers a non-invasive, label-free, bioelectrical imaging to decode simultaneous firing patterns from thousands of connected neuronal ensembles in acute OB slices. The platform can measure spontaneous and drug-induced extracellular field potential activity with substantially improved spatiotemporal resolution over conventional OB-based biosensors. Also, we employ our OB-CMOS recordings to perform multidimensional analysis to instantiate specific neurophysiological metrics underlying the olfactory spatiotemporal coding that emerged from the OB interconnected layers. Our results delineate the computational implications of large-scale activity patterns in functional olfactory processing. The systematic interplay of the experimental CMOS-base platform architecture and the high-content characterization of the olfactory circuit with various computational analyses endow significant functional interrogations of the OB information processing, high-spatiotemporal connectivity mapping, and global circuit dynamics. Thus, our study can inspire the design of advanced biomimetic olfactory-based biosensors and neuromorphic approaches for diagnostic biomarkers and drug discovery applications.

Neurovascular coupling (NVC) is the process associating local cerebral blood flow (CBF) to neuronal activity (NA). Although NVC provides the basis for the blood oxygen level dependent (BOLD) effect used in functional MRI (fMRI), the relationship between NVC and NA is still unclear. Since recent studies reported cerebellar non-linearities in BOLD signals during motor tasks execution, we investigated the NVC/NA relationship using a range of input frequencies in acute mouse cerebellar slices of vermis and hemisphere. The capillary diameter increased in response to mossy fiber activation in the 6–300 Hz range, with a marked inflection around 50 Hz (vermis) and 100 Hz (hemisphere). The corresponding NA was recorded using high-density multi-electrode arrays and correlated to capillary dynamics through a computational model dissecting the main components of granular layer activity. Here, NVC is known to involve a balance between the NMDAR-NO pathway driving vasodilation and the mGluRs-20HETE pathway driving vasoconstriction. Simulations showed that the NMDAR-mediated component of NA was sufficient to explain the time course of the capillary dilation but not its non-linear frequency dependence, suggesting that the mGluRs-20HETE pathway plays a role at intermediate frequencies. These parallel control pathways imply a vasodilation–vasoconstriction competition hypothesis that could adapt local hemodynamics at the microscale bearing implications for fMRI signals interpretation.

‍

Highlights

• FeP, a fluorescent probe, is suitable for brain imaging ferrous iron flux in vivo

• 3D two-photon imaging reveals elevation of ferrous iron in the epileptic mouse brain

• Dihydroartemisinin (DHA) can modulate iron homeostasis in the epileptic mouse

• DHA exhibits a potential antiepileptic effect in the epileptic mouse model

Summary

Imbalanced iron homeostasis plays a crucial role in neurological diseases, yet direct imaging evidence revealing the distribution of active ferrous iron (Fe²⁺) in the living brain remains scarce. Here, we present a near-infrared excited two-photon fluorescent probe (FeP) for imaging changes of Fe²⁺ flux in the living epileptic mouse brain. In vivo 3D two-photon brain imaging with FeP directly revealed abnormal elevation of Fe²⁺ in the epileptic mouse brain. Moreover, we found that dihydroartemisinin (DHA), a lead compound discovered through probe-based high-throughput screening, plays a critical role in modulating iron homeostasis. In addition, we revealed that DHA might exert its antiepileptic effects by modulating iron homeostasis in the brain and finally inhibiting ferroptosis. This work provides a reliable chemical tool for assessing the status of ferrous iron in the living epileptic mouse brain and may aid the rapid discovery of antiepileptic drug candidates.

Imbalanced iron homeostasis plays a crucial role in neurological diseases, yet direct imaging evidence revealing the distribution of active ferrous iron (Fe²⁺) in the living brain remains scarce. Here, we present a near-infrared excited two-photon fluorescent probe (FeP) for imaging changes of Fe²⁺ flux in the living epileptic mouse brain. In vivo 3D two-photon brain imaging with FeP directly revealed abnormal elevation of Fe²⁺ in the epileptic mouse brain. Moreover, we found that dihydroartemisinin (DHA), a lead compound discovered through probe-based high-throughput screening, plays a critical role in modulating iron homeostasis. In addition, we revealed that DHA might exert its antiepileptic effects by modulating iron homeostasis in the brain and finally inhibiting ferroptosis. This work provides a reliable chemical tool for assessing the status of ferrous iron in the living epileptic mouse brain and may aid the rapid discovery of antiepileptic drug candidates.

Optical stimulation technologies are gaining great consideration in cardiology, neuroscience studies, and drug discovery pathways by providing control over cell activity with high spatio-temporal resolution. However, this high precision requires manipulation of biological processes at genetic level concealing its development from broad scale application. Therefore, translating these technologies into tools for medical or pharmacological applications remains a challenge. Here, an all-optical nongenetic method for the modulation of electrogenic cells is introduced. It is demonstrated that plasmonic metamaterials can be used to elicit action potentials by converting near infrared laser pulses into stimulatory currents. The suggested approach allows for the stimulation of cardiomyocytes and neurons directly on commercial complementary metal-oxide semiconductor microelectrode arrays coupled with ultrafast pulsed laser, providing both stimulation and network-level recordings on the same device.

Rho GTPases are a class of G-proteins involved in several aspects of cellular biology, including the regulation of actin cytoskeleton. The most studied members of this family are RHOA and RAC1 that act in concert to regulate actin dynamics. Recently, Rho GTPases gained much attention as synaptic regulators in the mammalian central nervous system (CNS). In this context, ARHGAP22 protein has been previously shown to specifically inhibit RAC1 activity thus standing as critical cytoskeleton regulator in cancer cell models; however, whether this function is maintained in neurons in the CNS is unknown. Here, we generated a knockout animal model for arhgap22 and provided evidence of its role in the hippocampus. Specifically, we found that ARHGAP22 absence leads to RAC1 hyperactivity and to an increase in dendritic spine density with defects in synaptic structure, molecular composition, and plasticity. Furthermore, arhgap22 silencing causes impairment in cognition and a reduction in anxiety-like behavior in mice. We also found that inhibiting RAC1 restored synaptic plasticity in ARHGAP22 KO mice. All together, these results shed light on the specific role of ARHGAP22 in hippocampal excitatory synapse formation and function as well as in learning and memory behaviors.

Human-induced pluripotent stem cells (hiPSCs) with their differentiation protocols, constitute a potential tool to investigate the various biological mechanisms of different human cells, such as those of the central nervous system. With the advent of such technique, we have increased the knowledge of biological mechanisms of neural diseases and, new therapies are now emerging. In particular, three-dimensional (3D) neural cell culture models including brain organoids and neurospheroids are increasingly used as in vitro platforms for studying human brain cell biology and drug screening, in genome engineering and transplantation as potential treatment for some neurodegenerative diseases. In this work, we exploited a particular differentiation protocol to generate engineered excitatory cortical neurospheroids of human origin. To assess functional network activity, we used standard Micro Electrodes Arrays (60 channels) and for the first time CMOS based devices (4096 channels). Sample cultures showed electrophysiological activity in 4 weeks and these first results suggest future possible applications for drug screening and transplantation.

Methods for the conversion of human induced pluripotent stem cells (hiPSCs) into motor neurons (MNs) have opened to the generation of patient-derived in vitro systems that can be exploited for MN disease modelling. However, the lack of simplified and consistent protocols and the fact that hiPSC-derived MNs are often functionally immature yet limit the opportunity to fully take advantage of this technology, especially in research aimed at revealing the disease phenotypes that are manifested in functionally mature cells. In this study, we present a robust, optimized monolayer procedure to rapidly convert hiPSCs into enriched populations of motor neuron progenitor cells (MNPCs) that can be further amplified to produce a large number of cells to cover many experimental needs. These MNPCs can be efficiently differentiated towards mature MNs exhibiting functional electrical and pharmacological neuronal properties. Finally, we report that MN cultures can be long-term maintained, thus offering the opportunity to study degenerative phenomena associated with pathologies involving MNs and their functional, networked activity. These results indicate that our optimized procedure enables the efficient and robust generation of large quantities of MNPCs and functional MNs, providing a valid tool for MNs disease modelling and for drug discovery applications.

Defects in DNA single‐strand break repair (SSBR) are linked with neurological dysfunction but the underlying mechanisms remain poorly understood. Here, we show that hyperactivity of the DNA strand break sensor protein Parp1 in mice in which the central SSBR protein Xrcc1 is conditionally deleted (Xrcc1^Nes‐Cre) results in lethal seizures and shortened lifespan. Using electrophysiological recording and synaptic imaging approaches, we demonstrate that aberrant Parp1 activation triggers seizure‐like activity in Xrcc1‐defective hippocampus ex vivo and deregulated presynaptic calcium signalling in isolated hippocampal neurons in vitro. Moreover, we show that these defects are prevented by Parp1 inhibition or deletion and, in the case of Parp1 deletion, that the lifespan of Xrcc1^Nes‐Cre mice is greatly extended. This is the first demonstration that lethal seizures can be triggered by aberrant Parp1 activity at unrepaired SSBs, highlighting PARP inhibition as a possible therapeutic approach in hereditary neurological disease.

‍

Drug-induced cardiotoxicity is a major barrier to drug development and a main cause of withdrawal of marketed drugs. Drugs can strongly alter the spontaneous functioning of the heart by interacting with the cardiac membrane ion channels. If these effects only surface during in vivo preclinical tests, clinical trials or worse after commercialization, the societal and economic burden will be significant and seriously hinder the efficient drug development process. Hence, cardiac safety pharmacology requires in vitro electrophysiological screening assays of all drug candidates to predict cardiotoxic effects before clinical trials. In the past 10 years, microelectrode array (MEA) technology began to be considered a valuable approach in pharmaceutical applications. However, an effective tool for high-throughput intracellular measurements, compatible with pharmaceutical standards, is not yet available. Here, we propose laser-induced optoacoustic poration combined with CMOS-MEA technology as a reliable and effective platform to detect cardiotoxicity. This approach enables the acquisition of high-quality action potential recordings from large numbers of cardiomyocytes within the same culture well, providing reliable data using single-well MEA devices and single cardiac syncytia per each drug. Thus, this technology could be applied in drug safety screening platforms reducing times and costs of cardiotoxicity assessments, while simultaneously improving the data reliability.

‍

‍

Fine control of protein stoichiometry at synapses underlies brain function and plasticity. How proteostasis is controlled independently for each type of synaptic protein in a synapse-specific and activity-dependent manner remains unclear. Here, we show that Susd4, a gene coding for a complement-related transmembrane protein, is expressed by many neuronal populations starting at the time of synapse formation. Constitutive loss-of-function of Susd4 in the mouse impairs motor coordination adaptation and learning, prevents long-term depression at cerebellar synapses, and leads to misregulation of activity-dependent AMPA receptor subunit GluA2 degradation. We identified several proteins with known roles in the regulation of AMPA receptor turnover, in particular ubiquitin ligases of the NEDD4 subfamily, as SUSD4 binding partners. Our findings shed light on the potential role of SUSD4 mutations in neurodevelopmental diseases.

Retinal dystrophies often lead to blindness. Developing therapeutic interventions to restore vision is therefore of paramount importance. Here we demonstrate the ability of pluripotent stem cell-derived cone precursors to engraft and restore light responses in the Pde6brd1 mouse, an end-stage photoreceptor degeneration model. Our data show that up to 1.5% of precursors integrate into the host retina, differentiate into cones, and engraft in close apposition to the host bipolar cells. Half of the transplanted mice exhibited visual behavior and of these 33% showed binocular light sensitivity. The majority of retinal ganglion cells exhibited contrast-sensitive ON, OFF or ON-OFF light responses and even motion sensitivity; however, quite a few exhibited unusual responses (eg, light-induced suppression), presumably reflecting remodeling of the neural retina. Our data indicate that despite relatively low engraftment yield, pluripotent stem cell-derived cone precursors can elicit light responsiveness even at advanced degeneration stages. Further work is needed to improve engraftment yield and counteract retinal remodeling to achieve useful clinical applications.

The retina is a complex circuit of the central nervous system whose aim is to encode visual stimuli prior the higher order processing performed in the visual cortex. Due to the importance of its role, modeling the retina to advance in interpreting its spiking activity output is a well studied problem. In particular, it has been shown that latent variable models can be used to model the joint distribution of Retinal Ganglion Cells (RGCs). In this work, we validate the applicability of Restricted Boltzmann Machines to model the spiking activity responses of a large a population of RGCs recorded with high-resolution electrode arrays. In particular, we show that latent variables can encode modes in the RGC activity distribution that are closely related to the visual stimuli. In contrast to previous work, we further validate our findings by comparing results associated with recordings from retinas under normal and altered encoding conditions obtained by pharmacological manipulation. In these conditions, we observe that the model reflects well-known physiological behaviors of the retina. Finally, we show that we can also discover temporal patterns, associated with distinct dynamics of the stimuli.

It is widely assumed that electrical synapses in the mammalian brain, especially between interneurons, underlie neuronal synchrony. In the hippocampus, principal cells also establish electrical synapses with each other and have also been implicated in network oscillations, whereby the origin of fast electrical activity has been attributed to ectopic spikelets and dendro-dendritic or axo-axonal gap junctions. However, if electrical synapses were in axo-dendritic connections, where chemical synapses occur, the synaptic events would be mixed, having an electrical component preceding the chemical one. This type of communication is less well studied, mainly because it is not easily detected. Moreover, a possible scenario could be that an electrical synapse coexisted with a chemical one, but in a nonconductive state; hence, it would be considered inexistent. Could chemical synapses have a quiescent electrical component? If so, can silent electrical synapses be activated to be detected? We addressed this possibility, and we here report that, indeed, the connexin-36-containing glutamatergic mossy fiber synapses of the rat hippocampus express previously unrecognized electrical synapses, which are normally silent. We reveal that these synapses are pH sensitive, actuate in vitro and in vivo, and that the electrical signaling is bidirectional. With the simultaneous recording of hundreds of cells, we could reveal the existence of an electrical circuit in the hippocampus of adult rats of either sex consisting of principal cells where the nodes are interregional glutamatergic synapses containing silent but ready-to-use gap junctions.

The generation of laminated and light responsive retinal organoids from induced pluripotent stem cells (iPSCs) provides a powerful tool for the study of retinal diseases and drug discovery and a robust platform for cell-based therapies. The aim of this study is to investigate whether retinal organoids can retain their morphological and functional characteristics upon storage at room temperature (RT) conditions and shipment by air using a commercially available container that maintains the environment at ambient temperature. Morphological analysis and measurements of neuroepithelial thickness revealed no differences between control, RT incubated and shipped organoids. Similarly immunohistochemical analysis showed no differences in cell type composition and position within the laminated retinal structure. All groups showed a similar response to light, suggesting that the biological function of retinal organoids was not affected by RT storage or shipment. These findings provide an advance in transport of ready-made retinal organoids, increasing their availability to many research and pharma labs worldwide and facilitating cross-collaborative research.

The mammalian claustrum, owing to its widespread connectivity with other forebrain structures, has been hypothesized to mediate functions that range from decision-making to consciousness¹. Here we report that a homologue of the claustrum, identified by single-cell transcriptomics and viral tracing of connectivity, also exists in a reptile—the Australian bearded dragon Pogona vitticeps. In Pogona, the claustrum underlies the generation of sharp waves during slow-wave sleep. The sharp waves, together with superimposed high-frequency ripples², propagate to the entire neighbouring pallial dorsal ventricular ridge (DVR). Unilateral or bilateral lesions of the claustrum suppress the production of sharp-wave ripples during slow-wave sleep in a unilateral or bilateral manner, respectively, but do not affect the regular and rapidly alternating sleep rhythm that is characteristic of sleep in this species³. The claustrum is thus not involved in the generation of the sleep rhythm itself. Tract tracing revealed that the reptilian claustrum projects widely to a variety of forebrain areas, including the cortex, and that it receives converging inputs from, among others, areas of the mid- and hindbrain that are known to be involved in wake–sleep control in mammals^4,5,6. Periodically modulating the concentration of serotonin in the claustrum, for example, caused a matching modulation of sharp-wave production there and in the neighbouring DVR. Using transcriptomic approaches, we also identified a claustrum in the turtle Trachemys scripta, a distant reptilian relative of lizards. The claustrum is therefore an ancient structure that was probably already present in the brain of the common vertebrate ancestor of reptiles and mammals. It may have an important role in the control of brain states owing to the ascending input it receives from the mid- and hindbrain, its widespread projections to the forebrain and its role in sharp-wave generation during slow-wave sleep.

Electrical coupling has been reported to occur only between homotypic retinal ganglion cells, in line with the concept of parallel processing in the early visual system. Here, however, we show reciprocal correlated firing between heterotypic ganglion cells in multielectrode array recordings during light stimulation in retinas of adult guinea pigs of either sex. Heterotypic coupling was further confirmed via tracer spread after intracellular injections of single cells with neurobiotin. Both electrically coupled cell types were sustained ON center ganglion cells but showed distinct light response properties and receptive field sizes. We identified one of the involved cell types as sustained ON α-ganglion cells. The presence of electrical coupling between heterotypic ganglion cells introduces a network motif in which the signals of distinct ganglion cell types are partially mixed at the output stage of the retina.

Large-scale neural recording methods now allow us to observe large populations of identified single neurons simultaneously, opening a window into neural population dynamics in living organisms. However, distilling such large-scale recordings to build theories of emergent collective dynamics remains a fundamental statistical challenge. The neural field models of Wilson, Cowan, and colleagues remain the mainstay of mathematical population modeling owing to their interpretable, mechanistic parameters and amenability to mathematical analysis. Inspired by recent advances in biochemical modeling, we develop a method based on moment closure to interpret neural field models as latent state-space point-process models, making them amenable to statistical inference. With this approach we can infer the intrinsic states of neurons, such as active and refractory, solely from spiking activity in large populations. After validating this approach with synthetic data, we apply it to high-density recordings of spiking activity in the developing mouse retina. This confirms the essential role of a long lasting refractory state in shaping spatiotemporal properties of neonatal retinal waves. This conceptual and methodological advance opens up new theoretical connections between mathematical theory and point-process state-space models in neural data analysis.