Pitch discrimination in optimal and suboptimal acoustic environments: electroencephalographic, magnetoencephalographic, and behavioral evidence

Pitch discrimination is a fundamental property of the human auditory system. Our understanding of pitch-discrimination mechanisms is important from both theoretical and clinical perspectives. The discrimination of spectrally complex sounds is crucial in the processing of music and speech. Current methods of cognitive neuroscience can track the brain processes underlying sound processing either with precise temporal (EEG and MEG) or spatial resolution (PET and fMRI). A combination of different techniques is therefore required in contemporary auditory research. One of the problems in comparing the EEG/MEG and fMRI methods, however, is the fMRI acoustic noise. In the present thesis, EEG and MEG in combination with behavioral techniques were used, first, to define the ERP correlates of automatic pitch discrimination across a wide frequency range in adults and neonates and, second, they were used to determine the effect of recorded acoustic fMRI noise on those adult ERP and ERF correlates during passive and active pitch discrimination. Pure tones and complex 3-harmonic sounds served as stimuli in the oddball and matching-to-sample paradigms. The results suggest that pitch discrimination in adults, as reflected by MMN latency, is most accurate in the 1000-2000 Hz frequency range, and that pitch discrimination is facilitated further by adding harmonics to the fundamental frequency. Newborn infants are able to discriminate a 20% frequency change in the 250-4000 Hz frequency range, whereas the discrimination of a 5% frequency change was unconfirmed. Furthermore, the effect of the fMRI gradient noise on the automatic processing of pitch change was more prominent for tones with frequencies exceeding 500 Hz, overlapping with the spectral maximum of the noise. When the fundamental frequency of the tones was lower than the spectral maximum of the noise, fMRI noise had no effect on MMN and P3a, whereas the noise delayed and suppressed N1 and exogenous N2. Noise also suppressed the N1 amplitude in a matching-to-sample working memory task. However, the task-related difference observed in the N1 component, suggesting a functional dissociation between the processing of spatial and non-spatial auditory information, was partially preserved in the noise condition. Noise hampered feature coding mechanisms more than it hampered the mechanisms of change detection, involuntary attention, and the segregation of the spatial and non-spatial domains of working-memory. The data presented in the thesis can be used to develop clinical ERP-based frequency-discrimination protocols and combined EEG and fMRI experimental paradigms.

Technical efficiency of blood component preparation in blood centres of 10 European countries

Various reasons, such as ethical issues in maintaining blood resources, growing costs, and strict requirements for safe blood, have increased the pressure for efficient use of resources in blood banking. The competence of blood establishments can be characterized by their ability to predict the volume of blood collection to be able to provide cellular blood components in a timely manner as dictated by hospital demand. The stochastically varying clinical need for platelets (PLTs) sets a specific challenge for balancing supply with requests. Labour has been proven a primary cost-driver and should be managed efficiently. International comparisons of blood banking could recognize inefficiencies and allow reallocation of resources. Seventeen blood centres from 10 countries in continental Europe, Great Britain, and Scandinavia participated in this study. The centres were national institutes (5), parts of the local Red Cross organisation (5), or integrated into university hospitals (7). This study focused on the departments of blood component preparation of the centres. The data were obtained retrospectively by computerized questionnaires completed via Internet for the years 2000-2002. The data were used in four original articles (numbered I through IV) that form the basis of this thesis. Non-parametric data envelopment analysis (DEA, II-IV) was applied to evaluate and compare the relative efficiency of blood component preparation. Several models were created using different input and output combinations. The focus of comparisons was on the technical efficiency (II-III) and the labour efficiency (I, IV). An empirical cost model was tested to evaluate the cost efficiency (IV). Purchasing power parities (PPP, IV) were used to adjust the costs of the working hours and to make the costs comparable among countries. The total annual number of whole blood (WB) collections varied from 8,880 to 290,352 in the centres (I). Significant variation was also observed in the annual volume of produced red blood cells (RBCs) and PLTs. The annual number of PLTs produced by any method varied from 2,788 to 104,622 units. In 2002, 73% of all PLTs were produced by the buffy coat (BC) method, 23% by aphaeresis and 4% by the platelet-rich plasma (PRP) method. The annual discard rate of PLTs varied from 3.9% to 31%. The mean discard rate (13%) remained in the same range throughout the study period and demonstrated similar levels and variation in 2003-2004 according to a specific follow-up question (14%, range 3.8%-24%). The annual PLT discard rates were, to some extent, associated with production volumes. The mean RBC discard rate was 4.5% (range 0.2%-7.7%). Technical efficiency showed marked variation (median 60%, range 41%-100%) among the centres (II). Compared to the efficient departments, the inefficient departments used excess labour resources (and probably) production equipment to produce RBCs and PLTs. Technical efficiency tended to be higher when the (theoretical) proportion of lost WB collections (total RBC+PLT loss) from all collections was low (III). The labour efficiency varied remarkably, from 25% to 100% (median 47%) when working hours were the only input (IV). Using the estimated total costs as the input (cost efficiency) revealed an even greater variation (13%-100%) and overall lower efficiency level compared to labour only as the input. In cost efficiency only, the savings potential (observed inefficiency) was more than 50% in 10 departments, whereas labour and cost savings potentials were both more than 50% in six departments. The association between department size and efficiency (scale efficiency) could not be verified statistically in the small sample. In conclusion, international evaluation of the technical efficiency in component preparation departments revealed remarkable variation. A suboptimal combination of manpower and production output levels was the major cause of inefficiency, and the efficiency did not directly relate to production volume. Evaluation of the reasons for discarding components may offer a novel approach to study efficiency. DEA was proven applicable in analyses including various factors as inputs and outputs. This study suggests that analytical models can be developed to serve as indicators of technical efficiency and promote improvements in the management of limited resources. The work also demonstrates the importance of integrating efficiency analysis into international comparisons of blood banking.

