Sugarcane genome architecture decrypted with chromosome-specific oligo probes

Identification and validation of intergeneric hybrids between Saccharum officinarum and Erianthus rockii using molecular and cytogenetic tools

Sugarcane (Saccharum spp.), a vital crop for sugar and bioenergy production, faces challenges in breeding due to its narrow genetic base and susceptibility to environmental stresses. To enhance genetic diversity, distant hybridization with wild relatives such as Erianthus rockii, known for its drought tolerance and strong ratooning ability, offers a promising strategy. However, identifying true intergeneric hybrids remains a critical challenge. This study identified true intergeneric hybrids from crosses between S. officinarum and E. rockii using tetra-primer ARMS-PCR and HRM curve analysis targeting SNPs in the nrDNA-ITS region, with genomic in situ hybridization (GISH) validating chromosome composition. The results confirmed 13 true hybrids among 16 progeny, showing a chromosome inheritance pattern of 40 chromosomes derived from S. officinarum and 15 from E. rockii. Hybrid plants, while exhibiting traits from both parents, displayed inferior yield and sugar content compared to the maternal parent, indicating the necessity for further backcrossing to improve agronomic performance. This study pioneers the application of ARMS-PCR and HRM in identifying E. rockii hybrids and underscores their potential in advancing sugarcane breeding by facilitating the introgression of beneficial traits from wild relatives.

Gibberellin biosynthesis gene ScGA20 oxidase enhances sugarcane growth by modulating genes associated with phytohormone and growth processes

Sugarcane is a globally significant crop for sugar and energy production, yet its breeding potential is limited by its complex genetic background and restricted genetic diversity. Molecular breeding presents a promising approach to boost sugarcane productivity. Gibberellins (GA) plays a critical role in plant growth and development, with GA20-oxidase being the most crucial enzyme in GA biosynthesis. In this study, we isolated the sugarcane gene encoding GA20-oxidase (ScGA20ox, OR283803) enzyme, a 43.88 kDa hydrophilic protein, and introduced it into sugarcane variety GT42 via biolistics transformation. Transgenic sugarcane expressing Ubi-driven ScGA20ox exhibited elevated GA levels and accelerated growth. Transcriptome analysis revealed that genes associated with plant hormone metabolism and growth-related pathways, including photosynthesis, plant hormone signal transduction, and starch and sucrose metabolism were involved in the GA20ox-mediated growth in sugarcane. This work underscores the potential of GA20ox for enhancing sugarcane yields through transgenic approaches, while advancing our understanding of GA's regulatory role in the growth of this crucial sugar crop.

Functional divergence of CAD-like family genes in Saccharum complex under biotic and abiotic stress

KEY MESSAGE

A total of 54 genes of membrane attack complex/perforin (MACPF) superfamily were identified in Saccharum complex and function divergence among SsaCAD-like genes were present in plant against stressors. The membrane attack complex/perforin (MACPF) superfamily belongs to pore-forming proteins involving in innate and adaptive immunity in eukaryotes. The constitutively activated cell death (CAD) proteins contained the MACPF domain participate in plant defense responses under adverse conditions. However, the characteristics and functions of CAD-like genes in sugarcane are still poorly understood. In this study, 54 CAD-like genes were identified in three genomes from Saccharum complex including two clones (Np-X and AP85-441) of Saccharum spontaneum and a clone (Yunnan2009-3) of Erianthus rufipilus. All CAD-like genes were categorized into five phylogenetic groups (I-V). Various cis-acting elements related to stress responses, such as phytohormone response elements, were found in promoter regions. Transcriptome and RT-qPCR analysis demonstrated these genes possessing diverse expression profiles. The SsaCAD1-like1 gene was upregulated in sugarcane cultivars after cold treatment and infection by Xanthomonas albilineans (Xa) causing leaf scald. Meanwhile, this gene was downregulated under drought and ABA treatments but was upregulated and then downregulated across time-points of SA treatment. The SsaCAD4-like1 gene was downregulated under five abiotic stressors. Expression levels of two alleles (SsaCAD2-like1/2) were significantly decreased under all abiotic stressors except for salinity treatment. Similar expression patterns of three alleles (SsaCAD3-like1/2/4) were found under abiotic stress. The SsaCAD1-like1 exhibited a negative role but SsaCAD3-like2/4 acted as positive roles in transgenic Arabidopsis lines against bacterial pathogen infection. Our results provide novel gene resources for developing disease-resistant cultivars in sugarcane.

Exploiting historical agronomic data to develop genomic prediction strategies for early clonal selection in the Louisiana sugarcane variety development program

Genomic selection can enhance the rate of genetic gain of cane and sucrose yield in sugarcane (Saccharum L.), an important industrial crop worldwide. We assessed the predictive ability (PA) for six traits, such as theoretical recoverable sugar (TRS), number of stalks (NS), stalk weight (SW), cane yield (CY), sugar yield (SY), and fiber content (Fiber) using 20,451 single nucleotide polymorphisms (SNPs) with 22 statistical models based on the genomic estimated breeding values of 567 genotypes within and across five stages of the Louisiana sugarcane breeding program. TRS and SW with high heritability showed higher PA compared to other traits, while NS had the lowest. Machine learning (ML) methods, such as random forest and support vector machine (SVM), outperformed others in predicting traits with low heritability. ML methods predicted TRS and SY with the highest accuracy in cross-stage predictions, while Bayesian models predicted NS and CY with the highest accuracy. Extended genomic best linear unbiased prediction models accounting for dominance and epistasis effects showed a slight improvement in PA for a few traits. When both NS and TRS, which can be available as early as stage 2, were considered in a multi-trait selection model, the PA for SY in stage 5 could increase up to 0.66 compared to 0.30 with a single-trait model. Marker density assessment suggested 9091 SNPs were sufficient for optimal PA of all traits. The study demonstrated the potential of using historical data to devise genomic prediction strategies for clonal selection early in sugarcane breeding programs.

Oligo-FISH-Based Analysis of the Mechanisms Underlying Chromosome Number Variation in Saccharum spontaneum

