# Epigenetic control of gene function in schistosomes: a source of therapeutic targets?

Alejandro Cabezas Cruz, Julien Lancelot, Stéphanie Caby, Guilherme Oliveira, Raymond J. Pierce

► To cite this version:

Alejandro Cabezas Cruz, Julien Lancelot, Stéphanie Caby, Guilherme Oliveira, Raymond J. Pierce. Epigenetic control of gene function in schistosomes: a source of therapeutic targets?. Frontiers in Genetics, 2014, 5 (317). hal-02630213

HAL Id: hal-02630213

https://hal.inrae.fr/hal-02630213

Submitted on 27 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Epigenetic control of gene function in schistosomes: a source of therapeutic targets?

Alejandro Cabezas-Cruz<fnref n="1" />, Julien Lancelot<fnref n="1" />, Stéphanie Caby<fnref n="1" />, Guilherme Oliveira<fnref n="2" /> and Raymond J. Pierce<fnref n="1" />*

**Edited by:**

Paul J. Brindley, The George Washington University, USA

**Reviewed by:**

Amparo Latorre, University of Valencia, Spain
Abd El-Latif Hesham, Assiut University, Egypt

**\*Correspondence:**

Raymond J. Pierce, INSERM U1019 – CNRS UMR 8204, Center for Infection and Immunity of Lille, Université de Lille, Institut Pasteur de Lille, 1 rue du Professeur A. Calmette, 59019 Lille, France
e-mail: raymond.pierce@pasteur-lille.fr

The discovery of the epigenetic regulation of gene expression has revolutionized both our understanding of how genomes function and approaches to the therapy of numerous pathologies. Schistosomes are metazoan parasites and as such utilize most, if not all the epigenetic mechanisms in play in their vertebrate hosts: histone variants, histone tail modifications, non-coding RNA and, perhaps, DNA methylation. Moreover, we are acquiring an increasing understanding of the ways in which these mechanisms come into play during the complex schistosome developmental program. In turn, interest in the actors involved in epigenetic mechanisms, particularly the enzymes that carry out epigenetic modifications of histones or nucleic acid, as therapeutic targets has been stimulated by the finding that their inhibitors exert profound effects, not only on survival, but also on the reproductive function of *Schistosoma mansoni*. Here, we review our current knowledge, and what we can infer, about the role of epigenetic mechanisms in schistosome development, differentiation and survival. We will consider which epigenetic actors can be targeted for drug discovery and what strategies can be employed to develop potent, selective inhibitors as drugs to cure schistosomiasis.

**Keywords:** epigenetics, schistosome, drug discovery, histone modifying enzymes, DNA methylation, microRNAs

## INTRODUCTION

Schistosomiasis is caused by flatworm parasites of the genus *Schistosoma*, five species of which infect humans in 74 tropical and sub-tropical countries. It is estimated that more than 230 million people are infected, of which 90% are in sub-Saharan Africa (Colley et al., 2014 for review). Nearly 30 years after its introduction the treatment and control of schistosomiasis relies almost exclusively on praziquantel, the only drug effective against all schistosome species infecting humans. Its use has been and remains an unquestionable success story; mass treatments of school-age children in sub-Saharan Africa under the Schistosomiasis Control Initiative (Fenwick et al., 2009) hold the promise of a marked continent-wide reduction in disease morbidity and mortality. Nevertheless, the massive use of this drug may well-lead to the selection of resistant/tolerant parasite strains. Episodes of drug tolerance have been reported (Doenhoff et al., 2008; Melman et al., 2009) and can be induced readily in the laboratory (Fallon and Doenhoff, 1994). In addition, changes to the local genetic polymorphism of parasites following treatment of the population have been detected (Norton et al., 2010; Coeli et al., 2013), suggesting an effect of drug selection pressure. Although the precise mechanism of action of praziquantel in killing schistosomes is unknown, its initial effects include the rapid influx of $Ca^{2+}$ ions and calcium-dependent muscle contraction and paralysis (Day et al., 1992) and this may be mediated via its interaction with a calcium channel beta subunit

(Kohn et al., 2001). However, resistance to praziquantel may be mediated by increased expression of the P-glycoprotein efflux pump, which is often involved in drug resistance mechanisms, following exposure to praziquantel (Messerli et al., 2009). Whether or not such reports are the harbingers of the development of resistance by schistosomes toward praziquantel, reliance on a single drug is patently untenable in the medium to long term.

Most of the current efforts to identify new drug leads for schistosomiasis and other neglected parasitic diseases rely on the screening of random compound libraries directly on the parasite maintained in culture (phenotypic screening). The recent publication of the genome sequences of a variety of parasites including the three main species of schistosomes that infect humans (Berriman et al., 2009; Zhou et al., 2009; Young et al., 2012) now means that approaches targeting specific gene products or pathways can be envisaged. These can include enzymes with activities specific to the parasite, or at least not found in the human host (e.g., Sayed et al., 2008), metabolic bottlenecks, or molecules that are targeted in other pathologies. For these a wide knowledge base and extensive libraries of inhibitors may already exist that can be exploited as starting points for the development of parasite-selective compounds. This type of approach also has the advantage that the molecular mechanism of action of a given compound, which is required for any new drug, is much easier to determine than with the random screening approach. However, both strategies are still used and both have proved fruitful sources

<footnote n="1">¹ Institut National de la Santé et de la Recherche Médicale U1019 – Centre National de la Recherche Scientifique UMR 8204, Center for Infection and Immunity of Lille, Institut Pasteur de Lille, Université de Lille, Lille, France</footnote>

<footnote n="2">² Genomics and Computational Biology Group, Fundação Oswaldo Cruz, Center for Excellence in Bioinformatics, Centro de Pesquisas René Rachou, National Institute of Science and Technology in Tropical Diseases, Belo Horizonte, Brazil</footnote>
of new drugs (Swinney and Anthony, 2011) although a more recent analysis of discovery of first-in-class drugs suggests a growing predominance of target-based approaches (Eder et al., 2014). Drug discovery is not a zero-sum game, but more a Nash equilibrium (Nash, 1950; Holt and Roth, 2004) in which the coexistence of strategies is not only possible but can also be highly productive.

Schistosomes are digenean parasites that successively infect freshwater snails (the intermediate host) and the vertebrate definitive host. They reproduce both asexually (within the snail host) and sexually (vertebrate host) and their life-cycle includes four distinct morphological forms and separate sexes at the adult worm stage (Colley et al., 2014). The complexity of schistosome development and differentiation implies a tight control of gene transcription at all stages of the life-cycle and that epigenetic mechanisms are likely to play a crucial role in these processes, suggesting that they are viable drug targets. In other pathologies, but most intensively in cancer, the targeting of epigenetic processes is increasingly exploited. Indeed, two histone deacetylase (HDAC) inhibitors have already been approved for use and a number of other candidate drugs are undergoing clinical trials (Arrowsmith et al., 2012). Moreover, large libraries of compounds that affect epigenetic actors are available for testing against parasites. Here, we will consider which epigenetic mechanisms can be targeted in schistosomes and what methodologies can be used to develop parasite-selective drug leads.

## EPIGENETIC MECHANISMS AS DRUG TARGETS

The term “epigenetics” envelops a variety of heritable changes in gene expression that are linked to structural modifications of the chromatin, without changes to the DNA sequence. These include DNA methylation, reversible post-translational modifications of histones, histone variants, chromatin remodeling factors and non-coding RNAs. Viewed as potential targets, the most readily “druggable” are the enzymes that carry out DNA methylation and histone modifications, and increasingly, micro-RNAs (miRNAs) among the non-coding RNA categories (Figure 1).

The investigation of the role of epigenetic mechanisms in the control of gene transcription in schistosomes, and hence in biological processes like development and reproduction, is in its early stages. Nevertheless, the knowledge so far acquired, or inferred from the nature of schistosomes as invertebrate metazoan organisms and from a detailed analysis of the epigenetic actors encoded in their genomes, can be exploited to develop novel therapeutic strategies. Moreover, insights into schistosome epigenetic mechanisms has been gained from studies aimed at developing such strategies, including for example the characterization of the actions of inhibitors of histone modifying enzymes (HMEs), or from transcript knockdown studies. Here we will review the current state of knowledge of the epigenetic apparatus in schistosomes, including the still disputed significance of DNA methylation, the miRNA repertoire, the histone modifying enzyme complement and the potential for the development of novel drug treatments targeting these elements.

## DNA METHYLATION

DNA methylation encompasses the methylation or hydroxymethylation of cytosine residues, mainly, but not exclusively

within CpG dinucleotides (Baubec and Schübeler, 2014) and is an important epigenetic mark associated with gene repression. In vertebrates three DNA methyltransferases (Dnmts) establish (Dnmt3a and Dnmt3b) and maintain (Dnmt1) DNA methylation marks. A further member of this family, Dnmt2 is primarily a tRNA methyltransferase with only weak DNA methyltransferase activity (Schaefer and Lyko, 2010). Disruption of DNA methylation patterns is present in a variety of diseases, particularly in cancer in which many oncogenic pathways lead to Dnmt1 overexpression, an overall DNA hypomethylation concomitant with hypermethylation of tumor suppressor genes at CpG islands in the promoter regions. Agents that provoke DNA demethylation, such as 5-azacytidine and 5-aza-2'-deoxycytidine (decitabine) have been approved for use in myelodysplastic syndrome (Yoo and Jones, 2006), but their mode of action is complex. It involves conversion to a triphosphate metabolite, incorporation into DNA provoking a DNA damage response and covalent trapping of Dnmt isoforms, followed by proteolysis of the Dnmts, demethylation and the reactivation of the hypermethylated genes (Streseman and Lyko, 2008). The cytotoxicity of the 5-aza nucleosides and the lack of a direct inhibitory effect on Dnmts have led to the search for leads for new drug development. One example is laccaic acid A, a recently developed direct DNA-competitive inhibitor of Dnmt1 (Fagan et al., 2013).

