Numerical analyses of Runge–Kutta implicit–explicit schemes for horizontally explicit, vertically implicit solutions of atmospheric models

With the prospect of exascale computing, computational methods requiring only local data become especially attractive. Consequently, the typical domain decomposition of atmospheric models means horizontally-explicit vertically-implicit (HEVI) time-stepping schemes warrant further attention. In this analysis, Runge-Kutta implicit-explicit schemes from the literature are analysed for their stability and accuracy using a von Neumann stability analysis of two linear systems. Attention is paid to the numerical phase to indicate the behaviour of phase and group velocities. Where the analysis is tractable, analytically derived expressions are considered. For more complicated cases, amplification factors have been numerically generated and the associated amplitudes and phase diagnosed. Analysis of a system describing acoustic waves has necessitated attributing the three resultant eigenvalues to the three physical modes of the system. To do so, a series of algorithms has been devised to track the eigenvalues across the frequency space. The result enables analysis of whether the schemes exactly preserve the non-divergent mode; and whether there is evidence of spurious reversal in the direction of group velocities or asymmetry in the damping for the pair of acoustic modes. Frequency ranges that span next-generation high-resolution weather models to coarse-resolution climate models are considered; and a comparison is made of errors accumulated from multiple stability-constrained shorter time-steps from the HEVI scheme with a single integration from a fully implicit scheme over the same time interval. Two schemes, “Trap2(2,3,2)” and “UJ3(1,3,2)”, both already used in atmospheric models, are identified as offering consistently good stability and representation of phase across all the analyses. Furthermore, according to a simple measure of computational cost, “Trap2(2,3,2)” is the least expensive.

Identification and expression analyses of Cytosolic Glutamine Synthetase genes in barley (Hordeum vulgare L.)

Glutamine synthetase (GS) is a key enzyme in nitrogen (N) assimilation, particularly during seed development. Three cytosolic GS isoforms (HvGS1) were identified in barley (Hordeum vulgare L. cv Golden Promise). Quantitation of gene expression, localization and response to N supply revealed that each gene plays a non-redundant role in different tissues and during development. Localization of HvGS1_1 in vascular cells of different tissues, combined with its abundance in the stem and its response to changes in N supply, indicate that it is important in N transport and remobilization. HvGS1_1 is located on chromosome 6H at 72.54 cM, close to the marker HVM074 which is associated with a major quantitative trait locus (QTL) for grain protein content (GPC). HvGS1_1 may be a potential candidate gene to manipulate barley GPC. HvGS1_2 mRNA was localized to the leaf mesophyll cells, in the cortex and pericycle of roots, and was the dominant HvGS1 isoform in these tissues. HvGS1_2 expression increased in leaves with an increasing supply of N, suggesting its role in the primary assimilation of N. HvGS1_3 was specifically and predominantly localized in the grain, being highly expressed throughout grain development. HvGS1_3 expression increased specifically in the roots of plants grown on high NH+4, suggesting that it has a primary role in grain N assimilation and also in the protection against ammonium toxicity in roots. The expression of HvGS1 genes is directly correlated with protein and enzymatic activity, indicating that transcriptional regulation is of prime importance in the control of GS activity in barley.

High-resolution melt and morphological analyses of mealybugs (Hemiptera: Pseudococcidae) from cacao: tools for the control of Cacao Swollen Shoot Virus spread

BACKGROUND: Mealybugs (Hemiptera: Coccoidea: Pseudococcidae) are key vectors of badnaviruses, including Cacao Swollen Shoot Virus (CSSV) the most damaging virus affecting cacao (Theobroma cacao L.). The effectiveness of mealybugs as virus vectors is species dependent and it is therefore vital that CSSV resistance breeding programmes in cacao incorporate accurate mealybug identification. In this work the efficacy of a CO1-based DNA barcoding approach to species identification was evaluated by screening a range of mealybugs collected from cacao in seven countries. RESULTS: Morphologically similar adult females were characterised by scanning electron microscopy and then, following DNA extraction, were screened with CO1 barcoding markers. A high degree of CO1 sequence homology was observed for all 11 individual haplotypes including those accessions from distinct geographical regions. This has allowed for the design of a High Resolution Melt (HRM) assay capable of rapid identification of the commonly encountered mealybug pests of cacao. CONCLUSIONS: HRM Analysis (HRMA) readily differentiated between mealybug pests of cacao that can not necessarily be identified by conventional morphological analysis. This new approach, therefore, has potential to facilitate breeding for resistance to CSSV and other mealybug transmitted diseases.

