Evolutionary Genomics of Prokaryotic Viruses

Evolutionary history of biological entities is recorded within their nucleic acid sequences and can (sometimes) be deciphered by thorough genomic analysis. In this study we sought to gain insights into the diversity and evolution of bacterial and archaeal viruses. Our primary interest was pointed towards those virus groups/families for which comprehensive genomic analysis was not previously possible due to the lack of sufficient amount of genomic data. During the course of this work twenty-five putative proviruses integrated into various prokaryotic genomes were identified, enabling us to undertake a comparative genomics approach. This analysis allowed us to test the previously formulated evolutionary hypotheses and also provided valuable information on the molecular mechanisms behind the genome evolution of the studied virus groups.

Mu in vitro Transposition Technology in Functional Genetics and Genomics: Applications on Mouse and Bacteriophages

Transposons are mobile elements of genetic material that are able to move in the genomes of their host organisms using a special form of recombination called transposition. Bacteriophage Mu was the first transposon for which a cell-free in vitro transposition reaction was developed. Subsequently, the reaction has been refined and the minimal Mu in vitro reaction is useful in the generation of comprehensive libraries of mutant DNA molecules that can be used in a variety of applications. To date, the functional genetics applications of Mu in vitro technology have been subjected to either plasmids or genomic regions and entire genomes of viruses cloned on specific vectors. This study expands the use of Mu in vitro transposition in functional genetics and genomics by describing novel methods applicable to the targeted transgenesis of mouse and the whole-genome analysis of bacteriophages. The methods described here are rapid, efficient, and easily applicable to a wide variety of organisms, demonstrating the potential of the Mu transposition technology in the functional analysis of genes and genomes. First, an easy-to-use, rapid strategy to generate construct for the targeted mutagenesis of mouse genes was developed. To test the strategy, a gene encoding a neuronal K+/Cl- cotransporter was mutagenised. After a highly efficient transpositional mutagenesis, the gene fragments mutagenised were cloned into a vector backbone and transferred into bacterial cells. These constructs were screened with PCR using an effective 3D matrix system. In addition to traditional knock-out constructs, the method developed yields hypomorphic alleles that lead into reduced expression of the target gene in transgenic mice and have since been used in a follow-up study. Moreover, a scheme is devised to rapidly produce conditional alleles from the constructs produced. Next, an efficient strategy for the whole-genome analysis of bacteriophages was developed based on the transpositional mutagenesis of uncloned, infective virus genomes and their subsequent transfer into susceptible host cells. Mutant viruses able to produce viable progeny were collected and their transposon integration sites determined to map genomic regions nonessential to the viral life cycle. This method, applied here to three very different bacteriophages, PRD1, ΦYeO3 12, and PM2, does not require the target genome to be cloned and is directly applicable to all DNA and RNA viruses that have infective genomes. The method developed yielded valuable novel information on the three bacteriophages studied and whole-genome data can be complemented with concomitant studies on individual genes. Moreover, end-modified transposons constructed for this study can be used to manipulate genomes devoid of suitable restriction sites.

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.

Integrative functional genomics analysis of sustained polyploidy phenotypes in breast cancer cells identifies an oncogenic profile for GINS2

Aneuploidy is among the most obvious differences between normal and cancer cells. However, mechanisms contributing to development and maintenance of aneuploid cell growth are diverse and incompletely understood. Functional genomics analyses have shown that aneuploidy in cancer cells is correlated with diffuse gene expression signatures and that aneuploidy can arise by a variety of mechanisms, including cytokinesis failures, DNA endoreplication and possibly through polyploid intermediate states. Here, we used a novel cell spot microarray technique to identify genes with a loss-of-function effect inducing polyploidy and/or allowing maintenance of polyploid cell growth of breast cancer cells. Integrative genomics profiling of candidate genes highlighted GINS2 as a potential oncogene frequently overexpressed in clinical breast cancers as well as in several other cancer types. Multivariate analysis indicated GINS2 to be an independent prognostic factor for breast cancer outcome (p = 0.001). Suppression of GINS2 expression effectively inhibited breast cancer cell growth and induced polyploidy. In addition, protein level detection of nuclear GINS2 accurately distinguished actively proliferating cancer cells suggesting potential use as an operational biomarker.

