Structural organization of the barley D-hordein locus in comparison with its orthologous regions of wheat genomes

Beneficial glycaemic effects of high-amylose barley bread compared to wheat bread in type 2 diabetes

BACKGROUND

Cereals foods with a high content of dietary fibres or amylose have potential to lower postprandial glucose levels. Optimisation of cereal foods may improve management of type 2 diabetes (T2D).

METHODS

We investigated the impact on 4 h postprandial glucose responses given as incremental area under curve (iAUC) of bread made of either 50% RNAi-based (genetically modified) amylose-only barley flour (AmOn) (and 50% wheat flour), 50% hulless barley flour (and 50% wheat flour) or 75% hulless barley (and 25% wheat flour) in subjects with T2D compared with 100% wheat flour bread.

DESIGN

Twenty adults with T2D were randomly allocated to one of four breads at four separate visits. We measured fasting and 4 h postprandial responses of glucose, insulin, glucagon, triacylglycerol (TG), free fatty acids (FFA), glucagon-like peptide-1 (GLP-1) and gastric inhibitory polypeptide (GIP). Mixed model ANOVA was used to examine the differences.

RESULTS

Bread made from 50% AmOn lowered the 4 h postprandial glucose by 34%, 27%, 23% (P < 0.05) compared with 100% wheat, 50% or 75% hulless barley, respectively. Bread made from 75% hulless barley reduced the postprandial glucose response (iAUC) by 11% (P < 0.05) compared to 100% wheat bread. Postprandial insulin responses (iAUC) were reduced for 50% AmOn compared with 100% wheat and 50% hulless barley and for 75% hulless compared to 50% hulless barley bread (P < 0.05). 4 h postprandial glucagon (tAUC) did not differ between the four bread types (P > 0.05). Lower postprandial GIP (iAUC) was observed after all barley breads compared to 100% wheat (P < 0.05), whereas no difference was seen in postprandial GLP-1. Postprandial TG and FFA (tAUC) were difficult to judge due to differences in fasting values.

CONCLUSIONS

Bread made by replacing wheat flour with either 50% high-amylose or 75% hulless barley flour lowered postprandial glucose responses compared to 100% wheat bread indicating a beneficial impact on glucose regulation in T2D subjects. This trial was registered at clinicaltrials.gov as NCT04646746.

Low Gluten Beers Contain Variable Gluten and Immunogenic Epitope Content

Gluten content labels inform food choice and people practicing a gluten-free diet rely upon them to avoid illness. The regulations differ between jurisdictions, especially concerning fermented foodstuffs such as beer. Gluten abundance is typically measured using ELISAs, which have come into question when testing fermented or hydrolysed foodstuffs such as beer. Mass spectrometry can be used to directly identify gluten peptides and reveal false negatives recorded by ELISA. In this survey of gluten in control and gluten-free beers, gluten protein fragments that contain known immunogenic epitopes were detected using liquid chromatography-mass spectrometry in multiple beers that claim to be gluten-free and have sufficiently low gluten content, as measured by ELISA, to qualify as being gluten-free in some jurisdictions. In fact, several purportedly gluten-free beers showed equivalent or higher hordein content than some of the untreated, control beers. The shortcomings of ELISAs for beer gluten testing are summarised, the mismatch between ELISA and mass spectrometry results are explored, and the suitability of existing regulations as they pertain to the gluten content in fermented foods in different jurisdictions are discussed.

Tracing 100 million years of grass genome evolutionary plasticity

Grasses derive from a family of monocotyledonous plants that includes crops of major economic importance such as wheat, rice, sorghum and barley, sharing a common ancestor some 100 million years ago. The genomic attributes of plant adaptation remain obscure and the consequences of recurrent whole genome duplications (WGD) or polyploidization events, a major force in plant evolution, remain largely speculative. We conducted a comparative analysis of omics data from ten grass species to unveil structural (inversions, fusions, fissions, duplications, substitutions) and regulatory (expression and methylation) basis of genome plasticity, as possible attributes of plant long lasting evolution and adaptation. The present study demonstrates that diverged polyploid lineages sharing a common WGD event often present the same patterns of structural changes and evolutionary dynamics, but these patterns are difficult to generalize across independent WGD events as a result of non-WGD factors such as selection and domestication of crops. Polyploidy is unequivocally linked to the evolutionary success of grasses during the past 100 million years, although it remains difficult to attribute this success to particular genomic consequences of polyploidization, suggesting that polyploids harness the potential of genome duplication, at least partially, in lineage-specific ways. Overall, the present study clearly demonstrates that post-polyploidization reprogramming is more complex than traditionally reported in investigating single species and calls for a critical and comprehensive comparison across independently polyploidized lineages.

The CC-NB-LRR protein BSR1 from Brachypodium confers resistance to Barley stripe mosaic virus in gramineous plants by recognising TGB1 movement protein

Although some nucleotide binding, leucine-rich repeat immune receptor (NLR) proteins conferring resistance to specific viruses have been identified in dicot plants, NLR proteins involved in viral resistance have not been described in monocots. We have used map-based cloning to isolate the CC-NB-LRR (CNL) Barley stripe mosaic virus (BSMV) resistance gene barley stripe resistance 1 (BSR1) from Brachypodium distachyon Bd3-1 inbred line. Stable BSR1 transgenic Brachypodium line Bd21-3, barley (Golden Promise) and wheat (Kenong 199) plants developed resistance against BSMV ND18 strain. Allelic variation analyses indicated that BSR1 is present in several Brachypodium accessions collected from countries in the Middle East. Protein domain swaps revealed that the intact LRR domain and the C-terminus of BSR1 are required for resistance. BSR1 interacts with the BSMV ND18 TGB1 protein in planta and shows temperature-sensitive antiviral resistance. The R³⁹⁰ and T³⁹² residues of TGB1_ND (ND18 strain) and the G¹⁹⁶ and K¹⁹⁷ residues within the BSR1 P-loop motif are key amino acids required for immune activation. BSR1 is the first cloned virus resistance gene encoding a typical CNL protein in monocots, highlighting the utility of the Brachypodium model for isolation and analysis of agronomically important genes for crop improvement.

