Mitosis is a complex self-organising process that achieves high fidelity separation of duplicated chromosomes into two daughter cells through capture and alignment of chromosomes to the spindle mid-plane. Chromosome movements are driven by kinetochores (KTs), multi-protein machines that attach chromosomes to microtubules (MTs), and through those attachments both control and generate directional forces. Using lattice light sheet microscopy imaging and automated near-complete tracking of kinetochores at fine spatio-temporal resolution, we produce a detailed atlas of kinetochore metaphase-anaphase dynamics in untransformed human cells (RPE1). Such data allows dynamic models to be reverse engineered and biological hypotheses to be addressed. We determined the support from this dataset for 17 models of metaphase dynamics using Bayesian inference, demonstrating (1) substantial sister asymmetry that generates transverse organisation of the metaphase plate (MPP), (2) substantial spatial organisation of KT dynamic properties within the MPP, and (3) significant time dependence of the K-fiber mechanical parameters whereby K-fiber forces tune over the last 5 mins of metaphase towards a set point, referred to as the anaphase ready state. These spatio-temporal trends are robust to perturbation of the spindle assembly pathway (nocodazole washout treatment), suggesting that the underlying processes generating kinetochore heterogeneity are intrinsic to mitosis and possibly play a role in ensuring high-fidelity segregation.

In this paper, we present an innovative application of the Constant False Alarm Rate (CFAR) algorithm, traditionally utilized in radar signal processing, to enhance the accuracy and reliability of kinetochore (KT) tracking in live-cell lattice light-sheet microscopy (LLSM) imaging of human cells during the mitotic phase of cell division. Fluorescently labelled KTs appear as spots in diffraction-limited light microscopy. Traditional KT detection methods, such as adaptive thresholding, often struggle with the dynamic and noisy backgrounds of cells, leading to less efficient KT identification. By adapting the CFAR algorithm to the specific challenges of $K T$ detection in 3D, we present a method that offers improved precision and stability in detecting KTs across different stages of mitosis. The performance of the CFAR-KT method was rigorously compared to the adaptive thresholding approach across a cohort of 31 cells, with results highlighting CFAR-KT’s enhanced detection efficiency. Despite a slightly lower mean detection count compared to the adaptive method, the CFAR-KT method achieved lower false positives and a higher success rate in tracking KTs over the duration of the cell division process, underscoring its effectiveness in capturing the dynamics of KTs.

MOTIVATION: Lattice light-sheet microscopy (LLSM) is revolutionizing cell biology since it enables fast, high-resolution extended imaging in three dimensions combined with a drastic reduction in photo-toxicity and bleaching. However, analysis of such datasets still remains a major challenge. RESULTS: Automated tracking of kinetochores, the protein complex facilitating and controlling microtubule attachment of the chromosomes within the mitotic spindle, provides quantitative assessment of chromosome dynamics in mitosis. Here, we extend existing open-source kinetochore tracking software (KiT) to track (and pair) kinetochores throughout prometaphase to anaphase in LLSM data. One of the key improvements is a regularization term in the objective function to enforce biological information about the number of kinetochores in a human mitotic cell, as well as improved diagnostic tools. This software provides quantitative insights into how kinetochores robustly ensure congression and segregation of chromosomes during mitosis. AVAILABILITY AND IMPLEMENTATION: KiT is free, open-source software implemented in MATLAB and can be downloaded as a package from https://github.com/cmcb-warwick/KiT. The source repository is available at https://bitbucket.org/jarmond/kit (tag v2.4.0) and under continuing development. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Gaussian spot fitting methods have significantly extended the spatial range where fluorescent microscopy can be used, with recent techniques approaching nanometre (nm) resolutions. However, small inter-fluorophore distances are systematically over-estimated for typical molecular scales. This bias can be corrected computationally, but current algorithms are limited to correcting distances between pairs of fluorophores. Here we present a flexible Bayesian computational approach that infers the distances and angles between multiple fluorophores and has several advantages over these previous methods. Specifically it improves confidence intervals for small lengths, estimates measurement errors of each fluorophore individually and infers the correlations between polygon lengths. The latter is essential for determining the full multi-fluorophore 3D architecture. We further developed the algorithm to infer the mixture composition of a heterogeneous population of multiple polygon states. We use our algorithm to analyse the 3D architecture of the human kinetochore, a macro-molecular complex that is essential for high fidelity chromosome segregation during cell division. Using triple fluorophore image data we unravel the mixture of kinetochore states during human mitosis, inferring the conformation of microtubule attached and unattached kinetochores and their proportions across mitosis. We demonstrate that the attachment conformation correlates with intersister tension and sister alignment to the metaphase plate.

