Signatures of natural selection in tree topology shape of serially sampled viral phylogenies.

Tree shape metrics can be computed fast for trees of any size, which makes them promising alternatives to intensive statistical methods and parameter-rich evolutionary models in the era of massive data availability. Previous studies have demonstrated their effectiveness in unveiling important parameters in viral evolutionary dynamics, although the impact of natural selection on the shape of tree topologies has not been thoroughly investigated. We carried out a forward-time and individual-based simulation to investigate whether tree shape metrics of several kinds could predict the selection regime employed to generate the data. To examine the impact of the genetic diversity of the founder viral population, simulations were run under two opposing starting configurations of the genetic diversity of the infecting viral population. We found that four evolutionary regimes, namely, negative, positive, and frequency-dependent selection, as well as neutral evolution, were successfully distinguished by tree topology shape metrics. Two metrics from the Laplacian spectral density profile (principal eigenvalue and peakedness) and the number of cherries were the most informative for indicating selection type. The genetic diversity of the founder population had an impact on differentiating evolutionary scenarios. Tree imbalance, which has been frequently associated with the action of natural selection on intrahost viral diversity, was also characteristic of neutrally evolving serially sampled data. Metrics calculated from empirical analysis of HIV datasets indicated that most tree topologies exhibited shapes closer to the frequency-dependent selection or neutral evolution regimes.

Assessing the relative performance of fast molecular dating methods for phylogenomic data.

Advances in genome sequencing techniques produced a significant growth of phylogenomic datasets. This massive amount of data represents a computational challenge for molecular dating with Bayesian approaches. Rapid molecular dating methods have been proposed over the last few decades to overcome these issues. However, a comparative evaluation of their relative performance on empirical data sets is lacking. We analyzed 23 empirical phylogenomic datasets to investigate the performance of two commonly employed fast dating methodologies: penalized likelihood (PL), implemented in treePL, and the relative rate framework (RRF), implemented in RelTime. They were compared to Bayesian analyses using the closest possible substitution models and calibration settings. We found that RRF was computationally faster and generally provided node age estimates statistically equivalent to Bayesian divergence times. PL time estimates consistently exhibited low levels of uncertainty. Overall, to approximate Bayesian approaches, RelTime is an efficient method with significantly lower computational demand, being more than 100 times faster than treePL. Thus, to alleviate the computational burden of Bayesian divergence time inference in the era of massive genomic data, molecular dating can be facilitated using the RRF, allowing evolutionary hypotheses to be tested more quickly and efficiently.

The performance of outgroup-free rooting under evolutionary radiations.

Tree rooting implies a temporal dimension to phylogenies. Only after defining the position of the root node is that the ancestral-descendant relationship between branches can be fully deduced. Rooting has been usually carried out by employing evolutionarily close outgroup lineages, which is a drawback when these lineages are unavailable or unknown. Alternatively, outgroup-free rooting methods were proposed, which rely on the constancy of evolutionary rates to varying degrees. In this work we analyzed the performance of two of these methods, the midpoint rooting (MPR) and the minimal ancestor deviation (MAD), in rooting topologies evolved under challenging scenarios of fast evolutionary radiations derived from empirical data, characterized by short internal branches near the crown node. Considering all branch length combinations investigated, both methods exhibited average success rates below 50%, although MAD slightly outperformed MPR. Moreover, tree balance significantly impacted the relative performance of the methods. We found that, in four-taxa unrooted trees, the outcome of whether both methodologies will correctly root the tree can be roughly predicted by two simple dimensionless metrics: the coefficient of variation of the external branch lengths, and the ratio between the internal branch length to the total sum of branch lengths, which were employed to devise a general linear model that allowed calculating the probability of correct placing the root node for any four-taxa tree. We predicted that the performance of both outgroup-free rooting methods on loci representing the placental mammal radiation ranged between 50% and 75%.

Challenges in estimating virus divergence times in short epidemic timescales with special reference to the evolution of SARS-CoV-2 pandemic.

