Hypoxia in tomato (Solanum lycopersicum) fruit during ripening: Biophysical elucidation by a 3D reaction-diffusion model.

Respiration provides energy, substrates, and precursors to support physiological changes of the fruit during climacteric ripening. A key substrate of respiration is oxygen that needs to be supplied to the fruit in a passive way by gas transfer from the environment. Oxygen gradients may develop within the fruit due to its bulky size and the dense fruit tissues, potentially creating hypoxia that may have a role in the spatial development of ripening. This study presents a 3D reaction-diffusion model using tomato (Solanum lycopersicum) fruit as a test subject, combining the multiscale fruit geometry generated from magnetic resonance imaging and microcomputed tomography with varying respiration kinetics and contrasting boundary resistances obtained through independent experiments. The model predicted low oxygen levels in locular tissue under atmospheric conditions, and the oxygen level was markedly lower upon scar occlusion, aligning with microsensor profiling results. The locular region was in a hypoxic state, leading to its low aerobic respiration with high CO2 accumulation by fermentative respiration, while the rest of the tissues remained well oxygenated. The model further revealed that the hypoxia is caused by a combination of diffusion resistances and respiration rates of the tissue. Collectively, this study reveals the existence of the respiratory gas gradients and its biophysical causes during tomato fruit ripening, providing richer information for future studies on localized endogenous ethylene biosynthesis and fruit ripening.

Causes and consequences of endogenous hypoxia on growth and metabolism of developing maize kernels.

Maize (Zea mays) kernels are the largest cereal grains, and their endosperm is severely oxygen deficient during grain fill. The causes, dynamics, and mechanisms of acclimation to hypoxia are minimally understood. Here, we demonstrate that hypoxia develops in the small, growing endosperm, but not the nucellus, and becomes the standard state, regardless of diverse structural and genetic perturbations in modern maize (B73, popcorn, sweet corn), mutants (sweet4c, glossy6, waxy), and non-domesticated wild relatives (teosintes and Tripsacum species). We also uncovered an interconnected void space at the chalazal pericarp, providing superior oxygen supply to the placental tissues and basal endosperm transfer layer. Modeling indicated a very high diffusion resistance inside the endosperm, which, together with internal oxygen consumption, could generate steep oxygen gradients at the endosperm surface. Manipulation of oxygen supply induced reciprocal shifts in gene expression implicated in controlling mitochondrial functions (23.6 kDa Heat-Shock Protein, Voltage-Dependent Anion Channel 2) and multiple signaling pathways (core hypoxia genes, cyclic nucleotide metabolism, ethylene synthesis). Metabolite profiling revealed oxygen-dependent shifts in mitochondrial pathways, ascorbate metabolism, starch synthesis, and auxin degradation. Long-term elevated oxygen supply enhanced the rate of kernel development. Altogether, evidence here supports a mechanistic framework for the establishment of and acclimation to hypoxia in the maize endosperm.

Spatio-temporal dynamics of the metabolome of climacteric fruit during ripening and post-harvest storage.

Fruit quality traits are determined to a large extent by their metabolome. The metabolite content of climacteric fruit changes drastically during ripening and post-harvest storage, and has been investigated extensively. However, the spatial distribution of metabolites and how it changes in time has received much less attention as fruit are usually considered as homogenous plant organs. Yet, spatio-temporal changes of starch, which is hydrolyzed during ripening, has been used for a long time as a ripening index. As vascular transport of water, and hence convective transport of metabolites, slows down in mature fruit and even stalls after detachment, spatio-temporal changes in their concentration are probably affected by diffusive transport of gaseous molecules that act as substrate (O2), inhibitor (CO2), or regulator (ethylene and NO) of the metabolic pathways that are active during climacteric ripening. In this review, we discuss such spatio-temporal changes of the metabolome and how they are affected by transport of metabolic gases and gaseous hormones. As there are currently no techniques available to measure the metabolite distribution repeatedly by non-destructive means, we introduce reaction-diffusion models as an in silico tool to compute it. We show how the different components of such a model can be integrated and used to better understand the role of spatio-temporal changes of the metabolome in ripening and post-harvest storage of climacteric fruit that is detached from the plant, and discuss future research needs.

The role of chloroplast movement in C4 photosynthesis: a theoretical analysis using a three-dimensional reaction-diffusion model for maize.

