People who use drugs (PWUD) are at a high risk of contracting and developing severe coronavirus disease 2019 (COVID-19) and other infectious diseases due to their lifestyle, comorbidities, and the detrimental effects of opioids on cellular immunity. However, there is limited research on vaccine responses in PWUD, particularly regarding the role that T cells play in the immune response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Here, we show that before vaccination, PWUD did not exhibit an increased frequency of preexisting cross-reactive T cells to SARS-CoV-2 and that, despite the inhibitory effects that opioids have on T-cell immunity, standard vaccination can elicit robust polyfunctional CD4+ and CD8+ T-cell responses that were similar to those found in controls. Our findings indicate that vaccination stimulates an effective immune response in PWUD and highlight targeted vaccination as an essential public health instrument for the control of COVID-19 and other infectious diseases in this group of high-risk patients.

Poor overall survival of hematopoietic stem cell transplantation (HSCT) recipients who developed COVID-19 underlies the importance of SARS-CoV-2 vaccination. Previous studies of vaccine efficacy have reported weak humoral responses but conflicting results on T cell immunity. Here, we have examined the relationship between humoral and T cell response in 48 HSCT recipients who received two doses of Moderna’s mRNA-1273 or Pfizer/BioNTech’s BNT162b2 vaccines. Nearly all HSCT patients had robust T cell immunity regardless of protective humoral responses, with 18/48 (37%, IQR 8.679-5601 BAU/mL) displaying protective IgG anti-receptor binding domain (RBD) levels (>2000 BAU/mL). Flow cytometry analysis of activation induced markers (AIMs) revealed that 90% and 74% of HSCT patients showed reactivity towards immunodominant spike peptides in CD8 + and CD4 + T cells, respectively. The response rate increased to 90% for CD4 + T cells as well when we challenged the cells with a complete set of overlapping peptides spanning the entire spike protein. T cell response was detectable as early as 3 months after transplant, but only CD4 + T cell reactivity correlated with IgG anti-RBD level and time after transplantation. Boosting increased seroconversion rate, while only one patient developed COVID-19 requiring hospitalization. Our data suggest that HSCT recipients with poor serological responses were protected from severe COVID-19 by vaccine-induced T cell responses.

During the COVID-19 pandemic we utilized an AI-driven T cell epitope prediction tool, the NEC Immune Profiler (NIP) to scrutinize and predict regions of T cell immunogenicity (hotspots) from the entire SARS-CoV-2 viral proteome. These immunogenic regions offer potential for the development of universally protective T cell vaccine candidates. Here, we validated and characterized T cell responses to a set of minimal epitopes from these AI-identified universal hotspots. Utilizing a flow cytometry-based T cell activation-induced marker (AIM) assay, we identified 59 validated screening hits, of which 56% (33 peptides) have not been previously reported. Notably, we found that most of these novel epitopes were derived from the non-spike regions of SARS-CoV-2 (Orf1ab, Orf3a, and E). In addition, ex vivo stimulation with NIP-predicted peptides from the spike protein elicited CD8+ T cell response in PBMC isolated from most vaccinated donors. Our data confirm the predictive accuracy of AI platforms modelling bona fide immunogenicity and provide a novel framework for the evaluation of vaccine-induced T cell responses.

During the COVID-19 pandemic, several SARS-CoV-2 variants of concern (VOC) emerged, bringing with them varying degrees of health and socioeconomic burdens. In particular, the Omicron VOC displayed distinct features of increased transmissibility accompanied by antigenic drift in the spike protein that partially circumvented the ability of pre-existing antibody responses in the global population to neutralize the virus. However, T cell immunity has remained robust throughout all the different VOC transmission waves and has emerged as a critically important correlate of protection against SARS-CoV-2 and its VOCs, in both vaccinated and infected individuals. Therefore, as SARS-CoV-2 VOCs continue to evolve, it is crucial that we characterize the correlates of protection and the potential for immune escape for both B cell and T cell human immunity in the population. Generating the insights necessary to understand T cell immunity, experimentally, for the global human population is at present a critical but a time consuming, expensive, and laborious process. Further, it is not feasible to generate global or universal insights into T cell immunity in an actionable time frame for potential future emerging VOCs. However, using computational means we can expedite and provide early insights into the correlates of T cell protection. In this study, we generated and revealed insights on the T cell epitope landscape for the five main SARS-CoV-2 VOCs observed to date. We demonstrated using a unique AI prediction platform, a significant conservation of presentable T cell epitopes across all mutated peptides for each VOC. This was modeled using the most frequent HLA alleles in the human population and covers the most common HLA haplotypes in the human population. The AI resource generated through this computational study and associated insights may guide the development of T cell vaccines and diagnostics that are even more robust against current and future VOCs, and their emerging subvariants.

