Modeling ecological and evolutionary processes in spatially structured populations

Many species inhabit fragmented landscapes, resulting either from anthropogenic or from natural processes. The ecological and evolutionary dynamics of spatially structured populations are affected by a complex interplay between endogenous and exogenous factors. The metapopulation approach, simplifying the landscape to a discrete set of patches of breeding habitat surrounded by unsuitable matrix, has become a widely applied paradigm for the study of species inhabiting highly fragmented landscapes. In this thesis, I focus on the construction of biologically realistic models and their parameterization with empirical data, with the general objective of understanding how the interactions between individuals and their spatially structured environment affect ecological and evolutionary processes in fragmented landscapes. I study two hierarchically structured model systems, which are the Glanville fritillary butterfly in the Åland Islands, and a system of two interacting aphid species in the Tvärminne archipelago, both being located in South-Western Finland. The interesting and challenging feature of both study systems is that the population dynamics occur over multiple spatial scales that are linked by various processes. My main emphasis is in the development of mathematical and statistical methodologies. For the Glanville fritillary case study, I first build a Bayesian framework for the estimation of death rates and capture probabilities from mark-recapture data, with the novelty of accounting for variation among individuals in capture probabilities and survival. I then characterize the dispersal phase of the butterflies by deriving a mathematical approximation of a diffusion-based movement model applied to a network of patches. I use the movement model as a building block to construct an individual-based evolutionary model for the Glanville fritillary butterfly metapopulation. I parameterize the evolutionary model using a pattern-oriented approach, and use it to study how the landscape structure affects the evolution of dispersal. For the aphid case study, I develop a Bayesian model of hierarchical multi-scale metapopulation dynamics, where the observed extinction and colonization rates are decomposed into intrinsic rates operating specifically at each spatial scale. In summary, I show how analytical approaches, hierarchical Bayesian methods and individual-based simulations can be used individually or in combination to tackle complex problems from many different viewpoints. In particular, hierarchical Bayesian methods provide a useful tool for decomposing ecological complexity into more tractable components.

Coastal environmental gradients – Key to reproduction habitat mapping of freshwater fish in the Baltic Sea

Habitat requirements of fish are most strict during the early life stages, and the quality and quantity of reproduction habitats lays the basis for fish production. A considerable number of fish species in the northern Baltic Sea reproduce in the shallow coastal areas, which are also the most heavily exploited parts of the brackish marine area. However, the coastal fish reproduction habitats in the northern Baltic Sea are poorly known. The studies presented in this thesis focused on the influence of environmental conditions on the distribution of coastal reproduction habitats of freshwater fish. They were conducted in vegetated littoral zone along an exposure and salinity gradient extending from the innermost bays to the outer archipelago on the south-western and southern coasts of Finland, in the northern Baltic Sea. Special emphasis was placed on reed-covered Phragmites australis shores, which form a dominant vegetation type in several coastal archipelago areas. The main aims of this research were to (1) develop and test new survey and mapping methods, (2) investigate the environmental requirements that govern the reproduction of freshwater fish in the coastal area and (3) survey, map and model the distribution of the reproduction habitats of pike (Esox lucius) and roach (Rutilus rutilus). The white plate and scoop method with a standardized sampling time and effort was demonstrated to be a functional method for sampling the early life stages of fish in dense vegetation and shallow water. Reed-covered shores were shown to form especially important reproduction habitats for several freshwater fish species, such as pike, roach, other cyprinids and burbot, in the northern Baltic Sea. The reproduction habitats of pike were limited to sheltered reed- and moss-covered shores of the inner and middle archipelago, where suitable zooplankton prey were available and the influence of the open sea was low. The reproduction habitats of roach were even more limited and roach reproduction was successful only in the very sheltered reed-covered shores of the innermost bay areas, where salinity remained low (< 4‰) during the spawning season due to freshwater inflow. After identifying the critical factors restricting the reproduction of pike and roach, the spatial distribution of their reproduction habitats was successfully mapped and modelled along the environmental gradients using only a few environmental predictor variables. Reproduction habitat maps are a valuable tool promoting the sustainable use and management of exploited coastal areas and helping to maintain the sustainability of fish populations. However, the large environmental gradients and the extensiveness of the archipelago zone in the northern Baltic Sea demand an especially high spatial resolution of the coastal predictor variables. Therefore, the current lack of accurate large-scale, high-resolution spatial data gathered at exactly the right time is a considerable limitation for predictive modelling of shallow coastal waters.

