The Challenge of Analyzing the Sugarcane Genome

Identification and expression analysis of nuclear factor Y transcription factor genes under drought, cold and Eldana infestation in sugarcane (Saccharum spp. hybrid)

BACKGROUND

The Nuclear Factor Y (NF-Y) transcription factor (TF) gene family plays a crucial role in plant development and response to stress. Limited information is available on this gene family in sugarcane.

OBJECTIVES

To identify sugarcane NF-Y genes through bioinformatic analysis and phylogenetic association and investigate the expression of these genes in response to abiotic and biotic stress.

METHODS

Sugarcane NF-Y genes were identified using comparative genomics from functionally annotated Poaceae and Arabidopsis species. Quantitative PCR and transcriptome analysis assigned preliminary functional roles to these genes in response to water deficit, cold and African sugarcane borer (Eldana saccharina) infestation.

RESULTS

We identify 21 NF-Y genes in sugarcane. Phylogenetic analysis revealed three main branches representing the subunits with potential discrepancies present in the assignment of numerical names of some NF-Y putative orthologs across the different species. Gene expression analysis indicated that three genes, ShNF-YA1, A3 and B3 were upregulated and two genes, NF-YA4 and A7 were downregulated, while three genes were upregulated, ShNF-YB2, B3 and C4, in the plants exposed to water deficit and cold stress, respectively. Functional involvement of NF-Y genes in the biotic stress response were also detected where three genes, ShNF-YA6, A3 and A7 were downregulated in the early resistant (cv. N33) response to Eldana infestation whilst only ShNF-YA6 was downregulated in the susceptible (cv. N11) early response.

CONCLUSIONS

Our research findings establish a foundation for investigating the function of ShNF-Ys and offer candidate genes for stress-resistant breeding and improvement in sugarcane.

A Comprehensive Analysis of Transcriptomics and Metabolomics Revealed Key Pathways Involved in Saccharum spontaneum Defense against Sporisorium scitamineum

Sugarcane smut, caused by Sporisorium scitamineum, poses a severe threat to sugarcane production. The genetic basis of sugarcane resistance to S. scitamineum remains elusive. A comparative transcriptomic and metabolomic study was conducted on two wild Saccharum species of S. spontaneum with contrast smut resistance. Following infection, the resistant line exhibited greater down-regulation of genes and metabolites compared to the susceptible line, indicating distinct biological processes. Lignan and lignin biosynthesis and SA signal transduction were activated in the resistant line, while flavonoid biosynthesis and auxin signal transduction were enhanced in the susceptible line. TGA2.2 and ARF14 were identified as playing positive and negative roles, respectively, in plant defense. Exogenous auxin application significantly increased the susceptibility of S. spontaneum to S. scitaminum. This study established the significant switching of defense signaling pathways in contrast-resistant S. spontaneum following S. scitamineum infection, offering a hypothetical model and candidate genes for further research into sugarcane smut disease.

A SNP variation in the Sucrose synthase (SoSUS) gene associated with sugar-related traits in sugarcane

Background

Sugarcane (Saccharum spp.) is an economically significant crop for both the sugar and biofuel industries. Breeding sugarcane cultivars with high-performance agronomic traits is the most effective approach for meeting the rising demand for sugar and biofuels. Molecular markers associated with relevant agronomic traits could drastically reduce the time and resources required to develop new sugarcane varieties. Previous sugarcane candidate gene association analyses have found single nucleotide polymorphism (SNP) markers associated with sugar-related traits. This study aims to validate these associated SNP markers of six genes, including Lesion simulating disease 1 (LSD), Calreticulin (CALR), Sucrose synthase 1 (SUS1), DEAD-box ATP-dependent RNA helicase (RH), KANADI1 (KAN1), and Sodium/hydrogen exchanger 7 (NHX7), in a diverse population in 2-year and two-location evaluations.

Methods

After genotyping of seven targeted SNP markers was performed by PCR Allelic Competitive Extension (PACE) SNP genotyping, the association with sugar-related traits and important cane yield component traits was determined on a set of 159 sugarcane genotypes. The marker-trait relationships were validated and identified by both t-test analysis and an association analysis based on the general linear model.

Results

The mSoSUS1_SNPCh10.T/C and mSoKAN1_SNPCh7.T/C markers that were designed from the SUS1 and KAN1 genes, respectively, showed significant associations with different amounts of sugar-related traits and yield components. The mSoSUS1_SNPCh10.T/C marker was found to have more significant association with sugar-related traits, including pol, CCS, brix, fiber and sugar yield, with p values of 6.08 × 10-6 to 4.35 × 10-2, as well as some cane yield component traits with p values of 1.61 × 10-4 to 3.35 × 10-2. The significant association is consistent across four environments.

Conclusion

Sucrose synthase (SUS) is considered a crucial enzyme involved in sucrose metabolism. This marker is a high potential functional marker that may be used in sugarcane breeding programs to select superior sugarcane with good fiber and high sugar contents.

Genomic prediction with machine learning in sugarcane, a complex highly polyploid clonally propagated crop with substantial non-additive variation for key traits

Sugarcane has a complex, highly polyploid genome with multi-species ancestry. Additive models for genomic prediction of clonal performance might not capture interactions between genes and alleles from different ploidies and ancestral species. As such, genomic prediction in sugarcane presents an interesting case for machine learning (ML) methods, which are purportedly able to deal with high levels of complexity in prediction. Here, we investigated deep learning (DL) neural networks, including multilayer networks (MLP) and convolution neural networks (CNN), and an ensemble machine learning approach, random forest (RF), for genomic prediction in sugarcane. The data set used was 2912 sugarcane clones, scored for 26,086 genome wide single nucleotide polymorphism markers, with final assessment trial data for total cane harvested (TCH), commercial cane sugar (CCS), and fiber content (Fiber). The clones in the latest trial (2017) were used as a validation set. We compared prediction accuracy of these methods to genomic best linear unbiased prediction (GBLUP) extended to include dominance and epistatic effects. The prediction accuracies from GBLUP models were up to 0.37 for TCH, 0.43 for CCS, and 0.48 for Fiber, while the optimized ML models had prediction accuracies of 0.35 for TCH, 0.38 for CCS, and 0.48 for Fiber. Both RF and DL neural network models have comparable predictive ability with the additive GBLUP model but are less accurate than the extended GBLUP model.

