The complex polyploid genome architecture of sugarcane.

Sugarcane, the world's most harvested crop by tonnage, has shaped global history, trade and geopolitics, and is currently responsible for 80% of sugar production worldwide1. While traditional sugarcane breeding methods have effectively generated cultivars adapted to new environments and pathogens, sugar yield improvements have recently plateaued2. The cessation of yield gains may be due to limited genetic diversity within breeding populations, long breeding cycles and the complexity of its genome, the latter preventing breeders from taking advantage of the recent explosion of whole-genome sequencing that has benefited many other crops. Thus, modern sugarcane hybrids are the last remaining major crop without a reference-quality genome. Here we take a major step towards advancing sugarcane biotechnology by generating a polyploid reference genome for R570, a typical modern cultivar derived from interspecific hybridization between the domesticated species (Saccharum officinarum) and the wild species (Saccharum spontaneum). In contrast to the existing single haplotype ('monoploid') representation of R570, our 8.7 billion base assembly contains a complete representation of unique DNA sequences across the approximately 12 chromosome copies in this polyploid genome. Using this highly contiguous genome assembly, we filled a previously unsized gap within an R570 physical genetic map to describe the likely causal genes underlying the single-copy Bru1 brown rust resistance locus. This polyploid genome assembly with fine-grain descriptions of genome architecture and molecular targets for biotechnology will help accelerate molecular and transgenic breeding and adaptation of sugarcane to future environmental conditions.

Genetic control of yield related stalk traits in sugarcane.

A major focus of sugarcane variety improvement programs is to increase sugar yield, which can be accomplished by either increasing the sugar content of the cane or by increasing cane yield, as the correlation between these traits is low. We used a cross between an Australian sugarcane variety Q165, and a Saccharum officinarum accession, IJ76-514, to dissect the inheritance of yield-related traits in the complex polyploid sugarcane. A population of 227 individuals was grown in a replicated field trial and evaluated over 3 years for stalk weight, stalk diameter, stalk number, stalk length and total biomass. Over 1,000 AFLP and SSR markers were scored across the population and used to identify quantitative trait loci (QTL). In total, 27 regions were found that were significant at the 5% threshold using permutation tests with at least one trait; individually, they explained from 4 to 10% of the phenotypic variation and up to 46% were consistent across years. With the inclusion of digeneic interactions, from 28 to 60% of the variation was explained for these traits. The 27 genomic regions were located on 22 linkage groups (LGs) in six of the eight homology groups (HGs) indicating that a number of alleles or quantitative trait alleles (QTA) at each QTL contribute to the trait; from one to three alleles had an effect on the traits for each QTL identified. Alleles of a candidate gene, TEOSINTE BRANCHED 1 (TB1), the major gene controlling branching in maize, were mapped in this population using either an SSR or SNP markers. Two alleles showed some association with stalk number, but unlike maize, TB1 is not a major gene controlling branching in sugarcane but only has a minor and variable effect.

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.

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.

Identification of QTL for reaction to three races of Colletotrichum trifolii and further analysis of inheritance of resistance in autotetraploid lucerne.

Anthracnose, caused by Colletotrichum trifolii, is one of the most serious diseases of lucerne worldwide. The disease is managed through deployment of resistant cultivars, but new pathotypes present a challenge to the successful implementation of this strategy. This paper reports the genetic map locations of quantitative trait loci (QTL) for reaction to races 1, 2 and 4 of C. trifolii in a single autotetraploid lucerne clone, designated W126 from the Australian cv. Trifecta. Resistance was mapped in a backcross population of 145 individuals, and reaction was assessed both by spray and injection inoculation of stems. Resistance to injection inoculation with races 1 and 4 was incompletely dominant and closely linked (phenotypic markers 2.2 cM apart); these resistances mapped to a linkage group homologous to Medicago truncatula linkage group 8. When the spray inoculation data were subjected to QTL analysis, the strongest QTL for resistance was located on linkage group 8; six QTL were identified for race 1 and four for race 4. Resistance to race 2 was incompletely recessive; four QTL were identified and these include one QTL on linkage group 4 that was also identified for race 1. Modelling of the interactions between individual QTL and marker effects allowed a total of 52-63% of the phenotypic variation to be described for each of the different races. These markers will have value in breeding lucerne, carrying multiple sources of resistance to the three known races of C. trifolii.

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

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

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 Denotes variety covered by Australian plant breeding rights., 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.