A Robust Structural PGN Model for Control of Cell-Cycle Progression Stabilized by Negative Feedbacks

The cell division cycle comprises a sequence of phenomena controlled by a stable and robust genetic network. We applied a probabilistic genetic network (PGN) to construct a hypothetical model with a dynamical behavior displaying the degree of robustness typical of the biological cell cycle. The structure of our PGN model was inspired in well-established biological facts such as the existence of integrator subsystems, negative and positive feedback loops, and redundant signaling pathways. Our model represents genes interactions as stochastic processes and presents strong robustness in the presence of moderate noise and parameters fluctuations. A recently published deterministic yeast cell-cycle model does not perform as well as our PGN model, even upon moderate noise conditions. In addition, self stimulatory mechanisms can give our PGN model the possibility of having a pacemaker activity similar to the observed in the oscillatory embryonic cell cycle.

Treatment of rice straw hemicellulosic hydrolysates with advanced oxidative processes: a new and promising detoxification method to improve the bioconversion process

Abstract Background The use of lignocellulosic constituents in biotechnological processes requires a selective separation of the main fractions (cellulose, hemicellulose and lignin). During diluted acid hydrolysis for hemicellulose extraction, several toxic compounds are formed by the degradation of sugars and lignin, which have ability to inhibit microbial metabolism. Thus, the use of a detoxification step represents an important aspect to be considered for the improvement of fermentation processes from hydrolysates. In this paper, we evaluated the application of Advanced Oxidative Processes (AOPs) for the detoxification of rice straw hemicellulosic hydrolysate with the goal of improving ethanol bioproduction by Pichia stipitis yeast. Aiming to reduce the toxicity of the hemicellulosic hydrolysate, different treatment conditions were analyzed. The treatments were carried out according to a Taguchi L16 orthogonal array to evaluate the influence of Fe+2, H2O2, UV, O3 and pH on the concentration of aromatic compounds and the fermentative process. Results The results showed that the AOPs were able to remove aromatic compounds (furan and phenolic compounds derived from lignin) without affecting the sugar concentration in the hydrolysate. Ozonation in alkaline medium (pH 8) in the presence of H2O2 (treatment A3) or UV radiation (treatment A5) were the most effective for hydrolysate detoxification and had a positive effect on increasing the yeast fermentability of rice straw hemicellulose hydrolysate. Under these conditions, the higher removal of total phenols (above 40%), low molecular weight phenolic compounds (above 95%) and furans (above 52%) were observed. In addition, the ethanol volumetric productivity by P. stipitis was increased in approximately twice in relation the untreated hydrolysate. Conclusion These results demonstrate that AOPs are a promising methods to reduce toxicity and improve the fermentability of lignocellulosic hydrolysates.

Consequences of plant population size for pollinator visitation and plant reproductive success

Habitat loss and fragmentation have a prominent role in determining the size of plant populations, and can affect plant-pollinator interactions. It is hypothesized that in small plant populations the ability to set seeds can be reduced due to limited pollination services, since individuals in small populations can receive less quantity or quality of visits. In this study, I investigated the effect of population size on plant reproductive success and insect visitation in 8 populations of two common species in the island of Lesvos, Greece (Mediterranean Sea), Echium plantagineum and Ballota acetabulosa, and of a rare perennial shrub endemic to north-central Italy, Ononis masquillierii. All the three species depended on insect pollinators for sexual reproduction. For each species, pollen limitation was present in all or nearly all populations, but the relationship between pollen limitation and population size was only present in Ononis masquillierii. However, in Echium plantagineum, significant relationships between both open-pollinated and handcrossed-pollinated seed sets and population size were found, being small populations comparatively less productive than large ones. Additionally, for this species, livestock grazing intensity was greater for small populations and for sparse patches, and had a negative influence on productivity of the remnant plants. Both Echium plantagineum and Ballota acetabulosa attracted a great number of insects, representing a wide spectrum of pollinators, thereby can be considered as generalist species. For Ballota acetabulosa, the most important pollinators were megachilid female bees, and insect diversity didn’t decrease with decreasing plant population size. By contrast, Ononis masquillierii plants generally received few visits, with flowers specialized on small bees (Lasioglossum spp.), representing the most important insect guild. In Echium plantagineum and Ballota acetabulosa, plants in small and large populations received the same amount of visits per flower, and no differences in the number of intraplant visited flowers were detected. On the contrary, large Ononis populations supported higher amounts of pollinators than small ones. At patch level, high Echium flower density was associated with more and higher quality pollinators. My results indicate that small populations were not subject to reduced pollination services than large ones in Echium plantagineum and Ballota acetabulosa, and suggest that grazing and resource limitation could have a major impact on population fitness in Echium plantagineum. The absence of any size effects in these two species can be explained in the light of their high local abundance, wide habitat specificity, and ability to compete with other co-flowering species for pollinators. By contrast, size represents a key characteristic for both pollination and reproduction in Ononis masquillierii populations, as an increase in size could mitigate the negative effects coming from the disadvantageous reproductive traits of the species. Finally, the widespread occurrence of pollen limitation in the three species may be the result of 1) an ongoing weakening or disruption of plantpollinator interactions derived from ecological perturbations, 2) an adaptive equilibrium in response to stochastic processes, and 3) the presence of unfavourable reproductive traits (for Ononis masquillierii).

