Computational Methods for Reconstruction and Analysis of Genome-Scale Metabolic Networks

Metabolism is the cellular subsystem responsible for generation of energy from nutrients and production of building blocks for larger macromolecules. Computational and statistical modeling of metabolism is vital to many disciplines including bioengineering, the study of diseases, drug target identification, and understanding the evolution of metabolism. In this thesis, we propose efficient computational methods for metabolic modeling. The techniques presented are targeted particularly at the analysis of large metabolic models encompassing the whole metabolism of one or several organisms. We concentrate on three major themes of metabolic modeling: metabolic pathway analysis, metabolic reconstruction and the study of evolution of metabolism. In the first part of this thesis, we study metabolic pathway analysis. We propose a novel modeling framework called gapless modeling to study biochemically viable metabolic networks and pathways. In addition, we investigate the utilization of atom-level information on metabolism to improve the quality of pathway analyses. We describe efficient algorithms for discovering both gapless and atom-level metabolic pathways, and conduct experiments with large-scale metabolic networks. The presented gapless approach offers a compromise in terms of complexity and feasibility between the previous graph-theoretic and stoichiometric approaches to metabolic modeling. Gapless pathway analysis shows that microbial metabolic networks are not as robust to random damage as suggested by previous studies. Furthermore the amino acid biosynthesis pathways of the fungal species Trichoderma reesei discovered from atom-level data are shown to closely correspond to those of Saccharomyces cerevisiae. In the second part, we propose computational methods for metabolic reconstruction in the gapless modeling framework. We study the task of reconstructing a metabolic network that does not suffer from connectivity problems. Such problems often limit the usability of reconstructed models, and typically require a significant amount of manual postprocessing. We formulate gapless metabolic reconstruction as an optimization problem and propose an efficient divide-and-conquer strategy to solve it with real-world instances. We also describe computational techniques for solving problems stemming from ambiguities in metabolite naming. These techniques have been implemented in a web-based sofware ReMatch intended for reconstruction of models for 13C metabolic flux analysis. In the third part, we extend our scope from single to multiple metabolic networks and propose an algorithm for inferring gapless metabolic networks of ancestral species from phylogenetic data. Experimenting with 16 fungal species, we show that the method is able to generate results that are easily interpretable and that provide hypotheses about the evolution of metabolism.

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.

The efficacy and potential risks of controlling sprouting in Finnish birches (Betula spp.) with the fungal decomposer Chondrostereum purpureum

Sprouting of fast-growing broad-leaved trees causes problems in young coniferous stands, under power transmission lines and along roads and railways. Public opinion and the Finnish Forest Certification System oppose the use of chemical herbicides to control sprouting, which means that most areas with problems rely on mechanical cutting. However, cutting is a poor control method for many broad-leaved species because the removal of leaders can stimulate the sprouting of side branches and cut stumps quickly re-sprout. In order to be effective, cutting must be carried out frequently but each cut increases the costs, making this control method increasingly difficult and expensive once begun. As such, alternative methods for sprout control that are both effective and environmentally sound represent a continuing challenge to managers and research biologists. Using biological control agents to prevent sprouting has been given serious consideration recently. Dutch and Canadian researchers have demonstrated the potential of the white-rot fungus Chondrostereum purpureum (Pers. ex Fr.) Pouzar as a control agent of stump sprouting in many hardwoods. These findings have focused the attention of the Finnish forestry community on the utilization of C. purpureum for biocontrol purposes. Primarily, this study sought determines the efficacy of native C. purpureum as an inhibitor of birch stump sprouting in Finland and to clarify its mode of action. Additionally, genotypic variation in Finnish C. purpureum was examined and the environmental risks posed by a biocontrol program using this fungus were assessed. Experimental results of the study demonstrated that C. purpureum clearly affects the sprouting of birch: both the frequency of living stumps and the number of living sprouts per stump were effectively reduced by the treatment. However, the treatment had no effect on the maximum height of new sprouts. There were clear differences among fungal isolates in preventing sprouting and those that possessed high oxidative activities as measured in the laboratory inhibited sprouting most efficiently in the field. The most effective treatment time during the growing season was in early and mid summer (May July). Genetic diversity in Nordic and Baltic populations of C. purpureum was found to be high at the regional scale but locally homogeneous. This natural distribution of diversity means that using local genotypes in biocontrol programs would effectively prevent the introduction of novel genes or genotypes. While a biocontrol program using local strains of C. purpureum would be environmentally neutral, pruned birches that are close to the treatment site would have a high susceptibility to infect by the fungus during the early spring.

