Microcollinearity between autopolyploid sugarcane and diploid sorghum genomes

Identification of candidate single-nucleotide polymorphisms (SNPs) and genes associated with sugarcane leaf scald disease

Leaf scald, caused by Xanthomonas albilineans, is a severe disease affecting sugarcane worldwide. One of the most practical ways to control it is by developing resistant sugarcane cultivars. It is essential to identify genes associated with the response to leaf scald. A panel of 170 sugarcane genotypes was evaluated for resistance to leaf scald in field conditions for 2 years, followed by a 1-year greenhouse experiment. The phenotypic evaluation data showed a wide continuous distribution, with heritability values ranging from 0.58 to 0.84. Thirteen single nucleotide polymorphisms (SNPs) were identified, significantly associated with leaf scald resistance. Among these, eight were stable across multiple environments and association models. The candidate genes identified and validated based on RNA-seq and qRT-PCR included two genes that encode NB-ARC leucine-rich repeat (LRR)-containing domain disease-resistance protein. These findings provide a basis for developing marker-assisted selection strategies in sugarcane breeding programs.

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.

New methods for sorghum transformation in temperate climates

Sorghum (Sorghum bicolor) is an emerging cereal crop in temperate climates due to its high drought tolerance and other valuable traits. Genetic transformation is an important tool for the improvement of cereals. However, sorghum is recalcitrant to genetic transformation which is almost only successful in warmer climates. Here, we test the application of two new techniques for sorghum transformation in temperate climates, namely transient transformation by Agrobacterium tumefaciens-mediated agroinfiltration and stable transformation using gold particle bombardment and leaf whorls as explants. We optimized the transient transformation method, including post-infiltration incubation of plants in the dark and using Agrobacterium grown on plates with a high cell density (OD₆₀₀ = 2.0). Expression of the green fluorescence protein (GFP)-tagged endogenous sorghum gene SbDHR2 was achieved with low transformation efficiency, and our results point out a potential weakness in using this approach for localization studies. Furthermore, we succeeded in the production of callus and somatic embryos from leaf whorls, although no genetic transformation was accomplished with this method. Both methods show potential, even if they seem to be influenced by climatic conditions and therefore need further optimization to be applied routinely in temperate climates.

Transcriptomic and Proteomic Landscape of Sugarcane Response to Biotic and Abiotic Stressors

Sugarcane, a C₄ plant, provides most of the world's sugar, and a substantial amount of renewable bioenergy, due to its unique sugar-accumulating and feedstock properties. Brazil, India, China, and Thailand are the four largest sugarcane producers worldwide, and the crop has the potential to be grown in arid and semi-arid regions if its stress tolerance can be improved. Modern sugarcane cultivars which exhibit a greater extent of polyploidy and agronomically important traits, such as high sugar concentration, biomass production, and stress tolerance, are regulated by complex mechanisms. Molecular techniques have revolutionized our understanding of the interactions between genes, proteins, and metabolites, and have aided in the identification of the key regulators of diverse traits. This review discusses various molecular techniques for dissecting the mechanisms underlying the sugarcane response to biotic and abiotic stresses. The comprehensive characterization of sugarcane's response to various stresses will provide targets and resources for sugarcane crop improvement.

A Genome-Wide Association Study and Genomic Prediction for Fiber and Sucrose Contents in a Mapping Population of LCP 85-384 Sugarcane

Sugarcane (Saccharum spp. hybrids) is an economically important crop for both sugar and biofuel industries. Fiber and sucrose contents are the two most critical quantitative traits in sugarcane breeding that require multiple-year and multiple-location evaluations. Marker-assisted selection (MAS) could significantly reduce the time and cost of developing new sugarcane varieties. The objectives of this study were to conduct a genome-wide association study (GWAS) to identify DNA markers associated with fiber and sucrose contents and to perform genomic prediction (GP) for the two traits. Fiber and sucrose data were collected from 237 self-pollinated progenies of LCP 85-384, the most popular Louisiana sugarcane cultivar from 1999 to 2007. The GWAS was performed using 1310 polymorphic DNA marker alleles with three models of TASSEL 5, single marker regression (SMR), general linear model (GLM) and mixed linear model (MLM), and the fixed and random model circulating probability unification (FarmCPU) of R package. The results showed that 13 and 9 markers were associated with fiber and sucrose contents, respectively. The GP was performed by cross-prediction with five models, ridge regression best linear unbiased prediction (rrBLUP), Bayesian ridge regression (BRR), Bayesian A (BA), Bayesian B (BB) and Bayesian least absolute shrinkage and selection operator (BL). The accuracy of GP varied from 55.8% to 58.9% for fiber content and 54.6% to 57.2% for sucrose content. Upon validation, these markers can be applied in MAS and genomic selection (GS) to select superior sugarcane with good fiber and high sucrose contents.

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.

A draft chromosome-scale genome assembly of a commercial sugarcane

Sugarcane accounts for a large portion of the worlds sugar production. Modern commercial cultivars are complex hybrids of S. officinarum, S. spontaneum, and several other Saccharum species, resulting in an auto-allopolyploid with 8-12 copies of each chromosome. The current genome assembly gold standard is to generate a long read assembly followed by chromatin conformation capture sequencing to scaffold. We used the PacBio RSII and chromatin conformation capture sequencing to sequence and assemble the genome of a South East Asian commercial sugarcane cultivar, known as Khon Kaen 3. The Khon Kaen 3 genome assembled into 104,477 contigs totalling 7 Gb, which scaffolded into 56 pseudochromosomes containing 5.2 Gb of sequence. Genome annotation produced 242,406 genes from 30,927 orthogroups. Aligning the Khon Kaen 3 genome sequence to S. officinarum and S. spontaneum revealed a high level of apparent recombination, indicating a chimeric assembly. This assembly error is explained by high nucleotide identity between S. officinarum and S. spontaneum, where 91.8% of S. spontaneum aligns to S. officinarum at 94% identity. Thus, the subgenomes of commercial sugarcane are so similar that using short reads to correct long PacBio reads produced chimeric long reads. Future attempts to sequence sugarcane must take this information into account.

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

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

