A whole-genome snapshot of 454 sequences exposes the composition of the barley genome and provides evidence for parallel evolution of genome size in wheat and barley

The genomes of barley and wheat, two of the world's most important crops, are very large and complex due to their high content of repetitive DNA. In order to obtain a whole-genome sequence sample, we performed two runs of 454 (GS20) sequencing on genomic DNA of barley cv. Morex, which yielded approximately 1% of a haploid genome equivalent. Almost 60% of the sequences comprised known transposable element (TE) families, and another 9% represented novel repetitive sequences. We also discovered high amounts of low-complexity DNA and non-genic low-copy DNA. We identified almost 2300 protein coding gene sequences and more than 660 putative conserved non-coding sequences. Comparison of the 454 reads with previously published genomic sequences suggested that TE families are distributed unequally along chromosomes. This was confirmed by in situ hybridizations of selected TEs. A comparison of these data for the barley genome with a large sample of publicly available wheat sequences showed that several TE families that are highly abundant in wheat are absent from the barley genome. This finding implies that the TE composition of their genomes differs dramatically, despite their very similar genome size and their close phylogenetic relationship.     

A structural and evolutionary analysis of a dispersed repetitive sequence

A family of dispersed repetitive sequences (Hch1) which is present in the genome of the wild barley Hordeum chilense was studied in detail. Hch1 sequences are found both as part of short tandem arrays and dispersed throughout the H. chilense chromosomes. Subcloning of sections of the sequence reveals that it is composed of unrelated classes of sequences which can also be found separately in other genomic locations. Analysis of these sequences in the genomes of wheat and two other wild barley species strongly suggests that specific amplifications and arrangements of the repeated sequences have taken place during speciation. Nucleotide sequence analysis fails to detect, in their entirity, the features shown by plant transposons.     

Molecular evolution of receptor-like kinase genes in hexaploid wheat. Independent evolution of orthologs after polyploidization and mechanisms of local rearrangements at paralogous loci

Hexaploid wheat is a young polyploid species and represents a good model to study mechanisms of gene evolution after polyploidization. Recent studies at the scale of the whole genome have suggested rapid genomic changes after polyploidization but so far the rearrangements that have occurred in terms of gene content and organization have not been analyzed at the microlevel in wheat. Here, we have isolated members of a receptor kinase (Lrk) gene family in hexaploid and diploid wheat, Aegilops tauschii, and barley (Hordeum vulgare). Phylogenetic analysis has allowed us to establish evolutionary relationships (orthology versus paralogy) between the different members of this gene family in wheat as well as with Lrk genes from barley. It also demonstrated that the sequences of the homoeologous Lrk genes evolved independently after polyploidization. In addition, we found evidence for gene loss during the evolution of wheat and barley. Analysis of large genomic fragments isolated from nonorthologous Lrk loci showed a high conservation of the gene content and gene organization at these loci on the homoeologous group 1 chromosomes of wheat and barley. Finally, sequence comparison of two paralogous fragments of chromosome 1B showed a large number of local events (sequence duplications, deletions, and insertions), which reveal rearrangements and mechanisms for genome enlargement at the microlevel.     

Updating of transposable element annotations from large wheat genomic sequences reveals diverse activities and gene associations

