Comparison and validation of tissue modelization and statistical classification methods in T1-weighted MR brain images.

This paper presents a validation study on statistical nonsupervised brain tissue classification techniques in magnetic resonance (MR) images. Several image models assuming different hypotheses regarding the intensity distribution model, the spatial model and the number of classes are assessed. The methods are tested on simulated data for which the classification ground truth is known. Different noise and intensity nonuniformities are added to simulate real imaging conditions. No enhancement of the image quality is considered either before or during the classification process. This way, the accuracy of the methods and their robustness against image artifacts are tested. Classification is also performed on real data where a quantitative validation compares the methods' results with an estimated ground truth from manual segmentations by experts. Validity of the various classification methods in the labeling of the image as well as in the tissue volume is estimated with different local and global measures. Results demonstrate that methods relying on both intensity and spatial information are more robust to noise and field inhomogeneities. We also demonstrate that partial volume is not perfectly modeled, even though methods that account for mixture classes outperform methods that only consider pure Gaussian classes. Finally, we show that simulated data results can also be extended to real data.

The JNK inhibitor XG-102 protects from ischemic damage with delayed intravenous administration also in the presence of recombinant tissue plasminogen activator.

BACKGROUND: XG-102 (formerly D-JNKI1), a TAT-coupled dextrogyre peptide which selectively inhibits the c-Jun N-terminal kinase, is a powerful neuroprotectant in mouse models of middle cerebral artery occlusion (MCAo) with delayed intracerebroventricular injection. We aimed to determine whether this neuroprotection could also be achieved by intravenous injection of XG-102, which is a more feasible approach for future use in stroke patients. We also tested the compatibility of the compound with recombinant tissue plasminogen activator (rtPA), commonly used for intravenous thrombolysis and known to enhance excitotoxicity. METHODS: Male ICR-CD1 mice were subjected to a 30-min-suture MCAo. XG-102 was injected intravenously in a single dose, 6 h after ischemia. Hippocampal slice cultures were subjected to oxygen (5%) and glucose (1 mM) deprivation for 30 min. rtPA was added after ischemia and before XG-102 administration, both in vitro and in vivo. RESULTS: The lowest intravenous dose achieving neuroprotection was 0.0003 mg/kg, which reduced the infarct volume after 48 h from 62 +/- 19 mm(3) (n = 18) for the vehicle-treated group to 18 +/- 9 mm(3) (n = 5, p < 0.01). The behavioral outcome was also significantly improved at two doses. Addition of rtPA after ischemia enhanced the ischemic damage both in vitro and in vivo, but XG-102 was still able to induce a significant neuroprotection. CONCLUSIONS: A single intravenous administration of XG-102 several hours after ischemia induces a powerful neuroprotection. XG-102 protects from ischemic damage in the presence of rtPA. The feasibility of systemic administration of this promising compound and its compatibility with rtPA are important steps for its development as a drug candidate in ischemic stroke.

High Molecular Weight Methacrylate Sealing of a Bridge Deck, HR-2031, Construction Report, 1988

The Iowa Department of Transportation used a high molecular weight methacrylate (HMWM) resin to seal a 3,340 ft. x 64 ft. bridge deck in October 1986. The sealing was necessary to prevent deicing salt brine from entering a substantial number of transverse cracks that coincided with the epoxy coated top steel and unprotected bottom steel. HMWM resin is a three component product composed of a monomer, a curnene hydroperoxide initiator and a cobalt naphthenate promoter. The HMWM was applied with a dual spray bar system and flat-fan nozzles. Initiated monomer delivered through one spray bar was mixed in the air with promoted monomer from the other spray bar. The application rate averaged 0.956 gallons per 100 square feet for the tined textured driving lanes. Dry sand was broadcast on the surface at an average coverage of 0.58 lbs. per square yard to maintain friction. Coring showed that the H.MWM resin penetrated the cracks more than two inches deep. Testing of the treated deck yielded Friction Numbers averaging 33 with a treaded tire compared to 36 prior to treatment. An inspection soon after treatment found five leaky cracks in one of the 15 spans. One inspection during a steady rain showed no leakage, but leakage from numerous cracks occurred during a subsequent rain. A second HMWM application was made on two spans to determine if a double application would prevent leakage. This evaluation has not been completed.

