Analysis and functional annotation of an expressed sequence tag collection for tropical crop sugarcane

Long read transcriptome sequencing of a sugarcane hybrid and its progenitors, Saccharum officinarum and S. spontaneum

Commercial sugarcane hybrids are derivatives from Saccharum officinarum and Saccharum spontaneum hybrids containing the full complement of S. officinarum and a few S. spontaneum chromosomes and recombinants with favorable agronomic characters from both the species. The combination of the two sub-genomes in varying proportions in addition to the recombinants presents a challenge in the study of gene expression and regulation in the hybrid. We now report the transcriptome analysis of the two progenitor species and a modern commercial sugarcane hybrid through long read sequencing technology. Transcripts were profiled in the two progenitor species S. officinarum (Black Cheribon), and S. spontaneum (Coimbatore accession) and a recent high yielding, high sugar variety Co 11015. The composition and contribution of the progenitors to a hybrid with respect to sugar, biomass, and disease resistance were established. Sugar related transcripts originated from S. officinarum while several stress and senescence related transcripts were from S. spontaneum in the hybrid. The hybrid had a higher number of transcripts related to sugar transporters, invertases, transcription factors, trehalose, UDP sugars, and cellulose than the two progenitor species. Both S. officinarum and the hybrid had an abundance of novel genes like sugar phosphate translocator, while S. spontaneum had just one. In general, the hybrid shared a larger number of transcripts with S. officinarum than with S. spontaneum, reflecting the genomic contribution, while the progenitors shared very few transcripts between them. The common isoforms among the three genotypes and unique isoforms specific to each genotype indicate that there is a high scope for improvement of the modern hybrids by utilizing novel gene isoforms from the progenitor species.

Functional and comparative analysis of THI1 gene in grasses with a focus on sugarcane

De novo synthesis of thiamine (vitamin B1) in plants depends on the action of thiamine thiazole synthase, which synthesizes the thiazole ring, and is encoded by the THI1 gene. Here, we investigated the evolution and diversity of THI1 in Poaceae, where C4 and C3 photosynthetic plants co-evolved. An ancestral duplication of THI1 is observed in Panicoideae that remains in many modern monocots, including sugarcane. In addition to the two sugarcane copies (ScTHI1-1 and ScTHI1-2), we identified ScTHI1-2 alleles showing differences in their sequence, indicating divergence between ScTHI1-2a and ScTHI1-2b. Such variations are observed only in the Saccharum complex, corroborating the phylogeny. At least five THI1 genomic environments were found in Poaceae, two in sugarcane, M. sinensis, and S. bicolor. The THI1 promoter in Poaceae is highly conserved at 300 bp upstream of the start codon ATG and has cis-regulatory elements that putatively bind to transcription factors associated with development, growth, development and biological rhythms. An experiment set to compare gene expression levels in different tissues across the sugarcane R570 life cycle showed that ScTHI1-1 was expressed mainly in leaves regardless of age. Furthermore, ScTHI1 displayed relatively high expression levels in meristem and culm, which varied with the plant age. Finally, yeast complementation studies with THI4-defective strain demonstrate that only ScTHI1-1 and ScTHI1-2b isoforms can partially restore thiamine auxotrophy, albeit at a low frequency. Taken together, the present work supports the existence of multiple origins of THI1 harboring genomic regions in Poaceae with predicted functional redundancy. In addition, it questions the contribution of the levels of the thiazole ring in C4 photosynthetic plant tissues or potentially the relevance of the THI1 protein activity.

Antisense Transcription in Plants: A Systematic Review and an Update on cis-NATs of Sugarcane

Initially, natural antisense transcripts (NATs, natRNAs, or asRNAs) were considered repressors; however, their functions in gene regulation are diverse. Positive, negative, or neutral correlations to the cognate gene expression have been noted. Although the first studies were published about 50 years ago, there is still much to be investigated regarding antisense transcripts in plants. A systematic review of scientific publications available in the Web of Science databases was conducted to contextualize how the studying of antisense transcripts has been addressed. Studies were classified considering three categories: “Natural antisense” (208), artificial antisense used in “Genetic Engineering” (797), or “Natural antisense and Genetic Engineering”-related publications (96). A similar string was used for a systematic search in the NCBI Gene database. Of the 1132 antisense sequences found for plants, only 0.8% were cited in PubMed and had antisense information confirmed. This value was the lowest when compared to fungi (2.9%), bacteria (2.3%), and mice (54.1%). Finally, we present an update for the cis-NATs identified in Saccharum spp. Of the 1413 antisense transcripts found in different experiments, 25 showed concordant expressions, 22 were discordant, 1264 did not correlate with the cognate genes, and 102 presented variable results depending on the experiment.

A short review on sugarcane: its domestication, molecular manipulations and future perspectives

Sugarcane (Saccharum spp.) is a special crop plant that underwent anthropogenic evolution from a wild grass species to an important food, fodder, and energy crop. Unlike any other grass species which were selected for their kernels, sugarcane was selected for its high stem sucrose accumulation. Flowering in sugarcane is not favored since flowering diverts the stored sugar resources for the reproductive and developmental energy needs. Cultivars are vegetatively propagated and sugarcane breeding is still essentially focused on conventional methods, since the knowledge of sugarcane genetics has lagged that of other major crops. Cultivar improvement has been extremely challenging due to its polyploidy and aneuploidy nature derived from a few interspecific hybridizations between Saccharum officinarum and Saccharum spontaneum, revealing the coexistence of two distinct genome organization modes in the modern variety. Alongside implementation of modern agricultural techniques, generation of hybrid clones, transgenics and genome edited events will help to meet the ever-growing bioenergy needs. Additionally, there are two common biotechnological approaches to improve plant stress tolerance, which includes marker-assisted selection (MAS) and genetic transformation. During the past two decades, the use of molecular approaches has contributed greatly to a better understanding of the genetic and biochemical basis of plant stress-tolerance and in some cases, it led to the development of plants with enhanced tolerance to abiotic stress. Hence, this review mainly intends on the events that shaped the sugarcane as what it is now and what challenges ahead and measures taken to further improve its yield, production and maximize utilization to beat the growing demands.

Root nitrate uptake in sugarcane (Saccharum spp.) is modulated by transcriptional and presumably posttranscriptional regulation of the NRT2.1/NRT3.1 transport system

KEY MESSAGE

Nitrate uptake in sugarcane roots is regulated at the transcriptional and posttranscriptional levels based on the physiological status of the plant and is likely a determinant mechanism for discrimination against nitrate. Sugarcane (Saccharum spp.) is one of the most suitable energy crops for biofuel feedstock, but the reduced recovery of nitrogen (N) fertilizer by sugarcane roots increases the crop carbon footprint. The low nitrogen use efficiency (NUE) of sugarcane has been associated with the significantly low nitrate uptake, which limits the utilization of the large amount of nitrate available in agricultural soils. To understand the regulation of nitrate uptake in sugarcane roots, we identified the major canonical nitrate transporter genes (NRTs-NITRATE TRANSPORTERS) and then determined their expression profiles in roots under contrasting N conditions. Correlation of gene expression with ¹⁵N-nitrate uptake revealed that under N deprivation or inorganic N (ammonium or nitrate) supply in N-sufficient roots, the regulation of ScNRT2.1 and ScNRT3.1 expression is the predominant mechanism for the modulation of the activity of the nitrate high-affinity transport system. Conversely, in N-deficient roots, the induction of ScNRT2.1 and ScNRT3.1 transcription is not correlated with the marked repression of nitrate uptake in response to nitrate resupply or high N provision, which suggested the existence of a posttranscriptional regulatory mechanism. Our findings suggested that high-affinity nitrate uptake is regulated at the transcriptional and presumably at the posttranscriptional levels based on the physiological N status and that the regulation of NRT2.1 and NRT3.1 activity is likely a determinant mechanism for the discrimination against nitrate uptake observed in sugarcane roots, which contributes to the low NUE in this crop species.

