DNA rearrangement in orthologous orp regions of the maize, rice and sorghum genomes

The homeologous Orp1 and Orp2 regions of maize and the orthologous regions in sorghum and rice were compared by generating sequence data for >486 kb of genomic DNA. At least three genic rearrangements differentiate the maize Orp1 and Orp2 segments, including an insertion of a single gene and two deletions that removed one gene each, while no genic rearrangements were detected in the maize Orp2 region relative to sorghum. Extended comparison of the orthologous Orp regions of sorghum and japonica rice uncovered numerous genic rearrangements and the presence of a transposon-rich region in rice. Only 11 of 27 genes (40%) are arranged in the same order and orientation between sorghum and rice. Of the 8 genes that are uniquely present in the sorghum region, 4 were found to have single-copy homologs in both rice and Arabidopsis, but none of these genes are located near each other, indicating frequent gene movement. Further comparison of the Orp segments from two rice subspecies, japonica and indica, revealed that the transposon-rich region is both an ancient and current hotspot for retrotransposon accumulation and genic rearrangement. We also identify unequal gene conversion as a mechanism for maize retrotransposon rearrangement.     

A highly conserved, small LTR retrotransposon that preferentially targets genes in grass genomes

LTR retrotransposons are often the most abundant components of plant genomes and can impact gene and genome evolution. Most reported LTR retrotransposons are large elements (>4 kb) and are most often found in heterochromatic (gene poor) regions. We report the smallest LTR retrotransposon found to date, only 292 bp. The element is found in rice, maize, sorghum and other grass genomes, which indicates that it was present in the ancestor of grass species, at least 50-80 MYA. Estimated insertion times, comparisons between sequenced rice lines, and mRNA data indicate that this element may still be active in some genomes. Unlike other LTR retrotransposons, the small LTR retrotransposons (SMARTs) are distributed throughout the genomes and are often located within or near genes with insertion patterns similar to MITEs (miniature inverted repeat transposable elements). Our data suggests that insertions of SMARTs into or near genes can, in a few instances, alter both gene structures and gene expression. Further evidence for a role in regulating gene expression, SMART-specific small RNAs (sRNAs) were identified that may be involved in gene regulation. Thus, SMARTs may have played an important role in genome evolution and genic innovation and may provide a valuable tool for gene tagging systems in grass.     

Retrotransposon insertion polymorphisms in six rice genes and their evolutionary history

Retrotransposons are abundant in higher plant genomes. Although retrotransposons associated with plant genes have been identified, little is known about their evolutionary conservation at the level of species and subspecies. In the present study, we investigated the phylogenetic distribution of long terminal repeat (LTR) retrotransposon, long interspersed nuclear element (LINE) and short interspersed nuclear element (SINE) insertions in six genes in 95 cultivated and wild rice genotypes. These six genes are likely to be functional based on nonsynonymous (Ka) to synonymous (Ks) substitution ratios which were found to be significantly <1. Different conservation patterns of these retrotransposons in genes were observed in cultivated and wild rice species. Four out of seven retrotransposon insertions appear to predate the ancestral Oryza AA genome. Two of these insertions in genes 4 and 5 occurred early in the evolutionary history of Oryza. Two retrotransposon insertions in gene 1 arose after the divergence of Asian cultivated rice from its wild ancestor. Furthermore, the retrotransposon insertion in gene 3 appears to have occurred in the ancestral lineage leading to temperate japonicas. Conservation of retrotransposon insertions in genes in specific groups, species, and lineages might be related to their specific function.     

The colinearity of the Sh2/A1 orthologous region in rice, sorghum and maize is interrupted and accompanied by genome expansion in the triticeae

The Sh2/A1 orthologous region of maize, rice, and sorghum contains five genes in the order Sh2, X1, X2, and two A1 homologs in tandem duplication. The Sh2 and A1 homologs are separated by approximately 20 kb in rice and sorghum and by approximately 140 kb in maize. We analyzed the fate of the Sh2/A1 region in large-genome species of the Triticeae (wheat, barley, and rye). In the Triticeae, synteny in the Sh2/A1 region was interrupted by a break between the X1 and X2 genes. The A1 and X2 genes remained colinear in homeologous chromosomes as in other grasses. The Sh2 and X1 orthologs also remained colinear but were translocated to a nonhomeologous chromosome. Gene X1 was duplicated on two nonhomeologous chromosomes, and surprisingly, a paralog shared homology much higher than that of the orthologous copy to the X1 gene of other grasses. No tandem duplication of A1 homologs was detected but duplication of A1 on a nonhomeologous barley chromosome 6H was observed. Intergenic distances expanded greatly in wheat compared to rice. Wheat and barley diverged from each other 12 million years ago and both show similar changes in the Sh2/A1 region, suggesting that the break in colinearity as well as X1 duplications and genome expansion occurred in a common ancestor of the Triticeae species.     

