Altered expression of a raspberry homologue of VRN1 is associated with disruption of dormancy induction and misregulation of subsets of dormancy-associated genes.

Winter dormancy is a key process in the phenology of temperate perennials. Climate change is severely impacting its course leading to economic losses in agriculture. A better understanding of the underlying mechanisms, as well as the genetic basis of the different responses, is necessary for the development of climate-resilient cultivars. This study aims to provide an insight into winter dormancy in red raspberry (Rubus idaeus L). We report the transcriptomic profiles during dormancy in two raspberry cultivars with contrasting responses. The cultivar 'Glen Ample' showed a typical perennial phenology, whereas 'Glen Dee' registered consistent dormancy dysregulation, exhibiting active growth and flowering out of season. RNA-seq combined with weighted gene co-expression network analysis identified gene clusters in both genotypes that exhibited time-dependent expression profiles. Functional analysis of 'Glen Ample' gene clusters highlighted the significance of the cell and structural development prior to dormancy entry as well the role of genetic and epigenetic processes such as RNAi and DNA methylation in regulating gene expression. Dormancy release in 'Glen Ample' was associated with up-regulation of transcripts associated with the resumption of metabolism, nucleic acid biogenesis, and processing signal response pathways. Many of the processes occurring in 'Glen Ample' were dysregulated in 'Glen Dee' and 28 transcripts exhibiting time-dependent expression in 'Glen Ample' that also had an Arabidopsis homologue were not found in 'Glen Dee'. These included a gene with homology to Arabidopsis VRN1 (RiVRN1.1) that exhibited a sharp decline in expression following dormancy induction in 'Glen Ample'. Characterization of the gene region in the 'Glen Dee' genome revealed two large insertions upstream of the ATG start codon. We propose that expression below detection level of a specific VRN1 homologue in 'Glen Dee' causes dormancy misregulation as a result of inappropriate expression of a subset of genes that are directly or indirectly regulated by RiVRN1.1.

A chromosome conformation capture ordered sequence of the barley genome.

Cereal grasses of the Triticeae tribe have been the major food source in temperate regions since the dawn of agriculture. Their large genomes are characterized by a high content of repetitive elements and large pericentromeric regions that are virtually devoid of meiotic recombination. Here we present a high-quality reference genome assembly for barley (Hordeum vulgare L.). We use chromosome conformation capture mapping to derive the linear order of sequences across the pericentromeric space and to investigate the spatial organization of chromatin in the nucleus at megabase resolution. The composition of genes and repetitive elements differs between distal and proximal regions. Gene family analyses reveal lineage-specific duplications of genes involved in the transport of nutrients to developing seeds and the mobilization of carbohydrates in grains. We demonstrate the importance of the barley reference sequence for breeding by inspecting the genomic partitioning of sequence variation in modern elite germplasm, highlighting regions vulnerable to genetic erosion.

APETALA2 control of barley internode elongation.

Many plants dramatically elongate their stems during flowering, yet how this response is coordinated with the reproductive phase is unclear. We demonstrate that microRNA (miRNA) control of APETALA2 (AP2) is required for rapid, complete elongation of stem internodes in barley, especially of the final 'peduncle' internode directly underneath the inflorescence. Disrupted miR172 targeting of AP2 in the Zeo1.b barley mutant caused lower mitotic activity, delayed growth dynamics and premature lignification in the peduncle leading to fewer and shorter cells. Stage- and tissue-specific comparative transcriptomics between Zeo1.b and its parent cultivar showed reduced expression of proliferation-associated genes, ectopic expression of maturation-related genes and persistent, elevated expression of genes associated with jasmonate and stress responses. We further show that applying methyl jasmonate (MeJA) phenocopied the stem elongation of Zeo1.b, and that Zeo1.b itself was hypersensitive to inhibition by MeJA but less responsive to promotion by gibberellin. Taken together, we propose that miR172-mediated restriction of AP2 may modulate the jasmonate pathway to facilitate gibberellin-promoted stem growth during flowering.