Statistical properties of Radio Halos and the re-acceleration model

Galaxy clusters occupy a special position in the cosmic hierarchy as they are the largest bound structures in the Universe. There is now general agreement on a hierarchical picture for the formation of cosmic structures, in which galaxy clusters are supposed to form by accretion of matter and merging between smaller units. During merger events, shocks are driven by the gravity of the dark matter in the diffuse barionic component, which is heated up to the observed temperature. Radio and hard-X ray observations have discovered non-thermal components mixed with the thermal Intra Cluster Medium (ICM) and this is of great importance as it calls for a “revision” of the physics of the ICM. The bulk of present information comes from the radio observations which discovered an increasing number of Mpcsized emissions from the ICM, Radio Halos (at the cluster center) and Radio Relics (at the cluster periphery). These sources are due to synchrotron emission from ultra relativistic electrons diffusing through µG turbulent magnetic fields. Radio Halos are the most spectacular evidence of non-thermal components in the ICM and understanding the origin and evolution of these sources represents one of the most challenging goal of the theory of the ICM. Cluster mergers are the most energetic events in the Universe and a fraction of the energy dissipated during these mergers could be channelled into the amplification of the magnetic fields and into the acceleration of high energy particles via shocks and turbulence driven by these mergers. Present observations of Radio Halos (and possibly of hard X-rays) can be best interpreted in terms of the reacceleration scenario in which MHD turbulence injected during these cluster mergers re-accelerates high energy particles in the ICM. The physics involved in this scenario is very complex and model details are difficult to test, however this model clearly predicts some simple properties of Radio Halos (and resulting IC emission in the hard X-ray band) which are almost independent of the details of the adopted physics. In particular in the re-acceleration scenario MHD turbulence is injected and dissipated during cluster mergers and thus Radio Halos (and also the resulting hard X-ray IC emission) should be transient phenomena (with a typical lifetime <» 1 Gyr) associated with dynamically disturbed clusters. The physics of the re-acceleration scenario should produce an unavoidable cut-off in the spectrum of the re-accelerated electrons, which is due to the balance between turbulent acceleration and radiative losses. The energy at which this cut-off occurs, and thus the maximum frequency at which synchrotron radiation is produced, depends essentially on the efficiency of the acceleration mechanism so that observations at high frequencies are expected to catch only the most efficient phenomena while, in principle, low frequency radio surveys may found these phenomena much common in the Universe. These basic properties should leave an important imprint in the statistical properties of Radio Halos (and of non-thermal phenomena in general) which, however, have not been addressed yet by present modellings. The main focus of this PhD thesis is to calculate, for the first time, the expected statistics of Radio Halos in the context of the re-acceleration scenario. In particular, we shall address the following main questions: • Is it possible to model “self-consistently” the evolution of these sources together with that of the parent clusters? • How the occurrence of Radio Halos is expected to change with cluster mass and to evolve with redshift? How the efficiency to catch Radio Halos in galaxy clusters changes with the observing radio frequency? • How many Radio Halos are expected to form in the Universe? At which redshift is expected the bulk of these sources? • Is it possible to reproduce in the re-acceleration scenario the observed occurrence and number of Radio Halos in the Universe and the observed correlations between thermal and non-thermal properties of galaxy clusters? • Is it possible to constrain the magnetic field intensity and profile in galaxy clusters and the energetic of turbulence in the ICM from the comparison between model expectations and observations? Several astrophysical ingredients are necessary to model the evolution and statistical properties of Radio Halos in the context of re-acceleration model and to address the points given above. For these reason we deserve some space in this PhD thesis to review the important aspects of the physics of the ICM which are of interest to catch our goals. In Chapt. 1 we discuss the physics of galaxy clusters, and in particular, the clusters formation process; in Chapt. 2 we review the main observational properties of non-thermal components in the ICM; and in Chapt. 3 we focus on the physics of magnetic field and of particle acceleration in galaxy clusters. As a relevant application, the theory of Alfv´enic particle acceleration is applied in Chapt. 4 where we report the most important results from calculations we have done in the framework of the re-acceleration scenario. In this Chapter we show that a fraction of the energy of fluid turbulence driven in the ICM by the cluster mergers can be channelled into the injection of Alfv´en waves at small scales and that these waves can efficiently re-accelerate particles and trigger Radio Halos and hard X-ray emission. The main part of this PhD work, the calculation of the statistical properties of Radio Halos and non-thermal phenomena as expected in the context of the re-acceleration model and their comparison with observations, is presented in Chapts.5, 6, 7 and 8. In Chapt.5 we present a first approach to semi-analytical calculations of statistical properties of giant Radio Halos. The main goal of this Chapter is to model cluster formation, the injection of turbulence in the ICM and the resulting particle acceleration process. We adopt the semi–analytic extended Press & Schechter (PS) theory to follow the formation of a large synthetic population of galaxy clusters and assume that during a merger a fraction of the PdV work done by the infalling subclusters in passing through the most massive one is injected in the form of magnetosonic waves. Then the processes of stochastic acceleration of the relativistic electrons by these waves and the properties of the ensuing synchrotron (Radio Halos) and inverse Compton (IC, hard X-ray) emission of merging clusters are computed under the assumption of a constant rms average magnetic field strength in emitting volume. The main finding of these calculations is that giant Radio Halos are naturally expected only in the more massive clusters, and that the expected fraction of clusters with Radio Halos is consistent with the observed one. In Chapt. 6 we extend the previous calculations by including a scaling of the magnetic field strength with cluster mass. The inclusion of this scaling allows us to derive the expected correlations between the synchrotron radio power of Radio Halos and the X-ray properties (T, LX) and mass of the hosting clusters. For the first time, we show that these correlations, calculated in the context of the re-acceleration model, are consistent with the observed ones for typical µG strengths of the average B intensity in massive clusters. The calculations presented in this Chapter allow us to derive the evolution of the probability to form Radio Halos as a function of the cluster mass and redshift. The most relevant finding presented in this Chapter is that the luminosity functions of giant Radio Halos at 1.4 GHz are expected to peak around a radio power » 1024 W/Hz and to flatten (or cut-off) at lower radio powers because of the decrease of the electron re-acceleration efficiency in smaller galaxy clusters. In Chapt. 6 we also derive the expected number counts of Radio Halos and compare them with available observations: we claim that » 100 Radio Halos in the Universe can be observed at 1.4 GHz with deep surveys, while more than 1000 Radio Halos are expected to be discovered in the next future by LOFAR at 150 MHz. This is the first (and so far unique) model expectation for the number counts of Radio Halos at lower frequency and allows to design future radio surveys. Based on the results of Chapt. 6, in Chapt.7 we present a work in progress on a “revision” of the occurrence of Radio Halos. We combine past results from the NVSS radio survey (z » 0.05 − 0.2) with our ongoing GMRT Radio Halos Pointed Observations of 50 X-ray luminous galaxy clusters (at z » 0.2−0.4) and discuss the possibility to test our model expectations with the number counts of Radio Halos at z » 0.05 − 0.4. The most relevant limitation in the calculations presented in Chapt. 5 and 6 is the assumption of an “averaged” size of Radio Halos independently of their radio luminosity and of the mass of the parent clusters. This assumption cannot be released in the context of the PS formalism used to describe the formation process of clusters, while a more detailed analysis of the physics of cluster mergers and of the injection process of turbulence in the ICM would require an approach based on numerical (possible MHD) simulations of a very large volume of the Universe which is however well beyond the aim of this PhD thesis. On the other hand, in Chapt.8 we report our discovery of novel correlations between the size (RH) of Radio Halos and their radio power and between RH and the cluster mass within the Radio Halo region, MH. In particular this last “geometrical” MH − RH correlation allows us to “observationally” overcome the limitation of the “average” size of Radio Halos. Thus in this Chapter, by making use of this “geometrical” correlation and of a simplified form of the re-acceleration model based on the results of Chapt. 5 and 6 we are able to discuss expected correlations between the synchrotron power and the thermal cluster quantities relative to the radio emitting region. This is a new powerful tool of investigation and we show that all the observed correlations (PR − RH, PR − MH, PR − T, PR − LX, . . . ) now become well understood in the context of the re-acceleration model. In addition, we find that observationally the size of Radio Halos scales non-linearly with the virial radius of the parent cluster, and this immediately means that the fraction of the cluster volume which is radio emitting increases with cluster mass and thus that the non-thermal component in clusters is not self-similar.

