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

The physical map of wheat chromosome 1BS provides insights into its gene space organization and evolution

Background

The wheat genome sequence is an essential tool for advanced genomic research and improvements. The generation of a high-quality wheat genome sequence is challenging due to its complex 17 Gb polyploid genome. To overcome these difficulties, sequencing through the construction of BAC-based physical maps of individual chromosomes is employed by the wheat genomics community. Here, we present the construction of the first comprehensive physical map of chromosome 1BS, and illustrate its unique gene space organization and evolution.

Results

Fingerprinted BAC clones were assembled into 57 long scaffolds, anchored and ordered with 2,438 markers, covering 83% of chromosome 1BS. The BAC-based chromosome 1BS physical map and gene order of the orthologous regions of model grass species were consistent, providing strong support for the reliability of the chromosome 1BS assembly. The gene space for chromosome 1BS spans the entire length of the chromosome arm, with 76% of the genes organized in small gene islands, accompanied by a two-fold increase in gene density from the centromere to the telomere.

Conclusions

This study provides new evidence on common and chromosome-specific features in the organization and evolution of the wheat genome, including a non-uniform distribution of gene density along the centromere-telomere axis, abundance of non-syntenic genes, the degree of colinearity with other grass genomes and a non-uniform size expansion along the centromere-telomere axis compared with other model cereal genomes. The high-quality physical map constructed in this study provides a solid basis for the assembly of a reference sequence of chromosome 1BS and for breeding applications.

In silico resurrection of the major vault protein suggests it is ancestral in modern eukaryotes

Vaults are very large oligomeric ribonucleoproteins conserved among a variety of species. The rat vault 3D structure shows an ovoid oligomeric particle, consisting of 78 major vault protein monomers, each of approximately 861 amino acids. Vaults are probably the largest ribonucleoprotein structures in eukaryote cells, being approximately 70 nm in length with a diameter of 40 nm—the size of three ribosomes and with a lumen capacity of 50 million ų. We use both protein sequences and inferred ancestral sequences for in silico virtual resurrection of tertiary and quaternary structures to search for vaults in a wide variety of eukaryotes. We find that the vault’s phylogenetic distribution is widespread in eukaryotes, but is apparently absent in some notable model organisms. Our conclusion from the distribution of vaults is that they were present in the last eukaryote common ancestor but they have apparently been lost from a number of groups including fungi, insects, and probably plants. Our approach of inferring ancestral 3D and quaternary structures is expected to be useful generally.

Characterization of phosphorus-regulated miR399 and miR827 and their isomirs in barley under phosphorus-sufficient and phosphorus-deficient conditions

BACKGROUND

miR399 and miR827 are both involved in conserved phosphorus (P) deficiency signalling pathways. miR399 targets the PHO2 gene encoding E2 enzyme that negatively regulates phosphate uptake and root-to-shoot allocation, while miR827 targets SPX-domain-containing genes that negatively regulate other P-responsive genes. However, the response of miR399 and miR827 to P conditions in barley has not been investigated.

RESULTS

In this study, we investigated the expression profiles of miR399 and miR827 in barley (Hordeum vulagre L.) under P-deficient and P-sufficient conditions. We identified 10 members of the miR399 family and one miR827 gene in barley, all of which were significantly up-regulated under deficient P. In addition, we found many isomirs of the miR399 family and miR827, most of which were also significantly up-regulated under deficient P. Several isomirs of miR399 members were found to be able to cleave their predicted targets in vivo. Surprisingly, a few small RNAs (sRNAs) derived from the single-stranded loops of the hairpin structures of MIR399b and MIR399e-1 were also found to be able to cleave their predicted targets in vivo. Many antisense sRNAs of miR399 and a few for miR827 were also detected, but they did not seem to be regulated by P. Intriguingly, the lowest expressed member, hvu-miR399k, had four-fold more antisense sRNAs than sense sRNAs, and furthermore under P sufficiency, the antisense sRNAs are more frequent than the sense sRNAs. We identified a potential regulatory network among miR399, its target HvPHO2 and target mimics HvIPS1 and HvIPS2 in barley under P-deficient and P-sufficient conditions.

CONCLUSIONS

Our data provide an important insight into the mechanistic regulation and function of miR399, miR827 and their isomirs in barley under different P conditions.

Resequencing data indicate a modest effect of domestication on diversity in barley: a cultigen with multiple origins

The levels of diversity and extent of linkage disequilibrium in cultivated species are largely determined by diversity in their wild progenitors. We report a comparison of nucleotide sequence diversity in wild and cultivated barley (Hordeum vulgare ssp. spontaneum and ssp. vulgare) at 7 nuclear loci totaling 9296bp, using sequence from Hordeum bulbosum to infer the ancestral state of mutations. The sample includes 36 accessions of cultivated barley, including 23 landraces (cultivated forms not subject to modern breeding) and 13 cultivated lines and genetic stocks compared to either 25 or 45 accessions of wild barley for the same loci. Estimates of nucleotide sequence diversity indicate that landraces retain >80% of the diversity in wild barley. The primary population structure in wild barley, which divides the species into eastern and western populations, is reflected in significant differentiation at all loci in wild accessions and at 3 of 7 loci in landraces. "Oriental" landraces have slightly higher diversity than "Occidental" landraces. Genetic assignment suggests more admixture from Occidental landraces into Oriental landraces than the converse, which may explain this difference. Based on θπ for silent sites, modern western cultivars have ~73% of the diversity found in landraces and ~71% of the diversity in wild barley.

A detailed gene expression study of the Miscanthus genus reveals changes in the transcriptome associated with the rejuvenation of spring rhizomes

BACKGROUND

The Miscanthus genus of perennial C4 grasses contains promising biofuel crops for temperate climates. However, few genomic resources exist for Miscanthus, which limits understanding of its interesting biology and future genetic improvement. A comprehensive catalog of expressed sequences were generated from a variety of Miscanthus species and tissue types, with an emphasis on characterizing gene expression changes in spring compared to fall rhizomes.

RESULTS

Illumina short read sequencing technology was used to produce transcriptome sequences from different tissues and organs during distinct developmental stages for multiple Miscanthus species, including Miscanthus sinensis, Miscanthus sacchariflorus, and their interspecific hybrid Miscanthus × giganteus. More than fifty billion base-pairs of Miscanthus transcript sequence were produced. Overall, 26,230 Sorghum gene models (i.e., ~ 96% of predicted Sorghum genes) had at least five Miscanthus reads mapped to them, suggesting that a large portion of the Miscanthus transcriptome is represented in this dataset. The Miscanthus × giganteus data was used to identify genes preferentially expressed in a single tissue, such as the spring rhizome, using Sorghum bicolor as a reference. Quantitative real-time PCR was used to verify examples of preferential expression predicted via RNA-Seq. Contiguous consensus transcript sequences were assembled for each species and annotated using InterProScan. Sequences from the assembled transcriptome were used to amplify genomic segments from a doubled haploid Miscanthus sinensis and from Miscanthus × giganteus to further disentangle the allelic and paralogous variations in genes.

CONCLUSIONS

