Novel Insights into the Molecular Genetic Background of Colorectal Tumors

Many of the genes predisposing to highly penetrant colorectal cancer (CRC) syndromes, including hereditary non-polyposis colorectal cancer (MLH1, MSH2, MSH6, PMS2), familial adenomatous polyposis (APC), Peutz-Jeghers syndrome (LKB1), juvenile polyposis (SMAD4, BMPR1A), MYH-associated polyposis (MYH), and Cowden syndrome (PTEN) have already been discovered. Identification of these genes has allowed a more precise classification of the hereditary CRC syndromes and provided a means for predictive genetic testing and surveillance. Some of the genes are also involved in sporadic cancer forms, and therefore the investigation of the rare CRC syndromes has been a breakthrough for general cancer research. Despite the accumulating knowledge on hereditary cancer syndromes, a significant number of familial CRCs remain molecularly unexplained after genetic testing, reflecting the possibility of other predisposing genes or existence of novel syndromes. Moreover, genetic variants conferring low-penetrance risk are still largely unknown. In this study, we examined the role of some new high- and low-penetrance alleles on CRC predisposition. We identified disease causing MYH mutations in a subset (9%) of patients with APC and AXIN2 mutation negative adenomatous polyposis. Due to differences in the pattern of inheritance and clinical manifestation, screening for mutations in MYH is beneficial in view of genetic counselling and surveillance. A novel functionally deficient MYH founder mutation A459D was identified in the Finnish population, and this finding had immediate clinical implications for genetic counselling of at risk families. Many patients with hamartomatous polyposis remain without molecular diagnosis due to atypical phenotypes. We therefore sought to classify 49 patients with unexplained hamartomatous or hyperplastic/mixed polyposis by extensive molecular analyses of PTEN, LKB1, BMPR1A, SMAD4, ENG, BRAF, MYH, and BHD along with revision of polyp histology. Mutations were identified in 11/49 (22%) of the patients. In 6 cases the molecular diagnosis was re-classified guiding surveillance and decisions for prophylactic surgery. Re-evaluation of polyp histology with subsequent more accurate selection of candidate gene analyses is beneficial and can be recommended for patients with unexplained polyposis. Furthermore, germline mutations in ENG underlying juvenile polyposis were described for the first time, characterizing a possible novel genetically defined form of hereditary CRC. Association analyses on two putative low-penetrance alleles, NOD2 3020insC and MDM2 SNP309 were performed in a population-based series of 1042 Finnish CRC patients and in cancer-free controls. In contrast to previous results, NOD2 3020insC did not associate with CRC or age at disease onset in the Finnish population. These data suggest that NOD2 3020insC alone might not be sufficient for CRC predisposition. MDM2 SNP309 was as common in the CRC cohort as in the healthy controls. Interesting trends, however, were observed, which after correction for multiple testing did not reach statistical significance. SNP309 was more common in female CRC patients and a trend towards an earlier age at disease onset was observed in women with SNP309. Subsequent studies have supported this observation and SNP309 could affect gender- or hormone-related tumorigenesis. Finally, a large-scale unbiased effort was designed to characterize the complete mutatome of CRC with microsatellite instability (MSI). Using an approach combining expression microarray and genome database searches, we were able to identify putative MSI target genes. Further characterization of one of the genes suggested that it might play a role also in microsatellite stable CRC and Peutz-Jeghers syndrome pathogenesis.

Novel prognostic factors in chronic myeloid leukemia

