Gene Expression : From Microarrays to Functional Genomics

The time of the large sequencing projects has enabled unprecedented possibilities of investigating more complex aspects of living organisms. Among the high-throughput technologies based on the genomic sequences, the DNA microarrays are widely used for many purposes, including the measurement of the relative quantity of the messenger RNAs. However, the reliability of microarrays has been strongly doubted as robust analysis of the complex microarray output data has been developed only after the technology had already been spread in the community. An objective of this study consisted of increasing the performance of microarrays, and was measured by the successful validation of the results by independent techniques. To this end, emphasis has been given to the possibility of selecting candidate genes with remarkable biological significance within specific experimental design. Along with literature evidence, the re-annotation of the probes and model-based normalization algorithms were found to be beneficial when analyzing Affymetrix GeneChip data. Typically, the analysis of microarrays aims at selecting genes whose expression is significantly different in different conditions followed by grouping them in functional categories, enabling a biological interpretation of the results. Another approach investigates the global differences in the expression of functionally related groups of genes. Here, this technique has been effective in discovering patterns related to temporal changes during infection of human cells. Another aspect explored in this thesis is related to the possibility of combining independent gene expression data for creating a catalog of genes that are selectively expressed in healthy human tissues. Not all the genes present in human cells are active; some involved in basic activities (named housekeeping genes) are expressed ubiquitously. Other genes (named tissue-selective genes) provide more specific functions and they are expressed preferably in certain cell types or tissues. Defining the tissue-selective genes is also important as these genes can cause disease with phenotype in the tissues where they are expressed. The hypothesis that gene expression could be used as a measure of the relatedness of the tissues has been also proved. Microarray experiments provide long lists of candidate genes that are often difficult to interpret and prioritize. Extending the power of microarray results is possible by inferring the relationships of genes under certain conditions. Gene transcription is constantly regulated by the coordinated binding of proteins, named transcription factors, to specific portions of the its promoter sequence. In this study, the analysis of promoters from groups of candidate genes has been utilized for predicting gene networks and highlighting modules of transcription factors playing a central role in the regulation of their transcription. Specific modules have been found regulating the expression of genes selectively expressed in the hippocampus, an area of the brain having a central role in the Major Depression Disorder. Similarly, gene networks derived from microarray results have elucidated aspects of the development of the mesencephalon, another region of the brain involved in Parkinson Disease.

A promoter-level mammalian expression atlas

Regulated transcription controls the diversity, developmental pathways and spatial organization of the hundreds of cell types that make up a mammal. Using single-molecule cDNA sequencing, we mapped transcription start sites (TSSs) and their usage in human and mouse primary cells, cell lines and tissues to produce a comprehensive overview of mammalian gene expression across the human body. We find that few genes are truly 'housekeeping', whereas many mammalian promoters are composite entities composed of several closely separated TSSs, with independent cell-type-specific expression profiles. TSSs specific to different cell types evolve at different rates, whereas promoters of broadly expressed genes are the most conserved. Promoter-based expression analysis reveals key transcription factors defining cell states and links them to binding-site motifs. The functions of identified novel transcripts can be predicted by coexpression and sample ontology enrichment analyses. The functional annotation of the mammalian genome 5 (FANTOM5) project provides comprehensive expression profiles and functional annotation of mammalian cell-type-specific transcriptomes with wide applications in biomedical research.

FGFR1 regulated gene-expression, cell proliferation and differentiation in the developing midbrain and hindbrain

The neuroectodermal tissue close to the midbrain hindbrain boundary (MHB) is an important secondary organizer in the developing neural tube. This so-called isthmic organizer (IsO) regulates cellular survival, patterning and proliferation in the midbrain (Mb) and rhombomere 1 (R1) of the hindbrain. Signaling molecules of the IsO, such as fibroblast growth factor 8 (FGF8) and WNT1 are expressed in distinct bands of cells around the MHB. It has been previously shown that FGF-receptor 1 (FGFR1) is required for the normal development of this brain region in the mouse embryo. In the present study, we have compared the gene expression profiles of wild-type and Fgfr1 mutant embryos. We show that the loss of Fgfr1 results in the downregulation of several genes expressed close to the MHB and in the disappearance of gene expression gradients in the midbrain and R1. Our microarray screen identified several previously uncharacterized genes which may participate in the development of midbrain R1 region. Our results also show altered neurogenesis in the midbrain and R1 of the Fgfr1 mutants. Interestingly, the neuronal progenitors in midbrain and R1 show different responses to the loss of signaling through FGFR1. As Wnt1 expression at the MHB region requires the FGF signaling pathway, WNT target genes, including Drapc1, were also identified in our screen. The microarray data analysis also suggested that the cells next to the midbrain hindbrain boundary express distinct cell cycle regulators. We showed that the cells close to the border appeared to have unique features. These cells proliferate less rapidly than the surrounding cells. Unlike the cells further away from the boundary, these cells express Fgfr1 but not the other FGF receptors. The slowly proliferating boundary cells are necessary for development of the characteristic isthmic constriction. They may also contribute to compartmentalization of this brain region.

