De novo sequencing and transcriptome analysis of a low temperature tolerant Saccharum spontaneum clone IND 00-1037

Overexpression of EaALDH7, an aldehyde dehydrogenase gene from Erianthus arundinaceus enhances salinity tolerance in transgenic sugarcane (Saccharum spp. Hybrid)

Aldehyde Dehydrogenases (ALDH), a group of enzymes, are associated with the detoxification of aldehydes, produced in plants during abiotic stress conditions. Salinity remains a pivotal abiotic challenge that poses a significant threat to cultivation and yield of sugarcane. In this study, an Aldehyde dehydrogenase gene (EaALDH7) from Erianthus arundinaceus was overexpressed in the commercial sugarcane hybrid cultivar Co 86032. The transgenic lines were evaluated at different NaCl concentrations ranging from 0mM to 200mM for various morpho-physiological and biochemical parameters. The control plants, subjected to salinity stress condition, exhibited morphological changes in protoxylem, metaxylem, pericycle and pith whereas the transgenic events were on par with plants under regular irrigation. The overexpressing (OE) lines showed less cell membrane injury and improved photosynthetic rate, transpiration rate, and stomatal conductance than the untransformed control plants under stress conditions. Elevated proline content, higher activity of enzymatic antioxidants such as sodium dismutase (SOD), catalase (CAT), glutathione reductase (GR) and ascorbate peroxidase (APX) and low level of malondialdehyde MDA and hydrogen peroxide (H2O2) in the transgenic lines. The analysis of EaALDH7 expression revealed a significant upregulation in the transgenic lines compared to that of the untransformed control during salt stress conditions. The current study highlights the potentials of EaALDH7 gene in producing salinity-tolerant sugarcane cultivars.

Transgenic sugarcane overexpressing Glyoxalase III improved germination and biomass production at formative stage under salinity and water-deficit stress conditions

The glyoxalase system, involving Glyoxalase I (GlyI) and Glyoxalase II (Gly II), plays a vital role in abiotic stress tolerance in plants. A novel enzyme Glyoxalase III (Gly III) was found recently from bacteria, yeast, and plant species. This enzyme provides a new way to detoxify Methylglyoxal (MG), a cytotoxic α-oxoaldehyde, which, in excess, can cause complete cell destruction by forming Reactive Oxygen Species (ROS) and Advanced Glycation End products (AGEs) or DNA/RNA mutation. In this background, the current study examined sugarcane transgenic events that exhibit an increase in expression of EaGly III, to assess their performance in terms of germination and biomass production during formative stage under stress conditions. Southern blot analysis outcomes confirmed the integration of transgene in the transgenic plants. The results from quantitative RT-PCR analyses confirmed high expression levels of EaGly III in transgenic events compared to wild type (WT) under salinity (100 and 200 mM NaCl) and drought (withholding watering) conditions. Transgenic events exhibited enhanced biomass productivity ranged between 0.141 Kg/pot and 0.395 Kg/pot under 200 mM salinity and 0.262 Kg/pot and 0.666 Kg/pot under drought stress. Further, transgenic events observed significantly higher germination rates under salinity and drought conditions compared to that of WT. Subcellular localization prediction by EaGlyIII-GFP fusion expression in sugarcane callus showed that it is distributed across the cytoplasm, thus indicating its widespread activity within the cell. These results strongly suggest that enhancing EaGly III activity is a useful strategy to improve the salinity and drought-tolerance in sugarcane as well as other crops.

Comparative genome-wide characterization of salt responsive micro RNA and their targets through integrated small RNA and de novo transcriptome profiling in sugarcane and its wild relative Erianthus arundinaceus

