Glycogen debranching enzyme activity in the muscles of meat producing animals

Muscle glycogen exists in two forms: low molecular weight pro-glycogen and high molecular weight macro-glycogen. The degradation of glycogen to glucose 1 phosphate and free glucose is catalysed by glycogen phosphorylase together with glycogen debranching enzyme (GDE). The process in which glycogen is broken down via anaerobic pathways to lactate, results in the acidification of the muscles and has a great influence on meat quality. Thus, the overall aim of this thesis was to characterise the post mortem action of GDE in muscles of meat production animals (pigs, cattle and chickens). Interest was focused on the differences in GDE activity between fast twitch glycolytic muscles and slow twitch oxidative muscles. The effects of pH, temperature, RN genotype (PRKAG3 gene), and of time post mortem on GDE activity were also investigated. This thesis showed that there are differences in GDE activity between animal species and between different muscles of an animal. It was shown that in pigs and cattle, higher GDE activity and phosphorylase activity exists in the fast twitch glycolytic muscles than in slow twitch oxidative muscles of the same animal. Thus, the high activity of these enzymes enables a faster rate of glycogenolysis in glycolytic M. longissimus dorsi compared to oxidative M. masseter. In chicken muscles, the GDE activity was low compared to pig or cattle muscles. Furthermore, the GDE activity in the glycolytic M. pectoralis superficialis was lower than in more oxidative M. quadriceps femoris despite the high phosphorylase activity in the former. The relative ratios between phosphorylase and GDE activity were higher in fast twitch glycolytic muscles than in slow twitch oxidative muscles of all studied animals. This suggests that the relatively low GDE activity compared to the phosphorylase activity in fast twitch glycolytic muscles may be a protection mechanism in living muscle against a very fast pH decrease. Chilling significantly decreased GDE activity and below 15 C porcine GDE was almost inactive. The effect of pH on GDE activity was only minor at the range normally found in post mortem muscles (pH 7.4 to 5.0). The GDE activity remained level for several hours after slaughter. During the first hours post mortem, GDE activity was similar in RN- carrier pigs and in wild type pigs. However, the GDE activity declined faster in M. longissimus dorsi from wild type pigs than in the RN carrier pigs, the difference between genotypes was significant after 24 h post mortem. Pro-glycogen and macro-glycogen contents were higher, pH decrease was faster and ultimate pH was lower in RN- carrier pigs than in wild type pigs. In the RN- carriers, the prolonged high GDE activity level may enable an extended pH decrease and lower ultimate pH in their muscles. In conclusion, GDE is not the main factor determining the rate or the extent of post mortem glycogenolysis, but under certain conditions, such as in very fast chilling, the inhibition of GDE activity in meat may reduce the rate of pH decrease and result in higher ultimate pH. The rate and extent of pH decrease affects several meat quality traits.

Bacillus cereus spores and cereulide in food-borne illness

B. cereus is a gram-positive bacterium that possesses two different forms of life:the large, rod-shaped cells (ca. 0.002 mm by 0.004 mm) that are able to propagate and the small (0.001 mm), oval shaped spores. The spores can survive in almost any environment for up to centuries without nourishment or water. They are insensitive towards most agents that normally kill bacteria: heating up to several hours at 90 ºC, radiation, disinfectants and extreme alkaline (≥ pH 13) and acid (≤ pH 1) environment. The spores are highly hydrophobic and therefore make them tend to stick to all kinds of surfaces, steel, plastics and live cells. In favorable conditions the spores of B. cereus may germinate into vegetative cells capable of producing food poisoning toxins. The toxins can be heat-labile protein formed after ingestion of the contaminated food, inside the gastrointestinal tract (diarrhoeal toxins), or heat stable peptides formed in the food (emesis causing toxin, cereulide). Cereulide cannot be inactivated in foods by cooking or any other procedure applicable on food. Cereulide in consumed food causes serious illness in human, even fatalities. In this thesis, B. cereus strains originating from different kinds of foods and environments and 8 different countries were inspected for their capability of forming cereulide. Of the 1041 isolates from soil, animal feed, water, air, used bedding, grass, dung and equipment only 1.2 % were capable of producing cereulide, whereas of the 144 isolates originating from foods 24 % were cereulide producers. Cereulide was detected by two methods: by its toxicity towards mammalian cells (sperm assay) and by its peculiar chemical structure using liquid-chromatograph-mass spectrometry equipment. B. cereus is known as one of the most frequent bacteria occurring in food. Most foods contain more than one kind of B. cereus. When randomly selected 100 isolates of B. cereus from commercial infant foods (dry formulas) were tested, 11% of these produced cereulide. Considering a frequent content of 103 to 104 cfu (colony forming units) of B. cereus per gram of infant food formula (dry), it appears likely that most servings (200 ml, 30 g of the powder reconstituted with water) may contain cereulide producers. When a reconstituted infant formula was inoculated with >105 cfu of cereulide producing B. cereus per ml and left at room temperature, cereulide accumulated to food poisoning levels (> 0.1 mg of cereulide per serving) within 24 hours. Paradoxically, the amount of cereulide (per g of food) increased 10 to 50 fold when the food was diluted 4 - 15 fold with water. The amount of the produced cereulide strongly depended on the composition of the formula: most toxin was formed in formulas with cereals mixed with milk, and least toxin in formulas based on milk only. In spite of the aggressive cleaning practices executed by the modern dairy industry, certain genotypes of B. cereus appear to colonise the silos tanks. In this thesis four strategies to explain their survival of their spores in dairy silos were identified. First, high survival (log 15 min kill ≤ 1.5) in the hot alkaline (pH >13) wash liquid, used at the dairies for cleaning-in-place. Second, efficient adherence of the spores to stainless steel from cold water. Third, a cereulide producing group with spores characterized by slow germination in rich medium and well preserved viability when exposed to heating at 90 ºC. Fourth, spores capable of germinating at 8 ºC and possessing the psychrotolerance gene, cspA. There were indications that spores highly resistant to hot 1% sodium hydroxide may be effectively inactivated by hot 0.9% nitric acid. Eight out of the 14 dairy silo tank isolates possessing hot alkali resistant spores were capable of germinating and forming biofilm in whole milk, not previously reported for B. cereus. In this thesis it was shown that cereulide producing B. cereus was capable of inhibiting the growth of cereulide non-producing B. cereus occurring in the same food. This phenomenon, called antagonism, has long been known to exist between B. cereus and other microbial species, e.g. various species of Bacillus, gram-negative bacteria and plant pathogenic fungi. In this thesis intra-species antagonism of B. cereus was shown for the first time. This brother-killing did not depend on the cereulide molecule, also some of the cereulide non-producers were potent antagonists. Interestingly, the antagonistic clades were most frequently found in isolates from food implicated with human illness. The antagonistic property was therefore proposed in this thesis as a novel virulence factor that increases the human morbidity of the species B. cereus, in particular of the cereulide producers.