Chronic myeloid leukemia (CML) is a malignant clonal blood disease that originates from a pluripotent hematopoietic stem cell. The cytogenetic hallmark of CML, the Philadelphia chromosome (Ph), is formed as a result of reciprocal translocation between chromosomes 9 and 22, which leads to a formation of a chimeric BCR-ABL fusion gene. The BCR-ABL protein is a constitutively active tyrosine kinase that changes the adhesion properties of cells, constitutively activates mitogenic signaling, enhances cell proliferation and reduces apoptosis. This results in leukemic growth and the clinical disease, CML. With the advent of targeted therapies against the BCR-ABL fusion protein, the treatment of CML has changed considerably during the recent decade. In this thesis, the clinical significance of different diagnostic methods and new prognostic factors in CML have been assessed. First, the association between two different methods for measuring CML disease burden (the RQ-PCR and the high mitotic index metaphase FISH) was assessed in bone marrow and peripheral blood samples. The correlation between positive RQ-PCR and metaphase FISH samples was high. However, RQ-PCR was more sensitive and yielded measurable transcripts in 40% of the samples that were negative by metaphase FISH. The study established a laboratory-specific conversion factor for setting up the International Scale when standardizing RQ-PCR measurements. Secondly, the amount of minimal residual disease (MRD) after allogeneic hematopoietic stem cell transplantation (alloHSCT) was determined. For this, metaphase FISH was done for the bone marrow samples of 102 CML patients. Most (68%), had no residual cells during the entire follow-up time. Some (12 %) patients had minor (<1%) MRD which decreased even further with time, whereas 19% had a progressive rise in MRD that exceeded 1% or had more than 1% residual cells when first detected. Residual cells did not become eradicated spontaneously if the frequency of Ph+ cells exceeded 1% during follow-up. Next, the impact of deletions in the derivative chromosome 9, was examined. Deletions were observed in 15% of the CML patients who later received alloHSCT. After alloHSCT, there was no difference in the total relapse rate in patients with or without deletions. Nor did the estimates of overall survival, transplant-related mortality, leukemia-free survival and relapse-free time show any difference between these groups. When conventional treatment regimens are used, the der(9) status could be an important criterion, in conjunction with other prognostic factors, when allogeneic transplantation is considered. The significance of der(9) deletions for patients treated with tyrosine kinase inhibitors is not clear and requires further investigation. In addition to the der(9) status of the patient, the significance of bone marrow lymphocytosis as a prognostic factor in CML was assessed. Bone marrow lymphocytosis during imatinib therapy was a positive predictive factor and heralded optimal response. When combined with major cytogenetic response at three months of treatment, bone marrow lymphocytosis predicted a prognostically important major molecular response at 18 months of imatinib treatment. Although the validation of these findings is warranted, the determination of the bone marrow lymphocyte count could be included in the evaluation of early response to imatinib treatment already now. Finally, BCR-ABL kinase domain mutations were studied in CML patients resistant against imatinib treatment. Point mutations detected in the kinase domain were the same as previously reported, but other sequence variants, e.g. deletions or exon splicing, were also found. The clinical significance of the other variations remains to be determined.

Management of multi-scale forest resource data over time

During the last decades there has been a global shift in forest management from a focus solely on timber management to ecosystem management that endorses all aspects of forest functions: ecological, economic and social. This has resulted in a shift in paradigm from sustained yield to sustained diversity of values, goods and benefits obtained at the same time, introducing new temporal and spatial scales into forest resource management. The purpose of the present dissertation was to develop methods that would enable spatial and temporal scales to be introduced into the storage, processing, access and utilization of forest resource data. The methods developed are based on a conceptual view of a forest as a hierarchically nested collection of objects that can have a dynamically changing set of attributes. The temporal aspect of the methods consists of lifetime management for the objects and their attributes and of a temporal succession linking the objects together. Development of the forest resource data processing method concentrated on the extensibility and configurability of the data content and model calculations, allowing for a diverse set of processing operations to be executed using the same framework. The contribution of this dissertation to the utilisation of multi-scale forest resource data lies in the development of a reference data generation method to support forest inventory methods in approaching single-tree resolution.

Landscape-scale soil erosion modelling and ecological restoration for a mountainous watershed in Sichuan, China

