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

Genetic association analysis in sugarcane (Saccharum spp.) for sucrose accumulation in humid environments in Colombia

BACKGROUND

Sucrose accumulation in sugarcane is affected by several environmental and genetic factors, with plant moisture being of critical importance for its role in the synthesis and transport of sugars within the cane stalks, affecting the sucrose concentration. In general, rainfall and high soil humidity during the ripening stage promote plant growth, increasing the fresh weight and decreasing the sucrose yield in the humid region of Colombia. Therefore, this study aimed to identify markers associated with sucrose accumulation or production in the humid environment of Colombia through a genome-wide association study (GWAS).

RESULTS

Sucrose concentration measurements were taken in 220 genotypes from the Cenicaña's diverse panel at 10 (early maturity) and 13 (normal maturity) months after planting. For early maturity data was collected during plant cane and first ratoon, while at normal maturity it was during plant cane, first, and second ratoon. A total of 137,890 SNPs were selected after sequencing the 220 genotypes through GBS, RADSeq, and whole-genome sequencing. After GWAS analysis, a total of 77 markers were significantly associated with sucrose concentration at both ages, but only 39 were close to candidate genes previously reported for sucrose accumulation and/or production. Among the candidate genes, 18 were highlighted because they were involved in sucrose hydrolysis (SUS6, CIN3, CINV1, CINV2), sugar transport (i.e., MST1, MST2, PLT5, SUT4, ERD6 like), phosphorylation processes (TPS genes), glycolysis (PFP-ALPHA, HXK3, PHI1), and transcription factors (ERF12, ERF112). Similarly, 64 genes were associated with glycosyltransferases, glycosidases, and hormones.

CONCLUSIONS

These results provide new insights into the molecular mechanisms involved in sucrose accumulation in sugarcane and contribute with important genomic resources for future research in the humid environments of Colombia. Similarly, the markers identified will be validated for their potential application within Cenicaña's breeding program to assist the development of breeding populations.

Sugarcane breeding: a fantastic past and promising future driven by technology and methods

Sugarcane is the most important sugar and energy crop in the world. During sugarcane breeding, technology is the requirement and methods are the means. As we know, seed is the cornerstone of the development of the sugarcane industry. Over the past century, with the advancement of technology and the expansion of methods, sugarcane breeding has continued to improve, and sugarcane production has realized a leaping growth, providing a large amount of essential sugar and clean energy for the long-term mankind development, especially in the face of the future threats of world population explosion, reduction of available arable land, and various biotic and abiotic stresses. Moreover, due to narrow genetic foundation, serious varietal degradation, lack of breakthrough varieties, as well as long breeding cycle and low probability of gene polymerization, it is particularly important to realize the leapfrog development of sugarcane breeding by seizing the opportunity for the emerging Breeding 4.0, and making full use of modern biotechnology including but not limited to whole genome selection, transgene, gene editing, and synthetic biology, combined with information technology such as remote sensing and deep learning. In view of this, we focus on sugarcane breeding from the perspective of technology and methods, reviewing the main history, pointing out the current status and challenges, and providing a reasonable outlook on the prospects of smart breeding.

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

Background

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

Methods

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

Results

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

Conclusion

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

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

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

Validation of candidate gene-based EST-SSR markers for sugar yield in sugarcane

Sugarcane (Saccharum spp.) is a widely cultivated crop that fulfils approximately 75% of the sucrose demand worldwide. Owing to its polyploidy and complex genetic nature, it is difficult to identify and map genes related to complex traits, such as sucrose content. However, association mapping is one of the alternatives for identifying genes or markers for marker-assisted selection. In the present study, EST-SSR primers were obtained from in silico studies. The functionality of each primer was tested using Blast2Go software, and 30 EST-SSR primers related to sugar content were selected. These markers were validated using association analysis. A total of 70 F1 diverse genotypes for sugar content were phenotypes with two check lines. All parameters related to sugar content were recorded. The results showed a significant variation between the genotypes for sugar yield traits such as Brix value, purity, and sucrose content, etc. Correlation studies revealed that the Brix%, sucrose content, and sucrose recovery were significantly correlated. An association analysis was performed using mixed linear model to avoid false positive associations. The association analysis revealed that the SEM 407 marker was significantly associated with Brix% and sucrose content. The SEM 407 primers are putatively related to diphosphate-fructose-6-phosphate 1-phosphotransferase which is associated with Brix% and sucrose content. This functional marker can be used for marker-assisted selection for sugar yield traits in sugarcane that could accelerate the sugarcane breeding program.