Biodiversity and Phylogeny of Planktic Cyanobacteria in Temperate Freshwater Lakes

Currently, the classification used for cyanobacteria is based mainly on morphology. In many cases the classification is known to be incongruent with the phylogeny of cyanobacteria. The evaluation of this classification is complicated by the fact that numerous strains are only described morphologically and have not been isolated. Moreover, the phenotype of many cyanobacterial strains alters during prolonged laboratory cultivation. In this thesis, cyanobacterial strains were isolated from lakes (mainly Lake Tuusulanjärvi) and both morphology and phylogeny of the isolates were investigated. The cyanobacterial community composition in Lake Tuusulanjärvi was followed for two years in order to relate the success of cyanobacterial phenotypes and genotypes to environmental conditions. In addition, molecular biological methods were compared with traditional microscopic enumeration and their ability and usefulness in describing the cyanobacterial diversity was evaluated. The Anabaena, Aphanizomenon, and Trichormus strains were genetically heterogeneous and polyphyletic. The phylogenetic relationships of the heterocytous cyanobacteria were not congruent with their classification. In contrast to heterocytous cyanobacteria, the phylogenetic relationships of the Snowella and Woronichinia strains, which had not been studied before this thesis, reflected the morphology of strains and followed their current classification. The Snowella strains formed a monophyletic cluster, which was most closely related to the Woronichinia strain. In addition, a new cluster of thin, filamentous cyanobacterial strains identified as Limnothrix redekei was revealed. This cluster was not closely related to any other known cyanobacteria. The cyanobacterial community composition in Lake Tuusulanjärvi was studied with molecular methods [denaturant gradient gel electrophoresis (DGGE) and cloning of the 16S rRNA gene], through enumerations of cyanobacteria under microscope, and by strain isolations. Microcystis, Anabaena/Aphanizomenon, and Synechococcus were the major groups in the cyanobacterial community in Lake Tuusulanjärvi during the two-year monitoring period. These groups showed seasonal succession, and their success was related to different environmental conditions. The major groups of the cyanobacterial community were detected by all used methods. However, cloning gave higher estimates than microscopy for the proportions of heterocytous cyanobacteria and Synechococcus. The differences were probably caused by the high 16S rRNA gene copy numbers in heterotrophic cyanobacteria and by problems in the identification and detection of unicellular cyanobacteria.

Associations of gene polymorphisms and nutrition with calcium homeostasis and bone mineral density : Studies on skeletal nutrigenetics