The global population is at present suffering from a pandemic of Coronavirus disease 2019 (COVID-19), caused by the novel coronavirus Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). The goal of this study was to use artificial intelligence (AI) to predict blueprints for designing universal vaccines against SARS-CoV-2, that contain a sufficiently broad repertoire of T-cell epitopes capable of providing coverage and protection across the global population. To help achieve these aims, we profiled the entire SARS-CoV-2 proteome across the most frequent 100 HLA-A, HLA-B and HLA-DR alleles in the human population, using host-infected cell surface antigen presentation and immunogenicity predictors from the NEC Immune Profiler suite of tools, and generated comprehensive epitope maps. We then used these epitope maps as input for a Monte Carlo simulation designed to identify statistically significant "epitope hotspot" regions in the virus that are most likely to be immunogenic across a broad spectrum of HLA types. We then removed epitope hotspots that shared significant homology with proteins in the human proteome to reduce the chance of inducing off-target autoimmune responses. We also analyzed the antigen presentation and immunogenic landscape of all the nonsynonymous mutations across 3,400 different sequences of the virus, to identify a trend whereby SARS-COV-2 mutations are predicted to have reduced potential to be presented by host-infected cells, and consequently detected by the host immune system. A sequence conservation analysis then removed epitope hotspots that occurred in less-conserved regions of the viral proteome. Finally, we used a database of the HLA haplotypes of approximately 22,000 individuals to develop a "digital twin" type simulation to model how effective different combinations of hotspots would work in a diverse human population; the approach identified an optimal constellation of epitope hotspots that could provide maximum coverage in the global population. By combining the antigen presentation to the infected-host cell surface and immunogenicity predictions of the NEC Immune Profiler with a robust Monte Carlo and digital twin simulation, we have profiled the entire SARS-CoV-2 proteome and identified a subset of epitope hotspots that could be harnessed in a vaccine formulation to provide a broad coverage across the global population.

This protocol predicts blueprints for vaccine design that contain a broad repertoire of T-cell epitopes optimized for the global population. The protocol first requires a screening of the SARS-CoV-2 proteome using immunogenicity predictors to generate comprehensive epitope maps. Then, these epitope maps are used as input to Monte Carlo simulations designed to identify statistically significant “epitope hotspot” regions in the virus that are most likely to be immunogenic. The epitope hotspots that share significant homology with proteins in the human proteome are removed to reduce the chance of inducing off-target autoimmune responses. Finally, a database of the actual HLA genotypes of citizens is used to develop a “digital twin” type simulation to model how effective different combinations of hotspots would work in a diverse human population. The approach identifies an optimal constellation of epitope hotspots that could provide maximum coverage in the human population.

Pituitary endocrine cells fire action potentials (APs) to regulate their cytosolic Ca2+ concentration and hormone secretion rate. Depending on animal species, cell type, and biological conditions, pituitary APs are generated either by TTX-sensitive Na+ currents (INa), high-voltage activated Ca2+ currents (ICa), or by a combination of the two. Previous computational models of pituitary cells have mainly been based on data from rats, where INa is largely inactivated at the resting potential, and spontaneous APs are predominantly mediated by ICa. Unlike in rats, spontaneous INa-mediated APs are consistently seen in pituitary cells of several other animal species, including several species of fish. In the current work we develop a computational model of gonadotropin releasing cells in the teleost fish medaka (Oryzias latipes). The model stands out from previous modeling efforts by being (1) the first model of a pituitary cell in teleosts, (2) the first pituitary cell model that fires sponateous APs that are predominantly mediated by INa, and (3) the first pituitary cell model where the kinetics of the depolarizing currents, INa and ICa, are directly fitted to voltage-clamp data. We explore the firing properties of the model, and compare it to the properties of previous models that fire ICa-based APs. We put a particular focus on how the big conductance K+ current (IBK) modulates the AP shape. Interestingly, we find that IBK can prolong AP duration in models that fire ICa-based APs, while it consistently shortens the duration of the predominantly INa-mediated APs in the medaka gonadotroph model. Although the model is constrained to experimental data from gonadotroph cells in medaka, it may likely provide insights also into other pituitary cell types that fire INa-mediated APs.

