De novo assembly and characterizing of the culm-derived meta-transcriptome from the polyploid sugarcane genome based on coding transcripts

The subgenome Saccharum spontaneum contributes to sugar accumulation in sugarcane as revealed by full-length transcriptomic analysis

INTRODUCTION

Modern sugarcane cultivars (Saccharum spp. hybrids) derived from crosses between S. officinarum and S. spontaneum, with high-sugar traits and excellent stress tolerance inherited respectively. However, the contribution of the S. spontaneum subgenome to sucrose accumulation is still unclear.

OBJECTIVE

To compensate for the absence of a high-quality reference genome, a transcriptome analysis method is needed to analyze the molecular basis of differential sucrose accumulation in sugarcane hybrids and to find clues to the contribution of the S. spontaneum subgenome to sucrose accumulation.

METHODS

PacBio full-length sequencing was used to complement genome annotation, followed by the identification of differential genes between the high and low sugar groups using differential alternative splicing analysis and differential expression analysis. At the subgenomic level, the factors responsible for differential sucrose accumulation were investigated from the perspective of transcriptional and post-transcriptional regulation.

RESULTS

A full-length transcriptome annotated at the subgenomic level was provided, complemented by 263,378 allele-defined transcript isoforms and 139,405 alternative splicing (AS) events. Differential alternative splicing (DA) analysis and differential expression (DE) analysis identified differential genes between high and low sugar groups and explained differential sucrose accumulation factors by the KEGG pathways. In some gene models, different or even opposite expression patterns of alleles from the same gene were observed, reflecting the potential evolution of these alleles toward novel functions in polyploid sugarcane. Among DA and DE genes in the sucrose source-sink complex pathway, we found some alleles encoding sucrose accumulation-related enzymes derived from the S. spontaneum subgenome were differentially expressed or had DA events between the two contrasting sugarcane hybrids.

CONCLUSION

Full-length transcriptomes annotated at the subgenomic level could better characterize sugarcane hybrids, and the S. spontaneum subgenome was found to contribute to sucrose accumulation.

Comparative de novo transcriptome analysis identifies salinity stress responsive genes and metabolic pathways in sugarcane and its wild relative Erianthus arundinaceus [Retzius] Jeswiet

Erianthus arundinaceus [Retzius] Jeswiet, a wild relative of sugarcane has a high biomass production potential and a reservoir of many genes for superior agronomic traits and tolerance to biotic and abiotic stresses. A comparative physiological, anatomical and root transcriptome analysis were carried out to identify the salt-responsive genes and metabolic pathways associated with salt-tolerant E. arundinaceus genotype IND99-907 and salinity-sensitive sugarcane genotype Co 97010. IND99-907 recorded growth of young leaves, higher proline content, higher relative water content, intact root anatomical structures and lower Na+/K+, Ca2+/K+ and Mg2+/K+ ratio as compared to the sugarcane genotype Co 97010. We have generated four de novo transcriptome assemblies between stressed and control root samples of IND99-907 and Co 97010. A total of 649 and 501 differentially expressed genes (FDR<0.01) were identified from the stressed and control libraries of IND99-907 and Co 97010 respectively. Genes and pathways related to early stress-responsive signal transduction, hormone signalling, cytoskeleton organization, cellular membrane stabilization, plasma membrane-bound calcium and proton transport, sodium extrusion, secondary metabolite biosynthesis, cellular transporters related to plasma membrane-bound trafficking, nucleobase transporter, clathrin-mediated endocytosis were highly enriched in IND99-907. Whereas in Co 97010, genes related to late stress-responsive signal transduction, electron transport system, senescence, protein degradation and programmed cell death, transport-related genes associated with cellular respiration and mitochondrial respiratory chain, vesicular trafficking, nitrate transporter and fewer secondary metabolite biosynthetic genes were highly enriched. A total of 27 pathways, 24 biological processes, three molecular functions and one cellular component were significantly enriched (FDR≤ 0.05) in IND99-907 as compared to 20 pathways, two biological processes without any significant molecular function and cellular components in Co 97010, indicates the unique and distinct expression pattern of genes and metabolic pathways in both genotypes. The genomic resources developed from this study is useful for sugarcane crop improvement through development of genic SSR markers and genetic engineering approaches.

Pathways of Photosynthesis in Non-Leaf Tissues