Heredity explains a major part of the variation in calcium homeostasis and bone strength, and the susceptibility to osteoporosis is polygenetically regulated. Bone phenotype results from the interplay between lifestyle and genes, and several nutritional factors modulate bone health throughout life. Thus, nutrigenetics examining the genetic variation in nutrient intake and homeostatic control is an important research area in the etiology of osteoporosis. Despite continuing progress in the search for candidate genes for osteoporosis, the results thus far have been inconclusive. The main objective of this thesis was to investigate the associations of lactase, vitamin D receptor (VDR), calcium sensing receptor (CaSR) and parathyroid hormone (PTH) gene polymorphisms and lifestyle factors and their interactions with bone health in Finns at varying stages of the skeletal life span. Markers of calcium homeostasis and bone remodelling were measured from blood and urine samples. Bone strength was measured at peripheral and central bone sites. Lifestyle factors were assessed with questionnaires and interviews. Genetic lactase non-persistence (the C/C-13910 genotype) was associated with lower consumption of milk from childhood, predisposing females in particular to inadequate calcium intake. Consumption of low-lactose milk and milk products was shown to decrease the risk for inadequate calcium intake. In young adulthood, bone loss was more common in males than in females. Males with the lactase C/C-13910 genotype may be more susceptible to bone loss than males with the other lactase genotypes, although calcium intake predicts changes in bone mass more than the lactase genotype. The BsmI and FokI polymorphisms of the VDR gene were associated with bone mass in growing adolescents, but the associations weakened with age. In young adults, the A986S polymorphism of the calcium sensing receptor gene was associated with serum ionized calcium concentrations, and the BstBI polymorphism of the parathyroid gene was related to bone strength. The FokI polymorphism and sodium intake showed an interaction effect on urinary calcium excretion. A novel gene-gene interaction between the VDR FokI and PTH BstBI gene polymorphisms was found in the regulation of PTH secretion and urinary calcium excretion. Further research should be carried out with more number of Finns at varying stages of the skeletal life span and more detailed measurements of bone strength. Research should concern mechanisms by which genetic variants affect calcium homeostasis and bone strength, and the role of diet-gene and gene-gene interactions in the pathogenesis of osteoporosis.

A dynamic programming model for optimising feeding and slaughter decisions regarding fattening pigs

Costs of purchasing new piglets and of feeding them until slaughter are the main variable expenditures in pig fattening. They both depend on slaughter intensity, the nature of feeding patterns and the technological constraints of pig fattening, such as genotype. Therefore, it is of interest to examine the effect of production technology and changes in input and output prices on feeding and slaughter decisions. This study examines the problem by using a dynamic programming model that links genetic characteristics of a pig to feeding decisions and the timing of slaughter and takes into account how these jointly affect the quality-adjusted value of a carcass. The model simulates the growth mechanism of a pig under optional feeding and slaughter patterns and then solves the optimal feeding and slaughter decisions recursively. The state of nature and the genotype of a pig are known in the analysis. The main contribution of this study is the dynamic approach that explicitly takes into account carcass quality while simultaneously optimising feeding and slaughter decisions. The method maximises the internal rate of return to the capacity unit. Hence, the results can have vital impact on competitiveness of pig production, which is known to be quite capital-intensive. The results suggest that producer can significantly benefit from improvements in the pig's genotype, because they improve efficiency of pig production. The annual benefits from obtaining pigs of improved genotype can be more than €20 per capacity unit. The annual net benefits of animal breeding to pig farms can also be considerable. Animals of improved genotype can reach optimal slaughter maturity quicker and produce leaner meat than animals of poor genotype. In order to fully utilise the benefits of animal breeding, the producer must adjust feeding and slaughter patterns on the basis of genotype. The results suggest that the producer can benefit from flexible feeding technology. The flexible feeding technology segregates pigs into groups according to their weight, carcass leanness, genotype and sex and thereafter optimises feeding and slaughter decisions separately for these groups. Typically, such a technology provides incentives to feed piglets with protein-rich feed such that the genetic potential to produce leaner meat is fully utilised. When the pig approaches slaughter maturity, the share of protein-rich feed in the diet gradually decreases and the amount of energy-rich feed increases. Generally, the optimal slaughter weight is within the weight range that pays the highest price per kilogram of pig meat. The optimal feeding pattern and the optimal timing of slaughter depend on price ratios. Particularly, an increase in the price of pig meat provides incentives to increase the growth rates up to the pig's biological maximum by increasing the amount of energy in the feed. Price changes and changes in slaughter premium can also have large income effects. Key words: barley, carcass composition, dynamic programming, feeding, genotypes, lean, pig fattening, precision agriculture, productivity, slaughter weight, soybeans

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

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

Computational Techniques for Haplotype Inference and for Local Alignment Significance

