Conceptual and statistical modelling of environmental effects in population dynamics

Population dynamics are generally viewed as the result of intrinsic (purely density dependent) and extrinsic (environmental) processes. Both components, and potential interactions between those two, have to be modelled in order to understand and predict dynamics of natural populations; a topic that is of great importance in population management and conservation. This thesis focuses on modelling environmental effects in population dynamics and how effects of potentially relevant environmental variables can be statistically identified and quantified from time series data. Chapter I presents some useful models of multiplicative environmental effects for unstructured density dependent populations. The presented models can be written as standard multiple regression models that are easy to fit to data. Chapters II IV constitute empirical studies that statistically model environmental effects on population dynamics of several migratory bird species with different life history characteristics and migration strategies. In Chapter II, spruce cone crops are found to have a strong positive effect on the population growth of the great spotted woodpecker (Dendrocopos major), while cone crops of pine another important food resource for the species do not effectively explain population growth. The study compares rate- and ratio-dependent effects of cone availability, using state-space models that distinguish between process and observation error in the time series data. Chapter III shows how drought, in combination with settling behaviour during migration, produces asymmetric spatially synchronous patterns of population dynamics in North American ducks (genus Anas). Chapter IV investigates the dynamics of a Finnish population of skylark (Alauda arvensis), and point out effects of rainfall and habitat quality on population growth. Because the skylark time series and some of the environmental variables included show strong positive autocorrelation, the statistical significances are calculated using a Monte Carlo method, where random autocorrelated time series are generated. Chapter V is a simulation-based study, showing that ignoring observation error in analyses of population time series data can bias the estimated effects and measures of uncertainty, if the environmental variables are autocorrelated. It is concluded that the use of state-space models is an effective way to reach more accurate results. In summary, there are several biological assumptions and methodological issues that can affect the inferential outcome when estimating environmental effects from time series data, and that therefore need special attention. The functional form of the environmental effects and potential interactions between environment and population density are important to deal with. Other issues that should be considered are assumptions about density dependent regulation, modelling potential observation error, and when needed, accounting for spatial and/or temporal autocorrelation.

The role of cyclase-associated protein (CAP) in actin dynamics during cell motility and morphogenesis

The highly dynamic remodeling of the actin cytoskeleton is responsible for most motile and morphogenetic processes in all eukaryotic cells. In order to generate appropriate spatial and temporal movements, the actin dynamics must be under tight control of an array of actin binding proteins (ABPs). Many proteins have been shown to play a specific role in actin filament growth or disassembly of older filaments. Very little is known about the proteins affecting recycling i.e. the step where newly depolymerized actin monomers are funneled into new rounds of filament assembly. A central protein family involved in the regulation of actin turnover is cyclase-associated proteins (CAP, called Srv2 in budding yeast). This 50-60 kDa protein was first identified from yeast as a suppressor of an activated RAS-allele and a factor associated with adenylyl cyclase. The CAP proteins harbor N-terminal coiled-coil (cc) domain, originally identified as a site for adenylyl cyclase binding. In the N-terminal half is also a 14-3-3 like domain, which is followed by central proline-rich domains and the WH2 domain. In the C-terminal end locates the highly conserved ADP-G-actin binding domain. In this study, we identified two previously suggested but poorly characterized interaction partners for Srv2/CAP: profilin and ADF/cofilin. Profilins are small proteins (12-16 kDa) that bind ATP-actin monomers and promote the nucleotide exchange of actin. The profilin-ATP-actin complex can be directly targeted to the growth of the filament barbed ends capped by Ena/VASP or formins. ADF/cofilins are also small (13-19 kDa) and highly conserved actin binding proteins. They depolymerize ADP-actin monomers from filament pointed ends and remain bound to ADP-actin strongly inhibiting nucleotide exchange. We revealed that the ADP-actin-cofilin complex is able to directly interact with the 14-3-3 like domain at the N-terminal region of Srv2/CAP. The C-terminal high affinity ADP-actin binding site of Srv2/CAP competes with cofilin for an actin monomer. Cofilin can thus be released from Srv2/CAP for the subsequent round of depolymerization. We also revealed that profilin interacts with the first proline-rich region of Srv2/CAP and that the binding occurs simultaneously with ADP-actin binding to C-terminal domain of Srv2/CAP. Both profilin and Srv2/CAP can promote nucleotide exchange of actin monomer. Because profilin has much higher affinity to ATP-actin than Srv2/CAP, the ATP-actin-profilin complex is released for filament polymerization. While a disruption of cofilin binding in yeast Srv2/CAP produces a severe phenotype comparable to Srv2/CAP deletion, an impairment of profilin binding from Srv2/CAP results in much milder phenotype. This suggests that the interaction with cofilin is essential for the function of Srv2/CAP, whereas profilin can also promote its function without direct interaction with Srv2/CAP. We also show that two CAP isoforms with specific expression patterns are present in mice. CAP1 is the major isoform in most tissues, while CAP2 is predominantly expressed in muscles. Deletion of CAP1 from non-muscle cells results in severe actin phenotype accompanied with mislocalization of cofilin to cytoplasmic aggregates. Together these studies suggest that Srv2/CAP recycles actin monomers from cofilin to profilin and thus it plays a central role in actin dynamics in both yeast and mammalian cells.