## DNA METHYLATION IN SCHISTOSOMES

The presence of functional DNA methylation marks in schistosome genomes is controversial. Early work in which Southern blot analysis was carried out for selected genes after digestion with methylcytosine tolerant or sensitive restriction enzyme isoschizomers (*HpaII* and *MspI*) showed no differences in the restriction profiles for adult male or female *S. mansoni* DNA (*Fantappié* et al., 2001). Moreover, the methylcytosine-dependent restriction endonuclease McrBC failed to digest *S. mansoni* DNA. However, a more recent study (Geyer et al., 2011) in which a variety of more sensitive methods including GC-MS, antimethylcytosine antibodies and targeted bisulfite sequencing were used concluded that cytosine methylation was indeed present and a hypermethylated repetitive intron within a forkhead gene was characterized. The only DNA methyltransferase encoded in the schistosome genome is Dnmt2. The methyltransferases usually associated with DNA methylation; Dnmt1 and Dnmt3 orthologs are both absent. Dnmt2 has only weak DNA methyltransferase activity but has robust methyltransferase activity toward tRNAAsp and other tRNAs (Goll et al., 2006). The diverse group of animal species (including *S. mansoni* and *Drosophila melanogaster*) that express only Dnmt2 have very low DNA methylation levels (Kraus and Reuter, 2011). However, Dnmt2 does retain some cytosine methyltransferase activity (Hermann et al., 2003) and Geyer et al. (2011) showed that siRNA knockdown of *SmDnmt2* transcripts reduced overall methylcytosine levels in the schistosome genome. These authors have further suggested that cytosine methylation is conserved throughout the phylum Platyhelminthes (Geyer et al., 2013). Against this, a comprehensive study (Raddatz et al., 2013) using whole-genome bisulfite sequencing showed that the *S. mansoni* genome lacked a detectable DNA methylation pattern, even at the “hypermethylated” locus identified by Geyer et al. (2011).
<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Schematic representation of the major druggable epigenetic mechanisms. Histone modifications shown are limited to acetylation and methylation since the enzymes (writers and erasers)</figcaption></figure>

and recognition domains (readers: bromodomains) are the most studied for drug development. Also shown are DNA methylation and microRNAs.

Some clusters of incompletely converted cytosines were detected outside this region, but were consistent with bisulfite deamination artifacts (Warnecke et al., 2002). However, although these results strongly suggested that the S. mansoni genome is in fact unmethylated, the criticism has been leveled that the life-cycle stage analyzed, adult male worms, has the lowest level of DNA methylation measured using an ELISA method (Geyer et al., 2013). Notwithstanding this controversy, which will only be resolved by genome-wide bisulfite sequencing of other life-cycle stages, Dnmt inhibitors were found to strongly affect adult worms, particularly in terms of the morphology of the ovaries and *in vitro* egg-laying (Geyer et al., 2011). Whether or not this is due to the inhibition of DNA or tRNA methylation, it does suggest that Dnmt inhibitors such as 5-azacytidine may provide the basis for developing precursors of novel anti-schistosome drugs.

## MICRO-RNAs

Non-coding (nc) RNAs include many different classes of transcripts that do not code for proteins, but have various regulatory

roles in transcription, stability or translation of protein-coding genes. Of these, miRNAs are the best characterized in terms of their functional roles and pathological implications, as well as therapeutic strategies targeting them (Ling et al., 2013). They are generated from long, capped and polyadenylated transcripts that are processed by a nuclear complex containing RNase III (Drosha: canonical pathway) or by the mRNA splicing machinery (non-canonical pathway) (Li and Rana, 2014 for review) into 60-100 nucleotide precursors that are then transported into the cytoplasm where they are processed by the RNase Dicer into mature, double stranded miRNAs (Figure 1). Classically, miRNAs regulate transcript levels through binding to the 3'UTR regions of their target mRNAs, usually resulting in translational inhibition or mRNA destruction. However, it is now clear that miRNAs may have other mechanisms of action, for instance increasing translation via the recruitment of protein complexes to the mRNA or by binding proteins that block translation (Elring et al., 2010). Different miRNAs have been shown to have either tumor-suppressive (e.g., miR-15a-miR-16-1 cluster) or oncogenic (mIR-21, mIR-17-mIR-92 cluster, mIR-155) properties.
Indeed, miRNAs can drive cancer: miR 155 overexpression on its own provokes lymphoblastic leukemia or lymphoma in transgenic mice (Costinean et al., 2006). In cancer therapy, the upregulation of tumor-suppressive miRNAs has the advantage of simultaneously affecting a number of coding or non-coding genes that are targeted by the miRNA and that may be involved in the same or interacting pathways. A disadvantage is that a given miRNA may have different or even opposite effects in different cell types, depending on the expression patterns of its target genes. However, such considerations would be less of an obstacle in the therapy of parasitic diseases where it can be assumed that any disruption, positive or negative, of miRNA effects would be potentially deleterious to the parasite. Most current therapeutic strategies targeting miRNA in cancer are aimed at downregulating or blocking the function of oncogenic miRNA. One example consists in the use of antisense oligonucleotides, particularly those containing locked nucleic acids (LNA anti-miRs) which are bicyclic RNA analogs in a locked configuration. One such compound, an anti-viral, miravirsen, is in clinical trials for the treatment of hepatitis C viral infection (Janssen et al., 2013; Lieberman and Sarnow, 2013).

## SCHISTOSOME miRNAs

A survey of the available S. mansoni EST sequences (Oliveira et al., 2011) concluded that 10.3% (21,107 sequences) match the genome but have no protein coding potential and are therefore possible ncRNAs. This in turn suggests that the parasite may use a range of ncRNAs in transcriptional and translational regulation. Moreover, the presence of proteins involved in miRNA processing (Drosha, Dicer, and Argonaute) (reviewed in Oliveira et al., 2011) supports a role for miRNA regulation of development and differentiation of schistosomes and explains the effectiveness of RNAi and siRNA knockdown of transcription in the parasite (Boyle et al., 2003). miRNAs were first identified in S. japonicum (Xue et al., 2008; Huang et al., 2009) in two separate studies that demonstrated the existence of a limited number of miRNA that are conserved in other organisms including humans and several hundred novel miRNAs. In the second study all the novel miRNAs (172) were identified by an inferred RNA hairpin and many were differentially expressed during the life-cycle (Huang et al., 2009). Deep sequencing studies (Hao et al., 2010; Wang et al., 2010) further showed the presence of large numbers of siRNAs derived from transposable elements, but also identified 38 novel S. japonicum miRNAs. In S. mansoni, the sequencing of a small-RNA cDNA library yielded 211 novel miRNA candidates of which 11 were further verified by Northern blotting (Simoes et al., 2011). Therefore, although further work is necessary to validate the schistosome specific miRNAs and determine which of them are shared between schistosome species, it is clear that these miRNAs are potential therapeutic targets. It is of note that flatworms show a gradual loss of conserved miRNAs during evolution (Fromm et al., 2013), which is suggested to be due to morphological simplification. However, they have equally gained specific ncRNAs, including endogenous siRNAs, which are differentially expressed during development, notably during the sexual differentiation of female worms (Cai et al., 2011; Sun et al., 2014). A therapeutic strategy based on LNA anti-miRs would have the advantage of targeting parasite-specific

sequences and hence avoiding off-target effects, but it is not yet known whether individual miRNAs could be valid therapeutic targets. There are several additional challenges associated with such a strategy (Ling et al., 2013), the main ones being to ensure bioavailability to the parasite and oral delivery, which would require a significant effort to investigate appropriate chemical substitutions.

## POST-TRANSLATIONAL MODIFICATIONS OF HISTONES

Histone post-translational modifications are currently under the most intensive study for drug development. The “writers” that add groups to histone N-terminal tails, “erasers” that remove them, or “readers” that recognize and bind them, are all potential therapeutic targets (Figure 1). The increasing variety of possible modifications includes phosphorylation, ubiquitinylation and sumoylation, but acetylation and methylation are the most abundant, most studied and their activity is mediated by the largest number of druggable proteins (Arrowsmith et al., 2012). Gene regulation is effected by combinations of these histone marks, leading to the “histone code” hypothesis (Strahl and Allis, 2000) whereby different chromatin states are defined by specific repertoires of marks.

Histone acetylation is a dynamic process regulated by histone acetyltransferases (HAT) that use acetyl-CoA as a co-factor and transfer an acetyl residue to the e-amino group of lysines, particularly in the N-terminal tails of histones H3 and H4. The HDAC that remove this mark belong to four classes in mammals. Classes I, II, and IV have structurally-related catalytic domains and a Zn²⁺-dependent catalytic mechanism (Gregoretti et al., 2004). The class III HDACs, or sirtuins, are phylogenetically unrelated and rely on NAD⁺ as a co-factor (Greiss and Gartner, 2009). Histone acetylation neutralizes the positive charge of the lysine, leading to a more relaxed structure permitting recruitment of the transcriptional machinery and in consequence is associated with transcriptional activation.

Histone methyl marks are written on lysine or arginine residues in histone tails by S-adenosylmethionine-dependent methyltransferases and erased by two classes of demethylase, the Jumonji family of demethylases that are 2-oxoglutarate-dependent, or the flavin-dependent lysine-specific demethylase 1 (KDM1/LSD1) and 2 (KDM2/LSD2). Unlike acetylation, a methyl group has no effect on the overall charge of the lysine or arginine residue that carries it and the effects of the mark are mediated by “reader” proteins that either compact the nucleosomes or form complexes with other regulatory proteins. Moreover, lysine residues can react with different reader domains depending on their position and degree (mono-, di, -or tri-) of methylation and can consequently integrate signal platforms determining activation or repression of transcription (Badeaux and Shi, 2013).

