Smart design of scaffolds obtained by biofabrication for tissue engineering applications

A engenharia de tecidos é um domínio tecnológico emergente em rápido desenvolvimento que se destina a produzir substitutos viáveis para a restauração, manutenção ou melhoria da função dos tecidos ou órgãos humanos. Uma das estratégias mais predominantes em engenharia de tecidos envolve crescimento celular sobre matrizes de suporte (scaffolds), biocompatíveis e biodegradáveis. Estas matrizes devem possuir não só elevadas propriedades mecânicas e vasculares, mas também uma elevada porosidade. Devido à incompatibilidade destes dois parâmetros, é necessário desenvolver estratégias de simulação de forma a obter estruturas optimizadas. A previsão real das propriedades mecânicas, vasculares e topológicas das matrizes de suporte, produzidas por técnicas de biofabricação, é muito importante para as diversas aplicações em engenharia de tecidos. A presente dissertação apresenta o estado da arte da engenharia de tecidos, bem como as técnicas de biofabricação envolvidas na produção de matrizes de suporte. Para o design optimizado de matrizes de suporte foi adoptada uma metodologia de design baseada tanto em métodos de elementos finitos para o cálculo do comportamento mecânico, vascular e as optimizações topológicas, como em métodos analíticos para a validação das simulações estruturais utilizando dados experimentais. Considerando que as matrizes de suporte são estruturas elementares do tipo LEGO, dois tipos de famílias foram consideradas, superfícies não periódicas e as superfícies triplas periódicas que descrevem superfícies naturais. Os objectivos principais desta dissertação são: i) avaliar as técnicas existentes de engenharia de tecidos; ii) avaliar as técnicas existentes de biofabricação para a produção de matrizes de suporte; iii) avaliar o desempenho e comportamento das matrizes de suporte; iv) implementar uma metodologia de design de matrizes de suporte em variáveis tais como a porosidade, geometria e comportamento mecânico e vascular por forma a auxiliar o processo de design; e por fim, v) validar experimentalmente a metodologia adoptada.

Pregnancy and newborns disorders followed by urine metabolomics

Chapter 1 introduces the scope of the work by identifying the clinically relevant prenatal disorders and presently available diagnostic methods. The methodology followed in this work is presented, along with a brief account of the principles of the analytical and statistical tools employed. A thorough description of the state of the art of metabolomics in prenatal research concludes the chapter, highlighting the merit of this novel strategy to identify robust disease biomarkers. The scarce use of maternal and newborn urine in previous reports enlightens the relevance of this work. Chapter 2 presents a description of all the experimental details involved in the work performed, comprising sampling, sample collection and preparation issues, data acquisition protocols and data analysis procedures. The proton Nuclear Magnetic Resonance (NMR) characterization of maternal urine composition in healthy pregnancies is presented in Chapter 3. The urinary metabolic profile characteristic of each pregnancy trimester was defined and a 21-metabolite signature found descriptive of the metabolic adaptations occurring throughout pregnancy. 8 metabolites were found, for the first time to our knowledge, to vary in connection to pregnancy, while known metabolic effects were confirmed. This chapter includes a study of the effects of non-fasting (used in this work) as a possible confounder. Chapter 4 describes the metabolomic study of 2nd trimester maternal urine for the diagnosis of fetal disorders and prediction of later-developing complications. This was achieved by applying a novel variable selection method developed in the context of this work. It was found that fetal malformations (FM) (and, specifically those of the central nervous system, CNS) and chromosomal disorders (CD) (and, specifically, trisomy 21, T21) are accompanied by changes in energy, amino acids, lipids and nucleotides metabolic pathways, with CD causing a further deregulation in sugars metabolism, urea cycle and/or creatinine biosynthesis. Multivariate analysis models´ validation revealed classification rates (CR) of 84% for FM (87%, CNS) and 85% for CD (94%, T21). For later-diagnosed preterm delivery (PTD), preeclampsia (PE) and intrauterine growth restriction (IUGR), it is found that urinary NMR profiles have early predictive value, with CRs ranging from 84% for PTD (11-20 gestational weeks, g.w., prior to diagnosis), 94% for PE (18-24 g.w. pre-diagnosis) and 94% for IUGR (2-22 g.w. pre-diagnosis). This chapter includes results obtained for an ultraperformance liquid chromatography-mass spectrometry (UPLC-MS) study of pre-PTD samples and correlation with NMR data. One possible marker was detected, although its identification was not possible. Chapter 5 relates to the NMR metabolomic study of gestational diabetes mellitus (GDM), establishing a potentially predictive urinary metabolic profile for GDM, 2-21 g.w. prior to diagnosis (CR 83%). Furthermore, the NMR spectrum was shown to carry information on individual phenotypes, able to predict future insulin treatment requirement (CR 94%). Chapter 6 describes results that demonstrate the impact of delivery mode (CR 88%) and gender (CR 76%) on newborn urinary profile. It was also found that newborn prematurity, respiratory depression, large for gestational age growth and malformations induce relevant metabolic perturbations (CR 82-92%), as well as maternal conditions, namely GDM (CR 82%) and maternal psychiatric disorders (CR 91%). Finally, the main conclusions of this thesis are presented in Chapter 7, highlighting the value of maternal or newborn urine metabolomics for pregnancy monitoring and disease prediction, towards the development of new early and non-invasive diagnostic methods.

Metabolic signature of lung cancer: a metabolomic study of human tissues and biofluids

