Network analysis reveals age- and virus-specific circuits in nasal epithelial cells of extremely premature infants.

BACKGROUND AND OBJECTIVES: Viral respiratory infections significantly affect young children, particularly extremely premature infants, resulting in high hospitalization rates and increased health-care burdens. Nasal epithelial cells, the primary defense against respiratory infections, are vital for understanding nasal immune responses and serve as a promising target for uncovering underlying molecular and cellular mechanisms. METHODS: Using a trans-well pseudostratified nasal epithelial cell system, we examined age-dependent developmental differences and antiviral responses to influenza A and respiratory syncytial virus through systems biology approaches. RESULTS: Our studies revealed differences in innate-receptor repertoires, distinct developmental pathways, and differentially connected antiviral network circuits between neonatal and adult nasal epithelial cells. Consensus network analysis identified unique and shared cellular-viral networks, emphasizing highly relevant virus-specific pathways, independent of viral replication kinetics. CONCLUSION: This research highlights the importance of nasal epithelial cells in innate antiviral immune responses and offers crucial insights that allow for a deeper understanding of age-related differences in nasal epithelial cell immunity following respiratory virus infections.

Genetic variants and down-regulation of CACNA1H in pheochromocytoma.

Pheochromocytoma (PCC) and abdominal paraganglioma (aPGL) (together abbreviated PPGL) frequently present with an underlying genetic event in a PPGL driver gene, and additional susceptibility genes are anticipated. Here we re-analysed whole-exome sequencing data for PCC patients and identified two patients with rare missense variants in the calcium voltage-gated channel subunit 1H gene (CACNA1H). CACNA1H variants were also found in the clinical setting in PCC patients using targeted sequencing and from analysis of The Cancer Genome Atlas database. In total, CACNA1H variants were found in six PCC cases. Three of these were constitutional and two are known to have functional consequences on hormone production and gene expression in primary aldosteronism and aldosterone-producing adrenocortical adenoma. In general, PPGL exhibited reduced CACNA1H mRNA expression as compared to normal adrenal. Immunohistochemistry showed strong CACNA1H (CaV3.2) staining in adrenal medulla while PPGL typically had weak or negative staining. Reduced CACNA1H gene expression was especially pronounced in PCC as compared to aPGL and in PPGL with Cluster 2 kinase signalling phenotype. Furthermore, CACNA1H levels correlated with HIF1A and HIF2A. Moreover, TCGA data revealed a correlation between CACNA1H methylation density and gene expression. Expression of rCacna1h in PC12 cells induced differential protein expression profiles determined by mass spectrometry as well as a shift in the membrane potential where maximum calcium currents were observed as determined by electrophysiology. The findings suggest the involvement of CACNA1H/CaV3.2 in pheochromocytoma development and establish a potential link between the aetiology of adrenomedullary and adrenocortical tumor development.

A Sample Preparation Procedure for Isobaric Labeling-Based Single-Cell Proteomics.

Mass spectrometry-based proteomics has traditionally been limited by the amount of input material for analysis. Single-cell proteomics has emerged as a challenging discipline due to the ultra-high sensitivity required. Isobaric labeling-based multiplex strategies with a carrier proteome offer an approach to overcome the sensitivity limitations. Following this as the basic strategy, we show here the general workflow for preparing cells for single-cell mass spectrometry-based proteomics. This protocol can also be applied to manually isolated cells when large cells, such as cardiomyocytes, are difficult to isolate properly with conventional fluorescence-activated cell sorting (FACS) sorter methods.

EquiCP: Targeted Single-Cell Proteomics by Mass Spectrometry with Isobaric Labeled Multiplexing.

Nontargeted single-cell proteomics analysis by mass spectrometry with sample multiplexing utilizing isobaric labeling is often performed using a carrier proteome. The presented protocol describes a targeted approach that replaces the carrier proteome with a set of synthetic peptides from selected proteins, which improves the identification and quantification of these proteins in single human cells.

A multi-platform analysis of human gingival crevicular fluid reveals ferroptosis as a relevant regulated cell death mechanism during the clinical progression of periodontitis.