Genetic diversity and population structure of Saccharum hybrids

Sugarcane breeding programs incorporate foreign material to broaden the genetic base, expanding the gene pool. In South America, the Inter-university Network for the Development of the Sugarcane Industry (RIDESA) and Estación Experimental Agroindustrial Obispo Colombres (EEAOC) sugarcane breeding programs from Brazil and Argentina, respectively, have never exchanged materials. In that sense, the knowledge of the genetic diversity and population structure among sugarcane genotypes of both germplasm banks, determined in a reliable way through their molecular profiles, will provide valuable information to select the best parental accessions for crossing aimed at the efficient introgression of desirable alleles. For that, the aim was to determine the genetic diversity and population structure of 96 Saccharum commercial hybrids from RIDESA and EEAOC sugarcane breeding programs by using TRAP, SSR and markers related to disease resistance (e.g. Bru1 and G1). Genetic structure was determined through genetic similarity analysis, analysis of molecular variance (AMOVA), Multidimensional scaling (MDS), and a Bayesian method. Average PIC values were 0.25 and 0.26, Ho values were 0.24 and 0.28, and He values were 0.25 and 0.28, for TRAP and SSR primers, respectively. Genetic similarity, MDS, and analysis of structure revealed that Brazilian and Argentinean genotypes clustered in two groups clearly differentiated, whereas AMOVA suggested that there is more variability within programs than between them. Regarding Bru1 markers, Brazilian genotypes showed high frequency of haplotype 1 (71.4%) whereas Argentinean genotypes showed high frequency of haplotype 4 (80.8%); haplotypes 1 and 4 are indicated for the presence and absence of the brown rust resistance gene (Bru1), respectively. Respecting the G1 marker, most of the evaluated genotypes (60.4%) showed the presence of the fragment, in a similar proportion for genotypes of both programs. In conclusion, the exchange of materials, at least the most diverse genotypes, between RIDESA and EEAOC breeding programs will allow extending the genetic base of their germplasm banks, and the knowledge of genetic diversity will help breeders to better manage crosses, increasing the probability of obtaining more productive varieties.

A first look at the ability to use genomic prediction for improving the ratooning ability of sugarcane

The sugarcane ratooning ability (RA) is the most important target trait for breeders seeking to enhance the profitability of sugarcane production by reducing the planting cost. Understanding the genetics governing the RA could help breeders by identifying molecular markers that could be used for genomics-assisted breeding (GAB). A replicated field trial was conducted for three crop cycles (plant cane, first ratoon, and second ratoon) using 432 sugarcane clones and used for conducting genome-wide association and genomic prediction of five sugar and yield component traits of the RA. The RA traits for economic index (EI), stalk population (SP), stalk weight (SW), tonns of cane per hectare (TCH), and tonns of sucrose per hectare (TSH) were estimated from the yield and sugar data. A total of six putative quantitative trait loci and eight nonredundant single-nucleotide polymorphism (SNP) markers were associated with all five tested RA traits and appear to be unique. Seven putative candidate genes were colocated with significant SNPs associated with the five RA traits. The genomic prediction accuracies for those tested traits were moderate and ranged from 0.21 to 0.36. However, the models fitting fixed effects for the most significant associated markers for each respective trait did not give any advantages over the standard models without fixed effects. As a result of this study, more robust markers could be used in the future for clone selection in sugarcane, potentially helping resolve the genetic control of the RA in sugarcane.

Combination of GWAS and FST-based approaches identified loci associated with economic traits in sugarcane

