Algorithmic Techniques in Gene Expression Processing. From Imputation to Visualization

The amount of biological data has grown exponentially in recent decades. Modern biotechnologies, such as microarrays and next-generation sequencing, are capable to produce massive amounts of biomedical data in a single experiment. As the amount of the data is rapidly growing there is an urgent need for reliable computational methods for analyzing and visualizing it. This thesis addresses this need by studying how to efficiently and reliably analyze and visualize high-dimensional data, especially that obtained from gene expression microarray experiments. First, we will study the ways to improve the quality of microarray data by replacing (imputing) the missing data entries with the estimated values for these entries. Missing value imputation is a method which is commonly used to make the original incomplete data complete, thus making it easier to be analyzed with statistical and computational methods. Our novel approach was to use curated external biological information as a guide for the missing value imputation. Secondly, we studied the effect of missing value imputation on the downstream data analysis methods like clustering. We compared multiple recent imputation algorithms against 8 publicly available microarray data sets. It was observed that the missing value imputation indeed is a rational way to improve the quality of biological data. The research revealed differences between the clustering results obtained with different imputation methods. On most data sets, the simple and fast k-NN imputation was good enough, but there were also needs for more advanced imputation methods, such as Bayesian Principal Component Algorithm (BPCA). Finally, we studied the visualization of biological network data. Biological interaction networks are examples of the outcome of multiple biological experiments such as using the gene microarray techniques. Such networks are typically very large and highly connected, thus there is a need for fast algorithms for producing visually pleasant layouts. A computationally efficient way to produce layouts of large biological interaction networks was developed. The algorithm uses multilevel optimization within the regular force directed graph layout algorithm.

The functional and biological properties of whey proteins : prospects for the development of functional foods

Selostus: Heraproteiinit terveysvaikutteisten elintarvikkeiden kehittämisessä

Biological contamination of barley and carrots by pathogens in soil fertilized with anthropogenic nutrients

Conference publication in Kalmar ECO-TECH’07. International Conference on Technologies for Waste and Wastewater Treatment, Energy from Waste, Remediation of Contaminated Sites and Emissions Related to Climate November 26-28 2007, Kalmar, Sweden.

Food Safety Testing: Rapid Molecular Methods for Chemical and Biological Hazards

The central goal of food safety policy in the European Union (EU) is to protect consumer health by guaranteeing a high level of food safety throughout the food chain. This goal can in part be achieved by testing foodstuffs for the presence of various chemical and biological hazards. The aim of this study was to facilitate food safety testing by providing rapid and user-friendly methods for the detection of particular food-related hazards. Heterogeneous competitive time-resolved fluoroimmunoassays were developed for the detection of selected veterinary residues, that is coccidiostat residues, in eggs and chicken liver. After a simplified sample preparation procedure, the immunoassays were performed either in manual format with dissociation-enhanced measurement or in automated format with pre-dried assay reagents and surface measurement. Although the assays were primarily designed for screening purposes providing only qualitative results, they could also be used in a quantitative mode. All the developed assays had good performance characteristics enabling reliable screening of samples at concentration levels required by the authorities. A novel polymerase chain reaction (PCR)-based assay system was developed for the detection of Salmonella spp. in food. The sample preparation included a short non-selective pre-enrichment step, after which the target cells were collected with immunomagnetic beads and applied to PCR reaction vessels containing all the reagents required for the assay in dry form. The homogeneous PCR assay was performed with a novel instrument platform, GenomEra^™, and the qualitative assay results were automatically interpreted based on end-point time-resolved fluorescence measurements and cut-off values. The assay was validated using various food matrices spiked with sub-lethally injured Salmonella cells at levels of 1-10 colony forming units (CFU)/25 g of food. The main advantage of the system was the exceptionally short time to result; the entire process starting from the pre-enrichment and ending with the PCR result could be completed in eight hours. In conclusion, molecular methods using state-of-the-art assay techniques were developed for food safety testing. The combination of time-resolved fluorescence detection and ready-to-use reagents enabled sensitive assays easily amenable to automation. Consequently, together with the simplified sample preparation, these methods could prove to be applicable in routine testing.

Osteoclastic tartrate-resistant acid phosphatase 5b. Diagnostic use and biological significance in bone physiology