Investigation of Techniques for Accelerating the Construction of Bridge Deck Overlays, 2016

Use of bridge deck overlays is important in maximizing bridge service life. Overlays can replace the deteriorated part of the deck, thus extending the bridge life. Even though overlay construction avoids the construction of a whole new bridge deck, construction still takes significant time in re-opening the bridge to traffic. Current processes and practices are time-consuming and multiple opportunities may exist to reduce overall construction time by modifying construction requirements and/or materials utilized. Reducing the construction time could have an effect on reducing the socioeconomic costs associated with bridge deck rehabilitation and the inconvenience caused to travelers. This work included three major tasks with literature review, field investigation, and laboratory testing. Overlay concrete mix used for present construction takes long curing hours and therefore an investigation was carried out to find fast-curing concrete mixes that could reduce construction time. Several fast-cuing concrete mixes were found and suggested for further evaluation. An on-going overlay construction project was observed and documented. Through these observations, several opportunities were suggested where small modifications in the process could lead to significant time savings. With current standards of the removal depth of substrate concrete in Iowa, it takes long hours for the removal process. Four different laboratory tests were performed with different loading conditions to determine the necessary substrate concrete removal depth for a proper bond between the substrate concrete and the new overlay concrete. Several parameters, such as failure load, bond stress, and stiffness, were compared for four different concrete removal depths. Through the results and observations of this investigation several conclusions were made which could reduce bridge deck overlay construction time.

Preliminary Studies of Remedial Measures for the Prevention of Bridge Deck Deterioration, HR-146, 1970

Four series of five specimens each were investigated for static and fatigue strength. These four series differed in that there were two variables, the first being the subsidence of concrete around reinforcing bars and the second being shrinkage due to two different curing conditions. The combinations of these variables were then compared to each other by use of, load-deflection curves and S-H fatigue curves.

A Continuous Span Aluminum Girder Concrete Deck Bridge - Part 2 of 2: Fatigue Tests of Aluminum Girders, Final Report, 1997

In 1957, the Iowa State Highway Commission, with financial assistance from the aluminum industry, constructed a 220-ft (67-m) long, four-span continuous, aluminum girder bridge to carry traffic on Clive Road (86th Street) over Interstate 80 near Des Moines, Iowa. The bridge had four, welded I-shape girders that were fabricated in pairs with welded diaphragms between an exterior and an interior girder. The interior diaphragms between the girder pairs were bolted to girder brackets. A composite, reinforced concrete deck served as the roadway surface. The bridge, which had performed successfully for about 35 years of service, was removed in the fall of 1993 to make way for an interchange at the same location. Prior to the bridge demolition, load tests were conducted to monitor girder and diaphragm bending strains and deflections in the northern end span. Fatigue testing of the aluminum girders that were removed from the end spans were conducted by applying constant-amplitude, cyclic loads. These tests established the fatigue strength of an existing, welded, flange-splice detail and added, welded, flange-cover plates and horizontal web plate attachment details. This part, Part 2, of the final report focuses on the fatigue tests of the aluminum girder sections that were removed from the bridge and on the analysis of the experimental data to establish the fatigue strength of full-size specimens. Seventeen fatigue fractures that were classified as Category E weld details developed in the seven girder test specimens. Linear regression analyses of the fatigue test results established both nominal and experimental stress-range versus load cycle relationships (SN curves) for the fatigue strength of fillet-welded connections. The nominal strength SN curve obtained by this research essentially matched the SN curve for Category E aluminum weldments given in the AASHTO LRFD specifications. All of the Category E fatigue fractures that developed in the girder test specimens satisfied the allowable SN relationship specified by the fatigue provisions of the Aluminum Association. The lower-bound strength line that was set at two standard deviations below the least squares regression line through the fatigue fracture data points related well with the Aluminum Association SN curve. The results from the experimental tests of this research have provided additional information regarding behavioral characteristics of full-size, aluminum members and have confirmed that aluminum has the strength properties needed for highway bridge girders.

