Characterisation of the double genome structure of modern sugarcane cultivars (Saccharum spp.) by molecular cytogenetics

Whole-genome sequencing of a worldwide collection of sugarcane cultivars (Saccharum spp.) reveals the genetic basis of cultivar improvement

Sugarcane is the main source of sugar worldwide, and 80% of the sucrose production comes from sugarcane. However, the genetic differentiation and basis of agronomic traits remain obscure. Here, we sequenced the whole-genome of 219 elite worldwide sugarcane cultivar accessions. A total of approximately 6 million high-quality genome-wide single nucleotide polymorphisms (SNPs) were detected. A genome-wide association study identified a total of 2198 SNPs that were significantly associated with sucrose content, stalk number, plant height, stalk diameter, cane yield, and sugar yield. We observed homozygous tendency of favor alleles of these loci, and over 80% of cultivar accessions carried the favor alleles of the SNPs or haplotypes associated with sucrose content. Gene introgression analysis showed that the number of chromosome segments from Saccharum spontaneum decreased with the breeding time of cultivars, while those from S. officinarum increased in recent cultivars. A series of selection signatures were identified in sugarcane improvement procession, of which 104 were simultaneously associated with agronomic traits and 45 of them were mainly associated with sucrose content. We further proposed that as per sugarcane transgenic experiments, ShN/AINV3.1 plays a positive role in increasing stalk number, plant height, and stalk diameter. These findings provide comprehensive resources for understanding the genetic basis of agronomic traits and will be beneficial to germplasm innovation, screening molecular markers, and future sugarcane cultivar improvement.

First investigation into the genetic control of meiosis in sugarcane

The sugarcane (Saccharum spp.) genome is one of the most complex of all. Modern varieties are highly polyploid and aneuploid as a result of hybridization between Saccharum officinarum and S. spontaneum. Little research has been done on meiotic control in polyploid species, with the exception of the wheat Ph1 locus harboring the ZIP4 gene (TaZIP4-B2) which promotes pairing between homologous chromosomes while suppressing crossover between homeologs. In sugarcane, despite its interspecific origin, bivalent association is favored, and multivalents, if any, are resolved at the end of prophase I. Thus, our aim herein was to investigate the purported genetic control of meiosis in the parental species and in sugarcane itself. We investigated the ZIP4 gene and immunolocalized meiotic proteins, namely synaptonemal complex proteins Zyp1 and Asy1. The sugarcane ZIP4 gene is located on chromosome 2 and expressed more abundantly in flowers, a similar profile to that found for TaZIP4-B2. ZIP4 expression is higher in S. spontaneum a neoautopolyploid, with lower expression in S. officinarum, a stable octoploid species. The sugarcane Zip4 protein contains a TPR domain, essential for scaffolding. Its 3D structure was also predicted, and it was found to be very similar to that of TaZIP4-B2, reflecting their functional relatedness. Immunolocalization of the Asy1 and Zyp1 proteins revealed that S. officinarum completes synapsis. However, in S. spontaneum and SP80-3280 (a modern variety), no nuclei with complete synapsis were observed. Importantly, our results have implications for sugarcane cytogenetics, genetic mapping, and genomics.

Photoperiodic and lighting treatments for flowering control and its genetic regulation in sugarcane breeding

Improvement of sugarcane is hampered due to its narrow genetic base, and the difficulty in synchronizing flowering further hinders the exploitation of the genetic potential of available germplasm resources. Therefore, the continuous evaluation and optimization of flowering control and induction techniques are vital for sugarcane improvement. In view of this, the review was conducted to investigate the current understanding of photoperiodic and lighting treatment effects on sugarcane flowering and its genetic regulation. Photoperiod facilities have made a significant contribution to flowering control in sugarcane; however, inductive photoperiods are still unknown for some genotypes, and some intended crosses are still impossible to produce because of unresponsive varieties. The effectiveness of lower red/far-red ratios in promoting sugarcane flowering has been widely understood. Furthermore, there is vast potential for utilizing blue, red, and far-red light wavelengths in the flowering control of sugarcane. In this context, light-emitting diodes (LEDs) remain efficient sources of light. Therefore, the combined use of photoperiod regimes with different light wavelengths and optimization of such treatment combinations might help to control and induce flowering in sugarcane parental clones. In sugarcane, FLOWERING LOCUS T (ScFT) orthologues from ScFT1 to ScFT13 have been identified, and interestingly, ScFT3 has evidently been identified as a floral inducer in sugarcane. However, independent assessments of different FT-like gene family members are recommended to comprehensively understand their role in the regulation of flowering. Similarly, we believe this review provides substantial information that is vital for the manipulation of flowering and exploitation of germplasm resources in sugarcane breeding.

Transcriptome and small RNA analysis unveils novel insights into the C4 gene regulation in sugarcane

MAIN CONCLUSION

This study reveals miRNA indirect regulation of C4 genes in sugarcane through transcription factors, highlighting potential key regulators like SsHAM3a. C4 photosynthesis is crucial for the high productivity and biomass of sugarcane, however, the miRNA regulation of C4 genes in sugarcane remains elusive. We have identified 384 miRNAs along the leaf gradients, including 293 known miRNAs and 91 novel miRNAs. Among these, 86 unique miRNAs exhibited differential expression patterns, and we identified 3511 potential expressed targets of these differentially expressed miRNAs (DEmiRNAs). Analyses using Pearson correlation coefficient (PCC) and Gene Ontology (GO) enrichment revealed that targets of miRNAs with positive correlations are integral to chlorophyll-related photosynthetic processes. In contrast, negatively correlated pairs are primarily associated with metabolic functions. It is worth noting that no C4 genes were predicted as targets of DEmiRNAs. Our application of weighted gene co-expression network analysis (WGCNA) led to a gene regulatory network (GRN) suggesting miRNAs might indirectly regulate C4 genes via transcription factors (TFs). The GRAS TF SsHAM3a emerged as a potential regulator of C4 genes, targeted by miR171y and miR171am, and exhibiting a negative correlation with miRNA expression along the leaf gradient. This study sheds light on the complex involvement of miRNAs in regulating C4 genes, offering a foundation for future research into enhancing sugarcane's photosynthetic efficiency.

Unveiling the predominance of Saccharum spontaneum alleles for resistance to orange rust in sugarcane using genome-wide association

KEY MESSAGE

Six QTLs of resistance to sugarcane orange rust were identified in modern interspecific hybrids by GWAS. For five of them, the resistance alleles originated from S. spontaneum. Altogether, they efficiently predict disease resistance. Sugarcane orange rust (SOR) is a threatening emerging disease in many sugarcane industries worldwide. Improving the genetic resistance of commercial cultivars remains the most promising solution to control this disease. In this study, an association panel of 568 modern interspecific sugarcane hybrids (Saccharum officinarum x S. spontaneum) from Réunion's breeding program was evaluated for its resistance to SOR under natural conditions of infection. Two genome-wide association studies (GWAS) were conducted between disease reactions and 183,842 single nucleotide polymorphism (SNP) markers obtained by targeted genotyping-by-sequencing. Five resistance quantitative trait loci (QTLs), named Oru1, Oru2, Oru3, Oru4 and Oru5, were identified using a single-locus GWAS (SL-GWAS). These five QTLs all originated from the species S. spontaneum. A multi-locus GWAS (ML-GWAS) uncovered an additional but less significant resistance QTL named Oru6, which originated from S. officinarum. All six QTLs had a moderate to major phenotypic effect on disease resistance. Prediction accuracy estimated with linear regression models based on each of the five QTLs identified by SL-GWAS was between 0.16-0.41. Altogether, these five QTLs provided a relatively high prediction accuracy of 0.60. In comparison, accuracies obtained with six genome-wide prediction models (i.e., GBLUP, Bayes-A, Bayes-B, Bayes-C, Bayesian Lasso and RKHS) reached only 0.65. The good prediction accuracy of disease resistance provided by the QTLs and the predominant S. spontaneum origin of their resistance alleles pave the way for effective marker-assisted breeding strategies.