TERMINAL FLOWER-1/CENTRORADIALIS inhibits tuberisation via protein interaction with the tuberigen activation complex.

Potato tuber formation is a secondary developmental programme by which cells in the subapical stolon region divide and radially expand to further differentiate into starch-accumulating parenchyma. Although some details of the molecular pathway that signals tuberisation are known, important gaps in our knowledge persist. Here, the role of a member of the TERMINAL FLOWER 1/CENTRORADIALIS gene family (termed StCEN) in the negative control of tuberisation is demonstrated for what is thought to be the first time. It is shown that reduced expression of StCEN accelerates tuber formation whereas transgenic lines overexpressing this gene display delayed tuberisation and reduced tuber yield. Protein-protein interaction studies (yeast two-hybrid and bimolecular fluorescence complementation) demonstrate that StCEN binds components of the recently described tuberigen activation complex. Using transient transactivation assays, we show that the StSP6A tuberisation signal is an activation target of the tuberigen activation complex, and that co-expression of StCEN blocks activation of the StSP6A gene by StFD-Like-1. Transcriptomic analysis of transgenic lines misexpressing StCEN identifies early transcriptional events in tuber formation. These results demonstrate that StCEN suppresses tuberisation by directly antagonising the function of StSP6A in stolons, identifying StCEN as a breeding marker to improve tuber initiation and yield through the selection of genotypes with reduced StCEN expression.

Chitosan primes plant defence mechanisms against Botrytis cinerea, including expression of Avr9/Cf-9 rapidly elicited genes.

Current crop protection strategies against the fungal pathogen Botrytis cinerea rely on a combination of conventional fungicides and host genetic resistance. However, due to pathogen evolution and legislation in the use of fungicides, these strategies are not sufficient to protect plants against this pathogen. Defence elicitors can stimulate plant defence mechanisms through a phenomenon known as defence priming. Priming results in a faster and/or stronger expression of resistance upon pathogen recognition by the host. This work aims to study defence priming by a commercial formulation of the elicitor chitosan. Treatments with chitosan result in induced resistance (IR) in solanaceous and brassicaceous plants. In tomato plants, enhanced resistance has been linked with priming of callose deposition and accumulation of the plant hormone jasmonic acid (JA). Large-scale transcriptomic analysis revealed that chitosan primes gene expression at early time-points after infection. In addition, two novel tomato genes with a characteristic priming profile were identified, Avr9/Cf-9 rapidly elicited protein 75 (ACRE75) and 180 (ACRE180). Transient and stable over-expression of ACRE75, ACRE180 and their Nicotiana benthamiana homologs, revealed that they are positive regulators of plant resistance against B. cinerea. This provides valuable information in the search for strategies to protect Solanaceae plants against B. cinerea.

BaRTv1.0: an improved barley reference transcript dataset to determine accurate changes in the barley transcriptome using RNA-seq.

