In planta conditions induce genomic changes in Pseudomonas syringae pv. phaseolicola

The co-evolution of bacterial plant pathogens and their hosts is a complex and dynamic process. Plant resistance can impose stress on invading pathogens that can lead to, and select for, beneficial changes in the bacterial genome. The Pseudomonas syringae pv. phaseolicola (Pph) genomic island PPHGI-1 carries an effector gene, avrPphB (hopAR1), which triggers the hypersensitive reaction in bean plants carrying the R3 resistance gene. Interaction between avrPphB and R3 generates an antimicrobial environment within the plant, resulting in the excision of PPHGI-1 and its loss from the genome. The loss of PPHGI-1 leads to the generation of a Pph strain able to cause disease in the plant. In this study, we observed that lower bacterial densities inoculated into resistant bean (Phaseolus vulgaris) plants resulted in quicker PPHGI-1 loss from the population, and that loss of the island was strongly influenced by the type of plant resistance encountered by the bacteria. In addition, we found that a number of changes occurred in the bacterial genome during growth in the plant, whether or not PPHGI-1 was lost. We also present evidence that the circular PPHGI-1 episome is able to replicate autonomously when excised from the genome. These results shed more light onto the plasticity of the bacterial genome as it is influenced by in planta conditions.

Pseudomonas syringae pv. phaseolicola: from ‘has bean’ to supermodel

Pseudomonas syringae pv. phaseolicola causes halo blight of the common bean, Phaseolus vulgaris, worldwide and remains difficult to control. Races of the pathogen cause either disease symptoms or a resistant hypersensitive response on a series of differentially reacting bean cultivars. The molecular genetics of the interaction between P. syringae pv. phaseolicola and bean, and the evolution of bacterial virulence, have been investigated in depth and this research has led to important discoveries in the field of plant-microbe interactions. In this review, we discuss several of the areas of study that chart the rise of P. syringae pv. phaseolicola from a common pathogen of bean plants to a molecular plant-pathogen supermodel bacterium. Taxonomy: Bacteria; Proteobacteria, gamma subdivision; order Pseudomonadales; family Pseudomonadaceae; genus Pseudomonas; species Pseudomonas syringae; Genomospecies 2; pathogenic variety phaseolicola. Microbiological properties: Gram-negative, aerobic, motile, rod-shaped, 1.5 µm long, 0.7-1.2 µm in diameter, at least one polar flagellum, optimal temperatures for growth of 25-30 °C, oxidase negative, arginine dihydrolase negative, levan positive and elicits the hypersensitive response on tobacco. Host range: Major bacterial disease of common bean (Phaseolus vulgaris) in temperate regions and above medium altitudes in the tropics. Natural infections have been recorded on several other legume species, including all members of the tribe Phaseoleae with the exception of Desmodium spp. and Pisum sativum. Disease symptoms: Water-soaked lesions on leaves, pods, stems or petioles, that quickly develop greenish-yellow haloes on leaves at temperatures of less than 23 °C. Infected seeds may be symptomless, or have wrinkled or buttery-yellow patches on the seed coat. Seedling infection is recognized by general chlorosis, stunting and distortion of growth. Epidemiology: Seed borne and disseminated from exudation by water-splash and wind occurring during rainfall. Bacteria invade through wounds and natural openings (notably stomata). Weedy and cultivated alternative hosts may also harbour the bacterium. Disease control: Some measure of control is achieved with copper formulations and streptomycin. Pathogen-free seed and resistant cultivars are recommended. Useful websites: Pseudomonas-plant interaction http://www.pseudomonas-syringae.org/; PseudoDB http://xbase.bham.ac.uk/pseudodb/; Plant Associated and Environmental Microbes Database (PAMDB) http://genome.ppws.vt.edu/cgi-bin/MLST/home.pl; PseudoMLSA Database http://www.uib.es/microbiologiaBD/Welcome.html.

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.

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.

Potential influences on the United Kingdom's floods of winter 2013/14

During the winter of 2013/14, much of the UK experienced repeated intense rainfall events and flooding. This had a considerable impact on property and transport infrastructure. A key question is whether the burning of fossil fuels is changing the frequency of extremes, and if so to what extent. We assess the scale of the winter flooding before reviewing a broad range of Earth system drivers affecting UK rainfall. Some drivers can be potentially disregarded for these specific storms whereas others are likely to have increased their risk of occurrence. We discuss the requirements of hydrological models to transform rainfall into river flows and flooding. To determine any general changing flood risk, we argue that accurate modelling needs to capture evolving understanding of UK rainfall interactions with a broad set of factors. This includes changes to multiscale atmospheric, oceanic, solar and sea-ice features, and land-use and demographics. Ensembles of such model simulations may be needed to build probability distributions of extremes for both pre-industrial and contemporary concentration levels of atmospheric greenhouse gases.

