Eurasian large mammal dynamics in response to changing environments during the Late Neogene

Nisäkkäiden levinneisyyteen, niiden morfologisiin ja ekologisiin piirteisiin vaikuttavat ympäristön sekä lyhyet että pitkäkestoiset muutokset, etenkin ilmaston ja kasvillisuuden vaihtelut. Työssä tutkittiin nisäkkäiden sopeutumista ilmastonmuutoksiin Euraasiassa viimeisen 24 miljoonan vuoden aikana. Tutkimuksessa keskityttiin varsinkin viimeiseen kahteen miljoonaan vuoteen, jonka aikana ilmasto muuttui voimakkaasti ja ihmisen toiminta alkoi tulla merkittäväksi. Tämän takia on usein vaikea erottaa, kummasta em. seikasta jonkin nisäkäslajin sukupuutto tai häviäminen alueelta johtui. Aineistona käytettiin laajaa venäjänkielistä kirjallisuutta, josta löytyvät tiedot ovat kääntämättöminä jääneet aiemmin länsimaisen tutkimuksen ulkopuolelle. Työssä käytettiin myös NOW-tietokantaa, jossa on fossiilisten nisäkkäiden löytöpaikat sekä niiden iät.

Computational Methods for Detecting Large-Scale Chromosome Rearrangements in SNP Data

Large-scale chromosome rearrangements such as copy number variants (CNVs) and inversions encompass a considerable proportion of the genetic variation between human individuals. In a number of cases, they have been closely linked with various inheritable diseases. Single-nucleotide polymorphisms (SNPs) are another large part of the genetic variance between individuals. They are also typically abundant and their measuring is straightforward and cheap. This thesis presents computational means of using SNPs to detect the presence of inversions and deletions, a particular variety of CNVs. Technically, the inversion-detection algorithm detects the suppressed recombination rate between inverted and non-inverted haplotype populations whereas the deletion-detection algorithm uses the EM-algorithm to estimate the haplotype frequencies of a window with and without a deletion haplotype. As a contribution to population biology, a coalescent simulator for simulating inversion polymorphisms has been developed. Coalescent simulation is a backward-in-time method of modelling population ancestry. Technically, the simulator also models multiple crossovers by using the Counting model as the chiasma interference model. Finally, this thesis includes an experimental section. The aforementioned methods were tested on synthetic data to evaluate their power and specificity. They were also applied to the HapMap Phase II and Phase III data sets, yielding a number of candidates for previously unknown inversions, deletions and also correctly detecting known such rearrangements.

The Gerbera cDNA Microarray: A Tool for Large-Scale Identification of Genes Involved in Flower Development

During the past ten years, large-scale transcript analysis using microarrays has become a powerful tool to identify and predict functions for new genes. It allows simultaneous monitoring of the expression of thousands of genes and has become a routinely used tool in laboratories worldwide. Microarray analysis will, together with other functional genomics tools, take us closer to understanding the functions of all genes in genomes of living organisms. Flower development is a genetically regulated process which has mostly been studied in the traditional model species Arabidopsis thaliana, Antirrhinum majus and Petunia hybrida. The molecular mechanisms behind flower development in them are partly applicable in other plant systems. However, not all biological phenomena can be approached with just a few model systems. In order to understand and apply the knowledge to ecologically and economically important plants, other species also need to be studied. Sequencing of 17 000 ESTs from nine different cDNA libraries of the ornamental plant Gerbera hybrida made it possible to construct a cDNA microarray with 9000 probes. The probes of the microarray represent all different ESTs in the database. From the gerbera ESTs 20% were unique to gerbera while 373 were specific to the Asteraceae family of flowering plants. Gerbera has composite inflorescences with three different types of flowers that vary from each other morphologically. The marginal ray flowers are large, often pigmented and female, while the central disc flowers are smaller and more radially symmetrical perfect flowers. Intermediate trans flowers are similar to ray flowers but smaller in size. This feature together with the molecular tools applied to gerbera, make gerbera a unique system in comparison to the common model plants with only a single kind of flowers in their inflorescence. In the first part of this thesis, conditions for gerbera microarray analysis were optimised including experimental design, sample preparation and hybridization, as well as data analysis and verification. Moreover, in the first study, the flower and flower organ-specific genes were identified. After the reliability and reproducibility of the method were confirmed, the microarrays were utilized to investigate transcriptional differences between ray and disc flowers. This study revealed novel information about the morphological development as well as the transcriptional regulation of early stages of development in various flower types of gerbera. The most interesting finding was differential expression of MADS-box genes, suggesting the existence of flower type-specific regulatory complexes in the specification of different types of flowers. The gerbera microarray was further used to profile changes in expression during petal development. Gerbera ray flower petals are large, which makes them an ideal model to study organogenesis. Six different stages were compared and specifically analysed. Expression profiles of genes related to cell structure and growth implied that during stage two, cells divide, a process which is marked by expression of histones, cyclins and tubulins. Stage 4 was found to be a transition stage between cell division and expansion and by stage 6 cells had stopped division and instead underwent expansion. Interestingly, at the last analysed stage, stage 9, when cells did not grow any more, the highest number of upregulated genes was detected. The gerbera microarray is a fully-functioning tool for large-scale studies of flower development and correlation with real-time RT-PCR results show that it is also highly sensitive and reliable. Gene expression data presented here will be a source for gene expression mining or marker gene discovery in the future studies that will be performed in the Gerbera Laboratory. The publicly available data will also serve the plant research community world-wide.

