Construction and analysis of a BAC library in the grass<Emphasis Type="Italic"> Brachypodium sylvaticum</Emphasis>: its use as a tool to bridge the gap between rice and wheat in elucidating gene content

A BAC library of 30,228 clones with an average insert size of 102 kb was constructed in the grass Brachypodium sylvaticum. Brachypodium has a simple genome, similar in size and repetitive DNA content to that of rice, and is more closely related than rice both to the major temperate cereals wheat and barley, and to the forage grasses. The library represents 6.6 genome equivalents, implying a 99.9% probability of recovering any specific sequence. The library was arrayed onto two high-density colony filters, which were screened with heterologous DNA probes from rice chromosome nine and from syntenous regions of wheat, barley, maize and oat. The construction of Brachypodium BAC contigs revealed that synteny between rice, wheat and Brachypodium was largely maintained over several regions of rice chromosome nine. This suggests that Brachypodium will be a useful tool in the elucidation of gene content in agronomically important cereal crops, complementing rice as a grass genome model.     

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

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

Structure and dynamics of retrotransposons at wheat centromeres and pericentromeres

Little is known of the dynamics of centromeric DNA in polyploid plants. We report the sequences of two centromere-associated bacterial artificial chromosome clones from a Triticum boeoticum library. Both autonomous and non-autonomous wheat centromeric retrotransposons (CRWs) were identified, both being closely associated with the centromeres of wheat. Fiber-fluorescence in situ hybridization and chromatin immunoprecipitation analysis showed that wheat centromeric retrotransposons (CRWs) represent a dominant component of the wheat centromere and are associated with centromere function. CRW copy number showed variation among different genomes: the D genome chromosomes contained fewer copies than either the A or B genome chromosomes. The frequency of lengthy continuous CRW arrays was higher than that in either rice or maize. The dynamics of CRWs and other retrotransposons at centromeric and pericentromeric regions during diploid speciation and polyploidization of wheat and its related species are discussed. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Zhao Liu and Wei Yue made an equal contribution to this work.     

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     

Sequencing chromosome 5D of Aegilops tauschii and comparison with its allopolyploid descendant bread wheat (Triticum aestivum)

Flow cytometric sorting of individual chromosomes and chromosome‐based sequencing reduces the complexity of large, repetitive Triticeae genomes. We flow‐sorted chromosome 5D of Aegilops tauschii, the D genome donor of bread wheat and sequenced it by Roche 454 GS FLX platform to approximately 2.2x coverage. Repetitive sequences represent 81.09% of the survey sequences of this chromosome, and Class I retroelements are the prominent type, with a particular abundance of LTR/Gypsy superfamily. Nonrepetitive sequences were assembled to cover 17.76% of the total chromosome regions. Up to 6188 nonrepetitive gene loci were predicted to be encoded by the 5D chromosome. The numbers and chromosomal distribution patterns of tRNA genes suggest abundance in tRNA^Lys and tRNA^Met species, while the nonrepetitive assembly reveals tRNA^Ala species as the most abundant type. A comparative analysis of the genomic sequences of bread wheat and Aegilops chromosome 5D indicates conservation of gene content. Orthologous unique genes, matching Aegilops 5D sequences, numbered 3730 in barley, 5063 in Brachypodium, 4872 in sorghum and 4209 in rice. In this study, we provide a chromosome‐specific view into the structure and organization of the 5D chromosome of Ae. tauschii, the D genome ancestor of bread wheat. This study contributes to our understanding of the chromosome‐level evolution of the wheat genome and presents a valuable resource in wheat genomics due to the recent hybridization of Ae. tauschii genome with its tetraploid ancestor.     

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

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

HpaII library indicates ‘methylation-free islands’ in wheat and barley

