Modelling Time-Dependent Acquisition of Positional Information

Theoretical and computational modelling are crucial to understand dynamics of embryonic development. In this tutorial chapter, we describe two models of gene networks performing time-dependent acquisition of positional information under control of a dynamic morphogen: a toy-model of a bistable gene under control of a morphogen, allowing for the numerical computation of a simple Waddington's epigenetic landscape, and a recently published model of gap genes in Tribolium under control of multiple enhancers. We present detailed commented implementations of the models using python and jupyter notebooks.

Untangling the Hairball: Fitness-Based Asymptotic Reduction of Biological Networks

Complex mathematical models of interaction networks are routinely used for prediction in systems biology. However, it is difficult to reconcile network complexities with a formal understanding of their behavior. Here, we propose a simple procedure (called ϕ¯) to reduce biological models to functional submodules, using statistical mechanics of complex systems combined with a fitness-based approach inspired by in silico evolution. The ϕ¯ algorithm works by putting parameters or combination of parameters to some asymptotic limit, while keeping (or slightly improving) the model performance, and requires parameter symmetry breaking for more complex models. We illustrate ϕ¯ on biochemical adaptation and on different models of immune recognition by T cells. An intractable model of immune recognition with close to a hundred individual transition rates is reduced to a simple two-parameter model. The ϕ¯ algorithm extracts three different mechanisms for early immune recognition, and automatically discovers similar functional modules in different models of the same process, allowing for model classification and comparison. Our procedure can be applied to biological networks based on rate equations using a fitness function that quantifies phenotypic performance.

[Treatment of wrinkles of the upper lip by emulsified fat or "Nanofat": Biological and clinical study about 4 cases]

OBJECTIVE

Emulsified fat injection showed its interest in aesthetic facial surgery. The adipose tissue harvested is mechanically emulsified and filtered. The suspension obtained is injected into the dermis through small diameter needles (27 to 30 gauges). The objective of our study was to evaluate the biological composition of emulsified fat and its clinical effectiveness in the treatment of peri-oral wrinkles in 4 patients aged 50 to 59 years.

MATERIAL AND METHOD

Each patient received an intradermal injection of emulsified fat in the peri-oral wrinkles prepared from abdominal fat under local anesthesia. The cell viability, stromal vascular fraction (FVS) composition in emulsified fat and the adipocyte differentiation capacity of mesenchymal stem cells (MSC) were studied. The clinical results were evaluated by standardized photographs, 3D microphotography, confocal microscopy, and self-evaluation of patient satisfaction over a period of 4 months.

RESULTS

The biological study of the emulsified fat found a lysis of all the adipocytes. The mean number of FVS cells was 126,330±2758 cells by cc of emulsified fat with preserved cell viability (85.1±6.84 %) and a good proportion of regeneratives cells (18.77±6.2 %). The clinical study found a tendency to decrease the volume of wrinkles on standardized photography and 3D microphotography no significative. Patients were satisfied with treatment with an average score of 7±1.15/10 to 4 months.

CONCLUSION

Intradermal injection of emulsified fat seems to be an interesting treatment of face wrinkles. Our study has shown its safety, but additional studies seems necessary to confirm its clinical efficacy.

Analysis of the prophages carried by human infecting isolates provides new insight into the evolution of Group B Streptococcus species

OBJECTIVES

Group B Streptococcus (GBS) emerged in the 1970s as a major cause of neonatal infections, and has been increasingly associated with infections in adults since the 1990s. Prophages have been suspected to have driven these epidemiological trends. We have characterized the prophages harboured by 275 human GBS isolates belonging to the major lineages.

METHODS

We applied whole genome sequencing (WGS) to 14 isolates representative of the diversity within GBS species, located and identified their prophages. Using prediction tools, we searched for prophage elements potentially involved with the ability of GBS to infect humans. Using the data obtained by WGS, we designed a PCR-based tool and studied the prophage content of 275 isolates.

RESULTS

WGS of the 14 isolates revealed 22 prophages (i) distributed into six groups (A-F), (ii) similar to phages and prophages from GBS and non-GBS streptococci recovered from livestock, and (iii) carrying genes encoding factors previously associated with host adaptation and virulence. PCR-based detection of prophages revealed the presence of at least one prophage in 72.4% of the 275 isolates and a significant association between neonatal infecting isolates and prophages C, and between adult infecting isolates and prophages A.

CONCLUSIONS

Our results suggest that prophages (possibly animal-associated) have conditioned bacterial adaptation and ability to cause infections in neonates and adults, and support a role of lysogeny with the emergence of GBS as a pathogen in human.

