Sequencing wheat chromosome arm 7BS delimits the 7BS/4AL translocation and reveals homoeologous gene conservation

Complex Triticeae genomes pose a challenge to genome sequencing efforts due to their size and repetitive nature. Genome sequencing can reveal details of conservation and rearrangements between related genomes. We have applied Illumina second generation sequencing technology to sequence and assemble the low copy and unique regions of Triticum aestivum chromosome arm 7BS, followed by the construction of a syntenic build based on gene order in Brachypodium. We have delimited the position of a previously reported translocation between 7BS and 4AL with a resolution of one or a few genes and report approximately 13% genes from 7BS having been translocated to 4AL. An additional 13 genes are found on 7BS which appear to have originated from 4AL. The gene content of the 7DS and 7BS syntenic builds indicate a total of ~77,000 genes in wheat. Within wheat syntenic regions, 7BS and 7DS share 740 genes and a common gene conservation rate of ~39% of the genes from the corresponding regions in Brachypodium, as well as a common rate of colinearity with Brachypodium of ~60%. Comparison of wheat homoeologues revealed ~84% of genes previously identified in 7DS have a homoeologue on 7BS or 4AL. The conservation rates we have identified among wheat homoeologues and with Brachypodium provide a benchmark of homoeologous gene conservation levels for future comparative genomic analysis. The syntenic build of 7BS is publicly available at http://www.wheatgenome.info.     

Exploiting co-linearity among grass species to map the <Emphasis Type="Italic">Aegilops ventricosa-</Emphasis>derived <Emphasis Type="Italic">Pch1</Emphasis> eyespot resistance in wheat and establish its relationship to <Emphasis Type="Italic">Pch2</Emphasis>

Introgressions into wheat from related species have been widely used as a source of agronomically beneficial traits. One such example is the introduction of the potent eyespot resistance gene Pch1 from the wild relative Aegilops ventricosa onto chromosome 7DL of wheat. In common with genes carried on many other such introgressions, the use of Pch1 in commercial wheat varieties has been hindered by linkage drag with yield-limiting traits. Attempts to break this linkage have been frustrated by a lack of co-dominant PCR markers suitable for identifying heterozygotes in F₂ populations. We developed conserved orthologous sequence (COS) markers, utilising the Brachypodium distachyon (Brachypodium) genome sequence, to provide co-dominant markers in the Pch1 region. These were supplemented with previously developed sequence-tagged site (STS) markers and simple sequence repeat (SSR) markers. Markers were applied to a panel of varieties and to a BC₆ F₂ population, segregating between wheat and Ae. ventricosa over the distal portion of 7DL, to identify recombinants in the region of Pch1. By exploiting co-linearity between wheat chromosome 7D, Brachypodium chromosome 1, rice chromosome 6 and sorghum chromosome 10, Pch1 was located to an interval between the flanking markers Xwg7S and Xcos7-9. Furthermore candidate gene regions were identified in Brachypodium (364 Kb), rice (178 Kb) and sorghum (315 Kb) as a prelude to the map-based cloning of the gene. In addition, using homoeologue transferable markers, we obtained evidence that the eyespot resistances Pch1 and Pch2 on chromosomes 7D and 7A, respectively, are potentially homoeoloci. It is anticipated that the COS marker methodology could be used for the identification of recombinants in other introgressions into wheat from wild relatives. This would assist the mapping of genes of interest and the breaking of deleterious linkages to enable greater use of these introgressions in commercial varieties.     

The epitopes in wheat proteins for defining toxic units relevant to human health

Wheat-related disorders are well-studied health problems. Knowledge of the composition and amounts of epitopes present in a single wheat sample represents a significant gap, and the detailed wheat proteome datasets now available can provide the necessary information to carry out an estimation of allergen prediction for a single cultivar. The combined use of genome sequence and allergen databases, prediction methodology, and cereal chemistry results in better understanding of the level of toxicity present in the end-products produced from wheat flour. The workflow presented in this review provides information about the number and distribution of epitopes at single protein, or protein fraction, levels. In addition, epitopes present in the highest frequency and harmful proteins expressed in the highest amount can be identified. The “epitope toxicity” value obtained in this way is a significant research output from the analysis of large datasets that can be applied to the food industry.     

