IMILAST: a community effort to intercompare extratropical cyclone detection and tracking algorithms

The variability of results from different automated methods of detection and tracking of extratropical cyclones is assessed in order to identify uncertainties related to the choice of method. Fifteen international teams applied their own algorithms to the same dataset—the period 1989–2009 of interim European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis (ERAInterim) data. This experiment is part of the community project Intercomparison of Mid Latitude Storm Diagnostics (IMILAST; see www.proclim.ch/imilast/index.html). The spread of results for cyclone frequency, intensity, life cycle, and track location is presented to illustrate the impact of using different methods. Globally, methods agree well for geographical distribution in large oceanic regions, interannual variability of cyclone numbers, geographical patterns of strong trends, and distribution shape for many life cycle characteristics. In contrast, the largest disparities exist for the total numbers of cyclones, the detection of weak cyclones, and distribution in some densely populated regions. Consistency between methods is better for strong cyclones than for shallow ones. Two case studies of relatively large, intense cyclones reveal that the identification of the most intense part of the life cycle of these events is robust between methods, but considerable differences exist during the development and the dissolution phases.

Quantifying the climatological relationship between extratropical cyclone intensity and atmospheric precursors

We introduce a novel technique in which linear regression analysis is applied to clusters of tracked cyclones to statistically assess the factors controlling cyclone development. We illustrate this technique by evaluating the differences between cyclones forming in the west and east North Atlantic (herein termed west and east Atlantic cyclones). Enhanced cyclone intensity 2 days after genesis is found to be associated with deeper upper-level troughs upstream of the cyclone center at the genesis time in both west and east Atlantic cyclones. However, whilst west Atlantic cyclones are also enhanced by the presence of strong fronts, east Atlantic cyclones are not. Instead, east Atlantic cyclones exhibit an enhancement when diabatically generated midlevel potential vorticity is present (with the enhancement being of approximately equal magnitude to that associated with the potential vorticity in the upper-level trough). This is consistent with the paradigm of latent heat release in the warm conveyor belt region playing an important role in the development of east Atlantic cyclones.

Shoot calcium and magnesium concentrations differ between subtaxa, are highly heritable, and associate with potentially pleiotropic loci in Brassica oleracea

Calcium (Ca) and magnesium (Mg) are the most abundant group II elements in both plants and animals. Genetic variation in shoot Ca and shoot Mg concentration (shoot Ca and Mg) in plants can be exploited to biofortify food crops and thereby increase dietary Ca and Mg intake for humans and livestock. We present a comprehensive analysis of within-species genetic variation for shoot Ca and Mg, demonstrating that shoot mineral concentration differs significantly between subtaxa (varietas). We established a structured diversity foundation set of 376 accessions to capture a high proportion of species-wide allelic diversity within domesticated Brassica oleracea, including representation of wild relatives (C genome, 1n = 9) from natural populations. These accessions and 74 modern F-1 hybrid cultivars were grown in glasshouse and field environments. Shoot Ca and Mg varied 2- and 2.3-fold, respectively, and was typically not inversely correlated with shoot biomass, within most subtaxa. The closely related capitata (cabbage) and sabauda (Savoy cabbage) subtaxa consistently had the highest mean shoot Ca and Mg. Shoot Ca and Mg in glasshouse-grown plants was highly correlated with data from the field. To understand and dissect the genetic basis of variation in shoot Ca and Mg, we studied homozygous lines from a segregating B. oleracea mapping population. Shoot Ca and Mg was highly heritable (up to 40). Quantitative trait loci (QTL) for shoot Ca and Mg were detected on chromosomes C2, C6, C7, C8, and, in particular, C9, where QTL accounted for 14 to 55 of the total genetic variance. The presence of QTL on C9 was substantiated by scoring recurrent backcross substitution lines, derived from the same parents. This also greatly increased the map resolution, with strong evidence that a 4-cM region on C9 influences shoot Ca. This region corresponds to a 0.41-Mb region on Arabidopsis (Arabidopsis thaliana) chromosome 5 that includes 106 genes. There is also evidence that pleiotropic loci on C8 and C9 affect shoot Ca and Mg. Map-based cloning of these loci will reveal how shoot-level phenotypes relate to Ca 21 and Mg 21 uptake and homeostasis at the molecular level.

Shoot zinc (Zn) concentration varies widely within Brassica oleracea L. and is affected by soil Zn and phosphorus (P) levels