This thesis which consists of an introduction and four peer-reviewed original publications studies the problems of haplotype inference (haplotyping) and local alignment significance. The problems studied here belong to the broad area of bioinformatics and computational biology. The presented solutions are computationally fast and accurate, which makes them practical in high-throughput sequence data analysis. Haplotype inference is a computational problem where the goal is to estimate haplotypes from a sample of genotypes as accurately as possible. This problem is important as the direct measurement of haplotypes is difficult, whereas the genotypes are easier to quantify. Haplotypes are the key-players when studying for example the genetic causes of diseases. In this thesis, three methods are presented for the haplotype inference problem referred to as HaploParser, HIT, and BACH. HaploParser is based on a combinatorial mosaic model and hierarchical parsing that together mimic recombinations and point-mutations in a biologically plausible way. In this mosaic model, the current population is assumed to be evolved from a small founder population. Thus, the haplotypes of the current population are recombinations of the (implicit) founder haplotypes with some point--mutations. HIT (Haplotype Inference Technique) uses a hidden Markov model for haplotypes and efficient algorithms are presented to learn this model from genotype data. The model structure of HIT is analogous to the mosaic model of HaploParser with founder haplotypes. Therefore, it can be seen as a probabilistic model of recombinations and point-mutations. BACH (Bayesian Context-based Haplotyping) utilizes a context tree weighting algorithm to efficiently sum over all variable-length Markov chains to evaluate the posterior probability of a haplotype configuration. Algorithms are presented that find haplotype configurations with high posterior probability. BACH is the most accurate method presented in this thesis and has comparable performance to the best available software for haplotype inference. Local alignment significance is a computational problem where one is interested in whether the local similarities in two sequences are due to the fact that the sequences are related or just by chance. Similarity of sequences is measured by their best local alignment score and from that, a p-value is computed. This p-value is the probability of picking two sequences from the null model that have as good or better best local alignment score. Local alignment significance is used routinely for example in homology searches. In this thesis, a general framework is sketched that allows one to compute a tight upper bound for the p-value of a local pairwise alignment score. Unlike the previous methods, the presented framework is not affeced by so-called edge-effects and can handle gaps (deletions and insertions) without troublesome sampling and curve fitting.

Computational Identification of Recessive Mutations in Cancers using High Throughput SNP-arrays

This thesis presents a highly sensitive genome wide search method for recessive mutations. The method is suitable for distantly related samples that are divided into phenotype positives and negatives. High throughput genotype arrays are used to identify and compare homozygous regions between the cohorts. The method is demonstrated by comparing colorectal cancer patients against unaffected references. The objective is to find homozygous regions and alleles that are more common in cancer patients. We have designed and implemented software tools to automate the data analysis from genotypes to lists of candidate genes and to their properties. The programs have been designed in respect to a pipeline architecture that allows their integration to other programs such as biological databases and copy number analysis tools. The integration of the tools is crucial as the genome wide analysis of the cohort differences produces many candidate regions not related to the studied phenotype. CohortComparator is a genotype comparison tool that detects homozygous regions and compares their loci and allele constitutions between two sets of samples. The data is visualised in chromosome specific graphs illustrating the homozygous regions and alleles of each sample. The genomic regions that may harbour recessive mutations are emphasised with different colours and a scoring scheme is given for these regions. The detection of homozygous regions, cohort comparisons and result annotations are all subjected to presumptions many of which have been parameterized in our programs. The effect of these parameters and the suitable scope of the methods have been evaluated. Samples with different resolutions can be balanced with the genotype estimates of their haplotypes and they can be used within the same study.

Polymorphic low penetrance genes and breast cancer : The role of genes involved in metabolism of xenobiotics, estrogens and reactive oxygen species

Various endogenous and exogenous factors have been reported to increase the risk of breast cancer. Many of those are related to prolonged lifetime exposure to estrogens. Furthermore, a positive family history of breast cancer and certain benign breast diseases are known to increase the risk of breast cancer. The role of lifestyle factors, such as use of alcohol and smoking has been an area of intensive study. Alcohol has been found to increase the risk of breast cancer, whereas the role of smoking has remained obscure. A multitude of enzymes are involved in the metabolism of estrogens and xenobiotics including the carcinogens found in tobacco smoke. Many of the metabolic enzymes exhibit genetic polymorphisms that can lead to inter-individual differences in their abilities to modify hazardous substrates. Therefore, in presence of a given chemical exposure, one subgroup of women may be more susceptible to breast carcinogenesis, since they carry unfavourable forms of the polymorphic genes involved in the metabolism of the chemical. In this work, polymorphic genes encoding for cytochrome P450 (CYP) 1A1 and 1B1, N-acetyl transferase 2 (NAT2), sulfotransferase 1A1 (SULT1A1), manganese superoxide dismutase (MnSOD) and vitamin D receptor (VDR) were investigated in relation to breast cancer susceptibility in a Finnish population. CYP1A1, CYP1B1 and SULT1A1 are involved in the metabolism of both estrogens and xenobiotics, whereas NAT2 is involved only in the latter. MnSOD is an antioxidant enzyme protecting cells from oxidative damage. VDR, in turn, mediates the effects of the active form of vitamin D (1,25(OH)2D3, calcitriol) on maintenance of calcium homeostasis and it has anti-proliferative effects in many cancer cells. A 1.3-fold (95% CIs 1.01-1.73) increased risk of breast cancer was seen among women who carried the NAT2 slow acetylator genotype and a 1.5-fold (95% CI 1.1-2.0) risk was found in women with a MnSOD variant A allele containing genotypes compared to women with the NAT2 rapid acetylator genotype or to those with the MnSOD VV genotype, respectively. Instead, women with the VDR a allele containing genotypes were found to be at a decreased risk for breast cancer (OR 0.73; 95% CI 0.54-0.98) compared to women with the AA genotype. No significant overall associations were found between SULT1A1 or CYP genotypes and breast cancer risk, whereas a combination of the CYP1B1 432Val allele containing genotypes with the NAT2 slow acetylator genotypes posed a 1.5-fold (95% CI 1.03-2.24) increased risk. Moreover, NAT2 slow acetylator genotype was found to be confined to women with an advanced stage of breast cancer (stages III and IV). Further evidence for the association of xenobiotic metabolising genes with breast cancer risk was found when active smoking was taken into account. Women who smoked less than 10 cigarettes/day and carried at least one CYP1B1 432Val variant allele, were at 3.1-fold (95% CI 1.32-7.12) risk of breast cancer compared to women who smoked the same amount but did not carry the variant allele. Furthermore, the risk was significantly increased with increasing number of the CYP1B1 432Val alleles (p for trend 0.005). In addition, women who smoked less than 5 pack-years and carried the NAT2 slow acetylator genotype were at a 2.6-fold (95% CI 1.01-6.48) increased risk of breast cancer compared to women who smoked the same amount but carried the NAT2 rapid acetylator genotype. Furthermore, the combination of the CYP1B1 432Val allele and the NAT2 slow acetylator genotype increased the risk of breast cancer by 2.5-fold (95% CI 1.11-5.45) among ever smokers. Instead, the MnSOD A allele was found to be a risk factor among postmenopausal long-term smokers (>15 years of smoking) (OR 5.1; 95% CI 1.4-18.4) or among postmenopausal women who had smoked more than 10 cigarettes/day (OR 5.5; 95% CI 1.3-23.4) compared to women who had similar smoking habits but carried the MnSOD V/V genotype. Similarly, within subgroups of postmenopausal women who were using oral contraceptives, hormone replacement therapy or alcohol, women carrying the MnSOD A allele genotypes seemed to be at increased risk of breast cancer compared to women with the MnSOD V/V genotype. A positive family history of breast cancer and high parity were shown to be inversely associated with breast cancer risk among women carrying the VDR ApaI a allele or among premenopausal women carrying the SULT1A1*2 allele, respectively.