Interspecific hybridization serves as a crucial strategy for innovating sugarcane germplasms. Currently, nearly all modern sugarcane varieties that incorporate genetic material are derived from Saccharum spontaneum. The number of chromosomes in S. spontaneum ranges from 40 to 128, contributing significantly to the diversity of its genetic resources. However, the genetic mechanisms driving chromosome number variation in S. spontaneum remain to be fully elucidated. Here, oligonucleotide fluorescence in situ hybridization (Oligo-FISH) was conducted to identify individual chromosomes and explore chromosome transmission during the intraspecific hybridization of S. spontaneum. The results indicate that from the progenies generated from S. spontaneum Yunnan2017-22 (2n = 8x = 64) and Yunnan82-1 (2n = 8x = 64) emerged two distinct karyotypes, 2n = 12x = 96 (A1) and 2n = 8x = 64 (A2, A33-1, A18). This implies that the chromosome inheritances were 2n + n and n + n in the progenies. However, self-pollinated samples of A1 (2n = 12x = 96) produced normal offspring C1 (2n = 94) and C2 (2n = 96). The 2n + n inheritance pattern did not continue. In another cross, the progenies derived from S. spontaneum Yunnan2017-41 (2n = 8x = 64) and Yunnan8 (2n = 10x = 80) carried a karyotype of 2n = 9x = 72, with n + n inheritance mode. These findings highlight the existence of two chromosome inheritance modes, 2n + n and n + n, in the context of the intraspecific hybridization of S. spontaneum. Additionally, hybridization between different ploidy S. spontaneum was also accompanied by chromosomal translocations (A1, A2, A18, A18) and loss (A2, A33-1, AA-4, and C2) that further resulted in the complexity of the S. spontaneum genome. Together, these findings highlight diverse chromosome inheritance in S. spontaneum hybridization, and provide a theoretical foundation for the further utilization of S. spontaneum germplasm in sugarcane breeding.

The highly allo-autopolyploid modern sugarcane genome and very recent allopolyploidization in Saccharum

Modern sugarcane, a highly allo-autopolyploid organism, has a very complex genome. In the present study, the karyotype and genome architecture of modern sugarcane were investigated, resulting in a genome assembly of 97 chromosomes (8.84 Gb). The allopolyploid genome was divided into subgenomes from Saccharum officinarum (Soh) and S. spontaneum (Ssh), with Soh dominance in the Saccharum hybrid (S. hybrid). Genome shock affected transcriptome dynamics during allopolyploidization. Analysis of an inbreeding population with 192 individuals revealed the underlying genetic basis of transgressive segregation. Population genomics of 310 Saccharum accessions clarified the breeding history of modern sugarcane. Using the haplotype-resolved S. hybrid genome as a reference, genome-wide association studies identified a potential candidate gene for sugar content from S. spontaneum. These findings illuminate the complex genome evolution of allopolyploids, offering opportunities for genomic enhancements and innovative breeding strategies for sugarcane.

ShF5H1 overexpression increases syringyl lignin and improves saccharification in sugarcane leaves

The agricultural sugarcane residues, bagasse and straws, can be used for second-generation ethanol (2GE) production by the cellulose conversion into glucose (saccharification). However, the lignin content negatively impacts the saccharification process. This polymer is mainly composed of guaiacyl (G), hydroxyphenyl (H), and syringyl (S) units, the latter formed in the ferulate 5-hydroxylase (F5H) branch of the lignin biosynthesis pathway. We have generated transgenic lines overexpressing ShF5H1 under the control of the C4H (cinnamate 4-hydroxylase) rice promoter, which led to a significant increase of up to 160% in the S/G ratio and 63% in the saccharification efficiency in leaves. Nevertheless, the content of lignin was unchanged in this organ. In culms, neither the S/G ratio nor sucrose accumulation was altered, suggesting that ShF5H1 overexpression would not affect first-generation ethanol production. Interestingly, the bagasse showed a significantly higher fiber content. Our results indicate that the tissue-specific manipulation of the biosynthetic branch leading to S unit formation is industrially advantageous and has established a foundation for further studies aiming at refining lignin modifications. Thus, the ShF5H1 overexpression in sugarcane emerges as an efficient strategy to improve 2GE production from straw.

Identification and Characterization of Key Genes for Nitrogen Utilization from Saccharum spontaneum Sub-Genome in Modern Sugarcane Cultivar

Sugarcane (Saccharum spp.) is globally considered an important crop for sugar and biofuel production. During sugarcane production, the heavy reliance on chemical nitrogen fertilizer has resulted in low nitrogen use efficiency (NUE) and high loss. Up until now, there has been extensive research on the transcriptomic dynamics during sugarcane response to low nitrogen (LN) stress. However, the specific contribution of S. spontaneum to the NUE of modern sugarcane remains unclear. In the present study, the comparative transcriptome analysis of two contrasting sugarcane cultivars in response to nitrogen deficiency was performed via the combination of genomes of S. spontaneum and S. officinarum. Sub-genome analysis indicated that S. spontaneum supports the high NUE of modern sugarcane by providing genes related to nitrogen and carbohydrate metabolism, photosynthesis, and amino acid metabolism. Additionally, the key genes involved in nitrogen metabolism from the S. spontaneum were successfully identified through weighted gene co-expression network analyses (WGCNA), and a high-affinity nitrate transporter named ScNRT2.3 was subsequently cloned. Heterogeneous expression of the ScNRT2.3, a cell membrane-localized protein, could enhance the growth of Arabidopsis under low nitrate conditions. Furthermore, a conserved protein module known as NAR2.1/NRT2.3 was shown to regulate the response to LN stress in sugarcane roots through molecular interaction. This work helps to clarify the contribution of S. spontaneum to the NUE in modern sugarcane, and the function of the ScNRT2.3 in sugarcane.

Comparison of genotyping assays for detection of targeted CRISPR/Cas mutagenesis in highly polyploid sugarcane

Sugarcane (Saccharum spp.) is an important biofuel feedstock and a leading source of global table sugar. Saccharum hybrid cultivars are highly polyploid (2n = 100-130), containing large numbers of functionally redundant hom(e)ologs in their genomes. Genome editing with sequence-specific nucleases holds tremendous promise for sugarcane breeding. However, identification of plants with the desired level of co-editing within a pool of primary transformants can be difficult. While DNA sequencing provides direct evidence of targeted mutagenesis, it is cost-prohibitive as a primary screening method in sugarcane and most other methods of identifying mutant lines have not been optimized for use in highly polyploid species. In this study, non-sequencing methods of mutant screening, including capillary electrophoresis (CE), Cas9 RNP assay, and high-resolution melt analysis (HRMA), were compared to assess their potential for CRISPR/Cas9-mediated mutant screening in sugarcane. These assays were used to analyze sugarcane lines containing mutations at one or more of six sgRNA target sites. All three methods distinguished edited lines from wild type, with co-mutation frequencies ranging from 2% to 100%. Cas9 RNP assays were able to identify mutant sugarcane lines with as low as 3.2% co-mutation frequency, and samples could be scored based on undigested band intensity. CE was highlighted as the most comprehensive assay, delivering precise information on both mutagenesis frequency and indel size to a 1 bp resolution across all six targets. This represents an economical and comprehensive alternative to sequencing-based genotyping methods which could be applied in other polyploid species.

Comparative karyotype analysis provides cytogenetic evidence for the origin of sweetpotato