This thesis reports the application of metabolomics to human tissues and biofluids (blood plasma and urine) to unveil the metabolic signature of primary lung cancer. In Chapter 1, a brief introduction on lung cancer epidemiology and pathogenesis, together with a review of the main metabolic dysregulations known to be associated with cancer, is presented. The metabolomics approach is also described, addressing the analytical and statistical methods employed, as well as the current state of the art on its application to clinical lung cancer studies. Chapter 2 provides the experimental details of this work, in regard to the subjects enrolled, sample collection and analysis, and data processing. In Chapter 3, the metabolic characterization of intact lung tissues (from 56 patients) by proton High Resolution Magic Angle Spinning (HRMAS) Nuclear Magnetic Resonance (NMR) spectroscopy is described. After careful assessment of acquisition conditions and thorough spectral assignment (over 50 metabolites identified), the metabolic profiles of tumour and adjacent control tissues were compared through multivariate analysis. The two tissue classes could be discriminated with 97% accuracy, with 13 metabolites significantly accounting for this discrimination: glucose and acetate (depleted in tumours), together with lactate, alanine, glutamate, GSH, taurine, creatine, phosphocholine, glycerophosphocholine, phosphoethanolamine, uracil nucleotides and peptides (increased in tumours). Some of these variations corroborated typical features of cancer metabolism (e.g., upregulated glycolysis and glutaminolysis), while others suggested less known pathways (e.g., antioxidant protection, protein degradation) to play important roles. Another major and novel finding described in this chapter was the dependence of this metabolic signature on tumour histological subtype. While main alterations in adenocarcinomas (AdC) related to phospholipid and protein metabolisms, squamous cell carcinomas (SqCC) were found to have stronger glycolytic and glutaminolytic profiles, making it possible to build a valid classification model to discriminate these two subtypes. Chapter 4 reports the NMR metabolomic study of blood plasma from over 100 patients and near 100 healthy controls, the multivariate model built having afforded a classification rate of 87%. The two groups were found to differ significantly in the levels of lactate, pyruvate, acetoacetate, LDL+VLDL lipoproteins and glycoproteins (increased in patients), together with glutamine, histidine, valine, methanol, HDL lipoproteins and two unassigned compounds (decreased in patients). Interestingly, these variations were detected from initial disease stages and the magnitude of some of them depended on the histological type, although not allowing AdC vs. SqCC discrimination. Moreover, it is shown in this chapter that age mismatch between control and cancer groups could not be ruled out as a possible confounding factor, and exploratory external validation afforded a classification rate of 85%. The NMR profiling of urine from lung cancer patients and healthy controls is presented in Chapter 5. Compared to plasma, the classification model built with urinary profiles resulted in a superior classification rate (97%). After careful assessment of possible bias from gender, age and smoking habits, a set of 19 metabolites was proposed to be cancer-related (out of which 3 were unknowns and 6 were partially identified as N-acetylated metabolites). As for plasma, these variations were detected regardless of disease stage and showed some dependency on histological subtype, the AdC vs. SqCC model built showing modest predictive power. In addition, preliminary external validation of the urine-based classification model afforded 100% sensitivity and 90% specificity, which are exciting results in terms of potential for future clinical application. Chapter 6 describes the analysis of urine from a subset of patients by a different profiling technique, namely, Ultra-Performance Liquid Chromatography coupled to Mass Spectrometry (UPLC-MS). Although the identification of discriminant metabolites was very limited, multivariate models showed high classification rate and predictive power, thus reinforcing the value of urine in the context of lung cancer diagnosis. Finally, the main conclusions of this thesis are presented in Chapter 7, highlighting the potential of integrated metabolomics of tissues and biofluids to improve current understanding of lung cancer altered metabolism and to reveal new marker profiles with diagnostic value.

Matérias-primas da plataforma do Mondego para cerâmica