Effects of large-scale human land use on Capercaillie (Tetrao urogallus L.) populations in Finland

The Capercaillie (Tetrao urogallus L.) is often used as a focal species for landscape ecological studies: the minimum size for its lekking area is 300 ha, and the annual home range for an individual may cover 30 80 km2. In Finland, Capercaillie populations have decreased by approximately 40 85%, with the declines likely to have started in the 1940s. Although the declines have partly stabilized from the 1990s onwards, it is obvious that the negative population trend was at least partly caused by changes in human land use. The aim of this thesis was to study the connections between human land use and Capercaillie populations in Finland, using several spatial and temporal scales. First, the effect of forest age structure on Capercaillie population trends was studied in 18 forestry board districts in Finland, during 1965 1988. Second, the abundances of Capercaillie and Moose (Alces alces L.) were compared in terms of several land-use variables on a scale of 50 × 50 km grids and in five regions in Finland. Third, the effects of forest cover and fine-grain forest fragmentation on Capercaillie lekking area persistence were studied in three study locations in Finland, on 1000 and 3000 m spatial scales surrounding the leks. The analyses considering lekking areas were performed with two definitions for forest: > 60 and > 152 m3ha 1 of timber volume. The results show that patterns and processes at large spatial scales strongly influence Capercaillie in Finland. In particular, in southwestern and eastern Finland, high forest cover and low human impact were found to be beneficial for this species. Forest cover (> 60 m3ha 1 of timber) surrounding the lekking sites positively affected lekking area persistence only at the larger landscape scale (3000 m radius). The effects of older forest classes were hard to assess due to scarcity of older forests in several study areas. Young and middle-aged forest classes were common in the vicinity of areas with high Capercaillie abundances especially in northern Finland. The increase in the amount of younger forest classes did not provide a good explanation for Capercaillie population decline in 1965 1988. In addition, there was no significant connection between mature forests (> 152 m3ha 1 of timber) and lekking area persistence in Finland. It seems that in present-day Finnish landscapes, area covered with old forest is either too scarce to efficiently explain the abundance of Capercaillie and the persistence of the lekking areas, or the effect of forest age is only important when considering smaller spatial scales than the ones studied in this thesis. In conclusion, larger spatial scales should be considered for assessing the future Capercaillie management. According to the proposed multi-level planning, the first priority should be to secure the large, regional-scale forest cover, and the second priority should be to maintain fine-grained, heterogeneous structure within the separate forest patches. A management unit covering hundreds of hectares, or even tens or hundreds of square kilometers, should be covered, which requires regional-level land-use planning and co-operation between forest owners.

Prognostic molecular factors and algorithms in diffuse large B-cell lymphoma