CRISPR technology towards genome editing of the perennial and semi-perennial crops citrus, coffee and sugarcane

Gene editing technologies have opened up the possibility of manipulating the genome of any organism in a predicted way. CRISPR technology is the most used genome editing tool and, in agriculture, it has allowed the expansion of possibilities in plant biotechnology, such as gene knockout or knock-in, transcriptional regulation, epigenetic modification, base editing, RNA editing, prime editing, and nucleic acid probing or detection. This technology mostly depends on in vitro tissue culture and genetic transformation/transfection protocols, which sometimes become the major challenges for its application in different crops. Agrobacterium-mediated transformation, biolistics, plasmid or RNP (ribonucleoprotein) transfection of protoplasts are some of the commonly used CRISPR delivery methods, but they depend on the genotype and target gene for efficient editing. The choice of the CRISPR system (Cas9, Cas12), CRISPR mechanism (plasmid or RNP) and transfection technique (Agrobacterium spp., PEG solution, lipofection) directly impacts the transformation efficiency and/or editing rate. Besides, CRISPR/Cas technology has made countries rethink regulatory frameworks concerning genetically modified organisms and flexibilize regulatory obstacles for edited plants. Here we present an overview of the state-of-the-art of CRISPR technology applied to three important crops worldwide (citrus, coffee and sugarcane), considering the biological, methodological, and regulatory aspects of its application. In addition, we provide perspectives on recently developed CRISPR tools and promising applications for each of these crops, thus highlighting the usefulness of gene editing to develop novel cultivars.

Assessing the genetic integrity of sugarcane germplasm in the USDA-ARS National Plant Germplasm System collection using single-dose SNP markers

The World Collection of Sugarcane and Related Grasses, maintained at the USDA-ARS in Miami, FL, is one of the largest sugarcane germplasm repositories in the world. However, the genetic integrity of the Saccharum spp. germplasm in this collection has not been fully analyzed. In this study, we employed a single-dose SNP panel to genotype 901 sugarcane accessions, representing six Saccharum species and various hybrids. Our analysis uncovered a high rate of clone mislabeling in the collection. Specifically, we identified 86 groups of duplicates, characterized by identical SNP genotypes, which encompassed 211 accessions (23% of the total clones), while 135 groups, constituting 471 clones (52% of the total), exhibited near-identical genotypes. In addition, twenty-seven homonymous groups were detected, which shared the same clone name but differed in SNP genotypes. Hierarchical analysis of population structure partitioned the Saccharum germplasm into five clusters, corresponding to S. barberi, S. sinense, S. officinarum, S. spontaneum and S. robustum/S. edule. An assignment test, based on the five Saccharum species, enabled correcting 141 instances of mislabeled species memberships and inaccuracies. Moreover, we clarified the species membership and parentage of 298 clones that had ambiguous passport records (e.g., 'Saccharum spp', 'unknown', and 'hybrid'). Population structure and genetic diversity in these five species were further supported by Principal Coordinate Analysis and neighbor-joining clustering analysis. Analysis of Molecular Variance revealed that within-species genetic variations accounted for 85% of the total molecular variance, with the remaining 15% attributed to among-species genetic variations. The single-dose SNP markers developed in this study offer a robust tool for characterizing sugarcane germplasm worldwide. These findings have important implications for sugarcane genebank management, germplasm exchange, and crop genetic improvement.

Genetic diversity and population structure of sugarcane introgressed hybrids by SSR markers

One hundred and seventy two introgressed hybrids from four different introgression groups and in three nobilized stages (F1, BC1,and BC2) involving two genera viz., Sugarcane (Saccharum spp) and Erianthusl (E.arundinaceus),and Coimbatore (Commercials) clones were studied for yield and quality traits, and as well with 30 SSR markers (simple sequence repeats) for diversity and population structure. The yield and quality components of the introgressed hybrids from improved S. spontaneum, improved S. robustum, S. barberi (Pathri) and E. arundinaceus (IK-76-76) showed significant variation. The traits viz., cane height, single cane weight and cane yield showed the highest level of phenotypic and genotypic variation. 456 polymorphic alleles in all with values ranging from 6 (SMC 863-GC) to 24 (NKS 43) with an average of 15 alleles per primer were studied. Polymorphism information content (PIC) generated by 30 SSR primers varied from 0.57 (NKS 1) to 0.94 (NKS 34) and the highest average PIC (0.93) was observed in hybrids of S. spontaneum (SSH). The highest number of observed and effective alleles, and Shannon's index were also registered by introgressed hybrids of SSH and the lowest variability was observed with hybrids of Co canes x S. barberi. SSR-based population genetic measures identified three mating groups involving improved S.officinarum, improved S.robustum and Coimbatore canes with maximum variability. Structure analysis identified 51 introgressed hybrids with abundant genetic variation and low gene flow as a new germplasm source for future sugarcane breeding programmes. Fifty-nine rare alleles were detected, and structure analysis revealed a low degree of admixture between the generated hybrids. Low gene flow estimates from group I (Nm: 1.052, where mostly F1 and BC1 hybrids were grouped together) explained that the maximum variability from parental clones was inherited until the BC1 generation, and high values of Nm in group III also indicated that most of the wild alleles were fixed in the BC1 population. Structure and cluster analysis (PCoA, UPGMA) explained that the grouping of introgressed hybrids is highly correlated to different stages of nobilization. UPGMA and PCoA methods grouped the individuals into three major clusters, with a clear differentiation of F1 and backcross hybrids grouped under different internal sub-clusters. The occurrence of wide genetic variability in the population and elite hybrids with a unique genetic base can be exploited in breeding programmes to develop trait specific donors for genetic enhancement and to mitigate climate change.

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

INTRODUCTION

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Integrated transcriptomic and proteomic analyses of two sugarcane (Saccharum officinarum Linn.) varieties differing in their lodging tolerance

BACKGROUND

Lodging seriously affects sugarcane stem growth and sugar accumulation, reduces sugarcane yield and sucrose content, and impedes mechanization. However, the molecular mechanisms underlying sugarcane lodging tolerance remain unclear. In this study, comprehensive transcriptomic and proteomic analyses were performed to explore the differential genetic regulatory mechanisms between upright (GT42) and lodged (GF98-296) sugarcane varieties.

RESULTS

The stain test showed that GT42 had more lignin and vascular bundles in the stem than GF98-296. The gene expression analysis revealed that the genes that were differentially expressed between the two varieties were mainly involved in the phenylpropanoid pathway at the growth stage. The protein expression analysis indicated that the proteins that were differentially expressed between the two varieties were related to the synthesis of secondary metabolites, the process of endocytosis, and the formation of aminoacyl-tRNA. Time-series analysis revealed variations in differential gene expression patterns between the two varieties, whereas significant protein expression trends in the two varieties were largely consistent, except for one profile. The expression of CYP84A, 4CL, and CAD from the key phenylpropanoid biosynthetic pathway was enhanced in GT42 at stage 2 but suppressed in GF98-296 at the growth stage. Furthermore, the expression of SDT1 in the nicotinate and nicotinamide metabolism was enhanced in GT42 cells but suppressed in GF98-296 cells at the growth stage.

CONCLUSION

Our findings provide reference data for mining lodging tolerance-related genes that are expected to facilitate the selective breeding of sugarcane varieties with excellent lodging tolerance.

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.