O presente trabalho divulga os resultados dos estudos levados a efeito nas matérias-primas argilosas que se integram tipologicamente na argila comum, ocorrentes na designada Plataforma do Mondego, Centro de Portugal, na região entre Miranda do Corvo, a SW, e Tábua, a NE. Para tal realizou-se a cartografia superficial dos sedimentos continentais do Cretácico e do Terciário, aí preservados, estabeleceram-se as colunas sedimentares, a nível regional e a nível local, nas áreas de maior potencial reconhecido e, caracterizaram-se as matérias-primas argilosas amostradas, em termos de composição, textura e aptidão cerâmica, objectivando a definição das zonas das jazidas com maior interesse industrial. A informação obtida e compilada, relativa ao caulino e ao feldspato é também apresentada, com vista a uma percepção da potencialidade destas matériasprimas cerâmicas ocorrentes na área. O estudo de caracterização dos materiais argilosos investigados teve como base 53 amostras obtidas em seis regiões da Plataforma do Mondego aqui definidas por motivo de organização do trabalho, segundo os critérios geográfico e geotectónico. As formações de Côja e de Campelo, do Terciário, constituem as unidades onde ocorrem litótipos produtivos, em termos de matéria-prima para a Cerâmica de Barro Vermelho ou Cerâmica de Construção. Relativamente às características texturais e composicionais, em síntese, apresentam-se os factos relevantes seguintes: A matéria-prima argilosa existente na área estudada materializa, granulometricamente, na maioria das amostras, silte, caindo os níveis amostrados com maior percentagem de argila, no domínio do silte argiloso. O défice em fracção argila implica limitações quanto à possibilidade de diversificação de produtos cerâmicos fabricados com estas matérias-primas tal-qual. Os diferentes métodos analíticos utilizados na caracterização mineralógica dos materiais amostrados confirmaram uma composição em termos dos minerais argilosos, consistindo de ilite/mica (em geral, o mais abundante), caulinite e esmectite, interestratificados e clorite. Os minerais não argilosos são quartzo (predominante), feldspato (sobretudo potássico) e hematite, com uma representação baixa. A mineralogia da fracção inferior a 2μm das amostras, não difere das amostras totais, salvo no teor mais elevado dos minerais argilosos e acentuada redução dos minerais não argilosos. Os resultados da análise química por fluorescência de raios X das amostras integrais correlacionam-se com as características mineralógicas observadas através das técnicas analíticas utilizadas. No respeitante às propriedades e comportamento cerâmico verifica-se: Os parâmetros relacionados com a plasticidade indicam que parte das pastas elaboradas com estes materiais argilosos têm uma trabalhabilidade aceitável, mas existem problemas de conformação e acentuada retracção num número significativo de amostras, devido à elevada plasticidade da maioria das amostras. A extrusão é satisfatória a óptima. A RMF e a retracção em seco assumem valores, respectivamente, moderados a baixos e moderados, embora seja necessário ter em conta o procedimento de extrusão dos provetes, sem dispositivo de vácuo. Todas as amostras foram sujeitas a cozedura a 900ºC, e um conjunto seleccionado foi cozido a 1000ºC e a 1100ºC. As fases mineralógicas ocorrentes após cozeduras a 900ºC e 1100ºC foram identificadas num conjunto de amostras, tendo-se evidenciado a coerência dessas fases, com a mineralogia das amostras em seco. Após cozedura a 900ºC, os valores de RMF das amostras satisfazem geralmente os valores mínimos, exigidos para o fabrico de tijolo, abobadilha e, com alguma frequência, de telha, como já se verificava com os valores daquela propriedade em seco. Os valores de retracção seco-cozido são em geral, modestos. A capacidade de absorção de água é maioritariamente elevada. A formação de vidro, sobretudo, condiciona o comportamento destas propriedades por cozedura dos provetes a 1000ºC e a 1100ºC. A coloração predominante em cru das matérias-primas argilosas amostradas é amarelo acastanhado a castanho avermelhado. Após cozedura a 900ºC, há um acentuado escurecimento e incremento no grau de vermelho. As cozeduras a 1000ºC e 1100ºC promovem escurecimento gradual, com ligeira influência na cor. O comportamento dos provetes após as cozeduras cerâmicas revelou-se homogéneo a cada uma das respectivas temperaturas, não se registando também defeitos significativos, nem eflorescências. A análise das amostras em termos composicionais e tecnológicos permitiu destacar as principais características e aspectos distintivos das matériasprimas, nas diferentes regiões definidas objectivando as diferentes potencialidades cerâmicas. Nesta abordagem comparativa foram consideradas só as amostras dos campos silte e silte arenoso, por serem aquelas com maior interesse para a Cerâmica de Construção. Nas colunas sedimentares das regiões de Tábua e Santa Quitéria constata-se a ocorrência de dois ritmos de sedimentação, que embora assumam características específicas em cada região, têm aspectos composicionais e tecnológicos em comum, traduzindo melhor aptidão cerâmica os ritmos inferiores. A conjugação da cartografia realizada com os estudos laboratoriais permite concluir que as regiões de Tábua e de Santa Quitéria serão as que têm maior potencial por explorar, em matéria-prima para Cerâmica de Construção, apesar da primeira já ser intensamente explorada. Na região de Tábua, as amostras têm como fases mineralógicas principais ilite e quartzo na mesma proporção média (35%) e caulinite (média=19%) que regista enriquecimento significativo na fracção argila (média=38%). As argilas desta região registam a cor em cru mais vermelha e pH mais ácido observados. O ritmo de sedimentação inferior, com esmectite e interestratificados e ligeiramente menos quartzoso, apresenta melhores propriedades cerâmicas. A composição mineralógica média das amostras da região de Santa Quitéria é próxima daquela da região de Tábua, mas menos caulinítica, em especial na sequência inferior, na qual o teor médio de caulinite na fracção argila (7%) é o mais baixo observado. Na região de Côja – Arganil, a actividade extractiva é significativa na bacia de Côja. Aqui, a exploração de novas áreas potenciais é condicionada pela cobertura conglomerática e por estruturas tectónicas e não tectónicas relacionadas com comportamento plástico. Na restante área desta região, os recursos argilosos são penalizados por material areno-conglomerático. As amostras desta região distinguem-se das restantes a nível textural pela maior fracção areia e mineralogicamente pela presença de clorite, teor reduzido de caulinite e elevado de feldspato. As potencialidades em barro vermelho na região de Sanguinheda não serão significativas, pois a Formação de Côja é predominantemente arcósica e a Formação de Campelo pouco espessa e conglomerática. As argilas com melhor aptidão cerâmica foram amostradas na região de Miranda do Corvo – Lousã. Na composição, estas amostras são as que contêm maior fracção de argila, maior teor de ilite e caulinite e menor de argilas expansivas. Registam os melhores valores nas características tecnológicas, nomeadamente a RMF e absorção de água. A consistência dos grupos amostrais e ritmos definidos em termos composicionais e tecnológicos é corroborada pelas técnicas de análise estatística multivariada aplicadas, que os identificam. As condições de amostragem na região de Tábua, permitiram a elaboração de uma coluna tipológica, na qual, a partir de uma caracterização expedita de amostras é possível identificar a sua afinidade com os ritmos argilosos definidos e, consequentemente, a sua situação na coluna sedimentar regional e aptidão cerâmica. A cartografia dos recursos argilosos potenciais, elaborada à escala 1:25.000, constitui um dos objectivos principais deste trabalho e nela constam unidades litológicas, nas quais são diferenciadas unidades argilosas tendo também em consideração a tipologia e guias mineralógicos resultantes da caracterização das amostras. As características composicionais destas unidades denotam diferente aptidão cerâmica e, portanto, permitem salientar as zonas com maior interesse económico.

Analytical methods to study fate of organic pollutants in environment

Environmental transport of pollutants comprises distinct processes such as volatilization, leaching and surface runoff. Sorption is one of the most important phenomena that affects leaching, and thus the fate of hydrophobic organic pollutants in soils and also control their distribution in the soil/water environment. The work developed focuses the optimization of analytical techniques for monitoring the sorption behaviour of organic pollutants, 17α- ethinylestradiol (EE2) and atrazine, and their fate in aqueous environment. Initially, the development of several analytical techniques, such as micellar electrokinetic chromatography, spectral deconvolution, using UV-Vis and fluorescence spectroscopy, and also enzyme linked immunosorbent assay was performed. Optimization, method performance and recovery tests are described and results discussed. Moreover, in order to evaluate the applicability of the previously optimized method, atrazine and EE2 sorption to soil samples was performed. The work developed provide several options, in terms of methodology to follow sorption of atrazine onto soils, however the choice depends on the laboratory conditions and on the analyst preferences. The advantages and disadvantages of each methodology should be evaluated first. The second part of this work consisted in the sorption behaviour study of those two different hydrophobic organic pollutants onto different soil samples. Soil organic matter chemical characterization, being essential to understand the binding mechanism responsible for the interactions, was made. The results of atrazine binding to organic matter pointed out that carboxyl units and aromaticrich organic matter are the most efficient binding agents for atrazine. EE2 adsorbs strongly to soil organic matter and is mainly stabilized by hydrophobic interactions, through aromatic nuclei face to face with surface and/or another EE2 molecule association. Farmyard manure soil contains higher aromatic and carboxyl units, indicating that this type of manure can be effectively used to minimize the residual toxicity of EE2 and atrazine present in soils, increasing the sorption and reducing leaching onto water resources. Since the final destination of organic pollutants can be ground, surface and/or waste water, atrazine and 17α-ethinylestradiol were quantified in several water samples.

Computational methods in the fractional calculus of variations and optimal control

The fractional calculus of variations and fractional optimal control are generalizations of the corresponding classical theories, that allow problem modeling and formulations with arbitrary order derivatives and integrals. Because of the lack of analytic methods to solve such fractional problems, numerical techniques are developed. Here, we mainly investigate the approximation of fractional operators by means of series of integer-order derivatives and generalized finite differences. We give upper bounds for the error of proposed approximations and study their efficiency. Direct and indirect methods in solving fractional variational problems are studied in detail. Furthermore, optimality conditions are discussed for different types of unconstrained and constrained variational problems and for fractional optimal control problems. The introduced numerical methods are employed to solve some illustrative examples.