Proteolytic activities and profiles as useful traits to select barley cultivars for beer production

Barley malting depends on hydrolytic enzymes that degrade storage macromolecules. Identifying barley cultivars with proteolytic activity that guarantees appropriate foaming, flavor, and aroma in the beer is of great importance. In this work, the proteolytic activity and profiles of brewing malt from Mexican barley cultivars were analyzed. Data showed that Cys- (at 50°C) and Ser-proteases (at 70°C) are the major contributors to proteolytic activity during mashing. Essential amino acids, necessary for fermentation and production of good flavor and aroma in beer, were detected at the end of mashing. According to our results, Mexican cultivar HV2005-19 exhibits similar proteolytic activities as those from cultivar Metcalfe, which is one of the most utilized for the brewing industry. Moreover, we propose Cys- and Ser-proteases as biochemical markers during mashing at 50 and 70°C, respectively, to select barley cultivars for beer production. PRACTICAL APPLICATIONS: Proteolytic activity, which depends on activation and de novo synthesis of proteases in the aleurone layer of barley seeds, is crucial in beer production. Identifying new barley varieties that have optimal proteolytic activities is of great interest for genetic improvement programs. In this study, we propose the variety HV2005-19 as a genotype with Cys- and Ser-proteases activity similar to that from Metcalfe, which is a top variety in the brewing industry.

Hordein Accumulation in Developing Barley Grains

The temporal pattern of accumulation of hordein storage proteins in developing barley grains was studied by enzyme-linked immunosorbent assay (ELISA), western blot and liquid chromatography tandem mass spectrometry (LC-MS/MS). Hordein accumulation was compared to the pattern seen for two abundant control proteins, serpin Z4 (an early accumulator) and lipid transferase protein (LTP1, a late accumulator). Hordeins were detected from 6 days post-anthesis (DPA) and peaked at 30 DPA. Changes in fresh weight indicate that desiccation begins at 20 DPA and by 37 DPA fresh weight had decreased by 35%. ELISA analysis of hordein content, expressed on a protein basis, increased to a maximum at 30 DPA followed by a 17% decrease by 37 DPA. The accumulation of 39 tryptic and 29 chymotryptic hordein peptides representing all classes of hordein was studied by LC-MS/MS. Most peptides increased to a maximum at 30 DPA, and either remained at the maximum or did not decrease significantly. Only five tryptic peptides, members of the related B1- and γ1-hordeins decreased significantly by 21-51% at 37 DPA. Thus, the concentration of some specific peptides was reduced while remaining members of the same family were not affected. The N-terminal signal region was removed by proteolysis during co-translation. In addition to a suite of previously characterized hordeins, two novel barley B-hordein isoforms mapping to wheat low molecular weight glutenins (LMW-GS-like B-hordeins), and two avenin-like proteins (ALPs) sharing homology with wheat ALPs, were identified. These identified isoforms have not previously been mapped in the barley genome. Cereal storage proteins provide significant nutritional content for human consumption and seed germination. In barley, the bulk of the storage proteins comprise the hordein family and the final hordein concentration affects the quality of baked and brewed products. It is therefore important to study the accumulation of hordeins as this knowledge may assist plant breeding for improved health outcomes (by minimizing triggering of detrimental immune responses), nutrition and food processing properties.

Analysis of novel high-molecular-weight prolamins from Leymus multicaulis (Kar. et Kir.) Tzvelev and L. chinensis (Trin. ex Bunge) Tzvelev

Nine novel high-molecular-weight prolamins (HMW-prolamins) were isolated from Leymus multicaulis and L. chinensis. Based on the structure of the repetitive domains, all nine genes were classified as D-hordeins but not high-molecular-weight glutenin subunits (HMW-GSs) that have been previously isolated in Leymus spp. Four genes, Lmul 1.2, 2.4, 2.7, and Lchi 2.5 were verified by bacterial expression, whereas the other five sequences (1.3 types) were classified as pseudogenes. The four Leymus D-hordein proteins had longer N-termini than those of Hordeum spp. [116/118 vs. 110 amino acid (AA) residues], whereas three (Lmul 1.2, 2.4, and 2.7) contained shorter N-termini than those of the Ps. juncea (116 vs. 118 AA residues). Furthermore, Lmul 1.2 was identified as the smallest D-hordein, and Lmul 1.2 and 2.7 had an additional cysteines. Phylogenetic analysis supported that the nine D-hordeins of Leymus formed two independent clades, with all the 1.3 types clustered with Ps. juncea Ns 1.3, whereas the others were clustered together with the D-hordeins from Hordeum and Ps. juncea and the HMW-GSs from Leymus. Within the clade of four D-hordein genes and HMW-GSs, the HMW-GSs of Leymus formed a separated branch that served as an intermediate between the D-hordeins of Ps. juncea and Leymus. These novel D-hordeins may be potentially utilized in the improvement of food processing properties particularly those relating to extra cysteine residues. The findings of the present study also provide basic information for understanding the HMW-prolamins among Triticeae species, as well as expand the sources of D-hordeins from Hordeum to Leymus.

Comparing Multiple Reaction Monitoring and Sequential Window Acquisition of All Theoretical Mass Spectra for the Relative Quantification of Barley Gluten in Selectively Bred Barley Lines

Celiac disease (CD) is a disease of the small intestine that occurs in genetically susceptible subjects triggered by the ingestion of cereal gluten proteins for which the only treatment is strict adherence to a life-long gluten-free diet. Barley contains four gluten protein families, and the existence of barley genotypes that do not accumulate the B-, C-, and D-hordeins paved the way for the development of an ultralow gluten phenotype. Using conventional breeding strategies, three null mutations behaving as recessive alleles were combined to create a hordein triple-null barley variety. Proteomics has become an invaluable tool for characterization and quantification of the protein complement of cereal grains. In this study multiple reaction monitoring (MRM) mass spectrometry, viewed as the gold standard for peptide quantification, was compared to the data-independent acquisition strategy known as SWATH-MS (sequential window acquisition of all theoretical mass spectra). SWATH-MS was comparable (p < 0.001) to MRM-MS for 32/33 peptides assessed across the four families of hordeins (gluten) in eight barley lines. The results of SWATH-MS analysis further confirmed the absence of the B-, C-, and D-hordeins in the triple-null barley line and showed significantly reduced levels ranging from <1% to 16% relative to wild-type (WT) cv Sloop for the minor γ-hordein class. SWATH-MS represents a valuable tool for quantitative proteomics based on its ability to generate reproducible data comparable with MRM-MS, but has the added benefits of allowing reinterrogation of data to improve analytical performance, ask new questions, and in this case perform quantification of trypsin-resistant proteins (C-hordeins) through analysis of their semi- or nontryptic fragments.