Monitoring of Fungal Growth and Degradation of Wood

Valko- ja ruskolahosienet tunnetaan luonnossa tehokkaimpina puun ja karikkeen lignoselluloosan lahottajina. Valkolahosienet pystyvät hajottamaan kaikkia puun osia: ligniiniä, selluloosaa ja hemiselluloosaa. Selektiivisesti ligniiniä hajottavat sienet lahottavat puusta suhteessa enemmän vaikeasti hajoavaa ligniiniä kuin selluloosaa tai hemiselluloosaa, jolloin jäljelle jää valkoista ja miltei puhdasta selluloosaa. Bioteknisissä sovelluksissa juuri selektiviiviset valkolahottajat ovat kiinnostavia. Niiden avulla voidaan puuhaketta esikäsitellä esimerkiksi paperinvalmistuksessa haitallisen ligniinin poistamiseksi. Ruskolahosienet ovat huomattavia puun, puutavaran ja puisten rakenteiden lahottajia, kuten tässä työssä käytetty Gloeophyllum trabeum (saunasieni ) ja Poria (Postia) placenta (istukkakääpä). Ruskolahosienet hajottavat puusta hemiselluloosan lisäksi selluloosaa, jolloin jää jäljelle ruskea ja jauhomaiseksi mureneva ligniini. Ruskolahosienet muovaavat ligniiniä jonkin verran. Kahden ruskolahosienen G. trabeumin ja P. placentan lisäksi tutkittiin valkolahosieniä, joista Ceriporiopsis subvermispora (karstakääpä) ja harvinainen Physisporinus rivulosus -sieni (talikääpä) hajottavat ligniiniä erittäin selektiivisesti. Phanerochaete chrysosporium on kaikkialla paljon tutkittu sieni, ja Phlebia radiata valkolahosientä (rusorypykkä) on tutkittu paljon mikrobiologian osastolla. Lisäksi tutkittiin Phlebia tremellosa -sienten (hytyrypykkä) ligninolyyttisten entsyymien tuottoa ja 14C-leimatun synteettisen ligniinin (DHP) hajotusta. P. radiata ja P. tremellosa -sienten on todettu aiemmin hajottavan ligniiniä selektiivisesti. Työssä selvitettiin miten sienten kasvua voi mitata, miten vertailukelpoisia eri mittaamismenetelmillä saadut tulokset ovat ja ilmenevätkö sienten aktiivisimmat kasvuvaiheet samaan aikaan eri menetelmillä mitattuna. Tärkeimmät tulokset olivat seuraavat havainnot: (i) P. radiata ja P. tremellosa -sienikannat tuottivat ligniini- ja mangaaniperoksidaasientsyymejä (LiP ja MnP) sekä lakkaasia, ja sienistä puhdistettiin 2-3 LiP- ja P. radiatasta yksi MnP-entsyymi; (ii) P. tremellosa -sienet hajottivat leimattua synteettistä ligniiniä (DHP) yhtä hyvin kuin paljon tutkitut P. chrysosporium ja P. radiata -sienet; (iii) puu, sienen luonnollinen kasvualusta, lisäsi valkolaho- ja ruskolahosienten demetoksylaatiota [O14CH3]-leimatusta ligniinin malliyhdisteestä 14CO2:ksi ilman puuta olleeseen alustaan verrattuna; (iv) demetoksylaatio (14CO2:n tuotto) oli normaalissa ilma-atmosfäärissä useimmiten parempi happeen verrattuna; (v) hapessa paras 14CO2:n tuotto saatiin puupalakasvatuksissa, joihin oli lisätty ravinnetyppeä tai typen lisäksi glukoosia sekä valkolaho- että ruskolahosienillä; (vi) ilmassa 14CO2:n tuotto oli puulla voimakkainta valkolahosienillä ilman lisäravinteita, kun taas G. trabeum -sienellä se oli yhtä hyvä eri alustoissa; (vii) biomassan muodostuminen rihmastojen ergosterolipitoisuuksista mitattuna oli ruskolahosienillä parempi kuin valkolahosienillä; (viii) ja biomassojen huippupitoisuudet olivat 6:lla sienellä eri suuruisia ja niiden maksimimäärien ajankohdat vaihtelivat viiden viikon kasvatusten kuluessa. Mikrobiologian osastolla Viikissä eristetty ja paljon tutkittu P. radiata -valkolahosieni oli mukana kaikissa tehdyissä kokeissa. Sienen LiP-aktiivisuus ja 14CO2:n tuotto 14C-rengas-leimatusta synteettisestä ligniinistä (DHP) korreloivat erittäin hyvin. Biomassan muodostuminen ergosterolilla määritettynä tuki hyvin entsyymiaktiivisuusmittauksilla ja isotooppikasvatuksilla saatuja tuloksia.