Phase transitions and nonlinear phenomena in neuronal network models

Communication and cooperation between billions of neurons underlie the power of the brain. How do complex functions of the brain arise from its cellular constituents? How do groups of neurons self-organize into patterns of activity? These are crucial questions in neuroscience. In order to answer them, it is necessary to have solid theoretical understanding of how single neurons communicate at the microscopic level, and how cooperative activity emerges. In this thesis we aim to understand how complex collective phenomena can arise in a simple model of neuronal networks. We use a model with balanced excitation and inhibition and complex network architecture, and we develop analytical and numerical methods for describing its neuronal dynamics. We study how interaction between neurons generates various collective phenomena, such as spontaneous appearance of network oscillations and seizures, and early warnings of these transitions in neuronal networks. Within our model, we show that phase transitions separate various dynamical regimes, and we investigate the corresponding bifurcations and critical phenomena. It permits us to suggest a qualitative explanation of the Berger effect, and to investigate phenomena such as avalanches, band-pass filter, and stochastic resonance. The role of modular structure in the detection of weak signals is also discussed. Moreover, we find nonlinear excitations that can describe paroxysmal spikes observed in electroencephalograms from epileptic brains. It allows us to propose a method to predict epileptic seizures. Memory and learning are key functions of the brain. There are evidences that these processes result from dynamical changes in the structure of the brain. At the microscopic level, synaptic connections are plastic and are modified according to the dynamics of neurons. Thus, we generalize our cortical model to take into account synaptic plasticity and we show that the repertoire of dynamical regimes becomes richer. In particular, we find mixed-mode oscillations and a chaotic regime in neuronal network dynamics.

Concentration of tumor biomarkers using aqueous biphasic systems

According to the World Health Organization, around 8.2 million people die each year with cancer. Most patients do not perform routine diagnoses and the symptoms, in most situations, occur when the patient is already at an advanced stage of the disease, consequently resulting in a high cancer mortality. Currently, prostate cancer is the second leading cause of death among males worldwide. In Portugal, this is the most diagnosed type of cancer and the third that causes more deaths. Taking into account that there is no cure for advanced stages of prostate cancer, the main strategy comprises an early diagnosis to increase the successful rate of the treatment. The prostate specific antigen (PSA) is an important biomarker of prostate cancer that can be detected in biological fluids, including blood, urine and semen. However, the commercial kits available are addressed for blood samples and the commonly used analytical methods for their detection and quantification requires specialized staff, specific equipment and extensive sample processing, resulting in an expensive process. Thus, the aim of this MSc thesis consisted on the development of a simple, efficient and less expensive method for the extraction and concentration of PSA from urine samples using aqueous biphasic systems (ABS) composed of ionic liquids. Initially, the phase diagrams of a set of aqueous biphasic systems composed of an organic salt and ionic liquids were determined. Then, their ability to extract PSA was ascertained. The obtained results reveal that in the tested systems the prostate specific antigen is completely extracted to the ionic-liquid-rich phase in a single step. Subsequently, the applicability of the investigated ABS for the concentration of PSA was addressed, either from aqueous solutions or urine samples. The low concentration of this biomarker in urine (clinically significant below 150 ng/mL) usually hinders its detection by conventional analytical techniques. The obtained results showed that it is possible to extract and concentrate PSA, up to 250 times in a single-step, so that it can be identified and quantified using less expensive techniques.

Concentration of tumor biomarkers using aqueous biphasic systems

According to the World Health Organization, around 8.2 million people die each year with cancer. Most patients do not perform routine diagnoses and the symptoms, in most situations, occur when the patient is already at an advanced stage of the disease, consequently resulting in a high cancer mortality. Currently, prostate cancer is the second leading cause of death among males worldwide. In Portugal, this is the most diagnosed type of cancer and the third that causes more deaths. Taking into account that there is no cure for advanced stages of prostate cancer, the main strategy comprises an early diagnosis to increase the successful rate of the treatment. The prostate specific antigen (PSA) is an important biomarker of prostate cancer that can be detected in biological fluids, including blood, urine and semen. However, the commercial kits available are addressed for blood samples and the commonly used analytical methods for their detection and quantification requires specialized staff, specific equipment and extensive sample processing, resulting in an expensive process. Thus, the aim of this MSc thesis consisted on the development of a simple, efficient and less expensive method for the extraction and concentration of PSA from urine samples using aqueous biphasic systems (ABS) composed of ionic liquids. Initially, the phase diagrams of a set of aqueous biphasic systems composed of an organic salt and ionic liquids were determined. Then, their ability to extract PSA was ascertained. The obtained results reveal that in the tested systems the prostate specific antigen is completely extracted to the ionic-liquid-rich phase in a single step. Subsequently, the applicability of the investigated ABS for the concentration of PSA was addressed, either from aqueous solutions or urine samples. The low concentration of this biomarker in urine (clinically significant below 150 ng/mL) usually hinders its detection by conventional analytical techniques. The obtained results showed that it is possible to extract and concentrate PSA, up to 250 times in a single-step, so that it can be identified and quantified using less expensive techniques.

Population genetics of C. Maenas : oceanography and larval dispersal

