Roma Sinha1†‡
Mathu Malar C2,3‡
Raghwan1†‡Subhadeep Das2,3
Sonali Das1
Mohammad Shadab1†Rukhsana Chowdhury1
Sucheta Tripathy2,3*
Nahid Ali1*- 1Infectious Diseases and Immunology Division, CSIR-Indian Institute of Chemical Biology, Kolkata, India
- 2Structural Biology and Bioinformatics Division, CSIR-Indian Institute of Chemical Biology, Kolkata, India
- 3Academy of Scientific and Innovative Research (AcSIR), New Delhi, India
Leishmania donovani possesses a complex heteroxenic life cycle where infective metacyclic promastigotes are pre-adapted to infect their host and cope up with intracellular stress. Exploiting the similarities between cultured and sandfly derived promastigotes, we used early and late passage cultured promastigotes to show specific changes at genome level which compromise pathogen fitness reflected in gene expression and infection studies. The pathogen loses virulence mostly via transcriptional and translational regulations and long-time cultivation makes them struggle to convert to virulent metacyclics. At the genomic level very subtle plasticity was observed between the early and the late passages mostly in defense-related, nutrient acquisition and signal transduction genes. Chromosome Copy number variation is seen in the early and late passages involving several genes that may be playing a role in pathogenicity. Our study highlights the importance of ABC transporters and calpain like cysteine proteases in parasite virulence in cultured promastigotes. Interestingly, these proteins are emerging as important patho-adaptive factors in clinical isolates of Leishmania. We found that the currently available genome of Leishmania in the NCBI database are from late passages. Our early passage genome can act as a reference for future studies on virulent isolates of Leishmania. The annotated leads from this study can be used for virulence surveillance and therapeutic studies in the Indian subcontinent.
Introduction
Leishmania donovani is an intracellular obligate parasite of mammalian macrophages and causes visceral leishmaniasis or kala-azar, which is a fatal disease if not treated on time. The disease severity varies from host to host, region to region as well as between parasite strains. Extensive research by different groups exploiting genomic, proteomic, metabolomic, immunologic, and animal models have pointed toward tripartite determinants mediating the development of this disease viz. the vector, host and pathogen. Although vector and host characteristics are important for symptomatic disease, parasite characteristics majorly determine the fate of the disease (McCall et al., 2013). Promastigotes are the culturable form of Leishmania. Stationary phase culture of Leishmania is expected to contain a large number of metacyclic parasites and have been routinely used for experimental infections. Metacyclic promastigotes are injected into the host during blood meal and are the infective stage of the parasite. We and others have observed that continuous axenic cultivation of Leishmania leads to loss of virulence over a period of time (Dey et al., 2002; Ali et al., 2013). The relative virulence of stationary phase promastigotes is proportional to peanut agglutinin negative promastigotes contained within these populations (da Silva and Sacks, 1987) which is attributable to exposed surface carbohydrates. It has also been observed that the in vitro maintenance of Leishmania promastigotes by cultivation over longer periods may reduce their ability to differentiate into amastigote forms (Moreira et al., 2012; Ali et al., 2013). Studies at the mRNA and proteomic levels indicated differential expression of virulence related genes in both promastigotes and amastigotes cultured axenically for long periods compared to freshly isolated/transformed or early passage parasites (Lei et al., 2010; Pescher et al., 2011; Ali et al., 2013; Magalhaes et al., 2014). These studies were done with the aim of standardizing protocols for the use of axenic Leishmania cultures in infection studies and to identify virulence factors. Studies at genomic level mostly identified parasite evolution in the Indian subcontinent under drug pressure or disease phenotype (Zhang et al., 2014; Imamura et al., 2016). As promastigotes are the infective form of Leishmania, and cultured stationary phase promastigotes show remarkable consistency with sandfly derived metacyclic promastigotes (Inbar et al., 2017), we chose the former system to unravel the genetic mechanism of parasite pre-adaptation to host environment and how it is lost with continuous in vitro passages. We report here the global changes in the genome and transcriptome of serially passaged L. donovani AG83 strain which cause specific alterations in few genes and their expression that lead to loss of virulence reflected in the infection studies. We also try to link these adaptive features with parasite’s nutritional environment and complex life cycle, and how closely they relate to the global evolution of L. donovani clinical isolates as available through published work. This study also throws light on the pathologically significant genes common to clinical and laboratory strains which must be considered for virulence surveillance at least in the Indian subcontinent.
Materials and Methods
Ethics Statement
All animal experiment protocols adhered to the guidelines of the Committee for the Purpose of Control and Supervision on Experimental Animals (CPCSEA), Ministry of Environment and Forest, Government of India, and were approved by the Animal Ethics Committee (147/1999/CPSCEA) of CSIR-IICB.
Parasite Culture and Maintenance in Animals
Leishmania donovani (MHOM/IN/1983/AG83; ATCC repository number PRA®-413TM) amastigotes from infected hamster spleen were transformed to promastigotes in Schneider’s Drosophila medium (Sigma-Aldrich) at 22°C and sub-cultured as promastigotes in M 199 (Sigma-Aldrich) both supplemented with 100 U/ml penicillin, 100 μg/ml streptomycin, and 10% FCS (Sigma-Aldrich) (Banerjee et al., 2008). Parasites were sub-cultured till 25th passage for some experiments. Golden Syrian Hamsters, 5–6 weeks old, reared in institute facilities were used for the purpose of parasite maintenance.
Macrophage Infection in Vitro
Macrophages isolated from the peritoneal cavity of hamsters (12–16 weeks old) by injecting chilled RPMI-1640-10% FCS were pooled and cultured overnight on glass cover slips as described elsewhere (Bhowmick et al., 2010). AG83 promastigotes of different passages were allowed to infect peritoneal macrophages (MoI of 10:1) for 3 h, washed with warm 20 mM PBS, pH 7.2 and further incubated in fresh medium till 72 h. At 24 and 72 h post-infection, cells were fixed in methanol and stained with Giemsa for microscopic determination of intracellular parasite numbers per 100 host cells. Data were analyzed using GraphPad Prism 5 software.
Determination of Splenic and Hepatic Parasite Burden
Hamsters (5–6 weeks) were infected by injecting 2 × 107 stationary phase promastigotes by intracardiac injection suspended in 20 mM PBS. Animals were sacrificed 8 weeks post-infection for determination of parasite burden by Leishman Donovan Units (LDU) in Giemsa stained impression smears and limiting dilution assay in fivefold serial dilutions of homogenized organs as described previously (Banerjee et al., 2008). Data analysis was performed on GraphPad Prism 5 software.
High Throughput Genome Sequencing
High quality genomic DNA was extracted from stationary phase promastigotes of the early (HTI4) and late (HTI5) passage of L. donovani AG83 using a commercial procedure as recommended by the manufacturer (Qiagen, Germany). Size check, integrity and presence of contaminants in the DNA samples were assessed through gel electrophoresis. DNA purity was measured using a Nano Drop 2000 spectrophotometer (Thermo Scientific, Waltham, MA, United States). Two separate libraries of L. donovani AG83 were prepared one each for early and late passages. De novo paired end sequencing was done on Illumina HiSeq 2500 with the read length of 125 base pairs with an insert size of 250 bp. A total of 66 and 44 million raw reads were generated from early and late passages, respectively.
Transcriptome Sequencing
For each passage, parasites were harvested by chilling on ice, spun down, and washed once in cold PBS solution, pH 7.4, and suspended in RNA Later. Total RNA was isolated from promastigotes of early (2nd), intermediate (11th), and late (25th) passage (referred to as q1, q2, and q3, respectively, in the figures) by using “Roche High Pure RNA Isolation Kit,” (Product no.11828665001) according to manufacturer’s protocol. The purity and concentration of each RNA sample was checked by using the Agilent 2100 bioanalyzer (Agilent Technologies, CA, United States) before proceeding for further downstream analyses. Paired end transcriptome sequencing was carried out on Illumina HiSeq platform, with the read length of 125 base pairs with an insert size of 250 bp generating 37–42 million raw reads (approximately). The library has been sequenced following manufacturer’s instructions using the HiSeq SBS Kit v4 (Part # 15034097Rev.B) to generate paired-end reads. Additional quality control of raw data using FastQC (11) was performed. The reads were preprocessed using Trimmomatic and the poor quality sequences, contaminated sequences were removed using the blast and blat based searches using the reference (LdBPK82A1; Bioproject: PRJNA171503) from Genbank. The high quality filtered reads were used for downstream analyses.
Genome and Transcriptome Assembly and Functional Annotation
The early and late passage genomes were assembled using Allpaths-LG assembler (Ribeiro et al., 2012). We simulated reads with 6 K and 20 K insert sizes from the existing genomes of LdBPK282A1 retrieved from Genbank (PRJNA171503) with WGsim (Li et al., 2008).
The assembly from Allpaths was resolved up to chromosome level reference based genome assembly by nucleotide alignment using Nucmer and further using our in-house scripts1. RNAseq reads from all passages were aligned using blat (Kent, 2002) with the Genbank data of LdBPK282A1 for removal of contaminants. The cleaned reads were used for transcriptome assembly using Trinity transcriptome assembler (Grabherr et al., 2011; Haas et al., 2013) (Supplementary Table S1). We used BUSCO (Simão et al., 2015) for predicting the genome completeness and the core ortholog analysis was computed using eukaryotic genomes (Supplementary Table S2).
We predicted protein coding genes of the assembled genomes using Augustus (Stanke et al., 2008) and Scipio pipeline with LdBPK282A1 coding sequences as training material (Bioproject: PRJNA171503). Predicted genes were annotated using BLAST (Altschul et al., 2009) as well as InterProScan (Quevillon et al., 2005) searches. Whole genome protein sequences were submitted to the gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. GO ids were mapped back to their ontology functions using in house scripts. GO enrichment was analyzed in the categories of molecular function, cellular component and biological process. Metabolic pathways analysis was done using the KASS server (Moriya et al., 2007).
Additionally, we used companion server (Steinbiss et al., 2016) for predicting orthologous genes from the two culture passages using Genbank Leishmania major as the reference. This server predicted many interesting features including the pseudogenes present in the genomes that were not predicted by Augustus. Genomes of two passages were compared using MUMmer, NUCmer, and PROmer (Kurtz et al., 2004).
Differential Expression Analysis of Transcripts
We carried out differential expression analysis of transcripts using Tuxedo suit tools (Langmead et al., 2009; Trapnell et al., 2012; Kim et al., 2013) by aligning the RNAseq reads of three passages with the reference genome of LdBPK282A1. The pattern of expression across samples was calculated by using the normalized Fragments per Kilo base of transcript per Million (FPKM) values. R scripts were used to construct box plots and perform hierarchical clustering (Supplementary Figure S1). Putative secretary functions were assigned to the DEGs by using Secretome tool (Bendtsen et al., 2004).
Total RNA was isolated from 4 × 107 parasite promastigotes of different passages, using TriZOL-chloroform method. Pure RNA pellet was suspended in 50 μl DEPC treated RNase free water and concentration was measured using Nanodrop 2000 (Thermo Fisher Scientific). Almost 2 μg pure RNA was used for cDNA synthesis following manufacturer’s protocol using Revert Aid First Strand cDNA synthesis kit (Thermo Fisher, MA, United States). Differential expression analysis of parasite genes was performed using SYBR green master mix (Roche). Ld GAPDH was used as reference gene. Reaction was run in Roche light cycler 96 (Roche) and fold change of different genes was expressed as 2-ΔΔCt, where ΔCt = gene Ct- reference gene Ct and Δ(ΔCt) = ΔCt test – ΔCt control. RNA was collected from n = 3 hamster and proceeded for further long passage preparation. Also data per group (A2, A11, and A25th passage) represents n = 3 experiments performed for each group.
SNP Identification
Sequenced illumina genomic reads of early and late passages were aligned with the reference genome of LdBPK282A1 using Burrows-Wheeler Aligner (Li and Durbin, 2009), After marking the duplicates using PICARD tools the variants were called using the Genome Analysis toolkit (GATK) (Van der Auwera et al., 2013), obtained variants were filtered with the cutoff of QD < 2.0 ||MQ < 40||FS > 60.0|| ReadPosRankSum < -8.0.
Copy Number Variant Detection
The reads mapped by using BWA were used to estimate the read depth for each sample at each chromosome using the coverage values from Samtools v0.1.18 (Li et al., 2009). After removing the duplicates using Picard the CNVnator (Abyzov et al., 2011) was used to identify copy number variants. We used 100 bp as the window size for creating the histogram bins of read depth in CNVnator. For CNVnator, we removed calls which were less than the cutoff of t-test calculated e-value of 1e-5.
Cell Viability and Apoptotic Cell Death Assay
Promastigotes of early (2nd) and late (25th) passages were treated with graduated doses of miltefosine for 24 h and the percent viability was determined by enzymatic reduction of 3-[4,5-dimethylthiazole-2-yl]-2,5-diphenyltetrazolium bromide (MTT) to MTT-formazan. In parallel experiments the miltefosine treated promastigotes were analyzed for apoptotic cell death by AnnexinV-FITC and propidium iodide staining acquired in FITC and PE channels, respectively, on FACS LSR Fortessa and analyzed for percent FITC positive and FITC-PI double positive cells on FACS Diva software.
Preparation of Leishmanial Antigen
The promastigotes (2nd–5th and 22nd–25th passage) were washed three times in cold 20 mM PBS, pH 7.2 and re-suspended at a concentration of 1.0 g cell pellet in 50 ml of cold 5 mMTris–HCl buffer, pH 7.6 and Leishmanial antigen (Lag) was prepared by ultrasonication as described earlier (Bhowmick et al., 2010). LAg was precipitated in 90% acetone-10% TCA, solubilised in rehydration buffer (7 M urea, 2 M thiourea, 2% Nonidet P-40, 2% dithiothreitol, and 2% pharmalyte, all from GE Healthcare LifeSciences, Little Chalfont, United Kingdom except NP-40 from Sigma) overnight and estimated by Bradford method (Bradford, 1976). Two independent batches were prepared from each passage for analysis.
Isoelectric Focusing
The isoelectric focusing (IEF) was performed using the PROTEAN IEF Cell system (Bio-Rad). 500–600 μg of whole cell protein was added to 300 μl of rehydration buffer for 11 cm IPG strip (pH 4–7 non-linear, Bio-Rad). Protein sample was applied to the IPG strips and kept at room temperature (16–18 h) for passive rehydration. IEF was performed at 250 V for 20 min; 10,000 V for 2 h 50 min; 10,000–40,000 V/h and an optional step of 500–1000 V 10–15 h.
Second Dimension-SDS-PAGE
Before the second dimension, proteins in the IEF strips were reduced in equilibration buffer (6 M urea, 2% SDS, 300 mM Tris–HCl pH 8.8, 20% glycerol) containing 20 mg/ml DTT and alkylated in dark with 25 mg/ml iodoacetamide, again dissolved in equilibration buffer, for 15–20 min each. Strips were washed with SDS-PAGE running buffer and separated across 12% SDS-PAGE gels (30% acrylamide, 0.8% bis-acrylamide) using a vertical system (Bio-Rad) and standard Tris/glycine/SDS buffer. Gels were run at 30 mA/gel until the tracking dye left the gel. The protein standard was purchased from Bio-Rad (pre-stained SDS-PAGE broad range). For western blot analysis, total cell lysates of early and late passage promastigotes were run on SDS-PAGE followed by incubation with antibodies for gp63 and HDAC. The images were developed using Luminata Forte (Thermo Fisher, MA, United States) on Gel Doc XR (Bio-Rad). The anti-gp63 antibody raised in rabbit (in house) and anti-mouse HDAC antibody (Cell Signaling Technology) were used for the assay.
Spot Handling and Tryptic Digestion for Matrix-Assisted Laser Desorption/Ionization-Time of Flight Mass Spectrometry (MALDI-ToF MS/MS)
The 2-DE gels were stained with Coomassie Brilliant Blue G-250. Gel images were taken on the Gel Doc XR+ (Bio-Rad) and analysis was done using the PD Quest software version 8.0.1 (Bio-Rad Laboratories) and manual checking. The gel having a higher number of spots was used to locate the corresponding protein spots within different gels. In this study, we have applied a cut-off of at least 1.5-fold up-regulation/down-regulation of protein for the study. Experimental molecular weights and isoelectric points were calibrated according to Bio-Rad standard 2-DE PAGE marker. Spots were manually excised and processed according to the In-Gel Tryptic Digestion Kit (Thermo Fisher, MA, United States). The trypsin digested samples were concentrated using speed-vac before loading the samples for MALDI-TOF Mass Spectrometry analysis on 4800 MALDI ToF/ToF analyzer (Applied Biosystems, CA, United States). Zip Tip protein concentrator (Millipore, MA, United States) was used for concentrating and de-salting of the low abundance proteins. The samples were dissolved in a solvent consisting of 0.1% trifluoroacetate and 50% acetonitrile in MilliQ Water. Then, 0.5 μl of sample solution was mixed with 0.5 μl of matrix solution (1 mg/ml α-cyano-4-hydroxycinnamic acid dissolved in the aforementioned solvent), applied to a 384-MALDI sample target plate, and dried in air. Peptides were evaporated with a ND: YAG laser at 355 nm, using a delayed extraction approach. They were accelerated with 25 kV injection pulse for ToF analysis. Each spectrum was the cumulative average of 1000 laser shots. The MS/MS spectrum was collected in MS/MS 1 kV positive reflectron mode with fragments generated by post source decay. The MS/MS mass tolerance was set to ±20 ppm. After processing, 10 MS/MS precursors were selected (Minimum signal to noise ratio-50). Before each analysis, the instrument was calibrated with the Applied Biosystems 4700 Proteomics Analyzer Calibration Mixture. Data interpretation was performed using the GPS Explorer Software (Applied Biosystems, CA, United States), and an automated database search was performed using the MASCOT program (Matrix Science Ltd., London, United Kingdom).
Data Availability
The Genome sequences of Early and Late passages of L. donovani MHOM/IN/1983/AG83 is available at NCBI with submission numbers GCA_001989955.1 and GCA_001989975.1, respectively.
Results and Discussion
Continuous in Vitro Passaging Leads to Loss of Virulence of Promastigotes
We first compared the infection causing capacity of promastigotes from different passages in vitro and in vivo. Stationary promastigotes recovered from different passages were quantified and employed in the experiments. Hamsters and hamster derived macrophages were highly susceptible to infection, with infection persisting up to 25th passage (Figures 1, 2) and displayed a decline in infectivity with more time spent in the culture. The average number of amastigotes per macrophage was 14.4 ± 3.2 at 24 h, which increased to 21.6 ± 3.9 at 72 h (Figure 1A) post-infection with early passage. There was over a 50-fold reduction in the number of intra-macrophagic amastigotes compared to the 2nd passage when parasites of 25th passage were used for infection (Figure 1, p < 0.0001, One way Anova). Parasites from the 2nd passage could infect 96.83 ± 1.4% of macrophages at 72 h while only 8 ± 1.6% cells were infected with promastigotes from 25th passage (Figure 1B). From these experiments it seemed the infectivity is affected at two stages, first at the level of host-parasite interaction evident from lower infection at 24 h, and second, parasite sustenance inside macrophages, as seen in clearance of parasites after 72 h. Hamsters infected with 25th passage promastigotes presented over 6.5-fold less parasite burden in both liver (LDU of 515.6 ± 82.44) and spleen (LDU of 211.5 ± 42.21) as compared to those infected with the 2nd passage (LDU of 3358 ± 132.9 and 1416 ± 126.6, respectively, in liver and spleen) (Figures 2A,B, p < 0.0001, One Way Anova) at 8 week post-infection. As the number of live amastigotes inside the organs gives an idea of the disease severity, we additionally assessed the parasite burden in the infected liver and spleen of differentially infected hamsters by serial dilution method. The number of parasites showed a gradual decline with passage number thus reinstating the fact that continuous axenic culture of promastigotes leads to loss of virulence of the parasite. As shown in Figure 2C, the number of live amastigotes came down from 1016 in the liver of hamsters infected with virulent parasites (2nd passage) to as low as 100 parasites in the animals with attenuated infection (25th passage) (p < 0.0001, One Way Anova). The splenic infection level was very high in hamsters infected with early passage promastigotes, and reached levels as high as 1021 (Figure 2D, p < 0.0001, One Way Anova) whereas the parasites failed to survive and multiply in the animals infected with late passage, and only about 100 parasites could be detected in the piece of organ examined by us. The hamsters with mild infection also demonstrated reduced liver and spleen size (Figure 2E) unlike the severely infected groups. Thus early passage promastigotes displayed increased fitness inside host cells reflected in higher parasite burden while late passage promastigotes were successfully cleared by macrophages.
FIGURE 1. In vitro hamster derived macrophage infection with different passages. Hamster macrophages were infected with stationary promastigotes of Leishmania donovani AG83, as described in Section “Materials and Methods” and the cells were incubated for 24 h and 72 h at 37°C, 5% CO2. (A) The number of amastigotes in each passage and (B) the percentage of infected cells were analyzed by counting 100 cells. Data shown are mean ± SE of experiment performed in quadruplet.
FIGURE 2. Infection in hamsters. Hamsters were infected by intracardiac injection of 2 × 107 L. donovani AG83 promastigotes. Animals were sacrificed at 8 week post-infection and (A,C) liver and (B,D) spleen parasite burden was determined by (A,B) LDU and (C,D) LDA. Data represent mean ± SE of two independent experiments (n = 4–6). (E) The mean weights of liver and spleen were compared with respective organ weights of uninfected animals.
Late Passage Genome Has Higher Number of Pseudogenes
Genome plasticity in Leishmania has been linked to its ability to adapt to different environments (Sterkers et al., 2012). This prompted us to check for the genetic adaptations which may account for adaptability to growth in culture medium and the associated loss of infectivity as a result of continuous passaging. The genomes of early and late passages were assembled into complete 36 chromosomes with a genome size of 32.2 and 32.1 MB, respectively (Table 1A). The genomes of these two passages, however, had about 11.3% gaps compared to 3.67% gaps in LdBPK282A1 from Genbank [Accession number: GCA_000227135.2]. A chromosome wise comparison of early and late passages with LdBPK282A1 (Supplementary Figure S2) didn’t reveal much difference, whereas the whole genome comparisons clearly pointed out the small but prominent changes in the genome of the late passage, which was comparable to the NCBI reference LdBPK282A1 strain (Figure 3, upper panel). We used RATT and companion for transferring annotation from LdBPK282A1 strain (Downing et al., 2011) that resulted in 7563 and 7552 predicted protein coding genes, respectively, in early and late passages, compared to 7967 protein coding genes in LdBPK282A1 (Table 1B). The less number of protein coding genes from our cultures are mostly attributed to the presence of assembly gaps. The numbers of predicted pseudogenes in early and late passages were 73 and 82, respectively, and this number is significantly higher than 54 pseudogenes predicted in LdBPK282A1 strain, but similar to what is reported in the Srilankan L. donovani (70) strains (Zhang et al., 2014). Presence of large protein families such as calpain-like proteases and amastins are found to occur in large directional clusters in Genbank strain as well as our strains. We noticed a large amastin cluster in chromosome 8 in Genbank that are present in chromosome 10 in our strain (Table 2).
TABLE 1A. Summary of assembly statistics for early passage (HTI4), late passage (HTI5) and Genbank Leishmania donovani LdBPK282A1 genomes.
TABLE 1B. Gene prediction statistics for early (HTI4) and late (HTI5) passages using companion server.
FIGURE 3. Mummerplot comparisons and Gene Ontology analysis of early and late passage genomes to the NCBI reference genomes. (Upper) The y-axis represents HTI4 (left) or HTI5 (middle) whole genome assemblies compared to x-axis representingLdBPK282A1 assembly. The HTI4 (x-axis) versus HTI5 (y-axis) chromosomal assembly comparison is given in the right panel. The dots represent the positions of conserved DNA sequences on the genomes. (Lower) REVIGO was used to visualize the summary of significantly enriched GO terms. The scatterplots show the cluster representatives in a two dimensional space derived by applying multidimensional scaling to a matrix of the GO terms’ semantic similarities for early (left panel) and late (right panel) passages. Clustal differences are represented in boxes. Bubble color indicates the user-provided p-value; size indicates the frequency of the GO term in the underlying GOA database (bubbles of more general terms are larger).
TABLE 2. Gene Copy number in early, late passages and Genbank strain LdBPK282A1.
The GO Biological Process annotation of early and late passages indicates there are no changes in the gene numbers of processes in early and late passages. However, genes responsible for defense system (GO: 0006952) and positive regulation of cell proliferation (GO: 0008284) are missing in later passages (Figure 3, lower panel and Supplementary Data Sheet S1). Thus repeated ex host passages have made the late passage parasites suited to extracellular life where selection pressure against host defense system is not needed. Further, the cellular machinery needed for transformation to infective form is also perturbed which may impair their ability to infect the host.
Chromosome Copy Number Variation Drives Gene Expression in Leishmania
The previous studies on L. donovani (Leprohon et al., 2009; Dumetz et al., 2017) reported that aneuploidy and Copy number variants regulate the gene expression. In our study we identified SNPs from early and late passages. Total 4390 heterozygous loci were found in early passage and 4356 heterozygous loci were found in late passage of the genome.
The local copy number variant detection was performed on the basis of read depth binning ratio approach. It was already known that deletion and duplication event plays an important role in these genomes to maintain the virulence in the different environmental conditions (Dumetz et al., 2017).
The copy Number variation on our passages was done using LdBPK282A1 as reference, with an e-value cut-off of 1e-5 for filtering. Three hundred and eighty CNV events were qualified in late passage and 365 CNV events were qualified in early passage (Supplementary Data Sheet S2). The size of CNV events ranged from 0.7 to 271 kb. We checked for the regions which were overlapping with the protein coding genes. A total of 230 deletion events and 135 duplication events were identified in the early passage while 234 deletions events and 146 duplication events were identified in the late passage genome (Supplementary Data Sheet S2). Uniq gene lists that had undergone changes are listed in Supplementary Data Sheet S2: uniq_CNV. The genes undergoing structural changes in later passages compared to earlier passages comprised of several ABC Transporters (10), Amino acid transporters (7), Amastins (3), GP63, calpain like cysteine proteases to name a few. Interestingly, we have also reported changes in gene expression profiling of GP63 and Calpain like cysteine proteases between the late and early passages. We found more CNV events in chromosomes 5, 6, 8, 15, and 31 in the genome of early passage as per earlier reports (Laffitte et al., 2016; Iantorno et al., 2017). High number of CNV events play a major part in the gene expression regulation of in vitro cultured promastigotes (Dumetz et al., 2017). Interestingly we observed an increased frequency of CNVs in the protein coding regions of the late passage particularly in chromosome 3, 4, 13, 16, and 20. Protein coding regions corresponding to calpain like cysteine protease genes on chromosome 20 displayed multiple duplication events. Phosphoglycerate kinase B, cytosolic fragment on chromosome 20 showed an insertion event in the late passage which was not found in the earlier passage. This gene was reported to be over expressed in Srilankan L. donovani strains causing visceral disease versus those causing cutaneous lesions (Zhang et al., 2014). Noticeable protein coding genes falling in the variant regions on chromosome 13 of late passage genome were some acetyl transferases including histone acetyltrasferase, N-acetyl transferase subunit ARD1, RAS-related protein RAB5, mitogen activated protein kinase 2, etc. which may have a role in cell cycle progression and morphogenesis of Leishmania (Wiese, 1998; Yadav et al., 2016).
Single Nucleotide Polymorphisms in the ABC Transporter Coding Genes
The evolution of pathogenicity in microbes has been attributed to ordered changes in the functionality of the genes as a result of physiological constraints encountered by the organism in their immediate environment. A KEGG pathway analysis didn’t show major differences in the number of genes involved (Supplementary Data Sheet S3) although detailed mutational analysis revealed interesting changes linked directly or indirectly to loss of virulence in the later passage. Among the defense gene products, ABC transporters are important proteins involved in drug resistance, nutrient acquisition and pathogen virulence (Glavinas et al., 2004) and have expanded in the pathogenic protists as a parasitic adaptation (Zhang et al., 2014; Jackson et al., 2016). Genomes of L. donovani AG83 in early and late passages contain 42 copies of the genes spread over the genome. However, there are at least two instances where the ABC transporters are undergoing polymorphism that could result in reduced pathogen fitness. ABC transporter at LDON_230007700.1; HTI5:chr23: 91219–96174 in late passage has undergone several changes at the nucleotide level leading to an inactive ABC transporter in the late passages. Our Copy Number variant analysis also detected CNV duplication event of 12000 bp in the region of 86801–98800 in chromosome 23 in late passage. In earlier passages, a functional ABC transporter [HTI4:Chr23:91908–92025] resided inside a larger gene locus at HTI4:chr23:89807–101617 (Figure 4). We found a CNV duplication event as well in chromosome 23: 86801–98700 in early passage that corroborates this finding.
FIGURE 4. Nucleotide and protein sequence comparison in L. donovani ABC transporter gene in Chr 23 of early and late passages. (Upper) Shows multiple nucleotide polymorphisms in HTI4 and HTI5 including indels and nucleotide substitutions. (Lower) Describes the protein sequence polymorphism as a result of sequence changes. The differences leading to functional loss are marked as boxes.
Interestingly, there is another ABC transporter gene in chromosome 31 with gene id LDON_310017500.1 in HTI4 and gene id LDON_310017000.1 in HTI5 where a substitution at 5012th position (GGG ->GCA a non-synonymous substitution which leads to a change from G->A amino acid at 1671st position in HTI5 (Figure 5A). The same gene undergoes two substitutions at the 5010th position (from AGC ->AGT) (Figure 5B) leading to synonymous substitution. The Genbank L. donovani strain LdBPK282A1 gene corresponding to this gene has the ‘A’ variant in this locus (Figure 5C). This observation is also supported by our CNV analysis in chromosome 31 as well as by Imamura et al. (2016). This gene codes for pentamidine resistance transporter protein. Major facilitator Protein (MFS class) is responsible for solute transport via membrane. Recently their role in stress responses and virulence has been proposed by various groups (Shah et al., 2014; Zhang et al., 2014). One such gene in HTI4 has been pseudogenized and this gene is absent in HTI5 (assembly gaps can’t be ruled out as one of the plausible causes though). A small segment of the gene (297–360 bp) is copied at two places both in HTI5 and HTI4 with one point non-synonymous mutation (AAC->AGC) (gene id: LDON_290021200 in HTI4) resulting in N ->S (Figure 6). The presence of this duplicated domain and its role in gene regulation and pathogenicity may be intriguing. Functional members of this family is present in duplicates in HTI5 as well as HTI4 (HTI5: Chr29: 711070–711429; HTI5: Chr29: 710625–710970 and HTI4: chr29: 711414–711772; HTI4:chr29: 710661–711007). This indicates selection pressure is working on MFS classes of proteins. Studies on pathogenic organisms demonstrate that the acquisition of iron is very important in the intracellular pathogenesis process and Leishmania possesses molecular machinery for iron regulation (Huynh and Andrews, 2008). Interestingly, these transporter proteins have links to the MFS superfamily (Huynh et al., 2006) and insight into genetic structure, mechanism and regulation via domain duplication in this class of proteins appears intriguing in designing targeted therapies.
FIGURE 5. Synonymous and non-synonymous substitution in early and late passages leading to altered ABC transporter genes in chromosome 31. (A) A pair wise comparison between the protein coding genes in HTI4 and HTI5 showing single amino acid mismatches. (B) Nucleotide comparison in two genes indicates few substitutions, out of which the first one is a synonymous substitution and the second one is a non-synonymous substitution. A comparison between HTI4, HTI5 and Genbank strain (C) clearly indicates change of G–A in late passage same as the reference genome.
FIGURE 6. Schematic showing domain copy of MFS gene in Chromosome 29. A pseudo gene of MFS class of proteins present in HTI4 is missing in HTI5. However, a small domain from region 297–360 is being copied at two different locations in chromosome 29 with a non-synonymous substitution (AAC->AGC).
A fully functional acetyl-CoA synthatase gene in HTI4 has undergone modification [chr23: 199275–199598 in HTI5 and chr23: 198391–19851 in HTI4] leading to loss of function (Figure 7). There is also domain duplication in HTI4: LDON_230010500.1 from 1 to 118th position. The stringency in nutrient acquisition is key to preparedness of the parasite to intracellular environment (Carman et al., 2008) and loss of function of acetyl-CoA synthetase may restrict the parasite’s capacity to thrive in nutrient poor conditions using alternative carbon sources (McConville et al., 2015). Domain duplication in the early passage promastigotes (HTI4) indicate enhanced intracellular survival strategy.
FIGURE 7. Comparison of Acetyl Co-A synthetase genes in early and late passage. Gene coding for Acetyl CoA Synthetase undergoes several substitutions in the late passage which makes it non-functional.
Calpain like proteases play a very important role in infection process. A single insertion in 172519th position in HTI5 on chromosome 27 (Chr27: 172458–188744) leads to frameshift mutation causing this gene to become a pseudogene. The corresponding functional gene in early passage is present at Chr27: 172465–188750 (Figure 8). In case of Genbank L. donovani genome, the protein sequence is more identical to the HTI5 protein sequence. Recent reports have pointed toward the regulatory role played by expressed pseudogenes in cancer cells and parasites (Wen et al., 2011). To check the downstream impact of mutations in ABC transporters and CALPs, we checked the relative sensitivity of early and late passage promastigotes to the drug miltefosine. Early passage parasites displayed an IC50 of 30 μM whereas the late passage promastigotes showed increased apoptotic death under drug pressure and an IC50 of 18 μM (Supplementary Figure S3). Interestingly, parasite CALPs may serve other functions in the intracellular form which determine disease outcome and host responses (Branquinha et al., 2013) and are thus potent drug targets. This may open up new avenues in understanding Leishmania biology. The implications of these modifications need further investigation.
FIGURE 8. Frameshift mutation in Calpain like cysteine protease gene due to a single insertion in late passage. A single insertion in HTI5 (as well as Genbank LDBPK 282 A1) completely changes the coding frame of the gene (Upper panel) making it inactive. The changes in late passage and Genbank strain is consistent.
Comparative Transcriptome Profiles Between Passages Reveal Calpain Proteins Play an Important Role in Maintaining Virulence
The information from genomic studies done in this work as well as previous studies on axenic promastigotes indicated a strong link between parasite pre-adaptation and virulence. We hypothesized that continuous axenic cultivation may lead to altered transcript expression. To check for specific changes we analyzed the differential levels of transcript expression in three different culture stages (early, intermediate, and late passages) using RNAseq. Overall, there was no massive transcript expression switch between the early and late passages (Figure 9 and Supplementary Figure S1) although specific genes presented differential expression (Supplementary Data Sheet S4) which correlated to certain extent with the CNV. The relative gene expression of a few genes was validated by RT-PCR (Supplementary Figure S4). Among the genes coding for surface active proteins, there was down regulation in membrane bound acid phosphatase and surface antigen like protein transcript levels (Supplementary Data Sheet S4) in the later passages. Acid phosphatase activity is needed for virulence and is also involved in endosome sorting (Katakura and Kobayashi, 1988). This may also be an adaptive response to the hydrolytic environment encountered inside sandfly gut and mammalian cells (Papadaki et al., 2015). Cyclin-dependent kinase pho85-like protein is a morphogenesis related protein having strong link to virulence in Ustilago (Castillo-Lluva et al., 2007) although its function in Leishmania has not been elucidated to date; it is implicated in environmental signaling in yeast (Carroll and O’Shea, 2002). This gene was under-expressed in the later passages. Among the cytoskeletal and flagellar proteins, a paraflagellar rod protein 1D (PFR 1D) was over expressed in early passage compared to the intermediate passage consistent with previous proteomic data (Magalhaes et al., 2014) on differentially passaged promastigotes. PFR 1D is essential for proper flagellar motility, a marker of promastigote virulence (Ginger et al., 2008) and it seems that flagellar morphology gets affected after few passages and ultimately leads to loss of virulence. Interestingly, multiple members of calpain like cysteine proteases (CALPs) (XLOC_001409, XLOC_000608, XLOC_001023; Supplementary Data Sheet 4) presented a typical pattern of expression where they were down-regulated between early and intermediate passages and ultimately had increased expression in the final passages. Some of the calpains have lost protease domain to gain microtubule organization function like SMP-1 (Branquinha et al., 2013). SMP-1 is required in all promastigote stages and also for transformation of amastigotes to promastigotes for the development of flagella (Tull et al., 2010). Thus over expression of these atypical CALPs in the later passages indicate struggle of the promastigotes to convert to infective forms. The role of calpains in signal transduction and cytoskeletal remodeling is well documented in trypanosomatids (Galetovic et al., 2011), and structural deformities may affect the virulence function of the parasite.
FIGURE 9. Hierarchical Clustering Heatmap of differentially expressed genes in early, intermediate and late passages. Heat maps were generated to show the comparative log 2 abundance ratios between early, intermediate, and late passages using RNAseq data. Upregulated proteins are depicted by blue bars, downregulated proteins by red bars.
Metabolic reprogramming is a hallmark of Leishmania life-cycle where the parasite switches between two extreme environmental conditions. In the sandfly gut glucose and amino acids are the primary source of carbon while once inside the host, fatty acids and amino acids provide the necessary carbon for metabolism. Increased transcript abundance of 3-hydroxyisobutyryl-coenzyme a hydrolase-like protein in the early passages relative to the late is in consistency with the preparedness of the infection ready parasites to use beta-oxidation of fatty acids, as observed by others (Alcolea et al., 2009) (Supplementary Data Sheet S4). As opposed to the earlier observations (Moreno et al., 2014), level of tyrosine amino transferase presented a reverse trend, as it was slightly downregulated in the early passage, although protein expression may reveal different results as observed previously for Trypanosoma cruzi (dos Santos et al., 2012). Purine acquisition and inter-conversion is an important determinant of environmental sensing and morphogenesis in Leishmania. We observed a decrease in expression of a putative ribonucleoside-diphosphate reductase small chain, involved in this pathway in the late passage promastigotes in sync with earlier observation in L. donovani under purine stressed condition (Martin et al., 2014). On the other hand an ATP diphosphohydrolase was upregulated in the virulent parasites indicating preparedness for host-cell invasion and immunomodulation (Figueiredo et al., 2016) reflected in the differential host response in our infection studies. Amino acid transporter (aATP11) which has also been associated with purine starved condition (Martin et al., 2014) triggering parasite differentiation to metacyclic form (Inbar et al., 2017) was upregulated in the late passage, again reinstating the fact that these parasites lag in morphological transformation and try hard to attain infective status when growing in culture for a long time. Prostaglandin f2-alpha synthase involved in prostaglandin biosynthetic pathway is reduced in expression in the more infective promastigotes while the gene expression is more pronounced in the later passages indicating non-completion of metacyclogenesis as reported earlier (Araujo-Santos et al., 2014). Leishmania salvage folate and pteridines from their host or vector as they cannot synthesize them. The expression of these genes is higher in infective promastigotes and is known to aid in survival under extreme intra macrophagic conditions (Papadopoulou et al., 2002; Alcolea et al., 2010). As expected a gene for this protein on chromosome 6 was upregulated till 11th passage followed by a significant under expression in the late passage (Figure 9 and Supplementary Data Sheet S4). Altogether it signifies that nutrient availability is strongly linked to gene expression in Leishmania and that many metabolites and classical metabolic enzymes are capable of altering the state of cellular differentiation and development as has been extensively studied in cancer cells (van der Knaap and Verrijzer, 2016). Metallopeptidases are important virulence factors in Leishmania and we already observed differential regulation of the major surface peptidase gp63 at protein level. Here, we observed slightly upregulated expression of two peptidases a metallopeptidase and a thimet oligopeptidase in the late passage. Transcriptomic data from late logarithmic phase promastigotes of Trypanosoma brucei growing in normal and purine supplemented medium displayed under expression of thimet oligopeptidase A in the former system linking it to purine stressed condition when approaching metacyclic state (Fernandez-Moya et al., 2014) corroborating our finding.
As pointed out above, genetic and epigenetic control of gene expression is central to adaptation of the parasite to its immediate and prospective environment. Consistent with that we observed an increased abundance of a histone-like transcription factor (CBF/NF-Y) and archaeal histone transcript in early passage relative to the late passage. The histone modifying enzyme histone deacetylase was also over expressed in the early passage. This protein has been associated with the disease causing stage of the parasite (Vergnes et al., 2005) and its higher activity may be linked to more invasiveness in the parasite. The higher protein level expression of HDAC in early passage promastigotes was also confirmed by western blot on whole cell lysates (Supplementary Figure S5). Together with this, protein modifying enzyme farnesyltransferase presented a downward trend in the late passage (Figure 9 and Supplementary Data Sheet S4). This gene product has been associated with severe growth impairment in trypanosomatids (Gillespie et al., 2007). A putative pre-mRNA splicing factor ATP-dependent RNA helicase, member of DEAD box RNA helicase family presented increased expression till 11th passage. These proteins were shown to be strongly associated with mRNA of differentiating promastigotes as well as interactome of in vitro cultured amastigotes (Saxena et al., 2007). These proteins are mostly transiently upregulated or downregulated controlling the translational machinery. One of them has been associated with growth defects in L. infantum (Barhoumi et al., 2006). DNA repair related MutS-like proteins play important role in preventing genotoxicity and oxidative damage, a survival strategy adapted by many prokaryotic and eukaryotic pathogens (Genois et al., 2014), was upregulated in the virulent passage. Among the other differentially expressed transcripts are non-protein coding rRNA genes and conserved proteins with unknown functions. Non-coding RNAs along with proteins in the secreted exosomal milieu of Leishmania are known to modulate host–parasite interaction and changes in the composition of these exosomes may affect infectivity (Lambertz et al., 2015).
The results indicate that subtle genomic and transcriptomic variations drive the adaptation of promastigotes to culture condition and these changes mostly lead to loss of virulence in addition to morphological and metabolic modifications. Of particular interest is the role played by peptidases of different classes mainly metallopeptidases and calpain-like cysteine proteases in host–parasite interaction. Another important aspect is the importance of purine stress and modes of purine acquisition which modulate parasite and host responses to infection.
Significant Regulation of Virulence Proteins in Comparative Protein Expression Studies
Till date the loss of infectivity in late passage culture has been explained by reduced expression of virulence factors, particularly those which can interact with and modulate the host cells. Proteomic comparison of crude membrane extracts of L. donovani AG83 promastigotes (LAg) revealed altered protein expression, particularly a decrease in expression of some known virulence factors after continuous passaging (Table 3 and Supplementary Figure S6). There was very little overlap between transcriptomic and proteomic data which reinstates the fact that leishmanial gene expression is majorly controlled at post-transcriptional level (Lahav et al., 2011). A total of 35 differentially expressed spots were identified, out of which 31 spots presented higher expression in the early passage while four spots were highly expressed in the late passage. For example, gp63 (Olivier et al., 2012), HSP 70 (Ramirez et al., 2013a), elongation factor 1 alpha (Nandan et al., 2003), etc. have all been reported as virulence factors and were over expressed in the early passage. Western blot analysis on whole cell lysates of early and late passage promastigotes confirmed increased expression of gp63 in the virulent parasites (Supplementary Figure S5). The pyruvate dehydrogenase complex and its component dihydrolipoamide acetyltransferase has been associated with multiple functions including oxidative defense and regulation of gene expression in various pathogenic microbes (Spalding and Prigge, 2010), and presented an increased expression in the early passage. Apart from these, many hypothetical proteins also displayed differential protein expression and may determine the invasive capacity of the parasite. Transcript expression of each gene is also presented in Table 3 and presented a similar trend. Most of the proteins which displayed differential expression in our study form a part of the infectious exosomal proteins (Silverman et al., 2010). The changes in exosomal milieu probably changes the modulatory capacity of the parasites of late passage inside macrophages, which in turn make them more susceptible to microbicidal attack. The ultimate protein repertoire required for successful host–parasite interaction is the result of ordered changes in the genetic make-up and/or gene expression in the promastigotes.
TABLE 3. List of proteins differentially expressed (≥1.5-folds) between early and late passages identified by MALDI-ToF MS/MS.
An upregulated calpain like protein related to small myristoylated protein-1 (SMP-1) may also signify the struggle of late passage promastigotes to transform into infective metacyclic forms and as already mentioned in the previous section. It is also interesting to note that some calpain-like proteins, particularly those lacking the protease domain, are upregulated in the late passages while some others with protease activity are downregulated or psedogenized in the avirulent late passage parasites. Another consistently upregulated protein in the late passage was a putative ATPase beta subunit. Although there is no report in protozoans, a report linked increased acid sensitivity of bacterial pathogen to ATPase overexpression (McEntire et al., 2004). This may explain partly the non-viability of late passage parasites inside acidic host phagolysosomes. Beta tubulin was significantly up regulated in the non-infectious stage of the pathogen (Coulson et al., 1996). On the contrary alpha tubulin was mostly down regulated in the later passages. Tubulin plays a very important role in cytoskeleton formation and its expression is post-transcriptionally controlled and tightly linked to parasite transformation (Ramirez et al., 2013b). During the infection cycle, the pathogen utilizes hosts cytoskeletal machinery for pathogenicity (Gruenheid and Finlay, 2003) probably leading to less expression of these proteins.
Our study indicates that pathogenicity in L. donovani is a complex mechanism wherein the parasites are pre-adapted to a pathogenic lifestyle where nutritional stress and environmental sensing lead to global changes in both genes and gene expressions, ultimately leading to the production of effector molecules inside the host culminating in infection. The heteroxenic lifecycle and pleomorphic nature of the parasite adds another dimension to this process as successful transformation to metacyclic promastigotes determines virulence phenotype both in in vivo (other studies) in sandfly and in vitro (others and our study) in culture medium. It is re-established from these studies that long term cultivation of Leishmania promastigotes in axenic environment diminishes their adaptability toward intracellular life, and many of the factors involved in this process are markers of natural virulence attenuation. Leishmania have evolved several strategies to sense their immediate environment individually as well as collectively and respond to the same by controlling gene expression, and ultimately cell division and transformation. Transporter proteins are known to be involved in this pathway in other pathogens (Hlavacek et al., 2009). Continuous axenic cultivation possibly reverses this pathoadaptive response by inducing mutations in these genes which leads to loss of virulence. Further, pseudogenization of calpain like proteins involved in signal transduction and parasite transformation results in growth arrest at non-infective stages. Studies on in vitro cultured drug resistant and sensitive clinical isolates of L. donovani provided evidence that the former system has added advantage of increased metacyclogenesis in culture and higher infectivity index in host cells (Vanaerschot et al., 2010), which is also reflected in the miltefosine sensitivity assay done on the early and late passage promastigotes. Moreover, early and late passage promastigotes display differential infectivity (this work) because interaction of virulent L. donovani with host macrophages triggers differential transcriptional modulation in the latter as compared to avirulent strains (manuscript communicated elsewhere). A parallel gene expression study in the intracellular amastigotes also pointed toward higher expression of transporter and calcium-dependent cysteine-type endo peptidase activity genes (GO: 0005215; GO: 0004198) in more virulent forms (manuscript communicated elsewhere), strengthening the relevance of pre-adaptation in virulent parasites. The most important factor responsible for adaptive changes is the nutritional status of the immediate environment, sensing which the promastigotes undergo morphological and metabolic switch. This phenomenon is also strongly connected with transporter system in parasitic protists (Dean et al., 2014). In our genetic studies, we found domain duplication in acetyl CoA synthetase gene in the early passage promastigotes. Recent works have highlighted the importance of this enzyme in L. donovani infectivity owing to lipid and ergosterol biosynthesis (Soumya et al., 2017). On a broader note, this enzyme and its product acetyl CoA has been linked to global gene expression and chromatin regulation (van der Knaap and Verrijzer, 2016). This supports the theory of differential gene expression, and epigenetic changes in the virulent versus avirulent parasites, and reiterates the fact that nutrition and transcription are tightly linked. Another facet of the study relates to the importance of DNA repair and recombination activity prevalent in these parasites which promote antigenic variation suitable for intracellular adaptability and pathogenicity and needs further delving. This was indicated by differential regulation of MutS like protein at transcriptional level which has been implicated in quorum sensing independent developmental regulation in Trypanosoma recently (Zimmermann et al., 2017). Nevertheless, this in vitro cultured parasite model of gene regulation and virulence provides a basis for understanding not merely parasite adaptation to culture conditions and its genetic basis but also underlines the fact that common pathways are involved in intraspecific and intrastrain variations observed in Leishmania. Most of the analysis done so far in literature is based on comparison of several clinical isolates in Leishmania. However, comparison of various cultured passages of promastigotes are reported here.
Interestingly, in this study, we report a large number of ABC transporters undergoing structural changes in later passages, thus their inactivation may be a possible cause of loss of virulence. We have also found that in addition to several important genes such as GP63 and heat shock proteins which are already known to play major role in virulence, calpain like cysteine proteases may have significant role in loss of virulence since the genes undergo alteration and the expression levels change.
This work re-establishes the findings of other studies at transcriptomic and proteomic level to exercise caution while using serially passaged promastigotes for infectivity and therapeutic studies but also adds genetic analysis to this list. The revelation from our study that the NCBI reference strain was derived from clinical isolate which was cultured in vitro for multiple passages reiterates that subtle changes in the genome due to axenic cultivation may have important implications on the virulence phenotype. Moreover, similar genomic modifications may come into play when Leishmania infects different hosts, different vectors or even different sites of the same host which leads to differential host responses (Elso et al., 2004). Thus within-host or even within-vector selection pressure appears to be the major factor driving the modifications in the parasite which make them adaptable to their environment (Shaw, 1997), and is accompanied by no major gene losses but rather subtle polymorphisms that possibly alter the functionality of the genes. It is further evidenced from the higher expression of transporter and calpain-like proteins in visceralizing L. donovani strains from Sri Lanka (Zhang et al., 2014) compared to strains causing cutaneous lesions. Virulence is a function of nutritional status of the parasites and more time spent in the nutrient rich culture medium diminishes their pathoadaptive characteristics. Significant is the fact that changes in genes and gene expression occur in same/similar set of genes both in culture and in nature for virulence attenuation and cultured promastigotes can well mimic their sandfly counterparts but long-term cultures should be avoided. Changes in transporters and certain house-keeping genes can be markers of virulence. Additionally the new class of calpain peptidases are emerging as important virulence determinants along with more established metallopeptidases and cathepsin- like cysteine peptidases.
Conclusion
The work presented here highlights the importance of newly emerging pathoadaptive factors like transporters and calpain-like proteins in Leishmania virulence. Our work can add to the pool of information on adaptive genomics in L. donovani and careful mining of these genes can be used for virulence surveillance at least throughout the Indian subcontinent. The genome sequences of early and late passage promastigotes can act as references for future genomic studies on Leishmania, particularly related to virulence and drug resistance. Annotated leads from this study and further annotation and functional characterization of hypothetical proteins can provide novel drug/vaccine targets for disease management.
Author Contributions
RS and NA conceived the study. RS, R, NA, and ST designed the experiments. RS, ST, MC, SdD, and R analyzed the data. RS, ST, MC, R, and NA wrote the manuscript. RS, R, SD, and MS performed and optimized the experiments. RS, MC, SD, and ST prepared figures. RC contributed equipment and reagents. All authors critically revised and approved the manuscript.
Funding
This work was funded by Council of Scientific and Industrial Research (CSIR), Government of India Network Project BSC0114. RS, R, MS, and SD received fellowships from CSIR. MC received fellowship from Department of Science and Technology and SdD from Department of Biotechnology, Government of India.
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.
The handling Editor declared his shared affiliation with the reviewers MB and NM.
Acknowledgments
We acknowledge Mr. Sandip Chakraborty for Mass Spectrometry support and Mr. Janmenjay Midya for assistance in animal experiments from CSIR-IICB. We are thankful to Mr. Samrat Ghosh and Mr. Arijit Panda from Computational Biology and Genomics Lab, Structural Biology and Bioinformatics Division, CSIR-IICB for setting up the internal Blast Page. The assistance of Ms. Shreya Dev during experiments is duly acknowledged.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.01279/full#supplementary-material
Footnotes
References
Abyzov, A., Urban, A. E., Snyder, M., and Gerstein, M. (2011). CNVnator: an approach to discover, genotype, and characterize typical and atypical CNVs from family and population genome sequencing. Genome Res. 21, 974–984. doi: 10.1101/gr.114876.110
Alcolea, P. J., Alonso, A., Gomez, M. J., Moreno, I., Dominguez, M., Parro, V., et al. (2010). Transcriptomics throughout the life cycle of Leishmania infantum: high down-regulation rate in the amastigote stage. Int. J. Parasitol. 40, 1497–1516. doi: 10.1016/j.ijpara.2010.05.013
Alcolea, P. J., Alonso, A., Sanchez-Gorostiaga, A., Moreno-Paz, M., Gomez, M. J., Ramos, I., et al. (2009). Genome-wide analysis reveals increased levels of transcripts related with infectivity in peanut lectin non-agglutinated promastigotes of Leishmania infantum. Genomics 93, 551–564. doi: 10.1016/j.ygeno.2009.01.007
Ali, K. S., Rees, R. C., Terrell-Nield, C., and Ali, S. A. (2013). Virulence loss and amastigote transformation failure determine host cell responses to Leishmania mexicana. Parasite Immunol. 35, 441–456. doi: 10.1111/pim.12056
Altschul, S. F., Gertz, E. M., Agarwala, R., Schaffer, A. A., and Yu, Y. K. (2009). PSI-BLAST pseudocounts and the minimum description length principle. Nucleic Acids Res. 37, 815–824. doi: 10.1093/nar/gkn981
Araujo-Santos, T., Rodriguez, N. E., Moura-Pontes, S., Dixt, U. G., Abanades, D. R., Bozza, P. T., et al. (2014). Role of prostaglandin F2alpha production in lipid bodies from Leishmania infantum chagasi: insights on virulence. J. Infect. Dis. 210, 1951–1961. doi: 10.1093/infdis/jiu299
Banerjee, A., De, M., and Ali, N. (2008). Complete cure of experimental visceral leishmaniasis with amphotericin B in stearylamine-bearing cationic liposomes involves down-regulation of IL-10 and favorable T cell responses. J. Immunol. 181, 1386–1398. doi: 10.4049/jimmunol.181.2.1386
Barhoumi, M., Tanner, N. K., Banroques, J., Linder, P., and Guizani, I. (2006). Leishmania infantum LeIF protein is an ATP-dependent RNA helicase and an eIF4A-like factor that inhibits translation in yeast. FEBS J. 273, 5086–5100. doi: 10.1111/j.1742-4658.2006.05506.x
Bendtsen, J. D., Jensen, L. J., Blom, N., Von Heijne, G., and Brunak, S. (2004). Feature-based prediction of non-classical and leaderless protein secretion. Protein Eng. Des. Sel. 17, 349–356. doi: 10.1093/protein/gzh037
Bhowmick, S., Mazumdar, T., Sinha, R., and Ali, N. (2010). Comparison of liposome based antigen delivery systems for protection against Leishmania donovani. J. Control. Release 141, 199–207. doi: 10.1016/j.jconrel.2009.09.018
Bradford, M. M. (1976). A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem. 72, 248–254. doi: 10.1016/0003-2697(76)90527-3
Branquinha, M. H., Marinho, F. A., Sangenito, L. S., Oliveira, S. S., Goncalves, K. C., Ennes-Vidal, V., et al. (2013). Calpains: potential targets for alternative chemotherapeutic intervention against human pathogenic trypanosomatids. Curr. Med. Chem. 20, 3174–3185. doi: 10.2174/0929867311320250010
Carman, A. J., Vylkova, S., and Lorenz, M. C. (2008). Role of acetyl coenzyme A synthesis and breakdown in alternative carbon source utilization in Candida albicans. Eukaryot. Cell 7, 1733–1741. doi: 10.1128/EC.00253-08
Carroll, A. S., and O’Shea, E. K. (2002). Pho85 and signaling environmental conditions. Trends Biochem. Sci. 27, 87–93. doi: 10.1016/S0968-0004(01)02040-0
Castillo-Lluva, S., Alvarez-Tabares, I., Weber, I., Steinberg, G., and Perez-Martin, J. (2007). Sustained cell polarity and virulence in the phytopathogenic fungus Ustilago maydis depends on an essential cyclin-dependent kinase from the Cdk5/Pho85 family. J. Cell Sci. 120, 1584–1595. doi: 10.1242/jcs.005314
Coulson, R. M., Connor, V., Chen, J. C., and Ajioka, J. W. (1996). Differential expression of Leishmania major beta-tubulin genes during the acquisition of promastigote infectivity. Mol. Biochem. Parasitol. 82, 227–236. doi: 10.1016/0166-6851(96)02739-9
da Silva, R., and Sacks, D. L. (1987). Metacyclogenesis is a major determinant of Leishmania promastigote virulence and attenuation. Infect. Immun. 55, 2802–2806.
Dean, P., Major, P., Nakjang, S., Hirt, R. P., and Embley, T. M. (2014). Transport proteins of parasitic protists and their role in nutrient salvage. Front. Plant Sci. 5:153. doi: 10.3389/fpls.2014.00153
Dey, T., Afrin, F., Anam, K., and Ali, N. (2002). Infectivity and virulence of Leishmania donovani promastigotes: a role for media, source, and strain of parasite. J. Eukaryot. Microbiol. 49, 270–274. doi: 10.1111/j.1550-7408.2002.tb00369.x
dos Santos, P. F., Ruiz, J. C., Soares, R. P., Moreira, D. S., Rezende, A. M., Folador, E. L., et al. (2012). Molecular characterization of the hexose transporter gene in benznidazole resistant and susceptible populations of Trypanosoma cruzi. Parasit. Vectors 5:161. doi: 10.1186/1756-3305-5-161
Downing, T., Imamura, H., Decuypere, S., Clark, T. G., Coombs, G. H., Cotton, J. A., et al. (2011). Whole genome sequencing of multiple Leishmania donovani clinical isolates provides insights into population structure and mechanisms of drug resistance. Genome Res. 21, 2143–2156. doi: 10.1101/gr.123430.111
Dumetz, F., Imamura, H., Sanders, M., Seblova, V., Myskova, J., Pescher, P., et al. (2017). Modulation of aneuploidy in Leishmania donovani during adaptation to different in vitro and in vivo environments and its impact on gene expression. mBio 8:e00599-17. doi: 10.1128/mBio.00599-17
Elso, C. M., Roberts, L. J., Smyth, G. K., Thomson, R. J., Baldwin, T. M., Foote, S. J., et al. (2004). Leishmaniasis host response loci (lmr1-3) modify disease severity through a Th1/Th2-independent pathway. Genes Immun. 5, 93–100. doi: 10.1038/sj.gene.6364042
Fernandez-Moya, S. M., Carrington, M., and Estevez, A. M. (2014). A short RNA stem-loop is necessary and sufficient for repression of gene expression during early logarithmic phase in trypanosomes. Nucleic Acids Res. 42, 7201–7209. doi: 10.1093/nar/gku358
Figueiredo, A. B., Souza-Testasicca, M. C., and Afonso, L. C. (2016). Purinergic signaling and infection by Leishmania: a new approach to evasion of the immune response. Biomed. J. 39, 244–250. doi: 10.1016/j.bj.2016.08.004
Galetovic, A., Souza, R. T., Santos, M. R., Cordero, E. M., Bastos, I. M., Santana, J. M., et al. (2011). The repetitive cytoskeletal protein H49 of Trypanosoma cruzi is a calpain-like protein located at the flagellum attachment zone. PLoS One 6:e27634. doi: 10.1371/journal.pone.0027634
Genois, M. M., Paquet, E. R., Laffitte, M. C., Maity, R., Rodrigue, A., Ouellette, M., et al. (2014). DNA repair pathways in trypanosomatids: from DNA repair to drug resistance. Microbiol. Mol. Biol. Rev. 78, 40–73. doi: 10.1128/MMBR.00045-13
Gillespie, J. R., Yokoyama, K., Lu, K., Eastman, R. T., Bollinger, J. G., Van Voorhis, W. C., et al. (2007). C-terminal proteolysis of prenylated proteins in trypanosomatids and RNA interference of enzymes required for the post-translational processing pathway of farnesylated proteins. Mol. Biochem. Parasitol. 153, 115–124. doi: 10.1016/j.molbiopara.2007.02.009
Ginger, M. L., Portman, N., and Mckean, P. G. (2008). Swimming with protists: perception, motility and flagellum assembly. Nat. Rev. Microbiol. 6, 838–850. doi: 10.1038/nrmicro2009
Glavinas, H., Krajcsi, P., Cserepes, J., and Sarkadi, B. (2004). The role of ABC transporters in drug resistance, metabolism and toxicity. Curr. Drug Deliv. 1, 27–42. doi: 10.2174/1567201043480036
Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Trinity: reconstructing a full-length transcriptome without a genome from RNA-Seq data. Nat. Biotechnol. 29:644. doi: 10.1186/1471-2105-12-S14-S2
Gruenheid, S., and Finlay, B. B. (2003). Microbial pathogenesis and cytoskeletal function. Nature 422, 775–781. doi: 10.1038/nature01603
Haas, B. J., Papanicolaou, A., Yassour, M., Grabherr, M., Blood, P. D., Bowden, J., et al. (2013). De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat. Protoc. 8, 1494–1512. doi: 10.1038/nprot.2013.084
Hlavacek, O., Kucerova, H., Harant, K., Palkova, Z., and Vachova, L. (2009). Putative role for ABC multidrug exporters in yeast quorum sensing. FEBS Lett. 583, 1107–1113. doi: 10.1016/j.febslet.2009.02.030
Huynh, C., and Andrews, N. W. (2008). Iron acquisition within host cells and the pathogenicity of Leishmania. Cell. Microbiol. 10, 293–300. doi: 10.1111/j.1462-5822.2007.01095.x
Huynh, C., Sacks, D. L., and Andrews, N. W. (2006). A Leishmania amazonensis ZIP family iron transporter is essential for parasite replication within macrophage phagolysosomes. J. Exp. Med. 203, 2363–2375. doi: 10.1084/jem.20060559
Iantorno, S. A., Durrant, C., Khan, A., Sanders, M. J., Beverley, S. M., Warren, W. C., et al. (2017). Gene expression in Leishmania is regulated predominantly by gene dosage. mBio 8:e01393-17. doi: 10.1128/mBio.01393-17
Imamura, H., Downing, T., Van Den Broeck, F., Sanders, M. J., Rijal, S., Sundar, S., et al. (2016). Evolutionary genomics of epidemic visceral leishmaniasis in the Indian subcontinent. eLife 5:e12613. doi: 10.7554/eLife.12613
Inbar, E., Hughitt, V. K., Dillon, L. A., Ghosh, K., El-Sayed, N. M., and Sacks, D. L. (2017). The transcriptome of Leishmania major developmental stages in their natural sand fly vector. mBio 8:e00029-17. doi: 10.1128/mBio.00029-17
Jackson, A. P., Otto, T. D., Aslett, M., Armstrong, S. D., Bringaud, F., Schlacht, A., et al. (2016). Kinetoplastid phylogenomics reveals the evolutionary innovations associated with the origins of parasitism. Curr. Biol. 26, 161–172. doi: 10.1016/j.cub.2015.11.055
Katakura, K., and Kobayashi, A. (1988). Acid phosphatase activity of virulent and avirulent clones of Leishmania donovani promastigotes. Infect. Immun. 56, 2856–2860.
Kent, W. J. (2002). BLAT–the BLAST-like alignment tool. Genome Res. 12, 656–664. doi: 10.1101/gr.229202
Kim, D., Pertea, G., Trapnell, C., Pimentel, H., Kelley, R., and Salzberg, S. L. (2013). TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 14, 1–13. doi: 10.1186/gb-2013-14-4-r36
Kurtz, S., Phillippy, A., Delcher, A. L., Smoot, M., Shumway, M., Antonescu, C., et al. (2004). Versatile and open software for comparing large genomes. Genome Biol. 5:R12. doi: 10.1186/gb-2004-5-2-r12
Laffitte, M. N., Leprohon, P., Papadopoulou, B., and Ouellette, M. (2016). Plasticity of the Leishmania genome leading to gene copy number variations and drug resistance. F1000Res. 5:2350. doi: 10.12688/f1000research.9218.1
Lahav, T., Sivam, D., Volpin, H., Ronen, M., Tsigankov, P., Green, A., et al. (2011). Multiple levels of gene regulation mediate differentiation of the intracellular pathogen Leishmania. FASEB J. 25, 515–525. doi: 10.1096/fj.10-157529
Lambertz, U., Oviedo Ovando, M. E., Vasconcelos, E. J., Unrau, P. J., Myler, P. J., and Reiner, N. E. (2015). Small RNAs derived from tRNAs and rRNAs are highly enriched in exosomes from both old and new world Leishmania providing evidence for conserved exosomal RNA Packaging. BMC Genomics 16:151. doi: 10.1186/s12864-015-1260-7
Langmead, B., Trapnell, C., Pop, M., and Salzberg, S. L. (2009). Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 10, 1–10. doi: 10.1186/gb-2009-10-3-r25
Lei, S. M., Romine, N. M., and Beetham, J. K. (2010). Population changes in Leishmania chagasi promastigote developmental stages due to serial passage. J. Parasitol. 96, 1134–1138. doi: 10.1645/GE-2566.1
Leprohon, P., Legare, D., Raymond, F., Madore, E., Hardiman, G., Corbeil, J., et al. (2009). Gene expression modulation is associated with gene amplification, supernumerary chromosomes and chromosome loss in antimony-resistant Leishmania infantum. Nucleic Acids Res. 37, 1387–1399. doi: 10.1093/nar/gkn1069
Li, H., and Durbin, R. (2009). Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics 25, 1754–1760. doi: 10.1093/bioinformatics/btp324
Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., et al. (2009). The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078–2079. doi: 10.1093/bioinformatics/btp352
Li, H., Ruan, J., and Durbin, R. (2008). Mapping short DNA sequencing reads and calling variants using mapping quality scores. Genome Res. 18, 1851–1858. doi: 10.1101/gr.078212.108
Magalhaes, R. D., Duarte, M. C., Mattos, E. C., Martins, V. T., Lage, P. S., Chavez-Fumagalli, M. A., et al. (2014). Identification of differentially expressed proteins from Leishmania amazonensis associated with the loss of virulence of the parasites. PLoS Negl. Trop. Dis. 8:e2764. doi: 10.1371/journal.pntd.0002764
Martin, J. L., Yates, P. A., Soysa, R., Alfaro, J. F., Yang, F., Burnum-Johnson, K. E., et al. (2014). Metabolic reprogramming during purine stress in the protozoan pathogen Leishmania donovani. PLoS Pathog. 10:e1003938. doi: 10.1371/journal.ppat.1003938
McCall, L. I., Zhang, W. W., and Matlashewski, G. (2013). Determinants for the development of visceral leishmaniasis disease. PLoS Pathog. 9:e1003053. doi: 10.1371/journal.ppat.1003053
McConville, M. J., Saunders, E. C., Kloehn, J., and Dagley, M. J. (2015). Leishmania carbon metabolism in the macrophage phagolysosome- feast or famine? F1000Res. 4:938. doi: 10.12688/f1000research.6724.1
McEntire, J. C., Carman, G. M., and Montville, T. J. (2004). Increased ATPase activity is responsible for acid sensitivity of nisin-resistant Listeria monocytogenes ATCC 700302. Appl. Environ. Microbiol. 70, 2717–2721. doi: 10.1128/AEM.70.5.2717-2721.2004
Moreira, D., Santarem, N., Loureiro, I., Tavares, J., Silva, A. M., Amorim, A. M., et al. (2012). Impact of continuous axenic cultivation in Leishmania infantum virulence. PLoS Negl. Trop. Dis. 6:e1469. doi: 10.1371/journal.pntd.0001469
Moreno, M. A., Alonso, A., Alcolea, P. J., Abramov, A., De Lacoba, M. G., Abendroth, J., et al. (2014). Tyrosine aminotransferase from Leishmania infantum: a new drug target candidate. Int. J. Parasitol. Drugs Drug Resist. 4, 347–354. doi: 10.1016/j.ijpddr.2014.06.001
Moriya, Y., Itoh, M., Okuda, S., Yoshizawa, A. C., and Kanehisa, M. (2007). KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 35, W182–W185. doi: 10.1093/nar/gkm321
Nandan, D., Cherkasov, A., Sabouti, R., Yi, T., and Reiner, N. E. (2003). Molecular cloning, biochemical and structural analysis of elongation factor-1 alpha from Leishmania donovani: comparison with the mammalian homologue. Biochem. Biophys. Res. Commun. 302, 646–652. doi: 10.1016/S0006-291X(03)00216-X
Olivier, M., Atayde, V. D., Isnard, A., Hassani, K., and Shio, M. T. (2012). Leishmania virulence factors: focus on the metalloprotease GP63. Microbes Infect. 14, 1377–1389. doi: 10.1016/j.micinf.2012.05.014
Papadaki, A., Politou, A. S., Smirlis, D., Kotini, M. P., Kourou, K., Papamarcaki, T., et al. (2015). The Leishmania donovani histidine acid ecto-phosphatase LdMAcP: insight into its structure and function. Biochem. J. 467, 473–486. doi: 10.1042/BJ20141371
Papadopoulou, B., Roy, G., Breton, M., Kundig, C., Dumas, C., Fillion, I., et al. (2002). Reduced infectivity of a Leishmania donovani biopterin transporter genetic mutant and its use as an attenuated strain for vaccination. Infect. Immun. 70, 62–68. doi: 10.1128/IAI.70.1.62-68.2002
Pescher, P., Blisnick, T., Bastin, P., and Spath, G. F. (2011). Quantitative proteome profiling informs on phenotypic traits that adapt Leishmania donovani for axenic and intracellular proliferation. Cell. Microbiol. 13, 978–991. doi: 10.1111/j.1462-5822.2011.01593.x
Quevillon, E., Silventoinen, V., Pillai, S., Harte, N., Mulder, N., Apweiler, R., et al. (2005). InterProScan: protein domains identifier. Nucleic Acids Res. 33, W116–W120. doi: 10.1093/nar/gki442
Ramirez, C. A., Dea-Ayuela, M. A., Gutierrez-Blazquez, M. D., Bolas-Fernandez, F., Requena, J. M., and Puerta, C. J. (2013a). Identification of proteins interacting with HSP70 mRNAs in Leishmania braziliensis. J. Proteomics 94, 124–137. doi: 10.1016/j.jprot.2013.09.008
Ramirez, C. A., Requena, J. M., and Puerta, C. J. (2013b). Alpha tubulin genes from Leishmania braziliensis: genomic organization, gene structure and insights on their expression. BMC Genomics 14:454. doi: 10.1186/1471-2164-14-454
Ribeiro, F. J., Przybylski, D., Yin, S., Sharpe, T., Gnerre, S., Abouelleil, A., et al. (2012). Finished bacterial genomes from shotgun sequence data. Genome Res. 22, 2270–2277. doi: 10.1101/gr.141515.112
Saxena, A., Lahav, T., Holland, N., Aggarwal, G., Anupama, A., Huang, Y., et al. (2007). Analysis of the Leishmania donovani transcriptome reveals an ordered progression of transient and permanent changes in gene expression during differentiation. Mol. Biochem. Parasitol. 152, 53–65. doi: 10.1016/j.molbiopara.2006.11.011
Shah, A. H., Singh, A., Dhamgaye, S., Chauhan, N., Vandeputte, P., Suneetha, K. J., et al. (2014). Novel role of a family of major facilitator transporters in biofilm development and virulence of Candida albicans. Biochem. J. 460, 223–235. doi: 10.1042/BJ20140010
Shaw, J. (1997). Ecological and evolutionary pressures on leishmanial parasites. Braz. J. Genet. 20, 123–128. doi: 10.1590/S0100-84551997000100021
Silverman, J. M., Clos, J., De’oliveira, C. C., Shirvani, O., Fang, Y., Wang, C., et al. (2010). An exosome-based secretion pathway is responsible for protein export from Leishmania and communication with macrophages. J. Cell Sci. 123, 842–852. doi: 10.1242/jcs.056465
Simão, F. A., Waterhouse, R. M., Ioannidis, P., Kriventseva, E. V., and Zdobnov, E. M. (2015). BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31, 3210–3212. doi: 10.1093/bioinformatics/btv351
Soumya, N., Panara, M. N., Neerupudi, K. B., and Singh, S. (2017). Functional analysis of an AMP forming acetyl CoA synthetase from Leishmania donovani by gene overexpression and targeted gene disruption approaches. Parasitol. Int. 66, 992–1002. doi: 10.1016/j.parint.2016.11.001
Spalding, M. D., and Prigge, S. T. (2010). Lipoic acid metabolism in microbial pathogens. Microbiol. Mol. Biol. Rev. 74, 200–228. doi: 10.1128/MMBR.00008-10
Stanke, M., Diekhans, M., Baertsch, R., and Haussler, D. (2008). Using native and syntenically mapped cDNA alignments to improve de novo gene finding. Bioinformatics 24, 637–644. doi: 10.1093/bioinformatics/btn013
Steinbiss, S., Silva-Franco, F., Brunk, B., Foth, B., Hertz-Fowler, C., Berriman, M., et al. (2016). Companion: a web server for annotation and analysis of parasite genomes. Nucleic Acids Res. 44, W29–W34. doi: 10.1093/nar/gkw292
Sterkers, Y., Lachaud, L., Bourgeois, N., Crobu, L., and Bastien, P. (2012). Novel insights into genome plasticity in Eukaryotes: mosaic aneuploidy in Leishmania. Mol. Microbiol. 86, 15–23. doi: 10.1111/j.1365-2958.2012.08185.x
Trapnell, C., Roberts, A., Goff, L., Pertea, G., Kim, D., Kelley, D. R., et al. (2012). Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 7, 562–578. doi: 10.1038/nprot.2012.016
Tull, D., Naderer, T., Spurck, T., Mertens, H. D., Heng, J., Mcfadden, G. I., et al. (2010). Membrane protein SMP-1 is required for normal flagellum function in Leishmania. J. Cell Sci. 123, 544–554. doi: 10.1242/jcs.059097
Van der Auwera, G. A., Carneiro, M. O., Hartl, C., Poplin, R., Del Angel, G., Levy-Moonshine, A., et al. (2013). From FastQ data to high-confidence variant calls: the genome analysis toolkit best practices pipeline. Curr. Protoc. Bioinformatics 43, 11.10.11–11.10.33. doi: 10.1002/0471250953.bi1110s43
van der Knaap, J. A., and Verrijzer, C. P. (2016). Undercover: gene control by metabolites and metabolic enzymes. Genes Dev. 30, 2345–2369. doi: 10.1101/gad.289140.116
Vanaerschot, M., Maes, I., Ouakad, M., Adaui, V., Maes, L., De Doncker, S., et al. (2010). Linking in vitro and in vivo survival of clinical Leishmania donovani strains. PLoS One 5:e12211. doi: 10.1371/journal.pone.0012211
Vergnes, B., Sereno, D., Tavares, J., Cordeiro-Da-Silva, A., Vanhille, L., Madjidian-Sereno, N., et al. (2005). Targeted disruption of cytosolic SIR2 deacetylase discloses its essential role in Leishmania survival and proliferation. Gene 363, 85–96. doi: 10.1016/j.gene.2005.06.047
Wen, Y. Z., Zheng, L. L., Liao, J. Y., Wang, M. H., Wei, Y., Guo, X. M., et al. (2011). Pseudogene-derived small interference RNAs regulate gene expression in African Trypanosoma brucei. Proc. Natl. Acad. Sci. U.S.A. 108, 8345–8350. doi: 10.1073/pnas.1103894108
Wiese, M. (1998). A mitogen-activated protein (MAP) kinase homologue of Leishmania mexicana is essential for parasite survival in the infected host. EMBO J. 17, 2619–2628. doi: 10.1093/emboj/17.9.2619
Yadav, A., Chandra, U., and Saha, S. (2016). Histone acetyltransferase HAT4 modulates navigation across G2/M and re-entry into G1 in Leishmania donovani. Sci. Rep. 6:27510. doi: 10.1038/srep27510
Zhang, W. W., Ramasamy, G., Mccall, L. I., Haydock, A., Ranasinghe, S., Abeygunasekara, P., et al. (2014). Genetic analysis of Leishmania donovani tropism using a naturally attenuated cutaneous strain. PLoS Pathog. 10:e1004244. doi: 10.1371/journal.ppat.1004244
Keywords: Leishmania donovani, genomics, in vitro passage, promastigotes, genome plasticity
Citation: Sinha R, C MM, Raghwan, Das S, Das S, Shadab M, Chowdhury R, Tripathy S and Ali N (2018) Genome Plasticity in Cultured Leishmania donovani: Comparison of Early and Late Passages. Front. Microbiol. 9:1279. doi: 10.3389/fmicb.2018.01279
Received: 19 March 2018; Accepted: 25 May 2018;
Published: 03 July 2018.
Edited by:
Renato A. Mortara, Federal University of São Paulo, BrazilReviewed by:
Nilmar Silvio Moretti, Federal University of São Paulo, BrazilMarcelo R. S. Briones, Federal University of São Paulo, Brazil
Lucile Maria Floeter-Winter, Universidade de São Paulo, Brazil
Copyright © 2018 Sinha, C, Raghwan, Das, Das, Shadab, Chowdhury, Tripathy and Ali. 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) and the copyright owner(s) 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.
*Correspondence: Sucheta Tripathy, tsucheta@iicb.res.in; tsucheta@gmail.com Nahid Ali, nali@iicb.res.in
†Present address: Roma Sinha, Mater Medical Research Institute, Translational Research Institute, Woolloongabba, Brisbane, QLD, Australia
Raghwan, Division of Bacteriology, National Institute of Cholera and Enteric Diseases, Kolkata, India
Mohammad Shadab, Department of Dermatology, School of Medicine, The University of Alabama, Birmingham, AL, United States
‡These authors have contributed equally to this work.