Characterization of small GTPase Cdc42 from the ectomycorrhizal fungus Suillus bovinus and Agrobacterium tumefaciens-mediated transformation of fungi

Ectomycorrhizal formation between the host tree, Pinus sylvestris and fungal symbiont, Suillus bovinus was investigated at the molecular level by isolating genes regulating the organization of the actin cytoskeleton in the fungal partner S. bovinus. An Agrobacterium tumefaciens mediated transformation (ATMT) system was developed for the ectomycorrhizal fungi in order to assign specific functions to the cloned molecules. The developed ATMT system was also used to transform a plant pathogenic fungus, Helminthosporium turcicum, to hygromycin B resistance. Small GTPases Cdc42 and Rac1, the regulators of actin cytoskeleton in eukaryotes were isolated from S. bovinus. Sbcdc42 and Sbrac1, are both expressed in vegetative and in the symbiotic hyphae of S. bovinus . Using IIF microscopy, Cdc42 and actin were co-localized at the tips of vegetative hyphae and were visualized in association with the plasma membrane in swollen cells typical to the symbiotic hyphae. These results suggest that the small GTPases Cdc42 may play a significant role in the polarized growth of S. bovinus hyphae and regulate fungal morphogenesis during ectomycorrhiza formation through reorganization of the actin cytoskeleton. The functional equality of Cdc42 was tested in yeast complementation experiments using a Saccharomyces cerevisiae temperature sensitive mutant, cdc42-1ts. The genomic clone of CDC42 was isolated from S. bovinus genomic DNA via specific primers for Cdc42. The analogous S. cerevisiae cdc42 mutations, dominant active G12V and dominant negative D118A, were generated in the Sbcdc42 gene by in-vitro mutagenesis. The ectomycorrhizal fungi, S. bovinus, P. involutus and H. cylindroporum were transformed using ATMT and phleomycin as a selectable marker. PCR screeing suggested that the T-DNA was inserted in all the three fungal genomes but the fate of integration could not be proved by Southern blot analysis. An alternative Agrobacterium strain, AGL-1 and selection marker, hygromycin was used to transform our model fungus S. bovinus. PCR and Southern analysis suggested an improved efficiency of transformation. All the transformed fungal colonies selected for hygromycin gave positives in PCR and the Southerns showed multiple or single copy T-DNA integrations into the S. bovinus genome. Using the same Agrobacterium strain and the selectable marker, a maize pathogen, H. turcicum was also subjected to ATMT. The H. turcicum transformation data suggested the single copy T-DNA integrations into the genome of the screened transformants that further confirms wider applicability of the ATMT. The plasmids carrying the wild-type (pHGCDC42) and the mutated Sbcdc42 alleles (pHGGV; pHGDA) under Agaricus bisporus gpd promoter were constructed in an A. tumefaciens vector. ATMT was used to transform S. bovinus with the plasmids carrying the wild-type and mutated Sbcdc42 alleles. The isolation of Sbcdc42 and Sbrac1 genes and some other functionally related genes from ectomycorrhizal fungus, S. bovinus will form the basis of future work to resolve the signalling pathway leading to ectomycorrhizal symbiosis. The development of ATMT system will be a valuable tool in analysing the exact function of signalling pathway components in ectomycorrhizal symbiosis or in plant pathogenic interactions. The transformation frequency and broad applicability along with the simplicity of T-DNA integration make Agrobacterium a valuable, new and a powerfull tool for targeted and insertional mutagenesis in these plant associated fungi. The developed ATMT systems should therefore make it possible to generate large number of transformants with tagged genes which could then be screened for their specific roles in symbiosis and pathogenecity, respectively.