GRID and docking analyses reveal a molecular basis for flavonoid inhibition of src-family kinase activity

Flavonoids reduce cardiovascular disease risk through anti-inflammatory, anti-coagulant and anti-platelet actions. One key flavonoid inhibitory mechanism is blocking kinase activity that drives these processes. Flavonoids attenuate activities of kinases including phosphoinositide-3-kinase (PI3K), Fyn, Lyn, Src, Syk, PKC, PIM1/2, ERK, JNK, and PKA. X-ray crystallographic analyses of kinase-flavonoid complexes show that flavonoid ring systems and their hydroxyl substitutions are important structural features for their binding to kinases. A clearer understanding of structural interactions of flavonoids with kinases is necessary to allow construction of more potent and selective counterparts. We examined flavonoid (quercetin, apigenin and catechin) interactions with Src-family kinases (Lyn, Fyn and Hck) applying the Sybyl docking algorithm and GRID. A homology model (Lyn) was used in our analyses to demonstrate that high quality predicted kinase structures are suitable for flavonoid computational studies. Our docking results revealed potential hydrogen bond contacts between flavonoid hydroxyls and kinase catalytic site residues. Identification of plausible contacts indicated that quercetin formed the most energetically stable interactions, apigenin lacked hydroxyl groups necessary for important contacts, and the non-planar structure of catechin could not support predicted hydrogen bonding patterns. GRID analysis using a hydroxyl functional group supported docking results. Based on these findings, we predicted that quercetin would inhibit activities of Src-family kinases with greater potency than apigenin and catechin. We validated this prediction using in vitro kinase assays. We conclude that our study can be used as a basis to construct virtual flavonoid interaction libraries to guide drug discovery using these compounds as molecular templates.

Time-dependent gel to gel transformation of a peptide based supramolecular gelator

A dipeptide with a long fatty acid chain at its N-terminus gives hydrogels in phosphate buffer in the pH range 7.0–8.5. The hydrogel with a gelator concentration of 0.45% (w/v) at pH 7.46 (physiological pH) provides a very good platform to study dynamic changes within a supramolecular framework as it exhibits remarkable change in its appearance with time. Interestingly, the first formed transparent hydrogel gradually transforms into a turbid gel within 2 days. These two forms of the hydrogel have been thoroughly investigated by using small angle X-ray scattering (SAXS), powder X-ray diffraction (PXRD), field emission scanning electron microscopic (FE-SEM) and high-resolution transmission electron microscopic (HR-TEM) imaging, FT-IR and rheometric analyses. The SAXS and low angle PXRD studies substantiate different packing arrangements for the gelator molecules for these two different gel states (the freshly prepared and the aged hydrogel). Moreover, rheological studies of these two gels reveal that the aged gel is stiffer than the freshly prepared gel.

Comparative in situ analyses of cell wall matrix polysaccharide dynamics in developing rice and wheat grain

Cell wall polysaccharides of wheat and rice endosperm are an important source of dietary fibre. Monoclonal antibodies specific to cell wall polysaccharides were used to determine polysaccharide dynamics during the development of both wheat and rice grain. Wheat and rice grain present near synchronous developmental processes and significantly different endosperm cell wall compositions, allowing the localisation of these polysaccharides to be related to developmental changes. Arabinoxylan (AX) and mixed-linkage glucan (MLG) have analogous cellular locations in both species, with deposition of AX and MLG coinciding with the start of grain filling. A glucuronoxylan (GUX) epitope was detected in rice, but not wheat endosperm cell walls. Callose has been reported to be associated with the formation of cell wall outgrowths during endosperm cellularisation and xyloglucan is here shown to be a component of these anticlinal extensions, occurring transiently in both species. Pectic homogalacturonan (HG) was abundant in cell walls of maternal tissues of wheat and rice grain, but only detected in endosperm cell walls of rice in an unesterified HG form. A rhamnogalacturonan-I (RG-I) backbone epitope was observed to be temporally regulated in both species, detected in endosperm cell walls from 12 DAA in rice and 20 DAA in wheat grain. Detection of the LM5 galactan epitope showed a clear distinction between wheat and rice, being detected at the earliest stages of development in rice endosperm cell walls, but not detected in wheat endosperm cell walls, only in maternal tissues. In contrast, the LM6 arabinan epitope was detected in both species around 8 DAA and was transient in wheat grain, but persisted in rice until maturity.