Mitochondrial Malfunction in Tumor Formation and Neuromuscular Diseases : Consequences of defects in fumarate hydratase and in the respiratory chain

The mitochondrion is an organelle of outmost importance, and the mitochondrial network performs an array of functions that go well beyond ATP synthesis. Defects in mitochondrial performance lead to diseases, often affecting nervous system and muscle. Although many of these mitochondrial diseases have been linked to defects in specific genes, the molecular mechanisms underlying the pathologies remain unclear. The work in this thesis aims to determine how defects in mitochondria are communicated within - and interpreted by - the cells, and how this contributes to disease phenotypes. Fumarate hydratase (FH) is an enzyme of the citrate cycle. Recessive defects in FH lead to infantile mitochondrial encephalopathies, while dominant mutations predispose to tumor formation. Defects in succinate dehydrogenase (SDH), the enzyme that precedes FH in the citrate cycle, have also been described. Mutations in SDH subunits SDHB, SDHC and SDHD are associated with tumor predisposition, while mutations in SDHA lead to a characteristic mitochondrial encephalopathy of childhood. Thus, the citrate cycle, via FH and SDH, seems to have essential roles in mitochondrial function, as well as in the regulation of processes such as cell proliferation, differentiation or death. Tumor predisposition is not a typical feature of mitochondrial energy deficiency diseases. However, defects in citrate cycle enzymes also affect mitochondrial energy metabolism. It is therefore necessary to distinguish what is specific for defects in citrate cycle, and thus possibly associated with the tumor phenotype, from the generic consequences of defects in mitochondrial aerobic metabolism. We used primary fibroblasts from patients with recessive FH defects to study the cellular consequences of FH-deficiency (FH-). Similarly to the tumors observed in FH- patients, these fibroblasts have very low FH activity. The use of primary cells has the advantage that they are diploid, in contrast with the aneuploid tumor cells, thereby enabling the study of the early consequences of FH- in diploid background, before tumorigenesis and aneuploidy. To distinguish the specific consequences of FH- from typical consequences of defects in mitochondrial aerobic metabolism, we used primary fibroblasts from patients with MELAS (mitochondrial encephalopathy with lactic acidosis and stroke-like episodes) and from patients with NARP (neuropathy, ataxia and retinitis pigmentosa). These diseases also affect mitochondrial aerobic metabolism but are not known to predispose to tumor formation. To study in vivo the systemic consequences of defects in mitochondrial aerobic metabolism, we used a transgenic mouse model of late-onset mitochondrial myopathy. The mouse contains a transgene with an in-frame duplication of a segment of Twinkle, the mitochondrial replicative helicase, whose defects underlie the human disease progressive external ophthalmoplegia. This mouse model replicates the phenotype in the patients, particularly neuronal degeneration, mitochondrial myopathy, and subtle decrease of respiratory chain activity associated with mtDNA deletions. Due to the accumulation of mtDNA deletions, the mouse was named deletor. We first studied the consequences of FH- and of respiratory chain defects for energy metabolism in primary fibroblasts. To further characterize the effects of FH- and respiratory chain malfunction in primary fibroblasts at transcriptional level, we used expression microarrays. In order to understand the in vivo consequences of respiratory chain defects in vivo, we also studied the transcriptional consequences of Twinkle defects in deletor mice skeletal muscle, cerebellum and hippocampus. Fumarate accumulated in the FH- homozygous cells, but not in the compound heterozygous lines. However, virtually all FH- lines lacked cytoplasmic FH. Induction of glycolysis was common to FH-, MELAS and NARP fibroblasts. In deletor muscle glycolysis seemed to be upregulated. This was in contrast with deletor cerebellum and hippocampus, where mitochondrial biogenesis was in progress. Despite sharing a glycolytic pattern in energy metabolism, FH- and respiratory chain defects led to opposite consequences in redox environment. FH- was associated with reduced redox environment, while MELAS and NARP displayed evidences of oxidative stress. The deletor cerebellum had transcriptional induction of antioxidant defenses, suggesting increased production of reactive oxygen species. Since the fibroblasts do not represent the tissues where the tumors appear in FH- patients, we compared the fibroblast array data with the data from FH- leiomyomas and normal myometrium. This allowed the determination of the pathways and networks affected by FH-deficiency in primary cells that are also relevant for myoma formation. A key pathway regulating smooth muscle differentiation, SRF (serum response factor)-FOS-JUNB, was found to be downregulated in FH- cells and in myomas. While in the deletor mouse many pathways were affected in a tissue-specific basis, like FGF21 induction in the deletor muscle, others were systemic, such as the downregulation of ALAS2-linked heme synthesis in all deletor tissues analyzed. However, interestingly, even a tissue-specific response of FGF21 excretion could elicit a global starvation response. The work presented in this thesis has contributed to a better understanding of mitochondrial stress signalling and of pathways interpreting and transducing it to human pathology.