Genetical and comparative genomics of Brassica under altered Ca supply identifies Arabidopsis Ca-transporter orthologs

Although Ca transport in plants is highly complex, the overexpression of vacuolar Ca2+ transporters in crops is a promising new technology to improve dietary Ca supplies through biofortification. Here, we sought to identify novel targets for increasing plant Ca accumulation using genetical and comparative genomics. Expression quantitative trait locus (eQTL) mapping to 1895 cis- and 8015 trans-loci were identified in shoots of an inbred mapping population of Brassica rapa (IMB211 × R500); 23 cis- and 948 trans-eQTLs responded specifically to altered Ca supply. eQTLs were screened for functional significance using a large database of shoot Ca concentration phenotypes of Arabidopsis thaliana. From 31 Arabidopsis gene identifiers tagged to robust shoot Ca concentration phenotypes, 21 mapped to 27 B. rapa eQTLs, including orthologs of the Ca2+ transporters At-CAX1 and At-ACA8. Two of three independent missense mutants of BraA.cax1a, isolated previously by targeting induced local lesions in genomes, have allele-specific shoot Ca concentration phenotypes compared with their segregating wild types. BraA.CAX1a is a promising target for altering the Ca composition of Brassica, consistent with prior knowledge from Arabidopsis. We conclude that multiple-environment eQTL analysis of complex crop genomes combined with comparative genomics is a powerful technique for novel gene identification/prioritization.

Infrared and millimetre-wave scintillometry in the suburban environment – Part 1: Structure parameters

Scintillometry, a form of ground-based remote sensing, provides the capability to estimate surface heat fluxes over scales of a few hundred metres to kilometres. Measurements are spatial averages, making this technique particularly valuable over areas with moderate heterogeneity such as mixed agricultural or urban environments. In this study, we present the structure parameters of temperature and humidity, which can be related to the sensible and latent heat fluxes through similarity theory, for a suburban area in the UK. The fluxes are provided in the second paper of this two-part series. A millimetre-wave scintillometer was combined with an infrared scintillometer along a 5.5 km path over northern Swindon. The pairing of these two wavelengths offers sensitivity to both temperature and humidity fluctuations, and the correlation between wavelengths is also used to retrieve the path-averaged temperature–humidity correlation. Comparison is made with structure parameters calculated from an eddy covariance station located close to the centre of the scintillometer path. The performance of the measurement techniques under different conditions is discussed. Similar behaviour is seen between the two data sets at sub-daily timescales. For the two summer-to-winter periods presented here, similar evolution is displayed across the seasons. A higher vegetation fraction within the scintillometer source area is consistent with the lower Bowen ratio observed (midday Bowen ratio < 1) compared with more built-up areas around the eddy covariance station. The energy partitioning is further explored in the companion paper.

The identification of genes important in pseudomonas syringae pv. phaseolicola plant colonisation using in vitro screening of transposon libraries

The bacterial plant pathogen Pseudomonas syringae pv. phaseolicola (Pph) colonises the surface of common bean plants before moving into the interior of plant tissue, via wounds and stomata. In the intercellular spaces the pathogen proliferates in the apoplastic fluid and forms microcolonies (biofilms) around plant cells. If the pathogen can suppress the plant’s natural resistance response, it will cause halo blight disease. The process of resistance suppression is fairly well understood, but the mechanisms used by the pathogen in colonisation are less clear. We hypothesised that we could apply in vitro genetic screens to look for changes in motility, colony formation, and adhesion, which are proxies for infection, microcolony formation and cell adhesion. We made transposon (Tn) mutant libraries of Pph strains 1448A and 1302A and found 106/1920 mutants exhibited alterations in colony morphology, motility and biofilm formation. Identification of the insertion point of the Tn identified within the genome highlighted, as expected, a number of altered motility mutants bearing mutations in genes encoding various parts of the flagellum. Genes involved in nutrient biosynthesis, membrane associated proteins, and a number of conserved hypothetical protein (CHP) genes were also identified. A mutation of one CHP gene caused a positive increase in in planta bacterial growth. This rapid and inexpensive screening method allows the discovery of genes important for in vitro traits that can be correlated to roles in the plant interaction