Bioinformatic analyses to uncover genes involved in trehalose metabolism in the polyploid sugarcane

Trehalose-6-phosphate (T6P) is an intermediate of trehalose biosynthesis that plays an essential role in plant metabolism and development. Here, we comprehensively analyzed sequences from enzymes of trehalose metabolism in sugarcane, one of the main crops used for bioenergy production. We identified protein domains, phylogeny, and in silico expression levels for all classes of enzymes. However, post-translational modifications and residues involved in catalysis and substrate binding were analyzed only in trehalose-6-phosphate synthase (TPS) sequences. We retrieved 71 putative full-length TPS, 93 trehalose-6-phosphate phosphatase (TPP), and 3 trehalase (TRE) of sugarcane, showing all their conserved domains, respectively. Putative TPS (Classes I and II) and TPP sugarcane sequences were categorized into well-known groups reported in the literature. We measured the expression levels of the sequences from one sugarcane leaf transcriptomic dataset. Furthermore, TPS Class I has specific N-glycosylation sites inserted in conserved motifs and carries catalytic and binding residues in its TPS domain. Some of these residues are mutated in TPS Class II members, which implies loss of enzyme activity. Our approach retrieved many homo(eo)logous sequences for genes involved in trehalose metabolism, paving the way to discover the role of T6P signaling in sugarcane.

Analysis of the PEBP gene family and identification of a novel FLOWERING LOCUS T orthologue in sugarcane

Sugarcane (Saccharum spp.) is an important economic crop for both sugar and biomass, the yields of which are negatively affected by flowering. The molecular mechanisms controlling flowering in sugarcane are nevertheless poorly understood. RNA-seq data analysis and database searches have enabled a comprehensive description of the PEBP gene family in sugarcane. It is shown to consist of at least 13 FLOWERING LOCUS T (FT)-like genes, two MOTHER OF FT AND TFL (MFT)-like genes, and four TERMINAL FLOWER (TFL)-like genes. As expected, these genes all show very high homology to their corresponding genes in Sorghum, and also to FT-like, MFT-like, and TFL-like genes in maize, rice, and Arabidopsis. Functional analysis in Arabidopsis showed that the sugarcane ScFT3 gene can rescue the late flowering phenotype of the Arabidopsis ft-10 mutant, whereas ScFT5 cannot. High expression levels of ScFT3 in leaves of short day-induced sugarcane plants coincided with initial stages of floral induction in the shoot apical meristem as shown by histological analysis of meristem dissections. This suggests that ScFT3 is likely to play a role in floral induction in sugarcane; however, other sugarcane FT-like genes may also be involved in the flowering process.

A Sugarcane G-Protein-Coupled Receptor, ShGPCR1, Confers Tolerance to Multiple Abiotic Stresses

Sugarcane (Saccharum spp.) is a prominent source of sugar and serves as bioenergy/biomass feedstock globally. Multiple biotic and abiotic stresses, including drought, salinity, and cold, adversely affect sugarcane yield. G-protein-coupled receptors (GPCRs) are components of G-protein-mediated signaling affecting plant growth, development, and stress responses. Here, we identified a GPCR-like protein (ShGPCR1) from sugarcane and energy cane (Saccharum spp. hybrids) and characterized its function in conferring tolerance to multiple abiotic stresses. ShGPCR1 protein sequence contained nine predicted transmembrane (TM) domains connected by four extracellular and four intracellular loops, which could interact with various ligands and heterotrimeric G proteins in the cells. ShGPCR1 sequence displayed other signature features of a GPCR, such as a putative guanidine triphosphate (GTP)-binding domain, as well as multiple myristoylation and protein phosphorylation sites, presumably important for its biochemical function. Expression of ShGPCR1 was upregulated by drought, salinity, and cold stresses. Subcellular imaging and calcium (Ca²⁺) measurements revealed that ShGPCR1 predominantly localized to the plasma membrane and enhanced intracellular Ca²⁺ levels in response to GTP, respectively. Furthermore, constitutive overexpression of ShGPCR1 in sugarcane conferred tolerance to the three stressors. The stress-tolerance phenotype of the transgenic lines corresponded with activation of multiple drought-, salinity-, and cold-stress marker genes, such as Saccharum spp. LATE EMBRYOGENESIS ABUNDANT, DEHYDRIN, DROUGHT RESPONSIVE 4, GALACTINOL SYNTHASE, ETHYLENE RESPONSIVE FACTOR 3, SALT OVERLY SENSITIVE 1, VACUOLAR Na+/H+ ANTIPORTER 1, NAM/ATAF1/2/CUC2, COLD RESPONSIVE FACTOR 2, and ALCOHOL DEHYDROGENASE 3. We suggest that ShGPCR1 plays a key role in conferring tolerance to multiple abiotic stresses, and the engineered lines may be useful to enhance sugarcane production in marginal environments with fewer resources.

Modern Biotechnologies: Innovative and Sustainable Approaches for the Improvement of Sugarcane Tolerance to Environmental Stresses

Sugarcane (Saccharum spp.) is one of the most important industrial cash crops, contributing to the world sugar industry and biofuel production. It has been cultivated and improved from prehistoric times through natural selection and conventional breeding and, more recently, using the modern tools of genetic engineering and biotechnology. However, the heterogenicity, complex poly-aneuploid genome and susceptibility of sugarcane to different biotic and abiotic stresses represent impediments that require us to pay greater attention to the improvement of the sugarcane crop. Compared to traditional breeding, recent advances in breeding technologies (molecular marker-assisted breeding, sugarcane transformation, genome-editing and multiple omics technologies) can potentially improve sugarcane, especially against environmental stressors. This article will focus on efficient modern breeding technologies, which provide crucial clues for the engineering of sugarcane cultivars resistant to environmental stresses.

Overexpression of a Sugarcane BAHD Acyltransferase Alters Hydroxycinnamate Content in Maize Cell Wall

The purification of hydroxycinnamic acids [p-coumaric acid (pCA) and ferulic acid (FA)] from grass cell walls requires high-cost processes. Feedstocks with increased levels of one hydroxycinnamate in preference to the other are therefore highly desirable. We identified and conducted expression analysis for nine BAHD acyltransferase ScAts genes from sugarcane. The high conservation of AT10 proteins, together with their similar gene expression patterns, supported a similar role in distinct grasses. Overexpression of ScAT10 in maize resulted in up to 75% increase in total pCA content. Mild hydrolysis and derivatization followed by reductive cleavage (DFRC) analysis showed that pCA increase was restricted to the hemicellulosic portion of the cell wall. Furthermore, total FA content was reduced up to 88%, resulting in a 10-fold increase in the pCA/FA ratio. Thus, we functionally characterized a sugarcane gene involved in pCA content on hemicelluloses and generated a C4 plant that is promising for valorizing pCA production in biorefineries.

Applying Molecular Phenotyping Tools to Explore Sugarcane Carbon Potential

