The DISC1 Pathway in the Genetic Etiology of Schizophrenia

Schizophrenia is a severe psychotic disorder affecting 0.5-1 % of the population. The disorder is characterized by hallucinations; delusions; disorganized behavior and speech; avolition; anhedonia; flattened affect and cognitive deficits. The etiology of the disorder is complex with evidence for multiple genes contributing to the onset of the disorder along with environmental factors. DISC1 is one of the most promising candidate genes for schizophrenia. It codes for a protein which takes part in numerous molecular interactions along several pathways. This network, termed as the DISC1 pathway, is evidently important for the development and maturation of the central nervous system from the embryo until young adulthood. Disruption at these pathways is thought to predispose schizophrenia. In the present study, we have studied the DISC1 pathway in the etiology of schizophrenia in the Finnish population. We have utilized large Finnish samples; the schizophrenia family sample where DISC1 was originally shown to associate with schizophrenia and the Northern Finland birth cohort 1966 (NFBC66). Several DISC1 binding partners displayed evidence for association in the family sample along with DISC1. Through a genome-wide linkage study, we found a significant linkage signal to a locus where a DISC1 binding partner NDE1 is located at the carriers of a certain DISC1 risk variant. In a follow-up study, genetic markers in NDE1 displayed significant evidence for association with schizophrenia. Further exploration of association between 11 genes of the DISC1 pathway and schizophrenia led to recognition of novel variants in NDEL1, PDE4B and PDE4D that significantly either increased or decreased the risk for schizophrenia. Further, we found evidence that DISC1 itself has a significant role in the human mental functioning even in the healthy population. Variants in DISC1 had a significant effect on anhedonia which is a trait present at everybody but is in its severe form one of the main symptoms of schizophrenia and correlates with the risk of developing the disorder. Further, utilizing genome-wide marker data, we recognized three genes; MIR620; CCDC141 and LCT; that are closely related to the DISC1 pathway but which effects on anhedonia were observable only at the individuals who carried these specific DISC1 variants. Our findings significantly add up to the previous evidence for the involvement of DISC1 and the DISC1 pathway in the etiology of schizophrenia and psychosis. Our results support the concept of a number of DISC1 pathway related genes contributing in the etiology of schizophrenia along with DISC1 and provide new candidates for the studies of schizophrenia. Our findings also significantly increase the importance of DISC1 itself as having a role in psychological functioning in the general population.

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 molecular basis of Marinesco-Sjögren syndrome

Marinesco-Sjögren syndrome (MSS) is a rare autosomal recessive neurodegenerative disorder characterized by cerebellar ataxia due to cerebellar cortical atrophy, infantile- or childhood-onset bilateral cataracts, progressive myopathy, and mild to severe mental retardation. Additional features include hypergonadotropic hypogonadism, various skeletal abnormalities, short stature, and strabismus. The neuroradiologic hallmarks are hypoplasia of both the vermis and cerebellar hemispheres. The histopathologic findings include severe cerebellar atrophy and loss of Purkinje and granule cells. The common pathologic findings in muscle biopsy are variation in muscle fiber size, atrophic fibers, fatty replacement, and rimmed vacuole formation. The presence of marked cerebellar atrophy with myopathy distinguishes MSS from another rare syndrome, the congenital cataracts, facial dysmorphism, and neuropathy syndrome (CCFDN). Previously, work by others had resulted in the identification of an MSS locus on chromosome 5q31. A subtype of MSS with myoglobinuria and neuropathy had been linked to the CCFDN locus on chromosome 18qter, at which mutations in the CTDP1 gene had been identified. We confirmed linkage to the previously identified locus on chromosome 5q31 in two Finnish families with eight affected individuals, reduced the critical region by fine-mapping, and identified SIL1 as a gene underlying MSS. We found a common homozygous founder mutation in all Finnish patients. The same mutation was also present in patient samples from Norway and Sweden. Altogether, we identified eight mutations in SIL1, including nonsense, frameshift, splice site alterations, and one missense mutation. SIL1 encodes a nucleotide exchange factor for the endoplasmic reticulum (ER) resident heat-shock protein 70 chaperone GRP78. GRP78 functions in protein synthesis and quality control of the newly synthesized polypeptides. It senses and responds to stressful cellular conditions. We showed that in mice, SIL1 and GRP78 show highly similar spatial and temporal tissue expression in developing and mature brain, eye, and muscle. Studying endogenous proteins in mouse primary hippocampal neurons, we found that SIL1 and GRP78 colocalize and that SIL1 localizes to the ER. We studied the subcellular localization of two mutant proteins, a missense mutant found in two patients and an artificial mutant lacking the ER retrieval signal, and found that both mutant proteins formed aggregates within the ER. Well in line with our findings and the clinical features of MSS, recent work by Zhao et al. showed that a truncation of SIL1 causes ataxia and cerebellar Purkinje cell loss in the naturally occurring woozy mutant mouse. Prior to Purkinje cell degeneration, the unfolded protein response is initiated and abnormal protein accumulations are present. MSS thus joins the group of protein misfolding and accumulation diseases. These findings highlight the importance of SIL1 and the role of the ER in neuronal function and survival. The results presented in this thesis provide tools for the molecular genetic diagnostics of MSS and give a basis for future studies on the molecular pathogenesis of MSS. Understanding the mechanisms behind this pleiotropic syndrome may provide insights into more common forms of ataxia, myopathy, and neurodegeneration.

