Mechanical constraints imposed by 3D cellular geometry and arrangement modulate growth patterns in the Arabidopsis embryo

Morphogenesis occurs in 3D space over time and is guided by coordinated gene expression programs. Here we use postembryonic development in Arabidopsis plants to investigate the genetic control of growth. We demonstrate that gene expression driving the production of the growth-stimulating hormone gibberellic acid and downstream growth factors is first induced within the radicle tip of the embryo. The center of cell expansion is, however, spatially displaced from the center of gene expression. Because the rapidly growing cells have very different geometry from that of those at the tip, we hypothesized that mechanical factors may contribute to this growth displacement. To this end we developed 3D finite-element method models of growing custom-designed digital embryos at cellular resolution. We used this framework to conceptualize how cell size, shape, and topology influence tissue growth and to explore the interplay of geometrical and genetic inputs into growth distribution. Our simulations showed that mechanical constraints are sufficient to explain the disconnect between the experimentally observed spatiotemporal patterns of gene expression and early postembryonic growth. The center of cell expansion is the position where genetic and mechanical facilitators of growth converge. We have thus uncovered a mechanism whereby 3D cellular geometry helps direct where genetically specified growth takes place.

Geometry and kinematics of the Roisan-Cignana Shear Zone, and the orogenic evolution of the Dent Blanche Tectonic System (Western Alps)

The Dent Blanche Tectonic System (DBTS) is a composite thrust sheet derived from the previously thinned passive Adriatic continental margin. A kilometric high-strain zone, the Roisan-Cignana Shear Zone (RCSZ) defines the major tectonic boundary within the DBTS and separates it into two subunits, the Dent Blanche s.s. nappe to the northwest and the Mont Mary nappe to the southeast. Within this shear zone, tectonic slices of Mesozoic and pre-Alpine meta-sediments became amalgamated with continental basement rocks of the Adriatic margin. The occurrence of high pressure assemblages along the contact between these tectonic slices indicates that the amalgamation occurred prior to or during the subduction process, at an early stage of the Alpine orogenic cycle. Detailed mapping, petrographic and structural analysis show that the Roisan-Cignana Shear Zone results from several superimposed Alpine structural and metamorphic stages. Subduction of the continental fragments is recorded by blueschist-facies deformation, whereas the Alpine collision is reflected by a greenschist facies overprint associated with the development of large-scale open folds. The postnappe evolution comprises the development of low-angle brittle faults, followed by large-scale folding (Vanzone phase) and finally brittle extensional faults. The RCSZ shows that fragments of continental crust had been torn off the passive continental margin prior to continental collision, thus recording the entire history of the orogenic cycle. The role of preceding Permo-Triassic lithospheric thinning, Jurassic rifting, and ablative subduction processes in controlling the removal of crustal fragments from the reactivated passive continental margin is discussed. Results of this study constrain the temporal sequence of the tectono-metamorphic processes involved in the assembly of the DBTS, but they also show limits on the interpretation. In particular it remains difficult to judge to what extent precollisional rifting at the Adriatic continental margin preconditioned the efficiency of convergent processes, i.e. accretion, subduction, and orogenic exhumation.

SrGAP2-Dependent Integration of Membrane Geometry and Slit-Robo-Repulsive Cues Regulates Fibroblast Contact Inhibition of Locomotion

Migrating fibroblasts undergo contact inhibition of locomotion (CIL), a process that was discovered five decades ago and still is not fully understood at the molecular level. We identify the Slit2-Robo4-srGAP2 signaling network as a key regulator of CIL in fibroblasts. CIL involves highly dynamic contact protrusions with a specialized actin cytoskeleton that stochastically explore cell-cell overlaps between colliding fibroblasts. A membrane curvature-sensing F-BAR domain pre-localizes srGAP2 to protruding edges and terminates their extension phase in response to cell collision. A FRET-based biosensor reveals that Rac1 activity is focused in a band at the tip of contact protrusions, in contrast to the broad activation gradient in contact-free protrusions. SrGAP2 specifically controls the duration of Rac1 activity in contact protrusions, but not in contact-free protrusions. We propose that srGAP2 integrates cell edge curvature and Slit-Robo-mediated repulsive cues to fine-tune Rac1 activation dynamics in contact protrusions to spatiotemporally coordinate CIL.

