Microcollinearity between autopolyploid sugarcane and diploid sorghum genomes
© Wang et al; licensee BioMed Central Ltd. 2010
Received: 2 December 2009
Accepted: 23 April 2010
Published: 23 April 2010
Sugarcane (Saccharum spp.) has become an increasingly important crop for its leading role in biofuel production. The high sugar content species S. officinarum is an octoploid without known diploid or tetraploid progenitors. Commercial sugarcane cultivars are hybrids between S. officinarum and wild species S. spontaneum with ploidy at ~12×. The complex autopolyploid sugarcane genome has not been characterized at the DNA sequence level.
The microsynteny between sugarcane and sorghum was assessed by comparing 454 pyrosequences of 20 sugarcane bacterial artificial chromosomes (BACs) with sorghum sequences. These 20 BACs were selected by hybridization of 1961 single copy sorghum overgo probes to the sugarcane BAC library with one sugarcane BAC corresponding to each of the 20 sorghum chromosome arms. The genic regions of the sugarcane BACs shared an average of 95.2% sequence identity with sorghum, and the sorghum genome was used as a template to order sequence contigs covering 78.2% of the 20 BAC sequences. About 53.1% of the sugarcane BAC sequences are aligned with sorghum sequence. The unaligned regions contain non-coding and repetitive sequences. Within the aligned sequences, 209 genes were annotated in sugarcane and 202 in sorghum. Seventeen genes appeared to be sugarcane-specific and all validated by sugarcane ESTs, while 12 appeared sorghum-specific but only one validated by sorghum ESTs. Twelve of the 17 sugarcane-specific genes have no match in the non-redundant protein database in GenBank, perhaps encoding proteins for sugarcane-specific processes. The sorghum orthologous regions appeared to have expanded relative to sugarcane, mostly by the increase of retrotransposons.
The sugarcane and sorghum genomes are mostly collinear in the genic regions, and the sorghum genome can be used as a template for assembling much of the genic DNA of the autopolyploid sugarcane genome. The comparable gene density between sugarcane BACs and corresponding sorghum sequences defied the notion that polyploidy species might have faster pace of gene loss due to the redundancy of multiple alleles at each locus.
Sugarcane (Saccharum spp. L., Poaceae) is a large, perennial, tropical or subtropical grass widely grown primarily for sugar production worldwide. It is a first generation biofuel crop used for ethanol and biomass production as an alternative source of energy . About 75% of the world's sugar (sucrose) supply is from sugarcane and the other 25% from sugar beet (Beta vulgaris L., Chenopodiaceae). As a C4 plant, sugarcane has been recognized as one of the world's most efficient crops in converting solar energy into chemical energy, specialized for sucrose production [2–4]. Sugarcane is also among the crops having the most favorable energy input/output ratio [5, 6].
The genus Saccharum includes six species based on morphology, chromosome numbers, and geographical distribution, and they are S. spontaneum, S. robustum, S. officinarum, S. barberi, S. sinense, and S. edule. Recent genomic and molecular cytogenetic data provided strong evidence that S. barberi and S. sinense were derived from interspecific hybridization between S. officinarum and S. spontaneum. Although it has not been proven, S. edule is thought to be an intergeneric hybrid between either S. officinarum or S. robustum with a related genus that might account for its aborted inflorescence . S. robustum and S. spontaneum are two wild species with different basic chromosomes, x = 10 and x = 8, respectively [9–12]. These two wild species have a wide range of chromosome numbers and ploidy levels with 2n = 60 - 170 for S. robustum and 2n = 40 - 128 for S. spontaneum. S. robustum has been postulated to be the progenitor of the high sugar content species, S. officinarum (2n = 8× = 80). The unique basic chromosome number and distinctive DNA fingerprints of S. spontaneum from the other species of Saccharum are the reasons for a proposal to divide this genus to only two species, S. spontaneum as traditionally defined and S. officinarum including all other species and interspecific hybrids .
Hybrids derived from crosses involving a female S. officinarum (2n = 80) and a male S. spontaneum exhibit 2n + n transmission, conserving the entire genome of S. officinarum, a phenomena known as female restitution . The female restitution remains true in the first backcross between female S. officinarum and the 2n + n F1, but breaks down at the subsequent backcross. Dutch breeders in early 1900 utilized this unusual phenomenon in sugarcane improvement to integrate resistance genes for biotic and abiotic stresses from the wild species S. spontaneum and quickly recover the high sugar content property by a few backcrosses to the high sugar content species S. officinarum. For that reason, all current sugarcane cultivars in production are hybrids with 80-90% of the genome from S. officinarum and 10-20% of the genome from S. spontaneum[15, 16].
The complexity of the autopolyploid genome and the interspecific hybridization of modern cultivars hindered progress in genetic/genomic research and the application of genomic tools in sugarcane breeding programs. The only sugarcane bacterial artificial chromosome (BAC) library that we are aware of was constructed from the hybrid cultivar R570 with 2n = 115 chromosomes. Given an estimated genome size of 10 Gb, this BAC library provided 1.3× coverage of the polyploid genome and 14× coverage of the basic chromosome set [17, 18]. Sugarcane cultivars used for genetic mapping often have more than 100 chromosomes, and all sugarcane genetic maps constructed to date appear to be incomplete, due to the large number of chromosomes to be mapped and the limited genomic sequences available for developing markers [16, 19–22]. This deficiency has restricted the application of marker-assisted selection, because much of the genome cannot yet be scanned for target traits. However, sugarcane is in the forefront as a source of biofuel, and this has stimulated investment from both private and public sectors in sugarcane research. Coupled with the decreased cost of DNA sequencing using the next generation sequencing technologies, the once daunting and prohibitive task of sequencing the autopolyploid sugarcane genome becomes a real possibility.
Sorghum (Sorghum bicolor (L.) Moench) is a major cereal crop that provides food, feed, fiber, and fuel. It is domesticated in northern Africa with the ability to be productive in adverse environment and is valued for its drought tolerance. In comparison to polyploid sugarcane, sorghum is a diploid with 10 chromosomes and has a genome of about 730 Mb . The size of the monoploid sugarcane genome is estimated to be of similar magnitude, approximately 930 Mb; the genome complexity in sugarcane comes from the ploidy level and the two genome sets presented in sugarcane hybrids. The sorghum genome has been sequenced because of its small genome size and its importance for food security and biofuel production in diverse environments, particularly developing countries in the tropics . Sorghum and sugarcane belong to the same subtribe, Saccharinae, within the grass family Poaceae  and are close relatives to each other, sharing a common ancestor about 8-9 million years ago . The completion of the sorghum genome sequence offered unprecedented opportunities for sugarcane genomic research .
The synteny between sugarcane and sorghum genome has been reported before using DNA markers , but no details of microsynteny are available except one pair of sugarcane bacterial artificial chromosomes (BACs) containing the Adh1 gene . The objectives of this project are: (1) to test the suitability of using the sorghum genome as a template to assist assembly of sugarcane sequences generated from 20 selected sugarcane BACs using 454 Flex; (2) to explore the sequence features of the autopolyploid by examining a large set of long contiguous genomic sequences; and (3) to test the hypothesis of microcollinearity between sorghum and sugarcane at 20 chromosomal locations across the respective genomes.
Selection and sequencing of 20 sugarcane BACs
Summary of sequencing results of the 20 sugarcane BACs using 454 Flex and the 2 BACs using Sanger method.
General BAC clone information
BAC sequence assembly
Sugarcane - sorghum sequence alignment
Insert size (kb)
No. of contigs
Assembled BAC (bp)
No. of contig >10 k
No. of un-ambigious contigsb
No. of ordered contigs
Total length of ordered contigs (kb)
Spanned in sugarcane (bp)
Spanned in sorghum (bp)
Initial assembly of BAC sequences based on the row and column pools resulted in a total of 287 contigs, ranging from 5 to 42 contigs per BAC with a combined length of 2.61 Mb (Table 1). These assembled contigs were examined to eliminate redundant small assemblies, and 218 unambiguous contigs with a total of 1.99 Mb sequences were sorted out to represent the sequence of 19 BACs with an estimated sum of 1.86 Mb (Table 1). To order the multiple contigs for each BAC, orthologous sorghum sequences were used as templates to orient these sugarcane contigs and fill the gaps between the contigs when possible. A sum of 1.45 Mb sugarcane contigs were unambiguously ordered, accounting for 78.2% of the 1.86 Mb BAC sequences (Table 1). Contigs of two BACs, SC0172L01 and SC0175O07, were not ordered because all the contigs were less than 9 kb with a few genes predicted but projected to different sorghum super contigs. Among the 1.45 Mb ordered sugarcane contig sequences, the sequence aligned with sorghum sequence spanned 0.99 Mb (Table 1), accounting for about 53.1% of the sugarcane BAC sequences.
Gene content and repetitive sequences
To annotate the sequenced sugarcane BACs, the repetitive sequences of the assembled BAC contigs were first masked by RepeatMasker using a repeat database combined from RepBase databases and TIGR Plant Repeat Databases. The remaining sequences were annotated using sugarcane ESTs, sorghum ESTs, and predicted sorghum gene models. In total, 209 protein coding genes were predicted in the 19 sugarcane BAC sequences, including 155 (74.2%) validated by sugarcane ESTs, 28 by sorghum ESTs, and additional 26 corresponding to sorghum predicted gene models (Additional File 3). In the corresponding sorghum regions, 202 genes were annotated, including 122 (60.4%) validated by sorghum ESTs, 31 by sugarcane ESTs, and additional 49 by sorghum gene models. A total of 171 (81.8%) genes from the 19 BACs are predicted to be true orthologs between sugarcane and sorghum. Sequence identity between these two species across 140 kb coding region revealed a 95.2% exon sequence identity with a range of 80.1 to 100%.
Among these annotated genes, 19 sugarcane genes have no counterpart on the corresponding sorghum region while 12 sorghum genes have no counterpart on the corresponding sugarcane regions. The 19 putative sugarcane-specific genes are distributed on 11 BACs and are all supported by sugarcane ESTs. Among them, 17 have no orthologs in the sorghum genome and 2 have orthologs in other parts of the sorghum genome rather than in the corresponding regions. For the 12 sorghum genes without sugarcane orthologs, one is supported by sorghum EST, and the other 11 are predicted by sorghum gene models. The 17 sugarcane-specific genes were blasted against the non-redundant protein database in GenBank. Only five of them presented significant match (e value < E-06), two with known function, encoding beta-galactosidase 6 and auxin efflux carrier, respectively, and the other three encoding hypothetical proteins. The remaining 12 ESTs have no match in the non-redundant protein database.
Summary of repetitive sequences in sugarcane BACs and the orthologous euchromatic regions in sorghum.
Sugarcane BAC sequences
Sorghum homologous region
Percentage of the sequence (%)
Percentage of the sequence (%)
In the approximately 2 Mb sugarcane sequences, about 5 kb were simple repeats. We have designed 44 pairs of primer from the sequences flanking the simple sequence repeats (SSRs) for genetics mapping. Thirty six SSRs were amplified successfully and six were polymorphic between S. officinarum LA Purple and S. robustum Molokai 5829, the parents of our sugarcane mapping population (Additional File 6).
Comparative analysis of sugarcane and sorghum homologous sequences
Chromosomal rearrangements in the aligned regions between sugarcane and sorghum genomic sequences.
Total length (bp)
Number of events
Average length (bp) per event
Duplication in sugarcane
Duplication in sorghum
Tandem gene duplication was observed in both sugarcane and sorghum genomes. These genes were identified by aligning genome sequence with sugarcane and sorghum ESTs. For example, a gene duplication was found on sorghum chromosome 7 but not in the corresponding region of sugarcane BAC 118L15 contig 12. Three tandem copies of the gene predicted to encode the 60S ribosomal L10A protein is found in sorghum and only one copy in sugarcane. Another gene encoding serine carboxypeptidase 2 was duplicated in sugarcane BAC 109D21 contig 11-13-14 but not in the corresponding region of sorghum chromosome 2. There were two copies of this gene in sugarcane and one in sorghum . In another case, a gene encoding for receptor kinase was duplicated in sugarcane BAC 108J08 but not in the corresponding region of sorghum chromosome 9 (Additional File 8).
Estimated divergence time among sugarcane, sorghum, and rice.
Validation of 454 sequence assembly using Sanger sequences
To assess the quality of 454 Flex sequence assembly, we selected two BACs SC172L01 and SC118L15 for Sanger sequencing. SC118L15 appeared to harbor a substantial amount of rearrangement between sugarcane and sorghum and SC172L01 had 19 relatively short contigs assembled from 454 reads with the longest contig of 7 kb and a few genes being annotated that scattered in different regions of the sorghum genome. The Sanger sequence reads of these two BACs, were assembled into 3 and 6 contigs, respectively, after multiple attempts to close the remaining gaps by primer walking, each showing a reduction from the 19 and 12 contigs of 454 sequences. The Sanger sequences of these two BACs also matched the estimated insert size of 80 kb for SC172L01 and 60 kb for SC118L15, while the assembled 454 contigs are 52.4 kb and 101.7 kb for these two BACs respectively (Table 1).
The previous 454 sequence of SC172L01 had 19 short contigs (each less than 8 kb) and six genes including four retrotransposase genes. Its alignment to sorghum sequence was puzzling because the six genes aligned to different regions of the sorghum genome. The more complete Sanger sequences of SC172L01 allowed us to align it with the sorghum genome, which revealed chromosomal rearrangements after the divergence of these two species. The BAC SC172L01 aligned largely to sorghum chromosome 10, but a portion aligned to chromosomes 8 and 9 (Figures 6c, d, and Additional File 9). Though the Sanger sequence of SC118L15 was less complete than its 454 sequence, the assembly from Sanger reads further confirmed the local chromosome rearrangements harbored by this BAC (Figure 6e).
Because of the large number of chromosomes (often >100) and the nature of autopolyploids, both high density genetic mapping and physical mapping have proven to be challenging tasks in sugarcane. Currently, there is no physical map and no saturated genetic map that covers all chromosomes. Alternative approaches would need to be tested for a potential sugarcane genome sequencing project. Our results showed that the sorghum genome is an excellent template for assembly of sugarcane euchromatic sequences. The initial assembly of pooled 454 BAC sequences showed 40% inflation compared with estimated insert sizes, which likely was caused by multiple assemblies of repetitive sequences. After aligning the sequences with the sorghum genome using orthologous genes as anchors, 78.2% of the sugarcane BAC contigs could be ordered unambiguously and 53.1% of the sugarcane BAC sequences aligned with the sorghum genic regions. Sequences that were not aligned consisted of repetitive and non-coding sequences.
The suitability of the sorghum genome as a template for sugarcane genomic sequence assembly, at least for the genic regions, will be critical for strategic planning to sequence the sugarcane genome. Current BAC by BAC or whole genome shotgun sequencing approach would require a high density genetic map ideally with a density at two markers per Mb and a physical map with a10× genome coverage. For sugarcane, the only BAC library available is constructed from commercial hybrid cultivar R570 with 1.3× genome coverage . A 10× coverage BAC library would require one million clones with an average insert size of 100 kb, an expensive and laborious task. A high density genetic map would require mapping 20,000 markers for the 115 chromosomes of R570, and these markers would have to be sequence tagged to be useful for sequence assembly, not anonymous markers such as amplified fragment length polymorphism (AFLP) markers. In the past 20 years, 13 sugarcane maps have been constructed, and each of them covers only a fraction of the genome with less than 2,000 markers, and the majority of the markers in recent maps are AFLP markers [16, 19–22]. Fortunately, the cost of sequencing is declining rapidly with increased throughput. Most likely, a draft of the sugarcane genome will be generated before an ultra high density (2 markers per Mb) genetic map and a physical map are available, using the sorghum genome as a template for sequence assembly.
The sugarcane genome has gone through at least two rounds of genome wide duplication events to become an octoploid since its divergence from a common ancestor shared with sorghum. The two rounds of duplications might have occurred after the speciation event separated the two wild species S. robustum (x = 10) and S. spontaneum (x = 8) since these two species has different basic chromosome number [9–12], within 2 million years . Although each octoploid has eight genomes, it is not possible to distinguish each individual genome and every genome is a mosaic of all eight genome segments, because every chromosome is free to pair and recombine with any one of the other seven homologous chromosomes during meiosis, although it should be noted that most genetic maps of sugarcane showed some evidence of preferential pairing [19, 27]. For this reason, it might not be possible to distinguish the two recent genome wide duplications, and a minimum tiling path of BAC clones would be as a good representative as any one single genome in the octoploid. The hybrid cultivar R570 has 2n = 115 chromosomes with potentially 12 genomes. We selected a single BAC from each of 20 euchromatic regions corresponding to 20 distinctive chromosome arms (ended up with 18 arms due to the empty clone and a misplaced BAC), representing one of the potential 12 genome. We found more genes in sugarcane sequenced fragments than in sorghum in the aligned homologous regions (209 vs. 202), and more putative sugarcane specific genes (17) than sorghum specific genes (12). Two of the 19 initially annotated sugarcane specific genes have orthologs in other part of the sorghum genome, which left 17 to be most likely sugarcane specific genes. All 17 putative sugarcane-specific genes were validated by sugarcane ESTs, while only one of the 12 putative sorghum-specific genes was validated by sorghum ESTs. Moreover, 12 of the 17 sugarcane specific genes have no match in the non-redundant protein database in GenBank, suggesting that they are likely involved in sugarcane-specific processes. Although we masked the repetitive sequences of the BACs using plant repeat database, it is possible that some of them could be low copy transposable elements since we don't have a sugarcane specific repeat database.
The sugarcane EST project (SUCEST) yielded a database containing 237,954 ESTs assembled into 33,620 unigenes from 26 different cDNA libraries . This EST database validated 74.2% of the 209 annotated genes on the 19 sugarcane BACs, while only 60.4% of the 202 sorghum annotated genes were validated by sorghum ESTs. It might be a general rule that the EST databases of polyploid organisms represent higher percentage of genes than their diploid counterparts, because the multiple (12 in the case of sugarcane hybrids) allelic forms of each gene would result in greater chance of a particular allelic form to be sequenced in a collection of a wide range of tissues and developmental stages. However, more alleles don't necessarily increase the chance of a particular gene to be expressed in any type of tissues or developmental stages, as we have discovered two developmental stage specific genes in our RT PCR experiment involving 47 predicted genes.
The subtribe Saccharinae includes three major biofuel crops, sugarcane, Miscanthus, and sorghum. Sugarcane and Miscanthus are closely related and belong to the Saccharum complex . Sorghum is their closest relative outside of the Saccharum complex. Our estimate of a common ancestor shared by sugarcane and sorghum about 7.7 million years ago is in line with the 8-9 million years estimated by Jannoo et al . This time frame should be also applied to Miscanthus as it is a member of the Saccharum complex.
Most of the BAC sequences aligned with the sorghum sequences collinearly. However, one of the BAC (172L01) aligned to multiple chromosomes of sorghum, indicating large scale chromosomal rearrangements between sugarcane and sorghum genomes. Numerous local small scale (within a BAC) rearrangements between sugarcane and sorghum genomes were also detected. These sequence arrangements at both intra- and inter-chromosomal scales between the two species reflect their evolutionary history after their divergence about 8 million years ago. Our sugarcane BAC sequences provide the view of a representative genome of the possible 12 genomes in R570. It would be more interesting to document the rearrangements among sugarcane homologs, which should be far fewer.
The 2C genome size of R570 is about 10 Gb with an average of 87 Mb per chromosome among its 115 chromosomes, larger than the ~73 Mb per chromosome in sorghum . However, our data suggest that the sorghum sequences appear to be expanded compare to the sugarcane orthologous sequences studied, due to accumulation of retroelements, contradicting the genome size estimates from flow cytometry. If what we observed truly reflect the features of these two genomes, the basic genome of sugarcane (x = 10 or x = 8) could be smaller than that of sorghum. The discrepancy between the direct sequence comparison and the genome size estimates could be due to tendency of overestimating genome size by flow cytometry, as demonstrated by the sequenced genomes of rice and poplar [29, 30]. It is also possible that the discrepancy is caused by inaccurate assembly of repetitive sequences of the sugarcane BACs generated by 454 Flex. Finally, the small sampling of sugarcane BACs that we studied may not be representative of the genome as a whole.
Sugarcane has been cultivated and improved over thousands of years, beginning in prehistoric times with selection initially on natural variations and continuing with the modern techniques of hybridization and genetic engineering. Enormous yield increase has been achieved in the last century by breeding for yield, disease and insect resistance, and stress tolerance. While sugarcane farmers throughout the world face constant challenges to sustain profitability and protect the environment , breeders face not only those challenges but also a biological constraint as the gap between average farm yield and genetic yield potential is narrowed through improved agronomic practices . Sequencing the complex genome of autopolyploid sugarcane will provide the genomic resources to study genes and gene interactions controlling sugar yield, biomass yield, and other agronomic traits. A sugarcane genome sequence has the potential to revolutionize sugarcane improvement programs by providing high throughput genome wide screening for genomic selection , and for mining promoters of specific alleles.
Sugarcane is an economically important tropic crop primarily for sugar production but increasingly for biofeul production. Its large polyploid genome coupled with interspecific hybridization and aneuploid hindered sugarcane genomic research. The available genome sequence of sorghum, a closely relative of sugarcane, provides an exceptional opportunity to unravel the complex sugarcane genome. In this study, we strategically selected 20 sugarcane BACs each corresponding to a sorghum chromosome arm for sequencing to study the genome structure and organization. Sequence comparisons revealed that sugarcane genome is mostly collinear in the genic regions with sorghum genome, and the coding region of sugarcane and sorghum shared an average of 95.2% sequence identity. The unaligned regions between sugarcane and sorghum sequences were occupied mostly by repeats. The sorghum genome is an excellent template for assembling the genic DNA of the autopolyploid sugarcane genome. The comparable gene density between sugarcane and sorghum and the high number of sugarcane specific genes indicated that sugarcane genome might have retained more, not less, genes after the divergence of these two genera. Polyploidy species might not have faster pace of gene loss despite the redundancy of multiple alleles at each locus.
Selection of sugarcane BACs and sequencing
The detailed procedures for cloned, large insert genomic DNA isolation entailed a modified cleared lysate procedure as described in detail earlier [34, 35]. BAC DNA was prepared for sequencing on the 454/Roche GS-FLX as described by the manufacturer . Briefly this entailed shearing the purified BAC DNA via nebulization and subsequent end repair, as described , followed by ligation of adapter sequences and a second round of end repair to yield a blunt ended DNA library that then was quantified and pooled into 5 horizontal and 4 vertical pools (9 pools for 20 BACs total) using the Clone-Array Pooled Shotgun Sequencing strategy [38–40]. After dilution and emPCR amplification , the DNAs were loaded onto a 454/Roche GE-FLX for massively parallel pyrosequencing. The resulting sequence data was deconvoluted to individual BAC shotgun reads that then were assembled, first using the manufacturer supplied Newbler assembler and then by assembly with Phrap .
The GenBank accession numbers of these 19 BACs are: FJ348715-FJ348733. One of the 20 BACs has no insert.
BAC Sequence Annotation
The assembled BAC contig sequences were aligned to each other to identify and exclude ambiguous sequences such as sequence duplication, overlapping, and imbedding. The repeat sequences were masked in the un-ambiguous assembled sequence by RepeatMasker using a known repeats database combined from RepBase databases, TIGR Plant Repeat Databases. Genes were firstly annotated based on the spliced alignment of unmasked sequence assembly to sugarcane expression sequence tag (EST) from TIGR Plant Transcript Assemblies (Saccharum officinarum release 2). The reference gene set was further enriched by alignment of the unmasked assembly to closely related sorghum orthologues identified from sorghum EST and annotated CDs.
Two BAC clones, SC172L01 (60 kb) and SC11815 (80 kb), were sequenced using the shotgun approach with at least 10× coverage using Sanger sequencers. Approximately 3 ug BAC DNA was randomly sheared by a nebulizer (Invitrogen Corp.Carlsbad, CA USA) to produce fragments of 2-4 kb and then precipitated with Pellet Paint Co-Precipitant (Novagen, Darmstadt, Germany). The blunt-ends of DNA fragments were repaired using DNA terminator end repair kit (Lucigen, WI 53562 USA). The fragments with approximately 3 kb size were selected by cutting the gel slide, which was separated on a 1% agarose gel in 1 × TAE buffer through electrophoresis and purified using QIAquick Gel Extraction kit (Qiagen, Germantown, MD). The ligation and transformation were conducted by using Cells Clone Smart Blunt Cloning Kits (Lucigen, WI 53562 USA) following the manufacturer's instruction. DNAs from the 3-kb libraries were isolated through high through-put plasmid DNA miniprep and then subjected to cycle-sequence with ABI BigDye Terminator v3.1 and analyzed on a 3730XL DNA Analyzer (Applied Biosystems, Foster City, CA, USA). These two BACs were sequenced with 10× coverage.
Phred/Phrap/Consed and CAP3 packages were used for sequence assembly. Gaps in assembly and regions of low-quality were resolved by resequencing subclones identified by Autofinish, sequencing PCR products, and/or additional random subclone sequencing. All BAC clones were manually examined for signs of mis-assembly. Suspect regions were clarified either by ambiguous read removal, PCR amplification and sequencing, and/or alignment with a neighboring BAC. A BAC was not considered complete until all inconsistent read pairs were resolved and Consed reported an error rate of less than 1/10,000 bases. The GenBank accession numbers of the two BAC are: GU207345 and GU207346.
Comparative analysis of sugarcane and sorghum sequence
To identify the corresponding sorghum super contigs of the sugarcane BAC sequences, genes on sugarcane were aligned to the sorghum genome sequence in the web based BLAST search engine, http://www.phytozome.net/search.php?show=blast. The corresponding sorghum super contigs were used as anchor sequences to arrange the order (flip if necessary) of the sugarcane BAC sequence assemblies. The pairwise sequence comparisons and alignments between the arrayed sugarcane assembly and corresponding sorghum sequence were carried out on a visualized sequence alignment program WebACT . The sequence collinearity and local rearrangements were identified based on the above alignments with a bit value > 200.
Divergence time estimation between species
Orthologous gene pairs were identified. Protein sequences of orthologous gene pairs were aligned by CLUSTALW  and the protein alignments were converted back to DNA alignments using PAL2NAL . A few alignments were not reliable and discarded from further analysis. Ks (synonymous substitutions per synonymous site) values for these gene pairs were calculated using Nei-Gojobori method  implemented in PAML  package. The median Ks value was taken. The species divergence time were estimated by this formula: T_div = Ks/(2*6.5e-9). We used the commonly accepted synonymous substitution rate for grass lineage, estimated by Gaut et al. .
Closing gaps between BAC contigs
The gap sizes between the ordered adjacent contigs of each sugarcane BACs were estimated based on the corresponding sorghum gapless sequence. For the gap sizes less than 2 kb, primers were designed from two borders of the flanked contigs to amplify the fragment for gap filling by sequencing the PCR products. Primers were synthesized by Bioneer Inc (Alameda, CA 94501).
The PCR condition were as following: 10- μl PCR mix contained 1-ul of glycerol BAC stock as template DNA, 1× PCR buffer, 0.15 mM of each dNTP, 2.0 mM MgCl2, 0.15 μM of reverse and forward primers, and 0.8 units of Taq polymerase. The PCR reactions were performed using a 2720 Thermal Cycler (Applied Biosystems, USA), in which the reaction mixture was incubated at 94°C for 6 min, then for 35 cycles of 30 s of denaturing at 94°C, 40 s of annealing at 55°C, and 55 s of extension at 72°C, and then with a final extension at 72°C for 7 min. PCR products were separated on 1.5% agarose gels. The PCR products with single band were purified by using QIAquick PCR purification kit (QIAGEN Science, Maryland, USA) and sequenced at an ABI 3730XL core facility at Biotechnology Center of University of Illinois at Urbana- Champaign.
At least one intron was covered by primers designed for RT-PCR experiments to control genomic DNA contamination. Total RNA was extracted from two different tissues, mature leaf and leaf roll, and two genotypes, LA Purple and 95-4655. Approximately 2 μg of total RNA was treated with RNase-free DNase I (Promega, WI, USA) and reverse transcribed using RETROscript kit (Invitrogen, CA, USA). The synthesized cDNAs served as templates for RT-PCR. Four cDNA samples along with BAC DNA as positive control, RNA mix without DNase treatment and RNA mix after DNase treatment, were used as templates for PCR amplification. PCR products were run in 1% agarose gel.
We acknowledge our colleagues at the University of Oklahoma's Advanced Center for Genome Technology, Chunmei Qu and Ping Wang for their assistance with 454 GS-FLX sequencing sample preparation and Steve Kenton for his help with deconvoluting the pooled BACs and their subsequent assembly. We also thank Eric Tang for assistance on sequencing two BACs using Sanger sequencers. This project is supported by start-up funds from the University of Illinois to RM and a grant from the Energy Bioscience Institute (EBI) to SPM, MEH, RM, and DSR.
- Lam E, Shine J, da Silva J, Lawton M, Bonos S, Calvino M, Carrer H, Silva-Filho MC, Glynn N, Helsel Z, Ma J, Richard E, Souza G, Ming R: Improving Sugarcane for Biofuel: Engineering for an even better feedstock. Global Change Biology Bioenergy. 2009, 1: 251-255.View ArticleGoogle Scholar
- Aragon C, Carvalho LC, Gonzalez J, Escalona M, Amancio S: Sugarcane (Saccharum sp. Hybrid) propagated in headspace renovating systems shows autotrophic characteristics and develop improved anti-oxidative response. Tropical Plant Biology. 2009, 2: 38-50. 10.1007/s12042-008-9026-x.View ArticleGoogle Scholar
- De Maria Felix J, Papini-Terzi FS, Rocha FR, Vencio RZN, Vicentini R, Nishiyama MY, Cesar Ulian E, Souza GM, Menossi M: Expression profiling of signal transduction components in a sugarcane population segregating for sugar content. Tropical Plant Biology. 2009, 2: 98-109. 10.1007/s12042-009-9031-8.View ArticleGoogle Scholar
- Rae AL, Jackson MA, Nguyen CH, Bonnett GD: Functional specialization of vacuoles in sugarcane leaf and stem. Tropical Plant Biology. 2009, 2: 13-22. 10.1007/s12042-008-9019-9.View ArticleGoogle Scholar
- Heichel GH: Comparative efficiency of energy use in crop production. Bull Conn Ag Expt Sta. 1974, 739: 1-26.Google Scholar
- Tew TL: Genetic engeneering for greater energy efficiency in sugarcane. Conf Renewable Energy Technol. Dec 5. 1980, ; HonoluluGoogle Scholar
- D'Hont A, Paulet F, Glaszmann JC: Oligoclonal interspecific origin of 'North Indian' and 'Chinese' sugarcanes. Chromosome Res. 2002, 10 (3): 253-262. 10.1023/A:1015204424287.PubMedView ArticleGoogle Scholar
- Daniels J, Roach BT: Taxonomy and evolution. Sugarcane Improvement Through Breeding. Edited by: Heinz DJ. 1987, Amsterdam: Elsevier Press, 7-84.View ArticleGoogle Scholar
- D'Hont A, Grivet L, Feldmann P, Rao PS, Berding N, Glaszmann JC: Identification and characterization of sugarcane intergeneric hybrids, Saccharum officinarum x Erianthus arundinaceus, with molecular markers and DNA in situ hybridization. Theor Appl Genet. 1995, 91: 320-326. 10.1007/BF00220894.PubMedView ArticleGoogle Scholar
- D'Hont A, Grivet L, Feldmann P, Rao S, Berding N, Glaszmann JC: Characterisation of the double genome structure of modern sugarcane cultivars (Saccharum spp.) by molecular cytogenetics. Mol Gen Genet. 1996, 250 (4): 405-413. 10.1007/BF02174028.PubMedView ArticleGoogle Scholar
- D'Hont A, Ison D, Alix K, Roux C, Glaszmann JC: Determination of basic chromosome numbers in the genus Saccharum by physical mapping of ribosomal RNA genes. Genome. 1998, 41: 221-225. 10.1139/gen-41-2-221.View ArticleGoogle Scholar
- Ha S, Moore PH, Heinz D, Kato S, Ohmido N, Fukui K: Quantitative chromosome map of the polyploid Saccharum spontaneum by multicolor fluorescence in situ hybridization and imaging methods. Plant Mol Biol. 1999, 39 (6): 1165-1173. 10.1023/A:1006133804170.PubMedView ArticleGoogle Scholar
- Irvine JE: Saccharum species as horticultural classes. Theor Appl Genet. 1999, 98: 186-194. 10.1007/s001220051057.View ArticleGoogle Scholar
- Bremer G: Problems in breeding and cytology of sugarcane. Euphytica. 1961, 10: 59-78. 10.1007/BF00037206.View ArticleGoogle Scholar
- Grivet L, D'Hont A, Roques D, Feldmann P, Lanaud C, Glaszmann JC: RFLP Mapping in Cultivated Sugarcane (Saccharum spp.): Genome Organization in a Highly Polyploid and Aneuploid Interspecific Hybrid. Genetics. 1996, 142: 987-1000.PubMed CentralPubMedGoogle Scholar
- Hoarau JY, Grivet L, Offmann B, Raboin LM, Diorflar JP, Payet J, Hellmann M, D'Hont A, Glaszmann JC: Genetic dissection of a modern sugarcane cultivar (Saccharum spp.). II. Detection of QTLs for yield components. Theor Appl Genet. 2002, 105 (6-7): 1027-1037. 10.1007/s00122-002-1047-5.PubMedView ArticleGoogle Scholar
- Tomkins JP, Yu Y, Miller-Smith H, Frisch DA, Woo SS, Wing R: A bacterial artificial chromosome library for sugarcane. Theor Appl Genet. 1999, 99: 419-424. 10.1007/s001220051252.PubMedView ArticleGoogle Scholar
- D'Hont A: Unraveling the genome structure of polyploids using FISH and GISH; examples of sugarcane and banana. Cytogenet Genome Res. 2005, 109 (1-3): 27-33. 10.1159/000082378.PubMedView ArticleGoogle Scholar
- Ming R, Liu SC, Lin YR, da Silva J, Wilson W, Braga D, van Deynze A, Wenslaff TF, Wu KK, Moore PH, Burnquist W, Sorrells ME, Irvine JE, Paterson AH: Detailed alignment of saccharum and sorghum chromosomes: comparative organization of closely related diploid and polyploid genomes. Genetics. 1998, 150 (4): 1663-1682.PubMed CentralPubMedGoogle Scholar
- Ming R, Moore PH, Wu KK, D'Hont A, Glaszmann JC, Tew TL, Mirkov TE, da Silva J, Jifon J, Rai M, Schnell RJ, Brumbley SM, Lakshmanan P, Comstock JC, Paterson AH: Sugarcane improvement through breeding and biotechnology. Plant Breeding Review. 2005, 27: 15-118.Google Scholar
- Aitken KS, Jackson PA, McIntyre CL: A combination of AFLP and SSR markers provides extensive map coverage and identification of homo(eo)logous linkage groups in a sugarcane cultivar. Theor Appl Genet. 2005, 110 (5): 789-801. 10.1007/s00122-004-1813-7.PubMedView ArticleGoogle Scholar
- Garcia AA, Kido EA, Meza AN, Souza HM, Pinto LR, Pastina MM, Leite CS, Silva JA, Ulian EC, Figueira A, Souza AP: Development of an integrated genetic map of a sugarcane (Saccharum spp.) commercial cross, based on a maximum-likelihood approach for estimation of linkage and linkage phases. Theor Appl Genet. 2006, 112 (2): 298-314. 10.1007/s00122-005-0129-6.PubMedView ArticleGoogle Scholar
- Paterson AH, Bowers JE, Bruggmann R, Dubchak I, Grimwood J, Gundlach H, Haberer G, Hellsten U, Mitros T, Poliakov A, Schmutz J, Spannagl M, Tang H, Wang X, Wicker T, Bharti AK, Chapman J, Feltus FA, Gowik U, Grigoriev IV, Lyons E, Maher CA, Martis M, Narechania A, Otillar RP, Penning BW, Salamov AA, Wang Y, Zhang L, Carpita NC, Freeling M, Gingle AR, Hash CT, Keller B, Klein P, Kresovich S, McCann MC, Ming R, Peterson DG, Mehboob-ur-Rahman , Ware D, Westhoff P, Mayer KF, Messing J, Rokhsar DS: The Sorghum bicolor genome and the diversification of grasses. Nature. 2009, 457 (7229): 551-556. 10.1038/nature07723.PubMedView ArticleGoogle Scholar
- Clayton WD: Andropogoneae. Grass Systematics and Evolution. Edited by: Soderstrom TR, Hilu KW, Campbell CS, Barkworth ME. 1987, Smithsonian Institution Press, Washington, DC, 307-309.Google Scholar
- Jannoo N, Grivet L, Chantret N, Garsmeur O, Glaszmann JC, Arruda P, D'Hont A: Orthologous comparison in a gene-rich region among grasses reveals stability in the sugarcane polyploid genome. Plant J. 2007, 50 (4): 574-585. 10.1111/j.1365-313X.2007.03082.x.PubMedView ArticleGoogle Scholar
- Bowers JE, Arias MA, Asher R, Avise JA, Ball RT, Brewer GA, Buss RW, Chen AH, Edwards TM, Estill JC, Exum HE, Goff VH, Herrick KL, Steele CL, Karunakaran S, Lafayette GK, Lemke C, Marler BS, Masters SL, McMillan JM, Nelson LK, Newsome GA, Nwakanma CC, Odeh RN, Phelps CA, Rarick EA, Rogers CJ, Ryan SP, Slaughter KA, Soderlund CA, Tang H, Wing RA, Paterson AH: Comparative physical mapping links conservation of microsynteny to chromosome structure and recombination in grasses. Proc Natl Acad Sci USA. 2005, 102 (37): 13206-13211. 10.1073/pnas.0502365102.PubMed CentralPubMedView ArticleGoogle Scholar
- da Silva JA, Honeycutt RJ, Burnquist W, Al-Janabi SM, Sorrells ME, Tanksley SD, Sobral WS: Saccharum spontaneum L. 'SES208' genetic linkage map containing RFLP and PCR-based markers. Mol Breeding. 1995, 1: 165-179. 10.1007/BF01249701.View ArticleGoogle Scholar
- Vettore AL, da Silva FR, Kemper EL, Souza GM, da Silva AM, Ferro MI, Henrique-Silva F, Giglioti EA, Lemos MV, Coutinho LL, Nobrega MP, Carrer H, Franca SC, Bacci Junior M, Goldman MH, Gomes SL, Nunes LR, Camargo LE, Siqueira WJ, Van Sluys MA, Thiemann OH, Kuramae EE, Santelli RV, Marino CL, Targon ML, Ferro JA, Silveira HC, Marini DC, Lemos EG, Monteiro-Vitorello CB, Tambor JH, Carraro DM, Roberto PG, Martins VG, Goldman GH, de Oliveira RC, Truffi D, Colombo CA, Rossi M, de Araujo PG, Sculaccio SA, Angella A, Lima MM, de Rosa Junior VE, Siviero F, Coscrato VE, Machado MA, Grivet L, Di Mauro SM, Nobrega FG, Menck CF, Braga MD, Telles GP, Cara FA, Pedrosa G, Meidanis J, Arruda P: Analysis and functional annotation of an expressed sequence tag collection for tropical crop sugarcane. Genome Res. 2003, 13 (12): 2725-2735. 10.1101/gr.1532103.PubMed CentralPubMedView ArticleGoogle Scholar
- International Rice Genome Sequencing Project: The map-based sequence of the rice genome. Nature. 2005, 436 (7052): 793-800. 10.1038/nature03895.View ArticleGoogle Scholar
- Tuskan GA, Difazio S, Jansson S, Bohlmann J, Grigoriev I, Hellsten U, Putnam N, Ralph S, Rombauts S, Salamov A, Schein J, Sterck L, Aerts A, Bhalerao RR, Bhalerao RP, Blaudez D, Boerjan W, Brun A, Brunner A, Busov V, Campbell M, Carlson J, Chalot M, Chapman J, Chen GL, Cooper D, Coutinho PM, Couturier J, Covert S, Cronk Q, Cunningham R, Davis J, Degroeve S, Dejardin A, Depamphilis C, Detter J, Dirks B, Dubchak I, Duplessis S, Ehlting J, Ellis B, Gendler K, Goodstein D, Gribskov M, Grimwood J, Groover A, Gunter L, Hamberger B, Heinze B, Helariutta Y, Henrissat B, Holligan D, Holt R, Huang W, Islam-Faridi N, Jones S, Jones-Rhoades M, Jorgensen R, Joshi C, Kangasjarvi J, Karlsson J, Kelleher C, Kirkpatrick R, Kirst M, Kohler A, Kalluri U, Larimer F, Leebens-Mack J, Leple JC, Locascio P, Lou Y, Lucas S, Martin F, Montanini B, Napoli C, Nelson DR, Nelson C, Nieminen K, Nilsson O, Pereda V, Peter G, Philippe R, Pilate G, Poliakov A, Razumovskaya J, Richardson P, Rinaldi C, Ritland K, Rouze P, Ryaboy D, Schmutz J, Schrader J, Segerman B, Shin H, Siddiqui A, Sterky F, Terry A, Tsai CJ, Uberbacher E, Unneberg P, Vahala J, Wall K, Wessler S, Yang G, Yin T, Douglas C, Marra M, Sandberg G, Peer Van de Y, Rokhsar D: The genome of black cottonwood, Populus trichocarpa (Torr. & Gray). Science. 2006, 313 (5793): 1596-1604. 10.1126/science.1128691.PubMedView ArticleGoogle Scholar
- Glaz B: Integrated crop management for sustainable crop production: Recent advances. International Sugar Journal. 2003, 105: 175-186.Google Scholar
- Cassman KG: Ecological intensification of cereal production systems: yield potential, soil quality, and precision agriculture. Proc Natl Acad Sci USA. 1999, 96 (11): 5952-5959. 10.1073/pnas.96.11.5952.PubMed CentralPubMedView ArticleGoogle Scholar
- Heffner EL, Sorrells ME, Jannink JL: Genomic selection for crop improvement. Crop Sci. 2009, 49: 1-12. 10.2135/cropsci2008.08.0512.View ArticleGoogle Scholar
- Bodenteich A, Chissoe S, Wang YF, Roe BA: Shotgun Cloning as the Strategy of Choice to Generate Templates for High-throughput Dideoxynucleotide Sequencing. Automated DNA Sequencing and Analysis Techniques. Edited by: Venter JC. 1993, London: Academic Press, 42-50.Google Scholar
- Roe B, Crabtree J, Khan A: DNA Isolation and Sequencing: Essential Techniques Series. Edited by: Rickwood D. 1996, New York: J. Wiley and SonsGoogle Scholar
- Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA, Berka J, Braverman MS, Chen YJ, Chen Z, Dewell SB, Du L, Fierro JM, Gomes XV, Godwin BC, He W, Helgesen S, Ho CH, Irzyk GP, Jando SC, Alenquer ML, Jarvie TP, Jirage KB, Kim JB, Knight JR, Lanza JR, Leamon JH, Lefkowitz SM, Lei M, Li J, Lohman KL, Lu H, Makhijani VB, McDade KE, McKenna MP, Myers EW, Nickerson E, Nobile JR, Plant R, Puc BP, Ronan MT, Roth GT, Sarkis GJ, Simons JF, Simpson JW, Srinivasan M, Tartaro KR, Tomasz A, Vogt KA, Volkmer GA, Wang SH, Wang Y, Weiner MP, Yu P, Begley RF, Rothberg JM: Genome sequencing in microfabricated high-density picolitre reactors. Nature. 2005, 437 (7057): 376-380.PubMed CentralPubMedGoogle Scholar
- Roe BA: Shotgun library construction for DNA sequencing. Methods Mol Biol. 2004, 255: 171-187.PubMedGoogle Scholar
- Cai WW, Chen R, Gibbs RA, Bradley A: A clone-array pooled shotgun strategy for sequencing large genomes. Genome Res. 2001, 11 (10): 1619-1623. 10.1101/gr.198101.PubMedView ArticleGoogle Scholar
- Csuros M, Li B, Milosavljevic A: Clone-array pooled shotgun mapping and sequencing: design and analysis of experiments. Genome Inform. 2003, 14: 186-195.PubMedGoogle Scholar
- Csuros M, Milosavljevic A: Pooled Genomic Indexing (PGI): analysis and design of experiments. J Comput Biol. 2004, 11 (5): 1001-1021.PubMedView ArticleGoogle Scholar
- Green P: Phrap. 1999, [http://phrap.org]Google Scholar
- Abbott JC, Aanensen DM, Bentley SD: WebACT: an online genome comparison suite. Methods Mol Biol. 2007, 395: 57-74.PubMedView ArticleGoogle Scholar
- 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: 4673-4680. 10.1093/nar/22.22.4673.PubMed CentralPubMedView ArticleGoogle Scholar
- Suyama M, Torrents D, Bork P: PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments. Nucleic Acids Res. 2006, 34: 609-612. 10.1093/nar/gkl315.View ArticleGoogle Scholar
- Nei M, Gojobori T: Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. Mol Biol Evol. 1986, 3 (5): 418-26.PubMedGoogle Scholar
- Yang Z: PAML: a program package for phylogenetic analysis by maximum likelihood. Comput Appl Biosci. 1997, 13: 555-556.PubMedGoogle Scholar
- Gaut BS, Morton BR, McCaig BC, Clegg MT: Substitution rate comparisons between grasses and palms: Synonymous ratedifferences at the nuclear gene Adh parallel rate differences at the plastid gene rbcL. Proc Natl Acad Sci USA. 1996, 93: 10274-10279. 10.1073/pnas.93.19.10274.PubMed CentralPubMedView ArticleGoogle Scholar