This large expressed sequence tag collection creates a valuable resource for the study of Miscanthus biology by providing detailed gene sequence information and tissue preferred expression patterns. We have successfully generated a database of transcriptome assemblies and demonstrated its use in the study of genes of interest. Analysis of gene expression profiles revealed biological pathways that exhibit altered regulation in spring compared to fall rhizomes, which are consistent with their different physiological functions. The expression profiles of the subterranean rhizome provides a better understanding of the biological activities of the underground stem structures that are essentials for perenniality and the storage or remobilization of carbon and nutrient resources.

Emerging technologies advancing forage and turf grass genomics

Grassland is of major importance for agricultural production and provides valuable ecosystem services. Its impact is likely to rise in changing socio-economic and climatic environments. High yielding forage grass species are major components of sustainable grassland production. Understanding the genome structure and function of grassland species provides opportunities to accelerate crop improvement and thus to mitigate the future challenges of increased feed and food demand, scarcity of natural resources such as water and nutrients, and high product qualities. In this review, we will discuss a selection of technological developments that served as main drivers to generate new insights into the structure and function of nuclear genomes. Many of these technologies were originally developed in human or animal science and are now increasingly applied in plant genomics. Our main goal is to highlight the benefits of using these technologies for forage and turf grass genome research, to discuss their potentials and limitations as well as their relevance for future applications.

Kmasker--a tool for in silico prediction of single-copy FISH probes for the large-genome species Hordeum vulgare

Specific localization of large genomic fragments by fluorescence in situ hybridization (FISH) is challenging in large- genome plant species due to the high content of repetitive sequences. We report the automated work flow (Kmasker) for in silico extraction of unique genomic sequences of large genomic fragments suitable for FISH in barley. This method can be widely used for the integration of genetic and cytogenetic maps in plants and other species with large and complex genomes if the probe sequence (e.g. BACs, sequence contigs) and a low coverage (8-fold) of unassembled sequences of the species of interest are available. Kmasker has been made publicly available as a web tool at http://webblast.ipk-gatersleben.de/kmasker.

Unlocking Triticeae genomics to sustainably feed the future

The tribe Triticeae includes the major crops wheat and barley. Within the last few years, the whole genomes of four Triticeae species-barley, wheat, Tausch's goatgrass (Aegilops tauschii) and wild einkorn wheat (Triticum urartu)-have been sequenced. The availability of these genomic resources for Triticeae plants and innovative analytical applications using next-generation sequencing technologies are helping to revitalize our approaches in genetic work and to accelerate improvement of the Triticeae crops. Comparative genomics and integration of genomic resources from Triticeae plants and the model grass Brachypodium distachyon are aiding the discovery of new genes and functional analyses of genes in Triticeae crops. Innovative approaches and tools such as analysis of next-generation populations, evolutionary genomics and systems approaches with mathematical modeling are new strategies that will help us discover alleles for adaptive traits to future agronomic environments. In this review, we provide an update on genomic tools for use with Triticeae plants and Brachypodium and describe emerging approaches toward crop improvements in Triticeae.

Applications and challenges of next-generation sequencing in Brassica species

Next-generation sequencing (NGS) produces numerous (often millions) short DNA sequence reads, typically varying between 25 and 400 bp in length, at a relatively low cost and in a short time. This revolutionary technology is being increasingly applied in whole-genome, transcriptome, epigenome and small RNA sequencing, molecular marker and gene discovery, comparative and evolutionary genomics, and association studies. The Brassica genus comprises some of the most agro-economically important crops, providing abundant vegetables, condiments, fodder, oil and medicinal products. Many Brassica species have undergone the process of polyploidization, which makes their genomes exceptionally complex and can create difficulties in genomics research. NGS injects new vigor into Brassica research, yet also faces specific challenges in the analysis of complex crop genomes and traits. In this article, we review the advantages and limitations of different NGS technologies and their applications and challenges, using Brassica as an advanced model system for agronomically important, polyploid crops. Specifically, we focus on the use of NGS for genome resequencing, transcriptome sequencing, development of single-nucleotide polymorphism markers, and identification of novel microRNAs and their targets. We present trends and advances in NGS technology in relation to Brassica crop improvement, with wide application for sophisticated genomics research into agronomically important polyploid crops.

Fine mapping of the Rrs1 resistance locus against scald in two large populations derived from Spanish barley landraces

In two Spanish barley landraces with outstanding resistance to scald, the Rrs1 Rh4 locus was fine mapped including all known markers used in previous studies and closely linked markers were developed. Scald, caused by Rhynchosporium commune, is one of the most prevalent barley diseases worldwide. A search for new resistance sources revealed that Spanish landrace-derived lines SBCC145 and SBCC154 showed outstanding resistance to scald. They were crossed to susceptible cultivar Beatrix to create large DH-mapping populations of 522 and 416 DH lines that were scored for disease resistance in the greenhouse using two R. commune isolates. To ascertain the pattern of resistance, parents and reference barley lines with known scald resistance were phenotyped with a panel of differential R. commune isolates. Subpopulations were genotyped with the Illumina GoldenGate 1,536 SNP Assay and a large QTL in the centromeric region of chromosome 3H, known to harbour several scald resistance genes and/or alleles, was found in both populations. Five SNP markers closest to the QTL were converted into CAPS markers. These CAPS markers, together with informative SSR markers used in other scald studies, confirmed the presence of the Rrs1 locus. The panel of differential scald isolates indicated that the allele carried by both donors was Rrs1 Rh4 . The genetic distance between Rrs1 and its flanking markers was 1.2 cM (11_0010) proximally and 0.9 cM (11_0823) distally, which corresponds to a distance of just below 9 Mbp. The number and nature of scald resistance genes on chromosome 3H are discussed. The effective Rrs1 allele found and the closely linked markers developed are already useful tools for molecular breeding programs and provide a good step towards the identification of candidate genes.

Fine mapping and chromosome walking towards the Ror1 locus in barley (Hordeum vulgare L.)

The Ror1 gene was fine-mapped to the pericentric region of barley chromosome 1HL. Recessively inherited loss-of-function alleles of the barley (Hordeum vulgare) Mildew resistance locus o (Mlo) gene confer durable broad-spectrum disease resistance against the obligate biotrophic fungal powdery mildew pathogen Blumeria graminis f.sp. hordei. Previous genetic analyses revealed two barley genes, Ror1 and Ror2, that are Required for mlo-specified resistance and basal defence. While Ror2 was cloned and shown to encode a t-SNARE protein (syntaxin), the molecular nature or Ror1 remained elusive. Ror1 was previously mapped to the centromeric region of the long arm of barley chromosome 1H. Here, we narrowed the barley Ror1 interval to 0.18 cM and initiated a chromosome walk using barley yeast artificial chromosome (YAC) clones, next-generation DNA sequencing and fluorescence in situ hybridization. Two non-overlapping YAC contigs containing Ror1 flanking genes were identified. Despite a high degree of synteny observed between barley and the sequenced genomes of the grasses rice (Oryza sativa), Brachypodium distachyon and Sorghum bicolor across the wider chromosomal area, the genes in the YAC contigs showed extensive interspecific rearrangements in orientation and order. Consequently, the position of a Ror1 homolog in these species could not be precisely predicted, nor was a barley gene co-segregating with Ror1 identified. These factors have prevented the molecular identification of the Ror1 gene for the time being.