The estimation of evolutionary parameters provides essential information for designing public health policies. In short time intervals, however, nucleotide substitutions are ineffective to record all complexities of virus population dynamics. In this sense, the current SARS-CoV-2 pandemic poses a challenge for evolutionary analysis. We used computer simulation to evolve populations in scenarios of varying temporal intervals to evaluate the impact of the age of an epidemic on estimates of time and geography. Before estimating virus timescales, the shape of tree topologies can be used as a proxy to assess the effectiveness of the virus phylogeny in providing accurate estimates of evolutionary parameters. In short timescales, estimates have larger uncertainty. We compared the predictions from simulations with empirical data. The tree shape of SARS-CoV-2 was closer to shorter timescales scenarios, which yielded parametric estimates with larger uncertainty, suggesting that estimates from these datasets should be evaluated cautiously. To increase the accuracy of the estimates of virus transmission times between populations, the uncertainties associated with the age estimates of both the crown and stem nodes should be communicated. We place the age of the common ancestor of the current SARS-CoV-2 pandemic in late September 2019, corroborating an earlier emergence of the virus.

Phylogenetic analysis of the ATP-binding cassette proteins suggests a new ABC protein subfamily J in Aedes aegypti (Diptera: Culicidae).

BACKGROUND: We performed an in-depth analysis of the ABC gene family in Aedes aegypti (Diptera: Culicidae), which is an important vector species of arthropod-borne viral infections such as chikungunya, dengue, and Zika. Despite its importance, previous studies of the Arthropod ABC family have not focused on this species. Reports of insecticide resistance among pests and vectors indicate that some of these ATP-dependent efflux pumps are involved in compound traffic and multidrug resistance phenotypes. RESULTS: We identified 53 classic complete ABC proteins annotated in the A. aegypti genome. A phylogenetic analysis of Aedes aegypti ABC proteins was carried out to assign the novel proteins to the ABC subfamilies. We also determined 9 full-length sequences of DNA repair (MutS, RAD50) and structural maintenance of chromosome (SMC) proteins that contain the ABC signature. CONCLUSIONS: After inclusion of the putative ABC proteins into the evolutionary tree of the gene family, we classified A. aegypti ABC proteins into the established subfamilies (A to H), but the phylogenetic positioning of MutS, RAD50 and SMC proteins among ABC subfamilies-as well as the highly supported grouping of RAD50 and SMC-prompted us to name a new J subfamily of A. aegypti ABC proteins.

Survey for positively selected coding regions in the genome of the hematophagous tsetse fly Glossina morsitans identifies candidate genes associated with feeding habits and embryonic development.

Tsetse flies are responsible for the transmission of Trypanossoma sp. to vertebrate animals in Africa causing huge health issues and economic loss. The availability of the genome sequence of Glossina morsitans enabled the discovery of several genes related to medically important phenotypes and novel physiological features. However, a genome-wide scan for coding regions that underwent positive selection is still missing, which is surprising given the evolution of traits associated with the hematophagy in this lineage. In this study, we employed an experimental design that controlled for the rate of false positives and we performed a scan of 3,318 G. morsitans genes. We found 145 genes with significant historical signal of positive selection. These genes were categorized into 18 functional classes after careful manual annotation. Based on their attributed functions, we identified candidate genes related with feeding habits and embryonic development. When our results were contrasted with gene expression data, we confirmed that most genes that underwent adaptive molecular evolution were frequently expressed in organs associated with key physiological evolutionary innovations in the G. morsitans lineage, namely, the salivary gland, the midgut, fat body tissue, and in the spermatophore.

A multigene timescale and diversification dynamics of Ciliophora evolution.

Ciliophora is one of the most diverse lineages of unicellular eukaryotes. Nevertheless, a robust timescale including all main lineages and employing properly identified ciliate fossils as primary calibrations is lacking. Here, we inferred a time-calibrated multigene phylogeny of Ciliophora evolution, and we used this timetree to investigate the rates and patterns of lineage diversification through time. We implemented a two-step analytical approach that favored both gene and taxon sampling, reducing the uncertainty of time estimates and yielding narrower credibility intervals on the ribosomal-derived chronogram. We estimate the origin of Ciliophora at 1143 Ma, which is substantially younger than previously proposed ages, and the huge diversity explosion occurred during the Paleozoic. Among the current groups recognized as classes, Spirotrichea diverged earlier, its origin was dated at ca. 850 Ma, and Protocruziea was the younger class, with crown age estimated at 56 Ma. Macroevolutionary analysis detected a significant rate shift in diversification dynamics in the spirotrichean clade Hypotrichia + Oligotrichia + Choreotrichia, which had accelerated speciation rate ca. 570 Ma, during the Ediacaran-Cambrian transition. For all crown lineages investigated, speciation rates declined through time, whereas extinction rates remained low and relatively constant throughout the evolutionary history of ciliates.

The range of sampling times affects Zika virus evolutionary rates and divergence times.