Ferroptosis is implicated in the pathogenesis of numerous chronic-inflammatory diseases, yet its association with progressive periodontitis remains unexplored. To investigate the involvement and significance of ferroptosis in periodontitis progression, we assessed sixteen periodontitis-diagnosed patients. Disease progression was clinically monitored over twelve weeks via weekly clinical evaluations and gingival crevicular fluid (GCF) collection was performed for further analyses. Clinical metrics, proteomic data, in silico methods, and bioinformatics tools were combined to identify protein profiles linked to periodontitis progression and to explore their potential connection with ferroptosis. Subsequent western blot analyses validated key findings. Finally, a single-cell RNA sequencing (scRNA-seq) dataset (GSE164241) for gingival tissues was analyzed to elucidate cellular dynamics during periodontitis progression. Periodontitis progression was identified as occurring at a faster rate than traditionally thought. GCF samples from progressing and non-progressing periodontal sites showed quantitative and qualitatively distinct proteomic profiles. In addition, specific biological processes and molecular functions during progressive periodontitis were revealed and a set of hub proteins, including SNCA, CA1, HBB, SLC4A1, and ANK1 was strongly associated with the clinical progression status of periodontitis. Moreover, we found specific proteins - drivers or suppressors - associated with ferroptosis (SNCA, FTH1, HSPB1, CD44, and GCLC), revealing the co-occurrence of this specific type of regulated cell death during the clinical progression of periodontitis. Additionally, the integration of quantitative proteomic data with scRNA-seq analysis suggested the susceptibility of fibroblasts to ferroptosis. Our analyses reveal proteins and processes linked to ferroptosis for the first time in periodontal patients, which offer new insights into the molecular mechanisms of progressive periodontal disease. These findings may lead to novel diagnostic and therapeutic strategies.

Proteomic analysis of the urothelial cancer landscape.

Urothelial bladder cancer (UC) has a wide tumor biological spectrum with challenging prognostic stratification and relevant therapy-associated morbidity. Most molecular classifications relate only indirectly to the therapeutically relevant protein level. We improve the pre-analytics of clinical samples for proteome analyses and characterize a cohort of 434 samples with 242 tumors and 192 paired normal mucosae covering the full range of UC. We evaluate sample-wise tumor specificity and rank biomarkers by target relevance. We identify robust proteomic subtypes with prognostic information independent from histopathological groups. In silico drug prediction suggests efficacy of several compounds hitherto not in clinical use. Both in silico and in vitro data indicate predictive value of the proteomic clusters for these drugs. We underline that proteomics is relevant for personalized oncology and provide abundance and tumor specificity data for a large part of the UC proteome ( www.cancerproteins.org ).

A cell autonomous regulator of neuronal excitability modulates tau in Alzheimer's disease vulnerable neurons.

Neurons from layer II of the entorhinal cortex (ECII) are the first to accumulate tau protein aggregates and degenerate during prodromal Alzheimer's disease. Gaining insight into the molecular mechanisms underlying this vulnerability will help reveal genes and pathways at play during incipient stages of the disease. Here, we use a data-driven functional genomics approach to model ECII neurons in silico and identify the proto-oncogene DEK as a regulator of tau pathology. We show that epigenetic changes caused by Dek silencing alter activity-induced transcription, with major effects on neuronal excitability. This is accompanied by the gradual accumulation of tau in the somatodendritic compartment of mouse ECII neurons in vivo, reactivity of surrounding microglia, and microglia-mediated neuron loss. These features are all characteristic of early Alzheimer's disease. The existence of a cell-autonomous mechanism linking Alzheimer's disease pathogenic mechanisms in the precise neuron type where the disease starts provides unique evidence that synaptic homeostasis dysregulation is of central importance in the onset of tau pathology in Alzheimer's disease.

Cancer cell immunity-related protein co-expression networks are associated with early-stage solid-predominant lung adenocarcinoma.