Creation of the first ultra-low gluten barley (Hordeum vulgare L.) for coeliac and gluten-intolerant populations

Coeliac disease is a well-defined condition that is estimated to affect approximately 1% of the population worldwide. Noncoeliac gluten sensitivity is a condition that is less well defined, but is estimated to affect up to 10% of the population, and is often self-diagnosed. At present, the only remedy for both conditions is a lifelong gluten-free diet. A gluten-free diet is often expensive, high in fat and low in fibre, which in themselves can lead to adverse health outcomes. Thus, there is an opportunity to use novel plant breeding strategies to develop alternative gluten-free grains. In this work, we describe the breeding and characterization of a novel ultra-low gluten (ULG) barley variety in which the hordein (gluten) content was reduced to below 5 ppm. This was achieved using traditional breeding strategies to combine three recessive alleles, which act independently of each other to lower the hordein content in the parental varieties. The grain of the initial variety was shrunken compared to wild-type barleys. We implemented a breeding strategy to improve the grain size to near wild-type levels and demonstrated that the grains can be malted and brewed successfully. The ULG barley has the potential to provide novel healthy foods and beverages for those who require a gluten-free diet.

Transcriptome assembly and analysis of Tibetan Hulless Barley (Hordeum vulgare L. var. nudum) developing grains, with emphasis on quality properties

Background

Hulless barley is attracting increasing attention due to its unique nutritional value and potential health benefits. However, the molecular biology of the barley grain development and nutrient storage are not well understood. Furthermore, the genetic potential of hulless barley has not been fully tapped for breeding.

Methodology/Principal Findings

In the present study, we investigated the transcriptome features during hulless barley grain development. Using Illumina paired-end RNA-Sequencing, we generated two data sets of the developing grain transcriptomes from two hulless barley landraces. A total of 13.1 and 12.9 million paired-end reads with lengths of 90 bp were generated from the two varieties and were assembled to 48,863 and 45,788 unigenes, respectively. A combined dataset of 46,485 All-Unigenes were generated from two transcriptomes with an average length of 542 bp, and 36,278 among were annotated with gene descriptions, conserved protein domains or gene ontology terms. Furthermore, sequences and expression levels of genes related to the biosynthesis of storage reserve compounds (starch, protein, and β-glucan) were analyzed, and their temporal and spatial patterns were deduced from the transcriptome data of cultivated barley Morex.

Conclusions/Significance

We established a sequences and functional annotation integrated database and examined the expression profiles of the developing grains of Tibetan hulless barley. The characterization of genes encoding storage proteins and enzymes of starch synthesis and (1–3;1–4)-β-D-glucan synthesis provided an overview of changes in gene expression associated with grain nutrition and health properties. Furthermore, the characterization of these genes provides a gene reservoir, which helps in quality improvement of hulless barley.

The barley Frost resistance-H2 locus

Frost resistance-H2 (Fr-H2) is a major QTL affecting freezing tolerance in barley, yet its molecular basis is still not clearly understood. To gain a better insight into the structural characterization of the locus, a high-resolution linkage map developed from the Nure × Tremois cross was initially implemented to map 13 loci which divided the 0.602 cM total genetic distance into ten recombination segments. A PCR-based screening was then applied to identify positive bacterial artificial chromosome (BAC) clones from two genomic libraries of the reference genotype Morex. Twenty-six overlapping BACs from the integrated physical-genetic map were 454 sequenced. Reads assembled in contigs were subsequently ordered, aligned and manually curated in 42 scaffolds. In a total of 1.47 Mbp, 58 protein-coding sequences were identified, 33 of which classified according to similarity with sequences in public databases. As three complete barley C-repeat Binding Factors (HvCBF) genes were newly identified, the locus contained13 full-length HvCBFs, four Related to AP2 Triticeae (RAPT) genes, and at least five CBF pseudogenes. The final overall assembly of Fr-H2 includes more than 90 % of target region: all genes were identified along the locus, and a general survey of Repetitive Elements obtained. We believe that this gold-standard sequence for the Morex Fr-H2 will be a useful genomic tool for structural and evolutionary comparisons with Fr-H2 in winter-hardy cultivars along with Fr-2 of other Triticeae crops.

Characterization of x-type high-molecular-weight glutenin promoters (x-HGP) from different genomes in Triticeae

The sequences of x-type high-molecular-weight glutenin promoter (x-HGP) from 21 diploid Triticeae species were cloned and sequenced. The lengths of x-HGP varied from 897 to 955 bp, and there are 329 variable sites including 105 singleton sites and 224 polymorphic sites. Genetic distances of pairwise X-HGP sequences ranged from 0.30 to 16.40% within 21 species and four outgroup species of Hordeum. All five recognized regulatory elements emerged and showed higher conservation in the x-HGP of 21 Triticeae species. Most variations were distributed in the regions among or between regulatory elements. A 22 bp and 50 bp insertions which were the copy of adjacent region with minor change, were found in the x-HGP of Ae. speltoides and Ps. Huashanica, and could be regarded as genome specific indels. The phylogeny of media-joining network and neighbour-joining tree both supported the topology were composed of three sperate clusters. Especially, the cluster I comprising the x-HGP sequences of Aegilops, Triticum, Henrardia, Agropyron and Taeniatherum was highly supporting by both network and NJ tree. As conferring to higher level and temporal and spatial expression, x-HGP can used as the source of promoter for constructing transgenic plants which allow endosperm-specific expression of exogenous gene on higher level. In addition, the x-HGP has enough conservation and variation; so it should be valuable in phylogenetic analyses of Triticeae family members.