Probabilistic approaches in civil engineering: generation of random fields and structural identification with genetic algorithms

The inherent stochastic character of most of the physical quantities involved in engineering models has led to an always increasing interest for probabilistic analysis. Many approaches to stochastic analysis have been proposed. However, it is widely acknowledged that the only universal method available to solve accurately any kind of stochastic mechanics problem is Monte Carlo Simulation. One of the key parts in the implementation of this technique is the accurate and efficient generation of samples of the random processes and fields involved in the problem at hand. In the present thesis an original method for the simulation of homogeneous, multi-dimensional, multi-variate, non-Gaussian random fields is proposed. The algorithm has proved to be very accurate in matching both the target spectrum and the marginal probability. The computational efficiency and robustness are very good too, even when dealing with strongly non-Gaussian distributions. What is more, the resulting samples posses all the relevant, welldefined and desired properties of “translation fields”, including crossing rates and distributions of extremes. The topic of the second part of the thesis lies in the field of non-destructive parametric structural identification. Its objective is to evaluate the mechanical characteristics of constituent bars in existing truss structures, using static loads and strain measurements. In the cases of missing data and of damages that interest only a small portion of the bar, Genetic Algorithm have proved to be an effective tool to solve the problem.

Study of the circulation processes in the northern Adriatic sea - coastal area and Venice lagoon inlets

