Full-Length Enriched cDNA Libraries and ORFeome Analysis of Sugarcane Hybrid and Ancestor Genotypes

<div><p>Sugarcane is a major crop used for food and bioenergy production. Modern cultivars are hybrids derived from crosses between <i>Saccharum officinarum</i> and <i>Saccharum spontaneum</i>. Hybrid cultivars combine favorable characteristics from ancestral species and contain a genome that is highly polyploid and aneuploid, containing 100–130 chromosomes. These complex genomes represent a huge challenge for molecular studies and for the development of biotechnological tools that can facilitate sugarcane improvement. Here, we describe full-length enriched cDNA libraries for <i>Saccharum officinarum</i>, <i>Saccharum spontaneum</i>, and one hybrid genotype (SP803280) and analyze the set of open reading frames (ORFs) in their genomes (i.e., their ORFeomes). We found 38,195 (19%) sugarcane-specific transcripts that did not match transcripts from other databases. Less than 1.6% of all transcripts were ancestor-specific (i.e., not expressed in SP803280). We also found 78,008 putative new sugarcane transcripts that were absent in the largest sugarcane expressed sequence tag database (SUCEST). Functional annotation showed a high frequency of protein kinases and stress-related proteins. We also detected natural antisense transcript expression, which mapped to 94% of all plant KEGG pathways; however, each genotype showed different pathways enriched in antisense transcripts. Our data appeared to cover 53.2% (17,563 genes) and 46.8% (937 transcription factors) of all sugarcane full-length genes and transcription factors, respectively. This work represents a significant advancement in defining the sugarcane ORFeome and will be useful for protein characterization, single nucleotide polymorphism and splicing variant identification, evolutionary and comparative studies, and sugarcane genome assembly and annotation.</p></div

Description of cloned full-length cDNA libraries.

a<p>Excluding low-quality sequences and repeated clones.</p>b<p>Excluding sugarcane-specific transcripts.</p><p>Description of cloned full-length cDNA libraries.</p

Data from Ion PGM and 454 sequencing (after trimming) and Trinity assembly output.

<p>Data from Ion PGM and 454 sequencing (after trimming) and Trinity assembly output.</p

Functional annotation of sugarcane full-length cDNA contigs.

a<p>For some databases, several contigs showed more than one annotation or categorization.</p>b<p>Based on annotation from <i>Sorghum bicolor, Z. mays</i>, <i>P. virgatum</i>, <i>Setaria italica</i>, <i>O. sativa</i>, and <i>B. distachyon</i>.</p><p>Functional annotation of sugarcane full-length cDNA contigs.</p

Number of grass transcripts mapping to sugarcane contigs.

<p>A, Percentage of total grass transcripts mapping to sugarcane contigs. B, Total grass transcripts mapping to sugarcane contigs (white bars) and total sugarcane contigs mapping to each grass database (black bars). C, Total sugarcane contigs mapping to grasses, Uniprot, and NR databases and total unmatched sugarcane contigs (putative sugarcane-specific transcripts).</p

Number of reads identified as natural antisense transcripts (NATs), based on their orientation on alignment to grass genes.

<p>Tissues are indicated as follows: In1, immature internodes; In5, intermediate internodes; L, leaves.</p

Top four categories of genotype-specific contigs based on Phytozome annotations.

<p>Only contigs from leaf samples in each genotype were considered.</p

Number of contigs and percentage of coverage by antisense reads.

<p>Non-overlapping contig coverage by antisense reads was calculated for all contigs showing antisense expression (28,844 contigs).</p

Full-length enrichment for library cloning and next generation sequencing (NGS).

<p>Full-length (blue line with 5′ cap) or truncated (short blue line without 5′ cap) mRNAs were reverse transcribed into first-strand cDNA using oligo-dT primers (red arrow). The mRNA:cDNA hybrid was treated with RNase I (scissor) to remove the single-stranded RNA that was not fully extended by the first-strand cDNA, followed by selection for full-length transcripts using Cap-antibody magnetic beads to enrich the full-length mRNA:cDNA. The full-length single-stranded DNA (FLssDNA) was eluted from beads and used for both cDNA library cloning (lower left) and NGS (lower right). For full-length library cloning, a double-stranded adaptor (green) was linked to the 5′ end of ssDNA. Second-strand cDNA synthesis was then carried out, followed by cloning into a vector. For NGS, the full-length enriched ssDNA was fragmented by sonication to target fragments in the range of 200–400 bp, followed by ligation of the double-stranded DNA sequencing adaptor mixture (purple) to 3′ and 5′ ends of ssDNA. To maintain the complexity of the library while enriching the full-length cDNA for NGS, the original polyA mRNA was also fragmented using RNAse III, followed by ligation of the double-stranded RNA sequencing adaptor mixture (brown) to 3′ and 5′ ends of mRNA. After first- and second-strand synthesis, the polyA and capped mRNA and polyA and non-capped mRNA samples were mixed in a 3∶1 ratio and applied to the downstream NGS procedure.</p

Length distribution of sugarcane ORFs, full-length transcripts (FL), and UTRs.

<p>Graphs denote the comparison of sugarcane length distribution (gray bars) of ORFs (A) and full-length transcripts (B) to other grasses (colored lines). Length distribution of 5′ and 3′ UTRs (C, black and white bars, respectively) of sugarcane full-length transcripts is shown as well.</p