The evolutionary patterns of barley pericentromeric chromosome regions, as shaped by linkage disequilibrium and domestication

The distribution of recombination events along large cereal chromosomes is uneven and is generally restricted to gene-rich telomeric ends. To understand how the lack of recombination affects diversity in the large pericentromeric regions, we analysed deep exome capture data from a final panel of 815 Hordeum vulgare (barley) cultivars, landraces and wild barleys, sampled from across their eco-geographical ranges. We defined and compared variant data across the pericentromeric and non-pericentromeric regions, observing a clear partitioning of diversity both within and between chromosomes and germplasm groups. Dramatically reduced diversity was found in the pericentromeres of both cultivars and landraces when compared with wild barley. We observed a mixture of completely and partially differentiated single-nucleotide polymorphisms (SNPs) between domesticated and wild gene pools, suggesting that domesticated gene pools were derived from multiple wild ancestors. Patterns of genome-wide linkage disequilibrium, haplotype block size and number, and variant frequency within blocks showed clear contrasts among individual chromosomes and between cultivars and wild barleys. Although most cultivar chromosomes shared a single major pericentromeric haplotype, chromosome 7H clearly differentiated the two-row and six-row types associated with different geographical origins. Within the pericentromeric regions we identified 22 387 non-synonymous SNPs, 92 of which were fixed for alternative alleles in cultivar versus wild accessions. Surprisingly, only 29 SNPs found exclusively in the cultivars were predicted to be 'highly deleterious'. Overall, our data reveal an unconventional pericentromeric genetic landscape among distinct barley gene pools, with different evolutionary processes driving domestication and diversification.

A reversible light- and genotype-dependent acquired thermotolerance response protects the potato plant from damage due to excessive temperature

A powerful acquired thermotolerance response in potato was demonstrated and characterised in detail, showing the time course required for tolerance, the reversibility of the process and requirement for light. Potato is particularly vulnerable to increased temperature, considered to be the most important uncontrollable factor affecting growth and yield of this globally significant crop. Here, we describe an acquired thermotolerance response in potato, whereby treatment at a mildly elevated temperature primes the plant for more severe heat stress. We define the time course for acquiring thermotolerance and demonstrate that light is essential for the process. In all four commercial tetraploid cultivars that were tested, acquisition of thermotolerance by priming was required for tolerance at elevated temperature. Accessions from several wild-type species and diploid genotypes did not require priming for heat tolerance under the test conditions employed, suggesting that useful variation for this trait exists. Physiological, transcriptomic and metabolomic approaches were employed to elucidate potential mechanisms that underpin the acquisition of heat tolerance. This analysis indicated a role for cell wall modification, auxin and ethylene signalling, and chromatin remodelling in acclimatory priming resulting in reduced metabolic perturbation and delayed stress responses in acclimated plants following transfer to 40 °C.

Correction to: A reversible light- and genotype-dependent acquired thermotolerance response protects the potato plant from damage due to excessive temperature