Animal model and molecular interactions of Cln5

Neuronal ceroid lipofuscinoses (NCLs) are a family of inherited pediatric neurodegenerative disorders, leading to retinal degeneration, death of selective neuronal populations and accumulation of autofluorscent ceroid-lipopigments. The clinical manifestations are generally similar in all forms. The Finnish variant late infantile neuronal ceroid lipofuscinosis (vLINCLFin) is a form of NCL, especially enriched in the Finnish population. The aim of this thesis was to analyse the brain pathology of vLINCLFin utilising the novel Cln5-/- mouse model. Gene expression profiling of the brains of already symptomatic Cln5-/- mice revealed that inflammation, neurodegeneration and defects in myelinization are the major characteristics of the later stages of the disease. Histological characterization of the brain pathology confirmed that the thalamocortical system is affected in Cln5-/- mice, similarly to the other NCL mouse models. However, whereas the brain pathology in all other analyzed NCL mice initiate in the thalamus and spread only months later to the cortex, we observed that the sequence of events is uniquely reversed in Cln5-/- mice; beginning in the cortex and spreading to the thalamus only months later. We could also show that even though neurodegeneration is inititated in the cortex, reactive gliosis and loss of myelin are evident in specific nuclei of the thalamus already in the 1 month old brain. To obtain a deeper insight into the disturbed metabolic pathways, we performed gene expression profiling of presymptomatic mouse brains. We validated these findings with immunohistological analyses, and could show that cytoskeleton and myelin were affected in Cln5-/- mice. Comparison of gene expression profiling results of Cln5-/- and Cln1-/- mice, further highlighted that these two NCL models share a common defective pathway, leading to disturbances in the neuronal growth cone and cytoskeleton. Encouraged by the evidence of this defected pathway, we analyzed the molecular interactions of NCL-proteins and observed that Cln5 and Cln1/Ppt1 proteins interact with each other. Furthermore, we demonstrated that Cln5 and Cln1/Ppt1 share an interaction partner, the F1-ATP synthase, potentially linking both vLINCLFIN and INCL diseases to disturbed lipid metabolism. In addition, Cln5 was shown to interact with other NCL proteins; Cln2, Cln3, Cln6 and Cln8, implicating a central role for Cln5 in the NCL pathophysiology. This study is the first to describe the brain pathology and gene expression changes in the Cln5-/- mouse. Together the findings presented in this thesis represent novel information of the disease processes and the molecular mechanisms behind vLINCLFin and have highlighted that vLINCLFin forms a very important model to analyze the pathophysiology of NCL diseases.

Molecular Pathways of Disturbed Sleep and Depression: Studies on Adenosine and Gene Expression Patterns