Triticeae species (including wheat, barley and rye) have huge and complex genomes due to polyploidization and a high content of transposable elements (TEs). TEs are known to play a major role in the structure and evolutionary dynamics of Triticeae genomes. During the last 5 years, substantial stretches of contiguous genomic sequence from various species of Triticeae have been generated, making it necessary to update and standardize TE annotations and nomenclature. In this study we propose standard procedures for these tasks, based on structure, nucleic acid and protein sequence homologies. We report statistical analyses of TE composition and distribution in large blocks of genomic sequences from wheat and barley. Altogether, 3.8 Mb of wheat sequence available in the databases was analyzed or re-analyzed, and compared with 1.3 Mb of re-annotated genomic sequences from barley. The wheat sequences were relatively gene-rich (one gene per 23.9 kb), although wheat gene-derived sequences represented only 7.8% (159 elements) of the total, while the remainder mainly comprised coding sequences found in TEs (54.7%, 751 elements). Class I elements [mainly long terminal repeat (LTR) retrotransposons] accounted for the major proportion of TEs, in terms of sequence length as well as element number (83.6% and 498, respectively). In addition, we show that the gene-rich sequences of wheat genome A seem to have a higher TE content than those of genomes B and D, or of barley gene-rich sequences. Moreover, among the various TE groups, MITEs were most often associated with genes: 43.1% of MITEs fell into this category. Finally, the TRIM and copia elements were shown to be the most active TEs in the wheat genome. The implications of these results for the evolution of diploid and polyploid wheat species are discussed. Electronic Supplementary Material Supplementary material is available for this article at     

Effective Isolation of Retrotransposons and Repetitive DNA Families from the Wheat Genome

New classes of repetitive DNA elements were effectively identified by isolating small fragments of the elements from the wheat genome. A wheat A genome library was constructed from Triticum monococcum by degenerate cleavage with EcoO109I, the recognition sites of which consisted of 5'-PuGGNCCPy-3'multi-sequences. Three novel repetitive sequences pTm6, pTm69 and pTm58 derived from the A genome were screened and tested for high copy number using a blotting approach. pTm6 showed identity with integrase domains of the barley Ty1-Copia-retrotransposon BARE-1 and pTm58 showed similarity to the barley Ty3-gypsy-like retrotransposon Romani. pTm69, however, constituted a tandem array with useful genomic specificities, but did not share any identity with known repetitive elements. This study also sought to isolate wheat D-genome-specific repetitive elements regardless of the level of methylation, by genomic subtraction. Total genomic DNA of Aegilops tauschii was cleaved into short fragments with a methylation-insensitive 4 bp cutter, Mbol, and then common DNA sequences between Ae. tauschii and Triticum turgidum were subtracted by annealing with excess T. turgidum genomic DNA. The D genome repetitive sequence pAt1 was isolated and used to identify an additional novel repetitive sequence family from wheat bacterial artificial chromosomes with a size range of 1 395-1 850 bp. The methods successfully led pathfinding of two unique repetitive families.     

Resistance gene analogs in barley and their relationship to rust resistance genes.

Regions of amino acid conservation in the NBS domain of NBS-LRR resistance proteins facilitated the PCR isolation of eight resistance gene analog (RGA) sequences from genomic DNA of rice, barley, and Aegilops tauschii. These clones and other RGAs previously isolated from maize, rice, and wheat were assigned to 13 classes by DNA-sequence comparison and by their patterns of hybridisation to restricted barley DNA. Using a doubled-haploid mapping population, probes from 12 RGA classes were used to map 17 loci in the barley genome. Many of these probes have been used for mapping in wheat, and the collective data indicate that the positions of orthologous RGAs are conserved between barley and wheat. RGA loci were identified in the vicinity of barley leaf rust resistance loci Rph4, Rph7, and Rph10. Recombinants were identified between RGA loci and Rph7 and Rph10, while a cluster of RGA sequences detected by probe 5.2 cosegregated with Rph4 in 55 F2 lines.     

Analysis of Intraspecies Diversity in Wheat and Barley Genomes Identifies Breakpoints of Ancient Haplotypes and Provides Insight into the Structure of Diploid and Hexaploid Triticeae Gene Pools