Construction and characterization of two BAC libraries from Brachypodium distachyon, a new model for grass genomics.

Brachypodium is well suited as a model system for temperate grasses because of its compact genome and a range of biological features. In an effort to develop resources for genome research in this emerging model species, we constructed 2 bacterial artificial chromosome (BAC) libraries from an inbred diploid Brachypodium distachyon line, Bd21, using restriction enzymes HindIII and BamHI. A total of 73,728 clones (36,864 per BAC library) were picked and arrayed in 192,384-well plates. The average insert size for the BamHI and HindIII libraries is estimated to be 100 and 105 kb, respectively, and inserts of chloroplast origin account for 4.4% and 2.4%, respectively. The libraries individually represent 9.4- and 9.9-fold haploid genome equivalents with combined 19.3-fold genome coverage, based on a genome size of 355 Mb reported for the diploid Brachypodium, implying a 99.99% probability that any given specific sequence will be present in each library. Hybridization of the libraries with 8 starch biosynthesis genes was used to empirically evaluate this theoretical genome coverage; the frequency at which these genes were present in the library clones gave an estimated coverage of 11.6- and 19.6-fold genome equivalents. To obtain a first view of the sequence composition of the Brachypodium genome, 2185 BAC end sequences (BES) representing 1.3 Mb of random genomic sequence were compared with the NCBI GenBank database and the GIRI repeat database. Using a cutoff expectation value of E<10-10, only 3.3% of the BESs showed similarity to repetitive sequences in the existing database, whereas 40.0% had matches to the sequences in the EST database, suggesting that a considerable portion of the Brachypodium genome is likely transcribed. When the BESs were compared with individual EST databases, more matches hit wheat than maize, although their EST collections are of a similar size, further supporting the close relationship between Brachypodium and the Triticeae. Moreover, 122 BESs have significant matches to wheat ESTs mapped to individual chromosome bin positions. These BACs represent colinear regions containing the mapped wheat ESTs and would be useful in identifying additional markers for specific wheat chromosome regions.     

Comparison of a high-density genetic linkage map to genome features in the model grass <Emphasis Type="Italic">Brachypodium distachyon</Emphasis>

The small annual grass Brachypodium distachyon (Brachypodium) is rapidly emerging as a powerful model system to study questions unique to the grasses. Many Brachypodium resources have been developed including a whole genome sequence, highly efficient transformation and a large germplasm collection. We developed a genetic linkage map of Brachypodium using single nucleotide polymorphism (SNP) markers and an F₂ mapping population of 476 individuals. SNPs were identified by targeted resequencing of single copy genomic sequences. Using the Illumina GoldenGate Genotyping platform we placed 558 markers into five linkage groups corresponding to the five chromosomes of Brachypodium. The unusually long total genetic map length, 1,598 centiMorgans (cM), indicates that the Brachypodium mapping population has a high recombination rate. By comparing the genetic map to genome features we found that the recombination rate was positively correlated with gene density and negatively correlated with repetitive regions and sites of ancestral chromosome fusions that retained centromeric repeat sequences. A comparison of adjacent genome regions with high versus low recombination rates revealed a positive correlation between interspecific synteny and recombination rate.     

Using genic sequence capture in combination with a syntenic pseudo genome to map a deletion mutant in a wheat species

Mapping‐by‐sequencing analyses have largely required a complete reference sequence and employed whole genome re‐sequencing. In species such as wheat, no finished genome reference sequence is available. Additionally, because of its large genome size (17 Gb), re‐sequencing at sufficient depth of coverage is not practical. Here, we extend the utility of mapping by sequencing, developing a bespoke pipeline and algorithm to map an early‐flowering locus in einkorn wheat (Triticum monococcum L.) that is closely related to the bread wheat genome A progenitor. We have developed a genomic enrichment approach using the gene‐rich regions of hexaploid bread wheat to design a 110‐Mbp NimbleGen SeqCap EZ in solution capture probe set, representing the majority of genes in wheat. Here, we use the capture probe set to enrich and sequence an F₂ mapping population of the mutant. The mutant locus was identified in T. monococcum, which lacks a complete genome reference sequence, by mapping the enriched data set onto pseudo‐chromosomes derived from the capture probe target sequence, with a long‐range order of genes based on synteny of wheat with Brachypodium distachyon. Using this approach we are able to map the region and identify a set of deleted genes within the interval.     