Comparative expression analysis of sucrose phosphate synthase gene family in a low and high sucrose Pakistani sugarcane cultivars

Sugarcane is the world's largest cultivated crop by biomass and is the main source of sugar and biofuel. Sucrose phosphate synthase (SPS) enzymes are directly involved in the synthesis of sucrose. Here, we analyzed and compared one of the important gene families involved in sucrose metabolism in a high and low sucrose sugarcane cultivar. A comprehensive in silico analysis of the SoSPS family displayed their phylogenetic relationship, gene and protein structure, miRNA targets, protein interaction network (PPI), gene ontology and collinearity. This was followed by a spatial expression analysis in two different sugarcane varieties. The phylogenetic reconstruction distributed AtSPS, ZmSPS, OsSPS, SoSPS and SbSPS into three main groups (A, B, C). The regulatory region of SoSPS genes carries ABRE, ARE, G-box, and MYC as the most dominant cis-regulatory elements. The PPI analysis predicted a total of 14 unique proteins interacting with SPS. The predominant expression of SPS in chloroplast clearly indicates that they are the most active in the organelle which is the hub of photosynthesis. Similarly, gene ontology attributed SPS to sucrose phosphate synthase and glucosyl transferase molecular functions, as well as sucrose biosynthetic and disaccharide biological processes. Overall, the expression of SPS in CPF252 (high sucrose variety) was higher in leaf and culm as compared to that of CPF 251 (low sucrose variety). In brief, this study adds to the present literature about sugarcane, sucrose metabolism and role of SPS in sucrose metabolism thereby opening up further avenues of research in crop improvement.

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.

Evaluation of sugarcane promising clones based on the morphophysiological traits developed from fuzz

Sugarcane is one of the critical commercial crops and principal sources of ethanol and sugar worldwide. Unfavorable conditions and poor seed setting rates hinder variety development in sugarcane. Countries like Pakistan directly import fuzz (true seed) and other propagation material from the USA, China, Brazil, etc. In this study, we imported fuzz from China, developed 29 genotypes germinating in the glasshouse, and evaluated at field conditions along with two local checks (CPF-251 and HSF-240). Morphophysiological data were recorded, including plant height (PH), cane length (CL), internodal length (IL), tiller number (TN), brix percentage (B), cane diameter (CD), chlorophyll a (Chl. a), chlorophyll b (Chl. b), and total chlorophyll (T. Chl). Results showed highly significant (p < 0.001) differences among the sugarcane accessions for all the studied traits. High broad-sense heritability (81.89% to 99.91%) was recorded for all the studied parameters. Genetic Advance (GA) ranges from 4.6% to 65.32%. The highest GA was observed for PH (65.32%), followed by CL (63.28%). Chlorophyll leaching assay was also performed at different time points (0, 50, 100, 150, and 200 min). All the genotypes showed the same leaching trend at all times, and better performing genotypes showed less leaching compared to poor performing, indicating the high amount of cutin and wax on the leaf surface. Correlation analysis showed that PH, CL, IL, and TN had significant associations. Principal components analysis (PCA) further confirms these results. Based on PCA and correlation results, PH, CL, IL, and TN can be utilized as a selection criterion for sugarcane improvement. Genotypes such as NS-4a, NS-5, NS-6, NS-8, NS-9, and NS-15 are recommended for future breeding programs related to sugarcane variety development.

Recent Advances in Assembly of Complex Plant Genomes

Over the past 20 years, tremendous advances in sequencing technologies and computational algorithms have spurred plant genomic research into a thriving era with hundreds of genomes decoded already, ranging from those of nonvascular plants to those of flowering plants. However, complex plant genome assembly is still challenging and remains difficult to fully resolve with conventional sequencing and assembly methods due to high heterozygosity, highly repetitive sequences, or high ploidy characteristics of complex genomes. Herein, we summarize the challenges of and advances in complex plant genome assembly, including feasible experimental strategies, upgrades to sequencing technology, existing assembly methods, and different phasing algorithms. Moreover, we list actual cases of complex genome projects for readers to refer to and draw upon to solve future problems related to complex genomes. Finally, we expect that the accurate, gapless, telomere-to-telomere, and fully phased assembly of complex plant genomes could soon become routine.

A complete gap-free diploid genome in Saccharum complex and the genomic footprints of evolution in the highly polyploid Saccharum genus

A diploid genome in the Saccharum complex facilitates our understanding of evolution in the highly polyploid Saccharum genus. Here we have generated a complete, gap-free genome assembly of Erianthus rufipilus, a diploid species within the Saccharum complex. The complete assembly revealed that centromere satellite homogenization was accompanied by the insertions of Gypsy retrotransposons, which drove centromere diversification. An overall low rate of gene transcription was observed in the palaeo-duplicated chromosome EruChr05 similar to other grasses, which might be regulated by methylation patterns mediated by homologous 24 nt small RNAs, and potentially mediating the functions of many nucleotide-binding site genes. Sequencing data for 211 accessions in the Saccharum complex indicated that Saccharum probably originated in the trans-Himalayan region from a diploid ancestor (x = 10) around 1.9-2.5 million years ago. Our study provides new insights into the origin and evolution of Saccharum and accelerates translational research in cereal genetics and genomics.

The evolution and expression of stomatal regulators in C3 and C4 crops: Implications on the divergent drought tolerance

