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.

Assessment of Simulated Surveillance Testing and Quarantine in a SARS-CoV-2-Vaccinated Population of Students on a University Campus.

Importance: The importance of surveillance testing and quarantine on university campuses to limit SARS-CoV-2 transmission needs to be reevaluated in the context of a complex and rapidly changing environment that includes vaccines, variants, and waning immunity. Also, recent US Centers for Disease Control and Prevention guidelines suggest that vaccinated students do not need to participate in surveillance testing. Objective: To evaluate the use of surveillance testing and quarantine in a fully vaccinated student population for whom vaccine effectiveness may be affected by the type of vaccination, presence of variants, and loss of vaccine-induced or natural immunity over time. Design Setting and Participants: In this simulation study, an agent-based Susceptible, Exposed, Infected, Recovered model was developed with some parameters estimated using data from the 2020 to 2021 academic year at Duke University (Durham, North Carolina) that described a simulated population of 5000 undergraduate students residing on campus in residential dormitories. This study assumed that 100% of residential undergraduates are vaccinated. Under varying levels of vaccine effectiveness (90%, 75%, and 50%), the reductions in the numbers of positive cases under various mitigation strategies that involved surveillance testing and quarantine were estimated. Main Outcomes and Measures: The percentage of students infected with SARS-CoV-2 each day for the course of the semester (100 days) and the total number of isolated or quarantined students were estimated. Results: A total of 5000 undergraduates were simulated in the study. In simulations with 90% vaccine effectiveness, weekly surveillance testing was associated with only marginally reduced viral transmission. At 50% to 75% effectiveness, surveillance testing was estimated to reduce the number of infections by as much as 93.6%. A 10-day quarantine protocol for exposures was associated with only modest reduction in infections until vaccine effectiveness dropped to 50%. Increased testing of reported contacts was estimated to be at least as effective as quarantine at limiting infections. Conclusions and Relevance: In this simulated modeling study of infection dynamics on a college campus where 100% of the student body is vaccinated, weekly surveillance testing was associated with a substantial reduction of campus infections with even a modest loss of vaccine effectiveness. Model simulations also suggested that an increased testing cadence can be as effective as a 10-day quarantine period at limiting infections. Together, these findings provide a potential foundation for universities to design appropriate mitigation protocols for the 2021 to 2022 academic year.

Implementation of a Pooled Surveillance Testing Program for Asymptomatic SARS-CoV-2 Infections on a College Campus - Duke University, Durham, North Carolina, August 2-October 11, 2020.

On university campuses and in similar congregate environments, surveillance testing of asymptomatic persons is a critical strategy (1,2) for preventing transmission of SARS-CoV-2, the virus that causes coronavirus disease 2019 (COVID-19). All students at Duke University, a private research university in Durham, North Carolina, signed the Duke Compact (3), agreeing to observe mandatory masking, social distancing, and participation in entry and surveillance testing. The university implemented a five-to-one pooled testing program for SARS-CoV-2 using a quantitative, in-house, laboratory-developed, real-time reverse transcription-polymerase chain reaction (RT-PCR) test (4,5). Pooling of specimens to enable large-scale testing while minimizing use of reagents was pioneered during the human immunodeficiency virus pandemic (6). A similar methodology was adapted for Duke University's asymptomatic testing program. The baseline SARS-CoV-2 testing plan was to distribute tests geospatially and temporally across on- and off-campus student populations. By September 20, 2020, asymptomatic testing was scaled up to testing targets, which include testing for residential undergraduates twice weekly, off-campus undergraduates one to two times per week, and graduate students approximately once weekly. In addition, in response to newly identified positive test results, testing was focused in locations or within cohorts where data suggested an increased risk for transmission. Scale-up over 4 weeks entailed redeploying staff members to prepare 15 campus testing sites for specimen collection, developing information management tools, and repurposing laboratory automation to establish an asymptomatic surveillance system. During August 2-October 11, 68,913 specimens from 10,265 graduate and undergraduate students were tested. Eighty-four specimens were positive for SARS-CoV-2, and 51% were among persons with no symptoms. Testing as a result of contact tracing identified 27.4% of infections. A combination of risk-reduction strategies and frequent surveillance testing likely contributed to a prolonged period of low transmission on campus. These findings highlight the importance of combined testing and contact tracing strategies beyond symptomatic testing, in association with other preventive measures. Pooled testing balances resource availability with supply-chain disruptions, high throughput with high sensitivity, and rapid turnaround with an acceptable workload.