Sichuanissa Tiibetin ylängön metsäkato on pysähtynyt mutta eroosio-ongelmat jatkuvat Viikin tropiikki-instituutin tutkija Ping ZHOU kartoitti trooppisen metsänhoidon alaan kuuluvassa väitöskirjatyössään maaperän eroosioalttiutta ja sen riippuvuutta metsäkasvillisuudesta Jangtsen tärkeää sivuhaaraa Min-jokea ympäröivällä n. 7400 neliökilometrin suuruisella valuma-alueella Sichuanin Aba-piirikunnassa. Aineistonaan hän käytti muun muassa satelliittikartoitustietoja ja mittaustuloksia yli 600 maastokoealalta. Tutkimuksen nimi suomeksi on "Maaperän eroosion mallinnus ja vuoristoisen valuma-alueen ekologinen ennallistaminen Sichuanissa Kiinassa". Aikaisempien tutkimusten perusteella oli tiedossa että metsien häviäminen tällä alueella pysähtyi jo 1980-luvun alussa. Sen jälkeen on metsien pinta-ala hitaasti kasvanut etupäässä sen vuoksi, että teollinen puunhakkuu luonnonmetsissä kiellettiin kokonaan v. 1998 ja 25 astetta jyrkemmillä rinteillä myös maatalouden harjoittaminen on saatu lopetetuksi viljelijöille tarjottujen taloudellisten houkuttimien avulla. Täten myös pelto- ja laidunmaata on voitu ennallistaa metsäksi. Ping Zhou pystyi jakamaan 5700 metrin korkeuteen saakka kohoavan vuoristoalueen eroosioalttiudeltaan erilaisiin vyöhykkeisiin rinteen kaltevuuden, sademäärän, kasvipeitteen ja maalajin perusteella. Noin 15 prosentilla tutkitun valuma-alueen pinta-alasta, lähinnä Min-joen pääuomaa ympäröivillä jyrkillä rinteillä, eroosioriski oli suuri tai erittäin suuri. Eri tyyppisellä kasvillisuudella oli hyvin erilainen vaikutus eroosioalttiuteen, ja myös alueen sijainti vuoriston eri korkeuksilla vaikutti eroosioon. Säästyneet lähes luonnontilaiset havumetsät, joita on etupäässä vuoriston ylimmissä osissa 2600-4000 metrin korkeudella, edistävät tehokkaasti metsän luontaista uudistumista ja levittäytymistä vaurioituneille alueille. Säilyneiden metsien puulajikoostumus antoi tutkimuksessa mahdollisuuden ennustaa metsien tulevaa kehitystä koko tutkitulla valuma-alueella sen eri korkeusvyöhykkeissä ja eri maaperätyypeillä. Ennallistamisen kannalta ongelmallisimpia olivat alueet joilta metsäpeite oli lähinnä puiden teollisen hakkuun vuoksi kokonaan hävinnyt ja joilla maaperä yleisesti oli eroosion pahoin kuluttama. Näillä alueilla ei ole tehty juuri mitään uudistamis- tai ennallistamistoimenpiteitä. Niillä metsien ennallistaminen vaatii myös puiden tai pensaiden istuttamista. Tähän sopivia ovat erityisesti ilmakehän typpeä sitovat lajit, joista alueella kasvaa luontaisena mm. sama tyrnilaji joka esiintyy myös Suomessa. Työssä tutkittiin yli kahdeksankymmenen paikallisen luontaisen puulajin (joista peräti noin kolmannes on havupuulajeja) ekologisia ominaisuuksia ja soveltuvuutta metsien ennallistamiseen. Avainasemassa työn onnistumisen kannalta ovat nyt paikalliset asukkaat, joiden maankäytön muutokset ovat jo selvästi edistänet luonnonmetsän ennalleen palautumista. Suomen Akatemia rahoitti vuosina 2004-2006 VITRI:n tutkimushanketta, josta Ping Zhou'n väitöskirjatyö muodosti keskeisen osan. Kenttätyö Sichuanissa avasi mahdollisuuden hedelmälliseen monitieteiseen yhteistyöhön ja tutkijavaihtoon Kiinan tiedeakatemian alaisen Chengdun biologiainstituutin (CIB) kanssa; tämä tieteellinen kanssakäyminen jatkuu edelleen.

Plant species diversity of buffer zones in agricultural landscapes: in search of determinants from the local to regional scale

