The tumour microenvironment shapes innate lymphoid cells in patients with hepatocellular carcinoma.

OBJECTIVE: Hepatocellular carcinoma (HCC) represents a typical inflammation-associated cancer. Tissue resident innate lymphoid cells (ILCs) have been suggested to control tumour surveillance. Here, we studied how the local cytokine milieu controls ILCs in HCC. DESIGN: We performed bulk RNA sequencing of HCC tissue as well as flow cytometry and single-cell RNA sequencing of enriched ILCs from non-tumour liver, margin and tumour core derived from 48 patients with HCC. Simultaneous measurement of protein and RNA expression at the single-cell level (AbSeq) identified precise signatures of ILC subgroups. In vitro culturing of ILCs was used to validate findings from in silico analysis. Analysis of RNA-sequencing data from large HCC cohorts allowed stratification and survival analysis based on transcriptomic signatures. RESULTS: RNA sequencing of tumour, non-tumour and margin identified tumour-dependent gradients, which were associated with poor survival and control of ILC plasticity. Single-cell RNA sequencing and flow cytometry of ILCs from HCC livers identified natural killer (NK)-like cells in the non-tumour tissue, losing their cytotoxic profile as they transitioned into tumour ILC1 and NK-like-ILC3 cells. Tumour ILC composition was mediated by cytokine gradients that directed ILC plasticity towards activated tumour ILC2s. This was liver-specific and not seen in ILCs from peripheral blood mononuclear cells. Patients with high ILC2/ILC1 ratio expressed interleukin-33 in the tumour that promoted ILC2 generation, which was associated with better survival. CONCLUSION: Our results suggest that the tumour cytokine milieu controls ILC composition and HCC outcome. Specific changes of cytokines modify ILC composition in the tumour by inducing plasticity and alter ILC function.

The effects of telomere shortening on cancer cells: a network model of proteomic and microRNA analysis.

Previously, we have shown that shortening of telomeres by telomerase inhibition sensitized cancer cells to cisplatinum, slowed their migration, increased DNA damage and impaired DNA repair. The mechanism behind these effects is not fully characterized. Its clarification could facilitate novel therapeutics development and may obviate the time consuming process of telomere shortening achieved by telomerase inhibition. Here we aimed to decipher the microRNA and proteomic profiling of cancer cells with shortened telomeres and identify the key mediators in telomere shortening-induced damage to those cells. Of 870 identified proteins, 98 were differentially expressed in shortened-telomere cells. 47 microRNAs were differentially expressed in these cells; some are implicated in growth arrest or act as oncogene repressors. The obtained data was used for a network construction, which provided us with nodal candidates that may mediate the shortened-telomere dependent features. These proteins' expression was experimentally validated, supporting their potential central role in this system.

Common and specific signatures of gene expression and protein-protein interactions in autoimmune diseases.

The aim of this study is to understand intracellular regulatory mechanisms in peripheral blood mononuclear cells (PBMCs), which are either common to many autoimmune diseases or specific to some of them. We incorporated large-scale data such as protein-protein interactions, gene expression and demographical information of hundreds of patients and healthy subjects, related to six autoimmune diseases with available large-scale gene expression measurements: multiple sclerosis (MS), systemic lupus erythematosus (SLE), juvenile rheumatoid arthritis (JRA), Crohn's disease (CD), ulcerative colitis (UC) and type 1 diabetes (T1D). These data were analyzed concurrently by statistical and systems biology approaches tailored for this purpose. We found that chemokines such as CXCL1-3, 5, 6 and the interleukin (IL) IL8 tend to be differentially expressed in PBMCs of patients with the analyzed autoimmune diseases. In addition, the anti-apoptotic gene BCL3, interferon-γ (IFNG), and the vitamin D receptor (VDR) gene physically interact with significantly many genes that tend to be differentially expressed in PBMCs of patients with the analyzed autoimmune diseases. In general, similar cellular processes tend to be differentially expressed in PBMC in the analyzed autoimmune diseases. Specifically, the cellular processes related to cell proliferation (for example, epidermal growth factor, platelet-derived growth factor, nuclear factor-κB, Wnt/ß-catenin signaling, stress-activated protein kinase c-Jun NH2-terminal kinase), inflammatory response (for example, interleukins IL2 and IL6, the cytokine granulocyte-macrophage colony-stimulating factor and the B-cell receptor), general signaling cascades (for example, mitogen-activated protein kinase, extracellular signal-regulated kinase, p38 and TRK) and apoptosis are activated in most of the analyzed autoimmune diseases. However, our results suggest that in each of the analyzed diseases, apoptosis and chemotaxis are activated via different subsignaling pathways. Analyses of the expression levels of dozens of genes and the protein-protein interactions among them demonstrated that CD and UC have relatively similar gene expression signatures, whereas the gene expression signatures of T1D and JRA relatively differ from the signatures of the other autoimmune diseases. These diseases are the only ones activated via the Fcɛ pathway. The relevant genes and pathways reported in this study are discussed at length, and may be helpful in the diagnoses and understanding of autoimmunity and/or specific autoimmune diseases.

