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.

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.

Xanthine oxidoreductase in essential hypertension and metabolic syndrome : Experimental studies on rodent models

Hypertension, obesity, dyslipidemia and dysglycemia constitute metabolic syndrome, a major public health concern, which is associated with cardiovascular mortality. High dietary salt (NaCl) is the most important dietary risk factor for elevated blood pressure. The kidney has a major role in salt-sensitive hypertension and is vulnerable to harmful effects of increased blood pressure. Elevated serum urate is a common finding in these disorders. While dysregulation of urate excretion is associated with cardiovascular diseases, present studies aimed to clarify the role of xanthine oxidoreductase (XOR), i.e. xanthine dehydrogenase (XDH) and its post-translational isoform xanthine oxidase (XO), in cardiovascular diseases. XOR yields urate from hypoxanthine and xanthine. Low oxygen levels upregulate XOR in addition to other factors. In present studies higher renal XOR activity was found in hypertension-prone rats than in the controls. Furthermore, NaCl intake increased renal XOR dose-dependently. To clarify whether XOR has any causal role in hypertension, rats were kept on NaCl diets for different periods of time, with or without a XOR inhibitor, allopurinol. While allopurinol did not alleviate hypertension, it prevented left ventricular and renal hypertrophy. Nitric oxide synthases (NOS) produce nitric oxide (NO), which mediates vasodilatation. A paucity of NO, produced by NOS inhibition, aggravated hypertension and induced renal XOR, whereas NO generating drug, alleviated salt-induced hypertension without changes in renal XOR. Zucker fa/fa rat is an animal model of metabolic syndrome. These rats developed substantial obesity and modest hypertension and showed increased hepatic and renal XOR activities. XOR was modified by diet and antihypertensive treatment. Cyclosporine (CsA) is a fungal peptide and one of the first-line immunosuppressive drugs used in the management of organ transplantation. Nephrotoxicity ensue high doses resulting in hypertension and limit CsA use. CsA increased renal XO substantially in salt-sensitive rats on a high NaCl diet, indicating a possible role for this reactive oxygen species generating isoform in CsA nephrotoxicity. Renal hypoxia, common to these rodent models of hypertension and obesity, is one of the plausible XOR inducing factors. Although XOR inhibition did not prevent hypertension, present experimental data indicate that XOR plays a role in the pathology of salt-induced cardiac and renal hypertrophy.

Serum uric acid and metabolic risk factors in three ethnic groups: Asian Indians and Creoles in Mauritius and Chinese in Qingdao, China