Decifrar a complexa interacção entre os ciclos de vida de espécies marinhas e a oceanografia revela-se fundamental para a compreensão do fluxo genético e da conectividade no meio marinho. Nas espécies marinhas com desenvolvimento indirecto o fluxo de genes entre populações depende da distância que separa as populações, bem como da interacção entre a duração do desenvolvimento larvar, do comportamento das larvas e dos padrões de circulação oceânica. A conectividade larvar influencia uma variedade de processos como a dinâmica de stocks e de populações, a distribuição e limites geográficos das espécies, a estrutura genética das populações e a dispersão de espécies invasivas e reveste-se consequentemente de uma importância fundamental na identificação das unidades populacionais evolucionariamente relevantes e para a gestão e conservação marinhas. Os marcadores genéticos e os Modelos Individuais Acoplados a Modelos Físico-Biológicos (“ICPBMs”) são actualmente ferramentas fundamentais para o estudo dos padrões de dispersão larvar e para avaliar o nível de conectividade populacional. A presente tese respeita à avaliação das escalas espaciais de conectividade de populações de uma espécie costeira, o caranguejo Carcinus maenas, e utiliza conjuntamente informação de marcadores genéticos, análise de séries temporais de fornecimento de larvas e um modelo numérico de circulação oceânica. O primeiro capítulo introduz a temática da conectividade em espécies marinhas e inclui algumas referências aos métodos moleculares, analíticos e de modelação seguidos ao longo da tese. Através da utilização de múltiplas ferramentas – avaliação da estrutura genética geográfica de C. maenas na sua distribuição nativa com recurso a marcadores de DNA (microssatélites) (Capítulo 2), avaliação da estrutura genética temporal das larvas que formam os eventos de fornecimento larvar à Ria de Aveiro, NW Portugal (Capítulo 3), descrição da variabilidade inter-anual do fornecimento larvar à Ria de Aveiro, NW Portugal (Capítulo 4) e validação de um modelo ICPBM que descreve os padrões observados de fornecimento (Capítulo 5) – esta tese espera poder contribuir para uma melhor compreensão dos mecanismos que regulam o fluxo de genes e a conectividade entre populações de organismos marinhos. No Capítulo 6 são apresentadas as principais conclusões da investigação. A análise genética com recurso a microssatélites indicou que as populações de C. maenas são geneticamente homogéneas ao longo de várias centenas de km, dentro da distribuição nativa da espécie. Paralelamente, não foram encontrados indícios da existência de reprodução por “sweepstakes” em C. maenas de populações da costa oeste da Península Ibérica, visto que não se obtiveram diferenças genéticas significativas entre os eventos larvares. Também não se encontrou qualquer estrutura familiar entre as larvas que formam cada episódio de fornecimento, e não houve nenhuma redução significativa da variabilidade genética das larvas quando comparada com a de caranguejos adultos. A análise de séries temporais de suprimento de larvas na Ria de Aveiro em cinco anos estudados indica que este é um fenómeno episódico e variável, sendo os maiores episódios de fornecimento coincidentes com as marés vivas e acentuados por fortes ventos de sul. O modelo ICPBM foi validado com sucesso e parece fornecer uma estimativa realística das escalas espaciais e temporais de dispersão larvar, de acordo com as observações da estrutura genética e da ausência de reprodução por “sweepstake” em C. maenas da costa oeste da Península Ibérica

Metabonomics of human amniotic fluid for prenatal diagnostics

O trabalho apresentado nesta tese teve como principais objectivos contribuir para o conhecimento da composição do líquido amniótico humano (LA), colhido no 2º trimestre de gravidez, assim como investigar possíveis alterações na sua composição devido à ocorrência de patologias pré-natais, recorrendo à metabonómica e procurando, assim, definir novos biomarcadores de doenças da grávida e do feto. Após uma introdução descrevendo o estado da arte relacionado com este trabalho (Capítulo 1) e os princípios das metodologias analíticas usadas (Capítulo 2), seguida de uma descrição dos aspectos experimentais associados a esta tese (Capítulo 3), apresentam-se os resultados da caracterização da composição química do LA (gravidez saudável) por espectroscopia de ressonância magnética nuclear (RMN), assim como da monitorização da sua estabilidade durante o armazenamento e após ciclos de congelamento-descongelamento (Capítulo 4). Amostras de LA armazenadas a -20°C registaram alterações significativas, tornando-se estas menos pronunciadas (mas ainda mensuráveis) a -70°C, temperatura recomendada para o armazenamento de LA. Foram também observadas alterações de composição após 1-2 ciclos de congelamento-descongelamento (a ter em conta aquando da reutilização de amostras), assim como à temperatura ambiente (indicando um período máximo de 4h para a manipulação e análise de LA). A aquisição de espectros de RMN de 1H de alta resolução e RMN acoplado (LC-NMR/MS) permitiu a detecção de 75 compostos no LA do 2º trimestre, 6 dos quais detectados pela primeira vez no LA. Experiências de difusão (DOSY) permitiram ainda a caracterização das velocidades de difusão e massas moleculares médias das proteínas mais abundantes. O Capítulo 5 descreve o estudo dos efeitos de malformações fetais (FM) e de cromossomopatias (CD) na composição do LA do 2º trimestre de gravidez. A extensão deste trabalho ao estudo dos efeitos de patologias no LA que ocorrem no 3º trimestre de gravidez é descrita no Capítulo 6, nomeadamente no que se refere ao parto pré-termo (PTD), pré-eclampsia (PE), restrição do crescimento intra-uterino (IUGR), ruptura prematura de membranas (PROM) e diabetes mellitus gestacional (GDM). Como complemento a estes estudos, realizou-se uma análise preliminar da urina materna do 2º trimestre para o estudo de FM e GDM, descrita no Capítulo 7. Para interpretação dos dados analíticos, obtidos por espectroscopia RMN de 1H, cromatografia líquida de ultra eficiência acoplada a espectrometria de massa (UPLC-MS) e espectroscopia do infravermelho médio (MIR), recorreu-se à análise discriminante pelos métodos dos mínimos quadrados parciais e o método dos mínimos quadrados parciais ortogonal (PLS-DA e OPLS-DA) e à correlação espectral. Após análise por validação cruzada de Monte-Carlo (MCCV), os modelos PLS-DA de LA permitiram distinguir as FM dos controlos (sensibilidades 69-85%, especificidades 80-95%, taxas de classificação 80-90%), revelando variações metabólicas ao nível do metabolismo energético, dos metabolismos dos aminoácidos e glícidos assim como possíveis alterações ao nível do funcionamento renal. Observou-se também um grande impacto das FM no perfil metabólico da urina materna (medido por UPLC-MS), tendo no entanto sido registados modelos PLS-DA com menor sensibilidade (40-60%), provavelmente devido ao baixo número de amostras e maior variabilidade da composição da urina (relativamente ao LA). Foram sugeridos possíveis marcadores relacionados com a ocorrência de FM, incluindo lactato, glucose, leucina, valina, glutamina, glutamato, glicoproteínas e conjugados de ácido glucurónico e/ou sulfato e compostos endógenos e/ou exógenos (<1 M) (os últimos visíveis apenas na urina). No LA foram também observadas variações metabólicas devido à ocorrência de vários tipos de cromossomopatias (CD), mas de menor magnitude. Os perfis metabólicos de LA associado a pré- PTD produziram modelos que, apesar do baixo poder de previsão, sugeriram alterações precoces no funcionamento da unidade fetoplacentária, hiperglicémia e stress oxidativo. Os modelos obtidos para os grupos pré- IUGR pré- PE, pré- PROM e pré-diagnóstico GDM (LA e urina materna) registaram baixo poder de previsão, indicando o pouco impacto destas condições na composição do LA e/ou urina do 2º trimestre. Os resultados obtidos demonstram as potencialidades da análise dos perfis metabólicos do LA (e, embora com base em menos estudos, da urina materna) do 2º trimestre para o desenvolvimento de novos e complementares métodos de diagnóstico, nomeadamente para FM e PTD.

