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

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

Genetic 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.

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.

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.

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.

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

Sugarcane (Saccharum spp.) has a complex genome with variable ploidy and frequent aneuploidy, which hampers the understanding of phenotype and genotype relations. Despite this complexity, genome-wide association studies (GWAS) may be used to identify favorable alleles for target traits in core collections and then assist breeders in better managing crosses and selecting superior genotypes in breeding populations. Therefore, in the present study, we used a diversity panel of sugarcane, called the Brazilian Panel of Sugarcane Genotypes (BPSG), with the following objectives: (i) estimate, through a mixed model, the adjusted means and genetic parameters of the five yield traits evaluated over two harvest years; (ii) detect population structure, linkage disequilibrium (LD) and genetic diversity using simple sequence repeat (SSR) markers; (iii) perform GWAS analysis to identify marker-trait associations (MTAs); and iv) annotate the sequences giving rise to SSR markers that had fragments associated with target traits to search for putative candidate genes. The phenotypic data analysis showed that the broad-sense heritability values were above 0.48 and 0.49 for the first and second harvests, respectively. The set of 100 SSR markers produced 1,483 fragments, of which 99.5% were polymorphic. These SSR fragments were useful to estimate the most likely number of subpopulations, found to be four, and the LD in BPSG, which was stronger in the first 15 cM and present to a large extension (65 cM). Genetic diversity analysis showed that, in general, the clustering of accessions within the subpopulations was in accordance with the pedigree information. GWAS performed through a multilocus mixed model revealed 23 MTAs, six, three, seven, four and three for soluble solid content, stalk height, stalk number, stalk weight and cane yield traits, respectively. These MTAs may be validated in other populations to support sugarcane breeding programs with introgression of favorable alleles and marker-assisted selection.

GBS-based single dosage markers for linkage and QTL mapping allow gene mining for yield-related traits in sugarcane.

BACKGROUND: Sugarcane (Saccharum spp.) is predominantly an autopolyploid plant with a variable ploidy level, frequent aneuploidy and a large genome that hampers investigation of its organization. Genetic architecture studies are important for identifying genomic regions associated with traits of interest. However, due to the genetic complexity of sugarcane, the practical applications of genomic tools have been notably delayed in this crop, in contrast to other crops that have already advanced to marker-assisted selection (MAS) and genomic selection. High-throughput next-generation sequencing (NGS) technologies have opened new opportunities for discovering molecular markers, especially single nucleotide polymorphisms (SNPs) and insertion-deletion (indels), at the genome-wide level. The objectives of this study were to (i) establish a pipeline for identifying variants from genotyping-by-sequencing (GBS) data in sugarcane, (ii) construct an integrated genetic map with GBS-based markers plus target region amplification polymorphisms and microsatellites, (iii) detect QTLs related to yield component traits, and (iv) perform annotation of the sequences that originated the associated markers with mapped QTLs to search putative candidate genes. RESULTS: We used four pseudo-references to align the GBS reads. Depending on the reference, from 3,433 to 15,906 high-quality markers were discovered, and half of them segregated as single-dose markers (SDMs) on average. In addition to 7,049 non-redundant SDMs from GBS, 629 gel-based markers were used in a subsequent linkage analysis. Of 7,678 SDMs, 993 were mapped. These markers were distributed throughout 223 linkage groups, which were clustered in 18 homo(eo)logous groups (HGs), with a cumulative map length of 3,682.04 cM and an average marker density of 3.70 cM. We performed QTL mapping of four traits and found seven QTLs. Our results suggest the presence of a stable QTL across locations. Furthermore, QTLs to soluble solid content (BRIX) and fiber content (FIB) traits had markers linked to putative candidate genes. CONCLUSIONS: This study is the first to report the use of GBS for large-scale variant discovery and genotyping of a mapping population in sugarcane, providing several insights regarding the use of NGS data in a polyploid, non-model species. The use of GBS generated a large number of markers and still enabled ploidy and allelic dosage estimation. Moreover, we were able to identify seven QTLs, two of which had great potential for validation and future use for molecular breeding in sugarcane.

De novo assembly and transcriptome analysis of contrasting sugarcane varieties.

Sugarcane is an important crop and a major source of sugar and alcohol. In this study, we performed de novo assembly and transcriptome annotation for six sugarcane genotypes involved in bi-parental crosses. The de novo assembly of the sugarcane transcriptome was performed using short reads generated using the Illumina RNA-Seq platform. We produced more than 400 million reads, which were assembled into 72,269 unigenes. Based on a similarity search, the unigenes showed significant similarity to more than 28,788 sorghum proteins, including a set of 5,272 unigenes that are not present in the public sugarcane EST databases; many of these unigenes are likely putative undescribed sugarcane genes. From this collection of unigenes, a large number of molecular markers were identified, including 5,106 simple sequence repeats (SSRs) and 708,125 single-nucleotide polymorphisms (SNPs). This new dataset will be a useful resource for future genetic and genomic studies in this species.