A genome-wide association and fine-mapping study of white rust resistance in hexaploid chrysanthemum cultivars with a wild diploid reference genome

White rust caused by Puccinia horiana is one of the most serious diseases of chrysanthemum (Chrysanthemum × morifolium). In this study, we report the DNA markers associated with resistance against P. horiana via a simple approach using the genome of a wild diploid relative, Chrysanthemum seticuspe. First, we identified the important region of the genome in the resistant cultivar "Ariesu" via a genome-wide association study. Simplex single nucleotide polymorphism (SNP) markers mined from ddRAD-Seq were used in a biparental population originating from crosses between resistant "Ariesu" and susceptible "Yellow Queen". The C. seticuspe genome was used as a reference. For the fine mapping of P. horiana resistance locus 2 (Phr2), a comparative whole genome sequencing study was conducted. Although the genome sequences of chrysanthemum cultivars assembled via the short-read approach were fragmented, reliable genome alignments were reconstructed by mapping onto the chromosome level of the C. seticuspe pseudomolecule. Base variants were then identified by comparing the assembled genome sequences of resistant "Ariesu" and susceptible "Yellow Queen". Consequently, SNP markers that were closer to Phr2 compared with ddRAD-Seq markers were obtained. These SNP markers co-segregated with resistance in F₁ progenies originating from resistant "Ariesu" and showed robust transferability for detecting Phr2-conferring resistance among chrysanthemum genetic resources. The wild C. seticuspe pseudomolecule, a de facto monoploid genome used for ddRAD-Seq analysis and assembled genome sequence comparison, demonstrated this method's utility as a model for developing DNA markers in hexaploid chrysanthemum cultivars.

Impact of Agroclimatic Variables on Proteogenomics in Sugar Cane (Saccharum spp.) Plant Productivity

Sugar cane (Saccharum spp. hybrids) is a major crop for sugar and renewable bioenergy worldwide, grown in arid and semiarid regions. China, the world's fourth-largest sugar producer after Brazil, India, and the European Union, all share ∼80% of the global production, and the remaining ∼20% of sugar comes from sugar beets, mostly grown in the temperate regions of the Northern Hemisphere, also used as a raw material in production of bioethanol for renewable energy. In view of carboxylation strategies, sugar cane qualifies as one of the best C₄ crop. It has dual CO₂ concentrating mechanisms located in its unique Krantz anatomy, having dimorphic chloroplasts located in mesophylls and bundle sheath cells for integrated operation of C₄ and C₃ carbon fixation cycles, regulated by enzymes to upgrade/sustain an ability for improved carbon assimilation to acquire an optimum carbon economy by producing enhanced plant biomass along with sugar yield under elevated temperature and strong irradiance with improved water-use efficiency. These superior intrinsic physiological carbon metabolisms encouraged us to reveal and recollect the facts for moving ahead with the molecular approaches to reveal the expression of proteogenomics linked with plant productivity under abiotic stress during its cultivation in specific agrizones globally.

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

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

Identification of potential MTAs and candidate genes for juice quality- and yield-related traits in Saccharum clones: a genome-wide association and comparative genomic study

Sugarcane is an economically important commercial crop which provides raw material for the production of sugar, jaggery, bioethanol, biomass and other by-products. Sugarcane breeding till today heavily relies on conventional breeding approaches which is time consuming, laborious and costly. Integration of marker-assisted selection (MAS) in sugarcane genetic improvement programs for difficult to select traits like sucrose content, resistance to pests and diseases and tolerance to abiotic stresses will accelerate varietal development. In the present study, association mapping approach was used to identify QTLs and genes associated with sucrose and other important yield-contributing traits. A mapping panel of 110 diverse sugarcane genotypes and 148 microsatellite primers were used for structured association mapping study. An optimal subpopulation number (ΔK) of 5 was identified by structure analysis. GWAS analysis using TASSEL identified a total of 110 MTAs which were localized into 27 QTLs by GLM and MLM (Q + K, PC + K) approaches. Among the 24 QTLs sequenced, 12 were able to identify potential candidate genes, viz., starch branching enzyme, starch synthase 4, sugar transporters and G3P-DH related to carbohydrate metabolism and hormone pathway-related genes ethylene insensitive 3-like 1, reversion to ethylene sensitive1-like, and auxin response factor associated to juice quality- and yield-related traits. Six markers, NKS 5_185, SCB 270_144, SCB 370_256, NKS 46_176 and UGSM 648_245, associated with juice quality traits and marker SMC31CUQ_304 associated with NMC were validated and identified as significantly associated to the traits by one-way ANOVA analysis. In conclusion, 24 potential QTLs identified in the present study could be used in sugarcane breeding programs after further validation in larger population. The candidate genes from carbohydrate and hormone response pathway presented in this study could be manipulated with genome editing approaches to further improve sugarcane crop.

Characterization of Repetitive DNA in Saccharum officinarum and Saccharum spontaneum by Genome Sequencing and Cytological Assays

In most plant species, DNA repeated elements such as satellites and retrotransposons are composing the majority of their genomes. Saccharum officinarum (2n = 8x = 80) and S. spontaneum (2n = 40-128) are the two fundamental donors of modern sugarcane cultivars. These two species are polyploids with large genome sizes and are enriched in repetitive elements. In this work, we adopted a de novo strategy to isolate highly repetitive and abundant sequences in S. officinarum LA Purple and S. spontaneum SES208. The findings obtained from alignment to the genome assemblies revealed that the vast majority of the repeats (97.9% in LA Purple and 96.5% in SES208) were dispersed in the respective genomes. Fluorescence in situ hybridization assays were performed on 27 representative repeats to investigate their distributions and abundances. The results showed that the copies of some highly repeated sequences, including rDNA and centromeric or telomeric repeats, were underestimated in current genome assemblies. The analysis of the raw read mapping strategy showed more copy numbers for all studied repeats, suggesting that copy number underestimation is common for highly repeated sequences in current genome assemblies of LA Purple and SES208. In addition, the data showed that the centromeric retrotransposons in all SES208 centromeres were absent in certain S. spontaneum clones with different ploidies. This rapid turnover of centromeric DNA in sugarcane provides new clues regarding the pattern of centromeric retrotransposon formation and accumulation.

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.

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.

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

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

