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

Characterization and Use in Wheat Breeding of Leaf Rust Resistance Genes from Durable Varieties

Simple Summary

Wheat leaf rust is one of the most significant diseases worldwide, incited by a parasitic fungus which infects leaves, affecting grain yield. This pathogen is spread by the wind over large areas through microscopic spores. This huge number of spores favors the selection of virulent forms; therefore, there is a continuous need for new resistance genes to control this disease without fungicides. These resistant genes are naturally found in resistant wheat varieties and can be introduced by standard crosses. In this work, seven resistant genes were introduced into several commercial susceptible varieties. The selection of resistance genes was assisted by DNA markers that are close to these genes on the chromosome. Additionally, the selection of desirable traits from the commercial variety was also assisted by DNA markers to accelerate the process. In field testing, the varieties developed here were resistant to leaf rust, and suitable for commercial use.

Abstract

Leaf rust is one of the most significant diseases of wheat worldwide. In Argentina, it is one of the main reasons for variety replacement that becomes susceptible after large-scale use. Some varieties showed durable resistance to this disease, including Buck Manantial and Sinvalocho MA. RILs (Recombinant Inbred Lines) were developed for each of these varieties and used in genetics studies to identify components of resistance, both in greenhouse inoculations using leaf rust races, and in field evaluations under natural population infections. In Buck Manantial, the APR gene LrBMP1 was associated with resistance in field tests. In crosses involving Sinvalocho MA, four genes were previously identified and associated with resistance in field testing: APR (Adult Plant Resistance) gene LrSV1, the APR genetic system LrSV2 + LrcSV2 and the ASR (All Stage Resistance) gene LrG6. Using backcrosses, LrBMP1 was introgressed in four commercial susceptible varieties and LrSV1, LrSV2 + LrcSV2 and LrG6 were simultaneously introgressed in three susceptible commercial varieties. The use of molecular markers for recurrent parent background selection allowed us to select resistant lines with more than 80% similarity to commercial varieties. Additionally, progress towards positional cloning of the genetic system LrSV2 + LrcSV2 for leaf rust APR is reported.

The NLR-Annotator Tool Enables Annotation of the Intracellular Immune Receptor Repertoire

Disease resistance genes encoding nucleotide-binding and leucine-rich repeat (NLR) intracellular immune receptor proteins detect pathogens by the presence of pathogen effectors. Plant genomes typically contain hundreds of NLR-encoding genes. The availability of the hexaploid wheat (Triticum aestivum) cultivar Chinese Spring reference genome allows a detailed study of its NLR complement. However, low NLR expression and high intrafamily sequence homology hinder their accurate annotation. Here, we developed NLR-Annotator, a software tool for in silico NLR identification independent of transcript support. Although developed for wheat, we demonstrate the universal applicability of NLR-Annotator across diverse plant taxa. We applied our tool to wheat and combined it with a transcript-validated subset of genes from the reference gene annotation to characterize the structure, phylogeny, and expression profile of the NLR gene family. We detected 3,400 full-length NLR loci, of which 1,560 were confirmed as expressed genes with intact open reading frames. NLRs with integrated domains mostly group in specific subclades. Members of another subclade predominantly locate in close physical proximity to NLRs carrying integrated domains, suggesting a paired helper function. Most NLRs (88%) display low basal expression (in the lower 10 percentile of transcripts). In young leaves subjected to biotic stress, we found up-regulation of 266 of the NLRs To illustrate the utility of our tool for the positional cloning of resistance genes, we estimated the number of NLR genes within the intervals of mapped rust resistance genes. Our study will support the identification of functional resistance genes in wheat to accelerate the breeding and engineering of disease-resistant varieties.

Cloning of the wheat Yr15 resistance gene sheds light on the plant tandem kinase-pseudokinase family

Yellow rust, caused by Puccinia striiformis f. sp. tritici (Pst), is a devastating fungal disease threatening much of global wheat production. Race-specific resistance (R)-genes are used to control rust diseases, but the rapid emergence of virulent Pst races has prompted the search for a more durable resistance. Here, we report the cloning of Yr15, a broad-spectrum R-gene derived from wild emmer wheat, which encodes a putative kinase-pseudokinase protein, designated as wheat tandem kinase 1, comprising a unique R-gene structure in wheat. The existence of a similar gene architecture in 92 putative proteins across the plant kingdom, including the barley RPG1 and a candidate for Ug8, suggests that they are members of a distinct family of plant proteins, termed here tandem kinase-pseudokinases (TKPs). The presence of kinase-pseudokinase structure in both plant TKPs and the animal Janus kinases sheds light on the molecular evolution of immune responses across these two kingdoms.

Advances in Wheat and Pathogen Genomics: Implications for Disease Control

The gene pool of wheat and its wild and domesticated relatives contains a plethora of resistance genes that can be exploited to make wheat more resilient to pathogens. Only a few of these genes have been isolated and studied at the molecular level. In recent years, we have seen a shift from classical breeding to genomics-assisted breeding, which makes use of the enormous advancements in DNA sequencing and high-throughput molecular marker technologies for wheat improvement. These genomic advancements have the potential to transform wheat breeding in the near future and to significantly increase the speed and precision at which new cultivars can be bred. This review highlights the genomic improvements that have been made in wheat and its pathogens over the past years and discusses their implications for disease-resistance breeding.

Impact of transposable elements on genome structure and evolution in bread wheat

BACKGROUND

Transposable elements (TEs) are major components of large plant genomes and main drivers of genome evolution. The most recent assembly of hexaploid bread wheat recovered the highly repetitive TE space in an almost complete chromosomal context and enabled a detailed view into the dynamics of TEs in the A, B, and D subgenomes.

RESULTS

The overall TE content is very similar between the A, B, and D subgenomes, although we find no evidence for bursts of TE amplification after the polyploidization events. Despite the near-complete turnover of TEs since the subgenome lineages diverged from a common ancestor, 76% of TE families are still present in similar proportions in each subgenome. Moreover, spacing between syntenic genes is also conserved, even though syntenic TEs have been replaced by new insertions over time, suggesting that distances between genes, but not sequences, are under evolutionary constraints. The TE composition of the immediate gene vicinity differs from the core intergenic regions. We find the same TE families to be enriched or depleted near genes in all three subgenomes. Evaluations at the subfamily level of timed long terminal repeat-retrotransposon insertions highlight the independent evolution of the diploid A, B, and D lineages before polyploidization and cases of concerted proliferation in the AB tetraploid.

CONCLUSIONS

Even though the intergenic space is changed by the TE turnover, an unexpected preservation is observed between the A, B, and D subgenomes for features like TE family proportions, gene spacing, and TE enrichment near genes.

Chromosome-scale comparative sequence analysis unravels molecular mechanisms of genome dynamics between two wheat cultivars

BACKGROUND

Recent improvements in DNA sequencing and genome scaffolding have paved the way to generate high-quality de novo assemblies of pseudomolecules representing complete chromosomes of wheat and its wild relatives. These assemblies form the basis to compare the dynamics of wheat genomes on a megabase scale.

RESULTS

Here, we provide a comparative sequence analysis of the 700-megabase chromosome 2D between two bread wheat genotypes-the old landrace Chinese Spring and the elite Swiss spring wheat line 'CH Campala Lr22a'. Both chromosomes were assembled into megabase-sized scaffolds. There is a high degree of sequence conservation between the two chromosomes. Analysis of large structural variations reveals four large indels of more than 100 kb. Based on the molecular signatures at the breakpoints, unequal crossing over and double-strand break repair were identified as the molecular mechanisms that caused these indels. Three of the large indels affect copy number of NLRs, a gene family involved in plant immunity. Analysis of SNP density reveals four haploblocks of 4, 8, 9 and 48 Mb with a 35-fold increased SNP density compared to the rest of the chromosome. Gene content across the two chromosomes was highly conserved. Ninety-nine percent of the genic sequences were present in both genotypes and the fraction of unique genes ranged from 0.4 to 0.7%.

CONCLUSIONS