Nephrin and the Pathogenesis of Nephropathy - Emphasis on Type I Diabetes

End-stage renal disease is an increasingly common pathologic condition, with a current incidence of 87 per million inhabitants in Finland. It is the end point of various nephropathies, most common of which is the diabetic nephropathy. This thesis focuses on exploring the role of nephrin in the pathogenesis of diabetic nephropathy. Nephrin is a protein of the glomerular epithelial cell, or podocyte, and it appears to have a crucial function as a component of the filtration slit diaphragm in the kidney glomeruli. Mutations in the nephrin gene NPHS1 lead to massive proteinuria. Along with the originally described location in the podocyte, nephrin has now been found to be expressed in the brain, testis, placenta and pancreatic beta cells. In type 1 diabetes, the fundamental pathologic event is the autoimmune destruction of the beta cells. Autoantibodies against various beta cell antigens are generated during this process. Due to the location of nephrin in the beta cell, we hypothesized that patients with type 1 diabetes may present with nephrin autoantibodies. We also wanted to test whether such autoantibodies could be involved in the pathogenesis of diabetic nephropathy. The puromycin aminonucleoside nephrosis model in the rat, the streptozotocin model in the rat, and the non-obese diabetic mice were studied by immunochemical techniques, in situ -hybridization and the polymerase chain reaction -based methods to resolve the expression of nephrin mRNA and protein in experimental nephropathies. To test the effect of antiproteinuric therapies, streptozotocin-treated rats were also treated with aminoguanidine or perindopril. To detect nephrin antibodies we developed a radioimmunoprecipitation assay and analyzed follow-up material of 66 patients with type 1 diabetes. In the puromycin aminonucleoside nephrosis model, the nephrin expression level was uniformly decreased together with the appearance of proteinuria. In the streptozotocin-treated rats and in non-obese diabetic mice, the nephrin mRNA and protein expression levels were seen to increase in the early stages of nephropathy. However, as observed in the streptozotocin rats, in prolonged diabetic nephropathy the expression level decreased. We also found out that treatment with perindopril could not only prevent proteinuria but also a decrease in nephrin expression in streptozotocin-treated rats. Aminoguanidine did not have an effect on nephrin expression, although it could attenuate the proteinuria. Circulating antibodies to nephrin in patients with type 1 diabetes were found, although there was no correlation with the development of diabetic nephropathy. At diagnosis, 24% of the patients had these antibodies, while at 2, 5 and 10 years of disease duration the respective proportions were 23%, 14% and 18%. During the total follow-up of 16 to 19 years after diagnosis of diabetes, 14 patients had signs of nephropathy and 29% of them tested positive for nephrin autoantibodies in at least one sample. In conclusion, this thesis work could show changes of nephrin expression along with the development of proteinuria. The autoantibodies against nephrin are likely generated in the autoimmune process leading to type 1 diabetes. However, according to the present work it is unlikely that these autoantibodies are contributing significantly to the development of diabetic nephropathy.

