Next Article in Journal
Host–Pathogen Interactions between Xanthomonas fragariae and Its Host Fragaria × ananassa Investigated with a Dual RNA-Seq Analysis
Next Article in Special Issue
Central Asian Rodents as Model Animals for Leishmania major and Leishmania donovani Research
Previous Article in Journal
First Cases of Natural Infections with Borrelia hispanica in Two Dogs and a Cat from Europe
Previous Article in Special Issue
Leishmania Immunity: Advancing Immunotherapy and Vaccine Development
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Absence of C-5 DNA Methylation in Leishmania donovani Allows DNA Enrichment from Complex Samples

by
Bart Cuypers
1,2,†,
Franck Dumetz
1,†,‡,
Pieter Meysman
2,
Kris Laukens
2,
Géraldine De Muylder
1,
Jean-Claude Dujardin
1,3 and
Malgorzata Anna Domagalska
1,*
1
Molecular Parasitology, Institute of Tropical Medicine, 2000 Antwerp, Belgium
2
ADReM Data Lab, Department of Computer Science, University of Antwerp, 2000 Antwerp, Belgium
3
Department of Biomedical Sciences, University of Antwerp, 2000 Antwerp, Belgium
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Present address: Department of Pathology, University of Cambridge, Cambridge CB2 1QP, UK.
Microorganisms 2020, 8(8), 1252; https://doi.org/10.3390/microorganisms8081252
Submission received: 11 July 2020 / Revised: 12 August 2020 / Accepted: 12 August 2020 / Published: 18 August 2020
(This article belongs to the Special Issue Leishmania and Leishmaniasis)

Abstract

:
Cytosine C5 methylation is an important epigenetic control mechanism in a wide array of eukaryotic organisms and generally carried out by proteins of the C-5 DNA methyltransferase family (DNMTs). In several protozoans, the status of this mechanism remains elusive, such as in Leishmania, the causative agent of the disease leishmaniasis in humans and a wide array of vertebrate animals. In this work, we showed that the Leishmania donovani genome contains a C-5 DNA methyltransferase (DNMT) from the DNMT6 subfamily, whose function is still unclear, and verified its expression at the RNA level. We created viable overexpressor and knock-out lines of this enzyme and characterized their genome-wide methylation patterns using whole-genome bisulfite sequencing, together with promastigote and amastigote control lines. Interestingly, despite the DNMT6 presence, we found that methylation levels were equal to or lower than 0.0003% at CpG sites, 0.0005% at CHG sites, and 0.0126% at CHH sites at the genomic scale. As none of the methylated sites were retained after manual verification, we conclude that there is no evidence for DNA methylation in this species. We demonstrated that this difference in DNA methylation between the parasite (no detectable DNA methylation) and the vertebrate host (DNA methylation) allowed enrichment of parasite vs. host DNA using methyl-CpG-binding domain columns, readily available in commercial kits. As such, we depleted methylated DNA from mixes of Leishmania promastigote and amastigote DNA with human DNA, resulting in average Leishmania:human enrichments from 62× up to 263×. These results open a promising avenue for unmethylated DNA enrichment as a pre-enrichment step before sequencing Leishmania clinical samples.

1. Introduction

DNA methylation is an epigenetic mechanism responsible for a diverse set of functions across the three domains of life; eubacteria, archaebacteria, and eukaryota. In prokaryotes, many DNA methylation enzymes are part of so-called restriction-modification systems, which play a crucial role in their defense against phages and viruses. Prokaryotic methylation typically occurs on the C5 position of cytosine (cytosine C5 methylation), the exocyclic amino groups of adenine (adenine-N6 methylation), or cytosine (cytosine-N4 methylation) [1]. In eukaryotic species, DNA methylation is mostly restricted to 5-methylcytosine (me5C) and best characterized in mammals, where 70–80% of the CpG motifs are methylated [2]. As such, DNA methylation controls a wide range of important cellular functions, such as genomic imprinting, X-chromosome inactivation (in humans), gene expression, and the repression of transposable elements. Consequently, defects in genetic imprinting are associated with a variety of human diseases, and changes in DNA methylation patterns are a common hallmark of cancer [3,4]. Eukaryotic DNA methylation can also occur at CHG and CHH (where H is A, C, or T) sites [5], which was considered to occur primarily in plants. However, studies from the past decade demonstrate that CHG and CHH methylations are also frequent in several mammalian cell types, such as embryonic stem cells, oocytes, and brain cells [5,6,7,8].
Me5C methylation is mediated by a group of enzymes called C-5 DNA methyltransferases (DNMTs). This ancient group of enzymes shares a common ancestry, and their core domains are conserved across prokaryotes and eukaryotes [1]. Different DNMT subfamilies have developed distinct roles within epigenetic control mechanisms. For example, in mammals, DNMT3a and DNMT3b are responsible for de novo methylation, such as during germ cell differentiation and early development, or in specific tissues undergoing dynamic methylation [9]. In contrast, DNMT1 is responsible for maintaining methylation patterns, particularly during the S phase of the cell cycle, where it methylates the newly generated hemimethylated sites on the DNA daughter strands [10]. Some DNMTs have also changed substrate over the course of evolution. A large family of DNMTs, called DNMT2, has been shown to methylate the 38th position of different tRNAs to yield ribo-5-methylcytidine (rm5C) in a range of eukaryotic organisms, including humans [11], mice [12], Arabidopsis thaliana [13], and Drosophila melanogaster [14]. Therefore, DNMT2s are now often referred to as ‘tRNA methyltransferases’ or trDNMT and are known to carry out diverse regulatory functions [15]. However, in other eukaryotic taxa, DNMT2 appears to be a genuine DNMT, as DNMT2 can catalyze DNA methylation in Plasmodium falciparum [16] and Schistosoma mansoni [17]. In Entamoeba histolytica, both DNA and RNA can be used as substrates for DNMT2 [18,19]. The increase in available reference genomes of non-model eukaryotic species has recently resulted in the discovery of new DNMTs, such as DNMT5, DNMT6, or even SymbioLINE-DNMT—a massive family of DNMTs, so far only found in the dinoflagellate Symbiodinium [20].
Indeed, DNMT-mediated C5 methylation has been shown to be of major functional importance in a wide array of eukaryotic species, including protozoans, such as Toxoplasma gondii and Plasmodium [16,21]. In contrast, studies have failed to detect any C5 DNA methylation in eukaryotic species, such as Caenorhabditis elegans, Saccharomyces cerevisiae, and Schizosaccharomyces pombe [22,23]. In many other protozoans, the presence and potential role of DNA methylation remain elusive. This is especially true for Leishmania, a trypanosomatid parasite (Phylum Euglenozoa), despite its medical and veterinary importance. Leishmania is the causative agent of leishmaniasis—a disease that ranges from self-healing cutaneous lesions to lethal visceral leishmaniasis—in humans and a wide variety of vertebrate animals.
Leishmania features molecular biology that is remarkably different from other eukaryotes. This includes a system of polycistronic transcription of functionally unrelated genes [24]. The successful transcription of these cistrons depends on several known epigenetic modifications at the transcription start sites (acetylated histone H3) and transcription termination sites (β-D-glucosyl-hydroxymethyluracil, also called ‘Base J’). However, little research has been done towards other epigenetic modifications [25]. We were, therefore, interested in the 5-C methylation status of Leishmania, which has been poorly explored to date. In this context, a single study on a wide range of eukaryotic species lacking DNMT1 has reported the absence of CG-specific methylation in Leishmania major, however, using only a single sample of an unspecified life stage [26]. The study also does not comment on CHH and CHG specific methylation, which can be relevant as well. Contrastingly, Militello et al. demonstrated Me5C methylation in Trypanosoma brucei, another trypanosomatid species, although at low levels (0.01%) [27]. To clarify the status of C-5 DNA methylation in trypanosomatids and Leishmania, in particular, we presented the first comprehensive study of genomic methylation in Leishmania across different parasite life stages, making use of high-resolution whole-genome bisulfite sequencing.

2. Materials and Methods

2.1. In Silico Identification and Phylogeny of Putative DNMTs