In humans with a loss of uricase the final oxidation product of purine catabolism is uric acid (UA). The prevalence of hyperuricemia has been increasing around the world accompanied by a rapid increase in obesity and diabetes. Since hyperuricemia was first described as being associated with hyperglycemia and hypertension by Kylin in 1923, there has been a growing interest in the association between elevated UA and other metabolic abnormalities of hyperglycemia, abdominal obesity, dyslipidemia, and hypertension. The direction of causality between hyperuricemia and metabolic disorders, however, is unceartain. The association of UA with metabolic abnormalities still needs to be delineated in population samples. Our overall aims were to study the prevalence of hyperuricemia and the metabolic factors clustering with hyperuricemia, to explore the dynamical changes in blood UA levels with the deterioration in glucose metabolism and to estimate the predictive capability of UA in the development of diabetes. Four population-based surveys for diabetes and other non-communicable diseases were conducted in 1987, 1992, and 1998 in Mauritius, and in 2001-2002 in Qingdao, China. The Qingdao study comprised 1 288 Chinese men and 2 344 women between 20-74, and the Mauritius study consisted of 3 784 Mauritian Indian and Mauritian Creole men and 4 442 women between 25-74. In Mauritius, re-exams were made in 1992 and/or 1998 for 1 941 men (1 409 Indians and 532 Creoles) and 2 318 non pregnant women (1 645 Indians and 673 Creoles), free of diabetes, cardiovascular diseases, and gout at baseline examinations in 1987 or 1992, using the same study protocol. The questionnaire was designed to collect demographic details, physical examinations and standard 75g oral glucose tolerance tests were performed in all cohorts. Fasting blood UA and lipid profiles were also determined. The age-standardized prevalence in Chinese living in Qingdao was 25.3% for hyperuricemia (defined as fasting serum UA > 420 μmol/l in men and > 360 μmol/l in women) and 0.36% for gout in adults between 20-74. Hyperuricemia was more prevalent in men than in women. One standard deviation increase in UA concentration was associated with the clustering of metabolic risk factors for both men and women in three ethnic groups. Waist circumference, body mass index, and serum triglycerides appeared to be independently associated with hyperuricemia in both sexes and in all ethnic groups except in Chinese women, in whom triglycerides, high-density lipoprotein cholesterol, and total cholesterol were associated with hyperuricemia. Serum UA increased with increasing fasting plasma glucose levels up to a value of 7.0 mmol/l, but significantly decreased thereafter in mainland Chinese. An inverse relationship occurred between 2-h plasma glucose and serum UA when 2-h plasma glucose higher than 8.0 mmol/l. In the prospective study in Mauritius, 337 (17.4%) men and 379 (16.4%) women developed diabetes during the follow-up. Elevated UA levels at baseline increased 1.14-fold in risk of incident diabetes in Indian men and 1.37-fold in Creole men, but no significant risk was observed in women. In conclusion, the prevalence of hyperuricemia was high in Chinese in Qingdao, blood UA was associated with the clustering of metabolic risk factors in Mauritian Indian, Mauritian Creole, and Chinese living in Qingdao, and a high baseline UA level independently predicted the development of diabetes in Mauritian men. The clinical use of UA as a marker of hyperglycemia and other metabolic disorders needs to be further studied. Keywords: Uric acid, Hyperuricemia, Risk factors, Type 2 Diabetes, Incidence, Mauritius, Chinese

Microbial quality of hemp (Cannabis sativa L.) and flax (Linum usitatissimum L.) from plants to thermal insulation

Flax and hemp have traditionally been used mainly for textiles, but recently interest has also been focused on non-textile applications. Microbial quality throughout the whole processing chain of bast fibres has not previously been studied. This study concentrates on the microbial quality and possible microbial risks in the production chain of hemp and flax fibres and fibrous thermal insulations. In order to be able to utilize hemp and flax fibres, the bast fibres must be separated from the rest of the plant. Non-cellulosic components can be removed with various pretreatment processes, which are associated with a certain risk of microbial contamination. In this study enzymatic retting and steam explosion (STEX) were examined as pretreatment processes. On the basis of the results obtained in this study, the microbial contents on stalks of both plants studied increased at the end of the growing season and during the winter. However, by processing and mechanical separation it is possible to produce fibres containing less moulds and bacteria than the whole stem. Enzymatic treatment encouraged the growth of moulds in fibres. Steam explosion reduced the amount of moulds in fibres. Dry thermal treatment used in this study did not markedly reduce the amount of microbes. In this project an emission measurement chamber was developed which was suitable for measurements of emissions from both mat type and loose fill type insulations, and capable of interdisciplinary sampling. In this study, the highest amounts of fungal emissions were in the range of 10^3 10^5 cfu/m^3 from the flax and hemp insulations at 90% RH of air. The fungal emissions from stone wool, glass wool and recycled paper insulations were below 10^2 cfu/m^3 even at 90% RH. Equally low values were obtained from bast fibrous materials in lower humidities (at 30% and 80% RH of air). After drying of moulded insulations at 30% RH, the amounts of emitted moulds were in all cases higher compared to the emissions at 90% RH before drying. The most common fungi in bast fibres were Penicillium and Rhizopus. The widest variety of different fungi was in the untreated hemp and linseed fibres and in the commercial loose-fill flax insulation. Penicillium, Rhizopus and Paecilomyces were the most tolerant to steam explosion. According to the literature, the most common fungi in building materials and indoor air are Penicillium, Aspergillus and Cladosporium, which were all found in some of the bast fibre materials in this study. As organic materials, hemp and flax fibres contain high levels of nutrients for microbial growth. The amount of microbes can be controlled and somewhat decreased by the processing methods presented.

Plant secondary compounds and soil microbial processes in carbon and nitrogen cycling in relation to tree species

The aim of this study was to explore soil microbial activities related to C and N cycling and the occurrence and concentrations of two important groups of plant secondary compounds, terpenes and phenolic compounds, under silver birch (Betula pendula Roth), Norway spruce (Picea abies (L.) Karst) and Scots pine (Pinus sylvestris L.) as well as to study the effects of volatile monoterpenes and tannins on soil microbial activities. The study site, located in Kivalo, northern Finland, included ca. 70-year-old adjacent stands dominated by silver birch, Norway spruce and Scots pine. Originally the soil was very probably similar in all three stands. All forest floor layers (litter (L), fermentation layer (F) and humified layer (H)) under birch and spruce showed higher rates of CO2 production, greater net mineralisation of nitrogen and higher amounts of carbon and nitrogen in microbial biomass than did the forest floor layers under pine. Concentrations of mono-, sesqui-, di- and triterpenes were higher under both conifers than under birch, while the concentration of total water-soluble phenolic compounds as well as the concentration of condensed tannins tended to be higher or at least as high under spruce as under birch or pine. In general, differences between tree species in soil microbial activities and in concentrations of secondary compounds were smaller in the H layer than in the upper layers. The rate of CO2 production and the amount of carbon in the microbial biomass correlated highly positively with the concentration of total water-soluble phenolic compounds and positively with the concentration of condensed tannins. Exposure of soil to volatile monoterpenes and tannins extracted and fractionated from spruce and pine needles affected carbon and nitrogen transformations in soil, but the effects were dependent on the compound and its molecular structure. Monoterpenes decreased net mineralisation of nitrogen and probably had a toxic effect on part of the microbial population in soil, while another part of the microbes seemed to be able to use monoterpenes as a carbon source. With tannins, low-molecular-weight compounds (also compounds other than tannins) increased soil CO2 production and nitrogen immobilisation by soil microbes while the higher-molecular-weight condensed tannins had inhibitory effects. In conclusion, plant secondary compounds may have a great potential in regulation of C and N transformations in forest soils, but the real magnitude of their significance in soil processes is impossible to estimate.

Responses of soil microbial communities to clonal variation of Norway spruce

In boreal forests, microorganisms have a pivotal role in nutrient and water supply of trees as well as in litter decomposition and nutrient cycling. This reinforces the link between above-ground and below-ground communities in the context of sustainable productivity of forest ecosystems. In northern boreal forests, the diversity of microbes associated with the trees is high compared to the number of distinct tree species. In this thesis, the aim was to study whether conspecific tree individuals harbour different soil microbes and whether the growth of the trees and the community structure of the associated microbes are connected. The study was performed in a clonal field trial of Norway spruce, which was established in a randomized block design in a clear-cut area. Since out-planting in 1994, the spruce clones showed two-fold growth differences. The fast-growing spruce clones were associated with a more diverse community of ectomycorrhizal fungi than the slow-growing spruce clones. These growth performance groups also differed with respect to other aspects of the associated soil microorganisms: the species composition of ectomycorrhizal fungi, in the amount of extraradical fungal mycelium, in the structure of bacterial community associated with the mycelium, and in the structure of microbial community in the organic layer. The communities of fungi colonizing needle litter of the spruce clones in the field did not differ and the loss of litter mass after two-years decomposition was equal. In vitro, needles of the slow-growing spruce clones were colonized by a more diverse community of endophytic fungi that were shown to be significant needle decomposers. This study showed a relationship between the growth of Norway spruce clones and the community structure of the associated soil microbes. Spatial heterogeneity in soil microbial community was connected with intraspecific variation of trees. The latter may therefore influence soil biodiversity in monospecific forests.

Microbial life and deposits in paper machine circuits

Väitöskirjani käsittele mikrobien ja erilaisten kemikaalien rooleja saostumien ja biofilmien muodostumisessa paperi- ja kartonkikoneilla. "Saostuma" tässä työssä tarkoittaa kiinteän aineen kertymää konepinnoille tai rajapinnoille konekierroissa, jotka on tarkoitettu massasulppujen, lietteiden, vesien tai ilman kuljetukseen. Saostumasta tulee "biofilmi" silloin kun sen oleellinen rakennekomponentti on mikrobisolut tai niiden tuotteet. Väitöstyöni työhypoteesina oli, että i. tietämys saostumien koostumuksesta, sekä ii. niiden rakenteesta, biologisista, fysikaalis-kemiallisista ja teknisistä ominaisuuksista ohjaavat tutkijaa löytämään ympäristöä säästäviä keinoja estää epätoivottujen saostumien muodostus tai purkaa jo muodostuneita saostumia. Selvittääkseni saostumien koostumista ja rakennetta käytin monia erilaisia analytiikan työkaluja, kuten elektronimikroskopiaa, konfokaali-laser mikroskopiaa (CLSM), energiadispersiivistä röntgenanalyysiä (EDX), pyrolyysi kaasukromatografiaa yhdistettynä massaspektrometriaan (Py-GCMS), joninvaihtokromatografiaa, kaasukromatografiaa ja mikrobiologisia analyysejä. Osallistuin aktiivisesti innovatiivisen, valon takaisinsirontaan perustuvan sensorin kehittämistyöhön, käytettäväksi biofilmin kasvun mittaukseen suoraan koneen vesikierroista ja säiliöistä. Työni osoitti, että monet paperinvalmistuksessa käytetyistä kemikaaleista reagoivat keskenään tuottaen orgaanisia tahmakerroksia konekiertojen teräspinnoille. Löysin myös kerrostumia, jotka valomikroskooppisessa tarkastelussa oli tulkittu mikrobeiksi, mutta jotka elektronimikroskopia paljasti alunasta syntyneiksi, alumiinihydroksidiksi joka saostui pH:ssa 6,8 kiertokuitua käyttävän koneen viiravesistä. Monet paperintekijät käyttävät vieläkin alunaa kiinnitysaineena vaikka prosessiolot ovat muuttuneet happamista neutraaleiksi. Sitä pidetään paperitekijän "aspiriinina", mutta väitöstutkimukseni osoitti sen riskit. Löysin myös orgaanisia saostumia, joiden alkuperä oli aineiden, kuten pihkan, saippuoituminen (kalsium saippuat) niin että muodostui tahmankasvua ylläpitävä alusta monilla paperi- ja kartonkikoneilla. Näin solumuodoiltaan Deinococcus geothermalista muistuttavia bakteereita kasvamassa lujasti teräskoepalojen pintaan kiinnittyneinä pesäkkeinä, kun koepaloja upotettiin paperikoneiden vesikiertoihin. Nämä deinokokkimaiset pesäkkeet voivat toimia jalustana, tarttumisalustana muiden mikrobien massoille, joka selittäisi miksi saostumat yleisesti sisältävät deinokokkeja pienenä, muttei koskaan pääasiallisena rakenneosana. Kun paperikoneiden käyttämien vesien (raakavedet, lämminvesi, biologisesti puhdistettu jätevesi) laatua tutkitaan, mittausmenetelmällä on suuri merkitys. Koepalan upotusmenetelmällä todettu biofilmikasvu ja viljelmenetelmällä mitattu bakteerisaastuneisuus korreloivat toisiinsa huonosti etenkin silloin kun likaantumisessa oli mukana rihmamaiseti kasvavia bakteereja. Huoli ympäristöstä on pakottanut paperi- ja kartonkikoneiden vesikiertojen sulkemiseen. Vesien kierrätys ja prosessivesien uudelleenkäyttö nostavat prosessilämpötilaa ja lisäävät koneella kiertävien kolloidisten ja liuenneiden aineiden määriä. Tutkin kiertovesien pitoisuuksia kolmessa eriasteisesti suljetussa tehtaassa, joiden päästöt olivat 0 m3, 0,5 m3 ja 4 m3 jätevettä tuotetonnia kohden, perustuen puhdistetun jäteveden uudelleen käyttöön. Nollapäästöisellä tehtaalla kiertovesiin kertyi paljon orgaanisesti sidottua hiiltä (> 10 g L-1), etenkin haihtuvina happoina (maito-, etikka-, propioni- ja voi-). Myös sulfaatteja, klorideja, natriumia ja kalsiumia kertyi paljon, > 1 g L-1 kutakin. Pääosa (>40%) kaikista bakteereista oli 16S rRNA geenisekvenssianalyysien tulosten perusteella sukua, joskin etäistä (< 96%) ainoastaan Enterococcus cecorum bakteerille. 4 m3 päästävältä tehtaalta löytyi lisäksi Bacillus thermoamylovorans ja Bacillus coagulans. Tehtaiden saostumat sisälsivät arkkeja suurina pitoisuuksina, ≥ 108 g-1, mutta tunnistukseen riittävää sekvenssisamanlaisuutta löytyi vain yhteen arkkisukuun, Methanothrix. Tutkimustulokset osoittivat että tehtaan vesikiertojen sulkeminen vähensi rajusti mikrobiston monimuotoisuutta, muttei estänyt liuenneen aineen ja kiintoaineen mineralisoitumista.

Algorithms for 13C metabolic flux analysis

The metabolism of an organism consists of a network of biochemical reactions that transform small molecules, or metabolites, into others in order to produce energy and building blocks for essential macromolecules. The goal of metabolic flux analysis is to uncover the rates, or the fluxes, of those biochemical reactions. In a steady state, the sum of the fluxes that produce an internal metabolite is equal to the sum of the fluxes that consume the same molecule. Thus the steady state imposes linear balance constraints to the fluxes. In general, the balance constraints imposed by the steady state are not sufficient to uncover all the fluxes of a metabolic network. The fluxes through cycles and alternative pathways between the same source and target metabolites remain unknown. More information about the fluxes can be obtained from isotopic labelling experiments, where a cell population is fed with labelled nutrients, such as glucose that contains 13C atoms. Labels are then transferred by biochemical reactions to other metabolites. The relative abundances of different labelling patterns in internal metabolites depend on the fluxes of pathways producing them. Thus, the relative abundances of different labelling patterns contain information about the fluxes that cannot be uncovered from the balance constraints derived from the steady state. The field of research that estimates the fluxes utilizing the measured constraints to the relative abundances of different labelling patterns induced by 13C labelled nutrients is called 13C metabolic flux analysis. There exist two approaches of 13C metabolic flux analysis. In the optimization approach, a non-linear optimization task, where candidate fluxes are iteratively generated until they fit to the measured abundances of different labelling patterns, is constructed. In the direct approach, linear balance constraints given by the steady state are augmented with linear constraints derived from the abundances of different labelling patterns of metabolites. Thus, mathematically involved non-linear optimization methods that can get stuck to the local optima can be avoided. On the other hand, the direct approach may require more measurement data than the optimization approach to obtain the same flux information. Furthermore, the optimization framework can easily be applied regardless of the labelling measurement technology and with all network topologies. In this thesis we present a formal computational framework for direct 13C metabolic flux analysis. The aim of our study is to construct as many linear constraints to the fluxes from the 13C labelling measurements using only computational methods that avoid non-linear techniques and are independent from the type of measurement data, the labelling of external nutrients and the topology of the metabolic network. The presented framework is the first representative of the direct approach for 13C metabolic flux analysis that is free from restricting assumptions made about these parameters.In our framework, measurement data is first propagated from the measured metabolites to other metabolites. The propagation is facilitated by the flow analysis of metabolite fragments in the network. Then new linear constraints to the fluxes are derived from the propagated data by applying the techniques of linear algebra.Based on the results of the fragment flow analysis, we also present an experiment planning method that selects sets of metabolites whose relative abundances of different labelling patterns are most useful for 13C metabolic flux analysis. Furthermore, we give computational tools to process raw 13C labelling data produced by tandem mass spectrometry to a form suitable for 13C metabolic flux analysis.

Below-ground processes in meadow soil under elevated ozone and carbon dioxide : - Greenhouse gas fluxes, N cycling and microbial communities

