Genotyping for genetic association studies: Methods and applications

In this thesis, two separate single nucleotide polymorphism (SNP) genotyping techniques were set up at the Finnish Genome Center, pooled genotyping was evaluated as a screening method for large-scale association studies, and finally, the former approaches were used to identify genetic factors predisposing to two distinct complex diseases by utilizing large epidemiological cohorts and also taking environmental factors into account. The first genotyping platform was based on traditional but improved restriction-fragment-length-polymorphism (RFLP) utilizing 384-microtiter well plates, multiplexing, small reaction volumes (5 µl), and automated genotype calling. We participated in the development of the second genotyping method, based on single nucleotide primer extension (SNuPeTM by Amersham Biosciences), by carrying out the alpha- and beta tests for the chemistry and the allele-calling software. Both techniques proved to be accurate, reliable, and suitable for projects with thousands of samples and tens of markers. Pooled genotyping (genotyping of pooled instead of individual DNA samples) was evaluated with Sequenom s MassArray MALDI-TOF, in addition to SNuPeTM and PCR-RFLP techniques. We used MassArray mainly as a point of comparison, because it is known to be well suited for pooled genotyping. All three methods were shown to be accurate, the standard deviations between measurements being 0.017 for the MassArray, 0.022 for the PCR-RFLP, and 0.026 for the SNuPeTM. The largest source of error in the process of pooled genotyping was shown to be the volumetric error, i.e., the preparation of pools. We also demonstrated that it would have been possible to narrow down the genetic locus underlying congenital chloride diarrhea (CLD), an autosomal recessive disorder, by using the pooling technique instead of genotyping individual samples. Although the approach seems to be well suited for traditional case-control studies, it is difficult to apply if any kind of stratification based on environmental factors is needed. Therefore we chose to continue with individual genotyping in the following association studies. Samples in the two separate large epidemiological cohorts were genotyped with the PCR-RFLP and SNuPeTM techniques. The first of these association studies concerned various pregnancy complications among 100,000 consecutive pregnancies in Finland, of which we genotyped 2292 patients and controls, in addition to a population sample of 644 blood donors, with 7 polymorphisms in the potentially thrombotic genes. In this thesis, the analysis of a sub-study of pregnancy-related venous thromboses was included. We showed that the impact of factor V Leiden polymorphism on pregnancy-related venous thrombosis, but not the other tested polymorphisms, was fairly large (odds ratio 11.6; 95% CI 3.6-33.6), and increased multiplicatively when combined with other risk factors such as obesity or advanced age. Owing to our study design, we were also able to estimate the risks at the population level. The second epidemiological cohort was the Helsinki Birth Cohort of men and women who were born during 1924-1933 in Helsinki. The aim was to identify genetic factors that might modify the well known link between small birth size and adult metabolic diseases, such as type 2 diabetes and impaired glucose tolerance. Among ~500 individuals with detailed birth measurements and current metabolic profile, we found that an insertion/deletion polymorphism of the angiotensin converting enzyme (ACE) gene was associated with the duration of gestation, and weight and length at birth. Interestingly, the ACE insertion allele was also associated with higher indices of insulin secretion (p=0.0004) in adult life, but only among individuals who were born small (those among the lowest third of birth weight). Likewise, low birth weight was associated with higher indices of insulin secretion (p=0.003), but only among carriers of the ACE insertion allele. The association with birth measurements was also found with a common haplotype of the glucocorticoid receptor (GR) gene. Furthermore, the association between short length at birth and adult impaired glucose tolerance was confined to carriers of this haplotype (p=0.007). These associations exemplify the interaction between environmental factors and genotype, which, possibly due to altered gene expression, predisposes to complex metabolic diseases. Indeed, we showed that the common GR gene haplotype associated with reduced mRNA expression in thymus of three individuals (p=0.0002).

Molecular Genetics of Schizophrenia and Related Intermediate Phenotypes in a Founder Population