Soil salinity and saline irrigation water are major constraints in sugarcane affecting the production of cane and sugar yield. To understand the salinity induced responses and to identify novel genomic resources, integrated de novo transcriptome and small RNA sequencing in sugarcane wild relative, Erianthus arundinaceus salt tolerant accession IND 99-907 and salt-sensitive sugarcane genotype Co 97010 were performed. A total of 362 known miRNAs belonging to 62 families and 353 miRNAs belonging to 63 families were abundant in IND 99-907 and Co 97010 respectively. The miRNA families such as miR156, miR160, miR166, miR167, miR169, miR171, miR395, miR399, miR437 and miR5568 were the most abundant with more than ten members in both genotypes. The differential expression analysis of miRNA reveals that 221 known miRNAs belonging to 48 families and 130 known miRNAs belonging to 42 families were differentially expressed in IND 99-907 and Co 97010 respectively. A total of 12,693 and 7982 miRNA targets against the monoploid mosaic genome and a total of 15,031 and 12,152 miRNA targets against the de novo transcriptome were identified for differentially expressed known miRNAs of IND 99-907 and Co 97010 respectively. The gene ontology (GO) enrichment analysis of the miRNA targets revealed that 24, 12 and 14 enriched GO terms (FDR < 0.05) for biological process, molecular function and cellular component respectively. These miRNAs have many targets that associated in regulation of biotic and abiotic stresses. Thus, the genomic resources generated through this study are useful for sugarcane crop improvement through biotechnological and advanced breeding approaches.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03867-7.

Isolation of 5' regulatory region of COLD1 gene and its functional characterization through transient expression analysis in tobacco and sugarcane

Chilling Tolerant Divergence 1 (COLD1) gene consists of Golgi pH Receptor (GPHR) as well as Abscisic Acid-linked G Protein-Coupled Receptor (ABA_GPCR), which are the major transmembrane proteins in plants. This gene expression has been found to be differentially regulated, under various stress conditions, in wild Saccharum-related genera, Erianthus arundinaceus, compared to commercial sugarcane variety. In this study, Rapid Amplification of Genomic Ends (RAGE) technique was employed to isolate the 5' upstream region of COLD1 gene to gain knowledge about the underlying stress regulatory mechanism. The current study established the cis-acting elements, main promoter regions, and Transcriptional Start Site (TSS) present within the isolated 5' upstream region (Cold1P) of COLD1, with the help of specific bioinformatics techniques. Phylogenetic analysis results revealed that the isolated Cold1P promoter is closely related to the species, Sorghum bicolor. Cold1P promoter-GUS gene construct was generated in pCAMBIA 1305.1 vector that displayed a constitutive expression of the GUS reporter gene in both monocot as well as dicot plants. The histochemical GUS assay outcomes confirmed that Cold1P can drive expression in both monocot as well as dicot plants. Cold1P's activities under several abiotic stresses such as cold, heat, salt, and drought, revealed its differential expression profile in commercial sugarcane variety. The highest activity of the GUS gene was found after 24 h of cold stress, driven by the isolated Cold1P promoter. The outcomes from GUS fluorimetric assay correlated with that of the GUS expression findings. This is the first report on Cold1P isolated from the species, E. arundinaceus.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03650-8.

Expression analysis of transcription factors in sugarcane during cold stress

Transcription factors (TF) are a wide class of genes in plants, and these can regulate the expression of other genes in response to various environmental stresses (biotic and abiotic). In the current study, transcription factor activity in sugarcane was examined during cold stress. Initially, RNA transcript reads of two sugarcane cultivars (ROC22 and GT08-1108) under cold stress were downloaded from SRA NCBI database. The reads were aligned into a reference genome and the differential expression analyses were performed with the R/Bioconductor edgeR package. Based on our analyses in the ROC22 cultivar, 963 TF genes were significantly upregulated under cold stress among a total of 5649 upregulated genes, while 293 TF genes were downregulated among a total of 3,289 downregulated genes. In the GT08-1108 cultivar, 974 TF genes were identified among 5,649 upregulated genes and 283 TF genes were found among 3,289 downregulated genes. Most transcription factors were annotated with GO categories related to protein binding, transcription factor binding, DNA-sequence-specific binding, transcription factor complex, transcription factor activity in RNA polymerase II, the activity of nucleic acid binding transcription factor, transcription corepressor activity, sequence-specific regulatory region, the activity of transcription factor of RNA polymerase II, transcription factor cofactor activity, transcription factor activity from plastid promoter, transcription factor activity from RNA polymerase I promoter, polymerase II and RNA polymerase III. The findings of above results will help to identify differentially expressed transcription factors during cold stress. It also provides a comprehensive analysis of the regulation of the transcription activity of many genes. Therefore, this study provides the molecular basis for improving cold tolerance in sugarcane and other economically important grasses.