Study of the impact of tissue density heterogeneities on 3-dimensional abdominal dosimetry: comparison between dose kernel convolution and direct monte carlo methods.

Dose kernel convolution (DK) methods have been proposed to speed up absorbed dose calculations in molecular radionuclide therapy. Our aim was to evaluate the impact of tissue density heterogeneities (TDH) on dosimetry when using a DK method and to propose a simple density-correction method. METHODS: This study has been conducted on 3 clinical cases: case 1, non-Hodgkin lymphoma treated with (131)I-tositumomab; case 2, a neuroendocrine tumor treatment simulated with (177)Lu-peptides; and case 3, hepatocellular carcinoma treated with (90)Y-microspheres. Absorbed dose calculations were performed using a direct Monte Carlo approach accounting for TDH (3D-RD), and a DK approach (VoxelDose, or VD). For each individual voxel, the VD absorbed dose, D(VD), calculated assuming uniform density, was corrected for density, giving D(VDd). The average 3D-RD absorbed dose values, D(3DRD), were compared with D(VD) and D(VDd), using the relative difference Δ(VD/3DRD). At the voxel level, density-binned Δ(VD/3DRD) and Δ(VDd/3DRD) were plotted against ρ and fitted with a linear regression. RESULTS: The D(VD) calculations showed a good agreement with D(3DRD). Δ(VD/3DRD) was less than 3.5%, except for the tumor of case 1 (5.9%) and the renal cortex of case 2 (5.6%). At the voxel level, the Δ(VD/3DRD) range was 0%-14% for cases 1 and 2, and -3% to 7% for case 3. All 3 cases showed a linear relationship between voxel bin-averaged Δ(VD/3DRD) and density, ρ: case 1 (Δ = -0.56ρ + 0.62, R(2) = 0.93), case 2 (Δ = -0.91ρ + 0.96, R(2) = 0.99), and case 3 (Δ = -0.69ρ + 0.72, R(2) = 0.91). The density correction improved the agreement of the DK method with the Monte Carlo approach (Δ(VDd/3DRD) < 1.1%), but with a lesser extent for the tumor of case 1 (3.1%). At the voxel level, the Δ(VDd/3DRD) range decreased for the 3 clinical cases (case 1, -1% to 4%; case 2, -0.5% to 1.5%, and -1.5% to 2%). No more linear regression existed for cases 2 and 3, contrary to case 1 (Δ = 0.41ρ - 0.38, R(2) = 0.88) although the slope in case 1 was less pronounced. CONCLUSION: This study shows a small influence of TDH in the abdominal region for 3 representative clinical cases. A simple density-correction method was proposed and improved the comparison in the absorbed dose calculations when using our voxel S value implementation.

Tissue-specific opposing functions of the inflammasome adaptor ASC in the regulation of epithelial skin carcinogenesis.

A chronic inflammatory microenvironment favors tumor progression through molecular mechanisms that are still incompletely defined. In inflammation-induced skin cancers, IL-1 receptor- or caspase-1-deficient mice, or mice specifically deficient for the inflammasome adaptor protein ASC (apoptosis-associated speck-like protein containing a CARD) in myeloid cells, had reduced tumor incidence, pointing to a role for IL-1 signaling and inflammasome activation in tumor development. However, mice fully deficient for ASC were not protected, and mice specifically deficient for ASC in keratinocytes developed more tumors than controls, suggesting that, in contrast to its proinflammatory role in myeloid cells, ASC acts as a tumor-suppressor in keratinocytes. Accordingly, ASC protein expression was lost in human cutaneous squamous cell carcinoma, but not in psoriatic skin lesions. Stimulation of primary mouse keratinocytes or the human keratinocyte cell line HaCaT with UVB induced an ASC-dependent phosphorylation of p53 and expression of p53 target genes. In HaCaT cells, ASC interacted with p53 at the endogenous level upon UVB irradiation. Thus, ASC in different tissues may influence tumor growth in opposite directions: it has a proinflammatory role in infiltrating cells that favors tumor development, but it also limits keratinocyte proliferation in response to noxious stimuli, possibly through p53 activation, which helps suppressing tumors.

