Geospatial environmental data modelling applications using remote sensing, GIS and spatial statistics

This thesis presents novel modelling applications for environmental geospatial data using remote sensing, GIS and statistical modelling techniques. The studied themes can be classified into four main themes: (i) to develop advanced geospatial databases. Paper (I) demonstrates the creation of a geospatial database for the Glanville fritillary butterfly (Melitaea cinxia) in the Åland Islands, south-western Finland; (ii) to analyse species diversity and distribution using GIS techniques. Paper (II) presents a diversity and geographical distribution analysis for Scopulini moths at a world-wide scale; (iii) to study spatiotemporal forest cover change. Paper (III) presents a study of exotic and indigenous tree cover change detection in Taita Hills Kenya using airborne imagery and GIS analysis techniques; (iv) to explore predictive modelling techniques using geospatial data. In Paper (IV) human population occurrence and abundance in the Taita Hills highlands was predicted using the generalized additive modelling (GAM) technique. Paper (V) presents techniques to enhance fire prediction and burned area estimation at a regional scale in East Caprivi Namibia. Paper (VI) compares eight state-of-the-art predictive modelling methods to improve fire prediction, burned area estimation and fire risk mapping in East Caprivi Namibia. The results in Paper (I) showed that geospatial data can be managed effectively using advanced relational database management systems. Metapopulation data for Melitaea cinxia butterfly was successfully combined with GPS-delimited habitat patch information and climatic data. Using the geospatial database, spatial analyses were successfully conducted at habitat patch level or at more coarse analysis scales. Moreover, this study showed it appears evident that at a large-scale spatially correlated weather conditions are one of the primary causes of spatially correlated changes in Melitaea cinxia population sizes. In Paper (II) spatiotemporal characteristics of Socupulini moths description, diversity and distribution were analysed at a world-wide scale and for the first time GIS techniques were used for Scopulini moth geographical distribution analysis. This study revealed that Scopulini moths have a cosmopolitan distribution. The majority of the species have been described from the low latitudes, sub-Saharan Africa being the hot spot of species diversity. However, the taxonomical effort has been uneven among biogeographical regions. Paper III showed that forest cover change can be analysed in great detail using modern airborne imagery techniques and historical aerial photographs. However, when spatiotemporal forest cover change is studied care has to be taken in co-registration and image interpretation when historical black and white aerial photography is used. In Paper (IV) human population distribution and abundance could be modelled with fairly good results using geospatial predictors and non-Gaussian predictive modelling techniques. Moreover, land cover layer is not necessary needed as a predictor because first and second-order image texture measurements derived from satellite imagery had more power to explain the variation in dwelling unit occurrence and abundance. Paper V showed that generalized linear model (GLM) is a suitable technique for fire occurrence prediction and for burned area estimation. GLM based burned area estimations were found to be more superior than the existing MODIS burned area product (MCD45A1). However, spatial autocorrelation of fires has to be taken into account when using the GLM technique for fire occurrence prediction. Paper VI showed that novel statistical predictive modelling techniques can be used to improve fire prediction, burned area estimation and fire risk mapping at a regional scale. However, some noticeable variation between different predictive modelling techniques for fire occurrence prediction and burned area estimation existed.

Characterisation and processing of amorphous binary mixtures with low glass transition temperature

The number of drug substances in formulation development in the pharmaceutical industry is increasing. Some of these are amorphous drugs and have glass transition below ambient temperature, and thus they are usually difficult to formulate and handle. One reason for this is the reduced viscosity, related to the stickiness of the drug, that makes them complicated to handle in unit operations. Thus, the aim in this thesis was to develop a new processing method for a sticky amorphous model material. Furthermore, model materials were characterised before and after formulation, using several characterisation methods, to understand more precisely the prerequisites for physical stability of amorphous state against crystallisation. The model materials used were monoclinic paracetamol and citric acid anhydrate. Amorphous materials were prepared by melt quenching or by ethanol evaporation methods. The melt blends were found to have slightly higher viscosity than the ethanol evaporated materials. However, melt produced materials crystallised more easily upon consecutive shearing than ethanol evaporated materials. The only material that did not crystallise during shearing was a 50/50 (w/w, %) blend regardless of the preparation method and it was physically stable at least two years in dry conditions. Shearing at varying temperatures was established to measure the physical stability of amorphous materials in processing and storage conditions. The actual physical stability of the blends was better than the pure amorphous materials at ambient temperature. Molecular mobility was not related to the physical stability of the amorphous blends, observed as crystallisation. Molecular mobility of the 50/50 blend derived from a spectral linewidth as a function of temperature using solid state NMR correlated better with the molecular mobility derived from a rheometer than that of differential scanning calorimetry data. Based on the results obtained, the effect of molecular interactions, thermodynamic driving force and miscibility of the blends are discussed as the key factors to stabilise the blends. The stickiness was found to be affected glass transition and viscosity. Ultrasound extrusion and cutting were successfully tested to increase the processability of sticky material. Furthermore, it was found to be possible to process the physically stable 50/50 blend in a supercooled liquid state instead of a glassy state. The method was not found to accelerate the crystallisation. This may open up new possibilities to process amorphous materials that are otherwise impossible to manufacture into solid dosage forms.