Brachypodium distachyon: making hay with a wild grass

Brachypodium distachyon is a wild grass with a short life cycle. Although it is related to small grain cereals such as wheat, its genome is only a fraction of the size. A draft genome sequence is currently available, and molecular and genetic tools have been developed for transformation, mutagenesis and gene mapping. Accessions collected from across its ancestral range show a surprising degree of phenotypic variation in many traits, including those implicated in domestication of the cereals. Thus, given its rapid cycling time and ease of cultivation, Brachypodium will be a useful model for investigating problems in grass biology.     

First survey of the wheat chromosome 5A composition through a next generation sequencing approach

Wheat is one of the world's most important crops and is characterized by a large polyploid genome. One way to reduce genome complexity is to isolate single chromosomes using flow cytometry. Low coverage DNA sequencing can provide a snapshot of individual chromosomes, allowing a fast characterization of their main features and comparison with other genomes. We used massively parallel 454 pyrosequencing to obtain a 2x coverage of wheat chromosome 5A. The resulting sequence assembly was used to identify TEs, genes and miRNAs, as well as to infer a virtual gene order based on the synteny with other grass genomes. Repetitive elements account for more than 75% of the genome. Gene content was estimated considering non-redundant reads showing at least one match to ESTs or proteins. The results indicate that the coding fraction represents 1.08% and 1.3% of the short and long arm respectively, projecting the number of genes of the whole chromosome to approximately 5,000. 195 candidate miRNA precursors belonging to 16 miRNA families were identified. The 5A genes were used to search for syntenic relationships between grass genomes. The short arm is closely related to Brachypodium chromosome 4, sorghum chromosome 8 and rice chromosome 12; the long arm to regions of Brachypodium chromosomes 4 and 1, sorghum chromosomes 1 and 2 and rice chromosomes 9 and 3. From these similarities it was possible to infer the virtual gene order of 392 (5AS) and 1,480 (5AL) genes of chromosome 5A, which was compared to, and found to be largely congruent with the available physical map of this chromosome.     

Genome-wide survey of alternative splicing in the grass Brachypodium distachyon: a emerging model biosystem for plant functional genomics

A draft sequence of the genome of Brachypodium distachyon, the emerging grass model, was recently released. This represents a unique opportunity to determine its functional diversity compared to the genomes of other model species. Using homology mapping of assembled expressed sequence tags with chromosome scale pseudomolecules, we identified 128 alternative splicing events in B. distachyon. Our study identified that retention of introns is the major type of alternative splicing events (53%) in this plant and highlights the prevalence of splicing site recognition for definition of introns in plants. We have analyzed the compositional profiles of exon-intron junctions by base-pairing nucleotides with U1 snRNA which serves as a model for describing the possibility of sequence conservation. The alternative splicing isoforms identified in this study are novel and represent one of the potentially biologically significant means by which B. distachyon controls the function of its genes. Our observations serve as a basis to understand alternative splicing events of cereal crops with more complex genomes, like wheat or barley.     

Intraspecific sequence comparisons reveal similar rates of non-collinear gene insertion in the B and D genomes of bread wheat

ABSTRACT: BACKGROUND: Polyploidization is considered one of the main mechanisms of plant genome evolution. The presence of multiple copies of the same gene reduces selection pressure and permits sub-functionalization and neo-functionalization leading to plant diversification, adaptation and speciation. In bread wheat, polyploidization and the prevalence of transposable elements resulted in massive gene duplication and movement. As a result, the number of genes which are non-collinear to genomes of related species seems markedly increased in wheat. RESULTS: We used new-generation sequencing (NGS) to generate sequence of a Mb-sized region from wheat chromosome arm 3DS. Sequence assembly of 24 BAC clones resulted in two scaffolds of 1,264,820 and 333,768 bases. The sequence was annotated and compared to the homoeologous region on wheat chromosome 3B and orthologous loci of Brachypodium distachyon and rice. Among 39 coding sequences in the 3DS scaffolds, 32 have a homoeolog on chromosome 3B. In contrast, only fifteen and fourteen orthologs were identified in the corresponding regions in rice and Brachypodium, respectively. Interestingly, five pseudogenes were identified among the non-collinear coding sequences at the 3B locus, while none was found at the 3DS locus. CONCLUSION: Direct comparison of two Mb-sized regions of the B and D genomes of bread wheat revealed similar rates of non-collinear gene insertion in both genomes with a majority of gene duplications occurring before their divergence. Relatively low proportion of pseudogenes was identified among non-collinear coding sequences. Our data suggest that the pseudogenes did not originate from insertion of non-functional copies, but were formed later during the evolution of hexaploid wheat. Some evidence was found for gene erosion along the B genome locus.     