Flight metabolic rate in the Glanville fritillary butterfly

Dispersal is a highly important life history trait. In fragmented landscapes the long-term persistence of populations depends on dispersal. Evolution of dispersal is affected by costs and benefits and these may differ between different landscapes. This results in differences in the strength and direction of natural selection on dispersal in fragmented landscapes. Dispersal has been shown to be a nonrandom process that is associated with traits such as flight ability in insects. This thesis examines genetic and physiological traits affecting dispersal in the Glanville fritillary butterfly (Melitaea cinxia). Flight metabolic rate is a repeatable trait representing flight ability. Unlike in many vertebrates, resting metabolic rate cannot be used as a surrogate of maximum metabolic rate as no strong correlation between the two was found in the Glanville fritillary. Resting and flight metabolic rate are affected by environmental variables, most notably temperature. However, only flight metabolic rate has a strong genetic component. Molecular variation in the much-studied candidate locus phosphoglucose isomerase (Pgi), which encodes the glycolytic enzyme PGI, has an effect on carbohydrate metabolism in flight. This effect is temperature dependent: in low to moderate temperatures individuals with the heterozygous genotype at the single nucleotide polymorphism (SNP) AA111 have higher flight metabolic rate than the common homozygous genotype. At high temperatures the situation is reversed. This finding suggests that variation in enzyme properties is indeed translated to organismal performance. High-resolution data on individual female Glanville fritillaries moving freely in the field were recorded using harmonic radar. There was a strong positive correlation between flight metabolic rate and dispersal rate. Flight metabolic rate explained one third of the observed variation in the one-hour movement distance. A fine-scaled analysis of mobility showed that mobility peaked at intermediate ambient temperatures but the two common Pgi genotypes differed in their reaction norms to temperature. As with flight metabolic rate, heterozygotes at SNP AA111 were the most active genotype in low to moderate temperatures. The results show that molecular variation is associated with variation in dispersal rate through the link of flight physiology under the influence of environmental conditions. The evolutionary pressures for dispersal differ between males and females. The effect of flight metabolic rate on dispersal was examined in both sexes in field and laboratory conditions. The relationship between flight metabolic rate and dispersal rate in the field and flight duration in the laboratory were found to differ between the two sexes. In females the relationship was positive, but in males the longest distances and flight durations were recorded for individuals with low flight metabolic rate. These findings may reflect male investment in mate locating. Instead of dispersing, males with high flight metabolic rate may establish territories and follow a perching strategy when locating females and hence move less on the landscape level. Males with low metabolic rate may be forced to disperse due to low competitive success or may show adaptations to an alternative strategy: patrolling. In the light of life history trade-offs and the rate of living theory having high metabolic rate may carry a cost in the form of shortened lifespan. Experiments relating flight metabolic rate to longevity showed a clear correlation in the opposite direction: high flight metabolic rate was associated with long lifespan. This suggests that individuals with high metabolic rate do not pay an extra physiological cost for their high flight capacity, rather there are positive correlations between different measures of fitness. These results highlight the importance of condition.