The ozone transfer between atmosphere and vegetation. A study on Scots pine in the field

Ozone (O3) is a reactive gas present in the troposphere in the range of parts per billion (ppb), i.e. molecules of O3 in 109 molecules of air. Its strong oxidative capacity makes it a key element in tropospheric chemistry and a threat to the integrity of materials, including living organisms. Knowledge and control of O3 levels are an issue in relation to indoor air quality, building material endurance, respiratory human disorders, and plant performance. Ozone is also a greenhouse gas and its abundance is relevant to global warming. The interaction of the lower troposphere with vegetated landscapes results in O3 being removed from the atmosphere by reactions that lead to the oxidation of plant-related components. Details on the rate and pattern of removal on different landscapes as well as the ultimate mechanisms by which this occurs are not fully resolved. This thesis analysed the controlling processes of the transfer of ozone at the air-plant interface. Improvement in the knowledge of these processes benefits the prediction of both atmospheric removal of O3 and its impact on vegetation. This study was based on the measurement and analysis of multi-year field measurements of O3 flux to Scots pine (Pinus sylvestris L.) foliage with a shoot-scale gas-exchange enclosure system. In addition, the analyses made use of simultaneous CO2 and H2O exchange, canopy-scale O3, CO2 and H2O exchange, foliage surface wetness, and environmental variables. All data was gathered at the SMEAR measuring station (southern Finland). Enclosure gas-exchange techniques such as those commonly used for the measure of CO2 and water vapour can be applied to the measure of ozone gas-exchange in the field. Through analysis of the system dynamics the occurring disturbances and noise can be identified. In the system used in this study, the possible artefacts arising from the ozone reactivity towards the system materials in combination with low background concentrations need to be taken into account. The main artefact was the loss of ozone towards the chamber walls, which was found to be very variable. The level of wall-loss was obtained from simultaneous and continuous measurements, and was included in the formulation of the mass balance of O3 concentration inside the chamber. The analysis of the field measurements in this study show that the flux of ozone to the Scots pine foliage is generated in about equal proportions by stomatal and non-stomatal controlled processes. Deposition towards foliage and forest is sustained also during night and winter when stomatal gas-exchange is low or absent. The non-stomatal portion of the flux was analysed further. The pattern of flux in time was found to be an overlap of the patterns of biological activity and presence of wetness in the environment. This was seen to occur both at the shoot and canopy scale. The presence of wetness enhanced the flux not only in the presence of liquid droplets but also during existence of a moisture film on the plant surfaces. The existence of these films and their relation to the ozone sinks was determined by simultaneous measurements of leaf surface wetness and ozone flux. The results seem to suggest ozone would be reacting at the foliage surface and the reaction rate would be mediated by the presence of surface wetness. Alternative mechanisms were discussed, including nocturnal stomatal aperture and emission of reactive volatile compounds. The prediction of the total flux could thus be based on a combination of a model of stomatal behaviour and a model of water absorption on the foliage surfaces. The concepts behind the division of stomatal and non-stomatal sinks were reconsidered. This study showed that it is theoretically possible that a sink located before or near the stomatal aperture prevents or diminishes the diffusion of ozone towards the intercellular air space of the mesophyll. This obstacle to stomatal diffusion happens only under certain conditions, which include a very low presence of reaction sites in the mesophyll, an extremely strong sink located on the outer surfaces or stomatal pore. The relevance, or existence, of this process in natural conditions would need to be assessed further. Potentially strong reactions were considered, including dissolved sulphate, volatile organic compounds, and apoplastic ascorbic acid. Information on the location and the relative abundance of these compounds would be valuable. The highest total flux towards the foliage and forest happens when both the plant activity and ambient moisture are high. The highest uptake into the interior of the foliage happens at large stomatal apertures, provided that scavenging reactions located near the stomatal pore are weak or non-existent. The discussion covers the methodological developments of this study, the relevance of the different controlling factors of ozone flux, the partition amongst its component, and the possible mechanisms of non-stomatal uptake.

