Construction of a global map of human gene expression : the process, tools and analysis

This thesis studies human gene expression space using high throughput gene expression data from DNA microarrays. In molecular biology, high throughput techniques allow numerical measurements of expression of tens of thousands of genes simultaneously. In a single study, this data is traditionally obtained from a limited number of sample types with a small number of replicates. For organism-wide analysis, this data has been largely unavailable and the global structure of human transcriptome has remained unknown. This thesis introduces a human transcriptome map of different biological entities and analysis of its general structure. The map is constructed from gene expression data from the two largest public microarray data repositories, GEO and ArrayExpress. The creation of this map contributed to the development of ArrayExpress by identifying and retrofitting the previously unusable and missing data and by improving the access to its data. It also contributed to creation of several new tools for microarray data manipulation and establishment of data exchange between GEO and ArrayExpress. The data integration for the global map required creation of a new large ontology of human cell types, disease states, organism parts and cell lines. The ontology was used in a new text mining and decision tree based method for automatic conversion of human readable free text microarray data annotations into categorised format. The data comparability and minimisation of the systematic measurement errors that are characteristic to each lab- oratory in this large cross-laboratories integrated dataset, was ensured by computation of a range of microarray data quality metrics and exclusion of incomparable data. The structure of a global map of human gene expression was then explored by principal component analysis and hierarchical clustering using heuristics and help from another purpose built sample ontology. A preface and motivation to the construction and analysis of a global map of human gene expression is given by analysis of two microarray datasets of human malignant melanoma. The analysis of these sets incorporate indirect comparison of statistical methods for finding differentially expressed genes and point to the need to study gene expression on a global level.

Computational Methods for Locating and Analyzing Conserved Gene Regulatory DNA Elements

This thesis presents methods for locating and analyzing cis-regulatory DNA elements involved with the regulation of gene expression in multicellular organisms. The regulation of gene expression is carried out by the combined effort of several transcription factor proteins collectively binding the DNA on the cis-regulatory elements. Only sparse knowledge of the 'genetic code' of these elements exists today. An automatic tool for discovery of putative cis-regulatory elements could help their experimental analysis, which would result in a more detailed view of the cis-regulatory element structure and function. We have developed a computational model for the evolutionary conservation of cis-regulatory elements. The elements are modeled as evolutionarily conserved clusters of sequence-specific transcription factor binding sites. We give an efficient dynamic programming algorithm that locates the putative cis-regulatory elements and scores them according to the conservation model. A notable proportion of the high-scoring DNA sequences show transcriptional enhancer activity in transgenic mouse embryos. The conservation model includes four parameters whose optimal values are estimated with simulated annealing. With good parameter values the model discriminates well between the DNA sequences with evolutionarily conserved cis-regulatory elements and the DNA sequences that have evolved neutrally. In further inquiry, the set of highest scoring putative cis-regulatory elements were found to be sensitive to small variations in the parameter values. The statistical significance of the putative cis-regulatory elements is estimated with the Two Component Extreme Value Distribution. The p-values grade the conservation of the cis-regulatory elements above the neutral expectation. The parameter values for the distribution are estimated by simulating the neutral DNA evolution. The conservation of the transcription factor binding sites can be used in the upstream analysis of regulatory interactions. This approach may provide mechanistic insight to the transcription level data from, e.g., microarray experiments. Here we give a method to predict shared transcriptional regulators for a set of co-expressed genes. The EEL (Enhancer Element Locator) software implements the method for locating putative cis-regulatory elements. The software facilitates both interactive use and distributed batch processing. We have used it to analyze the non-coding regions around all human genes with respect to the orthologous regions in various other species including mouse. The data from these genome-wide analyzes is stored in a relational database which is used in the publicly available web services for upstream analysis and visualization of the putative cis-regulatory elements in the human genome.

Novel CDNF/MANF protein family : Molecular structure, expression and neurotrophic activity