Sugarcane is a globally important plant for both sugar and biofuel production. Although conventional breeding has played an important role in increasing the productivity of sugarcane, it takes a long time to achieve breeding goals such as high yield and resistant to diseases. Molecular breeding, including marker-assisted breeding and genomic selection, can accelerate genetic improvement by selecting elites at the seedling stage with DNA markers. However, only a few DNA markers associated with important traits were identified in sugarcane. The purpose of this study was to identify DNA markers associated with sugar content, stalk diameter, and sugarcane top borer resistance. The sugarcane samples with trait records were genotyped using the restriction site-associated DNA sequencing (RADseq) technology. Using F_ST analysis and genome-wide association study (GWAS), a total of 9, 23 and 9 DNA variants (single nucleotide polymorphisms (SNPs)/insertions and deletions (indels)) were associated with sugar content, stalk diameter, and sugarcane top borer resistance, respectively. The identified genetic variants were on different chromosomes, suggesting that these traits are complex and determined by multiple genetic factors. These DNA markers identified by both approaches have the potential to be used in selecting elite clones at the seeding stage in our sugarcane breeding program to accelerate genetic improvement. Certainly, it is essential to verify the reliability of the identified DNA markers associated with traits before they are used in molecular breeding in other populations.

QTL mapping and identification of candidate genes linked to red rot resistance in sugarcane

Sugarcane (Saccharum species hybrid) is one of the most important commercial crops cultivated worldwide for products like white sugar, bagasse, ethanol, etc. Red rot is a major sugarcane disease caused by a hemi-biotrophic fungus, Colletotrichum falcatum Went., which can potentially cause a reduction in yield up to 100%. Breeding for red rot-resistant sugarcane varieties has become cumbersome due to its complex genome and frequent generation of new pathotypes of red rot fungus. In the present study, a genetic linkage map was developed using a selfed population of a popular sugarcane variety CoS 96268. A QTL linked to red rot resistance (qREDROT) was identified, which explained 26% of the total phenotypic variation for the trait. A genotype-phenotype network analysis performed to account for epistatic interactions, identified the key markers involved in red rot resistance. The differential expression of the genes located in the genomic region between the two flanking markers of the qREDROT as well as in the vicinity of the markers identified through the genotype-phenotype network analysis in a set of contrasting genotypes for red rot infection further confirmed the mapping results. Further, the expression analysis revealed that the plant defense-related gene coding 26S protease regulatory subunit is strongly associated with the red rot resistance. The findings can help in the screening of disease resistant genotypes for developing red rot-resistant varieties of sugarcane.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03481-7.

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.

Identification of QTLs and critical genes related to sugarcane mosaic disease resistance

Mosaic viral diseases affect sugarcane productivity worldwide. Mining disease resistance-associated molecular markers or genes is a key component of disease resistance breeding programs. In the present study, 285 F₁ progeny were produced from a cross between Yuetang 93-159, a moderately resistant variety, and ROC22, a highly susceptible variety. The mosaic disease symptoms of these progenies, with ROC22 as the control, were surveyed by natural infection under 11 different environmental conditions in the field and by artificial infections with a mixed sugarcane mosaic virus (SCMV) and sorghum mosaic virus (SrMV) inoculum. Analysis of consolidated survey data enabled the identification of 29 immune, 55 highly resistant, 70 moderately resistant, 62 susceptible, and 40 highly susceptible progenies. The disease response data and a high-quality SNP genetic map were used in quantitative trait locus (QTL) mapping. The results showed that the correlation coefficients (0.26~0.91) between mosaic disease resistance and test environments were significant (p< 0.001), and that mosaic disease resistance was a highly heritable quantitative trait (H² = 0.85). Seven mosaic resistance QTLs were located to the SNP genetic map, each QTL accounted for 3.57% ~ 17.10% of the phenotypic variation explained (PVE). Furthermore, 110 pathogen response genes and 69 transcription factors were identified in the QTLs interval. The expression levels of nine genes (Soffic.07G0015370-1P, Soffic.09G0015410-2T, Soffic.09G0016460-1T, Soffic.09G0016460-1P, Soffic.09G0017080-3C, Soffic.09G0018730-3P, Soffic.09G0018730-3C, Soffic.09G0019920-3C and Soffic.03G0019710-2C) were significantly different between resistant and susceptible progenies, indicating their key roles in sugarcane resistance to SCMV and SrMV infection. The seven QTLs and nine genes can provide a certain scientific reference to help sugarcane breeders develop varieties resistant to mosaic diseases.

