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 three-dimensional reaction-diffusion model using tomato (Solanum lycopersicum) fruit as a test subject, combining the multiscale fruit geometry generated from magnetic resonance imaging and micro-computed 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.

Automatic 3D cell segmentation of fruit parenchyma tissue from X-ray micro CT images using deep learning.

BACKGROUND: High quality 3D information of the microscopic plant tissue morphology-the spatial organization of cells and intercellular spaces in tissues-helps in understanding physiological processes in a wide variety of plants and tissues. X-ray micro-CT is a valuable tool that is becoming increasingly available in plant research to obtain 3D microstructural information of the intercellular pore space and individual pore sizes and shapes of tissues. However, individual cell morphology is difficult to retrieve from micro-CT as cells cannot be segmented properly due to negligible density differences at cell-to-cell interfaces. To address this, deep learning-based models were trained and tested to segment individual cells using X-ray micro-CT images of parenchyma tissue samples from apple and pear fruit with different cell and porosity characteristics. RESULTS: The best segmentation model achieved an Aggregated Jaccard Index (AJI) of 0.86 and 0.73 for apple and pear tissue, respectively, which is an improvement over the current benchmark method that achieved AJIs of 0.73 and 0.67. Furthermore, the neural network was able to detect other plant tissue structures such as vascular bundles and stone cell clusters (brachysclereids), of which the latter were shown to strongly influence the spatial organization of pear cells. Based on the AJIs, apple tissue was found to be easier to segment, as the porosity and specific surface area of the pore space are higher and lower, respectively, compared to pear tissue. Moreover, samples with lower pore network connectivity, proved very difficult to segment. CONCLUSIONS: The proposed method can be used to automatically quantify 3D cell morphology of plant tissue from micro-CT instead of opting for laborious manual annotations or less accurate segmentation approaches. In case fruit tissue porosity or pore network connectivity is too low or the specific surface area of the pore space too high, native X-ray micro-CT is unable to provide proper marker points of cell outlines, and one should rely on more elaborate contrast-enhancing scan protocols.

Unraveling the Kinetics of the 10-23 RNA-Cleaving DNAzyme.

DNA-based enzymes, or DNAzymes, are single-stranded DNA sequences with the ability to catalyze various chemical reactions, including the cleavage of the bond between two RNA nucleotides. Lately, an increasing interest has been observed in these RNA-cleaving DNAzymes in the biosensing and therapeutic fields for signal generation and the modulation of gene expression, respectively. Additionally, multiple efforts have been made to study the effects of the reaction environment and the sequence of the catalytic core on the conversion of the substrate into product. However, most of these studies have only reported alterations of the general reaction course, but only a few have focused on how each individual reaction step is affected. In this work, we present for the first time a mathematical model that describes and predicts the reaction of the 10-23 RNA-cleaving DNAzyme. Furthermore, the model has been employed to study the effect of temperature, magnesium cations and shorter substrate-binding arms of the DNAzyme on the different kinetic rate constants, broadening the range of conditions in which the model can be exploited. In conclusion, this work depicts the prospects of such mathematical models to study and anticipate the course of a reaction given a particular environment.

Increasing the Robustness of SIFT-MS Volatilome Fingerprinting by Introducing Notional Analyte Concentrations.

Selected ion flow tube-mass spectrometry (SIFT-MS) is an analytical technique for volatile detection and quantification. SIFT-MS can be applied in a "white box" approach, measuring concentrations of target compounds, or as a "black box" fingerprinting technique, scanning all product ions during a full scan. Combining SIFT-MS full scan data acquired from multibatches or large-scale experiments remains problematic due to signal fluctuation over time. The standard approach of normalizing full scan data to the total signal intensity was insufficient. This study proposes a new approach to correct SIFT-MS fingerprinting data. In this concept, all of the product ions from a full scan are considered individual compounds for which notional concentrations can be calculated. Converting ion count rates into notional analyte concentrations accounts for any changes in the instrument parameters. The benefits of the proposed approach are demonstrated on three years of data from both multibatches and long-term experiments showing a significant reduction of system-induced fluctuations providing a better focus on the changes of interest.

