An R package for simulating metapopulation dynamics and range expansion under environmental change

The metapopulation paradigm is central in ecology and conservation biology to understand the dynamics of spatially-structured populations in fragmented landscapes. Metapopulations are often studied using simulation modelling, and there is an increasing demand of user-friendly software tools to simulate metapopulation responses to environmental change. Here we describe the MetaLandSim R package, mwhich integrates ideas from metapopulation and graph theories to simulate the dynamics of real and virtual metapopulations. The package offers tools to (i) estimate metapopulation parameters from empirical data, (ii) to predict variation in patch occupancy over time in static and dynamic landscapes, either real or virtual, and (iii) to quantify the patterns and speed of metapopulation expansion into empty landscapes. MetaLandSim thus provides detailed information on metapopulation processes, which can be easily combined with land use and climate change scenarios to predict metapopulation dynamics and range expansion for a variety of taxa and ecological systems.

Using complex network metrics to predict the persistence of metapopulations with asymmetric connectivity patterns

Almost all metapopulation modelling assumes that connectivity between patches is only a function of distance, and is therefore symmetric. However, connectivity will not depend only on the distance between the patches, as some paths are easy to traverse, while others are difficult. When colonising organisms interact with the heterogeneous landscape between patches, connectivity patterns will invariably be asymmetric. There have been few attempts to theoretically assess the effects of asymmetric connectivity patterns on the dynamics of metapopulations. In this paper, we use the framework of complex networks to investigate whether metapopulation dynamics can be determined by directly analysing the asymmetric connectivity patterns that link the patches. Our analyses focus on “patch occupancy” metapopulation models, which only consider whether a patch is occupied or not. We propose three easily calculated network metrics: the “asymmetry” and “average path strength” of the connectivity pattern, and the “centrality” of each patch. Together, these metrics can be used to predict the length of time a metapopulation is expected to persist, and the relative contribution of each patch to a metapopulation’s viability. Our results clearly demonstrate the negative effect that asymmetry has on metapopulation persistence. Complex network analyses represent a useful new tool for understanding the dynamics of species existing in fragmented landscapes, particularly those existing in large metapopulations.

Optimising weed control by considering disturbance

Mounting scientific evidence suggests newly imposed disturbance and/or alterations to existing disturbances facilitate invasion. Several empirical studies have explored the role of disturbance in invasion, but little work has been done to fit current understanding into a format useful for practical control efforts. We are working towards addressing this shortcoming by developing a metapopulation model couched in a decision theory framework. This approach has allowed us to investigate how incorporating the negative effects of disturbance on native vegetation into decision-making can change optimal control measures. In this paper, we present some preliminary results.

Extinction risk in cloud forest fragments under climate change and habitat loss

Aim: To quantify the consequences of major threats to biodiversity, such as climate and land-use change, it is important to use explicit measures of species persistence, such as extinction risk. The extinction risk of metapopulations can be approximated through simple models, providing a regional snapshot of the extinction probability of a species. We evaluated the extinction risk of three species under different climate change scenarios in three different regions of the Mexican cloud forest, a highly fragmented habitat that is particularly vulnerable to climate change. Location: Cloud forests in Mexico. Methods: Using Maxent, we estimated the potential distribution of cloud forest for three different time horizons (2030, 2050 and 2080) and their overlap with protected areas. Then, we calculated the extinction risk of three contrasting vertebrate species for two scenarios: (1) climate change only (all suitable areas of cloud forest through time) and (2) climate and land-use change (only suitable areas within a currently protected area), using an explicit patch-occupancy approximation model and calculating the joint probability of all populations becoming extinct when the number of remaining patches was less than five. Results: Our results show that the extent of environmentally suitable areas for cloud forest in Mexico will sharply decline in the next 70 years. We discovered that if all habitat outside protected areas is transformed, then only species with small area requirements are likely to persist. With habitat loss through climate change only, high dispersal rates are sufficient for persistence, but this requires protection of all remaining cloud forest areas. Main conclusions: Even if high dispersal rates mitigate the extinction risk of species due to climate change, the synergistic impacts of changing climate and land use further threaten the persistence of species with higher area requirements. Our approach for assessing the impacts of threats on biodiversity is particularly useful when there is little time or data for detailed population viability analyses. © 2013 John Wiley & Sons Ltd.

Stochastic dynamics of interacting haematopoietic stem cell niche lineages