Concentrate feeding strategies for growing and finishing dairy bulls offered grass silage-based diets

The first aim of this thesis was to produce data for evaluating, developing and recommending biologically and economically efficient energy and protein feeding strategies for growing and finishing dairy bulls offered grass silage-based diets. The second aim was to calculate the energy and protein supplies of the dairy bulls fed different grass silage-cereal-based diets and, based on this, to estimate the possible need to revise the current Finnish energy and protein recommendations for growing dairy bulls. The third aim was to demonstrate the phosphorus supply of dairy bulls fed grass silage-cereal-based diets with or without protein supplementation in relation to current feeding recommendations for phosphorus. The results indicate that protein supplement is not needed for finishing dairy bulls (live weight more than 250 kg) when they are fed good-quality grass silage (digestible organic matter more than 650 g/kg dry matter, restricted fermentation with low concentrations of fermentation acids and ammonia N) and grain-based concentrate with a moderate (300-700 g/kg dry matter) concentrate level. The results also suggest that with total mixed ration feeding it is possible to use rather high concentrate proportions (700 g/kg dry matter) in feeding dairy bulls. According to this study, barley fibre is a suitable energy supplement with good-quality silage for growing dairy bulls. The results suggest that 50% of barley grain can be replaced with barley fibre without affecting growth. Also oats is a suitable energy supplement for dairy bulls. However, as a consequence of decreased energy intake, the gain and feed conversion of the bulls were slightly reduced in this study when barley grain was replaced by oats in the diet. Ultimately, the rationality of the use of barley fibre and oats in the future will depend on the price in relation to other concentrates. During the feeding experiments the calculated supply of energy was 10% higher than in the Finnish feeding recommendations for the present growth rate. This indicates that there is a need to update the Finnish feeding recommendations for dairy-breed growing bulls, and further calculations are needed for the energy supply of growing dairy bulls. The calculated supply of AAT (amino acids absorbed from the small intestine) was 38% higher than in the Finnish feeding recommendations for the present growth. Possibly, the present AAT-PBV system is not an optimal protein evaluation system for growing dairy bulls more than 250 kg live weight. The calculations based on the feeding experiments and the Finnish feeding recommendations indicate that in most cases the dairy bulls (live weight more than 250 kg) received enough P from the basic grass silage cereal-based diets without additional mineral feeds. Therefore there is no need to add P in the form of mineral mixtures.

Lignocellulose degradation and humus modification by the fungus Paecilomyces inflatus

Composting is the biological conversion of solid organic waste into usable end products such as fertilizers, substrates for mushroom production and biogas. Although composts are highly variable in their bulk composition, composting material is generally based on lignocellulose compounds derived from agricultural, forestry, fruit and vegetable processing, household and municipal wastes. Lignocellulose is very recalcitrant; however it is rich and abundant source of carbon and energy. Therefore lignocellulose degradation is essential for maintaining the global carbon cycle. In compost, the active component involved in the biodegradation and conversion processes is the resident microbial population, among which microfungi play a very important role. In composting pile the warm, humid, and aerobic environment provides the optimal conditions for their development. Microfungi use many carbon sources, including lignocellulosic polymers and can survive in extreme conditions. Typically microfungi are responsible for compost maturation. In order to improve the composting process, more information is needed about the microbial degradation process. Better knowledge on the lignocellulose degradation by microfungi could be used to optimize the composting process. Thus, this thesis focused on lignocellulose and humic compounds degradation by a microfungus Paecilomyces inflatus, which belongs to a flora of common microbial compost, soil and decaying plant remains. It is a very common species in Europe, North America and Asia. The lignocellulose and humic compounds degradation was studied using several methods including measurements of carbon release from 14C-labelled compounds, such as synthetic lignin (dehydrogenative polymer, DHP) and humic acids, as well as by determination of fibre composition using chemical detergents and sulphuric acid. Spectrophotometric enzyme assays were conducted to detect extracellular lignocellulose-degrading hydrolytic and oxidative enzymes. Paecilomyces inflatus secreted clearly extracellular laccase to the culture media. Laccase was involved in the degradation process of lignin and humic acids. In compost P. inflatus mineralised 6-10% of 14C-labelled DHP into carbon dioxide. About 15% of labelled DHP was converted into water-soluble compounds. Also humic acids were partly mineralised and converted into water-soluble material, such as low-molecular mass fulvic acid-like compounds. Although laccase activity in aromatics-rich compost media clearly is connected with the degradation process of lignin and lignin-like compounds, it may preferentially effect the polymerisation and/or detoxification of such aromatic compounds. P. inflatus can degrade lignin and carbohydrates also while growing in straw and in wood. The cellulolytic enzyme system includes endoglucanase and β-glucosidase. In P. inflatus the secretion of these enzymes was stimulated by low-molecular-weight aromatics, such as soil humic acid and veratric acid. When strains of P. inflatus from different ecophysiological origins were compared, indications were found that specific adaptation strategies needed for lignocellulosics degradation may operate in P. inflatus. The degradative features of these microfungi are on relevance for lignocellulose decomposition in nature, especially in soil and compost environments, where basidiomycetes are not established. The results of this study may help to understand, control and better design the process of plant polymer conversion in compost environment, with a special emphasis on the role of ubiquitous microfungi.