Adenoviral gene therapy for non-small cell lung cancer

Adenoviral gene therapy is an experimental approach to cancer refractory to standard cancer therapies. Adenoviruses can be utilized as vectors to deliver therapeutic transgenes into cancer cells, while gene therapy with oncolytic adenoviruses exploits the lytic potential of viruses to kill tumor cells. Although adenoviruses demonstrate several advantages over other vectors - such as the unparalleled transduction efficacy and natural tropism to a wide range of tissues - the gene transfer efficacy to cancer cells has been limited, consequently restricting the therapeutic effect. There are, however, several approaches to circumvent this problem. We utilized different modified adenoviruses to obtain information on adenovirus tropism towards non-small cell lung cancer (NSCLC) cells. To enhance therapeutic outcome, oncolytic adenoviruses were evaluated. Further, to enhance gene delivery to tumors, we used mesenchymal stem cells (MSCs) as carriers. To improve adenovirus specificity, we investigated whether widely used cyclooxygenase 2 (Cox-2) promoter is induced by adenovirus infection in nontarget cells and whether selectivity can be retained by the 3 untranslated region (UTR) AU-rich elements. In addition, we investigated whether switching adenovirus fiber can retain gene delivery in the presence of neutralizing antibodies. Our results show that adenoviruses, whose capsids were modified with arginine-glycine-aspartatic acid (RGD-4C), the serotype 3 knob, or polylysins displayed enhanced gene transfer into NSCLC cell lines and fresh clinical specimens from patients. The therapeutic efficacy was further improved by using respective oncolytic adenoviruses with isogenic 24bp deletion in the E1A gene. Cox-2 promoter was also shown to be induced in normal and tumor cells following adenovirus infection, but utilization of 3 UTR elements can increase the tumor specificity of the promoter. Further, the results suggested that use of MSCs could enhance the bioavailability and delivery of adenoviruses into human tumors, although cells had no tumor tropism per se. Finally, we demonstrated that changing adenovirus fiber can allow virus to escape from existing neutralizing antibodies when delivered systemically. In conclusion, these results reveal that adenovirus gene transfer and specificity can be increased by using modified adenoviruses and MSCs as carriers, and fiber modifications simultaneously decrease the effect of neutralizing antibodies. This promising data suggest that these approaches could translate into clinical testing in patients with NSCLC refractory to current modalities.

Studies on immune activation in rheumatoid arthritis and reactive arthritis

Rheumatoid arthritis (RA) is an autoimmune disease characterized by synovitis, progressive joint destruction, and disability. Reactive arthritis (ReA) is a sterile joint inflammation following a distant mucosal infection. The clinical course of these diseases is variable and cannot be predicted with reasonable accuracy by clinical and laboratory markers. The predictive value of circulating soluble interleukin-2 receptor (sIL-2R), a marker of lymphocyte activation, measured by Immulite® automated immunoassay analyzer, was evaluated in two cohorts of RA patients. In 175 patients with active early RA randomized to treatment with either on disease-modifying antirheumatic drug (DMARD) or a combination of 3 DMARDs and prednisolone, low baseline sIL-2R level predicted remission after 6 months in patients treated with a single DMARD. In 24 patients with active RA refractory to DMARDs, low baseline sIL-2R level predicted rapid clinical response to treatment with infliximab, an anti-tumour necrosis factor antibody. Furthermore, in a cohort of 26 patients with acute ReA, high baseline sIL-2R level predicted remission after 6 months. Levels of circulating soluble E-selectin (sE-selectin), a marker of endothelial activation, were measured annually by enzyme-linked immunosorbent assay (ELISA) in a cohort of 85 patients with early RA. During a five-year follow-up, sE-selectin levels were associated with activity and outcome of RA. The levels of neutrophil and monocyte CD11b/CD18 expression measured by flow cytometry, and circulating levels of sE-selectin measured by ELISA, and procalcitonin by immunoluminometric assay, were compared in 28 patients with acute ReA and 16 patients with early RA. The levels of the markers were comparable in ReA, RA, and healthy control subjects. In conlusion, sIL-2R may provide a new predictive marker in early RA treated with a single DMARD and refractory RA treated with infliximab. In addition, sIL-2R level predicts remission in acute ReA.