Since we still know very little about stem cells in their natural environment, it is useful to explore their dynamics through modelling and simulation, as well as experimentally. Most models of stem cell systems are based on deterministic differential equations that ignore the natural heterogeneity of stem cell populations. This is not appropriate at the level of individual cells and niches, when randomness is more likely to affect dynamics. In this paper, we introduce a fast stochastic method for simulating a metapopulation of stem cell niche lineages, that is, many sub-populations that together form a heterogeneous metapopulation, over time. By selecting the common limiting timestep, our method ensures that the entire metapopulation is simulated synchronously. This is important, as it allows us to introduce interactions between separate niche lineages, which would otherwise be impossible. We expand our method to enable the coupling of many lineages into niche groups, where differentiated cells are pooled within each niche group. Using this method, we explore the dynamics of the haematopoietic system from a demand control system perspective. We find that coupling together niche lineages allows the organism to regulate blood cell numbers as closely as possible to the homeostatic optimum. Furthermore, coupled lineages respond better than uncoupled ones to random perturbations, here the loss of some myeloid cells. This could imply that it is advantageous for an organism to connect together its niche lineages into groups. Our results suggest that a potential fruitful empirical direction will be to understand how stem cell descendants communicate with the niche and how cancer may arise as a result of a failure of such communication.

Will biological control of Lantana camara ever succeed? Patterns, processes & prospects

Despite biocontrol research spanning over 100 years, the hybrid weed, commonly referred to as Lantana camara, is not under adequate control. Host specificity and varietal preference of released agents, climatic suitability of a region for released agents, number of agents introduced and range or area of infestation appear to play a role in limiting biocontrol success. At least one of 41 species of mainly leaf- or flower-feeding insects has been introduced, or spread, to 41 of the 70 countries or regions where lantana occurs. Over half (26) of these species have established, achieving varying levels of herbivory and presumably some degree of control. Accurate taxonomy of the plant and adaptation of potential agents to the host plant are some of the better predictors of at least establishment success. Retrospective analysis of the hosts of introduced biocontrol agents for L. camara show that a greater proportion of agents that were collected from L. camara or Lantana urticifolia established, than agents that were collected from other species of Lantana. Of the introduced agents that had established and were oligophagous, 18 out of 22 established. The proportion of species establishing, declined with the number of species introduced. However, there was no trend when oceanic islands were treated separately from mainland areas and the result is likely an artefact of how introductions have changed over time. A calculated index of the degree of herbivory due to agents known to have caused some damage per country, was not related to land area infested with lantana for mainlands nor for oceanic islands. However, the degree of herbivory is much higher on islands than mainlands. This difference between island and mainland situations may reflect population dynamics in patchy or metapopulation landscapes. Basic systematic studies of the host remain crucial to successful biocontrol, especially of hybrid weeds like L. camara. Potential biocontrol agents should be monophages collected from the most closely related species to the target weed or be phytophages that attack several species of lantana. Suitable agents should be released in the most ideal ecoclimatic area. Since collection of biocontrol agents has been limited to a fraction of the known number of phytophagous species available, biocontrol may be improved by targeting insects that feed on stems and roots, as well as the agents that feed on leaves and flowers.

Spatial subdivision among assemblages of Spanish mackerel, Scomberomorus commerson (Pisces: Scombridae) across northern Australia: implications for fisheries management

Aim: This study investigated the use of stable δ13C and δ18O isotopes in the sagittal otolith carbonate of narrow-barred Spanish mackerel, Scomberomorus commerson, as indicators of population structure across Australia. Location: Samples were collected from 25 locations extending from the lower west coast of Western Australia (30°), across northern Australian waters, and to the east coast of Australia (18°) covering a coastline length of approximately 9500 km, including samples from Indonesia. Methods: The stable δ13C and δ18O isotopes in the sagittal otolith carbonate of S. commerson were analysed using standard mass spectrometric techniques. The isotope ratios across northern Australian subregions were subjected to an agglomerative hierarchical cluster analysis to define subregions. Isotope ratios within each of the subregions were compared to assess population structure across Australia. Results: Cluster analysis separated samples into four subregions: central Western Australia, north Western Australia, northern Australia and the Gulf of Carpentaria and eastern Australia. Isotope signatures for fish from a number of sampling sites from across Australia and Indonesia were significantly different, indicating population separation. No significant differences were found in otolith isotope ratios between sampling times (no temporal variation). Main conclusions: Significant differences in the isotopic signatures of S. commerson demonstrate that there is unlikely to be any substantial movement of fish among these spatially discrete adult assemblages. The lack of temporal variation among otolith isotope ratios indicates that S. commerson populations do not undergo longshore spatial shifts in distribution during their life history. The temporal persistence of spatially explicit stable isotopic signatures indicates that, at these spatial scales, the population units sampled comprise functionally distinct management units or separate ‘stocks’ for many of the purposes of fisheries management. The spatial subdivision evident among populations of S. commerson across northern and western Australia indicates that it may be advantageous to consider S. commerson population dynamics and fisheries management from a metapopulation perspective (at least at the regional level).

