An annotated CNS transcriptome of the medicinal leech, Hirudo verbana: De novo sequencing to characterize genes associated with nervous system activity

Molecular Insights Into the Sensory Adaption of the Cave-Dwelling Leech Sinospelaeobdella wulingensis to the Karst Cave Environment

Karst caves are a unique environment significantly different from the external environment; adaptation of cave-dwelling animals to the cave environment is often accompanied by shifts in the sensory systems. Aquatic and terrestrial leeches have been found in the karst caves. In this study, we conducted a transcriptome analysis on the cave-dwelling leech Sinospelaeobdella wulingensis. A total of 29,286 unigenes were obtained by assembling the clean reads, and only 395 genes are differentially expressed in winter and summer samples. Two piezo-type mechanosensitive ion channels (Piezos), eight transient receptor potential channels (TRPs), and six ionotropic glutamate receptors (iGluRs) were identified in the transcriptome. These channels/receptors are transmembrane proteins sharing conserved structural features in the respective protein families. SwPiezo1 shares high identity with Piezos in non-caving leeches. SwiGluRs are conserved in protein sequence and share high identities with homologous proteins in other leeches. In contrast, SwTRPs belong to different subfamilies and share diverse identities with TRPs in other species. Gene expression analysis showed that two SwPiezos, five SwTRPs, and one SwiGluR are abundantly expressed in both winter and summer samples. These results suggest that SwPiezos, SwTRPs, and SwiGluRs are candidate sensory channels/receptors that may have roles in mechanosensory and chemosensory systems. High expression levels of Piezo and TRP genes imply a mechanosensory adaptation of S. wulingensis to the hanging living style in caves. Furthermore, enrichment of sensory genes in the oral sucker indicates the important role of this tissue in response to environmental stimuli. Similar gene expression profiles in winter and summer samples imply a stable physiological status of S. wulingensis in the cave environment.

Acceleration of genome rearrangement in clitellate annelids

Comparisons of multiple metazoan genomes have revealed the existence of ancestral linkage groups (ALGs), genomic scaffolds sharing sets of orthologous genes that have been inherited from ancestral animals for hundreds of millions of years (Simakov et al. 2022; Schultz et al. 2023) These ALGs have persisted across major animal taxa including Cnidaria, Deuterostomia, Ecdysozoa and Spiralia. Notwithstanding this general trend of chromosome-scale conservation, ALGs have been obliterated by extensive genome rearrangements in certain groups, most notably including Clitellata (oligochaetes and leeches), a group of easily overlooked invertebrates that is of tremendous ecological, agricultural and economic importance (Charles 2019; Barrett 2016). To further investigate these rearrangements, we have undertaken a comparison of 12 clitellate genomes (including four newly sequenced species) and 11 outgroup representatives. We show that these rearrangements began at the base of the Clitellata (rather than progressing gradually throughout polychaete annelids), that the inter-chromosomal rearrangements continue in several clitellate lineages and that these events have substantially shaped the evolution of the otherwise highly conserved Hox cluster.

Transcriptomic profiling and the first spatial expression analysis of candidate genes in the salivary gland of the East Asian medicinal leech, Hirudo nipponia

Hirudo nipponia, a blood-sucking leech native to East Asia, possesses a rich repertoire of active ingredients in its saliva, showcasing significant medical potential due to its anticoagulant, anti-inflammatory, and antibacterial effects against human diseases. Despite previous studies on the transcriptomic and proteomic characteristics of leech saliva, which have identified medicinal compounds, our knowledge of tissue-specific transcriptomes and their spatial expression patterns remains incomplete. In this study, we conducted an extensive transcriptomic profiling of the salivary gland tissue in H. nipponia based on de novo assemblies of tissue-specific transcriptomes from the salivary gland, teeth, and general head region. Through gene ontology (GO) analysis and hierarchical clustering, we discovered a novel set of anti-coagulant factors-i.e., Hni-Antistasin, Hni-Ghilanten, Hni-Bdellin, Hni-Hirudin-as well as a previously unrecognized immune-related gene, Hni-GLIPR1 and uncharacterized salivary gland specific transcripts. By employing in situ hybridization, we provided the first visualization of gene expression sites within the salivary gland of H. nipponia. Our findings expand on our understanding of transcripts specifically expressed in the salivary gland of blood-sucking leeches, offering valuable resources for the exploration of previously unidentified substances with medicinal applications.

Cell anatomy and network input explain differences within but not between leech touch cells at two different locations