Farmland birds and habitat heterogeneity in intensively cultivated boreal agricultural landscapes

The biodiversity of farmland ecosystems has decreased remarkably during the latter half of the 20th century, and this development is due to intensive farming with its various environmental effects. In the countries of the EU the Common Agricultural Policy (CAP) is the main determinant affecting farmland biodiversity, since the agricultural policy defines guidelines of agricultural practices. In addition to policies promoting intensive farming, CAP also includes national agri-environment schemes (AES), in which a part of subsidies paid to farmers is directed to acts that are presumed to promote environmental protection and biodiversity. In order to shape AES into relevant and powerful tools for biodiversity protection, detailed studies on the effects of agriculture on species and species assemblages are needed. In my thesis I investigated the importance of habitat heterogeneity and effects of different habitat and landscape characteristics on farmland bird abundance and diversity in typical cereal cultivation-dominated southern Finnish agricultural environments. The extensive data used were collected by territory mapping. My two main study species were the drastically declined ortolan bunting (Emberiza hortulana) and the phenomenally increased tree sparrow (Passer montanus); in addition I studied assemblages of 20 species breeding in open arable and edge/bush habitats. In light of my results I discuss whether the Finnish AES take into account the habitat needs of farmland birds, and I provide suggestions for improvement of the future AES. My results show that heterogeneity of both uncultivated and cultivated habitats increases abundance and species richness among farmland birds, but in this respect the amount and diversity of uncultivated habitats are essential. Ditches in particular are a keystone structure for farmland birds in boreal landscapes. Ditches lined by trees or bushes increased ortolan bunting abundance. Loss of that kind of ditches (and clearance of forest and bush patches), reduced breeding ortolan buntings, mainly by decreasing availability of song-posts that are important for the breeding groups of the species. Heterogeneity of uncultivated habitats, most importantly open ditches and the habitat patch richness, increased densities and species richnesses of species assemblages of open arable and edge/bush habitats. Human impact (winter-feeding, nest-boxes) affected favourably the tree sparrow s rapid range expansion in southern Finland, but any habitat types had no significant effects. At the moment the Finnish agri-environmental policy does not conserve farmland ditches as a habitat type. Instead, sub-surface drainage is financially promoted. This is a fatal mistake as far as farmland biodiversity is concerned. In addition to the maintenance of ditches, at least the following aspects should be included more than is done previously in the measures of the future AES: 1) promotion of diverse crop rotation (especially by promoting animal husbandry), 2) maintenance of tree and bush vegetation in islets and along ditches, 3) promotion of organic farming.

The ecology of the mountain hare (Lepus timidus L.) in managed boreal forests: habitat associations at different forest scales

Aims of this thesis This study is part of a larger hare project in Finland, which provides answers to basic ecological questions regarding the mountain hare. This study of the ecology of the mountain hare focuses in particular on different levels of managed boreal forest. The feeding habits and intensity of mountain hares in winter are explored, and the connections between mountain hares versus the forest structure are also studied (e.g. habitat use and the importance of different forest layers for hares). The use of the environment by hares at the landscape level was examined (forest patch structures), and the home ranges of mountain hares were studied. Finally, the productivity and survival rate of mountain hare populations were also studied (discussion e.g. predator effects on hare populations). Conclusions Feeding intensity seemed to be highest in the spring-winter, when home ranges were also largest. Favourable food species are covered by snow in winter and the mobility of hares is highest during late winter. A shortage of suitable food species may be problematic for hares, especially during the winter period. In this study mountain hares preferred a dense shrub layer at local level and deciduous and mixed tree forest over coniferous forest at the landscape level. Food and shelter are vital for hares and the preference for particular habitats may also affect the population dynamics of the mountain hare. It would be possible to improve the quality of food and shelter or at least prevent the most negative habitat changes through forest management. At a local level it is also possible to add supplementary food for hares through the winter period. The intensive clearing of young sapling stands and especially the removal of deciduous shrubs and trees reduces the quality of habitats for the mountain hare. Mountain hares primarily live in forest habitat and it is possible that changes in the forest structure play a crucial role in mountain hare habitat preference. Ecological knowledge of the mountain hare is vital to create habitat structure more suitable for the species. More deciduous trees should be saved in managing forests and the mechanical clearing of the shrub layer should be done carefully.

Habitat use and population dynamics of brown bears (Ursus arctos) in Scandinavia