Dynamic evolution of rht-1 homologous regions in grass genomes

Hexaploid bread wheat contains A, B, and D three subgenomes with its well-characterized ancestral genomes existed at diploid and tetraploid levels, making the wheat act as a good model species for studying evolutionary genomic dynamics. Here, we performed intra- and inter-species comparative analyses of wheat and related grass genomes to examine the dynamics of homologous regions surrounding Rht-1, a well-known "green revolution" gene. Our results showed that the divergence of the two A genomes in the Rht-1 region from the diploid and tetraploid species is greater than that from the tetraploid and hexaploid wheat. The divergence of D genome between diploid and hexaploid is lower than those of A genome, suggesting that D genome diverged latter than others. The divergence among the A, B and D subgenomes was larger than that among different ploidy levels for each subgenome which mainly resulted from genomic structural variation of insertions and, perhaps deletions, of the repetitive sequences. Meanwhile, the repetitive sequences caused genome expansion further after the divergence of the three subgenomes. However, several conserved non-coding sequences were identified to be shared among the three subgenomes of wheat, suggesting that they may have played an important role to maintain the homolog of three subgenomes. This is a pilot study on evolutionary dynamics across the wheat ploids, subgenomes and differently related grasses. Our results gained new insights into evolutionary dynamics of Rht-1 region at sequence level as well as the evolution of wheat during the plolyploidization process.

The MS(E)-proteomic analysis of gliadins and glutenins in wheat grain identifies and quantifies proteins associated with celiac disease and baker's asthma

Precise content of gliadin (Glia) and glutenin (Glu) proteins in wheat grain are largely unknown despite their association with celiac disease, various allergies, and physical processing properties of wheat. Developing methods to quantitatively measure clinically relevant proteins could support advancement in understanding exposure thresholds and clinical study design. The aim of this study was to use a data-independent mass spectrometry (MS(E)) approach for quantifying gliadin and glutenin proteins in wheat grain. The biologically replicated analysis yielded concentrations for 34 gliadin and 22 glutenin proteins. The primary focus of this survey was on measuring celiac disease proteins and baker's asthma associated proteins along with the proteins associated with viscoelastic properties of wheat flour and grain texture. The technical coefficients of variation ranged from 0.12 to 1.39 and indicate that MS(E) proteomics is a reproducible quantitative method for the determination of gliadin and glutenin content in the highly complex matrix of protein extracts from wheat grain. This article is part of a Special Issue entitled: Translational Plant Proteomics.

Quantitative RT-PCR based platform for rapid quantification of the transcripts of highly homologous multigene families and their members during grain development

BACKGROUND

Cereal storage proteins represent one of the most important sources of protein for food and feed and they are coded by multigene families. The expression of the storage protein genes exhibits a temporal fluctuation but also a response to environmental stimuli. Analysis of temporal gene expression combined with genetic variation in large multigene families with high homology among the alleles is very challenging.

RESULTS

We designed a rapid qRT-PCR system with the aim of characterising the variation in the expression of hordein genes families. All the known D-, C-, B-, and γ-hordein sequences coding full length open reading frames were collected from commonly available databases. Phylogenetic analysis was performed and the members of the different hordein families were classified into subfamilies. Primer sets were designed to discriminate the gene expression level of whole families, subfamilies or individual members. The specificity of the primer sets was validated before successfully applying them to a cDNA population derived from developing grains of field grown Hordeum vulgare cv. Barke. The results quantify the number of moles of transcript contributed to a particular gene family and its subgroups. More over the results indicate the genotypic specific gene expression.

CONCLUSIONS

Quantitative RT-PCR with SYBR Green labelling can be a useful technique to follow gene expression levels of large gene families with highly homologues members. We showed variation in the temporal expression of genes coding for barley storage proteins. The results imply that our rapid qRT-PCR system was sensitive enough to identify the presence of alleles and their expression profiles. It can be used to check the temporal fluctuations in hordein expressions or to find differences in their response to environmental stimuli. The method could be extended for cultivar recognition as some of the sequences from the database originated from cv. Golden Promise were not expressed in the studied barley cultivar Barke although showed primer specificity with their cloned DNA sequences.

Mitotic illegitimate recombination is a mechanism for novel changes in high-molecular-weight glutenin subunits in wheat-rye hybrids

Wide hybrids can have novel traits or changed expression of a quantitative trait that their parents do not have. These phenomena have long been noticed, yet the mechanisms are poorly understood. High-molecular-weight glutenin subunits (HMW-GS) are seed storage proteins encoded by Glu-1 genes that only express in endosperm in wheat and its related species. Novel HMW-GS compositions have been observed in their hybrids. This research elucidated the molecular mechanisms by investigating the causative factors of novel HMW-GS changes in wheat-rye hybrids. HMW-GS compositions in the endosperm and their coding sequences in the leaves of F₁ and F₂ hybrids between wheat landrace Shinchunaga and rye landrace Qinling were investigated. Missing and/or additional novel HMW-GSs were observed in the endosperm of 0.5% of the 2078 F₁ and 22% of 36 F₂ hybrid seeds. The wildtype Glu-1Ax null allele was found to have 42 types of short repeat sequences of 3-60 bp long that appeared 2 to 100 times. It also has an in-frame stop codon in the central repetitive region. Analyzing cloned allele sequences of HMW-GS coding gene Glu-1 revealed that deletions involving the in-frame stop codon had happened, resulting in novel ∼1.8-kb Glu-1Ax alleles in some F₁ and F₂ plants. The cloned mutant Glu-1Ax alleles were expressed in Escherichia coli, and the HMW-GSs produced matched the novel HMW-GSs found in the hybrids. The differential changes between the endosperm and the plant of the same hybrids and the data of E. coli expression of the cloned deletion alleles both suggested that mitotic illegitimate recombination between two copies of a short repeat sequence had resulted in the deletions and thus the changed HMW-GS compositions. Our experiments have provided the first direct evidence to show that mitotic illegitimate recombination is a mechanism that produces novel phenotypes in wide hybrids.