Ecological communities in variable environments : dynamics and diversity under coloured environmental stochasticity

While environmental variation is an ubiquitous phenomenon in the natural world which has for long been appreciated by the scientific community recent changes in global climatic conditions have begun to raise consciousness about the economical, political and sociological ramifications of global climate change. Climate warming has already resulted in documented changes in ecosystem functioning, with direct repercussions on ecosystem services. While predicting the influence of ecosystem changes on vital ecosystem services can be extremely difficult, knowledge of the organisation of ecological interactions within natural communities can help us better understand climate driven changes in ecosystems. The role of environmental variation as an agent mediating population extinctions is likely to become increasingly important in the future. In previous studies population extinction risk in stochastic environmental conditions has been tied to an interaction between population density dependence and the temporal autocorrelation of environmental fluctuations. When populations interact with each other, forming ecological communities, the response of such species assemblages to environmental stochasticity can depend, e.g., on trophic structure in the food web and the similarity in species-specific responses to environmental conditions. The results presented in this thesis indicate that variation in the correlation structure between species-specific environmental responses (environmental correlation) can have important qualitative and quantitative effects on community persistence and biomass stability in autocorrelated (coloured) environments. In addition, reddened environmental stochasticity and ecological drift processes (such as demographic stochasticity and dispersal limitation) have important implications for patterns in species relative abundances and community dynamics over time and space. Our understanding of patterns in biodiversity at local and global scale can be enhanced by considering the relevance of different drift processes for community organisation and dynamics. Although the results laid out in this thesis are based on mathematical simulation models, they can be valuable in planning effective empirical studies as well as in interpreting existing empirical results. Most of the metrics considered here are directly applicable to empirical data.

Mycorrhizal fungi and nitrogen dynamics in drained peatland

Fungi have a fundamental role in carbon and nutrient transformations in the acids soils of boreal regions, such as peatlands, where high amounts of carbon (C) and nutrients are stored in peat, the pH is relatively low and the nutrient uptake of trees is highly dependent on mycorrhizae. In this thesis, the aim was to examine nitrogen (N) transformations and the availability of dissolved N compounds in forestry-drained peatlands, to compare the fungal community biomass and structure at various peat N levels, to investigate the growth of ectomycorrhizal fungi with variable P and K availability and to assess how the ectomycorrhizal fungi (ECM) affect N transformations. Both field and laboratory experiments were carried out. The peat N concentration did not affect the soil fungal community structure within a site. Phosphorus (P) and potassium (K) deficiency of the trees as well as the degree of decomposition and dissolved organic nitrogen (DON) concentration of the peat were shown to affect the fungal community structure and biomass of ECMs, highlighting the complexity of the below ground system on drained peatlands. The biomass of extrametrical mycorrhizal mycelia (EMM) was enhanced by P and/or K deficiency of the trees, and ECM biomass in the roots was increased by P deficiency. Thus, PK deficiency in drained peatlands may increase the allocation of C by the tree to ECMs. It was also observed that fungi can alter N mineralization processes in the rhizosphere but variously depending on fungal species and fertility level of peat. Gross N mineralization did not vary but the net N mineralization rate significantly increased along the N gradient in both field and laboratory experiments. Gross N immobilization also significantly increased when the peat N concentration increased. Nitrification was hardly detectable in either field or laboratory experiments. During the growing season, dissolved inorganic N (DIN) fluctuated much more than the relatively stable DON. Special methodological challenges associated with sampling and analysis in microbial studies on peatlands are discussed.