Essays on socially optimal phosphorus policies in crop production

Phosphorus is a nutrient needed in crop production. While boosting crop yields it may also accelerate eutrophication in the surface waters receiving the phosphorus runoff. The privately optimal level of phosphorus use is determined by the input and output prices, and the crop response to phosphorus. Socially optimal use also takes into account the impact of phosphorus runoff on water quality. Increased eutrophication decreases the economic value of surface waters by Deteriorating fish stocks, curtailing the potential for recreational activities and by increasing the probabilities of mass algae blooms. In this dissertation, the optimal use of phosphorus is modelled as a dynamic optimization problem. The potentially plant available phosphorus accumulated in soil is treated as a dynamic state variable, the control variable being the annual phosphorus fertilization. For crop response to phosphorus, the state variable is more important than the annual fertilization. The level of this state variable is also a key determinant of the runoff of dissolved, reactive phosphorus. Also the loss of particulate phosphorus due to erosion is considered in the thesis, as well as its mitigation by constructing vegetative buffers. The dynamic model is applied for crop production on clay soils. At the steady state, the analysis focuses on the effects of prices, damage parameterization, discount rate and soil phosphorus carryover capacity on optimal steady state phosphorus use. The economic instruments needed to sustain the social optimum are also analyzed. According to the results the economic incentives should be conditioned on soil phosphorus values directly, rather than on annual phosphorus applications. The results also emphasize the substantial effects the differences in varying discount rates of the farmer and the social planner have on optimal instruments. The thesis analyzes the optimal soil phosphorus paths from its alternative initial levels. It also examines how erosion susceptibility of a parcel affects these optimal paths. The results underline the significance of the prevailing soil phosphorus status on optimal fertilization levels. With very high initial soil phosphorus levels, both the privately and socially optimal phosphorus application levels are close to zero as the state variable is driven towards its steady state. The soil phosphorus processes are slow. Therefore, depleting high phosphorus soils may take decades. The thesis also presents a methodologically interesting phenomenon in problems of maximizing the flow of discounted payoffs. When both the benefits and damages are related to the same state variable, the steady state solution may have an interesting property, under very general conditions: The tail of the payoffs of the privately optimal path as well as the steady state may provide a higher social welfare than the respective tail of the socially optimal path. The result is formalized and an applied to the created framework of optimal phosphorus use.

Multiangular Spectrometry and Optical Properties of Debris Covered Surfaces