Structural variation and evolutionary relationship of novel HMW glutenin subunits from Elymus glaucus

High molecular weight (HMW) glutenin subunits (GS) are important seed storage proteins relevant to the end-use quality of wheat and other cereal crops. Here we report the isolation and characterization of two novel HMW-GS alleles (1St 1.4 and 1St1.1) from the perennial Triticeae species Elymus glaucus. The amino acid (aa) sequences of E. glaucus 1St1.4 and 1St1.1 were predicted as 434 aa and 358 aa, respectively. Both subunits comprise a signal peptide with a conserved N-terminal domain, a central repetitive domain and a C-terminal domain. Elymus glaucus 1St 1.4 and 1St1.1 exhibit several distinct characteristics different from other known HMW-GSs. The lengths of repetitive domains in E. glaucus 1St 1.4 and 1St1.1 are substantially smaller than those of other known HMW-GSs, in which 1St1.1 (only 358 aa) is the smallest subunit identified so far. The N-terminal domains of E. glaucus 1St 1.4 and 1St1.1 are homologous to y-type subunits, whereas their C-terminal domains are similar to x-type subunits. Our results indicate that E. glaucus 1St 1.4 and 1St1.1 are novel HMW-GS variants or isoforms, and the characterization of both subunits can enhance our understanding on the structural differentiation and evolutionary relationship of HMW-GSs in Triticeae species.

Conserved globulin gene across eight grass genomes identify fundamental units of the loci encoding seed storage proteins

The wheat high molecular weight (HMW) glutenins are important seed storage proteins that determine bread-making quality in hexaploid wheat (Triticum aestivum). In this study, detailed comparative sequence analyses of large orthologous HMW glutenin genomic regions from eight grass species, representing a wide evolutionary history of grass genomes, reveal a number of lineage-specific sequence changes. These lineage-specific changes, which resulted in duplications, insertions, and deletions of genes, are the major forces disrupting gene colinearity among grass genomes. Our results indicate that the presence of the HMW glutenin gene in Triticeae genomes was caused by lineage-specific duplication of a globulin gene. This tandem duplication event is shared by Brachypodium and Triticeae genomes, but is absent in rice, maize, and sorghum, suggesting the duplication occurred after Brachypodium and Triticeae genomes diverged from the other grasses ~35 Ma ago. Aside from their physical location in tandem, the sequence similarity, expression pattern, and conserved cis-acting elements responsible for endosperm-specific expression further support the paralogous relationship between the HMW glutenin and globulin genes. While the duplicated copy in Brachypodium has apparently become nonfunctional, the duplicated copy in wheat has evolved to become the HMW glutenin gene by gaining a central prolamin repetitive domain.

EST mining for structure and expression of genes in the region of the wheat high-molecular-weight glutenin loci

An in-depth analysis was carried out with expressed sequence tags (ESTs) for genes in and near the HMW-GS loci. Considerations for using ESTs are discussed, including the occurrence of chimeric and aberrant HMW-GS ESTs. Complete gene sequences demonstrated the feasibility of constructing accurate full-length coding regions from EST assemblies and found, or supported, errors in several previously reported HMW-GS gene sequences. New complete HMW-GS gene sequences are reported for the cultivars Chinese Spring and Glenlea. The Ay subunit gene, which is considered null in cultivated wheats, was shown to transcribe in at least two germplasms. Analyses support the conclusion that of the five known genes within this genomic region, the two HMW-GS genes and the globulin gene are highly expressed. The other two genes, encoding a receptor kinase and a protein kinase, have one and no identifiable wheat EST, respectively, although ESTs are found for the orthologous genes in barley. The ESTs of all five genes within the HMW-GS region are either definitely associated with the endosperm or possibly originate from imbibed seed, suggesting the four distinct gene classes in this region are part of a seed or endosperm chromatin domain. EST resources were also used to determine relative abundance of ESTs for all classes of wheat prolamines and indicated differential levels of expression both among germplasms and among the three genomes of hexaploid wheats.

Structural characterization of Brachypodium genome and its syntenic relationship with rice and wheat

Brachypodium distachyon (Brachypodium) has been recently recognized as an emerging model system for both comparative and functional genomics in grass species. In this study, 55,221 repeat masked Brachypodium BAC end sequences (BES) were used for comparative analysis against the 12 rice pseudomolecules. The analysis revealed that approximately 26.4% of BES have significant matches with the rice genome and 82.4% of the matches were homologous to known genes. Further analysis of paired-end BES and approximately 1.0 Mb sequences from nine selected BACs proved to be useful in revealing conserved regions and regions that have undergone considerable genomic changes. Differential gene amplification, insertions/deletions and inversions appeared to be the common evolutionary events that caused variations of microcolinearity at different orthologous genomic regions. It was found that approximately 17% of genes in the two genomes are not colinear in the orthologous regions. Analysis of BAC sequences also revealed higher gene density (approximately 9 kb/gene) and lower repeat DNA content (approximately 13.1%) in Brachypodium when compared to the orthologous rice regions, consistent with the smaller size of the Brachypodium genome. The 119 annotated Brachypodium genes were BLASTN compared against the wheat EST database and deletion bin mapped wheat ESTs. About 77% of the genes retrieved significant matches in the EST database, while 9.2% matched to the bin mapped ESTs. In some cases, genes in single Brachypodium BACs matched to multiple ESTs that were mapped to the same deletion bins, suggesting that the Brachypodium genome will be useful for ordering wheat ESTs within the deletion bins and developing specific markers at targeted regions in the wheat genome.

Low-pass shotgun sequencing of the barley genome facilitates rapid identification of genes, conserved non-coding sequences and novel repeats

Background

Barley has one of the largest and most complex genomes of all economically important food crops. The rise of new short read sequencing technologies such as Illumina/Solexa permits such large genomes to be effectively sampled at relatively low cost. Based on the corresponding sequence reads a Mathematically Defined Repeat (MDR) index can be generated to map repetitive regions in genomic sequences.

Results