Validation of the Swiss methane emission inventory by atmospheric observations and inverse modelling

Atmospheric inverse modelling has the potential to provide observation-based estimates of greenhouse gas emissions at the country scale, thereby allowing for an independent validation of national emission inventories. Here, we present a regional-scale inverse modelling study to quantify the emissions of methane (CH₄) from Switzerland, making use of the newly established CarboCount-CH measurement network and a high-resolution Lagrangian transport model. In our reference inversion, prior emissions were taken from the "bottom-up" Swiss Greenhouse Gas Inventory (SGHGI) as published by the Swiss Federal Office for the Environment in 2014 for the year 2012. Overall we estimate national CH₄ emissions to be 196 ± 18 Gg yr⁻¹ for the year 2013 (1σ uncertainty). This result is in close agreement with the recently revised SGHGI estimate of 206 ± 33 Gg yr⁻¹ as reported in 2015 for the year 2012. Results from sensitivity inversions using alternative prior emissions, uncertainty covariance settings, large-scale background mole fractions, two different inverse algorithms (Bayesian and extended Kalman filter), and two different transport models confirm the robustness and independent character of our estimate. According to the latest SGHGI estimate the main CH₄ source categories in Switzerland are agriculture (78 %), waste handling (15 %) and natural gas distribution and combustion (6 %). The spatial distribution and seasonal variability of our posterior emissions suggest an overestimation of agricultural CH₄ emissions by 10 to 20 % in the most recent SGHGI, which is likely due to an overestimation of emissions from manure handling. Urban areas do not appear as emission hotspots in our posterior results, suggesting that leakages from natural gas distribution are only a minor source of CH₄ in Switzerland. This is consistent with rather low emissions of 8.4 Gg yr⁻¹ reported by the SGHGI but inconsistent with the much higher value of 32 Gg yr⁻¹ implied by the EDGARv4.2 inventory for this sector. Increased CH₄ emissions (up to 30 % compared to the prior) were deduced for the north-eastern parts of Switzerland. This feature was common to most sensitivity inversions, which is a strong indicator that it is a real feature and not an artefact of the transport model and the inversion system. However, it was not possible to assign an unambiguous source process to the region. The observations of the CarboCount-CH network provided invaluable and independent information for the validation of the national bottom-up inventory. Similar systems need to be sustained to provide independent monitoring of future climate agreements.

Computational fluid dynamics modelling in cardiovascular medicine.

This paper reviews the methods, benefits and challenges associated with the adoption and translation of computational fluid dynamics (CFD) modelling within cardiovascular medicine. CFD, a specialist area of mathematics and a branch of fluid mechanics, is used routinely in a diverse range of safety-critical engineering systems, which increasingly is being applied to the cardiovascular system. By facilitating rapid, economical, low-risk prototyping, CFD modelling has already revolutionised research and development of devices such as stents, valve prostheses, and ventricular assist devices. Combined with cardiovascular imaging, CFD simulation enables detailed characterisation of complex physiological pressure and flow fields and the computation of metrics which cannot be directly measured, for example, wall shear stress. CFD models are now being translated into clinical tools for physicians to use across the spectrum of coronary, valvular, congenital, myocardial and peripheral vascular diseases. CFD modelling is apposite for minimally-invasive patient assessment. Patient-specific (incorporating data unique to the individual) and multi-scale (combining models of different length- and time-scales) modelling enables individualised risk prediction and virtual treatment planning. This represents a significant departure from traditional dependence upon registry-based, population-averaged data. Model integration is progressively moving towards 'digital patient' or 'virtual physiological human' representations. When combined with population-scale numerical models, these models have the potential to reduce the cost, time and risk associated with clinical trials. The adoption of CFD modelling signals a new era in cardiovascular medicine. While potentially highly beneficial, a number of academic and commercial groups are addressing the associated methodological, regulatory, education- and service-related challenges.