Sorting at embryonic boundaries requires high heterotypic interfacial tension

The establishment of sharp boundaries is essential for segregation of embryonic tissues during development, but the underlying mechanism of cell sorting has remained unclear. Opposing hypotheses have been proposed, either based on global tissue adhesive or contractile properties or on local signalling through cell contact cues. Here we use ectoderm-mesoderm separation in Xenopus to directly evaluate the role of these various parameters. We find that ephrin-Eph-based repulsion is very effective at inducing and maintaining separation, whereas differences in adhesion or contractility have surprisingly little impact. Computer simulations support and generalise our experimental results, showing that a high heterotypic interfacial tension between tissues is key to their segregation. We propose a unifying model, in which conditions of sorting previously considered as driven by differential adhesion/tension should be viewed as suboptimal cases of heterotypic interfacial tension.The mechanisms that cause different cells to segregate into distinct tissues are unclear. Here the authors show in Xenopus that formation of a boundary between two tissues is driven by local tension along the interface rather than by global differences in adhesion or cortical contractility.

Looking at hospitalized persons throughout the prism of the handicap

OBJECTIVE

To describe the disability status of non-selected hospitalized persons.

METHODS AND FINDINGS

We conducted a cross-sectional survey to assess activity limitations of every person older than 18 years hospitalized in a regional university hospital covering all medical fields. Evaluators rated, on a scale from 0 to 4, 22 selected items of the International Classification of Functioning (ICF), covering the 6 following domains: learning and applying knowledge, general tasks and demands, communication, mobility, self-care, and interpersonal interactions and relationships. Univariate and multivariate analyses were performed to analyze the prevalence, severity and profile of the handicap in terms of sociodemographic characteristics and care pathways.

RESULTS

Among 1572 eligible persons, 1267 (81%) were surveyed (mean age 62.7±20.4years; 655 males [51.7%]). Overall, 82% showed at least one activity limitation. For 52%, disability was severe or total for at least one ICF item. Prevalence of disabilities was higher for mobility (75%) and self-care domains (63%). Disability was strongly related to age: age older than 80years versus 18 to 44years (OR=12.8 95% CI 6.4-27.9]; P<0.01). Disability was associated with hospitalization in rehabilitation units (96%; OR=4.3 [95% CI 2.2-5.3]; P<0.01). Severe disability was associated with hospitalization in critical care units (OR=6.7 [CI 3.2-15.1]; P<0.001) and psychiatry units (OR=5.3 [CI 2.7-11.4]; P<0.001).

CONCLUSION

Handicap was common in hospitalized persons, involving all 6 tested ICF activity domains, particularly mobility and self-care. This study alerts care givers, hospital administrators, and in general, people influencing health policies about the need to plan actions to reduce activity limitations of hospitalized persons, whatever the cause of the hospitalization.

Antimicrobial activity of ceftaroline against methicillin-resistant Staphylococcus aureus (MRSA) isolates collected in 2013-2014 at the Geneva University Hospitals

Ceftaroline is a broad-spectrum antibiotic with activity against methicillin-resistant Staphylococcus aureus (MRSA) strains. Ceftaroline susceptibility of an MRSA set archived between 1994 and 2003 in the Geneva University Hospitals detected a high percentage (66 %) of ceftaroline resistance in clonotypes ST228 and ST247 and correlated with mutations in PBP2a. The ceftaroline mechanism of action is based on the inhibition of PBP2a; thus, the identification of PBP2a mutations of recently circulating clonotypes in our institution was investigated. We analyzed ceftaroline susceptibility in MRSA isolates (2013 and 2014) and established that resistant strains correlated with PBP2a mutations and specific clonotypes. Ninety-six MRSA strains were analyzed from independent patients and were isolated from blood cultures (23 %), deep infections (38.5 %), and superficial (skin or wound) infections (38.5 %). This sample showed a ceftaroline minimum inhibitory concentration (MIC) range between 0.25 and 2 μg/ml and disk diameters ranging from 10 to 30 mm, with a majority of strains showing diameters ≥20 mm. Based on the European Committee on Antimicrobial Susceptibility Testing (EUCAST) breakpoints, 76 % (73/96) of isolates showed susceptibility to ceftaroline. Nevertheless, we still observed 24 % (23/96) of resistant isolates (MIC = 2 μg/ml). All resistant isolates were assigned to clonotype ST228 and carried the N146K mutation in PBP2a. Only two ST228 isolates showed ceftaroline susceptibility. The decreasing percentage of ceftaroline-resistant isolates in our hospital can be explained by the decline of ST228 clonotype circulating in our hospital since 2008. We present evidence that ceftaroline is active against recent MRSA strains from our hospital; however, the presence of PBP2a variants in particular clonotypes may affect ceftaroline efficacy.