The current list of inhibitors of “Histone Modifying Enzymes” (HMEs) approved for use in humans or in clinical trials reflects the initial concentration in this field on the HDACs as therapeutic targets. Of 18 such compounds, 17 are HDAC inhibitors (including the two approved compounds, Vorinostat, and Romidepsin) and another is a sirtuin (Sirt) 1 inhibitor (Arrowsmith et al., 2012; West and Johnstone, 2014). Of the HDAC inhibitors in
clinical trials some inhibit class I and II enzymes indiscriminately, whilst others are more selective. Romidepsin preferentially inhibits the class I HDACs 1, 2, 3, and 8, whilst Vorinostat inhibits the class II HDAC6 and HDAC8 only poorly (Arrowsmith et al., 2012). Selectivity for a given HDAC or class may be of therapeutic importance as these enzymes have different targets. HDACs generally deacetylate both histones and other proteins. HDAC6 is not involved in epigenetic signaling at all but deacetylates tubulin and Hsp90 (Hubbert et al., 2002; Kovacs et al., 2005), while the only known HDAC8 substrate is SMC3, a component of the cohesin complex (Deardorff et al., 2012). HDAC inhibitors occupy the hydrophobic tunnel in the enzymes that accommodates the acetyllysine substrate and coordinate the zinc ion at the base of the tunnel, for example with a hydroxamate grouping as for Vorinostat. Selectivity can be based on differences between the make-up and architecture of the tunnel, or on the surface-accessible rim. However, within the cell HDACs are often part of multi-protein complexes that may alter substrate and inhibitor specificities compared to isolated recombinant proteins (Bantscheff et al., 2011) and this remark is likely even more pertinent in the case of other HMEs.

Of the seven mammalian NAD⁺-dependent class III deacetylases, or sirtuins (Sirt) Sirts 1, 2, 3, 6, and 7 have been shown to possess deacetylase activity (Feldman et al., 2012), although Sirt6 is also a fatty acylase (Jiang et al., 2013). The other two sirtuins, 4 and 5 are both predominately mitochondrial. Sirt5 is a demalonylase and desuccinylase (Du et al., 2011) and has recently been shown to regulate a novel lysine modification, glutarylation (Tan et al., 2014). Sirt4 exhibits ADP-ribosyltransferase activity (Haigis et al., 2006). Sirtuin inhibitors have been developed against Sirts 1, 2, and 3 that couple the deacetylation reaction with the cleavage of NAD⁺, liberating free nicotinamide. Inhibitors can bind either to the conserved NAD⁺-binding C-pocket, like nicotinamide itself, or the acetyllysine peptide-binding cleft between the large and small domains of the enzyme (Yuan and Mamorstein, 2012), or both. In addition, since Sirt1 expression has been associated with increased lifespan and memory, allosteric activators such as resveratrol have been explored as therapeutic agents Hubbard and Sinclair, 2014). Only one sirtuin inhibitor, selisistat, is currently in clinical trials, however, for Huntington's disease (Arrowsmith et al., 2012).

HAT generally lack obvious druggable sites and few selective inhibitors are currently available. The HAT catalytic domains have a conserved organization around a central fold where the acetyl-CoA cofactor binds. The peptide substrate binding site in the only solved structure is shallow and solvent accessible, reducing its capacity to be targeted by drugs (Arrowsmith et al., 2012). Among the inhibitors so far described are natural substances that promiscuously bind a variety of targets (Piaz et al., 2011), or isothiazolone covalent modifiers (Ghizzoni et al., 2009). These latter include the more recently developed pyridoisothiazolones that effectively inhibit cancer cell proliferation (Furdas et al., 2011). However, a potent, selective inhibitor of the HAT EP300, C646, has been developed that binds at the cofactor pocket and has pro-apoptotic effects on prostate cancer cells (Bowers et al., 2010). This indicates that at least certain HATs are valid, standalone therapeutic targets, but effective screening may depend on

the reconstitution of multi-protein complexes in which they are active in the cell and which may modulate their enzymatic activity (Arrowsmith et al., 2012).

The protein methyltransferases include both lysine (KMT) and arginine (PRMT) methyltransferases that are phylogenetically unrelated but share the requirement for S-adenosylmethionine as a cofactor and a cofactor binding site adjacent to the channel that binds the peptide substrate (Arrowsmith et al., 2012). Both these sites can be used to generate selective and potent inhibitors for both PRMTs and KMTs (Spannhoff et al., 2009; Dowden et al., 2010). However, in order to screen certain of the latter enzymes, the reconstitution of protein complexes is a prerequisite. The KMT component of the PRC2 transcription repression complex, EZH2, which methylates H3K27, is inactive on its own and minimally requires the presence of at least two members of the complex, EED and SUZ12 (Helin and Dhanak, 2013). High-throughput screening has been carried out on a complex additionally containing AEBP2 and RbAp48 proteins and has yielded highly selective inhibitors such as GSK126, which is a promising lead for the treatment of lymphoma (McCabe et al., 2012). Complex recomposition is not always necessary; the H3K79 methyltransferase DOT1L requires no partner proteins and a highly selective inhibitor, EPZ-5676, has been developed using a structure-guided strategy, with significant activity in a rat xenograft model of MLL-rearranged leukemia (Daigle et al., 2013).

Histone demethylases are under increasing scrutiny as drug targets (Hojfeldt et al., 2013). Two unrelated families of proteins exert demethylase activity, the LSD family and the JMJC domain-containing demethylases. In mammals the LSD (for Lysine-Specific Demethylase) family comprises only two members, but they share an amine oxidase-like domain with several metabolic enzymes in addition to a SWIRM (SW13, RSC8, and Moira) domain that is only found in some chromatin-associated enzymes. These enzymes use FAD as a cofactor that binds to one of two folded subdomains of the amine oxidase-like domain, the other binding the substrate. However, substrate specificity may be regulated through interactions of the SWIRM domain with protein partners (Metzger et al., 2005). The conservation of the amine oxidase-like domain means that some inhibitors of monoamine oxidases, such as tranylcypromine also inhibit LSD1 (Schenk et al., 2012), however, selectivity for the latter can be improved in derivatives. One such inhibitor, ORY-1001, has greatly improved selectivity for LSD1 and is entering clinical trials (Maes et al., 2013). Few selective or potent inhibitors of JMJC (Jumonji C) demethylases have yet been reported. In these enzymes (31 family members in humans) the JMJC domain is catalytic and the enzymatic mechanism involves two cofactors, Fe(II) and 2-oxoglutarate, which are bound to it. Members of this family can demethylate mono-, di- or tri-methylated lysines. Most reported inhibitors contain metal chelating groups (analogous to HDAC inhibitors, but the latter are poor JMJC inhibitors) that compete with the 2-oxoglutarate cofactor. One remarkably selective inhibitor of JMJD3, an H3K27me3-specific demethylase, has been developed using a structure-guided approach (Kruidenier et al., 2012). This inhibitor, GSK-J1, interacts with critical amino
acids involved in the binding of both the 2-oxoglutarate cofactor and the histone peptide substrate, and is competitive with the former but non-competitive with the latter. It is selective for both JMJD3 and Utx (a closely related JMJC demethylase active on the same substrate as JMJD3) but inactive against other members of the JMJC family. Moreover, GSK-J1 inhibits pro-inflammatory functions of human primary macrophages, indicating a possible therapeutic role. Moreover, a related JMJD3 inhibitor, GSK-J4, has recently been shown to have anti-tumorigenic activity against T-cell acute lymphoblastic leukemia (Ntziachristos et al., 2014).

Of the “readers” of epigenetic marks, bromodomain proteins, which read acetylated lysine residues, have so far attracted the most attention as drug targets. Bromodomains are composed of a characteristic antiparallel bundle of four α-helices that binds acetyllysine in a pocket at one extremity. They were first identified in the *Drosophila* gene brahma, hence the name, and in humans 61 bromodomains have been identified in 46 proteins, some containing more than one, and belong to eight distinct families (Hewings et al., 2012). Although the biological roles of most bromodomain proteins remain unknown, the connection of some of them with diseases such as cancer is becoming clearer (Arrowsmith et al., 2012; Filippakopoulos and Knapp, 2014). Small molecule inhibitors have been principally been developed against the BET bromodomain family. A triazolobenzidine (Nicodeme et al., 2010) was initially isolated as an inducer of ApoA1 expression and a derivative with enhanced activity in a reporter gene assay was only subsequently shown to interact with the bromodomains of BET family members BRD2, 3, and 4 using a chemoproteomic approach (Chung et al., 2011). Following an observation that thieodiazepines could bind to BRD4 a structure-based approach led to the development of a novel thieno-triazolo-1'4-diazepine named JQ1 (Filippakopoulos et al., 2011). This molecule binds to the BRD4 acetyllysine site and has reduced activity against other BRDs, but little or none toward other bromodomains. JQ1 induced a differentiation phenotype and growth arrest in a cell line derived from human squamous carcinoma. Since these founder studies, further selective inhibitors of BET bromodomains have been identified and clinical trials are ongoing, notably for JQ1.

## SCHISTOSOME HISTONE MODIFICATIONS

An *in silico* analysis of the schistosome genome predictions and EST libraries (Anderson et al., 2012) showed that 21 of the 29 histone genes predicted in the genome are expressed in *S. mansoni*, the remainder being either unexpressed or having divergent sequences. Importantly, the N-terminal tails of the nucleosomal histones H3 and H4 are highly conserved, suggesting the functional conservation of the histone marks found in mammalian histones. The study of histone marks and their role in schistosome development and differentiation is still in its infancy, but studies involving inhibitors of HDACs and HATs have demonstrated the importance of histone acetylation and the interest of these enzymes as potential therapeutic targets. HDAC inhibitors including Trichostatin A (TSA) blocked the *in vitro* transformation of *S. mansoni* miracidia into primary

sporocysts in a dose-dependent manner (Azzi et al., 2009) and this correlated with an increase in histone H4 acetylation. More recently, the same authors showed that differences in the levels of histone H3K9 acetylation on the promoters of genes encoding polymorphic mucins correlated with their differential expression in parasite strains compatible or incompatible with a given strain of the intermediate host, the freshwater snail *Biomphalaria glabrata*. TSA treatment ablated these strain-specific differences in expression (Perrin et al., 2013).

