Behavioural, population, and genetic processes affecting metapopulation dynamics of the Glanville fritillary butterfly

In my thesis I have been studying the effects of population fragmentation and extinction-recolonization dynamics on genetic and evolutionary processes in the Glanville fritillary butterfly (Melitaea cinxia). By conducting crosses within and among newly-colonized populations and using several fitness measures, I found a strong decrease in fitness following colonization by a few related individuals, and a strong negative relationship between parental relatedness and offspring fitness. Thereafter, I was interested in determining the number and relatedness of individuals colonizing new populations, which I did using a set of microsatellites I had previously developed for this species. Additionally, I am interested in the evolution of key life-history traits. By following the lifetime reproductive success of males emerging at different times in a semi-natural setup, I demonstrated that protandry is adaptive in males, and I was able to rule out, for M. cinxia, alternative incidental hypotheses evoked to explain the evolution of protandry in insects. Finally, in work I did together with Prof. Hanna Kokko, I am proposing bet-hedging as a new mechanism that could explain the evolution of polyandry in M. cinxia.

The life cycle and genetic structure of the red alga Furcellaria lumbricalis on a salinity gradient

The life cycle and genetic diversity of the red alga Furcellaria lumbricalis (Hudson) Lamouroux were investigated in 15 populations in northern Europe. The occurrence of different life cycle phases and seasonality of reproduction were studied in four brackish populations in the northern Baltic Sea. Furthermore, a new method, based on genome screening with ISSR markers combined with a restriction-ligation method, was developed to discover microsatellite markers for population genetic analyses. The mitochondrial DNA cox2-3 spacer sequence and four microsatellite markers were used to examine the genetic diversity and differentiation of red algal populations in northern Europe. In addition, clonality and small-scale genetic structure of one Irish and four Baltic Sea populations were studied with microsatellite markers. It was discovered that at the low salinities of the northern Baltic Sea, only tetrasporophytes and males were present in the populations of F. lumbricalis and that winter was the main season for tetrasporangial production. Furthermore, the population occurring at the lowest salinity (3.6 practical salinity units, psu) did not produce spores. The size of the tetraspores was smaller in the Baltic Sea populations than that in the Irish population, and there were more deformed spores in the Baltic Sea populations than in the Irish populations. Studies with microsatellite markers indicated that clonality is a common phenomenon in the Baltic Sea populations of F. lumbricalis, although the proportion of clonal individuals varied among populations. Some genetic divergence occurred within locations both in Ireland and in the northern Baltic Sea. Even though no carpogonia were detected in the field samples during the study, the microsatellite data indicated that sexual reproduction occurs at least occasionally in the northern Baltic Sea. The genetic diversity of F. lumbricalis was highest in Brittany, France. Since no variation was discovered in the mtDNA cox2-3 spacer sequence, which is generally regarded as an informative phylogeographic marker in red algae, it can be assumed that the studied populations probably share the same origin.

Should I stay or should I go? : Reproductive and dispersal strategies of site-specific liverwort species