Higher PUFA and omega-3 PUFA, CLA, α-tocopherol and iron, but lower iodine and selenium concentrations in organic milk: a systematic literature review and meta- and redundancy analyses

Demand for organic milk is partially driven by consumer perceptions that it is more nutritious. However, there is still considerable uncertainty over whether the use of organic production standards affects milk quality. Here we report results of meta-analyses based on 170 published studies comparing the nutrient content of organic and conventional bovine milk. There were no significant differences in total SFA and MUFA concentrations between organic and conventional milk. However, concentrations of total PUFA and n-3 PUFA were significantly higher in organic milk, by an estimated 7 (95 % CI −1, 15) % and 56 (95 % CI 38, 74) %, respectively. Concentrations of α-linolenic acid (ALA), very long-chain n-3 fatty acids (EPA+DPA+DHA) and conjugated linoleic acid were also significantly higher in organic milk, by an 69 (95 % CI 53, 84) %, 57 (95 % CI 27, 87) % and 41 (95 % CI 14, 68) %, respectively. As there were no significant differences in total n-6 PUFA and linoleic acid (LA) concentrations, the n-6:n-3 and LA:ALA ratios were lower in organic milk, by an estimated 71 (95 % CI −122, −20) % and 93 (95 % CI −116, −70) %. It is concluded that organic bovine milk has a more desirable fatty acid composition than conventional milk. Meta-analyses also showed that organic milk has significantly higher α-tocopherol and Fe, but lower I and Se concentrations. Redundancy analysis of data from a large cross-European milk quality survey indicates that the higher grazing/conserved forage intakes in organic systems were the main reason for milk composition differences.

Biases in comparative analyses of extinction risk: mind the gap

1. Comparative analyses are used to address the key question of what makes a species more prone to extinction by exploring the links between vulnerability and intrinsic species’ traits and/or extrinsic factors. This approach requires comprehensive species data but information is rarely available for all species of interest. As a result comparative analyses often rely on subsets of relatively few species that are assumed to be representative samples of the overall studied group. 2. Our study challenges this assumption and quantifies the taxonomic, spatial, and data type biases associated with the quantity of data available for 5415 mammalian species using the freely available life-history database PanTHERIA. 3. Moreover, we explore how existing biases influence results of comparative analyses of extinction risk by using subsets of data that attempt to correct for detected biases. In particular, we focus on links between four species’ traits commonly linked to vulnerability (distribution range area, adult body mass, population density and gestation length) and conduct univariate and multivariate analyses to understand how biases affect model predictions. 4. Our results show important biases in data availability with c.22% of mammals completely lacking data. Missing data, which appear to be not missing at random, occur frequently in all traits (14–99% of cases missing). Data availability is explained by intrinsic traits, with larger mammals occupying bigger range areas being the best studied. Importantly, we find that existing biases affect the results of comparative analyses by overestimating the risk of extinction and changing which traits are identified as important predictors. 5. Our results raise concerns over our ability to draw general conclusions regarding what makes a species more prone to extinction. Missing data represent a prevalent problem in comparative analyses, and unfortunately, because data are not missing at random, conventional approaches to fill data gaps, are not valid or present important challenges. These results show the importance of making appropriate inferences from comparative analyses by focusing on the subset of species for which data are available. Ultimately, addressing the data bias problem requires greater investment in data collection and dissemination, as well as the development of methodological approaches to effectively correct existing biases.

Interim analyses and sequential designs in phase III studies