CD1 expression and CD1-restricted T cell activity in normal and tumour-bearing human liver

CD1d-restricted natural killer T (NKT) cells expressing invariant Valpha14Jalpha18 T cell receptor alpha-chains are abundant in murine liver and are implicated in the control of malignancy, infection and autoimmunity. Invariant NKT cells have potent anti-metastatic effects in mice and phase I clinical trials involving their homologues in humans are ongoing. However, invariant NKT cells are less abundant in human liver ( approximately 0.5% of hepatic T cells) than in murine liver (up to 50%) and it is not known if other hepatic T cells are CD1-restricted. We have examined expression of CD1a, CD1b, CD1c and CD1d mRNA and protein in human liver and evaluated the reactivity of mononuclear cells (MNC) from histologically normal and tumour-bearing human liver specimens against these CD1 isoforms. Messenger RNA for all CD1 isotypes was detectable in all liver samples. CD1c and CD1d were expressed at the protein level by hepatic MNC. CD1d, only, was detectable at the cell surface, but CD1c and CD1d were found at an intracellular location in significant numbers of liver MNC. CD1b was not expressed by MNC from healthy livers but was detectable within MNC in all tumour samples tested. Hepatic T cells exhibited reactivity against C1R cells expressing transfected CD1c and CD1d, but neither CD1a nor CD1b. These cells secreted interferon-gamma (IFN-gamma) but not interleukin-4 (IL-4) upon stimulation. In contrast, similar numbers of peripheral T cells released 13- and 16-fold less IFN-gamma in response to CD1c and CD1d, respectively. CD1c and CD1d expression and T cell reactivity were not altered in tumour-bearing liver specimens compared to histologically normal livers. These data suggest that, in addition to invariant CD1d-restricted NKT cells, autoreactive T cells that recognise CD1c and CD1d and release inflammatory cytokines are abundant in human liver.

Genetic aspects of outcome in kidney transplantation: cytokine and thrombosis associated candidate genes and gene expression biomarkers

Kidney transplantation (Tx) is the treatment of choice for end stage renal disease. Immunosuppressive medications are given to prevent an immunological rejection of the transplant. However, immunosuppressive drugs increase e.g. the risk of infection, cancer or nephrotoxicity. A major genetic contributors to immunological acceptance of the graft are human leukocyte antigen (HLA) genes. Also other non-HLA gene polymorphisms may predict the future risk of complications before Tx, possibly enabling individualised immunotherapy. Graft function after Tx is monitored using non-specific clinical symptoms and laboratory markers. The definitive diagnosis of graft rejection however relies on a biopsy of the graft. In the acute rejection (AR) diagnostics there is a need for an alternative to biopsy that would be an easily repeatable and simple method for regular use. Frequent surveillance of acute or subclinical rejection (SCR) may improve long-term function. In this thesis, associations between cytokine and thrombosis associated candidate genes and the outcome of kidney Tx were studied. Cytotoxic and co-stimulatory T lymphocyte molecule gene expression biomarkers for the diagnosis of the AR and the SCR were also investigated. We found that polymorphisms in the cytokine genes tumor necrosis factor and interleukin 10 (IL10) of the recipients were associated with AR. In addition, certain IL10 gene polymorphisms of the donors were associated with the incidence of cytomegalovirus infection and occurrence of later infection in a subpopulation of recipients. Further, polymorphisms in genes related to the risk of thrombosis and those of certain cytokines were not associated with the occurrence of thrombosis, infarction, AR or graft survival. In the study of biomarkers for AR, whole blood samples were prospectively collected from adult kidney Tx patients. With real-time quantitative PCR (RT-QPCR) gene expression quantities of CD154 and ICOS differentiated the patients with AR from those without, but not from the patients with other causes of graft dysfunction. Biomarkers for SCR were studied in paediatric kidney Tx patients. We used RT-QPCR to quantify the gene expression of immunological candidate genes in a low-density array format. In addition, we used RT-QPCR to validate the results of the microarray analysis. No gene marker differentiated patients with SCR from those without SCR. This research demonstrates the lack of robust markers among polymorphisms or biomarkers in investigated genes that could be included in routine analysis in a clinical laboratory. In genetic studies, kidney Tx can be regarded as a complex trait, i.e. several environmental and genetic factors may determine its outcome. A number of currently unknown genetic factors probably influence the results of Tx.