Large carnivore populations are currently recovering from past extirpation efforts and expanding back into their original habitats. At the same time human activities have resulted in very few wilderness areas left with suitable habitats and size large enough to maintain populations of large carnivores without human contact. Consequently the long-term future of large carnivores depends on their successful integration into landscapes where humans live. Thus, understanding their behaviour and interaction with surrounding habitats is of utmost importance in the development of management strategies for large carnivores. This applies also to brown bears (Ursus arctos) that were almost exterminated from Scandinavia and Finland at the turn of the century, but are now expanding their range with the current population estimates being approximately 2600 bears in Scandinavia and 840 in Finland. This thesis focuses on the large-scale habitat use and population dynamics of brown bears in Scandinavia with the objective to develop modelling approaches that support the management of bear populations. Habitat analysis shows that bear home ranges occur mainly in forested areas with a low level of human influence relative to surrounding areas. Habitat modelling based on these findings allows identification and quantification of the potentially suitable areas for bears in Scandinavia. Additionally, this thesis presents novel improvements to home range estimation that enable realistic estimates of the effective area required for the bears to establish a home range. This is achieved through fitting to the radio-tracking data to establish the amount of temporal autocorrelation and the proportion of time spent in different habitat types. Together these form a basis for the landscape-level management of the expanding population. Successful management of bears requires also assessment of the consequences of harvest on the population viability. An individual-based simulation model, accounting for the sexually selected infanticide, was used to investigate the possibility of increasing the harvest using different hunting strategies, such as trophy harvest of males. The results indicated that the population can sustain twice the current harvest rate. However, harvest should be changed gradually while carefully monitoring the population growth as some effects of increased harvest may manifest themselves only after a time-delay. The results and methodological improvements in this thesis can be applied to the Finnish bear population and to other large carnivores. They provide grounds for the further development of spatially-realistic management-oriented models of brow bear dynamics that can make projections of the future distribution of bears while accounting for the development of human activities.

Co-Signaling by Neurotrophic Factors and the Extracellular Matrix for Axonal Growth and Neuronal Survival

Neurotrophic factors (NTFs) and the extracellular matrix (ECM) are important regulators of axonal growth and neuronal survival in mammalian nervous system. Understanding of the mechanisms of this regulation is crucial for the development of posttraumatic therapies and drug intervention in the injured nervous system. NTFs act as soluble, target-derived extracellular regulatory molecules for a wide range of physiological functions including axonal guidance and the regulation of programmed cell death in the nervous system. The ECM determines cell adhesion and regulates multiple physiological functions via short range cell-matrix interactions. The present work focuses on the mechanisms of the action of NTFs and the ECM on axonal growth and survival of cultured sensory neurons from dorsal root ganglia (DRG). We first examined signaling mechanisms of the action of the glial cell line-derived neurotrophic factor (GDNF) family ligands (GFLs) on axonal growth. GDNF, neurturin (NRTN) and artemin (ART) but not persephin (PSPN) promoted axonal initiation in cultured DRG neurons from young adult mice. This effect required Src family kinase (SFK) activity. In neurons from GFRalpha2-deficient mice, NRTN did not significantly promote axonal initiation. GDNF and NRTN induced extensive lamellipodia formation on neuronal somata and growth cones. This study suggested that GDNF, NRTN and ARTN may serve as stimulators of nerve regeneration under posttraumatic conditions. Consequently we studied the convergence of signaling pathways induced by NTFs and the ECM molecule laminin in the intracellular signaling network that regulates axonal growth. We demonstrated that co-stimulation of DRG neurons with NTFs (GDNF, NRTN or nerve growth factor (NGF)) and laminin leads to axonal growth that requires activation of SFKs. A different, SFK-independent signaling pathway evoked axonal growth on laminin in the absence of the NTFs. In contrast, axonal branching was regulated by SFKs both in the presence and in the absence of NGF. We proposed and experimentally verified a Boolean model of the signaling network triggered by NTFs and laminin. Our results put forward an approach for predictable, Boolean logics-driven pharmacological manipulation of a complex signaling network. Finally we found that N-syndecan, the receptor for the ECM component HB-GAM was required for the survival of neonatal sensory neurons in vitro. We demonstrated massive cell death of cultured DRG neurons from mice deficient in the N-syndecan gene as compared to wild type controls. Importantly, this cell death could not be prevented by NGF the neurotrophin which activates multiple anti-apoptotic cascades in DRG neurons. The survival deficit was observed during first postnatal week. By contrast, DRG neurons from young adult N-syndecan knock-out mice exhibited normal survival. This study identifies a completely new syndecan-dependent type of signaling that regulates cell death in neurons.