This comparative analysis of two high-quality chromosome assemblies enabled a comprehensive assessment of large structural variations and gene content. The insight obtained from this analysis will form the basis of future wheat pan-genome studies.

Impact of transposable elements on polyploid plant genomes

BACKGROUND

The growing wealth of knowledge on whole-plant genome sequences is highlighting the key role of transposable elements (TEs) in plant evolution, as a driver of drastic changes in genome size and as a source of an important number of new coding and regulatory sequences. Together with polyploidization events, TEs should thus be considered the major players in evolution of plants.

SCOPE

This review outlines the major mechanisms by which TEs impact plant genome evolution and how polyploidy events can affect these impacts, and vice versa. These include direct effects on genes, by providing them with new coding or regulatory sequences, an effect on the epigenetic status of the chromatin close to genes, and more subtle effects by imposing diverse evolutionary constraints to different chromosomal regions. These effects are particularly relevant after polyploidization events. Polyploidization often induces bursts of transposition probably due to a relaxation in their epigenetic control, and, in the short term, this can increase the rate of gene mutations and changes in gene regulation due to the insertion of TEs next to or into genes. Over longer times, TE bursts may induce global changes in genome structure due to inter-element recombination including losses of large genome regions and chromosomal rearrangements that reduce the genome size and the chromosome number as part of a process called diploidization.

CONCLUSIONS

TEs play an essential role in genome and gene evolution, in particular after polyploidization events. Polyploidization can induce TE activity that may explain part of the new phenotypes observed. TEs may also play a role in the diploidization that follows polyploidization events. However, the extent to which TEs contribute to diploidization and fractionation bias remains unclear. Investigating the multiple factors controlling TE dynamics and the nature of ancient and recent polyploid genomes may shed light on these processes.

Effects of ploidy level and haplotype on variation of photosynthetic traits: Novel evidence from two Fragaria species

To reveal the effects of ploidy level and haplotype on photosynthetic traits, we chose 175 genotypes of wild strawberries belonging to two haplotypes at two types of ploidy levels (diploidy and tetraploidy) and measured photosynthetic traits. Our results revealed that ploidy significantly affected the characteristics of light-response curves, CO2-response curves, and leaf gas exchange parameters, except intercellular CO2 concentration (Ci). Tetraploid species had a lower light saturation point (LSP) and CO2 saturation point (CSP), higher light compensation point (LCP), dark respiration (Rd), and CO2 compensation point (CCP) than diploid species. Furthermore, tetraploid species have lower photosynthetic capacity than diploid species, including net photosynthetic rate (Pn), stomatal conductivity (Gs), and transpiration rate (Tr). In addition, haplotype had a significant effect on LSP, CSP, Tr, and Ci as well as a significant interactive effect between ploidy and haplotype on the maximal photosynethic rate of the light-response curve and Rd. Most of the variance existed within haplotypes among individuals. These results suggest that polyploidization was the main driver for the evolution of photosynthesis with increasing ploidy level (i.e. from diploidy to tetraploidy in Fragaria species), while the origin of a chromosome could also affect the photosynthetic traits and the polyploidization effect on photosynthetic traits.

Simulating Next-Generation Sequencing Datasets from Empirical Mutation and Sequencing Models

An obstacle to validating and benchmarking methods for genome analysis is that there are few reference datasets available for which the “ground truth” about the mutational landscape of the sample genome is known and fully validated. Additionally, the free and public availability of real human genome datasets is incompatible with the preservation of donor privacy. In order to better analyze and understand genomic data, we need test datasets that model all variants, reflecting known biology as well as sequencing artifacts. Read simulators can fulfill this requirement, but are often criticized for limited resemblance to true data and overall inflexibility. We present NEAT (NExt-generation sequencing Analysis Toolkit), a set of tools that not only includes an easy-to-use read simulator, but also scripts to facilitate variant comparison and tool evaluation. NEAT has a wide variety of tunable parameters which can be set manually on the default model or parameterized using real datasets. The software is freely available at github.com/zstephens/neat-genreads.

Molecular genetics and evolution of disease resistance in cereals

Contents 320 I. 320 II. 321 III. 321 IV. 322 V. 324 VI. 328 VII. 329 330 References 330 SUMMARY: Cereal crops produce a large part of the globally consumed food and feed. Because of the constant presence of devastating pathogens, the molecular characterization of disease resistance is a major research area and highly relevant for breeding. There has been recent and accelerating progress in the understanding of three distinct resistance mechanisms in cereals: resistance conferred by plasma membrane-localized receptor proteins; race-specific resistance conferred by intracellular immune receptors; and quantitative disease resistance. Intracellular immune receptors provide a particularly rich source for evolutionary studies, and have, for example, resulted in the recent discovery of a novel detection mechanism based on integrated decoy domains. Evolutionary studies have also revealed the origins of active resistance genes in both wild progenitors of today's cereals as well as in cultivated forms. In addition, independent evolution of orthologous genes in related cereals has resulted in resistance to different pathogen species. Quantitative resistance genes have been best characterized in wheat. The quantitative resistance genes identified so far in wheat encode transporter proteins or unusual kinase proteins. The recent discoveries in these three different resistance mechanisms have contributed to the basic molecular understanding of cereal immunity against pathogens and have suggested novel applications for resistance breeding.

Suppressed recombination and unique candidate genes in the divergent haplotype encoding Fhb1, a major Fusarium head blight resistance locus in wheat

Fine mapping and sequencing revealed 28 genes in the non-recombining haplotype containing Fhb1 . Of these, only a GDSL lipase gene shows a pathogen-dependent expression pattern. Fhb1 is a prominent Fusarium head blight resistance locus of wheat, which has been successfully introgressed in adapted breeding material, where it confers a significant increase in overall resistance to the causal pathogen Fusarium graminearum and the fungal virulence factor and mycotoxin deoxynivalenol. The Fhb1 region has been resolved for the susceptible wheat reference genotype Chinese Spring, yet the causal gene itself has not been identified in resistant cultivars. Here, we report the establishment of a 1 Mb contig embracing Fhb1 in the donor line CM-82036. Sequencing revealed that the region of Fhb1 deviates from the Chinese Spring reference in DNA size and gene content, which explains the repressed recombination at the locus in the performed fine mapping. Differences in genes expression between near-isogenic lines segregating for Fhb1 challenged with F. graminearum or treated with mock were investigated in a time-course experiment by RNA sequencing. Several candidate genes were identified, including a pathogen-responsive GDSL lipase absent in susceptible lines. The sequence of the Fhb1 region, the resulting list of candidate genes, and near-diagnostic KASP markers for Fhb1 constitute a valuable resource for breeding and further studies aiming to identify the gene(s) responsible for F. graminearum and deoxynivalenol resistance.

Zn2+ blocks annealing of complementary single-stranded DNA in a sequence-selective manner

Zinc is the second most abundant trace element essential for all living organisms. In human body, 30–40% of the total zinc ion (Zn²⁺) is localized in the nucleus. Intranuclear free Zn²⁺ sparks caused by reactive oxygen species have been observed in eukaryotic cells, but question if these free Zn²⁺ outrages could have affected annealing of complementary single-stranded (ss) DNA, a crucial step in DNA synthesis, repair and recombination, has never been raised. Here the author reports that Zn²⁺ blocks annealing of complementary ssDNA in a sequence-selective manner under near-physiological conditions as demonstrated in vitro using a low-temperature EDTA-free agarose gel electrophoresis (LTEAGE) procedure. Specifically, it is shown that Zn²⁺ does not block annealing of repetitive DNA sequences lacking CG/GC sites that are the major components of junk DNA. It is also demonstrated that Zn²⁺ blocks end-joining of double-stranded (ds) DNA fragments with 3′ overhangs mimicking double-strand breaks, and prevents renaturation of long stretches (>1 kb) of denatured dsDNA, in which Zn²⁺-tolerant intronic DNA provides annealing protection on otherwise Zn²⁺-sensitive coding DNA. These findings raise a challenging hypothesis that Zn²⁺-ssDNA interaction might be among natural forces driving eukaryotic genomes to maintain the Zn²⁺-tolerant repetitive DNA for adapting to the Zn²⁺-rich nucleus.