Background: Adenosine is a potent sleep-promoting substance, and one of its targets is the basal forebrain. Fairly little is known about its mechanism of action in the basal forebrain and about the receptor subtype mediating its regulating effects on sleep homeostasis. Homeostatic deficiency might be one of the causes of the profoundly disturbed sleep pattern in major depressive disorder, which could explain the reduced amounts of delta-activity-rich stages 3 and 4. Since major depression has a relatively high heritability, and on the other hand adenosine regulates sleep homeostasis and might also be involved in mood modulation, adenosine-related genes should be considered for their possible contribution to a predisposition for depression and disturbed sleep in humans. Depression is a complex disorder likely involving the abnormal functioning of several genes. Novel target genes which could serve as the possible common substrates for depression and comorbid disturbed sleep should be identified. In this way specific brain areas related to sleep regulation should be studied by using animal model of depression which represents more homogenous phenotype as compared to humans. It is also important to study these brain areas during the development of depressive-like features to understand how early changes could facilitate pathophysiological changes in depression. Aims and methods: We aimed to find out whether, in the basal forebrain, adenosine induces recovery non-rapid eye movement (NREM) sleep after prolonged waking through the A1 or/and A2A receptor subtype. A1 and A2A receptor antagonists were perfused into the rat basal forebrain during 3 h of sleep deprivation, and the amount of NREM sleep and delta power during recovery NREM sleep were analyzed. We then explored whether polymorphisms in genes related to the metabolism, transport and signaling of adenosine could predispose to depression accompanied by signs of disturbed sleep. DNA from 1423 individuals representative of the Finnish population and including controls and cases with depression, depression accompanied by early morning awakenings and depression accompanied by fatigue, was used in the study to investigate the possible association between polymorphisms from adenosine-related genes and cases. Finally to find common molecular substrates of depression and disturbed sleep, gene expression changes were investigated in specific brain areas in the rat clomipramine model of depression. We focused on the basal forebrain of 3-week old clomipramine-treated rats which develop depressive-like symptoms later in adulthood and on the hypothalamus of adult female clomipramine-treated rats. Results: Blocking of the A1 receptor during sleep deprivation resulted in a reduction of the recovery NREM sleep amount and delta power, whereas A2A receptor antagonism had no effect. Polymorphisms in adenosine-related genes SLC29A3 (equilibrative nucleoside transporter type 3) in women and SLC28A1 (concentrative nucleoside transporter type 1) in men associated with depression alone as well as when accompanied by early morning awakenings and fatigue. In Study III the basal forebrain of postnatal rats treated with clomipramine displayed disturbances in gamma-aminobutyric acid (GABA) receptor type A signaling, in synaptic transmission and possible epigenetic changes. CREB1 was identified as a common transcription denominator which also mediates epigenetic regulation. In the hypothalamus the major changes included the expression of genes in GABA-A receptor pathway, K+ channel-related, glutamatergic and mitochondrial genes, as well as an overexpression of genes related to RNA and mRNA processing. Conclusions: Adenosine plays an important role in sleep homeostasis by promoting recovery NREM sleep via the A1 receptor subtype in the basal forebrain. Also adenosine levels might contribute to the risk of depression with disturbed sleep, since the genes encoding nucleoside transporters showed the strongest associations with depression alone and when accompanied by signs of disturbed sleep in both women and men. Sleep and mood abnormalities in major depressive disorder could be a consequence of multiple changes at the transcriptional level, GABA-A receptor signaling and synaptic transmission in sleep-related basal forebrain and the hypothalamus.

Methods to improve gene signal: Application to cDNA microarrays

Microarrays are high throughput biological assays that allow the screening of thousands of genes for their expression. The main idea behind microarrays is to compute for each gene a unique signal that is directly proportional to the quantity of mRNA that was hybridized on the chip. A large number of steps and errors associated with each step make the generated expression signal noisy. As a result, microarray data need to be carefully pre-processed before their analysis can be assumed to lead to reliable and biologically relevant conclusions. This thesis focuses on developing methods for improving gene signal and further utilizing this improved signal for higher level analysis. To achieve this, first, approaches for designing microarray experiments using various optimality criteria, considering both biological and technical replicates, are described. A carefully designed experiment leads to signal with low noise, as the effect of unwanted variations is minimized and the precision of the estimates of the parameters of interest are maximized. Second, a system for improving the gene signal by using three scans at varying scanner sensitivities is developed. A novel Bayesian latent intensity model is then applied on these three sets of expression values, corresponding to the three scans, to estimate the suitably calibrated true signal of genes. Third, a novel image segmentation approach that segregates the fluorescent signal from the undesired noise is developed using an additional dye, SYBR green RNA II. This technique helped in identifying signal only with respect to the hybridized DNA, and signal corresponding to dust, scratch, spilling of dye, and other noises, are avoided. Fourth, an integrated statistical model is developed, where signal correction, systematic array effects, dye effects, and differential expression, are modelled jointly as opposed to a sequential application of several methods of analysis. The methods described in here have been tested only for cDNA microarrays, but can also, with some modifications, be applied to other high-throughput technologies. Keywords: High-throughput technology, microarray, cDNA, multiple scans, Bayesian hierarchical models, image analysis, experimental design, MCMC, WinBUGS.

Endoplasmic reticulum stress and cell death mechanisms in the cell model of Huntington’s disease