Language in Acquisition : Early Lexical Development and Associations between Lexicon and Grammar Findings from Full-Term and Very-Low-Birth-Weight Finnish Children

The aim was to analyse the growth and compositional development of the receptive and expressive lexicons between the ages 0,9 and 2;0 in the full-term (FT) and the very-low-birth-weight (VLBW) children who are acquiring Finnish. The associations between the expressive lexicon and grammar at 1;6 and 2;0 in the FT children were also studied. In addition, the language skills of the VLBW children at 2;0 were analysed, as well as the predictive value of early lexicon to the later language performance. Four groups took part in the studies: the longitudinal (N = 35) and cross-sectional (N = 146) samples of the FT children, and the longitudinal (N = 32) and cross-sectional (N = 66) samples of VLBW children. The data was gathered by applying of the structured parental rating method (the Finnish version of the Communicative Development Inventory), through analysis of the children´s spontaneous speech and by administering a a formal test (Reynell Developmental Language Scales). The FT children acquired their receptive lexicons earlier, at a faster rate and with larger individual variation than their expressive lexicons. The acquisition rate of the expressive lexicon increased from slow to faster in most children (91%). Highly parallel developmental paths for lexical semantic categories were detected in the receptive and expressive lexicons of the Finnish children when they were analysed in relation to the growth of the lexicon size, as described in the literature for children acquiring other languages. The emergence of grammar was closely associated with expressive lexical growth. The VLBW children acquired their receptive lexicons at a slower rate and had weaker language skills at 2;0 than the full-term children. The compositional development of both lexicons happened at a slower rate in the VLBW children when compared to the FT controls. However, when the compositional development was analysed in relation to the growth of lexicon size, this development occurred qualitatively in a nearly parallel manner in the VLBW children as in the FT children. Early receptive and expressive lexicon sizes were significantly associated with later language skills in both groups. The effect of the background variables (gender, length of the mother s basic education, birth weight) on the language development in the FT and the VLBW children differed. The results provide new information of early language acquisition by the Finnish FT and VLBW children. The results support the view that the early acquisition of the semantic lexical categories is related to lexicon growth. The current findings also propose that the early grammatical acquisition is closely related to the growth of expressive vocabulary size. The language development of the VLBW children should be followed in clinical work.

Epidemiological data of seasonal variation in mood and behaviour

