Uncovering functionally relevant signaling pathways using microarray-based expression profiling

The introduction of microarray technology to the scientific and medical communities has fundamentally altered the way in which we now address basic biomedical questions. Microarrays technology facilitates a more complete and inclusive experimental approach where alterations in the transcript level of entire genomes can be simultaneously assayed in response to a variety of stimuli. Conceptually different approaches to the development of microarray technology have resulted in the generation of two different array formats: oligonucleotide arrays and cDNA arrays. The application of microarray and related technologies to identify specific targets of defined genes that have clearly been implicated in cancer progression requires a specific experimental approach. The objective of tiffs approach is to define changes in transcriptional profile that occur in response to modulating the expression level of the gene to be studied. The resulting altered expression profile can then be viewed as a blueprint by which that gene effects its cellular function. We have used oligonucleotide array-based expression profiling in collaboration with Affymetrix to identify downstream transcriptional targets of the BRCA1 tumor-suppressor gene as a means of defining its function. BRCA1 has been implicated in at least three functional pathways, namely, mediating the cellular response to DNA damage, as a cell cycle checkpoint protein and in the regulation of transcription. The physiological significance of these properties and their implications for the function of BRCA1 as a tumor-suppressor gene remain to be defined.

The androgen receptor controls expression of the cancer-associated sTn antigen and cell adhesion through induction of ST6GalNAc1 in prostate cancer

Patterns of glycosylation are important in cancer, but the molecular mechanisms that drive changes are often poorly understood. The androgen receptor drives prostate cancer (PCa) development and progression to lethal metastatic castration-resistant disease. Here we used RNA-Seq coupled with bioinformatic analyses of androgen-receptor (AR) binding sites and clinical PCa expression array data to identify ST6GalNAc1 as a direct and rapidly activated target gene of the AR in PCa cells. ST6GalNAc1 encodes a sialytransferase that catalyses formation of the cancer-associated sialyl-Tn antigen (sTn), which we find is also induced by androgen exposure. Androgens induce expression of a novel splice variant of the ST6GalNAc1 protein in PCa cells. This splice variant encodes a shorter protein isoform that is still fully functional as a sialyltransferase and able to induce expression of the sTn-antigen. Surprisingly, given its high expression in tumours, stable expression of ST6GalNAc1 in PCa cells reduced formation of stable tumours in mice, reduced cell adhesion and induced a switch towards a more mesenchymal-like cell phenotype in vitro. ST6GalNAc1 has a dynamic expression pattern in clinical datasets, beingsignificantly up-regulated in primary prostate carcinoma but relatively down-regulated in established metastatic tissue. ST6GalNAc1 is frequently upregulated concurrently with another important glycosylation enzyme GCNT1 previously associated with prostate cancer progression and implicated in Sialyl Lewis X antigen synthesis. Together our data establishes an androgen-dependent mechanism for sTn antigen expression in PCa, and are consistent with a general role for the androgen receptor in driving important coordinate changes to the glycoproteome during PCa progression.

Regulation of human dUTPase gene expression and p53-mediated transcriptional repression in response to oxaliplatin-induced DNA damage

Deoxyuridine triphosphate nucleotidohydrolase (dUTPase) catalyzes the hydrolysis of dUTP to dUMP and PPi. Although dUTP is a normal intermediate in DNA synthesis, its accumulation and misincorporation into DNA is lethal. Importantly, uracil misincorporation is a mechanism of cytotoxicity induced by fluoropyrimidine chemotherapeutic agents including 5-fluorouracil (5-FU) and elevated expression of dUTPase is negatively correlated with clinical response to 5-FU-therapy. In this study we performed the first functional characterization of the dUTPase promoter and demonstrate a role for E2F-1 and Sp1 in driving dUTPase expression. We establish a direct role for both mutant and wild-type forms of p53 in modulating dUTPase promoter activity. Treatment of HCT116 p53(+/+) cells with the DNA-damaging agent oxaliplatin induced a p53-dependent transcriptional downregulation of dUTPase not observed in the isogenic null cell line. Oxaliplatin treatment induced enrichment of p53 at the dUTPase promoter with a concomitant reduction in Sp1. The suppression of dUTPase by oxaliplatin promoted increased levels of dUTP that was enhanced by subsequent addition of fluoropyrimidines. The novel observation that oxaliplatin downregulates dUTPase expression may provide a mechanistic basis contributing to the synergy observed between 5-FU and oxaliplatin in the clinic. Furthermore, these studies provide the first evidence of a direct transcriptional link between the essential enzyme dUTPase and the tumor suppressor p53.