Dynamics of Finnish starlings in 1951-2005 : from monitoring to population modelling

Farmland bird species have been declining in Europe. Many declines have coincided with general intensification of farming practices. In Finland, replacement of mixed farming, including rotational pastures, with specialized cultivation has been one of the most drastic changes from the 1960s to the 1990s. This kind of habitat deterioration limits the persistence of populations, as has been previously indicated from local populations. Integrated population monitoring, which gathers species-specific information of population size and demography, can be used to assess the response of a population to environment changes also at a large spatial scale. I targeted my analysis at the Finnish starling (Sturnus vulgaris). Starlings are common breeders in farmland habitats, but severe declines of local populations have been reported from Finland in the 1970s and 1980s and later from other parts of Europe. Habitat deterioration (replacement of pasture and grassland habitats with specialized cultivation areas) limits reproductive success of the species. I analysed regional population data in order to exemplify the importance of agricultural change to bird population dynamics. I used nestling ringing and nest-card data from 1951 to 2005 in order to quantify population trends and per capita reproductive success within several geographical regions (south/north and west/east aspects). I used matrix modelling, acknowledging age-specific survival and fecundity parameters and density-dependence, to model population dynamics. Finnish starlings declined by 80% from the end of the 1960s up to the end of the 1980s. The observed patterns and the model indicated that the population decline was due to the decline of the carrying capacity of farmland habitats. The decline was most severe in north Finland where populations largely become extinct. However, habitat deterioration was most severe in the southern breeding areas. The deteriorations in habitat quality decreased reproduction, which finally caused the decline. I suggest that poorly-productive northern populations have been partly maintained by immigration from the highly-productive southern populations. As the southern populations declined, ceasing emigration caused the population extinction in north. This phenomenon was explained with source sink population dynamics, which I structured and verified on the basis of a spatially explicit simulation model. I found that southern Finnish starling population exhibits ten-year cyclic regularity, a phenomenon that can be explained with delayed density-dependence in reproduction.

The effects of habitat edges and trampling intensity on vegetation in urban forests

Although changes in urban forest vegetation have been documented in previous Finnish studies, the reasons for these changes have not been studied explicitly. Especially, the consequences of forest fragmentation, i.e. the fact that forest edges receive more solar radiation, wind and air-borne nutrients than interiors have been ignored. In order to limit the change in urban forest vegetation we need to know why it occurs. Therefore, the effects of edges and recreational use of urban forests on vegetation were investigated together in this thesis to reveal the relative strengths of these effects and to provide recommendations for forest management. Data were collected in the greater Helsinki area (in the cities of Helsinki, Vantaa and Espoo, and in the municipalities of Sipoo and Tuusula) and in the Lahti region (in the city of Lahti and in the municipality of Hollola) by means of systematic and randomized vegetation and soil sampling and tree measurements. Sample plots were placed from the forest edges to the interiors to investigate the effects of forest edges, and on paths of different levels of wear and off these paths to investigate the effects of trampling. The natural vegetation of mesic and sub-xeric forest site types studied was sensitive both to the effects of the edge and to trampling. The abundances of dwarf shrubs and bryophytes decreased, while light- and nitrogen-demanding herbs and grasses - and especially Sorbus aucuparia – were favoured at the edges and next to the paths. Results indicated that typical forest site types at the edges are changing toward more nitrophilic vegetation communities. Covers of the most abundant forest species decreased considerably – even tens of percentages – from interiors to the edges indicating strong edge effects. These effects penetrated at least up to 50 m from the forest edges into the interiors, especially at south to west facing open edges. The effects of trampling were pronounced on paths and even low levels of trampling decreased the abundances of certain species considerably. The effects of trampling extended up to 8 m from path edges. Results showed that the fragmentation of urban forest remnants into small and narrow patches should be avoided in order to maintain natural forest understorey vegetation in the urban setting. Thus, urban forest fragments left within urban development should be at least 3 ha in size, and as circular as possible. Where the preservation of representative original forest interior vegetation is a management aim, closed edges with conifers can act as an effective barrier against solar radiation, wind and urban load, thereby restricting the effects of the edge. Tree volume at the edge should be at least 225-250 m3 ha-1 and the proportion of conifers (especially spruce) 80% or more of the tree species composition. Closed, spruce-dominated edges may also prevent the excessive growth of S. aucuparia saplings at urban forest edges. In addition, closed edges may guide people’s movements to the maintained paths, thus preventing the spontaneous creation of dense path networks. In urban areas the effects of edges and trampling on biodiversity may be considerable, and are important to consider when the aim of management is to prevent the development of homogeneous herb-grass dominated vegetation communities, as was observed at the investigated edges.

