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.

Sampling strategies for sugarcane using either clonal replicates or diverse genotypes can bias the conclusions of RNA-Seq studies.

A key procedure for ensuring statistical confidence in differential gene expression analyses is to use biological replicates to compare distinct groups. Biological replicates allow the estimation of the residual variation in the gene expression levels among samples of a given experimental condition. In sugarcane, it is possible to obtain an estimate of residual variability at two levels: among samples of distinct genotypes of the same experimental treatment, or clonal replicates of the same genotype. The sequencing costs are often a limitation to leveraging both these levels in the same study, stressing the relevance of efforts to determine an appropriate experimental design. We aim to investigate this question by comparing the transcriptional profiles of young sugarcane culms with different sucrose levels using both sampling strategies. Our results show that clonal replicates provided enough statistical power to identify nearly three times more deferentially expressed genes than the more diverse strategy. However, it resulted in potentially less meaningful biological results, because many of the significant genes were likely related to the particular genotype of choice, rather than representing a common expression profile for the compared groups. This study supports the development of sound experimental designs in new studies regarding differential expression for sugarcane.

Planting Season Impacts Sugarcane Stem Development, Secondary Metabolite Levels, and Natural Antisense Transcription.

To reduce the potentially irreversible environmental impacts caused by fossil fuels, the use of renewable energy sources must be increased on a global scale. One promising source of biomass and bioenergy is sugarcane. The study of this crop's development in different planting seasons can aid in successfully cultivating it in global climate change scenarios. The sugarcane variety SP80-3280 was field grown under two planting seasons with different climatic conditions. A systems biology approach was taken to study the changes on physiological, morphological, agrotechnological, transcriptomics, and metabolomics levels in the leaf +1, and immature, intermediate and mature internodes. Most of the variation found within the transcriptomics and metabolomics profiles is attributed to the differences among the distinct tissues. However, the integration of both transcriptomics and metabolomics data highlighted three main metabolic categories as the principal sources of variation across tissues: amino acid metabolism, biosynthesis of secondary metabolites, and xenobiotics biodegradation and metabolism. Differences in ripening and metabolite levels mainly in leaves and mature internodes may reflect the impact of contrasting environmental conditions on sugarcane development. In general, the same metabolites are found in mature internodes from both "one-year" and "one-and-a-half-year sugarcane", however, some metabolites (i.e., phenylpropanoids with economic value) and natural antisense transcript expression are only detected in the leaves of "one-year" sugarcane.

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.

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.

Field microenvironments regulate crop diel transcript and metabolite rhythms.

Most research in plant chronobiology has been done in laboratory conditions. However, laboratories usually fail to mimic natural conditions and their slight fluctuations, highlighting or obfuscating rhythmicity. High-density crops, such as sugarcane (Saccharum hybrid), generate field microenvironments with specific light and temperature regimes resulting from mutual shading. We measured the metabolic and transcriptional rhythms in the leaves of 4-month-old (4 mo) and 9 mo field-grown sugarcane. Most of the assayed rhythms in 9 mo sugarcane peaked >1 h later than in 4 mo sugarcane, including rhythms of the circadian clock gene, LATE ELONGATED HYPOCOTYL (LHY). We hypothesized that older sugarcane perceives dawn later than younger sugarcane as a consequence of self-shading. As a test, we measured LHY rhythms in plants on the east and the west sides of a field. We also tested if a wooden wall built between lines of sugarcane plants changed their rhythms. The LHY peak was delayed in the plants in the west of the field or beyond the wall; both shaded at dawn. We conclude that plants in the same field may have different phases resulting from field microenvironments, impacting important agronomical traits, such as flowering time, stalk weight and number.

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.

Assessment of Gene Flow to Wild Relatives and Nutritional Composition of Sugarcane in Brazil.

The commercial release of genetically modified organisms (GMO) requires a prior environmental and human/animal health risk assessment. In Brazil, the National Biotechnology Technical Commission (CTNBio) requires a survey of the area of natural occurrence of wild relatives of the GMO in the Brazilian ecosystems to evaluate the possibility of introgressive hybridization between sexually compatible species. Modern sugarcane cultivars, the focus of this study, derive from a series of hybridization and backcrossing events among Saccharum species. The so-called "Saccharum broad sense" group includes around 40 species from a few genera, including Erianthus, found in various tropical regions, particularly South-Eastern Asia. In Brazil, three native species, originally considered to belong to Erianthus, were reclassified as S. angustifolium (Nees) Trin., S. asperum (Nees) Steud., and S. villosum Steud., based on inflorescence morphology. Thus, we have investigated the potential occurrence of gene flow among the Brazilian Saccharum native species and commercial hybrids as a requisite for GMO commercial release. A comprehensive survey was carried out to map the occurrence of the three native Saccharum species in Brazil, concluding that they are sympatric with sugarcane cultivation only from around 14°S southwards, which precludes most Northeastern sugarcane-producing states from undergoing introgression. Based on phenology, we concluded that the Brazilian Saccharum species are unable to outcross naturally with commercial sugarcane since the overlap between the flowering periods of sugarcane and the native species is limited. A phylogenomic reconstruction based on the full plastid genome sequence showed that the three native Saccharum species are the taxa closest to sugarcane in Brazil, being closer than introduced Erianthus or Miscanthus. A 2-year study on eight nutritional composition traits of the 20 main sugarcane cultivars cultivated in Brazil was carried out in six environments. The minimum and maximum values obtained were, in percent: moisture (62.6-82.5); sucrose (9.65-21.76); crude fiber (8.06-21.03); FDN (7.20-20.68); FDA (4.55-16.90); lipids (0.06-1.59); ash (0.08-2.67); and crude protein (0.18-1.18). Besides a considerable amount of genetic variation and plastic responses, many instances of genotype-by-environment interaction were detected.

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.

Time-series expression profiling of sugarcane leaves infected with Puccinia kuehnii reveals an ineffective defense system leading to susceptibility.

KEY MESSAGE: Successful orange rust development on sugarcane can potentially be explained as suppression of the plant immune system by the pathogen or delayed plant signaling to trigger defense responses. Puccinia kuehnii is an obligate biotrophic fungus that infects sugarcane leaves causing a disease called orange rust. It spread out to other countries resulting in reduction of crop yield since its first outbreak. One of the knowledge gaps of that pathosystem is to understand the molecular mechanisms altered in susceptible plants by this biotic stress. Here, we investigated the changes in temporal expression of transcripts in pathways associated with the immune system. To achieve this purpose, we used RNA-Seq to analyze infected leaf samples collected at five time points after inoculation. Differential expression analyses of adjacent time points revealed substantial changes at 12, 48 h after inoculation and 12 days after inoculation, coinciding with the events of spore germination, haustoria post-penetration and post-sporulation, respectively. During the first 24 h, a lack of transcripts involved with resistance mechanisms was revealed by underrepresentation of hypersensitive and defense response related genes. However, two days after inoculation, upregulation of genes involved with immune response regulation provided evidence of some potential defense response. Events related to biotic stress responses were predominantly downregulated in the initial time points, but expression was later restored to basal levels. Genes involved in carbohydrate metabolism showed evidence of repression followed by upregulation, possibly to ensure the pathogen nutritional requirements were met. Our results support the hypothesis that P. kuehnii initially suppressed sugarcane genes involved in plant defense systems. Late overexpression of specific regulatory pathways also suggests the possibility of an inefficient recognition system by a susceptible sugarcane genotype.

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.

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

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

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.

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.

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.

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.

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.

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.