ABSTRACTWhile leeches in the genusHirudohave 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. Bioinformatic analyses identified 2,812 putative genes whose expression differed significantly among the samples. These genes clustered into 7 groups which could be associated with one or more of the identified cell types. We verified predicted expression patterns throughin situhybridization 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 twopiezogenes, two of ~65deg/enacgenes, and one of at least 16transient receptor potential(trp) genes as prime candidates for involvement in sensory transduction in the three distinct classes of leech mechanosensory neurons.