Monte Carlo Dose Calculation on Deforming Anatomy

This article presents the implementation and validation of a dose calculation approach for deforming anatomical objects. Deformation is represented by deformation vector fields leading to deformed voxel grids representing the different deformation scenarios. Particle transport in the resulting deformed voxels is handled through the approximation of voxel surfaces by triangles in the geometry implementation of the Swiss Monte Carlo Plan framework. The focus lies on the validation methodology which uses computational phantoms representing the same physical object through regular and irregular voxel grids. These phantoms are chosen such that the new implementation for a deformed voxel grid can be compared directly with an established dose calculation algorithm for regular grids. Furthermore, separate validation of the aspects voxel geometry and the density changes resulting from deformation is achieved through suitable design of the validation phantom. We show that equivalent results are obtained with the proposed method and that no statistically significant errors are introduced through the implementation for irregular voxel geometries. This enables the use of the presented and validated implementation for further investigations of dose calculation on deforming anatomy.

Screening for Chlamydia trachomatis: a systematic review of the economic evaluations and modelling

OBJECTIVE: To review systematically and critically, evidence used to derive estimates of costs and cost effectiveness of chlamydia screening. METHODS: Systematic review. A search of 11 electronic bibliographic databases from the earliest date available to August 2004 using keywords including chlamydia, pelvic inflammatory disease, economic evaluation, and cost. We included studies of chlamydia screening in males and/or females over 14 years, including studies of diagnostic tests, contact tracing, and treatment as part of a screening programme. Outcomes included cases of chlamydia identified and major outcomes averted. We assessed methodological quality and the modelling approach used. RESULTS: Of 713 identified papers we included 57 formal economic evaluations and two cost studies. Most studies found chlamydia screening to be cost effective, partner notification to be an effective adjunct, and testing with nucleic acid amplification tests, and treatment with azithromycin to be cost effective. Methodological problems limited the validity of these findings: most studies used static models that are inappropriate for infectious diseases; restricted outcomes were used as a basis for policy recommendations; and high estimates of the probability of chlamydia associated complications might have overestimated cost effectiveness. Two high quality dynamic modelling studies found opportunistic screening to be cost effective but poor reporting or uncertainty about complication rates make interpretation difficult. CONCLUSION: The inappropriate use of static models to study interventions to prevent a communicable disease means that uncertainty remains about whether chlamydia screening programmes are cost effective or not. The results of this review can be used by health service managers in the allocation of resources, and health economists and other researchers who are considering further research in this area.

A Closed-Form, Consistent and Robust Solution to Uncalibrated Photometric Stereo Via Local Diffuse Reflectance Maxima

Images of an object under different illumination are known to provide strong cues about the object surface. A mathematical formalization of how to recover the normal map of such a surface leads to the so-called uncalibrated photometric stereo problem. In the simplest instance, this problem can be reduced to the task of identifying only three parameters: the so-called generalized bas-relief (GBR) ambiguity. The challenge is to find additional general assumptions about the object, that identify these parameters uniquely. Current approaches are not consistent, i.e., they provide different solutions when run multiple times on the same data. To address this limitation, we propose exploiting local diffuse reflectance (LDR) maxima, i.e., points in the scene where the normal vector is parallel to the illumination direction (see Fig. 1). We demonstrate several noteworthy properties of these maxima: a closed-form solution, computational efficiency and GBR consistency. An LDR maximum yields a simple closed-form solution corresponding to a semi-circle in the GBR parameters space (see Fig. 2); because as few as two diffuse maxima in different images identify a unique solution, the identification of the GBR parameters can be achieved very efficiently; finally, the algorithm is consistent as it always returns the same solution given the same data. Our algorithm is also remarkably robust: It can obtain an accurate estimate of the GBR parameters even with extremely high levels of outliers in the detected maxima (up to 80 % of the observations). The method is validated on real data and achieves state-of-the-art results.