Protease activated receptor 2 expression and functionality in caries

Introduction: Protease activated receptors (PARs) are G-protein-coupled transmembrane receptors that are expressed on many cell types and implicated in various inflammatory processes in vivo. The induction of PAR2 as a result of the inflammatory response associated with dental caries remains to be determined. Objectives: The aim was to localise the expression of PAR2 in human dental pulp from carious teeth and to confirm receptor functionality using an in vitro assay. Methods: Dental pulp sections from decalcified carious teeth were examined by immunocytochemsitry. Membrane preparations from cultured pulp fibroblasts were subject to SDS-PAGE and immunoblotting to confirm fibroblast-associated immunoreactivity. The functionality of PAR2 on dental pulp fibroblasts was studied using calcium imaging in the presence of several potential activators including a PAR2 agonist (PAR2-AP), trypsin and pulpal enzymes from a carious tooth. Results: Immunocytochemistry revealed intense PAR2 immunoreactivity on pulpal fibroblasts subjacent to carious lesions but not in surrounding regions of the dental pulp. Pulp specimens from a dental injury model showed no expression of PAR2, suggesting its expression was related to cellular changes associated with ongoing caries. The localisation of PAR2 staining to pulpal fibroblasts in carious teeth was confirmed by Western blotting which revealed PAR2 immunoreactive bands in membrane fractions prepared from pulp fibroblasts. In functional studies, challenge of cultured pupal fibroblasts with PAR2-AP, trypsin and an extract of proteolytic enzymes from a carious dental pulp, showed specific activation of PAR2. Conclusions: This work demonstrates that PAR2 is functional and inducible in human dental pulp fibroblasts in response to caries and that endogenous pulpal enzymes can activate PAR2.

Sequence and expression of complement factor H gene cluster variants and their roles in age-related macular degeneration risk

Purpose: To investigate how potentially functional genetic variants are coinherited on each of four common complement factor H (CFH) and CFH-related gene haplotypes and to measure expression of these genes in eye and liver tissues.

Methods: We sequenced the CFH region in four individuals (one homozygote for each of four common CFH region haplotypes) to identify all genetic variants. We studied associations between the haplotypes and AMD phenotypes in 2157 cases and 1150 controls. We examined RNA-seq profiles in macular and peripheral retina and retinal pigment epithelium/choroid/sclera (RCS) from eight eye donors and three liver samples.

Results: The haplotypic coinheritance of potentially functional variants (including missense variants, novel splice sites, and the CFHR3–CFHR1 deletion) was described for the four common haplotypes. Expression of the short and long CFH transcripts differed markedly between the retina and liver. We found no expression of any of the five CFH-related genes in the retina or RCS, in contrast to the liver, which is the main source of the circulating proteins.

Conclusions: We identified all genetic variants on common CFH region haplotypes and described their coinheritance. Understanding their functional effects will be key to developing and stratifying AMD therapies. The small scale of our expression study prevented us from investigating the relationships between CFH region haplotypes and their expression, and it will take time and collaboration to develop epidemiologic-scale studies. However, the striking difference between systemic and ocular expression of complement regulators shown in this study suggests important implications for the development of intraocular and systemic treatments.

Network signatures based on gene pair expression ratios improve classification and the analysis of muscle-invasive urothelial cancer