This work is a detailed study of hydrodynamic processes in a defined area, the littoral in front of the Venice Lagoon and its inlets, which are complex morphological areas of interconnection. A finite element hydrodynamic model of the Venice Lagoon and the Adriatic Sea has been developed in order to study the coastal current patterns and the exchanges at the inlets of the Venice Lagoon. This is the first work in this area that tries to model the interaction dynamics, running together a model for the lagoon and the Adriatic Sea. First the barotropic processes near the inlets of the Venice Lagoon have been studied. Data from more than ten tide gauges displaced in the Adriatic Sea have been used in the calibration of the simulated water levels. To validate the model results, empirical flux data measured by ADCP probes installed inside the inlets of Lido and Malamocco have been used and the exchanges through the three inlets of the Venice Lagoon have been analyzed. The comparison between modelled and measured fluxes at the inlets outlined the efficiency of the model to reproduce both tide and wind induced water exchanges between the sea and the lagoon. As a second step, also small scale processes around the inlets that connect the Venice lagoon with the Northern Adriatic Sea have been investigated by means of 3D simulations. Maps of vorticity have been produced, considering the influence of tidal flows and wind stress in the area. A sensitivity analysis has been carried out to define the importance of the advection and of the baroclinic pressure gradients in the development of vortical processes seen along the littoral close to the inlets. Finally a comparison with real data measurements, surface velocity data from HF Radar near the Venice inlets, has been performed, which allows for a better understanding of the processes and their seasonal dynamics. The results outline the predominance of wind and tidal forcing in the coastal area. Wind forcing acts mainly on the mean coastal current inducing its detachment offshore during Sirocco events and an increase of littoral currents during Bora events. The Bora action is more homogeneous on the whole coastal area whereas the Sirocco strengthens its impact in the South, near Chioggia inlet. Tidal forcing at the inlets is mainly barotropic. The sensitivity analysis shows how advection is the main physical process responsible for the persistent vortical structures present along the littoral between the Venice Lagoon inlets. The comparison with measurements from HF Radar not only permitted a validation the model results, but also a description of different patterns in specific periods of the year. The success of the 2D and the 3D simulations on the reproduction both of the SSE, inside and outside the Venice Lagoon, of the tidal flow, through the lagoon inlets, and of the small scale phenomena, occurring along the littoral, indicates that the finite element approach is the most suitable tool for the investigation of coastal processes. For the first time, as shown by the flux modeling, the physical processes that drive the interaction between the two basins were reproduced.

Stochastic models and dynamic measures for the characterization of bistable circuits