A systems level strategy for analyzing the cell death network: implication in exploring the apoptosis/autophagy connection.

The mammalian cell death network comprises three distinct functional modules: apoptosis, autophagy and programmed necrosis. Currently, the field lacks systems level approaches to assess the extent to which the intermodular connectivity affects cell death performance. Here, we developed a platform that is based on single and double sets of RNAi-mediated perturbations targeting combinations of apoptotic and autophagic genes. The outcome of perturbations is measured both at the level of the overall cell death responses, using an unbiased quantitative reporter, and by assessing the molecular responses within the different functional modules. Epistatic analyses determine whether seemingly unrelated pairs of proteins are genetically linked. The initial running of this platform in etoposide-treated cells, using a few single and double perturbations, identified several levels of connectivity between apoptosis and autophagy. The knock down of caspase3 turned on a switch toward autophagic cell death, which requires Atg5 or Beclin-1. In addition, a reciprocal connection between these two autophagic genes and apoptosis was identified. By applying computational tools that are based on mining the protein-protein interaction database, a novel biochemical pathway connecting between Atg5 and caspase3 is suggested. Scaling up this platform into hundreds of perturbations potentially has a wide, general scope of applicability, and will provide the basis for future modeling of the cell death network.

The effect of phenotypic plasticity on evolution in multipeaked fitness landscapes.

When facing the challenge of developing an individual that best fits its environment, nature demonstrates an interesting combination of two fundamentally different adaptive mechanisms: genetic evolution and phenotypic plasticity. Following numerous computational models, it has become the accepted wisdom that lifetime acclimation (e.g. via learning) smooths the fitness landscape and consequently accelerates evolution. However, analytical studies, focusing on the effect of phenotypic plasticity on evolution in simple unimodal landscapes, have often found that learning hinders the evolutionary process rather than accelerating it. Here, we provide a general framework for studying the effect of plasticity on evolution in multipeaked landscapes and introduce a rigorous mathematical analysis of these dynamics. We show that the convergence rate of the evolutionary process in a given arbitrary one-dimensional fitness landscape is dominated by the largest descent (drawdown) in the landscape and provide numerical evidence to support an analogous dominance also in multidimensional landscapes. We consider several schemes of phenotypic plasticity and examine their effect on the landscape drawdown, identifying the conditions under which phenotypic plasticity is advantageous. The lack of such a drawdown in unimodal landscapes vs. its dominance in multipeaked landscapes accounts for the seemingly contradictory findings of previous studies.

Distributed synchrony in a cell assembly of spiking neurons.

We investigate the formation of a Hebbian cell assembly of spiking neurons, using a temporal synaptic learning curve that is based on recent experimental findings. It includes potentiation for short time delays between pre- and post-synaptic neuronal spiking, and depression for spiking events occurring in the reverse order. The coupling between the dynamics of synaptic learning and that of neuronal activation leads to interesting results. One possible mode of activity is distributed synchrony, implying spontaneous division of the Hebbian cell assembly into groups, or subassemblies, of cells that fire in a cyclic manner. The behavior of distributed synchrony is investigated both by simulations and by analytic calculations of the resulting synaptic distributions.

Cortical spreading depression and the pathogenesis of brain disorders: a computational and neural network-based investigation.

This paper reviews our recent studies of the role of cortical spreading depression (CSD) in the pathogenesis of brain disorders. Our investigation is a computational one, involving the development and utilization of a complex neuro-metabolic model of the interactions assumed to occur in the cortex during the passage of multiple CSD waves. Incorporating these neuro-metabolic changes of CSD within a neural network model of normoxic cortex produces cortical activation patterns during the passage of a CSD wave that, projected onto the visual fields, resemble the visual hallucinations observed during the migraine aura. When focal ischemia is simulated with the model, the evoked CSD waves are found to affect the expansion of the infarction into the ischemic penumbra. Our findings support the hypothesis that CSD does play an important pathogenic role in these and other neurological disorders, and suggest additional experimental studies that may further substantiate it.

