Root growth models: towards a new generation of continuous approaches

Models of root system growth emerged in the early 1970s, and were based on mathematical representations of root length distribution in soil. The last decade has seen the development of more complex architectural models and the use of computer-intensive approaches to study developmental and environmental processes in greater detail. There is a pressing need for predictive technologies that can integrate root system knowledge, scaling from molecular to ensembles of plants. This paper makes the case for more widespread use of simpler models of root systems based on continuous descriptions of their structure. A new theoretical framework is presented that describes the dynamics of root density distributions as a function of individual root developmental parameters such as rates of lateral root initiation, elongation, mortality, and gravitropsm. The simulations resulting from such equations can be performed most efficiently in discretized domains that deform as a result of growth, and that can be used to model the growth of many interacting root systems. The modelling principles described help to bridge the gap between continuum and architectural approaches, and enhance our understanding of the spatial development of root systems. Our simulations suggest that root systems develop in travelling wave patterns of meristems, revealing order in otherwise spatially complex and heterogeneous systems. Such knowledge should assist physiologists and geneticists to appreciate how meristem dynamics contribute to the pattern of growth and functioning of root systems in the field.

Effects of oxidised low density lipoprotein on dendritic cells: a possible immunoregulatory component of the atherogenic micro-environment?

Objective: The objective of this study was to explore the relationship between low density lipoprotein (LDL) and dendritic cell (DC) activation, based upon the hypothesis that reactive oxygen species (ROS)-mediated modification of proteins that may be present in local DC microenvironments could be important as mediators of this activation. Although LDL are known to be oxidised in vivo, and taken up by macrophages during atherogenesis; their effect on DC has not been explored previously. Methods: Human DCs were prepared from peripheral blood monocytes using GM-CSF and IL-4. Plasma LDLs were isolated by sequential gradient centrifugation, oxidised in CuSO4, and oxidation arrested to yield mild, moderate and highly oxidised LDL forms. DCs exposed to these LDLs were investigated using combined phenotypic, functional (autologous T cell activation), morphological and viability assays. Results: Highly-oxidised LDL increased DC HLA-DR, CD40 and CD86 expression, corroborated by increased DC-induced T cell proliferation. Both native and oxidised LDL induced prominent DC clustering. However, high concentrations of highly-oxidised LDL inhibited DC function, due to increased DC apoptosis. Conclusions: This study supports the hypothesis that oxidised LDL are capable of triggering the transition from sentinel to messenger DC. Furthermore, the DC clustering–activation–apoptosis sequence in the presence of different LDL forms is consistent with a regulatory DC role in immunopathogenesis of atheroma. A sequence of initial accumulation of DC, increasing LDL oxidation, and DC-induced T cell activation, may explain why local breach of tolerance can occur. Above a threshold level, however, supervening DC apoptosis limits this, contributing instead to the central plaque core.

Phylogenies reveal predictive power of traditional medicine in bioprospecting

There is controversy about whether traditional medicine can guide drug discovery, and investment in ethnobotanically led research has fluctuated. One view is that traditionally used plants are not necessarily efficacious and there are no robust methods for distinguishing the ones that are most likely to be bioactive when selecting species for further testing. Here, we reconstruct a genus-level molecular phylogeny representing the 20,000 species found in the floras of three disparate biodiversity hotspots: Nepal, New Zealand and the Cape of South Africa. Borrowing phylogenetic methods from community ecology, we reveal significant clustering of the 1,500 traditionally used species, and provide a direct measure of the relatedness of the three medicinal floras. We demonstrate shared phylogenetic patterns across the floras: related plants from these regions are used to treat medical conditions in the same therapeutic areas. This strongly suggests independent discovery of plant efficacy, an interpretation corroborated by the presence of a significantly greater proportion of known bioactive species in these plant groups than in a random sample. Phylogenetic cross-cultural comparison can focus screening efforts on a subset of traditionally used plants that are richer in bioactive compounds, and could revitalise the use of traditional knowledge in bioprospecting.