Solute transport in crystalline rocks at Äspö — I: Geological basis and model calibration

Water-conducting faults and fractures were studied in the granite-hosted A¨ spo¨ Hard Rock Laboratory (SE Sweden). On a scale of decametres and larger, steeply dipping faults dominate and contain a variety of different fault rocks (mylonites, cataclasites, fault gouges). On a smaller scale, somewhat less regular fracture patterns were found. Conceptual models of the fault and fracture geometries and of the properties of rock types adjacent to fractures were derived and used as input for the modelling of in situ dipole tracer tests that were conducted in the framework of the Tracer Retention Understanding Experiment (TRUE-1) on a scale of metres. After the identification of all relevant transport and retardation processes, blind predictions of the breakthroughs of conservative to moderately sorbing tracers were calculated and then compared with the experimental data. This paper provides the geological basis and model calibration, while the predictive and inverse modelling work is the topic of the companion paper [J. Contam. Hydrol. 61 (2003) 175]. The TRUE-1 experimental volume is highly fractured and contains the same types of fault rocks and alterations as on the decametric scale. The experimental flow field was modelled on the basis of a 2D-streamtube formalism with an underlying homogeneous and isotropic transmissivity field. Tracer transport was modelled using the dual porosity medium approach, which is linked to the flow model by the flow porosity. Given the substantial pumping rates in the extraction borehole, the transport domain has a maximum width of a few centimetres only. It is concluded that both the uncertainty with regard to the length of individual fractures and the detailed geometry of the network along the flowpath between injection and extraction boreholes are not critical because flow is largely one-dimensional, whether through a single fracture or a network. Process identification and model calibration were based on a single uranine breakthrough (test PDT3), which clearly showed that matrix diffusion had to be included in the model even over the short experimental time scales, evidenced by a characteristic shape of the trailing edge of the breakthrough curve. Using the geological information and therefore considering limited matrix diffusion into a thin fault gouge horizon resulted in a good fit to the experiment. On the other hand, fresh granite was found not to interact noticeably with the tracers over the time scales of the experiments. While fracture-filling gouge materials are very efficient in retarding tracers over short periods of time (hours–days), their volume is very small and, with time progressing, retardation will be dominated by altered wall rock and, finally, by fresh granite. In such rocks, both porosity (and therefore the effective diffusion coefficient) and sorption Kds are more than one order of magnitude smaller compared to fault gouge, thus indicating that long-term retardation is expected to occur but to be less pronounced.

DYPTOP: a cost-efficient TOPMODEL implementation to simulate sub-grid spatio-temporal dynamics of global wetlands and peatlands

Simulating the spatio-temporal dynamics of inundation is key to understanding the role of wetlands under past and future climate change. Earlier modelling studies have mostly relied on fixed prescribed peatland maps and inundation time series of limited temporal coverage. Here, we describe and assess the the Dynamical Peatland Model Based on TOPMODEL (DYPTOP), which predicts the extent of inundation based on a computationally efficient TOPMODEL implementation. This approach rests on an empirical, grid-cell-specific relationship between the mean soil water balance and the flooded area. DYPTOP combines the simulated inundation extent and its temporal persistency with criteria for the ecosystem water balance and the modelled peatland-specific soil carbon balance to predict the global distribution of peatlands. We apply DYPTOP in combination with the LPX-Bern DGVM and benchmark the global-scale distribution, extent, and seasonality of inundation against satellite data. DYPTOP successfully predicts the spatial distribution and extent of wetlands and major boreal and tropical peatland complexes and reveals the governing limitations to peatland occurrence across the globe. Peatlands covering large boreal lowlands are reproduced only when accounting for a positive feedback induced by the enhanced mean soil water holding capacity in peatland-dominated regions. DYPTOP is designed to minimize input data requirements, optimizes computational efficiency and allows for a modular adoption in Earth system models.