Germin-like proteins (GLPs) in cereal genomes: gene clustering and dynamic roles in plant defence

The recent release of the genome sequences of a number of crop and model plant species has made it possible to define the genome organisation and functional characteristics of specific genes and gene families of agronomic importance. For instance, Sorghum bicolor, maize (Zea mays) and Brachypodium distachyon genome sequences along with the model grass species rice (Oryza sativa) enable the comparative analysis of genes involved in plant defence. Germin-like proteins (GLPs) are a small, functionally and taxonomically diverse class of cupin-domain containing proteins that have recently been shown to cluster in an area of rice chromosome 8. The genomic location of this gene cluster overlaps with a disease resistance QTL that provides defence against two rice fungal pathogens (Magnaporthe oryzae and Rhizoctonia solani). Studies showing the involvement of GLPs in basal host resistance against powdery mildew (Blumeria graminis ssp.) have also been reported in barley and wheat. In this mini-review, we compare the close proximity of GLPs in publicly available cereal crop genomes and discuss the contribution that these proteins, and their genome sequence organisation, play in plant defence.     

Analysis of the wheat endosperm transcriptome

Among the cereals, wheat is the most widely grown geographically and is part of the staple diet in much of the world. Understanding how the cereal endosperm develops and functions will help generate better tools to manipulate grain qualities important to end-users. We used a genomics approach to identify and characterize genes that are expressed in the wheat endosperm. We analyzed the 17,949 publicly available wheat endosperm EST sequences to identify genes involved in the biological processes that occur within this tissue. Clustering and assembly of the ESTs resulted in the identification of 6,187 tentative unique genes, 2,358 of which formed contigs and 3,829 remained as singletons. A BLAST similarity search against the NCBI non-redundant sequence database revealed abundant messages for storage proteins, putative defense proteins, and proteins involved in starch and sucrose metabolism. The level of abundance of the putatively identified genes reflects the physiology of the developing endosperm. Half of the identified genes have unknown functions. Approximately 61% of the endosperm ESTs has been tentatively mapped in the hexaploid wheat genome. Using microarrays for global RNA profiling, we identified endosperm genes that are specifically up regulated in the developing grain.     

Sequencing and assembly of low copy and genic regions of isolated Triticum aestivum chromosome arm 7DS

The genome of bread wheat (Triticum aestivum) is predicted to be greater than 16 Gbp in size and consist predominantly of repetitive elements, making the sequencing and assembly of this genome a major challenge. We have reduced genome sequence complexity by isolating chromosome arm 7DS and applied second‐generation technology and appropriate algorithmic analysis to sequence and assemble low copy and genic regions of this chromosome arm. The assembly represents approximately 40% of the chromosome arm and all known 7DS genes. Comparison of the 7DS assembly with the sequenced genomes of rice (Oryza sativa) and Brachypodium distachyon identified large regions of conservation. The syntenic relationship between wheat, B. distachyon and O. sativa, along with available genetic mapping data, has been used to produce an annotated draft 7DS syntenic build, which is publicly available at http://www.wheatgenome.info . Our results suggest that the sequencing of isolated chromosome arms can provide valuable information of the gene content of wheat and is a step towards whole‐genome sequencing and variation discovery in this important crop.     

Strong population structure characterizes weediness gene evolution in the invasive grass species Brachypodium distachyon