Enhanced sugar accumulation and regulated plant hormone signalling genes contribute to cold tolerance in hypoploid Saccharum spontaneum

BACKGROUND

Wild sugarcane Saccharum spontaneum plants vary in ploidy, which complicates the utilization of its germplasm in sugarcane breeding. Investigations on cold tolerance in relation to different ploidies in S. spontaneum may promote the exploitation of its germplasm and accelerate the improvement of sugarcane varieties.

RESULTS

A hypoploid clone 12-23 (2n = 54) and hyperploid clone 15-28 (2n = 92) of S. spontaneum were analysed under cold stress from morphological, physiological, and transcriptomic perspectives. Compared with clone 15-28, clone 12-23 plants had lower plant height, leaf length, internode length, stem diameter, and leaf width; depressed stomata and prominent bristles and papillae; and thick leaves with higher bulliform cell groups and thicker adaxial epidermis. Compared with clone 15-28, clone 12-23 showed significantly lower electrical conductivity, significantly higher water content, soluble protein content, and superoxide dismutase activity, and significantly higher soluble sugar content and peroxidase activity. Under cold stress, the number of upregulated genes and downregulated genes of clone 12-23 was higher than clone 15-28, and many stress response genes and pathways were affected and enriched to varying degrees, particularly sugar and starch metabolic pathways and plant hormone signalling pathways. Under cold stress, the activity of 6-phosphate glucose trehalose synthase, trehalose phosphate phosphatase, and brassinosteroid-signalling kinase and the content of trehalose and brassinosteroids of clone 12-23 increased.

CONCLUSIONS

Compared with hyperploid clone 15-28, hypoploid clone 12-23 maintained a more robust osmotic adjustment system through sugar accumulation and hormonal regulation, which resulted in stronger cold tolerance.

Chromosome Pairing in Polyploid Grasses

Polyploids are species in which three or more sets of chromosomes coexist. Polyploidy frequently occurs in plants and plays a major role in their evolution. Based on their origin, polyploid species can be divided into two groups: autopolyploids and allopolyploids. The autopolyploids arise by multiplication of the chromosome sets from a single species, whereas allopolyploids emerge from the hybridization between distinct species followed or preceded by whole genome duplication, leading to the combination of divergent genomes. Having a polyploid constitution offers some fitness advantages, which could become evolutionarily successful. Nevertheless, polyploid species must develop mechanism(s) that control proper segregation of genetic material during meiosis, and hence, genome stability. Otherwise, the coexistence of more than two copies of the same or similar chromosome sets may lead to multivalent formation during the first meiotic division and subsequent production of aneuploid gametes. In this review, we aim to discuss the pathways leading to the formation of polyploids, the occurrence of polyploidy in the grass family (Poaceae), and mechanisms controlling chromosome associations during meiosis, with special emphasis on wheat.

Evaluation of the cold tolerance of Saccharum spontaneum L. clones with different ploidy levels on the basis of morphological and physiological indices

Saccharum spontaneum L. is one of the most important germplasm resources for modern sugarcane breeding. Exploring the cold tolerance of S. spontaneum clones with different ploidy levels and screening for cold-tolerant material can be helpful in parent selection for breeding cold-tolerant sugarcane. Morphological indices, leaf ultrastructure and physiological indices were used to evaluate the cold tolerance of 36 S. spontaneum clones with different ploidy levels (2n = 40, 48, 54, 60, 64, 78, 80, 88, 92 and 96). The morphological indices of S. spontaneum clones with different ploidy levels were positively correlated with ploidy. Under low-temperature stress, the chloroplast and mitochondrial structures of the clones with high ploidy were more severely damaged than were those of clones with low ploidy. A comprehensive evaluation of the physiological indices showed that the 36 S. spontaneum clones could be divided into four categories: strongly cold tolerant, cold tolerant, moderately cold tolerant and cold sensitive. Correlation analysis of the morphological indices and cold tolerance revealed a significant negative correlation between cold tolerance and ploidy. On the basis of the morphological and physiological indices, optimal stepwise regression equations that can be used for the selection of cold-tolerant S. spontaneum resources were established. The S. spontaneum clones with low ploidy are more cold tolerant than those with high ploidy. Clones 12-37, 13-10 and 12-23 are strongly cold-tolerant germplasm resources, which suggests these germplasm sources have high potential for use in breeding cold-tolerant sugarcane.

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.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Linkage Analysis and Haplotype Phasing in Experimental Autopolyploid Populations with High Ploidy Level Using Hidden Markov Models

Modern SNP genotyping technologies allow measurement of the relative abundance of different alleles for a given locus and consequently estimation of their allele dosage, opening a new road for genetic studies in autopolyploids. Despite advances in genetic linkage analysis in autotetraploids, there is a lack of statistical models to perform linkage analysis in organisms with higher ploidy levels. In this paper, we present a statistical method to estimate recombination fractions and infer linkage phases in full-sib populations of autopolyploid species with even ploidy levels for a set of SNP markers using hidden Markov models. Our method uses efficient two-point procedures to reduce the search space for the best linkage phase configuration and reestimate the final parameters by maximizing the likelihood of the Markov chain. To evaluate the method, and demonstrate its properties, we rely on simulations of autotetraploid, autohexaploid and autooctaploid populations and on a real tetraploid potato data set. The results show the reliability of our approach, including situations with complex linkage phase scenarios in hexaploid and octaploid populations.

Linkage Analysis and Haplotype Phasing in Experimental Autopolyploid Populations with High Ploidy Level Using Hidden Markov Models

Modern SNP genotyping technologies allow measurement of the relative abundance of different alleles for a given locus and consequently estimation of their allele dosage, opening a new road for genetic studies in autopolyploids. Despite advances in genetic linkage analysis in autotetraploids, there is a lack of statistical models to perform linkage analysis in organisms with higher ploidy levels. In this paper, we present a statistical method to estimate recombination fractions and infer linkage phases in full-sib populations of autopolyploid species with even ploidy levels for a set of SNP markers using hidden Markov models. Our method uses efficient two-point procedures to reduce the search space for the best linkage phase configuration and reestimate the final parameters by maximizing the likelihood of the Markov chain. To evaluate the method, and demonstrate its properties, we rely on simulations of autotetraploid, autohexaploid and autooctaploid populations and on a real tetraploid potato data set. The results show the reliability of our approach, including situations with complex linkage phase scenarios in hexaploid and octaploid populations.

