1. Introduction
Leishmania parasites are trypanosomatid protozoans transmitted through the bite of sandflies. During their life cycle, these protozoans alternate between extracellular flagellated promastigotes in sandfly gut and intracellular non-flagellated amastigotes in mammalian macrophages.
Leishmania parasites cause a group of diseases in humans with symptoms ranging from self-limiting skin lesions (cutaneous leishmaniasis) to life-threating infections of the bone marrow, liver, and spleen (visceral leishmaniasis) [
1]. Current treatments are plagued with high toxicity, low efficacy, and increasing resistance [
2]. Thus, new drugs are needed to control leishmaniasis.
The sterol biosynthesis pathway presents multiple targets for antileishmanial drugs. Among them is sterol C14-demethylase (C14DM, also known as CYP51) which is the primary target of azoles [
3,
4]. C14DM is a heme-containing enzyme catalyzing the oxidative removal of the C14-methyl group from sterol intermediates (
Figure 1)[
5]. Genetic or chemical inhibition of C14DM in
Leishmania leads to increased plasma membrane fluidity, hypersensitivity to heat, and compromised mitochondrial functions [
6,
7]. In addition to C14-demethylation, two more methyl groups need to be removed from the C4 position during sterol synthesis (
Figure 1). In
Saccharomyces cerevisiae, three enzymes catalyze sterol C4-demethylation: an O
2-dependent C4 sterol methyl oxidase (ERG25) that sequentially generates hydroxy, aldehyde and carboxylate intermediates; a C4 decarboxylase (ERG26) removing CO
2 from carboxylate intermediates to yield 3-ketosterol intermediates; and a 3-ketosterol reductase (ERG27) which regenerates the C3-hydroxy group [
8,
9,
10]. While
S. cerevisiae and vertebrates use the same C4 sterol methyl oxidase (ERG25) on both C4 methyl groups, plants have two distinct ERG25 homologs that oxidize 4,4-dimethylsterols and 4α-methylsterols, respectively [
11,
12] (
Figure 1).
A recent study has identified
Leishmania donovani CYP5122A1 (XP_003861867.1) as a non-canonical sterol C4-methyl oxidase [
13]. CYP5122A1 is well conserved among multiple trypanosomatid species. Recombinant LdCYP5122A1 protein binds and converts C4-methylated sterols into oxidized metabolites in the presence of NADPH and oxygen. Importantly, CYP5122A1 is essential for both promastigotes and intracellular amastigotes in
L. donovani. Heterozygous CYP5122A1 mutants (half knockout) in which one of the two
CYP5122A1 alleles is deleted showed reduced infectivity in mice [
13,
14]. In addition, inhibitors of both CYP5122A1 and C14DM display superior antileishmanial activity against
L. donovani in comparison to C14DM-selective inhibitors [
13,
15]. These findings suggest that C4DM is a promising drug target.
In addition to CYP5122A1, Leishmania parasites possess a putative ERG25 in their genomes. It is not known whether the leishmanial ERG25 protein is involved in the C4-demethylation reaction or whether its function overlaps with CYP5122A1. In this study, we investigated the roles of ERG25 in growth, virulence and sterol synthesis by characterizing genetic mutants of ERG25 in Leishmania major. We also explored the essentiality of CYP5122A1 in L. major. Our findings argue that CYP5122A1 is the sole sterol C4-methyl oxidase in Leishmania.
3. Discussion
Sterol biosynthesis requires the formation of lanosterol followed by the removal of three methyl groups: two from the C4 position and one from the C14 position (
Figure 1). Compared to the C14DM/CYP51 which is well characterized biochemically and genetically as the main target for azole drugs, the C4DM is less studied. In
L. donovani (the etiological agent for visceral leishmaniasis), we recently identified CYP5122A1 as a bona fide sterol C4-methyl oxidase capable of converting C4-methylated sterols into mono-hydroxy, aldehyde and carboxylate intermediates in vitro [
13]. Chromosomal CYP5122A1-null mutants cannot survive without a complementing CYP5122A1-episome, indicating that it is required for
L. donovani promastigotes and amastigotes. In this report, we investigated whether the essentiality of CYP5122A1 can be extended to
L. major which causes cutaneous leishmaniasis. Previous studies from us and others have revealed differences in the requirement of sterol and sphingolipid metabolism among
Leishmania species [
17,31–33]. For example, while C14DM/CYP51 can be removed from
L. major, it is indispensable in
L. donovani [
17,32]. In addition, the SL-degrading enzyme ISCL is required for the virulence of
L. major but not
L. amazonensis [31,33]. Thus, it is important to determine if different
Leishmania species have distinct requirements for C4DM.
Using the same plasmid segregation method as described for
L. donovani CYP5122A1 [
13,
20], we showed that
L. major 22A1‒ +pXNG4-
22A1 promastigotes retained a high level of episome even when facing strong negative selection from GCV treatment, whereas
22A1+/‒ +pXNG4-
22A1 lost the episome after 6 passages. When single clones isolated from
22A1‒ +pXNG4-
22A1 GCV treatment group were allowed to proliferate, they showed a mix of GFP-low and GFP-high cells, suggesting that GFP-low cells (lacking pXNG4-
22A1) could not perpetuate by themselves. Thus, like in
L. donovani, 22A1 is essential for
L. major promastigotes. Given the fact that
L. major C14DM/CYP51-null mutants are viable [
6], these findings suggest that removal of the two C4 methyl groups is more crucial than the removal of C14-methyl group. It is possible that the accumulation of C4-methyl intermediates is more harmful to
Leishmania membrane than the accumulation of C14-methyl intermediates [
22].
In the intracellular mammalian stage,
L. major amastigotes are less dependent on
de novo lipid synthesis and can salvage and remodel host lipids to support their activity [
6,
23,
24,
25]. To determine whether 22A1 is dispensable for amastigotes, BALB/c mice infected with
22A1‒ +pXNG4-
22A1 were treated with GCV or PBS for 14 days. As illustrated in
Figure 7, GCV treatment dramatically reduced the virulence and growth of
22A1‒ +pXNG4-
22A1 in mice, yet these parasites retained the pXNG4-
22A1 plasmid at 16 weeks post infection. These results are consistent with the notion that 22A1 is required during the amastigote stage. For mice infected by
22A1+/‒ +pXNG4-
22A1, GCV treatment caused an initial delay in footpad lesion formation which might arise from the elevated level of pXNG4-
22A1 in the early stage of infection. After 7 weeks, these parasites lost the episome and displayed normal levels of virulence and growth in mice, indicating that the plasmid is not essential for these parasites which contain one of chromosomal
22A1. These findings are largely consistent with the data from
L. donovani CYP5122A1. However, one distinction is that
L. major 22A1+/‒ (half knockout) shows WT-level of infectivity while
L. donovani half knockout has attenuated virulence [
13,
14]. Also,
L. major 22A1-overexpression and half knockout did not affect stress response or bulk sterol composition. These results suggest that
L. major parasites can tolerate the impact of 22A1 over- or under-expression better than
L. donovani. We did detect a slight accumulation of squalene in 22A1-overexpressing cells, suggesting that its overproduction can affect the early steps of sterol synthesis.
Overall, our results demonstrate that 22A1 is essential for
L. major amastigotes just like it is for
L. donovani amastigotes. This is notable because most sterols in intracellular amastigotes consist of host-derived cholesterol and only a small portion was synthesized
de novo [
6,
26]. These findings mark the potential of 22A1 as a pan-
Leishmania drug target. In culture promastigotes, the overexpression of 22A1 led to 2-3-fold resistance to several 22A1 and CYP51 inhibitors including DB766 and posaconazole (
Figure 8). Future studies will aim to identify more potent and selective inhibitors of 22A1.
Despite these findings, it is not clear whether CYP5122A1 is the sole sterol C4 methyl oxidase in
Leishmania. In
S. cerevisiae and animals, the two C-4 methyl groups are removed consecutively, and a single sterol C4 methyl oxidase catalyzes the initial oxygenation reactions [
8,
9,
10,
27]. The
S. cerevisiae ERG25-null mutants are ergosterol auxotroph and accumulate C4,4-dimethylzymosterol [
8]. Unlike
S. cerevisiae,
Aspergillus fumigatus has the two functional sterol C4 methyl oxidases (ERG25A and ERG25B) with ERG25A being the dominant enzyme. ERG25A mutants are viable but show susceptibility to hypoxia and ER stress [
28].
Candida albicans also has two sterol C4 methyl oxidase genes and both appear to be essential suggesting that they work in tandem for efficient activity [
29]. Similarly, plant cells possess two sterol C4 methyl oxidases each removing one C4-methyl group in a subsequent, but biochemically distinct mechanism [
11,
30,
31]. In
Arabidopsis thaliana, removal of each of the two C4-methyl groups requires a different sterol C4 methyl oxidase and mutants of the individual enzymes primarily accumulate their respective sterol intermediates [
28]. Finally, the methanotrophic gram-negative bacterium
Methylococcus capsulatus contains two C-4 demethylases that are distinct from the eukaryotic enzymes [
32]. Overall, these studies demonstrate diverse mechanisms involved in the C-4 demethylation step in different organisms.
In addition to
22A1, there are other genes in
Leishmania parasites encoding proteins homologous to
S. cerevisiae sterol C4 methyl oxidase: the previously characterized lathosterol oxidase (LSO, LmjF.23.1300) and a putative ERG25 (LmjF.36.2540). These orthologs show 29% and 26% identity to
S. cerevisiae sterol C4 methyl oxidase respectively.
L. major LSO-null mutants synthesize sterols without the C-5-C-6 double bond but show no defects in C4-demthylation [
33]. In this study, we characterized the putative ERG25 in
L. major. Like other sterol biosynthetic enzymes, ERG25 is largely localized in the ER.
L. major ERG25-null mutants were fully viable with growth rate and virulence similar to WT parasites (
Figure 4). These mutants synthesize ergosterol and 5-dehydroepisterol as their main sterols and their overall sterol composition is very similar to WT and add-back parasites (
Figure 5). These observations suggest that ERG25 is not involved in the oxidative removal of C4-methyl groups. Unlike other sterol biosynthetic enzymes, the deletion or overexpression of ERG25 had no impact on the cellular levels or localization of LPG and GP63 [
6,
18,
33].
While the precise function of ERG25 in
Leishmania awaits further investigation, it is important to note that ERG25 is not well conserved among trypanosomatids. According to Tritrypdb, ERG25 orthologs are present in
Leishmania spp.,
Angomonas deanei,
Bodo saltans,
Crithidia fasciculata,
Endotrypanum monterogeii,
Leptomonas spp.,
Paratrypanosoma confusum, and
Porcisia hertigi, but absent in
Trypanosoma spp. In contrast, CYP5122A1 orthologs are identified in all of the above organisms plus
Trypanosoma spp. including
T. brucei and
T. cruzi [
34]. Given the dispensable nature of
ERG25 in
L. major sterol synthesis, it is possible that
ERG25 has evolved to diverge from the sterol C4 methyl oxidase function and was lost in some trypanosomatid species. Meanwhile, CYP5122A1 was conserved by most trypanosomatids as the only or dominant sterol C4 methyl oxidase.
Figure 1.
Demethylation steps catalyzed by C4DM and C14DM in Leishmania sterol synthesis pathway. I: lanosterol; II: 4,14-dimethylzymosterol; III: 14-methylzymosterol; IV: 4,4-dimethylzymosterol; V: 4-methylzymosterol; VI: zymosterol.
Figure 1.
Demethylation steps catalyzed by C4DM and C14DM in Leishmania sterol synthesis pathway. I: lanosterol; II: 4,14-dimethylzymosterol; III: 14-methylzymosterol; IV: 4,4-dimethylzymosterol; V: 4-methylzymosterol; VI: zymosterol.
Figure 2.
Cellular localization of Cyp5122A1 and Erg25 in L. major. (A-B) Log phase promastigotes of WT +pXG-22A1-GFP (A) or WT + pXG-GFP-ERG25 (B) were labeled with an anti-T. brucei BiP antibody (ER marker) and processed for immunofluorescence microscopy. Merge images were the overlay of GFP and BiP. Overlap coefficients were determined from 30 cells each (average ± standard deviations). Scale bars: 10 μm. (C) Log phase promastigotes of WT, WT +pXG-22A1-GFP, WT + pXG-C14DM-GFP or WT + pXG-GFP-ERG25 were analyzed by western blot using antibodies against GFP (top) or α-tubulin.
Figure 2.
Cellular localization of Cyp5122A1 and Erg25 in L. major. (A-B) Log phase promastigotes of WT +pXG-22A1-GFP (A) or WT + pXG-GFP-ERG25 (B) were labeled with an anti-T. brucei BiP antibody (ER marker) and processed for immunofluorescence microscopy. Merge images were the overlay of GFP and BiP. Overlap coefficients were determined from 30 cells each (average ± standard deviations). Scale bars: 10 μm. (C) Log phase promastigotes of WT, WT +pXG-22A1-GFP, WT + pXG-C14DM-GFP or WT + pXG-GFP-ERG25 were analyzed by western blot using antibodies against GFP (top) or α-tubulin.
Figure 3.
Southern blots to confirm the deletions of chromosomal ERG25 and 22A1 alleles. Genomic DNA samples from WT,
ERG25+/-,
erg25¯ (three clones),
22A1+/‒, and
22A1¯+pXNG4-
22A1 (two clones) parasites were digested with restriction enzymes and hybridized with radiolabeled probes for the open reading frames (top) or upstream flanking regions (bottom) of
22A1 (A) or
Erg25 (B). Bands corresponding to
22A1, pXNG4-
22A1, Erg25, and antibiotic-resistant genes (
BSD/HYG/PAC) are indicated. Details of the Southern blots were included in
supplemental Figure. S3-S4.
Figure 3.
Southern blots to confirm the deletions of chromosomal ERG25 and 22A1 alleles. Genomic DNA samples from WT,
ERG25+/-,
erg25¯ (three clones),
22A1+/‒, and
22A1¯+pXNG4-
22A1 (two clones) parasites were digested with restriction enzymes and hybridized with radiolabeled probes for the open reading frames (top) or upstream flanking regions (bottom) of
22A1 (A) or
Erg25 (B). Bands corresponding to
22A1, pXNG4-
22A1, Erg25, and antibiotic-resistant genes (
BSD/HYG/PAC) are indicated. Details of the Southern blots were included in
supplemental Figure. S3-S4.
Figure 4.
ERG25 is not required for proliferation or virulence in L. major. (A-C) WT, erg25– and erg25–/+ERG25 promastigotes were inoculated in complete M199 media at 1.0 x 105 cells/ml and culture densities were determined daily (A). After cultures entered the stationary stage (when densities reached 2.0 x 107 cells/ml), percentages of round cells (B) and dead cells (C) were determined daily. (D-E) Day 3 stationary phase promastigotes were injected into the footpads of female BALB/c mice. The development of footpad lesions was monitored weekly (D) and parasite numbers in infected footpads were determined at 4- and 7-weeks (one mouse each) post infection via limiting dilution assay (E). Error bars represent standard deviations.
Figure 4.
ERG25 is not required for proliferation or virulence in L. major. (A-C) WT, erg25– and erg25–/+ERG25 promastigotes were inoculated in complete M199 media at 1.0 x 105 cells/ml and culture densities were determined daily (A). After cultures entered the stationary stage (when densities reached 2.0 x 107 cells/ml), percentages of round cells (B) and dead cells (C) were determined daily. (D-E) Day 3 stationary phase promastigotes were injected into the footpads of female BALB/c mice. The development of footpad lesions was monitored weekly (D) and parasite numbers in infected footpads were determined at 4- and 7-weeks (one mouse each) post infection via limiting dilution assay (E). Error bars represent standard deviations.
Figure 5.
ERG25 deletion does not affect sterol composition in L. major. (A-D) Partial GC chromatograms of free sterols from WT (A), erg25– (B), erg25–/+ERG25 (C) and erg25–/+GFP-ERG25 (D) promastigotes. The peaks represented ergosterol (1), 5-dehydroepisterol (2), cholesterol (3), cholesta-5,7,24-trienol (4) and cholesta-3,5-diene standard (std). Sterol compositions were summarized in E.
Figure 5.
ERG25 deletion does not affect sterol composition in L. major. (A-D) Partial GC chromatograms of free sterols from WT (A), erg25– (B), erg25–/+ERG25 (C) and erg25–/+GFP-ERG25 (D) promastigotes. The peaks represented ergosterol (1), 5-dehydroepisterol (2), cholesterol (3), cholesta-5,7,24-trienol (4) and cholesta-3,5-diene standard (std). Sterol compositions were summarized in E.
Figure 6.
22A1 is essential for L. major promastigotes in culture. (A) Promastigotes were continuously cultivated in the presence or absence of GCV or nourseothricin (SAT) and percentages of GFP-high cells were determined by flow cytometry for each passage. Error bars represent standard deviations from three replicates. (B–D) After 6 passages, single clones were isolated from 22A1¯+pXNG4-22A1 cells grown in the presence of GCV by serial dilution and allowed to proliferate. Flow cytometry analyses were performed on the original 22A1¯+pXNG4-22A1 GCV population (B) and two representative clones (C and D). In B–D, percentages of GFP-high cells were indicated as R3. (E) Promastigotes were inoculated at 1.0 x 105 cells/ml in the presence or absence of GCV or SAT and culture densities were determined daily (***: p < 0.001).
Figure 6.
22A1 is essential for L. major promastigotes in culture. (A) Promastigotes were continuously cultivated in the presence or absence of GCV or nourseothricin (SAT) and percentages of GFP-high cells were determined by flow cytometry for each passage. Error bars represent standard deviations from three replicates. (B–D) After 6 passages, single clones were isolated from 22A1¯+pXNG4-22A1 cells grown in the presence of GCV by serial dilution and allowed to proliferate. Flow cytometry analyses were performed on the original 22A1¯+pXNG4-22A1 GCV population (B) and two representative clones (C and D). In B–D, percentages of GFP-high cells were indicated as R3. (E) Promastigotes were inoculated at 1.0 x 105 cells/ml in the presence or absence of GCV or SAT and culture densities were determined daily (***: p < 0.001).
Figure 7.
22A1 is indispensable for L. major amastigotes. BALB/c mice were infected in the footpads with stationary phase promastigotes and treated with either GCV or PBS daily for 14 days. (A) Footpad lesions were recorded weekly. (B) Average parasite numbers per infected footpad were determined at the indicated times post infection by qPCR. (C) Average pXNG4-22A1 copy numbers in amastigotes were determined at the indicated times post infection by qPCR. ND: not determined. Error bars represent standard deviations from 4-5 mice (A) or 2-3 mice (B and C). **: p < 0.01.
Figure 7.
22A1 is indispensable for L. major amastigotes. BALB/c mice were infected in the footpads with stationary phase promastigotes and treated with either GCV or PBS daily for 14 days. (A) Footpad lesions were recorded weekly. (B) Average parasite numbers per infected footpad were determined at the indicated times post infection by qPCR. (C) Average pXNG4-22A1 copy numbers in amastigotes were determined at the indicated times post infection by qPCR. ND: not determined. Error bars represent standard deviations from 4-5 mice (A) or 2-3 mice (B and C). **: p < 0.01.
Figure 8.
Effects of 22A1 overexpression on promastigotes’ sensitivity to sterol synthesis inhibitors. Log phase promastigotes were inoculated in various concentrations of ketoconazole (A), posaconazole (B) or DB766 (C). Cell densities were determined after 48 hours, and percentages of growth were calculated using cells grown in the absence of inhibitors as controls. Error bars represent standard deviations from three experiments.
Figure 8.
Effects of 22A1 overexpression on promastigotes’ sensitivity to sterol synthesis inhibitors. Log phase promastigotes were inoculated in various concentrations of ketoconazole (A), posaconazole (B) or DB766 (C). Cell densities were determined after 48 hours, and percentages of growth were calculated using cells grown in the absence of inhibitors as controls. Error bars represent standard deviations from three experiments.
Figure 9.
Expression levels of 22A1 and C14DM in L. major promastigotes. (A) Whole cell lysates from log phase promastigotes were denatured at 65 ℃ and analyzed by western blot using antibodies against C14DM (top), 22A1 (middle), or α-tubulin (bottom). (B) The relative abundances of 22A1 and C14DM were determined after three experiments using α-tubulin as the loading control. Error bars represent standard deviations. ***p < 0.001.
Figure 9.
Expression levels of 22A1 and C14DM in L. major promastigotes. (A) Whole cell lysates from log phase promastigotes were denatured at 65 ℃ and analyzed by western blot using antibodies against C14DM (top), 22A1 (middle), or α-tubulin (bottom). (B) The relative abundances of 22A1 and C14DM were determined after three experiments using α-tubulin as the loading control. Error bars represent standard deviations. ***p < 0.001.
Table 1.
EC25 and EC50 values for posaconazole and DB766.
Table 1.
EC25 and EC50 values for posaconazole and DB766.
Inhibitor |
Effective concentration (EC) |
WT |
22A1+/‒ |
22A1+/‒ +pXNG4-22A1 |
22A1‒ +pXNG4-22A1
|
Posaconazole |
EC25 (µM ± SD) |
0.29 ± 0.089 |
0.21 ± 0.041 |
0.48 ± 0.15 |
0.47 ± 0.14 |
EC50 (µM ± SD) |
1.0 ± 0.096 |
0.74 ± 0.18 |
3.0 ± 1.0* |
1.7 ± 0.38* |
DB766 |
EC25 (nM ± SD) |
0.90 ± 0.12 |
0.96 ± 0.16 |
1.6 ± 0.33 |
1.9 ± 0.094 |
EC50 (nM ± SD) |
1.9 ± 0.31 |
2.04 ± 0.41 |
3.8 ± 0.68** |
3.5 ± 0.19** |