Karen S. Aitken

Sugarcane biotechnology: The challenges and opportunities

SummaryCommercial sugarcane, belonging to the genus Saccharum (Poaceae), is an important industrial crop accounting for nearly 70% of sugar produced worldwide. Compared to other major crops, efforts to improve sugarcane are limited and relatively recent, with the first introduction of interspecific hybrids about 80 yr ago. Progress in traditional breeding of sugareane, a highly polyploid and frequently aneuploid plant, is impeded by its narrow gene pool, complex genome, poor fertility, and the long breeding/selection cycle. These constraints, however, make sugarcane a good candidate for molecular breeding. In the past decade considerable progress has been made in understanding and manipulating the sugarcane genome using various biotechnological and cell biological approaches. Notable among them are the creation of transgenic plants with improved agronomic or other important traits, advances in genomics and molecular markers, and progress in understanding the molecular aspects of sucrose transport and accumulation. More recently, substantial effort has been directed towards developing sugarcane as a biofactory for high-value products. While these achievements are commendable, a greater understanding of the sugarcane genome, and cell and whole plant physiology, will accelerate the implementation of commercially significant biotechnology outcomes. We anticipate that the rapid advancements in molecular biology and emerging biotechnology innovations would play a significant role in the future sugarcane crop improvement programs and offer many new opportunities to develop it as a new-generation industrial crop.

A combination of AFLP and SSR markers provides extensive map coverage and identification of homo(eo)logous linkage groups in a sugarcane cultivar

Sugarcane varieties are complex polyploids carrying in excess of 100 chromosomes and are derived from interspecific hybridisation between the domesticated Saccharum officinarum and the wild relative S. spontaneum. A map was constructed in , an Australian cultivar, from a segregating F1 population, using 40 amplified fragment length polymorphism (AFLP) primer combinations, five randomly amplified DNA fingerprints (RAF) primers and 72 simple sequence repeat (SSR) primers. Using these PCR-based marker systems, we generated 1,365 polymorphic markers, of which 967 (71%) were single-dose (SD) markers. Of these SD 967 markers, 910 were distributed on 116 linkage groups (LGs) with a total map length of 9,058.3 cM. Genome organisation was significantly greater than observed in previously reported maps for Saccharum spp. With the addition of 123 double-dose markers, 36 (3:1) segregating markers and a further five SD markers, 1,074 markers were mapped onto 136 LGs. Repulsion phase linkage detected preferential pairing for 40 LGs, which formed 11 LG pairs and three multi-chromosome pairing groups. Using SSRs, double-dose markers and repulsion phase linkage, we succeeded in forming 127 of the 136 LGs into eight homo(eo)logy groups (HG). Two HGs were each represented by two sets of LGs. These sets of LGs potentially correspond to S. officinarum chromosomes, with each set aligning to either end of one or two larger LGs. The larger chromosomes in the two HGs potentially correspond to S. spontaneum chromosomes. This suggestion is consistent with the different basic chromosome number of the two species that are hybridised to form sugarcane cultivars, S. spontaneum (x=8) and S. officinarum (x=10), and illustrates the structural relationship between the genomes of these two species. The discrepancy of coverage between HGs highlights the difficulty in mapping large parts of the genome.

Targeted single nucleotide polymorphism (SNP) discovery in a highly polyploid plant species using 454 sequencing

Discovering single nucleotide polymorphisms (SNPs) in specific genes in a heterozygous polyploid plant species, such as sugarcane, is challenging because of the presence of a large number of homologues. To discover SNPs for mapping genes of interest, 454 sequencing of 307 polymerase chain reaction (PCR) amplicons (> 59 kb of sequence) was undertaken. One region of a four-gasket sequencing run, on a 454 Genome Sequencer FLX, was used for pooled PCR products amplified from each parent of a quantitative trait locus (QTL) mapping population (IJ76-514 x Q165). The sequencing yielded 96,755 (IJ76-514) and 86,241 (Q165) sequences with perfect matches to a PCR primer used in amplification, with an average sequence depth of approximately 300 and an average read length of 220 bases. Further analysis was carried out on amplicons whose sequences clustered into a single contig using an identity of 80% with the program cap3. In the more polymorphic sugarcane parent (Q165), 94% of amplicons (227/242) had evidence of a reliable SNP--an average of one every 35 bases. Significantly fewer SNPs were found in the pure Saccharum officinarum parent--with one SNP every 58 bases and SNPs in 86% (213/247) of amplicons. Using automatic SNP detection, 1632 SNPs were detected in Q165 sequences and 1013 in IJ76-514. From 225 candidate SNP sites tested, 209 (93%) were validated as polymorphic using the Sequenom MassARRAY system. Amplicon re-sequencing using the 454 system enables cost-effective SNP discovery that can be targeted to genes of interest and is able to perform in the highly challenging area of polyploid genomes.

Associations between DNA markers and resistance to diseases in sugarcane and effects of population substructure.