Sugarcane (Saccharum spp.), a C₄ grass, has a peculiar feature: it accumulates, gradient-wise, large amounts of carbon (C) as sucrose in its culms through a complex pathway. Apart from being a sustainable crop concerning C efficiency and bioenergetic yield per hectare, sugarcane is used as feedstock for producing ethanol, sugar, high-value compounds, and products (e.g., polymers and succinate), and bioelectricity, earning the title of the world's leading biomass crop. Commercial cultivars, hybrids bearing high levels of polyploidy, and aneuploidy, are selected from a large number of crosses among suitable parental genotypes followed by the cloning of superior individuals among the progeny. Traditionally, these classical breeding strategies have been favoring the selection of cultivars with high sucrose content and resistance to environmental stresses. A current paradigm change in sugarcane breeding programs aims to alter the balance of C partitioning as a means to provide more plasticity in the sustainable use of this biomass for metabolic engineering and green chemistry. The recently available sugarcane genetic assemblies powered by data science provide exciting perspectives to increase biomass, as the current sugarcane yield is roughly 20% of its predicted potential. Nowadays, several molecular phenotyping tools can be applied to meet the predicted sugarcane C potential, mainly targeting two competing pathways: sucrose production/storage and biomass accumulation. Here we discuss how molecular phenotyping can be a powerful tool to assist breeding programs and which strategies could be adopted depending on the desired final products. We also tackle the advances in genetic markers and mapping as well as how functional genomics and genetic transformation might be able to improve yield and saccharification rates. Finally, we review how "omics" advances are promising to speed up plant breeding and reach the unexplored potential of sugarcane in terms of sucrose and biomass production.

The Wild Sugarcane and Sorghum Kinomes: Insights Into Expansion, Diversification, and Expression Patterns

The protein kinase (PK) superfamily is one of the largest superfamilies in plants and the core regulator of cellular signaling. Despite this substantial importance, the kinomes of sugarcane and sorghum have not been profiled. Here, we identified and profiled the complete kinomes of the polyploid Saccharum spontaneum (Ssp) and Sorghum bicolor (Sbi), a close diploid relative. The Sbi kinome was composed of 1,210 PKs; for Ssp, we identified 2,919 PKs when disregarding duplications and allelic copies, and these were related to 1,345 representative gene models. The Ssp and Sbi PKs were grouped into 20 groups and 120 subfamilies and exhibited high compositional similarities and evolutionary divergences. By utilizing the collinearity between the species, this study offers insights into Sbi and Ssp speciation, PK differentiation and selection. We assessed the PK subfamily expression profiles via RNA-Seq and identified significant similarities between Sbi and Ssp. Moreover, coexpression networks allowed inference of a core structure of kinase interactions with specific key elements. This study provides the first categorization of the allelic specificity of a kinome and offers a wide reservoir of molecular and genetic information, thereby enhancing the understanding of Sbi and Ssp PK evolutionary history.

Sugarcane cystatins: From discovery to biotechnological applications

Phytocystatins are tight-binding cysteine protease inhibitors produced by plants. The first phytocystatin described was isolated from Oryza sativa and, since then, cystatins from several plant species were reported, including from sugarcane. Sugarcane cystatins were unraveled in Sugarcane EST project database, after sequencing of cDNA libraries from various sugarcane tissues at different developmental stages and six sugarcane cystatins were cloned, expressed and characterized (CaneCPI-1 to CaneCPI-6). These recombinant proteins were produced in different expression systems and inhibited several cysteine proteases, including human cathepsins B and L, which can be involved in pathologies, such as cancer. In this review, we summarize a comprehensive history of all sugarcane cystatins, presenting an updated phylogenetic analysis; chromosomal localization, and genomic organization. We also present protein docking of CaneCPI-5 in the active site of human cathepsin B, insights about canecystatins structures; recombinant expression in different systems, comparison of their inhibitory activities against human cysteine cathepsins B, K, L, S, V, falcipains from Plasmodium falciparum and a cathepsin L-like from the sugarcane weevil Sphenophorus levis; and enlighten their potential and current applications in agriculture and health.

Amino Acid and Carbohydrate Metabolism Are Coordinated to Maintain Energetic Balance during Drought in Sugarcane

The ability to expand crop plantations without irrigation is a major goal to increase agriculture sustainability. To achieve this end, we need to understand the mechanisms that govern plant growth responses under drought conditions. In this study, we combined physiological, transcriptomic, and genomic data to provide a comprehensive picture of drought and recovery responses in the leaves and roots of sugarcane. Transcriptomic profiling using oligoarrays and RNA-seq identified 2898 (out of 21,902) and 46,062 (out of 373,869) transcripts as differentially expressed, respectively. Co-expression analysis revealed modules enriched in photosynthesis, small molecule metabolism, alpha-amino acid metabolism, trehalose biosynthesis, serine family amino acid metabolism, and carbohydrate transport. Together, our findings reveal that carbohydrate metabolism is coordinated with the degradation of amino acids to provide carbon skeletons to the tricarboxylic acid cycle. This coordination may help to maintain energetic balance during drought stress adaptation, facilitating recovery after the stress is alleviated. Our results shed light on candidate regulatory elements and pave the way to biotechnology strategies towards the development of drought-tolerant sugarcane plants.

Sugarcane mosaic virus mediated changes in cytosine methylation pattern and differentially transcribed fragments in resistance-contrasting sugarcane genotypes

Sugarcane mosaic virus (SCMV) is the causal agent of sugarcane mosaic disease (SMD) in Brazil; it is mainly controlled by using resistant cultivars. Studies on the changes in sugarcane transcriptome provided the first insights about the molecular basis underlying the genetic resistance to SMD; nonetheless, epigenetic modifications such as cytosine methylation is also informative, considering its roles in gene expression regulation. In our previous study, differentially transcribed fragments (DTFs) were obtained using cDNA-amplified fragment length polymorphism by comparing mock- and SCMV-inoculated plants from two sugarcane cultivars with contrasting responses to SMD. In this study, the identification of unexplored DTFs was continued while the same leaf samples were used to evaluate SCMV-mediated changes in the cytosine methylation pattern by using methylation-sensitive amplification polymorphism. This analysis revealed minor changes in cytosine methylation in response to SCMV infection, but distinct changes between the cultivars with contrasting responses to SMD, with higher hypomethylation events 24 and 72 h post-inoculation in the resistant cultivar. The differentially methylated fragments (DMFs) aligned with transcripts, putative promoters, and genomic regions, with a preponderant distribution within CpG islands. The transcripts found were associated with plant immunity and other stress responses, epigenetic changes, and transposable elements. The DTFs aligned with transcripts assigned to stress responses, epigenetic changes, photosynthesis, lipid transport, and oxidoreductases, in which the transcriptional start site is located in proximity with CpG islands and tandem repeats. Real-time quantitative polymerase chain reaction results revealed significant upregulation in the resistant cultivar of aspartyl protease and VQ protein, respectively, selected from DMF and DTF alignments, suggesting their roles in genetic resistance to SMD and supporting the influence of cytosine methylation in gene expression. Thus, we identified new candidate genes for further validation and showed that the changes in cytosine methylation may regulate important mechanisms underlying the genetic resistance to SMD.

Genome biology of the paleotetraploid perennial biomass crop Miscanthus

Miscanthus is a perennial wild grass that is of global importance for paper production, roofing, horticultural plantings, and an emerging highly productive temperate biomass crop. We report a chromosome-scale assembly of the paleotetraploid M. sinensis genome, providing a resource for Miscanthus that links its chromosomes to the related diploid Sorghum and complex polyploid sugarcanes. The asymmetric distribution of transposons across the two homoeologous subgenomes proves Miscanthus paleo-allotetraploidy and identifies several balanced reciprocal homoeologous exchanges. Analysis of M. sinensis and M. sacchariflorus populations demonstrates extensive interspecific admixture and hybridization, and documents the origin of the highly productive triploid bioenergy crop M. × giganteus. Transcriptional profiling of leaves, stem, and rhizomes over growing seasons provides insight into rhizome development and nutrient recycling, processes critical for sustainable biomass accumulation in a perennial temperate grass. The Miscanthus genome expands the power of comparative genomics to understand traits of importance to Andropogoneae grasses.

