Gary J. Muehlbauer

A physical, genetic and functional sequence assembly of the barley genome

Barley (Hordeum vulgare L.) is among the world’s earliest domesticated and most important crop plants. It is diploid with a large haploid genome of 5.1 gigabases (Gb). Here we present an integrated and ordered physical, genetic and functional sequence resource that describes the barley gene-space in a structured whole-genome context. We developed a physical map of 4.98 Gb, with more than 3.90 Gb anchored to a high-resolution genetic map. Projecting a deep whole-genome shotgun assembly, complementary DNA and deep RNA sequence data onto this framework supports 79,379 transcript clusters, including 26,159 ‘high-confidence’ genes with homology support from other plant genomes. Abundant alternative splicing, premature termination codons and novel transcriptionally active regions suggest that post-transcriptional processing forms an important regulatory layer. Survey sequences from diverse accessions reveal a landscape of extensive single-nucleotide variation. Our data provide a platform for both genome-assisted research and enabling contemporary crop improvement.

The Fusarium graminearum Genome Reveals a Link Between Localized Polymorphism and Pathogen Specialization

We sequenced and annotated the genome of the filamentous fungus Fusarium graminearum, a major pathogen of cultivated cereals. Very few repetitive sequences were detected, and the process of repeat-induced point mutation, in which duplicated sequences are subject to extensive mutation, may partially account for the reduced repeat content and apparent low number of paralogous (ancestrally duplicated) genes. A second strain of F. graminearum contained more than 10,000 single-nucleotide polymorphisms, which were frequently located near telomeres and within other discrete chromosomal segments. Many highly polymorphic regions contained sets of genes implicated in plant-fungus interactions and were unusually divergent, with higher rates of recombination. These regions of genome innovation may result from selection due to interactions of F. graminearum with its plant hosts.

RNA-Seq Atlas of Glycine max : A guide to the soybean transcriptome

BackgroundNext generation sequencing is transforming our understanding of transcriptomes. It can determine the expression level of transcripts with a dynamic range of over six orders of magnitude from multiple tissues, developmental stages or conditions. Patterns of gene expression provide insight into functions of genes with unknown annotation.ResultsThe RNA Seq-Atlas presented here provides a record of high-resolution gene expression in a set of fourteen diverse tissues. Hierarchical clustering of transcriptional profiles for these tissues suggests three clades with similar profiles: aerial, underground and seed tissues. We also investigate the relationship between gene structure and gene expression and find a correlation between gene length and expression. Additionally, we find dramatic tissue-specific gene expression of both the most highly-expressed genes and the genes specific to legumes in seed development and nodule tissues. Analysis of the gene expression profiles of over 2,000 genes with preferential gene expression in seed suggests there are more than 177 genes with functional roles that are involved in the economically important seed filling process. Finally, the Seq-atlas also provides a means of evaluating existing gene model annotations for the Glycine max genome.ConclusionsThis RNA-Seq atlas extends the analyses of previous gene expression atlases performed using Affymetrix GeneChip technology and provides an example of new methods to accommodate the increase in transcriptome data obtained from next generation sequencing. Data contained within this RNA-Seq atlas of Glycine max can be explored at http://www.soybase.org/soyseq.

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.

Fungal Development and Induction of Defense Response Genes During Early Infection of Wheat Spikes by Fusarium graminearum

Fusarium head blight (FHB) of wheat is a crippling disease that causes severe economic losses in many of the wheat-growing regions of the world. Temporal patterns of fungus development and transcript accumulation of defense response genes were studied in Fusarium graminearum-inoculated wheat spikes within the first 48 to 76 h after inoculation (hai). Microscopy of inoculated glumes revealed that the fungus appeared to penetrate through stomata, exhibited subcuticular growth along stomatal rows, colonized glume parenchyma cells, and sporulated within 48 to 76 hai. No major differences in the timing of these events were found between Sumai 3 (resistant) and Wheaton (susceptible) genotypes. In complementary experiments, RNA was extracted from spikes at several time intervals up to 48 hai and temporal expression patterns were determined for defense response genes encoding peroxidase, PR-1, PR-2 (beta-1,3-glucanase), PR-3 (chitinase), PR-4, and PR-5 (thaumatin-like protein). In both genotypes, transcripts for the six defense response genes accumulated as early as 6 to 12 hai during F. graminearum infection and peaked at 36 to 48 hai. Greater and earlier PR-4 and PR-5 transcript accumulation was observed in Sumai 3, compared with Wheaton. Our results show that the timing of defense response gene induction is correlated with F. graminearum infection.

A bacterial artificial chromosome library for barley (Hordeum vulgare L.) and the identification of clones containing putative resistance genes

