Abstract
The stimulant methylphenidate (MPH) and the non-stimulant atomoxetine (ATX) are frequently used for the treatment of attention-deficit/hyperactivity disorder (ADHD); however, the function of these drugs in different types of brain cells and their effects on related genes remain largely unknown. To address these questions, we built a pipeline for the simultaneous examination of the activity behavior and transcriptional responses of Drosophila melanogaster at single-cell resolution following drug treatment. We selected the Drosophila with significantly increased locomotor activities (hyperactivity-like behavior) following the administration of each drug in comparison with the control (same food as the drug-treated groups with 5% sucrose, yeast, and blue food dye solution) using EasyFlyTracker. Subsequently, single cell RNA sequencing (scRNASEQ) was used to capture the transcriptome of 82,917 cells, unsupervised clustering analysis of which yielded 28 primary cell clusters representing the major cell types in adult Drosophila brain. Indeed, both neuronal and glial cells responded to MPH and ATX. Further analysis of differentially expressed genes (DEGs) revealed distinct transcriptional changes associated with these two drugs, such as two well-studied dopamine receptor genes (Dop2R and DopEcR) were responsive to MPH but not to ATX at their optimal doses, in addition to genes involved in dopamine metabolism pathways such as Syt1, Sytalpha, Syt7, and Ih in different cell types. More importantly, MPH also suppressed the expression of genes encoding other neurotransmitter receptors and synaptic signaling molecules in many cell types, especially those for Glu and GABA, while the responsive effects of ATX were much weaker. In addition to monoaminergic neuronal transmitters, other neurotransmitters have also shown a similar pattern with respect to a stronger effect associated with MPH than with ATX. Moreover, we identified four distinct glial cell subtypes responsive to the two drugs and detected a greater number of differentially expressed genes associated with ensheathing and astrocyte-like glia. Furthermore, our study provides a rich resource of candidate target genes, supported by drug set enrichment analysis (P = 2.10E-4; hypergeometric test), for the further exploration of drug repurposing. The whole list of candidates can be found at ADHDrug ( http://adhdrug.cibr.ac.cn/ ). In conclusion, we propose a fast and cost-efficient pipeline to explore the underlying molecular mechanisms of ADHD drug treatment in Drosophila brain at single-cell resolution, which may further facilitate drug repurposing applications.