During the last few years, a great deal of interest has risen concerning the applications of stochastic methods to several biochemical and biological phenomena. Phenomena like gene expression, cellular memory, bet-hedging strategy in bacterial growth and many others, cannot be described by continuous stochastic models due to their intrinsic discreteness and randomness. In this thesis I have used the Chemical Master Equation (CME) technique to modelize some feedback cycles and analyzing their properties, including experimental data. In the first part of this work, the effect of stochastic stability is discussed on a toy model of the genetic switch that triggers the cellular division, which malfunctioning is known to be one of the hallmarks of cancer. The second system I have worked on is the so-called futile cycle, a closed cycle of two enzymatic reactions that adds and removes a chemical compound, called phosphate group, to a specific substrate. I have thus investigated how adding noise to the enzyme (that is usually in the order of few hundred molecules) modifies the probability of observing a specific number of phosphorylated substrate molecules, and confirmed theoretical predictions with numerical simulations. In the third part the results of the study of a chain of multiple phosphorylation-dephosphorylation cycles will be presented. We will discuss an approximation method for the exact solution in the bidimensional case and the relationship that this method has with the thermodynamic properties of the system, which is an open system far from equilibrium.In the last section the agreement between the theoretical prediction of the total protein quantity in a mouse cells population and the observed quantity will be shown, measured via fluorescence microscopy.

Statistical problems related to excitation threshold and reset value of membrane potentials

Die vorliegende Arbeit ist motiviert durch biologische Fragestellungen bezüglich des Verhaltens von Membranpotentialen in Neuronen. Ein vielfach betrachtetes Modell für spikende Neuronen ist das Folgende. Zwischen den Spikes verhält sich das Membranpotential wie ein Diffusionsprozess X der durch die SDGL dX_t= beta(X_t) dt+ sigma(X_t) dB_t gegeben ist, wobei (B_t) eine Standard-Brown'sche Bewegung bezeichnet. Spikes erklärt man wie folgt. Sobald das Potential X eine gewisse Exzitationsschwelle S überschreitet entsteht ein Spike. Danach wird das Potential wieder auf einen bestimmten Wert x_0 zurückgesetzt. In Anwendungen ist es manchmal möglich, einen Diffusionsprozess X zwischen den Spikes zu beobachten und die Koeffizienten der SDGL beta() und sigma() zu schätzen. Dennoch ist es nötig, die Schwellen x_0 und S zu bestimmen um das Modell festzulegen. Eine Möglichkeit, dieses Problem anzugehen, ist x_0 und S als Parameter eines statistischen Modells aufzufassen und diese zu schätzen. In der vorliegenden Arbeit werden vier verschiedene Fälle diskutiert, in denen wir jeweils annehmen, dass das Membranpotential X zwischen den Spikes eine Brown'sche Bewegung mit Drift, eine geometrische Brown'sche Bewegung, ein Ornstein-Uhlenbeck Prozess oder ein Cox-Ingersoll-Ross Prozess ist. Darüber hinaus beobachten wir die Zeiten zwischen aufeinander folgenden Spikes, die wir als iid Treffzeiten der Schwelle S von X gestartet in x_0 auffassen. Die ersten beiden Fälle ähneln sich sehr und man kann jeweils den Maximum-Likelihood-Schätzer explizit angeben. Darüber hinaus wird, unter Verwendung der LAN-Theorie, die Optimalität dieser Schätzer gezeigt. In den Fällen OU- und CIR-Prozess wählen wir eine Minimum-Distanz-Methode, die auf dem Vergleich von empirischer und wahrer Laplace-Transformation bezüglich einer Hilbertraumnorm beruht. Wir werden beweisen, dass alle Schätzer stark konsistent und asymptotisch normalverteilt sind. Im letzten Kapitel werden wir die Effizienz der Minimum-Distanz-Schätzer anhand simulierter Daten überprüfen. Ferner, werden Anwendungen auf reale Datensätze und deren Resultate ausführlich diskutiert.

Dynamics of epithelial monolayers assessed by acoustic and impedimetric whole cell biosensors