Challenges in mercury speciation and fractionation in soil and sediment

This investigation focused on the development, test and validation of methodologies for mercury fractionation and speciation in soil and sediment. After an exhaustive review of the literature, several methods were chosen and tested in well characterised soil and sediment samples. Sequential extraction procedures that divide mercury fractions according to their mobility and potential availability in the environment were investigated. The efficiency of different solvents for fractionation of mercury was evaluated, as well as the adequacy of different analytical instruments for quantification of mercury in the extracts. Kinetic experiments to establish the equilibrium time for mercury release from soil or sediment were also performed. It was found that in the studied areas, only a very small percentage of mercury is present as mobile species and that mobility is associated to higher aluminium and manganese contents, and that high contents of organic matter and sulfur result in mercury tightly bound to the matrix. Sandy soils tend to release mercury faster that clayey soils, and therefore, texture of soil or sediment has a strong influence on the mobility of mercury. It was also understood that analytical techniques for quantification of mercury need to be further developed, with lower quantification limits, particularly for mercury quantification of less concentrated fractions: water-soluble e exchangeable. Although the results provided a better understanding of the distribution of mercury in the sample, the complexity of the procedure limits its applicability and robustness. A proficiency-testing scheme targeting total mercury determination in soil, sediment, fish and human hair was organised in order to evaluate the consistency of results obtained by different laboratories, applying their routine methods to the same test samples. Additionally, single extractions by 1 mol L-1 ammonium acetate solution, 0.1 mol L-1 HCl and 0.1 mol L-1 CaCl2, as well as extraction of the organometallic fraction were proposed for soil; the last was also suggested for sediment and fish. This study was important to update the knowledge on analytical techniques that are being used for mercury quantification, the associated problems and sources of error, and to improve and standardize mercury extraction techniques, as well as to implement effective strategies for quality control in mercury determination. A different, “non chemical-like” method for mercury species identification was developed, optimised and validated, based on the thermo-desorption of the different mercury species. Compared to conventional extraction procedures, this method has advantages: it requires little to no sample treatment; a complete identification of species present is obtained in less than two hours; mercury losses are almost neglectable; can be considered “clean”, as no residues are produced; the worldwide comparison of results obtained is easier and reliable, an important step towards the validation of the method. Therefore, the main deliverables of this PhD thesis are an improved knowledge on analytical procedures for identification and quantification of mercury species in soils and sediments, as well as a better understanding of the factors controlling the behaviour of mercury in these matrices.

Quantificação de mercúrio em águas de consumo: validação do método na Luságua

Para o ser humano, a água sempre foi um recurso essencial ao longo da sua evolução. Hoje em dia, com todo o avanço tecnológico, a água é um bem muito vulnerável às diversas atividades antropogénicas. Entre os muitos contaminantes que podem afetar a qualidade da água para consumo humano, os metais causam grande preocupação devido à sua elevada toxicidade. O mercúrio é um dos contaminantes que deve ser devidamente controlado devido ao seu elevado grau de toxicidade. Este Estágio foi realizado no Laboratório Luságua e o seu principal objetivo foi avaliar se era possível baixar o limite de quantificação da técnica existente na empresa para quantificar o mercúrio em águas de consumo e validar o método através do cálculo de vários parâmetros de controlo de qualidade. Foi ainda feita a comparação entre a técnica de quantificação de mercúrio existente na Luságua (CV-AAS) e a existente na Universidade de Aveiro (CVAFS) para avaliar se havia uma mais-valia para a Luságua se adquirisse um novo equipamento para analisar o mercúrio em águas. Os parâmetros de validação avaliados nos dois métodos derem resultados semelhantes, não identificando a necessidade atual de substituir o equipamento existente na Luságua, até porque se conseguiu baixar o limite de quantificação, atingindo assim o objetivo estabelecido no início para este Estágio.

Determinant factors for the use of imposex in TBT pollution monitoring

