Constraint-based modeling of minimum set covering: application to species differentation

Work presented in the context of the European Master in Computational Logics, as partial requisit for the graduation as Master in Computational Logics

Design and development of a microfluidic platform for use with colorimetric gold nanoprobe assays

Due to the importance and wide applications of the DNA analysis, there is a need to make genetic analysis more available and more affordable. As such, the aim of this PhD thesis is to optimize a colorimetric DNA biosensor based on gold nanoprobes developed in CEMOP by reducing its price and the needed volume of solution without compromising the device sensitivity and reliability, towards the point of care use. Firstly, the price of the biosensor was decreased by replacing the silicon photodetector by a low cost, solution processed TiO2 photodetector. To further reduce the photodetector price, a novel fabrication method was developed: a cost-effective inkjet printing technology that enabled to increase TiO2 surface area. Secondly, the DNA biosensor was optimized by means of microfluidics that offer advantages of miniaturization, much lower sample/reagents consumption, enhanced system performance and functionality by integrating different components. In the developed microfluidic platform, the optical path length was extended by detecting along the channel and the light was transmitted by optical fibres enabling to guide the light very close to the analysed solution. Microfluidic chip of high aspect ratio (~13), smooth and nearly vertical sidewalls was fabricated in PDMS using a SU-8 mould for patterning. The platform coupled to the gold nanoprobe assay enabled detection of Mycobacterium tuberculosis using 3 8l on DNA solution, i.e. 20 times less than in the previous state-of-the-art. Subsequently, the bio-microfluidic platform was optimized in terms of cost, electrical signal processing and sensitivity to colour variation, yielding 160% improvement of colorimetric AuNPs analysis. Planar microlenses were incorporated to converge light into the sample and then to the output fibre core increasing 6 times the signal-to-losses ratio. The optimized platform enabled detection of single nucleotide polymorphism related with obesity risk (FTO) using target DNA concentration below the limit of detection of the conventionally used microplate reader (i.e. 15 ng/μl) with 10 times lower solution volume (3 μl). The combination of the unique optical properties of gold nanoprobes with microfluidic platform resulted in sensitive and accurate sensor for single nucleotide polymorphism detection operating using small volumes of solutions and without the need for substrate functionalization or sophisticated instrumentation. Simultaneously, to enable on chip reagents mixing, a PDMS micromixer was developed and optimized for the highest efficiency, low pressure drop and short mixing length. The optimized device shows 80% of mixing efficiency at Re = 0.1 in 2.5 mm long mixer with the pressure drop of 6 Pa, satisfying requirements for the application in the microfluidic platform for DNA analysis.

Behavioral and synaptic circuit analysis in models of neuropsychiatric disorders: Dissecting the in vivo role of the postsynaptic density proteins nArgBP2 and Shank3 using genetic engineered mice

Dissertation presented to obtain the Ph.D degree in Biology

Genetic and functional analysis of the bacillus subtilis BSP1 gam cluster

Mannans (linear mannan, glucomannan, galactomannan and galactoglucomannan) are the major constituents of the hemicellulose fraction in softwoods and show great importance as a renewable resource for fuel or feedstock applications. As complex polysaccharides, mannans can only be degraded through a synergistic action of different mannan-degrading enzymes, mannanases. Microbial mannanases are mainly extracellular enzymes that can act in wide range of pH and temperature, contributing to pulp and paper, pharmaceutical, food and feed, oil and textile successful industrial applications. Knowing and controlling these microbial mannan-degrading enzymes are essential to take advantage of their great biotechnological potential. The genome of the laboratory 168 strain of Bacillus subtilis carries genes gmuA-G dedicated to the degradation and utilization of glucomannan, including an extracellular -mannanase. Recently, the genome sequence of an undomesticated strain of B. subtilis, BSP1, was determined. In BSP1, the gmuA-G operon is maintained, interestingly, however, a second cluster of genes was found (gam cluster), which comprise a second putative extracellular β-mannanase, and most likely specify a system for the degradation and utilization of a different mannan polymer, galactoglucomannan. The genetic organization and function of the gam cluster, and whether its presence in BSP1 strain results in new hemicellulolytic capabilities, compared to those of the laboratory strain, was address in this work. In silico and in vivo mRNA analyses performed in this study revealed that the gam cluster, comprising nine genes, is organized and expressed in at least six different transcriptional units. Furthermore, cloning, expression, and production of Bbsp2923 in Escherichia coli was achieved and preliminary characterization shows that the enzyme is indeed a β-mannanase. Finally, the high hemicellulolytic capacity of the undomesticated B. subtilis BSP1, demonstrated in this work by qualitative analyses, suggests potential to be used in the food and feed industries.