Sequencing of chloroplast genomes from wheat, barley, rye and their relatives provides a detailed insight into the evolution of the Triticeae tribe

Using Roche/454 technology, we sequenced the chloroplast genomes of 12 Triticeae species, including bread wheat, barley and rye, as well as the diploid progenitors and relatives of bread wheat Triticum urartu, Aegilops speltoides and Ae. tauschii. Two wild tetraploid taxa, Ae. cylindrica and Ae. geniculata, were also included. Additionally, we incorporated wild Einkorn wheat Triticum boeoticum and its domesticated form T. monococcum and two Hordeum spontaneum (wild barley) genotypes. Chloroplast genomes were used for overall sequence comparison, phylogenetic analysis and dating of divergence times. We estimate that barley diverged from rye and wheat approximately 8–9 million years ago (MYA). The genome donors of hexaploid wheat diverged between 2.1–2.9 MYA, while rye diverged from Triticum aestivum approximately 3–4 MYA, more recently than previously estimated. Interestingly, the A genome taxa T. boeoticum and T. urartu were estimated to have diverged approximately 570,000 years ago. As these two have a reproductive barrier, the divergence time estimate also provides an upper limit for the time required for the formation of a species boundary between the two. Furthermore, we conclusively show that the chloroplast genome of hexaploid wheat was contributed by the B genome donor and that this unknown species diverged from Ae. speltoides about 980,000 years ago. Additionally, sequence alignments identified a translocation of a chloroplast segment to the nuclear genome which is specific to the rye/wheat lineage. We propose the presented phylogeny and divergence time estimates as a reference framework for future studies on Triticeae.

Sequence-based mapping of the polyploid wheat genome

The emergence of new sequencing technologies has provided fast and cost-efficient strategies for high-resolution mapping of complex genomes. Although these approaches hold great promise to accelerate genome analysis, their application in studying genetic variation in wheat has been hindered by the complexity of its polyploid genome. Here, we applied the next-generation sequencing of a wheat doubled-haploid mapping population for high-resolution gene mapping and tested its utility for ordering shotgun sequence contigs of a flow-sorted wheat chromosome. A bioinformatical pipeline was developed for reliable variant analysis of sequence data generated for polyploid wheat mapping populations. The results of variant mapping were consistent with the results obtained using the wheat 9000 SNP iSelect assay. A reference map of the wheat genome integrating 2740 gene-associated single-nucleotide polymorphisms from the wheat iSelect assay, 1351 diversity array technology, 118 simple sequence repeat/sequence-tagged sites, and 416,856 genotyping-by-sequencing markers was developed. By analyzing the sequenced megabase-size regions of the wheat genome we showed that mapped markers are located within 40-100 kb from genes providing a possibility for high-resolution mapping at the level of a single gene. In our population, gene loci controlling a seed color phenotype cosegregated with 2459 markers including one that was located within the red seed color gene. We demonstrate that the high-density reference map presented here is a useful resource for gene mapping and linking physical and genetic maps of the wheat genome.

Separating homeologs by phasing in the tetraploid wheat transcriptome

Background

The high level of identity among duplicated homoeologous genomes in tetraploid pasta wheat presents substantial challenges for de novo transcriptome assembly. To solve this problem, we develop a specialized bioinformatics workflow that optimizes transcriptome assembly and separation of merged homoeologs. To evaluate our strategy, we sequence and assemble the transcriptome of one of the diploid ancestors of pasta wheat, and compare both assemblies with a benchmark set of 13,472 full-length, non-redundant bread wheat cDNAs.

Results

A total of 489 million 100 bp paired-end reads from tetraploid wheat assemble in 140,118 contigs, including 96% of the benchmark cDNAs. We used a comparative genomics approach to annotate 66,633 open reading frames. The multiple k-mer assembly strategy increases the proportion of cDNAs assembled full-length in a single contig by 22% relative to the best single k-mer size. Homoeologs are separated using a post-assembly pipeline that includes polymorphism identification, phasing of SNPs, read sorting, and re-assembly of phased reads. Using a reference set of genes, we determine that 98.7% of SNPs analyzed are correctly separated by phasing.

Conclusions

Our study shows that de novo transcriptome assembly of tetraploid wheat benefit from multiple k-mer assembly strategies more than diploid wheat. Our results also demonstrate that phasing approaches originally designed for heterozygous diploid organisms can be used to separate the close homoeologous genomes of tetraploid wheat. The predicted tetraploid wheat proteome and gene models provide a valuable tool for the wheat research community and for those interested in comparative genomic studies.

Genome-wide analysis of Stowaway-like MITEs in wheat reveals high sequence conservation, gene association, and genomic diversification

The diversity and evolution of wheat (Triticum-Aegilops group) genomes is determined, in part, by the activity of transposable elements that constitute a large fraction of the genome (up to 90%). In this study, we retrieved sequences from publicly available wheat databases, including a 454-pyrosequencing database, and analyzed 18,217 insertions of 18 Stowaway-like miniature inverted-repeat transposable element (MITE) families previously characterized in wheat that together account for approximately 1.3 Mb of sequence. All 18 families showed high conservation in length, sequence, and target site preference. Furthermore, approximately 55% of the elements were inserted in transcribed regions, into or near known wheat genes. Notably, we observed significant correlation between the mean length of the MITEs and their copy number. In addition, the genomic composition of nine MITE families was studied by real-time quantitative polymerase chain reaction analysis in 40 accessions of Triticum spp. and Aegilops spp., including diploids, tetraploids, and hexaploids. The quantitative polymerase chain reaction data showed massive and significant intraspecific and interspecific variation as well as genome-specific proliferation and nonadditive quantities in the polyploids. We also observed significant differences in the methylation status of the insertion sites among MITE families. Our data thus suggest a possible role for MITEs in generating genome diversification and in the establishment of nascent polyploid species in wheat.

Arabidopsis Fused kinase and the Kinesin-12 subfamily constitute a signalling module required for phragmoplast expansion

The conserved Fused kinase plays vital but divergent roles in many organisms from Hedgehog signalling in Drosophila to polarization and chemotaxis in Dictyostelium. Previously we have shown that Arabidopsis Fused kinase termed TWO-IN-ONE (TIO) is essential for cytokinesis in both sporophytic and gametophytic cell types. Here using in vivo imaging of GFP-tagged microtubules in dividing microspores we show that TIO is required for expansion of the phragmoplast. We identify the phragmoplast-associated kinesins, PAKRP1/Kinesin-12A and PAKRP1L/Kinesin-12B, as TIO-interacting proteins and determine TIO-Kinesin-12 interaction domains and their requirement in male gametophytic cytokinesis. Our results support the role of TIO as a functional protein kinase that interacts with Kinesin-12 subfamily members mainly through the C-terminal ARM repeat domain, but with a contribution from the N-terminal kinase domain. The interaction of TIO with Kinesin proteins and the functional requirement of their interaction domains support the operation of a Fused kinase signalling module in phragmoplast expansion that depends upon conserved structural features in diverse Fused kinases.

Marker utility of miniature inverted-repeat transposable elements for wheat biodiversity and evolution

Transposable elements (TEs) account for up to 80% of the wheat genome and are considered one of the main drivers of wheat genome evolution. However, the contribution of TEs to the divergence and evolution of wheat genomes is not fully understood. In this study, we have developed 55 miniature inverted-repeat transposable element (MITE) markers that are based on the presence/absence of an element, with over 60% of these 55 MITE insertions associated with wheat genes. We then applied these markers to assess genetic diversity among Triticum and Aegilops species, including diploid (AA, BB and DD genomes), tetraploid (BBAA genome) and hexaploid (BBAADD genome) species. While 18.2% of the MITE markers showed similar insertions in all species indicating that those are fossil insertions, 81.8% of the markers showed polymorphic insertions among species, subspecies, and accessions. Furthermore, a phylogenetic analysis based on MITE markers revealed that species were clustered based on genus, genome composition, and ploidy level, while 47.13% genetic divergence was observed between the two main clusters, diploids versus polyploids. In addition, we provide evidence for MITE dynamics in wild emmer populations. The use of MITEs as evolutionary markers might shed more light on the origin of the B-genome of polyploid wheat.