The low availability of zinc (Zn) in soils and crops affects dietary Zn intake worldwide. This study sought to determine if the natural genetic variation in shoot Zn concentrations (Zn(shoot)) is sufficient to pursue a crop improvement breeding strategy in a leafy vegetable crop. The gene-pool of Brassica oleracea L. was sampled using a large (n = 376) diversity foundation set (DFS), representing almost all species-wide common allelic variation, and 74 commercial varieties (mostly F(1)). The DFS genotypes were grown at low and high soil phosphorus (P) levels under glasshouse and field conditions, and also in a Zn-deficient soil, with or without Zn-fertilisation, in a glasshouse. Despite the large variation in Zn(shoot) among genotypes, environment had a profound effect on Zn(shoot) The heritability of Zn(shoot) was significant, but relatively low, among 90 doubled-haploid (DH) lines from a mapping population. While several quantitative trait loci (QTL) associated with Zn(shoot) occurred on chromosomes C2, C3, C5, C7, and C9, these were generally weak and conditional upon growth conditions. Breeding for Zn(shoot) in B. oleracea is therefore likely to be challenging. Shoot P concentrations increased substantially in all genotypes under low soil Zn conditions. Conversely, only some genotypes had increased Zn(shoot) at low soil P levels. Sufficient natural genetic variation may therefore exist to study some of the interactions between Zn and P nutrition.

Error, reproducibility and sensitivity: a pipeline for data processing of Agilent oligonucleotide expression arrays

Background: Expression microarrays are increasingly used to obtain large scale transcriptomic information on a wide range of biological samples. Nevertheless, there is still much debate on the best ways to process data, to design experiments and analyse the output. Furthermore, many of the more sophisticated mathematical approaches to data analysis in the literature remain inaccessible to much of the biological research community. In this study we examine ways of extracting and analysing a large data set obtained using the Agilent long oligonucleotide transcriptomics platform, applied to a set of human macrophage and dendritic cell samples. Results: We describe and validate a series of data extraction, transformation and normalisation steps which are implemented via a new R function. Analysis of replicate normalised reference data demonstrate that intrarray variability is small (only around 2 of the mean log signal), while interarray variability from replicate array measurements has a standard deviation (SD) of around 0.5 log(2) units (6 of mean). The common practise of working with ratios of Cy5/Cy3 signal offers little further improvement in terms of reducing error. Comparison to expression data obtained using Arabidopsis samples demonstrates that the large number of genes in each sample showing a low level of transcription reflect the real complexity of the cellular transcriptome. Multidimensional scaling is used to show that the processed data identifies an underlying structure which reflect some of the key biological variables which define the data set. This structure is robust, allowing reliable comparison of samples collected over a number of years and collected by a variety of operators. Conclusions: This study outlines a robust and easily implemented pipeline for extracting, transforming normalising and visualising transcriptomic array data from Agilent expression platform. The analysis is used to obtain quantitative estimates of the SD arising from experimental (non biological) intra- and interarray variability, and for a lower threshold for determining whether an individual gene is expressed. The study provides a reliable basis for further more extensive studies of the systems biology of eukaryotic cells.

Diversity in expression of phosphorus (P) responsive genes in Cucumis melo L.

Background: Phosphorus (P) is a major limiting nutrient for plant growth in many soils. Studies in model species have identified genes involved in plant adaptations to low soil P availability. However, little information is available on the genetic bases of these adaptations in vegetable crops. In this respect, sequence data for melon now makes it possible to identify melon orthologues of candidate P responsive genes, and the expression of these genes can be used to explain the diversity in the root system adaptation to low P availability, recently observed in this species. Methodology and Findings: Transcriptional responses to P starvation were studied in nine diverse melon accessions by comparing the expression of eight candidate genes (Cm-PAP10.1, Cm-PAP10.2, Cm-RNS1, Cm-PPCK1, Cm-transferase, Cm-SQD1, Cm-DGD1 and Cm-SPX2) under P replete and P starved conditions. Differences among melon accessions were observed in response to P starvation, including differences in plant morphology, P uptake, P use efficiency (PUE) and gene expression. All studied genes were up regulated under P starvation conditions. Differences in the expression of genes involved in P mobilization and remobilization (Cm-PAP10.1, Cm-PAP10.2 and Cm-RNS1) under P starvation conditions explained part of the differences in P uptake and PUE among melon accessions. The levels of expression of the other studied genes were diverse among melon accessions, but contributed less to the phenotypical response of the accessions. Conclusions: This is the first time that these genes have been described in the context of P starvation responses in melon. There exists significant diversity in gene expression levels and P use efficiency among melon accessions as well as significant correlations between gene expression levels and phenotypical measurements.