Background: Solid-predominant lung adenocarcinoma (SPA), which is one of the high-risk subtypes with poor prognosis and unsatisfactory response to chemotherapy and targeted therapy in lung adenocarcinoma, remains molecular profile unclarified. Weighted correlation network analysis (WGCNA) was used for data mining, especially for studying biological networks based on pairwise correlations between variables. This study aimed to identify disease-related protein co-expression networks associated with early-stage SPA. Methods: We assessed cancerous cells laser-microdissected from formalin-fixed paraffin-embedded (FFPE) tissues of a SPA group (n = 5), referencing a low-risk subtype, a lepidic predominant subtype group (LPA) (n = 4), and another high-risk subtype, micropapillary predominant subtype (MPA) group (n = 3) and performed mass spectrometry-based proteomic analysis. Disease-related co-expression networks associated with the SPA subtype were identified by WGCNA and their upstream regulators and causal networks were predicted by Ingenuity Pathway Analysis. Results: Among the forty WGCNA network modules identified, two network modules were found to be associated significantly with the SPA subtype. Canonical enriched pathways were highly associated with cellular growth, proliferation, and immune response. Upregulated HLA class I molecules HLA-G and HLA-B implicated high mutation burden and T cell activation in the SPA subtype. Upstream analysis implicated the involvement of highly activated oncogenic regulators, MYC, MLXIPL, MYCN, the redox master regulator NFE2L2, and the highly inhibited LARP1, leading to oncogenic IRES-dependent translation, and also regulators of the adaptive immune response, including highly activated IFNG, TCRD, CD3-TCR, CD8A, CD8B, CD3, CD80/CD86, and highly inhibited LILRB2. Interestingly, the immune checkpoint molecule HLA-G, which is the counterpart of LILRB2, was highly expressed characteristically in the SPA subtype and might be associated with antitumor immunity. Conclusion: Our findings provide a disease molecular profile based on protein co-expression networks identified for the high-risk solid predominant adenocarcinoma, which will help develop future therapeutic strategies.

Proteomic workflows for deep phenotypic profiling of 3D organotypic liver models.

Organotypic human tissue models constitute promising systems to facilitate drug discovery and development. They allow to maintain native cellular phenotypes and functions, which enables long-term pharmacokinetic and toxicity studies, as well as phenotypic screening. To trace relevant phenotypic changes back to specific targets or signaling pathways, comprehensive proteomic profiling is the gold-standard. A multitude of proteomic workflows have been applied on 3D tissue models to quantify their molecular phenotypes; however, their impact on analytical results and biological conclusions in this context has not been evaluated. The performance of twelve mass spectrometry-based global proteomic workflows that differed in the amount of cellular input, lysis protocols and quantification methods was compared for the analysis of primary human liver spheroids. Results differed majorly between protocols in the total number and subcellular compartment bias of identified proteins, which is particularly relevant for the reliable quantification of transporters and drug metabolizing enzymes. Using a model of metabolic dysfunction-associated steatotic liver disease, we furthermore show that critical disease pathways are robustly identified using a standardized high throughput-compatible workflow based on thermal lysis, even using only individual spheroids (1500 cells) as input. The results increase the applicability of proteomic profiling to phenotypic screens in organotypic microtissues and provide a scalable platform for deep phenotyping from limited biological material.

Epigenetic modulators link mitochondrial redox homeostasis to cardiac function in a sex-dependent manner.

While excessive production of reactive oxygen species (ROS) is a characteristic hallmark of numerous diseases, clinical approaches that ameliorate oxidative stress have been unsuccessful. Here, utilizing multi-omics, we demonstrate that in cardiomyocytes, mitochondrial isocitrate dehydrogenase (IDH2) constitutes a major antioxidative defense mechanism. Paradoxically reduced expression of IDH2 associated with ventricular eccentric hypertrophy is counterbalanced by an increase in the enzyme activity. We unveil redox-dependent sex dimorphism, and extensive mutual regulation of the antioxidative activities of IDH2 and NRF2 by a feedforward network that involves 2-oxoglutarate and L-2-hydroxyglutarate and mediated in part through unconventional hydroxy-methylation of cytosine residues present in introns. Consequently, conditional targeting of ROS in a murine model of heart failure improves cardiac function in sex- and phenotype-dependent manners. Together, these insights may explain why previous attempts to treat heart failure with antioxidants have been unsuccessful and open new approaches to personalizing and, thereby, improving such treatment.

Ultralight Ultrafast Enzymes.

Inorganic materials depleted of heavy stable isotopes are known to deviate strongly in some physicochemical properties from their isotopically natural counterparts. Here we explored for the first time the effect of simultaneous depletion of the heavy carbon, hydrogen, oxygen and nitrogen isotopes on the bacterium E. coli and the enzymes expressed in it. Bacteria showed faster growth, with most proteins exhibiting higher thermal stability, while for recombinant enzymes expressed in depleted media, faster kinetics was discovered. At room temperature, luciferase, thioredoxin and dihydrofolate reductase and Pfu DNA polymerase showed up to a 250 % increase in activity compared to the native counterparts, with an additional ∼50 % increase at 10 °C. Diminished conformational and vibrational entropy is hypothesized to be the cause of the accelerated kinetics. Ultralight enzymes may find an application where extreme reaction rates are required.