The perennial grass Miscanthus is a promising biomass crop. Here, via genomics and transcriptomics, the authors reveal its allotetraploid origin, characterize gene expression associated with rhizome development and nutrient recycling, and describe the hybrid origin of the triploid M. x giganteus.

Molecular Mechanisms Underlying Sugarcane Response to Aluminum Stress by RNA-Seq

Some metals are beneficial to plants and contribute to critical physiological processes. Some metals, however, are not. The presence of aluminum ions (Al³⁺) can be very toxic, especially in acidic soils. Considerable parts of the world's arable land are acidic in nature; mechanistically elucidating a plant's response to aluminum stress is critical to mitigating this stress and improving the quality of plants. To identify the genes involved in sugarcane response to aluminum stress, we generated 372 million paired-end RNA sequencing reads from the roots of CTC-2 and RB855453, which are two contrasting cultivars. Data normalization resulted in 162,161 contigs (contiguous sequences) and 97,335 genes from a de novo transcriptome assembly (trinity genes). A total of 4858 and 1307 differently expressed genes (DEGs) for treatment versus control were identified for the CTC-2 and RB855453 cultivars, respectively. The DEGs were annotated into 34 functional categories. The majority of the genes were upregulated in the CTC-2 (tolerant cultivar) and downregulated in RB855453 (sensitive cultivar). Here, we present the first root transcriptome of sugarcane under aluminum stress. The results and conclusions of this study are a crucial launch pad for future genetic and genomic studies of sugarcane. The transcriptome analysis shows that sugarcane tolerance to aluminum may be explained by an efficient detoxification mechanism combined with lateral root formation and activation of redox enzymes. We also present a hypothetical model for aluminum tolerance in the CTC-2 cultivar.

Role of NGS and SNP genotyping methods in sugarcane improvement programs

Sugarcane (Saccharum spp.) is one of the most economically significant crops because of its high sucrose content and it is a promising biomass feedstock for biofuel production. Sugarcane genome sequencing and analysis is a difficult task due to its heterozygosity and polyploidy. Long sequence read technologies, PacBio Single-Molecule Real-Time (SMRT) sequencing, the Illumina TruSeq, and the Oxford Nanopore sequencing could solve the problem of genome assembly. On the applications side, next generation sequencing (NGS) technologies played a major role in the discovery of single nucleotide polymorphism (SNP) and the development of low to high throughput genotyping platforms. The two mainstream high throughput genotyping platforms are the SNP microarray and genotyping by sequencing (GBS). This paper reviews the NGS in sugarcane genomics, genotyping methodologies, and the choice of these methods. Array-based SNP genotyping is robust, provides consistent SNPs, and relatively easier downstream data analysis. The GBS method identifies large scale SNPs across the germplasm. A combination of targeted GBS and array-based genotyping methods should be used to increase the accuracy of genomic selection and marker-assisted breeding.

Comparative analysis of drought-responsive transcriptomes of sugarcane genotypes with differential tolerance to drought

Water stress causes considerable yield losses in sugarcane. To investigate differentially expressed genes under water stress, two sugarcane genotypes were subjected to three water-deficit levels (mild, moderate, and severe) and subsequent recovery and leaf transcriptome was generated using Illumina NextSeq sequencing. Among the differentially expressed genes, the tolerant genotype Co 06022 generated 2970 unigenes (p ≤ 0.05, functionally known, non-redundant DEGs) at 2-day stress, and there was a progressive decrease in the expressed genes as the stress period increased with 2109 unigenes at 6-day stress and 2307 unigenes at 10-day stress. There was considerable reduction at recovery with 1334 unigenes expressed at 10 days after recovery. However, in the susceptible genotype Co 8021, the number of unigenes expressed at 2 days was lower (2025) than the tolerant genotype and a further reduction was seen at 6-day stress (1552). During recovery, more differentially expressed genes were observed in the susceptible cultivar indicating that the cultivar has to activate more functions/processes to recover from the damage caused by stress. Comparison of DEGs between all stages of stress and recovery in both genotypes revealed that, the commonly up- and down-regulated genes across different stages were approximately double in the tolerant genotype. The most enriched gene ontology classes were heme binding, peroxidase activity and metal ion binding in the biological process and response to oxidative stress, hydrogen peroxide catabolic process and response to stress in the molecular function category. The cellular component was enriched with DEGs involved in extracellular region followed by integral component of membrane. The KEGG pathway analysis revealed important metabolic activities and functionally important genes involved in mitigating water-deficit stress in both the varieties. In addition, several unannotated genes in important pathways were detected and together may provide novel insights into water-deficit tolerance mechanisms in sugarcane. The reliability of the observed expression patterns was confirmed by qRT-PCR. The results of this study will help to identify useful genes for improving drought tolerance in sugarcane.

Bringing to light the molecular evolution of GUX genes in plants

Hemicellulose and cellulose are essential polysaccharides for plant development and major components of cell wall. They are also an important energy source for the production of ethanol from plant biomass, but their conversion to fermentable sugars is hindered by the complex structure of cell walls. The glucuronic acid substitution of xylan (GUX) enzymes attach glucuronic acid to xylan, a major component of hemicellulose, decreasing the efficiency of enzymes used for ethanol production. Since loss-of-function gux mutants of Arabidopsis thaliana enhance enzyme accessibility and cell wall digestion without adverse phenotypes, GUX genes are potential targets for genetically improving energy crops. However, comprehensive identification of GUX in important species and their evolutionary history are largely lacking. Here, we identified putative GUX proteins using hidden Markov model searches with the GT8 domain and a GUX-specific motif, and inferred the phylogenetic relationship of 18 species with Maximum likelihood and Bayesian approaches. Each species presented a variable number of GUX, and their evolution can be explained by a mixture of divergent, concerted and birth-and-death evolutionary models. This is the first broad insight into the evolution of GUX gene family in plants and will potentially guide genetic and functional studies in species used for biofuel production.

Assembly of the 373k gene space of the polyploid sugarcane genome reveals reservoirs of functional diversity in the world's leading biomass crop

BACKGROUND

Sugarcane cultivars are polyploid interspecific hybrids of giant genomes, typically with 10-13 sets of chromosomes from 2 Saccharum species. The ploidy, hybridity, and size of the genome, estimated to have >10 Gb, pose a challenge for sequencing.

RESULTS

Here we present a gene space assembly of SP80-3280, including 373,869 putative genes and their potential regulatory regions. The alignment of single-copy genes in diploid grasses to the putative genes indicates that we could resolve 2-6 (up to 15) putative homo(eo)logs that are 99.1% identical within their coding sequences. Dissimilarities increase in their regulatory regions, and gene promoter analysis shows differences in regulatory elements within gene families that are expressed in a species-specific manner. We exemplify these differences for sucrose synthase (SuSy) and phenylalanine ammonia-lyase (PAL), 2 gene families central to carbon partitioning. SP80-3280 has particular regulatory elements involved in sucrose synthesis not found in the ancestor Saccharum spontaneum. PAL regulatory elements are found in co-expressed genes related to fiber synthesis within gene networks defined during plant growth and maturation. Comparison with sorghum reveals predominantly bi-allelic variations in sugarcane, consistent with the formation of 2 "subgenomes" after their divergence ∼3.8-4.6 million years ago and reveals single-nucleotide variants that may underlie their differences.

