Synchronous versus asynchronous modeling of gene regulatory networks.

MOTIVATION: In silico modeling of gene regulatory networks has gained some momentum recently due to increased interest in analyzing the dynamics of biological systems. This has been further facilitated by the increasing availability of experimental data on gene-gene, protein-protein and gene-protein interactions. The two dynamical properties that are often experimentally testable are perturbations and stable steady states. Although a lot of work has been done on the identification of steady states, not much work has been reported on in silico modeling of cellular differentiation processes. RESULTS: In this manuscript, we provide algorithms based on reduced ordered binary decision diagrams (ROBDDs) for Boolean modeling of gene regulatory networks. Algorithms for synchronous and asynchronous transition models have been proposed and their corresponding computational properties have been analyzed. These algorithms allow users to compute cyclic attractors of large networks that are currently not feasible using existing software. Hereby we provide a framework to analyze the effect of multiple gene perturbation protocols, and their effect on cell differentiation processes. These algorithms were validated on the T-helper model showing the correct steady state identification and Th1-Th2 cellular differentiation process. AVAILABILITY: The software binaries for Windows and Linux platforms can be downloaded from http://si2.epfl.ch/~garg/genysis.html.

Repeated injections of 131I-rituximab show patient-specific stable biodistribution and tissue kinetics.

PURPOSE: It is generally assumed that the biodistribution and pharmacokinetics of radiolabelled antibodies remain similar between dosimetric and therapeutic injections in radioimmunotherapy. However, circulation half-lives of unlabelled rituximab have been reported to increase progressively after the weekly injections of standard therapy doses. The aim of this study was to evaluate the evolution of the pharmacokinetics of repeated 131I-rituximab injections during treatment with unlabelled rituximab in patients with non-Hodgkin's lymphoma (NHL). METHODS: Patients received standard weekly therapy with rituximab (375 mg/m2) for 4 weeks and a fifth injection at 7 or 8 weeks. Each patient had three additional injections of 185 MBq 131I-rituximab in either treatment weeks 1, 3 and 7 (two patients) or weeks 2, 4 and 8 (two patients). The 12 radiolabelled antibody injections were followed by three whole-body (WB) scintigraphic studies during 1 week and blood sampling on the same occasions. Additional WB scans were performed after 2 and 4 weeks post 131I-rituximab injection prior to the second and third injections, respectively. RESULTS: A single exponential radioactivity decrease for WB, liver, spleen, kidneys and heart was observed. Biodistribution and half-lives were patient specific, and without significant change after the second or third injection compared with the first one. Blood T(1/2)beta, calculated from the sequential blood samples and fitted to a bi-exponential curve, was similar to the T(1/2) of heart and liver but shorter than that of WB and kidneys. Effective radiation dose calculated from attenuation-corrected WB scans and blood using Mirdose3.1 was 0.53+0.05 mSv/MBq (range 0.48-0.59 mSv/MBq). Radiation dose was highest for spleen and kidneys, followed by heart and liver. CONCLUSION: These results show that the biodistribution and tissue kinetics of 131I-rituximab, while specific to each patient, remained constant during unlabelled antibody therapy. RIT radiation doses can therefore be reliably extrapolated from a preceding dosimetry study.

Oligo-Miocene extensional tectonics and fluid flow across the Northern Snake Range detachment system, Nevada

The Northern Snake Range (Nevada) represents a spectacular example of a metamorphic core complex and exposes a complete section from the mylonitic footwall into the hanging wall of a fossil detachment system. Paired geochronological and stable isotopic data of mylonitic quartzite within the detachment footwall reveal that ductile deformation and infiltration of meteoric fluids occurred between 27 and 23 Ma. Ar-40/Ar-39 ages display complex recrystallization-cooling relationships but decrease systematically from 26.9 +/- 0.2 Ma at the top to 21.3 +/- 0.2 Ma at the bottom of footwall mylonite. Hydrogen isotope (delta D) values in white mica are very low (-150 to -145 %) within the top 80-90 m of detachment footwall, in contrast to values obtained from the deeper part of the section where values range from -77 to -64 %, suggesting that time-integrated interaction between rock and meteoric fluid was restricted to the uppermost part of the mylonitic footwall. Pervasive mica-water hydrogen isotope exchange is difficult to reconcile with models of Ar-40 loss during mylonitization solely by volume diffusion. Rather, we interpret the Ar-40/Ar-39 ages of white mica with low-delta D values to date syn-mylonitic hydrogen and argon isotope exchange, and we conclude that the hydrothermal system of the Northern Snake Range was active during late Oligocene (27-23 Ma) and has been exhumed by the combined effects of ductile strain, extensional detachment faulting, and erosion.