The origin of hexaploid sweetpotato [Ipomoea batatas (L.) Lam.] remains controversial. Comparative karyotype analysis is particularly useful in determining species relationships and the origin of polyploid species. In previous study, we developed a set of oligo probes and identified all chromosomes of Ipomoea nil, a model diploid Ipomoea species. Here, we found that this set of oligo probes could be used to identify all chromosomes of sweetpotato and its wild relatives with different ploidy. Karyotypes based on individually identified chromosomes were established and the number and position of 5S and 35S rDNA loci were determined for these Ipomoea species. Comparison of their karyotypes revealed distinct variations in the karyotypic parameters. Karyological relationships among these species were revealed by principal coordinate analysis (PCoA) based on six quantitative parameters (x, 2n, TCL, MCA, CVCL and CVCI). These results show that I. trifida is the most closely related diploid species to sweetpotato, and other diploid species could be excluded from consideration as its possible diploid ancestor. In addition, our study also provides cytogenetic evidence for the segmental allopolyploid hypothesis of sweetpotato origin.

Developing Adventitious Root Meristems Induced by Layering for Plant Chromosome Preparation

Chromosome numbers and morphology are important characteristics of a species and its evolution. Root tips are the most commonly used tissue as a source of actively dividing cells for chromosome visualization in plants. Previously, rapidly growing root tips were collected from germinating kernels or from seedlings growing in pots or fields. However, the use of adventitious roots (ARs) derived from aerial tissue as meristems for chromosome visualization has always been overlooked. Here, we successfully induced ARs in 12 materials that were investigated, with the exception of Sorghum nitidum. Using ARs meristem we obtained high-quality chromosome spreads for Morus alba, Broussonetia papyrifera, Lolium multiflorum, Sorghum sudanense, S. propinquum, S. bicolor × S. sudanense, Zea mays, Z. mexicana, Glycine max, Medicago sativa, and Brassica napus. The results reported here demonstrate that layering is an alternative and effective method for producing meristematic cells for high-quality chromosome preparation in plant species producing ARs. For species that produce ARs by layering, this protocol is particularly valuable for the development of cost-effective and high-throughput non-invasive cytogenetic studies.

The extent of multiallelic, co-editing of LIGULELESS1 in highly polyploid sugarcane tunes leaf inclination angle and enables selection of the ideotype for biomass yield

Sugarcane (Saccharum spp. hybrid) is a prime feedstock for commercial production of biofuel and table sugar. Optimizing canopy architecture for improved light capture has great potential for elevating biomass yield. LIGULELESS1 (LG1) is involved in leaf ligule and auricle development in grasses. Here, we report CRISPR/Cas9-mediated co-mutagenesis of up to 40 copies/alleles of the putative LG1 in highly polyploid sugarcane (2n = 100-120, x = 10-12). Next generation sequencing revealed co-editing frequencies of 7.4%-100% of the LG1 reads in 16 of the 78 transgenic lines. LG1 mutations resulted in a tuneable leaf angle phenotype that became more upright as co-editing frequency increased. Three lines with loss of function frequencies of ~12%, ~53% and ~95% of lg1 were selected following a randomized greenhouse trial and grown in replicated, multi-row field plots. The co-edited LG1 mutations were stably maintained in vegetative progenies and the extent of co-editing remained constant in field tested lines L26 and L35. Next generation sequencing confirmed the absence of potential off targets. The leaf inclination angle corresponded to light transmission into the canopy and tiller number. Line L35 displaying loss of function in ~12% of the lg1 NGS reads exhibited an 18% increase in dry biomass yield supported by a 56% decrease in leaf inclination angle, a 31% increase in tiller number, and a 25% increase in internode number. The scalable co-editing of LG1 in highly polyploid sugarcane allows fine-tuning of leaf inclination angle, enabling the selection of the ideotype for biomass yield.

In vivo induction of sugarcane (Saccharum spp.) haploids by genome editing

Mutation of the MATRILINEAL gene in sugarcane induces in vivo maternal haploids and suggests prospects for sugarcane breeding.

Sequencing vs. amplification for the estimation of allele dosages in sugarcane (Saccharum spp.)

Premise

Detecting single-nucleotide polymorphisms (SNPs) in a cost-effective way is fundamental in any plant breeding pipeline. Here, we compare three genotyping techniques for their ability to reproduce the allele dosage of SNPs of interest in sugarcane (Saccharum spp.).

Methods

To identify a reproducible technique to estimate allele dosage for the validation of SNP markers, the correlation between Flex-Seq, kompetitive allele-specific PCR (KASP), and genotyping-by-sequencing and restriction site-associated DNA sequencing (GBS+RADseq) was determined for a set of 76 SNPs. To find alternative methodologies for allele dosage estimation, the KASP and Flex-Seq techniques were compared for the same set of SNPs. For the three techniques, a population of 53 genotypes from the diverse sugarcane panel of the Centro de Investigación de la Caña de Azúcar (Cenicaña), Colombia, was selected.

Results

The average Pearson correlation coefficients between GBS+RADseq and Flex-Seq, GBS+RADseq and KASP, and Flex-Seq and KASP were 0.62 ± 0.27, 0.38 ± 0.27, and 0.38 ± 0.30, respectively.

Discussion

Flex-Seq reproduced the allele dosages determined using GBS+RADseq with good levels of precision because of its depth of sequencing and ability to target specific positions in the genome. Additionally, Flex-Seq outperformed KASP by allowing the conversion of a higher number of SNPs and a more accurate estimation of the allele dosage. Flex-Seq has therefore become the genotyping methodology of choice for marker validation at Cenicaña.

Molecular evaluation and phenotypic screening of brown and orange rust in Saccharum germplasm

Brazil is the largest global producer of sugarcane and plays a significant role-supplier of sugar and bioethanol. However, diseases such as brown and orange rust cause substantial yield reductions and economic losses, due decrease photosynthesis and biomass in susceptible cultivars. Molecular markers associated with resistance genes, such as Bru1 (brown rust) and G1 (orange rust), could aid in predicting resistant genotypes. In this study, we sought to associate the phenotypic response of 300 sugarcane accessions with the genotypic response of Bru1 and G1 markers. The field trials were conducted in a randomized block design, and five six-month-old plants per plot were evaluated under natural disease conditions. Genotypic information about the presence or absence of Bru1 (haplotype 1) and G1 gene was obtained after extraction of genomic DNA and conventional PCR. Of the total accessions evaluated, 60.3% (181) showed resistance to brown rust in the field, and of these, 70.7% (128) had the Bru1 gene present. Considering the field-resistant accessions obtained from Brazilian breeding programs (116), the Bru1 was present in 77,6% of these accessions. While alternative resistance sources may exist, Bru1 likely confers enduring genetic resistance in current Brazilian cultivars. Regarding the phenotypic reaction to orange rust, the majority of accessions, 96.3% (288), were field resistant, and of these, 52.7% (152) carried the G1 marker. Although less efficient for predicting resistance when compared to Bru1, the G1 marker could be part of a quantitative approach when new orange rust resistance genes are described. Therefore, these findings showed the importance of Bru1 molecular markers for the early selection of resistant genotypes to brown rust by genetic breeding programs.