A large number of wheat (Triticum aestivum) and barley (Hordeum vulgare) varieties have evolved in agricultural ecosystems since domestication. Because of the large, repetitive genomes of these Triticeae crops, sequence information is limited and molecular differences between modern varieties are poorly understood. To study intraspecies genomic diversity, we compared large genomic sequences at the Lr34 locus of the wheat varieties Chinese Spring, Renan, and Glenlea, and diploid wheat Aegilops tauschii. Additionally, we compared the barley loci Vrs1 and Rym4 of the varieties Morex, Cebada Capa, and Haruna Nijo. Molecular dating showed that the wheat D genome haplotypes diverged only a few thousand years ago, while some barley and Ae. tauschii haplotypes diverged more than 500,000 years ago. This suggests gene flow from wild barley relatives after domestication, whereas this was rare or absent in the D genome of hexaploid wheat. In some segments, the compared haplotypes were very similar to each other, but for two varieties each at the Rym4 and Lr34 loci, sequence conservation showed a breakpoint that separates a highly conserved from a less conserved segment. We interpret this as recombination breakpoints of two ancient haplotypes, indicating that the Triticeae genomes are a heterogeneous and variable mosaic of haplotype fragments. Analysis of insertions and deletions showed that large events caused by transposable element insertions, illegitimate recombination, or unequal crossing over were relatively rare. Most insertions and deletions were small and caused by template slippage in short homopolymers of only a few base pairs in size. Such frequent polymorphisms could be exploited for future molecular marker development.     

Organization of the histone H3 genes in soybean, barley and wheat

Several variants of the replacement histone H3 genes from soybean, barley and wheat have been cloned and sequenced. Analysis of segregating populations in barley and soybean, as well as analysis of clones isolated from a soybean genomic library, suggested that these genes are dispersed throughout the genome. Several genes contain introns located in similar positions, but of different lengths and sequence. Comparison of mRNA levels in different tissues revealed that the intron-containing and intronless genes have different expression patterns. The distribution of the introns in the histone H3 genes across several plant species suggests that some of the introns might have been lost during the evolution of the gene family. Sequence divergence among introns and gene-flanking sequences in cloned gene variants allowed us to use them as specific probes for localizing individual gene copies and analyzing the genomic distribution of these variants across a range of genotypes.Journal paper No. J-16127 of the Iowa Agriculture and Home Economics Experiment Station, Ames, IowaMention of a trademark or proprietary product does not constitute a guarantee or warranty of the product by the United States Department of Agriculture and does not imply its approval to the exclusion of other products that may be suitable     

Analysis of a contiguous 211 kb sequence in diploid wheat (Triticum monococcum L.) reveals multiple mechanisms of genome evolution

In plant species with large genomes such as wheat or barley, genome organization at the level of DNA sequence is largely unknown. The largest sequences that are publicly accessible so far from Triticeae genomes are two 60 kb and 66 kb intervals from barley. Here, we report on the analysis of a 211 kb contiguous DNA sequence from diploid wheat (Triticum monococcum L.). Five putative genes were identified, two of which show similarity to disease resistance genes. Three of the five genes are clustered in a 31 kb gene-enriched island while the two others are separated from the cluster and from each other by large stretches of repetitive DNA. About 70% of the contig is comprised of several classes of transposable elements. Ten different types of retrotransposons were identified, most of them forming a pattern of nested insertions similar to those found in maize and barley. Evidence was found for major deletion, insertion and duplication events within the analysed region, suggesting multiple mechanisms of genome evolution in addition to retrotransposon amplification. Seven types of foldback transposons, an element class previously not described for wheat genomes, were characterized. One such element was found to be closely associated with genes in several Triticeae species and may therefore be of use for the identification of gene-rich regions in these species.     

Integrating cereal genomics to support innovation in the Triticeae

The genomic resources of small grain cereals that include some of the most important crop species such as wheat, barley, and rye are attaining a level of completion that now is contributing to new structural and functional studies as well as refining molecular marker development and mapping strategies for increasing the efficiency of breeding processes. The integration of new efforts to obtain reference sequences in bread wheat and barley, in particular, is accelerating the acquisition and interpretation of genome-level analyses in both of these major crops.     

Rapid Elimination of Low-Copy DNA Sequences in Polyploid Wheat: A Possible Mechanism for Differentiation of Homoeologous Chromosomes

