Molecular basis of evolutionary events that shaped the hardness locus in diploid and polyploid wheat species (Triticum and Aegilops)

Template switching can create complex LTR retrotransposon insertions in Triticeae genomes

Background

The LTR (long terminal repeat) retrotransposons of higher plants are replicated by a mutagenic life cycle containing transcription and reverse transcription steps. The DNA copies are often subject to recombination once integrated into the genome. Complex elements, where two elements share an LTR, are not uncommon. They are thought to result from heterologous recombination between two adjacent elements that occurs following their integration.

Results

Here, we present evidence for another potential mechanism for the creation of complex elements, involving abnormal template switching during reverse transcription. The template switching creates a large, complex daughter element, formed by the fusion of two parent sequences, which is then inserted into the genome.

Conclusion

Those complex elements are part of the genome structure of plants in the Poaceae, especially in the Triticeae, but not of Arabidopsis. Hence, retrotransposon dynamics shaping the genome are lineage-specific.

Genome plasticity a key factor in the success of polyploid wheat under domestication

Wheat was domesticated about 10,000 years ago and has since spread worldwide to become one of the major crops. Its adaptability to diverse environments and end uses is surprising given the diversity bottlenecks expected from recent domestication and polyploid speciation events. Wheat compensates for these bottlenecks by capturing part of the genetic diversity of its progenitors and by generating new diversity at a relatively fast pace. Frequent gene deletions and disruptions generated by a fast replacement rate of repetitive sequences are buffered by the polyploid nature of wheat, resulting in subtle dosage effects on which selection can operate.

Differential impact of retrotransposon populations on the genome of allotetraploid tobacco (Nicotiana tabacum)

LTR-retrotransposons contribute substantially to the structural diversity of plant genomes. Recent models of genome evolution suggest that retrotransposon amplification is offset by removal of retrotransposon sequences, leading to a turnover of retrotransposon populations. While bursts of amplification have been documented, it is not known whether removal of retrotransposon sequences occurs continuously, or is triggered by specific stimuli over short evolutionary periods. In this work, we have characterized the evolutionary dynamics of four populations of copia-type retrotransposons in allotetraploid tobacco (Nicotiana tabacum) and its two diploid progenitors Nicotiana sylvestris and Nicotiana tomentosiformis. We have used SSAP (Sequence-Specific Amplification Polymorphism) to evaluate the contribution retrotransposons have made to the diversity of tobacco and its diploid progenitor species, to quantify the contribution each diploid progenitor has made to tobacco's retrotransposon populations, and to estimate losses or amplifications of retrotransposon sequences subsequent to tobacco's formation. Our results show that the tobacco genome derives from a turnover of retrotransposon sequences with removals concomitant with new insertions. We have detected unique behaviour specific to each retrotransposon population, with differences likely reflecting distinct evolutionary histories and activities of particular elements. Our results indicate that the retrotransposon content of a given plant species is strongly influenced by the host evolutionary history, with periods of rapid turnover of retrotransposon sequences stimulated by allopolyploidy.

Genetic and epigenetic alteration among three homoeologous genes of a class E MADS box gene in hexaploid wheat

Bread wheat (Triticum aestivum) is a hexaploid species with A, B, and D ancestral genomes. Most bread wheat genes are present in the genome as triplicated homoeologous genes (homoeologs) derived from the ancestral species. Here, we report that both genetic and epigenetic alterations have occurred in the homoeologs of a wheat class E MADS box gene. Two class E genes are identified in wheat, wheat SEPALLATA (WSEP) and wheat LEAFY HULL STERILE1 (WLHS1), which are homologs of Os MADS45 and Os MADS1 in rice (Oryza sativa), respectively. The three wheat homoeologs of WSEP showed similar genomic structures and expression profiles. By contrast, the three homoeologs of WLHS1 showed genetic and epigenetic alterations. The A genome WLHS1 homoeolog (WLHS1-A) had a structural alteration that contained a large novel sequence in place of the K domain sequence. A yeast two-hybrid analysis and a transgenic experiment indicated that the WLHS1-A protein had no apparent function. The B and D genome homoeologs, WLHS1-B and WLHS1-D, respectively, had an intact MADS box gene structure, but WLHS1-B was predominantly silenced by cytosine methylation. Consequently, of the three WLHS1 homoeologs, only WLHS1-D functions in hexaploid wheat. This is a situation where three homoeologs are differentially regulated by genetic and epigenetic mechanisms.