Histone acetyltransferase inhibition also has developmental consequences in schistosomes, particularly in egg maturation. The schistosome ortholog of the HAT GCN5 has been shown to acetylate H3 and H2A, and in particular H3K14 (de Moraes Maciel et al., 2008) and the CBP/P300 ortholog SmCBP1 primarily acetylates H4 (Bertin et al., 2006; Fantappié et al., 2008). Knockdown of either or both of these HATs in adult schistosomes has been shown to markedly reduce the transcription levels of the major eggshell protein p14 and to affect egg development. Moreover, these effects are reproduced by treating adult worm pairs with an HAT inhibitor, PU139 (Carneiro et al., 2014). After both inhibitor treatment and RNAi to knock down transcripts of the HATs, the phenotypic effects on egg laying and development were correlated with decreased acetylation of H3 and H4, increased methylation at H3K27, a marker of transcriptional repression, on the p14 proximal promoter.

The effects of both HDAC and HAT inhibitors on schistosomes suggest that histone acetylation may be a legitimate therapeutic target and this was further supported by a preliminary study showing that HDAC inhibitors like TSA and valproic acid could induce time and dose-dependent death of schistosomes (adult worms or schistosomula larvae) in culture (Dubois et al., 2009). Parasite death was associated with the induction of apoptosis in schistosomula shown by both TUNEL staining and the activation of the effector caspases 3/7. Once more, the molecular basis of these effects was evidenced by a global increase in histone H3 and H4 acetylation and significantly increased H4 acetylation at the proximal promoters of HDAC target genes, correlated with an increase in their transcription. More recently, similar effects have been observed using inhibitors of the class III HDAC, the sirtuins (Lancelot et al., 2013). Inhibitors of Sirtuins 1 and 2 such as salermide also induce apoptosis and death of schistosomula in culture. Moreover, salermide induces marked morphological alterations to the female worm genital apparatus, the arrest of egg-laying and the separation of worm pairs.

## SCHISTOSOME HISTONE MODIFYING ENZYMES: WHICH ARE THE BEST TARGETS?

Schistosome HDACs and HATs are clearly potential drug targets, as are probably other HMEs, but there are several challenges that could potentially impede drug development:

- These enzymes are evolutionarily conserved, particularly their catalytic domains, and in order to avoid potential side effects, inhibitors that are selective for the schistosome enzyme have to be developed.
<table><caption>Table 1 | Identity and characteristics of *Schistosoma mansoni* histone modifying enzymes.</caption><thead><tr><th>HME type</th><th>Class</th><th>Closest human ortholog</th><th>Size (aa)</th><th>Substrate specificity</th><th>Gene Id<sup>a</sup></th></tr></thead><tbody><tr><td rowspan="10">HDAC</td><td>I</td><td>HDAC1</td><td>517*</td><td></td><td>Smp_005210</td></tr><tr><td>I</td><td>HDAC3</td><td>418*</td><td></td><td>Smp_093280</td></tr><tr><td>I</td><td>HDAC8</td><td>440*</td><td></td><td>Smp_091990</td></tr><tr><td>II</td><td>HDAC4</td><td>291</td><td></td><td>Smp_191310</td></tr><tr><td>II</td><td>HDAC5</td><td>701</td><td></td><td>Smp_069380</td></tr><tr><td>II</td><td>HDAC6</td><td>1132</td><td></td><td>Smp_138770</td></tr><tr><td>III (Sirtuin)</td><td>Sirt1</td><td>568*</td><td>H1 – H3 – H4</td><td>Smp_138640</td></tr><tr><td>III (Sirtuin)</td><td>Sirt2</td><td>337*</td><td>H4K16</td><td>Smp_084140</td></tr><tr><td>III (Sirtuin)</td><td>Sirt5</td><td>305*</td><td></td><td>Smp_055090</td></tr><tr><td>III (Sirtuin)</td><td>Sirt6</td><td>386*</td><td>H3K9 – H3K56</td><td>Smp_134630</td></tr><tr><td>III (Sirtuin)</td><td>Sirt7</td><td>517*</td><td></td><td>Smp_024670</td></tr><tr><td rowspan="10">HAT</td><td>GNAT</td><td>GCN5 (KAT2A)</td><td>899*</td><td>H3K9 – H3K14 – H3K18 H2B</td><td>Smp_070190</td></tr><tr><td>GNAT</td><td>HAT1 (KAT1)</td><td>435</td><td>H4K5 – H4K12</td><td>Smp_178700</td></tr><tr><td>MYST</td><td>Tip60 (KAT5)</td><td>463</td><td>H2AK5 – H3K14 – H4K5 – H4K8 – H4K12 – H4K16</td><td>Smp_053140</td></tr><tr><td>MYST</td><td>MYST1 (KAT8)</td><td>496</td><td>H4K16</td><td>Smp_194520</td></tr><tr><td>MYST</td><td>MYST2 (KAT7)</td><td>400</td><td>H4K5 – H4K8 – H4K12 – H3</td><td>Smp_171700</td></tr><tr><td>MYST</td><td>MYST3 (KAT6A)</td><td>971</td><td>H3K14</td><td>Smp_131320</td></tr><tr><td>CBP/p300</td><td>CBP/SmCBP1 (KAT3A)</td><td>2093*</td><td>H2AK5 – H2BK15 – H3K14 – H3K18 – H4K5 – H4K8</td><td>Smp_105910</td></tr><tr><td>CBP/p300</td><td>CBP/SmCBP2 (KAT3A)</td><td>1892</td><td>H2AK5 – H2BK15 – H3K14 – H3K18 – H4K5 – H4K8</td><td>Smp_127010</td></tr><tr><td>TAFII250</td><td>TFIID subunit 1</td><td>2241</td><td>H3 – H4</td><td>Smp_166840</td></tr><tr><td rowspan="20">HMT</td><td>SET</td><td>EZH1</td><td>1026</td><td>H3K27</td><td>Smp_078900</td></tr><tr><td>SET</td><td>MLL3 (KMT2C)</td><td>399</td><td>H3K4</td><td>Smp_070210</td></tr><tr><td>SET</td><td>MLL3 (KMT2C)</td><td>1560</td><td>H3K4</td><td>Smp_138030</td></tr><tr><td>SET</td><td>MLL1/4 (KMT2D)</td><td>3002</td><td>H3K4</td><td>Smp_144180</td></tr><tr><td>SET</td><td>MLL5 (KMT2E)</td><td>751</td><td>H3K4</td><td>Smp_161010</td></tr><tr><td>SET</td><td>C20orf11/MLL5/Ranbp9</td><td>1305</td><td></td><td>Smp_009980</td></tr><tr><td>SET</td><td>NSD2/WHSC1</td><td>1746</td><td>H3K4 – H4K20</td><td>Smp_160700</td></tr><tr><td>SET</td><td>NSD1/2 (KMT3B)</td><td>1343</td><td>H3K36 – H4K44</td><td>Smp_137060</td></tr><tr><td>SET</td><td>SET8 (KMT5A)</td><td>409</td><td>H4K20</td><td>Smp_055310</td></tr><tr><td>SET</td><td>SUV 39H2 (KMT1B)</td><td>586</td><td>H3K9</td><td>Smp_027300</td></tr><tr><td>SET</td><td>SUV4-20H1 (KMT5C)</td><td>613</td><td>H4K20</td><td>Smp_062530</td></tr><tr><td>SET</td><td>SETD2</td><td>1575</td><td>H3K36</td><td>Smp_133910</td></tr><tr><td>SET</td><td>SETD1B</td><td>1720/1822</td><td>H3K4</td><td>Smp_140390</td></tr><tr><td>SET</td><td>SETDB</td><td>918/1032</td><td></td><td>Smp_150850</td></tr><tr><td>SET</td><td>SETMAR</td><td>250</td><td>H3K9</td><td>Smp_043580</td></tr><tr><td>SET</td><td>SET/MYND4</td><td>782</td><td></td><td>Smp_000700</td></tr><tr><td>SET</td><td>SET/MYND4</td><td>527</td><td></td><td>Smp_124950</td></tr><tr><td>SET</td><td>SET/MYND5</td><td>423/429/433</td><td></td><td>Smp_121610</td></tr><tr><td>DOT1</td><td>DOT1L (KMT4)</td><td></td><td>H3K79</td><td>Smp_165000</td></tr><tr><td>PRMT</td><td>PRMT1</td><td>252/359/334</td><td>H4R3</td><td>Smp_029240</td></tr><tr><td>PRMT</td><td>PRMT3</td><td>1564</td><td></td><td>Smp_127950</td></tr><tr><td>PRMT</td><td>PRMT4/CARM1</td><td>737</td><td>H3R2 – H3R17 – H3R26</td><td>Smp_070340</td></tr><tr><td>PRMT</td><td>PRMT5</td><td>630</td><td>H2A – H4</td><td>Smp_171150</td></tr><tr><td>PRMT</td><td>PRMT7</td><td>755</td><td></td><td>Smp_025550</td></tr><tr><td rowspan="5">HDM</td><td>KDM1</td><td>LSD1A</td><td>1043</td><td></td><td>Smp_150560</td></tr><tr><td>KDM1</td><td>LSD1A</td><td>916</td><td></td><td>Smp_160810</td></tr><tr><td>KDM1</td><td>LSD1 (KDM1)</td><td>1073</td><td>H3K4 – H3K9</td><td>Smp_162940</td></tr><tr><td>JmjC</td><td>JMJD1B (KDM3)</td><td>273</td><td>H3K9</td><td>Smp_161410</td></tr><tr><td>JmjC</td><td>JMJD2C (KDM4C)</td><td>1136</td><td>H3K9 – H3K36</td><td>Smp_132170</td></tr></tbody></table>