The skeleton undergoes continuous turnover throughout life. In women, an increase in bone turnover is pronounced during childhood and puberty and after menopause. Bone turnover can be monitored by measuring biochemical markers of bone resorption and bone formation. Tartrate-resistant acid phosphatase (TRACP) is an enzyme secreted by osteoclasts, macrophages and dendritic cells. The secreted enzyme can be detected from the blood circulation by recently developed immunoassays. In blood circulation, the enzyme exists as two isoforms, TRACP 5a with an intact polypeptide chain and TRACP 5b in which the polypeptide chain consists of two subunits. The 5b form is predominantly secreted by osteoclasts and is thus associated with bone turnover. The secretion of TRACP 5b is not directly related to bone resorption; instead, the levels are shown to be proportional to the number of osteoclasts. Therefore, the combination of TRACP 5b and a marker reflecting bone degradation, such as C-terminal cross-linked telopeptides of type I collagen (CTX), enables a more profound analysis of the changes in bone turnover. In this study, recombinant TRACP 5a-like protein was proteolytically processed into TRACP 5b-like two subunit form. The 5b-like form was more active both as an acid phosphatase and in producing reactive oxygen species, suggesting a possible function for TRACP 5b in osteoclastic bone resorption. Even though both TRACP 5a and 5b were detected in osteoclasts, serum TRACP 5a levels demonstrated no change in response to alendronate treatment of postmenopausal women. However, TRACP 5b levels decreased substantially, demonstrating that alendronate decreases the number of osteoclasts. This was confirmed in human osteoclast cultures, showing that alendronate decreased the number of osteoclats by inducing osteoclast apoptosis, and TRACP 5b was not secreted as an active enzyme from the apoptotic osteoclasts. In peripubertal girls, the highest levels of TRACP 5b and other bone turnover markers were observed at the time of menarche, whereas at the same time the ratio of CTX to TRACP 5b was lowest, indicating the presence of a high number of osteoclasts with decreased resorptive activity. These results support the earlier findings that TRACP 5b is the predominant form of TRACP secreted by osteoclasts. The major source of circulating TRACP 5a remains to be established, but is most likely other cells of the macrophage-monocyte system. The results also suggest that bone turnover can be differentially affected by both osteoclast number and their resorptive activity, and provide further support for the possible clinical use of TRACP 5b as a marker of osteoclast number.

The oxygenation in biological wastewater treatment of pulp and paper industry

The UPM-Kymmene Oyj Pietarsaari pulp and paper Mill biological wastewater treatment plant was built in the 1980's and the plant has been in use ever since. During the past years there have been problems with deviations. The wastewater treatment plant needs update, especially the aeration basin, where the old surface aerators cannot produce enough mixing and indroduce oxygen enough to the wastewater. In this thesis how extra aeration with oxygen affects the wastewater treatment plant effluent was studied. In the literature part the main focus is in aeration devices, which can be used in biological wastewater treatment. The target is to compare different kind of aerators, which are suitable for pulp and paper wastewater treatment. Studies show, that EDI-aerators are commonly used and also most suitable. In the experimental part, the focus is on the Pietarsaari Mills wastewater treatment plant and oxygen aeration during autumn 2008. This thesis presents the results of the trial run. Studies show, that extra oxygen devices can produce lot a of mixing and the oxygenation capacity was more than what the micro-organisms needed. The effect on sludge quality could not been seen during the trial runs.

Visualization of communication system testing logs

During the past decades testing has matured from ad-hoc activity into being an integral part of the development process. The benefits of testing are obvious for modern communication systems, which operate in heterogeneous environments amongst devices from various manufacturers. The increased demand for testing also creates demand for tools and technologies that support and automate testing activities. This thesis discusses applicability of visualization techniques in the result analysis part of the testing process. Particularly, the primary focus of this work is visualization of test execution logs produced by a TTCN-3 test system. TTCN-3 is an internationally standardized test specification and implementation language. The TTCN-3 standard suite includes specification of a test logging interface and a graphical presentation format, but no immediate relationship between them. This thesis presents a technique for mapping the log events to the graphical presentation format along with a concrete implementation, which is integrated with the Eclipse Platform and the OpenTTCN Tester toolchain. Results of this work indicate that for majority of the log events, a visual representation may be derived from the TTCN-3 standard suite. The remaining events were analysed and three categories relevant in either log analysis or implementation of the visualization tool were identified: events indicating insertion of something into the incoming queue of a port, events indicating a mismatch and events describing the control flow during the execution. Applicability of the results is limited into the domain of TTCN-3, but the developed mapping and the implementation may be utilized with any TTCN-3 tool that is able to produce the execution log in the standardized XML format.

Evaluating multidimensional visualization techniques in data mining tasks

