Dynamical modeling and analysis of large cellular regulatory networks.

The dynamical analysis of large biological regulatory networks requires the development of scalable methods for mathematical modeling. Following the approach initially introduced by Thomas, we formalize the interactions between the components of a network in terms of discrete variables, functions, and parameters. Model simulations result in directed graphs, called state transition graphs. We are particularly interested in reachability properties and asymptotic behaviors, which correspond to terminal strongly connected components (or "attractors") in the state transition graph. A well-known problem is the exponential increase of the size of state transition graphs with the number of network components, in particular when using the biologically realistic asynchronous updating assumption. To address this problem, we have developed several complementary methods enabling the analysis of the behavior of large and complex logical models: (i) the definition of transition priority classes to simplify the dynamics; (ii) a model reduction method preserving essential dynamical properties, (iii) a novel algorithm to compact state transition graphs and directly generate compressed representations, emphasizing relevant transient and asymptotic dynamical properties. The power of an approach combining these different methods is demonstrated by applying them to a recent multilevel logical model for the network controlling CD4+ T helper cell response to antigen presentation and to a dozen cytokines. This model accounts for the differentiation of canonical Th1 and Th2 lymphocytes, as well as of inflammatory Th17 and regulatory T cells, along with many hybrid subtypes. All these methods have been implemented into the software GINsim, which enables the definition, the analysis, and the simulation of logical regulatory graphs.

Large multi-gene phylogenetic trees of the grasses (Poaceae): progress towards complete tribal and generic level sampling.

In this paper we included a very broad representation of grass family diversity (84% of tribes and 42% of genera). Phylogenetic inference was based on three plastid DNA regions rbcL, matK and trnL-F, using maximum parsimony and Bayesian methods. Our results resolved most of the subfamily relationships within the major clades (BEP and PACCMAD), which had previously been unclear, such as, among others the: (i) BEP and PACCMAD sister relationship, (ii) composition of clades and the sister-relationship of Ehrhartoideae and Bambusoideae + Pooideae, (iii) paraphyly of tribe Bambuseae, (iv) position of Gynerium as sister to Panicoideae, (v) phylogenetic position of Micrairoideae. With the presence of a relatively large amount of missing data, we were able to increase taxon sampling substantially in our analyses from 107 to 295 taxa. However, bootstrap support and to a lesser extent Bayesian inference posterior probabilities were generally lower in analyses involving missing data than those not including them. We produced a fully resolved phylogenetic summary tree for the grass family at subfamily level and indicated the most likely relationships of all included tribes in our analysis.

Serrated polyps of the large intestine: a molecular study comparing sessile serrated adenomas and hyperplastic polyps.

AIMS: To compare the molecular profile of a series of sessile serrated adenomas (SSAs) and hyperplastic polyps (HPs), in order to distinguish these lesions, SSAs having a potential role in the genesis of serrated adenocarcinomas through a serrated pathway in which methylation plays a key role. METHODS AND RESULTS: Twelve HPs and sixteen SSAs of the right and left colon were investigated for microsatellite instability, DNA mismatch repair genes, p53, p16, and beta-catenin expression, MLH1 and p16 (CDKN2A) gene methylation, and KRAS and BRAF mutations. Both SSAs and HPs were microsatellite stable. MLH1 and MSH2 protein silencing, aberrant cytoplasmic expression and methylation of p16 were found to be exclusive to right-sided SSAs. The MLH1 promoter gene was frequently methylated in right-sided SSAs in contrast with HPs. Abnormal p53 and beta-catenin expression was present in both SSAs and HPs. BRAF and KRAS mutation were mutually exclusive, but KRAS mutation was present only in left-sided SSAs and HPs. CONCLUSIONS: HPs and SSAs may be related lesions. However, at least right-sided SSAs differ from left-sided SSAs and HPs in the occurrence of MLH1 and p16 methylation, supporting the hypothesis that SSAs could be precursors of serrated adenocarcinomas.

On problems of calculating energy expenditure and substrate utilization from respiratory exchange data.