(Continued)
<table><caption>**Table 1 | Continued**</caption><thead><tr><th>HME type</th><th>Class</th><th>Closest human ortholog</th><th>Size (aa)</th><th>Substrate specificity</th><th>Gene Id<sup>a</sup></th></tr></thead><tbody><tr><td></td><td>JmjC</td><td>JMJD4</td><td>809</td><td></td><td>Smp_147870</td></tr><tr><td></td><td>JmjC</td><td>JMJD6</td><td>839</td><td></td><td>Smp_137240</td></tr><tr><td></td><td>JmjC</td><td>JHDM1D (KDM7)</td><td>653</td><td>H3K36</td><td>Smp_127230</td></tr><tr><td></td><td>JmjC</td><td>Jarid (KDM5)</td><td>2372</td><td>H3K4</td><td>Smp_156290</td></tr><tr><td></td><td>JmjC</td><td>jarid (KDM5)</td><td>1639</td><td>H3K4</td><td>Smp_019170</td></tr><tr><td></td><td>JmjC</td><td>UTX (KDM6A)</td><td>1137</td><td>H3K27</td><td>Smp_034000</td></tr></tbody></table>

- If selective inhibitors can be developed, the targeted HME has to be essential to the parasite so that its inhibition is lethal. Many HMEs have overlapping specificities (Table 1) and the inhibition of one may be compensated by another. However, in the case of enzymes methylating (EZH) or demethylating (Utx) H3K27, only one isoform is present in schistosomes (Table 1), against two in humans, suggesting that these enzymes may represent particularly sensitive targets.

In consequence, for these particular therapeutic targets, an approach involving both target validation, notably by transcript knockdown (RNAi), and structural studies to determine specificities in the structure of the catalytic pocket, is essential.

The S. mansoni genome encodes 55 HMEs involved in protein acetylation/deacetylation or methylation/demethylation (Table 1) (Pierce et al., 2012). Some of these, including the class I HDACs (Oger et al., 2008) and the sirtuins (Lancelot et al., 2013) have been cloned and characterized and preliminary choices of targets can be made based on their degree of sequence conservation. In addition, the human orthologs of several of these enzymes are known to function only as part of a multiprotein complex as is the case for the H3K27 methyltransferase EZH2 mentioned previously (McCabe et al., 2012). High-throughput screening of the human enzyme has been done, using a five protein complex, but the resources devoted to carry out such a strategy for the development of anti-cancer therapies are not available in the case of neglected tropical diseases. Therefore, although the schistosome ortholog of EZH2 is a unique target, its screening will pose difficulties and it is therefore preferable to target an enzyme that is active on its own (like the H3K79 methyltransferase DOT1L for example Daigle et al., 2013). In addition, certain HMEs are very large proteins (Table 1) and although the production of truncated proteins containing the catalytic domain can be envisaged, this may affect both enzyme activity and the conformation of the catalytic pocket, limiting the relevance of screening or structural data.

These considerations reduce the choice of viable HME targets in schistosomes, but a large number remain. Two filters can be used to further limit the choice: the use of HME class inhibitors to determine whether enzyme families contain potential targets and transcript knockdown using RNAi to validate individual HMEs as stand-alone therapeutic targets. RNAi is still the only available means to achieve targeted knockdown of gene function in schistosomes, but its efficacy is transcript-dependent and phenotypes are not always observed (Stefanic et al., 2010).

## DEVELOPMENT OF SELECTIVE INHIBITORS AS DRUGS: THE CHALLENGES

An illustration of the strategy that can be employed to designate therapeutic targets is provided by the S. mansoni class I histone deacetylase 8 (SmHDAC8). Of the class I HDACs, SmHDAC8 was initially designated as a potential target for two reasons. First, transcript expression levels of SmHDAC8 are higher than those of SmHDAC1 and SmHDAC3 throughout the life cycle, notably in adult female worms (Oger et al., 2008). It is notable that HDAC 8 transcript levels are generally much lower than those of HDAC1 and 3 in normal human cells, but are markedly upregulated in some cancer cell lines and tissues (Nakagawa et al., 2007). Second, the analysis of the primary sequence of SmHDAC8 showed that it is less well-conserved compared to its human ortholog than the other two class I enzymes. This is demonstrated by the sequence alignment and phylogenetic analysis shown in Figure 2. The alignment shows that the essential residues for HDAC activity are conserved, but that the catalytic domain sequence contains insertions and substitutions that might indicate a change in architecture of the catalytic pocket, notably the replacement of a methionine (M274) in the human HDAC8 by a histidine (H292) in SmHDAC8. The status of SmHDAC8 as a stand-alone therapeutic target was enhanced by the use of generic HDAC8 inhibitors (unpublished results) and particularly by transcript knockdown using RNA interference (Marek et al., 2013). The latter showed that treatment of schistosomula with double-stranded RNA, followed by their injection i.v. into mice and harvesting of surviving worms 35 days later, led to a reduced worm recovery compared to mice treated with dsRNA encoding green fluorescent protein as a control.

Whilst molecular modeling seemed to show that only the charge difference within the catalytic pocket provided by His292 differentiated SmHDAC8 from its human counterpart, structural analysis by X-ray crystallography demonstrated a further important difference (Marek et al., 2013). In the schistosome enzyme, amino acid substitutions surrounding the catalytic pocket allow a change in the configuration of the side chain of phenylalanine 151. The side chain of the equivalent residue in human HDAC8, Phe152, adopts an obligatory flipped-in conformation that contributes to the narrow, hydrophobic tunnel accommodating the substrate or inhibitors. In contrast, in SmHDAC8 this side chain is free to adopt a flipped-out configuration, allowing the pocket to accommodate more bulky substrates or inhibitors. This difference, together with the charge difference, allows the possibility of identifying selective inhibitors for the schistosome

<footnote>*Validated by cDNA cloning.</footnote>

<footnote><sup>a</sup>Gene ID according to the genome annotation.</footnote>
<figure><img src="image_4.png" /><figcaption>**FIGURE 2 |** Phylogenetic relationship amongst human and Schistosome class I HDACs. Panel (A) displays a phylogenetic tree [ML, consensus from maximum likelihood; NJ, neighbor joining; and MB, Mr Bayes methods) built using the amino acid sequences from the catalytic domains of class I HDACs (HDAC1, 3, and 8) present in *S. mansoni* (Sm), *S. japonicum* (Sj), and *S. haematobium* (Sh) and their human (h) orthologs. Panel (B) shows an amino acid sequence alignment constructed with MAFFT with the class I HDACs used above. Conserved (black background and white letters), conservative changes (gray</figcaption></figure>

background and black letters) and less conservative changes (gray background and white letters) in amino acid positions are highlighted. Black (F151) and brown (Y99) triangles show residues that despite being conserved present conformational differences between human and Schistosomes HDAC8 (Marek et al., 2013). The blue triangle shows the position of the replacement of a methionine (M274) in the human HDAC8 by a histidine (H292) in SmHDAC8 (boxed amino acids below blue triangle). Residues that interact with the SmHDAC8 specific inhibitor J1075 are highlighted with blue lines.

enzyme. Indeed, an *in silico* screen, based on the crystal structure of SmHDAC8 and involving the docking of a large number of potential inhibitors, led to the identification of an inhibitor, J1075 (Marek et al., 2013), which had greatly improved selectivity for SmHDAC8 compared to human class I and II HDACs. Moreover, this inhibitor caused dose and time-dependent death of schistosome larvae in culture via the induction of apoptosis. Optimization of this and other inhibitors identified by this strategy is ongoing, and further potential drug precursors have

been identified (Stolfa et al., 2014). It is also notable that the structural specificities of the SmHDAC8 enzyme compared to the human ortholog are shared with HDAC8 in other flatworms, including other schistosome species (Figure 2), *Echinococcus* sp. and *Clonorchis sinensis* (Marek et al., 2013). Therefore, drugs effective against *S. mansoni* may well be applicable to these other species.

This example proves the concept that individual epigenetic enzymes can be valid therapeutic targets, and that, even though
these enzymes generally have conserved catalytic domains, suffi-
cient differences in structure can exist to allow the development
of selective inhibitors that are drug precursors. It remains to be
seen, however, whether these inhibitors can be made sufficiently
selective to preclude potentially harmful side effects and whether
they can be developed into drugs useable in a single oral dose in
humans.

CONCLUSIONS

Epigenetic processes provide a wealth of potential therapeutic
targets for the development of novel therapies against schistoso-
miasis and other parasitic diseases. The most readily exploitable of
these targets are the HMEs, as well as perhaps Dnmt2, which lend
themselves to target-based drug discovery strategies, necessary to
ensure the development of parasite-selective drugs. A structure-
based strategy has been initiated for S. mansoni HDAC8, involving
the solution of the 3D structure of the catalytic domain and
*in silico* docking of potential inhibitors. However, such enzymes
can also be screened directly using random compound libraries
and high-throughput methodologies using enzyme inhibition as
the read-out. Moreover, the existing extensive libraries of HME
inhibitors can be used for phenotypic screening of compound
libraries for lethal effects on the parasite itself. In all these cases,
since the molecular mechanism of action of the drug precur-
sor is known, time will be saved in the process of optimizing
selectivity.

The remaining potential targets discussed above, notably the
histone modification readers, like bromodomain proteins, and
miRNAs, require validation as effective targets and pose greater
challenges for drug development. In order to identify the schis-
tosome bromodomain proteins, an exhaustive analysis of the
genomic data has still to be done. However, since both bro-
modomain proteins and miRNAs are under intensive investiga-
tion, particularly for the development of anti-cancer therapies,
methodologies will be developed that could be exploited and
adapted for the treatment of parasitic diseases.

Work on all the potential targets discussed here can bene-
fit from the increasing knowledge base and compound libraries
accrued, notably in the development of anti-cancer therapies.
This "piggy-backing" approach (Dissous and Grevelding, 2010)
holds great promise and can in part mitigate the relative lack of
investment in efforts to improve the control and treatment of
schistosomiasis and the other neglected parasitic diseases.