Ancient diversity of splicing motifs and protein surfaces in the wild emmer wheat (Triticum dicoccoides) LR10 coiled coil (CC) and leucine-rich repeat (LRR) domains

In this study, we explore the diversity and its distribution along the wheat leaf rust resistance protein LR10 three-dimensional structure. Lr10 is a leaf rust resistance gene encoding a coiled coil-nucleotide-binding site-leucine-rich repeat (CC-NBS-LRR) class of protein. Lr10 was cloned and sequenced from 58 accessions representing diverse habitats of wild emmer wheat in Israel. Nucleotide diversity was very high relative to other wild emmer wheat genes (π= 0.029). The CC domain was found to be the most diverse domain and subject to positive selection. Superimposition of the diversity on the CC three-dimensional structure showed that some of the variable and positively selected residues were solvent exposed and may interact with other proteins. The LRR domain was relatively conserved, but showed a hotspot of amino acid variation between two haplotypes in the ninth repeat. This repeat was longer than the other LRRs, and three-dimensional modelling suggested that an extensive α helix structure was formed in this region. The two haplotypes also differed in splicing regulation motifs. In genotypes with one haplotype, an intron was alternatively spliced in this region, whereas, in genotypes with the other haplotype, this intron did not splice at all. The two haplotypes are proposed to be ancient and maintained by balancing selection.

The Genomes of All Angiosperms: A Call for a Coordinated Global Census

Recent advances in biological instrumentation and associated experimental technologies now permit an unprecedented efficiency and scale for the acquisition of genomic data, at ever-decreasing costs. Further advances, with accompanying decreases in cost, are expected in the very near term. It now becomes appropriate to discuss the best uses of these technologies in the context of the angiosperms. This white paper proposes a complete genomic census of the approximately 500,000 species of flowering plants, outlines the goals of this census and their value, and provides a road map towards achieving these goals in a timely manner.

Duplication and partitioning in evolution and function of homoeologous Q loci governing domestication characters in polyploid wheat

The Q gene encodes an AP2-like transcription factor that played an important role in domestication of polyploid wheat. The chromosome 5A Q alleles (5AQ and 5Aq) have been well studied, but much less is known about the q alleles on wheat homoeologous chromosomes 5B (5Bq) and 5D (5Dq). We investigated the organization, evolution, and function of the Q/q homoeoalleles in hexaploid wheat (Triticum aestivum L.). Q/q gene sequences are highly conserved within and among the A, B, and D genomes of hexaploid wheat, the A and B genomes of tetraploid wheat, and the A, S, and D genomes of the diploid progenitors, but the intergenic regions of the Q/q locus are highly divergent among homoeologous genomes. Duplication of the q gene 5.8 Mya was likely followed by selective loss of one of the copies from the A genome progenitor and the other copy from the B, D, and S genomes. A recent V(329)-to-I mutation in the A lineage is correlated with the Q phenotype. The 5Bq homoeoalleles became a pseudogene after allotetraploidization. Expression analysis indicated that the homoeoalleles are coregulated in a complex manner. Combined phenotypic and expression analysis indicated that, whereas 5AQ plays a major role in conferring domestication-related traits, 5Dq contributes directly and 5Bq indirectly to suppression of the speltoid phenotype. The evolution of the Q/q loci in polyploid wheat resulted in the hyperfunctionalization of 5AQ, pseudogenization of 5Bq, and subfunctionalization of 5Dq, all contributing to the domestication traits.

The impact of Ty3-gypsy group LTR retrotransposons Fatima on B-genome specificity of polyploid wheats

BACKGROUND

Transposable elements (TEs) are a rapidly evolving fraction of the eukaryotic genomes and the main contributors to genome plasticity and divergence. Recently, occupation of the A- and D-genomes of allopolyploid wheat by specific TE families was demonstrated. Here, we investigated the impact of the well-represented family of gypsy LTR-retrotransposons, Fatima, on B-genome divergence of allopolyploid wheat using the fluorescent in situ hybridisation (FISH) method and phylogenetic analysis.

RESULTS

FISH analysis of a BAC clone (BAC_2383A24) initially screened with Spelt1 repeats demonstrated its predominant localisation to chromosomes of the B-genome and its putative diploid progenitor Aegilops speltoides in hexaploid (genomic formula, BBAADD) and tetraploid (genomic formula, BBAA) wheats as well as their diploid progenitors. Analysis of the complete BAC_2383A24 nucleotide sequence (113,605 bp) demonstrated that it contains 55.6% TEs, 0.9% subtelomeric tandem repeats (Spelt1), and five genes. LTR retrotransposons are predominant, representing 50.7% of the total nucleotide sequence. Three elements of the gypsy LTR retrotransposon family Fatima make up 47.2% of all the LTR retrotransposons in this BAC. In situ hybridisation of the Fatima_2383A24-3 subclone suggests that individual representatives of the Fatima family contribute to the majority of the B-genome specific FISH pattern for BAC_2383A24. Phylogenetic analysis of various Fatima elements available from databases in combination with the data on their insertion dates demonstrated that the Fatima elements fall into several groups. One of these groups, containing Fatima_2383A24-3, is more specific to the B-genome and proliferated around 0.5-2.5 MYA, prior to allopolyploid wheat formation.

CONCLUSION

The B-genome specificity of the gypsy-like Fatima, as determined by FISH, is explained to a great degree by the appearance of a genome-specific element within this family for Ae. speltoides. Moreover, its proliferation mainly occurred in this diploid species before it entered into allopolyploidy.Most likely, this scenario of emergence and proliferation of the genome-specific variants of retroelements, mainly in the diploid species, is characteristic of the evolution of all three genomes of hexaploid wheat.

Mechanisms of haplotype divergence at the RGA08 nucleotide-binding leucine-rich repeat gene locus in wild banana (Musa balbisiana)

BACKGROUND

Comparative sequence analysis of complex loci such as resistance gene analog clusters allows estimating the degree of sequence conservation and mechanisms of divergence at the intraspecies level. In banana (Musa sp.), two diploid wild species Musa acuminata (A genome) and Musa balbisiana (B genome) contribute to the polyploid genome of many cultivars. The M. balbisiana species is associated with vigour and tolerance to pests and disease and little is known on the genome structure and haplotype diversity within this species. Here, we compare two genomic sequences of 253 and 223 kb corresponding to two haplotypes of the RGA08 resistance gene analog locus in M. balbisiana "Pisang Klutuk Wulung" (PKW).

RESULTS

Sequence comparison revealed two regions of contrasting features. The first is a highly colinear gene-rich region where the two haplotypes diverge only by single nucleotide polymorphisms and two repetitive element insertions. The second corresponds to a large cluster of RGA08 genes, with 13 and 18 predicted RGA genes and pseudogenes spread over 131 and 152 kb respectively on each haplotype. The RGA08 cluster is enriched in repetitive element insertions, in duplicated non-coding intergenic sequences including low complexity regions and shows structural variations between haplotypes. Although some allelic relationships are retained, a large diversity of RGA08 genes occurs in this single M. balbisiana genotype, with several RGA08 paralogs specific to each haplotype. The RGA08 gene family has evolved by mechanisms of unequal recombination, intragenic sequence exchange and diversifying selection. An unequal recombination event taking place between duplicated non-coding intergenic sequences resulted in a different RGA08 gene content between haplotypes pointing out the role of such duplicated regions in the evolution of RGA clusters. Based on the synonymous substitution rate in coding sequences, we estimated a 1 million year divergence time for these M. balbisiana haplotypes.

CONCLUSIONS

A large RGA08 gene cluster identified in wild banana corresponds to a highly variable genomic region between haplotypes surrounded by conserved flanking regions. High level of sequence identity (70 to 99%) of the genic and intergenic regions suggests a recent and rapid evolution of this cluster in M. balbisiana.