Chloroplasts movement within mesophyll cells in C4 plants is hypothesized to enhance the CO2 concentrating mechanism, but this is difficult to verify experimentally. A three-dimensional (3D) leaf model can help analyse how chloroplast movement influences the operation of the CO2 concentrating mechanism. The first volumetric reaction-diffusion model of C4 photosynthesis that incorporates detailed 3D leaf anatomy, light propagation, ATP and NADPH production, and CO2, O2 and bicarbonate concentration driven by diffusional and assimilation/emission processes was developed. It was implemented for maize leaves to simulate various chloroplast movement scenarios within mesophyll cells: the movement of all mesophyll chloroplasts towards bundle sheath cells (aggregative movement) and movement of only those of interveinal mesophyll cells towards bundle sheath cells (avoidance movement). Light absorbed by bundle sheath chloroplasts relative to mesophyll chloroplasts increased in both cases. Avoidance movement decreased light absorption by mesophyll chloroplasts considerably. Consequently, total ATP and NADPH production and net photosynthetic rate increased for aggregative movement and decreased for avoidance movement compared with the default case of no chloroplast movement at high light intensities. Leakiness increased in both chloroplast movement scenarios due to the imbalance in energy production and demand in mesophyll and bundle sheath cells. These results suggest the need to design strategies for coordinated increases in electron transport and Rubisco activities for an efficient CO2 concentrating mechanism at very high light intensities.

Microstructural changes enhance oxygen transport in tomato (Solanum lycopersicum) fruit during maturation and ripening.

Climacteric ripening of tomato fruit is initiated by a characteristic surge of the production rate of ethylene, accompanied by an increase in respiration rate. As both activities consume O2 and produce CO2 , gas concentration gradients develop in the fruit that cause diffusive transport. This may, in turn, affect respiration and ethylene biosynthesis. Gas diffusion in fruit depends on the amount and connectivity of cells and intercellular spaces in 3D. We investigated micromorphological changes in different tomato tissues during development and ripening by visualizing cells and pores based on high-resolution micro-computed tomography, and computed effective O2 diffusivity coefficients based on microstructural features of the tissues. We demonstrated that mesocarp and septa tissues have larger cells but small and more disconnected pores than the placenta and columella, resulting in relatively lower effective O2 diffusivity coefficients. Cell disintegration occurred in the mesocarp and septa during ripening, indicating lysigenous air pore formation and resulting in a gradual increase of the effective O2 diffusivity. The results suggest that hypoxic conditions caused by the increasing size and, hence, diffusion resistance of the growing fruit may induce an increase of tissue porosity that results in a greatly enhanced O2 diffusivity and, thus, helps to alleviate them.

The initiation of bud burst in grapevine features dynamic regulation of the apoplastic pore size.

The physiological constraints on bud burst in woody perennials, including vascular development and oxygenation, remain unresolved. Both light and tissue oxygen status have emerged as important cues for vascular development in other systems; however, grapevine buds have only a facultative light requirement, and data on the tissue oxygen status have been confounded by the spatial variability within the bud. Here, we analysed apoplastic development at early stages of grapevine bud burst and combined molecular modelling with histochemical techniques to determine the pore size of cell walls in grapevine buds. The data demonstrate that quiescent grapevine buds were impermeable to apoplastic dyes (acid fuchsin and eosin Y) until after bud burst was established. The molecular exclusion size was calculated to be 2.1 nm, which would exclude most macromolecules except simple sugars and phytohormones until after bud burst. We used micro-computed tomography to demonstrate that tissue oxygen partial pressure data correlated well with structural heterogeneity of the bud and differences in tissue density, confirming that the primary bud complex becomes rapidly and preferentially oxygenated during bud burst. Taken together, our results reveal that the apoplastic porosity is highly regulated during the early stages of bud burst, suggesting a role for vascular development in the initial, rapid oxygenation of the primary bud complex.

In silico study of the role of cell growth factors in photosynthesis using a virtual leaf tissue generator coupled to a microscale photosynthesis gas exchange model.