ACKNOWLEDGMENTS

This work and the authors of this manuscript received funding
from the European Union's Seventh Framework Programme for
research, technological development and demonstration under
grant agreements nos. 241865 (SEtTReND) and 602080 (A-
ParaDDisE). The work of Julien Lancelot, Stéphanie Caby,
Raymond J. Pierce, and Alejandro Cabezas-Cruz was also sup-
ported by institutional funds from the Centre National de la
Recherche Scientifique (CNRS), the Institut National de la Santé
et de la Recherche Médicale (INSERM) and the Université de
Lille 2. GO was also supported by contracts CNPq 573839/2008-
5, 403049/20121, and 309312/2012-4), CAPES (070/13), and
FAPEMIG (PPM-00439-10).

REFERENCES

Anderson, L., Pierce, R. J., and Verjovski-Almeida, S. (2012). *Schistosoma mansoni* histones: from transcription to chromatin regulation; an *in silico* analysis. *Mol. Biochem. Parasitol.* 183, 105–114. doi: 10.1016/j.molbiopara.2012.03.001

Arrowsmith, C. H., Bountra, C., Fish, P. V., Lee, K., and Schapira, M. (2012). Epigenetic protein families: a new frontier for drug discovery. *Nat. Rev. Drug Discov.* 11, 384–400. doi: 10.1038/nrd3674

Azzi, A., Cosseau, C., and Grunau, C. (2009). *Schistosoma mansoni*: developmental arrest of miracidia treated with histone deacetylase inhibitors. *Exp. Parasitol.* 121, 88–291. doi: 10.1016/j.exppara.2008.11.010

Badeaux, A. I., and Shi, Y. (2013). Emerging roles for chromatin as a signal integration and storage platform. *Nat. Rev. Mol. Cell Biol.* 14, 211–224. doi: 10.1038/nrm3545

Bantscheff, M., Hopf, C., Savitski, M. M., Dittmann, A., Grandi, P., Michon, A. M., et al. (2011). Chemoproteomics profiling of HDAC inhibitors reveals selective targeting of HDAC complexes. *Nat. Biotechnol.* 29, 255–265. doi: 10.1038/nbt.1759

Baubec, T., and Schübeler, D. (2014). Genomic patterns and context-
specific interpretation of DNA methylation. *Genet. Dev.* 25, 85–92. doi:
10.1016/j.gde.2013.11.015

Berriman, M., Haas, B., Loverde, P. T., Wilson, R. A., Dillon, G. P., Cerqueira, G. C., et al. (2009). The Genome of the blood fluke *Schistosoma mansoni*. *Nature* 462, 352–358. doi: 10.1038/nature08160

Bertin, B., Oger, F., Cornette, J., Caby, S., Noël, C., Capron, M., et al. (2006). *Schistosoma mansoni* CBP/P300 has a conserved domain structure and interacts functionally with the nuclear receptor SmFtz-F1. *Mol. Biochem. Parasitol.* 146, 180–191. doi: 10.1016/j.molbiopara.2005.12.006

Bowers, E. M., Yan, G., Mukherjee, C., Orry, A., Wang, L., Holbert, M. A., et al. (2010). Virtual ligand screening of the p300/CBP histone acetyltransferase: identification of a selective small molecule inhibitor. *Chem. Biol.* 17, 471–482. doi: 10.1016/j.chembiol.2010.03.006

Boyle, J. P., Shoemaker, C. B., and Yoshino, T. P. (2003). Using RNA interference to manipulate endogenous gene expression in *Schistosoma mansoni* sporocysts. *Mol. Biochem. Parasitol.* 128, 205–215. doi: 10.1016/S0166-6851(03)00078-1

Cai, P., Hou, N., Piao, X., Liu, S., Liu, H., Yang, F., et al. (2011). Profiles of small non-coding RNAs in *Schistosoma japonicum* during development. *PLoS Negl. Trop. Dis.* 5:e1256. doi: 10.1371/journal.pntd.0001256

Carneiro, V. C., de Abreu da Silva, I. C., Torres, E. J., Caby, S., Lancelot, J., Vanderstraete, M., et al. (2014). Epigenetic changes modulate schistosome egg formation and are a novel target for reducing transmission of schistosomiasis. *PLoS Pathog.* 8:e1004116. doi: 10.1371/journal.ppat.1004116

Chung, C.-W., Coste, H., White, J. H., Miguez, O., Wilde, J., Gosmini, R. L., et al. (2011). Discovery and characterization of small molecule inhibitors of the BET family bromodomains. *J. Med. Chem.* 54, 3827–3838. doi: 10.1021/jm200108t

Coeli, R., Baba, E. H., Araujo, N., Coelho, P. M. C. Z., and Oliveira, G. (2013). Praziquantel Treatment Decreases *Schistosoma mansoni* genetic diversity in experimental infections. *PLoS Negl. Trop. Dis.* 7:e2596. doi: 10.1371/journal.pntd.0002596

Colley, D. G., Bustinduy, A. L., Secor, W. E., and King, C. H. (2014). Human
schistosomiasis. *Lancet* 383, 2253–2264. doi: 10.1016/SD140-6736(13)61949-2

Costinean, S., Zanesi, N., Pekarsky, Y., Tili, E., Volinia, S., Heerema, N., et al. (2006). Pre-B cell proliferation and lymphoblastic leukemia /high grade lymphoma in Eμ-miR155 transgenic mice. *Proc. Natl. Acad. Sci. U.S.A.* 103, 7024–7029. doi: 10.1073/pnas.0602266103

Daigle, S. R., Olhava, E. J., Therkelsen, C. A., Basavapathruni, A., Jin, L., Boriack-Sjodin, P. A., et al. (2013). Potent inhibition of DOT1L as treatment for MLL-fusion leukemia. *Blood* 122, 1017–1025. doi: 10.1182/blood-2013-04-497644

Day, T. A., Bennet, J. L., and Pax, R. A. (1992). Praziquantel: the enigmatic antiparasitic. *Parasitol. Today* 8, 342–344.

Deardorff, M. A., Bando, M., Nakato, R., Watrin, E., Itoh, T., Minamino, M., et al. (2012). HDAC8 mutations in Cornelia de Lange syndrome affect the cohesin acetylation cycle. *Nature* 489, 313–317. doi: 10.1038/nature11316

de Moraes Maciel, R., da Costa, R. F., de Oliveira, F. M., Rumjanek, F. D., and Fantappié, M. R. (2008). Protein acetylation sites mediated by *Schistosoma mansoni* GCN5. *Biochem. Biophys. Res. Commun.* 370, 53–56. doi: 10.1016/j.bbrc.2008.03.022

Dissous, C., and Grevelding, C. G. (2010). Piggy-backing the concept of cancer drugs for schistosomiasis treatment: a tangible perspective? *Trends Parasitol.* 27, 59–66. doi: 10.1016/j.pt.2010.09.001
Doenhoff, M. J., Cioli, D., and Utzinger, J. (2008). Praziquantel: mechanisms of action, resistance and new derivatives for schistosomiasis. *Curr. Opin. Infect. Dis.* 21, 659–667. doi: 10.1097/QCO.0b013e328318978f

Dowden, J., Hong, W., Parry, R. V., Pike, R. A., and Ward, S. G. (2010). Toward the development of potent and selective bisubstrate inhibitors of protein arginine methyltransferases. *Bioorg. Med. Chem. Lett.* 20, 2103–2105. doi: 10.1016/j.bmcl.2010.02.069

Du, J., Zhou, Y., Su, X., Yu, J. J., Khan, S., Jiang, H., et al. (2011). Sirt5 is a NAD-dependent protein lysine demalonylase and desuccinylase. *Science* 334, 806–809. doi: 10.1126/science.1207861

Dubois, F., Caby, S., Oger, F., Cosseau, C., Capron, M., Grunau, C., et al. (2009). Histone deacetylase inhibitors induce apoptosis, histone hyperacetylation and up-regulation of gene expression in *Schistosoma mansoni*. *Mol. Biochem. Parasitol.* 168, 7–15. doi: 10.1016/j.molbiopara.2009.06.001

Eder, J., Sedrani, R., and Wiesmann, C. (2014). The discovery of first-in-class drugs: origins and evolution. *Nat. Rev. Drug Discov.* 13, 577–587. doi: 10.1038/nrd4336

Elring, A. M., Harb, J. G., Neviani, P., Garton, C., Oaks, J. J., Spizzo, R., et al. (2010). miR-328 functions as an RNA decoy to modulate hnRNP E2 regulation of mRNA translation in leukemic blasts. *Cell* 140, 652–665. doi: 10.1016/j.cell.2010.01.007

Fagan, R. L., Crydermand, D. E., Kopelovich, L., Wallraty, L. L., and Brenner, C. (2013). Laccaic acid A is a direct, DNA-competitive inhibitor of DNA methyl-transferase 1. *J. Biol. Chem.* 288, 23858–23867. doi: 10.1074/jbc.M113.480517

Fallon, P. G., and Doenhoff, M. J. (1994). Drug-resistant schistosomiasis: resistance to praziquantel and oxamniquine induced in *Schistosoma mansoni* in mice is drug specific. *Am. J. Trop. Med. Hyg.* 51, 83–88.

Fantappié, M. R., Bastos de Oliveira, F. M., de Moraes Maciel, R., Rumjanek, F. D., Wu, W., and LoVerde, P. T. (2008). Cloning of SmNCoA-62, a novel nuclear receptor co-activator from *Schistosoma mansoni*: assembly of a complex with a SmRXR1/SmNR1 heterodimer, SmGCN5 and SmCBP1. *Int. J. Parasitol.* 38, 1133–1147. doi: 10.1016/j.ijpara.2008.02.003

Fantappié, M. R., Rodrigues, E., Gimba, P., and Rumjanek, F. D. (2001). Lack of DNA methylation in *Schistosoma mansoni*. *Exp. Parasitol.* 98, 162–166. doi: 10.1006/expr.2001.4630