PacBio Single-Molecule Long-Read Sequencing Reveals Genes Tolerating Manganese Stress in Schima superba Saplings

Schima superba (Theaceae) is a subtropical evergreen tree and is used widely for forest firebreaks and gardening. It is a plant that tolerates salt and typically accumulates elevated amounts of manganese in the leaves. With large ecological amplitude, this tree species grows quickly. Due to its substantial biomass, it has a great potential for soil remediation. To evaluate the thorough framework of the mRNA, we employed PacBio sequencing technology for the first time to generate S. Superba transcriptome. In this analysis, overall, 511,759 full length non-chimeric reads were acquired, and 163,834 high-quality full-length reads were obtained. Overall, 93,362 open reading frames were obtained, of which 78,255 were complete. In gene annotation analyses, the Kyoto Encyclopedia of Genes and Genomes (KEGG), Clusters of Orthologous Genes (COG), Gene Ontology (GO), and Non-Redundant (Nr) databases were allocated 91,082, 71,839, 38,914, and 38,376 transcripts, respectively. To identify long non-coding RNAs (lncRNAs), we utilized four computational methods associated with protein families (Pfam), Cooperative Data Classification (CPC), Coding Assessing Potential Tool (CPAT), and Coding Non-Coding Index (CNCI) databases and observed 8,551, 9,174, 20,720, and 18,669 lncRNAs, respectively. Moreover, nine genes were randomly selected for the expression analysis, which showed the highest expression of Gene 6 (Na_Ca_ex gene), and CAX (CAX-interacting protein 4) was higher in manganese (Mn)-treated group. This work provided significant number of full-length transcripts and refined the annotation of the reference genome, which will ease advanced genetic analyses of S. superba.

Single-Molecule Long-Read Sequencing of Purslane (Portulaca oleracea) and Differential Gene Expression Related with Biosynthesis of Unsaturated Fatty Acids

This study aimed to obtain the full-length transcriptome of purslane (Portulaca oleracea); assorted plant samples were used for single-molecule real-time (SMRT) sequencing. Based on SMRT, functional annotation of transcripts, transcript factors (TFs) analysis, simple sequence repeat analysis and long non-coding RNAs (LncRNAs) prediction were accomplished. Total 15.33-GB reads were produced; with 9,350,222 subreads and the average length of subreads, 1640 bp was counted. With 99.99% accuracy, after clustering, 132,536 transcripts and 78,559 genes were detected. All unique SMART transcripts were annotated in seven functional databases. 4180 TFs (including transcript regulators) and 7289 LncRNAs were predicted. The results of RNA-seq were confirmed with qRT–PCR analysis. Illumina sequencing of leaves and roots of two purslane genotypes was carried out. Amounts of differential expression genes and related KEGG pathways were found. The expression profiles of related genes in the biosynthesis of unsaturated fatty acids pathway in leaves and roots of two genotypes of purslane were analyzed. Differential expression of genes in this pathway built the foundation of ω-3 fatty acid accumulation in different organs and genotypes of purslane. The aforementioned results provide sequence information and may be a valuable resource for whole-genome sequencing of purslane in the future.

Physiological changes and transcriptome profiling in Saccharum spontaneum L. leaf under water stress and re-watering conditions

As the polyploidy progenitor of modern sugarcane, Saccharum spontaneum is considered to be a valuable resistance source to various biotic and abiotic stresses. However, little has been reported on the mechanism of drought tolerance in S. spontaneum. Herein, the physiological changes of S. spontaneum GXS87-16 at three water-deficit levels (mild, moderate, and severe) and after re-watering during the elongation stage were investigated. RNA sequencing was utilized for global transcriptome profiling of GXS87-16 under severe drought and re-watered conditions. There were significant alterations in the physiological parameters of GXS87-16 in response to drought stress and then recovered differently after re-watering. A total of 1569 differentially expressed genes (DEGs) associated with water stress and re-watering were identified. Notably, the majority of the DEGs were induced by stress. GO functional annotations and KEGG pathway analysis assigned the DEGs to 47 GO categories and 93 pathway categories. The pathway categories were involved in various processes, such as RNA transport, mRNA surveillance, plant hormone signal transduction, and plant-pathogen interaction. The reliability of the RNA-seq results was confirmed by qRT-PCR. This study shed light on the regulatory processes of drought tolerance in S. spontaneum and identifies useful genes for genetic improvement of drought tolerance in sugarcane.