CONCLUSIONS

This assembly represents a large step towards a whole-genome assembly of a commercial sugarcane cultivar. It includes a rich diversity of genes and homo(eo)logous resolution for a representative fraction of the gene space, relevant to improve biomass and food production.

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.

Overexpression of an evolutionarily conserved drought-responsive sugarcane gene enhances salinity and drought resilience

BACKGROUND AND AIMS

Improving drought adaptation is more pressing for crops such as sugarcane, rice, wheat and maize, given the high dependence of these crops on irrigation. One option for enhancing adaptation to water limitation in plants is by transgenic approaches. An increasing number of genes that are associated with mechanisms used by plants to cope with water scarcity have been discovered. Genes encoding proteins with unknown functions comprise a relevant fraction of the genes that are modulated by drought. We characterized a gene in response to environmental stresses to gain insight into the unknown fraction of the sugarcane genome. Scdr2 (Sugarcane drought-responsive 2) encodes a small protein and shares highly conserved sequences within monocots, dicots, algae and fungi.

METHODS

Plants overexpressing the Scdr2 sugarcane gene were examined in response to salinity and drought. Measurements of the gas exchange parameters, germination rate, water content, dry mass and oxidative damage were performed. Seeds as well as juvenile plants were used to explore the resilience level of the transgenic plants when compared with wild-type plants.

KEY RESULTS

Overexpression of Scdr2 enhanced germination rates in tobacco seeds under drought and salinity conditions. Juvenile transgenic plants overexpressing Scdr2 and subjected to drought and salinity stresses showed higher photosynthesis levels, internal CO2 concentration and stomatal conductance, reduced accumulation of hydrogen peroxide in the leaves, no penalty for photosystem II and faster recovery after submission to both stress conditions. Respiration was not strongly affected by both stresses in the Scdr2 transgenic plants, whereas wild-type plants exhibited increased respiration rates.

CONCLUSIONS

Scdr2 is involved in the response mechanism to abiotic stresses. Higher levels of Scdr2 enhanced resilience to salinity and drought, and this protection correlated with reduced oxidative damage. Scdr2 confers, at the physiological level, advantages to climate limitations. Therefore, Scdr2 is a potential target for improving sugarcane resilience to abiotic stress.

The C-Terminal Domains SnRK2 Box and ABA Box Have a Role in Sugarcane SnRK2s Auto-Activation and Activity

Resistance to drought stress is fundamental to plant survival and development. Abscisic acid (ABA) is one of the major hormones involved in different types of abiotic and biotic stress responses. ABA intracellular signaling has been extensively explored in Arabidopsis thaliana and occurs via a phosphorylation cascade mediated by three related protein kinases, denominated SnRK2s (SNF1-related protein kinases). However, the role of ABA signaling and the biochemistry of SnRK2 in crop plants remains underexplored. Considering the importance of the ABA hormone in abiotic stress tolerance, here we investigated the regulatory mechanism of sugarcane SnRK2s-known as stress/ABA-activated protein kinases (SAPKs). The crystal structure of ScSAPK10 revealed the characteristic SnRK2 family architecture, in which the regulatory SnRK2 box interacts with the kinase domain αC helix. To study sugarcane SnRK2 regulation, we produced a series of mutants for the protein regulatory domains SnRK2 box and ABA box. Mutations in ScSAPK8 SnRK2 box aimed at perturbing its interaction with the protein kinase domain reduced protein kinase activity in vitro. On the other hand, mutations to ScSAPK ABA box did not impact protein kinase activity but did alter the protein autophosphorylation pattern. Taken together, our results demonstrate that both SnRK2 and ABA boxes might play a role in sugarcane SnRK2 function.

Sugarcane Omics: An Update on the Current Status of Research and Crop Improvement

Sugarcane is an important crop from Poaceae family, contributing about 80% of the total world's sucrose with an annual value of around US$150 billion. In addition, sugarcane is utilized as a raw material for the production of bioethanol, which is an alternate source of renewable energy. Moving towards sugarcane omics, a remarkable success has been achieved in gene transfer from a wide variety of plant and non-plant sources to sugarcane, with the accessibility of efficient transformation systems, selectable marker genes, and genetic engineering gears. Genetic engineering techniques make possible to clone and characterize useful genes and also to improve commercially important traits in elite sugarcane clones that subsequently lead to the development of an ideal cultivar. Sugarcane is a complex polyploidy crop, and hence no single technique has been found to be the best for the confirmation of polygenic and phenotypic characteristics. To better understand the application of basic omics in sugarcane regarding agronomic characters and industrial quality traits as well as responses to diverse biotic and abiotic stresses, it is important to explore the physiology, genome structure, functional integrity, and collinearity of sugarcane with other more or less similar crops/plants. Genetic improvements in this crop are hampered by its complex genome, low fertility ratio, longer production cycle, and susceptibility to several biotic and abiotic stresses. Biotechnology interventions are expected to pave the way for addressing these obstacles and improving sugarcane crop. Thus, this review article highlights up to date information with respect to how advanced data of omics (genomics, transcriptomic, proteomics and metabolomics) can be employed to improve sugarcane crops.

Label-Free Quantitative Proteomics of Enriched Nuclei from Sugarcane (Saccharum ssp) Stems in Response to Drought Stress

Drought is considered the major abiotic stress limiting crop productivity. This study seeks to identify proteins involved in the drought response in sugarcane stems submitted to drought stress. The integration of nuclei enrichment sample preparation with the shotgun proteomic approach results in great coverage of the sugarcane stem proteome with 5381 protein groups identified. A total of 1204 differentially accumulated proteins are detected in response to drought, among which 586 and 618 are increased and reduced in abundance, respectively. A total of 115 exclusive proteins are detected, being 41 exclusives of drought-stressed plants and 74 exclusives of control plants. In the control plants, most of these proteins are related to cell wall metabolism, indicating that drought affects negatively the cell wall metabolism. Also, 37 transcription factors (TFs) are identified, which are low abundant nuclear proteins and are differentially accumulated in response to drought stress. These TFs are associated to protein domains such as leucine-rich (bZIP), C2H2, NAC, C3H, LIM, Myb-related, heat shock factor (HSF) and auxin response factor (ARF). Increased abundance of chromatin remodeling and RNA processing proteins are also observed. It is suggested that these variations result from an imbalance of protein synthesis and degradation processes induced by drought.

Gene Duplication in the Sugarcane Genome: A Case Study of Allele Interactions and Evolutionary Patterns in Two Genic Regions

