Genome-Wide Association Study Unravels Quantitative Trait Loci and Genes Associated with Yield-Related Traits in Sugarcane.

Sugarcane, a major sugar and energy crop worldwide faces an increasing demand for higher yields. Identifying yield-related markers and candidate genes is valuable for breeding high-yield varieties using molecular techniques. In this work, seven yield-related traits were evaluated in a diversity panel of 159 genotypes, derived from Tripidium arundinaceum, Saccharum spontaneum, and modern sugarcane genotypes. All traits exhibited significant genetic variance with high heritability and high correlations. Genetic diversity analysis reveals a genomic decay of 23 kb and an average single nucleotide polymorphism (SNP) number of 25,429 per genotype. These 159 genotypes were divided into 4 subgroups. Genome-wide association analysis identified 47 SNPs associated with brix, spanning 36 quantitative trait loci (QTLs), and 138 SNPs for other traits across 104 QTLs, covering all 32 chromosomes. Interestingly, 12 stable QTLs associated with yield-related traits were identified, which contained 35 candidate genes. This work provides markers and candidate genes for marker-assisted breeding to improve sugarcane yields.

Cold-Induced Physiological and Biochemical Alternations and Proteomic Insight into the Response of Saccharum spontaneum to Low Temperature.

Sugarcane, a cash crop, is easily affected by low temperature, which results in a decrease in yield and sugar production. Breeding a new variety with cold tolerance is an essential strategy to reduce loss from cold stress. The identification of germplasms and genes/proteins with cold tolerance is a vital step in breeding sugarcane varieties with cold tolerance via a conventional program and molecular technology. In this study, the physiological and biochemical indices of 22 genotypes of S. spontaneum were measured, and the membership function analysis method was used to comprehensively evaluate the cold tolerance ability of these genotypes. The physiological and biochemical indices of these S. spontaneum genotypes showed a sophisticated response to low temperature. On the basis of the physiological and chemical indices, the genotypes were classified into different cold tolerance groups. Then, the high-tolerance genotype 1027 and the low-tolerance genotype 3217 were selected for DIA-based proteomic analysis by subjecting them to low temperature. From the four comparison groups, 1123, 1341, 751, and 1693 differentially abundant proteins (DAPs) were identified, respectively. The DAPs based on genotypes or treatments participated in distinct metabolic pathways. Through detailed analysis of the DAPs, some proteins related to protein homeostasis, carbohydrate and energy metabolism, amino acid transport and metabolism, signal transduction, and the cytoskeleton may be involved in sugarcane tolerance to cold stress. Furthermore, five important proteins related to cold tolerance were discovered for the first time in this study. This work not only provides the germplasms and candidate target proteins for breeding sugarcane varieties with cold tolerance via a conventional program and molecular breeding, but also helps to accelerate the determination of the molecular mechanism underlying cold tolerance in sugarcane.

MicroSugar: A database of comprehensive miRNA target prediction framework for sugarcane (Saccharum officinarum L.).

