Renata C. Ferreira

The Genome Sequence of Leishmania (Leishmania) amazonensis: Functional Annotation and Extended Analysis of Gene Models

We present the sequencing and annotation of the Leishmania (Leishmania) amazonensis genome, an etiological agent of human cutaneous leishmaniasis in the Amazon region of Brazil. L. (L.) amazonensis shares features with Leishmania (L.) mexicana but also exhibits unique characteristics regarding geographical distribution and clinical manifestations of cutaneous lesions (e.g. borderline disseminated cutaneous leishmaniasis). Predicted genes were scored for orthologous gene families and conserved domains in comparison with other human pathogenic Leishmania spp. Carboxypeptidase, aminotransferase, and 3′-nucleotidase genes and ATPase, thioredoxin, and chaperone-related domains were represented more abundantly in L. (L.) amazonensis and L. (L.) mexicana species. Phylogenetic analysis revealed that these two species share groups of amastin surface proteins unique to the genus that could be related to specific features of disease outcomes and host cell interactions. Additionally, we describe a hypothetical hybrid interactome of potentially secreted L. (L.) amazonensis proteins and host proteins under the assumption that parasite factors mimic their mammalian counterparts. The model predicts an interaction between an L. (L.) amazonensis heat-shock protein and mammalian Toll-like receptor 9, which is implicated in important immune responses such as cytokine and nitric oxide production. The analysis presented here represents valuable information for future studies of leishmaniasis pathogenicity and treatment.

Molecular Identification, Antifungal Susceptibility Profile, and Biofilm Formation of Clinical and Environmental Rhodotorula Species Isolates

ABSTRACT Rhodotorula species are emergent fungal pathogens capable of causing invasive infections, primarily fungemia. They are particularly problematic in immunosuppressed patients when using a central venous catheter. In this study, we evaluated the species distribution of 51 clinical and 8 environmental Rhodotorula species isolates using the ID32C system and internal transcribed spacer (ITS) sequencing. Antifungal susceptibility testing and biofilm formation capability using a crystal violet staining assay were performed. Using ITS sequencing as the gold standard, the clinical isolates were identified as follows: 44 R. mucilaginosa isolates, 2 R. glutinis isolates, 2 R. minuta isolates, 2 R. dairenensis isolates, and 1 Rhodosporidium fluviale isolate. The environmental isolates included 7 R. mucilaginosa isolates and 1 R. slooffiae isolate. Using the ID32C system, along with a nitrate assimilation test, only 90.3% of the isolates tested were correctly identified. In the biofilm formation assay, R. mucilaginosa and R. minuta exhibited greater biofilm formation ability compared to the other Rhodotorula species; the clinical isolates of R. mucilaginosa showed greater biofilm formation compared to the environmental isolates (P = 0.04). Amphotericin B showed good in vitro activity (MIC ≤ 1 μg/ml) against planktonic cells, whereas voriconazole and posaconazole showed poor activity (MIC50/MIC90, 2/4 μg/ml). Caspofungin and fluconazole MICs were consistently high for all isolates tested (≥64 μg/ml and ≥ 4 μg/ml, respectively). In this study, we emphasized the importance of molecular methods to correctly identify Rhodotorula species isolates and non-R. mucilaginosa species in particular. The antifungal susceptibility profile reinforces amphotericin B as the antifungal drug of choice for the treatment of Rhodotorula infections. To our knowledge, this is the first study evaluating putative differences in the ability of biofilm formation among different Rhodotorula species.

Genetic diversity of medically important and emerging Candida species causing invasive infection