Target Cell Topography and Cytoskeletal Reorganization in Natural Killer Cell Activity

The immune system has to recognize and destroy abnormal or infected cells to maintain homeostasis. Natural killer (NK) cells directly recognize and kill transformed or virus-infected cells without prior sensitization. We have studied both virus-infected and tumor cells in order to identify the target structures involved in triggering NK activity. Mouse/human cell hybrids containing various human chromosomes were used as targets. The human chromosome responsible for activating NK cell killing was identified to chromosome number 6. The results suggest that activated NK cells recognize ligands that are encoded on human chromosome 6. We showed that the ligand on the target cell side was intercellular adhesion molecule 2 (ICAM-2). There was no difference in the level of expression of ICAM-2, however, but a drastic difference was seen in the distribution of the molecule: ICAM-2 was evenly distributed on the surface of the NK-resistant cells, but almost totally redistributed to the tip of uropods, bud-like extensions, which were absent from the parental cells. Interestingly, the gene coding for cytoskeletal linker protein ezrin has been localized to human chromosome 6, and there was a colocalization of ezrin and ICAM-2 in the uropods. Furthermore, the transfected human ezrin into NK cell-resistant cells induced uropod formation, ICAM-2 and ezrin redistribution to newly formed uropods, and sensitized target cells to NK cell killing. These data reveal a novel form of NK cell recognition: target structures are already present on normal cells; they become detectable only after abnormal redistribution into hot spots on the target cell membrane. NK cells are central players in the defence against virus infections. They inhibit the spread of infection, allowing time for specific immune responses to develop. The virus-proteins that directly activate human NK cell killing are largely unknown. We studied the sensitivity of virus-specific early proteins of Semliki Forest virus (SFV) to NK killing. The viral non-structural proteins (nsP1-4) translated early in the virus cycle were transfected in NK-resistant cells. Viral early gene nsP1 alone efficiently sensitized target cells to NK activity, and the tight membrane association of nsP1 seems to be critical in the triggering of NK killing. NsP1 protein colocalized with (redistributed) ezrin in filopodia-like structures to which the NK cells were bound. The results suggest that also in viral infections NK cells react to rapid changes in membrane topography. Based on the results of this thesis, a new model of target cell recognition of NK cells can be suggested: reorganization of the cytoskeleton induces alterations in cell surface topography, and this new pattern of surface molecules is recognized as "altered-self".

Characterization of the TRIM37 gene and mutations underlying mulibrey nanism