Delineating the structural, functional and evolutionary relationships of sucrose phosphate synthase gene family II in wheat and related grasses

BACKGROUND

Sucrose phosphate synthase (SPS) is an important component of the plant sucrose biosynthesis pathway. In the monocotyledonous Poaceae, five SPS genes have been identified. Here we present a detailed analysis of the wheat SPSII family in wheat. A set of homoeologue-specific primers was developed in order to permit both the detection of sequence variation, and the dissection of the individual contribution of each homoeologue to the global expression of SPSII.

RESULTS

The expression in bread wheat over the course of development of various sucrose biosynthesis genes monitored on an Affymetrix array showed that the SPS genes were regulated over time and space. SPSII homoeologue-specific assays were used to show that the three homoeologues contributed differentially to the global expression of SPSII. Genetic mapping placed the set of homoeoloci on the short arms of the homoeologous group 3 chromosomes. A resequencing of the A and B genome copies allowed the detection of four haplotypes at each locus. The 3B copy includes an unspliced intron. A comparison of the sequences of the wheat SPSII orthologues present in the diploid progenitors einkorn, goatgrass and Triticum speltoides, as well as in the more distantly related species barley, rice, sorghum and purple false brome demonstrated that intronic sequence was less well conserved than exonic. Comparative sequence and phylogenetic analysis of SPSII gene showed that false purple brome was more similar to Triticeae than to rice. Wheat - rice synteny was found to be perturbed at the SPS region.

CONCLUSION

The homoeologue-specific assays will be suitable to derive associations between SPS functionality and key phenotypic traits. The amplicon sequences derived from the homoeologue-specific primers are informative regarding the evolution of SPSII in a polyploid context.

Cloning and phylogenetic analysis of phytoene synthase 1 (Psy1) genes in common wheat and related species

Cloning and phylogenetic analysis of Psy1 genes in common wheat and its relatives would help to understand the genetic diversity and evolution of Psy1 gene in common wheat and its related species. In the present study, common wheat (AABBDD) and eight relative species, including T. urartu (A(u)A(u)), T. boeoticum (A(m)A(m)), T. monococcum (A(m)A(m)), Ae. speltoides (SS), Ae. tauschii (DD), T. dicoccoides (AABB), T. dicoccum (AABB) and T. spelta (AABBDD), were sampled for the isolation of novel alleles at Psy1-A1, Psy1-B1/Psy1-S1 and Psy1-D1 loci corresponding to common wheat Psy1 genes, and 27 new alleles were identified at these loci, designated Psy1-A1f through Psy1-A1k, Psy1-A1m and Psy1-A1n, Psy1-B1h through Psy1-B1m, Psy1-S1a through Psy1-S1c, Psy1-D1a through Psy1-D1j, respectively. The genes contained six exons and five introns, and the sequences of exons were more conserved compared with those of introns. The Psy1-A1 genes encoded a polypeptide of 428 aminoacid residues, with one residue longer than those encoded by Psy1-D1 genes. The Psy1-B1/Psy1-S1 genes encoded four types of polypeptides, with 421 (Psy1-B1h through Psy1-B1j, Psy1-B1l), 427 (Psy1-B1k, Psy1-S1a and Psy1-S1c), 428 (Psy1-B1m), and 429 (Psy1-S1b) aminoacid residues, respectively. Neighbor joining tree was generated based on the gene sequences of the 27 novel alleles and those of the 13 alleles reported previously in common wheat and its relatives. The phylogenetic tree consisted of two subtrees. The subtree I comprised 11 of 14 alleles at Psy1-A1 locus, nine of 16 alleles at Psy1-B1/Psy1-S1 locus, and ten novel alleles at Psy1-D1 locus, while the subtree II included the other three alleles at Psy1-A1 locus, the remaining four Psy1-B1 alleles and three Psy1-S1 alleles. The alleles from different clusters showed high sequence divergences, indicated by various SNPs and InDels (insertion/deletion). The phylogenetic relationships of these allelic variants at the three loci in common wheat and its relatives also supported the hypothesis that common wheat was originated by recurrent hexaploidization events. In addition, 193 Chinese wheat cultivars with different yellow pigment contents were genotyped with two novel STS markers YP7D-1 and YP7D-2. The results indicated that 191 cultivars contained the allele Psy1-D1a, and two had Psy1-D1g.

Recruitment of closely linked genes for divergent functions: the seed storage protein (Glu-3) and powdery mildew (Pm3) genes in wheat (Triticum aestivum L.)

Wheat seed storage protein gene loci (Glu-3) and powdery mildew resistance gene loci (Pm3 and Pm3-like) are closely linked on the short arms of homoeologous group 1 chromosomes. To study the structural organization of the Glu-3/Pm3 loci, three bacterial artificial chromosome clones were sequenced from the A, B, and D genomes of hexaploid wheat. The A and B genome clones contained a Glu-3 adjacent to a Pm3-like gene organized in a conserved Glu-3/SFR159/Pm3-like structure. The D genome clone contained clusters of resistance gene analogs but no Pm3. Its similarity to the A and B genome was limited to the Glu-3/SFR159 region. Comparison of the B genome PM3-like deduced amino acid sequence with known PM3 functional isotypes reinforced the hypothesis of allelic evolution via block exchange by gene conversion/recombination. The advent of glutenin genes and the formation of the Glu-3/SFR159/Pm3 locus occurred after divergence of wheat from rice and Brachypodium. Comparison of the A genome homologous sequences permitted an estimate of time of divergence of approximately 0.3 million years ago. The B genome sequences were not colinear indicating that they could either be paralogs or represent different B genome progenitors. Analysis of the 11 complete retrotransposons indicated a time of divergence ranging from 0.29 to 5.62 million years ago, consistent with their complex nested structure.

Two different CC-NBS-LRR genes are required for Lr10-mediated leaf rust resistance in tetraploid and hexaploid wheat

Comparative study of disease resistance genes in crop plants and their relatives provides insight on resistance gene function, evolution and diversity. Here, we studied the allelic diversity of the Lr10 leaf rust resistance gene, a CC-NBS-LRR coding gene originally isolated from hexaploid wheat, in 20 diploid and tetraploid wheat lines. Besides a gene in the tetraploid wheat variety 'Altar' that is identical to the hexaploid wheat Lr10, two additional, functional resistance alleles showing sequence diversity were identified by virus-induced gene silencing in tetraploid wheat lines. In contrast to most described NBS-LRR proteins, the N-terminal CC domain of LR10 was found to be under strong diversifying selection. A second NBS-LRR gene at the Lr10 locus, RGA2, was shown through silencing to be essential for Lr10 function. Interestingly, RGA2 showed much less sequence diversity than Lr10. These data demonstrate allelic diversity of functional genes at the Lr10 locus in tetraploid wheat, and these new genes can now be analyzed for agronomic relevance. Lr10-based resistance is highly unusual both in its dependence on two, only distantly, related CC-NBS-LRR proteins, as well as in the pattern of diversifying selection in the N-terminal domain. This indicates a new and complex molecular mechanism of pathogen detection and signal transduction.

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.

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.

New insights into the origin of the B genome of hexaploid wheat: evolutionary relationships at the SPA genomic region with the S genome of the diploid relative Aegilops speltoides

Background

Several studies suggested that the diploid ancestor of the B genome of tetraploid and hexaploid wheat species belongs to the Sitopsis section, having Aegilops speltoides (SS, 2n = 14) as the closest identified relative. However molecular relationships based on genomic sequence comparison, including both coding and non-coding DNA, have never been investigated. In an attempt to clarify these relationships, we compared, in this study, sequences of the Storage Protein Activator (SPA) locus region of the S genome of Ae. speltoides (2n = 14) to that of the A, B and D genomes co-resident in the hexaploid wheat species (Triticum aestivum, AABBDD, 2n = 42).

Results

Four BAC clones, spanning the SPA locus of respectively the A, B, D and S genomes, were isolated and sequenced. Orthologous genomic regions were identified as delimited by shared non-transposable elements and non-coding sequences surrounding the SPA gene and correspond to 35 268, 22 739, 43 397 and 53 919 bp for the A, B, D and S genomes, respectively. Sequence length discrepancies within and outside the SPA orthologous regions are the result of non-shared transposable elements (TE) insertions, all of which inserted after the progenitors of the four genomes divergence.