Indirect calorimetry based on respiratory exchange measurement has been successfully used from the beginning of the century to obtain an estimate of heat production (energy expenditure) in human subjects and animals. The errors inherent to this classical technique can stem from various sources: 1) model of calculation and assumptions, 2) calorimetric factors used, 3) technical factors and 4) human factors. The physiological and biochemical factors influencing the interpretation of calorimetric data include a change in the size of the bicarbonate and urea pools and the accumulation or loss (via breath, urine or sweat) of intermediary metabolites (gluconeogenesis, ketogenesis). More recently, respiratory gas exchange data have been used to estimate substrate utilization rates in various physiological and metabolic situations (fasting, post-prandial state, etc.). It should be recalled that indirect calorimetry provides an index of overall substrate disappearance rates. This is incorrectly assumed to be equivalent to substrate "oxidation" rates. Unfortunately, there is no adequate golden standard to validate whole body substrate "oxidation" rates, and this contrasts to the "validation" of heat production by indirect calorimetry, through use of direct calorimetry under strict thermal equilibrium conditions. Tracer techniques using stable (or radioactive) isotopes, represent an independent way of assessing substrate utilization rates. When carbohydrate metabolism is measured with both techniques, indirect calorimetry generally provides consistent glucose "oxidation" rates as compared to isotopic tracers, but only when certain metabolic processes (such as gluconeogenesis and lipogenesis) are minimal or / and when the respiratory quotients are not at the extreme of the physiological range. However, it is believed that the tracer techniques underestimate true glucose "oxidation" rates due to the failure to account for glycogenolysis in the tissue storing glucose, since this escapes the systemic circulation. A major advantage of isotopic techniques is that they are able to estimate (given certain assumptions) various metabolic processes (such as gluconeogenesis) in a noninvasive way. Furthermore when, in addition to the 3 macronutrients, a fourth substrate is administered (such as ethanol), isotopic quantification of substrate "oxidation" allows one to eliminate the inherent assumptions made by indirect calorimetry. In conclusion, isotopic tracers techniques and indirect calorimetry should be considered as complementary techniques, in particular since the tracer techniques require the measurement of carbon dioxide production obtained by indirect calorimetry. However, it should be kept in mind that the assessment of substrate oxidation by indirect calorimetry may involve large errors in particular over a short period of time. By indirect calorimetry, energy expenditure (heat production) is calculated with substantially less error than substrate oxidation rates.

Anesthésie locale par tumescence [Tumescent local anesthesia].

Tumescent anesthesia is a local anesthesia produced by the infiltration of a large volume of very dilute anesthetic solution into the subcutaneous tissue. Many surgical procedures (liposuction, facelift, varicose vein surgery, etc.), which were previously performed under general or locoregional anesthesia, can now be achieved with local tumescent anesthesia. This type of anesthesia has many advantages: reduction of both anesthetic risks and surgical complications (bleeding, hematoma), prolonged anesthesia reducing the need for postoperative analgesia, and costs reduction because all these surgical procedures can be performed on an outpatient basis.

Clinical and genetic findings in a large cohort of patients with congenital myopathies due to mutations in the skeletal muscle ryanodine receptor (RYR1) gene

RYR1 mutations are the most common cause of structural congenital myopathies and may exhibit both dominant and recessive inheritance. Histopathological findings are variable and include central cores, multi-minicores, type 1 predominance/ uniformity, fibre type disproportion, increased internal nucleation and fatty and connective tissue. Until recently, diagnostic RYR1 sequencing was limited to mutational hotspots due to the large size of the gene. Since the introduction of full RYR1 sequencing in 2007 we have detected pathogenic mutations in 77 families: 39 had dominant inheritance and 38 recessive inheritance. In some cases with presumably recessive inheritance, only one heterozygous mutation inherited from an asymptomatic parent was identified. Of 28 dominant mutations, 6 were novel; 37 of the 59 recessive mutations were also novel. Dominant mutations were more frequently in recognized hotspot regions, while recessive mutations were distributed throughout the coding sequence. Dominant mutations were predominantly missense, whereas recessive mutations included many nonsense and splice mutations expected to result in reduced RyR1 protein. There was wide clinical variability in patients with both dominant and recessive inheritance. As a group, those with dominant mutations were generally more mildly affected than those with recessive inheritance, who had earlier onset and were weaker with more functional limitations. Extraocular muscle involvement was almost exclusively observed in the recessive group. Bulbar involvement was also more prominent in this group, resulting in a larger number requiring gastrostomy insertion. In conclusion, genomic sequencing of the entire RYR1 leads to the detection of many novel mutations, but may miss large genetic rearrangements in some cases. Assigning pathogenicity to novel mutations is often difficult and interpretation of genetic results in the context of clinical, histological and, increasingly, muscle MRI findings is essential.