Effective neuronal learning with ineffective Hebbian learning rules.

In this article we revisit the classical neuroscience paradigm of Hebbian learning. We find that it is difficult to achieve effective associative memory storage by Hebbian synaptic learning, since it requires network-level information at the synaptic level or sparse coding level. Effective learning can yet be achieved even with nonsparse patterns by a neuronal process that maintains a zero sum of the incoming synaptic efficacies. This weight correction improves the memory capacity of associative networks from an essentially bounded one to a memory capacity that scales linearly with network size. It also enables the effective storage of patterns with multiple levels of activity within a single network. Such neuronal weight correction can be successfully carried out by activity-dependent homeostasis of the neuron's synaptic efficacies, which was recently observed in cortical tissue. Thus, our findings suggest that associative learning by Hebbian synaptic learning should be accompanied by continuous remodeling of neuronally driven regulatory processes in the brain.

Emergence of memory-driven command neurons in evolved artificial agents.

Using evolutionary simulations, we develop autonomous agents controlled by artificial neural networks (ANNs). In simple lifelike tasks of foraging and navigation, high performance levels are attained by agents equipped with fully recurrent ANN controllers. In a set of experiments sharing the same behavioral task but differing in the sensory input available to the agents, we find a common structure of a command neuron switching the dynamics of the network between radically different behavioral modes. When sensory position information is available, the command neuron reflects a map of the environment, acting as a location-dependent cell sensitive to the location and orientation of the agent. When such information is unavailable, the command neuron's activity is based on a spontaneously evolving short-term memory mechanism, which underlies its apparent place-sensitive activity. A two-parameter stochastic model for this memory mechanism is proposed. We show that the parameter values emerging from the evolutionary simulations are near optimal; evolution takes advantage of seemingly harmful features of the environment to maximize the agent's foraging efficiency. The accessibility of evolved ANNs for a detailed inspection, together with the resemblance of some of the results to known findings from neurobiology, places evolved ANNs as an excellent candidate model for the study of structure and function relationship in complex nervous systems.

Similarity in perception: a window to brain organization.

This paper presents a neural model of similarity perception in identification tasks. It is based on self-organizing maps and population coding and is examined through five different identification experiments. Simulating an identification task, the neural model generates a confusion matrix that can be compared directly with that of human subjects. The model achieves a fairly accurate match with the pertaining experimental data both during training and thereafter. To achieve this fit, we find that the entire activity in the network should decline while learning the identification task, and that the population encoding of the specific stimuli should become sparse as the network organizes. Our results, thus, suggest that a self-organizing neural model employing population coding can account for identification processing while suggesting computational constraints on the underlying cortical networks.

NMDA receptor delayed maturation and schizophrenia.

This paper presents the hypothesis that NMDA receptor delayed maturation (NRDM) may lead to the pathogenesis of schizophrenic psychotic symptoms. This hypothesis is further analyzed in the language of a neural modeling formulation. This formulation points to a possible chain of pathological events, leading from molecular-level NRDM to over-increased synaptic plasticity, and to the formation of pathological attractors, a putative macroscopic-level correlate of schizophrenic positive symptoms. The relations of the NRDM hypothesis to other alterations which are assumed to take place in schizophrenia are discussed, together with possible ways to test this hypothesis.

Reinforcement-driven dimensionality reduction--a model for information processing in the basal ganglia.

Although anatomical studies of the basal ganglia show the existence of extensive convergence and lateral inhibitory connections, physiological studies failed to show correlated neural activity or lateral interaction in these nuclei. These seemingly contradictory results could be explained with a model in which the basal ganglia reduce the dimensionality of cortical information using optimal extraction methods. Simulations of this model predict a transient change in the efficacy of the feed-forward and lateral synapses following changes in reinforcement signal, causing an increase in correlated firing rates. This process ultimately restores the steady-state situation with diminished efficacy of lateral inhibition and no correlation of firing. Our experimental results confirm the model's predictions: rate correlations show a drastic decrease between the input stage (cortex) and output stage (pallidum). Moreover, preliminary analysis revealed that pallidal correlations show a transient increase following discrepancies between the animal's predictions and reality. We therefore propose that by using a reinforcement-driven dimensionality reduction process the basal ganglia achieve efficient extraction of cortical salient information that may then be used by the frontal cortex for execution and planning of forthcoming actions.