Development of simple sequence repeat markers for sugarcane from data mining of expressed sequence tags

Sugarcane (Saccharum spp. hybrids) is a worldwide acclaimed important agricultural crop used primarily for sugar production and biofuel. Sugarcane's genetic complexity, aneuploidy, and extreme heterozygosity make it a challenging crop in developing improved varieties. The molecular breeding programs promise to develop nutritionally improved varieties for both direct consumption and commercial application. Therefore, to address these challenges, the development of simple sequence repeats (SSRs) has been proven to be a powerful molecular tool in sugarcane. This study involved the collection of 285216 expressed sequence tags (ESTs) from sugarcane, resulting in 23666 unigenes, including 4547 contigs. Our analysis identified 4120 unigenes containing a total of 4960 SSRs, with the most abundant repeat types being monomeric (44.33%), dimeric (13.10%), and trimeric (39.68%). We further chose 173 primers to analyze the banding pattern in 10 sugarcane accessions by PAGE analysis. Additionally, functional annotation analysis showed that 71.07%, 53.6%, and 10.3% unigenes were annotated by Uniport, GO, and KEGG, respectively. GO annotations and KEGG pathways were distributed across three functional categories: molecular (46.46%), cellular (33.94%), and biological pathways (19.6%). The cluster analysis indicated the formation of four distinct clusters among selected sugarcane accessions, with maximum genetic distance observed among the varieties. We believe that these EST-SSR markers will serve as valuable references for future genetic characterization, species identification, and breeding efforts in sugarcane.

All nonhomologous chromosomes and rearrangements in Saccharum officinarum × Saccharum spontaneum allopolyploids identified by oligo-based painting

Modern sugarcane cultivars (Saccharum spp., 2n = 100~120) are complex polyploids primarily derived from interspecific hybridization between S. officinarum and S. spontaneum. Nobilization is the theory of utilizing wild germplasm in sugarcane breeding, and is the foundation for utilizing S. spontaneum for stress resistance. However, the exact chromosomal transmission remains elusive due to a lack of chromosome-specific markers. Here, we applied chromosome-specific oligonucleotide (oligo)-based probes for identifying chromosomes 1-10 of the F1 hybrids between S. officinarum and S. spontaneum. Then, S. spontaneum-specific repetitive DNA probes were used to distinguish S. spontaneum in these hybrids. This oligo- fluorescence in situ hybridization (FISH) system proved to be an efficient tool for revealing individual chromosomal inheritance during nobilization. We discovered the complete doubling of S. officinarum-derived chromosomes in most F1 hybrids. Notably, we also found defective S. officinarum-derived chromosome doubling in the F1 hybrid Yacheng75-4191, which exhibited 1.5n transmission for all nonhomologous chromosomes. Altogether, these results highlight the presence of variable chromosome transmission in nobilization between S. officinarum and S. spontaneum, including 1.5n + n and 2n + n. These findings provide robust chromosome markers for in-depth studies into the molecular mechanism underlying chromosome doubling during the nobilization, as well as tracing chromosomal inheritance for sugarcane breeding.

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.

Flowering in sugarcane-insights from the grasses

Flowering is a crucial phase for angiosperms to continue their species propagation and is highly regulated. In the current review, flowering in sugarcane and the associated mechanisms are elaborately presented. In sugarcane, flowering has two effects, wherein it is a beneficial factor from the breeder's perspective and crucial for crop improvement, but commercially, it depletes the sucrose reserves from the stalks; hence, less value is assigned. Different species of Saccharum genus are spread across geographical latitudes, thereby proving their ability to grow in multiple inductive daylengths of different locations according in the habituated zone. In general, sugarcane is termed an intermediate daylength plant with quantitative short-day behaviour as it requires reduction in daylength from 12 h 55 min to 12 h or 12 h 30 min. The prime concern in sugarcane flowering is its erratic flowering nature. The transition to reproductive stage which reverts to vegetative stage if there is any deviation from ambient temperature and light is also an issue. Spatial and temporal gene expression patterns during vegetative to reproductive stage transition and after reverting to vegetative state could possibly reveal how the genetic circuits are being governed. This review will also shed a light on potential roles of genes and/or miRNAs in flowering in sugarcane. Knowledge of transcriptomic background of circadian, photoperiod, and gibberellin pathways in sugarcane will enable us to better understand of variable response in floral development.

Transcriptome dynamics provides insights into divergences of the photosynthesis pathway between Saccharum officinarum and Saccharum spontaneum

Saccharum spontaneum and Saccharum officinarum contributed to the genetic background of modern sugarcane cultivars. Saccharum spontaneum has shown a higher net photosynthetic rate and lower soluble sugar than S. officinarum. Here, we analyzed 198 RNA-sequencing samples to investigate the molecular mechanisms for the divergences of photosynthesis and sugar accumulation between the two Saccharum species. We constructed gene co-expression networks based on differentially expressed genes (DEGs) both for leaf developmental gradients and diurnal rhythm. Our results suggested that the divergence of sugar accumulation may be attributed to the enrichment of major carbohydrate metabolism and the oxidative pentose phosphate pathway. Compared with S. officinarum, S. spontaneum DEGs showed a high enrichment of photosynthesis and contained more complex regulation of photosynthesis-related genes. Noticeably, S. spontaneum lacked gene interactions with sulfur assimilation stimulated by photorespiration. In S. spontaneum, core genes related to clock and photorespiration displayed a sensitive regulation by the diurnal rhythm and phase-shift. Small subunit of Rubisco (RBCS) displayed higher expression in the source tissues of S. spontaneum. Additionally, it was more sensitive under a diurnal rhythm, and had more complex gene networks than that in S. officinarum. This indicates that the differential regulation of RBCS Rubisco contributed to photosynthesis capacity divergence in both Saccharum species.

Transcriptome analysis of sugarcane reveals rapid defense response of SES208 to Xanthomonas albilineans in early infection

BACKGROUND

Diseases are the major factor affecting the quality and yield of sugarcane during its growth and development. However, our knowledge about the factors regulating disease responses remain limited. The present study focuses on identifying genes regulating transcriptional mechanisms responsible for resistance to leaf scald caused by Xanthomonas albilineans in S. spontaneum and S. officinarum.

RESULTS

After inoculation of the two sugarcane varieties SES208 (S. spontaneum) and LA Purple (S. officinarum) with Xanthomonas albilineans, SES208 exhibited significantly greater resistance to leaf scald caused by X. albilineans than did LA Purple. Using transcriptome analysis, we identified a total of 4323 and 1755 differentially expressed genes (DEGs) in inoculated samples of SES208 and LA Purple, respectively. Significantly, 262 DEGs were specifically identified in SES208 that were enriched for KEGG pathway terms such as plant-pathogen interaction, MAPK signaling pathway, and plant hormone signal transduction. Furthermore, we built a transcriptional regulatory co-expression network that specifically identified 16 and 25 hub genes in SES208 that were enriched for putative functions in plant-pathogen interactions, MAPK signaling, and plant hormone signal transduction. All of these essential genes might be significantly involved in resistance-regulating responses in SES208 after X. albilineans inoculation. In addition, we found allele-specific expression in SES208 that was associated with the resistance phenotype of SES208 when infected by X. albilineans. After infection with X. albilineans, a great number of DEGs associated with the KEGG pathways 'phenylpropanoid biosynthesis' and 'flavonoid biosynthesis' exhibited significant expression changes in SES208 compared to LA Purple that might contribute to superior leaf scald resistance in SES208.

CONCLUSIONS

