Quasi-integrals and the origin of bursting pH oscillations in an enzyme model system

Bursting Oscillation, Mixed-Mode oscillation, Slow Manifold, Quasi-Integral, slow-fast analysis, QSSA

Model-based diagnosis system: application to constructed wetlands

Report for the scientific sojourn carried out at the Model-based Systems and Qualitative Reasoning Group (Technical University of Munich), from September until December 2005. Constructed wetlands (CWs), or modified natural wetlands, are used all over the world as wastewater treatment systems for small communities because they can provide high treatment efficiency with low energy consumption and low construction, operation and maintenance costs. Their treatment process is very complex because it includes physical, chemical and biological mechanisms like microorganism oxidation, microorganism reduction, filtration, sedimentation and chemical precipitation. Besides, these processes can be influenced by different factors. In order to guarantee the performance of CWs, an operation and maintenance program must be defined for each Wastewater Treatment Plant (WWTP). The main objective of this project is to provide a computer support to the definition of the most appropriate operation and maintenance protocols to guarantee the correct performance of CWs. To reach them, the definition of models which represent the knowledge about CW has been proposed: components involved in the sanitation process, relation among these units and processes to remove pollutants. Horizontal Subsurface Flow CWs are chosen as a case study and the filtration process is selected as first modelling-process application. However, the goal is to represent the process knowledge in such a way that it can be reused for other types of WWTP.

A vertebrate reproductive system involving three ploidy levels : hybrid origin of triploids in a contact zone of diploid and tetraploid green toads (Bufo viridis subgroup)

The rise and consequences of polyploidy in vertebrates, whose origin was associated with genome duplications, may be best studied in natural diploid and polyploid populations. In a diploid/tetraploid (2n/4n) geographic contact zone of Palearctic green toads in northern Kyrgyzstan, we examine 4ns and triploids (3n) of unknown genetic composition and origins. Using mitochondrial and nuclear sequence, and nuclear microsatellite markers in 84 individuals, we show that 4n (Bufo pewzowi) are allopolyploids, with a geographically proximate 2n species (B. turanensis) being their maternal ancestor and their paternal ancestor as yet unidentified. Local 3n forms arise through hybridization. Adult 3n mature males (B. turanensis mtDNA) have 2n mothers and 4n fathers, but seem distinguishable by nuclear profiles from partly aneuploid 3n tadpoles (with B. pewzowi mtDNA). These observations suggest multiple pathways to the formation of triploids in the contact zone, involving both reciprocal origins. To explain the phenomena in the system, we favor a hypothesis where 3n males (with B. turanensis mtDNA) backcross with 4n and 2n females. Together with previous studies of a separately evolved, sexually reproducing 3n lineage, these observations reveal complex reproductive interactions among toads of different ploidy levels and multiple pathways to the evolution of polyploid lineages.

A peripheral dose calculation model for estimating second cancer risk after breast radiotherapy