Corrigendum: Brain endothelial cells exposure to malaria parasites links type I interferon signalling to antigen presentation, immunoproteasome activation, endothelium disruption, and cellular metabolism.

[This corrects the article DOI: 10.3389/fimmu.2023.1149107.].

Autoreactive B cells against malondialdehyde-induced protein cross-links are present in the joint, lung, and bone marrow of rheumatoid arthritis patients.

Autoantibodies to malondialdehyde (MDA) proteins constitute a subset of anti-modified protein autoantibodies in rheumatoid arthritis (RA), which is distinct from citrulline reactivity. Serum anti-MDA IgG levels are commonly elevated in RA and correlate with disease activity, CRP, IL6, and TNF-α. MDA is an oxidation-associated reactive aldehyde that together with acetaldehyde mediates formation of various immunogenic amino acid adducts including linear MDA-lysine, fluorescent malondialdehyde acetaldehyde (MAA)-lysine, and intramolecular cross-linking. We used single-cell cloning, generation of recombinant antibodies (n = 356 from 25 donors), and antigen-screening to investigate the presence of class-switched MDA/MAA+ B cells in RA synovium, bone marrow, and bronchoalveolar lavage. Anti-MDA/MAA+ B cells were found in bone marrow plasma cells of late disease and in the lung of both early disease and risk-individuals and in different B cell subsets (memory, double negative B cells). These were compared with previously identified anti-MDA/MAA from synovial memory and plasma cells. Seven out of eight clones carried somatic hypermutations and all bound MDA/MAA-lysine independently of protein backbone. However, clones with somatic hypermutations targeted MAA cross-linked structures rather than MDA- or MAA-hapten, while the germline-encoded synovial clone instead bound linear MDA-lysine in proteins and peptides. Binding patterns were maintained in germline converted clones. Affinity purification of polyclonal anti-MDA/MAA from patient serum revealed higher proportion of anti-MAA versus anti-MDA compared to healthy controls. In conclusion, IgG anti-MDA/MAA show distinct targeting of different molecular structures. Anti-MAA IgG has been shown to promote bone loss and osteoclastogenesis in vivo and may contribute to RA pathogenesis.

Systems-level temporal immune-metabolic profile in Crimean-Congo hemorrhagic fever virus infection.

Crimean-Congo hemorrhagic fever (CCHF) caused by CCHF virus (CCHFV) is one of the epidemic-prone diseases prioritized by the World Health Organisation as public health emergency with an urgent need for accelerated research. The trajectory of host response against CCHFV is multifarious and remains unknown. Here, we reported the temporal spectrum of pathogenesis following the CCHFV infection using genome-wide blood transcriptomics analysis followed by advanced systems biology analysis, temporal immune-pathogenic alterations, and context-specific progressive and postinfection genome-scale metabolic models (GSMM) on samples collected during the acute (T0), early convalescent (T1), and convalescent-phase (T2). The interplay between the retinoic acid-inducible gene-I-like/nucleotide-binding oligomerization domain-like receptor and tumor necrosis factor signaling governed the trajectory of antiviral immune responses. The rearrangement of intracellular metabolic fluxes toward the amino acid metabolism and metabolic shift toward oxidative phosphorylation and fatty acid oxidation during acute CCHFV infection determine the pathogenicity. The upregulation of the tricarboxylic acid cycle during CCHFV infection, compared to the noninfected healthy control and between the severity groups, indicated an increased energy demand and cellular stress. The upregulation of glycolysis and pyruvate metabolism potentiated energy generation through alternative pathways associated with the severity of the infection. The downregulation of metabolic processes at the convalescent phase identified by blood cell transcriptomics and single-cell type proteomics of five immune cells (CD4+ and CD8+ T cells, CD14+ monocytes, B cells, and NK cells) potentially leads to metabolic rewiring through the recovery due to hyperactivity during the acute phase leading to post-viral fatigue syndrome.

Toward Single Bacterium Proteomics.