Due to the recent development in CCD technology aerial photography is now slowly changing from film to digital cameras. This new aspect in remote sensing allows and requires also new automated analysis methods. Basic research on reflectance properties of natural targets is needed so that computerized processes could be fully utilized. For this reason an instrument was developed at Finnish Geodetic Institute for measurement of multiangular reflectance of small remote sensing targets e.g. forest understorey or asphalt. Finnish Geodetic Institute Field Goniospectrometer (FiGIFiGo) is a portable device that is operated by 1 or 2 persons. It can be reassembled to a new location in 15 minutes and after that a target's multiangular reflectance can be measured in 10 - 30 minutes (with one illumination angle). FiGIFiGo has effective spectral range approximately from 400 nm to 2000 nm. The measurements can be made either outside with sunlight or in laboratory with 1000 W QTH light source. In this thesis FiGIFiGo is introduced and the theoretical basis of such reflectance measurements are discussed. A new method is introduced for extraction of subcomponent proportions from reflectance of a mixture sample, e.g. for retrieving proportion of lingonberry's reflectance in observation of lingonberry-lichen sample. This method was tested by conducting a series of measurements on reflectance properties of artificial samples. The component separation method yielded sound results and brought up interesting aspects in targets' reflectances. The method and the results still need to be verified with further studies, but the preliminary results imply that this method could be a valuable tool in analysis of such mixture samples.

Remote sensing of boreal land cover: estimation of forest attributes and extent

Remote sensing provides methods to infer land cover information over large geographical areas at a variety of spatial and temporal resolutions. Land cover is input data for a range of environmental models and information on land cover dynamics is required for monitoring the implications of global change. Such data are also essential in support of environmental management and policymaking. Boreal forests are a key component of the global climate and a major sink of carbon. The northern latitudes are expected to experience a disproportionate and rapid warming, which can have a major impact on vegetation at forest limits. This thesis examines the use of optical remote sensing for estimating aboveground biomass, leaf area index (LAI), tree cover and tree height in the boreal forests and tundra taiga transition zone in Finland. The continuous fields of forest attributes are required, for example, to improve the mapping of forest extent. The thesis focus on studying the feasibility of satellite data at multiple spatial resolutions, assessing the potential of multispectral, -angular and -temporal information, and provides regional evaluation for global land cover data. Preprocessed ASTER, MISR and MODIS products are the principal satellite data. The reference data consist of field measurements, forest inventory data and fine resolution land cover maps. Fine resolution studies demonstrate how statistical relationships between biomass and satellite data are relatively strong in single species and low biomass mountain birch forests in comparison to higher biomass coniferous stands. The combination of forest stand data and fine resolution ASTER images provides a method for biomass estimation using medium resolution MODIS data. The multiangular data improve the accuracy of land cover mapping in the sparsely forested tundra taiga transition zone, particularly in mires. Similarly, multitemporal data improve the accuracy of coarse resolution tree cover estimates in comparison to single date data. Furthermore, the peak of the growing season is not necessarily the optimal time for land cover mapping in the northern boreal regions. The evaluated coarse resolution land cover data sets have considerable shortcomings in northernmost Finland and should be used with caution in similar regions. The quantitative reference data and upscaling methods for integrating multiresolution data are required for calibration of statistical models and evaluation of land cover data sets. The preprocessed image products have potential for wider use as they can considerably reduce the time and effort used for data processing.

A Probabilistic Approach to the Primary Visual Cortex

What can the statistical structure of natural images teach us about the human brain? Even though the visual cortex is one of the most studied parts of the brain, surprisingly little is known about how exactly images are processed to leave us with a coherent percept of the world around us, so we can recognize a friend or drive on a crowded street without any effort. By constructing probabilistic models of natural images, the goal of this thesis is to understand the structure of the stimulus that is the raison d etre for the visual system. Following the hypothesis that the optimal processing has to be matched to the structure of that stimulus, we attempt to derive computational principles, features that the visual system should compute, and properties that cells in the visual system should have. Starting from machine learning techniques such as principal component analysis and independent component analysis we construct a variety of sta- tistical models to discover structure in natural images that can be linked to receptive field properties of neurons in primary visual cortex such as simple and complex cells. We show that by representing images with phase invariant, complex cell-like units, a better statistical description of the vi- sual environment is obtained than with linear simple cell units, and that complex cell pooling can be learned by estimating both layers of a two-layer model of natural images. We investigate how a simplified model of the processing in the retina, where adaptation and contrast normalization take place, is connected to the nat- ural stimulus statistics. Analyzing the effect that retinal gain control has on later cortical processing, we propose a novel method to perform gain control in a data-driven way. Finally we show how models like those pre- sented here can be extended to capture whole visual scenes rather than just small image patches. By using a Markov random field approach we can model images of arbitrary size, while still being able to estimate the model parameters from the data.