Unraveling the electronic structures of low-valent naphthalene and anthracene iron complexes: x-ray, spectroscopic, and density functional theory studies

Naphthalene and anthracene transition metalates are potent reagents, but their electronic structures have remained poorly explored. A study of four Cp*-substituted iron complexes (Cp* = pentamethylcyclopentadienyl) now gives rare insight into the bonding features of such species. The highly oxygen- and water-sensitive compounds [K(18-crown- 6){Cp*Fe(η4-C10H8)}] (K1), [K(18-crown-6){Cp*Fe(η4-C14H10)}] (K2), [Cp*Fe(η4-C10H8)] (1), and [Cp*Fe(η4-C14H10)] (2) were synthesized and characterized by NMR, UV−vis, and 57Fe Mössbauer spectroscopy. The paramagnetic complexes 1 and 2 were additionally characterized by electron paramagnetic resonance (EPR) spectroscopy and magnetic susceptibility measurements. The molecular structures of complexes K1, K2, and 2 were determined by single-crystal X-ray crystallography. Cyclic voltammetry of 1 and 2 and spectroelectrochemical experiments revealed the redox properties of these complexes, which are reversibly reduced to the monoanions [Cp*Fe(η4-C10H8)]− (1−) and [Cp*Fe(η4-C14H10)]− (2−) and reversibly oxidized to the cations [Cp*Fe(η6-C10H8)]+ (1+) and [Cp*Fe(η6-C14H10)]+ (2+). Reduced orbital charges and spin densities of the naphthalene complexes 1−/0/+ and the anthracene derivatives 2−/0/+ were obtained by density functional theory (DFT) methods. Analysis of these data suggests that the electronic structures of the anions 1− and 2− are best represented by low-spin FeII ions coordinated by anionic Cp* and dianionic naphthalene and anthracene ligands. The electronic structures of the neutral complexes 1 and 2 may be described by a superposition of two resonance configurations which, on the one hand, involve a low-spin FeI ion coordinated by the neutral naphthalene or anthracene ligand L, and, on the other hand, a low-spin FeII ion coordinated to a ligand radical L•−. Our study thus reveals the redox noninnocent character of the naphthalene and anthracene ligands, which effectively stabilize the iron atoms in a low formal, but significantly higher spectroscopic oxidation state.

Importance of density-compensated temperature change for deep North Atlantic Ocean heat uptake

The efficiency with which the oceans take up heat has a significant influence on the rate of global warming. Warming of the ocean above 700 m over the past few decades has been well documented. However, most of the ocean lies below 700 m. Here we analyse observations of heat uptake into the deep North Atlantic. We find that the extratropical North Atlantic as a whole warmed by 1.45±0.5×1022 J between 1955 and 2005, but Lower North Atlantic Deep Water cooled, most likely as an adjustment from an early twentieth-century warm period. In contrast, the heat content of Upper North Atlantic Deep Water exhibited strong decadal variability. We demonstrate and quantify the importance of density-compensated temperature anomalies for long-term heat uptake into the deep North Atlantic. These anomalies form in the subpolar gyre and propagate equatorwards. High salinity in the subpolar gyre is a key requirement for this mechanism. In the past 50 years, suitable conditions have occurred only twice: first during the 1960s and again during the past decade. We conclude that heat uptake through density-compensated temperature anomalies will contribute to deep ocean heat uptake in the near term. In the longer term, the importance of this mechanism will be determined by competition between the multiple processes that influence subpolar gyre salinity in a changing climate.

Synthesis, chemical, electrochemical characterization of oxomolybdenum(V) complexes of 2-(3,5-dimethyl pyrazol-1-yl) benzothiazole ligand: crystal structure of ligand and oxomolybdenum(VI) complex and density functional theory (DFT) calculation