Ecological impacts of Phlebiopsis gigantea biocontrol treatment against Heterobasidion spp. as revealed by fungal community profiling and population analyses

Wood decay fungi belonging to the species complex Heterobasidion annosum sensu lato are among the most common and economically important species causing root rot and stem decay in conifers of the northern temperate regions. New infections by these pathogens can be suppressed by tree stump treatments using chemical or biological control agents. In Finland, the corticiaceous fungus Phlebiopsis gigantea has been formulated into a commercial biocontrol agent called Rotstop (Verdera Ltd.). This thesis addresses the ecological impacts of Rotstop biocontrol treatment on the mycoflora of conifer stumps. Locally, fungal communities within Rotstop-treated and untreated stumps were analyzed using a novel method based on DGGE profiling of small subunit ribosomal DNA fragments amplified directly from wood samples. Population analyses for P. gigantea and H. annosum s.l. were conducted to evaluate possible risks associated with local and/or global distribution of the Rotstop strain. Based on molecular community profiling by DGGE, we detected a few individual wood-inhabiting fungal species (OTUs) that seemed to have suffered or benefited from the Rotstop biocontrol treatment. The DGGE analyses also revealed fungal diversity not retrieved by cultivation and some fungal sequence types untypical for decomposing conifer wood. However, statistical analysis of DGGE community profiles obtained from Rotstop-treated and untreated conifer stumps revealed that the Rotstop treatment had not caused a statistically significant reduction in the species diversity of wood-inhabiting fungi within our experimental forest plots. Locally, ISSR genotyping of cultured P. gigantea strains showed that the Rotstop biocontrol strain was capable of surviving up to six years within treated Norway spruce stumps, while in Scots pine stumps it was sooner replaced by successor fungal species. In addition, the spread of resident P. gigantea strains into Rotstop-treated forest stands seemed effective in preventing the formation of genetically monomorphic populations in the short run. On a global scale, we detected a considerable level of genetic differentiation between the interfertile European and North American populations of P. gigantea. These results strongly suggest that local biocontrol strains should be used in order to prevent global spread of P. gigantea and hybrid formation between geographically isolated populations. The population analysis for H. annosum s.l. revealed a collection of Chinese fungal strains that showed a high degree of laboratory fertility with three different allopatric H. annosum s.l. taxa. However, based on the molecular markers, the Chinese strains could be clearly affiliated with the H. parviporum taxonomical cluster, which thus appears to have a continuous distribution range from Europe through southern Siberia to northern China. Keywords: Rotstop, wood decay, DGGE, ISSR fingerprinting, ribosomal DNA

The genome packaging machinery of dsDNA bacteriophage PRD1

