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

Evidence of neutral transcriptome evolution in plants

The transcriptome of an organism is its set of gene transcripts (mRNAs) at a defined spatial and temporal locus. Because gene expression is affected markedly by environmental and developmental perturbations, it is widely assumed that transcriptome divergence among taxa represents adaptive phenotypic selection. This assumption has been challenged by neutral theories which propose that stochastic processes drive transcriptome evolution. To test for evidence of neutral transcriptome evolution in plants, we quantified 18 494 gene transcripts in nonsenescent leaves of 14 taxa of Brassicaceae using robust cross-species transcriptomics which includes a two-step physical and in silicobased normalization procedure based on DNA similarity among taxa. Transcriptome divergence correlates positively with evolutionary distance between taxa and with variation in gene expression among samples. Results are similar for pseudogenes and chloroplast genes evolving at different rates. Remarkably, variation in transcript abundance among root-cell samples correlates positively with transcriptome divergence among root tissues and among taxa. Because neutral processes affect transcriptome evolution in plants, many differences in gene expression among or within taxa may be nonfunctional, reflecting ancestral plasticity and founder effects. Appropriate null models are required when comparing transcriptomes in space and time

Using art to illuminate social workers' stress

Summary: This article aims to capture the self-defined holistic interaction between stressors, stress reactions and coping for social workers’ stress. We utilised an arts-based intervention in the form of a single drawing with 80 social workers, who were actively guided to explore their own stressors, stress reactions and coping. Findings: Our findings suggest that whilst social workers define their stressors as being related to a lack of social work professional and managerial support (a macro problem), they experience this stress as lack of personal efficacy and self-worth and expect to cope by drawing on their inner strengths rather than by challenging the ‘system’ (micro solutions). Thus, they tend not to utilise available systemic ideas and theories in social work to address their own problems. Application:These findings offer a way of exploring stress and coping as an interactive whole that helps to understand both systemic stressors and the gaps in social workers’ coping methods. It also explores the relationship between stressors, stress reactions and coping through the personal drawings and narratives of participating social workers. The findings are relevant for supervisors, and managers helping social workers to manage stress, and offer an example of how visual methods might be used as a pedagogic tool in social work education and practice

Intercomparison of modal and sectional aerosol microphysics representations within the same 3-D global chemical transport model

In the most advanced aerosol-climate models it is common to represent the aerosol particle size distribution in terms of several log-normal modes. This approach, motivated by computational efficiency, makes assumptions about the shape of the particle distribution that may not always capture the properties of global aerosol. Here, a global modal aerosol microphysics module (GLOMAP-mode) is evaluated and improved by comparing against a sectional version (GLOMAP-bin) and observations in the same 3-D global offline chemistry transport model. With both schemes, the model captures the main features of the global particle size distribution, with sub-micron aerosol approximately unimodal in continental regions and bi-modal in marine regions. Initial bin-mode comparisons showed that the current values for two size distribution parameter settings in the modal scheme (mode widths and inter-modal separation sizes) resulted in clear biases compared to the sectional scheme. By adjusting these parameters in the modal scheme, much better agreement is achieved against the bin scheme and observations. Annual mean surface-level mass of sulphate, sea-salt, black carbon (BC) and organic carbon (OC) are within 25% in the two schemes in nearly all regions. Surface level concentrations of condensation nuclei (CN), cloud condensation nuclei (CCN), surface area density and condensation sink also compare within 25% in most regions. However, marine CCN concentrations between 30° N and 30° S are systematically 25–60% higher in the modal model, which we attribute to differences in size-resolved particle growth or cloud-processing. Larger differences also exist in regions or seasons dominated by biomass burning and in free-troposphere and high-latitude regions. Indeed, in the free-troposphere, GLOMAP-mode BC is a factor 2–4 higher than GLOMAP-bin, likely due to differences in size-resolved scavenging. Nevertheless, in most parts of the atmosphere, we conclude that bin-mode differences are much less than model-observation differences, although some processes are missing in these runs which may pose a bigger challenge to modal schemes (e.g., boundary layer nucleation and ultra-fine sea-spray). The findings here underline the need for a spectrum of complexity in global models, with size-resolved aerosol properties predicted by modal schemes needing to be continually benchmarked and improved against freely evolving sectional schemes and observations

Aerosol microphysics simulations of the Mt.~Pinatubo eruption with the UM-UKCA composition-climate model