The genetic colinearity of rice and other cereals on the basis of genomic sequence analysis

Small segments of rice genome sequence have been compared with that of the model plant Arabidopsis thaliana and with several closer relatives, including the cereals maize, rice, sorghum, barley and wheat. The rice genome is relatively stable relative to those of other grasses. Nevertheless, comparisons with other cereals have demonstrated that the DNA between cereal genes is highly variable and evolves rapidly. Genic regions have undergone many more small rearrangements than have been revealed by recombinational mapping studies. Tandem gene duplication/deletion is particularly common, but other types of deletions, inversions and translocations also occur. The many thousands of small genic rearrangements within the rice genome complicate but do not negate its use as a model for larger cereal genomes.     

Transposable elements,genes and recombination in a 215-kb contig from wheat chromosome 5A(m)

Sequencing of a contiguous 215-kb interval of Triticum monococcum showed the presence of five genes in the same order as in previously sequenced colinear barley and rice BACs. Gene 2 was in the same orientation in wheat and rice but inverted in barley. Gene density in this region was 1 gene per 43 kb and the ratio of physical to genetic distance was estimated to be 2,700 kb cM^–1. Twenty more-or-less intact retrotransposons were found in the intergenic regions, covering at least 70% of the sequenced region. The insertion times of 11 retrotransposons were less than 5 million years ago and were consistent with their nested structure. Five new families of retroelements and the first full-length elements for two additional retrotransposon families were discovered in this region. Significantly higher values of GC content were observed for Triticeae BACs compared with rice BACs. Relative enrichment or depletion of certain dinucleotides was observed in the comparison of introns, exons and retrotransposons. A higher proportion of transitions in CG and CNG sites that are targets for cytosine methylation was observed in retrotransposons (76%) than in introns (37%). These results showed that the wheat genome is a complex mixture of different sequence elements, but with general patterns of content and interspersion that are similar to those seen in maize and barley. Electronic Publication     

Gene-containing regions of wheat and the other grass genomes

Deletion line-based high-density physical maps revealed that the wheat (Triticum aestivum) genome is partitioned into gene-rich and -poor compartments. Available deletion lines have bracketed the gene-containing regions to about 10% of the genome. Emerging sequence data suggest that these may further be partitioned into "mini" gene-rich and gene-poor regions. An average of about 10% of each gene-rich region seem to contain genes. Sequence analyses in various species suggest that uneven distribution of genes may be a characteristic of all grasses and perhaps all higher organisms. Comparison of the physical maps with genetic linkage maps showed that recombination in wheat and barley (Hordeum vulgare) is confined to the gene-containing regions. Number of genes, gene density, and the extent of recombination vary greatly among the gene-rich regions. The gene order, relative region size, and recombination are highly conserved within the tribe Triticeae and moderately conserved within the family. Gene-poor regions are composed of retrotransposon-like non-transcribing repeats and pseudogenes. Direct comparisons of orthologous regions indicated that gene density in wheat is about one-half compared with rice (Oryza sativa). Genome size difference between wheat and rice is, therefore, mainly because of amplification of the gene-poor regions. Presence of species-, genera-, and family-specific repeats reveal a repeated invasion of the genomes by different retrotransposons over time. Preferential transposition to adjacent locations and presence of vital genes flanking a gene-rich region may have restricted retrotransposon amplification to gene-poor regions, resulting into tandem blocks of non-transcribing repeats. Insertional inactivation by adjoining retro-elements and selection seem to have played a major role in stabilizing genomes.     

Analysis of genes associated with retrotransposons in the rice genome