Summary A library of wheat genomic DNA HpaII tiny fragments (HTF), sized below 500 bp, has been constructed. Of the clones in the library 80% belong to the single/low-copy category, while 12% of the clones are nuclear repetitive sequences and 8% originate from the chloroplast and mitochondrial DNA. This result shows a substantial enrichment in the single/low-copy sequences of the wheat genome, which contains at least 80% repetitive sequences. Twenty-nine random single/lowcopy clones were analysed further for wheat chromosome location, cross-hybridisation to barley DNA and their association with rare-cutting, C-methylation-sensitive restriction sites. The results show that the HTF clones are associated more frequently than expected with NotI, MluI, NruI and PstI sites in wheat and barley genomic DNA. The 12% repetitive fraction of the clones contain both moderately and highly repetitive sequences, but no tandemly repeated sequences. The level of enrichment for single/low-copy sequences indicates that libraries of this type are a valuable source of probes for RFLP mapping. In addition, the close association of the HTF clones with rare-cutting restriction enzyme sites ensures that HTF clones will have a useful role in the construction of long-range physical maps in wheat.     

A candidate autonomous version of the wheat MITE <Emphasis Type="Italic">Hikkoshi</Emphasis> is present in the rice genome

A miniature inverted-repeat transposable element (MITE), designated as Hikkoshi, was previously identified in the null Wx-A1 allele of Turkish bread wheat lines. This MITE is 165 bp in size and has 12-bp terminal inverted repeats (TIRs) flanked by 8-bp target site duplications (TSDs). Southern and PCR analyses demonstrated the presence of multiple copies of Hikkoshi in the wheat genome. Database searches indicated that Hikkoshi MITEs are also present in barley, rice and maize. A 3.4-kb element that has Hikkoshi-like TIRs flanked by 8-bp TSDs has now been identified in the rice genome. This element shows high similarity to the 5 subterminal region of the wheat Hikkoshi MITE and contains a transposase (TPase) coding region. The TPase has two conserved domains, ZnF_TTF and hATC, and its amino acid sequence shows a high degree of homology to TPases encoded by Tip100 transposable elements belonging to the hAT superfamily. We designated the 3.4-kb element as OsHikkoshi. Several wheat clones deposited in EST databases showed sequence similarity to the TPase ORF of OsHikkoshi. The sequence information from the TPase of OsHikkoshi will thus be useful in isolating the autonomous element of the Hikkoshi system from wheat.     

BAC-FISH in wheat identifies chromosome landmarks consisting of different types of transposable elements

Fluorescence in situ hybridization (FISH) has been widely used in the physical mapping of genes and chromosome landmarks in plants and animals. Bacterial artificial chromosomes (BACs) contain large inserts making them amenable for FISH mapping. We used BAC-FISH to study genome organization and evolution in hexaploid wheat and its relatives. We selected 56 restriction fragment length polymorphism (RFLP) locus-specific BAC clones from libraries of Aegilops tauschii (the D-genome donor of hexaploid wheat) and A-genome diploid Triticum monococcum. Different types of repetitive sequences were identified using BAC-FISH. Two BAC clones gave FISH patterns similar to the repetitive DNA family pSc119; one BAC clone gave a FISH pattern similar to the repetitive DNA family pAs1. In addition, we identified several novel classes of repetitive sequences: one BAC clone hybridized to the centromeric regions of wheat and other cereal species, except rice; one BAC clone hybridized to all subtelomeric chromosome regions in wheat, rye, barley and oat; one BAC clone contained a localized tandem repeat and hybridized to five D-genome chromosome pairs in wheat; and four BAC clones hybridized only to a proximal region in the long arm of chromosome 4A of hexaploid wheat. These repeats are valuable markers for defined chromosome regions and can also be used for chromosome identification. Sequencing results revealed that all these repeats are transposable elements (TEs), indicating the important role of TEs, especially retrotransposons, in genome evolution of wheat.Communicated by P.B. Moens     

Molecular and cytological analyses of large tracks of centromeric DNA reveal the structure and evolutionary dynamics of maize centromeres