Computationally Efficient Methods for MDL-Optimal Density Estimation and Data Clustering

The Minimum Description Length (MDL) principle is a general, well-founded theoretical formalization of statistical modeling. The most important notion of MDL is the stochastic complexity, which can be interpreted as the shortest description length of a given sample of data relative to a model class. The exact definition of the stochastic complexity has gone through several evolutionary steps. The latest instantation is based on the so-called Normalized Maximum Likelihood (NML) distribution which has been shown to possess several important theoretical properties. However, the applications of this modern version of the MDL have been quite rare because of computational complexity problems, i.e., for discrete data, the definition of NML involves an exponential sum, and in the case of continuous data, a multi-dimensional integral usually infeasible to evaluate or even approximate accurately. In this doctoral dissertation, we present mathematical techniques for computing NML efficiently for some model families involving discrete data. We also show how these techniques can be used to apply MDL in two practical applications: histogram density estimation and clustering of multi-dimensional data.

Computational Methods for Locating and Analyzing Conserved Gene Regulatory DNA Elements

This thesis presents methods for locating and analyzing cis-regulatory DNA elements involved with the regulation of gene expression in multicellular organisms. The regulation of gene expression is carried out by the combined effort of several transcription factor proteins collectively binding the DNA on the cis-regulatory elements. Only sparse knowledge of the 'genetic code' of these elements exists today. An automatic tool for discovery of putative cis-regulatory elements could help their experimental analysis, which would result in a more detailed view of the cis-regulatory element structure and function. We have developed a computational model for the evolutionary conservation of cis-regulatory elements. The elements are modeled as evolutionarily conserved clusters of sequence-specific transcription factor binding sites. We give an efficient dynamic programming algorithm that locates the putative cis-regulatory elements and scores them according to the conservation model. A notable proportion of the high-scoring DNA sequences show transcriptional enhancer activity in transgenic mouse embryos. The conservation model includes four parameters whose optimal values are estimated with simulated annealing. With good parameter values the model discriminates well between the DNA sequences with evolutionarily conserved cis-regulatory elements and the DNA sequences that have evolved neutrally. In further inquiry, the set of highest scoring putative cis-regulatory elements were found to be sensitive to small variations in the parameter values. The statistical significance of the putative cis-regulatory elements is estimated with the Two Component Extreme Value Distribution. The p-values grade the conservation of the cis-regulatory elements above the neutral expectation. The parameter values for the distribution are estimated by simulating the neutral DNA evolution. The conservation of the transcription factor binding sites can be used in the upstream analysis of regulatory interactions. This approach may provide mechanistic insight to the transcription level data from, e.g., microarray experiments. Here we give a method to predict shared transcriptional regulators for a set of co-expressed genes. The EEL (Enhancer Element Locator) software implements the method for locating putative cis-regulatory elements. The software facilitates both interactive use and distributed batch processing. We have used it to analyze the non-coding regions around all human genes with respect to the orthologous regions in various other species including mouse. The data from these genome-wide analyzes is stored in a relational database which is used in the publicly available web services for upstream analysis and visualization of the putative cis-regulatory elements in the human genome.

Carbon balance and component CO2 fluxes in boreal Scots pine stands