Gene expression changes in phosphorus deficient potato (Solanum tuberosum L.) leaves and the potential for diagnostic gene expression markers

Background: There are compelling economic and environmental reasons to reduce our reliance on inorganic phosphate (Pi) fertilisers. Better management of Pi fertiliser applications is one option to improve the efficiency of Pi fertiliser use, whilst maintaining crop yields. Application rates of Pi fertilisers are traditionally determined from analyses of soil or plant tissues. Alternatively, diagnostic genes with altered expression under Pi limiting conditions that suggest a physiological requirement for Pi fertilisation, could be used to manage Pifertiliser applications, and might be more precise than indirect measurements of soil or tissue samples. Results: We grew potato (Solanum tuberosum L.) plants hydroponically, under glasshouse conditions, to control their nutrient status accurately. Samples of total leaf RNA taken periodically after Pi was removed from the nutrient solution were labelled and hybridised to potato oligonucleotide arrays. A total of 1,659 genes were significantly differentially expressed following Pi withdrawal. These included genes that encode proteins involved in lipid, protein, and carbohydrate metabolism, characteristic of Pi deficient leaves and included potential novel roles for genes encoding patatin like proteins in potatoes. The array data were analysed using a support vector machine algorithm to identify groups of genes that could predict the Pi status of the crop. These groups of diagnostic genes were tested using field grown potatoes that had either been fertilised or unfertilised. A group of 200 genes could correctly predict the Pi status of field grown potatoes. Conclusions: This paper provides a proof-of-concept demonstration for using microarrays and class prediction tools to predict the Pi status of a field grown potato crop. There is potential to develop this technology for other biotic and abiotic stresses in field grown crops. Ultimately, a better understanding of crop stresses may improve our management of the crop, improving the sustainability of agriculture.

Shoot yield drives phosphorus use efficiency in Brassica oleracea and correlates with root architecture traits

The environmental and financial costs of using inorganic phosphate fertilizers to maintain crop yield and quality are high. Breeding crops that acquire and use phosphorus (P) more efficiently could reduce these costs. The variation in shoot P concentration (shoot-P) and various measures of P use efficiency (PUE) were quantified among 355 Brassica oleracea L. accessions, 74 current commercial cultivars, and 90 doubled haploid (DH) mapping lines from a reference genetic mapping population. Accessions were grown at two or more external P concentrations in glasshouse experiments; commercial and DH accessions were also grown in replicated field experiments. Within the substantial species-wide diversity observed for shoot-P and various measures of PUE in B. oleracea, current commercial cultivars have greater PUE than would be expected by chance. This may be a consequence of breeding for increased yield, which is a significant component of most measures of PUE, or early establishment. Root development and architecture correlate with PUE; in particular, lateral root number, length, and growth rate. Significant quantitative trait loci associated with shoot-P and PUE occur on chromosomes C3 and C7. These data provide information to initiate breeding programmes to improve PUE in B. oleracea.

Regulatory hotspots are associated with plant gene expression under varying soil phosphorus supply in Brassica rapa

Gene expression is a quantitative trait that can be mapped genetically in structured populations to identify expression quantitative trait loci (eQTL). Genes and regulatory networks underlying complex traits can subsequently be inferred. Using a recently released genome sequence, we have defined cis- and trans-eQTL and their environmental response to low phosphorus (P) availability within a complex plant genome and found hotspots of trans-eQTL within the genome. Interval mapping, using P supply as a covariate, revealed 18,876 eQTL. trans-eQTL hotspots occurred on chromosomes A06 and A01 within Brassica rapa; these were enriched with P metabolism-related Gene Ontology terms (A06) as well as chloroplast-and photosynthesis-related terms (A01). We have also attributed heritability components to measures of gene expression across environments, allowing the identification of novel gene expression markers and gene expression changes associated with low P availability. Informative gene expression markers were used to map eQTL and P use efficiency-related QTL. Genes responsive to P supply had large environmental and heritable variance components. Regulatory loci and genes associated with P use efficiency identified through eQTL analysis are potential targets for further characterization and may have potential for crop improvement.