Schizophrenia, affecting about 1% of population worldwide, is a severe mental disorder characterized by positive and negative symptoms, such as psychosis and anhedonia, as well as cognitive deficits. At present, schizophrenia is considered a complex disorder of neurodevelopmental origin with both genetic and environmental factors contributing to its onset. Although a number of candidate genes for schizophrenia have been highlighted, only very few schizophrenia patients are likely to share identical genetic liability. This study is based on the nation-wide schizophrenia family sample of the National Institute for Health and Welfare, and represents one of the largest and most well-characterized familial series in the world. In the first part of this study, we investigated the roles of the DTNBP1, NRG1, and AKT1 genes in the background of schizophrenia in Finland. Although these genes are associated with schizophrenia liability in several populations, any significant association with clinical diagnostic information of schizophrenia remained absent in our sample of 441 schizophrenia families. In the second part of this study, we first replicated schizophrenia linkage on the long arm of chromosome 7 in 352 schizophrenia families. In the following association analysis, we utilized additional clinical disorder features and intermediate phenotypes – endophenotypes - in addition to diagnostic information from altogether 290 neuropsychologically assessed schizophrenia families. An intragenic short tandem repeat allele of the regional RELN gene, supposed to play a role in the background of several neurodevelopmental disorders, showed significant association with poorer cognitive functioning and more severe schizophrenia symptoms. Additionally, this risk allele was significantly more prevalent among the individuals affected with schizophrenia spectrum disorders. We have previously identified linkage of schizophrenia and its cognitive endophenotypes on the long arms of chromosomes 2, 4, and 5. In the last part of this study, we selected altogether 104 functionally relevant candidate genes from the linked regions. We detected several promising associations, of which especially interesting are the ERBB4 gene, showing association with the severity of schizophrenia symptoms and impairments in traits related to verbal abilities, and the GRIA1 gene, showing association with the severity of schizophrenia symptoms. Our results extend the previous evidence that the genetic risk for schizophrenia is at least partially mediated via the effects of the candidate genes and their combinations on relevant brain systems, resulting in alterations in different disorder domains, such as the cognitive deficits.

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.

Molecular Genetic Background of Tumours in Hereditary Leiomyomatosis and Renal Cell Cancer Syndrome

Hereditary Leiomyomatosis and Renal Cell Cancer (HLRCC) is a hereditary tumour predisposition syndrome. Its phenotype includes benign cutaneous and uterine leiomyomas (CLM, ULM) with high penetrance and rarer renal cell cancer (RCC), most commonly of papillary type 2 subtype. Over 130 HLRCC families have been identified world-wide but the RCC phenotype seems to concentrate in families from Finland and North America for unknown reasons. HLRCC is caused by heterozygous germline mutations in the fumarate hydratase (FH) gene. FH encodes the enzyme fumarase from mitochondrial citric acid cycle. Fumarase enzyme activity or type or site of the FH mutation are unassociated with disease phenotype. The strongest evidence for tumourigenesis mechanism in HLRCC supports a hypoxia inducible factor driven process called pseudohypoxia resulting from accumulation of the fumarase substrate fumarate. In this study, to assess the importance of gene- or exon-level deletions or amplifications of FH in patients with HLRCC-associated phenotypes, multiplex ligation-dependent probe amplification (MLPA) method was used. One novel FH mutation, deletion of exon 1, was found in a Swedish male patient with an evident HLRCC phenotype with CLM, RCC, and a family history of ULM and RCC. Six other patients with CLM and 12 patients with only RCC or uterine leiomyosarcoma (ULMS) remained FH mutation-negative. These results suggest that copy number aberrations of FH or its exons are an infrequent cause of HLRCC and that only co-occurrence of benign tumour types justifies FH-mutation screening in RCC or ULMS patients. Determination of the genomic profile of 11 HLRCC-associated RCCs from Finnish patients was performed by array comparative genomic hybridization. The most common copy number aberrations were gains of 2, 7, and 17 and losses of 13q12.3-q21.1, 14, 18, and X. When compared to aberrations of sporadic papillary RCCs, HLRCC-associated RCCs harboured a distinct DNA copy number profile and lacked many of the changes characterizing the sporadic RCCs. The findings suggest a divergent molecular pathway for tumourigenesis of papillary RCCs in HLRCC. In order to find a genetic modifier of RCC risk in HLRCC, genome-wide linkage and identical by descent (IBD) analysis studies were performed in Finnish HLRCC families with microsatellite marker mapping and SNP-array platforms. The linkage analysis identified only one locus of interest, the FH gene locus in 1q43, but no mutations were found in the genes of the region. IBD analysis yielded no convincing haplotypes shared by RCC patients. Although these results do not exclude the existence of a genetic modifier for RCC risk in HLRCC, they emphasize the role of FH mutations in the malignant tumourigenesis of HLRCC. To study the benign tumours in HLRCC, genome-wide DNA copy number and gene expression profiles of sporadic and HLRCC ULMs were defined with modern SNP- and gene-expression array platforms. The gene expression array suggests novel genes involved in FH-deficient ULM tumourigenesis and novel genes with putative roles in propagation of sporadic ULM. Both the gene expression and copy number profiles of HLRCC ULMs differed from those of sporadic ULMs indicating distinct molecular basis of the FH-deficient HLRCC tumours.

