Search for genetic variants influencing human height

Human growth and attained height are determined by a combination of genetic and environmental effects and in modern Western societies > 80% of the observed variation in height is determined by genetic factors. Height is a fundamental human trait that is associated with many socioeconomic and psychosocial factors and health measures, however little is known of the identity of the specific genes that influence height variation in the general population. This thesis work aimed to identify the genetic variants that influence height in the general population by genome-wide linkage analysis utilizing large family samples. The study focused on analysis of three separate sets of families consisting of: 1) 1,417 individuals from 277 Finnish families (FinnHeight), 2) 8,450 individuals from 3,817 families from Australia and Europe (EUHeight) and 3) 9,306 individuals from 3,302 families from the United States (USHeight). The most significant finding in this study was found in the Finnish family sample where we a locus in the chromosomal region 1p21 was linked to adult height. Several regions showed evidence for linkage in the Australian, European and US families with 8q21 and 15q25 being the most significant. The region on 1p21 was followed up with further studies and we were able to show that the collagen 11-alpha-1 gene (COL11A1) residing at this location was associated with adult height. This association was also confirmed in an independent Finnish population cohort (Health 2000) consisting of 6,542 individuals. From this population sample, we estimated that homozygous males and females for this gene variant were 1.1 and 0.6 cm taller than the respective controls. In this thesis work we identified a gene variant in the COL11A1 gene that influences human height, although this variant alone explains only 0.1% of height variation in the Finnish population. We also demonstrated in this study that special stratification strategies such as performing sex-limited analyses, focusing on dizygous twin pairs, analyzing ethnic groups within a population separately and utilizing homogenous populations such as the Finns can improve the statistical power of finding QTL significantly. Also, we concluded from the results of this study that even though genetic effects explain a great proportion of height variance, it is likely that there are tens or even hundreds of genes with small individual effects underlying the genetic architecture of height.

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).

Novel Multiple Sclerosis Predisposing Genetic Variants Outside the HLA Region

Multiple sclerosis (MS) is a chronic inflammatory disease of the central nervous system (CNS). Both environmental factors and several predisposing genes are required to generate MS. Despite intensive research these risk factors are still largely unknown, the pathogenesis of MS demyelination is poorly understood, and no curative treatment exists. Both prevalence and familial occurrence of MS are exceptionally high in a Finnish population subisolate, Southern Ostrobothnia, presumably due to enrichment of predisposing genetic variants within this region. Previous linkage scan on MS pedigrees from Southern Ostrobothnia detected three main MS loci on chromosomes 5p, 6p (HLA) and 17q. Linkage studies in other populations have also provided independent evidence for the location of MS susceptibility genes in these regions. Further, these loci are syntenic to the experimental autoimmune encephalomyelitis (EAE) susceptibility loci of rodents. In this thesis work an effort was made to localize MS predisposing alleles of the linked loci outside the HLA region by studying familial MS cases from the Southern Ostrobothnia isolate. Analysis of the 5p locus revealed one region, flanking the complement component 7 (C7) gene. The identified relatively rare haplotype seems to have a fairly large effect on genetic susceptibility of MS (frequency MS 12%, controls 4%; p=0.000003, OR=2.73). Evidence for association with alleles of the region and MS was seen also in more heterogeneous populations. Convincingly, plasma C7 protein levels and complement activity correlated with the risk haplotype identified. The finding stimulated us to study other complement cascade genes in MS. No evidence for association could be observed with the complement component coding genes outside 5p. A scan of the 17q locus provided evidence for association with variants of the protein kinase C alpha (PRKCA) gene (p=0.0001). Modest evidence for association with PRKCA was observed also in Canadian MS families. Finally we used a candidate gene based approach to identify potential MS loci. Mutations of DAP12 and TREM2 cause a recessively inherited CNS white matter disease PLOSL. Interestingly, DAP12 and TREM2 are located in MS regions on 6p and 19q, and we tested them as potential candidate genes in the Finnish MS sample. No evidence for association with MS was observed. This thesis provides an example of how extended families from special populations can be utilized in fine-mapping of the linked loci. A first relatively rare MS variant was identified utilizing the strength of a Finnish population subisolate. This variant seems to have an effect on activity of the complement system, which has previously been suggested to have an important role in the pathogenesis of MS.