Mit der Zielsetzung der vorliegenden Arbeit wurde die detailierten Analyse von Migrationsdynamiken epithelilaler Monolayer anhand zweier neuartiger in vitro Biosensoren verfolgt, der elektrischen Zell-Substrat Impedanz Spektroskopie (electrical cell-substrate impedance sensing, ECIS) sowie der Quarz Kristall Mikrowaage (quartz crystal microbalance, QCM). Beide Methoden erwiesen sich als sensitiv gegenüber der Zellmotilität und der Nanozytotoxizität.rnInnerhalb des ersten Projektes wurde ein Fingerprinting von Krebszellen anhand ihrer Motilitätsdynamiken und der daraus generierten elektrischen oder akkustischen Fluktuationen auf ECIS oder QCM Basis vorgenommen; diese Echtzeitsensoren wurdene mit Hilfe klassicher in vitro Boyden-Kammer Migrations- und Invasions-assays validiert. Fluktuationssignaturen, also Langzeitkorrelationen oder fraktale Selbstähnlichkeit aufgrund der kollektiven Zellbewegung, wurden über Varianz-, Fourier- sowie trendbereinigende Fluktuationsanalyse quantifiziert. Stochastische Langzeitgedächtnisphänomene erwiesen sich als maßgebliche Beiträge zur Antwort adhärenter Zellen auf den QCM und ECIS-Sensoren. Des weiteren wurde der Einfluss niedermolekularer Toxine auf die Zytoslelettdynamiken verfolgt: die Auswirkungen von Cytochalasin D, Phalloidin und Blebbistatin sowie Taxol, Nocodazol und Colchicin wurden dabei über die QCM und ECIS Fluktuationsanalyse erfasst.rnIn einem zweiten Projektschwerpunkt wurden Adhäsionsprozesse sowie Zell-Zell und Zell-Substrat Degradationsprozesse bei Nanopartikelgabe charackterisiert, um ein Maß für Nanozytotoxizität in Abhangigkeit der Form, Funktionalisierung Stabilität oder Ladung der Partikel zu erhalten.rnAls Schlussfolgerung ist zu nennen, dass die neuartigen Echtzeit-Biosensoren QCM und ECIS eine hohe Zellspezifität besitzen, auf Zytoskelettdynamiken reagieren sowie als sensitive Detektoren für die Zellvitalität fungieren können.

The behaviour of osteocyte dendritic processes in bone under cyclic load (il comportamento dei processi dendritici degli osteociti nel tessuto osseo sottoposti a carico ciclico)

Bone is continually being removed and replaced through the actions of basic multicellular units (BMU). This constant upkeep is necessary to remove microdamage formed naturally due to fatigue and thus maintain the integrity of the bone. The repair process in bone is targeted, meaning that a BMU travels directly to the site of damage and repairs it. It is still unclear how targeted remodelling is stimulated and directed but it is highly likely that osteocytes play a role. A number of theories have been advanced to explain the microcrack osteocyte interaction but no complete mechanism has been demonstrated. Osteocytes are connected to each other by dendritic processes. The “scissors model" proposed that the rupture of these processes where they cross microcracks signals the degree of damage and the urgency of the necessary repair. In its original form it was proposed that under applied compressive loading, microcrack faces will be pressed together and undergo relative shear movement. If this movement is greater than the width of an osteocyte process, then the process will be cut in a “scissors like" motion, releasing RANKL, a cytokine known to be essential in the formation of osteoclasts from pre-osteoclasts. The main aim of this thesis was to investigate this theoretical model with a specific focus on microscopy and finite element modelling. Previous studies had proved that cyclic stress was necessary for osteocyte process rupture to occur. This was a divergence from the original “scissors model" which had proposed that the cutting of cell material occurred in one single action. The present thesis is the first study to show fatigue failure in cellular processes spanning naturally occurring cracks and it's the first study to estimate the cyclic strain range and relate it to the number of cycles to failure, for any type of cell. Rupture due to shear movement was ruled out as microcrack closing never occurred, as a result of plastic deformation of the bone. Fatigue failure was found to occur due to cyclic tensile stress in the locality of the damage. The strain range necessary for osteocyte process rupture was quantified. It was found that the lower the process strain range the greater the number of cycles to cell process failure. FEM modelling allowed to predict stress in the vicinity of an osteocyte process and to analyse its interaction with the bone surrounding it: simulations revealed evident creep effects in bone during cyclic loading. This thesis confirms and dismisses aspects of the “scissors model". The observations support the model as a viable mechanism of microcrack detection by the osteocyte network, albeit in a slightly modified form where cyclic loading is necessary and the method of rupture is fatigue failure due to cyclic tensile motion. An in depth study was performed focusing on microscopy analysis of naturally occurring cracks in bone and FEM simulation analysis of an osteocyte process spanning a microcrack in bone under cyclic load.

Impact of anthropogenic stressors on the population genetics of the common bioindicator Dreissena polymorpha (Pallas, 1771)