Dosing regimen of gentamicin in neonates: evaluation of current guidelines based on a large population pharmacokinetic analysis

Objectives: Gentamicin is among the most commonly prescribed antibiotics in newborns, but large interindividual variability in exposure levels exists. Based on a population pharmacokinetic analysis of a cohort of unselected neonates, we aimed to validate current dosing recommendations from a recent reference guideline (Neofax®). Methods: From 3039 concentrations collected in 994 preterm (median gestational age 32.3 weeks, range 24.2-36.5) and 455 term newborns, treated at the University Hospital of Lausanne between 2006 and 2011, a population pharmacokinetic analysis was performed with NONMEM®. Model-based simulations were used to assess the ability of dosing regimens to bring concentrations into targets: trough ≤ 1mg/L and peak ~ 8mg/L. Results: A two-compartment model best characterized gentamicin pharmacokinetics. Model parameters are presented in the table. Body weight, gestational age and postnatal age positively influence clearance, which decreases under dopamine administration. Body weight and gestational age influence the distribution volume. Model based simulations confirm that preterm infants need doses superior to 4 mg/kg, and extended dosage intervals, up to 48 hours for very preterm newborns, whereas most term newborns would achieve adequate exposure under 4 mg/kg q. 24 h. More than 90% of neonates would achieve trough concentrations below 2 mg/L and peaks above 6 mg/L following most recent guidelines. Conclusions: Simulated gentamicin exposure demonstrates good accordance with recent dosing recommendations for target concentration achievement.

Application of a roughness-length representation to parameterize energy loss in 3-D numerical simulations of large rivers

Recent technological advances in remote sensing have enabled investigation of the morphodynamics and hydrodynamics of large rivers. However, measuring topography and flow in these very large rivers is time consuming and thus often constrains the spatial resolution and reach-length scales that can be monitored. Similar constraints exist for computational fluid dynamics (CFD) studies of large rivers, requiring maximization of mesh-or grid-cell dimensions and implying a reduction in the representation of bedform-roughness elements that are of the order of a model grid cell or less, even if they are represented in available topographic data. These ``subgrid'' elements must be parameterized, and this paper applies and considers the impact of roughness-length treatments that include the effect of bed roughness due to ``unmeasured'' topography. CFD predictions were found to be sensitive to the roughness-length specification. Model optimization was based on acoustic Doppler current profiler measurements and estimates of the water surface slope for a variety of roughness lengths. This proved difficult as the metrics used to assess optimal model performance diverged due to the effects of large bedforms that are not well parameterized in roughness-length treatments. However, the general spatial flow patterns are effectively predicted by the model. Changes in roughness length were shown to have a major impact upon flow routing at the channel scale. The results also indicate an absence of secondary flow circulation cells in the reached studied, and suggest simpler two-dimensional models may have great utility in the investigation of flow within large rivers. Citation: Sandbach, S. D. et al. (2012), Application of a roughness-length representation to parameterize energy loss in 3-D numerical simulations of large rivers, Water Resour. Res., 48, W12501, doi: 10.1029/2011WR011284.

Determinants of generic drug substitution in Switzerland.