Development of an Axiom Sugarcane100K SNP array for genetic map construction and QTL identification

An Axiom Sugarcane100K SNP array has been designed and successfully utilized to construct the sugarcane genetic map and to identify the QTLs associated with SCYLV resistance. To accelerate genetic studies in sugarcane, an Axiom Sugarcane100K single-nucleotide polymorphism (SNP) array was designed and customized in this study. Target enrichment sequencing 300 sugarcane accessions selected from the world collection of sugarcane and related grass species yielded more than four million SNPs, from which a total of 31,449 single-dose (SD) SNPs and 68,648 low-dosage (33,277 SD and 35,371 double dose) SNPs from two datasets, respectively, were selected and tiled on Affymetrix Axiom SNP array. Most of selected SNPs (91.77%) were located within genic regions (12,935 genes), with an average of 7.1 SNPs/gene according to sorghum gene models. This array was used to genotype 469 sugarcane clones, including one F₁ population derived from the cross between Green German and IND81-146, one selfing population derived from CP80-1827, and 11 diverse sugarcane accessions as controls. Results of genotyping revealed a high polymorphic SNP rate (77.04%) among the 469 samples. Three linkage maps were constructed by using SD SNP markers, including a genetic map for Green German with 3482 SD SNP markers spanning 3336 cM, a map for IND81-146 with 1513 SD SNP markers spanning 2615 cM, and a map for CP80-1827 with 536 SD SNP markers spanning 3651 cM. Quantitative trait loci (QTL) analysis identified 18 QTLs controlling Sugarcane yellow leaf virus resistance segregating in the two mapping populations, harboring 27 disease-resistant genes. This study demonstrated the successful development and utilization of a SNP array as an efficient genetic tool for high-throughput genotyping in highly polyploid sugarcane.

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

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

Assembly of allele-aware, chromosomal-scale autopolyploid genomes based on Hi-C data

Construction of chromosome-level assembly is a vital step in achieving the goal of a 'Platinum' genome, but it remains a major challenge to assemble and anchor sequences to chromosomes in autopolyploid or highly heterozygous genomes. High-throughput chromosome conformation capture (Hi-C) technology serves as a robust tool to dramatically advance chromosome scaffolding; however, existing approaches are mostly designed for diploid genomes and often with the aim of reconstructing a haploid representation, thereby having limited power to reconstruct chromosomes for autopolyploid genomes. We developed a novel algorithm (ALLHiC) that is capable of building allele-aware, chromosomal-scale assembly for autopolyploid genomes using Hi-C paired-end reads with innovative 'prune' and 'optimize' steps. Application on simulated data showed that ALLHiC can phase allelic contigs and substantially improve ordering and orientation when compared to other mainstream Hi-C assemblers. We applied ALLHiC on an autotetraploid and an autooctoploid sugar-cane genome and successfully constructed the phased chromosomal-level assemblies, revealing allelic variations present in these two genomes. The ALLHiC pipeline enables de novo chromosome-level assembly of autopolyploid genomes, separating each allele. Haplotype chromosome-level assembly of allopolyploid and heterozygous diploid genomes can be achieved using ALLHiC, overcoming obstacles in assembling complex genomes.

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

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

Genome-Wide Association Studies Identified Resistance Loci to Orange Rust and Yellow Leaf Virus Diseases in Sugarcane (Saccharum spp.)

Sugarcane (Saccharum spp.) supplies globally ∼80% of table sugar and 60% of bioethanol. Sugarcane orange rust and Sugarcane yellow leaf virus (SCYLV) are major sugarcane diseases, causing up to 50 and 40% yield losses, respectively. Sugarcane cultivars resistant to these diseases are needed to sustain sugarcane production in several regions. Dissecting DNA sequence variants controlling disease resistance provides a valuable tool for fulfilling a breeding strategy to develop resistant cultivars. In this study, we evaluated disease reactions to orange rust and SCYLV of a sugarcane diversity panel in repeated trials. We conducted a genome-wide association study between high-density markers and disease resistance reactions. We identified 91 putative DNA markers and 82 candidate genes significantly associated with resistance to one of the two diseases. These provide an important genetic resource for finding genes and molecular markers for disease resistance. Our results emphasized the importance of utilizing a wide germplasm collection for breeding resistant sugarcane cultivars.

Target enrichment sequencing of 307 germplasm accessions identified ancestry of ancient and modern hybrids and signatures of adaptation and selection in sugarcane (Saccharum spp.), a 'sweet' crop with 'bitter' genomes

Sugarcane (Saccharum spp.) is a highly energy-efficient crop primarily for sugar and bio-ethanol production. Sugarcane genetics and cultivar improvement have been extremely challenging largely due to its complex genomes with high polyploidy levels. In this study, we deeply sequenced the coding regions of 307 sugarcane germplasm accessions. Nearly five million sequence variations were catalogued. The average of 98× sequence depth enabled different allele dosages of sequence variation to be differentiated in this polyploid collection. With selected high-quality genome-wide SNPs, we performed population genomic studies and environmental association analysis. Results illustrated that the ancient sugarcane hybrids, S. barberi and S. sinense, and modern sugarcane hybrids are significantly different in terms of genomic compositions, hybridization processes and their potential ancestry contributors. Linkage disequilibrium (LD) analysis showed a large extent of LD in sugarcane, with 962.4 Kbp, 2739.2 Kbp and 3573.6 Kbp for S. spontaneum, S. officinarum and modern S. hybrids respectively. Candidate selective sweep regions and genes were identified during domestication and historical selection processes of sugarcane in addition to genes associated with environmental variables at the original locations of the collection. This research provided an extensive amount of genomic resources for sugarcane community and the in-depth population genomic analyses shed light on the breeding and evolution history of sugarcane, a highly polyploid species.