Studies of the genetic epidemiology of cardiovascular disease: focus on inflammatory candidate genes

Cardiovascular disease (CVD) is a complex disease with multifactorial aetiology. Both genetic and environmental factors contribute to the disease risk. The lifetime risk for CVD differs markedly between men and women, men being at increased risk. Inflammatory reaction contributes to the development of the disease by promoting atherosclerosis in artery walls. In the first part of this thesis, we identified several inflammatory related CVD risk factors associating with the amount of DNA from whole blood samples, indicating a potential source of bias if a genetic study selects the participants based on the available amount of DNA. In the following studies, this observation was taken into account by applying whole genome amplification to samples otherwise subjected to exclusion due to very low DNA yield. We continued by investigating the contribution of inflammatory genes to the risk for CVD separately in men and women, and looked for sex-genotype interaction. In the second part, we explored a new candidate gene and its role in the risk for CVD. Selenoprotein S (SEPS1) is a membrane protein residing in the endoplasmic reticulum where it participates in retro-translocation of unfolded proteins to cytosolic protein degradation. Previous studies have indicated that SEPS1 protects cells from oxidative stress and that variations in the gene are associated with circulating levels of inflammatory cytokines. In our study, we identified two variants in the SEPS1 gene, which associated with coronary heart disease and ischemic stroke in women. This is, to our knowledge, the first study suggesting a role of SEPS1 in the risk for CVD after extensively examining the variation within the gene region. In the third part of this thesis, we focused on a set of seven genes (angiotensin converting enzyme, angiotensin II receptor type I, C-reactive protein (CRP), and fibrinogen alpha-, beta-, and gamma-chains (FGA, FGB, FGG)) related to inflammatory cytokine interleukin 6 (IL6) and their association with the risk for CVD. We identified one variant in the IL6 gene conferring risk for CVD in men and a variant pair from IL6 and FGA genes associated with decreased risk. Moreover, we identified and confirmed an association between a rare variant in the CRP gene and lower CRP levels, and found two variants in the FGA and FGG genes associating with fibrinogen. The results from this third study suggest a role for the interleukin 6 pathway genes in the pathogenesis of CVD and warrant further studies in other populations. In addition to the IL6 -related genes, we describe in this thesis several sex-specific associations in other genes included in this study. The majority of the findings were evident only in women encouraging other studies of cardiovascular disease to include and analyse women separately from men.

Novel prognostic factors in chronic myeloid leukemia

Chronic myeloid leukemia (CML) is a malignant clonal blood disease that originates from a pluripotent hematopoietic stem cell. The cytogenetic hallmark of CML, the Philadelphia chromosome (Ph), is formed as a result of reciprocal translocation between chromosomes 9 and 22, which leads to a formation of a chimeric BCR-ABL fusion gene. The BCR-ABL protein is a constitutively active tyrosine kinase that changes the adhesion properties of cells, constitutively activates mitogenic signaling, enhances cell proliferation and reduces apoptosis. This results in leukemic growth and the clinical disease, CML. With the advent of targeted therapies against the BCR-ABL fusion protein, the treatment of CML has changed considerably during the recent decade. In this thesis, the clinical significance of different diagnostic methods and new prognostic factors in CML have been assessed. First, the association between two different methods for measuring CML disease burden (the RQ-PCR and the high mitotic index metaphase FISH) was assessed in bone marrow and peripheral blood samples. The correlation between positive RQ-PCR and metaphase FISH samples was high. However, RQ-PCR was more sensitive and yielded measurable transcripts in 40% of the samples that were negative by metaphase FISH. The study established a laboratory-specific conversion factor for setting up the International Scale when standardizing RQ-PCR measurements. Secondly, the amount of minimal residual disease (MRD) after allogeneic hematopoietic stem cell transplantation (alloHSCT) was determined. For this, metaphase FISH was done for the bone marrow samples of 102 CML patients. Most (68%), had no residual cells during the entire follow-up time. Some (12 %) patients had minor (<1%) MRD which decreased even further with time, whereas 19% had a progressive rise in MRD that exceeded 1% or had more than 1% residual cells when first detected. Residual cells did not become eradicated spontaneously if the frequency of Ph+ cells exceeded 1% during follow-up. Next, the impact of deletions in the derivative chromosome 9, was examined. Deletions were observed in 15% of the CML patients who later received alloHSCT. After alloHSCT, there was no difference in the total relapse rate in patients with or without deletions. Nor did the estimates of overall survival, transplant-related mortality, leukemia-free survival and relapse-free time show any difference between these groups. When conventional treatment regimens are used, the der(9) status could be an important criterion, in conjunction with other prognostic factors, when allogeneic transplantation is considered. The significance of der(9) deletions for patients treated with tyrosine kinase inhibitors is not clear and requires further investigation. In addition to the der(9) status of the patient, the significance of bone marrow lymphocytosis as a prognostic factor in CML was assessed. Bone marrow lymphocytosis during imatinib therapy was a positive predictive factor and heralded optimal response. When combined with major cytogenetic response at three months of treatment, bone marrow lymphocytosis predicted a prognostically important major molecular response at 18 months of imatinib treatment. Although the validation of these findings is warranted, the determination of the bone marrow lymphocyte count could be included in the evaluation of early response to imatinib treatment already now. Finally, BCR-ABL kinase domain mutations were studied in CML patients resistant against imatinib treatment. Point mutations detected in the kinase domain were the same as previously reported, but other sequence variants, e.g. deletions or exon splicing, were also found. The clinical significance of the other variations remains to be determined.