Mediterranean annual grasses have invaded California and have replaced vast areas of native grassland. One of these invasive grasses is Brachypodium distachyon , a new model species for the grasses with extensive genomic resources and a nearly completed genome sequence. This study shows that the level of genetic variation in invaded California grasslands is lower compared to the native range in Eurasia. The invaded regions are characterized by highly differentiated populations of B. distachyon isolated by distance, most likely as a result of founder effects and a dearth of outcrossing events. EXP6 and EXP10 encoding α-expansins responsible for rapid growth, and AGL11 and AGL13 encoding proteins involved in vegetative phase regulation, appear to be under purifying selection with no evidence for local adaptation. Our data show that B. distachyon has diverged only recently from related Brachypodium species and that tetraploidization might have been as recent as a few thousand years ago. Observed low genetic variation in EXP10 and AGL13 appears to have been present in Eurasia before tetraploidization, potentially as a result of strong selective pressures on advantageous mutations, which are most likely responsible for its fast growth and rapid completion of its life cycle.     

Exploiting the Brachypodium Tool Box in cereal and grass research

It is now a decade since Brachypodium distachyon (Brachypodium) was suggested as a model species for temperate grasses and cereals. Since then transformation protocols, large expressed sequence tag (EST) databases, tools for forward and reverse genetic screens, highly refined cytogenetic probes, germplasm collections and, recently, a complete genome sequence have been generated. In this review, we will describe the current status of the Brachypodium Tool Box and how it is beginning to be applied to study a range of biological traits. Further, as genomic analysis of larger cereals and forage grasses genomes are becoming easier, we will re-evaluate Brachypodium as a model species. We suggest that there remains an urgent need to employ reverse genetic and functional genomic approaches to identify the functionality of key genetic elements, which could be employed subsequently in plant breeding programmes; and a requirement for a Pooideae reference genome to aid assembling large pooid genomes. Brachypodium is an ideal system for functional genomic studies, because of its easy growth requirements, small physical stature, and rapid life cycle, coupled with the resources offered by the Brachypodium Tool Box.     

Molecular cloning of a new receptor-like kinase gene encoded at the Lr10 disease resistance locus of wheat

More than 100 resistance genes against wheat rust pathogens have been described in wheat and its relatives. Although many of them have been extensively used in wheat resistance breeding, none of these resistance loci has yet been analyzed at the molecular level. By screening a set of near-isogenic lines carrying different leaf rust resistance genes with a wheat probe encoding a serine/threonine protein kinase, we detected a polymorphic DNA fragment in the line with the Lr10 resistance gene. This fragment mapped to the Lr10 disease resistance locus and encodes a receptor-like protein kinase which we called LRK10. LRK10 contains a new type of extracellular domain not found in known plant or animal receptor kinases. Several conserved amino acids in S-domain glycoproteins and receptor-like kinases were also found in LRK10, suggesting that LRK10 and S-domain proteins belong to the same superfamily of specific recognition molecules in plants. Lrk10 was expressed at low levels in young seedlings and belongs to a gene family. Analysis of wheat lines with and without the Lr10 gene demonstrated that Lrk10 and Lr10 belong to the same genetic locus. We conclude that gene isolation based on protein kinase homology can identify new receptor domains and provide candidates for disease resistance genes in the complex wheat genome.     

Sub-proteome analysis of novel IgE-binding proteins from Bermuda grass pollen

Bermuda grass (Cynodon dactylon) pollen (BGP) is one of the most common causes of airway allergic disease, and has been shown to contain over 12 allergenic proteins on 1-D immunoglobulin E (IgE) immunoblots. However, only a few allergens have been identified and characterized. Cyn d 1 is a major allergen and the most abundant protein in BGP, representing 15% of the whole-pollen extract. To investigate variability in the IgE-reactive patterns of BGP-sensitized patients and to identify other prevalent allergens, a BGP extract was passed through an affinity column to remove Cyn d 1, and the non-bound material was collected and analyzed by 2-DE. IgE-reactive proteins were subsequently characterized by immunoblotting using serum samples from ten BGP-allergic patients. The prevalent IgE-reactive proteins were identified by MALDI-TOF MS, N-terminal sequence similarity, and LC-MS/MS. Here, we present a sub-proteome approach for allergen investigation and its use for determining BGP 2-DE profiles and identifying six novel allergens.