The aim of this study was to measure seasonal variation in mood and behaviour. The dual vulnerability and latitude effect hypothesis, the risk of increased appetite, weight and other seasonal symptoms to develop metabolic syndrome, and perception of low illumination in quality of life and mental well-being were assessed. These variations are prevalent in persons who live in high latitudes and need balancing of metabolic processes to adapt to environmental changes due to seasons. A randomized sample of 8028 adults aged 30 and over (55% women) participated in an epidemiological health examination study, The Health 2000, applying the probability proportional to population size method for a range of socio-demographic characteristics. They were present in a face-to-face interview at home and health status examination. The questionnaires included the modified versions of the Seasonal Pattern Assessment Questionnaire (SPAQ) and Beck Depression Inventory (BDI), the Health Related Quality of Life (HRQoL) instrument 15D, and the General Health Questionnaire (GHQ). The structured and computerized Munich Composite International Diagnostic Interview (M-CIDI) as part of the interview was used to assess diagnoses of mental disorders, and, the National Cholesterol Education Program Adult Treatment Panel III (NCEP-ATPIII) criteria were assessed using all the available information to detect metabolic syndrome. A key finding was that 85% of this nationwide representative sample had seasonal variation in mood and behaviour. Approximately 9% of the study population presented combined seasonal and depressive symptoms with a significant association between their scores, and 2.6% had symptoms that corresponded to Seasonal Affective Disorder (SAD) in severity. Seasonal variations in weight and appetite are two important components that increase the risk of metabolic syndrome. Other factors such as waist circumference and major depressive disorder contributed to the metabolic syndrome as well. Persons reported of having seasonal symptoms were associated with a poorer quality of life and compromised mental well-being, especially if indoors illumination at home and/or at work was experienced as being low. Seasonal and circadian misalignments are suggested to associate with metabolic disorders, and could be remarked if individuals perceive low illumination levels at home and/or at work that affect the health-related quality of life and mental well-being. Keywords: depression, health-related quality of life, illumination, latitude, mental well-being, metabolic syndrome, seasonal variation, winter.

The effect of temperature on height growth of Scots pine in northern Finland

The effect of temperature on height growth of Scots pine in the northern boreal zone in Lapland was studied in two different time scales. Intra-annual growth was monitored in four stands in up to four growing seasons using an approximately biweekly measurement interval. Inter-annual growth was studied using growth records representing seven stands and five geographical locations. All the stands were growing on a dry to semi-dry heath that is a typical site type for pine stands in Finland. The applied methodology is based on applied time-series analysis and multilevel modelling. Intra-annual elongation of the leader shoot correlated with temperature sum accumulation. Height growth ceased when, on average, 41% of the relative temperature sum of the site was achieved (observed minimum and maximum were 38% and 43%). The relative temperature sum was calculated by dividing the actual temperature sum by the long-term mean of the total annual temperature sum for the site. Our results suggest that annual height growth ceases when a location-specific temperature sum threshold is attained. The positive effect of the mean July temperature of the previous year on annual height increment proved to be very strong at high latitudes. The mean November temperature of the year before the previous had a statistically significantly effect on height increment in the three northernmost stands. The effect of mean monthly precipitation on annual height growth was statistically insignificant. There was a non-linear dependence between length and needle density of annual shoots. Exceptionally low height growth results in high needle-density, but the effect is weaker in years of average or good height growth. Radial growth and next year s height growth are both largely controlled by current July temperature. Nevertheless, their growth variation in terms of minimum and maximum is not necessarily strongly correlated. This is partly because height growth is more sensitive to changes in temperature. In addition, the actual effective temperature period is not exactly the same for these two growth components. Yet, there is a long-term balance that was also statistically distinguishable; radial growth correlated significantly with height growth with a lag of 2 years. Temperature periods shorter than a month are more effective variables than mean monthly values, but the improvement is on the scale of modest to good when applying Julian days or growing-degree-days as pointers.

Postglacial climate changes and vegetation responses in northern Europe