We provided the first systematical transcriptome map that the higher resistance of SES208 is associated with and elicited by the rapid activation of multiple clusters of defense response genes after infection by X. albilineans and not merely due to changes in the expression of genes generically associated with stress resistance. These results will serve as the foundation for further understanding of the molecular mechanisms of resistance against X. albilineans in S. spontaneum.

Identification of Saccharum CaM gene family and function characterization of ScCaM1 during cold and oxidant exposure in Pichia pastoris

BACKGROUND

Calmodulin (CaM) plays an essential role in binding calcium ions and mediating the interpretation of Ca²⁺ signals in plants under various stresses. However, the evolutionary relationship of CaM family proteins in Saccharum has not been elucidated.

OBJECTIVE

To deduce and explore the evolution and function of Saccharum CaM family.

METHODS

A total of 104 typical CaMs were obtained from Saccharum spontaneum and other 18 plant species. The molecular characteristics and evolution of those CaM proteins were analyzed. A typical CaM gene, ScCaM1, was subsequently cloned from sugarcane (Saccharum spp. hybrid). Its expression patterns in different tissues and under various abiotic stresses were assessed by quantitative real-time PCR. Then the green fluorescent protein was used to determine the subcellular localization of ScCaM1. Finally, the function of ScCaM1 was evaluated via heterologous yeast expression systems.

RESULTS

Three typical CaM members (SsCaM1, SsCaM2, and SsCaM3) were identified from the S. spontaneum genome database. CaMs were originated from the two last common ancestors before the origin of angiosperms. The number of CaM family members did not correlate to the genome size but correlated with allopolyploidization events. The ScCaM1 was more highly expressed in buds and roots than in other tissues. The expression patterns of ScCaM1 suggested that it was involved in responses to various abiotic stresses in sugarcane via different hormonal signaling pathways. Noteworthily, its expression levels appeared relatively stable during the cold exposure in the cold-tolerant variety but significantly suppressed in the cold-susceptible variety. Moreover, the recombinant yeast (Pichia pastoris) overexpressing ScCaM1 grew better than the wild-type yeast strain under cold and oxidative stresses. It was revealed that the ScCaM1 played a positive role in reactive oxygen species scavenging and conferred enhanced cold and oxidative stress tolerance to cells.

CONCLUSION

This study provided comprehensive information on the CaM gene family in Saccharum and would facilitate further investigation of their functional characterization.

Genome-wide identification and expression analysis of the NAC transcription factor family in Saccharum spontaneum under different stresses

The NAC (NAM, ATAF1/2, and CUC2) transcription factor family is one of the largest families unique to plants and is involved in plant growth and development, organs, morphogenesis, and stress responses. The NAC family has been identified in many plants. As the main source of resistance genes for sugarcane breeding, the NAC gene family in the wild species Saccharum spontaneum has not been systematically studied. In this study, 115 SsNAC genes were identified in the S. spontaneum genome, and these genes were heterogeneously distributed on 25 chromosomes. Phylogenetic analysis divided the SsNAC family members into 18 subgroups, and the gene structure and conserved motif analysis further supported the phylogenetic classification. Four groups of tandemly duplicated genes and nine pairs of segmentally duplicated genes were detected. The SsNAC gene has different expression patterns at different developmental stages of stems and leaves. Further qRT-PCR analysis showed that drought, low-temperature, salinity, pathogenic fungi, and other stresses as well as abscisic acid (ABA) and methyl jasmonate (MeJA) treatments significantly induced the expression of 12 SsNAC genes, indicating that these genes may play a key role in the resistance of S. spontaneum to biotic and abiotic stresses. In summary, the results from this study provide comprehensive information on the NAC transcription factor family, providing a reference for further functional studies of the SsNAC gene.

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.

Genome-wide characterization of cys-tathionine-β-synthase domain-containing proteins in sugarcane reveals their role in defense responses under multiple stressors

Cys-tathionine-β-synthase (CBS) domain-containing proteins (CDCPs) are essential for regulating plant responses to various biotic and abiotic stressors. This study describes the systematic identification and characterization of CDCP family genes in Saccharum spontaneum. A total of 95 SsCDCP genes and eight phylogenetic groups were identified that were distributed over 29 chromosomes of the AP85-441 genome. Most (78/95) SsCDCPs underwent fragment duplication events, and 64 gene pairs were located in synteny blocks. Expression profiling of nine ShCDCPs was also carried out in the Saccharum spp. cultivars ROC22 and MT11-611 that are resistant and susceptible to red stripe, respectively, in response to: (i) Infection by the bacterial pathogen Acidovorax avenue subsp. avenae (Aaa); (ii) abiotic stressors (drought and salinity); and (iii) exogenous salicylic acid (SA) treatment. Members of one gene pair (ShCBSD-PB1-5A and ShCBSD-PB1-7A-1) with a fragment duplication event acted as negative regulators in sugarcane under four stresses, as supported by the significantly decreased expression levels of ShCBSD-PB1-5A (23-83%) and ShCBSD-PB1-7A-1 (15-75%) at all-time points, suggesting that they have functional redundancy. Genes in another pair, ShCBS-4C and ShCBS-4D-1, which have a fragment duplication event, play opposing regulatory roles in sugarcane exposed to multiple stresses, particularly Aaa and NaCl treatments. ShCBS-4C expression was significantly decreased by 32-77%, but ShCBS-4D-1 expression was dramatically upregulated by 1.2-6.2-fold in response to Aaa treatment of both cultivars across all-time points. This result suggested that both genes exhibited functional divergence. Meanwhile, the expression of SsCBSDCBS-5A was significantly upregulated in ROC22 by 1.4-4.6-fold in response to the four stressors. These findings provide important clues for further elucidating the function of ShCDCP genes in sugarcane responding to a diverse range of stresses.

Isolation and sequencing of a single copy of an introgressed chromosome from a complex genome for gene and SNP identification

This manuscript describes the identification, isolation and sequencing of a single chromosome containing high value resistance genes from a complex polyploid where sequencing the whole genome is too costly. The large complex genomes of many crops constrain the use of new technologies for genome-assisted selection and genetic improvement. One method to simplify a genome is to break it into individual chromosomes by flow cytometry; however, in many crop species most chromosomes cannot be isolated individually. Flow sorting of a single copy of a chromosome has been developed in wheat, and here we demonstrate its use to identify markers of interest in an Erianthus/Sacchurum hybrid. Erianthus/Saccharum hybrids are of interest because Erianthus is known to be highly resistant to soil borne diseases which cause extensive sugarcane yield losses in Australia. Sugarcane (Saccharum) cultivars are autopolyploids with a highly complex genome and over 100 chromosomes. Flow cytometry for sugarcane, as in most crops, does not resolve individual chromosomes to a karyotype peak for sorting. To isolate a single chromosome, we used genomic in situ hybridization (GISH) to identify the flow karyotype region containing the Erianthus chromosomes, flow sorted single chromosomes from this region, PCR screened for the Erianthus chromosomes and sequenced them. One Erianthus chromosome amplified and sequenced well, and from this data we could identify 57 resistant type genes and SNPs in nearly half of these genes. We developed KASP SNP assays and demonstrated that the identified SNP markers segregated as expected in a small introgression population. The pipeline we developed here to flow sort and sequence single chromosomes could be used in any crop with a large complex genome to rapidly discover and develop markers to important loci.

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.

Quantitative Trait Loci Mapping and Development of KASP Marker Smut Screening Assay Using High-Density Genetic Map and Bulked Segregant RNA Sequencing in Sugarcane (Saccharum spp.)