Retrotransposons comprise a significant fraction of the rice genome. Despite their prevalence, the effects of retrotransposon insertions are not well understood, especially with regard to how they affect the expression of genes. In this study, we identified one-sixth of rice genes as being associated with retrotransposons, with insertions either in the gene itself or within its putative promoter region. Among genes with insertions in the promoter region, the likelihood of the gene being expressed was shown to be directly proportional to the distance of the retrotransposon from the translation start site. In addition, retrotransposon insertions in the transcribed region of the gene were found to be positively correlated with the presence of alternative splicing forms. Furthermore, preferential association of retrotransposon insertions with genes in several functional classes was identified. Some of the retrotransposons that are part of full-length cDNA (fl-cDNA) contribute splice sites and give rise to novel exons. Several interesting trends concerning the effects of retrotransposon insertions on gene expression were identified. Taken together, our data suggests that retrotransposon association with genes have a role in gene regulation. The data presented in this study provides a foundation for experimental studies to determine the role of retrotransposons in gene regulation.     

Structure–function analysis of the barley genome: the gene-rich region of chromosome 2HL

A major gene-rich region on the end of the long arm of Triticeae group 2 chromosomes exhibits high recombination frequencies, making it an attractive region for positional cloning. Traits known to be controlled by this region include chasmogamy/cleistogamy, frost tolerance at flowering, grain yield, head architecture, and resistance to Fusarium head blight and rusts. To assist these cloning efforts, we constructed detailed genetic maps of barley chromosome 2H, including 61 polymerase chain reaction markers. Colinearity with rice occurred in eight distinct blocks, including five blocks in the terminal gene-rich region. Alignment of rice sequences from the junctions of colinear chromosome segments provided no evidence for the involvement of long (>2.5 kb) inverted repeats in generating inversions. However, reuse of some junction sequences in two or three separate evolutionary breakage/fusion events was implicated, suggesting the presence of fragile sites. Sequencing across 91 gene fragments totaling 107 kb from four barley genotypes revealed the highest single nucleotide substitution and insertion–deletion polymorphism levels in the terminal regions of the chromosome arms. The maps will assist in the isolation of genes from the chromosome 2L gene-rich region in barley and wheat by providing markers and accelerating the identification of the corresponding points in the rice genome sequence. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Comparative mapping of the barley Ppd-H1 photoperiod response gene region,which lies close to a junction between two rice linkage segments

Comparative mapping of cereals has shown that chromosomes of barley, wheat, and maize can be described in terms of rice "linkage segments." However, little is known about marker order in the junctions between linkage blocks or whether this will impair comparative analysis of major genes that lie in such regions. We used genetic and physical mapping to investigate the relationship between the distal part of rice chromosome 7L, which contains the Hd2 heading date gene, and the region of barley chromosome 2HS containing the Ppd-H1 photoperiod response gene, which lies near the junction between rice 7 and rice 4 linkage segments. RFLP markers were mapped in maize to identify regions that might contain Hd2 or Ppd-H1 orthologs. Rice provided useful markers for the Ppd-H1 region but comparative mapping was complicated by loss of colinearity and sequence duplications that predated the divergence of rice, maize, and barley. The sequences of cDNA markers were used to search for homologs in the Arabidopsis genome. Homologous sequences were found for 13 out of 16 markers but they were dispersed in Arabidopsis and did not identify any candidate equivalent region. The implications of the results for comparative trait mapping in junction regions are discussed.     

Whole chloroplast genome comparison of rice,maize, and wheat: implications for chloroplast gene diversification and phylogeny of cereals

The fully sequenced chloroplast genomes of maize (subfamily Panicoideae), rice (subfamily Bambusoideae), and wheat (subfamily Pooideae) provide the unique opportunity to investigate the evolution of chloroplast genes and genomes in the grass family (Poaceae) by whole-genome comparison. Analyses of nucleotide sequence variations in 106 cereal chloroplast genes with tobacco sequences as the outgroup suggested that (1) most of the genic regions of the chloroplast genomes of maize, rice, and wheat have evolved at similar rates; (2) RNA genes have highly conservative evolutionary rates relative to the other genes; (3) photosynthetic genes have been under strong purifying selection; (4) between the three cereals, 14 genes which account for about 28% of the genic region have evolved with heterogeneous nucleotide substitution rates; and (5) rice genes tend to have evolved more slowly than the others at loci where rate heterogeneity exists. Although the mechanism that underlies chloroplast gene diversification is complex, our analyses identified variation in nonsynonymous substitution rates as a genetic force that generates heterogeneity, which is evidence of selection in chloroplast gene diversification at the intrafamilial level. Phylogenetic trees constructed with the variable nucleotide sites of the chloroplast genes place maize basal to the rice-wheat clade, revealing a close relationship between the Bambusoideae and Pooideae.     