Vascular growth factors and progenitor cells in cardiac allograft arteriosclerosis

Heart transplantation is the only therapeutic modality for many end-stage heart diseases but poor long-term survival remains a challenging problem. This is mainly due to the development of cardiac allograft arteriosclerosis (TxCAD) that is an accelerated form of coronary artery disease. Both traditional cardiovascular and transplantation-related risk factors for TxCAD have been identified but options for therapy are limited. TxCAD involves dysfunction of cardiac allograft vascular cells. Activated endothelial cells (EC) regulate allograft inflammation and secrete smooth muscle cell (SMC) growth factors. In turn, SMC and their progenitors invade the intima of the injured vessels and occlude the affected coronary arteries. Different vascular growth factors have to be delicately regulated in normal vascular development. In the present study, experimental heterotopic transplantation models were used to study the role of angiogenic and pro-inflammatory vascular endothelial growth factor (VEGF), EC growth factor angiopoietin (Ang), and SMC mitogen platelet-derived growth factor (PDGF) in the development of TxCAD. Pharmacological and gene transfer approaches were used to target these growth factors and to assess their therapeutic potential. This study shows that alloimmune response in heart transplants upregulates VEGF expression, and induces allograft angiogenesis that involves donor-derived primitive EC. Intracoronary adenoviral VEGF gene transfer increased macrophage infiltration, intimal angiogenesis and TxCAD. VEGF inhibition with PTK787 decreased allograft inflammation and TxCAD, and simultaneous PDGF inhibition with imatinib further decreased TxCAD. Specific inhibition of two VEGF-receptors (VEGFR) decreased allograft inflammation and TxCAD, and VEGFR-2 inhibition normalized the density of primitive and mature capillaries in the allografts. Adenovirus-mediated transient Ang1 expression in the allograft had anti-inflammatory and anti-arteriosclerotic effects. Adeno-associated virus (AAV)-mediated prolonged Ang1 or Ang2 expression had similar anti-inflammatory effects. However, AAV-Ang1 activated allograft SMC whereas AAV-Ang2 had no effects on SMC activation and decreased the development of TxCAD. These studies indicate an interplay of inflammation, angiogenesis and arteriosclerosis in cardiac allografts, and show that vascular growth factors are important regulators in the process. Also, VEGF inhibition, PDGF inhibition and angiopoietin therapy with clinically-relevant pharmacological agents or novel gene therapy approaches may counteract vascular dysfunction in cardiac allografts, and have beneficial effects on the survival of heart transplant patients in the future.

The effects of environmental factors on biomass and microcystin production by the freshwater cyanobacterial genera Microcystis and Anabaena