Urothelial cancer (UC) is highly recurrent and can progress from non-invasive (NMIUC) to a more aggressive muscle-invasive (MIUC) subtype that invades the muscle tissue layer of the bladder. We present a proof of principle study that network-based features of gene pairs can be used to improve classifier performance and the functional analysis of urothelial cancer gene expression data. In the first step of our procedure each individual sample of a UC gene expression dataset is inflated by gene pair expression ratios that are defined based on a given network structure. In the second step an elastic net feature selection procedure for network-based signatures is applied to discriminate between NMIUC and MIUC samples. We performed a repeated random subsampling cross validation in three independent datasets. The network signatures were characterized by a functional enrichment analysis and studied for the enrichment of known cancer genes. We observed that the network-based gene signatures from meta collections of proteinprotein interaction (PPI) databases such as CPDB and the PPI databases HPRD and BioGrid improved the classification performance compared to single gene based signatures. The network based signatures that were derived from PPI databases showed a prominent enrichment of cancer genes (e.g., TP53, TRIM27 and HNRNPA2Bl). We provide a novel integrative approach for large-scale gene expression analysis for the identification and development of novel diagnostical targets in bladder cancer. Further, our method allowed to link cancer gene associations to network-based expression signatures that are not observed in gene-based expression signatures.

Evolution and functional relevance of Ataxin-3 paralogues

O gene ataxin-3 (ATXN3; 14q32.1) codifica uma proteína expressa ubiquamente, envolvida na via ubiquitina-proteassoma e na repressão da transcrição. Grande relevância tem sido dada ao gene ATXN3 após a identificação de uma expansão (CAG)n na sua região codificante, responsável pela ataxia mais comum em todo o mundo, SCA3 ou doença de Machado-Joseph (DMJ). A DMJ é uma doença neurodegenerativa, autossómica dominante, de início tardio. O tamanho do alelo expandido explica apenas uma parte do pleomorfismo da doença, evidenciando a importância do estudo de outros modificadores. Em doenças de poliglutaminas (poliQ), a toxicidade é causada por um ganho de função da proteína expandida; no entanto, a proteína normal parece ser, também, um dos agentes modificadores da patogénese. O gene ATXN3 possui dois parálogos humanos gerados por retrotransposição: ataxin-3 like (ATXN3L) no cromossoma X, e LOC100132280, ainda não caracterizado, no cromossoma 8. Estudos in vitro evidenciaram a capacidade da ATXN3L para clivar cadeias de ubiquitina, sendo o seu domínio proteolítico mais eficiente do que o domínio da ATXN3 parental. O objetivo deste estudo foi explorar a origem e a evolução das retrocópias ATXN3L e LOC100132280 (aqui denominadas ATXN3L1 e ATXN3L2), assim como testar a relevância funcional de ambas através de abordagens evolutivas e funcionais. Deste modo, para estudar a divergência evolutiva dos páralogos do gene ATXN3: 1) analisaram-se as suas filogenias e estimou-se a data de origem dos eventos de retrotransposição; 2) avaliaram-se as pressões seletivas a que têm sido sujeitos os três parálogos, ao longo da evolução dos primatas; e 3) explorou-se a evolução das repetições CAG, localizadas em três contextos genómicos diferentes, provavelmente sujeitos a diferentes pressões seletivas. Finalmente, para o retrogene que conserva uma open reading frame (ORF) intacta, ATXN3L1, analisou-se, in silico, a conservação dos locais e domínios proteicos da putativa proteína. Ademais, para este retrogene, foi estudado o padrão de expressão de mRNA, através da realização de PCR de Transcriptase Reversa, em 16 tecidos humanos. Os resultados obtidos sugerem que dois eventos independentes de retrotransposição estiveram na origem dos retrogenes ATXN3L1 e ATXN3L2, tendo o primeiro ocorrido há cerca de 63 milhões de anos (Ma) e o segundo após a divisão Platirrínios-Catarrínios, há cerca de 35 Ma. Adicionalmente, outras retrocópias foram encontradas em primatas e outros mamíferos, correspondendo, no entanto, a eventos mais recentes e independentes de retrotransposição. A abordagem evolutiva mostrou a existência de algumas constrições selectivas associadas à evolução do gene ATXN3L1, à semelhança do que acontece com ATXN3. Por outro lado, ATXN3L2 adquiriu codões stop prematuros que, muito provavelmente, o tornaram num pseudogene processado. Os resultados da análise de expressão mostraram que o gene ATXN3L1 é transcrito, pelo menos, em testículo humano; no entanto, a optimização final da amplificação específica dos transcriptos ATXN3L1 permitirá confirmar se a expressão se estende a outros tecidos. Relativamente ao mecanismo de mutação inerente à repetição CAG, os dois parálogos mostraram diferentes padrões de evolução: a retrocópia ATXN3L1 é altamente interrompida e pouco polimórfica, enquanto a ATXN3L2 apresenta tratos puros de (CAG)n em algumas espécies e tratos hexanucleotídicos de CGGCAG no homem e no chimpanzé. A recente aquisição da repetição CGGCAG pode ter resultado de uma mutação inicial de CAG para CGG, seguida de instabilidade que proporcionou a expansão dos hexanucleótidos.Estudos futuros poderão ser realizados no sentido de confirmar o padrão de expressão do gene ATXN3L1 e de detetar proteína endógena in vivo. Adicionalmente, a caracterização da proteina ataxina-3 like 1 e dos seus interatores moleculares poderá povidenciar informação acerca da sua relevância no estado normal e patológico.