Plants have leaves as specialised organs that capture light energy by photosynthesis. However, photosynthesis is also found in other plant organs. Photosynthesis may be found in the petiole, stems, flowers, fruits, and seeds. All photosynthesis can contribute to the capture of carbon and growth of the plant. The benefit to the plant of photosynthesis in these other tissues or organs may often be associated with the need to re-capture carbon especially in storage organs that have high respiration rates. Some plants that conduct C3 photosynthesis in the leaves have been reported to use C4 photosynthesis in petioles, stems, flowers, fruits, or seeds. These pathways of non-leaf photosynthesis may be especially important in supporting plant growth under stress and may be a key contributor to plant growth and survival. Pathways of photosynthesis have directionally evolved many times in different plant lineages in response to environmental selection and may also have differentiated in specific parts of the plant. This consideration may be useful in the breeding of crop plants with enhanced performance in response to climate change.

Comparative Transcriptome Profiling of Resistant and Susceptible Sugarcane Cultivars in Response to Infection by Xanthomonas albilineans

Sugarcane (Saccharum spp. hybrids) is a major source of sugar and renewable bioenergy crop worldwide and suffers serious yield losses due to many pathogen infections. Leaf scald caused by Xanthomonas albilineans is a major bacterial disease of sugarcane in most sugarcane-planting countries. The molecular mechanisms of resistance to leaf scald in this plant are, however, still unclear. We performed a comparative transcriptome analysis between resistant (LCP 85-384) and susceptible (ROC20) sugarcane cultivars infected by X. albilineans using the RNA-seq platform. 24 cDNA libraries were generated with RNA isolated at four time points (0, 24, 48, and 72 h post inoculation) from the two cultivars with three biological replicates. A total of 105,783 differentially expressed genes (DEGs) were identified in both cultivars and the most upregulated and downregulated DEGs were annotated for the processes of the metabolic and single-organism categories, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of the 7612 DEGs showed that plant-pathogen interaction, spliceosome, glutathione metabolism, protein processing in endoplasmic reticulum, and plant hormone signal transduction contributed to sugarcane's response to X. albilineans infection. Subsequently, relative expression levels of ten DEGs determined by quantitative reverse transcription-PCR (qRT-PCR), in addition to RNA-Seq data, indicated that different plant hormone (auxin and ethylene) signal transduction pathways play essential roles in sugarcane infected by X. albilineans. In conclusion, our results provide, for the first time, valuable information regarding the transcriptome changes in sugarcane in response to infection by X. albilineans, which contribute to the understanding of the molecular mechanisms underlying the interactions between sugarcane and this pathogen and provide important clues for further characterization of leaf scald resistance in sugarcane.

Genomic resources for energy cane breeding in the post genomics era

Sugarcane is one of the most sustainable energy crops among cultivated crops presenting the highest tonnage of cultivated plants. Its high productivity of sugar, bioethanol and bioelectricity make it a promising green alternative to petroleum. Furthermore, the myriad of products that can be derived from sugarcane biomass has been driving breeding programs towards varieties with a higher yield of fiber and a more vigorous and sustainable performance: the energy cane. Here we provide an overview of the energy cane including plant description, breeding efforts, types, and end-uses. In addition, we describe recently published genomic resources for the development of this crop, discuss current knowledge of cell wall metabolism, bioinformatic tools and databases available for the community.

The Impact of cDNA Normalization on Long-Read Sequencing of a Complex Transcriptome

Normalization of cDNA is widely used to improve the coverage of rare transcripts in analysis of transcriptomes employing next-generation sequencing. Recently, long-read technology has been emerging as a powerful tool for sequencing and construction of transcriptomes, especially for complex genomes containing highly similar transcripts and transcript-spliced isoforms. Here, we analyzed the transcriptome of sugarcane, a highly polyploidy plant genome, by PacBio isoform sequencing (Iso-Seq) of two different cDNA library preparations, with and without a normalization step. The results demonstrated that, while the two libraries included many of the same transcripts, many longer transcripts were removed, and many new generally shorter transcripts were detected by normalization. For the same input cDNA and data yield, the normalized library recovered more total transcript isoforms and number of predicted gene families and orthologous groups, resulting in a higher representation for the sugarcane transcriptome, compared to the non-normalized library. The non-normalized library, on the other hand, included a wider transcript length range with more longer transcripts above ∼1.25 kb and more transcript isoforms per gene family and gene ontology terms per transcript. A large proportion of the unique transcripts comprising ∼52% of the normalized library were expressed at a lower level than the unique transcripts from the non-normalized library, across three tissue types tested including leaf, stalk, and root. About 83% of the total 5,348 predicted long noncoding transcripts was derived from the normalized library, of which ∼80% was derived from the lowly expressed fraction. Functional annotation of the unique transcripts suggested that each library enriched different functional transcript fractions. This demonstrated the complementation of the two approaches in obtaining a complete transcriptome of a complex genome at the sequencing depth used in this study.