This thesis clarifies important molecular pathways that are activated during the cell death observed in Huntington’s disease. Huntington’s disease is one of the most common inherited neurodegenerative diseases, which is primarily inherited in an autosomal dominant manner. HD is caused by an expansion of CAG repeats in the first exon of the IT15 gene. IT15 encodes the production of a Huntington’s disease protein huntingtin. Mutation of the IT15 gene results in a long stretch of polyQ residues close to the amino-terminal region of huntingtin. Huntington’s disease is a fatal autosomal neurodegenerative disorder. Despite the current knowledge of HD, the precise mechanism behind the selective neuronal death, and how the disease propagates, still remains an enigma. The studies mainly focused on the control of endoplasmic reticulum (ER) stress triggered by the mutant huntingtin proteins. The ER is a delicate organelle having essential roles in protein folding and calcium regulation. Even the slightest perturbations on ER homeostasis are effective enough to trigger ER stress and its adaptation pathways, called unfolded protein response (UPR). UPR is essential for cellular homeostasis and it adapts ER to the changing environment and decreases ER stress. If adaptation processes fail and stress is excessive and prolonged; irreversible cell death pathways are engaged. The results showed that inhibition of ER stress with chemical agents are able to decrease cell death and formation of toxic cell aggregates caused by mutant huntingtin proteins. The study concentrated also to the NF-κB (nuclear factor-kappaB) pathway, which is activated during ER stress. NF-κB pathway is capable to regulate the levels of important cellular antioxidants. Cellular antioxidants provide a first line of defence against excess reactive oxygen species. Excess accumulation of reactive oxygen species and subsequent activation of oxidative stress damages motley of vital cellular processes and induce cell degeneration. Data showed that mutant huntingtin proteins downregulate the expression levels of NF-κB and vital antioxidants, which was followed by increased oxidative stress and cell death. Treatment with antioxidants and inhibition of oxidative stress were able to counteract these adverse effects. In addition, thesis connects ER stress caused by mutant huntingtin to the cytoprotective autophagy. Autophagy sustains cellular balance by degrading potentially toxic cell proteins and components observed in Huntington’s disease. The results revealed that cytoprotective autophagy is active at the early points (24h) of ER stress after expression of mutant huntingtin proteins. GADD34 (growth arrest and DNA damage-inducible gene 34), which is previously connected to the regulation of translation during cell stress, was shown to control the stimulation of autophagy. However, GADD34 and autophagy were downregulated at later time points (48h) during mutant huntingtin proteins induced ER stress, and subsequently cell survival decreased. Overexpression GADD34 enhanced autophagy and decreased cell death, indicating that GADD34 plays a critical role in cell protection. The thesis reveales new interesting data about the neuronal cell death pathways seen in Huntington’s disease, and how cell degeneration is partly counteracted by various therapeutic agents. Expression of mutant huntingtin proteins is shown to alter signaling events that control ER stress, oxidative stress and autophagy. Despite that Huntington’s disease is mainly an untreatable disorder; these findings offer potential targets and neuroprotective strategies in designing novel therapies for Huntington’s disease.

ORP1L, a new Rab7 effector and regulator of late endosome functions

