Spatial ecology of a specialist insect herbivore : the leaf-mining moth Tischeria ekebladella on the pedunculate oak Quercus robur

Herbivorous insects comprise a major part of terrestrial biodiversity, and their interactions with their host plants and natural enemies are of vast ecological importance. A large body of research demonstrates that the ecology and evolution of these insects may be affected by trophic interactions, by abiotic influences, and by intraspecific processes, but so far research on these individual aspects has rarely been combined. This thesis uses the leaf-mining moth Tischeria ekebladella and the pedunculate oak (Quercus robur) as a case study to assess how spatial variation in trophic interactions and the physical distribution of host trees jointly affect the distribution, dynamics and evolution of a host-specific herbivore. With respect to habitat quality, Tischeria ekebladella experiences abundant variation at several spatial scales. Most of this variation occurs at small scales notably among leaves and shoots within individual trees. While hypothetically this could cause moths to evolve an ability to select leaves and shoots of high quality, I did not find any coupling between female preference and offspring performance. Based on my studies on temporal variation in resource quality I therefore propose that unpredictable temporal changes in the relative rankings of individual resource units may render it difficult for females to predict the fate of their developing offspring. With respect to intraspecific processes, my results suggest that limited moth dispersal in relation to the spatial distribution of oak trees plays a key role in determining the regional distribution of Tischeria ekebladella. The distribution of the moth is aggregated at the landscape level, where local leaf miner populations are less likely to be present where oaks are scarce. A modelling exercise based on empirical dispersal estimates revealed that the moth population on Wattkast an island in south-western Finland is spatially structured overall, but that the relative importance of local and regional processes on tree-specific moth dynamics varies drastically across the landscape. To conclude, my work in the oak-Tischeria ekebladella system demonstrates that the local abundance and regional distribution of a herbivore may be more strongly influenced by the spatial location of host trees than by their relative quality. Hence, it reveals the importance of considering spatial context in the study of herbivorous insects, and forms a bridge between the classical fields of plant-insect interactions and spatial ecology.

Ecology, population genetics and conservation of the African violet

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

Spatial ecology of an oak-associated herbivore community : A metacommunity of herbivores

Habitat fragmentation is currently affecting many species throughout the world. As a consequence, an increasing number of species are structured as metapopulations, i.e. as local populations connected by dispersal. While excellent studies of metapopulations have accumulated over the past 20 years, the focus has recently shifted from single species to studies of multiple species. This has created the concept of metacommunities, where local communities are connected by the dispersal of one or several of their member species. To understand this higher level of organisation, we need to address not only the properties of single species, but also establish the importance of interspecific interactions. However, studies of metacommunities are so far heavily biased towards laboratory-based systems, and empirical data from natural systems are urgently needed. My thesis focuses on a metacommunity of insect herbivores on the pedunculate oak Quercus robur a tree species known for its high diversity of host-specific insects. Taking advantage of the amenability of this system to both observational and experimental studies, I quantify and compare the importance of local and regional factors in structuring herbivore communities. Most importantly, I contrast the impact of direct and indirect competition, host plant genotype and local adaptation (i.e. local factors) to that of regional processes (as reflected by the spatial context of the local community). As a key approach, I use general theory to generate testable hypotheses, controlled experiments to establish causal relations, and observational data to validate the role played by the pinpointed processes in nature. As the central outcome of my thesis, I am able to relegate local forces to a secondary role in structuring oak-based insect communities. While controlled experiments show that direct competition does occur among both conspecifics and heterospecifics, that indirect interactions can be mediated by both the host plant and the parasitoids, and that host plant genotype may affect local adaptation, the size of these effects is much smaller than that of spatial context. Hence, I conclude that dispersal between habitat patches plays a prime role in structuring the insect community, and that the distribution and abundance of the target species can only be understood in a spatial framework. By extension, I suggest that the majority of herbivore communities are dependent on the spatial structure of their landscape and urge fellow ecologists working on other herbivore systems to either support or refute my generalization.

Phytoplankton ecology of subarctic lakes in Finnish Lapland