Age-specific incidence of all neoplasms after colorectal cancer.

PURPOSE: Patients diagnosed with a specific neoplasm tend to have a subsequent excess risk of the same neoplasm. The age incidence of a second neoplasm at the same site is approximately constant with age, and consequently the relative risk is greater at younger age. It is unclear whether such a line of reasoning can be extended from a specific neoplasm to the incidence of all neoplasms in subjects diagnosed with a defined neoplasm. METHODS: We considered the age-specific incidence of all non-hormone-related epithelial neoplasms after a first primary colorectal cancer (n = 9542) in the Vaud Cancer Registry data set. RESULTS: In subjects with a previous colorectal cancer, the incidence rate of all other epithelial non-hormone-related cancers was stable around 800 per 100,000 between age 30 and 60 years, and rose only about twofold to reach 1685 at age 70 to 79 years and 1826 per 100,000 at age 80 years or older. After excluding synchronous cancers, the rise was only about 1.5-fold, that is, from about 700 to 1000. In the general population, the incidence rate of all epithelial non-hormone-related cancers was 29 per 100,000 at age 30 to 39 years, and rose 30-fold to 883 per 100,000 at age 70 to 79 years. Excluding colorectal cancers, the rise of all non-hormone-related cancers was from 360 per 100,000 at age 40 to 49 years to 940 at age 70 to 79 years after colorectal cancer, and from 90 to 636 per 100,000 in the general population (i.e., 2.6- vs. 7.1-fold). CONCLUSIONS: The rise of incidence with age of all epithelial non-hormone-related second cancers after colorectal cancer is much smaller than in the general population. This can possibly be related to the occurrence of a single mutational event in a population of susceptible individuals, although alternative models are plausible within the complexity of the process of carcinogenesis.

Developmental and pathological lymphangiogenesis: from models to human disease.

The lymphatic vascular system, the body's second vascular system present in vertebrates, has emerged in recent years as a crucial player in normal and pathological processes. It participates in the maintenance of normal tissue fluid balance, the immune functions of cellular and antigen trafficking and absorption of fatty acids and lipid-soluble vitamins in the gut. Recent scientific discoveries have highlighted the role of lymphatic system in a number of pathologic conditions, including lymphedema, inflammatory diseases, and tumor metastasis. Development of genetically modified animal models, identification of lymphatic endothelial specific markers and regulators coupled with technological advances such as high-resolution imaging and genome-wide approaches have been instrumental in understanding the major steps controlling growth and remodeling of lymphatic vessels. This review highlights the recent insights and developments in the field of lymphatic vascular biology.

Evaluating the ability of habitat suitability models to predict species presences.

Models predicting species spatial distribution are increasingly applied to wildlife management issues, emphasising the need for reliable methods to evaluate the accuracy of their predictions. As many available datasets (e.g. museums, herbariums, atlas) do not provide reliable information about species absences, several presence-only based analyses have been developed. However, methods to evaluate the accuracy of their predictions are few and have never been validated. The aim of this paper is to compare existing and new presenceonly evaluators to usual presence/absence measures. We use a reliable, diverse, presence/absence dataset of 114 plant species to test how common presence/absence indices (Kappa, MaxKappa, AUC, adjusted D-2) compare to presenceonly measures (AVI, CVI, Boyce index) for evaluating generalised linear models (GLM). Moreover we propose a new, threshold-independent evaluator, which we call "continuous Boyce index". All indices were implemented in the B10MAPPER software. We show that the presence-only evaluators are fairly correlated (p > 0.7) to the presence/absence ones. The Boyce indices are closer to AUC than to MaxKappa and are fairly insensitive to species prevalence. In addition, the Boyce indices provide predicted-toexpected ratio curves that offer further insights into the model quality: robustness, habitat suitability resolution and deviation from randomness. This information helps reclassifying predicted maps into meaningful habitat suitability classes. The continuous Boyce index is thus both a complement to usual evaluation of presence/absence models and a reliable measure of presence-only based predictions.