The area of intensively managed forests, in which required conditions for several liverwort species are seldom found, has expanded over the forest landscape during the last century. Liverworts are very sensitive to habitat changes, because they demand continuously moist microclimate. Consequently, about third of the forest liverworts have been classified as threatened or near threatened in Finland. The general objective of this thesis is to increase knowledge of the reproductive and dispersal strategies of the substrate-specific forest bryophytes. A further aim was to develop recommendations for conservation measures for species inhabiting unstable and stable habitats in forest landscape. Both population ecological and genetic methods have been applied in the research. Anastrophyllum hellerianum inhabits spatially and temporally limited substrate patches, decaying logs, which can be considered as unstable habitats. The results show that asexual reproduction by gemmae is the dominant mode of reproduction, whereas sexual reproduction is considerably infrequent. Unlike previously assumed, not only spores but also the asexual propagules may contribute to long-distance dispersal. The combination of occasional spore production and practically continuous, massive gemma production facilitates dispersal both on a local scale and over long distances, and it compensates for the great propagule losses that take place preceding successful establishment at suitable sites. However, establishment probability of spores may be restricted because of environmental and biological limitations linked to the low success of sexual reproduction. Long-lasting dry seasons are likely to result in a low success of sexual reproduction and decreased release rate of gemmae from the shoots, and consequent fluctuations in population sizes. In the long term, the substratum limitation is likely to restrict population sizes and cause local extinctions, especially in small-sized remnant populations. Contrastingly, larger forest fragments with more natural disturbance dynamics, to which the species is adapted, are pivotal to species survival. Trichocolea tomentella occupies stable spring and mesic habitats in woodland. The relatively small populations are increasingly fragmented with a high risk for extinction for extrinsic reasons. The results show that T. tomentella mainly invests in population persistence by effective clonal growth via forming independent ramets and in competitive ability, and considerably less in sexuality and dispersal potential. The populations possess relatively high levels of genetic diversity regardless of population size and of degree of isolation. Thus, the small-sized populations inhabiting stable habitats should not be neglected when establishing conservation strategies for the species and when considering the habitat protection of small spring sites. Restricted dispersal capacity, also on a relatively small spatial scale, is likely to prevent successful (re-)colonization in the potential habitat patches of recovering forest landscapes. By contrast, random short-range dispersal of detached vegetative fragments within populations at suitable habitat seems to be frequent. Thus, the restoration actions of spring and streamside habitats close to the populations of T. tomentella may contribute to population expansion. That, in turn, decreases the harmful effects of environmental stochasticity.

Phylogeography and Adaptive Divergence of Three-spined Stickleback Populations

Genetic studies on phylogeography and adaptive divergence in Northern Hemisphere fish species such as three-spined stickleback (Gasterosteus aculeatus) provide an excellent opportunity to investigate genetic mechanisms underlying population differentiation. According to the theory, the process of population differentiation results from a complex interplay between random and deterministic processes as well historical factors. The main scope in this thesis was to study how historical factors like the Pleistocene ice ages have shaped the patterns molecular diversity in three-spined stickleback populations in Europe and how this information could be utilized in the conservation genetic context. Furthermore, identifying footprints of natural selection at the DNA level might be used in identifying genes involved in evolutionary change. Overall, the results from phylogeographic studies indicate that the three-spined stickleback has colonized the Atlantic basin relatively recently but constitutes three major evolutionary lineages in Europe. In addition, the colonization of freshwater appears to result from multiple and independent invasions by the marine conspecifics. Molecular data together with morphology suggest that the most divergent freshwater populations are located in the Balkan Peninsula and these populations deserve a special conservation genetic status without warranting further taxonomical classification. In order to investigate the adaptive divergence in Fennoscandian three-spined stickleback populations several approaches were used. First, sequence variability in the Eda-gene, coding for the number of lateral plates, was concordant with the previously observed global pattern. Full plated allele is in high frequencies among marine populations whereas low plated allele dominates in the freshwater populations. Second, a microsatellite based genome scan identified both indications of balancing and directional selection in the three-spined stickleback genome, i.e. loci with unusually similar or unusually different allele frequencies over populations. The directionally selected loci were mainly associated with the adaptation to freshwater. A follow up study conducting a more detailed analysis in a chromosome region containing a putatively selected gene locus identified a fairly large genomic region affected by natural selection. However, this region contained several gene predictions, all of which might be the actual target of natural selection. All in all, the phylogeographic and adaptive divergence studies indicate that most of the genetic divergence has occurred in the freshwater populations whereas the marine populations have remained relatively uniform.

Human genetic variation in the Baltic Sea region: Features of population history and natural selection

