Rico A. Caldo

A New Resource for Cereal Genomics: 22K Barley GeneChip Comes of Age

In recent years, access to complete genomic sequences, coupled with rapidly accumulating data related to RNA and protein expression patterns, has made it possible to determine comprehensively how genes contribute to complex phenotypes. However, for major crop plants, publicly available, standard platforms for parallel expression analysis have been limited. We report the conception and design of the new publicly available, 22K Barley1 GeneChip probe array, a model for plants without a fully sequenced genome. Array content was derived from worldwide contribution of 350,000 high-quality ESTs from 84 cDNA libraries, in addition to 1,145 barley (Hordeum vulgare) gene sequences from the National Center for Biotechnology Information nonredundant database. Conserved sequences expressed in seedlings of wheat (Triticum aestivum), oat (Avena strigosa), rice (Oryza sativa), sorghum (Sorghum bicolor), and maize (Zea mays) were identified that will be valuable in the design of arrays across grasses. To enhance the usability of the data, BarleyBase, a MIAME-compliant, MySQL relational database, serves as a public repository for raw and normalized expression data from the Barley1 GeneChip probe array. Interconnecting links with PlantGDB and Gramene allow BarleyBase users to perform gene predictions using the 21,439 non-redundant Barley1 exemplar sequences or cross-species comparison at the genome level, respectively. We expect that this first generation array will accelerate hypothesis generation and gene discovery in disease defense pathways, responses to abiotic stresses, development, and evolutionary diversity in monocot plants.

Interaction-Dependent Gene Expression in Mla-Specified Response to Barley Powdery Mildew

Plant recognition of pathogen-derived molecules influences attack and counterattack strategies that affect the outcome of host–microbe interactions. To ascertain the global framework of host gene expression during biotrophic pathogen invasion, we analyzed in parallel the mRNA abundance of 22,792 host genes throughout 36 (genotype × pathogen × time) interactions between barley (Hordeum vulgare) and Blumeria graminis f. sp hordei (Bgh), the causal agent of powdery mildew disease. A split-split-plot design was used to investigate near-isogenic barley lines with introgressed Mla6, Mla13, and Mla1 coiled-coil, nucleotide binding site, Leu-rich repeat resistance alleles challenged with Bgh isolates 5874 (AvrMla6 and AvrMla1) and K1 (AvrMla13 and AvrMla1). A linear mixed model analysis was employed to identify genes with significant differential expression (P value < 0.0001) in incompatible and compatible barley-Bgh interactions across six time points after pathogen challenge. Twenty-two host genes, of which five were of unknown function, exhibited highly similar patterns of upregulation among all incompatible and compatible interactions up to 16 h after inoculation (hai), coinciding with germination of Bgh conidiospores and formation of appressoria. By contrast, significant divergent expression was observed from 16 to 32 hai, during membrane-to-membrane contact between fungal haustoria and host epidermal cells, with notable suppression of most transcripts identified as differentially expressed in compatible interactions. These findings provide a link between the recognition of general and specific pathogen-associated molecules in gene-for-gene specified resistance and support the hypothesis that host-specific resistance evolved from the recognition and prevention of the pathogens suppression of plant basal defense.

An atlas of gene expression from seed to seed through barley development.

Assaying relative and absolute levels of gene expression in a diverse series of tissues is a central step in the process of characterizing gene function and a necessary component of almost all publications describing individual genes or gene family members. However, throughout the literature, such studies lack consistency in genotype, tissues analyzed, and growth conditions applied, and, as a result, the body of information that is currently assembled is fragmented and difficult to compare between different studies. The development of a comprehensive platform for assaying gene expression that is available to the entire research community provides a major opportunity to assess whole biological systems in a single experiment. It also integrates detailed knowledge and information on individual genes into a unified framework that provides both context and resource to explore their contributions in a broader biological system. We have established a data set that describes the expression of 21,439 barley genes in 15 tissues sampled throughout the development of the barley cv. Morex grown under highly controlled conditions. Rather than attempting to address a specific biological question, our experiment was designed to provide a reference gene expression data set for barley researchers; a gene expression atlas and a comparative data set for those investigating genes or regulatory networks in other plant species. In this paper we describe the tissues sampled and their transcriptomes, and provide summary information on genes that are either specifically expressed in certain tissues or show correlated expression patterns across all 15 tissue samples. Using specific examples and an online tutorial, we describe how the data set can be interrogated for patterns and levels of barley gene expression and how the resulting information can be used to generate and/or test specific biological hypotheses.

Estimating the number of true null hypotheses from a histogram of p values