Association between markers and sugarcane diseases were investigated in a collection of 154 sugarcane clones, consisting of important ancestors or parents, and cultivars. 1,068 polymorphic AFLP and 141 SRR markers were scored across all clones. Data on the four most important diseases in the Australian sugarcane industry were obtained; these diseases being pachymetra root rot (Pachymetra chaunorhiza B.J. Croft & M.W. Dick), leaf scald (Xanthomonas albilineans Dowson), Fiji leaf gall (Fiji disease virus), and smut (Ustilago scitaminea H. & P. Sydow). By a simple regression analysis, association between markers and diseases could be readily detected. However, many of these associations were due to the effects of embedded population structure and random effects. After taking population structure into account, we found that 59% of the phenotypic variation in smut resistance ratings could be accounted for by 11 markers, 32% of variation for leaf scald and pachymetra root rot rating by 4 markers, and 26% of Fiji leaf gall by 5 markers. The results suggest that marker–trait associations can be readily detected in populations generated from modern sugarcane breeding programs. This may be due to special features of past sugarcane breeding programs leading to persistent linkage disequilibrium in modern parental populations.

Quantitative trait loci identified for sugar related traits in a sugarcane (Saccharum spp.) cultivar × Saccharum officinarum population

AbstractThe identification of markers linked to quantitative trait loci (QTLs) for increased sugar accumulation could improve the effectiveness of current breeding strategies in sugarcane. Progeny from a cross between a high sucrose producing cultivar, (denotes Australian plant breeding rights), and a Saccharum officinarum clone, IJ76-514 were grown in two field experiments in different years, and evaluated in the early and mid-season phases of crop maturity, to identify robust QTLs in affecting sucrose content in cane. Using an extensive genetic map constructed for with over 1,000 AFLP and SSR markers, a total of 37 QTLs were identified for brix and pol of which, 16 were detected in both experiments. Of these 37 QTL, 30 were clustered into 12 genomic regions in six of the eight homo(eo)logous groups. Each QTL explained from 3 to 9% of the phenotypic variation observed. Both positive and negative effects were identified and the location of the QTLs on linkage groups belonging to the same homo(eo)logy group suggested that a number of the QTLs were allelic forms of the same genes. Of the 37 QTLs identified, the majority were significant in both early and mature cane, but 8 were identified as early specific QTLs and 9 as mature cane QTLs. In total, 97 interactions were significant (P<10−5) and these were localised to 32 genomic regions of which 6 were detected with both years’ data. Models including all the QTLs explained from 37 to 66% of the total phenotypic variation, depending on the trait. The results will be subsequently applied in marker assisted breeding.

Diversity Arrays Technology effectively reveals DNA polymorphism in a large and complex genome of sugarcane

Diversity Arrays Technology (DArT) provides whole genome profiling for hundreds to thousands of polymorphic markers in a single assay using a high-throughput microarray platform. The presented work aimed to establish DArT genotyping for the genetically challenging genome of sugarcane. Due to the genome complexity of this sugar-producing crop of high economic importance, an application of DArT genotyping to this species required extensive testing and optimization. As the method of genome complexity reduction determines the efficiency of polymorphism identification in DArT, various approaches and several methods were tested, in order to establish the most optimal. The sugarcane DArT markers generated with these established methods identified high genetic differentiation of sugarcane ancestral species from modern cultivars, in agreement with the data available for other types of molecular markers for this crop. The majority of sugarcane DArT markers segregated in a Mendelian fashion and were readily incorporated into the framework genetic map. As the DArT markers are sequence-ready genomic clones, we sequenced 384 clones and found that one-third of sequenced markers came from the transcribed portion of the sugarcane genome. The presented results further validate the potential of DArT technology in providing cost-effective genetic profiles for plants, irrespective of their genome complexity, for effective applications in molecular-assisted breeding, diversity analysis or genetic identity testing.

Construction of a genetic linkage map for Saccharum officinarum incorporating both simplex and duplex markers to increase genome coverage

Saccharum officinarum L. is an octoploid with 80 chromosomes and a basic chromosome number of x = 10. It has high stem sucrose and contributes 80% of the chromosomes to the interspecific sugarcane cultivars that are grown commercially for sucrose. A genetic linkage map was developed for S. officinarum (clone IJ76-514) using a segregating population generated from a cross between Q165 (a commercial sugarcane cultivar) and IJ76-514. In total, 40 AFLP and 72 SSR primer pairs were screened across the population, revealing 595 polymorphic bands inherited from IJ76-514. These 595 markers displayed a frequency distribution different from all other sugarcane genetic maps produced, with only 40% being simplex markers (segregated 1:1). Of these 240 simplex markers, 178 were distributed on 47 linkage groups (LGs) and 62 remained unlinked. With the addition of 234 duplex markers and 80 biparental simplex markers (segregating 3:1), 534 markers formed 123 LGs. Using the multi-allelic SSR markers, repulsion phase linkage, and alignment with the Q165 linkage map, 105 of the 123 LGs could be grouped into 10 homology groups (HGs). These 10 HGs were further assigned to the 8 HGs observed in cultivated sugarcane and S. spontaneum. Analysis of repulsion phase linkage indicated that IJ76-514 is neither a complete autopolyploid nor an allopolyploid. Detection of 28 repulsion linkages that occurred between 6 pairs of LGs located in 4 HGs suggested the occurrence of limited preferential chromosome pairing in this species.