BACKGROUND: The time required to analyse RNA-seq data varies considerably, due to discrete steps for computational assembly, quantification of gene expression and splicing analysis. Recent fast non-alignment tools such as Kallisto and Salmon overcome these problems, but these tools require a high quality, comprehensive reference transcripts dataset (RTD), which are rarely available in plants. RESULTS: A high-quality, non-redundant barley gene RTD and database (Barley Reference Transcripts - BaRTv1.0) has been generated. BaRTv1.0, was constructed from a range of tissues, cultivars and abiotic treatments and transcripts assembled and aligned to the barley cv. Morex reference genome (Mascher et al. Nature; 544: 427-433, 2017). Full-length cDNAs from the barley variety Haruna nijo (Matsumoto et al. Plant Physiol; 156: 20-28, 2011) determined transcript coverage, and high-resolution RT-PCR validated alternatively spliced (AS) transcripts of 86 genes in five different organs and tissue. These methods were used as benchmarks to select an optimal barley RTD. BaRTv1.0-Quantification of Alternatively Spliced Isoforms (QUASI) was also made to overcome inaccurate quantification due to variation in 5' and 3' UTR ends of transcripts. BaRTv1.0-QUASI was used for accurate transcript quantification of RNA-seq data of five barley organs/tissues. This analysis identified 20,972 significant differentially expressed genes, 2791 differentially alternatively spliced genes and 2768 transcripts with differential transcript usage. CONCLUSION: A high confidence barley reference transcript dataset consisting of 60,444 genes with 177,240 transcripts has been generated. Compared to current barley transcripts, BaRTv1.0 transcripts are generally longer, have less fragmentation and improved gene models that are well supported by splice junction reads. Precise transcript quantification using BaRTv1.0 allows routine analysis of gene expression and AS.

Identification of TIMING OF CAB EXPRESSION 1 as a temperature-sensitive negative regulator of tuberization in potato.

For many potato cultivars, tuber yield is optimal at average daytime temperatures in the range 14-22 °C. Above this range, tuber yield is reduced for most cultivars. We previously reported that moderately elevated temperature increases steady-state expression of the core circadian clock gene TIMING OF CAB EXPRESSION 1 (StTOC1) in developing tubers, whereas expression of the StSP6A tuberization signal is reduced, along with tuber yield. In this study we provide evidence that StTOC1 links environmental signalling with potato tuberization by suppressing StSP6A autoactivation in the stolons. We show that transgenic lines silenced in StTOC1 expression exhibit enhanced StSP6A transcript levels and changes in gene expression in developing tubers that are indicative of an elevated sink strength. Nodal cuttings of StTOC1 antisense lines displayed increased tuber yields at moderately elevated temperatures, whereas tuber yield and StSP6A expression were reduced in StTOC1 overexpressor lines. Here we identify a number of StTOC1 binding partners and demonstrate that suppression of StSP6A expression is independent of StTOC1 complex formation with the potato homolog StPIF3. Down-regulation of StTOC1 thus provides a strategy to mitigate the effects of elevated temperature on tuber yield.

Ancient barley landraces adapted to marginal soils demonstrate exceptional tolerance to manganese limitation.

BACKGROUND AND AIMS: Micronutrient deficiency in cereals is a problem of global significance, severely reducing grain yield and quality in marginal soils. Ancient landraces represent, through hundreds of years of local adaptation to adverse soil conditions, a unique reservoir of genes and unexplored traits for enhancing yield and abiotic stress tolerance. Here we explored and compared the genetic variation in a population of Northern European barley landraces and modern elite varieties, and their tolerance to manganese (Mn) limitation. METHODS: A total of 135 barley accessions were genotyped and the genetic diversity was explored using Neighbor-Joining clustering. Based on this analysis, a sub-population of genetically diverse landraces and modern elite control lines were evaluated phenotypically for their ability to cope with Mn-deficient conditions, across three different environments increasing in complexity from hydroponics through pot experiments to regional field trials. KEY RESULTS: Genetically a group of Scottish barley landraces (Bere barley) were found to cluster according to their island of origin, and accessions adapted to distinct biogeographical zones with reduced soil fertility had particularly larger Mn, but also zinc (Zn) and copper (Cu) concentrations in the shoot. Strikingly, when grown in an alkaline sandy soil in the field, the locally adapted landraces demonstrated an exceptional ability to acquire and translocate Mn to developing leaves, maintain photosynthesis and generate robust grain yields, whereas modern elite varieties totally failed to complete their life cycle. CONCLUSIONS: Our results highlight the importance of gene pools of local adaptation and the value of ancient landrace material to identify and characterize genes that control nutrient use efficiency traits in adverse environments to raise future crop production and improve agricultural sustainability in marginal soils. We propose and discuss a model summarizing the physiological mechanisms involved in the complex trait of tolerance to Mn limitation.