Multiscale quantification of tissue behavior during amniote embryo axis elongation

Embryonic axis elongation is a complex multi-tissue morphogenetic process responsible for the formation of the posterior part of the amniote body. How movements and growth are coordinated between the different posterior tissues (e.g. neural tube, axial and paraxial mesoderm, lateral plate, ectoderm, endoderm) to drive axis morphogenesis remain largely unknown. Here, we use quail embryos to quantify cell behavior and tissue movements during elongation. We quantify the tissue-specific contribution to axis elongation by using 3D volumetric techniques, then quantify tissue-specific parameters such as cell density and proliferation. To study cell behavior at a multi-tissue scale, we used high-resolution 4D imaging of transgenic quail embryos expressing fluorescent proteins. We developed specific tracking and image analysis techniques to analyze cell motion and compute tissue deformations in 4D. This analysis reveals extensive sliding between tissues during axis extension. Further quantification of tissue tectonics showed patterns of rotations, contractions and expansions, which are coherent with the multi-tissue behavior observed previously. Our approach defines a quantitative and multiscale method to analyze the coordination between tissue behaviors during early vertebrate embryo morphogenetic events.

Phenotypic spandrel: absolute discrimination and ligand antagonism

We consider the general problem of sensitive and specific discrimination between biochemical species. An important instance is immune discrimination between self and not-self, where it is also observed experimentally that ligands just below the discrimination threshold negatively impact response, a phenomenon called antagonism. We characterize mathematically the generic properties of such discrimination, first relating it to biochemical adaptation. Then, based on basic biochemical rules, we establish that, surprisingly, antagonism is a generic consequence of any strictly specific discrimination made independently from ligand concentration. Thus antagonism constitutes a 'phenotypic spandrel': a phenotype existing as a necessary by-product of another phenotype. We exhibit a simple analytic model of discrimination displaying antagonism, where antagonism strength is linear in distance from the detection threshold. This contrasts with traditional proofreading based models where antagonism vanishes far from threshold and thus displays an inverted hierarchy of antagonism compared to simpler models. The phenotypic spandrel studied here is expected to structure many decision pathways such as immune detection mediated by TCRs and FCϵRIs, as well as endocrine signalling/disruption.

Oxidants at the Surface of Mars: A Review in Light of Recent Exploration Results

In 1976, the Viking landers carried out the most comprehensive search for organics and microbial life in the martian regolith. Their results indicate that Mars' surface is lifeless and, surprisingly, depleted in organics at part-per-billion levels. Several biology experiments on the Viking landers gave controversial results that have since been explained by the presence of oxidizing agents on the surface of Mars. These oxidants may degrade abiotic or biological organics, resulting in their nondetection in the regolith. As several exploration missions currently focus on the detection of organics on Mars (or will do so in the near future), knowledge of the oxidative state of the surface is fundamental. It will allow for determination of the capability of organics to survive on a geological timescale, the most favorable places to seek them, and the best methods to process the samples collected at the surface. With this aim, we review the main oxidants assumed to be present on Mars, their possible formation pathways, and those laboratory studies in which their reactivity with organics under Mars-like conditions has been evaluated. Among the oxidants assumed to be present on Mars, only four have been detected so far: perchlorate ions (ClO₄^-) in salts, hydrogen peroxide (H₂O₂) in the atmosphere, and clays and metal oxides composing surface minerals. Clays have been suggested as catalysts for the oxidation of organics but are treated as oxidants in the following to keep the structure of this article straightforward. This work provides an insight into the oxidizing potential of the surface of Mars and an estimate of the stability of organic matter in an oxidizing environment. Key Words: Mars surface-Astrobiology-Oxidant-Chemical reactions. Astrobiology 16, 977-996.

Imipenem heteroresistance in nontypeable Haemophilus influenzae is linked to a combination of altered PBP3, slow drug influx and direct efflux regulation

OBJECTIVE

To investigate the potential roles of PBPs, efflux pumps and slow drug influx for imipenem heteroresistance in nontypeable Haemophilus influenzae (NTHi).

METHODS