Novel matrix metalloproteinases in intestinal inflammation and in cancer of the gastrointestinal tract

Matrix metalloproteinases (MMPs) comprise a family of 23 zinc-dependent human endopeptidases that can degrade virtually all components of the extracellular matrix (ECM). They are classified into eight subgroups according to their structure and into six subgroups based on their substrate-specificity. MMPs have been implicated in inflammation, tissue destruction, cell migration, arthritis, vascular remodeling, angiogenesis, and tumor growth and invasion. MMPs are inhibited by their natural inhibitors, tissue inhibitors of metalloproteinases (TIMPs). Different MMPs function in the same tasks depending on the tissue or cancer subtype. I investigated the role of recently discovered MMPs, especially MMPs-19 and -26, in intestinal inflammation, in intestinal and cutaneous wound healing, and in intestinal cancer. Several MMPs and TIMPs were studied to determine their exact location at tissue level and to obtain information on possible functions of MMPs in such tissues and diseases as the healthy intestine, inflammatory bowel disease (IBD), neonatal necrotizing enterocolitis (NEC), pyoderma gangrenosum (PG), and colorectal as well as pancreatic cancers. In latent celiac disease (CD), I attempted to identify markers to predict later onset of CD in children and adolescents. The main methods used were immunohistochemistry, in situ hybridization, and Taqman RT-PCR. My results show that MMP-26 is important for re-epithelialization in intestinal and cutaneous wound healing. In colon and pancreatic cancers, MMP-26 seems to be a marker of invasive potential, although it is not itself expressed at the invasive front. MMP-21 is upregulated in pancreatic cancer and may be associated with tumor differentiation. MMPs-19 and -28 are associated with normal tissue turnover in the intestine, but they disappear in tumor progression as if they were protective markers . MMP-12 is an essential protease in intestinal inflammation and tissue destruction, as seen here in NEC and in previous CD studies. In patients with type 1 diabetes (T1D), MMPs-1, -3, and -12 were upregulated in the intestinal mucosa. Furthermore, MMP-7 was strongly elevated in NEC. In a model of aberrant wound repair, PG, MMPs-8, -9, and 10 and TNFα may promote ECM destruction, while absence of MMP-1 and MMP-26 from keratinocytes retards re-epithelialization. Based on my results, I suggest MMP-26 to be considered a putative marker for poor prognosis in pancreatic and colon cancer. However, since it functions differently in various tissues and tumor subtypes, this use cannot be generalized. Furthermore, MMP-26 is a beneficial marker for wound healing if expressed by migrating epithelial cells. MMP-12 expression in latent CD patients warrants research in a larger patient population to confirm its role as a specific marker for CD in pathologically indistinct cases. MMP-7 should be considered one of the most crucial proteases in NEC-associated tissue destruction; hence, specific inhibitors of this MMP are worth investigating. In PG, TNFα inhibitors are potential therapeutic agents, as shown already in clinical trials. In conclusion, studies of several MMPs in specific diseases and in healthy tissues are needed to elucidate their roles at the tissue level. MMPs and TIMPs are not exclusively destructive or reparative in tissues. They seem to function differently in different tissues. To identify selective MMP inhibitors, we must thoroughly understand the MMP profile (degradome) and their functions in various organs not to interfere with normal reparative functions during wound repair or beneficial host-response effects during cancer initiation and growth.

Matrix Metalloproteinases and their Tissue Inhibitors as Biomarkers in Ulcerative Colitis and Crohn s Disease