Microcolinearity and genome evolution in the AdhA region of diploid and polyploid cotton (Gossypium)

Genome sizes vary by several orders of magnitude, driven by mechanisms such as illegitimate recombination and transposable element proliferation. Prior analysis of the CesA region in two cotton genomes that diverged 5-10 million years ago (Ma), and acquired a twofold difference in genome size, revealed extensive local conservation of genic and intergenic regions, with no evidence of the global genome size difference. The present study extends the comparison to include BAC sequences surrounding the gene encoding alcohol dehydrogenase A (AdhA) from four cotton genomes: the two co-resident genomes (A(T) and D(T)) of the allotetraploid, Gossypium hirsutum, as well as the model diploid progenitors, Gossypium arboreum (A) and Gossypium raimondii (D). In contrast to earlier work, evolution in the AdhA region reflects, in a microcosm, the overall difference in genome size, with a nearly twofold difference in aligned sequence length. Most size differences may be attributed to differential accumulation of retroelements during divergence of the genome diploids from their common ancestor, but in addition there has been a biased accumulation of small deletions, such that those in the smaller D genome are on average twice as large as those in the larger A genome. The data also provide evidence for the global phenomenon of 'genomic downsizing' in polyploids shortly after formation. This in part reflects a higher frequency of small deletions post-polyploidization, and increased illegitimate recombination. In conjunction with previous work, the data here confirm the conclusion that genome size evolution reflects many forces that collectively operate heterogeneously among genomic regions.

Orthologous comparison in a gene-rich region among grasses reveals stability in the sugarcane polyploid genome

Modern sugarcane (Saccharum spp.) is an important grass that contributes 60% of the raw sugar produced worldwide and has a high biofuel production potential. It was created about a century ago through hybridization of two highly polyploid species, namely S. officinarum and S. spontaneum. We investigated genome dynamics in this highly polyploid context by analyzing two homoeologous sequences (97 and 126 kb) in a region that has already been studied in several cereals. Our findings indicate that the two Saccharum species diverged 1.5-2 million years ago from one another and 8-9 million years ago from sorghum. The two sugarcane homoeologous haplotypes show perfect colinearity as well as high gene structure conservation. Apart from the insertion of a few retrotransposable elements, high homology was also observed for the non-transcribed regions. Relative to sorghum, the sugarcane sequences displayed colinearity, with the exception of two genes present only in sorghum, and striking homology in most non-coding parts of the genome. The gene distribution highlighted high synteny and colinearity with rice, and partial colinearity with each homoeologous maize region, which became perfect when the sequences were combined. The haplotypes observed in sugarcane may thus closely represent the ancestral Andropogoneae haplotype. This analysis of sugarcane haplotype organization at the sequence level suggests that the high ploidy in sugarcane did not induce generalized reshaping of its genome, thus challenging the idea that polyploidy quickly induces generalized rearrangement of genomes. These results also confirm the view that sorghum is the model of choice for sugarcane.

Promoter analysis and immunolocalisation show that puroindoline genes are exclusively expressed in starchy endosperm cells of wheat grain

The purolindolines are small cysteine-rich proteins which are present in the grain of wheat. They have a major impact on the utilisation of the grain as they are the major determinants of grain texture, which affects both milling and baking properties. Bread and durum wheats were transformed with constructs comprising the promoter regions of the Puroindoline a (Pina) and Puroindoline b (Pinb) genes fused to the uidA (GUS) reporter gene. Nine lines showing 3:1 segregation for the transgene and comprising all transgene/species combinations were selected for detailed analysis of transgene expression during grain development. This showed that transgene expression occurred only in the starchy endosperm cells and was not observed in any other seed or vegetative tissues. The location of the puroindoline proteins in these cells was confirmed by tissue printing of developing grain, using a highly specific monoclonal antibody for detection and an antibody to the aleurone-localised 8S globulin as a control. This provides clear evidence that puroindolines are only synthesised and accumulated in the starchy endosperm cells of the wheat grain.

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.