An intrinsic oscillator drives the blood stage cycle of the malaria parasite Plasmodium falciparum.

The blood stage of the infection of the malaria parasite Plasmodium falciparum exhibits a 48-hour developmental cycle that culminates in the synchronous release of parasites from red blood cells, which triggers 48-hour fever cycles in the host. This cycle could be driven extrinsically by host circadian processes or by a parasite-intrinsic oscillator. To distinguish between these hypotheses, we examine the P. falciparum cycle in an in vitro culture system and show that the parasite has molecular signatures associated with circadian and cell cycle oscillators. Each of the four strains examined has a different period, which indicates strain-intrinsic period control. Finally, we demonstrate that parasites have low cell-to-cell variance in cycle period, on par with a circadian oscillator. We conclude that an intrinsic oscillator maintains Plasmodium's rhythmic life cycle.

Guidelines for Genome-Scale Analysis of Biological Rhythms.

Genome biology approaches have made enormous contributions to our understanding of biological rhythms, particularly in identifying outputs of the clock, including RNAs, proteins, and metabolites, whose abundance oscillates throughout the day. These methods hold significant promise for future discovery, particularly when combined with computational modeling. However, genome-scale experiments are costly and laborious, yielding "big data" that are conceptually and statistically difficult to analyze. There is no obvious consensus regarding design or analysis. Here we discuss the relevant technical considerations to generate reproducible, statistically sound, and broadly useful genome-scale data. Rather than suggest a set of rigid rules, we aim to codify principles by which investigators, reviewers, and readers of the primary literature can evaluate the suitability of different experimental designs for measuring different aspects of biological rhythms. We introduce CircaInSilico, a web-based application for generating synthetic genome biology data to benchmark statistical methods for studying biological rhythms. Finally, we discuss several unmet analytical needs, including applications to clinical medicine, and suggest productive avenues to address them.

The Local Edge Machine: inference of dynamic models of gene regulation.

We present a novel approach, the Local Edge Machine, for the inference of regulatory interactions directly from time-series gene expression data. We demonstrate its performance, robustness, and scalability on in silico datasets with varying behaviors, sizes, and degrees of complexity. Moreover, we demonstrate its ability to incorporate biological prior information and make informative predictions on a well-characterized in vivo system using data from budding yeast that have been synchronized in the cell cycle. Finally, we use an atlas of transcription data in a mammalian circadian system to illustrate how the method can be used for discovery in the context of large complex networks.

SW1PerS: Sliding windows and 1-persistence scoring; discovering periodicity in gene expression time series data.

BACKGROUND: Identifying periodically expressed genes across different processes (e.g. the cell and metabolic cycles, circadian rhythms, etc) is a central problem in computational biology. Biological time series may contain (multiple) unknown signal shapes of systemic relevance, imperfections like noise, damping, and trending, or limited sampling density. While there exist methods for detecting periodicity, their design biases (e.g. toward a specific signal shape) can limit their applicability in one or more of these situations. METHODS: We present in this paper a novel method, SW1PerS, for quantifying periodicity in time series in a shape-agnostic manner and with resistance to damping. The measurement is performed directly, without presupposing a particular pattern, by evaluating the circularity of a high-dimensional representation of the signal. SW1PerS is compared to other algorithms using synthetic data and performance is quantified under varying noise models, noise levels, sampling densities, and signal shapes. Results on biological data are also analyzed and compared. RESULTS: On the task of periodic/not-periodic classification, using synthetic data, SW1PerS outperforms all other algorithms in the low-noise regime. SW1PerS is shown to be the most shape-agnostic of the evaluated methods, and the only one to consistently classify damped signals as highly periodic. On biological data, and for several experiments, the lists of top 10% genes ranked with SW1PerS recover up to 67% of those generated with other popular algorithms. Moreover, the list of genes from data on the Yeast metabolic cycle which are highly-ranked only by SW1PerS, contains evidently non-cosine patterns (e.g. ECM33, CDC9, SAM1,2 and MSH6) with highly periodic expression profiles. In data from the Yeast cell cycle SW1PerS identifies genes not preferred by other algorithms, hence not previously reported as periodic, but found in other experiments such as the universal growth rate response of Slavov. These genes are BOP3, CDC10, YIL108W, YER034W, MLP1, PAC2 and RTT101. CONCLUSIONS: In biological systems with low noise, i.e. where periodic signals with interesting shapes are more likely to occur, SW1PerS can be used as a powerful tool in exploratory analyses. Indeed, by having an initial set of periodic genes with a rich variety of signal types, pattern/shape information can be included in the study of systems and the generation of hypotheses regarding the structure of gene regulatory networks.