We have generated 574 Mbp of Illumina/Solexa sequences from barley total genomic DNA, representing about 10% of a genome equivalent. From these sequences we generated an MDR index which was then used to identify and mark repetitive regions in the barley genome. Comparison of the MDR plots with expert repeat annotation drawing on the information already available for known repetitive elements revealed a significant correspondence between the two methods. MDR-based annotation allowed for the identification of dozens of novel repeat sequences, though, which were not recognised by hand-annotation. The MDR data was also used to identify gene-containing regions by masking of repetitive sequences in eight de-novo sequenced bacterial artificial chromosome (BAC) clones. For half of the identified candidate gene islands indeed gene sequences could be identified. MDR data were only of limited use, when mapped on genomic sequences from the closely related species Triticum monococcum as only a fraction of the repetitive sequences was recognised.

Conclusion

An MDR index for barley, which was obtained by whole-genome Illumina/Solexa sequencing, proved as efficient in repeat identification as manual expert annotation. Circumventing the labour-intensive step of producing a specific repeat library for expert annotation, an MDR index provides an elegant and efficient resource for the identification of repetitive and low-copy (i.e. potentially gene-containing sequences) regions in uncharacterised genomic sequences. The restriction that a particular MDR index can not be used across species is outweighed by the low costs of Illumina/Solexa sequencing which makes any chosen genome accessible for whole-genome sequence sampling.

D hordeins of Hordeum chilense: a novel source of variation for improvement of wheat

The high molecular weight subunits of wheat (Triticum aestivum L.) glutenin (HMW-GS) are important in determining the bread-making quality of flour and dough. There is therefore interest in transferring orthologous HMW-GS present in other grass species into wheat by wide crossing in order to extend the range of end use properties. In this work, we have isolated and characterized two genes encoding D hordeins from Hordeum chilense (Roem. et Schult.) lines H1 and H7, representing two ecotypes. The fragments were 4,305 bp for line H1 and 4,227 for line H7 and contained the promoter, coding and terminator regions. Both sequences differ in the presence of single base changes (SNPs) and insertions/deletions in the open reading frame (ORF). The encoded proteins comprise 870 and 896 amino acids for lines H1 and H7, respectively. The primary structure is similar to those of D hordeins of cultivated barley (H. vulgare L.) and HMW-GS from wheat. However, the D hordeins from H. chilense are significantly larger than those from cultivated barley due to the presence of longer repetitive regions. The H. chilense D hordeins also differ from those of cultivated barley in the distribution of the cysteine residues: whereas the D hordeins of cultivated barley contain ten cysteines with four in the repetitive domain, only nine are present in the H. chilense proteins with two in the repetitive domain. As in the HMW-GS, the central part of the D hordein proteins comprises repeated sequences based on short peptide motifs. The repetitive domain is divided in three regions named as R1 (N-terminal repeats), R2 (central degenerate repeats) and R3 (C-terminal repeats). Hexapeptide motifs are present throughout the repetitive domains of D hordeins with a consensus motif of PFQGQQ in R1 and R2 and PHQGQQ in R3. In addition, the tetrapeptide motif TTVS, which is characteristic of D hordeins of cultivated barley is present in the repetitive domain close to the protein C-terminus.

Triticeae genomics: advances in sequence analysis of large genome cereal crops

Whole genome sequencing provides direct access to all genes of an organism and represents an essential step towards a systematic understanding of (crop) plant biology. Wheat and barley, two of the most important crop species worldwide, have two- to five-fold larger genomes than human - too large to be completely sequenced at current costs. Nevertheless, significant progress has been made to unlock the gene contents of these species by sequencing expressed sequence tags (EST) for high-density mapping and as a basis for elucidating gene function on a large scale. Several megabases of genomic (BAC) sequences have been obtained providing a first insight into the complexity of these huge cereal genomes. However, to fully exploit the information of the wheat and barley genomes for crop improvement, sequence analysis of a significantly larger portion of the Triticeae genomes is needed. In this review an overview of the current status of Triticeae genome sequencing and a perspective concerning future developments in cereal structural genomics is provided.

Characterization of HMW prolamines and their coding sequences from Crithopsis delileana

The high-molecular-weight (HMW) prolamines subunits and their coding sequences from wheat-related diploid species Crithopsis delileana were investigated. Only one HMW prolamine subunit with the similar electrophoresis mobility to the y-type HMW glutenin subunit of hexaploid wheat was observed in two accessions of C. delileana by SDS-PAGE analyses of the total storage protein fractions. It was confirmed by sequencing and expression analysis that this prolamine subunit was an x-type subunit. The amino acid sequence of this subunit had the similar typical structure to those of x-type HMW glutenin genes previously described in wheat. An in-frame stop codon was found in the coding sequences of y-type prolamine subunits. It was found by specifically extraction of HMW prolamines and sequence analysis that the coding regions of Ky prolamine subunit gene is very likely to be not expressed as a full-length protein. Phylogenetic analysis indicated that the Kx subunit could be clustered together with 1Ax1 subunit by an interior paralleled branch, and Ky subunit (inactive) was most closely related to the 1Ay subunit. The coding sequences of Kx subunit could successfully be expressed in bacterial expression system, and the expressed protein had the same electrophoresis mobility as the Kx subunit from the seed of C. delileana. It was the first time that the HMW prolamines subunits encoded by K genome of C. delileana were characterized.

Types and rates of sequence evolution at the high-molecular-weight glutenin locus in hexaploid wheat and its ancestral genomes

The Glu-1 locus, encoding the high-molecular-weight glutenin protein subunits, controls bread-making quality in hexaploid wheat (Triticum aestivum) and represents a recently evolved region unique to Triticeae genomes. To understand the molecular evolution of this locus region, three orthologous Glu-1 regions from the three subgenomes of a single hexaploid wheat species were sequenced, totaling 729 kb of sequence. Comparing each Glu-1 region with its corresponding homologous region from the D genome of diploid wheat, Aegilops tauschii, and the A and B genomes of tetraploid wheat, Triticum turgidum, revealed that, in addition to the conservation of microsynteny in the genic regions, sequences in the intergenic regions, composed of blocks of nested retroelements, are also generally conserved, although a few nonshared retroelements that differentiate the homologous Glu-1 regions were detected in each pair of the A and D genomes. Analysis of the indel frequency and the rate of nucleotide substitution, which represent the most frequent types of sequence changes in the Glu-1 regions, demonstrated that the two A genomes are significantly more divergent than the two B genomes, further supporting the hypothesis that hexaploid wheat may have more than one tetraploid ancestor.