Neurotrophic factors (NTFs) are secreted proteins which promote the survival of neurons, formation and maintenance of neuronal contacts and regulate synaptic plasticity. NTFs are also potential drug candidates for the treatment of neurodegenerative diseases. Parkinson’s disease (PD) is mainly caused by the degeneration of midbrain dopaminergic neurons. Current therapies for PD do not stop the neurodegeneration or repair the affected neurons. Thus, search of novel neurotrophic factors for midbrain dopaminergic neurons, which could also be used as therapeutic proteins, is highly warranted. In the present study, we identified and characterized a novel protein named conserved dopamine neurotrophic factor (CDNF), a homologous protein to mesencephalic astrocyte-derived neurotrophic factor (MANF). Others have shown that MANF supports the survival of embryonic midbrain dopaminergic neurons in vitro, and protects cultured cells against endoplasmic reticulum (ER) stress. CDNF and MANF form a novel evolutionary conserved protein family with characteristic eight conserved cysteine residues in their primary structure. The vertebrates have CDNF and MANF encoding genes, whereas the invertebrates, including Drosophila and Caenorhabditis have a single homologous CDNF/MANF gene. In this study we show that CDNF and MANF are secreted proteins. They are widely expressed in the mammalian brain, including the midbrain and striatum, and in several non-neuronal tissues. We expressed and purified recombinant human CDNF and MANF proteins, and tested the neurotrophic activity of CDNF on midbrain dopaminergic neurons using a 6-hydroxydopamine (6-OHDA) rat model of PD. In this model, a single intrastriatal injection of CDNF protected midbrain dopaminergic neurons and striatal dopaminergic fibers from the 6-OHDA toxicity. Importantly, an intrastriatal injection of CDNF also restored the functional activity of the nigrostriatal dopaminergic system when given after the striatal 6-OHDA lesion. Thus, our study shows that CDNF is a potential novel therapeutic protein for the treatment of PD. In order to elucidate the molecular mechanisms of CDNF and MANF activity, we resolved their crystal structure. CDNF and MANF proteins have two domains; an amino (N)-terminal saposin-like domain and a presumably unfolded carboxy (C)-terminal domain. The saposin-like domain, which is formed by five α-helices and stabilized by three intradomain disulphide bridges, may bind to lipids or membranes. The C-terminal domain contains an internal cysteine bridge in a CXXC motif similar to that of thiol/disulphide oxidoreductases and isomerases, and may thus facilitate protein folding in the ER. Our studies suggest that CDNF and MANF are novel potential therapeutic proteins for the treatment of neurodegenerative diseases. Future studies will reveal the neurotrophic and cytoprotective mechanisms of CDNF and MANF in more detail.

Structure and function of GDNF receptor alpha splice variants

The growth factors of the glial cell line-derived neurotrophic factor (GDNF) family consisting of GDNF, neurturin (NRTN), artemin (ARTN) and persephin (PSPN), are involved in the development, differentiation and maintenance of many types of neurons. They also have important functions outside the nervous system in the development of kidney, testis and thyroid gland. Each of these GFLs preferentially binds to one of the glycosylphosphatidylinositol (GPI)-anchored GDNF family receptors α (GFRα). GDNF binds to GFRα1, NRTN to GFRα2, ARTN to GFRα3 and PSPN to GFRα4. The GFLs in the complex with their cognate GFRα receptors all bind to and signal through the receptor tyrosine kinase RET. Alternative splicing of the mouse GFRα4 gene yields three splice isoforms. These had been described as putative GPI-anchored, transmembrane and soluble forms. My goal was to characterise the function of the different forms of mouse GFRα4. I firstly found that the putative GPI-anchored GFRα4 (GFRα4-GPI) is glycosylated, membrane-bound, GPI-anchored and interacts with PSPN and RET. We also showed that mouse GFRα4-GPI mediates PSPN-induced phosphorylation of RET, promotes PSPN-dependent neuronal differentiation of the rat pheochromocytoma cell line PC6-3 and PSPN-dependent survival of cerebellar granule neurons (CGN). However, although this receptor can mediate PSPN-signalling and activate RET, GFRα4-GPI does not recruit RET into lipid rafts. The recruitment of RET into lipid rafts has previously been thought to be a crucial event for GDNF- and GFL-mediated signalling via RET. I secondly demonstrated that the putative transmembrane GFRα4 (GFRα4-TM) is indeed a real transmembrane GFRα4 protein. Although it has a weak binding capacity for PSPN, it can not mediate PSPN-dependent phosphorylation of RET, neuronal differentiation or survival. These data show that GFRα4-TM is inactive as a receptor for PSPN. Surprisingly, GFRα4-TM can negatively regulate PSPN-mediated signalling via GFRα4-GPI. GFRα4-TM interacts with GFRα4-GPI and blocks PSPN-induced phosphorylation of RET, neuronal differentiation as well as survival. Taken together, our data show that GFRα4-TM may act as a dominant negative inhibitor of PSPN-mediated signaling. The most exciting part of my work was the finding that the putative soluble GFRα4 (GFRα4-sol) can form homodimers and function as an agonist of the RET receptor. In the absence of PSPN, GFRα4-sol can promote the phosphorylation of RET, trigger the activation of the PI-3K/AKT pathway, induce neuronal differentiation and support the survival of CGN. Our findings are in line with a recent publication showing the GFRα4-sol might contribute to the inherited cancer syndrome multiple endocrine neoplasia type 2. Our data provide an explanation to how GFRα4-sol may cause or modify the disease. Mammalian GFRα4 receptors all lack the first Cys-rich domain which is present in other GFRα receptors. In the final part of my work I have studied the function of this particular domain. I created a truncated GFRα1 construct lacking the first Cys-rich domain. Using binding assays in both cellular and cell-free systems, phosphorylation assays with RET, as well as neurite outgrowth assays, we found that the first Cys-rich domain contributes to an optimal function of GFRα1, by stabilizing the interaction between GDNF and GFRα1.