Evaluation of a Sugarcane (Saccharum spp.) Hybrid F1 Population Phenotypic Diversity and Construction of a Rapid Sucrose Yield Estimation Model for Breeding

Sugarcane is the major sugar-producing crop worldwide, and hybrid F₁ populations are the primary populations used in breeding. Challenged by the sugarcane genome's complexity and the sucrose yield's quantitative nature, phenotypic selection is still the most commonly used approach for high-sucrose yield sugarcane breeding. In this study, a hybrid F₁ population containing 135 hybrids was constructed and evaluated for 11 traits (sucrose yield (SY) and its related traits) in a randomized complete-block design during two consecutive growing seasons. The results revealed that all the traits exhibited distinct variation, with the coefficient of variation (CV) ranging from 0.09 to 0.35, the Shannon-Wiener diversity index (H') ranging between 2.64 and 2.98, and the broad-sense heritability ranging from 0.75 to 0.84. Correlation analysis revealed complex correlations between the traits, with 30 trait pairs being significantly correlated. Eight traits, including stalk number (SN), stalk diameter (SD), internode length (IL), stalk height (SH), stalk weight (SW), Brix (B), sucrose content (SC), and yield (Y), were significantly positively correlated with sucrose yield (SY). Cluster analysis based on the 11 traits divided the 135 F₁ hybrids into three groups, with 55 hybrids in Group I, 69 hybrids in Group II, and 11 hybrids in Group III. The principal component analysis indicated that the values of the first four major components' vectors were greater than 1 and the cumulative contribution rate reached 80.93%. Based on the main component values of all samples, 24 F₁ genotypes had greater values than the high-yielding parent 'ROC22' and were selected for the next breeding stage. A rapid sucrose yield estimation equation was established using four easily measured sucrose yield-related traits through multivariable linear stepwise regression. The model was subsequently confirmed using 26 sugarcane cultivars and 24 F₁ hybrids. This study concludes that the sugarcane F₁ population holds great genetic diversity in sucrose yield-related traits. The sucrose yield estimation model, ySY=2.01xSN+8.32xSD+0.79xB+3.44xSH-47.64, can aid to breed sugarcane varieties with high sucrose yield.

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.

SSR-based genome-wide association study in turkish durum wheat germplasms revealed novel QTL of accumulated platinum

BACKGROUND

Durum wheat has a genetic capacity to accumulate toxic metals that can exceed the safety limit of the international standards, which may seriously affect human health. Identifying germplasms with low, nontoxic accumulated metal contents is important to select and develop new varieties. Thus, the objective of this study is to identify the levels of accumulated platinum in durum wheat and detect novel QTL.

METHODS AND RESULTS

Platinum contents were determined using 130 durum genotypes. Results generally showed low values of accumulated Pt and significantly less than the maximum grain's Pt content determined by international standards. Pt contents among genotypes varied from ≤ 0.001 to 0.72 µg/kg with an average of 0.02. Landraces showed the lowest average accumulated Pt. GWAS was then performed with 780 SSR markers. Five QTL were detected and explained 14.4-23.1% of the total phenotypic variation. Chromosomes 3 A, 3B, and 5B appear to be hotspots and may play a crucial role in accumulated Pt and were harbored in 1, 3, and 1 QTL, respectively.

CONCLUSIONS

This assessment of accumulated Pt within a unique panel included accessions mostly from Turkish regions, and GWAS used is the first study regarding accumulated Pt indices to reveal novel QTL. It will allow breeders to accelerate their selection of proper genotypes according to desired alleles and offer an opportunity to apply MAS to minimize Pt toxicity in durum wheat. Results indicated that the significance of genome (B) regions are likely related to the inheritance control of Pt content and may play a pivotal role regarding durum wheat's Pt contents. Nonetheless, these novel QTL should be validated in independent populations in numerous environments.

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.

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.

Genomic Selection in Sugarcane: Current Status and Future Prospects

Sugarcane is a C4 and agro-industry-based crop with a high potential for biomass production. It serves as raw material for the production of sugar, ethanol, and electricity. Modern sugarcane varieties are derived from the interspecific and intergeneric hybridization between Saccharum officinarum, Saccharum spontaneum, and other wild relatives. Sugarcane breeding programmes are broadly categorized into germplasm collection and characterization, pre-breeding and genetic base-broadening, and varietal development programmes. The varietal identification through the classic breeding programme requires a minimum of 12-14 years. The precise phenotyping in sugarcane is extremely tedious due to the high propensity of lodging and suckering owing to the influence of environmental factors and crop management practices. This kind of phenotyping requires data from both plant crop and ratoon experiments conducted over locations and seasons. In this review, we explored the feasibility of genomic selection schemes for various breeding programmes in sugarcane. The genetic diversity analysis using genome-wide markers helps in the formation of core set germplasm representing the total genomic diversity present in the Saccharum gene bank. The genome-wide association studies and genomic prediction in the Saccharum gene bank are helpful to identify the complete genomic resources for cane yield, commercial cane sugar, tolerances to biotic and abiotic stresses, and other agronomic traits. The implementation of genomic selection in pre-breeding, genetic base-broadening programmes assist in precise introgression of specific genes and recurrent selection schemes enhance the higher frequency of favorable alleles in the population with a considerable reduction in breeding cycles and population size. The integration of environmental covariates and genomic prediction in multi-environment trials assists in the prediction of varietal performance for different agro-climatic zones. This review also directed its focus on enhancing the genetic gain over time, cost, and resource allocation at various stages of breeding programmes.

Genome-wide approaches for the identification of markers and genes associated with sugarcane yellow leaf virus resistance

Sugarcane yellow leaf (SCYL), caused by the sugarcane yellow leaf virus (SCYLV) is a major disease affecting sugarcane, a leading sugar and energy crop. Despite damages caused by SCYLV, the genetic base of resistance to this virus remains largely unknown. Several methodologies have arisen to identify molecular markers associated with SCYLV resistance, which are crucial for marker-assisted selection and understanding response mechanisms to this virus. We investigated the genetic base of SCYLV resistance using dominant and codominant markers and genotypes of interest for sugarcane breeding. A sugarcane panel inoculated with SCYLV was analyzed for SCYL symptoms, and viral titer was estimated by RT-qPCR. This panel was genotyped with 662 dominant markers and 70,888 SNPs and indels with allele proportion information. We used polyploid-adapted genome-wide association analyses and machine-learning algorithms coupled with feature selection methods to establish marker-trait associations. While each approach identified unique marker sets associated with phenotypes, convergences were observed between them and demonstrated their complementarity. Lastly, we annotated these markers, identifying genes encoding emblematic participants in virus resistance mechanisms and previously unreported candidates involved in viral responses. Our approach could accelerate sugarcane breeding targeting SCYLV resistance and facilitate studies on biological processes leading to this trait.

Modern Biotechnologies: Innovative and Sustainable Approaches for the Improvement of Sugarcane Tolerance to Environmental Stresses

Sugarcane (Saccharum spp.) is one of the most important industrial cash crops, contributing to the world sugar industry and biofuel production. It has been cultivated and improved from prehistoric times through natural selection and conventional breeding and, more recently, using the modern tools of genetic engineering and biotechnology. However, the heterogenicity, complex poly-aneuploid genome and susceptibility of sugarcane to different biotic and abiotic stresses represent impediments that require us to pay greater attention to the improvement of the sugarcane crop. Compared to traditional breeding, recent advances in breeding technologies (molecular marker-assisted breeding, sugarcane transformation, genome-editing and multiple omics technologies) can potentially improve sugarcane, especially against environmental stressors. This article will focus on efficient modern breeding technologies, which provide crucial clues for the engineering of sugarcane cultivars resistant to environmental stresses.

Pairwise linkage disequilibrium estimation for polyploids

Many tasks in statistical genetics involve pairwise estimation of linkage disequilibrium (LD). The study of LD in diploids is mature. However, in polyploids, the field lacks a comprehensive characterization of LD. Polyploids also exhibit greater levels of genotype uncertainty than diploids, yet no methods currently exist to estimate LD in polyploids in the presence of such genotype uncertainty. Furthermore, most LD estimation methods do not quantify the level of uncertainty in their LD estimates. Our study contains three major contributions. (i) We characterize haplotypic and composite measures of LD in polyploids. These composite measures of LD turn out to be functions of common statistical measures of association. (ii) We derive procedures to estimate haplotypic and composite LD in polyploids in the presence of genotype uncertainty. We do this by estimating LD directly from genotype likelihoods, which may be obtained from many genotyping platforms. (iii) We derive standard errors of all LD estimators that we discuss. We validate our methods on both real and simulated data. Our methods are implemented in the R package ldsep, available on the Comprehensive R Archive Network https://cran.r-project.org/package=ldsep.

Novel Tools for Adjusting Spatial Variability in the Early Sugarcane Breeding Stage

The detection of spatial variability in field trials has great potential for accelerating plant breeding progress due to the possibility of better controlling non-genetic variation. Therefore, we aimed to evaluate a digital soil mapping approach and a high-density soil sampling procedure for identifying and adjusting spatial dependence in the early sugarcane breeding stage. Two experiments were conducted in regions with different soil classifications. High-density sampling of soil physical and chemical properties was performed in a regular grid to investigate the structure of spatial variability. Soil apparent electrical conductivity (ECa) was measured in both experimental areas with an EM38-MK2^® sensor. In addition, principal component analysis (PCA) was employed to reduce the dimensionality of the physical and chemical soil data sets. After conducting the PCA and obtaining different thematic maps, we determined each experimental plot's exact position within the field. Tons of cane per hectare (TCH) data for each experiment were obtained and analyzed using mixed linear models. When environmental covariates were considered, a previous forward model selection step was applied to incorporate the variables. The PCA based on high-density soil sampling data captured part of the total variability in the data for Experimental Area 1 and was suggested to be an efficient index to be incorporated as a covariate in the statistical model, reducing the experimental error (residual variation coefficient, CVe). When incorporated into the different statistical models, the ECa information increased the selection accuracy of the experimental genotypes. Therefore, we demonstrate that the genetic parameter increased when both approaches (spatial analysis and environmental covariates) were employed.

Applying Molecular Phenotyping Tools to Explore Sugarcane Carbon Potential

Sugarcane (Saccharum spp.), a C₄ grass, has a peculiar feature: it accumulates, gradient-wise, large amounts of carbon (C) as sucrose in its culms through a complex pathway. Apart from being a sustainable crop concerning C efficiency and bioenergetic yield per hectare, sugarcane is used as feedstock for producing ethanol, sugar, high-value compounds, and products (e.g., polymers and succinate), and bioelectricity, earning the title of the world's leading biomass crop. Commercial cultivars, hybrids bearing high levels of polyploidy, and aneuploidy, are selected from a large number of crosses among suitable parental genotypes followed by the cloning of superior individuals among the progeny. Traditionally, these classical breeding strategies have been favoring the selection of cultivars with high sucrose content and resistance to environmental stresses. A current paradigm change in sugarcane breeding programs aims to alter the balance of C partitioning as a means to provide more plasticity in the sustainable use of this biomass for metabolic engineering and green chemistry. The recently available sugarcane genetic assemblies powered by data science provide exciting perspectives to increase biomass, as the current sugarcane yield is roughly 20% of its predicted potential. Nowadays, several molecular phenotyping tools can be applied to meet the predicted sugarcane C potential, mainly targeting two competing pathways: sucrose production/storage and biomass accumulation. Here we discuss how molecular phenotyping can be a powerful tool to assist breeding programs and which strategies could be adopted depending on the desired final products. We also tackle the advances in genetic markers and mapping as well as how functional genomics and genetic transformation might be able to improve yield and saccharification rates. Finally, we review how "omics" advances are promising to speed up plant breeding and reach the unexplored potential of sugarcane in terms of sucrose and biomass production.

Candidate gene association of gene expression data in sugarcane contrasting for sucrose content

Association mapping of gene expression data, generated from transcriptome and proteome studies, provides a means of understanding the functional significance and trait association potential of candidate genes. In this study, we applied candidate gene association mapping to validate sugarcane genes, using data from the starch and sucrose metabolism pathway, transcriptome, and proteome. We performed multiplex PCR targeted amplicon sequencing of 109 candidate genes, using NGS technology. A range of statistical models, both single-locus and multi-locus, were compared for minimization of false positives in association mapping of four sugar-related traits with different heritability. The Fixed and random model Circulating Probability Unification model effectively suppressed false positives for both low- and high-heritability traits. We identified favorable alleles of the candidate genes involved in signalling and transcriptional regulation. The results will support genetic improvement of sugarcane and may help clarify the genetic architecture of sugar-related traits.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

An enriched sugarcane diversity panel for utilization in genetic improvement of sugarcane

Sugarcane crop is important for both sugar and biofuels. A world collection of sugarcane and related grasses (WCSRG) maintained at Miami, FL contains > 1,200 non-redundant clones of different species and genera within the Saccharum complex. However, linkage of undesirable alleles with useful genes in wild species has hindered its efficient utilization in sugarcane breeding. A core collection developed previously with smaller number of clones representing WCSRG did not take into account > 120 wild/exotic clones maintained at the USDA-ARS Sugarcane Research Unit in Houma, Louisiana. Moreover, the genome complexity and sub-tropical to temperate growing climate of Louisiana warrant a region-specific core collection that can be used for base-broadening breeding aimed at efficient introgression of desirable alleles. Genetic diversity of 1,485 clones within WCSRG and Louisiana (commercials, wild/exotic) using 423 SSR alleles showed an average gene diversity (h) at 0.208 among all species groups where Erianthus-like Saccharum species (ELSS), Miscanthus spp., and S. spontaneum each formed a distinct cluster, Saccharum robustum, S. officinarum, hybrid cultivars, and S. edule grouped together in a major cluster, and Saccharum sinense and S. barberi formed distinct grouping. A 309-clone diversity panel (SDP1) was developed that captured the genetic diversity based on the combination of maximum length subtree and manual selection to maximize representation of Louisiana clones and minimize import of clones from Miami. SDP1 shared 324 alleles out of the 423 alleles in the entire population of 1,485 clones and captured the genetic diversity of the entire collection with an average gene diversity (h) at 0.163. The variation within (11–17%) and among (83–89%) the populations in SDP1 were comparable with the entire population of 1,485 clones (9–15% and 85–91%, respectively). The breadth of the genetic variation of SDP1 was exemplified by the intra- and inter-specific diversity of a 190-clone mini-core collection with markers derived from known cold-responsive genes. SDP1 will facilitate genome-wide association studies for identification of trait-specific markers for use in marker-assisted breeding in Louisiana and elsewhere.

Molecular diversity and genetic structure of Saccharum complex accessions

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

Rhythms of Transcription in Field-Grown Sugarcane Are Highly Organ Specific

Circadian clocks improve plant fitness in a rhythmic environment. As each cell has its own circadian clock, we hypothesized that sets of cells with different functions would have distinct rhythmic behaviour. To test this, we investigated whether different organs in field-grown sugarcane follow the same rhythms in transcription. We assayed the transcriptomes of three organs during a day: leaf, a source organ; internodes 1 and 2, sink organs focused on cell division and elongation; and internode 5, a sink organ focused on sucrose storage. The leaf had twice as many rhythmic transcripts (>68%) as internodes, and the rhythmic transcriptomes of the internodes were more like each other than to those of the leaves. Among the transcripts expressed in all organs, only 7.4% showed the same rhythmic pattern. Surprisingly, the central oscillators of these organs - the networks that generate circadian rhythms - had similar dynamics, albeit with different amplitudes. The differences in rhythmic transcriptomes probably arise from amplitude differences in tissue-specific circadian clocks and different sensitivities to environmental cues, highlighted by the sampling under field conditions. The vast differences suggest that we must study tissue-specific circadian clocks in order to understand how the circadian clock increases the fitness of the whole plant.

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.