Mechanosensory cells in the leech share several common features with mechanoreceptors in the human glabrous skin. Previous studies showed that the six T (touch) cells in each body segment of the leech are highly variable in their responses to somatic current injection and change their excitability over time. Here, we investigate three potential reasons for this variability in excitability by comparing the responses of T cells at two soma locations (T2 and T3): (1) Differential effects of time-dependent changes in excitability, (2) divergent synaptic input from the network, and (3) different anatomical structures. These hypotheses were explored with a combination of electrophysiological double recordings, 3D reconstruction of neurobiotin-filled cells, and compartmental model simulations. Current injection triggered significantly more spikes with shorter latency and larger amplitudes in cells at soma location T2 than at T3. During longer recordings, cells at both locations increased their excitability over time in the same way. T2 and T3 cells received the same amount of synaptic input from the unstimulated network, and the polysynaptic connections between both T cells were mutually symmetric. However, we found a striking anatomical difference: While in our data set all T2 cells innervated two roots connecting the ganglion with the skin, 50% of the T3 cells had only one root process. The sub-sample of T3 cells with one root process was significantly less excitable than the T3 cells with two root processes and the T2 cells. To test if the additional root process causes higher excitability, we simulated the responses of 3D reconstructed cells of both anatomies with detailed multi-compartment models. The anatomical subtypes do not differ in excitability when identical biophysical parameters and a homogeneous channel distribution are assumed. Hence, all three hypotheses may contribute to the highly variable T cell responses, but none of them is the only factor accounting for the observed systematic difference in excitability between cells at T2 vs. T3 soma location. Therefore, future patch clamp and modeling studies are needed to analyze how biophysical properties and spatial distribution of ion channels on the cell surface contribute to the variability and systematic differences of electrophysiological phenotypes.

Neural circuit regulation by identified modulatory projection neurons

Rhythmic behaviors (e.g., walking, breathing, and chewing) are produced by central pattern generator (CPG) circuits. These circuits are highly dynamic due to a multitude of input they receive from hormones, sensory neurons, and modulatory projection neurons. Such inputs not only turn CPG circuits on and off, but they adjust their synaptic and cellular properties to select behaviorally relevant outputs that last from seconds to hours. Similar to the contributions of fully identified connectomes to establishing general principles of circuit function and flexibility, identified modulatory neurons have enabled key insights into neural circuit modulation. For instance, while bath-applying neuromodulators continues to be an important approach to studying neural circuit modulation, this approach does not always mimic the neural circuit response to neuronal release of the same modulator. There is additional complexity in the actions of neuronally-released modulators due to: (1) the prevalence of co-transmitters, (2) local- and long-distance feedback regulating the timing of (co-)release, and (3) differential regulation of co-transmitter release. Identifying the physiological stimuli (e.g., identified sensory neurons) that activate modulatory projection neurons has demonstrated multiple “modulatory codes” for selecting particular circuit outputs. In some cases, population coding occurs, and in others circuit output is determined by the firing pattern and rate of the modulatory projection neurons. The ability to perform electrophysiological recordings and manipulations of small populations of identified neurons at multiple levels of rhythmic motor systems remains an important approach for determining the cellular and synaptic mechanisms underlying the rapid adaptability of rhythmic neural circuits.

Spatio-Temporal Changes of Extracellular Matrix (ECM) Stiffness in the Development of the Leech Hirudo verbana

The invertebrate leech Hirudo verbana represents a powerful experimental animal model for improving the knowledge about the functional interaction between the extracellular matrix (ECM) and cells within the tissue microenvironment (TME), and the key role played by ECM stiffness during development and growth. Indeed, the medicinal leech is characterized by a simple anatomical organization reproducing many aspects of the basic biological processes of vertebrates and in which a rapid spatiotemporal development is well established and easily assessed. Our results show that ECM structural organization, as well as the amount of fibrillar and non-fibrillar collagen are deeply different from hatching leeches to adult ones. In addition, the changes in ECM remodelling occurring during the different leech developmental stages, leads to a gradient of stiffness regulating both the path of migratory cells and their fates. The ability of cells to perceive and respond to changes in ECM composition and mechanics strictly depend on nuclear or cytoplasmic expression of Yes-Associated Protein 1 (YAP1), a key mediator converting mechanical signals into transcriptional outputs, expression, and activation.

A review of the bioeffects of low-intensity focused ultrasound and the benefits of a cellular approach