Whole chloroplast genome and gene locus phylogenies reveal the taxonomic placement and relationship of Tripidium (Panicoideae: Andropogoneae) to sugarcane

Background

For over 50 years, attempts have been made to introgress agronomically useful traits from Erianthus sect. Ripidium (Tripidium) species into sugarcane based on both genera being part of the ‘Saccharum Complex’, an interbreeding group of species believed to be involved in the origins of sugarcane. However, recent low copy number gene studies indicate that Tripidium and Saccharum are more divergent than previously thought. The extent of genus Tripidium has not been fully explored and many species that should be included in Tripidium are still classified as Saccharum. Moreover, Tripidium is currently defined as incertae sedis within the Andropogoneae, though it has been suggested that members of this genus are related to the Germainiinae.

Results

Eight newly-sequenced chloroplasts from potential Tripidium species were combined in a phylogenetic study with 46 members of the Panicoideae, including seven Saccharum accessions, two Miscanthidium and three Miscanthus species. A robust chloroplast phylogeny was generated and comparison with a gene locus phylogeny clearly places a monophyletic Tripidium clade outside the bounds of the Saccharinae. A key to the currently identified Tripidium species is presented.

Conclusion

For the first time, we have undertaken a large-scale whole plastid study of eight newly assembled Tripidium accessions and a gene locus study of five Tripidium accessions. Our findings show that Tripidium and Saccharum are 8 million years divergent, last sharing a common ancestor 12 million years ago. We demonstrate that four species should be removed from Saccharum/Erianthus and included in genus Tripidium. In a genome context, we show that Tripidium evolved from a common ancestor with and extended Germainiinae clade formed from Germainia, Eriochrysis, Apocopis, Pogonatherum and Imperata. We re-define the ‘Saccharum complex’ to a group of genera that can interbreed in the wild and extend the Saccharinae to include Sarga along with Sorghastrum, Microstegium vimineum and Polytrias (but excluding Sorghum). Monophyly of genus Tripidium is confirmed and the genus is expanded to include Tripidium arundinaceum, Tripidium procerum, Tripidium kanashiroi and Tripidium rufipilum. As a consequence, these species are excluded from genus Saccharum. Moreover, we demonstrate that genus Tripidium is distinct from the Germainiinae.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1356-9) contains supplementary material, which is available to authorized users.

Recent polyploidization events in three Saccharum founding species

The complexity of polyploid Saccharum genomes hindered progress of genome research and crop improvement in sugarcane. To understand their genome structure, transcriptomes of 59 F₁ individuals derived from S. officinarumLA Purple and S. robustum Molokai 5829 (2n = 80, x = 10 for both) were sequenced, yielding 11 157 and 8998 SNPs and 83 and 105 linkage groups, respectively. Most markers in each linkage group aligned to single sorghum chromosome. However, 71 interchromosomal rearrangements were detected between sorghum and S. officinarum or S. robustum, and 24 (33.8%) of them were shared between S. officinarum and S. robustum, indicating their occurrence before the speciation event that separated these two species. More than 2000 gene pairs from S. spontaneum, S. officinarum and S. robustum were analysed to estimate their divergence time. Saccharum officinarum and S. robustum diverged about 385 thousand years ago, and the whole-genome duplication events occurred after the speciation event because of shared interchromosomal rearrangements. The ancestor of these two species diverged from S. spontaneum about 769 thousand years ago, and the reduction in basic chromosome number from 10 to 8 in S. spontaneum occurred after the speciation event but before the two rounds of whole-genome duplication. Our results proved that S. officinarum is a legitimate species in its own right and not a selection from S. robustum during the domestication process in the past 10 000 years. Our findings rejected a long-standing hypothesis and clarified the timing of speciation and whole-genome duplication events in Saccharum.

Genome-wide association mapping identifies markers associated with cane yield components and sucrose traits in the Louisiana sugarcane core collection

Sugarcane is an economically important crop for both food and biofuel industries. Marker-assisted breeding in sugarcane is becoming a reality with the recent development and deployment of markers linked with disease resistance genes. Large linkage disequilibrium in sugarcane makes genome-wide association studies (GWAS) a better alternative to biparental mapping to identify markers associated with agronomic traits. GWAS was conducted on a Louisiana core collection to identify marker-trait associations (MTA) for 11 cane yield and sucrose traits using single nucleotide polymorphism (SNP) and insertion-deletion (Indel) markers. Significant (P < .05) MTAs were identified for all traits where the top ranked markers explained up to 15% of the total phenotypic variation. High correlations (0.732 to 0.999) were observed between sucrose traits and 56 markers were found consistent across multiple traits. These markers following validation in more diverse populations could be used in marker-assisted selection of clones in sugarcane breeding program in Louisiana and elsewhere.

Sweet Sorghum Originated through Selection of Dry, a Plant-Specific NAC Transcription Factor Gene

Sorghum (Sorghum bicolor) is the fifth most popular crop worldwide and a C₄ model plant. Domesticated sorghum comes in many forms, including sweet cultivars with juicy stems and grain sorghum with dry, pithy stems at maturity. The Dry locus, which controls the pithy/juicy stem trait, was discovered over a century ago. Here, we found that Dry gene encodes a plant-specific NAC transcription factor. Dry was either deleted or acquired loss-of-function mutations in sweet sorghum, resulting in cell collapse and altered secondary cell wall composition in the stem. Twenty-three Dry ancestral haplotypes, all with dry, pithy stems, were found among wild sorghum and wild sorghum relatives. Two of the haplotypes were detected in domesticated landraces, with four additional dry haplotypes with juicy stems detected in improved lines. These results imply that selection for Dry gene mutations was a major step leading to the origin of sweet sorghum. The Dry gene is conserved in major cereals; fine-tuning its regulatory network could provide a molecular tool to control crop stem texture.

Comparative Analysis of Homologous Sequences of Saccharum officinarum and Saccharum spontaneum Reveals Independent Polyploidization Events