Sugarcane (Saccharum spp.) is highly polyploid and aneuploid. Modern cultivars are derived from hybridization between S. officinarum and S. spontaneum. This combination results in a genome exhibiting variable ploidy among different loci, a huge genome size (~10 Gb) and a high content of repetitive regions. An approach using genomic, transcriptomic, and genetic mapping can improve our knowledge of the behavior of genetics in sugarcane. The hypothetical HP600 and Centromere Protein C (CENP-C) genes from sugarcane were used to elucidate the allelic expression and genomic and genetic behaviors of this complex polyploid. The physically linked side-by-side genes HP600 and CENP-C were found in two different homeologous chromosome groups with ploidies of eight and ten. The first region (Region01) was a Sorghum bicolor ortholog region with all haplotypes of HP600 and CENP-C expressed, but HP600 exhibited an unbalanced haplotype expression. The second region (Region02) was a scrambled sugarcane sequence formed from different noncollinear genes containing partial duplications of HP600 and CENP-C (paralogs). This duplication resulted in a non-expressed HP600 pseudogene and a recombined fusion version of CENP-C and the orthologous gene Sobic.003G299500 with at least two chimeric gene haplotypes expressed. It was also determined that it occurred before Saccharum genus formation and after the separation of sorghum and sugarcane. A linkage map was constructed using markers from nonduplicated Region01 and for the duplication (Region01 and Region02). We compare the physical and linkage maps, demonstrating the possibility of mapping markers located in duplicated regions with markers in nonduplicated region. Our results contribute directly to the improvement of linkage mapping in complex polyploids and improve the integration of physical and genetic data for sugarcane breeding programs. Thus, we describe the complexity involved in sugarcane genetics and genomics and allelic dynamics, which can be useful for understanding complex polyploid genomes.

Cell Wall Proteome of Sugarcane Young and Mature Leaves and Stems

By characterizing the cell wall proteomes of different sugarcane organs (leaves and stems) at two developmental stages (young vs mature/apical vs basal), it is possible to address unique characteristics in each of them. Four-month-old leaves show a higher proportion of oxido-reductases and proteins related to lipid metabolism (LM), besides a lower proportion of proteins acting on polysaccharides, in comparison to 4-month-old internodes. It is possible to note that sugarcane leaves and young stems have the highest LM rate than all species, which is assumed to be linked to cuticle formation. The data generated enrich the number of cell wall proteins (CWPs) identified in sugarcane, reaching 277. To our knowledge, sugarcane has now the second higher coverage of monocot CWP in plants.

Integration of sugarcane production technologies for enhanced cane and sugar productivity targeting to increase farmers' income: strategies and prospects

The idea of doubling the farmers' income in next 5 years has been slated by the Government of India. The specific target of increasing sugarcane farmers' income could be achieved by developing cost-effective technologies, transferring them from laboratory to land, educating the farmers and creating a linkage between all stakeholders. Consistent efforts shall be required to harness all possible sources for increasing farmer's income in and outside the agriculture sector with respect to improvement in sugarcane and sugar productivity, enhancement in resource use efficiency and adopting various other ways and means including intercropping, management of pests and diseases, use of biotechnological tools and minimizing post-harvest deterioration. The advances in sugarcane biotechnology could become remarkable in the coming years, both in terms of improving productivity as well as increasing the value and utility of this crop substantially. In future, genetically modified sugarcane varieties with increased resistance to different biotic and abiotic stresses would serve more towards sugarcane crop improvement. Any possibility of enhancement in the income of sugarcane farmers shall also be dependent upon the profitability and sustainability of the sugar industry. Integration of sugarcane production technologies for improvement in farm productivity, diversified sugarcane production system, reduced cost of cultivation along with increased processing plant efficiency and diversification to produce value added products shall ensure smooth and higher payment to the farmers. Development of low-cost technologies to convert "waste to resource" on a smaller scale shall also help the farmers to increase their income further. This paper focuses on possible measures to be taken up in each aspects of sugarcane cultivation including biotechnological approaches to achieve the goal of enhancing the income of sugarcane farmers substantially, particularly in the sub-tropical region of India.

A miniature inverted-repeat transposable element, AddIn-MITE, located inside a WD40 gene is conserved in Andropogoneae grasses

Miniature inverted-repeat transposable elements (MITEs) have been associated with genic regions in plant genomes and may play important roles in the regulation of nearby genes via recruitment of small RNAs (sRNA) to the MITEs loci. We identified eight families of MITEs in the sugarcane genome assembly with MITE-Hunter pipeline. These sequences were found to be upstream, downstream or inserted into 67 genic regions in the genome. The position of the most abundant MITE (Stowaway-like) in genic regions, which we call AddIn-MITE, was confirmed in a WD40 gene. The analysis of four monocot species showed conservation of the AddIn-MITE sequence, with a large number of copies in their genomes. We also investigated the conservation of the AddIn-MITE’ position in the WD40 genes from sorghum, maize and, in sugarcane cultivars and wild Saccharum species. In all analyzed plants, AddIn-MITE has located in WD40 intronic region. Furthermore, the role of AddIn-MITE-related sRNA in WD40 genic region was investigated. We found sRNAs preferentially mapped to the AddIn-MITE than to other regions in the WD40 gene in sugarcane. In addition, the analysis of the small RNA distribution patterns in the WD40 gene and the structure of AddIn-MITE, suggests that the MITE region is a proto-miRNA locus in sugarcane. Together, these data provide insights into the AddIn-MITE role in Andropogoneae grasses.

Network Analyses and Data Integration of Proteomics and Metabolomics From Leaves of Two Contrasting Varieties of Sugarcane in Response to Drought

Uncovering the molecular mechanisms involved in the responses of crops to drought is crucial to understand and enhance drought tolerance mechanisms. Sugarcane (Saccharum spp.) is an important commercial crop cultivated mainly in tropical and subtropical areas for sucrose and ethanol production. Usually, drought tolerance has been investigated by single omics analysis (e.g. global transcripts identification). Here we combine label-free quantitative proteomics and metabolomics data (GC-TOF-MS), using a network-based approach, to understand how two contrasting commercial varieties of sugarcane, CTC15 (tolerant) and SP90-3414 (susceptible), adjust their leaf metabolism in response to drought. To this aim, we propose the utilization of regularized canonical correlation analysis (rCCA), which is a modification of classical CCA, and explores the linear relationships between two datasets of quantitative variables from the same experimental units, with a threshold set to 0.99. Light curves revealed that after 4 days of drought, the susceptible variety had its photosynthetic capacity already significantly reduced, while the tolerant variety did not show major reduction. Upon 12 days of drought, photosynthesis in the susceptible plants was completely reduced, while the tolerant variety was at a third of its rate under control conditions. Network analysis of proteins and metabolites revealed that different biological process had a stronger impact in each variety (e.g. translation in CTC15, generation of precursor metabolites, response to stress and energy in SP90-3414). Our results provide a reference data set and demonstrate that rCCA can be a powerful tool to infer experimentally metabolite-protein or protein-metabolite associations to understand plant biology.

Overexpression of ScMYBAS1 alternative splicing transcripts differentially impacts biomass accumulation and drought tolerance in rice transgenic plants

Drought is the most significant environmental stress for agricultural production worldwide, and tremendous efforts have been made to improve crop yield under the increasing water scarcity. Transcription factors are major players in the regulation of water stress-related genes in plants. Recently, different MYB transcription factors were characterized for their involvement in drought response. A sugarcane R2R3-MYB gene (ScMYBAS1) and its four alternative forms of transcript (ScMYAS1-2, ScMYBAS1-3, ScMYBAS1-4 and ScMYBAS1-5) were identified in this study. The subcellular localization, in Nicotiniana benthamiana, of the TFs fused in frame with GFP revealed that ScMYBAS1-2-GFP and ScMYBAS1-3-GFP were observed in the nucleus. The overexpression of ScMYBAS1-2 and ScMYBAS1-3 spliced transcripts in rice promoted change in plant growth under both well-watered and drought conditions. The ScMYBAS1-2 and ScMYBAS1-3 transgenic lines revealed a higher relative water content (RWC) compared to the wild type before maximum stress under drought conditions. The ScMYBAS1-2 transgenic lines showed a reduction in biomass (total dry weight). Conversely, ScMYBAS1-3 showed an increased biomass (total dry weight) relative to the wild-type. The overexpression of ScMYBAS1-3 in rice transgenic lines showed involvement with drought tolerance and biomass and, for this reason, was considered a good target for plant transformation, particularly for use in developing genotypes with drought tolerance and biomass accumulation.