Background: Since generic drugs have the same therapeutic effect as the original formulation but at generally lower costs, their use should be more heavily promoted. However, a considerable number of barriers to their wider use have been observed in many countries. The present study examines the influence of patients, physicians and certain characteristics of the generics' market on generic substitution in Switzerland.Methods: We used reimbursement claims' data submitted to a large health insurer by insured individuals living in one of Switzerland's three linguistic regions during 2003. All dispensed drugs studied here were substitutable. The outcome (use of a generic or not) was modelled by logistic regression, adjusted for patients' characteristics (gender, age, treatment complexity, substitution groups) and with several variables describing reimbursement incentives (deductible, co-payments) and the generics' market (prices, packaging, co-branded original, number of available generics, etc.).Results: The overall generics' substitution rate for 173,212 dispensed prescriptions was 31%, though this varied considerably across cantons. Poor health status (older patients, complex treatments) was associated with lower generic use. Higher rates were associated with higher out-of-pocket costs, greater price differences between the original and the generic, and with the number of generics on the market, while reformulation and repackaging were associated with lower rates. The substitution rate was 13% lower among hospital physicians. The adoption of the prescribing practices of the canton with the highest substitution rate would increase substitution in other cantons to as much as 26%.Conclusions: Patient health status explained a part of the reluctance to substitute an original formulation by a generic. Economic incentives were efficient, but with a moderate global effect. The huge interregional differences indicated that prescribing behaviours and beliefs are probably the main determinant of generic substitution.

Treatment of large bone defects with the Ilizarov technique.

Between 1985 and 1990 we treated 11 large segmental bone defects (average 6.7 cm) in ten patients with the Ilizarov technique. Open fractures, type III according to Gustilo, represented the largest group (8 of 11 cases). The average delay before the Ilizarov technique was initiated was 8.9 months. The external fixator was usually maintained for 1 year. Bone regeneration was obtained in every case. Consolidation was not fulfilled with this technique in three cases. The complications observed were one refracture, four leg-length discrepancies (average 1.5 cm), and five axial deformities exceeding 5 degrees. No pin-track infection was observed. In our limited series of four type IIIC open fractures treated by the Ilizarov technique, no patients required amputation. The Ilizarov technique is particularly useful in the treatment of large bone defects, without major complications, especially if there is an adequate initial debridement.

A modular approach for integrative analysis of large-scale gene-expression and drug-response data

High-throughput technologies are now used to generate more than one type of data from the same biological samples. To properly integrate such data, we propose using co-modules, which describe coherent patterns across paired data sets, and conceive several modular methods for their identification. We first test these methods using in silico data, demonstrating that the integrative scheme of our Ping-Pong Algorithm uncovers drug-gene associations more accurately when considering noisy or complex data. Second, we provide an extensive comparative study using the gene-expression and drug-response data from the NCI-60 cell lines. Using information from the DrugBank and the Connectivity Map databases we show that the Ping-Pong Algorithm predicts drug-gene associations significantly better than other methods. Co-modules provide insights into possible mechanisms of action for a wide range of drugs and suggest new targets for therapy

Des céphalées, cinq avis spécialisés: les pièges diagnostiques pour la maladie de Horton [Headaches, five expert opinions: the pitfalls involved in the diagnosis of giant cell arteritis].

Giant cell arteritis (GCA) (or Horton's disease) is a systemic disease affecting the vessels of medium and large sizes. The incidence increases with age (the disease develops rarely before age 50) and the etiology remains unknown. Clinical manifestations may vary (including asthenia, temporal headache, visual disturbances, etc.) and GCA can potentially lead to dramatic consequences (permanent loss of vision). Although some anomalies in the investigations may help in the diagnosis of GCA, research and confirmation of the diagnosis of GCA may be difficult, especially when the symptoms presented by patients are spread out in time and appear to be nonspecific at first.

Thymic selection generates a large T cell pool recognizing a self-peptide in humans.

The low frequency of self-peptide-specific T cells in the human preimmune repertoire has so far precluded their direct evaluation. Here, we report an unexpected high frequency of T cells specific for the self-antigen Melan-A/MART-1 in CD8 single-positive thymocytes from human histocompatibility leukocyte antigen-A2 healthy individuals, which is maintained in the peripheral blood of newborns and adults. Postthymic replicative history of Melan-A/MART-1-specific CD8 T cells was independently assessed by quantifying T cell receptor excision circles and telomere length ex vivo. We provide direct evidence that the large T cell pool specific for the self-antigen Melan-A/MART-1 is mostly generated by thymic output of a high number of precursors. This represents the only known naive self-peptide-specific T cell repertoire directly accessible in humans.