Postglacial climate changes and vegetation responses were studied using a combination of biological and physical indicators preserved in lake sediments. Low-frequency trends, high-frequency events and rapid shifts in temperature and moisture balance were probed using pollen-based quantitative temperature reconstructions and oxygen-isotopes from authigenic carbonate and aquatic cellulose, respectively. Pollen and plant macrofossils were employed to shed light on the presence and response rates of plant populations in response to climate changes, particularly focusing on common boreal and temperate tree species. Additional geochemical and isotopic tracers facilitated the interpretation of pollen- and oxygen-isotope data. The results show that the common boreal trees were present in the Baltic region (~55°N) during the Lateglacial, which contrasts with the traditional view of species refuge locations in the south-European peninsulas during the glacial/interglacial cycles. The findings of this work are in agreement with recent paleoecological and genetic evidence suggesting that scattered populations of tree species persisted at higher latitudes, and that these taxa were likely limited to boreal trees. Moreover, the results demonstrate that stepwise changes in plant communities took place in concert with major climate fluctuations of the glacial/interglacial transition. Postglacial climate trends in northern Europe were characterized by rise, maxima and fall in temperatures and related changes in moisture balance. Following the deglaciation of the Northern Hemisphere and the early Holocene reorganization of the ice-ocean-atmosphere system, the long-term temperature trends followed gradually decreasing summer insolation. The early Holocene (~11,700-8000 cal yr BP) was overall cool, moist and oceanic, although the earliest Holocene effective humidity may have been low particularly in the eastern part of northern Europe. The gradual warming trend was interrupted by a cold event ~8200 cal yr BP. The maximum temperatures, ~1.5-3.0°C above modern values, were attained ~8000-4000 cal yr BP. This mid-Holocene peak warmth was coupled with low lake levels, low effective humidity and summertime drought. The late Holocene (~4000 cal yr BP-present) was characterized by gradually decreasing temperatures, higher lake levels and higher effective humidity. Moreover, the gradual trends of the late Holocene were probably superimposed by higher-frequency variability. The spatial variability of the Holocene temperature and moisture balance patterns were tentatively attributed to the differing heat capacities of continents and oceans, changes in atmospheric circulation modes and position of sites and subregions with respect to large water bodies and topographic barriers. The combination of physical and biological proxy archives is a pivotal aspect of this work, because non-climatic factors, such as postglacial migration, disturbances and competitive interactions, can influence reshuffling of vegetation and hence, pollen-based climate reconstructions. The oxygen-isotope records and other physical proxies presented in this work manifest that postglacial climate changes were the main driver of the establishment and expansion of temperate and boreal tree populations, and hence, large-scale and long-term vegetation patterns were in dynamic equilibrium with climate. A notable exception to this pattern may be the postglacial invasion of Norway spruce and the related suppression of mid-Holocene temperate forest. This salient step in north-European vegetation history, the development of the modern boreal ecosystem, cannot be unambiguously explained by current evidence of postglacial climate changes. The results of this work highlight that plant populations, including long-lived trees, may be able to respond strikingly rapidly to changes in climate. Moreover, interannual and seasonal variation and extreme events can exert an important influence on vegetation reshuffling. Importantly, the studies imply that the presence of diffuse refuge populations or local stands among the prevailing vegetation may have provided the means for extraordinarily rapid vegetation responses. Hence, if scattered populations are not provided and tree populations are to migrate long distances, their capacity to keep up with predicted rates of future climate change may be lower than previously thought.

Remote sensing of boreal land cover: estimation of forest attributes and extent

Remote sensing provides methods to infer land cover information over large geographical areas at a variety of spatial and temporal resolutions. Land cover is input data for a range of environmental models and information on land cover dynamics is required for monitoring the implications of global change. Such data are also essential in support of environmental management and policymaking. Boreal forests are a key component of the global climate and a major sink of carbon. The northern latitudes are expected to experience a disproportionate and rapid warming, which can have a major impact on vegetation at forest limits. This thesis examines the use of optical remote sensing for estimating aboveground biomass, leaf area index (LAI), tree cover and tree height in the boreal forests and tundra taiga transition zone in Finland. The continuous fields of forest attributes are required, for example, to improve the mapping of forest extent. The thesis focus on studying the feasibility of satellite data at multiple spatial resolutions, assessing the potential of multispectral, -angular and -temporal information, and provides regional evaluation for global land cover data. Preprocessed ASTER, MISR and MODIS products are the principal satellite data. The reference data consist of field measurements, forest inventory data and fine resolution land cover maps. Fine resolution studies demonstrate how statistical relationships between biomass and satellite data are relatively strong in single species and low biomass mountain birch forests in comparison to higher biomass coniferous stands. The combination of forest stand data and fine resolution ASTER images provides a method for biomass estimation using medium resolution MODIS data. The multiangular data improve the accuracy of land cover mapping in the sparsely forested tundra taiga transition zone, particularly in mires. Similarly, multitemporal data improve the accuracy of coarse resolution tree cover estimates in comparison to single date data. Furthermore, the peak of the growing season is not necessarily the optimal time for land cover mapping in the northern boreal regions. The evaluated coarse resolution land cover data sets have considerable shortcomings in northernmost Finland and should be used with caution in similar regions. The quantitative reference data and upscaling methods for integrating multiresolution data are required for calibration of statistical models and evaluation of land cover data sets. The preprocessed image products have potential for wider use as they can considerably reduce the time and effort used for data processing.

Environmental change and anthropogenic impact on lake sediments during the Holocene in the Finnish - Karelian inland area

