Skip to main content

Chromosome arm-specific BAC end sequences permit comparative analysis of homoeologous chromosomes and genomes of polyploid wheat

Abstract

Background

Bread wheat, one of the world’s staple food crops, has the largest, highly repetitive and polyploid genome among the cereal crops. The wheat genome holds the key to crop genetic improvement against challenges such as climate change, environmental degradation, and water scarcity. To unravel the complex wheat genome, the International Wheat Genome Sequencing Consortium (IWGSC) is pursuing a chromosome- and chromosome arm-based approach to physical mapping and sequencing. Here we report on the use of a BAC library made from flow-sorted telosomic chromosome 3A short arm (t3AS) for marker development and analysis of sequence composition and comparative evolution of homoeologous genomes of hexaploid wheat.

Results

The end-sequencing of 9,984 random BACs from a chromosome arm 3AS-specific library (TaaCsp3AShA) generated 11,014,359 bp of high quality sequence from 17,591 BAC-ends with an average length of 626 bp. The sequence represents 3.2% of t3AS with an average DNA sequence read every 19 kb. Overall, 79% of the sequence consisted of repetitive elements, 1.38% as coding regions (estimated 2,850 genes) and another 19% of unknown origin. Comparative sequence analysis suggested that 70-77% of the genes present in both 3A and 3B were syntenic with model species. Among the transposable elements, gypsy/sabrina (12.4%) was the most abundant repeat and was significantly more frequent in 3A compared to homoeologous chromosome 3B. Twenty novel repetitive sequences were also identified using de novo repeat identification. BESs were screened to identify simple sequence repeats (SSR) and transposable element junctions. A total of 1,057 SSRs were identified with a density of one per 10.4 kb, and 7,928 junctions between transposable elements (TE) and other sequences were identified with a density of one per 1.39 kb. With the objective of enhancing the marker density of chromosome 3AS, oligonucleotide primers were successfully designed from 758 SSRs and 695 Insertion Site Based Polymorphisms (ISBPs). Of the 96 ISBP primer pairs tested, 28 (29%) were 3A-specific and compared to 17 (18%) for 96 SSRs.

Conclusion

This work reports on the use of wheat chromosome arm 3AS-specific BAC library for the targeted generation of sequence data from a particular region of the huge genome of wheat. A large quantity of sequences were generated from the A genome of hexaploid wheat for comparative genome analysis with homoeologous B and D genomes and other model grass genomes. Hundreds of molecular markers were developed from the 3AS arm-specific sequences; these and other sequences will be useful in gene discovery and physical mapping.

Background