Integration of Transcriptional and Post-transcriptional Analysis Revealed the Early Response Mechanism of Sugarcane to Cold Stress

Cold stress causes major losses to sugarcane production, yet the precise molecular mechanisms that cause losses due to cold stress are not well-understood. To survey miRNAs and genes involved in cold tolerance, RNA-seq, miRNA-seq, and integration analyses were performed on Saccharum spontaneum. Results showed that a total of 118,015 genes and 6,034 of these differentially expressed genes (DEGs) were screened. Protein–protein interaction (PPI) analyses revealed that ABA signaling via protein phosphatase 2Cs was the most important signal transduction pathway and late embryogenesis abundant protein was the hub protein associated with adaptation to cold stress. Furthermore, a total of 856 miRNAs were identified in this study and 109 of them were differentially expressed in sugarcane responding to cold stress. Most importantly, the miRNA–gene regulatory networks suggested the complex post-transcriptional regulation in sugarcane under cold stress, including 10 miRNAs−42 genes, 16 miRNAs−70 genes, and three miRNAs−18 genes in CT vs. LT0.5, CT vs. LT1, and CT0.5 vs. LT1, respectively. Specifically, key regulators from 16 genes encoding laccase were targeted by novel-Chr4C_47059 and Novel-Chr4A_40498, while five LRR-RLK genes were targeted by Novel-Chr6B_65233 and Novel-Chr5D_60023, 19 PPR repeat proteins by Novel-Chr5C_57213 and Novel-Chr5D_58065. Our findings suggested that these miRNAs and cell wall-related genes played vital regulatory roles in the responses of sugarcane to cold stress. Overall, the results of this study provide insights into the transcriptional and post-transcriptional regulatory network underlying the responses of sugarcane to cold stress.

Characterization and analysis of the transcriptome in Opisina arenosella from different developmental stages using single-molecule real-time transcript sequencing and RNA-seq

Opisina arenosella is one of the main pests harming coconut trees. To date, there have been few studies on the molecular genetics, biochemistry and physiology of O. arenosella at the transcriptional level, and there are no available reference genomes. Here, Illumina RNA sequencing combined with PacBio single-molecule real-time analysis was applied to study the transcriptome of this pest at different developmental stages, providing reference data for transcript expression analysis. Twelve samples of O. arenosella from different stages of development were sequenced using Illumina RNA sequencing, and the pooled RNA samples were sequenced with PacBio technology (Iso-Seq). A full-length transcriptome with 41,938 transcripts was captured, and the N50 and N90 lengths were 3543 bp and 1646 bp, respectively. A total of 36,925 transcripts were annotated in public databases, 6493 of which were long noncoding RNAs, while 2510 represented alternative splicing events. There were significant differences in the gene expression profiles at different developmental stages, with high levels of differential gene expression associated with growth, development, carbohydrate metabolism and immunity. This work provides resources and information for the study of the transcriptome and gene function of O. arenosella and provides a valuable foundation for understanding the changes in gene expression during development.

Current breeding and genomic approaches to enhance the cane and sugar productivity under abiotic stress conditions

Sugarcane (Saccharum spp.) crop is vulnerable to many abiotic stresses such as drought, salinity, waterlogging, cold and high temperature due to climate change. Over the past few decades new breeding and genomic approaches have been used to enhance the genotypic performance under abiotic stress conditions. In sugarcane, introgression of genes from wild species and allied genera for abiotic stress tolerance traits plays a significant role in the development of several stress-tolerant varieties. Moreover, the genomics and transcriptomics approaches have helped to elucidate the key genes/TFs and pathways involved in abiotic stress tolerance in sugarcane. Several novel miRNAs families /proteins or regulatory elements that are responsible for drought, salinity, and cold tolerance have been identified through high-throughput sequencing. The existing sugarcane monoploid genome sequence information opens new gateways and opportunities for researchers to improve the desired traits through efficient genome editing tools, such as the clustered regularly interspaced short palindromic repeat-Cas (CRISPR/Cas) system. TALEN mediated mutations in a highly conserved region of the caffeic acid O-methyltransferase (COMT) of sugarcane significantly reduces the lignin content in the cell wall which is amenable for biofuel production from lignocellulosic biomass. In this review, we focus on current breeding with genomic approaches and their substantial role in enhancing cane production under the abiotic stress conditions, which is expected to provide new insights to plant breeders and biotechnologists to modify their strategy in developing stress-tolerant sugarcane varieties, which can highlight the future demand of cane, bio-energy, and viability of sugar industries.