AbstractBreast cancer is one of the most common cancers affecting one in eight women during their lives. Survival rates have increased steadily thanks to early diagnosis with mammography screening and more efficient treatment strategies. Post-operative radiation therapy is a standard of care in the management of breast cancer and has been shown to reduce efficiently both local recurrence rate and breast cancer mortality. Radiation therapy is however associated with some late effects for long-term survivors. Radiation-induced secondary cancer is a relatively rare but severe late effect of radiation therapy. Currently, radiotherapy plans are essentially optimized to maximize tumor control and minimize late deterministic effects (tissue reactions) that are mainly associated with high doses (» 1 Gy). With improved cure rates and new radiation therapy technologies, it is also important to evaluate and minimize secondary cancer risks for different treatment techniques. This is a particularly challenging task due to the large uncertainties in the dose-response relationship.In contrast with late deterministic effects, secondary cancers may be associated with much lower doses and therefore out-of-field doses (also called peripheral doses) that are typically inferior to 1 Gy need to be determined accurately. Out-of-field doses result from patient scatter and head scatter from the treatment unit. These doses are particularly challenging to compute and we characterized it by Monte Carlo (MC) calculation. A detailed MC model of the Siemens Primus linear accelerator has been thoroughly validated with measurements. We investigated the accuracy of such a model for retrospective dosimetry in epidemiological studies on secondary cancers. Considering that patients in such large studies could be treated on a variety of machines, we assessed the uncertainty in reconstructed peripheral dose due to the variability of peripheral dose among various linac geometries. For large open fields (> 10x10 cm2), the uncertainty would be less than 50%, but for small fields and wedged fields the uncertainty in reconstructed dose could rise up to a factor of 10. It was concluded that such a model could be used for conventional treatments using large open fields only.The MC model of the Siemens Primus linac was then used to compare out-of-field doses for different treatment techniques in a female whole-body CT-based phantom. Current techniques such as conformai wedged-based radiotherapy and hybrid IMRT were investigated and compared to older two-dimensional radiotherapy techniques. MC doses were also compared to those of a commercial Treatment Planning System (TPS). While the TPS is routinely used to determine the dose to the contralateral breast and the ipsilateral lung which are mostly out of the treatment fields, we have shown that these doses may be highly inaccurate depending on the treatment technique investigated. MC shows that hybrid IMRT is dosimetrically similar to three-dimensional wedge-based radiotherapy within the field, but offers substantially reduced doses to out-of-field healthy organs.Finally, many different approaches to risk estimations extracted from the literature were applied to the calculated MC dose distribution. Absolute risks varied substantially as did the ratio of risk between two treatment techniques, reflecting the large uncertainties involved with current risk models. Despite all these uncertainties, the hybrid IMRT investigated resulted in systematically lower cancer risks than any of the other treatment techniques. More epidemiological studies with accurate dosimetry are required in the future to construct robust risk models. In the meantime, any treatment strategy that reduces out-of-field doses to healthy organs should be investigated. Electron radiotherapy might offer interesting possibilities with this regard.RésuméLe cancer du sein affecte une femme sur huit au cours de sa vie. Grâce au dépistage précoce et à des thérapies de plus en plus efficaces, le taux de guérison a augmenté au cours du temps. La radiothérapie postopératoire joue un rôle important dans le traitement du cancer du sein en réduisant le taux de récidive et la mortalité. Malheureusement, la radiothérapie peut aussi induire des toxicités tardives chez les patients guéris. En particulier, les cancers secondaires radio-induits sont une complication rare mais sévère de la radiothérapie. En routine clinique, les plans de radiothérapie sont essentiellement optimisées pour un contrôle local le plus élevé possible tout en minimisant les réactions tissulaires tardives qui sont essentiellement associées avec des hautes doses (» 1 Gy). Toutefois, avec l'introduction de différentes nouvelles techniques et avec l'augmentation des taux de survie, il devient impératif d'évaluer et de minimiser les risques de cancer secondaire pour différentes techniques de traitement. Une telle évaluation du risque est une tâche ardue étant donné les nombreuses incertitudes liées à la relation dose-risque.Contrairement aux effets tissulaires, les cancers secondaires peuvent aussi être induits par des basses doses dans des organes qui se trouvent hors des champs d'irradiation. Ces organes reçoivent des doses périphériques typiquement inférieures à 1 Gy qui résultent du diffusé du patient et du diffusé de l'accélérateur. Ces doses sont difficiles à calculer précisément, mais les algorithmes Monte Carlo (MC) permettent de les estimer avec une bonne précision. Un modèle MC détaillé de l'accélérateur Primus de Siemens a été élaboré et validé avec des mesures. La précision de ce modèle a également été déterminée pour la reconstruction de dose en épidémiologie. Si on considère que les patients inclus dans de larges cohortes sont traités sur une variété de machines, l'incertitude dans la reconstruction de dose périphérique a été étudiée en fonction de la variabilité de la dose périphérique pour différents types d'accélérateurs. Pour de grands champs (> 10x10 cm ), l'incertitude est inférieure à 50%, mais pour de petits champs et des champs filtrés, l'incertitude de la dose peut monter jusqu'à un facteur 10. En conclusion, un tel modèle ne peut être utilisé que pour les traitements conventionnels utilisant des grands champs.Le modèle MC de l'accélérateur Primus a été utilisé ensuite pour déterminer la dose périphérique pour différentes techniques dans un fantôme corps entier basé sur des coupes CT d'une patiente. Les techniques actuelles utilisant des champs filtrés ou encore l'IMRT hybride ont été étudiées et comparées par rapport aux techniques plus anciennes. Les doses calculées par MC ont été comparées à celles obtenues d'un logiciel de planification commercial (TPS). Alors que le TPS est utilisé en routine pour déterminer la dose au sein contralatéral et au poumon ipsilatéral qui sont principalement hors des faisceaux, nous avons montré que ces doses peuvent être plus ou moins précises selon la technTque étudiée. Les calculs MC montrent que la technique IMRT est dosimétriquement équivalente à celle basée sur des champs filtrés à l'intérieur des champs de traitement, mais offre une réduction importante de la dose aux organes périphériques.Finalement différents modèles de risque ont été étudiés sur la base des distributions de dose calculées par MC. Les risques absolus et le rapport des risques entre deux techniques de traitement varient grandement, ce qui reflète les grandes incertitudes liées aux différents modèles de risque. Malgré ces incertitudes, on a pu montrer que la technique IMRT offrait une réduction du risque systématique par rapport aux autres techniques. En attendant des données épidémiologiques supplémentaires sur la relation dose-risque, toute technique offrant une réduction des doses périphériques aux organes sains mérite d'être étudiée. La radiothérapie avec des électrons offre à ce titre des possibilités intéressantes.