Drought stress is a major environmental hazard. Stomatal development is highly responsive to abiotic stress and has been used as a cellular marker for drought-tolerant crop selection. C3 and C4 crops have evolved into different photosynthetic systems and physiological responses to water deficits. The genome sequences of maize, sorghum, and sugarcane make it possible to explore the association of the stomatal response to drought stress with the evolution of the key stomatal regulators. In this study, phylogenic analysis, gene expression analysis and stomatal assay under drought stress were used to investigate the drought tolerance of C3 and C4 plants. Our data shows that C3 and C4 plants exhibit different drought responses at the cellular level. Drought represses the growth and stomatal development of C3 crops but has little effect on that of C4 plants. In addition, stomatal development is unresponsive to drought in drought-tolerant C3 crops but is repressed in drought-tolerant C4 plants. The different developmental responses to drought in C3 and C4 plants might be associated with the divergent expression of their SPEECHLESS genes. In particular, C4 crops have evolved to generate multiple SPEECHLESS homologs with different genetic structure and expression levels. Our research provides not only molecular evidence that supports the evolutionary history of C4 from C3 plants but also a possible molecular model that controls the cellular response to abiotic stress in C3 and C4 crops.

Functional characterization of the sugarcane (Saccharum spp.) ammonium transporter AMT2;1 suggests a role in ammonium root-to-shoot translocation

AMMONIUM TRANSPORTER/METHYLAMMONIUM PERMEASE/RHESUS (AMT) family members transport ammonium across membranes in all life domains. Plant AMTs can be categorized into AMT1 and AMT2 subfamilies. Functional studies of AMTs, particularly AMT1-type, have been conducted using model plants but little is known about the function of AMTs from crops. Sugarcane (Saccharum spp.) is a major bioenergy crop that requires heavy nitrogen fertilization but depends on a low carbon-footprint for competitive sustainability. Here, we identified and functionally characterized sugarcane ScAMT2;1 by complementing ammonium uptake-defective mutants of Saccharomyces cerevisiae and Arabidopsis thaliana. Reporter gene driven by the ScAMT2;1 promoter in A. thaliana revealed preferential expression in the shoot vasculature and root endodermis/pericycle according to nitrogen availability and source. Arabidopsis quadruple mutant plants expressing ScAMT2;1 driven by the CaMV35S promoter or by a sugarcane endogenous promoter produced significantly more biomass than mutant plants when grown in NH₄ ⁺ and showed more ¹⁵N-ammonium uptake by roots and nitrogen translocation to shoots. In A. thaliana, ScAMT2;1 displayed a K_m of 90.17 µM and V_max of 338.99 µmoles h^-1 g^-1 root DW. Altogether, our results suggest that ScAMT2;1 is a functional high-affinity ammonium transporter that might contribute to ammonium uptake and presumably to root-to-shoot translocation under high NH₄ ⁺ conditions.

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.

A joint learning approach for genomic prediction in polyploid grasses

Poaceae, among the most abundant plant families, includes many economically important polyploid species, such as forage grasses and sugarcane (Saccharum spp.). These species have elevated genomic complexities and limited genetic resources, hindering the application of marker-assisted selection strategies. Currently, the most promising approach for increasing genetic gains in plant breeding is genomic selection. However, due to the polyploidy nature of these polyploid species, more accurate models for incorporating genomic selection into breeding schemes are needed. This study aims to develop a machine learning method by using a joint learning approach to predict complex traits from genotypic data. Biparental populations of sugarcane and two species of forage grasses (Urochloa decumbens, Megathyrsus maximus) were genotyped, and several quantitative traits were measured. High-quality markers were used to predict several traits in different cross-validation scenarios. By combining classification and regression strategies, we developed a predictive system with promising results. Compared with traditional genomic prediction methods, the proposed strategy achieved accuracy improvements exceeding 50%. Our results suggest that the developed methodology could be implemented in breeding programs, helping reduce breeding cycles and increase genetic gains.

Taxonomically Restricted Genes Are Associated With Responses to Biotic and Abiotic Stresses in Sugarcane (Saccharum spp.)

Orphan genes (OGs) are protein-coding genes that are restricted to particular clades or species and lack homology with genes from other organisms, making their biological functions difficult to predict. OGs can rapidly originate and become functional; consequently, they may support rapid adaptation to environmental changes. Extensive spread of mobile elements and whole-genome duplication occurred in the Saccharum group, which may have contributed to the origin and diversification of OGs in the sugarcane genome. Here, we identified and characterized OGs in sugarcane, examined their expression profiles across tissues and genotypes, and investigated their regulation under varying conditions. We identified 319 OGs in the Saccharum spontaneum genome without detected homology to protein-coding genes in green plants, except those belonging to Saccharinae. Transcriptomic analysis revealed 288 sugarcane OGs with detectable expression levels in at least one tissue or genotype. We observed similar expression patterns of OGs in sugarcane genotypes originating from the closest geographical locations. We also observed tissue-specific expression of some OGs, possibly indicating a complex regulatory process for maintaining diverse functional activity of these genes across sugarcane tissues and genotypes. Sixty-six OGs were differentially expressed under stress conditions, especially cold and osmotic stresses. Gene co-expression network and functional enrichment analyses suggested that sugarcane OGs are involved in several biological mechanisms, including stimulus response and defence mechanisms. These findings provide a valuable genomic resource for sugarcane researchers, especially those interested in selecting stress-responsive genes.

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.

Transcriptome changes in the developing sugarcane culm associated with high yield and early-season high sugar content

Sugarcane, with its exceptional carbon dioxide assimilation, biomass and sugar yield, has a high potential for the production of bio-energy, bio-plastics and high-value products in the food and pharmaceutical industries. A crucial challenge for long-term economic viability and environmental sustainability is also to optimize the production of biomass composition and carbon sequestration. Sugarcane varieties such as KQ228 and Q253 are highly utilized in the industry. These varieties are characterized by a high early-season sugar content associated with high yield. In order to investigate these correlations, 1,440 internodes were collected and combined to generate a set of 120 samples in triplicate across 24 sugarcane cultivars at five different development stages. Weighted gene co-expression network analysis (WGCNA) was used and revealed for the first time two sets of co-expressed genes with a distinct and opposite correlation between fibre and sugar content. Gene identification and metabolism pathways analysis was used to define these two sets of genes. Correlation analysis identified a large number of interconnected metabolic pathways linked to sugar content and fibre content. Unsupervised hierarchical clustering of gene expression revealed a stronger level of segregation associated with the genotypes than the stage of development, suggesting a dominant genetic influence on biomass composition and facilitating breeding selection. Characterization of these two groups of co-expressed key genes can help to improve breeding program for high fibre, high sugar species or plant synthetic biology.