Mulibrey nanism is a hereditary developmental disorder, characterized by prenatal onset growth failure without postnatal catch-up growth, distinctive craniofacial features, progressive cardiopathy and failure of sexual maturation. In addition, the patients develop insulin resistance syndrome and type 2 diabetes and they have an increased risk of developing tumors. The TRIM37 gene that underlies mulibrey nanism encodes for a member of the tripartite motif (TRIM) protein family. The physiological function of TRIM37 and the pathogenetic mechanisms leading from TRIM37 dysfunction to the mulibrey nanism phenotype are unknown. However, TRIM37 localizes at least partially to peroxisomes, and possesses ubiquitin E3-ligase activity. Thus, it may mediate ubiquitin dependent protein degradation, suggesting that accumulation of yet unknown substrate proteins may underlie the disease pathogenesis. In this study, the TRIM37 gene was characterized in detail. A transcription initiation window, with several separate transcription start sites, was identified and the putative promoter region immediately upstream from the transcription initiation window was shown to possess basal promoter activity. Further, several alternative splice variants of the gene were identified, including a highly expressed testis specific variant, encoding for an identical protein product with the main transcript. Expression of TRIM37 mRNA was detected in several different tissues, with highest expression seen in testis and in brain, when the expression patterns of the two major transcripts in different human tissues were studied by quantitative real-time PCR. Several mulibrey nanism patients were studied and thirteen novel mutations in TRIM37 were found, including three mutations (p.Gly322Val, p.Cys109Ser, p.Glu271_Ser287), that are likely to express mutant TRIM37 proteins. These mutations were further shown to alter the subcellular localization of the mutant proteins. Most of the mulibrey nanism associated mutations however, lead to premature termination codons and degradation of mRNA. All the TRIM37 mutations identified to date predict loss-of-function alleles, and thus no phenotype-genotype correlation is seen among the patients. In order to understand the pathogenetic mechanisms underlying mulibrey nanism, an animal model for the disorder is needed. For the development of a Trim37 knock-out mouse, the mouse Trim37 gene was characterized. Alternative splice variants, were identified, including a testis specific variant predicting a longer protein product. Further, a strictly tissue and cell-specific pattern of Trim37 expression was observed in developing and adult mouse tissues, when studied by immunohistochemical methods. This distribution of Trim37 expression in mouse tissues is in agreement with the clinical findings in human mulibrey nanism patients. This thesis work gives new tools for the diagnostics of mulibrey nanism as well as for studying the molecular pathogenesis behind this interesting disorder.

The role of AIP and CDKN1B/p27Kip1 in endocrine neoplasia

Identification of genes predisposing to tumor syndromes has raised general awareness of tumorigenesis. Genetic testing of tumor susceptibility genes aids the recognition of individuals at increased risk of tumors. Identification of novel predisposing genes enables further studies concerning the classification of potential associated tumors and the definition of target patient group. Pituitary adenomas are common, benign neoplasms accounting for approximately 15% of all intracranial tumors. Accurate incidence estimation is challenging since a great portion of these adenomas are small and asymptomatic. Clinically relevant adenomas, that cause symptoms due to the expansion of the cell mass or the over-secretion of normally produced hormones, occur in approximately one of 1 000 individuals. Although the majority of pituitary adenomas are sporadic, a minority occur as components of familial syndromes, such as Multiple Endocrine Neoplasia type 1 (MEN1) and Carney complex (CNC). MEN1 syndrome is caused by germ-line mutations in the MEN1 gene, whereas most of the CNC patients carry the mutated protein kinase A (PKA) regulatory subunit-1-α (PRKAR1A) gene. Recently, other conditions predisposing to endocrine tumors have been identified: Pituitary Adenoma Predisposition (PAP) and MEN type 4 (MEN4). PAP was originally identified in a genetically homogeneous Finnish population. In a population based cohort from Northern Finland, aryl hydrocarbon receptor-interacting protein (AIP) gene mutations were found in 16% of all patients diagnosed with growth hormone (GH) producing pituitary adenoma, and in 40% of the subset of patients who were diagnosed under the age of 35 years. Since AIP mutations were originally described in a defined, homogeneous population from Northern Finland, it was relevant to study whether mutations also occur in more heterogeneous populations. In patient cohorts with different ethnic origins and variable clinical phenotypes, germ-line AIP mutations were detectable at low frequencies (range 0.8-7.4%). AIP mutation-positive patients were often diagnosed with a GH-producing adenoma at a young age, and usually had no family history of endocrine tumors. The low frequency of AIP mutations in randomly selected patients, and the lack of any family history of pituitary adenomas create a challenge for the identification of PAP patients. Our preliminary study suggests that AIP immunohistochemistry may serve as a pre-screening tool to distinguish between the AIP mutation-negative and the mutation-positive tumors. Tumors of various endocrine glands are components of MEN1 and CNC syndromes. Somatic MEN1 and PRKAR1A mutations in sporadic pituitary adenomas are rare, but occur in some of the other tumors related to these syndromes. The role of AIP mutations in endocrine neoplasia was studied and our results indicated that somatic AIP mutations are rare or non-existent in sporadic tumors of endocrine glands (0 of 111). Furthermore, germ-line AIP mutations in prolactin producing adenomas (2 of 9) confirmed the role of this pituitary tumor type in the PAP phenotype. Thyroid disorders are common in the general population, and the majority of them are sporadic. Interestingly, it has been suggested that thyroid disorders might be more common in PAP families. For this reason we studied germ-line AIP mutations in 93 index cases from familial non-medullary thyroid cancer (NMTC) families. The underlying gene or genes for familial NMTC have not been identified yet. None of the patients had any potentially pathogenic AIP mutation. This suggests that AIP is unlikely to play a role in familial NMTCs. A novel multiple endocrine syndrome was originally described in rats with phenotypic features of human MEN type 1 and 2. Germ-line mutations of cyclin-dependent kinase inhibitor 1B (CDKN1B also known as p27Kip1) gene were reported later in these rats and a germ-line mutation was also identified in one human family with MEN1-like phenotype (later named MEN4). To confirm the importance of this gene’s mutations in humans, we performed a mutation screening in MEN-like patients and in patients with pituitary adenoma. Our results indicate that CDKN1B/p27Kip1 mutations appear in a small portion of MEN1-like patients (one of 36), and that such mutations are rare or non-existent in both familial (0 of 19) and sporadic pituitary adenoma patients (0 of 50). In conclusion, this work strengthens the tumor susceptibility role of AIP and CDKN1B/p27Kip1 in endocrine neoplasia. Clarifying the PAP phenotype facilitates the identification of potential AIP mutation carriers. Genetic counseling can be offered to the relatives and follow-up of the mutation carriers can be organized, hence an earlier diagnosis is feasible.