To study genome evolution in allopolyploid plants, we analyzed polyploid wheats and their diploid progenitors for the occurrence of 16 low-copy chromosome- or genome-specific sequences isolated from hexaploid wheat. Based on their occurrence in the diploid species, we classified the sequences into two groups: group I, found in only one of the three diploid progenitors of hexaploid wheat, and group II, found in all three diploid progenitors. The absence of group II sequences from one genome of tetraploid wheat and from two genomes of hexaploid wheat indicates their specific elimination from these genomes at the polyploid level. Analysis of a newly synthesized amphiploid, having a genomic constitution analogous to that of hexaploid wheat, revealed a pattern of sequence elimination similar to the one found in hexaploid wheat. Apparently, speciation through allopolyploidy is accompanied by a rapid, nonrandom elimination of specific, low-copy, probably noncoding DNA sequences at the early stages of allopolyploidization, resulting in further divergence of homoeologous chromosomes (partially homologous chromosomes of different genomes carrying the same order of gene loci). We suggest that such genomic changes may provide the physical basis for the diploid-like meiotic behavior of polyploid wheat.     

DNA condensation with polyamines. II. Electron microscopic studies

Approximately 75% of the wheat and rye genomes consist of repeated sequence DNA. Three-quarters of the non-repeated or few copy sequences in wheat are less than 1000 base-pairs long, whilst in rye approximately half of the non-repeated or few copy sequences are in this size class. Most of the remaining non-repeated or few copy sequences appear to be a few thousand base-pairs long.In this paper a somewhat novel approach has been used to quantitatively analyse the linear organisation of the large proportion of repeated sequence DNA as well as the non-repeated DNA in the wheat and rye genomes. Repeated sequences in the genomes of oats, barley, wheat and rye have been used as probes to distinguish and isolate four different groups of repeated sequences and their neighbouring sequences from the wheat and rye genomes. Radioactively labelled wheat or rye DNA fragments ranging from 200 to over 9000 nucleotides long were incubated separately with large excesses of denatured unlabelled oats, barley, wheat and rye DNAs to C_ot values which enable all the repeated sequences of the unlabelled DNA to renature. The following parameters were then determined from the proportions of total labelled DNA in fragments which had at least partially renatured. (1) The proportions of the repeated sequences in the labelled DNAs that were able to hybridise to each unlabelled DNA; (2) the mean distance apart of the hybridising sequences on the longer labelled fragments; and (3) the proportion of the genome in which the hybridising sequences were concentrated. Analysis of these results, together with those of separate experiments designed to quantitatively estimate the nature of sequences unable to reanneal with the repeated sequences of each of the probe DNAs, have enabled schematic maps to be drawn which show how the repeated and non-repeated sequences are arranged in the wheat and rye genomes.Both genomes are constructed from millions of relatively short sequences, most of them considerably shorter than 3000 base-pairs. This structure was recognised because adjacent sequences can be distinguished by their frequency of repetition (i.e. repeated or non-repeated) or by their evolutionary origin. Approximately 40 to 45% of the wheat genome and 30 to 35% of the rye genome consists of short non-repeated sequences interspersed between short repeated sequences. Approximately 50% of the wheat genome and 60% of the rye genome consists of tandemly arranged repeated sequences of different evolutionary origins. It is postulated that much of this complex repeated sequence DNA could have arisen from amplification of compound sequences, each containing repeated and non-repeated sequence DNA.Short repeated sequences with a number average length of around 200 base-pairs and which occupy about 20% of the wheat and rye genomes are related to repeated sequences also found in oats and barley. They are concentrated in 60 to 70% of the wheat and rye genomes, being interspersed with different short repeated sequences and a significant proportion of the short non-repeated sequences.Rye chromosomes contain more DNA than wheat chromosomes. This is principally, but not entirely, due to additional repeated sequence DNA. Many quantitative changes appear to have occurred in both genomes, possibly affecting most families of repeated sequences, since wheat and rye diverged from a common ancestor. Both species contain species-specific repeated sequences (24% of rye genome; 16% of wheat genome) but a large proportion of these are closely interspersed with repeated sequences found in both genomes.     