Natural sciences generate an increasing amount of data in a wide range of formats developed by different research groups and commercial companies. At the same time there is a growing desire to share data along with publications in order to enable reproducible research. Open formats have publicly available specifications which facilitate data sharing and reproducible research. Hierarchical Data Format 5 (HDF5) is a popular open format widely used in neuroscience, often as a foundation for other, more specialized formats. However, drawbacks related to HDF5's complex specification have initiated a discussion for an improved replacement. We propose a novel alternative, the Experimental Directory Structure (Exdir), an open specification for data storage in experimental pipelines which amends drawbacks associated with HDF5 while retaining its advantages. HDF5 stores data and metadata in a hierarchy within a complex binary file which, among other things, is not human-readable, not optimal for version control systems, and lacks support for easy access to raw data from external applications. Exdir, on the other hand, uses file system directories to represent the hierarchy, with metadata stored in human-readable YAML files, datasets stored in binary NumPy files, and raw data stored directly in subdirectories. Furthermore, storing data in multiple files makes it easier to track for version control systems. Exdir is not a file format in itself, but a specification for organizing files in a directory structure. Exdir uses the same abstractions as HDF5 and is compatible with the HDF5 Abstract Data Model. Several research groups are already using data stored in a directory hierarchy as an alternative to HDF5, but no common standard exists. This complicates and limits the opportunity for data sharing and development of common tools for reading, writing, and analyzing data. Exdir facilitates improved data storage, data sharing, reproducible research, and novel insight from interdisciplinary collaboration. With the publication of Exdir, we invite the scientific community to join the development to create an open specification that will serve as many needs as possible and as a foundation for open access to and exchange of data.

Computational models in neuroscience typically contain many parameters that are poorly constrained by experimental data. Uncertainty quantification and sensitivity analysis provide rigorous procedures to quantify how the model output depends on this parameter uncertainty. Unfortunately, the application of such methods is not yet standard within the field of neuroscience. Here we present Uncertainpy, an open-source Python toolbox, tailored to perform uncertainty quantification and sensitivity analysis of neuroscience models. Uncertainpy aims to make it quick and easy to get started with uncertainty analysis, without any need for detailed prior knowledge. The toolbox allows uncertainty quantification and sensitivity analysis to be performed on already existing models without needing to modify the model equations or model implementation. Uncertainpy bases its analysis on polynomial chaos expansions, which are more efficient than the more standard Monte-Carlo based approaches. Uncertainpy is tailored for neuroscience applications by its built-in capability for calculating characteristic features in the model output. The toolbox does not merely perform a point-to-point comparison of the "raw" model output (e.g., membrane voltage traces), but can also calculate the uncertainty and sensitivity of salient model response features such as spike timing, action potential width, average interspike interval, and other features relevant for various neural and neural network models. Uncertainpy comes with several common models and features built in, and including custom models and new features is easy. The aim of the current paper is to present Uncertainpy to the neuroscience community in a user-oriented manner. To demonstrate its broad applicability, we perform an uncertainty quantification and sensitivity analysis of three case studies relevant for neuroscience: the original Hodgkin-Huxley point-neuron model for action potential generation, a multi-compartmental model of a thalamic interneuron implemented in the NEURON simulator, and a sparsely connected recurrent network model implemented in the NEST simulator.

AIM:Cachexia is a severe wasting disorder involving loss of body- and muscle mass reducing survival and quality of life in cancer patients. We aim at determining if cachexia is a mere perturbation of the protein balance or if the condition also involves a degenerative loss of myonuclei within the fibre syncytia or loss of whole muscle fibres.METHODS:We induced cachexia by xenografting PC3 prostate cancer cells in nu/nu mice. Six weeks later, we counted myonuclei by in vivo microscopic imaging of single live fibres in the extensor digitorum longus muscle (EDL), and the EDL, soleus and tibialis anterior muscles were also harvested for ex vivo histology.RESULTS:The mice lost on average 15% of the whole-body wt. The muscle wet weight of the glycolytic, fast EDL was reduced by 14%, the tibialis anterior by 17%, and the slow, oxidative soleus by 6%. The fibre cross-sectional area in the EDL was reduced by 21% with no loss of myonuclei or any significant reduction in the number of muscle fibres. TUNEL-positive nuclei or fibres with embryonic myosin were rare both in cachectic and control muscles, and haematoxylin-eosin staining revealed no clear signs of muscle pathology.CONCLUSION:The data suggest that the cachexia induced by xenografted prostate tumours induces a pronounced atrophy not accompanied by a loss of myonuclei or a loss of muscle fibres. Thus, stem cell related treatment might be redundant, and the quest for treatment options should rather focus on intervening with intracellular pathways regulating muscle fibre size.