Amyloid diseases at old age : a pathological, epidemiological, and genetic study

Along with the increased life span of individuals, the burden of old age-associated diseases has inevitably increased. Alzheimer s disease (AD), probably the most well known geriatric disease, belongs to the old age-associated amyloid diseases. The purpose of this study was to investigate the frequency, genetic and health-associated risk factors, mutual association, and amyloid proteins in two old age-associated amyloid disorders senile systemic amyloidosis (SSA) and cerebral amyloid angiopathy (CAA) as part of the prospective population-based Vantaa 85+ autopsy study on a Finnish population aged 85 years or more (Studies I-III), completed with a case report on a patient with advanced AGel amyloidosis (Study IV). The numbers of patients investigated in the studies (I-III) were 256, 74, and 63, respectively. The diagnosis and grading of amyloid were based upon histological examination of tissue samples obtained post mortem and stained with Congo red. The amyloid fibril and associated proteins were characterized by immunohistochemical staining methods. The genotype frequencies of 20 polymorphisms in 9 genes and information on health-associated risk factors in subjects with and without SSA and CAA were compared. In a Finnish population ≥ 95 years of age, SSA and CAA occurred in 36% and 49% of the subjects, respectively. In total, two-thirds of these very elderly individuals had SSA, CAA, or both. However, in only 14% of the population these two conditions co-occurred. In subjects 85 years or older, the prevalence of SSA was 25%. In this population, SSA was associated with age at the time of death (p=0.002), myocardial infarctions (MIs; p=0.004), the G/G (Val/Val) genotype of the exon 24 polymorphism in the alpha2-macroglobulin (α2M) gene (p=0.042) and with the H2 haplotype of the tau gene (p=0.016). In contrast, the presence of CAA was strongly associated with APOE e4 (p=0.0003), with histopathological AD (p=0.0005), and with clinical dementia (p=0.01) in both e4+ (p=0.02) and e4- (p=0.06) individuals. Apart from demonstrating the amyloid fibril proteins, complement proteins 3d (C3d) and 9 (C9) were detected in the amyloid deposits of CAA and AGel amyloidosis, and α2M protein was found in fibrous scar tissue close to SSA. In conclusion, this first population based study on SSA shows that both SSA and CAA are common in very elderly individuals. Old age, MIs, the exon 24 polymorphism of the α2M gene, and H1/H2 polymorphism of the tau gene associate with SSA while clinical dementia and APOE ε4 genotype associate with CAA. The high prevalence of CAA, combined with its association with clinical dementia independent of APOE genotype, neuropathological AD, or SSA, also highlights its clinical significance in the very aged, among which the serious end stage complications of CAA, namely multiple infarctions and hemorrhages, are rare. The report on a patient having advanced AGel amyloidosis added knowledge on the disease and showed that this generally benign condition occasionally may lead to death. Further studies are warranted to confirm the findings in other populations. Also, the role of α2M and tau in the pathogenesis of SSA and the involvement of complement in the process of amyloid beta (Aβ) protein elimination from the brain remain to be clarified. Finally, the high prevalence of SSA in the elderly raises the need for prospective clinical studies to define its clinical significance.

Kohti hyvää suometsien hoitoa : harvennusten ja kunnostusojitusten vaikutus ojitusaluemänniköiden puuntuotokseen ja metsänkasvatuksen taloustulokseen