Sugarcane is one of the most important industrial crops globally. It is the second largest source of bioethanol, and a major crop for biomass-derived electricity and sugar worldwide. Smut, caused by Sporisorium scitamineum, is a major sugarcane disease in many countries, and is managed by smut-resistant varieties. In China, smut remains the single largest constraint for sugarcane production, and consequently it impacts the value of sugarcane as an energy feedstock. Quantitative trait loci (QTLs) associated with smut resistance and linked diagnostic markers are valuable tools for smut resistance breeding. Here, we developed an F1 population (192 progeny) by crossing two sugarcane varieties with contrasting smut resistance and used for genome-wide single nucleotide polymorphism (SNP) discovery and mapping, using a high-throughput genotyping method called "specific locus amplified fragment sequencing (SLAF-seq) and bulked-segregant RNA sequencing (BSR-seq). SLAF-seq generated 148,500 polymorphic SNP markers. Using SNP and previously identified SSR markers, an integrated genetic map with an average 1.96 cM marker interval was produced. With this genetic map and smut resistance scores of the F1 individuals from four crop years, 21 major QTLs were mapped, with a phenotypic variance explanation (PVE) > 8.0%. Among them, 10 QTLs were stable (repeatable) with PVEs ranging from 8.0 to 81.7%. Further, four QTLs were detected based on BSR-seq analysis. aligning major QTLs with the genome of a sugarcane progenitor Saccharum spontaneum, six markers were found co-localized. Markers located in QTLs and functional annotation of BSR-seq-derived unigenes helped identify four disease resistance candidate genes located in major QTLs. 77 SNPs from major QTLs were then converted to Kompetitive Allele-Specific PCR (KASP) markers, of which five were highly significantly linked to smut resistance. The co-localized QTLs, candidate resistance genes, and KASP markers identified in this study provide practically useful tools for marker-assisted sugarcane smut resistance breeding.

Utilization of a Sugarcane100K Single Nucleotide Polymorphisms Microarray-Derived High-Density Genetic Map in Quantitative Trait Loci Mapping and Function Role Prediction of Genes Related to Chlorophyll Content in Sugarcane

Chlorophyll is the most important pigment for plant photosynthesis that plays an important role in crop growth and production. In this study, the chlorophyll content trait was explored to improve sugarcane yield. Two hundred and eighty-five F₁ progenies from the cross YT93-159 × ROC22 with significantly different chlorophyll contents were included as test materials. The chlorophyll content of the +1 leaves during elongation phase was measured using a SPAD-502 meter through a three-crop cycle (plant cane, first ratoon, and second ratoon). Linkage analysis was conducted on a high-density genetic map constructed based on the sugarcane 100K SNP chip. In addition, Fv/Fm, plant height, stalk diameter, brix data were collected on plant cane during the elongation and maturation phases. The results showed that the +1 leaf SPAD values, which can be used as an important reference to evaluate the growth potential of sugarcane, were significantly and positively correlated with the Fv/Fm during elongation phase, as well as with plant height, stalk diameter, and brix during maturity phase (P < 0.01). The broad sense heritability (H ²) of the chlorophyll content trait was 0.66 for plant cane crop, 0.67 for first ratoon crop, and 0.73 for second ratoon crop, respectively, indicating that this trait was mainly controlled by genetic factors. Thirty-one quantitative trait loci (QTL) were detected by QTL mapping. Among them, a major QTL, qCC-R1, could account for 12.95% of phenotypic variation explained (PVE), and the other 30 minor QTLs explained 2.37-7.99% PVE. Twenty candidate genes related to chlorophyll content were identified in the QTLs plus a 200-Kb extension region within either sides, of which four were homologous genes involved in the chlorophyll synthesis process and the remaining 16 played a certain role in chlorophyll catabolic pathway, chloroplast organization, or photosynthesis. These results provide a theoretical reference for analyzing the genetic mechanism of chlorophyll synthesis and subsequent improvement of photosynthetic characteristics in sugarcane.

Sugarcane Ratooning Ability: Research Status, Shortcomings, and Prospects

Sugarcane is an important sugar crop and it can be subjected to ratooning for several years. The advantages of ratooning include quality improvement, efficiency enhancement, and reduced costs and energy use. The genotype, environment, cultivation management, and harvesting technology affect the productivity and longevity of ratoon cane, with the genetic basis being the most critical factor. However, the majority of research has been focused on only limited genotypes, and a few studies have evaluated up to 100 sugarcane germplasm resources. They mainly focus on the comparison among different genotypes or among plant cane, different selection strategies for the first and second ratoon crops, together with screening indicators for the selection of stronger ratooning ability. In this paper, previous studies are reviewed in order to analyze the importance of sugarcane ratooning, the indicative traits used to evaluate ratooning ability, the major factors influencing the productivity and longevity of ratooning, the genetic basis of variation in ratooning ability, and the underlying mechanisms. Furthermore, the shortcomings of the existing research on sugarcane ratooning are highlighted. We then discuss the focus of future ratoon sugarcane research and the technical methods that will shorten the selection cycle and increase the genetic gain of ratooning ability, particularly the development of linked markers. This review is expected to provide a reference for understanding the mechanisms underlying the formation of ratooning ability and for breeding sugarcane varieties with a strong ratooning ability.

Chromosome Painting Provides Insights Into the Genome Structure and Evolution of Sugarcane

The genus Saccharum is composed of species with high polyploidy and highly varied chromosome numbers, laying a challenge for uncovering its genomic structure and evolution. We developed a chromosome 2 painting (CP2) probe by designing oligonucleotides covering chromosome 2 of Saccharum spontaneum (2n = 8x = 64). Fluorescence in situ hybridization (FISH) using this CP2 probe revealed six types of ploidies from twenty S. spontaneum clones, including 6x, 8x, 10x, 11x, 12x, and 13x clones. The finding of S. spontaneum clones with uneven of ploid suggested that certain S. spontaneum clones come from hybridization. It renews our knowledge that S. spontaneum is derived from autopolyploidization. Combined with a S. spontaneum-specific probe, chromosome 2-derived chromosome or fragments from either S. spontaneum or Saccharum officinarum can be identified in sugarcane modern cultivars. We revealed unexpected high level of interspecific recombination from introgressive S. spontaneum chromosomes (>50.0%) in cultivars ROC22 and ZZ1, indicating frequent chromosome exchange in cultivars. Intriguingly, we observed interspecific recombination recurring among either homoeologous or non-homoeologous chromosomes in sugarcane cultivars. These results demonstrated that chromosome painting FISH is a powerful tool in the genome dissection of sugarcane and provide new insights into the genome structure and evolution of the complex genus Saccharum.

Comparative phylogenetic analysis of CBL reveals the gene family evolution and functional divergence in Saccharum spontaneum

BACKGROUND

The identification and functional analysis of genes that improve tolerance to low potassium stress in S. spontaneum is crucial for breeding sugarcane cultivars with efficient potassium utilization. Calcineurin B-like (CBL) protein is a calcium sensor that interacts with specific CBL-interacting protein kinases (CIPKs) upon plants' exposure to various abiotic stresses.

RESULTS

In this study, nine CBL genes were identified from S. spontaneum. Phylogenetic analysis of 113 CBLs from 13 representative plants showed gene expansion and strong purifying selection in the CBL family. Analysis of CBL expression patterns revealed that SsCBL01 was the most commonly expressed gene in various tissues at different developmental stages. Expression analysis of SsCBLs under low K⁺ stress indicated that potassium deficiency moderately altered the transcription of SsCBLs. Subcellular localization showed that SsCBL01 is a plasma membrane protein and heterologous expression in yeast suggested that, while SsCBL01 alone could not absorb K⁺, it positively regulated K⁺ absorption mediated by the potassium transporter SsHAK1.

CONCLUSIONS

This study provided insights into the evolution of the CBL gene family and preliminarily demonstrated that the plasma membrane protein SsCBL01 was involved in the response to low K⁺ stress in S. spontaneum.