Genetic studies on recurrent miscarriage

Recurrent miscarriage (RM) is defined as three consecutive pregnancy failures and is estimated to affect ~1% of couples trying to conceive. The cause of RM remains unknown in approximately 50% of cases. In this study, it was hypothesized that some of the underlying factors yet to be discovered are genetic. The aim was to search for mutations in genes AMN, EPCR, TM, and p53 known to cause miscarriage in mouse models and thereby find new genetic causes for unexplained miscarriages in humans. In addition, the mitochondrial genome was studied because mitochondria are involved in processes important in early development. Furthermore, sex chromosome characteristics suggested to underlie miscarriage were also studied. A total of 40 couples and 8 women with unexplained RM were collected for this study and screened for mutations in the candidate genes. Six interesting exonic or potential splice site disrupting variations were detected. However, their phenotypic effects cannot be determined without further investigations. Additionally, an association between the C11992A polymorphism of the p53 gene and RM was detected. The results indicate that women carrying the C/A or A/A genotype have a two-fold higher risk for RM than women with a C/C genotype. This strengthens the results of previous studies reporting that p53 sequence variations may cause miscarriage. The role of variation C11992A in embryonic development is, however, difficult to predict without further studies When screening the mitochondrial genome a heteroplasmic mtDNA variation was found in an unexpected high number of women, as heteroplasmic variations are reported to be rare. One novel variation and 18 previously reported polymorphisms were detected in the mitochondrial genome. Although the detected variations are likely to be neutral polymorphisms, a role in the aetiology of miscarriage cannot be excluded as some mtDNA variations may be pathogenic only when a threshold is reached. Recent publications have reported skewed X chromosome inactivation and Y chromosome microdeletions to be associated with RM. Therefore, these sex chromosome abnormalities in the context of RM were investigated. No associations between skewed X chromosome inactivation or Y chromosome microdeletions and RM in the Finnish patients were detected. Data on ancestral birthplaces of the patients were collected to study any possible geographic clustering, which would indicate a common predisposing factor. The results showed clustering of the birthplaces in eastern Finland in a subset of patients. This suggests a possibility of an enriched susceptibility gene which may contribute to RM.

Search For Susceptibility Genes In Schizophrenia

Schizophrenia is a severe mental disorder affecting 0.4-1% of the population worldwide. It is characterized by impairments in the perception of reality and by significant social or occupational dysfunction. The disorder is one of the major contributors to the global burden of diseases. Studies of twins, families, and adopted children point to strong genetic components for schizophrenia, but environmental factors also play a role in the pathogenesis of disease. Molecular genetic studies have identified several potential positional candidate genes. The strongest evidence for putative schizophrenia susceptibility loci relates to the genes encoding dysbindin (DTNBP1) and neuregulin (NRG1), but studies lack impressive consistency in the precise genetic regions and alleles implicated. We have studied the role of three potential candidate genes by genotyping 28 single nucleotide polymorphisms in the DNTBP1, NRG1, and AKT1 genes in a large schizophrenia family sample consisting of 441 families with 865 affected individuals from Finland. Our results do not support a major role for these genes in the pathogenesis of schizophrenia in Finland. We have previously identified a region on chromosome 5q21-34 as a susceptibility locus for schizophrenia in a Finnish family sample. Recently, two studies reported association between the γ-aminobutyric acid type A receptor cluster of genes in this region and one study showed suggestive evidence for association with another regional gene encoding clathrin interactor 1 (CLINT1, also called Epsin 4 and ENTH). To further address the significance of these genes under the linkage peak in the Finnish families, we genotyped SNPs of these genes, and observed statistically significant association of variants between GABRG2 and schizophrenia. Furthermore, these variants also seem to affect the functioning of the working memory. Fetal events and obstetric complications are associated with schizophrenia. Rh incompatibility has been implicated as a risk factor for schizophrenia in several epidemiological studies. We conducted a family-based candidate-gene study that assessed the role of maternal-fetal genotype incompatibility at the RhD locus in schizophrenia. There was significant evidence for an RhD maternal-fetal genotype incompatibility, and the risk ratio was estimated at 2.3. This is the first candidate-gene study to explicitly test for and provide evidence of a maternal-fetal genotype incompatibility mechanism in schizophrenia. In conclusion, in this thesis we found evidence that one GABA receptor subunit, GABRG2, is significantly associated with schizophrenia. Furthermore, it also seems to affect to the functioning of the working memory. In addition, an RhD maternal-fetal genotype incompatibility increases the risk of schizophrenia by two-fold.