Easily repairable networks: reconnecting nodes after damage.

We introduce a simple class of distribution networks that withstand damage by being repairable instead of redundant. Instead of asking how hard it is to disconnect nodes through damage, we ask how easy it is to reconnect nodes after damage. We prove that optimal networks on regular lattices have an expected cost of reconnection proportional to the lattice length, and that such networks have exactly three levels of structural hierarchy. We extend our results to networks subject to repeated attacks, in which the repairs themselves must be repairable. We find that, in exchange for a modest increase in repair cost, such networks are able to withstand any number of attacks.

Checkpoints couple transcription network oscillator dynamics to cell-cycle progression.

BACKGROUND: The coupling of cyclin dependent kinases (CDKs) to an intrinsically oscillating network of transcription factors has been proposed to control progression through the cell cycle in budding yeast, Saccharomyces cerevisiae. The transcription network regulates the temporal expression of many genes, including cyclins, and drives cell-cycle progression, in part, by generating successive waves of distinct CDK activities that trigger the ordered program of cell-cycle events. Network oscillations continue autonomously in mutant cells arrested by depletion of CDK activities, suggesting the oscillator can be uncoupled from cell-cycle progression. It is not clear what mechanisms, if any, ensure that the network oscillator is restrained when progression in normal cells is delayed or arrested. A recent proposal suggests CDK acts as a master regulator of cell-cycle processes that have the potential for autonomous oscillatory behavior. RESULTS: Here we find that mitotic CDK is not sufficient for fully inhibiting transcript oscillations in arrested cells. We do find that activation of the DNA replication and spindle assembly checkpoints can fully arrest the network oscillator via overlapping but distinct mechanisms. Further, we demonstrate that the DNA replication checkpoint effector protein, Rad53, acts to arrest a portion of transcript oscillations in addition to its role in halting cell-cycle progression. CONCLUSIONS: Our findings indicate that checkpoint mechanisms, likely via phosphorylation of network transcription factors, maintain coupling of the network oscillator to progression during cell-cycle arrest.

Design and analysis of large-scale biological rhythm studies: a comparison of algorithms for detecting periodic signals in biological data.

MOTIVATION: To discover and study periodic processes in biological systems, we sought to identify periodic patterns in their gene expression data. We surveyed a large number of available methods for identifying periodicity in time series data and chose representatives of different mathematical perspectives that performed well on both synthetic data and biological data. Synthetic data were used to evaluate how each algorithm responds to different curve shapes, periods, phase shifts, noise levels and sampling rates. The biological datasets we tested represent a variety of periodic processes from different organisms, including the cell cycle and metabolic cycle in Saccharomyces cerevisiae, circadian rhythms in Mus musculus and the root clock in Arabidopsis thaliana. RESULTS: From these results, we discovered that each algorithm had different strengths. Based on our findings, we make recommendations for selecting and applying these methods depending on the nature of the data and the periodic patterns of interest. Additionally, these results can also be used to inform the design of large-scale biological rhythm experiments so that the resulting data can be used with these algorithms to detect periodic signals more effectively.