Visual data mining (VDM) tools employ information visualization techniques in order to represent large amounts of high-dimensional data graphically and to involve the user in exploring data at different levels of detail. The users are looking for outliers, patterns and models – in the form of clusters, classes, trends, and relationships – in different categories of data, i.e., financial, business information, etc. The focus of this thesis is the evaluation of multidimensional visualization techniques, especially from the business user’s perspective. We address three research problems. The first problem is the evaluation of projection-based visualizations with respect to their effectiveness in preserving the original distances between data points and the clustering structure of the data. In this respect, we propose the use of existing clustering validity measures. We illustrate their usefulness in evaluating five visualization techniques: Principal Components Analysis (PCA), Sammon’s Mapping, Self-Organizing Map (SOM), Radial Coordinate Visualization and Star Coordinates. The second problem is concerned with evaluating different visualization techniques as to their effectiveness in visual data mining of business data. For this purpose, we propose an inquiry evaluation technique and conduct the evaluation of nine visualization techniques. The visualizations under evaluation are Multiple Line Graphs, Permutation Matrix, Survey Plot, Scatter Plot Matrix, Parallel Coordinates, Treemap, PCA, Sammon’s Mapping and the SOM. The third problem is the evaluation of quality of use of VDM tools. We provide a conceptual framework for evaluating the quality of use of VDM tools and apply it to the evaluation of the SOM. In the evaluation, we use an inquiry technique for which we developed a questionnaire based on the proposed framework. The contributions of the thesis consist of three new evaluation techniques and the results obtained by applying these evaluation techniques. The thesis provides a systematic approach to evaluation of various visualization techniques. In this respect, first, we performed and described the evaluations in a systematic way, highlighting the evaluation activities, and their inputs and outputs. Secondly, we integrated the evaluation studies in the broad framework of usability evaluation. The results of the evaluations are intended to help developers and researchers of visualization systems to select appropriate visualization techniques in specific situations. The results of the evaluations also contribute to the understanding of the strengths and limitations of the visualization techniques evaluated and further to the improvement of these techniques.

Disturbance and ecosystem functioning : the role of changes in benthic biological traits

Rapid changes in biodiversity are occurring globally, as a consequence of anthropogenic disturbance. This has raised concerns, since biodiversity is known to significantly contribute to ecosystem functions and services. Marine benthic communities participate in numerous functions provided by soft-sedimentary ecosystems. Eutrophication-induced oxygen deficiency is a growing threat against infaunal communities, both in open sea areas and in coastal zones. There is thus a need to understand how such disturbance affects benthic communities, and what is lost in terms of ecosystem functioning if benthic communities are harmed. In this thesis, the status of benthic biodiversity was assessed for the open Baltic Sea, a system severely affected by broad-scale hypoxia. Long-term monitoring data made it possible to establish quantitative biodiversity baselines against which change could be compared. The findings show that benthic biodiversity is currently severely impaired in large areas of the open Baltic Sea, from the Bornholm Basin to the Gulf of Finland. The observed reduction in biodiversity indicates that benthic communities are structurally and functionally impoverished in several of the sub-basins due to the hypoxic stress. A more detailed examination of disturbance impacts (through field studies and -experiments) on benthic communities in coastal areas showed that changes in benthic community structure and function took place well before species were lost from the system. The degradation of benthic community structure and function was directed by the type of disturbance, and its specific temporal and spatial characteristics. The observed shifts in benthic trait composition were primarily the result of reductions in species’ abundances, or of changes in demographic characteristics, such as the loss of large, adult bivalves. Reduction in community functions was expressed as declines in the benthic bioturbation potential and in secondary biomass production. The benthic communities and their degradation accounted for a substantial proportion of the changes observed in ecosystem multifunctionality. Individual ecosystem functions (i.e. measures of sediment ecosystem metabolism, elemental cycling, biomass production, organic matter transformation and physical structuring) were observed to differ in their response to increasing hypoxic disturbance. Interestingly, the results suggested that an impairment of ecosystem functioning could be detected at an earlier stage if multiple functions were considered. Importantly, the findings indicate that even small-scale hypoxic disturbance can reduce the buffering capacity of sedimentary ecosystem, and increase the susceptibility of the system towards further stress. Although the results of the individual papers are context-dependent, their combined outcome implies that healthy benthic communities are important for sustaining overall ecosystem functioning as well as ecosystem resilience in the Baltic Sea.

Interpreting marine benthic ecosystem functioning in coastal waters: validating the biological trait concept