Enhanced sugar accumulation and regulated plant hormone signalling genes contribute to cold tolerance in hypoploid Saccharum spontaneum

BACKGROUND

Wild sugarcane Saccharum spontaneum plants vary in ploidy, which complicates the utilization of its germplasm in sugarcane breeding. Investigations on cold tolerance in relation to different ploidies in S. spontaneum may promote the exploitation of its germplasm and accelerate the improvement of sugarcane varieties.

RESULTS

A hypoploid clone 12-23 (2n = 54) and hyperploid clone 15-28 (2n = 92) of S. spontaneum were analysed under cold stress from morphological, physiological, and transcriptomic perspectives. Compared with clone 15-28, clone 12-23 plants had lower plant height, leaf length, internode length, stem diameter, and leaf width; depressed stomata and prominent bristles and papillae; and thick leaves with higher bulliform cell groups and thicker adaxial epidermis. Compared with clone 15-28, clone 12-23 showed significantly lower electrical conductivity, significantly higher water content, soluble protein content, and superoxide dismutase activity, and significantly higher soluble sugar content and peroxidase activity. Under cold stress, the number of upregulated genes and downregulated genes of clone 12-23 was higher than clone 15-28, and many stress response genes and pathways were affected and enriched to varying degrees, particularly sugar and starch metabolic pathways and plant hormone signalling pathways. Under cold stress, the activity of 6-phosphate glucose trehalose synthase, trehalose phosphate phosphatase, and brassinosteroid-signalling kinase and the content of trehalose and brassinosteroids of clone 12-23 increased.

CONCLUSIONS

Compared with hyperploid clone 15-28, hypoploid clone 12-23 maintained a more robust osmotic adjustment system through sugar accumulation and hormonal regulation, which resulted in stronger cold tolerance.

Genome-Wide Analysis of the DREB Subfamily in Saccharum spontaneum Reveals Their Functional Divergence During Cold and Drought Stresses

Drought and cold stresses are the main environmental factors that affect the yield of sugarcane, and DREB genes play very important roles in tolerance to drought, cold, and other environmental stresses. In this study, bioinformatics analysis was performed to characterize Saccharum spontaneum SsDREB genes. RNA sequencing (RNA-seq) was used to detect the expression profiles of SsDREBs induced by cold and drought stresses. According to our results, there are 110 SsDREB subfamily proteins in S. spontaneum, which can be classified into six groups; 106 of these genes are distributed among 29 chromosomes. Inter- and intraspecies synteny analyses suggested that all DREB groups have undergone gene duplication, highlighting the polyploid events that played an important role in the expansion of the DREB subfamily. Furthermore, RNA-seq results showed that 45 SsDREBs were up- or downregulated under cold stress; 35 of them were found to be involved in responding to drought stress. According to protein–protein interaction analysis, SsDREB100, SsDREB102, and SsDREB105 play key roles during the response to cold stress. These results reveal that functional divergence exists between collinear homologous genes or among common origin genes in the DREB subfamily of S. spontaneum. This study presents a comprehensive analysis and systematic understanding of the precise mechanism of SsDREBs in response to abiotic stress and will lead to improvements in sugarcane.

Transcriptome analysis reveals the mechanism underlying improved glutathione biosynthesis and secretion in Candida utilis during selenium enrichment