Computational Comparison of Continuous and Discontinuous Galerkin Time-Stepping Methods for Nonlinear Initial Value Problems

This article centers on the computational performance of the continuous and discontinuous Galerkin time stepping schemes for general first-order initial value problems in R n , with continuous nonlinearities. We briefly review a recent existence result for discrete solutions from [6], and provide a numerical comparison of the two time discretization methods.

Computer Assisted Diagnosis and Treatment Planning of Femoroacetabular Impingement (FAI)

Femoroacetabular impingement (FAI) is a dynamic conflict of the hip defined by a pathological, early abutment of the proximal femur onto the acetabulum or pelvis. In the past two decades, FAI has received increasing focus in both research and clinical practice as a cause of hip pain and prearthrotic deformity. Anatomical abnormalities such as an aspherical femoral head (cam-type FAI), a focal or general overgrowth of the acetabulum (pincer-type FAI), a high riding greater or lesser trochanter (extra-articular FAI), or abnormal torsion of the femur have been identified as underlying pathomorphologies. Open and arthroscopic treatment options are available to correct the deformity and to allow impingement-free range of motion. In routine practice, diagnosis and treatment planning of FAI is based on clinical examination and conventional imaging modalities such as standard radiography, magnetic resonance arthrography (MRA), and computed tomography (CT). Modern software tools allow three-dimensional analysis of the hip joint by extracting pelvic landmarks from two-dimensional antero-posterior pelvic radiographs. An object-oriented cross-platform program (Hip2Norm) has been developed and validated to standardize pelvic rotation and tilt on conventional AP pelvis radiographs. It has been shown that Hip2Norm is an accurate, consistent, reliable and reproducible tool for the correction of selected hip parameters on conventional radiographs. In contrast to conventional imaging modalities, which provide only static visualization, novel computer assisted tools have been developed to allow the dynamic analysis of FAI pathomechanics. In this context, a validated, CT-based software package (HipMotion) has been introduced. HipMotion is based on polygonal three-dimensional models of the patient’s pelvis and femur. The software includes simulation methods for range of motion, collision detection and accurate mapping of impingement areas. A preoperative treatment plan can be created by performing a virtual resection of any mapped impingement zones both on the femoral head-neck junction, as well as the acetabular rim using the same three-dimensional models. The following book chapter provides a summarized description of current computer-assisted tools for the diagnosis and treatment planning of FAI highlighting the possibility for both static and dynamic evaluation, reliability and reproducibility, and its applicability to routine clinical use.

Reduced muscle mass and bone size in pediatric patients with inflammatory bowel disease

BACKGROUND: Decreased bone mineral density has been reported in children with inflammatory bowel disease (IBD). We used peripheral quantitative computed tomography (pQCT) to assess bone mineralization, geometry, and muscle cross-sectional area (CSA) in pediatric IBD. METHODS: In a cross-sectional study, pQCT of the forearm was applied in 143 IBD patients (mean age 13.9 +/- 3.5 years); 29% were newly diagnosed, 98 had Crohn's disease, and 45 had ulcerative colitis. Auxological data, cumulative glucocorticoid dose, disease activity indices, laboratory markers for inflammation, and bone metabolism were related to the results of pQCT. RESULTS: Patients were compromised in height (-0.82 +/- 1.1 SD), weight (-0.77 +/- 1.0 SD), muscle mass (-1.12 +/- 1.0 SD), and total bone cross-sectional area (-0.79 +/- 1.0 SD) compared to age- and sex-matched healthy controls (z-scores). In newly diagnosed patients, the ratio of bone mineral mass per muscle CSA was higher than in those with longer disease duration (1.00 versus 0.30, P = 0.007). Serum albumin level and disease activity correlated with muscle mass, accounting for 41.0% of variability in muscle mass (P < 0.01). The trabecular bone mineral density z-score was on average at the lower normal level (-0.40 +/- 1.3 SD, P < 0.05). CONCLUSIONS: Reduced bone geometry was explained only in part by reduced height. Bone disease in children with IBD seems to be secondary to muscle wasting, which is already present at diagnosis. With longer disease duration, bone adapts to the lower muscle CSA. Serum albumin concentration is a good marker for muscle wasting and abnormal bone development.