The article A reversible light- and genotype-dependent acquired thermotolerance response protects the potato plant from damage due to excessive temperature, written by Almudena Trapero-Mozos, Laurence J. M. Ducreux, Craita E. Bita, Wayne Morris, Cosima Wiese, Jenny A. Morris, Christy Paterson, Peter E. Hedley, Robert D. Hancock, and Mark Taylor, was originally published electronically on the publisher's internet portal (currently SpringerLink) on 8 March 2018 without open access. With the author(s)' decision to opt for Open Choice the copyright of the article changed on 30 April 2018 to © The Author(s) 2018 and the article is forthwith distributed under the terms of the Creative Commons Attribution 4.0 International License ( http://creativecommons.org/licenses/by/4.0/ ), which permits use, duplication, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license and indicate if changes were made.The ​original ​article ​has ​been ​corrected.

Redox Control of Aphid Resistance through Altered Cell Wall Composition and Nutritional Quality

The mechanisms underpinning plant perception of phloem-feeding insects, particularly aphids, remain poorly characterized. Therefore, the role of apoplastic redox state in controlling aphid infestation was explored using transgenic tobacco (Nicotiana tabacum) plants that have either high (PAO) or low (TAO) ascorbate oxidase (AO) activities relative to the wild type. Only a small number of leaf transcripts and metabolites were changed in response to genotype, and cell wall composition was largely unaffected. Aphid fecundity was decreased significantly in TAO plants compared with other lines. Leaf sugar levels were increased and maximum extractable AO activities were decreased in response to aphids in all genotypes. Transcripts encoding the Respiratory Burst Oxidase Homolog F, signaling components involved in ethylene and other hormone-mediated pathways, photosynthetic electron transport components, sugar, amino acid, and cell wall metabolism, were increased significantly in the TAO plants in response to aphid perception relative to other lines. The levels of galactosylated xyloglucan were decreased significantly in response to aphid feeding in all the lines, the effect being the least in the TAO plants. Similarly, all lines exhibited increases in tightly bound (1→4)-β-galactan. Taken together, these findings identify AO-dependent mechanisms that limit aphid infestation.

Assessing the genetic diversity of rice originating from Bangladesh, Assam and West Bengal

BACKGROUND

Genetic diversity among rice cultivars from Bangladesh and North East India was assessed using a custom 384-SNP microarray assay. A total of 511 cultivars were obtained from several sources, choosing landraces likely to be from the aus subpopulation and modern improved cultivars from Bangladesh. Cultivars from the OryzaSNP set and Rice Diversity Panel 1 (RDP1) were also included for reference.

RESULTS

The population analysis program STRUCTURE was used to infer putative population groups in the panel, revealing four groups: indica (76 cultivars), japonica (55) and two distinct groups within the aus subpopulation (aus-1 = 99, aus-2 = 151). Principal Component Analysis was used to confirm the four population groups identified by STRUCTURE. The analysis revealed cultivars that belonged to neither aus-1 nor aus-2 but which are clearly aus based on the combined probabilities of their membership of the two aus groups which have been termed aus-admix (96). Information obtained from the panel of 511 cultivars was used to assign rice groups to 74 additional landraces obtained from Assam and West Bengal. While both the aus-1 and aus-2 groups were represented approximately equally in India, aus-2 (which includes cultivar N 22) was more common in Bangladesh, but was not found at all in West Bengal.

CONCLUSIONS

Examining the distribution of landrace names within theaus-1 and aus-2 groups suggests that aus-1 is associated with the term "boro", a word used to describe a winter growing season in Bangladesh and Assam. The information described here has been used to select a population of 300 cultivars for Genome Wide Association studies of the aus rice subpopulation.

WHIRLY1 Functions in the Control of Responses to Nitrogen Deficiency But Not Aphid Infestation in Barley

WHIRLY1 is largely targeted to plastids, where it is a major constituent of the nucleoids. To explore WHIRLY1 functions in barley (Hordeum vulgare), RNA interference-knockdown lines (W1-1, W1-7, and W1-9) that have very low levels of HvWHIRLY1 transcripts were characterized in plants grown under optimal and stress conditions. The WHIRLY1-1 (W1-1), W1-7, and W1-9 plants were phenotypically similar to the wild type but produced fewer tillers and seeds. Photosynthesis rates were similar in all lines, but W1-1, W1-7, and W1-9 leaves had significantly more chlorophyll and less sucrose than the wild type. Transcripts encoding specific subsets of chloroplast-localized proteins, such as ribosomal proteins, subunits of the RNA polymerase, and thylakoid nicotinamide adenine dinucleotide (reduced) and cytochrome b6/f complexes, were much more abundant in the W1-7 leaves than the wild type. Although susceptibility of aphid (Myzus persicae) infestation was similar in all lines, the WHIRLY1-deficient plants showed altered responses to nitrogen deficiency, maintaining higher photosynthetic CO2 assimilation rates than the wild type under limiting nitrogen. Although all lines showed globally similar low nitrogen-dependent changes in transcripts and metabolites, the increased abundance of FAR-RED IMPAIRED RESPONSE1-like transcripts in nitrogen-deficient W1-7 leaves infers that WHIRLY1 has a role in communication between plastid and nuclear genes encoding photosynthetic proteins during abiotic stress.

Nitrogen deficiency in barley (Hordeum vulgare) seedlings induces molecular and metabolic adjustments that trigger aphid resistance

Agricultural nitrous oxide (N2O) pollution resulting from the use of synthetic fertilizers represents a significant contribution to anthropogenic greenhouse gas emissions, providing a rationale for reduced use of nitrogen (N) fertilizers. Nitrogen limitation results in extensive systems rebalancing that remodels metabolism and defence processes. To analyse the regulation underpinning these responses, barley (Horedeum vulgare) seedlings were grown for 7 d under N-deficient conditions until net photosynthesis was 50% lower than in N-replete controls. Although shoot growth was decreased there was no evidence for the induction of oxidative stress despite lower total concentrations of N-containing antioxidants. Nitrogen-deficient barley leaves were rich in amino acids, sugars and tricarboxylic acid cycle intermediates. In contrast to N-replete leaves one-day-old nymphs of the green peach aphid (Myzus persicae) failed to reach adulthood when transferred to N-deficient barley leaves. Transcripts encoding cell, sugar and nutrient signalling, protein degradation and secondary metabolism were over-represented in N-deficient leaves while those associated with hormone metabolism were similar under both nutrient regimes with the exception of mRNAs encoding proteins involved in auxin metabolism and responses. Significant similarities were observed between the N-limited barley leaf transcriptome and that of aphid-infested Arabidopsis leaves. These findings not only highlight significant similarities between biotic and abiotic stress signalling cascades but also identify potential targets for increasing aphid resistance with implications for the development of sustainable agriculture.

Towards an understanding of the control of 'crumbly' fruit in red raspberry

The genetic disorder known as ‘crumbly’ fruit is becoming a serious problem in the European raspberry industry. The study set out to examine the crumbly phenotype in a red raspberry mapping population under two environments (field and polytunnel) across six seasons in an effort to understand variability of the syndrome and to examine whether genetic factors were important and if so, whether QTL associated with the phenotype could be identified. This highlighted that seasonal, environmental (field or polytunnel) and genetic factors all influence the condition. Two QTL that are important for the genetic control of the condition have been located on linkage groups one and three, and an association with ripening time has been identified.

Characterization of Arabidopsis Transcriptional Responses to Different Aphid Species Reveals Genes that Contribute to Host Susceptibility and Non-host Resistance

Aphids are economically important pests that display exceptional variation in host range. The determinants of diverse aphid host ranges are not well understood, but it is likely that molecular interactions are involved. With significant progress being made towards understanding host responses upon aphid attack, the mechanisms underlying non-host resistance remain to be elucidated. Here, we investigated and compared Arabidopsis thaliana host and non-host responses to aphids at the transcriptional level using three different aphid species, Myzus persicae, Myzus cerasi and Rhopalosiphum pisum. Gene expression analyses revealed a high level of overlap in the overall gene expression changes during the host and non-host interactions with regards to the sets of genes differentially expressed and the direction of expression changes. Despite this overlap in transcriptional responses across interactions, there was a stronger repression of genes involved in metabolism and oxidative responses specifically during the host interaction with M. persicae. In addition, we identified a set of genes with opposite gene expression patterns during the host versus non-host interactions. Aphid performance assays on Arabidopsis mutants that were selected based on our transcriptome analyses identified novel genes contributing to host susceptibility, host defences during interactions with M. persicae as well to non-host resistance against R. padi. Understanding how plants respond to aphid species that differ in their ability to infest plant species, and identifying the genes and signaling pathways involved, is essential for the development of novel and durable aphid control in crop plants.

Aphids are phloem-feeding insects that cause feeding damage and transmit plant viruses to many crops. While most aphid species are restricted to one or few host plants, some aphids can infest a wide range of plant species. These insects spend a considerable time on non-hosts, where they probe the leaf tissue and secrete saliva, but for unknown reasons are unable to ingest phloem sap. This suggests that aphids interact with non-host plants at the molecular level, but potentially do not suppress plant defences and/or promote the release of nutrients. We compared gene expression of plants during host and non-host interactions with aphids to identify genes involved in immunity. We found significant overlap in the plant responses to aphids regardless of the type of interaction. Despite this, we identified a set of genes specifically affected during host or non-host interactions with specific aphid species. In addition, we showed that several of these genes contribute to host and/or non-host immunity. These findings are important, as they advance our understanding of the plant cellular processes involved in host and non-host responses against insect pests. Understanding mechanisms of host and non-host resistance to plant parasites is essential for development of novel control strategies.

tropiTree: an NGS-based EST-SSR resource for 24 tropical tree species

The development of genetic tools for non-model organisms has been hampered by cost, but advances in next-generation sequencing (NGS) have created new opportunities. In ecological research, this raises the prospect for developing molecular markers to simultaneously study important genetic processes such as gene flow in multiple non-model plant species within complex natural and anthropogenic landscapes. Here, we report the use of bar-coded multiplexed paired-end Illumina NGS for the de novo development of expressed sequence tag-derived simple sequence repeat (EST-SSR) markers at low cost for a range of 24 tree species. Each chosen tree species is important in complex tropical agroforestry systems where little is currently known about many genetic processes. An average of more than 5,000 EST-SSRs was identified for each of the 24 sequenced species, whereas prior to analysis 20 of the species had fewer than 100 nucleotide sequence citations. To make results available to potential users in a suitable format, we have developed an open-access, interactive online database, tropiTree (http://bioinf.hutton.ac.uk/tropiTree), which has a range of visualisation and search facilities, and which is a model for the efficient presentation and application of NGS data.

Vitamin C and the abscisic acid-insensitive 4 transcription factor are important determinants of aphid resistance in Arabidopsis

AIMS

Aphids, like other insects, are probably unable to synthesize vitamin C (ascorbic acid), which is therefore an essential dietary nutrient that has to be obtained from the host plant. Plant responses to aphids involve hormones such as salicylic acid (SA), jasmonic acid (JA), and abscisic acid (ABA), but hormone/redox interactions remain poorly characterized. We therefore investigated hormone/redox signaling in the response of Arabidopsis thaliana to infestation by the aphid Myzus persicae, focusing on the interactions between ascorbic acid and ABA, together with the influence of altered ascorbate and ABA signaling on the SA- and JA-dependent pathways.

RESULTS

Whole-genome microarray analysis revealed highly dynamic transcriptional responses to aphid infestation with extensive differences between transcript profiles of infested and systemic leaves, revealing aphid-dependent effects on the suites of transcripts involved in the redox, SA, and ABA responses. Central roles for ascorbate, ABA-insensitive 4 (ABI4), and oxidative signal-inducible 1 in plant resistance to aphids were demonstrated by altered fecundity on respective mutants. However, ABA had a negative effect on aphid resistance, as did ABI4 or redox-responsive transcription factor 1. The decrease in aphid fecundity observed in mutants defective in ascorbate accumulation (vtc2) was absent from abi4vtc2 double mutants that are also deficient in ABA signaling (abi4). Aphid-dependent transcriptome responses reveal a role for ascorbate-regulated receptor-like kinases in plant defenses against aphids.

INNOVATION

Vitamin C deficiency enhances plant resistance to aphids through redox signaling pathways rather than dietary requirements.

CONCLUSION

ABI4 is a linchpin of redox regulation of the innate immune response to aphids.

A metabolic regulator modulates virulence and quorum sensing signal production in Pectobacterium atrosepticum

Plant cell wall-degrading enzymes (PCWDE) are key virulence determinants in the pathogenesis of the potato pathogen Pectobacterium atrosepticum. In this study, we report the impact on virulence of a transposon insertion mutation in the metJ gene that codes for the repressor of the methionine biosynthesis regulon. In a mutant strain defective for the small regulatory RNA rsmB, PCWDE are not produced and virulence in potato tubers is almost totally abolished. However, when the metJ gene is disrupted in this background, the rsmB(-) phenotype is suppressed and virulence and PCWDE production are restored. Additionally, when metJ is disrupted, production of the quorum-sensing signal, N-(3-oxohexanoyl)-homoserine lactone, is increased. The metJ mutant strains showed pleiotropic transcriptional impacts affecting approximately a quarter of the genome. Genes involved in methionine biosynthesis were most highly upregulated but many virulence-associated transcripts were also upregulated. This is the first report of the impact of the MetJ repressor on virulence in bacteria.

An exceptionally high nucleotide and haplotype diversity and a signature of positive selection for the eIF4E resistance gene in barley are revealed by allele mining and phylogenetic analyses of natural populations

In barley, the eukaryotic translation initiation factor 4E (eIF4E) gene situated on chromosome 3H is recognized as an important source of resistance to the bymoviruses Barley yellow mosaic virus and Barley mild mosaic virus. In modern barley cultivars, two recessive eIF4E alleles, rym4 and rym5, confer different isolate-specific resistances. In this study, the sequence of eIF4E was analysed in 1090 barley landraces and noncurrent cultivars originating from 84 countries. An exceptionally high nucleotide diversity was evident in the coding sequence of eIF4E but not in either the adjacent MCT-1 gene or the sequence-related eIF(iso)4E gene situated on chromosome 1H. Surprisingly, all nucleotide polymorphisms detected in the coding sequence of eIF4E resulted in amino acid changes. A total of 47 eIF4E haplotypes were identified, and phylogenetic analysis using maximum likelihood provided evidence of strong positive selection acting on this barley gene. The majority of eIF4E haplotypes were found to be specific to distinct geographic regions. Furthermore, the eI4FE haplotype diversity (uh) was found to be considerably higher in East Asia, whereas SNP genotyping identified a comparatively low degree of genome-wide genetic diversity in 16 of 17 tested accessions (each carrying a different eIF4E haplotype) from this same region. In addition, selection statistic calculations using coalescent simulations showed evidence of non-neutral variation for eIF4E in several geographic regions, including East Asia, the region with a long history of the bymovirus-induced yellow mosaic disease. Together these findings suggest that eIF4E may play a role in barley adaptation to local habitats.

Highly parallel gene-to-BAC addressing using microarrays

Second-generation sequencing now provides the potential for low-cost generation of whole-genome sequences. However, for large-genome organisms with high repetitive DNA content, genome-wide short read sequence assembly is currently impossible, with accurate ordering and localization of genes still relying heavily on integration with physical and genetic maps. To facilitate this process, we have used Agilent microarrays to simultaneously address thousands of gene sequences to individual BAC clones and contiguous sequences that form part of an emerging physical map of the large and currently unsequenced 5.3-Gb barley genome. The approach represents a cost-effective, highly parallel alternative to traditional addressing methods. By coupling the gene-to-BAC address data with gene-based molecular markers, thousands of BACs can be anchored directly to the genetic map, thereby generating a framework for orientating and ordering genes, and providing direct links to phenotypic traits.

Combining genetical genomics and bulked segregant analysis-based differential expression: an approach to gene localization

Positional gene isolation in unsequenced species generally requires either a reference genome sequence or an inference of gene content based on conservation of synteny with a genomic model. In the large unsequenced genomes of the Triticeae cereals the latter, i.e. conservation of synteny with the rice and Brachypodium genomes, provides a powerful proxy for establishing local gene content and order. However, efficient exploitation of conservation of synteny requires 'homology bridges' between the model genome and the target region that contains a gene of interest. As effective homology bridges are generally the sequences of genetically mapped genes, increasing the density of these genes around a target locus is an important step in the process. We used bulked segregant analysis (BSA) of transcript abundance data to identify genes located in a specific region of the barley genome. The approach is valuable because only a relatively small proportion of barley genes are currently placed on a genetic map. We analyzed eQTL datasets from the reference Steptoe × Morex doubled haploid population and showed a strong association between differential gene expression and cis-regulation, with 83% of differentially expressed genes co-locating with their eQTL. We then performed BSA by assembling allele-specific pools based on the genotypes of individuals at the partial resistance QTL Rphq11. BSA identified a total of 411 genes as differentially expressed, including HvPHGPx, a gene previously identified as a promising candidate for Rphq11. The genetic location of 276 of these genes could be determined from both eQTL datasets and conservation of synteny, and 254 (92%) of these were located on the target chromosome. We conclude that the identification of differential expression by BSA constitutes a novel method to identify genes located in specific regions of interest. The datasets obtained from such studies provide a robust set of candidate genes for the analysis and serve as valuable resources for targeted marker development and comparative mapping with other grass species.

Candidate genes associated with bud dormancy release in blackcurrant (Ribes nigrum L.)

BACKGROUND

The detrimental effects of mild winter temperatures on the consistency of cropping of blackcurrant (Ribes nigrum L.) in parts of Europe have led to increasing interest in the genetic control of dormancy release in this species. This study examined patterns of gene expression in leaf buds of blackcurrant to identify key differential changes in these profiles around the time of budbreak.

RESULTS

Using leaf bud tissue of blackcurrant, a cDNA library was generated as a source of blackcurrant ESTs for construction of a custom microarray, which was used to identify differential gene expression during dormancy release. Gene activity was lowest in early stages of dormancy, increasing to reach a maximum around the time of budbreak. Genes with significantly changing expression profiles were clustered and evidence is provided for the transient activity of genes previously associated with dormancy processes in other species. Expression profiling identified candidate genes which were mapped onto a blackcurrant genetic linkage map containing budbreak-related QTL. Three genes, which putatively encode calmodulin-binding protein, beta tubulin and acetyl CoA carboxylase respectively, were found to co-localise with budbreak QTL.

CONCLUSIONS

This study provides insight into the genetic control of dormancy transition in blackcurrant, identifying key changes in gene expression around budbreak. Genetic mapping of ESTs enabled the identification of genes which co-localise with previously-characterised blackcurrant QTL, and it is concluded that these genes have probable roles in release of dormancy and can therefore provide a basis for the development of genetic markers for future breeding deployment.

Effect of salinity on water relations of wild barley plants differing in salt tolerance

BACKGROUND AND AIMS

Certain lines of wild barley (Hordeum spontaneum) are more tolerant of salinity than others. The physiological basis of this difference is examined in a comparative study of a saline-tolerant and saline-intolerant line that emphasizes plant water relations.

METHODOLOGY

Effects of salt-treatment (75 mM maximum) extending from a few hours to 3 weeks were quantified in 8-day-old seedlings of a saline-sensitive wild barley line ('T-1') and a less saline-sensitive line ('20-45'). Plants were grown in nutrient culture. Levels of mRNA of the HtPIP2;4 aquaporin (AQP) gene were determined together with a range of physiological responses including root hydraulic conductivity, osmotic potential of root xylem sap, transpiration, leaf relative water content, root water content, leaf water potential, leaf sap osmolality, leaf length, leaf area and chlorophyll content.

PRINCIPAL RESULTS

Salt treatment inhibited transpiration and hydraulic conductivity more in salt-tolerant '20-45' plants than in salt-sensitive 'T-1'. In '20-45', the effect was paralleled by a fast (within a few hours) and persistent (3 days) down-regulation of aquaporin. In salt-sensitive 'T-1' plants, aquaporin down-regulation was delayed for up to 24 h. Greater tolerance in '20-45' plants was characterized by less inhibition of leaf area, root fresh weight, leaf water content and chlorophyll concentration. Leaf water potentials were similar in both lines.

CONCLUSIONS

(i) Decline in hydraulic conductivity in salt-treated barley plants is important for stomatal closure, (ii) lowered transpiration rate is beneficial for salt tolerance, at least at the seedling stage and (iii) changes in AQP expression are implicated in the control of whole plant hydraulic conductivity and the regulation of shoot water relations.

An eQTL analysis of partial resistance to Puccinia hordei in barley

BACKGROUND

Genetic resistance to barley leaf rust caused by Puccinia hordei involves both R genes and quantitative trait loci. The R genes provide higher but less durable resistance than the quantitative trait loci. Consequently, exploring quantitative or partial resistance has become a favorable alternative for controlling disease. Four quantitative trait loci for partial resistance to leaf rust have been identified in the doubled haploid Steptoe (St)/Morex (Mx) mapping population. Further investigations are required to study the molecular mechanisms underpinning partial resistance and ultimately identify the causal genes.

METHODOLOGY/PRINCIPAL FINDINGS

We explored partial resistance to barley leaf rust using a genetical genomics approach. We recorded RNA transcript abundance corresponding to each probe on a 15K Agilent custom barley microarray in seedlings from St and Mx and 144 doubled haploid lines of the St/Mx population. A total of 1154 and 1037 genes were, respectively, identified as being P. hordei-responsive among the St and Mx and differentially expressed between P. hordei-infected St and Mx. Normalized ratios from 72 distant-pair hybridisations were used to map the genetic determinants of variation in transcript abundance by expression quantitative trait locus (eQTL) mapping generating 15685 eQTL from 9557 genes. Correlation analysis identified 128 genes that were correlated with resistance, of which 89 had eQTL co-locating with the phenotypic quantitative trait loci (pQTL). Transcript abundance in the parents and conservation of synteny with rice allowed us to prioritise six genes as candidates for Rphq11, the pQTL of largest effect, and highlight one, a phospholipid hydroperoxide glutathione peroxidase (HvPHGPx) for detailed analysis.

CONCLUSIONS/SIGNIFICANCE

The eQTL approach yielded information that led to the identification of strong candidate genes underlying pQTL for resistance to leaf rust in barley and on the general pathogen response pathway. The dataset will facilitate a systems appraisal of this host-pathogen interaction and, potentially, for other traits measured in this population.

Microarray comparative genomic hybridisation analysis incorporating genomic organisation, and application to enterobacterial plant pathogens

Microarray comparative genomic hybridisation (aCGH) provides an estimate of the relative abundance of genomic DNA (gDNA) taken from comparator and reference organisms by hybridisation to a microarray containing probes that represent sequences from the reference organism. The experimental method is used in a number of biological applications, including the detection of human chromosomal aberrations, and in comparative genomic analysis of bacterial strains, but optimisation of the analysis is desirable in each problem domain.We present a method for analysis of bacterial aCGH data that encodes spatial information from the reference genome in a hidden Markov model. This technique is the first such method to be validated in comparisons of sequenced bacteria that diverge at the strain and at the genus level: Pectobacterium atrosepticum SCRI1043 (Pba1043) and Dickeya dadantii 3937 (Dda3937); and Lactococcus lactis subsp. lactis IL1403 and L. lactis subsp. cremoris MG1363. In all cases our method is found to outperform common and widely used aCGH analysis methods that do not incorporate spatial information. This analysis is applied to comparisons between commercially important plant pathogenic soft-rotting enterobacteria (SRE) Pba1043, P. atrosepticum SCRI1039, P. carotovorum 193, and Dda3937.Our analysis indicates that it should not be assumed that hybridisation strength is a reliable proxy for sequence identity in aCGH experiments, and robustly extends the applicability of aCGH to bacterial comparisons at the genus level. Our results in the SRE further provide evidence for a dynamic, plastic 'accessory' genome, revealing major genomic islands encoding gene products that provide insight into, and may play a direct role in determining, variation amongst the SRE in terms of their environmental survival, host range and aetiology, such as phytotoxin synthesis, multidrug resistance, and nitrogen fixation.

Molecular regulation of somatic embryogenesis in potato: an auxin led perspective

Potato internodal segments (INS) treated with the auxin 2,4-dichlorophenoxyacetic acid can be induced to develop somatic embryos upon their transfer to an auxin-free medium, while the continuous presence of auxin in the medium suppresses the progression of embryogenically-induced somatic cells to embryos. We have employed these contrasting pathways, in combination with potato microarrays representing circa 10,000 genes, to profile global gene expression patterns during the progression of somatic embryogenesis in potato. The induction phase, characterised by the presence of auxin, was analysed by the direct comparison of RNA isolated from freshly excised (0 days) and embryogenically induced (14 days) INS explants. RNAs from embryo-forming (withdrawal of auxin after 14 days) and embryo-inhibitory (continuous presence of auxin) conditions, isolated over a range of time-points until the emergence of somatic embryos, were compared in a loop design to identify auxin responsive genes putatively involved in the process of somatic embryogenesis. A total of 402 transcripts were found to be showing significant differential expression patterns during somatic embryogenesis 'induction' phase, 524 during 'embryo-transition' phase, while 44 transcripts were common to both phases. Functional classification of these transcripts, using Gene Ontology vocabularies (molecular and biological), revealed that a significant proportion of transcripts were involved in processes which are more relevant to somatic embryogenesis such as apoptosis, development, reproduction, stress and signal transduction. This is the first study profiling global gene expression patterns during true somatic embryogenesis initiated from mature and completely differentiated explants and has enabled the description of stage-specific expression patterns of a large number of genes during potato somatic embryogenesis (PSE). The significance of the key identified genes during critical stages of somatic embryogenesis is discussed.

Quorum sensing coordinates brute force and stealth modes of infection in the plant pathogen Pectobacterium atrosepticum

Quorum sensing (QS) in vitro controls production of plant cell wall degrading enzymes (PCWDEs) and other virulence factors in the soft rotting enterobacterial plant pathogen Pectobacterium atrosepticum (Pba). Here, we demonstrate the genome-wide regulatory role of QS in vivo during the Pba-potato interaction, using a Pba-specific microarray. We show that 26% of the Pba genome exhibited differential transcription in a QS (expI-) mutant, compared to the wild-type, suggesting that QS may make a greater contribution to pathogenesis than previously thought. We identify novel components of the QS regulon, including the Type I and II secretion systems, which are involved in the secretion of PCWDEs; a novel Type VI secretion system (T6SS) and its predicted substrates Hcp and VgrG; more than 70 known or putative regulators, some of which have been demonstrated to control pathogenesis and, remarkably, the Type III secretion system and associated effector proteins, and coronafacoyl-amide conjugates, both of which play roles in the manipulation of plant defences. We show that the T6SS and a novel potential regulator, VirS, are required for full virulence in Pba, and propose a model placing QS at the apex of a regulatory hierarchy controlling the later stages of disease progression in Pba. Our findings indicate that QS is a master regulator of phytopathogenesis, controlling multiple other regulators that, in turn, co-ordinately regulate genes associated with manipulation of host defences in concert with the destructive arsenal of PCWDEs that manifest the soft rot disease phenotype.

DsbA plays a critical and multifaceted role in the production of secreted virulence factors by the phytopathogen Erwinia carotovora subsp. atroseptica

Erwinia carotovora subsp. atroseptica is an enterobacterial phytopathogen causing economically significant soft rot disease. Pathogenesis is mediated by multiple secreted virulence factors, many of which are secreted by the type II (Out) secretion system. DsbA catalyzes the introduction of disulfide bonds into periplasmic and secreted proteins. In this study, the extracellular proteome (secretome) of wild type E. carotovora subsp. atroseptica SCRI1043, and dsbA and out mutants, was analyzed by spectral counting mass spectrometry. This revealed that dsbA inactivation had a huge impact on the secretome and identified diverse DsbA- and Out-dependent secreted proteins, representing known, predicted, and novel candidate virulence factors. Further characterization of the dsbA mutant showed that secreted enzyme activities, motility, production of the quorum-sensing signal, and virulence were absent or substantially reduced. The impact of DsbA on secreted virulence factor production was mediated at multiple levels, including impacting on the Out secretion system and the virulence gene regulatory network. Transcriptome analyses revealed that the abundance of a broad, but defined, set of transcripts, including many virulence factors, was altered in the dsbA mutant, identifying a new virulence regulon responsive to extracytoplasmic conditions. In conclusion, DsbA plays a crucial, multifaceted role in the pathogenesis of E. carotovora subsp. atroseptica.

Efflux pump gene expression in Erwinia chrysanthemi is induced by exposure to phenolic acids

Salicylic acid (SA) is an important signaling molecule in local and systemic plant resistance. Following infection by microbial pathogens and the initial oxidative burst in plants, SA accumulation functions in the amplification of defense gene expression. Production of pathogenesis-related proteins and toxic antimicrobial chemicals serves to protect the plant from infection. Successful microbial pathogens utilize a variety of mechanisms to rid themselves of toxic antimicrobial compounds. Important among these mechanisms are multidrug-resistance pumps that bring about the active efflux of toxic compounds from microbial cells. Here, we show that a combination SA and its precursors, t-cinnamic acid and benzoic acid, can activate expression of specific multidrug efflux pump-encoding genes in the plant pathogen Erwinia chrysanthemi and enhance survival of the bacterium in the presence of model as well as plant-derived antimicrobial chemicals. This ability of plant-pathogenic bacteria to co-opt plant defense-signaling molecules to activate multidrug efflux pumps may have evolved to ensure bacterial survival in susceptible host plants.

Co-ordinated gene expression during phases of dormancy release in raspberry (Rubus idaeus L.) buds

Bud break in raspberry (Rubus idaeus L.) is often poor and uneven, with many of the subapical buds remaining in a dormant state. In order to determine the dormancy status of raspberry buds, an empirical measure of bud burst in a growth-permissive environment following exposure to chilling (4 degrees C cold storage) was developed. For cv. Glen Ample, percentage bud burst in intact canes and isolated nodes was recorded after 14 d. Isolated nodes (a measure of endodormancy) achieved 100% bud burst after approximately 1500 h chilling whereas buds on intact plants (combined endo- and paradormancy) required an additional 1000 h chilling. A microarray approach was used to follow changes in gene expression that occurred during dormancy transition. The probes for the microarrays were obtained from endodormant and paradormant raspberry bud cDNA libraries. The expression profiles of 5300 clones from these libraries were subjected to principal component analysis to determine the most significant expression patterns. Sequence analysis of these clones, in many cases, enabled their functional categorization and the development of hypotheses concerning the mechanisms of bud dormancy release. Thus a set of novel candidates for key dormancy-related genes from raspberry buds have been identified. Bud dormancy is fundamental to the study of plant developmental processes and, in addition, its regulation is of significant economic importance to fruit and horticultural industries.

Genome-wide SNP discovery and linkage analysis in barley based on genes responsive to abiotic stress

More than 2,000 genome-wide barley single nucleotide polymorphisms (SNPs) were developed by resequencing unigene fragments from eight diverse accessions. The average genome-wide SNP frequency observed in 877 unigenes was 1 SNP per 200 bp. However, SNP frequency was highly variable with the least number of SNP and SNP haplotypes observed within European cultivated germplasm reflecting effects of breeding history on genetic diversity. More than 300 SNP loci were mapped genetically in three experimental mapping populations which allowed the construction of an integrated SNP map incorporating a large number of RFLP, AFLP and SSR markers (1,237 loci in total). The genes used for SNP discovery were selected based on their transcriptional response to a variety of abiotic stresses. A set of known barley abiotic stress QTL was positioned on the linkage map, while the available sequence and gene expression information facilitated the identification of genes potentially associated with these traits. Comparison of the sequenced SNP loci to the rice genome sequence identified several regions of highly conserved gene order providing a framework for marker saturation in barley genomic regions of interest. The integration of genome-wide SNP and expression data with available genetic and phenotypic information will facilitate the identification of gene function in barley and other non-model organisms.

Genome-wide SNP discovery and linkage analysis in barley based on genes responsive to abiotic stress

More than 2,000 genome-wide barley single nucleotide polymorphisms (SNPs) were developed by resequencing unigene fragments from eight diverse accessions. The average genome-wide SNP frequency observed in 877 unigenes was 1 SNP per 200 bp. However, SNP frequency was highly variable with the least number of SNP and SNP haplotypes observed within European cultivated germplasm reflecting effects of breeding history on genetic diversity. More than 300 SNP loci were mapped genetically in three experimental mapping populations which allowed the construction of an integrated SNP map incorporating a large number of RFLP, AFLP and SSR markers (1,237 loci in total). The genes used for SNP discovery were selected based on their transcriptional response to a variety of abiotic stresses. A set of known barley abiotic stress QTL was positioned on the linkage map, while the available sequence and gene expression information facilitated the identification of genes potentially associated with these traits. Comparison of the sequenced SNP loci to the rice genome sequence identified several regions of highly conserved gene order providing a framework for marker saturation in barley genomic regions of interest. The integration of genome-wide SNP and expression data with available genetic and phenotypic information will facilitate the identification of gene function in barley and other non-model organisms.

Single-feature polymorphism discovery in the barley transcriptome

A probe level model for analysis of GeneChip gene expression data is presented which identified more than 10,000 single-feature polymorphisms between two barley genotypes, with a high sensitivity. This method is applicable to all oligonucleotide microarray data.

A probe-level model for analysis of GeneChip gene-expression data is presented which identified more than 10,000 single-feature polymorphisms (SFP) between two barley genotypes. The method has good sensitivity, as 67% of known single-nucleotide polymorphisms (SNP) were called as SFPs. This method is applicable to all oligonucleotide microarray data, accounts for SNP effects in gene-expression data and represents an efficient and versatile approach for highly parallel marker identification in large genomes.

A comparative analysis of transcript abundance using SAGE and Affymetrix arrays

A number of methods are currently used for gene expression profiling. They differ in scale, economy and sensitivity. We present the results of a direct comparison between serial analysis of gene expression (SAGE) and the Barley1 Affymetrix GeneChip. Both technology platforms were used to obtain quantitative measurements of transcript abundance using identical RNA samples and assessed for their ability to quantify differential gene expression. For SAGE, a total of 82,122 tags were generated from two independent libraries representing whole developing barley caryopsis and dissected embryos. The Barley1 GeneChip contains 22,791 probe sets. Results obtained from both methods are generally comparable, indicating that both will lead to similar conclusions regarding transcript levels and differential gene expression. However, excluding singletons, 24.4% of the unique SAGE tags had no corresponding probe set on the Barley1 array indicating that a broader snapshot of gene expression was obtained by SAGE. Discrepancies were observed for a number of "genes" and these are discussed.

Metabolic engineering of high carotenoid potato tubers containing enhanced levels of beta-carotene and lutein

In order to enhance the carotenoid content of potato tubers, transgenic potato plants have been produced expressing an Erwinia uredovora crtB gene encoding phytoene synthase, specifically in the tuber of Solanum tuberosum L. cultivar Desiree which normally produces tubers containing c. 5.6 microg carotenoid g(-1) DW and also in Solanum phureja L. cv. Mayan Gold which has a tuber carotenoid content of typically 20 microg carotenoid g(-1) DW. In developing tubers of transgenic crtB Desiree lines, carotenoid levels reached 35 microg carotenoid g(-1) DW and the balance of carotenoids changed radically compared with controls: beta-carotene levels in the transgenic tubers reached c. 11 microg g(-1) DW, whereas control tubers contained negligible amounts and lutein accumulated to a level 19-fold higher than empty-vector transformed controls. The crtB gene was also transformed into S. phureja (cv. Mayan Gold), again resulting in an increase in total carotenoid content to 78 microg carotenoid g(-1) DW in the most affected transgenic line. In these tubers, the major carotenoids were violaxanthin, lutein, antheraxanthin, and beta-carotene. No increases in expression levels of the major carotenoid biosynthetic genes could be detected in the transgenic tubers, despite the large increase in carotenoid accumulation. Microarray analysis was used to identify a number of genes that were consistently up- or down-regulated in transgenic crtB tubers compared with empty vector controls. The implications of these data from a nutritional standpoint and for further modifications of tuber carotenoid content are discussed.

Characterisation of early transcriptional changes involving multiple signalling pathways in the Mla13 barley interaction with powdery mildew ( Blumeria graminis f. sp. hordei)

Suppression subtractive hybridisation was used to isolate 21 cDNAs ( bmi1- bmi21) up-regulated 1-5 h post-inoculation (hpi) in a barley ( Hordeum vulgare L. cv. Pallas) near-isogenic line (NIL) P11 ( Mla13) challenged with either avirulent or virulent isolates of Blumeria graminis f. sp. hordei. Transcriptional changes at these time-points are crucial for the Mla-mediated hypersensitive response [W.R. Bushnell and Z. Liu (1994) Physiol Mol Plant Pathol 44:389-402]. Seven sequences were up-regulated by 1 hpi, when the pathogen has formed only the primary germ tube. Some transcripts were similar to genes with a role in regulating programmed cell death in animals, including NF kappaB and oxysterol-binding protein. Moreover, bmi7, similar to rice resistance gene Xa21, was rapidly up-regulated in both compatible and incompatible interactions, but was then down-regulated by 5 hpi in the virulent interaction. Only nine of the transcripts were up-regulated in mlo5 resistance in cv. Pallas NIL P22, confirming differential pathway induction between Mla13 and mlo5. However, eight sequences up-regulated in the Mla13 response in P11 were already highly elevated in uninoculated mlo5 mutant P22, suggesting that they may be negatively regulated by wild-type Mlo. Regulation of bmi sequences was investigated using salicylic acid, methyl jasmonate, ethylene, H(2)O(2), abscisic acid, wounding and a glucan elicitor. No single stimulus up-regulated all genes, suggesting either combinations of these stimuli, or additional stimuli, are involved in early Mla13 and mlo5 resistances. Whereas H(2)O(2) up- or down-regulated 17 of the transcripts detected in Northern analyses, salicylic acid stimulated only down-regulation of 5 transcripts.

Differential expression of invertase genes in internal and external phloem tissues of potato (Solanum tuberosum L.)

The cloning of promoter sequences of two invertase genes from potato (Solanum tuberosum L.) is described. Histochemical analysis of series of reporter transgenic lines reveals phloem-specific expression from both promoters, with one expressed preferentially in internal phloem and the other in external phloem of stem vascular bundles.

Requirements for mini-exon inclusion in potato invertase mRNAs provides evidence for exon-scanning interactions in plants

Invertases are responsible for the breakdown of sucrose to fructose and glucose. In all but one plant invertase gene, the second exon is only 9 nt in length and encodes three amino acids of a five-amino-acid sequence that is highly conserved in all invertases of plant origin. Sequences responsible for normal splicing (inclusion) of exon 2 have been investigated in vivo using the potato invertase, invGF gene. The upstream intron 1 is required for inclusion whereas the downstream intron 2 is not. Mutations within intron 1 have identified two sequence elements that are needed for inclusion: a putative branchpoint sequence and an adjacent U-rich region. Both are recognized plant intron splicing signals. The branchpoint sequence lies further upstream from the 3' splice site of intron 1 than is normally seen in plant introns. All dicotyledonous plant invertase genes contain this arrangement of sequence elements: a distal branchpoint sequence and adjacent, downstream U-rich region. Intron 1 sequences upstream of the branchpoint and sequences in exons 1, 2, or 3 do not determine inclusion, suggesting that intron or exon splicing enhancer elements seen in vertebrate mini-exon systems are absent. In addition, mutation of the 3' and 5' splice sites flanking the mini-exon cause skipping of the mini-exon, suggesting that both splice sites are required. The branchpoint/U-rich sequence is able to promote splicing of mini-exons of 6, 3, and 1 nt in length and of a chicken cTNT mini-exon of 6 nt. These sequence elements therefore act as a splicing enhancer and appear to function via interactions between factors bound at the branchpoint/U-rich region and at the 5' splice site of intron 2, activating removal of this intron followed by removal of intron 1. This first example of splicing of a plant mini-exon to be analyzed demonstrates that particular arrangement of standard plant intron splicing signals can drive constitutive splicing of a mini-exon.

Microsatellite repeats are not randomly distributed within Norway spruce (Picea abies K.) expressed sequences

A Norway spruce (Picea abies K.) cDNA library obtained from vegetative bud tissue was screened for the presence of (AG)n and (AC)n microsatellite repeats. Ten (AG)n and six (AC)n microsatellites were found, with an average length of 25.5 repeat units. Most of the microsatellites are simple perfect repeats. The microsatellite distribution within the clones is clearly non-random, with different classes of repeats lying in different positions relative to the coding region and in a highly conserved orientation. An estimate of the frequency of dinucleotide microsatellites in expressed regions was obtained, showing that SSRs (simple sequence repeats) are found in genes about 20 times less frequently than in random genomic clones, with (AG)n repeats more frequent than (AC)n repeats. Potential applications of these sequences as expressed region-based molecular markers are shown by developing six SSR markers for the detection of natural variation in Norway spruce populations and testing two of them for the identification of illegitimate progenies from a mapping population. Key words: Picea abies, microsatellites, SSRs, ESTs, population genetics, trees.

Microsatellite repeats are not randomly distributed within Norway spruce (Picea abies K.) expressed sequences

A Norway spruce (Picea abies K.) cDNA library obtained from vegetative bud tissue was screened for the presence of (AG)n and (AC)n microsatellite repeats. Ten (AG)n and six (AC)n microsatellites were found, with an average length of 25.5 repeat units. Most of the microsatellites are simple perfect repeats. The microsatellite distribution within the clones is clearly non-random, with different classes of repeats lying in different positions relative to the coding region and in a highly conserved orientation. An estimate of the frequency of dinucleotide microsatellites in expressed regions was obtained, showing that SSRs (simple sequence repeats) are found in genes about 20 times less frequently than in random genomic clones, with (AG)n repeats more frequent than (AC)n repeats. Potential applications of these sequences as expressed region-based molecular markers are shown by developing six SSR markers for the detection of natural variation in Norway spruce populations and testing two of them for the identification of illegitimate progenies from a mapping population.

Expression of tandem invertase genes associated with sexual and vegetative growth cycles in potato

The organisation of two invertase genes (invGE and invGF) linked in direct tandem repeat within the potato genome is detailed. The genes exhibit a similar intron/exon structure which differs from previously described plant invertase genes; while intron locations are conserved between the genes, minor differences in exon length are seen. Both genes encode enzymes with putative extracellular location. Biochemical analysis of gene expression showed expression in floral tissues for both genes, with expression of the upstream gene (invGE) also detected in leaf tissue. Promoter sequences from both genes have been fused to the beta-glucuronidase (GUS) reporter gene (uidA) and transformed into potato. One promoter-GUS reporter construct was also transformed into tobacco. Histochemical analysis of transgenic lines defined specific expression from the downstream (invGF) promoter in potato and tobacco pollen, with expression first detected in the late uninucleate stage of tobacco microspore development. The invGE promoter determined expression in pollen and other floral tissues, but also at lateral nodes in stem, root and tuber. An association of invertase expression with generative tissue, both in vegetative and sexual modes of growth, is indicated.

Exon skipping induced by cold stress in a potato invertase gene transcript

We show that two invertase genes in potato, like most other plant invertase genes, include a very short second exon of 9 bp which encodes the central three amino acids of a motif highly conserved in invertases of diverse origin. This mini-exon is one of the smallest known in plants and pre-mRNA from these genes may be susceptible to alternative splicing, because of a potential requirement for specialized interaction with the splicing machinery to ensure correct processing for the production of a mature mRNA. No evidence of aberrant post-transcriptional processing was observed during normal invertase gene expression in potato. The fidelity of post-transcriptional processing of the pre-mRNA from one of the genes was perturbed by cold stress, resulting in the deletion of the mini-exon from some transcripts. This alternative splicing event occurred under cold stress in both leaf and stem, but was not induced by wounding. This adds an example of exon skipping and the induction of alternative processing by cold stress to the small number of transcripts which have been shown to exhibit alternative splicing in plants. The differential sensitivity of post-transcriptional processing to cold stress observed for the two transcripts examined will permit further dissection of the nucleotide sequence requirements for their accurate splicing.

Characterisation of a complementary DNA encoding a novel plant enzyme with sucrolytic activity

The cloning of a 1332 bp cDNA from a potato (Solanum tuberosum L.) cv. Cara leaf cDNA expression library, using an antibody raised against a purified tuber protein preparation with sucrolytic activity, is described. The corresponding gene in potato is of low copy number, is expressed in a variety of tissues, and encodes a protein which includes several domains with similarity to database sequences, including ferredoxin from Clostridium pasteurianum. Expression of the cDNA in E. coli yields a fusion protein with sucrolytic activity.

Potato (Solanum tuberosum) invertase-encoding cDNAs and their differential expression

A full-length cDNA clone encoding a potato invertase (Inv) has been isolated. It is highly related (77% nucleotide identity) to a previously characterised potato cDNA clone encoding a putative extracellular Inv. These Inv genes encode a subfamily of apoplastic enzymes which are shown to be distinct, on the basis of sequence similarity, from the related subfamily of vacuolar enzymes. In order to differentiate between the expression of the two potato genes encoding apoplastic Inv, a single-stranded conformational polymorphism (SSCP) assay was developed for products generated by reverse transcription-polymerase chain reaction (RT-PCR) utilising primers designed to amplify both potato sequences. Using this approach, we have shown that these two identified Inv from potato are expressed in a tissue-specific and developmentally regulated manner.

cDNA cloning and expression of a potato (Solanum tuberosum) invertase

A cDNA clone encoding an invertase isoenzyme has been isolated from a potato leaf cDNA library. The deduced amino acid sequence shows significant similarities to previously characterised invertases. The highest degree of overall similarity, including the signal peptide sequence, is to carrot cell wall invertase, suggesting that the potato gene encodes an apoplastic enzyme. Expression of the gene, as determined by RT-PCR, is detected in stem and leaf tissue, and at lower levels in tuber, but is absent from roots.

Developmental regulation of insulin like growth factor II gene expression in the pig

The expression of the insulin like growth factor II gene has been examined in the developing porcine embryo. It was found that IGF II transcripts were present in abundant quantities in first trimester embryos as well as in the extraembryonic tissues amnion and allantochorion. The expression of the IGF II gene was high in the fetal liver where a prominent 2.3 kB transcript and a less abundant 5.4 kB transcript were produced. However, neither of these transcripts could be traced in the adult liver. Instead we found two different IGF II transcripts with the sizes of 4.7 and 1.2 kB in the adult liver. These findings indicate that the IGF II gene is under developmental control with the possible existence of different promoters.