Feldman, J. L., Dittenhafer-Reed, K. E., and Denu, J. M. (2012). Sirtuin catalysis and regulation. *J. Biol. Chem.* 287, 42419–42427. doi: 10.1074/jbc.R112.378877

Fenwick, A., Webster, J. P., Bosque-Oliva, E., Blair, L., Fleming, F. M., Zhang, Y., et al. (2009). The Schistosomiasis Control Initiative (SCI): rationale, development and implementation from 2002-2008. *Parasitology* 136, 1719–1730. doi: 10.1017/S0031182009990400

Filippakopoulos, P., and Knapp, S. (2014). Targeting bromodomains: epigenetic readers of lysine acetylation. *Nat. Rev. Drug. Discov.* 13, 337–356. doi: 10.1038/nrd4286

Filippakopoulos, P., Qi, J., Picaud, S., Shen, Y., Smith, W. B., Fedorov, O., et al. (2011). Selective inhibition of BET bromodomains. *Nature* 468, 1067–1073. doi: 10.1038/nature09504

Fromm, B., Worren, M. M., Hahn, C., Hovig, E., and Bachmann, L. (2013). Substantial loss of conserved and gain of novel microRNA families in flatworms. *Mol. Biol. Evol.* 30, 2619–2628. doi: 10.1093/molbev/mst155

Furdas, S. D., Shekfeh, S., Bissinger, E.-M., Wagner, J. M., Schlimme, S., Valkov, V., et al. (2011). Synthesis and biological testing of novel pyridoisothiazolones as histone acetyltransferase inhibitors. *Bioorg. Med. Chem.* 19, 3678–3689. doi: 10.1016/j.bmc.2011.01.063

Geyer, K. K., Chalmers, I. W., MacKintosh, N., Hirst, J. E., Geoghegan, R., Badets, M., et al. (2013). Cytosine methylation is a conserved epigenetic feature found throughout the phylum Platyhelminthes. *BMC Genomics* 14:462. doi: 10.1186/1471-2164-14-462

Geyer, K. K., Rodriguez Lopez, C. M., Chalmers, I. W., Munshi, S. E., Truscott, M., Heald, J., et al. (2011). Cytosine methylation regulates oviposition in the pathogenic blood fluke *Schistosoma mansoni*. *Nat. Commun.* 2, 424. doi: 10.1038/ncomms1433

Ghizzoni, M., Haisma, H. J., and Dekker, F. J. (2009). Reactivity of isothiazolones and isothiazolone-1-oxides in the inhibition of PCAF histone acetyl-transferase. *Eur. J. Med. Chem.* 44, 4855–4861. doi: 10.1016/j.ejmech.2009.07.025

Goll, M. G., Kirpekar, F., Maggert, K. A., Yoder, J. A., Hsieh, C. L., Zhang, X., et al. (2006). Methylation of tRNA Asp by the DNA methyltransferase homolog Dnmt2. *Science* 311, 395–398. doi: 10.1126/science.1120976

Gregoretti, I. V., Lee, Y. M., and Goodson, H. V. (2004). Molecular evolution of the histone deacetylase family: functional implications of phylogenetic analysis. *J. Mol. Biol.* 338, 17–31. doi: 10.1016/j.jmb.2004.02.006

Greiss, S., and Gartner, A. (2009). Sirtuin/Sir2 phylogeny, evolutionary considerations and structural conservation. *Mol. Cells* 28, 407–415. doi: 10.1007/s10059-009-0169-x

Haigis, M. C., Mostoslavsky, R., Haigis, K. M., Fahie, K., Christodoulou, D. C., Murphy, A. J., et al. (2006). SIRT4 inhibits glutamate dehydrogenase and opposes the effects of calorie restriction in pancreatic beta cells. *Cell* 126, 941–954. doi: 10.1016/j.cell.2006.06.057

Hao, L., Cai, P., Jiang, N., Wang, H., and Chen, Q. (2010). Identification and characterization of microRNAs and endogenous siRNAs in *Schistosoma japonicum*. *BMC Genomics* 11:55. doi: 10.1186/1471-2164-11-55

Helin, K., and Dhanak, D. (2013). Chromatin proteins and modifications as drug targets. *Nature* 502, 480–488. doi: 10.1038/nature12751

Hermann, A., Schmitt, S., and Jeltsch, A. (2003). The human Dnmt2 has residual DNA-(cytosine-C5) methyltransferase activity. *J. Biol. Chem.* 278, 31717–31721. doi: 10.1074/jbc.M30544820

Hewings, D. S., Rooney, T. P., Jennings, L. E., Hay, D. A., Schofield, C. J., Brennan, P. E., et al. (2012). Progress in the development and application of small molecule inhibitors of bromodomain-acetyl-lysine interactions. *J. Med. Chem.* 55, 9393–9413. doi: 10.1021/jm300915b

Hojfeldt, J. W., Agger, K., and Helin, K. (2013). Histone lysine demethylases as targets for anticancer therapy. *Nat. Rev. Drug. Discov.* 12, 917–930. doi: 10.1038/nrd4154

Holt, C. A., and Roth, A. E. (2004). The Nash equilibrium: a perspective. *Proc. Natl. Acad. Sci. U.S.A.* 101, 3999–4002. doi: 10.1073/pnas.0308738101

Huang, J., Hao, P., Chen, H., Hu, W., Yan, Q., Liu, F., et al. (2009). Genome-wide identification of *Schistosoma japonicum* microRNAs using a deep-sequencing approach. *PLoS ONE* 4:e8206. doi: 10.1371/journal.pone.0008206

Hubbard, B. P., and Sinclair, D. A. (2014). Small molecule Sirt1 activators for the treatment of aging and age-related diseases. *Trends Pharmacol. Sci.* 35, 146–154. doi: 10.1016/j.tips.2013.12.004

Hubbert, C., Guardiola, A., Shao, R., Kawaguchi, Y., Ito, A., Nixon, A., et al. (2002). HDAC6 is a microtubule-associated deacetylase. *Nature* 417, 455–458. doi: 10.1038/417455a

Janssen, H. L., Reesink, H. W., Lawitz, E. J., Zeuzem, S., Rodriguez-Torres, M., Patel, K., et al. (2013). Treatment of HCV infection by targeting microRNA. *N. Engl. J. Med.* 368, 1685–1694. doi: 10.1056/NEJMoa1209026

Jiang, H., Khan, S., Wang, Y., Charron, G., He, B., Sebastian, C., et al. (2013). Sirt6 regulates TNF- secretion through hydrolysis of long-chain fatty acyl lysine. *Nature* 496, 110–113. doi: 10.1038/nature12038

Kohn, A. B., Anderson, P. A., Roberts-Misterly, J. M., and Greenberg, R. M. (2001). Schistosome calcium channel beta subunits. Unusual modulatory effects and potential role in the action of the antischistosomal drug praziquantel. *J. Biol. Chem.* 276, 36873–36876. doi: 10.1074/jbc.C100273200

Kovacs, J. J., Murphy, P. J., Gaillard, S., Zhao, X., Wu, J. T., Nicchitta, C. V., et al. (2005). HDAC6 regulates Hsp90 acetylation and chaperone-dependent activation of glucocorticoid receptor. *Mol. Cell* 18, 601–607. doi: 10.1016/j.molcel.2005.04.021

Kraus, V., and Reuter, G. (2011). DNA methylation in *Drosophila* – a critical evaluation. *Prog. Mol. Biol. Transf. Sci.* 101, 177–191. doi: 10.1016/B978-0-12-387685-0-00003-2

Kruidenier, L., Chung, C.-W., Cheng, Z., Liddle, J., Che, K., Joberty, G., et al. (2012). A selective jumonji H3K27 demethylase inhibitor modulates the proinflammatory macrophage response. *Nature* 488, 404–408. doi: 10.1038/nature11262

Lieberman, J., and Sarnow, P. (2013). Micromanaging hepatitis C virus. *N. Engl. J. Med.* 368, 1741–1743. doi: 10.1056/NEJMMe1301348

Lancelot, J., Caby, S., Dubois-Abdesselem, F., Vanderstraete, M., Trolet, J., Oliveira, G., et al. (2013). *Schistosoma mansoni* sirtuins: characterization and potential as chemotherapeutic targets. *PLoS Negl. Trop. Dis.* 7:e2428. doi: 10.1371/journal.pntd.0002428

Li, Z., and Rana, T. M. (2014). Therapeutic targeting of microRNAs: current status and future challenges. *Nat. Rev. Drug Discov.* 13, 622–638. doi: 10.1038/nrd4359

Ling, H., Fabbri, M., and Calin, G. A. (2013). MicroRNAs and other non-coding RNAs as targets for anticancer drug development. *Nat. Rev. Drug. Discov.* 12, 847–865. doi: 10.1038/nrd4140
Maes, T., Tirapu, I., Mascaro, C., Ortega, A., Estiarte, A., and Valls, N. (2013). Preclinical characterization of a potent and selective inhibitor of the histone demethylase KDM1A for MLL leukemia. *J. Clin. Oncol.* 31(suppl.):e13543.

Marek, M., Kannan, S., Hauser, A. T., Moraes Mourão, M., Caby, S., Cura, V., et al. (2013). Structural basis for the inhibition of histone deacetylase 8 (HDAC8), a key epigenetic player in the blood fluke *Schistosoma mansoni*. *PLoS Pathog.* 9:e1003645. doi: 10.1371/journal.ppat.1003645

McCabe, M. T., Ott, H. M., Ganji, G., Korenchuk, S., Thompson, C., Van Aller, G. S., et al. (2012). EZH2 inhibition as a therapeutic strategy for lymphoma with EZH2-activating mutations. *Nature* 492, 108–112. doi: 10.1038/nature11606

Melman, S. D., Steinauer, M. L., Cunningham, C., Kubatko, L. S., Mwangi, I. N., Wynn, N. B., et al. (2009). Reduced susceptibility to praziquantel among naturally occurring Kenya isolates of *Schistosoma mansoni*. *PLoS Negl. Trop. Dis.* 3:e504. doi: 10.1371/journal.pntd.0000504

