A Finite Variable Difference Relaxation Scheme for hyperbolic-parabolic equations

Using the framework of a new relaxation system, which converts a nonlinear viscous conservation law into a system of linear convection-diffusion equations with nonlinear source terms, a finite variable difference method is developed for nonlinear hyperbolic-parabolic equations. The basic idea is to formulate a finite volume method with an optimum spatial difference, using the Locally Exact Numerical Scheme (LENS), leading to a Finite Variable Difference Method as introduced by Sakai [Katsuhiro Sakai, A new finite variable difference method with application to locally exact numerical scheme, journal of Computational Physics, 124 (1996) pp. 301-308.], for the linear convection-diffusion equations obtained by using a relaxation system. Source terms are treated with the well-balanced scheme of Jin [Shi Jin, A steady-state capturing method for hyperbolic systems with geometrical source terms, Mathematical Modeling Numerical Analysis, 35 (4) (2001) pp. 631-645]. Bench-mark test problems for scalar and vector conservation laws in one and two dimensions are solved using this new algorithm and the results demonstrate the efficiency of the scheme in capturing the flow features accurately.

Population dynamics of blue mussels in a variable environment at the edge of their range

Spatial and temporal variation in the abundance of species can often be ascribed to spatial and temporal variation in the surrounding environment. Knowledge of how biotic and abiotic factors operate over different spatial and temporal scales in determining distribution, abundance, and structure of populations lies at the heart of ecology. The major part of the current ecological theory stems from studies carried out in central parts of the distributional range of species, whereas knowledge of how marginal populations function is inadequate. Understanding how marginal populations, living at the edge of their range, function is however in a key position to advance ecology and evolutionary biology as scientific disciplines. My thesis focuses on the factors affecting dynamics of marginal populations of blue mussels (Mytilus edulis) living close to their tolerance limits with regard to salinity. The thesis aims to highlight the dynamics at the edge of the range and contrast these with dynamics in more central parts of the range in order to understand the potential interplay between the central and the marginal part in the focal system. The objectives of the thesis are approached by studies on: (1) factors affecting regional patterns of the species, (2) long-term temporal dynamics of the focal species spaced along a regional salinity gradient, (3) selective predation by increasing populations of roach (Rutilus rutilus) when feeding on their main food item, the blue mussel, (4) the primary and secondary effects of local wave exposure gradients and (5) the role of small-scale habitat heterogeneity as determinants of large-scale pattern. The thesis shows that populations of blue mussels are largely determined by large scale changes in sea water salinity, affecting mainly recruitment success and longevity of local populations. In opposite to the traditional view, the thesis strongly indicate that vertebrate predators strongly affect abundance and size structure of blue mussel populations, and that the role of these predators increases towards the margin where populations are increasingly top-down controlled. The thesis also indicates that the positive role of biogenic habitat modifiers increases towards the marginal areas, where populations of blue mussels are largely recruitment limited. Finally, the thesis shows that local blue mussel populations are strongly dependent on high water turbulence, and therefore, dense populations are constrained to offshore habitats. Finally, the thesis suggests that ongoing sedimentation of rocky shores is detrimental for the species, affecting recruitment success and post-recruit survival, pushing stable mussel beds towards offshore areas. Ongoing large scale changes in the Baltic Sea, especially dilution processes with attendant effects, are predicted to substantially contract the distributional range of the mussel, but also affect more central populations. The thesis shows that in order to understand the functioning of marginal populations, research should (1) strive for multi-scale approaches in order to link ecosystem patterns with ecosystem processes, and (2) challenge the prevailing tenets that origin from research carried out in central areas that may not be valid at the edge.

Alternate pathways of the early secretory route in yeast