Characterization of the gene Mre11 and evidence of silencing after polyploidization in Triticum

The MRE11 protein is a component of the highly conserved MRN complex, along with RAD50 and NBS1. This complex is crucial in the repair of breaks in double stranded DNA, and is involved in many other cell processes. The present paper reports the molecular characterization of Mre11 gene in all three genomes of wheat, making use of the diploid species Triticum monococcum (genome A) and Aegilops Tauschii (genome D), the tetraploid T. turgidum (genomes A and B), and the hexaploid T. aestivum (genomes A, B and D). The genomic sequences characterized ranged from 4,662 to 4,766 bp in length; the cDNA corresponding to the processed mRNA was 2,440-2,510 bp long. In all cases, Mre11 coded for a highly conserved protein of 699 amino acids with a structure involving 22 exons. Mre11 expression was determined by real-time PCR in all the species analysed. The tetraploid species showed an expression similar to that of the diploid Ae. tauschii and lower than that of T. monococcum. Stronger expression was detected in the hexaploid T. aestivum. The SSCP technique was modified by introducing fluorescent labelling to the procedure in order to analyse the expression of the different Mre11 genes (i.e., those belonging to the different genomes) in the polyploid species. In both polyploids, the Mre11 gene belonging to the B genome was the least expressed. This probably reflects a first step in the process of silencing duplicate genes after polyploidization.

Measuring global gene expression in polyploidy; a cautionary note from allohexaploid wheat

The number of global gene expression studies has increased significantly in recent years. It is assumed that the different techniques employed report similar levels of gene expression for each sequence type. While this may be true for many species, polyploids containing homoeologous and paralogous gene copies represent a unique situation. In this paper, we describe the comparison of the Affymetrix GeneChip Wheat Genome Array, an in-house custom-spotted complementary DNA array and quantitative reverse transcription-polymerase chain reaction (PCR) for the study of gene expression in hexaploid wheat. Analysis of the data generated from each platform revealed little concordance and suggested that global comparisons are not possible. Potential causes of these inter-platform discrepancies were investigated and revealed to be due to the inability of the platforms to discriminate between different but related transcripts. Our results also showed that the traditionally used array validation technique, quantitative reverse transcription PCR, differs in its discriminatory ability, resulting in the poor confirmation rates seen in previous polyploid studies. These findings have implications for gene expression studies in polyploid organisms and highlight the need for homoeologous- and paralogous-specific arrays when investigating polyploid gene expression.

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.

Capturing diversity in the cereals: many options but little promiscuity

It is generally recognized by geneticists and plant breeders alike that there is a need to further improve the ability to capture and manipulate genetic diversity. The effective harnessing of diversity in traditional breeding programmes is limited and, therefore, it is vital that meiotic recombination can be manipulated given that it plays a pivotal role in generating diversity. With the advent of a wider range of genomics technologies, our understanding of meiotic processes should increase rapidly. Although comparative genetics has been useful, particularly in the broader grass family, the development of physical maps, long-range sequencing and transcript profiles promises to unravel the complexities of genomes as large or larger than wheat. Highlighting the most significant findings to date, this review pools the knowledge on these tools and reproductive processes.

Morgane, a new LTR retrotransposon group, and its subfamilies in wheats

Transposable elements are the main components of grass genomes, especially in Triticeae species. In a previous analysis, we identified a very short element, Morgane_CR626934-1; here we describe more precisely this unusual element. Morgane_CR626934-1 shows high sequence identity (until 98%) with ESTs belonging to other possible small elements, expressed under abiotic and biotic stress conditions. No putative functional polyprotein could be identified in all of these different Morgane-like sequences. Moreover, elements from the Morgane_CR626934-1 subfamily are found only in wheats and Agropyrum genomes and among these species, only Ae. tauschii and T. aestivum present a high copy number of these elements. They are highly conserved in wheat genomes (95.5%). Based on the uncommon characteristics of the described Morgane-like elements, we proposed to classify them in a new group within the Class I LTR retrotransposon, the Morgane group.