A Reappraisal of Polyploidy Events in Grasses (Poaceae) in a Rapidly Changing World

Around 80% of megaflora species became extinct at the Cretaceous–Paleogene (K–Pg) boundary. Subsequent polyploidy events drove the survival of thousands of plant species and played a significant historical role in the development of the most successful modern cereal crops. However, current and rapid global temperature change poses an urgent threat to food crops worldwide, including the world’s big three cereals: rice, wheat, and maize, which are members of the grass family, Poaceae. Some minor cereals from the same family (such as teff) have grown in popularity in recent years, but there are important knowledge gaps regarding the similarities and differences between major and minor crops, including how polyploidy affects their biological processes under natural and (a)biotic stress conditions and thus the potential to harness polyploidization attributes for improving crop climate resilience. This review focuses on the impact of polyploidy events on the Poaceae family, which includes the world’s most important food sources, and discusses the past, present, and future of polyploidy research for major and minor crops. The increasing accessibility to genomes of grasses and their wild progenitors together with new tools and interdisciplinary research on polyploidy can support crop improvement for global food security in the face of climate change.

Limited allele-specific gene expression in highly polyploid sugarcane

Polyploidy is widespread in plants, allowing the different copies of genes to be expressed differently in a tissue-specific or developmentally specific way. This allele-specific expression (ASE) has been widely reported, but the proportion and nature of genes showing this characteristic have not been well defined. We now report an analysis of the frequency and patterns of ASE at the whole-genome level in the highly polyploid sugarcane genome. Very high depth whole-genome sequencing and RNA sequencing revealed strong correlations between allelic proportions in the genome and in expressed sequences. This level of sequencing allowed discrimination of each of the possible allele doses in this 12-ploid genome. Most genes were expressed in direct proportion to the frequency of the allele in the genome with examples of polymorphisms being found with every possible discrete level of dose from 1:11 for single-copy alleles to 12:0 for monomorphic sites. The rarer cases of ASE were more frequent in the expression of defense-response genes, as well as in some processes related to the biosynthesis of cell walls. ASE was more common in genes with variants that resulted in significant disruption of function. The low level of ASE may reflect the recent origin of polyploid hybrid sugarcane. Much of the ASE present can be attributed to strong selection for resistance to diseases in both nature and domestication.

An autopolyploid-suitable polyBSA-seq strategy for screening candidate genetic markers linked to leaf blight resistance in sugarcane

An autopolyploid-suitable polyBSA-seq strategy was developed for screening candidate genetic markers linked to leaf blight resistance in sugarcane. Due to the complex genome architecture, the quantitative trait loci mappings and linkage marker selections for agronomic traits of autopolyploid crops were mainly limited to the time-consuming and cost intensive construction of genetic maps. To map resistance-linked markers for sugarcane leaf blight (SLB) caused by Stagonospora tainanensis, the autopolyploid-suitable bulk-segregant analysis based on the sequencing (polyBSA-seq) strategy was successfully applied for the first time. Resistant- and susceptible-bulks (R- and S-bulks) constructed from the extreme-phenotypic sugarcane F₁ lines of YT93-159 × ROC22 were deep sequenced with 195.0 × for bulks and 74.4 × for parents. Informative single-dose variants (ISDVs) present as one copy in one parent and null in the other parent were detected based on the genome sequence of LA Purple, an autooctoploid Saccharum officinarum, to screen candidate linkage markers (CLMs). The proportion of the number of short reads harboring ISDVs in the total short reads covering a given genomic position was defined as ISDV index and the ISDVs with indices met the threshold set in this study (0.04-0.14) were selected as CLMs. In total, three resistance- and one susceptibility-related CLMs for SLB resistance were identified by the polyBSA-seq. Among them, two markers on chromosome 10 were less than 300 Kb apart. Furthermore, the RNA-seq was used to calculate the expression level of genes within 1.0 Mb from the aforementioned four CLMs, which demonstrated that twelve genes were differentially expressed between resistant and susceptible clones, including a receptor-like kinase and an ethylene-responsive transcription factor. This is the first reported polyBSA-seq in autopolyploid sugarcane, which specifically tailored for the fast selection of the CLMs and causal genes associated with important agronomic traits.

Gene Editing Technologies for Sugarcane Improvement: Opportunities and Limitations

Plant-based biofuels present a promising alternative to depleting non-renewable fuel resources. One of the benefits of biofuel is reduced environmental impact, including reduction in greenhouse gas emission which causes climate change. Sugarcane is one of the most important bioenergy crops. Sugarcane juice is used to produce table sugar and first-generation biofuel (e.g., bioethanol). Sugarcane bagasse is also a potential material for second-generation cellulosic biofuel production. Researchers worldwide are striving to improve sugarcane biomass yield and quality by a variety of means including biotechnological tools. This paper reviews the use of sugarcane as a feedstock for biofuel production, and gene manipulation tools and approaches, including RNAi and genome-editing tools, such as TALENs and CRISPR-Cas9, for improving its quality. The specific focus here is on CRISPR system because it is low cost, simple in design and versatile compared to other genome-editing tools. The advance of CRISPR-Cas9 technology has transformed plant research with its ability to precisely delete, insert or replace genes in recent years. Lignin is the primary material responsible for biomass recalcitrance in biofuel production. The use of genome editing technology to modify lignin composition and distribution in sugarcane cell wall has been realized. The current and potential applications of genome editing technology for sugarcane improvement are discussed. The advantages and limitations of utilizing RNAi and TALEN techniques in sugarcane improvement are discussed as well.

Genetic manipulation of trehalose-6-phosphate synthase results in changes in the soluble sugar profile in transgenic sugarcane stems

