Skip to main content

Effects of midazolam, pentobarbital and ketamine on the mRNA expression of ion channels in a model organism Daphnia pulex



Over the last few decades intensive studies have been carried out on the molecular targets mediating general anesthesia as well as the effects of general anesthetics. The γ-aminobutyric acid type A receptor (GABAAR) has been indicated as the primary target of general anaesthetics such as propofol, etomidate and isoflurane, and sedating drugs including benzodiazepines and barbiturates. The GABAAR is also involved in drug tolerance and dependence. However, the involvement of other ion channels is possible.


Using reverse transcription and quantitative PCR techniques, we systematically investigated changes in the mRNA levels of ion channel genes in response to exposure to midazolam, pentobarbital and ketamine in a freshwater model animal, Daphnia pulex. To retrieve the sequences of Daphnia ion channel genes, Blast searches were performed based on known human or Drosophila ion channel genes. Retrieved sequences were clustered with the maximum-likelihood method. To quantify changes in gene expression after the drug treatments for 4 hours, total RNA was extracted and reverse transcribed into cDNA and then amplified using quantitative PCR.


A total of 108 ion channel transcripts were examined, and 19, 11 and 11 of them are affected by midazolam (100 μM), pentobarbital (200 μM) and ketamine (100 μM), respectively, covering a wide variety of ion channel types. There is some degree of overlap with midazolam- and pentobarbital-induced changes in the mRNA expression profiles, but ketamine causes distinct changes in gene expression pattern.

In addition, flumazenil (10 μM) eliminates the effect of midazolam on the mRNA expression of the GABAA receptor subunit Rdl, suggesting a direct interaction between midazolam and GABAA receptors.


Recent research using high throughput technology suggests that changes in mRNA expression correlate with delayed protein expression. Therefore, the mRNA profile changes in our study may reflect the molecular targets not only in drug actions, but also in chronic drug addiction. Our data also suggest the possibility that hypnotic/anesthetic drugs are capable of altering the functions of the nervous system, as well as those non-nerve tissues with abundant ion channel expressions.

Peer Review reports


Midazolam is a benzodiazepine and widely used as an anxiolytic, anticonvulsant, sleep aid, muscle relaxant, and antipsychotic. Pentobarbital is a short-acting barbiturate that is used as a sedative and anesthetic agent. Like other barbiturates, pentobarbital produces a wide spectrum of dose-dependent effects, including sedation, hypnosis, anesthesia and finally coma. Although their use has decreased over the years because of high abuse potential, barbiturates are still being prescribed to many patients, such as epilepsy patients and people with sleep disorders. The principal mechanism of actions of benzodiazepines and barbiturates is believed to be positive allosteric modulation of the γ-aminobutyric acid (GABA) type A receptor (GABAAR) [13]. Ketamine, a rapid acting anesthetic agent and a popular drug of abuse, has diverse effects, including antidepressant action and analgesic effects on chronic pain [4]. These actions were traditionally believed to arise from the inhibition of NMDA receptors. However, NMDAR blockers, such as MK-801, fail to mimic the actions of ketamine [5, 6]. Furthermore, knockout of NR2A subunit reduces but does not eliminate the actions of ketamine [79]. Ketamine has also been shown to positively modulate the function of the cerebellar GABAA receptors containing α6 and δ subunits [10].

Chronic use of midazolam, phenobarbital and ketamine produces tolerance and physical dependence. The homeostatic theory of drug tolerance [11] claims in a modern form that the functional tolerance results from altered function or expression of proteins in a way to reduce the effects of the drugs. Based on this theory, for example, the GABAAR, which is positively modulated by drugs such as benzodiazepines and barbiturates, is expected to be downregulated after prolonged or repetitive drug exposure. This adaptive response may involve, although not necessarily, reduction in the mRNA levels of the GABAAR subunits. Therefore, comparison of changes in mRNA expression patterns in response to different drugs may offer clues to molecular targets involved in drug actions and drug tolerance/dependence.

The effects of benzodiazepines and barbiturates on the mRNA expression of different GABAAR subunits have been investigated over the last few decades [1217]. Studies of the transcriptional responses of the NMDAR genes to diazepam [1820] and flurazepam [21] are also available. Unfortunately, such studies have led to conflicting results. For example, diazepam was found to downregulate GABAAR α1 mRNA expression in the cerebral cortex by several research groups [1214], but others [17, 22] found diazepam ineffective. Similar situations were found for other GABAAR subunits [23]. Diazepam was found to increase the cortical mRNA contents of NR1 and NR2B subunits [1820]. However, another study showed decreases in hippocampal NR2B mRNA and protein after chronic flurazepam treatment [21]. Compared with the GABAAR and NMDAR, the effects of barbiturates, benzodiazepines and ketamine on the mRNA expression of other types of ion channels are largely unknown.

Ion channels are complex proteins forming ion-per-meable pathways through biological membranes. The ion channels tested in this study include the P-domain channels [24], the pentameric ligand-gated ion channels (pLGICs), the ENaC/Deg ion channels, the ATP-gated ion channels (P2X receptors), the calcium release-activated calcium (CRAC) channels, the inositol 1, 4, 5-trisphosphate receptors (IP3Rs) and ryanodine receptors (RyRs), and the chloride ion channels. The only known ion channel genes not included in this study were the invertebrate ionotropic receptors (IRs), a variant subfamily of iGluRs [25, 26].

The P-domain channels have a common pore architecture composed of four homologous pore-domains contributed by one, two or four subunits arranged in four-fold symmetry. This group of channels include voltage-gated potassium (Kv) channels, cyclic nucleotide-gated (CNG) channels, hyperpolarization-activated CNG (HCN) channels, voltage-gated calcium (CaV) channels, voltage-gated sodium (NaV) channels, sodium-leak channels (NALCN), two-pore channels (TPCs), transient receptor potential (TRP) channels, and glutamate-gated ion channels (or ionotropic glutamate receptors, iGluRs).

Most pLGICs are gated by extracellular ligands, and include the nicotinic acetylcholine receptors (nAChRs), 5-hydroxytryptamine type 3 receptors (5-HT3Rs) GABAARs and glycine receptors (GlyRs). The mammalian pLGIC superfamily also includes zinc-activated ion channels (ZACNs), the invertebrate pLGIC superfamily also includes the glutamate-gated chloride (GluCl) channels, histamine-gated chloride (HisCl) channels, and pH-sensitive chloride (pHCl) channels.