microRNA (miRNA) is a group of small non-coding RNA that plays important role in post-transcription of gene expression. With the studies about miRNA increase in sugarcane, the researchers lack an exhaustive resource to achieve the data. To fill this gap, we developed MicroSugar, a database that supported mRNA and miRNA annotation for sugarcane (http://suc.gene-db.com). MicroSugar is an integrated resource developed for 194,528 genes including 80,746 unigenes from long reads of Pacbio platform and 468 miRNAs from 72 samples. Internode elongation (jointing) is the key biological characteristic for the growth of sugarcane tillers into sugarcane stems. The present study combined the sequencing data from the different stages in internode elongation of stem and tiller. In total, the 14,300 3' untranslated region (UTR) sequences were extracted from the gene sequences and 3019 mRNAs as target of 327 miRNA were identified by miRanda algorithm and Spearman's Rho of expression levels. To determine the gene functions regulated by these miRNAs, the gene ontology enrichment analysis was performed and it confirmed that the over-represented Gene Ontology (GO) terms were associated with organism formation indicating the growth controlling function by miRNAs in sugarcane. Moreover, MicroSugar is a comprehensive and integrated database with a user-friendly responsive template. By browsing, searching and downloading of the nucleotide sequences, expression and miRNA targets, the user can retrieve information promptly. The database provides a valuable resource to facilitate the understanding of miRNA in sugarcane development and growth which will contribute to the study of sugarcane and other plants.

Transcriptome and metabolome analyses reveal new insights into chlorophyll, photosynthesis, metal ion and phenylpropanoids related pathways during sugarcane ratoon chlorosis.

BACKGROUND: Ratoon sugarcane is susceptible to chlorosis, characterized by chlorophyll loss, poor growth, and a multitude of nutritional deficiency mainly occurring at young stage. Chlorosis would significantly reduce the cane production. The molecular mechanism underlying this phenomenon remains unknown. We analyzed the transcriptome and metabolome of chlorotic and non-chlorotic sugarcane leaves of the same age from the same field to gain molecular insights into this phenomenon. RESULTS: The agronomic traits, such as plant height and the number of leaf, stalk node, and tillers declined in chlorotic sugarcane. Chlorotic leaves had substantially lower chlorophyll content than green leaves. A total of 11,776 differentially expressed genes (DEGs) were discovered in transcriptome analysis. In the KEGG enriched chlorophyll metabolism pathway, sixteen DEGs were found, eleven of which were down-regulated. Two photosynthesis pathways were also enriched with 32 genes downregulated and four genes up-regulated. Among the 81 enriched GO biological processes, there were four categories related to metal ion homeostasis and three related to metal ion transport. Approximately 400 metabolites were identified in metabolome analysis. The thirteen differentially expressed metabolites (DEMs) were all found down-regulated. The phenylpropanoid biosynthesis pathway was enriched in DEGs and DEMs, indicating a potentially vital role for phenylpropanoids in chlorosis. CONCLUSIONS: Chlorophyll production, metal ion metabolism, photosynthesis, and some metabolites in the phenylpropanoid biosynthesis pathway were considerably altered in chlorotic ratoon sugarcane leaves. Our finding revealed the relation between chlorosis and these pathways, which will help expand our mechanistic understanding of ratoon sugarcane chlorosis.

Characterization of full-length transcriptome in Saccharum officinarum and molecular insights into tiller development.

BACKGROUND: Although extensive breeding efforts are ongoing in sugarcane (Saccharum officinarum L.), the average yield is far below the theoretical potential. Tillering is an important component of sugarcane yield, however, the molecular mechanism underlying tiller development is still elusive. The limited genomic data in sugarcane, particularly due to its complex and large genome, has hindered in-depth molecular studies. RESULTS: Herein, we generated full-length (FL) transcriptome from developing leaf and tiller bud samples based on PacBio Iso-Seq. In addition, we performed RNA-seq from tiller bud samples at three developmental stages (T0, T1 and T2) to uncover key genes and biological pathways involved in sugarcane tiller development. In total, 30,360 and 20,088 high-quality non-redundant isoforms were identified in leaf and tiller bud samples, respectively, representing 41,109 unique isoforms in sugarcane. Likewise, we identified 1063 and 1037 alternative splicing events identified in leaf and tiller bud samples, respectively. We predicted the presence of coding sequence for 40,343 isoforms, 98% of which was successfully annotated. Comparison with previous FL transcriptomes in sugarcane revealed 2963 unreported isoforms. In addition, we characterized 14,946 SSRs from 11,700 transcripts and 310 lncRNAs. By integrating RNA-seq with the FL transcriptome, 468 and 57 differentially expressed genes (DEG) were identified in T1vsT0 and T2vsT0, respectively. Strong up-regulation of several pyruvate phosphate dikinase and phosphoenolpyruvate carboxylase genes suggests enhanced carbon fixation and protein synthesis to facilitate tiller growth. Similarly, up-regulation of linoleate 9S-lipoxygenase and lipoxygenase genes in the linoleic acid metabolism pathway suggests high synthesis of key oxylipins involved in tiller growth and development. CONCLUSIONS: Collectively, we have enriched the genomic data available in sugarcane and provided candidate genes for manipulating tiller formation and development, towards productivity enhancement in sugarcane.

Global transcriptome changes of elongating internode of sugarcane in response to mepiquat chloride.

BACKGROUND: Mepiquat chloride (DPC) is a chemical that is extensively used to control internode growth and create compact canopies in cultured plants. Previous studies have suggested that DPC could also inhibit gibberellin biosynthesis in sugarcane. Unfortunately, the molecular mechanism underlying the suppressive effects of DPC on plant growth is still largely unknown. RESULTS: In the present study, we first obtained high-quality long transcripts from the internodes of sugarcane using the PacBio Sequel System. A total of 72,671 isoforms, with N50 at 3073, were generated. These long isoforms were used as a reference for the subsequent RNA-seq. Afterwards, short reads generated from the Illumina HiSeq 4000 platform were used to compare the differentially expressed genes in both the DPC and the control groups. Transcriptome profiling showed that most significant gene changes occurred after six days post DPC treatment. These genes were related to plant hormone signal transduction and biosynthesis of several metabolites, indicating that DPC affected multiple pathways, in addition to suppressing gibberellin biosynthesis. The network of DPC on the key stage was illustrated by weighted gene co-expression network analysis (WGCNA). Among the 36 constructed modules, the top positive correlated module, at the stage of six days post spraying DPC, was sienna3. Notably, Stf0 sulfotransferase, cyclin-like F-box, and HOX12 were the hub genes in sienna3 that had high correlation with other genes in this module. Furthermore, the qPCR validated the high accuracy of the RNA-seq results. CONCLUSION: Taken together, we have demonstrated the key role of these genes in DPC-induced growth inhibition in sugarcane.

Enhanced Activity of Genes Associated With Photosynthesis, Phytohormone Metabolism and Cell Wall Synthesis Is Involved in Gibberellin-Mediated Sugarcane Internode Growth.

Internode elongation is an important trait in sugarcane as it affects the sugarcane yield. Gibberellin (GA) is a key modulator of internode elongation in sugarcane. Understanding the gene expression features of GA-mediated internode elongation has both scientific and practical significance. This study aimed to examine the transcriptomic changes in the internode elongation of sugarcane following GA treatment. Eighteen cDNA libraries from the internode tissues on days of 0, 3, and 6 in control and GA treatment groups were sequenced and their gene expression were studied. RNA-seq analysis revealed 1,338,723,248 reads and 70,821 unigenes from elongating internodes of sugarcane. Comparative studies discovered a large number of transcripts that were differentially expressed in GA-treated samples compared to the control. Further analysis revealed that the differentially expressed genes were enriched in the metabolic process, one-carbon compound transport, and single-organism process. Kyoto Encyclopedia of Genes and Genomes pathway annotation showed significant enrichment in photosynthesis and plant hormone signal transduction, indicating its involvement in internode elongation. The function analysis suggested that metabolic pathways and biosynthesis of secondary metabolites, plant hormones, and cell wall components were enriched in the internodes of the GA-treated plants. The hub genes were identified, with the function of cellulose synthesis. The results of this study provide a global view of mRNA changes during sugarcane internode elongation and extend our knowledge of the GA-mediated cellular processes involved in sugarcane stem growth.

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

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

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

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