Frequent gene movement and pseudogene evolution is common to the large and complex genomes of wheat, barley, and their relatives

All six arms of the group 1 chromosomes of hexaploid wheat (Triticum aestivum) were sequenced with Roche/454 to 1.3- to 2.2-fold coverage and compared with similar data sets from the homoeologous chromosome 1H of barley (Hordeum vulgare). Six to ten thousand gene sequences were sampled per chromosome. These were classified into genes that have their closest homologs in the Triticeae group 1 syntenic region in Brachypodium, rice (Oryza sativa), and/or sorghum (Sorghum bicolor) and genes that have their homologs elsewhere in these model grass genomes. Although the number of syntenic genes was similar between the homologous groups, the amount of nonsyntenic genes was found to be extremely diverse between wheat and barley and even between wheat subgenomes. Besides a small core group of genes that are nonsyntenic in other grasses but conserved among Triticeae, we found thousands of genic sequences that are specific to chromosomes of one single species or subgenome. By examining in detail 50 genes from chromosome 1H for which BAC sequences were available, we found that many represent pseudogenes that resulted from transposable element activity and double-strand break repair. Thus, Triticeae seem to accumulate nonsyntenic genes frequently. Since many of them are likely to be pseudogenes, total gene numbers in Triticeae are prone to pronounced overestimates.     

Comparative sequence analysis of colinear barley and rice bacterial artificial chromosomes

Colinearity of a large region from barley (Hordeum vulgare) chromosome 5H and rice (Oryza sativa) chromosome 3 has been demonstrated by mapping of several common restriction fragment-length polymorphism clones on both regions. One of these clones, WG644, was hybridized to rice and barley bacterial artificial chromosome (BAC) libraries to select homologous clones. One BAC from each species with the largest overlapping segment was selected by fingerprinting and blot hybridization with three additional restriction fragment-length polymorphism clones. The complete barley BAC 635P2 and a 50-kb segment of the rice BAC 36I5 were completely sequenced. A comparison of the rice and barley DNA sequences revealed the presence of four conserved regions, containing four predicted genes. The four genes are in the same orientation in rice, but the second gene is in inverted orientation in barley. The fourth gene is duplicated in tandem in barley but not in rice. Comparison of the homeologous barley and rice sequences assisted the gene identification process and helped determine individual gene structures. General gene structure (exon number, size, and location) was largely conserved between rice and barley and to a lesser extent with homologous genes in Arabidopsis. Colinearity of these four genes is not conserved in Arabidopsis compared with the two grass species. Extensive similarity was not found between the rice and barley sequences other than within the exons of the structural genes, and short stretches of homology in the promoters and 3' untranslated regions. The larger distances between the first three genes in barley compared with rice are explained by the insertion of different transposable retroelements.     

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

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

Conserved noncoding sequences among cultivated cereal genomes identify candidate regulatory sequence elements and patterns of promoter evolution

Surveys for conserved noncoding sequences (CNS) among genes from monocot cereal species were conducted to assess the general properties of CNS in grass genomes and their correlation with known promoter regulatory elements. Initial comparisons of 11 orthologous maize-rice gene pairs found that previously defined regulatory motifs could be identified within short CNS but could not be distinguished reliably from random sequence matches. Among the different phylogenetic footprinting algorithms tested, the VISTA tool yielded the most informative alignments of noncoding sequence. VISTA was used to survey for CNS among all publicly available genomic sequences from maize, rice, wheat, barley, and sorghum, representing >300 gene comparisons. Comparisons of orthologous maize-rice and maize-sorghum gene pairs identified 20 bp as a minimal length criterion for a significant CNS among grass genes, with few such CNS found to be conserved across rice, maize, sorghum, and barley. The frequency and length of cereal CNS as well as nucleotide substitution rates within CNS were consistent with the known phylogenetic distances among the species compared. The implications of these findings for the evolution of cereal gene promoter sequences and the utility of using the nearly completed rice genome sequence to predict candidate regulatory elements in other cereal genes by phylogenetic footprinting are discussed.     