This thesis examines the impacts of silvicultural activities on productivity and financial returns of Scots pine (Pinus sylvestris L.) stands on drained peatlands in Finland. The effects of ditch network maintenance operations (DNM) and thinnings, with different timings and intensities, were studied. Based on stand development simulations, the best regimes for different types of stands according to site type, climatic area, and stand silvicultural status were defined from the viewpoint of both wood production and financial profitability. Certain aspects affecting the management outcomes, such as the timing of the first thinning, were examined using data from thinning experiments. Long-term predictions of the impacts of different management regimes were carried out by simulating the development of well-representative model-stands which were composed from appropriate inventory data sets. The MOTTI stand simulator used to perform the simulations enables the predictions by utilizing specific models for drained peatland stands. In addition to natural stand dynamics, these models describe the effects of silvicultural treatments on the development of a given stand. The mean annual increment of merchantable wood (MAImerch) was used as the measure of wood productivity, and the financial feasibility of the regimes was compared using net present value (NPV) analysis. Silvicultural treatments, when applied to appropriately match stand condition, increased both the productivity and financial returns of stand management. Applying DNM resulted in a small increase in MAImerch. When thinning was introduced along with DNM, their combined effect on wood productivity was considerable. According to current operational practices, DNM is generally combined with thinning. In some cases, e.g., in sites of low productivity, the need for DNM may become apparent prior to the thinning stage. As for profitability, thinnings proved to be crucial. The regimes with heavy and late thinnings were generally more profitable than those with normal thinnings. Further, early thinning (relative to stand volume) lacked appeal when seeking a financially profitable removal from the first thinning. In young stands with an initially poor silvicultural condition, however, applying even a low-yielding first thinning considerably increased the NPV when compared to a regime with no thinning at all. Generally, the regimes resulting in the best profitability included heavier thinnings and fewer DNM and thinning treatments than did the regimes resulting in the best yield results. This study demonstrates considerable potential for profitable wood production-oriented management in pine stands on drained peatlands despite their challenging circumstances and long rotations. The results can be used for defining new and more site-specific silvicultural guidelines for various types of drained, pine-dominated peatland stands within the entire range of boreal conditions.

Design and testing of stand-specific bucking instructions for use on modern cut-to-length harvesters

This study addresses three important issues in tree bucking optimization in the context of cut-to-length harvesting. (1) Would the fit between the log demand and log output distributions be better if the price and/or demand matrices controlling the bucking decisions on modern cut-to-length harvesters were adjusted to the unique conditions of each individual stand? (2) In what ways can we generate stand and product specific price and demand matrices? (3) What alternatives do we have to measure the fit between the log demand and log output distributions, and what would be an ideal goodness-of-fit measure? Three iterative search systems were developed for seeking stand-specific price and demand matrix sets: (1) A fuzzy logic control system for calibrating the price matrix of one log product for one stand at a time (the stand-level one-product approach); (2) a genetic algorithm system for adjusting the price matrices of one log product in parallel for several stands (the forest-level one-product approach); and (3) a genetic algorithm system for dividing the overall demand matrix of each of the several log products into stand-specific sub-demands simultaneously for several stands and products (the forest-level multi-product approach). The stem material used for testing the performance of the stand-specific price and demand matrices against that of the reference matrices was comprised of 9 155 Norway spruce (Picea abies (L.) Karst.) sawlog stems gathered by harvesters from 15 mature spruce-dominated stands in southern Finland. The reference price and demand matrices were either direct copies or slightly modified versions of those used by two Finnish sawmilling companies. Two types of stand-specific bucking matrices were compiled for each log product. One was from the harvester-collected stem profiles and the other was from the pre-harvest inventory data. Four goodness-of-fit measures were analyzed for their appropriateness in determining the similarity between the log demand and log output distributions: (1) the apportionment degree (index), (2) the chi-square statistic, (3) Laspeyres quantity index, and (4) the price-weighted apportionment degree. The study confirmed that any improvement in the fit between the log demand and log output distributions can only be realized at the expense of log volumes produced. Stand-level pre-control of price matrices was found to be advantageous, provided the control is done with perfect stem data. Forest-level pre-control of price matrices resulted in no improvement in the cumulative apportionment degree. Cutting stands under the control of stand-specific demand matrices yielded a better total fit between the demand and output matrices at the forest level than was obtained by cutting each stand with non-stand-specific reference matrices. The theoretical and experimental analyses suggest that none of the three alternative goodness-of-fit measures clearly outperforms the traditional apportionment degree measure. Keywords: harvesting, tree bucking optimization, simulation, fuzzy control, genetic algorithms, goodness-of-fit