To identify putative C-5 cytosine-specific DNA methylases in Leishmania donovani, we obtained the hidden Markov model (hmm) for this protein family from PFAM version 32.0 (Accession number: PF00145) [28]. The hmm search tool of hmmer-3.2.1 (hmmer.org) was then used with default settings to screen the LdBPKV2 reference genome [29] for this hmm signature. The initial pairwise alignment between the identified L. donovani and T. brucei C5 DNA MTase was carried out with T-COFFEE V_11.00.d625267.
To construct a comprehensive phylogenetic tree of the C5 DNA MTase family, including members found in trypanosomatid species, we modified the approach from Ponts et al. [16]. Firstly, we downloaded the putative proteomes of a wide range of prokaryotic and eukaryotic species. These species were selected to cover the different C5 DNA MTase subfamilies [1]. Specifically, the following proteomes were obtained: Leishmania amazonensis MHOMBR71973M2269 [30], Leishmania braziliensis MHOMBR75M2904 [31], Leishmania major Friedlin [24], Leishmania tarentolae Parrot-TarII [32], Leptomonas seymouri ATCC30220 [33], Trypanosoma brucei TREU92 [34], Trypanosoma congolense IL3000_2019 [35], Trypanosoma cruzi CL Brener [20], Trypanosoma grayi ANR4 [36] and Trypanosoma vivax Y486 from TriTrypDB v46 [37], Plasmodium falciparum 3D7 and Plasmodium vivax P01 from PlasmodDB v46 [38,39,40], Cryptosporidium parvum Iowa II and Cryptosporidium hominis TU502 from CryptoDB v46 [41,42,43], Toxoplasma gondii ARI from ToxoDB v 46 [44,45], Agrobacterium tumefaciens (GCF_000971565.1) [46], Arabidopsis thaliana (GCF_000001735.4), Aspergillus fumigatus (GCF_000002655.1), Aspergillus nidulans FGSC A4 (GCF_000149205.2) [47], Bacillus subtilis 168 (GCF_000009045.1), Clostridium botulinum ATCC 3502 (GCF_000063585.1) [48], Drosophila melanogaster (GCF_000001215.4), Emiliana huxleyi (GCF_000372725.1) [49], Entamoeba histolytica HM-1:IMSS (GCF_000208925.1) [50], Escherichia coli K12 (GCF_000005845.2), Eubacterium rectale (GCF_000020605.1), Euglena gracilis Z1 (PRJNA298469) [51], Micromonas pusilla (GCF_000151265.2) [52], Neurospora crassa OR74A (GCF_000182925.2) [53], Nostoc punctiforme PCC 73102 (GCF_000020025.1), Homo sapiens GRCh38.p12 (GCF_000001405.38), Phaeodactylum tricornutum (GCF_000150955.2) [54], Saccharomyces cerevisiae S288C (GCF_000146045.2), Salmonella enterica CT18 (GCF_000195995.1) [55], Schizosaccharomyces pombe ASM294 (GCF_000002945.1) [56], Streptococcus pneumoniae R6 (GCF_000007045.1) [57], Thalassiosira pseudonana (GCF_000149405.2) and Treponema succinifaciens (GCF_000195275.1) [58] from NCBI, Ascobolus immersus RN42 [59] from the JGI Genome Portal (genome.jgi.doe.gov), and Danio rerio (GRCz11) from Ensembl (ensembl.org).
All obtained proteomes were then searched with the hmm signature for C5 DNA MTases, exactly as described above for L. donovani. All hits with an E-value <0.01 (i.e., 1 false-positive hit is expected in every 100 searches with different query sequences) were maintained, and all domains matching the query hmm were extracted and merged per protein. This set of sequences was aligned in Mega-X with the MUSCLE multiple sequence alignment algorithm [60,61] and converted to the PHYLIP format with the ALTER tool [62]. Phage sequences and closely related isoforms were removed for the phylogenetic tree to maintain visibility.
A maximum-likelihood tree of this alignment was generated with RAxML version 8.2.10 using the automatic protein model assignment algorithm (option: -m PROTGAMMAAUTO). RAxML was run in three steps: Firstly, 20 trees were generated, and only the one with the highest likelihood score was kept. Secondly, 1000 bootstrap replicates were generated. In a final step, the bootstrap bipartitions were drawn on the best tree from the first round. The tree was visualized in Figtree v1.4.4 (https://github.com/rambaut/figtree/).

2.2. Culturing and DNA Extraction for Bisulfite Sequencing

Promastigotes (extracellular life stage) of Leishmania donovani MHOM/NP/03/BPK282/0 cl4 (further called BPK282) and its genetically modified daughter lines (see below) were cultured in HOMEM (Gibco, Thermo Fisher Scientific, Waltham, MA, USA) supplemented with 20% (v/v) heat-inactivated fetal bovine serum at 26 °C. Amastigotes (intracellular life stage) of the same strain were obtained from three months’ infected golden Syrian hamster (Charles Rivers, Wilmington, MA, USA), as described in Dumetz et al. [29], based on Pescher et al. [63] and respecting BM2013-8 ethical clearance from Institute of Tropical Medicine (ITM) Animal Ethics Committee. Briefly, 5-week-old female golden hamsters were infected via intracardiac injection of 5 × 105 stationary phase promastigotes. Three months post-infection, hamsters were euthanized, and amastigotes were purified from the liver. The liver was homogenized in HOMEM, then cleared by centrifugation (130× g for 5 at room temperature). The supernatant was collected, and saponin lysed (1 mg × mL−1) under agitation for 5 min, then centrifuged for 10 min at 1800× g. The pellet was washed three times with ice-cold PBS. Remaining host cells were removed by Percoll centrifugation (GE Healthcare, Chicago, IL, USA) in a 15 mL Falcon tube: the parasites were resuspended in a 45% Percoll solution layered on a 90% Percoll cushion. The gradient interface was recovered after 30 min centrifugation at 3500× g, 15 °C, then washed 3 times with ice-cold PBS (1800× g, 10 min, 15 °C) prior to DNA extraction. T. brucei gambiense MBA bloodstream forms were obtained from OF-1 mice when the parasitemia was at its highest, according to ITM Animal Ethics Committee decision BM2013-7 (26/02/2013). Parasites were separated from the whole blood, as described in Tihon et al. [64]. Briefly, the parasites were separated from the blood by placing the whole blood on an anion exchanger diethylaminoethyl (DEAE)-cellulose resin (Whatman, Maidstone, UK) suspended in phosphate saline glucose (PSG) buffer, pH 8. After elution and two washes on PSG, DNA was extracted. DNA of L. donovani, both promastigotes and amastigotes, as well as T. brucei, was extracted using DNeasy Blood and Tissue Kit (Qiagen, Hilden, Germany) according to manufacturer instructions.
Arabidopsis thaliana Col-0 was grown for 21 days under long-day conditions, i.e., 16 h light and 8 h darkness. DNA was then extracted from the whole rosette leaves using the DNeasy Plant Mini Kit (Qiagen).

2.3. Genetic Engineering of L. donovani BPK282

We generated both an LdDNMT overexpressing (LdDNMT+) and the null mutant line (LdDNMT-/-) of L. donovani BPK282. All the PCR products generated to produce the constructs for LdDNMToverex and LdDNMTKO were sequenced at the VIB sequencing facility using the same primer as for the amplification. For LdDNMToverex, the overexpression construct, pLEXSY-DNMT, was generated by PCR amplification of LdBPK_251230 from BPK282 genomic DNA using Phusion (NEB, Ipswich, Massachusetts, US) and cloned inside the expression vector pLEXSY-Hyg2 (JENA Bioscience, Jena, Germany) using NEBuilder (NEB) according to manufacturer’s instruction for primer design and cloning instructions (Supplementary Table for primers list). Once generated, 10 µg of pLEXSY-DNMT was electroporated in 5 × 107 BPK282 promastigotes from logarithmic culture using cytomix on a GenePulserX (BioRad, Hercules, California, US), according to LeBowitz (1994) [65], and selected in vitro by adding 50 μg/mL hygromycin B (JENA Bioscience) until parasite growth [66]. Verification of overexpression was carried out by qPCR on a LightCycler480 (Roche, Basel, Switzerland) using SensiMix SYBR No-ROX (Bioline, London, UK) on cDNA. Briefly, 108 logarithmic-phase promastigotes were pelleted, RNA extraction was performed using RNAqueous-Micro Total RNA Isolation Kit (Thermo Fisher Scientific) and quantified by Qubit, and the Qubit RNA BR Assay (Thermo Fisher ScientificTranscriptor reverse transcriptase (Roche) was used to synthesize cDNA following manufacturer’s instructions. qPCRs were run on a LightCycler 480 (Roche) with a SensiMix SYBR No-ROX Kit (Bioline); primer sequences are available in Supplementary Table S1. Normalization was performed using two transcripts, previously described as stable in promastigotes and amastigotes by Dumetz et al. (2018) [67]—LdBPK_340035000 and LdBPK_240021200.
For the generation of LdDNMT-/-, a two-step gene replacement strategy was used: replacing the first allele of LdBPK_251230 by nourseothricin resistance gene (SAT) and the second allele by a puromycin resistance gene (Puro). Briefly, each drug resistance gene was PCR amplified from pCL3S and pCL3P using Phusion (NEB) and cloned between 300 bp of PCR amplified DNA fragments from the upstream and downstream region of LdBPK_251230 using NEBuilder (NEB) inside pUC19 for construct amplification in E. coli DH5α (Promega, Madison, WI, USA) (cf. primer list in Supplementary Table S1). Each replacement construct was excised from pUC19 using SmaI (NEB), dephosphorylated using Antarctic phosphatase (NEB), and 10 µg of DNA was used for the electroporation in the same conditions, as previously described, to insert the pLEXSY-DNMT. The knock-out was confirmed by whole-genome sequencing.

2.4. Bisulfite Sequencing and Data Analysis

For each sample, one microgram of genomic DNA was used for bisulfite conversion with innuCONVERT Bisulfite All-In-One Kit (Analytikjena, Jena, Germany). Sequencing libraries were prepared with the TruSeq DNA Methylation Kit according to the manufacturer’s instructions (Illumina, San Diego, CA, USA). The resulting libraries were paired-end (2 × 100 bp) and sequenced on the Illumina HiSeq 1500 platform of the University of Antwerp (Centre of Medical Genetics). The sequencing quality was first verified with FastQC v0.11.4. Raw reads generated for each sample were aligned to their respective reference genome with BSseeker 2-2.0.3 [68]: LdBPK282v2 [29] for L. donovani, TREU927 [34] for T. brucei, and Tair10 [69] for the A. thaliana positive control. Samtools fixmate (option -m) and samtools markdup (option -r) were then used to remove duplicate reads. CpG, CHG, and CHH methylation sites were subsequently called with the BS-Seeker2 ‘call’ tool using default settings and further filtered with our Python3 workflow called ‘Bisulfilter’ (available at https://github.com/CuypersBart/Bisulfilter). Genome-wide visualization of methylated regions was then carried out with ggplot2 in R [70]. In Leishmania, the positions that passed our detection thresholds (coverage >25, methylation percentage >0.8) were then manually inspected in IGV 2.5.0 [71].

2.5. Leishmania DNA Enrichment from A Mix of Human and Leishmania DNA

To check whether the lack of detectable DNA methylation in Leishmania can be used for the enrichment of Leishmania vs. (methylated) human DNA, we carried out methylated DNA removal on two types of samples: (1) An artificial mix of L. donovani BPK282/0 cl4 promastigote DNA with human DNA (Promega) from 1/15 to 1/150,000 (Leishmania:human) and (2) Linked promastigote and hamster-derived amastigote samples from 3 clinical Leishmania donovani strains (BPK275, BPK282, and BPK026), which were generated in previous work [29]. For this experiment, we used a 1/1500 artificial mix of promastigote DNA and human DNA (Promega) to reflect the median ratio found in clinical samples. For each of the three biological replicates (strains), we carried out the experiment in duplicate (technical replicates). All parasite DNA was extracted with the DNA (DNeasy Blood and Tissue Kit, Qiagen). Leishmania DNA (0.017 ng/μL) was then enriched from the human DNA (25 ng/μL) using NEBNext Microbiome DNA Enrichment Kit (NEB) according to manufacturer instructions. Evaluation of the ratio Leishmania/human DNA was performed by qPCR on LightCycler480 (Roche) using SensiMix SYBR No-ROX (Bioline) and RPL30 primers provided in the kit to measure human DNA and Leishmania CS primers (cysteine synthase) [72]. Cysteine synthase is a single copy gene localized on chromosome 36. Across all the different studies carried on L. donovani, chromosome 36 is shown to be disomic in all strains (Imamura et al. 2016). Leishmania DNA quantification using qPCR and targeting the CS gene was previously validated by sequencing by Domagalska et al. [73].

3. Results

3.1. The Leishmania Genome Contains a Putative C-5 DNA Methyltransferase (DNMT)

Eukaryotic DNA methylation typically requires the presence of a functional C-5 cytosine-specific DNA methylase (C5 DNA MTase). This type of enzymes specifically methylates the C-5 position of cytosines in DNA, using S-adenosyl methionine as a methyl-donor. To check for the presence of C5 DNA MTases in Leishmania donovani, we carried out a deep search of the parasite’s genome. In particular, we used the LdBPKv2 reference genome [29] and searched the predicted protein sequences of this assembly using the hidden Markov model (hmm) signature of the C5 DNA MTase protein family obtained from PFAM (PF00145) and obtained a single hit: the protein LdBPK_251230 (E-value: 2.7 × 10−40). LdBPK_251230 was already annotated as ‘modification methylase-like protein’ with a predicted length of 840 amino acids. In the following text, this protein has been referred to as LdDNMT. Interestingly, in another trypanosomatid species, Trypanosoma brucei, the homolog of this protein (Tb927.3.1360 or TbDNMT) has been previously studied in detail by Militello et al. [27]. Moreover, these authors showed that TbDNMT has all ten conserved domains that are present in functional DNMTs. We aligned TbDNMT with LdDNMT using T-Coffee (Figure 1) and found that these 10 domains are also present in LdDNMT, including the putative catalytic cysteine residue in domain IV.

3.2. Leishmania and Trypanosomatid C-5 DNA Methyltransferase Belongs to the Eukaryotic DNMT6 Family

To learn more about the putative function and evolutionary history of this protein, we wanted to characterize the position of LdDNMT and those of related trypanosomatid species within the DNMT phylogenetic tree. Consequently, we collected the publicly available, putative proteomes of a wide range of prokaryotic and eukaryotic species, searched them for the hmm signature of the C5-DNMT family, aligned the identified proteins, and generated a RAxML maximum likelihood tree. In total, we identified 211 putative family members in the genomes of 37 species (E-value < 0.01), including seven prokaryotic (Agrobacterium tumefaciens, Clostridium botulinum, Escherichia coli, Eubacterium rectale, Nostoc punctiforme, Salmonella enterica, and Treponema succinifaciens) and 30 eukaryotic species. These eukaryotic species were selected to contain organisms from the excavata phylum (of which Leishmania is part) and a range of other, often better-characterized phyla as a reference. The excavata species included 11 trypanosomatids (Trypanosoma grayi, Trypanosoma congolense, Trypanosoma vivax, Trypanosoma brucei, Trypanosoma cruzi, Leptomonas seymouri, Leishmania braziliensis, Leishmania tarentolae, Leishmania amazonensis, Leishmania major, and Leishmania donovani), 1 other non-trypanosomatid euglenozoid species (Euglena gracilis), and 1 other non-euglenozoid species (Naegleria gruberi). The other eukaryotic phyla included in the analysis were apicomplexa (Plasmodium vivax, Plasmodium falciparum, Cryptosporidium parvum, Cryptosporidium hominis), amoebozoa (Entamoeba histolytica), angiosperma (Arabidopsis thaliana, Oryza sativa), ascomycota (Ascobolus immerses, Aspergillus fumigatus, Aspergillus nidulans, Neurospora crassa), chlorophyta (Micromonas pusilla), chordata (Homo Sapiens, Danio rerio), haptophyta (Emiliania huxleyi), ochrophyta (Phaeodactylum tricornutum, Thalassiosira pseudonana) (Figure 2).
Our phylogenetic tree was able to clearly separate known DNMT subgroups, including DNMT1, DNMT2, DNMT3, DNMT4, DNMT5, and several groups of prokaryotic DNMTs [1,16,74]. Interestingly, the tree also showed that most trypanosomatids have exactly one DNMT, which is part of the much less-characterized DNMT6 group, as has been previously described for Leishmania major and Trypanosoma brucei [20]. This group of DNMTs has also been found in chlorophyta, haptophyta, ochrophyta (in our tree represented as Micromonas pusilla, Emiliania huxleyi, and Phaeodactylum tricornutum, respectively) and recently in dinoflagellates (e.g., Symbiodinium kawagutii and Symbiodinium minutum), but its function remains elusive [20,75]. Next to this array of unicellular eukaryotes, the most closely related branch to DNMT6 contains a series of bacterial DNMTs (in our tree represented by Agrobacterium tumefaciens, Eubacterium rectale, Escherichia coli, Salmonella enterica, and Treponema succinifaciens). This suggests that DNMT6 already diversified from the other DNMT groups within the pool of prokaryotic DNMTs. Interestingly, from our tree, the DNMT composition of other excavata species appeared to be very different. The euglenozoid, Euglena gracilis, has DNMT1, DNMT2, DNMT4, and DNMT5, while Naegleria gruberi has both a DNMT1 and a DNMT2. In contrast, we found that the majority of trypanosomatids have just one single DNMT6. This complex distribution of DNMT groups over the different excavata species is suggestive of extensive horizontal gene transfer, which is regarded upon as one of the main mechanisms for the acquisition of new MTases in eukaryotes [20].
Interestingly, Trypanosoma cruzi was the only trypanosomatid with more than one DNMT, having two proteins significantly match the hmm. These two DNMTs are not part of the DNMT6 cluster, but instead had no close relatives in the tree (the closest being prokaryotic DNMTs). Importantly, the protein database used to identify these two DNMTs originated from a six-frame translation of the T. cruzi Bug2148 draft assembly. Therefore, further improvement of this draft proteome is needed to confirm or reject our findings.

3.3. Whole-Genome Bisulfite Sequencing Reveals No Evidence for Functional C-5 Methylation

As (1) we identified LdBPK_251230 to be from the C5 DNA MTase family, (2) all 10 conserved domains were present. We decided to check for the presence and functional role of C5 DNA methylation in L. donovani. Therefore, we assessed the locations and degree of CpG, CHG, and CHH methylation across the entire Leishmania genome and within the two parasite life stages: amastigotes (intracellular mammalian life stage) and promastigotes (extracellular, insect life stage). Amastigotes were derived directly from infected hamsters, while promastigotes were obtained from axenic cultures. Promastigotes were divided into two batches, one passaged long-term in axenic culture, the other passaged once through a hamster and then sequenced at axenic passage 3; thus, allowing us to study also the effect of long vs. short term in vitro passaging. Arabidopsis thaliana and T. brucei were included as a positive control as the degree of CpG, CHG, and CHH methylation in A. thaliana is well known [76,77], while T. brucei is the only trypanosomatid in which (low) methylation levels were previously detected by mass spectrometry [27].
An overview of all sequenced samples can be found in Supplementary Table S2. All L. donovani samples were sequenced with at least 30 million 100 bp paired-end (PE) reads (60 million total) per sample, resulting in an average genomic coverage of at least 94X for the Leishmania samples. The T. brucei was sequenced with 69 million PE reads, resulting in 171X average coverage, and A. thaliana 27 million PE reads, resulting in 21X average coverage. Detailed mapping statistics can be found in Supplementary Table S3.
We first checked for global methylation patterns across the genome. Interestingly, we could not detect any methylated regions in Leishmania donovani promastigotes, both short (P3) and long-term in vitro passaged, nor in hamster-derived promastigotes or amastigotes (Figure 3). Minor increases in the CHH signal towards the start end of several chromosomes were manually checked in IGV and attributed to poor mapping in (low complexity) telomeric regions. This was in contrast to our positive control, Arabidopsis thaliana, which showed clear, highly methylated CpG, CHG, and CHH patterns across the genome. This distribution was consistent with prior results with MethGO observed on Arabidopsis thaliana, confirming that our methylation detection workflow was working [70].
In a second phase, we checked for individual sites that were fully methylated (>80% of the sequenced DNA at that site) using BS-Seeker2 and filtering the results with our automated Python3 workflow. CpG methylation in all three biological samples for L. donovani was lower than 0.0003%, CHG methylation lower than 0.0005%, and CHH methylation lower than 0.0126% (Table 1, Supplementary Table S4). However, when this low number of detected ‘methylated’ sites was manually verified in IGV, they could all clearly be attributed to regions where BS-Seeker2 wrongly called methylated bases, either because of poor mapping (often in repetitive, low complexity regions) or of strand biases. In reliably mapped regions, there was clearly no methylation. Similarly, we detected 0.0001% of CpG methylation, 0.0006% of CHG methylation, and 0.0040% of CHH methylation for T. brucei, which could all be attributed to mapping errors or strand biases. In A. thaliana, our positive control, we detected 21.05% of CpG methylation, 4.04% of CHG methylation, and 0.31% of CHH methylation, which is similar as the reported values in the literature [78,79], and demonstrated that our bioinformatic workflow could accurately detect methylated sites. We also checked sites with a lower methylation degree (>40%), which gave higher percentages, but this could be attributed to the increased noise level at this resolution (Supplementary Table S4). Indeed, even when applying stringent coverage criteria (>25×), this approach is susceptible to false-positive methylation calls, as we are checking millions of positions (in case of Leishmania, more than 5.8 million CG sites, 3.9 million CHG sites, and 9.3 million CHH sites).
To determine whether LdDNMT is essential, we generated an LdDNMT knock-out line (LdDNMT-/-) and validated it by calculating its LdDNMT copy number based on the sequencing coverage (Figure 4). Indeed, the copy number of the LdDNMT gene in LdDNMT-/- was reduced to zero. The fact that the LdDNMT-/- line was viable shows that LdDNMT is not an essential gene in promastigotes. Similar to the wild-type promastigotes and amastigotes described above, no CpG, CHG, or CHH methylation could be observed with bisulfite sequencing (Table 1, Supplementary Table S4). Vice versa, to check if overexpressing LdDNMT would affect the C-5 DNA-methylation patterns, we generated an overexpressor line (LdDNMT+). In this line, the LdDNMT genomic copy number was increased to 78 copies, and mRNA levels (Table 2) showed a 2.5-fold higher expression than the corresponding wild type. Although the LdDNMT+ initially seemed to have slightly higher methylation percentages (Table 1, Supplementary Table S4), none of the CpG, CHG, or CHH showed to be methylated after our manual validation in IGV, just as in the wild type lines and LdDNMT-/-. This suggested again that LdDNMT had no impact on DNA methylation, although we did not validate LdDNMT overexpression in LdDNMT+ on the protein level. Finally, neither the deletion nor the overexpression of LdDNMT had an obvious effect on the growth of the parasites, although this was not tested in a systematic way.

3.4. Absence of C5 DNA Methylation as a Leishmania vs. Host DNA Enrichment Strategy

The lack or low level of C5 DNA methylation opens the perspective for enriching Leishmania DNA in mixed parasite-host DNA samples, based on the difference in methylation status (the vertebrate host does show C5 DNA methylation). This could potentially be an interesting pre-enrichment step before the whole-genome sequencing analysis of clinical samples containing Leishmania. Furthermore, commercial kits for removing methylated DNA are readily available and typically contain a methyl-CpG-binding domain (MBD) column, which binds methylated DNA while allowing unmethylated DNA to flow through.
To test if these kits can be used for the relative enrichment of Leishmania DNA, we first generated artificially mixed samples using different ratios of L. donovani promastigote DNA with human DNA. Ratios were made starting from 1/15 to 1/15,000, which reflects the real ratio of Leishmania vs. human DNA in clinical samples [80]. From these mixes, Leishmania DNA was enriched using NEBNext Microbiome DNA Enrichment Kit (NEB) that specifically binds methylated DNA, while the non-methylated remains in the supernatant. Leishmania vs. human DNA ratios were determined by qPCR, targeting the RPL30 primers provided in the kit to measure human DNA and cysteine synthase for Leishmania. We observed an average of 263× enrichment of Leishmania vs. human DNA (Figure 5). This ranged between 378× for the lowest dilution (removing 99.8% of the human DNA) to 164× (removing 99.6% of the human DNA) in the highest diluted condition (1/15,000 Leishmania:human).
Secondly, we wanted to test if enrichment via MBD columns worked equally well on L. donovani amastigotes for (a) fundamental reasons, as the (indirect) second method to detect if there are any methylation differences between promastigotes and amastigotes, and (b) practical reasons, as it is the (intracellular) life stage encountered in clinical samples. Therefore, we also carried out this enrichment technique on three sets (three strains) of hamster-derived amastigotes and their promastigote controls. Similarly, as in the previous experiment, Leishmania-human DNA mixes were generated in a 1/1500 (Leishmania:human) ratio, after which enrichment was carried out with the NEBNext Microbiome DNA Enrichment Kit. The enrichment worked well for both life stages, the promastigote samples were, on average, 76.22 ± 14.28 times enriched, and the amastigote samples 61.68 ± 4.23 times (Table 3).

4. Discussion

With this work, we present the first comprehensive study addressing the status of DNA-methylation in Leishmania.
We demonstrated that the Leishmania genome contains a C5-DNMT (LdDNMT) that contains all 10 conserved DNMT domains. We also showed the gene is expressed at the RNA level. As the C5-DNMT family is diverse, and several family members are known to have adopted (partially) distinct functions during the course of evolution, we were particularly interested in the position of this DNMT within the evolutionary tree of this family, as it could direct hypotheses about the function of this protein. We found that LdDNMT is, in fact, a DNMT6 and that all studied Leishmania species and, more generally, the majority of trypanosomatids have exactly one copy of this DNMT6 in their genome (exc. T. cruzi) [20]. Interestingly, all other (non-trypanosomatid) species studied so far have either multiple DNMT6 copies and/or other DNMT subfamily members in their genomes [20,75]. Therefore, trypanosomatids might be a unique model species to further study the role of this elusive DNMT subfamily, as there can be no interaction with the effects of other DNMTs.
The fact that our LdDNMT knock-out line (verified by sequencing) was viable shows that DNMT6 is not essential for the survival of the parasite, at least in promastigotes and in our experimental conditions. However, at the same time, one might hypothesize that DNMT6 does offer a selective advantage to the parasite. First of all, the sequence of DNMT core domains is extremely conserved across the tree of life, and this is no different from those that we encountered in Leishmania. Secondly, Leishmania is characterized by high genome plasticity and features extensive gene copy number differences between strains [81,82]. Therefore, one might speculate that the parasite would have lost the gene a long time ago if it did not provide any selective advantage.
In addition, we aimed to characterize the DNA methylation patterns of the parasite’s genome. Therefore, we carried out the first multi-life stage whole-genome bisulfite sequencing experiment on Leishmania and trypanosomatids in general. We checked both the promastigote (both culture and amastigote life stage). Surprisingly, we did not find any evidence for DNA methylation in L. donovani even though we checked both for large, regional patterns (sensitive for low levels of methylation over longer distances), and site-specific analyses (sensitive for high levels of methylation at individual sites). This could either mean that there is indeed no DNA methylation in these species or that was below our detection threshold. Regarding this detection threshold, two factors should be considered. Firstly, bisulfite sequencing and analysis allows for the detection of specific sites that are consistently methylated across the genomes of a mix of cells. In our case, we looked for sites that are methylated in at least 80% or 40% of the cases. Thus, if Leishmania consistently methylates certain genomic positions, our pipeline would have uncovered this. However, if this methylation would be more random, or occurring in only a small subset of cells, we would not be able to distinguish this for random sequencing errors, and as such, we cannot exclude this possibility. Secondly, bisulfite sequencing typically suffers from poor genomic coverage due to the harsh BS treatment of the DNA [83]. In our L. donovani samples, we covered at least 30.14% of the CpG sites, 29.47% of the CHG sites, and 24.23% of the CHH sites (even though having more than 90× average coverage). However, as there are millions of CpG, CHG, and CHH sites in the genome, the chance is very small (0.75n, with n = number of methylated sites) that we would not have detected methylated sites, even if present in low numbers.
In any case, it is hard to imagine that any of the typical eukaryotic DNA methylation systems, such as genomic imprinting, chromosome inactivation, gene expression regulation, and/or the repression of transposable elements, could be of significance with such low methylation levels. On the other hand, given its phylogenetic position, it is perfectly possible that DNMT6 has changed its biological activity and now carries out another function. Indeed, as we described above, a similar phenomenon was observed with DNMT2 that switched its substrate from DNA to tRNA during the course of evolution [16,17].
Correspondingly, we did not observe any detectable DNA methylation for T. brucei. These findings are in contrast to what has been reported before by Militello et al., who detected 0.01% of 5MC in the T. brucei genome [27]. Besides, the methylated (orthologous) loci described in this paper could not be confirmed in the current work. However, this is maybe not surprising as the same authors reported later that TbDNMT might, in fact, methylate RNA, as they identified methylated sites in several tRNAs [71]. This would indeed explain why we did not observe C5-DNA methylation in T. brucei with high resolution, whole-genome bisulfite sequencing, and further suggest that a similar substrate switch to tRNA has occurred for DNMT6, just like has occurred for DNMT2. Further functional characterization of DNMT6 is required to verify this hypothesis.
From an applied perspective, this study opens new avenues for the enrichment of trypanosomatid DNA from clinical samples, which often have an abundance of host DNA. Indeed, depletion of methylated DNA could be included as a pre-enrichment step for existing enrichment approaches. For example, our group has recently obtained excellent sequencing results of clinical samples using SureSelect (97% of the samples for diagnostic SNPs, 83% for genome-wide information for sequenced samples), but was not able to sequence samples below 0.006% of Leishmania DNA content [73]. Perhaps the removal of methylated DNA could further enhance the sensitivity of this method. In the case of Leishmania, the technique could be useful for both enrichments from mammalian hosts and the insect vector, as it has been recently shown that the phlebotomine vector also carries Me5C in its genome [84]. The depletion of methylated DNA as a pre-enrichment step before whole-genome sequencing has also been successfully used before for the parasite Plasmodium falciparum (malaria) and shown to generate unbiased sequencing reads [85].
In conclusion, we demonstrated that the Leishmania genome encodes for a DNMT6, but DNA methylation is either absent or present in such a low proportion that it is unlikely to have a major functional role. Instead, we suggest that more investigation at the RNA level is required to address the function of DNMT6 in Leishmania. The absence of DNA methylation provides a new working tool for the enrichment of Leishmania DNA in clinical samples, thus facilitating future parasitological studies.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-2607/8/8/1252/s1; Table S1: Overview of all primers used in this work, Table S2: Overview of bisulfite converted and sequenced samples in this study with SRA accession numbers, Table S3: Sequencing statistics (raw reads, mapped reads and coverage per sample) for all bisulfite sequenced samples, Table S4: Methylation call statistics per bisulfite sequenced sample. TableRaw sequencing data is available in the Sequence Read Archive under project accession numbers PRJNA560731 and PRJNA560871.

Author Contributions

Designed the experiments: B.C., F.D., G.D.M., J.-C.D., M.A.D. Performed the experiments: B.C., F.D., M.A.D. Analyzed the data: B.C., F.D., P.M., K.L., M.A.D. Wrote the manuscript: B.C., F.D., P.M., K.L., J.-C.D., M.A.D. All authors have read and agreed to the published version of the manuscript.

Funding

The computational resources and services used in this work were provided by the VSC (Flemish Supercomputer Center), funded by the Research Foundation—Flanders (FWO) and the Flemish Government—department EWI. This work was also supported by the Department of Economy, Science, and Innovation in Flanders ITM-SOFIB (SINGLE project, to J.-C.D.). We thank the Center of Medical Genetics at the University of Antwerp for hosting the NGS facility. B.C. is a post-doctoral fellow funded from the FWO [12V5319N].

Acknowledgments

We thank G.Z. and G.B. for providing us with the A. thaliana DNA, as well as P.B. and N.B. for the bloodstream form of T. brucei gambiense MBA and J.C. for the Leishmania expression vectors pCL3S and pCL3P. This work was supported by the Interuniversity Attraction Poles Program of Belgian Science Policy [P7/41 to J.-C.D.] and by the organization “Les amis des Instituts Pasteur à Bruxelles, asbl” [F.D.].

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jurkowski, T.P.; Jeltsch, A. On the evolutionary origin of eukaryotic DNA methyltransferases and Dnmt2. PLoS ONE 2011, 6, e28104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Li, E.; Zhang, Y. DNA methylation in mammals. Cold Spring Harb Perspect. Biol. 2014, 6, a019133. [Google Scholar] [CrossRef] [PubMed]
  3. Vidal, E.; Sayols, S.; Moran, S.; Guillaumet-Adkins, A.; Schroeder, M.P.; Royo, R.; Orozco, M.; Gut, M.; Gut, I.; Lopez-Bigas, N.; et al. A DNA methylation map of human cancer at single base-pair resolution. Oncogene 2017, 36, 5648. [Google Scholar] [CrossRef]
  4. Jin, Z.; Liu, Y. DNA methylation in human diseases. Genes Dis. 2018, 5, 1–8. [Google Scholar] [CrossRef] [PubMed]
  5. Guo, W.; Chung, W.-Y.; Qian, M.; Pellegrini, M.; Zhang, M.Q. Characterizing the strand-specific distribution of non-CpG methylation in human pluripotent cells. Nucleic Acids Res. 2013, 42, 3009–3016. [Google Scholar] [CrossRef] [Green Version]
  6. Lister, R.; Pelizzola, M.; Dowen, R.H.; Hawkins, R.D.; Hon, G.; Tonti-Filippini, J.; Nery, J.R.; Lee, L.; Ye, Z.; Ngo, Q.M.; et al. Human DNA methylomes at base resolution show widespread epigenomic differences. Nature 2009, 462, 315–322. [Google Scholar] [CrossRef] [Green Version]
  7. Tomizawa, S.; Kobayashi, H.; Watanabe, T.; Andrews, S.; Hata, K.; Kelsey, G.; Sasaki, H. Dynamic stage-specific changes in imprinted differentially methylated regions during early mammalian development and prevalence of non-CpG methylation in oocytes. Development 2011, 138, 811–820. [Google Scholar] [CrossRef] [Green Version]
  8. Xie, W.; Barr, C.L.; Kim, A.; Yue, F.; Lee, A.Y.; Eubanks, J.; Dempster, E.L.; Ren, B. Base-resolution analyses of sequence and parent-of-origin dependent DNA methylation in the mouse genome. Cell 2012, 148, 816–831. [Google Scholar] [CrossRef] [Green Version]
  9. Chédin, F. Chapter 7-The DNMT3 family of mammalian de novo dna methyltransferases. In Progress in Molecular Biology and Translational Science; Cheng, X., Blumenthal, R.M., Eds.; Academic Press: Cambridge, MA, USA, 2011; Volume 101, pp. 255–285. [Google Scholar]
  10. Jones, P.A.; Liang, G. Rethinking how DNA methylation patterns are maintained. Nat. Rev. Genet. 2009, 10, 805–811. [Google Scholar] [CrossRef] [Green Version]
  11. Jurkowski, T.P.; Meusburger, M.; Phalke, S.; Helm, M.; Nellen, W.; Reuter, G.; Jeltsch, A. Human DNMT2 methylates tRNA(Asp) molecules using a DNA methyltransferase-like catalytic mechanism. RNA 2008, 14, 1663–1670. [Google Scholar] [CrossRef] [Green Version]
  12. Tuorto, F.; Liebers, R.; Musch, T.; Schaefer, M.; Hofmann, S.; Kellner, S.; Frye, M.; Helm, M.; Stoecklin, G.; Lyko, F. RNA cytosine methylation by Dnmt2 and NSun2 promotes tRNA stability and protein synthesis. Nat. Struct. Mol. Biol. 2012, 19, 900–905. [Google Scholar] [CrossRef]
  13. Burgess, A.L.; David, R.; Searle, I.R. Conservation of tRNA and rRNA 5-methylcytosine in the kingdom Plantae. BMC Plant Biol. 2015, 15, 199. [Google Scholar] [CrossRef] [Green Version]
  14. Schaefer, M.; Pollex, T.; Hanna, K.; Lyko, F. RNA cytosine methylation analysis by bisulfite sequencing. Nucleic Acids Res. 2009, 37, e12. [Google Scholar] [CrossRef] [Green Version]
  15. Jeltsch, A.; Ehrenhofer-Murray, A.; Jurkowski, T.P.; Lyko, F.; Reuter, G.; Ankri, S.; Nellen, W.; Schaefer, M.; Helm, M. Mechanism and biological role of Dnmt2 in Nucleic Acid Methylation. RNA Biol. 2016, 14, 1108–1123. [Google Scholar] [CrossRef] [Green Version]
  16. Ponts, N.; Fu, L.; Harris, E.Y.; Zhang, J.; Chung, D.-W.D.; Cervantes, M.C.; Prudhomme, J.; Atanasova-Penichon, V.; Zehraoui, E.; Bunnik, E.M.; et al. Genome-wide mapping of DNA methylation in the human malaria parasite Plasmodium falciparum. Cell Host Microbe 2013, 14, 696–706. [Google Scholar] [CrossRef] [Green Version]
  17. Motorin, Y.; Grosjean, H. Multisite-specific tRNA:m5C-methyltransferase (Trm4) in yeast Saccharomyces cerevisiae: Identification of the gene and substrate specificity of the enzyme. RNA 1999, 5, 1105–1118. [Google Scholar] [CrossRef] [Green Version]
  18. Jeltsch, A.; Nellen, W.; Lyko, F. Two substrates are better than one: Dual specificities for Dnmt2 methyltransferases. Trends Biochem. Sci. 2006, 31, 306–308. [Google Scholar] [CrossRef]
  19. Kaiser, S.; Jurkowski, T.P.; Kellner, S.; Schneider, D.; Jeltsch, A.; Helm, M. The RNA methyltransferase Dnmt2 methylates DNA in the structural context of a tRNA. RNA Biol. 2016, 14, 1241–1251. [Google Scholar] [CrossRef]
  20. Ponger, L.; Li, W.H. Evolutionary diversification of DNA methyltransferases in eukaryotic genomes. Mol. Biol. Evol. 2005, 22, 1119–1128. [Google Scholar] [CrossRef] [Green Version]
  21. Wei, H.; Jiang, S.; Chen, L.; He, C.; Wu, S.; Peng, H. Characterization of cytosine methylation and the DNA methyltransferases of toxoplasma gondii. Int. J. Biol. Sci. 2017, 13, 458–470. [Google Scholar] [CrossRef] [Green Version]
  22. Morselli, M.; Pastor, W.A.; Montanini, B.; Nee, K.; Ferrari, R.; Fu, K.; Bonora, G.; Rubbi, L.; Clark, A.T.; Ottonello, S.; et al. In vivo targeting of de novo DNA methylation by histone modifications in yeast and mouse. eLife 2015, 4, e06205. [Google Scholar] [CrossRef]
  23. Greer, E.L.; Blanco, M.A.; Gu, L.; Sendinc, E.; Liu, J.; Aristizábal-Corrales, D.; Hsu, C.-H.; Aravind, L.; He, C.; Shi, Y. DNA methylation on N6-Adenine in C. elegans. Cell 2015, 161, 868–878. [Google Scholar] [CrossRef] [Green Version]
  24. Ivens, A.C.; Peacock, C.S.; Worthey, E.A.; Murphy, L.; Aggarwal, G.; Berriman, M.; Sisk, E.; Rajandream, M.A.; Adlem, E.; Aert, R.; et al. The genome of the kinetoplastid parasite, Leishmania major. Science 2005, 309, 436–442. [Google Scholar] [CrossRef] [Green Version]
  25. Van Luenen, H.G.A.M.; Farris, C.; Jan, S.; Genest, P.-A.; Tripathi, P.; Velds, A.; Kerkhoven, R.M.; Nieuwland, M.; Haydock, A.; Ramasamy, G.; et al. Glucosylated hydroxymethyluracil, DNA base J, prevents transcriptional readthrough in Leishmania. Cell 2012, 150, 909–921. [Google Scholar] [CrossRef] [Green Version]
  26. Huff, J.T.; Zilberman, D. Dnmt1-independent CG methylation contributes to nucleosome positioning in diverse eukaryotes. Cell 2014, 156, 1286–1297. [Google Scholar] [CrossRef] [Green Version]
  27. Militello, K.T.; Wang, P.; Jayakar, S.K.; Pietrasik, R.L.; Dupont, C.D.; Dodd, K.; King, A.M.; Valenti, P.R. African trypanosomes contain 5-Methylcytosine in nuclear DNA. Eukaryot. Cell 2008, 7, 2012–2016. [Google Scholar] [CrossRef] [Green Version]
  28. Bateman, A.; Smart, A.; Luciani, A.; Salazar, G.A.; Mistry, J.; Richardson, L.J.; Qureshi, M.; El-Gebali, S.; Potter, S.C.; Finn, R.D.; et al. The Pfam protein families database in 2019. Nucleic Acids Res. 2018, 47, D427–D432. [Google Scholar] [CrossRef]
  29. Dumetz, F.; Imamura, H.; Sanders, M.; Seblova, V.; Myskova, J.; Pescher, P.; Vanaerschot, M.; Meehan, C.J.; Cuypers, B.; De Muylder, G.; et al. Modulation of aneuploidy in leishmania donovani during adaptation to different in vitro and in vivo environments and its impact on gene expression. MBio 2017, 8, e00599-17. [Google Scholar] [CrossRef] [Green Version]
  30. Real, F.; Vidal, R.O.; Carazzolle, M.F.; Mondego, J.M.C.; Costa, G.G.L.; Herai, R.H.; Würtele, M.; de Carvalho, L.M.; Carmona e Ferreira, R.; Mortara, R.A.; et al. The genome sequence of Leishmania (Leishmania) amazonensis: Functional annotation and extended analysis of gene models. DNA Res. 2013, 20, 567–581. [Google Scholar] [CrossRef] [Green Version]
  31. Peacock, C.S.; Seeger, K.; Harris, D.; Murphy, L.; Ruiz, J.C.; Quail, M.A.; Peters, N.; Adlem, E.; Tivey, A.; Aslett, M.; et al. Comparative genomic analysis of three Leishmania species that cause diverse human disease. Nat. Genet. 2007, 39, 839–847. [Google Scholar] [CrossRef] [Green Version]
  32. Raymond, F.; Boisvert, S.; Roy, G.; Ritt, J.-F.; Légaré, D.; Isnard, A.; Stanke, M.; Olivier, M.; Tremblay, M.J.; Papadopoulou, B.; et al. Genome sequencing of the lizard parasite Leishmania tarentolae reveals loss of genes associated to the intracellular stage of human pathogenic species. Nucleic Acids Res. 2012, 40, 1131–1147. [Google Scholar] [CrossRef] [Green Version]
  33. Kraeva, N.; Butenko, A.; Hlaváčová, J.; Kostygov, A.; Myškova, J.; Grybchuk, D.; Leštinová, T.; Votýpka, J.; Volf, P.; Opperdoes, F.; et al. Leptomonas seymouri: Adaptations to the dixenous life cycle analyzed by genome sequencing, transcriptome profiling and co-infection with leishmania donovani. PLoS Pathog. 2015, 11, e1005127. [Google Scholar] [CrossRef] [Green Version]
  34. Berriman, M.; Ghedin, E.; Hertz-Fowler, C.; Blandin, G.; Renauld, H.; Bartholomeu, D.C.; Lennard, N.J.; Caler, E.; Hamlin, N.E.; Haas, B.; et al. The genome of the African trypanosome Trypanosoma brucei. Science 2005, 309, 416–422. [Google Scholar] [CrossRef] [Green Version]
  35. Abbas, A.H.; Silva Pereira, S.; D’Archivio, S.; Wickstead, B.; Morrison, L.J.; Hall, N.; Hertz-Fowler, C.; Darby, A.C.; Jackson, A.P. The structure of a conserved telomeric region associated with variant antigen loci in the blood parasite trypanosoma congolense. Genome Biol. Evol. 2018, 10, 2458–2473. [Google Scholar] [CrossRef]
  36. Kelly, S.; Ivens, A.; Manna, P.T.; Gibson, W.; Field, M.C. A draft genome for the African crocodilian trypanosome Trypanosoma grayi. Sci. Data 2014, 1, 140024. [Google Scholar] [CrossRef] [Green Version]
  37. Aslett, M.; Aurrecoechea, C.; Berriman, M.; Brestelli, J.; Brunk, B.P.; Carrington, M.; Depledge, D.P.; Fischer, S.; Gajria, B.; Gao, X.; et al. TriTrypDB: A functional genomic resource for the Trypanosomatidae. Nucleic Acids Res. 2010, 38, D457–D462. [Google Scholar] [CrossRef] [Green Version]
  38. Aurrecoechea, C.; Brestelli, J.; Brunk, B.P.; Dommer, J.; Fischer, S.; Gajria, B.; Gao, X.; Gingle, A.; Grant, G.; Harb, O.S.; et al. PlasmoDB: A functional genomic database for malaria parasites. Nucleic Acids Res. 2009, 37, D539–D543. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Auburn, S.; Böhme, U.; Steinbiss, S.; Trimarsanto, H.; Hostetler, J.; Sanders, M.; Gao, Q.; Nosten, F.; Newbold, C.I.; Berriman, M.; et al. A new Plasmodium vivax reference sequence with improved assembly of the subtelomeres reveals an abundance of pir genes. Wellcome Open Res. 2016, 1, 4. [Google Scholar] [CrossRef]
  40. Gardner, M.J.; Hall, N.; Fung, E.; White, O.; Berriman, M.; Hyman, R.W.; Carlton, J.M.; Pain, A.; Nelson, K.E.; Bowman, S.; et al. Genome sequence of the human malaria parasite Plasmodium falciparum. Nature 2002, 419, 498–511. [Google Scholar] [CrossRef]
  41. Puiu, D.; Enomoto, S.; Buck, G.A.; Abrahamsen, M.S.; Kissinger, J.C. CryptoDB: The Cryptosporidium genome resource. Nucleic Acids Res. 2004, 32, D329–D331. [Google Scholar] [CrossRef]
  42. Abrahamsen, M.S.; Templeton, T.J.; Enomoto, S.; Abrahante, J.E.; Zhu, G.; Lancto, C.A.; Deng, M.; Liu, C.; Widmer, G.; Tzipori, S.; et al. Complete genome sequence of the apicomplexan, Cryptosporidium parvum. Science 2004, 304, 441–445. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Xu, P.; Widmer, G.; Wang, Y.; Ozaki, L.S.; Alves, J.M.; Serrano, M.G.; Puiu, D.; Manque, P.; Akiyoshi, D.; Mackey, A.J.; et al. The genome of Cryptosporidium hominis. Nature 2004, 431, 1107–1112. [Google Scholar] [CrossRef]
  44. Gajria, B.; Bahl, A.; Brestelli, J.; Dommer, J.; Fischer, S.; Gao, X.; Heiges, M.; Iodice, J.; Kissinger, J.C.; Mackey, A.J.; et al. ToxoDB: An integrated Toxoplasma gondii database resource. Nucleic Acids Res. 2008, 36, D553–D556. [Google Scholar] [CrossRef] [Green Version]
  45. Lorenzi, H.; Khan, A.; Behnke, M.S.; Namasivayam, S.; Swapna, L.S.; Hadjithomas, M.; Karamycheva, S.; Pinney, D.; Brunk, B.P.; Ajioka, J.W.; et al. Local admixture of amplified and diversified secreted pathogenesis determinants shapes mosaic Toxoplasma gondii genomes. Nat. Commun. 2016, 7, 10147. [Google Scholar] [CrossRef]
  46. Huang, Y.Y.; Cho, S.T.; Lo, W.S.; Wang, Y.C.; Lai, E.M.; Kuo, C.H. Complete genome sequence of agrobacterium tumefaciens Ach5. Genome Announc. 2015, 3, e00570-15. [Google Scholar] [CrossRef] [Green Version]
  47. Galagan, J.E.; Calvo, S.E.; Cuomo, C.; Ma, L.-J.; Wortman, J.R.; Batzoglou, S.; Lee, S.-I.; Baştürkmen, M.; Spevak, C.C.; Clutterbuck, J.; et al. Sequencing of Aspergillus nidulans and comparative analysis with A. fumigatus and A. oryzae. Nature 2005, 438, 1105–1115. [Google Scholar] [CrossRef] [PubMed]
  48. Sebaihia, M.; Peck, M.W.; Minton, N.P.; Thomson, N.R.; Holden, M.T.; Mitchell, W.J.; Carter, A.T.; Bentley, S.D.; Mason, D.R.; Crossman, L.; et al. Genome sequence of a proteolytic (Group I) Clostridium botulinum strain Hall A and comparative analysis of the clostridial genomes. Genome Res. 2007, 17, 1082–1092. [Google Scholar] [CrossRef] [Green Version]
  49. Puerta, M.V.S.; Bachvaroff, T.R.; Delwiche, C.F. The complete plastid genome sequence of the haptophyte emiliania huxleyi: A comparison to other plastid genomes. DNA Res. 2005, 12, 151–156. [Google Scholar] [CrossRef] [PubMed]
  50. Lorenzi, H.A.; Puiu, D.; Miller, J.R.; Brinkac, L.M.; Amedeo, P.; Hall, N.; Caler, E.V. New assembly, reannotation and analysis of the Entamoeba histolytica genome reveal new genomic features and protein content information. PLoS Negl. Trop. Dis. 2010, 4, e716. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Ebenezer, T.E.; Zoltner, M.; Burrell, A.; Nenarokova, A.; Novák Vanclová, A.M.G.; Prasad, B.; Soukal, P.; Santana-Molina, C.; O’Neill, E.; Nankissoor, N.N.; et al. Transcriptome, proteome and draft genome of Euglena gracilis. BMC Biol. 2019, 17, 11. [Google Scholar] [CrossRef] [Green Version]
  52. Worden, A.Z.; Lee, J.-H.; Mock, T.; Rouzé, P.; Simmons, M.P.; Aerts, A.L.; Allen, A.E.; Cuvelier, M.L.; Derelle, E.; Everett, M.V.; et al. Green evolution and dynamic adaptations revealed by genomes of the marine picoeukaryotes micromonas. Science 2009, 324, 268–272. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Galagan, J.E.; Calvo, S.E.; Borkovich, K.A.; Selker, E.U.; Read, N.D.; Jaffe, D.; FitzHugh, W.; Ma, L.J.; Smirnov, S.; Purcell, S.; et al. The genome sequence of the filamentous fungus Neurospora crassa. Nature 2003, 422, 859–868. [Google Scholar] [CrossRef] [PubMed]
  54. Bowler, C.; Allen, A.E.; Badger, J.H.; Grimwood, J.; Jabbari, K.; Kuo, A.; Maheswari, U.; Martens, C.; Maumus, F.; Otillar, R.P.; et al. The Phaeodactylum genome reveals the evolutionary history of diatom genomes. Nature 2008, 456, 239–244. [Google Scholar] [CrossRef] [PubMed]
  55. Parkhill, J.; Dougan, G.; James, K.D.; Thomson, N.R.; Pickard, D.; Wain, J.; Churcher, C.; Mungall, K.L.; Bentley, S.D.; Holden, M.T.; et al. Complete genome sequence of a multiple drug resistant Salmonella enterica serovar Typhi CT18. Nature 2001, 413, 848–852. [Google Scholar] [CrossRef] [PubMed]
  56. Wood, V.; Gwilliam, R.; Rajandream, M.A.; Lyne, M.; Lyne, R.; Stewart, A.; Sgouros, J.; Peat, N.; Hayles, J.; Baker, S.; et al. The genome sequence of Schizosaccharomyces pombe. Nature 2002, 415, 871–880. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Hoskins, J.; Alborn, W.E., Jr.; Arnold, J.; Blaszczak, L.C.; Burgett, S.; DeHoff, B.S.; Estrem, S.T.; Fritz, L.; Fu, D.J.; Fuller, W.; et al. Genome of the bacterium Streptococcus pneumoniae strain R6. J. Bacteriol. 2001, 183, 5709–5717. [Google Scholar] [CrossRef] [Green Version]
  58. Han, C.; Gronow, S.; Teshima, H.; Lapidus, A.; Nolan, M.; Lucas, S.; Hammon, N.; Deshpande, S.; Cheng, J.-F.; Zeytun, A.; et al. Complete genome sequence of Treponema succinifaciens type strain (6091). Stand. Genomic Sci. 2011, 4, 361–370. [Google Scholar] [CrossRef] [Green Version]
  59. Murat, C.; Payen, T.; Noel, B.; Kuo, A.; Morin, E.; Chen, J.; Kohler, A.; Krizsan, K.; Balestrini, R.; Da Silva, C.; et al. Pezizomycetes genomes reveal the molecular basis of ectomycorrhizal truffle lifestyle. Nat. Ecol. Evol. 2018, 2, 1956–1965. [Google Scholar] [CrossRef] [Green Version]
  60. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  61. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [Green Version]
  62. Glez-Peña, D.; Gomez-Blanco, D.; Reboiro-Jato, M.; Fdez-Riverola, F.; Posada, D. ALTER: Program-oriented conversion of DNA and protein alignments. Nucleic Acids Res. 2010, 38, W14–W18. [Google Scholar] [CrossRef] [PubMed]
  63. Pescher, P.; Blisnick, T.; Bastin, P.; Spath, G.F. Quantitative proteome profiling informs on phenotypic traits that adapt Leishmania donovani for axenic and intracellular proliferation. Cell. Microbiol. 2011, 13, 978–991. [Google Scholar] [CrossRef] [PubMed]
  64. Tihon, E.; Imamura, H.; Dujardin, J.C.; Van Den Abbeele, J.; Van den Broeck, F. Discovery and genomic analyses of hybridization between divergent lineages of Trypanosoma congolense, causative agent of Animal African Trypanosomiasis. Mol. Ecol. 2017, 26, 6524–6538. [Google Scholar] [CrossRef] [PubMed]
  65. LeBowitz, J.H. Transfection experiments with Leishmania. Methods Cell Biol. 1994, 45, 65–78. [Google Scholar] [PubMed]
  66. Morales, M.A.; Watanabe, R.; Dacher, M.; Chafey, P.; Osorio y Fortea, J.; Scott, D.A.; Beverley, S.M.; Ommen, G.; Clos, J.; Hem, S.; et al. Phosphoproteome dynamics reveal heat-shock protein complexes specific to the Leishmania donovani infectious stage. Proc. Natl. Acad. Sci. USA 2010, 107, 8381–8386. [Google Scholar] [CrossRef] [Green Version]
  67. Dumetz, F.; Cuypers, B.; Imamura, H.; Zander, D.; D’Haenens, E.; Maes, I.; Domagalska, M.A.; Clos, J.; Dujardin, J.C.; De Muylder, G. Molecular preadaptation to antimony resistance in leishmania donovani on the Indian subcontinent. mSphere 2018, 3, e00548-17. [Google Scholar] [CrossRef] [Green Version]
  68. Guo, W.; Fiziev, P.; Yan, W.; Cokus, S.; Sun, X.; Zhang, M.Q.; Chen, P.-Y.; Pellegrini, M. BS-Seeker2: A versatile aligning pipeline for bisulfite sequencing data. BMC Genom. 2013, 14, 774. [Google Scholar] [CrossRef] [Green Version]
  69. Berardini, T.Z.; Reiser, L.; Li, D.; Mezheritsky, Y.; Muller, R.; Strait, E.; Huala, E. The arabidopsis information resource: Making and mining the “gold standard” annotated reference plant genome. Genesis 2015, 53, 474–485. [Google Scholar] [CrossRef] [Green Version]
  70. Liao, W.W.; Yen, M.R.; Ju, E.; Hsu, F.M.; Lam, L.; Chen, P.Y. MethGo: A comprehensive tool for analyzing whole-genome bisulfite sequencing data. BMC Genomics 2015, 16 (Suppl. S12), 1–8. [Google Scholar] [CrossRef] [Green Version]
  71. Robinson, J.T.; Thorvaldsdóttir, H.; Wenger, A.M.; Zehir, A.; Mesirov, J.P. Variant review with the integrative genomics viewer. Cancer Res. 2017, 77, e31–e34. [Google Scholar] [CrossRef] [Green Version]
  72. Decuypere, S.; Vanaerschot, M.; Rijal, S.; Yardley, V.; Maes, L.; de Doncker, S.; Chappuis, F.; Dujardin, J.C. Gene expression profiling of Leishmania (Leishmania) donovani: Overcoming technical variation and exploiting biological variation. Parasitology 2008, 135, 183–194. [Google Scholar] [CrossRef] [Green Version]
  73. Domagalska, M.A.; Imamura, H.; Sanders, M.; Van den Broeck, F.; Bhattarai, N.R.; Vanaerschot, M.; Maes, I.; D’Haenens, E.; Rai, K.; Rijal, S.; et al. Genomes of intracellular Leishmania parasites directly sequenced from patients. bioRxiv 2019. [Google Scholar] [CrossRef] [Green Version]
  74. Cao, X.; Aufsatz, W.; Zilberman, D.; Mette, M.F.; Huang, M.S.; Matzke, M.; Jacobsen, S.E. Role of the DRM and CMT3 methyltransferases in RNA-Directed DNA methylation. Curr. Biol. 2003, 13, 2212–2217. [Google Scholar] [CrossRef] [PubMed]
  75. De Mendoza, A.; Bonnet, A.; Vargas-Landin, D.B.; Ji, N.; Li, H.; Yang, F.; Li, L.; Hori, K.; Pflueger, J.; Buckberry, S.; et al. Recurrent acquisition of cytosine methyltransferases into eukaryotic retrotransposons. Nat. Commun. 2018, 9, 1341. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Dubin, M.J.; Zhang, P.; Meng, D.; Remigereau, M.S.; Osborne, E.J.; Paolo Casale, F.; Drewe, P.; Kahles, A.; Jean, G.; Vilhjalmsson, B.; et al. DNA methylation in Arabidopsis has a genetic basis and shows evidence of local adaptation. eLife 2015, 4, e05255. [Google Scholar] [CrossRef] [PubMed]
  77. Cokus, S.J.; Feng, S.; Zhang, X.; Chen, Z.; Merriman, B.; Haudenschild, C.D.; Pradhan, S.; Nelson, S.F.; Pellegrini, M.; Jacobsen, S.E. Shotgun bisulphite sequencing of the Arabidopsis genome reveals DNA methylation patterning. Nature 2008, 452, 215–219. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Hsieh, T.F.; Ibarra, C.A.; Silva, P.; Zemach, A.; Eshed-Williams, L.; Fischer, R.L.; Zilberman, D. Genome-wide demethylation of Arabidopsis endosperm. Science 2009, 324, 1451–1454. [Google Scholar] [CrossRef] [Green Version]
  79. Zemach, A.; McDaniel, I.E.; Silva, P.; Zilberman, D. Genome-wide evolutionary analysis of eukaryotic DNA methylation. Science 2010, 328, 916–919. [Google Scholar] [CrossRef] [Green Version]
  80. Verma, S.; Kumar, R.; Katara, G.K.; Singh, L.C.; Negi, N.S.; Ramesh, V.; Salotra, P. Quantification of parasite load in clinical samples of leishmaniasis patients: IL-10 level correlates with parasite load in visceral leishmaniasis. PLoS ONE 2010, 5, e10107. [Google Scholar] [CrossRef]
  81. Cuypers, B.; Berg, M.; Imamura, H.; Dumetz, F.; De Muylder, G.; Domagalska, M.A.; Rijal, S.; Bhattarai, N.R.; Maes, I.; Sanders, M. Integrated genomic and metabolomic profiling of ISC1, an emerging Leishmania donovani population in the Indian subcontinent. Infect. Genet. Evol. 2018, 62, 170–178. [Google Scholar] [CrossRef]
  82. Imamura, H.; Downing, T.; Van den Broeck, F.; Sanders, M.J.; Rijal, S.; Sundar, S.; Mannaert, A.; Vanaerschot, M.; Berg, M.; De Muylder, G.; et al. Evolutionary genomics of epidemic visceral leishmaniasis in the Indian subcontinent. eLife 2016, 5, e12613. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Olova, N.; Krueger, F.; Andrews, S.; Oxley, D.; Berrens, R.V.; Branco, M.R.; Reik, W. Comparison of whole-genome bisulfite sequencing library preparation strategies identifies sources of biases affecting DNA methylation data. Genome Biol. 2018, 19, 33. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Bewick, A.J.; Vogel, K.J.; Moore, A.J.; Schmitz, R.J. The evolution of DNA methylation and its relationship to sociality in insects. bioRxiv 2016, 062455. [Google Scholar] [CrossRef]
  85. Oyola, S.O.; Gu, Y.; Manske, M.; Otto, T.D.; O’Brien, J.; Alcock, D.; Macinnis, B.; Berriman, M.; Newbold, C.I.; Kwiatkowski, D.P.; et al. Efficient depletion of host DNA contamination in malaria clinical sequencing. J. Clin. Microbiol. 2013, 51, 745–751. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Protein alignment of LdDNMT (LdBPK_251230) and TbDNMT generated with T-coffee, picturing the similarities between the 10 homologous domains of C5 DNA methyltransferases. Black highlights homology, and the red character displays the position of the catalytic cysteine residue.
Figure 1. Protein alignment of LdDNMT (LdBPK_251230) and TbDNMT generated with T-coffee, picturing the similarities between the 10 homologous domains of C5 DNA methyltransferases. Black highlights homology, and the red character displays the position of the catalytic cysteine residue.
Microorganisms 08 01252 g001
Figure 2. RAxML maximum likelihood tree showing the position of trypanosomatid DNMT (DNMT 6) within the DNMT family. Displayed branch bootstrap values are based on 1000 bootstraps. Line thickness is scaled for these bootstrap values.
Figure 2. RAxML maximum likelihood tree showing the position of trypanosomatid DNMT (DNMT 6) within the DNMT family. Displayed branch bootstrap values are based on 1000 bootstraps. Line thickness is scaled for these bootstrap values.
Microorganisms 08 01252 g002
Figure 3. CpG, CHG, and CHH genome-wide methylation patterns in (A) Leishmania donovani BPK282 P3 promastigotes (36 chromosomes), (B) Trypanosoma brucei brucei TREU927 (11 chromosomes), and (C) Arabidopsis thaliana Col-0 (5 chromosomes). Data was binned over 10,000 positions to remove local noise and variation.
Figure 3. CpG, CHG, and CHH genome-wide methylation patterns in (A) Leishmania donovani BPK282 P3 promastigotes (36 chromosomes), (B) Trypanosoma brucei brucei TREU927 (11 chromosomes), and (C) Arabidopsis thaliana Col-0 (5 chromosomes). Data was binned over 10,000 positions to remove local noise and variation.
Microorganisms 08 01252 g003
Figure 4. DNA/gene copy number based on genomic sequencing depth on chromosome 25 position 465,000–475,000. Both the LdDNMT knock-out (LdDNMT-/-) and LdDNMT overexpressor lines (LdDNMT+) were successful with, respectively, 0 and 64 copies of the gene. The plot shows also that the neighboring genes LdBPK_251220 and LdBPK_251240 were unaffected and had the standard disomic pattern.
Figure 4. DNA/gene copy number based on genomic sequencing depth on chromosome 25 position 465,000–475,000. Both the LdDNMT knock-out (LdDNMT-/-) and LdDNMT overexpressor lines (LdDNMT+) were successful with, respectively, 0 and 64 copies of the gene. The plot shows also that the neighboring genes LdBPK_251220 and LdBPK_251240 were unaffected and had the standard disomic pattern.
Microorganisms 08 01252 g004
Figure 5. Enrichment (X) of Leishmania DNA in artificial mixtures of Leishmania promastigote DNA and human DNA, with the mixtures ranging from 1:15 to 1:15,000 Leishmania:human DNA. Enrichments were carried out with the NEBNext Microbiome DNA Enrichment Kit (NEB), and the unmethylated Leishmania DNA was enriched on average 263 times.
Figure 5. Enrichment (X) of Leishmania DNA in artificial mixtures of Leishmania promastigote DNA and human DNA, with the mixtures ranging from 1:15 to 1:15,000 Leishmania:human DNA. Enrichments were carried out with the NEBNext Microbiome DNA Enrichment Kit (NEB), and the unmethylated Leishmania DNA was enriched on average 263 times.
Microorganisms 08 01252 g005
Table 1. CpG, CHG, and CHH methylation percentages in different Leishmania donovani lines (Ld), Trypanosoma brucei, and Arabidopsis thaliana (positive control).
Table 1. CpG, CHG, and CHH methylation percentages in different Leishmania donovani lines (Ld), Trypanosoma brucei, and Arabidopsis thaliana (positive control).
CpG (%)CHG (%)CHH (%)
LdPro0.00030.00050.0126
LdAmas0.00010.00030.0073
LdHamPro0.00020.00050.0113
LdDNMT+0.00130.00260.0627
LdDNMT-/-0.00020.00060.0079
Tbrucei0.00010.00060.0040
Athaliana21.04734.04010.3141
Table 2. qPCR estimation of LdBPK_251230 expression level (copy number) of Ldo-Pro and Ldo-DNMToverex.
Table 2. qPCR estimation of LdBPK_251230 expression level (copy number) of Ldo-Pro and Ldo-DNMToverex.
Ldo-ProLdDNMT+
RNA1.53 ± 0.23.78 ± 0.3
Table 3. Enrichment (X) of Leishmania DNA in artificial mixtures of Leishmania and human DNA for promastigotes and amastigotes of 3 clinical isolates (BPK026, BPK275, and BPK282). Enrichments were carried out with the NEBNext Microbiome DNA Enrichment Kit (NEB).
Table 3. Enrichment (X) of Leishmania DNA in artificial mixtures of Leishmania and human DNA for promastigotes and amastigotes of 3 clinical isolates (BPK026, BPK275, and BPK282). Enrichments were carried out with the NEBNext Microbiome DNA Enrichment Kit (NEB).
BPK026BPK275BPK282Average Enrichment (X)St.Dev
Promastigotes79.8588.3260.4776.2214.28
Amastigotes64.8356.8763.3361.684.23

Share and Cite

MDPI and ACS Style

Cuypers, B.; Dumetz, F.; Meysman, P.; Laukens, K.; De Muylder, G.; Dujardin, J.-C.; Domagalska, M.A. The Absence of C-5 DNA Methylation in Leishmania donovani Allows DNA Enrichment from Complex Samples. Microorganisms 2020, 8, 1252. https://doi.org/10.3390/microorganisms8081252

AMA Style

Cuypers B, Dumetz F, Meysman P, Laukens K, De Muylder G, Dujardin J-C, Domagalska MA. The Absence of C-5 DNA Methylation in Leishmania donovani Allows DNA Enrichment from Complex Samples. Microorganisms. 2020; 8(8):1252. https://doi.org/10.3390/microorganisms8081252

Chicago/Turabian Style

Cuypers, Bart, Franck Dumetz, Pieter Meysman, Kris Laukens, Géraldine De Muylder, Jean-Claude Dujardin, and Malgorzata Anna Domagalska. 2020. "The Absence of C-5 DNA Methylation in Leishmania donovani Allows DNA Enrichment from Complex Samples" Microorganisms 8, no. 8: 1252. https://doi.org/10.3390/microorganisms8081252

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop