Silvicultural decisions based on simulation-optimization systems

Forest management is facing new challenges under climate change. By adjusting thinning regimes, conventional forest management can be adapted to various objectives of utilization of forest resources, such as wood quality, forest bioenergy, and carbon sequestration. This thesis aims to develop and apply a simulation-optimization system as a tool for an interdisciplinary understanding of the interactions between wood science, forest ecology, and forest economics. In this thesis, the OptiFor software was developed for forest resources management. The OptiFor simulation-optimization system integrated the process-based growth model PipeQual, wood quality models, biomass production and carbon emission models, as well as energy wood and commercial logging models into a single optimization model. Osyczka s direct and random search algorithm was employed to identify optimal values for a set of decision variables. The numerical studies in this thesis broadened our current knowledge and understanding of the relationships between wood science, forest ecology, and forest economics. The results for timber production show that optimal thinning regimes depend on site quality and initial stand characteristics. Taking wood properties into account, our results show that increasing the intensity of thinning resulted in lower wood density and shorter fibers. The addition of nutrients accelerated volume growth, but lowered wood quality for Norway spruce. Integrating energy wood harvesting into conventional forest management showed that conventional forest management without energy wood harvesting was still superior in sparse stands of Scots pine. Energy wood from pre-commercial thinning turned out to be optimal for dense stands. When carbon balance is taken into account, our results show that changing carbon assessment methods leads to very different optimal thinning regimes and average carbon stocks. Raising the carbon price resulted in longer rotations and a higher mean annual increment, as well as a significantly higher average carbon stock over the rotation.

Gene flow from transgenic plant populations: models and applications for risk assessment

The future use of genetically modified (GM) plants in food, feed and biomass production requires a careful consideration of possible risks related to the unintended spread of trangenes into new habitats. This may occur via introgression of the transgene to conventional genotypes, due to cross-pollination, and via the invasion of GM plants to new habitats. Assessment of possible environmental impacts of GM plants requires estimation of the level of gene flow from a GM population. Furthermore, management measures for reducing gene flow from GM populations are needed in order to prevent possible unwanted effects of transgenes on ecosystems. This work develops modeling tools for estimating gene flow from GM plant populations in boreal environments and for investigating the mechanisms of the gene flow process. To describe spatial dimensions of the gene flow, dispersal models are developed for the local and regional scale spread of pollen grains and seeds, with special emphasis on wind dispersal. This study provides tools for describing cross-pollination between GM and conventional populations and for estimating the levels of transgenic contamination of the conventional crops. For perennial populations, a modeling framework describing the dynamics of plants and genotypes is developed, in order to estimate the gene flow process over a sequence of years. The dispersal of airborne pollen and seeds cannot be easily controlled, and small amounts of these particles are likely to disperse over long distances. Wind dispersal processes are highly stochastic due to variation in atmospheric conditions, so that there may be considerable variation between individual dispersal patterns. This, in turn, is reflected to the large amount of variation in annual levels of cross-pollination between GM and conventional populations. Even though land-use practices have effects on the average levels of cross-pollination between GM and conventional fields, the level of transgenic contamination of a conventional crop remains highly stochastic. The demographic effects of a transgene have impacts on the establishment of trangenic plants amongst conventional genotypes of the same species. If the transgene gives a plant a considerable fitness advantage in comparison to conventional genotypes, the spread of transgenes to conventional population can be strongly increased. In such cases, dominance of the transgene considerably increases gene flow from GM to conventional populations, due to the enhanced fitness of heterozygous hybrids. The fitness of GM plants in conventional populations can be reduced by linking the selectively favoured primary transgene to a disfavoured mitigation transgene. Recombination between these transgenes is a major risk related to this technique, especially because it tends to take place amongst the conventional genotypes and thus promotes the establishment of invasive transgenic plants in conventional populations.

Meadow plant growth and competition under elevated ozone and carbon dioxide

The main aim of my thesis project was to assess the impact of elevated ozone (O3) and carbon dioxide (CO2) on the growth, competition and community of meadow plants in northern Europe. The thesis project consisted of three separate O3 and CO2 exposure experiments that were conducted as open-top-chamber (OTC) studies at Jokioinen, SW Finland, and a smaller-scale experiment with different availabilities of resources in greenhouses in Helsinki. The OTC experiments included a competition experiment with two- and three-wise interactions, a mesocosm-scale meadow community with a large number of species, and a pot experiment that assessed intraspecific differences of Centaurea jacea ecotypes. The studied lowland hay meadow proved to be an O3-sensitive biotope, as the O3 concentrations used (40-50 ppb) were moderate, and yet, six out of nine species (Campanula rotundifolia, Centaurea jacea, Fragaria vesca, Ranunculus acris, Trifolium medium, Vicia cracca) showed either significant reductions in biomass or reproductive development, visible O3 injury or any two as a response to elevated O3. The plant species and ecotypes exhibited large intra- and interspecific variation in their response to O3, but O3 and CO2 concentrations did not cause changes in their interspecific competition or in community composition. However, the largest O3-induced growth reductions were seen in the least abundant species (C. rotundifolia and F. vesca), which may indicate O3-induced suppression of weak competitors. The overall effects of CO2 were relatively small and mainly restricted to individual species and several measured variables. Based on the present studies, most of the deleterious effects of tropospheric O3 are not diminished by a moderate increase in CO2 under low N availability, and variation exists between different species and variables. The present study indicates that the growth of several herb species decreases with increasing atmospheric O3 concentrations, and that these changes may pose a threat to the biodiversity of meadows. Ozone-induced reductions in the total community biomass production and N pool are likely to have important consequences for the nutrient cycling of the ecosystem.

The effects of environmental factors on biomass and microcystin production by the freshwater cyanobacterial genera Microcystis and Anabaena

Several cyanobacterial genera produce the hepatotoxins, microcystins. Microcystins are produced only in cells that have microcystin synthetase gene (mcy) clusters, which encode enzyme complexes involved in microcystin biosynthesis. Microcystin-producing and nonmicrocystin-producing genotypes of single cyanobacterial genus may occur simultaneously in situ. Previously, the effects of environmental factors on the growth and microcystin production of cyanobacteria have mainly been studied by means of isolated cyanobacteria cultures in the laboratory. Studies in the field have been difficult, owing to the lack of methods to identify and quantify the different genotypes. In this study, genus-specific microcystin synthetase E (mcyE) gene primers were designed and a method to identify and quantify the mcyE copy numbers was developed and used in situ. Microcystis and Anabaena mcyE genes were observed in two Finnish lakes. Microcystis appeared to be the most abundant microcystin producer in Lake Tuusulanjärvi and in one basin of Lake Hiidenvesi. Because the most potent microcystin-producing genus of a lake can be identified, it will be possible in the future to design genus-targeted strategies for lake restoration. Effects of P and N concentrations on the biomass of microcystin-producing and nonmicrocystin-producing Microcystis strains and an Anabaena strain were studied in cultures. P and N concentrations and their combined effect increased cyanobacterial biomass of all Microcystis strains. The biomass of microcystin-producing Microcystis was higher than that of nonmicrocystin-producing strains at high nutrient concentrations. The P concentration increased Anabaena biomass, but the effect of N concentration was statistically insignificant for growth yield, probably due to the ability of the genus to fix molecular N2. P and N concentrations and combined nutrients caused an increase in cellular microcystin concentrations of the Microcystis strain cultivated in chemostat cultures. Cyanobacteria are able to hydrolyse nutrients from organic matter through extracellular enzyme activities. Leucine aminopeptidase (LAP) activity was observed in an axenic N2-fixing Anabaena strain grown in batch cultures. The P concentration caused a statistically significant increase in LAP activity, whereas the effect of N concentration was insignificant. The highest LAP activities were observed in the most eutrophic basins of Lake Hiidenvesi. LAP activity probably originated mostly from attached heterotrophic bacteria and less from cyanobacteria.

Forest inventory-based large-scale forest biomass and carbon budget assessment: new enhanced methods and use of remote sensing for verification

In recent years, concern has arisen over the effects of increasing carbon dioxide (CO2) in the earth's atmosphere due to the burning of fossil fuels. One way to mitigate increase in atmospheric CO2 concentration and climate change is carbon sequestration to forest vegeta-tion through photosynthesis. Comparable regional scale estimates for the carbon balance of forests are therefore needed for scientific and political purposes. The aim of the present dissertation was to improve methods for quantifying and verifying inventory-based carbon pool estimates of the boreal forests in the mineral soils. Ongoing forest inventories provide a data based on statistically sounded sampling for estimating the level of carbon stocks and stock changes, but improved modelling tools and comparison of methods are still needed. In this dissertation, the entire inventory-based large-scale forest carbon stock assessment method was presented together with some separate methods for enhancing and comparing it. The enhancement methods presented here include ways to quantify the biomass of understorey vegetation as well as to estimate the litter production of needles and branches. In addition, the optical remote sensing method illustrated in this dis-sertation can be used to compare with independent data. The forest inventory-based large-scale carbon stock assessment method demonstrated here provided reliable carbon estimates when compared with independent data. Future ac-tivity to improve the accuracy of this method could consist of reducing the uncertainties regarding belowground biomass and litter production as well as the soil compartment. The methods developed will serve the needs for UNFCCC reporting and the reporting under the Kyoto Protocol. This method is principally intended for analysts or planners interested in quantifying carbon over extensive forest areas.

Radiative transfer, interception and scattering in coniferous forests : models and applications for production ecology and remote sensing