Genetic Engineering for Enhancing Sugarcane Tolerance to Biotic and Abiotic Stresses

Sugarcane, a vital cash crop, contributes significantly to the world's sugar supply and raw materials for biofuel production, playing a significant role in the global sugar industry. However, sustainable productivity is severely hampered by biotic and abiotic stressors. Genetic engineering has been used to transfer useful genes into sugarcane plants to improve desirable traits and has emerged as a basic and applied research method to maintain growth and productivity under different adverse environmental conditions. However, the use of transgenic approaches remains contentious and requires rigorous experimental methods to address biosafety challenges. Clustered regularly interspaced short palindromic repeat (CRISPR) mediated genome editing technology is growing rapidly and may revolutionize sugarcane production. This review aims to explore innovative genetic engineering techniques and their successful application in developing sugarcane cultivars with enhanced resistance to biotic and abiotic stresses to produce superior sugarcane cultivars.

Pathogenesis-Related 1 (PR1) Protein Family Genes Involved in Sugarcane Responses to Ustilago scitaminea Stress

Plant resistance against biotic stressors is significantly influenced by pathogenesis-related 1 (PR1) proteins. This study examines the systematic identification and characterization of PR1 family genes in sugarcane (Saccharum spontaneum Np-X) and the transcript expression of selected genes in two sugarcane cultivars (ROC22 and Zhongtang3) in response to Ustilago scitaminea pathogen infection. A total of 18 SsnpPR1 genes were identified at the whole-genome level and further categorized into four groups. Notably, tandem and segmental duplication occurrences were detected in one and five SsnpPR1 gene pairs, respectively. The SsnpPR1 genes exhibited diverse physio-chemical attributes and variations in introns/exons and conserved motifs. Notably, four SsnpPR1 (SsnpPR1.02/05/09/19) proteins displayed a strong protein-protein interaction network. The transcript expression of three SsnpPR1 (SsnpPR1.04/06/09) genes was upregulated by 1.2-2.6 folds in the resistant cultivar (Zhongtang3) but downregulated in the susceptible cultivar (ROC22) across different time points as compared to the control in response to pathogen infection. Additionally, SsnpPR1.11 was specifically upregulated by 1.2-3.5 folds at 24-72 h post inoculation (hpi) in ROC22, suggesting that this gene may play an important negative regulatory role in defense responses to pathogen infection. The genetic improvement of sugarcane can be facilitated by our results, which also establish the basis for additional functional characterization of SsnpPR1 genes in response to pathogenic stress.

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.

The complex polyploid genome architecture of sugarcane

Sugarcane, the world's most harvested crop by tonnage, has shaped global history, trade and geopolitics, and is currently responsible for 80% of sugar production worldwide1. While traditional sugarcane breeding methods have effectively generated cultivars adapted to new environments and pathogens, sugar yield improvements have recently plateaued2. The cessation of yield gains may be due to limited genetic diversity within breeding populations, long breeding cycles and the complexity of its genome, the latter preventing breeders from taking advantage of the recent explosion of whole-genome sequencing that has benefited many other crops. Thus, modern sugarcane hybrids are the last remaining major crop without a reference-quality genome. Here we take a major step towards advancing sugarcane biotechnology by generating a polyploid reference genome for R570, a typical modern cultivar derived from interspecific hybridization between the domesticated species (Saccharum officinarum) and the wild species (Saccharum spontaneum). In contrast to the existing single haplotype ('monoploid') representation of R570, our 8.7 billion base assembly contains a complete representation of unique DNA sequences across the approximately 12 chromosome copies in this polyploid genome. Using this highly contiguous genome assembly, we filled a previously unsized gap within an R570 physical genetic map to describe the likely causal genes underlying the single-copy Bru1 brown rust resistance locus. This polyploid genome assembly with fine-grain descriptions of genome architecture and molecular targets for biotechnology will help accelerate molecular and transgenic breeding and adaptation of sugarcane to future environmental conditions.

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.

Lignin reduction in sugarcane by performing CRISPR/Cas9 site-direct mutation of SoLIM transcription factor

Genetic engineering of plant cell walls is limited for reducing lignocellulose recalcitrance, so mild and/or green-like pretreatment is still required for sequential enzymatic saccharification. Here, we report a method to reduce lignin content in sugarcane stalks using the CRISPR/Cas 9 technique. Three target sequences of SoLIM were designed and fused to pRGEB32. The cassette constructs were introduced into sugarcane calli cv. KK3 through Agrobacterium-mediated transformation. We produced one base substitution and one insertion line for the 1st target site; two insertions, one deletion, and one base substitution for the 2nd target site; and one base substitution and insertion for the 3rd target site. qRT-PCR analysis of SoLIM, SoPAL, SoC4H, and SoCAD showeded that downregulation of SoLIM by single nucleotide insertions or deletions reduced the expression of SoPAL, SoC4H, and SoCAD. Consequently, the edited lines contained 9.74 to 51.46% less lignin content compared to that in the wild-type plants. The syringyl/guaiacyl (S/G) ratio of the edited lines ranged between 0.23 and 0.49, while the wild-type was 0.22. The histochemical evaluation and scanning electron microscopy of the cell walls supported this observation. A low lignin content sugarcane will provide a better feedstock for second-generation bioethanol production.

Translocations and inversions: major chromosomal rearrangements during Vigna (Leguminosae) evolution

KEY MESSAGE

Inversions and translocations are the major chromosomal rearrangements involved in Vigna subgenera evolution, being Vigna vexillata the most divergent species. Centromeric repositioning seems to be frequent within the genus. Oligonucleotide-based fluorescence in situ hybridization (Oligo-FISH) provides a powerful chromosome identification system for inferring plant chromosomal evolution. Aiming to understand macrosynteny, chromosomal diversity, and the evolution of bean species from five Vigna subgenera, we constructed cytogenetic maps for eight taxa using oligo-FISH-based chromosome identification. We used oligopainting probes from chromosomes 2 and 3 of Phaseolus vulgaris L. and two barcode probes designed from V. unguiculata (L.) Walp. genome. Additionally, we analyzed genomic blocks among the Ancestral Phaseoleae Karyotype (APK), two V. unguiculata subspecies (V. subg. Vigna), and V. angularis (Willd.) Ohwi & Ohashi (V. subg. Ceratotropis). We observed macrosynteny for chromosomes 2, 3, 4, 6, 7, 8, 9, and 10 in all investigated taxa except for V. vexillata (L.) A. Rich (V. subg. Plectrotropis), in which only chromosomes 4, 7, and 9 were unambiguously identified. Collinearity breaks involved with chromosomes 2 and 3 were revealed. We identified minor differences in the painting pattern among the subgenera, in addition to multiple intra- and interblock inversions and intrachromosomal translocations. Other rearrangements included a pericentric inversion in chromosome 4 (V. subg. Vigna), a reciprocal translocation between chromosomes 1 and 5 (V. subg. Ceratotropis), a potential deletion in chromosome 11 of V. radiata (L.) Wilczek, as well as multiple intrablock inversions and centromere repositioning via genomic blocks. Our study allowed the visualization of karyotypic patterns in each subgenus, revealing important information for understanding intrageneric karyotypic evolution, and suggesting V. vexillata as the most karyotypically divergent species.

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.