Tissue oxygenation and hemodynamic response to NO synthase inhibition in septic shock.

The objective of the study was to evaluate the tissue oxygenation and hemodynamic effects of NOS inhibition in clinical severe septic shock. Eight patients with septic shock refractory to volume loading and high level of adrenergic support were prospectively enrolled in the study. Increasing doses of NOS inhibitors [N(G)-nitro-L-arginine-methyl ester (L-NAME) or N(G)-monomethyl-L-arginine (L-NMMA)] were administered as i.v. bolus until a peak effect = 10 mmHg on mean blood pressure was obtained or until side effects occurred. If deemed clinically appropriate, a continuous infusion of L-NAME was instituted and adrenergic support weaning attempted. The bolus administration of NOS inhibitors transiently increased mean blood pressure by 10 mm Hg in all patients. Seven out of eight patients received an L-NAME infusion, associated over 24 h with a progressive decline in cardiac index (P < 0.001) and an increase in systemic vascular resistance (P < 0.01). Partial or total adrenergic support weaning was rapidly possible in 6/8 patients. Oxygen transport decreased (P < 0.001), but oxygen consumption remained unchanged in those patients in whom it could be measured by indirect calorimetry (5/8). Blood lactate and the difference between tonometric gastric and arterial PCO2 remained unchanged. There were 4/8 ICU survivors. We conclude that nitric oxide synthase inhibition in severe septic shock was followed with a progressive correction of the vasoplegic hemodynamic disturbances with finally normalization of cardiac output and systemic vascular resistances without any demonstrable deterioration in tissue oxygenation.

Local optima networks of hard combinatorial landscapes

Combinatorial optimization involves finding an optimal solution in a finite set of options; many everyday life problems are of this kind. However, the number of options grows exponentially with the size of the problem, such that an exhaustive search for the best solution is practically infeasible beyond a certain problem size. When efficient algorithms are not available, a practical approach to obtain an approximate solution to the problem at hand, is to start with an educated guess and gradually refine it until we have a good-enough solution. Roughly speaking, this is how local search heuristics work. These stochastic algorithms navigate the problem search space by iteratively turning the current solution into new candidate solutions, guiding the search towards better solutions. The search performance, therefore, depends on structural aspects of the search space, which in turn depend on the move operator being used to modify solutions. A common way to characterize the search space of a problem is through the study of its fitness landscape, a mathematical object comprising the space of all possible solutions, their value with respect to the optimization objective, and a relationship of neighborhood defined by the move operator. The landscape metaphor is used to explain the search dynamics as a sort of potential function. The concept is indeed similar to that of potential energy surfaces in physical chemistry. Borrowing ideas from that field, we propose to extend to combinatorial landscapes the notion of the inherent network formed by energy minima in energy landscapes. In our case, energy minima are the local optima of the combinatorial problem, and we explore several definitions for the network edges. At first, we perform an exhaustive sampling of local optima basins of attraction, and define weighted transitions between basins by accounting for all the possible ways of crossing the basins frontier via one random move. Then, we reduce the computational burden by only counting the chances of escaping a given basin via random kick moves that start at the local optimum. Finally, we approximate network edges from the search trajectory of simple search heuristics, mining the frequency and inter-arrival time with which the heuristic visits local optima. Through these methodologies, we build a weighted directed graph that provides a synthetic view of the whole landscape, and that we can characterize using the tools of complex networks science. We argue that the network characterization can advance our understanding of the structural and dynamical properties of hard combinatorial landscapes. We apply our approach to prototypical problems such as the Quadratic Assignment Problem, the NK model of rugged landscapes, and the Permutation Flow-shop Scheduling Problem. We show that some network metrics can differentiate problem classes, correlate with problem non-linearity, and predict problem hardness as measured from the performances of trajectory-based local search heuristics.

Real-time MR imaging of myocardial regional function using strain-encoding (SENC) with tissue through-plane motion tracking.