The promoter of the wheat puroindoline-a gene (PinA) exhibits a more complex pattern of activity than that of the PinB gene and is induced by wounding and pathogen attack in rice

Puroindolines form the molecular basis of wheat grain hardness. However, little is known about puroindoline gene regulation. We previously reported that the Triticum aestivum puroindoline-b gene (PinB) promoter directs beta-glucuronidase gene (uidA) seed-specific expression in transgenic rice. In this study, we isolated a puroindoline-a gene (PinA), analyzed PinA promoter activity by 5' deletions and compared PinA and PinB promoters in transgenic rice. Seeds of PinA-1214 and PinB-1063 transgenic plants strongly expressed uidA in endosperm, in the aleurone layer and in epidermis cells in a developmentally regulated manner. The GUS activity was also observed in PinA-1214 embryos. Whereas the PinB promoter is seed specific, the PinA promoter also directed, but to a lower level, uidA expression in roots of seedlings and in the vascular tissues of palea and pollen grains of dehiscent anthers during flower development. In addition, the PinA promoter was induced by wounding and by Magnaporthe grisea. By deletion analysis, we showed that the "390-bp" PinA promoter drives the same expression pattern as the "1214-bp" promoter. Moreover, the "214-bp" PinA promoter drives uidA expression solely in pollen grains of dehiscent anthers. The presence of putative cis-regulatory elements that may be related to PinA expression is discussed from an evolutionary point of view. By electrophoretic mobility shift assay, we showed that putative cis-elements (WUN-box, TCA motifs and as-1-like binding sites) whose presence in the PinA promoter may be related to wounding and/or the pathogen response form complexes with nuclear extracts isolated from wounded wheat leaves.

Homeologous recombination plays a major role in chromosome rearrangements that occur during meiosis of Brassica napus haploids

Chromosomal rearrangements can be triggered by recombination between distinct but related regions. Brassica napus (AACC; 2n = 38) is a recent allopolyploid species whose progenitor genomes are widely replicated. In this article, we analyze the extent to which chromosomal rearrangements originate from homeologous recombination during meiosis of haploid B. napus (n = 19) by genotyping progenies of haploid x euploid B. napus with molecular markers. Our study focuses on three pairs of homeologous regions selected for their differing levels of divergence (N1/N11, N3/N13, and N9/N18). We show that a high number of chromosomal rearrangements occur during meiosis of B. napus haploid and are transmitted by first division restitution (FDR)-like unreduced gametes to their progeny; half of the progeny of Darmor-bzh haploids display duplications and/or losses in the chromosomal regions being studied. We demonstrate that half of these rearrangements are due to recombination between regions of primary homeology, which represents a 10- to 100-fold increase compared to the frequency of homeologous recombination measured in euploid lines. Some of the other rearrangements certainly result from recombination between paralogous regions because we observed an average of one to two autosyndetic A-A and/or C-C bivalents at metaphase I of the B. napus haploid. These results are discussed in the context of genome evolution of B. napus.

Evolutionary dynamics of Waxy and the origin of hexaploid Spartina species (Poaceae)

We investigated the evolutionary dynamics of duplicated copies of the granule-bound starch synthase I gene (GBSSI or Waxy) within polyploid Spartina species. Molecular cloning, sequencing, and phylogenetic analyses revealed incongruences between the expected species phylogeny and the inferred gene trees. Some genes within species were more divergent than expected from ploidy level alone, suggesting the existence of paralogous sets of Waxy loci in Spartina. Phylogenetic analyses indicate that this paralogy originated from a duplication that occurred prior to the divergence of Spartina from other Chloridoideae. Gene tree topologies revealed three divergent homoeologous sequences in the hexaploid S. alterniflora that are consistent with the proposal of an allopolyploid origin of the hexaploid clade. Waxy sequences differ in insertion-deletion events in introns, which may be used to diagnose gene copies. Both paralogous and homoeologous coding regions appear to evolving under selective constraints.