Viral genomes are encapsidated within protective protein shells. This encapsidation can be achieved either by a co-condensation reaction of the nucleic acid and coat proteins, or by first forming empty viral particles which are subsequently packaged with nucleic acid, the latter mechanism being typical for many dsDNA bacteriophages. Bacteriophage PRD1 is an icosahedral, non-tailed dsDNA virus that has an internal lipid membrane, the hallmark of the Tectiviridae family. Although PRD1 has been known to assemble empty particles into which the genome is subsequently packaged, the mechanism for this has been unknown, and there has been no evidence for a separate packaging vertex, similar to the portal structures used for packaging in the tailed bacteriophages and herpesviruses. In this study, a unique DNA packaging vertex was identified for PRD1, containing the packaging ATPase P9, packaging factor P6 and two small membrane proteins, P20 and P22, extending the packaging vertex to the internal membrane. Lack of small membrane protein P20 was shown to totally abolish packaging, making it an essential part of the PRD1 packaging mechanism. The minor capsid proteins P6 was shown to be an important packaging factor, its absence leading to greatly reduced packaging efficiency. An in vitro DNA packaging mechanism consisting of recombinant packaging ATPase P9, empty procapsids and mutant PRD1 DNA with a LacZ-insert was developed for the analysis of PRD1 packaging, the first such system ever for a virus containing an internal membrane. A new tectiviral sequence, a linear plasmid called pBClin15, was identified in Bacillus cereus, providing material for sequence analysis of the tectiviruses. Analysis of PRD1 P9 and other putative tectiviral ATPase sequences revealed several conserved sequence motifs, among them a new tectiviral packaging ATPase motif. Mutagenesis studies on PRD1 P9 were used to confirm the significance of the motifs. P9-type putative ATPase sequences carrying a similar sequence motif were identified in several other membrane containing dsDNA viruses of bacterial, archaeal and eukaryotic hosts, suggesting that these viruses may have similar packaging mechanisms. Interestingly, almost the same set of viruses that were found to have similar putative packaging ATPases had earlier been found to share similar coat protein folds and capsid structures, and a common origin for these viruses had been suggested. The finding in this study of similar packaging proteins further supports the idea that these viruses are descendants of a common ancestor.

Identification of the fungal catabolic D-galacturonate pathway

Pectin is a natural polymer consisting mainly of D-galacturonic acid monomers. Microorganisms living on decaying plant material can use D-galacturonic acid for growth. Although bacterial pathways for D-galacturonate catabolism had been described previously, no eukaryotic pathway for D-galacturonate catabolism was known at the beginning of this work. The aim of this work was to identify such a pathway. In this thesis the pathway for D-galacturonate catabolism was identified in the filamentous fungus Trichoderma reesei. The pathway consisted of four enzymes: NADPH-dependent D-galacturonate reductase (GAR1), L-galactonate dehydratase (LGD1), L-threo-3-deoxy-hexulosonate aldolase (LGA1) and NADPH-dependent glyceraldehyde reductase (GLD1). In this pathway D-galacturonate was converted to pyruvate and glycerol via L-galactonate, L-threo-3-deoxy-hexulosonate and L-glyceraldehyde. The enzyme activities of GAR1, LGD1 and LGA1 were present in crude mycelial extract only when T. reesei was grown on D-galacturonate. The activity of GLD1 was equally present on all the tested carbon sources. The corresponding genes were identified either by purifying and sequencing the enzyme or by expressing genes with homology to other similar enzymes in a heterologous host and testing the activities. The new genes that were identified were expressed in Saccharomyces cerevisiae and resulted in active enzymes. The GAR1, LGA1 and GLD1 were also produced in S. cerevisiae as active enzymes with a polyhistidine-tag, and purified and characterised. GAR1 and LGA1 catalysed reversible reactions, whereas only the forward reactions were observed for LGD1 and GLD1. When gar1, lgd1 or lga1 was deleted in T. reesei the deletion strain was unable to grow with D-galacturonate as the only carbon source, demonstrating that all the corresponding enzymes were essential for D-galacturonate catabolism and that no alternative D-galacturonate pathway exists in T. reesei. A challenge for biotechnology is to convert cheap raw materials to useful and more valuable products. Filamentous fungi are especially useful for the conversion of pectin, since they are efficient producers of pectinases. Identification of the fungal D-galacturonate pathway is of fundamental importance for the utilisation of pectin and its conversion to useful products.