BackgroundGenetic variation in the ribosomal DNA (rDNA) internal transcribed spacer (ITS) region has been studied among fungi. However, the numbers of ITS sequence polymorphisms in the various Candida species and their associations with sources of invasive fungal infections remain poorly investigated. Here, we characterized the intraspecific and interspecific ITS diversity of Candida spp. strains collected from patients with bloodstream or oroesophageal candidiasis.MethodsWe selected cultures of representative medically important species of Candida as well as some rare and emerging pathogens. Identification was performed by micromorphology and by biochemical testing using an ID32C® system, as well as by the sequencing of rDNA ITS. The presence of intraspecific ITS polymorphisms was characterized based on haplotype networks, and interspecific diversity was characterized based on Bayesian phylogenetic analysis.ResultsAmong 300 Candida strains, we identified 76 C. albicans, 14 C. dubliniensis, 40 C. tropicalis, 47 C. glabrata, 34 C. parapsilosis (sensu stricto), 31 C. orthopsilosis, 3 C. metapsilosis, 21 Meyerozyma guilliermondii (C. guilliermondii), 12 Pichia kudriavzevii (C. krusei), 6 Clavispora lusitaniae (C. lusitaniae), 3 C. intermedia, 6 Wickerhamomyces anomalus (C. pelliculosa), and 2 C. haemulonii strains, and 1 C. duobushaemulonii, 1 Kluyveromyces marxianus (C. kefyr), 1 Meyerozyma caribbica (C. fermentati), 1 Pichia norvegensis (C. norvegensis), and 1 Lodderomyces elongisporus strain. Out of a total of seven isolates with inconsistent ID32C® profiles, ITS sequencing identified one C. lusitaniae strain, three C. intermedia strains, two C. haemulonii strains and one C. duobushaemulonii strain. Analysis of ITS variability revealed a greater number of haplotypes among C. albicans, C. tropicalis, C. glabrata and C. lusitaniae, which are predominantly related to endogenous sources of acquisition. Bayesian analysis confirmed the major phylogenetic relationships among the isolates and the molecular identification of the different Candida spp.ConclusionsMolecular studies based on ITS sequencing are necessary to identify closely related and emerging species. Polymorphism analysis of the ITS rDNA region demonstrated its utility as a genetic marker for species identification and phylogenetic relationships as well as for drawing inferences concerning the natural history of hematogenous infections caused by medically important and emerging Candida species.

Accurate Identification of Candida parapsilosis (Sensu Lato) by Use of Mitochondrial DNA and Real-Time PCR

ABSTRACT Candida parapsilosis is the Candida species isolated the second most frequently from blood cultures in South America and some European countries, such as Spain. Since 2005, this species has been considered a complex of 3 closely related species: C. parapsilosis, Candida metapsilosis, and Candida orthopsilosis. Here, we describe a real-time TaqMan-MGB PCR assay, using mitochondrial DNA (mtDNA) as the target, which readily distinguishes these 3 species. We first used comparative genomics to locate syntenic regions between these 3 mitochondrial genomes and then selected NADH5 as the target for the real-time PCR assay. Probes were designed to include a combination of different single-nucleotide polymorphisms that are able to differentiate each species within the C. parapsilosis complex. This new methodology was first tested using mtDNA and then genomic DNA from 4 reference and 5 clinical strains. For assay validation, a total of 96 clinical isolates and 4 American Type Culture Collection (ATCC) isolates previously identified by internal transcribed spacer (ITS) ribosomal DNA (rDNA) sequencing were tested. Real-time PCR using genomic DNA was able to differentiate the 3 species with 100% accuracy. No amplification was observed when DNA from other species was used as the template. We observed 100% congruence with ITS rDNA sequencing identification, including for 30 strains used in blind testing. This novel method allows a quick and accurate intracomplex identification of C. parapsilosis and saves time compared with sequencing, which so far has been considered the “gold standard” for Candida yeast identification. In addition, this assay provides a useful tool for epidemiological and clinical studies of these emergent species.

The hidden factors in impact factors: a perspective from Brazilian science