Several cyanobacterial genera produce the hepatotoxins, microcystins. Microcystins are produced only in cells that have microcystin synthetase gene (mcy) clusters, which encode enzyme complexes involved in microcystin biosynthesis. Microcystin-producing and nonmicrocystin-producing genotypes of single cyanobacterial genus may occur simultaneously in situ. Previously, the effects of environmental factors on the growth and microcystin production of cyanobacteria have mainly been studied by means of isolated cyanobacteria cultures in the laboratory. Studies in the field have been difficult, owing to the lack of methods to identify and quantify the different genotypes. In this study, genus-specific microcystin synthetase E (mcyE) gene primers were designed and a method to identify and quantify the mcyE copy numbers was developed and used in situ. Microcystis and Anabaena mcyE genes were observed in two Finnish lakes. Microcystis appeared to be the most abundant microcystin producer in Lake Tuusulanjärvi and in one basin of Lake Hiidenvesi. Because the most potent microcystin-producing genus of a lake can be identified, it will be possible in the future to design genus-targeted strategies for lake restoration. Effects of P and N concentrations on the biomass of microcystin-producing and nonmicrocystin-producing Microcystis strains and an Anabaena strain were studied in cultures. P and N concentrations and their combined effect increased cyanobacterial biomass of all Microcystis strains. The biomass of microcystin-producing Microcystis was higher than that of nonmicrocystin-producing strains at high nutrient concentrations. The P concentration increased Anabaena biomass, but the effect of N concentration was statistically insignificant for growth yield, probably due to the ability of the genus to fix molecular N2. P and N concentrations and combined nutrients caused an increase in cellular microcystin concentrations of the Microcystis strain cultivated in chemostat cultures. Cyanobacteria are able to hydrolyse nutrients from organic matter through extracellular enzyme activities. Leucine aminopeptidase (LAP) activity was observed in an axenic N2-fixing Anabaena strain grown in batch cultures. The P concentration caused a statistically significant increase in LAP activity, whereas the effect of N concentration was insignificant. The highest LAP activities were observed in the most eutrophic basins of Lake Hiidenvesi. LAP activity probably originated mostly from attached heterotrophic bacteria and less from cyanobacteria.

Factors affecting farmland prices in Finland

The purpose of the study was to analyse factors affecting the differences in land prices between regions. The key issue was to find out the policy effects on farmland prices. In addition to comprehensive literature review, a theoretical analysis as well as modern panel and spatial econometric techniques were utilized. The study clearly pointed out the importance of taking into account the possible spatial dependence. The data were exceptionally large, comprising more than 6 000 observations. Thus, it allowed a thorough econometric estimation including the possibility to take into account the spatial nature of the data. This study supports the view that there are many other factors that affect farmland prices besides pure agricultural returns. It was also found that the support clearly affects land prices. However, rather than assuming the discount rates for support and market returns to be similar, the rough analysis refers to the discount rate for support being a little lower. If this were true it would indicate that farmers rely more on support income than market returns. The results support the view presented in literature that land values are more responsive to government payments when these payments are perceived to be permanent. An important result of this study is that the structural differences between regions and the structural change in agriculture seemed to have a considerable role in affecting land prices. Firstly, the present structure affects the competition in the land market: the more dense farms are in the region the more there are potential buyers, and the land price increases. Secondly, the change in farm structure (especially in animal husbandry) connected to the policy changes that increase area-based support affects land prices. The effect comes from two sources. Growing farms need more land for the manure, and the proportion of retiring farmers may be lower. The introduction of the manure density variable proved to be an efficient way to aggregate the otherwise very difficult task of taking into account the environmental pressure caused by structural change in animal husbandry. Finally, infrastructure also has a very important role in determining the price level of agricultural land. If other industries are prospering in the surrounding area, agricultural viability also seems to improve. The non-farm opportunities offered to farm families make continuing and developing farming more tempting.

Environmental factors affecting the occurrence of periglacial landforms in Finnish Lapland: a numerical approach