Toxicant inputs from agriculture, industry and human settlements have been shown to severely affect freshwater ecosystems. Pollution can lead to changes in population genetic patterns through various genetic and stochastic processes. In my thesis, I investigated the impact of anthropogenic stressors on the population genetics of the zebra mussel Dreissena polymorpha. In order to analyze the genetics of zebra mussel populations, I isolated five new highly polymorphic microsatellite loci. Out of those and other already existing microsatellite markers for this species, I established a robust marker set of six microsatellite loci for D. polymorpha. rnMonitoring the biogeographical background is an important requirement when integrating population genetic measures into ecotoxicological studies. I analyzed the biogeographical background of eleven populations in a section of the River Danube (in Hungary and Croatia) and some of its tributaries, and another population in the River Rhine as genetic outgroup. Moreover, I measured abiotic water parameters at the sampling sites and analyzed if they were correlated with the genetic parameters of the populations. The genetic differentiation was basically consistent with the overall biogeographical history of the populations in the study region. However, the genetic diversity of the populations was not influenced by the geographical distance between the populations, but by the environmental factors oxygen and temperature and also by other unidentified factors. I found strong evidence that genetic adaptation of zebra mussel populations to local habitat conditions had influenced the genetic constitution of the populations. Moreover, by establishing the biogeographical baseline of molecular variance in the study area, I laid the foundation for interpreting population genetic results in ecotoxicological experiments in this region.rnIn a cooperation project with the Department of Zoology of the University of Zagreb, I elaborated an integrated approach in biomonitoring with D. polymorpha by combining the analysis techniques of microsatellite analysis, Comet assay and micronucleus test (MNT). This approach was applied in a case study on freshwater contamination by an effluent of a wastewater treatment plant (WWTP) in the River Drava (Croatia) and a complementary laboratory experiment. I assessed and compared the genetic status of two zebra mussel populations from a contaminated and a reference site. Microsatellite analysis suggested that the contaminated population had undergone a genetic bottleneck, caused by random genetic drift and selection, whereas a bottleneck was not detected in the reference population. The Comet assay did not indicate any difference in DNA damage between the two populations, but MNT revealed that the contaminated population had an increased percentage of micronuclei in hemocytes in comparison to the reference population. The laboratory experiment with mussels exposed to municipal wastewater revealed that mussels from the contaminated site had a lower percentage of tail DNA and a higher percentage of micronuclei than the reference population. These differences between populations were probably caused by an overall decreased fitness of mussels from the contaminated site due to genetic drift and by an enhanced DNA repair mechanism due to adaptation to pollution in the source habitat. Overall, the combination of the three biomarkers provided sufficient information on the impact of both treated and non-treated municipal wastewater on the genetics of zebra mussels at different levels of biological organization.rnIn my thesis, I could show that the newly established marker set of six microsatellite loci provided reliable and informative data for population genetic analyses of D. polymorpha. The adaptation of the analyzed zebra mussel populations to the local conditions of their habitat had a strong influence on their genetic constitution. We found evidence that the different genetic constitutions of two populations had influenced the outcome of our ecotoxicological experiment. Overall, the integrated approach in biomonitoring gave comprehensive information about the impact of both treated and non-treated municipal wastewater on the genetics of zebra mussels at different levels of biological organization and was well practicable in a first case study.

Active processes, morphology, and dynamics of icy debris fans: Landform evolution along rapidly degrading escarpments in alpine regions undergoing recent deglaciation

In the past few decades the impacts of climate warming have been significant in alpine glaciated regions. Many valley glaciers formerly linked as distributary glaciers to high-level icecaps have decoupled at their icefalls, exposing major escarpments and generating a suite of dynamic landforrns dominated by mass wasting. Ice-dominated landforms, here termed icy debris fans, develop rapidly by ice avalanching, rockfall, and icy debris flow. Field-based reconnaissance studies at two alpine settings, the Wrangell Mountains of Alaska and the Southern Alps of New Zealand, provide a preliminary morphogenetic model of spatial and temporal evolution of icy debris fans in a range of alpine settings. The influence of these processes on landform evolution is largely unrecognized in the literature dealing with post-glacial landform adjustment known as the paraglacial. A better understanding of these dynamic processes will be increasingly important because of the extreme geohazards characterizing these areas. Our field studies show that after glacier decoupling, icy debris fans begin to form along the base of bedrock escarpments at the mouths of catchments and prograde over valley glaciers. The presence of a distinct catchment, apex, and fan morphology distinguishes these landforms from other landforms common in periglacial hillslope settings receiving abundant clastic debris and ice. Ice avalanching is the most abundant process involved in icy debris fan formation. Fans developed below weakly incised catchments are dominated by ice avalanching and are composed primarily of ice with minor lithic detritus. Typically, avalanches fall into the fan catchments where sediments transform into grainflows that flow onto the fans. Once on the fans, avalanche deposits ablate rapidly, flattening and concentrating lithic fragments at the surface. Icy debris fans may become thick enough to become glaciers with splay crevasse systems. Fans developed below larger, more complex catchments are composed of higher proportions of lithic detritus resulting from temporary storage of ice and lithic detritus deposits within the catchment. Episodic outbursts of meltwater from the icecap may mix with the stored sediments and mobilize icy debris flows (mixture of ice and lithic clasts) onto the fans. Our observations indicate that the entire evolutionary cycle of icy debris fans probably occurs during an early paraglacial interval (i.e., decades to 100 years). Observations comparing avalanche frequency, volume, and fan morphologic evolution at the Alaska site between 2006 and 2010 illustrate complex response between icy debris fans even within the same cirque - where one fan may be growing while others are downwasting because of differences in ice supply controlled by their respective catchments and icecap contributions. As ice supply from the icecap diminishes through time, icy debris fans rapidly downwaste and eventually evolve into talus cones that receive occasional but ephemeral ice avalanches.

