Sistemas de informação para DNA microarrays

O projecto de sequenciação do genoma humano veio abrir caminho para o surgimento de novas áreas transdisciplinares de investigação, como a biologia computacional, a bioinformática e a bioestatística. Um dos resultados emergentes desde advento foi a tecnologia de DNA microarrays, que permite o estudo do perfil da expressão de milhares de genes, quando sujeitos a perturbações externas. Apesar de ser uma tecnologia relativamente consolidada, continua a apresentar um conjunto vasto de desafios, nomeadamente do ponto de vista computacional e dos sistemas de informação. São exemplos a optimização dos procedimentos de tratamento de dados bem como o desenvolvimento de metodologias de interpretação semi-automática dos resultados. O principal objectivo deste trabalho consistiu em explorar novas soluções técnicas para agilizar os procedimentos de armazenamento, partilha e análise de dados de experiências de microarrays. Com esta finalidade, realizou-se uma análise de requisitos associados às principais etapas da execução de uma experiência, tendo sido identificados os principais défices, propostas estratégias de melhoramento e apresentadas novas soluções. Ao nível da gestão de dados laboratoriais, é proposto um LIMS (Laboratory Information Management System) que possibilita a gestão de todos os dados gerados e dos procedimentos realizados. Este sistema integra ainda uma solução que permite a partilha de experiências, de forma a promover a participação colaborativa de vários investigadores num mesmo projecto, mesmo usando LIMS distintos. No contexto da análise de dados, é apresentado um modelo que facilita a integração de algoritmos de processamento e de análise de experiências no sistema desenvolvido. Por fim, é proposta uma solução para facilitar a interpretação biológica de um conjunto de genes diferencialmente expressos, através de ferramentas que integram informação existente em diversas bases de dados biomédicas.

Structural analysis of glycans and development of microarrays from Helcobacter pylori cell surface glycome