The rate of evolution of viral genomes is a fundamental parameter for understanding the origin and spread of epidemics. For instance, molecular dating is one of the many practical outcomes of evolutionary rate estimation. In this sense, the rate of evolution of ZIKV merits attention, because it has been shown to be higher than the average rate reported for other flaviviruses. It has been hypothesized that the higher rate of ZIKV evolution is due to a bias related to the analysis of sequences collected within a short time range, which would increase the chance of sampling slightly deleterious nucleotide polymorphisms. To investigate this hypothesis, we assembled datasets with different ranges of sampling times and also decomposed the ZIKV evolutionary rate into synonymous and non-synonymous rates. Our results demonstrated that the rate of ZIKV evolution is time dependent and that the observed increase in short-term rates is largely accounted for by a higher non-synonymous rate, suggesting the presence of slightly deleterious variants not yet eliminated by purifying selection. On the other hand, we show that synonymous rates were less impacted by the range of sampling times, generating timescales with reduced uncertainty. We conclude that, for inferring the ZIKV timescale and reconstructing the history of epidemics, synonymous changes are the most appropriate substitution type to be examined. We were able to obtain ZIKV divergence times that were time independent and exhibited greater precision than previous estimates. This observation should also hold for other serially sampled fast-evolving pathogens with evidence of time dependence of evolutionary rates.

Large ancestral effective population size explains the difficult phylogenetic placement of owl monkeys.

The phylogenetic position of owl monkeys, grouped in the genus Aotus, has been a controversial issue for understanding Neotropical primate evolution. Explanations of the difficult phylogenetic assignment of owl monkeys have been elusive, frequently relying on insufficient data (stochastic error) or scenarios of rapid speciation (adaptive radiation) events. Using a coalescent-based approach, we explored the population-level mechanisms likely explaining these topological discrepancies. We examined the topological variance of 2,192 orthologous genes shared between representatives of the three major Cebidae lineages and the outgroup. By employing a methodological framework that allows for reticulated tree topologies, our analysis explicitly tested for non-dichotomous evolutionary processes impacting the finding of the position of owl monkeys in the cebid phylogeny. Our findings indicated that Aotus is a sister lineage of the callitrichines. Most gene trees (>50%) failed to recover the species tree topology, although the distribution of gene trees mismatching the true species topology followed the standard expectation of the multispecies coalescent without reticulation. We showed that the large effective population size of the common ancestor of Aotus and callitrichines was the most likely factor responsible for generating phylogenetic uncertainty. On the other hand, fast speciation scenarios or introgression played minor roles. We propose that the difficult phylogenetic placement of Aotus is explained by population-level processes associated with the large ancestral effective size. These results shed light on the biogeography of the early cebid diversification in the Miocene, highlighting the relevance of evaluating phylogenetic relationships employing population-aware approaches.

Multispecies coalescent analysis confirms standing phylogenetic instability in Hexapoda.

The multispecies coalescent (MSC) has been increasingly used in phylogenomic analyses due to the accommodation of gene tree topological heterogeneity by taking into account population-level processes, such as incomplete lineage sorting. In this sense, the phylogeny of insect species, which are characterized by their large effective population sizes, is suitable for a coalescent-based analysis. Furthermore, studies so far recovered short internal branches at early divergences of the insect tree of life, indicating fast evolutionary radiations that increase the probability of incomplete lineage sorting in deep time. Here, we investigated the performance of the MSC for a phylogenomic data set of hexapods compiled by Misof et al. (2014, Science 346:763). Our analysis recovered the monophyly of most insect orders, and major phylogenetic relationships were in agreement with current insect systematics. We identified, however, some evolutionary associations that were consistently problematic. Most noticeable, Hexapod monophyly was disrupted by the sister group relationship between the remiped crustacean and Insecta. Additionally, the interordinal relationships within Polyneoptera and Neuropteroidea were found to be phylogenetically unstable. We show that these controversial phylogenetic arrangements were also poorly supported by previous analyses, and therefore, we evaluated their robustness to stochastic errors from sampling sites and terminals, confirming standing problems in hexapod phylogeny in the genomics age.

Comparative evaluation of maximum parsimony and Bayesian phylogenetic reconstruction using empirical morphological data.

