Extrinsic noise-induced transitions to bimodal dynamics have been largely investigated in a variety of chemical, physical, and biological systems. In the standard approach in physical and chemical ...systems, the key properties that make these systems mathematically tractable are that the noise appears linearly in the dynamical equations, and it is assumed Gaussian and white. In biology, the Gaussian approximation has been successful in specific systems, but the relevant noise being usually non-Gaussian, non-white, and nonlinear poses serious limitations to its general applicability. Here we revisit the fundamental features of linear Gaussian noise, pinpoint its limitations, and review recent new approaches based on nonlinear bounded noises, which highlight novel mechanisms to account for transitions to bimodal behaviour. We do this by considering a simple but fundamental gene expression model, the repressed gene, which is characterized by linear and nonlinear dependencies on external parameters. We then review a general methodology introduced recently, so-called nonlinear noise filtering, which allows the investigation of linear, nonlinear, Gaussian and non-Gaussian noises. We also present a derivation of it, which highlights its dynamical origin. Testing the methodology on the repressed gene confirms that the emergence of noise-induced transitions appears to be strongly dependent on the type of noise adopted, and on the degree of nonlinearity present in the system.
Multipotent neural crest (NC) progenitors generate an astonishing array of derivatives, including neuronal, skeletal components and pigment cells (chromatophores), but the molecular mechanisms ...allowing balanced selection of each fate remain unknown. In zebrafish, melanocytes, iridophores and xanthophores, the three chromatophore lineages, are thought to share progenitors and so lend themselves to investigating the complex gene regulatory networks (GRNs) underlying fate segregation of NC progenitors. Although the core GRN governing melanocyte specification has been previously established, those guiding iridophore and xanthophore development remain elusive. Here we focus on the iridophore GRN, where mutant phenotypes identify the transcription factors Sox10, Tfec and Mitfa and the receptor tyrosine kinase, Ltk, as key players. Here we present expression data, as well as loss and gain of function results, guiding the derivation of an initial iridophore specification GRN. Moreover, we use an iterative process of mathematical modelling, supplemented with a Monte Carlo screening algorithm suited to the qualitative nature of the experimental data, to allow for rigorous predictive exploration of the GRN dynamics. Predictions were experimentally evaluated and testable hypotheses were derived to construct an improved version of the GRN, which we showed produced outputs consistent with experimentally observed gene expression dynamics. Our study reveals multiple important regulatory features, notably a sox10-dependent positive feedback loop between tfec and ltk driving iridophore specification; the molecular basis of sox10 maintenance throughout iridophore development; and the cooperation between sox10 and tfec in driving expression of pnp4a, a key differentiation gene. We also assess a candidate repressor of mitfa, a melanocyte-specific target of sox10. Surprisingly, our data challenge the reported role of Foxd3, an established mitfa repressor, in iridophore regulation. Our study builds upon our previous systems biology approach, by incorporating physiologically-relevant parameter values and rigorous evaluation of parameter values within a qualitative data framework, to establish for the first time the core GRN guiding specification of the iridophore lineage.
The COVID-19 pandemic induced an extraordinary impact on public mental health to a degree not completely understood, especially in vulnerable populations such as breast cancer (BC) survivors. In this ...study, we described the short- (after 3-month) and long- (after 12-month) term effects of a multidisciplinary home-based lifestyle intervention in Italian women BC survivors during the first year of COVID-19 pandemic.
In total, 30 Italian BC survivors with risk factors for recurrence took part in the ongoing MoviS trial (protocol: NCT04818359). Between January 2020 and January 2021, a 3-month lifestyle intervention based on psychological counseling, nutrition, and exercise was carried out. Participants were asked to fill out psychological questionnaires for the assessment of quality of life (QoL) indicators (European Organization for Research and Treatment of Cancer QoL, EORTC-QLQ-C30) and psychological health measures such as fatigue (Brief Fatigue Inventory, BFI), distress (Distress Thermometer, DT and Psychological Distress Inventory, PDI), cancer-related fatigue (Verbal Rating Scale, VRS), and mood states (Profile of Mood States Questionnaire, POMS). IBM SPSS Statistical Software version 27.0 and R Project for Statistical Computing version 4.2.1 were used to process data. All participants were assessed at four time points: T0 (baseline), T1 (3-month), and follow-up at T2 and T3 (6- and 12-month, respectively) to measure primary (quality of life indicators) and secondary (psychological health) outcomes. Friedman non parametric test and Wilcoxon signed rank test (with Bonferroni correction) were conducted to investigate the statistically significant differences in psychometric scores and between assessment times.
Compared to baseline (T0), at T1 most of the QoL indicators (i.e., symptoms of fatigue and general health) were improved (p < 0.017) with the exception of a worsening in participants' social functioning ability. Also, perception of severity of fatigue, distress, cancer-related fatigue, depression, and anger enhanced. Compared to baseline (T0), at T3 we mainly observed a stable condition with T0-T1 pairwise comparison, however other secondary outcomes (i.e., fatigue mood state, confusion, and anxiety) significantly improved.
Our preliminary findings support the proposal of this lifestyle intervention for BC survivors. Despite the home-confinement due to the COVID-19 pandemic, the intervention surprisingly improved QoL indicators and psychological health of the participants.
Neural crest cells are highly multipotent stem cells, but it remains unclear how their fate restriction to specific fates occurs. The direct fate restriction model hypothesises that migrating cells ...maintain full multipotency, whilst progressive fate restriction envisages fully multipotent cells transitioning to partially-restricted intermediates before committing to individual fates. Using zebrafish pigment cell development as a model, we show applying NanoString hybridization single cell transcriptional profiling and RNAscope in situ hybridization that neural crest cells retain broad multipotency throughout migration and even in post-migratory cells in vivo, with no evidence for partially-restricted intermediates. We find that leukocyte tyrosine kinase early expression marks a multipotent stage, with signalling driving iridophore differentiation through repression of fate-specific transcription factors for other fates. We reconcile the direct and progressive fate restriction models by proposing that pigment cell development occurs directly, but dynamically, from a highly multipotent state, consistent with our recently-proposed Cyclical Fate Restriction model.
Sequence alignment underpins all of comparative genomics, yet it remains an incompletely solved problem. In particular, the statistical uncertainty within inferred alignments is often disregarded, ...while parametric or phylogenetic inferences are considered meaningless without confidence estimates. Here, we report on a theoretical and simulation study of pairwise alignments of genomic DNA at human-mouse divergence. We find that >15% of aligned bases are incorrect in existing whole-genome alignments, and we identify three types of alignment error, each leading to systematic biases in all algorithms considered. Careful modeling of the evolutionary process improves alignment quality; however, these improvements are modest compared with the remaining alignment errors, even with exact knowledge of the evolutionary model, emphasizing the need for statistical approaches to account for uncertainty. We develop a new algorithm, Marginalized Posterior Decoding (MPD), which explicitly accounts for uncertainties, is less biased and more accurate than other algorithms we consider, and reduces the proportion of misaligned bases by a third compared with the best existing algorithm. To our knowledge, this is the first nonheuristic algorithm for DNA sequence alignment to show robust improvements over the classic Needleman-Wunsch algorithm. Despite this, considerable uncertainty remains even in the improved alignments. We conclude that a probabilistic treatment is essential, both to improve alignment quality and to quantify the remaining uncertainty. This is becoming increasingly relevant with the growing appreciation of the importance of noncoding DNA, whose study relies heavily on alignments. Alignment errors are inevitable, and should be considered when drawing conclusions from alignments. Software and alignments to assist researchers in doing this are provided at http://genserv.anat.ox.ac.uk/grape/.
An understanding of the dynamics of the metabolic profile of a bacterial cell is sought from a dynamical systems analysis of kinetic models. This modelling formalism relies on a deterministic ...mathematical description of enzyme kinetics and their metabolite regulation. However, it is severely impeded by the lack of available kinetic information, limiting the size of the system that can be modelled. Furthermore, the subsystem of the metabolic network whose dynamics can be modelled is faced with three problems: how to parameterize the model with mostly incomplete steady state data, how to close what is now an inherently open system, and how to account for the impact on growth. In this study we address these challenges of kinetic modelling by capitalizing on multi-'omics' steady state data and a genome-scale metabolic network model. We use these to generate parameters that integrate knowledge embedded in the genome-scale metabolic network model, into the most comprehensive kinetic model of the central carbon metabolism of E. coli realized to date. As an application, we performed a dynamical systems analysis of the resulting enriched model. This revealed bistability of the central carbon metabolism and thus its potential to express two distinct metabolic states. Furthermore, since our model-informing technique ensures both stable states are constrained by the same thermodynamically feasible steady state growth rate, the ensuing bistability represents a temporal coexistence of the two states, and by extension, reveals the emergence of a phenotypically heterogeneous population.
The mechanisms generating stably differentiated cell-types from multipotent precursors are key to understanding normal development and have implications for treatment of cancer and the therapeutic ...use of stem cells. Pigment cells are a major derivative of neural crest stem cells and a key model cell-type for our understanding of the genetics of cell differentiation. Several factors driving melanocyte fate specification have been identified, including the transcription factor and master regulator of melanocyte development, Mitf, and Wnt signalling and the multipotency and fate specification factor, Sox10, which drive mitf expression. While these factors together drive multipotent neural crest cells to become specified melanoblasts, the mechanisms stabilising melanocyte differentiation remain unclear. Furthermore, there is controversy over whether Sox10 has an ongoing role in melanocyte differentiation. Here we use zebrafish to explore in vivo the gene regulatory network (GRN) underlying melanocyte specification and differentiation. We use an iterative process of mathematical modelling and experimental observation to explore methodically the core melanocyte GRN we have defined. We show that Sox10 is not required for ongoing differentiation and expression is downregulated in differentiating cells, in response to Mitfa and Hdac1. Unexpectedly, we find that Sox10 represses Mitf-dependent expression of melanocyte differentiation genes. Our systems biology approach allowed us to predict two novel features of the melanocyte GRN, which we then validate experimentally. Specifically, we show that maintenance of mitfa expression is Mitfa-dependent, and identify Sox9b as providing an Mitfa-independent input to melanocyte differentiation. Our data supports our previous suggestion that Sox10 only functions transiently in regulation of mitfa and cannot be responsible for long-term maintenance of mitfa expression; indeed, Sox10 is likely to slow melanocyte differentiation in the zebrafish embryo. More generally, this novel approach to understanding melanocyte differentiation provides a basis for systematic modelling of differentiation in this and other cell-types.
One of the most challenging problems in microbiology is to understand how a small fraction of microbes that resists killing by antibiotics can emerge in a population of genetically identical cells, ...the phenomenon known as persistence or drug tolerance. Its characteristic signature is the biphasic kill curve, whereby microbes exposed to a bactericidal agent are initially killed very rapidly but then much more slowly. Here we relate this problem to the more general problem of understanding the emergence of distinct growth phenotypes in clonal populations. We address the problem mathematically by adopting the framework of the phenomenon of so-called weak ergodicity breaking, well known in dynamical physical systems, which we extend to the biological context. We show analytically and by direct stochastic simulations that distinct growth phenotypes can emerge as a consequence of slow-down of stochastic fluctuations in the expression of a gene controlling growth rate. In the regime of fast gene transcription, the system is ergodic, the growth rate distribution is unimodal, and accounts for one phenotype only. In contrast, at slow transcription and fast translation, weakly non-ergodic components emerge, the population distribution of growth rates becomes bimodal, and two distinct growth phenotypes are identified. When coupled to the well-established growth rate dependence of antibiotic killing, this model describes the observed fast and slow killing phases, and reproduces much of the phenomenology of bacterial persistence. The model has major implications for efforts to develop control strategies for persistent infections.
Breast cancer (BC) is the most common invasive cancer in women, and exercise can significantly improve the outcomes of BC survivors. MoviS (Movement and Health Beyond Care) is a randomized controlled ...trial aimed to evaluate the potential health benefits of exercise and proper nutritional habits. This study aims to assess the efficacy of aerobic exercise training in improving quality of life (QoL) and health-related factors in high-risk BC.
One hundred seventy-two BC survivor women, aged 30-70 years, non-metastatic, stage 0-III, non-physically active, 6-12 months post-surgery, and post chemo- or radiotherapy, will be recruited in this study. Women will be randomly allocated to the intervention arm (lifestyle recommendations and MoviS Training) or control arm (lifestyle recommendations). The MoviS training consists of 12 weeks of aerobic exercise training (2 days/week of supervised and 1 day/week of unsupervised exercise) with a progressive increase in exercise intensity (40-70% of heart rate reserve) and duration (20-60 min). Both arms will receive counseling on healthy lifestyle habits (nutrition and exercise) based on the World Cancer Research Fund International (WCRF) 2018 guidelines. The primary outcome is the improvement of the QoL. The secondary outcomes are improvement of health-related parameters such as Mediterranean diet adherence, physical activity level, flexibility, muscular fitness, fatigue, cardiorespiratory fitness (estimated maximal oxygen uptake), echocardiographic parameters, heart rate variability (average of the standard deviations of all 5 min normal to normal intervals (ASDNN/5 min) and 24 h very low and low frequency), and metabolic, endocrine, and inflammatory serum biomarkers (glycemia, insulin resistance, progesterone, testosterone, and high-sensitivity C-reactive protein).
This trial aims to evaluate if supervised exercise may improve QoL and health-related factors of BC survivors with a high risk of recurrence. Findings from this project could provide knowledge improvement in the field of exercise oncology through the participation of a multidisciplinary team that will provide a coordinated program of cancer care to improve healthcare quality, improve prognosis, increase survival times and QoL, and reduce the risk of BC recurrence.
ClinicalTrials.gov NCT04818359 . Retrospectively registered on March 26, 2021.
Post-transcriptional gene regulation is driven by RNA-binding proteins (RBPs). Recent global approaches suggest widespread autoregulation of RBPs through binding to their own mRNA; however, little is ...known about the regulatory impact and quantitative models remain elusive. By integration of several independent kinetic parameters and abundance data, we modelled autoregulatory feedback loops for six canonical and non-canonical RBPs from the yeast
, namely Hrb1p, Hek2/Khd1p, Ski2p, Npl3p, Pfk2p, and Map1p. By numerically solving ordinary differential equations, we compared non-feedback models with models that considered the RPBs as post-transcriptional activators/repressors of their own expression. While our results highlight a substantial gap between predicted protein output and experimentally determined protein abundances applying a no-feedback model, addition of positive feedback loops are surprisingly versatile and can improve predictions towards experimentally determined protein levels, whereas negative feedbacks are particularly sensitive to cooperativity. Our data suggests that introduction of feedback loops supported by real data can improve models of post-transcriptional gene expression.