Actin Dynamics in Muscle Cells

In every cell, actin is a key component involved in migration, cytokinesis, endocytosis and generation of contraction. In non-muscle cells, actin filaments are very dynamic and regulated by an array of proteins that interact with actin filaments and/or monomeric actin. Interestingly, in non-muscle cells the barbed ends of the filaments are the predominant assembly place, whereas in muscle cells actin dynamics was reported to predominate at the pointed ends of thin filaments. The actin-based thin filament pointed (slow growing) ends extend towards the middle of the sarcomere's M-line where they interact with the thick filaments to generate contraction. The actin filaments in muscle cells are organized into a nearly crystalline array and are believed to be significantly less dynamic than the ones in other cell types. However, the exact mechanisms of the sarcomere assembly and turnover are largely unknown. Interestingly, although sarcomeric actin structures are believed to be relatively non-dynamic, many proteins promoting actin dynamics are expressed also in muscle cells (e.g ADF/cofilin, cyclase-associated protein and twinfilin). Thus, it is possible that the muscle-specific isoforms of these proteins promote actin dynamics differently from their non-muscle counterparts, or that actin filaments in muscle cells are more dynamic than previously thought. To study protein dynamics in live muscle cells, I used primary cell cultures of rat cardiomyocytes. My studies revealed that a subset of actin filaments in cardiomyocyte sarcomeres displays rapid turnover. Importantly, I discovered that the turnover of actin filaments depends on contractility of the cardiomyocytes and that the contractility-induced actin dynamics plays an important role in sarcomere maturation. Together with previous studies those findings suggest that sarcomeres undergo two types of actin dynamics: (1) contractility-dependent turnover of whole filaments and (2) regulatory pointed end monomer exchange to maintain correct thin filament length. Studies involving an actin polymerization inhibitor suggest that the dynamic actin filament pool identified here is composed of filaments that do not contribute to contractility. Additionally, I provided evidence that ADF/cofilins, together with myosin-induced contractility, are required to disassemble non-productive filaments in developing cardiomyocytes. In addition, during these studies we learned that isoforms of actin monomer binding protein twinfilin, Twf-1 and Twf-2a localise to myofibrils in cardiomyocytes and may thus contribute to actin dynamics in myofibrils. Finally, in collaboration with Roberto Dominguez s laboratory we characterized a new actin nucleator in muscle cells - leiomodin (Lmod). Lmod localises towards actin filament pointed ends and its depletion by siRNA leads to severe sarcomere abnormalities in cardiomyocytes. The actin filament nucleation activity of Lmod is enhanced by interactions with tropomyosin. We also revealed that Lmod expression correlates with the maturation of myofibrils, and that it associates with sarcomeres only at relatively late stages of myofibrillogenesis. Thus, Lmod is unlikely to play an important role in myofibril formation, but rather might be involved in the second step of the filament arrangement and/or maintenance through its ability to promote tropomyosin-induced actin filament nucleation occurring at the filament pointed ends. The results of these studies provide valuable new information about the molecular mechanisms underlying muscle sarcomere assembly and turnover. These data offer important clues to understanding certain physiological and pathological behaviours of muscle cells. Better understanding of the processes occurring in muscles might help to find strategies for determining, diagnosis, prognosis and therapy in heart and skeletal muscles diseases.

Structure, function and intracellular dynamics of alphavirus replication complexes