This thesis discusses the prehistoric human disturbance during the Holocene by means of case studies using detailed high-resolution pollen analysis from lake sediment. The four lakes studied are situated between 61o 40' and 61o 50' latitudes in the Finnish Karelian inland area and vary between 2.4 and 28.8 ha in size. The existence of Early Metal Age population was one important question. Another study question concerned the development of grazing, and the relationship between slash-and-burn cultivation and permanent field cultivation. The results were presented as pollen percentages and pollen concentrations (grains cm 3). Accumulation values (grains cm 2 yr 1) were calculated for Lake Nautajärvi and Lake Orijärvi sediment, where the sediment accumulation rate was precisely determined. Sediment properties were determined using loss-on-ignition (LOI) and magnetic susceptibility (k). Dating methods used include both conventional and AMS 14C determinations, paleomagnetic dating and varve choronology. The isolation of Lake Kirjavalampi on the northern shore of Lake Ladoga took place ca. 1460 1300 BC. The long sediment cores from Finland, Lake Kirkkolampi and Lake Orijärvi in southeastern Finland and Lake Nautajärvi in south central Finland all extended back to the Early Holocene and were isolated from the Baltic basin ca. 9600 BC, 8600 BC and 7675 BC, respectively. In the long sediment cores, the expansion of Alnus was visible between 7200 - 6840 BC. The spread of Tilia was dated in Lake Kirkkolampi to 6600 BC, in Lake Orijärvi to 5000 BC and at Lake Nautajärvi to 4600 BC. Picea is present locally in Lake Kirkkolampi from 4340 BC, in Lake Orijärvi from 6520 BC and in Lake Nautajärvi from 3500 BC onwards. The first modifications in the pollen data, apparently connected to anthropogenic impacts, were dated to the beginning of the Early Metal Period, 1880 1600 BC. Anthropogenic activity became clear in all the study sites by the end of the Early Metal Period, between 500 BC AD 300. According to Secale pollen, slash-and-burn cultivation was practised around the eastern study lakes from AD 300 600 onwards, and at the study site in central Finland from AD 880 onwards. The overall human impact, however, remained low in the studied sites until the Late Iron Age. Increasing human activity, including an increase in fire frequency was detected from AD 800 900 onwards in the study sites in eastern Finland. In Lake Kirkkolampi, this included cultivation on permanent fields, but in Lake Orijärvi, permanent field cultivation became visible as late as AD 1220, even when the macrofossil data demonstrated the onset of cultivation on permanent fields as early as the 7th century AD. On the northern shore of Lake Ladoga, local activity became visible from ca. AD 1260 onwards and at Lake Nautajärvi, sediment the local occupation was traceable from 1420 AD onwards. The highest values of Secale pollen were recorded both in Lake Orijärvi and Lake Kirjavalampi between ca. AD 1700 1900, and could be associated with the most intensive period of slash-and-burn from AD 1750 to 1850 in eastern Finland.

Paleomagnetic and rock magnetic study with emphasis on the Precambrian intrusions and impact structures in Fennoscandia and South Africa