The diversity of functions of eukaryotic cells is preserved by enclosing different enzymatic activities into membrane-bound organelles. Separation of exocytic proteins from those which remain in the endoplasmic reticulum (ER) casts the foundation for correct compartmentalization. The secretory pathway, starting from the ER membrane, operates by the aid of cytosolic coat proteins (COPs). In anterograde transport, polymerization of the COPII coat on the ER membrane is essential for the ER exit of proteins. Polymerization of the COPI coatomer on the cis-Golgi membrane functions for the retrieval of proteins from the Golgi for repeated use in the ER. The COPII coat is formed by essential proteins; Sec13/31p and Sec23/24p have been thought to be indispensable for the ER exit of all exocytic proteins. However, we found that functional Sec13p was not required for the ER exit of yeast endogenous glycoprotein Hsp150 in the yeast Saccharomyces cerevisiae. Hsp150 turned out to be an ATP phosphatase. ATP hydrolysis by a Walker motif located in the C-terminal domain of Hsp150 was an active mediator for the Sec13p and Sec24p independent ER exit. Our results suggest that in yeast cells a fast track transport route operates in parallel with the previously described cisternal maturation route of the Golgi. The fast track is used by Hsp150 with the aid of its C-terminal ATPase activity at the ER-exit. Hsp150 is matured with a half time of less than one minute. The cisternal maturation track is several-fold slower and used by other exocytic proteins studied so far. Operative COPI coat is needed for ER exit by a subset of proteins but not by Hsp150. We located a second active determinant to the Hsp150 polypeptide s N-terminal portion that guided also heterologous fusion proteins out of the ER in COPII coated vesicles under non-functional COPI conditions for several hours. Our data indicate that ER exit is a selective, receptor-mediated event, not a bulk flow. Furthermore, it suggests the existence of another retrieval pathway for essential reusable components, besides the COPI-operated retrotransport route. Additional experiments suggest that activation of the COPI primer, ADP ribosylation factor (ARF), is essential also for Hsp150 transport. Moreover, it seemed that a subset of proteins directly needed activated ARF in the anterograde transport to complete the ER exit. Our results indicate that coat structures and transport routes are more variable than it has been imagined.

U12-type Spliceosome: Localization and Effects of Splicing Efficiency on Gene Expression

The removal of non-coding sequences, introns, is an essential part of messenger RNA processing. In most metazoan organisms, the U12-type spliceosome processes a subset of introns containing highly conserved recognition sequences. U12-type introns constitute less than 0,5% of all introns and reside preferentially in genes related to information processing functions, as opposed to genes encoding for metabolic enzymes. It has previously been shown that the excision of U12-type introns is inefficient compared to that of U2-type introns, supporting the model that these introns could provide a rate-limiting control for gene expression. The low efficiency of U12-type splicing is believed to have important consequences to gene expression by limiting the production of mature mRNAs from genes containing U12-type introns. The inefficiency of U12-type splicing has been attributed to the low abundance of the components of the U12-type spliceosome in cells, but this hypothesis has not been proven. The aim of the first part of this work was to study the effect of the abundance of the spliceosomal snRNA components on splicing. Cells with a low abundance of the U12-type spliceosome were found to inefficiently process U12-type introns encoded by a transfected construct, but the expression levels of endogenous genes were not found to be affected by the abundance of the U12-type spliceosome. However, significant levels of endogenous unspliced U12-type intron-containing pre-mRNAs were detected in cells. Together these results support the idea that U12-type splicing may limit gene expression in some situations. The inefficiency of U12-type splicing has also promoted the idea that the U12-type spliceosome may control gene expression, limiting the mRNA levels of some U12-type intron-containing genes. While the identities of the primary target genes that contain U12-type introns are relatively well known, little has previously been known about the downstream genes and pathways potentially affected by the efficiency of U12-type intron processing. Here, the effects of U12-type splicing efficiency on a whole organism were studied in a Drosophila line with a mutation in an essential U12-type spliceosome component. Genes containing U12-type introns showed variable gene-specific responses to the splicing defect, which points to variation in the susceptibility of different genes to changes in splicing efficiency. Surprisingly, microarray screening revealed that metabolic genes were enriched among downstream effects, and that the phenotype could largely be attributed to one U12-type intron-containing mitochondrial gene. Gene expression control by the U12-type spliceosome could thus have widespread effects on metabolic functions in the organism. The subcellular localization of the U12-type spliceosome components was studied as a response to a recent dispute on the localization of the U12-type spliceosome. All components studied were found to be nuclear indicating that the processing of U12-type introns occurs within the nucleus, thus clarifying a question central to the field. The results suggest that the U12-type spliceosome can limit the expression of genes that contain U12-type introns in a gene-specific manner. Through its limiting role in pre-mRNA processing, the U12-type splicing activity can affect specific genetic pathways, which in the case of Drosophila are involved in metabolic functions.

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.

Differentiating between pollutants build-up on roads and roofs: Significance of roofs as a stormwater pollutant source