The effect of sodium selenite on batch culture of Candida utilis CCTCC M 209298 was investigated. Cell growth was inhibited while glutathione biosynthesis and secretion were improved during selenium enrichment. To reveal the mechanism underlying the decrease in biomass and the increase in glutathione, both metabolic flux analysis of key intermediates involved in glutathione metabolic pathway and transcriptome analysis of C. utilis by RNA-seq were carried out for selenized cells and the control without selenium enrichment. Results indicated that sodium selenite decreased carbon fluxes towards biomass but increased fluxes towards amino acids for the biosynthesis of glutathione and related amino acids. Selenium enrichment down-regulated a large number of genes involved in cell components and the cell cycle, resulting in decreased biomass as well as increased cell permeability. Moreover, several genes associated with transportation, binding, and mitochondrial and ribosomal functions for energy metabolism and protein synthesis were up-regulated in the presence of sodium selenite. All of these results disclosed the physiological response of C. utilis to sodium selenite.

Isolation and characterization of water-deficit stress-responsive α-expansin 1 (EXPA1) gene from Saccharum complex

In this study, full-length (1282-1330 bp) α-expansin 1 (EXPA1) gene from three different accessions belonging to Saccharum complex (Saccharum officinarum-SoEXPA1, Erianthus arundinaceus-EaEXPA1, and Saccharum spp. hybrid-ShEXPA1) was isolated using RAGE technique and characterized. The intronic and coding regions of isolated expansin genes ranged between 526-568 and 756-762 bp, respectively. An open reading frame encoding a polypeptide of 252 amino acids was obtained from S. officinarum and commercial sugarcane hybrid, whereas 254 amino acids were obtained in E. arundinaceus, a wild relative of Saccharum. Bioinformatics analysis of deduced protein revealed the presence of specific signature sequences and conserved amino acid residues crucial for the functioning of the protein. The predicted physicochemical characterization showed that the protein is stable in nature with instability index (II) value less than 40 and also clearly shown the dominance of random coil in the protein structure. Phylogenetic analysis revealed high conservation of EXPA1 among Saccharum complex and related crop species, Sorghum bicolor and Zea mays. The docking study of EXPA1 protein showed the interaction with xylose, which is present in xyloglucan of plant cell wall, elucidated the role of the expansin proteins in plant cell wall modification. This was further supported by the subcellular localization experiment in which it is clearly seen that the expansin protein localizes in the cell wall. Relative expression analysis of EXPA1 gene in Saccharum complex during drought stress showed high expression of the EaEXPA1 in comparison with SoEXPA1 and ShEXPA1 indicating possible role of EaEXPA1 in increased water-deficit stress tolerance in E. arundinaceus. These results suggest the potential use of EXPA1 for increasing the water-deficient stress tolerance levels in crop plants.

Uncovering full-length transcript isoforms of sugarcane cultivar Khon Kaen 3 using single-molecule long-read sequencing

Background

Sugarcane is an important global food crop and energy resource. To facilitate the sugarcane improvement program, genome and gene information are important for studying traits at the molecular level. Most currently available transcriptome data for sugarcane were generated using second-generation sequencing platforms, which provide short reads. The de novo assembled transcripts from these data are limited in length, and hence may be incomplete and inaccurate, especially for long RNAs.

Methods

We generated a transcriptome dataset of leaf tissue from a commercial Thai sugarcane cultivar Khon Kaen 3 (KK3) using PacBio RS II single-molecule long-read sequencing by the Iso-Seq method. Short-read RNA-Seq data were generated from the same RNA sample using the Ion Proton platform for reducing base calling errors.

Results

A total of 119,339 error-corrected transcripts were generated with the N50 length of 3,611 bp, which is on average longer than any previously reported sugarcane transcriptome dataset. 110,253 sequences (92.4%) contain an open reading frame (ORF) of at least 300 bp long with ORF N50 of 1,416 bp. The mean lengths of 5′ and 3′ untranslated regions in 73,795 sequences with complete ORFs are 1,249 and 1,187 bp, respectively. 4,774 transcripts are putatively novel full-length transcripts which do not match with a previous Iso-Seq study of sugarcane. We annotated the functions of 68,962 putative full-length transcripts with at least 90% coverage when compared with homologous protein coding sequences in other plants.

Discussion

The new catalog of transcripts will be useful for genome annotation, identification of splicing variants, SNP identification, and other research pertaining to the sugarcane improvement program. The putatively novel transcripts suggest unique features of KK3, although more data from different tissues and stages of development are needed to establish a reference transcriptome of this cultivar.

Transcriptome analysis highlights key differentially expressed genes involved in cellulose and lignin biosynthesis of sugarcane genotypes varying in fiber content