Sugarcane (Saccharum spp. hybrids) is an economically important crop widely grown in tropical and subtropical regions for sugar and ethanol production. However, the large genome size, high ploidy level, interspecific hybridization and aneuploidy make sugarcane one of the most complex genomes and have long hampered genome research in sugarcane. Modern sugarcane cultivars are derived from interspecific hybridization between S. officinarum and S. spontaneum with 80-90% of the genome from S. officinarum and 10-20% of the genome from S. spontaneum. We constructed bacterial artificial chromosome (BAC) libraries of S. officinarum variety LA Purple (2n = 8x = 80) and S. spontaneum haploid clone AP85-441 (2n = 4x = 32), and selected and sequenced 97 BAC clones from the two Saccharum BAC libraries. A total of 5,847,280 bp sequence from S. officinarum and 5,011,570 bp from S. spontaneum were assembled and 749 gene models were annotated in these BACs. A relatively higher gene density and lower repeat content were observed in S. spontaneum BACs than in S. officinarum BACs. Comparative analysis of syntenic regions revealed a high degree of collinearity in genic regions between Saccharum and Sorghum bicolor and between S. officinarum and S. spontaneum. In the syntenic regions, S. spontaneum showed expansion relative to S. officinarum, and both S. officinarum and S. spontaneum showed expansion relative to sorghum. Among the 75 full-length LTR retrotransposons identified in the Saccharum BACs, none of them are older than 2.6 mys and no full-length LTR elements are shared between S. officinarum and S. spontaneum. In addition, divergence time estimated using a LTR junction marker and a syntenic gene shared by 3 S. officinarum and 1 S. spontaneum BACs revealed that the S. spontaneum intergenic region was distant to those from the 3 homologous regions in S. officinarum. Our results suggested that S. officinarum and S. spontaneum experienced at least two rounds of independent polyploidization in each lineage after their divergence from a common ancestor.

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

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

Identification of Genomic Regions Controlling Leaf Scald Resistance in Sugarcane Using a Bi-parental Mapping Population and Selective Genotyping by Sequencing

Leaf scald, caused by Xanthomonas albilineans, is a major sugarcane disease worldwide. The disease is managed primarily with resistant cultivars obtained through classical breeding. However, erratic symptom expression hinders the reliability and reproducibility of selection for resistance. The development and use of molecular markers associated with incompatible/compatible reactions could overcome this limitation. The aim of the present work was to find leaf scald resistance-associated molecular markers in sugarcane to facilitate marker-assisted breeding. A genetic linkage map was constructed by selective genotyping of 89 pseudo F2 progenies of a cross between LCP 85-384 (resistant) and L 99-226 (susceptible) using 1,948 single dose (SD) markers generated from SSR, eSSR, and SNPs. Of these, 1,437 SD markers were mapped onto 294 linkage groups, which covered 19,464 cM with 120 and 138 LGs assigned to the resistant and susceptible parent, respectively. Composite interval mapping identified 8 QTLs associated with the disease response with LOD scores ranging from 3.0 to 7.6 and explained 5.23 to 16.93% of the phenotypic variance. Comparative genomics analysis with Sorghum bicolor allowed us to pinpoint three SNP markers that explained 16% phenotypic variance. In addition, representative stress-responsive genes close to the major effect QTLs showed upregulation in their expression in response to the bacterial infection in leaf/meristem tissue.

Transcriptomic analysis of transgressive segregants revealed the central role of photosynthetic capacity and efficiency in biomass accumulation in sugarcane

Sugarcane is among the most efficient crops in converting solar energy into chemical energy. However, due to its complex genome structure and inheritance, the genetic and molecular basis of biomass yield in sugarcane is still largely unknown. We created an F2 segregating population by crossing S. officinarum and S. spontaneum and evaluated the biomass yield of the F2 individuals. The F2 individuals exhibited clear transgressive segregation in biomass yield. We sequenced transcriptomes of source and sink tissues from 12 selected extreme segregants to explore the molecular basis of high biomass yield for future breeding of high-yielding energy canes. Among the 103,664 assembled unigenes, 10,115 and 728 showed significant differential expression patterns between the two extreme segregating groups in the top visible dewlap leaf and the 9^th culm internode, respectively. The most enriched functional categories were photosynthesis and fermentation in the high-biomass and the low-biomass groups, respectively. Our results revealed that high-biomass yield was mainly determined by assimilation of carbon in source tissues. The high-level expression of fermentative genes in the low-biomass group was likely induced by their low-energy status. Group-specific expression alleles which can be applied in the development of new high-yielding energy cane varieties via molecular breeding were identified.

Diversity of the Bacterial Microbiome in the Roots of Four Saccharum Species: S. spontaneum, S. robustum, S. barberi, and S. officinarum

Endophytic bacteria are nearly ubiquitously present in the internal tissues of plants, and some endophytes can promote plant growth. In this study, we sampled the roots of four ancestral species of sugarcane (two genotypes per species) and two sugarcane cultivars, and used 16S rRNA and nifH gene sequencing to characterize the root endophytic bacterial communities and diazotroph diversity. A total of 7,198 operational taxonomic units (OTUs) were detected for the endophytic bacteria community. The endophytic bacterial communities exhibited significantly different α- and β-diversities. From the 202 detected families in the sugarcane roots, a core microbiome containing 13 families was identified. The nifH gene was successfully detected in 9 of 30 samples from the four sugarcane species assayed, and 1,734 OTUs were merged for endophytic diazotrophs. In the tested samples, 43 families of endophytic diazotrophs were detected, and six families showed differences across samples. Among the 20 most abundant detected genera, 10 have been reported to be involved in nitrogen fixation in sugarcane. These findings demonstrate the diversity of the microbial communities in different sugarcane germplasms and shed light on the mechanism of biological nitrogen fixation in sugarcane.

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

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

Transcriptomic characterization and potential marker development of contrasting sugarcane cultivars