In this thesis, the genetic variation of human populations from the Baltic Sea region was studied in order to elucidate population history as well as evolutionary adaptation in this region. The study provided novel understanding of how the complex population level processes of migration, genetic drift, and natural selection have shaped genetic variation in North European populations. Results from genome-wide, mitochondrial DNA and Y-chromosomal analyses suggested that the genetic background of the populations of the Baltic Sea region lies predominantly in Continental Europe, which is consistent with earlier studies and archaeological evidence. The late settlement of Fennoscandia after the Ice Age and the subsequent small population size have led to pronounced genetic drift, especially in Finland and Karelia but also in Sweden, evident especially in genome-wide and Y-chromosomal analyses. Consequently, these populations show striking genetic differentiation, as opposed to much more homogeneous pattern of variation in Central European populations. Additionally, the eastern side of the Baltic Sea was observed to have experienced eastern influence in the genome-wide data as well as in mitochondrial DNA and Y-chromosomal variation – consistent with linguistic connections. However, Slavic influence in the Baltic Sea populations appears minor on genetic level. While the genetic diversity of the Finnish population overall was low, genome-wide and Y-chromosomal results showed pronounced regional differences. The genetic distance between Western and Eastern Finland was larger than for many geographically distant population pairs, and provinces also showed genetic differences. This is probably mainly due to the late settlement of Eastern Finland and local isolation, although differences in ancestral migration waves may contribute to this, too. In contrast, mitochondrial DNA and Y-chromosomal analyses of the contemporary Swedish population revealed a much less pronounced population structure and a fusion of the traces of ancient admixture, genetic drift, and recent immigration. Genome-wide datasets also provide a resource for studying the adaptive evolution of human populations. This study revealed tens of loci with strong signs of recent positive selection in Northern Europe. These results provide interesting targets for future research on evolutionary adaptation, and may be important for understanding the background of disease-causing variants in human populations.

Cytokinin signalling in the regulation of cambial development

Secondary growth of plants is of pivotal importance in terrestrial ecosystems, providing a significant carbon sink in the form of wood. As plant biomass accumulation results largely from the cambial growth, it is surprising that quite little is known about the hormonal or genetic control of this important process in any plant species. The central aim of my thesis studies was to explore the function of cytokinin in the regulation of cambial development. Since their discovery as regulators of plant cell divisions, cytokinins have been assumed to participate in the control of cambial development. Evidence for this action was deduced from hormone treatment experiments, where exogenously applied cytokinin was shown to enhance cambial cell divisions in diverse plant organs and species. In my thesis work, the conservation of cytokinin signalling and homeostasis genes between a herbaceous plant, Arabidopsis, and a hardwood tree species, Populus trichocarpa. Presumably reflecting the ancient origin of cytokinin signalling system, the Populus genome contains orthologs for all Arabidopsis cytokinin signalling and homeostasis genes. Thus, genes belonging to five main families of isopentenyl transferases (IPTs), cytokinin oxidases (CKXs), two-component receptors, histidine containing phosphotransmitters (HPts) and response regulators (RRs) were identified from the Populus genome. Three subfamilies associated with cytokinin signal transduction, the CKI1-like family of two-component receptors, the AHP4-like HPts, and the ARR22-like atypical RRs, were significantly larger in Populus genome than in Arabidopsis. Potential contribution to the extensive secondary development of Populus by the members of these considerably expanded gene families will be discussed. Representatives of all cytokinin signal transduction elements were expressed in the Populus cambial zone, and most of the expressed genes appeared to be slightly more abundant on the phloem side of the meristem. The abundance of cytokinin related genes in the cambium emphasizes the important role of this hormone in the regulation of the extensive secondary growth characteristic of tree species. The function of the pseudo HPts in primary vascular development was studied in Arabidopsis root vasculature. It was demonstrated that the pseudo HPt AHP6 has a role in locally inhibiting cytokinin signalling in the protoxylem position in the Arabidopsis root, thus enabling differentiation of the protoxylem cell file. The possible role of pseudo HPts in cambial development will be discussed. The expression peak of cytokinin signalling genes in the tree cambial zone strongly indicates that cytokinin has a role in the regulation of this meristem function. To address whether cytokinin signalling is required for cambial activity, transgenic Populus trees with modified cytokinin signalling were produced. These trees were expressing a cytokinin catabolic gene from Arabidopsis, CYTOKININ OXIDASE 2, (AtCKX2) under the promoter of a Betula CYTOKININ RECEPTOR 1 (BpCRE1). The pBpCRE1::CKX2 transgenic Populus trees showed a reduced concentration of a biologically active cytokinin, correlating with their impaired cytokinin response. Furthermore, the radial growth of these trees was compromised, as illustrated by a smaller stem diameter than in wild-type trees of the same height. Moreover, the level of cambial cytokinin signalling was down-regulated in these thin-stemmed trees. The reduced signalling correlated with a decreased number of meristematic cambial cells, implicating cytokinin activity as a direct regulator of cambial cell division activity. Together, the results of my study indicate that cytokinins are major hormonal regulators required for cambial development.

