Maximum-likelihood divergence date estimates based on rRNA gene sequences suggest two scenarios of Trypanosoma cruzi intraspecific evolution

Abstract

The phylogenetic relationships of Trypanosoma cruzi strains were inferred using maximum-likelihood from complete 18S rDNA sequences and D7-24Sα rDNA regions from 20 representative strains of T. cruzi. For this we sequenced the 18S rDNA of 14 strains and the D7-24Sα rDNA of four strains and aligned them to previously published sequences. Phylogenies inferred from these data sets identified four groups, named Riboclades 1, 2, 3, and 4, and a basal dichotomy that separated Riboclade 1 from Riboclades 2, 3, and 4. Substitution models and other parameters were optimized by hierarchical likelihood tests, and our analysis of the 18S rDNA molecular clock by the likelihood ratio test suggests that a taxa subset encompassing all 2,150 positions in the alignment supports rate constancy among lineages. The present analysis supports the notion that divergence dates of T. cruzi Riboclades can be estimated from 18S rDNA sequences and therefore, we present alternative evolutionary scenarios based on two different views of T. cruzi intraspecific divergence. The first assumes a faster evolutionary rate, which suggests that the divergence between T. cruzi I and II and the extant strains occurred in the Tertiary period (37-18 MYA). The other, which supports the hypothesis that the divergence between T. cruzi I and II occurred in the Cretaceous period (144-65 MYA) and the divergence of the extant strains occurred in the Tertiary period of the Cenozoic era (65-1.8 MYA), is consistent with our previously proposed hypothesis of divergence by geographical isolation and mammalian host coevolution.