Cloning and molecular characterization of three novel LMW-i glutenin subunit genes from cultivated einkorn (Triticum monococcum L.)

Three novel low molecular weight (LMW) glutenin subunits from cultivated einkorn (Triticum monococcum L., A(m)A(m), 2n = 2x = 14) were characterized by SDS-PAGE and molecular weights determined by MALDI-TOF-MS. Their coding genes were amplified and cloned with designed AS-PCR primers, revealing three complete gene sequences. All comprised upstream, open reading frame (ORF), downstream and no introns were present. The deduced amino acid sequences showed that all three genes, named as LMW-M1, LMW-M3 and LMW-M5, respectively, belonged to the LMW-i type subunits with the predicted molecular weight between 38.5206 and 38.7028 kDa. They showed high similarity with other LMW-i type genes from hexaploid bread wheats, but also displayed unique features. Particularly, LMW-M5 subunit contained an extra cysteine residue in the C-terminus except for eight conserved cysteines, which resulted from a single-nucleotide polymorphism (SNP) of the T-C transition, namely arginine --> cysteine substitution at position 242 from the N-terminal end. This is the first report that the LMW-i subunit contained nine cysteines residues that could result in a more highly cross-linked and more elastic glutenin suggesting that LMW-M5 gene may associates with good quality properties. In addition, a total of 25 SNPs and one insertions/deletions (InDels) were detected among three LMW-i genes, which could result in significant functional changes in polymer formation of gluten. It is anticipated that these SNPs could be used as reliable genetic markers during wheat quality improvement. The phylogenetic analysis indicated that LMW-i type genes apparently differed from LMW-m and LMW-s type genes and diverged early from the primitive LMW-GS gene family, at about 12.92 million years ago (MYA) while the differentiation of A(m) and A genomes was estimated at 3.98 MYA.

Quantification of the retrotransposon BARE-1 reveals the dynamic nature of the barley genome

We used quantitative real-time PCR analysis to measure the copy number of the BARE-1 retrotransposon in 5 cultivars of barley (Hordeum vulgare), as well as in samples from its wild relative, Hordeum spontaneum. Two sets of PCR primers were used to amplify regions within the long terminal repeat (LTR) and the reverse transcriptase (RT) gene of BARE-1 (GenBank accession Z17327). The LTR primers detected an average of 2.148 x 105 +/- 0.012 x 105 copies per haploid genome among barley samples, whereas the RT primers detected an average of 1.588 x 104 +/- 0.085 x 104 copies. The average ratio of LTR:RT was estimated to be 13.5:1. This finding indicates that more than 7% of the barley genome is occupied by BARE-1 elements in the form of solo LTRs and another 2.6% of the genome is occupied by the full-length element. Taken together, BARE-1 sequences represent approximately 9.6% of the barley genome among the barley plants used in this study. For the above estimation, a genome size of 5.44 x 103 Mb for H. vulgare and 5.39 x 103 Mb for H. spontaneum were assumed. Our study on quantification results of the BARE-1 for a small group of barley cultivars showed that there are significant differences among cultivars in terms of BARE-1 copy number, providing further evidence that BARE-1 is active and has a major role in shaping the barley genome as a result of breeding and selection. Quantification results also showed that most of the elements (> 90%) are present as truncated copies (solo LTRs). These results show that there is a high level of recombination leading to the formation of truncated elements and a subsequent DNA loss from the genome. Taken together, our study provides a glimpse into a dynamic micro-evolutionary process that is the by-product of genome reshuffling and directional selection in barley breeding programs.

Ancient haplotypes resulting from extensive molecular rearrangements in the wheat A genome have been maintained in species of three different ploidy levels

Plant genomes, in particular grass genomes, evolve very rapidly. The closely related A genomes of diploid, tetraploid, and hexaploid wheat are derived from a common ancestor that lived <3 million years ago and represent a good model to study molecular mechanisms involved in such rapid evolution. We have sequenced and compared physical contigs at the Lr10 locus on chromosome 1AS from diploid (211 kb), tetraploid (187 kb), and hexaploid wheat (154 kb). A maximum of 33% of the sequences were conserved between two species. The sequences from diploid and tetraploid wheat shared all of the genes, including Lr10 and RGA2 and define a first haplotype (H1). The 130-kb intergenic region between Lr10 and RGA2 was conserved in size despite its activity as a hot spot for transposon insertion, which resulted in >70% of sequence divergence. The hexaploid wheat sequence lacks both Lr10 and RGA2 genes and defines a second haplotype, H2, which originated from ancient and extensive rearrangements. These rearrangements included insertions of retroelements and transposons deletions, as well as unequal recombination within elements. Gene disruption in haplotype H2 was caused by a deletion and subsequent large inversion. Gene conservation between H1 haplotypes, as well as conservation of rearrangements at the origin of the H2 haplotype at three different ploidy levels indicate that the two haplotypes are ancient and had a stable gene content during evolution, whereas the intergenic regions evolved rapidly. Polyploidization during wheat evolution had no detectable consequences on the structure and evolution of the two haplotypes.

Large intraspecific haplotype variability at the Rph7 locus results from rapid and recent divergence in the barley genome