The decline of northern malaria and population dynamics of Plamodium vivax

Europe was declared malaria free in 1975. The disappearance of malaria has traditionally been attributed to numerous deliberate actions like vector control, the screening of houses, more efficient medication etc. Malaria, however, disappeared from many countries like Finland before any counter measures had even started. The aim of this thesis is to study the population ecology of P. vivax and its interaction with the human host and the vector. By finding the factors that attributed to the extinction of vivax malaria it might be possible to improve the modern strategy against P. vivax. The parasite was studied with data from Finland, which provides the longest time series (1749-2008) of malaria statistics in the world. The malaria vectors, Anopheles messeae and A. beklemishevi are still common species in the country. The eradication of vivax malaria is difficult because the parasite has a dormant stage that can cause a relapse long after a primary infection. It was now shown that P. vivax is able to detect the presence of a potential vector. A dormant stage is triggered even from a bite of an uninfected Anopheles mosquito. This optimizes the chances for the Plasmodium to reach a mosquito vector for sexual reproduction. The longevity of the dormant stage could be shown to be at least nine years. The parasite spends several years in its human host and the behaviour of the human carrier had a profound impact on the decline of the disease in Finland. Malaria spring epidemics could be explained by a previous warm summer. Neither annual nor summer mean temperature had any impact on the long term malaria trend. Malaria disappeared slowly from Finland without mosquito control. The sociological change from extended families to nuclear families led to decreased household size. The decreased household size correlated strongly with the decline of malaria. That led to an increased isolation of the subpopulations of P. vivax. Their habitat consisted of the bedrooms in which human carriers slept together with the overwintering vectors. The isolation of the parasite ultimately led to the extinction of vivax malaria. Metapopulation models adapted to local conditions should therefore be implemented as a tool for settlement planning and socio-economic development and become an integrated part of the fight against malaria.

Population Expansion and Genetic Structure in Carcharhinus brevipinna in the Southern Indo-Pacific

Background:Quantifying genetic diversity and metapopulation structure provides insights into the evolutionary history of a species and helps develop appropriate management strategies. We provide the first assessment of genetic structure in spinner sharks (Carcharhinus brevipinna), a large cosmopolitan carcharhinid, sampled from eastern and northern Australia and South Africa. Methods and Findings:Sequencing of the mitochondrial DNA NADH dehydrogenase subunit 4 gene for 430 individuals revealed 37 haplotypes and moderately high haplotype diversity (h = 0.6770 ±0.025). While two metrics of genetic divergence (ΦST and FST) revealed somewhat different results, subdivision was detected between South Africa and all Australian locations (pairwise ΦST, range 0.02717–0.03508, p values ≤ 0.0013; pairwise FST South Africa vs New South Wales = 0.04056, p = 0.0008). Evidence for fine-scale genetic structuring was also detected along Australia’s east coast (pairwise ΦST = 0.01328, p < 0.015), and between south-eastern and northern locations (pairwise ΦST = 0.00669, p < 0.04).Conclusions: The Indian Ocean represents a robust barrier to contemporary gene flow in C. brevipinna between Australia and South Africa. Gene flow also appears restricted along a continuous continental margin in this species, with data tentatively suggesting the delineation of two management units within Australian waters. Further sampling, however, is required for a more robust evaluation of the latter finding. Evidence indicates that all sampled populations were shaped by a substantial demographic expansion event, with the resultant high genetic diversity being cause for optimism when considering conservation of this commercially-targeted species in the southern Indo-Pacific.

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.

Individuals on the Move : Body Condition Dependent Dispersal and Quasi-Local Competition in Metapopulations