Climate is warming and it is especially seen in arctic areas, where the warming trend is expected to be greatest. Arctic freshwater ecosystems, which are a very characteristic feature of the arctic landscape, are especially sensitive to climate change. They could be used as early warning systems, but more information about the ecosystem functioning and responses are needed for proper interpretation of the observations. Phytoplankton species and assemblages could be especially suitable for climate-related studies, since they have short generation times and react rapidly to changes in the environment. In addition, phytoplankton provides a good tool for lake classifications, since different species have different requirements and tolerance ranges for various environmental factors. The use of biological indicators is especially useful in arctic areas, were many of the chemical factors commonly fall under the detection limit and therefore do not provide much information about the environment. This work brings new information about species distribution and dynamics of arctic freshwater phytoplankton in relation to environmental factors. The phytoplankton of lakes in Finnish Lapland and other European high-altitude or high-latitude areas were compared. Most lakes were oligotrophic and dominated by flagellated species belonging to chrysophytes, cryptophytes and dinoflagellates. In Finnish Lapland cryptophytes were of less importance, whereas desmids had high species richness in many of the lakes. In Pan-European scale, geographical and catchment-related factors were explaining most of the differences in species distributions between different districts, whereas lake water chemistry (especially conductivity, SiO2 and pH) was most important regionally. Seasonal and interannual variation of phytoplankton was studied in subarctic Lake Saanajärvi. Characteristic phytoplankton species in this oligotrophic, dimictic lake belonged mainly to chrysophytes and diatoms. The maximum phytoplankton biomass in Lake Saanajärvi occurs during autumn, while spring biomass is very low. During years with heavy snow cover the lake suffers from pH drop caused by melt waters, but the effects of this acid pulse are restricted to surface layers and last for a relatively short period. In addition to some chemical parameters (mainly Ca and nutrients), length of the mixing cycle and physical factors such as lake water temperature and thermal stability of water column had major impact on phytoplankton dynamics. During a year with long and strong thermal stability, the phytoplankton community developed towards an equilibrium state, with heavy dominance of only a few taxa for a longer period of time. During a year with higher windiness and less thermal stability, the species composition was more diverse and species with different functional strategies were able to occur simultaneously. The results of this work indicate that although arctic lakes in general share many common features concerning their catchment and water chemistry, large differences in biological features can be found even in a relatively small area. Most likely the lakes with very different algal flora do not respond in a similar way to differences in the environmental factors, and more information about specific arctic lake types is needed. The results also show considerable year to year differences in phytoplankton species distribution and dynamics, and these changes are most likely linked to climatic factors.

Designing Dense Radar Networks for Weather Surveillance in Tropical Regions.

Design considerations are presented for a dense weather radar network to support multiple services including aviation. Conflicts, tradeoffs and optimization issues in the context of operation in a tropical region are brought out. The upcoming Indian radar network is used as a case study. Algorithms for data mosaicing are briefly outlined.

Analyzing changes in macro-level drivers of country-specific emissions : The role of consumption, technology and trade