The use of discrete morphological data in Bayesian phylogenetics has increased significantly over the last years with the proposal of total evidence analysis and the treatment of fossils as terminal taxa in Bayesian molecular dating. Both approaches rely on the assumption that probabilistic Markov models reasonably accommodate all the complexity of morphological evolution of discrete traits. The performance of such morphological models used in Bayesian phylogenetics has been thoroughly investigated, but conclusions so far were based mostly on simulated data. In this study, we have surveyed MorphoBank and obtained a large number of morphological matrices to evaluate Bayesian phylogenetic inference (BI) under Lewis' Mk model in comparison with the maximum parsimony (MP) algorithm. We found that trees estimated by both methods frequently differed and that BI generated a larger amount of polytomic tree topologies. The number of trees contained in the 95% Bayesian credibility interval was significantly greater than the number of equally parsimonious trees. We also investigated which factors mostly influenced the topological difference between maximum parsimony and Bayesian tree topologies and found that the number of terminals in morphological matrices was the variable with the highest association with the topological distance between trees inferred by BI and MP. Surprisingly, we show that differences between both approaches were not influenced by increasing sample size. Our results, which were based on a large set of empirical matrices, corroborate recent findings that BI is less precise than MP.

Discovery and evolution of novel hemerythrin genes in annelid worms.

BACKGROUND: Despite extensive study on hemoglobins and hemocyanins, little is known about hemerythrin (Hr) evolutionary history. Four subgroups of Hrs have been documented, including: circulating Hr (cHr), myohemerythrin (myoHr), ovohemerythrin (ovoHr), and neurohemerythrin (nHr). Annelids have the greatest diversity of oxygen carrying proteins among animals and are the only phylum in which all Hr subgroups have been documented. To examine Hr diversity in annelids and to further understand evolution of Hrs, we employed approaches to survey annelid transcriptomes in silico. RESULTS: Sequences of 214 putative Hr genes were identified from 44 annelid species in 40 different families and Bayesian inference revealed two major clades with strong statistical support. Notably, the topology of the Hr gene tree did not mirror the phylogeny of Annelida as presently understood, and we found evidence of extensive Hr gene duplication and loss in annelids. Gene tree topology supported monophyly of cHrs and a myoHr clade that included nHrs sequences, indicating these designations are functional rather than evolutionary. CONCLUSIONS: The presence of several cHrs in early branching taxa suggests that a variety of Hrs were present in the common ancestor of extant annelids. Although our analysis was limited to expressed-coding regions, our findings demonstrate a greater diversity of Hrs among annelids than previously reported.

Multilocus phylogeny and statistical biogeography clarify the evolutionary history of major lineages of turtles.

Despite their complex evolutionary history and the rich fossil record, the higher level phylogeny and historical biogeography of living turtles have not been investigated in a comprehensive and statistical framework. To tackle these issues, we assembled a large molecular dataset, maximizing both taxonomic and gene sampling. As different models provide alternative biogeographical scenarios, we have explicitly tested such hypotheses in order to reconstruct a robust biogeographical history of Testudines. We scanned publicly available databases for nucleotide sequences and composed a dataset comprising 13 loci for 294 living species of Testudines, which accounts for all living genera and 85% of their extant species diversity. Phylogenetic relationships and species divergence times were estimated using a thorough evaluation of fossil information as calibration priors. We then carried out the analysis of historical biogeography of Testudines in a fully statistical framework. Our study recovered the first large-scale phylogeny of turtles with well-supported relationships following the topology proposed by phylogenomic works. Our dating result consistently indicated that the origin of the main clades, Pleurodira and Cryptodira, occurred in the early Jurassic. The phylogenetic and historical biogeographical inferences permitted us to clarify how geological events affected the evolutionary dynamics of crown turtles. For instance, our analyses support the hypothesis that the breakup of Pangaea would have driven the divergence between the cryptodiran and pleurodiran lineages. The reticulated pattern in the ancestral distribution of the cryptodiran lineage suggests a complex biogeographic history for the clade, which was supposedly related to the complex paleogeographic history of Laurasia. On the other hand, the biogeographical history of Pleurodira indicated a tight correlation with the paleogeography of the Gondwanan landmasses.

Evolutionary analysis of Chironius snakes unveils cryptic diversity and provides clues to diversification in the Neotropics.