Sugarcane (Saccharum officinarum L.) is an important crop for sugar production and bioenergy worldwide. In this study, we performed transcriptome sequencing for six contrasting sugarcane genotypes involved in leaf abscission, tolerance to pokkah boeng disease and drought stress. More than 465 million high-quality reads were generated, which were de novo assembled into 93,115 unigenes. Based on a similarity search, 43,526 (46.74%) unigenes were annotated against at least one of the public databases. Functional classification analyses showed that these unigenes are involved in a wide range of metabolic pathways. Comparative transcriptome analysis revealed that many unigenes involved in response to abscisic acid and ethylene were up-regulated in the easy leaf abscission genotype, and unigenes associated with response to jasmonic acid and salicylic acid were up-regulated in response to the pokkah boeng disease in the tolerance genotype. Moreover, unigenes related to peroxidase, antioxidant activity and signal transduction were up-regulated in response to drought stress in the tolerant genotype. Finally, we identified a number of putative markers, including 8,630 simple sequence repeats (SSRs) and 442,152 single-nucleotide polymorphisms (SNPs). Our data will be important resources for future gene discovery, molecular marker development, and genome studies in sugarcane.

The Challenge of Analyzing the Sugarcane Genome

Reference genome sequences have become key platforms for genetics and breeding of the major crop species. Sugarcane is probably the largest crop produced in the world (in weight of crop harvested) but lacks a reference genome sequence. Sugarcane has one of the most complex genomes in crop plants due to the extreme level of polyploidy. The genome of modern sugarcane hybrids includes sub-genomes from two progenitors Saccharum officinarum and S. spontaneum with some chromosomes resulting from recombination between these sub-genomes. Advancing DNA sequencing technologies and strategies for genome assembly are making the sugarcane genome more tractable. Advances in long read sequencing have allowed the generation of a more complete set of sugarcane gene transcripts. This is supporting transcript profiling in genetic research. The progenitor genomes are being sequenced. A monoploid coverage of the hybrid genome has been obtained by sequencing BAC clones that cover the gene space of the closely related sorghum genome. The complete polyploid genome is now being sequenced and assembled. The emerging genome will allow comparison of related genomes and increase understanding of the functioning of this polyploidy system. Sugarcane breeding for traditional sugar and new energy and biomaterial uses will be enhanced by the availability of these genomic resources.

Comprehensively Characterizing the Cytological Features of Saccharum spontaneum by the Development of a Complete Set of Chromosome-Specific Oligo Probes

Chromosome-specific identification is a powerful technique in the study of genome structure and evolution. However, there is no reliable cytogenetic marker to unambiguously identify each of the chromosomes in sugarcane (Saccharum spp., Poaceae), which has a complex genome with a high level of ploidy and heterozygosity. In this study, we developed a set of oligonucleotide (oligo)-based probes through bioinformatic design and massive synthetization. These probes produced a clear and bright single signal in each of the chromosomes and their eight homologous chromosomes in the ancient species Saccharum spontaneum (2n = 8x = 64). Thus, they can be used as reliable markers to robustly label each of the chromosomes in S. spontaneum. We then obtained the karyotype data and established a nomenclature based on chromosomal sizes for the eight chromosomes of the octoploid S. spontaneum. In addition, we also found that the 45S and 5S rDNAs demonstrated high copy number variations among different homologous chromosomes, indicating a rapid evolution of the highly repeated sequence after polyploidization. Our fluorescence in situ hybridization (FISH) assay also demonstrated that these probes could be used as cross-species markers between or within the genera of Sorghum and Saccharum. By comparing FISH analyses, we discovered that several chromosome rearrangement events occurred in S. spontaneum, which might have contributed to the basic chromosome number reduction from 10 in sorghum to 8 in sugarcane. Consistent identification of individual chromosomes makes molecular cytogenetic study possible in sugarcane and will facilitate fine chromosomal structure and karyotype evolution of the genus Saccharum.

Development and Applications of Chromosome-Specific Cytogenetic BAC-FISH Probes in S. spontaneum

Saccharum spontaneum is a major Saccharum species that contributed to the origin of modern sugarcane cultivars, and due to a high degree of polyploidy is considered to be a plant species with one of the most complex genetics. Fluorescence in situ hybridization (FISH) is a powerful and widely used tool in genome studies. Here, we demonstrated that FISH based on bacterial artificial chromosome (BAC) clones can be used as a specific cytological marker to identify S. spontaneum individual chromosomes and study the relationship between S. spontaneum and other related species. We screened low-copy BACs as probes from the sequences of a high coverage of S. spontaneum BAC library based on BLAST search of the sorghum genome. In total, we isolated 49 positive BAC clones, and identified 27 BAC clones that can give specific signals on the S. spontaneum chromosomes. Of the 27 BAC probes, 18 were confirmed to be able to discriminate the eight basic chromosomes of S. spontaneum. Moreover, BAC-24, BAC-66, BAC-78, BAC-69, BAC-71, BAC-73, and BAC-77 probes were used to construct physical maps of chromosome 1 and chromosome 2 of S. spontaneum, which indicated synteny in Sb01 between S. spontaneum and sorghum. Furthermore, we found that BAC-14 and BAC-19 probes, corresponding to the sorghum chromosomes 2 and 8, respectively, localized to different arms of the same S. spontaneum chromosome, suggesting that there was an inter-chromosomal rearrangement event between S. spontaneum and sorghum. Our study provides the first set of chromosome-specific cytogenetic markers in Saccharum and is critical for future advances in cytogenetics and genome sequencing studies in Saccharum.

Analysis of Three Sugarcane Homo/Homeologous Regions Suggests Independent Polyploidization Events of Saccharum officinarum and Saccharum spontaneum

Whole genome duplication has played an important role in plant evolution and diversification. Sugarcane is an important crop with a complex hybrid polyploid genome, for which the process of adaptation to polyploidy is still poorly understood. In order to improve our knowledge about sugarcane genome evolution and the homo/homeologous gene expression balance, we sequenced and analyzed 27 BACs (Bacterial Artificial Chromosome) of sugarcane R570 cultivar, containing the putative single-copy genes LFY (seven haplotypes), PHYC (four haplotypes), and TOR (seven haplotypes). Comparative genomic approaches showed that these sugarcane loci presented a high degree of conservation of gene content and collinearity (synteny) with sorghum and rice orthologous regions, but were invaded by transposable elements (TE). All the homo/homeologous haplotypes of LFY, PHYC, and TOR are likely to be functional, because they are all under purifying selection (dN/dS ≪ 1). However, they were found to participate in a nonequivalently manner to the overall expression of the corresponding gene. SNPs, indels, and amino acid substitutions allowed inferring the S. officinarum or S. spontaneum origin of the TOR haplotypes, which further led to the estimation that these two sugarcane ancestral species diverged between 2.5 and 3.5 Ma. In addition, analysis of shared TE insertions in TOR haplotypes suggested that two autopolyploidization may have occurred in the lineage that gave rise to S. officinarum, after its divergence from S. spontaneum.