Trehalose is a non-reducing disaccharide widely distributed in nature. The trehalose biosynthetic intermediate, trehalose 6-phosphate (Tre6P) is an essential regulatory and signaling molecule involved in both regulation of carbon metabolism and photosynthesis. To investigate the effect of altered trehalose synthesis on sucrose accumulation in sugarcane (Saccharum spp. hybrid), we independently overexpressed the Escherichia coli otsA (trehalose-6-phosphate synthase; TPS) and otsB (trehalose-6-phosphate phosphatase; TPP) genes and additionally partially silenced native TPS expression. In mature cane, sucrose levels in the otsA transgenic plants were lowered, whereas sucrose levels in the otsB transgenic plants were increased. Partial silencing of TPS expression in sugarcane transformed with a TPS-targeted microRNA recombinant construct was confirmed in leaf and mature internode tissue of transgenic plants. Most of the silencing transgenic lines accumulated trehalose at lower levels than the wild-type (WT) plants. The immature stalk tissue of these transgenic lines had lower levels of glucose and fructose, whereas the mature internode tissue had lower sucrose and glucose levels, when compared with the WT. Furthermore, various minor metabolites and sugars were detected in the sugarcane plants, which mostly decreased as the stalk tissue of the cane matured. The results demonstrate that manipulation of Tre6P/trehalose metabolism has the potential to modify the profile of soluble sugars accumulated in sugarcane stems.

Sugarcane monoculture drives microbial community composition, activity and abundance of agricultural-related microorganisms

Sugarcane monoculture (SM) often leads to soil problems, like soil acidification, degradation, and soil-borne diseases, which ultimately pose a negative impact on agricultural productivity and sustainability. Understanding the change in microbial communities' composition, activities, and functional microbial taxa associated with the plant and soil under SM is unclear. Using multidisciplinary approaches such as Illumina sequencing, measurements of soil properties, and enzyme activities, we analyzed soil samples from three sugarcane fields with different monoculture histories (1-, 2-, and 4-year cultivation times, respectively). We observed that SM induced soil acidity and had adverse effects on soil fertility, i.e., soil organic matter (OM), total nitrogen (TN), total carbon (TC), and available potassium (AK), as well as enzyme activities indicative for carbon, phosphorus, and nitrogen cycles. Non-metric multidimensional scaling (NMDS) analysis showed that SM time greatly affected soil attribute patterns. We observed strong correlation among soil enzymes activities and soil physiochemical properties (soil pH, OM, and TC). Alpha diversity analysis showed a varying response of the microbes to SM time. Bacterial diversity increased with increasing oligotrophs (e.g., Acidobacteria and Chloroflexi), while fungal diversity decreased with reducing copiotrophs (e.g., Ascomycota). β-Diversity analysis showed that SM time had a great influence on soil microbial structure and soil properties, which led to the changes in major components of microbial structure (soil pH, OM, TC, bacteria and soil pH; TC, fungi). Additionally, SM time significantly stimulated (four bacterial and ten fungal) and depleted (12 bacterial and three fungal) agriculturally and ecologically important microbial genera that were strongly and considerably correlated with soil characteristics (soil pH, OM, TC, and AK). In conclusion, SM induces soil acidity, reduces soil fertility, shifts microbial structure, and reduces its activity. Furthermore, most beneficial bacterial genera decreased significantly due to SM, while beneficial fungal genera showed a reverse trend. Therefore, mitigating soil acidity, improving soil fertility, and soil enzymatic activities, including improved microbial structure with beneficial service to plants and soil, can be an effective measure to develop a sustainable sugarcane cropping system.

Unraveling the Genome of a High Yielding Colombian Sugarcane Hybrid

Recent developments in High Throughput Sequencing (HTS) technologies and bioinformatics, including improved read lengths and genome assemblers allow the reconstruction of complex genomes with unprecedented quality and contiguity. Sugarcane has one of the most complicated genomes among grassess with a haploid length of 1Gbp and a ploidies between 8 and 12. In this work, we present a genome assembly of the Colombian sugarcane hybrid CC 01-1940. Three types of sequencing technologies were combined for this assembly: PacBio long reads, Illumina paired short reads, and Hi-C reads. We achieved a median contig length of 34.94 Mbp and a total genome assembly of 903.2 Mbp. We annotated a total of 63,724 protein coding genes and performed a reconstruction and comparative analysis of the sucrose metabolism pathway. Nucleotide evolution measurements between orthologs with close species suggest that divergence between Saccharum officinarum and Saccharum spontaneum occurred <2 million years ago. Synteny analysis between CC 01-1940 and the S. spontaneum genome confirms the presence of translocation events between the species and a random contribution throughout the entire genome in current sugarcane hybrids. Analysis of RNA-Seq data from leaf and root tissue of contrasting sugarcane genotypes subjected to water stress treatments revealed 17,490 differentially expressed genes, from which 3,633 correspond to genes expressed exclusively in tolerant genotypes. We expect the resources presented here to serve as a source of information to improve the selection processes of new varieties of the breeding programs of sugarcane.

Near-infrared spectroscopy outperforms genomics for predicting sugarcane feedstock quality traits

The main objectives of this study were to evaluate the prediction performance of genomic and near-infrared spectroscopy (NIR) data and whether the integration of genomic and NIR predictor variables can increase the prediction accuracy of two feedstock quality traits (fiber and sucrose content) in a sugarcane population (Saccharum spp.). The following three modeling strategies were compared: M1 (genome-based prediction), M2 (NIR-based prediction), and M3 (integration of genomics and NIR wavenumbers). Data were collected from a commercial population comprised of three hundred and eighty-five individuals, genotyped for single nucleotide polymorphisms and screened using NIR spectroscopy. We compared partial least squares (PLS) and BayesB regression methods to estimate marker and wavenumber effects. In order to assess model performance, we employed random sub-sampling cross-validation to calculate the mean Pearson correlation coefficient between observed and predicted values. Our results showed that models fitted using BayesB were more predictive than PLS models. We found that NIR (M2) provided the highest prediction accuracy, whereas genomics (M1) presented the lowest predictive ability, regardless of the measured traits and regression methods used. The integration of predictors derived from NIR spectroscopy and genomics into a single model (M3) did not significantly improve the prediction accuracy for the two traits evaluated. These findings suggest that NIR-based prediction can be an effective strategy for predicting the genetic merit of sugarcane clones.

Bracing for sustainable agriculture: the development and function of brace roots in members of Poaceae