Functional analysis of String/CDC25 phosphatase in post-mitotic neurons

Cell cycle and differentiation are two highly coordinated processes during organ development. Recent studies have demonstrated that core cell cycle regulators also play cell cycle-independent functions in post-mitotic neurons, and are essential for the maintenance of neuronal homeostasis. CDC25 phosphatases are well-established CDK activators and their activity is mainly associated to proliferating tissues. The expression and activity of mammalian CDC25s has been reported in adult brains. However, their physiological relevance and the potential substrates in a non-proliferative context have never been addressed. string (stg) encodes the Drosophila CDC25 homolog. Previous studies from our group showed that stg is expressed in photoreceptors (PRs) and in lamina neurons, which are two differentiated cell types that compose the fly visual system. The aims of this work are to uncover the function of stg and to identify its potential neuronal substrates, using the Drosophila visual system as a model. To gain insight into the function of stg in a non-dividing context we used the GAL4/UAS system to promote downregulation of stg in PR-neurons, through the use of an RNAi transgene. The defects caused by stg loss-of-function were evaluated in the developing eye imaginal disc by immunofluorescence, and during adult stages by scanning electron microscopy. This genetic approach was combined with a specific proteomic method, two-dimensional difference gel electrophoresis (2D-DIGE), to identify the potential substrates in PR-cells. Our results showed that stg downregulation in PRs affects the well-patterned retina organization, inducing the loss of apical maintenance of PR-nuclei on the eye disc, and ommatidia disorganization. We also detected an abnormal accumulation of cytoskeletal proteins and a disruption of the axon structure. As a consequence, the projection of PR-axons into the lamina and medulla neuropils of the optic lobe was impaired. Upon stg downregulation, we also detected that PR-cells accumulate Cyclin B. Although the rough eye phenotype observed upon stg downregulation suggests neurodegeneration, we did not detect neuronal death during larval stages, suggesting that it likely occurs during pupal stages or during adulthood. By 2D-DIGE, we identified seven proteins which were differentially expressed upon stg downregulation, and are potential neuronal substrates of Stg. Altogether, our observations suggest that Stg phosphatase plays an essential role in the Drosophila visual system neurons, regulating several cell components and processes in order to ensure their homeostasis.

Evolution, gene regulation and functional analysis of BMP2 in fish

