The parasite intraerythrocytic cycle and human circadian cycle are coupled during malaria infection.

During infections with the malaria parasites Plasmodium vivax, patients exhibit rhythmic fevers every 48 h. These fever cycles correspond with the time the parasites take to traverse the intraerythrocytic cycle (IEC). In other Plasmodium species that infect either humans or mice, the IEC is likely guided by a parasite-intrinsic clock [Rijo-Ferreiraet al., Science 368, 746-753 (2020); Smith et al., Science 368, 754-759 (2020)], suggesting that intrinsic clock mechanisms may be a fundamental feature of malaria parasites. Moreover, because Plasmodium cycle times are multiples of 24 h, the IECs may be coordinated with the host circadian clock(s). Such coordination could explain the synchronization of the parasite population in the host and enable alignment of IEC and circadian cycle phases. We utilized an ex vivo culture of whole blood from patients infected with P. vivax to examine the dynamics of the host circadian transcriptome and the parasite IEC transcriptome. Transcriptome dynamics revealed that the phases of the host circadian cycle and the parasite IEC are correlated across multiple patients, showing that the cycles are phase coupled. In mouse model systems, host-parasite cycle coupling appears to provide a selective advantage for the parasite. Thus, understanding how host and parasite cycles are coupled in humans could enable antimalarial therapies that disrupt this coupling.

Experimental guidance for discovering genetic networks through hypothesis reduction on time series.

Large programs of dynamic gene expression, like cell cyles and circadian rhythms, are controlled by a relatively small "core" network of transcription factors and post-translational modifiers, working in concerted mutual regulation. Recent work suggests that system-independent, quantitative features of the dynamics of gene expression can be used to identify core regulators. We introduce an approach of iterative network hypothesis reduction from time-series data in which increasingly complex features of the dynamic expression of individual, pairs, and entire collections of genes are used to infer functional network models that can produce the observed transcriptional program. The culmination of our work is a computational pipeline, Iterative Network Hypothesis Reduction from Temporal Dynamics (Inherent dynamics pipeline), that provides a priority listing of targets for genetic perturbation to experimentally infer network structure. We demonstrate the capability of this integrated computational pipeline on synthetic and yeast cell-cycle data.

Alignment of Synchronized Time-Series Data Using the Characterizing Loss of Cell Cycle Synchrony Model for Cross-Experiment Comparisons.

Investigating the cell cycle often depends on synchronizing cell populations to measure various parameters in a time series as the cells traverse the cell cycle. However, even under similar conditions, replicate experiments display differences in the time required to recover from synchrony and to traverse the cell cycle, thus preventing direct comparisons at each time point. The problem of comparing dynamic measurements across experiments is exacerbated in mutant populations or in alternative growth conditions that affect the synchrony recovery time and/or the cell-cycle period. We have previously published a parametric mathematical model named Characterizing Loss of Cell Cycle Synchrony (CLOCCS) that monitors how synchronous populations of cells release from synchrony and progress through the cell cycle. The learned parameters from the model can then be used to convert experimental time points from synchronized time-series experiments into a normalized time scale (lifeline points). Rather than representing the elapsed time in minutes from the start of the experiment, the lifeline scale represents the progression from synchrony to cell-cycle entry and then through the phases of the cell cycle. Since lifeline points correspond to the phase of the average cell within the synchronized population, this normalized time scale allows for direct comparisons between experiments, including those with varying periods and recovery times. Furthermore, the model has been used to align cell-cycle experiments between different species (e.g., Saccharomyces cerevisiae and Schizosaccharomyces pombe), thus enabling direct comparison of cell-cycle measurements, which may reveal evolutionary similarities and differences.

Inherent Dynamics Visualizer, an Interactive Application for Evaluating and Visualizing Outputs from a Gene Regulatory Network Inference Pipeline.

Developing gene regulatory network models is a major challenge in systems biology. Several computational tools and pipelines have been developed to tackle this challenge, including the newly developed Inherent Dynamics Pipeline. The Inherent Dynamics Pipeline consists of several previously published tools that work synergistically and are connected in a linear fashion, where the output of one tool is then used as input for the following tool. As with most computational techniques, each step of the Inherent Dynamics Pipeline requires the user to make choices about parameters that don't have a precise biological definition. These choices can substantially impact gene regulatory network models produced by the analysis. For this reason, the ability to visualize and explore the consequences of various parameter choices at each step can help increase confidence in the choices and the results.The Inherent Dynamics Visualizer is a comprehensive visualization package that streamlines the process of evaluating parameter choices through an interactive interface within a web browser. The user can separately examine the output of each step of the pipeline, make intuitive changes based on visual information, and benefit from the automatic production of necessary input files for the Inherent Dynamics Pipeline. The Inherent Dynamics Visualizer provides an unparalleled level of access to a highly intricate tool for the discovery of gene regulatory networks from time series transcriptomic data.