This review article highlights the historical developments and current state of knowledge of an important neuromodulation technology: low-intensity focused ultrasound. Because compelling studies have shown that focused ultrasound can modulate neuronal activity non-invasively, especially in deep brain structures with high spatial specificity, there has been a renewed interest in attempting to understand the specific bioeffects of focused ultrasound at the cellular level. Such information is needed to facilitate the safe and effective use of focused ultrasound to treat a number of brain and nervous system disorders in humans. Unfortunately, to date, there appears to be no singular biological mechanism to account for the actions of focused ultrasound, and it is becoming increasingly clear that different types of nerve cells will respond to focused ultrasound differentially based on the complement of their ion channels, other membrane biophysical properties, and arrangement of synaptic connections. Furthermore, neurons are apparently not equally susceptible to the mechanical, thermal and cavitation-related consequences of focused ultrasound application-to complicate matters further, many studies often use distinctly different focused ultrasound stimulus parameters to achieve a reliable response in neural activity. In this review, we consider the benefits of studying more experimentally tractable invertebrate preparations, with an emphasis on the medicinal leech, where neurons can be studied as unique individual cells and be synaptically isolated from the indirect effects of focused ultrasound stimulation on mechanosensitive afferents. In the leech, we have concluded that heat is the primary effector of focused ultrasound neuromodulation, especially on motoneurons in which we observed a focused ultrasound-mediated blockade of action potentials. We discuss that the mechanical bioeffects of focused ultrasound, which are frequently described in the literature, are less reliably achieved as compared to thermal ones, and that observations ascribed to mechanical responses may be confounded by activation of synaptically-coupled sensory structures or artifacts associated with electrode resonance. Ultimately, both the mechanical and thermal components of focused ultrasound have significant potential to contribute to the sculpting of specific neural outcomes. Because focused ultrasound can generate significant modulation at a temperature <5°C, which is believed to be safe for moderate durations, we support the idea that focused ultrasound should be considered as a thermal neuromodulation technology for clinical use, especially targeting neural pathways in the peripheral nervous system.

Transcriptional profiling of identified neurons in leech

Background

While leeches in the genus Hirudo have long been models for neurobiology, the molecular underpinnings of nervous system structure and function in this group remain largely unknown. To begin to bridge this gap, we performed RNASeq on pools of identified neurons of the central nervous system (CNS): sensory T (touch), P (pressure) and N (nociception) neurons; neurosecretory Retzius cells; and ganglia from which these four cell types had been removed.

Results

Bioinformatic analyses identified 3565 putative genes whose expression differed significantly among the samples. These genes clustered into 9 groups which could be associated with one or more of the identified cell types. We verified predicted expression patterns through in situ hybridization on whole CNS ganglia, and found that orthologous genes were for the most part similarly expressed in a divergent leech genus, suggesting evolutionarily conserved roles for these genes. Transcriptional profiling allowed us to identify candidate phenotype-defining genes from expanded gene families. Thus, we identified one of eight hyperpolarization-activated cyclic-nucleotide gated (HCN) channels as a candidate for mediating the prominent sag current in P neurons, and found that one of five inositol triphosphate receptors (IP3Rs), representing a sub-family of IP3Rs absent from vertebrate genomes, is expressed with high specificity in T cells. We also identified one of two piezo genes, two of ~ 65 deg/enac genes, and one of at least 16 transient receptor potential (trp) genes as prime candidates for involvement in sensory transduction in the three distinct classes of leech mechanosensory neurons.

Conclusions

Our study defines distinct transcriptional profiles for four different neuronal types within the leech CNS, in addition to providing a second ganglionic transcriptome for the species. From these data we identified five gene families that may facilitate the sensory capabilities of these neurons, thus laying the basis for future work leveraging the strengths of the leech system to investigate the molecular processes underlying and linking mechanosensation, cell type specification, and behavior.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07526-0.

Effects of calcium-activated potassium channel modulators on afterhyperpolarizing potentials in identified motor and mechanosensory neurons of the medicinal leech