Messerli, S. M., Kasinathan, R. S., Morgan, W., Spranger, S., and Greenberg, R. M. (2009). *Schistosoma mansoni* P-glycoprotein levels increase in response to praziquantel exposure and correlate with reduced praziquantel susceptibility. *Mol. Biochem. Parasitol.* 167, 54–59. doi: 10.1016/j.molbiopara.2009.04.007

Metzger, E., Wissmann, M., Yin, N., Müller, J. M., Schneider, R., Peters, A. H., et al. (2005). LSD1 demethylates repressive histone marks to promote androgen-receptor dependent transcription. *Nature* 437, 436–439. doi: 10.1038/nature04020

Nakagawa, M., Oda, Y., Eguchi, T., Aishima, S.-I., Yao, T., Hosoi, F., et al. (2007). Expression profile of class I histone deacetylases in human cancer tissues. *Oncol. Rep.* 18, 769–774. doi: 10.3892/or.18.4.769

Nash, J. F. (1950). Equilibrium points in N-person games. *Proc. Natl. Acad. Sci. U.S.A.* 36, 48–49.

Nicodeme, E., Jeffrey, K. L., Schaefer, U., Beinke, S., Dewell, S., Chung, C.-W., et al. (2010). Suppression of inflammation by a synthetic histone mimic. *Nature* 468, 1119–1123. doi: 10.1038/nature09589

Norton, A. J., Gower, C. M., Lamberton, P. H., Webster, B. L., Lwambo, N. J., Blair, L., et al. (2010). Genetic consequences of mass human chemotherapy for *Schistosoma mansoni*: population structure pre- and post-praziquantel treatment in Tanzania. *Am. J. Trop. Med. Hyg.* 83, 951–957. doi: 10.4269/ajtmh.2010.10-0283

Ntziachristos, P., Tsigirigos, A., Welstead, C. G., Trimarchi, T., Bakogianni, S., Xu, L., et al. (2014). Contrasting roles of histone 3 lysine 27 demethylases in acute lymphoblastic leukaemia. *Nature*. doi: 10.1074/jbc.C100273200. [Epub ahead of print].

Oger, F., Dubois, F., Caby, S., Noël, C., Cornette, J., Bertin, B., et al. (2008). The class I histone deacetylases of the platyhelminth parasite *Schistosoma mansoni*. *Biochem. Biophys. Res. Commun.* 377, 1079–1084. doi: 10.1016/j.bbrcc.2008

Oliveira, K. C., Carvalho, M. L. P., Maracaja-Coutinho, V., Kitajima, J. P., and Verjovski-Almeida, S. (2011). Non-coding RNAs in schistosomes: an unexplored world. *An. Acad. Bras. Cienc.* 83, 673–694. doi: 10.1590/S0001-37652011000200026

Perrin, C., Lepesant, J. M., Roger, E., Duval, D., Fneich, S., Thuillier, V., et al. (2013). *Schistosoma mansoni* mucin gene (SmPoMuc) expression: epigenetic control to shape adaptation to a new host. *PLoS Pathog.* 9:e1003571. doi: 10.1371/journal.ppat.1003571

Piaz, F. D., Vassallo, A., Rubio, O. C., Castellano, S., Sbardella, G., and De Tommasi, N. (2011). Chemical biology of histone acetyltransferase natural compounds modulators. *Mol. Divers.* 15, 401–416. doi: 10.1007/s11030-010-9299-5

Pierce, R. J., Dubois-Abdesselem, F., Lancelot, J., Andrade, L., and Oliveira, G. (2012). Targeting schistosome histone modifying enzymes for drug development. *Curr. Pharm. Des.* 18, 3567–3578. doi: 10.2174/138161212801327248

Raddatz, G., Guzzardo, P. M., Olova, N., Fantappié, M. R., Rampp, M., Schaeffer, M., et al. (2013). Dnmt2-dependent methylomes lack defined DNA methylation patterns. *Proc. Natl. Acad. Sci. U.S.A.* 110, 8627–8631. doi: 10.1073/pnas.1306723110

Sayed, A. A., Simeonov, A., Thomas, C. J., Inglese, J., Austin, C. P., and Williams, D. L. (2008). Identification of oxadiazoles as new drug leads for the control of schistosomiasis. *Nat. Med.* 14, 407–412. doi: 10.1038/nm1737

Schaefer, M., and Lyko, F. (2010). Solving the Dnmt2 enigma. *Chromosoma* 119, 35–40. doi: 10.1007/s00412-009-0240-6

Schenk, T., Chen, W. C., Göllner, S., Howell, L., Jin, L., Hebestreit, K., et al. (2012). Inhibition of the LSD1 (KDM1A) demethylase reactivates the all-trans-retinoic differentiation pathway in acute myeloid leukemia. *Nature Med.* 18, 605–611. doi: 10.1038/nm.2661

Simoes, M. C., Lee, J., Djikeng, A., Cerqueira, G. C., Zerlotini, A., Silva-Pereira, R. A., et al. (2011). Identification of *Schistosoma mansoni* microRNAs. *BMC Genomics* 12:47. doi: 10.1186/1471-2164-12-47

Spannhoff, A., Hauser, A. T., Heinke, R., Sippl, W., and Jung, M. (2009). The emerging therapeutic potential of histone methyltransferase and demethylase inhibitors. *CheMedChem* 4, 1568–1582. doi: 10.1002/cmdc.200900301

Stefanic, S., Dvorak, J., Horn, M., Braschi, S., Sojka, D., Ruelas, D. Q., et al. (2010). RNA interference in *Schistosoma mansoni* schistosomula: selectivity, sensitivity and operation for larger-scale screening. *PLoS Negl. Trop. Dis.* 4:e850. doi: 10.1371/journal.pntd.0000850

Stolfa, D. A., Marek, M., Lancelot, J., Hauser, A. T., Walter, A., Leproult, E., et al. (2014). Molecular basis for the antiparasitic activity of a mercaptoacetamide derivative that inhibits histone deacetylase 8 (HDAC8) from the human pathogen *Schistosoma mansoni*. *J. Mol. Biol.* doi: 10.1016/j.jmb.2014.03.007. [Epub ahead of print].

Strahl, B., and Allis, C. (2000). The language of covalent histone modifications. *Nature* 403, 41–45. doi: 10.1038/47412

Streseman, C., and Lyko, F. (2008). Modes of action of the DNA methyltransferase inhibitors azacytidine and decitabine. *Int. J. Cancer* 123, 8–13. doi: 10.1002/ijc.23607

Sun, J., Wang, S., Li, C., Ren, Y., and Wang, J. (2014). Novel expression profiles of microRNAs suggest that specific miRNAs regulate gene expression for the sexual maturation of female *Schistosoma japonicum* after pairing. *Parasit. Vectors* 7:177. doi: 10.1186/1756-3305-7-177

Swinney, D. C., and Anthony, J. (2011). How were new medicines discovered? *Nat. Rev. Drug. Discov.* 10, 507–519. doi: 10.1038/nrd3480

Tan, M., Peng, C., Anderson, K. A., Chhoy, P., Xie, Z., Dai, L., et al. (2014). Lysine glutarylation is a protein posttranslational modification regulated by Sirt5. *Cell Metab.* 19, 605–617. doi: 10.1016/j.cmet.2014.03.014

Wang, Z., Xue, X., Sun, J., Luo, R., Xu, X., Jiang, Y., et al. (2010). An ≪ in-depth ≫ description of the small non-coding RNA population of *Schistosoma japonicum* schistosomulum. *PLoS Negl. Trop. Dis.* 4:e596. doi: 10.1371/journal.pntd.0000596

Warnecke, P. M., Stirzaker, C., Song, J., Grunau, C., Melki, J. R., and Clark, S. J. (2002). Identification and resolution of artifacts in bisulfite sequencing. *Methods* 27, 101–107. doi: 10.1016/S1046-2023(02)00060-9

West, A. C., and Johnstone, R. W. (2014). New and emerging HDAC inhibitors for cancer treatment. *J. Clin. Invest.* 124, 30–39. doi: 10.1172/JCI69738

Xue, X., Sun, J., Zhang, Q., Wang, Z., Huang, Y., and Pan, W. (2008). Identification and characterization of novel microRNAs from *Schistosoma japonicum*. *PLoS ONE* 3:e4034. doi: 10.1371/journal.pone.0004034

Yoo, C. B., and Jones, P. A. (2006). Epigenetic therapy of cancer. Past, present, and future. *Nat. Rev. Drug Discov.* 5, 37–50. doi: 10.1038/nrd1930

Young, N. D., Jex, A. R., Li, B., Liu, S., Yang, L., Xiong, Z., et al. (2012). Whole-genome sequence of *Schistosoma haematobium*. *Nat. Genet.* 44, 221–225. doi: 10.1038/ng.1065

Yuan, H., and Mamorstein, R. (2012). Structural basis for sirtuin activity and inhibition. *J. Biol. Chem.* 287, 42428–42435. doi: 10.1074/jbc.R112.372300

Zhou, Y., Zheng, H., Chen, Y., Zhang, L., Wang, K., Guo, J., et al. (2009). The *Schistosoma japonicum* genome reveals features of host-parasite interplay. *Nature* 460, 345–351. doi: 10.1038/nature08140

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Received: 26 June 2014; accepted: 24 August 2014; published online: 10 September 2014.

Citation: Cabezas-Cruz A, Lancelot J, Caby S, Oliveira G and Pierce RJ (2014) *Epigenetic control of gene function in schistosomes: a source of therapeutic targets?* *Front. Genet.* 5:317. doi: 10.3389/fgene.2014.00317

This article was submitted to Evolutionary and Genomic Microbiology, a section of the journal Frontiers in Genetics.

Copyright © 2014 Cabezas-Cruz, Lancelot, Caby, Oliveira and Pierce. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.