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.

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.

The genetic associations of acute anterior uveitis and their overlap with the genetics of ankylosing spondylitis

Acute anterior uveitis (AAU) involves inflammation of the iris and ciliary body of the eye. It occurs both in isolation and as a complication of ankylosing spondylitis (AS). It is strongly associated with HLA-B*27, but previous studies have suggested that further genetic factors may confer additional risk. We sought to investigate this using the Illumina Exomechip microarray, to compare 1504 cases with AS and AAU, 1805 with AS but no AAU and 21 133 healthy controls. We also used a heterogeneity test to test the differences in effect size between AS with AAU and AS without AAU. In the analysis comparing AS+AAU+ cases versus controls, HLA-B*27 and HLA-A*02:01 were significantly associated with the presence of AAU (P<10−300 and P=6 × 10−8, respectively). Secondary independent association with PSORS1C3 (P=4.7 × 10−5) and TAP2 (P=1.1 × 10−5) were observed in the major histocompatibility complex. There was a new suggestive association with a low-frequency variant at zinc-finger protein 154 in the AS without AAU versus control analysis (zinc-finger protein 154 (ZNF154), P=2.2 × 10−6). Heterogeneity testing showed that rs30187 in ERAP1 has a larger effect on AAU compared with that in AS alone. These findings also suggest that variants in ERAP1 have a differential impact on the risk of AAU when compared with AS, and hence the genetic risk for AAU differs from AS.

trans gene regulation in adaptive evolution: A genetic algorithm model

This is a continuation of earlier studies on the evolution of infinite populations of haploid genotypes within a genetic algorithm framework. We had previously explored the evolutionary consequences of the existence of indeterminate—“plastic”—loci, where a plastic locus had a finite probability in each generation of functioning (being switched “on”) or not functioning (being switched “off”). The relative probabilities of the two outcomes were assigned on a stochastic basis. The present paper examines what happens when the transition probabilities are biased by the presence of regulatory genes. We find that under certain conditions regulatory genes can improve the adaptation of the population and speed up the rate of evolution (on occasion at the cost of lowering the degree of adaptation). Also, the existence of regulatory loci potentiates selection in favour of plasticity. There is a synergistic effect of regulatory genes on plastic alleles: the frequency of such alleles increases when regulatory loci are present. Thus, phenotypic selection alone can be a potentiating factor in a favour of better adaptation.

Genetic studies of the PRP17 gene of Saccharomyces cerevisiae: a domain essential for function maps to a nonconserved region of the protein

The PRP17 gene product is required for the second step of pre-mRNA splicing reactions. The C-terminal half of this protein bears four repeat units with homology to the beta transducin repeat. Missense mutations in three temperature-sensitive prp17 mutants map to a region in the N-terminal half of the protein. We have generated, in vitro, 11 missense alleles at the beta transducin repeat units and find that only one affects function in vivo. A phenotypically silent missense allele at the fourth repeat unit enhances the slow-growing phenotype conferred by an allele at the third repeat, suggesting an interaction between these domains. Although many missense mutations in highly conserved amino acids lack phenotypic effects, deletion analysis suggests an essential role for these units. Only mutations in the N-terminal nonconserved domain of PRP17 are synthetically lethal in combination with mutations in PRP16 and PRP18, two other gene products required for the second splicing reaction. A mutually allele-specific interaction between Prp17 and snr7, with mutations in U5 snRNA, was observed. We therefore suggest that the functional region of Prp17p that interacts with Prp18p, Prp16p, and U5 snRNA is the N terminal region of the protein.

Flight metabolic rate in the Glanville fritillary butterfly