Ecological consequences of genetic modifications - an invasion analysis approach

Biological invasions are considered as one of the greatest threats to biodiversity, as they may lead to disruption and homogenization of natural communities, and in the worst case, to native species extinctions. The introduction of gene modified organisms (GMOs) to agricultural, fisheries and forestry practices brings them into contact with natural populations. GMOs may appear as new invasive species if they are able to (1) invade into natural habitats or (2) hybridize with their wild relatives. The benefits of GMOs, such as increased yield or decreased use of insecticides or herbicides in cultivation, may thus be reduced due the potential risks they may cause. A careful ecological risk analysis therefore has to precede any responsible GMO introduction. In this thesis I study ecological invasion in relation to GMOs, and what kind of consequences invasion may have in natural populations. A set of theoretical models that combine life-history evolution, population dynamics, and population genetics were developed for the hazard identification part of ecological risks assessment of GMOs. In addition, the potential benefits of GMOs in management of an invasive pest were analyzed. In the first study I showed that a population that is fluctuating due to scramble-type density dependence (due to, e.g., nutrient competition in plants) may be invaded by a population that is relatively more limited by a resource (e.g., light in plants) that is a cause of contest-type density dependence. This result emphasises the higher risk of invasion in unstable environments. The next two studies focused on escape of a growth hormone (GH) transgenic fish into a natural population. The results showed that previous models may have given too pessimistic a view of the so called Trojan gene -effect, where the invading genotype is harmful for the population as a whole. The previously suggested population extinctions did not occur in my studies, since the changes in mating preferences caused by the GH-fish were be ameliorated by decreased level of competition. The GH-invaders may also have to exceed a threshold density before invasion can be successful. I also showed that the prevalence of mature parr (aka. sneaker) strategy among GH-fish may have clear effect on invasion outcome. The fourth study assessed the risks and developed methods against the invasion of the Colorado Potato Beetle (CPB, Leptinotarsa decemlineata). I showed that the eradication of CPB is most important for the prevention of their establishment, but the cultivation of transgenic Bt-potato could also be effective. In general, my results emphasise that invasion of transgenic species or genotypes to be possible under certain realistic conditions and resulting in competitive exclusion, population decline through outbreeding depression and genotypic displacement of native species. Ecological risk assessment should regard the decline and displacement of the wild genotype by an introduced one as a consequence that is as serious as the population extinction. It will also be crucial to take into account different kinds of behavioural differences among species when assessing the possible hazards that GMOs may cause if escaped. The benefits found of GMO crops effectiveness in pest management may also be too optimistic since CPB may evolve resistance to Bt-toxin. The models in this thesis could be further applied in case specific risk assessment of GMOs by supplementing them with detailed data of the species biology, the effect of the transgene introduced to the species, and also the characteristics of the populations or the environments in the risk of being invaded.

Human parvovirus B19: tissue persistence and prevalence of prototypic and new variants

Human parvovirus B19 is a minute ssDNA virus causing a wide variety of diseases, including erythema infectiosum, arthropathy, anemias, and fetal death. After primary infection, genomic DNA of B19 has been shown to persist in solid tissues of not only symptomatic but also of constitutionally healthy, immunocompetent individuals. In this thesis, the viral DNA was shown to persist as an apparently intact molecule of full length, and without persistence-specific mutations. Thus, although the mere presence of B19 DNA in tissue can not be used as a diagnostic criterion, a possible role in the pathogenesis of diseases e.g. through mRNA or protein production can not be excluded. The molecular mechanism, the host-cell type and the possible clinical significance of B19 DNA tissue persistence are yet to be elucidated. In the beginning of this work, the B19 genomic sequence was considered highly conserved. However, new variants were found: V9 was detected in 1998 in France, in serum of a child with aplastic crisis. This variant differed from the prototypic B19 sequences by ~10 %. In 2002 we found, persisting in skin of constitutionally healthy humans, DNA of another novel B19 variant, LaLi. Genetically this variant differed from both the prototypic sequences and the variant V9 also by ~10%. Simultaneously, B19 isolates with DNA sequences similar to LaLi were introduced by two other groups, in the USA and France. Based on phylogeny, a classification scheme based on three genotypes (B19 types 1-3) was proposed. Although the B19 virus is mainly transmitted via the respiratory route, blood and plasma-derived products contaminated with high levels of B19 DNA have also been shown to be infectious. The European Pharmacopoeia stipulates that, in Europe, from the beginning of 2004, plasma pools for manufacture must contain less than 104 IU/ml of B19 DNA. Quantitative PCR screening is therefore a prerequisite for restriction of the B19 DNA load and obtaining of safe plasma products. Due to the DNA sequence variation among the three B19 genotypes, however, B19 PCR methods might fail to detect the new variants. We therefore examined the suitability of the two commercially available quantitative B19 PCR tests, LightCycler-Parvovirus B19 quantification kit (Roche Diagnostics) and RealArt Parvo B19 LC PCR (Artus), for detection, quantification and differentiation of the three B19 types known, including B19 types 2 and 3. The former method was highly sensitive for detection of the B19 prototype but was not suitable for detection of types 2 and 3. The latter method detected and differentiated all three B19 virus types. However, one of the two type-3 strains was detected at a lower sensitivity. Then, we assessed the prevalence of the three B19 virus types among Finnish blood donors, by screening pooled plasma samples derived from >140 000 blood-donor units: none of the pools contained detectable levels of B19 virus types 2 or 3. According to the results of other groups, B19 type 2 was absent also among Danish blood-donors, and extremely rare among symptomatic European patients. B19 type 3 has been encountered endemically in Ghana and (apparently) in Brazil, and sporadical cases have been detected in France and the UK. We next examined the biological characteristics of these virus types. The p6 promoter regions of virus types 1-3 were cloned in front of a reporter gene, the constructs were transfected into different cell lines, and the promoter activities were measured. As a result, we found that the activities of the three p6 promoters, although differing in sequence by >20%, were of equal strength, and most active in B19-permissive cells. Furthermore, the infectivity of the three B19 types was examined in two B19-permissive cell lines. RT-PCR revealed synthesis of spliced B19 mRNAs, and immunofluorescence verified the production of NS1 and VP proteins in the infected cells. These experiments suggested similar host-cell tropism and showed that the three virus types are strains of the same species, i.e. human parvovirus B19. Last but not least, the sera from subjects infected in the past either with B19 type 1 or type 2 (as evidenced by tissue persistence of the respective DNAs), revealed in VP1/2- and VP2-EIAs a 100 % cross-reactivity between virus types 1 and 2. These results, together with similar studies by others, indicate that the three B19 genotypes constitute a single serotype.