We use a stratosphere–troposphere composition–climate model with interactive sulfur chemistry and aerosol microphysics, to investigate the effect of the 1991 Mount Pinatubo eruption on stratospheric aerosol properties. Satellite measurements indicate that shortly after the eruption, between 14 and 23 Tg of SO2 (7 to 11.5 Tg of sulfur) was present in the tropical stratosphere. Best estimates of the peak global stratospheric aerosol burden are in the range 19 to 26 Tg, or 3.7 to 6.7 Tg of sulfur assuming a composition of between 59 and 77 % H2SO4. In light of this large uncertainty range, we performed two main simulations with 10 and 20 Tg of SO2 injected into the tropical lower stratosphere. Simulated stratospheric aerosol properties through the 1991 to 1995 period are compared against a range of available satellite and in situ measurements. Stratospheric aerosol optical depth (sAOD) and effective radius from both simulations show good qualitative agreement with the observations, with the timing of peak sAOD and decay timescale matching well with the observations in the tropics and mid-latitudes. However, injecting 20 Tg gives a factor of 2 too high stratospheric aerosol mass burden compared to the satellite data, with consequent strong high biases in simulated sAOD and surface area density, with the 10 Tg injection in much better agreement. Our model cannot explain the large fraction of the injected sulfur that the satellite-derived SO2 and aerosol burdens indicate was removed within the first few months after the eruption. We suggest that either there is an additional alternative loss pathway for the SO2 not included in our model (e.g. via accommodation into ash or ice in the volcanic cloud) or that a larger proportion of the injected sulfur was removed via cross-tropopause transport than in our simulations. We also critically evaluate the simulated evolution of the particle size distribution, comparing in detail to balloon-borne optical particle counter (OPC) measurements from Laramie, Wyoming, USA (41° N). Overall, the model captures remarkably well the complex variations in particle concentration profiles across the different OPC size channels. However, for the 19 to 27 km injection height-range used here, both runs have a modest high bias in the lowermost stratosphere for the finest particles (radii less than 250 nm), and the decay timescale is longer in the model for these particles, with a much later return to background conditions. Also, whereas the 10 Tg run compared best to the satellite measurements, a significant low bias is apparent in the coarser size channels in the volcanically perturbed lower stratosphere. Overall, our results suggest that, with appropriate calibration, aerosol microphysics models are capable of capturing the observed variation in particle size distribution in the stratosphere across both volcanically perturbed and quiescent conditions. Furthermore, additional sensitivity simulations suggest that predictions with the models are robust to uncertainties in sub-grid particle formation and nucleation rates in the stratosphere

Epidemiological associations between brachycephaly and upper respiratory tract disorders in dogs attending veterinary practices in England

Background: Brachycephalic dog breeds are increasingly common. Canine brachycephaly has been associated with upper respiratory tract (URT) disorders but reliable prevalence data remain lacking. Using primary-care veterinary clinical data, this study aimed to report the prevalence and breed-type risk factors for URT disorders in dogs. Results: The sampling frame included 170,812 dogs attending 96 primary-care veterinary clinics participating within the VetCompass Programme. Two hundred dogs were randomly selected from each of three extreme brachycephalic breed types (Bulldog, French Bulldog and Pug) and three common small-to medium sized breed types (moderate brachycephalic: Yorkshire Terrier and non-brachycephalic: Border Terrier and West Highland White Terrier). Information on all URT disorders recorded was extracted from individual patient records. Disorder prevalence was compared between groups using the chi-squared test or Fisher’s test, as appropriate. Risk factor analysis used multivariable logistic regression modelling. During the study, 83 (6.9 %) study dogs died. Extreme brachycephalic dogs (median longevity: 8.6 years, IQR: 2.4-10.8) were significantly younger at death than the moderate and non-brachycephalic group of dogs (median 12.7 years, IQR 11.1-15.0) (P \u3c 0.001). A higher proportion of deaths in extreme brachycephalic breed types were associated with URT disorders (4/24 deaths, 16.7 %) compared with the moderate and non-brachycephalic group (0/59 deaths, 0.0 %) (P = 0.001). The prevalence of having at least one URT disorder in the extreme brachycephalic group was higher (22.0 %, 95 % confidence interval (CI): 18.0-26.0) than in the moderate and non-brachycephalic group (9.7 %, 95 % CI: 7.1-12.3, P \u3c 0.001). The prevalence of URT disorders varied significantly by breed type: Bulldogs 19.5 %, French Bulldogs 20.0 %, Pugs 26.5 %, Border Terriers 9.0 %, West Highland White Terriers 7.0 % and Yorkshire Terriers 13.0 % (P \u3c 0.001). After accounting for the effects of age, bodyweight, sex, neutering and insurance, extreme brachycephalic dogs had 3.5 times (95 % CI: 2.4-5.0, P \u3c 0.001) the odds of at least one URT disorder compared with the moderate and non-brachycephalic group. Conclusions: In summary, this study reports that URT disorders are commonly diagnosed in Bulldog, French Bulldog, Pug, Border Terrier, WHWT and Yorkshire Terrier dogs attending primary-care veterinary practices in England. The three extreme brachycephalic breed types (Bulldog, French Bulldog and Pug) were relatively short-lived and predisposed to URT disorders compared with three other small-to-medium size breed types that are commonly owned (moderate brachycephalic Yorkshire Terrier and non-brachycephalic: Border Terrier and WHWT). Conclusions: In summary, this study reports that URT disorders are commonly diagnosed in Bulldog, French Bulldog, Pug, Border Terrier, WHWT and Yorkshire Terrier dogs attending primary-care veterinary practices in England. The three extreme brachycephalic breed types (Bulldog, French Bulldog and Pug) were relatively short-lived and predisposed to URT disorders compared with three other small-to-medium size breed types that are commonly owned (moderate brachycephalic Yorkshire Terrier and non-brachycephalic: Border Terrier and WHWT)