Reproduction and population structure in European aspen

The European aspen (Populus tremula) is a keystone species for biodiversity in boreal forests. However, the future of aspen may be threatened, because large aspens have mostly been removed from managed forests, whereas regeneration and the long-term persistence of mature trees are subjects of concern in protected areas. Aspen is a pioneer tree, and it can reproduce both sexually by seed and asexually by root suckers. Through asexual reproduction aspen forms clones, groups of genetically identical trees (ramets). In my thesis, I have studied the structure of aspen populations in terms of number, size, clonal and demographic properties. Additionally, I have investigated the emergence and survival of seedlings as well as the seed quantity and quality in crosses between the European and hybrid aspen. To study the regeneration and population structure, mature aspens were recorded in old-growth and managed forests in eastern Finland based on a large-scale inventory (11 400 ha). In addition, small aspen trees were surveyed on sample plots. Clonal structure was investigated both by morphological characters and by DNA-based markers (microsatellites). Seedling emergence and survival was studied with two sowing experiments. With crosses between European and hybrid aspens we wanted to study whether elevated temperatures due to climate change would benefit the different crosses of European and hybrid aspen unequally and thus affect the gene flow between the two species. The average volumes of mature aspen were 5.3 m3/ha in continuous old-growth, and 0.8 m3/ha in managed forests. Results indicate also that large aspen trees in managed forests are a legacy of the past less intensively managed forest landscapes. Long-term persistence of aspen in protected areas can only be secured by restoration measures creating sufficiently large gaps for regeneration. More emphasis should be given to sparing aspens in thinnings and to retaining of mature aspens in regeneration cutting in managed forests. Aspen was found to be spatially aggregated in the landscape. This could be explained by site type, disturbance history and / or limitations in seed dispersal. Clonal structure does not explain the spatial aggregation, since average size of the clones was only 2.3 ramets, and most clones (70 %) consisted of just one ramet. The small size of the clones suggests that most of them are relatively young. Therefore, sexual reproduction may be more common than has previously been thought. Seedling emergence was most successful in mineral soil especially, when the site had been burned. Only few seedlings occurred on humus. Survival of the seedlings was low, and strongly dependent on moisture, but also on seedbed conditions. The seeds were found to maintain their germinability longer than has earlier been thought to be possible. Interspecific crosses produced more seeds with higher quality than intraspecific crosses. When temperature was elevated, germination of hybrid aspen seeds increased more than seeds from P. tremula x P. tremula crosses. These results suggest that hybrid aspen may have a significant genetic impact on the European aspen, and this effect may become strengthened by climate warming.