Genômica comparativa em ambiente computacional distribuído: aplicabilidade e potencial no estudo de homologia entre protozoários / Comparative genomics in a distributed computing scenario: applicability and potential on protozoa homology study

A inferência de homologia entre organismos é uma atividade da genômicacomparativa que possibilita compreender melhor a relação entre os mesmos e, porconseguinte, sua distância evolutiva. Especificamente, a identificação de genesortólogos, ou seja, aqueles que têm sua origem em um ancestral comum, permiteoferecer melhorias na anotação funcional de genes, uma vez que genes ortólogostendem a ter sua função conservada.Com a crescente disponibilidade de genomas através de técnicas de NGS, aconstrução e atualização de bases de dados de ortólogos representam um desafioconstante, pois demandam o estudo e identificação das relações entre os genes detais organismos, em um volume de dados cada vez mais extenso e a um custocomputacional cada vez mais elevado.Nesta tese propomos a solução para nuvem computacional elastic-OrthoSearch, umworkflow científico de genômica comparativa inspirado no OrthoSearch, responsávelpela inferência de homologia entre organismos com o uso de abordagem baseadaem melhores hits recíprocos e perfis de Markov.Também propomos uma metodologia para criação de bases de ortólogos construídaatravés do reuso do OrthoSearch. Esta metodologia mostrou-se capaz de alavancara oferta de grupos ortólogos e assim auxiliar, por exemplo, na identificação de alvosde protozoários...

Homology inference among organisms is a comparative genomics tasks which allowsfor a better understanding on how such organisms are related to each other and ontheir evolutionary distance. Specifically, the identification of orthologous genes –those who share a common ancestor – allows for functional gene annotationimprovements, as orthologous genes tend to preserve their functions.The increasing amount of genomic data provided by the NGS techniques makes theorthologous databases’ building and update processes a challenging task. It requiresthe identification and study of the organisms’ genes relationships, in an extensivedata volume and at an increasing computational cost.In this thesis we propose elastic-OrthoSearch, a cloud-enabled comparativegenomics scientific workflow, derived from OrthoSearch. It aims at providinghomology inference among organisms, in a reciprocal best hits and Markov profilesapproach.We also propose an improved orthologous database creation methodology built ontop of OrthoSearch. Such methodology has shown means to offer a broaderorthologous groups dataset, which could in turn aid on Protozoa target identification...

Genes de virulência e diversidade genética em Salmonella spp. isoladas de amostras de origem suína / Virulence genes and genetic diversity in Salmonella spp. isolated from samples of swine origin

A diversificação da produção industrial de alimentos de origem suína e o intercâmbio comercial de animais e seus derivados destinados ao consumo humano podem ser importantes disseminadores de sorovares de Salmonella spp. na cadeia alimentar. Objetivou-se avaliar em 86 cepas de Salmonella spp., isoladas em granja de terminação e no abate de suínos, a ocorrência de três genes de virulência (invA, agfA e lpfA), bem como a similaridade genética entre elas. A ocorrência do gene invA foi verificada em 100% das amostras. O gene lpfA foi detectado em 80,23% (69/86) das cepas, não foi detectado em S. Panama e estava presente em todas as cepas de S. Infantis. O gene agfA foi detectado em 63,95% (55/86) das amostras. S. Agona apresentou positividade para todos os genes de virulência estudados. A análise de homologia entre as cepas agrupou os diferentes sorovares em clusters. A similaridade foi independente do local de isolamento, o que demonstra a presença de clones ao longo da cadeia de produção e a existência de multiplicidade de fontes para a infecção dos animais, como a ração, e a contaminação cruzada das carcaças. A pesquisa de genes de virulência e a avaliação da proximidade gênica permitem a caracterização e um maior entendimento sobre cepas de Salmonella circulantes na cadeia produtiva de suínos e, assim, podem subsidiar medidas de controle durante o processo produtivo com o objetivo de garantir a saúde do consumidor.(AU)