Recent hypotheses to explain tropical diversity involves the Neogene and Quaternary geoclimatic dynamics, but the absence of unambiguous data permitting the choice between alternative hypotheses makes a general theory for the origin of tropical biodiversity far to be achieved. The occurrence of Chironius snakes in well-defined biogeographical regions led us to adopt Chironius as a model to unveil patterns of vertebrate diversification in the Neotropics. Here, we used molecular markers and records on geographic distribution to investigate Chironius evolution and, subsequently, providing hints on diversification in the Neotropics. To avoid analyzing nominal species that do not constitute exclusive evolutionary lineages, we firstly conducted a species delimitation study prior to carrying out the species distribution modeling analysis. We generated 161 sequences of 12S, 16S, c-mos and rag2 for 15 species and 50 specimens, and included additional data from GenBank yielding a matrix of 137 terminals, and performed the following evolutionary analyses: inference of a concatenated gene tree, estimation of gene divergence times, inference of the coalescent-based phylogeny of Chironius, estimation of effective population sizes and modeling potential distribution of species across the last millennia. We tested for species boundaries within Chironius by implementing a coalescent-based Bayesian species delimitation approach. Our analyses supported the monophyly of Chironius, although our findings underscored cryptic candidate species in C. flavolineatus and C. exoletus. The inferred timetree suggested that Chironius snakes have evolved in the early Miocene (ca. 20.2Mya) and began to diversify from the late Miocene to the early Pliocene, values that are much older than previously reported. Following genetic divergence of virtually all extant Chironius species investigated, the effective sizes of the populations have expanded when compared to their MRCAs. The evolutionary and SDM data from C. brazili and C. diamantina provided additional evidence of the origin of species in the Neotropics. We argue that temperature, instead of altitude, has been the major driving factor in the evolution of both species, and thus we present a case for the consequences of global warming.

Expanded phylogenetic analyses of the class Heterotrichea (Ciliophora, Postciliodesmatophora) using five molecular markers and morphological data.

Most studies of the molecular evolution of Heterotrichea have been based solely on the 18S-rDNA gene, which were inconsistent with morphological classification. Because of the limitations of single locus phylogenies and the recurring problem of lack of resolution of deeper nodes found in previous studies, we present hypotheses of the evolution of internal groups of the class Heterotrichea based on multi-loci analyses (18S-rDNA, 28S-rDNA, ITS1-5.8S-ITS2 region, COI and alpha-tubulin) and morphological data. Phylogenetic trees from protein coding gene data are presented for Heterotrichea for the first time. Phylogenetic analyses included Bayesian inference, maximum likelihood, maximum parsimony methods, and optimal trees were statistically compared to alternative topologies from the literature. Additionally, the Bayesian concordance approach (BCA algorithm) was used to assess the concordance factor between topologies obtained from isolated analyses. Because different loci may evolve at different rates, resulting in different gene topologies, we also estimated a species tree for Heterotrichea using the STAR coalescence-based method. The results show that: (1) single gene trees are inconsistent regarding the position of some heterotrichean families; (2) the concatenation of all data in a total-evidence tree improved the resolution of deep nodes among the heterotrichean families and genera; (3) the coalescent-based species tree is consistent with phylogenies based on the 18S-rDNA gene and shows Spirostomidae as the stem group of Heterotrichea; (4) however, the total-evidence tree suggests that the large Heterotrichea cluster is divided into nine lineages in which Peritromidae diverges at the base of the Heterotrichea tree.

Boronated tartrolon antibiotic produced by symbiotic cellulose-degrading bacteria in shipworm gills.

Shipworms are marine wood-boring bivalve mollusks (family Teredinidae) that harbor a community of closely related Gammaproteobacteria as intracellular endosymbionts in their gills. These symbionts have been proposed to assist the shipworm host in cellulose digestion and have been shown to play a role in nitrogen fixation. The genome of one strain of Teredinibacter turnerae, the first shipworm symbiont to be cultivated, was sequenced, revealing potential as a rich source of polyketides and nonribosomal peptides. Bioassay-guided fractionation led to the isolation and identification of two macrodioloide polyketides belonging to the tartrolon class. Both compounds were found to possess antibacterial properties, and the major compound was found to inhibit other shipworm symbiont strains and various pathogenic bacteria. The gene cluster responsible for the synthesis of these compounds was identified and characterized, and the ketosynthase domains were analyzed phylogenetically. Reverse-transcription PCR in addition to liquid chromatography and high-resolution mass spectrometry and tandem mass spectrometry revealed the transcription of these genes and the presence of the compounds in the shipworm, suggesting that the gene cluster is expressed in vivo and that the compounds may fulfill a specific function for the shipworm host. This study reports tartrolon polyketides from a shipworm symbiont and unveils the biosynthetic gene cluster of a member of this class of compounds, which might reveal the mechanism by which these bioactive metabolites are biosynthesized.