Buffer zones are vegetated strip-edges of agricultural fields along watercourses. As linear habitats in agricultural ecosystems, buffer strips dominate and play a leading ecological role in many areas. This thesis focuses on the plant species diversity of the buffer zones in a Finnish agricultural landscape. The main objective of the present study is to identify the determinants of floral species diversity in arable buffer zones from local to regional levels. This study was conducted in a watershed area of a farmland landscape of southern Finland. The study area, Lepsämänjoki, is situated in the Nurmijärvi commune 30 km to the north of Helsinki, Finland. The biotope mosaics were mapped in GIS. A total of 59 buffer zones were surveyed, of which 29 buffer strips surveyed were also sampled by plot. Firstly, two diversity components (species richness and evenness) were investigated to determine whether the relationship between the two is equal and predictable. I found no correlation between species richness and evenness. The relationship between richness and evenness is unpredictable in a small-scale human-shaped ecosystem. Ordination and correlation analyses show that richness and evenness may result from different ecological processes, and thus should be considered separately. Species richness correlated negatively with phosphorus content, and species evenness correlated negatively with the ratio of organic carbon to total nitrogen in soil. The lack of a consistent pattern in the relationship between these two components may be due to site-specific variation in resource utilization by plant species. Within-habitat configuration (width, length, and area) were investigated to determine which is more effective for predicting species richness. More species per unit area increment could be obtained from widening the buffer strip than from lengthening it. The width of the strips is an effective determinant of plant species richness. The increase in species diversity with an increase in the width of buffer strips may be due to cross-sectional habitat gradients within the linear patches. This result can serve as a reference for policy makers, and has application value in agricultural management. In the framework of metacommunity theory, I found that both mass effect(connectivity) and species sorting (resource heterogeneity) were likely to explain species composition and diversity on a local and regional scale. The local and regional processes were interactively dominated by the degree to which dispersal perturbs local communities. In the lowly and intermediately connected regions, species sorting was of primary importance to explain species diversity, while the mass effect surpassed species sorting in the highly connected region. Increasing connectivity in communities containing high habitat heterogeneity can lead to the homogenization of local communities, and consequently, to lower regional diversity, while local species richness was unrelated to the habitat connectivity. Of all species found, Anthriscus sylvestris, Phalaris arundinacea, and Phleum pretense significantly responded to connectivity, and showed high abundance in the highly connected region. We suggest that these species may play a role in switching the force from local resources to regional connectivity shaping the community structure. On the landscape context level, the different responses of local species richness and evenness to landscape context were investigated. Seven landscape structural parameters served to indicate landscape context on five scales. On all scales but the smallest scales, the Shannon-Wiener diversity of land covers (H') correlated positively with the local richness. The factor (H') showed the highest correlation coefficients in species richness on the second largest scale. The edge density of arable field was the only predictor that correlated with species evenness on all scales, which showed the highest predictive power on the second smallest scale. The different predictive power of the factors on different scales showed a scaledependent relationship between the landscape context and local plant species diversity, and indicated that different ecological processes determine species richness and evenness. The local richness of species depends on a regional process on large scales, which may relate to the regional species pool, while species evenness depends on a fine- or coarse-grained farming system, which may relate to the patch quality of the habitats of field edges near the buffer strips. My results suggested some guidelines of species diversity conservation in the agricultural ecosystem. To maintain a high level of species diversity in the strips, a high level of phosphorus in strip soil should be avoided. Widening the strips is the most effective mean to improve species richness. Habitat connectivity is not always favorable to species diversity because increasing connectivity in communities containing high habitat heterogeneity can lead to the homogenization of local communities (beta diversity) and, consequently, to lower regional diversity. Overall, a synthesis of local and regional factors emerged as the model that best explain variations in plant species diversity. The studies also suggest that the effects of determinants on species diversity have a complex relationship with scale.

Protein cross-linking with oxidative enzymes and transglutaminase : Effects in meat protein systems