Impervious surfaces in an urban catchment are the primary stormwater pollutant contributing areas. Appropriate treatment of stormwater runoff from these impervious surfaces is essential to safeguard the urban water environment. While urban roads have received significant research attention in this regard, roofs have not been well investigated. Key pollutant processes such as build-up on roads and roofs can be different due to the different surface characteristics. This entails different treatment strategies being needed for road and roofs. The research study characterized roof pollutants build-up by differentiating with road surfaces. It was noted that pollutants are more highly concentrated on particles and particularly finer particles in the case of roof surfaces, compared to road surfaces. Additionally, pollutants built-up on roof surfaces tend to be relatively more variable from one day to another in terms of pollutant loads. These results highlight the significance of roofs as a stormwater pollutant source and the important need for a specific stormwater treatment strategy rather than the application of a combined approach for treating stormwater runoff from both, roads and roofs.

Sostdc1 Plays an Essential Role in Mammalian Tooth Patterning : Insight into the Rodent Dental Evolution

This thesis work focuses on the role of TGF-beta family antagonists during the development of mouse dentition. Tooth develops through an interaction between the dental epithelium and underlying neural crest derived mesenchyme. The reciprocal signaling between these tissues is mediated by soluble signaling molecules and the balance between activatory and inhibitory signals appears to be essential for the pattern formation. We showed the importance of Sostdc1 in the regulation of tooth shape and number. The absence of Sostdc1 altered the molar cusp patterning and led to supernumerary tooth formation both in the molar and incisor region. We showed that initially, Sostdc1 expression is in the mesenchyme, suggesting that dental mesenchyme may limit supernumerary tooth induction. We tested this in wild-type incisors by minimizing the amount of mesenchymal tissue surrounding the incisor tooth germs prior to culture in vitro. The cultured teeth phenocopied the extra incisor phenotype of the Sostdc1-deficient mice. Furthermore, we showed that minimizing the amount of dental mesenchyme in cultured Sostdc1-deficient incisors caused the formation of additional de novo incisors that resembled the successional incisor development resulting from activated Wnt signaling. Sostdc1 seemed to be able to inhibit both mesenchymal BMP4 and epithelial canonical Wnt signaling, which thus allows Sostdc1 to restrict the enamel knot size and regulate the tooth shape and number. Our work emphasizes the dual role for the tooth mesenchyme as a suppressor as well as an activator during tooth development. We found that the placode, forming the thick mouse incisor, is prone to disintegration during initiation of tooth development. The balance between two mesenchymal TGF-beta family signals, BMP4 and Activin is essential in this regulation. The inhibition of BMP4 or increase in Activin signaling led to the splitting of the large incisor placode into two smaller placodes resulting in thin incisors. These two signals appeared to have different effects on tooth epithelium and the analysis of the double null mutant mice lacking Sostdc1 and Follistatin indicated that these TGF-beta inhibitors regulate the mutual balance of BMP and Activin in vivo. In addition, this work provides an alternative explanation for the issue of incisor identity published in Science by Tucker et al. in 1998 and proposes that the molar like morphology that can be obtained by inhibiting BMP signaling is due to partial splitting of the incisor placodes and not due to change in tooth identity from the incisor to the molar. This thesis work presents possible molecular mechanisms that may have modified the mouse dental pattern during evolution leading to the typical rodent dentition of modern mouse. The rodent dentition is specialized for gnawing and consists of two large continuously growing incisors and toothless diastema region separating the molars and incisors. The ancestors of rodents had higher number of more slender incisors together with canines and premolars. Additionally, murine rodents, which include the mouse, have lost their ability for tooth replacement. This work has revealed that the inhibitory molecules appear to play a role in the tooth number suppression by delineating the spatial and temporal action of the inductive signals. The results suggest that Sostdc1 plays an essential role in several stages of tooth development through the regulation of both the BMP and Wnt pathway. The work shows a dormant sequential tooth forming potential present in wild type mouse incisor region and gives a new perspective on tooth suppression by dental mesenchyme. It reveals as well a novel mechanism to create a large mouse incisor through the regulation of mesenchymal balance between inductive and inhibitory signals.

The assessment of ecological condition in south-east Queensland, Australia: An evaluation of reliability across variable environments and surrogate efficacy for biodiversity values

Multimetric ecological condition assessment has become an important biodiversity management tool. This study was the first to examine the reliability of these ecological surrogates across variable environments, and the implications for surrogate efficacy. It was demonstrated that through strategic application and design of the multimetric ecological condition index, the effects of environmental gradients and disturbance regimes can be mitigated, and that ecological condition assessment may serve as a scientifically rigorous approach for conservation planning.