Bacteria are orders of magnitude smaller than mammalian cells, and while single cell proteomics (SCP) currently detects and quantifies several thousands of proteins per mammalian cell, it is not clear whether conventional SCP methods will be suitable for bacteria. Here we report on the first successful attempt to detect proteins from individual Escherichia coli bacteria, with validation of our findings by comparison with two bacteria samples and bulk proteomics data. Data are available via ProteomeXchange with the identifier PXD043473.

3D microperfusion of mesoscale human microphysiological liver models improves functionality and recapitulates hepatic zonation.

Hepatic in vitro models that accurately replicate phenotypes and functionality of the human liver are needed for applications in toxicology, pharmacology and biomedicine. Notably, it has become clear that liver function can only be sustained in 3D culture systems at physiologically relevant cell densities. Additionally, drug metabolism and drug-induced cellular toxicity often follow distinct spatial micropatterns of the metabolic zones in the liver acinus, calling for models that capture this zonation. We demonstrate the manufacture of accurate liver microphysiological systems (MPS) via engineering of 3D stereolithography printed hydrogel chips with arrays of diffusion open synthetic vasculature channels at spacings approaching in vivo capillary distances. Chip designs are compatible with seeding of cell suspensions or preformed liver cell spheroids. Importantly, primary human hepatocytes (PHH) and hiPSC-derived hepatocyte-like cells remain viable, exhibit improved molecular phenotypes compared to isogenic monolayer and static spheroid cultures and form interconnected tissue structures over the course of multiple weeks in perfused culture. 3D optical oxygen mapping of embedded sensor beads shows that the liver MPS recapitulates oxygen gradients found in the acini, which translates into zone-specific acet-ami-no-phen toxicity patterns. Zonation, here naturally generated by high cell densities and associated oxygen and nutrient utilization along the flow path, is also documented by spatial proteomics showing increased concentration of periportal- versus perivenous-associated proteins at the inlet region and vice versa at the outlet region. The presented microperfused liver MPS provides a promising platform for the mesoscale culture of human liver cells at phenotypically relevant densities and oxygen exposures. STATEMENT OF SIGNIFICANCE: A full 3D tissue culture platform is presented, enabled by massively parallel arrays of high-resolution 3D printed microperfusion hydrogel channels that functionally mimics tissue vasculature. The platform supports long-term culture of liver models with dimensions of several millimeters at physiologically relevant cell densities, which is difficult to achieve with other methods. Human liver models are generated from seeded primary human hepatocytes (PHHs) cultured for two weeks, and from seeded spheroids of hiPSC-derived human liver-like cells cultured for two months. Both model types show improved functionality over state-of-the-art 3D spheroid suspensions cultured in parallel. The platform can generate physiologically relevant oxygen gradients driven by consumption rather than supply, which was validated by visualization of embedded oxygen-sensitive microbeads, which is exploited to demonstrate zonation-specific toxicity in PHH liver models.

Analysis of local extracellular matrix identifies different aetiologies behind bicuspid and tricuspid aortic valve degeneration and suggests therapies.

Aortic valve degeneration (AVD) is a life-threatening condition that has no medical treatment and lacks individual therapies. Although extensively studied with standard approaches, aetiologies behind AVD are unclear. We compared abundances of extracellular matrix (ECM) proteins from excised valve tissues of 88 patients with isolated AVD of normal tricuspid (TAV) and congenital bicuspid aortic valves (BAV), quantified more than 1400 proteins per ECM sample by mass spectrometry, and demonstrated that local ECM preserves molecular cues of the pathophysiological processes. The BAV ECM showed enrichment with fibrosis markers, namely Tenascin C, Osteoprotegerin, and Thrombospondin-2. The abnormal physical stress on BAV may cause a mechanical injury leading to a continuous Tenascin C-driven presence of myofibroblasts and persistent fibrosis. The TAV ECM exhibited enrichment with Annexin A3 (p = 1.1 × 10-16 and the fold change 6.5) and a significant deficit in proteins involved in high-density lipid metabolism. These results were validated by orthogonal methods. The difference in the ECM landscape suggests distinct aetiologies between AVD of BAV and TAV; warrants different treatments of the patients with BAV and TAV; elucidates the molecular basis of AVD; and implies possible new therapeutic approaches. Our publicly available database (human_avd_ecm.surgsci.uu.se) is a rich source for medical doctors and researchers who are interested in AVD or heart ECM in general. Systematic proteomic analysis of local ECM using the methods described here may facilitate future studies of various tissues and organs in development and disease.