The effects of tyrosinase, laccase and transglutaminase (TG) were studied in different meat protein systems. The study was focused on the effects of the enzymes on the gel formation properties of myofibrils, and on the textural and water-holding properties of the heated meat systems. The cross-linking efficiency of a novel Trichoderma reesei tyrosinase was compared to that of the commercial Agaricus bisporus tyrosinase. Trichoderma tyrosinase was found to be superior compared to the Agaricus enzyme in its protein cross-linking efficiency and in the incorporation of a small molecule into a complex proteinaceous substrate. Tyrosinase, laccase and TG all polymerised myofibrillar proteins, but laccase was also found to cause protein fragmentation. A positive connection between covalent cross-link and gel formation was observed with tyrosinase and TG. Laccase was able to increase the gel formation only slightly. With an excessive laccase dosage the gel formation declined due to protein fragmentation. Tyrosinase, laccase and TG had different effects on the texture and water-holding of the heated chicken breast meat homogenates. Tyrosinase improved the firmness of the homogenate gels free of phosphate and with a low amount of meat. TG improved the firmness of all studied homogenates. Laccase weakened the gel firmness of the low-meat, low-salt and low-salt/phosphate homogenates and maintained the firmness on the control level in the homogenate free of phosphate. Tyrosinase was the only enzyme capable of reducing the weight loss in the homogenates containing a low amount of meat and a low amount of NaCl. TG was the only enzyme that could positively affect the firmness of the homogenate gel containing both low NaCl and phosphate amounts. In pilot scale the test products were made of coarsely ground chicken breast fillet with a moderate amount of salt. Increasining the amount of meat, salt and TG contents favoured the development of firmness of the test products. The evaporation loss decreased slightly along with increasing TG and NaCl amounts in the experimental conditions used, indicating a positive interaction between these two factors. In this work it was shown that tyrosinase, laccase and TG affected the same myofibrillar proteins, i.e. myosin and troponin T. However, these enzymes had distinguishable effects on the gel formation of a myofibril system as well as on the textural and water-holding properties of the finely ground meat homogenates, reflecting distinctions at least in the reaction mechanisms and target amino acid availability in the protein substrates for these enzymes.

Technicolor

In technicolor theories the scalar sector of the Standard Model is replaced by a strongly interacting sector. Although the Standard Model has been exceptionally successful, the scalar sector causes theoretical problems that make these theories seem an attractive alternative. I begin my thesis by considering QCD, which is the known example of strong interactions. The theory exhibits two phenomena: confinement and chiral symmetry breaking. I find the low-energy dynamics to be similar to that of the sigma models. Then I analyze the problems of the Standard Model Higgs sector, mainly the unnaturalness and triviality. Motivated by the example of QCD, I introduce the minimal technicolor model to resolve these problems. I demonstrate the minimal model to be free of anomalies and then deduce the main elements of its low-energy particle spectrum. I find the particle spectrum contains massless or very light technipions, and also technibaryons and techni-vector mesons with a high mass of over 1 TeV. Standard Model fermions remain strictly massless at this stage. Thus I introduce the technicolor companion theory of flavor, called extended technicolor. I show that the Standard Model fermions and technihadrons receive masses, but that they remain too light. I also discuss flavor-changing neutral currents and precision electroweak measurements. I then show that walking technicolor models partly solve these problems. In these models, contrary to QCD, the coupling evolves slowly over a large energy scale. This behavior adds to the masses so that even the light technihadrons are too heavy to be detected at current particle accelerators. Also all observed masses of the Standard Model particles can be generated, except for the bottom and top quarks. Thus it is shown in this thesis that, excluding the masses of third generation quarks, theories based on walking technicolor can in principle produce the observed particle spectrum.

The kinetics and reactivity of several polyatomic free radicals in reactions with NO2, O2, Cl2, and HCl