Neuronal regulation: A mechanism for synaptic pruning during brain maturation.

Human and animal studies show that mammalian brains undergo massive synaptic pruning during childhood, losing about half of the synapses by puberty. We have previously shown that maintaining the network performance while synapses are deleted requires that synapses be properly modified and pruned, with the weaker synapses removed. We now show that neuronal regulation, a mechanism recently observed to maintain the average neuronal input field of a postsynaptic neuron, results in a weight-dependent synaptic modification. Under the correct range of the degradation dimension and synaptic upper bound, neuronal regulation removes the weaker synapses and judiciously modifies the remaining synapses. By deriving optimal synaptic modification functions in an excitatory-inhibitory network, we prove that neuronal regulation implements near-optimal synaptic modification and maintains the performance of a network undergoing massive synaptic pruning. These findings support the possibility that neural regulation complements the action of Hebbian synaptic changes in the self-organization of the developing brain.

From parallel to serial processing: a computational study of visual search.

A novel computational model of a preattentive system performing visual search is presented. The model processes displays of lines, reproduced from Wolfe, Friedman-Hill, Stewart, and O'Connell's (1992) and Treisman and Sato's (1990) visual-search experiments. The response times measured in these experiments suggest that some of the displays are searched serially, whereas others are scanned in parallel. Our neural network model operates in two phases. First, the visual displays are compressed via standard methods (principal component analysis), to overcome assumed biological capacity limitations. Second, the compressed representations are further processed to identify a target in the display. The model succeeds in fast detection of targets in experimentally labeled parallel displays, but fails with serial ones. Analysis of the compressed internal representations reveals that compressed parallel displays contain global information that enables instantaneous target detection. However, in representations of serial displays, this global information is obscure, and hence, a target detection system should resort to a serial, attentional scan of local features across the display. Our analysis provides a numerical criterion that is strongly correlated with the experimental response time slopes and enables us to reformulate Duncan and Humphreys's (1989) search surface, using precise quantitative measures. Our findings provide further insight into the important debate concerning the dichotomous versus continuous views of parallel/serial visual search.

Associative memory in a multi-modular network.

Recent imaging studies suggest that object knowledge is stored in the brain as a distributed network of many cortical areas. Motivated by these observations, we study a multimodular associative memory network, whose functional goal is to store patterns with different coding levels--patterns that vary in the number of modules in which they are encoded. We show that in order to accomplish this task, synaptic inputs should be segregated into intramodular projections and intermodular projections, with the latter undergoing additional nonlinear dendritic processing. This segregation makes sense anatomically if the intermodular projections represent distal synaptic connections on apical dendrites. It is then straightforward to show that memories encoded in more modules are more resilient to focal afferent damage. Further hierarchical segregation of intermodular connections on the dendritic tree improves this resilience, allowing memory retrieval from input to just one of the modules in which it is encoded.

Pathogenic mechanisms in ischemic damage: a computational study.

The pathogenesis of penumbral tissue infarction during acute ischemic stroke is controversial. This peri-infarct tissue may subsequently die, or survive and recuperate, and its preservation has been a prime goal of recent therapeutic trials in acute stroke. Two major hypotheses currently under consideration are that penumbral tissue is recruited into an infarct by cortical spreading depression (CSD) waves, or by a non-wave self-propagating process such as glutamate excitotoxicity (GE). Careful experimental attempts to discriminate between these two hypotheses have so far been quite ambiguous. Using a computational metabolic model of acute focal stroke we show here that the spatial patterns of tissue damage arising from artificially induced foci of infarction having specific geometric shapes are inherently different. This is due to the distinct propagation characteristics underlying self-regenerating waves and non-wave diffusional processes. The experimental testing of these predicted spatial patterns of damage may help determine the relative contributions of the two pathological mechanisms hypothesized for ischemic tissue damage.

Synaptic pruning in development: a computational account.

Research with humans and primates shows that the developmental course of the brain involves synaptic overgrowth followed by marked selective pruning. Previous explanations have suggested that this intriguing, seemingly wasteful phenomenon is utilized to remove, "erroneous" synapses. We prove that this interpretation is wrong if synapses are Hebbian. Under limited metabolic energy resources restricting the amount and strength of synapses, we show that memory performance is maximized if synapses are first overgrown and then pruned following optimal "minimal-value" deletion. This optimal strategy leads to interesting insights concerning childhood amnesia.