Optimization of crop production requires root systems to function in water uptake, nutrient use, and anchorage. In maize, two types of nodal roots-subterranean crown and aerial brace roots function in anchorage and water uptake and preferentially express multiple water and nutrient transporters. Brace root development shares genetic control with juvenile-to-adult phase change and flowering time. We present a comprehensive list of the genes known to alter brace roots and explore these as candidates for QTL studies in maize and sorghum. Brace root development and function may be conserved in other members of Poaceae, however research is limited. This work highlights the critical knowledge gap of aerial nodal root development and function and suggests new focus areas for breeding resilient crops.

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.

Shallower Root Spatial Distribution Induced by Phosphorus Deficiency Contributes to Topsoil Foraging and Low Phosphorus Adaption in Sugarcane (Saccharum officinarum L.)

Low phosphorus (P) availability in acid soils is one of the main limiting factors in sugarcane (Saccharum officinarum L.) production. Reconstruction of the root system architecture (RSA) is a vital mechanism for crop low P adaption, while the RSA of sugarcane has not been studied in detail because of its complex root system. In this study, reconstruction of the RSA and its relationship with P acquisition were investigated in a P-efficient sugarcane genotype ROC22 (R22) and two P-inefficient genotypes Yunzhe 03-103 (YZ) and Japan 2 (JP). An efficient dynamic observation room was developed to monitor the spatiotemporal alternation of sugarcane root length density (RLD) and root distribution in soil with heterogeneous P locations. The sugarcane RSA was reconstructed under P deficiency, and R22 had an earlier response than YZ and JP and presented an obvious feature of root shallowness. Compared with the normal P condition, the shallow RLD was increased by 112% in R22 under P deficiency while decreased by 26% in YZ and not modified in JP. Meanwhile, R22 exhibited a shallower root distribution than YZ and JP under P deficiency, supported by 51 and 24% greater shallow RLD, and 96 and 67% greater shallow root weight, respectively. The ratio of shallow RLD to total RLD in R22 was 91% greater than YZ, and the ratio of shallow root weight to total root weight in R22 was greater than that of YZ and JP by 94 and 30%, respectively. As a result, R22 had a higher shoot P accumulation than YZ and JP, which thereby increased the relative leaf sheath inorganic P concentration (RLPC) by 47 and 56%, relative shoot biomass (RSB) by 36 and 33%, and relative cane weight (RCW) by 31 and 36%, compared with YZ and JP under P deficiency, respectively. We verified the reliability and efficiency of a dynamic observation room and demonstrated that a shallower root distribution contributed to improving topsoil foraging, P acquisition, and low P adaption under P deficiency in sugarcane. Therefore, a shallower root distribution merits consideration as an evaluation trait for breeding P efficient sugarcane genotypes and genetic improvement.

Leaping into the Unknown World of Sporisorium scitamineum Candidate Effectors

Sporisorium scitamineum is a biotrophic fungus causing sugarcane smut disease. In this study, we set up a pipeline and used genomic and dual transcriptomic data previously obtained by our group to identify candidate effectors of S. scitamineum and their expression profiles in infected smut-resistant and susceptible sugarcane plants. The expression profile of different genes after infection in contrasting sugarcane genotypes assessed by RT-qPCR depended on the plant genotypes and disease progression. Three candidate effector genes expressed earlier only in resistant plants, four expressed in both genotypes, and three later in susceptible plants. Ten genes were cloned and transiently expressed in N. benthamiana leaves to determine their subcellular location, while four localized in more than one compartment. Two candidates, g3890 having a nucleoplasmic and mitochondrial location and g5159 targeting the plant cell wall, were selected to obtain their possible corresponding host targets using co-immunoprecipitation (CoIP) experiments and mass spectrometry. Various potential interactors were identified, including subunits of the protein phosphatase 2A and an endochitinase. We investigated the presence of orthologs in sugarcane and using transcriptome data present their expression profiles. Orthologs of sugarcane shared around 70% similarity. Identifying a set of putative fungal effectors and their plant targets provides a valuable resource for functional characterization of the molecular events leading to smut resistance in sugarcane plants and uncovers further opportunities for investigation.

Transcriptional reprogramming of major defense-signaling pathways during defense priming and sugarcane-Colletotrichum falcatum interaction

Red rot caused by Colletotrichum falcatum poses a serious threat to sugarcane cultivation in many tropical and sub-tropical countries. Deciphering the molecular network of major defense-signaling pathways in sugarcane cultivars with varying red rot resistance is essential to elucidate the phenomenon of defense priming exerted by resistance inducers. Therefore, in this study, expression pattern of transcripts coding for major defense-signaling pathway regulatory genes was profiled during compatible and incompatible interactions and in response to defense priming using qRT-PCR. Candidate genes that were profiled are involved in or related to hypersensitive response and reactive oxygen species production (HR/ROS), salicylic acid (SA), and jasmonic acid/ethylene (JA/ET) pathways. For compatible and incompatible interactions, susceptible (CoC 671), field tolerant (Co 86032) and resistant (Co 93009) sugarcane cultivars were used, whereas for defense priming, benzothiadiazole (BTH) and the pathogen-associated molecular patterns (PAMPs) of C. falcatum viz., CfEPL1 (eliciting plant response-like) and CfPDIP1 (plant defense inducing protein) were used in CoC 671 cultivar. Results indicated that the master regulator of defense pathways, nonexpressor of pathogenesis-related genes 1 (NPR1) was highly upregulated in incompatible interactions (in both Co 86032 and Co 93009) than the compatible interaction along with SA pathway-associated genes. Similarly, in response to defense priming with BTH, CfEPL1 and CfPDIP1, only the SA pathway-associated genes showed considerable upregulation at 0 h post inoculation (hpi) and other intermittent time points. Overall, this study showed that SA-mediated defense pathway is the most predominant pathway reprogrammed during priming with BTH, CfEPL1 and CfPDIP1 and substantiated the earlier findings that these agents indeed induce systemic resistance against red rot of sugarcane.

Differential expression in leaves of Saccharum genotypes contrasting in biomass production provides evidence of genes involved in carbon partitioning

BACKGROUND