Dispersal is a highly important life history trait. In fragmented landscapes the long-term persistence of populations depends on dispersal. Evolution of dispersal is affected by costs and benefits and these may differ between different landscapes. This results in differences in the strength and direction of natural selection on dispersal in fragmented landscapes. Dispersal has been shown to be a nonrandom process that is associated with traits such as flight ability in insects. This thesis examines genetic and physiological traits affecting dispersal in the Glanville fritillary butterfly (Melitaea cinxia). Flight metabolic rate is a repeatable trait representing flight ability. Unlike in many vertebrates, resting metabolic rate cannot be used as a surrogate of maximum metabolic rate as no strong correlation between the two was found in the Glanville fritillary. Resting and flight metabolic rate are affected by environmental variables, most notably temperature. However, only flight metabolic rate has a strong genetic component. Molecular variation in the much-studied candidate locus phosphoglucose isomerase (Pgi), which encodes the glycolytic enzyme PGI, has an effect on carbohydrate metabolism in flight. This effect is temperature dependent: in low to moderate temperatures individuals with the heterozygous genotype at the single nucleotide polymorphism (SNP) AA111 have higher flight metabolic rate than the common homozygous genotype. At high temperatures the situation is reversed. This finding suggests that variation in enzyme properties is indeed translated to organismal performance. High-resolution data on individual female Glanville fritillaries moving freely in the field were recorded using harmonic radar. There was a strong positive correlation between flight metabolic rate and dispersal rate. Flight metabolic rate explained one third of the observed variation in the one-hour movement distance. A fine-scaled analysis of mobility showed that mobility peaked at intermediate ambient temperatures but the two common Pgi genotypes differed in their reaction norms to temperature. As with flight metabolic rate, heterozygotes at SNP AA111 were the most active genotype in low to moderate temperatures. The results show that molecular variation is associated with variation in dispersal rate through the link of flight physiology under the influence of environmental conditions. The evolutionary pressures for dispersal differ between males and females. The effect of flight metabolic rate on dispersal was examined in both sexes in field and laboratory conditions. The relationship between flight metabolic rate and dispersal rate in the field and flight duration in the laboratory were found to differ between the two sexes. In females the relationship was positive, but in males the longest distances and flight durations were recorded for individuals with low flight metabolic rate. These findings may reflect male investment in mate locating. Instead of dispersing, males with high flight metabolic rate may establish territories and follow a perching strategy when locating females and hence move less on the landscape level. Males with low metabolic rate may be forced to disperse due to low competitive success or may show adaptations to an alternative strategy: patrolling. In the light of life history trade-offs and the rate of living theory having high metabolic rate may carry a cost in the form of shortened lifespan. Experiments relating flight metabolic rate to longevity showed a clear correlation in the opposite direction: high flight metabolic rate was associated with long lifespan. This suggests that individuals with high metabolic rate do not pay an extra physiological cost for their high flight capacity, rather there are positive correlations between different measures of fitness. These results highlight the importance of condition.

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.

Applications of gene sequence polymorphisms in evolutionary genetic studies of Atlantic salmon (Salmo salar) and other teleost fish species

Evolutionary genetics incorporates traditional population genetics and studies of the origins of genetic variation by mutation and recombination, and the molecular evolution of genomes. Among the primary forces that have potential to affect the genetic variation within and among populations, including those that may lead to adaptation and speciation, are genetic drift, gene flow, mutations and natural selection. The main challenges in knowing the genetic basis of evolutionary changes is to distinguish the adaptive selection forces that cause existent DNA sequence variants and also to identify the nucleotide differences responsible for the observed phenotypic variation. To understand the effects of various forces, interpretation of gene sequence variation has been the principal basis of many evolutionary genetic studies. The main aim of this thesis was to assess different forms of teleost gene sequence polymorphisms in evolutionary genetic studies of Atlantic salmon (Salmo salar) and other species. Firstly, the level of Darwinian adaptive evolution affected coding regions of the growth hormone (GH) gene during the teleost evolution was investigated based on the sequence data existing in public databases. Secondly, a target gene approach was used to identify within population variation in the growth hormone 1 (GH1) gene in salmon. Then, a new strategy for single nucleotide polymorphisms (SNPs) discovery in salmonid fishes was introduced, and, finally, the usefulness of a limited number of SNP markers as molecular tools in several applications of population genetics in Atlantic salmon was assessed. This thesis showed that the gene sequences in databases can be utilized to perform comparative studies of molecular evolution, and some putative evidence of the existence of Darwinian selection during the teleost GH evolution was presented. In addition, existent sequence data was exploited to investigate GH1 gene variation within Atlantic salmon populations throughout its range. Purifying selection is suggested to be the predominant evolutionary force controlling the genetic variation of this gene in salmon, and some support for gene flow between continents was also observed. The novel approach to SNP discovery in species with duplicated genome fragments introduced here proved to be an effective method, and this may have several applications in evolutionary genetics with different species - e.g. when developing gene-targeted markers to investigate quantitative genetic variation. The thesis also demonstrated that only a few SNPs performed highly similar signals in some of the population genetic analyses when compared with the microsatellite markers. This may have useful applications when estimating genetic diversity in genes having a potential role in ecological and conservation issues, or when using hard biological samples in genetic studies as SNPs can be applied with relatively highly degraded DNA.