Tandem quadruplication of HMA4 in the zinc (Zn) and cadmium (Cd) hyperaccumulator Noccaea caerulescens

Zinc (Zn) and cadmium (Cd) hyperaccumulation may have evolved twice in the Brassicaceae, in Arabidopsis halleri and in the Noccaea genus. Tandem gene duplication and deregulated expression of the Zn transporter, HMA4, has previously been linked to Zn/Cd hyperaccumulation in A. halleri. Here, we tested the hypothesis that tandem duplication and deregulation of HMA4 expression also occurs in Noccaea. A Noccaea caerulescens genomic library was generated, containing 36,864 fosmid pCC1FOS (TM) clones with insert sizes similar to 20-40 kbp, and screened with a PCR-generated HMA4 genomic probe. Gene copy number within the genome was estimated through DNA fingerprinting and pooled fosmid pyrosequencing. Gene copy numbers within individual clones was determined by PCR analyses with novel locus specific primers. Entire fosmids were then sequenced individually and reads equivalent to 20-fold coverage were assembled to generate complete whole contigs. Four tandem HMA4 repeats were identified in a contiguous sequence of 101,480 bp based on sequence overlap identities. These were flanked by regions syntenous with up and downstream regions of AtHMA4 in Arabidopsis thaliana. Promoter-reporter beta-glucuronidase (GUS) fusion analysis of a NcHMA4 in A. thaliana revealed deregulated expression in roots and shoots, analogous to AhHMA4 promoters, but distinct from AtHMA4 expression which localised to the root vascular tissue. This remarkable consistency in tandem duplication and deregulated expression of metal transport genes between N. caerulescens and A. halleri, which last shared a common ancestor > 40 mya, provides intriguing evidence that parallel evolutionary pathways may underlie Zn/Cd hyperaccumulation in Brassicaceae.

A Brassica exon array for whole-transcript gene expression profiling

Affymetrix GeneChip (R) arrays are used widely to study transcriptional changes in response to developmental and environmental stimuli. GeneChip (R) arrays comprise multiple 25-mer oligonucleotide probes per gene and retain certain advantages over direct sequencing. For plants, there are several public GeneChip (R) arrays whose probes are localised primarily in 39 exons. Plant whole-transcript (WT) GeneChip (R) arrays are not yet publicly available, although WT resolution is needed to study complex crop genomes such as Brassica, which are typified by segmental duplications containing paralogous genes and/or allopolyploidy. Available sequence data were sampled from the Brassica A and C genomes, and 142,997 gene models identified. The assembled gene models were then used to establish a comprehensive public WT exon array for transcriptomics studies. The Affymetrix GeneChip (R) Brassica Exon 1.0 ST Array is a 5 mu M feature size array, containing 2.4 million 25-base oligonucleotide probes representing 135,201 gene models, with 15 probes per gene distributed among exons. Discrimination of the gene models was based on an E-value cut-off of 1E(-5), with <= 98 sequence identity. The 135 k Brassica Exon Array was validated by quantifying transcriptome differences between leaf and root tissue from a reference Brassica rapa line (R-o-18), and categorisation by Gene Ontologies (GO) based on gene orthology with Arabidopsis thaliana. Technical validation involved comparison of the exon array with a 60-mer array platform using the same starting RNA samples. The 135 k Brassica Exon Array is a robust platform. All data relating to the array design and probe identities are available in the public domain and are curated within the BrassEnsembl genome viewer at http://www.brassica.info/BrassEnsembl/index.html.

Within-host evolution of Staphylococcus aureus during asymptomatic carriage