A New Class of Wheat Gliadin Genes and Proteins

The utility of mining DNA sequence data to understand the structure and expression of cereal prolamin genes is demonstrated by the identification of a new class of wheat prolamins. This previously unrecognized wheat prolamin class, given the name δ-gliadins, is the most direct ortholog of barley γ3-hordeins. Phylogenetic analysis shows that the orthologous δ-gliadins and γ3-hordeins form a distinct prolamin branch that existed separate from the γ-gliadins and γ-hordeins in an ancestral Triticeae prior to the branching of wheat and barley. The expressed δ-gliadins are encoded by a single gene in each of the hexaploid wheat genomes. This single δ-gliadin/γ3-hordein ortholog may be a general feature of the Triticeae tribe since examination of ESTs from three barley cultivars also confirms a single γ3-hordein gene. Analysis of ESTs and cDNAs shows that the genes are expressed in at least five hexaploid wheat cultivars in addition to diploids Triticum monococcum and Aegilops tauschii. The latter two sequences also allow assignment of the δ-gliadin genes to the A and D genomes, respectively, with the third sequence type assumed to be from the B genome. Two wheat cultivars for which there are sufficient ESTs show different patterns of expression, i.e., with cv Chinese Spring expressing the genes from the A and B genomes, while cv Recital has ESTs from the A and D genomes. Genomic sequences of Chinese Spring show that the D genome gene is inactivated by tandem premature stop codons. A fourth δ-gliadin sequence occurs in the D genome of both Chinese Spring and Ae. tauschii, but no ESTs match this sequence and limited genomic sequences indicates a pseudogene containing frame shifts and premature stop codons. Sequencing of BACs covering a 3 Mb region from Ae. tauschii locates the δ-gliadin gene to the complex Gli-1 plus Glu-3 region on chromosome 1.     

Development of 5006 Full-Length CDNAs in Barley: A Tool for Accessing Cereal Genomics Resources

A collection of 5006 full-length (FL) cDNA sequences was developed in barley. Fifteen mRNA samples from various organs and treatments were pooled to develop a cDNA library using the CAP trapper method. More than 60% of the clones were confirmed to have complete coding sequences, based on comparison with rice amino acid and UniProt sequences. Blastn homologies (E<1E-5) to rice genes and Arabidopsis genes were 89 and 47%, respectively. Of the 5028 possible amino acid sequences derived from the 5006 FLcDNAs, 4032 (80.2%) were classified into 1678 GreenPhyl multigenic families. There were 555 cDNAs showing low homology to both rice and Arabidopsis. Gene ontology annotation by InterProScan indicated that many of these cDNAs (71%) have no known molecular functions and may be unique to barley. The cDNAs showed high homology to Barley 1 GeneChip oligo probes (81%) and the wheat gene index (84%). The high homology between FLcDNAs (27%) and mapped barley expressed sequence tag enabled assigning linkage map positions to 151–233 FLcDNAs on each of the seven barley chromosomes. These comprehensive barley FLcDNAs provide strong platform to connect pre-existing genomic and genetic resources and accelerate gene identification and genome analysis in barley and related species.Key words: full-length cDNA, Hordeum vulgare, mRNA, gene ontology     

Sequencing of 15 622 gene‐bearing BACs clarifies the gene‐dense regions of the barley genome