Genome-wide analysis of R2R3-MYB transcription factors family in the autopolyploid Saccharum spontaneum: an exploration of dominance expression and stress response

BACKGROUND

Sugarcane (Saccharum) is the most critical sugar crop worldwide. As one of the most enriched transcription factor families in plants, MYB genes display a great potential to contribute to sugarcane improvement by trait modification. We have identified the sugarcane MYB gene family at a whole-genome level through systematic evolution analyses and expression profiling. R2R3-MYB is a large subfamily involved in many plant-specific processes.

RESULTS

A total of 202 R2R3-MYB genes (356 alleles) were identified in the polyploid Saccharum spontaneum genomic sequence and classified into 15 subgroups by phylogenetic analysis. The sugarcane MYB family had more members by a comparative analysis in sorghum and significant advantages among most plants, especially grasses. Collinearity analysis revealed that 70% of the SsR2R3-MYB genes had experienced duplication events, logically suggesting the contributors to the MYB gene family expansion. Functional characterization was performed to identify 56 SsR2R3-MYB genes involved in various plant bioprocesses with expression profiling analysis on 60 RNA-seq databases. We identified 22 MYB genes specifically expressed in the stem, of which RT-qPCR validated MYB43, MYB53, MYB65, MYB78, and MYB99. Allelic expression dominance analysis implied the differential expression of alleles might be responsible for the high expression of MYB in the stem. MYB169, MYB181, MYB192 were identified as candidate C4 photosynthetic regulators by C4 expression pattern and robust circadian oscillations. Furthermore, stress expression analysis showed that MYB36, MYB48, MYB54, MYB61 actively responded to drought treatment; 19 and 10 MYB genes were involved in response to the sugarcane pokkah boeng and mosaic disease, respectively.

CONCLUSIONS

This is the first report on genome-wide analysis of the MYB gene family in sugarcane. SsMYBs probably played an essential role in stem development and the adaptation of various stress conditions. The results will provide detailed insights and rich resources to understand the functional diversity of MYB transcription factors and facilitate the breeding of essential traits in sugarcane.

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.

CRISPR/Cas9-Mediated Multi-Allelic Gene Targeting in Sugarcane Confers Herbicide Tolerance

Sugarcane is the source of 80% of the sugar and 26% of the bioethanol produced globally. However, its complex, highly polyploid genome (2n = 100 - 120) impedes crop improvement. Here, we report efficient and reproducible gene targeting (GT) in sugarcane, enabling precise co-editing of multiple alleles via template-mediated and homology-directed repair (HDR) of DNA double strand breaks induced by the programmable nuclease CRISPR/Cas9. The evaluation of 146 independently transformed plants from five independent experiments revealed a targeted nucleotide replacement that resulted in both targeted amino acid substitutions W574L and S653I in the acetolactate synthase (ALS) in 11 lines in addition to single, targeted amino acid substitutions W574L or S653I in 25 or 18 lines, respectively. Co-editing of up to three ALS copies/alleles that confer herbicide tolerance was confirmed by Sanger sequencing of cloned long polymerase chain reaction (PCR) amplicons. This work will enable crop improvement by conversion of inferior alleles to superior alleles through targeted nucleotide substitutions.

Detection of a major QTL related to smut disease resistance inherited from a Japanese wild sugarcane using GRAS-Di technology

Smut disease of sugarcane causes considerable yield losses and the use of resistant varieties is the best control practice. Our group identified a Japanese wild sugarcane with highly smut disease resistance named 'Iriomote8'. In this study, we conducted QTL analysis for smut disease resistance using a mapping population derived from a resistant variety 'Yaenoushie', in which resistance is inherited from 'Iriomote8'. We identified 4813 non-redundant markers using GRAS-Di technology and developed a linkage map of mapping parents. We evaluated smut disease resistance of the mapping population by the inoculation test. Consequently, a large number of clones did not show the disease symptoms and the distribution of smut disease incidence tended to be "L shaped". Composite interval mapping detected an identical QTL for indices of smut disease incidence with a markedly high LOD score (26.6~45.6) at the end of linkage group 8 of 'Yaenoushie'. This QTL explained approximately 50% of the cases of smut disease incidence. In the mapping population, there were no correlations between the indices of smut disease incidence and other agronomic traits. In conclusion, this QTL could be used for marker-assisted selection to significantly improve smut disease resistance without negative effects on other agronomic traits.

Chromosome behavior during meiosis in pollen mother cells from Saccharum officinarum × Erianthus arundinaceus F1 hybrids

BACKGROUND

In recent years, sugarcane has attracted increasing attention as an energy crop. Wild resources are widely used to improve the narrow genetic base of sugarcane. However, the infertility of F₁ hybrids between Saccharum officinarum (S. officinarum) and Erianthus arundinaceus (E. arundinaceus) has hindered sugarcane breeding efforts. To discover the cause of this infertility, we studied the hybridization process from a cytological perspective.

RESULTS

We examined the meiotic process of pollen mother cells (PMCs) in three F₁ hybrids between S. officinarum and E. arundinaceus. Cytological analysis showed that the male parents, Hainan 92-77 and Hainan 92-105, had normal meiosis. However, the meiosis process in F₁ hybrids showed various abnormal phenomena, including lagging chromosomes, micronuclei, uneven segregation, chromosome bridges, and inability to form cell plates. Genomic in situ hybridization (GISH) showed unequal chromatin distribution during cell division. Interestingly, 96.70% of lagging chromosomes were from E. arundinaceus. Furthermore, fluorescence in situ hybridization (FISH) was performed using 45S rDNA and 5S rDNA as probes. Either 45S rDNA or 5S rDNA sites were lost during abnormal meiosis, and results of unequal chromosomal separation were also clearly observed in tetrads.

CONCLUSIONS

Using cytogenetic analysis, a large number of meiotic abnormalities were observed in F₁. GISH further confirmed that 96.70% of the lagging chromosomes were from E. arundinaceus. Chromosome loss was found by further investigation of repeat sequences. Our findings provide insight into sugarcane chromosome inheritance to aid innovation and utilization in sugarcane germplasm resources.

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

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

Genome-Wide Characterization of Lectin Receptor Kinases in Saccharum spontaneum L. and Their Responses to Stagonospora tainanensis Infection

Sugarcane is an important sugar and bioenergy ethanol crop, and the hyperploidy has led to stagnant progress in sugarcane genome decipherment, which also hindered the genome-wide analyses of versatile lectin receptor kinases (LecRKs). The published genome of Saccharum spontaneum, one of the two sugarcane ancestor species, enables us to study the characterization of LecRKs and their responses to sugarcane leaf blight (SLB) triggered by Stagonospora tainanensis. A total of 429 allelic and non-allelic LecRKs, which were classified into evolved independently three types according to signal domains and phylogeny, were identified based on the genome. Regarding those closely related LecRKs in the phylogenetic tree, their motifs and exon architectures of representative L- and G-types were similar or identical. LecRKs showed an unequal distribution on chromosomes and more G-type tandem repeats may come from the gene expansion. Comparing the differentially expressed LecRKs (DELs) in response to SLB in sugarcane hybrid and ancestor species S. spontaneum, we found that the DEL number in the shared gene sets was highly variable among each sugarcane accession, which indicated that the expression dynamics of LecRKs in response to SLB were quite different between hybrids and particularly between sugarcane hybrid and S. spontaneum. In addition, C-type LecRKs may participate in metabolic processes of plant–pathogen interaction, mainly including pathogenicity and plant resistance, indicating their putative roles in sugarcane responses to SLB infection. The present study provides a basic reference and global insight into the further study and utilization of LecRKs in plants.

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.