Achieving global food security for the estimated 9 billion people by 2050 is a major scientific challenge. Crop productivity is fundamentally restricted by the rate of fixation of atmospheric carbon. The dedicated enzyme, RubisCO, has a low turnover and poor specificity for CO2. This limitation of C3 photosynthesis (the basic carbon-assimilation pathway present in all plants) is alleviated in some lineages by use of carbon-concentrating-mechanisms, such as the C4 cycle-a biochemical pump that concentrates CO2 near RubisCO increasing assimilation efficacy. Most crops use only C3 photosynthesis, so one promising research strategy to boost their productivity focuses on introducing a C4 cycle. The simplest proposal is to use the cycle to concentrate CO2 inside individual chloroplasts. The photosynthetic efficiency would then depend on the leakage of CO2 out of a chloroplast. We examine this proposal with a 3D spatial model of carbon and oxygen diffusion and C4 photosynthetic biochemistry inside a typical C3-plant mesophyll cell geometry. We find that the cost-efficiency of C4 photosynthesis depends on the gas permeability of the chloroplast envelope, the C4 pathway having higher quantum efficiency than C3 for permeabilities below 300 μm/s. However, at higher permeabilities the C4 pathway still provides a substantial boost to carbon assimilation with only a moderate decrease in efficiency. The gains would be capped by the ability of chloroplasts to harvest light, but even under realistic light regimes a 100% boost to carbon assimilation is possible. This could be achieved in conjunction with lower investment in chloroplasts if their cell surface coverage is also reduced. Incorporation of this C4 cycle into C3 crops could thus promote higher growth rates and better drought resistance in dry, high-sunlight climates.

UNLABELLED: During mitosis, chromosomes are attached to the mitotic spindle via large protein complexes called kinetochores. The motion of kinetochores throughout mitosis is intricate and automated quantitative tracking of their motion has already revealed many surprising facets of their behaviour. Here, we present 'KiT' (Kinetochore Tracking)-an easy-to-use, open-source software package for tracking kinetochores from live-cell fluorescent movies. KiT supports 2D, 3D and multi-colour movies, quantification of fluorescence, integrated deconvolution, parallel execution and multiple algorithms for particle localization. AVAILABILITY AND IMPLEMENTATION: KiT is free, open-source software implemented in MATLAB and runs on all MATLAB supported platforms. KiT can be downloaded as a package from http://www.mechanochemistry.org/mcainsh/software.php The source repository is available at https://bitbucket.org/jarmond/kit and under continuing development. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online. CONTACT: [email protected].

Table of contents A1 Highlights from the eleventh ISCB Student Council Symposium 2015 Katie Wilkins, Mehedi Hassan, Margherita Francescatto, Jakob Jespersen, R. Gonzalo Parra, Bart Cuypers, Dan DeBlasio, Alexander Junge, Anupama Jigisha, Farzana Rahman O1 Prioritizing a drug’s targets using both gene expression and structural similarity Griet Laenen, Sander Willems, Lieven Thorrez, Yves Moreau O2 Organism specific protein-RNA recognition: A computational analysis of protein-RNA complex structures from different organisms Nagarajan Raju, Sonia Pankaj Chothani, C. Ramakrishnan, Masakazu Sekijima; M. Michael Gromiha O3 Detection of Heterogeneity in Single Particle Tracking Trajectories Paddy J Slator, Nigel J Burroughs O4 3D-NOME: 3D NucleOme Multiscale Engine for data-driven modeling of three-dimensional genome architecture Przemysław Szałaj, Zhonghui Tang, Paul Michalski, Oskar Luo, Xingwang Li, Yijun Ruan, Dariusz Plewczynski O5 A novel feature selection method to extract multiple adjacent solutions for viral genomic sequences classification Giulia Fiscon, Emanuel Weitschek, Massimo Ciccozzi, Paola Bertolazzi, Giovanni Felici O6 A Systems Biology Compendium for Leishmania donovani Bart Cuypers, Pieter Meysman, Manu Vanaerschot, Maya Berg, Hideo Imamura, Jean-Claude Dujardin, Kris Laukens O7 Unravelling signal coordination from large scale phosphorylation kinetic data Westa Domanova, James R. Krycer, Rima Chaudhuri, Pengyi Yang, Fatemeh Vafaee, Daniel J. Fazakerley, Sean J. Humphrey, David E. James, Zdenka Kuncic