To study genome evolution and diversity in barley (Hordeum vulgare), we have sequenced and compared more than 300 kb of sequence spanning the Rph7 leaf rust disease resistance gene in two barley cultivars. Colinearity was restricted to five genic and two intergenic regions representing <35% of the two sequences. In each interval separating the seven conserved regions, the number and type of repetitive elements were completely different between the two homologous sequences, and a single gene was absent in one cultivar. In both cultivars, the nonconserved regions consisted of approximately 53% repetitive sequences mainly represented by long-terminal repeat retrotransposons that have inserted <1 million years ago. PCR-based analysis of intergenic regions at the Rph7 locus and at three other independent loci in 41 H. vulgare lines indicated large haplotype variability in the cultivated barley gene pool. Together, our data indicate rapid and recent divergence at homologous loci in the genome of H. vulgare, possibly providing the molecular mechanism for the generation of high diversity in the barley gene pool. Finally, comparative analysis of the gene composition in barley, wheat (Triticum aestivum), rice (Oryza sativa), and sorghum (Sorghum bicolor) suggested massive gene movements at the Rph7 locus in the Triticeae lineage.

Comparative sequence analysis of the region harboring the hardness locus in barley and its colinear region in rice

The ancestral shared synteny concept has been advocated as an approach to positionally clone genes from complex genomes. However, the unified grass genome model and the study of grasses as a single syntenic genome is a topic of considerable controversy. Hence, more quantitative studies of cereal colinearity at the sequence level are required. This study compared a contiguous 300-kb sequence of the barley (Hordeum vulgare) genome with the colinear region in rice (Oryza sativa). The barley sequence harbors genes involved in endosperm texture, which may be the subject of distinctive evolutionary forces and is located at the extreme telomeric end of the short arm of chromosome 5H. Comparative sequence analysis revealed the presence of five orthologous genes and a complex, postspeciation evolutionary history involving small chromosomal rearrangements, a translocation, numerous gene duplications, and extensive transposon insertion. Discrepancies in gene content and microcolinearity indicate that caution should be exercised in the use of rice as a surrogate for map-based cloning of genes from large genome cereals such as barley.

Genomic organization of the complex alpha-gliadin gene loci in wheat

To better understand the molecular evolution of the large alpha-gliadin gene family, a half-million bacterial artificial chromosome (BAC) library clones from tetraploid durum wheat, Triticum turgidum ssp. durum (2n = 4x = 28, genome AB), were screened for large genomic segments carrying the alpha-gliadin genes of the Gli-2 loci on the group 6 homoeologous chromosomes. The resulting 220 positive BAC clones--each containing between one and four copies of alpha-gliadin sequences--were fingerprinted for contig assembly to produce contiguous chromosomal regions covering the Gli-2 loci. While contigs consisting of as many as 21 BAC clones and containing up to 17 alpha-gliadin genes were formed, many BAC clones remained as singletons. The accuracy of the order of BAC clones in the contigs was verified by Southern hybridization analysis of the BAC fingerprints using an alpha-gliadin probe. These results indicate that alpha-gliadin genes are not evenly dispersed in the Gli-2 locus regions. Hybridization of these BACs with probes for long terminal repeat retrotransposons was used to determine the abundance and distribution of repetitive DNA in this region. Sequencing of BAC ends indicated that 70% of the sequences were significantly similar to different classes of retrotransposons, suggesting that these elements are abundant in this region. Several mechanisms underlying the dynamic evolution of the Gli-2 loci are discussed.

Rapid genome evolution revealed by comparative sequence analysis of orthologous regions from four triticeae genomes

Bread wheat (Triticum aestivum) is an allohexaploid species, consisting of three subgenomes (A, B, and D). To study the molecular evolution of these closely related genomes, we compared the sequence of a 307-kb physical contig covering the high molecular weight (HMW)-glutenin locus from the A genome of durum wheat (Triticum turgidum, AABB) with the orthologous regions from the B genome of the same wheat and the D genome of the diploid wheat Aegilops tauschii (Anderson et al., 2003; Kong et al., 2004). Although gene colinearity appears to be retained, four out of six genes including the two paralogous HMW-glutenin genes are disrupted in the orthologous region of the A genome. Mechanisms involved in gene disruption in the A genome include retroelement insertions, sequence deletions, and mutations causing in-frame stop codons in the coding sequences. Comparative sequence analysis also revealed that sequences in the colinear intergenic regions of these different genomes were generally not conserved. The rapid genome evolution in these regions is attributable mainly to the large number of retrotransposon insertions that occurred after the divergence of the three wheat genomes. Our comparative studies indicate that the B genome diverged prior to the separation of the A and D genomes. Furthermore, sequence comparison of two distinct types of allelic variations at the HMW-glutenin loci in the A genomes of different hexaploid wheat cultivars with the A genome locus of durum wheat indicates that hexaploid wheat may have more than one tetraploid ancestor.

Sequencing of the Triticum monococcum Hardness locus reveals good microcolinearity with rice

The Hardness ( Ha) locus on chromosome 5D is the main determinant of grain texture in hexaploid wheat. The related genes Puroindoline-a and -b ( Pina-D1, Pinb-D1) and Grain Softness Protein ( Gsp-D1) are tightly linked at this locus. Mutations in the Pina-D1 and Pinb-D1 genes are associated with increased grain hardness. We report here the complete sequence of a 101-kb BAC clone from Triticum monococcum (A(m ) genome) which includes these three genes, and its comparison with the orthologous region in rice. The genes Gsp-A(m) 1, Pina-A(m) 1 and Pinb-A(m) 1 are separated by 37 kb and 32 kb, respectively, and are organized in the same transcriptional orientation. Four additional genes, including a pair of duplicated genes, were identified upstream of Gsp-A(m) 1 within a high-density gene island. These additional genes were found in the same order and orientation, and the same relative distances apart as similar genes previously annotated on rice chromosome 12. An interesting discovery was a small unannotated putative rice gene that was similar to the Gsp-A(m) 1 gene of T. monococcum (65% similarity at the protein level), and that was disposed in the same orientation, and located in the same position relative to the other orthologous genes. The high gene density observed in this BAC (1 gene per 14 kb) was expected for a distal chromosome region, but the level of microcolinearity with rice was higher than that reported in similar distal regions of other wheat chromosomes. Most of the BAC sequence (40%) was represented by repetitive elements, mainly concentrated in regions adjacent to the genes Pina-A(m) 1 and Pinb-A(m) 1. Rearrangements among these repetitive elements might provide an explanation for the frequent deletions observed at this locus in the genomes of the polyploid wheat species.