Coastal areas harbour high biodiversity, but are simultaneously affected by rapid degradations of species and habitats due to human interactions. Such alterations also affect the functioning of the ecosystem, which is primarily governed by the characteristics or traits expressed by the organisms present. Marine benthic fauna is nvolved in numerous functions such as organic matter transformation and transport, secondary production, oxygen transport as well as nutrient cycling. Approaches utilising the variety of faunal traits to assess benthic community functioning have rapidly increased and shown the need for further development of the concept. In this thesis, I applied biological trait analysis that allows for assessments of a multitude of categorical traits and thus evaluation of multiple functional aspects simultaneously. I determined the functional trait structure, diversity and variability of coastal zoobenthic communities in the Baltic Sea. The measures were related to recruitment processes, habitat heterogeneity, large-scale environmental and taxonomic gradients as well as anthropogenic impacts. The studies comprised spatial scales from metres to thousands of kilometres, and temporal scales spanning one season as well as a decade. The benthic functional structure was found to vary within and between seagrass landscape microhabitats and four different habitats within a coastal bay, in papers I and II respectively. Expressions of trait categories varied within habitats, while the density of individuals was found to drive the functional differences between habitats. The findings in paper III unveiled high trait richness of Finnish coastal benthos (25 traits and 102 cateogries) although this differed between areas high and low in salinity and human pressure. In paper IV, the natural reduction in taxonomic richness across the Baltic Sea led to an overall reduction in function. However, functional richness in terms of number of trait categories remained comparatively high at low taxon richness. Changes in number of taxa within trait categories were also subtle and some individual categories were maintained or even increased. The temporal analysis in papers I and III highlighted generalities in trait expressions and dominant trait categories in a seagrass landscape as well as a “type organism” for the northern Baltic Sea. Some initial findings were made in all four papers on the role of common and rare species and traits for benthic community functioning. The findings show that common and rare species may not always express the same trait categories in relation to each other. Rare species in general did not express unique functional properties. In order to advance the understanding of the approach, I also assessed some issues concerning the limitations of the concept. This was conducted by evaluating the link between trait category and taxonomic richness using especially univariate measures. My results also show the need to collaborate nationally and internationally on safeguarding the utility of taxonomic and trait data. The findings also highlight the importance of including functional trait information into current efforts in marine spatial planning and biomonitoring.

Physical Attractiveness as a Signal of Biological Quality

It has been commonly thought that standards of beauty are arbitrary cultural conventions that vary between cultures and time. In my thesis I found that it is not so. Instead, I show that attractiveness and preferred traits serve as cues to phenotypic qualities that provide selective benefits for those who choose their mates based on these criteria. In the first study I show that attractive men have a stronger antibody response to the hepatitis b vaccine and higher levels of testosterone than their less attractive peers. Men with low levels of testosterone also tend to have high levels of the stress hormone cortisol, suggesting that their immune responses may have been inhibited by stress hormones. Thus, facial attractiveness may serve as an honest cue of the strength of immune defence in men. In the second study, I show that the attractiveness of the male body is also a cue of better immunity. In addition, I show that adiposity, both in men’s faces and bodies, is a better cue of the strength of immunity and attractiveness than of masculinity. In the third study, I test the preferences of women from 13 countries for facial cues of testosterone and cortisol. I show that there is cross-cultural variation in women’s preference for cues of testosterone and cortisol in male faces. I found a relationship between the health of a nation and women’s preferences for cues of testosterone in the male face and the interaction between preferences for cues of testosterone and cortisol. I show also a relationship between preferences for cues of testosterone and a societal-level measure of parasite stress. Thus, it seems that societal-level ecological factors influence the relative value of traits as revealed by combinations of testosterone and stress hormones. In the fourth study, I show that women’s immune responsiveness (amount of antibodies produced) does not predict facial attractiveness. Instead, plasma cortisol level is negatively associated with attractiveness, indicating that stressed women look less attractive. Fat percentage is curvilinearly associated with facial attractiveness, indicating that being too thin or too fat reduces attractiveness. This study suggests that in contrast to men, facial attractiveness in women does not indicate the strength of immune defence, but is associated with other aspects of long-term health and fertility: circulating levels of the stress hormone cortisol and the percentage of body fat. In the last study I show that the attractiveness of men’s body odor is positively correlated with stress hormone levels, suggesting also that the attractiveness of body odors may signal the phenotypic quality of males to females. However, the attractiveness of men’s body odor was not associated with testosterone levels. My thesis suggests that the standard of beauty is not in the eye of the beholder. Instead, our standard of beauty is hardwired in our brains by genes that are selected by natural selection and also influenced by current environmental conditions.