Modeling of hybrid systems by piecewise decomposition

Piecewise linear models systems arise as mathematical models of systems in many practical applications, often from linearization for nonlinear systems. There are two main approaches of dealing with these systems according to their continuous or discrete-time aspects. We propose an approach which is based on the state transformation, more particularly the partition of the phase portrait in different regions where each subregion is modeled as a two-dimensional linear time invariant system. Then the Takagi-Sugeno model, which is a combination of local model is calculated. The simulation results show that the Alpha partition is well-suited for dealing with such a system

Hybrid consumption paths in the attribute space: A model and application with scanner data

This paper presents and estimates a dynamic choice model in the attribute space considering rational consumers. In light of the evidence of several state-dependence patterns, the standard attribute-based model is extended by considering a general utility function where pure inertia and pure variety-seeking behaviors can be explained in the model as particular linear cases. The dynamics of the model are fully characterized by standard dynamic programming techniques. The model presents a stationary consumption pattern that can be inertial, where the consumer only buys one product, or a variety-seeking one, where the consumer shifts among varied products.We run some simulations to analyze the consumption paths out of the steady state. Underthe hybrid utility assumption, the consumer behaves inertially among the unfamiliar brandsfor several periods, eventually switching to a variety-seeking behavior when the stationary levels are approached. An empirical analysis is run using scanner databases for three different product categories: fabric softener, saltine cracker, and catsup. Non-linear specifications provide the best fit of the data, as hybrid functional forms are found in all the product categories for most attributes and segments. These results reveal the statistical superiority of the non-linear structure and confirm the gradual trend to seek variety as the level of familiarity with the purchased items increases.

Mouse mammary tumor virus : a model system for regulation of gene expression

Mouse mammary tumor virus (MMTV) kommt prizipiell in zwei Formen vor. Erstens als integierte virale DNA (endogen vererbt), die in allen Zellen der Maus enthalten ist und zweitens als infektiöse Form, bei der sich die DNA nur im Kern von Brustdrüsenzellen integriert. Die erste Form verhält sich wie ein stummes Gen während die zweite Form aktiv ist, durch Glukocorticoide stimuliert wird und zum Mamma-Karzinom führt. Wir haben beide Typen von viralen Genen molekular geklont und durch Transfektion in verschiedene Zellen in Gewebekultur eingeführt. Wir konnten zeigen, dass sowohl die endogene DNA, wie dir infektiöse DNA in transfektieren Zellen aktiv ist und dass die Expression beider Gene durch Glukocorticoide stimuliert wird. Wir konnten die DNA Squenzen, die für dir Homonstimulierung nötig sind, in einem kleinen Fragment der viralen DNA lokalisieren. Bei der Sequenzanalyse dieses DNA-Stückes haben wir ein neues virales Gen entdeckt, das dir Information für ein Protein von ca. 40000 Moleklargewicht enthählt. Mit Hilfe eines Antikörpers suchen wir in verschiedenen Brustdrüsenzellen und Tumoren nach diesem Proetin, dessen Funktion noch nicht bekannt ist.

Prediction of 18-month survival in patients with primary myelodysplastic syndrome. A regression model and scoring system based on the combination of chromosome findings and the Bournemouth score.

The predictive potential of six selected factors was assessed in 72 patients with primary myelodysplastic syndrome using univariate and multivariate logistic regression analysis of survival at 18 months. Factors were age (above median of 69 years), dysplastic features in the three myeloid bone marrow cell lineages, presence of chromosome defects, all metaphases abnormal, double or complex chromosome defects (C23), and a Bournemouth score of 2, 3, or 4 (B234). In the multivariate approach, B234 and C23 proved to be significantly associated with a reduction in the survival probability. The similarity of the regression coefficients associated with these two factors means that they have about the same weight. Consequently, the model was simplified by counting the number of factors (0, 1, or 2) present in each patient, thus generating a scoring system called the Lausanne-Bournemouth score (LB score). The LB score combines the well-recognized and easy-to-use Bournemouth score (B score) with the chromosome defect complexity, C23 constituting an additional indicator of patient outcome. The predicted risk of death within 18 months calculated from the model is as follows: 7.1% (confidence interval: 1.7-24.8) for patients with an LB score of 0, 60.1% (44.7-73.8) for an LB score of 1, and 96.8% (84.5-99.4) for an LB score of 2. The scoring system presented here has several interesting features. The LB score may improve the predictive value of the B score, as it is able to recognize two prognostic groups in the intermediate risk category of patients with B scores of 2 or 3. It has also the ability to identify two distinct prognostic subclasses among RAEB and possibly CMML patients. In addition to its above-described usefulness in the prognostic evaluation, the LB score may bring new insights into the understanding of evolution patterns in MDS. We used the combination of the B score and chromosome complexity to define four classes which may be considered four possible states of myelodysplasia and which describe two distinct evolutional pathways.