Targeted system approach to ethylene biosynthesis and signaling of a heat tolerant tomato cultivar; the impact of growing season on fruit ripening.

Growing tomato in hot weather conditions is challenging for fruit production and yield. Tomato cv. Savior is a heat-tolerant cultivar which can be grown during both the Vietnamese winter (mild condition) and summer (hot condition) season. Understanding the mechanisms of ethylene biosynthesis and signaling are important for agriculture, as manipulation of these pathways can lead to improvements in crop yield, stress tolerance, and fruit ripening. The objective of this study was to investigate an overview of ethylene biosynthesis and signaling from target genes to proteins and metabolites and the impact of growing season on a heat tolerant tomato cultivar throughout fruit ripening and postharvest storage. This work also showed the feasibility of absolute protein quantification of ethylene biosynthesis enzymes. Summer fruit showed the delayed peak of ethylene production until the red ripe stage. The difference in postharvest ethylene production between winter and summer fruit appears to be regulated by the difference in accumulation of 1-aminocyclopropane-1-carboxylic acid (ACC) which depends on the putative up-regulation of SAM levels. The lack of differences in protein concentrations between winter and summer fruit indicate that heat stress did not alter the ethylene biosynthesis-related protein abundance in heat tolerant cultivar. The analysis results of enzymatic activity and proteomics showed that in both winter and summer fruit, the majority of ACO activity could be mainly contributed to the abundance of ACO5 and ACO6 isoforms, rather than ACO1. Likewise, ethylene signal transduction was largely controlled by the abundance of ethylene receptors ETR1, ETR3, ETR6, and ETR7 together with the constitute triple response regulator CTR1 for both winter and summer grown tomatoes. Altogether our results indicate that in the heat tolerant tomato cv. Savior, growing season mainly affects the ethylene biosynthesis pathway and leaves the signaling pathway relatively unaffected.

Early detection of Botrytis cinerea in strawberry fruit during quiescent infection using selected ion flow tube mass spectrometry (SIFT-MS).

Botrytis cinerea is a devastating pathogen that can cause huge postharvest losses of strawberry. Although this fungus usually infects strawberries through their flowers, symptoms mainly appear when fruit are fully mature. A fast and sensitive method to detect and quantify the fungal infection, prior to symptom development, is, therefore, needed. In this study, we explore the possibility of using the strawberry volatilome to identify biomarkers for B. cinerea infection. Strawberry flowers were inoculated with B. cinerea to mimic the natural infection. First, quantitative polymerase chain reaction (qPCR) was used to quantify B. cinerea in the strawberry fruit. The detection limit of qPCR for B. cinerea DNA extracted from strawberries was 0.01 ng. Subsequently, changes in the fruit volatilome at different fruit developmental stages were characterized using gas chromatography - mass spectrometry (GC-MS) and selected ion flow tube mass spectrometry (SIFT-MS). Based on GC-MS data, 1-octen-3-ol produced by B. cinerea was confirmed as a potential biomarker of B. cinerea infection. Moreover, the product ion NO+ 127, obtained by SIFT-MS measurements, was proposed as a potential biomarker for B. cinerea infection by comparing its relative level with that of 1-octen-3-ol (obtained by GC-MS) and B. cinerea (obtained by qPCR). Separate PLS regressions were carried out for each developmental stages, and 11 product ions were significantly altered at all developmental stages. Finally, PLS regressions using these 11 ions as variables allowed the discrimination between samples containing different amount of B. cinerea. This work showed that profiling the fruit's volatilome using SIFT-MS can be used as a potential alternative to detect B. cinerea during the quiescent stage of B. cinerea infection prior to symptom development. Moreover, the corresponding compounds of potential biomarkers suggest that the volatile changes caused by B. cinerea infection may contribute to strawberry defense.

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.