3D phenotyping and quantitative trait locus mapping identify core regions of the rice genome controlling root architecture.

Identification of genes that control root system architecture in crop plants requires innovations that enable high-throughput and accurate measurements of root system architecture through time. We demonstrate the ability of a semiautomated 3D in vivo imaging and digital phenotyping pipeline to interrogate the quantitative genetic basis of root system growth in a rice biparental mapping population, Bala × Azucena. We phenotyped >1,400 3D root models and >57,000 2D images for a suite of 25 traits that quantified the distribution, shape, extent of exploration, and the intrinsic size of root networks at days 12, 14, and 16 of growth in a gellan gum medium. From these data we identified 89 quantitative trait loci, some of which correspond to those found previously in soil-grown plants, and provide evidence for genetic tradeoffs in root growth allocations, such as between the extent and thoroughness of exploration. We also developed a multivariate method for generating and mapping central root architecture phenotypes and used it to identify five major quantitative trait loci (r(2) = 24-37%), two of which were not identified by our univariate analysis. Our imaging and analytical platform provides a means to identify genes with high potential for improving root traits and agronomic qualities of crops.

GiA Roots: software for the high throughput analysis of plant root system architecture.

BACKGROUND: Characterizing root system architecture (RSA) is essential to understanding the development and function of vascular plants. Identifying RSA-associated genes also represents an underexplored opportunity for crop improvement. Software tools are needed to accelerate the pace at which quantitative traits of RSA are estimated from images of root networks. RESULTS: We have developed GiA Roots (General Image Analysis of Roots), a semi-automated software tool designed specifically for the high-throughput analysis of root system images. GiA Roots includes user-assisted algorithms to distinguish root from background and a fully automated pipeline that extracts dozens of root system phenotypes. Quantitative information on each phenotype, along with intermediate steps for full reproducibility, is returned to the end-user for downstream analysis. GiA Roots has a GUI front end and a command-line interface for interweaving the software into large-scale workflows. GiA Roots can also be extended to estimate novel phenotypes specified by the end-user. CONCLUSIONS: We demonstrate the use of GiA Roots on a set of 2393 images of rice roots representing 12 genotypes from the species Oryza sativa. We validate trait measurements against prior analyses of this image set that demonstrated that RSA traits are likely heritable and associated with genotypic differences. Moreover, we demonstrate that GiA Roots is extensible and an end-user can add functionality so that GiA Roots can estimate novel RSA traits. In summary, we show that the software can function as an efficient tool as part of a workflow to move from large numbers of root images to downstream analysis.

Imaging and analysis platform for automatic phenotyping and trait ranking of plant root systems.

The ability to nondestructively image and automatically phenotype complex root systems, like those of rice (Oryza sativa), is fundamental to identifying genes underlying root system architecture (RSA). Although root systems are central to plant fitness, identifying genes responsible for RSA remains an underexplored opportunity for crop improvement. Here we describe a nondestructive imaging and analysis system for automated phenotyping and trait ranking of RSA. Using this system, we image rice roots from 12 genotypes. We automatically estimate RSA traits previously identified as important to plant function. In addition, we expand the suite of features examined for RSA to include traits that more comprehensively describe monocot RSA but that are difficult to measure with traditional methods. Using 16 automatically acquired phenotypic traits for 2,297 images from 118 individuals, we observe (1) wide variation in phenotypes among the genotypes surveyed; and (2) greater intergenotype variance of RSA features than variance within a genotype. RSA trait values are integrated into a computational pipeline that utilizes supervised learning methods to determine which traits best separate two genotypes, and then ranks the traits according to their contribution to each pairwise comparison. This trait-ranking step identifies candidate traits for subsequent quantitative trait loci analysis and demonstrates that depth and average radius are key contributors to differences in rice RSA within our set of genotypes. Our results suggest a strong genetic component underlying rice RSA. This work enables the automatic phenotyping of RSA of individuals within mapping populations, providing an integrative framework for quantitative trait loci analysis of RSA.