Soil hydrological properties and conditions, site preparation, and the long-term performance of planted Scots pine (Pinus sylvestris L.) on upland forest sites in Finnish Lapland

Scots pine (Pinus sylvestris L.) and Norway spruce (Picea abies (L.) Karst.) forests dominate in Finnish Lapland. The need to study the effect of both soil factors and site preparation on the performance of planted Scots pine has increased due to the problems encountered in reforestation, especially on mesic and moist, formerly spruce-dominated sites. The present thesis examines soil hydrological properties and conditions, and effect of site preparation on them on 10 pine- and 10 spruce-dominated upland forest sites. Finally, the effects of both the site preparation and reforestation methods, and soil hydrology on the long-term performance of planted Scots pine are summarized. The results showed that pine and spruce sites differ significantly in their soil physical properties. Under field capacity or wetter soil moisture conditions, planted pines presumably suffer from excessive soil water and poor soil aeration on most of the originally spruce sites, but not on the pine sites. The results also suggested that site preparation affects the soil-water regime and thus prerequisites for forest growth over two decades after site preparation. High variation in the survival and mean height of planted pine was found. The study suggested that on spruce sites, pine survival is the lowest on sites that dry out slowly after rainfall events, and that height growth is the fastest on soils that reach favourable aeration conditions for root growth soon after saturation, and/or where the average air-filled porosity near field capacity is large enough for good root growth. Survival, but not mean height can be enhanced by employing intensive site preparation methods on spruce sites. On coarser-textured pine sites, site preparation methods don t affect survival, but methods affecting soil fertility, such as prescribed burning and ploughing, seem to enhance the height growth of planted Scots pines over several decades. The use of soil water content in situ as the sole criterion for sites suitable for pine reforestation was tested and found to be a relatively uncertain parameter. The thesis identified new potential soil variables, which should be tested using other data in the future.

Digitoitujen väri-infrailmakuvien numeerinen monikuvatulkinta metsäninventoinnissa

The use of remote sensing imagery as auxiliary data in forest inventory is based on the correlation between features extracted from the images and the ground truth. The bidirectional reflectance and radial displacement cause variation in image features located in different segments of the image but forest characteristics remaining the same. The variation has so far been diminished by different radiometric corrections. In this study the use of sun azimuth based converted image co-ordinates was examined to supplement auxiliary data extracted from digitised aerial photographs. The method was considered as an alternative for radiometric corrections. Additionally, the usefulness of multi-image interpretation of digitised aerial photographs in regression estimation of forest characteristics was studied. The state owned study area located in Leivonmäki, Central Finland and the study material consisted of five digitised and ortho-rectified colour-infrared (CIR) aerial photographs and field measurements of 388 plots, out of which 194 were relascope (Bitterlich) plots and 194 were concentric circular plots. Both the image data and the field measurements were from the year 1999. When examining the effect of the location of the image point on pixel values and texture features of Finnish forest plots in digitised CIR photographs the clearest differences were found between front-and back-lighted image halves. Inside the image half the differences between different blocks were clearly bigger on the front-lighted half than on the back-lighted half. The strength of the phenomenon varied by forest category. The differences between pixel values extracted from different image blocks were greatest in developed and mature stands and smallest in young stands. The differences between texture features were greatest in developing stands and smallest in young and mature stands. The logarithm of timber volume per hectare and the angular transformation of the proportion of broadleaved trees of the total volume were used as dependent variables in regression models. Five different converted image co-ordinates based trend surfaces were used in models in order to diminish the effect of the bidirectional reflectance. The reference model of total volume, in which the location of the image point had been ignored, resulted in RMSE of 1,268 calculated from test material. The best of the trend surfaces was the complete third order surface, which resulted in RMSE of 1,107. The reference model of the proportion of broadleaved trees resulted in RMSE of 0,4292 and the second order trend surface was the best, resulting in RMSE of 0,4270. The trend surface method is applicable, but it has to be applied by forest category and by variable. The usefulness of multi-image interpretation of digitised aerial photographs was studied by building comparable regression models using either the front-lighted image features, back-lighted image features or both. The two-image model turned out to be slightly better than the one-image models in total volume estimation. The best one-image model resulted in RMSE of 1,098 and the two-image model resulted in RMSE of 1,090. The homologous features did not improve the models of the proportion of broadleaved trees. The overall result gives motivation for further research of multi-image interpretation. The focus may be improving regression estimation and feature selection or examination of stratification used in two-phase sampling inventory techniques. Keywords: forest inventory, digitised aerial photograph, bidirectional reflectance, converted image co­ordinates, regression estimation, multi-image interpretation, pixel value, texture, trend surface