Social network analysis applied to knowledge creation and transfer in the portuguese agricultural R&D field: an exploratory study

Paper presented at the 9th European Conference on Knowledge Management, Southampton Solent University, Southampton, UK, 4-5 Sep. 2008. URL: http://academic-conferences.org/eckm/eckm2008/eckm08-home.htm

Estrogen receptor dependent genetic and epigenetic factors of tamoxifen resistance

Resumo: A decisão da terapêutica hormonal no tratamento do cancro da mama baseiase na determinação do receptor de estrogénio alfa por imunohistoquímica (IHC). Contudo, a presença deste receptor não prediz a resposta em todas as situações, em parte devido a limitações do método IHC. Investigámos se a expressão dos genes ESR1 e ESR2, bem como a metilação dos respectivos promotores, pode estar relacionada com a evolução desfavorável de uma proporção de doentes tratados com tamoxifeno assim como com a perda dos receptores de estrogénio alfa (ERα) e beta (ERß). Amostras de 211 doentes com cancro da mama diagnosticado entre 1988 e 2004, fixadas em formalina e preservadas em parafina, foram utilizadas para a determinação por IHC da presença dos receptores ERα e ERß. O mRNA total do gene ESR1 e os níveis específicos do transcrito derivado do promotor C (ESR1_C), bem como dos transcritos ESR2_ß1, ESR2_ß2/cx, and ESR2_ß5 foram avaliados por Real-time PCR. Os promotores A e C do gene ESR1 e os promotores 0K e 0N do gene ESR2 foram investigados por análise de metilação dos dinucleotidos CpG usando bisulfite-PCR para análise com enzimas de restrição, ou para methylation specific PCR. Atendendo aos resultados promissores relacionados com a metilação do promotor do gene ESR1, complementamos o estudo com um método quantitativo por matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) suportado pelo software Epityper para a medição da metilação nos promotores A e C. Fez-se a avaliação da estabilidade do mRNA nas linhas celulares de cancro da mama MCF-7 e MDA-MB-231 tratadas com actinomicina D. Baixos níveis do transcrito ESR1_C associaram-se a uma melhor sobrevivência global (p = 0.017). Níveis elevados do transcrito ESR1_C associaram-se a uma resposta inferior ao tamoxifeno (HR = 2.48; CI 95% 1.24-4.99), um efeito mais pronunciado em doentes com tumores de fenótipo ERα/PgR duplamente positivo (HR = 3.41; CI 95% 1.45-8.04). A isoforma ESR1_C mostrou ter uma semi-vida prolongada, bem como uma estrutura secundária da região 5’UTR muito mais relaxada em comparação com a isoforma ESR1_A. A análise por Western-blot mostrou que ao nível da 21 proteína, a selectividade de promotores é indistinguivel. Não se detectou qualquer correlação entre os níveis das isoformas do gene ESR2 ou entre a metilação dos promotores do gene ESR2, e a detecção da proteína ERß. A metilação do promotor C do gene ESR1, e não do promotor A, foi responsável pela perda do receptor ERα. Estes resultados sugerem que os níveis do transcrito ESR1_C sejam usados como um novo potencial marcador para o prognóstico e predição de resposta ao tratamento com tamoxifeno em doentes com cancro da mama. Abstract: The decision of endocrine breast cancer treatment relies on ERα IHC-based assessment. However, ER positivity does not predict response in all cases in part due to IHC methodological limitations. We investigated whether ESR1 and ESR2 gene expression and respective promoter methylation may be related to non-favorable outcome of a proportion of tamoxifen treated patients as well as to ERα and ERß loss. Formalin-fixed paraffin-embedded breast cancer samples from 211 patients diagnosed between 1988 and 2004 were submitted to IHC-based ERα and ERß protein determination. ESR1 whole mRNA and promoter C specific transcript levels, as well as ESR2_ß1, ESR2_ß2/cx, and ESR2_ß5 transcripts were assessed by real-time PCR. ESR1 promoters A and C, and ESR2 promoters 0N and 0K were investigated by CpG methylation analysis using bisulfite-PCR for restriction analysis, or methylation specific PCR. Due to the promising results related to ESR1 promoter methylation, we have used a quantification method by matrixassisted laser desorption/ionization time-of-flight mass spectrometry (MALDITOF MS) together with Epityper software to measure methylation at promoters A and C. mRNA stability was assessed in actinomycin D treated MCF-7 and MDA-MB-231 cells. ERα protein was quantified using transiently transfected breast cancer cells. Low ESR1_C transcript levels were associated with better overall survival (p = 0.017). High levels of ESR1_C transcript were associated with non-favorable response in tamoxifen treated patients (HR = 2.48; CI 95% 1.24-4.99), an effect that was more pronounced in patients with ERα/PgR double-positive tumors (HR = 3.41; CI 95% 1.45-8.04). The ESR1_C isoform had a prolonged mRNA half-life and a more relaxed 5’UTR structure compared to ESR1_A isoform. Western-blot analysis showed that at protein level, the promoter selectivity is undistinguishable. There was no correlation between levels of ESR2 isoforms or ESR2 promoter methylation and ERß protein staining. ESR1 promoter C CpG methylation and not promoter A was responsible for ERα loss. We propose ESR1_C levels as a putative novel marker for breast cancer prognosis and prediction of tamoxifen response.

