Four-dimensional measurement of root system development using time-series three-dimensional volumetric data analysis by backward prediction

BACKGROUND: Root system architecture (RSA) is an essential characteristic for efficient water and nutrient absorption in terrestrial plants; its plasticity enables plants to respond to different soil environments. Better understanding of root plasticity is important in developing stress-tolerant crops. Non-invasive techniques that can measure roots in soils nondestructively, such as X-ray computed tomography (CT), are useful to evaluate RSA plasticity. However, although RSA plasticity can be measured by tracking individual root growth, only a few methods are available for tracking individual roots from time-series three-dimensional (3D) images. RESULTS: We developed a semi-automatic workflow that tracks individual root growth by vectorizing RSA from time-series 3D images via two major steps. The first step involves 3D alignment of the time-series RSA images by iterative closest point registration with point clouds generated by high-intensity particles in potted soils. This alignment ensures that the time-series RSA images overlap. The second step consists of backward prediction of vectorization, which is based on the phenomenon that the root length of the RSA vector at the earlier time point is shorter than that at the last time point. In other words, when CT scanning is performed at time point A and again at time point B for the same pot, the CT data and RSA vectors at time points A and B will almost overlap, but not where the roots have grown. We assumed that given a manually created RSA vector at the last time point of the time series, all RSA vectors except those at the last time point could be automatically predicted by referring to the corresponding RSA images. Using 21 time-series CT volumes of a potted plant of upland rice (Oryza sativa), this workflow revealed that the root elongation speed increased with age. Compared with a workflow that does not use backward prediction, the workflow with backward prediction reduced the manual labor time by 95%. CONCLUSIONS: We developed a workflow to efficiently generate time-series RSA vectors from time-series X-ray CT volumes. We named this workflow 'RSAtrace4D' and are confident that it can be applied to the time-series analysis of RSA development and plasticity. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s13007-022-00968-x.