Expression of glucose transporter isoforms and the insulin receptor during hamster preimplantation embryo development

During preimplantation development, embryos of many species are known to express up to five isoforms of the facilitative glucose transporter proteins (GLUT). Development of hamster blastocysts is inhibited by glucose. We therefore investigated GLUT isoform and insulin receptor (IR) expression in hamster preimplantation embryos cultured in glucose-free medium from the 8-cell stage onwards. We show that GLUT1, 3 and 8 mRNA are constitutively expressed from the 8-cell to the blastocyst stage. The IR is expressed from the morula stage onwards. Messenger RNA of the insulin-responsive GLUT4 was not detected at any stage. GLUT1 and 3 were localised by immunocytochemistry. GLUT1 was expressed in both embryoblast and trophoblast, in the latter, mainly in basal and lateral membranes directed towards the blastocoel. and embryoblast. GLUT3 was exclusively localised in the apical. membrane of trophoblast cells. We show that hamster preimplantation embryos express several GLUT isoforms thus closely resembling embryos of other mammalian species. Despite endogenous IR expression, the insulin-sensitive isoform GLUT4 was not expressed, indicating that the insulin-mediated glucose uptake known from classical insulin target cells may not be relevant for hamster blastocysts.

Expression, Enzymatic Activities and Subcellular Localization of Hepatitis E Virus and Semliki Forest Virus Replicase Proteins

Viruses are submicroscopic, infectious agents that are obligate intracellular parasites. They adopt various types of strategies for their parasitic replication and proliferation in infected cells. The nucleic acid genome of a virus contains information that redirects molecular machinery of the cell to the replication and production of new virions. Viruses that replicate in the cytoplasm and are unable to use the nuclear transcription machinery of the host cell have developed their own transcription and capping systems. This thesis describes replication strategies of two distantly related viruses, hepatitis E virus (HEV) and Semliki Forest virus (SFV), which belong to the alphavirus-like superfamily of positive-strand RNA viruses. We have demonstrated that HEV and SFV share a unique cap formation pathway specific for alphavirus-like superfamily. The capping enzyme first acts as a methyltransferase, catalyzing the transfer of a methyl group from S-adenosylmethionine to GTP to yield m7GTP. It then transfers the methylated guanosine to the end of viral mRNA. Both reactions are virus-specific and differ from those described for the host cell. Therefore, these capping reactions offer attractive targets for the development of antiviral drugs. Additionally, it has been shown that replication of SFV and HEV takes place in association with cellular membranes. The origin of these membranes and the intracellular localization of the components of the replication complex were studied by modern microscopy techniques. It was demonstrated that SFV replicates in cytoplasmic membranes that are derived from endosomes and lysosomes. According to our studies, site for HEV replication seems to be the intermediate compartment which mediates the traffic between endoplasmic reticulum and the Golgi complex. As a result of this work, a unique mechanism of cap formation for hepatitis E virus replicase has been characterized. It represents a novel target for the development of specific inhibitors against viral replication.

Nucleotide sequence and expression inE. coli of the complete P4 type VP4 from a G2 serotype human rotavirus

The complete sequence of a P4 type VP4 gene from a G2 serotype human rotavirus, IS2, isolated in India has been determined. Although the IS2 VP4 is highly homologous to the other P4 type alleles, it contained acidic amino acid substitutions at several positions that make it acidic among the P4 type alleles that are basic. Moreover, comparative sequence analysis revealed unusual polymorphism in members of the P4 type at amino acid position 393 which is highly conserved in members of other VP4 types. To date, expression of complete VP4 inE. coli has not been achieved. In this study we present successful expression inE. coli of the complete VP4 as well as VP8* and VP5* cleavage subunits in soluble form as fusion proteins of the maltose-binding protein (MBP) and their purification by single-step affinity chromatography. The hemagglutinating activity exhibited by the recombinant protein was specifically inhibited by the antiserum raised against it. Availability of pure VP4 proteins should facilitate development of polyclonal and monoclonal antibodies (MAbs) for P serotyping of rotaviruses.

Regulation of Osteoblast Differentiation and Gene Expression by Phosphodiesterases and Bioactive Peptides