Sensitivity of predictive species distribution models to change in grain size

Predictive species distribution modelling (SDM) has become an essential tool in biodiversity conservation and management. The choice of grain size (resolution) of environmental layers used in modelling is one important factor that may affect predictions. We applied 10 distinct modelling techniques to presence-only data for 50 species in five different regions, to test whether: (1) a 10-fold coarsening of resolution affects predictive performance of SDMs, and (2) any observed effects are dependent on the type of region, modelling technique, or species considered. Results show that a 10 times change in grain size does not severely affect predictions from species distribution models. The overall trend is towards degradation of model performance, but improvement can also be observed. Changing grain size does not equally affect models across regions, techniques, and species types. The strongest effect is on regions and species types, with tree species in the data sets (regions) with highest locational accuracy being most affected. Changing grain size had little influence on the ranking of techniques: boosted regression trees remain best at both resolutions. The number of occurrences used for model training had an important effect, with larger sample sizes resulting in better models, which tended to be more sensitive to grain. Effect of grain change was only noticeable for models reaching sufficient performance and/or with initial data that have an intrinsic error smaller than the coarser grain size.

Long-term health-related quality of life and walking capacity of lung recipients with and without bronchiolitis obliterans syndrome.

BACKGROUND: Outcome after lung transplantation (LTx) is affected by the onset of bronchiolitis obliterans syndrome (BOS) and lung function decline. Reduced health-related quality of life (HRQL) and physical mobility have been shown in patients developing BOS, but the impact on the capacity to walk is unknown. We aimed to compare the long-term HRQL and 6-minute walk test (6MWT) between lung recipients affected or not by BOS Grade > or =2. METHODS: Fifty-eight patients were prospectively followed for 5.6 +/- 2.9 years after LTx. Assessments included the St George's Respiratory Questionnaire (SGRQ) and the 6MWT, which were performed yearly. Moreover, clinical complications were recorded to estimate the proportion of the follow-up time lived without clinical intercurrences after transplant. Analyses were performed using adjusted linear regression and repeated-measures analysis of variance. RESULTS: BOS was a significant predictor of lower SGRQ scores (p < 0.01) and reduced time free of clinical complications (p = 0.001), but not of 6MWT distance (p = 0.12). At 7 years post-transplant, results were: 69.0 +/- 21.8% vs 86.9 +/- 5.6%, p < 0.05 (SGRQ); 58.5 +/- 21.6% vs 88.7 +/- 11.4%, p < 0.01 (proportion of time lived without clinical complications); and 82.2 +/- 10.9% vs 91.9 +/- 14.2%, p = 0.27 (percent of predicted 6MWT), respectively, for patients with BOS and without BOS. CONCLUSIONS: Despite significantly less time lived without clinical complications and progressive decline of self-reported health status, the capacity to walk of patients affected by BOS remained relatively stable over time. These findings may indicate that the development of moderate to severe BOS does not prevent lung recipients from walking independently and pursuing an autonomous life.

Evaluation of tumor vascularisation in two rat sarcoma models for studying isolated lung perfusion : Injection route determines the origin of tumour vessels