Two ex situ fungal technologies to treat contaminated soil

Wood-degrading fungi are able to degrade a large range of recalcitrant pollutants which resemble the lignin biopolymer. This ability is attributed to the production of lignin-modifying enzymes, which are extracellular and non-specific. Despite the potential of fungi in bioremediation, there is still an understanding gap in terms of the technology. In this thesis, the feasibility of two ex situ fungal bioremediation methods to treat contaminated soil was evaluated. Treatment of polycyclic aromatic hydrocarbons (PAHs)-contaminated marsh soil was studied in a stirred slurry-phase reactor. Due to the salt content in marsh soil, fungi were screened for their halotolerance, and the white-rot fungi Lentinus tigrinus, Irpex lacteus and Bjerkandera adusta were selected for further studies. These fungi degraded 40 - 60% of a PAH mixture (phenanthrene, fluoranthene, pyrene and chrysene) in a slurry-phase reactor (100 ml) during 30 days of incubation. Thereafter, B. adusta was selected to scale-up and optimize the process in a 5 L reactor. Maximum degradation of dibenzothiophene (93%), fluoranthene (82%), pyrene (81%) and chrysene (83%) was achieved with the free mycelium inoculum of the highest initial biomass (2.2 g/l). In autoclaved soil, MnP was the most important enzyme involved in PAH degradation. In non-sterile soil, endogenous soil microbes together with B. adusta also degraded the PAHs extensively, suggesting a synergic action between soil microbes and the fungus. A fungal solid-phase cultivation method to pretreat contaminated sawmill soil with high organic matter content was developed to enhance the effectiveness of the subsequent soil combustion. In a preliminary screening of 146 fungal strains, 28 out of 52 fungi, which extensively colonized non-sterile contaminated soil, were litter-decomposing fungi. The 18 strains further selected were characterized by their production of lignin-modifying and hydrolytic enzymes, of which MnP and endo-1,4-β-glucanase were the main enzymes during cultivation on Scots pine (Pinus sylvestris) bark. Of the six fungi selected for further tests, Gymnopilus luteofolius, Phanerochaete velutina, and Stropharia rugosoannulata were the most active soil organic matter degraders. The results showed that a six-month pretreatment of sawmill soil would result in a 3.5 - 9.5% loss of organic matter, depending on the fungus applied. The pretreatment process was scaled-up for a 0.56 m3 reactor, in which perforated plastic tubes filled with S. rugosoannulata growing on pine bark were introduced into the soil. The fungal pretreatment resulted in a soil mass loss of 30.5 kg, which represents 10% of the original soil mass (308 kg). Despite the fact that Scots pine bark contains several antimicrobial compounds, it was a suitable substrate for fungal growth and promoter of the production of oxidative enzymes, as well as an excellent and cheap natural carrier of fungal mycelium. This thesis successfully developed two novel fungal ex situ bioremediation technologies and introduce new insights for their further full-scale application. Ex situ slurry-phase fungal reactors might be applied in cases when the soil has a high water content or when the contaminant bioavailability is low; for example, in wastewater treatment plants to remove pharmaceutical residues. Fungal solid-phase bioremediation is a promising remediation technology to ex situ or in situ treat contaminated soil.

Optimization of sample preparation for next-generation sequencing with Illumina Genome Analyzer II