Comparative and functional genome analysis of fungi for development of the protein production host Trichoderma reesei

Filamentous fungi of the subphylum Pezizomycotina are well known as protein and secondary metabolite producers. Various industries take advantage of these capabilities. However, the molecular biology of yeasts, i.e. Saccharomycotina and especially that of Saccharomyces cerevisiae, the baker's yeast, is much better known. In an effort to explain fungal phenotypes through their genotypes we have compared protein coding gene contents of Pezizomycotina and Saccharomycotina. Only biomass degradation and secondary metabolism related protein families seem to have expanded recently in Pezizomycotina. Of the protein families clearly diverged between Pezizomycotina and Saccharomycotina, those related to mitochondrial functions emerge as the most prominent. However, the primary metabolism as described in S. cerevisiae is largely conserved in all fungi. Apart from the known secondary metabolism, Pezizomycotina have pathways that could link secondary metabolism to primary metabolism and a wealth of undescribed enzymes. Previous studies of individual Pezizomycotina genomes have shown that regardless of the difference in production efficiency and diversity of secreted proteins, the content of the known secretion machinery genes in Pezizomycotina and Saccharomycotina appears very similar. Genome wide analysis of gene products is therefore needed to better understand the efficient secretion of Pezizomycotina. We have developed methods applicable to transcriptome analysis of non-sequenced organisms. TRAC (Transcriptional profiling with the aid of affinity capture) has been previously developed at VTT for fast, focused transcription analysis. We introduce a version of TRAC that allows more powerful signal amplification and multiplexing. We also present computational optimisations of transcriptome analysis of non-sequenced organism and TRAC analysis in general. Trichoderma reesei is one of the most commonly used Pezizomycotina in the protein production industry. In order to understand its secretion system better and find clues for improvement of its industrial performance, we have analysed its transcriptomic response to protein secretion stress conditions. In comparison to S. cerevisiae, the response of T. reesei appears different, but still impacts on the same cellular functions. We also discovered in T. reesei interesting similarities to mammalian protein secretion stress response. Together these findings highlight targets for more detailed studies.

Molecular epidemiology of human rhinoviruses