Abstract Modern cultivated barley is an important cereal crop with an estimated genome size of 5000 Mb. To develop the resources for positional cloning and structural genomic analyses in barley, we constructed a bacterial artificial chromosome (BAC) library for the cultivar Morex using the cloning enzyme HindIII. The library contains 313344 clones (816 384-well plates). A random sampling of 504 clones indicated an average insert size of 106 kbp (range=30–195 kbp) and 3.4% empty vectors. Screening the colony filters for chloroplast DNA content indicated an exceptionally low 1.5% contamination with chloroplast DNA. Thus, the library provides 6.3 haploid genome equivalents allowing a >99% probability of recovering any specific sequence of interest. High-density filters were gridded robotically using a Genetix Q-BOT in a 4×4 double-spotted array on 22.5-cm2 filters. Each set of 17 filters allows the entire library to be screened with 18432 clones represented per filter. Screening the library with 40 single copy probes identified an average 6.4 clones per probe, with a range of 1–13 clones per probe. A set of resistance-gene analog (RGA) sequences identified 121 RGA-containing BAC clones representing 20 different regions of the genome with an average of 6.1 clones per locus. Additional screening of the library with a P-loop disease resistance primer probe identified 459 positive BAC clones. These data indicate that this library is a valuable resource for structural genomic applications in barley.

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.

INTERMEDIUM-C, a modifier of lateral spikelet fertility in barley, is an ortholog of the maize domestication gene TEOSINTE BRANCHED 1

The domestication of cereals has involved common changes in morphological features, such as seed size, seed retention and modification of vegetative and inflorescence architecture that ultimately contributed to an increase in harvested yield. In barley, this process has resulted in two different cultivated types, two-rowed and six-rowed forms, both derived from the wild two-rowed ancestor, with archaeo-botanical evidence indicating the origin of six-rowed barley early in the domestication of the species, some 8,600–8,000 years ago. Variation at SIX-ROWED SPIKE 1 (VRS1) is sufficient to control this phenotype. However, phenotypes imposed by VRS1 alleles are modified by alleles at the INTERMEDIUM-C (INT-C) locus. Here we show that INT-C is an ortholog of the maize domestication gene TEOSINTE BRANCHED 1 (TB1) and identify 17 coding mutations in barley TB1 correlated with lateral spikelet fertility phenotypes.

Transcriptome analysis of the barley-Fusarium graminearum interaction.

Fusarium head blight (FHB) of barley (Hordeum vulgare L.) is caused by Fusarium graminearum. FHB causes yield losses and reduction in grain quality primarily due to the accumulation of trichothecene mycotoxins such as deoxynivalenol (DON). To develop an understanding of the barley-F. graminearum interaction, we examined the relationship among the infection process, DON concentration, and host transcript accumulation for 22,439 genes in spikes from the susceptible cv. Morex from 0 to 144 h after F. graminearum and water control inoculation. We detected 467 differentially accumulating barley gene transcripts in the F. graminearum-treated plants compared with the water control-treated plants. Functional annotation of the transcripts revealed a variety of infection-induced host genes encoding defense response proteins, oxidative burst-associated enzymes, and phenylpropanoid pathway enzymes. Of particular interest was the induction of transcripts encoding potential trichothecene catabolic enzymes and transporters, and the induction of the tryptophan biosynthetic and catabolic pathway enzymes. Our results define three stages of E graminearum infection. An early stage, between 0 and 48 h after inoculation (hai), exhibited limited fungal development, low DON accumulation, and little change in the transcript accumulation status. An intermediate stage, between 48 and 96 hai, showed increased fungal development and active infection, higher DON accumulation, and increased transcript accumulation. A majority of the host gene transcripts were detected by 72 hai, suggesting that this is an important timepoint for the barley-F. graminearum interaction. A late stage also identified between 96 and 144 hai, exhibiting development of hyphal mats, high DON accumulation, and a reduction in the number of transcripts observed. Our study provides a baseline and hypothesis-generating dataset in barley during F. graminearum infection and in other grasses during pathogen infection.

Anchoring and ordering NGS contig assemblies by population sequencing (POPSEQ).

Next-generation whole-genome shotgun assemblies of complex genomes are highly useful, but fail to link nearby sequence contigs with each other or provide a linear order of contigs along individual chromosomes. Here, we introduce a strategy based on sequencing progeny of a segregating population that allows de novo production of a genetically anchored linear assembly of the gene space of an organism. We demonstrate the power of the approach by reconstructing the chromosomal organization of the gene space of barley, a large, complex and highly repetitive 5.1 Gb genome. We evaluate the robustness of the new assembly by comparison to a recently released physical and genetic framework of the barley genome, and to various genetically ordered sequence-based genotypic datasets. The method is independent of the need for any prior sequence resources, and will enable rapid and cost-efficient establishment of powerful genomic information for many species.