RENATA CARMONA E FERREIRA - Academia.edu (original) (raw)

Papers by RENATA CARMONA E FERREIRA

In this work we propose a model for gene expression based on the theory of random dynamical syste... more In this work we propose a model for gene expression based on the theory of random dynamical systems (RDS) and show that it has a "modularity property" in the following sense: given any collection of genes that are linked in a transcriptional network, if each of them is individually described by a certain class of RDS then there is a natural, and essentially unique, prescription for coupling them together, respecting the network topology, in such a way that the collective system formed by all genes is a RDS as well. Moreover, the class of RDS used to describe the individual genes is flexible enough to account for a wide range of stochastic behaviors within the realm of stationary processes.

Narcolepsy type 1 (NT1) is a rare and chronic neurological disease characterized by sudden sleep ... more Narcolepsy type 1 (NT1) is a rare and chronic neurological disease characterized by sudden sleep attacks, overwhelming daytime drowsiness, and cataplexy. To contribute to the understanding of NT1 genetic causes, here we describe a whole-genome analysis of a monozygotic twin pair discordant for NT1. Our study revealed that although both twins have the same pathogenic mutations in NT1 associated genes (such as HLA-DQB1*06:02:01, HLA-DRB1*11:01:02/*15:03:01) the unaffected twin has mutations in genes outside the HLA loci that could be suppressing the NT1 phenotype. These results support the notion that NT1 has an immunological basis but that protective mutations in non-HLA might interfere with the clinical manifestation of the disease.

The evolution of mitochondrial genomes is essential for the adaptation of yeasts to the variation... more The evolution of mitochondrial genomes is essential for the adaptation of yeasts to the variation of environmental levels of oxygen. Although Saccharomyces cerevisiae mitochondrial DNA lacks all complex I genes, respiration is possible because alternative NADH dehydrogenases are encoded by NDE1 and NDI1 nuclear genes. The proposed whole genome duplication (WGD) in the yeast ancestor at 150-100 million years ago caused nuclear gene duplications and secondary losses, although its relation to the loss of complex I mitocondrial is unknown. Here we present phylogenomic supertrees and supermatrix tree of 46 mitochondrial genomes showing that the loss of complex I predates WGD and occurred independently in the S. cerevisiae group and the fission yeast Schizosaccharomyces pombe. We also show that the branching patterns do not differ dramatically in supertrees and supermatrix phylogenies. Our inferences indicated consistent relations between conserved mitochondrial chromosomal gene order (sy...

Mathematical Biosciences, 2016

Here we propose a new approach to modeling gene expression based on the theory of random dynamica... more 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.

Genetics and molecular research : GMR, 2007

The analysis of transcriptional temporal noise could be an interesting means to study gene expres... more The analysis of transcriptional temporal noise could be an interesting means to study gene expression dynamics and stochasticity in eukaryotes. To study the statistical distributions of temporal noise in the eukaryotic model system Saccharomyces cerevisiae, we analyzed microarray data corresponding to one cell cycle for 6200 genes. We found that the temporal noise follows a lognormal distribution with scale invariance at the genome, chromosomal and sub-chromosomal levels. Correlation of temporal noise with the codon adaptation index suggests that at least 70% of all protein-coding genes are a noise minimization core of the genome. Accordingly, a mathematical model of individual gene expression dynamics was proposed, using an operator theoretical approach, which reveals strict conditions for noise variability and a possible global noise minimization/optimization strategy at the genome level. Our model and data show that minimal noise does not correspond to genes obeying a strictly de...

Medical mycology, Jan 22, 2015

Human and animal fungal pathogens are a growing threat worldwide leading to emerging infections a... more Human and animal fungal pathogens are a growing threat worldwide leading to emerging infections and creating new risks for established ones. There is a growing need for a rapid and accurate identification of pathogens to enable early diagnosis and targeted antifungal therapy. Morphological and biochemical identification methods are time-consuming and require trained experts. Alternatively, molecular methods, such as DNA barcoding, a powerful and easy tool for rapid monophasic identification, offer a practical approach for species identification and less demanding in terms of taxonomical expertise. However, its wide-spread use is still limited by a lack of quality-controlled reference databases and the evolving recognition and definition of new fungal species/complexes. An international consortium of medical mycology laboratories was formed aiming to establish a quality controlled ITS database under the umbrella of the ISHAM working group on "DNA barcoding of human and animal pa...

BMC Infectious Diseases, 2015

Background: Genetic variation in the ribosomal DNA (rDNA) internal transcribed spacer (ITS) regio... more Background: Genetic 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. Methods: We 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. Results: Among 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. Conclusions: Molecular 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.

Infection, Genetics and Evolution, 2015

Benznidazole (BZ) is one of the two drugs for Chagas disease treatment. In a previous study we sh... more Benznidazole (BZ) is one of the two drugs for Chagas disease treatment. In a previous study we showed that the Trypanosoma cruzi ABCG-like transporter gene, named TcABCG1, is over-expressed in parasite strains naturally resistant to BZ and that the gene of TcI BZ-resistant strains exhibited several single nucleotide polymorphisms (SNPs) as compared to the gene of CL Brener BZ-susceptible strain. Here we report the sequence of TcABCG1 gene of fourteen T. cruzi strains, with diverse degrees of BZ sensitivity and belonging to different discrete typing units (DTUs) and Tcbat group. Although DTU-specific SNPs and amino acid changes were identified, no direct correlation with BZ-resistance phenotype was found. Thus, it is plausible that the transporter abundance is a determinant factor for drug resistance, as pointed out above. Sequence data were used for Bayesian phylogenies and network genealogy analysis. The network showed a high degree of reticulation suggesting genetic exchange between the parasites. TcI and TcII clades were clearly separated. Tcbat sequences were close to TcI. A fourth clade clustered TcABCG1 haplotypes of TcV, TcVI and TcIII strains, with closer proximity to TcI. Analysis of the recombination patterns indicated that hybrid strains contain haplotypes that are mosaics most likely derived by intragenic recombination of parental sequences. The data confirm that TcII and TcIII as the parentals of TcV and TcVI DTUs. Since genetic fingerprint of TcI was found in TcIII, we sustain the previously proposed ''Two Hybridization model'' for the origin of hybrid strains. Among the twenty best BLASTP hits in databases, orthologues of TcABCG1 transporter were found in Leishmania spp. and African trypanosomes, though their function remains undescribed.

Frontiers in genetics, 2013

Examining marginal sequence similarities between bacterial type III secretion system components a... more Examining marginal sequence similarities between bacterial type III secretion system components and Trypanosoma cruzi surface proteins: horizontal gene transfer or convergent evolution? Front. Genet. 4:143.

Frontiers in Genetics, 2013

Frontiers in Cellular and Infection Microbiology, 2014

Enolase is secreted by Candida albicans and is present in its biofilms although its extracellular... more Enolase is secreted by Candida albicans and is present in its biofilms although its extracellular function is unknown. Here we show that extracellular enolase mediates the colonization of small intestine mucosa by C. albicans. Assays using intestinal mucosa disks show that C. albicans adhesion is inhibited, in a dose dependent mode, either by pretreatment of intestinal epithelium mucosa disks with recombinant C. albicans enolase (70% at 0.5 mg/ml enolase) or by pretreatment of C. albicans yeasts with anti-enolase antibodies (48% with 20 μg antiserum). Also using flow cytometry, immunoblots of conditioned media and confocal microscopy we demonstrate that enolase is present in biofilms and that the extracellular enolase is not an artifact due to cell lysis, but must represent functional secretion of a stable form. This is the first direct evidence that C. albicans' extracellular enolase mediates colonization on its primary translocation site. Also, because enolase is encoded by a single locus in C. albicans, its dual role peptide, as glycolytic enzyme and extracellular peptide, is a remarkable example of gene sharing in fungi.

Medical Mycology, 2013

in this region. On the other hand, C. glabrata is a highly prevalent etiological agent of candide... more in this region. On the other hand, C. glabrata is a highly prevalent etiological agent of candidemia in the USA and several European countries [3-8]. Candida rugosa has been extensively studied because of its biotechnological properties [9-13]. This species is able to produce extracellular lipases that have been used for the production of many products, such as fatty acids, synthesis of various esters, and kinetic resolution of racemic mixtures [14]. Candida rugosa was originally isolated from human faeces in 1917 and named Mycoderma rugosa [15]. It has been described as a colonising microorganism specifi cally associated with bovine mastitis [16]. Recently, C. rugosa has been recognised as a likely emergent human pathogen [17,18]. In the ARTEMIS DISK Antifungal Surveillance Program, Candida rugosa comprised 0.4% of 134,715 Candida spp. isolates from 127 different medical centres worldwide [18]. In 2003, we described an outbreak of C. rugosa among six hospitalized patients (fi ve out of six patients died) from a tertiary care hospital in Sao Paulo, Brazil [19]. The isolates exhibited a high degree of genotypic relatedness but were different from the

Journal of Clinical Microbiology, 2012

Candida parapsilosis is the Candida species isolated the second most frequently from blood cultur... more 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 Collec...

International Journal of Bioinformatics Research and Applications, 2009

Here we show that the transcriptional noise is an emergent property with scale invariance from ge... more Here we show that the transcriptional noise is an emergent property with scale invariance from genome level to the level of small Transcriptional Regulatory Genetic Networks (TRGN). We show that a small set of 9-12 genes reproduces the geometric mean value of transcriptional noise of the largest percolating networks and the whole 93-gene wide TRGN sub-network. Our results predict that the collapse of the standard deviation of the transcriptional noise as a function of gene sub-networks connectivity should occur for 1000 genes, the approximate size of the maximal interconnected percolating network cluster, which corresponds to the minimal genome size.

Infection, Genetics and Evolution, 2014

Chagas disease is caused by the protozoan Trypanosoma cruzi which affects 10 million people world... more Chagas disease is caused by the protozoan Trypanosoma cruzi which affects 10 million people worldwide. Very few kinases have been characterized in this parasite, including the phosphatidylinositol kinases (PIKs) that are at the heart of one of the major pathways of intracellular signal transduction. Recently, we have classified the PIK family in T. cruzi using five different models based on the presence of PIK conserved domains. In this study, we have mapped PIK genes to the chromosomes of two different T. cruzi lineages (G and CL Brener) and determined the cellular localization of two PIK members. The kinases have crucial roles in metabolism and are assumed to be conserved throughout evolution. For this reason, they should display a conserved localization within the same eukaryotic species. In spite of this, there is an extensive polymorphism regarding PIK localization at both genomic and cellular levels, among different T. cruzi isolates and between T. cruzi and Trypanosoma brucei, respectively. We showed in this study that the cellular localization of two PIK-related proteins (TOR1 and 2) in the T. cruzi lineage is distinct from that previously observed in T. brucei. In addition, we identified a new PIK gene with peculiar feature, that is, it codes for a FYVE domain at N-terminal position. FYVE-PIK genes are phylogenetically distant from the groups containing exclusively the FYVE or PIK domain. The FYVE-PIK architecture is only present in trypanosomatids and in virus such as Acanthamoeba mimivirus, suggesting a horizontal acquisition. Our Bayesian phylogenetic inference supports this hypothesis. The exact functions of this FYVE-PIK gene are unknown, but the presence of FYVE domain suggests a role in membranous compartments, such as endosome. Taken together, the data presented here strengthen the possibility that trypanosomatids are characterized by extensive genomic plasticity that may be considered in designing drugs and vaccines for prevention of Chagas disease.

Infection, Genetics and Evolution, 2013

The opportunistic fungal pathogen Candida albicans causes serious hematogenic hospital acquired c... more 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. 34 ptc Catheter ? ?

Infection, Genetics and Evolution, 2012

The diversity of Trypanosoma cruzi is categorized into six discrete typing units (DTUs) T. cruzi ... more 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.

Infection, Genetics and Evolution, 2010

Trypanosoma cruzi is highly diverse genetically and has been partitioned into six discrete typing... more 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 &amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;gt;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.

DNA Research, 2013

We present the sequencing and annotation of the Leishmania (Leishmania) amazonensis genome, an et... more 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 0-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 heatshock 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.

Antimicrobial Agents and Chemotherapy, 2012

Rhodotorula species are emergent fungal pathogens capable of causing invasive infections, primari... more 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 ...