Complete genome sequence of sporisorium scitamineum and biotrophic interaction transcriptome with sugarcane

Sporisorium scitamineum is a biotrophic fungus responsible for the sugarcane smut, a worldwide spread disease. This study provides the complete sequence of individual chromosomes of S. scitamineum from telomere to telomere achieved by a combination of PacBio long reads and Illumina short reads sequence data, as well as a draft sequence of a second fungal strain. Comparative analysis to previous available sequences of another strain detected few polymorphisms among the three genomes. The novel complete sequence described herein allowed us to identify and annotate extended subtelomeric regions, repetitive elements and the mitochondrial DNA sequence. The genome comprises 19,979,571 bases, 6,677 genes encoding proteins, 111 tRNAs and 3 assembled copies of rDNA, out of our estimated number of copies as 130. Chromosomal reorganizations were detected when comparing to sequences of S. reilianum, the closest smut relative, potentially influenced by repeats of transposable elements. Repetitive elements may have also directed the linkage of the two mating-type loci. The fungal transcriptome profiling from in vitro and from interaction with sugarcane at two time points (early infection and whip emergence) revealed that 13.5% of the genes were differentially expressed in planta and particular to each developmental stage. Among them are plant cell wall degrading enzymes, proteases, lipases, chitin modification and lignin degradation enzymes, sugar transporters and transcriptional factors. The fungus also modulates transcription of genes related to surviving against reactive oxygen species and other toxic metabolites produced by the plant. Previously described effectors in smut/plant interactions were detected but some new candidates are proposed. Ten genomic islands harboring some of the candidate genes unique to S. scitamineum were expressed only in planta. RNAseq data was also used to reassure gene predictions.

Comparative genetics in sugarcane enables structured map enhancement and validation of marker-trait associations

As sugarcane is a complex polyaneuploid with many chromosomes, large numbers of markers are required to generate genetic maps with reasonable levels of genome coverage. Comparative mapping was investigated as an approach for both quantitative trait loci (QTL) validation and genetic map enhancement in sugarcane. More than 1000 SSR and AFLP markers were scored in a bi-parental Australian sugarcane population (Q3) that was segregating widely for sugar content-related traits. Two maps were constructed, one for each parent. The Q117 (female) and MQ77-340 (male) maps each contained almost 400 markers distributed onto approximately 100 linkage groups (LGs), of which nearly half could be assigned to homology groups (HGs) on the basis of SSRs. Then, using common SSR and AFLP markers, the two Q3 parental maps were aligned with the maps of the French cultivar, R570, and of the Australian cultivar, Q165A (A denotes variety covered by Australian plant breeding rights). As a result of comparative mapping, all ten HGs in the Q117 map, and all eleven HGs in the MQ77-340 map could be re-assigned to seven of the expected eight sugarcane HGs, revealing that one sugarcane HG was not covered at all in either Q3 parental map, and that other HGs were poorly represented. QTL analysis in the Q3 population identified approximately 75 marker-trait associations (MTAs) from approximately 18 chromosomal regions or putative QTL in each map for three sugar content-related traits. QTL location appeared to be consistent between the 4 maps; two of the eight HGs were observed to contain MTAs for brix in two or three maps, strongly suggesting the location of sugar content-related trait loci in these HGs.

Identification of QTL for resistance and susceptibility to Stagonospora meliloti in autotetraploid lucerne

In eastern Australia and California, USA, one of the major lethal fungal diseases of lucerne (Medicago sativa) is Stagonospora root and crown rot, caused by Stagonospora meliloti. Quantitative trait loci (QTL) involved in resistance and susceptibility to S. meliloti were identified in an autotetraploid lucerne backcross population of 145 individuals. Using regression analysis and interval mapping, we detected one region each on linkage groups 2, 6 and 7 that were consistently associated with disease reaction to S. meliloti in two separate experiments. The largest QTL on linkage group 7, which is associated with resistance to S. meliloti, contributed up to 17% of the phenotypic variation. The QTL located on linkage group 2, which is potentially a resistance allele in repulsion to the markers for susceptibility to S. meliloti, contributed up to 8% of the phenotypic variation. The QTL located on linkage group 6, which is associated with susceptibility to S. meliloti, contributed up to 16% of the phenotypic variation. A further two unlinked markers contributed 5 and 8% of the phenotypic variation, and were detected in only one experiment. A total of 517 simple sequence repeat (SSR) markers from Medicago truncatula were screened on the parents of the mapping population. Only 27 (6%) SSR markers were polymorphic and could be incorporated into the autotetraploid map of M. sativa. This allowed alignment of our M. sativa linkage map with published M. truncatula maps. The markers linked to the QTL we have reported will be useful for marker assisted selection for partial resistance to S. meliloti in lucerne.