Independent Wheat B and G Genome Origins in Outcrossing Aegilops Progenitor Haplotypes

The origin of modern wheats involved alloploidization among related genomes. To determine if Aegilops speltoides was the donor of the B and G genomes in AABB and AAGG tetraploids, we used a 3-tiered approach. Using 70 amplified fragment length polymorphism (AFLP) loci, we sampled molecular diversity among 480 wheat lines from their natural habitats encompassing all S genome Aegilops, the putative progenitors of wheat B and G genomes. Fifty-nine Aegilops representatives for S genome diversity were compared at 375 AFLP loci with diploid, tetraploid, and 11 nulli-tetrasomic Triticum aestivum wheat lines. B genome-specific markers allowed pinning the origin of the B genome to S chromosomes of A. speltoides, while excluding other lineages. The outbreeding nature of A. speltoides influences its molecular diversity and bears upon inferences of B and G genome origins. Haplotypes at nuclear and chloroplast loci ACC1, G6PDH, GPT, PGK1, Q, VRN1, and ndhF for approximately 70 Aegilops and Triticum lines (0.73 Mb sequenced) reveal both B and G genomes of polyploid wheats as unique samples of A. speltoides haplotype diversity. These have been sequestered by the AABB Triticum dicoccoides and AAGG Triticum araraticum lineages during their independent origins.

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.

Doubling genome size without polyploidization: dynamics of retrotransposition-driven genomic expansions in Oryza australiensis, a wild relative of rice

Retrotransposons are the main components of eukaryotic genomes, representing up to 80% of some large plant genomes. These mobile elements transpose via a "copy and paste" mechanism, thus increasing their copy number while active. Their accumulation is now accepted as the main factor of genome size increase in higher eukaryotes, besides polyploidy. However, the dynamics of this process are poorly understood. In this study, we show that Oryza australiensis, a wild relative of the Asian cultivated rice O. sativa, has undergone recent bursts of three LTR-retrotransposon families. This genome has accumulated more than 90,000 retrotransposon copies during the last three million years, leading to a rapid twofold increase of its size. In addition, phenetic analyses of these retrotransposons clearly confirm that the genomic bursts occurred posterior to the radiation of the species. This provides direct evidence of retrotransposon-mediated variation of genome size within a plant genus.

Advanced resources for plant genomics: a BAC library specific for the short arm of wheat chromosome 1B

Common wheat (Triticum aestivum L., 2n = 6x = 42) is a polyploid species possessing one of the largest genomes among the cultivated crops (1C is approximately 17 000 Mb). The presence of three homoeologous genomes (A, B and D), and the prevalence of repetitive DNA make sequencing the wheat genome a daunting task. We have developed a novel 'chromosome arm-based' strategy for wheat genome sequencing to simplify this task; this relies on sub-genomic libraries of large DNA inserts. In this paper, we used a di-telosomic line of wheat to isolate six million copies of the short arm of chromosome 1B (1BS) by flow sorting. Chromosomal DNA was partially digested with HindIII and used to construct an arm-specific BAC library. The library consists of 65 280 clones with an average insert size of 82 kb. Almost half of the library (45%) has inserts larger than 100 kb, while 18% of the inserts range in size between 75 and 100 kb, and 37% are shorter than 75 kb. We estimated the chromosome arm coverage to be 14.5-fold, giving a 99.9% probability of identifying a clone corresponding to any sequence on the short arm of 1B. Each chromosome arm in wheat can be flow sorted from an appropriate cytogenetic stock, and we envisage that the availability of chromosome arm-specific BAC resources in wheat will greatly facilitate the development of ready-to-sequence physical maps and map-based gene cloning.

PLOTREP: a web tool for defragmentation and visual analysis of dispersed genomic repeats