Contrary to all current international recommendations on evaluation of academic achievement the evaluation of graduate programs in Brazil relies heavily on journal impact factors (Garfield, 2006; San Francisco Declaration on Research Assessment, 2013). The governmental agency CAPES from the Education Ministry monopolize this evaluation and pressure programs by the distribution of funding resources and departmental fellowships conditioned to adherence to a journal classification system called “Qualis” which is a discretization of the continuous distribution of journals ranking by their impact factors (Greenwood, 2007). In several institutions the graduate committee authorizes professors to act as thesis advisors only if in a certain period (e.g., 4 years) they publish at least one paper in a journal classified as “Qualis A2.” This classification has seven categories with decreasing impact factor ranges (A1, A2, B1, B2, B3, B4, B5, and C, http://www.capes.gov.br/avaliacao/qualis) and has certain percentile adjustments depending on the field of research. This has been hailed as a major cause for the enhancement of Brazilian scientific output although this system has several critics, demanding a profound review of evaluation criteria, and a proper adaptation to international guidelines (Rocha-e-Silva, 2009; Hermes-Lima, 2013). Most critics suggest that the excessive concern with publication in certain journals is in fact reducing the originality. Recently, to dismay of CAPES, a survey of the Brazilian science publication output, from 2001 to 2011, has shown that although Brazil has climbed from 17th place to 13th place in the total number of papers it has dropped from 31th to 40th place in citations (Rughetti, 2013).

Phylogenetic evidence based on Trypanosoma cruzi nuclear gene sequences and information entropy suggest that inter-strain intragenic recombination is a basic mechanism underlying the allele diversity of hybrid strains

The diversity of Trypanosoma cruzi is categorized into six discrete typing units (DTUs) T. cruzi I to VI. Several studies indicate that T. cruzi I and II are ancestors of T. cruzi III-VI which are considered products of independent hybridization events. The individual haplotypes or alleles of these hybrids cluster in three groups, either closer to T. cruzi I or T. cruzi II or forming a midpoint clade between T. cruzi I and II in network phylogenies. To understand the origins of these different sets of haplotypes and test the hypothesis of a direct correlation between high entropy and positive selection, we analyzed four nuclear protein coding genes. We show that hybrid strains contain haplotypes that are mosaics probably originated by intragenic recombination. Accordingly, in phylogenies, the hybrid haplotypes are closer to one or both parentals (T. cruzi I and II) depending on the proportion of parental sequences composing the mosaics. In addition, Shannon entropy, used to measure sequence diversity, is highly correlated with positive selection in the four genes here analyzed. Our data on recombination patterns also support the hypothesis of two hybridization events in the hybrid structures of T. cruzi III-VI. Data presented and discussed here are consistent with a scenario where TcI and TcII are phylogenetically divergent forming a hybrid zone in between (T. cruzi III-VI). We predict that because of the quasi-random nature of T. cruzi I and II hybridization more DTUs, with different haplotype combinations, will be discovered in the hybrid zone.

Network genealogy of 195-bp satellite DNA supports the superimposed hybridization hypothesis of Trypanosoma cruzi evolutionary pattern.

Trypanosoma cruzi is highly diverse genetically and has been partitioned into six discrete typing units (DTUs), recently re-named T. cruzi I-VI. Although T. cruzi reproduces predominantly by binary division, accumulating evidence indicates that particular DTUs are the result of hybridization events. Two major scenarios for the origin of the hybrid lineages have been proposed. It is accepted widely that the most heterozygous TcV and TcVI DTUs are the result of genetic exchange between TcII and TcIII strains. On the other hand, the participation of a TcI parental in the current genome structure of these hybrid strains is a matter of debate. Here, sequences of the T. cruzi-specific 195-bp satellite DNA of TcI, TcII, TcIII, TcV, and TcVI strains have been used for inferring network genealogies. The resulting genealogy showed a high degree of reticulation, which is consistent with more than one event of hybridization between the Tc DTUs. The data also strongly suggest that TcIII is a hybrid with two distinct sets of satellite sequences, and that genetic exchange between TcI and TcII parentals occurred within the pedigree of the TcV and TcVI DTUs. Although satellite DNAs belong to the fast-evolving portion of eukaryotic genomes, in >100 satellite units of nine T. cruzi strains we found regions that display 100% identity. No DTU-specific consensus motifs were identified, inferring species-wide conservation.

Candida mesorugosa sp. nov., a novel yeast species similar to Candida rugosa, isolated from a tertiary hospital in Brazil

Candida rugosa is a yeast species that is emerging as a causative agent of invasive infection, particularly in Latin America. Recently, C. pseudorugosa was proposed as a new species closely related to C. rugosa. We evaluated in this investigation the genetic heterogeneity within the C. rugosa species complex. All clinical isolates used in this study were identified phenotypically as C. rugosa but were genotypically different from the C. rugosa type, ATCC 10571. RAPD marker analysis revealed less than 83% similarity between our clinical isolates and the C. rugosa type strain. The D1/D2 region sequences of our clinical isolates showed 98% identity with C. rugosa but only 94-95% identity with C. pseudorugosa. The ITS rDNA sequences of the Brazilian isolates showed 91% identity with the C. rugosa ATCC 10571 ITS sequence. Network and Bayesian analyses of ITS and housekeeping gene sequences separated our clinical isolates into different branches from C. rugosa type strain. These differences are sufficient to reassign our isolates to a distinct species, named C. mesorugosa.

Intraspecific comparative genomics of Candida albicans mitochondria reveals non-coding regions under neutral evolution

The opportunistic fungal pathogen Candida albicans causes serious hematogenic hospital acquired candidiasis with worldwide impact on public health. Because of its importance as a nosocomial etiologic agent, C. albicans genome has been largely studied to identify intraspecific variation and several typing methods have been developed to distinguish closely related strains. Mitochondrial DNA can be useful for this purpose because, as compared to nuclear DNA, its higher mutational load and evolutionary rate readily reveals microvariants. Accordingly, we sequenced and assembled, with 8-fold coverage, the mitochondrial genomes of two C. albicans clinical isolates (L296 and L757) and compared these sequences with the genome sequence of reference strain SC5314. The genome alignment of 33,928 positions revealed 372 polymorphic sites being 230 in coding and 142 in non-coding regions. Three intergenic regions located between genes tRNAGly/COX1, NAD3/COB and ssurRNA/NAD4L, named IG1, IG2 and IG3, respectively, which showed high number of neutral substitutions, were amplified and sequenced from 18 clinical isolates from different locations in Latin America and 2 ATCC standard C. albicans strains. High variability of sequence and size were observed, ranging up to 56bp size difference and phylogenies based on IG1, IG2 and IG3 revealed three groups. Insertions of up to 49bp were observed exclusively in Argentinean strains relative to the other sequences which could suggest clustering by geographical polymorphism. Because of neutral evolution, high variability, easy isolation by PCR and full length sequencing these mitochondrial intergenic regions can contribute with a novel perspective in molecular studies of C. albicans isolates, complementing well established multilocus sequence typing methods.

A Model of Gene Expression Based on Random Dynamical Systems Reveals Modularity Properties of Gene Regulatory Networks

Here we propose a new approach to modeling gene expression based on the theory of random dynamical systems (RDS) that provides a general coupling prescription between the nodes of any given regulatory network given the dynamics of each node is modeled by a RDS. The main virtues of this approach are the following: (i) it provides a natural way to obtain arbitrarily large networks by coupling together simple basic pieces, thus revealing the modularity of regulatory networks; (ii) the assumptions about the stochastic processes used in the modeling are fairly general, in the sense that the only requirement is stationarity; (iii) there is a well developed mathematical theory, which is a blend of smooth dynamical systems theory, ergodic theory and stochastic analysis that allows one to extract relevant dynamical and statistical information without solving the system; (iv) one may obtain the classical rate equations form the corresponding stochastic version by averaging the dynamic random variables (small noise limit). It is important to emphasize that unlike the deterministic case, where coupling two equations is a trivial matter, coupling two RDS is non-trivial, specially in our case, where the coupling is performed between a state variable of one gene and the switching stochastic process of another gene and, hence, it is not a priori true that the resulting coupled system will satisfy the definition of a random dynamical system. We shall provide the necessary arguments that ensure that our coupling prescription does indeed furnish a coupled regulatory network of random dynamical systems. Finally, the fact that classical rate equations are the small noise limit of our stochastic model ensures that any validation or prediction made on the basis of the classical theory is also a validation or prediction of our model. We illustrate our framework with some simple examples of single-gene system and network motifs.