Bone morphogenetic proteins (BMPs) are multifunctional growth factors belonging to the transforming growth factor β (TGFβ) superfamily with a central role in bone formation and mineralization. BMP2, a founding member of this family, has demonstrated remarkable osteogenic properties and is clinically used to promote bone repair and fracture healing. Lack of basic data on factors regulating BMP2 expression and activity have hampered a better understanding of its role in bone formation and bone-related diseases. The objective of this work was to collect new functional data and determine spatiotemporal expression patterns in a fish system aiming towards a better understanding of BMP2 function and regulation. Transcriptional and post-transcriptional regulation of gilthead seabream BMP2 gene was inferred from luciferase reporter systems. Several bone- and cartilage-related transcription factors (e.g. RUNX3, MEF2c, SOX9 and ETS1) were found to regulate BMP2 transcription, while microRNA 20a was shown to affect stability of the BMP2 transcript and thus the mineralogenic capacity of fish bone-derived host cells. The regulation of BMP2 activity through an interaction with the matrix Gla protein (MGP) was investigated in vitro using BMP responsive elements (BRE) coupled to luciferase reporter gene. Although we demonstrated the functionality of the experimental system in a fish cell line and the activation of BMP signaling pathway by seabream BMP2, no conclusive evidence could be collected on a possible interaction beween MGP and BMP2. The evolutionary relationship among the members of BMP2/4/16 subfamily was inferred from taxonomic and phylogenetic analyses. BMP16 diverged prior to BMP2 and BMP4 and should be the result of an ancient genome duplication that occurred early in vertebrate evolution. Structural and functional data suggested that all three proteins are effectors of the BMP signaling pathway, but expression data revealed different spatiotemporal patterns in teleost fish suggesting distinct mechanisms of regulation. In this work, through the collection of novel data, we provide additional insight into the regulation, the structure and the phylogenetic relationship of BMP2 and its closely related family members.

Molecular characterization of Gla-rich protein (GRP) gene expression and function

Gla-rich protein (GRP) is a vitamin K-dependent protein related to bone and cartilage recently described. This protein is characterized by a large number of Gla (γ-carboxyglutamic acid) residues being the protein with the highest Gla content of any known protein. It was found in a widely variety of tissues but highest levels was found in skeletal and cartilaginous tissues. This small secreted protein was also expressed and accumulated in soft tissues and it was clearly associated with calcification pathologies in the same tissues. Although the biological importance of GRP remains to be elucidated, it was suggested a physiological role in cartilage development and calcification process during vertebrate skeleton formation. Using zebrafish, an accepted model to study skeletal development, we have described two grp paralog genes, grp1 and grp2, which exhibited distinct patterns of expression, suggesting different regulatory pathways for each gene. Gene synteny analysis showed that grp2 gene is more closely related to tetrapod grp, although grp1 gene was proposed to be the vertebrate ortholog by sequence comparison. In addition, we identified a functional promoter of grp2 gene and using a functional approach we confirmed the involvement of transcription factors from Sox family (Sox9b and Sox10) in the regulation of grp2 expression. In an effort to provide more information about the function of grp isoforms, we generated two zebrafish transgenic lines capable to overexpress conditionally grp genes and possible roles in the skeleton development were studied. To better understand GRP function a mammalian system was used and the analysis of knockout mice showed that GRP is involved in chondrocyte maturation and the absence of GRP is associated to proteoglycans loss in calcified articular cartilage. In addition, we detected differences in chondrogenesis markers in articular chondrocyte primary culture. Overall, our data suggest a main role for GRP on chondrocyte differentiation.

Molecular structure and functional analysis of runx3 in zebrafish

Tese de doutoramento, Ciências Biomédicas, Universidade do Algarve, Departamento de Ciências Biomédicas e Medicina, 2014

Identification and molecular characterization of bone-related micrornas: functional implications

Tese de doutoramento, Ciências Biomédicas, Departamento de Ciências Biomédicas e Medicina, Universidade do Algarve, 2014

The functional link between HFE and hepcidin- its contribution to iron metabolism regulation

Tese de doutoramento, Biologia (Biologia Molecular), Universidade de Lisboa, Faculdade de Ciências, 2013

The role of reciprocal regulation between TAL1 and miRNA expression in T-cell leukemia progression

Tese de doutoramento, Ciências Biomédicas (Biologia Celular e Molecular), Universidade de Lisboa, Faculdade de Medicina, 2015