The importance of supercontinents in our understanding of the geological evolution of the planet Earth has been recently emphasized. The role of paleomagnetism in reconstructing lithospheric blocks in their ancient paleopositions is vital. Paleomagnetism is the only quantitative tool for providing ancient latitudes and azimuthal orientations of continents. It also yields information of content of the geomagnetic field in the past. In order to obtain a continuous record on the positions of continents, dated intrusive rocks are required in temporal progression. This is not always possible due to pulse-like occurrences of dykes. In this work we demonstrate that studies of meteorite impact-related rocks may fill some gaps in the paleomagnetic record. This dissertation is based on paleomagnetic and rock magnetic data obtained from samples of the Jänisjärvi impact structure (Russian Karelia, most recent 40Ar-39Ar age of 682 Ma), the Salla diabase dyke (North Finland, U-Pb 1122 Ma), the Valaam monzodioritic sill (Russian Karelia, U-Pb 1458 Ma), and the Vredefort impact structure (South Africa, 2023 Ma). The paleomagnetic study of Jänisjärvi samples was made in order to obtain a pole for Baltica, which lacks paleomagnetic data from 750 to ca. 600 Ma. The position of Baltica at ca. 700 Ma is relevant in order to verify whether the supercontinent Rodinia was already fragmented. The paleomagnetic study of the Salla dyke was conducted to examine the position of Baltica at the onset of supercontinent Rodinia's formation. The virtual geomagnetic pole (VGP) from Salla dyke provides hints that the Mesoproterozoic Baltica - Laurentia unity in the Hudsonland (Columbia, Nuna) supercontinent assembly may have lasted until 1.12 Ga. Moreover, the new VGP of Salla dyke provides new constraint on the timing of the rotation of Baltica relative to Laurentia (e.g. Gower et al., 1990). A paleomagnetic study of the Valaam sill was carried out in order to shed light into the question of existence of Baltica-Laurentia unity in the supercontinent Hudsonland. Combined with results from dyke complex of the Lake Ladoga region (Schehrbakova et al., 2008) a new robust paleomagnetic pole for Baltica is obtained. This pole places Baltica on a latitude of 10°. This low latitude location is supported also by Mesoproterozoic 1.5 1.3 Ga red-bed sedimentation (for example the Satakunta sandstone). The Vredefort impactite samples provide a well dated (2.02 Ga) pole for the Kaapvaal Craton. Rock magnetic data reveal unusually high Koenigsberger ratios (Q values) in all studied lithologies of the Vredefort dome. The high Q values are now first time also seen in samples from the Johannesburg Dome (ca. 120 km away) where there is no impact evidence. Thus, a direct causative link of high Q values to the Vredefort impact event can be ruled out.

The high- and low-activity forms of human plasma phospholipid transfer protein (PLTP)

Plasma phospholipid transfer protein (PLTP) plays a crucial role in high-density lipoprotein (HDL) metabolism and reverse cholesterol transport (RCT). It mediates the generation of pre-beta-HDL particles, enhances the cholesterol efflux from peripheral cells to pre-beta-HDL, and metabolically maintains the plasma HDL levels by facilitating the transfer of post-lipolytic surface remnants of triglyceride-rich lipoproteins to HDL. In addition to the antiatherogenic properties, recent findings indicate that PLTP has also proatherogenic characteristics, and that these opposite characteristics of PLTP are dependent on the site of PLTP expression and action. In human plasma, PLTP exists in a high-activity (HA-PLTP) and a low-activity form (LA-PLTP), which are associated with macromolecular complexes of different size and composition. The aims of this thesis were to isolate the two PLTP forms from human plasma, to characterize the molecular complexes in which the HA- and LA-PLTP reside, and to study the interactions of the PLTP forms with apolipoproteins (apo) and the ability of apolipoproteins to regulate PLTP activity. In addition, we aimed to study the distribution of the two PLTP forms in a Finnish population sample as well as to find possible regulatory factors for PLTP by investigating the influence of lipid and glucose metabolism on the balance between the HA- and LA-PLTP. For these purposes, an enzyme-linked immunosorbent assay (ELISA) capable of determining the serum total PLTP concentration and quantitating the two PLTP forms separately was developed. In this thesis, it was demonstrated that the HA-PLTP isolated from human plasma copurified with apoE, whereas the LA-PLTP formed a complex with apoA-I. The separation of these two PLTP forms was carried out by a dextran sulfate (DxSO4)-CaCl2 precipitation of plasma samples before the mass determination. A similar immunoreactivity of the two PLTP forms in the ELISA could be reached after a partial sample denaturation by SDS. Among normolipidemic Finnish individuals, the mean PLTP mass was 6.6 +/- 1.5 mg/l and the mean PLTP activity 6.6 +/- 1.7 umol/ml/h. Of the serum PLTP concentration, almost 50% represented HA-PLTP. The results indicate that plasma HDL levels could regulate PLTP concentration, while PLTP activity could be regulated by plasma triglyceride-rich very low-density lipoprotein (VLDL) concentration. Furthermore, new evidence is presented that PLTP could also play a role in glucose metabolism. Finally, both PLTP forms were found to interact with apoA-I, apoA-IV, and apoE. In addition, both apoE and apoA-IV, but not apoA-I, were capable of activating the LA-PLTP. These findings suggest that the distribution of the HA- and LA-PLTP in human plasma is subject to dynamic regulation by apolipoproteins.

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

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