Novel Insights Into the Early Stages of Ratoon Stunting Disease of Sugarcane Inferred from Transcript and Protein Analysis

Despite of the importance of ratoon stunting disease, little is known on the responses of sugarcane to its causal agent, the vascular bacterial endophyte Leifsonia xyli subsp. xyli. The transcriptome and proteome of young plants of a susceptible cultivar with no symptoms of stunting but with relative low and high bacterial titers were compared at 30 and 60 days after inoculation. Increased bacterial titers were associated with alterations in the expression of 267 cDNAs and in the abundance of 150 proteins involved in plant growth, hormone metabolism, signal transduction and defense responses. Some alterations are predicted to benefit the pathogen, such as the up-regulation of genes involved in the synthesis of methionine. Also, genes and proteins of the cell division cycle were all down-regulated in plants with higher titers at both times. It is hypothesized that the negative effects on cell division related to increased bacterial titers is cumulative over time and its modulation by other host and environmental factors results in the stunting symptom.

A mosaic monoploid reference sequence for the highly complex genome of sugarcane

Sugarcane (Saccharum spp.) is a major crop for sugar and bioenergy production. Its highly polyploid, aneuploid, heterozygous, and interspecific genome poses major challenges for producing a reference sequence. We exploited colinearity with sorghum to produce a BAC-based monoploid genome sequence of sugarcane. A minimum tiling path of 4660 sugarcane BAC that best covers the gene-rich part of the sorghum genome was selected based on whole-genome profiling, sequenced, and assembled in a 382-Mb single tiling path of a high-quality sequence. A total of 25,316 protein-coding gene models are predicted, 17% of which display no colinearity with their sorghum orthologs. We show that the two species, S. officinarum and S. spontaneum, involved in modern cultivars differ by their transposable elements and by a few large chromosomal rearrangements, explaining their distinct genome size and distinct basic chromosome numbers while also suggesting that polyploidization arose in both lineages after their divergence.

Characterization of leaf transcriptome, development and utilization of unigenes-derived microsatellite markers in sugarcane (Saccharum sp. hybrid)

Sugarcane (Saccharum species hybrid) is the major source of sugar (> 80% sugar) in the world and is cultivated in more than 115 countries. It has recently gained attention as a source of biofuel (ethanol). Due to genomic complexity, the development of new genomic resources is imperative in understanding the gene regulation and function, and to fine tune the genetic improvement of sugarcane. In this study, a cDNA library was constructed from mature leaves so as to develop ESTs resources which were further compared with nucleotide and protein databases to explore the functional identity of sugarcane genes. The non-redundant ESTs (unigenes) were categorized into 18 metabolic functions. The major categories were bioenergetics and photosynthesis (4%), cell metabolism (5%), development related protein (3%), membrane-related, mobile genetic elements (5%), signal transduction (2%), DNA (1%), RNA (1%) and protein (2%) metabolism, other metabolic processes (3%), transcription factors (1%), transport (4%) and proteins related to stress/defense (4%). From 540 unique ESTs, 212 simple sequence repeats were identified, of which 206 were from 463 singlets and six were mined from 77 contig sequences. A total of 540 unique EST sequences were used for SSR search of which 97 (17.9%) contained specified SSR motifs, generating 212 unique SSRs. The genes characterized in this study and the EST-derived microsatellite markers identified from the cDNA library will enrich genomic resources for association- and linkage-mapping studies in sugarcane.

Detection and validation of EST-SSR markers associated with sugar-related traits in sugarcane using linkage and association mapping

Sugar-related traits are of great importance in sugarcane breeding. In the present study, quantitative trait loci (QTL) mapping validated with association mapping was used to identify expressed sequence tag-simple sequence repeats (EST-SSRs) associated with sugar-related traits. For linkage mapping, 524 EST-SSRs, 241 Amplified Fragment Length Polymorphisms, and 10 genomic SSR markers were mapped using 283 F₁ progenies derived from an interspecific cross. Six regions were identified using Multiple QTL Mapping, and 14 unlinked markers using single marker analysis. Association analysis was performed on a set of 200 accessions, based on the mixed linear model. Validation of the EST-SSR markers using association mapping within the target QTL genomic regions identified two EST-SSR markers showing a putative relationship with uridine diphosphate (UDP) glycosyltransferase, and beta-amylase, which are associated with pol and sugar yield. These functional markers can be used for marker-assisted selection of 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.

Potentials, Challenges, and Genetic and Genomic Resources for Sugarcane Biomass Improvement

Lignocellulosic biomass has become an emerging feedstock for second-generation bioethanol production. Sugarcane (Saccharum spp. hybrids), a very efficient perennial C4 plant with a high polyploid level and complex genome, is considered a top-notch candidate for biomass production due to its salient features viz. fast growth rate and abilities for high tillering, ratooning, and photosynthesis. Energy cane, an ideal type of sugarcane, has been bred specifically as a biomass crop. In this review, we described (1) biomass potentials of sugarcane and its underlying genetics, (2) challenges associated with biomass improvement such as large and complex genome, narrow gene pool in existing commercial cultivars, long breeding cycle, and non-synchronous flowering, (3) available genetic resources such as germplasm resources, and genomic and cell wall-related databases that facilitate biomass improvement, and (4) mining candidate genes controlling biomass in genomic databases. We extensively reviewed databases for biomass-related genes and their usefulness in biofuel generation. This review provides valuable resources for sugarcane breeders, geneticists, and broad scientific communities involved in bioenergy production.

Association of variation in the sugarcane transcriptome with sugar content

Background

Sugarcane is a major crop of the tropics cultivated mainly for its high sucrose content. The crop is genetically less explored due to its complex polyploid genome. Sucrose synthesis and accumulation are complex processes influenced by physiological, biochemical and genetic factors, and the growth environment. The recent focus on the crop for fibre and biofuel has led to a renewed interest on understanding the molecular basis of sucrose and biomass traits. This transcriptome study aimed to identify genes that are associated with and differentially regulated during sucrose synthesis and accumulation in the mature stage of sugarcane. Patterns of gene expression in high and low sugar genotypes as well as mature and immature culm tissues were studied using RNA-Seq of culm transcriptomes.

Results

In this study, 28 RNA-Seq libraries from 14 genotypes of sugarcane differing in their sucrose content were used for studying the transcriptional basis of sucrose accumulation. Differential gene expression studies were performed using SoGI (Saccharum officinarum Gene Index, 3.0), SAS (sugarcane assembled sequences) of sugarcane EST database (SUCEST) and SUGIT, a sugarcane Iso-Seq transcriptome database. In total, about 34,476 genes were found to be differentially expressed between high and low sugar genotypes with the SoGI database, 20,487 genes with the SAS database and 18,543 genes with the SUGIT database at FDR < 0.01, using the Baggerley’s test. Further, differential gene expression analyses were conducted between immature (top) and mature (bottom) tissues of the culm. The DEGs were functionally annotated using GO classification and the genes consistently associated with sucrose accumulation were identified.

Conclusions

The large number of DEGs may be due to the large number of genes that influence sucrose content or are regulated by sucrose content. These results indicate that apart from being a primary metabolite and storage and transport sugar, sucrose may serve as a signalling molecule that regulates many aspects of growth and development in sugarcane. Further studies are needed to confirm if sucrose regulates the expression of the identified DEGs or vice versa. The DEGs identified in this study may lead to identification of genes/pathways regulating sucrose accumulation and/or regulated by sucrose levels in sugarcane. We propose identifying the master regulators of sucrose if any in the future.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4302-5) contains supplementary material, which is available to authorized users.