Os compostos de tributilestanho (TBT) foram utilizados como biocidas em tintas antivegetativas (AV) e amplamente aplicados, durante décadas, de forma a evitar a bioincrustação em superfícies submersas, principalmente em cascos de embarcações. Os seus efeitos deletérios em organismos não-alvo tornaram-se evidentes após o aparecimento de gastrópodes prosobrâquios com imposex – sobreposição de caracteres sexuais masculinos sobre o tracto reprodutivo das fêmeas. Desde então, a expressão do imposex em prosobrânquios tem sido amplamente utilizada como biomarcador da poluição por TBT. Um dos objectivos da presente tese é avaliar se as mais recentes restrições legais na utilização das tintas AV com compostos organoestânicos (OTs) resultaram numa redução da poluição pelos mesmos na costa continental Portuguesa. Para tal, foi levada a cabo a análise da variação temporal do imposex como biomarcador dos níveis ambientais de TBT e a validação de procedimentos, de modo a seguir de forma precisa a evolução da intensidade deste fenómeno ao longo do tempo. O trabalho de investigação teve início no momento em que a ineficácia da legislação anterior (Directiva 89/677/CEE) na redução da poluição por TBT no litoral Português era reportada na literatura e quando estavam já agendados instrumentos decisivos para a diminuição definitiva deste tipo de poluição: o Regulamento (CE) N.º 782/2003 do Parlamento e do Conselho da União Europeia (UE) que bania as tintas AV baseadas em TBT na sua frota a partir de 1 de Julho de 2003; a "Convenção Internacional sobre o Controlo de Sistemas Antivegetativos Nocivos em Navios” (Convenção AFS), adoptada em 2001 pela Organização Marítima Internacional (OMI), que procurava erradicar os OTs da frota mundial até 2008. Os níveis de imposex e as concentrações de OTs nos tecidos de fêmeas de Nucella lapillus foram medidos em 17 locais de amostragem ao longo da costa Portuguesa em 2003, a fim de se avaliarem os impactos do TBT nas populações desta espécie e de se criar uma base de dados para seguir a sua evolução futura. O índice da sequência do vaso deferente (VDSI), o índice do tamanho relativo do pénis (RPSI) e a percentagem de fêmeas afectadas por imposex (%I) foram utilizados na medição da intensidade deste fenómeno em cada local de amostragem e os seus valores variaram entre 0,20-4,04, 0,0- 42,2% e 16,7-100,0%, respectivamente. Foram encontradas fêmeas estéreis em 3 locais de amostragem, com percentagens a variar entre 4,0 e 6,2%. As concentrações de TBT e dibutilestanho (DBT) nas fêmeas variaram respectivamente entre 23-138 e <10-62 ng Sn.g-1 de peso seco, e o conteúdo em TBT nos tecidos revelou-se significativamente correlacionado com o imposex (nomeadamente com o RPSI e o VDSI). Os níveis de expressão do fenómeno e o conteúdo em OTs nos tecidos, foram superiores na proximidade de portos, confirmando as conclusões obtidas anteriormente por outros autores de que os navios e a actividade dos estaleiros são as principais fontes destes compostos no litoral Português. As infra-estruturas associadas ás principais fontes de OTs – portos, estaleiros e marinas – estão geralmente localizadas no interior de estuários, motivo pelo qual estas áreas têm vindo a ser descritas como as mais afectadas por estes poluentes. Por esta razão, foi levado a cabo o estudo pormenorizado da poluição por TBT na Ria de Aveiro, como caso de estudo representativo da poluição por estes compostos num sistema estuarino em Portugal continental. N. lapillus foi usada como bioindicador para avaliar a tendência temporal da poluição por TBT nesta área entre 1997 e 2007. Foi registada uma diminuição da intensidade do imposex após 2003, embora as melhorias mais evidentes tenham sido observadas entre 2005 e 2007, provavelmente devido à implementação do Regulamento (CE) N.º 782/2003 que proibiu a aplicação de tintas AV com TBT em navios com a bandeira da UE. Apesar desses progressos, as análises ao conteúdo em OTs nos tecidos de fêmeas de N. lapillus e em amostras de água colhidas em 2006 indicaram contaminação recente por TBT na área de estudo, evidenciando assim a permanência de fontes de poluição. A utilização de N. lapillus como bioindicador da poluição por TBT na Ria de Aveiro apresenta algumas limitações uma vez que a espécie não ocorre nas áreas mais interiores da Ria e não vive em contacto com sedimentos. Assim, a informação obtida a partir da sua utilização como bioindicador é fundamentalmente relativa aos níveis de TBT na coluna de água. Foi então necessário recorrer a um bioindicador complementar – Hydrobia ulvae – para melhor avaliar a evolução temporal da poluição por TBT no interior deste sistema estuarino e estudar a persistência de TBT nos sedimentos. Não foi registada uma diminuição dos níveis de imposex em H. ulvae na Ria de Aveiro entre 1998 e 2007, apesar da aplicação do Regulamento (CE) N.º 782/2003. Pelo contrário, houve um aumento global significativo da percentagem de fêmeas afectadas por imposex e um ligeiro aumento do VDSI, contrastando com o que tem sido descrito para outros bioindicadores na Ria de Aveiro no mesmo período. Estes resultados mostram que a diferente biologia/ecologia das espécies indicadoras determina vias distintas de acumulação de TBT, apontando a importância da escolha do bioindicador dependendo do compartimento a ser monitorizado (sedimento versus água). A ingestão de sedimento como hábito alimentar em H. ulvae foi discutida como sendo a razão para a escolha da espécie como indicadora da contaminação dos sedimentos por TBT. Foram também estudados os métodos mais fiáveis para reduzir a influência de variáveis críticas na medição dos níveis de imposex em H. ulvae. As comparações de parâmetros do imposex baseados em medições do pénis devem ser sempre realizadas sob condições de narcotização bem standardizadas uma vez que este procedimento provoca um aumento significativo do comprimento do pénis (PL) em ambos os sexos. O VDSI, a % e o PL em ambos os sexos revelaram ser fortemente influenciados pelo tamanho dos espécimes: a utilização de fêmeas mais pequenas conduz à subestimação do VDSI, da %I e do PL, enquanto que diferenças no tamanho dos machos provocam variações no índice do comprimento relativo do pénis (RPLI), independentemente dos níveis de poluição por TBT. Existe, portanto, a necessidade de controlar algumas variáveis envolvidas na análise do imposex que mostraram afectar a fiabilidade dos resultados. Uma vez que N. lapillus é o principal bioindicador dos efeitos biológicos específicos do TBT para a área da OSPAR, foi também estudada a influência de algumas variáveis na avaliação dos níveis de imposex nesta espécie, especificamente as relacionadas com o ciclo reprodutor e o tamanho dos espécimes. O estudo do ciclo reprodutor e a variação sazonal/espacial do comprimento do pénis do macho (MPL) incidiu num único local no litoral Português (Areão – região de Aveiro) de forma a avaliar se o RPSI varia sazonal e espacialmente na mesma estação de amostragem e se tais variações influenciam os resultados obtidos em programas de monitorização do imposex. Nos meses de Dezembro de 2005 a Junho de 2007, foram encontrados espécimes de N. lapillus sexualmente maturos e potencialmente aptos para a reprodução. Contudo, foi também evidente um padrão sazonal do ciclo reprodutor – o estado de desenvolvimento da gametogénese nas fêmeas variou sazonalmente e ocorreu uma diminuição do volume da glândula da cápsula e do factor de condição no final do Verão / início do Outono. Contrariamente, a gametogénese nos machos não apresentou variação sazonal significativa, embora os valores mais baixos do factor de condição, do comprimento do pénis e dos volumes de esperma e da próstata tenham também sido registados no final do Verão / início do Outono. Além disso, o MPL mostrou variar, no mesmo local de amostragem, inversamente com a distância aos ninhos de cápsulas ovígeras; um aumento dos valores do MPL foi também observado em espécimes de maior tamanho. Todas estas variações no MPL introduzem desvios nos resultados da avaliação do imposex quando é usado o RPSI. A magnitude do erro envolvido foi quantificada e revelou-se superior em locais com níveis mais elevados de poluição por TBT. Apesar do RPSI ser um índice que fornece informação importante sobre os níveis de poluição por TBT, a sua interpretação deve ser cuidadosa e realizada complementarmente com os outros índices, destacando-se o VDSI como índice mais fiável e com significado biológico mais expressivo. Novas campanhas de monitorização do imposex em N. lapillus foram realizadas ao longo da costa Portuguesa em 2006 e 2008, e os resultados subsequentemente comparados com a base de dados criada em 2003,de forma a avaliar a evolução da poluição por TBT no litoral Português naquele período. Nestes estudos foram aplicados novos procedimentos na monitorização e tratamento dos dados, de forma a minimizar a variação nos índices de avaliação do imposex induzida pelos factores acima descritos, para seguir com maior consistência a tendência da poluição por TBT entre 2003 e 2008. Foi observado um declínio global significativo na intensidade de imposex na área de estudo durante o referido período e a qualidade ecológica da costa Portuguesa, segundo os termos definidos pela Comissão OSPAR, revelou melhorias notáveis após 2003 confirmando a eficácia do Regulamento (CE) N.º 782/2003. Não obstante, as populações de N. lapillus revelaram-se ainda amplamente afectadas por imposex, tendo sido detectadas emissões de TBT na água do mar ao longo da costa em 2006, apesar da restrição anteriormente referida. Estes inputs foram atribuídos principalmente aos navios que à data ainda circulavam com tintas AV com TBT aplicadas antes de 2003, uma vez que a sua utilização nas embarcações foi apenas proibida em 2008. Considerando que o Regulamento (CE) N.º 782/2003 constitui uma antecipação da proibição global da OMI que entrou em vigor em Setembro de 2008, prevê-se, por analogia, que haja uma rápida diminuição da poluição por TBT à escala mundial num futuro próximo. Na sequência desta previsão, é apresentada uma discussão teórica preliminar relativamente ás possíveis estratégias usadas por N. lapillus na recolonização de locais onde, no passado, as populações terão aparentemente sido extintas devido a níveis extremamente elevados de TBT. Estes locais são tipicamente zonas abrigadas junto de infra-estruturas portuárias, cuja recolonização por esta espécie será provavelmente muito lenta dada a mobilidade reduzida dos adultos e o ciclo de vida não apresentar fase larvar pelágica. Foram então equacionadas duas hipotéticas vias de recolonização de zonas abrigadas por espécimes provenientes de populações de costa aberta/exposta: a migração de adultos e/ou a dispersão de juvenis. No entanto, em ambos os casos, estaria implicada a transposição de um problema amplamente descrito na bibliografia: as populações de N. lapillus de costa aberta podem apresentar um fenótipo muito diferente das de zonas abrigadas, podendo inclusivamente variar no número de cromossomas. A recolonização pode portanto não ter sucesso pelo simples facto dos novos recrutas não estarem bem adaptados aos locais a recolonizar. Para analisar este problema, foram estudados tanto a forma da concha de espécimes de N. lapillus ao longo da costa Portuguesa como o respectivo número de cromossomas. Embora a forma da concha tenha revelado diferenças, de acordo com o grau de hidrodinamismo entre as populações de N. lapillus avaliadas, o cariótipo 2n = 26, típico de zonas expostas, foi registado em todos os locais amostrados. Por outro lado, foi também testada em laboratório a possível mortalidade de juvenis em dispersão no interior menos salino dos estuários. Foi então verificada a ocorrência de mortalidade elevada de juvenis expostos a salinidades baixas (=100% após 1 hora a salinidades ≤9 psu), o que também pode comprometer a recolonização de zonas estuarinas menos salinas por esta via. Mesmo assim, os juvenis mostraram um comportamento de flutuação à superfície da água em condições laboratoriais, que pode ser considerado um benefício específico na colonização de áreas mais internas dos estuários, se tal facto vier a ser confirmado em estudos de campo. As conclusões deste estudo contribuem certamente para a descrição do final da “história do TBT” dado que, uma vez controlados alguns factores determinantes no uso do imposex como biomarcador, a avaliação do declínio da poluição por estes compostos, agora esperado à escala global, se torna mais rigorosa.