Plant genomes may vary over 2,000-fold in genome size (Genlisea aurea, ~64 Mb to Paris japonica, ~149 Gb). Within the grass family, bread wheat, also known as common or hexaploid wheat (Triticum aestivum L., 2n = 6x = 42), has a large genome size of 17 Gb because of two allopolypoidization events and also a very high repetitive DNA content (reviewed in Gill et al. [1]). Each chromosome of wheat, at an average size of 809 Mb (17 Gb/21), is twice the size of the rice genome (390 Mb). Therefore, the physical mapping and sequencing of the wheat genome poses a technical challenge and is cost prohibitive, as discussed at an NSF-sponsored workshop [1]. The wheat genomics community decided on a chromosome-based approach to physical mapping and sequencing to overcome these problems and to provide an opportunity for division of labor and cost sharing by establishing an international consortium of scientists (http://www.wheatgenome.org/). This approach is feasible because of the ability to purify wheat chromosomes and/or chromosome arms by flow cytometric sorting. Particular chromosomes and chromosome arms can be isolated at over 90 % purity and in sufficient quantity for making chromosome and chromosome arm-specific BAC libraries [24].

Among the wheat chromosomes, 3B is the largest, sorts out as a single peak on a flow karyotype and its BAC library [5] was used to construct the first wheat chromosome physical map [6]. Another peak contains three of the smallest wheat chromosomes 1D, 4D and 6D, and their physical map construction is underway (http://wggrc.plantpath.ksu.edu/wheat/Dgenome/). For the remaining wheat chromosomes, lines carrying telocentric chromosomes, which represent individual arms, must be used to increase sorting resolution and library construction [2, 3]. Chromosome and chromosome arm-specific BAC libraries now have been constructed for most of the wheat chromosomes (http://lmcc.ieb.cz/dna-libraries/cereals)[4].

Flow-sorted chromosomes and libraries have been utilized for generating chromosome and chromosome arm-specific sequence information for gene discovery and comparative genome analysis. BAC end sequences (BESs) have been analyzed from BAC libraries and can be highly informative in determining genome content and organization, enhance the value of BACs as a genomic resource and also provide random sequence information. BESs have been used to estimate the distribution of the repetitive elements including the retrotransposons, DNA transposons and SSRs in several plant species [714]. Paux et al. [9] generated 11 Mb of random BESs from chromosome 3B of bread wheat and reported 86% of the sequences as repetitive elements, 1.2% as coding region and 13% was unknown. BESs can be an excellent source for marker development for genetic and physical mapping. BESs were used to develop simple sequence repeats, genic sequence-based markers and, more recently, inserted transposable elements junction based ISBP markers [9, 10, 1315].

Sequences from the whole-genome amplification of flow-sorted chromosomes or chromosome-arms have been published for several wheat and barley chromosomes [1622]. Wicker et al. [19] sequenced flow-sorted wheat group-1 chromosomes by Roche/454 technology at 1.3–2.2x coverage and used the sequence data to estimate gene syntenic relationships of the Triticeae with Brachypodium, rice and sorghum and suggested a large number of sequences that are nonsyntenic to model grasses are probably pseudogenes. Berkman et al. [17, 21] recently sequenced the flow sorted chromosomes 7BS and 7DS of wheat to isolate the low copy and genic sequences and observed ~ 60% rate of colinearity with Brachypodium. Recently, Rustenholnz et al. studied the gene space organization in wheat by mapping the expressed portion of chromosome 3B using barley microarrays [23] and a wheat unigene microarray [24] and observed a 2-fold increase in gene density from centromere towards the telomeres. A BES [9] and transcript map [24] of chromosome 3B provide an opportunity for comparative analysis of homoeologous chromosomes.

In this study, we used the first arm-specific library in wheat from chromosome 3AS [25] for which we are developing a physical map. The 3AS BAC library consists of 55, 296 clones (physical map 3ASv1.0; http://wggrc.plantpath.ksu.edu/wheat/3A/3A_index.html). We randomly selected nearly 10,000 BACs and end sequenced them to generate 11 Mb of sequence. These sequences were analyzed to obtain insight into the sequence composition of wheat chromosome arm 3AS and its evolutionary relationship with other homoeologous chromosomes and model grass genomes. Furthermore, we identified DNA motifs potentially useful as molecular markers to saturate chromosome 3AS.

Results and discussion

BAC-end sequencing and data processing

Chromosome 3A of bread wheat is among the largest wheat chromosomes with a metaphase size of 11.8 μm, is submetacentric with an arm ratio of 1.3 and accounts for 14.4% of the total A genome of hexaploid wheat [26]. Chromosome arm 3AS is estimated at ~355 Mb, equivalent to 0.8 times the size of the rice genome. We have fingerprinted and assembled contigs of 47,063 clones of the bread wheat cultivar Chinese Spring chromosome arm 3AS-specific BAC library (physical map 3ASv1.0; http://wggrc.plantpath.ksu.edu/wheat/3A/3A_index.html). The first phase of assembly of the 3AS physical map used FPC 8.5.2 [27], where 35,124 clones are assembled into 1,677 contigs with the remaining 11,939 clones existing as singletons. In parallel, a total of 9,984 random BACs were end-sequenced and, after preliminary screening for length and quality assessments, a total of 18,022 high-quality BESs were generated with success rate of over 90%. Among them, 431 (2.3%) sequences were eliminated after masking vector sequences and sequences contaminated with bacterial and organelle genomes using CROSS_ MATCH and BLASTN searches [28, 29]. In total, 17,591 BESs were generated with an average size of 626 bp. The total BESs length was 11,014,359 bp with a GC content of 44.5%, which suggests that wheat is AT rich, similar to other cereal genomes [10, 30, 31]. The BESs constituted 3.2% of chromosome arm 3AS. The distribution of forward and reverse BESs displayed a 1:1 ratio (9045 F and 8997R). Tracing back the 9,984 end-sequenced BAC clones to the first phase assembly of the 3AS fingerprint map revealed that 6,506 clones are located in 1,071 contigs and 3,688 singletons (Table 1). On an average, six BACs were end-sequenced per contig, theoretically providing coverage of one BES per 45.4 kb across the chromosome arm. All BESs generated herein have been deposited in the GenBank databases under the “GenBank Accession” EI: EI666997-EI676076 and ER772249-ER781190.

Table 1 Origin and distribution of 3AS BAC clones used in BAC end sequencing in the first FPC assembly of chromosome 3AS of bread wheat

Analysis of the repeat fraction and identification of novel repeat elements

The 17,591 BESs accounting for 11 Mb of t3AS were analyzed sequentially for their repeat and gene content using semi-automated pipeline. Based on similarity searches against a repeat database, 79.1% of the nucleotide sequence corresponded to repetitive sequences (Table 2). The class I TEs (retroelements) constitute 66.9% of the sequence, followed by 4.1% for class II TEs (DNA transposons). Within the class I TEs, long terminal repeat (LTR) retrotransposons were the most prominent elements accounting for 98.8%.The most common repeat families of class I were gypsy-like (12,150 reads), followed by copia-like LTR retrotransposons (3,368 reads). Non-LTR retrotransposons, such as LINE (222 reads, 0.7%) and SINE (1 read), were observed at a much lower frequency (Table 2). Clearly, LTR retrotransposons overshadow non-LTR retrotransposons both in number of matches in BESs and percent of sequence. The next most abundant repeats were the class II TEs DNA transposons with 1,626 BESs showing homology. The CACTA family was the most abundant type of DNA transposon accounting for 84.6 % of the class II elements.

Table 2 Occurrence and distribution of repetitive DNA in the BAC end sequences of bread wheat chromosome arm 3AS
Table 3 The frequency and distribution of SSRs in wheat homoeologous chromosomes The SSRs on chromosome arm 3AS were obtained from BESs analysis and SSRs from chromosome 3B similar analysis from BESs [9]

Public repeat databases were used to mask the BESs matching to known DNA repeat sequences. However, the repeat detection could be limited by the size and diversity of the repeat database, such as the absence of genome and species-specific elements. Furthermore, the rapid evolution of repetitive DNA elements may lead to the origin of several unknown repetitive DNA elements. To identify unknown and wheat-specific repetitive DNA elements, the BESs were masked with known repetitive elements and then compared to each other. The BESs with multiple hits (more than eight) at high stringency (at least 100 bp aligned contiguously) were identified as repeat sequences. Families of repeat sequences were obtained computationally from BESs and then subjected to multiple sequence alignment to obtain 122 consensus sequences. The consensus sequences were merged at overlapping regions and extended using CAP3 software [32] and detected 31 putative novel repeats with sizes ranging from 199 bp to 1,420 bp (Additional file 1). The putative repeat elements could be considerably longer but, because the average BES length was 623 bp, only a part of the repeats could be identified. The searching of all databases, including TREP, RepBase and Genbank, against putative repeats resulted in the identification of 11 protein coding genes representing large gene families, 10 of which also matched transposable proteins. The remaining 20 repeat elements showed no similarity in BLASTN or BLASTX [33] and were termed as Triticum novel repeat sequences (TRES). The novel, repeat sequences (TRES) were analyzed for their presence in wheat chromosome 3A by another BLAST against BESs. These TRES were repeated at a higher frequency compared to the initial search in the BESs because of the extended length of the repeat sequences. TRES2 had maximum number of hits in the BESs (84 at e <10-25), and the novel repeats totaled 84 matches in the BESs database (Additional file 2). Seventeen of the 20 TRESs had a hit on chromosome 3B BESs when BLASTN searched in NCBI nt/nr database.

Gene content and functional annotation

After masking the known and putatively new repetitive sequences, the remaining BESs were used to estimate the gene content of chromosome arm 3AS. The sequences were compared with wheat transcript assemblies (Triticum aestivum release 2) using TBLASTX and BLASTN algorithms [33]. We found 519 unique sequences that were similar (e-value 10-50) to an EST with a cumulative match length of 152,662 bp, representing about 1.38 % of the data set. Considering the size of 3AS as 355 Mb, the estimated transcribed fraction of 1.38 % corresponds to 4.8 Mb in 3AS and 11.4 Mb for the entire 3A chromosome. Analysis of full-length cDNAs in several grass species have shown the average lengths of FLcDNA varying between 1,000-2,000 bp. In 37,139 Oryza sativa and 7,455 maize FLcDNAs, the average length was found to be 1,746 ± 945 bp [34] and 1,442 bp [35], respectively. In the Triticeae tribe, the average length of 22,651 Uni-FLcDNAs in barley was 1,711 ± 863 bp [34]. Considering 1.7 kb as the average length of a gene based on the average in barley, our results suggest the presence of about 2,850 genes on the short arm of chromosome 3A and ~6,700 genes on chromosome 3A. This estimate is a little lower than the prediction of 7–8,000 genes in chromosome 3A based on preliminary results from Roche/454 sequences of sorted 3A chromosome arms (our unpublished results).

The gene estimates for the homoeologous 3B chromosome, which is 20% bigger than 3A, varies from 6,000 genes based on BES analysis [9] to over 8,400 to 9,500 genes, based on sequence analysis of 17 Mb of contiguous BAC sequences [36] and ~6,300 genes based on 1.2x coverage of sorted chromosome 3B by Solexa/Illumina sequencing [36]. The estimated number of genes on chromosomes 3A and 3B is much higher than the 4,856 genes predicted for syntenic rice chromosome 1 [30]. However, the predicted gene numbers may be overestimated, because it is difficult to distinguish bona fide genes from pseudogenes or gene fragments due to the incompleteness of genome or chromosome sequence data.

Repeat-masked BESs were functionally annotated with BLASTX searches against the NCBI non-redundant database and functionally classified using Blast2GO [37]. A total of 1,928 BESs had a hit in the database, of which 1,477 sequences (76%) had a definite hit, 109 sequences (5%) hit a transposon and another 342 sequence (18%) hit a hypothetical or predicted protein (Additional file 3). Of 1,928 sequences, only 864 were assigned to GO categories, and an enzyme code was identified for only 31 sequences (Additional file 3). Chromosome 3AS BESs were classified into three categories: biological processes, molecular functions and cellular components. A total of 241 sequences were classified as biological processes, 234 were cellular components and 71 were of metabolic function (Figure 1). Within the biological process classification, the largest categories were cellular processes (39.3% of the sequences) and metabolic processes (30.7% of the sequences), followed by biological regulation (11.6% of the sequences) (Figure 1A). For the cellular components (Figure 1B), 52.8% of the sequences were cell parts and 28.6% were membrane-bounded organelles. Within the molecular functions class (Figure 1C), 46%, 30%, and 22% of the sequences were categorized as transferase activity, hydrolase activity and nucleotide binding, respectively.

Figure 1
figure 1

Distribution of annotated sequences for the chromosome 3A BAC ends under the gene ontology classifications carried out with the automated software Blast2GO. GO-annotated blast hits (1,928) were assigned into the ontology classifications biological processes, cellular components and molecular functions. (A) Biological process. A total of 241 sequences were identified. Categories expressed as a percentage of the total are cellular process (39.4%), metabolic process (30.7%), biological regulation (11.6%), response to stimulus (7.8%) and localization (10.3%). (B) Cellular component. A total of 234 sequences were identified. Categories expressed as a percentage of the total are membrane-bounded organelle (28.6%), cell part (52.9%), organelle part (7.6%) and vesicle (10.6%). (C) Molecular function. A total of 71 sequences were identified. Categories expressed as a percentage of the total are nucleotide binding (22.5%), hydrolase activity (30.9%) and transferase activity (46.4%). The GO classification is based on the biological process class of an associate GO term with an alpha score of at least level 0.6 and ontology level 2 of Blast2GO [37].

Comparative sequence composition of homoeologous A, B and D genomes

The homoeologous genomes of wheat have undergone dramatic sequence diversification in the intergenic regions, but have still retained the backbone of gene synteny with their ancestral genomes. A significant amount of representative sequence data was generated from chromosome 3A in the present study and from chromosome 3B [9, 24] and the D genome [38]. Assuming that a chromosome arm could act as a good representative sample of the genome, we utilized 11 Mb BESs each from chromosome arm 3AS, chromosome 3B [9] and a 2.9 Mb of shotgun sequence from the D genome [38] to compare the composition of homoeologous genomes. Our estimate of the repeat content in chromosome 3AS (79.1%) is similar to the 73-82% repeat content of other A-genome chromosomes sequenced by next generation technology from flow-sorted chromosome-arm specific MDA DNA[19, 20, 22]. However, MDA is known to introduce a bias in the representation of individual amplified regions [39, 40], and a more accurate estimation must await quantitative analysis of the repeat content from flow-sorted chromosome arms DNA without the MDA step. However, BES information from chromosomes 3AS and 3B are the best available sequence samples to compare with homoeologous chromosomes and genomes and, together with next generation sequencing of wheat chromosomes 1A, 4A, 5A, 1B, 7BS, 1D and 7DS, provide a snapshot of the gene content of wheat.

The total repeat fraction in the A genome was 79.1%, which was less than that of the B genome (85.9%) and higher than that of the D genome (74%), indicating a positive correlation between TE content and genome size. The most abundant retrotransposons in all three genomes were gypsy-like, accounting for 33-53% in the A, B and D genomes, followed by copia-like (Figure 2). The A genome has a higher percentage of copia-like elements compared to the B and D genomes, however, the B genome is richest in gypsy-like elements. Interestingly, DNA transposons, such as CACTA elements, are in significantly higher proportion in the D genome compared to the A and B genomes (Figure 2). The Sabrina and Fatima elements seem to be the most abundant and have played a key role in expansion of the A genome along with other elements, such as Erika, Angela and Wilma, which form a considerable part of the retroelements in the genome (Additional file 4). For chromosome 5A, Sabrina (18%), Wilma (8.2%) and Fatima (6.3%) were the most abundant repeats [20]. The Sabrina, Angela, Barbara,Vagabond, Egug and Erika elements were significantly higher (p-value < 0.0001) in the A genome compared to the B genome based on BESs analysis. The ratio of class I elements to class II elements in the A (16.2) and B (12.3) genomes is quite similar, whereas it is significantly lower (4.1) in the D genome, indicating the significant role of DNA transposons in D-genome expansion.

Figure 2
figure 2

Histogram showing the comparative abundance of different retrotransposons and DNA transposon families (data expressed in per cent of genome) in the three homoeologous genomes of wheat. The BAC end sequences for chromosome arm 3AS (A genome) and comparison with the B and D genomes. The data representative for the B and D genomes was obtained in a similar way, using 10.8 Mb of BESs from wheat chromosome 3B [9], and the D genome was represented by 2.9 Mb sequence obtained from whole-genome, shotgun sequence form Ae. tauschii the D-genome donor of bread wheat [38].

Based on the analysis of 3B BESs and the reanalysis of the D-genome shotgun sequence [38], Paux et al.[9] hypothesized that the Triticeae genomes have evolved from an ancestral genome of approximately 1 Gb. The portion of the genome not related to the repetitive sequences in the A genome based on 3AS and 5A sequence data is estimated at 1.1 Gb (~20% of the 5.5 Gb A genome). However, most probably, this is an overestimate, because the wheat repeated sequence database is incomplete and the non-repetitive portion of the D genome based on evidence from bioinformatics and hybridization analyses was just 8.4% [38].

GC content is an important feature and was quite comparable between chromosomes 3AS (44.5%) and the other A-genome chromosomes 1AL (44.7%, [14]) 4AS (44.7%, [22]) however, the GC content of chromosome 4AL (41.4%, [22]) was considerably lower possibly because of biasness caused by MDA. The GC content in the B-genome chromosome 3B (44.6%, [9]) was quite similar to that of the A-genome chromosomes of wheat, and the GC content of wheat is comparable to other grasses, such as rye (45.9%, [10]), maize (46.5%, [31]) and rice (44.0%, [30]).

The BESs from chromosome 3AS were BLASTN against the 3,000-loci transcript map recently developed for chromosome 3B [24], and 100 BESs matched to 77 unigenes from chromosome 3B. Of these, 77 BESs (matching 77 unigenes), 45 had a hit in rice (40/77) or Brachypodium (40/77) with 70% (28/40) of the hits in rice and 77% (31/40) in Brachypodium in syntenic regions (Additional file 5). Rustenholz et al. [24] reported 35-42% of unigenes genes mapped on chromosome 3B were non-syntenic in model grasses genomes. Our results, based on the sequence from 3A together with 3B unigenes [24] showed that 23-30% of the genes are non-syntenic in wheat, which is consistent with data from chromosome arms 1AL, 1BL, and 1DL that showed 20-45% were non-syntenic [19]. A similar rate of colinearity was observed for chromosomes 7BS and 7DS with the corresponding region of the Brachypodium by sequence analysis of sorted chromosome arms [21]. This erosion of colinearity could be explained by a high rate of transposable element activity in hexaploid wheat, such as TE-driven gene movement [19]. Nearly, 40% (32/77) BESs having a hit to 3B unigenes had no hit in rice or Brachypodium. Rustenholz et al. [24] reported that 89% of these unigenes were expressed, and we found 11 of the 32 BESs had a hit with barley full-length cDNA suggesting these genes are less likely to be pseudogenes. Most probably, these genes are less conserved and, therefore, remain undetected in model genomes, however, a few of them could be novel wheat/Triticeae-specific genes. In the recent barley genome sequence analysis by Meyer et al. [18], similar results were observed where 23% of the barley genes had no sequence similarity in model genomes, although they were expressed in barley.

Based on our BES analysis, we estimated one gene per 123 kb and ~44,600 genes in the A genome of wheat. Vitulo et al. [20] recently sequenced the MDA from flow-sorted chromosome 5A of wheat and estimated 5,088 genes of chromosome 5A and with a gene density of one gene per 162 kb. The coding fraction for chromosome 3AS (1.38%) was comparable to those of chromosome 5A (1.08-1.30%, [20]) and 1AL (1.03%, [14]). In another sequence analysis and comparison between the chromosome 4A of wheat and model genome Brachypodium, Hernandez et al. [22] estimated 9,500 genes on chromosome 4A and ~61,500 in the A genome. However, taking into account the large number of pseudogenes and paralogs, they estimated a minimum number of 4,300 genes on chromosome 4A and 28,000 in the A genome of wheat [22]. The number of genes estimated in our analysis is higher than conservative estimates based on chromosomes 4A and 5A, however, it corresponds well with the recently estimated numbers of genes in the A genome based on BESs from chromosome 1AL (50,000 genes [14]) and the number of genes in the B genome based on 17 Mb of BAC sequence of chromosome 3B (41,000-50,000 genes [36]). These estimates are much lower than predicted in earlier studies based on the random sequencing of methylfiltrated libraries [41] and significantly higher than the estimated 36,000 genes for the D-genome donor Ae. tauschii based on BAC sequencing [42].

Berkman et al. [21] made a very stringent estimate of about 25,000-26,000 genes per genome of wheat based on sequence analysis of wheat chromosome arms 7BS and 7DS. However, the estimates made by Berkman et al. [21] and Hernandez et al. [22] are based on the genome zipper approach. These estimates do not take into account about 23-30% of the genes that lack sufficient sequence similarity to any gene in the three model grass genomes as suggested by Meyer et al. [18] in barley and our analysis of 3AS BESs and 3B unigenes in wheat.

On the other hand, the estimate of gene content by several other studies show that the relatively higher gene content in wheat (100,000-350,000) [9, 41] could be an inflation in the actual number of genes due to misannotations of pseudogenes as all these estimates are based on partial gene sequences or a relatively small data set compared to the large wheat genome. Furthermore, recent sequence analysis of the group-1 chromosomes of wheat showed that a large number of nonsysntenic genes are present in only one of three wheat genomes and only 60 % of them showed transcriptional evidence, suggesting that most of them are pseudogenes [19] making annotation of genes based on sequence similarity suboptimal. We may have a better estimate of the gene content once full-length gene sequences and a gold standard reference sequence of wheat becomes available allowing us to differentiate between genes and gene fragments.

Comparison to model grass genomes

Repeat-masked BESs were used for a comparative analysis with other model grass genomes. The 3AS BESs were subjected to a similarity search against rice chromosome 1, Brachypodium chromosome 2 and Sorghum chromosome 3 and also to chromosome 3B BESs using BLASTN. There were 193 hits to rice chromosome 1, 346 on Brachypodium chromosome 2 and 185 to Sorghum chromosome 3 (Figure 3, Additional file 6). There were 632 hits to 3B BESs demonstrating a high level of homology, as expected. The 3A BESs, when aligned to a model grass genome based on sequence similarity and ordered based on Brachypodium, showed a mosaic of synteny (Figure 3). The short arm of chromosome 3A showed a high level of synteny with short arm of rice chromosome 1, Sorghum chromosome 3 and Brachypodium chromosome 2 (Figure 3). Sequence hits were distributed all along chromosome 3AS with a mean of 0.9 hits/Mb of sequence and range of 0–27 hits/Mb of sequence, however, the density of hits were higher towards the distal end of the chromosome, which could be due to higher gene density in terminal part of the chromosome [24, 43]. A similar trend of increased gene density towards the distal end of the chromosome was seen in chromosome 3B [24]. The BESs hit 43 bins (bin size = 1 Mb) in rice chromosome 1 with a mean of 4.4 hits per bin and range of 0–21 hits per Mb of sequence. In Brachypodium chromosome 2, the mean number of hits per Mb of sequence was 6.9 with range of 0–31 hits/Mb of sequence. In Sorghum chromosome 3, the mean hits per Mb of sequence were 2.4. Our comparative analysis validates the high synteny (70-77%) between the short arm of chromosome 3A with the short arms of rice chromosome 1, Brachypodium chromosome 2 and Sorghum chromosome 3 (Figure 3).

Figure 3
figure 3

High-resolution comparative analysis between wheat chromosome 3AS and rice chromosome 1, Brachypodium Bd2 and Sorghum chromosome 3. High-density comparative analysis chromosome 3AS BAC end sequences versus the sequenced model grass genomes of rice, Brachypodium and Sorghum. The figure includes two concentric circles. The inner circle represents wheat chromosome 3AS (green) and orthologous chromosomes from rice (blue), Brachypodium (yellow) and Sorghum (red). The outer circle illustrates the heat map and colored bars indicate the density of BLAST hits per mega base of sequence (bin size = 1 Mb). Colored lines in the center represent the putative orthologous relationship between the wheat chromosome 3AS and rice chromosome 1, Brachypodium Bd2 and Sorghum chromosome 3.

Molecular markers from 3AS

In order to identify new putative markers for physical and genetic mapping, 17,591 BESs were analyzed for the presence of SSRs, which could be converted into microsatellite markers. A total of 1,057 perfect repeats of di-, tri- and tetranucleotide repeats were identified. The frequency of SSRs derived from the BESs was about one SSR per 10.4 kb (Table 3). The most abundant of these SSRs were dinucleotide 779 (73.8%) and trinucleotide 268 (25.2%); tetranucleotides were rare (1.0%). Nearly 28% of SSRs (both di- and trinucleotide) were AT rich (i.e., TA, AT and TAA). A comparison of the frequency and distribution of the SSRs between the A and B genomes of wheat indicates a higher density of SSRs on chromosome 3B (one SSR per 6.1 kb) compared to chromosome 3AS (one SSR per 10.4 kb). Furthermore, the dinucleotide motif (TA)n was most abundant (12.6%) on chromosome 3AS in contrast to the trinucleotide motif (AAG), which was the most abundant in wheat chromosome 3B [9]. Other species, such as Brachypodium, rice and maize, also have a higher frequency of trinucleotide motifs, however, species such as papaya show a higher frequency of dinucleotide motifs and soybean had higher frequency of tetranucleotide motifs [44]. In our study, the frequency of SSRs on chromosome 3AS is far lower than that of chromosome 3B. The higher frequency of AT-rich di- and trinucleotide motifs in wheat and maize suggests that AT-rich region are associated with repetitive sequences, in contrast to GC-rich regions, which are associated with the transcribed part of the genome.

Microsatellite primer pairs were successfully developed from 758 out of 1,057 SSRs. Of the 758 SSRs primer pairs, 155 were developed from class I microsatellites with more than nine motif repeats (Additional file 7). The BESs were also processed to identify junctions between TE and other sequences using ISBPfinder.pl [45]. In total, 7,928 insertion sites were identified in 17,591 BESs, and primer pairs were successfully designed from 695 suitable junctions (Additional file 7). A repeat was identified every 1.39 kb on t3AS, whereas a repeat was identified every 0.89 kb on t1AL [13], suggesting that repeat-based makers could play an important role in the saturation mapping of the wheat genome. A set of 96 class I SSRs and 96 ISBP primer pairs were tested by PCR amplification for 3AS specificity and utility across species.

For testing 3AS specificity, primer pairs were amplified on genomic DNA of a 3AS BAC pool, Chinese Spring (CS), and the CS aneuploids N3AT3B (nullisomic for 3A and tetrasomic for 3B), Dt3AS (nullisomic 3AL) and Dt3AL (nullisomic 3AS). An ISBP primer pair was considered 3AS specific if it did not amplify a product from N3AT3B or Dt3AL but amplified products from all other stocks (Figure 4). Nearly 18% of the SSRs did not produce amplicons from N3AT3B or Dt3AL suggesting these markers had 3AS specificity (Figure 4). Twenty-nine percent of the ISBP primer pairs did not amplify from N3ATB or Dt3AL indicating their specificity to 3AS (Figure 4). However, all the ISBP and ISSR markers that amplified from CS (93% SSRs and 98% ISBPs) also amplified from the CS t3AS BAC library super pools, suggesting that they are present on chromosome 3AS. The nonspecificity of nearly two-thirds of the ISBP and four-fifths of the SSRs markers may be attributed to a combination of several factors including amplification from the homoeologous chromosomes, t3AS has a portion of the long arm of chromosome 3A, flow-sorted chromosome arm-based libraries have a reported 10-13 % contamination [3]; duplications on other wheat chromosomes, because microsatellites have a high possibility of being moved by transposable elements; or bias in the sample tested.

Figure 4
figure 4

Example of the localization SSR and ISBP markers on chromosome 3AS and amplification across species. PCR products separated on 3 % agarose gels after amplification with one SSR (Xbem5) and one ISBP (Xisbp3A-14). Templates were DNA from (1) CS wheat, (2) a 3AS BAC library superpool, (3) N3AT3B, (4) Dt3AS, (5) Dt3AL, (6) T. monococcum subsp. monococcum TA4342-L96, (7) T. monococcum subsp. aegilopoides TA4342-L95, (8) T. turdigum subsp. durum cv. Langdon, (9) Langdon DIC-3A, (10) T. aestivum cv. Opata 85, (11) synthetic wheat W79849, (12) water, (L) 100-bp DNA ladder.

Nearly, 53% (23/44) of ISBPs designed from 1AL BESs were 1AL specific [13], suggesting that because the first primer is designed from repeat and the second one from unique sequences, the high specificity of the unique sequences becomes very significant for specific amplification.

The set of 96 class I SSRs and 96 ISBP primer pairs were also tested by PCR amplification on wheat genotypes, the parents of three mapping populations: cultivated diploid wheat (T. momococcum subsp. monococcum) and wild diploid wheat (T. monococcum subsp. aegilopoides); tetraploid wheat (T. turgidum subsp. durum) genotypes Langdon and Langdon DIC-3A; and bread wheat genotypes CS, Opata 85 and synthetic W7984 (Figure 4). More than 93% of the SSR and 97% of ISBP primer pairs gave amplification from diploid, tetraploid and hexaploid wheat. Twenty-six percent (25/96) ISBPs and 16% (15/96) SSRs were found polymorphic between parents of diploid mapping population T. monococcum and T. aegilopoides on 3% agarose gels. To eliminate the effect of amplification from 3B and 3D, we are mapping these markers in a diploid wheat (AA) T. monococcum subsp. monococcum/subsp. aegilopoides RIL population (data not shown). Amplification of the primer pairs from diploid, tetraploid and hexaploid wheat shows their reproducibility and wide adaptability of these markers.

In summary, a total of 1,453 new (758 SSR + 695 IBSP) putative markers have been developed spanning the 3AS arm, which adds over 200% more markers than mapped in all previous studies (<600 markers) on this arm. The universality and variability of SSRs make these sequences an attractive source of developing microsatellite markers. However, because most of the genome (>80%) is repetitive, the ISBP repeat insertion-based markers could be very important in saturating markers in this portion of the genome and, thus, could enhance the efficiency of breeding programs. The potential of such repeat-based markers in crop breeding for high-throughput, marker-assisted selection has been suggested by Paux et al. [45]. More than 30 genes and QTLs related to domestication traits, such as brittle rachis; tillering; disease resistance genes against stem rust, leaf rust, tanspot, and Septoria tritici blotch; and agronomic traits such as preharvest sprouting and root-shoot biomass, have been genetically mapped on chromosome 3AS of wheat (http://www.shigen.nig.ac.jp/wheat/komugi/genes/symbolClassList.jsp) [46]. In addition to the 1,453 molecular markers (SSR and ISBP) developed in the present study, we have obtained sequence tags of 161 of the estimated ~2,850 genes on chromosome 3AS. The 1,453 molecular markers plus 161 sequence tags of genes will aid in marker-assisted selection and the map-based cloning of above-mentioned agronomic traits on 3AS arm. The molecular markers will also help to align the physical and genetic maps of chromosome 3AS.

Conclusion

BESs are a valuable genomic resource for genome organization analysis, comparative genomics and marker development. Analysis of 17,591 3AS BESs revealed differences in the repeat landscape of the A- and B-genome homoeologous chromosomes and estimate of gene content, Triticum-specific and genome-specific genes. We also developed BESs-derived microsatellite and ISBP putative markers as a resource for integrating physical and genetic maps and for mapping and cloning of a large number of agronomically important genes on 3AS arm.

Methods

BAC library

We used the first of the two BAC libraries constructed from the short arm of chromosome 3A (3AS). The library (code TaaCsp3AShA) consists of 55,296 clones arranged in 144 384-well plates. With the average insert size of 80 kb and 11% contamination with other chromosomes, the library represents 10.9x equivalents of the chromosome arm [25].

BAC end sequencing

Twenty-six 384-well plates (9,984 BAC clones) randomly selected from the 3AS library were end sequenced (5’ & 3’) at the J. Craig Venter Institute (Rockville, MD). BAC clones of HindIII 96-deep-well blocks containing 1.2 ml/well of 2xYT medium were grown overnight at 37°C with shaking at 300 rpm. The cells were harvested by centrifugation, and the BAC DNA was purified using a REAL Prep 96 Plasmid Kit (Qiagen, Valencia, CA). For BAC-end sequencing, 5 μl of purified BAC DNA (~ 0.2 to 0.5 μg) was used in a sequencing reaction with ABI BigDye terminators (Applied Biosystems, Foster City, CA). Template DNA was sequenced from both directions with pCC1BAC/pIndigoBAC-5 Forward and Reverse End-Sequencing Primers (Epicentre, Madison, WI). Electrophoresis of the sequencing reaction was carried out with a 3730xl DNA Analyzer (Applied Biosystems, Foster City, CA).

Sequence data processing

The sequences and quality files from trace files were read by the Phred program for base calling and trimmed to remove vector and low quality bases [28, 29]. The high-quality sequence data were then filtered for sequences contaminated with Escherichia coli or with plant organelle genomes based on matches to the wheat mitochondria and chloroplast sequences.

Sequence analysis

The high-quality BESs were compared with several repeat databases, including the Triticeae repetitive (TREP) sequence database (http://wheat.pw.usda.gov/ITMI/Repeats/)[47, the TIGR plant repeat database (ftp://ftp.tigr.org/pub/data/TIGR_Plant_Repeats/), and RepBase at GIRI. (http://www.girinst.org/)[48]. Self-BLASTN and BLASTX at an e-value cutoff of 10−5 was used [33]. The composition and contents of repeat element in BESs was assessed using the RepeatMasker program (http://www.repeatmasker.org/) with the default settings using each of the above databases as the custom library option. Finally, a cutoff score of 250 was used to identify the repetitive sequences. The BESs were annotated based on their best match to the repeat database and categorized according to (http://wheat.pw.usda.gov/ITMI/Repeats/gene_annotation.pdf)[47]. Sequences matching known repeats were masked as N and repeat masked sequences were further used to identify novel repeats. Self-BLASTN was performed on repeat-masked BESs to identify sequences that had multiple strong matches to other BESs with an e-value <10−50. The blast output file was processed using RECON software (http://selab.janelia.org/recon.html) [49] to identify families with eight or more sequences. The sequences were extracted and aligned by ClustalW [50]. Consensus sequences for each family with a minimum 100-bp alignment and 80 % identity were selected. The consensus sequences were compared to each other and aligned using the CAP3 software [32] to merge overlapping regions and extend the sequences. The putative novel repeats were searched in BLAST against nonredundant nucleotide, EST and protein databases to remove any characterized sequences. The repeat-masked sequences were subjected to gene content analysis by a homology search using TBLASTX versus the T. aestivum Transcript assembly 2.0 (e-value 10−50) (http://www.jcvi.org/wheat/downloads.php). The cumulative match lengths were used to calculate the coding fraction, as described for repetitive elements. The repeat masked sequences were subjected to Blast2GO for functional annotation using default parameters [37].

Comparative sequence analysis

The rice (http://rice.plantbiology.msu.edu/), Brachypodium (http://ftp.brachypodium.org/) and Sorghum (http://www.phytozome.net/sorghum) sequences were downloaded from their respective databases. The 3B BESs were obtained from NCBI Genbank and 3B unigene sequences were requested from authors [24]. The repeat-masked BESs sequences were searched in model grass genomes by BLASTN against Brcahypodium chromosome 2, rice chromosome 1 and Sorghum chromosome 3 at an e-value of 1e-5 and a minimum hit length of 50 bp. The unique hits for each BESs were used for developing a circos map for visualization of synteny [51]. The 3AS BESs were compared with chromosome 3B unigenes, barley full length cDNA by BLASTN at an e-value of 1e-10.

Development of molecular markers from BESs

The BESs were screened computationally for the presence of di-, tri- and tetranucleotides SSRs using the SSRIT program (ftp://ftp.gramene.org/pub/gramene/software/scripts/ssr.pl). SSRs motifs that span ten or more nucleotides were recorded, and SSR location data in sequence were extracted from the SSR search output files. Sequences containing SSR were used to design primers using an in-house perl script and Pimer 3.0 with default settings, and ISBP primers were designed using IsbpFinder.pl [45]. All the SSR and ISBP primer pairs were amplified using the 0.8 μM of primer and 0.25 mM of dNTP, 1.25 mM MgCl2, 0.1 U Taq polymerase (Bioline, CA, USA) with 1× buffer and 50 ng template DNA. PCR amplification was conducted using the following procedure: nine cycles (95°C for 30 sec, 58°C for 60 sec minus 0.5°C each cycle, 72°C for 60 sec), 25 cycles (95°C for 30 sec, 53°C for 45 sec, 72°C for 60 sec) and one additional cycle (72°C for 7 minutes). Genomic DNA amplified with SSR and ISBP primers were separated on 3 % GenePure HiRes Agarose (ISC Bioexpress, Kaysville, UT, USA) gels.

Author’s contributions

BSG and WL: conceived the study; SKS and WL: carried out the sequence data analysis; PR and AC: BAC-end sequencing; HS and JD: provided the BAC library and participated in manuscript preparation; SKS, BSG and WL wrote the manuscript. All authors read and approved the final manuscript.

Abbreviations

BAC:

bacterial artificial chromosome

BESs:

BAC end sequences

EST:

expressed sequence tag

GO:

gene ontology

HSP:

high-scoring segment pair

IB:

isolation buffer

Gb:

Gigabase

kb:

kilobase

Mb:

Megabase

bp:

base pair

MDA:

multiple displacement amplification

SSR:

simple sequence repeat

ISBP:

insertion site based polymorphism

RIL:

recombinant inbred line..

References

  1. Gill BS, Appels R, Botha-Oberholster AM, Buell CR, Bennetzen JL, Chalhoub B, Chumley F, Dvorak J, Iwanaga M, Keller B, Li WL, McCombie WR, Ogihara Y, Quetier F, Sasaki T: A workshop report on wheat genome sequencing: international genome research on wheat consortium. Genetics. 2004, 168 (2): 1087-1096. 10.1534/genetics.104.034769.

    Article  PubMed  Google Scholar 

  2. Kubalakova M, Vrana J, Cihalikova J, Simkova H, Dolezel J: Flow karyotyping and chromosome sorting in bread wheat (Triticum aestivumL.). Theor Appl Genet. 2002, 104 (8): 1362-1372. 10.1007/s00122-002-0888-2.

    Article  PubMed  Google Scholar 

  3. Dolezel J, Kubalakova M, Paux E, Bartos J, Feuillet C: Chromosome-based genomics in the cereals. Chromosome Res. 2007, 15 (1): 51-66. 10.1007/s10577-006-1106-x.

    Article  PubMed  CAS  Google Scholar 

  4. Safar J, Simkova H, Kubalakova M, Cihalikova J, Suchankova P, Bartos J, Dolezel J: Development of chromosome-specific BAC resources for genomics of bread wheat. Cytogenet Genome Res. 2010, 129 (1–3): 211-223.

    Article  PubMed  CAS  Google Scholar 

  5. Safar J, Bartos J, Janda J, Bellec A, Kubalakova M, Valarik M, Pateyron S, Weiserova J, Tuskova R, Cihalikova J, Vrana J, Simkova H, Faivre-Rampant P, Sourdille P, Caboche M, Bernard M, Dolezel J, Chalhoub B: Dissecting large and complex genomes: flow sorting and BAC cloning of individual chromosomes from bread wheat. Plant J. 2004, 39 (6): 960-968. 10.1111/j.1365-313X.2004.02179.x.

    Article  PubMed  CAS  Google Scholar 

  6. Paux E, Sourdille P, Salse J, Saintenac C, Choulet F, Leroy P, Korol A, Michalak M, Kianian S, Spielmeyer W, Lagudah E, Somers D, Kilian A, Alaux M, Vautrin S, Berges H, Eversole K, Appels R, Safar J, Simkova H, Dolezel J, Bernard M, Feuillet C: A physical map of the 1-gigabase bread wheat chromosome 3B. Science. 2008, 322 (5898): 101-104. 10.1126/science.1161847.

    Article  PubMed  CAS  Google Scholar 

  7. Messing J, Bharti AK, Karlowski WM, Gundlach H, Kim HR, Yu Y, Wei FS, Fuks G, Soderlund CA, Mayer KFX, Wing RA: Sequence composition and genome organization of maize. Proc Natl Acad Sci USA. 2004, 101 (40): 14349-14354. 10.1073/pnas.0406163101.

    Article  PubMed  CAS  Google Scholar 

  8. Hong CP, Plaha P, Koo D-H, Yang T-J, Choi SR, Lee YK, Uhm T, Bang J-W, Edwards D, Bancrofts I, Park B-S, Lee J, Lim YP: A survey of theBrassica rapagenome by BAC-end sequence analysis and comparison withArabidopsis thaliana. Mol Cells. 2006, 22 (3): 300-307.

    PubMed  Google Scholar 

  9. Paux E, Roger D, Badaeva E, Gay G, Bernard M, Sourdille P, Feuillet C: Characterizing the composition and evolution of homoeologous genomes in hexaploid wheat through BAC-end sequencing on chromosome 3B. Plant J. 2006, 48 (3): 463-474. 10.1111/j.1365-313X.2006.02891.x.

    Article  PubMed  CAS  Google Scholar 

  10. Bartos J, Paux E, Kofler R, Havrankova M, Kopecky D, Suchankova P, Safar J, Simkova H, Town CD, Lelley T, Feuillet C, Dolezel J: A first survey of the rye (Secale cereale) genome composition through BAC end sequencing of the short arm of chromosome 1R. BMC Plant Biol. 2008, 8 (95): .

  11. Datema E, Mueller LA, Buels R, Giovannoni JJ, Visser RGF, Stiekema WJ, van Ham RCHJ: Comparative BAC end sequence analysis of tomato and potato reveals overrepresentation of specific gene families in potato. BMC Plant Biol. 2008, 8 (34): .

  12. Gonzalez VM, Rodriguez-Moreno L, Centeno E, Benjak A, Garcia-Mas J, Puigdomenech P, Aranda MA: Genome-wide BAC-end sequencing ofCucumis melousing two BAC libraries. BMC Genomics. 2010, 11 (618): .

  13. Dolezel J, Greilhuber J, Suda J: Estimation of nuclear DNA content in plants using flow cytometry. Nat Protoc. 2007, 2 (9): 2233-2244. 10.1038/nprot.2007.310.

    Article  PubMed  CAS  Google Scholar 

  14. Lucas S, Simkova H, Safar J, Jurman I, Cattonaro F, Vautrin S, Bellec A, Berges H, Dolezel J, Budak H: Functional features of a single chromosome arm in wheat (1AL) determined from its structure. Funct Integr Genomics. 2011, 12 (1): 173-182.

    Article  PubMed  Google Scholar 

  15. Shultz JL, Kazi S, Bashir R, Afzal JA, Lightfoot DA: The development of BAC-end sequence-based microsatellite markers and placement in the physical and genetic maps of soybean. Theor Appl Genet. 2007, 114 (6): 1081-1090. 10.1007/s00122-007-0501-9.

    Article  PubMed  CAS  Google Scholar 

  16. Mayer KFX, Taudien S, Martis M, Simkova H, Suchankova P, Gundlach H, Wicker T, Petzold A, Felder M, Steuernagel B, Scholz U, Graner A, Platzer M, Dolezel J, Stein N: Gene content and virtual gene order of barley chromosome 1 H. Plant Physiol. 2009, 151 (2): 496-505. 10.1104/pp.109.142612.

    Article  PubMed  CAS  Google Scholar 

  17. Berkman PJ, Skarshewski A, Lorenc MT, Lai K, Duran C, Ling EYS, Stiller J, Smits L, Imelfort M, Manoli S, McKenzie M, Kubalakova M, Simkova H, Batley J, Fleury D, Dolezel J, Edwards D: Sequencing and assembly of low copy and genic regions of isolatedTriticum aestivumchromosome arm 7DS. Plant Biotech J. 2011, 9 (7): 768-775. 10.1111/j.1467-7652.2010.00587.x.

    Article  CAS  Google Scholar 

  18. Mayer KFX, Martis M, Hedley PE, Simkova H, Liu H, Morris JA, Steuernagel B, Taudien S, Roessner S, Gundlach H, Kubalakova M, Suchankova P, Murat F, Felder M, Nussbaumer T, Graner A, Salse J, Endo T, Sakai H, Tanaka T, Itoh T, Sato K, Platzer M, Matsumoto T, Scholz U, Dolezel J, Waugh R, Stein N: Unlocking the barley genome by chromosomal and comparative genomics. Plant Cell. 2011, 23 (4): 1249-1263. 10.1105/tpc.110.082537.

    Article  PubMed  CAS  Google Scholar 

  19. Wicker T, Mayer KFX, Gundlach H, Martis M, Steuernagel B, Scholz U, Simkova H, Kubalakova M, Choulet F, Taudien S, Platzer M, Feuillet C, Fahima T, Budak H, Dolezel J, Keller B, Stein N: Frequent gene movement and pseudogene evolution is common to the large and complex genomes of wheat, barley, and their relatives. Plant Cell. 2011, 23 (5): 1706-1718. 10.1105/tpc.111.086629.

    Article  PubMed  CAS  Google Scholar 

  20. Vitulo N, Albiero A, Forcato C, Campagna D, Dal Pero F, Bagnaresi P, Colaiacovo M, Faccioli P, Lamontanara A, Simkova H, Kubalakova M, Perrotta G, Facella P, Lopez L, Pietrella M, Gianese G, Dolezel J, Giuliano G, Cattivelli L, Valle G, Stanca AM: First Survey of the wheat chromosome 5A composition through a next generation sequencing approach. PLoS One. 2011, 6 (10): e26421-10.1371/journal.pone.0026421.

    Article  PubMed  CAS  Google Scholar 

  21. Berkman PJ, Skarshewski A, Manoli S, Lorenc MT, Stiller J, Smits L, Lai K, Campbell E, Kubalakova M, Simkova H, Batley J, Dolezel J, Hernandez P, Edwards D: Sequencing wheat chromosome arm 7BS delimits the 7BS/4AL translocation and reveals homoeologous gene conservation. Theor Appl Genet. 2012, 124 (3): 423-432. 10.1007/s00122-011-1717-2.

    Article  PubMed  CAS  Google Scholar 

  22. Hernandez P, Martis M, Dorado G, Pfeifer M, Galvez S, Schaaf S, Jouve N, Simkova H, Valarik M, Dolezel J, Mayer KFX: Next-generation sequencing and syntenic integration of flow-sorted arms of wheat chromosome 4A exposes the chromosome structure and gene content. Plant J. 2012, 69 (3): 377-386. 10.1111/j.1365-313X.2011.04808.x.

    Article  PubMed  CAS  Google Scholar 

  23. Rustenholz C, Hedley PE, Morris J, Choulet F, Feuillet C, Waugh R, Paux E: Specific patterns of gene space organisation revealed in wheat by using the combination of barley and wheat genomic resources. BMC Genomics. 2010, 11 (714): .

  24. Rustenholz C, Choulet F, Laugier C, Safar J, Simkova H, Dolezel J, Magni F, Scalabrin S, Cattonaro F, Vautrin S, Bellec A, Berges H, Feuillet C, Paux E: A 3,000-loci transcription map of chromosome 3B unravels the structural and functional features of gene islands in hexaploid wheat. Plant Physiol. 2011, 157 (4): 1596-1608. 10.1104/pp.111.183921.

    Article  PubMed  CAS  Google Scholar 

  25. Safar J, Simkova H, Kubaláková M, Suchankova P, Cihalikova J, Bartos J, Fiocchetti F, Roselli M, Gill BS, Dolezel J, Lucretti S: Generating resources for genomics of wheat homoeologous chromosome group 3: 3AS- and 3DS-specific BAC libraries. J Genet Breed. 2007, 61: 151-160.

    Google Scholar 

  26. Gill BS, Friebe B, Endo TR: Standard karyotype and nomenclature system for description of chromosome bands and structural-aberrations in wheat (Triticum aestivum). Genome. 1991, 34 (5): 830-839. 10.1139/g91-128.

    Article  Google Scholar 

  27. Soderlund C, Longden I, Mott R: FPC: a system for building contigs from restriction fingerprinted clones. CABIOS. 1997, 13: 523-535.

    PubMed  CAS  Google Scholar 

  28. Lazo GR, Tong J, Miller R, Hsia C, Rausch C, Kang Y, Anderson OD: Software scripts for quality checking of high-throughput nucleic acid sequencers. Biotechniques. 2001, 30 (6): 1300.

    PubMed  CAS  Google Scholar 

  29. Chou HH, Holmes MH: DNA sequence quality trimming and vector removal. Bioinformatics. 2001, 17 (12): 1093-1104. 10.1093/bioinformatics/17.12.1093.

    Article  PubMed  CAS  Google Scholar 

  30. Goff SA, Ricke D, Lan TH, Presting G, Wang RL, Dunn M, Glazebrook J, Sessions A, Oeller P, Varma H, Hadley D, Hutchinson D, Martin C, Katagiri F, Lange BM, Moughamer T, Xia Y, Budworth P, Zhong JP, Miguel T, Paszkowski U, Zhang SP, Colbert M, Sun WL, Chen LL, Cooper B, Park S, Wood TC, Mao L, Quail P, Wing R, Dean R, Yu YS, Zharkikh A, Shen R, Sahasrabudhe S, Thomas A, Cannings R, Gutin A, Pruss D, Reid J, Tavtigian S, Mitchell J, Eldredge G, Scholl T, Miller RM, Bhatnagar S, Adey N, Rubano T, Tusneem N, Robinson R, Feldhaus J, Macalma T, Oliphant A, Briggs S: A draft sequence of the rice genome (Oryza sativa L. ssp japonica). Science. 2002, 296 (5565): 92-100. 10.1126/science.1068275.

    Article  PubMed  CAS  Google Scholar 

  31. Haberer G, Young S, Bharti AK, Gundlach H, Raymond C, Fuks G, Butler E, Wing RA, Rounsley S, Birren B, Nusbaum C, Mayer KFX, Messing J: Structure and architecture of the maize genome. Plant Physiol. 2005, 139 (4): 1612-1624. 10.1104/pp.105.068718.

    Article  PubMed  CAS  Google Scholar 

  32. Huang XQ, Madan A: CAP3: A DNA sequence assembly program. Genome Res. 1999, 9 (9): 868-877. 10.1101/gr.9.9.868.

    Article  PubMed  CAS  Google Scholar 

  33. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.

    Article  PubMed  CAS  Google Scholar 

  34. Matsumoto T, Tanaka T, Sakai H, Amano N, Kanamori H, Kurita K, Kikuta A, Kamiya K, Yamamoto M, Ikawa H, Fujii N, Hori K, Itoh T, Sato K: Comprehensive sequence analysis of 24,783 barley full-length cdnas derived from 12 clone libraries. Plant Physiol. 2011, 156 (1): 20-28. 10.1104/pp.110.171579.

    Article  PubMed  CAS  Google Scholar 

  35. Soderlund C, Descour A, Kudrna D, Bomhoff M, Boyd L, Currie J, Angelova A, Collura K, Wissotski M, Ashley E, Morrow D, Fernandes J, Walbot V, Yu Y: Sequencing, mapping, and analysis of 27,455 maize full-length cDNAs. PLoS Genetics. 2009, 5 (11): e1000740-10.1371/journal.pgen.1000740.

    Article  PubMed  Google Scholar 

  36. Choulet F, Wicker T, Rustenholz C, Paux E, Salse J, Leroy P, Schlub S, Le Paslier MC, Magdelenat G, Gonthier C, Couloux A, Budak H, Breen J, Pumphrey M, Liu SX, Kong XY, Jia JZ, Gut M, Brunel D, Anderson JA, Gill BS, Appels R, Keller B, Feuillet C: Megabase level sequencing reveals contrasted organization and evolution patterns of the wheat gene and transposable element spaces. Plant Cell. 2010, 22 (6): 1686-1701. 10.1105/tpc.110.074187.

    Article  PubMed  CAS  Google Scholar 

  37. Gotz S, Garcia-Gomez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, Robles M, Talon M, Dopazo J, Conesa A: High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008, 36 (10): 3420-3435. 10.1093/nar/gkn176.

    Article  PubMed  CAS  Google Scholar 

  38. Li WL, Zhang P, Fellers JP, Friebe B, Gill BS: Sequence composition, organization, and evolution of the core Triticeae genome. Plant J. 2004, 40 (4): 500-511. 10.1111/j.1365-313X.2004.02228.x.

    Article  PubMed  CAS  Google Scholar 

  39. Pinard R, de Winter A, Sarkis GJ, Gerstein MB, Tartaro KR, Plant RN, Egholm M, Rothberg JM, Leamon JH: Assessment of whole genome amplification-induced bias through high-throughput, massively parallel whole genome sequencing. BMC Genomics. 2006, 7 (216): .

  40. Pugh TJ, Delaney AD, Farnoud N, Flibotte S, Griffith M, Li HI, Qian H, Farinha P, Gascoyne RD, Marra MA: Impact of whole genome amplification on analysis of copy number variants. Nucleic Acids Res. 2008, 36 (13): e80-10.1093/nar/gkn378.

    Article  PubMed  CAS  Google Scholar 

  41. Rabinowicz PD, Citek R, Budiman MA, Nunberg A, Bedell JA, Lakey N, O'Shaughnessy AL, Nascimento LU, McCombie WR, Martienssen RA: Differential methylation of genes and repeats in land plants. Genome Res. 2005, 15 (10): 1431-1440. 10.1101/gr.4100405.

    Article  PubMed  CAS  Google Scholar 

  42. Massa AN, Wanjugi H, Deal KR, O'Brien K, You FM, Maiti R, Chan AP, Gu YQ, Luo MC, Anderson OD, Rabinowicz PD, Dvorak J, Devos KM: Gene space dynamics during the evolution ofAegilops tauschii,Brachypodium distachyon,Oryza sativa, andSorghum bicolorGenomes. Mol Biol Evol. 2011, 28 (9): 2537-2547. 10.1093/molbev/msr080.

    Article  PubMed  CAS  Google Scholar 

  43. Erayman M, Sandhu D, Sidhu D, Dilbirligi M, Baenziger PS, Gill KS: Demarcating the gene-rich regions of the wheat genome. Nucleic Acids Res. 2004, 32 (12): 3546-3565. 10.1093/nar/gkh639.

    Article  PubMed  CAS  Google Scholar 

  44. Huo N, Lazo GR, Vogel JP, You FM, Ma Y, Hayde DM, Coleman-Derr D, Hill TA, Dvorak J, Anderson OD, et al: The nuclear genome ofBrachypodium distachyon: analysis of BAC end sequences. Funct Integr Genomics. 2008, 8 (2): 135-147. 10.1007/s10142-007-0062-7.

    Article  PubMed  CAS  Google Scholar 

  45. Paux E, Faure S, Choulet F, Roger D, Gauthier V, Martinant J-P, Sourdille P, Balfourier F, Le Paslier M-C, Chauveau A, Cakir M, Gandon B, Feuillet C: Insertion site-based polymorphism markers open new perspectives for genome saturation and marker-assisted selection in wheat. Plant Biotech J. 2010, 8 (2): 196-210. 10.1111/j.1467-7652.2009.00477.x.

    Article  CAS  Google Scholar 

  46. McIntosh RA, Yamazaki Y, Devos KM, Dubcovsky J, Rogers WJ, Appels R: Catalogue of gene symbols for wheat. Wheat Inf Serv. 2003, 97: 27-37.

    Google Scholar 

  47. Wicker T, Matthews DE, Keller B: TREP: A database for Triticeae repetitive elements. TrendsPlant Sci. 2002, 7 (12): 561-562.

    Article  CAS  Google Scholar 

  48. Jurka J, Kapitonov VV, Pavlicek A, Klonowski P, Kohany O, Walichiewicz J: Repbase update, a database of eukaryotic repetitive elements. Cytogenet Genome Res. 2005, 110 (1–4): 462-467.

    Article  PubMed  CAS  Google Scholar 

  49. Bao ZR, Eddy SR: Automated de novo identification of repeat sequence families in sequenced genomes. Genome Res. 2002, 12 (8): 1269-1276. 10.1101/gr.88502.

    Article  PubMed  CAS  Google Scholar 

  50. Thompson JD, Higgins DG, Gibson TJ: CLUSTAL-W - Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22 (22): 4673-4680. 10.1093/nar/22.22.4673.

    Article  PubMed  CAS  Google Scholar 

  51. Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA: Circos: An information aesthetic for comparative genomics. Genome Res. 2009, 19 (9): 1639-1645. 10.1101/gr.092759.109.

    Article  PubMed  CAS  Google Scholar 

Download references

Acknowledgements

This work was supported by United States Department of Agriculture (USDA-NIFA grant nos. 2006-35604-17248 and 2008-35300-04588). HS and JD were supported by the Czech Ministry of Education, Youth and Sports and the European Regional Development Fund (Operational Program Research and Development for Innovations No. ED0007/01/01). We thank Jayachandran Pillamari for customized scripts and developing figures. This is contribution no. 12-240-J from the Kansas Agricultural Experiment Station, Manhattan, KS 66506–5502, USA.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Bikram S Gill.

Electronic supplementary material

Additional file 1

: Sequence for putative novel repeats identified from 3AS BAC end sequences. (FASTA 13 KB)

Additional file 2

: Novel repeats, length and number of hits identified in 3AS BAC end sequences. (XLS 26 KB)

Additional file 3

: Blast2GO annotations of repeat masked BESs. (XLS 446 KB)

Additional file 4

: Class wise distribution of repetitive elements in homoeologous genomes of hexaploid wheat based on BESs analysis. (XLS 53 KB)

Additional file 5

: Comparative analyses of wheat chromosomes 3AS BESs and 3B unigenes. (XLS 39 KB)

Additional file 6

: Microsynteny between wheat chromosome 3A and rice chromosome 1, Brachypodium chromosome 2 and Sorghum chromosome 3. (XLS 26 KB)

Additional file 7

: SSR and ISBP markers developed from BESs. (XLS 352 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Sehgal, S.K., Li, W., Rabinowicz, P.D. et al. Chromosome arm-specific BAC end sequences permit comparative analysis of homoeologous chromosomes and genomes of polyploid wheat. BMC Plant Biol 12, 64 (2012). https://doi.org/10.1186/1471-2229-12-64

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2229-12-64

Keywords