Identification, classification and transcriptional profiles of dirigent domain-containing proteins in sugarcane

Dirigent (DIR) proteins, encoded by DIR genes, are referred to as "dirigent" because they direct the outcome of the coupling of the monolignol coniferyl alcohol into (+) or (-) pinoresinol, the first intermediates in the enantiocomplementary pathways for lignan biosynthesis. DIR domain-containing or DIR-like proteins are, thus, termed for not having a clear characterization. A transcriptome- and genome-wide survey of DIR domain-containing proteins in sugarcane was carried out, in addition to phylogenetic, physicochemical and transcriptional analyses. A total of 120 non-redundant sequences containing the DIR domain were identified and classified into 64 groups according to phylogenetic and sequence alignment analyses. In silico analysis of transcript abundance showed that these sequences are expressed at low levels in leaves and genes in the same phylogenetic clade have similar expression patterns. Expression analysis of ShDIR1-like transcripts in the culm internodes of sugarcane demonstrates their abundance in mature internodes, their induction by nitrogen fertilization and their predominant expression in cells that have a lignified secondary cell wall, such as vascular bundles of young internodes and parenchymal cells of the pith of mature internodes. Due to the lack of information about the functional role of DIR in plants, a possible relationship is discussed between the ShDIR1-like transcriptional profile and cell wall development in parenchyma cells of sugarcane culm, which typically accumulates large amounts of sucrose. The number of genes encoding the DIR domain-containing proteins in sugarcane is intriguing and is an indication per se that these proteins may have an important metabolic role and thus deserve to be better studied.

Sugarcane Water Stress Tolerance Mechanisms and Its Implications on Developing Biotechnology Solutions

Sugarcane is a unique crop with the ability to accumulate high levels of sugar and is a commercially viable source of biomass for bioelectricity and second-generation bioethanol. Water deficit is the single largest abiotic stress affecting sugarcane productivity and the development of water use efficient and drought tolerant cultivars is an imperative for all major sugarcane producing countries. This review summarizes the physiological and molecular studies on water deficit stress in sugarcane, with the aim to help formulate more effective research strategies for advancing our knowledge on genes and mechanisms underpinning plant response to water stress. We also overview transgenic studies in sugarcane, with an emphasis on the potential strategies to develop superior sugarcane varieties that improve crop productivity in drought-prone environments.

A New Sugarcane Cystatin Strongly Binds to Dental Enamel and Reduces Erosion

Cystatin B was recently identified as an acid-resistant protein in acquired enamel pellicle; it could therefore be included in oral products to protect against caries and erosion. However, human recombinant cystatin is very expensive, and alternatives to its use are necessary. Phytocystatins are reversible inhibitors of cysteine peptidases that are found naturally in plants. In plants, they have several biological and physiological functions, such as the regulation of endogenous processes, defense against pathogens, and response to abiotic stress. Previous studies performed by our research group have reported high inhibitory activity and potential agricultural and medical applications of several sugarcane cystatins, including CaneCPI-1, CaneCPI-2, CaneCPI-3, and CaneCPI-4. In the present study, we report the characterization of a novel sugarcane cystatin, named CaneCPI-5. This cystatin was efficiently expressed in Escherichia coli, and inhibitory assays demonstrated that it was a potent inhibitor of human cathepsins B, K, and L ( Ki = 6.87, 0.49, and 0.34 nM, respectively). The ability of CaneCPI-5 to bind to dental enamel was evaluated using atomic force microscopy. Its capacity to protect against initial enamel erosion was also tested in vitro via changes in surface hardness. CaneCPI-5 showed a very large force of interaction with enamel (e.g., compared with mucin and casein) and significantly reduced initial enamel erosion. These results suggest that the inclusion of CaneCPIs in dental products might confer protection against enamel erosion.

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.

Ethylene-induced transcriptional and hormonal responses at the onset of sugarcane ripening

The effects of ethephon as a sugarcane ripener are attributed to ethylene. However, the role of this phytohormone at the molecular level is unknown. We performed a transcriptome analysis combined with the evaluation of sucrose metabolism and hormone profiling of sugarcane plants sprayed with ethephon or aminoethoxyvinylglycine (AVG), an ethylene inhibitor, at the onset of ripening. The differential response between ethephon and AVG on sucrose level and sucrose synthase activity in internodes indicates ethylene as a potential regulator of sink strength. The correlation between hormone levels and transcriptional changes suggests ethylene as a trigger of multiple hormone signal cascades, with approximately 18% of differentially expressed genes involved in hormone biosynthesis, metabolism, signalling, and response. A defence response elicited in leaves favoured salicylic acid over the ethylene/jasmonic acid pathway, while the upper internode was prone to respond to ethylene with strong stimuli on ethylene biosynthesis and signalling genes. Besides, ethylene acted synergistically with abscisic acid, another ripening factor, and antagonistically with gibberellin and auxin. We identified potential ethylene target genes and characterized the hormonal status during ripening, providing insights into the action of ethylene at the site of sucrose accumulation. A molecular model of ethylene interplay with other hormones is proposed.

GBS-based single dosage markers for linkage and QTL mapping allow gene mining for yield-related traits in sugarcane

BACKGROUND

Sugarcane (Saccharum spp.) is predominantly an autopolyploid plant with a variable ploidy level, frequent aneuploidy and a large genome that hampers investigation of its organization. Genetic architecture studies are important for identifying genomic regions associated with traits of interest. However, due to the genetic complexity of sugarcane, the practical applications of genomic tools have been notably delayed in this crop, in contrast to other crops that have already advanced to marker-assisted selection (MAS) and genomic selection. High-throughput next-generation sequencing (NGS) technologies have opened new opportunities for discovering molecular markers, especially single nucleotide polymorphisms (SNPs) and insertion-deletion (indels), at the genome-wide level. The objectives of this study were to (i) establish a pipeline for identifying variants from genotyping-by-sequencing (GBS) data in sugarcane, (ii) construct an integrated genetic map with GBS-based markers plus target region amplification polymorphisms and microsatellites, (iii) detect QTLs related to yield component traits, and (iv) perform annotation of the sequences that originated the associated markers with mapped QTLs to search putative candidate genes.

RESULTS

We used four pseudo-references to align the GBS reads. Depending on the reference, from 3,433 to 15,906 high-quality markers were discovered, and half of them segregated as single-dose markers (SDMs) on average. In addition to 7,049 non-redundant SDMs from GBS, 629 gel-based markers were used in a subsequent linkage analysis. Of 7,678 SDMs, 993 were mapped. These markers were distributed throughout 223 linkage groups, which were clustered in 18 homo(eo)logous groups (HGs), with a cumulative map length of 3,682.04 cM and an average marker density of 3.70 cM. We performed QTL mapping of four traits and found seven QTLs. Our results suggest the presence of a stable QTL across locations. Furthermore, QTLs to soluble solid content (BRIX) and fiber content (FIB) traits had markers linked to putative candidate genes.

CONCLUSIONS

This study is the first to report the use of GBS for large-scale variant discovery and genotyping of a mapping population in sugarcane, providing several insights regarding the use of NGS data in a polyploid, non-model species. The use of GBS generated a large number of markers and still enabled ploidy and allelic dosage estimation. Moreover, we were able to identify seven QTLs, two of which had great potential for validation and future use for molecular breeding in sugarcane.