RXLR Effector AVR2 Up-Regulates a Brassinosteroid-Responsive bHLH Transcription Factor to Suppress Immunity.

An emerging area in plant research focuses on antagonism between regulatory systems governing growth and immunity. Such cross talk represents a point of vulnerability for pathogens to exploit. AVR2, an RXLR effector secreted by the potato blight pathogen Phytophthora infestans, interacts with potato BSL1, a putative phosphatase implicated in growth-promoting brassinosteroid (BR) hormone signaling. Transgenic potato (Solanum tuberosum) plants expressing the effector exhibit transcriptional and phenotypic hallmarks of overactive BR signaling and show enhanced susceptibility to P. infestans Microarray analysis was used to identify a set of BR-responsive marker genes in potato, all of which are constitutively expressed to BR-induced levels in AVR2 transgenic lines. One of these genes was a bHLH transcription factor, designated StCHL1, homologous to AtCIB1 and AtHBI1, which are known to facilitate antagonism between BR and immune responses. Transient expression of either AVR2 or CHL1 enhanced leaf colonization by P. infestans and compromised immune cell death activated by perception of the elicitin Infestin1 (INF1). Knockdown of CHL1 transcript using Virus-Induced Gene Silencing (VIGS) reduced colonization of P. infestans on Nicotiana benthamiana Moreover, the ability of AVR2 to suppress INF1-triggered cell death was attenuated in NbCHL1-silenced plants, indicating that NbCHL1 was important for this effector activity. Thus, AVR2 exploits cross talk between BR signaling and innate immunity in Solanum species, representing a novel, indirect mode of innate immune suppression by a filamentous pathogen effector.

The redox state of the apoplast influences the acclimation of photosynthesis and leaf metabolism to changing irradiance.

The redox state of the apoplast is largely determined by ascorbate oxidase (AO) activity. The influence of AO activity on leaf acclimation to changing irradiance was explored in wild-type (WT) and transgenic tobacco (Nicotiana tobaccum) lines containing either high [pumpkin AO (PAO)] or low [tobacco AO (TAO)] AO activity at low [low light (LL); 250 µmol m-2  s-1 ] and high [high light (HL); 1600 µmol m-2  s-1 ] irradiance and following the transition from HL to LL. AO activities changed over the photoperiod, particularly in the PAO plants. AO activity had little effect on leaf ascorbate, which was significantly higher under HL than under LL. Apoplastic ascorbate/dehydroascorbate (DHA) ratios and threonate levels were modified by AO activity. Despite decreased levels of transcripts encoding ascorbate synthesis enzymes, leaf ascorbate increased over the first photoperiod following the transition from HL to LL, to much higher levels than LL-grown plants. Photosynthesis rates were significantly higher in the TAO leaves than in WT or PAO plants grown under HL but not under LL. Sub-sets of amino acids and fatty acids were lower in TAO and WT leaves than in the PAO plants under HL, and following the transition to LL. Light acclimation processes are therefore influenced by the apoplastic as well as chloroplastic redox state.

A comparative analysis of nonhost resistance across the two Triticeae crop species wheat and barley.