Helicobacter pylori is a bacterial pathogen that affects more than half of the world’s population with gastro-intestinal diseases and is associated with gastric cancer. The cell surface of H. pylori is decorated with lipopolysaccharides (LPSs) composed of three distinct regions: a variable polysaccharide moiety (O-chain), a structurally conserved core oligosaccharide, and a lipid A region that anchors the LPS to the cell membrane. The O-chain of H. pylori LPS, exhibits unique oligosaccharide structures, such as Lewis (Le) antigens, similar to those present in the gastric mucosa and are involved in interactions with the host. Glucan, heptoglycan, and riban domains are present in the outer core region of some H. pylori LPSs. Amylose-like glycans and mannans are also constituents of some H. pylori strains, possibly co-expressed with LPSs. The complexity of H. pylori LPSs has hampered the establishment of accurate structure-function relationships in interactions with the host, and the design of carbohydrate-based therapeutics, such as vaccines. Carbohydrate microarrays are recent powerful and sensitive tools for studying carbohydrate antigens and, since their emergence, are providing insights into the function of carbohydrates and their involvement in pathogen-host interactions. The major goals of this thesis were the structural analysis of LPSs from H. pylori strains isolated from gastric biopsies of symptomatic Portuguese patients and the construction of a novel pathogen carbohydrate microarray of these LPSs (H. pylori LPS microarray) for interaction studies with proteins. LPSs were extracted from the cell surface of five H. pylori clinical isolates and one NCTC strain (26695) by phenol/water method, fractionated by size exclusion chromatography and analysed by gas chromatography coupled to mass spectrometry. The oligosaccharides released after mild acid treatment of the LPS were analysed by electrospray mass spectrometry. In addition to the conserved core oligosaccharide moieties, structural analyses revealed the presence of type-2 Lex and Ley antigens and N-acetyllactosamine (LacNAc) sequences, typically found in H. pylori strains. Also, the presence of O-6 linked glucose residues, particularly in LPSs from strains 2191 and NCTC 26695, pointed out to the expression of a 6-glucan. Other structural domains, namely ribans, composed of O-2 linked ribofuranose residues were observed in the LPS of most of H. pylori clinical isolates. For the LPS from strain 14382, large amounts of O-3 linked galactose units, pointing to the occurrence of a galactan, a domain recently identified in the LPS of another H. pylori strain. A particular feature to the LPSs from strains 2191 and CI-117 was the detection of large amounts of O-4 linked N-acetylglucosamine (GlcNAc) residues, suggesting the presence of chitin-like glycans, which to our knowledge have not been described for H. pylori strains. For the construction of the H. pylori LPS microarray, the structurally analysed LPSs, as well as LPS-derived oligosaccharide fractions, prepared as neoglycolipid (NGL) probes were noncovalently immobilized onto nitrocellulosecoated glass slides. These were printed together with NGLs of selected sequence defined oligosaccharides, bacterial LPSs and polysaccharides. The H. pylori LPS microarray was probed for recognition with carbohydratebinding proteins (CBPs) of known specificity. These included Le and blood group-related monoclonal antibodies (mAbs), plant lectins, a carbohydratebinding module (CBM) and the mammalian immune receptors DC-SIGN and Dectin-1. The analysis of these CBPs provided new information that complemented the structural analyses and was valuable in the quality control of the constructed microarray. Microarray analysis revealed the occurrence of type-2 Lex and Ley, but not type-1 Lea or Leb antigens, supporting the results obtained in the structural analysis. Furthermore, the H. pylori LPSs were recognised by DC-SIGN, a mammalian lectin known to interact with this bacterium through fucosylated Le epitopes expressed in its LPSs. The -fucose-specific lectin UEA-I, showed restricted binding to probes containing type-2 blood group H sequence and to the LPSs from strains CI-117 and 14382. The presence of H-type-2, as well Htype- 1 in the LPSs from these strains, was confirmed using specific mAbs. Although H-type-1 determinant has been reported for H. pylori LPSs, this is the first report of the presence of H-type-2 determinant. Microarray analysis also revealed that plant lectins known to bind 4-linked GlcNAc chitin oligosaccharide sequences bound H. pylori LPSs. STL, which exhibited restricted and strong binding to 4GlcNAc tri- and pentasaccharides, differentially recognised the LPS from the strain CI-117. The chitin sequences recognised in the LPS could be internal, as no binding was detected to this LPS with WGA, known to be specific for nonreducing terminal of 4GlcNAc sequence. Analyses of the H. pylori LPSs by SDS-PAGE and Western blot with STL provided further evidence for the presence of these novel domains in the O-chain region of this LPS. H. pylori LPS microarray was also applied to analysis of two human sera. The first was from a case infected with H. pylori (H. pylori+ CI-5) and the second was from a non-infected control.The analysis revealed a higher IgG-reactivity towards H. pylori LPSs in the H. pylori+ serum, than the control serum. A specific IgG response was observed to the LPS isolated from the CI-5 strain, which caused the infection. The present thesis has contributed to extension of current knowledge on chemical structures of LPS from H. pylori clinical isolates. Furthermore, the H. pylori LPS microarray constructed enabled the study of interactions with host proteins and showed promise as a tool in serological studies of H. pyloriinfected individuals. Thus, it is anticipated that the use of these complementary approaches may contribute to a better understanding of the molecular complexity of the LPSs and their role in pathogenesis.

Nucella lapillus: imposex transcriptome analysis and phenotypic plasticity