In an earlier article, an intuitively appealing method for estimating the number of true null hypotheses in a multiple test situation was proposed. That article presented an iterative algorithm that relies on a histogram of observed p values to obtain the estimator. We characterize the limit of that iterative algorithm and show that the estimator can be computed directly without iteration. We compare the performance of the histogram-based estimator with other procedures for estimating the number of true null hypotheses from a collection of observed p values and find that the histogram-based estimator performs well in settings similar to those encountered in microarray data analysis. We demonstrate the approach using p values from a large microarray experiment aimed at uncovering molecular mechanisms of barley resistance to a fungal pathogen.

BarleyBase—an expression profiling database for plant genomics

BarleyBase (BB) (www.barleybase.org) is an online database for plant microarrays with integrated tools for data visualization and statistical analysis. BB houses raw and normalized expression data from the two publicly available Affymetrix genome arrays, Barley1 and Arabidopsis ATH1 with plans to include the new Affymetrix 61K wheat, maize, soybean and rice arrays, as they become available. BB contains a broad set of query and display options at all data levels, ranging from experiments to individual hybridizations to probe sets down to individual probes. Users can perform cross-experiment queries on probe sets based on observed expression profiles and/or based on known biological information. Probe set queries are integrated with visualization and analysis tools such as the R statistical toolbox, data filters and a large variety of plot types. Controlled vocabularies for gene and plant ontologies, as well as interconnecting links to physical or genetic map and other genomic data in PlantGDB, Gramene and GrainGenes, allow users to perform EST alignments and gene function prediction using Barley1 exemplar sequences, thus, enhancing cross-species comparison.

Stage-Specific Suppression of Basal Defense Discriminates Barley Plants Containing Fast- and Delayed-Acting Mla Powdery Mildew Resistance Alleles

Nonspecific recognition of pathogen-derived general elicitors triggers the first line of plant basal defense, which in turn, preconditions the host towards resistance or susceptibility. To elucidate how basal defense responses influence the onset of Mla (mildew resistance locus a)-specified resistance, we performed a meta-analysis of GeneChip mRNA expression for 155 basal defense-related genes of barley (Hordeum vulgare) challenged with Blumeria graminis f. sp. hordei, the causal agent of powdery mildew disease. In plants containing the fast-acting Mla1, Mla6, or Mla13 alleles, transcripts hyper-accumulated from 0 to 16 h after inoculation (hai) in both compatible and incompatible interactions. Suppression of basal defense-related transcripts was observed after 16 hai only in compatible interactions, whereas these transcripts were sustained or increased in incompatible interactions. By contrast, in plants containing wild-type and mutants of the delayed-acting Mla12 allele, an early hyper-induction of transcripts from 0 to 8 hai was observed, but the expression of many of these genes is markedly suppressed from 8 to 16 hai. These results suggest that the inhibition of basal defense facilitates the development of haustoria by the pathogen, consequently delaying the onset of host resistance responses. Thus, we hypothesize that the regulation of basal defense influences host-cell accessibility to the fungal pathogen and drives allelic diversification of gene-specific resistance phenotypes.

Code-Assisted Discovery of TAL Effector Targets in Bacterial Leaf Streak of Rice Reveals Contrast with Bacterial Blight and a Novel Susceptibility Gene

Bacterial leaf streak of rice, caused by Xanthomonas oryzae pv. oryzicola (Xoc) is an increasingly important yield constraint in this staple crop. A mesophyll colonizer, Xoc differs from X. oryzae pv. oryzae (Xoo), which invades xylem to cause bacterial blight of rice. Both produce multiple distinct TAL effectors, type III-delivered proteins that transactivate effector-specific host genes. A TAL effector finds its target(s) via a partially degenerate code whereby the modular effector amino acid sequence identifies nucleotide sequences to which the protein binds. Virulence contributions of some Xoo TAL effectors have been shown, and their relevant targets, susceptibility (S) genes, identified, but the role of TAL effectors in leaf streak is uncharacterized. We used host transcript profiling to compare leaf streak to blight and to probe functions of Xoc TAL effectors. We found that Xoc and Xoo induce almost completely different host transcriptional changes. Roughly one in three genes upregulated by the pathogens is preceded by a candidate TAL effector binding element. Experimental analysis of the 44 such genes predicted to be Xoc TAL effector targets verified nearly half, and identified most others as false predictions. None of the Xoc targets is a known bacterial blight S gene. Mutational analysis revealed that Tal2g, which activates two genes, contributes to lesion expansion and bacterial exudation. Use of designer TAL effectors discriminated a sulfate transporter gene as the S gene. Across all targets, basal expression tended to be higher than genome-average, and induction moderate. Finally, machine learning applied to real vs. falsely predicted targets yielded a classifier that recalled 92% of the real targets with 88% precision, providing a tool for better target prediction in the future. Our study expands the number of known TAL effector targets, identifies a new class of S gene, and improves our ability to predict functional targeting.