Industrial ecology is an important field of sustainability science. It can be applied to study environmental problems in a policy relevant manner. Industrial ecology uses ecosystem analogy; it aims at closing the loop of materials and substances and at the same time reducing resource consumption and environmental emissions. Emissions from human activities are related to human interference in material cycles. Carbon (C), nitrogen (N) and phosphorus (P) are essential elements for all living organisms, but in excess have negative environmental impacts, such as climate change (CO2, CH4 N2O), acidification (NOx) and eutrophication (N, P). Several indirect macro-level drivers affect emissions change. Population and affluence (GDP/capita) often act as upward drivers for emissions. Technology, as emissions per service used, and consumption, as economic intensity of use, may act as drivers resulting in a reduction in emissions. In addition, the development of country-specific emissions is affected by international trade. The aim of this study was to analyse changes in emissions as affected by macro-level drivers in different European case studies. ImPACT decomposition analysis (IPAT identity) was applied as a method in papers I III. The macro-level perspective was applied to evaluate CO2 emission reduction targets (paper II) and the sharing of greenhouse gas emission reduction targets (paper IV) in the European Union (EU27) up to the year 2020. Data for the study were mainly gathered from official statistics. In all cases, the results were discussed from an environmental policy perspective. The development of nitrogen oxide (NOx) emissions was analysed in the Finnish energy sector during a long time period, 1950 2003 (paper I). Finnish emissions of NOx began to decrease in the 1980s as the progress in technology in terms of NOx/energy curbed the impact of the growth in affluence and population. Carbon dioxide (CO2) emissions related to energy use during 1993 2004 (paper II) were analysed by country and region within the European Union. Considering energy-based CO2 emissions in the European Union, dematerialization and decarbonisation did occur, but not sufficiently to offset population growth and the rapidly increasing affluence during 1993 2004. The development of nitrogen and phosphorus load from aquaculture in relation to salmonid consumption in Finland during 1980 2007 was examined, including international trade in the analysis (paper III). A regional environmental issue, eutrophication of the Baltic Sea, and a marginal, yet locally important source of nutrients was used as a case. Nutrient emissions from Finnish aquaculture decreased from the 1990s onwards: although population, affluence and salmonid consumption steadily increased, aquaculture technology improved and the relative share of imported salmonids increased. According to the sustainability challenge in industrial ecology, the environmental impact of the growing population size and affluence should be compensated by improvements in technology (emissions/service used) and with dematerialisation. In the studied cases, the emission intensity of energy production could be lowered for NOx by cleaning the exhaust gases. Reorganization of the structure of energy production as well as technological innovations will be essential in lowering the emissions of both CO2 and NOx. Regarding the intensity of energy use, making the combustion of fuels more efficient and reducing energy use are essential. In reducing nutrient emissions from Finnish aquaculture to the Baltic Sea (paper III) through technology, limits of biological and physical properties of cultured fish, among others, will eventually be faced. Regarding consumption, salmonids are preferred to many other protein sources. Regarding trade, increasing the proportion of imports will outsource the impacts. Besides improving technology and dematerialization, other viewpoints may also be needed. Reducing the total amount of nutrients cycling in energy systems and eventually contributing to NOx emissions needs to be emphasized. Considering aquaculture emissions, nutrient cycles can be partly closed through using local fish as feed replacing imported feed. In particular, the reduction of CO2 emissions in the future is a very challenging task when considering the necessary rates of dematerialisation and decarbonisation (paper II). Climate change mitigation may have to focus on other greenhouse gases than CO2 and on the potential role of biomass as a carbon sink, among others. The global population is growing and scaling up the environmental impact. Population issues and growing affluence must be considered when discussing emission reductions. Climate policy has only very recently had an influence on emissions, and strong actions are now called for climate change mitigation. Environmental policies in general must cover all the regions related to production and impacts in order to avoid outsourcing of emissions and leakage effects. The macro-level drivers affecting changes in emissions can be identified with the ImPACT framework. Statistics for generally known macro-indicators are currently relatively well available for different countries, and the method is transparent. In the papers included in this study, a similar method was successfully applied in different types of case studies. Using transparent macro-level figures and a simple top-down approach are also appropriate in evaluating and setting international emission reduction targets, as demonstrated in papers II and IV. The projected rates of population and affluence growth are especially worth consideration in setting targets. However, sensitivities in calculations must be carefully acknowledged. In the basic form of the ImPACT model, the economic intensity of consumption and emission intensity of use are included. In seeking to examine consumption but also international trade in more detail, imports were included in paper III. This example demonstrates well how outsourcing of production influences domestic emissions. Country-specific production-based emissions have often been used in similar decomposition analyses. Nevertheless, trade-related issues must not be ignored.

Genetic analysis of chromosomal regions 2q33, 7q32 and 19q13 in multiple sclerosis susceptibility

Multiple sclerosis (MS) is an immune-mediated demyelinating disorder of the central nervous system (CNS) affecting 0.1-0.2% of Northern European descent population. MS is considered to be a multifactorial disease, both environment and genetics play a role in its pathogenesis. Despite several decades of intense research, the etiological and pathogenic mechanisms underlying MS remain still largely unknown and no curative treatment exists. The genetic architecture underlying MS is complex with multiple genes involved. The strongest and the best characterized predisposing genetic factors for MS are located, as in other immune-mediated diseases, in the major histocompatibility complex (MHC) on chromosome 6. In humans MHC is called human leukocyte antigen (HLA). Alleles of the HLA locus have been found to associate strongly with MS and remained for many years the only consistently replicable genetic associations. However, recently other genes located outside the MHC region have been proposed as strong candidates for susceptibility to MS in several studies. In this thesis a new genetic locus located on chromosome 7q32, interferon regulatory factor 5 (IRF5), was identified in the susceptibility to MS. In particular, we found that common variation of the gene was associated with the disease in three different populations, Spanish, Swedish and Finnish. We also suggested a possible functional role for one of the risk alleles with impact on the expression of the IRF5 locus. Previous studies have pointed out a possible role played by chromosome 2q33 in the susceptibility to MS and other autoimmune disorders. The work described here also investigated the involvement of this chromosomal region in MS predisposition. After the detection of genetic association with 2q33 (article-1), we extended our analysis through fine-scale single nucleotide polymorphism (SNP) mapping to define further the contribution of this genomic area to disease pathogenesis (article-4). We found a trend (p=0.04) for association to MS with an intronic SNP located in the inducible T-cell co-stimulator (ICOS) gene, an important player in the co-stimulatory pathway of the immune system. Expression analysis of ICOS revealed a novel, previously uncharacterized, alternatively spliced isoform, lacking the extracellular domain that is needed for ligand binding. The stability of the newly-identified transcript variant and its subcellular localization were analyzed. These studies indicated that the novel isoform is stable and shows different subcellular localization as compared to full-length ICOS. The novel isoform might have a regulatory function, but further studies are required to elucidate its function. Chromosome 19q13 has been previously suggested as one of the genomic areas involved in MS predisposition. In several populations, suggestive linkage signals between MS predisposition and 19q13 have been obtained. Here, we analysed the role of allelic variation in 19q13 by family based association analysis in 782 MS families collected from Finland. In this dataset, we were not able to detect any statistically significant associations, although several previously suggested markers were included to the analysis. Replication of the previous findings on the basis of linkage disequilibrium between marker allele and disease/risk allele appears notoriously difficult because of limitations such as allelic heterogeneity. Re-sequencing based approaches may be required for elucidating the role of chromosome 19q13 with MS. This thesis has resulted in the identification of a new MS susceptibility locus (IRF5) previously associated with other inflammatory or autoimmune disorders, such as SLE. IRF5 is one of the mediators of interferons biological function. In addition to providing new insight in the possible pathogenetic pathway of the disease, this finding suggests that there might be common mechanisms between different immune-mediated disorders. Furthermore the work presented here has uncovered a novel isoform of ICOS, which may play a role in regulatory mechanisms of ICOS, an important mediator of lymphocyte activation. Further work is required to uncover its functions and possible involvement of the ICOS locus in MS susceptibility.