O conhecimento de mecanismos de genómica funcional tem sido maioritariamente adquirido pela utilização de organismos modelo que são mantidos em condições laboratoriais. Contudo, estes organismos não reflectem as respostas a alterações ambientais. Por outro lado, várias espécies, ecologicamente bem estudadas, reflectem bem as interacções entre genes e ambiente mas que, das quais não existem recursos genéticos disponíveis. O imposex, caracterizado pela superimposição de caracteres sexuais masculinos em fêmeas, é induzido pelo tributilestanho (TBT) e trifenilestanho (TPT) e representa um dos melhores exemplos de disrupção endócrina com causas antropogénicas no ambiente aquático. Com o intuito de elucidar as bases moleculares deste fenómeno, procedeu-se à combinação das metodologias de pirosequenciação (sequenciação 454 da Roche) e microarrays (Agilent 4*180K) de forma a contribuir para um melhor conhecimento desta interacção gene-ambiente no gastrópode Nucella lapillus, uma espécie sentinela para imposex. O trancriptoma de N. lapillus foi sequenciado, reconstruído e anotado e posteriormente utilizado para a produção de um “array” de nucleótidos. Este array foi então utilizado para explorar níveis de expressão génica em resposta à contaminação por TBT. Os resultados obtidos confirmaram as hipóteses anteriormente propostas (esteróidica, neuroendócrina, retinóica) e adicionalmente revelou a existência de potenciais novos mecanismos envolvidos no fenómeno imposex. Evidência para alvos moleculares de disrupção endócrina não relacionados com funções reprodutoras, tais como, sistema imunitário, apoptose e supressores de tumores, foram identificados. Apesar disso, tendo em conta a forte componente reprodutiva do imposex, esta componente funcional foi a mais explorada. Assim, factores de transcrição e receptores nucleares lipofílicos, funções mitocondriais e actividade de transporte celular envolvidos na diferenciação de géneros estão na base de potenciais novos mecanismos associados ao imposex em N. lapillus. Em particular, foi identificado como estando sobre-expresso, um possível homólogo do receptor nuclear “peroxisome proliferator-activated receptor gamma” (PPARγ), cuja função na indução de imposex foi confirmada experimentalmente in vivo após injecção dos animais com Rosiglitazone, um conhecido ligando de PPARγ em vertebrados. De uma forma geral, os resultados obtidos mostram que o fenómeno imposex é um mecanismo complexo, que possivelmente envolve a cascata de sinalização envolvendo o receptor retinoid X (RXR):PPARγ “heterodimer” que, até à data não foi descrito em invertebrados. Adicionalmente, os resultados obtidos apontam para alguma conservação de mecanismos de acção envolvidos na disrupção endócrina em invertebrados e vertebrados. Finalmente, a informação molecular produzida e as ferramentas moleculares desenvolvidas contribuem de forma significativa para um melhor conhecimento do fenómeno imposex e constituem importantes recursos para a continuação da investigação deste fenómeno e, adicionalmente, poderão vir a ser aplicadas no estudo de outras respostas a alterações ambientais usando N. lapillus como organismo modelo. Neste sentido, N. lapillus foi também utilizada para explorar a adaptação na morfologia da concha em resposta a alterações naturais induzidas por acção das ondas e pelo risco de predação por caranguejos. O contributo da componente genética, plástica e da sua interacção para a expressão fenotípica é crucial para compreender a evolução de caracteres adaptativos a ambientes heterogéneos. A contribuição destes factores na morfologia da concha de N. lapillus foi explorada recorrendo a transplantes recíprocos e experiências laboratoriais em ambiente comum (com e sem influência de predação) e complementada com análises genéticas, utilizando juvenis provenientes de locais representativos de costas expostas e abrigadas da acção das ondas. As populações estudadas são diferentes geneticamente mas possuem o mesmo cariótipo. Adicionalmente, análises morfométricas revelaram plasticidade da morfologia da concha em ambas as direcções dos transplantes recíprocos e também a retenção parcial, em ambiente comum, da forma da concha nos indivíduos da F2, indicando uma correlação positiva (co-gradiente) entre heritabilidade e plasticidade. A presença de estímulos de predação por caranguejos estimulou a produção de conchas com labros mais grossos, de forma mais evidente em animais recolhidos de costas expostas e também provocou alterações na forma da concha em animais desta proveniência. Estes dados sugerem contra-gradiente em alterações provocadas por predação na morfologia da concha, na produção de labros mais grossos e em níveis de crescimento. O estudo das interacções gene-ambiente descritas acima demonstram a actual possibilidade de produzir recursos e conhecimento genómico numa espécie bem caracterizada ecologicamente mas com limitada informação genómica. Estes recursos permitem um maior conhecimento biológico desta espécie e abrirão novas oportunidades de investigação, que até aqui seriam impossíveis de abordar.