Conclusion

On the basis of conserved sequence length as well as identity of the shared non-TE regions and the SPA coding sequence, Ae speltoides appears to be more evolutionary related to the B genome of T. aestivum than the A and D genomes. However, the differential insertions of TEs, none of which are conserved between the two genomes led to the conclusion that the S genome of Ae. speltoides has diverged very early from the progenitor of the B genome which remains to be identified.

New insights into the origin of the B genome of hexaploid wheat: evolutionary relationships at the SPA genomic region with the S genome of the diploid relative Aegilops speltoides

BACKGROUND

Several studies suggested that the diploid ancestor of the B genome of tetraploid and hexaploid wheat species belongs to the Sitopsis section, having Aegilops speltoides (SS, 2n = 14) as the closest identified relative. However molecular relationships based on genomic sequence comparison, including both coding and non-coding DNA, have never been investigated. In an attempt to clarify these relationships, we compared, in this study, sequences of the Storage Protein Activator (SPA) locus region of the S genome of Ae. speltoides (2n = 14) to that of the A, B and D genomes co-resident in the hexaploid wheat species (Triticum aestivum, AABBDD, 2n = 42).

RESULTS

Four BAC clones, spanning the SPA locus of respectively the A, B, D and S genomes, were isolated and sequenced. Orthologous genomic regions were identified as delimited by shared non-transposable elements and non-coding sequences surrounding the SPA gene and correspond to 35,268, 22,739, 43,397 and 53,919 bp for the A, B, D and S genomes, respectively. Sequence length discrepancies within and outside the SPA orthologous regions are the result of non-shared transposable elements (TE) insertions, all of which inserted after the progenitors of the four genomes divergence.

CONCLUSION

On the basis of conserved sequence length as well as identity of the shared non-TE regions and the SPA coding sequence, Ae speltoides appears to be more evolutionary related to the B genome of T. aestivum than the A and D genomes. However, the differential insertions of TEs, none of which are conserved between the two genomes led to the conclusion that the S genome of Ae. speltoides has diverged very early from the progenitor of the B genome which remains to be identified.

Genome organisation and retrotransposon driven molecular evolution of the endosperm Hardness (Ha) locus in Triticum aestivum cv Glenlea

Wheat endosperm texture is controlled primarily by a locus (Ha), which comprises Gsp-1, Pina and Pinb genes encoding the so-called grain softness protein, puroindoline-a and puroindoline-b, respectively. Pina and Pinb were detected only on the D-genome of hexaploid wheat and its diploid progenitors while Gsp-1 was on all three homoeologous loci. Hexaploid cultivar Glenlea has a hard phenotype due to a null Pina genotype (D-genome) but the sequence organization is not reported. This study aimed at understanding the evolution of homoeologous Ha loci. Sequencing of three BAC clones from cv Glenlea was performed and sequence analyses delimited the Ha loci which spanned 3,925, 5,330 and 31,607 bp in the A-, B- and D-genomes, respectively. A solo LTR of Angela retroelement, downstream to Gsp-A1 and a fragment of Sabrina retroelement, downstream of Gsp-B1, were discovered. We propose that the insertion of these elements into the intergenic regions have driven the deletions of genomic segments harbouring Pina and Pinb genes in the A- and B-genomes of hexaploid wheat. Similarly, fragments of Romani and Vagabond retroelements were identified between truncated Pina and Pinb genes, indicating their role in the deletion of Pina in Glenlea, leading to its hard texture. Structural differences of the Ha locus region of the A-genome between two hexaploid wheat varieties namely Glenlea and Renan (CR626929), suggested the presence of more than one tetraploid ancestor in the origin of hexaploid wheat.

Dynamics and differential proliferation of transposable elements during the evolution of the B and A genomes of wheat

Transposable elements (TEs) constitute >80% of the wheat genome but their dynamics and contribution to size variation and evolution of wheat genomes (Triticum and Aegilops species) remain unexplored. In this study, 10 genomic regions have been sequenced from wheat chromosome 3B and used to constitute, along with all publicly available genomic sequences of wheat, 1.98 Mb of sequence (from 13 BAC clones) of the wheat B genome and 3.63 Mb of sequence (from 19 BAC clones) of the wheat A genome. Analysis of TE sequence proportions (as percentages), ratios of complete to truncated copies, and estimation of insertion dates of class I retrotransposons showed that specific types of TEs have undergone waves of differential proliferation in the B and A genomes of wheat. While both genomes show similar rates and relatively ancient proliferation periods for the Athila retrotransposons, the Copia retrotransposons proliferated more recently in the A genome whereas Gypsy retrotransposon proliferation is more recent in the B genome. It was possible to estimate for the first time the proliferation periods of the abundant CACTA class II DNA transposons, relative to that of the three main retrotransposon superfamilies. Proliferation of these TEs started prior to and overlapped with that of the Athila retrotransposons in both genomes. However, they also proliferated during the same periods as Gypsy and Copia retrotransposons in the A genome, but not in the B genome. As estimated from their insertion dates and confirmed by PCR-based tracing analysis, the majority of differential proliferation of TEs in B and A genomes of wheat (87 and 83%, respectively), leading to rapid sequence divergence, occurred prior to the allotetraploidization event that brought them together in Triticum turgidum and Triticum aestivum, <0.5 million years ago. More importantly, the allotetraploidization event appears to have neither enhanced nor repressed retrotranspositions. We discuss the apparent proliferation of TEs as resulting from their insertion, removal, and/or combinations of both evolutionary forces.

Lgals6, a 2-million-year-old gene in mice: a case of positive Darwinian selection and presence/absence polymorphism

Duplications of genes are widely considered to be a driving force in the evolutionary process. The fate of such duplicated genes (paralogs) depends mainly on the early stages of their evolution. Therefore, the study of duplications that have already started to diverge is useful to better understand their evolution. We present here the example of a 2-million-year-old segmental duplication at the origin of the Lgals4 and Lgals6 genes in the mouse genome. We analyzed the distribution of these genes in samples from 110 wild individuals and wild-derived inbred strains belonging to eight mouse species from Mus (Coelomys) pahari to M. musculus and 28 laboratory strains. Using a maximum-likelihood method, we show that the sequence of the Lgals6 gene has evolved under the influence of strong positive selection that is likely to result in its neofunctionalization. Surprisingly, despite this selection pressure, the Lgals6 gene is present in some mouse species, but not all. Furthermore, even within the species and populations where it is present, the Lgals6 gene is never fixed. To explain this paradox, we propose different hypotheses such as balanced selection and neutral retention of ancient polymophism and we discuss this unexpected result with regard to known galectin properties and response to infections by pathogens.

Evolutionary origin of the segmental duplication encompassing the wheat GLU-B1 locus encoding the overexpressed Bx7 (Bx7OE) high molecular weight glutenin subunit

Sequencing of a BAC clone encompassing the Glu-B1 locus in Glenlea, revealed a 10.3 Kb segmental duplication including the Bx7 gene and flanking an LTR retroelement. To better understand the evolution of this locus, two collections of wheat were surveyed. The first consisted of 96 diploid and tetraploid species accessions while the second consisted of 316 Triticum aestivum cultivars and landraces from 41 countries. The genotypes were first characterized by SDS-PAGE and a total of 40 of the 316 T. aestivum accessions were found to display the overexpressed Bx7 phenotype (Bx7OE). Three lines from the 96 diploid/tetraploid collection also displayed the stronger intensity staining characteristic of the Bx7(OE) subunit. The relative amounts of the Bx7 subunit to total HMW-GS were quantified by RP-HPLC for all Bx7OE accessions and a number of checks. The entire collection was assessed for the presence of four DNA markers namely an 18 bp indel of the coding region of Bx7 variant alleles, a 43 bp indel of the 5'-region and the left and right junctions of the LTR retrotransposon borders and the duplicated segment. All 43 accessions found to have the Bx7OE subunit by SDS-PAGE and RP-HPLC produced the four diagnostic PCR amplicons. None of the lines without the Bx7OE had the LTR retroelement/duplication genomic structure. However, the 18 and 43 bp indel were found in accessions other than Bx7OE. These results indicate that the overexpression of the Bx7 HMW-GS is likely the result of a single event, i.e., a gene duplication at the Glu-B1 locus mediated by the insertion of a retroelement. Also, the 18 and 43 bp indels pre-date the duplication event. Allelic variants Bx7*, Bx7 with and without 43 bp insert and Bx7OE were found in both tetraploid and hexaploid collections and shared the same genomic organization. Though the possibility of introgression from T. aestivum to T. turgidum cannot be ruled out, the three structural genomic changes of the B-genome taken together support the hypothesis of multiple polyploidization events involving different tetraploid progenitors.