In this thesis, the kinetics of several alkyl, halogenated alkyl, and alkenyl free radical reactions with NO2, O2, Cl2, and HCl reactants were studied over a wide temperature range in time resolved conditions. Laser photolysis photoionisation mass spectrometer coupled to a flow reactor was the experimental method employed and this thesis present the first measurements performed with the experimental system constructed. During this thesis a great amount of work was devoted to the designing, building, testing, and improving the experimental apparatus. Carbon-centred free radicals were generated by the pulsed 193 or 248 nm photolysis of suitable precursors along the tubular reactor. The kinetics was studied under pseudo-first-order conditions using either He or N2 buffer gas. The temperature and pressure ranges employed were between 190 and 500 K, and 0.5 45 torr, respectively. The possible role of heterogeneous wall reactions was investigated employing reactor tubes with different sizes, i.e. to significantly vary the surface to volume ratio. In this thesis, significant new contributions to the kinetics of carbon-centred free radical reactions with nitrogen dioxide were obtained. Altogether eight substituted alkyl (CH2Cl, CHCl2, CCl3, CH2I, CH2Br, CHBr2, CHBrCl, and CHBrCH3) and two alkenyl (C2H3, C3H3) free radical reactions with NO2 were investigated as a function of temperature. The bimolecular rate coefficients of all these reactions were observed to possess negative temperature dependencies, while pressure dependencies were not noticed for any of these reactions. Halogen substitution was observed to moderately reduce the reactivity of substituted alkyl radicals in the reaction with NO2, while the resonance stabilisation of the alkenyl radical lowers its reactivity with respect to NO2 only slightly. Two reactions relevant to atmospheric chemistry, CH2Br + O2 and CH2I + O2, were also investigated. It was noticed that while CH2Br + O2 reaction shows pronounced pressure dependence, characteristic of peroxy radical formation, no such dependence was observed for the CH2I + O2 reaction. Observed primary products of the CH2I + O2 reaction were the I-atom and the IO radical. Kinetics of CH3 + HCl, CD3 + HCl, CH3 + DCl, and CD3 + DCl reactions were also studied. While all these reactions possess positive activation energies, in contrast to the other systems investigated in this thesis, the CH3 + HCl and CD3 + HCl reactions show a non-linear temperature dependency on the Arrhenius plot. The reactivity of substituted methyl radicals toward NO2 was observed to increase with decreasing electron affinity of the radical. The same trend was observed for the reactions of substituted methyl radicals with Cl2. It is proposed that interactions of frontier orbitals are responsible to these observations and Frontier Orbital Theory could be used to explain the observed reactivity trends of these highly exothermic reactions having reactant-like transition states.

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.

Spatial variation of climate and the impact of disturbances on local climate and forest recovery in northern Finland

The structure and function of northern ecosystems are strongly influenced by climate change and variability and by human-induced disturbances. The projected global change is likely to have a pronounced effect on the distribution and productivity of different species, generating large changes in the equilibrium at the tree-line. In turn, movement of the tree-line and the redistribution of species produce feedback to both the local and the regional climate. This research was initiated with the objective of examining the influence of natural conditions on the small-scale spatial variation of climate in Finnish Lapland, and to study the interaction and feedback mechanisms in the climate-disturbances-vegetation system near the climatological border of boreal forest. The high (1 km) resolution spatial variation of climate parameters over northern Finland was determined by applying the Kriging interpolation method that takes into account the effect of external forcing variables, i.e., geographical coordinates, elevation, sea and lake coverage. Of all the natural factors shaping the climate, the geographical position, local topography and altitude proved to be the determining ones. Spatial analyses of temperature- and precipitation-derived parameters based on a 30-year dataset (1971-2000) provide a detailed description of the local climate. Maps of the mean, maximum and minimum temperatures, the frost-free period and the growing season indicate that the most favourable thermal conditions exist in the south-western part of Lapland, around large water bodies and in the Kemijoki basin, while the coldest regions are in highland and fell Lapland. The distribution of precipitation is predominantly longitudinally dependent but with the definite influence of local features. The impact of human-induced disturbances, i.e., forest fires, on local climate and its implication for forest recovery near the northern timberline was evaluated in the Tuntsa area of eastern Lapland, damaged by a widespread forest fire in 1960 and suffering repeatedly-failed vegetation recovery since that. Direct measurements of the local climate and simulated heat and water fluxes indicated the development of a more severe climate and physical conditions on the fire-disturbed site. Removal of the original, predominantly Norway spruce and downy birch vegetation and its substitution by tundra vegetation has generated increased wind velocity and reduced snow accumulation, associated with a large variation in soil temperature and moisture and deep soil frost. The changed structural parameters of the canopy have determined changes in energy fluxes by reducing the latter over the tundra vegetation. The altered surface and soil conditions, as well as the evolved severe local climate, have negatively affected seedling growth and survival, leading to more unfavourable conditions for the reproduction of boreal vegetation and thereby causing deviations in the regional position of the timberline. However it should be noted that other factors, such as an inadequate seed source or seedbed, the poor quality of the soil and the intensive logging of damaged trees could also exacerbate the poor tree regeneration. In spite of the failed forest recovery at Tunsta, the position and composition of the timberline and tree-line in Finnish Lapland may also benefit from present and future changes in climate. The already-observed and the projected increase in temperature, the prolonged growing season, as well as changes in the precipitation regime foster tree growth and new regeneration, resulting in an advance of the timberline and tree-line northward and upward. This shift in the distribution of vegetation might be decelerated or even halted by local topoclimatic conditions and by the expected increase in the frequency of disturbances.