Proteomic landscape of astrocytes and pericytes infected with HIV/SARS-CoV-2 mono/co-infection, impacting on neurological complications.

Background: Although most individuals recover from coronavirus disease 2019 (COVID-19) within a few weeks, some people continue to experience a wide range of symptoms known as post-acute sequelae of SARS-CoV-2 (PASC) or long COVID. Majority of patients with PASC develop neurological disorders like brain fog, fatigue, mood swings, sleep disorders, loss of smell and test among others collectively called neuro-PASC. While the people living with HIV (PWH) do not have a higher risk of developing severe disease and mortality/morbidity due to COVID-19. As a large section of PWH suffered from HIV-associated neurocognitive disorders (HAND), it is essential to understand the impact of neuro-PASC on people with HAND. In pursuit of this, we infected HIV/SARS-CoV-2 alone or together in primary human astrocytes and pericytes and performed proteomics to understand the impact of co-infection in the central nervous system. Methods: Primary human astrocytes and pericytes were infected with SARS-CoV-2 or HIV or HIV + SARS-CoV-2. The concentration of HIV and SARS-CoV-2 genomic RNA in the culture supernatant was quantified using reverse transcriptase quantitative real time polymerase chain reaction (RT-qPCR). This was followed by a quantitative proteomics analysis of mock, HIV, SARS-CoV-2, and HIV + SARS-CoV-2 infected astrocytes and pericytes to understand the impact of the virus in CNS cell types. Results: Both healthy and HIV-infected astrocytes and pericytes support abortive/low level of SARS-CoV-2 replication. In both mono-infected and co-infected cells, we observe a modest increase in the expression of SARS-CoV-2 host cell entry factors (ACE2, TMPRSS2, NRP1, and TRIM28) and inflammatory mediators (IL-6, TNF-α, IL-1ß and IL-18). Quantitative proteomic analysis has identified uniquely regulated pathways in mock vs SARS-CoV-2, mock vs HIV + SARS-CoV-2, and HIV vs HIV + SARS-CoV-2 infected astrocytes and pericytes. The gene set enrichment analysis revealed that the top ten enriched pathways are linked to several neurodegenerative disorders, including Alzheimer's disease, Parkinson's disease, Huntington's disease, and amyotrophic lateral sclerosis. Conclusions: Our study emphasizes the significance of long-term monitoring of patients co-infected with HIV and SARS-CoV-2 to detect and understand the development of neurological abnormalities. By unraveling the molecular mechanisms involved, we can identify potential targets for future therapeutic interventions.

Distribution, cellular localization, and colocalization of several peptide neurotransmitters in the central nervous system of Aplysia.

Neuropeptides are widely used as neurotransmitters in vertebrates and invertebrates. In vertebrates, a detailed understanding of their functions as transmitters has been hampered by the complexity of the nervous system. The marine mollusk Aplysia, with a simpler nervous system and many large, identified neurons, presents several advantages for addressing this question and has been used to examine the roles of tens of peptides in behavior. To screen for other peptides that might also play roles in behavior, we observed immunoreactivity in individual neurons in the central nervous system of adult Aplysia with antisera raised against the Aplysia peptide FMRFamide and two mammalian peptides that are also found in Aplysia, cholecystokinin (CCK) and neuropeptide Y (NPY), as well as serotonin (5HT). In addition, we observed staining of individual neurons with antisera raised against mammalian somatostatin (SOM) and peptide histidine isoleucine (PHI). However, genomic analysis has shown that these two peptides are not expressed in the Aplysia nervous system, and we have therefore labeled the unknown peptides stained by these two antibodies as XSOM and XPHI There was an area at the anterior end of the cerebral ganglion that had staining by antisera raised against many different transmitters, suggesting that this may be a modulatory region of the nervous system. There was also staining for XSOM and, in some cases, FMRFamide in the bag cell cluster of the abdominal ganglion. In addition, these and other studies have revealed a fairly high degree of colocalization of different neuropeptides in individual neurons, suggesting that the peptides do not just act independently but can also interact in different combinations to produce complex functions. The simple nervous system of Aplysia is advantageous for further testing these ideas.