The diversification of industrial food production of swine origin and trade of animals and their derivatives for human consumption may be important disseminators of serovars of Salmonella spp. in the food chain. This study aimed to evaluate 86 strains of Salmonella spp. isolated form in the finishing and slaughter of pigs, the occurrence of three virulence genes (invA, agfa and lpfA), as well as the genetic similarity between them. The occurrence of gene invA was observed in 100% of the samples. The gene lpfA was detected in 80.23% (69/86) strains and is not detected in S. Panama, but present in all strains of S. Infantis. The gene agfA was detected in 63.95% (55/86). S. Agona was positive for all virulence genes studied. The analysis of homology between the different serovars grouped the isolates in clusters. The similarity was regardless of the location of isolation, demonstrating the presence of clones along the production chain and that there are multiple sources for the infection of animals, such as feed, and cross-contamination of carcasses. A survey of virulence genes and evaluation of gene proximity allow characterization and better understanding of Salmonella strains circulating in the pig production chain, thus being able to support control measures during the production process in order to ensure consumer health.(AU)

Genes de virulência e diversidade genética em Salmonella spp. isoladas de amostras de origem suína / Virulence genes and genetic diversity in Salmonella spp. isolated from samples of swine origin

A diversificação da produção industrial de alimentos de origem suína e o intercâmbio comercial de animais e seus derivados destinados ao consumo humano podem ser importantes disseminadores de sorovares de Salmonella spp. na cadeia alimentar. Objetivou-se avaliar em 86 cepas de Salmonella spp., isoladas em granja de terminação e no abate de suínos, a ocorrência de três genes de virulência (invA, agfA e lpfA), bem como a similaridade genética entre elas. A ocorrência do gene invA foi verificada em 100% das amostras. O gene lpfA foi detectado em 80,23% (69/86) das cepas, não foi detectado em S. Panama e estava presente em todas as cepas de S. Infantis. O gene agfA foi detectado em 63,95% (55/86) das amostras. S. Agona apresentou positividade para todos os genes de virulência estudados. A análise de homologia entre as cepas agrupou os diferentes sorovares em clusters. A similaridade foi independente do local de isolamento, o que demonstra a presença de clones ao longo da cadeia de produção e a existência de multiplicidade de fontes para a infecção dos animais, como a ração, e a contaminação cruzada das carcaças. A pesquisa de genes de virulência e a avaliação da proximidade gênica permitem a caracterização e um maior entendimento sobre cepas de Salmonella circulantes na cadeia produtiva de suínos e, assim, podem subsidiar medidas de controle durante o processo produtivo com o objetivo de garantir a saúde do consumidor...

The diversification of industrial food production of swine origin and trade of animals and their derivatives for human consumption may be important disseminators of serovars of Salmonella spp. in the food chain. This study aimed to evaluate 86 strains of Salmonella spp. isolated form in the finishing and slaughter of pigs, the occurrence of three virulence genes (invA, agfa and lpfA), as well as the genetic similarity between them. The occurrence of gene invA was observed in 100% of the samples. The gene lpfA was detected in 80.23% (69/86) strains and is not detected in S. Panama, but present in all strains of S. Infantis. The gene agfA was detected in 63.95% (55/86). S. Agona was positive for all virulence genes studied. The analysis of homology between the different serovars grouped the isolates in clusters. The similarity was regardless of the location of isolation, demonstrating the presence of clones along the production chain and that there are multiple sources for the infection of animals, such as feed, and cross-contamination of carcasses. A survey of virulence genes and evaluation of gene proximity allow characterization and better understanding of Salmonella strains circulating in the pig production chain, thus being able to support control measures during the production process in order to ensure consumer health...

Overlapped sequence types (STs) and serogroups of avian pathogenic (APEC) and human extra-intestinal pathogenic (ExPEC) Escherichia coli isolated in Brazil.