Predicting the temporal response of seagrass meadows to dredging using Dynamic Bayesian Networks

Predicting temporal responses of ecosystems to disturbances associated with industrial activities is critical for their management and conservation. However, prediction of ecosystem responses is challenging due to the complexity and potential non-linearities stemming from interactions between system components and multiple environmental drivers. Prediction is particularly difficult for marine ecosystems due to their often highly variable and complex natures and large uncertainties surrounding their dynamic responses. Consequently, current management of such systems often rely on expert judgement and/or complex quantitative models that consider only a subset of the relevant ecological processes. Hence there exists an urgent need for the development of whole-of-systems predictive models to support decision and policy makers in managing complex marine systems in the context of industry based disturbances. This paper presents Dynamic Bayesian Networks (DBNs) for predicting the temporal response of a marine ecosystem to anthropogenic disturbances. The DBN provides a visual representation of the problem domain in terms of factors (parts of the ecosystem) and their relationships. These relationships are quantified via Conditional Probability Tables (CPTs), which estimate the variability and uncertainty in the distribution of each factor. The combination of qualitative visual and quantitative elements in a DBN facilitates the integration of a wide array of data, published and expert knowledge and other models. Such multiple sources are often essential as one single source of information is rarely sufficient to cover the diverse range of factors relevant to a management task. Here, a DBN model is developed for tropical, annual Halophila and temperate, persistent Amphibolis seagrass meadows to inform dredging management and help meet environmental guidelines. Specifically, the impacts of capital (e.g. new port development) and maintenance (e.g. maintaining channel depths in established ports) dredging is evaluated with respect to the risk of permanent loss, defined as no recovery within 5 years (Environmental Protection Agency guidelines). The model is developed using expert knowledge, existing literature, statistical models of environmental light, and experimental data. The model is then demonstrated in a case study through the analysis of a variety of dredging, environmental and seagrass ecosystem recovery scenarios. In spatial zones significantly affected by dredging, such as the zone of moderate impact, shoot density has a very high probability of being driven to zero by capital dredging due to the duration of such dredging. Here, fast growing Halophila species can recover, however, the probability of recovery depends on the presence of seed banks. On the other hand, slow growing Amphibolis meadows have a high probability of suffering permanent loss. However, in the maintenance dredging scenario, due to the shorter duration of dredging, Amphibolis is better able to resist the impacts of dredging. For both types of seagrass meadows, the probability of loss was strongly dependent on the biological and ecological status of the meadow, as well as environmental conditions post-dredging. The ability to predict the ecosystem response under cumulative, non-linear interactions across a complex ecosystem highlights the utility of DBNs for decision support and environmental management.

Allostery and conformational free energy changes in human tryptophanyl-tRNA synthetase from essential dynamics and structure networks

The interdependence of the concept of allostery and enzymatic catalysis, and they being guided by conformational mobility is gaining increased prominence. However, to gain a molecular level understanding of llostery and hence of enzymatic catalysis, it is of utter importance that the networks of amino acids participating in allostery be deciphered. Our lab has been exploring the methods of network analysis combined with molecular dynamics simulations to understand allostery at molecular level. Earlier we had outlined methods to obtain communication paths and then to map the rigid/flexible regions of proteins through network parameters like the shortest correlated paths, cliques, and communities. In this article, we advance the methodology to estimate the conformational populations in terms of cliques/communities formed by interactions including the side-chains and then to compute the ligand-induced population shift. Finally, we obtain the free-energy landscape of the protein in equilibrium, characterizing the free-energy minima accessed by the protein complexes. We have chosen human tryptophanyl-tRNA synthetase (hTrpRS), a protein esponsible for charging tryptophan to its cognate tRNA during protein biosynthesis for this investigation. This is a multidomain protein exhibiting excellent allosteric communication. Our approach has provided valuable structural as well as functional insights into the protein. The methodology adopted here is highly generalized to illuminate the linkage between protein structure networks and conformational mobility involved in the allosteric mechanism in any protein with known structure.

Essential Thrombocythaemia : Diagnosis, Prognostic Aspects, and the Outcome of Finnish patients