Individual movement is very versatile and inevitable in ecology. In this thesis, I investigate two kinds of movement body condition dependent dispersal and small-range foraging movements resulting in quasi-local competition and their causes and consequences on the individual, population and metapopulation level. Body condition dependent dispersal is a widely evident but barely understood phenomenon. In nature, diverse relationships between body condition and dispersal are observed. I develop the first models that study the evolution of dispersal strategies that depend on individual body condition. In a patchy environment where patches differ in environmental conditions, individuals born in rich (e.g. nutritious) patches are on average stronger than their conspecifics that are born in poorer patches. Body condition (strength) determines competitive ability such that stronger individuals win competition with higher probability than weak individuals. Individuals compete for patches such that kin competition selects for dispersal. I determine the evolutionarily stable strategy (ESS) for different ecological scenarios. My models offer explanations for both dispersal of strong individuals and dispersal of weak individuals. Moreover, I find that within-family dispersal behaviour is not always reflected on the population level. This supports the fact that no consistent pattern is detected in data on body condition dependent dispersal. It also encourages the refining of empirical investigations. Quasi-local competition defines interactions between adjacent populations where one population negatively affects the growth of the other population. I model a metapopulation in a homogeneous environment where adults of different subpopulations compete for resources by spending part of their foraging time in the neighbouring patches, while their juveniles only feed on the resource in their natal patch. I show that spatial patterns (different population densities in the patches) are stable only if one age class depletes the resource very much but mainly the other age group depends on it.

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.

Effects of environmental variation on ecological and evolutionary dynamics

Environmental variation is a fact of life for all the species on earth: for any population of any particular species, the local environmental conditions are liable to vary in both time and space. In today's world, anthropogenic activity is causing habitat loss and fragmentation for many species, which may profoundly alter the characteristics of environmental variation in remaining habitat. Previous research indicates that, as habitat is lost, the spatial configuration of remaining habitat will increasingly affect the dynamics by which populations are governed. Through the use of mathematical models, this thesis asks how environmental variation interacts with species properties to influence population dynamics, local adaptation, and dispersal evolution. More specifically, we couple continuous-time continuous-space stochastic population dynamic models to landscape models. We manipulate environmental variation via parameters such as mean patch size, patch density, and patch longevity. Among other findings, we show that a mixture of high and low quality habitat is commonly better for a population than uniformly mediocre habitat. This conclusion is justified by purely ecological arguments, yet the positive effects of landscape heterogeneity may be enhanced further by local adaptation, and by the evolution of short-ranged dispersal. The predicted evolutionary responses to environmental variation are complex, however, since they involve numerous conflicting factors. We discuss why the species that have high levels of local adaptation within their ranges may not be the same species that benefit from local adaptation during range expansion. We show how habitat loss can lead to either increased or decreased selection for dispersal depending on the type of habitat and the manner in which it is lost. To study the models, we develop a recent analytical method, Perturbation expansion, to enable the incorporation of environmental variation. Within this context, we use two methods to address evolutionary dynamics: Adaptive dynamics, which assumes mutations occur infrequently so that the ecological and evolutionary timescales can be separated, and via Genotype distributions, which assume mutations are more frequent. The two approaches generally lead to similar predictions yet, exceptionally, we show how the evolutionary response of dispersal behaviour to habitat turnover may qualitatively depend on the mutation rate.

Dispersal and related life history traits in the Glanville fritillary butterfly