Avian pathogenic Escherichia coli (APEC) strains belong to a category that is associated with colibacillosis, a serious illness in the poultry industry worldwide. Additionally, some APEC groups have recently been described as potential zoonotic agents. In this work, we compared APEC strains with extraintestinal pathogenic E. coli (ExPEC) strains isolated from clinical cases of humans with extra-intestinal diseases such as urinary tract infections (UTI) and bacteremia. PCR results showed that genes usually found in the ColV plasmid (tsh, iucA, iss, and hlyF) were associated with APEC strains while fyuA, irp-2, fepC sitDchrom, fimH, crl, csgA, afa, iha, sat, hlyA, hra, cnf1, kpsMTII, clpVSakai and malX were associated with human ExPEC. Both categories shared nine serogroups (O2, O6, O7, O8, O11, O19, O25, O73 and O153) and seven sequence types (ST10, ST88, ST93, ST117, ST131, ST155, ST359, ST648 and ST1011). Interestingly, ST95, which is associated with the zoonotic potential of APEC and is spread in avian E. coli of North America and Europe, was not detected among 76 APEC strains. When the strains were clustered based on the presence of virulence genes, most ExPEC strains (71.7%) were contained in one cluster while most APEC strains (63.2%) segregated to another. In general, the strains showed distinct genetic and fingerprint patterns, but avian and human strains of ST359, or ST23 clonal complex (CC), presented more than 70% of similarity by PFGE. The results demonstrate that some "zoonotic-related" STs (ST117, ST131, ST10CC, ST23CC) are present in Brazil. Also, the presence of moderate fingerprint similarities between ST359 E. coli of avian and human origin indicates that strains of this ST are candidates for having zoonotic potential.

Influence of overlapping genes on the evolution of human hepatitis B virus.

The aim of this work was to analyse the influence of overlapping genes on the evolution of hepatitis B virus (HBV). A differential evolutionary behaviour among genetic regions and clinical status was found. Dissimilar levels of conservation of the different protein regions could derive from alternative mechanisms to maintain functionality. We propose that, in overlapping regions, selective constraints on one of the genes could drive the substitution process. This would allow protein conservation in one gene by synonymous substitutions while mechanisms of tolerance to the change operate in the overlapping gene (e.g. usage of amino acids with high-degeneracy codons, differential codon usage and replacement by physicochemically similar amino acids). In addition, differential selection pressure according to the HBeAg status was found in all genes, suggesting that the immune response could be one of the factors that would constrain viral replication by interacting with different HBV proteins during the HBeAg(-) stage.

Modelagem conceitual do sistema de banco de dados ProteinWorldDB / Conceptual modeling of the database system ProteinWorldDB

Esta tese descreve o projeto conceitual do sistema de banco de dados ProteinWorldDB (PWDB). Um ponto importante da proposta do PWDB é permitir a construção de consultas e procedimentos no domínio da genômica comparativa sem a necessidade de comparação de sequências. Além disso, o PCG comparou milhões de sequências de proteína, incluindo o conjunto proteico total de centenas de genomas completos, utilizando programação dinâmica, e não um método heurístico, para os cálculos de similaridade. A estratégia do PCG, assim como a genômica, está fundamentada no conhecimento de que sequências biológicas por si só são pouco informativas; elas precisam ser analisadas a partir de um enfoque comparativo para a inferência de homologia. A comparação de sequências de diferentes organismos introduz uma perspectiva evolutiva ao processo, e o estudo comparativo de genomas completos pode ampliar a escala do conhecimento de um único processo biológico para o de sistemas biológicos complexos em células e organismos. Para responder eficientemente questões dessa natureza, o esquema conceitual apresentado associa bases de dados biológicos de referência aos índices de similaridade já pré-calculados e armazenados pelo PCGUtilizando um formato gráfico de fácil compreensão para representar conceitos e relacionamentos (diagrama ER), o esquema foi proposto para facilitar o planejamento de consultas e procedimentos por pesquisadores da área de genômica (sem conhecimento de linguagens de bancos de dados), assim como guiar o desenvolvimento e a implementação física do PWDB por profissionais da área de computação. Alguns exemplos são apresentados com o objetivo de demonstrar a utilização do esquema conceitual para a especificação de consultas e procedimentos, mesmo antes da existência de um esquema lógico...

