Molecular Cloning and Characterization of Novel Glutamate-Gated Chloride Channel Subunits from
Cys-loop ligand-gated ion channels (LGICs) mediate fast ionotropic neurotransmission. They are proven drug targets in nematodes and arthropods, but are poorly characterized in flatworms. In this study, we characterized the anion-selective, non-acetylcholine-gated Cys-loop LGICs from Schistosoma mansoni. Full-length cDNAs were obtained for SmGluCl-1 (Smp_096480), SmGluCl-2 (Smp_015630) and SmGluCl-3 (Smp_104890). A partial cDNA was retrieved for SmGluCl-4 (Smp_099500/Smp_176730). Phylogenetic analyses suggest that SmGluCl-1, SmGluCl-2, SmGluCl-3 and SmGluCl-4 belong to a novel clade of flatworm glutamate-gated chloride channels (GluCl) that includes putative genes from trematodes and cestodes. The flatworm GluCl clade was distinct from the nematode-arthropod and mollusc GluCl clades, and from all GABA receptors. We found no evidence of GABA receptors in S. mansoni. SmGluCl-1, SmGluCl-2 and SmGluCl-3 subunits were characterized by two-electrode voltage clamp (TEVC) in Xenopus oocytes, and shown to encode Cl−-permeable channels gated by glutamate. SmGluCl-2 and SmGluCl-3 produced functional homomers, while SmGluCl-1 formed heteromers with SmGluCl-2. Concentration-response relationships revealed that the sensitivity of SmGluCl receptors to L-glutamate is among the highest reported for GluCl receptors, with EC50 values of 7–26 µM. Chloride selectivity was confirmed by current-voltage (I/V) relationships. SmGluCl receptors are insensitive to 1 µM ivermectin (IVM), indicating that they do not belong to the highly IVM-sensitive GluClα subtype group. SmGluCl receptors are also insensitive to 10 µM meclonazepam, a schistosomicidal benzodiazepine. These results provide the first molecular evidence showing the contribution of GluCl receptors to L-glutamate signaling in S. mansoni, an unprecedented finding in parasitic flatworms. Further work is needed to elucidate the roles of GluCl receptors in schistosomes and to explore their potential as drug targets.
Published in the journal:
Molecular Cloning and Characterization of Novel Glutamate-Gated Chloride Channel Subunits from. PLoS Pathog 9(8): e32767. doi:10.1371/journal.ppat.1003586
Category:
Research Article
doi:
https://doi.org/10.1371/journal.ppat.1003586
Summary
Cys-loop ligand-gated ion channels (LGICs) mediate fast ionotropic neurotransmission. They are proven drug targets in nematodes and arthropods, but are poorly characterized in flatworms. In this study, we characterized the anion-selective, non-acetylcholine-gated Cys-loop LGICs from Schistosoma mansoni. Full-length cDNAs were obtained for SmGluCl-1 (Smp_096480), SmGluCl-2 (Smp_015630) and SmGluCl-3 (Smp_104890). A partial cDNA was retrieved for SmGluCl-4 (Smp_099500/Smp_176730). Phylogenetic analyses suggest that SmGluCl-1, SmGluCl-2, SmGluCl-3 and SmGluCl-4 belong to a novel clade of flatworm glutamate-gated chloride channels (GluCl) that includes putative genes from trematodes and cestodes. The flatworm GluCl clade was distinct from the nematode-arthropod and mollusc GluCl clades, and from all GABA receptors. We found no evidence of GABA receptors in S. mansoni. SmGluCl-1, SmGluCl-2 and SmGluCl-3 subunits were characterized by two-electrode voltage clamp (TEVC) in Xenopus oocytes, and shown to encode Cl−-permeable channels gated by glutamate. SmGluCl-2 and SmGluCl-3 produced functional homomers, while SmGluCl-1 formed heteromers with SmGluCl-2. Concentration-response relationships revealed that the sensitivity of SmGluCl receptors to L-glutamate is among the highest reported for GluCl receptors, with EC50 values of 7–26 µM. Chloride selectivity was confirmed by current-voltage (I/V) relationships. SmGluCl receptors are insensitive to 1 µM ivermectin (IVM), indicating that they do not belong to the highly IVM-sensitive GluClα subtype group. SmGluCl receptors are also insensitive to 10 µM meclonazepam, a schistosomicidal benzodiazepine. These results provide the first molecular evidence showing the contribution of GluCl receptors to L-glutamate signaling in S. mansoni, an unprecedented finding in parasitic flatworms. Further work is needed to elucidate the roles of GluCl receptors in schistosomes and to explore their potential as drug targets.
Introduction
Schistosomiasis, a disease caused by parasitic flatworms in the genus Schistosoma, is one of the most prevalent parasitic diseases in tropical and sub-tropical areas of the world. S. mansoni is responsible for the majority of schistosomiasis infections in sub-Saharan Africa, the Middle East, the Caribbean and South America [1]. Global statistics for 2003 revealed that an estimated 207 million people were infected, of whom >85% live in Africa; 120 million people suffered from clinical disease and 779 million were at risk of infection [2]. Schistosomiasis leads to a chronic, often debilitating disease that impairs growth, development and productivity in infected individuals, and is strongly linked to extreme poverty, particularly in sub-Saharan Africa [3], [4]. Moreover, in Africa alone, 280,000 deaths/year are attributed to the severe complications caused by schistosomiasis [5]. No vaccines are available against schistosome species, and schistosomiasis control relies almost entirely on a single drug, praziquantel. Growing concerns about sub-optimal efficacy of praziquantel and the prospect of drug resistance [6]–[10] highlight the need to identify targets for the discovery of new schistosomicidal drugs.
The nervous system of S. mansoni is an attractive target for development of new therapeutic drugs. It coordinates many functions vital to parasite survival and reproduction, including host attachment and penetration, motor activity and migration, feeding and excretion, pairing, and egg laying. The nervous system is also presumed to control long-distance signal transduction, through yet undefined mechanisms, since schistosomes lack a coelom and a proper circulatory system to support endocrine signaling [11], [12]. The S. mansoni genome predicts a rich diversity of neuroreceptors, including several members of the Cys-loop ligand-gated ion channel (Cys-loop LGIC) superfamily [13], [14]. The Cys-loop LGICs found in S. mansoni or in other schistosome species have not been characterized at the molecular level, and the physiological functions mediated by these receptors have not been elucidated. The only schistosome Cys-loop LGICs cloned to date are three putative nicotinic acetylcholine receptor subunits from S. haematobium, named ShAR1α, ShAR1β and ShAR2β [15], [16].
Cys-loop LGICs are neurotransmitter receptors instrumental in mediating fast ionotropic neurotransmission in vertebrate and invertebrate nervous systems. They are targets for an extensive range of active molecules, including insecticides [17], antiparasitic agents [18], anaesthetics, muscle relaxants and drugs for neurological disorders [19], [20]. Cys-loop LGICs are divided into excitatory and inhibitory receptors, based on permeability to cations or anions, respectively. Different families of Cys-loop LGICs are categorized based on ligand specificity. Vertebrate Cys-loop LGIC families include inhibitory γ-aminobutyric acid (GABA) and glycine (Gly) receptors, as well as excitatory nicotinic acetylcholine (nACh), serotonin (5-HT) and zinc receptors (R) [21]. Invertebrate Cys-loop LGIC families include not only the conventional inhibitory GABA receptors and excitatory nAChRs, but also a broad range of non-nicotinic inhibitory glutamate (GluCl), biogenic amines, ACh and pH-gated (pHCl) receptors, as well as excitatory GABA receptors [22]. All Cys-loop LGICs share a common pentameric structure arranged to form an ion-conducting pore. Channels can be homomeric or heteromeric. Cys-loop LGIC subunits are made up of an amino-terminal extracellular domain (ECD) flanked by a signal peptide (SP), and a carboxy-terminal transmembrane domain (TMD) [21]. The ECD contains the ligand-binding site and a 13 amino acid cysteine-loop signature motif (Cys-loop). The TMD is made up of four membrane spanning regions (M1-M4) and an intracellular domain (ICD) located between M3 and M4. A consensus motif on the cytoplasmic side of M2 determines whether channels are anion- or cation-selective [23], [24]. A second 12–14 amino acid Cys-loop motif located just upstream of M1 is found in a subset of inhibitory LGIC subunits which includes the glycine-, glutamate- and histamine-gated chloride channels (GlyRs, GluCls, HisCls) [22].
A potentially time-saving and cost-effective approach to address the need for new schistosomicidal drugs is to develop existing compounds with some antiparasitic activity. In this regard, the benzodiazepine (BZD) derivative meclonazepam – an anxiolytic BZD that is an allosteric modulator of mammalian Cys-loop GABA receptors [19] – has emerged as an interesting lead candidate. Meclonazepam exhibited potent therapeutic activity against both mature and immature stages of S. mansoni and S. haematobium. [25]. Unfortunately, the drug was later discarded as a lead candidate due to its lack of selectivity; sedative and hypnotic effects were seen in humans at therapeutic doses [26]. There is hope that medicinal chemistry efforts, combined with a better understanding of the molecular target of meclonazepam, could lead to the development of derivatives with greater selectivity towards schistosomes and minimal adverse effects in patients [27], [28].
We hypothesized that the molecular target of meclonazepam in S. mansoni could be a homolog of mammalian GABA receptors. Cys-loop LGICs are poorly characterized in flatworms, including S. mansoni. We cloned three of the four subunits identified by homology in S. mansoni. We characterized these subunits in Xenopus oocytes and demonstrated that they form functional GluCl receptors that are insensitive to ivermectin and meclonazepam. We also show that these GluCl subunits form, together with other putative subunits from trematodes and cestodes, a novel family of flatworm GluCl subunits distinct from their snail, nematode and arthropod counterparts.
Results and Discussion
A novel family of non-AChR-like, inhibitory Cys-loop LGICs in S. mansoni
A protein BLAST analysis of an early version of the S. mansoni genome database (v3.1), using vertebrate and invertebrate Cys-loop GABA receptor subunits as queries, identified five non-AChR-like, inhibitory Cys-loop LGIC subunit gene candidates. Smp_015630 (XM_002572982.1), Smp_096480 (XM_002580489.1), Smp_104890 (XM_002570508.1), Smp_099500 (XM_002570087.1) and Smp_176730 (XM_002570085.1) are predicted open reading frames (ORFs) that encode putative Cys-loop LGIC subunit genes related to GABA, glutamate and glycine receptors, but precise identities were not assigned to them in the first draft genome (v4.0) [13]. In the latest version of the genome (v5.0), Smp_015630 and Smp_096480 were mapped to chromosome W/Z and were tentatively annotated as glycine receptor subunit α1 and β, respectively. Smp_104890 was mapped to chromosome 1, and Smp_099500 and Smp_176730 were mapped to chromosome 6, but subtype predictions were not specified [14]. Interestingly, recent work pointed out that these putative genes possess a glycine residue in loop B of the ECD, a position where small residues are correlated with L-glutamate binding [29], [30]. Moreover, screening of the NCBI database (July 2012) also suggests that S. mansoni putative inhibitory Cys-loop LGIC subunits share closest homology to invertebrate GluCl subunits, rather than GABA or glycine receptors. This may be explained by recent increases in the number of invertebrate Cys-loop LGIC encoding gene sequences available in public databases.
SMART [31] InterproScan [32] and ScanProsite [33] analyses confirmed that, while some of these putative genes do not appear to encode full length subunits, they all possess core features characteristic of Cys-loop LGICs. In addition, four of these putative Cys-loop LGICs contain an additional Cys-loop motif equivalent to those found in GlyRs, GluCls, and HisCls. The 1,536 bp Smp_096480 predicted ORF encodes a full-length Cys-loop LGIC subunit. In contrast, the 1,269 bp Smp_015630 and 1,515 bp Smp_104890 predicted ORFs encode partial subunits lacking a SP and a M4. The 885 bp Smp_099500 predicted ORF corresponds to an ECD lacking a SP, while the 840 bp Smp_176730 predicted ORF encodes a TMD lacking M4 and is flanked at the N-terminus by a 14 amino acid Cys-loop motif corresponding to the second Cys-loop motif mentioned above.
RACE experiments generated full-length Smp_096480 (1,536 bp SmGluCl-1, KC861381; Figure 1A) and Smp_104890 (1,545 bp SmGluCl-3, KC861384) cDNAs, and two full-length Smp_015630 cDNA variants (1,659 bp SmGluCl-2.1, KC861382; and 1,638 bp SmGluCl-2.2, KC861383). Alignment of SmGluCl-2 and SmGluCl-3 cDNAs with their corresponding genomic DNAs shows that both are encoded in a single exon, contrasting with in silico predictions from the genome database (Figure 1B–C) [13]. Repeated attempts only retrieved a partial Smp_099500 cDNA which comprises Smp_176730 (1,986 bp SmGluCl-4, KC861385), but still lacks a SP and M4 (Figure 1D and Figure S1). In addition, the sequence corresponding to the last exon of the SmGluCl-4 cDNA did not correspond to any regions of the first draft genome (v4.0). The 3′RACE experiments that identified Smp_176730 as part of the Smp_099500 gene were further supported by PCR amplifications of adult S. mansoni cDNA using primer pairs combining a Smp_099500-specific sense primer with a Smp_176730-specific antisense primer. Sequencing of the resulting products confirmed that Smp_099500 and Smp_176730 belong to the same locus and encode a single Cys-loop LGIC subunit gene, SmGluCl-4. At odds with this finding, an 83.5 kb gap, unusually large for an intron, is found between the two Cys-loop motifs of the SmGluCl-4 cDNA when its sequence is aligned to genomic DNA. In comparison, the largest intron currently predicted in the S. mansoni genome is 33.8 kb [13]. The publication of chromosome 6 in the latest version of S. mansoni genome provided new data that allowed us to identify a putative 2,598 bp SmGluCl-4 ORF which included a SP and M4 (not shown). However, attempts to amplify the putative 5′end region of this SmGluCl-4 ORF from parasite cDNA were unsuccessful. PCRs using a sense primer targeting the predicted start codon in conjunction with an antisense primer targeting a known region of SmGluCl-4 failed to amplify SmGluCl-4 products. These results suggest that annotation of the region of chromosome 6 in which SmGluCl-4 is located may require further work.
SmGluCl-2 and SmGluCl-3 produce functional homomeric GluCl receptors
We functionally characterized full-length SmGluCl-1, SmGluCl-2, and SmGluCl-3 subunits. We did not pursue characterization of SmGluCl-4, due to the inability to obtain a full-length cDNA. We performed two-electrode voltage clamp (TEVC) experiments on Xenopus oocytes injected with in vitro-transcribed capped RNA (cRNA) encoding the SmGluCl-1, SmGluCl-2.1, SmGluCl-2.2, and SmGluCl-3 subunits. Application of 1 mM L-glutamate produced no detectable response in oocytes injected with SmGluCl-1 cRNA (Figure 2A), suggesting that SmGluCl-1 does not form a functional glutamate-gated homomeric receptor in Xenopus oocytes. In contrast, oocytes injected with SmGluCl-2.1 cRNA or SmGluCl-3 cRNA exhibited robust, rapidly activating and fully reversible currents upon application of 1 mM L-glutamate (Figure 2B–C), indicating that SmGluCl-2.1 and SmGluCl-3 form functional homomeric glutamate receptors. Glutamate-sensitive currents desensitized rapidly and with biphasic kinetics in the continued presence of L-glutamate. The L-glutamate response in oocytes injected with SmGluCl-2.2 cRNA was comparable to the response evoked in oocytes injected with the SmGluCl-2.1 isoform (data not shown). GABA, glycine, ACh, 5-HT, tyramine, histamine and L-aspartate (1 mM) elicited no detectable currents in oocytes injected with any of the subunits, either alone or in combination (data not shown). Control oocytes injected with water and non-injected oocytes showed no response to L-glutamate at concentrations up to 5 mM (data not shown).
Both SmGluCl-2.1 and SmGluCl-3 homomeric receptors were highly sensitive to L-glutamate and responded in a concentration-dependent manner. Oocytes injected with SmGluCl-2.1 cRNA responded to L-glutamate concentrations ≥2 µM; current amplitude saturated at 100 µM L-glutamate (Figure 3A). The effector concentration for half-maximum response (EC50) for L-glutamate on SmGluCl-2.1 receptors was 11.8±0.6 µM with a Hill coefficient of 2.1±0.2 (n = 8, Figure 3C). A Hill coefficient of >2 indicates at least two L-glutamate molecules must bind to activate the channels with positive cooperativity. Concentration-response curves for SmGluCl-2.2 receptors were indistinguishable from those obtained from SmGluCl-2.1 receptors (data not shown). Consequently, we pursued SmGluCl-2 characterization using SmGluCl-2.1. SmGluCl-3 receptors were more sensitive to L-glutamate than SmGluCl-2 receptors (P<0.05, paired-T test). Oocytes injected with SmGluCl-3 cRNA responded to L-glutamate concentrations ≥500 nM; current amplitude again saturated at 100 µM L-glutamate (Figure 3B). The EC50 for L-glutamate on SmGluCl-3 receptors was 6.9±0.5 µM and the Hill coefficient was 1.3±0.1 (n = 8, Figure 3C), showing no evidence of positive cooperativity.
SmGluCl-1 and SmGluCl-2 produce functional heteromeric GluCl receptors
To assess whether the SmGluCl-1 subunit could form a functional glutamate receptor, we co-expressed it with the other subunits. The underlying rationale was that, since a homomeric SmGluCl-1 receptor is not detectable in Xenopus oocytes, any shift in the concentration-response to L-glutamate relative to the homomeric receptors would be indicative of functional heteromeric glutamate receptors containing SmGluCl-1 subunits. We determined the L-glutamate concentration-response relationships of oocytes injected with a 1∶1 ratio of SmGluCl-1 cRNA and either SmGluCl-2.1 cRNA (SmGluCl-1:SmGluCl-2.1), or SmGluCl-3 cRNA (SmGluCl-1:SmGluCl-3). Co-injection of SmGluCl-1 and SmGluCl-2.1 cRNAs decreased sensitivity to L-glutamate compared to oocytes injected with SmGluCl-2.1 cRNA alone (P<0.05, paired-T test), confirming that SmGluCl-1 is a glutamate receptor subunit able to form heteromeric receptors with SmGluCl-2.1. The EC50 for L-glutamate on heteromeric SmGluCl-1:SmGluCl-2.1 receptors increased to 26.3±2.0 µM, while the Hill coefficient decreased to 1.2±0.1 (n = 8, Figure 4A), indicating significant loss of cooperativity in the heteromeric receptors compared to homomeric SmGluCl-2.1 receptors. It should be noted that the co-injection of a 1∶1 ratio of SmGluCl-1 and SmGluCl-2.1 cRNAs likely produced more than one population of heteromeric receptors with distinct subunit stoichiometries, along with a third population of homomeric SmGluCl-2.1 receptors. The concentration-response obtained reflects the contribution of the different populations of receptors present. It is not known whether the SmGluCl-1 and SmGluCl-2 combination occurs naturally in S. mansoni. In contrast, the evidence did not support the formation of heteromeric SmGluCl-1:SmGluCl-3 receptors. The EC50 and Hill coefficient for L-glutamate on SmGluCl-1:SmGluCl-3-injected oocytes were 7.0±0.7 µM and 1.0±0.1, respectively (n = 8, Figure 4B), indicating that co-injection of SmGluCl-1:SmGluCl-3 cRNAs did not alter sensitivity to L-glutamate compared to homomeric SmGluCl-3 receptors (P>0.05, paired-T test).
SmGluCl-2 and SmGluCl-3 receptors form Cl− permeable channels
In anionic LGICs, −2′ proline, −1′ alanine and 13′ threonine residues from the M2 transmembrane domain are the minimal determinants of anion selectivity [23], [24]. These three residues are conserved in the SmGluCl-1, SmGluCl-2, SmGluCl-3 and SmGluCl-4 subunits (Figure S1), which are therefore predicted to be permeable to Cl−. To confirm that the SmGluCls form Cl− channels, we performed current-voltage (I/V) relationship experiments on oocytes expressing homomeric SmGluCl-2.1 and SmGluCl-3 receptors. For SmGluCl-2.1, we compared the reversal potentials for glutamate-sensitive currents in normal ND96 solution containing 103.6 mM Cl− and 96 mM Na+, sodium gluconate-replaced ND96 with [Cl−] reduced to 45.6 mM, and choline chloride-replaced ND96 with [Na+] reduced to 38 mM. In oocytes injected with SmGluCl-2.1 cRNA, reduction of external [Cl−] shifted the reversal potential by 15.5±2.9 mV (n = 4), from −23.0±2.3 mV to −8.2±2.4 mV (P<0.0001, one-way ANOVA), reasonably close to the 20.7 mV shift predicted by the Nernst equation (Figure 5A), indicating that SmGluCl-2.1 receptors are Cl− selective. In contrast, reducing the extracellular [Na+] did not alter the reversal potential (−22.0±3.0 mV, n = 4) or the I/V relationship compared to normal ND96, indicating that SmGluCl-2.1 receptors are not permeable to Na+. We measured the reversal potentials of SmGluCl-3 receptors for glutamate-sensitive currents in normal ND96, sodium gluconate-replaced ND96 with [Cl−] reduced to 13.6 mM, and choline chloride-replaced ND96 containing 6 mM Na+. In oocytes injected with SmGluCl-3 cRNA, reduction of the external [Cl−] shifted the reversal potential by 37.4±2.8 mV (n = 4), from −15.6±9.8 mV to 21.9±9.2 mV (P<0.0001, one-way ANOVA), reasonably close to the 51.5 mV shift predicted by the Nernst equation (Figure 5B), indicating that SmGluCl-3 receptors are permeable to Cl−. Reducing the extracellular [Na+] had no influence on the I/V relationship or the reversal potential (−17.7±4.9 mV, n = 4) compared to normal ND96, indicating that SmGluCl-3 receptors are not permeable to Na+. Our electrophysiology data clearly demonstrate that SmGluCl-2.1 and SmGluCl-3 receptors are permeable to Cl−, as predicted from their M2 regions. These results can be extrapolated to all four S. mansoni GluCl subunits (SmGluCls), since they all possess the canonical molecular determinants for anion selectivity.
SmGluCl-1, SmGluCl-2 and SmGluCl-3 are insensitive to ivermectin
Macrocyclic lactones (MLs) such as ivermectin (IVM) are potent agonists of nematode and arthropod GluClα receptor subtypes, and also potentiate the response to L-glutamate [29]. The interaction of IVM with GluCl receptors is strongly correlated with its biological activity, and is therefore considered to be the main contributor to its anthelmintic and insecticidal effects [34]. In contrast, flatworms are insensitive to MLs and it has been postulated that they lack high-affinity ML binding sites [35], [36]. We examined whether IVM had agonist or modulatory effects on SmGluCl-1, SmGluCl-2.1, and SmGluCl-3 receptors. Application of 1 µM IVM induced no detectable channel activity, either by direct activation or by potentiation of L-glutamate, in oocytes injected with SmGluCl-1, SmGluCl-2.1, and SmGluCl-3 cRNAs alone or in combination (data not shown). The insensitivity of S. mansoni GluCls to IVM, even at a high concentration, confirms that they are pharmacologically distinct from the highly IVM-sensitive GluClα receptors found in nematodes and arthropods.
The crystal structure of the Caenorhabditis elegans GluClα receptor (GLC-1) complexed with IVM [37], together with site-directed mutagenesis and molecular modeling studies on human GlyRα1 [38], demonstrated that the ML binding site is located in the TMD interface of adjacent subunits. Analysis of the SmGluCl TMD residues, including SmGluCl-4, reveals that they lack key molecular determinants for ML binding. Most importantly, the combination of a proline residue in M1 and a glycine residue in M3 (Pro284 and Gly342 in C. elegans GLC-1, respectively) was shown to be essential for high IVM sensitivity, and these two residues are entirely conserved in all GluClα subunits [39]. In the SmGluCls, the M3-Gly residue is substituted by bulky residues (Figure S1), thereby disrupting a crucial interaction with IVM. Our electrophysiology results corroborate the prediction of Lynagh et al [29] that the large side chains replacing this M3-Gly residue in the schistosome putative inhibitory Cys-loop LGIC subunits would be sufficient to preclude IVM binding. The authors also pointed out that potent IVM activity requires polar residues at key positions in M1 and M2; hydrophobic residues at these positions could adversely affect IVM binding, even if an M3-Gly is present. More specifically, IVM appears to interact with a glutamine residue located in M1 (Gln280 in C. elegans GLC-1), and a polar residue at position M2–12′, -14′, -15′ or -19′ (Thr318, Gln320, Ser321 and Asn325 in C. elegans GLC-1, respectively). These interactions are likely to be disrupted in the SmGluCl subunits, since the amine group of the M1-Gln is substituted by hydroxylated residues, the M2–12′ Thr is replaced with hydrophobic residues, and the M2–19′ Asn residues is substituted by either hydrophobic or hydroxylated residues (Figure S1). These substitutions in M2 might create a more hydrophobic environment, which could be less favourable for IVM binding.
The electrophysiological data presented here, combined with current knowledge pertaining to the molecular basis of IMV sensitivity, provide further evidence that the insensitivity of schistosomes to ML is due to the lack of affinity of S. mansoni GluCls for IVM. To extend this observation to other parasitic flatworms, we performed a protein BLAST of the S. haematobium [40], S. japonicum, Clonorchis sinensis [41], Echinococcus multilocularis and Hymenolepis microstoma [42] genome databases, and identified 21 GluCl candidates homologous to the SmGluCls. We found that the same substitutions described above, likely responsible for the SmGluCls' lack of affinity for IVM, are conserved in all flatworm GluCl-like sequences examined (data not shown). This suggests that all these flatworm putative GluCls lack high affinity for IVM, providing a molecular explanation as to why flatworms are insensitive to MLs.
SmGluCl-1, SmGluCl-2, and SmGluCl-3 are insensitive to meclonazepam
As noted above, BZDs are classical allosteric modulators of mammalian inhibitory Cys-loop GABA receptors, and are extensively used in human medicine [19]. In mammals, BZDs also appear to have modulatory effects on GlyRs [43], [44]. In addition, BZDs can act as low-affinity allosteric inhibitors of α2-containing GlyRs in rats [45] and are active at insect GABA receptors[46]–[49]. Schistosomicidal activity against S. mansoni and S. haematobium has been demonstrated in a subset of BZDs at doses close to the therapeutic range, including meclonazepam and clonazepam [25], [27]. Schistosomicidal BZDs have physiological effects reminiscent of neuronal modulation, but their molecular target has yet to be identified. We hypothesized that meclonazepam might act on schistosome receptors related to mammalian Cys-loop GABA or Gly receptors. We asked whether SmGluCls, the only inhibitory LGIC found in S. mansoni, could be the molecular target of schistosomicidal BZDs. We therefore examined if meclonazepam had agonist or modulatory effects on SmGluCl-1, SmGluCl-2.1, and SmGluCl-3 receptors. Application of 10 µM meclonazepam had no direct effect and failed to potentiate L-glutamate-sensitive currents in oocytes injected with any of the subunits, either alone or in combination (data not shown). These results suggest that SmGluCls are not likely to be the molecular target of schistosomicidal BZDs. However, one cannot rule out the possibility that SmGluCl-4 is required to confer sensitivity to meclonazepam in SmGluCl assemblies as an alternative explanation for the lack of modulatory effects by the drug on Xenopus oocytes expressing the available subunits.
Alternatively, it is possible that meclonazepam targets previously uncharacterized proteins unrelated to mammalian BZD receptors. A previous study identified a low-affinity binding site for meclonazepam in S. mansoni membrane extracts, with Kd = 2 µM, but its pharmacological properties were considerably different from BZD binding sites in mammalian brain [50]. The nature of this binding site remains unknown. Noël et al [51] identified two BZD binding sites exhibiting pharmacological properties partly reminiscent of central BZD receptors (GABA-gated Cl− channel) and non-neuronal peripheral BZD receptor [52], respectively, in crude membrane extracts of S. mansoni. However, these BZD binding sites were discarded as putative molecular targets for meclonazepam, since neither exhibited significant affinity for meclonazepam [51], [53]. In agreement with the evidence for a low-affinity meclonazepam binding site, a recent study in in vitro cultures of S. mansoni reported half maximal lethal concentrations (LC50) of 3 µM for meclonazepam and 10 µM for clonazepam [28]. These LC50 values correlate with in vitro studies on S. mansoni adult males, which showed that 1–10 µM meclonazepam induces immediate, Ca2+-dependent spastic paralysis and extensive tegumental disruption, leading to parasite death [54]–[56]. Some classical Ca2+ channel blockers partially blocked the Ca2+-dependant paralysis and subsequent death induced by meclonazepam [57]. This raises the possibility that the molecular target of meclonazepam in schistosomes is involved, directly or indirectly, in Ca2+ channel functions or Ca2+ homeostasis. Interestingly, modulation of Ca2+ channel activity by BZDs has also been reported in mammals, also at µM concentrations [58]. This modulatory role does not seem to be associated with central or peripheral BZD receptors, and it is unclear whether it is mediated by a specific interaction with a low-affinity BZD receptor, or if it is a consequence of unspecific modulation of processes in excitable cells. [25], [59], [60].
However, discrepancies between in vitro and in vivo potencies of schistosomicidal BZDs suggest that the effect of meclonazepam on intracellular calcium in schistosomes cultured in vitro may not be the primary mode of action by which schistosomicidal BZDs exert their therapeutic effect in vivo. Baard et al [25] reported that a single dose of 0.2 to 0.3 mg/kg meclonazepam was curative in patients infected with S. haematobium or S. mansoni. Pharmacokinetic studies following single oral doses of meclonazepam and clonazepam show that plasma levels of the drugs corresponded to <3.5% of the original doses, and the proportion found in plasma tended to decrease as the dose increased [59]–[61]. Based on these figures, Cmax estimates for schistosomicidal doses of meclonazepam are likely to be < 600 nM, below the concentrations that affect schistosomes in culture. Although spastic paralysis may be due to binding to the low-affinity meclonazepam binding site in S. mansoni, schistosomicidal activity observed in vivo may be attributable to mechanisms mediated by unidentified receptors.
SmGluCl-1, SmGluCl-2, and SmGluCl-3 belong to a novel family of GluCl receptors pharmacologically distinct from their nematode, arthropod and mollusc counterparts
Heterologous expression of SmGluCls revealed pharmacological properties distinct from insect, nematode and snail GluCls. Our data indicate that SmGluCls exhibit much higher affinity for L-glutamate than their GluCl counterparts in free-living nematodes and molluscs. Interestingly, the SmGluCls resemble the insect and parasitic nematode GluClα clade in this regard, but are distinguished from them by their insensitivity to IVM. We measured EC50 values for L-glutamate of 6.9 µM and 11.8 µM for SmGluCl-3 and SmGluCl-2.1 homomeric receptors, respectively, and 26.3 µM for SmGluCl-1:SmGluCl-2.1 heteromeric receptors. Drosophila melanogaster GluClα EC50 values ranged from 19.5 to 23 µM [62], [63], while Musca domestica GluClα EC50 values ranged from 30 to 40.3 µM [64], [65]. In Haemonchus contortus, EC50 values of 8.4 µM and 27.6 µM were determined for homomeric GluClα (HcoGLC-5a) and GluClα3 (HcoGBR-2b, ortholog of C. elegans AVR-14b), respectively [66], [67]. Cooperia oncophora homomeric GluClα3 (ortholog of C. elegans AVR-14b) had an EC50 of 29.7 µM, whereas heteromers of GluClα3 and GluClβ had an EC50 of 13.4 µM [68]. However, 1 µM IVM had no effect on oocytes expressing any SmGluCl, indicating that they do not align with the high IVM affinity GluClα subtype group. On the other hand, SmGluCl are more similar to the GluCls from the snail Aplysia californica and GluClβ receptors from nematodes with regards to their insensitivity to IVM, but are 7- to 54- fold more sensitive to L-glutamate than these receptors. A. californica GluClAc1 and GluClAc2 EC50 values ranged from 196 to 499 µM [69]. C. elegans homomeric GluClβ had an EC50 of 380 µM, while C. oncophora homomeric GluClβ had an EC50 of 185.6 µM [68], [70].
The focus of the present study was to characterize a novel family of inhibitory GluCls in S. mansoni. Additional experiments in Xenopus oocytes are needed to expand the pharmacological profile of these receptors. For instance, unlike the nematode and arthropod GluCls previously characterized, no selective agonists of S. mansoni GluCls have been described. The SmGluCls lack high affinity for IVM, and activity of nodulisporic acid, a selective agonist of arthropod GluCls [71], has not been reported in S. mansoni. Identification of selective agonists of the SmGluCls would provide a very useful tool for RNAi-based studies examining the functions of these receptors in the worm and their potential as drug targets, in addition to potentially opening the way to the discovery of new schistosomicidal compounds.
The S. mansoni GluCl subunits belong to a previously uncharacterized clade of flatworm GluCl subunits
Seventeen Cys-loop LGIC subunits are predicted in the S. mansoni genome: the four inhibitory subunits described here and 13 other putative nAChR-like Cys-loop LGIC subunits [14]. We first examined the phylogenetic relationship between the four inhibitory Cys-loop LGIC subunits cloned from S. mansoni and representatives of vertebrate and invertebrate inhibitory and excitatory Cys-loop LGIC subunits (Figure 6A). The phylogenetic analysis included inhibitory subunits from C. elegans, H. contortus, D. melanogaster, the snail species A. californica, Haliotis asinina and Lymnaea stagnalis, and Rattus norvegicus. Several subunits homologous to the SmGluCls were identified by protein BLAST search against the genomes of related trematode and cestode species: S. haematobium [40], [72], S. japonicum [73], C. sinensis [41], E. multilocularis and H. microstoma [42]. We examined the phylogenetic relationship between the putative inhibitory Cys-loop LGIC subunits identified in these flatworms, the inhibitory Cys-loop LGIC subunits cloned from S. mansoni, and GluCl subunits from the snail, nematode and insect species mentioned above (Figure 6B).
The most striking finding was the evolutionary origin of the SmGluCls, which clearly belong to a novel clade of glutamate-gated Cl− channels distinct from GABA receptors, with branch reliability estimates >95% (aLRT test, Chi2-based parametric; Figure 6A). The separation from their counterparts in Ecdysozoa and Lophotrochozoa was significant, with branch reliability estimates of 94% and >95%, respectively (aLRT test, Chi2-based parametric; Figure 6A). In addition, SmGluCl homologs identified in other trematode and cestode species reveal an extensive distribution of putative GluCl-like subunits. Flatworm GluCls appear orthologous to the S. mansoni GluCls, with branch reliability estimates >95% (aLRT test, Chi2-based parametric; Figure 6B). This finding suggests a common ancestral origin of trematode and cestode GluCl subunits, and corroborates the separation of flatworm GluCls from Ecdysozoa and Lophotrochozoa (mollusc) GluCls. This conclusion is supported by branch reliability estimates >80% when bootstrapping (100 replicates) is used (Figure S2). The flatworm GluCl subunit clade appears to originate from an ancestral gene closely related to the Ecdysozoa GluCl clade, and distinct from the molluscan GluCl clade based on the GluClAc1 and GluClAc2 subunits from A. californica, and a putative GluCl subunit (CAA50492) from L. stagnalis, described by Kehoe et al [69]. The genome of the oyster Crassostrea gigas [74] encodes subunits only from the molluscan GluCl clade, suggesting conservation of glutamate signaling among molluscs. An explanation for why the lophotrochozoan S. mansoni GluCls would be more similar to the ecdysozoan GluCls than the lophotrochozoan mollusc GluCls is that the common ancestor encoded both avr-14-like (ecdysozoan) and GluClAc-like (molluscan) GluCls. In this model, the Ecdysozoa have lost the GluClAc-like class; within the Lophotrochozoa, the molluscs have lost the avr-14-like GluCl, whereas flatworms have lost the GluClAc-like GluCl.
Another important finding was that the SmGluCls are the only non-nAChR Cys-loop LGIC subunits predicted to be present in S. mansoni (Figure 6A). Likewise, in preliminary phylogenetic analyses of the LGIC subunits from other trematode and cestode species (not shown), all putative non-nAChR Cys-loop LGIC subunits identified were SmGluCl homologs that grouped with the flatworm GluCl clade, and were distinct from GABA receptors. Collectively, this suggests that parasitic flatworms lack Cys-loop GABA receptors, and raises the possibility that Cys-loop LGIC-mediated GABA signaling does not occur in cestodes and trematodes. The lack of bioinformatic evidence supporting the existence of GABA receptors is consistent with the fact that GABA-related physiological activity has not been confirmed experimentally in any of these species. In fact, chemotherapeutic effects of Cys-loop GABA receptor modulators have not been described in flatworms, including schistosomes.
The data raise questions about whether the canonical GABAergic signaling pathway exists in parasitic flatworms, including S. mansoni, and are difficult to reconcile with evidence for an extensive distribution of GABA in the central and peripheral nervous system of some trematode and cestode species. GABA immunoreactivity has been demonstrated in the nervous system of S. mansoni [75], as well as in the trematode Fasciola hepatica, the cestode Moniezia expansa, and the free-living turbellarians Polycelis nigra and Dugesia tigrina [76]–[78]. However, GABA distribution has not been documented in S. japonicum, S. haematobium, C. sinensis, E. multilocularis or H. microstoma. Glutamate decarboxylase, the enzyme that produces GABA from glutamate, is found in S. mansoni homogenates [75] and transcripts encoding putative GABA transporters and GABA receptor-associated protein are also present [79]. Evidence for GABA biosynthesis was also reported in D. tigrina homogenates [77], but has not been demonstrated in parasitic flatworms.
The SmGluCl family represents a novel, inhibitory component of the glutamatergic transmitter system in S. mansoni
Our findings constitute the first report of inhibitory GluCl receptors from platyhelminths, and give new insights into the glutamatergic transmitter system in S. mansoni. While the genome encodes enzymes involved in L-glutamate biosynthesis and transport, as well as several excitatory ionotropic glutamate-gated channels (iGluRs) and metabotropic glutamate receptors (mGluRs), no inhibitory GluCl subunits were predicted in this parasite [14]. The glutamate receptor subunits predicted from the S. mansoni genome include three putative kainate-like iGluR subunits, two putative NMDA-like iGluR subunits and two putative AMPA-like iGluR subunits, as well as three putative mGluR subunits.
Despite genomic evidence for glutamatergic neurotransmission, studies describing L-glutamate physiological actions are scarce in S. mansoni, and most of the putative glutamate receptors await functional characterization. L-glutamate-containing neurons have been identified in S. mansoni cercaria [80], but have not been reported in other life stages. L-glutamate has been implicated in the regulation of neuromuscular activity in schistosomes, but the mechanisms mediating these effects are unclear. Direct application of L-glutamate did not affect basal motor activity in intact worms [81], but a subsequent study provided evidence of glutamate-induced, concentration-dependent contractions mediated by a glutamate/Na+ co-transporter in isolated S. mansoni muscles fibers [82]. Kainate, an iGluR agonist, was subsequently shown to produce tonic muscle contraction and paralysis in whole worms [83]. In addition, a low-affinity kainate binding site was found to exhibit pharmacological properties consistent with excitatory glutamatergic receptors in a crude membrane extract of S. mansoni, providing indirect evidence for a kainate-like iGluR. No NMDA-like iGluRs have been reported in S. mansoni, but a phenotypic chemical screen identified the NMDA iGluR antagonist 7-chlorokynurenic acid as an inhibitor of in vitro miracidial transformation in schistosomes [84]. Two S. mansoni mGluR subunits, SmGluR (Smp_128940) and SmGBP (Smp_052660), have recently been characterized [85], [86]. Both subunits are distantly related to mGluRs from other species and have homologs in other flatworms [87]. Despite recent progress in this area, the functional roles of glutamate in schistosome neuromuscular systems remain ill-defined.
The rich diversity of glutamate receptors in S. mansoni argues for a more important role of the glutamatergic signaling pathway in parasitic flatworms than previously appreciated. The characterization of an inhibitory SmGluCl family is unprecedented in platyhelminths, and has important implications in S. mansoni biology, as well as in other parasitic flatworms encoding GluCl-like receptors. Further investigations are needed to clarify the functions of L-glutamate and the physiological relevance of SmGluCls, as well as their potential as a therapeutic target. Immunolocalization experiments and functional studies are in progress to examine the distribution and physiological roles of these receptors in S. mansoni.
Materials and Methods
Ethics statement
All animal procedures were approved by the McGill University Facility Animal Care Committee (FACC), in full compliance with the Canadian Council on Animal Care. Recovery of S. mansoni adult worms from mice was performed in accordance to the FACC animal protocol # 3346. Oocytes were harvested from mature X. laevis females in accordance to the FACC animal protocol # 5284.
Parasites
S. mansoni adult worms (Puerto Rican strain) were kindly provided by Dr. P. Ribeiro (McGill University). Sporocyst-infected Biomphalaria glabrata snails were obtained from Dr F. Lewis, Biomedical Research Institute (Bethesda, MD, USA) and cercarial shedding was induced by light exposure as previously described [88]. 28 day-old CD1 female mice were infected with approximately 150 cercaria/mouse and adult S. mansoni worms were recovered by portal perfusion 6 to 8 weeks later.
Identification of putative inhibitory Cys-loop LGIC in S. mansoni
The S. mansoni draft genome database was comprehensively searched to identify inhibitory Cys-loop LGIC gene candidates sharing homology with Cys-loop GABA receptor subunits. Smp_015630, Smp_096480, Smp_104890, Smp_099500 and Smp_176730 putative genes were identified by performing protein BLAST analyses of the S. mansoni genome database (version 3.1), using Cys-loop GABA receptor sequences from mammals and invertebrates as queries. Data retrieved from the S. mansoni draft genome database were produced at the Sanger Institute by the Schistosoma mansoni Sequencing Group. Smp_015630, Smp_096480, Smp_104890, Smp_099500 and Smp_176730 putative genes sequences are identical in later versions of the genome [13], [14].
The topology of Smp_015630, Smp_096480, Smp_104890, Smp_099500 and Smp_176730 putative amino acid sequences was analysed for functional domains, hydrophobic regions and signature motifs characteristic of the ligand-gated ion channel superfamily. SMART (http://smart.embl-heidelberg.de/) [31] and InterProScan (http://www.ebi.ac.uk/Tools/pfa/iprscan/) [32] analyses were used to predict the domain architecture and hydrophobic regions. The cysteine-loop signature motif was identified using ScanProsite (http://prosite.expasy.org/) [33]. Gene structures were drawn using the Exon-Intron Graphic Maker version 4 (wormweb.org/exonintron).
RACE experiments and cloning of SmGluCl-1, SmGluCl-2, SmGluCl-3 and SmGluCl-4
The sequences were extended to full length using 5′ and 3′RACE procedures (Rapid Amplification of cDNA Ends). For 5′RACE, a S. mansoni cDNA library (kindly provided by Dr P. LoVerde) was used as template in a PCR with a sense primer corresponding to the T3 promoter region in the pCMV-Script EX vector (5′-AATTAACCCTCACTAAAGGG-3′), used in conjunction with specific internal antisense primers (Table S1). Second rounds of nested PCRs were done, using the same T3 sense primer and nested specific internal antisense primers (Table S1). 3′RACE experiments were done using a commercial 3′RACE kit (Bioline), and first-strand cDNA obtained from reverse-transcribed adult S. mansoni total RNA was used as a template. Total RNA was extracted from adult S. mansoni worms using the TRIzol reagent (Invitrogen), and then purified using the RNeasy Plus Mini kit (Quiagen). Briefly, total RNA was extracted in TRIzol reagent according to manufacturers' instructions. The aqueous layer containing the total RNA fraction was recovered and mixed with the RLT buffer supplied whit the RNeasy kit, and the RNA purification proceeded according to manufacturers' instructions. RNA was reverse-transcribed using an oligo(dT)20 primer and Superscript Transcriptase III (Invitrogen). First-strand cDNA was synthesized at 55°C according to standard procedures, and used as a template for 3′RACE. A touch-up PCR was done, using the antisense 3′RACE adaptor supplied with the 3′RACE kit, in conjunction with specific internal sense primers (Table S2). Second rounds of nested PCR were done, using the 3′RACE outer primer supplied with the kit and nested specific internal sense primers (Table S2). The 5′ and 3′ RACE products were cloned into pJET1.2/blunt vector (Fermentas) and new ends were confirmed by DNA sequencing.
Once the full-length coding sequences were confirmed, SmGluCl-1 (Smp_096480), SmGluCl-2.1 (Smp_015630 isoform 1), SmGluCl-2.2 (Smp_015630 isoform 2), and SmGluCl-3 (Smp_104890) open reading frames (ORFs) were amplified directly from first-strand cDNA, and cloned into the pJET1.2/blunt vector (Fermentas). PCR reactions were done using specific primers targeting the predicted start and stop codons (Table S3). Primers were flanked with BglII, BamHI, EcoRV or SpeI restriction sites, to allow subsequent sub-cloning into the Xenopus expression vector pT7TS (AB255037; www.addgene.org/17091/) which flanks the insert with 5′ and 3′ untranslated region sequences from the X. laevis β-globin gene . The pT7TS multiple cloning site contains the following restriction sites: BglII, EcoRV and SpeI. BglII and BamHI have compatible ends. Product sequences were confirmed by DNA sequencing of pJET constructs for each gene.
SmGluCl-1, SmGluCl-2.1, SmGluCl-2.2, and SmGluCl-3 cloned sequences were transferred into the pT7TS Xenopus expression vector. The cloned products were digested out of the pJET1.2 constructs using the appropriate restriction enzymes (Table S3), and then sub-cloned into the pT7TS vector. Sub-cloning was confirmed by DNA sequencing of pT7TS constructs for each gene.
In vitro capped RNA synthesis
Capped RNA (cRNA) was synthesized in vitro using the mMESSAGE mMACHINE T7 polymerase kit (Ambion). The pT7TS constructs were used as templates. Briefly, the pT7TS constructs were linearized downstream of the X. laevis β-globin 3′UTR using XbaI and purified using the GeneJET PCR purification kit (Fermentas). The capped transcription reaction was assembled according to manufacturer's instructions and incubated at 37°C for 2 hrs. Template plasmid DNA was removed using the TURBO DNAse supplied with the kit, and capped RNA was recovered by LiCl precipitation and resuspended in nuclease-free H2O.
Oocyte expression and electrophysiology
Oocytes were harvested from mature X. laevis females using standard procedures [89] and maintained at 15°C in ND96 solution (96 mM NaCl, 2 mM KCl, 1.8 mM CaCl2, 1 mM MgCl2, and 5 mM Hepes, pH 7.5), as described [90]. Oocytes were injected with 35–50 ng cRNA in a total volume of 46 nl using the Nanoject system (Drummond Scientific, Broomall, PA) and incubated for 2 to 3 days in ND96 at 15°C before recordings were done. To express heteromeric receptors, subunit cRNAs were mixed at a ratio of 1∶1 before injection into oocytes. Uninjected oocytes and oocytes injected with nuclease-free H2O were used as negative controls.
Two-electrode voltage clamp (TEVC) experiments were carried out using a fast perfusion system and a Maltese Cross chamber (ALA Scientific Instruments, Westbury, NY). Microelectrode pipettes were filled with 3 M KCl and had a resistance between 1 and 3 MΩ. The bath was connected to the ground through a 3 M KCl agar bridge. Pharmacological compounds were purchased from Sigma, with the exception of meclonazepam, which was kindly provided by Roche. All experiments were carried out in ND96, and agonists and drugs were dissolved or diluted in ND96. In the case of IVM and meclonazepam, 10 mM stock solutions were prepared in DMSO, and diluted to the appropriate working concentrations in ND96. Final DMSO concentration in ND96 did not exceed 0.1%. Measurements were done using the AxoClamp 2B and Digidata 1322A 16-bit data acquisition system (Axon Instruments, Foster City, CA), and recordings were sampled at 1 kHz using Clampex 8.1 digital oscilloscope software (Axon Instruments, Foster City, CA). Data were filtered at 30 Hz and analysed using GraphPad Prism software (GraphPad Software, San Diego, CA, USA).
For L-glutamate concentration-response experiments, oocytes were clamped at −80 mV. Concentration-response curves were generated by applying increasing concentrations of L-glutamate, with a 2 to 3 min interval between applications to ensure full recovery from desensitization. For each oocyte, responses were normalized to a saturating concentration of L-glutamate. Data were fitted to a sigmoid concentration-response curve of variable slope using the equation:
where I is the response, Imax is the estimated maximal response, Imin is the estimated minimal response, EC50 is the concentration of agonist eliciting half-maximal response, [a] is the agonist concentration, and n is the Hill coefficient. EC50 values and Hill coefficients correspond to the mean ± SEM for 8 individual oocytes from 3 frogs. Data points on graphs are given as the mean ± SD. Statistical analyses were performed using paired two-tailed t-test, with a significance level of P<0.05.
SmGluCl-2.1 receptors desensitized almost completely in the continued presence of L-glutamate. As a consequence, I/V experiments could not be performed using a voltage ramp protocol. Instead, I/V relationships were determined by measuring the peak glutamate currents obtained at holding potentials ranging from −60 mV to 20 mV, and the ND96 Cl− and Na+ concentrations were adjusted accordingly for the reversal potentials to fall within that range. At holding potentials > 20 mV, the I/V relationship of SmGluCl-2.1 glutamate-sensitive currents no longer behaved in a linear fashion, due to the contribution from Ca2+-activated currents. Ion substitution was accomplished by replacing NaCl in standard ND96 with 58 mM choline chloride, or 58 mM sodium gluconate. Normal ND96 contained 103.6 mM Cl− and 96 mM Na+. Sodium gluconate-replaced ND96 had a Cl− concentration reduced to 45.6 mM, and choline chloride-replaced ND96 had a Na+ concentration reduced to 38 mM. Oocytes expressing SmGluCl-3 did not sustain voltage clamping for extended periods of time, preventing us from performing I/V experiments using the protocol applied for SmGluCl-2.1 receptors. SmGluCl-3 receptors desensitized at a slower rate than SmGluCl-2.1 in the continued presence of L-glutamate, allowing the use of a voltage ramp protocol to determine I/V relationships. For SmGluCl-3, I/V relationships were obtained using a 4 mV/s voltage ramps in the presence and absence of 5 µM L-glutamate. Currents obtained over a voltage range of −60–40 mV were generated by subtracting L-glutamate-free from L-glutamate-containing data. Ion substitution was accomplished by replacing NaCl in standard ND96 with 90 mM choline chloride, or 90 mM sodium gluconate. Normal ND96 contained 103.6 mM Cl− and 96 mM Na+. Sodium gluconate-replaced ND96 had a Cl− concentration reduced to 13.6 mM, while choline chloride-replaced ND96 saline contained 6 mM Na+. Data points on graphs are given as the mean ± SD of 4 oocytes. Statistical analyses of I/V relationships were performed using one-way ANOVA test, with a significance level of P<0.05.
Phylogeny
We analysed the evolutionary relationship between SmGluCl-1, SmGluCl-2, SmGluCl-3, and SmGluCl-4 from S. mansoni and other invertebrates and vertebrates inhibitory Cys-loop LGIC subunits. The analyses included subunits from rat (Rattus norvegicus, Rn), nematodes (Caenorhabditis elegans, Ce; Haemonchus contortus, Hco), the fruit fly Drosophila melanogaster (Dm), snails (Aplysia californica, Ac; Haliotis asinine, Ha; Lymnaea stagnalis, Ls), trematodes (Schistosoma haematobium, Sha; Schistosoma japonicum, Sjp; Clonorchis sinensis, Csi), and cestodes (Hymenolepis microstoma, HmN; Echinococcus multilocularis, EmuJ). In Figure 2A, reference was given to subunits for which the function had been confirmed by heterologous expression. Five subunits from D. melanogaster, three subunits from C. elegans and four subunits from snails that clearly group with the inhibitory Cys-loop LGIC family, but for which function has not yet been characterized, were also included. The cation-selective GABA-gated channel subunits from C. elegans and D. melanogaster were also included in the analysis. To build both trees, subunit amino acid sequences were first aligned using PROMALS-3D [91]. To optimize the accuracy of the alignment, alignment of sequences within groups in the first stage was done using the Promals algorithm. Default settings were used for all other parameters. Nonalignable, noninformative sites were removed manually using Jalview 2 (304 informative sites located in the LBD and M1, M2 and M3 were used for tree building) [92].
The phylogenetic analysis was performed on the Phylogeny.fr platform (http://www.phylogeny.fr/) [93], [94]. The phylogenetic tree was reconstructed using the maximum likelihood method implemented in the PhyML program (v3.0 aLRT) [95]. The WAG substitution model was selected assuming an estimated proportion of invariant sites (of 0.023 and 0.041 for tree in Figure 2A and 2B, respectively) and 8 gamma-distributed rate categories to account for rate heterogeneity across sites. The gamma shape parameter was estimated directly from the data (gamma = 1.563 and gamma = 1.048 for tree in Figure 2A and 2B, respectively). Reliability for internal branch was assessed using the aLRT test (Chi2-based parametric) [96], where values > 95% are considered reliable. Graphical representation and edition of the phylogenetic tree were performed with TreeDyn (v198.3) [97]. Alternatively, reliability for internal branch was assessed using the bootstrapping method (100 bootstrap replicates).
LGIC sequences and accession numbers
S. mansoni putative LGIC subunit sequences were acquired from the S. mansoni genome database (version 3.1), and are available in the current version of the S. mansoni genome (http://www.genedb.org/Homepage/Smansoni; v5.0). Nucleotide sequences of the cloned SmGluCl subunits were submitted to the GenBank database and were assigned the following accession numbers: SmGluCl-1 – KC861381, SmGluCl-2.1 – KC861382, SmGluCl-2.2 – KC861383, SmGluCl-3 – KC861384, and SmGluCl-4 (partial ORF) – KC861385.
C. elegans sequences were obtained from Wormbase (http://www.wormbase.org version WS230, march 2012) [98]. The D. melanogaster sequences were obtained from Flybase (http://flybase.org, database issue D706-14, march 2012) [99]. R. norvegicus, H. contortus, A. californica, H. asinina and L. stagnalis sequences with the following accession numbers were acquired from the NCBI Protein Database: RnGABRα1 – P62813, RnGABRβ2 – P63138, RnGlyRα1 – NP_037265, RnGlyRβ – P20781, HcoGLC-5a – AF076682, HcoGLC-6 – ABV68895, AcGABRα (putative) – AAL37252, AcGluClAc1 – GQ148562, AcGluClAc2 – GQ148563, HaGABR (putative) – ABN13868, LsGABRβ (putative) – P26714 and LsGABRζ (putative) – CAA50492.
Putative genes homologous to SmGluCls in parasitic flatworms were identified by performing protein BLAST analyses of the S. haematobium, S. japonicum, C. sinensis, H. microstoma and E. multilocularis genome databases. The translated sequence of SmGluCl-1, SmGluCl-2.1 and SmGluCl-3 full-length cDNAs, and SmGluCl-4 partial cDNA were used as queries. S. haematobium sequences were obtained from the genome database on SchistoDB 3.0 (http://schistodb.net/schisto/, database issue D579-82) [40], [72]. S. japonicum sequences were obtained from the genome database on GeneDB (http://www.genedb.org/Homepage/Sjaponicum) [73]. Sequences from C. sinensis [41] with the following accession numbers were acquired from the NCBI Protein Database: CsiGlyRα1 (putative) – GAA34110.2, CsiGlyRα2 (putative) – GAA51498.1, CsiGlyRαZ1 (putative) – GAA38381.2, CsiGlyRβ (putative) – GAA28305.1, and CsiGlyRβ (putative) – GAA36399.1. Sequences from H. microstoma were obtained from the draft genome databases on GeneDB (http://www.genedb.org/Homepage/Hmicrostoma) [42]. Sequences from E. multilocularis were obtained from the draft genome databases on GeneDB (http://www.genedb.org/Homepage/Emultilocularis). Data retrieved from the E. multilocularis draft genome database were produced at the Sanger Institute by the Echinococcus multilocularis Sequencing Group, in collaboration with Prof. Klaus Brehm of Institute for Hygiene and Microbiology, University of Wurzburg.
Supporting Information
Zdroje
1. GryseelsB, PolmanK, ClerinxJ, KestensL (2006) Human schistosomiasis. Lancet 368: 1106–1118.
2. SteinmannP, KeiserJ, BosR, TannerM, UtzingerJ (2006) Schistosomiasis and water resources development: systematic review, meta-analysis, and estimates of people at risk. Lancet Infect Dis 6: 411–425.
3. KingCH, DickmanK, TischDJ (2005) Reassessment of the cost of chronic helminthic infection: a meta-analysis of disability-related outcomes in endemic schistosomiasis. Lancet 365: 1561–1569.
4. KingCH, Dangerfield-ChaM (2008) The unacknowledged impact of chronic schistosomiasis. Chronic Illn 4: 65–79.
5. van der WerfM, VlasS, BrookerS, LoomanC, NagelkerkeN, et al. (2003) Quantification of clinical morbidity associated with schistosome infection in sub-Saharan Africa. Acta Trop 86: 125–139.
6. IsmailM, MetwallyA, FarghalyA, BruceJ, TaoL-F, et al. (1996) Characterization of isolates of Schistosoma mansoni from egyptian villagers that tolerate high doses of praziquantel. Am J Trop Med Hyg 55: 214–218.
7. GryseelsB, MbayeA, De VlasSJ, StelmaFF, GuisséF, et al. (2001) Are poor responses to praziquantel for the treatment of Schistosoma mansoni infections in Senegal due to resistance? An overview of the evidence. Trop Med Int Health 6: 864–873.
8. BotrosS, SayedH, AmerN, El-GhannamM, BennettJL, et al. (2005) Current status of sensitivity to praziquantel in a focus of potential drug resistance in Egypt. Int J Parasitol 35: 787–791.
9. MelmanSD, SteinauerML, CunninghamC, KubatkoLS, MwangiIN, et al. (2009) Reduced susceptibility to praziquantel among naturally occurring Kenyan isolates of Schistosoma mansoni. PLoS Negl Trop Dis 3: e504.
10. DoenhoffMJ, HaganP, CioliD, SouthgateV, Pica MattocciaL, et al. (2009) Praziquantel: its use in control of schistosomiasis in sub-Saharan Africa and current research needs. Parasitology 136: 1825–1835.
11. HaltonDW, GustafssonMKS (1996) Functional morphology of the platyhelminth nervous system. Parasitology 113: S47–S72.
12. RibeiroP, GearyTG (2010) Neuronal signaling in schistosomes: current status and prospects for postgenomics. Can J Zool 88: 1–22.
13. BerrimanM, HaasBJ, LoVerdePT, WilsonRA, DillonGP, et al. (2009) The genome of the blood fluke Schistosoma mansoni. Nature 460: 352–358.
14. ProtasioAV, TsaiIJ, BabbageA, NicholS, HuntM, et al. (2012) A systematically improved high quality genome and transcriptome of the human blood fluke Schistosoma mansoni. PLoS Negl Trop Dis 6: e1455.
15. BentleyGN, JonesAK, Oliveros ParraWG, AgnewA (2004) ShAR1α and ShAR1β: novel putative nicotinic acetylcholine receptor subunits from the platyhelminth blood fluke Schistosoma. Gene 329: 27–38.
16. BentleyGN, JonesAK, AgnewA (2007) ShAR2β, a divergent nicotinic acetylcholine receptor subunit from the blood fluke Schistosoma. Parasitology 134: 833–840.
17. Raymond-DelpechV, MatsudaK, SattelleB, RauhJ, SattelleD (2005) Ion channels: molecular targets of neuroactive insecticides. Invert Neurosci 5: 119–133.
18. van den EndenE (2009) Pharmacotherapy of helminth infection. Expert Opin Pharmacother 10: 435–451.
19. OlsenRW, SieghartW (2009) GABAA receptors: subtypes provide diversity of function and pharmacology. Neuropharmacology 56: 141–148.
20. LummisSCR (2012) 5-HT3 receptors. J Biol Chem 287: 40239–40245.
21. ThompsonAJ, LesterHA, LummisSC (2010) The structural basis of function in Cys-loop receptors. Q Rev Biophys 43: 449–499.
22. DentJ (2006) Evidence for a diverse cys-loop ligand-gated ion channel superfamily in early bilateria. J Mol Evol 62: 523–535.
23. GalziJ-L, Devillers-ThieryA, HussyN, BertrandS, ChangeuxJ-P, et al. (1992) Mutations in the channel domain of a neuronal nicotinic receptor convert ion selectivity from cationic to anionic. Nature 359: 500–505.
24. KeramidasA, MoorhouseAJ, FrenchCR, SchofieldPR, BarryPH (2000) M2 pore mutations convert the glycine receptor channel from being anion- to cation-selective. Biophys J 79: 247–259.
25. BaardAP, SommersDK, HoniballPJ, FourieED, du ToitLE (1979) Preliminary results in human schistosomiasis with Ro 11-3128. S Afr Med J 55: 617–618.
26. O'BoyleC, LambeR, DarraghA (1985) Central effects in man of the novel schistosomicidal benzodiazepine meclonazepam. Eur J Clin Pharmacol 29: 105–108.
27. MahajanA, KumarV, MansourN, BickleQ, ChibaleK (2008) Meclonazepam analogues as potential new antihelmintic agents. Bioorg Med Chem Lett 18: 2333–2336.
28. MenezesCMS, RiveraG, AlvesMA, do AmaralDN, ThibautJPB, et al. (2012) Synthesis, biological evaluation, and structure–activity relationship of clonazepam, meclonazepam, and 1,4-benzodiazepine compounds with schistosomicidal activity. Chem Biol Drug Des 79: 943–949.
29. LynaghT, LynchJW (2012) Ivermectin binding sites in human and invertebrate Cys-loop receptors. Trends Pharmacol Sci 33: 432–441.
30. CromerBA, MortonCJ, ParkerMW (2002) Anxiety over GABAA receptor structure relieved by AChBP. Trends Biochem Sci 27: 280–287.
31. SchultzJ, MilpetzF, BorkP, PontingCP (1998) SMART, a simple modular architecture research tool: Identification of signaling domains. Proc Natl Acad Sci U S A 95: 5857–5864.
32. QuevillonE, SilventoinenV, PillaiS, HarteN, MulderN, et al. (2005) InterProScan: protein domains identifier. Nucleic Acids Res 33: W116–120.
33. de CastroE, SigristCJA, GattikerA, BulliardV, Langendijk-GenevauxPS, et al. (2006) ScanProsite: detection of PROSITE signature matches and ProRule-associated functional and structural residues in proteins. Nucleic Acids Res 34: W362–W365.
34. ArenaJP, LiuKK, ParessPS, FrazierEG, CullyDF, et al. (1995) The mechanism of action of avermectins in Caenorhabditis elegans: correlation between activation of glutamate-sensitive chloride current, membrane binding, and biological activity. J Parasitol 81: 286–294.
35. ShoopWL, OstlindDA, RohrerSP, MickleG, HainesHW, et al. (1995) Avermectins and milbemycins against Fasciola hepatica: in vivo drug efficacy and in vitro receptor binding. Int J Parasitol 25: 923–927.
36. CampbellWC, FisherMH, StapleyEO, Albers-SchönbergG, JacobTA (1983) Ivermectin: a potent new antiparasitic agent. Science 221: 823–828.
37. HibbsRE, GouauxE (2011) Principles of activation and permeation in an anion-selective Cys-loop receptor. Nature 474: 54–60.
38. LynaghT, WebbTI, DixonCL, CromerBA, LynchJW (2011) Molecular determinants of Ivermectin sensitivity at the glycine receptor chloride channel. J Biol Chem 286: 43913–43924.
39. LynaghT, LynchJW (2010) A glycine residue essential for high ivermectin sensitivity in Cys-loop ion channel receptors. Int J Parasitol 40: 1477–1481.
40. YoungND, JexAR, LiB, LiuS, YangL, et al. (2012) Whole-genome sequence of Schistosoma haematobium. Nat Genet 44: 221–225.
41. WangX, ChenW, HuangY, SunJ, MenJ, et al. (2011) The draft genome of the carcinogenic human liver fluke Clonorchis sinensis. Genome Biol 12: R107.
42. OlsonPD, ZarowieckiM, KissF, BrehmK (2012) Cestode genomics – progress and prospects for advancing basic and applied aspects of flatworm biology. Parasite Immunol 34: 130–150.
43. LynchJW (2004) Molecular structure and function of the glycine receptor chloride channel. Physiol Rev 84: 1051–1095.
44. YoungAB, ZukinSR, SnyderSH (1974) Interaction of benzodiazepines with central nervous glycine receptors: possible mechanism of action. Proc Natl Acad Sci U S A 71: 2246–2250.
45. ThioLL, ShanmugamA, IsenbergK, YamadaK (2003) Benzodiazepines block α2-containing inhibitory glycine receptors in embryonic mouse hippocampal neurons. J Neurophysiol 90: 89–99.
46. LeesG, BeadleDJ, NeumannR, BensonJA (1987) Responses to GABA by isolated insect neuronal somata: pharmacology and modulation by a benzodiazepine and a barbiturate. Brain Res 401: 267–278.
47. SattelleDB, LummisSCR, WongJFH, RauhJJ (1991) Pharmacology of insect GABA receptors. Neurochem Res 16: 363–374.
48. BuckinghamSD, MatsudaK, HosieAM, BaylisHA, SquireMD, et al. (1996) Wild-type and insecticide-resistant homo-oligomeric GABA receptors of Drosophila melanogaster stably expressed in a Drosophila cell line. Neuropharmacology 35: 1393–1401.
49. HosieAM, SattelleDB (1996) Allosteric modulation of an expressed homo-oligomeric GABA-gated chloride channel of Drosophila melanogaster. Br J Pharmacol 117: 1229–1237.
50. BennettJL (1980) Characteristics of antischistosomal benzodiazepine binding sites in Schistosoma mansoni. J Parasitol 66: 742–747.
51. NoëlF, Mendonça-SilvaDL, ThibautJ-PB, LopesDVS (2007) Characterization of two classes of benzodiazepine binding sites in Schistosoma mansoni. Parasitology 134: 1003–1012.
52. GavishM, BachmanI, ShoukrunR, KatzY, VeenmanL, et al. (1999) Enigma of the peripheral benzodiazepine receptor. Pharmacol Rev 51: 629–650.
53. ThibautJPB, MonteiroLM, LeiteLCC, MenezesCMS, LimaLM, et al. (2009) The effects of 3-methylclonazepam on Schistosoma mansoni musculature are not mediated by benzodiazepine receptors. Eur J Pharmacol 606: 9–16.
54. PaxR, BennettJL, FettererR (1978) A benzodiazepine derivative and praziquantel: effects on musculature of Schistosoma mansoni and Schistosoma japonicum. Naunyn Schmiedebergs Arch Pharmacol 304: 309–315.
55. ForgetE, FélixH, NotteghemMJ, LégerN (1982) Appréciation de l'activité de dérivés schistosomicides en microscopie électronique. Bull Soc Pathol Exot Filiales 75: 192–200.
56. BrickerCS, DepenbuschJW, BennettJL, ThompsonDP (1983) The relationship between tegumental disruption and muscle contraction Schistosoma mansoni exposed to various compounds. Z Parasitenkd 69: 61–71.
57. Pica-MattocciaL, RuppelA, XiaCM, CioliD (2008) Praziquantel and the benzodiazepine Ro 11-3128 do not compete for the same binding sites in schistosomes. Parasitology 135: 47–54.
58. RampeD, TriggleDJ (1986) Benzodiazepines and calcium channel function. Trends Pharmacol Sci 7: 461–464.
59. BerlinA, DahlströmH (1975) Pharmacokinetics of the anticonvulsant drug clonazepam evaluated from single oral and intravenous doses and by repeated oral administration. Eur J Clin Pharmacol 9: 155–159.
60. CrevoisierC, DelisleMC, JosephI, FolettiG (2003) Comparative single-dose pharmacokinetics of clonazepam following intravenous, intramuscular and oral administration to healthy volunteers. Eur Neurol 49: 173–177.
61. CoassoloP, AubertC, CanoJ-P (1985) Plasma determination of 3-methylclonazepam by capillary gas chromatography. J Chromatogr 338: 347–355.
62. CullyDF, ParessPS, LiuKK, SchaefferJM, ArenaJP (1996) Identification of a Drosophila melanogaster glutamate-gated chloride channel sensitive to the antiparasitic agent avermectin. J Biol Chem 271: 20187–20191.
63. KaneNS, HirschbergB, QianS, HuntD, ThomasB, et al. (2000) Drug-resistant Drosophila indicate glutamate-gated chloride channels are targets for the antiparasitics nodulisporic acid and ivermectin. Proc Natl Acad Sci U S A 97: 13949–13954.
64. EguchiY, IharaM, OchiE, ShibataY, MatsudaK, et al. (2006) Functional characterization of Musca glutamate- and GABA-gated chloride channels expressed independently and coexpressed in Xenopus oocytes. Insect Mol Biol 15: 773–783.
65. HirataK, IshidaC, EguchiY, SakaiK, OzoeF, et al. (2008) Role of a serine residue (S278) in the pore-facing region of the housefly L-glutamate-gated chloride channel in determining sensitivity to noncompetitive antagonists. Insect Mol Biol 17: 341–350.
66. ForresterSG, PrichardRK, DentJA, BeechRN (2003) Haemonchus contortus: HcGluCla expressed in Xenopus oocytes forms a glutamate-gated ion channel that is activated by ibotenate and the antiparasitic drug ivermectin. Mol Biochem Parasitol 129: 115–121.
67. McCaveraS, RogersAT, YatesDM, WoodsDJ, WolstenholmeAJ (2009) An Ivermectin-sensitive glutamate-gated chloride channel from the parasitic nematode Hemonchus contortus. Mol Pharmacol 75: 1347–1355.
68. NjueAI, HayashiJ, KinneL, FengX-P, PrichardRK (2004) Mutation in the extracellular domains of glutamate-gated chloride channel a3 and b subunits from ivermectin-resistant Cooperia oncophora affect agonist sensitivity. J Neurochem 89: 1137–1147.
69. KehoeJ, BuldakovaS, AcherF, DentJ, BregestovskiP, et al. (2009) Aplysia cys-loop glutamate-gated chloride channels reveal convergent evolution of ligand specificity. J Mol Evol 69: 125–141.
70. CullyDF, VassilatisDK, LiuKK, ParessPS, Van der PloegLHT, et al. (1994) Cloning of an avermectin-sensitive glutamate-gated chloride channel from Caenorhabditis elegans. Nature 371: 707–711.
71. SmithMM, WarrenVA, ThomasBS, BrochuRM, ErtelEA, et al. (2000) Nodulisporic acid opens insect glutamate-gated chloride channels: identification of a new high affinity modulator. Biochemistry 39: 5543–5554.
72. ZerlotiniA, HeigesM, WangH, MoraesRLV, DominitiniAJ, et al. (2009) SchistoDB: a Schistosoma mansoni genome resource. Nucl Acids Res 37: D579–582.
73. Consortium TSjGSaFA (2009) The Schistosoma japonicum genome reveals features of host-parasite interplay. Nature 460: 345–351.
74. ZhangG, FangX, GuoX, LiL, LuoR, et al. (2012) The oyster genome reveals stress adaptation and complexity of shell formation. Nature 490: 49–54.
75. Mendonça-SilvaDL, GardinoPF, KubruslyRCC, De MelloFG, NoëlF (2004) Characterization of a GABAergic neurotransmission in adult Schistosoma mansoni. Parasitology 129: 137–146.
76. ErikssonKS, MauleAG, HaltonDW, PanulaPAJ, ShawC (1995) GABA in the nervous system of parasitic flatworms. Parasitology 110: 339–346.
77. ErikssonKS, PanulaP (1994) Gamma-aminobutyric acid in the nervous system of a planarian. J Comp Neurol 345: 528–536.
78. ErikssonK, PanulaP, ReuterM (1995) GABA in the nervous system of the planarian Polycelis nigra. Hydrobiologia 305: 285–289.
79. Verjovski-AlmeidaS, DeMarcoR, MartinsEA, GuimaraesPE, OjopiEP, et al. (2003) Transcriptome analysis of the acoelomate human parasite Schistosoma mansoni. Nat Genet 35: 148–157.
80. Solis-sotoJM, De Jong BrinkM (1994) Immunocytochemical study on biologically active neurosubstances in daughter sporocysts and cercariae of Trichobilharzia ocellata and Schistosoma mansoni. Parasitology 108: 301–311.
81. MellinTN, BuschRD, WangCC, KathG (1983) Neuropharmacology of the parasitic trematode, Schistosoma Mansoni. Am J Trop Med Hyg 32: 83–93.
82. MillerCL, DayTA, BennettJL, PaxRA (1996) Schistosoma mansoni: glutamate-induced contractions in isolated muscle fibers; evidence for a glutamate transporter. Exp Parasitol 84: 410–419.
83. Mendonça-SilvaDL, PessôaRF, NoëlF (2002) Evidence for the presence of glutamatergic receptors in adult Schistosoma mansoni. Biochem Pharmacol 64: 1337–1344.
84. TaftAS, NoranteFA, YoshinoTP (2010) The identification of inhibitors of Schistosoma mansoni miracidial transformation by incorporating a medium-throughput small-molecule screen. Exp Parasitol 125: 84–94.
85. TamanA, RibeiroP (2011) Glutamate-mediated signaling in Schistosoma mansoni: a novel glutamate receptor is expressed in neurons and the female reproductive tract. Mol Biochem Parasitol 176: 42–50.
86. TamanA, RibeiroP (2011) Characterization of a truncated metabotropic glutamate receptor in a primitive metazoan, the parasitic flatworm Schistosoma mansoni. PLoS ONE 6: e27119.
87. ZamanianM, KimberM, McVeighP, CarlsonS, MauleA, et al. (2011) The repertoire of G protein-coupled receptors in the human parasite Schistosoma mansoni and the model organism Schmidtea mediterranea.. BMC Genomics 12: 596.
88. El-ShehabiF, VermeireJJ, YoshinoTP, RibeiroP (2009) Developmental expression analysis and immunolocalization of a biogenic amine receptor in Schistosoma mansoni. Exp Parasitol 122: 17–27.
89. Goldin AL, Bernardo R (1992) Maintenance of Xenopus laevis and oocyte injection. Meth in Enzymol: Academic Press. pp. 266–279.
90. PutrenkoI, ZakikhaniM, DentJA (2005) A family of acetylcholine-gated chloride channel subunits in Caenorhabditis elegans. J Biol Chem 280: 6392–6398.
91. PeiJ, KimB-H, GrishinNV (2008) PROMALS3D: a tool for multiple protein sequence and structure alignments. Nucleic Acids Res 36: 2295–2300.
92. WaterhouseAM, ProcterJB, MartinDMA, ClampM, BartonGJ (2009) Jalview Version 2 — a multiple sequence alignment editor and analysis workbench. Bioinformatics 25: 1189–1191.
93. DereeperA, GuignonV, BlancG, AudicS, BuffetS, et al. (2008) Phylogeny.fr: robust phylogenetic analysis for the non-specialist. Nucleic Acids Res 36: W465–W469.
94. DereeperA, AudicS, ClaverieJ-M, BlancG (2010) BLAST-EXPLORER helps you building datasets for phylogenetic analysis. BMC Evol Biol 10: 1–6.
95. GuindonS, GascuelO (2003) A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol 52: 696–704.
96. AnisimovaM, GascuelO (2006) Approximate likelihood-ratio test for branches: a fast, accurate, and powerful alternative. Syst Biol 55: 539–552.
97. ChevenetF, BrunC, BañulsA-L, JacqB, ChristenR (2006) TreeDyn: towards dynamic graphics and annotations for analyses of trees. BMC Bioinformatics 7: 1–9.
98. SteinL, SternbergP, DurbinR, Thierry-MiegJ, SpiethJ (2001) WormBase: network access to the genome and biology of Caenorhabditis elegans.. Nucleic Acids Res 29: 82–86.
99. McQuiltonP, St. PierreSE, ThurmondJ, ConsortiumtF (2012) FlyBase 101 – the basics of navigating FlyBase. Nucleic Acids Res 40: D706–D714.
Štítky
Hygiena a epidemiológia Infekčné lekárstvo LaboratóriumČlánok vyšiel v časopise
PLOS Pathogens
2013 Číslo 8
- Parazitičtí červi v terapii Crohnovy choroby a dalších zánětlivých autoimunitních onemocnění
- Očkování proti virové hemoragické horečce Ebola experimentální vakcínou rVSVDG-ZEBOV-GP
- Koronavirus hýbe světem: Víte jak se chránit a jak postupovat v případě podezření?
Najčítanejšie v tomto čísle
- Host Immune Response to Intestinal Amebiasis
- Bed Bugs and Infectious Disease: A Case for the Arboviruses
- Discovery of Anthelmintic Drug Targets and Drugs Using Chokepoints in Nematode Metabolic Pathways
- Relevance of Trehalose in Pathogenicity: Some General Rules, Yet Many Exceptions