The development of biomass crops aims to meet industrial yield demands, in order to optimize profitability and sustainability. Achieving these goals in an energy crop like sugarcane relies on breeding for sucrose accumulation, fiber content and stalk number. To expand the understanding of the biological pathways related to these traits, we evaluated gene expression of two groups of genotypes contrasting in biomass composition.

RESULTS

First visible dewlap leaves were collected from 12 genotypes, six per group, to perform RNA-Seq. We found a high number of differentially expressed genes, showing how hybridization in a complex polyploid system caused extensive modifications in genome functioning. We found evidence that differences in transposition and defense related genes may arise due to the complex nature of the polyploid Saccharum genomes. Genotypes within both biomass groups showed substantial variability in genes involved in photosynthesis. However, most genes coding for photosystem components or those coding for phosphoenolpyruvate carboxylases (PEPCs) were upregulated in the high biomass group. Sucrose synthase (SuSy) coding genes were upregulated in the low biomass group, showing that this enzyme class can be involved with sucrose synthesis in leaves, similarly to sucrose phosphate synthase (SPS) and sucrose phosphate phosphatase (SPP). Genes in pathways related to biosynthesis of cell wall components and expansins coding genes showed low average expression levels and were mostly upregulated in the high biomass group.

CONCLUSIONS

Together, these results show differences in carbohydrate synthesis and carbon partitioning in the source tissue of distinct phenotypic groups. Our data from sugarcane leaves revealed how hybridization in a complex polyploid system resulted in noticeably different transcriptomic profiles between contrasting genotypes.

Innovations in plant genetics adapting agriculture to climate change

Developing new genotypes of plants is one of the key options for adaptation of agriculture to climate change. Plants may be required to provide resilience in changed climates or support the migration of agriculture to new regions. Very different genotypes may be required to perform in the modified environments of protected agriculture. Consumers will continue to demand taste, convenience, healthy and safe food and sustainably and ethically produced food, despite the greater challenges of climate in the future. Improving the nutritional value of foods in response to climate change is a significant challenge. Genomic sequences of relevant germplasm and an understanding of the functional role of alleles controlling key traits will be an enabling platform for this innovation.

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.

Molecular diversity and genetic structure of Saccharum complex accessions

Sugarcane is an important crop for food and energy security, providing sucrose and bioethanol from sugar content and bioelectricity from lignocellulosic bagasse. In order to evaluate the diversity and genetic structure of the Brazilian Panel of Sugarcane Genotypes (BPSG), a core collection composed by 254 accessions of the Saccharum complex, eight TRAP markers anchored in sucrose and lignin metabolism genes were evaluated. A total of 584 polymorphic fragments were identified and used to investigate the genetic structure of BPSG through analysis of molecular variance (AMOVA), principal components analysis (PCA), a Bayesian method using STRUCTURE software, genetic dissimilarity and phylogenetic tree. AMOVA showed a moderate genetic differentiation between ancestors and improved accessions, 0.14, and the molecular variance was higher within populations than among populations, with values of 86%, 95% and 97% when constrasting improved with ancestors, foreign with ancestors and improved with foreign, respectively. The PCA approach suggests clustering in according with evolutionary and Brazilian breeding sugarcane history, since improved accessions from older generations were positioned closer to ancestors than improved accessions from recent generations. This result was also confirmed by STRUCTURE analysis and phylogenetic tree. The Bayesian method was able to separate ancestors of the improved accessions while the phylogenetic tree showed clusters considering the family relatedness within three major clades; the first being composed mainly by ancestors and the other two mainly by improved accessions. This work can contribute to better management of the crosses considering functional regions of the sugarcane genome.

A flow cytometry-based analysis to establish a cell cycle synchronization protocol for Saccharum spp

Modern sugarcane is an unusually complex heteroploid crop, and its genome comprises two or three subgenomes. To reduce the complexity of sugarcane genome research, the ploidy level and number of chromosomes can be reduced using flow chromosome sorting. However, a cell cycle synchronization (CCS) protocol for Saccharum spp. is needed that maximizes the accumulation of metaphase chromosomes. For flow cytometry analysis in this study, we optimized the lysis buffer, hydroxyurea(HU) concentration, HU treatment time and recovery time for sugarcane. We determined the mitotic index by microscopic observation and calculation. We found that WPB buffer was superior to other buffers for preparation of sugarcane nuclei suspensions. The optimal HU treatment was 2 mM for 18 h at 25 °C, 28 °C and 30 °C. Higher recovery treatment temperatures were associated with shorter recovery times (3.5 h, 2.5 h and 1.5 h at 25 °C, 28 °C and 30 °C, respectively). The optimal conditions for treatment with the inhibitor of microtubule polymerization, amiprophos-methyl (APM), were 2.5 μM for 3 h at 25 °C, 28 °C and 30 °C. Meanwhile, preliminary screening of CCS protocols for Badila were used for some main species of genus Saccharum at 25 °C, 28 °C and 30 °C, which showed that the average mitotic index decreased from 25 °C to 30 °C. The optimal sugarcane CCS protocol that yielded a mitotic index of >50% in sugarcane root tips was: 2 mM HU for 18 h, 0.1 X Hoagland's Solution without HU for 3.5 h, and 2.5 μM APM for 3.0 h at 25 °C. The CCS protocol defined in this study should accelerate the development of genomic research and cytobiology research in sugarcane.

Flow cytometric characterisation of the complex polyploid genome of Saccharum officinarum and modern sugarcane cultivars

Sugarcane (Saccharum spp.) is a globally important crop for sugar and bioenergy production. Its highly polyploid, complex genome has hindered progress in understanding its molecular structure. Flow cytometric sorting and analysis has been used in other important crops with large genomes to dissect the genome into component chromosomes. Here we present for the first time a method to prepare suspensions of intact sugarcane chromosomes for flow cytometric analysis and sorting. Flow karyotypes were generated for two S. officinarum and three hybrid cultivars. Five main peaks were identified and each genotype had a distinct flow karyotype profile. The flow karyotypes of S. officinarum were sharper and with more discrete peaks than the hybrids, this difference is probably due to the double genome structure of the hybrids. Simple Sequence Repeat (SSR) markers were used to determine that at least one allelic copy of each of the 10 basic chromosomes could be found in each peak for every genotype, except R570, suggesting that the peaks may represent ancestral Saccharum sub genomes. The ability to flow sort Saccharum chromosomes will allow us to isolate and analyse chromosomes of interest and further examine the structure and evolution of the sugarcane genome.

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

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

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.