This thesis describes the conceptual design of the database system ProteinWorldDB(PWDB). An important point of the PWDB proposal is to allow the construction of queriesand procedures in the field of comparative genomics without the need for sequencecomparison. Moreover, the PCG compared millions of protein sequences, including theentire set of proteins from hundreds of complete genomes using dynamic programming,rather than a heuristic method, for calculating similarityPCG‘s strategy, like that of genomic studies in general, is grounded in the knowledgethat biological sequences alone are uninformative. They need to be analyzed from acomparative approach to infer homology. The comparison of sequences from differentorganisms introduces an evolutionary perspective to the process and the comparativestudy of complete genomes can expand our knowledge from a single biological processall the way to complex biological systems in cells and organisms. To efficiently answerquestions of this nature, the conceptual schema links selected international referencebiological databases to similarity indexes already precomputed and stored by the PCG.By using an easily understandable graphic format to represent concepts andrelationships (ER diagram), the schema was proposed to help the design of queries andprocedures by genomic researchers (who may not have knowledge of databaselanguages) as well as to guide the development and physical implementation of thesystem by developers. Some examples are presented to demonstrate the use of theconceptual schema for specifying queries and procedures, even before the existence ofa logical schema. The schema can be easily extended. Additional modules can be inserted/removed toinclude other protein sequences comparisons projects that may benefit from theinformation provided by the schema´s central module. Likewise, new databases specificto different areas (-omics, for example) can be cross-referenced to the schema...

Sigma70 promoters in Escherichia coli: specific transcription in dense regions of overlapping promoter-like signals.

We present here a computational analysis showing that sigma70 house-keeping promoters are located within zones with high densities of promoter-like signals in Escherichia coli, and we introduce strategies that allow for the correct computer prediction of sigma70 promoters. Based on 599 experimentally verified promoters of E.coli K-12, we generated and evaluated more than 200 weight matrices optimizing different criteria to obtain the best recognition matrices. The alignments generating the best statistical models did not fully correspond with the canonical sigma70 model. However, matrices that correspond to such a canonical model performed better as tools for prediction. We tested the predictive capacity of these matrices on 250 bp long regions upstream of gene starts, where 90% of the known promoters occur. The computational matrix models generated an average of 38 promoter-like signals within each 250 bp region. In more than 50% of the cases, the true promoter does not have the best score within the region. We observed, in fact, that real promoters occur mostly within regions with high densities of overlapping putative promoters. We evaluated several strategies to identify promoters. The best one uses an intrinsic score of the -10 and -35 hexamers that form the promoter as well as an extrinsic score that uses the distribution of promoters from the start of the gene. We were able to identify 86% true promoters correctly, generating an average of 4.7 putative promoters per region as output, of which 3.7, on average, exist in clusters, as a series of overlapping potentially competing RNA polymerase-binding sites. As far as we know, this is the highest predictive capability reported so far. This high signal density is found mainly within regions upstream of genes, contrasting with coding regions and regions located between convergently transcribed genes. These results are consistent with experimental evidence that show the existence of multiple overlapping promoter sites that become functional under particular conditions. This density is probably the consequence of a rich number of vestiges of promoters in evolution. We suggest that transcriptional regulators as well as other functional promoters play an important role in keeping these latent signals suppressed.

Homology among extra-cryptic DNA bands and the typical plasmids in brazilian Yersina pestis strains

Yersina pestis, the etiologic agent of plague, harbors three well-characterized plasmids: pFra (90-110kb), pYV (70kb) and pPst (9.5kb). Furthermore, some extra-cryptic DNA bands have been observed in a number of wild strains from several foci of the world. Additional bands have also been reported in Brazilian strains. Looking for any relationship among these cryptic DNA bands and the three-prototypical plasmids, we analyzed twelve strains displaying different plasmid content. The DNA bands were hybridized by southern blot with probes directed at the genes cafl, lcrV and 'pla' located respectively on the plasmids pFra, pYV and pPst. The probes were constructed by PCR amplification and labeled with digoxigenin. The Pla probe hybridized with its target (pPst) and with bands of about 35 kb suggesting some homology among them. The Cafl probe hybridized with the target (pFra) as well as with higher bands. The LcrV also hybridized with the target (pYV) and both with the bands higher than pFra and the bands between pFra and pYV. These results suggest that the large-cryptic bands could represent some rearrangement, open circular or linearized forms of the pFras and pYV plasmids.