Diffuse large B-cell lymphoma (DLBCL) is the most common of the non-Hodgkin lymphomas. As DLBCL is characterized by heterogeneous clinical and biological features, its prognosis varies. To date, the International Prognostic Index has been the strongest predictor of outcome for DLBCL patients. However, no biological characters of the disease are taken into account. Gene expression profiling studies have identified two major cell-of-origin phenotypes in DLBCL with different prognoses, the favourable germinal centre B-cell-like (GCB) and the unfavourable activated B-cell-like (ABC) phenotypes. However, results of the prognostic impact of the immunohistochemically defined GCB and non-GCB distinction are controversial. Furthermore, since the addition of the CD20 antibody rituximab to chemotherapy has been established as the standard treatment of DLBCL, all molecular markers need to be evaluated in the post-rituximab era. In this study, we aimed to evaluate the predictive value of immunohistochemically defined cell-of-origin classification in DLBCL patients. The GCB and non-GCB phenotypes were defined according to the Hans algorithm (CD10, BCL6 and MUM1/IRF4) among 90 immunochemotherapy- and 104 chemotherapy-treated DLBCL patients. In the chemotherapy group, we observed a significant difference in survival between GCB and non-GCB patients, with a good and a poor prognosis, respectively. However, in the rituximab group, no prognostic value of the GCB phenotype was observed. Likewise, among 29 high-risk de novo DLBCL patients receiving high-dose chemotherapy and autologous stem cell transplantation, the survival of non-GCB patients was improved, but no difference in outcome was seen between GCB and non-GCB subgroups. Since the results suggested that the Hans algorithm was not applicable in immunochemotherapy-treated DLBCL patients, we aimed to further focus on algorithms based on ABC markers. We examined the modified activated B-cell-like algorithm based (MUM1/IRF4 and FOXP1), as well as a previously reported Muris algorithm (BCL2, CD10 and MUM1/IRF4) among 88 DLBCL patients uniformly treated with immunochemotherapy. Both algorithms distinguished the unfavourable ABC-like subgroup with a significantly inferior failure-free survival relative to the GCB-like DLBCL patients. Similarly, the results of the individual predictive molecular markers transcription factor FOXP1 and anti-apoptotic protein BCL2 have been inconsistent and should be assessed in immunochemotherapy-treated DLBCL patients. The markers were evaluated in a cohort of 117 patients treated with rituximab and chemotherapy. FOXP1 expression could not distinguish between patients, with favourable and those with poor outcomes. In contrast, BCL2-negative DLBCL patients had significantly superior survival relative to BCL2-positive patients. Our results indicate that the immunohistochemically defined cell-of-origin classification in DLBCL has a prognostic impact in the immunochemotherapy era, when the identifying algorithms are based on ABC-associated markers. We also propose that BCL2 negativity is predictive of a favourable outcome. Further investigational efforts are, however, warranted to identify the molecular features of DLBCL that could enable individualized cancer therapy in routine patient care.

Iodine Brachytherapy for Large Uveal Melanomas

