Apparent versus true gene expression changes of three hypoxia-related genes in autopsy derived tissue and the importance of normalisation.

The aim of our work was to show how a chosen normal-isation strategy can affect the outcome of quantitative gene expression studies. As an example, we analysed the expression of three genes known to be upregulated under hypoxic conditions: HIF1A, VEGF and SLC2A1 (GLUT1). Raw RT-qPCR data were normalised using two different strategies: a straightforward normalisation against a single reference gene, GAPDH, using the 2(-ΔΔCt) algorithm and a more complex normalisation against a normalisation factor calculated from the quantitative raw data from four previously validated reference genes. We found that the two different normalisation strategies revealed contradicting results: normalising against a validated set of reference genes revealed an upregulation of the three genes of interest in three post-mortem tissue samples (cardiac muscle, skeletal muscle and brain) under hypoxic conditions. Interestingly, we found a statistically significant difference in the relative transcript abundance of VEGF in cardiac muscle between donors who died of asphyxia versus donors who died from cardiac death. Normalisation against GAPDH alone revealed no upregulation but, in some instances, a downregulation of the genes of interest. To further analyse this discrepancy, the stability of all reference genes used were reassessed and the very low expression stability of GAPDH was found to originate from the co-regulation of this gene under hypoxic conditions. We concluded that GAPDH is not a suitable reference gene for the quantitative analysis of gene expression in hypoxia and that validation of reference genes is a crucial step for generating biologically meaningful data.

Temporal evolution of MR perfusion in neonatal hypoxic-ischemic encephalopathy

PURPOSE: To illustrate the evolution of brain perfusion-weighted magnetic resonance imaging (PWI-MRI) in severe neonatal hypoxic-ischemic (HI) encephalopathy, and its possible relation to further neurodevelopmental outcome. MATERIALS AND METHODS: Two term neonates with HI encephalopathy underwent an early and a late MRI, including PWI. They were followed until eight months of age. A total of three "normal controls" were also included. Perfusion maps were obtained, and relative cerebral blood flow (rCBF) and cerebral blood volume (rCBV) values were measured. RESULTS: Compared to normal neonates, a hyperperfusion (increased rCBF and rCBV) was present on early scans in the whole brain. On late scans, hyperperfusion persisted in cortical gray matter (normalization of rCBF and rCBV ratios in white matter and basal ganglia, but not in cortical gray matter). Diffusion-weighted imaging (DWI) was normalized, and extensive lesions became visible on T2-weighted images. Both patients displayed very abnormal outcome: Patient 2 with the more abnormal early and late hyperperfusion being the worst. CONCLUSION: PWI in HI encephalopathy did not have the same temporal evolution as DWI, and remained abnormal for more than one week after injury. This could be a marker of an ongoing mechanism underlying severe neonatal HI encephalopathy. Evolution of PWI might help to predict further neurodevelopmental outcome.

Mutations in B3GALT6, which encodes a glycosaminoglycan linker region enzyme, cause a spectrum of skeletal and connective tissue disorders.

Proteoglycans (PGs) are a major component of the extracellular matrix in many tissues and function as structural and regulatory molecules. PGs are composed of core proteins and glycosaminoglycan (GAG) side chains. The biosynthesis of GAGs starts with the linker region that consists of four sugar residues and is followed by repeating disaccharide units. By exome sequencing, we found that B3GALT6 encoding an enzyme involved in the biosynthesis of the GAG linker region is responsible for a severe skeletal dysplasia, spondyloepimetaphyseal dysplasia with joint laxity type 1 (SEMD-JL1). B3GALT6 loss-of-function mutations were found in individuals with SEMD-JL1 from seven families. In a subsequent candidate gene study based on the phenotypic similarity, we found that B3GALT6 is also responsible for a connective tissue disease, Ehlers-Danlos syndrome (progeroid form). Recessive loss-of-function mutations in B3GALT6 result in a spectrum of disorders affecting a broad range of skeletal and connective tissues characterized by lax skin, muscle hypotonia, joint dislocation, and spinal deformity. The pleiotropic phenotypes of the disorders indicate that B3GALT6 plays a critical role in a wide range of biological processes in various tissues, including skin, bone, cartilage, tendon, and ligament.

Ectopic lymphoid-organ development occurs through interleukin 7-mediated enhanced survival of lymphoid-tissue-inducer cells.