Identification of dispersed or interspersed repeats, most of which are derived from transposons, retrotransposons or retrovirus-like elements, is an important step in genome annotation. Software tools that compare genomic sequences with precompiled repeat reference libraries using sensitive similarity-based methods provide reliable means of finding the positions of fragments homologous to known repeats. However, their output is often incomplete and fragmented owing to the mutations (nucleotide substitutions, deletions or insertions) that can result in considerable divergence from the reference sequence. Merging these fragments to identify the whole region that represents an ancient copy of a mobile element is challenging, particularly if the element is large and suffered multiple deletions or insertions. Here we report PLOTREP, a tool designed to post-process results obtained by sequence similarity search and merge fragments belonging to the same copy of a repeat. The software allows rapid visual inspection of the results using a dot-plot like graphical output. The web implementation of PLOTREP is available at .

Ogre elements--a distinct group of plant Ty3/gypsy-like retrotransposons

Ogre elements are a group of LTR retrotransposons recently discovered in legume plants, where they constitute almost 40% of the genome in some species. They are exceptional in their size (reaching 25 kb) and possess several specific features, including an intron within a polyprotein-coding region, and an extra open reading frame (ORF1) encoding a protein of unknown function located upstream of the gag gene. Although these features make Ogres interesting for further research, identification of additional elements from a broader range of plant taxa has been complicated by the divergence of their sequences, preventing their detection using similarity-based searches. Here we report the results of structure-based computational searches for Ogre elements in available plant genomic sequences, which proved to be more efficient and revealed occurrences of Ogres in three families of dicot plants (Leguminosae, Solanaceae and Salicaceae). In addition, a representative set of 85 elements was retrieved from a model legume species Medicago truncatula. All identified full-length elements were used for comparative analysis, which showed that in spite of only little conservation of their nucleotide sequences, their protein domains were highly conserved, including several regions within ORF1. Further, the elements shared the same functional regions, including a primer binding site complementary to tRNA(arg), a conserved motif within a polypurine tract, and a putative intron between the pro and rt/rh coding domains. These findings, together with analysis of their phylogenetic relationship to other retrotransposons based on similarities of rt domains suggest that Ogre elements from different plant taxa have a common origin and thus constitute a distinct group of Ty3/gypsy retrotransposons.

Molecular evolution of the puroindoline-a, puroindoline-b, and grain softness protein-1 genes in the tribe Triticeae

The genome organization of the Hardness locus in the tribe Triticeae constitutes an excellent model for studying the mechanisms of evolution that played a role in the preservation and potential functional innovations of duplicate genes. Here we applied the nonsynonymous-synonymous rate ratio (d ( N )/d ( S ) or omega) to measure the selective pressures at the paralogous puroindoline-a (Pina), puroindoline-b (Pinb), and grain softness protein-1 (Gsp-1) genes located at this locus. Puroindolines represent the molecular-genetic basis of grain texture. In addition, the puroindoline gene products have antimicrobial properties with potential role in plant defense. We document the complete coding sequences from the Triticum/Aegilops taxa, rye and barley including the A, D, C, H, M, N, R, S, and U genomes of the Triticeae. Maximum likelihood analyses performed on Bayesian phylogenetic trees showed distinct evolutionary patterns among Pina, Pinb, and Gsp-1. Positive diversifying selection appeared to drive the evolution of at least one of the three genes examined, suggesting that adaptive forces have operated at this locus. Results evidenced positive selection (omega > 4) at Pina and detected amino acid residues along the mature PIN-a protein with a high probability (>95%) of having evolved under adaptation. We hypothesized that positive selection at the Pina region is congruent with its role as a plant defense gene.

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.

Significant expansion of Vicia pannonica genome size mediated by amplification of a single type of giant retroelement

Amplification and eventual elimination of dispersed repeats, especially those of the retroelement origin, account for most of the profound size variability observed among plant genomes. In most higher plants investigated so far, differential accumulation of various families of elements contributes to these differences. Here we report the identification of giant Ty3/gypsy-like retrotransposons from the legume plant Vicia pannonica, which alone make up approximately 38% of the genome of this species. These retrotransposons have structural features of the Ogre elements previously identified in the genomes of pea and Medicago. These features include extreme size (25 kb), the presence of an extra ORF upstream of the gag-pol region, and a putative intron dividing the prot and rt coding sequences. The Ogre elements are evenly dispersed on V. pannonica chromosomes except for terminal regions containing satellite repeats, their individual copies show extraordinary sequence similarity, and at least part of them are transcriptionally active, which suggests their recent amplification. Similar elements were also detected in several other Vicia species but in most cases in significantly lower numbers. However, there was no obvious correlation of the abundance of Ogre sequences with the genome size of these species.