Reduced trabecular bone mineral density and cortical thickness accompanied by increased outer bone circumference in metacarpal bone of rheumatoid arthritis patients: a cross-sectional study

Introduction The objective of this study was to assess three-dimensional bone geometry and density at the epiphysis and shaft of the third meta-carpal bone of rheumatoid arthritis (RA) patients in comparison to healthy controls with the novel method of peripheral quantitative computed tomography (pQCT). Methods PQCT scans were performed in 50 female RA patients and 100 healthy female controls at the distal epiphyses and shafts of the third metacarpal bone, the radius and the tibia. Reproducibility was determined by coefficient of varia-tion. Bone densitometric and geometric parameters were compared between the two groups and correlated to disease characteristics. Results Reproducibility of different pQCT parameters was between 0.7% and 2.5%. RA patients had 12% to 19% lower trabecular bone mineral density (BMD) (P ≤ 0.001) at the distal epiphyses of radius, tibia and metacarpal bone. At the shafts of these bones RA patients had 7% to 16% thinner cortices (P ≤ 0.03). Total cross-sectional area (CSA) at the metacarpal bone shaft of pa-tients was larger (between 5% and 7%, P < 0.02), and relative cortical area was reduced by 13%. Erosiveness by Ratingen score correlated negatively with tra-becular and total BMD at the epiphyses and shaft cortical thickness of all measured bones (P < 0.04). Conclusions Reduced trabecular BMD and thinner cortices at peripheral bones, and a greater bone shaft diameter at the metacarpal bone suggest RA spe-cific bone alterations. The proposed pQCT protocol is reliable and allows measuring juxta-articular trabecular BMD and shaft geometry at the metacarpal bone.

The binding mode of side chain- and C3-modified epothilones to tubulin

The tubulin-binding mode of C3- and C15-modified analogues of epothilone A (Epo A) was determined by NMR spectroscopy and computational methods and compared with the existing structural models of tubulin-bound natural Epo A. Only minor differences were observed in the conformation of the macrocycle between Epo A and the C3-modified analogues investigated. In particular, 3-deoxy- (compound 2) and 3-deoxy-2,3-didehydro-Epo A (3) were found to adopt similar conformations in the tubulin-binding cleft as Epo A, thus indicating that the 3-OH group is not essential for epothilones to assume their bioactive conformation. None of the available models of the tubulin-epothilone complex is able to fully recapitulate the differences in tubulin-polymerizing activity and microtubule-binding affinity between C20-modified epothilones 6 (C20-propyl), 7 (C20-butyl), and 8 (C20-hydroxypropyl). Based on the results of transferred NOE experiments in the presence of tubulin, the isomeric C15 quinoline-based Epo B analogues 4 and 5 show very similar orientations of the side chain, irrespective of the position of the nitrogen atom in the quinoline ring. The quinoline side chain stacks on the imidazole moiety of beta-His227 with equal efficiency in both cases, thus suggesting that the aromatic side chain moiety in epothilones contributes to tubulin binding through strong van der Waals interactions with the protein rather than hydrogen bonding involving the heteroaromatic nitrogen atom. These conclusions are in line with existing tubulin polymerization and microtubule-binding data for 4, 5, and Epo B.