This work reports the ligational behavior of the neutral bidentate chelating molecule 2-(3,5-dimethyl pyrazol-1-yl) benzothiazole towards the oxomolybdenum(V) center. Both mononuclear complexes of the type (MoOX3L)-O-V and binuclear complexes of the formula (Mo2O4X2L2)-O-V (where X = Cl, Br) are isolated in the solid state. The complexes are characterized by elemental analyses, various spectroscopic techniques (UV-Vis IR), magnetic susceptibility measurement at room temperature, and cyclic voltammetry for their redox behavior at a platinum electrode in CH3CN. The mononuclear complexes (MoOX3L)-O-V are found to be paramagnetic while the binuclear complexes Mo2O4X2L2 are diamagnetic. Crystal and molecular structure of the ligand and the dioxomolybdenum complex (MoO2Br2L)-O-VI (obtained from the complex MoOBr3L during crystallization) have been solved by single crystal X-ray diffraction technique. Relevant DFT calculations of the ligand and the complex (MoO2Br2L)-O-VI are also carried out.

From ensemble forecasts to predictive distribution functions

The translation of an ensemble of model runs into a probability distribution is a common task in model-based prediction. Common methods for such ensemble interpretations proceed as if verification and ensemble were draws from the same underlying distribution, an assumption not viable for most, if any, real world ensembles. An alternative is to consider an ensemble as merely a source of information rather than the possible scenarios of reality. This approach, which looks for maps between ensembles and probabilistic distributions, is investigated and extended. Common methods are revisited, and an improvement to standard kernel dressing, called ‘affine kernel dressing’ (AKD), is introduced. AKD assumes an affine mapping between ensemble and verification, typically not acting on individual ensemble members but on the entire ensemble as a whole, the parameters of this mapping are determined in parallel with the other dressing parameters, including a weight assigned to the unconditioned (climatological) distribution. These amendments to standard kernel dressing, albeit simple, can improve performance significantly and are shown to be appropriate for both overdispersive and underdispersive ensembles, unlike standard kernel dressing which exacerbates over dispersion. Studies are presented using operational numerical weather predictions for two locations and data from the Lorenz63 system, demonstrating both effectiveness given operational constraints and statistical significance given a large sample.

Intermittently-visual tracking experiments reveal the roles of error-correction and predictive mechanisms in the human visual-motor control system

Prediction mechanism is necessary for human visual motion to compensate a delay of sensory-motor system. In a previous study, “proactive control” was discussed as one example of predictive function of human beings, in which motion of hands preceded the virtual moving target in visual tracking experiments. To study the roles of the positional-error correction mechanism and the prediction mechanism, we carried out an intermittently-visual tracking experiment where a circular orbit is segmented into the target-visible regions and the target-invisible regions. Main results found in this research were following. A rhythmic component appeared in the tracer velocity when the target velocity was relatively high. The period of the rhythm in the brain obtained from environmental stimuli is shortened more than 10%. The shortening of the period of rhythm in the brain accelerates the hand motion as soon as the visual information is cut-off, and causes the precedence of hand motion to the target motion. Although the precedence of the hand in the blind region is reset by the environmental information when the target enters the visible region, the hand motion precedes the target in average when the predictive mechanism dominates the error-corrective mechanism.

Evaluation of the use of high-density SNP genotyping to implement UPOV Model 2 for DUS testing in barley

Developments in high-throughput genotyping provide an opportunity to explore the application of marker technology in distinctness, uniformity and stability (DUS) testing of new varieties. We have used a large set of molecular markers to assess the feasibility of a UPOV Model 2 approach: “Calibration of threshold levels for molecular characteristics against the minimum distance in traditional characteristics”. We have examined 431 winter and spring barley varieties, with data from UK DUS trials comprising 28 characteristics, together with genotype data from 3072 SNP markers. Inter varietal distances were calculated and we found higher correlations between molecular and morphological distances than have been previously reported. When varieties were grouped by kinship, phenotypic and genotypic distances of these groups correlated well. We estimated the minimum marker numbers required and showed there was a ceiling after which the correlations do not improve. To investigate the possibility of breaking through this ceiling, we attempted genomic prediction of phenotypes from genotypes and higher correlations were achieved. We tested distinctness decisions made using either morphological or genotypic distances and found poor correspondence between each method.