Time-domain optimization of amplifiers based on distributed genetic algorithms

Thesis presented in partial fulfillment of the requirements for the degree of Doctor of Philosophy in the subject of Electrical and Computer Engineering

Genetic characterization of Portuguese Fasciola hepatica isolates

Dissertation presented to obtain the Master Degree in Molecular, Genetics and Biomedicine

Analysis of PAX2 gene alterations in renal cell tumors

Dissertation for applying to a Master’s Degree in Molecular Genetics and Biomedicine submitted to the Sciences and Technology Faculty of New University of Lisbon

Functional genomic analysis of heat stress in Vitis vinifera

Grapevine (Vitis vinifera) is one of most agro-economically important fruit crops worldwide, with a special relevance in Portugal where over 300 varieties are used for wine production. Due to global warming, temperature stress is currently a serious issue affecting crop production especially in temperate climates. Mobile genetic elements such as retrotransposons have been shown to be involved in environmental stress induced genetic and epigenetic modifications. In this study, sequences related to Grapevine Retrotransposon 1 (Gret1) were utilized to determine heat induced genomic and transcriptomic modifications in Touriga Nacional, a traditional Portuguese grapevine variety. For this purpose, growing canes were treated to 42 oC for four hours and leaf genomic DNA and RNA was utilized for various techniques to observe possible genomic alterations and variation in transcription levels of coding and non-coding sequences between non-treated plants and treated plants immediately after heat stress (HS-0 h) or after a 24 hour recovery period (HS-24 h). Heat stress was found to induce a significant decrease in Gret1 related sequences in HS-24 h leaves, indicating an effect of heat stress on genomic structure. In order to identify putative heat induced DNA modifications, genome wide approaches such as Amplified Fragment Length Polymorphism were utilized. This resulted in the identification of a polymorphic DNA fragment in HS-0 h and HS-24 h leaves whose sequence mapped to a genomic region flanking a house keeping gene (NADH) that is represented in multiple copies in the Vitis vinifera genome. Heat stress was also found to affect the transcript levels of various non-coding and gene coding sequences. Accordingly, quantitative real time PCR results established that Gret1 related sequences are up regulated immediately after heat stress whereas the level of transcript of genes involved in identification and repair of double strand breaks are significantly down regulated in HS-0 h plants. Taken together, the results of this work demonstrated heat stress affects both genomic integrity and transcription levels.