A study on the variations of optical and physical properties of aerosols over a tropical semi-arid station during grassland fire

The present paper records the results of a case study on the impact of an extensive grassland fire on the physical and optical properties of aerosols at a semi-arid station in southern India for the first time from ground based measurements using a MICROTOPS-II sunphotometer, an aethalometer and a quartz crystal microbalance impactor (QCM). Observations revealed a substantial increase in aerosol optical depth (AOD) at all wavelengths during burning days compared to normal days. High AOD values observed at shorter wavelengths suggest the dominance of accumulation mode particle loading over the study area. Daily mean aerosol size spectra shows, most of the time, power-law distribution. To characterize AOD, the Angstrom parameters (i.e., alpha and beta) were used. Wavelength exponent (1.38) and turbidity coefficient (0.21) are high during burning days compared to normal days, thereby suggesting an increase in accumulation mode particle loading. Aerosol size distribution suggested dominance of accumulation mode particle loading during burning days compared to normal days. A significant positive correlation was observed between AOD at 500 mn and water vapour and negative correlation between AOD at 500 nm and wind speed for burning and non-burning days. Diurnal variations of black carbon (BC) aerosol mass concentrations increased by a factor of similar to 2 in the morning and afternoon hours during burning period compared to normal days.

Spatial and vertical heterogeneities in aerosol properties over oceanic regions around India: Implications for radiative forcing

The influence of atmospheric aerosols on Earth's radiation budget and hence climate, though well recognized and extensively investigated in recent years, remains largely uncertain mainly because of the large spatio-temporal heterogeneity and the lack of data with adequate resolution. To characterize this diversity, a major multi-platform field campaign ICARB (Integrated Campaign for Aerosols, gases and Radiation Budget) was carried out during the pre-monsoon period of 2006 over the Indian landmass and surrounding oceans, which was the biggest such campaign ever conducted over this region. Based on the extensive and concurrent measurements of the optical and physical properties of atmospheric aerosols during ICARB, the spatial distribution of aerosol radiative forcing was estimated over the entire Bay of Bengal (BoB), northern Indian Ocean and Arabian Sea (AS) as well as large spatial variations within these regions. Besides being considerably lower than the mean values reported earlier for this region, our studies have revealed large differences in the forcing components between the BoB and the AS. While the regionally averaged aerosol-induced atmospheric forcing efficiency was 31 +/- 6 W m(-2) tau(-1) for the BoB, it was only similar to 18 +/- 7 W m(-2) tau(-1) for the AS. Airborne measurements revealed the presence of strong, elevated aerosol layers even over the oceans, leading to vertical structures in the atmospheric forcing, resulting in significant warming in the lower troposphere. These observations suggest serious climate implications and raise issues ranging from the impact of aerosols on vertical thermal structure of the atmospheric and hence cloud formation processes to monsoon circulation.

Water balance modelling in a tropical watershed under deciduous forest (Mule Hole, India): Regolith matric storage buffers the groundwater recharge process