Sugarcane (Saccharum spp. hybrids) is a potential lignocellulosic feedstock for biofuel production due to its exceptional biomass accumulation ability, high convertible carbohydrate content and a favorable energy input/output ratio. Genetic modification of biofuel traits to improve biomass conversion requires an understanding of the regulation of carbohydrate and lignin biosynthesis. RNA-Seq was used to investigate the transcripts differentially expressed between the immature and mature tissues of the sugarcane genotypes varying in fiber content. Most of the differentially expressed transcripts were found to be down-regulated during stem maturation, highlighting their roles in active secondary cell-wall development in the younger tissues of both high and low fiber genotypes. Several cellulose synthase genes (including CesA2, CesA4, CesA7 and COBRA-like protein), lignin biosynthesis-related genes (ρ-coumarate 3-hydroxylase, ferulate 5-hydroxylase, cinnamyl alcohol dehydrogenase and gentiobiase) and transcription regulators for the secondary cell-wall synthesis (including LIM, MYB, PLATZ, IAA24, C2H2 and C2C2 DOF zinc finger gene families) were exclusively differentially expressed between immature and mature tissues of high fiber genotypes. These findings reveal target genes for subsequent research on the regulation of cellulose and lignin metabolism.

Differential gene expression profiling through transcriptome approach of Saccharum spontaneum L. under low temperature stress reveals genes potentially involved in cold acclimation

Sugarcane (Saccharum sp.) is predominantly grown in both tropics and subtropics in India, and the subtropics alone contribute more than half of sugarcane production. Sugarcane active growth period in subtropics is restricted to 8-9 months mainly due to winter's low temperature stress prevailing during November to February every year. Being a commercial crop, tolerance to low temperature is important in sugarcane improvement programs. Development of cold tolerant sugarcane varieties require a deep knowledge on molecular mechanism naturally adapted by cold tolerant genotypes during low temperature stress. To understand gene regulation under low temperature stress, control and stressed (10 °C, 24 h) leaf samples of cold tolerant S. spontaneum IND 00-1037 collected from high altitude region in Arunachal Pradesh were used for transcriptome analysis using the Illumina NextSeq 500 platform with paired-end sequencing method. Raw reads of 5.1 GB (control) and 5.3 GB (stressed) obtained were assembled using trinity and annotated with UNIPROT, KEGG, GO, COG and SUCEST databases, and transcriptome was validated using qRT-PCR. The differential gene expression (DGE) analysis showed that 2583 genes were upregulated and 3302 genes were down-regulated upon low temperature stress. A total of 170 cold responsive transcriptional factors belonging to 30 families were differentially regulated. CBF6 (C-binding factor), a DNA binding transcriptional activation protein associated with cold acclimation and freezing tolerance was differentially upregulated. Many low temperature responsive genes involved in various metabolic pathways, viz. cold sensing through membrane fluidity, calcium and lipid signaling genes, MAP kinases, phytohormone signaling and biosynthetic genes, antioxidative enzymes, membrane and cellular stabilizing genes, genes involved in biosynthesis of polyunsaturated fatty acids, chaperones, LEA proteins, soluble sugars, osmoprotectants, lignin and pectin biosynthetic genes were also differentially upregulated. Potential cold responsive genes and transcriptional factors involved in cold tolerance mechanism in cold tolerant S. spontaneum IND 00-1037 were identified. Together, this study provides insights into the cold tolerance to low temperature stress in S. spontaneum, thus opening applications in the genetic improvement of cold stress tolerance in sugarcane.

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