Genetic alterations in microsatellite-unstable colorectal cancer

Colorectal cancer (CRC) is the third most common cancer in Finland. Of all CRC tumors, 15% display microsatellite-instability (MSI) caused by defective cellular mismatch repair. Cells displaying MSI accumulate a high number of mutations genome-wide, especially in short repeat areas, microsatellites. When targeting genes essential for cell growth or death, MSI can promote tumorigenesis. In non-coding areas, microsatellite mutations are generally considered as passenger events. Since the discovery of MSI and its linkage to cancer, more that 200 genes have been investigated for a role in MSI tumorigenesis. Although various criteria have been suggested for MSI target gene identification, the challenge has been to distinguish driver mutations from passenger mutations. This study aimed to clarify these key issues in the research field of MSI cancer. Prior to this, background mutation rate in MSI cancer has not been studied in a large-scale. We investigated the background mutation rate in MSI CRC by analyzing the spectrum of microsatellite mutations in non-coding areas. First, semenogelin I was studied for a possible role in MSI carcinogenesis. The intronic T9 repeat of semenogelin I was frequently mutated but no evidence for selection during tumorigenesis was obtained. Second, a sequencing approach was utilized to evaluate the general background mutation rate in MSI CRC. Both intronic and intergenic repeats harbored extremely high mutation rates of ≤ 87% and intergenic repeats were more unstable than the intronic repeats. As mutation rates of presumably neutral microsatellites can be high in MSI CRC in the absence of apparent selection pressure, high mutation frequency alone is not sufficient evidence for identification of driver MSI target genes. Next, an unbiased approach was designed to identify the mutatome of MSI CRC. By combining expression array data and a database search we identified novel genes possibly related to MSI CRC carcinogenesis. One of the genes was studied further. In the functional analysis this gene was observed to cause an abnormal cancer-prone cellular phenotype, possibly through altered responses to DNA damage. In our recent study, smooth muscle myosin heavy chain 11 (MYH11) was identified as a novel MSI CRC gene. Additionally, MYH11 has a well established role in acute myeloid leukemia (AML) through an oncogenic fusion protein CBFB-MYH11. We investigated further the role of MYH11 in AML by sequencing. Three novel missense variants of MYH11 were identified. None of the variants were present in the population-based control material. One of the identified variants, V71A, lies in the N-terminal SH3-like domain of MYH11 of unknown function. The other two variants, K1059E and R1792Q are located in the coil-coiled myosin rod essential for the regulation and filament formation of MYH11. The variant K1059E lies in the close proximity of the K1044N that has been functionally assessed in our earlier work of CRC and has been reported to cause total loss of MYH11 protein regulation. As the functional significance of the three novel variants examined in this work remains unknown, future studies should clarify the further role of MYH11 in AML leukaemogenesis and in other malignancies.

Human Torque teno virus: epidemiology, cell biology and immunology