Uveal melanoma is the most common primary intraocular malignancy in adults. Vision in the affected eye is threatened by both the tumor and side-effects from the treatments currently available. Poor prognosis for saving vision increases with tumor size and, consequently, enucleation has been the treatment of choice for large uveal melanomas in most centers. However, increasing evidence suggests that no survival benefit is gained (nor lost) by enucleation as compared to eye-conserving methods. The Helsinki University Eye Hospital has since 1990 offered episcleral iodine-125 plaque brachytherapy (IBT) for all patients unwilling to undergo enucleation for a large uveal melanoma. The primary aim of this study was to assess survival, local tumor recurrence and preservation of the eye and vision after IBT in a population-based series of 97 patients with uveal melanomas classified as large by the Collaborative Ocular Melanoma Study (COMS) criteria. Further aims included reporting the incidence of side-effects and assessing the role of intraocular dose distribution and clinical risk factors in their development. Finally, means to improve the current treatment were investigated by using computer models to compare existing plaques with collimating ones and by comparing the outcome of a subgroup of 54 IBT patients with very thick tumors with 33 patients with similarly-sized tumors managed with transscleral local resection (TSR) in Liverpool, United Kingdom. Kaplan-Meier estimates of all-cause and melanoma-specific survival at 5 years after IBT were 62% and 65%, respectively, and visually comparable with the survival experience of patients reported after enucleation by the COMS. Local recurrence developed in 6% of eyes and 84% of eyes were conserved at 5 years. Visual prognosis was guarded with 11% avoiding loss of 20/70 vision and 26% avoiding loss of 20/400 vision in the tumor eye at 2 years. Large tumor height and short distance from the posterior pole were independently associated with loss of vision. Using cumulative incidence analysis to account for competing risks, such as enucleation and metastatic death, the 5-year incidence of cataract after IBT was 79%, glaucoma 60%, optic neuropathy 46%, maculopathy 52%, persistent or recurring retinal detachment (RD) 25%, and vitreous hemorrhage 36%. In multivariate competing risks regression models, increasing tumor height was associated with cataract, iris neovascularization and RD. Maculopathy and optic neuropathy were associated with distance from the tumor to the respective structure. Median doses to the tumor apex, macula and optic disc were 81 Gy (range, 40-158), 79 Gy (range, 12-632), and 83 Gy (range, 10-377), respectively. Dose to the optic disc was independently associated with optic neuropathy, and both dose to the optic disc and dose to the macula predicted vision loss after IBT. Simulated treatment using collimating plaques resulted in clinically meaningful reduction in both optic disc (median reduction, 30 Gy) and macular (median reduction, 36 Gy) doses as compared to the actual treatment with standard plaques. In the subgroup of patients with uveal melanomas classified as large because of tumor height, cumulative incidence analysis revealed that while long-term preservation of 20/70 vision was rare after both IBT and TSR, preservation of 20/400 vision was better after TSR (32% vs. 5% at 5 years). In multivariate logistic regression models, TSR was independently associated with better preservation of 20/400 vision (OR 0.03 at 2 years, P=0.005) No cases of secondary glaucoma were observed after TSR and optic neuropathy was rare. However, local tumor recurrence was more common after TSR than it was after IBT (Cumulative incidence 41% vs. 7% at 5 years, respectively). In terms of survival, IBT seems to be a safe alternative to enucleation in managing large uveal melanomas. Local tumor control is no worse than with medium-sized tumors and the chances of avoiding secondary enucleation are good. Unfortunately, side-effects from radiotherapy are frequent, especially in thick tumors, and long-term prognosis of saving vision is consequently guarded. Some complications can be limited by using collimating plaques and by managing uveal melanomas that are large because of tumor height with TSR instead of IBT. However, the patient must be willing to accept a substantial risk of local tumor recurrence after TSR and it is best suited for cases in which the preservation of vision in the tumor eye is critical.

Asteroid identification using statistical orbital inversion methods

An efficient and statistically robust solution for the identification of asteroids among numerous sets of astrometry is presented. In particular, numerical methods have been developed for the short-term identification of asteroids at discovery, and for the long-term identification of scarcely observed asteroids over apparitions, a task which has been lacking a robust method until now. The methods are based on the solid foundation of statistical orbital inversion properly taking into account the observational uncertainties, which allows for the detection of practically all correct identifications. Through the use of dimensionality-reduction techniques and efficient data structures, the exact methods have a loglinear, that is, O(nlog(n)), computational complexity, where n is the number of included observation sets. The methods developed are thus suitable for future large-scale surveys which anticipate a substantial increase in the astrometric data rate. Due to the discontinuous nature of asteroid astrometry, separate sets of astrometry must be linked to a common asteroid from the very first discovery detections onwards. The reason for the discontinuity in the observed positions is the rotation of the observer with the Earth as well as the motion of the asteroid and the observer about the Sun. Therefore, the aim of identification is to find a set of orbital elements that reproduce the observed positions with residuals similar to the inevitable observational uncertainty. Unless the astrometric observation sets are linked, the corresponding asteroid is eventually lost as the uncertainty of the predicted positions grows too large to allow successful follow-up. Whereas the presented identification theory and the numerical comparison algorithm are generally applicable, that is, also in fields other than astronomy (e.g., in the identification of space debris), the numerical methods developed for asteroid identification can immediately be applied to all objects on heliocentric orbits with negligible effects due to non-gravitational forces in the time frame of the analysis. The methods developed have been successfully applied to various identification problems. Simulations have shown that the methods developed are able to find virtually all correct linkages despite challenges such as numerous scarce observation sets, astrometric uncertainty, numerous objects confined to a limited region on the celestial sphere, long linking intervals, and substantial parallaxes. Tens of previously unknown main-belt asteroids have been identified with the short-term method in a preliminary study to locate asteroids among numerous unidentified sets of single-night astrometry of moving objects, and scarce astrometry obtained nearly simultaneously with Earth-based and space-based telescopes has been successfully linked despite a substantial parallax. Using the long-term method, thousands of realistic 3-linkages typically spanning several apparitions have so far been found among designated observation sets each spanning less than 48 hours.