Determination of the environmental factors controlling earth surface processes and landform patterns is one of the central themes in physical geography. However, the identification of the main drivers of the geomorphological phenomena is often challenging. Novel spatial analysis and modelling methods could provide new insights into the process-environment relationships. The objective of this research was to map and quantitatively analyse the occurrence of cryogenic phenomena in subarctic Finland. More precisely, utilising a grid-based approach the distribution and abundance of periglacial landforms were modelled to identify important landscape scale environmental factors. The study was performed using a comprehensive empirical data set of periglacial landforms from an area of 600 km2 at a 25-ha resolution. The utilised statistical methods were generalized linear modelling (GLM) and hierarchical partitioning (HP). GLMs were used to produce distribution and abundance models and HP to reveal independently the most likely causal variables. The GLM models were assessed utilising statistical evaluation measures, prediction maps, field observations and the results of HP analyses. A total of 40 different landform types and subtypes were identified. Topographical, soil property and vegetation variables were the primary correlates for the occurrence and cover of active periglacial landforms on the landscape scale. In the model evaluation, most of the GLMs were shown to be robust although the explanation power, prediction ability as well as the selected explanatory variables varied between the models. The great potential of the combination of a spatial grid system, terrain data and novel statistical techniques to map the occurrence of periglacial landforms was demonstrated in this study. GLM proved to be a useful modelling framework for testing the shapes of the response functions and significances of the environmental variables and the HP method helped to make better deductions of the important factors of earth surface processes. Hence, the numerical approach presented in this study can be a useful addition to the current range of techniques available to researchers to map and monitor different geographical phenomena.

Geological factors affecting on after-use of Finnish cut-over peatlands : with implications on the carbon accumulation

In Finland, peat harvesting sites are utilized down almost to the mineral soil. In this situation the properties of mineral subsoil are likely to have considerable influence on the suitability for the various after-use forms. The aims of this study were to recognize the chemical and physical properties of mineral subsoils possibly limiting the after-use of cut-over peatlands, to define a minimum practice for mineral subsoil studies and to describe the role of different geological areas. The future percentages of the different after-use forms were predicted, which made it possible to predict also carbon accumulation in this future situation. Mineral subsoils of 54 different peat production areas were studied. Their general features and grain size distribution was analysed. Other general items studied were pH, electrical conductivity, organic matter, water soluble nutrients (P, NO3-N, NH4-N, S and Fe) and exchangeable nutrients (Ca, Mg and K). In some cases also other elements were analysed. In an additional case study carbon accumulation effectiveness before the intervention was evaluated on three sites in Oulu area (representing sites typically considered for peat production). Areas with relatively sulphur rich mineral subsoil and pool-forming areas with very fine and compact mineral subsoil together covered approximately 1/5 of all areas. These areas were unsuitable for commercial use. They were recommended for example for mire regeneration. Another approximate 1/5 of the areas included very coarse or very fine sediments. Commercial use of these areas would demand special techniques - like using the remaining peat layer for compensating properties missing from the mineral subsoil. One after-use form was seldom suitable for one whole released peat production area. Three typical distribution patterns (models) of different mineral subsoils within individual peatlands were found. 57 % of studied cut-over peatlands were well suited for forestry. In a conservative calculation 26% of the areas were clearly suitable for agriculture, horticulture or energy crop production. If till without large boulders was included, the percentage of areas suitable to field crop production would be 42 %. 9-14 % of all areas were well suitable for mire regeneration or bird sanctuaries, but all areas were considered possible for mire regeneration with correct techniques. Also another 11 % was recommended for mire regeneration to avoid disturbing the mineral subsoil, so total 20-25 % of the areas would be used for rewetting. High sulphur concentrations and acidity were typical to the areas below the highest shoreline of the ancient Litorina sea and Lake Ladoga Bothnian Bay zone. Also differences related to nutrition were detected. In coarse sediments natural nutrient concentration was clearly higher in Lake Ladoga Bothnian Bay zone and in the areas of Svecokarelian schists and gneisses, than in Granitoid area of central Finland and in Archaean gneiss areas. Based on this study the recommended minimum analysis for after-use planning was for pH, sulphur content and fine material (<0.06 mm) percentage. Nutrition capacity could be analysed using the natural concentrations of calcium, magnesium and potassium. Carbon accumulation scenarios were developed based on the land-use predictions. These scenarios were calculated for areas in peat production and the areas released from peat production (59300 ha + 15 671 ha). Carbon accumulation of the scenarios varied between 0.074 and 0.152 million t C a-1. In the three peatlands considered for peat production the long term carbon accumulation rates varied between 13 and 24 g C m-2 a-1. The natural annual carbon accumulation had been decreasing towards the time of possible intervention.