The first part of this work investigates the molecular epidemiology of a human enterovirus (HEV), echovirus 30 (E-30). This project is part of a series of studies performed in our research team analyzing the molecular epidemiology of HEV-B viruses. A total of 129 virus strains had been isolated in different parts of Europe. The sequence analysis was performed in three different genomic regions: 420 nucleotides (nt) in the VP4/VP2 capsid protein coding region, the entire VP1 capsid protein coding gene of 876 nt, and 150 nt in the VP1/2A junction region. The analysis revealed a succession of dominant sublineages within a major genotype. The temporally earlier genotypes had been replaced by a genetically homogenous lineage that has been circulating in Europe since the late 1970s. The same genotype was found by other research groups in North America and Australia. Globally, other cocirculating genetic lineages also exist. The prevalence of a dominant genotype makes E-30 different from other previously studied HEVs, such as polioviruses and coxsackieviruses B4 and B5, for which several coexisting genetic lineages have been reported. The second part of this work deals with molecular epidemiology of human rhinoviruses (HRVs). A total of 61 field isolates were studied in the 420-nt stretch in the capsid coding region of VP4/VP2. The isolates were collected from children under two years of age in Tampere, Finland. Sequences from the clinical isolates clustered in the two previously known phylogenetic clades. Seasonal clustering was found. Also, several distinct serotype-like clusters were found to co-circulate during the same epidemic season. Reappearance of a cluster after disappearing for a season was observed. The molecular epidemiology of the analyzed strains turned out to be complex, and we decided to continue our studies of HRV. Only five previously published complete genome sequences of HRV prototype strains were available for analysis. Therefore, all designated HRV prototype strains (n=102) were sequenced in the VP4/VP2 region, and the possibility of genetic typing of HRV was evaluated. Seventy-six of the 102 prototype strains clustered in HRV genetic group A (HRV-A) and 25 in group B (HRV-B). Serotype 87 clustered separately from other HRVs with HEV species D. The field strains of HRV represented as many as 19 different genotypes, as judged with an approximate demarcation of a 20% nt difference in the VP4/VP2 region. The interserotypic differences of HRV were generally similar to those reported between different HEV serotypes (i.e. about 20%), but smaller differences, less than 10%, were also observed. Because some HRV serotypes are genetically so closely related, we suggest that the genetic typing be performed using the criterion "the closest prototype strain". This study is the first systematic genetic characterization of all known HRV prototype strains, providing a further taxonomic proposal for classification of HRV. We proposed to divide the genus Human rhinoviruses into HRV-A and HRV-B. The final part of the work comprises a phylogenetic analysis of a subset (48) of HRV prototype strains and field isolates (12) in the nonstructural part of the genome coding for the RNA-dependent RNA polymerase (3D). The proposed division of the HRV strains in the species HRV-A and HRV-B was also supported by 3D region. HRV-B clustered closer to HEV species B, C, and also to polioviruses than to HRV-A. Intraspecies variation within both HRV-A and HRV-B was greater in the 3D coding region than in the VP4/VP2 coding region, in contrast to HEV. Moreover, the diversity of HRV in 3D exceeded that of HEV. One group of HRV-A, designated HRV-A', formed a separate cluster outside other HRV-A in the 3D region. It formed a cluster also in the capsid region, but located within HRV-A. This may reflect a different evolutionary history of distinct genomic regions among HRV-A. Furthermore, the tree topology within HRV-A in the 3D region differed from that in the VP4/VP2, suggesting possible recombination events in the evolution of the strains. No conflicting phylogenies were observed in any of the 12 field isolates. Possible recombination was further studied using the Similarity and Bootscanning analyses of the complete genome sequences of HRV available in public databases. Evidence for recombination among HRV-A was found, as HRV2 and HRV39 showed higher similarity in the nonstructural part of the genome. Whether HRV2 and HRV39 strains - and perhaps also some other HRV-A strains not yet completely sequenced - are recombinants remains to be determined.

Herpes Simplex Virus Infection, Pathological Pain and Recurrent Lymphocytic Meningitis

Background: Aims of the study were: (i) to characterise the clinical picture, immunological features and changes in brain morphology and function in patients with widespread unilateral pain and HSV-infections, and (ii) to analyse the prevalence, clinical symptoms and immunological predisposing factors of HSV-2 induced recurrent lymphocytic meningitis (RLM) in Southern Finland. Patients and methods: Patients for the studies were recruited from the Pain Clinic, and from the Department of Neurology, at Helsinki University Central Hospital. Plasma concentrations of IgM, IgA, IgG, and IgG1-4, and serum concentrations of C3, C4 were measured. Serological anti-HSV-1 and -2 antibody status was tested. C4 genotyping, HLA-A, HLA-B and HLA-DRB1 typing, MBL2 genotyping, and IgG1 and IgG3 allotyping (Gm) were performed. Clinical neurological examination, quantitative sensory testing, skin biopsy, and functional magnetic resonance imaging were also performed. Results: HSV probably has a role in the generation of a pathological pain state. Low serum IgG1 and IgG3 levels, made the patients vulnerable for recurring HSV infections. Both functional and structural changes were observed in the brain pain-processing areas in the patients: they had less pain-related activity in the insular cortices bilaterally, in the anterior cingular cortex (ACC), and in the thalamus, and the gray matter density was lower in the ACC, in the frontal and prefrontal cortices. In the meningitis studies it was shown that RLM is more common and less benign than previously reported, and that neuropathic pain is frequently present both during and after meningitis episodes. HLA-DRB1*01, HLA-B*27, and low IgG1 levels are predisposing factors for RLM. Conclusions: Patients are vulnerable to recurrent HSV infections because of subtle immunological abnormalities. HSV causes diverse clinical manifestations. First, the herpes simplex virus, or the inflammatory process triggered by it, may cause pathological widespread pain probably by activating glial cells in the CNS. In these patients, signs of alterations in the brain pain-processing areas can be demonstrated by functional brain imaging methods. Secondly, HSV-2 induced RLM is a rare complication of HSV-2 virus. The predisposing factors include low IgG1 subclass levels, HLA-DRB1*01 and HLA –B*27 genotypes. Neuropathic pain is frequently associated with RLM.