Innovative Fabrication of Hollow Microneedle Arrays Enabling Blood Sampling with a Self-Powered Microfluidic Patch.

Microneedles are gaining a lot of attention in the context of sampling cutaneous biofluids such as capillary blood. Their minimal invasiveness and user-friendliness make them a prominent substitute for venous puncture or finger-pricking. Although the latter is suitable for self-sampling, the impracticality of manual handling and the difficulty of obtaining enough qualitative sample is driving the search for better solutions. In this context, hollow microneedle arrays (HMNAs) are particularly interesting for completely integrating sample-to-answer solutions as they create a duct between the skin and the sampling device. However, the fabrication of sharp-tipped HMNAs with a high aspect ratio (AR) is challenging, especially since a length of ≥1500 µm is desired to reach the blood capillaries. In this paper, we first described a novel two-step fabrication protocol for HMNAs in stainless steel by percussion laser drilling and subsequent micro-milling. The HMNAs were then integrated into a self-powered microfluidic sampling patch, containing a capillary pump which was optimized to generate negative pressure differences up to 40.9 ± 1.8 kPa. The sampling patch was validated in vitro, showing the feasibility of sampling 40 µL of liquid. It is anticipated that our proof-of-concept is a starting point for more sophisticated all-in-one biofluid sampling and point-of-care testing systems.

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.

Mechanical damages and packaging methods along the fresh fruit supply chain: A review.

Mechanical damage of fresh fruit occurs throughout the postharvest supply chain leading to poor consumer acceptance and marketability. In this review, the mechanisms of damage development are discussed first. Mathematical modeling provides advanced ways to describe and predict the deformation of fruit with arbitrary geometry, which is important to understand their mechanical responses to external forces. Also, the effects of damage at the cellular and molecular levels are discussed as this provides insight into fruit physiological responses to damage. Next, direct measurement methods for damage including manual evaluation, optical detection, magnetic resonance imaging, and X-ray computed tomography are examined, as well as indirect methods based on physiochemical indexes. Also, methods to measure fruit susceptibility to mechanical damage based on the bruise threshold and the amount of damage per unit of impact energy are reviewed. Further, commonly used external and interior packaging and their applications in reducing damage are summarized, and a recent biomimetic approach for designing novel lightweight packaging inspired by the fruit pericarp. Finally, future research directions are provided.HIGHLIGHTSMathematical modeling has been increasingly used to calculate damage to fruit.Cell and molecular mechanisms response to fruit damage is an under-explored area.Susceptibility measurement of different mechanical forces has received attention.Customized design of reusable and biodegradable packaging is a hot topic of research.

Sexual repurposing of juvenile aposematism in locusts.

Adaptive plasticity requires an integrated suite of functional responses to environmental variation, which can include social communication across life stages. Desert locusts (Schistocerca gregaria) exhibit an extreme example of phenotypic plasticity called phase polyphenism, in which a suite of behavioral and morphological traits differ according to local population density. Male and female juveniles developing at low population densities exhibit green- or sand-colored background-matching camouflage, while at high densities they show contrasting yellow and black aposematic patterning that deters predators. The predominant background colors of these phenotypes (green/sand/yellow) all depend on expression of the carotenoid-binding "Yellow Protein" (YP). Gregarious (high-density) adults of both sexes are initially pinkish, before a YP-mediated yellowing reoccurs upon sexual maturation. Yellow color is especially prominent in gregarious males, but the reason for this difference has been unknown since phase polyphenism was first described in 1921. Here, we use RNA interference to show that gregarious male yellowing acts as an intrasexual warning signal, which forms a multimodal signal with the antiaphrodisiac pheromone phenylacetonitrile (PAN) to prevent mistaken sexual harassment from other males during scramble mating in a swarm. Socially mediated reexpression of YP thus adaptively repurposes a juvenile signal that deters predators into an adult signal that deters undesirable mates. These findings reveal a previously underappreciated sexual dimension to locust phase polyphenism, and promote locusts as a model for investigating the relative contributions of natural versus sexual selection in the evolution of phenotypic plasticity.