Most studies of life history evolution are based on the assumption that species exist at equilibrium and spatially distinct separated populations. In reality, this is rarely the case, as populations are often spatially structured with ephemeral local populations. Therefore, the characteristics of metapopulations should be considered while studying factors affecting life history evolution. Theoretical studies have examined spatial processes shaping the evolution of life history traits to some extent, but there is little empirical data and evidence to investigate model predictions. In my thesis I have tried to bridge the gap between theoretical and empirical studies by using the well-known Glanville fritillary (Melitaea cinxia) metapopulation as a model system. The long-term persistence of classic metapopulations requires sufficient dispersal to establish new local populations to compensate for local extinctions. Previous studies on the Glanville fritillary have shown that females establishing new populations are not a random sample from the metapopulation, but they are in fact more dispersive than females in old populations. Many other life-history traits, such as body size, fecundity and lifespan, may be related to dispersal rate. Therefore, I examined a range of correlated traits for their evolutionary and ecological consequences. I was particularly interested in how the traits vary under natural environmental conditions, hence all studies were conducted in a large (32 x 26 m) outdoor population cage built upon a natural habitat patch. Individuals for the experiments were sampled from newly-established and old populations within a large metapopulation. Results show that females originating from newly-established populations had higher within-habitat patch mobility than females from old populations. I showed that dispersal rate is heritable and that flight activity is related to variation in a gene encoding the glycolytic enzyme phosphoglucose isomerase. Both among-individual and among-population variation in dispersal are correlated with the reproductive performance of females, though I found no evidence for a trade-off between dispersal and fecundity in terms of lifetime egg production or clutch size. Instead, the results suggest that highly dispersive females from newly-established populations have a shorter lifespan than females from old populations, and that dispersive females may pay a cost in terms of reduced lifetime reproductive success due to increased time spent outside habitat patches. In summary, the results of this thesis show that genotype-dependent dispersal rate correlates with other life history traits in the Glanville fritillary, and that the rapid turnover of local populations (extinctions and re-colonisations) is likely to be the mechanism that maintains phenotypic variation in many life history traits at the metapopulation level.

Ecology, population genetics and conservation of the African violet

The ongoing rapid fragmentation of tropical forests is a major threat to global biodiversity. This is because many of the tropical forests are so-called biodiversity 'hotspots', areas that host exceptional species richness and concentrations of endemic species. Forest fragmentation has negative ecological and genetic consequences for plant survival. Proposed reasons for plant species' loss in forest fragments are, e.g., abiotic edge effects, altered species interactions, increased genetic drift, and inbreeding depression. To be able to conserve plants in forest fragments, the ecological and genetic processes that threaten the species have to be understood. That is possible only after obtaining adequate information on their biology, including taxonomy, life history, reproduction, and spatial and genetic structure of the populations. In this research, I focused on the African violet (genus Saintpaulia), a little-studied conservation flagship from the Eastern Arc Mountains and Coastal Forests hotspot of Tanzania and Kenya. The main objective of the research was to increase understanding of the life history, ecology and population genetics of Saintpaulia that is needed for the design of appropriate conservation measures. A further aim was to provide population-level insights into the difficult taxonomy of Saintpaulia. Ecological field work was conducted in a relatively little fragmented protected forest in the Amani Nature Reserve in the East Usambara Mountains, in northeastern Tanzania, complemented by population genetic laboratory work and ecological experiments in Helsinki, Finland. All components of the research were conducted with Saintpaulia ionantha ssp. grotei, which forms a taxonomically controversial population complex in the study area. My results suggest that Saintpaulia has good reproductive performance in forests with low disturbance levels in the East Usambara Mountains. Another important finding was that seed production depends on sufficient pollinator service. The availability of pollinators should thus be considered in the in situ management of threatened populations. Dynamic population stage structures were observed suggesting that the studied populations are demographically viable. High mortality of seedlings and juveniles was observed during the dry season but this was compensated by ample recruitment of new seedlings after the rainy season. Reduced tree canopy closure and substrate quality are likely to exacerbate seedling and juvenile mortality, and, therefore, forest fragmentation and disturbance are serious threats to the regeneration of Saintpaulia. Restoration of sufficient shade to enhance seedling establishment is an important conservation measure in populations located in disturbed habitats. Long-term demographic monitoring, which enables the forecasting of a population s future, is also recommended in disturbed habitats. High genetic diversities were observed in the populations, which suggest that they possess the variation that is needed for evolutionary responses in a changing environment. Thus, genetic management of the studied populations does not seem necessary as long as the habitats remain favourable for Saintpaulia. The observed high levels of inbreeding in some of the populations, and the reduced fitness of the inbred progeny compared to the outbred progeny, as revealed by the hand-pollination experiment, indicate that inbreeding and inbreeding depression are potential mechanisms contributing to the extinction of Saintpaulia populations. The relatively weak genetic divergence of the three different morphotypes of Saintpaulia ionantha ssp. grotei lend support to the hypothesis that the populations in the Usambara/lowlands region represent a segregating metapopulation (or metapopulations), where subpopulations are adapting to their particular environments. The partial genetic and phenological integrity, and the distinct trailing habit of the morphotype 'grotei' would, however, justify its placement in a taxonomic rank of its own, perhaps in a subspecific rank.