Toxicogenomics of natural and anthropogenic effectors in eleutheroembryos

The introduction of chemicals into the environment by human activities may represent a serious risk to environmental and human health. Environmental risk assessment requires the use of efficient and sensitive tools to determine the impact of contaminants on the ecosystems. The use of zebrafish for the toxicity assessment of pharmaceuticals, drugs, and pollutants, is becoming well accepted due to zebrafish unique advantages for the screening of compounds for hazard identification. The aim of the present work is to apply toxicogenomic approaches to identify novel biomarkers and uncovered potential modes of action of classic and emergent contaminants able to disrupt endocrine systems, such as the Retinoic Acid Receptor, Retinoid X Receptor and the Aryl Hydrocarbon Receptor. This study relies on different nuclear and cytosolic protein receptors and other conditional (ligand- or stress- activated) transcriptional factors that are intimately involved in the regulation of defensome genes and in mechanisms of chemical toxicity. The transcriptomic effects of organic compounds, endogenous compounds, and nanoparticles were analysed during the early stages of zebrafish development. Studying the gene expression profiles of exposed and unexposed organisms to pollutants using microarrays allowed the identification of specific gene markers and to establish a "genetic code" for the tested compounds. Changes in gene expression were observed at toxicant concentrations that did not cause morphological effects. Even at low toxicant concentrations, the observed changes in transcript levels were robust for some target genes. Microarray responses of selected genes were further complemented by the real time quantitative polymerase chain reaction (qRT-PCR) methodology. The combination of bio-informatic, toxicological analyses of differential gene expression profiles, and biochemical and phenotypic responses across the treatments allowed the identification of uncovered potential mechanisms of action. In addition, this work provides an integrated set of tools that can be used to aid management-decision making by improving the predictive capability to measure environmental stress of contaminants in freshwater ecosystems. This study also illustrates the potential of zebrafish embryos for the systematic, large-scale analysis of chemical effects on developing vertebrates.

Effect assessment of nanoparticles toxicity in the terrestrial compartment

Over 11 million tons of nanomaterials (NMs) have been produced in 2012 and predictions point the increase in production. Despite predictions and extended usage via consumer products and industry, the understanding of the potential impact of these materials on the environment is virtually absent. The main aim of this thesis is to understand how a selected group of nanomaterials (metal based particles) may impact soil invertebrates, with special focus on the mechanisms of response. Since a case-by-case Environmental Risk Assessment (ERA) of all the emerging contaminants (particularly NMs) is impossible, among others due to time and cost reasons, to gain understanding on the mechanism of action and response is very important to reach a common paradigm. Understanding the modes of action provides predictive characters in cross particle extrapolation. Besides, it also provides insight for the production of new and sustainable materials. Overall, the effects of the selected NMs (Copper and Silver, Titanium and Zirconium oxides) and the respective salt forms, were investigated at the gene expression (using high-throughput tools, microarray and qPCR technology), biochemical (using enzymatic assays for analysis of oxidative stress markers) and organism (survival and reproduction as in OECD test guidelines) levels, this using standard soil species (Enchytraeus albidus, Enchytraeus crypticus, Eisenia fetida). Gene expression analysis provided valuable information on the mechanisms affected by each of the NMs. The gene expression profile highlighted a (nano)material signature and the effect of the duration of exposure. The functional analyses integrated with the biochemical and organism data, revealed a good understanding power. The biochemical parameters (oxidative stress related) were distinct across the materials and also influenced by duration of exposure and concentration. The standardized organismal responses differed the least between the various materials. The overall outcome is that, in this context of NMs effect assessment, gene expression and enzymatic assays introduced a very important knowledge gap, which could not had been achieved by the standard organismal effects alone. A reoccurring issue with some metal based NMs is the possible dissolution and subsequent release of ions that then causes toxicity e.g. Cu-NPs or Ag-NPs release Cu2+ or Ag+. The oxidation state of the particles was investigated, although this was not the focus of the thesis. The study of fate, e.g. dissolution of NPs, is also only in its beginning and the appropriate techniques are currently being developed. The results showed a specific nanoparticle effect. The UV exposure with titanium dioxide nanoparticles increased its effect.