Constitutive Defense Mechanisms Have a Major Role in the Resistance of Woodland Strawberry Leaves Against Botrytis cinerea.

The necrotrophic fungus Botrytis cinerea is a major threat to strawberry cultivation worldwide. By screening different Fragaria vesca genotypes for susceptibility to B. cinerea, we identified two genotypes with different resistance levels, a susceptible genotype F. vesca ssp. vesca Tenno 3 (T3) and a moderately resistant genotype F. vesca ssp. vesca Kreuzkogel 1 (K1). These two genotypes were used to identify the molecular basis for the increased resistance of K1 compared to T3. Fungal DNA quantification and microscopic observation of fungal growth in woodland strawberry leaves confirmed that the growth of B. cinerea was restricted during early stages of infection in K1 compared to T3. Gene expression analysis in both genotypes upon B. cinerea inoculation suggested that the restricted growth of B. cinerea was rather due to the constitutive resistance mechanisms of K1 instead of the induction of defense responses. Furthermore, we observed that the amount of total phenolics, total flavonoids, glucose, galactose, citric acid and ascorbic acid correlated positively with higher resistance, while H2O2 and sucrose correlated negatively. Therefore, we propose that K1 leaves are more resistant against B. cinerea compared to T3 leaves, prior to B. cinerea inoculation, due to a lower amount of innate H2O2, which is attributed to a higher level of antioxidants and antioxidant enzymes in K1. To conclude, this study provides important insights into the resistance mechanisms against B. cinerea, which highly depend on the innate antioxidative profile and specialized metabolites of woodland strawberry leaves.

Time Is of the Essence-Early Activation of the Mevalonate Pathway in Apple Challenged With Gray Mold Correlates With Reduced Susceptibility During Postharvest Storage.

Apple is typically stored under low temperature and controlled atmospheric conditions to ensure a year round supply of high quality fruit for the consumer. During storage, losses in quality and quantity occur due to spoilage by postharvest pathogens. One important postharvest pathogen of apple is Botrytis cinerea. The fungus is a broad host necrotroph with a large arsenal of infection strategies able to infect over 1,400 different plant species. We studied the apple-B. cinerea interaction to get a better understanding of the defense response in apple. We conducted an RNAseq experiment in which the transcriptome of inoculated and non-inoculated (control and mock) apples was analyzed at 0, 1, 12, and 28 h post inoculation. Our results show extensive reprogramming of the apple's transcriptome with about 28.9% of expressed genes exhibiting significant differential regulation in the inoculated samples. We demonstrate the transcriptional activation of pathogen-triggered immunity and a reprogramming of the fruit's metabolism. We demonstrate a clear transcriptional activation of secondary metabolism and a correlation between the early transcriptional activation of the mevalonate pathway and reduced susceptibility, expressed as a reduction in resulting lesion diameters. This pathway produces the building blocks for terpenoids, a large class of compounds with diverging functions including defense. 1-MCP and hot water dip treatment are used to further evidence the key role of terpenoids in the defense and demonstrate that ethylene modulates this response.

Quality Evolution and Aroma Profile of Pointed Cabbage in Different Storage Regimes.