Both ENaC/Deg channels and P2X receptors are trimeric and share similar transmembrane topology. However, they are not homologous in amino acid sequences. The CRAC channels are hexameric plasma membrane proteins mediating the entry of extracellular Ca2+ when the intracellular Ca2+ stores are depleted. The IP3Rs and RyRs, on the other hand, are intracellular membrane proteins important to intracellular Ca2+ signaling. The chloride channels are anion permeable protein complexes (excluding the GABAAR and GlyR) and are less-well understood.

The aim of our study is to measure the effects of midazolam, pentobarbital and ketamine on the ion channel mRNA expression in the water flea Daphnia pulex, a freshwater crustacean with great potential for biomedical research (


Daphniacultures and treatments

The method for culturing D. pulex was the same as previously described [27]. Briefly, daphnids from a single clone were kept in a temperature- and photocycle-controlled tank (20 ± 1°C, 16:8 h light-dark cycle) and fed daily with a mixture of Saccharomycetes and Spirulina. One-fifth of the medium was renewed every second day. Midazolam, pentobarbital sodium and ketamine stock solutions of 10, 100 and 100 mM, respectively, were made in distilled water, and diluted to final concentrations before use. Since the constant jumping movements of Daphnia may complicate experimental results if the animals are immobilized by drugs while the animals in the control group are free to swim, drugs of subanesthetic concentrations based on the EC50 values (Figure 1) were used to investigate the transcriptional responses of ion channel genes. For this purpose, daphnids were exposed to midazolam (100 μM), pentobarbital (200 μM), and ketamine (100 μM) for 4 hours, and killed by crushing immediately for mRNA quantification by RT-qPCR. Protocols were approved by the Scientific and ethics committee of Sichuan university and adhered to the international guidelines for animal care.

Figure 1

Daphnia immobility dose responses of midazolam and midazolam + flumazenil (Panel A), pentobarbital and ketamine (Panel B). 20 daphnids were treated with midazolam, midazolam + flumazenil, pentobarbital and ketamine at series of concentrations. The numbers of immobilized daphnids were observed 4 hours after each treatment. The immobile ratio were presented as means ± SE (n = 3).

Sequence analysis and phylogenetic inference

The Daphnia genome (, NCBI and Uniprot protein database were searched for the Daphnia ion channel gene models based on known human or Drosophila genes. Sequences were analyzed with the maximum-likelihood (ML) method using MEGA 5.10 program [28]. The sequences used are designated in succession by the abbreviation of the species (Hs for Homo sapiens, Dm for Drosophila melanogaster, Ce for C. elegans, Am for Apis mellifera, and Dpul for Daphnia pulex) and the gene name.

RNA extraction, reverse transcription and polymerase chain reaction (PCR)

Total RNA extraction, reverse transcription and PCR were performed as previously described [27]. Briefly, RNA was harvested from 50 crushed daphnids for phase separation, precipitation, and quantification. cDNA was generated using PrimeScriptTM RT reagent Kit DRR037A (TaKaRa), and amplified first by regular PCR to screen primers, which were designed using Primer 3 software [29] based on the scaffold sequences ( Successful primers (Additional file 1: Table S1) were then used for qPCR to quantify changes in gene expression after drug treatment. qPCR was carried out on an iQ5 system (Bio-Rad) using SYBR Premix Ex TaqTM II KIT DRR081A (TaKaRa). Each reaction was run in triplicate and contained 2 μl of cDNA template along with 0.8 μM primers in a total volume of 20 μl. Cycling parameters were 95°C for 30 s to active the DNA polymerase, then 40 cycles of 95°C for 5 s, 55°C for 30 s and 72°C for 30 s. Melting-curves were performed to verify only a single product without primer-dimers. Data were normalized against a house-keeping reference gene β-actin, and were analyzed using the 2-ΔΔCT method [30].

Statistical analysis

All data were presented as means ± SE. qPCR data represent the average of 5 replicate experiments; all results were normalized to β-actin, an internal control, and then to control group. Differences in relative expression of genes were assessed using paired-t test (n = 5). Statistical significance was set at a level of P < 0.05 (*) and P < 0.01 (**).


Phylogenetic analysis

More than 120 gene models of ion channel proteins have been predicted in the Daphnia genome and 108 of them were successfully amplified by RT-PCR (Additional file 2: Table S2). For comparison purposes, these proteins are classified into different categories based on sequence homology molecular structure, and ion selectivity in the case of chloride channels.

The four-fold symmetric P-domain ion channels

A total of 53 P-domain channel genes are detected at transcript level in Daphnia (Additional file 2: Table S2), and they are classified into four superfamilies: Kv/CNG, CaV/NaV, TRP and iGluR. The Kv/CGN members are further classified into three families based on subunit transmembrane topology: the two transmembrane-helix (2TM) family (Additional file 3: Figure S1), the 6TM family (Additional file 4: Figure S2), and the 4TM (K2P) family (Additional file 5: Figure S3). The 4TM potassium genes resulted from duplication of 2TM genes during the evolution, thus the name K2P, while the 6TM domain consists of a 2TM domain and a voltage sensor domain. The CNG and HCN channels are homologous to the 6TM Kv channels (Additional file 4: Figure S2). The CaV or NaV channel contains a single principal subunit with four 6TM domains. The NALCN and TPC channels also belong to the NaV/CaV superfamily. The NALCN is a voltage-independent, TTX-insensitive, and nonselective cation channel underlining the background Na+ leak current [31]. TPC channels are intracellular ion channels mediating the second messenger NAADP-regulated Ca2+ release. The Daphnia NaV/CaV superfamily contains 7 members: 3 CaV, 2 NaV, 1 NALCN and 1 TPC genes (Additional file 6: Figure S4).

TRP channels are highly diverse in function, structure and distribution, with 28 mammalian TRP genes classified into six subfamilies: TRPC, TRPV, TRPM, TRPML, TRPA and TRPP. In addition, the invertebrates have a group known as TRPN, which is also found in zebrafish. The Daphnia genome contains 13 TRP subunit genes (Additional file 7: Figure S5).

The iGluRs are tetrameric and can be divided into three subfamilies based on pharmacology and homology: AMPA, NMDA and kainate (KA) receptors. The Daphnia genome contains eight iGluR homologs: Dpul_Glu-RI, Dpul_Nmdar1-3, Dpul_KaiR1-4 (Additional file 8: Figure S6). A variant iGluR subfamily, the ionotropic receptors (IRs), was not included in this study.

The pentameric ligand-gated ion channels

The Daphnia genome contains 20 pLGIC genes: 12 nAChR (Additional file 9: Figure S7), 5 GABAAR (Grd, Rdl, RdlL, Lcch3 and CG8916), 1 GluCl, and 2 HisCl genes (Additional file 10: Figure S8). There is no RdlL counterpart in Drosophila. The classification of the Daphnia pLGICs into nAChR and GABAAR groups is simply based on homology.

The ENaC/Deg channels and P2X receptors

The ENaC/Deg channels and P2X receptors are distinct classes of trimeric protein complexes. The Daphnia ENaC/Deg homologs are quite diverse, with 14 members detected at the transcript level (Additional file 11: Figure S9). Two Daphnia P2X genes are detected at the transcript level (Dpul_P2XL1 and Dpul_P2XL2 (Additional file 12: Figure S10).

The CRAC channels, IP3Rs, RyRs and the chloride channels

The Daphnia genome contains one CRAC channel gene, one IP3R gene, and one RyR gene (Additional file 12: Figure S10). There are at least five distinct classes of Cl channels, including the ClC channels, the CLIC proteins, bestrophin, the tweety chloride channels, and anoctamin/TMEM16. The Daphnia genome contains 7 ClC genes, 2 CLIC genes, 4 bestrophin genes, 1 tty gene and 3 TMEM16 genes (Additional file 13: Figure S11).

Dose-dependent immobility upon midazolam, pentobarbital, and ketamine treatment

In order to determine the subanesthetic concentrations for chronic treatment in our mRNA assay, daphnids were exposed to aquarium water containing a series of concentrations of midazolam, pentobarbital, or ketamine for four hours (Figure 1). Midazolam, pentobarbital and ketamine immobilized the daphnids at EC50 values of 0.65, 0.92 and 0.84 mM, respectively. Flumazenil (10 μM), a competitive antagonist of benzodiazepines and used to reverse the actions of benzodiazepines in clinical settings, shifted the dose-response curve of midazolam to the right (EC50 = 0.86 mM , Figure 1A), but had no effects of its own on the immobility of the daphnids up to 100 μM. Interestingly, ketamine at concentrations of >200 μM produced a consistent circling behavioral phenotype, mimicking the core behavior aspects of rodents and fish administered with ketamine [32, 33]. This aberrant behavior was completely absent in the control daphnids. Ketamine-induced immobility started at higher concentration (~400 μM). Based on these dose-dependent responses, we exposed daphnias to 100-μM midazolam, 200-μM pentobarbital and 100-μM ketamine for 4 hours for our mRNA assays. At these concentrations, even longer treatment (10 hours) did not result in death of the animals. The use of subanesthetic concentrations instead of anesthetic ones were due to the fact that daphnids jump constantly and this behavior may complicate the experimental results if the animals were immobilized by drugs while the animals in the drug-free (control) group were free to swim.

Effects of midazolam on the transcription of daphniaion channel genes

Out of the 108 genes tested, the transcription of 6 genes was upregulated by midazolam treatment, including the 6TM Kv channel shal, the CNG channel cngl, the nAChR subunit ACHA4, the GABAAR RdlL, DEG/ENaC-11 and P2XL1 (Figure 2). Meanwhile, midazolam downregulated transcription of 13 genes, including the 6TM Kv channels shaker and shawl1, the CNG channel CNGA1, the TRP channels nan and iav, the NMDA receptor Nmdar1, the GABAAR Rdl and GRD, the HisCl hclB, the CLC Clc-c1, Clc-c2, DEG/ENaC-10 and IP3R (Figure 2A). Flumazenil (10 μM) abolished the effects of midazolam on the mRNA expression of Rdl, Shawl1 and Shaker (Figure 2B). Flumazenil alone, however, did not affect the Rdl, Shawl1 and Shaker mRNA expression. In addition, we found that flumazenil (10 μM) alone downregulated Nmdar1 and IP3R transcription (Figure 2C).

Figure 2

Effects of midazolam and flumazenil on the expression of Daphnia ion channel genes. (A). Daphnids were exposed to aquarium water containing 100 μM midazolam for 4 hours and killed for mRNA quantification. Out of 108 ion channel genes, 19 genes are affected by midazolam. (B). Daphnids were co-treated with flumazenil (10 μM) along with midazolam (100 μM) for 4 hours and killed for mRNA quantification. Flumazenil abolishes the effects of midazolam on the mRNA expression of Rdl, Shawl1 and Shaker. (C) Daphnids were exposed to 10 μM flumazenil alone for 4 hours before mRNA quantification. The Nmdar1 and IP3R transcripts were significantly reduced by flumazenil. Statistical significance was set at a level of P < 0.05 (*) and P < 0.01 (**).

Effects of pentobarbital on the transcription of daphniaion channel genes

Pentobarbital (200 μM) downregulated the expression of the 6TM Kv channel genes Shawl1 and Shawl2, the CaV channel gene Ca-alpha1D, the NMDAR Nmdar1, the GABAARs Rdl and GRD, the CLC Clc-c2, while the KAR gene KaiR1 and KaiR4, the AMPAR gene Glu-RI and the nAChR ACHA4 were upregulated by pentobarbital (Figure 3).

Figure 3

Effects of pentobarbital on the expression of the Daphnia ion channel genes. Daphnids were exposed to 200 μM pentobarbital for 4 hours and killed for mRNA quantification. Seven genes are downregulated and four genes are upregulated at transcriptional level by pentobarbital. Statistical significance was set at a level of P < 0.05 (*) and P < 0.01 (**).

Effects of ketamine on the transcription of daphniaion channel genes

Ketamine (100 μM) downregulated transcription of 11 genes, including the 2TM Kir channels Ir and Irk2, the K2P channels Task6, Ork1 and TRESK, the TRP channel Trpgamma, the iGluRs Glu-RI and KaiR1, the CLC CLC-c1 and CLC-c2, and tty (Figure 4).

Figure 4

Effects of ketamine on the expression of Daphnia ion channel genes. Daphnids were exposed to 100 μM ketamine for 4 hours and killed for mRNA quantification. Among 108 ion channel genes, ketamine affects 11 genes at transcriptional level. Surprisingly, all of the 11 genes are downregulated by ketamine. Statistical significance was set at a level of P < 0.05 (*) and P < 0.01 (**).


Invertebrate model organisms, such as C. elegans[34] and Drosophila[35], are desirable models for the studies of the mechanisms of drug tolerance/dependence. These simple animals share many complex traits with mammals, even at behavioral levels [36]. Furthermore, the major groups of ion channels are conserved, simplifying the dissection of core molecular machinery responsible for drug addiction.

The choice of midazolam concentration (100 μM) for the mRNA assays was based on the dose-response curve (Figure 1). Compared to clinic concentrations, which are around 1 μM [37], the high midazolam concentration may result in nonspecific effects of the drug and consequently in changes in gene expression that are not related to midazolam action on GABAA receptors. Nonetheless, the down-regulation of the Rdl mRNA by midazolam appears to be directly mediated by GABAA receptors, since this effect is eliminated by flumazenil (10 μM), a competitive antagonist of benzodiazepines. In addition, the direction of Rdl mRNA regulation is consistent with the prediction by the homeostasis theory of drug tolerance [11]. Despite the role of mRNA regulation is not explicitly implied in the theory, recent studies suggest that protein and mRNA expression levels are correlated [3840]. Interestingly, flumazenil also eliminates midazolam-induced reduction in the mRNA levels of Shawl1 and Shaker (Figure 2B), suggesting a possibility of crosstalk between GABAAR and potassium channel signalings.

The short time course of drug treatment in the current study may undermine the relevance of the results to questions of addiction and drug dependence. For the development of tolerance and dependence in mammals, chronic treatments lasting weeks to months are required. In Drosophila, the lifespan of which is around 30 days under common culture conditions, the time courses of chronic treatments are much shorter, but vary greatly depending on experiment scenarios and the purposes of the study. For example, induction of neurodegeneration mimicking Parkinson's disease requires days of treatments [41, 42]. Alcohol addiction in flies, on the other hand, develops more rapidly. Two types of addiction have been identified [35]: rapid tolerance can be induced by a single brief (less than 60 min) exposure to ethanol, while chronic tolerance requires prolonged (~24 h) ones. Chronic alcohol tolerance depends on protein synthesis, while molecular events downstream of protein synthesis contribute to acute tolerance. However, in general the emerging picture of the regulation of mRNA and protein expression is perhaps more complex than initially thought [43]. Accumulating evidence shows that mRNA expression correlate best with delayed protein expression. For example, the mRNA abundance changes in yeasts occur in a time window of 20 to 240 min after rapamycin treatment, while the protein abundance changes mostly occur at 4 and 6 h of the treatment [44]. This temporal correlation between mRNA and protein expression is not specific to rapamycin, since oxidative stress also induces rapid transcriptional changes which peak at 60 min after the treatment and decay quickly, while the protein expression response is much slower [38]. Nothing is known about the temporal relationship between mRNA and protein abundances in the development of chronic drug tolerance. However, it is highly possible in our case that changes in ion channel mRNA levels may lead to delayed alteration in protein abundance and chronic response to the drug treatments.

Our results suggest that Daphnia are feasible model animals for the investigation of the role of GABAA receptors in drug addiction, including addiction to benzodiazepines. Many early studies suggest that insect GABAA receptors were relatively insensitive to benzodiazepines. These studies involved the use of recombinantly expressed GABAA receptors, such as the homo-oligomeric Rdl receptor [45, 46] and the Grd/LCCH3 receptor [47]. There is, however, evidence that native insect GABAA receptors are modulated by some benzodiazepines. For example, Lees et al. (1987) [48] showed that the 7-nitro benzodiazepine flunitrazepam enhanced the amplitude of GABA-induced currents by up to 70%. More recently, Buckingham et al. (2009) [49] observed that GABA-induced currents in acutely dissociated insect motor neurons were enhanced by both Ro5-4864 and diazepam, whereas clonazepam was ineffective. The pharmacological differences between recombinant and native GABAA receptors suggest more complex subunit combinations in native GABAA receptors. The involvement of GABAA receptors in benzodiazepine dependence appears to be evolutionarily conserved. It is observed in flatworms, the simplest bilaterian animals, that flumazenil (10 μM) antagonizes the abstinence-induced withdrawal from midazolam (10 μM), but has no effect of its own on the behavior of the animals [50].

In Daphnia, ketamine at lower concentrations produces increased movements in the form of circling swimming, while immobility starts at higher concentrations. The dissociation of hyperactivity from immobility is interesting and raises the possibility that the molecular targets responsible for ketamine-induced hyperactivity are different from those for immobility. Currently, we found that ketamine-induced changes in the mRNA expression profile are quite different from those seen with midazolam and pentobarbital treatments. It is worth noting that most of the Daphnia genes affected by ketamine are from the P-domain channel group, including the 2TM (Ir, Irk2), 4TM (Task6, ork1, TRESK), TRP (Trp-gamma), iGluR (Glu-R1, KaiR1) families. A P-domain is basically the 2TM domain and composed of two transmembrane helices connected by a loop region referred to as the P-loop. Four P-loops in a ion channel form the pore selectivity filter. It is believed that the open-channel blockers of NMDARs, such as MK-801 and ketamine, bind to the P-loop amino acid residues. It is therefore highly possible that ketamine may interact with the other P-domain channels in a similar way. The K2P channel, ork1, has been implicated in the regulation of cardiac automatic activity in Drosophila[51], as well as in mammals [52]. In mammals, the kainate receptors are implicated in chronic pain regulation [53]. The Daphnia KaiR1 gene is homologous to the mammalian KARs and was downregulated at the transcript level. Unfortunately, data regarding the physiology of KaiR1 is unavailable.

Compared with the P-domain channels and pLGICs, other ion channels are less well understood. Nonetheless, recent studies suggest that the CLC chloride channels play a role in synaptic transmission and plasticity [54, 55], as well as in neuronal excitability [56]. We found in Daphnia that the Clc-c2 mRNA expression was downregulated by midazolam, pentobarbital and ketamine, while the CLC-c1 mRNA expression was also downregulated by midazolam and ketamine, suggesting a new type of molecular targets possibly involved in drug addiction.


The major groups of ion channels are highly conserved across the animal kingdom. In addition, simple animals, such as C. elegans and Drosophila, share many complex traits with mammals, even at behavioral levels [36]. Therefore, the use of invertebrate model organisms may greatly simplify the dissection of core molecular machinery responsible for drug actions and addiction. In Daphnia, midazolam, pentobarbital and ketamine cause distinct mRNA expression profiles and this observation provide insights into potential novel molecular targets involved in drug actions and addiction.


  1. 1.

    Miller LG, Roy RB, Weill CL: Chronic clonazepam administration decreases gamma-aminobutyric acid A receptor function in cultured cortical neurons. Mol Pharmacol. 1989, 36 (5): 796-802.

    CAS  PubMed  Google Scholar 

  2. 2.

    Loscher W, Rogawski MA: How theories evolved concerning the mechanism of action of barbiturates. Epilepsia. 2012, 53 (Suppl 8): 12-25.

    Article  PubMed  Google Scholar 

  3. 3.

    Costa E, Auta J, Grayson DR, Matsumoto K, Pappas GD, Zhang X, Guidotti A: GABAA receptors and benzodiazepines: a role for dendritic resident subunit mRNAs. Neuropharmacology. 2002, 43 (6): 925-937. 10.1016/S0028-3908(02)00199-5.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Morgan CJ, Curran HV: Ketamine use: a review. Addiction. 2012, 107 (1): 27-38. 10.1111/j.1360-0443.2011.03576.x.

    Article  PubMed  Google Scholar 

  5. 5.

    Kelland MD, Soltis RP, Boldry RC, Walters JR: Behavioral and electrophysiological comparison of ketamine with dizocilpine in the rat. Physiol Behav. 1993, 54 (3): 547-554. 10.1016/0031-9384(93)90248-E.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Irifune M, Katayama S, Takarada T, Shimizu Y, Endo C, Takata T, Morita K, Dohi T, Sato T, Kawahara M: MK-801 enhances gabaculine-induced loss of the righting reflex in mice, but not immobility. Can J Anaesth. 2007, 54 (12): 998-1005. 10.1007/BF03016634.

    Article  PubMed  Google Scholar 

  7. 7.

    Petrenko AB, Yamakura T, Fujiwara N, Askalany AR, Baba H, Sakimura K: Reduced sensitivity to ketamine and pentobarbital in mice lacking the N-methyl-D-aspartate receptor GluRepsilon1 subunit. Anesth Analg. 2004, 99 (4): 1136-1140. 10.1213/01.ANE.0000131729.54986.30. table of contents

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Petrenko AB, Yamakura T, Askalany AR, Kohno T, Sakimura K, Baba H: Effects of ketamine on acute somatic nociception in wild-type and N-methyl-D-aspartate (NMDA) receptor epsilon1 subunit knockout mice. Neuropharmacology. 2006, 50 (6): 741-747. 10.1016/j.neuropharm.2005.11.019.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Sato Y, Kobayashi E, Murayama T, Mishina M, Seo N: Effect of N-methyl-D-aspartate receptor epsilon1 subunit gene disruption of the action of general anesthetic drugs in mice. Anesthesiology. 2005, 102 (3): 557-561. 10.1097/00000542-200503000-00013.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Hevers W, Hadley SH, Luddens H, Amin J: Ketamine, but not phencyclidine, selectively modulates cerebellar GABA(A) receptors containing alpha6 and delta subunits. J Neurosci. 2008, 28 (20): 5383-5393. 10.1523/JNEUROSCI.5443-07.2008.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Martin WR: XVI. A homeostatic and redundancy theory of tolerance to and dependence on narcotic analgesics. Res Publ Assoc Res Nerv Ment Dis. 1968, 46: 206-225.

    CAS  PubMed  Google Scholar 

  12. 12.

    Heninger C, Saito N, Tallman JF, Garrett KM, Vitek MP, Duman RS, Gallager DW: Effects of continuous diazepam administration on GABAA subunit mRNA in rat brain. J Mol Neurosci. 1990, 2 (2): 101-107. 10.1007/BF02876917.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Impagnatiello F, Pesold C, Longone P, Caruncho H, Fritschy JM, Costa E, Guidotti A: Modifications of gamma-aminobutyric acidA receptor subunit expression in rat neocortex during tolerance to diazepam. Mol Pharmacol. 1996, 49 (5): 822-831.

    CAS  PubMed  Google Scholar 

  14. 14.

    Longone P, Impagnatiello F, Guidotti A, Costa E: Reversible modification of GABAA receptor subunit mRNA expression during tolerance to diazepam-induced cognition dysfunction. Neuropharmacology. 1996, 35 (9–10): 1465-1473.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    Pesold C, Caruncho HJ, Impagnatiello F, Berg MJ, Fritschy JM, Guidotti A, Costa E: Tolerance to diazepam and changes in GABA(A) receptor subunit expression in rat neocortical areas. Neuroscience. 1997, 79 (2): 477-487. 10.1016/S0306-4522(96)00609-4.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Chen S, Huang X, Zeng XJ, Sieghart W, Tietz EI: Benzodiazepine-mediated regulation of alpha1, alpha2, beta1-3 and gamma2 GABA(A) receptor subunit proteins in the rat brain hippocampus and cortex. Neuroscience. 1999, 93 (1): 33-44. 10.1016/S0306-4522(99)00118-9.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Wu Y, Rosenberg HC, Chiu TH, Zhao TJ: Subunit- and brain region-specific reduction of GABAA receptor subunit mRNAs during chronic treatment of rats with diazepam. J Mol Neurosci. 1994, 5 (2): 105-120. 10.1007/BF02736752.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Tsuda M, Chiba Y, Suzuki T, Misawa M: Upregulation of NMDA receptor subunit proteins in the cerebral cortex during diazepam withdrawal. Eur J Pharmacol. 1998, 341 (2–3): R1-R2.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Perez MF, Salmiron R, Ramirez OA: NMDA-NR1 and -NR2B subunits mRNA expression in the hippocampus of rats tolerant to Diazepam. Behav Brain Res. 2003, 144 (1–2): 119-124.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Almiron RS, Perez MF, Ramirez OA: MK-801 prevents the increased NMDA-NR1 and NR2B subunits mRNA expression observed in the hippocampus of rats tolerant to diazepam. Brain Res. 2004, 1008 (1): 54-60. 10.1016/j.brainres.2004.01.080.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Van Sickle BJ, Cox AS, Schak K, Greenfield LJ, Tietz EI: Chronic benzodiazepine administration alters hippocampal CA1 neuron excitability: NMDA receptor function and expression(1). Neuropharmacology. 2002, 43 (4): 595-606. 10.1016/S0028-3908(02)00152-1.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Holt RA, Bateson AN, Martin IL: Chronic treatment with diazepam or abecarnil differently affects the expression of GABAA receptor subunit mRNAs in the rat cortex. Neuropharmacology. 1996, 35 (9–10): 1457-1463.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Uusi-Oukari M, Korpi ER: Regulation of GABA(A) receptor subunit expression by pharmacological agents. Pharmacol Rev. 2010, 62 (1): 97-135. 10.1124/pr.109.002063.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Zhorov BS, Tikhonov DB: Potassium, sodium, calcium and glutamate-gated channels: pore architecture and ligand action. J Neurochem. 2004, 88 (4): 782-799. 10.1111/j.1471-4159.2004.02261.x.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Benton R, Vannice KS, Gomez-Diaz C, Vosshall LB: Variant ionotropic glutamate receptors as chemosensory receptors in Drosophila. Cell. 2009, 136 (1): 149-162. 10.1016/j.cell.2008.12.001.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Croset V, Rytz R, Cummins SF, Budd A, Brawand D, Kaessmann H, Gibson TJ, Benton R: Ancient protostome origin of chemosensory ionotropic glutamate receptors and the evolution of insect taste and olfaction. PLoS Genet. 2010, 6 (8): e1001064-10.1371/journal.pgen.1001064.

    Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Sheng B, Liu J, Li GH: Metformin preconditioning protects Daphnia pulex from lethal hypoxic insult involving AMPK, HIF and mTOR signaling. Comp Biochem Physiol B Biochem Mol Biol. 2012, 163 (1): 51-58. 10.1016/j.cbpb.2012.04.009.

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28 (10): 2731-2739. 10.1093/molbev/msr121.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Meth Mol Biol. 2000, 132: 365-386.

    CAS  Google Scholar 

  30. 30.

    Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    Lu B, Su Y, Das S, Liu J, Xia J, Ren D: The neuronal channel NALCN contributes resting sodium permeability and is required for normal respiratory rhythm. Cell. 2007, 129 (2): 371-383. 10.1016/j.cell.2007.02.041.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Myslobodsky MS, Ackermann RF, Golovchinsky V, Engel J: Ketamine-induced rotation: interaction with GABA-transaminase inhibitors and picrotoxin. Pharmacol Biochem Behav. 1979, 11 (5): 483-486. 10.1016/0091-3057(79)90029-7.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Zakhary SM, Ayubcha D, Ansari F, Kamran K, Karim M, Leheste JR, Horowitz JM, Torres G: A behavioral and molecular analysis of ketamine in zebrafish. Synapse. 2011, 65 (2): 160-167. 10.1002/syn.20830.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  34. 34.

    Schafer WR: Addiction research in a simple animal model: the nematode Caenorhabditis elegans. Neuropharmacology. 2004, 47 (Suppl 1): 123-131.

    CAS  Article  PubMed  Google Scholar 

  35. 35.

    Kaun KR, Devineni AV, Heberlein U: Drosophila melanogaster as a model to study drug addiction. Hum Genet. 2012, 131 (6): 959-975. 10.1007/s00439-012-1146-6.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  36. 36.

    O'Kane CJ: Drosophila as a model organism for the study of neuropsychiatric disorders. Curr Top Behav Neurosci. 2011, 7: 37-60. 10.1007/7854_2010_110.

    Article  PubMed  Google Scholar 

  37. 37.

    Crevat-Pisano P, Dragna S, Granthil C, Coassolo P, Cano JP, Francois G: Plasma concentrations and pharmacokinetics of midazolam during anaesthesia. J Pharm Pharmacol. 1986, 38 (8): 578-582. 10.1111/j.2042-7158.1986.tb03084.x.

    CAS  Article  PubMed  Google Scholar 

  38. 38.

    Vogel C, Abreu Rde S, Ko D, Le SY, Shapiro BA, Burns SC, Sandhu D, Boutz DR, Marcotte EM, Penalva LO: Sequence signatures and mRNA concentration can explain two-thirds of protein abundance variation in a human cell line. Mol Syst Biol. 2010, 6: 400-

    Article  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Schwanhausser B, Busse D, Li N, Dittmar G, Schuchhardt J, Wolf J, Chen W, Selbach M: Global quantification of mammalian gene expression control. Nature. 2011, 473 (7347): 337-342. 10.1038/nature10098.

    Article  PubMed  Google Scholar 

  40. 40.

    Lackner DH, Schmidt MW, Wu S, Wolf DA, Bahler J: Regulation of transcriptome, translation, and proteome in response to environmental stress in fission yeast. Genome Biol. 2012, 13 (4): R25-10.1186/gb-2012-13-4-r25.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  41. 41.

    Coulom H, Birman S: Chronic exposure to rotenone models sporadic Parkinson's disease in Drosophila melanogaster. J Neurosci. 2004, 24 (48): 10993-10998. 10.1523/JNEUROSCI.2993-04.2004.

    CAS  Article  PubMed  Google Scholar 

  42. 42.

    Bonilla-Ramirez L, Jimenez-Del-Rio M, Velez-Pardo C: Acute and chronic metal exposure impairs locomotion activity in Drosophila melanogaster: a model to study Parkinsonism. Biometals. 2011, 24 (6): 1045-1057. 10.1007/s10534-011-9463-0.

    CAS  Article  PubMed  Google Scholar 

  43. 43.

    Gygi SP, Rochon Y, Franza BR, Aebersold R: Correlation between protein and mRNA abundance in yeast. Mol Cell Biol. 1999, 19 (3): 1720-1730.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  44. 44.

    Fournier ML, Paulson A, Pavelka N, Mosley AL, Gaudenz K, Bradford WD, Glynn E, Li H, Sardiu ME, Fleharty B, et al: Delayed correlation of mRNA and protein expression in rapamycin-treated cells and a role for Ggc1 in cellular sensitivity to rapamycin. Mol Cell Proteomics. 2010, 9 (2): 271-284. 10.1074/mcp.M900415-MCP200.

    CAS  Article  PubMed  Google Scholar 

  45. 45.

    Millar NS, Buckingham SD, Sattelle DB: Stable expression of a functional homo-oligomeric Drosophila GABA receptor in a Drosophila cell line. Proc Biol Sci. 1994, 258 (1353): 307-314. 10.1098/rspb.1994.0178.

    CAS  Article  PubMed  Google Scholar 

  46. 46.

    Hosie AM, Sattelle DB: Allosteric modulation of an expressed homo-oligomeric GABA-gated chloride channel of Drosophila melanogaster. Br J Pharmacol. 1996, 117 (6): 1229-1237. 10.1111/j.1476-5381.1996.tb16720.x.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  47. 47.

    Gisselmann G, Plonka J, Pusch H, Hatt H: Drosophila melanogaster GRD and LCCH3 subunits form heteromultimeric GABA-gated cation channels. Br J Pharmacol. 2004, 142 (3): 409-413. 10.1038/sj.bjp.0705818.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  48. 48.

    Lees G, Beadle DJ, Neumann R, Benson JA: Responses to GABA by isolated insect neuronal somata: pharmacology and modulation by a benzodiazepine and a barbiturate. Brain Res. 1987, 401 (2): 267-278. 10.1016/0006-8993(87)91411-9.

    CAS  Article  PubMed  Google Scholar 

  49. 49.

    Buckingham SD, Higashino Y, Sattelle DB: Allosteric modulation by benzodiazepines of GABA-gated chloride channels of an identified insect motor neurone. Invert Neurosci. 2009, 9 (2): 85-89. 10.1007/s10158-009-0091-0.

    CAS  Article  PubMed  Google Scholar 

  50. 50.

    Raffa RB, Cavallo F, Capasso A: Flumazenil-sensitive dose-related physical dependence in planarians produced by two benzodiazepine and one non-benzodiazepine benzodiazepine-receptor agonists. Eur J Pharmacol. 2007, 564 (1–3): 88-93.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  51. 51.

    Lalevee N, Monier B, Senatore S, Perrin L, Semeriva M: Control of cardiac rhythm by ORK1, a Drosophila two-pore domain potassium channel. Curr Biol. 2006, 16 (15): 1502-1508. 10.1016/j.cub.2006.05.064.

    CAS  Article  PubMed  Google Scholar 

  52. 52.

    Gurney A, Manoury B: Two-pore potassium channels in the cardiovascular system. Eur Biophys J. 2009, 38 (3): 305-318. 10.1007/s00249-008-0326-8.

    CAS  Article  PubMed  Google Scholar 

  53. 53.

    Bhangoo SK, Swanson GT: Kainate receptor signaling in pain pathways. Mol Pharmacol. 2013, 83 (2): 307-315. 10.1124/mol.112.081398.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  54. 54.

    Riazanski V, Deriy LV, Shevchenko PD, Le B, Gomez EA, Nelson DJ: Presynaptic CLC-3 determines quantal size of inhibitory transmission in the hippocampus. Nat Neurosci. 2011, 14 (4): 487-494. 10.1038/nn.2775.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  55. 55.

    Farmer LM, Le BN, Nelson DJ: CLC-3 chloride channels moderate long-term potentiation at Schaffer collateral-CA1 synapses. J Physiol. 2013, 591 (Pt 4): 1001-1015.

    CAS  Article  PubMed  Google Scholar 

  56. 56.

    Chen TT, Klassen TL, Goldman AM, Marini C, Guerrini R, Noebels JL: Novel brain expression of ClC-1 chloride channels and enrichment of CLCN1 variants in epilepsy. Neurology. 2013, 80 (12): 1078-1085. 10.1212/WNL.0b013e31828868e7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

Pre-publication history

  1. The pre-publication history for this paper can be accessed here:

Download references


This work was supported by a grant to GHL from the National Natural Science Foundation of China (Grant No. 30972861).

Author information



Corresponding author

Correspondence to Guo Hua Li.

Additional information

Competing interests

The authors declare that they have no competing interest.

Authors’ contributions

CD carried out the pharmacological and qPCR assays, performed the statistical analysis and helped draft the manuscript. AH carried out sequence analysis and screened the primers. YN was involved in the pharmacological studies. YZ contributed to the study design and helped draft manuscript. GHL designed the study and drafted the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Table S1: Gene name, Gene ID, Genome Map, Primer, product size and qPCR efficiency. The forward and reverse primers, gene ID, product size and qPCR efficiency for indicated genes were given, based on the scaffold sequences of the Daphnia Genomics Consortium. All primers were designed using Primer 3 (version 0.4.0). The expression and length of each gene were confirmed by reverse transcription PCR. qPCR efficiency of target genes was calculated by Bio-Rad iQ5 software. (DOC 216 KB)

Additional file 2: Table S2: Gene name, synonyms of the gene name and Uniprot accession numbers used in phylogenies. (DOC 580 KB)

Additional file 3: Figure S1: Phylogenetic tree of the 2TM-KV (Kir) channels. The Kir subunit possesses a P-domain composed of two transmembrane helices connected by a P-loop. Four Kir subunits assemble to for a channel. The mammalian Kir family comprises the classical inward-rectifying K+ channels (Kir2), the KATP channels (Kir6), G protein-activated K+ channels (Kir3), and K+-transport channels. The invertebrate Kir genes form distinct clusters. The Daphnia genome contains two Kir genes: Dpul_Ir and Dpul_Irk2. The homolog of the Drosophila Irk3 is absent in Daphnia genome. Abbreviation: Hs, Homo sapiens; Dm, Drosophila melanogaster; Ce, Caenorhabditis elegans; Dpul, Daphnia pulex. (TIFF 2 MB)

Additional file 4: Figure S2: Phylogenetic tree of the 6TM-Kv, CNG and HCN channels. The 6TM-Kv channels are highly diverse in sequence, structure and function. In addition to the P-domain, each 6TM-KV subunit obtains a voltage sensor domain (VSD) composed of four transmembrane helices. The CNG and HCN channels, although not classified as Kv channels, are homologous to the KCNH KV channels. Daphnia and Drosophila share similar gene sets for the 6TM-KV/CNG/HCN channel group, with subtle differences. The homolog of the Drosophila elk is absent in the Daphnia genome. Daphnia have two homologs of the Drosophila Shawl genes, namely Shawl1 and Shawl2, while Drosophila have one. (TIFF 2 MB)

Additional file 5: Figure S3: Phylogenetic tree of the 4TM-KV (K2P) channels. The K2P subunit consists of two P-domains, and two such subunits assemble to form a channel. Daphnia have five K2P members: Dpul_TWIK, Dpul_Task6, Dpul_Task7, Dpul_Ork1 and Dpul_TRESK. The Daphnia homolog of the mammalian TWIK is detected at transcript level, but the Drosophila counterpart is absent. As in Drosophila, the TREK member of the K2P channels is also absent in Daphnia. Abbreviation: Hs, Homo sapiens; Dm, Drosophila melanogaster; Ce, Caenorhabditis elegans; Dpul, Daphnia pulex. (TIFF 2 MB)

Additional file 6: Figure S4: Phylogenetic tree of the CaV, NaV, NALCN and TPC channels. The CaV and NaV channels are thought to arise from the potassium channels during evolution. The CaV/NaV channel contains a single principal (α1) subunit with four 6TM domains. The human genome contains ten α1 genes, which can be further clustered in three groups: CaV1 (L-type), CaV2 (P/Q, N, and R type), and CaV3 (T-type). The Daphnia genome predicts three genes: Dpul_Ca-alpha1D (CaV1), Dpul_Ca-alpha1T (CaV2) and Dpul_cac (CaV3). NALCN is represented by a single gene in human, Daphnia and Drosophila. Unlike the CaV, NaV and NALCN, the TPC subunit contains two P-domains and assembles as a dimer. The transcript of the TPC encoding gene is detected in Daphnia. TPC Gene models are also available for many insects, such as Apis mellifera, but absent in Drosophila. Abbreviation: Hs, Homo sapiens; Dm, Drosophila melanogaster; Am, Apis mellifera; Caenorhabditis elegans; Dpul, Daphnia pulex. (TIFF 2 MB)

Additional file 7: Figure S5: Phylogenetic tree of the TRP channels. Despite of the structure and function diversities, all TRP subunits have a common design for the transmembrane domain, which basically is a 6TM domain. So far 28 mammalian TRP subunit encoding genes are known and belong to six subfamilies known as the TRPC, TRPV, TRPM, TRPML, TRPP and the TRPA, and all these groups have invertebrate orthologs. The Daphnia genome predicts 13 TRP subunit genes and all were detected at transcript level, including a 3 TRPC-like genes (Dpul_TRP, Dpul_TRPL and Dpul_Trpgamma), 2 TRPM-like gene (Dpul_TRPM1, Dpul_TRPM2), 4 TRPA-like genes (Dpul_TRPA5, Dpul_pain, Dpul_pyx1, Dpul_pyx2), 3TRPV genes (Dpul_nan, Dpul_iav, TRPV_nompc), 1TRPML gene (TRPML). TRPP genes are absent in Daphnia. Abbreviation: Hs, Homo sapiens; Dm, Drosophila melanogaster; Ce, Caenorhabditis elegans; Am, Apis mellifera; Dr, Danio rerio; Dpul, Daphnia pulex. (TIFF 2 MB)

Additional file 8: Figure S6: Phylogenetic tree of the iGluRs. The iGluR is a heterotetramer containing four P-domains, each in one subunit. The mammalian iGluRs are divided into three groups: AMPA, NMDA and KA receptors. There are two other classes with respect to sequence similarity known as the delta class (GRID) and the kainate-binding proteins (KBP). The invertebrate iGluRs display considerable variations, especially in the KAR branch. In addition, a GRIN3-like gene is detected in Daphnia at the transcript level, but it is absent in Drosophila. Abbreviation: Hs, Homo sapiens; Dm, Drosophila melanogaster; Xl, Xenopus laevis; Gg, Gallus gallus; Dpul, Daphnia pulex. (TIFF 2 MB)

Additional file 9: Figure S7: Phylogenetic tree of the nAChRs. The nAChRs belong to the Cys-loop receptor superfamily. The nAChRs are highly diverse in both vertebrates and invertebrates. Daphnia genome encodes 12 putative nAChR genes. Abbreviation: Hs, Homo sapiens; Dm, Drosophila melanogaster; Gg, Gallus gallus; Am, Apis mellifera; Dpul, Daphnia pulex. (TIFF 2 MB)

Additional file 10: Figure S8: Phylogenetic tree of non-nAChR cys-loop receptors. In addition to the nAChRs, the Cys-loop receptor superfamily also includes 5-hydroxytryptamine type 3 receptors (5-HT3Rs), zinc-activated ion channels (ZACNs), GABAARs and glycine receptors (GlyRs). The Daphnia genome lacks the counterparts of the Drosophila genes CG11340, CG6927 and CG7589, which form a separate branch. The pHCl gene seen in Drosophila is also absent in Daphnia. (TIFF 4 MB)

Additional file 11: Figure S9: Phylogenetic tree of the Deg/ENaC channels. DEG/ENaC channels have been found in nematodes, insects and vertebrates and implicated in a broad spectrum of cellular functions. Mammalian DEG/ENaC channels fall into two major groups (EnaC and ASIC) with 9 members. The invertebrate DEG/ENaC members are highly diverse. Fourteen Daphnia ENaC/Deg homologs are detected at transcript level. (TIFF 1 MB)

Additional file 12: Figure S10: Phylogenetic trees of the P2X receptor, ORAI proteins, IP3Rs and RyRs. Two P2X gene transcripts are detected in Daphnia, meanwhile, one ORAI, one IP3R and one RyR are detected in Daphnia. (TIFF 2 MB)

Additional file 13: Figure S11: Phylogenetic trees of the CLC, CLIC, BEST, tweety and ANO proteins. Chloride channels are a functionally and structurally diverse group of anion selective channels. Chloride channels are poorly understood and classified into ClC, Clic, Best, Tweety and Ano families. There is no sequence homology between these families. (TIFF 2 MB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Dong, C., Hu, A., Ni, Y. et al. Effects of midazolam, pentobarbital and ketamine on the mRNA expression of ion channels in a model organism Daphnia pulex . BMC Anesthesiol 13, 32 (2013).

Download citation


  • Midazolam
  • Pentobarbital
  • Ketamine
  • mRNA
  • Ion channel
  • Daphnia pulex