Mining sequence variations in representative polyploid sugarcane germplasm accessions

Background

Sugarcane (Saccharum spp.) is one of the most important economic crops because of its high sugar production and biofuel potential. Due to the high polyploid level and complex genome of sugarcane, it has been a huge challenge to investigate genomic sequence variations, which are critical for identifying alleles contributing to important agronomic traits. In order to mine the genetic variations in sugarcane, genotyping by sequencing (GBS), was used to genotype 14 representative Saccharum complex accessions. GBS is a method to generate a large number of markers, enabled by next generation sequencing (NGS) and the genome complexity reduction using restriction enzymes.

Results

To use GBS for high throughput genotyping highly polyploid sugarcane, the GBS analysis pipelines in 14 Saccharum complex accessions were established by evaluating different alignment methods, sequence variants callers, and sequence depth for single nucleotide polymorphism (SNP) filtering. By using the established pipeline, a total of 76,251 non-redundant SNPs, 5642 InDels, 6380 presence/absence variants (PAVs), and 826 copy number variations (CNVs) were detected among the 14 accessions. In addition, non-reference based universal network enabled analysis kit and Stacks de novo called 34,353 and 109,043 SNPs, respectively. In the 14 accessions, the percentages of single dose SNPs ranged from 38.3% to 62.3% with an average of 49.6%, much more than the portions of multiple dosage SNPs. Concordantly called SNPs were used to evaluate the phylogenetic relationship among the 14 accessions. The results showed that the divergence time between the Erianthus genus and the Saccharum genus was more than 10 million years ago (MYA). The Saccharum species separated from their common ancestors ranging from 0.19 to 1.65 MYA.

Conclusions

The GBS pipelines including the reference sequences, alignment methods, sequence variant callers, and sequence depth were recommended and discussed for the Saccharum complex and other related species. A large number of sequence variations were discovered in the Saccharum complex, including SNPs, InDels, PAVs, and CNVs. Genome-wide SNPs were further used to illustrate sequence features of polyploid species and demonstrated the divergence of different species in the Saccharum complex. The results of this study showed that GBS was an effective NGS-based method to discover genomic sequence variations in highly polyploid and heterozygous species.

Evolution and expression of the fructokinase gene family in Saccharum

BACKGROUND

Sugarcane is an important sugar crop contributing up to about 80% of the world sugar production. Efforts to characterize the genes involved in sugar metabolism at the molecular level are growing since increasing sugar content is a major goal in the breeding of new sugarcane varieties. Fructokinases (FRK) are the main fructose phosphorylating enzymes with high substrate specificity and affinity.

RESULTS

In this study, by combining comparative genomics approaches with BAC resources, seven fructokinase genes were identified in S. spontaneum. Phylogenetic analysis based on representative monocotyledon and dicotyledon plant species suggested that the FRK gene family is ancient and its evolutionary history can be traced in duplicated descending order: SsFRK4, SsFRK6/SsFRK7,SsFRK5, SsFRK3 and SsFRK1/SsFRK2. Among the close orthologs, the number and position of exons in FRKs were conserved; in contrast, the size of introns varied among the paralogous FRKs in Saccharum. Genomic constraints were analyzed within the gene alleles and between S. spontaneum and Sorghum bicolor, and gene expression analysis was performed under drought stress and with exogenous applications of plant hormones. FRK1, which was under strong functional constraint selection, was conserved among the gene allelic haplotypes, and displayed dominant expression levels among the gene families in the control conditions, suggesting that FRK1 plays a major role in the phosphorylation of fructose. FRK3 and FRK5 were dramatically induced under drought stress, and FRK5 was also found to increase its expression levels in the mature stage of Saccharum. Similarly, FRK3 and FRK5 were induced in response to drought stress in Saccharum. FRK2 and FRK7 displayed lower expression levels than the other FRK family members; FRK2 was under strong genomic selection constraints whereas FRK7 was under neutral selection. FRK7 may have become functionally redundant in Saccharum through pseudogenization. FRK4 and FRK6 shared the most similar expression pattern: FRK4 was revealed to have higher expression levels in mature tissues than in premature tissues of Saccharum, and FRK6 presented a slight increase under drought stress.

CONCLUSIONS

Our study presents a comprehensive genomic study of the entire FRK gene family in Saccharum, providing the foundations for approaches to characterize the molecular mechanism regulated by the SsFRK family in sugarcane.

Isolation and characterization of centromeric repetitive DNA sequences in Saccharum spontaneum

Sugarcane (Saccharum hybrids spp.) is the most important sugar crop that accounts for ~75% of the world’s sugar production. Recently, a whole-genome sequencing project was launched on the wild species S. spontaneum. To obtain information on the DNA composition of the repeat-enriched region of the centromere, we conducted a genome-wide analysis of the DNA sequences associated with CenH3 (a mutant of histone H3 located in eukaryote centromeres) using chromatin immunoprecipitation followed by sequencing (ChIP-seq) method. We demonstrate that the centromeres contain mainly SCEN-like single satellite repeat (Ss1) and several Ty3/gypsy retrotransposon-related repeats (Ss166, Ss51, and Ss68). Ss1 dominates in the centromeric regions and spans up to 500 kb. In contrast, the Ty3/gypsy retrotransposon-related repeats are either clustered spanning over a short range, or dispersed in the centromere regions. Interestingly, Ss1 exhibits a chromosome-specific enrichment in the wild species S. spontaneum and S. robustum, but not in the domesticated species S. officinarum and modern sugarcane cultivars. This finding suggests an autopolyploid genome identity of S. spontaneum with a high level of homology among its eight sub-genomes. We also conducted a genome-wide survey of the repetitive DNAs in S. spontaneum following a similarity-based sequence clustering strategy. These results provide insight into the composition of sugarcane genome as well as the genome assembly of S. spontaneum.