Oxysterol binding protein (OSBP) homologues have been found in eukaryotic organisms ranging from yeast to humans. These evolutionary conserved proteins have in common the presence of an OSBP-related domain (ORD) which contains the fully conserved EQVSHHPP sequence motif. The ORD forms a barrel structure that binds sterols in its interior. Other domains and sequence elements found in OSBP-homologues include pleckstrin homology domains, ankyrin repeats and two phenylalanines in an acidic tract (FFAT) motifs, which target the proteins to distinct subcellular compartments. OSBP homologues have been implicated in a wide range of intracellular processes, including vesicle trafficking, lipid metabolism and cell signaling, but little is known about the functional mechanisms of these proteins. The human family of OSBP homologues consists of twelve OSBP-related proteins (ORP). This thesis work is focused on one of the family members, ORP1, of which two variants were found to be expressed tissue-specifically in humans. The shorter variant, ORP1S contains an ORD only. The N-terminally extended variant, ORP1L, comprises a pleckstrin homology domain and three ankyrin repeats in addition to the ORD. The two ORP1 variants differ in intracellular localization. ORP1S is cytosolic, while the ankyrin repeat region of ORP1L targets the protein to late endosomes/lysosomes. This part of ORP1L also has profound effects on late endosomal morphology, inducing perinuclear clustering of late endosomes. A central aim of this study was to identify molecular interactions of ORP1L on late endosomes. The morphological changes of late endosomes induced by overexpressed ORP1L implies involvement of small Rab GTPases, regulators of organelle motility, tethering, docking and/or fusion, in generation of the phenotype. A direct interaction was demonstrated between ORP1L and active Rab7. ORP1L prolongs the active state of Rab7 by stabilizing its GTP-bound form. The clustering of late endosomes/lysosomes was also shown to be linked to the minus end-directed microtubule-based dynein-dynactin motor complex through the ankyrin repeat region of ORP1L. ORP1L, Rab7 and the Rab7-interacting lysosomal protein (RILP) were found to be part of the same effector complex recruiting the dynein-dynactin complex to late endosomes, thereby promoting minus end-directed movement. The proteins were found to be physically close to each other on late endosomes and RILP was found to stabilize the ORP1L-Rab7 interaction. It is possible that ORP1L and RILP bind to each other through their C-terminal and N-terminal regions, respectively, when they are bridged by Rab7. With the results of this study we have been able to place a member of the uncharacterized OSBP-family, ORP1L, in the endocytic pathway, where it regulates motility and possibly fusion of late endosomes through interaction with the small GTPase Rab7.

The Gerbera cDNA Microarray: A Tool for Large-Scale Identification of Genes Involved in Flower Development

During the past ten years, large-scale transcript analysis using microarrays has become a powerful tool to identify and predict functions for new genes. It allows simultaneous monitoring of the expression of thousands of genes and has become a routinely used tool in laboratories worldwide. Microarray analysis will, together with other functional genomics tools, take us closer to understanding the functions of all genes in genomes of living organisms. Flower development is a genetically regulated process which has mostly been studied in the traditional model species Arabidopsis thaliana, Antirrhinum majus and Petunia hybrida. The molecular mechanisms behind flower development in them are partly applicable in other plant systems. However, not all biological phenomena can be approached with just a few model systems. In order to understand and apply the knowledge to ecologically and economically important plants, other species also need to be studied. Sequencing of 17 000 ESTs from nine different cDNA libraries of the ornamental plant Gerbera hybrida made it possible to construct a cDNA microarray with 9000 probes. The probes of the microarray represent all different ESTs in the database. From the gerbera ESTs 20% were unique to gerbera while 373 were specific to the Asteraceae family of flowering plants. Gerbera has composite inflorescences with three different types of flowers that vary from each other morphologically. The marginal ray flowers are large, often pigmented and female, while the central disc flowers are smaller and more radially symmetrical perfect flowers. Intermediate trans flowers are similar to ray flowers but smaller in size. This feature together with the molecular tools applied to gerbera, make gerbera a unique system in comparison to the common model plants with only a single kind of flowers in their inflorescence. In the first part of this thesis, conditions for gerbera microarray analysis were optimised including experimental design, sample preparation and hybridization, as well as data analysis and verification. Moreover, in the first study, the flower and flower organ-specific genes were identified. After the reliability and reproducibility of the method were confirmed, the microarrays were utilized to investigate transcriptional differences between ray and disc flowers. This study revealed novel information about the morphological development as well as the transcriptional regulation of early stages of development in various flower types of gerbera. The most interesting finding was differential expression of MADS-box genes, suggesting the existence of flower type-specific regulatory complexes in the specification of different types of flowers. The gerbera microarray was further used to profile changes in expression during petal development. Gerbera ray flower petals are large, which makes them an ideal model to study organogenesis. Six different stages were compared and specifically analysed. Expression profiles of genes related to cell structure and growth implied that during stage two, cells divide, a process which is marked by expression of histones, cyclins and tubulins. Stage 4 was found to be a transition stage between cell division and expansion and by stage 6 cells had stopped division and instead underwent expansion. Interestingly, at the last analysed stage, stage 9, when cells did not grow any more, the highest number of upregulated genes was detected. The gerbera microarray is a fully-functioning tool for large-scale studies of flower development and correlation with real-time RT-PCR results show that it is also highly sensitive and reliable. Gene expression data presented here will be a source for gene expression mining or marker gene discovery in the future studies that will be performed in the Gerbera Laboratory. The publicly available data will also serve the plant research community world-wide.