Rapid evolution and complex structural organization in genomic regions harboring multiple prolamin genes in the polyploid wheat genome

Genes encoding wheat prolamins belong to complicated multi-gene families in the wheat genome. To understand the structural complexity of storage protein loci, we sequenced and analyzed orthologous regions containing both gliadin and LMW-glutenin genes from the A and B genomes of a tetraploid wheat species, Triticum turgidum ssp. durum. Despite their physical proximity to one another, the gliadin genes and LMW-glutenin genes are organized quite differently. The gliadin genes are found to be more clustered than the LMW-glutenin genes which are separated from each other by much larger distances. The separation of the LMW-glutenin genes is the result of both the insertion of large blocks of repetitive DNA owing to the rapid amplification of retrotransposons and the presence of genetic loci interspersed between them. Sequence comparisons of the orthologous regions reveal that gene movement could be one of the major factors contributing to the violation of microcolinearity between the homoeologous A and B genomes in wheat. The rapid sequence rearrangements and differential insertion of repetitive DNA has caused the gene islands to be not conserved in compared regions. In addition, we demonstrated that the i-type LMW-glutenin originated from a deletion of 33-bps in the 5' coding region of the m-type gene. Our results show that multiple rounds of segmental duplication of prolamin genes have driven the amplification of the omega-gliadin genes in the region; such segmental duplication could greatly increase the repetitive DNA content in the genome depending on the amount of repetitive DNA present in the original duplicate region.

Contrasting rates of evolution in Pm3 loci from three wheat species and rice

The Pm3 gene from wheat confers resistance against powdery mildew and recent studies have shown that it is a member of a multigene family in the wheat genome. We compared genomic sequences ranging from 178 to 332 kb containing six Pm3-like genes and five gene fragments from orthologous loci in the A genome of wheat at three different ploidy levels. We found that the wheat Pm3 loci display an extremely dynamic evolution where sequence conservation is minimal between species and basically limited to very short sequences containing the genetic markers that define the orthology. The Pm3-like genes and their up- and downstream regions were reshuffled by multiple rearrangements, resulting in a complex mosaic of conserved and unique sequences. Comparison with rice showed that the known wheat Pm3-like genes represent only one branch of a large superfamily with several clusters in rice and suggests the presence of additional similar genes in the wheat genome. Estimates of divergence times and transposable-element insertions indicate that the Pm3 locus in wheat has undergone more drastic changes in its recent evolution than its counterpart in rice. This indicates that loci containing homologous resistance gene analogs can evolve at highly variable speeds in different species.

New insights into Oryza genome evolution: high gene colinearity and differential retrotransposon amplification

An approximately 247-kb genomic region from FF genome of wild rice Oryza brachyantha, possessing the smallest Oryza genome, was compared to the orthologous approximately 450-kb region from AA genome, O. sativa L. ssp. japonica. 37 of 38 genes in the orthologous regions are shared between japonica and O. brachyantha. Analyses of nucleotide substitution in coding regions suggest the two genomes diverged approximately 10 million years ago. Comparisons of transposable elements (TEs) reveal that the density of DNA TEs in O. brachyantha is comparable to O. sativa; however, the density of RNA TEs is dramatically lower. The genomic fraction of RNA TEs in japonica is two times greater than in O. brachyantha. Differences, particularly in RNA TEs, in this region and in BAC end sequences from five wild and two cultivated Oryza species explain major genome size differences between sativa and brachyantha. Gene expression analyses of three ObDREB1 genes in the sequenced region indicate orthologous genes retain similar expression patterns following cold stress. Our results demonstrate that size and number of RNA TEs play a major role in genomic differentiation and evolution in Oryza. Additionally, distantly related O. brachyantha shares colinearity with O. sativa, offering opportunities to use comparative genomics to explore the genetic diversity of wild species to improve cultivated rice.

Grinding up wheat: a massive loss of nucleotide diversity since domestication

Several demographic and selective events occurred during the domestication of wheat from the allotetraploid wild emmer (Triticum turgidum ssp. dicoccoides). Cultivated wheat has since been affected by other historical events. We analyzed nucleotide diversity at 21 loci in a sample of 101 individuals representing 4 taxa corresponding to representative steps in the recent evolution of wheat (wild, domesticated, cultivated durum, and bread wheats) to unravel the evolutionary history of cultivated wheats and to quantify its impact on genetic diversity. Sequence relationships are consistent with a single domestication event and identify 2 genetically different groups of bread wheat. The wild group is not highly polymorphic, with only 212 polymorphic sites among the 21,720 bp sequenced, and, during domestication, diversity was further reduced in cultivated forms--by 69% in bread wheat and 84% in durum wheat--with considerable differences between loci, some retaining no polymorphism at all. Coalescent simulations were performed and compared with our data to estimate the intensity of the bottlenecks associated with domestication and subsequent selection. Based on our 21-locus analysis, the average intensity of domestication bottleneck was estimated at about 3--giving a population size for the domesticated form about one third that of wild dicoccoides. The most severe bottleneck, with an intensity of about 6, occurred in the evolution of durum wheat. We investigated whether some of the genes departed from the empirical distribution of most loci, suggesting that they might have been selected during domestication or breeding. We detected a departure from the null model of demographic bottleneck for the hypothetical gene HgA. However, the atypical pattern of polymorphism at this locus might reveal selection on the linked locus Gsp1A, which may affect grain softness--an important trait for end-use quality in wheat.

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.

In vivo participation of red chlorophyll catabolite reductase in chlorophyll breakdown

A central reaction of chlorophyll breakdown, porphyrin ring opening of pheophorbide a to the primary fluorescent chlorophyll catabolite (pFCC), requires pheophorbide a oxygenase (PAO) and red chlorophyll catabolite reductase (RCCR), with red chlorophyll catabolite (RCC) as a presumably PAO-bound intermediate. In subsequent steps, pFCC is converted to different fluorescent chlorophyll catabolites (FCCs) and nonfluorescent chlorophyll catabolites (NCCs). Here, we show that RCCR-deficient Arabidopsis thaliana accumulates RCC and three RCC-like pigments during senescence, as well as FCCs and NCCs. We also show that the stereospecificity of Arabidopsis RCCR is defined by a small protein domain and can be reversed by a single Phe-to-Val exchange. Exploiting this feature, we prove the in vivo participation of RCCR in chlorophyll breakdown. After complementation of RCCR mutants with RCCRs exhibiting alternative specificities, patterns of chlorophyll catabolites followed the specificity of complementing RCCRs. Light-dependent leaf cell death observed in different RCCR-deficient lines strictly correlated with the accumulation of RCCs and the release of singlet oxygen, and PAO induction preceded lesion formation. These findings suggest that RCCR absence causes leaf cell death as a result of the accumulation of photodynamic RCC. We conclude that RCCR (together with PAO) is required for the detoxification of chlorophyll catabolites and discuss the biochemical role(s) for this enzyme.

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.

Development of simple sequence repeat markers specific for the Lr34 resistance region of wheat using sequence information from rice and Aegilops tauschii