Fifty-nine NTHi clinical isolates examined in this study were collected at Geneva University Hospitals between 2009 and 2014. Alterations in PBPs were investigated by gene sequencing. To evaluate the affinities of the PBPs to imipenem, steady-state concentration-response experiments were carried out using imipenem in a competition assay with Bocillin-FL. The effect of the carbonyl cyanide m-chlorophenylhydrazone (CCCP) on imipenem susceptibility was assessed using broth dilution and viable cell counting. Using whole-genome sequencing, we explored the potential roles of outer membrane protein P2 (OmpP2), LytM proteins and the dcw gene cluster in imipenem heteroresistance.

RESULTS

All 46 imipenem-heteroresistant isolates (IMI^hR) harboured amino acid substitutions in the ftsI gene, which encodes PBP3, corresponding to 25 different mutation patterns that varied from the ftsI gene mutation patterns found in imipenem-susceptible isolates. Among all PBPs, the highest affinity to imipenem was documented for PBP3 (IC₅₀, 0.004 μg/mL). Different amino acid substitutions and insertions were noted in OmpP2, suggesting a relationship with imipenem heteroresistance. The IMI^hR isolates were affected by CCCP differently and displayed a higher percentage of killing by imipenem in CCCP-treated cells at concentrations ranging between 0.5 and 8 μg/mL.

CONCLUSIONS

The present study provides robust evidence indicating that in combination with the altered PBP3, the slowed drug influx and its enhanced efflux due to the loss of regulation led to the development of imipenem heteroresistance in NTHi.

Critical Timing without a Timer for Embryonic Development

Timing of embryonic development is precisely controlled, but the mechanisms underlying biological timers are still unclear. Here, a validated model for timing under control of Sonic Hedgehog is revisited and generalized to an arbitrary number of genes. The developmental dynamics where a temporal sequence of gene expression recapitulates a steady-state spatial pattern can be realized through a simple network close to criticality, controlled by the duration of exposure to a morphogen. Criticality simultaneously accounts for many observed biological properties, such as timing, multistability, and canalization of genetic expression. This process can be parsimoniously generalized in many dimensions with a minimum number of genes, all repressing each other with asymmetrical strengths, which also explains sequential activation of different fates. Separation of timescales allows for a simple analytical interpretation. Finally, it is shown that even in the presence of noise, coupling between cells preserves criticality and robust patterning. The model offers a simple theoretical framework for the study of emergent developmental timers.

Positional information from oscillatory phase shifts : insights from in silico evolution

Complex cellular decisions are based on temporal dynamics of pathways, including genetic oscillators. In development, recent works on vertebrae formation have suggested that relative phase of genetic oscillators encode positional information, including differentiation front defining vertebrae positions. Precise mechanisms for this are still unknown. Here, we use computational evolution to find gene network topologies that can compute the phase difference between oscillators and convert it into a decoder morphogen concentration. Two types of networks are discovered, based on symmetry properties of the decoder gene. So called asymmetric networks are studied, and two submodules are identified converting phase information into an amplitude variable. Those networks naturally display a 'shock' for a well defined phase difference, that can be used to define a wavefront of differentiation. We show how implementation of these ideas reproduce experimental features of vertebrate segmentation.

Predicting Ancestral Segmentation Phenotypes from Drosophila to Anopheles Using In Silico Evolution

Molecular evolution is an established technique for inferring gene homology but regulatory DNA turns over so rapidly that inference of ancestral networks is often impossible. In silico evolution is used to compute the most parsimonious path in regulatory space for anterior-posterior patterning linking two Dipterian species. The expression pattern of gap genes has evolved between Drosophila (fly) and Anopheles (mosquito), yet one of their targets, eve, has remained invariant. Our model predicts that stripe 5 in fly disappears and a new posterior stripe is created in mosquito, thus eve stripe modules 3+7 and 4+6 in fly are homologous to 3+6 and 4+5 in mosquito. We can place Clogmia on this evolutionary pathway and it shares the mosquito homologies. To account for the evolution of the other pair-rule genes in the posterior we have to assume that the ancestral Dipterian utilized a dynamic method to phase those genes in relation to eve.

The last common ancestor of the fruit fly (Drosophila) and mosquito (Anopheles) lived more than 200 Million years ago. Can we use available data on insects alive today to infer what their ancestor looked like? In this manuscript, we focus on early embryonic development, when stripes of genetic expression appear and define the location of insect segments (“segmentation”). We use an evolutionary algorithm to reconstruct and predict dynamics of genes controlling stripes in the last common ancestor of fly and mosquito. We predict a new and different combinatorial logic of stripe formation in mosquito compared to fly, which is fully consistent with development of intermediate species such as moth-fly (Clogmia). Our simulations further suggest that the dynamics of gene expression in this last common ancestor were similar to other insects, such as wasps (Nasonia). Our method illustrates how computational methods inspired by machine learning and non-linear physics can be used to infer gene dynamics in species that disappeared millions of years ago.