What Has Been the Focus of Sugarcane Research? A Bibliometric Overview

Sugarcane is one of the main crops worldwide, and it has an important impact on environmental issues. A bibliometric mapping analysis of the research on sugarcane was carried out, using data on the titles, abstracts, and keywords of articles published in leading journals and other peer-reviewed documents available in the SCOPUS database from 1858 to 2019 (27 August), and this was subsequently analyzed with the software VOSviewer. The three most important countries that publish research and were most-cited regarding sugarcane were Brazil, the USA, and India. The analysis of the co-occurrence of terms shows that the main research areas were sugarcane bagasse and terms related to bioenergy, and on a second level of relevance agronomy topics related to increasing crop yields. This first attempt to visualize the abundance of publications regarding sugarcane in their totality is in itself a good starting point for further scientific discussion.

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

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

A genome-wide association study identified loci for yield component traits in sugarcane (Saccharum spp.)

Sugarcane (Saccharum spp.) has a complex genome with variable ploidy and frequent aneuploidy, which hampers the understanding of phenotype and genotype relations. Despite this complexity, genome-wide association studies (GWAS) may be used to identify favorable alleles for target traits in core collections and then assist breeders in better managing crosses and selecting superior genotypes in breeding populations. Therefore, in the present study, we used a diversity panel of sugarcane, called the Brazilian Panel of Sugarcane Genotypes (BPSG), with the following objectives: (i) estimate, through a mixed model, the adjusted means and genetic parameters of the five yield traits evaluated over two harvest years; (ii) detect population structure, linkage disequilibrium (LD) and genetic diversity using simple sequence repeat (SSR) markers; (iii) perform GWAS analysis to identify marker-trait associations (MTAs); and iv) annotate the sequences giving rise to SSR markers that had fragments associated with target traits to search for putative candidate genes. The phenotypic data analysis showed that the broad-sense heritability values were above 0.48 and 0.49 for the first and second harvests, respectively. The set of 100 SSR markers produced 1,483 fragments, of which 99.5% were polymorphic. These SSR fragments were useful to estimate the most likely number of subpopulations, found to be four, and the LD in BPSG, which was stronger in the first 15 cM and present to a large extension (65 cM). Genetic diversity analysis showed that, in general, the clustering of accessions within the subpopulations was in accordance with the pedigree information. GWAS performed through a multilocus mixed model revealed 23 MTAs, six, three, seven, four and three for soluble solid content, stalk height, stalk number, stalk weight and cane yield traits, respectively. These MTAs may be validated in other populations to support sugarcane breeding programs with introgression of favorable alleles and marker-assisted selection.

Analysis of the diversity and tissue specificity of sucrose synthase genes in the long read transcriptome of sugarcane

BACKGROUND

Sugarcane accumulates very high levels of sucrose in the culm. Elucidation of the molecular mechanisms that allows such high sucrose synthesis and accumulation (up to 650 mM) is made difficult by the complexity of the highly polyploid genome. Here we report the use of RNA Seq data to characterize the sucrose synthase (SuSy) genes expressed in the transcriptome of the mature sugarcane plant.

RESULTS

Four SuSy gene families were identified in the sugarcane Iso-Seq long read transcriptome (SUGIT) through gene annotation of transcripts that mapped to reference SuSy genes from sorghum and maize. In total, 38, 19, 14, and 2 transcripts were identified for the four corresponding SuSy genes 1, 2, 4 and 7, respectively. Comparative studies using available SuSy genes from sorghum (1, 2, 4, 6, 7) and maize (1-7) revealed that the sugarcane SuSy genes were interrupted by multiple introns and that they share a highly conserved gene structure. Spatial expression of the four SuSy genes in sugarcane genotypes and in the progenitor species, Saccharum spontaneum and Saccharum officinarum, was studied in the leaf and root tissues and also in three regions of the culm tissue; top, middle and bottom internodes. Expression profiles indicated that all SuSy transcripts were differentially expressed between the top and bottom tissues, with high expression in the top tissues, lower expression in the bottom and moderate expression in the middle, indicating a gradient of SuSy activity in the sugarcane culm. Further, the root tissue had similar expression levels to that of the top internodes while leaf tissues showed lower expression. In the progenitors, SuSy7 was found to be highly expressed in S. officinarum while the other three SuSy genes had moderate expression in both the progenitors.

CONCLUSIONS

The high expression of the SuSy genes in sink tissues, the top internodes and the roots suggests functional roles in sucrose utilization to support growth. The SuSy7 gene has not been previously reported in sugarcane. As sugarcane is unique in storing such high amounts of sucrose, it is possible that there are more SuSy genes/isoforms with specific expression patterns to be discovered in this complex system.

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 and miRNAs analyses enhance our understanding of the evolutionary advantages of polyploidy

Polyploid organisms have more than two sets of chromosomes, including autopolyploid via intraspecific genome doubling, and allopolyploid via merging genomes of distinct species by hybridization. Polyploid organisms are widespread in plants, indicating that polyploidy has some evolutionary advantages over its diploid ancestor. Actually, polyploidy is always tightly associated with hybrid vigor and adaptation to adverse environmental conditions. However, why polyploidy can develop such advantages is poorly known. MicroRNAs (miRNAs) are endogenous ∼21 nt small RNAs which can play important regulatory roles in animals and plants by targeting mRNAs for cleavage or translational repression. MicroRNAs are essential for cell development, differentiation, signal transduction, and show an adaptive response to biotic and abiotic stresses. Environmental stresses cause plants to over- or under-express certain miRNAs or synthesize new miRNAs to cope with stress. We have here reviewed our current knowledge on the molecular mechanisms, which can account for the evolutionary advantages of polyploidy over its diploid ancestor from genome-wide gene expression and microRNAs expression perspectives.