Genome-wide analysis of mitogen-activated protein (MAP) kinase gene family expression in response to biotic and abiotic stresses in sugarcane

To systematically analyze mitogen-activated protein (MAP) kinase gene families and their expression profiles in sugarcane (Saccharum spp. hybrids; Sh) under diverse biotic and abiotic stresses, we identified 15 ShMAPKs, 6 ShMAPKKs and 16 ShMAPKKKs genes in the sugarcane cultivar R570 genome. These were also confirmed in one S. spontaneum genome and two transcriptome datasets of sugarcane trigged by Acidovorax avenae subsp. avenae (Aaa) and Xanthomonas albilineans (Xa) infections. Phylogenetic analysis revealed that four subgroups were present in each ShMAPK and ShMAPKK family and three sub-families (RAF, MEKK and ZIK) presented in the ShMAPKKK family. Conserved protein motif and gene structure analyses supported the evolutionary relationships of the three families inferred from the phylogenetic analysis. All of the ShMAPK, ShMAPKK and ShMAPKKK genes identified in Saccharum spp. R570 were distributed on chromosomes 1-7 and 9-10. RNA-seq and qRT-PCR analyses indicated that ShMAPK07 and ShMAPKKK02 were defense-responsive genes in sugarcane challenged by both Aaa and Xa stimuli, while some genes were upregulated specifically by Aaa and Xa infection. Additionally, ShMAPK05 acted as a negative regulator under drought and salinity stress, but served as a positive regulator under salicylic acid (SA) treatment. ShMAPK07 plays a positive role under drought stress, but a negative role under SA treatment. ShMAPKKK01 was negatively modulated by both salinity stress and SA treatment, whereas ShMAPKKK06 was positively regulated by both of the two stress stimuli. Our results suggest that members of MAPK cascade gene families regulate adverse stress responses through multiple signal transduction pathways in sugarcane.

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.

Expression profiling of MADS-box gene family revealed its role in vegetative development and stem ripening in S. spontaneum

Sugarcane is the most important sugar and biofuel crop. MADS-box genes encode transcription factors that are involved in developmental control and signal transduction in plants. Systematic analyses of MADS-box genes have been reported in many plant species, but its identification and characterization were not possible until a reference genome of autotetraploid wild type sugarcane specie, Saccharum spontaneum is available recently. We identified 182 MADS-box sequences in the S. spontaneum genome, which were annotated into 63 genes, including 6 (9.5%) genes with four alleles, 21 (33.3%) with three, 29 (46%) with two, 7 (11.1%) with one allele. Paralogs (tandem duplication and disperse duplicated) were also identified and characterized. These MADS-box genes were divided into two groups; Type-I (21 Mα, 4 Mβ, 4 Mγ) and Type-II (32 MIKCc, 2 MIKC*) through phylogenetic analysis with orthologs in Arabidopsis and sorghum. Structural diversity and distribution of motifs were studied in detail. Chromosomal localizations revealed that S. spontaneum MADS-box genes were randomly distributed across eight homologous chromosome groups. The expression profiles of these MADS-box genes were analyzed in leaves, roots, stem sections and after hormones treatment. Important alleles based on promoter analysis and expression variations were dissected. qRT-PCR analysis was performed to verify the expression pattern of pivotal S. spontaneum MADS-box genes and suggested that flower timing genes (SOC1 and SVP) may regulate vegetative development.

Sugarcane genome architecture decrypted with chromosome-specific oligo probes

Sugarcane (Saccharum spp.) is probably the crop with the most complex genome. Modern cultivars (2n = 100-120) are highly polyploids and aneuploids derived from interspecific hybridization between Saccharum officinarum (2n = 80) and Saccharum spontaneum (2n = 40-128). Chromosome-specific oligonucleotide probes were used in combination with genomic in situ hybridization to analyze the genome architecture of modern cultivars and representatives of their parental species. The results validated a basic chromosome number of x = 10 for S. officinarum. In S. spontaneum, rearrangements occurred from a basic chromosome of x = 10, probably in the Northern part of India, in two steps leading to x = 9 and then x = 8. Each step involved three chromosomes that were rearranged into two. Further polyploidization led to the wide geographical extension of clones with x = 8. We showed that the S. spontaneum contribution to modern cultivars originated from cytotypes with x = 8 and varied in proportion between cultivars (13-20%). Modern cultivars had mainly 12 copies for each of the first four basic chromosomes, and a more variable number for those basic chromosomes whose structure differs between the two parental species. One-four of these copies corresponded to entire S. spontaneum chromosomes or interspecific recombinant chromosomes. In addition, a few inter-chromosome translocations were revealed. The new information and cytogenetic tools described in this study substantially improve our understanding of the extreme level of complexity of modern sugarcane cultivar genomes.

Genome-wide identification and characterization of DCL, AGO and RDR gene families in Saccharum spontaneum

RNA silencing is a conserved mechanism in eukaryotic organisms to regulate gene expression. Argonaute (AGO), Dicer-like (DCL) and RNA-dependent RNA polymerase (RDR) proteins are critical components of RNA silencing, but how these gene families’ functions in sugarcane were largely unknown. Most stress-resistance genes in modern sugarcane cultivars (Saccharum spp.) were originated from wild species of Saccharum, for example S. spontaneum. Here, we used genome-wide analysis and a phylogenetic approach to identify four DCL, 21 AGO and 11 RDR genes in the S. spontaneum genome (termed SsDCL, SsAGO and SsRDR, respectively). Several genes, particularly some of the SsAGOs, appeared to have undergone tandem or segmental duplications events. RNA-sequencing data revealed that four SsAGO genes (SsAGO18c, SsAGO18b, SsAGO10e and SsAGO6b) and three SsRDR genes (SsRDR2b, SsRDR2d and SsRDR3) tended to have preferential expression in stem tissue, while SsRDR5 was preferentially expressed in leaves. qRT-PCR analysis showed that SsAGO10c, SsDCL2 and SsRDR6b expressions were strongly upregulated, whereas that of SsAGO18b, SsRDR1a, SsRDR2b/2d and SsRDR5 was significantly depressed in S. spontaneum plants exposed to PEG-induced dehydration stress or infected with Xanthomonas albilineans, causal agent of leaf scald disease of sugarcane, suggesting that these genes play important roles in responses of S. spontaneum to biotic and abiotic stresses.

Differential expression between drought-tolerant and drought-sensitive sugarcane under mild and moderate water stress as revealed by a comparative analysis of leaf transcriptome

Sugarcane contributes 80% of global sugar production and to bioethanol generation for the bioenergy industry. Its productivity is threatened by drought that can cause up to 60% yield loss. This study used RNA-Seq to gain a better understanding of the underlying mechanism by which drought-tolerant sugarcane copes with water stress. We compared gene expression in KPS01-12 (drought-tolerant genotype) and UT12 (drought-sensitive genotype) that have significantly different yield loss rates under drought conditions. We treated KPS01-12 and UT12 with mild and moderate water stress and found differentially expressed genes in various biological processes. KPS01-12 had higher expression of genes that were involved in water retention, antioxidant secondary metabolite biosynthesis, and oxidative and osmotic stress response than UT12. In contrast, the sensitive genotype had more down-regulated genes that were involved in photosynthesis, carbon fixation and Calvin cycle than the tolerant genotype. Our obtained expression profiles suggest that the tolerant sugarcane has a more effective genetic response than the sensitive genotype at the initiation of drought stress. The knowledge gained from this study may be applied in breeding programs to improve sugarcane production in drought conditions.

Isolation, characterization and expression analysis of stress responsive plant nuclear transcriptional factor subunit (NF-YB2) from commercial Saccharum hybrid and wild relative Erianthus arundinaceus