Sugarcane biomass has been used for sugar, bioenergy and biomaterial production. The majority of the sugarcane biomass comes from the culm, which makes it important to understand the genetic control of biomass production in this part of the plant. A meta-transcriptome of the culm was obtained in an earlier study by using about one billion paired-end (150 bp) reads of deep RNA sequencing of samples from 20 diverse sugarcane genotypes and combining de novo assemblies from different assemblers and different settings. Although many genes could be recovered, this resulted in a large combined assembly which created the need for clustering to reduce transcript redundancy while maintaining gene content. Here, we present a comprehensive analysis of the effect of different assembly settings and clustering methods on de novo assembly, annotation and transcript profiling focusing especially on the coding transcripts from the highly polyploid sugarcane genome. The new coding sequence-based transcript clustering resulted in a better representation of transcripts compared to the earlier approach, having 121,987 contigs, which included 78,052 main and 43,935 alternative transcripts. About 73%, 67%, 61% and 10% of the transcriptome was annotated against the NCBI NR protein database, GO terms, orthologous groups and KEGG orthologies, respectively. Using this set for a differential gene expression analysis between the young and mature sugarcane culm tissues, a total of 822 transcripts were found to be differentially expressed, including key transcripts involved in sugar/fiber accumulation in sugarcane. In the context of the lack of a whole genome sequence for sugarcane, the availability of a well annotated culm-derived meta-transcriptome through deep sequencing provides useful information on coding genes specific to the sugarcane culm and will certainly contribute to understanding the process of carbon partitioning, and biomass accumulation in the sugarcane culm.

Transcriptomic characterization and potential marker development of contrasting sugarcane cultivars

Sugarcane (Saccharum officinarum L.) is an important crop for sugar production and bioenergy worldwide. In this study, we performed transcriptome sequencing for six contrasting sugarcane genotypes involved in leaf abscission, tolerance to pokkah boeng disease and drought stress. More than 465 million high-quality reads were generated, which were de novo assembled into 93,115 unigenes. Based on a similarity search, 43,526 (46.74%) unigenes were annotated against at least one of the public databases. Functional classification analyses showed that these unigenes are involved in a wide range of metabolic pathways. Comparative transcriptome analysis revealed that many unigenes involved in response to abscisic acid and ethylene were up-regulated in the easy leaf abscission genotype, and unigenes associated with response to jasmonic acid and salicylic acid were up-regulated in response to the pokkah boeng disease in the tolerance genotype. Moreover, unigenes related to peroxidase, antioxidant activity and signal transduction were up-regulated in response to drought stress in the tolerant genotype. Finally, we identified a number of putative markers, including 8,630 simple sequence repeats (SSRs) and 442,152 single-nucleotide polymorphisms (SNPs). Our data will be important resources for future gene discovery, molecular marker development, and genome studies in sugarcane.

A survey of the complex transcriptome from the highly polyploid sugarcane genome using full-length isoform sequencing and de novo assembly from short read sequencing

BACKGROUND

Despite the economic importance of sugarcane in sugar and bioenergy production, there is not yet a reference genome available. Most of the sugarcane transcriptomic studies have been based on Saccharum officinarum gene indices (SoGI), expressed sequence tags (ESTs) and de novo assembled transcript contigs from short-reads; hence knowledge of the sugarcane transcriptome is limited in relation to transcript length and number of transcript isoforms.

RESULTS

The sugarcane transcriptome was sequenced using PacBio isoform sequencing (Iso-Seq) of a pooled RNA sample derived from leaf, internode and root tissues, of different developmental stages, from 22 varieties, to explore the potential for capturing full-length transcript isoforms. A total of 107,598 unique transcript isoforms were obtained, representing about 71% of the total number of predicted sugarcane genes. The majority of this dataset (92%) matched the plant protein database, while just over 2% was novel transcripts, and over 2% was putative long non-coding RNAs. About 56% and 23% of total sequences were annotated against the gene ontology and KEGG pathway databases, respectively. Comparison with de novo contigs from Illumina RNA-Sequencing (RNA-Seq) of the internode samples from the same experiment and public databases showed that the Iso-Seq method recovered more full-length transcript isoforms, had a higher N50 and average length of largest 1,000 proteins; whereas a greater representation of the gene content and RNA diversity was captured in RNA-Seq. Only 62% of PacBio transcript isoforms matched 67% of de novo contigs, while the non-matched proportions were attributed to the inclusion of leaf/root tissues and the normalization in PacBio, and the representation of more gene content and RNA classes in the de novo assembly, respectively. About 69% of PacBio transcript isoforms and 41% of de novo contigs aligned with the sorghum genome, indicating the high conservation of orthologs in the genic regions of the two genomes.

CONCLUSIONS

The transcriptome dataset should contribute to improved sugarcane gene models and sugarcane protein predictions; and will serve as a reference database for analysis of transcript expression in sugarcane.