With its increasing popularity, the need for optimal storage conditions of pointed cabbages becomes more important to meet the year-round demand. Storage of the pointed varieties, however, is more difficult compared to the traditional, round varieties and is limited to a few weeks in normal air. Pointed cabbages are more susceptible to quality loss (shriveling, yellowing of leaves, weight loss, fungal, and bacterial infections) and tend to spoil much faster. In order to provide a year-round availability of the fresh product, storage under controlled atmosphere (CA) could offer a solution. In this study, pointed, white cabbage heads (Brassica oleracea var. capitata for. alba L. subv. Conica cv. 'Caraflex') were stored at 1°C from November 2018 to May 2019 under four different CA conditions (1 kPa O2 + 1.5 kPa CO2, 1 kPa O2 + 5 kPa CO2, 3 kPa O2 + 1.5 kPa CO2, and 3 kPa O2 + 5 kPa CO2), and compared to storage under normal air. Results showed that CA storage resulted in a prolonged storage life with a good quality retention for both texture and aroma. CA-stored cabbages showed less weight loss, shriveling, and yellowing. Internal quality parameters [color, soluble solids content (SSC)] were stable over the whole storage period for all objects. The aroma profiles of both the storage atmosphere and cabbage samples were impacted by storage duration. The aroma of cabbage juice was also affected by the storage regime. A clear separation was found for cabbage stored under CA compared to the reference group. From the CA-treatments studied, a combination of low oxygen (1 kPa O2) and elevated carbon dioxide levels (5 kPa CO2) showed the best results maintaining quality. Storage under CA resulted in a better resemblance to the aroma of freshly, harvested produce compared to cabbages stored in normal air.

Robust dynamic experiments for the precise estimation of respiration and fermentation parameters of fruit and vegetables.

Dynamic models based on non-linear differential equations are increasingly being used in many biological applications. Highly informative dynamic experiments are valuable for the identification of these dynamic models. The storage of fresh fruit and vegetables is one such application where dynamic experimentation is gaining momentum. In this paper, we construct optimal O2 and CO2 gas input profiles to estimate the respiration and fermentation kinetics of pear fruit. The optimal input profiles, however, depend on the true values of the respiration and fermentation parameters. Locally optimal design of input profiles, which uses a single initial guess for the parameters, is the traditional method to deal with this issue. This method, however, is very sensitive to the initial values selected for the model parameters. Therefore, we present a robust experimental design approach that can handle uncertainty on the model parameters.

Computer-aided food engineering.

Computer-aided food engineering (CAFE) can reduce resource use in product, process and equipment development, improve time-to-market performance, and drive high-level innovation in food safety and quality. Yet, CAFE is challenged by the complexity and variability of food composition and structure, by the transformations food undergoes during processing and the limited availability of comprehensive mechanistic frameworks describing those transformations. Here we introduce frameworks to model food processes and predict physiochemical properties that will accelerate CAFE. We review how investments in open access, such as code sharing, and capacity-building through specialized courses could facilitate the use of CAFE in the transformation already underway in digital food systems.

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.

X-ray computed tomography for 3D plant imaging.

X-ray computed tomography (CT) is a valuable tool for 3D imaging of plant tissues and organs. Applications include the study of plant development and organ morphogenesis, as well as modeling of transport processes in plants. Some challenges remain, however, including attaining higher contrast for easier quantification, increasing the resolution for imaging subcellular features, and decreasing image acquisition and processing time for high-throughput phenotyping. In addition, phase contrast, multispectral, dark-field, soft X-ray, and time-resolved imaging are emerging. At the same time, a large amount of 3D image data are becoming available, posing challenges for data management. We review recent advances in the area of X-ray CT for plant imaging, and describe opportunities for using such images for studying transport processes in plants.

3D Printing of Monolithic Capillarity-Driven Microfluidic Devices for Diagnostics.

Rapid diagnostic testing at the site of the patient is essential when a fully equipped laboratory is not accessible. To maximize the impact of this approach, low-cost, disposable tests that require minimal user-interference and external equipment are desired. Fluid transport by capillary wicking removes the need for bulky ancillary equipment to actuate and control fluid flow. Nevertheless, current microfluidic paper-based analytical devices based on this principle struggle with the implementation of multistep diagnostic protocols because of fabrication-related issues. Here, 3D-printed microfluidic devices are demonstrated in a proof-of-concept enzyme-linked immunosorbent assay in which a multistep assay timeline is completed by precisely engineering capillary wetting within printed porous bodies. 3D printing provides a scalable route to low-cost microfluidic devices and obviates the assembly of discrete components. The resulting rapid and seamless transition between digital data and physical objects allows for rapid design iterations, and opens up perspectives on distributed manufacturing.