Computational tools that allow in silico analysis of the role of cell growth and division on photosynthesis are scarce. We present a freely available tool that combines a virtual leaf tissue generator and a two-dimensional microscale model of gas transport during C3 photosynthesis. A total of 270 mesophyll geometries were generated with varying degrees of growth anisotropy, growth extent, and extent of schizogenous airspace formation in the palisade mesophyll. The anatomical properties of the virtual leaf tissue and microscopic cross-sections of actual leaf tissue of tomato (Solanum lycopersicum L.) were statistically compared. Model equations for transport of CO2 in the liquid phase of the leaf tissue were discretized over the geometries. The virtual leaf tissue generator produced a leaf anatomy of tomato that was statistically similar to real tomato leaf tissue. The response of photosynthesis to intercellular CO2 predicted by a model that used the virtual leaf tissue geometry compared well with measured values. The results indicate that the light-saturated rate of photosynthesis was influenced by interactive effects of extent and directionality of cell growth and degree of airspace formation through the exposed surface of mesophyll per leaf area. The tool could be used further in investigations of improving photosynthesis and gas exchange in relation to cell growth and leaf anatomy.

3-D microstructural changes in relation to the evolution of quality during ripening of mango (Mangifera indica L. cv. Carabao).

BACKGROUND: The ripening of mango involves changes in texture, flavor, and color, affecting the quality of the fruit. Previous studies have investigated the physiology on the evolution of quality during ripening but only a few have looked at microstructural changes during ripening. None of them has provided an insight into the relationhip between 3-D microstructure and the evolution of quality during ripening. As the 3-D microstructure of fruit tissue determines its mechanical and gas-transport properties, it is likely to affect fruit texture, respiratory metabolism, and other ripening processes. RESULTS: The present study focuses on the role of 3-D microstructural changes in relation to quality changes during mango ripening. Microstructural imaging using X-ray micro-computed tomography suggested the incidence of cell leakage, which was confirmed by the measurement of electrolyte leakage from the fruit peel. Due to cell leakage, porosity, pore connectivity, and pore local diameter were decreased whereas the tissue local diameter and pore specific area were increased. The decline in respiration and respiratory quotient during ripening followed the microstructural changes observed. Meanwhile, changes in aroma were observed such as a decrease in monoterpenes and an increase in esters and other fermentative metabolites. CONCLUSION: Overall, the results provide a complete, integrated picture of microstructural changes during ripening accompanying the evolution of fruit quality, suggesting functional relationships between the two. © 2020 Society of Chemical Industry.

Down-regulation of respiration in pear fruit depends on temperature.

The respiration rate of plant tissues decreases when the amount of available O2 is reduced. There is, however, a debate on whether the respiration rate is controlled either by diffusion limitation of oxygen or through regulatory processes at the level of the transcriptome. We used experimental and modelling approaches to demonstrate that both diffusion limitation and metabolic regulation affect the response of respiration of bulky plant organs such as fruit to reduced O2 levels in the surrounding atmosphere. Diffusion limitation greatly affects fruit respiration at high temperature, but at low temperature respiration is reduced through a regulatory process, presumably a response to a signal generated by a plant oxygen sensor. The response of respiration to O2 is time dependent and is highly sensitive, particularly at low O2 levels in the surrounding atmosphere. Down-regulation of the respiration at low temperatures may save internal O2 and relieve hypoxic conditions in the fruit.

Modelling postmortem evolution of pH in beef M. biceps femoris under two different cooling regimes.

A model based on enzyme kinetics was developed to predict differences in postmortem pH change in beef muscles as affected by cooling rate. For the calibration and validation of the model, pH and temperature measurements were conducted at different positions in M. biceps femoris following conventional carcass cooling or faster cooling of the muscle after hot boning. The glycogen conversion, and, hence, the pH fall, was observed to significantly vary with position and cooling regime but only during the initial hours of cooling. Comparison of the cooling regimes indicated that fast cooling following hot boning avoids heat shortening induced by the combined effect of high temperature and low pH.

Synchrotron X-ray computed laminography of the three-dimensional anatomy of tomato leaves.

Synchrotron radiation computed laminography (SR-CL) is presented as an imaging method for analyzing the three-dimensional (3D) anatomy of leaves. The SR-CL method was used to provide 3D images of 1-mm² samples of intact leaves at a pixel resolution of 750 nm. The method allowed visualization and quantitative analysis of palisade and spongy mesophyll cells, and showed local venation patterns, aspects of xylem vascular structure and stomata. The method failed to image subcellular organelles such as chloroplasts. We constructed 3D computer models of leaves that can provide a basis for calculating gas exchange, light penetration and water and solute transport. The leaf anatomy of two different tomato genotypes grown in saturating light conditions was compared by 3D analysis. Differences were found in calculated values of tissue porosity, cell number density, cell area to volume ratio and cell volume and cell shape distributions of palisade and spongy cell layers. In contrast, the exposed cell area to leaf area ratio in mesophyll, a descriptor that correlates to the maximum rate of photosynthesis in saturated light conditions, was no different between spongy and palisade cells or between genotypes. The use of 3D image processing avoids many of the limitations of anatomical analysis with two-dimensional sections.