Matrix metalloproteinases (MMPs) represent a family of 23 metalloendopeptidases, collectively capable of degrading all components of the extracellular matrix. MMPs have been implicated in several inflammatory processes such as arthritis, atherosclerosis, and even carcinomas. They are also involved in several beneficial activities such as epithelial repair. MMPs are inhibited by endogenous tissue inhibitors of matrix metalloproteinases (TIMP). In this study, MMPs were investigated in intestinal mucosa of inflammatory bowel diseases (IBD), chronic intestinal disorders. The main focus was to characterize mucosal inflammation in the intestine, but also cutaneous pyoderma gangrenosum (PG), to assess similarites with IBD inflammation. MMPs and TIMPs were mainly examined in colonic mucosa, in adult Crohn s disease (CD), and paediatric CD, ulcerative colitis (UC), and indeterminate colitis (IC). Ileal pouch mucosa of proctocolectomized paediatric onset IBD patients was also investigated to characterize pouch mucosa. The focus was on finding specific MMPs that could act as markers to differentiate between different IBD disorders, and MMPs that could be implied as markers for tissue injury, potentially serving as targets for MMP-inhibitors. All examinations were performed using immunohistochemistry. The results show that immunosuppressive agents decrease stromal expression of MMP-9 and -26 that could serve as specific targets for MMP-inhibitors in treating CD. In paediatric colonic inflammation, MMP-10 and TIMP-3 present as molecular markers for IBD inflammation, and MMP-7 for CD. MMP expression in the the pouch mucosa could not be classified as strictly IBD- or non-IBD-like. For the first time, this study describes the expression of MMP-3, -7, -9, -12, and TIMP-2 and -3 in pouch mucosa. The MMP profile in PG bears resemblance to both intestinal IBD inflammation and cutaneous inflammation. Based on the results, MMPs and their inhibitors emerge as promising tools in the differential diagnosis of IBD and characterization of the disease subtype, although further research is necessary. Furthermore, the expression of several MMPs in pouch has been described for the first time. While further research is warranted, the findings contribute to a better understanding of events occurring in IBD mucosa, as well as pyoderma gangrenosum.

Matrix metalloproteinases as biomarkers in premalignant and malignant tumors of the human skin

The incidence of non-melanoma skin cancer is increasing worldwide. Basal cell carcinoma followed by squamous cell carcinoma and malignant melanoma are the most frequent skin tumors. Immunosuppressed patients have an increased risk of neoplasia, of which non-melanoma skin cancer is the most common. Matrix metalloproteinases (MMPs) are proteolytic enzymes that collectively are capable of degrading virtually all components of the extracellular matrix. MMPs can also process substrates distinct from extracellular matrix proteins and influence cell proliferation, differentiation, angiogenesis, and apoptosis. MMP activity is regulated by their natural inhibitors, tissue inhibitors of metallopro-teinases (TIMPs). In this study, the expression patterns of MMPs, TIMPs, and certain cancer-related molecules were investigated in premalignant and malignant lesions of the human skin. As methods were used immunohistochemisty, in situ hybridization, and reverse transcriptase polymerase chain reaction (RT-PCR) from the cell cultures. Our aim was to evaluate the expression pattern of MMPs in extramammary Paget's disease in order to find markers for more advanced tumors, as well as to shed light on the origin of this rare neoplasm. Novel MMPs -21, -26, and -28 were studied in melanoma cell culture, in primary cutaneous melanomas, and their sentinel nodes. The MMP expression profile in keratoacanthomas and well-differentiated squamous cell carcinomas was analyzed to find markers to differentiate benign keratinocyte hyperproliferation from malignantly transformed cells. Squamous cell carcinomas of immunosuppressed organ transplant recipients were compared to squamous cell carcinomas of matched immunocompetent controls to investigate the factors explaining their more aggressive nature. We found that MMP-7 and -19 proteins are abundant in extramammary Paget's disease and that their presence may predict an underlying adenocarcinoma in these patients. In melanomas, MMP-21 was upregulated in early phases of melanoma progression, but disappeared from the more aggressive tumors with lymph node metastases. The presence of MMP-13 in primary melanomas and lymph node metastases may relate to more aggressive disease. In keratoacanthomas, the expression of MMP-7 and -9 is rare and therefore should raise a suspicion of well-differentiated squamous cell carcinomas. Furthermore, MMP-19 and p16 were observed in benign keratinocyte hyperproliferation of keratoacanthomas, whereas they were generally lost from malignant keratinocytes of SCCs. MMP-26 staining was significantly stronger in squamous cell carcinomas and Bowen s disease samples of organ transplant recipients and it may contribute to the more aggressive nature of squamous cell carcinomas in immunosuppressed patients. In addition, the staining for MMP-9 was significantly stronger in macrophages surrounding the tumors of the immunocompetent group and in neutrophils of those patients on cyclosporin medication. In conclusion, based on our studies, MMP-7 and -19 might serve as biomarkers for more aggressive extramammary Paget's disease and MMP-21 for malignant transformation of melanocytes. MMP -7, -9, and -26, however, could play an important role in the pathobiology of keratinocyte derived malignancies.