Analysis of outcome and cancer stem cells status in patients with resectable pancreatic adenocarcinoma

RESUMO: Actualmente, a única possibilidade de cura para doentes com adenocarcinoma do pâncreas (PDAC) é a ressecção cirúrgica, no início deste estudo, perguntamo-nos se os predictores clínico-patológicos clássicos de prognostico poderiam ser validados em uma grande cohort de doentes com cancro do pâncreas ressecável e se outros predictores clínicos poderiam ter um papel na decisão de que doentes beneficiariam de ressecção cirúrgica. No capítulo 2, observamos que até 30% dos doentes morrem no primeiro ano após a ressecção cirúrgica, pelo que o nosso objectivo foi determinar factores pré-operatórios que se correlacionam com mortalidade precoce após ressecação cirúrgica com recurso a um instrumento estatisticamente validado, o Charlson-Age Comorbidity Index (CACI), determinamos que um CACI score superior a 4 foi preditivo de internamentos prolongados (p <0,001), complicações pós-operatórias (p = 0,042), e mortalidade em 1 ano pós- ressecção cirúrgica (p <0,001). Um CACI superior a 6 triplicou a mortalidade no primeiro ano pós-cirurgia e estes doentes têm menos de 50% de probabilidade de estarem vivos um ano após a cirurgia. No capítulo 3, o nosso objectivo foi identificar uma proteína de superfície que se correlacionasse estatisticamente com o prognostico de doentes com adenocarcinoma do pâncreas e permitisse a distinção de subgrupos de doentes de acordo com as suas diferenças moleculares, perguntamo-nos ainda se essa proteína poderia ser um marcador de células-estaminais. No nosso trabalho anterior observamos que as células tumorais na circulação sanguínea apresentavam genes com características bifenotípica epitelial e mesenquimal, enriquecimento para genes de células estaminais (ALDH1A1 / ALDH1A2 e KLF4), e uma super-expressão de genes da matriz extracelular (colagénios, SPARC, e DCN) normalmente identificados no estroma de PDAC. Após a avaliação dos tumores primários com RNA-ISH, muitos dos genes identificados, foram encontrados co-localizando em uma sub-população de células na região basal dos ductos pancreáticos malignos. Além disso, observamos que estas células expressam o marcador SV2A neuroendócrino, e o marcador de células estaminais ALDH1A1/2. Em comparação com tumores negativos para SV2, os doentes com tumores SV2 positivos apresentaram níveis mais baixos de CA 19-9 (69% vs. 52%, p = 0,012), tumores maiores (> 4 cm, 23% vs. 10%, p = 0,0430), menor invasão de gânglios linfáticos (69% vs. 86%, p = 0,005) e tumores mais diferenciados (69% vs. 57%, p = 0,047). A presença de SV2A foi associada com uma sobrevida livre de doença mais longa (HR: 0,49 p = 0,009) bem como melhor sobrevida global (HR: 0,54 p = 0,018). Em conjunto, esta informação aponta para dois subtipos diferentes de adenocarcinoma do pâncreas, e estes subtipos co-relacionam estatisticamente com o prognostico de doentes, sendo este subgrupo definido pela presença do clone celular SV2A / ALDH1A1/2 positivo com características neuroendócrinas. No Capítulo 4, a expressão de SV2A no cancro do pâncreas foi validado em linhas celulares primárias. Demonstramos a heterogeneidade do adenocarcinoma do pâncreas de acordo com características clonais neuroendócrinas. Ao comparar as linhas celulares expressando SV2 com linhas celulares negativas, verificamos que as linhas celulares SV2+ eram mais diferenciadas, diferindo de linhas celulares SV2 negativas no que respeita a mutação KRAS, proliferação e a resposta à quimioterapia. No capítulo 5, perguntamo-nos se o clone celular SV2 positivo poderia explicar a resistência a quimioterapia observada em doentes. Observamos um aumento absoluto de clones celulares expressando SV2A, em múltiplas linhas de evidência - doentes, linhas de células primárias e xenotransplantes. Embora, tenhamos sido capazes de demonstrar que o adenocarcinoma do pâncreas é uma doença heterogénea, consideramos que a caracterização genética destes clones celulares expressando SV2A é de elevada importância. Pretendemos colmatar esta limitação com as seguintes estratégias: Após o tratamento com quimioterapia neoadjuvante na nossa coorte, realizamos microdissecação a laser das amostras primarias em parafina, de forma a analisar mutações genéticas observadas no adenocarcinoma pancreático; em segundo lugar, pretendemos determinar consequências de knockdown da expressão de SV2A em nossas linhas celulares seguindo-se o tratamento com gemicitabina para determinação do papel funcional de SV2A; finalmente, uma vez que os nossos esforços anteriores com um promotor - repórter e SmartFlare ™ falharam, o próximo passo será realizar RNA-ISH PrimeFlow™ seguido de FACS e RNA-seq para caracterização deste clone celular. Em conjunto, conseguimos provar com várias linhas de evidência, que o adenocarcinoma pancreático é uma doença heterogénea, definido por um clone de células que expressam SV2A, com características neuroendócrinas. A presença deste clone no tecido de doentes correlaciona-se estatisticamente com o prognostico da doença, incluindo sobrevida livre de doença e sobrevida global. Juntamente com padrões de proliferação e co-expressão de ALDH1A1/2, este clone parece apresentar um comportamento de células estaminais e está associado a resistência a quimioterapia, uma vez que a sua expressão aumenta após agressão química, quer em doentes, quer em linhas de células primárias.----------------------------- ABSTRACT: Currently, the only chance of cure for patients with pancreatic adenocarcinoma is surgical resection, at the beginning of my thesis studies, we asked if the classical clinicopathologic predictors of outcome could be validated in a large cohort of patients with early stage pancreatic cancer and if other clinical predictors could have a role on deciding which patients would benefit from surgery. In chapter 2, we found that up to 30% of patients die within the first year after curative intent surgery for pancreatic adenocarcinoma. We aimed at determining pre-operative factors that would correlate with early mortality following resection for pancreatic cancer using a statistically validated tool, the Charlson-Age Comorbidity Index (CACI). We found that a CACI score greater than 4 was predictive of increased length of stay (p<0.001), post-operative complications (p=0.042), and mortality within 1-year of pancreatic resection (p<0.001). A CACI score of 6 or greater increased 3-fold the odds of death within the first year. Patients with a high CACI score have less than 50% likelihood of being alive 1 year after surgery. In chapter 3 we aimed at identifying a surface protein that correlates with patient’s outcome and distinguishes sub-groups of patients according to their molecular differences and if this protein could be a cancer stem cell marker. The most abundant class of circulating tumor cells identified in our previous work was found to have biphenotypic features of epithelial to mesenchymal transition, enrichment for stem-cell associated genes (ALDH1A1/ALDH1A2 and KLF4), and an overexpression of extracellular matrix genes (Collagens, SPARC, and DCN) normally found in the stromal microenvironment of PDAC primary tumors. Upon evaluation of matched primary tumors with RNA-ISH, many of the genes identified were found to co-localize in a sub-population of cells at the basal region of malignant pancreatic ducts. In addition, these cells expressed the neuroendocrine marker SV2A, and the stem cell marker ALDH1A1/2. Compared to SV2 negative tumors, patients with SV2 positive tumors were more likely to present with lower CA 19-9 (69% vs. 52%, p = 0.012), bigger tumors (size > 4 cm, 23% vs. 10%, p= 0.0430), less nodal involvement (69% vs. 86%, p = 0.005) and lower histologic grade (69% vs. 57%, p = 0.047). The presence of SV2A expressing cells was associated with an improved disease free survival (HR: 0.49 p=0.009) and overall survival (HR: 0.54 p=0.018) and correlated linearly with ALDH1A2. Together, this information points to two different sub-types of pancreatic adenocarcinoma, and these sub-types correlated with patients’ outcome and were defined by the presence of a SV2A/ ALDH1A1/2 expressing clone with neuroendocrine features. In Chapter 4, SV2A expression in cancer was validated in primary cell lines. We were able to demonstrate pancreatic adenocarcinoma heterogeneity according to neuroendocrine clonal features. When comparing SV2 expressing cell lines with SV2 negative cell lines, we found that SV2+ cell lines were more differentiated and differ from SV2 negative cell lines regarding KRAS mutation, proliferation and response to chemotherapy. In Chapter 5 we aimed at determining if this SV2 positive clone could explain chemoresistance observed in patients. We found an absolute increase in SV2A expressing cells, with multiple lines of evidence, in patients, primary cell lines and xenografts. Although, we have been able to show evidence that pancreatic adenocarcinoma is a heterogeneous disease, our findings warrant further investigation. To further characterize SV2A expressing clones after treatment with neoadjuvant chemotherapy in our cohort, we have performed laser capture microdissection of the paraffin embedded tissue in this study and will analyze the tissue for known genetic mutations in pancreatic adenocarcinoma; secondly, we want to know what will happen after knocking down SV2A expression in our cell lines followed by treatment with gemcitabine to determine if SV2A is functionally important; finally, since our previous efforts with a promoter – reporter and SmartFlare™ have failed, we will utilize a novel PrimeFlow™ RNA-ISH assay followed by FACS and RNA sequencing to further characterize this cellular clone. Overall our data proves, with multiple lines of evidence, that pancreatic adenocarcinoma is a heterogeneous disease, defined by a clone of SV2A expressing cells, with neuroendocrine features. The presence of this clone in patients’ tissue correlates with patient’s disease free survival and overall survival. Together with patterns of proliferation and ALDH1A1/2 co-expression, this clone seems to present a stem-cell-like behavior and is associated with chemoresistance, since it increases after chemotherapy, both in patients and primary cell lines.

Genetic algorithm for shipping route estimation with long-range tracking data : automatic reconstruction of shipping routes based on the historical ship positions for maritime safety applications

Ship tracking systems allow Maritime Organizations that are concerned with the Safety at Sea to obtain information on the current location and route of merchant vessels. Thanks to Space technology in recent years the geographical coverage of the ship tracking platforms has increased significantly, from radar based near-shore traffic monitoring towards a worldwide picture of the maritime traffic situation. The long-range tracking systems currently in operations allow the storage of ship position data over many years: a valuable source of knowledge about the shipping routes between different ocean regions. The outcome of this Master project is a software prototype for the estimation of the most operated shipping route between any two geographical locations. The analysis is based on the historical ship positions acquired with long-range tracking systems. The proposed approach makes use of a Genetic Algorithm applied on a training set of relevant ship positions extracted from the long-term storage tracking database of the European Maritime Safety Agency (EMSA). The analysis of some representative shipping routes is presented and the quality of the results and their operational applications are assessed by a Maritime Safety expert.