Three-dimensional microscale modelling of CO2 transport and light propagation in tomato leaves enlightens photosynthesis.

We present a combined three-dimensional (3-D) model of light propagation, CO2 diffusion and photosynthesis in tomato (Solanum lycopersicum L.) leaves. The model incorporates a geometrical representation of the actual leaf microstructure that we obtained with synchrotron radiation X-ray laminography, and was evaluated using measurements of gas exchange and leaf optical properties. The combination of the 3-D microstructure of leaf tissue and chloroplast movement induced by changes in light intensity affects the simulated CO2 transport within the leaf. The model predicts extensive reassimilation of CO2 produced by respiration and photorespiration. Simulations also suggest that carbonic anhydrase could enhance photosynthesis at low CO2 levels but had little impact on photosynthesis at high CO2 levels. The model confirms that scaling of photosynthetic capacity with absorbed light would improve efficiency of CO2 fixation in the leaf, especially at low light intensity.

Automatic analysis of the 3-D microstructure of fruit parenchyma tissue using X-ray micro-CT explains differences in aeration.

BACKGROUND: 3D high-resolution X-ray imaging methods have emerged over the last years for visualising the anatomy of tissue samples without substantial sample preparation. Quantitative analysis of cells and intercellular spaces in these images has, however, been difficult and was largely based on manual image processing. We present here an automated procedure for processing high-resolution X-ray images of parenchyma tissues of apple (Malus × domestica Borkh.) and pear (Pyrus communis L.) as a rapid objective method for characterizing 3D plant tissue anatomy at the level of single cells and intercellular spaces. RESULTS: We isolated neighboring cells in 3D images of apple and pear cortex tissues, and constructed a virtual sieve to discard incorrectly segmented cell particles or unseparated clumps of cells. Void networks were stripped down until their essential connectivity features remained. Statistical analysis of structural parameters showed significant differences between genotypes in the void and cell networks that relate to differences in aeration properties of the tissues. CONCLUSIONS: A new model for effective oxygen diffusivity of parenchyma tissue is proposed that not only accounts for the tortuosity of interconnected voids, but also for significant diffusion across cells where the void network is not connected. This will significantly aid interpretation and analysis of future tissue aeration studies. The automated image analysis methodology will also support pheno- and genotyping studies where the 3D tissue anatomy plays a role.

Modeling the propagation of light in realistic tissue structures with MMC-fpf: a meshed Monte Carlo method with free phase function.

Monte Carlo methods commonly used in tissue optics are limited to a layered tissue geometry and thus provide only a very rough approximation for many complex media such as biological structures. To overcome these limitations, a Meshed Monte Carlo method with flexible phase function choice (fpf-MC) has been developed to function in a mesh. This algorithm can model the light propagation in any complexly shaped structure, by attributing optical properties to the different mesh elements. Furthermore, this code allows the use of different discretized phase functions for each tissue type, which can be simulated from the microstructural properties of the tissue, in combination with a tool for simulating the bulk optical properties of polydisperse suspensions. As a result, the scattering properties of tissues can be estimated from information on the microstructural properties of the tissue. This is important for the estimation of the bulk optical properties that can be used for the light propagation model, since many types of tissue have never been characterized in literature. The combination of these contributions, made it possible to use the MMC-fpf for modeling the light porapagation in plant tissue. The developed Meshed Monte Carlo code with flexible phase function choice (MMC-fpf) was successfully validated in simulation through comparison with the Monte Carlo code in Multi-Layered tissues (R2 > 0.9999) and experimentally by comparing the measured and simulated reflectance (RMSE = 0.015%) and transmittance (RMSE = 0.0815%) values for tomato leaves.

Quantitative 3D shape description of dust particles from treated seeds by means of X-ray micro-CT.