The growing interest for sequencing with higher throughput in the last decade has led to the development of new sequencing applications. This thesis concentrates on optimizing DNA library preparation for Illumina Genome Analyzer II sequencer. The library preparation steps that were optimized include fragmentation, PCR purification and quantification. DNA fragmentation was performed with focused sonication in different concentrations and durations. Two column based PCR purification method, gel matrix method and magnetic bead based method were compared. Quantitative PCR and gel electrophoresis in a chip were compared for DNA quantification. The magnetic bead purification was found to be the most efficient and flexible purification method. The fragmentation protocol was changed to produce longer fragments to be compatible with longer sequencing reads. Quantitative PCR correlates better with the cluster number and should thus be considered to be the default quantification method for sequencing. As a result of this study more data have been acquired from sequencing with lower costs and troubleshooting has become easier as qualification steps have been added to the protocol. New sequencing instruments and applications will create a demand for further optimizations in future.

Discovery of oxidative enzymes for food engineering : Tyrosinase and sulfhydryl oxidase

Enzymes offer many advantages in industrial processes, such as high specificity, mild treatment conditions and low energy requirements. Therefore, the industry has exploited them in many sectors including food processing. Enzymes can modify food properties by acting on small molecules or on polymers such as carbohydrates or proteins. Crosslinking enzymes such as tyrosinases and sulfhydryl oxidases catalyse the formation of novel covalent bonds between specific residues in proteins and/or peptides, thus forming or modifying the protein network of food. In this study, novel secreted fungal proteins with sequence features typical of tyrosinases and sulfhydryl oxidases were iden-tified through a genome mining study. Representatives of both of these enzyme families were selected for heterologous produc-tion in the filamentous fungus Trichoderma reesei and biochemical characterisation. Firstly, a novel family of putative tyrosinases carrying a shorter sequence than the previously characterised tyrosinases was discovered. These proteins lacked the whole linker and C-terminal domain that possibly play a role in cofactor incorporation, folding or protein activity. One of these proteins, AoCO4 from Aspergillus oryzae, was produced in T. reesei with a production level of about 1.5 g/l. The enzyme AoCO4 was correctly folded and bound the copper cofactors with a type-3 copper centre. However, the enzyme had only a low level of activity with the phenolic substrates tested. Highest activity was obtained with 4-tert-butylcatechol. Since tyrosine was not a substrate for AoCO4, the enzyme was classified as catechol oxidase. Secondly, the genome analysis for secreted proteins with sequence features typical of flavin-dependent sulfhydryl oxidases pinpointed two previously uncharacterised proteins AoSOX1 and AoSOX2 from A. oryzae. These two novel sulfhydryl oxidases were produced in T. reesei with production levels of 70 and 180 mg/l, respectively, in shake flask cultivations. AoSOX1 and AoSOX2 were FAD-dependent enzymes with a dimeric tertiary structure and they both showed activity on small sulfhydryl compounds such as glutathione and dithiothreitol, and were drastically inhibited by zinc sulphate. AoSOX2 showed good stabil-ity to thermal and chemical denaturation, being superior to AoSOX1 in this respect. Thirdly, the suitability of AoSOX1 as a possible baking improver was elucidated. The effect of AoSOX1, alone and in combi-nation with the widely used improver ascorbic acid was tested on yeasted wheat dough, both fresh and frozen, and on fresh water-flour dough. In all cases, AoSOX1 had no effect on the fermentation properties of fresh yeasted dough. AoSOX1 nega-tively affected the fermentation properties of frozen doughs and accelerated the damaging effects of the frozen storage, i.e. giving a softer dough with poorer gas retention abilities than the control. In combination with ascorbic acid, AoSOX1 gave harder doughs. In accordance, rheological studies in yeast-free dough showed that the presence of only AoSOX1 resulted in weaker and more extensible dough whereas a dough with opposite properties was obtained if ascorbic acid was also used. Doughs containing ascorbic acid and increasing amounts of AoSOX1 were harder in a dose-dependent manner. Sulfhydryl oxidase AoSOX1 had an enhancing effect on the dough hardening mechanism of ascorbic acid. This was ascribed mainly to the produc-tion of hydrogen peroxide in the SOX reaction which is able to convert the ascorbic acid to the actual improver dehydroascorbic acid. In addition, AoSOX1 could possibly oxidise the free glutathione in the dough and thus prevent the loss of dough strength caused by the spontaneous reduction of the disulfide bonds constituting the dough protein network. Sulfhydryl oxidase AoSOX1 is therefore able to enhance the action of ascorbic acid in wheat dough and could potentially be applied in wheat dough baking.