Developing new sugarcane varieties suitable for mechanized production in China: principles, strategies and prospects

The sugar industry, which relates to people's livelihood, is strategic and fundamental in the development of agricultural economy. In China, sugar derived from sugarcane accounts for approximately 85% of total sugar production. Mechanization is the "flower" of sugarcane industry. As the saying goes "when there are blooming flowers, there will be sweet honey." However, due to limitations in land resources, technology, equipment, organization, and management, mechanization throughout the sugarcane production process has not yet brought about the economic benefits that a mechanized system should provide and has not reached an ideal yield through the integration of agricultural machinery and agronomic practice. This paper briefly describes how to initiate the mechanization of Chinese sugarcane production to promote the sound, healthy, and rapid development of the sugarcane industry, and how to ultimately achieve the transformation of sugarcane breeding in China and the modernization of the sugarcane industry from three perspectives, namely, requirements of mechanized production for sugarcane varieties, breeding strategies for selecting new sugarcane varieties suitable for mechanized production, and screening for sugarcane varieties that are suitable for mechanization and diversification in variety distribution or arrangement in China. We also highlight the current challenges surrounding this topic and look forward to its bright prospects.

Optimising clonal performance in sugarcane: leveraging non-additive effects via mate-allocation strategies

Mate-allocation strategies in breeding programs can improve progeny performance by harnessing non-additive genetic effects. These approaches prioritise predicted progeny merit over parental breeding value, making them particularly appealing for clonally propagated crops such as sugarcane. We conducted a comparative analysis of mate-allocation strategies, exploring utilising non-additive and heterozygosity effects to maximise clonal performance with schemes that solely consider additive effects to optimise breeding value. Using phenotypic and genotypic data from a population of 2,909 clones evaluated in final assessment trials of Australian sugarcane breeding programs, we focused on three important traits: tonnes of cane per hectare (TCH), commercial cane sugar (CCS), and Fibre. By simulating families from all possible crosses (1,225) with 50 progenies each, we predicted the breeding and clonal values of progeny using two models: GBLUP (considering additive effects only) and extended-GBLUP (incorporating additive, non-additive, and heterozygosity effects). Integer linear programming was used to identify the optimal mate-allocation among selected parents. Compared to breeding value-based approaches, mate-allocation strategies based on clonal performance yielded substantial improvements, with predicted progeny values increasing by 57% for TCH, 12% for CCS, and 16% for fibre. Our simulation study highlights the effectiveness of mate-allocation approaches that exploit non-additive and heterozygosity effects, resulting in superior clonal performance. However, there was a notable decline in additive gain, particularly for TCH, likely due to significant epistatic effects. When selecting crosses based on clonal performance for TCH, the inbreeding coefficient of progeny was significantly lower compared to random mating, underscoring the advantages of leveraging non-additive and heterozygosity effects in mitigating inbreeding depression. Thus, mate-allocation strategies are recommended in clonally propagated crops to enhance clonal performance and reduce the negative impacts of inbreeding.

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.

High rates of structural rearrangements have shaped the chromosome evolution in dysploid Phaseolus beans

KEY MESSAGE

Karyotypes evolve through numerical and structural chromosome rearrangements. We show that Phaseolus leptostachyus, a wild bean, underwent a rapid genome reshuffling associated with the reduction from 11 to 10 chromosome pairs, but without whole genome duplication, the highest chromosome evolution rate known for plants. Plant karyotypes evolve through structural rearrangements often associated with polyploidy or dysploidy. The genus Phaseolus comprises ~ 90 species, five of them domesticated due to their nutritional relevance. Most of the species have 2n = 22 karyotypes and are highly syntenic, except for three dysploid karyotypes of species from the Leptostachyus group (2n = 20) that have accumulated several rearrangements. Here, we investigated the degrees of structural rearrangements among Leptostachyus and other Phaseolus groups by estimating their chromosomal evolution rates (CER). For this, we combined our oligo-FISH barcode system for beans and chromosome-specific painting probes for chromosomes 2 and 3, with rDNA and a centromeric probe to establish chromosome orthologies and identify structural rearrangements across nine Phaseolus species. We also integrated the detected rearrangements with a phylogenomic approach to estimate the CERs for each Phaseolus lineage. Our data allowed us to identify translocations, inversions, duplications and deletions, mostly in species belonging to the Leptostachyus group. Phaseolus leptostachyus showed the highest CER (12.31 rearrangements/My), a tenfold increase in contrast to the 2n = 22 species analysed. This is the highest rate known yet for plants, making it a model species for investigating the mechanisms behind rapid genome reshuffling in early species diversification.

Identification of ARF transcription factor gene family and its defense responses to bacterial infection and salicylic acid treatment in sugarcane

Auxin response factor (ARF) is a critical regulator in the auxin signaling pathway, involved in a variety of plant biological processes. Here, gene members of 24 SpapARFs and 39 SpnpARFs were identified in two genomes of Saccharum spontaneum clones AP85-441 and Np-X, respectively. Phylogenetic analysis showed that all ARF genes were clustered into four clades, which is identical to those ARF genes in maize (Zea mays) and sorghum (Sorghum bicolor). The gene structure and domain composition of this ARF family are conserved to a large degree across plant species. The SpapARF and SpnpARF genes were unevenly distributed on chromosomes 1-8 and 1-10 in the two genomes of AP85-441 and Np-X, respectively. Segmental duplication events may also contribute to this gene family expansion in S. spontaneum. The post-transcriptional regulation of ARF genes likely involves sugarcane against various stressors through a miRNA-medicated pathway. Expression levels of six representative ShARF genes were analyzed by qRT-PCR assays on two sugarcane cultivars [LCP85-384 (resistant to leaf scald) and ROC20 (susceptible to leaf scald)] triggered by Acidovorax avenae subsp. avenae (Aaa) and Xanthomonas albilineans (Xa) infections and salicylic acid (SA) treatment. ShARF04 functioned as a positive regulator under Xa and Aaa stress, whereas it was a negative regulator under SA treatment. ShARF07/17 genes played positive roles against both pathogenic bacteria and SA stresses. Additionally, ShARF22 was negatively modulated by Xa and Aaa stimuli in both cultivars, particularly LCP85-384. These findings imply that sugarcane ARFs exhibit functional redundancy and divergence against stressful conditions. This work lays the foundation for further research on ARF gene functions in sugarcane against diverse environmental stressors.

Oligonucleotide Fluorescence In Situ Hybridization: An Efficient Chromosome Painting Method in Plants