Métodos robustos em geoestatística

O objectivo principal da presente tese consiste no desenvolvimento de estimadores robustos do variograma com boas propriedades de eficiência. O variograma é um instrumento fundamental em Geoestatística, pois modela a estrutura de dependência do processo em estudo e influencia decisivamente a predição de novas observações. Os métodos tradicionais de estimação do variograma não são robustos, ou seja, são sensíveis a pequenos desvios das hipóteses do modelo. Essa questão é importante, pois as propriedades que motivam a aplicação de tais métodos, podem não ser válidas nas vizinhanças do modelo assumido. O presente trabalho começa por conter uma revisão dos principais conceitos em Geoestatística e da estimação tradicional do variograma. De seguida, resumem-se algumas noções fundamentais sobre robustez estatística. No seguimento, apresenta-se um novo método de estimação do variograma que se designou por estimador de múltiplos variogramas. O método consiste em quatro etapas, nas quais prevalecem, alternadamente, os critérios de robustez ou de eficiência. A partir da amostra inicial, são calculadas, de forma robusta, algumas estimativas pontuais do variograma; com base nessas estimativas pontuais, são estimados os parâmetros do modelo pelo método dos mínimos quadrados; as duas fases anteriores são repetidas, criando um conjunto de múltiplas estimativas da função variograma; por fim, a estimativa final do variograma é definida pela mediana das estimativas obtidas anteriormente. Assim, é possível obter um estimador que tem boas propriedades de robustez e boa eficiência em processos Gaussianos. A investigação desenvolvida revelou que, quando se usam estimativas discretas na primeira fase da estimação do variograma, existem situações onde a identificabilidade dos parâmetros não está assegurada. Para os modelos de variograma mais comuns, foi possível estabelecer condições, pouco restritivas, que garantem a unicidade de solução na estimação do variograma. A estimação do variograma supõe sempre a estacionaridade da média do processo. Como é importante que existam procedimentos objectivos para avaliar tal condição, neste trabalho sugere-se um teste para validar essa hipótese. A estatística do teste é um estimador-MM, cuja distribuição é desconhecida nas condições de dependência assumidas. Tendo em vista a sua aproximação, apresenta-se uma versão do método bootstrap adequada ao estudo de observações dependentes de processos espaciais. Finalmente, o estimador de múltiplos variogramas é avaliado em termos da sua aplicação prática. O trabalho contém um estudo de simulação que confirma as propriedades estabelecidas. Em todos os casos analisados, o estimador de múltiplos variogramas produziu melhores resultados do que as alternativas usuais, tanto para a distribuição assumida, como para distribuições contaminadas.