Sparse probability density function estimation using the minimum integrated square error

We develop a new sparse kernel density estimator using a forward constrained regression framework, within which the nonnegative and summing-to-unity constraints of the mixing weights can easily be satisfied. Our main contribution is to derive a recursive algorithm to select significant kernels one at time based on the minimum integrated square error (MISE) criterion for both the selection of kernels and the estimation of mixing weights. The proposed approach is simple to implement and the associated computational cost is very low. Specifically, the complexity of our algorithm is in the order of the number of training data N, which is much lower than the order of N2 offered by the best existing sparse kernel density estimators. Numerical examples are employed to demonstrate that the proposed approach is effective in constructing sparse kernel density estimators with comparable accuracy to those of the classical Parzen window estimate and other existing sparse kernel density estimators.

Optimising the analysis of transcript data using high density oligonucleotide arrays and genomic DNA-based probe selection

Background: Affymetrix GeneChip arrays are widely used for transcriptomic studies in a diverse range of species. Each gene is represented on a GeneChip array by a probe- set, consisting of up to 16 probe-pairs. Signal intensities across probe- pairs within a probe-set vary in part due to different physical hybridisation characteristics of individual probes with their target labelled transcripts. We have previously developed a technique to study the transcriptomes of heterologous species based on hybridising genomic DNA (gDNA) to a GeneChip array designed for a different species, and subsequently using only those probes with good homology. Results: Here we have investigated the effects of hybridising homologous species gDNA to study the transcriptomes of species for which the arrays have been designed. Genomic DNA from Arabidopsis thaliana and rice (Oryza sativa) were hybridised to the Affymetrix Arabidopsis ATH1 and Rice Genome GeneChip arrays respectively. Probe selection based on gDNA hybridisation intensity increased the number of genes identified as significantly differentially expressed in two published studies of Arabidopsis development, and optimised the analysis of technical replicates obtained from pooled samples of RNA from rice. Conclusion: This mixed physical and bioinformatics approach can be used to optimise estimates of gene expression when using GeneChip arrays.

Using genomic DNA-based probe-selection to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species

High-density oligonucleotide (oligo) arrays are a powerful tool for transcript profiling. Arrays based on GeneChip® technology are amongst the most widely used, although GeneChip® arrays are currently available for only a small number of plant and animal species. Thus, we have developed a method to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species and tested the method by analysing the transcriptome of Brassica oleracea L., a species for which no GeneChip® array is available, using a GeneChip® array designed for Arabidopsis thaliana (L.) Heynh. Genomic DNA from B. oleracea was labelled and hybridised to the ATH1-121501 GeneChip® array. Arabidopsis thaliana probe-pairs that hybridised to the B. oleracea genomic DNA on the basis of the perfect-match (PM) probe signal were then selected for subsequent B. oleracea transcriptome analysis using a .cel file parser script to generate probe mask files. The transcriptional response of B. oleracea to a mineral nutrient (phosphorus; P) stress was quantified using probe mask files generated for a wide range of gDNA hybridisation intensity thresholds. An example probe mask file generated with a gDNA hybridisation intensity threshold of 400 removed > 68 % of the available PM probes from the analysis but retained >96 % of available A. thaliana probe-sets. Ninety-nine of these genes were then identified as significantly regulated under P stress in B. oleracea, including the homologues of P stress responsive genes in A. thaliana. Increasing the gDNA hybridisation intensity thresholds up to 500 for probe-selection increased the sensitivity of the GeneChip® array to detect regulation of gene expression in B. oleracea under P stress by up to 13-fold. Our open-source software to create probe mask files is freely available http://affymetrix.arabidopsis.info/xspecies/ webcite and may be used to facilitate transcriptomic analyses of a wide range of plant and animal species in the absence of custom arrays.