Energiapuuharvennus osana metsänkasvatusta ja ilmastonmuutoksen hillitsemistä

The aim of this study was to compare the differences between forest management incorporating energy wood thinning and forest management based on silvicultural recommendations (baseline). Energy wood thinning was substituted for young stand thinning and the first commercial thinning of industrial wood. The study was based on the forest stand data from Southern Finland, which were simulated by the MOTTI-simulator. The main interest was to find out the climatic benefits resulting from carbon sequestration and energy substitution. The value of energy wood was set to substitute it for coal as an alternative energy fuel (emission trade). Other political instruments (Kemera subsidies) were also analysed. The largest carbon dioxide emission reductions were achieved as a combination of carbon sequestration and energy substitution (on average, a 26-90 % increase in discounted present value in the beginning of rotation) compared to the baseline. Energy substitution increased emission reductions more effectively than carbon sequestration, when maintaining dense young stands. According to the study, energy wood thinning as a part of forest management was more profitable than the baseline when the value of carbon dioxide averaged more than 15 €/CO2 and other political subsidies were unchanged. Alternatively, the price of energy wood should on average exceed 21 €/m3 on the roadside in order to be profitable in the absence of political instruments. The most cost-efficient employment of energy wood thinning occured when the dominant height was 12 meters, when energy substitution was taken into account. According to alternative forest management, thinning of sapling stands could be done earlier or less intensely than thinning based on silvicultural recommendations and the present criteria of subsidies. Consequently, the first commercial thinning could be profitable to carry out either as harvesting of industrial wood or energy wood, or as integrated harvesting depending on the costs of the harvesting methods available and the price level of small-size industrial wood compared to energy wood.

Genetic diversity of Atlantic salmon (Salmo salar) and European whitefish (Coregonus lavaretus) in the Baltic Sea

We described the patterns and extent of microsatellite DNA variation in historical and present-day Atlantic salmon (Salmo salar L.) stocks in the Baltic Sea and neighbouring areas, and in European whitefish (Coregonus lavaretus) ecotypes, populations and run-timing types in Finland. Moreover, the amount and pattern of genetic diversity in historical salmon populations before human impact were described, and the proportion of diversity maintained in the present hatchery stocks evaluated. Salmon populations in the Baltic Sea were, on average, significantly less variable than eastern Atlantic populations, and the diversity of landlocked populations (Lakes Vänern, Saimaa, Onega and Ladoga) was in turn significantly lower than that of anadromous salmon populations in the Baltic Sea populations. Within the Baltic Sea, the anadromous populations of Atlantic salmon formed three clear groups, corresponding to the northern (Gulf of Bothnia), eastern (Gulf of Finland and eastern Baltic Main Basin) and southern (western Baltic Main Basin) regions. Based on microsatellite data, three salmon population groups in the Baltic Sea were considered potentially different colonization lineages. In short- and long-term breeding programmes of Atlantic salmon, the average observed rate of loss of alleles was 4.9% and 2.0% per generation and the average rate of loss of heterozygosity was 1.4% and 1% per generation, respectively. When comparing the genetic parameters of stocks before and after hatchery breeding of several successive generations (Rivers Iijoki and Oulujoki), statistically significant changes in allele frequencies were common, while large wild stock in the Teno River has remained temporally very stable over 56 years. Despite the observed losses of genetic diversity in broodstock breeding, a large proportion of the genetic resources of the extirpated stocks are still conserved in the broodstocks. Genetic differentiation among European whitefish ecotypes was generally low, thus giving support to the hypothesis of one native European whitefish species in Fennoscandia. Among the ecotypes, the northern, large sparsely rakered, bottom-dwelling whitefish was the most unique. The known genetic differences in quantitative traits have thus either developed independently of potential phylogenetic lineages, or the lineages have mixed and the quantitative traits of the ecotypes, like gill-raker number, have later changed according to environment and selection pressures. Overall, genetic distances between the anadromous whitefish populations along the Finnish coast, especially in the Bothnian Bay area, were small. Wild whitefish populations studied had slightly higher allelic diversity than hatchery-reared populations in corresponding rivers.