Bone is a mineralized tissue that enables multiple mechanical and metabolic functions to be carried out in the skeleton. Bone contains distinct cell types: osteoblasts (bone-forming cells), osteocytes (mature osteoblast that embedded in mineralized bone matrix) and the osteoclasts (bone-resorbing cells). Remodelling of bone begins early in foetal life, and once the skeleton is fully formed in young adults, almost all of the metabolic activity is in this form. Bone is constantly destroyed or resorbed by osteoclasts and then replaced by osteoblasts. Many bone diseases, i.e. osteoporosis, also known as bone loss, typically reflect an imbalance in skeletal turnover. The cyclic adenosine monophosphate (cAMP) and the cyclic guanosine monophosphate (cGMP) are second messengers involved in a variety of cellular responses to such extracellular agents as hormones and neurotransmitters. In the hormonal regulation of bone metabolism, i.e. via parathyroid hormone (PTH), parathyroid hormone-related peptide (PTHrp) and prostaglandin E2 signal via cAMP. cAMP and cGMP are formed by adenylate and guanylate cyclases and are degraded by phosphodiesterases (PDEs). PDEs determine the amplitudes of cyclic nucleotide-mediated hormonal responses and modulate the duration of the signal. The activities of the PDEs are regulated by multiple inputs from other signalling systems and are crucial points of cross-talk between the pathways. Food-derived bioactive peptides are reported to express a variety of functions in vivo. The angiotensin-converting enzymes (ACEs) are involved in the regulation of the specific maturation or degradation of a number of mammalian bioactive peptides. The bioactive peptides offer also a nutriceutical and a nutrigenomic aspect to bone cell biology. The aim of this study was to investigate the influence of PDEs and bioactive peptides on the activation and the differentiation of human osteoblast cells. The profile of PDEs in human osteoblast-like cells and the effect of glucocorticoids on the function of cAMP PDEs, were investigated at the mRNA and enzyme levels. The effects of PDEs on bone formation and osteoblast gene expression were determined with chemical inhibitors and siRNAs (short interfering RNAs). The influence of bioactive peptides on osteoblast gene expression and proliferation was studied at the mRNA and cellular levels. This work provides information on how PDEs are involved in the function and the differentiation of osteoblasts. The findings illustrate that gene-specific silencing with an RNA interference (RNAi) method is useful in inhibiting, the gene expression of specific PDEs and further, PDE7 inhibition upregulates several osteogenic genes and increases bALP activity and mineralization in human mesenchymal stem cells-derived osteoblasts. PDEs appear to be involved in a mechanism by which glucocorticoids affect cAMP signaling. This may provide a potential route in the formation of glucocorticoid-induced bone loss, involving the down-regulation of cAMP-PDE. PDEs may play an important role in the regulation of osteoblastic differentiation. Isoleucine-proline-proline (IPP), a bioactive peptide, possesses the potential to increase osteoblast proliferation, differentiation and signalling.

Pheo-type: A diagnostic gene-expression assay for the classification of pheochromocytoma and paraganglioma

Context: Pheochromocytomas and paragangliomas (PPGLs) are heritable neoplasms that can be classified into gene-expression subtypes corresponding to their underlying specific genetic drivers. Objective: This study aimed to develop a diagnostic and research tool (Pheo-type) capable of classifying PPGL tumors into gene-expression subtypes that could be used to guide and interpret genetic testing, determine surveillance programs, and aid in elucidation of PPGL biology. Design: A compendium of published microarray data representing 205 PPGL tumors was used for the selection of subtype-specific genes that were then translated to the Nanostring gene-expression platform. A support vector machine was trained on the microarray dataset and then tested on an independent Nanostring dataset representing 38 familial and sporadic cases of PPGL of known genotype (RET, NF1, TMEM127, MAX, HRAS, VHL, and SDHx). Different classifier models involving between three and six subtypes were compared for their discrimination potential. Results: A gene set of 46 genes and six endogenous controls was selected representing six known PPGL subtypes; RTK1–3 (RET, NF1, TMEM127, and HRAS), MAX-like, VHL, and SDHx. Of 38 test cases, 34 (90%) were correctly predicted to six subtypes based on the known genotype to gene-expression subtype association. Removal of the RTK2 subtype from training, characterized by an admixture of tumor and normal adrenal cortex, improved the classification accuracy (35/38). Consolidation of RTK and pseudohypoxic PPGL subtypes to four- and then three-class architectures improved the classification accuracy for clinical application. Conclusions: The Pheo-type gene-expression assay is a reliable method for predicting PPGL genotype using routine diagnostic tumor samples.