Response of fungal and actinobacterial communities to water-level drawdown in boreal peatland sites

"We used PCR-DGGE fingerprinting and direct sequencing to analyse the response of fungal and actinobacterial communities to changing hydrological conditions at 3 different sites in a boreal peatland complex in Finland. The experimental design involved a short-term (3 years; STD) and a long-term (43 years; LTD) water-level drawdown. Correspondence analyses of DGGE bands revealed differences in the communities between natural sites representing the nutrient-rich mesotrophic fen, the nutrient-poorer oligotrophic fen, and the nutrient-poor ombrotrophic bog. Still, most fungi and actinobacteria found in the pristine peatland seemed robust to the environmental variables. Both fungal and actinobacterial diversity was higher in the fens than in the bog. Fungal diversity increased significantly after STD whereas actinobacterial diversity did not respond to hydrology. Both fungal and actinobacterial communities became more similar between peatland types after LTD, which was not apparent after STD. Most sequences clustered equally between the two main fungal phyla Ascomycota and Basidiomycota. Sequencing revealed that basidiomycetes may respond more (either positively or negatively) to hydrological changes than ascomycetes. Overall, our results suggest that fungal responses to water-level drawdown depend on peatland type. Actinobacteria seem to be less sensitive to hydrological changes, although the response of some may similarly depend on peatland type. (C) 2009 Elsevier Ltd. All rights reserved."

Microbial identification by detection of ligation probes on DNA microarray

Microbes in natural and artificial environments as well as in the human body are a key part of the functional properties of these complex systems. The presence or absence of certain microbial taxa is a correlate of functional status like risk of disease or course of metabolic processes of a microbial community. As microbes are highly diverse and mostly notcultivable, molecular markers like gene sequences are a potential basis for detection and identification of key types. The goal of this thesis was to study molecular methods for identification of microbial DNA in order to develop a tool for analysis of environmental and clinical DNA samples. Particular emphasis was placed on specificity of detection which is a major challenge when analyzing complex microbial communities. The approach taken in this study was the application and optimization of enzymatic ligation of DNA probes coupled with microarray read-out for high-throughput microbial profiling. The results show that fungal phylotypes and human papillomavirus genotypes could be accurately identified from pools of PCR amplicons generated from purified sample DNA. Approximately 1 ng/μl of sample DNA was needed for representative PCR amplification as measured by comparisons between clone sequencing and microarray. A minimum of 0,25 amol/μl of PCR amplicons was detectable from amongst 5 ng/μl of background DNA, suggesting that the detection limit of the test comprising of ligation reaction followed by microarray read-out was approximately 0,04%. Detection from sample DNA directly was shown to be feasible with probes forming a circular molecule upon ligation followed by PCR amplification of the probe. In this approach, the minimum detectable relative amount of target genome was found to be 1% of all genomes in the sample as estimated from 454 deep sequencing results. Signal-to-noise of contact printed microarrays could be improved by using an internal microarray hybridization control oligonucleotide probe together with a computational algorithm. The algorithm was based on identification of a bias in the microarray data and correction of the bias as shown by simulated and real data. The results further suggest semiquantitative detection to be possible by ligation detection, allowing estimation of target abundance in a sample. However, in practise, comprehensive sequence information of full length rRNA genes is needed to support probe design with complex samples. This study shows that DNA microarray has the potential for an accurate microbial diagnostic platform to take advantage of increasing sequence data and to replace traditional, less efficient methods that still dominate routine testing in laboratories. The data suggests that ligation reaction based microarray assay can be optimized to a degree that allows good signal-tonoise and semiquantitative detection.