Pituitary endocrine cells fire action potentials (APs) to regulate their cytosolic Ca2+ concentration and hormone secretion rate. Depending on animal species, cell type, and biological conditions, pituitary APs are generated either by TTX-sensitive Na+ currents ( I Na ), high-voltage activated Ca2+ currents ( I Ca ), or by a combination of the two. Previous computational models of pituitary cells have mainly been based on data from rats, where I Na is largely inactivated at the resting potential, and spontaneous APs are exclusively mediated by I Ca . As a part of the previous modeling studies, a paradoxical role was identified for the big conductance K+ current ( I BK ), which was found to prolong the duration of I Ca -mediated APs, and sometimes give rise to pseudo-plateau bursts, contrary to what one would expect from a hyperpolarizing current. Unlike in rats, spontaneous I Na -mediated APs are consistently seen in pituitary cells of several other animal species, including several species of fish. In the current work we develop the, to our knowledge, first computational model of a pituitary cell that fires I Na -mediated APs. Although we constrain the model to experimental data from gonadotrope cells in the teleost fish medaka ( Oryzias latipes ), it may likely provide insights also into other pituitary cell types that fire I Na -mediated APs. In the current work, we use the model to explore how the effect of I BK depends on the AP generating mechanisms of pituitary cells. We do this by comparing simulations on the medaka gonadotrope model (two versions thereof) with simulations on a previously developed model of a rat pituitary cell. Interestingly, we find that I BK has the opposite effect on APs in the two models, i.e. it reduces the duration of already fast I Na -mediated APs in the medaka model, and prolongs the duration of already slow I Ca -mediated APs in the rat model. Author summary Excitable cells elicit electrical pulses called action potentials (APs), which are generated and shaped by a combination of ion channels in the cell membrane. While neurons use APs for interneuronal communication and heart cells use them to generate heart-beats, pituitary cells use APs to regulate their cytosolic Ca2+ concentration, which in turn controls their hormone secretion rate. The amount of Ca2+ that enters the pituitary cell during an AP depends strongly on how long it lasts, and it is therefore important to understand the mechanisms that control this. Depending on animal species and biological conditions, pituitary APs may be initiated either by Ca2+ channels or Na+ channels. Here, we explore the differences between the two scenarios by comparing simulations on two different computer models: (i) a previously developed model which fires Na+-based APs, adapted to data from pituitary cells in rats, and (ii) a novel model that fires Ca2+-based APs, adapted to data from pituitary cells in the fish medaka. Interestingly, we find that the role of big conductance K+ (BK) channels, which are known to affect the duration of the AP, are opposite in the two models, i.e., they act to prolong Ca2+-based APs while they act to shorten Na+-based APs.

Educational software (apps) can improve science education by providing an interactive way of learning about complicated topics that are hard to explain with text and static illustrations. However, few educational apps are available for simulation of neural networks. Here, we describe an educational app, Neuronify, allowing the user to easily create and explore neural networks in a plug-and-play simulation environment. The user can pick network elements with adjustable parameters from a menu, i.e., synaptically connected neurons modelled as integrate-and-fire neurons and various stimulators (current sources, spike generators, visual, and touch) and recording devices (voltmeter, spike detector, and loudspeaker). We aim to provide a low entry point to simulation-based neuroscience by allowing students with no programming experience to create and simulate neural networks. To facilitate the use of Neuronify in teaching, a set of premade common network motifs is provided, performing functions such as input summation, gain control by inhibition, and detection of direction of stimulus movement. Neuronify is developed in C++ and QML using the cross-platform application framework Qt and runs on smart phones (Android, iOS) and tablet computers as well personal computers (Windows, Mac, Linux).

Perineuronal nets (PNNs) are specialized extracellular matrix (ECM) structures that condense around the soma and proximal dendrites of subpopulations of neurons. Emerging evidence suggests that they are involved in regulating brain plasticity. However, the expression of PNNs varies between and within brain areas. A lack of quantitative studies describing the distribution and cell-specificity of PNNs makes it difficult to reveal the functional roles of PNNs. In the current study, we examine the distribution of PNNs and the identity of PNN-enwrapped neurons in three brain areas with different cognitive functions: the dorsal hippocampus, medial entorhinal cortex (mEC) and primary visual cortex (V1). We compared rats and mice as knowledge from these species are often intermingled. The most abundant expression of PNNs was found in the mEC and V1, while dorsal hippocampus showed strikingly low levels of PNNs, apart from dense expression in the CA2 region. In hippocampus we also found apparent species differences in expression of PNNs. While we confirm that the PNNs enwrap parvalbumin-expressing (PV+) neurons in V1, we found that they mainly colocalize with excitatory CamKII-expressing neurons in CA2. In mEC, we demonstrate that in addition to PV+ cells, the PNNs colocalize with reelin-expressing stellate cells. We also show that the maturation of PNNs in mEC coincides with the formation of grid cell pattern, while PV+ cells, unlike in other cortical areas, are present from early postnatal development. Finally, we demonstrate considerable effects on the number of PSD-95-gephyrin puncta after enzymatic removal of PNNs.