Modelling of weather radar echoes from anomalous propagation using a hybrid parabolic equation method and NWP model data

Contamination of weather radar echoes by anomalous propagation (anaprop) mechanisms remains a serious issue in quality control of radar precipitation estimates. Although significant progress has been made identifying clutter due to anaprop there is no unique method that solves the question of data reliability without removing genuine data. The work described here relates to the development of a software application that uses a numerical weather prediction (NWP) model to obtain the temperature, humidity and pressure fields to calculate the three dimensional structure of the atmospheric refractive index structure, from which a physically based prediction of the incidence of clutter can be made. This technique can be used in conjunction with existing methods for clutter removal by modifying parameters of detectors or filters according to the physical evidence for anomalous propagation conditions. The parabolic equation method (PEM) is a well established technique for solving the equations for beam propagation in a non-uniformly stratified atmosphere, but although intrinsically very efficient, is not sufficiently fast to be practicable for near real-time modelling of clutter over the entire area observed by a typical weather radar. We demonstrate a fast hybrid PEM technique that is capable of providing acceptable results in conjunction with a high-resolution terrain elevation model, using a standard desktop personal computer. We discuss the performance of the method and approaches for the improvement of the model profiles in the lowest levels of the troposphere.

The wound watch: an objective staging system for wounds in the diabetic (db/db) mouse model.

As in cancer biology, in wound healing there is a need for objective staging systems to decide for the best treatment and predictors of outcome. We developed in the diabetic (db/db) wound healing model, a staging system, the "wound watch," based on the quantification of angiogenesis and cell proliferation in open wounds. In chronic wounds, there is often a lack of cellular proliferation and angiogenesis that leads to impaired healing. The wound watch addresses this by quantifying the proliferative phase of wound healing in two dimensions (cellular division and angiogenesis). The results are plotted in a two-dimensional graph to monitor the course of healing and compare the response to different treatments.

Current status of a model system: the gene Gp-9 and its association with social organization in fire ants.

The Gp-9 gene in fire ants represents an important model system for studying the evolution of social organization in insects as well as a rich source of information relevant to other major evolutionary topics. An important feature of this system is that polymorphism in social organization is completely associated with allelic variation at Gp-9, such that single-queen colonies (monogyne form) include only inhabitants bearing B-like alleles while multiple-queen colonies (polygyne form) additionally include inhabitants bearing b-like alleles. A recent study of this system by Leal and Ishida (2008) made two major claims, the validity and significance of which we examine here. After reviewing existing literature, analyzing the methods and results of Leal and Ishida (2008), and generating new data from one of their study sites, we conclude that their claim that polygyny can occur in Solenopsis invicta in the U.S.A. in the absence of expression of the b-like allele Gp-9(b) is unfounded. Moreover, we argue that available information on insect OBPs (the family of proteins to which GP-9 belongs), on the evolutionary/population genetics of Gp-9, and on pheromonal/behavioral control of fire ant colony queen number fails to support their view that GP-9 plays no role in the chemosensory-mediated communication that underpins regulation of social organization. Our analyses lead us to conclude that there are no new reasons to question the existing consensus view of the Gp-9 system outlined in Gotzek and Ross (2007).

Hepatitis B and C virus coinfection: a novel model system reveals the absence of direct viral interference.

Coinfection with hepatitis B virus (HBV) and hepatitis C virus (HCV) has been associated with severe liver disease and frequent progression to cirrhosis and hepatocellular carcinoma. Clinical evidence suggests reciprocal replicative suppression of the two viruses, or viral interference. However, interactions between HBV and HCV have been difficult to study due to the lack of appropriate model systems. We have established a novel model system to investigate interactions between HBV and HCV. Stable Huh-7 cell lines inducibly replicating HBV were transfected with selectable HCV replicons or infected with cell culture-derived HCV. In this system, both viruses were found to replicate in the same cell without overt interference. Specific inhibition of one virus did not affect the replication and gene expression of the other. Furthermore, cells harboring replicating HBV could be infected with cell culture-derived HCV, arguing against superinfection exclusion. Finally, cells harboring replicating HBV supported efficient production of infectious HCV. Conclusion: HBV and HCV can replicate in the same cell without evidence for direct interference in vitro. Therefore, the viral interference observed in coinfected patients is probably due to indirect mechanisms mediated by innate and/or adaptive host immune responses. These findings provide new insights into the pathogenesis of HBV-HCV coinfection and may contribute to its clinical management in the future.