Crop seeds are often treated with pesticides before planting. Pesticide-laden dust particles can be abraded from the seed coating during planting and expelled into the environment, damaging nontarget organisms. Drift of these dust particles depends on their size, shape and density. In this work, we used X-ray micro-CT to examine the size, shape (sphericity) and porosity of dust particles from treated seeds of various crops. The dust properties quantified in this work were very variable in different crops. This variability may be a result of seed morphology, seed batch, treatment composition, treatment technology, seed cleaning or an interaction of these factors. The intraparticle porosity of seed treatment dust particles varied from 0.02 to 0.51 according to the crop and generally increased with particle size. Calculated settling velocities demonstrated that accounting for particle shape and porosity is important in drift studies. For example, the settling velocity of dust particles with an equivalent diameter of 200 µm may vary between 0.1 and 1.2 m s(-1), depending on their shape and density. Our analysis shows that in a wind velocity of 5 m s(-1), such particles ejected at 1 m height may travel between 4 and 50 m from the source before settling. Although micro-CT is a valuable tool to characterize dust particles, the current image processing methodology limits the number of particles that can be analyzed.

Quantitative neutron imaging of water distribution, venation network and sap flow in leaves.

MAIN CONCLUSION: Quantitative neutron imaging is a promising technique to investigate leaf water flow and transpiration in real time and has perspectives towards studies of plant response to environmental conditions and plant water stress. The leaf hydraulic architecture is a key determinant of plant sap transport and plant-atmosphere exchange processes. Non-destructive imaging with neutrons shows large potential for unveiling the complex internal features of the venation network and the transport therein. However, it was only used for two-dimensional imaging without addressing flow dynamics and was still unsuccessful in accurate quantification of the amount of water. Quantitative neutron imaging was used to investigate, for the first time, the water distribution in veins and lamina, the three-dimensional venation architecture and sap flow dynamics in leaves. The latter was visualised using D2O as a contrast liquid. A high dynamic resolution was obtained by using cold neutrons and imaging relied on radiography (2D) as well as tomography (3D). The principle of the technique was shown for detached leaves, but can be applied to in vivo leaves as well. The venation network architecture and the water distribution in the veins and lamina unveiled clear differences between plant species. The leaf water content could be successfully quantified, though still included the contribution of the leaf dry matter. The flow measurements exposed the hierarchical structure of the water transport pathways, and an accurate quantification of the absolute amount of water uptake in the leaf was possible. Particular advantages of neutron imaging, as compared to X-ray imaging, were identified. Quantitative neutron imaging is a promising technique to investigate leaf water flow and transpiration in real time and has perspectives towards studies of plant response to environmental conditions and plant water stress.

The mechanism of improved aeration due to gas films on leaves of submerged rice.

Some terrestrial wetland plants, such as rice, have super-hydrophobic leaf surfaces which retain a gas film when submerged. O2 movement through the diffusive boundary layer (DBL) of floodwater, gas film and stomata into leaf mesophyll was explored by means of a reaction-diffusion model that was solved in a three-dimensional leaf anatomy model. The anatomy and dark respiration of leaves of rice (Oryza sativa L.) were measured and used to compute O2 fluxes and partial pressure of O2 (pO2 ) in the DBL, gas film and leaf when submerged. The effects of floodwater pO2 , DBL thickness, cuticle permeability, presence of gas film and stomatal opening were explored. Under O2 -limiting conditions of the bulk water (pO2 < 10 kPa), the gas film significantly increases the O2 flux into submerged leaves regardless of whether stomata are fully or partly open. With a gas film, tissue pO2 substantially increases, even for the slightest stomatal opening, but not when stomata are completely closed. The effect of gas films increases with decreasing cuticle permeability. O2 flux and tissue pO2 decrease with increasing DBL thickness. The present modelling analysis provides a mechanistic understanding of how leaf gas films facilitate O2 entry into submerged plants.

Cross-scale modelling of transpiration from stomata via the leaf boundary layer.