Accurate estimations of water balance are needed in semi-arid and sub-humid tropical regions, where water resources are scarce compared to water demand. Evapotranspiration plays a major role in this context, and the difficulty to quantify it precisely leads to major uncertainties in the groundwater recharge assessment, especially in forested catchments. In this paper, we propose to assess the importance of deep unsaturated regolith and water uptake by deep tree roots on the groundwater recharge process by using a lumped conceptual model (COMFORT). The model is calibrated using a 5 year hydrological monitoring of an experimental watershed under dry deciduous forest in South India (Mule Hole watershed). The model was able to simulate the stream discharge as well as the contrasted behaviour of groundwater table along the hillslope. Water balance simulated for a 32 year climatic time series displayed a large year-to-year variability, with alternance of dry and wet phases with a time period of approximately 14 years. On an average, input by the rainfall was 1090 mm year(-1) and the evapotranspiration was about 900 mm year(-1) out of which 100 mm year(-1) was uptake from the deep saprolite horizons. The stream flow was 100 mm year(-1) while the groundwater underflow was 80 mm year(-1). The simulation results suggest that (i) deciduous trees can uptake a significant amount of water from the deep regolith, (ii) this uptake, combined with the spatial variability of regolith depth, can account for the variable lag time between drainage events and groundwater rise observed for the different piezometers and (iii) water table response to recharge is buffered due to the long vertical travel time through the deep vadose zone, which constitutes a major water reservoir. This study stresses the importance of long term observations for the understanding of hydrological processes in tropical forested ecosystems. (C) 2009 Elsevier B.V. All rights reserved.

A theoretical model for the magnetic helicity of solar active regions

Active regions on the solar surface are known to possess magnetic helicity, which is predominantly negative in the northern hemisphere and positive in the southern hemisphere. Choudhuri et al. [Choudhuri, A.R. On the connection between mean field dynamo theory and flux tubes. Solar Phys. 215, 31–55, 2003] proposed that the magnetic helicity arises due to the wrapping up of the poloidal field of the convection zone around rising flux tubes which form active regions. Choudhuri [Choudhuri, A.R., Chatterjee, P., Nandy, D. Helicity of solar active regions from a dynamo model. ApJ 615, L57–L60, 2004] used this idea to calculate magnetic helicity from their solar dynamo model. Apart from getting broad agreements with observational data, they also predict that the hemispheric helicity rule may be violated at the beginning of a solar cycle. Chatterjee et al. [Chatterjee, P., Choudhuri, A.R., Petrovay, K. Development of twist in an emerging magnetic flux tube by poloidal field accretion. A&A 449, 781–789, 2006] study the penetration of the wrapped poloidal field into the rising flux tube due to turbulent diffusion using a simple 1-d model. They find that the extent of penetration of the wrapped field will depend on how weak the magnetic field inside the rising flux tube becomes before its emergence. They conclude that more detailed observational data will throw light on the physical conditions of flux tubes just before their emergence to the photosphere.

Systematic study on crystal-contact engineering of diphthine synthase: influence of mutations at crystal-packing regions on X-ray diffraction quality

It is well known that protein crystallizability can be influenced by site-directed mutagenesis of residues on the molecular surface of proteins, indicating that the intermolecular interactions in crystal-packing regions may play a crucial role in the structural regularity at atomic resolution of protein crystals. Here, a systematic examination was made of the improvement in the diffraction resolution of protein crystals on introducing a single mutation of a crystal-packing residue in order to provide more favourable packing interactions, using diphthine synthase from Pyrococcus horikoshii OT3 as a model system. All of a total of 21 designed mutants at 13 different crystal-packing residues yielded almost isomorphous crystals from the same crystallization conditions as those used for the wild-type crystals, which diffracted X-rays to 2.1 angstrom resolution. Of the 21 mutants, eight provided crystals with an improved resolution of 1.8 angstrom or better. Thus, it has been clarified that crystal quality can be improved by introducing a suitable single mutation of a crystal-packing residue. In the improved crystals, more intimate crystal-packing interactions than those in the wild-type crystal are observed. Notably, the mutants K49R and T146R yielded crystals with outstandingly improved resolutions of 1.5 and 1.6 angstrom, respectively, in which a large-scale rearrangement of packing interactions was unexpectedly observed despite the retention of the same isomorphous crystal form. In contrast, the mutants that provided results that were in good agreement with the designed putative structures tended to achieve only moderate improvements in resolution of up to 1.75 angstrom. These results suggest a difficulty in the rational prediction of highly effective mutations in crystal engineering.