Torque teno virus (TTV) was discovered in 1997 in the serum of a Japanese patient who had a post-transfusion hepatitis of unknown etiology. It is a small virus containing a circular single-stranded DNA genome which is unique among human viruses. Within a few years after its discovery, the TTVs were noted to form a large family of viruses with numerous genotypes. TTV is highly prevalent among the general population throughout the world, and persistent infections and co-infections with several genotypes occur frequently. However, the pathogenicity and the mechanism for the sustained occurrence of the virus in blood are at present unclear. To determine the prevalence of TTV in Finland, we set up PCR methods and examined the sera of asymptomatic subjects for the presence of TTV DNA and for genotype-6 DNA. TTV was found to be highly prevalent also in Finland; 85% of adults harbored TTV in their blood, and 4% were infected with genotype-6. In addition, TTV DNA was detected in a number of different tissues, with no tissue-type or symptom specificity. Most cell-biological events during TTV infections are at the moment unknown. Replicating TTV DNA has, however, been detected in liver and the hematopoietic compartment, and three mRNAs are known to be generated. To characterize TTV cell biology in more detail, we cloned in full length the genome of TTV genotype 6. We showed that in human kidney-derived cells TTV produces altogether six proteins with distinct subcellular localizations. TTV mRNA transcription was detected in all cell lines transfected with the full-length clone, and TTV DNA replicated in several of them, including those of erythroid, kidney, and hepatic origin. Furthermore, the viral DNA replication was shown to utilize the cellular DNA polymerases. Diagnoses of TTV infections have been based almost solely on PCR, whereas serological tests, measuring antibody responses, would give more information on many aspects of these infections. To investigate the TTV immunology in more detail, we produced all six TTV proteins for use as antigens in serological tests. We detected in human sera IgM and IgG antibodies to occur simultaneously with TTV DNA, and observed appearance of TTV DNA regardless of pre-existing antibodies, and disappearance of TTV DNA after antibody appearance. The genotype-6 nucleotide sequence remained stable for years within the infected subjects, suggesting that some mechanism other than mutations is used by this minute virus to evade our immune system and to establish chronic infections in immunocompetent subjects.

Genetic Heterogeneity in Autism Spectrum Disorders in a Population Isolate

Positional cloning has enabled hypothesis-free, genome-wide scans for genetic factors contributing to disorders or traits. Traditionally linkage analysis has been used to identify regions of interest, followed by meticulous fine mapping and candidate gene screening using association methods and finally sequencing of regions of interest. More recently, genome-wide association analysis has enabled a more direct approach to identify specific genetic variants explaining a part of the variance of the phenotype of interest. Autism spectrum disorders (ASDs) are a group of childhood onset neuropsychiatric disorders with shared core symptoms but varying severity. Although a strong genetic component has been established in ASDs, genetic susceptibility factors have largely eluded characterization. Here, we have utilized modern molecular genetic methods combined with the advantages provided by the special population structure in Finland to identify genetic risk factors for ASDs. The results of this study show that numerous genetic risk factors exist for ASDs even within a population isolate. Stratification based on clinical phenotype resulted in encouraging results, as previously identified linkage to 3p14-p24 was replicated in an independent family set of families with Asperger syndrome, but no other ASDs. Fine-mapping of the previously identified linkage peak for ASDs at 3q25-q27 revealed association between autism and a subunit of the 5-hydroxytryptamine receptor 3C (HTR3C). We also used dense, genome-wide single nucleotide polymorphism (SNP) data to characterize the population structure of Finns. We observed significant population substructure which correlates with the known history of multiple consecutive bottle-necks experienced by the Finnish population. We used this information to ascertain a genetically homogenous subset of autism families to identify possible rare, enriched risk variants using genome-wide SNP data. No rare enriched genetic risk factors were identified in this dataset, although a subset of families could be genealogically linked to form two extended pedigrees. The lack of founder mutations in this isolated population suggests that the majority of genetic risk factors are rare, de novo mutations unique to individual nuclear families. The results of this study are consistent with others in the field. The underlying genetic architecture for this group of disorders appears highly heterogeneous, with common variants accounting for only a subset of genetic risk. The majority of identified risk factors have turned out to be exceedingly rare, and only explain a subset of the genetic risk in the general population in spite of their high penetrance within individual families. The results of this study, together with other results obtained in this field, indicate that family specific linkage, homozygosity mapping and resequencing efforts are needed to identify these rare genetic risk factors.

Molecular genetics and evolution of Puumala hantavirus