Regeneration ecology and woody species diversity in the remnant natural forest at Wondo Genet, Ethiopia – Potential for landscape restoration

Regeneration ecology, diversity of native woody species and its potential for landscape restoration was studied in the remnant natural forest at the College of Forestry and Natural Resources at Wondo Genet, Ethiopia. The type of forest is Afromontane rainforest , with many valuable tree species like Aningeria adolfi-friederici, and it is an important provider of ecological, social and economical services for the population that lives in this area. The study contains two parts, natural regeneration studies (at the natural forest) and interviews with farmers in the nearby village of the remnant patch. The objective of the first part was to investigate the floristic composition, densitiy and regeneration profiles of native woody species in the forest, paying special attention to woody species that are considered the most relevant (socio-economic). The second part provided information on woody species preferred by the farmers and on multiple uses of the adjacent natural forest, it also provided information and analysed perceptions on forest degradation. Systematic plot sampling was used in the forest inventory. Twenty square plots of 20 x 20 m were assessed, with 38 identified woody species (the total number of species was 45), representing 26 families. Of these species 61% were trees, 13% shrubs, 11% lianas and 16% species that could have both life forms. An analysis of natural regeneration of five important tree species in the natural forest showed that Aningeria adolfi-friederici had the best regeneration results. An analysis of population structure (as determined by height classes) of two commercially important woody species in the forest, Aningeria adolfi-friederici and Podocarpus falcatus, showed a marked difference: Aningeria had a typical “reversed J” frequency distribution, while Podocarpus showed very low values in all height classes. Multi dimensional scaling (MDS) was used to map the sample plots according to their similarity in species composition, using the Sørensen quantitative index, coupled with indicator species analysis .Three groups were identified with respective indicator species: Group 1 – Adhatoda schimperiana, Group 2 – Olea hochstetteri , Group 3 – Acacia senegal and Aningeria adolfi-friederici. Thirty questionnaire interviews were conducted with farmers in the village of Gotu Onoma that use the nearby remant forest patch. Their tree preferences were exotic species such as Eucalyptus globulus for construction and fuelwood and Grevillea robusta for shade and fertility. Considering forest land degradation farmers were aware of the problem and suggested that the governmental institutions address the problem by planting more Eucalyptus globulus. The natural forest seemed to have moderate levels of disturbance and it was still floristically diverse. However, the low rate of natural regeneration of Podocarpus falcatus suggested that this species is threatened and must be a priority in conservation actions. Plantations and agroforestry seem to be possible solutions for rehabilitation of the surrounding degraded lands, thereby decreasing the existent pressure in the remnant natural forest.

Open and closed boundaries in large-scale convective dynamos

Earlier work has suggested that large-scale dynamos can reach and maintain equipartition field strengths on a dynamical time scale only if magnetic helicity of the fluctuating field can be shed from the domain through open boundaries. To test this scenario in convection-driven dynamos by comparing results for open and closed boundary conditions. Three-dimensional numerical simulations of turbulent compressible convection with shear and rotation are used to study the effects of boundary conditions on the excitation and saturation level of large-scale dynamos. Open (vertical field) and closed (perfect conductor) boundary conditions are used for the magnetic field. The contours of shear are vertical, crossing the outer surface, and are thus ideally suited for driving a shear-induced magnetic helicity flux. We find that for given shear and rotation rate, the growth rate of the magnetic field is larger if open boundary conditions are used. The growth rate first increases for small magnetic Reynolds number, Rm, but then levels off at an approximately constant value for intermediate values of Rm. For large enough Rm, a small-scale dynamo is excited and the growth rate in this regime increases proportional to Rm^(1/2). In the nonlinear regime, the saturation level of the energy of the mean magnetic field is independent of Rm when open boundaries are used. In the case of perfect conductor boundaries, the saturation level first increases as a function of Rm, but then decreases proportional to Rm^(-1) for Rm > 30, indicative of catastrophic quenching. These results suggest that the shear-induced magnetic helicity flux is efficient in alleviating catastrophic quenching when open boundaries are used. The horizontally averaged mean field is still weakly decreasing as a function of Rm even for open boundaries.