Saccharomycotin transcriptomics by next-generation sequencing

The non-standard decoding of the CUG codon in Candida cylindracea raises a number of questions about the evolutionary process of this organism and other species Candida clade for which the codon is ambiguous. In order to find some answers we studied the transcriptome of C. cylindracea, comparing its behavior with that of Saccharomyces cerevisiae (standard decoder) and Candida albicans (ambiguous decoder). The transcriptome characterization was performed using RNA-seq. This approach has several advantages over microarrays and its application is booming. TopHat and Cufflinks were the software used to build the protocol that allowed for gene quantification. About 95% of the reads were mapped on the genome. 3693 genes were analyzed, of which 1338 had a non-standard start codon (TTG/CTG) and the percentage of expressed genes was 99.4%. Most genes have intermediate levels of expression, some have little or no expression and a minority is highly expressed. The distribution profile of the CUG between the three species is different, but it can be significantly associated to gene expression levels: genes with fewer CUGs are the most highly expressed. However, CUG content is not related to the conservation level: more and less conserved genes have, on average, an equal number of CUGs. The most conserved genes are the most expressed. The lipase genes corroborate the results obtained for most genes of C. cylindracea since they are very rich in CUGs and nothing conserved. The reduced amount of CUG codons that was observed in highly expressed genes may be due, possibly, to an insufficient number of tRNA genes to cope with more CUGs without compromising translational efficiency. From the enrichment analysis, it was confirmed that the most conserved genes are associated with basic functions such as translation, pathogenesis and metabolism. From this set, genes with more or less CUGs seem to have different functions. The key issues on the evolutionary phenomenon remain unclear. However, the results are consistent with previous observations and shows a variety of conclusions that in future analyzes should be taken into consideration, since it was the first time that such a study was conducted.

Study of tRNA modifying enzymes and codon usage bias in cancer

Recent evidences indicate that tRNA modifications and tRNA modifying enzymes may play important roles in complex human diseases such as cancer, neurological disorders and mitochondrial-linked diseases. We postulate that expression deregulation of tRNA modifying enzymes affects the level of tRNA modifications and, consequently, their function and the translation efficiency of their tRNA corresponding codons. Due to the degeneracy of the genetic code, most amino acids are encoded by two to six synonymous codons. This degeneracy and the biased usage of synonymous codons cause alterations that can span from protein folding to enhanced translation efficiency of a select gene group. In this work, we focused on cancer and performed a meta-analysis study to compare microarray gene expression profiles, reported by previous studies and evaluate the codon usage of different types of cancer where tRNA modifying enzymes were found de-regulated. A total of 36 different tRNA modifying enzymes were found de-regulated in most cancer datasets analyzed. The codon usage analysis revealed a preference for codons ending in AU for the up-regulated genes, while the down-regulated genes show a preference for GC ending codons. Furthermore, a PCA biplot analysis showed this same tendency. We also analyzed the codon usage of the datasets where the CTU2 tRNA modifying enzyme was found deregulated as this enzyme affects the wobble position (position 34) of specific tRNAs. Our data points to a distinct codon usage pattern between up and downregulated genes in cancer, which might be caused by the deregulation of specific tRNA modifying enzymes. This codon usage bias may augment the transcription and translation efficiency of some genes that otherwise, in a normal situation, would be translated less efficiently.