bex-db: Bioinformatics workbench for comprehensive analysis of barley-expressed genes

Barley (Hordeum vulgare) is one of the world's most important cereal crops. Although its large and complex genome has held back barley genomics for quite a while, the whole genome sequence was released in 2012 by the International Barley Genome Sequencing Consortium (IBSC). Moreover, more than 30,000 barley full-length cDNAs (FLcDNAs) are now available in the public domain. Here we present the Barley Gene Expression Database (bex-db: http://barleyflc.dna.affrc.go.jp/bexdb/index.html) as a repository of transcriptome data including the sequences and the expression profiles of barley genes resulting from microarray analysis. In addition to FLcDNA sequences, bex-db also contains partial sequences of more than 309,000 novel expressed sequence tags (ESTs). Users can browse the data via keyword, sequence homology and expression profile search options. A genome browser was also developed to display the chromosomal locations of barley FLcDNAs and wheat (Triticum aestivum) transcripts as well as Aegilops tauschii gene models on the IBSC genome sequence for future comparative analysis of orthologs among Triticeae species. The bex-db should provide a useful resource for further genomics studies and development of genome-based tools to enhance the progress of the genetic improvement of cereal crops.

Next generation characterisation of cereal genomes for marker discovery

Cereal crops form the bulk of the world’s food sources, and thus their importance cannot be understated. Crop breeding programs increasingly rely on high-resolution molecular genetic markers to accelerate the breeding process. The development of these markers is hampered by the complexity of some of the major cereal crop genomes, as well as the time and cost required. In this review, we address current and future methods available for the characterisation of cereal genomes, with an emphasis on faster and more cost effective approaches for genome sequencing and the development of markers for trait association and marker assisted selection (MAS) in crop breeding programs.

A physical map of the short arm of wheat chromosome 1A

Bread wheat (Triticum aestivum) has a large and highly repetitive genome which poses major technical challenges for its study. To aid map-based cloning and future genome sequencing projects, we constructed a BAC-based physical map of the short arm of wheat chromosome 1A (1AS). From the assembly of 25,918 high information content (HICF) fingerprints from a 1AS-specific BAC library, 715 physical contigs were produced that cover almost 99% of the estimated size of the chromosome arm. The 3,414 BAC clones constituting the minimum tiling path were end-sequenced. Using a gene microarray containing ∼40 K NCBI UniGene EST clusters, PCR marker screening and BAC end sequences, we arranged 160 physical contigs (97 Mb or 35.3% of the chromosome arm) in a virtual order based on synteny with Brachypodium, rice and sorghum. BAC end sequences and information from microarray hybridisation was used to anchor 3.8 Mbp of Illumina sequences from flow-sorted chromosome 1AS to BAC contigs. Comparison of genetic and synteny-based physical maps indicated that ∼50% of all genetic recombination is confined to 14% of the physical length of the chromosome arm in the distal region. The 1AS physical map provides a framework for future genetic mapping projects as well as the basis for complete sequencing of chromosome arm 1AS.

SNPMeta: SNP annotation and SNP metadata collection without a reference genome

The increase in availability of resequencing data is greatly accelerating SNP discovery and has facilitated the development of SNP genotyping assays. This, in turn, is increasing interest in annotation of individual SNPs. Currently, these data are only available through curation, or comparison to a reference genome. Many species lack a reference genome, but are still important genetic models or are significant species in agricultural production or natural ecosystems. For these species, it is possible to annotate SNPs through comparison with cDNA, or data from well-annotated genes in public repositories. We present SNPMeta, a tool which gathers information about SNPs by comparison with sequences present in GenBank databases. SNPMeta is able to annotate SNPs from contextual sequence in SNP assay designs, and SNPs discovered through genotyping by sequencing (GBS) approaches. However, SNPs discovered through GBS occur throughout the genome, rather than only in gene space, and therefore do not annotate at high rates. SNPMeta can therefore be used to annotate SNPs in nonmodel species or species that lack a reference genome. Annotations generated by SNPMeta are highly concordant with annotations that would be obtained from a reference genome.

Overexpression of cytokinin dehydrogenase genes in barley (Hordeum vulgare cv. Golden Promise) fundamentally affects morphology and fertility

Barley is one of the most important cereal crops grown worldwide. It has numerous applications, but its utility could potentially be extended by genetically manipulating its hormonal balances. To explore some of this potential we identified gene families of cytokinin dehydrogenases (CKX) and isopentenyl transferases, enzymes that respectively irreversibly degrade and synthesize cytokinin (CK) plant hormones, in the raw sequenced barley genome. We then examined their spatial and temporal expression patterns by immunostaining and qPCR. Two CKX-specific antibodies, anti-HvCKX1 and anti-HvCKX9, predominantly detect proteins in the aleurone layer of maturing grains and leaf vasculature, respectively. In addition, two selected CKX genes were used for stable, Agrobacterium tumefaciens-mediated transformation of the barley cultivar Golden Promise. The results show that constitutive overexpression of CKX causes morphological changes in barley plants and prevents their transition to flowering. In all independent transgenic lines roots proliferated more rapidly and root-to-shoot ratios were higher than in wild-type plants. Only one transgenic line, overexpressing CKX under the control of a promoter from a phosphate transporter gene, which is expressed more strongly in root tissue than in aerial parts, yielded progeny. Analysis of several T1-generation plants indicates that plants tend to compensate for effects of the transgene and restore CK homeostasis later during development. Depleted CK levels during early phases of development are restored by down-regulation of endogenous CKX genes and reinforced de novo biosynthesis of CKs.

Gramene 2013: comparative plant genomics resources