A fast algorithm for sparse probability density function construction

A new sparse kernel density estimator is introduced. Our main contribution is to develop a recursive algorithm for the selection of significant kernels one at time using the minimum integrated square error (MISE) criterion for both kernel selection. The proposed approach is simple to implement and the associated computational cost is very low. Numerical examples are employed to demonstrate that the proposed approach is effective in constructing sparse kernel density estimators with competitive accuracy to existing kernel density estimators.

Oxidised low-density lipoproteins induce rapid platelet activation and shape change through tyrosine kinase and Rho kinase-signaling pathways

Oxidized low-density lipoproteins (oxLDL) generated in the hyperlipidemic state may contribute to unregulated platelet activation during thrombosis. Although the ability of oxLDL to activate platelets is established, the underlying signaling mechanisms remain obscure. Weshow that oxLDL stimulate platelet activation through phosphorylation of the regulatory light chains of the contractile protein myosin IIa (MLC). oxLDL, but not native LDL, induced shape change, spreading, and phosphorylation of MLC (serine 19) through a pathway that was ablated under conditions that blocked CD36 ligation or inhibited Src kinases, suggesting a tyrosine kinase–dependent mechanism. Consistent with this, oxLDL induced tyrosine phosphorylation of a number of proteins including Syk and phospholipase C g2. Inhibition of Syk, Ca21 mobilization, and MLC kinase (MLCK) only partially inhibited MLC phosphorylation, suggesting the presence of a second pathway. oxLDL activated RhoA and RhoA kinase (ROCK) to induce inhibitory phosphorylation of MLC phosphatase (MLCP). Moreover, inhibition of Src kinases prevented the activation of RhoA and ROCK, indicating that oxLDL regulates contractile signaling through a tyrosine kinase–dependent pathway that induces MLC phosphorylation through the dual activation of MLCK and inhibition of MLCP. These data reveal new signaling events downstream of CD36 that are critical in promoting platelet aggregation by oxLDL.

Available potential energy density for a multicomponent Boussinesq fluid with arbitrary nonlinear equation of state

In this paper, the concept of available potential energy (APE) density is extended to a multicomponent Boussinesq fluid with a nonlinear equation of state. As shown by previous studies, the APE density is naturally interpreted as the work against buoyancy forces that a parcel needs to perform to move from a notional reference position at which its buoyancy vanishes to its actual position; because buoyancy can be defined relative to an arbitrary reference state, so can APE density. The concept of APE density is therefore best viewed as defining a class of locally defined energy quantities, each tied to a different reference state, rather than as a single energy variable. An important result, for which a new proof is given, is that the volume integrated APE density always exceeds Lorenz’s globally defined APE, except when the reference state coincides with Lorenz’s adiabatically re-arranged reference state of minimum potential energy. A parcel reference position is systematically defined as a level of neutral buoyancy (LNB): depending on the nature of the fluid and on how the reference state is defined, a parcel may have one, none, or multiple LNB within the fluid. Multiple LNB are only possible for a multicomponent fluid whose density depends on pressure. When no LNB exists within the fluid, a parcel reference position is assigned at the minimum or maximum geopotential height. The class of APE densities thus defined admits local and global balance equations, which all exhibit a conversion with kinetic energy, a production term by boundary buoyancy fluxes, and a dissipation term by internal diffusive effects. Different reference states alter the partition between APE production and dissipation, but neither affect the net conversion between kinetic energy and APE, nor the difference between APE production and dissipation. We argue that the possibility of constructing APE-like budgets based on reference states other than Lorenz’s reference state is more important than has been previously assumed, and we illustrate the feasibility of doing so in the context of an idealised and realistic oceanic example, using as reference states one with constant density and another one defined as the horizontal mean density field; in the latter case, the resulting APE density is found to be a reasonable approximation of the APE density constructed from Lorenz’s reference state, while being computationally cheaper.