The forest vegetation takes up atmospheric carbon dioxide (CO2) in photosynthesis. Part of the fixed carbon is released back into the atmosphere during plant respiration but a substantial part is stored as plant biomass, especially in the stems of trees. Carbon also accumulates in the soil as litter and via the roots. CO2 is released into the atmosphere from these carbon stocks in the decomposition of dead biomass. Carbon balance of a forest stand is the difference between the CO2 uptake and CO2 efflux. This study quantifies and analyses the dynamics of carbon balance and component CO2 fluxes in four Southern Finnish Scots pine stands that covered the typical economic rotation time of 80 years. The study was based on direct flux measurements with chambers and eddy covariance (EC), and modelling of component CO2 fluxes. The net CO2 exchange of the stand was partitioned into component fluxes: photosynthesis of trees and ground vegetation, respiration of tree foliage and stems, and CO2 efflux from the soil. The relationships between the component fluxes and the environmental factors (light, temperature, atmospheric CO2, air humidity and soil moisture) were studied with mathematical modelling. The annual CO2 balance varied from a source of about 400 g C/m2 at a recently clearcut site to net CO2 uptake of 200 300 g C/m2 in a middle-aged (40-year-old) and a mature (75-year-old) stand. A 12-year-old sapling site was at the turning point from source to a sink of CO2. In the middle-aged stand, photosynthetic production was dominated by trees. Under closed pine canopies, ground vegetation accounted for 10 20% of stand photosynthesis whereas at the open sites the proportion and also the absolute photosynthesis of ground vegetation was much higher. The aboveground respiration was dominated by tree foliage which accounted for one third of the ecosystem respiration. Rate of wood respiration was in the order of 10% of total ecosystem respiration. CO2 efflux from the soil dominated the ecosystem respiratory fluxes in all phases of stand development. Instantaneous and delayed responses to the environmental driving factors could predict well within-year variability in photosynthetic production: In the short term and during the growing season photosynthesis follows primarily light while the seasonal variation is more strongly connected to temperature. The temperature relationship of the annual cycle of photosynthesis was found to be almost equal in the southern boreal zone and at the timberline in the northern boreal zone. The respiratory fluxes showed instantaneous and seasonal temperature relationships but they could also be connected to photosynthesis at an annual timescale.

Barrett's Esophagus and Associated Adenocarcinoma : studies on pathogenesis, clinical staging and cost-utility of cancer treatment

Esophageal and gastroesophageal junction (GEJ) adenocarcinoma is rapidly increasing disease with a pathophysiology connected to oxidative stress. Exact pre-treatment clinical staging is essential for optimal care of this lethal malignancy. The cost-effectiviness of treatment is increasingly important. We measured oxidative metabolism in the distal and proximal esophagus by myeloperoxidase activity (MPA), glutathione content (GSH), and superoxide dismutase (SOD) in 20 patients operated on with Nissen fundoplication and 9 controls during a 4-year follow-up. Further, we assessed the oxidative damage of DNA by 8-hydroxydeoxyguanosine (8-OHdG) in esophageal samples of subjects (13 Barrett s metaplasia, 6 Barrett s esophagus with high-grade dysplasia, 18 adenocarcinoma of the distal esophagus/GEJ, and 14 normal controls). We estimated the accuracy (42 patients) and preoperative prognostic value (55 patients) of PET compared with computed tomography (CT) and endoscopic ultrasound (EUS) in patients with adenocarcinoma of the esophagus/GEJ. Finally, we clarified the specialty-related costs and the utility of either radical (30 patients) or palliative (23 patients) treatment of esophageal/GEJ carcinoma by the 15 D health-related quality-of-life (HRQoL) questionnaire and the survival rate. The cost-utility of radical treatment of esophageal/GEJ carcinoma was investigated using a decision tree analysis model comparing radical, palliative, and hypothetical new treatment. We found elevated oxidative stress ( measured by MPA) and decreased antioxidant defense (measured by GSH) after antireflux surgery. This indicates that antireflux surgery is not a perfect solution for oxidative stress of the esophageal mucosa. Elevated oxidative stress in turn may partly explain why adenocarcinoma of the distal esophagus is found even after successful fundoplication. In GERD patients, proximal esophageal mucosal anti-oxidative defense seems to be defective before and even years after successful antireflux surgery. In addition, antireflux surgery apparently does not change the level of oxidative stress in the proximal esophagus, suggesting that defective mucosal anti-oxidative capacity plays a role in development of oxidative damage to the esophageal mucosa in GERD. In the malignant transformation of Barrett s esophagus an important component appears to be oxidative stress. DNA damage may be mediated by 8-OHdG, which we found to be increased in Barrett s epithelium and in high-grade dysplasia as well as in adenocarcinoma of the esophagus/GEJ compared with controls. The entire esophagus of Barrett s patients suffers from increased oxidative stress ( measured by 8-OhdG). PET is a useful tool in the staging and prognostication of adenocarcinoma of the esophagus/GEJ detecting organ metastases better than CT, although its accuracy in staging of paratumoral and distant lymph nodes is limited. Radical surgery for esophageal/GEJ carcinoma provides the greatest benefit in terms of survival, and its cost-utility appears to be the best of currently available treatments.