Gramene (http://www.gramene.org) is a curated online resource for comparative functional genomics in crops and model plant species, currently hosting 27 fully and 10 partially sequenced reference genomes in its build number 38. Its strength derives from the application of a phylogenetic framework for genome comparison and the use of ontologies to integrate structural and functional annotation data. Whole-genome alignments complemented by phylogenetic gene family trees help infer syntenic and orthologous relationships. Genetic variation data, sequences and genome mappings available for 10 species, including Arabidopsis, rice and maize, help infer putative variant effects on genes and transcripts. The pathways section also hosts 10 species-specific metabolic pathways databases developed in-house or by our collaborators using Pathway Tools software, which facilitates searches for pathway, reaction and metabolite annotations, and allows analyses of user-defined expression datasets. Recently, we released a Plant Reactome portal featuring 133 curated rice pathways. This portal will be expanded for Arabidopsis, maize and other plant species. We continue to provide genetic and QTL maps and marker datasets developed by crop researchers. The project provides a unique community platform to support scientific research in plant genomics including studies in evolution, genetics, plant breeding, molecular biology, biochemistry and systems biology.

Arabidopsis meiotic crossover hot spots overlap with H2A.Z nucleosomes at gene promoters

PRDM9 directs human meiotic crossover hot spots to intergenic sequence motifs, whereas budding yeast hot spots overlap regions of low nucleosome density (LND) in gene promoters. To investigate hot spots in plants, which lack PRDM9, we used coalescent analysis of genetic variation in Arabidopsis thaliana. Crossovers increased toward gene promoters and terminators, and hot spots were associated with active chromatin modifications, including H2A.Z, histone H3 Lys4 trimethylation (H3K4me3), LND and low DNA methylation. Hot spot-enriched A-rich and CTT-repeat DNA motifs occurred upstream and downstream, respectively, of transcriptional start sites. Crossovers were asymmetric around promoters and were most frequent over CTT-repeat motifs and H2A.Z nucleosomes. Pollen typing, segregation and cytogenetic analysis showed decreased numbers of crossovers in the arp6 H2A.Z deposition mutant at multiple scales. During meiosis, H2A.Z forms overlapping chromosomal foci with the DMC1 and RAD51 recombinases. As arp6 reduced the number of DMC1 or RAD51 foci, H2A.Z may promote the formation or processing of meiotic DNA double-strand breaks. We propose that gene chromatin ancestrally designates hot spots within eukaryotes and PRDM9 is a derived state within vertebrates.

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.

Making open data work for plant scientists

Despite the clear demand for open data sharing, its implementation within plant science is still limited. This is, at least in part, because open data-sharing raises several unanswered questions and challenges to current research practices. In this commentary, some of the challenges encountered by plant researchers at the bench when generating, interpreting, and attempting to disseminate their data have been highlighted. The difficulties involved in sharing sequencing, transcriptomics, proteomics, and metabolomics data are reviewed. The benefits and drawbacks of three data-sharing venues currently available to plant scientists are identified and assessed: (i) journal publication; (ii) university repositories; and (iii) community and project-specific databases. It is concluded that community and project-specific databases are the most useful to researchers interested in effective data sharing, since these databases are explicitly created to meet the researchers' needs, support extensive curation, and embody a heightened awareness of what it takes to make data reuseable by others. Such bottom-up and community-driven approaches need to be valued by the research community, supported by publishers, and provided with long-term sustainable support by funding bodies and government. At the same time, these databases need to be linked to generic databases where possible, in order to be discoverable to the majority of researchers and thus promote effective and efficient data sharing. As we look forward to a future that embraces open access to data and publications, it is essential that data policies, data curation, data integration, data infrastructure, and data funding are linked together so as to foster data access and research productivity.

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.

Quantitative trait loci-dependent analysis of a gene co-expression network associated with Fusarium head blight resistance in bread wheat (Triticum aestivum L.)

Background

Fusarium head blight (FHB) caused by Fusarium graminearum Schwabe is one of the most prevalent diseases of wheat (Triticum aestivum L.) and other small grain cereals. Resistance against the fungus is quantitative and more than 100 quantitative trait loci (QTL) have been described. Two well-validated and highly reproducible QTL, Fhb1 and Qfhs.ifa-5A have been widely investigated, but to date the underlying genes have not been identified.

Results

We have investigated a gene co-expression network activated in response to F. graminearum using RNA-seq data from near-isogenic lines, harboring either the resistant or the susceptible allele for Fhb1 and Qfhs.ifa-5A. The network identified pathogen-responsive modules, which were enriched for differentially expressed genes between genotypes or different time points after inoculation with the pathogen. Central gene analysis identified transcripts associated with either QTL within the network. Moreover, we present a detailed gene expression analysis of four gene families (glucanases, NBS-LRR, WRKY transcription factors and UDP-glycosyltransferases), which take prominent roles in the pathogen response.

Conclusions

A combination of a network-driven approach and differential gene expression analysis identified genes and pathways associated with Fhb1 and Qfhs.ifa-5A. We find G-protein coupled receptor kinases and biosynthesis genes for jasmonate and ethylene earlier induced for Fhb1. Similarly, we find genes involved in the biosynthesis and metabolism of riboflavin more abundant for Qfhs.ifa-5A.

Large-scale collection and analysis of full-length cDNAs from Brachypodium distachyon and integration with Pooideae sequence resources

A comprehensive collection of full-length cDNAs is essential for correct structural gene annotation and functional analyses of genes. We constructed a mixed full-length cDNA library from 21 different tissues of Brachypodium distachyon Bd21, and obtained 78,163 high quality expressed sequence tags (ESTs) from both ends of ca. 40,000 clones (including 16,079 contigs). We updated gene structure annotations of Brachypodium genes based on full-length cDNA sequences in comparison with the latest publicly available annotations. About 10,000 non-redundant gene models were supported by full-length cDNAs; ca. 6,000 showed some transcription unit modifications. We also found ca. 580 novel gene models, including 362 newly identified in Bd21. Using the updated transcription start sites, we searched a total of 580 plant cis-motifs in the −3 kb promoter regions and determined a genome-wide Brachypodium promoter architecture. Furthermore, we integrated the Brachypodium full-length cDNAs and updated gene structures with available sequence resources in wheat and barley in a web-accessible database, the RIKEN Brachypodium FL cDNA database. The database represents a “one-stop” information resource for all genomic information in the Pooideae, facilitating functional analysis of genes in this model grass plant and seamless knowledge transfer to the Triticeae crops.

Application of genotyping-by-sequencing on semiconductor sequencing platforms: a comparison of genetic and reference-based marker ordering in barley

The rapid development of next-generation sequencing platforms has enabled the use of sequencing for routine genotyping across a range of genetics studies and breeding applications. Genotyping-by-sequencing (GBS), a low-cost, reduced representation sequencing method, is becoming a common approach for whole-genome marker profiling in many species. With quickly developing sequencing technologies, adapting current GBS methodologies to new platforms will leverage these advancements for future studies. To test new semiconductor sequencing platforms for GBS, we genotyped a barley recombinant inbred line (RIL) population. Based on a previous GBS approach, we designed bar code and adapter sets for the Ion Torrent platforms. Four sets of 24-plex libraries were constructed consisting of 94 RILs and the two parents and sequenced on two Ion platforms. In parallel, a 96-plex library of the same RILs was sequenced on the Illumina HiSeq 2000. We applied two different computational pipelines to analyze sequencing data; the reference-independent TASSEL pipeline and a reference-based pipeline using SAMtools. Sequence contigs positioned on the integrated physical and genetic map were used for read mapping and variant calling. We found high agreement in genotype calls between the different platforms and high concordance between genetic and reference-based marker order. There was, however, paucity in the number of SNP that were jointly discovered by the different pipelines indicating a strong effect of alignment and filtering parameters on SNP discovery. We show the utility of the current barley genome assembly as a framework for developing very low-cost genetic maps, facilitating high resolution genetic mapping and negating the need for developing de novo genetic maps for future studies in barley. Through demonstration of GBS on semiconductor sequencing platforms, we conclude that the GBS approach is amenable to a range of platforms and can easily be modified as new sequencing technologies, analysis tools and genomic resources develop.

Development of SNP markers for genes of the phenylpropanoid pathway and their association to kernel and malting traits in barley

BACKGROUND

Flavonoids are an important class of secondary compounds in angiosperms. Next to certain biological functions in plants, they play a role in the brewing process and have an effect on taste, color and aroma of beer. The aim of this study was to reveal the haplotype diversity of candidate genes involved in the phenylpropanoid biosynthesis pathway in cultivated barley varieties (Hordeum vulgare L.) and to determine associations to kernel and malting quality parameters.

RESULTS

Five genes encoding phenylalanine ammonia-lyase (PAL), cinnamate 4-hydroxylase (C4H), chalcone synthase (CHS), flavanone 3-hydroxylase (F3H) and dihydroflavonol reductase (DFR) of the phenylpropanoid biosynthesis pathway were partially resequenced in 16 diverse barley reference genotypes. Their localization in the barley genome, their genetic structure, and their genetic variation e.g. single nucleotide polymorphism (SNP) and Insertion/Deletion (InDel) patterns were revealed. In total, 130 SNPs and seven InDels were detected. Of these, 21 polymorphisms were converted into high-throughput pyrosequencing markers. The resulting SNP and haplotype patterns were used to calculate associations with kernel and malting quality parameters.

CONCLUSIONS

SNP patterns were found to be highly variable for the investigated genes. The developed high-throughput markers are applicable for assessing the genetic variability and for the determination of haplotype patterns in a set of barley accessions. The candidate genes PAL, C4H and F3H were shown to be associated to several malting properties like glassiness (PAL), viscosity (C4H) or to final attenuation (F3H).

Next-generation survey sequencing and the molecular organization of wheat chromosome 6B

Common wheat (Triticum aestivum L.) is one of the most important cereals in the world. To improve wheat quality and productivity, the genomic sequence of wheat must be determined. The large genome size (∼17 Gb/1 C) and the hexaploid status of wheat have hampered the genome sequencing of wheat. However, flow sorting of individual chromosomes has allowed us to purify and separately shotgun-sequence a pair of telocentric chromosomes. Here, we describe a result from the survey sequencing of wheat chromosome 6B (914 Mb/1 C) using massively parallel 454 pyrosequencing. From the 4.94 and 5.51 Gb shotgun sequence data from the two chromosome arms of 6BS and 6BL, 235 and 273 Mb sequences were assembled to cover ∼55.6 and 54.9% of the total genomic regions, respectively. Repetitive sequences composed 77 and 86% of the assembled sequences on 6BS and 6BL, respectively. Within the assembled sequences, we predicted a total of 4798 non-repetitive gene loci with the evidence of expression from the wheat transcriptome data. The numbers and chromosomal distribution patterns of the genes for tRNAs and microRNAs in wheat 6B were investigated, and the results suggested a significant involvement of DNA transposon diffusion in the evolution of these non-protein-coding RNA genes. A comparative analysis of the genomic sequences of wheat 6B and monocot plants clearly indicated the evolutionary conservation of gene contents.

TAL effector nucleases induce mutations at a pre-selected location in the genome of primary barley transformants

Transcription activator-like effector nucleases (TALENs) enable targeted mutagenesis in a variety of organisms. The primary advantage of TALENs over other sequence-specific nucleases, namely zinc finger nucleases and meganucleases, lies in their ease of assembly, reliability of function, and their broad targeting range. Here we report the assembly of several TALENs for a specific genomic locus in barley. The cleavage activity of individual TALENs was first tested in vivo using a yeast-based, single-strand annealing assay. The most efficient TALEN was then selected for barley transformation. Analysis of the resulting transformants showed that TALEN-induced double strand breaks led to the introduction of short deletions at the target site. Additional analysis revealed that each barley transformant contained a range of different mutations, indicating that mutations occurred independently in different cells.

Major QTL for Fusarium crown rot resistance in a barley landrace

Fusarium crown rot (FCR) is a serious cereal disease in semi-arid regions worldwide. In assisting the effort of breeding cultivars with enhanced resistance, we identified several barley genotypes with high levels of FCR resistance. One of these genotypes, AWCS079 which is a barley landrace originating from Japan, was investigated by developing and assessing three populations of recombinant inbred lines. Two QTL, one located on the long arm of chromosome 1H (designated as Qcrs.cpi-1H) and the other on 3HL (designated as Qcrs.cpi-3H), were found to be responsible for the FCR resistance of this genotype. Qcrs.cpi-1H is novel as no other FCR loci have been reported on this chromosome arm. Qcrs.cpi-3H co-located with a reduced height (Rht) locus and the effectiveness of the former was significantly affected by the latter. The total phenotypic variance explained by these two QTL was over 60 %. Significant effects were detected for each of the QTL in each of the three populations assessed. The existence of these loci with major effects should not only facilitate breeding and exploitation of FCR-resistant barley cultivars but also their further characterization based on fine mapping and map-based gene cloning.

Modelling shifts in agroclimate and crop cultivar response under climate change

This paper aims: (i) to identify at national scale areas where crop yield formation is currently most prone to climate-induced stresses, (ii) to evaluate how the severity of these stresses is likely to develop in time and space, and (iii) to appraise and quantify the performance of two strategies for adapting crop cultivation to a wide range of (uncertain) climate change projections. To this end we made use of extensive climate, crop, and soil data, and of two modelling tools: N-AgriCLIM and the WOFOST crop simulation model. N-AgriCLIM was developed for the automatic generation of indicators describing basic agroclimatic conditions and was applied over the whole of Finland. WOFOST was used to simulate detailed crop responses at four representative locations. N-AgriCLIM calculations have been performed nationally for 3829 grid boxes at a 10 × 10 km resolution and for 32 climate scenarios. Ranges of projected shifts in indicator values for heat, drought and other crop-relevant stresses across the scenarios vary widely – so do the spatial patterns of change. Overall, under reference climate the most risk-prone areas for spring cereals are found in south-west Finland, shifting to south-east Finland towards the end of this century. Conditions for grass are likely to improve. WOFOST simulation results suggest that CO₂ fertilization and adjusted sowing combined can lead to small yield increases of current barley cultivars under most climate scenarios on favourable soils, but not under extreme climate scenarios and poor soils. This information can be valuable for appraising alternative adaptation strategies. It facilitates the identification of regions in which climatic changes might be rapid or otherwise notable for crop production, requiring a more detailed evaluation of adaptation measures. The results also suggest that utilizing the diversity of cultivar responses seems beneficial given the high uncertainty in climate change projections.

Transgenic barley: a prospective tool for biotechnology and agriculture

Barley (Hordeum vulgare L.) is one of the founder crops of agriculture, and today it is the fourth most important cereal grain worldwide. Barley is used as malt in brewing and distilling industry, as an additive for animal feed, and as a component of various food and bread for human consumption. Progress in stable genetic transformation of barley ensures a potential for improvement of its agronomic performance or use of barley in various biotechnological and industrial applications. Recently, barley grain has been successfully used in molecular farming as a promising bioreactor adapted for production of human therapeutic proteins or animal vaccines. In addition to development of reliable transformation technologies, an extensive amount of various barley genetic resources and tools such as sequence data, microarrays, genetic maps, and databases has been generated. Current status on barley transformation technologies including gene transfer techniques, targets, and progeny stabilization, recent trials for improvement of agricultural traits and performance of barley, especially in relation to increased biotic and abiotic stress tolerance, and potential use of barley grain as a protein production platform have been reviewed in this study. Overall, barley represents a promising tool for both agricultural and biotechnological transgenic approaches, and is considered an ancient but rediscovered crop as a model industrial platform for molecular farming.

Variation in the interaction between alleles of HvAPETALA2 and microRNA172 determines the density of grains on the barley inflorescence

Within the cereal grasses, variation in inflorescence architecture results in a conspicuous morphological diversity that in crop species influences the yield of cereal grains. Although significant progress has been made in identifying some of the genes underlying this variation in maize and rice, in the temperate cereals, a group that includes wheat, barley, and rye, only the dosage-dependent and highly pleiotropic Q locus in hexaploid wheat has been molecularly characterized. Here we show that the characteristic variation in the density of grains along the inflorescence, or spike, of modern cultivated barley (Hordeum vulgare) is largely the consequence of a perturbed interaction between microRNA172 and its corresponding binding site in the mRNA of an APELATA2 (AP2)-like transcription factor, HvAP2. We used genome-wide association and biparental mapping to identify HvAP2. By comparing inflorescence development and HvAP2 transcript abundance in an extreme dense-spike mutant and its nearly isogenic WT line, we show that HvAP2 turnover driven by microRNA 172 regulates the length of a critical developmental window that is required for elongation of the inflorescence internodes. Our data indicate that this heterochronic change, an altered timing of developmental events caused by specific temporal variation in the efficiency of HvAP2 turnover, leads to the striking differences in the size and shape of the barley spike.

Plant genome sequencing - applications for crop improvement

It is over 10 years since the genome sequence of the first crop was published. Since then, the number of crop genomes sequenced each year has increased steadily. The amazing pace at which genome sequences are becoming available is largely due to the improvement in sequencing technologies both in terms of cost and speed. Modern sequencing technologies allow the sequencing of multiple cultivars of smaller crop genomes at a reasonable cost. Though many of the published genomes are considered incomplete, they nevertheless have proved a valuable tool to understand important crop traits such as fruit ripening, grain traits and flowering time adaptation.

Identification and distribution of Puroindoline b-2 variant gene homologs in Hordeum

The barley hordoindoline genes (Hina and Hinb) are homologous to the wheat puroindoline genes (Pina and Pinb). These genes are involved in grain hardness, which is an important quality for barley processing. We identified novel variants of Hina and Hinb in 10 wild Hordeum species (H. bogdanii, H. brachyantherum, H. bulbosum, H. chilense, H. comosum, H. marinum, H. murinum, H. patagonicum, H. pusillum, and H. roshevitzii) covering all Hordeum genomes and preliminarily named them Hinc. These nucleotide sequences were highly similar to those of Puroindoline b-2 variant genes (Pinb-2v) and were located on chromosome 7I in H. chilense. The Hinc genes in H. bogdanii, H. bulbosum, H. patagonicum, and H. roshevitzii were pseudogenes possessing in-frame stop codons. We also found a partial Hinc sequence in H. murinum. This gene was not found in cultivated barley and H. vulgare subsp. spontaneum. The phylogenetic tree of Gsp-1, Hin, and Pin genes demonstrates that Hinc and Pinb-2v genes formed one cluster. Therefore, we considered that Hinc and Pinb-2v genes shared a common ancestral gene and were homologous to each other. We also studied the evolutional process of Gsp-1, Hin, and Pin genes. Our results suggested that Gsp-1 might be the most closely related to a putative ancestral gene on Ha locus.

Glutathione S-transferase (GST) family in barley: identification of members, enzyme activity, and gene expression pattern

Barley (Hordeum vulgare) is one of the most important cereals in many developing countries where drought stress considerably diminishes agricultural production. Glutathione S-transferases (GSTs EC 2.5.1.18) are multifunctional enzymes which play a crucial role in cellular detoxification and oxidative stress tolerance. In this study, 84 GST genes were identified in barley by a comprehensive in silico approach. Sequence alignment and phylogenetic analysis grouped these HvGST proteins in eight classes. The largest numbers of the HvGST genes (50) were included in the Tau class followed by 21 genes in Phi, five in Zeta, two in DHAR, two in EF1G, two in Lambda, and one each in TCHQD and Theta classes. Phylogenetic analysis of the putative GSTs from Arabidopsis, rice, and barley indicated that major functional diversification within the GST family predated the monocot/dicot divergence. However, intra-specious duplication seems to be common. Expression patterns of five GST genes from Phi and Tau classes were investigated in three barley genotypes (Yusof [drought-tolerant], Moroc9-75 [drought-sensitive], and HS1 [wild ecotype]) under control and drought-stressed conditions, during the vegetative stage. All investigated genes were up-regulated significantly under drought stress and/or showed a higher level of transcripts in the tolerant cultivar. Additionally, GST enzyme activity was superior in Yusof and induced in the extreme-drought-treated leaves, while it was not changed in Moroc9-75 under drought conditions. Moreover, the lowest and highest levels of lipid peroxidation were observed in the Yusof and Moroc9-75 cultivars, respectively. Based on the achieved results, detoxification and antioxidant activity of GSTs might be considered an important factor in the drought tolerance of barley genotypes for further investigations.

Retrotransposon replication in plants

Retrotransposons comprise the bulk of large plant genomes, replicating via an RNA intermediate whereby the original, integrated element remains in place. Of the two main orders, the LTR retrotransposons considerably outnumber the LINEs. LINEs integrate into target sites simultaneously with the RNA transcript being copied into cDNA by target-primed reverse transcription. LTR retrotransposon replication is basically equivalent to the intracellular phase of retroviral life cycles. The envelope gene giving extracellular mobility to retroviruses is in fact widespread in plants and their retrotransposons. Evolutionary analyses of the retrotransposons and retroviruses suggest that both form an ancient monophyletic group. The particular adaptations of LTR retrotransposons to plant life cycles enabling their success remain to be clarified.

A wheat homologue of PHYTOCLOCK 1 is a candidate gene conferring the early heading phenotype to einkorn wheat

An X-ray mutant showing an early flowering phenotype has been identified in einkorn wheat (Triticum monococcum L.), for which a major QTL for heading time was previously mapped in the telomeric region on the long arm of chromosome 3A. Recent advances in Triticeae genomics revealed that the gene order in this region is highly conserved between wheat and barley. Thus, we adopted a hypothetical gene order in barley, the so-called GenomeZipper, to develop DNA markers for fine mapping the target gene in wheat. We identified three genes tightly linked to the early heading phenotype. PCR analysis revealed that early-flowering is associated with the deletion of two genes in the mutant. Of the two deleted genes, one is an ortholog of the LUX ARRHYTHMO (LUX)/PHYTOCLOCK 1 (PCL1) gene found in Arabidopsis, which regulates the circadian clock and flowering time. We found distorted expression patterns of two clock genes (TOC1 and LHY) in the einkorn pcl1 deletion mutant as was reported for the Arabidopsis lux mutant. Transcript accumulation levels of photoperiod-response related genes, a photoperiod sensitivity gene (Ppd-1) and two wheat CONSTANS-like genes (WCO1 and TaHd1), were significantly higher in the einkorn wheat mutant. In addition, transcripts of the wheat florigen gene (WFT) accumulated temporally under short-day conditions in the einkorn wheat mutant. These results suggest that deletion of WPCL1 leads to abnormally higher expression of Ppd-1, resulting in the accumulation of WFT transcripts that triggers flowering even under short-day conditions. Our observations from gene mapping, gene deletions, and expression levels of flowering related genes strongly suggest that WPCL1 is the most likely candidate gene for controlling the early flowering phenotype in the einkorn wheat mutant.

Nitrogen-metabolism related genes in barley - haplotype diversity, linkage mapping and associations with malting and kernel quality parameters

Background

Several studies report about intra-specific trait variation of nitrogen-metabolism related traits, such as N(itrogen)-use efficiency, protein content, N-storage and remobilization in barley and related grass species. The goal of this study was to assess the intra-specific genetic diversity present in primary N-metabolism genes of barley and to investigate the associations of the detected haplotype diversity with malting and kernel quality related traits.

Results

Partial sequences of five genes related to N-metabolism in barley (Hordeum vulgare L.) were obtained, i.e. nitrate reductase 1, glutamine synthetase 2, ferredoxin-dependent glutamate synthase, aspartate aminotransferase and asparaginase. Two to five haplotypes in each gene were discovered in a set of 190 various varieties. The development of 33 SNP markers allowed the genotyping of all these barley varieties consisting of spring and winter types. Furthermore, these markers could be mapped in several doubled haploid populations. Cluster analysis based on haplotypes revealed a more uniform pattern of the spring barleys as compared to the winter barleys. Based on linear model approaches associations to several malting and kernel quality traits including soluble N and protein were identified.

Conclusions

A study was conducted to investigate the presence of sequence variation of several genes related to the primary N-metabolism in barley. The detected diversity could be related to particular phenotypic traits. Specific differences between spring and winter barleys most likely reflect different breeding aims. The developed markers can be used as tool for further genetic studies and marker-assisted selection during breeding of barley.

Insights into the loblolly pine genome: characterization of BAC and fosmid sequences

Despite their prevalence and importance, the genome sequences of loblolly pine, Norway spruce, and white spruce, three ecologically and economically important conifer species, are just becoming available to the research community. Following the completion of these large assemblies, annotation efforts will be undertaken to characterize the reference sequences. Accurate annotation of these ancient genomes would be aided by a comprehensive repeat library; however, few studies have generated enough sequence to fully evaluate and catalog their non-genic content. In this paper, two sets of loblolly pine genomic sequence, 103 previously assembled BACs and 90,954 newly sequenced and assembled fosmid scaffolds, were analyzed. Together, this sequence represents 280 Mbp (roughly 1% of the loblolly pine genome) and one of the most comprehensive studies of repetitive elements and genes in a gymnosperm species. A combination of homology and de novo methodologies were applied to identify both conserved and novel repeats. Similarity analysis estimated a repetitive content of 27% that included both full and partial elements. When combined with the de novo investigation, the estimate increased to almost 86%. Over 60% of the repetitive sequence consists of full or partial LTR (long terminal repeat) retrotransposons. Through de novo approaches, 6,270 novel, full-length transposable element families and 9,415 sub-families were identified. Among those 6,270 families, 82% were annotated as single-copy. Several of the novel, high-copy families are described here, with the largest, PtPiedmont, comprising 133 full-length copies. In addition to repeats, analysis of the coding region reported 23 full-length eukaryotic orthologous proteins (KOGS) and another 29 novel or orthologous genes. These discoveries, along with other genomic resources, will be used to annotate conifer genomes and address long-standing questions about gymnosperm evolution.

Effects of barley chromosome addition to wheat on behavior and development of Locusta migratoria nymphs

Locusta migratoria feeds on various Poaceae plants but barley. Barley genes related to feeding deterrence may be useful for developing novel resistant crops. We investigated the effects of barley cultivar Betzes, wheat cultivar Chinese Spring (CS), and six barley chromosome disomic addition lines of wheat (2H-7H) on locomotor activity, feeding behavior, survival and development of L. migratoria nymphs. Locomotor activity was similar in nymphs kept with wheat and 2H-7H in an actograph, whereas it was generally high in those kept with barely. No-choice and choice feeding tests suggested that barley genes related to inhibition of feeding by L. migratoria are located on barley chromosomes 5H and 6H and those related to the palatability of plants on chromosomes 2H, 5H and 6H. Rearing experiments suggested the presence of barley genes negatively affecting the survival and growth of locust nymphs on chromosomes 5H and 2H, respectively, and the effects are phase-dependent.

HvLUX1 is a candidate gene underlying the early maturity 10 locus in barley: phylogeny, diversity, and interactions with the circadian clock and photoperiodic pathways

Photoperiodic flowering is a major factor determining crop performance and is controlled by interactions between environmental signals and the circadian clock. We proposed Hvlux1, an ortholog of the Arabidopsis circadian gene LUX ARRHYTHMO, as a candidate underlying the early maturity 10 (eam10) locus in barley (Hordeum vulgare L.). The link between eam10 and Hvlux1 was discovered using high-throughput sequencing of enriched libraries and segregation analysis. We conducted functional, phylogenetic, and diversity studies of eam10 and HvLUX1 to understand the genetic control of photoperiod response in barley and to characterize the evolution of LUX-like genes within barley and across monocots and eudicots. We demonstrate that eam10 causes circadian defects and interacts with the photoperiod response gene Ppd-H1 to accelerate flowering under long and short days. The results of phylogenetic and diversity analyses indicate that HvLUX1 was under purifying selection, duplicated at the base of the grass clade, and diverged independently of LUX-like genes in other plant lineages. Taken together, these findings contribute to improved understanding of the barley circadian clock, its interaction with the photoperiod pathway, and evolution of circadian systems in barley and across monocots and eudicots.

Analysis of chromosomal polymorphism in barley (Hordeum vulgare L. ssp. vulgare) and between H. vulgare and H. chilense using three-color fluorescence in situ hybridization (FISH)

The aim of the present work was to study chromosomal polymorphism within cultivated barley (Hordeum vulgare ssp. vulgare) using three-color fluorescence in situ hybridization (FISH). The physical distribution of the most frequently used, highly repetitive DNA sequences (GAA)7 specific for pericentromeric heterochromatic regions, the ribosomal DNA clone pTa71, specific for the 45S rDNA, and the barley-specific telomere-associated sequence HvT01, was investigated to reveal genetic diversity in metaphase spreads of ten barley genotypes with diverse geographical origin, growth habit and row number. A wild relative of barley, Hordeum chilense was also studied in order to compare the polymorphism between and within Hordeum species. Significant differences in the hybridization patterns of all three DNA probes could be detected between the two related species, but only probes pTa71 and HvT01 showed variation in the intensity and/or position of hybridization sites among genotypes of H. vulgare ssp. vulgare. The extent of polymorphism was less than that earlier reported for molecular markers and was restricted to the long chromosome arms, with differences between the chromosomes. 1H and 3H proved to be the most variable chromosomes and 4H and 6H the most conserved.

Characterisation and analysis of the Aegilops sharonensis transcriptome, a wild relative of wheat in the Sitopsis section

Aegilops sharonensis Eig (Sharon goatgrass) is a wild diploid relative of wheat within the Sitopsis section of Aegilops. This species represents an untapped reservoir of genetic diversity for traits of agronomic importance, especially as a source of novel disease resistance. To gain a foothold in this genetic resource, we sequenced the cDNA from leaf tissue of two geographically distinct Ae. sharonensis accessions (1644 and 2232) using the 454 Life Sciences platform. We compared the results of two different assembly programs using different parameter sets to generate 13 distinct assemblies in an attempt to maximize representation of the gene space in de novo transcriptome assembly. The most sensitive assembly (71,029 contigs; N50 674 nts) retrieved 18,684 unique best reciprocal BLAST hits (BRBH) against six previously characterised grass proteomes while the most specific assembly (30,609 contigs; N50 815 nts) retrieved 15,687 BRBH. We combined these two assemblies into a set of 62,243 non-redundant sequences and identified 139 belonging to plant disease resistance genes of the nucleotide binding leucine-rich repeat class. Based on the non-redundant sequences, we predicted 37,743 single nucleotide polymorphisms (SNP), equivalent to one per 1,142 bp. We estimated the level of heterozygosity as 1.6% in accession 1644 and 30.1% in 2232. The Ae. sharonensis leaf transcriptome provides a rich source of sequence and SNPs for this wild wheat relative. These sequences can be used with existing monocot genome sequences and EST sequence collections (e.g. barley, Brachypodium, wheat, rice, maize and Sorghum) to assist with genetic and physical mapping and candidate gene identification in Ae. sharonensis. These resources provide an initial framework to further build on and characterise the genetic and genomic structure of Ae. sharonensis.

Leaf proteome alterations in the context of physiological and morphological responses to drought and heat stress in barley (Hordeum vulgare L.)

The objective of this study was to identify barley leaf proteins differentially regulated in response to drought and heat and the combined stresses in context of the morphological and physiological changes that also occur. The Syrian landrace Arta and the Australian cultivar Keel were subjected to drought, high temperature, or a combination of both treatments starting at heading. Changes in the leaf proteome were identified using differential gel electrophoresis and mass spectrometry. The drought treatment caused strong reductions of biomass and yield, while photosynthetic performance and the proteome were not significantly changed. In contrast, the heat treatment and the combination of heat and drought reduced photosynthetic performance and caused changes of the leaf proteome. The proteomic analysis identified 99 protein spots differentially regulated in response to heat treatment, 14 of which were regulated in a genotype-specific manner. Differentially regulated proteins predominantly had functions in photosynthesis, but also in detoxification, energy metabolism, and protein biosynthesis. The analysis indicated that de novo protein biosynthesis, protein quality control mediated by chaperones and proteases, and the use of alternative energy resources, i.e. glycolysis, play important roles in adaptation to heat stress. In addition, genetic variation identified in the proteome, in plant growth and photosynthetic performance in response to drought and heat represent stress adaption mechanisms to be exploited in future crop breeding efforts.

Homoeologous relationship of rye chromosome arms as detected with wheat PLUG markers

Based on the similarity in gene structure between rice and wheat, the polymerase chain reaction (PCR)-based landmark unique gene (PLUG) system enabled us to design primer sets that amplify wheat genic sequences including introns. From the previously reported wheat PLUG markers, we chose 144 markers that are distributed on different chromosomes and in known chromosomal regions (bins) to obtain rye-specific PCR-based markers. We conducted PCR with the 144 primer sets and the template of the Imperial rye genomic DNA and found that 131 (91.0%) primer sets successfully amplified PCR products. Of the 131 PLUG markers, 110 (76.4%) markers showed rye-specific PCR amplification with or without restriction enzyme digestion. We assigned 79 of the 110 markers to seven rye chromosomes (1R to 7R) using seven wheat-rye (cv. Imperial) chromosome addition and substitution lines: 12 to 1R, 8 to 2R, 11 to 3R, 8 to 4R, 16 to 5R, 12 to 6R, and 12 to 7R. Furthermore, we located their positions on the short or long (L) chromosome arm, using 13 Imperial rye telosomic lines of common wheat (except for 3RL). Referring to the chromosome bin locations of the 79 PLUG markers in wheat, we deduced the syntenic relationships between rye and wheat chromosomes. We also discussed chromosomal rearrangements in the rye genome with reference to the cytologically visible chromosomal gaps.

Conserved synteny-based anchoring of the barley genome physical map

Gene order is largely collinear in the small-grained cereals, a feature which has proved helpful in both marker development and positional cloning. The accuracy of a virtual gene order map ("genome zipper") for barley (Hordeum vulgare), developed by combining a genetic map of this species with a large number of gene locations obtained from the maps constructed in other grass species, was evaluated here both at the genome-wide level and at the fine scale in a representative segment of the genome. Comparing the whole genome "genome zipper" maps with a genetic map developed by using transcript-derived markers, yielded an accuracy of >94 %. The fine-scale comparison involved a 14 cM segment of chromosome arm 2HL. One hundred twenty-eight genes of the "genome zipper" interval were analysed. Over 95 % (45/47) of the polymorphic markers were genetically mapped and allocated to the expected region of 2HL, following the predicted order. A further 80 of the 128 genes were assigned to the correct chromosome arm 2HL by analysis of wheat-barley addition lines. All 128 gene-based markers developed were used to probe a barley bacterial artificial chromosome (BAC) library, delivering 26 BAC contigs from which all except two were anchored to the targeted zipper interval. The results demonstrate that the gene order predicted by the "genome zipper" is remarkably accurate and that the "genome zipper" represents a highly efficient informational resource for the systematic identification of gene-based markers and subsequent physical map anchoring of the barley genome.

Conserved synteny-based anchoring of the barley genome physical map

Gene order is largely collinear in the small-grained cereals, a feature which has proved helpful in both marker development and positional cloning. The accuracy of a virtual gene order map ("genome zipper") for barley (Hordeum vulgare), developed by combining a genetic map of this species with a large number of gene locations obtained from the maps constructed in other grass species, was evaluated here both at the genome-wide level and at the fine scale in a representative segment of the genome. Comparing the whole genome "genome zipper" maps with a genetic map developed by using transcript-derived markers, yielded an accuracy of >94 %. The fine-scale comparison involved a 14 cM segment of chromosome arm 2HL. One hundred twenty-eight genes of the "genome zipper" interval were analysed. Over 95 % (45/47) of the polymorphic markers were genetically mapped and allocated to the expected region of 2HL, following the predicted order. A further 80 of the 128 genes were assigned to the correct chromosome arm 2HL by analysis of wheat-barley addition lines. All 128 gene-based markers developed were used to probe a barley bacterial artificial chromosome (BAC) library, delivering 26 BAC contigs from which all except two were anchored to the targeted zipper interval. The results demonstrate that the gene order predicted by the "genome zipper" is remarkably accurate and that the "genome zipper" represents a highly efficient informational resource for the systematic identification of gene-based markers and subsequent physical map anchoring of the barley genome.

Arguments for standardizing transposable element annotation in plant genomes

Whole genome sequence assemblies have been generated for many plants. Annotation of transposable elements (TEs), which constitute the major proportion of genomes and play a significant role in epigenome alterations under stress, has not been given equal importance to that of genes. In this opinion article, we argue that the lack of focus dedicated to the fine-scale characterization of repeat fractions and the absence of consistent methods for their annotation impede our ability to critically understand the influence of TEs on the epigenome with implications in gene expression and non-Mendelian inheritance. Major structural changes occur over an evolutionary time scale. However, epigenetic regulation mediated by TEs can happen in a single generation, thus emphasizing the need for their standardized annotation.