Processes Regulating the Initiation and Postejaculatory Resumption of Copulatory Behavior in Male Hamsters

Studies using factor analysis have helped describe the organization of copulatory behavior in male rodents. However, the focus of these studies on a few traditional measures may have limited their results. To test this possibility, 74 sexually-experienced male hamsters were observed as they copulated with stimulus females. The measures collected exceeded the conventional ones in number, variety and independence. The factor analysis of these data revealed a structure with seven factors collectively accounting for 80% of the variance. Most resembled the factors in previous reports, reinforcing the contributions that the processes suggested by these factors make to the organization,of male behavior. But several other factors were more novel, possibly reflecting the use of measures that were novel or revised for greater independence. The most interesting of these were two factors focusing on early steps in the progression leading to ejaculation. Importantly, both incorporated measures from each of the three copulatory series that were observed. Past work suggests that independent processes control the times required to initiate copulation and later resume it after an ejaculation. In contrast, these results suggest the existence of two processes, each of which contributes to both the initiation and reinitiation of copulation. (C) 2014 Elsevier B.V. All rights reserved.

Generalizations of Current Status Data With Applications

In estimation of a survival function, current status data arises when the only information available on individuals is their survival status at a single monitoring time. Here we briefly review extensions of this form of data structure in two directions: (i) doubly censored current status data, where there is incomplete information on the origin of the failure time random variable, and (ii) current status information on more complicated stochastic processes. Simple examples of these data forms are presented for motivation.

IDENTIFICATION OF GENES CONTROLLING BIOLOGICAL PROCESSES AND PATHWAYS THROUGH STATISTICAL ANALYSIS AND NETWORK RECONSTRUCTION

The developmental processes and functions of an organism are controlled by the genes and the proteins that are derived from these genes. The identification of key genes and the reconstruction of gene networks can provide a model to help us understand the regulatory mechanisms for the initiation and progression of biological processes or functional abnormalities (e.g. diseases) in living organisms. In this dissertation, I have developed statistical methods to identify the genes and transcription factors (TFs) involved in biological processes, constructed their regulatory networks, and also evaluated some existing association methods to find robust methods for coexpression analyses. Two kinds of data sets were used for this work: genotype data and gene expression microarray data. On the basis of these data sets, this dissertation has two major parts, together forming six chapters. The first part deals with developing association methods for rare variants using genotype data (chapter 4 and 5). The second part deals with developing and/or evaluating statistical methods to identify genes and TFs involved in biological processes, and construction of their regulatory networks using gene expression data (chapter 2, 3, and 6). For the first part, I have developed two methods to find the groupwise association of rare variants with given diseases or traits. The first method is based on kernel machine learning and can be applied to both quantitative as well as qualitative traits. Simulation results showed that the proposed method has improved power over the existing weighted sum method (WS) in most settings. The second method uses multiple phenotypes to select a few top significant genes. It then finds the association of each gene with each phenotype while controlling the population stratification by adjusting the data for ancestry using principal components. This method was applied to GAW 17 data and was able to find several disease risk genes. For the second part, I have worked on three problems. First problem involved evaluation of eight gene association methods. A very comprehensive comparison of these methods with further analysis clearly demonstrates the distinct and common performance of these eight gene association methods. For the second problem, an algorithm named the bottom-up graphical Gaussian model was developed to identify the TFs that regulate pathway genes and reconstruct their hierarchical regulatory networks. This algorithm has produced very significant results and it is the first report to produce such hierarchical networks for these pathways. The third problem dealt with developing another algorithm called the top-down graphical Gaussian model that identifies the network governed by a specific TF. The network produced by the algorithm is proven to be of very high accuracy.