BACKGROUND: Nonhost resistance (NHR) protects plants against a vast number of non-adapted pathogens which implicates a potential exploitation as source for novel disease resistance strategies. Aiming at a fundamental understanding of NHR a global analysis of transcriptome reprogramming in the economically important Triticeae cereals wheat and barley, comparing host and nonhost interactions in three major fungal pathosystems responsible for powdery mildew (Blumeria graminis ff. ssp.), cereal blast (Magnaporthe sp.) and leaf rust (Puccinia sp.) diseases, was performed. RESULTS: In each pathosystem a significant transcriptome reprogramming by adapted- or non-adapted pathogen isolates was observed, with considerable overlap between Blumeria, Magnaporthe and Puccinia. Small subsets of these general pathogen-regulated genes were identified as differentially regulated between host and corresponding nonhost interactions, indicating a fine-tuning of the general pathogen response during the course of co-evolution. Additionally, the host- or nonhost-related responses were rather specific for each pair of adapted and non-adapted isolates, indicating that the nonhost resistance-related responses were to a great extent pathosystem-specific. This pathosystem-specific reprogramming may reflect different resistance mechanisms operating against non-adapted pathogens with different lifestyles, or equally, different co-option of the hosts by the adapted isolates to create an optimal environment for infection. To compare the transcriptional reprogramming between wheat and barley, putative orthologues were identified. Within the wheat and barley general pathogen-regulated genes, temporal expression profiles of orthologues looked similar, indicating conserved general responses in Triticeae against fungal attack. However, the comparison of orthologues differentially expressed between host and nonhost interactions revealed fewer commonalities between wheat and barley, but rather suggested different host or nonhost responses in the two cereal species. CONCLUSIONS: Taken together, our results suggest independent co-evolutionary forces acting on host pathosystems mirrored by barley- or wheat-specific nonhost responses. As a result of evolutionary processes, at least for the pathosystems investigated, NHR appears to rely on rather specific plant responses.

The Dynamics of Transcript Abundance during Cellularization of Developing Barley Endosperm.