Hexaploid wheat (Triticum aestivum L.) originated about 8,000 years ago from the hybridization of tetraploid wheat with diploid Aegilops tauschii Coss. containing the D-genome. Thus, the bread wheat D-genome is evolutionary young and shows a low degree of polymorphism in the bread wheat gene pool. To increase marker density around the durable leaf rust resistance gene Lr34 located on chromosome 7DS, we used molecular information from the orthologous region in rice. Wheat expressed sequence tags (wESTs) were identified by homology with the rice genes in the interval of interest, but were monomorphic in the 'Arina' x 'Forno' mapping population. To derive new polymorphic markers, bacterial artificial chromosome (BAC) clones representing a total physical size of approximately 1 Mb and belonging to four contigs were isolated from Ae. tauschii by hybridization screening with wheat ESTs. Several BAC clones were low-pass sequenced, resulting in a total of approximately 560 kb of sequence. Ten microsatellite sequences were found, and three of them were polymorphic in our population and were genetically mapped close to Lr34. Comparative analysis of marker order revealed a large inversion between the rice genome and the wheat D-genome. The SWM10 microsatellite is closely linked to Lr34 and has the same allele in the three independent sources of Lr34: 'Frontana', 'Chinese Spring', and 'Forno', as well in most of the genotypes containing Lr34. Therefore, SWM10 is a highly useful marker to assist selection for Lr34 in breeding programs worldwide.

Rapid generation of new powdery mildew resistance genes after wheat domestication

Plant defence against pathogens is controlled by disease resistance (R) gene products that directly or indirectly detect specific pathogen effectors. Plant-pathogen interactions have been proposed to follow a co-evolutionary arms-race model where R genes are recent and evolve rapidly in response to structural changes in matching pathogen effectors. However, the longevity and extensive polymorphism of R genes studied were more consistent with balancing selection maintaining ancient and diverse R genes or alleles. In bread wheat (Triticum aestivum), the Pm3 locus confers race-specific resistance to wheat powdery mildew (Blumeria graminis f.sp. triticii). Here we describe recently generated Pm3 resistance alleles that all derive from one susceptible allele, Pm3CS, which is widespread among hexaploid bread-wheat lines. One group of four Pm3 resistance alleles shows few, clearly delimited, polymorphic sequence blocks of ancient origin, embedded in sequences identical to Pm3CS and possibly derived from gene conversion. A second group of three alleles differs from Pm3CS by only two to five mutations, all non-synonymous, and all in the leucine-rich repeat-encoding region. Transient transformation experiments confirmed that Pm3 resistance specificities are based on one or few amino acid changes. The Pm3CS allele was found in wild tetraploid wheat, the ancestor of hexaploid bread wheat, specifically from southern Turkey, a region proposed to be the site of wheat domestication. Based on these data, we propose that the Pm3 resistance alleles were generated in agricultural ecosystems after domestication of wheat 10,000 years ago. The evolution of Pm3 alleles in wheat is best described by the model of evolved recycling, where novel genetic variation is integrated in plant populations together with recycling of old variation.

Mechanisms of genomic rearrangements and gene expression changes in plant polyploids

Polyploidy is produced by multiplication of a single genome (autopolyploid) or combination of two or more divergent genomes (allopolyploid). The available data obtained from the study of synthetic (newly created or human-made) plant allopolyploids have documented dynamic and stochastic changes in genomic organization and gene expression, including sequence elimination, inter-chromosomal exchanges, cytosine methylation, gene repression, novel activation, genetic dominance, subfunctionalization and transposon activation. The underlying mechanisms for these alterations are poorly understood. To promote a better understanding of genomic and gene expression changes in polyploidy, we briefly review origins and forms of polyploidy and summarize what has been learned from genome-wide gene expression analyses in newly synthesized auto-and allopolyploids. We show transcriptome divergence between the progenitors and in the newly formed allopolyploids. We propose models for transcriptional regulation, chromatin modification and RNA-mediated pathways in establishing locus-specific expression of orthologous and homoeologous genes during allopolyploid formation and evolution.

Single nucleotide polymorphism, genetic mapping, and expression of genes coding for the DOF wheat prolamin-box binding factor

Wheat prolamin-box binding factor (WPBF) was shown to be an activator of Triticum aestivum L. storage protein genes. Three homoeologous genes encoding this transcription factor were isolated from a bacterial artificial chromosome genomic library and sequenced. The genes all have two exons separated by an intron of approximately 1,000 bp where the second exon contains the entire coding sequence. Many differences were found between homoeologous sequences, but none of them is predicted to significantly alter the sequence of the putative encoded protein. The three homoeologous genes are specifically expressed in grain from 3 to 39 days after anthesis. The allelic variation of a genetically diverse collection of 27 bread wheat lines was assessed. One, five, and one single-nucleotide polymorphisms (SNPs) were detected in the wPbf genes for the A, B, and D genomes, respectively. Physical and genetic mapping utilizing some of the SNPs identified confirmed that wPbf genes are located close to the centromeres on the homoeologous group 5 chromosomes. The low level of allelic diversity found in wPbf genes may suggest that these genes play a key role and are thus constrained by selection.

Structural variation and evolution of a defense-gene cluster in natural populations of Aegilops tauschii

Genetic mapping and sequencing of plant genomes have been useful for investigating eukaryotic chromosome structural organization. In many cases, analyses have been limited in the number of representatives sampled from specific groups. The degree of intraspecific genome diversity remains in question. The possibility exists that a single model genome may have limited utility for identifying genes in related members of the species or genus. Crop improvement programs have particular interests in disease resistance genes that are harbored by wild relatives of modern cultivated crops. These genes are evolutionarily dynamic and under selective pressure by a broad range of pathogenic organisms. Using resistance gene analogs as models for gene evolution, intraspecific genome comparisons were made among populations of wild diploid wheat (Aegilops tauschii). We observed that deletion haplotypes are occurring frequently and independently in the genome. Haplotypes are geographically correlated and maintenance of gene complements in localized populations indicates selective advantage. Furthermore, deletion haplotypes are not detrimental to plant health, since genes without adaptive value in alternate environments are eliminated from the genome. Deletion haplotypes appear to be a common form of allelic variation in plants, and we address the consequences on genome restructuring and gene evolution.

Alpha-gliadin genes from the A, B, and D genomes of wheat contain different sets of celiac disease epitopes

Background

Bread wheat (Triticum aestivum) is an important staple food. However, wheat gluten proteins cause celiac disease (CD) in 0.5 to 1% of the general population. Among these proteins, the α-gliadins contain several peptides that are associated to the disease.

Results

We obtained 230 distinct α-gliadin gene sequences from severaldiploid wheat species representing the ancestral A, B, and D genomes of the hexaploid bread wheat. The large majority of these sequences (87%) contained an internal stop codon. All α-gliadin sequences could be distinguished according to the genome of origin on the basis of sequence similarity, of the average length of the polyglutamine repeats, and of the differences in the presence of four peptides that have been identified as T cell stimulatory epitopes in CD patients through binding to HLA-DQ2/8. By sequence similarity, α-gliadins from the public database of hexaploid T. aestivum could be assigned directly to chromosome 6A, 6B, or 6D. T. monococcum (A genome) sequences, as well as those from chromosome 6A of bread wheat, almost invariably contained epitope glia-α9 and glia-α20, but never the intact epitopes glia-α and glia-α2. A number of sequences from T. speltoides, as well as a number of sequences fromchromosome 6B of bread wheat, did not contain any of the four T cell epitopes screened for. The sequences from T. tauschii (D genome), as well as those from chromosome 6D of bread wheat, were found to contain all of these T cell epitopes in variable combinations per gene. The differences in epitope composition resulted mainly from point mutations. These substitutions appeared to be genome specific.

Conclusion

Our analysis shows that α-gliadin sequences from the three genomes of bread wheat form distinct groups. The four known T cell stimulatory epitopes are distributed non-randomly across the sequences, indicating that the three genomes contribute differently to epitope content. A systematic analysis of all known epitopes in gliadins and glutenins will lead to better understanding of the differences in toxicity among wheat varieties. On the basis of such insight, breeding strategies can be designed to generate less toxic varieties of wheat which may be tolerated by at least part of the CD patient population.