A hAT superfamily transposase recruited by the cereal grass genome

Transposable elements are ubiquitous genomic parasites with an ancient history of coexistence with their hosts. A few cases have emerged recently where these genetic elements have been recruited for normal function in the host organism. We have identified an expressed hobo/Ac/Tam (hAT) family transposase-like gene in cereal grasses which appears to represent such a case. This gene, which we have called gary, is found in one or two copies in barley, two diverged copies in rice and two very similar copies in hexaploid wheat. No gary homologues are found in Arabidopsis. In all three cereal species, an apparently complete 2.5 kb transposase-like open reading frame is present and nucleotide substitution data show evidence for positive selection, yet the predicted gary protein is probably not an active transposase, as judged by the absence of key amino acids required for transposase function. Gary is expressed in wheat and barley spikes and gary cDNA sequences are also found in rice, oat, rye, maize, sorghum and sugarcane. The short inverted terminal repeats, flanked by an eight-nucleotide host sequence duplication, which are characteristic of a hAT transposon are absent. Genetic mapping in barley shows that gary is located on the distal end of the long arm of chromosome 2H. Wheat homologues of gary map to the same approximate location on the wheat group 2 chromosomes by physical bin-mapping and the more closely related of the two rice garys maps to the syntenic location near the bottom of rice chromosome 4. These data suggest that gary has resided in a single genomic location for at least 60 Myr and has lost the ability to transpose, yet expresses a transposase-related protein that is being conserved under host selection. We propose that the gary transposase-like gene has been recruited by the cereal grasses for an unknown function.

Advent of a New Retrotransposon Structure: The Long Form of the Veju Elements

Transposable elements are the main component of plant genomes, especially in grass species. In a previous analysis, we have identified two unusual types of Class I elements, two homologous Veju TRIM elements, but with an unusual long structure. They are formed by the junction of a yet unidentified segment labelled unknown DNA, flanked by the borders of the classical Veju element. Here, we show that the long (Veju_L) and the short forms (Veju_S) coexist within wheat genomes. The associated unknown DNA had always the same origin, and the Veju_L came probably from either illegitimate recombinations or 'template switching' between the Veju_S and another unique unknown DNA sequence. This junction then evolved differently within wheat genomes.

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

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

Insertional polymorphism and antiquity of PDR1 retrotransposon insertions in pisum species

Sequences flanking 73 insertions of the retrotransposon PDR1 have been characterized, together with an additional 270 flanking regions from one side alone, from a diverse collection of Pisum germ plasm. Most of the identified flanking sequences are repetitious DNAs but more than expected (7%) lie within nuclear gene protein-coding regions. The approximate age of 52 of the PDR1 insertions has been determined by measuring sequence divergence among LTR pairs. These data show that PDR1 transpositions occurred within the last 5 MY, with a peak at 1-2.5 MYA. The insertional polymorphism of 68 insertions has been assessed across 47 selected Pisum accessions, representing the diversity of the genus. None of the insertions are fixed, showing that PDR1 insertions can persist in a polymorphic state for millions of years in Pisum. The insertional polymorphism data have been compared with the age estimations to ask what rules control the proliferation of PDR1 insertions in Pisum. Relatively recent insertions (< approximately 1.5 MYA) tend to be found in small subsets of the Pisum accessions set, "middle-aged" insertions (between approximately 1.5 and 2.5 MYA) vary greatly in their occurrence, and older insertions (> approximately 2.5 MYA) are mostly found in small subsets of Pisum. Finally, the average age estimate for PDR1 insertions, together with an existing data set for PDR1 retrotransposon SSAP markers, has been used to derive an estimate of the effective population size for Pisum of approximately 7.5 x 10(5).