Barley (Hordeum vulgare L.) possesses a large and highly repetitive genome of 5.1 Gb that has hindered the development of a complete sequence. In 2012, the International Barley Sequencing Consortium released a resource integrating whole‐genome shotgun sequences with a physical and genetic framework. However, because only 6278 bacterial artificial chromosome (BACs) in the physical map were sequenced, fine structure was limited. To gain access to the gene‐containing portion of the barley genome at high resolution, we identified and sequenced 15 622 BACs representing the minimal tiling path of 72 052 physical‐mapped gene‐bearing BACs. This generated ~1.7 Gb of genomic sequence containing an estimated 2/3 of all Morex barley genes. Exploration of these sequenced BACs revealed that although distal ends of chromosomes contain most of the gene‐enriched BACs and are characterized by high recombination rates, there are also gene‐dense regions with suppressed recombination. We made use of published map‐anchored sequence data from Aegilops tauschii to develop a synteny viewer between barley and the ancestor of the wheat D‐genome. Except for some notable inversions, there is a high level of collinearity between the two species. The software HarvEST:Barley provides facile access to BAC sequences and their annotations, along with the barley–Ae. tauschii synteny viewer. These BAC sequences constitute a resource to improve the efficiency of marker development, map‐based cloning, and comparative genomics in barley and related crops. Additional knowledge about regions of the barley genome that are gene‐dense but low recombination is particularly relevant.     

Organization and primary nucleotide sequence of 5S rRNA genes in barley Hordeum vulgare

A BamHI DNA fragment of 301 bp corresponding to the main repeating unit of 5S rRNA was isolated from barley genomic DNA. The primary nucleotide sequence of this fragment was determined and a high level of homology was found between coding sequences of 5S rRNA genes of barley, wheat and rye. At the same time, spacer's nucleotide sequences of different species of cereals were changed dramatically. At least two types of 5S rRNA tandem repeats of 301 and 450 bp were found in barley genome. Polymorphism for restriction fragment length in 5S rRNA repeats allowed to discriminate between all barley varieties used in this work.     

Sorting the wheat from the chaff: identifying miRNAs in genomic survey sequences of Triticum aestivum chromosome 1AL

Individual chromosome-based studies of bread wheat are beginning to provide valuable structural and functional information about one of the world's most important crops. As new genome sequences become available, identifying miRNA coding sequences is arguably as important a task as annotating protein coding sequences, but one that is not as well developed. We compared conservation-based identification of conserved miRNAs in 1.5× coverage survey sequences of wheat chromosome 1AL with a predictive method based on pre-miRNA hairpin structure alone. In total, 42 sequences expected to encode conserved miRNAs were identified on chromosome 1AL, including members of several miRNA families that have not previously been reported to be expressed in T. aestivum. In addition, we demonstrate that a number of sequences previously annotated as novel wheat miRNAs are closely related to transposable elements, particularly Miniature Inverted Terminal repeat Elements (MITEs). Some of these TE-miRNAs may well have a functional role, but separating true miRNA coding sequences from TEs in genomic sequences is far from straightforward. We propose a strategy for annotation to minimize the risk of mis-identifying TE sequences as miRNAs.     

Assembly and analysis of a qingke reference genome demonstrate its close genetic relation to modern cultivated barley

Qingke, the local name of hulless barley in the Tibetan Plateau, is a staple food for Tibetans. The availability of its reference genome sequences could be useful for studies on breeding and molecular evolution. Taking advantage of the third‐generation sequencer (PacBio), we de novo assembled a 4.84‐Gb genome sequence of qingke, cv. Zangqing320 and anchored a 4.59‐Gb sequence to seven chromosomes. Of the 46,787 annotated ‘high‐confidence’ genes, 31 564 were validated by RNA‐sequencing data of 39 wild and cultivated barley genotypes with wide genetic diversity, and the results were also confirmed by nonredundant protein database from NCBI. As some gaps in the reference genome of Morex were covered in the reference genome of Zangqing320 by PacBio reads, we believe that the Zangqing320 genome provides the useful supplements for the Morex genome. Using the qingke genome as a reference, we conducted a genome comparison, revealing a close genetic relationship between a hulled barley (cv. Morex) and a hulless barley (cv. Zangqing320), which is strongly supported by the low‐diversity regions in the two genomes. Considering the origin of Morex from its breeding pedigree, we then demonstrated a close genomic relationship between modern cultivated barley and qingke. Given this genomic relationship and the large genetic diversity between qingke and modern cultivated barley, we propose that qingke could provide elite genes for barley improvement.