Species-based and community-level approaches to conservation prioritization

One major reason for the global decline of biodiversity is habitat loss and fragmentation. Conservation areas can be designed to reduce biodiversity loss, but as resources are limited, conservation efforts need to be prioritized in order to achieve best possible outcomes. The field of systematic conservation planning developed as a response to opportunistic approaches to conservation that often resulted in biased representation of biological diversity. The last two decades have seen the development of increasingly sophisticated methods that account for information about biodiversity conservation goals (benefits), economical considerations (costs) and socio-political constraints. In this thesis I focus on two general topics related to systematic conservation planning. First, I address two aspects of the question about how biodiversity features should be valued. (i) I investigate the extremely important but often neglected issue of differential prioritization of species for conservation. Species prioritization can be based on various criteria, and is always goal-dependent, but can also be implemented in a scientifically more rigorous way than what is the usual practice. (ii) I introduce a novel framework for conservation prioritization, which is based on continuous benefit functions that convert increasing levels of biodiversity feature representation to increasing conservation value using the principle that more is better. Traditional target-based systematic conservation planning is a special case of this approach, in which a step function is used for the benefit function. We have further expanded the benefit function framework for area prioritization to address issues such as protected area size and habitat vulnerability. In the second part of the thesis I address the application of community level modelling strategies to conservation prioritization. One of the most serious issues in systematic conservation planning currently is not the deficiency of methodology for selection and design, but simply the lack of data. Community level modelling offers a surrogate strategy that makes conservation planning more feasible in data poor regions. We have reviewed the available community-level approaches to conservation planning. These range from simplistic classification techniques to sophisticated modelling and selection strategies. We have also developed a general and novel community level approach to conservation prioritization that significantly improves on methods that were available before. This thesis introduces further degrees of realism into conservation planning methodology. The benefit function -based conservation prioritization framework largely circumvents the problematic phase of target setting, and allowing for trade-offs between species representation provides a more flexible and hopefully more attractive approach to conservation practitioners. The community-level approach seems highly promising and should prove valuable for conservation planning especially in data poor regions. Future work should focus on integrating prioritization methods to deal with multiple aspects in combination influencing the prioritization process, and further testing and refining the community level strategies using real, large datasets.

Ecological consequences of genetic modifications - an invasion analysis approach