Fluorescence in situ hybridization (FISH) is an indispensable technique for studying chromosomes in plants. However, traditional FISH methods, such as BAC, rDNA, tandem repeats, and distributed repetitive sequence probe-based FISH, have certain limitations, including difficulties in probe synthesis, low sensitivity, cross-hybridization, and limited resolution. In contrast, oligo-based FISH represents a more efficient method for chromosomal studies in plants. Oligo probes are computationally designed and synthesized for any plant species with a sequenced genome and are suitable for single and repetitive DNA sequences, entire chromosomes, or chromosomal segments. Furthermore, oligo probes used in the FISH experiment provide high specificity, resolution, and multiplexing. Moreover, oligo probes made from one species are applicable for studying other genetically and taxonomically related species whose genome has not been sequenced yet, facilitating molecular cytogenetic studies of non-model plants. However, there are some limitations of oligo probes that should be considered, such as requiring prior knowledge of the probe design process and FISH signal issues with shorter probes of background noises during oligo-FISH experiments. This review comprehensively discusses de novo oligo probe synthesis with more focus on single-copy DNA sequences, preparation, improvement, and factors that affect oligo-FISH efficiency. Furthermore, this review highlights recent applications of oligo-FISH in a wide range of plant chromosomal studies.

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

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

Reconstruction of karyotypic evolution in Saccharum spontaneum species by comparative oligo-FISH mapping

BACKGROUND

Karyotype dynamics driven by chromosomal rearrangements has long been considered as a fundamental question in the evolutionary genetics. Saccharum spontaneum, the most primitive and complex species in the genus Saccharum, has reportedly undergone at least two major chromosomal rearrangements, however, its karyotypic evolution remains unclear.

RESULTS

In this study, four representative accessions, i.e., hypothetical diploid sugarcane ancestor (sorghum, x = 10), Sa. spontaneum Np-X (x = 10, tetraploid), 2012-46 (x = 9, hexaploid) and AP85-441 (x = 8, tetraploid), were selected for karyotype evolution studies. A set of oligonucleotide (oligo)-based barcode probes was developed based on the sorghum genome, which allowed universal identification of all chromosomes from sorghum and Sa. spontaneum. By comparative FISH assays, we reconstructed the karyotype evolutionary history and discovered that although chromosomal rearrangements resulted in greater variation in relative lengths of some chromosomes, all chromosomes maintained a conserved metacentric structure. Additionally, we found that the barcode oligo probe was not applicable for chromosome identification in both Sa. robustum and Sa. officinarum species, suggesting that sorghum is more distantly related to Sa. robustum and Sa. officinarum compared with Sa. spontaneum species.

CONCLUSIONS

Our study demonstrated that the barcode oligo-FISH is an efficient tool for chromosome identification and karyotyping research, and expanded our understanding of the karyotypic and chromosomal evolution in the genus Saccharum.

Screening of Candidate Genes Associated with Brown Stripe Resistance in Sugarcane via BSR-seq Analysis

Sugarcane brown stripe (SBS), caused by the fungal pathogen Helminthosporium stenospilum, is one of the most serious threats to sugarcane production. However, its outbreaks and epidemics require suitable climatic conditions, resulting in the inefficient improvement of the SBS resistance by phenotype selection. The sugarcane F1 population of SBS-resistant YT93-159 × SBS-susceptible ROC22 was used for constructing the bulks. Bulked segregant RNA-seq (BSR-seq) was then performed on the parents YT93-159 (T01) and ROC22 (T02), and the opposite bulks of 30 SBS-susceptible individuals mixed bulk (T03) and 30 SBS-resistant individuals mixed bulk (T04) collected from 287 F1 individuals. A total of 170.00 Gb of clean data containing 297,921 SNPs and 70,426 genes were obtained. Differentially expressed genes (DEGs) analysis suggested that 7787 and 5911 DEGs were identified in the parents (T01 vs. T02) and two mixed bulks (T03 vs. T04), respectively. In addition, 25,363 high-quality and credible SNPs were obtained using the genome analysis toolkit GATK for SNP calling. Subsequently, six candidate regions with a total length of 8.72 Mb, which were located in the chromosomes 4B and 7C of sugarcane wild species Saccharum spontaneum, were identified, and 279 genes associated with SBS-resistance were annotated by ED algorithm and ΔSNP-index. Furthermore, the expression profiles of candidate genes were verified by quantitative real-time PCR (qRT-PCR) analysis, and the results showed that eight genes (LRR-RLK, DHAR1, WRKY7, RLK1, BLH4, AK3, CRK34, and NDA2) and seven genes (WRKY31, CIPK2, CKA1, CDPK6, PFK4, CBL2, and PR2) of the 20 tested genes were significantly up-regulated in YT93-159 and ROC22, respectively. Finally, a potential molecular mechanism of sugarcane response to H. stenospilum infection is illustrate that the activations of ROS signaling, MAPK cascade signaling, Ca2+ signaling, ABA signaling, and the ASA-GSH cycle jointly promote the SBS resistance in sugarcane. This study provides abundant gene resources for the SBS resistance breeding in sugarcane.

Combining genomic selection with genome-wide association analysis identified a large-effect QTL and improved selection for red rot resistance in sugarcane

Red rot caused by the fungus Colletotrichum falcatum is the main disease limiting sugarcane productivity in several countries including the major producer India. The genetic basis for red rot resistance is unclear. We studied a panel of 305 sugarcane clones from the Australian breeding program for disease response phenotype and genotype using an Affymetrix® Axiom® array, to better understand the genetic basis of red rot resistance. SNP markers highly significantly associated with red rot response (≤ 10-8) were identified. Markers with largest effect were located in a single 14.6 Mb genomic region of sorghum (the closest diploid relative of sugarcane with a sequenced genome) suggesting the presence of a major-effect QTL. By genomic selection, the estimated selection accuracy was ~0.42 for red rot resistance. This was increased to ~0.5 with the addition of 29 highly significant SNPs as fixed effects. Analysis of genes nearby the markers linked to the QTL revealed many biotic stress responsive genes within this QTL, with the most significant SNP co-locating with a cluster of four chitinase A genes. The SNP markers identified here could be used to predict red rot resistance with high accuracy at any stage in the sugarcane breeding program.

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

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

Recent Advances in Sugarcane Genomics, Physiology, and Phenomics for Superior Agronomic Traits