Kinetochores are multi-protein complexes that mediate the physical coupling of sister chromatids to spindle microtubule bundles (called kinetochore (K)-fibres) from respective poles. These kinetochore-attached K-fibres generate pushing and pulling forces, which combine with polar ejection forces (PEF) and elastic inter-sister chromatin to govern chromosome movements. Classic experiments in meiotic cells using calibrated micro-needles measured an approximate stall force for a chromosome, but methods that allow the systematic determination of forces acting on a kinetochore in living cells are lacking. Here we report the development of mathematical models that can be fitted (reverse engineered) to high-resolution kinetochore tracking data, thereby estimating the model parameters and allowing us to indirectly compute the (relative) force components (K-fibre, spring force and PEF) acting on individual sister kinetochores in vivo. We applied our methodology to thousands of human kinetochore pair trajectories and report distinct signatures in temporal force profiles during directional switches. We found the K-fibre force to be the dominant force throughout oscillations, and the centromeric spring the smallest although it has the strongest directional switching signature. There is also structure throughout the metaphase plate, with a steeper PEF potential well towards the periphery and a concomitant reduction in plate thickness and oscillation amplitude. This data driven reverse engineering approach is sufficiently flexible to allow fitting of more complex mechanistic models; mathematical models of kinetochore dynamics can therefore be thoroughly tested on experimental data for the first time. Future work will now be able to map out how individual proteins contribute to kinetochore-based force generation and sensing.

MOTIVATION: There are a number of algorithms to infer causal regulatory networks from time series (gene expression) data. Here we analyse the phenomena of regulator interference, where regulators with similar dynamics mutually suppress both the probability of regulating a target and the associated link strength; for instance, interference between two identical strong regulators reduces link probabilities by ∼50%. RESULTS: We construct a robust method to define an interference-corrected causal network based on an analysis of the conditional link probabilities that recovers links lost through interference. On a large real network (Streptomyces coelicolor, phosphate depletion), we demonstrate that significant interference can occur between regulators with a correlation as low as 0.865, losing an estimated 34% of links by interference. However, levels of interference cannot be predicted from the correlation between regulators alone and are data specific. Validating against known networks, we show that high numbers of functional links are lost by regulator interference. Performance against other methods on DREAM4 data is excellent. AVAILABILITY AND IMPLEMENTATION: The method is implemented in R and is publicly available as the NIACS package at http://www2.warwick.ac.uk/fac/sci/systemsbiology/research/software.

Immune synapses formed by T and NK cells both show segregation of the integrin ICAM1 from other proteins such as CD2 (T cell) or KIR (NK cell). However, the mechanism by which these proteins segregate remains unclear; one key hypothesis is a redistribution based on protein size. Simulations of this mechanism qualitatively reproduce observed segregation patterns, but only in certain parameter regimes. Verifying that these parameter constraints in fact hold has not been possible to date, this requiring a quantitative coupling of theory to experimental data. Here, we address this challenge, developing a new methodology for analysing and quantifying image data and its integration with biophysical models. Specifically we fit a binding kinetics model to 2 colour fluorescence data for cytoskeleton independent synapses (2 and 3D) and test whether the observed inverse correlation between fluorophores conforms to size dependent exclusion, and further, whether patterned states are predicted when model parameters are estimated on individual synapses. All synapses analysed satisfy these conditions demonstrating that the mechanisms of protein redistribution have identifiable signatures in their spatial patterns. We conclude that energy processes implicit in protein size based segregation can drive the patternation observed in individual synapses, at least for the specific examples tested, such that no additional processes need to be invoked. This implies that biophysical processes within the membrane interface have a crucial impact on cell:cell communication and cell signalling, governing protein interactions and protein aggregation.

MOTIVATION: Gene expression measurements are the most common data source for reverse engineering gene interaction networks. When dealing with destructive sampling in time course experiments, it is common to average any available measurements for each time point and to treat this as the actual time series data for fitting the network, neglecting the variability contained in the repeated measurements. Proceeding in such a way can affect the retrieved network topology. RESULTS: We propose a fully Bayesian method for reverse engineering a gene interaction network, based on time course data with repeated measurements. The observations are treated as surrogate measurements of the underlying gene expression. As these measurements often contain outliers, we use a non-Gaussian specification for dealing with measurement error. The network interactions are assumed linear and an autoregressive model is specified, augmented with indicator variables that allow inference on the topology of the network. We analyse two in silico and one in vivo experiments, the latter dealing with the circadian clock in Arabidopsis thaliana. A systematic attenuation of the estimated regulation strengths and a concomitant overestimation of their precision is demonstrated when measurement error is disregarded. Thus, a clear improvement in the inferred topology for the synthetic datasets is demonstrated when this is included. Also, the influence of outliers in the retrieved network is demonstrated when using the in vivo data. AVAILABILITY: Matlab code and data used in the article are available from http://go.warwick.ac.uk/majuarez/home/materials.