PURPOSE: To implement real-time myocardial strain-encoding (SENC) imaging in combination with tracking the tissue displacement in the through-plane direction. MATERIALS AND METHODS: SENC imaging was combined with the slice-following technique by implementing three-dimensional (3D) selective excitation. Certain adjustments were implemented to reduce scan time to one heartbeat. A total of 10 volunteers and five pigs were scanned on a 3T MRI scanner. Spatial modulation of magnetization (SPAMM)-tagged images were acquired on planes orthogonal to the SENC planes for comparison. Myocardial infarction (MI) was induced in two pigs and the resulting SENC images were compared to standard delayed-enhancement (DE) images. RESULTS: The strain values computed from SENC imaging with slice-following showed significant difference from those acquired without slice-following, especially during systole (P < 0.01). The strain curves computed from the SENC images with and without slice-following were similar to those computed from the orthogonal SPAMM images, with and without, respectively, tracking the tag line displacement in the strain direction. The resulting SENC images showed good agreement with the DE images in identifying MI in infarcted pigs. CONCLUSION: Correction of through-plane motion in real-time cardiac functional imaging is feasible using slice-following. The strain measurements are more accurate than conventional SENC measurements in humans and animals, as validated with conventional MRI tagging.

Superior flow for bridge to life with self-expanding venous cannulas.

BACKGROUND: Recently, a compact cardiopulmonary support (CPS) system designed for quick set-up for example, during emergency cannulation, has been introduced. Traditional rectilinear percutaneous cannulas are standard for remote vascular access with the original design. The present study was designed to assess the potential of performance increase by the introduction of next-generation, self-expanding venous cannulas, which can take advantage of the luminal width of the venous vasculature despite a relatively small access orifice. METHODS: Veno-arterial bypass was established in three bovine experiments (69+/-10 kg). The Lifebridge (Lifebridge GmbH, Munich, Germany) system was connected to the right atrium in a trans-jugular fashion with various venous cannulas; and the oxygenated blood was returned through the carotid artery with a 17 F percutaneous cannula. Two different venous cannulas were studied, and the correlation between the centrifugal pump speed (1500-3900 RPM), flow and the required negative pressure on the venous side was established: (A) Biomedicus 19 F (Medtronic, Tolochenaz, Switzerland); (B) Smart canula 18 F/36 F (Smartcanula LLC, Lausanne, Switzerland). RESULTS: At 1500 RPM, the blood flow was 0.44+/-0.26 l min(-1) for the 19 F rectilinear cannula versus 0.73+/-0.34 l min(-1) for the 18/36 F self-expanding cannula. At 2500 RPM the blood flow was 1.63+/-0.62 l min(-1) for the 19F rectilinear cannula versus 2.13+/-0.34 l min(-1) for the 18/36 F self-expanding cannula. At 3500 RPM, the blood flow was 2.78+/-0.47 l min(-1) for the 19 F rectilinear cannula versus 3.64+/-0.39 l min(-1) for the 18/36 F self-expanding cannula (p<0.01 for 18/36 F vs 19 F). At 1500 RPM, the venous line pressure was 18+/-8 mmHg for the 19F rectilinear cannula versus 19+/-5 mmHg for the 18/36 F self-expanding cannula. At 2500 RPM the venous line pressure accounted for -22+/-32 mmHg for the 19 F rectilinear cannula versus 2+/-5 mmHg for the 18/36 F self-expanding cannula. At 3500 RPM, the venous line pressure was -112+/-42 mmHg for the rectilinear cannula versus 28+/-7 mmHg for the 18/36 F self-expanding cannula (p<0.01 for 18 F/36 F vs 19 F). Conclusions: The negative pressure required to achieve adequate venous drainage with the self-expanding venous cannula accounts for approximately 31% of the pressure necessary with the 19 F rectilinear cannula. In addition, a pump flow of more than 4 l min(-1) can be achieved with the self-expanding design and a well-accepted negative inlet pressure for minimal blood trauma of less than 50 mmHg.

Multi-scale systems analysis of plant development : beyond cellular signaling and tissue morphodynamics