This thesis focuses on how elevated CO2 and/or O3 affect the below-ground processes in semi-natural vegetation, with an emphasis on greenhouse gases, N cycling and microbial communities. Meadow mesocosms mimicking lowland hay meadows in Jokioinen, SW Finland, were enclosed in open-top chambers and exposed to ambient and elevated levels of O3 (40-50 ppb) and/or CO2 (+100 ppm) for three consecutive growing season, while chamberless plots were used as chamber controls. Chemical and microbiological analyses as well as laboratory incubations of the mesocosm soils under different treatments were used to study the effects of O3 and/or CO2. Artificially constructed mesocosms were also compared with natural meadows with regards to GHG fluxes and soil characteristics. In addition to research conducted at the ecosystem level (i.e. the mesocosm study), soil microbial communities were also examined in a pot experiment with monocultures of individual species. By comparing mesocosms with similar natural plant assemblage, it was possible to demonstrate that artificial mesocosms simulated natural habitats, even though some differences were found in the CH4 oxidation rate, soil mineral N, and total C and N concentrations in the soil. After three growing seasons of fumigations, the fluxes of N2O, CH4, and CO2 were decreased in the NF+O3 treatment, and the soil NH4+-N and mineral N concentrations were lower in the NF+O3 treatment than in the NF control treatment. The mesocosm soil microbial communities were affected negatively by the NF+O3 treatment, as the total, bacterial, actinobacterial, and fungal PLFA biomasses as well as the fungal:bacterial biomass ratio decreased under elevated O3. In the pot survey, O3 decreased the total, bacterial, actinobacterial, and mycorrhizal PLFA biomasses in the bulk soil and affected the microbial community structure in the rhizosphere of L. pratensis, whereas the bulk soil and rhizosphere of the other monoculture, A. capillaris, remained unaffected by O3. Elevated CO2 caused only minor and insignificant changes in the GHG fluxes, N cycling, and the microbial community structure. In the present study, the below-ground processes were modified after three years of moderate O3 enhancement. A tentative conclusion is that a decrease in N availability may have feedback effects on plant growth and competition and affect the N cycling of the whole meadow ecosystem. Ecosystem level changes occur slowly, and multiplication of the responses might be expected in the long run.

Effects of recreational use and fragmentation on the understorey vegetation and soil microbial communities of urban forests in southern Finland

The impacts of fragmentation and recreational use on the hemiboreal urban forest understorey vegetation and the microbial community of the humus layer (the phospholipid fatty acid (PLFA) pattern, microbial biomass and microbial activity, measured as basal respiration) were examined in the greater Helsinki area, southern Finland. Trampling tolerance of 1) herb-rich OMT, 2) mesic MT, and 3) sub-xeric VT forests (in decreasing order of fertility) was studied by comparing relative understorey vegetation cover (urban/untrampled reference ratio) of the three forest types. The trampling tolerance of forest vegetation increased with the productivity of the site (sub-xeric < mesic < herb-rich). Wear of understorey vegetation correlated positively with the number of residents (i.e., recreational pressure) around the forest patch. An increase of 15000 residents within a radius of 1 km around a forest patch was associated with ca. 30% decrease in the relative understorey vegetation cover. The cover of dwarf shrub Vaccinium myrtillus in particular decreased with increasing levels of wear. The cover of mosses in urban forests was less than half of that in untrampled reference areas. Cover of tree saplings, mainly Sorbus aucuparia, and some resilient herbs was higher than in the reference areas. In small urban forest fragments, broad-leaved trees, grasses and herbs were more abundant and mosses were scarcer than in larger urban forest areas. Thus, due to trampling and edge effects, resilient herb and grass species are replacing sensitive dwarf shrubs, mosses and lichens in urban forests. Differences in the soil microbial community structure were found between paths and untrampled areas and the effects of paths extended more than one meter from the paths. Paths supported approximately 25-30% higher microbial biomass with a transition zone of at least 1 m from the path edge. However, microbial activity per unit of biomass was lower on paths than in untrampled areas. Furthermore, microbial biomass and activity were 30-45% lower at the first 20 m into the forest fragments, due to low moisture content of humus near the edge. The decreased microbial activity detected at forest edges and paths implies decreased litter decomposition rates, and thus, a change in nutrient cycling. Changes in the decomposition and nutrient supply may in turn affect the diversity and function of plant communities in urban forests. Keywords: boreal forest vegetation, edge effects, phospholipid fatty acids, trampling, urban woodlands, wear