The effective population sizes of the anthropoid ancestors of the human-chimpanzee lineage provide insights on the historical biogeography of the great apes.

The recent development of methods that apply coalescent theory to phylogenetic problems has enabled the study of the population-level phenomena that drove the diversification of anthropoid primates. Effective population size, Ne, is one of the main parameters that constitute the theoretical underpinning of these new analytical approaches. For this reason, the ancestral N(e) of selected primate lineages has been thoroughly investigated. However, for some of these lineages, the estimates of ancestral N(e) reported in several studies present significant variation. This is the case for the common ancestor of humans and chimpanzees. Moreover, several ancestral anthropoid lineages have been ignored in the studies conducted so far. Because N(e) is fundamental to understand historic species demography, it is a crucial component of a complete description of the historical scenario of primate evolution. It also provides information that is helpful for differentiating between competing biogeographical hypotheses. In this study, the effective population sizes of the anthropoid ancestors of the human-chimp lineage are inferred using data sets of coding and noncoding sequences. A general pattern of a serial decline of population sizes is found between the ancestral lineage of Anthropoidea and that of Homo and Pan. When the theoretical distribution of gene trees was derived from the parametric estimates obtained, it closely corresponded to the empirical frequency of inferred gene trees along the genome. The most abrupt decrease of N(e) was found between the ancestors of all great apes and those of the African great apes alone. This suggests the occurrence of a genetic bottleneck during the evolution of Homininae, which corroborates the origin of African apes from a Eurasian ancestor.

The influence of taxon sampling on Bayesian divergence time inference under scenarios of rate heterogeneity among lineages.

Although taxon sampling is commonly considered an important issue in phylogenetic inference, it is rarely considered in the Bayesian estimation of divergence times. In fact, the studies conducted to date have presented ambiguous results, and the relevance of taxon sampling for molecular dating remains unclear. In this study, we developed a series of simulations that, after six hundred Bayesian molecular dating analyses, allowed us to evaluate the impact of taxon sampling on chronological estimates under three scenarios of among-lineage rate heterogeneity. The first scenario allowed us to examine the influence of the number of terminals on the age estimates based on a strict molecular clock. The second scenario imposed an extreme example of lineage specific rate variation, and the third scenario permitted extensive rate variation distributed along the branches. We also analyzed empirical data on selected mitochondrial genomes of mammals. Our results showed that in the strict molecular-clock scenario (Case I), taxon sampling had a minor impact on the accuracy of the time estimates, although the precision of the estimates was greater with an increased number of terminals. The effect was similar in the scenario (Case III) based on rate variation distributed among the branches. Only under intensive rate variation among lineages (Case II) taxon sampling did result in biased estimates. The results of an empirical analysis corroborated the simulation findings. We demonstrate that taxonomic sampling affected divergence time inference but that its impact was significant if the rates deviated from those derived for the strict molecular clock. Increased taxon sampling improved the precision and accuracy of the divergence time estimates, but the impact on precision is more relevant. On average, biased estimates were obtained only if lineage rate variation was pronounced.

Redescription and Phylogenetic Position of Condylostoma arenarium Spiegel, 1926 (Ciliophora, Heterotrichea) from Guanabara Bay, Brazil.

Details on Condylostoma arenarium infraciliature have not been described; therefore, it is considered a poorly known species. The lack of detailed description on C. arenarium morphology caused several misidentifications that have accumulated in the literature. In this study, we present the first complete description of C. arenarium infraciliature based on protargol-impregnated organisms and scanning electron microscopy. We also have inferred the phylogenetic position of this species based on 18S rRNA sequences. The main characteristics of C. arenarium population from Guanabara Bay are as follows: in vivo elongated body shape with 350-600 µm length × 70-220 µm width, they are highly contractile when subjected to disturbances, green-yellowish cortical granules are present, contractile vacuoles absent, V-shaped peristome comprises approximately 1/5 of the total length, adoral zone with 83-145 membranelles, 1-2 small frontal cirrus observed only in impregnated specimens, 10-15 fiber-like stripes arranged transversely on the inner wall of the oral cavity, 30-45 somatic kineties, moniliform macronucleus with 15-20 nodules. Some observations on morphogenesis of C. arenarium were also included. In phylogenetic analyses, C. arenarium clustered with Condylostoma sp. within a clade composed of three C. curva sequences with high support values.