Advances in sugarcane breeding have contributed significantly to improvements in agronomic traits and crop yield. However, the growing global demand for sugar and biofuel in the context of climate change requires further improvements in cane and sugar yields. Attempts to achieve the desired rates of genetic gain in sugarcane by conventional breeding means are difficult as many agronomic traits are genetically complex and polygenic, with each gene exerting small effects. Unlike those of many other crops, the sugarcane genome is highly heterozygous due to its autopolyploid nature, which further hinders the development of a comprehensive genetic map. Despite these limitations, many superior agronomic traits/genes for higher cane yield, sugar production, and disease/pest resistance have been identified through the mapping of quantitative trait loci, genome-wide association studies, and transcriptome approaches. Improvements in traits controlled by one or two loci are relatively easy to achieve; however, this is not the case for traits governed by many genes. Many desirable phenotypic traits are controlled by quantitative trait nucleotides (QTNs) with small and variable effects. Assembling these desired QTNs by conventional breeding methods is time consuming and inefficient due to genetic drift. However, recent developments in genomics selection (GS) have allowed sugarcane researchers to select and accumulate desirable alleles imparting superior traits as GS is based on genomic estimated breeding values, which substantially increases the selection efficiency and genetic gain in sugarcane breeding programs. Next-generation sequencing techniques coupled with genome-editing technologies have provided new vistas in harnessing the sugarcane genome to look for desirable agronomic traits such as erect canopy, leaf angle, prolonged greening, high biomass, deep root system, and the non-flowering nature of the crop. Many desirable cane-yielding traits, such as single cane weight, numbers of tillers, numbers of millable canes, as well as cane quality traits, such as sucrose and sugar yield, have been explored using these recent biotechnological tools. This review will focus on the recent advances in sugarcane genomics related to genetic gain and the identification of favorable alleles for superior agronomic traits for further utilization in sugarcane breeding programs.

Genomic insights into the recent chromosome reduction of autopolyploid sugarcane Saccharum spontaneum

Saccharum spontaneum is a founding Saccharum species and exhibits wide variation in ploidy levels. We have assembled a high-quality autopolyploid genome of S. spontaneum Np-X (2n = 4x = 40) into 40 pseudochromosomes across 10 homologous groups, that better elucidates recent chromosome reduction and polyploidization that occurred circa 1.5 million years ago (Mya). One paleo-duplicated chromosomal pair in Saccharum, NpChr5 and NpChr8, underwent fission followed by fusion accompanied by centromeric split around 0.80 Mya. We inferred that Np-X, with x = 10, most likely represents the ancestral karyotype, from which x = 9 and x = 8 evolved. Resequencing of 102 S. spontaneum accessions revealed that S. spontaneum originated in northern India from an x = 10 ancestor, which then radiated into four major groups across the Indian subcontinent, China, and Southeast Asia. Our study suggests new directions for accelerating sugarcane improvement and expands our knowledge of the evolution of autopolyploids.

Rapid and visual monitoring of alien sequences using crop wild relatives specific oligo-painting: The case of cucumber chromosome engineering

Wild species related to domesticated crops (crop wild relatives, or CWRs) represent a high level of genetic diversity that provides a practical gene pool for crop pre-breeding employed to address climate change and food demand challenges globally. Nevertheless, rapid identifying and visual tracking of alien chromosomes and sequences derived from CWRs have been a technical challenge for crop chromosome engineering. Here, a species-specific oligonucleotide (oligo) pool was developed by using the reference genome of Cucumis hystrix (HH, 2n = 2x = 24), a wild species carrying many favorable traits and interspecific compatibility with cultivated cucumber (C. sativus, CC, 2n = 2x = 14). These synthetic double-stranded oligo probes were applied to validate the assembly and characterize the chromosome architectures of C. hystrix, as well as to rapidly identify C. hystrix-chromosomes in diverse C. sativus-hystrix chromosome-engineered germplasms, including interspecific hybrid F1 (HC), synthetic allopolyploids (HHCC, CHC, and HCH) and alien additional lines (CC-H). Moreover, a ∼2Mb of C. hystrix-specific sequences, introduced into cultivated cucumber, were visualized by CWR-specific oligo-painting. These results demonstrate that the CWR-specific oligo-painting technique holds broad applicability for chromosome engineering of numerous crops, as it allows rapid identification of alien chromosomes, reliable detection of homoeologous recombination, and visual tracking of the introgression process. It is promising to achieve directed and high-precision crop pre-breeding combined with other breeding techniques, such as CRISPR/Cas9-mediated chromosome engineering.

Allele expression biases in mixed-ploid sugarcane accessions

Allele-specific expression (ASE) represents differences in the magnitude of expression between alleles of the same gene. This is not straightforward for polyploids, especially autopolyploids, as knowledge about the dose of each allele is required for accurate estimation of ASE. This is the case for the genomically complex Saccharum species, characterized by high levels of ploidy and aneuploidy. We used a Beta-Binomial model to test for allelic imbalance in Saccharum, with adaptations for mixed-ploid organisms. The hierarchical Beta-Binomial model was used to test if allele expression followed the expectation based on genomic allele dosage. The highest frequencies of ASE occurred in sugarcane hybrids, suggesting a possible influence of interspecific hybridization in these genotypes. For all accessions, genes showing ASE (ASEGs) were less frequent than those with balanced allelic expression. These genes were related to a broad range of processes, mostly associated with general metabolism, organelles, responses to stress and responses to stimuli. In addition, the frequency of ASEGs in high-level functional terms was similar among the genotypes, with a few genes associated with more specific biological processes. We hypothesize that ASE in Saccharum is largely a genotype-specific phenomenon, as a large number of ASEGs were exclusive to individual accessions.

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

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

Chromosome-specific painting unveils chromosomal fusions and distinct allopolyploid species in the Saccharum complex

Karyotypes provide key cytogenetic information on the phylogenetic relationships and evolutionary origins in related eukaryotic species. Despite our knowledge of the chromosome numbers of sugarcane and its wild relatives, the chromosome composition and evolution among the species in the Saccharum complex have been elusive owing to the complex polyploidy and the large numbers of chromosomes of these species. Oligonucleotide-based chromosome painting has become a powerful tool of cytogenetic studies especially for plant species with large numbers of chromosomes. We developed oligo-based chromosome painting probes for all 10 chromosomes in Saccharum officinarum (2n = 8x = 80). The 10 painting probes generated robust fluorescence in situ hybridization signals in all plant species within the Saccharum complex, including species in the genera Saccharum, Miscanthus, Narenga and Erianthus. We conducted comparative chromosome analysis using the same set of probes among species from four different genera within the Saccharum complex. Excitingly, we discovered several novel cytotypes and chromosome rearrangements in these species. We discovered that fusion from two different chromosomes is a common type of chromosome rearrangement associated with the species in the Saccharum complex. Such fusion events changed the basic chromosome number and resulted in distinct allopolyploids in the Saccharum complex.

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.

A comprehensive molecular cytogenetic analysis of the genome architecture in modern sugarcane cultivars