Golgi proteomics : Identification of a novel cartilage-specific Golgi protein GoPro49

The Golgi complex is a central organelle of the secretory pathway, responsible for a range of post-translational modifications, as well as for membrane traffic to the plasma membrane and to the endosomal-lysosomal pathway. In addition, this organelle has roles in cell migration, in the regulation of traffic, and as a mitotic check point. The structure of the Golgi complex is highly dynamic and able to respond to the amount of cargo being transported and the stage of the cell cycle. The Golgi proteome reflects the functions and structure of this organelle, and can be divided into three major groups: the Golgi resident proteins (e.g. modification enzymes), the Golgi matrix proteins (involved in structure and tethering events), and trafficking proteins (e.g. vesicle coat proteins and Rabs). The Golgi proteome has been studied on several occasions, from both rat liver and mammary gland Golgi membranes using proteomic approaches, but still little more than half of the estimated Golgi proteome is known. Nevertheless, methodological improvements and introduction of shotgun proteomics have increased the number of identified proteins, and especially the number of identified transmembrane proteins. Cartilage, even though not a typical tissue in which to study membrane traffic, secretes large amounts of extracellular matrix proteins that are extensively modified, especially by amino acid hydroxylation, glycosylation and sulfation. Furthermore, the cartilage ECM contains several, large oligomeric proteins (such as collagen II) that are difficult to assemble and transport. Indeed, cartilage has been shown to be susceptible to changes both in secretory pathway (e.g. the COPII coat assembly) and in post-translational modifications (e.g. heparan sulfate formation). Dental follicle, and the periodontal ligament (PDL) that it forms, are another type of connective tissue, and they have a role in anchoring teeth to bone. This anchorage is achieved by numerous matrix fibres that connect the bone matrix with the cementum. These tissues have in common the secretion of large matrix molecules. In this study the Golgi proteome was analysed from purified, stacked Golgi membranes isolated from rat liver. The identified, extensive proteome included a protein similar to Ab2-095, or Golgi protein 49kDa (GoPro49), which was shown to localise to the Golgi complex as an EGFP fusion protein. Surprisingly, in situ hybridisation showed the GoPro49 expression to be highly restricted to different mesenchymal tissues, especially in cartilage, and this expression pattern was clearly developmentally regulated. In addition to cartilage, GoPro49 was also expressed in the dental follicle, but was not observed in the mature PDL. Importantly, GoPro49 is the first specific marker for the dental follicle. Endogenous GoPro49 protein co-localised with β-COP in both chondrosarcoma and primary dental follicle cell lines. The COPI staining in these cells was highly dynamic, showing a number of tubules. This may reflect the type of secretory cargo they secrete. Currently GoPro49 is the only Golgi protein with such a restricted expression pattern.

Comparative and functional genome analysis of fungi for development of the protein production host Trichoderma reesei

Filamentous fungi of the subphylum Pezizomycotina are well known as protein and secondary metabolite producers. Various industries take advantage of these capabilities. However, the molecular biology of yeasts, i.e. Saccharomycotina and especially that of Saccharomyces cerevisiae, the baker's yeast, is much better known. In an effort to explain fungal phenotypes through their genotypes we have compared protein coding gene contents of Pezizomycotina and Saccharomycotina. Only biomass degradation and secondary metabolism related protein families seem to have expanded recently in Pezizomycotina. Of the protein families clearly diverged between Pezizomycotina and Saccharomycotina, those related to mitochondrial functions emerge as the most prominent. However, the primary metabolism as described in S. cerevisiae is largely conserved in all fungi. Apart from the known secondary metabolism, Pezizomycotina have pathways that could link secondary metabolism to primary metabolism and a wealth of undescribed enzymes. Previous studies of individual Pezizomycotina genomes have shown that regardless of the difference in production efficiency and diversity of secreted proteins, the content of the known secretion machinery genes in Pezizomycotina and Saccharomycotina appears very similar. Genome wide analysis of gene products is therefore needed to better understand the efficient secretion of Pezizomycotina. We have developed methods applicable to transcriptome analysis of non-sequenced organisms. TRAC (Transcriptional profiling with the aid of affinity capture) has been previously developed at VTT for fast, focused transcription analysis. We introduce a version of TRAC that allows more powerful signal amplification and multiplexing. We also present computational optimisations of transcriptome analysis of non-sequenced organism and TRAC analysis in general. Trichoderma reesei is one of the most commonly used Pezizomycotina in the protein production industry. In order to understand its secretion system better and find clues for improvement of its industrial performance, we have analysed its transcriptomic response to protein secretion stress conditions. In comparison to S. cerevisiae, the response of T. reesei appears different, but still impacts on the same cellular functions. We also discovered in T. reesei interesting similarities to mammalian protein secretion stress response. Together these findings highlight targets for more detailed studies.