Les plantes sont essentielles pour les sociétés humaines. Notre alimentation quotidienne, les matériaux de constructions et les sources énergétiques dérivent de la biomasse végétale. En revanche, la compréhension des multiples aspects développementaux des plantes est encore peu exploitée et représente un sujet de recherche majeur pour la science. L'émergence des technologies à haut débit pour le séquençage de génome à grande échelle ou l'imagerie de haute résolution permet à présent de produire des quantités énormes d'information. L'analyse informatique est une façon d'intégrer ces données et de réduire la complexité apparente vers une échelle d'abstraction appropriée, dont la finalité est de fournir des perspectives de recherches ciblées. Ceci représente la raison première de cette thèse. En d'autres termes, nous appliquons des méthodes descriptives et prédictives combinées à des simulations numériques afin d'apporter des solutions originales à des problèmes relatifs à la morphogénèse à l'échelle de la cellule et de l'organe. Nous nous sommes fixés parmi les objectifs principaux de cette thèse d'élucider de quelle manière l'interaction croisée des phytohormones auxine et brassinosteroïdes (BRs) détermine la croissance de la cellule dans la racine du méristème apical d'Arabidopsis thaliana, l'organisme modèle de référence pour les études moléculaires en plantes. Pour reconstruire le réseau de signalement cellulaire, nous avons extrait de la littérature les informations pertinentes concernant les relations entre les protéines impliquées dans la transduction des signaux hormonaux. Le réseau a ensuite été modélisé en utilisant un formalisme logique et qualitatif pour pallier l'absence de données quantitatives. Tout d'abord, Les résultats ont permis de confirmer que l'auxine et les BRs agissent en synergie pour contrôler la croissance de la cellule, puis, d'expliquer des observations phénotypiques paradoxales et au final, de mettre à jour une interaction clef entre deux protéines dans la maintenance du méristème de la racine. Une étude ultérieure chez la plante modèle Brachypodium dystachion (Brachypo- dium) a révélé l'ajustement du réseau d'interaction croisée entre auxine et éthylène par rapport à Arabidopsis. Chez ce dernier, interférer avec la biosynthèse de l'auxine mène à la formation d'une racine courte. Néanmoins, nous avons isolé chez Brachypodium un mutant hypomorphique dans la biosynthèse de l'auxine qui affiche une racine plus longue. Nous avons alors conduit une analyse morphométrique qui a confirmé que des cellules plus anisotropique (plus fines et longues) sont à l'origine de ce phénotype racinaire. Des analyses plus approfondies ont démontré que la différence phénotypique entre Brachypodium et Arabidopsis s'explique par une inversion de la fonction régulatrice dans la relation entre le réseau de signalisation par l'éthylène et la biosynthèse de l'auxine. L'analyse morphométrique utilisée dans l'étude précédente exploite le pipeline de traitement d'image de notre méthode d'histologie quantitative. Pendant la croissance secondaire, la symétrie bilatérale de l'hypocotyle est remplacée par une symétrie radiale et une organisation concentrique des tissus constitutifs. Ces tissus sont initialement composés d'une douzaine de cellules mais peuvent aisément atteindre des dizaines de milliers dans les derniers stades du développement. Cette échelle dépasse largement le seuil d'investigation par les moyens dits 'traditionnels' comme l'imagerie directe de tissus en profondeur. L'étude de ce système pendant cette phase de développement ne peut se faire qu'en réalisant des coupes fines de l'organe, ce qui empêche une compréhension des phénomènes cellulaires dynamiques sous-jacents. Nous y avons remédié en proposant une stratégie originale nommée, histologie quantitative. De fait, nous avons extrait l'information contenue dans des images de très haute résolution de sections transverses d'hypocotyles en utilisant un pipeline d'analyse et de segmentation d'image à grande échelle. Nous l'avons ensuite combiné avec un algorithme de reconnaissance automatique des cellules. Cet outil nous a permis de réaliser une description quantitative de la progression de la croissance secondaire révélant des schémas développementales non-apparents avec une inspection visuelle classique. La formation de pôle de phloèmes en structure répétée et espacée entre eux d'une longueur constante illustre les bénéfices de notre approche. Par ailleurs, l'exploitation approfondie de ces résultats a montré un changement de croissance anisotropique des cellules du cambium et du phloème qui semble en phase avec l'expansion du xylème. Combinant des outils génétiques et de la modélisation biomécanique, nous avons démontré que seule la croissance plus rapide des tissus internes peut produire une réorientation de l'axe de croissance anisotropique des tissus périphériques. Cette prédiction a été confirmée par le calcul du ratio des taux de croissance du xylème et du phloème au cours de développement secondaire ; des ratios élevés sont effectivement observés et concomitant à l'établissement progressif et tangentiel du cambium. Ces résultats suggèrent un mécanisme d'auto-organisation établi par un gradient de division méristématique qui génèrent une distribution de contraintes mécaniques. Ceci réoriente la croissance anisotropique des tissus périphériques pour supporter la croissance secondaire. - Plants are essential for human society, because our daily food, construction materials and sustainable energy are derived from plant biomass. Yet, despite this importance, the multiple developmental aspects of plants are still poorly understood and represent a major challenge for science. With the emergence of high throughput devices for genome sequencing and high-resolution imaging, data has never been so easy to collect, generating huge amounts of information. Computational analysis is one way to integrate those data and to decrease the apparent complexity towards an appropriate scale of abstraction with the aim to eventually provide new answers and direct further research perspectives. This is the motivation behind this thesis work, i.e. the application of descriptive and predictive analytics combined with computational modeling to answer problems that revolve around morphogenesis at the subcellular and organ scale. One of the goals of this thesis is to elucidate how the auxin-brassinosteroid phytohormone interaction determines the cell growth in the root apical meristem of Arabidopsis thaliana (Arabidopsis), the plant model of reference for molecular studies. The pertinent information about signaling protein relationships was obtained through the literature to reconstruct the entire hormonal crosstalk. Due to a lack of quantitative information, we employed a qualitative modeling formalism. This work permitted to confirm the synergistic effect of the hormonal crosstalk on cell elongation, to explain some of our paradoxical mutant phenotypes and to predict a novel interaction between the BREVIS RADIX (BRX) protein and the transcription factor MONOPTEROS (MP),which turned out to be critical for the maintenance of the root meristem. On the same subcellular scale, another study in the monocot model Brachypodium dystachion (Brachypodium) revealed an alternative wiring of auxin-ethylene crosstalk as compared to Arabidopsis. In the latter, increasing interference with auxin biosynthesis results in progressively shorter roots. By contrast, a hypomorphic Brachypodium mutant isolated in this study in an enzyme of the auxin biosynthesis pathway displayed a dramatically longer seminal root. Our morphometric analysis confirmed that more anisotropic cells (thinner and longer) are principally responsible for the mutant root phenotype. Further characterization pointed towards an inverted regulatory logic in the relation between ethylene signaling and auxin biosynthesis in Brachypodium as compared to Arabidopsis, which explains the phenotypic discrepancy. Finally, the morphometric analysis of hypocotyl secondary growth that we applied in this study was performed with the image-processing pipeline of our quantitative histology method. During its secondary growth, the hypocotyl reorganizes its primary bilateral symmetry to a radial symmetry of highly specialized tissues comprising several thousand cells, starting with a few dozens. However, such a scale only permits observations in thin cross-sections, severely hampering a comprehensive analysis of the morphodynamics involved. Our quantitative histology strategy overcomes this limitation. We acquired hypocotyl cross-sections from tiled high-resolution images and extracted their information content using custom high-throughput image processing and segmentation. Coupled with an automated cell type recognition algorithm, it allows precise quantitative characterization of vascular development and reveals developmental patterns that were not evident from visual inspection, for example the steady interspace distance of the phloem poles. Further analyses indicated a change in growth anisotropy of cambial and phloem cells, which appeared in phase with the expansion of xylem. Combining genetic tools and computational modeling, we showed that the reorientation of growth anisotropy axis of peripheral tissue layers only occurs when the growth rate of central tissue is higher than the peripheral one. This was confirmed by the calculation of the ratio of the growth rate xylem to phloem throughout secondary growth. High ratios are indeed observed and concomitant with the homogenization of cambium anisotropy. These results suggest a self-organization mechanism, promoted by a gradient of division in the cambium that generates a pattern of mechanical constraints. This, in turn, reorients the growth anisotropy of peripheral tissues to sustain the secondary growth.