Intracellular membrane alterations are hallmarks of positive-sense RNA (+RNA) virus replication. Strong evidence indicates that within these exotic compartments, viral replicase proteins engage in RNA genome replication and transcription. To date, fundamental questions such as the origin of altered membranes, mechanisms of membrane deformation and topological distribution and function of viral components, are still waiting for comprehensive answers. This study addressed some of the above mentioned questions for the membrane alterations induced during Semliki Forest virus (SFV) infection of mammalian cells. With the aid of electron and fluorescence microscopy coupled with radioactive labelling and immuno-cytochemistry techniques, our group and others showed that few hours after infection the four non structural proteins (nsP1-4) and newly synthesized RNAs of SFV colocalized in close proximity of small membrane invaginations, designated as spherules . These 50-70 nm structures were mainly detected in the perinuclear area, at the limiting membrane of modified endosomes and lysosomes, named CPV-I (cytopathic vacuoles type I). More rarely, spherules were also found at the plasma membrane (PM). In the first part of this study I present the first three-dimensional reconstruction of the CPV-I and the spherules, obtained by electron tomography after chemical or cryo-fixation. Different approaches for imaging these macromolecular assemblies to obtain better structure preservation and higher resolution are presented as unpublished data. This study provides insights into spherule organization and distribution of viral components. The results of this and other experiments presented in this thesis will challenge currently accepted models for virus replication complex formation and function. In a revisitation of our previous models, the second part of this work provides the first complete description of the biogenesis of the CPV-I. The results demonstrate that these virus-induced vacuoles, where hundreds of spherules accumulate at late stages during infection, represent the final phase of a journey initiated at the PM, which apparently serves as a platform for spherule formation. From the PM spherules were internalized by an endocytic event that required the activity of the class I PI3K, caveolin-1, cellular cholesterol and functional actin-myosin network. The resulting neutral endocytic carrier vesicle delivered the spherules to the membrane of pre-existing acidic endosomes via multiple fusion events. Microtubule based transport supported the vectorial transfer of these intermediates to the pericentriolar area where further fusions generated the CPV-I. A signal for spherule internalization was identified in one of the replicase proteins, nsP3. Infections of cells with viruses harbouring a deletion in a highly phosphorylated region of nsP3 did not result in the formation of CPV-Is. Instead, thousands of spherules remained at the PM throughout the infection cycle. Finally, the role of the replicase protein nsP2 during viral RNA replication and transcription was investigated. Three enzymatic activities, protease, NTPase and RNA-triphosphatase were studied with the aid of temperature sensitive mutants in vitro and, when possible, in vivo. The results highlighted the interplay of the different nsP2 functions during different steps of RNA replication and sub-genomic promoter regulation, and suggest that the protein could have different activities when participating in the replication complex or as a free enzyme.

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.

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.

Mechanism of RNA Translocation by a Viral Packaging Motor

Molecular motors are proteins that convert chemical energy into mechanical work. The viral packaging ATPase P4 is a hexameric molecular motor that translocates RNA into preformed viral capsids. P4 belongs to the ubiquitous class of hexameric helicases. Although its structure is known, the mechanism of RNA translocation remains elusive. Here we present a detailed kinetic study of nucleotide binding, hydrolysis, and product release by P4. We propose a stochastic-sequential cooperative model to describe the coordination of ATP hydrolysis within the hexamer. In this model the apparent cooperativity is a result of hydrolysis stimulation by ATP and RNA binding to neighboring subunits rather than cooperative nucleotide binding. Simultaneous interaction of neighboring subunits with RNA makes the otherwise random hydrolysis sequential and processive. Further, we use hydrogen/deuterium exchange detected by high resolution mass spectrometry to visualize P4 conformational dynamics during the catalytic cycle. Concerted changes of exchange kinetics reveal a cooperative unit that dynamically links ATP binding sites and the central RNA binding channel. The cooperative unit is compatible with the structure-based model in which translocation is effected by conformational changes of a limited protein region. Deuterium labeling also discloses the transition state associated with RNA loading which proceeds via opening of the hexameric ring. Hydrogen/deuterium exchange is further used to delineate the interactions of the P4 hexamer with the viral procapsid. P4 associates with the procapsid via its C-terminal face. The interactions stabilize subunit interfaces within the hexamer. The conformation of the virus-bound hexamer is more stable than the hexamer in solution, which is prone to spontaneous ring openings. We propose that the stabilization within the viral capsid increases the packaging processivity and confers selectivity during RNA loading. Finally, we use single molecule techniques to characterize P4 translocation along RNA. While the P4 hexamer encloses RNA topologically within the central channel, it diffuses randomly along the RNA. In the presence of ATP, unidirectional net movement is discernible in addition to the stochastic motion. The diffusion is hindered by activation energy barriers that depend on the nucleotide binding state. The results suggest that P4 employs an electrostatic clutch instead of cycling through stable, discrete, RNA binding states during translocation. Conformational changes coupled to ATP hydrolysis modify the electrostatic potential inside the central channel, which in turn biases RNA motion in one direction. Implications of the P4 model for other hexameric molecular motors are discussed.