Genetic Mechanisms Underlying Autism Spectrum Disorders

Autism is a childhood-onset developmental disorder characterized by deficits in reciprocal social interaction, verbal and non-verbal communication, and dependence on routines and rituals. It belongs to a spectrum of disorders (autism spectrum disorders, ASDs) which share core symptoms but show considerable variation in severity. The whole spectrum affects 0.6-0.7% of children worldwide, inducing a substantial public health burden and causing suffering to the affected families. Despite having a very high heritability, ASDs have shown exceptional genetic heterogeneity, which has complicated the identification of risk variants and left the etiology largely unknown. However, recent studies suggest that rare, family-specific factors contribute significantly to the genetic basis of ASDs. In this study, we investigated the role of DISC1 (Disrupted-in-schizophrenia-1) in ASDs, and identified association with markers and haplotypes previously associated with psychiatric phenotypes. We identified four polymorphic micro-RNA target sites in the 3 UTR of DISC1, and showed that hsa-miR-559 regulates DISC1 expression in vitro in an allele-specific manner. We also analyzed an extended autism pedigree with genealogical roots in Central Finland reaching back to the 17th century. To take advantage of the beneficial characteristics of population isolates to gene mapping and reduced genetic heterogeneity observed in distantly related individuals, we performed a microsatellite-based genome-wide screen for linkage and linkage disequilibrium in this pedigree. We identified a putative autism susceptibility locus on chromosome 19p13.3 and obtained further support for previously reported loci at 1q23 and 15q11-q13. To follow-up these findings, we extended our study sample from the same sub-isolate and initiated a genome-wide analysis of homozygosity and allelic sharing using high-density SNP markers. We identified a small number of haplotypes shared by different subsets of the genealogically connected cases, along with convergent biological pathways from SNP and gene expression data, which highlighted axon guidance molecules in the pathogenesis of ASDs. In conclusion, the results obtained in this thesis show that multiple distinct genetic variants are responsible for the ASD phenotype even within single pedigrees from an isolated population. We suggest that targeted resequencing of the shared haplotypes, linkage regions, and other susceptibility loci is essential to identify the causal variants. We also report a possible micro-RNA mediated regulatory mechanism, which might partially explain the wide-range neurobiological effects of the DISC1 gene.

Data integration from two microarray platforms identifies bi-allelic genetic inactivation of RIC8A in a breast cancer cell line

Background: Using array comparative genomic hybridization (aCGH), a large number of deleted genomic regions have been identified in human cancers. However, subsequent efforts to identify target genes selected for inactivation in these regions have often been challenging. Methods: We integrated here genome-wide copy number data with gene expression data and non-sense mediated mRNA decay rates in breast cancer cell lines to prioritize gene candidates that are likely to be tumour suppressor genes inactivated by bi-allelic genetic events. The candidates were sequenced to identify potential mutations. Results: This integrated genomic approach led to the identification of RIC8A at 11p15 as a putative candidate target gene for the genomic deletion in the ZR-75-1 breast cancer cell line. We identified a truncating mutation in this cell line, leading to loss of expression and rapid decay of the transcript. We screened 127 breast cancers for RIC8A mutations, but did not find any pathogenic mutations. No promoter hypermethylation in these tumours was detected either. However, analysis of gene expression data from breast tumours identified a small group of aggressive tumours that displayed low levels of RIC8A transcripts. qRT-PCR analysis of 38 breast tumours showed a strong association between low RIC8A expression and the presence of TP53 mutations (P = 0.006). Conclusion: We demonstrate a data integration strategy leading to the identification of RIC8A as a gene undergoing a classical double-hit genetic inactivation in a breast cancer cell line, as well as in vivo evidence of loss of RIC8A expression in a subgroup of aggressive TP53 mutant breast cancers.