Résumé Introduction: La perfusion isolée cytostatique du poumon est une technique attractive qui permet l'administration des doses élevées d'un agent cytostatique tout en épargnant dans la mesure du possible la circulation systémique. Cependant, la perfusion de l'artère pulmonaire risque d'épargner le territoire pulmonaire vascularisé par l'intermédiaire des artères bronchiques, ce qui pourrait diminuer l'efficacité de ce traitement au cas où la lésion ciblée est vascularisée par les artères bronchiques. Ce travail est destiné au développement d'un modèle tumoral au niveau des poumons de rongeur (rat) porteur d'un sarcome pulmonaire afin de déterminer si la voie d'injection des cellules tumorales (intraveineuse, versus intratumorale) influencera la vascularisation des tumeurs (provenant du système artères pulmonaires ou artères bronchiques). Méthod: Des tumeurs de sarcomes pulmonaires ont été générées par injection d'une suspension cellulaire de sarcome, soit par injection intraveineuse, soit directement dans le parenchyme pulmonaire par thoracotomie. Ensuite, une perfusion isolée du poumon porteur de la tumeur à l'aide de l'encre a été effectuée, soit par l'artère pulmonaire, soit par le système des artères bronchiques. La distribution de l'encre dans les vaisseaux tumoraux ainsi que dans les vaisseaux non tumoraux du poumon adjacent a été investiguée à l'aide d'une analyse histologique des poumons perfusés. Résultat: L'administration intraveineuse et intratumorale de la suspension de cellules tumorales résulte en des tumeurs similaires sur le plan histologique. Néanmoins, l'injection intra-parenchymateuse démontre des tumeurs plus homogènes et avec un développement plus prédictible, était associée à une survie plus longue qu'après injection intraveineuse. Les analyses histologiques après perfusion isolée à l'aide de l'encre démontre que les tumeurs résultant de l'injection intraveineuse ont développé une vascularisation se basant sur le système d'artères pulmonaires tandis que les tumeurs émergeant après injection intraparenchymateuse ont développé une vascularisation provenant du système des artères bronchiques. Conclusion: Ce travail démontre pour la première fois l'importance du mode de génération de tumeurs pulmonaires en ce qui concerne leur future vascularisation, ce qui pourrait avoir un impact sur leur traitement par perfusion isolée du poumon. Abstract Isolated cytostatic lung perfusion (ILP) is an attractive technique allowing delivery of a high-dose of cytostatic agents to the lungs while limiting systemic toxicity. In developing a rat model of ILP, we have analysed the effect of the route of tumour cell injection on the source of tumour vessels. Pulmonary sarcomas were estab¬lished by injecting a sarcoma cell suspension either by the intravenous (i.v.) route or directly into the lung paren¬chyma. Ink perfusion through either pulmonary artery (PA) or bronchial arteries (BA) was performed and the characteristics of the tumour deposits defined. i.v. and direct injection methods induced pulmonary sarcoma nodules, with similar histological features. The intraparenchymal injection of tumour cells resulted in more reli¬able and reproducible tumour growth and was associat¬ed with a longer survival of the animals. i.v. injected tumours developed a PA-derived vascular tree whereas directly injected tumours developed a BA-derived vasculature.

Avian erythrocytes have functional mitochondria, opening novel perspectives for birds as animal models in the study of ageing.

BACKGROUND: In contrast to mammalian erythrocytes, which have lost their nucleus and mitochondria during maturation, the erythrocytes of almost all other vertebrate species are nucleated throughout their lifespan. Little research has been done however to test for the presence and functionality of mitochondria in these cells, especially for birds. Here, we investigated those two points in erythrocytes of one common avian model: the zebra finch (Taeniopygia guttata). RESULTS: Transmission electron microscopy showed the presence of mitochondria in erythrocytes of this small passerine bird, especially after removal of haemoglobin interferences. High-resolution respirometry revealed increased or decreased rates of oxygen consumption by erythrocytes in response to the addition of respiratory chain substrates or inhibitors, respectively. Fluorometric assays confirmed the production of mitochondrial superoxide by avian erythrocytes. Interestingly, measurements of plasmatic oxidative markers indicated lower oxidative stress in blood of the zebra finch compared to a size-matched mammalian model, the mouse. CONCLUSIONS: Altogether, those findings demonstrate that avian erythrocytes possess functional mitochondria in terms of respiratory activities and reactive oxygen species (ROS) production. Interestingly, since blood oxidative stress was lower for our avian model compared to a size-matched mammalian, our results also challenge the idea that mitochondrial ROS production could have been one actor leading to this loss during the course of evolution. Opportunities to assess mitochondrial functioning in avian erythrocytes open new perspectives in the use of birds as models for longitudinal studies of ageing via lifelong blood sampling of the same subjects.

The MIGCLIM R package - seamless integration of dispersal constraints into projections of species distribution models

The MIGCLIM R package is a function library for the open source R software that enables the implementation of species-specific dispersal constraints into projections of species distribution models under environmental change and/or landscape fragmentation scenarios. The model is based on a cellular automaton and the basic modeling unit is a cell that is inhabited or not. Model parameters include dispersal distance and kernel, long distance dispersal, barriers to dispersal, propagule production potential and habitat invasibility. The MIGCLIM R package has been designed to be highly flexible in the parameter values it accepts, and to offer good compatibility with existing species distribution modeling software. Possible applications include the projection of future species distributions under environmental change conditions and modeling the spread of invasive species.