Modern sugarcane cultivars are derived from the hybridization of Saccharum officinarum (2n = 80) and S. spontaneum (2n = 40-128), leading to a variety of complex genomes with highly polyploid and varied chromosome structures. These complex genomes have hindered deciphering the genome structure and marker-assisted selection in sugarcane breeding. Ten cultivars were analyzed by fluorescence in situ hybridization adopting chromosome painting and S. spontaneum-specific probes. The results showed six types of chromosomes in the studied cultivars, including S. spontaneum or S. officinarum chromosomes, interspecific recombinations from homoeologous or nonhomoeologous chromosomes, and translocations of S. spontaneum or S. officinarum chromosomes. The results showed unexpectedly high proportions of interspecific recombinations in these cultivars (11.9-40.9%), which renew our knowledge that less than 13% of chromosomes result from interspecific exchanges. Also, the results showed a high frequency of translocations (an average of 2.15 translocations per chromosome) between S. officinarum chromosomes. The diverse types of chromosomes in cultivars imply that hybrid gametes of S. spontaneum and S. officinarum may form unusual chromosome pairs, including homoeologous or nonhomoeologous chromosomes either between or within S. spontaneum and S. officinarum. Moreover, we consistently observed 11 or 12 copies for the four studied chromosomes, i.e., chromosomes 1, 2, 7, and 8, suggesting steady transmission during the breeding program. By comparison, we found a relatively fewer copies of S. spontaneum chromosome 1 than those of S. spontaneum chromosomes 2, 7, and 8. These results provide deep insights into the structure of cultivars and may facilitate chromosome-assisted selection in sugarcane breeding.

Temporal Gene Expression in Apical Culms Shows Early Changes in Cell Wall Biosynthesis Genes in Sugarcane

Multiple genes in sugarcane control sucrose accumulation and the biosynthesis of cell wall components; however, it is unclear how these genes are expressed in its apical culms. To better understand this process, we sequenced mRNA from +1 stem internodes collected from four genotypes with different concentrations of soluble solids. Culms were collected at four different time points, ranging from six to 12-month-old plants. Here we show differentially expressed genes related to sucrose metabolism and cell wall biosynthesis, including genes encoding invertases, sucrose synthase and cellulose synthase. Our results showed increased expression of invertases in IN84-58, the genotype with lower sugar and higher fiber content, as well as delayed expression of secondary cell wall-related cellulose synthase for the other genotypes. Interestingly, genes involved with hormone metabolism were differentially expressed across time points in the three genotypes with higher soluble solids content. A similar result was observed for genes controlling maturation and transition to reproductive stages, possibly a result of selection against flowering in sugarcane breeding programs. These results indicate that carbon partitioning in apical culms of contrasting genotypes is mainly associated with differential cell wall biosynthesis, and may include early modifications for subsequent sucrose accumulation. Co-expression network analysis identified transcription factors related to growth and development, showing a probable time shift for carbon partitioning occurred in 10-month-old plants.

Oligo-FISH barcode in beans: a new chromosome identification system

An Oligo-FISH barcode system was developed for two model legumes, allowing the identification of all cowpea and common bean chromosomes in a single FISH experiment, and revealing new chromosome rearrangements. The FISH barcode system emerges as an effective tool to understand the chromosome evolution of economically important legumes and their related species. Current status on plant cytogenetic and cytogenomic research has allowed the selection and design of oligo-specific probes to individually identify each chromosome of the karyotype in a target species. Here, we developed the first chromosome identification system for legumes based on oligo-FISH barcode probes. We selected conserved genomic regions between Vigna unguiculata (Vu, cowpea) and Phaseolus vulgaris (Pv, common bean) (diverged ~ 9.7-15 Mya), using cowpea as a reference, to produce a unique barcode pattern for each species. We combined our oligo-FISH barcode pattern with a set of previously developed FISH probes based on BACs and ribosomal DNA sequences. In addition, we integrated our FISH maps with genome sequence data. Based on this integrated analysis, we confirmed two translocation events (involving chromosomes 1, 5, and 8; and chromosomes 2 and 3) between both species. The application of the oligo-based probes allowed us to demonstrate the participation of chromosome 5 in the translocation complex for the first time. Additionally, we detailed a pericentric inversion on chromosome 4 and identified a new paracentric inversion on chromosome 10. We also detected centromere repositioning associated with chromosomes 2, 3, 5, 7, and 9, confirming previous results for chromosomes 2 and 3. This first barcode system for legumes can be applied for karyotyping other Phaseolinae species, especially non-model, orphan crop species lacking genomic assemblies and cytogenetic maps, expanding our understanding of the chromosome evolution and genome organization of this economically important legume group.

Chromosome identification in oil palm (Elaeis guineensis) using in situ hybridization with massive pools of single copy oligonucleotides and transferability across Arecaceae species

Chromosome identification is essential for linking sequence and chromosomal maps, verifying sequence assemblies, showing structural variations and tracking inheritance or recombination of chromosomes and chromosomal segments during evolution and breeding programs. Unfortunately, identification of individual chromosomes and chromosome arms has been a major challenge for some economically important crop species with a near-continuous chromosome size range and similar morphology. Here, we developed oligonucleotide-based chromosome-specific probes that enabled us to establish a reference chromosome identification system for oil palm (Elaeis guineensis Jacq., 2n = 32). Massive oligonucleotide sequence pools were anchored to individual chromosome arms using dual and triple fluorescent in situ hybridization (EgOligoFISH). Three fluorescently tagged probe libraries were developed to contain, in total 52,506 gene-rich single-copy 47-mer oligonucleotides spanning each 0.2-0.5 Mb across strategically placed chromosome regions. They generated 19 distinct FISH signals and together with rDNA probes enabled identification of all 32 E. guineensis chromosome arms. The probes were able to identify individual homoeologous chromosome regions in the related Arecaceae palm species: American oil palm (Elaeis oleifera), date palm (Phoenix dactylifera) and coconut (Cocos nucifera) showing the comparative organization and concerted evolution of genomes in the Arecaceae. The oligonucleotide probes developed here provide a valuable approach to chromosome arm identification and allow tracking chromosome transfer in hybridization and breeding programs in oil palm, as well as comparative studies within Arecaceae.

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.

Arabidopsis-Based Dual-Layered Biological Network Analysis Elucidates Fully Modulated Pathways Related to Sugarcane Resistance on Biotrophic Pathogen Infection

We assembled a dual-layered biological network to study the roles of resistance gene analogs (RGAs) in the resistance of sugarcane to infection by the biotrophic fungus causing smut disease. Based on sugarcane-Arabidopsis orthology, the modeling used metabolic and protein-protein interaction (PPI) data from Arabidopsis thaliana (from Kyoto Encyclopedia of Genes and Genomes (KEGG) and BioGRID databases) and plant resistance curated knowledge for Viridiplantae obtained through text mining of the UniProt/SwissProt database. With the network, we integrated functional annotations and transcriptome data from two sugarcane genotypes that differ significantly in resistance to smut and applied a series of analyses to compare the transcriptomes and understand both signal perception and transduction in plant resistance. We show that the smut-resistant sugarcane has a larger arsenal of RGAs encompassing transcriptionally modulated subnetworks with other resistance elements, reaching hub proteins of primary metabolism. This approach may benefit molecular breeders in search of markers associated with quantitative resistance to diseases in non-model systems.

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.