Within the cereal grain, the endosperm and its nutrient reserves are critical for successful germination and in the context of grain utilization. The identification of molecular determinants of early endosperm development, particularly regulators of cell division and cell wall deposition, would help predict end-use properties such as yield, quality, and nutritional value. Custom microarray data have been generated using RNA isolated from developing barley grain endosperm 3 d to 8 d after pollination (DAP). Comparisons of transcript abundance over time revealed 47 gene expression modules that can be clustered into 10 broad groups. Superimposing these modules upon cytological data allowed patterns of transcript abundance to be linked with key stages of early grain development. Here, attention was focused on how the datasets could be mined to explore and define the processes of cell wall biosynthesis, remodeling, and degradation. Using a combination of spatial molecular network and gene ontology enrichment analyses, it is shown that genes involved in cell wall metabolism are found in multiple modules, but cluster into two main groups that exhibit peak expression at 3 DAP to 4 DAP and 5 DAP to 8 DAP. The presence of transcription factor genes in these modules allowed candidate genes for the control of wall metabolism during early barley grain development to be identified. The data are publicly available through a dedicated web interface (https://ics.hutton.ac.uk/barseed/), where they can be used to interrogate co- and differential expression for any other genes, groups of genes, or transcription factors expressed during early endosperm development.

Molecular and Biochemical Examination of Spraing Disease in Potato Tuber in Response to Tobacco rattle virus Infection.

Field-grown tubers of potato were examined for infection by Tobacco rattle virus (TRV) and consequent production of corky ringspot or spraing symptoms. A microarray study identified genes that are differentially expressed in tuber tissue in response to TRV infection and to spraing production, suggesting that hypersensitive response (HR) pathways are activated in spraing-symptomatic tubers. This was confirmed by quantitative reverse transcription-polymerase chain reaction (qRT-PCR) of a selected group of HR-related genes and by histochemical staining of excised tuber tissue with spraing symptoms. qRT-PCR of TRV in different regions of the same tuber slice showed that nonsymptomatic areas contained higher levels of virus relative to spraing-symptomatic areas. This suggests that spraing formation is associated with an active plant defense that reduces the level of virus in the infected tuber. Expression of two of the same plant defense genes was similarly upregulated in tubers that were infected with Potato mop-top virus, a virus that also induces spraing formation.

The low-recombining pericentromeric region of barley restricts gene diversity and evolution but not gene expression.

The low-recombining pericentromeric region of the barley genome contains roughly a quarter of the genes of the species, embedded in low-recombining DNA that is rich in repeats and repressive chromatin signatures. We have investigated the effects of pericentromeric region residency upon the expression, diversity and evolution of these genes. We observe no significant difference in average transcript level or developmental RNA specificity between the barley pericentromeric region and the rest of the genome. In contrast, all of the evolutionary parameters studied here show evidence of compromised gene evolution in this region. First, genes within the pericentromeric region of wild barley show reduced diversity and significantly weakened purifying selection compared with the rest of the genome. Second, gene duplicates (ohnolog pairs) derived from the cereal whole-genome duplication event ca. 60MYa have been completely eliminated from the barley pericentromeric region. Third, local gene duplication in the pericentromeric region is reduced by 29% relative to the rest of the genome. Thus, the pericentromeric region of barley is a permissive environment for gene expression but has restricted gene evolution in a sizeable fraction of barley's genes.

A sequence-ready physical map of barley anchored genetically by two million single-nucleotide polymorphisms.

Barley (Hordeum vulgare) is an important cereal crop and a model species for Triticeae genomics. To lay the foundation for hierarchical map-based sequencing, a genome-wide physical map of its large and complex 5.1 billion-bp genome was constructed by high-information content fingerprinting of almost 600,000 bacterial artificial chromosomes representing 14-fold haploid genome coverage. The resultant physical map comprises 9,265 contigs with a cumulative size of 4.9 Gb representing 96% of the physical length of the barley genome. The reliability of the map was verified through extensive genetic marker information and the analysis of topological networks of clone overlaps. A minimum tiling path of 66,772 minimally overlapping clones was defined that will serve as a template for hierarchical clone-by-clone map-based shotgun sequencing. We integrated whole-genome shotgun sequence data from the individuals of two mapping populations with published bacterial artificial chromosome survey sequence information to genetically anchor the physical map. This novel approach in combination with the comprehensive whole-genome shotgun sequence data sets allowed us to independently validate and improve a previously reported physical and genetic framework. The resources developed in this study will underpin fine-mapping and cloning of agronomically important genes and the assembly of a draft genome sequence.

The MORPH-R web server and software tool for predicting missing genes in biological pathways.

A biological pathway is the set of molecular entities involved in a given biological process and the interrelations among them. Even though biological pathways have been studied extensively, discovering missing genes in pathways remains a fundamental challenge. Here, we present an easy-to-use tool that allows users to run MORPH (MOdule-guided Ranking of candidate PatHway genes), an algorithm for revealing missing genes in biological pathways, and demonstrate its capabilities. MORPH supports the analysis in tomato, Arabidopsis and the two new species: rice and the newly sequenced potato genome. The new tool, called MORPH-R, is available both as a web server (at http://bioinformatics.psb.ugent.be/webtools/morph/) and as standalone software that can be used locally. In the standalone version, the user can apply the tool to new organisms using any proprietary and public data sources.

Barley whole exome capture: a tool for genomic research in the genus Hordeum and beyond.

Advanced resources for genome-assisted research in barley (Hordeum vulgare) including a whole-genome shotgun assembly and an integrated physical map have recently become available. These have made possible studies that aim to assess genetic diversity or to isolate single genes by whole-genome resequencing and in silico variant detection. However such an approach remains expensive given the 5 Gb size of the barley genome. Targeted sequencing of the mRNA-coding exome reduces barley genomic complexity more than 50-fold, thus dramatically reducing this heavy sequencing and analysis load. We have developed and employed an in-solution hybridization-based sequence capture platform to selectively enrich for a 61.6 megabase coding sequence target that includes predicted genes from the genome assembly of the cultivar Morex as well as publicly available full-length cDNAs and de novo assembled RNA-Seq consensus sequence contigs. The platform provides a highly specific capture with substantial and reproducible enrichment of targeted exons, both for cultivated barley and related species. We show that this exome capture platform provides a clear path towards a broader and deeper understanding of the natural variation residing in the mRNA-coding part of the barley genome and will thus constitute a valuable resource for applications such as mapping-by-sequencing and genetic diversity analyzes.

An evaluation of genotyping by sequencing (GBS) to map the Breviaristatum-e (ari-e) locus in cultivated barley.

UNLABELLED: We explored the use of genotyping by sequencing (GBS) on a recombinant inbred line population (GPMx) derived from a cross between the two-rowed barley cultivar 'Golden Promise' (ari-e.GP/Vrs1) and the six-rowed cultivar 'Morex' (Ari-e/vrs1) to map plant height. We identified three Quantitative Trait Loci (QTL), the first in a region encompassing the spike architecture gene Vrs1 on chromosome 2H, the second in an uncharacterised centromeric region on chromosome 3H, and the third in a region of chromosome 5H coinciding with the previously described dwarfing gene Breviaristatum-e (Ari-e). BACKGROUND: Barley cultivars in North-western Europe largely contain either of two dwarfing genes; Denso on chromosome 3H, a presumed ortholog of the rice green revolution gene OsSd1, or Breviaristatum-e (ari-e) on chromosome 5H. A recessive mutant allele of the latter gene, ari-e.GP, was introduced into cultivation via the cv. 'Golden Promise' that was a favourite of the Scottish malt whisky industry for many years and is still used in agriculture today. RESULTS: Using GBS mapping data and phenotypic measurements we show that ari-e.GP maps to a small genetic interval on chromosome 5H and that alternative alleles at a region encompassing Vrs1 on 2H along with a region on chromosome 3H also influence plant height. The location of Ari-e is supported by analysis of near-isogenic lines containing different ari-e alleles. We explored use of the GBS to populate the region with sequence contigs from the recently released physically and genetically integrated barley genome sequence assembly as a step towards Ari-e gene identification. CONCLUSIONS: GBS was an effective and relatively low-cost approach to rapidly construct a genetic map of the GPMx population that was suitable for genetic analysis of row type and height traits, allowing us to precisely position ari-e.GP on chromosome 5H. Mapping resolution was lower than we anticipated. We found the GBS data more complex to analyse than other data types but it did directly provide linked SNP markers for subsequent higher resolution genetic analysis.

Proteomics and transcriptomics of the BABA-induced resistance response in potato using a novel functional annotation approach.

BACKGROUND: Induced resistance (IR) can be part of a sustainable plant protection strategy against important plant diseases. ß-aminobutyric acid (BABA) can induce resistance in a wide range of plants against several types of pathogens, including potato infected with Phytophthora infestans. However, the molecular mechanisms behind this are unclear and seem to be dependent on the system studied. To elucidate the defence responses activated by BABA in potato, a genome-wide transcript microarray analysis in combination with label-free quantitative proteomics analysis of the apoplast secretome were performed two days after treatment of the leaf canopy with BABA at two concentrations, 1 and 10 mM. RESULTS: Over 5000 transcripts were differentially expressed and over 90 secretome proteins changed in abundance indicating a massive activation of defence mechanisms with 10 mM BABA, the concentration effective against late blight disease. To aid analysis, we present a more comprehensive functional annotation of the microarray probes and gene models by retrieving information from orthologous gene families across 26 sequenced plant genomes. The new annotation provided GO terms to 8616 previously un-annotated probes. CONCLUSIONS: BABA at 10 mM affected several processes related to plant hormones and amino acid metabolism. A major accumulation of PR proteins was also evident, and in the mevalonate pathway, genes involved in sterol biosynthesis were down-regulated, whereas several enzymes involved in the sesquiterpene phytoalexin biosynthesis were up-regulated. Interestingly, abscisic acid (ABA) responsive genes were not as clearly regulated by BABA in potato as previously reported in Arabidopsis. Together these findings provide candidates and markers for improved resistance in potato, one of the most important crops in the world.