Recruitment of patients to a clinical trial usually occurs over a period of time, resulting in the steady accumulation of data throughout the trial's duration. Yet, according to traditional statistical methods, the sample size of the trial should be determined in advance, and data collected on all subjects before analysis proceeds. For ethical and economic reasons, the technique of sequential testing has been developed to enable the examination of data at a series of interim analyses. The aim is to stop recruitment to the study as soon as there is sufficient evidence to reach a firm conclusion. In this paper we present the advantages and disadvantages of conducting interim analyses in phase III clinical trials, together with the key steps to enable the successful implementation of sequential methods in this setting. Examples are given of completed trials, which have been carried out sequentially, and references to relevant literature and software are provided.

Steric sea level variability (1993-2010) in an ensemble of ocean reanalyses and objective analyses

Quantifying the effect of the seawater density changes on sea level variability is of crucial importance for climate change studies, as the sea level cumulative rise can be regarded as both an important climate change indicator and a possible danger for human activities in coastal areas. In this work, as part of the Ocean Reanalysis Intercomparison Project, the global and regional steric sea level changes are estimated and compared from an ensemble of 16 ocean reanalyses and 4 objective analyses. These estimates are initially compared with a satellite-derived (altimetry minus gravimetry) dataset for a short period (2003–2010). The ensemble mean exhibits a significant high correlation at both global and regional scale, and the ensemble of ocean reanalyses outperforms that of objective analyses, in particular in the Southern Ocean. The reanalysis ensemble mean thus represents a valuable tool for further analyses, although large uncertainties remain for the inter-annual trends. Within the extended intercomparison period that spans the altimetry era (1993–2010), we find that the ensemble of reanalyses and objective analyses are in good agreement, and both detect a trend of the global steric sea level of 1.0 and 1.1 ± 0.05 mm/year, respectively. However, the spread among the products of the halosteric component trend exceeds the mean trend itself, questioning the reliability of its estimate. This is related to the scarcity of salinity observations before the Argo era. Furthermore, the impact of deep ocean layers is non-negligible on the steric sea level variability (22 and 12 % for the layers below 700 and 1500 m of depth, respectively), although the small deep ocean trends are not significant with respect to the products spread.

Genetic diversity at the Dhn3 locus in Turkish Hordeum spontaneum populations with comparative structural analyses

We analysed Hordeum spontaneum accessions from 21 different locations to understand the genetic diversity of HsDhn3 alleles and effects of single base mutations on the intrinsically disordered structure of the resulting polypeptide (HsDHN3). HsDHN3 was found to be YSK2-type with a low-frequency 6-aa deletion in the beginning of Exon 1. There is relatively high diversity in the intron region of HsDhn3 compared to the two exon regions. We have found subtle differences in K segments led to changes in amino acids chemical properties. Predictions for protein interaction profiles suggest the presence of a protein-binding site in HsDHN3 that coincides with the K1 segment. Comparison of DHN3 to closely related cereals showed that all of them contain a nuclear localization signal sequence flanking to the K1 segment and a novel conserved region located between the S and K1 segments [E(D/T)DGMGGR]. We found that H. vulgare, H. spontaneum, and Triticum urartu DHN3s have a greater number of phosphorylation sites for protein kinase C than other cereal species, which may be related to stress adaptation. Our results show that the nature and extent of mutations in the conserved segments of K1 and K2 are likely to be key factors in protection of cells.

abctools: an R package for tuning approximate Bayesian computation analyses

Approximate Bayesian computation (ABC) is a popular family of algorithms which perform approximate parameter inference when numerical evaluation of the likelihood function is not possible but data can be simulated from the model. They return a sample of parameter values which produce simulations close to the observed dataset. A standard approach is to reduce the simulated and observed datasets to vectors of summary statistics and accept when the difference between these is below a specified threshold. ABC can also be adapted to perform model choice. In this article, we present a new software package for R, abctools which provides methods for tuning ABC algorithms. This includes recent dimension reduction algorithms to tune the choice of summary statistics, and coverage methods to tune the choice of threshold. We provide several illustrations of these routines on applications taken from the ABC literature.