Organism-sediment interactions govern post-hypoxia recovery of ecosystem functioning

Hypoxia represents one of the major causes of biodiversity and ecosystem functioning loss for coastal waters. Since eutrophication-induced hypoxic events are becoming increasingly frequent and intense, understanding the response of ecosystems to hypoxia is of primary importance to understand and predict the stability of ecosystem functioning. Such ecological stability may greatly depend on the recovery patterns of communities and the return time of the system properties associated to these patterns. Here, we have examined how the reassembly of a benthic community contributed to the recovery of ecosystem functioning following experimentally-induced hypoxia in a tidal flat. We demonstrate that organism-sediment interactions that depend on organism size and relate to mobility traits and sediment reworking capacities are generally more important than recovering species richness to set the return time of the measured sediment processes and properties. Specifically, increasing macrofauna bioturbation potential during community reassembly significantly contributed to the recovery of sediment processes and properties such as denitrification, bedload sediment transport, primary production and deep pore water ammonium concentration. Such bioturbation potential was due to the replacement of the small-sized organisms that recolonised at early stages by large-sized bioturbating organisms, which had a disproportionately stronger influence on sediment. This study suggests that the complete recovery of organism-sediment interactions is a necessary condition for ecosystem functioning recovery, and that such process requires long periods after disturbance due to the slow growth of juveniles into adult stages involved in these interactions. Consequently, repeated episodes of disturbance at intervals smaller than the time needed for the system to fully recover organism-sediment interactions may greatly impair the resilience of ecosystem functioning.

Modelling chemistry in the nocturnal boundary layer above tropical rainforest and a generalised effective nocturnal ozone deposition velocity for sub-ppbv NOx conditions

Measurements of atmospheric composition have been made over a remote rainforest landscape. A box model has previously been demonstrated to model the observed daytime chemistry well. However the box model is unable to explain the nocturnal measurements of relatively high [NO] and [O3], but relatively low observed [NO2]. It is shown that a one-dimensional (1-D) column model with simple O3 -NOx chemistry and a simple representation of vertical transport is able to explain the observed nocturnal concentrations and predict the likely vertical profiles of these species in the nocturnal boundary layer (NBL). Concentrations of tracers carried over from the end of the night can affect the atmospheric chemistry of the following day. To ascertain the anomaly introduced by using the box model to represent the NBL, vertically-averaged NBL concentrations at the end of the night are compared between the 1-D model and the box model. It is found that, under low to medium [NOx] conditions (NOx <1 ppbv), a simple parametrisation can be used to modify the box model deposition velocity of ozone, in order to achieve good agreement between the box and 1-D models for these end-of-night concentrations of NOx and O3. This parametrisation would could also be used in global climate-chemistry models with limited vertical resolution near the surface. Box-model results for the following day differ significantly if this effective nocturnal deposition velocity for ozone is implemented; for instance, there is a 9% increase in the following day’s peak ozone concentration. However under medium to high [NOx] conditions (NOx > 1 ppbv), the effect on the chemistry due to the vertical distribution of the species means no box model can adequately represent chemistry in the NBL without modifying reaction rate constants

Obinutuzumab as consolidation after chemo-immunotherapy: Results of the UK National Cancer Research Institute phase II/III GALACTIC trial

The GA101 (obinutuzumab) monocLonal Antibody as Consolidation Therapy In chronic lymphocytic leukaemia (CLL) (GALACTIC) was a seamless phase II/III trial designed to test whether consolidation with obinutuzumab is safe and eradicates minimal residual disease (MRD) and, subsequently, whether this leads to prolonged progression-free survival (PFS) in patients with CLL who have recently responded to chemo-immunotherapy. Patients with a response 3–24 months after chemotherapy were assessed for MRD. MRD-positive patients were randomised to receive consolidation therapy with obinutuzumab or no consolidation. The trial closed after the phase II part due to slow recruitment. In all, 48 patients enrolled of whom 19 were MRD negative and were monitored. Of the 29 MRD-positive patients, 14 were randomised to receive consolidation and 15 to no consolidation. At 6 months after randomisation, 10 and 13 consolidated patients achieved MRD negativity by flow cytometry (sensitivity 10−4) in bone marrow and peripheral blood respectively. PFS was significantly better in consolidated patients compared to non-consolidated patients (p = 0.001). No difference was observed in PFS, overall survival or duration of MRD negativity when comparing the 10 MRD-negative patients after consolidation with the 19 MRD-negative patients in the monitoring group. Common adverse events in the consolidation arm were thrombocytopenia, infection, and cough. Only 1% of events were infusion-related reactions. This observation provides further evidence that consolidation to achieve MRD negativity improves outcomes in CLL and that obinutuzumab is well tolerated in patients with low levels of disease