Biological invasions are considered as one of the greatest threats to biodiversity, as they may lead to disruption and homogenization of natural communities, and in the worst case, to native species extinctions. The introduction of gene modified organisms (GMOs) to agricultural, fisheries and forestry practices brings them into contact with natural populations. GMOs may appear as new invasive species if they are able to (1) invade into natural habitats or (2) hybridize with their wild relatives. The benefits of GMOs, such as increased yield or decreased use of insecticides or herbicides in cultivation, may thus be reduced due the potential risks they may cause. A careful ecological risk analysis therefore has to precede any responsible GMO introduction. In this thesis I study ecological invasion in relation to GMOs, and what kind of consequences invasion may have in natural populations. A set of theoretical models that combine life-history evolution, population dynamics, and population genetics were developed for the hazard identification part of ecological risks assessment of GMOs. In addition, the potential benefits of GMOs in management of an invasive pest were analyzed. In the first study I showed that a population that is fluctuating due to scramble-type density dependence (due to, e.g., nutrient competition in plants) may be invaded by a population that is relatively more limited by a resource (e.g., light in plants) that is a cause of contest-type density dependence. This result emphasises the higher risk of invasion in unstable environments. The next two studies focused on escape of a growth hormone (GH) transgenic fish into a natural population. The results showed that previous models may have given too pessimistic a view of the so called Trojan gene -effect, where the invading genotype is harmful for the population as a whole. The previously suggested population extinctions did not occur in my studies, since the changes in mating preferences caused by the GH-fish were be ameliorated by decreased level of competition. The GH-invaders may also have to exceed a threshold density before invasion can be successful. I also showed that the prevalence of mature parr (aka. sneaker) strategy among GH-fish may have clear effect on invasion outcome. The fourth study assessed the risks and developed methods against the invasion of the Colorado Potato Beetle (CPB, Leptinotarsa decemlineata). I showed that the eradication of CPB is most important for the prevention of their establishment, but the cultivation of transgenic Bt-potato could also be effective. In general, my results emphasise that invasion of transgenic species or genotypes to be possible under certain realistic conditions and resulting in competitive exclusion, population decline through outbreeding depression and genotypic displacement of native species. Ecological risk assessment should regard the decline and displacement of the wild genotype by an introduced one as a consequence that is as serious as the population extinction. It will also be crucial to take into account different kinds of behavioural differences among species when assessing the possible hazards that GMOs may cause if escaped. The benefits found of GMO crops effectiveness in pest management may also be too optimistic since CPB may evolve resistance to Bt-toxin. The models in this thesis could be further applied in case specific risk assessment of GMOs by supplementing them with detailed data of the species biology, the effect of the transgene introduced to the species, and also the characteristics of the populations or the environments in the risk of being invaded.

Genetic aspects of outcome in kidney transplantation: cytokine and thrombosis associated candidate genes and gene expression biomarkers

Kidney transplantation (Tx) is the treatment of choice for end stage renal disease. Immunosuppressive medications are given to prevent an immunological rejection of the transplant. However, immunosuppressive drugs increase e.g. the risk of infection, cancer or nephrotoxicity. A major genetic contributors to immunological acceptance of the graft are human leukocyte antigen (HLA) genes. Also other non-HLA gene polymorphisms may predict the future risk of complications before Tx, possibly enabling individualised immunotherapy. Graft function after Tx is monitored using non-specific clinical symptoms and laboratory markers. The definitive diagnosis of graft rejection however relies on a biopsy of the graft. In the acute rejection (AR) diagnostics there is a need for an alternative to biopsy that would be an easily repeatable and simple method for regular use. Frequent surveillance of acute or subclinical rejection (SCR) may improve long-term function. In this thesis, associations between cytokine and thrombosis associated candidate genes and the outcome of kidney Tx were studied. Cytotoxic and co-stimulatory T lymphocyte molecule gene expression biomarkers for the diagnosis of the AR and the SCR were also investigated. We found that polymorphisms in the cytokine genes tumor necrosis factor and interleukin 10 (IL10) of the recipients were associated with AR. In addition, certain IL10 gene polymorphisms of the donors were associated with the incidence of cytomegalovirus infection and occurrence of later infection in a subpopulation of recipients. Further, polymorphisms in genes related to the risk of thrombosis and those of certain cytokines were not associated with the occurrence of thrombosis, infarction, AR or graft survival. In the study of biomarkers for AR, whole blood samples were prospectively collected from adult kidney Tx patients. With real-time quantitative PCR (RT-QPCR) gene expression quantities of CD154 and ICOS differentiated the patients with AR from those without, but not from the patients with other causes of graft dysfunction. Biomarkers for SCR were studied in paediatric kidney Tx patients. We used RT-QPCR to quantify the gene expression of immunological candidate genes in a low-density array format. In addition, we used RT-QPCR to validate the results of the microarray analysis. No gene marker differentiated patients with SCR from those without SCR. This research demonstrates the lack of robust markers among polymorphisms or biomarkers in investigated genes that could be included in routine analysis in a clinical laboratory. In genetic studies, kidney Tx can be regarded as a complex trait, i.e. several environmental and genetic factors may determine its outcome. A number of currently unknown genetic factors probably influence the results of Tx.