Calcium-activated potassium (K_Ca) channels contribute to multiple neuronal properties including spike frequency and afterhyperpolarizing potentials (AHPs). K_Ca channels are classified as K_Ca1.1, K_Ca2, or K_Ca3.1 based on single-channel conductance and pharmacology. Ca²⁺-dependent AHPs in vertebrates are categorized as fast, medium, or slow. Fast and medium AHPs are generated by K_Ca1.1 and K_Ca2 channels, respectively. The K_Ca subtype responsible for slow AHPs is unclear. Prolonged, Ca²⁺-dependent AHPs have been described in several leech neurons. Unfortunately, apamin and other K_Ca blockers often prove ineffective in the leech. An alternative approach is to utilize K_Ca modulators, which alter channel sensitivity to Ca²⁺. Vertebrate K_Ca2 channels are targeted selectively by the positive modulator CyPPA and the negative modulator NS8593. Here we show that AHPs in identified motor and mechanosensory leech neurons are enhanced by CyPPA and suppressed by NS8593. Our results indicate that K_Ca2 channels underlie prolonged AHPs in these neurons and suggest that K_Ca2 modulators may serve as effective tools to explore the role of K_Ca channels in leech physiology.

Feeding, hunger, satiety and serotonin in invertebrates

The serotonergic modulation of feeding behaviour has been intensively studied in several invertebrate groups, including Arthropoda, Annelida, Nematoda and Mollusca. These studies offer comparative information on feeding regulation across divergent phyla and also provide general insights into the neural control of feeding. Specifically, model invertebrates are ideal for parsing feeding behaviour into component parts and examining the underlying mechanisms at the levels of biochemical pathways, single cells and identified neural circuitry. Research has found that serotonin is crucial during certain phases of feeding behaviour, especially movements directly underlying food intake, but inessential during other phases. In addition, while the serotonin system can be manipulated systemically in many animals, invertebrate model organisms also allow manipulations at the level of single cells and molecules, revealing limited and precise serotonergic actions. The latter highlight the importance of local versus global modulatory effects of serotonin, a potentially significant consideration for drug and pesticide design.

Purinergic modulation of neuronal gap junction circuits in the CNS of the leech

Gap junctions (GJs) are widely distributed in brains across the animal kingdom. To visualize the GJ- coupled networks of two major mechanosensory neurons in the ganglia of medicinal leeches, we injected these cells with the GJ-permeable tracer Neurobiotin. When diffusion time was limited to only 30 min, tracer coupling was highly variable for both cells, suggesting a possible modulation of GJ permeability. In invertebrates the innexins (homologs of vertebrate pannexins) form the GJs. Because extracellular adenosine triphosphate (ATP) modulates pannexin and leech innexin hemichannel permeability and is released by leech glial cells following injury, we tested the effects of bath application of ATP after the injection of Neurobiotin and observed a significant increase in the number of neurons tracer coupled to the sensory neurons. This effect required the elevation of intracellular Ca²⁺ and could be produced by bath application of caffeine. Conversely, scavenging endogenous extracellular ATP with the ATPase apyrase decreased the number of coupled cells. ATP also increased electrical conductance and tracer permeability between the bilateral Retzius neurons. This modulatory effect of ATP on GJ coupling was blocked by siRNA knockdown of a P1-like adenosine receptor. Finally, exposure of leech ganglia to extracellular ATP induced a characteristic low frequency (<0.3 Hz) rhythmic bursting activity that was roughly synchronous among multiple neurons, a behavior that was significantly attenuated by the GJ blocker octanol. These findings highlight the mediation by ATP of a robust physiological mechanism for modifying neuronal circuits by rapidly recruiting neurons into active networks and entraining synchronized bursting activity.

The medicinal leech as a valuable model for better understanding the role of a TLR4-like receptor in the inflammatory process

Despite extensive investigation focused on both the molecular characteristics and the expression level of Toll-like receptors (TLRs) during the inflammatory response in vertebrates, few data are available in the literature on the role of these proteins in invertebrate's immune response. Here, we propose the medicinal leech as a valuable model to better elucidate the role of TLR4 and its related products, such as tumor necrosis factor (TNF-α), after activation of the leech peripheral immune system with the endogenous medicinal leech recombinant allograft inflammatory factor-1 (rHmAIF-1) or with an exogenous stimulus, such as lipopolysaccharide (LPS). Our results indicate that activated macrophages (HmAIF-1⁺) and granulocytes (CD11b⁺) express both TLR4 and its coreceptor CD14. Moreover, functional studies performed by injecting a cyanobacterium selective TLR4 antagonist CyP demonstrated that only the TLR4 pathway was blocked, while the immune response caused by lipoteichoic acid (LTA) treatment is not affected. These results are consistent with literature on vertebrates, indicating that TLR4 functions as a LPS receptor while the recognition of LTA may involve other pathways.