Plant nuclear factor (NF-Y) is a transcription activating factor, consisting of three subunits, and plays a key regulatory role in many stress-responsive mechanisms including drought and salinity stresses. NF-Ys function both as complex and individual subunits. Considering the importance of sugarcane as a commercial crop with high socio-economic importance and the crop being affected mostly by water deficit stress and salinity stress causing significant yield loss, nuclear transcriptional factor NF-YB2 was focused in this study. Plant nuclear factor subunit B2 from Erianthus arundinaceus (EaNF-YB2), a wild relative of sugarcane which is known for its drought and salinity stress tolerance, and commercial Saccharum hybrid Co 86032 (ShNF-YB2) was isolated and characterized. Both EaNF-YB2 and ShNF-YB2 genes are 543 bp long that encodes for a polypeptide of 180 amino acid residues. Comparison of EaNF-YB2 and ShNF-YB2 gene sequences revealed nucleotide substitutions at nine positions corresponding to three synonymous and six nonsynonymous amino acid substitutions that resulted in variations in physiochemical properties. However, multiple sequence alignment (MSA) of NF-YB2 proteins showed conservation of functionally important amino acid residues. In silico analysis revealed NF-YB2 to be a hydrophilic and intracellular protein, and EaNF-YB2 is thermally more stable than that of ShNF-YB2. Phylogenetic analysis suggested the lower rate of evolution of NF-YB2. Subcellular localization in sugarcane callus revealed NF-YB2 localization at nucleus that further evidenced it to be a transcription activation factor. Comparative RT-qPCR experiments showed a significantly higher level of NF-YB2 expression in E. arundinaceus when compared to that in the commercial Saccharum hybrid Co 86032 under drought and salinity stresses. Hence, EaNF-YB2 could be an ideal candidate gene, and its overexpression in sugarcane through genetic engineering approach might enhance tolerance to drought and salinity stresses.

Genetic analysis of agronomic traits in elite sugarcane (Saccharum spp.) germplasm

Sugarcane (Saccharum spp.) is an important economic crop, supplying up to 80% of the table sugar and ~60% of bio-ethanol worldwide. Due to population growth and dwindling fossil-fuel reserves, the demand for sugar and bio-ethanol requires significant improvement in sugarcane production. Breeding sugarcane cultivars with high-performance agronomic traits is undoubtedly the most efficient way to achieve this goal. Therefore, evaluating agronomic traits and dissecting underlying loci are critically important for this aim steps in providing genetic resources and molecular markers for selection. In this study, we assembled a diversity panel of 236 elite sugarcane germplasms originally collected from 12 countries. We evaluated 28 agronomic traits in the diversity panel with three replicates. The diversity panel was genotyped using amplified fragment length polymorphism markers, and a total of 1,359 markers were generated. Through the genome-wide association study, we identified three markers significantly associated with three traits evaluated at a stringent threshold (P < 0.05 after Bonferroni correction). The genotypes of the three associated markers grouped respective trait values into two distinct groups, supporting the reliability of these markers for breeding selection. Our study provides putative molecular markers linked to agronomic traits for breeding robust sugarcane cultivars. Additionally, this study emphasized the importance of sugarcane germplasm introduced from other countries and suggested that the use of these germplasms in breeding programs depends on local industrial needs.

Genome-wide association study of multiple yield traits in a diversity panel of polyploid sugarcane (Saccharum spp.)

Sugarcane (Saccharum spp.) is an important economic crop, contributing up to 80% of sugar and approximately 60% of biofuel globally. To meet the increased demand for sugar and biofuel supplies, it is critical to breed sugarcane cultivars with robust performance in yield traits. Therefore, dissection of causal DNA sequence variants is of great importance, as it provides genetic resources and fundamental information for crop improvement. In this study, we analyzed nine yield traits in a sugarcane diversity panel consisting of 308 accessions primarily selected from the World Collection of Sugarcane and Related Grasses. By genotyping the diversity panel via target enrichment sequencing, we identified a large number of sequence variants. Genome-wide association studies between the markers and traits were conducted, taking dosages and gene actions into consideration. In total, 217 nonredundant markers and 225 candidate genes were identified to be significantly associated with the yield traits, which can serve as a comprehensive genetic resource database for future gene identification, characterization, and selection for sugarcane improvement. We further investigated runs of homozygosity (ROH) in the sugarcane diversity panel. We characterized 282 ROHs and found that the occurrence of ROHs in the genome were nonrandom and probably under selection. The ROHs were associated with total weight and dry weight, and high ROHs resulted in a decrease in the two traits. This study suggests that genomic inbreeding has led to negative impacts on sugarcane yield.

Identification of Differentially Expressed Proteins in Sugarcane in Response to Infection by Xanthomonas albilineans Using iTRAQ Quantitative Proteomics

Sugarcane can suffer severe yield losses when affected by leaf scald, a disease caused by Xanthomonas albilineans. This bacterial pathogen colonizes the vascular system of sugarcane, which can result in reduced plant growth and plant death. In order to better understand the molecular mechanisms involved in the resistance of sugarcane to leaf scald, a comparative proteomic study was performed with two sugarcane cultivars inoculated with X. albilineans: one resistant (LCP 85-384) and one susceptible (ROC20) to leaf scald. The iTRAQ (isobaric tags for relative and absolute quantification) approach at 0 and 48 h post-inoculation (hpi) was used to identify and annotate differentially expressed proteins (DEPs). A total of 4295 proteins were associated with 1099 gene ontology (GO) terms by GO analysis. Among those, 285 were DEPs during X. albilineans infection in cultivars LCP 85-384 and ROC20. One hundred seventy-two DEPs were identified in resistant cultivar LCP 85-384, and 113 of these proteins were upregulated and 59 were downregulated. One hundred ninety-two DEPs were found in susceptible cultivar ROC20 and half of these (92) were upregulated, whereas the other half corresponded to downregulated proteins. The significantly upregulated DEPs in LCP 85-384 were involved in metabolic pathways, the biosynthesis of secondary metabolites, and the phenylpropanoid biosynthesis pathway. Additionally, the expression of seven candidate genes related to photosynthesis and glycolytic pathways, plant innate immune system, glycosylation process, plant cytochrome P450, and non-specific lipid transfer protein was verified based on transcription levels in sugarcane during infection by X. albilineans. Our findings shed new light on the differential expression of proteins in sugarcane cultivars in response to infection by X. albilineans. The identification of these genes provides important information for sugarcane variety improvement programs using molecular breeding strategies.

Alternative Splicing of Circadian Clock Genes Correlates With Temperature in Field-Grown Sugarcane

Alternative Splicing (AS) is a mechanism that generates different mature transcripts from precursor mRNAs (pre-mRNAs) of the same gene. In plants, a wide range of physiological and metabolic events are related to AS, as well as fast responses to changes in temperature. AS is present in around 60% of intron-containing genes in Arabidopsis, 46% in rice, and 38% in maize and it is widespread among the circadian clock genes. Little is known about how AS influences the circadian clock of C4 plants, like commercial sugarcane, a C4 crop with a complex hybrid genome. This work aims to test if the daily dynamics of AS forms of circadian clock genes are regulated by environmental factors, such as temperature, in the field. A systematic search for AS in five sugarcane clock genes, ScLHY, ScPRR37, ScPRR73, ScPRR95, and ScTOC1 using different organs of sugarcane sampled during winter, with 4 months old plants, and during summer, with 9 months old plants, revealed temperature- and organ-dependent expression of at least one alternatively spliced isoform in all genes. Expression of AS isoforms varied according to the season. Our results suggest that AS events in circadian clock genes are correlated with temperature.

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.