Essential thrombocythaemia (ET) is a myeloproliferative disease (MPD) characterized by thrombocytosis, i.e. a constant elevation of platelet count. Thrombocytosis may appear in MPDs (ET, polycythaemia vera, chronic myeloid leukaemia, myelofibrosis) and as a reactive phenomenon. The differential diagnosis of thrombocytosis is important, because the clinical course, need of therapy, and prognosis are different in patients with MPDs and in those with reactive thrombocytosis. ET patients may remain asymptomatic for years, but serious thrombohaemorrhagic and pregnancy-related complications may occur. The complications are difficult to predict. The aims of the present study were to evaluate the diagnostic findings, clinical course, and prognostic factors of ET. The present retrospective study consists of 170 ET patients. Two thirds had a platelet count < 1000 x 109/l. The diagnosis was supported by an increased number of megakaryocytes with an abnormal morphology in a bone marrow aspirate, aggregation defects in platelet function studies, and the presence of spontaneous erythroid and/or megakaryocytic colony formation in in vitro cultures of haematopoietic progenitors. About 70 % of the patients had spontaneous colony formation, while about 30 % had a normal growth pattern. Only a fifth of the patients remained asymptomatic. Half had a major thrombohaemorrhagic complication. The proportion of the patients suffering from thrombosis was as high as 45 %. About a fifth had major bleedings. Half of the patients had microvascular symptoms. Age over 60 years increased the risk of major bleedings, but the occurrence of thrombotic complications was similar in all age groups. Male gender, smoking in female patients, the presence of any spontaneous colony formation, and the presence of spontaneous megakaryocytic colony formation in younger patients were identified as risk factors for thrombosis. Pregnant ET patients had an increased risk of complications. Forty-five per cent of the pregnancies were complicated and 38 % of them ended in stillbirth. Treatment with acetylsalicylic acid alone or in combination with platelet lowering drugs improved the outcome of the pregnancy. The present findings about risk factors in ET as well as treatment outcome in the pregnancies of ET patients should be taken into account when planning treatment strategies for Finnish patients.

Characterization of single-stranded DNA-binding proteins from Mycobacteria. The carboxyl-terminal of domain of SSB is essential for stable association with its cognate RecA protein

Single-stranded DNA-binding proteins (SSB) play an important role in most aspects of DNA metabolism including DNA replication, repair, and recombination. We report here the identification and characterization of SSB proteins of Mycobacterium smegmatis and Mycobacterium tuberculosis. Sequence comparison of M. smegmatis SSB revealed that it is homologous to M. tuberculosis SSB, except for a small spacer connecting the larger amino-terminal domain with the extreme carboxyl-terminal tail. The purified SSB proteins of mycobacteria bound single-stranded DNA with high affinity, and the association and dissociation constants were similar to that of the prototype SSB. The proteolytic signatures of free and bound forms of SSB proteins disclosed that DNA binding was associated with structural changes at the carboxyl-terminal domain. Significantly, SSB proteins from mycobacteria displayed high affinity for cognate RecA, whereas Escherichia coli SSB did not under comparable experimental conditions. Accordingly, SSB and RecA were coimmunoprecipitated from cell lysates, further supporting an interaction between these proteins in vivo. The carboxyl-terminal domain of M. smegmatis SSB, which is not essential for interaction with ssDNA, is the site of binding of its cognate RecA. These studies provide the first evidence for stable association of eubacterial SSB proteins with their cognate RecA, suggesting that these two proteins might function together during DNA repair and/or recombination.

Apparent shortening of the Csp(3)-Csp(3) bond analysed via a variable-temperature X-ray diffraction study in racemic 1,1 '-binaphthalene-2,2 '-diyl diethyl bis(carbonate)

Crystal structure determination at room temperature [292 (2) K] of racemic 1,1'-binaphthalene-2,2'-diyl diethyl bis(carbonate), C26H22O6, showed that one of the terminal carbon-carbon bond lengths is very short [Csp(3)-Csp(3) = 1.327 (6) angstrom]. The reason for such a short bond length has been analysed by collecting data sets on the same crystal at 393, 150 and 90 K. The values of the corrected bond lengths clearly suggest that the shortening is mainly due to positional disorder at two sites, with minor perturbations arising as a result of thermal vibrations. The positional disorder has been resolved in the analysis of the 90 K data following the changes in the unit-cell parameters for the data sets at 150 and 90 K, which appear to be an artifact of a near centre of symmetry relationship between the two independent molecules in the space group P (1) over bar at these temperatures. Indeed, the unit cell at low temperature (150 and 90 K) is a supercell of the room-temperature unit cell.