We sequenced two maize bacterial artificial chromosome (BAC) clones anchored by the centromere-specific satellite repeat CentC. The two BACs, consisting of approximately 200 kb of cytologically defined centromeric DNA, are composed exclusively of satellite sequences and retrotransposons that can be classified as centromere specific or noncentromere specific on the basis of their distribution in the maize genome. Sequence analysis suggests that the original maize sequences were composed of CentC arrays that were expanded by retrotransposon invasions. Seven centromere-specific retrotransposons of maize (CRM) were found in BAC 16H10. The CRM elements inserted randomly into either CentC monomers or other retrotransposons. Sequence comparisons of the long terminal repeats (LTRs) of individual CRM elements indicated that these elements transposed within the last 1.22 million years. We observed that all of the previously reported centromere-specific retrotransposons in rice and barley, which belong to the same family as the CRM elements, also recently transposed with the oldest element having transposed approximately 3.8 million years ago. Highly conserved sequence motifs were found in the LTRs of the centromere-specific retrotransposons in the grass species, suggesting that the LTRs may be important for the centromere specificity of this retrotransposon family.     

TEnest: automated chronological annotation and visualization of nested plant transposable elements

Organisms with a high density of transposable elements (TEs) exhibit nesting, with subsequent repeats found inside previously inserted elements. Nesting splits the sequence structure of TEs and makes annotation of repetitive areas challenging. We present TEnest, a repeat identification and display tool made specifically for highly repetitive genomes. TEnest identifies repetitive sequences and reconstructs separated sections to provide full-length repeats and, for long-terminal repeat (LTR) retrotransposons, calculates age since insertion based on LTR divergence. TEnest provides a chronological insertion display to give an accurate visual representation of TE integration history showing timeline, location, and families of each TE identified, thus creating a framework from which evolutionary comparisons can be made among various regions of the genome. A database of repeats has been developed for maize (Zea mays), rice (Oryza sativa), wheat (Triticum aestivum), and barley (Hordeum vulgare) to illustrate the potential of TEnest software. All currently finished maize bacterial artificial chromosomes totaling 29.3 Mb were analyzed with TEnest to provide a characterization of the repeat insertions. Sixty-seven percent of the maize genome was found to be made up of TEs; of these, 95% are LTR retrotransposons. The rate of solo LTR formation is shown to be dissimilar across retrotransposon families. Phylogenetic analysis of TE families reveals specific events of extreme TE proliferation, which may explain the high quantities of certain TE families found throughout the maize genome. The TEnest software package is available for use on PlantGDB under the tools section (http://www.plantgdb.org/prj/TE_nest/TE_nest.html); the source code is available from (http://wiselab.org).     

Characterization of a deep-coverage carrot (<Emphasis Type="Italic">Daucus carota</Emphasis> L.) BAC library and initial analysis of BAC-end sequences

Carrot is the most economically important member of the Apiaceae family and a major source of provitamin A carotenoids in the human diet. However, carrot molecular resources are relatively underdeveloped, hampering a number of genetic studies. Here, we report on the synthesis and characterization of a bacterial artificial chromosome (BAC) library of carrot. The library is 17.3-fold redundant and consists of 92,160 clones with an average insert size of 121 kb. To provide an overview of the composition and organization of the carrot nuclear genome we generated and analyzed 2,696 BAC-end sequences (BES) from nearly 2,000 BACs, totaling 1.74 Mb of BES. This analysis revealed that 14% of the BES consists of known repetitive elements, with transposable elements representing more than 80% of this fraction. Eleven novel carrot repetitive elements were identified, covering 8.5% of the BES. Analysis of microsatellites showed a comparably low frequency for these elements in the carrot BES. Comparisons of the translated BES with protein databases indicated that approximately 10% of the carrot genome represents coding sequences. Moreover, among eight dicot species used for comparison purposes, carrot BES had highest homology to protein-coding sequences from tomato. This deep-coverage library will aid carrot breeding and genetics. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Nucleotide sequence data reported are available in the DDBJ/EMBL/GenBank databases under the accession numbers FJ147695–FJ150390.     

<Emphasis Type="Italic">Mutator</Emphasis>-like elements identified in melon,Arabidopsis and rice contain ULP1 protease domains

The transposon Mutator was first identified in maize, and is one of the most active mobile elements in plants. The Arabidopsis thaliana genome contains at least 200 Mutator-like elements (MULEs), which contain the Mutator-like transposase gene, and often additional genes. We have detected a novel type of MULEs in melon (CUMULE), which, besides the transposase, contains two ubiquitin-like specific protease-like sequences (ULP1). This element is not present in the observed location in some melon cultivars. Multiple copies of this element exist in the Cucumis melo genome, and it has been detected in other Cucurbitaceae species. Analysis of the A. thaliana genome revealed more than 90 CUMULE-like elements, containing one or two Ulp1-like sequences, although no evidence of mobility exists for these elements. We detected various putative transposable elements containing ULP1-like sequences in rice. The discovery of these MULEs in melon and Arabidopsis, and the existence of similar elements in rice and maize, suggest that a proteolytic function may be important for this subset of the MULE transposable elements. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users. Nucleotide sequence data reported are available in the GenBank database under the accession number AY524004.     

The occurrence of Ds-like sequences in cereal genomes

Summary The occurrence of DNA sequences similar to the Ds-element of sh-m5933 maize (Ds-like sequences) was studied in other representatives of the Gramineae. The approximate number of copies of such sequences found under gentle and stringent conditions of washing was determined by dot-hybridization. It was shown that in the maize genome the number of copies of Ds-like sequences exceeds about ten-fold the content of such sequences found in wheat, rye and barley genomes. Quantitative differences in Ds-like sequences between wheat species with various genomes and ploidies (when estimated per genome) as well as between different H. vulgare varieties was not determined. The various melting points (T_m) of DNA-duplexes formed when the Ds-element is hybridized with wheat, rye and barley DNA respectively do not show significant differences and are essentially lower than the T_m of the Ds-element (by 8°–9°C). Thus, these duplexes have 9–11% of nucleotide substitutions in comparison to Ds sh-m5933. The data obtained permit one to suppose the presence of a series of Ds-like sequences heterogenous for the length and degree of homology to the Ds-element isolated from the shrunken locus (sh-m5933) of maize DNA.     

A structural and evolutionary analysis of a dispersed repetitive sequence

A family of dispersed repetitive sequences (Hch1) which is present in the genome of the wild barley Hordeum chilense was studied in detail. Hch1 sequences are found both as part of short tandem arrays and dispersed throughout the H. chilense chromosomes. Subcloning of sections of the sequence reveals that it is composed of unrelated classes of sequences which can also be found separately in other genomic locations. Analysis of these sequences in the genomes of wheat and two other wild barley species strongly suggests that specific amplifications and arrangements of the repeated sequences have taken place during speciation. Nucleotide sequence analysis fails to detect, in their entirity, the features shown by plant transposons.     

Genome size and sequence composition of moso bamboo: A comparative study

Moso bamboo (Phyllostachys pubescens) is one of the world’s most important bamboo species. It has the largest area of all planted bamboo—over two-thirds of the total bamboo forest area—and the highest economic value in China. Moso bamboo is a tetraploid (4x=48) and a special member of the grasses family. Although several genomes have been sequenced or are being sequenced in the grasses family, we know little about the genome of the bambusoids (bamboos). In this study, the moso bamboo genome size was estimated to be about 2034 Mb by flow cytometry (FCM), using maize (cv. B73) and rice (cv. Nipponbare) as internal references. The rice genome has been sequenced and the maize genome is being sequenced. We found that the size of the moso bamboo genome was similar to that of maize but significantly larger than that of rice. To determine whether the bamboo genome had a high proportion of repeat elements, similar to that of the maize genome, approximately 1000 genome survey sequences (GSS) were generated. Sequence analysis showed that the proportion of repeat elements was 23.3% for the bamboo genome, which is significantly lower than that of the maize genome (65.7%). The bamboo repeat elements were mainly Gypsy/DIRS1 and Ty1/Copia LTR retrotransposons (14.7%), with a few DNA transposons. However, more genomic sequences are needed to confirm the above results due to several factors, such as the limitation of our GSS data. This study is the first to investigate sequence composition of the bamboo genome. Our results are valuable for future genome research of moso and other bamboos.     

Analysis of a contiguous 211 kb sequence in diploid wheat (Triticum monococcum L.) reveals multiple mechanisms of genome evolution

In plant species with large genomes such as wheat or barley, genome organization at the level of DNA sequence is largely unknown. The largest sequences that are publicly accessible so far from Triticeae genomes are two 60 kb and 66 kb intervals from barley. Here, we report on the analysis of a 211 kb contiguous DNA sequence from diploid wheat (Triticum monococcum L.). Five putative genes were identified, two of which show similarity to disease resistance genes. Three of the five genes are clustered in a 31 kb gene-enriched island while the two others are separated from the cluster and from each other by large stretches of repetitive DNA. About 70% of the contig is comprised of several classes of transposable elements. Ten different types of retrotransposons were identified, most of them forming a pattern of nested insertions similar to those found in maize and barley. Evidence was found for major deletion, insertion and duplication events within the analysed region, suggesting multiple mechanisms of genome evolution in addition to retrotransposon amplification. Seven types of foldback transposons, an element class previously not described for wheat genomes, were characterized. One such element was found to be closely associated with genes in several Triticeae species and may therefore be of use for the identification of gene-rich regions in these species.     

A genome-wide survey of switchgrass genome structure and organization

The perennial grass, switchgrass (Panicum virgatum L.), is a promising bioenergy crop and the target of whole genome sequencing. We constructed two bacterial artificial chromosome (BAC) libraries from the AP13 clone of switchgrass to gain insight into the genome structure and organization, initiate functional and comparative genomic studies, and assist with genome assembly. Together representing 16 haploid genome equivalents of switchgrass, each library comprises 101,376 clones with average insert sizes of 144 (HindIII-generated) and 110 kb (BstYI-generated). A total of 330,297 high quality BAC-end sequences (BES) were generated, accounting for 263.2 Mbp (16.4%) of the switchgrass genome. Analysis of the BES identified 279,099 known repetitive elements, >50,000 SSRs, and 2,528 novel repeat elements, named switchgrass repetitive elements (SREs). Comparative mapping of 47 full-length BAC sequences and 330K BES revealed high levels of synteny with the grass genomes sorghum, rice, maize, and Brachypodium. Our data indicate that the sorghum genome has retained larger microsyntenous regions with switchgrass besides high gene order conservation with rice. The resources generated in this effort will be useful for a broad range of applications.     

Comparative analysis in cereals of a key proline catabolism gene.

Proline accumulation and catabolism play significant roles in adaptation to a variety of plant stresses including osmotic stress, drought, temperature, freezing, UV irradiation, heavy metals and pathogen infection. In this study, the gene Δ¹ -pyrroline-5-carboxylate dehydrogenase (P5CDH), which catalyzes the second step in the conversion of proline to glutamate, is characterized in a number of cereal species. P5CDH genes from hexaploid wheat, Triticum turgidum (durum wheat), Aegilops tauschii, Triticum monococcum, barley, maize and rice were shown to be conserved in terms of gene structure and sequence, present as a single copy per haploid, non-polyploid genome and located in evolutionarily conserved linkage groups. A wheat cDNA sequence was shown by yeast complementation to encode a functional P5CDH activity. A divergently-transcribed rab7 gene was identified immediately 5′ of P5CDH in all grasses examined, except rice. The rab7/P5CDH intergenic region in these species, which presumably encompasses 5′ regulatory elements of both genes, showed a distinct pattern of sequence evolution with sequences in juxtaposition to each ORF conserved between barley, wheat, A. tauschii and T. monococcum. More distal 5′ sequence in this intergenic region showed a higher rate of divergence, with no homology observed between these regions in the wheat and barley genomes. Maize and rice showed no similarity in regions 5′ of P5CDH when compared with wheat, barley, and each other, apart from a 22 bp region of conserved non-coding sequence (CNS) that is similar to a proline response element identified in the promoter of the Arabidopsis proline dehydrogenase gene. A palindromic motif similar to this cereal CNS was also identified 5′ of the Arabidopsis AtP5CDH gene showing conservation of this sequence in monocot and dicot lineages.