Gene Expression Biomarkers Provide Sensitive Indicators of in Planta Nitrogen Status in Maize

Over the last several decades, increased agricultural production has been driven by improved agronomic practices and a dramatic increase in the use of nitrogen-containing fertilizers to maximize the yield potential of crops. To reduce input costs and to minimize the potential environmental impacts of nitrogen fertilizer that has been used to optimize yield, an increased understanding of the molecular responses to nitrogen under field conditions is critical for our ability to further improve agricultural sustainability. Using maize (Zea mays) as a model, we have characterized the transcriptional response of plants grown under limiting and sufficient nitrogen conditions and during the recovery of nitrogen-starved plants. We show that a large percentage (approximately 7%) of the maize transcriptome is nitrogen responsive, similar to previous observations in other plant species. Furthermore, we have used statistical approaches to identify a small set of genes whose expression profiles can quantitatively assess the response of plants to varying nitrogen conditions. Using a composite gene expression scoring system, this single set of biomarker genes can accurately assess nitrogen responses independently of genotype, developmental stage, tissue type, or environment, including in plants grown under controlled environments or in the field. Importantly, the biomarker composite expression response is much more rapid and quantitative than phenotypic observations. Consequently, we have successfully used these biomarkers to monitor nitrogen status in real-time assays of field-grown maize plants under typical production conditions. Our results suggest that biomarkers have the potential to be used as agronomic tools to monitor and optimize nitrogen fertilizer usage to help achieve maximal crop yields.

Quantitative and temporal definition of the Mla transcriptional regulon during barley-powdery mildew interactions.

Barley Mildew resistance locus a (Mla) is a major determinant of immunity to the powdery mildew pathogen, Blumeria graminis f. sp. hordei. Alleles of Mla encode cytoplasmic- and membrane-localized coiled-coil, nucleotide binding site, leucine-rich repeat proteins that mediate resistance when complementary avirulence effectors (AVR(a)) are present in the pathogen. Presence of an appropriate AVR(a) protein triggers nuclear relocalization of MLA, in which MLA binds repressing host transcription factors. Timecourse expression profiles of plants harboring Mla1, Mla6, and Mla12 wild-type alleles versus paired loss-of-function mutants were compared to discover conserved transcriptional targets of MLA and downstream signaling cascades. Pathogen-dependent gene expression was equivalent or stronger in susceptible plants at 20 h after inoculation (HAI) and was attenuated at later timepoints, whereas resistant plants exhibited a time-dependent strengthening of the transcriptional response, increasing in both fold change and the number of genes differentially expressed. Deregulation at 20 HAI implicated 16 HAI as a crucial point in determining the future trajectory of this interaction and was interrogated by quantitative analysis. In total, 28 potential transcriptional targets of the MLA regulon were identified. These candidate targets possess a diverse set of predicted functions, suggesting that multiple pathways are required to mediate the hypersensitive reaction.

Transcript-Based Cloning of RRP46, a Regulator of rRNA Processing and R Gene–Independent Cell Death in Barley–Powdery Mildew Interactions

Programmed cell death (PCD) plays a pivotal role in plant development and defense. To investigate the interaction between PCD and R gene–mediated defense, we used the 22K Barley1 GeneChip to compare and contrast time-course expression profiles of Blumeria graminis f. sp hordei (Bgh) challenged barley (Hordeum vulgare) cultivar C.I. 16151 (harboring the Mla6 powdery mildew resistance allele) and its fast neutron–derived Bgh-induced tip cell death1 mutant, bcd1. Mixed linear model analysis identified genes associated with the cell death phenotype as opposed to R gene–mediated resistance. One-hundred fifty genes were found at the threshold P value < 0.0001 and a false discovery rate <0.6%. Of these, 124 were constitutively overexpressed in the bcd1 mutant. Gene Ontology and rice (Oryza sativa) alignment-based annotation indicated that 68 of the 124 overexpressed genes encode ribosomal proteins. A deletion harboring six genes on chromosome 5H cosegregates with bcd1-specified cell death and is associated with misprocessing of rRNAs but segregates independent of R gene–mediated resistance. Barley stripe mosaic virus-induced gene silencing of one of the six deleted genes, RRP46 (rRNA-processing protein 46), phenocopied bcd1-mediated tip cell death. These findings suggest that RRP46, a critical component of the exosome core, mediates RNA processing and degradation involved in cell death initiation as a result of attempted penetration by Bgh during the barley–powdery mildew interaction but is independent of gene-for-gene resistance.