Ecological processes and large-scale climate relationships in northern coniferous forests

Ilmasto vaikuttaa ekologisiin prosesseihin eri tasoilla. Suuren mittakaavan ilmastoprosessit, yhdessä ilmakehän ja valtamerien kanssa, säätelevät paikallisia sääilmiöitä suurilla alueilla (mantereista pallopuoliskoihin). Tämä väistöskirja pyrkii selittämään kuinka suuren mittakaavan ilmasto on vaikuttanut tiettyihin ekologisiin prosesseihin pohjoisella havumetsäalueella. Valitut prosessit olivat puiden vuosilustojen kasvu, metsäpalojen esiintyminen ja vuoristomäntykovakuoriaisen aiheuttamat puukuolemat. Suuren mittakaavan ilmaston löydettiin vaikuttaneen näiden prosessien esiintymistiheyteen, kestoon ja levinneisyyteen keskeisten sään muuttujien välityksellä hyvin laajoilla alueilla. Tutkituilla prosesseilla oli vahva yhteys laajan mittakaavan ilmastoon. Yhteys on kuitenkin ollut hyvin dynaaminen ja muuttunut 1900-luvulla ilmastonmuutoksen aiheuttaessa muutoksia suuren mittakaavan ja alueellisten ilmastoprosessien välisiin sisäisiin suhteisiin.

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.

Computational Methods for Detecting Large-Scale Chromosome Rearrangements in SNP Data

Large-scale chromosome rearrangements such as copy number variants (CNVs) and inversions encompass a considerable proportion of the genetic variation between human individuals. In a number of cases, they have been closely linked with various inheritable diseases. Single-nucleotide polymorphisms (SNPs) are another large part of the genetic variance between individuals. They are also typically abundant and their measuring is straightforward and cheap. This thesis presents computational means of using SNPs to detect the presence of inversions and deletions, a particular variety of CNVs. Technically, the inversion-detection algorithm detects the suppressed recombination rate between inverted and non-inverted haplotype populations whereas the deletion-detection algorithm uses the EM-algorithm to estimate the haplotype frequencies of a window with and without a deletion haplotype. As a contribution to population biology, a coalescent simulator for simulating inversion polymorphisms has been developed. Coalescent simulation is a backward-in-time method of modelling population ancestry. Technically, the simulator also models multiple crossovers by using the Counting model as the chiasma interference model. Finally, this thesis includes an experimental section. The aforementioned methods were tested on synthetic data to evaluate their power and specificity. They were also applied to the HapMap Phase II and Phase III data sets, yielding a number of candidates for previously unknown inversions, deletions and also correctly detecting known such rearrangements.