Polymorphisms and evolutionary history of retrotransposon insertions in rice promoters

Retrotransposons are ubiquitous in higher plant genomes. The presence or absence of retrotransposons in whole genome and high throughput genomic sequence (HTGS) from cultivated and wild rice was investigated to understand the organization and evolution of retrotransposon insertions in promoter regions. Approximately half of the Oryza sativa subsp. japonica 'Nipponbare' promoters with retrotransposons conserved in Oryza sativa subsp. indica '93-11' and four wild rice species showed higher sequence conservation in retrotransposon than nonretrotransposon regions. We further investigated, in detail, the evolutionary dynamics of five retrotransposons in the promoter regions of 95 rice genotypes. Our data suggest that four of five insertions (Rp2-Rp5) occurred in the ancestor of AA genome, while the other insertion (Rp1) predates the ancestral divergence of Oryza officinalis (CC genome). Four retrotransposons (Rp2-Rp5) were present in 52% (Rp2), 29% (Rp3), 53% (Rp4), and 43% (Rp5) of the rice genotypes with AA genome type, and the fifth retrotransposon (Rp1) was present in 95% of the rice genotypes with AA, BBCC, or CC genome types. Furthermore, most of these retrotransposons were found to evolve slower than flanking promoter regions, suggesting a role in promoter function for regulating downstream genes.     

Genes controlling seed dormancy and pre-harvest sprouting in a rice-wheat-barley comparison

Pre-harvest sprouting results in significant economic loss for the grain industry around the world. Lack of adequate seed dormancy is the major reason for pre-harvest sprouting in the field under wet weather conditions. Although this trait is governed by multiple genes it is also highly heritable. A major QTL controlling both pre-harvest sprouting and seed dormancy has been identified on the long arm of barley chromosome 5H, and it explains over 70% of the phenotypic variation. Comparative genomics approaches among barley, wheat and rice were used to identify candidate gene(s) controlling seed dormancy and hence one aspect of pre-harvest sprouting. The barley seed dormancy/pre-harvest sprouting QTL was located in a region that showed good synteny with the terminal end of the long arm of rice chromosome 3. The rice DNA sequences were annotated and a gene encoding GA20-oxidase was identified as a candidate gene controlling the seed dormancy/pre-harvest sprouting QTL on 5HL. This chromosomal region also shared synteny with the telomere region of wheat chromosome 4AL, but was located outside of the QTL reported for seed dormancy in wheat. The wheat chromosome 4AL QTL region for seed dormancy was syntenic to both rice chromosome 3 and 11. In both cases, corresponding QTLs for seed dormancy have been mapped in rice.C. Li and P. Ni contributed equally to this work     

The wheat D-genome HMW-glutenin locus: BAC sequencing,gene distribution,and retrotransposon clusters

A bacterial-artificial-chromosome (BAC) clone from the genome of Triticum tauschii, the D-genome ancestor of hexaploid bread wheat, was sequenced and the presence of the two paralogous x- and y-type high-molecular-weight (HMW) glutenin genes of the Glu-D1 locus was confirmed. These two genes occur in the same orientation, are 51,893 bp apart, and the separating DNA includes a 31,000-bp cluster of retrotransposons. A second retrotransposon cluster of 32,000 bp follows the x-type HMW-glutenin gene region. Each HMW-glutenin gene is found within a region of mainly unique DNA sequence which includes multiple additional genes including an active endosperm globulin gene not previously reported in the Triticeae family, a leucine-rich-repeat (LRR) type gene truncated at the 5′ end of the BAC, a kinase gene of unknown activity, remnants of a paralogous second globulin gene, and genes similar to two hypothetical rice genes. The newly identified globulin genes are assigned to a locus designated Glo-2. Comparison to available orthologous regions of the wheat A and B genomes show rapid sequence divergences flanking the HMW-glutenin genes, and the absence of two hypothetical and unknown genes found 5′ to the B-genome x-type ortholog. The region surrounding the Glu-D1 locus is similar to other reported Triticeae BAC sequences; i.e. small gene islands separated by retrotransposon clusters. Electronic Publication