This work develops methods to account for shoot structure in models of coniferous canopy radiative transfer. Shoot structure, as it varies along the light gradient inside canopy, affects the efficiency of light interception per unit needle area, foliage biomass, or foliage nitrogen. The clumping of needles in the shoot volume also causes a notable amount of multiple scattering of light within coniferous shoots. The effect of shoot structure on light interception is treated in the context of canopy level photosynthesis and resource use models, and the phenomenon of within-shoot multiple scattering in the context of physical canopy reflectance models for remote sensing purposes. Light interception. A method for estimating the amount of PAR (Photosynthetically Active Radiation) intercepted by a conifer shoot is presented. The method combines modelling of the directional distribution of radiation above canopy, fish-eye photographs taken at shoot locations to measure canopy gap fraction, and geometrical measurements of shoot orientation and structure. Data on light availability, shoot and needle structure and nitrogen content has been collected from canopies of Pacific silver fir (Abies amabilis (Dougl.) Forbes) and Norway spruce (Picea abies (L.) Karst.). Shoot structure acclimated to light gradient inside canopy so that more shaded shoots have better light interception efficiency. Light interception efficiency of shoots varied about two-fold per needle area, about four-fold per needle dry mass, and about five-fold per nitrogen content. Comparison of fertilized and control stands of Norway spruce indicated that light interception efficiency is not greatly affected by fertilization. Light scattering. Structure of coniferous shoots gives rise to multiple scattering of light between the needles of the shoot. Using geometric models of shoots, multiple scattering was studied by photon tracing simulations. Based on simulation results, the dependence of the scattering coefficient of shoot from the scattering coefficient of needles is shown to follow a simple one-parameter model. The single parameter, termed the recollision probability, describes the level of clumping of the needles in the shoot, is wavelength independent, and can be connected to previously used clumping indices. By using the recollision probability to correct for the within-shoot multiple scattering, canopy radiative transfer models which have used leaves as basic elements can use shoots as basic elements, and thus be applied for coniferous forests. Preliminary testing of this approach seems to explain, at least partially, why coniferous forests appear darker than broadleaved forests in satellite data.

Engineering the pentose phosphate pathway of Saccharomyces cerevisiae for production of ethanol and xylitol

The baker s yeast Saccharomyces cerevisiae has a long tradition in alcohol production from D-glucose of e.g. starch. However, without genetic modifications it is unable to utilise the 5-carbon sugars D-xylose and L arabinose present in plant biomass. In this study, one key metabolic step of the catabolic D-xylose pathway in recombinant D-xylose-utilising S. cerevisiae strains was studied. This step, carried out by xylulokinase (XK), was shown to be rate-limiting, because overexpression of the xylulokinase-encoding gene XKS1 increased both the specific ethanol production rate and the yield from D xylose. In addition, less of the unwanted side product xylitol was produced. Recombinant D-xylose-utilizing S. cerevisiae strains have been constructed by expressing the genes coding for the first two enzymes of the pathway, D-xylose reductase (XR) and xylitol dehydrogenase (XDH) from the D-xylose-utilising yeast Pichia stipitis. In this study, the ability of endogenous genes of S. cerevisiae to enable D-xylose utilisation was evaluated. Overexpression of the GRE3 gene coding for an unspecific aldose reductase and the ScXYL2 gene coding for a xylitol dehydrogenase homologue enabled growth on D-xylose in aerobic conditions. However, the strain with GRE3 and ScXYL2 had a lower growth rate and accumulated more xylitol compared to the strain with the corresponding enzymes from P. stipitis. Use of the strictly NADPH-dependent Gre3p instead of the P. stipitis XR able to utilise both NADH and NADPH leads to a more severe redox imbalance. In a S. cerevisiae strain not engineered for D-xylose utilisation the presence of D-xylose increased xylitol dehydrogenase activity and the expression of the genes SOR1 or SOR2 coding for sorbitol dehydrogenase. Thus, D-xylose utilisation by S. cerevisiae with activities encoded by ScXYL2 or possibly SOR1 or SOR2, and GRE3 is feasible, but requires efficient redox balance engineering. Compared to D-xylose, D-glucose is a cheap and readily available substrate and thus an attractive alternative for xylitol manufacture. In this study, the pentose phosphate pathway (PPP) of S. cerevisiae was engineered for production of xylitol from D-glucose. Xylitol was formed from D-xylulose 5-phosphate in strains lacking transketolase activity and expressing the gene coding for XDH from P. stipitis. In addition to xylitol, ribitol, D-ribose and D-ribulose were also formed. Deletion of the xylulokinase-encoding gene increased xylitol production, whereas the expression of DOG1 coding for sugar phosphate phosphatase increased ribitol, D-ribose and D-ribulose production. Strains lacking phosphoglucose isomerase (Pgi1p) activity were shown to produce 5 carbon compounds through PPP when DOG1 was overexpressed. Expression of genes encoding glyceraldehyde 3-phosphate dehydrogenase of Bacillus subtilis, GapB, or NAD-dependent glutamate dehydrogenase Gdh2p of S. cerevisiae, altered the cellular redox balance and enhanced growth of pgi1 strains on D glucose, but co-expression with DOG1 reduced growth on higher D-glucose concentrations. Strains lacking both transketolase and phosphoglucose isomerase activities tolerated only low D-glucose concentrations, but the yield of 5-carbon sugars and sugar alcohols on D-glucose was about 50% (w/w).

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.