Development of Peyer's patches and lymph nodes requires the interaction between CD4+ CD3- IL-7Ralpha+ lymphoid-tissue inducer (LTi) and VCAM-1+ organizer cells. Here we showed that by promoting their survival, enhanced expression of interleukin-7 (IL-7) in transgenic mice resulted in accumulation of LTi cells. With increased IL-7 availability, de novo formation of VCAM-1+ Peyer's patch anlagen occurred along the entire fetal gut resulting in a 5-fold increase in Peyer's patch numbers. IL-7 overexpression also led to formation of multiple organized ectopic lymph nodes and cecal patches. After immunization, ectopic lymph nodes developed normal T cell-dependent B cell responses and germinal centers. Mice overexpressing IL-7 but lacking either RORgamma, a factor required for LTi cell generation, or lymphotoxin alpha1beta2 had neither Peyer's patches nor ectopic lymph nodes. Therefore, by controlling LTi cell numbers, IL-7 can regulate the formation of both normal and ectopic lymphoid organs.

Blockade of the renin-angiotensin system and renal tissue oxygenation as measured with BOLD-MRI in patients with type 2 diabetes.

AIM: To assess whether blockade of the renin-angiotensin system (RAS), a recognized strategy to prevent the progression of diabetic nephropathy, affects renal tissue oxygenation in type 2 diabetes mellitus (T2DM) patients. METHODS: Prospective randomized 2-way cross over study; T2DM patients with (micro)albuminuria and/or hypertension underwent blood oxygenation level-dependent magnetic resonance imaging (BOLD-MRI) at baseline, after one month of enalapril (20mgqd), and after one month of candesartan (16mgqd). Each BOLD-MRI was performed before and after the administration of furosemide. The mean R2* (=1/T2*) values in the medulla and cortex were calculated, a low R2* indicating high tissue oxygenation. RESULTS: Twelve patients (mean age: 60±11 years, eGFR: 62±22ml/min/1.73m(2)) completed the study. Neither chronic enalapril nor candesartan intake modified renal cortical or medullary R2* levels. Furosemide significantly decreased cortical and medullary R2* levels suggesting a transient increase in renal oxygenation. Medullary R2* levels correlated positively with urinary sodium excretion and systemic blood pressure, suggesting lower renal oxygenation at higher dietary sodium intake and blood pressure; cortical R2* levels correlated positively with glycemia and HbA1c. CONCLUSION: RAS blockade does not seem to increase renal tissue oxygenation in T2DM hypertensive patients. The response to furosemide and the association with 24h urinary sodium excretion emphasize the crucial role of renal sodium handling as one of the main determinants of renal tissue oxygenation.

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.

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.

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 diagnostic performance of perfusion-cardiovascular magnetic resonance versus SPECT to detect coronary artery disease: The secondary endpoints of the multicenter multivendor MR-IMPACT II (Magnetic Resonance Imaging for Myocardial Perfusion Assessment in Coronary Artery Disease Trial).

ABSTRACT: BACKGROUND: Perfusion-cardiovascular magnetic resonance (CMR) is generally accepted as an alternative to SPECT to assess myocardial ischemia non-invasively. However its performance vs gated-SPECT and in sub-populations is not fully established. The goal was to compare in a multicenter setting the diagnostic performance of perfusion-CMR and gated-SPECT for the detection of CAD in various populations using conventional x-ray coronary angiography (CXA) as the standard of reference. METHODS: In 33 centers (in US and Europe) 533 patients, eligible for CXA or SPECT, were enrolled in this multivendor trial. SPECT and CXA were performed within 4 weeks before or after CMR in all patients. Prevalence of CAD in the sample was 49% and 515 patients received MR contrast medium. Drop-out rates for CMR and SPECT were 5.6% and 3.7%, respectively (ns). The study was powered for the primary endpoint of non-inferiority of CMR vs SPECT for both, sensitivity and specificity for the detection of CAD (using a single-threshold reading), the results for the primary endpoint were reported elsewhere. In this article secondary endpoints are presented, i.e. the diagnostic performance of CMR versus SPECT in subpopulations such as multi-vessel disease (MVD), in men, in women, and in patients without prior myocardial infarction (MI). For diagnostic performance assessment the area under the receiver-operator-characteristics-curve (AUC) was calculated. Readers were blinded versus clinical data, CXA, and imaging results. RESULTS: The diagnostic performance (= area under ROC = AUC) of CMR was superior to SPECT (p = 0.0004, n = 425) and to gated-SPECT (p = 0.018, n = 253). CMR performed better than SPECT in MVD (p = 0.003 vs all SPECT, p = 0.04 vs gated-SPECT), in men (p = 0.004, n = 313) and in women (p = 0.03, n = 112) as well as in the non-infarct patients (p = 0.005, n = 186 in 1-3 vessel disease and p = 0.015, n = 140 in MVD). CONCLUSION: In this large multicenter, multivendor study the diagnostic performance of perfusion-CMR to detect CAD was superior to perfusion SPECT in the entire population and in sub-groups. Perfusion-CMR can be recommended as an alternative for SPECT imaging. TRIAL REGISTRATION: ClinicalTrials.gov, Identifier: NCT00977093.

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.