Background Staphylococcus aureus is a major cause of healthcare associated mortality, but like many important bacterial pathogens, it is a common constituent of the normal human body flora. Around a third of healthy adults are carriers. Recent evidence suggests that evolution of S. aureus during nasal carriage may be associated with progression to invasive disease. However, a more detailed understanding of within-host evolution under natural conditions is required to appreciate the evolutionary and mechanistic reasons why commensal bacteria such as S. aureus cause disease. Therefore we examined in detail the evolutionary dynamics of normal, asymptomatic carriage. Sequencing a total of 131 genomes across 13 singly colonized hosts using the Illumina platform, we investigated diversity, selection, population dynamics and transmission during the short-term evolution of S. aureus. Principal Findings We characterized the processes by which the raw material for evolution is generated: micro-mutation (point mutation and small insertions/deletions), macro-mutation (large insertions/deletions) and the loss or acquisition of mobile elements (plasmids and bacteriophages). Through an analysis of synonymous, non-synonymous and intergenic mutations we discovered a fitness landscape dominated by purifying selection, with rare examples of adaptive change in genes encoding surface-anchored proteins and an enterotoxin. We found evidence for dramatic, hundred-fold fluctuations in the size of the within-host population over time, which we related to the cycle of colonization and clearance. Using a newly-developed population genetics approach to detect recent transmission among hosts, we revealed evidence for recent transmission between some of our subjects, including a husband and wife both carrying populations of methicillin-resistant S. aureus (MRSA). Significance This investigation begins to paint a picture of the within-host evolution of an important bacterial pathogen during its prevailing natural state, asymptomatic carriage. These results also have wider significance as a benchmark for future systematic studies of evolution during invasive S. aureus disease.

Longitudinal investigation of the role of temperament and stressful life events in childhood anxiety

The current study investigated the longitudinal relationships between BI, life events, and anxiety in a sample of 102 behaviourally inhibited (BI) and 100 uninhibited (BUI) children aged 3 to 4 years. Children’s parents completed questionnaires on BI, stressful life events, and anxiety symptoms, and were administered a diagnostic interview three times in a 5-year period. In line with our hypotheses, negative life events, and negative behaviour- dependent life events (i.e. life events that are related to the children’s own behaviours) in particular, and the impact of negative life events, were predictive of increases in subsequent anxiety symptoms, the likelihood of having an anxiety disorder, and increased number of anxiety diagnoses over the five year follow-up period. Experiencing more positive, behaviour-independent life events decreased the risk of being diagnosed with an anxiety disorder. Furthermore, differences were found in life events between BI and BUI children. That is, BI children experienced fewer positive and specifically positive behaviour-dependent life events, and the impact of these positive life events was also lower in BI children than in BUI children. However, BI did not interact with life events in the prediction of anxiety problems as hypothesized. Therefore, this study seems to indicate that BI and life events act as additive risk factors in the development of anxiety problems.

Informing early intervention: preschool predictors of anxiety disorders in middle childhood

Background: To inform early intervention practice, the present research examines how child anxiety, behavioural inhibition, maternal overinvolvement, maternal negativity, mother-child attachment and maternal anxiety, as assessed at age four, predict anxiety at age nine. Method: 202 children (102 behaviourally inhibited and 100 behaviourally uninhibited) aged 3–4 years were initially recruited and the predictors outlined above were assessed. Diagnostic assessments, using the Anxiety Disorders Interview Schedule, were then conducted five years later. Results: Behavioural inhibition, maternal anxiety, and maternal overinvolvement were significant predictors of clinical anxiety, even after controlling for baseline anxiety (p,.05). No significant effect of negativity or attachment security was found over and above baseline anxiety (p..1). Conclusions: Preschool children who show anxiety, are inhibited, have overinvolved mothers and mothers with anxiety disorders are at increased risk for anxiety in middle childhood. These factors can be used to identify suitable participants for early intervention and can be targeted within intervention programs.

Interpretation bias in preschool children at risk for anxiety: a prospective study

A story-stem paradigm was used to assess interpretation bias in preschool children. Data were available for 131 children. Interpretation bias, behavioural inhibition (BI) and anxiety were assessed when children were aged between 3 years 2 months and 4 years 5 months. Anxiety was subsequently assessed 12-months, 2-years and 5-years later. A significant difference in interpretation bias was found between participants who met criteria for an anxiety diagnosis at baseline, with clinically anxious participants more likely to complete the ambiguous story-stems in a threat-related way. Threat interpretations significantly predicted anxiety symptoms at 12-month follow-up, after controlling for baseline symptoms, but did not predict anxiety symptoms or diagnoses at either 2-year or 5- year follow-up. There was little evidence for a relationship between BI and interpretation bias. Overall, the pattern of results was not consistent with the hypothesis that interpretation bias plays a role in the development of anxiety. Instead, some evidence for a role in the maintenance of anxiety over relatively short periods of time was found. The use of a story-stem methodology to assess interpretation bias in young children is discussed along with the theoretical and clinical implications of the findings.