BACKGROUND AND AIMS: Leaf transpiration is a key parameter for understanding land surface-climate interactions, plant stress and plant structure­function relationships. Transpiration takes place at the microscale level, namely via stomata that are distributed discretely over the leaf surface with a very low surface coverage (approx. 0·2-5%). The present study aims to shed more light on the dependency of the leaf boundary-layer conductance (BLC) on stomatal surface coverage and air speed. METHODS: An innovative three-dimensional cross-scale modelling approach was applied to investigate convective mass transport from leaves, using computational fluid dynamics. The gap between stomatal and leaf scale was bridged by including all these scales in the same computational model (10⁻5-10⁻¹ m), which implies explicitly modelling individual stomata. KEY RESULTS: BLC was strongly dependent on stomatal surface coverage and air speed. Leaf BLC at low surface coverage ratios (CR), typical for stomata, was still relatively high, compared with BLC of a fully wet leaf (hypothetical CR of 100%). Nevertheless, these conventional BLCs (CR of 100%), as obtained from experiments or simulations on leaf models, were found to overpredict the convective exchange. In addition, small variations in stomatal CR were found to result in large variations in BLCs. Furthermore, stomata of a certain size exhibited a higher mass transfer rate at lower CRs. CONCLUSIONS: The proposed cross-scale modelling approach allows us to increase our understanding of transpiration at the sub-leaf level as well as the boundary-layer microclimate in a way currently not feasible experimentally. The influence of stomatal size, aperture and surface density, and also flow-field parameters can be studied using the model, and prospects for further improvement of the model are presented. An important conclusion of the study is that existing measures of conductances (e.g. from artificial leaves) can be significantly erroneous because they do not account for microscopic stomata, but instead assume a uniform distribution of evaporation such as found for a fully-wet leaf. The model output can be used to correct or upgrade existing BLCs or to feed into higher-scale models, for example within a multiscale framework.

A plant cell division algorithm based on cell biomechanics and ellipse-fitting.

BACKGROUND AND AIMS: The importance of cell division models in cellular pattern studies has been acknowledged since the 19th century. Most of the available models developed to date are limited to symmetric cell division with isotropic growth. Often, the actual growth of the cell wall is either not considered or is updated intermittently on a separate time scale to the mechanics. This study presents a generic algorithm that accounts for both symmetrically and asymmetrically dividing cells with isotropic and anisotropic growth. Actual growth of the cell wall is simulated simultaneously with the mechanics. METHODS: The cell is considered as a closed, thin-walled structure, maintained in tension by turgor pressure. The cell walls are represented as linear elastic elements that obey Hooke's law. Cell expansion is induced by turgor pressure acting on the yielding cell-wall material. A system of differential equations for the positions and velocities of the cell vertices as well as for the actual growth of the cell wall is established. Readiness to divide is determined based on cell size. An ellipse-fitting algorithm is used to determine the position and orientation of the dividing wall. The cell vertices, walls and cell connectivity are then updated and cell expansion resumes. Comparisons are made with experimental data from the literature. KEY RESULTS: The generic plant cell division algorithm has been implemented successfully. It can handle both symmetrically and asymmetrically dividing cells coupled with isotropic and anisotropic growth modes. Development of the algorithm highlighted the importance of ellipse-fitting to produce randomness (biological variability) even in symmetrically dividing cells. Unlike previous models, a differential equation is formulated for the resting length of the cell wall to simulate actual biological growth and is solved simultaneously with the position and velocity of the vertices. CONCLUSIONS: The algorithm presented can produce different tissues varying in topological and geometrical properties. This flexibility to produce different tissue types gives the model great potential for use in investigations of plant cell division and growth in silico.

Cyclist drag in team pursuit: influence of cyclist sequence, stature, and arm spacing.

In team pursuit, the drag of a group of cyclists riding in a pace line is dependent on several factors, such as anthropometric characteristics (stature) and position of each cyclist as well as the sequence in which they ride. To increase insight in drag reduction mechanisms, the aerodynamic drag of four cyclists riding in a pace line was investigated, using four different cyclists, and for four different sequences. In addition, each sequence was evaluated for two arm spacings. Instead of conventional field or wind tunnel experiments, a validated numerical approach (computational fluid dynamics) was used to evaluate cyclist drag, where the bicycles were not included in the model. The cyclist drag was clearly dependent on his position in the pace line, where second and subsequent positions experienced a drag reduction up to 40%, compared to an individual cyclist. Individual differences in stature and position on the bicycle led to an intercyclist variation of this drag reduction at a specific position in the sequence, but also to a variation of the total drag of the group for different sequences. A larger drag area for the group was found when riding with wider arm spacing. Such numerical studies on cyclists in a pace line are useful for determining the optimal cyclist sequence for team pursuit.