An empirical model for calculation of the collimator contamination dose in therapeutic proton beams

Collimators are used as lateral beam shaping devices in proton therapy with passive scattering beam lines. The dose contamination due to collimator scattering can be as high as 10% of the maximum dose and influences calculation of the output factor or monitor units (MU). To date, commercial treatment planning systems generally use a zero-thickness collimator approximation ignoring edge scattering in the aperture collimator and few analytical models have been proposed to take scattering effects into account, mainly limited to the inner collimator face component. The aim of this study was to characterize and model aperture contamination by means of a fast and accurate analytical model. The entrance face collimator scatter distribution was modeled as a 3D secondary dose source. Predicted dose contaminations were compared to measurements and Monte Carlo simulations. Measurements were performed on two different proton beam lines (a fixed horizontal beam line and a gantry beam line) with divergent apertures and for several field sizes and energies. Discrepancies between analytical algorithm dose prediction and measurements were decreased from 10% to 2% using the proposed model. Gamma-index (2%/1 mm) was respected for more than 90% of pixels. The proposed analytical algorithm increases the accuracy of analytical dose calculations with reasonable computation times.

Randomized, placebo-controlled, double-blind clinical trial to evaluate the efficacy of polyhexanide for topical decolonization of MRSA carriers

OBJECTIVES

The objective of this study was to evaluate the efficacy of polyhexanide (Prontoderm(®)) in eliminating MRSA carriage.

METHODS

In a 1900 bed teaching hospital, MRSA-colonized patients were randomized into a double-blind, placebo-controlled superiority trial between January 2011 and July 2014. Patients were treated with either polyhexanide or placebo applied to the anterior nares (thrice daily) and skin (once daily) for 10 days. The primary outcome was MRSA decolonization at day 28 (D28) after the end of treatment assessed by ITT responder and PP analyses (microbiological follow-up ± 7 days and topical treatment ≥ 5 days). Secondary outcomes included safety, emergence of resistance and MRSA genotype changes. Registered trial number ISRCTN02288276.

RESULTS

Of 2590 patients screened, 146 (polyhexanide group, 71; placebo group, 75) were included. ITT analysis showed that 24/71 (33.8%) patients in the polyhexanide group versus 22/75 (29.3%) in the placebo group were MRSA-free at D28 (risk difference, 4.5%; 95% CI, -10.6% to 19.5%; P = 0.56). PP analysis confirmed the results with 19/53 (35.8%) decolonized polyhexanide-treated patients versus 17/56 (30.4%) in the placebo arm (risk difference, 5.5%; 95% CI, -12.2% to 23%; P = 0.54). Nine serious adverse events occurred in the polyhexanide group versus 12 in the placebo group; none was attributable to study medication. Emergence of polyhexanide resistance or cross-resistance between polyhexanide and chlorhexidine was not observed. No case of exogenous recolonization by a genotypically different MRSA strain was documented.

CONCLUSIONS

This study suggests that under real-life conditions, a single polyhexanide decolonization course is not effective in eradicating MRSA carriage.

Dynamics of the slowing segmentation clock reveal alternating two-segment periodicity

The formation of reiterated somites along the vertebrate body axis is controlled by the segmentation clock, a molecular oscillator expressed within presomitic mesoderm (PSM) cells. Although PSM cells oscillate autonomously, they coordinate with neighboring cells to generate a sweeping wave of cyclic gene expression through the PSM that has a periodicity equal to that of somite formation. The velocity of each wave slows as it moves anteriorly through the PSM, although the dynamics of clock slowing have not been well characterized. Here, we investigate segmentation clock dynamics in the anterior PSM in developing zebrafish embryos using an in vivo clock reporter, her1:her1-venus. The her1:her1-venus reporter has single-cell resolution, allowing us to follow segmentation clock oscillations in individual cells in real-time. By retrospectively tracking oscillations of future somite boundary cells, we find that clock reporter signal increases in anterior PSM cells and that the periodicity of reporter oscillations slows to about ∼1.5 times the periodicity in posterior PSM cells. This gradual slowing of the clock in the anterior PSM creates peaks of clock expression that are separated at a two-segment periodicity both spatially and temporally, a phenomenon we observe in single cells and in tissue-wide analyses. These results differ from previous predictions that clock oscillations stop or are stabilized in the anterior PSM. Instead, PSM cells oscillate until they incorporate into somites. Our findings suggest that the segmentation clock may signal somite formation using a phase gradient with a two-somite periodicity.