Background. Sterols such as cholesterol, are important components of cellular membranes. But unlike mammalian cells, the main sterols found in the membranes of trypanosomes and fungi are ergosterol, and other 24-methyl sterols, which are required for growth and viability. In spite of this strict requirement, this group of organisms have evolved different strategies to produce and/or obtain sterols. Trypanosoma cruzi is the causative agent of Chagas Disease. In this parasite, one of the few validated targets for chemotherapeutic intervention is the sterol biosynthesis pathway. In this work we present a study of the genetic diversity observed in genes of the isoprenoid and sterol biosynthesis pathways in T. cruzi, and a comparative analysis of the diversity found in other trypanosomatids.
Methodology/Principal Findings. Using a number of bioinformatic strategies, we first completed a number of holes in the pathway by identifying the sequences of genes that were missing and/or were truncated in the draft T. cruzi genome. Based on this analysis we identified a non-orthologous homolog of the yeast ERG25 gene (sterol methyl oxidase, SMO) and propose that the orthologs of ERG25 have been lost in trypanosomes (but not in leishmanias). Next, starting from a set of 16 T. cruzi strains representative of six major evolutionary lineages, we have amplified and sequenced ~ 24Kbp from 18 genes of the pathway, and identified a total of 975 SNPs or fixed differences, of which 28% represent nonsynonymous changes. We observed different patterns of accumulation of nucleotide changes for different genes of the pathway, from genes with a density of substitutions ranging from those close to the average (~2.5/100 bp) to some showing a high number of changes (11.4/100 bp, for a putative lathosterol oxidase gene). The majority of genes are under apparent purifying selection. However, two genes (TcPMK, TcSMO-like) have a ratio of nonsynonymous to synonymous changes that is close to neutrality. None of the nonsynonymous changes identified affect a catalytic or a ligand binding site residue. However, after mapping these changes on top of available structural data, we identified a number of changes that are in the close vicinity (7 Angstrom) of key residues, and that could therefore be functionally important. A comparative analysis of the corresponding T. brucei and Leishmania genes, obtained from available complete genomes highlights a high degree of conservation of the pathway, but with differences in the genes that are under apparent purifying selection in each case.
Conclusions/Significance. We have identified a number of genes of the sterol biosynthesis pathway that were missing from the T. cruzi genome assembly. Also, we have identified unequal apparent selection acting on these genes, which may provide essential information for the future of drug development studies focused on this pathway.