Puumala virus (PUUV) is the causative agent of nephropathia epidemica (NE), a mild form of hemorrhagic fever with renal syndrome. Finland has the highest documented incidence of NE with around 1000 cases diagnosed annually. PUUV is also found in other Scandinavian countries, Central Europe and the European part of Russia. PUUV belongs to the genus Hantavirus in the family Bunyaviridae. Hantaviruses are rodent-borne viruses each carried by a specific host that is persistently and asymptomatically infected by the virus. PUUV is carried by the bank voles (Myodes glareolus, previously known as Clethrionomys glareolus). Hantaviruses have co-evolved with their carrier rodents for millions of years and these host animals are the evolutionary scene of hantaviruses. In this study, PUUV sequences were recovered from bank voles captured in Denmark and Russian Karelia to study the evolution of PUUV in Scandinavia. Phylogenetic analysis of these strains showed a geographical clustering of genetic variants following the presumable migration pattern of bank voles during the recolonization of Scandinavia after the last ice age approximately 10 000 years ago. The currently known PUUV genome sequences were subjected to in-depth phylogenetic analyses and the results showed that genetic drift seems to be the major mechanism of PUUV evolution. In general, PUUV seems to evolve quite slowly following a molecular clock. We also found evidence for recombination in the evolution of some genetic lineages of PUUV. Viral microevolution was studied in controlled virus transmission in colonized bank voles and changes in quasispecies dynamics were recorded as the virus was transmitted from one animal to another. We witnessed PUUV evolution in vivo, as one synonymous mutation became repeatedly fixed in the viral genome during the experiment. The detailed knowledge on the PUUV diversity was used to establish new sensitive and specific detection methods for this virus. Direct viral invasion of the hypophysis was demonstrated for the first time in a lethal case of NE. PUUV detection was done by immunohistochemistry, in situ hybridization and RT-nested-PCR of the autopsy tissue samples.

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.

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.

Nuorten männiköiden kehitys - taimikonhoito ja ensiharvennus

The aim of this study was to investigate the effects of location, site type, regeneration method and precommercial thinning on the characteristics and development of young, even-aged, pure Scots pine stands. In addition, the effects of timing and intensity of first commercial thinning on the yield and profitability during the rotation period were also studied. The stand characteristics and external quality of young Scots pine stands and stand-level growth models were based on extensive inventory data of the Finnish Forest Research Institute for young Scots pine stands (3 measurement times, 192 stands). The effect of precommercial thinning on stand development was examined on the basis of long-term experiments (13 stands, 169 plots). The effect of timing and intensity of the first commercial thinning on yield and profitability were based on measurements made in first commercial thinnings (27 stands of Metsähallitus), and the further stand development was modeled using the MOTTI simulator. The thesis was based on four articles and a summary. Stand level growth models were developed for young, even-aged Scots pine stands. The models reliably predicted the development up until the first commercial thinning stage. The stand density of young Scots pine stands in Finland was moderately low compared to the target values. In addition, the external quality of pines was low on average. The low stand density and poor external quality will result in the need for quality tree selection in thinnings, if high quality sawn timber is required. In Northern Finland, only 20% of the dominant trees were classified as normal. This will lead to the situation where external quality will remain relatively poor up until the end of rotation. Early and light precommercial thinning (Hdom 3 m, to a density of 3000 trees per hectare) increased the thinning removal by 40% compared to late and more intensive precommercial thinning (at 7 meters to a density of 2000 trees per hectare). A model for the effect of precommercial thinning on merchantable thinning removal at the first commercial thinning was developed for forest management planning purposes. When the recommended time of first commercial thinning was delayed from a dominant height of 12 m to 16 m, or by ten years, the yield of merchantable wood was doubled. Simultaneously, the current value of the stumpage revenues (with 4% interest rate) was increased on the average by 65% (330 € per hectare). Variation in stumpage prices or interest rates did not have any effect on the final results. Without exception, delaying the first commercial thinning by ten years seemed to be the most profitable method. This presupposes that precommercial thinning has been carried out at the right time and that tree quality aspects do not be specially considered. Furthermore, the wood yield and economic outcome from the entire rotation were similar regardless of whether the first thinning was performed at the time currently recommended or ten years later.

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

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